基于不同構(gòu)型輔助梁的民機尾翼前緣設(shè)計與抗鳥撞性能研究

高 俊， 吳志斌，2， 孔令勇， 李 琪

(1.上海飛機設(shè)計研究院，上海 201210； 2.西北工業(yè)大學(xué) 航空學(xué)院，西安 710072)

鳥撞是飛機飛行中經(jīng)常遇到的沖擊動力學(xué)問題。僅2016年，中國境內(nèi)民航客機就發(fā)生鳥撞事件4 681起，構(gòu)成事故癥候202起，經(jīng)濟損失超過1.1 億元，其中尾翼鳥撞事故癥候率占比超過30%，遠高于其他部位[1]。CCAR25部規(guī)定，民用飛機應(yīng)在與規(guī)定質(zhì)量的鳥相撞后，仍能夠繼續(xù)安全飛行和著陸，其中尾翼結(jié)構(gòu)鳥體質(zhì)量規(guī)定為3.6 kg，其他部位鳥體質(zhì)量規(guī)定為1.8 kg[2]。尾翼鳥撞工況相比其他部位更為嚴酷。

尾翼抗鳥撞性能主要取決于前緣結(jié)構(gòu)設(shè)計[3]。尾翼前緣抗鳥撞設(shè)計有多種形式，包括鋁合金蒙皮及FMLs蒙皮[4-6]、前緣隔板或輔助梁[7-8]、蜂窩夾層[9]、三角板[10]、泡沫[11]等。其中輔助梁適用于前梁站位距離前緣較遠的結(jié)構(gòu)，在實際飛機尾翼中應(yīng)用極為廣泛。

研究飛機結(jié)構(gòu)抗鳥撞問題，通常有數(shù)值計算和試驗兩種方法。其中鳥撞試驗是驗證結(jié)構(gòu)抗鳥撞能力的最直接有效的手段而被廣泛使用。不過試驗方法周期長、成本高，相比而言數(shù)值計算具有迭代速度快、成本低、環(huán)保的優(yōu)勢，可以有效減少試驗數(shù)量并指導(dǎo)試驗件設(shè)計。數(shù)值計算時鳥體的模擬有Lagrange[12-13],Euler[14]，ALE(arbitrary lagrangian eulerian)[15]、SPH(smoothed particle hydrodynamics)[16-18]等多種方法，其中SPH方法是一種基于粒子的無網(wǎng)格方法，該方法把鳥體看作一組相互作用的離散粒子，通過內(nèi)插值方法控制各個SPH單元，是一種具有無網(wǎng)格性、自適應(yīng)性、穩(wěn)定性以及Lagrange性的動力學(xué)方法[19]，克服了Lagrange方法中的網(wǎng)格畸變、Euler方法不適合低速鳥體模擬和ALE方法中的重復(fù)計算問題,具有明顯的優(yōu)勢。

某型民機尾翼原始前緣使用了一種機加工的輔助梁，但是其抗鳥撞性能較差。本文提出了兩種形式的鈑金輔助梁結(jié)構(gòu)作為優(yōu)化設(shè)計方案。通過數(shù)值計算及試驗結(jié)合的方法，評估不同構(gòu)型的抗鳥撞性能，為尾翼前緣設(shè)計提供方案。

1 前緣輔助梁結(jié)構(gòu)設(shè)計

典型帶輔助梁的尾翼前緣在遭受鳥撞時，鳥體穿透前緣后會撞擊輔助梁，輔助梁通過變形、破壞等形式吸收能量，以保護后面的前梁等主承力結(jié)構(gòu)及作動器等關(guān)鍵系統(tǒng)和設(shè)備，因此在抗鳥撞過程中，輔助梁至關(guān)重要。傳統(tǒng)尾翼前緣使用的輔助梁通常是鋁合金機加件，這種形式的輔助梁無需模具，加工周期短、精度高、易裝配。但是需要意識到民用飛機尾翼前緣設(shè)計的關(guān)鍵因素是其抗鳥撞性能。

