李廣超，何 江，林貴平

電脈沖除冰系統(tǒng)的結(jié)構(gòu)力學(xué)性能分析

李廣超1，*，何 江2，林貴平1

（1.北京航空航天大學(xué)航空科學(xué)與工程學(xué)院，北京100083； 2.中國空間技術(shù)研究院，北京100081）

電脈沖除冰（EIDI）作為一種電動-機(jī)械式除冰系統(tǒng)，具有廣闊的應(yīng)用前景。本文采用線圈和金屬平板組成的電脈沖系統(tǒng)作為EIDI結(jié)構(gòu)動力學(xué)模型的研究對象，基于有限元分析軟件MSC／NASTRAN，建立了一套針對典型平板-線圈結(jié)構(gòu)的EIDI系統(tǒng)結(jié)構(gòu)動力學(xué)仿真模型，研究了模型的固有特性（模態(tài)與頻率）與時域響應(yīng)分析，并開展了仿真模型的實驗驗證。研究了基于模型的EIDI系統(tǒng)結(jié)構(gòu)動力學(xué)特性，重點(diǎn)關(guān)注最大應(yīng)變和位移在空間的分布規(guī)律，并討論了電容、電壓等放電參數(shù)的影響。研究表明，電脈沖作用下，平板最大位移分布特性近似“馬鞍型”，提升電容和電壓可以增大脈沖力峰值和結(jié)構(gòu)變形量，同時減小固定位置處的應(yīng)變。

飛機(jī)防冰除冰；電脈沖除冰（EIDI）；有限元分析；模態(tài)分析；應(yīng)力-應(yīng)變

在結(jié)冰條件下，飛機(jī)機(jī)翼前緣和進(jìn)氣道口易發(fā)生結(jié)冰現(xiàn)象，這將導(dǎo)致升力下降、飛行阻力增大，進(jìn)而引起飛機(jī)的操縱性和穩(wěn)定性品質(zhì)惡化，嚴(yán)重時甚至引起飛機(jī)失事［1-2］。因此，飛機(jī)結(jié)冰防護(hù)系統(tǒng)是現(xiàn)代飛機(jī)中必不可少的功能性組件。電脈沖除冰（EIDI）作為一種電動-機(jī)械式除冰系統(tǒng)，主要原理是通過機(jī)翼內(nèi)部線圈產(chǎn)生的電脈沖激勵機(jī)翼蒙皮高頻振蕩來除冰，具有耗能低、重量小、易維護(hù)和可靠性高等優(yōu)點(diǎn)，具有廣闊的應(yīng)用前景［3-5］。

對EIDI系統(tǒng)的研究，可分為電動力學(xué)研究、結(jié)構(gòu)動力學(xué)研究和應(yīng)用試驗研究等方向。在結(jié)構(gòu)動力學(xué)方面，主要研究系統(tǒng)結(jié)構(gòu)對脈沖力的動力響應(yīng)特性，如加速度、應(yīng)變等。在脈沖力作用下，蒙皮的真實響應(yīng)是沿三維空間分布且隨時間變化的。Bernhart和Gien［6］采用Love-Timoshenko運(yùn)動關(guān)系式建立了閉合的狀態(tài)方程組。該方法考慮了橫向剪切力的合成以及縱向債性對于結(jié)構(gòu)的影響，能夠以較高精度預(yù)估蒙皮的響應(yīng)能力。然而，真實的EIDI系統(tǒng)具有極為復(fù)雜的邊界條件，其中絕大多數(shù)無法利用此方法直接獲得精確解。因此該方法僅具有理論分析的價值，無法處理實際問題。

隨著計算機(jī)技術(shù)的不斷發(fā)展，采用有限元方法仿真求解結(jié)構(gòu)動力學(xué)問題已成為當(dāng)前的主要研究手段。Zumwalt［7］率先對鋁質(zhì)矩形平板在電脈沖作用下的結(jié)構(gòu)響應(yīng)展開研究，獲得了線圈附近法向脈沖力的分布形式；此外，計算了鋁板受上述脈沖力作用下的動態(tài)響應(yīng)（應(yīng)力和加速度隨時間的變化），與實驗情況進(jìn)行對比，計算得到的正加速度峰值比實驗值小12%，應(yīng)力最大值比實驗值小5%。

