• <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
    空中滑翔大比拼(下)——滑翔傘
    在线观看免费视频网站a站| 久久午夜亚洲精品久久| 亚洲国产欧美在线一区| 91九色精品人成在线观看| 日韩一区二区三区影片| h视频一区二区三区| 精品国产乱码久久久久久男人| 欧美精品一区二区免费开放| 精品少妇久久久久久888优播| 少妇 在线观看| 亚洲成人免费av在线播放| 成年人免费黄色播放视频| 久久久久久久精品吃奶| 可以免费在线观看a视频的电影网站| 午夜免费鲁丝| 精品国产超薄肉色丝袜足j| 在线观看免费视频网站a站| 18禁美女被吸乳视频| 精品一区二区三区四区五区乱码| 色播在线永久视频| 国产在线一区二区三区精| 亚洲第一青青草原| 久久午夜亚洲精品久久| 中文字幕av电影在线播放| 亚洲精品在线美女| 日韩欧美三级三区| 成年人黄色毛片网站| 成人影院久久| 不卡一级毛片| 精品人妻熟女毛片av久久网站| 又大又爽又粗| 免费av中文字幕在线| 在线十欧美十亚洲十日本专区| 亚洲少妇的诱惑av| 狠狠婷婷综合久久久久久88av| 99热国产这里只有精品6| 精品一区二区三区av网在线观看 | 亚洲国产看品久久| av网站免费在线观看视频| 老司机福利观看| 国产成+人综合+亚洲专区| 自线自在国产av| 大香蕉久久网| 两人在一起打扑克的视频| 9191精品国产免费久久| 新久久久久国产一级毛片| 久久久国产精品麻豆| av福利片在线| 精品国产国语对白av| 免费在线观看影片大全网站| 国产精品二区激情视频| 日本一区二区免费在线视频| 久久国产亚洲av麻豆专区| 少妇粗大呻吟视频| 一个人免费在线观看的高清视频| 亚洲精品乱久久久久久| 9热在线视频观看99| 国产亚洲精品第一综合不卡| 亚洲av电影在线进入| 女警被强在线播放| 午夜福利影视在线免费观看| 国产淫语在线视频| 欧美精品一区二区免费开放| 看免费av毛片| 精品亚洲成国产av| 欧美在线一区亚洲| 亚洲欧洲精品一区二区精品久久久| 亚洲五月色婷婷综合| 国产又色又爽无遮挡免费看| 黄片大片在线免费观看| 亚洲 欧美一区二区三区| 国产精品久久久久久精品古装| 十八禁网站网址无遮挡| 亚洲成a人片在线一区二区| 99九九在线精品视频| 日韩欧美国产一区二区入口| 免费女性裸体啪啪无遮挡网站| 国产成人免费观看mmmm| 夜夜骑夜夜射夜夜干| 久久久国产欧美日韩av| 久久精品91无色码中文字幕| 怎么达到女性高潮| 国产日韩欧美亚洲二区| 最近最新免费中文字幕在线| 久久热在线av| 亚洲av成人不卡在线观看播放网| 国产av精品麻豆| 日韩大片免费观看网站| 亚洲午夜精品一区,二区,三区| 精品一区二区三卡| 91国产中文字幕| 黄色a级毛片大全视频| 亚洲精品国产区一区二| 精品久久久久久久毛片微露脸| 女人精品久久久久毛片| 亚洲五月色婷婷综合| 欧美激情 高清一区二区三区| 久久性视频一级片| 啦啦啦视频在线资源免费观看| 男人舔女人的私密视频| 亚洲欧美一区二区三区黑人| 人妻久久中文字幕网| 一级a爱视频在线免费观看| 欧美午夜高清在线| 久久香蕉激情| 夜夜骑夜夜射夜夜干| 久久狼人影院| 国产av又大| 亚洲精品av麻豆狂野| 在线亚洲精品国产二区图片欧美| 一夜夜www| av网站在线播放免费| 又紧又爽又黄一区二区| 99精品欧美一区二区三区四区| 在线观看免费午夜福利视频| 99九九在线精品视频| 天天躁狠狠躁夜夜躁狠狠躁| 又大又爽又粗| 午夜91福利影院| 少妇精品久久久久久久| a在线观看视频网站| 国产片内射在线| 在线观看www视频免费| 黄频高清免费视频| 最新美女视频免费是黄的| 在线观看www视频免费| 国产在线免费精品| 国产激情久久老熟女| 中文亚洲av片在线观看爽 | 久9热在线精品视频| 在线观看免费高清a一片| 国产野战对白在线观看| 大香蕉久久成人网| 午夜91福利影院| 精品视频人人做人人爽| 69精品国产乱码久久久| 老司机深夜福利视频在线观看| 亚洲国产毛片av蜜桃av| 亚洲人成电影观看| 美国免费a级毛片| 精品视频人人做人人爽| 国产欧美日韩综合在线一区二区| 色老头精品视频在线观看| 亚洲第一青青草原| 亚洲精品乱久久久久久| 日本黄色视频三级网站网址 | 色精品久久人妻99蜜桃| 两个人看的免费小视频| 自线自在国产av| 18禁美女被吸乳视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品中文字幕在线视频| 亚洲七黄色美女视频| av一本久久久久| 久久久久久人人人人人| 亚洲av第一区精品v没综合| 国产精品 欧美亚洲| 久久婷婷成人综合色麻豆| 亚洲中文av在线| 亚洲成a人片在线一区二区| 欧美日韩精品网址| 免费高清在线观看日韩| 国产精品一区二区在线不卡| 国产日韩欧美亚洲二区| 亚洲精品自拍成人| 9191精品国产免费久久| 免费在线观看完整版高清| 日韩有码中文字幕| 久久国产精品影院| 精品国产一区二区三区久久久樱花| 美女扒开内裤让男人捅视频| 婷婷成人精品国产| 久久99一区二区三区| 91大片在线观看| 99久久人妻综合| 久久精品国产亚洲av高清一级| 久久影院123| 美女视频免费永久观看网站| 大型av网站在线播放| 欧美日韩一级在线毛片| 免费不卡黄色视频| 亚洲专区国产一区二区| 90打野战视频偷拍视频| 精品第一国产精品| 丰满迷人的少妇在线观看| 国产无遮挡羞羞视频在线观看| 国产三级黄色录像| 中文字幕最新亚洲高清| 亚洲全国av大片| 在线看a的网站| 我的亚洲天堂| 人妻 亚洲 视频| 日本一区二区免费在线视频| 麻豆乱淫一区二区| 免费观看人在逋| 国产日韩一区二区三区精品不卡| 亚洲少妇的诱惑av| 久久久国产成人免费| 欧美国产精品va在线观看不卡| 1024香蕉在线观看| 成人三级做爰电影| 十八禁高潮呻吟视频| 久久精品91无色码中文字幕| 夜夜骑夜夜射夜夜干| 欧美黄色片欧美黄色片| 欧美中文综合在线视频| 精品亚洲成国产av| 女警被强在线播放| 动漫黄色视频在线观看| 老鸭窝网址在线观看| 亚洲色图综合在线观看| 国产91精品成人一区二区三区 | 精品亚洲成国产av| 操出白浆在线播放| 午夜两性在线视频| 日韩视频一区二区在线观看| av天堂久久9| 黑人巨大精品欧美一区二区蜜桃| 国产黄色免费在线视频| 在线 av 中文字幕| 国产精品自产拍在线观看55亚洲 | 欧美成人午夜精品| 午夜福利免费观看在线| 91老司机精品| 欧美精品高潮呻吟av久久| 久久午夜亚洲精品久久| 在线观看66精品国产| 国产免费现黄频在线看| 久久人人爽av亚洲精品天堂| 国产免费av片在线观看野外av| 国产一区二区在线观看av| 在线观看免费视频日本深夜| 久久99一区二区三区| 啦啦啦视频在线资源免费观看| 国产主播在线观看一区二区| 精品一区二区三区四区五区乱码| av免费在线观看网站| 91麻豆精品激情在线观看国产 | 脱女人内裤的视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲男人天堂网一区| 久9热在线精品视频| 大片免费播放器 马上看| 国产一区有黄有色的免费视频| 女人久久www免费人成看片| 成人国产av品久久久| 宅男免费午夜| 黄片大片在线免费观看| 久久青草综合色| 999精品在线视频| 18禁国产床啪视频网站| 纵有疾风起免费观看全集完整版| 午夜激情久久久久久久| 在线 av 中文字幕| 高清视频免费观看一区二区| 美女主播在线视频| 丁香六月欧美| 女人高潮潮喷娇喘18禁视频| 亚洲性夜色夜夜综合| 亚洲成a人片在线一区二区| 久久精品国产综合久久久| 国产真人三级小视频在线观看| 建设人人有责人人尽责人人享有的| 久久亚洲精品不卡| 亚洲精华国产精华精| 一本色道久久久久久精品综合| 国产亚洲精品一区二区www | 12—13女人毛片做爰片一| 日本黄色日本黄色录像| 无人区码免费观看不卡 | 亚洲性夜色夜夜综合| 十分钟在线观看高清视频www| 精品高清国产在线一区| 每晚都被弄得嗷嗷叫到高潮| 免费在线观看完整版高清| 欧美乱妇无乱码| 天堂8中文在线网| 一区二区三区激情视频| 午夜福利一区二区在线看| 日本a在线网址| 亚洲avbb在线观看| 亚洲欧美激情在线| 国产精品免费一区二区三区在线 | 免费在线观看完整版高清| 丝袜人妻中文字幕| 午夜日韩欧美国产| avwww免费| 欧美成狂野欧美在线观看| 亚洲精品久久午夜乱码| 91字幕亚洲| 天天添夜夜摸| 黑人猛操日本美女一级片| 精品国产超薄肉色丝袜足j| 亚洲精品粉嫩美女一区| 日本一区二区免费在线视频| 亚洲欧美一区二区三区久久| 日韩欧美一区视频在线观看| 国产色视频综合| 精品第一国产精品| 丝瓜视频免费看黄片| 免费观看av网站的网址| 成人国语在线视频| 建设人人有责人人尽责人人享有的| 欧美黄色片欧美黄色片| 精品熟女少妇八av免费久了| 啦啦啦在线免费观看视频4| 99久久人妻综合| 免费女性裸体啪啪无遮挡网站| 亚洲国产欧美在线一区| av免费在线观看网站| 啦啦啦 在线观看视频| 久9热在线精品视频| 咕卡用的链子| 波多野结衣av一区二区av| 免费看十八禁软件| 久久青草综合色| 97在线人人人人妻| 久久人妻福利社区极品人妻图片| 可以免费在线观看a视频的电影网站| 97人妻天天添夜夜摸| 美女高潮喷水抽搐中文字幕| 精品一区二区三区四区五区乱码| 国产在视频线精品| 少妇被粗大的猛进出69影院| 国产精品成人在线| 国产亚洲精品一区二区www | 亚洲精品美女久久久久99蜜臀| 欧美另类亚洲清纯唯美| 国产精品98久久久久久宅男小说| 99re6热这里在线精品视频| 大片电影免费在线观看免费| 在线观看www视频免费| 欧美日韩亚洲综合一区二区三区_| 亚洲综合色网址| 成人特级黄色片久久久久久久 | 色婷婷av一区二区三区视频| 美国免费a级毛片| 精品视频人人做人人爽| av超薄肉色丝袜交足视频| 黑人欧美特级aaaaaa片| 丰满饥渴人妻一区二区三| 成年人免费黄色播放视频| 桃花免费在线播放| 又紧又爽又黄一区二区| 国产一区二区三区在线臀色熟女 | 天天添夜夜摸| 亚洲avbb在线观看| 精品一区二区三区视频在线观看免费 | kizo精华| 国产极品粉嫩免费观看在线| 精品卡一卡二卡四卡免费| 国产精品香港三级国产av潘金莲| 后天国语完整版免费观看| www日本在线高清视频| 一进一出好大好爽视频| e午夜精品久久久久久久| 每晚都被弄得嗷嗷叫到高潮| 狠狠婷婷综合久久久久久88av| 淫妇啪啪啪对白视频| 日韩成人在线观看一区二区三区| 一边摸一边抽搐一进一小说 | 久久 成人 亚洲| 国产一区二区三区视频了| 久久中文字幕一级| www.熟女人妻精品国产| 国产欧美日韩一区二区三区在线| 满18在线观看网站| 国产av国产精品国产| 日韩成人在线观看一区二区三区| 亚洲精品自拍成人| 国产免费视频播放在线视频| 久久中文看片网| 欧美一级毛片孕妇| 亚洲国产欧美日韩在线播放| 久久亚洲真实| 99精品在免费线老司机午夜| 亚洲成人免费av在线播放| 国产精品偷伦视频观看了| 一级毛片精品| 丁香欧美五月| 日本av免费视频播放| av不卡在线播放| 另类亚洲欧美激情| 久久午夜亚洲精品久久| 在线观看舔阴道视频| 久久亚洲精品不卡| 一区二区av电影网| 黑人巨大精品欧美一区二区mp4| 成年人免费黄色播放视频| 亚洲精品乱久久久久久| 不卡av一区二区三区| 99精品欧美一区二区三区四区| 日韩大片免费观看网站| 国产精品98久久久久久宅男小说| 亚洲人成77777在线视频| 国产人伦9x9x在线观看| avwww免费| 中文字幕色久视频| 欧美亚洲日本最大视频资源| 久久99一区二区三区| 精品高清国产在线一区| 成年动漫av网址| bbb黄色大片| 热re99久久精品国产66热6| 国产免费福利视频在线观看| 中文字幕av电影在线播放| 宅男免费午夜| 色综合欧美亚洲国产小说| 亚洲精品中文字幕在线视频| 满18在线观看网站| 王馨瑶露胸无遮挡在线观看| 国产免费av片在线观看野外av| 夜夜夜夜夜久久久久| 久久婷婷成人综合色麻豆| 狠狠狠狠99中文字幕| 2018国产大陆天天弄谢| 国产精品久久久久久精品古装| 精品免费久久久久久久清纯 | svipshipincom国产片| 一个人免费看片子| 男人操女人黄网站| 老司机深夜福利视频在线观看| 久久久久精品人妻al黑| 国产男女内射视频| 18禁观看日本| 99国产精品一区二区蜜桃av | 欧美日韩国产mv在线观看视频| 黄片小视频在线播放| 免费av中文字幕在线| 亚洲国产毛片av蜜桃av| 国产欧美日韩综合在线一区二区| 不卡一级毛片| 交换朋友夫妻互换小说| 中文字幕av电影在线播放| 亚洲第一av免费看| 视频区图区小说| 麻豆成人av在线观看| 性色av乱码一区二区三区2| 亚洲av成人一区二区三| 1024香蕉在线观看| 人妻久久中文字幕网| 国产精品影院久久| 成人手机av| 久久久精品区二区三区| 日韩大码丰满熟妇| 水蜜桃什么品种好| 叶爱在线成人免费视频播放| 最新在线观看一区二区三区| 无限看片的www在线观看| 人人妻人人爽人人添夜夜欢视频| 国产区一区二久久| 一级黄色大片毛片| 中文字幕另类日韩欧美亚洲嫩草| 91老司机精品| 午夜福利视频在线观看免费| aaaaa片日本免费| 老司机靠b影院| 亚洲国产精品一区二区三区在线| 美女福利国产在线| 99在线人妻在线中文字幕 | 国产成人啪精品午夜网站| 久久久国产欧美日韩av| av福利片在线| 亚洲午夜精品一区,二区,三区| 岛国毛片在线播放| 久久人人爽av亚洲精品天堂| 人妻一区二区av| 久久精品国产99精品国产亚洲性色 | 国产aⅴ精品一区二区三区波| h视频一区二区三区| 久久久精品国产亚洲av高清涩受| 不卡一级毛片| 麻豆国产av国片精品| 国产无遮挡羞羞视频在线观看| 久久久久久久久久久久大奶| 欧美性长视频在线观看| 亚洲av第一区精品v没综合| av不卡在线播放| 中文字幕另类日韩欧美亚洲嫩草| 黑人操中国人逼视频| 亚洲精品粉嫩美女一区| 国产在线精品亚洲第一网站| 午夜福利一区二区在线看| 91麻豆av在线| 国产日韩欧美在线精品| 99久久国产精品久久久| 999精品在线视频| 丝袜人妻中文字幕| 国产精品一区二区精品视频观看| 69av精品久久久久久 | 国产成人精品在线电影| 中文字幕av电影在线播放| 国产欧美亚洲国产| 一级毛片精品| 岛国毛片在线播放| 国产99久久九九免费精品| 国产伦理片在线播放av一区| 男女高潮啪啪啪动态图| 午夜福利一区二区在线看| 国产激情久久老熟女| 久久国产精品男人的天堂亚洲| 日本vs欧美在线观看视频| 欧美 日韩 精品 国产| 国产精品偷伦视频观看了| 亚洲熟女毛片儿| 满18在线观看网站| 色婷婷久久久亚洲欧美| 久久中文字幕人妻熟女| 亚洲欧美一区二区三区久久| 老司机影院毛片| 一本大道久久a久久精品| 精品免费久久久久久久清纯 | 国产精品一区二区精品视频观看| 啪啪无遮挡十八禁网站| 9热在线视频观看99| 久久久久久久精品吃奶| 无限看片的www在线观看| 99精国产麻豆久久婷婷| 亚洲国产欧美网| 国产三级黄色录像| 亚洲 欧美一区二区三区| 18禁国产床啪视频网站| 91大片在线观看| 国产欧美日韩一区二区三| 欧美久久黑人一区二区| 欧美黄色片欧美黄色片| 精品一区二区三区av网在线观看 | 超色免费av| 动漫黄色视频在线观看| 香蕉国产在线看| 国产无遮挡羞羞视频在线观看| 大香蕉久久成人网| 999久久久精品免费观看国产| 国产色视频综合| 在线观看66精品国产| 这个男人来自地球电影免费观看| 精品国产亚洲在线| 日韩有码中文字幕| 韩国精品一区二区三区| 欧美日韩福利视频一区二区| 纯流量卡能插随身wifi吗| 在线观看免费午夜福利视频| 国产一区二区三区视频了| 69精品国产乱码久久久| 飞空精品影院首页| 国产欧美日韩一区二区精品| 老熟妇仑乱视频hdxx| 免费在线观看黄色视频的| 黄色 视频免费看| 久9热在线精品视频| 色老头精品视频在线观看| 日本a在线网址| 老司机亚洲免费影院| 久久天躁狠狠躁夜夜2o2o| 黄色成人免费大全| 久久天躁狠狠躁夜夜2o2o| 黄色成人免费大全| 最近最新免费中文字幕在线| 欧美成人午夜精品| 嫁个100分男人电影在线观看| 亚洲精品成人av观看孕妇| 波多野结衣一区麻豆| 美女视频免费永久观看网站| 欧美精品一区二区免费开放| 亚洲国产欧美日韩在线播放| 女警被强在线播放| 99国产极品粉嫩在线观看| 一本色道久久久久久精品综合| 欧美黑人精品巨大| 人人妻人人添人人爽欧美一区卜| 欧美黑人欧美精品刺激| 精品少妇内射三级| 免费高清在线观看日韩| 日韩免费高清中文字幕av| 男女免费视频国产| 亚洲专区国产一区二区| 亚洲天堂av无毛| 成人国产一区最新在线观看| 一本—道久久a久久精品蜜桃钙片| 国产又色又爽无遮挡免费看| 制服人妻中文乱码| 美女高潮到喷水免费观看| 亚洲国产欧美网| 男女下面插进去视频免费观看| 亚洲专区国产一区二区| 麻豆国产av国片精品| 欧美+亚洲+日韩+国产| 亚洲第一欧美日韩一区二区三区 | 这个男人来自地球电影免费观看| 国产无遮挡羞羞视频在线观看| 热99re8久久精品国产| 亚洲第一av免费看| 国产男女内射视频| 一本一本久久a久久精品综合妖精| 国产免费视频播放在线视频| 久久人妻福利社区极品人妻图片| 老司机在亚洲福利影院| 成人手机av| 大型黄色视频在线免费观看| 色婷婷av一区二区三区视频| 国产xxxxx性猛交| 黄色怎么调成土黄色| 黄色视频,在线免费观看| 亚洲自偷自拍图片 自拍| 久久久久久久久久久久大奶| 亚洲三区欧美一区| 精品欧美一区二区三区在线| videosex国产| 精品国产超薄肉色丝袜足j| 国产激情久久老熟女| 亚洲九九香蕉| 中文字幕制服av| 成人18禁高潮啪啪吃奶动态图| 大型av网站在线播放| 十八禁高潮呻吟视频|