某型民機原始前緣構(gòu)型如圖1所示，該構(gòu)型使用了機加輔助梁及肋。此構(gòu)型稱之為構(gòu)型1，其輔助梁和肋的材料是7050-T7451，輔助梁的名義厚度為2.54 mm，即輔助梁大部分區(qū)域厚度為2.54 mm，但是加強筋、局部加強區(qū)及R區(qū)附近等區(qū)域厚度可能大于2.54 mm。

圖1 構(gòu)型1前緣結(jié)構(gòu)

構(gòu)型2的前緣結(jié)構(gòu)通過剛度匹配的理念，設(shè)計了鈑金結(jié)構(gòu)形式的輔助梁及肋，其示意圖如圖2所示。輔助梁和肋的材料是2024-T42，輔助梁厚度為2.54 mm。

圖2 構(gòu)型2前緣結(jié)構(gòu)

構(gòu)型3的輔助梁相比構(gòu)型2增加了兩條“波紋加強筋”，此設(shè)計有兩個目的，一是引導(dǎo)輔助梁在鳥撞過程中通過拉伸波紋變形吸收鳥體能量，二是提升腹板縱向剛度，防止局部屈曲。構(gòu)型3的結(jié)構(gòu)示意圖如圖3所示，其輔助梁和肋的材料是2024-T42，輔助梁厚度為2.54 mm。構(gòu)型2和構(gòu)型3的輔助梁截面對比如圖4所示。

圖3 構(gòu)型3前緣結(jié)構(gòu)

圖4 構(gòu)型2和構(gòu)型3輔助梁截面圖對比

3種不同輔助梁所用材料的性能參數(shù)，如表1所示。由表可知，用于機加工件的7050-T7451強度更高，但是用于鈑金件的2024-T42具有更好的延展性。

表1 兩種材料力學(xué)性能表

由表2可知，雖然3種構(gòu)型輔助梁的名義厚度相同，但是構(gòu)型1的機加工件的厚度有不同分區(qū)，且設(shè)計有加強筋條，使其質(zhì)量較構(gòu)型2和構(gòu)型3多2 kg。

表2 輔助梁質(zhì)量對比表

2 鳥撞數(shù)值計算

鳥撞是高度非線性的過程，本文中涉及的結(jié)構(gòu)動力學(xué)數(shù)值計算利用顯示動態(tài)分析軟件PAM-CRASH完成。

2.1 數(shù)值計算模型

2.1.1 鳥體模型

鳥體離散為SPH單元，外形是兩端為半球形的圓柱體，其結(jié)構(gòu)示意圖如圖5所示。

圖5 鳥體SPH模型尺寸示意圖

鳥體的本構(gòu)模型通過Mornaghan EOS定義。

(1)

式中：P0為初始壓力；ρ0為鳥體初始密度，本文取ρ0=9.6×10-7kg/mm3;B為體積模量;λ為狀態(tài)參數(shù)。這兩個參數(shù)通過分析與試驗對標得到，本文取B=0.128 GPa，λ=7.98。

2.1.2 尾翼結(jié)構(gòu)模型

建模時為了合理權(quán)衡計算精度和計算速度的關(guān)系，前緣、輔助梁和肋等鳥撞區(qū)域的建模尺寸選取為7 mm，其他區(qū)域的建模尺寸選取為20 mm。

鋁合金輔助梁采用Belytschko-Wong-Chiang(BWC)殼單元對其進行模擬，其材料采用Johnson-Cook本構(gòu)模型，不考慮材料的應(yīng)變率及溫度效應(yīng)，其材料模型本構(gòu)方程為

(2)

式中:σ為應(yīng)力；εp為塑性應(yīng)變；a,b和n分別為材料的屈服應(yīng)力、硬化系數(shù)、硬化指數(shù)，通過設(shè)置失效應(yīng)變來控制單元的失效刪除。

前緣結(jié)構(gòu)中的金屬玻纖混雜層板結(jié)構(gòu)同樣利用BWC殼單元對其進行模擬。蒙皮中金屬本構(gòu)采用式(2)方程構(gòu)建；對于蒙皮中的玻纖材料，選用單向復(fù)材整體鋪層模型，定義其各個方向力學(xué)參數(shù)。前緣結(jié)構(gòu)中的蜂窩采用縮減積分的六面體的模擬，使用改進的實體單元蜂窩材料，定義其各個方向力學(xué)參數(shù)。

