復(fù)合材料層合板鳥彈撞擊損傷分析

王 浪，孫有朝，周亞?wèn)|，楊 坤，曾海軍，陳 川

（1.南京航空航天大學(xué)民航學(xué)院，南京 210016；2.中國(guó)航發(fā)商用航空發(fā)動(dòng)機(jī)有限責(zé)任公司，上海 200241；3.中國(guó)南方航空公司股份有限公司，廣州 510470）

0 引言

隨著中國(guó)經(jīng)濟(jì)和科技的快速發(fā)展，飛機(jī)的數(shù)量和班次都在不斷增加，因鳥撞而發(fā)生的飛機(jī)安全事故也越來(lái)越多，據(jù)統(tǒng)計(jì)，在飛機(jī)的地面安全事故中，鳥類碰撞事故占21%。其中大量事故的起因是鳥類被吸入航空發(fā)動(dòng)機(jī)，造成發(fā)動(dòng)機(jī)葉片變形，進(jìn)而導(dǎo)致發(fā)動(dòng)機(jī)推力降低、葉片斷裂等嚴(yán)重后果，嚴(yán)重威脅民用和軍用航空的安全，引起國(guó)際社會(huì)的高度重視。因此，設(shè)計(jì)和研究航空發(fā)動(dòng)機(jī)葉片防鳥撞擊保護(hù)裝置對(duì)于提高航空飛行的安全性、減少財(cái)產(chǎn)損失、保障乘客生命安全具有極其重要的意義。

為避免鳥撞飛機(jī)事件的發(fā)生，提高航空器結(jié)構(gòu)的抗鳥撞擊能力，迫切需要通過(guò)開(kāi)發(fā)高性能材料和更好的耐撞性設(shè)計(jì)來(lái)提高飛機(jī)風(fēng)扇葉片結(jié)構(gòu)的抗鳥撞性能。一般來(lái)說(shuō)，鳥撞是持續(xù)時(shí)間短（一般在毫秒范圍內(nèi))的高強(qiáng)度載荷加載事件，具有復(fù)雜的耦合相互作用。在鳥彈的軟體撞擊作用下，鳥彈的幾何效應(yīng)、材料模型，以及初始速度和落角等關(guān)鍵參數(shù)都會(huì)影響風(fēng)扇葉片結(jié)構(gòu)的變形和損傷機(jī)理。為了解決復(fù)合材料風(fēng)扇葉片在鳥彈高速?zèng)_擊下變形和損壞的問(wèn)題，Hedayati等研究了平板不同厚度、鋪層方式對(duì)平板鳥撞總能量的影響；Di等研究了鳥撞風(fēng)扇葉片層合板厚度對(duì)飛機(jī)的影響，計(jì)算出鳥撞層合板的臨界厚度；Zhou等在不改變整體幾何形狀的情況下，考慮了不同材料體系（不同鋪層順序和總厚度)的鳥擊耐撞性；Kim等、Siddens等分別采用任意拉格朗日-歐拉（Arbitrary Lagrange-Euler，ALE)方法和無(wú)網(wǎng)格拉格朗日粒子方法研究了被鳥擊中的單個(gè)風(fēng)扇葉片的初始退化，并評(píng)估了對(duì)其他發(fā)動(dòng)機(jī)部件的后續(xù)損壞；Heimbs等介紹了在軟體沖擊下承受拉伸和壓縮預(yù)載荷的層壓復(fù)合板的動(dòng)態(tài)行為，發(fā)現(xiàn)拉伸和壓縮預(yù)加載會(huì)對(duì)損壞模式產(chǎn)生影響。上述研究表明，影響疊層復(fù)合材料鳥撞損傷的關(guān)鍵因素包括層合板的總厚度、鋪層角度、轉(zhuǎn)速和撞擊位置等。

本文進(jìn)行了鳥彈在4個(gè)角度下撞擊復(fù)合材料層合板的試驗(yàn)，并進(jìn)行了數(shù)值模擬，建立了復(fù)合材料層合板靶板，分析了復(fù)合材料層合板與彈體發(fā)生碰撞后的受力情況，建立風(fēng)扇真實(shí)模型模擬鳥彈撞擊風(fēng)扇。

1 復(fù)合材料層合板鳥撞數(shù)值模擬

