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

    納秒激光等離子體減阻數(shù)值模擬

    2014-06-07 07:15:52馬麗華盧朝梁陳永安
    激光與紅外 2014年2期
    關(guān)鍵詞:納秒來流馬赫數(shù)

    陳 豪,石 磊,馬麗華,徐 景,盧朝梁,陳永安

    (1.空軍工程大學(xué)信息與導(dǎo)航學(xué)院,陜西西安710077;2.95561部隊(duì),西藏拉薩850002)

    納秒激光等離子體減阻數(shù)值模擬

    陳 豪1,石 磊1,馬麗華1,徐 景2,盧朝梁1,陳永安1

    (1.空軍工程大學(xué)信息與導(dǎo)航學(xué)院,陜西西安710077;2.95561部隊(duì),西藏拉薩850002)

    針對激光等離子體減阻技術(shù)機(jī)理,采用納秒激光,數(shù)值模擬激光等離子體在流場中的演化過程,分析關(guān)鍵參數(shù)對納秒激光等離子體減阻性能的影響。結(jié)果表明:納秒激光能最大程度地提高激光等離子體減阻性能,阻力減小的百分比達(dá)99%,低阻力維持時(shí)間是入射激光持續(xù)時(shí)間的103倍;隨著來流馬赫數(shù)的增大,空氣來流強(qiáng)度增強(qiáng),導(dǎo)致減小阻力的百分比減小;隨著激光能量的增加,激光引致的沖擊波強(qiáng)度增大,使得減小阻力的百分比增加;隨著激光聚焦擊穿位置的增大,減小阻力的百分比減小,但低阻力維持的時(shí)間明顯增加。

    納秒激光;等離子體減阻;減小阻力的百分比;來流馬赫數(shù);激光能量;激光聚焦位置

    1 引 言

    1978年,美國倫瑟勒工學(xué)院的Myrabo率先提出利用激光在飛行器前端聚焦擊穿空氣形成高溫高壓等離子體,減小超聲速飛行器氣動阻力的新概念[1],其基本原理是在飛行器頭部通過激光等離子體形成高壓擾動,改變超聲速來流形成的弓形脫體激波波形,變?yōu)樾奔げǎ瑥亩鴾p小飛行器受到的阻力[2]。由于具有可控性好,減阻效率高等特點(diǎn),受到國內(nèi)外學(xué)者的廣泛關(guān)注[3]。近年來,基于微秒級CO2氣體激光,相繼開展了減阻機(jī)理、實(shí)驗(yàn)和數(shù)值模擬研究[4-7],討論分析了激光脈沖能量[8]、激光聚焦擊穿位置[9]及來流馬赫數(shù)[10]對減阻性能的影響,目前利用激光等離子體技術(shù)可以減小70%以上的波阻。

    隨著激光器水平的發(fā)展,體積更小、脈沖更窄、峰值功率更高的納秒固體激光器在航空航天應(yīng)用中的優(yōu)勢得以顯現(xiàn)。本文采用納秒激光,開展激光等離子體減阻技術(shù)數(shù)值模擬,分析納秒激光等離子體減阻性能。

    2 計(jì)算模型

    2.1 計(jì)算構(gòu)型及計(jì)算區(qū)域劃分

    數(shù)值計(jì)算中采用2維軸對稱計(jì)算模型,使用拋物形飛行器模型,滿足方程y2+z2=2px,p=10為拋物參數(shù),底部直徑為D=30 mm,高為H=11.25 mm,如圖1所示。為降低計(jì)算區(qū)域復(fù)雜度,基于計(jì)算分區(qū)思想,將計(jì)算區(qū)域劃分為3個區(qū)域,通過Gambit軟件生成計(jì)算網(wǎng)格。

    圖1 計(jì)算模型Fig.1 Computation model

    2.2 初邊值條件

    數(shù)值模擬中邊界條件如圖1所示,其中inflow為強(qiáng)制高超聲速入口邊界條件,wall為滑移固壁邊界條件,outflow為出口邊界條件,axis為對稱軸。計(jì)算中假設(shè)環(huán)境氣體是距離地球表面30 km高空的大氣,具體參數(shù)如表1所示。

    表1 環(huán)境氣體參數(shù)Tab.1 Parameters of ambient air

    2.3 數(shù)值模擬方案

    不考慮粘性和質(zhì)量力,控制方程采用Navier-Stokes方程。數(shù)值模擬計(jì)算基于有限體積法,采用Roe-FDS通量格式,通過基于Van Albada通量限制器的單調(diào)逆風(fēng)格式(MUSCL)方法獲得三階空間精度,時(shí)間積分采用二階精度的耦合隱式算法;用離散轉(zhuǎn)移方法(DTRM)分析等離子體輻射;使用標(biāo)準(zhǔn)κε雙方程模型進(jìn)行湍流計(jì)算。

    3 計(jì)算結(jié)果與分析

    3.1 納秒激光等離子體減阻流場演變

    首先以自由來流馬赫數(shù)Ma=5.5為例,分析納秒激光等離子體減阻流場演變。計(jì)算中采用單脈沖激光能量E=2J,脈沖寬度tp=20ns的納秒激光,并假設(shè)激光聚焦位置L=0.6D(L為激光聚焦點(diǎn)到飛行器頭部的距離,D為飛行器的直徑)。圖2顯示了計(jì)算時(shí)間t=0~50μs內(nèi)飛行器的阻力曲線。

    圖2 阻力曲線Fig.2 Air drag curve in the case of Ma=5.5

    從圖中可知,來流馬赫數(shù)Ma=5.5條件下,飛行器受到的穩(wěn)態(tài)阻力為57.86N,激光作用后的最小阻力僅有0.53N,減小阻力的百分比ηD達(dá)99%,明顯高于微秒激光70%的減阻效果[3],另一方面低阻力(小于穩(wěn)態(tài)阻力)維持時(shí)間達(dá)23μs,是入射激光持續(xù)時(shí)間20 ns的1150倍。

    分析曲線變化趨勢可知,納秒激光經(jīng)聚焦擊穿后產(chǎn)生的等離子體沖擊波在4.3μs時(shí)傳播到飛行器外壁面,并開始產(chǎn)生作用力,導(dǎo)致阻力迅速升高,在5.8μs時(shí)達(dá)到峰值,最大值為215 N。此后阻力迅速減小,在t=14μs時(shí)下降趨勢變緩,在21μs時(shí)達(dá)到最小值0.53 N。在t=22μs時(shí),阻力再次迅速升高,然而在t=23μs后,納秒激光等離子體的影響幾乎消失,在自由空氣來流作用下,阻力緩慢升高,而后下降,收斂至穩(wěn)態(tài)阻力。

    圖3給出了t=0~50μs之間流場內(nèi)的壓強(qiáng)等值線。

    在t=4.3μs時(shí),激光等離子體沖擊波剛剛傳播到飛行器外壁面,波形呈非對稱性傳播,前端波陣面(逆空氣來流方向傳播)受空氣來流影響傳播較慢,后端波陣面?zhèn)鞑ポ^快。

    在t=5.8μs時(shí),阻力迅速升高達(dá)到最大值。原因是納秒激光引致的沖擊波強(qiáng)波陣面與飛行器相互作用面積增加,使得飛行器外壁面壓力增大,導(dǎo)致出現(xiàn)阻力峰值。

    在t=10μs時(shí),飛行器受到的阻力迅速下降。原因是沖擊波與飛行器頭部弓形激波相互作用,改變激波結(jié)構(gòu),形成較弱的斜激波,同時(shí)伴隨沖擊波波陣面強(qiáng)度的衰減,飛行器受到的阻力迅速減小。

    在t=14μs時(shí),沖擊波傳播到飛行器外側(cè),同時(shí)沖擊波經(jīng)外壁面發(fā)射,在飛行器頭部形成低壓區(qū),使得阻力下降趨勢變緩。在t=21μs時(shí)達(dá)到最小值,僅為0.53 N。

    圖3 不同時(shí)間壓強(qiáng)等值線Fig.3 Pressure contours at the different time

    在t=22μs時(shí),阻力迅速升高,原因是受空氣來流影響,沖擊波前端波陣面及反射波陣面被吹向下游,并傳播到飛行器外壁面,產(chǎn)生作用力,導(dǎo)致阻力瞬時(shí)升高。

    在t=23μs瞬時(shí)達(dá)到31 N后,伴隨沖擊波強(qiáng)度的衰減,激光引致的沖擊波對飛行器的影響幾乎消失,在自由空氣來流的作用下,阻力緩慢升高,在t=36μs時(shí)達(dá)到第二阻力峰值,隨后逐漸減小,趨于穩(wěn)態(tài)。

    3.2 來流馬赫數(shù)的影響

    在相同激光參數(shù)及聚焦點(diǎn)位置條件下,計(jì)算比較分析來流馬赫數(shù)Ma=5~9范圍內(nèi)高超音速來流對激光等離子體減阻性能的影響。計(jì)算結(jié)果如表2所示。

    表2 不同來流條件下阻力減小情況Tab.2 Air drag reduction with different inflow

    從表中可以看出,隨著來流馬赫數(shù)的增大,飛行器受到的穩(wěn)態(tài)阻力快速增大,激光作用后的最小阻力增大,但減小的阻力也快速增大。

    圖4給出了不同來流馬赫數(shù)條件下飛行器受到的阻力隨時(shí)間變化曲線。

    圖4 不同來流條件下阻力曲線Fig.4 Air drag curveswith different Mach number

    從圖中可知,穩(wěn)態(tài)條件下,飛行器受到的阻力隨著來流馬赫數(shù)的增大明顯增大。在納秒激光作用下,隨著來流馬赫數(shù)的增大,空氣來流加速等離子體沖擊波吹向下游,使得沖擊波到達(dá)飛行器外壁面的時(shí)間提前,阻力峰值顯著增加,最小阻力增大。從低阻力維持時(shí)間看,隨來流馬赫數(shù)的增大,等離子體沖擊波的作用時(shí)間減小,導(dǎo)致低阻力維持時(shí)間縮短,阻力收斂時(shí)間提前。

    圖5給出了減小阻力的百分比ηD與來流馬赫數(shù)的關(guān)系。

    圖5 減小阻力的百分比隨來流馬赫數(shù)的變化關(guān)系Fig.5 The percentage of air drag reduction vs.Mach number

    從圖中擬合曲線可以看出,隨著來流馬赫數(shù)的增大,阻力減小的百分比呈下降趨勢。

    3.3 納秒激光能量的影響

    以來流馬赫數(shù)Ma=5.5為例,假設(shè)激光聚焦位置L=0.6D,分析納秒激光能量為0~2J范圍內(nèi)對減阻性能的影響。圖6給出了不同激光能量條件下,飛行器受到的阻力隨時(shí)間變化曲線。

    圖6 不同激光能量條件下阻力曲線Fig.6 Air drag curves with different laser energy

    從圖中可以看出,不同激光能量條件下,阻力曲線的變化規(guī)律是一致的。區(qū)別在于,隨著激光能量的增加,激光引致的沖擊波強(qiáng)度增大,沖擊波到達(dá)飛行器外壁面的時(shí)間縮短,阻力峰值增大,阻力最小值減小,低阻力維持的時(shí)間延長,阻力收斂時(shí)間推遲。

    圖7 減小阻力的百分比隨激光能量的變化關(guān)系Fig.7 The percentage of air drag reduction vs.laser enengy

    圖7 給出了減小阻力的百分比ηD與納秒激光能量之間的變化關(guān)系。

    從圖中可以看出,在納秒激光能量E<0.6J時(shí),減小阻力的百分比ηD隨著激光能量的增加而快速增加;而當(dāng)激光能量E>0.6J時(shí),減小阻力的百分比ηD達(dá)到97%以上,變化緩慢。

    3.4 聚焦位置的影響

    以來流馬赫數(shù)Ma=5.5為例,采用激光能量為E=2 J,脈沖寬度為tp=20 ns的納秒激光,分析激光聚焦擊穿位置對減阻性能的影響。圖8給出了不同點(diǎn)火位置的飛行器受到的阻力隨時(shí)間變化曲線。

    圖8 不同激光聚焦位置條件下阻力曲線Fig.8 Air drag curveswith different focusing location

    從圖中可知,相同來流環(huán)境及激光參數(shù)條件下,聚焦位置離飛行器頭部越近,等離子體沖擊

    圖9 減小阻力的百分比隨激光聚焦位置的變化Fig.9 The percentage of air drag reduction vs.laser focusing location

    波到達(dá)飛行器外壁面的時(shí)間越早,達(dá)到的阻力峰值越大,最小值越小。但從低阻力持續(xù)時(shí)間上看,聚焦位置越遠(yuǎn),低阻力維持時(shí)間越長。

    圖9給出了減小阻力的百分比隨點(diǎn)火位置的變化規(guī)律。

    從圖中擬合曲線可知,隨著激光聚焦位置的增大,減小阻力的百分比ηD減小,當(dāng)L/D<0.8時(shí),ηD變化較為緩慢;當(dāng)L/D>0.8時(shí),ηD快速下降。

    4 結(jié) 論

    以實(shí)驗(yàn)室已有的飛行器為例,采用納秒激光,開展激光等離子體減阻數(shù)值模擬仿真,簡要分析了來流馬赫數(shù)、脈沖激光能量和激光聚焦點(diǎn)對減小阻力的百分比的影響。主要結(jié)論如下:

    (1)納秒激光等離子體能夠很大程度提高減小阻力的百分比,在來流馬赫數(shù)為Ma=5.5條件下,達(dá)99.08%,并能夠增加低阻力維持的時(shí)間。

    (2)隨著來流馬赫數(shù)的增大,空氣來流對飛行器的阻力增大,對沖擊波的作用增大,導(dǎo)致阻力最大值升高,將沖擊波加速吹向下游,使得低阻力維持時(shí)間縮短,減小阻力的百分比減小。

    (3)隨著納秒激光能量的增加,激光引致的沖擊波強(qiáng)度增大,阻力最小值減小,減小阻力的百分比ηD增加,當(dāng)E<0.6J時(shí),ηD迅速增加,而當(dāng)E>0.6J時(shí),ηD達(dá)97%以上,變化緩慢。

    (4)隨著聚焦位置的增大,減小阻力的百分比ηD減小,但低阻力維持的時(shí)間明顯增加。

    [1] Myrabo.Solar-Powered global air transportation[R]. AIAA,1978:78-0689.

    [2] Sun Zhongxiang.Progress in plasma assisted drag reduction technology[J].Advances in Mechanics,2003,33(1):87-94.(in Chinese)孫宗祥.等離子體減阻技術(shù)的研究進(jìn)展[J].力學(xué)進(jìn)展,2003,33(1):87-94.

    [3] Hong Yanji,Li Qian,F(xiàn)ang Juan,et al.Advances in study of laser plasma drag reduction technology[J].Acta aeronautica et astronautica sinica,2010,31(1):93-101.(in Chinese)洪延姬,李倩,方娟,等.激光等離子體減阻技術(shù)研究進(jìn)展[J].航空學(xué)報(bào),2010,31(1):93-101.

    [4] Fang Juan,Hong Yanji,LiQian.Numerical analysis of interaction between single-pulse laser-induced plasma and bow shock in a supersonic flow[J].Plasma Science and Technology,2012,14(8):741-746.

    [5] David Sperber,F(xiàn)abian Schmid,Hans-Albert Eckel.Wave drag reduction of blunt bodies using laser sustained energy deposition in argon atmosphere[C].6thAIAA Flow Control Conference,2012:2012-2815.

    [6] HyoungJin Lee,InSeuck Jeung,SangHun Lee,et al.Experiment of flow control using laser energy deposition around high speed propulsion system[C].Beamed Energy Propulsion AIPConf.Proc.,2011,1402:416-423.

    [7] Fang Juan,Hong Yanji,Li Qian,et al.The influence of flight altitude on supersonic drag reduction with laser energy depositions[C].Beamed Energy Propulsion AIP Conf.Proc.,2011,1402:430-436.

    [8] Li Qian,Jin Xing,Cao Zhengrui,et al.Effects on aerodynamic drag of incident laser energy in technology of laser plasma point source drag reduction[J].Journal of propulsion technology,2010,31(3):377-380.(in Chinese)李倩,金星,曹正蕊等.激光等離子體點(diǎn)源減阻技術(shù)中入射能量對氣動阻力的影響[J].推進(jìn)技術(shù),2010,31(3):377-380.

    [9] Fang Juan,Hong Yanji,Huang Hui,et al.Influence of ignition location on drag reduction by laser plasma[J]. High power laser and particle beams,2010,22(9):2059-2062.(in Chinese)方娟,洪延姬,黃輝等.點(diǎn)火位置對激光等離子體減阻效能的影響[J].強(qiáng)激光與粒子束,2010,22(9):2059-2062.

    [10]Liu Zhun,Dou Zhiguo,Huang Hui,et al.The influence of the inflow mach number on laser plasma drag reduction[J].Development&Innovation ofMachinery&Electrical Products,2010,23(4):18-20.(in Chinese)劉準(zhǔn),竇志國,黃輝,等.來流馬赫數(shù)對激光等離子體減阻性能影響[J].機(jī)電產(chǎn)品開發(fā)與創(chuàng)新,2010,23(4):18-20.

    Numerical simulation of nanosecond laser plasma drag reduction

    CHEN Hao1,SHILei1,MA Li-hua1,XU Jing2,LU Chao-liang1,CHEN Yong-an1
    (1.School of Information and Navigation,Air Force Engineering University,Xi′an 710077,China;2.Army 95561 of PLA,Lhasa 850002,China)

    Aiming at the mechanism of laser plasma drag reduction technology,the propagation of the shock wave caused by the nanosecond laser is numerically simulated.The influence of the key parameters on the air drag reduction performance is analyzed.It can be concluded from the simulation results that the performance of the plasma drag reduction is improved greatly by the nanosecond laser,where the percentage of air drag reduction is close to 99%and the time of low air drag is103times of the incident laser time.The percentage of air drag reduction decreaseswith the Mach number increasing because the intensity of the air inflow increased.The percentage increases with the laser energy increasing because the shock wave caused by the nanosecond laser increases.The percentage decreaseswith the focusing location being far,but the time of low air drag obviously increases.

    nanosecond laser;plasma drag reduction;percentage of drag reduction;Mach number;laser energy;focusing location

    TN249;V439

    A

    10.3969/j.issn.1001-5078.2014.02.0

    1001-5078(2014)02-0131-05

    國家自然科學(xué)基金資助項(xiàng)目(No.61108068)。

    陳 豪(1989-),男,碩士研究生,主要從事激光等離子體減阻技術(shù)研究。E-mail:chenhao_8889@163.com

    2013-07-08;

    2013-07-27

    猜你喜歡
    納秒來流馬赫數(shù)
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    兩種典型來流條件下風(fēng)力機(jī)尾跡特性的數(shù)值研究
    能源工程(2022年2期)2022-05-23 13:51:48
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    載荷分布對可控?cái)U(kuò)散葉型性能的影響
    不同來流條件對溢洪道過流能力的影響
    亞納秒前沿脈沖強(qiáng)電場集成光學(xué)測試技術(shù)
    電子測試(2018年15期)2018-09-26 06:01:26
    基于高速數(shù)字電路的PD UHF信號納秒級陡脈沖源研制
    電測與儀表(2016年2期)2016-04-12 00:24:46
    空氣中納秒脈沖均勻DBD增加聚合物的表面親水性
    彈發(fā)匹配驗(yàn)證試驗(yàn)系統(tǒng)來流快速啟動技術(shù)研究
    玻璃氣體放電管與陶瓷氣體放電管的納秒脈沖響應(yīng)特性比較
    日韩av在线大香蕉| 天美传媒精品一区二区| 欧美性猛交黑人性爽| 黑人高潮一二区| 免费电影在线观看免费观看| 国产成年人精品一区二区| 丰满人妻一区二区三区视频av| 一进一出好大好爽视频| 黑人高潮一二区| 色综合站精品国产| 18禁在线播放成人免费| 偷拍熟女少妇极品色| 天天躁日日操中文字幕| 亚洲av电影不卡..在线观看| 淫妇啪啪啪对白视频| 国产伦在线观看视频一区| 亚洲成人av在线免费| 色在线成人网| 国产成人一区二区在线| 国内精品美女久久久久久| 国产av在哪里看| 亚洲一级一片aⅴ在线观看| 午夜亚洲福利在线播放| 国产单亲对白刺激| 中国美白少妇内射xxxbb| 一区二区三区高清视频在线| 麻豆国产97在线/欧美| 午夜老司机福利剧场| 欧美成人免费av一区二区三区| 国产激情偷乱视频一区二区| 99久久久亚洲精品蜜臀av| 国产蜜桃级精品一区二区三区| 日韩成人av中文字幕在线观看 | 欧美不卡视频在线免费观看| 精品人妻偷拍中文字幕| 国产一区二区在线av高清观看| 亚洲av熟女| 国产日本99.免费观看| 精品久久久噜噜| av福利片在线观看| 欧美日韩精品成人综合77777| 亚洲精品粉嫩美女一区| 久久久色成人| 一级毛片aaaaaa免费看小| 一级黄色大片毛片| 国产伦一二天堂av在线观看| 精品少妇黑人巨大在线播放 | 国产91av在线免费观看| 午夜福利成人在线免费观看| 国产av一区在线观看免费| 天堂网av新在线| 91在线观看av| 婷婷精品国产亚洲av| 99视频精品全部免费 在线| 亚洲av中文av极速乱| 精品久久久久久久久亚洲| 国产精品,欧美在线| 久久精品影院6| 日本免费a在线| 最后的刺客免费高清国语| 老司机福利观看| 特级一级黄色大片| 亚洲精品成人久久久久久| 中国国产av一级| 好男人在线观看高清免费视频| 又爽又黄a免费视频| 色综合色国产| 草草在线视频免费看| 69av精品久久久久久| 日本一二三区视频观看| 男人的好看免费观看在线视频| 99久久精品一区二区三区| 国产伦精品一区二区三区四那| 国产在视频线在精品| 校园春色视频在线观看| 国产av麻豆久久久久久久| 国产精品乱码一区二三区的特点| 婷婷亚洲欧美| 国产三级中文精品| 97超视频在线观看视频| 成人漫画全彩无遮挡| 日本爱情动作片www.在线观看 | 欧美3d第一页| 乱码一卡2卡4卡精品| 成人亚洲精品av一区二区| 午夜精品国产一区二区电影 | 五月玫瑰六月丁香| 欧美性感艳星| 成人无遮挡网站| 桃色一区二区三区在线观看| 乱码一卡2卡4卡精品| 国产精品99久久久久久久久| 高清午夜精品一区二区三区 | www日本黄色视频网| 亚洲综合色惰| 成人无遮挡网站| 99久国产av精品| 亚洲欧美成人精品一区二区| 亚洲av中文av极速乱| 久久久久久国产a免费观看| 1024手机看黄色片| 波多野结衣巨乳人妻| 亚洲一区二区三区色噜噜| 黑人高潮一二区| 最后的刺客免费高清国语| 日韩欧美国产在线观看| 色哟哟哟哟哟哟| 国产精品免费一区二区三区在线| 亚洲欧美日韩高清专用| 精品国内亚洲2022精品成人| 自拍偷自拍亚洲精品老妇| 亚州av有码| 在线观看午夜福利视频| 成人综合一区亚洲| 蜜臀久久99精品久久宅男| 欧美性猛交╳xxx乱大交人| 欧美xxxx性猛交bbbb| 国产欧美日韩精品一区二区| 村上凉子中文字幕在线| 亚洲va在线va天堂va国产| 十八禁国产超污无遮挡网站| ponron亚洲| 欧美一区二区亚洲| 国产成人aa在线观看| 中文在线观看免费www的网站| 国内精品一区二区在线观看| 欧美日本亚洲视频在线播放| 免费搜索国产男女视频| 成人欧美大片| 2021天堂中文幕一二区在线观| av中文乱码字幕在线| 狂野欧美激情性xxxx在线观看| 久久人人爽人人片av| 熟妇人妻久久中文字幕3abv| 亚洲美女视频黄频| 国产精品一区www在线观看| 亚洲精华国产精华液的使用体验 | 欧美另类亚洲清纯唯美| 99九九线精品视频在线观看视频| 搡老妇女老女人老熟妇| 成人无遮挡网站| 在现免费观看毛片| 国产麻豆成人av免费视频| 国产精品永久免费网站| 国产成人精品久久久久久| 国产一区二区亚洲精品在线观看| 99热这里只有精品一区| 永久网站在线| 最近中文字幕高清免费大全6| 亚洲美女搞黄在线观看 | 久久韩国三级中文字幕| 欧美最黄视频在线播放免费| 欧美日韩在线观看h| 国产精品久久久久久av不卡| 成人欧美大片| 精品国产三级普通话版| 简卡轻食公司| 特级一级黄色大片| 国产亚洲91精品色在线| 啦啦啦啦在线视频资源| 午夜福利在线观看免费完整高清在 | 久久亚洲精品不卡| 男人和女人高潮做爰伦理| 看黄色毛片网站| 黄片wwwwww| 国产亚洲精品久久久com| 小蜜桃在线观看免费完整版高清| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩高清专用| 国产精品一区二区三区四区免费观看 | 亚洲人与动物交配视频| 大香蕉久久网| 五月玫瑰六月丁香| 精品99又大又爽又粗少妇毛片| 99视频精品全部免费 在线| 女生性感内裤真人,穿戴方法视频| 成人二区视频| 99视频精品全部免费 在线| 淫秽高清视频在线观看| 久久人妻av系列| 久久精品人妻少妇| 97热精品久久久久久| a级一级毛片免费在线观看| 五月伊人婷婷丁香| 日韩 亚洲 欧美在线| 国产精品一区二区三区四区免费观看 | 如何舔出高潮| 国产免费一级a男人的天堂| 免费大片18禁| 一本精品99久久精品77| 18禁黄网站禁片免费观看直播| 国产精华一区二区三区| 久久久久国产精品人妻aⅴ院| 日本精品一区二区三区蜜桃| 国产单亲对白刺激| 久久热精品热| 日韩av在线大香蕉| 91狼人影院| 草草在线视频免费看| 久99久视频精品免费| 亚洲欧美日韩东京热| 国产精品人妻久久久影院| 日本黄色片子视频| 日韩欧美 国产精品| 国产高清有码在线观看视频| 国产伦精品一区二区三区四那| 一a级毛片在线观看| 春色校园在线视频观看| 精品日产1卡2卡| 亚洲人成网站高清观看| 久久精品综合一区二区三区| 在线观看午夜福利视频| 人妻少妇偷人精品九色| 一本精品99久久精品77| 在线观看免费视频日本深夜| 精品国产三级普通话版| 国产高清有码在线观看视频| 国产中年淑女户外野战色| 日日摸夜夜添夜夜添小说| 久久久久久九九精品二区国产| 高清日韩中文字幕在线| 久久久久久大精品| 亚洲四区av| 亚洲成a人片在线一区二区| 亚洲第一电影网av| 国产精品电影一区二区三区| av福利片在线观看| 亚洲美女视频黄频| 国产精品久久久久久av不卡| 欧洲精品卡2卡3卡4卡5卡区| 免费无遮挡裸体视频| 日本 av在线| 日日摸夜夜添夜夜爱| 久久精品国产亚洲av天美| 一级黄色大片毛片| 亚洲综合色惰| 日韩高清综合在线| 国产大屁股一区二区在线视频| 99视频精品全部免费 在线| 亚洲av免费在线观看| 亚洲精品国产av成人精品 | 不卡一级毛片| 天堂√8在线中文| 可以在线观看的亚洲视频| av在线老鸭窝| 欧美激情久久久久久爽电影| 蜜臀久久99精品久久宅男| 91av网一区二区| 欧美bdsm另类| 少妇裸体淫交视频免费看高清| 中文字幕精品亚洲无线码一区| 中文在线观看免费www的网站| 成年版毛片免费区| 国产 一区精品| 久久久久久大精品| 免费搜索国产男女视频| 午夜福利在线观看吧| 国产欧美日韩一区二区精品| 欧美国产日韩亚洲一区| 亚洲专区国产一区二区| 男女啪啪激烈高潮av片| 国产精品av视频在线免费观看| 成人漫画全彩无遮挡| 久久天躁狠狠躁夜夜2o2o| 色视频www国产| 不卡视频在线观看欧美| 精华霜和精华液先用哪个| 国产高清有码在线观看视频| 高清毛片免费看| 亚洲精品一卡2卡三卡4卡5卡| 久久久午夜欧美精品| 日韩三级伦理在线观看| 超碰av人人做人人爽久久| 免费观看的影片在线观看| 麻豆精品久久久久久蜜桃| 久久久欧美国产精品| 亚洲激情五月婷婷啪啪| 哪里可以看免费的av片| 国产美女午夜福利| 亚洲av中文字字幕乱码综合| 欧美日韩综合久久久久久| 一级a爱片免费观看的视频| 日韩强制内射视频| 久久人人爽人人片av| 国产av麻豆久久久久久久| 99久久久亚洲精品蜜臀av| 亚洲内射少妇av| 国产中年淑女户外野战色| 搞女人的毛片| 婷婷六月久久综合丁香| 亚洲七黄色美女视频| 亚洲国产精品国产精品| 寂寞人妻少妇视频99o| 亚洲无线在线观看| 日本黄色片子视频| 亚洲精品一卡2卡三卡4卡5卡| 国产精品人妻久久久影院| 美女 人体艺术 gogo| 日本熟妇午夜| 久久久久久久久大av| 国产欧美日韩一区二区精品| 在线国产一区二区在线| 22中文网久久字幕| 国产伦精品一区二区三区视频9| 岛国在线免费视频观看| a级毛色黄片| 不卡视频在线观看欧美| 亚洲欧美日韩无卡精品| 最近中文字幕高清免费大全6| 麻豆成人午夜福利视频| 国产欧美日韩精品一区二区| 国产探花极品一区二区| 小蜜桃在线观看免费完整版高清| 午夜福利在线观看免费完整高清在 | 亚洲性久久影院| 三级国产精品欧美在线观看| 午夜a级毛片| videossex国产| 99久国产av精品| 免费av观看视频| av在线观看视频网站免费| 丝袜美腿在线中文| 啦啦啦观看免费观看视频高清| 自拍偷自拍亚洲精品老妇| 国产精品久久视频播放| 中文字幕av在线有码专区| 一区福利在线观看| 国产真实乱freesex| 亚洲欧美日韩东京热| 久久人人爽人人爽人人片va| 久久精品国产99精品国产亚洲性色| 女同久久另类99精品国产91| 狂野欧美激情性xxxx在线观看| 99久久成人亚洲精品观看| 亚洲精品日韩在线中文字幕 | 国产一区亚洲一区在线观看| 毛片女人毛片| 中文在线观看免费www的网站| 国产高清不卡午夜福利| 无遮挡黄片免费观看| 网址你懂的国产日韩在线| 男女啪啪激烈高潮av片| 亚洲无线在线观看| 三级毛片av免费| 一级毛片aaaaaa免费看小| 国产一级毛片七仙女欲春2| 黄片wwwwww| 一边摸一边抽搐一进一小说| 亚洲国产精品久久男人天堂| 黄色一级大片看看| 黄片wwwwww| 久久午夜福利片| 欧美最黄视频在线播放免费| 少妇人妻一区二区三区视频| 亚洲va在线va天堂va国产| 精品久久久久久久人妻蜜臀av| 少妇丰满av| 好男人在线观看高清免费视频| 内射极品少妇av片p| 日韩欧美一区二区三区在线观看| 三级毛片av免费| 欧美日本视频| 99国产极品粉嫩在线观看| 国产男人的电影天堂91| 一本久久中文字幕| 成年女人毛片免费观看观看9| 看片在线看免费视频| 人人妻人人澡欧美一区二区| 精品久久久久久久久av| 久99久视频精品免费| 国产av麻豆久久久久久久| a级毛色黄片| 一级毛片久久久久久久久女| 亚洲七黄色美女视频| 免费看av在线观看网站| 欧美+日韩+精品| 久久国产乱子免费精品| av视频在线观看入口| 中出人妻视频一区二区| 黄色配什么色好看| 久久久成人免费电影| 免费观看的影片在线观看| 国产色爽女视频免费观看| 成年女人看的毛片在线观看| 亚洲美女视频黄频| 日本-黄色视频高清免费观看| 99久久中文字幕三级久久日本| 69av精品久久久久久| 国产 一区 欧美 日韩| 色哟哟哟哟哟哟| 一夜夜www| 青春草视频在线免费观看| 国产单亲对白刺激| av黄色大香蕉| 永久网站在线| 亚洲av免费在线观看| 久久国内精品自在自线图片| 日本在线视频免费播放| 国产爱豆传媒在线观看| 在线观看午夜福利视频| 直男gayav资源| 老熟妇仑乱视频hdxx| 乱人视频在线观看| 亚洲国产精品国产精品| 国产女主播在线喷水免费视频网站 | 五月伊人婷婷丁香| 国产成人aa在线观看| 在线观看66精品国产| 亚洲av一区综合| 免费人成在线观看视频色| 久久精品国产亚洲av香蕉五月| 亚洲成a人片在线一区二区| 国内精品美女久久久久久| 日日干狠狠操夜夜爽| 看黄色毛片网站| 亚洲av美国av| 欧美另类亚洲清纯唯美| 亚洲精品在线观看二区| 久久久色成人| 亚洲经典国产精华液单| 中国国产av一级| 免费人成在线观看视频色| 香蕉av资源在线| 国产69精品久久久久777片| 欧美最新免费一区二区三区| 亚洲国产欧美人成| 久久久色成人| 91在线观看av| 搡老岳熟女国产| 国产淫片久久久久久久久| 国产午夜精品论理片| 国产伦精品一区二区三区视频9| 最近2019中文字幕mv第一页| 97超碰精品成人国产| 九九久久精品国产亚洲av麻豆| 国产黄色视频一区二区在线观看 | 一本一本综合久久| 亚洲无线观看免费| 免费电影在线观看免费观看| 51国产日韩欧美| 少妇高潮的动态图| 精品一区二区免费观看| 久久久色成人| 国内精品一区二区在线观看| 亚洲18禁久久av| 最新在线观看一区二区三区| av福利片在线观看| av在线天堂中文字幕| 乱码一卡2卡4卡精品| 久久国产乱子免费精品| 亚洲欧美日韩高清专用| 免费av不卡在线播放| 搡老妇女老女人老熟妇| 日韩欧美在线乱码| 亚洲av成人av| 老师上课跳d突然被开到最大视频| 日本与韩国留学比较| 最后的刺客免费高清国语| 国产探花在线观看一区二区| 亚洲精品日韩av片在线观看| 午夜免费激情av| 少妇的逼水好多| 日本 av在线| 99热6这里只有精品| 在线a可以看的网站| 俄罗斯特黄特色一大片| 国产精品,欧美在线| 性色avwww在线观看| 日韩大尺度精品在线看网址| 精品久久久久久久久亚洲| 三级毛片av免费| 在线天堂最新版资源| 亚洲av熟女| 国产伦在线观看视频一区| 一级毛片电影观看 | 日本与韩国留学比较| 国产成人aa在线观看| av天堂中文字幕网| 欧美日韩乱码在线| 成人一区二区视频在线观看| 日韩欧美一区二区三区在线观看| 能在线免费观看的黄片| 亚洲第一区二区三区不卡| 亚洲欧美日韩卡通动漫| 免费人成在线观看视频色| 久久婷婷人人爽人人干人人爱| 热99re8久久精品国产| 全区人妻精品视频| 午夜亚洲福利在线播放| 波多野结衣巨乳人妻| 哪里可以看免费的av片| 国产成人a区在线观看| 免费电影在线观看免费观看| 别揉我奶头 嗯啊视频| 国产一区二区在线观看日韩| 国产淫片久久久久久久久| 欧美激情久久久久久爽电影| 久久久久久大精品| 午夜久久久久精精品| 久久这里只有精品中国| 别揉我奶头~嗯~啊~动态视频| 在线天堂最新版资源| 毛片一级片免费看久久久久| 大香蕉久久网| 亚洲aⅴ乱码一区二区在线播放| 伦精品一区二区三区| a级毛片a级免费在线| 99精品在免费线老司机午夜| 亚州av有码| 校园春色视频在线观看| 一区二区三区高清视频在线| 嫩草影院入口| 国产精品99久久久久久久久| 日韩强制内射视频| 天堂网av新在线| 91麻豆精品激情在线观看国产| 99热这里只有是精品在线观看| 亚洲一区二区三区色噜噜| 男人的好看免费观看在线视频| 熟女电影av网| 插阴视频在线观看视频| 99视频精品全部免费 在线| 一区二区三区免费毛片| 黄色日韩在线| 97热精品久久久久久| 在线观看免费视频日本深夜| 亚洲丝袜综合中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 久久人人爽人人片av| 国产爱豆传媒在线观看| 一本久久中文字幕| 国产精品乱码一区二三区的特点| 久久久久久国产a免费观看| 日韩大尺度精品在线看网址| 久久久久久久久大av| av黄色大香蕉| 亚洲av.av天堂| 99久久成人亚洲精品观看| 久久久成人免费电影| 亚洲第一电影网av| 精品国产三级普通话版| 成人亚洲精品av一区二区| 99国产极品粉嫩在线观看| 好男人在线观看高清免费视频| av福利片在线观看| 精品免费久久久久久久清纯| 国产成人一区二区在线| 日本免费a在线| 欧美激情国产日韩精品一区| 亚洲七黄色美女视频| 精品福利观看| av在线播放精品| 亚洲av二区三区四区| 99riav亚洲国产免费| 国产亚洲精品久久久久久毛片| 给我免费播放毛片高清在线观看| av在线蜜桃| 99热6这里只有精品| 欧美日韩在线观看h| 国产真实乱freesex| 欧美+日韩+精品| 丝袜美腿在线中文| 婷婷精品国产亚洲av| 午夜免费激情av| 全区人妻精品视频| 美女内射精品一级片tv| 免费无遮挡裸体视频| 岛国在线免费视频观看| 亚洲色图av天堂| 亚洲七黄色美女视频| 伦精品一区二区三区| 国产精品女同一区二区软件| 特级一级黄色大片| 日日干狠狠操夜夜爽| h日本视频在线播放| 亚洲精品久久国产高清桃花| 亚洲电影在线观看av| 麻豆成人午夜福利视频| 久久欧美精品欧美久久欧美| 一级av片app| 日本 av在线| 国产在线精品亚洲第一网站| 天堂√8在线中文| 亚洲精华国产精华液的使用体验 | 精品久久国产蜜桃| 久久久欧美国产精品| 欧美一区二区亚洲| 免费av观看视频| 偷拍熟女少妇极品色| 欧美激情国产日韩精品一区| 免费看美女性在线毛片视频| 亚洲人成网站在线观看播放| 最近视频中文字幕2019在线8| 国产伦一二天堂av在线观看| 久久6这里有精品| 国产亚洲精品久久久久久毛片| 亚洲最大成人手机在线| 亚洲一区二区三区色噜噜| 国产黄色视频一区二区在线观看 | 最后的刺客免费高清国语| 人妻夜夜爽99麻豆av| 午夜激情福利司机影院| 久久国产乱子免费精品| 精品人妻一区二区三区麻豆 | 在线免费观看不下载黄p国产| 国产一区二区三区av在线 | 一级毛片aaaaaa免费看小| 黄色视频,在线免费观看| 久久久久国产精品人妻aⅴ院| 在线观看午夜福利视频| 国产精品亚洲一级av第二区| 久久九九热精品免费| 日韩欧美 国产精品| 欧美xxxx性猛交bbbb| 欧美日本亚洲视频在线播放|