對于結(jié)構(gòu)中的緊固件，采用PLNK單元來模擬，用以定義其與結(jié)構(gòu)間載荷的傳遞及自身的失效。

2.1.3 接觸算法

鳥體和結(jié)構(gòu)體之間發(fā)生沖擊的過程中，必須遵循碰撞方程，在整體坐標系中，碰撞方程[20]可以表示為

[M]{a}+[C]{v}+[K]j5i0abt0b={Fex}

(3)

式中：[M]為質(zhì)量矩陣；{a}為加速度向量;[C]為阻尼矩陣;{v}為速度向量;[K]為剛度矩陣;j5i0abt0b為位移向量;{Fex}為外力向量。

定義[C]{v}+[K]j5i0abt0b={Fin}，則可以得到

[M]{a}={Fex}-{Fin}

(4)

本文采用非對稱主從接觸算法定義鳥體和結(jié)構(gòu)之間的接觸；采用自適應(yīng)接觸算法定義結(jié)構(gòu)之間的自接觸。

2.2 邊界條件及載荷施加

SPH鳥體模型及3種前緣結(jié)構(gòu)的有限元模型分別如圖6～圖9所示，有限元模型規(guī)模見表3。如圖10所示，將前緣及過渡段模型通過夾具固定到承力墻上，并承力墻底部所有節(jié)點固支約束；鳥體質(zhì)量為3.6 kg，速度為186 m/s，方向為逆航向。

圖6 鳥體SPH模型

圖7 構(gòu)型1前緣有限元模型

圖8 構(gòu)型2前緣有限元模型

圖9 構(gòu)型3前緣有限元模型

表3 有限元模型概況

圖10 邊界條件及載荷

2.3 數(shù)值計算

2.3.1 典型分析工況

定義典型的分析工況MID-A、MID-C、MID-AC如圖11所示，其中MID-A點位于前緣邊界和肋中間，MIC-C位于輔助梁與肋連接區(qū)，MID-AC則位于兩點之間。

圖11 鳥撞工況位置

2.3.2 構(gòu)型1分析結(jié)果

構(gòu)型1的前緣結(jié)構(gòu)鳥撞數(shù)值計算結(jié)果中各工況輔助梁的塑性應(yīng)變云圖，如圖12～圖14所示，分析結(jié)果中各工況的輔助梁均發(fā)生穿透性破壞，未穿透區(qū)域變形較小，整體而言該構(gòu)型抗鳥撞性能較差，難以保證輔助梁后方的結(jié)構(gòu)及系統(tǒng)設(shè)備的安全性。

圖12 構(gòu)型1 MID-A工況輔助梁塑性應(yīng)變云圖

圖13 構(gòu)型1中MID-AC工況輔助梁塑性應(yīng)變云圖

圖14 構(gòu)型1中MID-C工況輔助梁塑性應(yīng)變云圖

2.3.3 構(gòu)型2分析結(jié)果

構(gòu)型2的前緣結(jié)構(gòu)鳥撞數(shù)值計算結(jié)果的各工況輔助梁的塑性應(yīng)變云圖，如圖15～圖17所示，即構(gòu)型2的輔助梁在MID-A和MIC-C工況中僅發(fā)生變形，而MID-AC工況中輔助梁產(chǎn)生裂紋，整體而言構(gòu)型2的抗鳥撞性能較構(gòu)型1有明顯改善，這是由于該構(gòu)型鈑金輔助梁以抗鳥撞為主要設(shè)計因素，沒有加強筋條等硬點，避免了剛度突變，且材料延展性好，在發(fā)生鳥撞時，輔助梁在很大跨度范圍內(nèi)發(fā)生變形，吸收鳥撞過程中的能量，避免局部沖擊載荷過大而發(fā)生穿透。MID-AC工況中輔助梁產(chǎn)生裂紋，說明此構(gòu)型的抗鳥撞能力已經(jīng)達到極限。

圖15 構(gòu)型2中MID-A工況輔助梁損塑性應(yīng)變云圖

圖16 構(gòu)型2中MID-AC工況輔助梁應(yīng)變云圖

圖17 構(gòu)型2中MID-C工況輔助梁應(yīng)變云圖

2.3.4 構(gòu)型3分析結(jié)果

