王澤峰 ，李清龍

中國飛行試驗研究院，陜西 西安 710089

俯沖拉起是直升機(jī)典型的飛行譜設(shè)計動作之一，在實際使用中具有非常重要的戰(zhàn)術(shù)意義。在型號定型試飛過程中，法向過載包線驗證試飛屬于必須完成的內(nèi)容，也是通過直升機(jī)俯沖拉起試飛來進(jìn)行的。在執(zhí)行該動作過程中，直升機(jī)可能會出現(xiàn)極限的法向過載，此時旋翼槳轂系統(tǒng)部件尤其是直升機(jī)槳葉，都會承受大的載荷，嚴(yán)重者可能會出現(xiàn)結(jié)構(gòu)損壞的情況。因此，無論從直升機(jī)旋翼部件結(jié)構(gòu)強(qiáng)度的角度，還是從直升機(jī)試飛安全的角度考慮，都非常有必要研究直升機(jī)俯沖拉起過程中旋翼載荷的變化規(guī)律。

本文基于某型直升機(jī)旋翼載荷試飛，實測了直升機(jī)主槳葉在俯沖拉起試飛科目中的揮舞、擺振彎矩，分析其變化規(guī)律。

1 理論分析

以主槳葉揮舞為例進(jìn)行說明。圖1所示為主槳葉運(yùn)動過程中的受力分析[1-6]。由圖中可以看出，槳葉微段在運(yùn)動過程中主要受升力ΔT，重力ΔG以及離心力ΔF。三個力均對揮舞方向構(gòu)成彎矩，最終在三個方向力的合成作用下，反映在結(jié)構(gòu)上就是主槳葉揮舞方向的彎矩Mf，即：

式中：升力構(gòu)成的彎矩為重力構(gòu)成的彎矩為離心力構(gòu)成的彎矩為

圖1 槳葉揮舞方向受力分析Fig.1 A typical helicopter dive and pull-up motion

對于MG、MF而言，主要受質(zhì)量分布和旋翼轉(zhuǎn)速的影響，升力彎矩項的影響因素與槳葉升力的影響因素相關(guān)，可由葉素理論進(jìn)行分析。葉素示意圖如圖2所示。

圖2 直升機(jī)在垂直下降時的葉素分析Fig.2 Blade model while helicopter vertical drop

圖2中，V為直升機(jī)垂向速度；v1為旋翼誘導(dǎo)速度；φ為槳葉安裝角；ε為來流角；α為槳葉迎角?？紤]俯沖拉起動作過程中的旋翼迎角αs及前飛速度V0的影響，在槳葉方位角ψ處的槳葉上，徑向r處的相對氣流速度為：

式中：Vr為軸向相對氣流速度分量，Vc為徑向相對氣流分量。

又槳葉迎角

由升力計算公式可得：

由式（3）可以看出，在直升機(jī)俯沖拉起動作過程中，槳葉升力主要受直升機(jī)前飛速度、旋翼迎角、誘導(dǎo)速度、旋翼轉(zhuǎn)速等因素的影響，而在整個俯沖拉起動作過程中，前飛速度、旋翼迎角、誘導(dǎo)速度、旋翼轉(zhuǎn)速均為變量。因此，主槳葉升力也是一個動態(tài)變化過程，升力構(gòu)成的彎矩也是動態(tài)變化的，最終體現(xiàn)在主槳葉上的結(jié)構(gòu)揮舞彎矩[7～9]載荷值上也是動態(tài)變化的，且同樣受上述因素的影響。擺振方向彎矩結(jié)果類似。

2 試飛方法

依據(jù)GJB720.7A中3.2.10條款，直升機(jī)俯沖拉起試飛方法為：試驗機(jī)應(yīng)在要求的高度上，在不同重量、重心位置、旋翼轉(zhuǎn)速組合下，以不同速度進(jìn)行穩(wěn)定平飛、下降或俯沖，在規(guī)定的時間內(nèi)使操縱機(jī)構(gòu)勻速移動，直到產(chǎn)生規(guī)定的法向過載。在達(dá)到規(guī)定的過載值后，在規(guī)定的時間內(nèi)勻速返回，直到平飛所要求的位置。

