基于真實(shí)形狀鳥(niǎo)體的撞擊方向?qū)B(niǎo)撞分析影響研究

基于真實(shí)形狀鳥(niǎo)體的撞擊方向?qū)B(niǎo)撞分析影響研究

張大海1,2，費(fèi)慶國(guó)1,2，劉宏月1,2

(1.東南大學(xué)工程力學(xué)系，南京210096； 2.江蘇省工程力學(xué)分析重點(diǎn)實(shí)驗(yàn)室，南京210096)

摘要：鳥(niǎo)體形狀及撞擊方向?qū)B(niǎo)撞分析影響較大。真實(shí)鳥(niǎo)體形狀復(fù)雜，不同鳥(niǎo)類外形、大小各異，研究中常用簡(jiǎn)化替代模型模擬。4種常見(jiàn)鳥(niǎo)體替代模型為球體、圓柱體、兩端半球-中間圓柱體及橢球體。實(shí)際鳥(niǎo)撞事故中，鳥(niǎo)可能從頭、尾、翼及腹部等不同方向撞擊飛機(jī)，所致影響不盡相同。對(duì)此，建立真實(shí)形狀鳥(niǎo)體撞擊剛性靶有限元模型，研究不同撞擊方向?qū)B(niǎo)撞分析影響；進(jìn)行4種簡(jiǎn)化鳥(niǎo)體替代模型分析，并與真實(shí)形狀鳥(niǎo)體分析結(jié)果對(duì)比，研究替代模型對(duì)真實(shí)形狀不同撞擊方向的適用性。結(jié)果表明，不同撞擊方向?qū)Ψ治鼋Y(jié)果影響較大，腹部沖擊力極值最大，沖擊持時(shí)最短；頭、尾部沖擊極值較小。腹部沖擊動(dòng)能衰減最快、最多，尾、翼部次之，頭部最慢。尾部與兩端半球-中間圓柱體撞擊計(jì)算結(jié)果吻合較好；腹部與橢球體長(zhǎng)邊側(cè)撞擊計(jì)算結(jié)果吻合較好；頭、翼部與替代模型計(jì)算結(jié)果吻合較差。

關(guān)鍵詞：鳥(niǎo)撞；鳥(niǎo)體形狀；真實(shí)鳥(niǎo)體模型；撞擊方向；替代模型

中圖分類號(hào):O347文獻(xiàn)標(biāo)志碼：A

基金項(xiàng)目：中國(guó)民航局科技項(xiàng)目(MHRD20140207)；中央高校項(xiàng)目(ZXp012B004)

收稿日期：2015-03-20修改稿收到日期：2015-05-29

Effects of bird’s striking orientation on bird impact analysis based on a realistic bird model

ZHANGDa-hai1,2,FEIQing-guo1,2,LIUHong-yue1,2(1. Department of Engineering Mechanics, Southeast University, Nanjing 210096, China; 2. Jiangsu Provincial Key Laboratory of Engineering Mechanics, Nanjing 210096, China)

Abstract:Bird shape and strike orientation are two very important factors affecting the accuracy of bird strike analysis. The realistic bird has a complex configuration and there are different bird types with different body shapes, so researchers often use simple configurations to model birds in bird strike events. There are four substitute bird models most-frequently employed in literatures, they are sphere, straight-ended cylinder, hemispherical-ended cylinder and ellipsoid. During a bird strike event, a bird can strike an airplane from its head, tail, abdomen or wings. These orientations have different effects on the response of an airplane part. Firstly, a realistic bird model was established to study the effect of its strike orientation. The results obtained by using the four substitute bird models were compared with those using a realistic bird model in order to determine the best substitute bird model for modeling the force exerted by a real bird with different strike orientations. It was shown that the striking orientation of a bird has a big influence on the impact force; the impact force obtained by striking from abdomen is the largest and the impact duration is the shortest while those obtained by striking from head and tail are relatively less; meanwhile, the kinetic energy loss of bird tends to be most and fastest for striking from abdomen, for strikings from wing and tail they are second, for striking from head it is the smallest; besides, for string from tail, a hemispherical-ended cylinder model has the best results and for striking from abdomen, an ellipsoid model striking from its long side is the best substitute model, while for striking from head and wing, no substitute models are better.

