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

    基于蒙特卡洛仿真的高空尾渦運動特性

    2021-12-22 05:58:12魏志強李曉晨
    空軍工程大學(xué)學(xué)報 2021年5期
    關(guān)鍵詞:環(huán)量尾渦飛行高度

    魏志強,李曉晨

    (中國民航大學(xué)空中交通管理學(xué)院,天津,300300)

    飛機機翼在產(chǎn)生升力的過程中,上翼面及下翼面之間的壓強差將會導(dǎo)致在左右翼尖處形成轉(zhuǎn)向相反的翼尖渦流場[1]。當(dāng)后方飛機不慎進入前機的尾渦流場時,可能會發(fā)生失速、滾轉(zhuǎn)、急劇俯仰等危險情況,嚴重威脅飛行安全[2]。近年來,空域繁忙程度的增加以及空域資源的日趨緊張使得充分利用高空資源越來越迫切,其中尾渦遭遇問題是高空飛行高度層垂直間隔縮減的重要限制因素。

    空中交通流量的持續(xù)增長對空域利用效率提出更高要求。我國在8 900~12 500 m的平飛巡航高度區(qū)間內(nèi)實施縮小最低垂直間隔(reduced vertical separation minimum, RVSM),將飛行高度層之間的垂直間隔標準由600 m縮小到300 m,增加了空域容量。近年來,國際民航組織(international civil aviation organization,ICAO)開始著手研究在RVSM空域上擴高度層中實施縮小垂直間隔的可行性[3],但尚未考慮飛機尾流的消散和運動對下方飛機飛行安全的影響問題。

    目前針對尾渦流場參數(shù)的仿真計算在技術(shù)上可分為基于流體動力學(xué)方法的數(shù)值模擬技術(shù)和基于渦流演變機理的快速仿真計算技術(shù)。飛機尾渦流場的數(shù)值模擬主要包括大渦模擬、雷諾平均及分離渦模擬方法[4-6]。受計算機計算方法及運算能力的限制,在對飛機尾渦進行數(shù)值模擬時,有限的網(wǎng)格數(shù)量使得對尾渦形成階段的模擬較為清晰,但對于尾渦的消散和運動階段的數(shù)值模擬效果尚不夠理想。在尾渦流場參數(shù)的仿真計算方面,國外研究者通過建立數(shù)學(xué)模型描述尾渦消散機制。Holz?pfel等建立兩階段尾渦消散模型(two-phase wake vortex decay,D2P)來計算尾渦強度的衰減情況[7]。Sarpkaya等基于數(shù)值模擬數(shù)據(jù)的分析,認為尾渦的消散速率與雷諾數(shù)關(guān)系不大,而主要取決于大氣湍流度[8-9]。Proctor等建立了尾渦流場快速預(yù)測模型,并應(yīng)用于終端區(qū)尾流仿真系統(tǒng)平臺。這些數(shù)學(xué)模型由于響應(yīng)快速、運算高效、相對簡化而被廣泛應(yīng)用于尾渦流場參數(shù)的仿真計算中[10-11]。

    國內(nèi)研究者主要致力于尾渦流場建模及尾渦參數(shù)計算方面的研究。魏志強等利用“天河一號”超級計算機進行數(shù)值模擬實驗,計算分析了不同側(cè)風(fēng)下的渦量衰減、渦心速度等參數(shù)的變化情況[12]。沈淳等基于雷達探測反演方法預(yù)測飛機尾流行為及尾渦危險區(qū)域,為尾流動態(tài)安全間隔標準制定提供支撐[13]。谷潤平等建立多目標優(yōu)化模型,以誘導(dǎo)滾轉(zhuǎn)力矩系數(shù)反映尾渦遭遇嚴重度,應(yīng)用于飛機編隊飛行的前后機相對位置優(yōu)化[14]。魏志強等分析了飛機加裝翼尖小翼后尾流安全間隔的變化情況,以及航空器重新分類時的尾流遭遇風(fēng)險問題[15-16]。國內(nèi)外研究者在尾渦流場參數(shù)仿真計算和安全間隔標準的制定方面開展的研究主要針對中低空及飛機的起降階段,未對12 500 m以上高空空域的尾渦流場特性進行分析。

    為解決上述問題,文中首先建立了飛機尾渦物理模型。然后以某重型飛機為例分析12 500 m飛行高度以上的高空尾渦流場運動特性,根據(jù)蒙特卡洛方法仿真計算高空與中低空的尾渦下沉高度差異;最后分析高空尾渦渦核下沉運動的影響因素。研究結(jié)果為在RVSM空域上擴高度層中實施縮小垂直間隔標準提供參考,以期提高空域利用效率,降低高空尾渦遭遇風(fēng)險。

    1 飛機尾渦物理模型

    1.1 尾渦初始強度計算模型

    飛機尾渦一般使用渦旋環(huán)量Γ表示其強度。根據(jù)圓柱繞流氣動力知識,機翼微元上的升力計算公式為:

    l(x)=ρ∞v∞Γ(x)

    (1)

    式中:ρ∞為大氣密度;v∞為飛機真空速;Γ(x)為翼剖面附著渦環(huán)量;x為翼剖面與飛機縱軸之間的距離。對整個翼展積分,飛機升力L的計算公式為:

    (2)

    式中:B為飛機翼展。對于后掠翼飛機,翼剖面環(huán)量的計算公式為:

    (3)

    (4)

    所以,得尾渦初始強度計算公式如下:

    (5)

    式中:n為飛機法向過載;W為飛機質(zhì)量。由上式可以看出,尾渦初始環(huán)量的大小與飛機質(zhì)量、空氣密度、飛機速度、翼展等因素相關(guān)。

    1.2 尾渦消散過程

    綜合考慮尾渦消散的隨機特性,Holz?pfel和Robins應(yīng)用激光雷達開展多次尾渦探測試驗[17],在D2P模型基礎(chǔ)上添加隨機擾動項,形成尾渦隨機兩階段消散模型(probabilistic two-phase wake vortex decay,P2P)。尾渦的消散主要受大氣參數(shù)及飛機特性影響,P2P模型可以用來描述不同飛行高度處的尾渦消散情況。

    為了計算方便,通常使用無量綱基準參數(shù)。常用的無量綱基準參數(shù)包括參考下沉速度w0和尾渦參考時間t0,具體計算公式如下:

    (6)

    (7)

    式中:b0為尾渦初始渦核間距,通常為πB/4。

    根據(jù)P2P模型,尾渦的消散可分為2個階段:擴散階段及快速衰減階段。在擴散階段,尾渦消散的速率較為緩慢,尾渦環(huán)量計算公式如下:

    (8)

    (9)

    ε*=(εb0)1/3/w0

    (10)

    用浮力頻率N表征大氣層結(jié)穩(wěn)定度,無因次浮力頻率N*的計算公式如下:

    (11)

    (12)

    式中:hf為飛機的飛行高度;T為所在飛行高度的溫度;Cp為定壓比熱。與平流層相比,通常對流層內(nèi)大氣層結(jié)穩(wěn)定度較低。

    1.3 尾渦下沉模型

    尾渦形成后由于相互誘導(dǎo)作用而向下運動,渦核位置發(fā)生變化。渦核下沉速度vs及經(jīng)時間t渦核下沉高度hs的計算公式如下:

    (13)

    (14)

    在實際大氣中,尾渦的強度消散和渦核下沉運動會對下方飛機的飛行安全造成不利影響。通過計算并對比高空和中低空的尾渦運動差異,可以定量地評估在高空實施縮減尾流垂直間隔的可行性及安全性。

    2 尾渦運動特性計算分析

    2.1 計算模型的驗證分析

    SESAR和EUROCONTROL針對歐洲空域內(nèi)飛機巡航階段的尾渦遭遇安全問題聯(lián)合開展R-WAKE項目[18]。分別以形成尾渦流場的前機為A320、A330、B777,飛行高度為39 500 ft,在渦流耗散率ε及浮力頻率N均為0的條件下,參考文獻[18]給出R-WAKE項目的尾渦遭遇模擬數(shù)據(jù)。統(tǒng)計前機后方不同縱向位置處的尾渦渦核下沉高度,與第1節(jié)中模型的計算結(jié)果進行對比,如表1所示。相對誤差基本在4.2%范圍內(nèi),說明文中模型的計算精度可以接受。

    表1 尾渦渦核下沉高度數(shù)據(jù)對比

    2.2 仿真實驗條件

    考慮在飛機飛行過程中飛機的重量、速度及大氣參數(shù)的不確定性,采用蒙特卡洛方法(Monte-Carlo approach)對尾渦的運動情況進行仿真模擬。蒙特卡洛模擬屬于概率分析法,使用隨機變量代替常量,同時保證隨機變量符合一定的概率分布,使最終結(jié)果更加逼近使用常量時所得到的計算結(jié)果。蒙特卡洛仿真的具體實驗條件如表2所示。

    表2 蒙特卡洛仿真實驗條件

    2.3 尾渦渦核下沉計算分析

    根據(jù)式(1)~(5),使用蒙特卡洛方法對某重型機在不同飛行高度處的尾渦初始環(huán)量進行10 000次仿真計算,得到不同飛行高度處的尾渦初始環(huán)量均值如圖1所示。由圖1可以看出,高空中的尾渦初始環(huán)量較大,15 000 m飛行高度處的尾渦初始環(huán)量約為中低空尾渦初始環(huán)量最小值的3.1倍。因此在高空飛行范圍內(nèi),仍有尾渦遭遇不安全事件的發(fā)生和報道。

    圖1 蒙特卡洛仿真模擬的尾渦初始環(huán)量Γ0均值

    前機下方形成的尾渦危險區(qū)與后機所能承受的尾渦強度Γa有關(guān),尾渦危險區(qū)在垂直范圍上表示為尾渦環(huán)量消散到Γa時的渦核下沉高度。為反映不同飛機所能承受的尾渦強度不同,令Γa分別為150 m2/s、100 m2/s、50 m2/s,計算尾渦渦核下沉高度。依據(jù)機型性能數(shù)據(jù)和仿真實驗條件,使用蒙特卡洛方法對某重型機在不同飛行高度處的尾渦渦核下沉高度進行10 000次仿真計算,得到尾渦下沉高度平均值和標準差,仿真結(jié)果如圖2和圖3所示。

    從圖2可以看出,12 500 m以上的高空中,尾渦環(huán)量消散到Γa時渦核下沉高度隨飛行高度的增加而增加。這是因為飛行高度增加后,尾渦初始環(huán)量增大,因此由渦核的下沉高度計算公式可知,飛行高度增加,高空渦核的下沉速度相對較大,導(dǎo)致高空渦核下沉高度增加。從圖3可以看出,12 500 m以上的高空中,隨飛行高度增加,渦核下沉高度的標準差降低,因此外界因素的不確定性變化對高空渦核下沉運動的干擾較小。

    圖2 蒙特卡洛仿真模擬的渦核下沉高度均值

    圖3 蒙特卡洛仿真模擬的渦核下沉高度標準差

    與中低空相比,12 500 m以上高空中的渦核下沉高度在增大。為了計算12 500 m以上高空尾渦渦核下沉高度增大的幅度,設(shè)尾渦環(huán)量消散到Γa時,高空尾渦渦核最大下沉高度為hmax,中低空尾渦渦核最小下沉高度為hmin,則尾渦渦核下沉高度差值Δh的計算公式如下:

    Δh=hmax-hmin

    (15)

    2.4 高空與中低空的渦核下沉運動差異性分析

    計算不同蒙特卡洛仿真實驗次數(shù)下的Δh,統(tǒng)計平均值與相對標準差,結(jié)果如圖4及圖5所示??梢钥闯?,高空尾渦渦核下沉高度平均增加量為42.4~49.7 m,相對標準差為14.4%~17.6%。當(dāng)實驗次數(shù)超過10 000次后,計算結(jié)果與仿真實驗次數(shù)基本無關(guān)聯(lián)。

    圖4 不同實驗次數(shù)下的Δh平均值

    圖5 不同實驗次數(shù)下的Δh相對標準差

    當(dāng)環(huán)量消散到100 m2/s,統(tǒng)計50 000次蒙特卡洛仿真實驗下的Δh分布情況,繪制頻數(shù)分布直方圖及對應(yīng)的概率密度函數(shù)曲線,如圖6所示。由圖6可以看出,Δh的仿真值大部分分布在35~55 m范圍內(nèi),呈正態(tài)分布規(guī)律。

    圖6 Δh的分布概率密度函數(shù)

    表3為Δh的具體區(qū)間分布情況。由表3可以看出,Δh最大值所在區(qū)間為80~85 m,最小值所在區(qū)間為15~20 m。渦核下沉高度差值在45~50 m區(qū)間內(nèi)出現(xiàn)次數(shù)最多,占比約為28.3%。根據(jù)蒙特卡洛仿真模擬得到的計算數(shù)據(jù),當(dāng)飛機參數(shù)及外界因素發(fā)生不確定性變化時,高空與中低空尾渦渦核下沉高度存在較大差異。

    表3 Δh的區(qū)間分布

    3 影響因素分析

    飛機特性及大氣條件均會對尾渦的消散及運動造成影響。此節(jié)分析飛機質(zhì)量以及大氣湍流度的變化對RVSM空域上擴高度層中的尾渦渦核下沉高度的影響。

    3.1 飛機質(zhì)量對尾渦下沉高度的影響

    飛機的形狀和質(zhì)量影響形成尾渦的初始強度,進而影響尾渦的消散和渦核的運動,其中飛機質(zhì)量是主要的影響因素[19]。為了分析飛機質(zhì)量對尾渦下沉高度的影響,按照飛機質(zhì)量350~390 t,其余初始條件同2.2節(jié)的條件,由尾渦物理模型計算15 000 m飛行高度處當(dāng)尾渦環(huán)量消散到100 m2/s時對應(yīng)的尾渦下沉高度。圖7為不同飛機質(zhì)量下進行10 000次蒙特卡洛仿真實驗數(shù)據(jù)的統(tǒng)計值,包括中位數(shù)、60%~90%概率值等。

    圖7 不同飛機質(zhì)量下尾渦下沉高度變化趨勢

    從圖7可以看出,高空尾渦渦核下沉高度隨飛機質(zhì)量的增加而增加。這是因為飛機質(zhì)量增加后,導(dǎo)致尾渦初始環(huán)量增加,渦核的下沉速度相對較大,因此由渦核下沉高度的計算公式可知,飛機質(zhì)量增加,導(dǎo)致渦核下沉高度增加。飛機質(zhì)量從350 t增加到390 t的過程中,渦核下沉高度最大增加4.4 m,變化范圍為1.7%~2.1%。

    3.2 大氣湍流度對尾渦下沉高度的影響

    高空中的大氣湍流度較低,為了分析大氣湍流水平對尾渦渦核下沉高度的影響,按照渦流耗散率10-5~10-3m2/s3,其余初始條件同2.2節(jié)的條件計算15 000 m飛行高度處當(dāng)尾渦環(huán)量消散到100 m2/s時對應(yīng)的尾渦下沉高度。圖8為不同大氣湍流度下進行10 000次蒙特卡洛仿真實驗數(shù)據(jù)的統(tǒng)計值,包括中位數(shù)、60%~90%概率值等。

    圖8 不同大氣湍流度下尾渦下沉高度變化趨勢

    從圖8可以看出,隨著渦流耗散率增大,高空尾渦渦核下沉高度減小。這是因為渦流耗散率增大,意味著大氣紊亂程度增加,使得尾渦消散進入快速衰減階段的時間提前[20],從而加快尾渦消散過程,渦核下沉高度降低。

    4 結(jié)語

    1)隨著飛行高度的增加,12 500 m以上高空空域內(nèi)的渦核下沉高度在增加。與中低空相比,使用蒙特卡洛模擬得到的渦核下沉高度增加量均值為42.4~49.7 m。

    2)飛機質(zhì)量通過改變尾渦初始強度進而影響尾渦渦核運動;大氣湍流度通過改變尾渦消散過程進而影響尾渦渦核運動。飛機質(zhì)量的增加會導(dǎo)致高空尾渦渦核下沉高度增加,大氣湍流度的增加會導(dǎo)致尾渦渦核下沉高度減小。因此飛行參數(shù)及大氣湍流條件的改變可以降低尾渦影響范圍,提高尾渦遭遇安全性。

    3)空域資源的日趨緊張使得進一步提升飛機升限、充分利用高空資源越來越迫切,在12 500 m以上空域內(nèi)實施縮小垂直間隔時,應(yīng)警惕前機尾渦所帶來的不利影響。下一步擬通過流場數(shù)值模擬方法或探測系統(tǒng)獲得大量的高空尾渦衰變及運動數(shù)據(jù),進一步研究高空飛行中的尾渦流場演化特性。

    猜你喜歡
    環(huán)量尾渦飛行高度
    不同B-V頻率下的飛機尾渦數(shù)值模擬研究
    高空巡航階段的飛機尾渦流場演化特性研究
    飛行參數(shù)對六旋翼植保無人機霧滴在荔枝樹冠層沉積分布的影響
    葉輪出口環(huán)量非線性分布條件下混流泵性能研究
    等-變環(huán)量設(shè)計葉片軸流風(fēng)機性能研究
    流體機械(2020年7期)2020-09-10 10:00:18
    基于模式函數(shù)和變分法的螺旋槳最佳環(huán)量計算方法
    簡析二次雷達高度信息與飛機實際高度的關(guān)系
    基于激光雷達回波的動態(tài)尾渦特征參數(shù)計算
    干擾板作用下飛機尾渦流場近地演變機理研究
    FAA頒新政使小型無人機飛行高度翻倍
    航空模型(2016年5期)2016-07-25 08:59:26
    99在线人妻在线中文字幕| 久久6这里有精品| 日韩欧美一区二区三区在线观看| 精品久久国产蜜桃| 九九热线精品视视频播放| a级毛片a级免费在线| 中文字幕熟女人妻在线| 在线观看美女被高潮喷水网站| 最近手机中文字幕大全| 成人二区视频| 成熟少妇高潮喷水视频| 直男gayav资源| 日韩欧美国产在线观看| 18+在线观看网站| 欧美激情国产日韩精品一区| 成人特级黄色片久久久久久久| 亚洲精品亚洲一区二区| 免费av毛片视频| 亚洲最大成人av| 毛片一级片免费看久久久久| 国产真实乱freesex| 国产精品伦人一区二区| 久久人妻av系列| 91精品一卡2卡3卡4卡| 少妇熟女aⅴ在线视频| 亚洲国产高清在线一区二区三| 久久99精品国语久久久| 美女黄网站色视频| 日韩强制内射视频| 国产精品久久视频播放| 小说图片视频综合网站| 国产精品伦人一区二区| 日韩国内少妇激情av| 国内精品美女久久久久久| 亚洲美女视频黄频| 亚洲av一区综合| 美女内射精品一级片tv| 免费看av在线观看网站| 成人午夜精彩视频在线观看| 插逼视频在线观看| 亚洲真实伦在线观看| 国产精品1区2区在线观看.| 天天躁夜夜躁狠狠久久av| 亚洲欧洲国产日韩| 超碰av人人做人人爽久久| 91久久精品国产一区二区成人| 99在线视频只有这里精品首页| 一级av片app| 真实男女啪啪啪动态图| 不卡视频在线观看欧美| 欧美三级亚洲精品| 深夜精品福利| 欧美bdsm另类| 在现免费观看毛片| 中文字幕精品亚洲无线码一区| 97超碰精品成人国产| 少妇丰满av| 亚洲国产精品国产精品| 男女啪啪激烈高潮av片| 久久九九热精品免费| 日本免费一区二区三区高清不卡| 欧美xxxx性猛交bbbb| 12—13女人毛片做爰片一| 啦啦啦啦在线视频资源| 亚洲真实伦在线观看| 美女高潮的动态| 日韩欧美精品v在线| 人妻少妇偷人精品九色| av在线蜜桃| 美女cb高潮喷水在线观看| 色哟哟哟哟哟哟| 1024手机看黄色片| 成人美女网站在线观看视频| 亚洲av免费在线观看| 日韩大尺度精品在线看网址| 少妇被粗大猛烈的视频| 精品日产1卡2卡| 亚洲欧美成人综合另类久久久 | 一区二区三区高清视频在线| 日韩av在线大香蕉| 国产毛片a区久久久久| 99riav亚洲国产免费| 在线观看66精品国产| 精品午夜福利在线看| 午夜精品国产一区二区电影 | 一进一出抽搐gif免费好疼| 日韩高清综合在线| 国产成人精品久久久久久| 亚洲精华国产精华液的使用体验 | 亚洲内射少妇av| 国产又黄又爽又无遮挡在线| 日韩欧美一区二区三区在线观看| 联通29元200g的流量卡| 美女脱内裤让男人舔精品视频 | 中出人妻视频一区二区| 国产亚洲5aaaaa淫片| 亚洲图色成人| 美女xxoo啪啪120秒动态图| 久久久欧美国产精品| 午夜福利在线观看吧| 国产爱豆传媒在线观看| 成人午夜高清在线视频| 不卡一级毛片| 国产美女午夜福利| 观看美女的网站| 亚洲人与动物交配视频| 国产高清激情床上av| 直男gayav资源| 观看美女的网站| 色综合亚洲欧美另类图片| 亚洲在久久综合| 少妇人妻精品综合一区二区 | 最近视频中文字幕2019在线8| 欧美区成人在线视频| 亚洲精品456在线播放app| 国产精品日韩av在线免费观看| 日韩欧美 国产精品| av专区在线播放| 麻豆国产av国片精品| 亚州av有码| www.av在线官网国产| 国产精品一区www在线观看| 国产成人精品婷婷| 免费看av在线观看网站| 国产精品一区二区三区四区久久| 久久久久网色| 久久久久久久久久久免费av| 国产精品久久久久久精品电影| 国产精品人妻久久久久久| 亚洲七黄色美女视频| 五月玫瑰六月丁香| 久久精品影院6| 国产成人精品婷婷| 精品久久国产蜜桃| 91在线精品国自产拍蜜月| 亚洲aⅴ乱码一区二区在线播放| 久久久久九九精品影院| 日本一二三区视频观看| 乱系列少妇在线播放| 午夜免费激情av| 美女cb高潮喷水在线观看| 最后的刺客免费高清国语| 亚洲国产精品成人综合色| 男女视频在线观看网站免费| 一区二区三区四区激情视频 | 不卡视频在线观看欧美| 婷婷亚洲欧美| h日本视频在线播放| 黄色一级大片看看| 熟女电影av网| 三级国产精品欧美在线观看| 日日啪夜夜撸| 免费观看a级毛片全部| 色噜噜av男人的天堂激情| 国产精品国产三级国产av玫瑰| 精品久久久久久成人av| 亚洲人成网站在线播| 亚洲av成人精品一区久久| 特级一级黄色大片| 黄色一级大片看看| 亚洲欧美精品自产自拍| 国内精品美女久久久久久| 亚洲欧美日韩无卡精品| 国产高清激情床上av| 欧美一区二区国产精品久久精品| а√天堂www在线а√下载| 在线观看午夜福利视频| 亚洲va在线va天堂va国产| 老司机影院成人| 久久99热6这里只有精品| 国产av麻豆久久久久久久| 欧美激情国产日韩精品一区| 久久鲁丝午夜福利片| .国产精品久久| 国产免费一级a男人的天堂| 村上凉子中文字幕在线| 人人妻人人澡欧美一区二区| 人人妻人人澡欧美一区二区| 99热只有精品国产| 一夜夜www| 亚洲成人久久爱视频| 九九爱精品视频在线观看| 国产亚洲av片在线观看秒播厂 | 国产老妇女一区| 国产精品无大码| 日本色播在线视频| 午夜a级毛片| 久久精品国产自在天天线| 丰满的人妻完整版| 国产蜜桃级精品一区二区三区| 精品午夜福利在线看| 亚洲av熟女| 亚洲精品久久国产高清桃花| 狂野欧美白嫩少妇大欣赏| 直男gayav资源| 久久久久久久久中文| 熟女电影av网| 国产 一区精品| 亚洲内射少妇av| 一级毛片电影观看 | .国产精品久久| 蜜臀久久99精品久久宅男| 亚洲内射少妇av| 精品久久久久久久久av| 91久久精品国产一区二区三区| 好男人视频免费观看在线| 美女高潮的动态| 国产成人午夜福利电影在线观看| 九九在线视频观看精品| av在线蜜桃| 亚洲欧美成人综合另类久久久 | 亚洲一级一片aⅴ在线观看| 九九热线精品视视频播放| 亚洲婷婷狠狠爱综合网| 午夜精品国产一区二区电影 | 亚洲五月天丁香| 亚洲va在线va天堂va国产| 久久久久性生活片| 成年av动漫网址| 国产精品99久久久久久久久| 永久网站在线| 99热只有精品国产| 国产精品99久久久久久久久| 欧美区成人在线视频| 亚洲欧美成人综合另类久久久 | 色哟哟哟哟哟哟| 中文欧美无线码| 亚洲国产欧美人成| 国产老妇伦熟女老妇高清| 蜜桃亚洲精品一区二区三区| 婷婷色av中文字幕| 国产综合懂色| h日本视频在线播放| av天堂在线播放| 国产精品免费一区二区三区在线| 99久久无色码亚洲精品果冻| 久99久视频精品免费| 欧美极品一区二区三区四区| 最近的中文字幕免费完整| 日韩一区二区视频免费看| 国产精品综合久久久久久久免费| 国产精品,欧美在线| 午夜福利视频1000在线观看| 欧美极品一区二区三区四区| 又黄又爽又刺激的免费视频.| 亚洲最大成人中文| 一区二区三区高清视频在线| 国产精品国产三级国产av玫瑰| 性色avwww在线观看| 精品久久久久久成人av| 天堂av国产一区二区熟女人妻| 三级国产精品欧美在线观看| 日韩三级伦理在线观看| 卡戴珊不雅视频在线播放| 国产精品久久久久久久久免| 欧美潮喷喷水| 成人毛片60女人毛片免费| 亚洲18禁久久av| 亚洲av免费高清在线观看| 国产91av在线免费观看| 久久精品国产亚洲av涩爱 | 国产成人影院久久av| 国产精品蜜桃在线观看 | 久久精品人妻少妇| 美女国产视频在线观看| 91在线精品国自产拍蜜月| 韩国av在线不卡| 联通29元200g的流量卡| 免费看光身美女| 国产真实伦视频高清在线观看| 精品久久久久久久久av| 最近手机中文字幕大全| 午夜爱爱视频在线播放| 可以在线观看毛片的网站| 日韩欧美 国产精品| 色吧在线观看| 欧美在线一区亚洲| 亚洲最大成人手机在线| 级片在线观看| 成人亚洲欧美一区二区av| 国产午夜精品久久久久久一区二区三区| 长腿黑丝高跟| 一本久久中文字幕| 老女人水多毛片| 日韩在线高清观看一区二区三区| 黄色日韩在线| 18禁黄网站禁片免费观看直播| 精品免费久久久久久久清纯| 一级av片app| www.色视频.com| 丰满的人妻完整版| 一区二区三区免费毛片| 九九爱精品视频在线观看| 国产精品1区2区在线观看.| 99久久中文字幕三级久久日本| 国产亚洲av嫩草精品影院| 小蜜桃在线观看免费完整版高清| 国产亚洲精品久久久久久毛片| 亚洲真实伦在线观看| 校园人妻丝袜中文字幕| 嘟嘟电影网在线观看| 人妻制服诱惑在线中文字幕| 3wmmmm亚洲av在线观看| 免费观看人在逋| 女的被弄到高潮叫床怎么办| 久久久久国产网址| 久久精品夜夜夜夜夜久久蜜豆| 夜夜爽天天搞| 蜜桃亚洲精品一区二区三区| 人妻制服诱惑在线中文字幕| 99久国产av精品国产电影| 成熟少妇高潮喷水视频| 少妇高潮的动态图| 中文欧美无线码| 天美传媒精品一区二区| 亚洲一区高清亚洲精品| 天天一区二区日本电影三级| 婷婷色综合大香蕉| 两个人视频免费观看高清| 桃色一区二区三区在线观看| 国产大屁股一区二区在线视频| 久久久久久久久久成人| 亚洲中文字幕日韩| 九色成人免费人妻av| 国产麻豆成人av免费视频| 欧美日本视频| av天堂在线播放| 国产一区二区三区在线臀色熟女| 久久人人爽人人片av| 国产探花极品一区二区| 村上凉子中文字幕在线| 一区二区三区四区激情视频 | 日韩国内少妇激情av| 搡女人真爽免费视频火全软件| 亚洲精品国产av成人精品| 国产亚洲精品av在线| 日韩欧美一区二区三区在线观看| 久久九九热精品免费| 国产探花在线观看一区二区| 日韩精品青青久久久久久| 久久久精品欧美日韩精品| 蜜桃亚洲精品一区二区三区| 国产精品一区二区在线观看99 | 久久久久九九精品影院| 又爽又黄a免费视频| 欧美日韩一区二区视频在线观看视频在线 | 欧美日韩综合久久久久久| 久久99热这里只有精品18| 国产淫片久久久久久久久| 亚洲欧美成人精品一区二区| 少妇被粗大猛烈的视频| 亚洲精品影视一区二区三区av| 亚洲美女搞黄在线观看| 寂寞人妻少妇视频99o| 亚洲国产欧美人成| 久久久国产成人精品二区| 网址你懂的国产日韩在线| 国产精品,欧美在线| 国产精品久久久久久精品电影小说 | 国产精品综合久久久久久久免费| 国产成人福利小说| 国产亚洲欧美98| 人体艺术视频欧美日本| 我的女老师完整版在线观看| 久久久色成人| 亚洲欧美日韩高清在线视频| 最好的美女福利视频网| 久久6这里有精品| 久久久a久久爽久久v久久| 日韩欧美 国产精品| 我要搜黄色片| 中文字幕人妻熟人妻熟丝袜美| 两个人的视频大全免费| 亚洲精品国产成人久久av| 亚洲av电影不卡..在线观看| 全区人妻精品视频| 又粗又硬又长又爽又黄的视频 | 麻豆一二三区av精品| 天天躁夜夜躁狠狠久久av| 中国国产av一级| 又爽又黄a免费视频| 国产精品一区二区三区四区免费观看| 少妇人妻一区二区三区视频| 69av精品久久久久久| 亚洲成a人片在线一区二区| 亚洲性久久影院| 国产私拍福利视频在线观看| 中文字幕精品亚洲无线码一区| av在线亚洲专区| 免费不卡的大黄色大毛片视频在线观看 | 亚洲国产日韩欧美精品在线观看| 哪个播放器可以免费观看大片| 在线观看午夜福利视频| 成人欧美大片| 亚洲成人精品中文字幕电影| 亚洲五月天丁香| 亚洲av中文字字幕乱码综合| 亚洲自拍偷在线| 日韩亚洲欧美综合| 亚洲,欧美,日韩| 久久精品国产清高在天天线| 国产视频内射| 欧美区成人在线视频| 久久久精品欧美日韩精品| 如何舔出高潮| 中文亚洲av片在线观看爽| 激情 狠狠 欧美| 99久久人妻综合| 亚洲国产高清在线一区二区三| 亚洲aⅴ乱码一区二区在线播放| 午夜福利在线观看吧| 中文字幕免费在线视频6| 国产极品精品免费视频能看的| 成人漫画全彩无遮挡| 又爽又黄a免费视频| 国产爱豆传媒在线观看| 日日摸夜夜添夜夜爱| 青春草视频在线免费观看| 精品人妻视频免费看| 国产精华一区二区三区| 麻豆一二三区av精品| 少妇人妻一区二区三区视频| 久久久国产成人免费| 久久久久久久久久成人| 黑人高潮一二区| 欧美3d第一页| 精品不卡国产一区二区三区| 在线免费观看的www视频| 直男gayav资源| 亚洲欧美日韩无卡精品| videossex国产| 亚洲国产日韩欧美精品在线观看| 亚洲欧美日韩无卡精品| 男插女下体视频免费在线播放| 欧美高清成人免费视频www| avwww免费| 桃色一区二区三区在线观看| 日本-黄色视频高清免费观看| 嫩草影院精品99| 中文在线观看免费www的网站| 亚洲激情五月婷婷啪啪| 深爱激情五月婷婷| 日韩大尺度精品在线看网址| 日韩精品有码人妻一区| 日韩一区二区三区影片| 亚洲经典国产精华液单| 国产精品人妻久久久影院| 午夜福利视频1000在线观看| 老师上课跳d突然被开到最大视频| 午夜福利在线在线| 久久午夜福利片| 日韩欧美在线乱码| 亚洲激情五月婷婷啪啪| 色播亚洲综合网| 国内久久婷婷六月综合欲色啪| 国产精品1区2区在线观看.| 午夜福利成人在线免费观看| 国产毛片a区久久久久| 欧美性猛交黑人性爽| 男女那种视频在线观看| av免费观看日本| 热99re8久久精品国产| 欧美不卡视频在线免费观看| 欧美性猛交╳xxx乱大交人| 成人美女网站在线观看视频| 久久精品夜夜夜夜夜久久蜜豆| 久久99热6这里只有精品| 悠悠久久av| 18禁裸乳无遮挡免费网站照片| 亚洲欧美日韩东京热| 有码 亚洲区| 久久久精品大字幕| 97人妻精品一区二区三区麻豆| 久久精品夜色国产| 全区人妻精品视频| 免费人成视频x8x8入口观看| 久久亚洲国产成人精品v| 99在线人妻在线中文字幕| 又爽又黄a免费视频| avwww免费| 午夜福利成人在线免费观看| 只有这里有精品99| 亚洲人与动物交配视频| 深夜精品福利| 亚洲不卡免费看| 少妇被粗大猛烈的视频| 麻豆精品久久久久久蜜桃| 12—13女人毛片做爰片一| 国产亚洲av嫩草精品影院| 在线观看免费视频日本深夜| 综合色av麻豆| 亚洲在线自拍视频| 色吧在线观看| 内地一区二区视频在线| 精品久久久久久久久亚洲| 午夜激情福利司机影院| 悠悠久久av| 禁无遮挡网站| 熟女电影av网| 亚洲婷婷狠狠爱综合网| 精品少妇黑人巨大在线播放 | 91久久精品国产一区二区三区| 欧美一级a爱片免费观看看| 亚州av有码| 丝袜喷水一区| 国内揄拍国产精品人妻在线| 国产精品久久久久久久久免| 久久久午夜欧美精品| 午夜激情欧美在线| 欧美精品一区二区大全| 在现免费观看毛片| 寂寞人妻少妇视频99o| 最后的刺客免费高清国语| 亚洲经典国产精华液单| av视频在线观看入口| 久久久久免费精品人妻一区二区| 哪里可以看免费的av片| 又黄又爽又刺激的免费视频.| 黄片wwwwww| 美女高潮的动态| 亚洲精品乱码久久久久久按摩| 99久国产av精品国产电影| 国产熟女欧美一区二区| 欧美xxxx黑人xx丫x性爽| 色综合色国产| 国产黄色小视频在线观看| 欧美人与善性xxx| 一级毛片久久久久久久久女| 午夜免费男女啪啪视频观看| 变态另类丝袜制服| 联通29元200g的流量卡| 亚洲精品色激情综合| 国产综合懂色| 国产精品av视频在线免费观看| 欧美极品一区二区三区四区| 日产精品乱码卡一卡2卡三| 少妇熟女欧美另类| 欧美+亚洲+日韩+国产| 久99久视频精品免费| 91久久精品电影网| 亚洲国产色片| 欧美一区二区精品小视频在线| 免费观看的影片在线观看| 老熟妇乱子伦视频在线观看| 国产高清有码在线观看视频| 亚洲国产精品合色在线| 一级毛片久久久久久久久女| 91狼人影院| 国产黄片视频在线免费观看| 尾随美女入室| 中文字幕熟女人妻在线| 国产在视频线在精品| 国产蜜桃级精品一区二区三区| 亚洲久久久久久中文字幕| 男女做爰动态图高潮gif福利片| 色尼玛亚洲综合影院| 99久久精品热视频| 国产91av在线免费观看| 婷婷色综合大香蕉| 欧美另类亚洲清纯唯美| 成人三级黄色视频| 晚上一个人看的免费电影| 九九在线视频观看精品| 18禁在线无遮挡免费观看视频| 亚洲av二区三区四区| 性色avwww在线观看| 日本黄大片高清| 国产黄a三级三级三级人| 日韩一区二区视频免费看| 成年女人永久免费观看视频| 亚洲人成网站在线观看播放| 嫩草影院精品99| 日韩成人伦理影院| 亚洲av第一区精品v没综合| 伦理电影大哥的女人| 夜夜爽天天搞| 女同久久另类99精品国产91| 97热精品久久久久久| 欧美三级亚洲精品| 久久久久网色| 边亲边吃奶的免费视频| 日韩 亚洲 欧美在线| 日日啪夜夜撸| 国产精品一二三区在线看| www日本黄色视频网| 午夜福利成人在线免费观看| 久久亚洲精品不卡| 中文字幕人妻熟人妻熟丝袜美| 爱豆传媒免费全集在线观看| 国内久久婷婷六月综合欲色啪| 国产精品国产三级国产av玫瑰| 成人综合一区亚洲| 久久久久久久午夜电影| 欧美日韩国产亚洲二区| 两个人视频免费观看高清| 国语自产精品视频在线第100页| 高清毛片免费观看视频网站| 高清午夜精品一区二区三区 | 看片在线看免费视频| www日本黄色视频网| 亚洲四区av| 夫妻性生交免费视频一级片| 欧美日韩乱码在线| 日韩视频在线欧美| 婷婷六月久久综合丁香| 久久精品国产99精品国产亚洲性色| 身体一侧抽搐| 一级毛片我不卡| 伊人久久精品亚洲午夜| 色哟哟·www| 成人一区二区视频在线观看| 久久精品综合一区二区三区| 亚洲国产日韩欧美精品在线观看| 最近中文字幕高清免费大全6| 一个人看的www免费观看视频|