以某直升機(jī)載荷譜試飛為例，由于該型機(jī)旋翼轉(zhuǎn)速為恒定轉(zhuǎn)速，故旋翼轉(zhuǎn)速不作變量考慮。在試驗機(jī)選取起飛重量m1、m2兩個重量，每個重量選取縱向前重心、縱向正常重心、縱向后重心，共計6種構(gòu)型，每種構(gòu)型下選取低、中、高（壓力高度1000m，2000m，3000m）三個高度，以規(guī)定表速進(jìn)行俯沖拉起。典型的俯沖拉起動作[10]過程如圖3所示?？梢钥闯觯囷w動作分為以下幾個步驟：穩(wěn)定平飛段→進(jìn)入俯沖段→穩(wěn)定俯沖段→拉起段→改出平飛段。其中，在該型機(jī)試飛過程中，由進(jìn)入俯沖段開始到改出平飛段動作結(jié)束，直升機(jī)總距全程固定，通過縱向推拉桿和航向修正來完成該俯沖拉起動作。

圖3 俯沖拉起典型動作示意圖Fig.3 A typical helicopter dive and pull-up motion

3 主槳葉載荷實測與分析

3.1 主槳葉載荷實測

主槳葉作為主要升力面和操縱面，在實際飛行中主要承受揮舞彎矩和擺振彎矩，根據(jù)主槳葉靜力計算和試驗結(jié)果，確定強(qiáng)度較弱的剖面。本次主槳葉實測載荷，選取三個關(guān)鍵剖面（分別為距離槳葉根部110mm、1600mm、3639mm三個剖面），采用電阻應(yīng)變計法進(jìn)行揮舞彎矩和擺振彎矩的飛行實測。

由圖4、圖5可以看出，槳葉110mm、1600mm、3639mm剖面實際承受的是交變載荷，周期性強(qiáng)，從頻譜分析結(jié)果可以看出，揮舞彎矩和擺振彎矩的基頻為4.3Hz，等于該型機(jī)旋翼轉(zhuǎn)速，其他各階頻率為基頻的整數(shù)倍。從動力學(xué)角度考慮，符合槳葉揮舞運(yùn)動和擺振規(guī)律。

圖4 揮舞彎矩時間歷程Fig.4 Time history curve of fl apping moment

圖5 1600mm剖面擺振彎矩頻譜分析Fig.5 Spectrum analysis of 1600mm position lagging moment

3.2 載荷分析

從交變載荷的角度入手，統(tǒng)計各動作段的動載荷，對其進(jìn)行規(guī)律分析。

式中：Smax、Smin為旋翼轉(zhuǎn)一圈內(nèi)的載荷最大值、最小值，Sd為旋翼轉(zhuǎn)一圈的動載荷。

式中：Sdi為旋轉(zhuǎn)第i圈內(nèi)的動載荷，N為動作段旋轉(zhuǎn)圈數(shù)，S為動作段內(nèi)的動載荷均值。

選取主槳葉1600mm剖面揮舞、擺振彎矩，繪制俯沖拉起動作全程時間歷程曲線，如圖6和圖7所示。

圖6 槳葉1600剖面揮舞彎矩俯沖拉起全程變化趨勢Fig.6The whole change trend of 1600mm position flapping moment while helicopter dive and pull-up

由圖6和圖7可以看出，在轉(zhuǎn)入俯沖的過程中，主槳葉揮舞、擺振彎矩動載荷都有減小趨勢；在穩(wěn)定俯沖過程中，主槳葉揮舞彎矩動載荷有明顯增加趨勢；在拉起段過程中，主槳葉揮舞、擺振彎矩的動載荷增大，基本達(dá)到整個俯沖拉起過程中的最大值；在改平飛過程中，主槳葉揮舞、擺振彎矩動載荷減小。即穩(wěn)定俯沖和拉起段時為整個俯沖拉起過程中槳葉的嚴(yán)重受載段。