構(gòu)型3的前緣結(jié)構(gòu)鳥撞數(shù)值計算結(jié)果的各工況輔助梁的塑性應(yīng)變云圖，如圖18～圖20所示。相比構(gòu)型2，構(gòu)型3增加了兩條“波紋加強筋”，各工況均僅發(fā)生塑性變形，沒有發(fā)生斷裂，整體抗鳥撞性能較構(gòu)型2進一步提升，說明通過引導(dǎo)輔助梁在鳥撞過程中通過拉伸波紋變形吸收鳥體能量的思路是有效的。

圖18 構(gòu)型3中MID-A工況輔助梁應(yīng)變云圖

圖19 構(gòu)型3中MID-AC工況輔助梁損傷云圖

圖20 構(gòu)型3中MID-C工況輔助梁損傷云圖

2.4 分析結(jié)果

表4將2.3節(jié)的分析結(jié)果進行了梳理。由表可知，構(gòu)型2中的鈑金輔助梁的抗鳥撞性能較構(gòu)型1有明顯提升，而構(gòu)型3在鈑金輔助梁腹板上增加了“波紋加強筋”后，能夠進一步提升結(jié)構(gòu)的抗鳥撞能力。

表4 數(shù)值計算結(jié)果匯總

3 鳥撞試驗驗證及試后計算分析

根據(jù)分析結(jié)果，選取構(gòu)型1試驗件1件和構(gòu)型3試驗件2件，開展鳥撞試驗驗證。試驗中構(gòu)型1作為對比基準，主要考察構(gòu)型3的抗鳥撞性能。

鳥撞試驗系統(tǒng)原理圖如圖21所示，系統(tǒng)主要由鳥彈發(fā)射系統(tǒng)、激光測速系統(tǒng)、高速攝像系統(tǒng)、信號采集處理系統(tǒng)組成，鳥體速度由氣罐壓力大小控制。

圖21 鳥撞試驗系統(tǒng)原理圖

試驗件安裝狀態(tài)如圖22所示，試驗件通過過渡段與夾具相連，夾具安裝在承力墻上；在試驗件之后布置一塊鋁合金平板，稱之為“后板”，若試驗時輔助梁發(fā)生穿透，鳥體會撞擊后板使其變形，故可以通過后板的變形情況評估鳥體的剩余能量。試驗前在每件試驗件的輔助梁上粘貼應(yīng)變片，應(yīng)變片粘貼位置如圖23所示，其中LG表示應(yīng)變片,BG代表應(yīng)邊花。

試驗實施的概況如表5所示。試驗完成后，依據(jù)試驗實際測速調(diào)整鳥體速度、對照高速攝像微調(diào)鳥體姿態(tài)，調(diào)整模型后完成了試驗后鳥撞分析。

圖22 試驗件安裝狀態(tài)

圖23 應(yīng)變片在輔助梁上粘貼位置示意圖

表5 鳥撞試驗概況

3.1 構(gòu)型1的MID-A工況試驗及試后分析結(jié)果

構(gòu)型1的MID-A工況試驗結(jié)果，如圖24～圖26所示，試驗后分析結(jié)果如圖27～圖28所示，分析結(jié)果與試驗結(jié)果中前緣和輔助梁的失效形式吻合得較好，驗證了該構(gòu)型的輔助梁難以滿足抗鳥撞要求的事實。圖26的后板變形明顯，表明穿透輔助梁后鳥體仍有較大能量，分析和試驗的應(yīng)變對比如圖29所示，LG01，BG03和LG02的分析與試驗的應(yīng)變情況一致性較好，但是由于鳥體撞擊輔助梁前，已經(jīng)穿透了一層蜂窩夾層結(jié)構(gòu)前緣，故兩者仍有一定偏差，特別是6 ms之后；而試驗中靠近撞擊點的BG01和BG02應(yīng)變片的數(shù)值已經(jīng)超過了量程，是由于試驗過程中應(yīng)變片已經(jīng)損壞。

圖24 構(gòu)型1的MID-A工況試驗前緣損傷

圖25 構(gòu)型1的MID-A工況試驗輔助梁損傷

圖26 構(gòu)型1的MID-A工況試后分析板變形

