胡 偶，陳平劍

(中國直升機(jī)設(shè)計(jì)研究所，江西 景德鎮(zhèn) 333001)

旋翼槳葉載荷與槳葉結(jié)構(gòu)參數(shù)相關(guān)性研究

胡 偶，陳平劍

(中國直升機(jī)設(shè)計(jì)研究所，江西 景德鎮(zhèn) 333001)

基于CAMRAD II軟件，建立自由尾跡耦合柔性多體動力學(xué)的彈性槳葉結(jié)構(gòu)載荷分析模型，針對某直升機(jī)旋翼槳葉結(jié)構(gòu)設(shè)計(jì)與載荷分析問題，開展槳葉結(jié)構(gòu)參數(shù)與槳葉載荷水平的相關(guān)性研究。通過典型飛行狀態(tài)，包括超扭狀態(tài)和大前進(jìn)比前飛狀態(tài)，研究槳葉剖面揮舞剛度、擺振剛度、扭轉(zhuǎn)剛度等參數(shù)的變化對槳葉揮舞與擺振方向載荷水平的影響，并分析由結(jié)構(gòu)共振引起的槳葉載荷突增現(xiàn)象。

載荷；結(jié)構(gòu)參數(shù)；共振；旋翼槳葉

0 引言

直升機(jī)旋翼槳葉載荷分析是旋翼氣動設(shè)計(jì)與結(jié)構(gòu)設(shè)計(jì)的重要支撐，是旋翼系統(tǒng)方案設(shè)計(jì)與工程詳細(xì)設(shè)計(jì)階段旋翼槳葉強(qiáng)度與壽命評估的重要依據(jù)。旋翼槳葉載荷分析是典型的氣動/結(jié)構(gòu)耦合問題，且是一種綜合性的系統(tǒng)工作：首先需建立直升機(jī)模型、旋翼和尾槳模型以及槳轂-操縱系統(tǒng)模型(包括各種鉸結(jié)構(gòu)、變距拉桿等)，再將模型通過各種連接形式連接，然后基于氣動模型(如升力線方法)和多體動力學(xué)模型，開展旋翼槳葉載荷分析。載荷分析工作在旋翼系統(tǒng)設(shè)計(jì)過程中起著承上啟下的重要作用，是結(jié)構(gòu)設(shè)計(jì)、載荷分析、壽命評估這一設(shè)計(jì)流程必不可少的關(guān)鍵環(huán)節(jié)。

旋翼槳葉振動載荷評估是直升機(jī)空氣動力學(xué)與動力學(xué)領(lǐng)域面臨的難點(diǎn)。受限于復(fù)雜的氣流環(huán)境、非線性結(jié)構(gòu)彈性變形等耦合因素，采用數(shù)值分析方法很難獲得精確的振動載荷值。但從工程設(shè)計(jì)的角度來看，通過數(shù)值分析與解析分析方法，構(gòu)建高效、準(zhǔn)確的載荷工程計(jì)算方法，能夠滿足不同飛行狀態(tài)的載荷計(jì)算分析工作要求，形成可用于旋翼系統(tǒng)壽命與強(qiáng)度評估的載荷譜，從而有效縮短直升機(jī)旋翼系統(tǒng)研制周期，提高研制效率。

在國外，隨著計(jì)算技術(shù)的高速發(fā)展，形成了以CAMRAD II[1]為代表的直升機(jī)綜合氣彈分析軟件，在直升機(jī)工程設(shè)計(jì)領(lǐng)域得到了廣泛的應(yīng)用。該軟件是由Wayne Johnson[1-2]開發(fā)的一款直升機(jī)氣動與動力學(xué)綜合分析軟件，具備多體動力學(xué)、非線性有限元以及直升機(jī)空氣動力學(xué)耦合分析能力，可用于預(yù)先研究、方案設(shè)計(jì)、詳細(xì)設(shè)計(jì)等階段直升機(jī)的性能、載荷、振動、響應(yīng)及穩(wěn)定性分析。CAMRAD II軟件在旋翼槳葉載荷計(jì)算方面得到了廣泛的應(yīng)用驗(yàn)證[3-7]，其精度能夠滿足工程設(shè)計(jì)的要求。