Gien［8］采用有限元模型研究理想半圓柱前緣結(jié)構(gòu)在脈沖力作用下的結(jié)構(gòu)變形。為了驗證模型的正確性，首先在簡單邊界條件下對比模型結(jié)果與閉合狀態(tài)方程的精確解。結(jié)果顯示，對于正加速度的峰值二者誤差約為7%，負(fù)加速度峰值誤差約為4%，從而驗證了模型的正確性。在此基礎(chǔ)上，對半圓柱試驗結(jié)果與該有限元方法的計算結(jié)果進(jìn)行對比，發(fā)現(xiàn)計算與試驗的誤差約為7%。

Khatkhate等［3］采用MSC／NASTRAN有限元軟件，對鋁質(zhì)蒙皮和積冰層建立模型，重點(diǎn)研究了冰層與蒙皮之間的剪切應(yīng)力。同時，將計算結(jié)果與NASA Lewis研究中心開展的一系列針對EIDI系統(tǒng)的冰風(fēng)洞試驗結(jié)果進(jìn)行對比，驗證了該方法可用來預(yù)測除冰區(qū)域。Scavuzzo和Chu［9］在此基礎(chǔ)上進(jìn)行了改進(jìn)，考慮了蒙皮振動時的阻尼現(xiàn)象。他們的計算結(jié)果顯示，蒙皮與冰層之間的剪切應(yīng)力在線圈附近較大，從脈沖位置處沿弦向迅速衰減。

Labeas等［10］在前人研究基礎(chǔ)上，利用ANSYS有限元仿真軟件對EIDI系統(tǒng)的除冰過程進(jìn)行了仿真研究，分析了冰層厚度、線圈位置、線圈數(shù)、脈沖振幅和前緣曲率半徑對除冰效果的影響。

國內(nèi)李清英等［11-12］搭建了包括電路系統(tǒng)、除冰實驗機(jī)架與結(jié)除冰環(huán)境實驗艙的EIDI系統(tǒng)實驗臺，完成了脈沖放電電路的性能測試，同時也完成了振動實驗如加速度、位移、模態(tài)實驗，進(jìn)行了EIDI實驗。張文杰等［13］研究了冰層失效的范圍，運(yùn)用動力學(xué)分析了冰-鋁板界面之間的應(yīng)力狀態(tài)。

本文采用線圈和金屬平板組成的電脈沖系統(tǒng)作為EIDI結(jié)構(gòu)動力學(xué)模型的研究對象，使用商業(yè)有限元軟件MSC／NASTRAN，對實驗對象進(jìn)行固有特性（模態(tài)與頻率）與時域響應(yīng)分析，并利用實驗數(shù)據(jù)進(jìn)行對比和驗證，最后重點(diǎn)研究了基于模型的EIDI系統(tǒng)結(jié)構(gòu)動力學(xué)特性，蒙皮表面應(yīng)變與位移隨時間變化的規(guī)律及最大應(yīng)變和位移在空間的分布規(guī)律。

1 研究方法

1．1 有限元建模

本文選用一塊長方形鋁板作為飛機(jī)蒙皮的簡化計算／實驗?zāi)Ｐ停鐖D1所示，鋁板長300 mm、寬200mm、厚度2.5mm。4個邊角及長邊中點(diǎn)位置共有6個直徑為10mm的固定圓孔，實驗中通過螺栓固定。

利用同軸雙圓環(huán)電流電磁理論建立了EIDI電動力學(xué)數(shù)學(xué)模型，數(shù)值求解獲得脈沖載荷在時間域和空間域上的分布。脈沖力 F和通電線圈電流I隨放電時間t變化的示例見圖2（a）；不同時刻，脈沖力 F（r）沿線圈半徑方向 r分布見圖2（b）（r0為線圈半徑），本文所使用的EIDI電動力學(xué)數(shù)學(xué)模型詳見文獻(xiàn)［14］。

