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

    航空發(fā)動機露天試車臺側(cè)風(fēng)模擬研究

    2022-01-18 07:31:06艾延廷朱建勇
    燃氣渦輪試驗與研究 2021年4期
    關(guān)鍵詞:風(fēng)筒總壓臺風(fēng)

    汪 才,艾延廷,朱建勇,張 巍,林 山,喬 黎

    (1.沈陽航空航天大學(xué),沈陽 110136;2.中國航發(fā)沈陽發(fā)動機研究所,沈陽 110015)

    1 引言

    航空發(fā)動機在地面露天試車臺試車期間常會受到側(cè)風(fēng)影響,導(dǎo)致發(fā)動機工作不穩(wěn)定,影響發(fā)動機推力測量和性能評定。發(fā)動機側(cè)風(fēng)試驗的目的是驗證側(cè)風(fēng)條件下進氣道與發(fā)動機的兼容性[1-3]。自然空間中的強側(cè)風(fēng)可能會造成發(fā)動機進口氣流畸變,導(dǎo)致壓氣機工作不穩(wěn)定甚至喘振。因此,有必要進行側(cè)風(fēng)試驗研究,以獲得不同風(fēng)向及風(fēng)速下發(fā)動機的響應(yīng)和進氣道的流動特性,同時確定發(fā)動機地面工作側(cè)風(fēng)邊界,進而優(yōu)化進氣道與發(fā)動機的兼容性,并為飛機在側(cè)風(fēng)條件下的安全運行提供保障[4-6]。

    截至目前,國外就側(cè)風(fēng)對露天試車臺的影響做了大量的試驗和研究工作。英國R·R 公司、美國P&W 公司和美國GE 公司各自建造了風(fēng)速可達30 m/s、空氣流量可達2 000 kg/s 的大型側(cè)風(fēng)風(fēng)源設(shè)備。其中,R·R公司和P&W公司設(shè)計了全方位的風(fēng)源設(shè)備,分析了各個角度對發(fā)動機形成的側(cè)風(fēng)條件對發(fā)動機進氣畸變的影響[7];GE 公司設(shè)計了從0°、45°、90°三個方位進行側(cè)風(fēng)試驗的風(fēng)源設(shè)備,給出了發(fā)動機若能順利通過這三個方向的側(cè)風(fēng)試驗,則其他方向的側(cè)風(fēng)試驗都能通過的結(jié)論。Tourrette 等[8]發(fā)展了一套數(shù)值求解Navier-Stokes方程的軟件,并通過與試驗結(jié)果對比,驗證了數(shù)值方法和湍流模型的可靠性;Brix等[9]開展了側(cè)風(fēng)條件下由于短艙進氣道的吸入作用產(chǎn)生地面渦的風(fēng)洞試驗研究;Trapp等[10]就側(cè)風(fēng)條件下飛機起飛時產(chǎn)生的進氣道渦流對發(fā)動機進口流場的影響進行了分析與研究。相比之下,國內(nèi)對航空發(fā)動機露天試驗臺側(cè)風(fēng)影響研究起步較晚,目前的側(cè)風(fēng)試驗主要是在發(fā)動機進口安裝多點壓力測頭,來測量有無側(cè)風(fēng)條件下發(fā)動機進口壓力場的變化。劉永泉等[11]通過仿真流場壓力分布模擬了某型發(fā)動機的進氣畸變;李志平等[12]研究了發(fā)動機受進氣畸變導(dǎo)致其工作穩(wěn)定性變化的規(guī)律;劉浩等[13]采用數(shù)值模擬方法,研究了側(cè)風(fēng)來流條件對地面渦形成和發(fā)展的影響;王寶坤等[14]提出了風(fēng)速標定方法和開展側(cè)風(fēng)試驗的程序,給出了側(cè)風(fēng)裝置工作狀態(tài)與試驗所需風(fēng)速之間的映射關(guān)系。由于露天試車臺側(cè)風(fēng)模擬的復(fù)雜性和有效性,我國對側(cè)風(fēng)模擬裝置試驗的研究深度仍有待提高。

    本文基于流體力學(xué)理論,運用ANSYS Workbench 對側(cè)風(fēng)裝置試驗過程進行仿真,研究露天試車臺側(cè)風(fēng)裝置內(nèi)部流場和下游出口處的壓力、流量、速度的分布及演化規(guī)律。此外,采用Realizablek-ε湍流模型進行數(shù)值模擬,獲得了風(fēng)機壓力與流量的對應(yīng)關(guān)系,揭示了側(cè)風(fēng)裝置中風(fēng)機使用數(shù)量對流場穩(wěn)定性的影響。

    2 計算方法

    2.1 流體控制方程

    在本文側(cè)風(fēng)數(shù)值模擬穩(wěn)態(tài)計算工況下,需要同時滿足的控制方程包括質(zhì)量守恒方程、動量守恒方程和能量方程,各方程的詳細推導(dǎo)過程可參考文獻[15]。

    (1) 質(zhì)量守恒方程

    式中:ρ為氣流密度,u,v,ω分別為x、y、z方向的速度分量。

    (2) 動量守恒方程

    x、y、z三個方向的動量守恒方程為:

    式中:p為微元體壓力,τxx、τxy、τxz分別為因黏性作用在微元體表面上的黏性應(yīng)力分量,F(xiàn)x、Fy、Fz分別為微元體上的體力。

    對于牛頓流體,黏性應(yīng)力與流體的變形率成比例,有:

    式中:μ為動力黏度;λ為第二黏度,一般取-。

    將(3)式代入(2)式得,

    式中:Su、Sv、Sw為動量方程廣義源項,其中Su=Fx+sx,Sv=Fy+sy,Sw=Fz+sz。

    上式中sx、sy、sz的表達式如下:

    一般來講sx、sy、sz是小量,對于黏性為常數(shù)的不可壓流體,sx=sy=sz=0。

    (3) 能量守恒方程

    式中:e為物體內(nèi)能,T為溫度,q為熱流量,k為熱傳導(dǎo)系數(shù)。

    2.2 航空發(fā)動機側(cè)風(fēng)裝置模型

    不同露天試車臺所使用的側(cè)風(fēng)裝置結(jié)構(gòu)不盡相同,側(cè)風(fēng)裝置的模擬風(fēng)源可以由一臺渦軸發(fā)動機在減速裝置作用下連帶的一組飛機螺旋槳組成,也可以由安裝在風(fēng)源筒體同一平面上的電機風(fēng)扇組組成。本文的側(cè)風(fēng)裝置模型參照了美國GE公司的發(fā)動機側(cè)風(fēng)裝置,如圖1所示。

    圖1 美國GE公司的發(fā)動機側(cè)風(fēng)裝置Fig.1 Aero-engine crosswind device of GE company in American

    對發(fā)動機側(cè)風(fēng)裝置進行合理簡化,省略無關(guān)部分,得到如圖2所示的發(fā)動機側(cè)風(fēng)裝置模型。裝置全長91 500 mm,寬30 000 mm,高18 890 mm。其中前端進氣部分長13 000 mm,試驗段殼體長67 840 mm。風(fēng)源部分由19 個風(fēng)筒及安裝在其內(nèi)部的電機風(fēng)扇組組成,風(fēng)筒長5 660 mm,風(fēng)筒前端進氣口直徑1 800 mm,尾部出氣口直徑2 160 mm,風(fēng)扇組的排列方式如圖3所示。風(fēng)扇模型是在沒有真實模擬風(fēng)扇的三維模型條件下的簡化,為無厚度平面,通過設(shè)置經(jīng)過該平面的壓力增量近似模擬風(fēng)扇的增壓效果。

    圖2 發(fā)動機側(cè)風(fēng)裝置模型Fig.2 Model of engine crosswind device

    圖3 風(fēng)扇排列方式及編號Fig.3 Fan arrangement and number

    2.3 網(wǎng)格劃分

    在側(cè)風(fēng)裝置前10 倍動力段距離和側(cè)風(fēng)裝置后30倍動力段距離處建立延伸段,以便能更好地觀察氣流流出下游出口之后的跡線,如圖4 所示。為提高側(cè)風(fēng)模擬計算精度,分區(qū)進行網(wǎng)格劃分,同時經(jīng)數(shù)值仿真對網(wǎng)格無關(guān)性進行驗證,網(wǎng)格如圖5 所示。風(fēng)源筒體及計算域網(wǎng)格適當加密,網(wǎng)格數(shù)量140萬。風(fēng)扇動力段網(wǎng)格數(shù)量50萬,如圖6所示。

    圖4 完整流場網(wǎng)格Fig.4 Complete flow field grid

    圖5 風(fēng)源筒體及計算域網(wǎng)格Fig.5 Fan barrel and computational domain grid

    圖6 風(fēng)扇動力段網(wǎng)格Fig.6 Fan power section grid

    2.4 邊界條件設(shè)置

    分別計算在中央1臺風(fēng)機、中間7臺風(fēng)機與全部19臺風(fēng)機工作(圖7)下,每臺風(fēng)扇進口分別提供200,500,750,1 000,1 500,2 000 Pa壓強時,側(cè)風(fēng)裝置下游出口處氣流的壓力、流量和速度。計算域邊界條件設(shè)置如圖8所示,紫色區(qū)域設(shè)置為壓力入口,黃色區(qū)域設(shè)置為壓力出口,粉色區(qū)域設(shè)置為無滑移、絕熱壁面。

    圖7 風(fēng)機工作情況Fig.7 Working condition of fan

    圖8 邊界條件設(shè)置Fig.8 Boundary condition setting

    運用ANSYS Workbench Fluent 18.2進行數(shù)值模擬,計算時的求解器選擇密度基求解器,流體設(shè)定為理想氣體,計算模型的速度特征選擇相對速度求解。計算過程中,采用的控制方程為三維雷諾平均Navier-Stokes方程,空間離散采用二階迎風(fēng)格式,湍流模型選用Realizablek-ε湍流模型。

    3 側(cè)風(fēng)物理量分析

    側(cè)風(fēng)裝置出口氣流流動會受到側(cè)風(fēng)裝置不同數(shù)量風(fēng)機工作的影響,下面從側(cè)風(fēng)裝置出口壓力、流量和速度三個方面分析側(cè)風(fēng)裝置出口氣流的流動情況。

    3.1 側(cè)風(fēng)裝置出口壓力變化規(guī)律

    采用數(shù)值模擬方法對側(cè)風(fēng)裝置出口處的壓力進行研究,對比分析不同數(shù)量風(fēng)機工作時的出口壓力變化規(guī)律。圖9為各個工況下出口壓力的分布云圖(從左到右分別對應(yīng)風(fēng)扇進口總壓200,500,750,1 000,1 500,2000 Pa)??梢?,1 臺、7 臺、19 臺風(fēng)機工作時,出口壓力均隨進口總壓改變,且最大壓力位置并不處于出口截面中心位置。1 臺風(fēng)機工作時,最大壓力出現(xiàn)在出口截面上部,且進口總壓越大越明顯;7臺風(fēng)機和19臺風(fēng)機工作時,高進口總壓的出口截面壓力分布比1 臺風(fēng)機工作時的更均勻,最大壓力出現(xiàn)在出口截面下部和上部的邊角。

    圖9 不同數(shù)量風(fēng)機工作時出口截面的壓力分布Fig.9 Pressure distribution of outlet section when different numbers of fans are working

    進一步對風(fēng)機工作時的跡線圖進行分析。進口總壓2 000 Pa,1臺、7臺、19臺風(fēng)機工作時側(cè)風(fēng)裝置的氣流跡線分別如圖10~圖12 所示。1 臺風(fēng)機工作時,氣流跡線在出口截面處小幅度偏向上方,導(dǎo)致氣流在出口前一較短距離內(nèi)向上方聚集,形成高壓區(qū),使得出口截面處高壓區(qū)域整體偏向上方。此時,在收斂段會產(chǎn)生回流區(qū),如圖10(b)所示。風(fēng)機所在風(fēng)筒為進氣狀態(tài),而其余未工作風(fēng)機所在風(fēng)筒為開口風(fēng)筒,氣流可以自由進出。出口截面之后會有1 000 mm左右的空間處于負壓狀態(tài),即該區(qū)域絕對壓強低于外界大氣壓強,造成除10號風(fēng)機所在風(fēng)筒外的其余風(fēng)筒處于從外界吸氣的狀態(tài)。但是由于該負壓相對較小(>-10 Pa),因此氣流在這些風(fēng)筒內(nèi)的速度并不是很高,進入收斂段后不能及時從下游出口排出,在負壓區(qū)形成回流區(qū)。

    圖10 1臺風(fēng)機工作時側(cè)風(fēng)裝置的跡線圖(進口總壓2 000 Pa)Fig.10 Trace diagram of working device of 1 fan with 2 000 Pa total inlet pressure

    圖11 7臺風(fēng)機工作時側(cè)風(fēng)裝置的跡線圖(進口總壓2 000 Pa)Fig.11 Trace diagram of working device of 7 fans with 2 000 Pa total inlet pressure

    圖12 19臺風(fēng)機工作時側(cè)風(fēng)裝置的跡線圖(進口總壓2 000 Pa)Fig.12 Trace diagram of working device of 19 fans with 2 000 Pa total inlet pressure

    7臺風(fēng)機工作時,在收斂段同樣會產(chǎn)生回流區(qū),但其形成機理與1 臺風(fēng)機工作時的形成機理不同。此時,工作風(fēng)機向內(nèi)吸氣會在風(fēng)筒截面后形成高壓區(qū),使被吸入空氣因為內(nèi)外壓差被強行從工作風(fēng)機所在風(fēng)筒擠出收斂段內(nèi)部空間,從而在內(nèi)部形成回流區(qū)。

    19臺風(fēng)機工作時,在收斂段內(nèi)部沒有明顯的回流區(qū)產(chǎn)生,流場均勻穩(wěn)定。

    3.2 側(cè)風(fēng)裝置出口流量變化規(guī)律

    上述從壓強方面解釋了1 臺風(fēng)機與7 臺風(fēng)機工作時收斂段內(nèi)部產(chǎn)生回流區(qū)的原因,本節(jié)通過比較進、出口流量變化規(guī)律對產(chǎn)生回流區(qū)的原因進行分析說明。不同數(shù)量風(fēng)機工作時的氣體流量見表1??梢钥闯?,1 臺風(fēng)機工作時,所有風(fēng)筒處于進氣狀態(tài),進氣流量等于出氣流量。7臺風(fēng)機工作時,只有工作風(fēng)機所在風(fēng)筒進氣,而氣流從非工作風(fēng)機所在風(fēng)筒與下游出口流出,使得進氣流量等于出氣流量。所以1臺風(fēng)機與7臺風(fēng)機工作時產(chǎn)生回流區(qū)的方式不同。

    表1 不同數(shù)量風(fēng)機工作時的氣體流量Table 1 Gas flow during operation of different number of fans

    圖13為進口總壓2 000 Pa時,1臺、7臺、19臺風(fēng)機工作時每臺風(fēng)機對進氣流量的貢獻曲線圖。從圖中可清晰看出不同數(shù)量風(fēng)機工作時每臺風(fēng)機輸送的流量:1臺風(fēng)機工作時,幾乎所有進口流量由10號工作風(fēng)機提供,其余18 臺未工作風(fēng)機也有少量進氣,但對總體進口流量影響不大,只是會導(dǎo)致側(cè)風(fēng)裝置收斂延長段內(nèi)部產(chǎn)生輕微的紊流;7臺風(fēng)機工作時,進口流量由5、6、9、10、11、14、15 號工作風(fēng)機提供,其余12 臺未工作風(fēng)機不僅沒有向裝置內(nèi)部提供流量,反而造成吸入流量損失,進而造成收斂延伸段內(nèi)部產(chǎn)生嚴重的回流區(qū),對出口流量影響較大;19 臺風(fēng)機工作時,由于全部風(fēng)機都在工作,所以每臺風(fēng)機承擔的流量近似,這與1臺風(fēng)機和7臺風(fēng)機工作時工作風(fēng)機與非工作風(fēng)機提供的流量有明顯不同。圖中也進一步說明,1臺風(fēng)機和7臺風(fēng)機工作時產(chǎn)生回流區(qū)的機理完全不同。

    圖13 進口氣壓2 000 Pa時各風(fēng)機流量貢獻曲線Fig.13 Contribution curve of fan flow at inlet pressure of 2 000 Pa

    對比三種數(shù)量風(fēng)機工作情況,1 臺風(fēng)機工作時效率最高,且未工作風(fēng)機也能提供一定的流量;7臺風(fēng)機工作時,流場相對穩(wěn)定,但未工作風(fēng)機會引起一部分流量損失;19 臺風(fēng)機工作時,流場穩(wěn)定均勻。因此,試驗過程中應(yīng)盡量選用全部19臺風(fēng)機工作。

    3.3 側(cè)風(fēng)裝置出口速度變化規(guī)律

    1 臺、7 臺、19 臺風(fēng)機工作時各進口壓力下的出口速度分布云圖如圖14所示??梢钥闯?,隨著進口總壓增大,出口速度也隨之增大,且其在單一條件中的分布也具有規(guī)律性,即越靠近壁面速度越低。1臺風(fēng)機工作時,出口截面的最大速度位于中心點右上方位置,這是受未工作風(fēng)機影響,使得流場紊亂所致。工作風(fēng)機數(shù)量增多,流場更加均勻。

    圖14 不同數(shù)量風(fēng)機工作時的出口速度變化Fig.14 Variation diagram of outlet speed when different numbers of fans are working

    以出口截面中心為原點建立坐標系,水平方向每間隔200 mm 取一個坐標點,原點左右取±80 mm的兩點,則出口截面處水平方向的速度分布如圖15所示??梢姡隹谒俣仍诳拷诿嫣帟E然下降,越靠近壁面處速度越低,且近壁面處速度梯度較大。這是因為緊貼壁面處有一層因壁面限制而導(dǎo)致脈動消失的層流薄層,其黏滯力使得流速急劇下降。

    圖15 不同工況下的出口速度分布曲線Fig.15 Outlet velocity distribution curve under different working conditions

    4 結(jié)論

    利用數(shù)值模擬方法,對某發(fā)動機露天試車臺側(cè)風(fēng)裝置不同工況進行了仿真計算分析,主要得出以下結(jié)論:

    (1) 相比全部19 臺風(fēng)機工作,中央1 臺風(fēng)機工作時的總壓損失大,且氣流流動紊亂;相比中間7臺風(fēng)機和全部19臺風(fēng)機工作,中央1臺風(fēng)機工作時的流量、速度受總壓的影響小。

    (2) 中央1 臺風(fēng)機單獨工作時效率最高,比19臺風(fēng)機一起工作時每臺風(fēng)機提供流量的能力都強,其周圍不工作的風(fēng)機受壓差作用會從外界吸氣,使得側(cè)風(fēng)裝置出口流量高于工作風(fēng)機吸入流量,但未工作風(fēng)機吸入的氣流流速低且紊亂。

    (3) 19臺風(fēng)機工作時,在壁面和中心區(qū)之間的區(qū)域流動較為穩(wěn)定;中央1臺風(fēng)機工作時,流速和壓力較高區(qū)域都偏離中心區(qū);中間7臺風(fēng)機工作時,流動相對穩(wěn)定,但流量會受未工作風(fēng)機影響而降低。試驗時,應(yīng)盡量避免中央1臺風(fēng)機和中間7臺風(fēng)機工作引起的試驗誤差。

    猜你喜歡
    風(fēng)筒總壓臺風(fēng)
    臺風(fēng)過韓
    總壓探針性能結(jié)構(gòu)敏感性分析
    可調(diào)式總壓耙設(shè)計及應(yīng)用
    離心式壓縮機異型風(fēng)筒法蘭車削簡圖的計算模型研究
    臺風(fēng)來了
    小讀者(2020年4期)2020-06-16 03:33:46
    亞聲速條件下總壓探針臨壁效應(yīng)的數(shù)值研究
    2 m超聲速風(fēng)洞流場變速壓控制方法研究
    臺風(fēng)愛搗亂
    機械通風(fēng)冷卻塔風(fēng)筒內(nèi)流場數(shù)值模擬研究
    高瓦斯長大隧道大直徑風(fēng)筒安裝施工技術(shù)
    精品国产一区二区久久| 国产人伦9x9x在线观看| 99久久99久久久精品蜜桃| 又紧又爽又黄一区二区| 热99久久久久精品小说推荐| 狠狠狠狠99中文字幕| 午夜精品久久久久久毛片777| 欧美一级毛片孕妇| 一本大道久久a久久精品| 亚洲伊人色综图| 中文字幕人妻丝袜一区二区| 亚洲天堂av无毛| 免费在线观看黄色视频的| 国产精品久久久久久人妻精品电影 | kizo精华| 国产精品久久久久久精品古装| 午夜福利影视在线免费观看| 9色porny在线观看| 高清av免费在线| 国产极品粉嫩免费观看在线| 精品国产乱码久久久久久小说| 波多野结衣一区麻豆| 大香蕉久久成人网| 国产熟女午夜一区二区三区| 脱女人内裤的视频| 少妇粗大呻吟视频| 国产欧美日韩精品亚洲av| 久久精品成人免费网站| 大片电影免费在线观看免费| 精品国产国语对白av| 欧美日韩黄片免| 亚洲美女黄片视频| 丝袜美腿诱惑在线| 色婷婷av一区二区三区视频| 亚洲专区字幕在线| 国产精品香港三级国产av潘金莲| 国产福利在线免费观看视频| 夫妻午夜视频| 亚洲精品一卡2卡三卡4卡5卡| 色在线成人网| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品 国内视频| 18禁裸乳无遮挡动漫免费视频| 久久影院123| 成年人午夜在线观看视频| 99久久人妻综合| 精品乱码久久久久久99久播| 欧美成狂野欧美在线观看| 亚洲色图 男人天堂 中文字幕| 欧美成狂野欧美在线观看| 亚洲第一av免费看| 天天躁夜夜躁狠狠躁躁| 免费日韩欧美在线观看| 日本vs欧美在线观看视频| 精品人妻熟女毛片av久久网站| 一区二区三区激情视频| 色视频在线一区二区三区| 久久99一区二区三区| 两个人免费观看高清视频| 亚洲专区国产一区二区| 一边摸一边抽搐一进一小说 | 91成人精品电影| 9191精品国产免费久久| 久久久水蜜桃国产精品网| 欧美人与性动交α欧美软件| 国产单亲对白刺激| 精品福利永久在线观看| 一区二区日韩欧美中文字幕| 久久国产精品影院| 最近最新中文字幕大全电影3 | 777米奇影视久久| 国产有黄有色有爽视频| 在线永久观看黄色视频| 好男人电影高清在线观看| 国产日韩欧美视频二区| 黄色毛片三级朝国网站| 老司机午夜福利在线观看视频 | 国产极品粉嫩免费观看在线| 国产成人精品久久二区二区91| 蜜桃在线观看..| videos熟女内射| 黄色视频在线播放观看不卡| 热re99久久国产66热| 国产精品亚洲av一区麻豆| 一边摸一边抽搐一进一小说 | 国产一区二区三区在线臀色熟女 | 国产亚洲精品久久久久5区| 啦啦啦免费观看视频1| 日本黄色视频三级网站网址 | 91九色精品人成在线观看| 两个人看的免费小视频| 最新在线观看一区二区三区| 成年人黄色毛片网站| 国产男女超爽视频在线观看| 亚洲成人国产一区在线观看| 国产不卡一卡二| 欧美亚洲日本最大视频资源| 99国产精品一区二区三区| 午夜免费鲁丝| 亚洲国产欧美在线一区| bbb黄色大片| 久久精品国产亚洲av高清一级| 亚洲精品国产精品久久久不卡| 日韩欧美国产一区二区入口| 另类精品久久| 老司机在亚洲福利影院| 欧美乱码精品一区二区三区| 高潮久久久久久久久久久不卡| 熟女少妇亚洲综合色aaa.| 亚洲人成电影免费在线| 亚洲性夜色夜夜综合| tube8黄色片| 久久久久久人人人人人| netflix在线观看网站| 久久久久精品国产欧美久久久| 欧美日韩福利视频一区二区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲久久久国产精品| 亚洲成人免费av在线播放| 国产精品98久久久久久宅男小说| 色尼玛亚洲综合影院| 十八禁人妻一区二区| 亚洲精品一卡2卡三卡4卡5卡| 9色porny在线观看| 亚洲人成电影免费在线| 亚洲精品自拍成人| 精品免费久久久久久久清纯 | aaaaa片日本免费| 国产成人影院久久av| 欧美日韩中文字幕国产精品一区二区三区 | 精品少妇久久久久久888优播| 久久久国产成人免费| 天天添夜夜摸| 国产欧美日韩一区二区三| 天堂俺去俺来也www色官网| 精品欧美一区二区三区在线| 中文亚洲av片在线观看爽 | 亚洲人成电影观看| 黄片小视频在线播放| 自线自在国产av| 天天躁日日躁夜夜躁夜夜| 国产单亲对白刺激| 国产色视频综合| 性色av乱码一区二区三区2| 十八禁网站免费在线| 国产真人三级小视频在线观看| 纯流量卡能插随身wifi吗| 国产麻豆69| 国产区一区二久久| 亚洲精华国产精华精| 99国产精品一区二区蜜桃av | 新久久久久国产一级毛片| 女同久久另类99精品国产91| 亚洲专区字幕在线| xxxhd国产人妻xxx| 人人澡人人妻人| 男女高潮啪啪啪动态图| 法律面前人人平等表现在哪些方面| 亚洲成a人片在线一区二区| 国产欧美日韩综合在线一区二区| 超色免费av| 大型av网站在线播放| 国产免费福利视频在线观看| 久久免费观看电影| 国产色视频综合| 久久精品国产亚洲av高清一级| 亚洲性夜色夜夜综合| av有码第一页| 亚洲精品美女久久av网站| 国产伦理片在线播放av一区| 国产精品国产av在线观看| kizo精华| 久久毛片免费看一区二区三区| 两个人免费观看高清视频| 欧美精品高潮呻吟av久久| 新久久久久国产一级毛片| 女人爽到高潮嗷嗷叫在线视频| 激情在线观看视频在线高清 | 亚洲精品久久午夜乱码| 午夜精品国产一区二区电影| 亚洲午夜理论影院| 亚洲欧美日韩高清在线视频 | 精品一品国产午夜福利视频| 成人国产一区最新在线观看| 三上悠亚av全集在线观看| 亚洲精品成人av观看孕妇| 国产精品国产高清国产av | 黑人欧美特级aaaaaa片| 亚洲综合色网址| 国产97色在线日韩免费| 中文字幕色久视频| 国产男女内射视频| 制服人妻中文乱码| 人人澡人人妻人| 日日夜夜操网爽| 色播在线永久视频| 在线观看免费高清a一片| 国产成人精品久久二区二区免费| 大片免费播放器 马上看| 亚洲va日本ⅴa欧美va伊人久久| 手机成人av网站| 日韩一区二区三区影片| 黄网站色视频无遮挡免费观看| 国产无遮挡羞羞视频在线观看| 狠狠精品人妻久久久久久综合| 国产精品影院久久| 精品亚洲成国产av| 亚洲国产看品久久| 一进一出好大好爽视频| 精品少妇黑人巨大在线播放| 国产成人精品在线电影| 国产有黄有色有爽视频| 国产精品 国内视频| 国产精品偷伦视频观看了| 午夜精品久久久久久毛片777| 日本撒尿小便嘘嘘汇集6| 熟女少妇亚洲综合色aaa.| 免费在线观看日本一区| 一区二区三区国产精品乱码| 无限看片的www在线观看| 97人妻天天添夜夜摸| 侵犯人妻中文字幕一二三四区| 999精品在线视频| 午夜免费成人在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 少妇被粗大的猛进出69影院| 欧美黑人欧美精品刺激| 日韩熟女老妇一区二区性免费视频| 午夜福利免费观看在线| a级毛片在线看网站| 操出白浆在线播放| 一本色道久久久久久精品综合| 男女高潮啪啪啪动态图| 久久国产亚洲av麻豆专区| 欧美精品啪啪一区二区三区| 动漫黄色视频在线观看| 久久久国产欧美日韩av| 精品乱码久久久久久99久播| 精品国产一区二区三区久久久樱花| 亚洲一区二区三区欧美精品| 久久香蕉激情| 涩涩av久久男人的天堂| 欧美av亚洲av综合av国产av| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人欧美| 精品一区二区三卡| 人人妻人人澡人人爽人人夜夜| 麻豆国产av国片精品| 色综合欧美亚洲国产小说| 精品卡一卡二卡四卡免费| 日韩一区二区三区影片| 亚洲国产看品久久| 国产免费现黄频在线看| 久久九九热精品免费| 久久热在线av| 欧美亚洲日本最大视频资源| 在线天堂中文资源库| 日韩一卡2卡3卡4卡2021年| 精品少妇一区二区三区视频日本电影| 最近最新中文字幕大全免费视频| 咕卡用的链子| 国产午夜精品久久久久久| 多毛熟女@视频| 三上悠亚av全集在线观看| 免费观看人在逋| 成人三级做爰电影| 免费看a级黄色片| 亚洲人成电影观看| 下体分泌物呈黄色| 91av网站免费观看| 黄色毛片三级朝国网站| 中文字幕人妻丝袜一区二区| 国产片内射在线| 国产亚洲精品第一综合不卡| 亚洲av国产av综合av卡| 欧美国产精品va在线观看不卡| 黄色视频不卡| 黄片播放在线免费| 日韩欧美一区二区三区在线观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 啦啦啦视频在线资源免费观看| 精品国产乱码久久久久久小说| √禁漫天堂资源中文www| 在线十欧美十亚洲十日本专区| 国产精品香港三级国产av潘金莲| 久久这里只有精品19| 久久精品亚洲熟妇少妇任你| 成年动漫av网址| 欧美精品人与动牲交sv欧美| 色综合欧美亚洲国产小说| 国产老妇伦熟女老妇高清| 丝瓜视频免费看黄片| 成人亚洲精品一区在线观看| 最黄视频免费看| 亚洲色图综合在线观看| 欧美日韩亚洲高清精品| 少妇被粗大的猛进出69影院| 18禁美女被吸乳视频| 色婷婷av一区二区三区视频| 日日夜夜操网爽| 女人被躁到高潮嗷嗷叫费观| 考比视频在线观看| 1024视频免费在线观看| 最新美女视频免费是黄的| 婷婷成人精品国产| 免费高清在线观看日韩| 午夜久久久在线观看| 十八禁高潮呻吟视频| 亚洲色图 男人天堂 中文字幕| 97人妻天天添夜夜摸| 国内毛片毛片毛片毛片毛片| 男男h啪啪无遮挡| 久久中文看片网| 777米奇影视久久| 乱人伦中国视频| 一区在线观看完整版| 9热在线视频观看99| 黄色视频,在线免费观看| 看免费av毛片| 日韩欧美三级三区| 制服人妻中文乱码| 国产91精品成人一区二区三区 | 18禁观看日本| 后天国语完整版免费观看| 又紧又爽又黄一区二区| 亚洲七黄色美女视频| 黑人操中国人逼视频| 捣出白浆h1v1| 淫妇啪啪啪对白视频| 久久中文看片网| 国产片内射在线| 久久亚洲精品不卡| av一本久久久久| 一区二区日韩欧美中文字幕| 亚洲国产欧美在线一区| 国产精品香港三级国产av潘金莲| 午夜激情久久久久久久| 午夜福利乱码中文字幕| 日韩大片免费观看网站| 国产日韩欧美亚洲二区| 女人精品久久久久毛片| 精品亚洲乱码少妇综合久久| 99久久精品国产亚洲精品| 少妇粗大呻吟视频| 人人妻人人澡人人看| 丰满饥渴人妻一区二区三| 久久久久久久久久久久大奶| 国产精品一区二区精品视频观看| 视频在线观看一区二区三区| 少妇 在线观看| av网站免费在线观看视频| 一本大道久久a久久精品| 日韩视频一区二区在线观看| 在线观看免费高清a一片| 捣出白浆h1v1| 丁香六月欧美| 国产一区二区三区在线臀色熟女 | 黄片播放在线免费| 天堂动漫精品| 淫妇啪啪啪对白视频| 欧美老熟妇乱子伦牲交| 91大片在线观看| av视频免费观看在线观看| 成人18禁在线播放| 国产不卡av网站在线观看| 午夜激情久久久久久久| av在线播放免费不卡| 自拍欧美九色日韩亚洲蝌蚪91| 午夜福利视频在线观看免费| 首页视频小说图片口味搜索| 天天影视国产精品| 自拍欧美九色日韩亚洲蝌蚪91| 欧美av亚洲av综合av国产av| 国产男女内射视频| 久久精品aⅴ一区二区三区四区| 亚洲欧洲日产国产| 日韩欧美一区视频在线观看| 一级片'在线观看视频| 无遮挡黄片免费观看| 少妇的丰满在线观看| videos熟女内射| 菩萨蛮人人尽说江南好唐韦庄| 在线亚洲精品国产二区图片欧美| 一进一出抽搐动态| av线在线观看网站| av片东京热男人的天堂| 在线观看人妻少妇| 精品欧美一区二区三区在线| 在线观看免费视频网站a站| 午夜久久久在线观看| 搡老岳熟女国产| 国产单亲对白刺激| 男女午夜视频在线观看| 成人国产av品久久久| 亚洲av成人不卡在线观看播放网| 免费看十八禁软件| 亚洲专区字幕在线| 狠狠婷婷综合久久久久久88av| 天天躁日日躁夜夜躁夜夜| 国产一区二区激情短视频| 久久午夜综合久久蜜桃| 精品乱码久久久久久99久播| 制服人妻中文乱码| 1024香蕉在线观看| 黄色丝袜av网址大全| kizo精华| 久久午夜综合久久蜜桃| 亚洲,欧美精品.| 免费av中文字幕在线| 亚洲欧美激情在线| 最新的欧美精品一区二区| 曰老女人黄片| 国产成人精品无人区| 欧美日韩国产mv在线观看视频| 五月开心婷婷网| 我的亚洲天堂| 欧美激情 高清一区二区三区| 国产精品免费大片| avwww免费| 高清在线国产一区| 热99久久久久精品小说推荐| 天堂8中文在线网| 国产精品国产av在线观看| 成人18禁高潮啪啪吃奶动态图| 精品一品国产午夜福利视频| 国产aⅴ精品一区二区三区波| 国产精品香港三级国产av潘金莲| 三级毛片av免费| 日韩免费av在线播放| 汤姆久久久久久久影院中文字幕| 美女视频免费永久观看网站| 国产欧美日韩一区二区精品| 天天躁夜夜躁狠狠躁躁| 热99re8久久精品国产| 午夜两性在线视频| 亚洲一区二区三区欧美精品| 亚洲色图av天堂| 久久久久国内视频| 国产视频一区二区在线看| 亚洲五月婷婷丁香| 久久久精品区二区三区| 久久精品成人免费网站| 国产亚洲av高清不卡| 国产成人精品无人区| 大香蕉久久成人网| 亚洲一区中文字幕在线| cao死你这个sao货| 国产精品1区2区在线观看. | 男女免费视频国产| 最新美女视频免费是黄的| 美女高潮喷水抽搐中文字幕| 两个人看的免费小视频| 女性生殖器流出的白浆| 亚洲国产欧美网| 极品教师在线免费播放| 一本综合久久免费| 天天添夜夜摸| 下体分泌物呈黄色| 久久亚洲真实| 丰满饥渴人妻一区二区三| 高清毛片免费观看视频网站 | 三级毛片av免费| 亚洲精品在线美女| 亚洲国产av影院在线观看| 亚洲少妇的诱惑av| av又黄又爽大尺度在线免费看| videos熟女内射| 自拍欧美九色日韩亚洲蝌蚪91| 十八禁网站网址无遮挡| 国产aⅴ精品一区二区三区波| 久久精品aⅴ一区二区三区四区| 最新美女视频免费是黄的| 久久中文看片网| 久久热在线av| 美女高潮到喷水免费观看| av国产精品久久久久影院| 99re6热这里在线精品视频| 50天的宝宝边吃奶边哭怎么回事| 天堂动漫精品| 午夜91福利影院| 黄色视频,在线免费观看| 亚洲国产欧美网| 亚洲精品中文字幕一二三四区 | 欧美黑人欧美精品刺激| 淫妇啪啪啪对白视频| 亚洲久久久国产精品| 亚洲国产欧美网| 90打野战视频偷拍视频| 激情在线观看视频在线高清 | 2018国产大陆天天弄谢| 国产精品香港三级国产av潘金莲| 9色porny在线观看| 在线天堂中文资源库| 一本综合久久免费| 成人免费观看视频高清| 一本综合久久免费| 少妇裸体淫交视频免费看高清 | 欧美人与性动交α欧美精品济南到| 激情视频va一区二区三区| 制服诱惑二区| 狠狠精品人妻久久久久久综合| 少妇的丰满在线观看| 午夜福利欧美成人| 天堂8中文在线网| 黄色视频不卡| 黑人欧美特级aaaaaa片| 国产精品自产拍在线观看55亚洲 | 日韩制服丝袜自拍偷拍| 亚洲精品国产一区二区精华液| 国产真人三级小视频在线观看| 精品第一国产精品| 美女高潮到喷水免费观看| 老司机亚洲免费影院| 国产精品偷伦视频观看了| 超碰97精品在线观看| 精品一区二区三区四区五区乱码| 嫩草影视91久久| 老司机午夜十八禁免费视频| 91成人精品电影| 久久久久精品国产欧美久久久| 91精品三级在线观看| 欧美日韩国产mv在线观看视频| 十八禁网站免费在线| 99精品在免费线老司机午夜| 精品久久久久久久毛片微露脸| 精品少妇内射三级| cao死你这个sao货| 午夜91福利影院| 一级,二级,三级黄色视频| 国产高清视频在线播放一区| 老司机午夜十八禁免费视频| 另类亚洲欧美激情| 黄色a级毛片大全视频| 欧美日韩亚洲综合一区二区三区_| 国产精品久久久久成人av| 亚洲av日韩精品久久久久久密| 午夜91福利影院| 天堂中文最新版在线下载| 少妇 在线观看| 青草久久国产| 新久久久久国产一级毛片| 国产激情久久老熟女| 国产在线观看jvid| 亚洲欧洲精品一区二区精品久久久| 国产高清激情床上av| 一级毛片精品| 欧美激情高清一区二区三区| 欧美性长视频在线观看| 大型av网站在线播放| 国产精品免费视频内射| 国产色视频综合| 久久久久久免费高清国产稀缺| 成人国产av品久久久| 亚洲性夜色夜夜综合| 亚洲自偷自拍图片 自拍| 免费看十八禁软件| 免费在线观看影片大全网站| 一本色道久久久久久精品综合| 国产精品久久久av美女十八| 精品国内亚洲2022精品成人 | 一个人免费在线观看的高清视频| 大片免费播放器 马上看| 亚洲精品国产精品久久久不卡| 高清欧美精品videossex| 纵有疾风起免费观看全集完整版| 欧美精品啪啪一区二区三区| 国产伦人伦偷精品视频| 成年版毛片免费区| 国产国语露脸激情在线看| 国产亚洲精品第一综合不卡| 黑人欧美特级aaaaaa片| 免费高清在线观看日韩| 一边摸一边做爽爽视频免费| 亚洲人成电影免费在线| 国产精品香港三级国产av潘金莲| 久久人人爽av亚洲精品天堂| 日本a在线网址| 亚洲精品成人av观看孕妇| 亚洲成a人片在线一区二区| 美女午夜性视频免费| 啦啦啦视频在线资源免费观看| 国产有黄有色有爽视频| 黄色 视频免费看| 中文欧美无线码| 悠悠久久av| 欧美中文综合在线视频| 国产一区二区三区综合在线观看| 国产精品1区2区在线观看. | 国产黄色免费在线视频| 男女午夜视频在线观看| 我的亚洲天堂| 少妇裸体淫交视频免费看高清 | 亚洲色图 男人天堂 中文字幕| 99久久人妻综合| av天堂久久9| 19禁男女啪啪无遮挡网站| 97人妻天天添夜夜摸| 中文欧美无线码| 天天影视国产精品| 91大片在线观看| kizo精华| 美女福利国产在线| 青草久久国产| 乱人伦中国视频| 国产在线精品亚洲第一网站| 黄频高清免费视频| 黑人巨大精品欧美一区二区mp4| 日本vs欧美在线观看视频| av视频免费观看在线观看| 香蕉久久夜色| 亚洲一区二区三区欧美精品| 亚洲欧美色中文字幕在线| 国产在线精品亚洲第一网站| 啪啪无遮挡十八禁网站| 每晚都被弄得嗷嗷叫到高潮|