Key words:bird strike; bird shape; realistic bird model; striking orientation; substitute bird models

鳥(niǎo)撞指高速飛行器與鳥(niǎo)類相撞的航空事故[1]， 僅1990～2007年美國(guó)境內(nèi)發(fā)生約82000起，且呈逐年遞增趨勢(shì)[2]。而我國(guó)則由2002年94起增加到2011年1538起[3]，造成巨大人員傷亡及經(jīng)濟(jì)損失。鳥(niǎo)撞事故嚴(yán)重威脅航空安全，因此飛機(jī)結(jié)構(gòu)設(shè)計(jì)中鳥(niǎo)撞問(wèn)題已成必須考慮的重要內(nèi)容之一。

鳥(niǎo)撞問(wèn)題為高度非線性沖擊動(dòng)力學(xué)問(wèn)題，鳥(niǎo)體與靶體間存在極強(qiáng)的耦合作用[4]。早期鳥(niǎo)撞研究主要以試驗(yàn)為主，數(shù)值模擬則采用解耦算法。由于計(jì)算機(jī)及有限元技術(shù)發(fā)展，使鳥(niǎo)撞耦合分析更真實(shí)模擬整個(gè)撞擊過(guò)程成為可能。

由于鳥(niǎo)類身體結(jié)構(gòu)復(fù)雜，研究中常用簡(jiǎn)化的替代模型模擬，主要有球體、圓柱體、兩端半球-中間圓柱體及橢球體4種，見(jiàn)圖1。

圖1　四種簡(jiǎn)化鳥(niǎo)體替代模型 Fig.1 Four substitute bird models

對(duì)鳥(niǎo)體替代模型，Richard[5]認(rèn)為替代鳥(niǎo)體應(yīng)標(biāo)準(zhǔn)化，并推薦使用圖1的圓柱體、兩端半球-中間圓柱體及橢球體模擬鳥(niǎo)體。國(guó)標(biāo)[6]規(guī)定鳥(niǎo)撞試驗(yàn)的鳥(niǎo)彈外形為圓柱體，長(zhǎng)徑比為2∶1。Meguid等[7-8]對(duì)文獻(xiàn)[5]的3種形狀鳥(niǎo)體撞擊剛性靶及風(fēng)機(jī)葉片進(jìn)行數(shù)值計(jì)算影響研究，認(rèn)為鳥(niǎo)體、靶體初始接觸面積對(duì)沖擊力影響較大，鳥(niǎo)體長(zhǎng)徑比對(duì)沖擊壓力影響不大。

圖2為一架客機(jī)遭遇鳥(niǎo)群“襲擊”示意圖。實(shí)際鳥(niǎo)撞事故中，因鳥(niǎo)類飛行軌跡具有隨機(jī)性，鳥(niǎo)體可能從頭、尾、翼、腹部等不同方向撞擊飛機(jī)，所致影響不相同。由于鳥(niǎo)體結(jié)構(gòu)的復(fù)雜性，簡(jiǎn)化鳥(niǎo)體替代模型已遠(yuǎn)不能滿足研究需要，故嘗試建立真實(shí)形狀的鳥(niǎo)體進(jìn)行研究。McCallum等[9-10]用加拿大天鵝實(shí)際物理尺寸建立鳥(niǎo)體研究表明，含頭、長(zhǎng)頸的鳥(niǎo)體模型能正確預(yù)測(cè)沖擊載荷及靶體破壞，且計(jì)算結(jié)果沖擊時(shí)程更長(zhǎng)，沖擊力極值更大。Lakshmi[11]通過(guò)建立真實(shí)鳥(niǎo)體模型，對(duì)不同部位取不同密度研究表明，分析小鳥(niǎo)撞問(wèn)題時(shí)可用傳統(tǒng)的簡(jiǎn)化鳥(niǎo)體軀干模型，大鳥(niǎo)撞時(shí)鳥(niǎo)體中其它部位不可忽視。Hedayati等[12-13]建立幾何形狀更準(zhǔn)確精細(xì)的鳥(niǎo)體模型，研究真實(shí)形狀鳥(niǎo)體從各部位對(duì)結(jié)構(gòu)沖擊的靶體損傷，并將真實(shí)形狀鳥(niǎo)體計(jì)算結(jié)果與替代模型計(jì)算結(jié)果及實(shí)驗(yàn)結(jié)果對(duì)比表明，真實(shí)鳥(niǎo)體模型計(jì)算結(jié)果與理論計(jì)算值差別較大，與實(shí)驗(yàn)值吻合較好。