1.1 鳥彈幾何模型的建立及網(wǎng)格劃分

利用Ironpython語(yǔ)言對(duì)彈體尺寸以及復(fù)合材料層合板靶板尺寸進(jìn)行參數(shù)化，定義鳥彈以及復(fù)合材料層合板的尺寸與實(shí)際試驗(yàn)狀態(tài)下的尺寸一致，程序化方法有利于模型生成，同時(shí)縮短工作時(shí)間。將鳥彈形狀簡(jiǎn)化為圓柱體，鳥彈幾何尺寸如圖1所示，靶板為長(zhǎng)方形。對(duì)鳥彈模型進(jìn)行全六面體網(wǎng)格劃分，網(wǎng)格數(shù)量為9483，鳥彈模型網(wǎng)格劃分如圖2所示。

1.2 復(fù)合材料層合板幾何模型的建立

對(duì)復(fù)合材料層合板殼體進(jìn)行建模，將網(wǎng)格模型生成K文件導(dǎo)入Isprepost進(jìn)行TSHELL單元的鋪層，為便于定義每層材料，分別對(duì)每層的積分點(diǎn)進(jìn)行層數(shù)定義。從板上方到板下方共分5層，每層的鋪層角度分別為0°、45°、90°、-45°、0°，網(wǎng)格數(shù)量為3960。復(fù)合材料層合板尺寸、模型網(wǎng)格如圖3、4所示。

圖1 鳥彈幾何尺寸

圖2 鳥彈模型網(wǎng)格

圖3 復(fù)合材料層合板尺寸

圖4 復(fù)合材料層合板模型網(wǎng)格

圖5 在不同入射角度下鳥彈撞擊的組裝彈體-靶板模型

鳥彈在不同入射角度下撞擊靶板時(shí)的組裝彈體-靶板模型如圖5所示。圖中紅色矩形為鳥彈，綠色矩形為靶板，鳥彈入射角度為靶板與重力方向夾角，分別為0°、15°、30°、45°，鳥彈撞擊在矩形靶板一條長(zhǎng)邊的中部。

考慮到邊界效應(yīng)的影響，試驗(yàn)采用全模型進(jìn)行仿真計(jì)算。

1.3 仿真計(jì)算約束和材料模型的確立

1.3.1 仿真計(jì)算約束條件的設(shè)置

該仿真計(jì)算模型采用ANSYS/LSDYNA中的拉格朗日算法，彈體與復(fù)合材料層合板之間的接觸為侵蝕接 觸（*CONTACT_ERODING_SURFACE_TO_SURFACE，簡(jiǎn)稱為ESTS），建立的全模型無(wú)需設(shè)置對(duì)稱邊界約束，同時(shí)根據(jù)試驗(yàn)條件，在目標(biāo)板件單短邊施加固定邊界約束，彈體速度設(shè)置為140 m/s。具體約束條件如圖6所示，在試驗(yàn)中平板采用一端懸臂的方式夾持，如圖7所示。

圖6 目標(biāo)板件的約束條件

圖7 試驗(yàn)平板一端懸臂夾持

1.3.2 仿真計(jì)算材料模型的設(shè)置

數(shù)值模擬面臨的挑戰(zhàn)之一是建立一個(gè)合適的彈體和復(fù)合材料層合板材料模型，本文考慮到試驗(yàn)中鳥彈出現(xiàn)的流體分散現(xiàn)象，采用*MAT_ELASTIC_PLASTIC_HYDRO材料模型。

靶板模擬采用LS-DYNA中的COMPOSITE損傷材料模型，該模型適合描述復(fù)合材料等脆性材料在大應(yīng)變、高應(yīng)變率和高壓下的動(dòng)態(tài)響應(yīng)，鳥彈和靶板材料模型參數(shù)分別見(jiàn)表1、2，其中靶板、方向分別為靶板長(zhǎng)邊、短邊方向，方向?yàn)榕c靶板垂直的方向。

表1 鳥彈材料模型參數(shù)

表2 靶板材料模型參數(shù)

2 數(shù)值模擬結(jié)果分析

