• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于數(shù)值虛擬飛行技術(shù)的飛行器動(dòng)態(tài)特性分析

    2016-11-14 00:58:02黃宇閻超席柯王文
    航空學(xué)報(bào) 2016年8期
    關(guān)鍵詞:飛行技術(shù)返回艙迎角

    黃宇, 閻超,*, 席柯, 王文

    1.北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院, 北京 100083 2.中國兵器工業(yè)導(dǎo)航與控制技術(shù)研究所, 北京 100089

    ?

    基于數(shù)值虛擬飛行技術(shù)的飛行器動(dòng)態(tài)特性分析

    黃宇1, 閻超1,*, 席柯2, 王文1

    1.北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院, 北京100083 2.中國兵器工業(yè)導(dǎo)航與控制技術(shù)研究所, 北京100089

    基于結(jié)合結(jié)構(gòu)重疊網(wǎng)格、閉環(huán)PID控制器、舵偏控制律、剛體六自由度運(yùn)動(dòng)和非定常Navier-Stokes方程求解等模塊的數(shù)值虛擬飛行技術(shù),對“起源號”返回艙、基本帶翼導(dǎo)彈外形的多自由度非定常運(yùn)動(dòng)、受控特性及控制參數(shù)的整定開展了模擬。分析了不同自由度(DOF)下飛行器的運(yùn)動(dòng)特性,飛行器受擾動(dòng)后的穩(wěn)定性及控制參數(shù)的整定。計(jì)算結(jié)果表明:利用數(shù)值虛擬飛行技術(shù)可有效地開展復(fù)雜流動(dòng)下飛行物體非線性運(yùn)動(dòng)問題的研究,對研究飛行器在非線性流動(dòng)下的動(dòng)態(tài)特性、受控特性、流動(dòng)機(jī)理研究以及控制律的設(shè)計(jì)檢驗(yàn)具有較高的工程價(jià)值和應(yīng)用前景。

    計(jì)算流體力學(xué)(CFD); 重疊網(wǎng)格; PID控制; 虛擬飛行; 動(dòng)導(dǎo)數(shù); 控制參數(shù)整定; 動(dòng)態(tài)特性; 非定常流動(dòng)

    飛行器在大迎角、過失速等狀態(tài)下飛行時(shí),會(huì)出現(xiàn)復(fù)雜的分離渦等強(qiáng)非線性、非定常流動(dòng),與飛行器運(yùn)動(dòng)耦合將會(huì)使飛行器產(chǎn)生一系列的非定常運(yùn)動(dòng)現(xiàn)象,如極限環(huán)搖滾、俯仰振蕩等,從而引起飛行器氣動(dòng)力的突變、遲滯及非對稱特性,極易導(dǎo)致飛行器運(yùn)動(dòng)進(jìn)入發(fā)散狀態(tài),嚴(yán)重影響飛行器的安定性、穩(wěn)定性及受控性。由此,現(xiàn)代先進(jìn)空戰(zhàn)武器設(shè)計(jì)對飛行器的動(dòng)態(tài)特性、穩(wěn)定性及控制效率等的預(yù)測技術(shù)提出了要求[1]。

    飛行器出現(xiàn)非定常、非線性流動(dòng)情況下,線性化數(shù)據(jù)不能正確地刻畫飛行器非線性的氣動(dòng)力和運(yùn)動(dòng)規(guī)律[2-3]。為了能有效地研究飛行器非定常運(yùn)動(dòng)中的氣動(dòng)、控制等問題,國內(nèi)外研究人員從20世紀(jì)80年代開始發(fā)展了一種新型的風(fēng)洞試驗(yàn)技術(shù),即風(fēng)洞虛擬飛行技術(shù)[4-6](Wind Tunnel Based Virtual Flight Testing)。通過將實(shí)驗(yàn)縮比模型安裝有若干個(gè)自由度的支撐機(jī)構(gòu)上,在風(fēng)洞中開展運(yùn)動(dòng)模型的動(dòng)力特性、控制特性的綜合實(shí)驗(yàn),在一定程度上模擬飛行器的真實(shí)飛行狀態(tài),獲得飛行器動(dòng)態(tài)品質(zhì)及控制特性。但由于風(fēng)洞試驗(yàn)的固有缺陷,如支架/洞壁干擾、機(jī)械阻尼、相似準(zhǔn)則等難以完全滿足等情況,導(dǎo)致難以準(zhǔn)確地測量動(dòng)態(tài)參數(shù)。此外,由于風(fēng)洞一般只能提供直勻流,且模型支撐機(jī)構(gòu)并不能在大范圍內(nèi)實(shí)現(xiàn)真實(shí)飛行自由度,因而模擬的仿真程度也受到制約。

    近年來,計(jì)算流體力學(xué)(Computational Fluid Dynamics, CFD)及大規(guī)模高性能計(jì)算技術(shù)發(fā)展迅速,CFD在基礎(chǔ)理論和工程實(shí)踐上都取得了令人矚目的成就,以CFD技術(shù)為基礎(chǔ)實(shí)現(xiàn)的基于方程的數(shù)值虛擬飛行技術(shù)已經(jīng)成為可能[7-8]。數(shù)值虛擬飛行技術(shù)將飛行動(dòng)力學(xué)方程/流體力學(xué)方程/控制律耦合實(shí)時(shí)求解,可開展全尺寸,真實(shí)運(yùn)動(dòng)/氣流參數(shù)的非定常模擬,避免了風(fēng)洞虛擬飛行技術(shù)的固有缺點(diǎn),運(yùn)動(dòng)仿真程度更高。此外數(shù)值虛擬飛行技術(shù)靈活性更好、耗費(fèi)低、周期短,在機(jī)理研究中對流場、運(yùn)動(dòng)參數(shù)的獲取更為方便和豐富,在實(shí)際工程應(yīng)用中可在各個(gè)階段介入研究,為設(shè)計(jì)及驗(yàn)證飛行器氣動(dòng)、控制等系統(tǒng)的特性提供有力的依據(jù),填補(bǔ)風(fēng)洞虛擬飛行及飛行試驗(yàn)的不足[9]。

    本文基于數(shù)值虛擬飛行技術(shù),開展了返回艙、帶翼導(dǎo)彈外形的動(dòng)態(tài)特性研究。研究結(jié)果表明:數(shù)值虛擬飛行技術(shù)對復(fù)雜運(yùn)動(dòng)下的飛行器動(dòng)態(tài)特性、飛行器的操縱特性、控制系統(tǒng)參數(shù)的整定等問題的研究具有較高的工程價(jià)值和應(yīng)用前景。

    1 數(shù)值虛擬飛行方法

    1.1流動(dòng)主控方程

    主控方程為非定常、可壓縮、雷諾平均的三維Navier-Stokes方程。采用來流參數(shù)為無量綱特征變量,對方程無量綱化,并將方程轉(zhuǎn)換到一般曲線坐標(biāo)系下,可得如下形式[8]:

    (1)

    式中:Q為守恒變量向量;F、G和H為對流通量;Fv、Gv和Hv為黏性通量;Re為來流雷諾數(shù)。

    采用基于結(jié)構(gòu)網(wǎng)格的有限體積法求解,對流通量計(jì)算采用帶熵修正的Roe格式,采用MUSCL插值和vanAlbada限制器達(dá)到二階空間精度。時(shí)間推進(jìn)格式在定常計(jì)算時(shí)選用隱式LU-SGS方法,在非定常計(jì)算時(shí)選用雙時(shí)間步LU-SGS格式。采用一方程的Spalart-Allmaras湍流模型計(jì)算獲得湍流渦黏性。利用基于共享內(nèi)存的Open-MP并行技術(shù)加速計(jì)算。

    1.2結(jié)構(gòu)重疊網(wǎng)格技術(shù)

    網(wǎng)格生成是數(shù)值模擬的重要部分之一,好的網(wǎng)格質(zhì)量是獲得合理計(jì)算結(jié)果的前提條件。結(jié)構(gòu)重疊網(wǎng)格邏輯關(guān)系簡單,魯棒性好,適合開展工程計(jì)算。在大部分計(jì)算區(qū)域,結(jié)構(gòu)重疊網(wǎng)格計(jì)算方法同結(jié)構(gòu)網(wǎng)格完全一樣,并能保證和結(jié)構(gòu)網(wǎng)格相同的計(jì)算精度、計(jì)算效率及流動(dòng)特性捕捉能力,而網(wǎng)格間的重疊技術(shù)降低了復(fù)雜拓?fù)渫庑蔚木W(wǎng)格生成難度[10-12]。另一方面,虛擬模擬實(shí)際問題時(shí),被模擬物體子部件間往往會(huì)出現(xiàn)相對運(yùn)動(dòng),如舵面/襟翼等部件受控偏轉(zhuǎn)、飛行器脫離掛架等,結(jié)構(gòu)重疊網(wǎng)格將各子部件劃為獨(dú)立的子區(qū)域,可有效地處理這些多體運(yùn)動(dòng)、分離問題[12]。求解多體、部件運(yùn)動(dòng)等位置或形狀隨時(shí)間變化的問題時(shí),網(wǎng)格無需重新生成,各子區(qū)域網(wǎng)格隨部件作剛性運(yùn)動(dòng),每一時(shí)間步只需要進(jìn)行一次重疊過程以建立插值信息關(guān)系,相比網(wǎng)格重構(gòu)、變形等動(dòng)網(wǎng)格方法實(shí)現(xiàn)更為簡單高效,且適用于任意幅度的運(yùn)動(dòng)[12-13]。

    1.3飛行器剛體6DOF運(yùn)動(dòng)方程組

    飛行器或運(yùn)動(dòng)部件作為剛體處理,在空間中剛體運(yùn)動(dòng)有6個(gè)自由度(DOF),剛體在空間的平動(dòng)、轉(zhuǎn)動(dòng)、位置及姿態(tài)的變化可以用6個(gè)動(dòng)力學(xué)方程和6個(gè)運(yùn)動(dòng)學(xué)方程充分描述[14]。

    剛體質(zhì)心動(dòng)力學(xué)方程為

    (2)

    式中:a和V分別為慣性系下加速度和速度矢量;P和F分別為推力和氣動(dòng)力在體軸坐標(biāo)系下矢量,下標(biāo)b表示體軸系;G為慣性系下重力矢量;L為坐標(biāo)轉(zhuǎn)換矩陣;m為質(zhì)量。

    慣性坐標(biāo)系下質(zhì)心運(yùn)動(dòng)學(xué)方程為

    (3)

    體軸坐標(biāo)系下剛體繞質(zhì)心轉(zhuǎn)動(dòng)的動(dòng)力學(xué)方程為

    (4)

    (5)

    由式(4)可得體軸系下繞質(zhì)心轉(zhuǎn)動(dòng)的動(dòng)力學(xué)方程為

    (6)

    再結(jié)合體軸系下角速度分量同姿態(tài)角變化關(guān)系式(7),可獲得關(guān)于俯仰角θ、偏航角φ及滾轉(zhuǎn)角γ的剛體動(dòng)力學(xué)和運(yùn)動(dòng)學(xué)聯(lián)立的方程組為

    (7)

    1.4控制系統(tǒng)?;?/p>

    飛行控制系統(tǒng)一般由舵回路、穩(wěn)定回路和制導(dǎo)回路這3個(gè)反饋回路構(gòu)成。這3個(gè)回路主要實(shí)現(xiàn)氣動(dòng)力(氣動(dòng)面)控制、擔(dān)負(fù)穩(wěn)定與控制功能,以及按照一定導(dǎo)引律給出控制指令[15]。

    在數(shù)值虛擬飛行中可以使用任意的控制器,復(fù)雜的控制器的實(shí)現(xiàn)、考察不在本文研究范圍中。為簡化起見,本文采用PID控制器,實(shí)現(xiàn)舵偏反饋控制,并形成穩(wěn)定的制導(dǎo)回路來跟蹤指定迎角變化規(guī)律。PID控制器如圖1所示。

    圖1 PID控制器框圖Fig.1 Block diagram of PID controller

    圖中:R(t)為設(shè)定值,C(t)為檢測值,E(t)為設(shè)定值和檢測值間的偏差,可見PID控制器是通過對偏差進(jìn)行運(yùn)算獲得控制值U(t)。PID運(yùn)算包含3部分,第1部分為比例控制Proportional,比例控制器為一個(gè)增益(放大)器,其將輸入信號放大Kp倍,Kp稱為增益(放大)系數(shù),增加Kp可提高系統(tǒng)準(zhǔn)確性,降低系統(tǒng)慣性,但會(huì)影響控制穩(wěn)定性。第2部分為積分控制器Integral,增加積分控制器參數(shù)Ki可以減小控制系統(tǒng)的穩(wěn)態(tài)誤差,但有可能引起系統(tǒng)振蕩。第3部分為比例微分控制器Derivative,改變Kd可影響系統(tǒng)響應(yīng)速度,即影響控制系統(tǒng)的快速性[16]。

    在模擬中,考慮到非定常CFD計(jì)算的物理時(shí)間步長遠(yuǎn)遠(yuǎn)小于舵機(jī)的響應(yīng)時(shí)間,因此假設(shè)舵面每隔Td接受一次控制信號。在Td時(shí)間內(nèi),舵偏規(guī)律設(shè)為

    (8)

    式中:Δd為需要接受控制指令后舵面需要偏轉(zhuǎn)的總舵偏角;d(t)為舵偏角關(guān)于時(shí)間的函數(shù);int為取整函數(shù)。

    1.5氣動(dòng)/運(yùn)動(dòng)/控制集成流程

    數(shù)值虛擬飛行的核心流程是氣動(dòng)/運(yùn)動(dòng)/控制(CFD/RBD/FCS)耦合求解。其中,CFD模塊采用非定常方法求解獲得每一物理時(shí)間步飛行器受到的氣動(dòng)力和力矩。剛體運(yùn)動(dòng)(RigidBodyDynamics,RBD)模塊根據(jù)力和力矩情況,求解剛體6DOF方程,更新物體位置(網(wǎng)格坐標(biāo))。

    本文使用松耦合算法[17]求解剛體動(dòng)力學(xué)方程與流動(dòng)控制方程的耦合問題,采用龍格庫塔法求解動(dòng)力學(xué)方程,采用雙歐法求解飛行器姿態(tài)變化。在獲得剛體下一時(shí)刻的位移和姿態(tài)角后,由式(9)獲得n+1時(shí)刻的網(wǎng)格點(diǎn)坐標(biāo)為

    Rn+1=Ln+1·R0

    (9)

    式中:Ln+1為由n+1時(shí)刻的姿態(tài)角確定的轉(zhuǎn)換矩陣;R0為初始的網(wǎng)格坐標(biāo)。

    控制系統(tǒng)(FlightControlSystem,F(xiàn)CS)模塊根據(jù)飛行器姿態(tài)與期望的偏差形成反饋控制,結(jié)合重疊網(wǎng)格等動(dòng)網(wǎng)格技術(shù)控制氣動(dòng)操縱面,控制飛行器姿態(tài)。然后再反復(fù)上述過程,直到運(yùn)動(dòng)結(jié)束。流程如下:

    1) 為減小計(jì)算時(shí)間,采用定常計(jì)算獲得初始流場。

    2) 以初場為基礎(chǔ)開展非定常流動(dòng)計(jì)算。

    3) 每推進(jìn)一個(gè)物理時(shí)間步,求解飛行器的氣動(dòng)力和力矩,代入RBD求解獲得物體位移,更新整體網(wǎng)格位置。

    4) 根據(jù)FCS,控制操縱面運(yùn)動(dòng),更新操縱面物體(網(wǎng)格)位置。

    5) 開展下一物理時(shí)間步流場計(jì)算,重復(fù)步驟2~5,直至物理時(shí)間推進(jìn)結(jié)束。

    2 數(shù)值方法驗(yàn)證

    用三維的超聲速彈道外形[18](HyperBallisticShape,HBS)對程序和方法的正確性和可信度進(jìn)行驗(yàn)證。HBS外形如圖2所示,具有試驗(yàn)數(shù)據(jù)可供對比,是一個(gè)廣泛用于超聲速動(dòng)態(tài)非定常計(jì)算程序驗(yàn)證的標(biāo)準(zhǔn)算例。

    圖2 超聲速彈道模型外形及尺寸Fig.2 Shape and size of HBS model

    采用文獻(xiàn)[18]中來流參數(shù)開展計(jì)算,具體為:來流馬赫數(shù)6.85,質(zhì)心在彈長0.72倍處,以前段直徑d為特征長度的雷諾數(shù)為0.72×106。對比在0°、5°、10°、15°和20° 迎角下計(jì)算和試驗(yàn)獲得的俯仰導(dǎo)數(shù),俯仰動(dòng)導(dǎo)數(shù)采用振幅為1° 的正弦強(qiáng)迫運(yùn)動(dòng)計(jì)算獲得,強(qiáng)迫運(yùn)動(dòng)的減縮頻率取0.01,計(jì)算物理時(shí)間步長為2ms。圖3(a)為計(jì)算獲得的靜導(dǎo)數(shù)Cmα同試驗(yàn)值對比,程序準(zhǔn)確計(jì)算出了約在8° 迎角附近導(dǎo)數(shù)符號的改變。圖3(b)為俯仰動(dòng)導(dǎo)數(shù)同試驗(yàn)值對比,計(jì)算同試驗(yàn)數(shù)據(jù)符合較好,表明本文程序和非定常方法是可行的。

    圖3 HBS俯仰靜、動(dòng)導(dǎo)數(shù)計(jì)算與試驗(yàn)數(shù)據(jù)對比Fig.3 Comparison of computational static and dynamic pitch derivative with test data of HBS

    3 數(shù)值模擬結(jié)果及分析

    3.1返回艙外形平面運(yùn)動(dòng)動(dòng)態(tài)特性研究

    起源號返回艙網(wǎng)格如圖4所示,具體外形尺寸見文獻(xiàn)[19],模型物理參數(shù)見表1。

    圖4 起源號返回艙表面及對稱面網(wǎng)格Fig.4 Surface and symmetry plane grid of Genesis capsule

    表1 起源號返回艙模型物理屬性

    Table 1 Physical property of Genesis capsule model

    ItemValued/mm69.42M/g633Izz/kg·m23.01×10-4Iyy/kg·m21.99×10-4L/mm44.33Xcg/mm23.83

    圖5為初始馬赫數(shù)Ma0為4.5,釋放迎角為0.5° 時(shí)返回艙總迎角隨射程變化規(guī)律同文獻(xiàn)數(shù)據(jù)對比。由圖可見,該狀態(tài)彈道的前段總迎角較小,在120 m后迎角開始發(fā)散,計(jì)算也出現(xiàn)了相同的趨勢。需要說明的是,文獻(xiàn)中并未說明返回艙釋放的角速度,計(jì)算時(shí)設(shè)初始角速度為零,從結(jié)果來看這不會(huì)對后段彈道迎角發(fā)散這一規(guī)律產(chǎn)生較大影響。圖6為Ma=3.5,大氣參數(shù)為海平面時(shí),采用CFD自由振蕩計(jì)算獲得的動(dòng)導(dǎo)數(shù)同文獻(xiàn)[19]中數(shù)據(jù)對比。計(jì)算結(jié)果同文獻(xiàn)數(shù)據(jù)接近,符合該返回艙在馬赫數(shù)大于2.5時(shí)在大迎角下為動(dòng)穩(wěn)定的結(jié)論。

    圖5 Ma0=4.5,釋放迎角為0.5° 時(shí),返回艙總迎角隨彈道變化規(guī)律計(jì)算與試驗(yàn)對比Fig.5 Comparison of computational trajectory history of total angle of attack with test data at Ma0=4.5 and release α=0.5°

    圖6 Ma=3.5時(shí),返回艙俯仰動(dòng)導(dǎo)數(shù)計(jì)算同文獻(xiàn)對比Fig.6 Comparison of computational capsule dynamic pitch derivative with reference data at Ma=3.5

    以上計(jì)算結(jié)果同試驗(yàn)結(jié)果吻合良好,可確定本文方法的正確性。下面的研究主要關(guān)注不同自由度運(yùn)動(dòng)對返回艙運(yùn)動(dòng)特性的影響,因此本外形算例來流參數(shù)統(tǒng)一為:來流馬赫數(shù)為3.5,飛行高度為25 km,來流壓力、溫度根據(jù)地球大氣參數(shù)表獲得,計(jì)算皆采用層流假設(shè),物理時(shí)間步長統(tǒng)一取為0.3 ms。由于自由俯仰運(yùn)動(dòng)時(shí),該系統(tǒng)是自治系統(tǒng),其極限環(huán)狀態(tài)(振幅和頻率)與初值無關(guān),因此,為了較快出現(xiàn)極限環(huán)運(yùn)動(dòng),計(jì)算預(yù)置10° 迎角,并且設(shè)置返回艙的質(zhì)量為表1中的1/4,以減小慣量,使運(yùn)動(dòng)迅速出現(xiàn)振蕩,以節(jié)約計(jì)算時(shí)間。

    3.1.1平面單自由度俯仰運(yùn)動(dòng)

    當(dāng)返回艙繞質(zhì)心做單自由度俯仰振動(dòng)時(shí),其動(dòng)力學(xué)方程為

    (10)

    在某迎角的小的鄰近范圍內(nèi),可以假設(shè)俯仰靜、動(dòng)導(dǎo)數(shù)為常量,則式(10)為一典型的二階阻尼系統(tǒng),其解析解為

    (11)

    式中:ξl為阻尼;ωl為振動(dòng)頻率,具體表達(dá)式為

    (12)

    圖7為返回艙的靜、動(dòng)導(dǎo)數(shù)隨迎角變化情況。由圖可見,各迎角下該外形靜導(dǎo)數(shù)變化不大,而動(dòng)導(dǎo)數(shù)在小迎角處是正值,表現(xiàn)為動(dòng)不穩(wěn)定;迎角大于5° 后,動(dòng)導(dǎo)數(shù)變?yōu)樨?fù)值,表現(xiàn)為動(dòng)穩(wěn)定。因此,返回艙在小迎角下受到擾動(dòng),由于動(dòng)不穩(wěn)定性,擾動(dòng)振幅將增大,在大迎角時(shí),氣流起阻尼作用,從而衰減振幅。在一定條件下,這樣的動(dòng)導(dǎo)數(shù)變化規(guī)律可使得飛行器出現(xiàn)俯仰極限環(huán)振蕩狀態(tài)。

    圖7 “起源號”返回艙各迎角下俯仰靜、動(dòng)導(dǎo)數(shù)Fig.7 Static and dynamic pitch derivative of Genesiscapsule under different angles of attack

    圖8為自由俯仰振蕩迎角隨時(shí)間的變化規(guī)律??梢?,由前述動(dòng)導(dǎo)數(shù)特性,從10° 迎角釋放外形后,飛行器振幅逐漸增加,當(dāng)振幅達(dá)到15° 左右時(shí),飛行器出現(xiàn)等周期的極限環(huán)振蕩。

    圖8 自由俯仰運(yùn)動(dòng)迎角隨時(shí)間變化規(guī)律Fig.8 Time history of angle of attack of free pitching motion

    3.1.2平面俯仰+沉浮雙自由度運(yùn)動(dòng)

    下面研究引入質(zhì)心沉浮運(yùn)動(dòng)后,返回艙的運(yùn)動(dòng)特性。圖 9為雙自由度和單自由度運(yùn)動(dòng)下飛行器迎角隨時(shí)間變化規(guī)律。由圖可見,引入沉浮運(yùn)動(dòng)后,返回艙仍形成等周期和振幅極限環(huán)俯仰振動(dòng),但相比單自由度振動(dòng),其振幅加大,頻率減小。

    圖9 雙自由度和單自由度運(yùn)動(dòng)下迎角隨時(shí)間變化規(guī)律的對比Fig.9 Comparison of time histories of angle of attack of one and two DOF motion

    圖10為返回艙質(zhì)心沉浮速度V隨時(shí)間變化規(guī)律,其與迎角隨時(shí)間變化有1/4周期的相位差。此外,注意到,由沉浮運(yùn)動(dòng)引起的迎角變化不會(huì)超過0.02°,這表明引起俯仰振動(dòng)頻率和振幅變化的主要原因并不是沉浮速度,而是引入沉浮運(yùn)動(dòng)后,系統(tǒng)阻尼變小導(dǎo)致的,這跟文獻(xiàn)中理論分析一致[20-21]。返回艙平面雙自由度運(yùn)動(dòng)下的俯仰靜、動(dòng)導(dǎo)數(shù)及失穩(wěn)特性和單自由度運(yùn)動(dòng)的類似,在此不再贅述。

    圖10 返回艙雙自由度運(yùn)動(dòng)質(zhì)心沉浮速度隨時(shí)間的變化規(guī)律Fig.10 Time history of centroid velocity of two DOF motion of return capsule

    3.1.3平面自由飛三自由度運(yùn)動(dòng)

    返回艙返回時(shí),在氣動(dòng)阻力的作用下減速運(yùn)動(dòng),引入流向運(yùn)動(dòng)自由度,即可考察返回艙平面自由飛行的動(dòng)態(tài)特性,這種自由飛更加接近真實(shí)飛行情況。圖11為平面自由飛運(yùn)動(dòng)下飛行器迎角隨時(shí)間的變化規(guī)律。由圖可見,考慮氣動(dòng)阻力后,返回艙振幅更大,這表明飛行器變得更為動(dòng)不穩(wěn)定。此外,可以看出,振蕩頻率逐漸減小,振蕩周期變長,這是由于氣動(dòng)阻力使飛行器速度降低造成的。

    圖11 平面自由飛和雙自由度下迎角隨時(shí)間的變化規(guī)律對比Fig.11 Comparison of time histories of angle of attack of planar free and two DOF motion

    圖12為平面自由飛運(yùn)動(dòng)下返回艙質(zhì)心沉浮速度隨時(shí)間的變化規(guī)律??梢姡杂娠w下沉浮速度振幅比雙自由度運(yùn)動(dòng)下的振幅有所增加,且周期逐漸變長。圖13為質(zhì)心流向速度隨時(shí)間的變化規(guī)律,空氣阻力的存在使飛行器流向速度由3.5馬赫降到了2.2馬赫左右。同樣,基于自由飛計(jì)算開展參數(shù)辨識,也可獲得相應(yīng)的靜、動(dòng)導(dǎo)數(shù),在此不再累述。

    圖12 平面自由飛質(zhì)心沉浮速度隨時(shí)間變化規(guī)律Fig.12 Time history of centroid velocity of planar free motion

    圖13 平面自由飛質(zhì)心流向速度隨時(shí)間變化規(guī)律Fig.13 Time history of flow direction centroid velocity of planar free motion

    3.2受控導(dǎo)彈縱向數(shù)值虛擬飛行模擬

    計(jì)算外形為尾部安裝有4個(gè)翼的導(dǎo)彈。該外形具有大量的靜、動(dòng)態(tài)試驗(yàn)數(shù)據(jù)及計(jì)算對比數(shù)據(jù)[22-23]。該導(dǎo)彈外形涵蓋了大部分戰(zhàn)術(shù)導(dǎo)彈的典型氣動(dòng)布局特征。導(dǎo)彈為十字翼布局,尾舵為均勻厚度。計(jì)算外形尺寸及舵號定義見圖14(圖中CG為模型重心),模型物理參數(shù)見表 2。圖15為Ma=4.42,Red=2.96×107時(shí)采用不同時(shí)間步長計(jì)算獲得的俯仰角隨射程變化規(guī)律同試驗(yàn)值之間的對比,其中Red是以彈身直徑d為特征長度的雷諾數(shù);圖 16為Red=8.6×106,Ma=1.96,無量綱質(zhì)心位置X/d=6.1時(shí)導(dǎo)彈各迎角下動(dòng)導(dǎo)數(shù)計(jì)算與試驗(yàn)對比。以上計(jì)算結(jié)果同試驗(yàn)結(jié)果吻合良好,可確定本文方法的正確性。此外由圖 15可見,采用0.5 ms和1 ms時(shí)間步長的計(jì)算結(jié)果同試驗(yàn)值得差異皆不大,因此為減少計(jì)算機(jī)時(shí),之后計(jì)算中物理時(shí)間步長均取1 ms。本文只考慮俯仰通道模擬,因此只考慮1號和2號舵實(shí)現(xiàn)舵偏功能,舵后緣向下偏轉(zhuǎn)定義為正舵偏,采用重疊網(wǎng)格實(shí)現(xiàn)舵偏功能。此外,在下述計(jì)算皆采用一方程的Spalart-Allmaras湍流模型模擬湍流黏性,子迭代收斂殘差取為0.01。

    圖14 基本帶翼導(dǎo)彈模型外形及尺寸Fig.14 Shape and size of basic finner

    表2 基本帶翼導(dǎo)彈計(jì)算模型物理屬性

    Table 2Physical property of calculation model of basic finner

    ItemValueL/m0.3M/kg1.58Ixx/kg·m21.92×10-4Iyy/kg·m297.85×10-4d/m0.03Xcg/m0.165

    圖15 導(dǎo)彈俯仰角隨射程變化計(jì)算同試驗(yàn)值對比Fig.15 Comparison of computational projectile pitch angle variation with range and test data

    圖16 Red=8.6×106時(shí)動(dòng)導(dǎo)數(shù)計(jì)算試驗(yàn)值對比Fig.16 Comparison of computational dynamic derivative with test data at Red=8.6×106

    3.2.1飛行器姿態(tài)穩(wěn)定特性研究

    飛行器的穩(wěn)定特性包括動(dòng)穩(wěn)定性和動(dòng)操縱性。當(dāng)飛行器受外界擾動(dòng)后的運(yùn)動(dòng)體現(xiàn)其動(dòng)穩(wěn)定性,當(dāng)飛行器受操縱力控制后的運(yùn)動(dòng)體現(xiàn)其動(dòng)操縱性。設(shè)計(jì)擾動(dòng)規(guī)律,開展飛行器姿態(tài)穩(wěn)定研究,可以檢驗(yàn)飛行器不同的運(yùn)動(dòng)模態(tài),測試飛行器在各種極端情況下的靜、動(dòng)穩(wěn)定性,且可為氣動(dòng)參數(shù)辨識提供足夠的信息。

    采用時(shí)域法分析一個(gè)系統(tǒng)的性能,一般采用一些典型的輸入作為測試信號,如階躍輸入,脈沖輸入等[16],結(jié)合導(dǎo)彈飛行的特點(diǎn),本文選取如下4種舵偏規(guī)律作為擾動(dòng)輸入,測試基本外形導(dǎo)彈在不同舵偏擾動(dòng)下的姿態(tài)穩(wěn)定性:

    1) 階躍型舵偏輸入

    d(t)=5°t≥0

    (13)

    2) 方形脈沖舵偏輸入

    (14)

    3) 梯形舵偏輸入

    (15)

    4) 正弦舵偏輸入

    (16)

    式中:d(t)為舵偏角隨時(shí)間的函數(shù);t為物理時(shí)間。下面計(jì)算中來流條件為Ma=2,高度為10 km,來流壓力及密度根據(jù)大氣參數(shù)表獲得,飛行器初始迎角為0°,無姿態(tài)角速度及姿態(tài)角加速度。

    基于小擾動(dòng)理論,基本外形導(dǎo)彈帶舵面控制的無量綱縱向運(yùn)動(dòng)方程可寫為

    (17)

    式中:Iz為轉(zhuǎn)動(dòng)慣量;θ為俯仰角。

    式(17)的特征方程為

    (18)

    圖17為階躍舵偏下導(dǎo)彈俯仰角響應(yīng)曲線及俯仰角速度響應(yīng)曲線。階躍舵偏反映了舵偏(擾動(dòng))突變的一種極限情況,而飛行受擾動(dòng)后收斂于平衡位置的過程反映了飛行的動(dòng)穩(wěn)定性。由圖可見,階躍舵偏(擾動(dòng))引起的角速度達(dá)200°/s左右,飛行器俯仰角振蕩收斂到-5.5° 左右,振蕩次數(shù)、幅度衰減情況反映了飛行器的操縱品質(zhì)。式(17)在階躍舵偏下的解為

    (19)

    式中:αeq為配平攻角;φ為相位裕度,具體表達(dá)式為

    (20)

    結(jié)合式(18)與式(19),即可辨識出俯仰阻尼導(dǎo)數(shù)等參數(shù)。

    圖18為偶極脈沖舵偏下導(dǎo)彈俯仰角響應(yīng)曲線及俯仰角速度響應(yīng)曲線。在t=0.2 s時(shí),此時(shí)導(dǎo)彈具有正俯仰角速度,0.2 s時(shí)刻后,舵偏變負(fù),產(chǎn)生的抬頭力矩加劇了抬頭的趨勢。當(dāng)干擾頻率與飛行器固有頻率接近時(shí),可能引起飛行器共振。偶極脈沖舵偏可反映飛行器受到瞬時(shí)、隨機(jī)且頻率接近飛行器固有頻率的干擾的動(dòng)態(tài)性能。

    圖19為梯形脈沖舵偏下導(dǎo)彈俯仰角響應(yīng)曲線及俯仰角速度響應(yīng)曲線。梯形脈沖舵偏反映了等速舵偏與階躍舵偏的組合情況。由圖可見,此時(shí)超調(diào)量小于階躍型舵偏,且振動(dòng)頻率也小于階躍型舵偏,最大俯仰角速度約為100°/s,這種舵偏規(guī)律代表了緩慢的機(jī)動(dòng)飛行。

    圖17 階躍舵偏下的俯仰響應(yīng)曲線Fig.17 Pitch response curve of phase step impulse rudder deflexion

    圖18 偶極脈沖舵偏下的俯仰響應(yīng)曲線Fig.18 Pitch response curve of dipole impulse rudder deflexion

    圖19 梯形脈沖舵偏下的俯仰響應(yīng)曲線Fig.19 Pitch response curve of echelon impulse rudder deflexion

    圖20為諧波脈沖舵偏下導(dǎo)彈俯仰角響應(yīng)曲線及俯仰角速度響應(yīng)曲線。由圖可見,此種情況下,舵偏(擾動(dòng))引起的最大俯仰角速度及俯仰力矩相對前3種情況是最小的,此時(shí)飛行器實(shí)施精確跟蹤,處于理想化的反復(fù)修正過程。

    圖20 諧波脈沖舵偏下的俯仰響應(yīng)曲線Fig.20 Pitch response curve of harmonic wave impulse rudder deflexion

    3.2.2基于數(shù)值虛擬飛行對控制參數(shù)的整定

    飛行器在飛行過程中,控制系統(tǒng)會(huì)根據(jù)飛行狀態(tài)發(fā)出實(shí)時(shí)控制指令,修正飛行姿態(tài)。系統(tǒng)受控后,要么導(dǎo)致受控發(fā)散,要么趨于受控狀態(tài)。飛行器是否會(huì)出現(xiàn)受控發(fā)散,即控制是否具有穩(wěn)定性,可以采用3.3.1節(jié)中的方法檢驗(yàn)。而對于受控穩(wěn)定情況,則要考慮控制系統(tǒng)的準(zhǔn)確性,快速性及平穩(wěn)性。

    控制系統(tǒng)的準(zhǔn)確性,快速性及平穩(wěn)性與控制系統(tǒng)的參數(shù)有很強(qiáng)的關(guān)系。同一種控制方式,面對各種各樣的飛行器,滿足上述要求的控制參數(shù)也大不相同。即使對于同一個(gè)飛行器,由于飛行過程中燃料消耗導(dǎo)致的質(zhì)量、慣性矩及重心的物理參數(shù)的改變,也會(huì)使飛行器受控后響應(yīng)發(fā)生改變。而對于處于非定常/非線性流動(dòng)下的飛行器,氣動(dòng)力、力矩的突變、遲滯特性更會(huì)使得基于線化氣動(dòng)導(dǎo)數(shù)等參數(shù)獲得控制參數(shù)不能在全范圍內(nèi)保證控制的穩(wěn)定、準(zhǔn)確、快速及平穩(wěn)性。因此,整定控制參數(shù),是數(shù)值虛擬飛行一個(gè)重要的應(yīng)用,也是風(fēng)洞等常規(guī)手段難以實(shí)現(xiàn)的。

    本文采用PID控制器,實(shí)現(xiàn)舵偏反饋控制,跟蹤指定迎角規(guī)律。測試了不同PID參數(shù)(如表3 所示)下導(dǎo)彈的縱向受控性能,以獲得較好的控制參數(shù)。

    模擬的來流條件為:Ma=2,高度為10 km,初始迎角為a=0°,初始俯仰角速度為零,角加速度也為零。導(dǎo)彈姿態(tài)角變化要求為:機(jī)動(dòng)到5° 迎角飛行,超調(diào)不超過20%,并維持姿態(tài)穩(wěn)定。按上述要求及計(jì)算模型物理屬性,根據(jù)PID控制理論[16],選取如下范圍的PID控制參數(shù)及舵受控時(shí)間Td,對該范圍的參數(shù)整定,以達(dá)到最佳控制效果。

    表3 PID控制參數(shù)

    圖21~圖24為各控制參數(shù)下導(dǎo)彈的響應(yīng)時(shí)間歷程。從圖中可以看到,控制系統(tǒng)根據(jù)飛行器姿態(tài)角的誤差發(fā)出控制指令,控制舵面,飛行器實(shí)施機(jī)動(dòng),然后反復(fù)調(diào)節(jié),直至飛行器達(dá)到期望姿態(tài)。由相圖可見,初始階段,飛行器的俯仰角速度增加,迎角也隨之增加,當(dāng)迎角接近指定值5° 時(shí),微分控制器開始制動(dòng),減小俯仰角速度,當(dāng)姿態(tài)角度超過5° 后,在控制器作用下,舵面反復(fù)調(diào)整,使得俯仰角速度變負(fù),最后逐漸收斂于相圖上的螺旋點(diǎn)。

    對比Case 1和Case 2可以看出,Case 1超調(diào)量約20%,Case 2超調(diào)量約10%,Case 1的最大俯仰角速率約為33°/s,而Case 2為30°/s。這表明增加Kp,會(huì)提高系統(tǒng)的響應(yīng)速度,但會(huì)使得超調(diào)增加,振蕩過渡過程變長。對比Case 3和Case 4相圖可以看出,在迎角達(dá)到2° 左右時(shí),控制系統(tǒng)開始制動(dòng),即減小俯仰角速率。制動(dòng)后Case 3的最大俯仰角速率為39°/s,Case 4的最大俯仰角速率為42°/s。這表明增加Kd,使得系統(tǒng)的制動(dòng)提前,降低俯仰角速度峰值,能減少過渡振蕩過程。對比Case 2和Case 3,Td減小加快了俯仰角速率,提高了系統(tǒng)響應(yīng)速度。

    圖21 Case 1的俯仰響應(yīng)曲線Fig.21 Pitch response curve of Case 1

    圖22 Case 2的俯仰響應(yīng)曲線Fig.22 Pitch response curve of Case 2

    圖23 Case 3的俯仰響應(yīng)曲線Fig.23 Pitch response curve of Case 3

    上述算例表明,Case 3的控制參數(shù)有較好的特性。下面采用Case 3控制參數(shù)對有初始俯仰角速度的狀態(tài)進(jìn)行控制,驗(yàn)證控制效果,初始俯仰角速度設(shè)為50°/s。

    圖25為采用Case 3的控制參數(shù)下的俯仰角控制曲線和相圖。導(dǎo)彈在初始俯仰角速度的作用下,迎角增加。由于導(dǎo)彈是靜、動(dòng)穩(wěn)定,且此時(shí)舵偏角較小,俯仰角速度迅速減小變?yōu)樨?fù)值。當(dāng)舵偏角負(fù)向變大時(shí),產(chǎn)生抬頭力矩,迎角增加,然后逐漸調(diào)整到預(yù)設(shè)迎角。初始俯仰角速度的存在使得彈體動(dòng)阻尼效應(yīng)較為明顯,且初始角速度越大,動(dòng)態(tài)效應(yīng)表現(xiàn)越強(qiáng)。

    圖24 Case 4的俯仰響應(yīng)曲線Fig.24 Pitch response curve of Case 4

    圖25 有初始俯仰角速度時(shí)的俯仰響應(yīng)曲線Fig.25 Pitch response curve of state with initial pitching angular speed

    4 結(jié) 論

    本文利用數(shù)值虛擬飛行技術(shù),對返回艙外形及受控基本帶翼導(dǎo)彈外形開展了動(dòng)態(tài)特性研究。

    對返回艙外形的研究表明:返回艙外形具有靜態(tài)穩(wěn)定性和動(dòng)不穩(wěn)定性。對于平面三自由度俯仰運(yùn)動(dòng),俯仰阻尼導(dǎo)數(shù)決定了返回艙的動(dòng)穩(wěn)定性。引入沉浮運(yùn)動(dòng)使得返回艙動(dòng)穩(wěn)定性變差。引入流向自由度后,氣動(dòng)阻尼使飛行器速度降低,動(dòng)穩(wěn)定性變差,且極限環(huán)周期逐漸變長。

    對受控基本外形導(dǎo)彈的研究表明:通過輸入典型的擾動(dòng)(舵偏)信號,可以檢測導(dǎo)彈系統(tǒng)的穩(wěn)定、受控性能。不同控制參數(shù)對系統(tǒng)的控制性能有較大影響。增大增益系數(shù),可以提高響應(yīng)速度,但會(huì)增加超調(diào)量,引起振蕩。增加微分系數(shù),可以控制系統(tǒng)制動(dòng)幅度及時(shí)間,減小超調(diào)量,使振蕩過渡過程盡快結(jié)束。

    上述研究表明,結(jié)合結(jié)構(gòu)重疊網(wǎng)格技術(shù)、閉環(huán)控制律的數(shù)值虛擬飛行技術(shù)能夠有效地研究常規(guī)手段無法研究的飛行器動(dòng)態(tài)特性問題。其應(yīng)用靈活性好,周期短,對探索飛行器氣動(dòng)/運(yùn)動(dòng)耦合問題,控制參數(shù)的整定,控制律的研究和檢測,獲取飛行器非定常的運(yùn)動(dòng)導(dǎo)數(shù)以及飛行的安全性測試具有較好的工程應(yīng)用價(jià)值及應(yīng)用前景。

    [1]呂光男. 風(fēng)洞虛擬飛行實(shí)驗(yàn)中的飛行力學(xué)與控制問題研究[D]. 南京: 南京航空航天大學(xué), 2009: 1-13

    LYU G N. The problem of flight mechanics and control in wind tunnel based virtual flight testing[D]. Nanjing: Nanjing University of Aeronautics and Astonautics, 2009:1-13 (in Chinese).

    [2]李周復(fù). 風(fēng)洞特種實(shí)驗(yàn)技術(shù)[M]. 北京: 航空工業(yè)出版社, 2010: 210-250.

    LI Z F. Special wind tunnel testing technology[M]. Beijing: Aviation Industry Press, 2010: 210-250 (in Chinese).

    [3]劉偉, 楊小亮, 趙云飛. 高超聲速飛行器加速度導(dǎo)數(shù)數(shù)值模擬[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2010, 28(4): 426-429.

    LIU W, YANG X L, ZHAO Y F. Numerical simulation of acceleration derivative of hypersonic aircraft[J]. Acta Aerodynamica Sinica, 2010, 28(4): 426-429 (in Chinese).

    [4]RATLIFF C L, MARQUART E J. An assessment of a potential test technique: virtual flight testing (VFT): AIAA-1995-3415[R]. Reston: AIAA, 1995.

    [5]GEBERT G, KELLY J, LOPEZ J. Wind tunnel based virtual flight testing: AIAA-2000-0829[R], Reston: AIAA, 2000.

    [6]MAGILL J C, CATALDI P, MORENCY J R, et al. Demonstration of a wire suspension for wind-tunnel virtual flight testing[J]. Journal of Spacecraft and Rockets, 2009, 46(3): 624-633.

    [7]張來平, 馬戎, 常興華, 等. 虛擬飛行中氣動(dòng)運(yùn)動(dòng)和控制耦合的數(shù)值模擬技術(shù)[J]. 力學(xué)進(jìn)展, 2014, 44: 201410.

    ZHANG L P, MA R, CHANG X H, et al. Review of aerodynamics/kinematics/flight-control coupling methods in virtual flight simulations[J]. Advances in Mechanics, 2014, 44: 201410.

    [8]閻超. 計(jì)算流體力學(xué)方法及應(yīng)用[M]. 北京: 北京航空航天大學(xué)出版社, 2006: 18-24.

    YAN C. Computational fluid dynamics methods and applications. Beijing[M]. Beijing: Beihang University Press 2006: 18-24 (in Chinese).

    [9]SALAS M D. Digital flight: The last CFD aeronautical grand challenge[J]. Journal of Scientific Computing, 2006, 28(2-3): 479-505.

    [10]ROGERS S E,ROTH K,NASH S M,et al.Advances in overset CFD processes applied to subsonic high-lift aircraft: AIAA-2000-4216[R]. Reston: AIAA, 2000.

    [11]GOETZ H K, JAMES E K, HENRY C L, et al. Validation of overflow for computing plume effects during the ARES I stage separation process: AIAA-2011-0170[R]. Reston: AIAA, 2011

    [12]范晶晶, 閻超, 張輝. 重疊網(wǎng)格洞面優(yōu)化技術(shù)的改進(jìn)與應(yīng)用[J]. 航空學(xué)報(bào), 2010, 31(6): 1127-1133.

    FAN J J, YAN C, ZHANG H. Improvement of hole-surface optimization technique in overset grids and its application[J]. Acta Aeronauticaet Astronautica sinica, 2010, 31(6): 1127-1133 (in Chinese).

    [13]袁武. 新型重疊網(wǎng)格方法研究及在復(fù)雜多體氣動(dòng)問題中的應(yīng)用[D]. 北京: 北京航空航天大學(xué), 2012: 24-35.

    YUAN W. Investigations on novel chimera grid methods and its applications to complex multibody aerodynamic problems[D]. Beijing: Beihang University, 2012: 24-35 (in Chinese).

    [14]肖業(yè)倫, 航空航天器運(yùn)動(dòng)的建?!w行動(dòng)力學(xué)的理論基礎(chǔ)[M]. 北京: 北京航空航天大學(xué)出版社, 2003: 64-78.

    XIAO Y L. Modeling of the motion of space vehicle—the basic of flight mechnics[M]. Beijing: Beihang University Press 2003: 64-78 (in Chinese).

    [15]吳森堂, 費(fèi)玉華. 飛行控制系統(tǒng)[M]. 北京: 北京航空航天大學(xué)出版社, 2009: 85-93.

    WU S T, FEI Y H, Flight control system[M]. Beijing: Beihang University Press, 2009: 85-93 (in Chinese).

    [16]MEI X R, ZHUANG X Y. Principle of auto control[M]. Beijing: Science Press, 2009: 54-61.

    [17]栗長江. 高超聲速飛行器多自由度動(dòng)態(tài)特性研究[D]. 長沙: 國防科學(xué)技術(shù)大學(xué), 2010: 51-54.

    LI C J. Study of mutidimensional motion of hypersonic flight vehicle[D]. Changsha: National University of Science and Technology, 2010: 51-54 (in Chinese).

    [18]EAST R A, HUTT G R. Comparison of predictions and experimental data for hypersonic pitching motion stability[J]. Journal of Spacecraft, 1988, 25(3): 225-233

    [19]CHEATWOOD F, WINCHENBACH G, HATHAWAY W, et al. Dynamic stability testing of the genesis sample return capsule: AIAA-2000-1009[R]. Reston: AIAA, 2000.

    [20]袁先旭. 非定常流動(dòng)數(shù)值模擬及飛行器動(dòng)態(tài)特性分析研究[D]. 綿陽: 中國空氣動(dòng)力研究與發(fā)展中心研究生部, 2002: 101-117.

    YUAN X X. Numerical simulation for unsteady flowsand research on dynamic characteristicsof vehicle[D]. Mianyang: China Aerodynamics Research and Development Center. 2002: 101-117 (in Chinese).

    [21]席柯. 虛擬飛行數(shù)值模擬方法及飛行器動(dòng)態(tài)特性研究[D]. 北京: 北京航空航天大學(xué), 2015.

    XI K. Virtual flight numerical simulation method and research on dynamic characteristics of vehicle[D]. Beijing, Beihang University, 2015 (in Chinese).

    [22]DUPUIS A D, HATHAWAY W. Aeroballistic range tests of the basic finner reference projectile at supersonic velocities: DREV-TM-9703[R]. Valcartier: Defense Research Establishment, 1997.

    [23]BHAGWANDIN V A, SAHU J. Numerical prediction of pitch damping stability derivatives for finned projectiles: AIAA-2011-3028[R]. Reston: AIAA, 2011.

    黃宇男, 博士研究生。主要研究方向: 計(jì)算流體力學(xué)、 混合網(wǎng)格數(shù)值方法。

    Tel: 010-82338071

    E-mail: huangyu@buaa.edu.cn

    閻超男, 博士, 教授, 博士生導(dǎo)師。主要研究方向: 計(jì)算流體力學(xué)、 空氣動(dòng)力學(xué)。

    Tel: 010-82317019

    E-mail: yanchao@buaa.edu.cn

    Analysis of flying vehicle’s dynamic characteristics based onnumerical virtual flight technology

    HUANG Yu1, YAN Chao1,*, XI Ke2, WANG Wen1

    1. School of Aeronautic Science and Engineering, Beihang University, Beijing100083, China 2. Institute of Ordnance Industry Navigation and Control Technology, Beijing100089, China

    With the numerical virtual flight technology based on structured overlapping grid, closed PID controller, rudder control, rigid body motion and unsteady N-S equation solver, the unsteady motion characteristic, dynamic stability and adjustment of control parameter of Genesis capsule and basic finner projectile have been simulated. The motion characteristic with different degree of freedom (DOF) motion, the motion stability and the adjustment of control parameter with disturbance of those flying vehicles have been analyzed. The calculation results show that the study of flying vehicle’s nonlinear motion under the condition of complex fluid flows can be effectively carried out by numerical virtual flight technology, which has practical value and application prospect in the area of simulation and prediction of flying vehicle’s motion and controlling characteristic, fluid mechanics study and design of control rule under the condition of unsteady and nonlinear aerodynamics.

    CFD; overlapping grid; PID control; virtual flight; dynamic derivative; control parameters adjustment; dynamic characteristic; unsteady flow

    2016-01-11; Revised: 2016-02-17; Accepted: 2016-04-06; Published online: 2016-04-2015:26

    . Tel.: 010-82317019E-mail: yanchao@buaa.edu.cn

    2016-01-11; 退修日期: 2016-02-17; 錄用日期: 2016-04-06;

    時(shí)間: 2016-04-2015:26

    www.cnki.net/kcms/detail/11.1929.V.20160420.1526.004.html

    .Tel.: 010-82317019E-mail:yanchao@buaa.edu.cn

    10.7527/S1000-6893.2016.0120

    V211.3

    A

    1000-6893(2016)08-2525-14

    引用格式: 黃宇, 閻超, 席柯, 等. 基于數(shù)值虛擬飛行技術(shù)的飛行器動(dòng)態(tài)特性分析[J]. 航空學(xué)報(bào), 2016, 37(8): 2525-2538. HUANG Y, YAN C, XI K, et al. Analysis of flying vehicle’s dynamic characteristics based on numerical virtual flight technology[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(8): 2525-2538.

    http://hkxb.buaa.edu.cnhkxb@buaa.edu.cn

    URL: www.cnki.net/kcms/detail/11.1929.V.20160420.1526.004.html

    猜你喜歡
    飛行技術(shù)返回艙迎角
    “ 神舟十三號”返回艙安全著陸
    軍事文摘(2022年12期)2022-07-13 03:11:16
    連續(xù)變迎角試驗(yàn)數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    The innovation path of Chinese pilot training model
    飛行技術(shù)專業(yè)“機(jī)載設(shè)備”課程教學(xué)改革與研究
    中國市場(2016年36期)2016-10-19 04:56:47
    實(shí)踐十號返回艙回家
    太空探索(2016年5期)2016-07-12 15:17:55
    多用途飛船縮比返回艙成功著陸
    太空探索(2016年7期)2016-07-10 12:10:15
    民航特色專業(yè)學(xué)風(fēng)建設(shè)的實(shí)踐與探索
    失速保護(hù)系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    高超音速技術(shù)將用于未來客機(jī)中
    免费看a级黄色片| 18在线观看网站| 亚洲精品国产精品久久久不卡| 免费在线观看亚洲国产| 极品少妇高潮喷水抽搐| 人妻久久中文字幕网| 精品乱码久久久久久99久播| 超碰成人久久| 日本精品一区二区三区蜜桃| 国产精品乱码一区二三区的特点 | 在线观看www视频免费| 男人的好看免费观看在线视频 | 成人黄色视频免费在线看| 婷婷精品国产亚洲av在线 | 一边摸一边抽搐一进一小说 | 丰满饥渴人妻一区二区三| 欧美成狂野欧美在线观看| 黄色 视频免费看| 日本撒尿小便嘘嘘汇集6| 精品久久久久久久毛片微露脸| 一本大道久久a久久精品| 青草久久国产| 久久久国产欧美日韩av| 欧美另类亚洲清纯唯美| 午夜福利在线免费观看网站| 久久国产精品男人的天堂亚洲| 精品少妇久久久久久888优播| 黄片大片在线免费观看| 搡老岳熟女国产| 成年人黄色毛片网站| 黑人操中国人逼视频| 免费观看a级毛片全部| 国产三级黄色录像| 三上悠亚av全集在线观看| 伦理电影免费视频| 成人手机av| 国产精品.久久久| 妹子高潮喷水视频| 国精品久久久久久国模美| 成人国语在线视频| 国产精品久久视频播放| 美女高潮喷水抽搐中文字幕| 亚洲国产欧美日韩在线播放| 亚洲午夜理论影院| 国产成人免费观看mmmm| 俄罗斯特黄特色一大片| 午夜影院日韩av| 老司机在亚洲福利影院| xxx96com| 亚洲欧美激情综合另类| 91av网站免费观看| 搡老岳熟女国产| 天堂俺去俺来也www色官网| 欧美乱码精品一区二区三区| xxx96com| 亚洲avbb在线观看| 99精品久久久久人妻精品| 99在线人妻在线中文字幕 | 精品一区二区三卡| 亚洲午夜理论影院| 99久久精品国产亚洲精品| 色94色欧美一区二区| 精品电影一区二区在线| 成人影院久久| 老鸭窝网址在线观看| 亚洲av美国av| 一区二区三区激情视频| 在线观看舔阴道视频| 国产成人精品久久二区二区免费| 精品一区二区三区av网在线观看| 欧美精品人与动牲交sv欧美| 日韩欧美一区视频在线观看| 婷婷成人精品国产| 精品国内亚洲2022精品成人 | 在线观看免费午夜福利视频| 精品久久久久久久毛片微露脸| 精品一区二区三区视频在线观看免费 | 老司机影院毛片| 一二三四在线观看免费中文在| 精品午夜福利视频在线观看一区| 久久久水蜜桃国产精品网| 麻豆乱淫一区二区| 淫妇啪啪啪对白视频| 中文字幕制服av| videosex国产| 亚洲三区欧美一区| 天天操日日干夜夜撸| 欧美黑人欧美精品刺激| 欧美午夜高清在线| 精品一品国产午夜福利视频| 国产野战对白在线观看| 黑人巨大精品欧美一区二区蜜桃| 丝袜在线中文字幕| 国产无遮挡羞羞视频在线观看| 欧美日韩av久久| 99国产精品一区二区三区| 亚洲av片天天在线观看| xxxhd国产人妻xxx| 国产日韩一区二区三区精品不卡| 狂野欧美激情性xxxx| 国产97色在线日韩免费| 久久午夜亚洲精品久久| 中亚洲国语对白在线视频| 国产男女超爽视频在线观看| 视频区欧美日本亚洲| 黄色视频不卡| 亚洲成av片中文字幕在线观看| 欧美国产精品va在线观看不卡| 久久影院123| 成熟少妇高潮喷水视频| 99久久综合精品五月天人人| 人人妻人人澡人人看| 后天国语完整版免费观看| a级毛片在线看网站| 日韩欧美免费精品| 亚洲九九香蕉| 欧美日韩乱码在线| 91成年电影在线观看| videosex国产| 久久精品人人爽人人爽视色| 捣出白浆h1v1| 亚洲成人国产一区在线观看| av在线播放免费不卡| 亚洲国产精品sss在线观看 | 久久久精品国产亚洲av高清涩受| 亚洲人成电影免费在线| 国产熟女午夜一区二区三区| 一区二区三区精品91| 在线观看免费日韩欧美大片| 亚洲五月天丁香| 这个男人来自地球电影免费观看| 国产成人精品久久二区二区91| 欧美在线一区亚洲| 国产色视频综合| 日本wwww免费看| 国产男女内射视频| 国产高清国产精品国产三级| 免费观看精品视频网站| 亚洲一区二区三区不卡视频| 老熟妇乱子伦视频在线观看| 一边摸一边做爽爽视频免费| 每晚都被弄得嗷嗷叫到高潮| 国产野战对白在线观看| 一a级毛片在线观看| 亚洲久久久国产精品| 露出奶头的视频| 亚洲av熟女| 色精品久久人妻99蜜桃| 日韩有码中文字幕| 无人区码免费观看不卡| 精品午夜福利视频在线观看一区| 久久久国产一区二区| 国产精品 欧美亚洲| 黄色a级毛片大全视频| 高清欧美精品videossex| 欧美激情久久久久久爽电影 | 女人久久www免费人成看片| 五月开心婷婷网| 亚洲av日韩在线播放| 日韩大码丰满熟妇| 99精品久久久久人妻精品| 国产区一区二久久| 国产精品免费大片| 国产片内射在线| 人妻丰满熟妇av一区二区三区 | 亚洲欧美日韩另类电影网站| 欧美亚洲日本最大视频资源| 老司机靠b影院| av天堂在线播放| 精品一品国产午夜福利视频| 亚洲视频免费观看视频| 捣出白浆h1v1| 久久中文字幕人妻熟女| 久久精品成人免费网站| 久久香蕉激情| 99久久精品国产亚洲精品| 亚洲国产中文字幕在线视频| 国产精品一区二区免费欧美| 在线国产一区二区在线| 久久精品亚洲熟妇少妇任你| 欧美精品亚洲一区二区| 黄色 视频免费看| 在线视频色国产色| 在线免费观看的www视频| 露出奶头的视频| 交换朋友夫妻互换小说| 欧美乱妇无乱码| 99热国产这里只有精品6| 人人妻人人澡人人看| 成熟少妇高潮喷水视频| 久久国产精品人妻蜜桃| 天天躁日日躁夜夜躁夜夜| 在线观看舔阴道视频| 国产视频一区二区在线看| 国产欧美日韩一区二区三区在线| 丁香欧美五月| www.999成人在线观看| 久久人妻福利社区极品人妻图片| 最近最新免费中文字幕在线| 精品熟女少妇八av免费久了| 狂野欧美激情性xxxx| 欧美久久黑人一区二区| 无限看片的www在线观看| 日本黄色视频三级网站网址 | 人成视频在线观看免费观看| 免费观看a级毛片全部| www.精华液| 日日爽夜夜爽网站| 久久国产精品男人的天堂亚洲| 1024香蕉在线观看| 欧美日韩亚洲综合一区二区三区_| 日韩大码丰满熟妇| 亚洲成a人片在线一区二区| 国产精品永久免费网站| 成人精品一区二区免费| 下体分泌物呈黄色| 自线自在国产av| 在线永久观看黄色视频| 日韩熟女老妇一区二区性免费视频| 天天操日日干夜夜撸| 欧美人与性动交α欧美精品济南到| 中文字幕色久视频| 91老司机精品| 99国产精品一区二区三区| 老汉色∧v一级毛片| 日韩精品免费视频一区二区三区| 91精品三级在线观看| 亚洲国产精品一区二区三区在线| 亚洲人成电影免费在线| 亚洲国产中文字幕在线视频| 午夜91福利影院| 999久久久精品免费观看国产| 精品一区二区三卡| 国产一区二区激情短视频| 夜夜爽天天搞| 久9热在线精品视频| 欧洲精品卡2卡3卡4卡5卡区| 看片在线看免费视频| 亚洲欧美激情在线| 天天影视国产精品| 国产亚洲欧美在线一区二区| 亚洲一区二区三区欧美精品| 欧美人与性动交α欧美精品济南到| 国产精品久久电影中文字幕 | 欧美亚洲日本最大视频资源| 国产主播在线观看一区二区| 国产精品一区二区在线观看99| 91国产中文字幕| 成熟少妇高潮喷水视频| 亚洲中文日韩欧美视频| 亚洲熟女毛片儿| 一a级毛片在线观看| 国产精品自产拍在线观看55亚洲 | 日韩免费高清中文字幕av| 老熟妇仑乱视频hdxx| 他把我摸到了高潮在线观看| 51午夜福利影视在线观看| 国产极品粉嫩免费观看在线| 久久影院123| 欧美不卡视频在线免费观看 | 亚洲熟女精品中文字幕| 国产精品美女特级片免费视频播放器 | 国产欧美日韩一区二区三| 国产av精品麻豆| 欧美日韩一级在线毛片| 欧美日韩亚洲国产一区二区在线观看 | av网站在线播放免费| 亚洲午夜理论影院| 身体一侧抽搐| 国产精品久久久人人做人人爽| 最近最新中文字幕大全免费视频| 极品人妻少妇av视频| 我的亚洲天堂| 欧美激情久久久久久爽电影 | 在线看a的网站| 精品一区二区三区四区五区乱码| 国产在线精品亚洲第一网站| 精品少妇一区二区三区视频日本电影| 日本欧美视频一区| 少妇裸体淫交视频免费看高清 | 99久久99久久久精品蜜桃| 黄色 视频免费看| 久热这里只有精品99| 国产伦人伦偷精品视频| 欧美日韩成人在线一区二区| 国产成人啪精品午夜网站| 搡老熟女国产l中国老女人| 妹子高潮喷水视频| 亚洲,欧美精品.| 亚洲五月天丁香| 国产单亲对白刺激| 天天躁日日躁夜夜躁夜夜| 成人国产一区最新在线观看| 99精品欧美一区二区三区四区| 久久久久久久久久久久大奶| 91麻豆av在线| 亚洲自偷自拍图片 自拍| 亚洲成人免费av在线播放| 成人亚洲精品一区在线观看| 如日韩欧美国产精品一区二区三区| svipshipincom国产片| 9热在线视频观看99| 天天影视国产精品| 成年人黄色毛片网站| 国产极品粉嫩免费观看在线| 亚洲在线自拍视频| 日本黄色视频三级网站网址 | 欧美日韩亚洲国产一区二区在线观看 | 午夜免费鲁丝| 亚洲中文字幕日韩| 亚洲黑人精品在线| 亚洲av日韩精品久久久久久密| 成年人黄色毛片网站| 人人妻人人添人人爽欧美一区卜| 欧美日韩亚洲高清精品| 视频区欧美日本亚洲| 成人三级做爰电影| 亚洲欧美日韩另类电影网站| www.999成人在线观看| 日韩一卡2卡3卡4卡2021年| 色婷婷久久久亚洲欧美| 亚洲av美国av| 久99久视频精品免费| 欧美国产精品va在线观看不卡| 国产精品免费一区二区三区在线 | 中文字幕高清在线视频| 美女视频免费永久观看网站| 天天添夜夜摸| 成人18禁在线播放| 精品电影一区二区在线| 国产精品 国内视频| 不卡av一区二区三区| 国内毛片毛片毛片毛片毛片| 久久精品国产综合久久久| 国产免费现黄频在线看| netflix在线观看网站| 国产亚洲欧美在线一区二区| 丝袜美足系列| 久久精品人人爽人人爽视色| 99久久99久久久精品蜜桃| 久久精品国产清高在天天线| 国产97色在线日韩免费| 亚洲少妇的诱惑av| 精品久久久久久电影网| 精品少妇一区二区三区视频日本电影| 黄频高清免费视频| a级毛片黄视频| 在线看a的网站| 免费在线观看视频国产中文字幕亚洲| 狠狠狠狠99中文字幕| 欧美性长视频在线观看| 人成视频在线观看免费观看| 黄色成人免费大全| 嫁个100分男人电影在线观看| 久久国产精品影院| 国产精华一区二区三区| 免费在线观看影片大全网站| 18禁黄网站禁片午夜丰满| 中文字幕制服av| 在线观看日韩欧美| 99久久人妻综合| 免费在线观看亚洲国产| 久久久国产欧美日韩av| 欧美成人免费av一区二区三区 | 国产精品秋霞免费鲁丝片| 久久人人97超碰香蕉20202| 国产精品美女特级片免费视频播放器 | 国产成人免费无遮挡视频| 岛国毛片在线播放| 久久精品国产综合久久久| 成在线人永久免费视频| 欧美 亚洲 国产 日韩一| 黄频高清免费视频| 老司机在亚洲福利影院| 18在线观看网站| 黄网站色视频无遮挡免费观看| 男人的好看免费观看在线视频 | 19禁男女啪啪无遮挡网站| 国精品久久久久久国模美| 美女高潮到喷水免费观看| 日日夜夜操网爽| 国产亚洲欧美在线一区二区| 在线国产一区二区在线| 啦啦啦 在线观看视频| 国产精品免费大片| 91麻豆精品激情在线观看国产 | 身体一侧抽搐| 国产av又大| 成人av一区二区三区在线看| 精品久久久久久电影网| 国产在线观看jvid| 国产精品永久免费网站| 亚洲第一青青草原| 国精品久久久久久国模美| 女性生殖器流出的白浆| 在线观看舔阴道视频| 少妇裸体淫交视频免费看高清 | 黄色 视频免费看| 午夜精品久久久久久毛片777| ponron亚洲| 欧美日韩一级在线毛片| 极品少妇高潮喷水抽搐| 人妻久久中文字幕网| 一区福利在线观看| 欧美精品av麻豆av| 亚洲熟妇熟女久久| 欧美乱码精品一区二区三区| 国产成人精品无人区| 国产精品98久久久久久宅男小说| 91国产中文字幕| 飞空精品影院首页| 亚洲成a人片在线一区二区| 亚洲国产精品sss在线观看 | 久久国产亚洲av麻豆专区| 国产一区二区三区在线臀色熟女 | 天天影视国产精品| 男女高潮啪啪啪动态图| 亚洲av成人av| 91在线观看av| 交换朋友夫妻互换小说| 国产亚洲一区二区精品| 黄色片一级片一级黄色片| 极品少妇高潮喷水抽搐| 美女扒开内裤让男人捅视频| 欧美老熟妇乱子伦牲交| 国产1区2区3区精品| 99精国产麻豆久久婷婷| 曰老女人黄片| 欧美精品一区二区免费开放| 韩国精品一区二区三区| 国产成人欧美在线观看 | 啦啦啦在线免费观看视频4| 五月开心婷婷网| 香蕉久久夜色| a级毛片在线看网站| 激情视频va一区二区三区| 在线观看免费午夜福利视频| 天天躁狠狠躁夜夜躁狠狠躁| 高清在线国产一区| 在线天堂中文资源库| 精品国产国语对白av| 亚洲精品中文字幕一二三四区| av免费在线观看网站| 中文欧美无线码| 伊人久久大香线蕉亚洲五| 首页视频小说图片口味搜索| 少妇被粗大的猛进出69影院| 精品少妇久久久久久888优播| 成人亚洲精品一区在线观看| 岛国在线观看网站| 成人手机av| 韩国精品一区二区三区| 久久久国产一区二区| 国产一卡二卡三卡精品| 看免费av毛片| 中文字幕人妻丝袜制服| 女同久久另类99精品国产91| 一进一出抽搐动态| 日韩 欧美 亚洲 中文字幕| 欧美精品人与动牲交sv欧美| av免费在线观看网站| 制服诱惑二区| 国产在线一区二区三区精| 99久久综合精品五月天人人| 亚洲精品一二三| 真人做人爱边吃奶动态| 亚洲精品中文字幕一二三四区| 亚洲五月色婷婷综合| 久久午夜亚洲精品久久| 十八禁高潮呻吟视频| 国产精品免费大片| 天堂动漫精品| 咕卡用的链子| av超薄肉色丝袜交足视频| 身体一侧抽搐| 亚洲免费av在线视频| 久久久精品免费免费高清| 人妻丰满熟妇av一区二区三区 | 99精品在免费线老司机午夜| 欧美激情高清一区二区三区| 又紧又爽又黄一区二区| 国产精品美女特级片免费视频播放器 | 欧美日韩av久久| 久久午夜亚洲精品久久| 女人高潮潮喷娇喘18禁视频| 日韩欧美在线二视频 | 老司机深夜福利视频在线观看| 麻豆av在线久日| 纯流量卡能插随身wifi吗| 狂野欧美激情性xxxx| 宅男免费午夜| 亚洲av成人一区二区三| 一级片免费观看大全| 12—13女人毛片做爰片一| 首页视频小说图片口味搜索| 久久天堂一区二区三区四区| 美女福利国产在线| 亚洲欧美一区二区三区黑人| 精品国产超薄肉色丝袜足j| 免费不卡黄色视频| 丰满的人妻完整版| 亚洲中文日韩欧美视频| 啪啪无遮挡十八禁网站| 国产高清视频在线播放一区| 99久久综合精品五月天人人| 男女床上黄色一级片免费看| av电影中文网址| 18在线观看网站| 老熟女久久久| 男女高潮啪啪啪动态图| 久久久久久久午夜电影 | 日韩欧美一区二区三区在线观看 | 男女床上黄色一级片免费看| av电影中文网址| 日韩大码丰满熟妇| av天堂在线播放| 亚洲精品国产色婷婷电影| 视频区欧美日本亚洲| 亚洲色图av天堂| 精品熟女少妇八av免费久了| 午夜福利视频在线观看免费| 男女高潮啪啪啪动态图| 人人妻人人澡人人爽人人夜夜| 成年版毛片免费区| 国产欧美日韩综合在线一区二区| 日韩精品免费视频一区二区三区| 国产精品电影一区二区三区 | 欧美亚洲日本最大视频资源| 狂野欧美激情性xxxx| 欧美黑人欧美精品刺激| 国产精品久久久久成人av| 欧美精品一区二区免费开放| 亚洲国产毛片av蜜桃av| 亚洲欧美日韩高清在线视频| 亚洲七黄色美女视频| 免费在线观看影片大全网站| 人妻久久中文字幕网| 成年人免费黄色播放视频| 美女国产高潮福利片在线看| 在线观看日韩欧美| 午夜福利乱码中文字幕| 一区二区三区国产精品乱码| 麻豆av在线久日| 精品人妻1区二区| 久久午夜综合久久蜜桃| 香蕉丝袜av| 精品国产亚洲在线| 中出人妻视频一区二区| 啦啦啦 在线观看视频| 老熟妇乱子伦视频在线观看| 国产精品美女特级片免费视频播放器 | 国产精品99久久99久久久不卡| bbb黄色大片| 少妇被粗大的猛进出69影院| 欧美日韩国产mv在线观看视频| 最新的欧美精品一区二区| 涩涩av久久男人的天堂| 精品一区二区三区av网在线观看| 亚洲av日韩在线播放| 三上悠亚av全集在线观看| 国产精品 欧美亚洲| 天天添夜夜摸| 激情视频va一区二区三区| 亚洲第一欧美日韩一区二区三区| 日韩三级视频一区二区三区| 九色亚洲精品在线播放| av片东京热男人的天堂| 人人妻人人爽人人添夜夜欢视频| 亚洲va日本ⅴa欧美va伊人久久| 久久久国产成人免费| 久久亚洲精品不卡| 黑人操中国人逼视频| 国产色视频综合| 午夜福利欧美成人| 色综合婷婷激情| 国产精华一区二区三区| 精品一区二区三卡| 看黄色毛片网站| 男人的好看免费观看在线视频 | 久久精品91无色码中文字幕| 国产日韩欧美亚洲二区| 精品国产乱子伦一区二区三区| 99国产精品一区二区蜜桃av | 国产精品 欧美亚洲| 99精国产麻豆久久婷婷| 身体一侧抽搐| 国产精品欧美亚洲77777| 大陆偷拍与自拍| 婷婷成人精品国产| √禁漫天堂资源中文www| 国产av精品麻豆| 丝袜在线中文字幕| 丰满人妻熟妇乱又伦精品不卡| 久久中文字幕人妻熟女| 12—13女人毛片做爰片一| 国产精品永久免费网站| 国产亚洲精品久久久久5区| 中文字幕制服av| 中文字幕最新亚洲高清| 亚洲欧美一区二区三区黑人| 亚洲 国产 在线| 欧美久久黑人一区二区| 交换朋友夫妻互换小说| 免费看a级黄色片| 欧美最黄视频在线播放免费 | 国产精品秋霞免费鲁丝片| 人妻一区二区av| 手机成人av网站| 夜夜夜夜夜久久久久| 在线天堂中文资源库| 欧美成人午夜精品| 久久久国产一区二区| 亚洲性夜色夜夜综合| 12—13女人毛片做爰片一| 中文字幕最新亚洲高清|