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

    復雜約束條件下再入高超聲速滑翔飛行器軌跡快速優(yōu)化①

    2017-01-05 09:40:20邢清華毛藝帆
    固體火箭技術 2016年6期
    關鍵詞:偽譜禁飛區(qū)約束條件

    王 路,邢清華,毛藝帆

    (空軍工程大學 防空反導學院,西安 710051)

    復雜約束條件下再入高超聲速滑翔飛行器軌跡快速優(yōu)化①

    王 路,邢清華,毛藝帆

    (空軍工程大學 防空反導學院,西安 710051)

    作為直接法的典型代表,高斯偽譜方法在處理復雜約束條件(含路徑點或禁飛區(qū)約束)下再入高超聲速滑翔飛行器軌跡優(yōu)化問題時,僅能保證所得軌跡在各高斯節(jié)點處嚴格滿足各項約束,而節(jié)點之間軌跡的可行性無法保證,為解決這一問題,文章提出改進多階段高斯偽譜方法。該方法首先使用新定義的兩類節(jié)點(固定節(jié)點和虛擬節(jié)點)將軌跡分段,其中固定節(jié)點是用來保證路徑點與高斯節(jié)點重合,虛擬節(jié)點則是用來保證禁飛區(qū)附近分布更多的高斯節(jié)點,此分段方式能保證所得軌跡在任意位置可行;然后,向各分段軌跡插入指定數(shù)目的高斯節(jié)點;最后,使用序列二次規(guī)劃方法對此多階段軌跡優(yōu)化模型進行求解。仿真結果表明,改進多階段高斯偽譜方法在處理復雜約束條件下再入高超聲速滑翔飛行器軌跡優(yōu)化問題時快速有效。

    再入高超聲速滑翔飛行器;軌跡優(yōu)化;路徑點;禁飛區(qū);高斯偽譜方法

    0 引言

    高超聲速滑翔飛行器因其高升阻比外形設計具有良好的氣動性能,能憑借空氣動力的作用,實現(xiàn)遠距離無動力滑翔飛行。由于飛行器運動方程形式復雜、再入滑翔軌跡對控制量敏感,且飛行過程中受過載、熱流密度等約束影響較大,其軌跡優(yōu)化問題一直是研究的熱點。

    通常軌跡優(yōu)化的數(shù)值解法從本質(zhì)上可分為直接法和間接法[1]。間接法根據(jù)一階最優(yōu)必要條件和龐特里亞金最小化原理,將最優(yōu)化問題轉化為兩點邊值問題[2],缺點是對初值高度敏感,很難收斂,時效性不高。而直接法將原始最優(yōu)控制問題轉化為非線性規(guī)劃問題,在求解過程中具有較大的收斂半徑,且不受一階必要條件的限制[3]。

    作為直接法中的典型代表,高斯偽譜方法使用少量的節(jié)點,并具有較快的收斂速度,目前已被廣泛應用于普通約束條件下高超聲速飛行器上升段[4]或再入段[5-11]的軌跡優(yōu)化問題,這里的普通約束主要包括動壓、過載、熱流密度等過程約束、控制約束及終端狀態(tài)約束。然而,在飛行器實際再入過程中,不得不考慮某些真實約束,如路徑點和禁飛區(qū)約束。路徑點是指偵察或載荷投送時,飛行器必須經(jīng)過的空間位置,禁飛區(qū)是指受地理政治或敵方武器系統(tǒng)威脅,飛行器應嚴格避讓的某些特殊區(qū)域。

    高斯偽譜法將狀態(tài)變量和控制變量在一系列高斯節(jié)點上離散,所有約束條件在高斯節(jié)點上嚴格滿足,但高斯節(jié)點的分布呈中間疏松、兩端密集的特點。此時,考慮路徑點和禁飛區(qū)約束條件,可能會帶來以下問題:

    (1)相鄰高斯節(jié)點之間的軌跡不一定滿足禁飛區(qū)約束或過程約束;

    (2)路徑點約束條件不一定能夠滿足高精度。

    針對問題1和問題2,Jiang[12]提出了基于高斯偽譜方法的多階段軌跡優(yōu)化技術,通過所提出的3類插入點,將整個再入軌跡劃分為不同階段,由于其所設計的3類插入點都是固定的,所以最終得到的軌跡必須嚴格經(jīng)過各插入點。但通過研究發(fā)現(xiàn),在軌跡生成之前,只有路徑點是確定的,而禁飛區(qū)附近的飛行軌跡是未知的,所以只有路徑點可提前設置為固定插入點,Jorris[13-14]和Hu[15]所研究的軌跡優(yōu)化問題考慮了路徑點和禁飛區(qū)等復雜約束條件,但沒有具體討論高斯偽譜方法存在的上述兩類問題。針對以上問題,還有一類解決方法,就是使用大量的高斯節(jié)點對問題進行求解,但代價是十分耗時,實用價值不強。

    為克服傳統(tǒng)高斯偽譜方法在處理復雜約束條件下軌跡優(yōu)化問題時存在的不足,本文提出改進多階段高斯偽譜方法,通過使用新定義的兩類節(jié)點(固定節(jié)點和虛擬節(jié)點)將軌跡分段,在不增加總節(jié)點數(shù)目的同時,改變節(jié)點的分布。一方面,使路徑點與高斯節(jié)點重合,從而保證路徑點約束被嚴格滿足;另一方面,是禁飛區(qū)附近分布高密度的高斯節(jié)點,從而保證禁飛區(qū)約束滿足要求,又由于節(jié)點數(shù)目并沒有增加。所以,可保證所提算法的快速性。仿真結果驗證了其有效性。

    1 問題描述

    本節(jié)主要介紹高超聲速飛行器再入軌跡優(yōu)化問題模型,包括再入動態(tài)方程、約束條件和目標函數(shù)。

    1.1 再入動態(tài)方程

    在高超聲速飛行器的再入階段,假設目標為質(zhì)點,并將地球視為非旋轉勻質(zhì)球體,則飛行器再入運動方程如下[16]:

    (1)

    式中V為飛行器速度,m/s;θ為航跡傾角,rad;ψ為航跡偏角,rad;r為飛行器質(zhì)心至地心的距離,m;λ為地理經(jīng)度,rad;φ為地理緯度,rad;γc為滾轉角,rad;m為飛行器質(zhì)量,kg;D為阻力,D=CDρSV2/2;L為升力,L=CLρSV2/2;CD為阻力系數(shù);CL為升力系數(shù);S為氣動參考面積。

    在高超聲速條件下,升力、阻力系數(shù)CL、CD可近似表示為攻角α(°)與速度V(Ma)的函數(shù)[17]。

    (2)

    式中CLi和CDi(i=0,…,3)為常系數(shù)。

    1.2 約束條件

    1.2.1 過程約束

    高超聲速飛行器再入滑翔階段是一個高動態(tài)的過程,為滿足系統(tǒng)熱防護和結構穩(wěn)定性的需要,再入過程需滿足一系列苛刻約束條件的限制,如動壓限制、熱流密度限制、過載限制。

    (3)

    1.2.2 控制約束

    為了保證系統(tǒng)的穩(wěn)定性與可控性,攻角與滾轉角應滿足如下約束:

    (4)

    式中 下標“min”與“max”為對應變量的下限與上限。

    1.2.3 終端約束

    為使飛行器具備一定的末段打擊能力,對終端約束條件也應進行一定的限制。

    (5)

    式中 下標“f”為終端狀態(tài)。

    1.2.4 路徑點與禁飛區(qū)約束

    路徑點約束:為了完成諸如偵察或載荷投擲等特定任務,飛行器需精準的飛過指定路徑點上空,對應路徑點約束表示如下:

    (6)

    式中Npoint為路徑點的總個數(shù)。

    禁飛區(qū)約束:受地理政治或敵方武器系統(tǒng)威脅的限制,飛行器應嚴格避讓某些特殊區(qū)域,在此假設禁飛區(qū)為不限高度的圓柱體[13],第j個禁飛區(qū)的中心點坐標用(λj,φj)表示,對應的禁飛區(qū)半徑用Rj表示,則禁飛區(qū)約束表示為

    (7)

    式中Nzone為禁飛區(qū)的個數(shù)。

    1.3 目標函數(shù)

    通常情況下,不同的任務對應不同的目標函數(shù)。如果是為了維持飛行器較好的熱防護,目標函數(shù)可設置為最小熱流密度或最小峰值熱流密度;如果打擊或終點目標是固定的,目標函數(shù)可設置為達到時間最短;而如果是為了測試飛行器的性能,目標函數(shù)可設置為橫向或縱向航程最大等等。

    1.4 軌跡優(yōu)化

    對于高超聲速飛行器,其再入軌跡優(yōu)化問題就是尋找控制量α和γ,在滿足各項約束條件的同時,使得目標函數(shù)最優(yōu)。

    2 改進多階段高斯偽譜方法

    2.1 高斯偽譜法及其分析

    2.1.1 高斯偽譜法

    可用一般最優(yōu)控制問題描述高超聲速飛行器再入軌跡優(yōu)化問題,這里采用Bolza問題形式[18]進行表述。連續(xù)Bolza問題可描述為:求解控制變量u(t)∈Rm和狀態(tài)變量x(t)∈Rn,使如下目標函數(shù)[7]最小化。

    (8)

    同時滿足如下運動方程、過程約束以及終端約束:

    (9)

    式(8)、式(9)被稱為連續(xù)Bolza 問題的數(shù)學描述。高斯偽譜法在高斯節(jié)點上將狀態(tài)變量和控制變量進行離散,并利用這些離散點構造全局拉格朗日插值多項式來近似系統(tǒng)的動力學方程,從而將連續(xù)的無限維最優(yōu)控制問題轉換為非線性規(guī)劃問題進行求解[19]。對于高超聲速飛行器,飛行持續(xù)時間為從t0到tf,而高斯節(jié)點的分布區(qū)間為[-1,1]。因此,需將時間區(qū)域從t∈[t0,tf]轉化至τ∈[-1,1],轉化公式如式(10):

    (10)

    需要注意的是,N個離散的節(jié)點包含了初始節(jié)點τ0=-1、終端節(jié)點τf=1以及K=N-2個內(nèi)部節(jié)點,且這K個節(jié)點為K階拉格朗日插值多項式的根,則狀態(tài)和控制變量可近似表示如下:

    (11)

    其中,Li(τ)和Li*(τ)分別為N+1和N階拉格朗日插值多項式的基函數(shù),計算公式如下:

    (12)

    為獲得動態(tài)運動方程的約束,對式(11)求導可得

    (13)

    高斯偽譜方法中每個節(jié)點處目標狀態(tài)變量的倒數(shù)都應滿足式(13)。微分矩陣D∈RN×N+1可離線生成,其各元素計算公式如下[20]:

    式中PN(τ)為N階拉格朗日插值多項式;i=1,2,…,N;k=1,2,…,N。

    然后,將式(9)所述動態(tài)方程約束轉化為如下等式約束:

    其中,k=1,2,…,N。由于式(12)只表示高斯節(jié)點的目標狀態(tài),不包括初始時刻和終端時刻的目標狀態(tài),因此應增加2個端點處的等式約束如下[21]:

    (16)

    式中ωk為高斯積分的權重。

    然后,目標函數(shù)中的積分項通過高斯積分近似,則目標函數(shù)可表示為

    (17)

    同時,可將終端約束和過程約束表示為

    至此,將最優(yōu)控制問題轉化為非線性規(guī)劃問題,描述如下:求解非線性規(guī)劃問題的變量:

    (19)

    在滿足式(15)、式(16)和(式18)所述約束條件的同時,使得式(17)所述目標函數(shù)達到最小。

    2.1.2 分析

    如圖1所示,高斯節(jié)點的分布呈中間疏松兩端密集的特點。

    圖1 拉格朗日高斯節(jié)點分布示意圖Fig.1 Distribution of LG nodes

    高斯節(jié)點的這種分布特性會帶來以下幾方面的問題:

    (1)當高斯節(jié)點數(shù)目確定時,節(jié)點的分布位置在區(qū)間[-1,1]是確定的,如圖2所示,由于約束條件只在節(jié)點處嚴格滿足,所以很難保證所生成軌跡滿足路徑點約束,除非路徑點和某一高斯節(jié)點恰好重合。

    (2)如圖3所示,盡管禁飛區(qū)約束在所有高斯節(jié)點均嚴格滿足,但高斯節(jié)點之間的飛行器軌跡穿過了禁飛區(qū)。

    (3)過程約束條件不能保證嚴格滿足,如圖4所示。盡管某項過程約束指標在所有高斯節(jié)點均嚴格滿足,但該指標在多處超出了允許范圍,這會為飛行器維持其結構的穩(wěn)定性帶來嚴重挑戰(zhàn)。

    圖2 路徑點附近軌跡示意圖Fig.2 Example of trajectory near the waypoints

    圖3 禁飛區(qū)附近軌跡示意圖Fig.3 Example of trajectory near the no-fly zone

    圖4 過程約束變化示意圖Fig.4 Change flow of path constraint

    2.2 改進多階段高斯偽譜算法

    為克服前文所述傳統(tǒng)高斯偽譜方法在處理復雜約束條件下的軌跡優(yōu)化問題時存在的問題,本節(jié)提出改進多階段高斯偽譜算法。

    2.2.1 相關定義

    定義1 固定節(jié)點:只包含位置信息,不包含時間信息的節(jié)點,如路徑點。

    定義2 虛擬節(jié)點:只含時間信息,不包含位置信息的節(jié)點。

    2.2.2 算法流程

    (1)步驟1 使用固定節(jié)點將軌跡劃分為多階段

    如圖5所示,將路徑點設置為固定節(jié)點,使用固定節(jié)點將軌跡劃分為多階段,如此以來,路徑點約束可嚴格滿足。轉步驟2。

    圖5 軌跡分段示意圖(使用固定節(jié)點)Fig.5 Trajectory divided by fixed nodes

    (2)步驟2 模型求解

    假設軌跡被劃分為Z段,每階段的高斯節(jié)點數(shù)目為N,則第z階段的目標函數(shù)可表示為

    (20)

    則總目標函數(shù)可表示為

    (21)

    對應的終端和過程約束表示為

    (22)

    式(15)轉變?yōu)?/p>

    (23)

    式(16)轉變?yōu)?/p>

    (24)

    為保證各階段連接點處的連續(xù)性,增加如下約束:

    (25)

    至此,式(20)~式(25)將傳統(tǒng)高斯偽譜法轉化為多階段高斯偽譜法。模型的處理方式同傳統(tǒng)高斯偽譜法一致,本文選用適用于求解大規(guī)模非線性規(guī)劃問題的序列二次規(guī)劃方法。

    模型求解后,得到當前最優(yōu)軌跡,轉步驟3。

    (3)步驟3 軌跡可行性判別

    (26)

    經(jīng)檢測如果各約束條件在所有檢測點處嚴格滿足,則停止,所得軌跡即為最優(yōu)軌跡;否則,轉步驟4。

    (4)步驟4 增加虛擬節(jié)點將軌跡劃分為多階段

    步驟3中不滿足約束條件的檢測點由于只包含時間信息,因此將其視為虛擬節(jié)點,如圖7所示。軌跡被固定節(jié)點和新增虛擬節(jié)點劃分為多階段,由于高斯節(jié)點的分布特性,在虛擬節(jié)點的周圍會有較多的高斯節(jié)點分布,通過這種方式,可保證過程約束或路徑約束得到滿足。然后,轉步驟2。

    圖6 高斯節(jié)點和檢測點分布示意圖Fig.6 Distribution of LG nodes and checkpoints

    圖7 虛擬節(jié)點插入后禁飛區(qū)約束得到滿足示意圖Fig.7 No-fly zone constraint being satisfied by inserting virtual node

    2.2.3 初值生成算法

    本文使用傳統(tǒng)的初始化策略:先采用較少的高斯節(jié)點,求解出一條低精度的軌跡;然后,通過對此低精度軌跡進行插值,得到理想數(shù)目高斯節(jié)點所對應的各節(jié)點初始值。

    2.2.4 改進多階段高斯偽譜法分析

    改進多階段高斯偽譜法主要思想為:首先將路徑點視為固定節(jié)點,并用固定節(jié)點將軌跡分段,進行求解,從而使得到的軌跡嚴格滿足路徑點約束條件;為使全局軌跡嚴格滿足禁飛區(qū)和路經(jīng)約束,需要對所生成軌跡進行可行性驗證,若不滿足要求,則需增加虛擬節(jié)點,將軌跡重新分段,并進行求解,直至所生成軌跡嚴格滿足所有約束條件后即可。

    相比傳統(tǒng)高斯偽譜方法,改進多階段高斯偽譜方法在處理復雜約束條件下的軌跡優(yōu)化問題時,具有如下優(yōu)勢:

    (1)通過使用固定節(jié)點,將軌跡進行分段,所得軌跡嚴格滿足路徑點約束條件。

    (2)通過使用虛擬節(jié)點,將軌跡進行分段,所得軌跡不僅在高斯節(jié)點處嚴格滿足禁飛區(qū)以及過程約束,高斯節(jié)點之間的可行性同樣可保證。

    (3)固定節(jié)點和虛擬節(jié)點的加入,沒有改變高斯節(jié)點的數(shù)目,而是改變了高斯節(jié)點的分布。所以,算法的速度沒有受到影響。

    3 仿真及分析

    表1 初始時刻及終端時刻飛行器狀態(tài)Table 1 Initial and terminal condition of vehicle

    表2 路徑點與禁飛區(qū)參數(shù)設置Table 2 Parameters of waypoint and no-fly zone constraints

    3.1 粒子群算法仿真結果

    粒子的數(shù)目為20,其他參數(shù)設置與文獻[22]一致,所得最優(yōu)軌跡如圖8所示。

    (a) 三維軌跡 (b) 俯視軌跡 (c) 禁飛區(qū)1附近軌跡

    (d) 禁飛區(qū)2附近軌跡 (e) 路徑點附近軌跡圖8 粒子群算法仿真結果Fig.8 Results of optimal trajectory by PSO

    3.2 傳統(tǒng)高斯偽譜方法仿真結果

    高斯偽譜方法選用的高斯節(jié)點數(shù)目為30,仿真結果如圖9所示。

    (a) 三維軌跡

    (b) 俯視軌跡

    (d) 禁飛區(qū)2附近軌跡圖9 高斯偽譜方法仿真結果Fig.9 Results of optimal trajectory by GPM

    3.3 改進多階段高斯偽譜方法仿真結果

    為與傳統(tǒng)高斯偽譜方法進行比較,改進多階段高斯偽譜方法所選取的高斯節(jié)點數(shù)目同樣為30,仿真結果如圖10所示。

    (a) 三維軌跡

    (b) 俯視軌跡

    (c) 禁飛區(qū)1附近軌跡

    (d) 禁飛區(qū)2附近軌跡圖10 改進多階段高斯偽譜方法仿真結果Fig.10 Results of optimal trajectory by improved GPM

    3.4 仿真分析

    由圖8(a)、(b)可看出,所得軌跡較為平滑,圖8(c)、(d)說明禁飛區(qū)約束能夠嚴格滿足,圖8(e)說明路徑點約束不能滿足。

    圖9(a)給出了飛行器的三維運動軌跡。可見,所得軌跡較為平滑。圖9(b)~(d)顯示出了傳統(tǒng)高斯偽譜方法的不足,即只能保證所得軌跡在高斯節(jié)點處嚴格滿足約束條件,而無法保證其他部分的軌跡是否滿足禁飛區(qū)等復雜約束。

    由圖10給出的仿真結果可見,通過使用固定節(jié)點和虛擬節(jié)點將軌跡分為了4段,而總的高斯節(jié)點的數(shù)目沒有增加,只是改變了高斯節(jié)點的分布,所得軌跡嚴格滿足路徑點與禁飛區(qū)約束。

    此外,粒子群算法求得最終飛行器飛行時間為1 592 s,傳統(tǒng)高斯偽譜方法和改進高斯偽譜方法的最終結果分別為1 540 s和1 563 s。由此可見,高斯偽譜方法和改進高斯偽譜方法的精度都優(yōu)于粒子群算法,相比改進高斯偽譜方法,傳統(tǒng)高斯偽譜法雖然求得的目標函數(shù)值更小,但由于所得軌跡穿越了禁飛區(qū),因此軌跡無法使用。

    最后,分析不同仿真方法的耗時。粒子群算法的平均耗時約為2 min,但所得軌跡也只是接近全局最優(yōu);高斯偽譜方法的耗時約為1 min,但30個高斯節(jié)點不能保證節(jié)點之間的軌跡滿足約束條件,如果將高斯節(jié)點的數(shù)目增加,雖然所得軌跡精度會提高,但仿真耗時同樣會提高;相比之下,同樣是30個高斯節(jié)點,但本文所提算法的耗時為30 s。這是由于固定節(jié)點的插入將路徑點約束轉化成了終端約束,而虛擬節(jié)點的插入將禁飛區(qū)約束只在某分段模型中進行求解,而傳統(tǒng)高斯偽譜方法是將路徑點和禁飛區(qū)約束在整個模型中進行求解。因此,本文所提算法耗時較小,且精度更高。

    從以上仿真結果,可得出以下結論:

    (1)使用固定節(jié)點將軌跡進行分段,所生成軌跡能夠嚴格滿足路徑點約束;

    (2)使用虛擬節(jié)點將軌跡進行分段,會有更多的高斯節(jié)點分布于虛擬節(jié)點附近,從而實現(xiàn)所生成軌跡能嚴格滿足禁飛區(qū)約束和過程約束;

    (3)所提出的改進多階段高斯偽譜算法只是改變了高斯節(jié)點的分布,而沒有增加高斯節(jié)點的數(shù)目。因此,相比其他類型多階段偽譜方法,本文所提算法復雜度明顯降低。

    4 結論

    為解決復雜約束條件下再入高超聲速滑翔飛行器軌跡優(yōu)化問題,本文提出改進多階段高斯偽譜方法。通過使用固定節(jié)點和虛擬節(jié)點,將軌跡進行分段,在不增加高斯節(jié)點數(shù)目的前提下,改變了高斯節(jié)點的分布,以實現(xiàn)路徑點和禁飛區(qū)附近有較多的高斯節(jié)點分布。仿真結果表明,所得軌跡具有較高的精度,各項約束條件能夠嚴格滿足。

    [1] Huntington G T, Benson D, Rao A V. A comparison of accuracy and computational efficiency of three pseudospectral methods[C]//Proc. of the AIAA Guidance, Navigation, and Control Conference and Exhibit.2007, AIAA 20-23.

    [2] Guo Xiao, Zhu Ming. Direct trajectory optimization based on a mapped Chebyshev pseudospectral method[J]. Chinese Journal of Aeronautics, 2013, 26(2):401-412.

    [3] Ross I M, Fahroo F B. A perspective on methods for trajectory optimization[R].AIAA 2002-4727.

    [4] Bailing Tian. 3DOF ascent phase trajectory optimization for aircraft based on adaptive gauss pseudospectral method[C]//2012 Third International Conference on Intelligent Control and Information Processing, 2012: 431-453.

    [5] 孫勇,張卯瑞,梁曉玲. 求解含負責約束非線性最優(yōu)控制問題的改進Gauss偽譜法[J]. 自動化學報,2013,39(5):672-678.

    [6] 陳琦,王中原,常思江. 基于Gauss偽譜法的滑翔彈道快速優(yōu)化[J]. 彈道學報,2013,21(15):105-108.

    [7] 李柯,聶萬勝,馮必鳴. 助推-滑翔飛行器彈道分段優(yōu)化研究[J]. 指揮控制與仿真,2012,34(5):21-25.

    [8] 湯亮,楊建民,陳風雨. 多約束條件下升力滑翔式再入軌跡優(yōu)化[J].導彈與航天運載技術,2013(1):1-5.

    [9] 曹紅錦,葛致磊. 滑翔增程彈彈道優(yōu)化設計研究[J].四川兵工學報,2013,34(11):8-10.

    [10] 張鼎逆,劉毅. 基于改進遺傳算法和序列二次 規(guī)劃的再入軌跡優(yōu)化[J]. 浙江大學學報,2014,48(1):161-167.

    [11] 雍恩米,唐國金,陳磊. 基于Gauss偽普方法的高超聲速飛行器再入軌跡快速優(yōu)化[J].宇航學報,2008,29(6):1766-1771.

    [12] Jiang Z, Rui Z. Reentry trajectory optimization for hypersonic vehicle satisfying complex constraints[J]. Chinese Journal of Aeronautics, 2013, 26(6): 1544-1553.

    [13] Jorris T R, Cobb R G. Three-dimensional trajectory optimization satisfying waypoint and no-fly zone constraints[J]. Journal of Guidance, Control, and Dynamics, 2009, 32(2): 551-572.

    [14] Jorris T R. Common aero vehicle autonomous reentry trajectory optimization satisfying waypoint and no-fly zone constraints[D]. Ohio, USA: Air Force Institute of Technology, 2007.

    [15] Hu Z D. Research on trajectory planning and guidance for space-based strike weapon[D]. Changsha, China: Graduate School of National University of Defense Technology, 2009.

    [16] Vinh N X. Optimal trajectories in atmospheric flight[D]. New York:Elsevier Scientific Publishing Company,1981.

    [17] 孫勇,段廣仁,張卯瑞,等. 高超聲速飛行器再入過程改進氣動系數(shù)模型[J]. 系統(tǒng)工程與電子技術,2011,33(1):134-137.

    [18] Benson D A, Huntington G T, Thorvaldsen T P, et al. Direct trajectory optimization and costate estimation via an orthogonal collocation method[J]. Journal of Guidance, Control, and Dynamics, 2006,29(6): 1435-1440.

    [19] Huntington G T. Advancement and analysis of a Gauss pseudospectral transcription for optimal control problems[D]. Cambridge, United States: Massachusetts Institute of Technology, 2007.

    [20] Liu Y, Qian Y J, Li J Q, et al. Mars exploring trajectory optimization using gauss pseudo-spectral method[C]//Proc. of the International Conference on Mechatronics and Automation.2012: 2371-2377.

    [21] Yong E M, Chen L, Tang G J. Trajectory optimization of hypersonic gliding reentry vehicle based on the physical programming[J]. Acta Aeronautica et Astronautica Sinica, 2008, 29(5):1092-1097.

    [22] 謝富強,吳浩,唐靈靈. 基于粒子群算法的飛行器再入軌跡優(yōu)化[J].計算技術與自動化,2008,27(4):72-75.

    (編輯:薛永利)

    Trajectory rapid optimization for reentry hypersonic glide vehicle satisfying complex constraints

    WANG Lu,XING Qing-hua,MAO Yi-fan

    (School of Air and Missile Defense, Air Force Engineering University, Xi'an 710051, China)

    As a typical representation of direct method, in dealing with reentry hypersonic glide vehicle trajectory optimization problem under complex constraints including waypoint or no-fly zone constraints, gauss pseudo-spectral method (GPM) can merely ensure that all the constraints can be satisfied strictly at each Legendre-Gauss (LG) node, but cannot ensure the feasibility of the trajectory between each two neighbor nodes. To overcome this problem, an improved multi-segment gauss pseudo-spectral method was proposed. In this method, first of all, the trajectory was divided into multiple segments by two kinds of proposed nodes, fixed node and virtual node, the first one is used to ensure that each waypoint locates at one of LG nodes, and the second one is used to ensure that there will be more LG nodes locating near the no-fly zones, such that all constraints for whole trajectory can be satisfied; then the number of LG nodes for each segment was given; finally, the sequential quadratic programming method was used to solve the model. Simulation results demonstrate that the proposed optimization method can rapidly generate a reentry trajectory satisfying complex constraints.

    reentry hypersonic glide vehicle;trajectory optimization;waypoint;no-fly zone;gauss pseudospectral method

    2015-07-18;

    2015-09-02。

    國家自然科學基金(61272011);全軍軍事類研究生資助課題(2014JY526)。

    王路(1987—),男,博士生,研究方向為地基反臨近空間飛行器指揮決策問題。E-mail:408191081@qq.com

    V448.2

    A

    1006-2793(2016)06-0839-08

    10.7673/j.issn.1006-2793.2016.06.018

    猜你喜歡
    偽譜禁飛區(qū)約束條件
    基于一種改進AZSVPWM的滿調(diào)制度死區(qū)約束條件分析
    矩陣偽譜的新定位集及其在土壤生態(tài)系統(tǒng)的應用
    A literature review of research exploring the experiences of overseas nurses in the United Kingdom (2002–2017)
    大疆更新多邊形禁飛區(qū)策略
    航空模型(2017年7期)2017-08-15 10:24:44
    紊流環(huán)境下四維軌跡優(yōu)化的偽譜方法研究
    飛行力學(2017年1期)2017-02-15 03:00:58
    線性規(guī)劃的八大妙用
    偽譜法及其在飛行器軌跡優(yōu)化設計領域的應用綜述*
    基于Radau偽譜方法的軌跡優(yōu)化
    人臉禁飛區(qū)
    自然與科技(2013年6期)2013-04-29 00:44:03
    一類禁飛區(qū)后方安全撤離軌跡的設計方法研究
    bbb黄色大片| 97碰自拍视频| 精品一区二区三区av网在线观看| 亚洲va在线va天堂va国产| 丰满的人妻完整版| 精品一区二区三区视频在线观看免费| 精品久久久久久久久久免费视频| 久久国内精品自在自线图片| 亚洲av免费在线观看| 日韩人妻高清精品专区| 在线国产一区二区在线| 婷婷丁香在线五月| 久久久久久久精品吃奶| 美女黄网站色视频| 久久欧美精品欧美久久欧美| 欧美精品啪啪一区二区三区| 最新中文字幕久久久久| 全区人妻精品视频| 久久久久久久久中文| 老司机福利观看| 日韩高清综合在线| 亚洲精品粉嫩美女一区| 91av网一区二区| 全区人妻精品视频| 在线播放无遮挡| 亚洲国产精品久久男人天堂| 好男人在线观看高清免费视频| 少妇裸体淫交视频免费看高清| 国产蜜桃级精品一区二区三区| 成人av一区二区三区在线看| 日韩,欧美,国产一区二区三区 | 亚洲欧美清纯卡通| 精品久久国产蜜桃| 日本 av在线| 波多野结衣高清无吗| 国产精品一区二区三区四区免费观看 | 国产精品一区二区性色av| 国产亚洲精品av在线| 99久久成人亚洲精品观看| 精品欧美国产一区二区三| 最近视频中文字幕2019在线8| 免费不卡的大黄色大毛片视频在线观看 | 最新中文字幕久久久久| 悠悠久久av| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲天堂国产精品一区在线| 99热6这里只有精品| 一个人观看的视频www高清免费观看| 日韩欧美三级三区| 亚洲精品一区av在线观看| 婷婷精品国产亚洲av| av专区在线播放| 欧美成人免费av一区二区三区| 一本久久中文字幕| 亚洲精品国产成人久久av| 他把我摸到了高潮在线观看| 99国产极品粉嫩在线观看| 免费av不卡在线播放| 国产亚洲精品综合一区在线观看| 日本三级黄在线观看| 美女 人体艺术 gogo| 中亚洲国语对白在线视频| 日韩强制内射视频| 亚洲av.av天堂| 2021天堂中文幕一二区在线观| 成年版毛片免费区| 国产精品人妻久久久久久| 亚洲男人的天堂狠狠| 国产成人aa在线观看| 99久国产av精品| 亚洲欧美日韩东京热| 久久人妻av系列| 一级av片app| 国产麻豆成人av免费视频| 免费搜索国产男女视频| 狂野欧美激情性xxxx在线观看| 黄色女人牲交| 国产精品久久久久久亚洲av鲁大| 欧美最新免费一区二区三区| 亚洲国产日韩欧美精品在线观看| 九色国产91popny在线| 日本黄色片子视频| 国产淫片久久久久久久久| 亚洲七黄色美女视频| 成人特级av手机在线观看| 久久这里只有精品中国| 国产午夜精品久久久久久一区二区三区 | 亚洲国产欧美人成| 他把我摸到了高潮在线观看| 亚洲精品在线观看二区| 国产老妇女一区| 国产精品99久久久久久久久| av在线老鸭窝| 午夜福利成人在线免费观看| 精品人妻熟女av久视频| 春色校园在线视频观看| 午夜福利欧美成人| 午夜影院日韩av| 国产高潮美女av| 亚洲成人免费电影在线观看| 内射极品少妇av片p| 麻豆一二三区av精品| 国产aⅴ精品一区二区三区波| 一级av片app| 国内精品久久久久久久电影| 久久午夜亚洲精品久久| 欧美成人一区二区免费高清观看| 一区二区三区四区激情视频 | 久久久久久久亚洲中文字幕| 精品久久久久久成人av| 99热6这里只有精品| 午夜精品久久久久久毛片777| 国产成人影院久久av| 一区二区三区四区激情视频 | 免费人成在线观看视频色| 国产大屁股一区二区在线视频| 99热这里只有精品一区| 久久久国产成人精品二区| 欧美xxxx性猛交bbbb| 久久九九热精品免费| 日韩高清综合在线| 成人无遮挡网站| av视频在线观看入口| 亚洲精品日韩av片在线观看| 欧美日韩综合久久久久久 | 波多野结衣高清无吗| 1000部很黄的大片| 三级国产精品欧美在线观看| 亚洲第一区二区三区不卡| 久久香蕉精品热| 欧美精品啪啪一区二区三区| 亚洲,欧美,日韩| 看黄色毛片网站| 九色国产91popny在线| 97热精品久久久久久| 国产精品av视频在线免费观看| 熟妇人妻久久中文字幕3abv| 三级毛片av免费| 国产精品精品国产色婷婷| 亚洲专区国产一区二区| 欧美成人免费av一区二区三区| 乱人视频在线观看| 22中文网久久字幕| 午夜视频国产福利| 真人一进一出gif抽搐免费| 色吧在线观看| 九九在线视频观看精品| 亚洲国产精品久久男人天堂| 欧美色欧美亚洲另类二区| av在线老鸭窝| 九九爱精品视频在线观看| 亚洲第一电影网av| 亚洲国产精品久久男人天堂| 如何舔出高潮| 麻豆av噜噜一区二区三区| 一进一出抽搐动态| 日韩亚洲欧美综合| 国产亚洲av嫩草精品影院| 久久久久九九精品影院| 一区二区三区免费毛片| 久久久久久久精品吃奶| 欧美zozozo另类| 精品免费久久久久久久清纯| 一级毛片久久久久久久久女| 美女xxoo啪啪120秒动态图| 亚洲最大成人av| 中文字幕熟女人妻在线| 噜噜噜噜噜久久久久久91| 国国产精品蜜臀av免费| 丰满乱子伦码专区| 九九热线精品视视频播放| 国产一区二区三区在线臀色熟女| 噜噜噜噜噜久久久久久91| 日韩欧美免费精品| 日韩欧美免费精品| 亚洲七黄色美女视频| 欧美激情在线99| 嫩草影院入口| 国产精品久久视频播放| 午夜免费激情av| 午夜免费激情av| 亚洲国产色片| 亚洲人与动物交配视频| 非洲黑人性xxxx精品又粗又长| 99热6这里只有精品| 欧美丝袜亚洲另类 | av天堂中文字幕网| 在线天堂最新版资源| 国产极品精品免费视频能看的| 不卡一级毛片| 日本免费一区二区三区高清不卡| 亚洲自拍偷在线| 97人妻精品一区二区三区麻豆| 少妇人妻一区二区三区视频| 亚洲自拍偷在线| 少妇丰满av| 精品一区二区免费观看| 国产精品免费一区二区三区在线| 97超视频在线观看视频| 亚洲av成人av| a级毛片a级免费在线| 天美传媒精品一区二区| 少妇猛男粗大的猛烈进出视频 | 免费大片18禁| 免费搜索国产男女视频| 丝袜美腿在线中文| 99久久中文字幕三级久久日本| 国产高潮美女av| 国产伦一二天堂av在线观看| 日韩欧美在线二视频| 亚洲四区av| 国产一区二区三区在线臀色熟女| 18禁黄网站禁片免费观看直播| 国产精品98久久久久久宅男小说| 老师上课跳d突然被开到最大视频| 99精品久久久久人妻精品| 观看美女的网站| 精品一区二区免费观看| 精品人妻偷拍中文字幕| 黄色配什么色好看| 欧美一级a爱片免费观看看| 在线看三级毛片| 成年免费大片在线观看| 亚洲av电影不卡..在线观看| 又黄又爽又刺激的免费视频.| 日日夜夜操网爽| 国产精品亚洲美女久久久| 黄色配什么色好看| 国产一级毛片七仙女欲春2| 日韩欧美国产一区二区入口| 少妇被粗大猛烈的视频| 亚州av有码| 黄色视频,在线免费观看| 男女视频在线观看网站免费| 精品欧美国产一区二区三| 韩国av在线不卡| 内地一区二区视频在线| 成人一区二区视频在线观看| 国产成人一区二区在线| 亚洲av五月六月丁香网| 精品一区二区三区视频在线观看免费| 自拍偷自拍亚洲精品老妇| 国产蜜桃级精品一区二区三区| 床上黄色一级片| 天美传媒精品一区二区| x7x7x7水蜜桃| 久久亚洲真实| 日本一二三区视频观看| 久久午夜福利片| 欧美xxxx性猛交bbbb| av在线亚洲专区| 99久久久亚洲精品蜜臀av| 2021天堂中文幕一二区在线观| 搡女人真爽免费视频火全软件 | 亚洲专区国产一区二区| 日韩一本色道免费dvd| 一级a爱片免费观看的视频| 午夜a级毛片| 欧美一区二区国产精品久久精品| 99国产精品一区二区蜜桃av| av女优亚洲男人天堂| 色av中文字幕| 白带黄色成豆腐渣| 成人国产一区最新在线观看| 国产伦精品一区二区三区视频9| 国产精品自产拍在线观看55亚洲| 午夜老司机福利剧场| 久久久久性生活片| 亚洲成av人片在线播放无| 三级男女做爰猛烈吃奶摸视频| 国产三级中文精品| 香蕉av资源在线| 国产精品自产拍在线观看55亚洲| 久久国内精品自在自线图片| 国产亚洲欧美98| 又黄又爽又刺激的免费视频.| 国产伦一二天堂av在线观看| 国产极品精品免费视频能看的| 美女免费视频网站| 久久人人爽人人爽人人片va| 岛国在线免费视频观看| 午夜激情欧美在线| 亚洲天堂国产精品一区在线| 韩国av一区二区三区四区| 欧美区成人在线视频| 午夜福利欧美成人| 日韩强制内射视频| 久久久色成人| 亚洲美女搞黄在线观看 | 国产 一区 欧美 日韩| 内射极品少妇av片p| 一区二区三区高清视频在线| 久久精品国产自在天天线| 国产成人福利小说| 国模一区二区三区四区视频| 日本黄色片子视频| 联通29元200g的流量卡| 黄色日韩在线| 精品人妻视频免费看| 男女边吃奶边做爰视频| 在线观看午夜福利视频| 蜜桃亚洲精品一区二区三区| 亚洲精华国产精华液的使用体验 | bbb黄色大片| 日本一本二区三区精品| 深夜a级毛片| 91久久精品国产一区二区成人| 久久人妻av系列| 亚洲精品456在线播放app | 午夜福利高清视频| 91麻豆精品激情在线观看国产| 日本a在线网址| 三级毛片av免费| 亚洲欧美日韩高清在线视频| 国产精品1区2区在线观看.| 亚洲精品在线观看二区| 99在线人妻在线中文字幕| 久久久久国内视频| 午夜福利成人在线免费观看| 韩国av一区二区三区四区| ponron亚洲| 五月玫瑰六月丁香| 亚洲五月天丁香| 免费看日本二区| 亚洲午夜理论影院| 免费在线观看影片大全网站| 国产精品一及| 亚洲欧美日韩卡通动漫| xxxwww97欧美| 免费一级毛片在线播放高清视频| 日韩中文字幕欧美一区二区| 国产高清激情床上av| 噜噜噜噜噜久久久久久91| 麻豆成人午夜福利视频| 久久久久久久亚洲中文字幕| 欧美区成人在线视频| 久久婷婷人人爽人人干人人爱| 两个人视频免费观看高清| 亚洲成a人片在线一区二区| 人妻丰满熟妇av一区二区三区| 久久香蕉精品热| 日韩欧美一区二区三区在线观看| 亚洲aⅴ乱码一区二区在线播放| 干丝袜人妻中文字幕| 日日摸夜夜添夜夜添av毛片 | 在线观看午夜福利视频| 99久国产av精品| 变态另类成人亚洲欧美熟女| 亚洲自偷自拍三级| 天天躁日日操中文字幕| 亚洲精品456在线播放app | 在线免费观看的www视频| 午夜爱爱视频在线播放| 国产精品一区二区三区四区免费观看 | 麻豆精品久久久久久蜜桃| 波多野结衣巨乳人妻| 午夜a级毛片| 免费搜索国产男女视频| 亚洲av不卡在线观看| 少妇熟女aⅴ在线视频| 97人妻精品一区二区三区麻豆| 深爱激情五月婷婷| 两人在一起打扑克的视频| 午夜日韩欧美国产| 波多野结衣高清无吗| a级毛片a级免费在线| 无人区码免费观看不卡| 午夜视频国产福利| 亚洲 国产 在线| a级毛片免费高清观看在线播放| 国产淫片久久久久久久久| 中文字幕久久专区| 91久久精品电影网| 色哟哟·www| 亚洲人与动物交配视频| 中文资源天堂在线| 村上凉子中文字幕在线| 欧美性猛交黑人性爽| 精品无人区乱码1区二区| 亚洲精华国产精华液的使用体验 | 国产亚洲精品久久久com| 国产色爽女视频免费观看| 成年女人看的毛片在线观看| 免费不卡的大黄色大毛片视频在线观看 | 一级av片app| 99久久九九国产精品国产免费| 久久精品久久久久久噜噜老黄 | 校园人妻丝袜中文字幕| 国产v大片淫在线免费观看| 看免费成人av毛片| 久久亚洲精品不卡| 一边摸一边抽搐一进一小说| 国产大屁股一区二区在线视频| 日本-黄色视频高清免费观看| 欧美精品国产亚洲| 好男人在线观看高清免费视频| 欧美最黄视频在线播放免费| 久久热精品热| 黄色视频,在线免费观看| 国产单亲对白刺激| 色综合色国产| 久久久久久伊人网av| 最后的刺客免费高清国语| 国产精品一区二区免费欧美| 身体一侧抽搐| 淫秽高清视频在线观看| 成人二区视频| 免费av不卡在线播放| 97超视频在线观看视频| 精品一区二区三区视频在线| 麻豆一二三区av精品| 2021天堂中文幕一二区在线观| 97热精品久久久久久| 赤兔流量卡办理| 国产精品一区二区性色av| 亚洲精品国产成人久久av| 亚洲一区高清亚洲精品| 婷婷亚洲欧美| 久久久久久久久久黄片| 成人亚洲精品av一区二区| 久久香蕉精品热| 麻豆国产97在线/欧美| 美女高潮的动态| 国产白丝娇喘喷水9色精品| 黄色女人牲交| 好男人在线观看高清免费视频| 特大巨黑吊av在线直播| 国产亚洲精品综合一区在线观看| 99久久久亚洲精品蜜臀av| 国产主播在线观看一区二区| 久久九九热精品免费| 欧美高清成人免费视频www| av国产免费在线观看| 一本久久中文字幕| 国产av麻豆久久久久久久| 午夜影院日韩av| 亚洲精华国产精华精| 亚洲人成网站高清观看| 老熟妇仑乱视频hdxx| 噜噜噜噜噜久久久久久91| 特大巨黑吊av在线直播| av.在线天堂| 噜噜噜噜噜久久久久久91| 精品99又大又爽又粗少妇毛片 | 亚洲熟妇熟女久久| 美女高潮的动态| 波野结衣二区三区在线| 九色成人免费人妻av| 成人鲁丝片一二三区免费| 国产高清不卡午夜福利| 99在线人妻在线中文字幕| 性欧美人与动物交配| 国产精品亚洲一级av第二区| 琪琪午夜伦伦电影理论片6080| 亚洲欧美精品综合久久99| 小说图片视频综合网站| 亚洲精品成人久久久久久| 我要搜黄色片| 精品人妻一区二区三区麻豆 | 亚洲欧美日韩东京热| 国产不卡一卡二| 麻豆国产av国片精品| 一个人看视频在线观看www免费| 一进一出好大好爽视频| 无人区码免费观看不卡| 99久久九九国产精品国产免费| 亚洲国产欧洲综合997久久,| 成年免费大片在线观看| 女同久久另类99精品国产91| 淫秽高清视频在线观看| 别揉我奶头 嗯啊视频| 久久久久久九九精品二区国产| 最近中文字幕高清免费大全6 | 欧美激情久久久久久爽电影| 精品午夜福利在线看| 免费看日本二区| 国内精品美女久久久久久| 一个人免费在线观看电影| 国产高清三级在线| 久久人人精品亚洲av| 嫁个100分男人电影在线观看| 中文资源天堂在线| 免费电影在线观看免费观看| 亚洲成a人片在线一区二区| 国产亚洲精品久久久久久毛片| 国产一区二区三区在线臀色熟女| 黄色丝袜av网址大全| 日日啪夜夜撸| 久久久久九九精品影院| 国产精品亚洲美女久久久| 亚洲精品色激情综合| 91精品国产九色| 国产三级中文精品| 午夜精品久久久久久毛片777| 亚洲一区二区三区色噜噜| 精品久久久久久久末码| 成年免费大片在线观看| 成年女人看的毛片在线观看| 欧美日韩乱码在线| 国产精品电影一区二区三区| 日韩精品有码人妻一区| 午夜激情欧美在线| 99视频精品全部免费 在线| 国内毛片毛片毛片毛片毛片| 最近最新中文字幕大全电影3| 深夜a级毛片| 可以在线观看的亚洲视频| 波野结衣二区三区在线| 国产精品精品国产色婷婷| 国产在线精品亚洲第一网站| 久久精品久久久久久噜噜老黄 | 男人和女人高潮做爰伦理| 日韩国内少妇激情av| 有码 亚洲区| 又爽又黄a免费视频| 少妇人妻一区二区三区视频| 成人av在线播放网站| av天堂中文字幕网| 欧美日韩乱码在线| 国产伦精品一区二区三区四那| 黄色女人牲交| 成年女人毛片免费观看观看9| 日韩欧美 国产精品| 中国美女看黄片| 不卡视频在线观看欧美| 欧美日韩中文字幕国产精品一区二区三区| 高清在线国产一区| 91av网一区二区| 三级毛片av免费| 成年女人看的毛片在线观看| 免费人成视频x8x8入口观看| 在线国产一区二区在线| 精品久久国产蜜桃| 又黄又爽又刺激的免费视频.| 在线观看一区二区三区| 97超级碰碰碰精品色视频在线观看| 身体一侧抽搐| 午夜老司机福利剧场| 久久久久国产精品人妻aⅴ院| 91久久精品国产一区二区三区| 国产精品一区二区免费欧美| 久久精品夜夜夜夜夜久久蜜豆| 成年免费大片在线观看| 亚洲专区中文字幕在线| 午夜福利视频1000在线观看| 亚洲色图av天堂| 日本免费a在线| 国产精品人妻久久久久久| 天堂√8在线中文| 国产aⅴ精品一区二区三区波| 午夜福利欧美成人| 国产69精品久久久久777片| 亚洲专区中文字幕在线| 亚洲av免费高清在线观看| 小说图片视频综合网站| 亚洲欧美日韩卡通动漫| 国产伦在线观看视频一区| 亚洲精品国产成人久久av| 免费在线观看日本一区| 久久中文看片网| 久久久久久久久久久丰满 | 舔av片在线| 在线观看舔阴道视频| 一夜夜www| 日韩精品中文字幕看吧| 日日干狠狠操夜夜爽| 国产成人一区二区在线| 露出奶头的视频| 看黄色毛片网站| 久久精品综合一区二区三区| 国产探花在线观看一区二区| 美女高潮喷水抽搐中文字幕| 级片在线观看| 丰满的人妻完整版| 国产精品永久免费网站| avwww免费| 色av中文字幕| www日本黄色视频网| 99久国产av精品| 日本色播在线视频| 女人被狂操c到高潮| 欧美日本亚洲视频在线播放| 久久久久性生活片| 亚洲精品在线观看二区| 一级黄片播放器| 国产精品久久久久久亚洲av鲁大| 精品人妻熟女av久视频| 午夜福利高清视频| 熟女人妻精品中文字幕| 国产高清视频在线播放一区| 99在线视频只有这里精品首页| 精品一区二区三区人妻视频| 久久久色成人| 一区二区三区高清视频在线| 韩国av在线不卡| 3wmmmm亚洲av在线观看| 又爽又黄a免费视频| 色哟哟·www| 国产伦在线观看视频一区| 久久6这里有精品| 欧美三级亚洲精品| 韩国av在线不卡| 亚洲精品成人久久久久久| 亚洲国产精品久久男人天堂| 88av欧美| 看黄色毛片网站| 岛国在线免费视频观看| 精品不卡国产一区二区三区| 真人一进一出gif抽搐免费| 成人毛片a级毛片在线播放| 国产精品嫩草影院av在线观看 | 亚洲熟妇熟女久久| 18+在线观看网站| 亚洲aⅴ乱码一区二区在线播放| 在线观看66精品国产|