• <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
    耀變體噴流高能電子譜的形成機制
    “平流層”是個啥——話說飛行高度
    国产 精品1| 日韩不卡一区二区三区视频在线| 精品国产露脸久久av麻豆| 亚洲精品美女久久av网站| 一区二区日韩欧美中文字幕| 成年av动漫网址| 777久久人妻少妇嫩草av网站| 免费在线观看视频国产中文字幕亚洲 | 18禁动态无遮挡网站| 国产高清国产精品国产三级| 亚洲综合色惰| 国产精品无大码| 国产av码专区亚洲av| 日韩在线高清观看一区二区三区| 久热久热在线精品观看| 久久午夜福利片| 精品第一国产精品| 国产高清国产精品国产三级| 汤姆久久久久久久影院中文字幕| 大陆偷拍与自拍| 亚洲在久久综合| av视频免费观看在线观看| 久久精品熟女亚洲av麻豆精品| 国产欧美日韩一区二区三区在线| 精品少妇一区二区三区视频日本电影 | 国产精品久久久久久精品电影小说| www.精华液| 中文字幕人妻熟女乱码| 汤姆久久久久久久影院中文字幕| 亚洲激情五月婷婷啪啪| 老鸭窝网址在线观看| 纵有疾风起免费观看全集完整版| 精品国产露脸久久av麻豆| 国产av精品麻豆| 波野结衣二区三区在线| 十八禁网站网址无遮挡| 少妇被粗大的猛进出69影院| 精品一品国产午夜福利视频| 国产亚洲精品第一综合不卡| 久久久久久久国产电影| 麻豆乱淫一区二区| 国产成人aa在线观看| 男人爽女人下面视频在线观看| 9191精品国产免费久久| 亚洲精品av麻豆狂野| 国产综合精华液| 国产片内射在线| 亚洲精品国产色婷婷电影| www.精华液| 亚洲精品国产一区二区精华液| 成年女人毛片免费观看观看9 | av有码第一页| 深夜精品福利| 捣出白浆h1v1| 亚洲经典国产精华液单| 日本黄色日本黄色录像| 精品人妻熟女毛片av久久网站| 国产不卡av网站在线观看| 久久精品国产亚洲av涩爱| 国产欧美亚洲国产| 十八禁网站网址无遮挡| 人妻 亚洲 视频| 亚洲精品久久成人aⅴ小说| 亚洲第一av免费看| 高清在线视频一区二区三区| 在线亚洲精品国产二区图片欧美| 午夜av观看不卡| 男人操女人黄网站| 亚洲伊人久久精品综合| 亚洲综合色惰| 大陆偷拍与自拍| 久久毛片免费看一区二区三区| 少妇精品久久久久久久| 丝袜喷水一区| 在线观看人妻少妇| 久久精品国产a三级三级三级| 国产精品久久久久久久久免| 国产精品熟女久久久久浪| 国产1区2区3区精品| 中文字幕av电影在线播放| 久久午夜综合久久蜜桃| av一本久久久久| 午夜激情久久久久久久| 久久久久视频综合| 欧美日韩一级在线毛片| 免费高清在线观看视频在线观看| 人人妻人人澡人人爽人人夜夜| 不卡视频在线观看欧美| 久久久精品免费免费高清| 亚洲国产毛片av蜜桃av| 免费黄网站久久成人精品| 狂野欧美激情性bbbbbb| 国产精品久久久久成人av| 中文字幕精品免费在线观看视频| 秋霞在线观看毛片| 一级黄片播放器| 天美传媒精品一区二区| 美女xxoo啪啪120秒动态图| 一边亲一边摸免费视频| 2018国产大陆天天弄谢| 少妇人妻久久综合中文| 欧美bdsm另类| 人人妻人人澡人人看| av福利片在线| 最新的欧美精品一区二区| 久久99热这里只频精品6学生| 亚洲综合色惰| 一区在线观看完整版| 一级片免费观看大全| 在线看a的网站| 亚洲第一区二区三区不卡| 九草在线视频观看| 十八禁网站网址无遮挡| 亚洲精品一二三| 91精品三级在线观看| kizo精华| 18+在线观看网站| 在线观看免费视频网站a站| 美女午夜性视频免费| 精品视频人人做人人爽| 新久久久久国产一级毛片| 另类精品久久| 秋霞在线观看毛片| 麻豆精品久久久久久蜜桃| 国产乱来视频区| 久热久热在线精品观看| 波多野结衣av一区二区av| h视频一区二区三区| 国产精品不卡视频一区二区| 女性生殖器流出的白浆| 性高湖久久久久久久久免费观看| videosex国产| 亚洲欧美精品自产自拍| 亚洲第一区二区三区不卡| 亚洲av中文av极速乱| 免费高清在线观看视频在线观看| 极品人妻少妇av视频| 亚洲av福利一区| 青春草亚洲视频在线观看| 我的亚洲天堂| 久久午夜福利片| 欧美日韩av久久| 一级黄片播放器| 纯流量卡能插随身wifi吗| 精品99又大又爽又粗少妇毛片| 亚洲精品一区蜜桃| 一级爰片在线观看| 少妇的逼水好多| 不卡视频在线观看欧美| 岛国毛片在线播放| 亚洲精品自拍成人| 亚洲av.av天堂| 大片电影免费在线观看免费| 2018国产大陆天天弄谢| 久久久久久久久免费视频了| 久久精品国产综合久久久| 久久久久国产精品人妻一区二区| 欧美 日韩 精品 国产| 久久久国产欧美日韩av| 18+在线观看网站| 精品国产露脸久久av麻豆| 免费高清在线观看视频在线观看| 又大又黄又爽视频免费| 欧美日韩亚洲国产一区二区在线观看 | 亚洲国产毛片av蜜桃av| 亚洲精品一二三| 一级毛片黄色毛片免费观看视频| 超色免费av| 久久久精品国产亚洲av高清涩受| 激情五月婷婷亚洲| 好男人视频免费观看在线| 少妇 在线观看| 久久久久久人妻| 少妇的逼水好多| 亚洲经典国产精华液单| 18在线观看网站| 免费不卡的大黄色大毛片视频在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲国产最新在线播放| 女性被躁到高潮视频| 天天躁夜夜躁狠狠躁躁| 亚洲人成网站在线观看播放| 女性生殖器流出的白浆| 人成视频在线观看免费观看| 欧美日韩一级在线毛片| 久久久精品94久久精品| 国产成人a∨麻豆精品| 狂野欧美激情性bbbbbb| 乱人伦中国视频| 另类亚洲欧美激情| 国产成人精品在线电影| 亚洲激情五月婷婷啪啪| xxx大片免费视频| 久久国产亚洲av麻豆专区| 欧美精品亚洲一区二区| 黑人巨大精品欧美一区二区蜜桃| 国产午夜精品一二区理论片| 午夜福利在线观看免费完整高清在| 久久久亚洲精品成人影院| 在线观看免费日韩欧美大片| 菩萨蛮人人尽说江南好唐韦庄| 久久久精品94久久精品| 久久精品久久精品一区二区三区| 中文乱码字字幕精品一区二区三区| 久久久久国产精品人妻一区二区| 日日啪夜夜爽| 高清视频免费观看一区二区| 国产有黄有色有爽视频| freevideosex欧美| 午夜影院在线不卡| av线在线观看网站| 国产亚洲av片在线观看秒播厂| 欧美激情 高清一区二区三区| 99精国产麻豆久久婷婷| 精品福利永久在线观看| 青草久久国产| 中国三级夫妇交换| 激情五月婷婷亚洲| 久久影院123| 在现免费观看毛片| 蜜桃国产av成人99| 亚洲国产毛片av蜜桃av| 国产成人午夜福利电影在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩,欧美,国产一区二区三区| 你懂的网址亚洲精品在线观看| a级毛片在线看网站| 精品国产乱码久久久久久小说| 色哟哟·www| 在线观看免费高清a一片| 久久精品夜色国产| 日韩视频在线欧美| 国产精品.久久久| 欧美亚洲 丝袜 人妻 在线| 各种免费的搞黄视频| 777米奇影视久久| 精品卡一卡二卡四卡免费| 曰老女人黄片| 日日啪夜夜爽| 女人久久www免费人成看片| 街头女战士在线观看网站| 国产精品三级大全| 午夜福利视频精品| 国产精品一二三区在线看| 精品一品国产午夜福利视频| 精品国产乱码久久久久久男人| 日韩中字成人| 国产在线免费精品| 不卡视频在线观看欧美| 大香蕉久久成人网| 亚洲欧洲精品一区二区精品久久久 | 91成人精品电影| 国产精品一区二区在线观看99| 国产 一区精品| 伦理电影免费视频| 亚洲四区av| 亚洲欧美色中文字幕在线| 99国产综合亚洲精品| 欧美日韩视频精品一区| 亚洲欧洲日产国产| 超碰成人久久| 日韩制服丝袜自拍偷拍| 欧美激情 高清一区二区三区| 日韩人妻精品一区2区三区| 日韩中文字幕视频在线看片| 啦啦啦在线免费观看视频4| 麻豆精品久久久久久蜜桃| 欧美在线黄色| 精品国产一区二区三区久久久樱花| 免费在线观看黄色视频的| 人人妻人人添人人爽欧美一区卜| 国产精品99久久99久久久不卡 | 亚洲图色成人| 久久鲁丝午夜福利片| 欧美亚洲 丝袜 人妻 在线| 欧美中文综合在线视频| 亚洲欧美清纯卡通| 精品亚洲乱码少妇综合久久| 国产一区二区三区综合在线观看| 香蕉国产在线看| 黑人欧美特级aaaaaa片| 九色亚洲精品在线播放| 精品国产一区二区三区四区第35| 熟女电影av网| 国产一区二区在线观看av| 亚洲成人手机| 纯流量卡能插随身wifi吗| 国产一区二区三区综合在线观看| 国产黄频视频在线观看| 国产精品 欧美亚洲| 亚洲欧美中文字幕日韩二区| 日韩精品免费视频一区二区三区| 天天操日日干夜夜撸| 色婷婷av一区二区三区视频| 夜夜骑夜夜射夜夜干| 欧美精品高潮呻吟av久久| 夫妻性生交免费视频一级片| 日本午夜av视频| 欧美日韩视频高清一区二区三区二| 曰老女人黄片| 国产色婷婷99| 午夜福利网站1000一区二区三区| av又黄又爽大尺度在线免费看| 免费日韩欧美在线观看| 欧美在线黄色| 天天操日日干夜夜撸| 国产精品欧美亚洲77777| 人妻少妇偷人精品九色| 日本欧美视频一区| av.在线天堂| 在线观看三级黄色| 亚洲国产最新在线播放| 国产亚洲av片在线观看秒播厂| 亚洲国产精品一区二区三区在线| 十分钟在线观看高清视频www| 亚洲国产欧美网| 久久精品国产a三级三级三级| 国产精品一区二区在线观看99| 欧美日韩精品成人综合77777| 精品少妇内射三级| 自线自在国产av| 色94色欧美一区二区| av片东京热男人的天堂| 男男h啪啪无遮挡| 精品亚洲成国产av| videosex国产| 99久久中文字幕三级久久日本| 亚洲国产最新在线播放| 高清视频免费观看一区二区| 老鸭窝网址在线观看| 亚洲欧洲精品一区二区精品久久久 | 亚洲男人天堂网一区| 国产精品偷伦视频观看了| 久久午夜综合久久蜜桃| 成人二区视频| 日韩熟女老妇一区二区性免费视频| 精品一区二区三区四区五区乱码 | 欧美成人午夜精品| av视频免费观看在线观看| 男女边摸边吃奶| 欧美中文综合在线视频| 人妻少妇偷人精品九色| 国产精品香港三级国产av潘金莲 | 99久国产av精品国产电影| 伊人久久国产一区二区| 大陆偷拍与自拍| 在线观看三级黄色| 丝袜脚勾引网站| 蜜桃在线观看..| 亚洲精品久久久久久婷婷小说| 国产日韩一区二区三区精品不卡| 人妻少妇偷人精品九色| 国产精品香港三级国产av潘金莲 | 中国国产av一级| 久久精品久久久久久久性| 国产色婷婷99| 日韩一本色道免费dvd| 夫妻午夜视频| 亚洲精品成人av观看孕妇| 亚洲人成网站在线观看播放| 一区二区三区精品91| 国产av国产精品国产| 在线观看人妻少妇| 精品国产露脸久久av麻豆| 在线观看三级黄色| 国产黄色视频一区二区在线观看| 久久这里有精品视频免费| 亚洲美女黄色视频免费看| 免费观看无遮挡的男女| 日日啪夜夜爽| 亚洲欧美一区二区三区黑人 | 国产在线一区二区三区精| 成人亚洲欧美一区二区av| 熟女少妇亚洲综合色aaa.| 一本大道久久a久久精品| 亚洲国产毛片av蜜桃av| 日韩人妻精品一区2区三区| 一级毛片电影观看| 亚洲欧美一区二区三区国产| 少妇被粗大的猛进出69影院| 麻豆av在线久日| 91在线精品国自产拍蜜月| 黑人巨大精品欧美一区二区蜜桃| 精品国产一区二区三区久久久樱花| 中文字幕色久视频| 亚洲欧美一区二区三区国产| 成人午夜精彩视频在线观看| 国产精品偷伦视频观看了| 美女福利国产在线| 天天躁狠狠躁夜夜躁狠狠躁| 丝袜人妻中文字幕| 青青草视频在线视频观看| 99久久人妻综合| 高清黄色对白视频在线免费看| 黑人猛操日本美女一级片| 超碰97精品在线观看| 99国产精品免费福利视频| 少妇人妻 视频| 一本大道久久a久久精品| 亚洲精华国产精华液的使用体验| 婷婷色综合大香蕉| 大片免费播放器 马上看| 国产免费又黄又爽又色| 久久99蜜桃精品久久| 男男h啪啪无遮挡| 亚洲三级黄色毛片| 制服诱惑二区| 一级毛片我不卡| 午夜免费观看性视频| 国产精品成人在线| 国产在视频线精品| 亚洲精品国产色婷婷电影| 久久人人爽人人片av| www.av在线官网国产| 欧美精品国产亚洲| 国产无遮挡羞羞视频在线观看| 波野结衣二区三区在线| 777米奇影视久久| 免费看av在线观看网站| 最近最新中文字幕大全免费视频 | 一区二区日韩欧美中文字幕| 久久精品亚洲av国产电影网| 国产免费一区二区三区四区乱码| 国产又色又爽无遮挡免| 久久精品国产自在天天线| 亚洲国产精品成人久久小说| 在线精品无人区一区二区三| 精品久久蜜臀av无| 国产精品一二三区在线看| 叶爱在线成人免费视频播放| 两个人看的免费小视频| 国产亚洲av片在线观看秒播厂| 一区二区三区精品91| 免费在线观看视频国产中文字幕亚洲 | 一二三四在线观看免费中文在| 夫妻性生交免费视频一级片| av天堂久久9| 国产又爽黄色视频| 午夜日本视频在线| 大片免费播放器 马上看| 欧美在线黄色| 一区二区av电影网| 97人妻天天添夜夜摸| 国产精品.久久久| 久久久久国产一级毛片高清牌| 久久女婷五月综合色啪小说| 建设人人有责人人尽责人人享有的| 18禁裸乳无遮挡动漫免费视频| 亚洲国产欧美在线一区| 一二三四中文在线观看免费高清| 亚洲人成网站在线观看播放| 国产精品秋霞免费鲁丝片| 免费黄色在线免费观看| 色网站视频免费| 久久精品久久精品一区二区三区| 国产免费又黄又爽又色| 最新中文字幕久久久久| 黄色视频在线播放观看不卡| 亚洲成人av在线免费| 国产亚洲最大av| 国产一区有黄有色的免费视频| 精品国产乱码久久久久久小说| 激情视频va一区二区三区| 日本-黄色视频高清免费观看| 国产成人av激情在线播放| 极品少妇高潮喷水抽搐| 欧美亚洲 丝袜 人妻 在线| 人妻人人澡人人爽人人| 国产淫语在线视频| 久久这里只有精品19| 一级a爱视频在线免费观看| 国产高清不卡午夜福利| 人人澡人人妻人| 亚洲第一av免费看| 亚洲美女视频黄频| 国产亚洲一区二区精品| 精品福利永久在线观看| 午夜日韩欧美国产| 赤兔流量卡办理| 久久久久国产一级毛片高清牌| 美女视频免费永久观看网站| 男的添女的下面高潮视频| 免费观看无遮挡的男女| 亚洲精华国产精华液的使用体验| 日本免费在线观看一区| av国产精品久久久久影院| 欧美激情 高清一区二区三区| 中文精品一卡2卡3卡4更新| 日韩制服丝袜自拍偷拍| 欧美亚洲 丝袜 人妻 在线| 国产老妇伦熟女老妇高清| 亚洲成av片中文字幕在线观看 | 国语对白做爰xxxⅹ性视频网站| 亚洲第一青青草原| 国产成人精品婷婷| 国产野战对白在线观看| 婷婷色综合大香蕉| 一边摸一边做爽爽视频免费| 精品酒店卫生间| 综合色丁香网| 宅男免费午夜| 精品久久久久久电影网| 人妻系列 视频| 精品福利永久在线观看| 精品卡一卡二卡四卡免费| 日本黄色日本黄色录像| 国产精品 欧美亚洲| 在线观看一区二区三区激情| 男人操女人黄网站| 建设人人有责人人尽责人人享有的| 日本免费在线观看一区| av卡一久久| 欧美激情极品国产一区二区三区| 中文字幕最新亚洲高清| 亚洲国产看品久久| 久久久久国产一级毛片高清牌| 伊人久久国产一区二区| 黄片无遮挡物在线观看| 国产片特级美女逼逼视频| 成年美女黄网站色视频大全免费| 亚洲国产欧美网| 日韩熟女老妇一区二区性免费视频| 久久青草综合色| 精品少妇黑人巨大在线播放| 成人漫画全彩无遮挡| 久久毛片免费看一区二区三区| 国产精品人妻久久久影院| av在线app专区| 国产国语露脸激情在线看| 中文字幕亚洲精品专区| 哪个播放器可以免费观看大片| 日韩一区二区三区影片| 亚洲第一区二区三区不卡| videossex国产| 国产在线免费精品| av国产精品久久久久影院| 亚洲,欧美,日韩| av天堂久久9| 免费看不卡的av| 亚洲欧洲精品一区二区精品久久久 | 精品一区二区三区四区五区乱码 | 亚洲欧美一区二区三区黑人 | 午夜福利在线免费观看网站| 国产1区2区3区精品| 欧美+日韩+精品| 在线观看国产h片| 日韩在线高清观看一区二区三区| 男女国产视频网站| 久久久精品免费免费高清| 国产激情久久老熟女| 久久国产精品男人的天堂亚洲| 日韩三级伦理在线观看| 波多野结衣av一区二区av| 99久久中文字幕三级久久日本| 十八禁高潮呻吟视频| 日韩一卡2卡3卡4卡2021年| 亚洲伊人色综图| 亚洲伊人久久精品综合| 久久精品国产鲁丝片午夜精品| 国产精品秋霞免费鲁丝片| 日日啪夜夜爽| 精品国产乱码久久久久久小说| 纯流量卡能插随身wifi吗| 99热国产这里只有精品6| 伊人久久大香线蕉亚洲五| 最新中文字幕久久久久| 国产欧美日韩一区二区三区在线| 精品国产乱码久久久久久小说| 啦啦啦啦在线视频资源| 欧美日韩精品成人综合77777| 在线亚洲精品国产二区图片欧美| 欧美日韩一级在线毛片| www.自偷自拍.com| 亚洲av男天堂| 日本91视频免费播放| 国产成人精品婷婷| 久久精品aⅴ一区二区三区四区 | 欧美人与性动交α欧美软件| 寂寞人妻少妇视频99o| 人妻少妇偷人精品九色| 亚洲av欧美aⅴ国产| 少妇人妻精品综合一区二区| 亚洲av国产av综合av卡| 天堂俺去俺来也www色官网| 91aial.com中文字幕在线观看| 国产一级毛片在线| 最新中文字幕久久久久| 丝袜喷水一区| 久久久精品94久久精品| 精品人妻在线不人妻| 久久久精品94久久精品| 欧美变态另类bdsm刘玥| 亚洲av电影在线观看一区二区三区| 蜜桃在线观看..| 人人妻人人添人人爽欧美一区卜| 欧美激情 高清一区二区三区| 91在线精品国自产拍蜜月| 色94色欧美一区二区| 9191精品国产免费久久| 熟女电影av网| 下体分泌物呈黄色| 亚洲美女黄色视频免费看| 狠狠精品人妻久久久久久综合| 久久精品久久久久久噜噜老黄| 欧美成人午夜免费资源| 菩萨蛮人人尽说江南好唐韦庄| 天天影视国产精品| 欧美国产精品一级二级三级| 亚洲国产精品999| 亚洲精品av麻豆狂野| 久久久久久久久久久免费av| 最黄视频免费看| 伊人久久大香线蕉亚洲五| 色婷婷av一区二区三区视频| 亚洲三级黄色毛片|