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

    小容積流量下汽輪機(jī)調(diào)節(jié)級動葉三維流場與級后工質(zhì)溫度的研究*

    2016-12-25 02:06:16王涵北
    化工機(jī)械 2016年4期
    關(guān)鍵詞:動葉工質(zhì)容積

    李 勇 王涵北

    (東北電力大學(xué)能源與動力工程學(xué)院)

    小容積流量下汽輪機(jī)調(diào)節(jié)級動葉三維流場與級后工質(zhì)溫度的研究*

    李 勇**王涵北

    (東北電力大學(xué)能源與動力工程學(xué)院)

    以某300MW汽輪機(jī)調(diào)節(jié)級為研究對象,介紹了比例系數(shù)法和數(shù)值模擬法的基本原理。對不同相對容積流量比和不同工況下的調(diào)節(jié)級動葉三維流場進(jìn)行了分析。將數(shù)值模擬方法得到的調(diào)節(jié)級級后工質(zhì)溫度與比例系數(shù)法、速度三角形法的計(jì)算結(jié)果進(jìn)行對比,得到不同計(jì)算方法間的相對誤差,并對不同相對容積流量下級后工質(zhì)溫度的計(jì)算方法進(jìn)行了探討。

    汽輪機(jī) 動葉三維流場 級后工質(zhì)溫度 相對容積流量比

    隨著風(fēng)能、太陽能等清潔能源在電力能源系統(tǒng)中所占比例的不斷增加,要求參與調(diào)峰的大容量火電機(jī)組必須通過降低負(fù)荷,甚至是啟停調(diào)峰的方式來適應(yīng)新能源發(fā)電的介入。因此,火電機(jī)組的快速啟停工作受到國內(nèi)外學(xué)者的廣泛關(guān)注。在確定火電機(jī)組啟停速度和壽命損耗中,汽輪機(jī)調(diào)節(jié)級后的溫度變化率對優(yōu)化火電機(jī)組啟停具有重要影響。為確定主蒸汽流量較小時的級后工質(zhì)參數(shù),文獻(xiàn)[1~4]提出當(dāng)汽輪機(jī)負(fù)荷低于80%額定負(fù)荷時,采用比例系數(shù)法進(jìn)行參數(shù)計(jì)算,然而這種估算方法并未考慮汽輪機(jī)沖轉(zhuǎn)流量變化帶來的影響,必然引起調(diào)節(jié)級級后溫度變化率產(chǎn)生較大誤差。

    隨著CFD技術(shù)的興起與日漸成熟,數(shù)值模擬方法在汽輪機(jī)調(diào)節(jié)級和小容積流量研究中獲得了廣泛應(yīng)用。文獻(xiàn)[5~7]通過模擬調(diào)節(jié)級三維流場,分析了工質(zhì)經(jīng)過調(diào)節(jié)級堵塞區(qū)域時通道中靜壓的變化情況。謝金偉等通過建立完整的調(diào)節(jié)閥至調(diào)節(jié)級段模型,對不同工況下該段的三維流動與壓力損失特性進(jìn)行了研究[8]。文獻(xiàn)[9~12]對汽輪機(jī)調(diào)節(jié)級氣動性能和氣流激振突變方面進(jìn)行了研究,提出了優(yōu)化調(diào)節(jié)級結(jié)構(gòu)的可行方案。文獻(xiàn)[13,14]通過對低負(fù)荷運(yùn)行的汽輪機(jī)末幾級進(jìn)行數(shù)值模擬,得到當(dāng)容積流量小到一定程度時末級動葉會做負(fù)功的結(jié)論。張仲彬等對多級流動進(jìn)行了數(shù)值分析,結(jié)果表明,小容積流量工況時上級流場的分離流動會導(dǎo)致下級流場內(nèi)分離區(qū)域增多[15]。

    目前國內(nèi)外學(xué)者普遍關(guān)注的是小容積流量工況下汽輪機(jī)的末級流動狀態(tài),而調(diào)節(jié)級小容積流量方面的研究卻鮮有報道。在此,筆者對不同相對容積流量比和不同工況下的調(diào)節(jié)級動葉三維流場進(jìn)行分析,并將采用數(shù)值模擬方法得到的調(diào)節(jié)級級后工質(zhì)溫度與比例系數(shù)法、速度三角形法計(jì)算結(jié)果進(jìn)行對比,得到不同計(jì)算方法間的相對誤差,以指導(dǎo)不同容積流量下級后工質(zhì)溫度的計(jì)算。

    1 比例系數(shù)法

    比例系數(shù)法中首先根據(jù)熱力特性說明書查出汽輪機(jī)額定負(fù)荷下的主蒸汽溫度、壓力,再進(jìn)行計(jì)算:

    (1)

    (2)

    式中p0、p01——額定工況和變工況下的主蒸汽壓力,MPa;

    pi、pi1——額定工況和變工況下某一級特征面的蒸汽壓力,MPa;

    t0、t01——額定工況和變工況下的主蒸汽溫度,℃;

    ti、ti1——額定工況和變工況下某一級特征面的蒸汽溫度,℃。

    由式(1)、(2)可知,比例系數(shù)法存在兩方面問題。首先,對于汽輪機(jī)沖轉(zhuǎn)階段而言,主蒸汽參數(shù)基本不變,因此該方法計(jì)算出的調(diào)節(jié)級級后工質(zhì)參數(shù)也不隨時間改變。而事實(shí)上,汽輪機(jī)沖轉(zhuǎn)時需要調(diào)節(jié)閥全開來保證汽缸和轉(zhuǎn)子受熱均勻,因此多數(shù)電廠采用高壓主汽閥控制汽輪機(jī)轉(zhuǎn)速以便進(jìn)行低速檢查和中速暖機(jī),這就使主汽閥開度在沖轉(zhuǎn)階段隨時間變化,進(jìn)而引起調(diào)節(jié)級入口、出口工質(zhì)參數(shù)隨時間改變。其次,式(1)、(2)中沒有體現(xiàn)沖轉(zhuǎn)流量改變對調(diào)節(jié)級級后工質(zhì)參數(shù)的影響。

    2 物理模型和數(shù)值方法

    2.1 物理模型

    筆者對某300MW汽輪機(jī)調(diào)節(jié)級葉片及其流域進(jìn)行建模,葉片幾何尺寸見表1。

    表1 某300MW汽輪機(jī)調(diào)節(jié)級葉片幾何尺寸

    采用Gambit軟件建模時,由表1可知,調(diào)節(jié)級動、靜葉片數(shù)不同,同一工況下單位靜葉流道與單位動葉流道的質(zhì)量流量不同,因此不能采用單葉片模型進(jìn)行研究,而建立整周模型所需網(wǎng)格數(shù)量大,容易導(dǎo)致計(jì)算機(jī)內(nèi)存超限并延長計(jì)算時間,經(jīng)觀察發(fā)現(xiàn),5個靜葉與4個動葉構(gòu)成的流域具有相同質(zhì)量流量,遂建立如圖1所示的物理模型。筆者選用的調(diào)節(jié)級模型部分進(jìn)汽度高達(dá)0.96,建模時可忽略其影響。

    在調(diào)節(jié)閥全開的情況下,過熱蒸汽通過汽輪機(jī)蒸汽室進(jìn)入調(diào)節(jié)級靜葉,經(jīng)噴嘴加速、動葉做功后流入調(diào)節(jié)級汽室進(jìn)行充分混合。為保證建立的

    圖1 某300MW汽輪機(jī)調(diào)節(jié)級物理模型

    模型與實(shí)際情況相似,同時考慮到網(wǎng)格劃分的需要,設(shè)置動葉入口截面距動葉前緣點(diǎn)0.2mm,動葉出口截面距動葉尾緣點(diǎn)0.2mm,并在動葉后建立軸向距離等于真實(shí)長度的調(diào)節(jié)級汽室模型,在保證動葉出口流體充分發(fā)展的同時使流動更加真實(shí)。汽室部分在動葉生成流域的基礎(chǔ)上切割而成,不但提高了動葉與汽室連接處的順滑性,也為劃分高質(zhì)量網(wǎng)格提供了保證。葉柵通道網(wǎng)格(結(jié)構(gòu)化網(wǎng)格)劃分情況如圖2所示。

    圖2 葉柵通道網(wǎng)格劃分情況

    2.2 網(wǎng)格無關(guān)性討論

    為兼顧計(jì)算速度和準(zhǔn)確度,筆者對模型進(jìn)行了網(wǎng)格無關(guān)性討論。在保證計(jì)算方法、邊界條件不變的前提下,分別對網(wǎng)格數(shù)為184萬、292萬、346萬、422萬的模型進(jìn)行數(shù)值模擬。結(jié)果表明,當(dāng)網(wǎng)格數(shù)大于346萬后,隨著網(wǎng)格數(shù)的增加模擬結(jié)果幾乎不變,故綜合考慮后取346萬網(wǎng)格的方案。

    2.3 數(shù)值方法

    采用Fluent軟件進(jìn)行計(jì)算,選擇適用于參數(shù)研究的標(biāo)準(zhǔn)k-ε湍流模型,應(yīng)用Simplec算法和二階迎風(fēng)離散格式,采用旋轉(zhuǎn)周期性邊界條件,動葉流域及其固體壁面均繞x軸以給定轉(zhuǎn)速旋轉(zhuǎn),靜葉、汽室流域及其固體壁面保持靜止。靜葉、動葉、調(diào)節(jié)級汽室三者交界面采用滑移網(wǎng)格技術(shù)進(jìn)行數(shù)據(jù)交換。入口邊界條件為壓力入口,出口邊界條件為質(zhì)量流量出口(mass-flow inlet中的outward normals),過熱蒸汽采用相較于理想氣體更接近真實(shí)水蒸氣性質(zhì)的Redlich-Kwong氣體模型,能量方程收斂精度為10-6,速度方程收斂精度為10-3。

    3 計(jì)算結(jié)果分析

    3.1 動葉三維流場的分析

    由文獻(xiàn)[15]可知,當(dāng)本級靜葉入口工質(zhì)均勻時,容積流量減小很難引起靜葉流域發(fā)生明顯的流動分離,故主要分析動葉三維流場的情況。為了排除汽輪機(jī)轉(zhuǎn)速對流場的影響,將動葉及其流域轉(zhuǎn)速均設(shè)為3kr/min。相對容積流量比k為:

    k=G1v1/(Gv)

    (3)

    式中G、G1——閥門全開和變工況下調(diào)節(jié)級動葉出口蒸汽質(zhì)量流量,kg/s;

    v、v1——閥門全開和變工況下調(diào)節(jié)級動葉出口蒸汽比容,m3/kg。

    以閥門全開工況為基準(zhǔn),圖3為不同k下同一動葉流道50%葉高截面處的流線圖,圖4為k為0.290時同一工況下不同葉高截面處的流線圖。

    圖3 不同k下同一動葉流道50%葉高截面處的流線圖

    圖4 同一工況下不同葉高截面處的流線圖

    由圖3可知,隨著k值的減小,調(diào)節(jié)級動葉壓力面前緣點(diǎn)附近會出現(xiàn)流動分離,分離區(qū)與主流之間有一條明顯的分離線。出現(xiàn)流動分離的原因是容積流量較小時調(diào)節(jié)級靜葉出口汽流絕對速度較小,由于動葉高速旋轉(zhuǎn),導(dǎo)致汽流以較大負(fù)攻角進(jìn)入動葉流域并被吸力面分為兩支,其中一支在旋轉(zhuǎn)葉片的作用下沿本流道繼續(xù)向下游運(yùn)動,另一分支流入相鄰流道并在動葉壓力面附近產(chǎn)生漩渦。當(dāng)k>0.026時,壓力面漩渦軸向影響范圍隨k值的減小而增大,這是因?yàn)榇藭r負(fù)攻角較大,而且容積流量相對較高,工質(zhì)經(jīng)堵塞段加速后具有較大的軸向分速度,使再附點(diǎn)逐漸遠(yuǎn)離分離點(diǎn);當(dāng)k<0.026時,動葉出口吸力面附近出現(xiàn)了流動分離且分離區(qū)隨k值的減小而增大,此分離區(qū)的出現(xiàn)使動葉內(nèi)工質(zhì)橫向流動明顯加劇,加之此時容積流量過小,汽流軸向分速度較低,使壓力面回流區(qū)軸向影響范圍稍有減小。

    由圖4可知,同一工況下不同葉高截面處的流動狀態(tài)存在差異,壓力面流動分離現(xiàn)象首先在動葉根部前緣點(diǎn)附近產(chǎn)生,且當(dāng)k=0.290時回流區(qū)范圍隨葉高的增加而減小直至消失。為了進(jìn)一步說明壓力面發(fā)生流動分離時的動葉三維流場情況,引入圖5所示的不同k值時動葉壓力面和下端壁極限流線圖。

    a. k=0.290

    b. k=0.330

    圖5a中k=0.290時的工況與文獻(xiàn)[16]類似,此時在動葉壓力面與下端壁附近均發(fā)生了流動分離。k=0.290時動葉壓力面前緣點(diǎn)附近沿葉高方向形成一條分離線,靠近分離線的下游流體在40%左右軸向弦長內(nèi)形成了大面積回流區(qū)且漩渦占據(jù)了整個葉高,導(dǎo)致近壁區(qū)域出現(xiàn)強(qiáng)烈的回流和徑向流動。40%軸向弦長及其后區(qū)域徑向流動逐步減弱,但整體流線呈斜下方趨勢。此外,k=0.290時下端壁附近也出現(xiàn)了明顯的流動分離,分離區(qū)內(nèi)反向流體沿斜上方運(yùn)動并向分離線靠攏。對比發(fā)現(xiàn),k=0.330時壓力面與下端壁附近均未出現(xiàn)回流,流動效果較好。綜上所述,k=0.290為調(diào)節(jié)級動葉流域出現(xiàn)流動分離現(xiàn)象的起始點(diǎn)。

    3.2 級后工質(zhì)溫度的分析

    速度三角形法和比例系數(shù)法都是計(jì)算級后工

    質(zhì)參數(shù)的傳統(tǒng)方法,由于筆者主要介紹數(shù)值模擬方法,因此只將新方法與傳統(tǒng)方法所得結(jié)果進(jìn)行比較。假設(shè)此時主汽閥和調(diào)節(jié)閥全開,則3種方法計(jì)算結(jié)果之間的相對誤差如圖6所示。

    圖6 不同級后工質(zhì)溫度計(jì)算方法間

    由圖6可知,在k>0.026范圍內(nèi)數(shù)值模擬法與速度三角形法計(jì)算結(jié)果的相對誤差小于2%,且基本保持不變,這是因?yàn)檎{(diào)節(jié)級徑高比較大,流場沿半徑方向雖有變化但變化程度較小,平均直徑處汽流參數(shù)仍適用于變工況的計(jì)算,雖然k≤0.290時動葉壓力面附近已出現(xiàn)流動分離,但此時仍可用速度三角形法估算級后工質(zhì)溫度;當(dāng)k≤0.026時,由于動葉出口吸力面附近出現(xiàn)了流動分離,使平均直徑處汽流出汽角相比于設(shè)計(jì)工況發(fā)生了較大變化,導(dǎo)致兩種方法計(jì)算結(jié)果的相對誤差呈大幅上升趨勢,此時速度三角形法得到的級后工質(zhì)溫度存在較大偏差,應(yīng)采用數(shù)值模擬方法進(jìn)行計(jì)算。比例系數(shù)法計(jì)算結(jié)果與數(shù)值模擬結(jié)果的相對誤差在k=0.330時已大于2%,且隨k值的減小而增大,在k<0.110范圍內(nèi)增幅更為明顯,可見比例系數(shù)法所得結(jié)果在容積流量較小時并不理想。綜上所述,為得到準(zhǔn)確的調(diào)節(jié)級級后工質(zhì)溫度,當(dāng)k>0.026時可采用速度三角形法進(jìn)行計(jì)算,當(dāng)k≤0.026時應(yīng)采用數(shù)值模擬方法。

    4 結(jié)論

    4.1對動葉三維流場的分析表明,當(dāng)k=0.290時,由于入口汽流負(fù)攻角較大,下端壁和壓力面附近首次出現(xiàn)流動分離。當(dāng)0.026

    4.2對數(shù)值模擬結(jié)果與速度三角形法計(jì)算結(jié)果的分析表明,當(dāng)k>0.026時,兩者的相對誤差小于2%,雖然在0.026

    4.3對數(shù)值模擬結(jié)果與比例系數(shù)法計(jì)算結(jié)果的分析表明,當(dāng)k≤0.330時,兩者相對誤差大于2%,可見容積流量較小時,采用比例系數(shù)法所得結(jié)果存在較大誤差,相比之下,數(shù)值模擬結(jié)果準(zhǔn)確度更高。

    [1] 陳堅(jiān)紅,阮偉,任浩仁,等.冷態(tài)啟動溫升率對汽輪機(jī)轉(zhuǎn)子壽命影響的研究[J].動力工程,1997,17(6):11~15.

    [2] 袁鵬飛,盛德仁,陳堅(jiān)紅,等.汽輪機(jī)轉(zhuǎn)子熱疲勞壽命損耗監(jiān)測面的有限元分析[J].電站系統(tǒng)工程,2004,20(2):43~44.

    [3] 付忠廣,杜志鋒,王夢嬌.660MW超超臨界汽輪機(jī)轉(zhuǎn)子熱應(yīng)力的研究[J].現(xiàn)代電力,2012,29(1):60~64.

    [4] 袁利軍.300MW機(jī)組中壓缸啟動熱應(yīng)力計(jì)算與啟動優(yōu)化[D].北京:華北電力大學(xué),2005.

    [5] 訾宏達(dá),方韌,徐星仲,等.汽輪機(jī)部分進(jìn)汽調(diào)節(jié)級三維流場的非定常數(shù)值模擬[J].動力工程,2007,27(5):721~724.

    [6]Lampart P,Szymaniak M,Rzadkowski R.Unsteady Load of Partial Admission Control Stage Rotor of a Lager Power Steam Turbine[C].ASME Turbo Expo 2004:Power for Land,Sea and Air.Glasgow:American Society of Mechanical Engineers, 2004:237~246.

    [7] Hills N J.Whole Turbine CFD Modelling[C].ASME Turbo Expo 2007:Power for Land,Sea and Air.Montreal:American Society of Mechanical Engineers,2007:817~824.

    [8] 謝金偉,王新軍,廖高良,等.汽輪機(jī)調(diào)節(jié)閥-調(diào)節(jié)級段三維流動與壓力損失數(shù)值研究[J].中國電機(jī)工程學(xué)報,2014,34(20): 3376~3383.

    [9] Hushmandi N B,Hu J,Fridh J,et al.Numerical Study of Unsteady Flow Phenomena in a Partial Admission Axial Steam Yurbine[C].ASME Turbo Expo 2008:Power for Land,Sea and Air.Berlin:American Society of Mechanical Engineers,2008: 713~722.

    [10] He L.Computation of Unsteady Flow through Steam Turbine Blade Row at Partial Adimission[J].Proceeding of the Institution of Mechanical Engineers Part A:Journal of Power and Energy,1997,211(3):197~205.

    [11] 陽虹,楊建道,李軍,等.汽輪機(jī)調(diào)節(jié)級氣動性能分析和結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J].熱力透平,2011,40(1):5~10.

    [12] 李建蘭,郭緒宏,余辰,等.汽輪機(jī)調(diào)節(jié)級汽流激振突變模型研究[J].中國電機(jī)工程學(xué)報,2013,33(11):39~46.

    [13] 楊建道.空冷50MW機(jī)組低壓缸變工況三維流場數(shù)值分析[J].熱力透平,2007,36(1):53~56.

    [14] Senoo S,Segawa K,Hamatake H,et al.Computation for Unsteady Compressible Flows in a Multi-Stage Steam Turbine with Steam Properties at Low Load Operations[J].Journal of Engineering for Gas Turbines and Power,2011,133(10): 711~722.

    [15] 張仲彬,田艷靜,曹麗華,等.小容積流量工況下汽輪機(jī)級內(nèi)流場的數(shù)值分析[J].化工機(jī)械,2013,40(1):94~97.

    [16] 朱光宇,俞茂錚,戴麗萍,等.汽輪機(jī)調(diào)節(jié)級動葉柵大負(fù)沖角工況下的三維分離流動研究[J].動力工程學(xué)報,2004,24(2): 189~194.

    Studyof3DFlowFieldsinGoverningStageandSteamTemperatureafterStageofSteamTurbineatSmallVolumetricFlow

    LI Yong, WANG Han-bei

    (SchoolofEnergyandPowerEngineering,NortheastDianliUniversity,Jilin132012,China)

    Taking the governing stage of a 300MW steam turbine as the object of study, the basic principle of proportional coefficient method and numerical simulation method was described; and the 3D flow fields in moving blades at governing stage with different volume flow ratios and at various working conditions were analyzed. Through comparing the steam temperature after governing stage and from the numerical simulation method with the results from the proportional coefficient method and velocity triangle method, the relative error between various calculation methods can be reached. These different calculation methods for steam temperature at different volumetric flows were discussed.

    steam turbine, 3D flow field in moving blade, steam temperature after stage, relative volume flow ratio

    *國家自然科學(xué)基金項(xiàng)目(51376041)。

    **李 勇,男,1964年1月生,教授。吉林省吉林市,132012。

    TQ053

    A

    0254-6094(2016)04-0489-06

    2015-09-21)

    猜你喜歡
    動葉工質(zhì)容積
    海洋溫差能發(fā)電熱力循環(huán)系統(tǒng)的工質(zhì)優(yōu)選
    怎樣求醬油瓶的容積
    M701DA燃機(jī)壓氣機(jī)動葉鎖鍵的配合方式研究
    采用R1234ze(E)/R245fa的非共沸混合工質(zhì)有機(jī)朗肯循環(huán)系統(tǒng)實(shí)驗(yàn)研究
    發(fā)電廠汽輪機(jī)振動異常增大的原因分析
    采用二元非共沸工質(zhì)的有機(jī)朗肯循環(huán)熱力學(xué)分析
    巧求容積
    試析引風(fēng)機(jī)動葉故障原因及處理措施
    若干低GWP 純工質(zhì)在空調(diào)系統(tǒng)上的應(yīng)用分析
    截斷的自適應(yīng)容積粒子濾波器
    亚洲精品中文字幕在线视频 | 美女内射精品一级片tv| 亚洲自偷自拍三级| 亚洲伊人久久精品综合| 18禁在线播放成人免费| 六月丁香七月| 一夜夜www| 男女边摸边吃奶| 精品久久久久久久久久久久久| 丰满乱子伦码专区| 大香蕉久久网| 国产久久久一区二区三区| 两个人视频免费观看高清| 亚洲精品日韩在线中文字幕| 99久久精品一区二区三区| 久久精品久久精品一区二区三区| 日本免费在线观看一区| 精品国内亚洲2022精品成人| 国产精品美女特级片免费视频播放器| 麻豆av噜噜一区二区三区| 久久草成人影院| videos熟女内射| 免费人成在线观看视频色| 欧美激情久久久久久爽电影| 欧美日本视频| 超碰97精品在线观看| 国产黄片视频在线免费观看| www.色视频.com| 国产成人午夜福利电影在线观看| 最近2019中文字幕mv第一页| 69人妻影院| av在线天堂中文字幕| 国产 一区精品| 国产伦理片在线播放av一区| 青春草国产在线视频| 能在线免费观看的黄片| 久久久久久伊人网av| 午夜福利在线在线| 99久久精品一区二区三区| 狠狠精品人妻久久久久久综合| 午夜福利在线观看吧| 久久久成人免费电影| 国产熟女欧美一区二区| 人妻一区二区av| 成人亚洲欧美一区二区av| 精品一区二区三区人妻视频| 国产亚洲5aaaaa淫片| av在线蜜桃| 国产大屁股一区二区在线视频| 国产精品久久久久久精品电影小说 | 青春草国产在线视频| 亚洲成人久久爱视频| 最后的刺客免费高清国语| 欧美成人一区二区免费高清观看| 久久久成人免费电影| 亚洲国产高清在线一区二区三| 赤兔流量卡办理| 亚洲精品乱码久久久v下载方式| 18禁在线无遮挡免费观看视频| 国产精品日韩av在线免费观看| 成人特级av手机在线观看| 亚洲熟妇中文字幕五十中出| 亚洲在久久综合| 国产亚洲一区二区精品| 在线 av 中文字幕| 赤兔流量卡办理| 免费av不卡在线播放| 亚洲精品久久久久久婷婷小说| 97超视频在线观看视频| 免费观看性生交大片5| 大片免费播放器 马上看| 久久精品综合一区二区三区| 91狼人影院| 国产v大片淫在线免费观看| 亚洲一级一片aⅴ在线观看| 白带黄色成豆腐渣| 亚洲av中文字字幕乱码综合| 日韩强制内射视频| 精品一区二区三卡| 好男人在线观看高清免费视频| 亚洲在线观看片| 国产精品1区2区在线观看.| 中文欧美无线码| 七月丁香在线播放| 国产精品国产三级专区第一集| 亚洲电影在线观看av| 国内精品美女久久久久久| 久久精品国产自在天天线| 亚洲欧洲国产日韩| 欧美三级亚洲精品| 大香蕉久久网| 人体艺术视频欧美日本| 国产高潮美女av| 国产成年人精品一区二区| 精品一区二区三区视频在线| 小蜜桃在线观看免费完整版高清| 欧美日韩精品成人综合77777| 亚洲aⅴ乱码一区二区在线播放| 韩国高清视频一区二区三区| 直男gayav资源| 欧美 日韩 精品 国产| 国产精品一二三区在线看| 99热6这里只有精品| 国产黄色小视频在线观看| 一个人免费在线观看电影| 两个人视频免费观看高清| 嫩草影院新地址| 人妻少妇偷人精品九色| 天堂影院成人在线观看| 中文字幕亚洲精品专区| 久久精品夜色国产| 久热久热在线精品观看| 女人被狂操c到高潮| 伊人久久国产一区二区| 天天躁夜夜躁狠狠久久av| 国产高清国产精品国产三级 | 色综合站精品国产| 美女高潮的动态| 亚洲性久久影院| 亚洲精品自拍成人| 亚洲国产精品专区欧美| 国产不卡一卡二| 国产一区二区三区av在线| 亚洲真实伦在线观看| 亚洲精品久久久久久婷婷小说| 国产精品福利在线免费观看| 97人妻精品一区二区三区麻豆| 99久久精品热视频| 性色avwww在线观看| 色尼玛亚洲综合影院| 欧美变态另类bdsm刘玥| 亚洲图色成人| 国产色婷婷99| 高清视频免费观看一区二区 | 久久久久久国产a免费观看| 一个人观看的视频www高清免费观看| 亚洲国产最新在线播放| 国产乱来视频区| 日本黄色片子视频| av线在线观看网站| 亚洲精品中文字幕在线视频 | 久久99蜜桃精品久久| 日本-黄色视频高清免费观看| 中文天堂在线官网| 日本wwww免费看| 最近手机中文字幕大全| 亚洲精品一二三| 2022亚洲国产成人精品| 性色avwww在线观看| av在线老鸭窝| 天堂中文最新版在线下载 | 亚洲精华国产精华液的使用体验| 亚洲,欧美,日韩| 日韩欧美一区视频在线观看 | 九草在线视频观看| 国产精品人妻久久久久久| 国产精品1区2区在线观看.| av在线观看视频网站免费| 久久久精品欧美日韩精品| 特大巨黑吊av在线直播| 亚洲电影在线观看av| 国产视频首页在线观看| 国产人妻一区二区三区在| 欧美日韩亚洲高清精品| 亚洲av不卡在线观看| 极品少妇高潮喷水抽搐| 久久亚洲国产成人精品v| 人妻夜夜爽99麻豆av| 国内精品美女久久久久久| 亚洲精品,欧美精品| 日韩制服骚丝袜av| 日本av手机在线免费观看| 高清在线视频一区二区三区| 国产av在哪里看| 毛片一级片免费看久久久久| 熟妇人妻不卡中文字幕| 人妻一区二区av| 美女脱内裤让男人舔精品视频| 午夜福利在线在线| 男的添女的下面高潮视频| 在线免费观看的www视频| 欧美精品国产亚洲| 久久久精品94久久精品| 色播亚洲综合网| 在现免费观看毛片| 久久久欧美国产精品| 日韩中字成人| 97精品久久久久久久久久精品| 欧美激情久久久久久爽电影| 久久久久久国产a免费观看| 晚上一个人看的免费电影| 国产黄片美女视频| 美女主播在线视频| 少妇丰满av| 国产亚洲午夜精品一区二区久久 | 精品国产三级普通话版| 精品久久久久久电影网| 一夜夜www| 欧美xxxx性猛交bbbb| 精品亚洲乱码少妇综合久久| 22中文网久久字幕| 久久午夜福利片| 波多野结衣巨乳人妻| 国内揄拍国产精品人妻在线| 最近2019中文字幕mv第一页| 亚洲最大成人中文| 亚洲欧洲国产日韩| 国产成人a∨麻豆精品| 又爽又黄无遮挡网站| 国产 一区 欧美 日韩| av在线老鸭窝| 不卡视频在线观看欧美| 久久国产乱子免费精品| 精品一区二区免费观看| 午夜激情久久久久久久| 性插视频无遮挡在线免费观看| 久久精品国产亚洲网站| 日韩欧美一区视频在线观看 | 麻豆av噜噜一区二区三区| 97超碰精品成人国产| 午夜福利网站1000一区二区三区| 亚洲av国产av综合av卡| 久久久久精品久久久久真实原创| 久久精品综合一区二区三区| 午夜福利在线观看吧| 青春草国产在线视频| 午夜精品国产一区二区电影 | av女优亚洲男人天堂| 真实男女啪啪啪动态图| 久久精品国产自在天天线| 成人性生交大片免费视频hd| 国产高清不卡午夜福利| av在线播放精品| 青春草视频在线免费观看| 中文欧美无线码| 国产精品久久视频播放| 九九爱精品视频在线观看| 久久精品国产鲁丝片午夜精品| 国产精品伦人一区二区| 国产高清三级在线| 欧美日本视频| 在线观看av片永久免费下载| 91午夜精品亚洲一区二区三区| 最近最新中文字幕免费大全7| 好男人在线观看高清免费视频| 国产欧美日韩精品一区二区| 婷婷六月久久综合丁香| 又爽又黄无遮挡网站| 三级国产精品欧美在线观看| 欧美成人一区二区免费高清观看| 精品酒店卫生间| 国产亚洲91精品色在线| 久久精品国产亚洲网站| 高清视频免费观看一区二区 | av在线观看视频网站免费| 久久久久久久久大av| 国产黄色小视频在线观看| av女优亚洲男人天堂| 中国美白少妇内射xxxbb| 欧美日韩综合久久久久久| 国产精品无大码| 日本猛色少妇xxxxx猛交久久| 99热全是精品| 成人高潮视频无遮挡免费网站| 草草在线视频免费看| 久99久视频精品免费| 白带黄色成豆腐渣| 日韩av免费高清视频| 少妇熟女欧美另类| 国产高清三级在线| 久久这里有精品视频免费| 国产高清不卡午夜福利| 精品久久久久久电影网| 国产免费视频播放在线视频 | 久久久久久久久久成人| 麻豆久久精品国产亚洲av| 亚洲av.av天堂| 久久99热这里只频精品6学生| 欧美 日韩 精品 国产| 少妇高潮的动态图| 搡女人真爽免费视频火全软件| 亚洲精品国产成人久久av| 18禁在线无遮挡免费观看视频| 亚洲精品456在线播放app| 久久久久久久久久久免费av| 日韩中字成人| 最近中文字幕高清免费大全6| 欧美 日韩 精品 国产| 成人高潮视频无遮挡免费网站| 自拍偷自拍亚洲精品老妇| av一本久久久久| 99热6这里只有精品| 国产成人一区二区在线| 我的老师免费观看完整版| 国产av在哪里看| 色综合站精品国产| av女优亚洲男人天堂| 免费看不卡的av| 只有这里有精品99| 国产人妻一区二区三区在| av.在线天堂| 亚洲av电影在线观看一区二区三区 | 亚洲人成网站在线观看播放| 人妻制服诱惑在线中文字幕| 高清午夜精品一区二区三区| 亚洲自拍偷在线| 久久综合国产亚洲精品| 国产一级毛片在线| 亚洲av中文字字幕乱码综合| 婷婷色麻豆天堂久久| 免费观看性生交大片5| 国产91av在线免费观看| 91狼人影院| 亚洲精品aⅴ在线观看| 一级片'在线观看视频| 国产人妻一区二区三区在| 中文字幕久久专区| 全区人妻精品视频| 亚洲av免费在线观看| 成年人午夜在线观看视频 | 永久免费av网站大全| 最近视频中文字幕2019在线8| 国产精品嫩草影院av在线观看| 久久精品久久久久久久性| 日韩一区二区视频免费看| 天堂中文最新版在线下载 | 特级一级黄色大片| 国产午夜精品久久久久久一区二区三区| 久久久精品94久久精品| 亚洲国产成人一精品久久久| 白带黄色成豆腐渣| 亚洲真实伦在线观看| 舔av片在线| 又爽又黄无遮挡网站| 久久久欧美国产精品| 国产黄频视频在线观看| 国产高清有码在线观看视频| 久久草成人影院| 99久久九九国产精品国产免费| 丰满人妻一区二区三区视频av| av女优亚洲男人天堂| 午夜久久久久精精品| 久久久久久久久久久免费av| 国内揄拍国产精品人妻在线| 久久久久九九精品影院| 免费黄网站久久成人精品| 国产精品1区2区在线观看.| 99久久中文字幕三级久久日本| 99热全是精品| 午夜免费男女啪啪视频观看| 美女主播在线视频| 在线免费十八禁| 精品酒店卫生间| 亚洲av中文av极速乱| 精品午夜福利在线看| 在线免费十八禁| 国产在线男女| 亚洲av国产av综合av卡| 久久人人爽人人爽人人片va| 亚洲最大成人手机在线| 久久久欧美国产精品| 亚洲成人精品中文字幕电影| 欧美成人a在线观看| 啦啦啦韩国在线观看视频| 蜜臀久久99精品久久宅男| 久久久久久久午夜电影| 少妇熟女欧美另类| 91狼人影院| 一级毛片 在线播放| 亚洲欧美中文字幕日韩二区| 国产探花极品一区二区| 99久久精品热视频| 五月天丁香电影| 午夜免费观看性视频| av播播在线观看一区| 性色avwww在线观看| 久久久久久久久中文| 菩萨蛮人人尽说江南好唐韦庄| 好男人在线观看高清免费视频| 少妇熟女aⅴ在线视频| 超碰av人人做人人爽久久| 亚洲,欧美,日韩| 校园人妻丝袜中文字幕| 亚洲色图av天堂| 亚洲精品影视一区二区三区av| 黄色日韩在线| 91久久精品国产一区二区成人| 亚洲人成网站在线观看播放| 99热网站在线观看| 国产乱来视频区| 肉色欧美久久久久久久蜜桃 | eeuss影院久久| 一本久久精品| av在线亚洲专区| 亚洲欧美成人综合另类久久久| 精品人妻熟女av久视频| 青春草视频在线免费观看| 91aial.com中文字幕在线观看| 亚洲精品影视一区二区三区av| 成人综合一区亚洲| 婷婷色综合大香蕉| 又爽又黄无遮挡网站| 亚洲激情五月婷婷啪啪| 亚洲欧美清纯卡通| av天堂中文字幕网| 国产成年人精品一区二区| 亚洲精品第二区| 色综合色国产| 日韩三级伦理在线观看| 亚洲一区高清亚洲精品| av.在线天堂| 视频中文字幕在线观看| 特大巨黑吊av在线直播| 日本-黄色视频高清免费观看| 在线观看美女被高潮喷水网站| 欧美丝袜亚洲另类| 人人妻人人澡人人爽人人夜夜 | 成人特级av手机在线观看| 蜜桃亚洲精品一区二区三区| 成人毛片60女人毛片免费| 91在线精品国自产拍蜜月| 国产精品人妻久久久久久| 国产男人的电影天堂91| 又爽又黄无遮挡网站| 寂寞人妻少妇视频99o| 精品久久国产蜜桃| 日本黄色片子视频| 99视频精品全部免费 在线| 亚洲欧洲日产国产| 久久99蜜桃精品久久| 好男人视频免费观看在线| 一级毛片aaaaaa免费看小| 国产精品久久久久久精品电影| 在现免费观看毛片| 欧美激情久久久久久爽电影| 91久久精品国产一区二区三区| 男女那种视频在线观看| 免费电影在线观看免费观看| 中文字幕亚洲精品专区| 国产精品国产三级专区第一集| 亚洲四区av| 国产黄片美女视频| 亚洲欧美日韩卡通动漫| 免费av观看视频| 婷婷色综合大香蕉| 男人爽女人下面视频在线观看| 伦理电影大哥的女人| 久久久久久久午夜电影| 亚洲精品乱码久久久v下载方式| 日韩一本色道免费dvd| 美女内射精品一级片tv| 免费看光身美女| 成人美女网站在线观看视频| 国产白丝娇喘喷水9色精品| 能在线免费看毛片的网站| 亚洲成人av在线免费| 欧美zozozo另类| 91狼人影院| 午夜激情久久久久久久| av黄色大香蕉| 看非洲黑人一级黄片| 国产美女午夜福利| 夜夜看夜夜爽夜夜摸| 亚洲精品久久午夜乱码| 欧美潮喷喷水| 日韩三级伦理在线观看| 超碰97精品在线观看| 91久久精品电影网| 成人无遮挡网站| 日韩国内少妇激情av| 噜噜噜噜噜久久久久久91| 可以在线观看毛片的网站| 97超视频在线观看视频| 色尼玛亚洲综合影院| 亚洲国产精品sss在线观看| 三级毛片av免费| 亚洲欧美精品自产自拍| 97热精品久久久久久| 国产一区二区三区综合在线观看 | 亚洲综合色惰| 美女高潮的动态| 91午夜精品亚洲一区二区三区| 亚洲欧美清纯卡通| 免费大片黄手机在线观看| av在线天堂中文字幕| 联通29元200g的流量卡| 日韩不卡一区二区三区视频在线| 啦啦啦啦在线视频资源| 国产精品综合久久久久久久免费| 成人无遮挡网站| 人人妻人人澡人人爽人人夜夜 | 精品久久久精品久久久| 一个人看的www免费观看视频| 亚洲最大成人av| 秋霞伦理黄片| 高清日韩中文字幕在线| 午夜久久久久精精品| 黄片wwwwww| 久久精品国产亚洲网站| 国产高潮美女av| 2021少妇久久久久久久久久久| 国产精品久久久久久久电影| 欧美成人一区二区免费高清观看| 夫妻性生交免费视频一级片| 水蜜桃什么品种好| 婷婷色综合大香蕉| 欧美成人午夜免费资源| 亚洲怡红院男人天堂| 日日啪夜夜爽| 亚洲精品日韩av片在线观看| 三级国产精品片| 两个人的视频大全免费| 99久国产av精品国产电影| 亚洲国产高清在线一区二区三| 午夜福利视频精品| 亚洲最大成人中文| 校园人妻丝袜中文字幕| 国产探花极品一区二区| 五月玫瑰六月丁香| 99热这里只有是精品50| 国产不卡一卡二| av线在线观看网站| 在线观看一区二区三区| 超碰97精品在线观看| 亚洲欧美日韩无卡精品| 免费看日本二区| 成人国产麻豆网| 久久久久久久亚洲中文字幕| 国产av不卡久久| 午夜视频国产福利| 国产av不卡久久| 国产探花在线观看一区二区| 亚洲精品色激情综合| 观看免费一级毛片| 亚洲国产精品专区欧美| 久久久久久久久久久丰满| av国产免费在线观看| 欧美 日韩 精品 国产| 国产av在哪里看| 狂野欧美白嫩少妇大欣赏| 身体一侧抽搐| 国产在视频线精品| 国产精品一区二区性色av| 国产麻豆成人av免费视频| 美女脱内裤让男人舔精品视频| 少妇猛男粗大的猛烈进出视频 | 午夜视频国产福利| 简卡轻食公司| 日韩精品青青久久久久久| 少妇人妻精品综合一区二区| 永久免费av网站大全| 麻豆成人午夜福利视频| 亚洲精品aⅴ在线观看| 亚洲国产精品成人久久小说| 精品人妻偷拍中文字幕| 国产乱来视频区| 97人妻精品一区二区三区麻豆| 国产在线一区二区三区精| or卡值多少钱| av.在线天堂| 国产精品久久久久久久久免| 日韩av免费高清视频| 国产白丝娇喘喷水9色精品| 亚洲精品国产av成人精品| 成人毛片a级毛片在线播放| 日韩一区二区三区影片| 国产黄色免费在线视频| 男人狂女人下面高潮的视频| 2022亚洲国产成人精品| 干丝袜人妻中文字幕| 国产黄色免费在线视频| 成人二区视频| 婷婷色综合大香蕉| 伊人久久国产一区二区| 欧美bdsm另类| 午夜久久久久精精品| 久久精品国产自在天天线| 中文字幕制服av| 国产精品久久久久久精品电影| 老师上课跳d突然被开到最大视频| 色吧在线观看| 可以在线观看毛片的网站| 舔av片在线| 欧美不卡视频在线免费观看| 成人性生交大片免费视频hd| 亚洲国产精品专区欧美| 亚洲性久久影院| av女优亚洲男人天堂| 日韩制服骚丝袜av| 一级毛片aaaaaa免费看小| 美女xxoo啪啪120秒动态图| 日韩,欧美,国产一区二区三区| 两个人视频免费观看高清| 三级男女做爰猛烈吃奶摸视频| videossex国产| 日本-黄色视频高清免费观看| 一区二区三区乱码不卡18| 国内精品美女久久久久久| 熟女人妻精品中文字幕| av在线观看视频网站免费| 精品99又大又爽又粗少妇毛片| 在线观看一区二区三区| 亚洲精品国产av蜜桃| 久久精品夜夜夜夜夜久久蜜豆| 国产色爽女视频免费观看| 精品久久久久久久久av| 人人妻人人看人人澡| 久久99热这里只频精品6学生| 色网站视频免费| 成年人午夜在线观看视频 | 国产麻豆成人av免费视频| 爱豆传媒免费全集在线观看| 国产成人一区二区在线| 18禁在线无遮挡免费观看视频| 国产免费一级a男人的天堂|