將圖2分析獲得的脈沖載荷的作用區(qū)域簡化為一個受力中心點(diǎn)與平板中心同軸的圓形。受力區(qū)域的直徑為63.5mm，在整個受力區(qū)內(nèi)，沿半徑方向等間距劃分為10個同心的圓環(huán)。假定在每一個圓環(huán)內(nèi)，所受的力是平均作用在各圓環(huán)面積上的。在有限元模擬中，以壓強(qiáng)方式進(jìn)行載荷加載。

圖1 機(jī)翼蒙皮簡化模型Fig.1 Simplified model of wing skin

圖2 脈沖載荷的時空的分布Fig.2 Impulse load distribution in space and time

本文有限元建模采用平面模型。4個邊角處的圓孔相對于整個平板的尺寸不可忽略，因此有限元建模時必須保留其存在。為了方便載荷加載，在平板的中心建立了10個同心圓環(huán)。在網(wǎng)格劃分的時候，將圓環(huán)之間共用的結(jié)點(diǎn)重合。模型材料為硬鋁，彈性模量為70 GPa，泊松比為0.33，密度為2 700 kg／m3。通過與實驗數(shù)據(jù)的對比，將結(jié)構(gòu)阻尼系數(shù)設(shè)置為0.1。模型的約束方式和載荷加載區(qū)域如圖 3所示?？紤]到實際物理模型是通過埋螺帽固定平板的4個圓孔，在有限元模型中，約束方式是固定圓孔邊上的結(jié)點(diǎn)的x、y、z 3個方向的位移Ux、Uy、Uz。

圖3 模型的約束方式和載荷加載區(qū)域Fig.3 Model constraint condition and load adding region

1．2 實驗方法與設(shè)備

為驗證有限元計算的可靠性，本文還在已有EIDI實驗臺［15］上開展了鋁板結(jié)構(gòu)響應(yīng)的實驗測量，測量內(nèi)容包括鋁板表面應(yīng)變與位移。實驗件與計算模型相同，電脈沖發(fā)生裝置的充電電壓實驗在600～800 V范圍內(nèi)可調(diào)，充電電容在200～400μF范圍內(nèi)可調(diào)，鋁板與線圈間距為2.5 mm。在圖1所標(biāo)示的1～8號測點(diǎn)位置粘貼應(yīng)變片。EIDI實驗臺與應(yīng)變測量系統(tǒng)見圖4。位移測量采用激光非接觸測試系統(tǒng)。

圖4 EIDI實驗臺與應(yīng)變測量系統(tǒng)Fig.4 EIDI experimental platform and strain measurement system

2 結(jié)果與分析

2．1 模態(tài)分析

為了了解模型的振動以哪階模態(tài)為主，首先對有限元模型進(jìn)行模態(tài)分析，得出模型的固有頻率。表1給出了模型的固有頻率。

采用模態(tài)疊加法對模型進(jìn)行瞬態(tài)響應(yīng)分析，得出中心點(diǎn)的位移如圖5（a）所示。最大位移為1.619mm，發(fā)生在時間t＝2.39×10-3s。圖5（b）給出了中心點(diǎn)位移數(shù)據(jù)的頻譜分析?？梢钥闯?，中心點(diǎn)的位移主要以149.92Hz以及599.7 Hz的振動為主。對比模型的固有頻率可知，模型受到?jīng)_擊載荷作用下，主要以第1階模態(tài)和第7階模態(tài)振動為主。

圖5 模型中心點(diǎn)的瞬態(tài)響應(yīng)Fig.5 Transient response ofmodel centroid

2．2 實驗驗證

根據(jù)仿真模型，可求解獲得蒙皮表面任意位置處在整個脈沖過程中的應(yīng)變 ε。為了驗證仿真結(jié)果的正確性，將其與實驗測量的結(jié)果進(jìn)行對比。圖6對比了圖1中5號測點(diǎn)處實驗測量與有限元仿真所獲得的應(yīng)變的時間歷程。如圖6所示，在實驗參數(shù)范圍內(nèi)，仿真模型能夠準(zhǔn)確預(yù)測出應(yīng)變峰值及達(dá)到峰值的時間。在峰值之后的過程中，仿真結(jié)果與實驗數(shù)據(jù)的變化趨勢基本相同。但仿真分析的數(shù)值普遍低于實驗，其原因可能是有限元模型中結(jié)構(gòu)振動的阻尼特性設(shè)置偏大，導(dǎo)致蒙皮在經(jīng)過第1次振動后，振動幅度快速降低。

