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

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

?

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

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

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

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

0　引言

結(jié)構(gòu)所受載荷是結(jié)構(gòu)設(shè)計(jì)的依據(jù),識(shí)別結(jié)構(gòu)在服役過(guò)程所受的動(dòng)態(tài)載荷是結(jié)構(gòu)優(yōu)化設(shè)計(jì)的前提?；趯?shí)際工程結(jié)構(gòu)受力的復(fù)雜性,通常行之有效的方法是根據(jù)結(jié)構(gòu)系統(tǒng)的實(shí)測(cè)動(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]在服役過(guò)程中受底部推力、發(fā)射筒約束與彈體周圍流體動(dòng)力共同作用,因此彈體結(jié)構(gòu)受力非常復(fù)雜。通常導(dǎo)彈彈體為細(xì)長(zhǎng)體結(jié)構(gòu),為了簡(jiǎn)化受力分析,可以將彈體結(jié)構(gòu)假定為非均勻梁。彈體截面受到的內(nèi)力載荷主要為彎矩與軸力,扭矩與剪力較小可以不予考慮。如何利用模型或?qū)崗椩囼?yàn)手段,通過(guò)彈體響應(yīng)測(cè)量數(shù)據(jù)來(lái)識(shí)別彈體梁截面所受的彎矩與軸力是型號(hào)設(shè)計(jì)部門面臨的研究課題。

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

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

1　載荷識(shí)別方法

1.1彈體截面彎矩與軸力識(shí)別

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

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

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

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

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

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

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

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

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

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

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

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

由于彎矩引起的應(yīng)變?cè)谙嗉訒r(shí)互相抵消,則軸力可以由(9)式求得:

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

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

由(8)式和(9)式可知,本方法彎矩、軸力識(shí)別精度取決于彎矩-應(yīng)變轉(zhuǎn)換系數(shù)以及軸力-應(yīng)變轉(zhuǎn)換系數(shù)的計(jì)算精度。通常由于實(shí)際導(dǎo)彈結(jié)構(gòu)的復(fù)雜性,很難通過(guò)理論計(jì)算精確確定。為確保載荷識(shí)別精度,應(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),采用逐漸遞增(或遞減)純彎彎矩(或軸力)的加載方式,獲取一系列截面彎矩(或軸力)載荷及與相對(duì)應(yīng)的應(yīng)變數(shù)據(jù);3)通過(guò)線性回歸方法擬合得到彎矩-應(yīng)變載荷轉(zhuǎn)換系數(shù)以及軸力-應(yīng)變載荷轉(zhuǎn)換系數(shù)。

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

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

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

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

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

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

為了驗(yàn)證識(shí)別方法的可靠性與識(shí)別精度,在實(shí)驗(yàn)室進(jìn)行了上述載荷識(shí)別方法的驗(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.將圓筒一端通過(guò)固定支座固定于剪力墻上,一端自由并采用集中力方式加載。載荷識(shí)別加載方式如圖3所示。

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

由于剪力對(duì)軸向應(yīng)變影響很小,在進(jìn)行彎矩-應(yīng)變轉(zhuǎn)換系數(shù)標(biāo)定試驗(yàn)時(shí)采用了較為容易實(shí)現(xiàn)的等效彎矩加載方式。如圖4所示,在筒體的懸臂端增加一段彎矩加載梁,在加載梁的兩個(gè)相距為L(zhǎng)的作用點(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)共測(cè)量了模型梁的8個(gè)截面。每個(gè)截面分別有8個(gè)單向應(yīng)變片,測(cè)點(diǎn)位置分別為周向角(測(cè)點(diǎn)徑向與y軸的夾角)0°、30°、45°、90°、180°、210°、225°、270°.應(yīng)變片周向分布如圖5所示,黑色菱形表示應(yīng)變片位置,8個(gè)截面應(yīng)變片周向分布位置完全相同。將所有的應(yīng)變片進(jìn)行編號(hào)為xxyyy,其中xx表示為截面編號(hào),而yyy表示應(yīng)變片周向角度,如第2個(gè)截面30°方向的應(yīng)變片編號(hào)為02030.

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

根據(jù)1.2節(jié)提出的逐級(jí)加載方法,按照?qǐng)D4的加載方案與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

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

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

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

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

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

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

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

3　結(jié)論

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

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

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

[1]劉玉標(biāo),劉維瑋,程載斌,等.通過(guò)縮比試驗(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)載荷識(shí)別方法[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)截面彎矩識(shí)別方法[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)洞測(cè)力天平研制[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)力天平測(cè)力系統(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ī)尾翼載荷飛行測(cè)量研究[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ī)飛行載荷測(cè)量載荷方程的優(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

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

收稿日期:2015-04-06

DOI:10.3969/ j.issn.1000-1093.2016.02.020

中圖分類號(hào):TB122

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

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