應(yīng)用簡化準(zhǔn)則與應(yīng)用一般方程計算載荷的對比研究

李淑萍，岳海龍，高 峰(.沈陽航空航天大學(xué) 航空航天工程學(xué)部(院)，沈陽 036；.遼寧通用航空研究院 設(shè)計與制造部，沈陽036)

應(yīng)用簡化準(zhǔn)則與應(yīng)用一般方程計算載荷的對比研究

李淑萍1,2，岳海龍1，高 峰2
(1.沈陽航空航天大學(xué) 航空航天工程學(xué)部(院)，沈陽 110136；2.遼寧通用航空研究院 設(shè)計與制造部，沈陽110136)

飛機(jī)載荷的計算關(guān)系到飛機(jī)的安全與重量。針對某雙座電動飛機(jī)的載荷，分別采用ASTM F2245-11中的簡化準(zhǔn)則以及一般方程求解的方法對載荷進(jìn)行了計算。結(jié)果表明采用一般方程求解的方法計算機(jī)翼的機(jī)動載荷，在法向方向的載荷比簡化方法求解的結(jié)果小10%以上，在切向方向上小50%以上。平尾的機(jī)動載荷在采用一般方程后的計算結(jié)果要比采用簡化方法獲得的結(jié)果小20%以上。平尾的嚴(yán)重載荷主要來自陣風(fēng)。偏航載荷的計算結(jié)果表明，兩種方法相差不多。相應(yīng)地，飛機(jī)結(jié)構(gòu)重量可以減輕10%以上。將該重量用于增加電池的重量，可以提高飛機(jī)的續(xù)航能力28%以上。

載荷；一般方程；簡化準(zhǔn)則

隨著低空空域的開放，我國進(jìn)入了通用航空的快速發(fā)展時期。RX1E飛機(jī)是全球首個獲得適航許可的雙座電動飛機(jī),該飛機(jī)也是我國自主研制的飛機(jī)，采用了流線型的外觀以及大展弦比的設(shè)計。該飛機(jī)是按照ASTM F2245-11[1]進(jìn)行設(shè)計的，該規(guī)范對于飛行載荷的計算給出了一般方法和簡化方法兩種計算方法供用戶選擇。為了對比這兩種方法，本文開展了相關(guān)研究。文獻(xiàn)[2-5]研究了現(xiàn)代飛機(jī)高速度、低重量和高機(jī)動性性能飛機(jī)連續(xù)紊流的影響，文獻(xiàn)[6-10]研究了高空、長航時的無人偵察機(jī)等的氣動彈性問題；文獻(xiàn)[11]是基于Euler方法和有限體積法計算的氣動力，耗時長，費(fèi)用高；文獻(xiàn)[12]討論了CCAR23部的載荷。對于常規(guī)布局的低速飛機(jī)來說，飛機(jī)上的載荷可以采用線性平衡方程獲得。RX1E飛機(jī)屬于600 kg以下輕型飛機(jī)，它的載荷計算是按照ASTM F2245-11中簡化原則計算的，可以不考慮全機(jī)的平衡。本文對比了ASTM F2245-11中簡化方法和列平衡方程[13]所得結(jié)果，可以給將來的設(shè)計做參考。

1 所用理論

(1)先求n=1部分的平衡載荷

飛機(jī)初始水平飛行時，發(fā)動機(jī)推力等于阻力。法向過載為1,假設(shè)推力線沒有明顯的偏轉(zhuǎn)，見圖1所示，則：

圖1 飛機(jī)的受力圖

T=D

LWB+LT-mg=0

-M0+LWBlW-LTlT-Tt=0

(1)

m:飛機(jī)質(zhì)量

g:重力加速度

LWB:飛機(jī)上的升力，不包括平尾的貢獻(xiàn)

T:發(fā)動機(jī)凈推力

D:總阻力

LT:平尾上的升力

MD:機(jī)翼的零升俯仰力矩

lW:沿機(jī)身方向，重心到機(jī)翼氣動力作用點的距離