圖2　客機(jī)鳥(niǎo)撞事故 Fig.2 A flock of birds striking an airplane

雖鳥(niǎo)類不同，體型、大小相差較大，但其結(jié)構(gòu)類似，均由橢球狀軀干、細(xì)長(zhǎng)頸、一對(duì)翅膀、遠(yuǎn)小于軀干的頭部及附屬結(jié)構(gòu)組成。本文用真實(shí)帆背潛鴨幾何尺寸建立三維鳥(niǎo)體模型，利用非線性有限元數(shù)值模擬方法建立耦合鳥(niǎo)撞分析模型。分析鳥(niǎo)體從頭、尾、翼、腹部4個(gè)不同方向撞擊靶體，研究撞擊方向影響；并與4種簡(jiǎn)化的替代模型計(jì)算結(jié)果對(duì)比，研究簡(jiǎn)化模型在替代真實(shí)形狀鳥(niǎo)體不同方向撞擊的適用性，分析結(jié)果可為相似鳥(niǎo)體撞擊提供參考。

1鳥(niǎo)體模型

1.1真實(shí)鳥(niǎo)體幾何模型

一架波音767起飛階段遭遇一群帆背潛鴨“襲擊”，造成一臺(tái)發(fā)動(dòng)機(jī)起火，風(fēng)擋、剎車系統(tǒng)不同程度受損。帆背潛鴨大小、質(zhì)量與國(guó)標(biāo)規(guī)定的鳥(niǎo)體類似，本文真實(shí)鳥(niǎo)體模型用帆背潛鴨幾何尺寸建立，三維模型見(jiàn)圖3，幾何尺寸為：鳥(niǎo)體長(zhǎng)290 mm，翼展450 mm，質(zhì)量1.484 kg。

圖3　帆背潛鴨及三維模型 Fig.3 A flying canvasback and its 3D model

1.2替代鳥(niǎo)體幾何模型

建立4種簡(jiǎn)化鳥(niǎo)體替代模型，幾何尺寸見(jiàn)表1。按國(guó)標(biāo)規(guī)定，圓柱體、兩端半球-中間圓柱體及橢球體長(zhǎng)徑比為2∶1。

表1　四種鳥(niǎo)體替代模型幾何尺寸

1.3鳥(niǎo)體材料模型

真實(shí)鳥(niǎo)體結(jié)構(gòu)較復(fù)雜，其力學(xué)特性不均勻、不相似，較難按真實(shí)本構(gòu)進(jìn)行數(shù)值模擬，多采用不同方法近似，如用彈塑性本構(gòu)[14-16]及流體本構(gòu)進(jìn)行模擬[17-19]。研究表明，高速撞擊時(shí)鳥(niǎo)體力學(xué)特性不均勻性、不相似性可忽略；撞擊速度大于100 m/s時(shí)鳥(niǎo)體可視為均勻流體射流撞擊靶體[20-21]。均勻流體材料本構(gòu)方程可表示為

(1)

狀態(tài)方程能成功描述壓縮的流體力學(xué)特性，一般形式為

P=f(ρ,Em)

(2)

(3)

式中：PH為Hugoniot壓力；Γ0為材料常數(shù)，稱為Grüneisen比；ρ0為參考密度；EH為Hugoniot能。

聯(lián)立消去Γ、EH，即

(4)

對(duì)線性US-UP形式，有

(5)

式中：co為參考聲速；s為US-UP曲線斜率。本文狀態(tài)方程中ρ0=980 kg/m3，c0=1 483 m/s，Γ0=0，s=0。

1.4鳥(niǎo)體模型及數(shù)值模擬方法驗(yàn)證

