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

    新一代四機并聯(lián)火箭發(fā)動機噴流熱環(huán)境數(shù)值研究

    2023-04-06 07:30:32趙一霖來霄毅馬祿創(chuàng)侯凌霄葉哲霄
    空天防御 2023年1期
    關(guān)鍵詞:箭體飛行高度噴流

    趙一霖,嚴 立,來霄毅,馬祿創(chuàng),侯凌霄,葉哲霄

    (上海宇航系統(tǒng)工程研究所,上海 201109)

    0 引 言

    隨著航天產(chǎn)業(yè)的快速發(fā)展,對運載火箭運載能力的要求越來越高,而采用多噴管火箭發(fā)動機是提高運載火箭運載能力的主要途徑之一。多噴管火箭發(fā)動機工作時會產(chǎn)生多股噴流相互干擾的現(xiàn)象,同時噴流的引射和相互干擾的共同作用形成回流,底部復(fù)雜的熱環(huán)境及流動特性與火箭的飛行參數(shù)、推進劑種類、發(fā)動機出口的燃氣參數(shù)、發(fā)動機噴管結(jié)構(gòu)參數(shù)及其排列方式等因素相關(guān),多因素的影響對運載火箭底部結(jié)構(gòu)防熱設(shè)計帶來風(fēng)險。

    目前國內(nèi)的現(xiàn)役常溫推進劑火箭發(fā)動機推力較小,其噴流溫度和速度也較小,同時發(fā)動機噴管尺寸較小,噴管之間的距離在500 mm 左右,這使得現(xiàn)役常溫推進劑的火箭底部最大熱流密度在200 kW/m2左右。現(xiàn)役新一代捆綁運載火箭助推器與芯級的間距較大,噴流的交互作用較弱,噴流輻射作用對底部加熱影響更大,火箭底部最大熱流密度在500 kW/m2左右。對比現(xiàn)役火箭,新一代四機并聯(lián)大推力液氧煤油運載火箭的特點為:發(fā)動機噴管出口的噴流溫度和速度高,噴流中CO2和H2O 的含量高,發(fā)動機噴管尺寸大且噴管的間距小于400 mm。這些特點導(dǎo)致噴流之間形成的相互干擾作用更強,噴流的輻射加熱和對流加熱的共同作用使得火箭底部熱環(huán)境更加惡劣。因此,開展對新一代四機并聯(lián)大推力液氧煤油發(fā)動機噴流特性及噴流干擾對底部熱環(huán)境影響的研究是必要的。

    多噴管火箭發(fā)動機是由單噴管火箭發(fā)動機并聯(lián)組成的,研究單噴管噴流流場可為火箭底部熱環(huán)境的研究打下堅實的基礎(chǔ)。文獻[1?4]針對單噴管火箭發(fā)動機,主要研究了不考慮噴流的二次燃燒時,不同飛行馬赫數(shù)、飛行高度對尾流場的膨脹壓縮波系個數(shù)、膨脹半徑、核心區(qū)長度的影響。文獻[5?8]研究了尾焰流場的復(fù)燃下不同飛行條件對單噴管尾焰軸向和徑向溫度及組分的影響。文獻[9]采用LU 隱式算法數(shù)值模擬了有無噴流干擾對軸對稱超聲速底部流場結(jié)構(gòu)的影響。文獻[10]采用RANS/LES 混合方法研究了超聲速底部流場的流動機理,得出IDDES方法預(yù)測的下游渦結(jié)構(gòu)更加清晰。

    單噴管火箭發(fā)動機的研究結(jié)果不能簡單地疊加之后運用于多噴管火箭發(fā)動機,因此文獻[11?12]通過試驗測量研究了單機和雙機氫氧發(fā)動機噴流的激波總壓力場,結(jié)果表明在相同的軸向位置,雙機干擾噴流的最高壓力是單機噴流的5 倍,噴流的干擾范圍隨軸向距離的增大而增大。文獻[13]數(shù)值仿真研究了雙機并聯(lián)氫氧火箭發(fā)動機尾焰流場,發(fā)現(xiàn)兩臺發(fā)動機尾焰流場的相互擠壓,形成撞擊激波,使內(nèi)區(qū)燃氣溫度、壓力升高,速度降低。文獻[14?15]對四噴管火箭發(fā)動機進行三維數(shù)值模擬,研究了不同飛行高度對尾流流場和火箭底部輻射熱流密度的影響,結(jié)果表明:隨飛行高度的增加,噴流擴張角變大,四噴管噴流之間形成干擾激波,并產(chǎn)生回流,在飛行高度為20 km左右,底板的中心位置出現(xiàn)最大熱流密度。

    火箭底部防熱問題是火箭研制階段應(yīng)重點關(guān)注的技術(shù)問題。研制階段對火箭底部熱環(huán)境評估過低會導(dǎo)致防熱措施不足,降低火箭的可靠性;評估過高則使防熱結(jié)構(gòu)質(zhì)量增加,降低運載能力。因此研究火箭發(fā)動機噴流對箭體底部的加熱和箭體底部流場結(jié)構(gòu)有利于了解火箭底部熱環(huán)境。本文基于標準κ-ε湍流模型與離散坐標輻射(discrete ordinates,DO)模型,采用FLUENT 軟件對四機并聯(lián)火箭發(fā)動機整體建模計算,綜合比較分析了不同飛行高度下噴流流場及噴流干擾形成的回流對火箭底部熱流密度的影響,可為合理評估火箭底部熱環(huán)境,進而設(shè)計出合理的熱防護結(jié)構(gòu)提供理論依據(jù)。

    1 物理模型及數(shù)值計算方法

    1.1 物理模型及網(wǎng)格化分

    火箭一級采用四個液氧煤油的液體發(fā)動機,火箭幾何模型如圖1 所示。模型網(wǎng)格采用ICEM 軟件完成,采用結(jié)構(gòu)化網(wǎng)格,針對不同的計算區(qū)域采用不同大小的網(wǎng)格,對發(fā)動機推力室及噴管出口附近進行網(wǎng)格加密處理,如圖2所示,共生成1 500萬六面體網(wǎng)格。發(fā)動機的參數(shù)如表1所示。

    表1 發(fā)動機主要參數(shù)Tab.1 Parameters of rocket Engine

    圖1 運載火箭幾何模型Fig.1 Geometry model of launch vehicle

    圖2 模型網(wǎng)格和尾部網(wǎng)格劃分Fig.2 Model meshes and afterbody meshes

    1.2 控制方程

    1.2.1多組分N-S方程

    基于液體運載火箭多組分燃氣滿足連續(xù)介質(zhì)假設(shè)和理想氣體狀態(tài)方程的特點,建立三維多組分可壓縮Navier-Stokes 方程,包括質(zhì)量、動量與能量輸運的控制方程,形式如下。

    多組分輸運方程為

    式中:t為時間;ρ為密度;ν為速度向量;l為多組分中的一組分;Yl為組分l的質(zhì)量分數(shù);Rl為組分l在化學(xué)反應(yīng)后的凈生成率;Sl為自定義源項的離散相引起的生成率。

    組分擴散通量Jl可表示為

    式中:T為環(huán)境溫度;Dl,m為組份l的質(zhì)量耗散系數(shù);DT,l為組份l的熱擴散系數(shù)。

    在直角坐標系中,單一組分l的可壓縮N-S 方程的守恒形式為

    式(3)~(6)中:U為流動變量;F、G、H為氣流通量矢量,F(xiàn)v、Gv、Hv為黏性通量矢量;K為熱傳導(dǎo)系數(shù);T為環(huán)境溫度;p、ρ、e、τ、μ分別為壓力、密度、比動能、應(yīng)力、黏性系數(shù);u、v、w分別為速度在x、y、z方向上的分量。

    注:式(6)中,ux為u在x方向上的偏導(dǎo)。

    1.2.2湍流模型

    標準κ-ε湍流模型,近壁面為標準壁面函數(shù),此湍流模型本身具有較好的穩(wěn)定性、經(jīng)濟性和比較高的計算精度。

    湍流輸運方程為

    式(7)~(9)中:k為湍流動能;ε為湍流耗散率;Gk為由平均速度梯度引起的湍動能k的產(chǎn)生項;Gb為由浮力引起的湍動能k的產(chǎn)生項;YM為可壓湍流中的脈動擴張項;σk、σε為與湍動能k和耗散率ε二者相對應(yīng)的Prandtl數(shù);Sk、Sε為用戶定義的源項;C1ε、C2ε、C3ε、Cμ是經(jīng)驗系數(shù),通常為常數(shù)。

    1.2.3輻射傳輸模型

    離散坐標輻射(DO)模型,DO 模型使用范圍最大的輻射模型,可以計算所有光學(xué)厚度的輻射問題。輻射模型中計算氣體的吸收系數(shù)時采用了FLUENT 軟件提供的灰氣體加權(quán)平均模型(WSGGM)。

    輻射方程為

    式中:r為位置向量;s為方向向量;a為氣體吸收系數(shù);σs為氣體散射系數(shù);n為氣體折射率;σ為黑體輻射常數(shù);I(r,s)為位置r處的微元體在s方向上的輻射出射度;s′為散射方向向量;Φ(s,s′)為散射相函數(shù);Ω′為空間立體角。

    1.2 邊界條件

    噴管喉部設(shè)置為壓力進口邊界,采用發(fā)動機穩(wěn)態(tài)工作時的總壓和總溫;箭體頭部和周圍外部進口邊界為壓力遠場邊界,計算考慮外界環(huán)境壓力和溫度的變化對噴流的影響;箭體尾部邊界為壓力出口邊界,邊界設(shè)置與飛行高度有關(guān);箭體、噴管外壁和箭體底部為300 K定溫邊界。

    根據(jù)海拔高度不同設(shè)置4 個工況,對底部熱環(huán)境進行仿真分析,工況設(shè)置如表2所示。

    表2 自由來流參數(shù)Tab.2 Parameters of free flow

    2 計算結(jié)果分析

    2.1 不同飛行高度馬赫數(shù)分布分析

    圖3 是Z=0 截面馬赫數(shù)等值線圖,飛行高度為7 km 時,在箭體頂端沒有形成激波,噴流之間幾乎無交互作用。隨著飛行時間的增加,火箭進入超音速飛行,箭體頭部形成弓形激波,噴流的擴張角變大,自由來流與噴流相互作用形成羽流激波、桶形激波和自由剪切層即膨脹邊界。在噴流內(nèi)部即噴流相互作用的區(qū)域,噴流之間形成干擾激波,從圖4 可以看出,速度會出現(xiàn)負值,這表明底部出現(xiàn)回流,回流的軸向速度隨飛行高度的增加而升高。速度為零的位置是噴流之間發(fā)生交匯并產(chǎn)生返流的位置,不同飛行高度返流點位置如表3所示。當飛行高度高于13 km,從X軸向速度線圖可以看出,開始回流點的軸向位置隨飛行高度的增加而減小。

    表3 不同飛行高度下馬赫數(shù)為零的軸向位置Tab.3 The axial position where Mach number is zero at different flight altitudes

    圖3 不同飛行高度馬赫數(shù)等值線圖Fig.3 Isograms of Mach number at different flight altitudes

    圖4 不同飛行高度X軸向馬赫數(shù)分布Fig.4 Axial distribution curve of the Mach number at different flight altitudes

    2.2 高低空底部流場對比分析

    圖5(a)、(b)是飛行高度為7 km 時底部參數(shù)分布圖。從圖可以看出,在飛行高度為7 km 時,火箭飛行馬赫數(shù)為0.72,發(fā)動機噴流主要集中在噴管中心區(qū)域,噴流膨脹角較小,噴流之間因膨脹產(chǎn)生相互作用而形成一個壓力較高的區(qū)域,噴流的引射作用將空氣吸入底部區(qū)域,被吸入底部的自由流空氣向下游流動時不能通過噴流相互作用區(qū),在發(fā)動機出口附近形成腎形渦對,增強了與噴流摻混。噴流相互作用區(qū)與底部的壓差較小,返流以亞音速沖擊底板中心,并在底板附近形成沿噴管間隙流動的壁面射流,此時壁面射流的馬赫數(shù)小于0.1。圖5(c)、(d)給出了飛行高度為21 km 時底部參數(shù)分布圖。在飛行高度為21 km 時,飛行馬赫數(shù)為1.89,環(huán)境壓力降低,噴流膨脹角增大,四噴管噴流的相互作用在四噴管出口附近形成一個高壓噴流相互作用區(qū)。噴管內(nèi)壁面的噴流邊界層隨噴流的膨脹波流動到噴流相互作用區(qū),但低動能的噴流邊界層[16]與被吸入底部的自由流空氣不能通過高壓的噴流相互作用區(qū),因此出現(xiàn)返流現(xiàn)象。噴流相互作用區(qū)與底部之間形成很大的壓差,返流逐漸加速到超音速,形成返流激波,最后在距離底板中心0.2 m 左右形成一個正激波,超音速返流穿過正激波后,導(dǎo)致返流降為亞音速,壓力和溫度急劇升高。相對于低空,高空的壁面射流的馬赫數(shù)增大到0.6 左右。從圖6(a)、(b)可以看出,低空時返流的溫度較低,高空時返流溫度較高,致使底板的熱環(huán)境更加惡劣。

    圖5 底部截面馬赫數(shù)流線圖Fig.5 Mach number flow field at the bottom of rocket

    圖6 溫度分布云圖Fig.6 Temperature contour

    從圖7(a)、(b)可知低空時返流主要為空氣,高空時返流主要集中在軸線附近的區(qū)域,影響區(qū)域變大且高溫噴流氣體摩爾分數(shù)占比增大。從圖8 可以看出,在X=2.042 m 時X軸向速度為零,此位置是開始返流的位置,在X=0~2.042 m 速度會出現(xiàn)負值,這表明底部出現(xiàn)返流。從開始返流到底板,返流空氣摩爾分數(shù)從90%左右逐漸升高到99%左右,這說明返流氣體主要是被吸入底部的空氣,只有少量被卷吸的噴流氣體。返流溫度從900 K左右逐漸下降到357 K左右,如圖9 所示。飛行高度為21 km 時開始返流的位置為X=2.105 m。從開始返流到底板,返流空氣摩爾分數(shù)從19%左右逐漸升高到42%左右,如圖10所示,空氣摩爾分數(shù)逐漸升高是因為返流過程中不斷地與被吸入底部區(qū)域的空氣摻混,但返流氣體組分中噴流氣體占50%以上。返流溫度從3 200 K 左右波動下降到2 780 K左右,如圖11所示。

    圖7 空氣摩爾分數(shù)流線圖Fig.7 Mole fraction of air

    圖8 7 km X軸向馬赫數(shù)和空氣摩爾分數(shù)分布曲線Fig.8 Axial distribution curve of the Mach number and mole fraction of air at 7 km

    圖10 21 km X軸向馬赫數(shù)和空氣摩爾分數(shù)分布曲線Fig.10 Axial distribution curve of the Mach number and mole fraction of air at 21 km

    圖11 21 km X軸向壓力和溫度分布曲線Fig.11 Axial distribution curve of the pressure and temperature at 21 km

    2.3 不同飛行高度底部熱流分布分析

    圖12 為底部熱流密度分布云圖。從圖中可以看出,防熱底板的最大總熱流密度出現(xiàn)在中心位置,高熱流密度區(qū)主要分布在中心區(qū)域和靠近邊緣的各噴管之間狹長區(qū)域。這一是因為這兩區(qū)域受輻射的面積較大;二是因為高溫反流直接沖擊防熱底板后沿著各噴管之間狹長區(qū)域流出。噴管周圍區(qū)域的熱流密度都較小,主要是噴流作用相當于引射器,將自由流空氣吸入底部區(qū)域,起到冷卻和減弱熱輻射的作用。隨著飛行高度的升高,防熱底板熱流密度呈現(xiàn)先升高后降低的趨勢,如圖13 所示,最大總熱流密度出現(xiàn)在飛行高度21 km 處,總熱流密度達到748 kW/m2。飛行高度為30 km 到60 km 期間,噴流返流的速度和溫度增大幅度較小,但噴流返流的密度下降,導(dǎo)致總熱流密度呈下降趨勢。

    圖12 熱流密度分布云圖Fig.12 Heat flux contour

    圖13 防熱底板最大熱流密度隨時間變化圖Fig.13 Max.heat flux vs time on the bottom of the rocket

    發(fā)動機的總熱流密度隨高度的增加呈現(xiàn)先升高后下降的趨勢,如圖14 所示,在飛行高度為7 km 時總熱流密度達到最大,為297 kW/m2。發(fā)動機大熱流密度區(qū)主要分布在發(fā)動機出口靠近箭體軸線的位置,這主要是因為四個發(fā)動機的高溫交匯區(qū)距離發(fā)動機出口較近,輻射作用對發(fā)動機的加熱占主要地位。尾段的熱流密度隨飛行高度的增大,從負熱流密度向正熱流密度變化,總熱流密度是先升高后降低,如圖15 所示。隨著飛行高度的繼續(xù)增大,尾段主要受輻射加熱的影響,尾段總熱流密度增大,總熱流密度在飛行高度42 km時達到最大,為18 kW/m2。

    圖14 發(fā)動機最大熱流密度隨時間變化圖Fig.14 Max.heat flux vs time of the rocket engine

    圖15 尾段最大熱流密度隨時間變化圖Fig.15 Max.heat flux vs time of the afterbody

    3 結(jié)束語

    本文通過對四機并聯(lián)火箭發(fā)動機噴流流場進行數(shù)值仿真,分析了不同飛行高度下噴流干擾形成的回流特性和箭體底部熱流分布。分析結(jié)果表明,四臺發(fā)動機噴流在不同飛行高度下,在發(fā)動機出口附近產(chǎn)生一個交互回流區(qū)。通過分析高低空底部流場發(fā)現(xiàn):隨飛行高度的增加,返流向底板流動速度從亞音速加速到超音速,返流溫度也逐漸升高,防熱底板熱流密度呈現(xiàn)先升高后降低的趨勢,最大熱流密度點出現(xiàn)在飛行高度21 km 處。本文對流場的結(jié)構(gòu)、返流的溫度和組分等方面的研究,可以較為清晰地了解不同飛行高度的熱流密度變化原因,對箭體底部防熱區(qū)域細化,進而減輕防熱結(jié)構(gòu)重量具有較強的指導(dǎo)意義。

    猜你喜歡
    箭體飛行高度噴流
    航天箭體用薄壁結(jié)構(gòu)機械加工噪聲分析及防治
    飛行參數(shù)對六旋翼植保無人機霧滴在荔枝樹冠層沉積分布的影響
    “慧眼”發(fā)現(xiàn)迄今距離黑洞最近的高速噴流
    簡析二次雷達高度信息與飛機實際高度的關(guān)系
    不簡單的小旋翼
    航空模型(2017年3期)2017-07-28 20:56:40
    噴流干擾氣動熱數(shù)值模擬的若干影響因素
    FAA頒新政使小型無人機飛行高度翻倍
    航空模型(2016年5期)2016-07-25 08:59:26
    骨骼清奇 合理最美——淺析長征五號箭體結(jié)構(gòu)
    太空探索(2016年12期)2016-07-18 11:13:43
    耀變體噴流高能電子譜的形成機制
    “平流層”是個啥——話說飛行高度
    午夜免费观看性视频| 18在线观看网站| 又大又黄又爽视频免费| 丰满迷人的少妇在线观看| 色婷婷久久久亚洲欧美| 天天躁狠狠躁夜夜躁狠狠躁| 狠狠婷婷综合久久久久久88av| 五月伊人婷婷丁香| av在线播放精品| 久久精品久久久久久久性| 两个人免费观看高清视频| av.在线天堂| 精品国产一区二区久久| 久久韩国三级中文字幕| 久久精品国产亚洲av天美| 老汉色av国产亚洲站长工具| 91aial.com中文字幕在线观看| 久久午夜综合久久蜜桃| 免费看av在线观看网站| 在线天堂中文资源库| 18在线观看网站| 欧美少妇被猛烈插入视频| 夜夜骑夜夜射夜夜干| 亚洲国产欧美日韩在线播放| 女人久久www免费人成看片| 国产 一区精品| 亚洲av电影在线观看一区二区三区| 亚洲国产色片| 国产精品国产av在线观看| 在线看a的网站| 十分钟在线观看高清视频www| 嫩草影院入口| 伊人久久大香线蕉亚洲五| 国产一区二区三区综合在线观看| 日韩精品有码人妻一区| 日韩欧美精品免费久久| 天天影视国产精品| 人成视频在线观看免费观看| 欧美少妇被猛烈插入视频| 97人妻天天添夜夜摸| 精品一区二区三卡| 视频在线观看一区二区三区| 在线观看美女被高潮喷水网站| 麻豆乱淫一区二区| 亚洲在久久综合| 久久99精品国语久久久| 亚洲 欧美一区二区三区| 啦啦啦啦在线视频资源| 大香蕉久久网| 少妇人妻久久综合中文| 久久97久久精品| 精品少妇一区二区三区视频日本电影 | 欧美 亚洲 国产 日韩一| 亚洲少妇的诱惑av| 纵有疾风起免费观看全集完整版| 亚洲欧洲精品一区二区精品久久久 | 国产精品久久久久久av不卡| 国产日韩欧美亚洲二区| 在线天堂最新版资源| 日日撸夜夜添| 亚洲一区中文字幕在线| 国产成人午夜福利电影在线观看| 免费人妻精品一区二区三区视频| 啦啦啦中文免费视频观看日本| 国产成人欧美| 久久女婷五月综合色啪小说| 一本色道久久久久久精品综合| 一本久久精品| 午夜免费观看性视频| 日韩一本色道免费dvd| 大话2 男鬼变身卡| 菩萨蛮人人尽说江南好唐韦庄| 国产综合精华液| 人妻系列 视频| 国产成人午夜福利电影在线观看| 青草久久国产| 国产成人a∨麻豆精品| 999精品在线视频| 国产精品国产av在线观看| 边亲边吃奶的免费视频| 天堂8中文在线网| 亚洲成人一二三区av| 日本黄色日本黄色录像| 人人妻人人澡人人看| 国产精品久久久久久av不卡| 女人高潮潮喷娇喘18禁视频| 边亲边吃奶的免费视频| 两个人看的免费小视频| 亚洲美女搞黄在线观看| 欧美黄色片欧美黄色片| 亚洲熟女精品中文字幕| 大片免费播放器 马上看| 国产乱人偷精品视频| 国产成人精品在线电影| 亚洲av福利一区| 9热在线视频观看99| 少妇的逼水好多| 亚洲精品久久久久久婷婷小说| 国产日韩一区二区三区精品不卡| 久久久久久人妻| 国产高清不卡午夜福利| √禁漫天堂资源中文www| 一级黄片播放器| 欧美黄色片欧美黄色片| 亚洲av男天堂| 99久国产av精品国产电影| 精品一品国产午夜福利视频| 午夜av观看不卡| 国精品久久久久久国模美| 久久精品国产a三级三级三级| 叶爱在线成人免费视频播放| 国产精品欧美亚洲77777| 人妻系列 视频| 日韩,欧美,国产一区二区三区| 久久精品国产自在天天线| 99久久综合免费| 国产亚洲精品第一综合不卡| 国产精品 国内视频| 欧美精品av麻豆av| 丝瓜视频免费看黄片| 叶爱在线成人免费视频播放| 国产高清国产精品国产三级| 老熟女久久久| 亚洲欧美一区二区三区黑人 | 亚洲精品美女久久久久99蜜臀 | 国产精品偷伦视频观看了| 日本av手机在线免费观看| 成年av动漫网址| 国产精品av久久久久免费| av又黄又爽大尺度在线免费看| 9热在线视频观看99| 1024香蕉在线观看| 日韩制服骚丝袜av| 少妇 在线观看| 在线精品无人区一区二区三| 亚洲精品乱久久久久久| 性色av一级| av有码第一页| 国产精品三级大全| 老汉色∧v一级毛片| 国产白丝娇喘喷水9色精品| 最近最新中文字幕免费大全7| av女优亚洲男人天堂| 国产伦理片在线播放av一区| 日日撸夜夜添| 热re99久久精品国产66热6| 性高湖久久久久久久久免费观看| 国产精品三级大全| 99久久人妻综合| 制服丝袜香蕉在线| 伦理电影免费视频| 精品人妻一区二区三区麻豆| 亚洲一区二区三区欧美精品| 日日爽夜夜爽网站| 一区二区日韩欧美中文字幕| 亚洲美女视频黄频| 日韩一区二区三区影片| av.在线天堂| 伦理电影免费视频| 在线观看人妻少妇| 街头女战士在线观看网站| 丝袜美足系列| kizo精华| 成人国语在线视频| 国产精品久久久久久精品古装| 亚洲成色77777| 高清av免费在线| 三级国产精品片| 性高湖久久久久久久久免费观看| 激情五月婷婷亚洲| 国产日韩欧美亚洲二区| 国产 一区精品| 亚洲综合精品二区| 精品少妇黑人巨大在线播放| av又黄又爽大尺度在线免费看| 国产乱来视频区| 精品酒店卫生间| av又黄又爽大尺度在线免费看| 麻豆精品久久久久久蜜桃| 亚洲中文av在线| 自拍欧美九色日韩亚洲蝌蚪91| 麻豆精品久久久久久蜜桃| 国产免费现黄频在线看| 欧美国产精品一级二级三级| 国产在线一区二区三区精| 欧美精品人与动牲交sv欧美| 最近中文字幕高清免费大全6| 日本91视频免费播放| 日韩av不卡免费在线播放| 亚洲 欧美一区二区三区| 少妇精品久久久久久久| 最新的欧美精品一区二区| 午夜福利视频精品| 18禁国产床啪视频网站| 亚洲精品乱久久久久久| 三上悠亚av全集在线观看| 免费观看在线日韩| 久久99热这里只频精品6学生| 午夜福利视频在线观看免费| 一区在线观看完整版| 国语对白做爰xxxⅹ性视频网站| 国产av一区二区精品久久| 日韩伦理黄色片| 亚洲精品日本国产第一区| 欧美黄色片欧美黄色片| 2018国产大陆天天弄谢| 丝瓜视频免费看黄片| a级毛片在线看网站| 美女福利国产在线| 免费在线观看完整版高清| 免费播放大片免费观看视频在线观看| 久久久久国产精品人妻一区二区| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久久精品精品| 国产精品久久久av美女十八| 久久精品久久久久久噜噜老黄| 男人爽女人下面视频在线观看| 久久鲁丝午夜福利片| 日韩中文字幕视频在线看片| 久久久久久人妻| 90打野战视频偷拍视频| 又粗又硬又长又爽又黄的视频| 超碰97精品在线观看| 成年美女黄网站色视频大全免费| 久久 成人 亚洲| 午夜免费观看性视频| 亚洲欧美成人综合另类久久久| 国产 精品1| 欧美人与善性xxx| 各种免费的搞黄视频| 免费看av在线观看网站| 亚洲一区二区三区欧美精品| 亚洲综合精品二区| 性高湖久久久久久久久免费观看| 黄片播放在线免费| 久久久精品免费免费高清| 国产欧美日韩综合在线一区二区| 黄网站色视频无遮挡免费观看| 国产av码专区亚洲av| 亚洲国产精品一区三区| 成年动漫av网址| 久久久久精品久久久久真实原创| 黄网站色视频无遮挡免费观看| 免费观看在线日韩| videossex国产| 亚洲第一青青草原| 丰满乱子伦码专区| 精品国产一区二区三区四区第35| 性少妇av在线| 亚洲av综合色区一区| 大话2 男鬼变身卡| 色播在线永久视频| 欧美激情极品国产一区二区三区| 青春草亚洲视频在线观看| 91成人精品电影| 黄片无遮挡物在线观看| 日韩 亚洲 欧美在线| 久久久久久久久久人人人人人人| 最近中文字幕高清免费大全6| 日韩欧美一区视频在线观看| 亚洲国产欧美在线一区| 精品亚洲成a人片在线观看| 街头女战士在线观看网站| 亚洲精品日本国产第一区| 久久久久久久久久久免费av| 免费观看性生交大片5| 80岁老熟妇乱子伦牲交| 欧美精品国产亚洲| 欧美精品高潮呻吟av久久| 免费在线观看黄色视频的| 少妇熟女欧美另类| 一级片免费观看大全| 精品国产国语对白av| 国产成人一区二区在线| 午夜久久久在线观看| 乱人伦中国视频| 只有这里有精品99| 久久精品夜色国产| 天天躁夜夜躁狠狠久久av| 咕卡用的链子| 国产欧美亚洲国产| 综合色丁香网| 亚洲欧美成人精品一区二区| 久久久国产一区二区| 国产乱人偷精品视频| 哪个播放器可以免费观看大片| 侵犯人妻中文字幕一二三四区| 熟女电影av网| 天天躁夜夜躁狠狠久久av| 日本wwww免费看| 最近中文字幕高清免费大全6| 9热在线视频观看99| 亚洲精品美女久久久久99蜜臀 | 在线观看免费视频网站a站| 日韩中字成人| 亚洲欧美精品综合一区二区三区 | 啦啦啦在线观看免费高清www| 又黄又粗又硬又大视频| 伊人久久大香线蕉亚洲五| 免费女性裸体啪啪无遮挡网站| 一二三四中文在线观看免费高清| 国产成人一区二区在线| 亚洲情色 制服丝袜| 日韩精品免费视频一区二区三区| 亚洲综合色网址| 丁香六月天网| 国产麻豆69| 久久久久久久大尺度免费视频| 欧美日韩成人在线一区二区| 国产精品99久久99久久久不卡 | 精品午夜福利在线看| 亚洲精品一区蜜桃| 国产欧美日韩综合在线一区二区| 三上悠亚av全集在线观看| 两性夫妻黄色片| 最近2019中文字幕mv第一页| 日韩人妻精品一区2区三区| 免费日韩欧美在线观看| 免费女性裸体啪啪无遮挡网站| 午夜日韩欧美国产| 老汉色av国产亚洲站长工具| 青青草视频在线视频观看| 免费看不卡的av| 国产激情久久老熟女| 免费少妇av软件| 成人国语在线视频| 久久久久久久亚洲中文字幕| 成年女人毛片免费观看观看9 | 一本久久精品| 日本免费在线观看一区| 人妻人人澡人人爽人人| 精品国产乱码久久久久久男人| 男男h啪啪无遮挡| 国产欧美亚洲国产| 在线观看免费高清a一片| 久久精品国产亚洲av涩爱| 国产亚洲最大av| 久久久精品免费免费高清| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av福利一区| tube8黄色片| av免费在线看不卡| 一个人免费看片子| 日韩中文字幕视频在线看片| 国产xxxxx性猛交| 成人国语在线视频| 久久国产精品男人的天堂亚洲| 成人国产av品久久久| 免费在线观看完整版高清| 国产视频首页在线观看| 人妻 亚洲 视频| 有码 亚洲区| 国产又爽黄色视频| 亚洲精品国产一区二区精华液| 国产又爽黄色视频| 久久久国产一区二区| 巨乳人妻的诱惑在线观看| 久久久久视频综合| 精品国产一区二区三区四区第35| www.精华液| 日韩一区二区视频免费看| 欧美日韩亚洲国产一区二区在线观看 | 欧美精品av麻豆av| 亚洲国产精品999| 精品第一国产精品| 亚洲国产日韩一区二区| 久久 成人 亚洲| 伦理电影大哥的女人| 日日爽夜夜爽网站| 免费高清在线观看日韩| 99久久精品国产国产毛片| 香蕉丝袜av| 我要看黄色一级片免费的| 香蕉丝袜av| 2021少妇久久久久久久久久久| 午夜福利在线观看免费完整高清在| 春色校园在线视频观看| 久久99一区二区三区| www日本在线高清视频| 国产成人精品无人区| 久久国内精品自在自线图片| 边亲边吃奶的免费视频| 秋霞在线观看毛片| 亚洲色图 男人天堂 中文字幕| 热re99久久精品国产66热6| 国产精品无大码| 婷婷色av中文字幕| 免费在线观看黄色视频的| 亚洲精品国产av蜜桃| 欧美日韩一级在线毛片| 久久精品国产亚洲av天美| 人人妻人人澡人人看| 亚洲av在线观看美女高潮| 亚洲国产av影院在线观看| 亚洲av在线观看美女高潮| 七月丁香在线播放| 国产又爽黄色视频| 亚洲欧洲日产国产| 久久国产精品男人的天堂亚洲| 2018国产大陆天天弄谢| 蜜桃在线观看..| 久久久久人妻精品一区果冻| 国产探花极品一区二区| 亚洲国产欧美网| 久久久久久久久久人人人人人人| 男男h啪啪无遮挡| 国产av精品麻豆| 国产精品久久久久成人av| 成年美女黄网站色视频大全免费| 国产精品99久久99久久久不卡 | 国产精品女同一区二区软件| 国产成人精品福利久久| 人妻人人澡人人爽人人| 人成视频在线观看免费观看| 黄网站色视频无遮挡免费观看| 亚洲熟女精品中文字幕| 亚洲精品一区蜜桃| 午夜激情久久久久久久| 免费看不卡的av| 色视频在线一区二区三区| 美女高潮到喷水免费观看| 王馨瑶露胸无遮挡在线观看| 最近最新中文字幕大全免费视频 | 欧美 亚洲 国产 日韩一| 亚洲精品aⅴ在线观看| 一级a爱视频在线免费观看| 亚洲人成电影观看| 久久久久久久大尺度免费视频| 久久这里只有精品19| 少妇精品久久久久久久| www.熟女人妻精品国产| 精品少妇黑人巨大在线播放| 久久精品国产综合久久久| 日本-黄色视频高清免费观看| 久久久久精品性色| 青草久久国产| 久久精品国产a三级三级三级| 十分钟在线观看高清视频www| 美女午夜性视频免费| 欧美少妇被猛烈插入视频| 伦精品一区二区三区| 免费观看无遮挡的男女| 免费观看在线日韩| 韩国av在线不卡| 午夜av观看不卡| 国产深夜福利视频在线观看| 啦啦啦在线观看免费高清www| 久久毛片免费看一区二区三区| 性色av一级| 这个男人来自地球电影免费观看 | 亚洲三区欧美一区| 肉色欧美久久久久久久蜜桃| 久久精品国产综合久久久| 国产免费福利视频在线观看| 国产精品国产三级国产专区5o| 亚洲一区二区三区欧美精品| 男人爽女人下面视频在线观看| 久久精品国产亚洲av高清一级| 97人妻天天添夜夜摸| 91aial.com中文字幕在线观看| 丝瓜视频免费看黄片| 国产精品亚洲av一区麻豆 | 亚洲av中文av极速乱| 精品国产超薄肉色丝袜足j| 99热国产这里只有精品6| 国产野战对白在线观看| 欧美少妇被猛烈插入视频| 99国产精品免费福利视频| 亚洲婷婷狠狠爱综合网| 美女大奶头黄色视频| 啦啦啦啦在线视频资源| 国产成人免费观看mmmm| 久久毛片免费看一区二区三区| av国产精品久久久久影院| 国产熟女欧美一区二区| 亚洲欧美精品综合一区二区三区 | 91在线精品国自产拍蜜月| 一级爰片在线观看| av福利片在线| 欧美亚洲日本最大视频资源| 日本午夜av视频| 男人舔女人的私密视频| 一区二区av电影网| 日本91视频免费播放| 欧美精品国产亚洲| 久久精品熟女亚洲av麻豆精品| 国产精品亚洲av一区麻豆 | 日韩视频在线欧美| 夫妻午夜视频| 免费看不卡的av| 一级毛片电影观看| 欧美97在线视频| 欧美少妇被猛烈插入视频| 亚洲,欧美精品.| 亚洲精品久久午夜乱码| 精品酒店卫生间| 一区二区三区乱码不卡18| 永久网站在线| 久久久久网色| 久久精品人人爽人人爽视色| 亚洲国产最新在线播放| 永久网站在线| 国产精品久久久久成人av| av网站免费在线观看视频| 侵犯人妻中文字幕一二三四区| 久久精品亚洲av国产电影网| 久久久精品区二区三区| 亚洲精品日韩在线中文字幕| 男人爽女人下面视频在线观看| 欧美精品人与动牲交sv欧美| 麻豆乱淫一区二区| 如日韩欧美国产精品一区二区三区| 久久久久久免费高清国产稀缺| 免费在线观看黄色视频的| 91成人精品电影| 日韩中文字幕视频在线看片| av在线老鸭窝| 亚洲综合色网址| 国产精品不卡视频一区二区| 久久精品国产自在天天线| 三上悠亚av全集在线观看| 亚洲国产av影院在线观看| 久久久久久久亚洲中文字幕| 国产精品 欧美亚洲| 深夜精品福利| 日日摸夜夜添夜夜爱| 国产在线免费精品| 久久韩国三级中文字幕| 你懂的网址亚洲精品在线观看| 国产精品嫩草影院av在线观看| 国产福利在线免费观看视频| 久久人人爽人人片av| 美女主播在线视频| 欧美成人精品欧美一级黄| 国产 一区精品| 国产一级毛片在线| 性少妇av在线| 亚洲精品久久成人aⅴ小说| 搡老乐熟女国产| 欧美激情极品国产一区二区三区| 天天躁日日躁夜夜躁夜夜| 精品福利永久在线观看| 一级a爱视频在线免费观看| 欧美亚洲 丝袜 人妻 在线| 一二三四在线观看免费中文在| 亚洲色图 男人天堂 中文字幕| 久热久热在线精品观看| 国产精品偷伦视频观看了| 欧美亚洲日本最大视频资源| 国产成人91sexporn| 精品亚洲乱码少妇综合久久| 免费不卡的大黄色大毛片视频在线观看| av视频免费观看在线观看| 免费观看性生交大片5| 99国产综合亚洲精品| 777久久人妻少妇嫩草av网站| 国产精品一国产av| www.自偷自拍.com| 色婷婷av一区二区三区视频| 久久精品国产亚洲av天美| 欧美bdsm另类| 成年av动漫网址| 亚洲综合色惰| 国产又爽黄色视频| 男女国产视频网站| 亚洲成av片中文字幕在线观看 | 纯流量卡能插随身wifi吗| 成人国产麻豆网| 在线观看人妻少妇| www.精华液| 亚洲男人天堂网一区| 国产成人免费无遮挡视频| 亚洲av电影在线进入| 亚洲精品国产av成人精品| 亚洲一码二码三码区别大吗| 亚洲欧美一区二区三区国产| 如日韩欧美国产精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 一本色道久久久久久精品综合| 午夜福利乱码中文字幕| 丰满迷人的少妇在线观看| 国产一区二区激情短视频 | 国产福利在线免费观看视频| 18禁观看日本| www.熟女人妻精品国产| 中文字幕人妻丝袜制服| 欧美亚洲日本最大视频资源| 欧美中文综合在线视频| 亚洲av电影在线进入| 在线观看www视频免费| 亚洲av福利一区| 日韩熟女老妇一区二区性免费视频| 韩国av在线不卡| xxx大片免费视频| 女人久久www免费人成看片| 一边摸一边做爽爽视频免费| av又黄又爽大尺度在线免费看| 美女高潮到喷水免费观看| 99热全是精品| 成年女人毛片免费观看观看9 | 乱人伦中国视频| 天堂8中文在线网| 久久久久久久久久久久大奶| 涩涩av久久男人的天堂| av国产精品久久久久影院| 亚洲精品视频女| 午夜影院在线不卡| 精品少妇内射三级| 免费人妻精品一区二区三区视频| 美女午夜性视频免费| videossex国产| av免费观看日本| 亚洲视频免费观看视频| 老司机亚洲免费影院| 菩萨蛮人人尽说江南好唐韦庄|