• <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型氣動司控道岔的改造
    色综合站精品国产| 亚洲人与动物交配视频| 又爽又黄无遮挡网站| 色尼玛亚洲综合影院| 亚洲国产精品久久男人天堂| 亚洲 国产 在线| 2021天堂中文幕一二区在线观| 亚洲片人在线观看| 精品无人区乱码1区二区| 亚洲人成电影免费在线| 99精品久久久久人妻精品| 青草久久国产| 国产亚洲av嫩草精品影院| 亚洲色图 男人天堂 中文字幕| 别揉我奶头~嗯~啊~动态视频| 欧美中文日本在线观看视频| 一本精品99久久精品77| 亚洲在线观看片| 亚洲成av人片免费观看| 亚洲精品美女久久av网站| 欧美午夜高清在线| 亚洲性夜色夜夜综合| 少妇熟女aⅴ在线视频| 久久午夜亚洲精品久久| 在线永久观看黄色视频| 午夜两性在线视频| 亚洲av电影不卡..在线观看| 狠狠狠狠99中文字幕| 国产一区在线观看成人免费| 久久久久久久久中文| 亚洲欧美激情综合另类| 成人无遮挡网站| 亚洲狠狠婷婷综合久久图片| 欧美成狂野欧美在线观看| 国产91精品成人一区二区三区| 久久香蕉精品热| 久久久久久久久中文| 国产真实乱freesex| 此物有八面人人有两片| 免费观看精品视频网站| 可以在线观看的亚洲视频| 哪里可以看免费的av片| 波多野结衣高清无吗| 真人做人爱边吃奶动态| 天天躁日日操中文字幕| 他把我摸到了高潮在线观看| 国产三级中文精品| av天堂中文字幕网| 91在线观看av| 亚洲一区高清亚洲精品| 69av精品久久久久久| 88av欧美| 人人妻人人澡欧美一区二区| 香蕉久久夜色| 制服人妻中文乱码| 精品无人区乱码1区二区| 欧美性猛交╳xxx乱大交人| 九九久久精品国产亚洲av麻豆 | 成熟少妇高潮喷水视频| 免费看美女性在线毛片视频| 香蕉av资源在线| 日本黄色视频三级网站网址| 窝窝影院91人妻| 欧美黄色淫秽网站| 亚洲五月婷婷丁香| 老司机福利观看| 亚洲九九香蕉| 亚洲乱码一区二区免费版| 亚洲午夜精品一区,二区,三区| 又黄又粗又硬又大视频| 国产精品亚洲美女久久久| 色吧在线观看| 色老头精品视频在线观看| 欧美另类亚洲清纯唯美| 日日夜夜操网爽| 午夜福利在线观看吧| 91av网一区二区| 国产一区二区三区视频了| 亚洲最大成人中文| 啦啦啦韩国在线观看视频| 欧美极品一区二区三区四区| 三级毛片av免费| 国产亚洲精品久久久com| 嫩草影院精品99| 国产av不卡久久| 男女做爰动态图高潮gif福利片| 久久国产乱子伦精品免费另类| 国内少妇人妻偷人精品xxx网站 | 亚洲激情在线av| 99国产极品粉嫩在线观看| 亚洲人成网站在线播放欧美日韩| 99国产精品一区二区三区| xxx96com| 一卡2卡三卡四卡精品乱码亚洲| 久久这里只有精品19| 久久天躁狠狠躁夜夜2o2o| 视频区欧美日本亚洲| 国产精品,欧美在线| 一级a爱片免费观看的视频| 999精品在线视频| 少妇熟女aⅴ在线视频| 亚洲国产日韩欧美精品在线观看 | 啪啪无遮挡十八禁网站| 国产一区二区激情短视频| 久久国产乱子伦精品免费另类| 久久久久亚洲av毛片大全| av天堂中文字幕网| 国产精品一区二区三区四区免费观看 | 午夜精品在线福利| 一个人免费在线观看电影 | 欧美另类亚洲清纯唯美| 夜夜夜夜夜久久久久| 麻豆国产av国片精品| 亚洲五月婷婷丁香| 两性夫妻黄色片| 后天国语完整版免费观看| 男人和女人高潮做爰伦理| 男插女下体视频免费在线播放| 高清毛片免费观看视频网站| 欧美成人一区二区免费高清观看 | 看片在线看免费视频| 啦啦啦观看免费观看视频高清| 男女那种视频在线观看| 亚洲欧洲精品一区二区精品久久久| 一卡2卡三卡四卡精品乱码亚洲| 日本三级黄在线观看| 制服人妻中文乱码| 99热精品在线国产| 免费一级毛片在线播放高清视频| 看免费av毛片| 久久精品国产综合久久久| 亚洲国产日韩欧美精品在线观看 | 国产亚洲精品综合一区在线观看| 91在线观看av| 日本熟妇午夜| 又黄又爽又免费观看的视频| 又紧又爽又黄一区二区| 亚洲成人精品中文字幕电影| 亚洲片人在线观看| 国产一区在线观看成人免费| 久久精品国产亚洲av香蕉五月| 国产欧美日韩一区二区三| 丰满的人妻完整版| 宅男免费午夜| 亚洲精品国产精品久久久不卡| av在线天堂中文字幕| 中国美女看黄片| 男插女下体视频免费在线播放| 最近在线观看免费完整版| 久久这里只有精品19| 精品午夜福利视频在线观看一区| 色尼玛亚洲综合影院| 真实男女啪啪啪动态图| 亚洲在线观看片| 少妇的逼水好多| 亚洲va日本ⅴa欧美va伊人久久| 日韩精品中文字幕看吧| 在线a可以看的网站| 12—13女人毛片做爰片一| 国产高潮美女av| 99国产极品粉嫩在线观看| 国产又色又爽无遮挡免费看| 成人三级黄色视频| 久久精品夜夜夜夜夜久久蜜豆| 他把我摸到了高潮在线观看| 国产高清视频在线观看网站| 女生性感内裤真人,穿戴方法视频| 俄罗斯特黄特色一大片| 狠狠狠狠99中文字幕| 99riav亚洲国产免费| 精品欧美国产一区二区三| 国产精品美女特级片免费视频播放器 | 最新在线观看一区二区三区| 亚洲av电影不卡..在线观看| 国产人伦9x9x在线观看| 亚洲人成网站高清观看| 亚洲av成人一区二区三| 亚洲国产精品合色在线| 国产乱人伦免费视频| 亚洲av五月六月丁香网| 亚洲精品在线美女| 2021天堂中文幕一二区在线观| 国产一区在线观看成人免费| 午夜免费观看网址| 国产精品爽爽va在线观看网站| 亚洲狠狠婷婷综合久久图片| 天天躁日日操中文字幕| 欧美成狂野欧美在线观看| 很黄的视频免费| 精品99又大又爽又粗少妇毛片 | 欧洲精品卡2卡3卡4卡5卡区| 精品久久久久久久久久久久久| 国产av麻豆久久久久久久| 亚洲无线观看免费| 99国产精品一区二区三区| 免费av不卡在线播放| 久久香蕉精品热| 久久久久久九九精品二区国产| 人人妻人人澡欧美一区二区| 亚洲成人中文字幕在线播放| 亚洲,欧美精品.| 亚洲国产欧美人成| av在线天堂中文字幕| 一个人观看的视频www高清免费观看 | 亚洲av成人不卡在线观看播放网| 精品电影一区二区在线| 精华霜和精华液先用哪个| 在线观看66精品国产| 亚洲精品一卡2卡三卡4卡5卡| 淫妇啪啪啪对白视频| netflix在线观看网站| 又黄又爽又免费观看的视频| 黑人巨大精品欧美一区二区mp4| 12—13女人毛片做爰片一| 亚洲精品一区av在线观看| 欧美av亚洲av综合av国产av| 美女黄网站色视频| 亚洲无线在线观看| 一本久久中文字幕| 女人被狂操c到高潮| 在线国产一区二区在线| 久久香蕉精品热| 日韩欧美一区二区三区在线观看| 香蕉丝袜av| 亚洲欧美激情综合另类| 日韩大尺度精品在线看网址| 婷婷精品国产亚洲av| 久久国产乱子伦精品免费另类| 麻豆国产97在线/欧美| 十八禁人妻一区二区| 色综合欧美亚洲国产小说| 熟妇人妻久久中文字幕3abv| 此物有八面人人有两片| 变态另类丝袜制服| 精品久久久久久久末码| 精品午夜福利视频在线观看一区| 一区二区三区国产精品乱码| 久久精品91蜜桃| 精品国产美女av久久久久小说| 免费看美女性在线毛片视频| 日韩有码中文字幕| 久久久久久久久中文| 亚洲一区二区三区色噜噜| 97超视频在线观看视频| 日本免费a在线| 狂野欧美白嫩少妇大欣赏| 久久午夜综合久久蜜桃| 成人高潮视频无遮挡免费网站| 亚洲精品在线观看二区| 国产午夜精品久久久久久| 国产淫片久久久久久久久 | 国产激情偷乱视频一区二区| 老司机在亚洲福利影院| 久久亚洲精品不卡| 母亲3免费完整高清在线观看| www国产在线视频色| 欧美高清成人免费视频www| 999精品在线视频| 欧美午夜高清在线| 丰满人妻一区二区三区视频av | 男人的好看免费观看在线视频| 丝袜人妻中文字幕| 高清毛片免费观看视频网站| 精品久久久久久久末码| 一边摸一边抽搐一进一小说| 两个人视频免费观看高清| 狠狠狠狠99中文字幕| 老鸭窝网址在线观看| 欧美在线一区亚洲| 99riav亚洲国产免费| 又大又爽又粗| 亚洲国产精品成人综合色| 99国产精品一区二区蜜桃av| 色播亚洲综合网| 国产精品 国内视频| a在线观看视频网站| 日本 欧美在线| 国产成人精品久久二区二区91| 两个人视频免费观看高清| 91在线精品国自产拍蜜月 | 成年女人看的毛片在线观看| 精品国产亚洲在线| 他把我摸到了高潮在线观看| 欧美在线一区亚洲| 99久久精品国产亚洲精品| 亚洲av中文字字幕乱码综合| av在线蜜桃| 国内揄拍国产精品人妻在线| 亚洲欧美日韩卡通动漫| 日本撒尿小便嘘嘘汇集6| 怎么达到女性高潮| 99久久久亚洲精品蜜臀av| 国产精品一区二区三区四区久久| 69av精品久久久久久| 又爽又黄无遮挡网站| 99国产精品一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品自产拍在线观看55亚洲| 国产真实乱freesex| 91老司机精品| 日韩中文字幕欧美一区二区| 色吧在线观看| 国产亚洲av嫩草精品影院| xxx96com| 九九在线视频观看精品| 国产私拍福利视频在线观看| 午夜福利在线观看吧| 一卡2卡三卡四卡精品乱码亚洲| 国产精品国产高清国产av| 国产精品久久久av美女十八| 国产一区二区三区在线臀色熟女| 免费在线观看亚洲国产| 少妇的逼水好多| 人人妻,人人澡人人爽秒播| 日日干狠狠操夜夜爽| 法律面前人人平等表现在哪些方面| 国产精品亚洲一级av第二区| 亚洲av美国av| 欧美色视频一区免费| 精品一区二区三区av网在线观看| 国产精品综合久久久久久久免费| 国产欧美日韩精品亚洲av| 男人和女人高潮做爰伦理| 波多野结衣巨乳人妻| 亚洲精品在线观看二区| 国内毛片毛片毛片毛片毛片| 欧美中文日本在线观看视频| 日韩国内少妇激情av| 欧美在线黄色| 99国产精品一区二区蜜桃av| 精品久久久久久久久久久久久| 一区二区三区激情视频| 国产成年人精品一区二区| 久9热在线精品视频| 国产欧美日韩精品一区二区| 最好的美女福利视频网| 精品免费久久久久久久清纯| 人妻久久中文字幕网| 国产精品一区二区三区四区免费观看 | 黄频高清免费视频| 欧美中文日本在线观看视频| 亚洲av五月六月丁香网| 婷婷六月久久综合丁香| 精品久久久久久久末码| svipshipincom国产片| 手机成人av网站| 在线免费观看的www视频| 99在线视频只有这里精品首页| 精品久久久久久久人妻蜜臀av| 人人妻人人澡欧美一区二区| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产中文字幕在线视频| 午夜两性在线视频| 日本精品一区二区三区蜜桃| 午夜精品久久久久久毛片777| 久久精品国产综合久久久| 欧美在线黄色| 精品熟女少妇八av免费久了| 亚洲无线观看免费| 国产成人系列免费观看| 91久久精品国产一区二区成人 | a级毛片在线看网站| 日本一本二区三区精品| 精品电影一区二区在线| 香蕉丝袜av| 丰满人妻熟妇乱又伦精品不卡| 成年女人看的毛片在线观看| 波多野结衣高清作品| e午夜精品久久久久久久| 99国产精品一区二区三区| 亚洲精品在线美女| 久久精品aⅴ一区二区三区四区| 99久久无色码亚洲精品果冻| 久久久精品大字幕| 丁香欧美五月| 国产精品久久视频播放| 日韩大尺度精品在线看网址| 此物有八面人人有两片| 免费搜索国产男女视频| 一级毛片女人18水好多| 91av网一区二区| 国产精华一区二区三区| 人人妻人人看人人澡| 婷婷六月久久综合丁香| 欧美日韩一级在线毛片| 在线观看舔阴道视频| 国产真人三级小视频在线观看| 成年人黄色毛片网站| 国产成+人综合+亚洲专区| 色综合亚洲欧美另类图片| 久9热在线精品视频| 欧美绝顶高潮抽搐喷水| 18禁黄网站禁片免费观看直播| 黄片大片在线免费观看| 99久久精品一区二区三区| 啦啦啦免费观看视频1| 后天国语完整版免费观看| 在线国产一区二区在线| 黄色女人牲交| 精品无人区乱码1区二区| 欧美日韩精品网址| 久久久久久久午夜电影| 美女扒开内裤让男人捅视频| 国产99白浆流出| 亚洲av第一区精品v没综合| 久久人妻av系列| 巨乳人妻的诱惑在线观看| av女优亚洲男人天堂 | 国产一区二区三区在线臀色熟女| 国产精品久久电影中文字幕| 欧美三级亚洲精品| 制服丝袜大香蕉在线| 丰满人妻熟妇乱又伦精品不卡| www.自偷自拍.com| 午夜福利在线观看吧| 亚洲国产日韩欧美精品在线观看 | 亚洲精品色激情综合| av片东京热男人的天堂| 在线观看免费视频日本深夜| 国产精品爽爽va在线观看网站| 国产精品永久免费网站| 国产三级中文精品| 欧美成人一区二区免费高清观看 | 国产一区二区在线av高清观看| 老司机福利观看| 国产一级毛片七仙女欲春2| 极品教师在线免费播放| 国产高清三级在线| 麻豆久久精品国产亚洲av| 两性午夜刺激爽爽歪歪视频在线观看| 精品国产美女av久久久久小说| 国产真实乱freesex| 国产精品99久久99久久久不卡| 最近视频中文字幕2019在线8| 日本在线视频免费播放| 午夜福利高清视频| 搡老妇女老女人老熟妇| 免费观看精品视频网站| av欧美777| 国内精品久久久久久久电影| av福利片在线观看| 亚洲av日韩精品久久久久久密| 午夜久久久久精精品| 日本五十路高清| 两性夫妻黄色片| 曰老女人黄片| 免费高清视频大片| 精品日产1卡2卡| 男人舔女人的私密视频| 亚洲国产欧美一区二区综合| 脱女人内裤的视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲人与动物交配视频| 欧美日韩综合久久久久久 | 国产精品一区二区三区四区久久| 精品久久蜜臀av无| 久久久久久九九精品二区国产| 成人无遮挡网站| 久久精品亚洲精品国产色婷小说| 亚洲av成人av| 亚洲国产中文字幕在线视频| 午夜福利在线在线| 免费人成视频x8x8入口观看| 窝窝影院91人妻| 黄片小视频在线播放| 午夜福利成人在线免费观看| 嫩草影视91久久| 国产成人影院久久av| 免费看光身美女| 亚洲精品粉嫩美女一区| 伦理电影免费视频| 亚洲国产精品999在线| 在线观看日韩欧美| 国产精品98久久久久久宅男小说| 久久中文字幕人妻熟女| 成人一区二区视频在线观看| 国产免费av片在线观看野外av| 又爽又黄无遮挡网站| 男人舔女人的私密视频| 美女cb高潮喷水在线观看 | 少妇裸体淫交视频免费看高清| 一级毛片高清免费大全| 国产av不卡久久| 99久久成人亚洲精品观看| 国内揄拍国产精品人妻在线| 美女扒开内裤让男人捅视频| 一级a爱片免费观看的视频| 精品无人区乱码1区二区| 国产亚洲av高清不卡| 岛国在线免费视频观看| 国产亚洲精品综合一区在线观看| 男人舔女人下体高潮全视频| 亚洲五月婷婷丁香| 日日摸夜夜添夜夜添小说| 午夜免费成人在线视频| 亚洲国产欧美人成| 欧美日韩一级在线毛片| 亚洲真实伦在线观看| 热99在线观看视频| 一个人观看的视频www高清免费观看 | 宅男免费午夜| 最新在线观看一区二区三区| 一进一出抽搐gif免费好疼| 给我免费播放毛片高清在线观看| 久久国产精品影院| 免费人成视频x8x8入口观看| 香蕉久久夜色| 视频区欧美日本亚洲| 怎么达到女性高潮| 国产高清激情床上av| 午夜福利成人在线免费观看| 精品久久久久久,| 日本黄色片子视频| 日本五十路高清| 狂野欧美白嫩少妇大欣赏| 亚洲成人久久性| 性色avwww在线观看| 黄色成人免费大全| 欧美在线一区亚洲| 在线十欧美十亚洲十日本专区| 露出奶头的视频| 麻豆一二三区av精品| 亚洲激情在线av| 日本三级黄在线观看| 国产精品香港三级国产av潘金莲| 99久久精品国产亚洲精品| 日韩欧美国产一区二区入口| 一进一出抽搐动态| 不卡一级毛片| 国产亚洲av嫩草精品影院| 久久国产精品影院| www.www免费av| 免费无遮挡裸体视频| 日本 欧美在线| 国产黄片美女视频| 亚洲黑人精品在线| 在线十欧美十亚洲十日本专区| 我要搜黄色片| 欧美乱码精品一区二区三区| 俄罗斯特黄特色一大片| 欧美乱码精品一区二区三区| 久久久精品欧美日韩精品| 免费一级毛片在线播放高清视频| 欧美黑人欧美精品刺激| 久久人妻av系列| 久久久久久久久久黄片| 亚洲在线自拍视频| 亚洲熟妇中文字幕五十中出| 亚洲人成网站高清观看| 中文字幕人妻丝袜一区二区| 午夜福利成人在线免费观看| 99re在线观看精品视频| 熟女少妇亚洲综合色aaa.| av国产免费在线观看| 亚洲欧美精品综合久久99| 91av网一区二区| 精品午夜福利视频在线观看一区| 美女免费视频网站| 成年女人毛片免费观看观看9| 国产野战对白在线观看| 夜夜夜夜夜久久久久| 欧美又色又爽又黄视频| 在线观看免费午夜福利视频| 老汉色∧v一级毛片| 亚洲色图 男人天堂 中文字幕| 国产午夜福利久久久久久| 美女cb高潮喷水在线观看 | 日本免费a在线| 超碰成人久久| 亚洲欧美日韩卡通动漫| 欧美日韩瑟瑟在线播放| 亚洲国产精品成人综合色| 老司机在亚洲福利影院| 亚洲精品久久国产高清桃花| 国产亚洲精品久久久久久毛片| 中国美女看黄片| 全区人妻精品视频| 国产亚洲av嫩草精品影院| 亚洲中文字幕一区二区三区有码在线看 | 欧洲精品卡2卡3卡4卡5卡区| 在线永久观看黄色视频| 在线播放国产精品三级| 99热只有精品国产| 可以在线观看毛片的网站| 亚洲狠狠婷婷综合久久图片| 免费在线观看视频国产中文字幕亚洲| 日韩av在线大香蕉| 精品不卡国产一区二区三区| 国产亚洲av嫩草精品影院| 色哟哟哟哟哟哟| 成人性生交大片免费视频hd| 久久久久久九九精品二区国产| 免费在线观看亚洲国产| 村上凉子中文字幕在线| 精品午夜福利视频在线观看一区| 国产精品香港三级国产av潘金莲| 天堂影院成人在线观看| 老汉色av国产亚洲站长工具| 色综合婷婷激情| 制服丝袜大香蕉在线| 婷婷亚洲欧美| 18禁裸乳无遮挡免费网站照片| 国产极品精品免费视频能看的| 后天国语完整版免费观看| 日本黄色视频三级网站网址| 精品熟女少妇八av免费久了| 免费看a级黄色片| 老司机福利观看| 欧美午夜高清在线| 中文字幕久久专区| 天天添夜夜摸| 日韩精品中文字幕看吧| 嫁个100分男人电影在线观看| 12—13女人毛片做爰片一| 日本撒尿小便嘘嘘汇集6|