為驗(yàn)證本文鳥(niǎo)體模型及數(shù)值模擬方法的準(zhǔn)確性，對(duì)文獻(xiàn)[7]中鳥(niǎo)撞剛性靶進(jìn)行數(shù)值計(jì)算。鳥(niǎo)體模型采用兩端半球-中間圓柱體模擬，撞擊速度200 m/s。圖4為鳥(niǎo)撞分析中鳥(niǎo)體變形。由圖4看出，鳥(niǎo)體表現(xiàn)出明顯的流體特性，隨撞擊繼續(xù)，鳥(niǎo)體逐步向外擴(kuò)展，動(dòng)能逐步轉(zhuǎn)化為熱能及靶體內(nèi)能。

圖4　兩端半球-中間圓柱型鳥(niǎo)體撞擊剛性靶 Fig.4 A hemispherical-ended cylinder substitute bird model striking a rigid target

將沖擊力計(jì)算結(jié)果與仿真結(jié)果及實(shí)驗(yàn)結(jié)果[22]對(duì)比，見(jiàn)圖5。由圖5看出，本文計(jì)算值與實(shí)驗(yàn)值及文獻(xiàn)中計(jì)算值均吻合較好，從而驗(yàn)證所用鳥(niǎo)體模型及數(shù)值模擬方法的有效性。

圖5　替代鳥(niǎo)體模型沖擊剛性靶歸一化沖擊力隨其時(shí)間變化 Fig.5 Variation of normalized impact force versus normalized time for a hemispherical-ended cylindrical bird striking a rigid target

2鳥(niǎo)撞方向影響

鳥(niǎo)類飛行軌跡隨機(jī)性較大，鳥(niǎo)體可能從任何方向撞擊飛機(jī)，而鳥(niǎo)體結(jié)構(gòu)不對(duì)稱性導(dǎo)致不同方向撞擊靶體結(jié)果亦不相同。本文重點(diǎn)分析鳥(niǎo)體從頭、尾、翼及腹部4個(gè)不同方向撞擊靶體，研究鳥(niǎo)撞方向?qū)Ψ治鲇绊?。圖6為鳥(niǎo)體從4個(gè)不同方向撞擊剛性靶體示意圖。

圖6　鳥(niǎo)體從各方向撞擊靶體 Fig.6 Deformation of birds striking the target form different sides

圖7為鳥(niǎo)體從不同方向、部位撞擊剛性靶體所得沖擊力及鳥(niǎo)體動(dòng)能隨時(shí)間變化曲線。由于計(jì)算采用特定鳥(niǎo)體(帆背潛鴨)幾何參數(shù)，主要對(duì)結(jié)果進(jìn)行定性分析，因此可擴(kuò)展到其它鳥(niǎo)體撞擊，對(duì)沖擊力及時(shí)間進(jìn)行歸一化處理。由圖7看出，①較圖5替代鳥(niǎo)體模型結(jié)果，真實(shí)形狀鳥(niǎo)體所得沖擊力出現(xiàn)明顯次波峰。因鳥(niǎo)體與靶體作用時(shí)前端單元迅速減速產(chǎn)生峰值，此時(shí)后端單元并未減速。隨撞擊后端單元及靶體繼續(xù)作用，產(chǎn)生新峰值。某種程度上可認(rèn)為沖擊力大小、單位時(shí)間與靶體產(chǎn)生作用的鳥(niǎo)體單元成一定正比例關(guān)系。因真實(shí)形狀鳥(niǎo)體結(jié)構(gòu)不規(guī)則性，會(huì)產(chǎn)生明顯次波峰。②撞擊方向?qū)_擊力、鳥(niǎo)體動(dòng)能變化均影響較大，不同方向沖擊力大小、極值出現(xiàn)時(shí)間、撞擊持時(shí)及鳥(niǎo)體動(dòng)能衰減均有較大差別。腹部撞擊沖擊力極值最大，翼部次之，為腹部撞擊的65.8%，頭、尾部最小，分別為腹部的36.8%、34.2%。腹部撞擊持時(shí)遠(yuǎn)小于其它三方向計(jì)算結(jié)果，因真實(shí)鳥(niǎo)體飛行狀態(tài)下腹部方向尺寸遠(yuǎn)小于鳥(niǎo)體長(zhǎng)及翼展。③腹、尾部撞擊力極值出現(xiàn)時(shí)間最早，腹部撞擊力極值出現(xiàn)后迅速衰減，尾部則持時(shí)較長(zhǎng)；翼部初始時(shí)較小，中間出現(xiàn)峰值后慢慢衰減，因翼部撞擊時(shí)開(kāi)始、結(jié)束階段為薄弱的鳥(niǎo)翼與靶體作用，中間為鳥(niǎo)體軀干與靶體作用；頭部撞擊力極值出現(xiàn)時(shí)間最晚，因頭部撞擊時(shí)先與靶體作用的為相對(duì)薄弱的頭、頸部，最后為軀干。④腹部撞擊鳥(niǎo)體動(dòng)能衰減最快、最多，衰減約70%，因鳥(niǎo)體質(zhì)量主要集中在軀干部分，腹部撞擊時(shí)鳥(niǎo)體軀干先接觸靶體，鳥(niǎo)體動(dòng)能會(huì)迅速衰減；而為尾、翼、頭部沖擊鳥(niǎo)體衰減最緩、最少，衰減約29%。翼部撞擊開(kāi)始時(shí)動(dòng)能衰減較緩慢，軀干與靶體作用時(shí)衰減急劇加速，鳥(niǎo)體動(dòng)能-時(shí)間曲線出現(xiàn)明顯拐點(diǎn)。

