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

    某超燃沖壓發(fā)動(dòng)機(jī)尾噴管的三維流場(chǎng)特性

    2021-07-29 02:09:18王建明齊曉航欒思琦王成軍陳寶嶺
    科學(xué)技術(shù)與工程 2021年14期
    關(guān)鍵詞:羽流馬赫數(shù)激波

    王建明,齊曉航,欒思琦,王成軍,陳寶嶺

    (1.沈陽(yáng)航空航天大學(xué)航空發(fā)動(dòng)機(jī)學(xué)院,遼寧省航空推進(jìn)系統(tǒng)先進(jìn)測(cè)試技術(shù)重點(diǎn)實(shí)驗(yàn)室,沈陽(yáng) 110136;2.沈陽(yáng)航天新樂(lè)有限責(zé)任公司,沈陽(yáng) 110034)

    超燃沖壓發(fā)動(dòng)機(jī)作為高超聲速飛行器的動(dòng)力裝置,相比于傳統(tǒng)的渦噴和渦扇發(fā)動(dòng)機(jī),它在高馬赫數(shù)下能夠擁有較大的推重比[1]。它的結(jié)構(gòu)簡(jiǎn)單,由進(jìn)氣道、隔離段、燃燒室和尾噴管四部分構(gòu)成,構(gòu)型上采用飛行器機(jī)體一體化設(shè)計(jì),將飛行器的下表面作為進(jìn)氣道和尾噴管的上壁面,以減少迎風(fēng)面積和重量。超燃沖壓發(fā)動(dòng)機(jī)尾噴管是提供發(fā)動(dòng)機(jī)推力的主要部件,它產(chǎn)生推力的原理是:燃燒室高溫高壓氣體在噴管入口處膨脹加速,利用尾噴管上壁面作為外膨脹面以增大燃?xì)獾呐蛎浢娣e[2],讓出口處氣流和進(jìn)口處產(chǎn)生動(dòng)量差,從而使發(fā)動(dòng)機(jī)獲得凈推力。

    影響超燃沖壓發(fā)動(dòng)機(jī)尾噴管性能的因素眾多,對(duì)尾噴管性能影響規(guī)律的研究是近年來(lái)該領(lǐng)域的重要研究方向。Huang等[3]做了氫燃料超燃沖壓發(fā)動(dòng)機(jī)噴管的化學(xué)非平衡流動(dòng)分析,方法是利用CHEMKIN軟件對(duì)Hyshot型超燃沖壓發(fā)動(dòng)機(jī)不同工況下噴管內(nèi)的化學(xué)非平衡現(xiàn)象進(jìn)行了數(shù)值分析。通過(guò)對(duì)凍結(jié)流、平衡流和非平衡流三種化學(xué)狀態(tài)的測(cè)試和比較,論證了化學(xué)反應(yīng)對(duì)噴嘴流場(chǎng)的影響。結(jié)果表明:真實(shí)非平衡流模擬計(jì)算結(jié)果介于凍結(jié)流和平衡流兩種極限情況之間,更接近于凍結(jié)流。在噴管入口燃?xì)馔耆紵那闆r下,自由基復(fù)合反應(yīng)會(huì)釋放出大量熱量,這些熱量主要用于提高氣體溫度,對(duì)推力作用很小。在噴管入口不完全燃燒情況下,噴管內(nèi)既有燃燒反應(yīng),又有自由基復(fù)合,反應(yīng)放熱對(duì)推力的影響取決于噴管進(jìn)口處燃燒的完全程度,推力相比于凍結(jié)流則大幅度提高。Mo等[4]做了非均勻進(jìn)氣超燃沖壓發(fā)動(dòng)機(jī)噴管設(shè)計(jì)與冷態(tài)流動(dòng)的試驗(yàn)。提出了一種考慮非均勻入流影響的非對(duì)稱(chēng)超燃沖壓發(fā)動(dòng)機(jī)噴管設(shè)計(jì)計(jì)算程序。給出了給定非均勻馬赫數(shù)分布和相應(yīng)的質(zhì)量加權(quán)平均均勻馬赫數(shù)分布的典型設(shè)計(jì)實(shí)例。然后,通過(guò)計(jì)算流體力學(xué)分析和冷態(tài)流動(dòng)實(shí)驗(yàn)測(cè)量,對(duì)非均勻設(shè)計(jì)的噴管性能改善進(jìn)行了量化。計(jì)算和實(shí)驗(yàn)的結(jié)果表明,非均勻馬赫數(shù)噴管設(shè)計(jì)法設(shè)計(jì)出的噴管較均勻馬赫數(shù)方法下整體性能更優(yōu)。Lü等[5]做了考慮橫向展開(kāi)和幾何約束的三維超燃沖壓發(fā)動(dòng)機(jī)噴管設(shè)計(jì)的研究,提出了一種基于準(zhǔn)二維超音速流動(dòng)和最大推力理論的考慮橫向展開(kāi)和幾何約束的三維噴管設(shè)計(jì)方法。為了生成三維噴管的外形,采用特征線法計(jì)算無(wú)粘流場(chǎng)以及參考溫度法修正附面層厚度,并用計(jì)算流體力學(xué)方法計(jì)算了噴管的氣動(dòng)性能。結(jié)果表明:初始圓弧半徑對(duì)軸向推力系數(shù)影響較小,而橫向輪廓、下整流罩長(zhǎng)度和初始膨脹角的變化對(duì)軸向推力系數(shù)有顯著影響。并和用流線追蹤法設(shè)計(jì)的三維噴管進(jìn)行了對(duì)比研究,驗(yàn)證了新方法的優(yōu)越性。黃偉等[6]也用數(shù)值模擬的方法論述了尾噴管構(gòu)型對(duì)性能的影響,研究表明對(duì)于俯仰力矩,上壁面切線角越小越好,外罩長(zhǎng)度越長(zhǎng)越好,而對(duì)于升力則剛好相反,所以在設(shè)計(jì)尾噴管構(gòu)型時(shí)要充分平衡升力和俯仰力矩的需求。Zhang等[7]做了非均勻入口馬赫數(shù)的超燃沖壓發(fā)動(dòng)機(jī)噴管三維流場(chǎng)仿真實(shí)驗(yàn),分析某種構(gòu)型的噴管內(nèi)的馬赫數(shù)和壓力等值線,結(jié)果表明相比于均勻馬赫數(shù)入口條件,非均勻入口條件下的推力和陀螺力矩都會(huì)降低。

    目前關(guān)于超燃沖壓發(fā)動(dòng)機(jī)尾噴管流場(chǎng)的數(shù)值模擬大多都是基于二維流場(chǎng),討論的內(nèi)容也是以入口參數(shù)對(duì)發(fā)動(dòng)機(jī)推力、升力等性能參數(shù)的影響為主,很少有對(duì)噴管流場(chǎng)三維結(jié)構(gòu)的研究,而噴管流場(chǎng)的結(jié)構(gòu)特性也是影響噴管性能的重要因素?,F(xiàn)著重于研究噴管三維流場(chǎng)的結(jié)構(gòu),通過(guò)觀察各個(gè)截面上密度、馬赫數(shù)等值線、渦量及速度矢量線圖來(lái)分析三維流場(chǎng)中的波系和渦系的分布特性。

    1 計(jì)算模型與結(jié)果驗(yàn)證

    1.1 計(jì)算模型和網(wǎng)格

    噴管計(jì)算模型選自文獻(xiàn)[8],三維模型及結(jié)構(gòu)組成如圖1(a)所示,二維結(jié)構(gòu)示意圖如圖1(b)所示,主要設(shè)計(jì)參數(shù)及尺寸如表1所示,x是主流方向。

    計(jì)算網(wǎng)格如圖2所示,圖2(a)是整個(gè)計(jì)算域的網(wǎng)格,圖2(b)是噴管處的網(wǎng)格。在噴管近壁面處的網(wǎng)格均進(jìn)行了加密處理,網(wǎng)格原點(diǎn)取在噴管入口的左下角。網(wǎng)格大部分由O型網(wǎng)格組成,局部位置采用非結(jié)構(gòu)網(wǎng)格,近壁面第一層網(wǎng)格高度設(shè)置為0.005 mm,y+<1,網(wǎng)格總數(shù)約230萬(wàn)。

    圖1 噴管模型圖Fig.1 Model of nozzle

    表1 噴管主要設(shè)計(jì)參數(shù)Table 1 Main design parameters of nozzle

    1.2 數(shù)值方法

    流動(dòng)可視為理想氣體可壓縮定常流動(dòng)。在三維坐標(biāo)系下微分形式的N-S方程如下:

    (1)

    (2)

    (3)

    式(1)~式(3)分別是質(zhì)量、動(dòng)量和能量的守恒形式方程[9-10]。式中:t為時(shí)間;i、j代表方向;ρ為流體密度;p為靜壓;u為流體速度;E為流體微團(tuán)總能量;T為溫度;gi和Fi分別為i方向上控制體所受的重力體積力和外部體積力;τij為應(yīng)力張量;k為有效熱傳導(dǎo)系數(shù);Sh為其他熱源項(xiàng)。

    圖2 計(jì)算域與網(wǎng)格Fig.2 Computing domain and grid

    現(xiàn)選用基于SST(shear-stress transport)k-ω湍流模型。其中湍動(dòng)能k的輸運(yùn)方程為

    (4)

    特定湍流耗散率ω的輸運(yùn)方程為

    (5)

    式中:Γk和Γω分別為k和ω的有效耗散系數(shù);Gk和Gω分別為k和ω的產(chǎn)生項(xiàng);Yk和Yω分別為k和ω由于湍流引起的耗散項(xiàng);Xω是交叉擴(kuò)散項(xiàng)。湍動(dòng)能的耗散項(xiàng)Yk可表示為

    Yk=ρβkωFDES

    (6)

    式(6)中:FDES為湍流強(qiáng)度項(xiàng)。

    數(shù)值模擬用FLUENT軟件,求解器選擇壓力基,計(jì)算方法采用Coupled算法,選取Green-Gause節(jié)點(diǎn)格式,密度、壓力等參數(shù)均選取二階迎風(fēng)格式做定常計(jì)算,共迭代14 000步。設(shè)置流體為理想氣體,流體黏性采用Sutherland公式。所有壁面皆為絕熱、無(wú)滑移壁面。根據(jù)參考文獻(xiàn)設(shè)定邊界條件為噴管入口靜壓35 574 Pa,靜溫298 K,馬赫數(shù)2.5,自由流靜壓2 940 Pa,總溫673 K,馬赫數(shù)7.1。

    1.3 網(wǎng)格穩(wěn)定性和計(jì)算結(jié)果驗(yàn)證

    為了保證計(jì)算精度的基礎(chǔ)上加快計(jì)算速度,對(duì)網(wǎng)格的獨(dú)立性進(jìn)行驗(yàn)證,分別計(jì)算了180萬(wàn)、200萬(wàn)、230萬(wàn)和270萬(wàn)的網(wǎng)格,并對(duì)計(jì)算結(jié)果進(jìn)行無(wú)量綱化,計(jì)算結(jié)果如圖3所示,橫坐標(biāo)為噴管上壁面中心線位置坐標(biāo),縱坐標(biāo)是噴管上壁面中心線上壓力無(wú)量綱化值,根據(jù)文獻(xiàn)[8],P是上壁面中心線的壓力值,Pc是試驗(yàn)臺(tái)模擬的燃燒室出口氣體的壓力值,按照文獻(xiàn)[8]中無(wú)量綱化的方法取(P/Pc)以10為底的對(duì)數(shù)值。如圖3可知當(dāng)網(wǎng)格在230萬(wàn)時(shí)結(jié)果已經(jīng)不再隨網(wǎng)格量的增加而改變,因此選取230萬(wàn)的網(wǎng)格進(jìn)行計(jì)算。

    仿真結(jié)果與文獻(xiàn)[8]的實(shí)驗(yàn)結(jié)果的對(duì)比圖如圖4所示,圖4中的橫坐標(biāo)為噴管上壁面中心線坐標(biāo),縱坐標(biāo)是噴管上壁面中心線處壓力的無(wú)量綱化的數(shù)值,圖4表明仿真的曲線結(jié)果與實(shí)驗(yàn)測(cè)得的九個(gè)測(cè)量點(diǎn)的結(jié)果基本吻合,進(jìn)一步說(shuō)明了仿真結(jié)果的準(zhǔn)確性和可信性。

    圖3 網(wǎng)格獨(dú)立性驗(yàn)證Fig.3 Grid independence verification

    圖4 數(shù)值結(jié)果與實(shí)驗(yàn)值對(duì)比Fig.4 Comparison of numerical and experimental values

    2 三維流場(chǎng)結(jié)構(gòu)的結(jié)果與分析

    2.1 三維流場(chǎng)中膨脹波分布

    為了觀察三維流場(chǎng)結(jié)構(gòu),除了截取噴管對(duì)稱(chēng)面(z=0.025 m)以外,另取y=0.005 7 m、x=0.02 m、x=0.13 m三處截面,截面位置示意圖如圖5所示。

    流場(chǎng)各個(gè)截面的密度分布圖如圖6所示。根據(jù)流場(chǎng)中密度分布情況可以確定膨脹波的位置。圖6(a)中,發(fā)現(xiàn)噴管內(nèi)存在沿流向密度下降明顯的邊界,流場(chǎng)中密度下降、壓力必然下降,因此該處即為管內(nèi)膨脹波的位置。圖6(a)中可以看到這兩道膨脹波的起始位置分別在下壁面噴管入口處和上壁面內(nèi)噴管出口處,這些位置產(chǎn)生膨脹波的原因是這兩處都是超聲速氣流發(fā)生外折的位置。根據(jù)圖6(a)中所示兩道膨脹波會(huì)在內(nèi)噴管出口處附近交匯,管內(nèi)氣流在經(jīng)過(guò)這兩道膨脹波時(shí)膨脹加速,進(jìn)而產(chǎn)生推力和升力。圖6(b)中可以看到圖6(a)中噴管內(nèi)膨脹波在y方向上的位置,同時(shí)可以觀察到產(chǎn)生于兩側(cè)壁面出口處的膨脹波,側(cè)壁面產(chǎn)生的兩道膨脹波均存在內(nèi)折角。根據(jù)圖6(c),可以得知膨脹波之所以發(fā)生內(nèi)折,是因?yàn)閭?cè)壁尾端的射流邊界與側(cè)壁形成了外鈍角,而且射流邊界處的密度小于噴管出口的密度,即射流邊界處壓強(qiáng)小于入口處壓強(qiáng),因此會(huì)在兩個(gè)側(cè)壁尾端產(chǎn)生內(nèi)折的膨脹波。

    圖5 截面位置示意圖Fig.5 Schematic diagram of section position

    2.2 三維流場(chǎng)激波和剪切層分布

    圖7是流場(chǎng)中各個(gè)截面的馬赫數(shù)等值線,由于激波處馬赫數(shù)等值線會(huì)發(fā)生匯聚,因此根據(jù)馬赫數(shù)等值線匯聚的位置可以確定激波的大致位置。圖7(a)中可以觀察到上、下壁面出口處的羽流激波,羽流激波的存在是因?yàn)榇嬖诟呖沼鹆餍?yīng)。高空羽流效應(yīng)的產(chǎn)生是因?yàn)槌曀賴娏麟x開(kāi)噴管時(shí)會(huì)在外部的超聲速介質(zhì)中流動(dòng),二者的氣流參數(shù)不同,此時(shí)噴管處的超聲速氣流在模擬狀態(tài)下的超聲速介質(zhì)中擴(kuò)散,就會(huì)產(chǎn)生類(lèi)似于發(fā)動(dòng)機(jī)在高空飛行環(huán)境下所產(chǎn)生的高空羽流效應(yīng),高空羽流效應(yīng)會(huì)導(dǎo)致飛行器產(chǎn)生附體斜激波,這類(lèi)激波被稱(chēng)為羽流激波[11-12]。在上壁面尾部,由于上方自由流穿過(guò)羽流激波,導(dǎo)致上壁面尾部附近的壓力明顯增高,進(jìn)而導(dǎo)致此處附近的自由流對(duì)主流產(chǎn)生壓縮,進(jìn)而形成一道由壓力決定的管內(nèi)斜激波。圖7(b)中,可以看到噴管側(cè)壁面同樣存在羽流激波。因此可以推斷,外噴管流場(chǎng)周?chē)挠鹆骷げㄊ浅虱h(huán)狀分布在噴管流場(chǎng)周?chē)?。并且根?jù)圖中馬赫數(shù)等值線的分布情況可以看出隨著流動(dòng)的發(fā)展,激波處馬赫數(shù)等值線變得越來(lái)越稀疏,這表明激波會(huì)隨著流動(dòng)的發(fā)展強(qiáng)度逐漸減弱。

    觀察馬赫數(shù)等值線圖,發(fā)現(xiàn)除了激波位置外,仍存在其他的等值線發(fā)生匯聚的位置。實(shí)際上流場(chǎng)中除了激波處,剪切層處的馬赫數(shù)等值線也會(huì)發(fā)生匯聚,因此需要判定這些位置上是否存在剪切層。剪切層最明顯的特征是它的內(nèi)部渦量很高,大多數(shù)研究者根據(jù)這個(gè)特征來(lái)判定流場(chǎng)中是否存在剪切層。渦量云圖如圖8所示,根據(jù)圖8(a),發(fā)現(xiàn)了噴管上、下壁面尾端的剪切層包含z方向的渦。圖8(b)中可以觀察到噴管主流兩側(cè)剪切層則包含y方向的渦。

    圖6 流場(chǎng)密度云圖Fig.6 Cloud map of density

    圖7 流場(chǎng)馬赫數(shù)等值線圖Fig.7 Contour map of Mach number

    2.3 三維流場(chǎng)渦流分布

    流向的渦流比較難捕捉,為了確定流場(chǎng)中流向渦的具體位置,做出了各個(gè)截面上的速度矢量線圖,通過(guò)觀察各個(gè)截面的速度矢量線圖來(lái)確定渦流的位置,矢量圖中速度的大小值取基于該截面下的y方向和z方向的合成速度的值。根據(jù)圖9(a),發(fā)現(xiàn)流場(chǎng)中有渦流的存在,由于該截面方向垂直于主流,因此這些渦流為流向渦,根據(jù)矢量線的分布,發(fā)現(xiàn)渦產(chǎn)生于內(nèi)噴管出口附近。渦流的形成原因是由于內(nèi)噴管出口處產(chǎn)生的高空羽流效應(yīng),使得主流的超聲速氣流受到羽流的影響,氣流會(huì)沿橫向分離并發(fā)生卷折,進(jìn)而形成流向渦。圖9(b)中,上壁面出口處也觀察到了流向渦的存在。它的形成原因是由于上壁面的邊界層存在橫向分離的趨勢(shì),這一趨勢(shì)導(dǎo)致氣流離開(kāi)上壁面后,進(jìn)而在尾流中生成流向渦[13-14]。

    圖8 流場(chǎng)渦量云圖Fig.8 Vorticity cloud map of flow field

    2.4 三維流場(chǎng)主要結(jié)構(gòu)示意圖

    圖10是綜合前面結(jié)果得到了三維流場(chǎng)結(jié)構(gòu)示意圖。圖10結(jié)果與文獻(xiàn)[8]的實(shí)驗(yàn)結(jié)果相比,文獻(xiàn)[8]所得到的流場(chǎng)結(jié)構(gòu)示意圖并沒(méi)有討論流場(chǎng)中的膨脹波,討論的流場(chǎng)結(jié)構(gòu)也都是基于二維截面下的。本文揭示的流動(dòng)結(jié)構(gòu)則是基于三維流場(chǎng),比文獻(xiàn)[8]中的更全面,并且對(duì)所有物理現(xiàn)象的產(chǎn)生原因加以論述。

    3 結(jié)論

    通過(guò)數(shù)值模擬的方法對(duì)超燃沖壓發(fā)動(dòng)機(jī)三維定常流場(chǎng)結(jié)構(gòu)進(jìn)行分析,結(jié)論如下。

    圖9 速度矢量線圖Fig.9 Velocity vector line diagram

    (1)噴管內(nèi)存在膨脹波,管內(nèi)的膨脹波會(huì)發(fā)生交匯形成膨脹波面。同時(shí)側(cè)壁也產(chǎn)生內(nèi)折的膨脹波,該處膨脹波會(huì)隨著氣流膨脹,強(qiáng)度逐漸變?nèi)酢?/p>

    (2)噴管外流場(chǎng)周?chē)嬖谟鹆骷げê图羟袑?,羽流激波?huì)形成一個(gè)激波面分布在外流場(chǎng)周?chē)?,剪切層則分布在羽流激波和主流之間。在上壁面出口處存在羽流激波和因邊界層分離產(chǎn)生的斜激波,并且在這兩道激波之間也存在剪切層。

    (3)內(nèi)噴管出口處氣流會(huì)在噴管主流兩側(cè)形成流向渦。外噴管出口處氣流也會(huì)在上壁面出口兩側(cè)形成渦流。內(nèi)噴管出口附近的流向渦主要是受到羽流效應(yīng)的影響,尾流中的流向渦則是由于上壁面氣流沿橫向分離造成的。

    猜你喜歡
    羽流馬赫數(shù)激波
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    水下羽流追蹤方法研究進(jìn)展
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    載荷分布對(duì)可控?cái)U(kuò)散葉型性能的影響
    一種基于聚類(lèi)分析的二維激波模式識(shí)別算法
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    斜激波入射V形鈍前緣溢流口激波干擾研究
    適于可壓縮多尺度流動(dòng)的緊致型激波捕捉格式
    水下管道向下泄漏的羽/射流特性
    室內(nèi)多股羽流混合運(yùn)動(dòng)機(jī)理模型研究進(jìn)展分析
    欧美黑人欧美精品刺激| 亚洲国产精品成人综合色| 老鸭窝网址在线观看| 国产精品嫩草影院av在线观看 | 国产一级毛片七仙女欲春2| 成人国产一区最新在线观看| 亚洲av成人av| 国产欧美日韩一区二区精品| 蜜桃亚洲精品一区二区三区| 看黄色毛片网站| 久久精品91蜜桃| avwww免费| 波多野结衣高清作品| 欧美成人性av电影在线观看| 性欧美人与动物交配| 国产精品野战在线观看| 日韩人妻高清精品专区| 亚洲18禁久久av| 亚洲aⅴ乱码一区二区在线播放| 欧美色欧美亚洲另类二区| 亚洲天堂国产精品一区在线| 国产视频一区二区在线看| 国产午夜精品久久久久久一区二区三区 | 欧美+亚洲+日韩+国产| 欧美黑人欧美精品刺激| 久久午夜亚洲精品久久| 怎么达到女性高潮| 黄色日韩在线| 中文字幕人成人乱码亚洲影| 欧美在线一区亚洲| 中出人妻视频一区二区| 搞女人的毛片| 午夜福利18| 欧美日韩中文字幕国产精品一区二区三区| 日韩亚洲欧美综合| 日韩大尺度精品在线看网址| 欧美一区二区精品小视频在线| 在线免费观看的www视频| 99久久无色码亚洲精品果冻| 99久久成人亚洲精品观看| 熟妇人妻久久中文字幕3abv| 中文字幕av成人在线电影| 国产高清视频在线观看网站| 日韩欧美精品免费久久 | 99国产精品一区二区三区| 人人妻,人人澡人人爽秒播| 老汉色av国产亚洲站长工具| 国产欧美日韩一区二区精品| 国产高清三级在线| 身体一侧抽搐| 亚洲av美国av| 亚洲最大成人手机在线| 国产aⅴ精品一区二区三区波| 韩国av一区二区三区四区| ponron亚洲| 丰满人妻一区二区三区视频av | 中文资源天堂在线| 国产精品影院久久| 亚洲狠狠婷婷综合久久图片| 成人鲁丝片一二三区免费| 欧美高清成人免费视频www| 国内少妇人妻偷人精品xxx网站| 久久中文看片网| 18美女黄网站色大片免费观看| 亚洲一区高清亚洲精品| 老司机午夜福利在线观看视频| 国产精品野战在线观看| 此物有八面人人有两片| 国产精品精品国产色婷婷| 操出白浆在线播放| a级一级毛片免费在线观看| 国产精品98久久久久久宅男小说| 99久久精品热视频| 欧美乱妇无乱码| 国产高清videossex| 日韩精品中文字幕看吧| 精品乱码久久久久久99久播| 1024手机看黄色片| 色综合婷婷激情| 97碰自拍视频| 天美传媒精品一区二区| e午夜精品久久久久久久| 91在线观看av| eeuss影院久久| 成人午夜高清在线视频| 亚洲国产精品久久男人天堂| 欧美三级亚洲精品| 男女床上黄色一级片免费看| bbb黄色大片| 一区二区三区国产精品乱码| 久久国产乱子伦精品免费另类| АⅤ资源中文在线天堂| 99久久无色码亚洲精品果冻| 深爱激情五月婷婷| 国产成年人精品一区二区| 国产免费av片在线观看野外av| 男人舔女人下体高潮全视频| 美女大奶头视频| 国产91精品成人一区二区三区| 床上黄色一级片| 天天添夜夜摸| 女人十人毛片免费观看3o分钟| 久久久久九九精品影院| 别揉我奶头~嗯~啊~动态视频| 亚洲精品456在线播放app | 99国产极品粉嫩在线观看| 精品国产三级普通话版| 国产国拍精品亚洲av在线观看 | 婷婷丁香在线五月| 久久久精品大字幕| 成人特级黄色片久久久久久久| 九九久久精品国产亚洲av麻豆| 91久久精品国产一区二区成人 | 亚洲 国产 在线| 久久香蕉国产精品| 中文字幕精品亚洲无线码一区| 国产激情欧美一区二区| 久久欧美精品欧美久久欧美| 久久6这里有精品| 国产成人啪精品午夜网站| 三级男女做爰猛烈吃奶摸视频| 亚洲激情在线av| 久久九九热精品免费| 制服人妻中文乱码| 久久精品人妻少妇| 超碰av人人做人人爽久久 | 日韩欧美国产一区二区入口| 一夜夜www| 亚洲国产精品久久男人天堂| 欧美一区二区精品小视频在线| or卡值多少钱| 蜜桃亚洲精品一区二区三区| 国产极品精品免费视频能看的| 天堂√8在线中文| 少妇人妻精品综合一区二区 | 欧美区成人在线视频| 亚洲,欧美精品.| 日韩亚洲欧美综合| 久久精品国产亚洲av涩爱 | 欧美黑人巨大hd| 国产高清videossex| 999久久久精品免费观看国产| АⅤ资源中文在线天堂| 免费在线观看亚洲国产| 激情在线观看视频在线高清| 99热只有精品国产| 亚洲人成伊人成综合网2020| 国产黄片美女视频| 深夜精品福利| 91麻豆av在线| 久久香蕉精品热| 亚洲熟妇中文字幕五十中出| 精品熟女少妇八av免费久了| 又黄又爽又免费观看的视频| 天堂√8在线中文| 亚洲人成伊人成综合网2020| 精品不卡国产一区二区三区| 淫秽高清视频在线观看| 欧美中文综合在线视频| 久久精品91蜜桃| 欧美日韩中文字幕国产精品一区二区三区| 精品99又大又爽又粗少妇毛片 | 欧美绝顶高潮抽搐喷水| 欧美在线黄色| 在线观看日韩欧美| 国产精品乱码一区二三区的特点| 国产精品一区二区三区四区久久| 最新在线观看一区二区三区| 日本成人三级电影网站| 国产精品免费一区二区三区在线| 久久久色成人| 国产精品一区二区免费欧美| 日本五十路高清| 亚洲欧美日韩卡通动漫| 国产成人a区在线观看| 亚洲精品亚洲一区二区| 91在线精品国自产拍蜜月 | 亚洲熟妇熟女久久| 欧美极品一区二区三区四区| 每晚都被弄得嗷嗷叫到高潮| 51午夜福利影视在线观看| 黄色丝袜av网址大全| 色哟哟哟哟哟哟| 亚洲无线在线观看| 老司机午夜福利在线观看视频| 女同久久另类99精品国产91| 国产精品女同一区二区软件 | 欧美色欧美亚洲另类二区| 日韩精品中文字幕看吧| 国产成人福利小说| 老司机福利观看| 亚洲中文字幕一区二区三区有码在线看| 中文字幕熟女人妻在线| 亚洲国产中文字幕在线视频| 国产男靠女视频免费网站| 免费搜索国产男女视频| 国产亚洲精品综合一区在线观看| 国产私拍福利视频在线观看| 国产69精品久久久久777片| 精品福利观看| 成人特级av手机在线观看| 老司机午夜福利在线观看视频| 夜夜看夜夜爽夜夜摸| a级毛片a级免费在线| 亚洲第一电影网av| 女人高潮潮喷娇喘18禁视频| 嫁个100分男人电影在线观看| 一级毛片高清免费大全| 国产精品一及| 亚洲第一欧美日韩一区二区三区| 国产午夜福利久久久久久| 深爱激情五月婷婷| 精品国内亚洲2022精品成人| 国产免费av片在线观看野外av| 女警被强在线播放| 在线天堂最新版资源| 偷拍熟女少妇极品色| 国产精品久久电影中文字幕| 亚洲一区二区三区色噜噜| 叶爱在线成人免费视频播放| 高潮久久久久久久久久久不卡| 亚洲美女黄片视频| 最近在线观看免费完整版| 老司机深夜福利视频在线观看| 香蕉av资源在线| 国产真实乱freesex| 日本撒尿小便嘘嘘汇集6| 国产精品1区2区在线观看.| 99在线人妻在线中文字幕| 性欧美人与动物交配| 亚洲电影在线观看av| 最近视频中文字幕2019在线8| 日韩有码中文字幕| 小说图片视频综合网站| 久久香蕉国产精品| 男女那种视频在线观看| 国产真人三级小视频在线观看| 成年人黄色毛片网站| 国产伦在线观看视频一区| 91麻豆精品激情在线观看国产| 97超级碰碰碰精品色视频在线观看| 成年版毛片免费区| 成人特级av手机在线观看| 国产精品永久免费网站| 亚洲av成人不卡在线观看播放网| 亚洲最大成人手机在线| www日本在线高清视频| 国产私拍福利视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 97碰自拍视频| 最后的刺客免费高清国语| 国产亚洲精品av在线| 一进一出好大好爽视频| 国产成人欧美在线观看| 老汉色∧v一级毛片| 亚洲黑人精品在线| 成人18禁在线播放| 国产精品1区2区在线观看.| 国产一区二区在线av高清观看| 在线观看日韩欧美| 欧美乱码精品一区二区三区| 91麻豆精品激情在线观看国产| 最新在线观看一区二区三区| 动漫黄色视频在线观看| 99久久无色码亚洲精品果冻| 亚洲无线观看免费| 中文字幕人妻丝袜一区二区| 亚洲国产中文字幕在线视频| 嫩草影视91久久| 亚洲精品久久国产高清桃花| 免费看光身美女| 亚洲,欧美精品.| 国内精品久久久久久久电影| 在线天堂最新版资源| 亚洲av熟女| 国产av麻豆久久久久久久| 国产伦精品一区二区三区视频9 | 欧美区成人在线视频| 亚洲性夜色夜夜综合| 亚洲av成人精品一区久久| 欧美av亚洲av综合av国产av| 无人区码免费观看不卡| 国产一区二区三区视频了| 狠狠狠狠99中文字幕| 99热精品在线国产| 亚洲18禁久久av| 亚洲国产中文字幕在线视频| 99国产极品粉嫩在线观看| 久久午夜亚洲精品久久| 午夜两性在线视频| 69人妻影院| 在线国产一区二区在线| 国产一区二区激情短视频| 国产综合懂色| 热99在线观看视频| 亚洲avbb在线观看| 身体一侧抽搐| 国内毛片毛片毛片毛片毛片| 真人做人爱边吃奶动态| 亚洲国产欧洲综合997久久,| 国产成人a区在线观看| 日韩欧美一区二区三区在线观看| svipshipincom国产片| 国产野战对白在线观看| 欧美色欧美亚洲另类二区| 国产爱豆传媒在线观看| 久久性视频一级片| 九九在线视频观看精品| 亚洲精华国产精华精| 老司机深夜福利视频在线观看| 国产精品爽爽va在线观看网站| 久久国产精品人妻蜜桃| 九九久久精品国产亚洲av麻豆| 性欧美人与动物交配| 国产乱人视频| 亚洲欧美日韩高清专用| 999久久久精品免费观看国产| 在线观看午夜福利视频| 亚洲av五月六月丁香网| 免费看美女性在线毛片视频| 18+在线观看网站| 美女高潮喷水抽搐中文字幕| 日本三级黄在线观看| 日韩 欧美 亚洲 中文字幕| 97超级碰碰碰精品色视频在线观看| 亚洲成人久久性| 18美女黄网站色大片免费观看| 欧美精品啪啪一区二区三区| 亚洲av第一区精品v没综合| 亚洲无线观看免费| 99国产综合亚洲精品| 精品久久久久久久人妻蜜臀av| 婷婷六月久久综合丁香| 国产久久久一区二区三区| 啦啦啦免费观看视频1| 狂野欧美白嫩少妇大欣赏| 全区人妻精品视频| 久久性视频一级片| 国产99白浆流出| 色精品久久人妻99蜜桃| 精品福利观看| 国产精品久久久久久人妻精品电影| 尤物成人国产欧美一区二区三区| 午夜a级毛片| 欧美3d第一页| 91av网一区二区| 美女 人体艺术 gogo| 日韩中文字幕欧美一区二区| 中亚洲国语对白在线视频| 老熟妇乱子伦视频在线观看| 亚洲七黄色美女视频| 欧美日本视频| 俄罗斯特黄特色一大片| 国产熟女xx| 哪里可以看免费的av片| 免费人成在线观看视频色| 美女黄网站色视频| 久久久久精品国产欧美久久久| 五月伊人婷婷丁香| 精品99又大又爽又粗少妇毛片 | 99视频精品全部免费 在线| 51午夜福利影视在线观看| 国产成+人综合+亚洲专区| 变态另类丝袜制服| 99热精品在线国产| 欧美性猛交黑人性爽| 怎么达到女性高潮| 国产欧美日韩精品亚洲av| 观看免费一级毛片| 欧美黄色片欧美黄色片| 精品人妻一区二区三区麻豆 | 免费av毛片视频| 亚洲精品乱码久久久v下载方式 | www.熟女人妻精品国产| 国内精品美女久久久久久| 精品久久久久久久末码| 国产精品亚洲一级av第二区| 18美女黄网站色大片免费观看| 久久久久免费精品人妻一区二区| 综合色av麻豆| 欧美乱码精品一区二区三区| 99精品欧美一区二区三区四区| 久久人妻av系列| 久久人人精品亚洲av| 女警被强在线播放| 日韩人妻高清精品专区| 欧美激情在线99| 精品久久久久久久久久免费视频| 神马国产精品三级电影在线观看| 美女黄网站色视频| 欧美一级毛片孕妇| 老熟妇仑乱视频hdxx| 麻豆国产av国片精品| 757午夜福利合集在线观看| 91在线精品国自产拍蜜月 | 在线a可以看的网站| 叶爱在线成人免费视频播放| 国产精品野战在线观看| 黄色视频,在线免费观看| 欧美高清成人免费视频www| 亚洲av美国av| 99热这里只有是精品50| 男女床上黄色一级片免费看| 精品人妻一区二区三区麻豆 | 国产成人a区在线观看| 色综合亚洲欧美另类图片| 欧美丝袜亚洲另类 | 99久久99久久久精品蜜桃| 高清日韩中文字幕在线| 国产在线精品亚洲第一网站| 99精品欧美一区二区三区四区| 国产精品久久电影中文字幕| 少妇丰满av| 国产日本99.免费观看| 我的老师免费观看完整版| 日韩欧美精品免费久久 | 国产精品98久久久久久宅男小说| or卡值多少钱| 啦啦啦韩国在线观看视频| 亚洲精品粉嫩美女一区| 国产伦一二天堂av在线观看| 香蕉丝袜av| 免费av观看视频| 国产精品女同一区二区软件 | 中文字幕熟女人妻在线| 国产欧美日韩精品一区二区| 国产真实伦视频高清在线观看 | 精品免费久久久久久久清纯| 亚洲午夜理论影院| 男女之事视频高清在线观看| 亚洲在线观看片| АⅤ资源中文在线天堂| 国产激情偷乱视频一区二区| www.999成人在线观看| 成人欧美大片| 成年女人永久免费观看视频| 丁香六月欧美| 欧美性猛交黑人性爽| 真实男女啪啪啪动态图| 国产精品三级大全| 久久精品国产综合久久久| 色综合亚洲欧美另类图片| xxx96com| www.色视频.com| 一本精品99久久精品77| 国产伦精品一区二区三区视频9 | 成人国产综合亚洲| 成人无遮挡网站| 老汉色∧v一级毛片| bbb黄色大片| 欧美av亚洲av综合av国产av| 日本三级黄在线观看| 欧美黑人欧美精品刺激| 高清在线国产一区| 免费高清视频大片| 久久精品亚洲精品国产色婷小说| 久久久久久久精品吃奶| 亚洲色图av天堂| 午夜免费激情av| 久久久久久人人人人人| 18禁黄网站禁片午夜丰满| 日韩成人在线观看一区二区三区| 校园春色视频在线观看| 热99在线观看视频| 夜夜躁狠狠躁天天躁| 亚洲欧美日韩无卡精品| www.熟女人妻精品国产| 波多野结衣高清作品| 天堂网av新在线| 国产精品日韩av在线免费观看| 免费观看精品视频网站| 极品教师在线免费播放| 变态另类成人亚洲欧美熟女| 欧美成人一区二区免费高清观看| 色哟哟哟哟哟哟| 亚洲成a人片在线一区二区| 人人妻,人人澡人人爽秒播| 欧美另类亚洲清纯唯美| 欧美在线黄色| 日本一本二区三区精品| 欧美性感艳星| av福利片在线观看| 日韩高清综合在线| 每晚都被弄得嗷嗷叫到高潮| 蜜桃亚洲精品一区二区三区| www日本黄色视频网| 国产在视频线在精品| 亚洲内射少妇av| 18+在线观看网站| 超碰av人人做人人爽久久 | 国产不卡一卡二| svipshipincom国产片| 精品久久久久久,| 亚洲第一电影网av| 亚洲精品国产精品久久久不卡| 免费观看精品视频网站| 免费观看人在逋| 三级男女做爰猛烈吃奶摸视频| 岛国在线免费视频观看| 亚洲精品在线美女| 又爽又黄无遮挡网站| 欧美又色又爽又黄视频| 亚洲国产中文字幕在线视频| 十八禁网站免费在线| 亚洲一区二区三区不卡视频| 亚洲国产精品成人综合色| 国产午夜福利久久久久久| 国产视频内射| 国产精品av视频在线免费观看| 最近最新中文字幕大全电影3| or卡值多少钱| 黄色日韩在线| 国产成人av教育| 国产男靠女视频免费网站| 高清在线国产一区| 国产一区二区亚洲精品在线观看| 久久精品国产综合久久久| 天天躁日日操中文字幕| 麻豆成人av在线观看| 亚洲国产精品成人综合色| 日本黄大片高清| 熟女人妻精品中文字幕| 亚洲精品日韩av片在线观看 | 欧美国产日韩亚洲一区| 变态另类丝袜制服| www日本在线高清视频| 欧美色视频一区免费| 免费看a级黄色片| 男插女下体视频免费在线播放| 黄色成人免费大全| 亚洲最大成人手机在线| 亚洲真实伦在线观看| 国产麻豆成人av免费视频| 欧美+亚洲+日韩+国产| xxxwww97欧美| 成人三级黄色视频| 亚洲av成人精品一区久久| av在线天堂中文字幕| 小蜜桃在线观看免费完整版高清| 露出奶头的视频| 国产精品98久久久久久宅男小说| 男女视频在线观看网站免费| 久久国产精品影院| 非洲黑人性xxxx精品又粗又长| 欧美一区二区精品小视频在线| 两人在一起打扑克的视频| 他把我摸到了高潮在线观看| 午夜免费男女啪啪视频观看 | 午夜日韩欧美国产| 亚洲精华国产精华精| 亚洲精品亚洲一区二区| 亚洲国产精品久久男人天堂| 久久香蕉国产精品| 老司机午夜十八禁免费视频| 午夜a级毛片| 少妇的逼水好多| 久久久久久久精品吃奶| 欧美三级亚洲精品| 99久国产av精品| 国产黄a三级三级三级人| 99久久99久久久精品蜜桃| 俺也久久电影网| 亚洲精品一区av在线观看| 亚洲精品456在线播放app | 麻豆久久精品国产亚洲av| 成人鲁丝片一二三区免费| 成人三级黄色视频| 高潮久久久久久久久久久不卡| 一个人免费在线观看电影| 99久久九九国产精品国产免费| 亚洲人成电影免费在线| av欧美777| 国产在线精品亚洲第一网站| 男女之事视频高清在线观看| 久久久久久人人人人人| 两人在一起打扑克的视频| 日本精品一区二区三区蜜桃| 色老头精品视频在线观看| 俺也久久电影网| 人人妻人人澡欧美一区二区| 真实男女啪啪啪动态图| 麻豆国产97在线/欧美| 国产精品1区2区在线观看.| 午夜激情福利司机影院| 亚洲国产精品sss在线观看| 成人一区二区视频在线观看| 51午夜福利影视在线观看| 亚洲精品一卡2卡三卡4卡5卡| 非洲黑人性xxxx精品又粗又长| 国产精品免费一区二区三区在线| 国产精品一区二区免费欧美| 51午夜福利影视在线观看| 亚洲专区中文字幕在线| 精品午夜福利视频在线观看一区| 欧美在线黄色| 国产综合懂色| 90打野战视频偷拍视频| 激情在线观看视频在线高清| 国产成人影院久久av| 老熟妇仑乱视频hdxx| 精品日产1卡2卡| 男人舔女人下体高潮全视频| 日本三级黄在线观看| 中文字幕熟女人妻在线| 亚洲 欧美 日韩 在线 免费| 校园春色视频在线观看| 午夜福利免费观看在线| 熟女电影av网| 欧美3d第一页| 国产精品久久久人人做人人爽| 欧美日本视频| 国产淫片久久久久久久久 | 嫩草影视91久久| 久久久久精品国产欧美久久久| 床上黄色一级片|