在國內(nèi)，習(xí)娟等人[8]針對某無人直升機(jī)開展了直升機(jī)旋翼系統(tǒng)載荷分析技術(shù)研究，王浩文等[9]建立了旋翼系統(tǒng)綜合氣彈分析方法，孫韜等[10]針對典型前飛狀態(tài)開展旋翼振動載荷分析研究，吳杰等[11]比較了力積分法、反力法和曲率法三類旋翼槳葉結(jié)構(gòu)載荷計(jì)算方法的預(yù)測精度和適用范圍。

本文基于CAMRAD II軟件，建立耦合自由尾跡模型和柔性多體力學(xué)思想的彈性槳葉結(jié)構(gòu)載荷分析模型。選取典型飛行狀態(tài)，開展旋翼槳葉結(jié)構(gòu)參數(shù)與槳葉載荷水平相關(guān)性分析，通過槳葉剖面揮舞剛度、擺振剛度、扭轉(zhuǎn)剛度等參數(shù)的變化，研究槳葉揮舞與擺振方向載荷隨槳葉結(jié)構(gòu)參數(shù)的變化關(guān)系。最后，從槳葉模態(tài)分析和載荷諧波分析的角度，揭示了結(jié)構(gòu)共振引起的槳葉載荷突增現(xiàn)象。

1 旋翼槳葉載荷計(jì)算模型

直升機(jī)旋翼槳葉載荷包括靜載和動載，計(jì)算狀態(tài)包括超轉(zhuǎn)、超扭、水平飛行、水平轉(zhuǎn)彎、螺旋轉(zhuǎn)彎、自轉(zhuǎn)下滑等狀態(tài)[8]。本文分別選取具有代表性的超扭狀態(tài)和最大速度平飛狀態(tài)開展槳葉靜載荷與動載荷分析。超扭狀態(tài)模擬的是一種極限懸停狀態(tài)，旋翼總距設(shè)為旋翼失速前的最大允許值，周期變距為零，旋翼轉(zhuǎn)速為正常轉(zhuǎn)速，該狀態(tài)考核槳葉的靜載荷的極限值。

采用CAMRAD II軟件，建立彈性槳葉載荷分析模型，主要計(jì)算流程如圖1所示。

動載荷計(jì)算過程中采用的是全機(jī)配平模型，槳葉為彈性槳葉模型，取6階模態(tài)和6階諧波，自由尾跡選取Johnson模型[4]；靜載荷采用孤立旋翼模型計(jì)算，取6階模態(tài)和0階諧波，自由尾跡選取Johnson模型[4]。

2 旋翼槳葉載荷計(jì)算結(jié)果及分析

2.1 旋翼槳葉結(jié)構(gòu)參數(shù)

以某樣例直升機(jī)作為分析對象，其旋翼槳葉片數(shù)為5片，基準(zhǔn)轉(zhuǎn)速(Ω)為258 RPM，最大平飛速度(VMAX)為260km/h。為研究槳葉結(jié)構(gòu)參數(shù)對槳葉載荷的影響，首先，設(shè)計(jì)了一副基礎(chǔ)模型槳葉(方案 A)，圖2給出了計(jì)算得到的基礎(chǔ)槳葉(方案A)的固有特性。然后，在方案A槳葉的基礎(chǔ)上，通過調(diào)整槳葉翼型起始段之后剖面的揮舞剛度(EIFALP)、擺振剛度(EILAG)和扭轉(zhuǎn)剛度(GJ)等結(jié)構(gòu)參數(shù)分布，形成如表1所示的5個備選槳葉結(jié)構(gòu)方案，其中：在方案A結(jié)構(gòu)參數(shù)的基礎(chǔ)上，(-)表示結(jié)構(gòu)參數(shù)未發(fā)生明顯變化，+表示參數(shù)值提高，-表示參數(shù)值降低。

圖3展示了各方案槳葉段剖面特性參數(shù)的變化，包括槳葉剖面揮舞剛度分布與擺振剛度分布。

表1 槳葉剖面結(jié)構(gòu)參數(shù)的變化

2.2 槳葉載荷與槳葉結(jié)構(gòu)相關(guān)性分析