圖7 鳥(niǎo)體從不同方向撞擊靶體歸一化沖擊力及鳥(niǎo)體動(dòng)能隨歸一化時(shí)間變化Fig.7Variationofnormalizedimpactforceandkineticenergyofbirdversusnormalizedtimeforahemispherical-endedcylindricalbirdstrikingarigidtarget圖8 真實(shí)形狀鳥(niǎo)體頭部撞擊與4種替代模型撞擊歸一化沖擊力及鳥(niǎo)體動(dòng)能隨歸一化時(shí)間變化Fig.8Variationofnormalizedimpactforceandkineticenergyofbirdversusnormalizedtimefordifferentbirdsubstitutesandtherealisticbirdstrikingfromheadside圖9 真實(shí)形狀鳥(niǎo)體尾部撞擊與4種替代模型撞擊歸一化沖擊力及鳥(niǎo)體動(dòng)能隨歸一化時(shí)間變化Fig.9Variationofnormalizedimpactforceandkineticenergyofbirdversusnormalizedtimefordifferentbirdsubstitutesandtherealisticbirdstrikingfromtailside

3鳥(niǎo)體替代模型分析

用1.2節(jié)建立的鳥(niǎo)體替代模型進(jìn)行鳥(niǎo)撞分析，將計(jì)算結(jié)果與真實(shí)形狀鳥(niǎo)撞分析結(jié)果對(duì)比，研究替代模型的適用性，選取可替代真實(shí)形狀鳥(niǎo)體不同撞擊方向的簡(jiǎn)化模型。

3.1頭部撞擊

將真實(shí)形狀鳥(niǎo)體從頭部撞擊所得沖擊力及鳥(niǎo)體動(dòng)能與替代模型對(duì)比，對(duì)結(jié)果進(jìn)行歸一化處理，見(jiàn)圖8。由圖8看出，真實(shí)形狀鳥(niǎo)體頭部沖擊力極值出現(xiàn)時(shí)間約在整個(gè)分析的0.9處，而替代鳥(niǎo)體模型全部出現(xiàn)在0.5前；球體沖擊力極值計(jì)算結(jié)果與頭部最接近，球體沖擊力極值為真實(shí)鳥(niǎo)體頭部撞擊的1.09倍。由于鳥(niǎo)體長(zhǎng)遠(yuǎn)大于替代模型長(zhǎng)度且質(zhì)量主要集中在中后段軀干部，故頭部撞擊時(shí)鳥(niǎo)體動(dòng)能衰減緩慢，只衰減約29%，而替代模型衰減均在50%以上。

綜合考慮沖擊力與鳥(niǎo)體動(dòng)能衰減，4種替代模型均不能較好與真實(shí)鳥(niǎo)體頭部撞擊吻合，若只關(guān)注沖擊力極值大小，可選球體進(jìn)行替代。

3.2尾部撞擊