在上述模型的基礎(chǔ)上，對(duì)數(shù)值模擬鳥撞復(fù)合材料層合板損傷情況進(jìn)行了分析。在0°、15°、30°、45°4個(gè)角度下鳥彈與復(fù)合材料層合板接觸的運(yùn)動(dòng)過(guò)程如圖8所示。由于復(fù)合材料層合板角度發(fā)生變化，初始撞擊時(shí)間不相同。在4個(gè)角度下鳥彈的速度時(shí)程曲線如圖9所示。從圖中可見(jiàn)，在4個(gè)角度下復(fù)合材料層合板的鳥彈初始速度均為140 m/s，當(dāng)復(fù)合材料層合板與鳥彈接觸時(shí)，鳥彈的速度開(kāi)始減小，45°復(fù)合材料層合板最先與鳥彈接觸，其次分別是30°、15°、0°復(fù)合材料層合板，而后隨著時(shí)間的延長(zhǎng)，鳥彈速度不斷減小。45°復(fù)合材料層合板的鳥彈速度達(dá)到平衡的時(shí)間最短，但30°復(fù)合材料層合板的殘余速度最小。

圖8 在不同入射角度下鳥彈彈體與復(fù)合材料層合板接觸的運(yùn)動(dòng)過(guò)程

圖9 在不同入射角度下鳥彈的速度時(shí)程曲線

鳥彈沿風(fēng)扇進(jìn)氣方向撞擊風(fēng)扇葉片時(shí)，切鳥角度會(huì)改變，不同角度的撞擊對(duì)葉片造成的損傷也不同。為了分析復(fù)合材料層合板與鳥彈碰撞后的受力情況，對(duì)4個(gè)角度下復(fù)合材料層合板與彈體碰撞進(jìn)行了數(shù)值模擬，不同角度下復(fù)合材料層合板的應(yīng)力如圖10所示。從圖10（a）中可見(jiàn)，0°復(fù)合材料層合板在第1.11 ms時(shí)產(chǎn)生應(yīng)力，在第1.19 ms時(shí)開(kāi)始出現(xiàn)裂紋，裂紋從斷裂位置呈折線形延展，在第1.77 ms時(shí)復(fù)合材料層合板全部斷裂；從圖10（b）中可見(jiàn)，15°復(fù)合材料層合板在第0.93 ms時(shí)產(chǎn)生應(yīng)力，在第1.04 ms時(shí)出現(xiàn)裂紋，裂紋從斷裂位置呈折線形延展，在第1.63 ms時(shí)復(fù)合材料層合板全部斷裂；從圖10（c）中可見(jiàn)，30°復(fù)合材料層合板在第0.76 ms時(shí)產(chǎn)生應(yīng)力，在第0.91 ms時(shí)開(kāi)始出現(xiàn)裂紋，與0°和15°復(fù)合材料層合板不同，裂紋呈直線形延展，在第1.61 ms時(shí)層合板完全斷裂；從圖10（d）中可見(jiàn)，45°復(fù)合材料層合板在第0.61 ms時(shí)產(chǎn)生應(yīng)力，在第0.82 ms時(shí)出現(xiàn)裂紋，與30°復(fù)合材料層合板破壞裂縫相同，裂紋呈直線形延展，在第1.57 ms時(shí)平板完全斷裂?？梢?jiàn)由于角度不同，鳥彈距復(fù)合材料層合板的距離也不同，角度越大，距離越近，鳥彈與復(fù)合材料層合板接觸的時(shí)間越短，但觀察到一個(gè)有趣的現(xiàn)象，角度越大，復(fù)合材料層合板從開(kāi)始產(chǎn)生應(yīng)力到出現(xiàn)裂紋的時(shí)間越長(zhǎng)，說(shuō)明角度越大，復(fù)合材料層合板的抗鳥撞效果越好。綜合比較鳥彈的殘余速度和復(fù)合材料層合板的損壞情況可知，鳥彈入射角度為30°時(shí)復(fù)合材料層合板的抗鳥撞效果最好。

圖10 在不同入射角度下鳥彈撞擊時(shí)復(fù)合材料層合板的應(yīng)力分析

