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

    再生冷卻結(jié)構(gòu)非定常換熱計算研究

    2023-05-30 05:05:28張若凌張磊王想義蔣勁顧磊
    實驗流體力學(xué) 2023年2期
    關(guān)鍵詞:發(fā)動機

    張若凌,張磊,王想義,蔣勁,顧磊

    中國空氣動力研究與發(fā)展中心 高超聲速沖壓發(fā)動機技術(shù)重點實驗室,綿陽 621000

    0 引 言

    碳氫燃料超燃沖壓發(fā)動機工作時,燃氣溫度很高,即使采用復(fù)合材料也無法承受燃燒室內(nèi)巨大的熱載荷。在發(fā)動機中,一般采用燃油再生冷卻方法解決熱防護問題[1-3]。采用燃油再生冷卻既可以利用燃料的吸熱性質(zhì)降低發(fā)動機壁面溫度,又可以減小系統(tǒng)重量。受燃料熱沉能力限制和積碳影響,目前國際上普遍認為在飛行馬赫數(shù)低于8 時可采用碳氫燃料[4-6]進行再生冷卻。

    中國空氣動力研究與發(fā)展中心(CARDC)從2003年開始進行超燃沖壓發(fā)動機再生冷卻傳熱的計算和試驗研究,基于國內(nèi)某航空煤油的組分測量建立了冷卻劑的物性模型,結(jié)合大規(guī)模并行CFD 流場計算工具AHL3D[7-9],建立了定常換熱分析工具和初步的冷卻通道設(shè)計方法,進行了冷卻劑分別為水和航空煤油的發(fā)動機試驗研究[10-13]。

    目前針對超燃沖壓發(fā)動機再生冷卻結(jié)構(gòu)非定常換熱的研究未見報道。進行非定常換熱分析,可以獲得發(fā)動機的熱平衡特性,對發(fā)動機再生冷卻結(jié)構(gòu)設(shè)計和試驗研究十分重要,可為熱結(jié)構(gòu)風(fēng)洞建設(shè)提供依據(jù)。

    本文首先建立發(fā)動機再生冷卻結(jié)構(gòu)的非定常換熱計算模型;然后,通過對比發(fā)動機(冷卻劑分別為水和航空煤油)計算與試驗數(shù)據(jù),驗證非定常換熱計算模型,分析非正規(guī)階段溫升時間;最后,計算分析發(fā)動機初始壁面溫度對達到熱平衡時間的影響。

    1 非定常換熱計算模型

    在壁面溫度趨近于平衡值的過程中,存在2 個階段:在第一階段里,溫度分布主要受初始溫度分布影響,即非正規(guī)階段;在第二階段里,不同時刻的溫度分布主要取決于邊界條件和物性,即正規(guī)階段[14]。若超燃沖壓發(fā)動機的壁面較薄,則非正規(guī)階段溫升的時間遠小于達到熱平衡的時間,下文的分析基于發(fā)動機壁面的溫度分布處于正規(guī)階段的假設(shè)。

    發(fā)動機壁面由圖1 所示的再生冷卻面板組成。假設(shè)發(fā)動機固壁內(nèi)在沿著冷卻劑流動的方向上沒有傳熱。一個通道在流向某位置上的截面如圖2 所示(y 軸沿冷卻通道的高度方向),燃氣垂直紙面向內(nèi)流動,冷卻劑垂直紙面向外流動,冷卻通道寬Lc、高Hc,肋片厚度D,一個冷卻通道的寬度為Lc+D,襯層厚度b1,涂層(TBC)厚度b2。

    圖1 再生冷卻面板Fig. 1 Regenerative cooling panel

    圖2 再生冷卻通道截面Fig. 2 Regenerative channel section

    燃氣側(cè)壁面熱流密度q1和襯層冷卻劑側(cè)熱流密度q2的計算公式[12]如下:

    式中:h1為燃氣與壁面的換熱系數(shù),T01為燃氣恢復(fù)總溫,T 為襯層的平均溫度,?T 為襯層兩側(cè)溫差的一半,?T1為涂層兩側(cè)溫差,h2為冷卻劑與襯層壁面的換熱系數(shù),T2為冷卻劑的平均溫度。

    發(fā)動機壁面選用不銹鋼或高溫合金材料,其導(dǎo)熱系數(shù)比冷卻劑大很多,熱量經(jīng)冷卻通道四周傳向冷卻劑。燃氣側(cè)壁面熱流密度q1可以分為4 部分:一部分熱流密度使冷卻劑正下方的襯層和涂層升溫,一部分熱流密度q2通過襯層傳向冷卻劑,一部分熱流密度q4傳給肋基,還有一部分熱流密度q5使肋基正下方的襯層和涂層升溫。假設(shè)冷卻面板外側(cè)絕熱,把肋片上端半寬為Lc/2 的外側(cè)金屬層部分也看作厚度為D/2 的肋片,則整個肋片的高度為Hc+Lc/2。

    假設(shè)整個襯層的溫度分布相同,通過肋基的熱流密度q4可以分成2 部分:一部分使整個肋片升溫,一部分是通過肋片傳向冷卻劑的熱流密度q3。q3計算公式如下:

    式中:h3為肋片與冷卻劑的平均換熱系數(shù),Tfin為肋片平均溫度。

    由于金屬的熱阻遠小于流體,肋基處熱流密度比冷卻劑側(cè)熱流密度大,即q4>q2。當發(fā)動機壁面溫度達到平衡分布時,肋基處熱流密度可以通過肋片效率η定義:

    式中:k1為襯層熱導(dǎo)率,下標“eq”表示平衡狀態(tài)。

    為計算非定常階段的溫升,令每一時刻使襯層、涂層升溫和通過襯層傳向冷卻劑的總熱流密度為φq1(φ為熱流比例系數(shù)),并假設(shè)在每一時刻均滿足:

    由式(6)、(7)可推導(dǎo)出:

    襯層和肋片升溫過程的控制方程分別為:

    式中:ρ1、c1分別為襯層的密度和比熱;ρ2、c2分別為涂層的密度和比熱,t 為時間。由圖2、式(10)、(11)可知,單位時間的溫升等于流入的熱流密度與流出的熱流密度之差除以單位面積單位時間溫升所需的熱量。式(10)可用來確定q5,同時根據(jù)式(9)和(10)可以確定q4。式(11)可以改寫為:

    考察襯層處情形,認為在每一時刻襯層和涂層的溫度分布處于正規(guī)階段,假設(shè):

    式中:k2為涂層熱導(dǎo)率。將式(1)、(2)代入式(13),得到:

    將式(1)、(2)、(8)代入式(10),并用式(14)替換?T、?T1,得到:

    其中,C 和λ表達式分別為:

    式(15)就是襯層溫升過程的控制方程,初值條件為襯層平均溫度T(t = 0) =T0。若近似認為物性和換熱系數(shù)為常數(shù),其解為:

    式中:Teq為平衡溫度, (Teq?T0)e?λt表征非定常過程,即襯層溫度向平衡溫度趨近的過程。需要指出的是,h1隨燃氣和壁面溫度變化不大,h2隨冷卻劑溫度變化較大。

    定義襯層非定常溫升時間τ為襯層平均溫度達到0.99 倍平衡溫度的時間,則:

    可以看出,τ除了與材料物性、結(jié)構(gòu)和冷卻劑的換熱性質(zhì)相關(guān)外,還與襯層初始溫度T0相關(guān),T0與Teq越接近(相對差別須大于1%),τ越小。

    肋片處的溫升過程根據(jù)式(12)求解,本文不做進一步推導(dǎo)。在計算中冷卻劑流動采用一維方法處理,沿著流動方向步進,就可以得到沿冷卻通道壁面的溫度分布。換熱系數(shù)h1采用Eckert 參考焓方法[15]計算,考慮進口效應(yīng)的換熱系數(shù)h2采用修正Bartz方程[16-17]計算。

    2 非定常換熱計算模型驗證

    2.1 發(fā)動機熱平衡時間分析

    CARDC 于2004—2005 年在電弧加熱器設(shè)備(圖3)上開展了水冷卻和航空煤油冷卻直連式超燃沖壓發(fā)動機(無燃燒)換熱試驗。發(fā)動機采用不銹鋼材料(1Cr18Ni9Ti)制成,燃氣側(cè)壁面沒有涂層。發(fā)動機氣體內(nèi)流道入口高32 mm、寬100 mm,冷卻通道寬3 mm、高3 mm,肋片寬3 mm,襯層厚3 mm。設(shè)備來流總溫1100 K,來流總壓1.6 MPa,空氣質(zhì)量流量3.2 kg/s,燃燒室入口馬赫數(shù)為2,冷卻劑入口溫度297 K、壓力7 MPa。

    圖3 再生冷卻發(fā)動機測試設(shè)備Fig. 3 Regeneratively cooled scramjet test facility

    圖4 和5 分別是水冷卻和航空煤油冷卻試驗條件下的典型外壁溫度歷程曲線??梢钥闯?,水冷卻試驗中模型的外壁溫度基本達到了穩(wěn)定狀態(tài),非定常溫升時間約23 s;航空煤油冷卻試驗中模型的外壁溫度未能達到穩(wěn)態(tài)分布,有效試驗時間約28.5 s。

    圖4 水冷卻試驗中外壁溫度歷程Fig. 4 Histories of the outer wall temperature in water cooling test

    圖5 航空煤油冷卻試驗中外壁溫度歷程Fig. 5 Histories of the outer wall temperature in kerosene cooling test

    采用非定常換熱計算模型,在相同工況下,對試驗中的超燃沖壓發(fā)動機結(jié)構(gòu)開展水冷和航空煤油冷卻換熱特性計算。單條冷卻通道的水和航空煤油質(zhì)量流量分別為3.5 和3.0 g/s。計算得到一條冷卻通道的冷卻劑側(cè)壁面不同點的溫升曲線分別如圖6和7 所示。由圖可知,水冷時發(fā)動機結(jié)構(gòu)達到熱平衡的時間約為20 s,航空煤油冷卻時達到熱平衡的時間約為45 s。

    圖6 水冷卻試驗條件下壁溫計算Fig. 6 Calculated wall temperature under water cooling condition

    圖7 航空煤油冷卻試驗條件下壁溫計算Fig. 7 Calculated wall temperature under kerosene cooling condition

    需要指出的是,在圖4、6 和7 中,確定模型外壁溫度基本達到穩(wěn)態(tài)的時間的方法具有一定的主觀性。更加嚴謹?shù)姆椒ㄒ蠼o出溫度達到穩(wěn)態(tài)的一致性判據(jù),例如在穩(wěn)態(tài)溫度未知時,采用溫度對時間導(dǎo)數(shù)為該導(dǎo)數(shù)最大值的5%或10%所需的時間,或者在穩(wěn)態(tài)溫度值已知時采用溫度達到0.99 倍平衡溫度的時間。本文意在說明相同的試驗參數(shù)下,航空煤油冷卻試驗中模型外壁溫度基本達到穩(wěn)態(tài)的時間明顯比水冷卻試驗更長,因此本文采用的方法也是可行的。

    從式(16)、(17)和(19)可以看出,采用不同冷卻劑得到發(fā)動機達到熱平衡的時間差別較大,主要原因是冷卻劑換熱系數(shù)不同。在所討論的條件下,水的換熱系數(shù)比航空煤油大,即水的冷卻性能比航空煤油好,使發(fā)動機達到熱平衡的時間更短。

    對比圖4 和6 可知,在水冷條件下,計算和試驗得到的發(fā)動機達到熱平衡時間差別約為3 s(在計算中假定肋片效率η在溫升過程中保持不變并忽略相鄰單元沿冷卻劑流向的傳熱,這些是產(chǎn)生誤差的因素),二者吻合較好,對本文建立的再生冷卻結(jié)構(gòu)非定常換熱計算模型進行了初步驗證。

    2.2 非正規(guī)階段溫升時間分析

    研究非正規(guī)階段溫升時間時,可采用平板模型近似替代帶冷卻通道的燃燒室壁板。該階段冷卻劑吸熱較少,簡化計算中可不考慮。一側(cè)加熱、一側(cè)絕熱的平板溫度分布的精確解為[14]:

    式中:Tp0為平板初始溫度,qp為加熱的熱流密度,ρp、cp、bp、kp分別為平板的密度、比熱、厚度和熱傳導(dǎo)系數(shù),熱擴散系數(shù)α=kp/(ρpcp),x = 0 為絕熱側(cè)壁面位置,x=bp為給定熱流密度側(cè)壁面位置。

    式(20)右端第四項表征非正規(guī)階段的溫度分布,隨著時間增長呈指數(shù)衰減,當(即傅立葉數(shù)Fo)較大時,可忽略此項。選取3.1 節(jié)不銹鋼材料的物性參數(shù),取ρp=8000 kg/m3、kp= 16 W/(m·K)(300 K 對應(yīng)的數(shù)據(jù))、bp= 3 mm、cp= 500 J/(K·kg)。當指數(shù)為–3 時(e–3= 5.0%),F(xiàn)o = 0.304,非正規(guī)階段的溫升時間為0.68 s。

    由圖4 和5 可知,水冷卻發(fā)動機非正規(guī)階段溫升時間約為2.5 s,航空煤油冷卻時該時間約為1 s,比式(20)計算得到的時間長。這是因為式(20)適用于一側(cè)絕熱平板,實際上再生冷卻超燃沖壓發(fā)動機的壁板不能簡單假設(shè)為平板,且冷卻劑還要從襯層吸熱。由于溫度較低時航空煤油的換熱系數(shù)比水小,因此航空煤油冷卻發(fā)動機非正規(guī)階段溫升時間較短。

    2.3 壁面預(yù)加熱對熱平衡影響分析

    在建設(shè)超燃沖壓發(fā)動機熱結(jié)構(gòu)風(fēng)洞時,需要掌握發(fā)動機達到熱平衡的時間,采用傳熱性能良好的材料和通道構(gòu)型可以縮短達到熱平衡時間。由式(19)可知,在發(fā)動機試車前對發(fā)動機壁面進行預(yù)加熱處理(即提高T0初始值),可縮短發(fā)動機達到熱平衡的時間,從而節(jié)約熱結(jié)構(gòu)風(fēng)洞的運行時間。

    采用本文的非定常換熱計算模型,計算發(fā)動機一條冷卻通道壁面初始溫度分別為300 和800 K 時的非定常階段溫升過程,得到不同時刻冷卻劑側(cè)壁面沿程溫度分布如圖8 和9 所示。

    圖8 初始溫度300 K 時壁面溫度沿通道方向的分布Fig. 8 Distributions of wall temperature along tube under initial wall temperature 300 K

    圖9 初始溫度800 K 時壁面溫度沿通道方向的分布Fig. 9 Distributions of wall temperature along tube under initial wall temperature 800 K

    計算采用的發(fā)動機壁面材料為不銹鋼(1Cr18Ni9Ti),燃氣側(cè)壁面有0.1 mm 的ZrO2涂層。來流條件對應(yīng)飛行馬赫數(shù)6,發(fā)動機冷卻和燃燒油氣比均為1.1(燃氣總溫最高約2700 K)。冷卻劑為航空煤油,入口壓力7 MPa,從發(fā)動機出口向入口方向流動,每根通道流量1 g/s。由圖8 和9 可知,初始壁面溫度分別為300 和800 K 時,發(fā)動機達到熱平衡時間分別約為23 和13 s。由此可知,對發(fā)動機壁面進行預(yù)加熱可以縮短達到熱平衡的時間。需要指出的是,此處計算得到的達到熱平衡的時間,未考慮冷卻劑流量的影響。

    進行非定常階段溫升時間的初步估算對于一些實際應(yīng)用意義重大。在包含冷卻流動的傳熱試驗中,傳熱過程是十分復(fù)雜的。壁溫升高特性主要取決于:1)固體的傳熱過程,非常接近線性;2)在冷卻通道里的流動中,存在復(fù)雜的非線性特征,即冷卻劑物性參數(shù)發(fā)生很大變化以及實際的對流邊界層存在非線性特征;3)燃燒室熱流密度分布的不均勻特性,主要是總溫和激波等干擾帶來的加熱條件分布不均。

    若要采用某個簡單公式來進行非定常階段溫升時間的初步估算,必須采用一些簡單假設(shè)。通過反復(fù)驗算發(fā)現(xiàn),對于典型的再生冷卻超燃沖壓發(fā)動機而言,冷卻劑物性取300 K 對應(yīng)的數(shù)據(jù),燃氣恢復(fù)溫度取1650 K(對應(yīng)飛行馬赫數(shù)6 ),換熱系數(shù)和壁面物性取壁溫900 K 對應(yīng)的數(shù)據(jù),利用式(19)進行非定常階段溫升時間初步估算是合適的。以圖8 和9 研究的再生冷卻超燃沖壓發(fā)動機結(jié)構(gòu)(1/λ=5.34 s,C/λ=892.2 K)為例:當初始壁面溫度300 K時,得到非定常階段溫升時間22.39 s,與圖8 結(jié)果一致;初始壁面溫度800 K 時得到非定常階段溫升時間12.46 s,與圖9 結(jié)果一致。

    3 結(jié) 論

    本文建立了超燃沖壓發(fā)動機再生冷卻結(jié)構(gòu)的一種非定常換熱計算模型,進行了發(fā)動機換熱特性分析,得到以下結(jié)論:

    1)在所討論的水冷卻條件下,計算和試驗得到的發(fā)動機達到熱平衡時間分別為23 s 和20 s,二者較為吻合,初步驗證了本文建立的非定常換熱計算模型。

    2)給出了一種發(fā)動機非正規(guī)階段溫升時間的計算方法,且溫升時間的計算與測量結(jié)果吻合較好。

    3)計算結(jié)果表明,試驗前對發(fā)動機壁面進行預(yù)加熱可以縮短發(fā)動機達到熱平衡的時間,從而節(jié)約熱結(jié)構(gòu)風(fēng)洞的運行時間。

    4)本文建立的再生冷卻結(jié)構(gòu)非定常換熱計算模型具有一定的適用性,能夠應(yīng)用于超燃沖壓發(fā)動機的傳熱特性分析。

    猜你喜歡
    發(fā)動機
    元征X-431實測:奔馳發(fā)動機編程
    2015款寶馬525Li行駛中發(fā)動機熄火
    2012年奔馳S600發(fā)動機故障燈偶爾點亮
    發(fā)動機空中起動包線擴展試飛組織與實施
    RD-180超級火箭發(fā)動機的興衰
    太空探索(2016年8期)2016-07-10 09:21:58
    奔馳E200車發(fā)動機故障燈常亮
    奔馳E260冷車時發(fā)動機抖動
    新一代MTU2000發(fā)動機系列
    2013年車用發(fā)動機排放控制回顧(下)
    2013年車用發(fā)動機排放控制回顧(上)
    男的添女的下面高潮视频| 国产午夜福利久久久久久| 欧美日韩在线观看h| 97在线视频观看| 亚洲av免费高清在线观看| 麻豆乱淫一区二区| 超碰97精品在线观看| 亚洲精品一区蜜桃| 久久这里有精品视频免费| 人人妻人人澡人人爽人人夜夜 | 亚洲av二区三区四区| 国产精品无大码| АⅤ资源中文在线天堂| 日本一二三区视频观看| 亚洲人成网站在线播| 免费人成在线观看视频色| .国产精品久久| 亚洲精品自拍成人| 国产黄a三级三级三级人| 嫩草影院入口| 欧美另类亚洲清纯唯美| 日韩欧美在线乱码| 小说图片视频综合网站| 亚洲欧美日韩东京热| 欧美人与善性xxx| 亚洲电影在线观看av| 国产精品无大码| 97在线视频观看| 日本欧美国产在线视频| 精品一区二区三区视频在线| 亚洲综合精品二区| 2021少妇久久久久久久久久久| 国产精品久久视频播放| 青青草视频在线视频观看| 免费观看人在逋| 视频中文字幕在线观看| 美女内射精品一级片tv| 日韩亚洲欧美综合| 熟女人妻精品中文字幕| 水蜜桃什么品种好| 成年免费大片在线观看| 亚洲国产精品国产精品| 欧美日韩国产亚洲二区| 免费av不卡在线播放| 日本五十路高清| 国产私拍福利视频在线观看| 听说在线观看完整版免费高清| 看片在线看免费视频| 中文资源天堂在线| 欧美日韩一区二区视频在线观看视频在线 | 能在线免费观看的黄片| 一级av片app| 亚洲三级黄色毛片| 国产黄片美女视频| 熟女电影av网| 久久精品国产自在天天线| 日本午夜av视频| 噜噜噜噜噜久久久久久91| 岛国在线免费视频观看| 亚洲一区高清亚洲精品| 麻豆一二三区av精品| 网址你懂的国产日韩在线| 成人毛片a级毛片在线播放| 亚洲自拍偷在线| 欧美又色又爽又黄视频| 国产国拍精品亚洲av在线观看| av线在线观看网站| 99热全是精品| 六月丁香七月| 亚洲最大成人手机在线| 水蜜桃什么品种好| 99热全是精品| 亚洲图色成人| 精品久久久久久久人妻蜜臀av| 特级一级黄色大片| 三级男女做爰猛烈吃奶摸视频| 欧美精品一区二区大全| 黄片wwwwww| 亚洲内射少妇av| 欧美成人a在线观看| 免费在线观看成人毛片| 亚洲精品一区蜜桃| av在线亚洲专区| 国产午夜福利久久久久久| 亚洲精品成人久久久久久| 天堂√8在线中文| 国产成人a∨麻豆精品| 久99久视频精品免费| av黄色大香蕉| 亚洲性久久影院| 天天躁夜夜躁狠狠久久av| 亚洲怡红院男人天堂| 日日摸夜夜添夜夜添av毛片| 久久久国产成人精品二区| 国产人妻一区二区三区在| 欧美成人免费av一区二区三区| 可以在线观看毛片的网站| 日本熟妇午夜| 国产乱人视频| 亚洲怡红院男人天堂| 久久综合国产亚洲精品| 1024手机看黄色片| 老司机福利观看| 亚洲人成网站在线观看播放| 美女脱内裤让男人舔精品视频| 国产av在哪里看| 亚洲av成人av| 亚洲婷婷狠狠爱综合网| 成人鲁丝片一二三区免费| 国产成人a区在线观看| 我的老师免费观看完整版| 精品久久国产蜜桃| 中文字幕制服av| 色网站视频免费| 国产精品一区二区三区四区久久| 亚洲国产精品成人久久小说| 欧美变态另类bdsm刘玥| 中文欧美无线码| 精品免费久久久久久久清纯| 亚洲不卡免费看| 久久精品91蜜桃| 男女边吃奶边做爰视频| 欧美另类亚洲清纯唯美| 大话2 男鬼变身卡| 最近视频中文字幕2019在线8| 嘟嘟电影网在线观看| .国产精品久久| 亚洲熟妇中文字幕五十中出| 成人二区视频| 午夜福利在线观看免费完整高清在| 成年女人永久免费观看视频| 女的被弄到高潮叫床怎么办| 最近中文字幕2019免费版| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产三级中文精品| 国产黄片视频在线免费观看| 中文字幕久久专区| 少妇人妻一区二区三区视频| 大香蕉97超碰在线| 国产激情偷乱视频一区二区| 亚洲欧美成人综合另类久久久 | 男人舔奶头视频| 亚洲av一区综合| 欧美成人精品欧美一级黄| 亚洲经典国产精华液单| 欧美精品国产亚洲| 亚洲人成网站在线播| 天堂中文最新版在线下载 | 国产精品国产三级专区第一集| 中文字幕熟女人妻在线| 18禁在线无遮挡免费观看视频| 狂野欧美白嫩少妇大欣赏| 乱码一卡2卡4卡精品| 纵有疾风起免费观看全集完整版 | 亚洲国产最新在线播放| 中文资源天堂在线| 最新中文字幕久久久久| 欧美色视频一区免费| 97超碰精品成人国产| 简卡轻食公司| 国产免费男女视频| 亚洲国产欧美在线一区| 日韩一区二区三区影片| 日韩,欧美,国产一区二区三区 | 亚洲人与动物交配视频| 国产乱人偷精品视频| 国产老妇女一区| 国产一区亚洲一区在线观看| 国产亚洲精品久久久com| 欧美xxxx性猛交bbbb| 久久精品熟女亚洲av麻豆精品 | 干丝袜人妻中文字幕| 老师上课跳d突然被开到最大视频| 中文字幕人妻熟人妻熟丝袜美| 18+在线观看网站| 中文资源天堂在线| 天堂中文最新版在线下载 | 欧美色视频一区免费| 老司机福利观看| 国产成人a区在线观看| 精品少妇黑人巨大在线播放 | 久久6这里有精品| 欧美色视频一区免费| 秋霞在线观看毛片| 久久久久久久久久久丰满| 成人综合一区亚洲| 看免费成人av毛片| 日韩,欧美,国产一区二区三区 | 美女cb高潮喷水在线观看| 免费黄色在线免费观看| 非洲黑人性xxxx精品又粗又长| 麻豆精品久久久久久蜜桃| 99热精品在线国产| 亚洲欧美日韩东京热| 美女大奶头视频| 一级毛片久久久久久久久女| 狠狠狠狠99中文字幕| 国内精品宾馆在线| 欧美极品一区二区三区四区| 日本熟妇午夜| 22中文网久久字幕| 国产不卡一卡二| 有码 亚洲区| 日本一二三区视频观看| 欧美成人午夜免费资源| 日韩国内少妇激情av| 日韩欧美三级三区| 亚洲久久久久久中文字幕| 男的添女的下面高潮视频| 91久久精品国产一区二区成人| 青春草视频在线免费观看| 亚洲国产精品sss在线观看| 午夜福利视频1000在线观看| 伦理电影大哥的女人| 国产亚洲91精品色在线| 精品久久久久久久久久久久久| 又粗又爽又猛毛片免费看| 亚洲婷婷狠狠爱综合网| 日本熟妇午夜| 久久久久久大精品| 成人综合一区亚洲| 免费播放大片免费观看视频在线观看 | 欧美成人精品欧美一级黄| 我的老师免费观看完整版| 91在线精品国自产拍蜜月| 国产淫语在线视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 美女xxoo啪啪120秒动态图| 日韩欧美国产在线观看| av免费观看日本| 亚洲欧美中文字幕日韩二区| 亚洲国产高清在线一区二区三| 舔av片在线| 极品教师在线视频| 黄色日韩在线| 久久综合国产亚洲精品| 国产色爽女视频免费观看| 欧美区成人在线视频| av在线天堂中文字幕| 亚洲av成人精品一区久久| 婷婷色综合大香蕉| 在线免费观看不下载黄p国产| 在线观看66精品国产| 中文字幕熟女人妻在线| 国产精品野战在线观看| 爱豆传媒免费全集在线观看| 色尼玛亚洲综合影院| 日本五十路高清| 男人和女人高潮做爰伦理| 国内精品一区二区在线观看| 99九九线精品视频在线观看视频| 午夜视频国产福利| 亚洲av免费在线观看| 欧美另类亚洲清纯唯美| 亚洲精品日韩av片在线观看| 久久99热6这里只有精品| 大话2 男鬼变身卡| 男人的好看免费观看在线视频| 亚洲最大成人中文| 偷拍熟女少妇极品色| 夜夜爽夜夜爽视频| 七月丁香在线播放| 99热这里只有精品一区| www日本黄色视频网| 乱系列少妇在线播放| 在线观看66精品国产| 午夜亚洲福利在线播放| av在线亚洲专区| 长腿黑丝高跟| 国产一区有黄有色的免费视频 | 变态另类丝袜制服| 国产一区有黄有色的免费视频 | 国产成人福利小说| 国产高清三级在线| 久久亚洲精品不卡| 欧美成人精品欧美一级黄| 国产三级中文精品| 熟女人妻精品中文字幕| 午夜福利网站1000一区二区三区| 毛片女人毛片| 熟女电影av网| 亚洲国产最新在线播放| 夫妻性生交免费视频一级片| 亚洲欧美日韩卡通动漫| 少妇熟女aⅴ在线视频| or卡值多少钱| 国产成人午夜福利电影在线观看| 床上黄色一级片| 国产亚洲av片在线观看秒播厂 | 久久人人爽人人爽人人片va| 日本一本二区三区精品| 美女黄网站色视频| 亚洲美女搞黄在线观看| 精品久久久久久久久亚洲| 99热全是精品| 国产爱豆传媒在线观看| 中文乱码字字幕精品一区二区三区 | 国产欧美日韩精品一区二区| 久久婷婷人人爽人人干人人爱| 久久99热6这里只有精品| 欧美一区二区精品小视频在线| 久久久久九九精品影院| 欧美区成人在线视频| 国产免费视频播放在线视频 | 黄片wwwwww| 色视频www国产| 91午夜精品亚洲一区二区三区| 国产在视频线精品| 精品国产露脸久久av麻豆 | 成人性生交大片免费视频hd| 亚洲国产高清在线一区二区三| 性插视频无遮挡在线免费观看| 免费看光身美女| 一级二级三级毛片免费看| 天堂av国产一区二区熟女人妻| 黑人高潮一二区| 久久久久久久久久久丰满| 国产一区二区在线观看日韩| 亚洲av不卡在线观看| 欧美日韩在线观看h| 99热全是精品| 国产亚洲5aaaaa淫片| 亚洲欧美一区二区三区国产| 亚洲激情五月婷婷啪啪| 亚洲四区av| 在线a可以看的网站| 亚洲国产精品成人综合色| 日韩国内少妇激情av| 99九九线精品视频在线观看视频| 欧美日韩国产亚洲二区| 国产精品精品国产色婷婷| 日韩一区二区视频免费看| 亚洲国产欧洲综合997久久,| 久久久久久久亚洲中文字幕| 久久久久久伊人网av| 国产久久久一区二区三区| 亚洲成人中文字幕在线播放| 久热久热在线精品观看| 日本色播在线视频| 中文字幕亚洲精品专区| 精品国产三级普通话版| av线在线观看网站| 国产精品一及| 观看免费一级毛片| 亚洲av福利一区| 最近最新中文字幕大全电影3| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品福利在线免费观看| 男人狂女人下面高潮的视频| 精华霜和精华液先用哪个| 欧美丝袜亚洲另类| 婷婷色麻豆天堂久久 | 在线观看一区二区三区| 亚洲精品乱码久久久v下载方式| 男人舔女人下体高潮全视频| 色尼玛亚洲综合影院| 日韩一本色道免费dvd| 一二三四中文在线观看免费高清| 欧美高清性xxxxhd video| 伦精品一区二区三区| 亚洲av中文av极速乱| 久久6这里有精品| 尾随美女入室| 在线免费观看的www视频| 国产69精品久久久久777片| 亚洲中文字幕日韩| 偷拍熟女少妇极品色| 成人美女网站在线观看视频| 美女大奶头视频| 99在线视频只有这里精品首页| 国产探花极品一区二区| 天堂√8在线中文| 免费黄网站久久成人精品| 中文字幕av成人在线电影| 久久午夜福利片| 欧美日本视频| 亚洲伊人久久精品综合 | 欧美变态另类bdsm刘玥| 色综合站精品国产| 中文字幕av在线有码专区| 男女国产视频网站| 在线播放无遮挡| 亚洲最大成人av| 中国美白少妇内射xxxbb| 一级二级三级毛片免费看| 国产成人一区二区在线| 久久精品久久精品一区二区三区| 看十八女毛片水多多多| 一区二区三区乱码不卡18| 熟女电影av网| 亚洲欧美中文字幕日韩二区| 2022亚洲国产成人精品| 国产精品伦人一区二区| 成人国产麻豆网| 99久久成人亚洲精品观看| 亚洲av.av天堂| 又黄又爽又刺激的免费视频.| 亚洲18禁久久av| 免费人成在线观看视频色| 国产精品日韩av在线免费观看| 久久久欧美国产精品| 美女xxoo啪啪120秒动态图| 国产探花极品一区二区| 亚洲婷婷狠狠爱综合网| 国产成人精品久久久久久| 亚洲欧美中文字幕日韩二区| 免费看a级黄色片| 免费看日本二区| 少妇被粗大猛烈的视频| 天堂影院成人在线观看| 国产高清视频在线观看网站| 亚洲人成网站在线观看播放| 校园人妻丝袜中文字幕| av福利片在线观看| 欧美色视频一区免费| 天天躁夜夜躁狠狠久久av| 日韩一区二区三区影片| 欧美bdsm另类| 久久久a久久爽久久v久久| 国产一区亚洲一区在线观看| 国产午夜精品一二区理论片| 免费观看人在逋| 国产毛片a区久久久久| 精品久久久久久久末码| 天堂中文最新版在线下载 | 久久久久久大精品| 国产成人freesex在线| 天天躁日日操中文字幕| 身体一侧抽搐| 国产视频内射| 搡女人真爽免费视频火全软件| 国产亚洲av片在线观看秒播厂 | 精品免费久久久久久久清纯| av福利片在线观看| 麻豆成人午夜福利视频| 国产精品国产三级专区第一集| 国产精品伦人一区二区| 久久久久久伊人网av| 大又大粗又爽又黄少妇毛片口| 久久久久久大精品| www.av在线官网国产| 丰满少妇做爰视频| 日韩欧美 国产精品| 人妻系列 视频| 精品久久久久久久久亚洲| 久久精品国产鲁丝片午夜精品| 亚洲精品色激情综合| 国产在线男女| 亚洲av熟女| 国产av码专区亚洲av| 成人性生交大片免费视频hd| 国产真实伦视频高清在线观看| 亚洲国产精品久久男人天堂| 长腿黑丝高跟| 欧美激情国产日韩精品一区| 最后的刺客免费高清国语| 亚洲欧美精品自产自拍| 最近的中文字幕免费完整| 成人三级黄色视频| 国产乱人偷精品视频| 久久人人爽人人爽人人片va| 麻豆成人av视频| 最近中文字幕高清免费大全6| 成人亚洲欧美一区二区av| 精品免费久久久久久久清纯| av在线亚洲专区| 免费搜索国产男女视频| 麻豆久久精品国产亚洲av| 中文字幕精品亚洲无线码一区| av国产免费在线观看| 中文在线观看免费www的网站| 精品久久久久久久久av| 搡老妇女老女人老熟妇| 夜夜看夜夜爽夜夜摸| 国产精品一区www在线观看| 国产精品无大码| 亚洲精品一区蜜桃| 国产精品三级大全| 一个人看视频在线观看www免费| 欧美xxxx黑人xx丫x性爽| 欧美精品一区二区大全| 亚洲欧美精品综合久久99| 日韩大片免费观看网站 | 久久久成人免费电影| 成人一区二区视频在线观看| 最近中文字幕高清免费大全6| 九九在线视频观看精品| 一级黄片播放器| 3wmmmm亚洲av在线观看| 久久久精品大字幕| 一级毛片久久久久久久久女| 熟女电影av网| 日韩一区二区三区影片| 亚洲伊人久久精品综合 | 久久久久久久久久黄片| 尾随美女入室| 亚洲va在线va天堂va国产| 午夜日本视频在线| 99热这里只有精品一区| 国产黄a三级三级三级人| 两性午夜刺激爽爽歪歪视频在线观看| 国产精华一区二区三区| 日本五十路高清| 国产精品永久免费网站| 男女边吃奶边做爰视频| 国产精品一区二区三区四区免费观看| 国产精品熟女久久久久浪| 汤姆久久久久久久影院中文字幕 | 一级毛片aaaaaa免费看小| 视频中文字幕在线观看| 国产午夜福利久久久久久| 日本黄色片子视频| 亚洲综合色惰| 欧美日本视频| 久久久久九九精品影院| 2021天堂中文幕一二区在线观| 亚洲不卡免费看| 有码 亚洲区| 国产毛片a区久久久久| 国产老妇伦熟女老妇高清| 国产 一区精品| 99热这里只有是精品50| 一个人免费在线观看电影| 人体艺术视频欧美日本| 国产精品一二三区在线看| 亚洲伊人久久精品综合 | 长腿黑丝高跟| 欧美激情久久久久久爽电影| 日韩亚洲欧美综合| 黄片wwwwww| 一级毛片aaaaaa免费看小| 亚洲欧美日韩卡通动漫| 成人毛片a级毛片在线播放| 中文字幕亚洲精品专区| 精品久久久久久久久av| 久久久久久久久久久免费av| 国产精品美女特级片免费视频播放器| 亚洲欧美成人综合另类久久久 | 免费搜索国产男女视频| 日韩一本色道免费dvd| 国产久久久一区二区三区| 欧美性猛交黑人性爽| 亚洲成色77777| 非洲黑人性xxxx精品又粗又长| 色综合亚洲欧美另类图片| 亚洲精华国产精华液的使用体验| 18禁裸乳无遮挡免费网站照片| 亚洲精品久久久久久婷婷小说 | 黄片wwwwww| 成人欧美大片| 亚洲五月天丁香| 日日啪夜夜撸| av在线老鸭窝| 激情 狠狠 欧美| 国产成人一区二区在线| 狠狠狠狠99中文字幕| 亚洲经典国产精华液单| 亚洲高清免费不卡视频| 亚洲中文字幕日韩| 国产老妇伦熟女老妇高清| 中文天堂在线官网| 最近的中文字幕免费完整| 欧美日本亚洲视频在线播放| 午夜精品国产一区二区电影 | 欧美丝袜亚洲另类| 日韩欧美精品v在线| 国产久久久一区二区三区| 美女黄网站色视频| 精品国内亚洲2022精品成人| 国产视频首页在线观看| 少妇熟女aⅴ在线视频| 午夜福利在线观看吧| 国内精品美女久久久久久| 亚洲人与动物交配视频| 精品人妻视频免费看| 亚洲不卡免费看| 亚洲国产欧美在线一区| 波多野结衣巨乳人妻| 日本免费a在线| 黄色日韩在线| 国产精品熟女久久久久浪| 日韩一本色道免费dvd| 又爽又黄a免费视频| 最近视频中文字幕2019在线8| 日本-黄色视频高清免费观看| 女的被弄到高潮叫床怎么办| 丰满乱子伦码专区| 欧美日韩一区二区视频在线观看视频在线 | 午夜a级毛片| 人人妻人人看人人澡| 一级毛片我不卡| 91av网一区二区| 成人午夜精彩视频在线观看| 性色avwww在线观看| 日本午夜av视频| 秋霞伦理黄片| 久久久久九九精品影院| 婷婷色综合大香蕉| 国产成人午夜福利电影在线观看| 国产亚洲精品久久久com| 国产不卡一卡二| 亚洲欧美成人精品一区二区| 在现免费观看毛片| 日本wwww免费看| 欧美+日韩+精品| 国产成人aa在线观看| 丰满人妻一区二区三区视频av| 九草在线视频观看| 秋霞在线观看毛片| 中文精品一卡2卡3卡4更新| 亚洲av二区三区四区| 国产69精品久久久久777片| 麻豆国产97在线/欧美| 小说图片视频综合网站| 欧美成人精品欧美一级黄|