將真實(shí)形狀鳥(niǎo)體從尾部撞擊所得沖擊力及鳥(niǎo)體動(dòng)能與替代模型進(jìn)行對(duì)比，并對(duì)結(jié)果歸一化處理，見(jiàn)圖9。由圖9看出，尾部沖擊力極值與兩端半球-中間圓柱體計(jì)算結(jié)果吻合較好，后者沖擊力極值為真實(shí)鳥(niǎo)體尾部撞擊的61.5%，兩者出現(xiàn)時(shí)間均在分析的0.1前；球體沖擊力極值為真實(shí)鳥(niǎo)體尾部撞擊的95.5%，但沖擊力持時(shí)遠(yuǎn)小于真實(shí)鳥(niǎo)體尾部。由于鳥(niǎo)體質(zhì)量主要集中在軀干部，尾部沖擊時(shí)動(dòng)能衰減較快，衰減約49%，與除球體外的3替代模型結(jié)果近似。

綜合考慮沖擊力與鳥(niǎo)體動(dòng)能衰減，兩端半球-中間圓柱體與真實(shí)形狀鳥(niǎo)體尾部撞擊吻合最好，為最優(yōu)替代模型選擇。

3.3翼部撞擊

將真實(shí)形狀鳥(niǎo)體翼部撞擊所得沖擊力及鳥(niǎo)體動(dòng)能與替代模型對(duì)比，對(duì)結(jié)果歸一化處理，見(jiàn)圖10。由圖10看出，真實(shí)形狀鳥(niǎo)體翼部撞擊沖擊力極值出現(xiàn)時(shí)間約在整個(gè)撞擊分析的0.58處，最接近的鳥(niǎo)體替代模型為橢球體(0.45)，但橢球體所得沖擊力極值僅為翼部撞擊的32.5%；圓柱體沖擊力極值為翼部撞擊的1.13倍，但出現(xiàn)時(shí)間相差較大。由于鳥(niǎo)體質(zhì)量主要集中在軀干部，翼部撞擊開(kāi)始時(shí)薄弱的鳥(niǎo)翼先與靶體作用，動(dòng)能衰減緩慢，軀干與靶體作用時(shí)衰減急劇加速，鳥(niǎo)體動(dòng)能-時(shí)間曲線出現(xiàn)明顯拐點(diǎn)，動(dòng)能最終衰減約49%。

綜合考慮沖擊力與鳥(niǎo)體動(dòng)能衰減，4種替代模型均不能較好與真實(shí)鳥(niǎo)體翼部撞擊吻合。如只關(guān)注沖擊力極值大小，可選圓柱體替代；若只關(guān)注沖擊力極值出現(xiàn)時(shí)間，可選橢球體模擬；若只關(guān)注鳥(niǎo)體最終衰減，可選除球體外的3種替代模型進(jìn)行模擬。

圖10　真實(shí)形狀鳥(niǎo)體翼部撞擊與4種替代模型撞擊歸一化 沖擊力及鳥(niǎo)體動(dòng)能隨歸一化時(shí)間變化 Fig.10 Variation of normalized impact force and kinetic energy of bird versus normalized time for different bird substitutes and the realistic bird striking from wing side

3.4腹部撞擊

圖11　真實(shí)形狀鳥(niǎo)體腹部與4種替代鳥(niǎo)體模型撞擊 歸一化沖擊力及鳥(niǎo)體動(dòng)能隨歸一化時(shí)間變化 Fig.11 Variation of normalized impact force and kinetic energy of bird versus normalized time for different bird substitutes and the realistic bird striking from abdomen side

將真實(shí)形狀鳥(niǎo)體從腹部撞擊所得沖擊力及鳥(niǎo)體動(dòng)能與替代模型對(duì)比，對(duì)結(jié)果歸一化處理，見(jiàn)圖11。由圖11看出，真實(shí)形狀鳥(niǎo)體腹部撞擊力時(shí)程與鳥(niǎo)體動(dòng)能衰減時(shí)程均遠(yuǎn)小于4種替代模型計(jì)算值，原因?yàn)檎鎸?shí)鳥(niǎo)體腹部撞擊時(shí)腹部方向尺寸遠(yuǎn)小于鳥(niǎo)體長(zhǎng)及翼展。

進(jìn)一步分析發(fā)現(xiàn)，真實(shí)鳥(niǎo)體腹部撞擊作用時(shí)間約為替代模型的一半，而替代模型長(zhǎng)徑比為2∶1，因此對(duì)4種鳥(niǎo)體替代模型長(zhǎng)邊側(cè)撞擊進(jìn)行分析，橢球體替代模型撞擊模型見(jiàn)圖12。