圖7 槳葉1600剖面擺振彎矩俯沖拉起全程變化趨勢Fig.7 The whole change trend of 1600mm position lagging moment while helicopter dive and pull-up

這主要是由于在轉(zhuǎn)入俯沖過程中，直升機(jī)低頭，旋翼槳盤前傾，旋翼升力減小，故動載荷減?。辉诜€(wěn)定俯沖過程中，槳葉迎角一定的情況下，直升機(jī)速度增大，槳葉相對來流速度也增大，旋翼升力增大，故動載荷增大；在拉起過程中，槳葉迎角增大，旋翼升力增大，達(dá)到最大值，故動載荷也相應(yīng)增大，達(dá)到最大值；在改平飛過程中，旋翼槳盤再次前傾，旋翼升力減小，故動載荷再次減小。

統(tǒng)計不同法向過載下的主槳葉揮舞、擺振彎矩動載荷值，如圖8、圖9所示。

從圖8、圖9可以看出，在同一重量、重心構(gòu)型下，110mm、1600mm兩個剖面的實測揮舞彎矩、擺振彎矩動載荷值隨法向過載的增大，也呈增大趨勢。這主要是由于大的法向過載意味著大的旋翼拉力，也即大的主槳葉升力，故動載荷值也增大。在包線范圍內(nèi)的最大法向過載俯沖拉起狀態(tài)下，槳葉揮擺彎矩實測動載荷值接近最大限制值，即1h疲勞限制值。

圖8 同一構(gòu)型下，110mm剖面實測動載荷值隨過載的變化Fig.8Under the same configuration， the dynamic load measured in the 110mm section varies with the overload

圖9 同一構(gòu)型下，1600mm剖面實測動載荷值隨過載的變化Fig.9Under the same configuration， the dynamic load measured in the 1600mm section varies with the overload

4 結(jié)論

通過飛行試驗，可以得到以下結(jié)論：

（1）在整個俯沖拉起動作過程中，主槳葉在穩(wěn)定俯沖和拉起段時受載嚴(yán)重，在試飛過程中，應(yīng)重點(diǎn)監(jiān)控這兩個時間段的槳葉載荷值。

（2）主槳葉剖面的揮舞、擺振彎矩動載荷值隨法向過載增加而增大。對于直升機(jī)大法向過載的俯沖拉起應(yīng)關(guān)注槳葉的揮舞、擺振彎矩。

[1] Ormiston R A. Comparison of several methods for predicting loads on hypothetical helicopter rotor [C]// AHS/NASA-Ames Specialist’s Meeting on Rotorcraft Dynamics， 1974.

[2] Charies B D， Hassan A A. A correlation study of rotor inf l ow in forward fl ight [C]// The 44th AHS Forum， 1988.

[3] Huber H， Polz G. Studies on blade to blade and Rotor-fuselage tail interferences [C]// AGARD Conference Proceeding， 1982.

[4] Johson W. Helicopter theory [M]. Princeton: Princeton University Press， 1980.

[5] Hansford R E. A unified formulation of rotor load prediction Methods [C]// The 41th AHS Forum，1985.

[6] Landgrebe A J. An analytical method for predicting rotor wake geometry [J]. JAHS， 1969， 14（4）：20-32.

[7] 朱世晉，李南慧.直升機(jī)旋翼氣動載荷的計算與分析 [J].南京航空航天大學(xué)學(xué)報，1979（3）：103-118.ZHU Shijin， LI Nanhui. The analysis and calculation of aerodynamic load of helicopter rotor [J]. Journal of Nanjing University of Aeronautics and Astronautics， 1979（3）：103-118.（in Chinese）

[8] David J. Statistical analysis of structural flight test data [R].A.H.S， 1988.

[9] 孫之釗，蕭秋庭，徐桂祺. 直升機(jī)強(qiáng)度 [M].北京：航空工業(yè)出版社，1990.SUN Zhizhao， XIAO Qiuting， XU Guiqi. Helicopter structure strength [M]. Beijing： Aviation Industry Press， 1990. （in Chinese）

[10] Robert T N， Chen J. Kinematic properties of the helicopter in coordinated turns [R]. NASA， 1973.