在4個(gè)不同入射角度下鳥彈撞擊復(fù)合材料層合板時(shí)，復(fù)合材料層合板應(yīng)變隨時(shí)間變化的試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果對(duì)比如圖11所示。圖中實(shí)線為仿真結(jié)果曲線，虛線為試驗(yàn)結(jié)果曲線；應(yīng)變仿真結(jié)果最大值與試驗(yàn)結(jié)果最大值的比較見(jiàn)表3。仿真結(jié)果與試驗(yàn)結(jié)果的誤差

從表3中可見(jiàn)，應(yīng)變仿真結(jié)果與試驗(yàn)結(jié)果的誤差在10%以內(nèi)，驗(yàn)證了數(shù)值模型的正確性。

圖11 在不同入射角下鳥彈撞擊時(shí)復(fù)合材料層合板應(yīng)變隨時(shí)間變化的試驗(yàn)與數(shù)值模擬結(jié)果對(duì)比

表3 復(fù)合材料層合板應(yīng)變最大值仿真結(jié)果與試驗(yàn)結(jié)果的比較

為模擬真實(shí)環(huán)境下鳥彈撞擊風(fēng)扇葉片的情況，建立風(fēng)扇葉片的計(jì)算模型，如圖12所示。

風(fēng)扇轉(zhuǎn)速為6000 r/min時(shí)，在不同入射角度（通過(guò)投鳥角度改變）下鳥彈撞擊的速度時(shí)程曲線如圖13所示。鳥彈初始速度均為140 m/s，不同入射角度下鳥彈與風(fēng)扇葉片碰撞的時(shí)刻分別為：0°時(shí)為第0.79 ms，15°時(shí)為第0.62 ms，30°時(shí)為第0.47 ms，45°時(shí)為第0.79 ms。從圖中可見(jiàn)，入射角度越大，鳥彈與風(fēng)扇葉片的接觸時(shí)間就越短，這與試驗(yàn)結(jié)果一致，且角度為30°時(shí)，復(fù)合材料層合板的殘余速度最小，這在工程上十分有利。

圖12 風(fēng)扇葉片計(jì)算模型

圖13 風(fēng)扇轉(zhuǎn)速為6000 r/min時(shí)，在不同入射角度下鳥彈撞擊的速度時(shí)程曲線

在不同入射角度下鳥彈位移隨時(shí)間變化曲線如圖14所示。從圖中可見(jiàn)，各角度下鳥彈撞擊初始時(shí)刻的位移均為0，而后隨著時(shí)間的延長(zhǎng)，各角度下鳥彈位移不斷增大，呈線性增長(zhǎng)趨勢(shì)，可以看到風(fēng)扇葉片位移在0°時(shí)最大，在15°時(shí)次之，在30°和45°時(shí)幾乎相等。

對(duì)鳥彈與風(fēng)扇葉片靶板接觸的單元進(jìn)行應(yīng)力分析，各角度下接觸單元的應(yīng)力最大值如圖15所示。從圖中可見(jiàn)，風(fēng)扇葉片靶板的應(yīng)力在鳥彈入射角度為0°時(shí)最小，隨著入射角度的增大，風(fēng)扇葉片靶板的應(yīng)力也隨之增大，在入射角度為15°和30°時(shí)應(yīng)力值相差不多，在45°時(shí)應(yīng)力最大。綜上所述，在風(fēng)扇葉片靶板為30°時(shí)，鳥彈殘余速度最小，風(fēng)扇葉片靶板的位移最小，承受應(yīng)力較小，靶板的抗鳥撞效果也最好。

3 結(jié)論

圖14 在不同入射角度下鳥彈撞擊時(shí)風(fēng)扇葉片位移隨時(shí)間變化的曲線

圖15 在不同入射角度下風(fēng)扇葉片靶板接觸單元的應(yīng)力最大值

（1）復(fù)合材料層合板鳥撞數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果一致，綜合比較鳥彈的殘余速度和復(fù)合材料層合板的損壞情況，發(fā)現(xiàn)在鳥彈入射角度為30°時(shí)復(fù)合材料層合板的抗鳥撞效果最好。

（2）風(fēng)扇轉(zhuǎn)速為6000 r/min時(shí)，鳥彈入射角度越大，鳥彈與風(fēng)扇葉片接觸的時(shí)間就越短，接觸碰撞后鳥彈的速度減小越多、殘余速度越小。鳥彈入射角度為30°時(shí)，復(fù)合材料層合板的抗鳥撞效果仍然最好。