lT:沿機(jī)身方向，重心到平尾氣動力作用點的距離

t:推力軸線到重心的距離(重心以上)

通過公式(1)求出LT，進(jìn)而求LWB，計算相應(yīng)的迎角，利用該迎角，計算更精確的尾翼升力，直到獲得收斂的解。

(2)再求Δn部分的平衡載荷

(2)

(3)

(4)

(5)

LTα=m(n-1)g(ST/S)(α1T/a1WB)[(1-dε/dα)+ρSa1WBlT/(2m)]

(6)

m:飛機(jī)質(zhì)量

Hm:縱向機(jī)動裕度

n:過載

g:重力加速度

ρ:當(dāng)?shù)卮髿饷芏?/p>

S:機(jī)翼面積

V0:飛機(jī)飛行方向的速度

a2:操縱面傾轉(zhuǎn)角產(chǎn)生的升力線斜率

ηss:舵偏角

LTSS:考慮ηss的水平安定面載荷

ST:平尾面積

a1WB:翼-身有效升力曲線斜率

Δα:迎角的變化

LWα:機(jī)翼載荷增量

LTα:平尾由于迎角變化產(chǎn)生的增量

a1T:基于平尾面積ST的平尾升力線斜率

dε/dα:機(jī)翼傾轉(zhuǎn)誘導(dǎo)的尾部下洗變化率

lT:沿機(jī)身方向，重心到平尾氣動力作用點的距離

利用公式(2)計算舵偏角，利用公式(3)計算考慮舵偏角的水平安定面載荷，利用公式(4)計算迎角的變化，利用公式(5)計算機(jī)翼載荷的增量，利用公式(6)計算平尾由于迎角變化產(chǎn)生的增量。平尾的總載荷可由LTα與LTSS的和加上初始配平載荷獲得。機(jī)翼的總載荷可由LWα加上初始配平載荷獲得。

2 RX1E飛機(jī)數(shù)據(jù)

進(jìn)行載荷計算的前提是需要給出有關(guān)RX1E飛機(jī)的參數(shù)以及載荷的工況。只有有關(guān)飛機(jī)的參數(shù)以及載荷的工況確定以后，才能開展載荷的計算研究。表1至表4的數(shù)據(jù)為有關(guān)飛機(jī)的參數(shù)以及載荷的工況。

表1 RX1E飛機(jī)的慣性參數(shù)

表2 RX1E飛機(jī)的幾何參數(shù)

表3 氣動參數(shù)

表4 載荷工況

表3中帶*的數(shù)據(jù)由于沒有風(fēng)洞吹風(fēng)數(shù)據(jù)，其中舵面偏轉(zhuǎn)引起的升力線斜率是參照文獻(xiàn)[13]中第6.2.1(1)a襟翼偏轉(zhuǎn)引起的升力線斜率進(jìn)行的，對于大舵面偏角時根據(jù)文獻(xiàn)[13]中圖6-46進(jìn)行了非線性修正。由于舵偏引起的升力系數(shù)與載荷呈線性關(guān)系，建議通過風(fēng)洞實驗獲得準(zhǔn)確值。

3 機(jī)翼對稱載荷計算對比

按照ASTM F2245-11的要求，機(jī)翼對稱載荷的最低要求見表5所示：

表5 機(jī)翼對稱載荷的最低要求

注：W為飛機(jī)最大起飛重量500 kg。

圖2為該機(jī)1∶6縮比模型在中航工業(yè)空氣動力研究院FL-8風(fēng)洞實驗所得。從圖2中可以看出，最大升力系數(shù)為1.57，最小升力系數(shù)為-0.3。因此按照氣動特性分析，過載無法達(dá)到-2,不需對比G點。

按照表5的要求，對于給定的飛機(jī)參數(shù)(參見第1節(jié))，可以得到A點和D點的載荷(見表6中的第4列)。按照文獻(xiàn)[13]中第AD1.5.3.3中提供的方法可以獲得相關(guān)載荷(見表6中第5列)，該方法是采用平衡方程進(jìn)行計算。

