升力體高超聲速飛行器非定常滾轉力矩建模研究

高 清，李建華，李 潛

(中國航天空氣動力技術研究院，北京 100074)

升力體高超聲速飛行器非定常滾轉力矩建模研究

高 清，李建華*，李 潛

(中國航天空氣動力技術研究院，北京 100074)

為了深入研究升力體高超聲速飛行器相對薄弱的橫側向穩定性問題，對滾轉自由振蕩風洞試驗數據進行了譜分析，建立了多個氣動頻率余弦和形式的非定常滾轉力矩模型。由于升力體高超聲速飛行器的滾轉自由振蕩曲線呈現非定常、非線性和一定的周期性特征，且對試驗結果的譜分析發現，在多種氣動狀態下，都存在除機械振動頻率外的三個振動頻率，將滾轉力矩表達為此三個頻率余弦和的形式。該滾轉非定常氣動力矩數學模型捕捉了試驗的基本趨勢涵蓋了其主要的量值范圍，反映了升力體高超聲速飛行器橫向流場擾流的多尺度和周期性特征。

升力體;高超聲速;滾轉力矩;非定常氣動力;建模

0 引 言

飛行器空氣動力學問題的復雜性之一就在于繞流結構的多尺度特征。飛行器繞流易出現分離和分離渦的周期性脫落。繞飛行器的來流即使是定常，也會因分離渦的周期性脫落使流動呈現非定常和非對稱特征。在高超聲速條件下飛行器繞流的周期性、非定常過程更加強烈，剪切流和渦結構由于軸向動量的增大被拉長，此時橫向擾動極易使流動的軸向動能轉化為橫側向動能。

通常在亞跨和低超聲速，繞流中低頻的大尺度結構能量遠大于高頻的小尺度結構能量，使流體動力學參數通常表現為單頻、線性。但對高超聲速，隨馬赫數的增加，繞流中高頻、小尺度結構的能量及其對飛行器的影響越來越大，流體動力學參數無法繼續由線性系數描述。小尺度流動結構的重要特征是流動的三維性，高超聲速的橫流擾動強烈，尤其是滾轉力矩呈現強烈的非定常和非線性特征，必須尋求一種新的參數形式，進行非定常氣動力建模。

根據泛函分析理論，廣義非定常氣動力可寫為Duhamel積分形式［1］:

它表示飛行器繞流引起的氣動力和力矩不僅與流動變量的瞬時值有關，且與其變化歷程相關，但要確定機動飛行器上瞬時氣動載荷和運動變量間的關系非常困難。為此，人們提出了各種假定和近似，以期找到既符合基本原理又適用于具體實踐的數學表達式。

非定?？諝鈩恿Ｗ?0世紀初在國外開始發展以來，陸續出現了多種非定?？諝鈩恿Φ慕７椒ǎ瑲w納起來大致有:積分模型［2-5］、微分模型［6-9］、狀態空間模型［10-11］、神經網絡模糊邏輯模型［10-19］和代數模型［20-22］。積分模型在理論上是完備的，但過于復雜，且模型中所含未知函數很難通過數值計算、風洞試驗和飛行試驗確定，因此很難應用;微分方程模型存在非線性項物理意義不直觀、非線性項近似表達式中的人為加權函數不易確定等缺點;代數模型基于準定常假設，是最簡單、最方便的非線性模型。常用的代數模型有級數函數模型和樣條函數模型等。采用代數模型時，須根據物理問題選擇適當的狀態參數、階次及組合項，用最少的項來逼近物理本質。

研究通過對升力體高超聲速飛行器非定常氣動力風洞試驗數據的分析，試圖建立客觀直接、有明確物理意義、并具有一定工程實用指導意義的滾轉力矩數學模型。

1 建模原理

研究基于風洞試驗獲取的非定常氣動力數據，根據數據的頻譜特性，分析建立一種級數形式的非定常氣動力表達式。級數的項數、各項的系數、頻率由風洞試驗結果確定。

由于本研究中升力體高超聲速飛行器的滾轉自由振動曲線與Hopf分叉后的振蕩曲線類似(圖1)，所以借鑒了Hopf分叉［23-25］后氣動力的表達式，將滾轉力矩表達為振動頻率余弦和的形式。

圖1 Hopf分叉Fig.1 Hopf bifurcation

2 建模方法

對滾轉非定常氣動力試驗數據的譜分析發現，滾轉方向的氣動力矩存在多個頻譜，流場特征可近似由幾個主要振動頻率共同描述。將滾轉力矩表達為mx的形式，其中mx表示滾轉力矩，an表示振幅，ωn表示根據試驗結果的頻譜分析得到的主要振動頻率，λn表示相位。

