劉玉標(biāo),邢云林,申仲翰(中國科學(xué)院力學(xué)研究所流固耦合系統(tǒng)力學(xué)重點(diǎn)實(shí)驗(yàn)室,北京100190)

?

基于應(yīng)變測量的彈體截面載荷識別

劉玉標(biāo),邢云林,申仲翰
(中國科學(xué)院力學(xué)研究所流固耦合系統(tǒng)力學(xué)重點(diǎn)實(shí)驗(yàn)室,北京100190)

摘要:如何通過實(shí)彈試驗(yàn)中測量的導(dǎo)彈彈體應(yīng)變數(shù)據(jù)獲得彈體截面等效載荷是型號設(shè)計(jì)部門面臨的研究課題。引入一種數(shù)值反演方法,利用導(dǎo)彈試驗(yàn)過程中實(shí)測彈體表面應(yīng)變數(shù)據(jù)實(shí)現(xiàn)彈體截面總體載荷識別。為解決實(shí)際測量過程中可能出現(xiàn)的數(shù)據(jù)壞點(diǎn)對識別結(jié)果的干擾,提高載荷識別的可靠性,提出基于線性相關(guān)性理論的測量數(shù)據(jù)壞點(diǎn)剔除方法,設(shè)計(jì)驗(yàn)證試驗(yàn),檢驗(yàn)識別方法的精度和可靠性。試驗(yàn)結(jié)果表明:該識別方法簡單易行,且具有很高的識別精度與工程應(yīng)用價(jià)值。

關(guān)鍵詞:兵器科學(xué)與技術(shù);彈體;載荷識別;線性相關(guān)性;最小二乘法

0　引言

結(jié)構(gòu)所受載荷是結(jié)構(gòu)設(shè)計(jì)的依據(jù),識別結(jié)構(gòu)在服役過程所受的動(dòng)態(tài)載荷是結(jié)構(gòu)優(yōu)化設(shè)計(jì)的前提?；趯?shí)際工程結(jié)構(gòu)受力的復(fù)雜性,通常行之有效的方法是根據(jù)結(jié)構(gòu)系統(tǒng)的實(shí)測動(dòng)態(tài)響應(yīng)和結(jié)構(gòu)系統(tǒng)本身固有動(dòng)態(tài)特性反演結(jié)構(gòu)所受的動(dòng)態(tài)載荷。

潛射導(dǎo)彈[1]在服役過程中受底部推力、發(fā)射筒約束與彈體周圍流體動(dòng)力共同作用,因此彈體結(jié)構(gòu)受力非常復(fù)雜。通常導(dǎo)彈彈體為細(xì)長體結(jié)構(gòu),為了簡化受力分析,可以將彈體結(jié)構(gòu)假定為非均勻梁。彈體截面受到的內(nèi)力載荷主要為彎矩與軸力,扭矩與剪力較小可以不予考慮。如何利用模型或?qū)崗椩囼?yàn)手段,通過彈體響應(yīng)測量數(shù)據(jù)來識別彈體梁截面所受的彎矩與軸力是型號設(shè)計(jì)部門面臨的研究課題。

