郭軍 符偉超 李旭陽 盧麗穎

摘要：通過對某民用飛機天線罩結(jié)構(gòu)特點的分析及鳥撞計算分析所需條件，截取了天線罩及機體結(jié)構(gòu)的幾何模型；根據(jù)該天線罩幾何模型的各零部件特點，采用了不同的單元模式劃分有限元網(wǎng)格；針對連接接頭的受力特性，利用軸連接副和球連接副實現(xiàn)了對上下接頭連接處運動副的模擬；通過光滑粒子流體動力學（SPH）鳥撞擊天線罩有限元模型，模擬了天線罩鳥撞過程；最后，通過對不同鳥撞工況下天線罩結(jié)構(gòu)的受力情況分析，為民機天線罩的鳥撞試驗驗證方案制訂提供參考。

關(guān)鍵詞：天線罩；鳥撞；運動副；有限元；SPH

中圖分類號：V226文獻標識碼：ADOI：10.19452/j.issn1007-5453.2021.12.014

鳥撞事故是飛行器常見的事故之一。對于載人飛行器來說，鳥撞事故有可能造成機毀人亡的嚴重危害[1]。根據(jù)美國聯(lián)邦航空局（FAA）統(tǒng)計，1988—2021年，飛機因與野生動物相撞而發(fā)生損毀的數(shù)量高達263架，2017年野生動物與飛機相撞事故相較2016年增加了1069起，其中鳥撞比例高達95%[2]。民用航空對飛機的結(jié)構(gòu)抗鳥撞設(shè)計有明確的要求，對飛機的機頭、風擋、機翼及尾翼均有具體的條款規(guī)定。對于暴露在機體結(jié)構(gòu)外面，且有可能發(fā)生鳥撞并引發(fā)飛行安全的部件，民航適航管理部門對其都有抗鳥撞要求。鳥類撞擊試驗提供了一種檢測鳥類撞擊力的直接方法。然而，飛機結(jié)構(gòu)的設(shè)計通常涉及設(shè)計—制造—測試的多次迭代，并且進行鳥類撞擊試驗不僅耗時而且成本高。此外，試驗數(shù)據(jù)往往高度分散，給結(jié)構(gòu)設(shè)計帶來了障礙。而鳥撞試驗的仿真模擬減少了試驗耗時長、成本高的缺點，近年來幾乎成為學術(shù)研究與工程應(yīng)用中必備的方法[3-4]。

國內(nèi)外在鳥撞模擬中常使用圓柱體、半球端圓柱體、橢圓體和球體來替代真實的鳥體，該方法可以反映鳥軀干的主要質(zhì)量和形狀對碰撞過程的影響[5]。Frederik等[6]通過試驗與仿真相結(jié)合的方法對鳥類在撞擊過程中鳥類材料的沖擊與穩(wěn)態(tài)壓力進行了分析，揭示了傾斜角等因素對沖擊狀態(tài)的影響。賈建東等[7]采用SPH方法對圓弧風擋受鳥體撞擊的過程進行了數(shù)值模擬，得到了結(jié)構(gòu)的受載與損傷情況，為風擋安全的設(shè)計提供了參考。楊瑞進等[8]基于光滑粒子流體動力學（SPH）方法對航空發(fā)動機進氣風扇抗鳥撞性能進行了研究，仿真結(jié)果與試驗結(jié)果吻合良好。Grimaldi等[9]采用SPH方法研究了風擋幾何形狀、撞擊角度和擋板曲率對結(jié)構(gòu)撞擊響應(yīng)的影響。閻軍等[10]針對寬弦風扇葉片空腔結(jié)構(gòu)采用加權(quán)求和的方式實現(xiàn)了不同位置鳥撞工況下風扇結(jié)構(gòu)性能的統(tǒng)一評價和設(shè)計。飛機用于導(dǎo)航及通信的天線，基本安裝在機體結(jié)構(gòu)外部，通過加裝天線罩來保持飛機的氣動外形，同時避免天線受自然環(huán)境的鹽類物質(zhì)侵蝕破壞[11]。天線罩在受到鳥撞后，天線及機體結(jié)構(gòu)均有可能出現(xiàn)嚴重損傷，危及飛機的安全飛行。Sebastian等[12]應(yīng)用積木式驗證方法，從材料級到天線罩部件級對天線罩結(jié)構(gòu)性能進行了表征，為天線罩抗鳥撞的仿真模擬提供了參考。白杰等對CFM56-5B型發(fā)動機風扇葉片抗鳥撞性能進行了有限元仿真，得出鳥撞葉片的機械損傷程度與鳥的速度、受撞擊位置之間的關(guān)系。綜上所述，目前對飛機天線罩在受到鳥撞載荷下的響應(yīng)相關(guān)分析較少，有必要對其進行鳥撞分析，以評估飛鳥對飛機安全飛行的影響。