建模過程:首先，對升力體高超聲速飛行器模型的地面無風自由滾轉振動數據進行譜分析(圖2)，獲得機械阻尼的頻率(10.74 Hz)。然后對四種狀態下的滾轉非定常氣動力試驗結果進行譜分析(圖3)，確定主要的振動頻率及其幅值和相位，發現在這些狀態下，試驗曲線都存在除機械阻尼頻率外的三個振動頻率，按振動頻率對應的幅值從大到小的順序依次為:10.25 Hz、9.77 Hz、11.23 Hz。而且根據頻譜圖，機械阻尼對應頻率10.74 Hz的幅值小于三個氣動阻尼頻率對應的幅值，說明在研究的氣動狀態下，氣動阻尼在總阻尼中所占的比例較高，機械阻尼干擾相對較小。

圖2 地面無風自由滾轉振動數據及其譜分析圖Fig.2 Ground-tested roll-free-oscillation result and its spectrum analysis

圖3 滾轉非定常氣動力試驗數據及其譜分析圖Fig.3 Wind-tunnel-tested roll-free-oscillation result and its spectrum analysis

滾轉力矩可近似由三個振動頻率的余弦函數表達為:

其中ai表示氣動主頻對應的幅值，λi表示相位。對于不同的氣動狀態，三個氣動主頻對應的幅值和相位稍有不同(圖4)。

由圖4可見，在試驗研究的四種狀態下，滾轉力矩模型中頻率項幅值存在較明顯的規律性，頻率項幅值隨馬赫數的增加而減小，隨攻角的增加而略有增加，相位角隨馬赫數和攻角的變化也呈現一定的規律性。

本節通過對非定常氣動力試驗數據的譜分析，獲得了模型系統的三個振動頻率、幅值和相位，建立了這三個氣動振動頻率余弦和形式的滾轉力矩數學模型。

3 建模模型的物理意義

根據動態試驗相似條件，試驗模型振動頻率的無量綱值與真實飛行器在大氣中的無量綱運動頻率應相同［26］。因為只有兩者相同時，真實飛行器和試驗模型的繞流才是相似的。以Ma=5，飛行高度30 km配平攻角α=1.2°狀態為例，真實飛行器在大氣中的滾轉振動的無量綱頻率可表示為:

圖4 主要振動頻率及相位幅值Fig.4 Main vibration frequencies and their phases and amplitudes

其中L表示試驗模型的參考長度。

試驗模型的無量綱振動頻率可表示為:

根據動態試驗相似關系，有:

即，橫向參考長度和試驗模型的振動頻率之間存在式(6)的線性比例關系。對非定常氣動力風洞試驗數據的譜分析結果表明，模型的振動曲線存在除機械振動頻率外的三個較為明顯的振動頻率，對應飛行器橫向繞流有三個特征尺度，即橫向分離或轉捩流動有三個不同的尺度。這三個特征長度可根據式(6)由試驗振動頻率直接獲得(表1)。

表1 試驗振動頻率與特征長度Table 1 Main vibration frequencies and their characteristic lengths

如選取飛行器的橫向跨度為特征長度，則該特征長度在飛行器上的分布見圖5，流場特性在這三個位置上可能出現非對稱分離、轉捩等變化。

圖5 橫向參考長度在飛行器上的位置Fig.5 Distribution of lateral reference lengths on aircraft

從流場角度分析，在高超聲速飛行時，飛行器繞流中湍流與轉捩等小尺度流動的貢獻明顯增加，不同尺度流動的頻率隨馬赫數增加也越來越近，高超聲速繞流的多尺度效應是高超聲速流動有別于低速流動最重要的特征，也是高超聲速飛行器氣動力出現強烈非線性的根本原因之一。滾轉力矩多頻、非線性特征使高超聲速飛行與控制面臨較大困難。

4 模型的有效性

為了驗證建立的三個振動頻率余弦和形式滾轉力矩模型模擬試驗數據的效果，將各狀態非定常滾轉力矩數學模型數據與相應試驗數據進行對比，如果趨勢和量值基本一致，說明這三個氣動振動主頻基本可代表該氣動狀態主要特征，反之，則可能是氣動振動主頻選擇不合適，或三個氣動振動主頻不足以描述該狀態的氣動特征。

圖6 三個氣動主頻模擬數據與相應試驗數據的對比Fig.6 Comparison of the model of roll moment constructed in this paper with the corresponding aerodynamic test data

