• <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)特性比較
    汤姆久久久久久久影院中文字幕 | 搡老妇女老女人老熟妇| 男人的好看免费观看在线视频| 亚洲高清免费不卡视频| 18禁裸乳无遮挡免费网站照片| 18禁动态无遮挡网站| 少妇人妻精品综合一区二区| 老女人水多毛片| 伦理电影大哥的女人| 成人毛片a级毛片在线播放| 波多野结衣高清无吗| 亚洲av福利一区| 成人无遮挡网站| 久久久久久久国产电影| 国产色婷婷99| 久久久久久九九精品二区国产| 国产一区亚洲一区在线观看| 18禁在线无遮挡免费观看视频| av国产久精品久网站免费入址| 中文字幕制服av| 欧美激情国产日韩精品一区| 亚洲精品成人久久久久久| 深爱激情五月婷婷| 免费在线观看成人毛片| 欧美日韩国产亚洲二区| 国产精品综合久久久久久久免费| 国产精品日韩av在线免费观看| 国产精品熟女久久久久浪| 天堂影院成人在线观看| 亚洲欧美精品自产自拍| 日本猛色少妇xxxxx猛交久久| 在现免费观看毛片| 亚洲丝袜综合中文字幕| 欧美又色又爽又黄视频| 91久久精品国产一区二区三区| 日韩av在线大香蕉| 一边摸一边抽搐一进一小说| 久久99热这里只频精品6学生 | av线在线观看网站| 国产成人a区在线观看| 亚洲无线观看免费| 欧美日韩一区二区视频在线观看视频在线 | 在线观看66精品国产| 国产一级毛片七仙女欲春2| 热99在线观看视频| 丰满乱子伦码专区| 日韩欧美精品免费久久| 免费黄色在线免费观看| 别揉我奶头 嗯啊视频| 黄色欧美视频在线观看| 国产精品.久久久| 国内少妇人妻偷人精品xxx网站| 亚洲精品久久久久久婷婷小说 | 久久久色成人| 建设人人有责人人尽责人人享有的 | 欧美xxxx黑人xx丫x性爽| 美女cb高潮喷水在线观看| 简卡轻食公司| 国产一级毛片在线| 最近的中文字幕免费完整| 晚上一个人看的免费电影| 国产色婷婷99| 亚洲在久久综合| 亚洲中文字幕日韩| 26uuu在线亚洲综合色| 国产成人91sexporn| 淫秽高清视频在线观看| 国产免费视频播放在线视频 | 国国产精品蜜臀av免费| 亚洲内射少妇av| 美女xxoo啪啪120秒动态图| 一夜夜www| 中文字幕熟女人妻在线| 亚洲综合精品二区| 亚洲18禁久久av| 久久精品夜夜夜夜夜久久蜜豆| 欧美一区二区国产精品久久精品| 一区二区三区乱码不卡18| 欧美成人一区二区免费高清观看| 女人被狂操c到高潮| 久久人妻av系列| 高清视频免费观看一区二区 | 久久久久免费精品人妻一区二区| 国产又黄又爽又无遮挡在线| 亚洲美女搞黄在线观看| 黄色欧美视频在线观看| 在线播放无遮挡| 国产av码专区亚洲av| 永久网站在线| 久久久午夜欧美精品| 黑人高潮一二区| 欧美激情在线99| 99热这里只有是精品在线观看| 精品国内亚洲2022精品成人| 一卡2卡三卡四卡精品乱码亚洲| 极品教师在线视频| 国产91av在线免费观看| 成年av动漫网址| 99久久成人亚洲精品观看| 亚洲欧美成人精品一区二区| 在线a可以看的网站| 黄片无遮挡物在线观看| 亚洲18禁久久av| 18禁在线无遮挡免费观看视频| 久久人人爽人人片av| 中文字幕久久专区| 午夜免费男女啪啪视频观看| 国产熟女欧美一区二区| 久久热精品热| 日韩制服骚丝袜av| 菩萨蛮人人尽说江南好唐韦庄 | 中国国产av一级| 青春草亚洲视频在线观看| 嫩草影院新地址| 亚洲成人久久爱视频| 日韩一区二区三区影片| 99久久成人亚洲精品观看| 国产成人freesex在线| 美女脱内裤让男人舔精品视频| 身体一侧抽搐| 亚洲成人精品中文字幕电影| 菩萨蛮人人尽说江南好唐韦庄 | 国产真实乱freesex| 蜜桃亚洲精品一区二区三区| 国产伦理片在线播放av一区| 午夜福利在线观看免费完整高清在| 男人舔女人下体高潮全视频| 中国国产av一级| 亚洲欧美日韩高清专用| 国产精品久久电影中文字幕| 国产av一区在线观看免费| 日本一本二区三区精品| 啦啦啦韩国在线观看视频| 黄色一级大片看看| 蜜臀久久99精品久久宅男| 亚洲精品色激情综合| 欧美性猛交╳xxx乱大交人| 2021少妇久久久久久久久久久| 午夜激情福利司机影院| 亚洲美女搞黄在线观看| 精品久久久久久久人妻蜜臀av| 精品一区二区三区人妻视频| 久久精品久久久久久噜噜老黄 | 91精品一卡2卡3卡4卡| 亚洲电影在线观看av| 亚洲成人中文字幕在线播放| 日韩三级伦理在线观看| 久久综合国产亚洲精品| 国产免费视频播放在线视频 | 中文字幕精品亚洲无线码一区| 老司机福利观看| 精品久久久久久久人妻蜜臀av| 精品不卡国产一区二区三区| 亚洲美女视频黄频| 国产精品女同一区二区软件| 欧美色视频一区免费| 乱系列少妇在线播放| 国产乱来视频区| 亚洲不卡免费看| 可以在线观看毛片的网站| 麻豆成人av视频| 精品酒店卫生间| 男人舔奶头视频| 99热网站在线观看| av在线播放精品| 国产高清有码在线观看视频| 亚洲va在线va天堂va国产| 秋霞在线观看毛片| 老女人水多毛片| 久久久久久久久久久免费av| 男人舔女人下体高潮全视频| 久久99热这里只频精品6学生 | 亚洲在久久综合| 国产 一区精品| 激情 狠狠 欧美| 亚洲最大成人手机在线| 中文字幕人妻熟人妻熟丝袜美| 最近中文字幕2019免费版| 成人国产麻豆网| 国产高清有码在线观看视频| 久久草成人影院| 国产成人freesex在线| 欧美日韩精品成人综合77777| 国产一区二区在线av高清观看| 午夜老司机福利剧场| 国产综合懂色| 国产熟女欧美一区二区| 亚洲精品色激情综合| 亚洲国产精品成人久久小说| 秋霞伦理黄片| 午夜老司机福利剧场| 中文字幕免费在线视频6| 内射极品少妇av片p| av在线观看视频网站免费| 欧美一区二区国产精品久久精品| 国产片特级美女逼逼视频| 日韩欧美国产在线观看| 岛国毛片在线播放| 亚洲国产精品sss在线观看| 偷拍熟女少妇极品色| 国产麻豆成人av免费视频| 国产精品电影一区二区三区| 午夜福利网站1000一区二区三区| 在线免费观看的www视频| 午夜精品在线福利| 大又大粗又爽又黄少妇毛片口| 亚洲av日韩在线播放| 亚洲熟妇中文字幕五十中出| 国内精品宾馆在线| 亚洲三级黄色毛片| 亚洲av日韩在线播放| 久久久久久伊人网av| 国语对白做爰xxxⅹ性视频网站| 亚洲人成网站高清观看| 日本黄色视频三级网站网址| 国产欧美另类精品又又久久亚洲欧美| 在现免费观看毛片| 麻豆一二三区av精品| 自拍偷自拍亚洲精品老妇| 色播亚洲综合网| 欧美高清性xxxxhd video| 熟女电影av网| 精品久久国产蜜桃| 久久精品国产亚洲网站| 亚洲成人久久爱视频| 一卡2卡三卡四卡精品乱码亚洲| 国产精品伦人一区二区| 久久精品综合一区二区三区| 三级毛片av免费| 尾随美女入室| 成人三级黄色视频| av在线播放精品| 国产综合懂色| 汤姆久久久久久久影院中文字幕 | 午夜视频国产福利| 我的老师免费观看完整版| 久久久精品94久久精品| 免费看av在线观看网站| 一个人看的www免费观看视频| 日韩欧美三级三区| 在线免费十八禁| 久久鲁丝午夜福利片| 99久久中文字幕三级久久日本| 99久久无色码亚洲精品果冻| 亚洲人成网站在线播| 日本黄色视频三级网站网址| 亚洲国产精品专区欧美| 久久99热这里只频精品6学生 | 大话2 男鬼变身卡| 亚洲精品国产成人久久av| 欧美xxxx性猛交bbbb| 色综合色国产| 五月玫瑰六月丁香| 日韩亚洲欧美综合| 老女人水多毛片| h日本视频在线播放| 黄片wwwwww| 国产黄片视频在线免费观看| 女人被狂操c到高潮| 国产白丝娇喘喷水9色精品| 精品久久久久久久久亚洲| 欧美性感艳星| 蜜桃亚洲精品一区二区三区| 日韩在线高清观看一区二区三区| 欧美又色又爽又黄视频| 综合色丁香网| 亚洲国产成人一精品久久久| 国产毛片a区久久久久| 日韩制服骚丝袜av| 午夜福利视频1000在线观看| 又黄又爽又刺激的免费视频.| 国产午夜福利久久久久久| 少妇的逼好多水| 欧美丝袜亚洲另类| 亚洲性久久影院| 少妇丰满av| 性色avwww在线观看| 黄色欧美视频在线观看| 午夜福利在线在线| 69人妻影院| 亚洲图色成人| 你懂的网址亚洲精品在线观看 | 插逼视频在线观看| 好男人视频免费观看在线| 欧美激情国产日韩精品一区| 18禁在线播放成人免费| 看黄色毛片网站| 黑人高潮一二区| av线在线观看网站| 亚洲av成人精品一二三区| 国产视频内射| 亚洲精品久久久久久婷婷小说 | 国产老妇女一区| 久久6这里有精品| 久久亚洲国产成人精品v| 精品一区二区免费观看| 久久精品国产亚洲av天美| 国产精品无大码| 亚洲人与动物交配视频| 久久久亚洲精品成人影院| 国内精品宾馆在线| 国产极品精品免费视频能看的| 美女内射精品一级片tv| 国产白丝娇喘喷水9色精品| 观看美女的网站| 少妇的逼水好多| 亚洲av男天堂| videos熟女内射| 久久精品国产自在天天线| 国产69精品久久久久777片| 中文字幕av成人在线电影| 一本一本综合久久| 欧美三级亚洲精品| 午夜免费激情av| 人人妻人人澡人人爽人人夜夜 | 一本一本综合久久| 国内精品美女久久久久久| 亚洲国产最新在线播放| 99视频精品全部免费 在线| 99久久无色码亚洲精品果冻| 国产黄片美女视频| 丰满乱子伦码专区| 在线播放国产精品三级| 亚洲精品自拍成人| 久久综合国产亚洲精品| 欧美另类亚洲清纯唯美| 久久久久国产网址| 亚洲精品乱久久久久久| 极品教师在线视频| 人人妻人人澡人人爽人人夜夜 | 国内精品一区二区在线观看| 在线观看av片永久免费下载| 精品不卡国产一区二区三区| 久久亚洲精品不卡| 国产精品三级大全| 高清毛片免费看| 日本-黄色视频高清免费观看| 日韩强制内射视频| 欧美高清性xxxxhd video| 午夜福利高清视频| 亚洲av熟女| 久久精品熟女亚洲av麻豆精品 | 国产又色又爽无遮挡免| 自拍偷自拍亚洲精品老妇| 两个人视频免费观看高清| 女人被狂操c到高潮| 成人亚洲欧美一区二区av| 亚洲国产精品专区欧美| 欧美xxxx黑人xx丫x性爽| 日韩在线高清观看一区二区三区| 亚洲国产精品专区欧美| 国产乱来视频区| 插逼视频在线观看| 免费看美女性在线毛片视频| 十八禁国产超污无遮挡网站| 国产一区二区三区av在线| 午夜精品在线福利| 国内精品一区二区在线观看| 国产一区二区在线av高清观看| 99热这里只有精品一区| 日韩一本色道免费dvd| 国产成人a∨麻豆精品| 看黄色毛片网站| 色吧在线观看| 又爽又黄无遮挡网站| 日本黄大片高清| 99国产精品一区二区蜜桃av| 国产精品电影一区二区三区| 久久亚洲精品不卡| 日韩一区二区三区影片| 国产成人精品久久久久久| 日日干狠狠操夜夜爽| 人人妻人人看人人澡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | .国产精品久久| 日本爱情动作片www.在线观看| 久久久欧美国产精品| 国产一区亚洲一区在线观看| 亚洲精品456在线播放app| 寂寞人妻少妇视频99o| 男女视频在线观看网站免费| 久久精品熟女亚洲av麻豆精品 | 精品一区二区三区人妻视频| 欧美zozozo另类| 亚洲欧美成人精品一区二区| 精华霜和精华液先用哪个| 简卡轻食公司| 色网站视频免费| 免费观看人在逋| 国产日韩欧美在线精品| 国产真实伦视频高清在线观看| 99久久精品国产国产毛片| 色综合色国产| 乱人视频在线观看| 久久久久久久亚洲中文字幕| 国产不卡一卡二| 日韩一区二区三区影片| 久久精品91蜜桃| 日韩av在线大香蕉| kizo精华| 亚洲精品乱久久久久久| 国产91av在线免费观看| 日本五十路高清| 久久久久网色| 国产精品1区2区在线观看.| 18禁在线播放成人免费| 亚洲人成网站高清观看| 久久精品夜夜夜夜夜久久蜜豆| 麻豆乱淫一区二区| 天堂中文最新版在线下载 | 伦精品一区二区三区| 成年女人永久免费观看视频| av在线老鸭窝| 亚洲怡红院男人天堂| 国产一级毛片在线| 黄片无遮挡物在线观看| 国产中年淑女户外野战色| 国产精品久久久久久久电影| 久久国内精品自在自线图片| 国产免费福利视频在线观看| 久久人人爽人人片av| 亚洲乱码一区二区免费版| 国产精品久久久久久久久免| 久久久久久伊人网av| 天堂网av新在线| 色综合站精品国产| 欧美+日韩+精品| 欧美zozozo另类| 99视频精品全部免费 在线| 精品欧美国产一区二区三| 国产高潮美女av| 麻豆乱淫一区二区| 最新中文字幕久久久久| 国产一区有黄有色的免费视频 | av视频在线观看入口| 一个人看的www免费观看视频| 亚洲欧美清纯卡通| 又粗又爽又猛毛片免费看| 免费看av在线观看网站| 久久久成人免费电影| 天堂网av新在线| 久久久久久大精品| 欧美激情久久久久久爽电影| 一区二区三区乱码不卡18| 边亲边吃奶的免费视频| 三级国产精品片| 日韩中字成人| 国产 一区 欧美 日韩| 国产成人福利小说| 亚洲在久久综合| 一区二区三区免费毛片| 成人毛片60女人毛片免费| 在线免费十八禁| 国产精品一区二区性色av| 美女cb高潮喷水在线观看| 老司机影院毛片| 天天躁日日操中文字幕| 内射极品少妇av片p| 亚洲欧洲日产国产| 午夜久久久久精精品| 看非洲黑人一级黄片| 国产一级毛片在线| 久久精品91蜜桃| 日本猛色少妇xxxxx猛交久久| 国产亚洲av片在线观看秒播厂 | 欧美成人精品欧美一级黄| 久久久亚洲精品成人影院| 超碰97精品在线观看| 九九久久精品国产亚洲av麻豆| 久久久a久久爽久久v久久| 男女边吃奶边做爰视频| 看免费成人av毛片| 亚洲av免费高清在线观看| 亚洲av中文av极速乱| 18禁动态无遮挡网站| 久久久久久久久大av| 国产精品野战在线观看| 两个人视频免费观看高清| 全区人妻精品视频| 午夜激情欧美在线| 看非洲黑人一级黄片| 国产高清三级在线| 深爱激情五月婷婷| 精品国产露脸久久av麻豆 | 午夜精品国产一区二区电影 | 久久久久久九九精品二区国产| 少妇熟女欧美另类| 搡女人真爽免费视频火全软件| 久久久久久久亚洲中文字幕| 舔av片在线| 亚洲成色77777| 99在线视频只有这里精品首页| 午夜福利在线在线| 亚洲国产精品国产精品| 丰满乱子伦码专区| 一个人看的www免费观看视频| 久久久国产成人免费| av在线观看视频网站免费| 晚上一个人看的免费电影| 免费观看a级毛片全部| 久久精品国产99精品国产亚洲性色| 草草在线视频免费看| 尤物成人国产欧美一区二区三区| 国产私拍福利视频在线观看| av在线老鸭窝| 中国国产av一级| 免费av不卡在线播放| 能在线免费观看的黄片| 午夜a级毛片| 99热精品在线国产| 99国产精品一区二区蜜桃av| 亚洲精品aⅴ在线观看| 免费看日本二区| 亚洲av免费在线观看| 亚洲无线观看免费| 美女cb高潮喷水在线观看| 国产亚洲5aaaaa淫片| 国产乱人偷精品视频| 国产精品爽爽va在线观看网站| 久久久久久久久久久免费av| 最近2019中文字幕mv第一页| 午夜a级毛片| av专区在线播放| 91久久精品国产一区二区成人| 日本黄大片高清| 久久草成人影院| 亚洲自偷自拍三级| 久久精品国产亚洲av涩爱| 国产在视频线精品| 亚洲乱码一区二区免费版| 国产亚洲av片在线观看秒播厂 | 国产精品久久视频播放| 欧美又色又爽又黄视频| 日韩欧美国产在线观看| 久久综合国产亚洲精品| 天堂中文最新版在线下载 | 99国产精品一区二区蜜桃av| 91久久精品电影网| 一级黄片播放器| 国产午夜精品久久久久久一区二区三区| 建设人人有责人人尽责人人享有的 | 国产真实乱freesex| 国产爱豆传媒在线观看| 免费黄色在线免费观看| 日本与韩国留学比较| 如何舔出高潮| 99热网站在线观看| 国产日韩欧美在线精品| 精品久久国产蜜桃| kizo精华| 精品免费久久久久久久清纯| 大香蕉久久网| 亚洲国产最新在线播放| 日本黄大片高清| 久久精品夜色国产| 久久99精品国语久久久| 久久草成人影院| 日韩制服骚丝袜av| 国产成人a区在线观看| 色哟哟·www| 国产午夜精品论理片| 一区二区三区四区激情视频| 免费av不卡在线播放| 国内精品宾馆在线| 精品人妻偷拍中文字幕| 国产精品久久电影中文字幕| 国产精品爽爽va在线观看网站| 人体艺术视频欧美日本| 婷婷色av中文字幕| 床上黄色一级片| 欧美最新免费一区二区三区| 国产精品国产三级国产专区5o | 在线播放无遮挡| 99久久精品国产国产毛片| 五月伊人婷婷丁香| 国产伦精品一区二区三区四那| 久久午夜福利片| 中文欧美无线码| 亚洲三级黄色毛片| 亚洲精品影视一区二区三区av| 国产av在哪里看| 天堂网av新在线| 99久久无色码亚洲精品果冻| 亚洲精品久久久久久婷婷小说 | 99久久精品热视频| 国产伦一二天堂av在线观看| 亚洲国产色片| 日韩制服骚丝袜av| 国产精品99久久久久久久久| 成人av在线播放网站| 搡女人真爽免费视频火全软件| 91精品一卡2卡3卡4卡| 在现免费观看毛片| 亚洲美女视频黄频| 麻豆成人av视频| 日韩高清综合在线| 日韩av不卡免费在线播放| 国产老妇伦熟女老妇高清| 国产大屁股一区二区在线视频| 国产成人精品婷婷| 午夜免费男女啪啪视频观看| 精品免费久久久久久久清纯| 亚洲av二区三区四区| 成人性生交大片免费视频hd| 91aial.com中文字幕在线观看| 国产亚洲av片在线观看秒播厂 | 午夜福利在线在线| 久久精品国产亚洲av天美| 国产一区二区在线av高清观看| 黄色一级大片看看| 老司机影院成人| 自拍偷自拍亚洲精品老妇| 国产大屁股一区二区在线视频| 欧美性猛交╳xxx乱大交人| 99视频精品全部免费 在线|