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

    跨流區(qū)超低軌航天器快速氣動力計算方法

    2022-12-03 16:38:20郭晨林陳方趙艷彬尤超藍錢勇李文龍楊麗麗
    上海航天 2022年5期
    關鍵詞:氣動力飛行器氣動

    郭晨林,陳方,趙艷彬,尤超藍,錢勇,李文龍,楊麗麗

    (1.上海交通大學 航空航天學院,上海 200240;2.上海衛(wèi)星工程研究所研發(fā)中心,上海 201109)

    0 引言

    超低軌飛行器由于其重訪周期短、對地成像觀測性能強的特點,在通信、對地探測等領域展現(xiàn)出了巨大的潛力[1],但超低軌飛行器受到的氣動攝動大,造成軌道維持與姿態(tài)控制的困難。同時,由于其運行高度在120~300 km 之間,屬于過渡流區(qū)與自由分子流區(qū),對于氣動力的估算缺少經驗理論與工具,造成了總體設計與動力選型上的困難。因此,當前急需有關超低軌衛(wèi)星的快速氣動力計算工具,以實現(xiàn)合理的動力選型,延長衛(wèi)星工作壽命[2]與加速超低軌飛行器的初步設計迭代。

    對于工作在自由分子流區(qū)中的飛行器的氣動力特性求解,BIRD[3]提出的直接模擬蒙特卡洛法(Direct Simulation Monte Carlo,DSMC)求解分子之間的碰撞效應高精度求解其氣動力特性,DAVIS[4]提出的測試粒子蒙特卡洛法(Test Particle Monte Carlo,TPMC)省略了分子間的碰撞效應,雖然降低求解精度,但是大幅減小了計算量,提升了計算效率。Maxwell 通過建立簡單的面元與氣體粒子的碰撞模型,求解自由分子流區(qū)中的面元受力狀況,進而積分求解飛行器受力[5]。該方法成為自由分子流區(qū)中快速氣動力計算的基本方法,后續(xù)許多學者根據分子與面元碰撞時的動量損失與漫反射效應對Maxwell 提出的模型進行優(yōu)化,擴展了該模型的精度與適用范圍[6-7]。

    在過渡流區(qū)中,氣體的運動規(guī)律更為復雜,難以抽象出簡化的物理模型。當前求解過渡流區(qū)中的流動特性多依賴于DSMC 方法[8]。由于在該區(qū)域中大氣密度較高,意味著DSMC 模擬的分子數(shù)量多,所需要的計算量大,因此,求解過渡流區(qū)氣動特性的方法為:首先對連續(xù)流區(qū)的氣動特性與自由分子流區(qū)中的氣動特性進行求解,再使用橋函數(shù)加權來近似過渡流區(qū)的氣動特性。常用的橋函數(shù)有sine-square 橋函數(shù)[9]與erf-log 橋函數(shù)[10],通過假定多個待定系數(shù)函數(shù)擬合連續(xù)流區(qū)-過渡流區(qū)-自由分子流區(qū)的部分計算結果,實現(xiàn)對其他工況的結果進行估計。以上橋函數(shù)的擬合只對特定的氣流來流條件與當前幾何外形的飛行器適用,當外部條件發(fā)生變化時,需要新的DSMC 或試驗結果進行標定。

    目前,過渡流區(qū)的計算高度依賴于DSMC,TPMC 等高精度計算工具,快速氣動力計算工具十分有限,常用的氣動估算方法存在一定的局限性。同時,面元法應用的前景廣泛,通過對面元法進行進一步改進,以上問題能被較好的解決。本文旨在解決超低軌飛行器的快速氣動計算需求,提出一套適用于超低軌衛(wèi)星從過渡流區(qū)到自由流區(qū)的快速氣動估算方法。在面元法的基礎上,對迎風單元篩選算法與橋函數(shù)進行了設計,達到了提高精度,增加適用范圍的效果。

    1 超低軌衛(wèi)星快速氣動力計算方法

    本文提出的超低軌衛(wèi)星快速氣動力計算方法主要分為3 個部分,如圖1 所示。包括輸入數(shù)據、數(shù)據的預處理方式和氣動力計算方法。其中,有2 個過程顯著影響計算精度和計算效率:①數(shù)據預處理與迎風單元與背風單元的劃分算法;②氣動力的計算方法和橋函數(shù)模型。本文將對此進行較詳細的闡述。

    圖1 超低軌衛(wèi)星快速氣動力計算方法流程Fig.1 Flow chart of the rapid aerodynamic calculation method for ultra-low orbit spacecrafts

    1.1 數(shù)據預處理

    在輸入的模型STL 格式文件中,包含著衛(wèi)星的幾何尺寸信息,首先需要對其表面單元進行離散,得到其表面的三角形單元如圖2 所示。在建模的過程中,常以飛行器的體坐標系進行建模,為方便后續(xù)的迎風背風單元算法及氣動力計算,需要對離散單元在風軸系和體軸系下進行自由轉換。通過以下旋轉矩陣可以實現(xiàn)轉換過程:

    圖2 迎風單元與背風單元篩選Fig.2 Screening diagram of windward and leeward units

    式中:Cbo為體坐標系ObXbYbZb至風軸坐標系OwXwYwZw下的旋轉矩陣;α、β分別為當前飛行器的迎角、側滑角。

    進行坐標轉換后,需對表面的單元進行分類,分為迎風單元與背風單元。根據單元屬性的不同,在后續(xù)的氣動力計算中使用不同的計算方式。本文設計了可用于非凸表面的迎風單元篩選算法。通過驗證發(fā)現(xiàn),隨著網格密度達到一定要求時,其具有較高的計算精度,滿足對復雜外形的超低軌衛(wèi)星氣動力估算的需求。迎風單元和背風單元篩選算法的輸入包含如下數(shù)據:

    1)節(jié)點位置Pi=[xi,yi,zi],為編號為i節(jié)點其對應在風軸系下的空間坐標位置;

    2)單元列表Un=[u,v,w],n為單元編號,u,v,w為三角形單元的3 個頂點節(jié)點編號。

    算法的輸出包含以下幾個部分:

    1)迎風單元列表Wn=[u',v',w'],u',v',w'為三角形單元的3 個頂點節(jié)點編號;

    2)背風單元列表Ln=[u',v',w'],u',v',w'為三角形單元的3 個頂點節(jié)點的編號;

    3)迎風節(jié)點列表PWi=[xi,yi,zi],xi,yi,zi為節(jié)點在風軸系下的坐標位置。

    具體算法見表1。

    表1 迎風單元篩選算法流程Tab.1 Flow chart of the screening algorithm for windward units

    從迎風方向依次選取節(jié)點,判斷其是否在迎風單元列表的投影內,若不在任意投影內,則將該節(jié)點加入迎風節(jié)點列表中,并將與其連接的單元加入迎風單元列表。最后將所有總單元列表Un減去篩選出的迎風單元列表,則可得到背風單元列表。具體篩選過程如圖2 所示,在對表面進行離散后的衛(wèi)星網格中,黑色單元為篩選出的迎風單元,灰色的是待篩選的單元列表。紅色節(jié)點與藍色節(jié)點為兩個待篩選節(jié)點。在圖示的沿著來流方向的投影平面上,紅色節(jié)點在黑色單元投影外,所以其為新的迎風單元;而藍色節(jié)點在當黑色單元的投影內,其為背風節(jié)點。

    為驗證該算法的準確性,以及提出在后續(xù)計算中劃分網格的精細程度,本文設計了如下算例,如圖3(a)所示,該圖所示的梨狀外形,為一個非凸的幾何結構,其在軸向投影方向的最大投影半徑R=0.5 以此作為模型的特征尺度,因此可以得到其準確的投影面積Sproj=πR2=0.25π。本文依次增加特征尺度與三角形網格大小R/δx之間的比值,觀察通過迎風單元與背風單元算法計算出的投影面積與實際投影面積的差距,驗證算法精度及求解合適的特征尺度與網格大小的比值關系R/δx。計算的結果如圖3(b)所示,發(fā)現(xiàn)在比值為10 時,算法計算出的投影面積與實際的投影面積十分接近,進一步縮小網格尺度,增加特征尺度與網格大小的比值,精度不再改善。因此,本文認為該算法的準確性良好,在對后續(xù)的計算模型進行表面單元的離散過程中,將選取特征尺度與網格大小的比值R/δx=10。

    圖3 迎風單元與背風單元篩選算例模型Fig.3 Example model for screening windward and leeward units

    1.2 氣動力計算

    在氣動力的計算過程中,需要從大氣模型中根據低軌衛(wèi)星的飛行高度和經緯度獲得當前的流動條件,獲得準確的流動條件對算法的精度與可靠性也有較大的影響。因此,選取合適的大氣模型也至關重要。當前常用的大氣模型有以下幾種,US Standard 1976、MSISe-90、MSISe-00 和Jacchia-71等[11-13]。本文選取NRLMSISE-00 作為計算使用的大氣模型,主要基于其擁有以下良好性質[14]:①使用廣泛,有豐富的程序接口可調用;②能夠實現(xiàn)努森數(shù)(Knudsen Number,Kn)與海拔高度的轉換;③其數(shù)據豐富,滿足本文的計算需求。

    本文選取的模型計算條件見表2。

    表2 大氣模型參數(shù)Tab.2 Parameters of the atmospheric model

    為方便后續(xù)計算,本文固定部分大氣模型參數(shù)(見表2),只保留高度? 作為變量。以上數(shù)據取太陽豐度為平均值時的大氣系數(shù),能反映當前的大氣較為普遍的狀態(tài)。選取以上參數(shù)后,本文假設大氣模型為函數(shù)Fa,其輸入變量為高度?,則其可輸出Kn、T、m、ρ分別為當前高度下的努森數(shù)、大氣溫度、大氣的摩爾質量與大氣密度,即:

    根據不同高度對應的流區(qū)的不同,本文將計算分為3 個部分,見表3。

    表3 對應不同流區(qū)適用的不同氣動力計算方式Tab.3 Aerodynamic calculation methods for different flow regions

    1.2.1 自由分子流模型

    在自由分子流中,分子之間的碰撞概率很小,經碰撞后反射分子對流動的影響可以忽略不計[6]。因此,主要關注的過程是氣體分子與衛(wèi)星表面發(fā)生碰撞(見圖4),分子與衛(wèi)星表面發(fā)生能量與動量的交換情況。因此,自由分子流動時衛(wèi)星的受力狀況是由氣體-表面相互作用性質決定的。在假設氣體僅與衛(wèi)星表面發(fā)生動量交換的前提下,作用在衛(wèi)星表面上的力等于氣體動量的變化率。對于一個外凸形狀的衛(wèi)星,其單位面積受力大小與氣體的入射動量與反射動量的大小有關。以一個衛(wèi)星的表面單元為例,其受到分子的作用力可表示為

    圖4 氣體分子與單元表面碰撞模型Fig.4 Model for the collision of gas molecule and unit surface(a)Maxwell 模型(b)牛頓理論

    式中:f為分子的作用力;A為衛(wèi)星表面的單位面積;p為法向的動量通量;τ為切向的動量通量。下標a、b 分別對應著入射通量和反射通量,法向和切向入射動量流pa、τa取決于入射速度(V)和質量流(dQ);其反射動量通量為-pbn+τbt;n為在物面表面的法向分量;t為在物面表面的切向分量。

    在自由分子流中,可以認為

    式中:dQ=ρVcosδ,δ為粒子入射方向與平面之間的夾角,ρ為分子流密度;V為分子流速度大小。

    在流動中氣體分子的速度服從Maxwell 速度分布函數(shù)[15]。則入射動量通量可以表示為

    式中:m為單個氣體分子的質量;F(u)為麥克斯韋分布函數(shù)。

    式中:s為分子速度比;V為在速度分布函數(shù)下最可能的速度大?。籖c為氣體常數(shù);T∞為當前高度下的來流溫度。

    目前,如何求解確定反射通量是個難點。有許多采用簡化的氣體-表面作用模型來對反射動量通量進行建模,如Maxwell[15]、Schamberg、Schaaf&Chambre[6]等模型。本 文選取Sentman 模型[16]進行建模,假設一個粒子在撞擊后會依據一定的概率分布向不同方向進行漫反射,定義分子與平面發(fā)生碰撞之后,平面的法向力系數(shù)與切向力系數(shù)分別為

    式中:Tk,b、Tk,a為

    式中:αacc為能量調節(jié)系數(shù),根據飛行器的高度確定[16],隨著飛行高度的增加,能量調節(jié)系數(shù)的取值不斷減??;Tw為飛行器的表面壁溫。

    1.2.2 連續(xù)流區(qū)模型

    牛頓理論(牛頓正弦平方律)闡述了一種推導平面與流體相互作用的理論模型,其假設氣體分子為質點且互相孤立,氣體與平面撞擊后,其動量在物面的法向上發(fā)生完全的動量交換,分子沿著物面的切向動量完全保留。其理論在低速流動條件下偏差較大,但物體在高速連續(xù)流區(qū)運動中,其對氣動力的估算效果良好。物面受到的法向力系數(shù)CP為

    牛頓理論提出的模型僅僅與氣流方向與物面的夾角有關,在馬赫數(shù)不夠大時,會出現(xiàn)明顯的偏差。為了對該偏差進行修正,Lester Lee 提出了新的計算方式:

    其中:

    式中:Ma∞為自由來流的馬赫數(shù)。

    在牛頓理論中,背風單元不受氣動力的作用,但在連續(xù)流區(qū)中,忽略背風單元的氣動力會導致計算出現(xiàn)明顯的誤差。因此,本文使用以下普朗特-邁耶膨脹波理論估算背風單元的氣動系數(shù)[17]

    式中:γ為絕熱氣體指數(shù),在量熱完全氣體的假設下γ為1.4。

    在實際估算中,以述方法能夠較好地估計連續(xù)流區(qū)中的氣動力系數(shù)值,但由于缺少黏性力的估算,依然會有一定的偏差,為做進一步修正,本文使用可壓縮流的摩擦系數(shù)計算公式對摩擦力系數(shù)進行估算。

    1.2.3 過渡流區(qū)模型

    在過渡流區(qū)的氣動力計算中,本文使用橋函數(shù)的方式進行計算,即在計算出連續(xù)流區(qū)與自由分子流區(qū)的氣動特性參數(shù)后,根據在過渡流區(qū)中的高度h由大氣模型推導出當前高度下的努森數(shù)Kn對連續(xù)流區(qū)與自由分子流區(qū)的氣動特性參數(shù)進行權重分配,得到當前位置的氣動特性參數(shù)。

    式中:Cfm為在自由分子流中求得的氣動力參數(shù);Ccont為在連續(xù)流區(qū)中求得的氣動力參數(shù);Pb為橋函數(shù)。當前常用的橋函數(shù)有sine-squared 橋函數(shù)與log-erf 橋函數(shù)。

    在設定合理的參數(shù)條件下,2 種函數(shù)均表現(xiàn)出良好的擬合性能。以上2 種函數(shù)都需要根據求解物體的幾何外形的變化與來流條件的改變對相關參數(shù)進行修正。在實際使用過程中,需要使用部分實驗數(shù)據或高精度的求解結果對參數(shù)進行調整,這增加了工作量。

    為了對上述問題進行改善,本文提出了以logistic 函數(shù)為基礎的橋函數(shù),并在函數(shù)中加入了形狀修正因子,拓展了橋函數(shù)的適用范圍。本文設計的橋函數(shù)如下:

    式中:q為形狀因子,其與氣體來流方向與物面間的夾角有關;橋函數(shù)為與當?shù)豄n數(shù)和形狀有關的函數(shù)。

    2 模型驗證

    為驗證提出的快速算法的有效性及正確性,本文結合DSMC 的計算結果,對以上方法進行驗證。吳子牛等[10]曾對圓柱體與鈍頭雙楔體進行DSMC與橋函數(shù)的研究,其研究對象如圖5 所示,2 個模型展向上的長度均為1 000 mm。

    圖5 幾何外形Fig.5 Geometric shape

    計算條件見表4,給定模型壁溫分別500 K、300 K;自由來流溫度條件給定為300 K。選取圓柱模型的速度工況為馬赫數(shù)4 與馬赫數(shù)16,鈍頭雙楔體模型速度工況為馬赫數(shù)4 與馬赫數(shù)8;選取計算努森數(shù)Kn為10-3~102,在這個范圍內大氣特性變化顯著。其余輸入條件由Nrlmsise-00 大氣模型給出。

    表4 模型驗證的計算條件Tab.4 Calculation conditions for model verification

    為了進行對比,本文對DSMC 結果及使用本文提出的logistic-log 函數(shù)與sine-squared 函數(shù)進行比較。圖6 表示了在本節(jié)所研究兩種不同模型中,本文所提出的logistic-log 橋函數(shù)與sine-squared 橋函數(shù)在2 種不同模型中不同條件下的其努森數(shù)(Kn)與阻力系數(shù)(CD)之間的關系如圖6 所示。在圖6(a)中,紅色線代表的是Ma=16 結果,藍色線代表的是Ma=4 結果。其中,粗虛線為logistic-log 橋函數(shù)的計算結果,而細虛線為sine-squared 橋函數(shù)的計算結果,三角符號代表著使用自由分子流算法(Free Molecular Fuction,F(xiàn)MF)的計算結果,其余圖像中的散點代表的是DSMC 計算數(shù)據。由圖6 可知,對于圓柱的計算工況,在Ma=16 的高速工況下,logisticlog 橋函數(shù)與sine-squared 橋函數(shù)與DSMC 結果吻合都較好。在Ma=4 的工況下,sine-squared 橋函數(shù)的計算結果在Kn數(shù)較高時與實驗數(shù)據產生了較為明顯的差異。在圖6(b)中,紅色線代表的是Ma=8結果,藍色線代表的是Ma=4 結果。其中,粗虛線為logistic-log 橋函數(shù)的計算結果,而細虛線為sinesquared 橋函數(shù)的計算結果,三角符號代表著使用自由分子流算法(Free Molecular Fuction,F(xiàn)MF)的計算結果,其余圖像中的散點代表的是DSMC 計算數(shù)據。同樣可以發(fā)現(xiàn),對于鈍頭雙楔體的計算工況,在速度較高的情況下(Ma=8),2 種橋函數(shù)與實驗結果基本吻合。當Ma=4 時,使用sine-squared 橋函數(shù)在Kn=12 時出現(xiàn)了最大值,與DSMC 數(shù)據相差較大,也違背了隨著Kn數(shù)的增加,氣動阻力系數(shù)也隨之增加的基本規(guī)律。以上結果說明了在固定系數(shù)的情況下,sine-squared 橋函數(shù)無法滿足在不同外形,多種工況下均保持良好的估算精度;而logistic-log 橋函數(shù)可以實現(xiàn)在不改變系數(shù)的情況下滿足對多種工況的氣動系數(shù)估算。

    圖6 不同橋函數(shù)對應的氣動特性估算結果Fig.6 Results of the aerodynamic characteristics corresponding to different bridge functions

    為驗證本文提出的形狀因子對其結果的影響,本 文對q=1+0.632 6cosδ+6.124 9與q=1 進行對比,計算結果如圖7 所示。

    在圖7(a)中,紅色線圖代表的是Ma=16 結果,藍色線代表的是Ma=4 結果。粗虛線為形狀因子q=1+0.632 6cosδ+6.124 9 的計算結果,而細虛線為形狀因子q=1 的計算結果,圖像中的散點代表的是DSMC 計算數(shù)據。由圖7 可知,對于圓柱的計算工況,在Ma=16 的情況下,2 種形狀因子與實驗結果的吻合情況較好。在Ma=4 的計算條件下,不使用q=1+0.632 6cosδ+6.124 9 的形狀因子則會出現(xiàn)明顯的誤差。對于鈍頭雙楔體的計算工況,發(fā)現(xiàn)形狀因子在不同工況下的影響不大。因為在鈍頭雙楔體中,表面外形的傾角相對固定,形狀因子對其影響則不明顯;而在圓柱工況中,圓柱各個位置的傾角均有變化,則形狀因子在速度較低的情況下起到了明顯的修正作用。在速度較低時,形狀因子起到的修正效果明顯,在高速情況下形狀因子的修正作用不顯著。為了進一步量化計算結果與DSMC 數(shù)據的吻合程度,本文引入了均方誤差對計算結果進行了比較,其計算公式如下:

    圖7 形狀因子對計算結果的影響Fig.7 Effects of the shape factor on the calculation results

    式中:ε為誤差值為DSMC 的采樣數(shù)據為使用不同橋函數(shù)得到的數(shù)值結果數(shù)據;k為DSMC 的數(shù)據維度。

    通過對均方誤差進行比較,計算出不同方法的誤差,見表5。可見本文提出的計算方式誤差最小。

    表5 不同橋函數(shù)計算結果均方誤差(ε)Tab.5 Mean square errors of the calculation results of different bridge functions

    綜上所述,本文提出logistic-log 橋函數(shù)與當前常用的橋函數(shù)相比,具有適用范圍廣,估算精度高的特點。本文提出的形狀因子q能夠根據計算外形對計算結果進行修正。

    3 典型超低軌飛行器氣動力計算

    歐洲空間局(European Space Agency,ESA)的地球重力與海洋環(huán)流探測衛(wèi)星GOCE 為典型的超低軌飛行器[18],有關學者也使用了DSMC 或TPMC方法對其氣動特性進行了研究[19-20],本文第2 章提出的方法對其在工作高度范圍內的阻力系數(shù)進行測算。

    在不同側滑角下,不同攻角下隨著工作高度的變化GOCE 的阻力系數(shù)與阻力大小的變化趨勢如圖8 所示。隨著工作高度的增加,阻力系數(shù)的變化范圍越小,GOCE 在俯仰與偏航時所受的氣動力變化相對較低。但是隨著運行高度降低,大氣密度顯著增加,阻力呈幾何量級增大。在平飛狀態(tài)下,阻力的最大值為42.526 N,最小值為0.005 41 N。因此,超低軌衛(wèi)星在進行工作時,需要及時進行姿軌控制維持其高度。

    圖8 不同高度下GOCE 的阻力與阻力系數(shù)大小與其姿態(tài)角度關系Fig.8 Relationships between the drag or drag coefficient of GOCE and its attitude angle at different altitudes

    續(xù)圖8 不同高度下GOCE 的阻力與阻力系數(shù)大小與其姿態(tài)角度關系ContinuedFig.8 Relationships between the drag or drag coefficient of GOCE and its attitude angle at different altitudes

    GOCE 形態(tài)細長,姿態(tài)變化導致的氣動力變化明顯,在運行高度在100 km 時,其阻力系數(shù)變化范圍在3.2~9.8;在高度為260 km 時其阻力系數(shù)的變化范圍在4.1~9.5。其太陽能帆板占主要浸濕面積,在偏航時,迎風面積變化明顯,導致了其在側滑時的阻力系數(shù)與阻力變化明顯高于俯仰時的情況。

    4 結束語

    1)本文提出了一種超低軌飛行器的快速氣動力計算工具方法。使用面元法進行求解,通過將飛行器進行便面單元離散為迎風單元與背風單元,分別求解不同種類單元的氣動力系數(shù),通過對離散單元的數(shù)值結果進行積分得到了飛行器的氣動力特性。通過重新設計橋函數(shù),降低計算方法的數(shù)據依賴程度,提高了應用范圍。

    2)使用2 種橋函數(shù)計算,并對比鈍頭雙楔體與圓柱體的DSMC 結果,說明了算法的可靠性和準確性,解釋了本文提出的在計算不同工況氣動力時,形狀因子對計算結果有修正作用。

    3)分析闡述了典型超低軌飛行器GOCE 在100~260 km 工作范圍內隨著俯仰角與偏航角變化氣動特性的變化趨勢。發(fā)現(xiàn)在發(fā)生姿態(tài)變化時的阻力系數(shù)變化顯著,其變化范圍在3.2~9.8,隨著高度的降低阻力系數(shù)變化的范圍越大。因此,飛行器在低軌工作時,需要為飛行器設計舵面或陀螺儀對其姿態(tài)進行控制。在過渡流區(qū)內,隨著高度增加,氣動阻力減小,阻力最大值為42.526 N,最小值為0.005 41 N。由此可見,超低軌飛行器需要在寬域過渡流區(qū)內運行,需要額外的動力裝置。

    猜你喜歡
    氣動力飛行器氣動
    中寰氣動執(zhí)行機構
    高超聲速飛行器
    基于NACA0030的波紋狀翼型氣動特性探索
    飛行載荷外部氣動力的二次規(guī)劃等效映射方法
    基于反饋線性化的RLV氣動控制一體化設計
    測控技術(2018年9期)2018-11-25 07:44:24
    復雜飛行器的容錯控制
    電子制作(2018年2期)2018-04-18 07:13:25
    側風對拍動翅氣動力的影響
    神秘的飛行器
    高速鐵路接觸線覆冰后氣動力特性的風洞試驗研究
    KJH101-127型氣動司控道岔的改造
    一边摸一边做爽爽视频免费| 久久精品国产清高在天天线| 亚洲色图 男人天堂 中文字幕| 免费在线观看影片大全网站| 国产精品永久免费网站| 最近最新中文字幕大全免费视频| 91字幕亚洲| 在线观看舔阴道视频| 女警被强在线播放| 久久中文看片网| 精品电影一区二区在线| 精品久久久久久久人妻蜜臀av | 欧美日韩精品网址| 久久国产精品人妻蜜桃| 日韩一卡2卡3卡4卡2021年| 精品久久久久久久毛片微露脸| 丁香欧美五月| 欧美日韩精品网址| 九色国产91popny在线| 18禁裸乳无遮挡免费网站照片 | 国产成人啪精品午夜网站| 免费观看人在逋| 日韩欧美免费精品| 97人妻天天添夜夜摸| 搞女人的毛片| 日韩高清综合在线| 在线视频色国产色| 国产精品99久久99久久久不卡| 如日韩欧美国产精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 一级a爱片免费观看的视频| 老司机靠b影院| 国产高清视频在线播放一区| 亚洲国产欧美日韩在线播放| 色哟哟哟哟哟哟| 麻豆久久精品国产亚洲av| 动漫黄色视频在线观看| 国产亚洲精品av在线| 性色av乱码一区二区三区2| 免费在线观看影片大全网站| 亚洲专区字幕在线| 99久久精品国产亚洲精品| 一区二区三区精品91| 满18在线观看网站| 19禁男女啪啪无遮挡网站| 午夜福利欧美成人| 国产精品久久久av美女十八| 国产一级毛片七仙女欲春2 | 午夜福利高清视频| 中文字幕av电影在线播放| 色综合欧美亚洲国产小说| 亚洲成人免费电影在线观看| 精品一区二区三区视频在线观看免费| 亚洲va日本ⅴa欧美va伊人久久| 亚洲自偷自拍图片 自拍| 岛国在线观看网站| 天堂影院成人在线观看| √禁漫天堂资源中文www| 亚洲欧洲精品一区二区精品久久久| 18禁黄网站禁片午夜丰满| 久久香蕉激情| 亚洲欧美激情综合另类| 级片在线观看| 国产伦一二天堂av在线观看| 日韩中文字幕欧美一区二区| 男女下面进入的视频免费午夜 | 午夜亚洲福利在线播放| 日韩高清综合在线| 欧美色欧美亚洲另类二区 | 一边摸一边抽搐一进一出视频| 一边摸一边做爽爽视频免费| 欧美日本视频| 国产一级毛片七仙女欲春2 | 在线观看66精品国产| 色哟哟哟哟哟哟| 亚洲少妇的诱惑av| 在线观看一区二区三区| 国产一区二区三区在线臀色熟女| 国产成人影院久久av| 久久久精品国产亚洲av高清涩受| xxx96com| 成人国产综合亚洲| 妹子高潮喷水视频| 男人舔女人下体高潮全视频| videosex国产| 午夜精品久久久久久毛片777| 99久久精品国产亚洲精品| 日本三级黄在线观看| 欧美国产日韩亚洲一区| 久久精品亚洲精品国产色婷小说| 老司机深夜福利视频在线观看| 亚洲专区字幕在线| 999久久久国产精品视频| 黄频高清免费视频| 免费在线观看亚洲国产| av电影中文网址| 麻豆久久精品国产亚洲av| 天堂动漫精品| 一边摸一边抽搐一进一小说| 久久性视频一级片| 欧美色欧美亚洲另类二区 | 好男人电影高清在线观看| 精品不卡国产一区二区三区| 成人特级黄色片久久久久久久| 欧美精品亚洲一区二区| 最好的美女福利视频网| 亚洲精品国产精品久久久不卡| 亚洲av第一区精品v没综合| 国产av又大| 精品午夜福利视频在线观看一区| 可以在线观看毛片的网站| 欧美不卡视频在线免费观看 | 久久人妻福利社区极品人妻图片| 欧美人与性动交α欧美精品济南到| 久久精品国产99精品国产亚洲性色 | 一区二区三区高清视频在线| 国产不卡一卡二| 搡老岳熟女国产| 免费高清视频大片| av视频在线观看入口| av福利片在线| 亚洲在线自拍视频| 女同久久另类99精品国产91| 久久人人精品亚洲av| 亚洲成人国产一区在线观看| 高清黄色对白视频在线免费看| 一区在线观看完整版| 成熟少妇高潮喷水视频| 成人永久免费在线观看视频| 国产精品影院久久| 一二三四社区在线视频社区8| 久久久水蜜桃国产精品网| 高清黄色对白视频在线免费看| 欧美日韩福利视频一区二区| 欧美一级a爱片免费观看看 | 国内精品久久久久久久电影| 首页视频小说图片口味搜索| 国产亚洲精品综合一区在线观看 | 啦啦啦 在线观看视频| 免费人成视频x8x8入口观看| 日韩欧美国产在线观看| www日本在线高清视频| 动漫黄色视频在线观看| av欧美777| 深夜精品福利| a级毛片在线看网站| 国产精品秋霞免费鲁丝片| 最近最新中文字幕大全电影3 | 国产成人精品无人区| 亚洲色图 男人天堂 中文字幕| aaaaa片日本免费| 午夜久久久久精精品| 免费在线观看视频国产中文字幕亚洲| 日韩欧美国产一区二区入口| 色综合站精品国产| 欧美日韩福利视频一区二区| 成人欧美大片| 99香蕉大伊视频| 成人亚洲精品一区在线观看| 精品国产一区二区久久| 麻豆一二三区av精品| 国产成人影院久久av| 国产精品九九99| 18禁观看日本| 日韩欧美国产在线观看| 国产精品亚洲av一区麻豆| 91大片在线观看| 熟妇人妻久久中文字幕3abv| 日韩免费av在线播放| 亚洲av美国av| 亚洲成人免费电影在线观看| 两性夫妻黄色片| 亚洲中文字幕一区二区三区有码在线看 | 国产欧美日韩一区二区三区在线| 亚洲午夜理论影院| 亚洲第一欧美日韩一区二区三区| 黑丝袜美女国产一区| 免费av毛片视频| 嫩草影视91久久| 精品国产亚洲在线| 欧美黄色片欧美黄色片| 嫩草影院精品99| 免费看美女性在线毛片视频| 欧美色欧美亚洲另类二区 | 亚洲美女黄片视频| 亚洲无线在线观看| 欧美色视频一区免费| 三级毛片av免费| 国产av精品麻豆| 丝袜人妻中文字幕| 欧美一级a爱片免费观看看 | 嫁个100分男人电影在线观看| 日韩大码丰满熟妇| 中亚洲国语对白在线视频| 免费看十八禁软件| 如日韩欧美国产精品一区二区三区| 国产又色又爽无遮挡免费看| 久久婷婷人人爽人人干人人爱 | 国产精品精品国产色婷婷| 亚洲第一青青草原| 美女高潮到喷水免费观看| 99热只有精品国产| 女人精品久久久久毛片| 欧美日韩一级在线毛片| 久久午夜综合久久蜜桃| 又紧又爽又黄一区二区| 日韩欧美三级三区| 久久中文字幕一级| 一卡2卡三卡四卡精品乱码亚洲| 国产成+人综合+亚洲专区| 亚洲熟女毛片儿| 亚洲国产精品999在线| 亚洲第一av免费看| 老司机靠b影院| 久久国产乱子伦精品免费另类| 宅男免费午夜| 美女高潮到喷水免费观看| 国产精品久久久久久人妻精品电影| tocl精华| 久久久精品国产亚洲av高清涩受| √禁漫天堂资源中文www| 成人欧美大片| 精品不卡国产一区二区三区| 波多野结衣av一区二区av| avwww免费| 人人妻人人爽人人添夜夜欢视频| 久久精品影院6| 成在线人永久免费视频| 久久精品91无色码中文字幕| 国产欧美日韩综合在线一区二区| 欧美一级毛片孕妇| 日韩精品中文字幕看吧| 久久精品国产亚洲av香蕉五月| 亚洲精华国产精华精| 女警被强在线播放| 久久久久国产精品人妻aⅴ院| 神马国产精品三级电影在线观看 | 国产91精品成人一区二区三区| 国产精品久久视频播放| 精品久久久久久成人av| av免费在线观看网站| 成人三级黄色视频| 久久国产乱子伦精品免费另类| 成人亚洲精品一区在线观看| 操美女的视频在线观看| 黄片播放在线免费| 男人舔女人的私密视频| 黄色视频,在线免费观看| 欧美中文日本在线观看视频| 岛国在线观看网站| 欧美在线黄色| 88av欧美| 欧美日韩福利视频一区二区| 色播在线永久视频| 亚洲成人精品中文字幕电影| а√天堂www在线а√下载| 久久精品亚洲精品国产色婷小说| 一区福利在线观看| 日韩精品青青久久久久久| 欧美色欧美亚洲另类二区 | 亚洲欧美一区二区三区黑人| 91av网站免费观看| 亚洲欧美一区二区三区黑人| 日韩精品免费视频一区二区三区| 午夜免费观看网址| 国产免费av片在线观看野外av| 一a级毛片在线观看| 久久精品国产亚洲av香蕉五月| 国产一区二区三区综合在线观看| 宅男免费午夜| 在线免费观看的www视频| 久久精品人人爽人人爽视色| 国产高清有码在线观看视频 | 老司机午夜福利在线观看视频| 高潮久久久久久久久久久不卡| 国产精品野战在线观看| 午夜成年电影在线免费观看| 婷婷丁香在线五月| 91在线观看av| 亚洲熟女毛片儿| 深夜精品福利| 亚洲黑人精品在线| 国产人伦9x9x在线观看| 亚洲色图综合在线观看| 高潮久久久久久久久久久不卡| av天堂久久9| 久久久久久久久免费视频了| 久久国产精品人妻蜜桃| 国产成+人综合+亚洲专区| 女同久久另类99精品国产91| 少妇粗大呻吟视频| 亚洲国产毛片av蜜桃av| 男女午夜视频在线观看| 免费av毛片视频| 国产精品免费一区二区三区在线| 90打野战视频偷拍视频| 少妇裸体淫交视频免费看高清 | 婷婷精品国产亚洲av在线| 精品第一国产精品| 禁无遮挡网站| 桃红色精品国产亚洲av| 99久久久亚洲精品蜜臀av| 久久欧美精品欧美久久欧美| 99re在线观看精品视频| 一本久久中文字幕| 日日摸夜夜添夜夜添小说| 天天躁夜夜躁狠狠躁躁| 免费女性裸体啪啪无遮挡网站| 法律面前人人平等表现在哪些方面| 激情视频va一区二区三区| 久久久久久大精品| 中文字幕色久视频| 国产视频一区二区在线看| 亚洲中文av在线| 色av中文字幕| 色哟哟哟哟哟哟| 久久中文字幕人妻熟女| av天堂久久9| 人人妻,人人澡人人爽秒播| 国产精品久久久av美女十八| 成在线人永久免费视频| 国产亚洲欧美精品永久| 老熟妇仑乱视频hdxx| 国产精品香港三级国产av潘金莲| 91国产中文字幕| 女人爽到高潮嗷嗷叫在线视频| 欧美av亚洲av综合av国产av| 亚洲色图 男人天堂 中文字幕| 国产亚洲精品综合一区在线观看 | aaaaa片日本免费| 成年版毛片免费区| 男女之事视频高清在线观看| 熟妇人妻久久中文字幕3abv| 久久精品国产亚洲av高清一级| 亚洲一区二区三区不卡视频| 黄频高清免费视频| 欧美乱色亚洲激情| 国产高清有码在线观看视频 | 欧美日韩亚洲国产一区二区在线观看| 美女免费视频网站| 老司机午夜福利在线观看视频| 国产三级在线视频| 国产区一区二久久| 搡老岳熟女国产| 性欧美人与动物交配| 涩涩av久久男人的天堂| 国产区一区二久久| 一级a爱片免费观看的视频| 免费不卡黄色视频| 国产亚洲精品第一综合不卡| 亚洲欧美激情综合另类| 老汉色∧v一级毛片| 最近最新中文字幕大全免费视频| 亚洲精品av麻豆狂野| 久久精品国产清高在天天线| 男人操女人黄网站| 国产伦一二天堂av在线观看| 亚洲第一青青草原| 精品熟女少妇八av免费久了| 淫妇啪啪啪对白视频| 涩涩av久久男人的天堂| 亚洲七黄色美女视频| 亚洲情色 制服丝袜| 神马国产精品三级电影在线观看 | ponron亚洲| 69av精品久久久久久| 亚洲在线自拍视频| 黄色 视频免费看| 俄罗斯特黄特色一大片| 级片在线观看| 国产精品1区2区在线观看.| 久久久国产精品麻豆| 熟女少妇亚洲综合色aaa.| 久久国产精品影院| 男女床上黄色一级片免费看| x7x7x7水蜜桃| 亚洲专区国产一区二区| 长腿黑丝高跟| 精品日产1卡2卡| 天天躁夜夜躁狠狠躁躁| 久久午夜亚洲精品久久| 18禁观看日本| 国产亚洲欧美98| 欧美一区二区精品小视频在线| 操出白浆在线播放| 人人妻人人澡欧美一区二区 | 91麻豆精品激情在线观看国产| 欧美国产日韩亚洲一区| 美女高潮喷水抽搐中文字幕| 久久 成人 亚洲| 国产精品九九99| 一级毛片高清免费大全| 国产精品自产拍在线观看55亚洲| 成人永久免费在线观看视频| 免费人成视频x8x8入口观看| 在线观看免费视频日本深夜| 久久九九热精品免费| 免费在线观看影片大全网站| 国产极品粉嫩免费观看在线| 亚洲精品在线观看二区| 中出人妻视频一区二区| 国产蜜桃级精品一区二区三区| 久久久久精品国产欧美久久久| 波多野结衣av一区二区av| 久久精品成人免费网站| 无限看片的www在线观看| 女性生殖器流出的白浆| 国产成人精品久久二区二区免费| 久久中文字幕一级| 9热在线视频观看99| 18美女黄网站色大片免费观看| 午夜两性在线视频| 欧美日韩福利视频一区二区| 国产亚洲精品第一综合不卡| 老司机福利观看| 黄色毛片三级朝国网站| 国产一区二区在线av高清观看| 国产激情欧美一区二区| 免费看美女性在线毛片视频| 精品福利观看| 亚洲五月婷婷丁香| 亚洲 国产 在线| 性少妇av在线| 99re在线观看精品视频| 国产aⅴ精品一区二区三区波| www.999成人在线观看| 99精品在免费线老司机午夜| 国内毛片毛片毛片毛片毛片| 国产精品98久久久久久宅男小说| 午夜精品在线福利| 午夜免费成人在线视频| 91麻豆av在线| 他把我摸到了高潮在线观看| 日日夜夜操网爽| 成人亚洲精品一区在线观看| 欧美精品啪啪一区二区三区| 精品国产超薄肉色丝袜足j| 中国美女看黄片| 女同久久另类99精品国产91| 99香蕉大伊视频| 午夜免费成人在线视频| 成年女人毛片免费观看观看9| 国产男靠女视频免费网站| 狂野欧美激情性xxxx| 亚洲三区欧美一区| 青草久久国产| www.自偷自拍.com| 亚洲第一欧美日韩一区二区三区| 国产私拍福利视频在线观看| 日日爽夜夜爽网站| 久久中文字幕人妻熟女| 国产精品免费视频内射| av福利片在线| 老司机靠b影院| 亚洲成av人片免费观看| 亚洲欧美日韩另类电影网站| 777久久人妻少妇嫩草av网站| 大陆偷拍与自拍| 免费观看精品视频网站| 少妇被粗大的猛进出69影院| 国产精品亚洲av一区麻豆| 麻豆久久精品国产亚洲av| 久久精品成人免费网站| 色婷婷久久久亚洲欧美| 国产精品九九99| 91大片在线观看| 国产精品香港三级国产av潘金莲| 无人区码免费观看不卡| 国产乱人伦免费视频| 69av精品久久久久久| 一区二区三区国产精品乱码| 一进一出好大好爽视频| 人妻久久中文字幕网| 国产亚洲精品第一综合不卡| 18禁裸乳无遮挡免费网站照片 | 美女扒开内裤让男人捅视频| 少妇被粗大的猛进出69影院| 91精品国产国语对白视频| 欧美日韩亚洲综合一区二区三区_| 亚洲一区中文字幕在线| 欧美色视频一区免费| 最好的美女福利视频网| 正在播放国产对白刺激| 久久久久久大精品| 午夜激情av网站| 久久久国产成人免费| 欧美日本亚洲视频在线播放| 天堂影院成人在线观看| 精品久久久久久久毛片微露脸| 曰老女人黄片| 99香蕉大伊视频| 97人妻精品一区二区三区麻豆 | 黑人操中国人逼视频| 大型av网站在线播放| 不卡一级毛片| 日本vs欧美在线观看视频| 精品人妻在线不人妻| 国产激情欧美一区二区| 国产成人精品无人区| 国产一级毛片七仙女欲春2 | 久久久久精品国产欧美久久久| 日韩欧美国产在线观看| 免费少妇av软件| 欧美另类亚洲清纯唯美| 日韩精品免费视频一区二区三区| 国产精品 欧美亚洲| 日韩欧美一区视频在线观看| 精品无人区乱码1区二区| 午夜福利高清视频| 亚洲一卡2卡3卡4卡5卡精品中文| 999久久久精品免费观看国产| 国产成人影院久久av| 亚洲欧美精品综合一区二区三区| 69精品国产乱码久久久| 国产精品久久久av美女十八| 亚洲欧美日韩高清在线视频| 51午夜福利影视在线观看| 中文亚洲av片在线观看爽| 在线观看免费午夜福利视频| 男女下面进入的视频免费午夜 | 91老司机精品| 香蕉丝袜av| 亚洲成人免费电影在线观看| 最新在线观看一区二区三区| 亚洲,欧美精品.| 亚洲午夜理论影院| 日韩高清综合在线| 国产三级在线视频| 九色国产91popny在线| 久久中文字幕人妻熟女| 夜夜看夜夜爽夜夜摸| 日韩一卡2卡3卡4卡2021年| 日韩欧美一区视频在线观看| 看黄色毛片网站| 久久亚洲真实| 亚洲av成人不卡在线观看播放网| 黄色丝袜av网址大全| 成人亚洲精品一区在线观看| 亚洲五月色婷婷综合| 黄色片一级片一级黄色片| 两个人免费观看高清视频| 色综合婷婷激情| 在线观看午夜福利视频| 国产亚洲精品第一综合不卡| 欧洲精品卡2卡3卡4卡5卡区| 国产精品免费一区二区三区在线| 啦啦啦观看免费观看视频高清 | 看免费av毛片| 女性生殖器流出的白浆| 一级毛片女人18水好多| 长腿黑丝高跟| 麻豆一二三区av精品| 黑丝袜美女国产一区| 搡老熟女国产l中国老女人| 少妇裸体淫交视频免费看高清 | 久久中文字幕人妻熟女| 亚洲自偷自拍图片 自拍| 满18在线观看网站| 一二三四在线观看免费中文在| 视频在线观看一区二区三区| 国产午夜福利久久久久久| 免费在线观看黄色视频的| av中文乱码字幕在线| 婷婷丁香在线五月| 日韩欧美一区二区三区在线观看| 50天的宝宝边吃奶边哭怎么回事| 天堂√8在线中文| 一级毛片高清免费大全| 精品国产乱码久久久久久男人| 亚洲一区中文字幕在线| 亚洲成av人片免费观看| 午夜福利欧美成人| 日日摸夜夜添夜夜添小说| 久久天躁狠狠躁夜夜2o2o| 久久欧美精品欧美久久欧美| 久久热在线av| 亚洲一码二码三码区别大吗| 自线自在国产av| 久久午夜综合久久蜜桃| 国产亚洲精品久久久久5区| 亚洲成人久久性| 多毛熟女@视频| 俄罗斯特黄特色一大片| 黄色 视频免费看| √禁漫天堂资源中文www| 国产精品精品国产色婷婷| 身体一侧抽搐| √禁漫天堂资源中文www| 极品教师在线免费播放| 老汉色∧v一级毛片| 多毛熟女@视频| 久久中文字幕一级| 成人av一区二区三区在线看| 在线播放国产精品三级| 国产伦人伦偷精品视频| 18禁观看日本| 久久国产精品人妻蜜桃| 大型av网站在线播放| 亚洲伊人色综图| 精品第一国产精品| 男女下面插进去视频免费观看| 亚洲全国av大片| 欧美久久黑人一区二区| 日日干狠狠操夜夜爽| 亚洲专区中文字幕在线| 高清在线国产一区| 久久人妻熟女aⅴ| 自线自在国产av| 国产在线精品亚洲第一网站| 亚洲七黄色美女视频| 波多野结衣高清无吗| 亚洲欧美日韩高清在线视频| a级毛片在线看网站| 亚洲国产毛片av蜜桃av| 成人精品一区二区免费|