圖2 縱向氣動特性(Cy-α)

D點計算的載荷的負(fù)號表示向前。規(guī)范規(guī)定的為切向向后0.2 W，是比較保守的。A點計算的切向向前的載荷是比較小的，因為有阻力的抵消(零升阻力系數(shù)取0.02，K取0.0208)，規(guī)范規(guī)定的是比較保守的。機(jī)翼的剪切載荷是在此基礎(chǔ)上進(jìn)行計算的。簡化算法中沒有滾轉(zhuǎn)載荷的計算，無法進(jìn)行對比，也缺少計算所需的滾轉(zhuǎn)力矩系數(shù)。

4 平尾對稱機(jī)動載荷計算對比

對于舵面載荷，ASTM F2245-11在X1.進(jìn)行了規(guī)定，按照該規(guī)定進(jìn)行計算，所得結(jié)果如表7中第2列所示。按照文獻(xiàn)[13]中第AD1.5.3.3中提供的方法可以獲得相關(guān)載荷(見表7中第3列)，該方法是采用平衡方程進(jìn)行計算。

表6 機(jī)翼對稱載荷情況計算對比

表7 平尾對稱機(jī)動載荷比較

從上述計算中可以看出規(guī)范規(guī)定的是偏保守的，這里計算的是機(jī)動載荷。

5 平尾陣風(fēng)載荷計算

平尾的陣風(fēng)載荷也是很大的，不容忽略。按照ASTM2245-11規(guī)范第5.4.3條水平安定面的突風(fēng)載荷以及X4.1的要求計算如下表8所示。

從表8看出，平尾的嚴(yán)重工況主要來自向上的突風(fēng)載荷，在設(shè)計時需要注意。

表8 平尾突風(fēng)載荷計算

6 偏航載荷

偏航載荷屬于非對稱載荷。對于舵面載荷，ASTM F2245-11在X1.進(jìn)行了規(guī)定，按照該規(guī)定進(jìn)行計算，所得結(jié)果如表9中第2列所示。按照文獻(xiàn)[13]中第AD1.7.3.3中提供的方法可以獲得相關(guān)載荷(見表9中第3列)，該方法是采用平衡方程進(jìn)行計算。

表9 偏航機(jī)動載荷比較

在這里，A點的舵偏為20°，按照文獻(xiàn)[14]圖6-46進(jìn)行了非線性修正，K取0.7。偏航機(jī)動載荷計算的結(jié)果與規(guī)范的規(guī)定相差不大。

對于超輕型飛機(jī)的垂尾來說，陣風(fēng)載荷同樣是不可忽視的。按照ASTM F2245-11中X3.節(jié)所提供的方法，對垂尾的橫向陣風(fēng)進(jìn)行了計算(見表10)。

表10 橫向陣風(fēng)載荷

經(jīng)過計算，陣風(fēng)載荷的計算(見表10)與規(guī)定的機(jī)動載荷以及計算的結(jié)果(見表9)很相近，這主要是由于計算方法的一致性。

7 結(jié)論

采用一般方程求解的方法與采用簡化準(zhǔn)則的方法相比較，得到以下結(jié)論：

(1)在法向方向的載荷比簡化方法求解的結(jié)果小10%以上。在切向方向上小50%以上；

(2)平尾的機(jī)動載荷在采用一般方程后的計算結(jié)果要比采用簡化方法獲得的結(jié)果小20% 以上；

(3)偏航載荷的計算結(jié)果表明，兩種方法相差不多。

該飛機(jī)的結(jié)構(gòu)重量為202 kg,如果采用該方法進(jìn)行計算，飛機(jī)結(jié)構(gòu)重量可以減輕10%以上，即20 kg。該飛機(jī)的電池重量為72.5 kg,將該重量用于增加電池的重量，可以提高飛機(jī)的續(xù)航能力28%以上。此外，各個翼面的陣風(fēng)載荷是很大的，不能忽視。

[1]ASTM F2245-11 Standard Specification for Design and Performance of a Light Sport Airplane[S].