圖6給出Ma=5、α=0°狀態下，建立的滾轉力矩模型數據與試驗數據的對比?？梢?，滾轉力矩數學模型基本捕捉了試驗曲線的基本趨勢，也涵蓋了主要的量值范圍，而且周期性的表達式反映了橫向力矩周期性的物理特征，可初步認定，采用三個振動主頻余弦和模擬滾轉力矩是適當的。

5 結 論

升力體高超聲速飛行器滾轉自由振蕩試驗發現，模型的振動具有非線性、非定常、多頻、周期性特征，飛行器的滾轉力矩已不能由線性系數繼續描述。對該試驗模型多種氣動狀態試驗數據的譜分析表明，存在除機械振動頻率外的三個振動頻率。此三個振動頻率對應飛行器橫向繞流有三個特征尺度，即橫向分離或轉捩流動有三個不同的尺度。研究表明，多尺度特征是高超聲速繞流有別于低速繞流最重要的特征，也是高超聲速飛行器氣動力出現強烈非線性的根本原因之一。

建立了這三個頻率的余弦和形式的滾轉非定常氣動力模型。從滾轉力矩建模值與相應氣動數據的對比來看，數學模型捕捉了試驗曲線的基本趨勢，也涵蓋了主要的量值范圍，而且周期性的表達形式反映了橫向力矩周期性的物理特征。初步認為本文建立的非定常滾轉力矩數學模型是適當的。

［1］ Nixon David，Rodman Laura C.A study of supermanuever aerodynamics［R］.ADA 218378，1990.

［2］ Tobak M，Chapman G T，Schiff L B.Mathematical modeling of the aerodynamic characteristics in flight dynamics［R］.NASA TM-85880，1985.

［3］ Tobak M，Chapman G T，Uenal A.Modeling aerodynamic discontinuities and the onset of chaos in flight dynamical systems［R］.NASA TM-89420，1986.

［4］ Gupta N K，Iliff K W.Identification of integro-differential systems for application to unsteady aerodynamics and aeroelasticity［R］.AIAA 85-1763，1985.

［5］ Gupta N K，Iliff K W.Identification of unsteady aerodynamics and aeroelastic integro-differential systems［R］.NASA TM-86749，1985.

［6］ Bhaskar T G，Devi J V.Nonuniform stability and boundedness criteria for set differential equations［J］.Applicable Analysis，2005，84 (2):131-142.

［7］ Bhaskar T G，Devi J V.Stability criteria for set differential equations［J］.Mathematical and Computer Modelling，2005，41(1):1371-1378.

［8］ Lakshmikantham V，Bhaskar T G，Devi J V.Theory of set differential equations in metric spaces［M］.Cambridge Scientific Publisher，2006:56-122.

［9］ Feng Pan，Tao Cai，Jie Chen，et al.The stability analysis of particle swarm optimization without lipschitz condition constraint［J］.Control Theory and Application，2004，1(2):46-57.

［10］Shi Zhiwei，Ni Fangyuan，Chen Yongliang，The state-space models based on two-step linear regression method［J］.Acta Aerodynamica Sinica，2013，31(6):699-703.(in Chinese)

史志偉，倪芳原，陳永亮.基于兩步線性回歸的狀態空間模型建立與驗證［J］.空氣動力學學報，2013，31(6):699-703.

［11］Sun Haisheng，Zhang Haiyou，Liu Zhitao.Comparative evaluation of unsteady aerodynamics modeling approaches at high angle of attack［J］.Acta Aerodynamica Sinica，2011，29(6):733-737.(in Chinese)

孫海生，張海酉，劉志濤.大迎角非定常氣動力建模方法研究［J］.空氣動力學學報，2011，29(6):733-737.

［12］Yang Meng，Huang Da.The effects of variables on unsteady aerodynamic modeling［J］.Acta Aerodynamica Sinica，2011，29(3): 355-359.(in Chinese)

楊勐，黃達，建模變量對非定?？諝鈩恿祵W模型的影響［J］.空氣動力學學報，2011，29(3):355-359.

［13］Shi Zhiwei，Wang Zhenghua，Li Juncheng.The research of RBFNN in modeling of nonlinear unsteady aerodynamics［J］.Acta Aerodynamica Sinica，2012，30(1):108-112.(in Chinese)

史志偉，王崢華，李俊成.徑向基神經網絡在非線性非定常氣動力建模中的應用研究［J］.空氣動力學學報，2012，30(1):108-112.