此外，還將模型中心點(diǎn)處沿外法線方向的位移進(jìn)行了對比驗證，如圖7所示。仿真結(jié)果和實驗數(shù)據(jù)的中心點(diǎn)位移時間歷程幾乎相同，尤其是峰值（或谷值）出現(xiàn)的時刻以及幅值能夠較好的匹配。綜上，仿真模型的計算結(jié)果與實驗測量結(jié)果在應(yīng)變率與法向位移兩方面均能較好吻合，由此驗證了本文有限元仿真的可靠性。

圖6 5號測點(diǎn)處應(yīng)變的時間歷程Fig.6 Time histories of strain on No.5 station

圖7 模型中心點(diǎn)位移的時間歷程Fig.7 Time histories of displacement ofmodel centroid

2．3 最大應(yīng)變和位移的空間分布

最大位移（或應(yīng)變），是指在一個電脈沖信號周期內(nèi)蒙皮表明所達(dá)到的最大變化值zmax。該參數(shù)是一個與時間無關(guān)的物理量，僅取決于位置點(diǎn)坐標(biāo)（x，y），即zmax＝f（x，y），便于理解蒙皮表面的結(jié)構(gòu)變形特性。

圖8給出用有限元仿真得到的最大位移和最大應(yīng)變在蒙皮表面的分布（EIDI系統(tǒng)參數(shù)：電容200μF、電壓600 V、間距2.5mm）。此工況下，脈沖力的最大值為730N。如圖8所示，最大位移的分布特性近似“馬鞍型”。此外，在線圈對應(yīng)位置處的應(yīng)變較大，且沿半徑方向逐步遞減，而在長邊處應(yīng)變再次增大。這是由邊界條件決定的：4個頂點(diǎn)固定的約束條件，決定了在較長的邊界上將會產(chǎn)生較大的變形。在圖8（b）中，4個固定點(diǎn)附近的應(yīng)變遠(yuǎn)大于平板其他位置。如果不考慮這些固定點(diǎn)，則線圈對應(yīng)位置處、靠近y軸與蒙皮交匯處、靠近x軸與蒙皮交匯處的地方，應(yīng)變較大。

進(jìn)一步的仿真結(jié)果表明，增大EIDI系統(tǒng)參數(shù)至電容400μF、電壓700 V，脈沖力的峰值將增大至1800N。此工況下，最大應(yīng)變、最大位移的數(shù)值（如圖9所示）均大于圖8所示工況，最大位移（見圖9（a））在蒙皮表面的分布規(guī)律與圖8（a）無異。但線圈對應(yīng)位置處的最大應(yīng)變明顯高于其他位置的結(jié)果，甚至高于固定螺栓位置。導(dǎo)致這一結(jié)果的原因，主要是脈沖力的絕對值改變了蒙皮表面的振動特性，此外，脈沖力的作用時間也發(fā)現(xiàn)了顯著改變。在該工況下，蒙皮的結(jié)構(gòu)變形量得到顯著提升，有利于除冰效果；另一方面，固定位置處的變形量反而減小，避免了疲勞損傷。

圖8 蒙皮表面在1個電脈沖周期內(nèi)的最大響應(yīng)（EIDI系統(tǒng)參數(shù)：電容200μF、電壓600 V、間距2.5mm）Fig.8 Maximum skin response within one electro-impulse period（EIDI system parameters：capacitance equals to 200μF，voltage equals to 600 V，spacing equals to 2.5mm）

圖9 蒙皮表面在1個電脈沖周期內(nèi)的最大響應(yīng)（EIDI系統(tǒng)參數(shù)：電容400μF、電壓700 V、間距2.5mm）Fig.9 Maximum skin response within one electro-impulse period（EIDI system parameters：capacitance equals to 400μF，voltage equals to 700 V，spacing equals to 2.5mm）

3 結(jié) 論