在固體火箭截面彎矩識別研究方面,尹云玉[2-3]采用模態(tài)疊加原理,通過模態(tài)彎矩來識別截面橫向響應(yīng)與彎矩。其中彈體結(jié)構(gòu)位移數(shù)據(jù)是利用固體火箭自由飛行狀態(tài)時(shí)的低頻振動(dòng)加速度測量數(shù)據(jù)進(jìn)行兩次積分得到,但是該方法在低頻時(shí)位移積分誤差較大,即使采用積分補(bǔ)償也難以解決。另外該方法只適用于自由飛行的固體火箭橫向載荷識別,尚不適用于發(fā)射過程帶有約束飛行的導(dǎo)彈彈體載荷識別。在結(jié)構(gòu)受力測量技術(shù)方面,測力天平技術(shù)的發(fā)展非常迅速,在風(fēng)洞試驗(yàn)中已作為常規(guī)測量設(shè)備得到廣泛的應(yīng)用。應(yīng)變天平[4-6]是通過天平結(jié)構(gòu)的應(yīng)變測量組成電橋,通過標(biāo)定試驗(yàn)建立內(nèi)力載荷與應(yīng)變電橋之間的關(guān)系,可以直接測量出飛行器在風(fēng)洞中受到的6分量外力的大小和方向。測力天平既需要有較高的靈敏度與信噪比以保證測量結(jié)果準(zhǔn)確,又需要有較大的剛度和強(qiáng)度以保證被測結(jié)構(gòu)的安全[7]。如采用天平作為測量儀器測量結(jié)構(gòu)內(nèi)力,需要將天平安裝在被測結(jié)構(gòu)物內(nèi)部,這將改變被測結(jié)構(gòu)物的靜動(dòng)態(tài)響應(yīng)特性和風(fēng)洞流場,對于潛射導(dǎo)彈結(jié)構(gòu)實(shí)施難度較大。在飛機(jī)機(jī)翼內(nèi)力識別研究方面,使用比較多的是應(yīng)變測量法。其測量原理與應(yīng)變天平基本相同,只是將應(yīng)變片直接粘貼在機(jī)翼測點(diǎn)位置組成電橋,使用回歸分析法[8-9]和神經(jīng)網(wǎng)絡(luò)法[10-11]通過由地面校準(zhǔn)試驗(yàn)建立內(nèi)力載荷與應(yīng)變電橋之間的關(guān)系,即載荷方程。然后將飛行試驗(yàn)測得的電橋響應(yīng)數(shù)據(jù)帶入載荷方程中得到內(nèi)力載荷。載荷方程的確定方法可以等。應(yīng)變測量法中由于一種電橋主要對一種載荷敏感,因此必須針對不同的內(nèi)力載荷組成多種電橋進(jìn)行測量,且測量電橋靈敏度對結(jié)構(gòu)載荷識別精度較為敏感[8]。

