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

    基于導(dǎo)航點改進Gauss偽譜法規(guī)劃滑翔導(dǎo)彈航跡

    2016-12-22 00:57:02牛曉潔李邦杰舒健生潘樂飛楊奇松
    彈道學(xué)報 2016年4期
    關(guān)鍵詞:偽譜經(jīng)度滑翔

    牛曉潔,李邦杰,舒健生,潘樂飛,楊奇松

    (火箭軍工程大學(xué) 初級指揮學(xué)院,西安 710025)

    ?

    基于導(dǎo)航點改進Gauss偽譜法規(guī)劃滑翔導(dǎo)彈航跡

    牛曉潔,李邦杰,舒健生,潘樂飛,楊奇松

    (火箭軍工程大學(xué) 初級指揮學(xué)院,西安 710025)

    為了研究滑翔導(dǎo)彈的航跡規(guī)劃問題,利用優(yōu)化理論,規(guī)劃了多約束條件下導(dǎo)航點之間的航跡?;趯?dǎo)航點位置改進Gauss偽譜法(Gauss pseudospectral method,GPM),利用改進GPM離散控制變量和狀態(tài)變量,將最優(yōu)控制問題轉(zhuǎn)化為非線性規(guī)劃問題,利用改進序列二次規(guī)劃算法求解。GPM求解分段航跡規(guī)劃問題,需要迭代獲取下段航跡的飛行時間,降低了算法效率。改進GPM能夠有效彌補單純以時間為自變量帶來的諸多規(guī)劃中的不足,具有更為寬泛的優(yōu)化目標適應(yīng)能力。算例結(jié)果表明,改進后的算法能夠準確、高效地規(guī)劃一條合理的航跡,滿足飛行約束條件。

    滑翔導(dǎo)彈;航跡規(guī)劃;改進高斯偽譜法

    助推-滑翔導(dǎo)彈是高科技戰(zhàn)爭中常實施的殺手锏之一,其航跡規(guī)劃問題作為關(guān)鍵技術(shù)之一值得研究[1]。早期的研究采用解析法求解軌跡優(yōu)化問題,但由于在處理復(fù)雜的非線性系統(tǒng)時力不從心,學(xué)者們開始把研究焦點轉(zhuǎn)向數(shù)值方法,一般分為間接法和直接法[2]。近些年,直接法中的偽譜法由于其高精度、高效率,在軌跡優(yōu)化方面得到了廣泛的應(yīng)用[3]。高斯偽譜法[4]是2005年由麻省理工的Benson給出的一種新的求解最優(yōu)控制問題的直接方法,并從理論上證明了高斯偽譜法與最優(yōu)控制理論中的一階必要條件是一致的。GPM是將狀態(tài)變量和控制變量在一系列Gauss點上進行離散,最終可將最優(yōu)控制問題轉(zhuǎn)化為受一系列代數(shù)約束的參數(shù)優(yōu)化問題,即非線性規(guī)劃問題[5]。

    許多文章都對GPM進行了研究,主要分為3大類:一是給定相應(yīng)的指標泛函設(shè)計航跡,如文獻[6]以最大橫程為優(yōu)化指標,利用高斯偽譜法求解了升力式飛行器的再入軌跡,文獻[7]基于高斯偽譜法以飛行時間為優(yōu)化指標對滑翔彈道進行了優(yōu)化設(shè)計;二是建立模型解決軌跡優(yōu)化問題,如文獻[8]基于高斯偽譜法研究了助推-滑翔飛行器多階段軌跡優(yōu)化問題,文獻[9]針對瞄準對抗問題建立了多階最優(yōu)控制模型,基于GPM對飛機瞄準目標問題進行了研究;三是利用GPM對控制系統(tǒng)進行設(shè)計,如文獻[10]利用Gauss偽譜法對滑翔段到巡航段的過渡段控制率進行了設(shè)計。

    GPM對時間變量進行歸一化處理,求解分段航跡規(guī)劃問題時,需要迭代獲取下段航跡的飛行時間,增加了計算量,降低了算法效率。提出了一種利用導(dǎo)航點位置改進GPM的新思路,以經(jīng)度或者緯度為自變量進行歸一化處理。由于導(dǎo)航點位置可以提前確定,在求解分段航跡規(guī)劃問題時,避免了因為迭代時間導(dǎo)致的算法效率降低的問題。本文中以經(jīng)度或者緯度為自變量的GPM方法,可以有效地彌補單純的以時間為自變量帶來的諸多規(guī)劃中的不便,具有更為寬泛的優(yōu)化目標適應(yīng)能力,以滿足諸如速度控制、能量管理等優(yōu)化目標。文中以經(jīng)度為例進行算法設(shè)計,并進行了仿真驗證。

    1 滑翔段軌跡優(yōu)化問題

    1.1 無動力滑翔運動模型

    考慮地球為旋轉(zhuǎn)圓球時,高超聲速飛行器無動力再入運動方程為[11]

    (1)

    式中:ω,r,v,γ,ψ,λ,B,γc,Fx,Fy分別為地球自轉(zhuǎn)角速度、地心距、速度、飛行路徑角、航向角、經(jīng)度、緯度、傾側(cè)角、氣動升力和阻力。航跡角是速度向量與當(dāng)?shù)厮矫娴膴A角,向上為正。航向角是速度向量在當(dāng)?shù)厮矫嫱队芭c正北方向的夾角,順時針旋轉(zhuǎn)為正。

    文中研究導(dǎo)航點之間的航跡規(guī)劃問題均在滑翔段完成,在北天東坐標系下確定導(dǎo)航點位置后,滑翔導(dǎo)彈飛行的縱向二維平面也就確定了。為了便于研究,該段航跡規(guī)劃時,假設(shè)忽略射面外的運動,僅以二維即只關(guān)注射面內(nèi)的飛行參數(shù),則規(guī)劃式的動力學(xué)模型與動力學(xué)模型都如下式所示:

    (2)

    式中:u,vr,Θ,h分別為切向速度、法向速度、速度傾角和高程。

    1.2 約束條件

    為保證飛行器在滑翔段飛行的穩(wěn)定性,本文考慮了滑翔過程中的相關(guān)約束條件,包括控制量約束和彈道約束:熱流密度、動壓和過載等[12]。

    1.2.1 控制約束

    考慮導(dǎo)彈滑翔過程的實際情況,控制量應(yīng)滿足一定約束,即

    αmin≤α≤αmax

    (3)

    設(shè)定攻角的變化范圍為10°~20°,變化率不超過5(°)/s,基于三自由度模型,假定控制系統(tǒng)能夠保證導(dǎo)彈在豎直平面內(nèi)滑翔飛行,文中取導(dǎo)航點之間導(dǎo)彈飛行的縱向二維平面與東天面的夾角為定值。

    1.2.2 熱流約束

    (4)

    1.2.3 動壓約束

    q=ρv2/2≤qmax

    (5)

    式中:ρ為大氣密度,文中取飛行器可承受的最大動壓qmax=60kPa。

    1.2.4 過載約束

    (6)

    式中:CD,CL,S分別為氣動阻力系數(shù)、升力系數(shù)和氣動參考面積。左邊表達式為法向氣動過載,右邊為飛行器可以承受的最大過載,文中取nmax=4。

    2 改進思路及模型

    2.1 改進思路

    GPM改進背景及思路:在滑翔導(dǎo)彈滑翔過程中,處于某一位置時相對于開始時刻的時間不能很好、快速地獲取,但所處位置的經(jīng)緯度能夠通過計算得出,若是能對經(jīng)緯度進行歸一化處理,將原有GPM中速度對時間求導(dǎo)的約束條件轉(zhuǎn)化為速度對經(jīng)緯度求導(dǎo)的約束條件,將會簡化問題的復(fù)雜程度,提高計算效率。

    2.2 經(jīng)緯度歸一化模型

    已知兩導(dǎo)航點經(jīng)、緯度坐標,可以求得相對于北天東坐標系下兩導(dǎo)航點間的弧線段與Oxz平面之間的夾角σ,以及弧線段在經(jīng)、緯度方向的跨度大小。根據(jù)跨度大小決定以緯度B或經(jīng)度λ為基準進行歸算處理,即若經(jīng)度λ方向跨度大就將經(jīng)度λ歸算到τ。

    (7)

    2.2.1 經(jīng)度歸一化

    改進GPM將2個導(dǎo)航點之間對應(yīng)的經(jīng)度λ標準化為區(qū)間τ∈[-1,1]之間的量,因此對經(jīng)度作如下變換:

    (8)

    可以看出經(jīng)度λ與τ呈現(xiàn)一一對應(yīng)的關(guān)系,并且單調(diào)遞增。

    對應(yīng)的高程變量有:

    (9)

    其他狀態(tài)變量同理可得。

    2.2.2 高斯偽譜法配點選擇

    綜合考慮航跡規(guī)劃問題,程序中采用K階Legendre多項式的根以及τ0=-1作為離散節(jié)點,Legendre多項式為

    (10)

    2.2.3 狀態(tài)變量和控制變量選擇

    在航跡規(guī)劃時,選擇攻角作為控制量。由于K階Legendre多項式共有K+1個根,則作為控制量的攻角可表示為

    α=(α1α2…αK+1)

    (11)

    狀態(tài)量在對應(yīng)時域下的值為

    (12)

    2.2.4 拉格朗日插值函數(shù)

    狀態(tài)變量為

    (13)

    2.2.5 構(gòu)造等式約束

    將式(7)與經(jīng)狀態(tài)變量表達式求導(dǎo)后的函數(shù)式聯(lián)立,可得等式約束式為

    (14)

    2.2.6 目標函數(shù)

    為了提高航跡規(guī)劃的效率,文中選擇最短滑翔時間作為性能指標泛函:

    (15)

    3 求解優(yōu)化問題的算法流程

    3.1 航跡規(guī)劃問題流程分析

    首先利用改進GPM離散最優(yōu)控制問題,利用二次規(guī)劃算法,求解近似原問題的二次規(guī)劃子問題來近似求解原問題。將典型二次規(guī)劃問題進行簡化分析,轉(zhuǎn)化為只含等式約束的最優(yōu)控制問題,利用起作用約束集法調(diào)用拉格朗日算法,迭代更新起作用集和可行解,最終求得航跡規(guī)劃問題的最優(yōu)解。具體流程圖如圖1所示。

    3.2 C++環(huán)境下程序構(gòu)架

    在C++環(huán)境下編寫程序?qū)崿F(xiàn)改進GPM離散最優(yōu)控制問題,使用起作用集法求解二次規(guī)劃子問題。

    ①構(gòu)造等式約束。

    程序中以K階Legendre多項式的根為配點以及τ0=-1作為插值節(jié)點,利用拉格朗日插值函數(shù)求微分來近似式(7)運動模型中的左函數(shù),之后與函數(shù)聯(lián)立,構(gòu)造等式約束式。相對于狀態(tài)變量而言,每個節(jié)點對應(yīng)3個等式約束,K+1個節(jié)點就對應(yīng)3(K+1)個等式約束式。相對于控制變量而言,每個節(jié)點對應(yīng)一個等式約束,所以二次規(guī)劃問題總共對應(yīng)3(K+1)+K+1個等式約束。

    利用高斯積分簡化目標函數(shù),求解等式約束式以及目標函數(shù)的梯度。

    ②求取函數(shù)初值。

    以往SQP法解二次規(guī)劃子問題時,首先需要求解線性規(guī)劃式,判斷給定的初始解是否可行,之后進行下一步的計算。這種方法對初值較為敏感,計算較為復(fù)雜。

    文中在C環(huán)境下編寫程序時,只是給出了正定矩陣,利用拉格朗日函數(shù)求解下降方向向量,同時以此作為二次規(guī)劃子問題的初始可行解。與經(jīng)典序列二次規(guī)劃算法相比,簡化了算法流程,通過仿真驗證了算法的高效性和合理性。

    ③二次規(guī)劃尋優(yōu)。

    首先根據(jù)等式約束,確定初始起作用集和非起作用集。使用拉格朗日法求解只含有起作用集約束的線性規(guī)劃問題,根據(jù)相關(guān)判定條件迭代更新起作用集和問題可行解,最終求出二次規(guī)劃子問題的最優(yōu)解,也即是原先二次規(guī)劃問題的最優(yōu)解。

    4 算例分析

    以遠程高超聲速滑翔式再入飛行器為仿真對象,氣動參數(shù)采用美國波音公司98年設(shè)計的CAV-H相關(guān)參數(shù)擬合得到[15],其最大升阻比為3.5;氣動參考面積為0.484m2;質(zhì)量為907.2kg;最大飛行攻角和最大升阻比攻角分別為20°和10°;起始導(dǎo)航點經(jīng)度、緯度、高程分別為29°,72°和60.37km;目標導(dǎo)航點經(jīng)緯度坐標、高程分別為87°,74°和48.8km;滑翔開始時刻的速度為6 606m/s,最大飛行高度為70km,約束條件如1.2節(jié)所示。

    取5階LG點以及τ0=-1作為插值節(jié)點,優(yōu)化計算在CPU為3.4GHz/Corei3,操作系統(tǒng)為WindowsXP系統(tǒng)的微機上實現(xiàn),VisualStudio2010環(huán)境下編寫程序,以時域變換為基礎(chǔ)的Gauss偽譜法求解上述算例需要0.078s,使用文中改進后的算法規(guī)劃航跡的時間為0.016s。兩者對比可知,后者在滿足飛行約束條件的基礎(chǔ)上大大縮短了航跡規(guī)劃所需的時間。

    優(yōu)化計算結(jié)果見圖2。圖2(a)中攻角變化曲線表明,滑翔導(dǎo)彈長時間保持最大升阻比攻角飛行,在起始導(dǎo)航點和目標導(dǎo)航點位置確定的情況下保證了最短飛行時間這一優(yōu)化指標;圖2(b)中速度傾角始終為負值,說明滑翔導(dǎo)彈始終朝向目標點飛行,保證了最短飛行時間的優(yōu)化指標;圖2(c)、圖2(d)中速度、飛行高度都能夠滿足終端約束條件限制;圖2(e)、圖2(f)、圖2(g)中過載、動壓和熱流密度的變化均在給定的變化范圍內(nèi),滿足飛行過程中的約束條件。仿真結(jié)果表明,在滿足相同約束條件的前提下,與以往相比,利用改進后的GPM能夠準確、高效地規(guī)劃出一條合理的航跡。

    圖2 航跡規(guī)劃結(jié)果

    5 結(jié)束語

    文中基于導(dǎo)航點位置改進GPM,利用改進后的算法對滑翔導(dǎo)彈導(dǎo)航點之間的航跡進行了優(yōu)化求解,并選取最短飛行時間作為優(yōu)化指標,在滿足各種飛行約束的條件下,將攻角作為控制變量進行優(yōu)化,得到了滑翔段導(dǎo)航點之間的最優(yōu)航跡。從本文的研究中可以看出,文中改進后的算法與以往GPM相比,跳出了飛行時間的限制,在已知導(dǎo)航點位置后,滑翔導(dǎo)彈的縱向飛行平面也就隨之確定,簡化了運動模型,減少了計算量,縮短了求解航跡規(guī)劃問題所用的時間。文中的研究成果可為滑翔導(dǎo)彈在線航跡規(guī)劃提供一定的技術(shù)支撐,但文中改進算法是在地球為旋轉(zhuǎn)圓球的基礎(chǔ)上提出的,與實際不太相符,下一步將以地球橢球為背景設(shè)計研究相應(yīng)的改進算法。

    [1] 孟中杰,黃攀峰,閆杰.基于改進稀疏A*算法的高超聲速飛行器航跡規(guī)劃技術(shù)[J].西北工業(yè)大學(xué)學(xué)報,2010,28(2):182-186. MENG Zhong-jie,HUANG Pan-feng,YAN Jie.Exploring trajectory planning for hypersonic vehicle using improved sparse A* algorithm[J].Journal of Northwestern Polytechnical University,2010,28(2):182-186.(in Chinese)

    [2]BETTS J T.Survey numerical methods for trajectory optimiza-tion[J].AIAA Journal of Guidance Control and Dynamics,1998,21(2):193-207.

    [3]胡松啟,陳雨.偽譜法在飛行器軌跡優(yōu)化中應(yīng)用分析[J].火箭推進,2014,40(5):61-68. HU Song-qi,CHEN Yu.Analysis of pseudo-spectral methods applied to aircraft trajectory optimization[J].Journal of Rocket Propulsion,2014,40(5):61-68.(in Chinese)

    [4]BENSON D A.A Gauss pseudospectral transcription for optimal control[D].Massachusetts:Massachusetts Institute of Technology,2005.

    [5]唐國金,羅亞中,雍恩米.航天器軌跡優(yōu)化理論、方法及應(yīng)用[M].北京:科學(xué)出版社,2011. TANG Guo-jin,LUO Ya-zhong,YONG En-mi.Theory,method and application of spacecraft trajectory optimization[M].Beijing:Science Press,2011.(in Chinese)

    [6]周文雅,楊滌,李順利.利用高斯偽譜法求解具有最大橫程的再入軌跡[J].系統(tǒng)工程與電子技術(shù),2010,32(5):1 038-1 042. ZHOU Wen-ya,YANG Di,LI Shun-li.Solution of reentry trajectory with maximum cross range by using Gauss pseudospectral method[J].Systems Engineering and Electronics,2010,32(5):1 038-1 042.(in Chinese)

    [7]李鐵鵬,王中原,陳琦,等.基于高斯偽譜法的滑翔彈道優(yōu)化算法研究[J].彈箭與制導(dǎo)學(xué)報,2014,34(2):113-116. LI Tie-peng,WANG Zhong-yuan,CHEN Qi,et al.Study of gliding trajectory optimization algorithm based on Gauss pseudospectral method[J].Journal of Projectiles,Rockets,Missiles and Guidance,2014,34(2):113-116.(in Chinese)

    [8]GUAN Cheng-qi,CHEN Cong.Multiphase path-constrained trajectory optimization for the boost-glide vehicle via the Gauss pseudospectral method[J].Journal of Astronautics,2010,31(11):2 512-2 518.

    [9]程建鋒,董新民,薛建平,等.基于Gauss偽譜法的飛機最優(yōu)目標瞄準控制[J].計算機應(yīng)用,2013,33(11):3 291-3 295. CHENG Jian-feng,DONG Xin-min,XUE Jian-Ping,et al.Aircraft optimal target aiming control based on Gauss pseudospectral method[J].Journal of Computer Applications,2013,33(11):3 291-3 295.(in Chinese)

    [10]汪立新,冷杉,王建華.基于Gauss偽譜法的導(dǎo)彈過渡段控制律設(shè)計[J].控制工程,2013,20(5):906-909. WANG Li-xin,LENG Shan,WANG Jian-hua.Design of control law of transition zone of missile based on Gauss pseudospectral method[J].Control Engineering of China,2013,20(5):906-909.(in Chinese)

    [11]BOLLINO K P,ROSS I M.A pseudospectral feedback method for real-time optimal guidance of reentry vehicles[C]//Proceedings of the American Control Conference.New York:American Control Conference,2007:3 861-3 867.

    [12]馮林爽,劉磊,王永驥.基于PSO和SQP混合優(yōu)化算法的高超聲速飛行器軌跡優(yōu)化[C]//2015 IEEE Chinese Control and Decision Conference.Qingdao:IEEE,2015:4 518-4 522. FENG Lin-shuang,LIU Lei,WANG Yong-ji.Trajectories optimization of hypersonic vehicle based on a hybrid optimization algorithm of PSO and SQP[C]//2015 IEEE Chinese Control and Decision Conference.Qingdao:IEEE,2015:4 518-4 522.(in Chinese)

    [13]張科南,周浩,陳萬春.高超聲速飛行器多約束多種機動突防模式彈道規(guī)劃[J].彈道學(xué)報,2012,24(3):85-90. ZHANG Ke-nan,ZHOU Hao,CHEN Wan-chun.Trajectory planning for hypersonic vehicle with multiple constraints and multiple maneuvering penetration strategies[J].Journal of Ballistics,2012,24(3):85-90.(in Chinese)

    [14]劉建業(yè),曾慶化,趙偉,等.導(dǎo)航系統(tǒng)理論與應(yīng)用[M].西安:西北工業(yè)大學(xué)出版社,2010. LIU Jian-ye,ZENG Qing-hua,ZHAO Wei,et al.Navigation system theory and application[M].Xi’an:Northwestern Polytechnic University Press,2010.(in Chinese)

    [14]DARBY C L,HAGER W W,RAO A V.An hp-adaptive pseudospectral method for solving optimal control problems[J].Optimal Control Applications and Methods,2011,32(4):476-502.

    Trajectory Planning for Gliding Missile Via Improved Gauss Pseudospectral Method Based on Navigation Points

    NIU Xiao-jie,LI Bang-jie,SHU Jian-sheng,PAN Le-fei,YANG Qi-song

    (Primary Command College,Rocket Force University,Xi’an 710025,China)

    In order to study the trajectory planning problem of glide missile,the trajectory of the navigation points under the multi-constraint conditions was planned by using optimization theory.The improved Gauss pseudospectral method(GPM)was based on navigation points.The control variables and state variables were discretized by improved GPM,and the optimal control problem was converted into a nonlinear programming problem,and the problem was solved by using the improved sequential quadratic programming(SQP)method.When the segment path is planned by GPM,the flight time of the next trajectory is obtained via iterative method,and the efficiency of the algorithm is reduced.The improved GPM can effectively compensate the lack of time being taken as the independent variable in many trajectory planning,and it has a wider adaptability to optimize the target.The results show that the improved algorithm can accurately and efficiently plan a reasonable trajectory to meet the flight constraints.

    glide missile;trajectory planning;improved Gauss pseudospectral method(GPM)

    2016-06-30

    牛曉潔(1992- ),男,碩士研究生,研究方向為飛行器總體、結(jié)構(gòu)分析與飛行力學(xué)。E-mail:18706766973@163.com。

    V412.4

    A

    1004-499X(2016)04-0036-06

    猜你喜歡
    偽譜經(jīng)度滑翔
    蘭花螳螂會滑翔
    巧用四步法 妙解地方時
    巧用四步法 妙解地方時
    攻天掠地的先鋒武器——滑翔導(dǎo)彈
    矩陣偽譜的新定位集及其在土壤生態(tài)系統(tǒng)的應(yīng)用
    讓你心跳加速的“滑翔飛板”
    紊流環(huán)境下四維軌跡優(yōu)化的偽譜方法研究
    偽譜法及其在飛行器軌跡優(yōu)化設(shè)計領(lǐng)域的應(yīng)用綜述*
    空中滑翔大比拼(下)——滑翔傘
    基于Radau偽譜方法的軌跡優(yōu)化
    亚洲第一电影网av| tocl精华| 久久精品国产清高在天天线| 欧美bdsm另类| 男女下面进入的视频免费午夜| 欧美黄色片欧美黄色片| 嫁个100分男人电影在线观看| 桃色一区二区三区在线观看| 99精品在免费线老司机午夜| 午夜福利在线观看免费完整高清在 | 亚洲av二区三区四区| 亚洲成人精品中文字幕电影| 久久久久久久午夜电影| 国产中年淑女户外野战色| 真实男女啪啪啪动态图| 精品国产超薄肉色丝袜足j| 久久亚洲精品不卡| 夜夜爽天天搞| 日本在线视频免费播放| 九九在线视频观看精品| 中文字幕人妻丝袜一区二区| 午夜精品在线福利| 美女被艹到高潮喷水动态| 欧美区成人在线视频| 国产一区二区在线观看日韩 | 午夜精品一区二区三区免费看| 变态另类成人亚洲欧美熟女| 99久久精品热视频| 黄片大片在线免费观看| 高潮久久久久久久久久久不卡| 欧美区成人在线视频| www.色视频.com| 亚洲黑人精品在线| av黄色大香蕉| 日本黄大片高清| 丰满人妻一区二区三区视频av | 搞女人的毛片| 无遮挡黄片免费观看| 噜噜噜噜噜久久久久久91| 操出白浆在线播放| 人人妻人人看人人澡| 三级国产精品欧美在线观看| 天堂√8在线中文| 久久性视频一级片| av专区在线播放| 三级国产精品欧美在线观看| 90打野战视频偷拍视频| 成人鲁丝片一二三区免费| 欧美日本视频| av视频在线观看入口| 夜夜夜夜夜久久久久| 亚洲午夜理论影院| 亚洲成人中文字幕在线播放| 免费大片18禁| 男女视频在线观看网站免费| 亚洲人成网站在线播| 欧美一级a爱片免费观看看| 一个人看视频在线观看www免费 | 国产 一区 欧美 日韩| 国产精品一区二区三区四区免费观看 | 亚洲最大成人手机在线| 一级毛片高清免费大全| 亚洲美女黄片视频| 最近最新中文字幕大全免费视频| 久久久久精品国产欧美久久久| 国产老妇女一区| www.色视频.com| 一夜夜www| 午夜精品久久久久久毛片777| a级毛片a级免费在线| 午夜免费成人在线视频| av福利片在线观看| 中文字幕高清在线视频| 亚洲av中文字字幕乱码综合| 黄色丝袜av网址大全| 午夜精品一区二区三区免费看| 日本熟妇午夜| 12—13女人毛片做爰片一| 国产黄片美女视频| 伊人久久精品亚洲午夜| 又黄又粗又硬又大视频| 我的老师免费观看完整版| 人人妻,人人澡人人爽秒播| 亚洲美女视频黄频| 少妇的逼好多水| 中文字幕人妻丝袜一区二区| 国产私拍福利视频在线观看| 国产精品免费一区二区三区在线| АⅤ资源中文在线天堂| 一个人免费在线观看的高清视频| 听说在线观看完整版免费高清| 免费大片18禁| 亚洲av日韩精品久久久久久密| 麻豆国产97在线/欧美| 免费观看人在逋| 日本一二三区视频观看| 99久久综合精品五月天人人| 亚洲在线自拍视频| 变态另类丝袜制服| 午夜a级毛片| 丰满人妻一区二区三区视频av | 久久99热这里只有精品18| 88av欧美| 一夜夜www| 男人的好看免费观看在线视频| 亚洲精品国产精品久久久不卡| 丰满乱子伦码专区| 国产一区二区三区在线臀色熟女| 偷拍熟女少妇极品色| 成人午夜高清在线视频| 色播亚洲综合网| 真人一进一出gif抽搐免费| 一本一本综合久久| 国产真实乱freesex| 桃色一区二区三区在线观看| 首页视频小说图片口味搜索| 首页视频小说图片口味搜索| 亚洲性夜色夜夜综合| 最近最新中文字幕大全免费视频| 少妇人妻精品综合一区二区 | 男人舔女人下体高潮全视频| 国产精品,欧美在线| 色精品久久人妻99蜜桃| 免费在线观看成人毛片| 国产成人a区在线观看| 成人国产综合亚洲| 12—13女人毛片做爰片一| 一进一出好大好爽视频| 舔av片在线| 久久性视频一级片| 国产精品久久电影中文字幕| 18+在线观看网站| 久久精品91蜜桃| 国产一区二区在线观看日韩 | 老司机午夜福利在线观看视频| 国内精品久久久久精免费| 手机成人av网站| 国产男靠女视频免费网站| 亚洲18禁久久av| 免费搜索国产男女视频| 免费在线观看成人毛片| 一级毛片女人18水好多| 免费av不卡在线播放| 国产真实伦视频高清在线观看 | 一个人看的www免费观看视频| 日韩精品青青久久久久久| 免费电影在线观看免费观看| 69人妻影院| 国产97色在线日韩免费| 国产亚洲精品av在线| 亚洲成人久久性| 搡女人真爽免费视频火全软件 | 国产精品久久久久久久久免 | 日日干狠狠操夜夜爽| 国产精品久久久久久人妻精品电影| 免费人成视频x8x8入口观看| 国产成人影院久久av| 国产精品美女特级片免费视频播放器| 亚洲国产日韩欧美精品在线观看 | 亚洲在线自拍视频| 91av网一区二区| 欧美xxxx黑人xx丫x性爽| 精品午夜福利视频在线观看一区| 免费无遮挡裸体视频| 亚洲熟妇中文字幕五十中出| 欧美成人免费av一区二区三区| 亚洲人成网站在线播放欧美日韩| 噜噜噜噜噜久久久久久91| 欧美另类亚洲清纯唯美| 国内精品一区二区在线观看| 欧美av亚洲av综合av国产av| 欧洲精品卡2卡3卡4卡5卡区| 亚洲男人的天堂狠狠| 国产激情偷乱视频一区二区| 在线国产一区二区在线| 国产精品98久久久久久宅男小说| 每晚都被弄得嗷嗷叫到高潮| 国产伦精品一区二区三区视频9 | 日本在线视频免费播放| 午夜福利在线在线| 最近最新免费中文字幕在线| 国产探花在线观看一区二区| 国产高清激情床上av| 1024手机看黄色片| 免费无遮挡裸体视频| 亚洲内射少妇av| 免费看日本二区| 成年女人看的毛片在线观看| 成人午夜高清在线视频| 欧美大码av| 欧美区成人在线视频| 一区二区三区高清视频在线| 一本综合久久免费| 亚洲 欧美 日韩 在线 免费| 97人妻精品一区二区三区麻豆| a级毛片a级免费在线| e午夜精品久久久久久久| 欧美绝顶高潮抽搐喷水| 热99在线观看视频| 色精品久久人妻99蜜桃| 亚洲天堂国产精品一区在线| 久久久久免费精品人妻一区二区| 麻豆成人av在线观看| 美女免费视频网站| 长腿黑丝高跟| 国产高清有码在线观看视频| 好男人在线观看高清免费视频| 国产精品一区二区免费欧美| 亚洲成人久久性| 最近最新中文字幕大全免费视频| 亚洲国产欧美网| 国产蜜桃级精品一区二区三区| 在线观看日韩欧美| 国产精品亚洲一级av第二区| 黄色视频,在线免费观看| 精品人妻一区二区三区麻豆 | 精品免费久久久久久久清纯| 中文字幕av在线有码专区| 日本 av在线| 少妇人妻精品综合一区二区 | 久久久久亚洲av毛片大全| 免费看十八禁软件| 免费看光身美女| 天堂动漫精品| 亚洲第一欧美日韩一区二区三区| 日韩人妻高清精品专区| 桃红色精品国产亚洲av| 精品福利观看| 乱人视频在线观看| 日日夜夜操网爽| eeuss影院久久| av女优亚洲男人天堂| 国内精品一区二区在线观看| 极品教师在线免费播放| 精品久久久久久久人妻蜜臀av| 99久久精品热视频| 国产日本99.免费观看| 国产精品亚洲一级av第二区| 亚洲aⅴ乱码一区二区在线播放| av专区在线播放| 草草在线视频免费看| 午夜激情欧美在线| 国产伦一二天堂av在线观看| 国产亚洲精品久久久久久毛片| 怎么达到女性高潮| 麻豆久久精品国产亚洲av| 麻豆国产97在线/欧美| 免费在线观看成人毛片| www日本在线高清视频| 丰满的人妻完整版| 成年人黄色毛片网站| 免费看日本二区| 成人国产综合亚洲| 久久精品综合一区二区三区| 国产在线精品亚洲第一网站| 色综合亚洲欧美另类图片| 最新在线观看一区二区三区| 国内少妇人妻偷人精品xxx网站| 中亚洲国语对白在线视频| 欧美日韩国产亚洲二区| 色综合站精品国产| 怎么达到女性高潮| 国产精品,欧美在线| 特大巨黑吊av在线直播| 亚洲成人精品中文字幕电影| 97超级碰碰碰精品色视频在线观看| 亚洲国产高清在线一区二区三| 精品免费久久久久久久清纯| 国产成+人综合+亚洲专区| 老汉色av国产亚洲站长工具| 法律面前人人平等表现在哪些方面| 午夜免费成人在线视频| 国产中年淑女户外野战色| 精品一区二区三区视频在线观看免费| 女人十人毛片免费观看3o分钟| 国产伦精品一区二区三区四那| 免费电影在线观看免费观看| av黄色大香蕉| 国产一区在线观看成人免费| 性欧美人与动物交配| 日韩人妻高清精品专区| 精品一区二区三区视频在线 | 欧美在线一区亚洲| 99久久99久久久精品蜜桃| 免费av毛片视频| 少妇的丰满在线观看| 午夜福利高清视频| 欧美日韩福利视频一区二区| 听说在线观看完整版免费高清| 黄色丝袜av网址大全| 丁香欧美五月| 成年女人毛片免费观看观看9| 日本免费一区二区三区高清不卡| 一级作爱视频免费观看| 国内精品美女久久久久久| 18美女黄网站色大片免费观看| 国产精品99久久久久久久久| 超碰av人人做人人爽久久 | 精品国产美女av久久久久小说| 亚洲欧美精品综合久久99| 亚洲真实伦在线观看| 国产淫片久久久久久久久 | 日本五十路高清| 成年人黄色毛片网站| 一级黄色大片毛片| 一边摸一边抽搐一进一小说| www日本在线高清视频| 99国产综合亚洲精品| 天堂av国产一区二区熟女人妻| 亚洲18禁久久av| 精品人妻一区二区三区麻豆 | 99久久综合精品五月天人人| 久久久久亚洲av毛片大全| 亚洲aⅴ乱码一区二区在线播放| 国产97色在线日韩免费| 亚洲精品国产精品久久久不卡| 国产亚洲av嫩草精品影院| 亚洲av中文字字幕乱码综合| 亚洲国产精品sss在线观看| 久久久久久久精品吃奶| 国产色爽女视频免费观看| 1000部很黄的大片| 欧美性猛交╳xxx乱大交人| 熟女少妇亚洲综合色aaa.| 国产精品1区2区在线观看.| 久久这里只有精品中国| 国产亚洲精品久久久com| 中文字幕av在线有码专区| 午夜福利18| 亚洲av熟女| 亚洲在线观看片| 久久国产乱子伦精品免费另类| 一个人免费在线观看电影| 观看免费一级毛片| 神马国产精品三级电影在线观看| 午夜激情福利司机影院| 一区二区三区高清视频在线| 国产精品美女特级片免费视频播放器| 亚洲专区中文字幕在线| 国产美女午夜福利| 国内少妇人妻偷人精品xxx网站| 欧美色视频一区免费| 国产精品久久久久久久电影 | 激情在线观看视频在线高清| 久久久久久久精品吃奶| 亚洲av二区三区四区| 亚洲国产精品成人综合色| 757午夜福利合集在线观看| 叶爱在线成人免费视频播放| 国产精品99久久久久久久久| av视频在线观看入口| 天天一区二区日本电影三级| 69av精品久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 欧美性猛交╳xxx乱大交人| 欧美高清成人免费视频www| 亚洲精品国产精品久久久不卡| 脱女人内裤的视频| 女生性感内裤真人,穿戴方法视频| 99热精品在线国产| 久久久久久国产a免费观看| 看免费av毛片| 麻豆国产97在线/欧美| 美女被艹到高潮喷水动态| 亚洲精品粉嫩美女一区| 亚洲国产欧美网| 亚洲狠狠婷婷综合久久图片| 在线国产一区二区在线| 国内毛片毛片毛片毛片毛片| 欧美精品啪啪一区二区三区| 国产精品 国内视频| 亚洲色图av天堂| 欧美激情在线99| 99久久九九国产精品国产免费| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区三区四区免费观看 | 黑人欧美特级aaaaaa片| 天堂√8在线中文| 脱女人内裤的视频| a级毛片a级免费在线| 日日夜夜操网爽| 男女做爰动态图高潮gif福利片| 在线观看午夜福利视频| 日本免费a在线| 精品一区二区三区视频在线 | 国产一区在线观看成人免费| 国产真实乱freesex| 日韩精品中文字幕看吧| 国产av不卡久久| 老司机午夜福利在线观看视频| 悠悠久久av| 欧美在线一区亚洲| 看黄色毛片网站| www.色视频.com| 成年免费大片在线观看| 成年女人看的毛片在线观看| 51午夜福利影视在线观看| 成人一区二区视频在线观看| 五月伊人婷婷丁香| 18禁黄网站禁片午夜丰满| 日本成人三级电影网站| 国内毛片毛片毛片毛片毛片| 免费观看精品视频网站| 久久国产乱子伦精品免费另类| 18禁黄网站禁片午夜丰满| 动漫黄色视频在线观看| 欧美激情在线99| 手机成人av网站| 18美女黄网站色大片免费观看| 中文亚洲av片在线观看爽| 高清毛片免费观看视频网站| 老熟妇仑乱视频hdxx| 黄色片一级片一级黄色片| 欧美又色又爽又黄视频| 日韩亚洲欧美综合| 少妇熟女aⅴ在线视频| 黄色成人免费大全| 亚洲精品国产精品久久久不卡| www国产在线视频色| 亚洲中文字幕一区二区三区有码在线看| 最近视频中文字幕2019在线8| 国产精品永久免费网站| 国产三级中文精品| 性色avwww在线观看| 99国产精品一区二区三区| 麻豆国产97在线/欧美| 丰满人妻熟妇乱又伦精品不卡| 熟女电影av网| svipshipincom国产片| 婷婷精品国产亚洲av| 亚洲性夜色夜夜综合| 制服人妻中文乱码| 国内精品美女久久久久久| 亚洲黑人精品在线| 在线十欧美十亚洲十日本专区| 18禁在线播放成人免费| 国产伦精品一区二区三区四那| 亚洲,欧美精品.| 亚洲av不卡在线观看| 黄色日韩在线| 天堂√8在线中文| 成年女人毛片免费观看观看9| 国产乱人视频| 男插女下体视频免费在线播放| 手机成人av网站| а√天堂www在线а√下载| 久久精品国产综合久久久| 哪里可以看免费的av片| 色av中文字幕| 国产一级毛片七仙女欲春2| 亚洲av第一区精品v没综合| 久久精品国产综合久久久| 麻豆成人av在线观看| 男人舔女人下体高潮全视频| 日韩亚洲欧美综合| 国产视频一区二区在线看| 亚洲熟妇熟女久久| 熟女人妻精品中文字幕| 亚洲人成网站在线播放欧美日韩| 日韩av在线大香蕉| 熟妇人妻久久中文字幕3abv| 国产精品亚洲美女久久久| 亚洲国产欧美人成| 高潮久久久久久久久久久不卡| av欧美777| 亚洲av不卡在线观看| 亚洲真实伦在线观看| 午夜精品久久久久久毛片777| 国产毛片a区久久久久| 最近视频中文字幕2019在线8| 国产成人av激情在线播放| 久久久成人免费电影| 亚洲精品乱码久久久v下载方式 | 免费大片18禁| 在线播放国产精品三级| 亚洲国产精品sss在线观看| 亚洲人成电影免费在线| 熟女电影av网| 免费看十八禁软件| 欧美日韩精品网址| 亚洲国产精品成人综合色| 亚洲av五月六月丁香网| 男女视频在线观看网站免费| 日韩精品中文字幕看吧| 偷拍熟女少妇极品色| 久久精品亚洲精品国产色婷小说| 天美传媒精品一区二区| 又黄又爽又免费观看的视频| 国产黄色小视频在线观看| 日韩国内少妇激情av| 69av精品久久久久久| 岛国在线观看网站| 噜噜噜噜噜久久久久久91| 精品国内亚洲2022精品成人| 美女高潮喷水抽搐中文字幕| 亚洲av电影不卡..在线观看| 婷婷精品国产亚洲av在线| 精品熟女少妇八av免费久了| 两个人的视频大全免费| 亚洲人成电影免费在线| 午夜福利高清视频| 国产成人aa在线观看| 日韩欧美 国产精品| 99精品久久久久人妻精品| 99视频精品全部免费 在线| a在线观看视频网站| 国产黄片美女视频| 国产精品久久久久久人妻精品电影| 伊人久久精品亚洲午夜| 日日干狠狠操夜夜爽| 色老头精品视频在线观看| 18禁在线播放成人免费| 国产精品99久久久久久久久| 搡老妇女老女人老熟妇| 少妇高潮的动态图| 一区二区三区免费毛片| 亚洲人成网站在线播放欧美日韩| 黄片大片在线免费观看| 精品国产超薄肉色丝袜足j| 美女cb高潮喷水在线观看| 精品国产亚洲在线| 久久午夜亚洲精品久久| 亚洲精品在线美女| 一本一本综合久久| 国产高清视频在线播放一区| bbb黄色大片| 嫩草影视91久久| 在线看三级毛片| 亚洲自拍偷在线| 免费在线观看成人毛片| 午夜福利在线在线| 可以在线观看的亚洲视频| 国产免费一级a男人的天堂| 久久久国产成人免费| 亚洲国产高清在线一区二区三| 黄片大片在线免费观看| 成年女人看的毛片在线观看| 真人一进一出gif抽搐免费| 三级国产精品欧美在线观看| 青草久久国产| 免费在线观看影片大全网站| 久久久国产成人免费| 99久久久亚洲精品蜜臀av| 国产精品美女特级片免费视频播放器| 亚洲一区二区三区色噜噜| 毛片女人毛片| 91久久精品电影网| 少妇人妻精品综合一区二区 | 狂野欧美激情性xxxx| 在线观看一区二区三区| www.熟女人妻精品国产| 久久久久久久久久黄片| 亚洲欧美日韩高清在线视频| 少妇的逼好多水| 亚洲成av人片在线播放无| 男女之事视频高清在线观看| 网址你懂的国产日韩在线| 十八禁人妻一区二区| 国内精品美女久久久久久| 久久久国产成人免费| 午夜免费成人在线视频| www国产在线视频色| 国产高清激情床上av| 精品一区二区三区av网在线观看| 亚洲精品美女久久久久99蜜臀| 真人做人爱边吃奶动态| 超碰av人人做人人爽久久 | 欧美bdsm另类| 国产av一区在线观看免费| АⅤ资源中文在线天堂| 欧美高清成人免费视频www| 日本免费a在线| 麻豆成人av在线观看| 亚洲国产色片| 村上凉子中文字幕在线| 99热这里只有精品一区| 国产私拍福利视频在线观看| 日日摸夜夜添夜夜添小说| 国内毛片毛片毛片毛片毛片| 免费无遮挡裸体视频| 日本一二三区视频观看| 日韩精品青青久久久久久| 最新在线观看一区二区三区| 国产真人三级小视频在线观看| 亚洲天堂国产精品一区在线| 国产av不卡久久| 91在线观看av| 亚洲av熟女| 熟女电影av网| 一级作爱视频免费观看| 日本五十路高清| 欧美一级毛片孕妇| 久久婷婷人人爽人人干人人爱| av黄色大香蕉| 哪里可以看免费的av片| 国产午夜精品论理片| 啪啪无遮挡十八禁网站| 欧美三级亚洲精品| 真人一进一出gif抽搐免费| 男人舔奶头视频| 特大巨黑吊av在线直播| 丰满的人妻完整版| 国产在线精品亚洲第一网站| 99在线人妻在线中文字幕| 波多野结衣巨乳人妻| av视频在线观看入口| 亚洲成人免费电影在线观看| 99精品久久久久人妻精品| 日韩成人在线观看一区二区三区| 高清在线国产一区| 有码 亚洲区| 日本熟妇午夜| 亚洲成av人片免费观看| 久久久久亚洲av毛片大全| 亚洲第一欧美日韩一区二区三区|