圖12　橢球體型鳥(niǎo)體替代模型撞擊剛性靶 Fig.12 An ellipsoid substitute bird model striking a rigid target

將真實(shí)形狀鳥(niǎo)體從腹部撞擊靶體所得沖擊力及鳥(niǎo)體動(dòng)能與替代模型長(zhǎng)邊側(cè)撞擊計(jì)算結(jié)果進(jìn)行對(duì)比，結(jié)果進(jìn)行歸一化處理，見(jiàn)圖13。由圖13看出，腹部撞擊力與橢球體計(jì)算結(jié)果吻合最好，橢球體沖擊力極值為真實(shí)鳥(niǎo)體腹部撞擊的89.4%，兩者出現(xiàn)時(shí)間均在撞擊分析的0.1前，且均在撞擊分析總時(shí)間0.8處基本衰減為0。鳥(niǎo)體動(dòng)能衰減與橢球體結(jié)果吻合最好。

圖13　真實(shí)形狀鳥(niǎo)體腹部撞擊與4種替代鳥(niǎo)體模型長(zhǎng)邊側(cè) 撞擊歸一化沖擊力及鳥(niǎo)體動(dòng)能隨歸一化時(shí)間變化 Fig.13 Variation of normalized impact force and kinetic energy of bird versus normalized time for different bird substitutes’ long sides and the realistic bird striking from abdomen side

綜合考慮沖擊力與鳥(niǎo)體動(dòng)能衰減，橢球體長(zhǎng)邊側(cè)撞擊計(jì)算結(jié)果與真實(shí)形狀鳥(niǎo)體尾部撞擊吻合最好，為最優(yōu)替代模型選擇。

4結(jié)論

通過(guò)真實(shí)帆背潛鴨幾何尺寸建立鳥(niǎo)體模型，研究其不同撞擊方向?qū)Ψ治鲇绊?；并將?jì)算結(jié)果與4種常見(jiàn)鳥(niǎo)體替代模型計(jì)算結(jié)果對(duì)比，研究替代模型的適用性，結(jié)論如下：

(1)鳥(niǎo)撞方向?qū)_擊力影響較大，不同撞擊方向沖擊力大小、極值出現(xiàn)時(shí)間及沖擊力持時(shí)均差別較大。

(2)腹部撞擊沖擊力持時(shí)最短，沖擊力極值遠(yuǎn)大于其它方向，且撞擊鳥(niǎo)體動(dòng)能衰減最快、最多，危害最嚴(yán)重。

(3)尾部撞擊力與兩端半球-中間圓柱體替代模型吻合較好，腹部與橢球體替代模型長(zhǎng)邊側(cè)撞擊吻合較好；頭、翼部撞擊力與4種鳥(niǎo)體替代模型吻合不好。

參考文獻(xiàn)

[1]劉軍,李玉龍,郭偉國(guó),等.鳥(niǎo)體本構(gòu)模型參數(shù)反演I：鳥(niǎo)撞平板試驗(yàn)研究[J].航空學(xué)報(bào),2011, 32(5):802-811.

LIU Jun, LI Yu-long, GUO Wei-guo, et al. Parameters inversion on bird constitutive model part I: study on experiment of bird striking on plate[J]. Acta Aeronautica et Astronautica Sinica,2011, 32(5): 802-811.

[2]Dolbeer R A, Wright S E. Wildlife strikes to civil aircraft in the united states 1990-2007[R]. FAA National Wildlife Strike Database, Serial Report, 2008.

[3]中國(guó)民用航空局機(jī)場(chǎng)司. 2011年度中國(guó)民航鳥(niǎo)擊航空器事件分析報(bào)告[R].中國(guó)民航科學(xué)技術(shù)研究院,2012.

[4]劉富,張嘉振,童明波,等.2024-T3鋁合金動(dòng)力學(xué)實(shí)驗(yàn)及其平板鳥(niǎo)撞動(dòng)態(tài)響應(yīng)分析[J].振動(dòng)與沖擊,2014,33(4):113-118.