圖4展示了超扭狀態(tài)槳葉靜載荷與槳葉結(jié)構(gòu)參數(shù)之間的關(guān)系，表2給出了槳葉0.5R剖面處的靜載荷對比。數(shù)值模擬結(jié)果表明，槳葉段的剖面揮舞彎矩主要受槳葉剖面揮舞剛度的影響，隨著揮舞剛度的增加，槳葉載荷呈線性增加趨勢，且剖面擺振剛度變化對揮舞彎矩沒有影響；在擺振方向上，槳葉剖面擺振彎矩受擺振剛度和揮舞剛度耦合影響，比較方案C與方案D可以發(fā)現(xiàn)，雖然方案D剖面擺振剛度比方案C高7.5%，但由于其揮舞剛度低20%，反而導(dǎo)致方案D的0.5R剖面擺振要小于方案C槳葉。綜上，在將葉結(jié)構(gòu)設(shè)計(jì)時，必須重點(diǎn)考慮槳葉揮舞方向剛度的設(shè)計(jì)，因?yàn)閾]舞剛度的變化能夠同時引起揮舞與擺振方向載荷的變化。

圖5展示了最大速度平飛狀態(tài)槳葉動載荷與槳葉結(jié)構(gòu)參數(shù)之間的關(guān)系，表2給出了槳葉0.5R剖面處動載荷對比。數(shù)值模擬結(jié)果表明，槳葉段動載荷受結(jié)構(gòu)剛度的影響較為復(fù)雜，排除方案D槳葉結(jié)構(gòu)，可以發(fā)現(xiàn)隨著剖面揮舞剛度的增加，剖面揮舞彎矩增大；隨著剖面擺振剛度增加，擺振彎矩增大。但是方案D槳葉結(jié)構(gòu)，其載荷變化不滿足這一規(guī)律，主要表現(xiàn)為擺振彎矩增加顯著。

2.3 結(jié)構(gòu)共振引起的載荷突增現(xiàn)象分析

為分析方案D槳葉擺振彎矩突增的原因，對比方案D與方案E槳葉結(jié)構(gòu)，兩者的差別在于扭轉(zhuǎn)剛度不一樣，方案D的扭轉(zhuǎn)剛度相比于方案E有15%的降低，因此初步確認(rèn)是由于扭轉(zhuǎn)剛度變化導(dǎo)致載荷發(fā)生變化的。

表2 不同槳葉0.5R剖面處槳葉載荷比較

首先，通過模態(tài)分析法，取最大速度平飛狀態(tài)對應(yīng)的總距值，進(jìn)行方案D與方案E槳葉結(jié)構(gòu)的固有特性分析，得到如表3所示槳葉固有頻率特性。計(jì)算結(jié)果表明，方案D與方案E槳葉結(jié)構(gòu)前3階固有頻率值一致，但該總距條件下，方案D槳葉結(jié)構(gòu)的扭轉(zhuǎn)1階頻率相比于方案E槳葉更靠近6Ω。因此，從槳葉固有特性分析，導(dǎo)致載荷突增的原因可能是方案D槳葉6階模態(tài)與6Ω靠得近。

表3 方案D與方案E槳葉固有頻率比較

取槳葉0.5R處剖面動載荷隨方位角變化的關(guān)系，如圖6所示。揮舞彎矩方案D與方案E的載荷變化趨勢一致，而擺振彎矩則兩者有顯著的差別，方案D槳葉擺振方向載荷在周期內(nèi)幅值變化劇烈。

為此，對圖6中方案D和方案E槳葉的擺振彎矩進(jìn)行諧波分析[12]，獲得不同階次(Ω)下的擺振載荷幅值，如圖7所示。通過諧波分析可以看出，對于方案E槳葉結(jié)構(gòu)，隨著階次的增加，擺振彎矩呈減小的趨勢，6階項(xiàng)稍有增大；而方案D槳葉，其擺振6階項(xiàng)載荷突然增加，載荷呈發(fā)散趨勢。

因而從槳葉固有頻率和載荷諧波分析兩個方面可以得出，由于存在擺振-扭轉(zhuǎn)的耦合作用，方案D槳葉結(jié)構(gòu)由于扭轉(zhuǎn)剛度配置不合理，導(dǎo)致直升機(jī)以最大速度平飛時，出現(xiàn)結(jié)構(gòu)共振響應(yīng)，從而引起槳葉載荷高階項(xiàng)增加，增大了振動載荷水平。

3 結(jié)論