1幾何模型選取

本文的研究目的是分析鳥撞后天線罩的響應(yīng)以及飛機機體結(jié)構(gòu)承受載荷的情況。所以幾何模型必須包含天線罩、機體結(jié)構(gòu)、天線罩與機體結(jié)構(gòu)的連接結(jié)構(gòu)及天線假件。

對于機體結(jié)構(gòu)的選取，因在鳥撞仿真計算時，邊界條件采用四邊固支，為消除邊界固支的剛度影響，集合模型最終選取了9個隔框（10跨）、10個長桁（11跨）范圍內(nèi)的機體結(jié)構(gòu)。其他主要結(jié)構(gòu)包括天線罩、天線假件、天線安裝底座、連接機體與底座的7組接頭以及加強耳片等，如圖1、圖2所示。

2有限元模型

用于鳥撞仿真計算的有限元模型，選取了主要結(jié)構(gòu)進行了網(wǎng)格劃分。根據(jù)結(jié)構(gòu)的實際尺寸，天線安裝底盤采用體單元，其他結(jié)構(gòu)均采用殼單元，各結(jié)構(gòu)之間的連接采用一維剛性連接單元，鳥體采用SPH鳥，總體有限元模型如圖3所示。

2.1材料參數(shù)

各部件材料參數(shù)如下：蒙皮、長桁及隔采用鋁2024-T351。托盤和接頭采用鋁7075-T651。天線罩為復(fù)合材料。具體參數(shù)見表1。其中，ρ為密度，λ為泊松比，E為彈性模量，Бt為切向模量，бy為屈服應(yīng)變。

2.2單元網(wǎng)格劃分

根據(jù)不同的部件及實際功能，對其進行了單元網(wǎng)格劃分。天線安裝底座和7組接頭劃分為體單元，其中底座采用六面體單元，接頭采用為四面體單元。其他結(jié)構(gòu)抽取 CATIA模型中的部件中面劃分為殼單元，主要以四邊形單元為主，極少部分為三角形單元。

2.3接頭連接形式

該天線罩及天線裝在底座上，底座與機體結(jié)構(gòu)通過7組接頭連接。考慮到機體結(jié)構(gòu)的受力問題，設(shè)計接頭連接時，不同位置接頭的承載方式不一樣。要求1、5接頭承受YZ方向的力，3、4接頭承受XZ方向的力，2、6、7承受Z方向的力。接頭位置與形式如圖4所示，其中7號接頭靠近機尾方向。

本文利用軸連接副和球連接副（見圖5），釋放了連接單元的特定自由度，實現(xiàn)了接頭設(shè)計所要求的承載方式。有限元分析軟件LS_DYNA中的軸連接副用于模擬力接頭設(shè)計的軸運動副。下部接頭的連接設(shè)計使用了球運動副和柱運動副，在DYNA中沒有現(xiàn)成的組合連接副，本文采用釋放連接副與單元連接點相應(yīng)自由度的方法實現(xiàn)組合運動副的模擬。計算結(jié)果顯示，該方法成功模擬了接頭的設(shè)計運動副狀態(tài)。

2.4鳥體選擇