本文基于有限元分析軟件MSC／NASTRAN，建立了一套針對典型平板-線圈結(jié)構(gòu)的EIDI系統(tǒng)結(jié)構(gòu)動力學(xué)仿真模型。模型計算的應(yīng)變和位移結(jié)果與實驗測量數(shù)據(jù)吻合較好，模型的正確性得以驗證。

根據(jù)有限元仿真結(jié)果，提出了最大應(yīng)變和最大位移的概念，分析了這2個參數(shù)在蒙皮表面的分布規(guī)律和特性：

1）本文給定的邊界條件下，最大位移的分布特性近似“馬鞍型”。此外，在線圈對應(yīng)位置處的應(yīng)較大，且沿半徑方向逐步遞減，而在長邊處應(yīng)變再次增大。

2）電容和電壓由200μF和600 V提升至400μF和700V后，脈沖力峰值增大，且振動特性發(fā)生變化。一方面使得結(jié)構(gòu)變形量顯著提升，另一方面削弱了固定位置處的應(yīng)變。

本文研究結(jié)果為后續(xù)除冰準(zhǔn)則關(guān)系的建立提供了數(shù)據(jù)支撐。

（References）

［1］鐘長生，洪冠新.飛機(jī)翼面結(jié)冰對飛行特性影響的研究［J］.航空科學(xué)技術(shù)，2004（3）：32-34.

ZHONG C S，HONG G X.Review about effects of lifting surface icing on aircraft flight characteristics［J］.Aeronautical Science＆Technology，2004（3）：32-34（in Chinese）.

［2］LYNCH F T，KHODADOUST A.Effects of ice accretions on aircraft aerodynamics［J］.Progress in Aerospace Sciences，2001，37（8）：669-767.

［3］KHATKHATE A A，SCAVUZZO R J，CHU M L.A finite element study of the EIDI system［C］∥Proceedings of the 26 th AIAA Aerospace Sciences Meeting.Reston：AIAA，1988：1-8.

［4］THOMAS S K，CASSONIR P，MACARTHUR C D.Aircraft anti-icing and de-icing techniques and modeling［J］.Journal of Aircraft，1996，33（5）：841-854.

［5］周莉，徐浩軍，龔勝科.飛機(jī)結(jié)冰特性及防除冰技術(shù)研究［J］.中國安全科學(xué)學(xué)報，2010，20（6）：105-110.

ZHOU L，XU H J，GONG S K.Research of aircraft icing characteristics and anti-icing and de-icing technology［J］.China Safety Science Journal，2010，20（6）：105-110（in Chinese）.

［6］BERNHARTW D，GIEN P H.A structural dynamics investigation related to EIDI applications：AIAA-1986-0550［R］.Reston：AIAA，1986.

［7］ZUMWALT GW.Analysis and tests for design of an electro-impulse de-icing system：NASA／CR-1987-174919［R］.Washington，D.C.：NASA，1987.

［8］GIEN PH.Free vibration and transient analysis of a semi-cylinder subjected to electro-impulse de-icing forces［M］.W ichita，Kansas：W ichita State University，1985：19-37.

［9］SCAVUZZO R J，CHU M L.Finite element studies of the electro impulse de-icing system［J］.Journal of Aircraft，1990，27（9）：757-763.

［10］LABEASG N，DIAMANTAKOS I D，SUNARIC M M.Simulation of the electroimpulse de-icing process of aircraftwings［J］. Journal of Aircraft，2006，43（6）：1876-1885.

［11］李清英，朱春玲，白天.電脈沖除冰系統(tǒng)的除冰實驗與數(shù)值模擬［J］.航空動力學(xué)報，2012，27（2）：350-356.

LIQ Y，ZHU C L，BAI T.De-icing experiment and numerical simulation of the electro-impulse de-icing system［J］.Journal of Aerospace Power，2012，27（2）：350-356（in Chinese）.

［12］李清英，朱春玲，白天.電脈沖除冰系統(tǒng)激勵簡化與影響因素［J］.航空學(xué)報，2012，25（8）：1384-1393.

LIQ Y，ZHU C L，BAI T.Simplification of de-icing excitation and influential factors of the electro-impulse de-icing system［J］.Acta Aeronautica et Astronautica Sinica，2012，25（8）：1384-1393（in Chinese）.