LIU Fu, ZHANG Jia-zhen, TONG Ming-bo, et al. Dynamic tests and bird impact dynamic response analysis for a 2024-T3 aluminum alloy plate [J]. Journal of Vibration and Shock, 2014,33(4):113-118.

[5]Richard B. The development of a substitute artificial bird by the International Bird Strike Research Group for use in aircraft component testing[C]//International Bird Strike Committee, ISBC25/WP-IE3,Amsterdam, 2000.

[6]GJB 2464A-2007,飛機(jī)透明件鳥(niǎo)撞試驗(yàn)方法[S].

[7]Meguid S A, Mao R H, Ng T Y. FE analysis of geometry effects of an artificial bird striking an aero engine fan blade [J]. International Journal of Impact Engineering, 2008,35:487-497.

[8]Mao R H, Meguid S A, Ng T Y. Transient three dimensional finite element analysis of a bird striking a fan blade[J]. International Journal of Mechanics and Materials in Design, 2008, 4:79-96.

[9]McCallum S C, Constantinou C. The influence of bird-shape in bird-strike analysis[C]//5thEuropean LS-DYNA Users Conference,2005.

[10]McCallum S, Shoji H, Akiyama H. Development of an advanced multi-material bird-strike model using the smoothed particle hydrodynamics method[J]. International Journal of Crashworthiness,2013, 18(6): 579-596.

[11]Lakshmi S N. Models and methods for bird strike load predictions [D].Wichita : Wichita State University,2007.

[12]Hedayati R, Ziaei-Rad S. A new bird model and the effect of bird geometry in impacts from various orientatons[J]. Aerospace Science and Technology,2013(28):9-20.

[13]Hedayati R, Sadighi M,Mohammadi-Aghdam M. On the difference of pressure readings from the numerical, experimental and theoretical results in different bird strike studies[J]. Aerospace Science and Technology,2014(32): 260-266.

[14]Wang F S, Yue Z F. Numerical simulation of damage and failure in aircraft windshield structure against bird strike[J]. Materials and Desigh,2010, 31:687-695.

[15]王富生,李立州,王新軍,等.鳥(niǎo)體材料參數(shù)的一種反演方法[J]. 航空學(xué)報(bào),2007, 28(2):344-347.

WANG Fu-sheng, LI Li-zhou, WANG Xin-jun, et al. A method to identify bird’s material parameters[J]. Acta Aeronautica et Astronautica Sinica,2007, 28(2):344-347.

[16]杜龍.基于歐拉-拉格朗日方法的復(fù)合材料機(jī)翼前緣鳥(niǎo)撞模擬[J].振動(dòng)與沖擊,2012, 31(7):137-141.

DU Long. Finite element analysis of bird striking on a composite wing leading edge based on eulerian-lagrangian method[J].Journal of Vibration and Shock,2012,31(7):137-141.

[17]Dar U A, Zhang W, Xu Y. FEanalysis of dynamic response of aircraft windshield against bird impact[J].International Journal of Aerospace Engineering,2013.

[18]林長(zhǎng)亮,王益鋒,王浩文,等.直升機(jī)旋翼槳葉鳥(niǎo)撞動(dòng)態(tài)響應(yīng)計(jì)算[J].振動(dòng)與沖擊,2013, 32(10):62-68.

LIN Chang-liang,WANG Yi-feng,WANG Hao-wen, et al. Computation of dynamic response of helicopter blade under bird impact[J].Journal of Vibration and Shock,2013,32(10): 62-68.

[19]Smojver I, Ivancevic D. Bird strike damage analysis in aircraft structure using abaqus/explict and coupled eulerian lagrangian approach[J].Composites Science and Technology, 2011, 71:489-498.

[20]Liu Jun, Li Yu-long, Gao Xiao-sheng. Bird strike on a flat plate: experiments and numerical simulations[J].International Journal of Impact Engineering 2014, 70:21-37.

[21]Lavoie M A,Gakwaya A, Ensan M N, et al. Bird’s substitute tests results and evaluation of available numerical methods [J].International Journal of Impact Engineering,2009,36: 1276-1287.

[22]Willows M, Driffill B. GARTEUR (group for aeronautical research and technology in europe) bird strike group, round robin work package: rigid wall phase1 and task 1[M]. Hants: DERA Farnborough, 1999.

第一作者解江男，博士，講師，1982年1月生