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

    高超聲速滑翔飛行器滑翔段初始狀態(tài)的唯一性和最優(yōu)性分析*

    2020-10-28 06:18:20孟繁卿田康生
    國防科技大學學報 2020年5期
    關鍵詞:滑翔超聲速彈道

    孟繁卿,田康生

    (1. 空軍預警學院 研究生大隊, 湖北 武漢 430019; 2. 空軍預警學院 預警情報系, 湖北 武漢 430019)

    臨近空間高超聲速滑翔飛行器通常是指飛行馬赫數(shù)大于5、在距地面20~100 km的高度依靠空氣動力滑翔的一種飛行器。臨近空間高超聲速滑翔飛行器的顯著特點是飛行高度比彈道導彈低,飛行速度遠高于普通航空飛行器。因此高超聲速滑翔飛行器兼具彈道導彈的高速、航空飛行器的靈活機動的特點,被視為改變戰(zhàn)爭規(guī)則的“撒手锏”武器。高超聲速滑翔飛行器可多次變軌,其彈道大部分為滑翔段,而滑翔段的彈道特性除受控制律控制外,主要受滑翔段初始狀態(tài)影響。

    文獻[1]針對平衡滑翔彈道對初始狀態(tài)精度要求高的特點,設計了基于預測校正的制導律,增強了平衡滑翔彈道對初始狀態(tài)的魯棒性。文獻[2]利用正則攝動的方法得出了高超聲速滑翔飛行器平穩(wěn)滑翔的解析解,分析了平穩(wěn)滑翔彈道的動態(tài)特性。文獻[3-5]以射程最大、能量最省、總氣動加熱最小等為優(yōu)化目標,分析了不同飛行器初始狀態(tài)以及飛行器升阻比等氣動參數(shù)對彈道的影響,指出飛行器彈道特性主要受飛行器初始狀態(tài)影響。文獻[6-7]利用高斯偽譜法,分別從燃料最省和航路規(guī)劃的角度對高超聲速飛行器的軌跡進行了優(yōu)化。文獻[8]利用物理規(guī)劃和遺傳算法,求解了射程最大、約束變量最小、彈道最穩(wěn)定的優(yōu)化彈道。文獻[9]利用模式搜索算法對高超聲速飛行器軌跡進行優(yōu)化,分析約束條件對飛行器軌跡和彈道的影響,為確定飛行器的初始狀態(tài)提供了參考。文獻[10-11]指出跳躍滑翔彈道的初始狀態(tài)決定了飛行器的飛行能力,用分段高斯偽譜法對多約束條件下的滑翔彈道進行了優(yōu)化。文獻[12]從雷達跟蹤精度的角度對跳躍滑翔彈道的突防性能進行了分析,得出了彈道的機動性能越強、飛行器的突防性能越優(yōu)的結論。文獻[13]利用正交試驗分析了初始狀態(tài)對跳躍滑翔彈道機動性能的影響,通過回歸分析估計了機動性能最優(yōu)時的滑翔段初始狀態(tài)。

    以上文獻的研究主要集中于平衡滑翔彈道的彈道特性分析以及跳躍滑翔彈道的控制律優(yōu)化設計,尚缺乏滑翔段初始狀態(tài)對彈道特性影響規(guī)律的研究。因為高超聲速滑翔飛行器是一種無動力飛行器,滑翔段初始狀態(tài)對飛行器滑翔段的飛行能力和彈道特性至關重要,所以有必要研究滑翔段初始狀態(tài)對高超聲速滑翔飛行器彈道特性的影響規(guī)律。

    1 高超聲速滑翔飛行器滑翔段運動模型

    1.1 飛行器滑翔段運動方程

    根據(jù)高超聲速滑翔飛行器在滑翔段的受力情況,建立速度和位置的四自由度運動方程,如式(1)所示。式中,V是飛行器飛行速度,γ是速度傾角,(x,z)是飛行器在縱向平面內(nèi)的位置坐標,m是飛行器質(zhì)量,g0是重力加速度常數(shù),L是飛行器所受升力,D是飛行器所受阻力。

    (1)

    其中升力、阻力的計算可由式(2)求得。ρ為大氣密度,如式(3)所示,在此采用指數(shù)模型計算飛行器所在飛行高度的大氣密度,ρ0=1.225 kg/m3為海平面大氣密度常數(shù),ξ=8434 m為高度常數(shù);S是飛行器的參考面積;CL和CD為飛行器升力系數(shù)、阻力系數(shù)[14]。

    (2)

    (3)

    1.2 滑翔段彈道控制參數(shù)模型

    1.2.1 縱向平衡滑翔彈道控制參數(shù)模型

    飛行器縱向平衡滑翔時,在縱向平面內(nèi)受力平衡,速度傾角變化率為零(如式(4)所示),又因為飛行器在縱向平面內(nèi)的主要控制參數(shù)是攻角α,所以令飛行器攻角保持不變、同時滿足速度傾角變化率為零的限制條件[8],即可實現(xiàn)臨近空間高超聲速滑翔飛行器的縱向平衡滑翔。

    (4)

    1.2.2 縱向跳躍滑翔彈道控制參數(shù)模型

    控制臨近空間高超聲速滑翔飛行器在縱向平面內(nèi)跳躍運動的主要參數(shù)是攻角α,縱向跳躍滑翔彈道控制參數(shù)模型即攻角α參數(shù)模型[15],如式(5)所示。

    (5)

    式中,αmax為最大飛行攻角,αKmax為最大升阻比攻角,V1和V2分別為攻角參數(shù)變化的兩個臨界速度。

    2 縱向平衡滑翔初始狀態(tài)唯一性分析

    2.1 飛行器狀態(tài)解析解

    由平衡滑翔條件,將式(2)代入式(4)整理可得:

    (6)

    所以可推導出飛行器速度表達式為:

    (7)

    將式(3)代入式(7)整理可得:

    (8)

    式(8)對z求微分可得:

    (9)

    將式(8)代入式(9)可得:

    (10)

    將式(1)代入式(10)可得:

    (11)

    由式(10)和式(11)可得:

    (12)

    將式(2)、式(7)代入式(12)整理可得:

    (13)

    由式(13)可得飛行器速度V的解析式:

    (14)

    由式(7)和式(14)可得:

    (15)

    由式(3)和式(15)可得:

    (16)

    2.2 初始狀態(tài)唯一性分析

    令V=f(γ),對式(14)求導可得:

    (17)

    由式(17)可知f′(γ)>0,所以f(γ)為γ上的單調(diào)遞增函數(shù)。因為f(γ)為單調(diào)遞增函數(shù),所以f-1(V)一定存在,且和f(γ)單調(diào)性一致。

    (18)

    由f-1(V)與f(γ)單調(diào)性一致可得,γ>f-1(Vmin),當取Vmin=1700 m/s時,可得f-1(Vmin)=-2.084°,又因為文獻[15]已經(jīng)證明γ為負的小量,所以可將γ的取值范圍縮小為γ∈(-2.084°,0°)。

    令z=g(γ),C1=ρ0ξSCL,C2=ρ0ξSCD,對式(16)求導并整理可得:

    (19)

    令h(γ)=C1sinγ+C2cosγ可得:

    (20)

    將C1、C2代入式(20)可得:

    (21)

    由文獻[15]可知5°≤α≤25°,所以可得:

    (22)

    因為γ∈(-2.084°,0°),所以可得h(γ)>0。將式(20)代入式(19)可得:

    (23)

    由式(23)可知g′(γ)<0,所以g(γ)為γ上的單調(diào)遞減函數(shù)。又因為f(γ)為γ上的單調(diào)遞增函數(shù),所以對初始狀態(tài)(V0,γ0,z0),只要確定三個初始量中的一個量,其余兩個也唯一確定。當飛行器滑翔段初始狀態(tài)(V0,γ0,z0)滿足解析式(14)和式(16)時,飛行器才會以平衡滑翔彈道滑翔,否則只能以近似平衡滑翔彈道滑翔,或產(chǎn)生大幅度震蕩以跳躍滑翔彈道滑翔。

    2.3 仿真校驗與分析

    2.3.1 仿真設置

    以美國洛克希德·馬丁公司開發(fā)的高超聲速通用氣動飛行器CAV-H為例,對平衡滑翔彈道滑翔段初始狀態(tài)分析進行驗證。CAV-H飛行器質(zhì)量為907.2 kg,參考面積為0.483 9 m2。取目標點位置(2×106m,-100 m),αKmax=12°,x0=0 m。由式(14)和式(16)確定飛行器在45 km高度的平衡滑翔初始狀態(tài)V0=3563 m/s,γ0=-0.247 2°,z0=-45 000 m。

    為分析縱向平衡滑翔時滑翔段初始狀態(tài)中各變量對彈道特性的影響,利用單因素敏感性分析方法對平衡滑翔彈道滑翔段的初始狀態(tài)進行敏感性分析。如表1所示,設計了三組仿真實驗進行對比:①保持(V0,γ0)不變,改變初始高度z0;②保持(γ0,z0)不變,改變初始速度V0;③保持(V0,z0)不變,改變初始速度傾角γ0。

    表1 彈道0和三組對比彈道的滑翔段初始狀態(tài)

    2.3.2 彈道軌跡分析

    圖1是彈道0和三組對比彈道的彈道軌跡仿真。對比圖1(a)、圖1(b)、圖1(c)可知,當初始高度、初始速度傾角偏離初始狀態(tài)真值時,相應彈道都會圍繞平衡滑翔彈道上下小幅度跳躍,跳躍幅度越來越小、逐漸收斂于平衡滑翔彈道,且圖1(c)中三個彈道的波動幅度小于圖1(a)中三個彈道的波動幅度。當初始速度偏離平衡滑翔初始速度真值時,第二組三個彈道在平衡滑翔彈道上方跳躍,且與平衡滑翔彈道的偏差越來越大。從三組對比仿真可知,初始速度偏離平衡滑翔初始速度真值對平衡滑翔彈道軌跡影響最大,初始速度傾角偏離平衡滑翔初始速度傾角真值對平衡滑翔彈道軌跡影響最小。

    (a) 彈道0和第一組彈道的彈道軌跡仿真圖

    2.3.3 飛行速度誤差分析

    圖2是三組對比仿真彈道的飛行速度與彈道0飛行速度的誤差曲線圖。對比圖2(a)、圖2(b)、圖2(c)可知,第二組三個彈道的初始速

    (a) 彈道0和第一組彈道的飛行速度誤差

    度偏離平衡滑翔初始速度真值時,相應彈道的速度誤差最大,同時飛行過程中,速度誤差曲線降低的趨勢并不明顯。第三組三個彈道的初始速度傾角偏離平衡滑翔初始速度傾角真值時,對彈道的速度誤差影響最小,雖然第三組的速度誤差曲線呈現(xiàn)收斂于非零值的趨勢,但其速度誤差數(shù)值在三組對比仿真中最小。

    2.3.4 速度傾角誤差分析

    圖3是三組對比仿真彈道的速度傾角與彈道0速度傾角的誤差曲線圖。對比圖3(a)、圖3(b)、圖3(c)可知,第二組三個彈道的初始速度偏離平衡滑翔初始速度真值時,相應彈道的速度傾角誤差介于第一組和第三組之間,但第二組三個彈道的速度傾角誤差曲線呈現(xiàn)發(fā)散、不斷變大的趨勢。第三組三個彈道的初始速度傾角偏離平衡滑翔初始速度傾角真值時,對彈道的速度傾角誤差影響最小,其速度傾角誤差曲線呈現(xiàn)不斷下降的趨勢,且速度傾角誤差數(shù)值最小。

    (a) 彈道0和第一組彈道的速度傾角誤差

    2.3.5 飛行高度誤差分析

    圖4是三組對比仿真彈道的飛行高度與彈道0飛行高度的誤差曲線圖。對比圖4(a)、圖4(b)、圖4(c)可知,第二組三個彈道的初始速度偏離平衡滑翔初始速度真值時,相應彈道的高度誤差與第一組大致相當,均高于第三組,但第二組三個彈道的高度誤差曲線振蕩發(fā)散、誤差數(shù)值不斷變大。第三組三個彈道的初始速度傾角偏離平衡滑翔初始速度傾角真值時,其高度誤差曲線振蕩收斂,彈道高度誤差數(shù)值在三組對比仿真中最小,對彈道高度誤差影響最小。

    (a) 彈道0和第一組彈道的飛行高度誤差

    3 縱向跳躍滑翔初始狀態(tài)最優(yōu)性分析

    3.1 最優(yōu)性問題描述

    由前文分析可知,飛行器滑翔彈道產(chǎn)生跳躍的原因是初始狀態(tài)不滿足平衡滑翔條件,即跳躍滑翔彈道滑翔段的初始狀態(tài)不唯一。正是因為跳躍滑翔彈道滑翔段初始狀態(tài)不能唯一確定,所以有必要對其初始狀態(tài)的最優(yōu)性進行研究,即尋找跳躍滑翔彈道滑翔段的最優(yōu)初始狀態(tài)。在此以飛行彈道的彈道性能最優(yōu)為目標,初始狀態(tài)最優(yōu)性的數(shù)學模型如下:

    彈道性能:f(V0,γ0,z0)

    約束條件:

    由初始狀態(tài)最優(yōu)性的數(shù)學模型可知,求解最優(yōu)的初始狀態(tài)為多約束條件下的優(yōu)化問題。為便于問題求解,引入懲罰函數(shù),將有約束優(yōu)化問題轉化為無約束優(yōu)化問題,如式(24)所示,為彈道性能函數(shù)f(V0,γ0,z0)引入懲罰因子pi(i=1,2,…,5)。當滑翔彈道不滿足約束條件時,需要在相應彈道的性能函數(shù)上疊加懲罰因子。

    (24)

    3.2 評價指標

    為了評價以不同初始狀態(tài)進行跳躍滑翔時的彈道特性,特選擇以下幾個評價指標。

    1)機動變軌次數(shù):機動變軌次數(shù)Nm定義為飛行器縱向加速度符號發(fā)生變化的次數(shù)。高超聲速滑翔飛行器的彈道上下跳躍,所以無法像彈道導彈一樣反推導彈的發(fā)射點和預測導彈的落點。機動變軌次數(shù)是高超聲速滑翔飛行器機動性能的重要體現(xiàn),機動變軌次數(shù)越多,飛行器彈道越靈活多變。

    2)首跳高度:首跳高度ΔH是指飛行器第一次跳躍時的高度大小,其數(shù)值等于滑翔段初始位置與彈道軌跡曲線第一個極值點(極大值或極小值)之間的高度差。因為無法準確定義跳躍滑翔彈道的平衡位置,且隨著飛行速度的下降,飛行器的跳躍幅度越來越小,所以采用首跳高度來衡量高超聲速滑翔飛行器彈道的跳躍強度。

    3)飛行時間:飛行時間ΔT是指飛行器滑翔段初始時刻與滑翔段結束時刻的差值。飛行時間越短,飛行器暴露在雷達、衛(wèi)星等探測裝備探測區(qū)的時間越短,越不容易被發(fā)現(xiàn)。即便被探測裝備發(fā)現(xiàn),因為飛行時間短,留給攔截武器的時間將被大大壓縮,也將極大提高飛行器的生存能力。

    4)剩余機械能:剩余機械能Ef是指滑翔段結束時刻高超聲速滑翔飛行器總的機械能。高超聲速滑翔飛行器在滑翔段做無動力滑翔,所以飛行過程中其機械能是不斷減小的?;瓒问S鄼C械能即飛行器末段的最大機械能。飛行器的剩余機械能越大,則為末段彈道提供的可利用能量空間越大。

    綜合以上評價指標以及約束條件,可確定求解最優(yōu)初始狀態(tài)時的目標函數(shù),如式(25)所示。其中ω1、ω2、ω3、ω4分別為四個評價指標的指標權重,可通過本征向量法確定各評價指標的指標權重值為(0.472 9,0.169 9,0.284 4,0.072 8)。

    f(V0,γ0,z0)=max(ω1Nm+ω2ΔH+

    (25)

    3.3 粒子群算法

    粒子群算法是一種全局迭代搜索算法,該算法將優(yōu)化問題的解看作是搜索空間中的粒子。算法起始時,隨機產(chǎn)生種群大小為Npso的初始粒子群,即初始解集。在每次迭代中,利用目標函數(shù)計算粒子的適應度值,并更新粒子個體最優(yōu)解和粒子群全局最優(yōu)解。兩個最優(yōu)解更新完畢后,按照式(26)更新各粒子的速度和位置,實現(xiàn)粒子在解空間的局部搜索和全局搜索。當算法達到終止條件時,迭代結束并輸出搜索結果。

    (26)

    式中,vj(t)為第t次迭代后粒子的速度,wpso為搜索慣性權重,C1pso和C2pso為學習因子,r1(t)、r2(t)為介于0和1之間的隨機數(shù),pj(t)為第t次迭代后粒子的個體最優(yōu)解,g(t)為第t次迭代后粒子群的全局最優(yōu)解,zj(t)為第t次迭代后粒子的位置。

    3.4 仿真校驗與分析

    3.4.1 仿真校驗

    通過粒子群算法可求得最優(yōu)的初始狀態(tài)(V0,γ0,z0)為(4864.718 m/s,-4.727 8°,-93 691.965 m),機動變軌次數(shù)Nm=12,首跳高度ΔH=53.181 km,飛行時間ΔT=1085 s,剩余機械能Ef=4.711×108J,目標函數(shù)值f(V0,γ0,z0)=0.568 3。

    為分析縱向跳躍滑翔彈道滑翔段最優(yōu)初始狀態(tài)中各變量對彈道特性的影響,利用單因素敏感性分析方法對跳躍滑翔彈道滑翔段的最優(yōu)初始狀態(tài)進行敏感性分析。如表2所示,設計以下三組對比實驗:①在最優(yōu)的初始速度和速度傾角(V0,γ0)條件下,改變初始高度z0;②在最優(yōu)的初始速度傾角和高度(γ0,z0)條件下,改變初始速度V0;③在最優(yōu)的初始速度和高度(V0,z0)條件下,改變初始速度傾角γ0。

    表2 彈道10和三組對比彈道的滑翔段初始狀態(tài)

    三組對比實驗以及最優(yōu)初始狀態(tài)對應彈道的性能指標值和目標函數(shù)值如表3所示。由表3可知,最優(yōu)初始狀態(tài)的目標函數(shù)值最大,第五組三個彈道的目標函數(shù)值最小。下面從彈道軌跡和縱向加速度兩個方面對三組彈道和最優(yōu)初始狀態(tài)對應彈道的彈道特性進行分析。

    表3 彈道10和三組對比彈道的評價指標值和目標函數(shù)值

    3.4.2 彈道軌跡分析

    圖5是三組對比彈道和彈道10的彈道軌跡仿真圖。從圖5(a)中可以看出,當初始高度偏離最優(yōu)初始高度時,彈道的首跳高度會相應降低,且偏離程度越大,首跳高度越小。從圖5(b)中可以看出第五組的三個彈道,滑翔段終點距離目標點的終端距離大于200 km,不滿足終端約束。雖然第五組三個彈道的性能指標,除機動變軌次數(shù)外均優(yōu)于彈道10的性能指標,但因為第五組的三個彈道不滿足終端約束條件,所以第五組三個彈道初始狀態(tài)的目標函數(shù)值明顯低于其他兩組彈道和彈道10初始狀態(tài)的目標函數(shù)值。從圖5(c)中可以看出,當初始速度傾角偏離最優(yōu)初始速度傾角時,彈道的首跳高度會相應降低,且偏離程度越大,首跳高度越小。

    對比圖5(a)、圖5(b)、圖5(c)可知,當彈道的初始狀態(tài)偏離最優(yōu)初始狀態(tài)時,彈道軌跡均會偏離最優(yōu)彈道軌跡,相應初始狀態(tài)的目標函數(shù)值都低于最優(yōu)初始狀態(tài)的目標函數(shù)值。當初始速度偏離最優(yōu)初始速度時,相應彈道軌跡偏離最優(yōu)彈道的程度越來越大,且不能滿足終端約束條件,對彈道性能影響最大。當初始速度傾角偏離最優(yōu)初始速度傾角時,彈道軌跡偏離最優(yōu)彈道的程度最低,首跳高度下降在幾十米的量級。當初始高度偏離最優(yōu)初始高度時,彈道軌跡偏離最優(yōu)彈道的程度介于第六組和第五組之間,其首跳高度下降在幾百米的量級。

    (a) 彈道10和第四組彈道的彈道軌跡仿真圖

    3.4.3 縱向加速度分析

    圖6是兩組對比彈道和彈道10的縱向加速度曲線,因為第五組的三個彈道不滿足終端約束條件,所以主要對第四組和第六組彈道的縱向加速度進行分析。從圖6(a)中可以看出,當初始高度偏離最優(yōu)初始高度時,縱向加速度曲線會偏離最優(yōu)初始狀態(tài)對應彈道的縱向加速度曲線。在飛行結束之前,彈道11、彈道12、彈道13均能完成最后一次縱向加速度符號翻轉,所以彈道11、彈道12和彈道13的縱向機動變軌次數(shù)與彈道10相同。從圖6(b)中可以看出,當初始速度傾角偏離最優(yōu)初始速度傾角時,縱向加速度曲線會偏離最優(yōu)初始狀態(tài)對應彈道的縱向加速度曲線。在飛行結束之前,彈道17、彈道18、彈道19均未能完成最后一次縱向加速度符號翻轉,所以彈道17、彈道18和彈道19的縱向機動變軌次數(shù)少于彈道10的縱向機動變軌次數(shù)。

    (a) 彈道10和第四組彈道的縱向加速度

    4 結論

    針對臨近空間高超聲速滑翔飛行器滑翔段初始狀態(tài)影響其彈道特性的規(guī)律問題,在平衡滑翔條件下,從理論上證明了縱向平衡滑翔彈道滑翔段初始狀態(tài)的唯一性。由敏感性分析可知,當初始速度偏離平衡滑翔初始速度真值時,隨著飛行過程的進行,其彈道逐漸發(fā)散,對彈道特性的影響最大;當初始速度傾角偏離平衡滑翔初始速度傾角真值時,其彈道軌跡、狀態(tài)變量均趨向于收斂到平衡滑翔的彈道軌跡和狀態(tài)變量,對彈道特性的影響最小。在跳躍滑翔條件下,利用粒子群算法求解了彈道性能最優(yōu)時的滑翔段初始狀態(tài)。由敏感性分析可知,當初始速度偏離跳躍滑翔最優(yōu)初始速度時,其彈道不滿足約束條件,相應彈道初始狀態(tài)的目標函數(shù)值最低;當初始速度傾角偏離最優(yōu)初始速度傾角時,從彈道軌跡的角度分析,其彈道軌跡最接近最優(yōu)的彈道軌跡。

    猜你喜歡
    滑翔超聲速彈道
    蘭花螳螂會滑翔
    彈道——打勝仗的奧秘
    高超聲速出版工程
    高超聲速飛行器
    攻天掠地的先鋒武器——滑翔導彈
    一維彈道修正彈無線通信系統(tǒng)研制
    電子制作(2019年7期)2019-04-25 13:17:48
    超聲速旅行
    讓你心跳加速的“滑翔飛板”
    小學時代(2017年16期)2017-06-19 19:33:19
    基于PID控制的二維彈道修正彈仿真
    制導與引信(2016年3期)2016-03-20 16:02:02
    空中滑翔大比拼(下)——滑翔傘
    色哟哟哟哟哟哟| 国产精品无大码| 国产精品人妻久久久久久| 久久久久性生活片| 日本一二三区视频观看| 亚州av有码| 亚洲美女黄片视频| 22中文网久久字幕| 一级黄色大片毛片| 成人永久免费在线观看视频| 一级毛片aaaaaa免费看小| 国产久久久一区二区三区| 亚洲国产色片| 卡戴珊不雅视频在线播放| 超碰av人人做人人爽久久| 亚洲国产欧洲综合997久久,| 亚洲国产欧洲综合997久久,| 日韩欧美三级三区| 国内揄拍国产精品人妻在线| 精品人妻偷拍中文字幕| 精品一区二区三区人妻视频| 国产视频一区二区在线看| 日韩,欧美,国产一区二区三区 | .国产精品久久| 欧美激情久久久久久爽电影| 日韩欧美在线乱码| 亚洲av成人精品一区久久| 久久精品国产99精品国产亚洲性色| 99热这里只有是精品在线观看| 精品欧美国产一区二区三| 成人无遮挡网站| 欧美成人一区二区免费高清观看| 两个人视频免费观看高清| 日韩欧美免费精品| 亚洲性夜色夜夜综合| a级毛片免费高清观看在线播放| 精品午夜福利在线看| 99国产极品粉嫩在线观看| 22中文网久久字幕| 亚洲激情五月婷婷啪啪| 国内精品美女久久久久久| 国产免费一级a男人的天堂| 熟女电影av网| 三级经典国产精品| 国产精品一区二区三区四区免费观看 | 精品一区二区三区视频在线| 亚洲国产色片| 国产人妻一区二区三区在| 69av精品久久久久久| 岛国在线免费视频观看| 日韩强制内射视频| 一区福利在线观看| 美女免费视频网站| 亚洲精品一区av在线观看| av天堂在线播放| 色在线成人网| 亚洲欧美日韩无卡精品| 人人妻,人人澡人人爽秒播| 人妻夜夜爽99麻豆av| 蜜臀久久99精品久久宅男| 婷婷色综合大香蕉| 一进一出好大好爽视频| 日韩精品有码人妻一区| 男人舔女人下体高潮全视频| 成人三级黄色视频| av国产免费在线观看| 村上凉子中文字幕在线| ponron亚洲| 国产午夜福利久久久久久| 成人特级黄色片久久久久久久| 亚洲精品一卡2卡三卡4卡5卡| 搞女人的毛片| 伊人久久精品亚洲午夜| 色5月婷婷丁香| 可以在线观看毛片的网站| 99久久精品一区二区三区| 成人永久免费在线观看视频| 看黄色毛片网站| 看黄色毛片网站| 国产精品99久久久久久久久| 不卡视频在线观看欧美| 黄色配什么色好看| 在线天堂最新版资源| 亚洲欧美精品自产自拍| 亚洲国产高清在线一区二区三| 不卡视频在线观看欧美| 真人做人爱边吃奶动态| 麻豆一二三区av精品| 日韩成人伦理影院| 国产亚洲精品av在线| 国产伦在线观看视频一区| 国产精品乱码一区二三区的特点| 三级经典国产精品| 久久精品国产亚洲av天美| 麻豆乱淫一区二区| 午夜精品国产一区二区电影 | 免费av不卡在线播放| 成人特级黄色片久久久久久久| 日日啪夜夜撸| 亚洲成人久久性| .国产精品久久| 日本-黄色视频高清免费观看| 亚洲成人久久性| av天堂在线播放| 美女免费视频网站| 日韩制服骚丝袜av| 波多野结衣高清作品| 欧美一区二区国产精品久久精品| 久久久精品94久久精品| 嫩草影院精品99| 国产黄色小视频在线观看| 日韩强制内射视频| 一夜夜www| 久久午夜亚洲精品久久| a级毛色黄片| 熟女人妻精品中文字幕| 人妻少妇偷人精品九色| 三级毛片av免费| 成年免费大片在线观看| 中文字幕人妻熟人妻熟丝袜美| 精品国内亚洲2022精品成人| 日日摸夜夜添夜夜添av毛片| 少妇高潮的动态图| 成人亚洲欧美一区二区av| 国产一级毛片七仙女欲春2| 日本与韩国留学比较| 欧美日韩精品成人综合77777| 久久亚洲精品不卡| 在线免费观看的www视频| 国产精品三级大全| 国产精品女同一区二区软件| 一本久久中文字幕| 免费搜索国产男女视频| 中国美女看黄片| 美女 人体艺术 gogo| 人妻久久中文字幕网| 成人美女网站在线观看视频| 人人妻人人看人人澡| 不卡视频在线观看欧美| av福利片在线观看| 亚洲av中文av极速乱| 午夜a级毛片| 国产精品三级大全| 女人十人毛片免费观看3o分钟| 99riav亚洲国产免费| 久久精品夜色国产| 内地一区二区视频在线| 欧美另类亚洲清纯唯美| 一级av片app| 欧洲精品卡2卡3卡4卡5卡区| 国产在线精品亚洲第一网站| 赤兔流量卡办理| 白带黄色成豆腐渣| 嫩草影视91久久| 久久精品久久久久久噜噜老黄 | 简卡轻食公司| 国产麻豆成人av免费视频| 亚洲美女搞黄在线观看 | 无遮挡黄片免费观看| 婷婷六月久久综合丁香| 淫秽高清视频在线观看| 99九九线精品视频在线观看视频| 国产一区亚洲一区在线观看| 97在线视频观看| 亚洲aⅴ乱码一区二区在线播放| 国产亚洲av嫩草精品影院| 国产中年淑女户外野战色| 亚洲国产精品久久男人天堂| 国产亚洲精品综合一区在线观看| h日本视频在线播放| 中文字幕av成人在线电影| 搡老岳熟女国产| 九九热线精品视视频播放| 97热精品久久久久久| 精品无人区乱码1区二区| 精品99又大又爽又粗少妇毛片| 青春草视频在线免费观看| 黄色配什么色好看| 午夜福利在线观看免费完整高清在 | 久久久久国内视频| 麻豆成人午夜福利视频| 中文字幕免费在线视频6| 看黄色毛片网站| 亚洲成人av在线免费| av天堂在线播放| 亚州av有码| 我要搜黄色片| 麻豆国产av国片精品| 成人二区视频| 亚洲欧美日韩高清在线视频| 亚洲av二区三区四区| 亚洲激情五月婷婷啪啪| 国产精品国产三级国产av玫瑰| 国产黄色小视频在线观看| 熟女人妻精品中文字幕| 亚洲精品影视一区二区三区av| 18禁在线无遮挡免费观看视频 | 精品久久久久久久久久久久久| 十八禁国产超污无遮挡网站| 久久鲁丝午夜福利片| 老熟妇仑乱视频hdxx| 久久精品国产自在天天线| 露出奶头的视频| 欧美xxxx黑人xx丫x性爽| 别揉我奶头~嗯~啊~动态视频| 日本三级黄在线观看| 亚洲精华国产精华液的使用体验 | 三级经典国产精品| 91午夜精品亚洲一区二区三区| 亚洲综合色惰| 日本爱情动作片www.在线观看 | 国产极品精品免费视频能看的| 欧美日韩在线观看h| 日韩三级伦理在线观看| 日韩成人av中文字幕在线观看 | 亚洲三级黄色毛片| 国产黄片美女视频| 91av网一区二区| 久久久久久大精品| 美女cb高潮喷水在线观看| 啦啦啦啦在线视频资源| 久久久久久久亚洲中文字幕| 老熟妇乱子伦视频在线观看| 国产精品久久久久久久电影| 国内揄拍国产精品人妻在线| 俄罗斯特黄特色一大片| 午夜爱爱视频在线播放| 久久久久久久亚洲中文字幕| av福利片在线观看| 村上凉子中文字幕在线| 成人特级av手机在线观看| 久久久久国内视频| 性欧美人与动物交配| 国产在线男女| 日本-黄色视频高清免费观看| 精品久久久久久久久亚洲| av在线天堂中文字幕| 免费搜索国产男女视频| 夜夜夜夜夜久久久久| 国产激情偷乱视频一区二区| 国产亚洲精品久久久久久毛片| 热99在线观看视频| 日日干狠狠操夜夜爽| 国产视频内射| 精品久久久久久久久久免费视频| 熟女电影av网| 99久久中文字幕三级久久日本| 天天躁日日操中文字幕| 欧美成人免费av一区二区三区| 特级一级黄色大片| 婷婷精品国产亚洲av在线| 国产欧美日韩精品亚洲av| 乱人视频在线观看| 国产精品美女特级片免费视频播放器| 亚洲在线观看片| 亚洲经典国产精华液单| 国产一区二区在线观看日韩| 亚洲精品国产成人久久av| 简卡轻食公司| 最近最新中文字幕大全电影3| 免费无遮挡裸体视频| 日本与韩国留学比较| 我要看日韩黄色一级片| 国产麻豆成人av免费视频| 国产探花在线观看一区二区| 九九在线视频观看精品| 男女下面进入的视频免费午夜| 午夜亚洲福利在线播放| 国产淫片久久久久久久久| 99久国产av精品国产电影| 男人狂女人下面高潮的视频| 国产av一区在线观看免费| 欧美高清性xxxxhd video| 国产精品野战在线观看| 我的老师免费观看完整版| 日本精品一区二区三区蜜桃| 一边摸一边抽搐一进一小说| 中文字幕免费在线视频6| 国产精品99久久久久久久久| 日韩av不卡免费在线播放| 亚洲精品粉嫩美女一区| 日韩亚洲欧美综合| 亚洲不卡免费看| 免费看光身美女| 亚洲精品日韩av片在线观看| 一级毛片久久久久久久久女| 精品无人区乱码1区二区| 少妇的逼水好多| 此物有八面人人有两片| 九色成人免费人妻av| 一本一本综合久久| 欧美国产日韩亚洲一区| 亚洲国产日韩欧美精品在线观看| 日本与韩国留学比较| 黄色视频,在线免费观看| 非洲黑人性xxxx精品又粗又长| 精品国内亚洲2022精品成人| 人妻夜夜爽99麻豆av| 精品午夜福利在线看| 人妻少妇偷人精品九色| 99精品在免费线老司机午夜| 一进一出抽搐gif免费好疼| 黑人高潮一二区| 国内精品宾馆在线| 淫妇啪啪啪对白视频| 女同久久另类99精品国产91| 国产蜜桃级精品一区二区三区| 丝袜美腿在线中文| 成熟少妇高潮喷水视频| 网址你懂的国产日韩在线| 国产高清激情床上av| 美女内射精品一级片tv| 女人被狂操c到高潮| 一个人看视频在线观看www免费| 日本与韩国留学比较| 国产亚洲精品综合一区在线观看| 成年女人看的毛片在线观看| 日韩欧美三级三区| 麻豆成人午夜福利视频| 精品国内亚洲2022精品成人| 一区二区三区高清视频在线| 亚洲国产欧美人成| 三级毛片av免费| 亚洲国产日韩欧美精品在线观看| 人人妻人人澡欧美一区二区| 亚洲最大成人av| 韩国av在线不卡| 内射极品少妇av片p| 一级a爱片免费观看的视频| 欧美性猛交黑人性爽| 少妇人妻精品综合一区二区 | 日本精品一区二区三区蜜桃| 国产亚洲av嫩草精品影院| 联通29元200g的流量卡| 国产三级中文精品| 深爱激情五月婷婷| 一级黄片播放器| 国产伦一二天堂av在线观看| 国产精品无大码| 69人妻影院| 蜜桃久久精品国产亚洲av| 日本成人三级电影网站| 国产亚洲91精品色在线| 国产亚洲av嫩草精品影院| 国产又黄又爽又无遮挡在线| 亚洲一区二区三区色噜噜| 亚洲av成人av| 狂野欧美白嫩少妇大欣赏| 中文资源天堂在线| 精品人妻视频免费看| 久久综合国产亚洲精品| 欧美极品一区二区三区四区| 欧美最新免费一区二区三区| 日韩制服骚丝袜av| 男人的好看免费观看在线视频| 一级毛片aaaaaa免费看小| 成年女人看的毛片在线观看| 亚洲最大成人av| 性欧美人与动物交配| 亚洲欧美日韩东京热| 久久久成人免费电影| 男人和女人高潮做爰伦理| 成人一区二区视频在线观看| av中文乱码字幕在线| 精品免费久久久久久久清纯| 国内久久婷婷六月综合欲色啪| av.在线天堂| 国产一区二区三区在线臀色熟女| 色av中文字幕| 在线观看66精品国产| 国产亚洲91精品色在线| 午夜爱爱视频在线播放| 日韩欧美在线乱码| 一级a爱片免费观看的视频| 一边摸一边抽搐一进一小说| av在线老鸭窝| 中文字幕av成人在线电影| 午夜亚洲福利在线播放| 久久精品国产鲁丝片午夜精品| 精品熟女少妇av免费看| 久久久久久久午夜电影| 老司机福利观看| 日本黄色视频三级网站网址| 99在线人妻在线中文字幕| 国产日本99.免费观看| 日韩精品中文字幕看吧| 美女黄网站色视频| 能在线免费观看的黄片| 尾随美女入室| 一区福利在线观看| 春色校园在线视频观看| 男女那种视频在线观看| 97碰自拍视频| 中文字幕熟女人妻在线| 老司机福利观看| 日本成人三级电影网站| 国产黄a三级三级三级人| 女人被狂操c到高潮| 国产精品久久久久久久电影| 国产精品爽爽va在线观看网站| 99久国产av精品| av.在线天堂| or卡值多少钱| 亚洲精品乱码久久久v下载方式| 成人毛片a级毛片在线播放| 超碰av人人做人人爽久久| 日韩欧美国产在线观看| 亚洲精品在线观看二区| 久久韩国三级中文字幕| a级一级毛片免费在线观看| 欧美日本亚洲视频在线播放| 男女之事视频高清在线观看| 又爽又黄a免费视频| 22中文网久久字幕| 精品免费久久久久久久清纯| 久久人妻av系列| 日本黄色视频三级网站网址| 有码 亚洲区| 女人被狂操c到高潮| 九九在线视频观看精品| 哪里可以看免费的av片| 成年免费大片在线观看| 午夜福利在线观看免费完整高清在 | 国语自产精品视频在线第100页| 精品午夜福利在线看| 成人鲁丝片一二三区免费| 插逼视频在线观看| 丰满乱子伦码专区| 久久久久精品国产欧美久久久| 成人av一区二区三区在线看| 女人十人毛片免费观看3o分钟| 国产激情偷乱视频一区二区| 少妇丰满av| 久99久视频精品免费| 黄色日韩在线| 国产一区二区三区在线臀色熟女| 国产精品久久久久久久久免| 婷婷精品国产亚洲av| 中国国产av一级| 国产高清三级在线| 午夜影院日韩av| 日本欧美国产在线视频| 搡老熟女国产l中国老女人| 国产精品久久久久久精品电影| a级毛色黄片| 亚洲电影在线观看av| 天美传媒精品一区二区| 日韩大尺度精品在线看网址| 97人妻精品一区二区三区麻豆| 久久国内精品自在自线图片| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲自偷自拍三级| 极品教师在线视频| 中文字幕av成人在线电影| 热99re8久久精品国产| 久久99热6这里只有精品| 亚洲第一区二区三区不卡| 亚洲av熟女| 在线播放无遮挡| 深夜精品福利| 男人狂女人下面高潮的视频| 婷婷精品国产亚洲av| 成人漫画全彩无遮挡| 日本一本二区三区精品| 看非洲黑人一级黄片| 亚洲四区av| 久久鲁丝午夜福利片| 狂野欧美激情性xxxx在线观看| 天堂影院成人在线观看| 国产 一区精品| 1024手机看黄色片| 在线a可以看的网站| 欧美日韩乱码在线| 精品久久久久久久人妻蜜臀av| 国产蜜桃级精品一区二区三区| 国产精品一二三区在线看| 欧美三级亚洲精品| 国产精品久久久久久av不卡| 伦精品一区二区三区| 国产国拍精品亚洲av在线观看| 99久久成人亚洲精品观看| 国产男靠女视频免费网站| 亚洲国产精品成人久久小说 | 看免费成人av毛片| 成人特级av手机在线观看| 成人美女网站在线观看视频| 国产不卡一卡二| 精品无人区乱码1区二区| 夜夜看夜夜爽夜夜摸| 一夜夜www| 亚洲在线观看片| 亚洲美女视频黄频| 国产黄色视频一区二区在线观看 | 国产精品久久久久久av不卡| 免费av观看视频| 亚洲18禁久久av| 别揉我奶头 嗯啊视频| 午夜影院日韩av| 精品一区二区三区视频在线| 国产精品1区2区在线观看.| 国产精品女同一区二区软件| 黄色欧美视频在线观看| 看非洲黑人一级黄片| 欧美色欧美亚洲另类二区| 国产一级毛片七仙女欲春2| 欧美成人a在线观看| 国产伦在线观看视频一区| 久久久久久久午夜电影| 一本精品99久久精品77| av国产免费在线观看| 在线免费观看不下载黄p国产| 国产乱人视频| 日韩 亚洲 欧美在线| 12—13女人毛片做爰片一| 国产亚洲av嫩草精品影院| 寂寞人妻少妇视频99o| 亚洲精品影视一区二区三区av| 国产免费男女视频| 久久人人爽人人爽人人片va| 色尼玛亚洲综合影院| 神马国产精品三级电影在线观看| 99久国产av精品国产电影| 亚洲av二区三区四区| 91在线精品国自产拍蜜月| 日韩欧美在线乱码| 国产国拍精品亚洲av在线观看| 中文在线观看免费www的网站| 久久亚洲精品不卡| 免费看av在线观看网站| 91久久精品电影网| 国内精品宾馆在线| 国产老妇女一区| av视频在线观看入口| 亚洲,欧美,日韩| 欧美一区二区国产精品久久精品| 又黄又爽又刺激的免费视频.| 亚洲欧美成人综合另类久久久 | 18禁黄网站禁片免费观看直播| 一级黄片播放器| 成人毛片a级毛片在线播放| 国产亚洲精品久久久com| 99视频精品全部免费 在线| 狠狠狠狠99中文字幕| 欧美人与善性xxx| 亚洲无线在线观看| 黄色欧美视频在线观看| 久久久久久久久大av| 久久久久精品国产欧美久久久| 99久久中文字幕三级久久日本| 最近最新中文字幕大全电影3| 精品午夜福利在线看| 成人一区二区视频在线观看| 精品久久久久久久久av| 青春草视频在线免费观看| 免费不卡的大黄色大毛片视频在线观看 | 日本a在线网址| 亚洲性久久影院| 日韩强制内射视频| 午夜免费男女啪啪视频观看 | 成人美女网站在线观看视频| 一进一出好大好爽视频| 国产精品永久免费网站| 婷婷六月久久综合丁香| 欧洲精品卡2卡3卡4卡5卡区| 男女做爰动态图高潮gif福利片| 别揉我奶头 嗯啊视频| 久久久久久伊人网av| 老师上课跳d突然被开到最大视频| 久久久久久久久大av| 国产精品永久免费网站| 免费av毛片视频| 又黄又爽又刺激的免费视频.| 亚洲国产欧洲综合997久久,| 欧美最黄视频在线播放免费| 国产高清视频在线观看网站| 国产aⅴ精品一区二区三区波| 日韩制服骚丝袜av| 亚洲第一电影网av| 亚洲av.av天堂| 男女那种视频在线观看| 国产精品一区二区三区四区久久| 国产色婷婷99| 男女边吃奶边做爰视频| 99久久精品一区二区三区| 国产精品三级大全| 97超碰精品成人国产| 超碰av人人做人人爽久久| 少妇熟女aⅴ在线视频| av在线播放精品| 波野结衣二区三区在线| 午夜老司机福利剧场| 小说图片视频综合网站| 国产亚洲av嫩草精品影院| 国内精品美女久久久久久| 成人av在线播放网站| 亚洲精品亚洲一区二区| 国产中年淑女户外野战色| 午夜影院日韩av| 特大巨黑吊av在线直播| avwww免费| 国内久久婷婷六月综合欲色啪| av免费在线看不卡| 观看免费一级毛片| 国产欧美日韩精品亚洲av| 久久人妻av系列| 长腿黑丝高跟| 大又大粗又爽又黄少妇毛片口| 91久久精品国产一区二区三区| 欧美高清成人免费视频www| 乱人视频在线观看| 欧美激情久久久久久爽电影| 天堂√8在线中文| 黄色欧美视频在线观看| a级毛片a级免费在线| 中文字幕av在线有码专区|