本文通過應(yīng)變片直接測量彈體某一截面上多個(gè)角度的軸向應(yīng)變數(shù)據(jù),利用地面試驗(yàn)預(yù)先標(biāo)定的彈體截面彎矩、軸力-應(yīng)變轉(zhuǎn)換系數(shù)矩陣,在時(shí)域內(nèi)通過數(shù)值反演得到導(dǎo)彈彈體截面上受到的彎矩和軸力載荷。該方法僅使用應(yīng)變片作為測量設(shè)備,不會(huì)改變潛射導(dǎo)彈試驗(yàn)時(shí)的外部流場特征及被測結(jié)構(gòu)的受力狀態(tài)。各個(gè)截面的應(yīng)變片無需組成電橋,結(jié)構(gòu)載荷識別精度不受測量電橋靈敏度影響。且該方法適用性不受導(dǎo)彈飛行過程中有無約束條件的限制,識別結(jié)果無模態(tài)截?cái)嗾`差[2-3]。

1　載荷識別方法

1.1彈體截面彎矩與軸力識別

通常導(dǎo)彈彈體為圓筒形細(xì)長體結(jié)構(gòu),本文假定導(dǎo)彈彈體為歐拉梁。根據(jù)材料力學(xué)理論,對于一條僅受軸力、剪力、彎矩載荷,且符合歐拉假定的梁,梁截面上任意一點(diǎn)的軸向應(yīng)變,由兩部分組成:

式中:My、Mz分別為繞y和z軸的彎矩;y和z為應(yīng)變測點(diǎn)的坐標(biāo)值;E為材料彈性模量;F為軸力;S為截面面積;Iy、Iz為截面慣性矩。

如在梁上某一截面任意周向角θ上布置兩個(gè)應(yīng)變片測點(diǎn)(測點(diǎn)1和測點(diǎn)2)測量梁的軸向應(yīng)變,且假定測量應(yīng)變片已進(jìn)行了溫度補(bǔ)償,測量應(yīng)變數(shù)據(jù)中不含溫度引起的應(yīng)變。如圖1所示,對于圓筒形截面:

圖1　彈體梁截面應(yīng)變測點(diǎn)Fig.1 Strain measuring points on the cross-section

由彎矩引起的軸向應(yīng)變εM(θ)可根據(jù)(5)式獲得:

式中:ε1(θ)、ε2(θ)分別表示測點(diǎn)1與測點(diǎn)2測量的軸向應(yīng)變,測點(diǎn)1和測點(diǎn)2稱之為與周向角θ對應(yīng)的應(yīng)變測量對。如在一截面上相同半徑r位置兩個(gè)不同周向角上(θ1,θ2)布置軸向應(yīng)變測量對,則有

式中:A為彎矩-應(yīng)變轉(zhuǎn)換系數(shù)矩陣。理論上講,每一截面僅需要2個(gè)軸向應(yīng)變測量對,共4個(gè)測點(diǎn)的應(yīng)變值就可以解出截面彎矩My與Mz.

彈體飛行過程中,應(yīng)變片可能損壞,2個(gè)軸向應(yīng)變測量對的載荷識別方案的精度與可靠性都無法保證。為了提高載荷識別結(jié)果的精度與可靠性,在實(shí)際應(yīng)用中可布置多個(gè)角度軸向應(yīng)變測量對。圖2表示4個(gè)角度的軸向應(yīng)變測量方案的測點(diǎn)位置示意圖,A-C、B-D、E-G、F-H分別為4個(gè)角度上的應(yīng)變測點(diǎn)。如實(shí)際測量方案采用n個(gè)角度軸向應(yīng)變測量對,共2n個(gè)測點(diǎn)。根據(jù)彎矩與應(yīng)變的關(guān)系式,則有

式中:Ay(θ)、Az(θ)為對應(yīng)測點(diǎn)的彎矩-應(yīng)變轉(zhuǎn)換系數(shù)。

圖2　截面上4個(gè)角度測量對應(yīng)的應(yīng)變片位置Fig.2 Strain measuring point pairs on the cross- section

(7)式為線性方程組,由于方程的個(gè)數(shù)大于未知數(shù)的個(gè)數(shù),只能通過矩陣偽逆得出最小二乘誤差意義下方程的解[12]:

式中:AL += [ATA]-1AT是A的左偽逆。

由于彎矩引起的應(yīng)變在相加時(shí)互相抵消,則軸力可以由(9)式求得:

式中:Tχ=1/ ES為軸力-應(yīng)變轉(zhuǎn)換系數(shù)。

1.2載荷-應(yīng)變轉(zhuǎn)換系數(shù)標(biāo)定

由(8)式和(9)式可知,本方法彎矩、軸力識別精度取決于彎矩-應(yīng)變轉(zhuǎn)換系數(shù)以及軸力-應(yīng)變轉(zhuǎn)換系數(shù)的計(jì)算精度。通常由于實(shí)際導(dǎo)彈結(jié)構(gòu)的復(fù)雜性,很難通過理論計(jì)算精確確定。為確保載荷識別精度,應(yīng)預(yù)先進(jìn)行標(biāo)定試驗(yàn)獲得載荷-應(yīng)變轉(zhuǎn)換系數(shù)。

標(biāo)定試驗(yàn)步驟如下:1)將貼好應(yīng)變片的導(dǎo)彈彈體一端固定,一端加載純彎彎矩(或軸力);2)在彈體結(jié)構(gòu)強(qiáng)度容許的范圍內(nèi),采用逐漸遞增(或遞減)純彎彎矩(或軸力)的加載方式,獲取一系列截面彎矩(或軸力)載荷及與相對應(yīng)的應(yīng)變數(shù)據(jù);3)通過線性回歸方法擬合得到彎矩-應(yīng)變載荷轉(zhuǎn)換系數(shù)以及軸力-應(yīng)變載荷轉(zhuǎn)換系數(shù)。

1.3壞點(diǎn)剔除方法

在實(shí)際試驗(yàn)測量過程中,不可避免會(huì)出現(xiàn)應(yīng)變片損壞的情況。應(yīng)變片損壞后,測量得到的應(yīng)變結(jié)果難以預(yù)料。如果載荷識別時(shí)利用了此壞點(diǎn)的測量數(shù)據(jù)將對結(jié)果帶來難以預(yù)測的誤差,降低載荷識別結(jié)果的可靠性。所以在進(jìn)行載荷識別之前,對壞點(diǎn)數(shù)據(jù)進(jìn)行剔除是非常必要的。工程上一般是通過查看測點(diǎn)數(shù)據(jù)的幅值、頻率等識別壞點(diǎn),具有較強(qiáng)的主觀性。本文基于測量數(shù)據(jù)的線性相關(guān)性提出如下的測量數(shù)據(jù)壞點(diǎn)剔除原則,可自動(dòng)識別出壞點(diǎn)。其基本思路為:由于同一截面應(yīng)變數(shù)據(jù)是由同一截面彎矩、軸力載荷所形成,它們之間應(yīng)是線性相關(guān)的。假定同一截面大部分測量通道數(shù)據(jù)是可靠的,則可通過進(jìn)行互相關(guān)分析,確定各數(shù)據(jù)通道之間互相關(guān)系數(shù)。

任意兩測量通道數(shù)據(jù)X、Y的互相關(guān)系數(shù)定義為RXY(τ).m個(gè)測量通道數(shù)據(jù)Y1,Y2,…,Ym的互相關(guān)矩陣為

因此,對每個(gè)截面各個(gè)測點(diǎn)的應(yīng)變數(shù)據(jù)做互相關(guān)分析。如果同截面數(shù)據(jù)完全線性相關(guān),則互相關(guān)系數(shù)都是1.如果某一測點(diǎn)應(yīng)變片損壞,則測量的應(yīng)變數(shù)據(jù)與其他大部分測點(diǎn)應(yīng)變數(shù)據(jù)相關(guān)性會(huì)很差,互相關(guān)系數(shù)將小于1.互相關(guān)系數(shù)數(shù)值越小,相關(guān)性會(huì)越差,說明是壞點(diǎn)的可能性越大。如果互相關(guān)系數(shù)小于設(shè)定的容許閾值,那么就認(rèn)為該點(diǎn)為壞點(diǎn)。容許閾值視具體試驗(yàn)情況及應(yīng)變片的測量信噪比而定。

由(5)式和(6)式可知,采用本文識別方法進(jìn)行載荷識別時(shí),一個(gè)角度上同時(shí)需要兩個(gè)測點(diǎn)的應(yīng)變數(shù)據(jù)。因此,若某個(gè)角度軸向應(yīng)變測量對中一個(gè)測點(diǎn)被剔除,則其對應(yīng)的另一個(gè)測點(diǎn)也相應(yīng)地被剔除。

2　驗(yàn)證試驗(yàn)

為了驗(yàn)證識別方法的可靠性與識別精度,在實(shí)驗(yàn)室進(jìn)行了上述載荷識別方法的驗(yàn)證試驗(yàn)。

驗(yàn)證試驗(yàn)坐標(biāo)系設(shè)定為:模擬彈體梁固定端面為坐標(biāo)原點(diǎn),梁軸向?yàn)棣?豎直方向?yàn)閥軸,z軸按右手法則生成。驗(yàn)證試驗(yàn)?zāi)M彈體梁采用圓筒形截面鋁材加工,外徑534 mm,壁厚7 mm.將圓筒一端通過固定支座固定于剪力墻上,一端自由并采用集中力方式加載。載荷識別加載方式如圖3所示。

圖3　載荷識別驗(yàn)證試驗(yàn)加載方式Fig.3 Loading mode of load identification test

由于剪力對軸向應(yīng)變影響很小,在進(jìn)行彎矩-應(yīng)變轉(zhuǎn)換系數(shù)標(biāo)定試驗(yàn)時(shí)采用了較為容易實(shí)現(xiàn)的等效彎矩加載方式。如圖4所示,在筒體的懸臂端增加一段彎矩加載梁,在加載梁的兩個(gè)相距為L的作用點(diǎn)分別施加方向相反的作用力F,那么作用在筒體端面上的彎矩為M = FL.

圖4　載荷-應(yīng)變轉(zhuǎn)換系數(shù)標(biāo)定試驗(yàn)彎矩加載方式Fig.4 Moment loading mode of loading calibration test

試驗(yàn)共測量了模型梁的8個(gè)截面。每個(gè)截面分別有8個(gè)單向應(yīng)變片,測點(diǎn)位置分別為周向角(測點(diǎn)徑向與y軸的夾角)0°、30°、45°、90°、180°、210°、225°、270°.應(yīng)變片周向分布如圖5所示,黑色菱形表示應(yīng)變片位置,8個(gè)截面應(yīng)變片周向分布位置完全相同。將所有的應(yīng)變片進(jìn)行編號為xxyyy,其中xx表示為截面編號,而yyy表示應(yīng)變片周向角度,如第2個(gè)截面30°方向的應(yīng)變片編號為02030.

圖5　梁截面上的應(yīng)變片周向分布圖Fig.5 Strain measuring point pairs on the cross-section of beam

根據(jù)1.2節(jié)提出的逐級加載方法,按照圖4的加載方案與1.2節(jié)試驗(yàn)步驟進(jìn)行載荷-應(yīng)變轉(zhuǎn)換系數(shù)標(biāo)定試驗(yàn),得到各個(gè)截面的彎矩(軸力)-應(yīng)變轉(zhuǎn)換系數(shù)如表1所示。

表1　驗(yàn)證試驗(yàn)?zāi)Ｐ蛷澗?應(yīng)變轉(zhuǎn)換系數(shù)Tab.1 The moment-strain conversion factor

在載荷識別試驗(yàn)中,采用圖3所示加載方式,在模擬彈體梁上的自由端施加集中力,作用力為

式中:集中力Fχ、Fy、Fz的單位為kN;時(shí)間t的單位為s.加載過程中同時(shí)測量試驗(yàn)梁上各點(diǎn)的應(yīng)變數(shù)據(jù)。其中截面4上45°角應(yīng)變數(shù)據(jù)如圖6所示。為了驗(yàn)證數(shù)據(jù)剔除方法的可行性,試驗(yàn)前將測點(diǎn)04030的應(yīng)變片人為損壞,故該點(diǎn)測量值中含有較大的噪聲。通過與該截面其他各測點(diǎn)數(shù)據(jù)的互相關(guān)分析,該測點(diǎn)與該截面其他各測點(diǎn)數(shù)據(jù)的互相關(guān)系數(shù)小于0.3,根據(jù)1.3節(jié)提出的壞點(diǎn)剔除原則,此測點(diǎn)及對應(yīng)的軸向應(yīng)變測量對中另一測點(diǎn)04210的數(shù)據(jù)被剔除。

圖6　應(yīng)變片04045的應(yīng)變Fig.6 Strain at measuring point 04045

圖7和圖8分別給出了模擬彈體梁的截面4軸力與彎矩載荷識別值與實(shí)際加載值的比較。由圖7和圖8可以看出,如果將壞點(diǎn)04030點(diǎn)測量數(shù)據(jù)剔除,則識別值與實(shí)際加載值基本一致,平均相對誤差小于1%,在峰值處相對誤差小于2%.圖9表示模擬彈體梁截面4彎矩載荷識別值在數(shù)據(jù)壞點(diǎn)剔除前后的對比。如果不將壞點(diǎn)剔除,則識別值與實(shí)際加載值存在較大的誤差,平均相對誤差超過10%,在峰值處相對誤差大于50%.由此表明,在載荷識別過程中,數(shù)據(jù)壞點(diǎn)剔除的是十分必要的。

圖7　截面4軸力Fχ時(shí)域曲線對比Fig.7 Axial force Fχon cross-section 4

圖8　截面4彎矩My時(shí)域曲線對比Fig.8 Moment Myon cross-section 4

圖9　截面4彎矩Mz時(shí)域曲線對比Fig.9 Moment Mzon cross-section 4

3　結(jié)論

本文引入了一種數(shù)值反演方法,利用導(dǎo)彈試驗(yàn)過程中實(shí)測彈體表面應(yīng)變數(shù)據(jù)實(shí)現(xiàn)了彈體截面總體載荷的識別,同時(shí)為了降低實(shí)際測量過程可能出現(xiàn)的數(shù)據(jù)壞點(diǎn)對識別結(jié)果的干擾,提出了基于線性相關(guān)性理論的測量數(shù)據(jù)壞點(diǎn)剔除方法,并在實(shí)驗(yàn)室構(gòu)建模擬試驗(yàn)對方法進(jìn)行了驗(yàn)證。

本文方法由于只需測量彈體表面應(yīng)變,測量方案和反演計(jì)算程序簡單易于實(shí)現(xiàn),無需構(gòu)建復(fù)雜的測量橋路,識別方法精度不受測量橋路靈敏度影響,且不受彈體局部結(jié)構(gòu)形式復(fù)雜性及試驗(yàn)中彈體所受約束條件的限制,同時(shí)測量設(shè)備不會(huì)影響潛射導(dǎo)彈試驗(yàn)時(shí)的外流場特征及被測結(jié)構(gòu)的受力狀態(tài)。通過型號設(shè)計(jì)部門實(shí)彈試驗(yàn)數(shù)據(jù)分析處理的實(shí)際應(yīng)用表明,該方法具有很高的識別精度與工程實(shí)用價(jià)值,且本文提出的數(shù)據(jù)壞點(diǎn)剔除方法可以快速辨別壞點(diǎn)數(shù)據(jù),進(jìn)一步提高載荷識別結(jié)果的精度與可靠性。

參考文獻(xiàn)(References)

[1]劉玉標(biāo),劉維瑋,程載斌,等.通過縮比試驗(yàn)數(shù)據(jù)計(jì)算潛射導(dǎo)彈的載荷[J].兵工學(xué)報(bào), 2012.33(8): 916-920.LIU Yu-biao, LIU Wei-wei, CHENG Zai-bin, et al.Load calculation for submarine launched missile based on scaled model test [J].Acta Armamentarii, 2012.33(8): 916-920.(in Chinese)

[2]尹云玉.固體火箭橫向響應(yīng)載荷識別方法[J].宇航學(xué)報(bào), 2009, 30(6): 2135-2139.YIN Yun-yu.The identification method of transverse responding load for solid rocket [J].Journal of Astronautics, 2009, 30(6): 2135-2139.(in Chinese)

[3]尹云玉.固體火箭自由運(yùn)行狀態(tài)截面彎矩識別方法[J].宇航學(xué)報(bào), 2010, 31(7): 1706-1710.YIN Yun-yu.An identification method for transverse response moment of solid rocket during free flight state[J].Journal of Astronautics, 2010, 31(7):1706-1710.(in Chinese)

[4]彭超,史玉杰,陳竹.某飛機(jī)部件高速風(fēng)洞測力天平研制[J].實(shí)驗(yàn)流體力學(xué), 2009, 23(4):98-101.PENG Chao, SHI Yu-jie, CHEN Zhu.Forcing balance development of a certain airplane components in high speed wind tunnel [J].Journal of Experiments in Fluid Mechanics, 2009, 23(4): 98-101.(in Chinese)

[5]胡凱征,劉獻(xiàn)偉,黃帥軍,等.導(dǎo)彈掛飛應(yīng)力天平測力系統(tǒng)設(shè)計(jì)[J].科學(xué)技術(shù)與工程, 2008, 8(20):5716-5719.HU Kai-zheng, LIU Xian-wei, HUANG Shuai-jun, et al.Design and development of system for the measurement of missile’s embarkation flight loads based on strain-gauge balance [ J].Science Technology and Engineering, 2008,8(20):5716-5719.(in Chinese)

[6]楊雙龍.風(fēng)洞應(yīng)變天平動(dòng)態(tài)特性與動(dòng)態(tài)校正方法研究[D].合肥:合肥工業(yè)大學(xué), 2014.YANG Shuang-long.Studies on dynamic characteristics and dynamic correction methods for wind tunnel strain gauge balance [D].Hefei: Hefei University of Technology, 2014.(in Chinese)

[7]解亞君,宋筆鋒,郭琦,等.桿式應(yīng)變天平的優(yōu)化設(shè)計(jì)[J].傳感器技術(shù), 2002, 21(12):18-20.XIE Ya-jun, SONG Bi-feng, GUO Qi, et al.Optimizing design of beam strain-gauge balance[J].Journal of Transducer Technology, 2002, 21(12):18-20.(in Chinese)

[8]湯吉晨.飛機(jī)尾翼載荷飛行測量研究[J].航空學(xué)報(bào), 1989, 10(10):B474-B478.TANG Ji-chen.Investigation on the determination of airplane tail loads by flight tests[J].Acta Aeronautica et Astronautica Sinica, 1989, 10(10):B474-B478.(in Chinese)

[9]胡昌榮.飛機(jī)飛行載荷測量載荷方程的優(yōu)選[J].航空學(xué)報(bào), 1994,15(1):102-105.HU Chang-rong.The optimization of load equation during aircraft flight load measurement[J].Acta Aeronautica et Astronautica Sinica, 1994, 15(1):102-105.(in Chinese)

[10]尚琳,何發(fā)東,楊華保,等.基于神經(jīng)網(wǎng)絡(luò)的垂尾飛行載荷模型研究[J].強(qiáng)度與環(huán)境, 2011, 38(5):49-53.SHANG Lin, HE Fa-dong, YANG Hua-bao, et al.Research on flight loads model of vertical tail based on neural networks[J].Structure & Environment Engineering, 2011, 38(5):49-53.(in Chinese)

[11]Allen M J, Dibley R P.Modeling aircraft wing loads from flight data using neural networks[R].Edwards, CA:SAE, 2003.

[12]方保镕,周繼東,李醫(yī)民.矩陣論[M].北京:清華大學(xué)出版社, 2009.FANG Bao-rong, ZHOU Ji-dong, LI Yi-min.Matrix theory[M].Beijing:Tsinghua University Press, 2009.(in Chinese)

Identification of Cross-section Loads Based on Measured Strain of Missile Body

LIU Yu-biao, XING Yun-lin, SHEN Zhong-han
(Key Laboratory for Mechanics in Fluid Solid Coupling Systems, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China)

Abstract:A pressing problem facing to missile designer is to identify the cross section loads of missile load by measured strains data which would be conveniently obtained in the flight test.A load identification method based on the measured strain data of the surface of the missile body is introduced to determine the equivalent bending moment and the axial force of the cross section.In order to eliminate the interference of the fault points of the measured data, a defect rejection principle based on linear correlation theory is presented to improve the reliability of the identification method.A verification test is designed to verify the accuracy and reliability of the identification method.Results show that the method has the high identification accuracy and can be applied to the engineering practice.

Key words:ordnance science and technology; missile body; load identification; linear correlation; least square method

作者簡介:劉玉標(biāo)(1961—),男,副研究員,碩士生導(dǎo)師。E-mail:yubiao@ imech.ac.cn

收稿日期:2015-04-06

DOI:10.3969/ j.issn.1000-1093.2016.02.020

中圖分類號:TB122

文獻(xiàn)標(biāo)志碼:A

文章編號:1000-1093(2016)02-0332-06