圖27 構(gòu)型1的MID-A工況試后分析前緣損傷

圖28 構(gòu)型1的MID-A工況試后分析輔助梁損傷

圖29 構(gòu)型1MID-A工況試驗與分析應(yīng)變對比

3.2 構(gòu)型3的MID-A工況試驗及試后分析結(jié)果

構(gòu)型3的MID-A工況試驗結(jié)果，如圖30～圖31所示，試驗后分析結(jié)果如圖32～圖33所示，分析和試驗中輔助梁的變形形式和范圍基本一致，即輔助梁僅發(fā)生大范圍變形而未發(fā)生失效，與3.1節(jié)的試驗結(jié)果對比可知，構(gòu)型2的輔助梁抗鳥撞能力明顯提升。分析和試驗的應(yīng)變對比如圖34所示，BG03和LG02的分析與試驗結(jié)果一致性較好，不過由于鳥體撞擊輔助梁前已經(jīng)穿透前緣，導(dǎo)致兩者之間仍存在一定差異，特別是6 ms之后；試驗中靠近撞擊點的LG01，BG01和BG02的數(shù)值已經(jīng)超過了量程，是由于試驗過程中應(yīng)變片已經(jīng)損壞。

圖30 構(gòu)型3的MID-A工況試驗前緣損傷

圖31 構(gòu)型3的MID-A工況試驗輔助梁損傷

圖32 構(gòu)型3的MID-A工況試后分析前緣損傷

圖33 構(gòu)型3的MID-A工況試后分析輔助梁損傷

圖34 構(gòu)型3MID-A工況試驗與分析應(yīng)變對比

3.3 構(gòu)型3的MID-AC工況試驗

構(gòu)型3的MID-AC工況試驗結(jié)果，如圖35～圖37所示，試驗后分析結(jié)果如圖38～圖40所示，分析與試驗結(jié)果吻合較好，較好地驗證了靠近肋的MID-AC工況的抗鳥撞性能，表明了該型輔助梁較好的抗鳥撞性；圖37的肋上有少量鳥血，是因為前緣和輔助梁連接的部分緊固件在試驗中發(fā)生失效，導(dǎo)致鳥血從兩個零件的縫隙中滲入，輔助梁本身并無裂紋。分析與試驗應(yīng)變對比如圖41所示，LG01和LG02的分析與試驗一致性較好，不過由于鳥體撞擊輔助梁前已經(jīng)穿透前緣，導(dǎo)致兩者之間仍存在一定差異，特別是6 ms之后；試驗中靠近撞擊點的BG01,BG02和BG03的數(shù)值已經(jīng)超過了量程，是由于試驗過程中應(yīng)變片已經(jīng)損壞。

圖35 構(gòu)型3的MID-AC工況試驗前緣損傷

圖36 構(gòu)型3的MID-AC工況試驗輔助梁損傷

圖37 構(gòu)型3的MID-AC工況結(jié)構(gòu)損傷(側(cè)后方)

圖38 構(gòu)型3的MID-AC工況試后分析前緣損傷

圖39 構(gòu)型3的MID-AC工況試后分析輔助梁損傷

圖40 構(gòu)型3的MID-AC工況試后分析(側(cè)后方)

圖41 構(gòu)型3MID-AC試驗與分析應(yīng)變對比

3.4 結(jié)果總結(jié)

對3.1節(jié)～3.3節(jié)的試驗及對應(yīng)試驗后分析結(jié)果進行梳理，可得結(jié)果匯總于表6。

表6 試驗結(jié)果及試驗后分析結(jié)果

將3.1節(jié)～3.3節(jié)的試驗結(jié)果進行匯總，可以得到圖42中的損傷范圍示意圖。

圖42 鳥撞試驗前緣蒙皮及輔助梁損傷范圍

3.5 結(jié)構(gòu)抗鳥撞機理分析

由3.3節(jié)和3.4節(jié)可知，3次試驗和試驗后分析的前緣蒙皮的損傷形式基本一致，因為兩個構(gòu)型試驗件的差異主要在輔助梁及對應(yīng)的肋，對前緣基本無影響。

但是構(gòu)型1(MID-A工況)和構(gòu)型3(MID-A和MID-AC工況)的輔助梁損傷形式完全不同：