［14］Shi Zhiwei，Wu Genxing，Huang Da.The validation of unsteady aerodynamic model based on the high amplitude harmonic rolling motion test［J］.Acta Aerodynamica Sinica，2010，28(6):650-654.(in Chinese)

史志偉，吳根興，黃達，基于大振幅諧波運動的非定常氣動模型風洞實驗驗證［J］.空氣動力學學報，2010，28(6):650-654.

［15］Kumar R，Ganguli R，Omkar S N，et al.Rotorcraft parameter estimation from real time flight data using radial basis function networks［J］.Journal of Aircraft，2008，45(1):333-341.

［16］Xia Y S，Feng G，Wang J.A novel recurrent neural network for solving nonlinear optimization problems with inequality constraints［J］.IEEE Trans on Neural Networks，2008，19(8):1340-1353.

［17］Peyada N K，Ghosh A K.Longitudinal aerodynamic parameter estimation using neural network and Gauss Newton method［J］.Journal of Aerospace Sciences and Technology，2009，61(2):295-304.

［18］Singh S，Ghosh A K.Parameter estimation from flight data of a missile using maximum likelihood and neural network method［C］// 2006 AIAA Flight Mechanics Conference and Exhibit，2006:21-24.

［19］ Malmathanraj R，Thamarai Selvi S E.Prediction of aerodynamic characteristics using neural networks［J］.Asian Journal of Information Technology，2008，7(1):19-26.

［20］蔡金獅.動力學系統辨識與建模［M］.北京:國防工業出版社，1991.

［21］Klein V，Batterson J G，Murphy P C.Determination of airplane model structure from flight data by using modified stepwise regression［R］.NASA TP-1916，1981.

［22］Klein V，Batterson J G.Determination of airplane model structure from flight data using splines and stepwise regression［R］.NASA TP-2126，1983.

［23］Librescu L，Marzocca P，Silva W A，et al.Supersonic/hypersonic flutter and postflutter of geometrically imperfect circular cylindrical panels［J］.Journal of Spacecraft and Rockets，2002，39(5):802-812.

［24］Chaitanya V S K.Stability analysis of structurally unstable man-machine system involving time delays［J］.Nonlinear Analysis:Real Word Applications，2005，6(5):845-857.

［25］Hassard B，Kazarinoff N，Wan Y H.Theory and application of hopf bifurcation［M］.Cambridge:Cambridge University Press，1981: 14-35.

［26］李周復，主編.風洞特種試驗技術［M］.北京:航空工業出版社，2010.

Non-steady roll moment modeling for hypersonic lifting configuration

Gao Qing，Li Jianhua*，Li Qian
(China Academy of Aerospace Aerodynamic，Beijing 100074，China)

Thorough study of lateral stability of hypersonic lifting aircrafts needs exact expression of aerodynamic roll moment.Spectrum analysis is carried out using roll free-oscillation wind-tunnel test data，and a mathematical model of roll moment is proposed as an expression by the sum of the cosine of the main frequencies from spectrum analysis.The roll test curves of hypersonic lifting model present non-linearity，non-steady and periodical characteristics，and the results of spectrum analysis of roll vibration data under several different flow conditions show that they share three main vibration frequencies apart from the mechanical vibration frequency，which indicates that there are three scales of lateral separation or transition in the flow field.Comparison of the mathematical model of roll moment constructed with the corresponding aerodynamic test data indicates，this model can capture the primary tendency of test curves，and cover the main magnitude domain of roll aerodynamic moment.The mathematical model gives prominence to the multi-scale and periodical characteristics of lateral flow field of hypersonic lifting model.

lifting-configuration;hypersonic;roll moment;non-steady aerodynamic;aerodynamic modeling

V211.3;O177.7

A

10.7638/kqdlxxb-2013.0061

0258-1825(2015)03-0392-05

2013-06-04;

2013-10-16

高清(1983-)，女，山東曹縣人，博士，高級工程師，研究方向:高超聲速空氣動力學和飛行力學.E-mail:1034505942@qq.com

李建華*(1977-)，男，江西廣昌人，高級工程師，研究方向:飛行器總體設計和氣動布局設計.E-mail:ssanliren@126.com

高清，李建華，李潛.升力體高超聲速飛行器非定常滾轉力矩建模研究［J］.空氣動力學學報，2015，33(3):392-396.

10.7638/kqdlxxb-2013.0061 Gao Q，Li J H，Li Q.Non-steady roll moment modeling for hypersonic lifting configuration［J］.Acta Aerodynamica Sinica，2015，33(3):392-396.