以某直升機(jī)旋翼為研究對象，開展了槳葉結(jié)構(gòu)設(shè)計(jì)參數(shù)與載荷的相關(guān)性分析研究?；贑AMRAD II分析軟件，建立了彈性槳葉結(jié)構(gòu)載荷計(jì)算方法，并針對槳葉結(jié)構(gòu)設(shè)計(jì)中的關(guān)鍵參數(shù)：剖面揮舞剛度、擺振剛度和扭轉(zhuǎn)剛度的分布等，研究結(jié)構(gòu)參數(shù)對槳葉揮舞彎矩與擺振彎矩的影響。研究結(jié)果表明，槳葉揮舞與擺振方向載荷主要受揮舞剛度與擺振剛度影響；由于擺振-扭轉(zhuǎn)耦合作用，當(dāng)扭轉(zhuǎn)剛度布置不合理時，可能引起載荷共振響應(yīng)現(xiàn)象的產(chǎn)生。通過直升機(jī)旋翼槳葉結(jié)構(gòu)與槳葉載荷的相關(guān)性分析，摸清槳葉結(jié)構(gòu)參數(shù)對槳葉載荷水平的影響，可為槳葉結(jié)構(gòu)設(shè)計(jì)提供技術(shù)支撐，并避免可能產(chǎn)生的結(jié)構(gòu)共振響應(yīng)現(xiàn)象。

[1] Johnson W. CAMRAD II, Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics[J]. Johnson Aeronautics, 1992.

[2] Johnson W. CAMRAD II, Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics, Rotorcraft Applications[J]. Johnson Aeronautics, 1993.

[3] Johnson, W. Rotorcraft aerodynamics models for a comprehensive analysis[C].Presented at the 54th Annual Forum of the AHS, Washington, DC, May 20-22, 1998.

[4] Johnson, W. A General Free Wake Geometry Calculation for Wings and Rotors[C]. American Helicopter Society 51st Annual Forum, Ft. Worth, TX, May 1995.

[5] Nguyen K, Johnson W. Evaluation of Dynamic Stall Models with UH-60A Airloads Flight Test Data[C]. American Helicopter Society 54th Annual Forum, Washington, D.C., May 1998.

[6] Jones H E, Kunz D L. Comprehensive Modeling of the Apache in CAMRAD II[C]. American Helicopter Society Structure Specialists' Meeting, Williamsburg, VA,October 2001.

[7] Acree C W , Jr. A CAMRAD II Model of the V-22 for Whirl-Flutter Analysis[R]. NASA TM 2004-212801, July 2004.

[8] 習(xí) 娟，吳裕平，陳平劍. 直升機(jī)旋翼系統(tǒng)載荷分析技術(shù)[J].直升機(jī)技術(shù)，2007, 3: 52-54.

[9] 王浩文，高 正．采用綜合氣彈分析方法的旋翼非定常氣動載荷計(jì)算[J].南京航空航天大學(xué)學(xué)報(bào)，2003，35(3): 268-272．

[10] 孫 韜，譚劍鋒，王浩文. 典型飛行狀態(tài)下的旋翼振動載荷計(jì)算與分析[J].南京航空航天大學(xué)學(xué)報(bào)，2011,43(3): 302-307.

[11] 吳 杰，楊衛(wèi)東，虞志浩. 旋翼槳葉結(jié)構(gòu)載荷計(jì)算方法比較研究[J].振動與沖擊，2014, 33(7): 210-214.

[12] Johnson W. Helicopter Theory[M]. Princeton: Princeton University Press, 1980.

Correlation Study for Blade Loads and Blade Structure Properties of Rotor

HU Ou, CHEN Pingjian

(China Helicopter Research and Development Institute, Jingdezhen 333001, China)

This paper presented predictions of blade structure loads using free wake model and flexible multibody dynamics method based on the CAMRAD II software. The research was focused on the problem of structure design and predictions of blade structure loads, and the relationship of blade structure properties and structure loads was studied. Two typical flight conditions were investigated including over-torque case and high speed case. The blade section properties, such as flap bending stiffness, lag bending stiffness and torsion stiffness were changed to observe how the loads changed. Meanwhile, the resonance phenomenon of blade loads was revealed.

loads； structure properties；resonance；blade

2016-09-23

胡 偶(1984-)，男，安徽池州人，博士，工程師，主要從事旋翼系統(tǒng)載荷分析、計(jì)算流體力學(xué)、網(wǎng)格生成技術(shù)等研究。

1673-1220(2017)01-008-05

V211.52

A