構(gòu)型1的MID-A工況中，輔助梁發(fā)生穿透，長度約610 mm，未發(fā)生穿透區(qū)域塑性變形很小。

而構(gòu)型3的MID-A工況和MID-AC工況中，輔助梁發(fā)生了非常大的變形，但是沒有發(fā)生穿透。這是因為相比構(gòu)型1，構(gòu)型3的輔助梁為鈑金件，沒有機加輔助梁中的加強筋條及各種硬點，避免了剛度突變；同時構(gòu)型3使用的2024-T42材料雖然強度較低，但是其延展性明顯優(yōu)于用于制造機加件的7050-T7451，使得其可以發(fā)生很大變形而不失效，并因此吸收大量能量。

由圖30和圖35還可以看出，構(gòu)型3的兩次試驗中，連接前緣蒙皮和輔助梁的緊固件發(fā)生了大量失效，而圖24的構(gòu)型1中相應(yīng)緊固件失效的數(shù)量要少得多，其原因如圖43所示，構(gòu)型3的前緣發(fā)生鳥撞時，輔助梁發(fā)生很大變形，從而使得其與前緣蒙皮連接的緊固件處于剪切與拉伸的復(fù)合受載中，嚴酷的受載環(huán)境使得這些緊固件大量失效，并因此吸收更多能量，這也是此構(gòu)型抗鳥撞性能更好的重要原因；而構(gòu)型1的機加輔助梁在試驗中發(fā)生的變形小得多，對應(yīng)的緊固件的所受載荷以剪切為主，載荷水平較低，在緊固件被大量剪斷前，輔助梁已經(jīng)發(fā)生穿透，隨后緊固件載荷迅速下降，因此試驗中緊固件失效較少。

圖43 前緣蒙皮與鈑金輔助梁連接緊固件失效原理示意圖

從圖43還可以看出，橫截面帶有“波紋加強筋”的構(gòu)型3，在鳥撞過程中，“波紋”被拉伸變平，因此和構(gòu)型2的相比，此構(gòu)型將吸收更多能量。此特征也可以從圖31和圖36中觀察到。

綜上所述，本文設(shè)計的構(gòu)型3中的輔助梁，會通過輔助梁發(fā)生大變形、緊固件大量失效、“波紋”被拉伸變形等形式吸能，提升結(jié)構(gòu)的抗鳥撞性能。

4 結(jié) 論

本文根據(jù)某型民機使用的帶機加輔助梁的尾翼前緣原始結(jié)構(gòu)，提出了兩種使用鈑金輔助梁的尾翼前緣新方案，并通過數(shù)值計算選擇了抗鳥撞性能更好的方案作為優(yōu)化設(shè)計構(gòu)型。對原構(gòu)型和優(yōu)化設(shè)計構(gòu)型進行了抗鳥撞性能驗證及機理分析。根據(jù)本文研究，可以得到如下結(jié)論：

(1)數(shù)值計算與試驗結(jié)果相似度很好，結(jié)構(gòu)的失效形式和破壞程度接近，說明本文采用的分析方法合理，可以用于抗鳥撞性能評估，并指導(dǎo)結(jié)構(gòu)設(shè)計。

(2)采用機加工形式的輔助梁的前緣結(jié)構(gòu)，在遭受鳥撞時，輔助梁沒有產(chǎn)生太大塑性變形，而是發(fā)生穿透性破壞，其整體抗鳥撞性能較差。

(3)采用鈑金輔助梁的前緣結(jié)構(gòu)，在遭受鳥撞時，輔助梁會發(fā)生很大的塑性變形，緊固件大量失效，從而吸收大量鳥體能量，表現(xiàn)出較好的抗鳥撞性能；在輔助梁增加“波紋加強筋”設(shè)計后，其抗鳥撞能力進一步提升。

(4)相比機加工輔助梁，鈑金形式輔助梁質(zhì)量更輕，在實際應(yīng)用中具有明顯優(yōu)勢。

總之，本文提出的使用鈑金輔助梁的尾翼前緣結(jié)構(gòu)抗鳥撞性能較好，且具有質(zhì)量優(yōu)勢，可以應(yīng)用到民機尾翼前緣抗鳥撞設(shè)計中。