[2]嚴(yán)德，楊超，肖志鵬.彈性飛機(jī)平衡的陣風(fēng)外載荷計算與分析[J].北京航空航天大學(xué)學(xué)報，2012，38(10)：1321-1325.

[3]嚴(yán)德，楊超.基于試驗氣動力的縱向機(jī)動飛行載荷分析[J].北京航空航天大學(xué)學(xué)報，2007，33(3)：253-256.

[4]鄧立東，李天，薛曉春.飛機(jī)非線性飛行載荷計算方法研究[J].航空學(xué)報，2002，23(4)：317-320.

[5]曹奇凱.戰(zhàn)斗機(jī)機(jī)動過程與飛行載荷綜合設(shè)計[J].北京航空航天大學(xué)學(xué)報，2011,39(9)：1105-1109.

[6]孫本華.軍用飛機(jī)飛行載荷計算方法研究[J].空氣動力學(xué)學(xué)報 2006,24(2)：238-242.

[7]潘登，吳志剛，楊超，等.大柔性飛機(jī)非線性飛行載荷分析及優(yōu)化[J].航空學(xué)報，2010,31(11)：2146-2151.

[8]潘登，吳志剛，楊超.柔性飛機(jī)大變形狀態(tài)飛行載荷分析[J].北京航空航天大學(xué)學(xué)報，2010,31(10)：1158-1161.

[9]肖志鵬，萬志強(qiáng)，梁猛，等.雙尾撐布局飛機(jī)飛行載荷分析[J].北京航空航天大學(xué)學(xué)報，2008,34(4)：456-459.

[10]劉湘寧,向錦武.大展弦比柔性復(fù)合材料飛機(jī)機(jī)翼的非線性顫振分析[J].航空學(xué)報，2006，27(2)：213-218.

[11]萬志強(qiáng)，鄧立東，楊超，等.基于非線性試驗氣動力的飛機(jī)靜氣動彈性響應(yīng)分析[J].航空學(xué)報，2005,26(4)：2146-2151.

[12]王世安，尹貴魯，周喜軍.基于數(shù)值模擬的氣動載荷計算方法研究[J].空氣動力學(xué)學(xué)報,2004,22(2):160-164.

[13]丁亞修.民用飛機(jī)載荷計算研究與軟件研制[D].西安：西北工業(yè)大學(xué),2003.

[14]Denis Howe.飛機(jī)載荷與結(jié)構(gòu)布局[M].北京：航空工業(yè)出版社，2014.

[15]飛機(jī)設(shè)計手冊總編委會.飛機(jī)設(shè)計手冊第六冊 氣動手冊[M].北京：航空工業(yè)出版社，2002.

(責(zé)任編輯：吳萍 英文審校：趙歡)

Comparison of loading calculation results between simplified design load criteria and general equations

LI Shu-ping1,2,YUE Hai-long1,GAO Feng2
(1.Faculty of Aerospace Engineering,Shenyang Aerospace University,Shenyang 110136 ,China;2.Design and Manufacture Department,Liaoning General Aviation Institute,Shenyang 110136,China)

Loading calculation of airplane involves its safety and weight.The loading of electric airplane with some double-seat was calculated using both simplified design load criteria(ASTM F2245-11)and general equations.The results show the normal load of wing resulting from general equations decreased by more than 10% comparing the simplified criteria,the tangential load reduced by more than 50%,but loading of horizontal wing decreased by more than 20%.The overload of horizontal wing came from gust mainly.The results from loading calculation of yaw indicated the two calculation methods were similar.Consequently,the structure weight of airplane decreased by more than 10%.If the reduced weight is substituted by batteries,the flight endurance of airplane will increase by more than 28%.

loading;general equations;simplified criteria

2016-10-28

李淑萍(1968-)，女，河北河間人，副教授，主要研究方向：飛機(jī)結(jié)構(gòu)設(shè)計，E-mail:lsp@sau.edu.cn。

2095-1248(2017)02-0033-05

V211.47

A

10.3969/j.issn.2095-1248.2017.02.006