［13］張文杰，董文俊，王斌團(tuán).電脈沖除冰仿真冰層松脫準(zhǔn)則研究［J］.計算機(jī)工程與應(yīng)用，2012，48（3）：232-233.

ZHANG W J，DONGW J，WANG B T.Study on de-icing criterion of electro-impulse de-icing simulation［J］.Computer Engineering and Applications，2012，48（3）：232-233（in Chinese）.

［14］李廣超，何江，林貴平.電脈沖除冰（EIDI）系統(tǒng)電動力學(xué)模型分析［J］.航空動力學(xué)報，2011，26（1）：54-59.

LIG C，HE J，LIN G P.Research of electrodynamics in electroimpulse de-icing（EIDI）system［J］.Journal of Aerospace Power，2011，26（1）：54-59（in Chinese）.

［15］李廣超，何江，林貴平.電脈沖除冰（EIDI）技術(shù)研究［J］.航空動力學(xué)報，2011，26（8）：1718-1735.

LIG C，HE J，LIN G P.Electro-impulse de-icing（EIDI）technology study［J］.Journal of Aerospace Power，2011，26（8）：1718-1735（in Chinese）.

Tel.：010-82339235

E-mail：liguangchao＠buaa.edu.cn

Perform ance analysis on structural dynam ics of EIDIsystem

LIGuangchao1，*，HE Jiang2，LIN Guiping1

（1.School of Aeronautic Science and Engineering，Beijing University of Aeronautics and Astronautics，Beijing 100083，China；2.China Academy of Space Technology，Beijing 100081，China）

As an electromechanical de-icing system，electro-impulse de-icing（EIDI）technology has broad application prospect.In the present paper，structural dynamics characteristics of a simplified EIDI system consisting of coil and metal p late were numerically and experimentally studied.A finite element simulation model of an EIDI system was constructed，and the structural dynamics was analyzed by FEM software MSC／NASTRAN.The inherent characteristics（modal and frequency）as well as the transient response of themodel were analyzed.The simulation was then validated by the experimentalmeasurement.The structure dynamics characteristics of an EIDImodel were studied.Themaximum strain and disp lacement of the skin surface characteristic locations were evaluated with different capacitance，voltage and clearance between the skin and the coils.It is found that under the electro-impulse，themaximum disp lacement of the metal plate presents a saddle-like distribution.Increasing the capacitance and voltage will increase the impulse peak value and the structural deformation，and meanwhile，the strain of fixed location will be decreased.

aircraft anti-icing and de-icing；electro-impulse de-icing（EIDI）；finite element analysis； modal analysis；stress-strain

2015-09-17；Accep ted：2015-12-21；Pub lished online：2016-01-20 16：40

Aeronautical Science Foundation of China（2013ZC51030）

V244.1+5

A

1001-5965（2016）10-2069-06

李廣超 男，博士，副研究員。主要研究方向：流動控制、飛機(jī)除防冰。

http：∥bhxb.buaa.edu.cn jbuaa＠buaa.edu.cn

DO I：10.13700／j.bh.1001-5965.2015.0609

2015-09-17；錄用日期：2015-12-21；網(wǎng)絡(luò)出版時間：2016-01-20 16：40

www.cnki.net／kcms／detail／11.2625.V.20160120.1640.004.htm l

航空科學(xué)基金（2013ZC51030）

*通訊作者：Tel.：010-82339235 E-mail：liguangchao＠buaa.edu.cn

李廣超，何江，林貴平.電脈沖除冰系統(tǒng)的結(jié)構(gòu)力學(xué)性能分析［J］.北京航空航天大學(xué)學(xué)報，2016，42（10）：2069-2074. LIG C，HE J，LIN G P.Performance analysis on structural dynamics ofEIDIsystem［J］.Journal ofBeijing University ofAeronautics and Astronautics，2016，42（10）：2069-2074（in Chinese）.

URL：www.cnki.net／kcms／detail／11.2625.V.20160120.1640.004.htm l

*Correspond ing au thor.Tel.：010-82339235 E-mail：liguangchao＠buaa.edu.cn