鳥撞過程是兩個強度、剛度相差很大的物體發(fā)生碰撞的過程。強度小的鳥體通常會發(fā)生大變形。因此，在進行仿真模擬時，如果采用Lagrange網(wǎng)格會發(fā)生扭曲問題，同時網(wǎng)格會出現(xiàn)負體積，導(dǎo)致計算終止。SPH算法可以克服Euler方法難跟蹤變形和不識別材料界面位形的問題，同時解決了Lagrange方法在大變形下狀態(tài)下的網(wǎng)格扭曲問題。本文鳥撞分析采用SPH模型鳥體。鳥體為兩端半球體，直徑112.3mm，中間圓柱體長度112.3mm，質(zhì)量1.8kg，單元為質(zhì)量元。SPH鳥體參數(shù)見表2。

3計算分析

天線鳥撞仿真計算的初始條件為：鳥撞速度188.3m/s，方向沿航向。計算工況4個：（1）鳥體中心距離蒙皮高度125mm；（2）距離蒙皮高度225mm；（3）鳥體中心距離蒙皮高度125mm，角度30°；（4）鳥體中心距離蒙皮高度225mm，角度30°。計算采用LS-DYNA求解器，利用10核并行計算。

3.1第1工況

在工況1中，7個接頭的連接點力—時間歷程曲線如圖6所示。為直觀展示各個接頭所承受載荷情況，對接頭三個方向的力進行了統(tǒng)計，見表3。鳥撞后結(jié)構(gòu)整體的變形情況如圖7所示（剖視圖）。從圖6和表3可知，鳥撞過程4號接頭承受的瞬時載荷最大。

3.2第2工況

對照3.1節(jié)，計算結(jié)果給出了工況2下7個接頭的連接點力。其曲線如圖8所示。表4列出了工況2下各個接頭三個方向的載荷瞬時最大值。從圖8和表4可知，該工況鳥撞過程也是4號接頭承受的瞬時載荷最大。

3.3第3工況

工況3下7個接頭的連接點力時間歷程曲線如圖9所示。從圖9和表5可知，該工況鳥撞過程也是4號接頭承受的瞬時載荷最大。

3.4第4工況

工況4下7個接頭的連接點力時間歷程曲線如圖10所示。從圖10和表6可知，該工況鳥撞過程也是4號接頭承受的瞬時載荷最大。

3.5計算結(jié)果分析

在工況1和工況2下，第4個接頭的X方向瞬時力最大，鳥體距離蒙皮125mm的工況下力的最大值為60080N，距離225mm的工況下力的最大值為56450N。將工況1和工況2結(jié)果與工況3和工況4對比可以看出，改變角度后依然是第四接頭承受的瞬時載荷最大，且改變角度后，第四接頭的峰值載荷基本不變，說明鳥體撞擊角度對接頭承力的影響較小。這表明，鳥體距離蒙皮越遠接頭4所承受X方向的力越小。也就是說，該天線罩在鳥撞過程中，鳥體距離蒙皮越近，傳遞給機體結(jié)構(gòu)的力越大。

經(jīng)過對該天線罩鳥撞仿真結(jié)果分析可知，鳥撞位置距離蒙皮越近，其對機體結(jié)構(gòu)造成的危害越大，因此該天線罩的抗鳥撞驗證區(qū)域應(yīng)以選擇距離蒙皮較近的位置。

4結(jié)論

通過本文的研究，得到了某天線罩遭受不同工況下鳥撞后，通過天線底座傳遞給機體結(jié)構(gòu)的力大小規(guī)律：鳥體距離蒙皮越近，天線罩鳥撞后傳給機體結(jié)構(gòu)的力越大；相同速度下，鳥體撞擊角度對接頭承力峰值影響較小。該仿真分析結(jié)果為天線罩的鳥撞試驗設(shè)計提供了工況選擇指明了趨勢。后續(xù)研究將根據(jù)試驗結(jié)果修正鳥撞有限元模型，找出影響計算結(jié)果的因素，為其他鳥撞仿真計算的建模提供參考。

參考文獻

[1]劉小川，王計真，白春玉.人工鳥研究進展及在飛機結(jié)構(gòu)抗鳥撞中的應(yīng)用[J].振動與沖擊，2021，40（12）：80-89. Liu Xiaochuan， Wang Jizhen， Bai Chunyu. Overview on aircraft bird and application on the structural bird-strike[J]. Journal of Vibration and Shock，2021，40（12）：80-89. （in Chinese）

[2]Dobeeer R A，Begier M J. Wildlife strikes to civil aircraft in the united states，1990-2017[R]. The FAA National Wildlife Strike Database，2019.

[3]Liu J，Li Y，Gao X. Bird strike on a flat plate：experiments and numerical simulations[R]. International Journal of Impact Engineering，2014，70：21-37.

[4]高俊，耿立超，吳志斌.民用飛機方向舵抗鳥撞分析研究[J].民用飛機設(shè)計與研究，2021（1）：117-122. Gao Jun， Geng Lichao， Wu Zhibin. Bird strike resistance analysis and research on civil aircraft rudder[J]. Civil Aircraft Design & Research， 2021（1）：117-122. （in Chinese）

[5]SebastianH.Computationalmethodsforbirdstrike simulations：A review[J]. Computers and Structures，2011，89（23）：2093-2112.

[6]Frederik A，Geert L，Wim V P，et al. Numerical and experimental investigation of the shock and steady state pressuresinthebirdmaterialduringbirdstrike[J]. International Journal of Impact Engineering，2017，107：12-22.

[7]賈建東，李志強，楊建林.用SPH和有限元方法研究鳥撞飛機風擋問題[J].航空學報，2010，31（1）：136-142. Jia Jiandong， Li Zhiqiang， Yang Jianlin. A study of bird impact on aircraft windshield using SPH and finite element method[J]. Acta Aeronautica et Astronautica Sinica， 2010， 31（1）：136-142.（in Chinese）

[8]楊瑞進，姜楠，池劍虹.基于SPH方法鳥撞航空發(fā)動機進氣風扇的數(shù)值分析[J].航空科學技術(shù)，2012（6）：35-37. Yang Ruijin， Jiang Nan， Chi Jianhong. Numerical analysis based on the SPH for bird strike aeroengine intake fan[J]. Aeronautical Science & Technology， 2012（6）：35-37.（in Chinese）

[9]Grimaldi A，Sollo A，Guida M，et al. Parametric study of a SPH high velocity impact analysis-A bird strike windshield application-science direct[J]. Composite Structures，2013，96（4）：616-630.

[10]閻軍，張晨光，霍思旭.寬弦風扇葉片空腔結(jié)構(gòu)多目標輕量化設(shè)計[J].航空科學技術(shù)，2021，32（4）：60-64. Yan Jun， Zhang Chenguang， Huo Sixu. Multi-objective lightweight design of the cavity structure of wide-chord fan blades[J]. Aeronautical Science & Technology， 2021，32（4）：60-64. （in Chinese）

[11]Rafieipour H，Setoodeh A R，Lau K T. Mechanical and electromagnetic behavior of fabricated hybrid composite sandwich radome with a new optimized frequency-selective surface[J]. Composite Structures，2021，273：54-65.

[12]Sebastian H，Tim W H，Clemens B，et al. Bird strike on aircraft radome：Dynamic characterization of quartz fiber composite sandwich for accurate，predictive impact simulations[C]// EPJ Web of Conferences，2018.

Simulation Analysis of Bird Strike Calculation for Civil Aircraft Radome

Guo Jun1，F(xiàn)u Weichao1，Li Xuyang1，Lu Liying2

1. Aviation Key Laboratory of Science and Technology on Structures Impact Dynamics，Aircraft Strength Research Institute of China，Xian 710065，China

2. Aircraft Maintenance and Engineering Corporation，Beijing 100621，China

Abstract： Through the analysis of the structural characteristics of a civil aircraft radome and the conditions required for bird strike calculation and analysis， the geometric model of the radome and the body structure is intercepted. According to the characteristics of the components of the radome geometric model， different unit modes are used to divide the finite element mesh. According to the force characteristics of the connection joints， the shaft connection pair and the ball connection pair are used to simulate the motion pair of the upper and lower joints. The SPH bird impact radome finite element model is used to simulate the radome bird strike process. Finally， through the analysis of the force situation of the radome structure under different bird strike conditions， it provides reference for the development of the bird strike experimental verification scheme of the civil aircraft radome.

Key Words： radome； bird strike； motion pair； finite element； SPH