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

    湍流強(qiáng)度對水平軸風(fēng)力機(jī)氣動性能的影響

    2016-07-08 06:23:45李仁年李德順
    甘肅科學(xué)學(xué)報 2016年3期
    關(guān)鍵詞:數(shù)值模擬

    李仁年,任 鵬,李德順

    (1.蘭州理工大學(xué) 能源與動力工程學(xué)院,甘肅 蘭州 730050;2.甘肅省風(fēng)力機(jī)技術(shù)工程研究中心,甘肅 蘭州 730050)

    ?

    湍流強(qiáng)度對水平軸風(fēng)力機(jī)氣動性能的影響

    李仁年1,2,任鵬1,李德順1,2

    (1.蘭州理工大學(xué) 能源與動力工程學(xué)院,甘肅 蘭州730050;2.甘肅省風(fēng)力機(jī)技術(shù)工程研究中心,甘肅 蘭州730050)

    摘要為了研究水平軸風(fēng)力機(jī)氣動性能隨湍流強(qiáng)度的影響,基于CFD軟件對不同來流風(fēng)速工況下的33kW水平軸風(fēng)力機(jī)風(fēng)輪模型進(jìn)行三維數(shù)值模擬,對比分析風(fēng)力機(jī)在湍流強(qiáng)度I為0.1%、14%、25%時的氣動性能。結(jié)果表明:隨著來流湍流強(qiáng)度的增加,水平軸風(fēng)力機(jī)葉片表面壓差會有一定程度的減小,從而導(dǎo)致風(fēng)力機(jī)風(fēng)輪轉(zhuǎn)矩減小,風(fēng)力機(jī)風(fēng)能利用效率明顯降低。

    關(guān)鍵詞數(shù)值模擬;水平軸風(fēng)力機(jī);湍流強(qiáng)度;氣動性能

    近年來,化石能源的急速消耗與其產(chǎn)生的大氣污染問題受到了世界各國政府的高度重視。風(fēng)能作為一種清潔環(huán)保的可再生能源,它的開發(fā)利用成為解決這一世界難題不可或缺的途徑。風(fēng)能開發(fā)利用包括了風(fēng)場風(fēng)資源評估和風(fēng)力機(jī)設(shè)計選型等重要環(huán)節(jié)。其中,湍流強(qiáng)度是風(fēng)場風(fēng)資源評估的重要指標(biāo)之一,同時也是決定風(fēng)力機(jī)設(shè)計標(biāo)準(zhǔn)和安全性能的重要參數(shù)[1]。湍流強(qiáng)度反映了風(fēng)的湍流特性,風(fēng)的湍流擾動會在風(fēng)力機(jī)葉片上產(chǎn)生非定常的載荷波動,從而影響風(fēng)力機(jī)的氣動性能[2-5]。

    目前,研究風(fēng)力機(jī)氣動性能的主要方法已由傳統(tǒng)的二維葉素理論方法發(fā)展為三維流場數(shù)值模擬方法[6]。王強(qiáng)[7]基于CFD軟件,采用SST k-ω湍流模型對低風(fēng)速來流工況下的NREL PHASE VI風(fēng)力機(jī)葉片進(jìn)行全三維數(shù)值模擬,并將模擬結(jié)果與美國國家可再生能源實驗室的實驗數(shù)據(jù)進(jìn)行對比,結(jié)果表明風(fēng)力機(jī)葉片氣動性能可以通過三維數(shù)值模擬方法相對準(zhǔn)確地反映出來。黃宸武等[8,9]針對水平軸風(fēng)力機(jī)及其縮比模型,在相同雷諾數(shù)和尖速比的條件下利用CFD方法計算得到各模型的氣動性能,通過比較和分析各模型的流場信息,揭示了引起氣動性能差異的主要原因。王琦峰等[10]采用數(shù)值模擬的方法,分別使用S-A、標(biāo)準(zhǔn)k-ε、RNG k-ε和SST k-ω四種湍流模型,對水平軸風(fēng)力機(jī)的氣動性能進(jìn)行研究,并將計算結(jié)果與實驗結(jié)果進(jìn)行對比,認(rèn)為SST k-ω湍流模型可以很好地描述流場流動情況。

    我們以水平軸風(fēng)力機(jī)為研究對象,根據(jù)相似理論及風(fēng)力機(jī)的設(shè)計參數(shù),運用Pro/E軟件建立對應(yīng)的三維風(fēng)輪縮比模型,并運用FLUENT軟件對三維風(fēng)輪縮比模型進(jìn)行流場分析,計算不同湍流強(qiáng)度工況下的風(fēng)力機(jī)氣動性能,得到相應(yīng)的氣動性能曲線,對比分析湍流強(qiáng)度對水平軸風(fēng)力機(jī)性能的影響規(guī)律。

    1湍流強(qiáng)度

    湍流在空氣動力學(xué)中指的是短時間(一般少于10 min)內(nèi)的風(fēng)速波動。湍流產(chǎn)生的原因主要有兩個:一個是當(dāng)氣流流動時,由于地形差異(如山峰)造成的與地表的摩擦或者阻滯作用;另一個是因大氣溫度差異和空氣密度差異引起的氣流垂直流動。通常這兩種原因彼此影響[11]。

    湍流是一個復(fù)雜的隨機(jī)過程,可以用其統(tǒng)計規(guī)律進(jìn)行描述。湍流強(qiáng)度用來描述風(fēng)速隨時間變化的程度,反映脈動風(fēng)速的相對強(qiáng)度,是描述大氣湍流運動特性的最重要的特征量。湍流強(qiáng)度定義為標(biāo)準(zhǔn)風(fēng)速偏差與平均風(fēng)速的比率,即

    (1)

    地表的粗糙度和高度決定湍流強(qiáng)度的大小。但是,湍流強(qiáng)度通常會受到地貌特征的影響,如丘陵和山脈,以及位于上風(fēng)向的樹和建筑等障礙物。同時,它還受大氣熱運動的影響。例如,當(dāng)氣溫升高時,接近地面的空氣會穿過大氣上升,這將引起大規(guī)模的湍流對流。

    用于設(shè)計計算的湍流強(qiáng)度的確切值可以由其他的一些用于風(fēng)力機(jī)設(shè)計計算的標(biāo)準(zhǔn)來進(jìn)行規(guī)定,其縱向風(fēng)速分量標(biāo)準(zhǔn)偏差和湍流強(qiáng)度分別為

    σ=I15(15+aVhub)/(a+1),

    (2)

    (3)

    其中:I15為15 m/s風(fēng)速時計算的湍流強(qiáng)度特征值;a為斜度系數(shù);Vhub為輪轂高度平均風(fēng)速。對于較高湍流強(qiáng)度等級A,取I15=0.18,a=2;對于較低湍流強(qiáng)度等級B,取I15=0.16,a=3。

    IEC61400-1標(biāo)準(zhǔn)下風(fēng)場湍流強(qiáng)度曲線如圖1所示。等級A的湍流強(qiáng)度曲線高于等級B的湍流強(qiáng)度曲線,風(fēng)場湍流強(qiáng)度隨著風(fēng)速的增加而減小。

    2風(fēng)輪三維數(shù)值計算

    2.1計算模型與網(wǎng)格劃分

    選取某實驗風(fēng)場33kW兩葉片水平軸風(fēng)力機(jī)為原型建立幾何縮比模型。原型風(fēng)輪直徑為14.8 m,輪轂中心高16.1 m,切入風(fēng)速4 m/s,切出風(fēng)速為23 m/s,額定轉(zhuǎn)速為85 r/min,額定風(fēng)速為11 m/s,額定功率為33kW,葉片采用NACA4415-4430系列翼型[13]。根據(jù)相似理論,在保證模型和原型風(fēng)力機(jī)幾何相似的條件下,取模型風(fēng)輪直徑為500 mm。

    流場計算域為圓柱形,圓柱半徑為4倍風(fēng)輪半徑,上游延伸6倍風(fēng)輪半徑,下游延伸10倍風(fēng)輪半徑[14]。整體采用結(jié)構(gòu)化網(wǎng)格,網(wǎng)格單元數(shù)為845萬,網(wǎng)格質(zhì)量滿足計算要求。

    2.2計算方法和邊界條件

    數(shù)值模擬采用SST k-ω湍流模型,采用標(biāo)準(zhǔn)壁面函數(shù)求解近壁流動,壓力和速度耦合采用SIMPLE算法,相關(guān)變量均采用二階迎風(fēng)離散格式處理,旋轉(zhuǎn)區(qū)域采用多重參考系模型進(jìn)行計算。設(shè)置入口邊界條件為速度進(jìn)口,出口邊界條件為自由出流,葉片和輪轂表面均為無滑移壁面。

    在設(shè)置速度進(jìn)口邊界條件時,通過給定湍流強(qiáng)度和湍流粘度比來定義流場的湍流。在現(xiàn)代低湍流度風(fēng)洞中,自由流的湍流強(qiáng)度可由風(fēng)洞的特征長度估算出來,其湍流強(qiáng)度值約為0.1%。根據(jù)標(biāo)準(zhǔn),風(fēng)速在4~23 m/s之間時,風(fēng)場湍流強(qiáng)度大體處于14%~25%的范圍內(nèi),故分別取湍流強(qiáng)度為0.1%、14%、25%。在大部分外部流動的自由流邊界上,湍流粘度比的數(shù)值很小。一般情況下,湍流粘度比通常在1~10之間,故取其值為10。

    3模擬結(jié)果對比分析

    3.1葉片表面壓力分布

    分別取來流速度為7 m/s、11 m/s、23 m/s,風(fēng)輪轉(zhuǎn)速為2 516 r/min,對應(yīng)的葉尖速比分別為9.41、5.98、2.86。沿葉片展向在0.2R、0.6R和0.8R處做葉片截面翼型,計算得到各截面翼型在湍流強(qiáng)度分別為0.1%、14%和25%時的壓力分布曲線,如圖2~圖4所示。

    由圖2~圖4可以看出,在同一來流風(fēng)速下,從0.2R~0.8R截面,各截面翼型的壓力面壓力增大且吸力面壓力減小,翼型表面壓差增大;隨著來流風(fēng)速的增加,葉片壓力面壓力增大且吸力面壓力減小,葉片表面壓差增大。當(dāng)來流風(fēng)速為7 m/s時,沿展向0.2R處截面翼型吸力面在I=0.1%下的壓力值略小于I=14%和I=25%下的壓力值,同時,I=14%下的值略小于I=25%下的值;沿展向0.6R和0.8R處截面翼型分別位于葉片中部和葉片尖部,此時三種湍流強(qiáng)度對表面壓力影響不是很明顯。當(dāng)來流風(fēng)速為11 m/s時,湍流強(qiáng)度對葉片表面壓力的影響非常明顯,I=14%和I=25%時三處截面翼型的吸力面壓力值均比I=0.1%時的大,而壓力面的壓力值基本保持一致,導(dǎo)致葉片表面壓差變小,影響風(fēng)輪出力。當(dāng)來流風(fēng)速為23m/s時,沿展向0.2R處截面翼型吸力面靠近尾緣區(qū)域的壓力存在不規(guī)律性分布,說明翼型工作在大攻角下存在劇烈的流動分離,不同湍流強(qiáng)度的來流引起的分離區(qū)內(nèi)流動比較復(fù)雜,但總體上I=0.1%時該翼型表面的壓差要大于I=14%和I=25%時的壓差。

    綜上所述,湍流強(qiáng)度的變化會引起葉片表面壓力的變化,尤其在葉片根部表現(xiàn)最明顯,且在風(fēng)速為11 m/s時不同湍流強(qiáng)度對葉片表面壓力的影響最敏感。當(dāng)來流湍流強(qiáng)度增加時,截面翼型的吸力面壓力值增大,壓力面壓力值無明顯變化,翼型表面壓差減小,從而引起翼型升力的減小,降低風(fēng)輪的氣動性能。

    3.2轉(zhuǎn)矩性能

    圖5為湍流強(qiáng)度分別為0.1%、14%、25%時風(fēng)輪轉(zhuǎn)矩隨風(fēng)速變化的曲線圖,圖6為湍流強(qiáng)度分別為0.1%、14%、25%時風(fēng)輪風(fēng)能利用系數(shù)隨葉尖速比的變化曲線。

    由圖5可以看出,風(fēng)輪轉(zhuǎn)矩隨風(fēng)速增大而增大。風(fēng)速在4~7 m/s時三種湍流強(qiáng)度下的風(fēng)輪轉(zhuǎn)矩?zé)o明顯差異,說明在低風(fēng)速工況下風(fēng)輪轉(zhuǎn)矩對湍流強(qiáng)度的變化并不敏感。當(dāng)風(fēng)速超過7 m/s時風(fēng)輪轉(zhuǎn)矩會隨著湍流強(qiáng)度的增加而明顯減小,風(fēng)輪在高湍流強(qiáng)度I=25%時轉(zhuǎn)矩最小,這說明在中、高風(fēng)速工況下風(fēng)輪轉(zhuǎn)矩對湍流強(qiáng)度的變化很敏感,且湍流強(qiáng)度越大,風(fēng)輪出力越小。

    由圖6可以看出,風(fēng)輪的風(fēng)能利用系數(shù)隨葉尖速比的增大而先增大后減小。隨著來流湍流強(qiáng)度的增加,風(fēng)能利用系數(shù)會明顯遞減。當(dāng)葉尖速比為2.86,I=25%時風(fēng)輪風(fēng)能利用系數(shù)最低,比I=14%和I=0.1%時分別低5.4%和10.9%。風(fēng)輪隨著葉尖速比的增大受到來流湍流強(qiáng)度的影響越明顯,當(dāng)葉尖速比為5.98時,風(fēng)能利用系數(shù)達(dá)到最大值,此時不同湍流強(qiáng)度對風(fēng)輪風(fēng)能利用系數(shù)的影響最明顯,I=14%和I=25%時的風(fēng)輪風(fēng)能利用系數(shù)比I=0.1%時分別降低了12.4%和22.7%。

    綜上所述,來流湍流強(qiáng)度對水平軸風(fēng)力機(jī)的氣動性能影響明顯,隨著湍流強(qiáng)度的增加,風(fēng)力機(jī)風(fēng)輪轉(zhuǎn)矩有一定程度的降低,導(dǎo)致風(fēng)力機(jī)風(fēng)能利用效率降低。

    4結(jié)論

    通過對水平軸風(fēng)力機(jī)風(fēng)輪模型進(jìn)行三維定常湍流數(shù)值模擬,對比分析風(fēng)輪分別在湍流強(qiáng)度I為0.1%、14%、25%時的氣動性能,得到以下結(jié)論:

    (1)隨著來流湍流強(qiáng)度的增加,葉片吸力面壓力增大,壓力面壓力基本保持不變,導(dǎo)致葉片表面壓差減小,降低了風(fēng)輪的氣動性能,其中,葉根附近的表面壓力變化較為明顯;

    (2)隨著來流湍流強(qiáng)度的增加,風(fēng)輪轉(zhuǎn)矩和風(fēng)能利用系數(shù)明顯減小,水平軸風(fēng)力機(jī)氣動性能會發(fā)生明顯的降低,高湍流強(qiáng)度I=25%時風(fēng)輪的最大風(fēng)能利用性能比低湍流強(qiáng)度I=0.1%時低22.7%。

    參考文獻(xiàn)

    [1]王承凱.風(fēng)場湍流強(qiáng)度的計算及其對風(fēng)電機(jī)組選型的影響[C]//2008年中國電機(jī)工程學(xué)會年會論文集.西安:中國電機(jī)工程學(xué)會,2008.

    [2]宋麗莉,毛慧琴,錢光明,等.熱帶氣旋對風(fēng)力發(fā)電的影響分析[J].太陽能學(xué)報,2006,27(9):961-965.

    [3]歐陽華,巫發(fā)明,王靛,等.湍流強(qiáng)度對兆瓦級風(fēng)電機(jī)組的影響[J].大功率變流技術(shù),2013,36(3):82-86.

    [4]楊佳元,歐陽華,盛科,等.WT計算湍流強(qiáng)度用于風(fēng)電場安全性評估的研究[J].大功率變流技術(shù),2013,36(5):14-16.

    [5]???李國文,王虎彬,等.湍流強(qiáng)度對S型風(fēng)力機(jī)氣動性能影響的研究[J].沈陽航空航天大學(xué)學(xué)報,2012,29(4):25-28.

    [6]李宇紅,張慶麟.風(fēng)力機(jī)葉片三維流動特性與氣動性能的數(shù)值分析[J].太陽能學(xué)報,2008,29(9):1 172-1 176.

    [7]王強(qiáng).水平軸風(fēng)力機(jī)氣動性能數(shù)值模擬研究[D].南京:南京航空航天大學(xué),2010.

    [8]黃宸武.基于相似理論風(fēng)力機(jī)氣動性能預(yù)測研究[D].北京:中國科學(xué)院大學(xué),2012.

    [9]黃宸武,楊科,劉強(qiáng),等.相同雷諾數(shù)和尖速比下縮比風(fēng)力機(jī)流動特性分析[J].工程熱物理學(xué)報,2013,34(3):440-444.

    [10]王琦峰,姜哲,安亦然.基于NREL PHASE VI風(fēng)力機(jī)數(shù)值模擬的湍流模式研究[J].水動力學(xué)研究與進(jìn)展,2013,28(5):606-614.

    [11]Tony Burton,Nick Jenkins,David Sharpe.風(fēng)能技術(shù)[M].第2版.北京:科學(xué)出版社,2014.

    [12]胡隱樵,陳晉北,左洪超.湍流強(qiáng)度定理和湍流發(fā)展的宏觀機(jī)制[J].中國科學(xué),2007,37(2):272-281.

    [13]李德順,李仁年,王秀勇.基于外場實驗的風(fēng)力機(jī)葉片三維效應(yīng)研究[J].應(yīng)用數(shù)學(xué)和力學(xué),2013,34(10):1 073-1 082.

    [14]范忠瑤,康順,王建錄.風(fēng)力機(jī)葉片三維數(shù)值計算方法確認(rèn)研究[J].太陽能學(xué)報,2010,31(3):279-285.

    he Influence of Turbulence Intensity on Aerodynamic Performance of Horizontal Axis Wind Turbine

    Li Rennian1,2,Ren Peng1,Li Deshun1,2

    (1.School of Energy and Power Engineering,Lanzhou University of Technology,Lanzhou 730050,China;2.Gansu Provincial Engineering Center for Wind Turbines,Lanzhou 730050,China)

    Key wordsNumerical simulation;Horizontal axis wind turbine;Turbulence intensity;Aerodynamic performance

    AbstractIn order to study the influence of turbulence intensity on aerodynamic performance of horizontal axis wind turbine,a three-dimensionally numerical simulation of wind wheel model of 33kWhorizontal axis wind turbine under different working conditions of wind speed of incoming flow is conducted based on CFD software by comparing and analyzing aerodynamic performance of wind turbine when the turbulence intensity(I) is 0.1%,14% and 25%.The results show that differential pressure on the surface of horizontal axis wind turbine blade decreases on a certain degree with the increase of turbulence intensity,which causes that torque of wind wheel of wind turbine decreases and wind power utilization efficiency of wind turbine is clearly reduced.

    doi:10.16468/j.cnki.issn1004-0366.2016.03.010.

    收稿日期:2015-04-08;修回日期:2015-05-11.

    基金項目:國家自然科學(xué)基金(51166009).

    作者簡介:李仁年(1963-),男,甘肅民勤人,教授,博士生導(dǎo)師,研究方向為流體機(jī)械.E-mail:fox925@yeah.net. 通訊作者:任鵬.E-mail:rp72433@163.com.

    中圖分類號:TK83

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1004-0366(2016)03-0044-05

    引用格式:Li Rennian,Ren Peng,Li Deshun.The Influence of Turbulence Intensity on Aerodynamic Performance of Horizontal Axis Wind Turbine[J].Journal of Gansu Sciences,2016,28(3):44-48.[李仁年,任鵬,李德順.湍流強(qiáng)度對水平軸風(fēng)力機(jī)氣動性能的影響[J].甘肅科學(xué)學(xué)報,2016,28(3):44-48.]

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測
    鉆孔灌注樁樁底沉渣對樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測試、預(yù)處理和加工研究報告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    黑丝袜美女国产一区| 国产成人精品婷婷| 美女国产视频在线观看| 一级,二级,三级黄色视频| 超色免费av| 日本与韩国留学比较| 黄色配什么色好看| 成人影院久久| 亚洲av成人精品一二三区| 国产一区二区三区av在线| 蜜臀久久99精品久久宅男| 一区在线观看完整版| 久久久久精品久久久久真实原创| 中文欧美无线码| 五月天丁香电影| 91在线精品国自产拍蜜月| 久久久久人妻精品一区果冻| 国产黄色免费在线视频| 久久久久久久亚洲中文字幕| 18禁动态无遮挡网站| 菩萨蛮人人尽说江南好唐韦庄| 黄色欧美视频在线观看| 99热这里只有是精品在线观看| 亚洲,一卡二卡三卡| 国产免费现黄频在线看| 久久久欧美国产精品| 亚洲精品456在线播放app| 久久久久人妻精品一区果冻| 欧美精品一区二区免费开放| 人体艺术视频欧美日本| 纯流量卡能插随身wifi吗| 最近手机中文字幕大全| 色哟哟·www| 日本wwww免费看| 美女大奶头黄色视频| 免费高清在线观看日韩| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 少妇的逼好多水| 3wmmmm亚洲av在线观看| 国产女主播在线喷水免费视频网站| 国产成人午夜福利电影在线观看| 午夜福利网站1000一区二区三区| 色婷婷av一区二区三区视频| 国产亚洲最大av| 精品少妇内射三级| 亚洲av不卡在线观看| 日韩在线高清观看一区二区三区| 欧美亚洲 丝袜 人妻 在线| 美女国产视频在线观看| 国模一区二区三区四区视频| 亚洲av在线观看美女高潮| 久久av网站| 交换朋友夫妻互换小说| 又粗又硬又长又爽又黄的视频| 丰满少妇做爰视频| 久久精品国产鲁丝片午夜精品| 十八禁网站网址无遮挡| 国产一区二区三区综合在线观看 | 国产在线视频一区二区| 亚洲欧美色中文字幕在线| 久久av网站| 满18在线观看网站| 成人国产av品久久久| 久久国产精品男人的天堂亚洲 | 日韩伦理黄色片| 成人国产av品久久久| 免费观看在线日韩| 久久人人爽人人爽人人片va| 国产精品一区www在线观看| 国产一区有黄有色的免费视频| 黄色一级大片看看| 国产伦精品一区二区三区视频9| 美女国产高潮福利片在线看| 精品亚洲成a人片在线观看| 日韩成人伦理影院| 国产一区二区三区av在线| 99热国产这里只有精品6| 国产亚洲欧美精品永久| 亚洲成人手机| 99热这里只有是精品在线观看| 久久毛片免费看一区二区三区| 中文天堂在线官网| 国产亚洲最大av| 在现免费观看毛片| 两个人免费观看高清视频| 日韩av不卡免费在线播放| 国产精品三级大全| a 毛片基地| 99久久精品一区二区三区| 久久久久久久久久成人| 日韩av在线免费看完整版不卡| 自拍欧美九色日韩亚洲蝌蚪91| 一级毛片我不卡| 高清在线视频一区二区三区| 狠狠婷婷综合久久久久久88av| 日韩大片免费观看网站| 女性被躁到高潮视频| 男人添女人高潮全过程视频| 午夜日本视频在线| 日韩成人av中文字幕在线观看| 丝袜在线中文字幕| 国产av精品麻豆| 女人精品久久久久毛片| 曰老女人黄片| 欧美激情国产日韩精品一区| 日韩制服骚丝袜av| 男人爽女人下面视频在线观看| 亚洲三级黄色毛片| 亚洲激情五月婷婷啪啪| 最近中文字幕高清免费大全6| 久热久热在线精品观看| 最近中文字幕2019免费版| 亚洲国产av新网站| 美女国产高潮福利片在线看| 亚洲第一av免费看| 免费高清在线观看视频在线观看| 精品少妇内射三级| 51国产日韩欧美| 国产成人aa在线观看| 免费播放大片免费观看视频在线观看| 亚洲欧美精品自产自拍| 男女免费视频国产| 国产极品天堂在线| 十分钟在线观看高清视频www| 九色成人免费人妻av| 久久人人爽人人爽人人片va| 国产男人的电影天堂91| 国产伦精品一区二区三区视频9| 免费观看a级毛片全部| 国产熟女午夜一区二区三区 | 91在线精品国自产拍蜜月| 亚洲国产欧美在线一区| 日韩一本色道免费dvd| 亚洲,欧美,日韩| 久久精品熟女亚洲av麻豆精品| 蜜桃久久精品国产亚洲av| 国产亚洲精品第一综合不卡 | 午夜影院在线不卡| 日韩在线高清观看一区二区三区| 女人久久www免费人成看片| 天堂中文最新版在线下载| 22中文网久久字幕| 国产午夜精品久久久久久一区二区三区| 搡女人真爽免费视频火全软件| 免费久久久久久久精品成人欧美视频 | 男女免费视频国产| 青青草视频在线视频观看| 欧美激情极品国产一区二区三区 | 国模一区二区三区四区视频| 80岁老熟妇乱子伦牲交| 日韩精品有码人妻一区| 亚洲国产精品成人久久小说| 日韩,欧美,国产一区二区三区| 日韩不卡一区二区三区视频在线| 久久久久精品性色| www.av在线官网国产| 亚洲综合色惰| 麻豆乱淫一区二区| 人人妻人人澡人人爽人人夜夜| 色吧在线观看| 建设人人有责人人尽责人人享有的| 精品少妇黑人巨大在线播放| 99久久中文字幕三级久久日本| a级毛片黄视频| 免费看不卡的av| 亚洲天堂av无毛| 一区二区三区免费毛片| 免费看不卡的av| 国产极品粉嫩免费观看在线 | 国产免费现黄频在线看| 男女高潮啪啪啪动态图| av线在线观看网站| 久久免费观看电影| 美女内射精品一级片tv| 在线播放无遮挡| 日韩伦理黄色片| 欧美 亚洲 国产 日韩一| h视频一区二区三区| 国产精品偷伦视频观看了| 嫩草影院入口| 97超碰精品成人国产| av福利片在线| 婷婷色综合www| 精品久久久噜噜| 国产精品国产三级专区第一集| 欧美日韩视频高清一区二区三区二| 乱人伦中国视频| 寂寞人妻少妇视频99o| 狂野欧美激情性xxxx在线观看| 国产精品蜜桃在线观看| 纵有疾风起免费观看全集完整版| 中文字幕av电影在线播放| 久久久国产欧美日韩av| 观看av在线不卡| 中文字幕精品免费在线观看视频 | 国产精品国产三级专区第一集| 少妇猛男粗大的猛烈进出视频| 国产色婷婷99| 日本午夜av视频| 午夜福利影视在线免费观看| 欧美人与善性xxx| 国产成人免费观看mmmm| 超碰97精品在线观看| 亚洲成色77777| 九色成人免费人妻av| 美女视频免费永久观看网站| 交换朋友夫妻互换小说| 99久久综合免费| 一区二区日韩欧美中文字幕 | 我的女老师完整版在线观看| 欧美亚洲日本最大视频资源| 亚洲精品aⅴ在线观看| 校园人妻丝袜中文字幕| 视频区图区小说| 国产老妇伦熟女老妇高清| 欧美人与善性xxx| 中国三级夫妇交换| 91精品国产九色| 成年美女黄网站色视频大全免费 | 亚洲人成77777在线视频| 精品久久久久久久久av| 精品亚洲乱码少妇综合久久| 日本av免费视频播放| 美女cb高潮喷水在线观看| 在线播放无遮挡| 日日撸夜夜添| 国产欧美日韩一区二区三区在线 | 欧美亚洲日本最大视频资源| 精品人妻偷拍中文字幕| 精品久久久精品久久久| 精品少妇黑人巨大在线播放| av一本久久久久| 久久久久久人妻| 男女啪啪激烈高潮av片| 日韩欧美精品免费久久| 97超碰精品成人国产| 99久久中文字幕三级久久日本| 天堂俺去俺来也www色官网| 国产有黄有色有爽视频| 高清不卡的av网站| 亚洲人成网站在线观看播放| 亚洲一级一片aⅴ在线观看| 一级黄片播放器| 欧美最新免费一区二区三区| 国产亚洲最大av| 涩涩av久久男人的天堂| 精品久久久久久久久亚洲| 欧美日韩av久久| 国产精品.久久久| 精品久久久久久久久av| 国产高清三级在线| 日韩精品有码人妻一区| 99热这里只有是精品在线观看| 国产成人精品在线电影| 中文精品一卡2卡3卡4更新| 夜夜骑夜夜射夜夜干| 麻豆成人av视频| 国产一区二区三区av在线| 高清不卡的av网站| 在线精品无人区一区二区三| 国产高清国产精品国产三级| 久久狼人影院| 各种免费的搞黄视频| 婷婷成人精品国产| 日韩一本色道免费dvd| 国产在视频线精品| 日韩熟女老妇一区二区性免费视频| 搡老乐熟女国产| 久久国产精品男人的天堂亚洲 | 高清不卡的av网站| 国产av一区二区精品久久| 日韩伦理黄色片| 麻豆乱淫一区二区| 亚洲精品乱码久久久久久按摩| 午夜91福利影院| 91精品三级在线观看| 午夜久久久在线观看| 精品亚洲成a人片在线观看| 国产熟女午夜一区二区三区 | 欧美人与性动交α欧美精品济南到 | 久久人妻熟女aⅴ| 久久97久久精品| 欧美变态另类bdsm刘玥| 有码 亚洲区| 免费黄频网站在线观看国产| 26uuu在线亚洲综合色| 久久精品夜色国产| 国产成人精品无人区| 美女福利国产在线| 十八禁网站网址无遮挡| av一本久久久久| 青春草视频在线免费观看| 极品人妻少妇av视频| 亚洲av综合色区一区| 日韩一本色道免费dvd| 国产亚洲精品久久久com| 岛国毛片在线播放| 国产欧美亚洲国产| 最近中文字幕高清免费大全6| 国产成人午夜福利电影在线观看| 国国产精品蜜臀av免费| 久久久a久久爽久久v久久| 99热国产这里只有精品6| 性色av一级| 男女边摸边吃奶| 中文乱码字字幕精品一区二区三区| 人体艺术视频欧美日本| 精品国产国语对白av| 欧美日本中文国产一区发布| 国产综合精华液| 午夜视频国产福利| 国产精品欧美亚洲77777| 久久久a久久爽久久v久久| 亚洲怡红院男人天堂| 美女视频免费永久观看网站| 免费观看a级毛片全部| 天堂俺去俺来也www色官网| 一级,二级,三级黄色视频| 成年人午夜在线观看视频| 中文字幕人妻丝袜制服| 2021少妇久久久久久久久久久| 久久久久久久久久人人人人人人| 男男h啪啪无遮挡| 亚洲久久久国产精品| 国产视频首页在线观看| 国产成人午夜福利电影在线观看| 大片电影免费在线观看免费| 麻豆精品久久久久久蜜桃| 在线精品无人区一区二区三| 欧美激情国产日韩精品一区| 十分钟在线观看高清视频www| 亚洲欧美一区二区三区国产| 亚洲中文av在线| 啦啦啦在线观看免费高清www| 天堂俺去俺来也www色官网| 插逼视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 日日爽夜夜爽网站| 在线观看免费视频网站a站| 3wmmmm亚洲av在线观看| 在现免费观看毛片| 人妻制服诱惑在线中文字幕| 国产国拍精品亚洲av在线观看| 国产精品久久久久久久久免| 国产成人精品久久久久久| 欧美另类一区| 欧美日韩亚洲高清精品| 日韩一本色道免费dvd| 亚洲精品国产色婷婷电影| 狂野欧美激情性xxxx在线观看| 国产爽快片一区二区三区| 这个男人来自地球电影免费观看 | 久久97久久精品| 精品一区二区三卡| 3wmmmm亚洲av在线观看| 黄色视频在线播放观看不卡| 大香蕉久久网| 黑人猛操日本美女一级片| 18禁观看日本| 22中文网久久字幕| 毛片一级片免费看久久久久| 秋霞伦理黄片| 视频区图区小说| 亚洲欧美日韩另类电影网站| 久久女婷五月综合色啪小说| 在线观看国产h片| 中文字幕制服av| 成人毛片a级毛片在线播放| 日本与韩国留学比较| 国产免费一区二区三区四区乱码| 国产精品不卡视频一区二区| 一本—道久久a久久精品蜜桃钙片| 欧美丝袜亚洲另类| 人妻夜夜爽99麻豆av| 高清不卡的av网站| 午夜视频国产福利| 七月丁香在线播放| 热re99久久精品国产66热6| 大话2 男鬼变身卡| av国产精品久久久久影院| 欧美97在线视频| 看非洲黑人一级黄片| 少妇熟女欧美另类| 国产有黄有色有爽视频| 国产成人av激情在线播放 | 国产男女超爽视频在线观看| 多毛熟女@视频| 晚上一个人看的免费电影| 精品国产露脸久久av麻豆| 观看美女的网站| 精品人妻在线不人妻| www.色视频.com| 51国产日韩欧美| 亚洲国产日韩一区二区| 大片电影免费在线观看免费| av线在线观看网站| 男人操女人黄网站| 日日摸夜夜添夜夜爱| 亚洲丝袜综合中文字幕| 又大又黄又爽视频免费| 久久久久国产精品人妻一区二区| 久久久久国产网址| 亚洲精品色激情综合| a级毛色黄片| 欧美最新免费一区二区三区| 激情五月婷婷亚洲| 九色亚洲精品在线播放| 国产精品99久久久久久久久| 少妇人妻 视频| 欧美日韩国产mv在线观看视频| 亚洲丝袜综合中文字幕| 亚洲伊人久久精品综合| 亚洲成人一二三区av| 视频在线观看一区二区三区| 一本久久精品| 亚洲国产日韩一区二区| 搡老乐熟女国产| 男女无遮挡免费网站观看| 亚洲伊人久久精品综合| 欧美亚洲日本最大视频资源| 亚洲av国产av综合av卡| 夜夜骑夜夜射夜夜干| 成人亚洲精品一区在线观看| 国产免费视频播放在线视频| 99国产综合亚洲精品| av国产久精品久网站免费入址| 岛国毛片在线播放| 国产精品一区二区三区四区免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 熟妇人妻不卡中文字幕| 日产精品乱码卡一卡2卡三| 久久国内精品自在自线图片| 九草在线视频观看| 美女xxoo啪啪120秒动态图| h视频一区二区三区| 日韩三级伦理在线观看| 特大巨黑吊av在线直播| 日产精品乱码卡一卡2卡三| 亚洲国产精品成人久久小说| 我要看黄色一级片免费的| 国产精品一区二区在线观看99| 夜夜骑夜夜射夜夜干| 天天操日日干夜夜撸| a级毛色黄片| 精品国产乱码久久久久久小说| av在线app专区| 大片电影免费在线观看免费| 亚洲性久久影院| 亚洲色图综合在线观看| 中国美白少妇内射xxxbb| 国产精品一区二区在线不卡| 午夜视频国产福利| 97在线视频观看| 秋霞伦理黄片| 亚洲欧洲日产国产| 免费人成在线观看视频色| 亚洲伊人久久精品综合| 丝袜在线中文字幕| 国产黄频视频在线观看| a级毛片黄视频| 黑丝袜美女国产一区| 制服人妻中文乱码| 日韩精品免费视频一区二区三区 | 肉色欧美久久久久久久蜜桃| 91精品国产国语对白视频| 熟女av电影| 日韩av免费高清视频| 精品午夜福利在线看| 在线亚洲精品国产二区图片欧美 | 欧美精品亚洲一区二区| 天天操日日干夜夜撸| 日本欧美国产在线视频| 欧美日韩一区二区视频在线观看视频在线| 国产免费视频播放在线视频| 99视频精品全部免费 在线| 亚洲情色 制服丝袜| 成人午夜精彩视频在线观看| 日韩三级伦理在线观看| 国产在线一区二区三区精| 婷婷色综合大香蕉| 国产69精品久久久久777片| 亚洲第一区二区三区不卡| 亚洲,一卡二卡三卡| 欧美变态另类bdsm刘玥| 色哟哟·www| 国产免费又黄又爽又色| 国产精品人妻久久久影院| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久久大av| 久久精品国产a三级三级三级| 精品国产一区二区久久| 免费黄网站久久成人精品| 视频区图区小说| 成年av动漫网址| 少妇被粗大猛烈的视频| 久久久久久久大尺度免费视频| 国产熟女午夜一区二区三区 | 美女大奶头黄色视频| 亚洲国产精品成人久久小说| 国产黄色免费在线视频| av线在线观看网站| 久久久精品94久久精品| 精品酒店卫生间| 色婷婷久久久亚洲欧美| 成人手机av| 精品亚洲成国产av| 国产日韩一区二区三区精品不卡 | 校园人妻丝袜中文字幕| 夫妻性生交免费视频一级片| 国产 精品1| 少妇被粗大猛烈的视频| 久久国产精品男人的天堂亚洲 | videossex国产| 国产免费一级a男人的天堂| 日本黄色片子视频| 成人国产av品久久久| 亚洲,欧美,日韩| 午夜福利影视在线免费观看| 日韩一本色道免费dvd| 女的被弄到高潮叫床怎么办| 大码成人一级视频| 在线精品无人区一区二区三| 亚洲四区av| 成人手机av| 一个人免费看片子| 久久综合国产亚洲精品| 国产有黄有色有爽视频| 啦啦啦啦在线视频资源| 十分钟在线观看高清视频www| 日韩成人伦理影院| 亚洲综合精品二区| 超色免费av| 精品亚洲成国产av| 人人澡人人妻人| 插逼视频在线观看| 青春草国产在线视频| 久久久久久伊人网av| 在线观看免费视频网站a站| 97在线视频观看| 人人妻人人澡人人爽人人夜夜| 国产精品嫩草影院av在线观看| 亚洲精品久久久久久婷婷小说| 免费黄网站久久成人精品| 国产精品欧美亚洲77777| 久久精品人人爽人人爽视色| 午夜激情福利司机影院| 久久久久久久久大av| 欧美激情 高清一区二区三区| 国产免费一区二区三区四区乱码| 蜜臀久久99精品久久宅男| 亚洲国产精品国产精品| 久久久久人妻精品一区果冻| 精品国产露脸久久av麻豆| 久久精品久久久久久久性| 久久久久视频综合| 国产成人freesex在线| 在线 av 中文字幕| 日韩精品有码人妻一区| 18+在线观看网站| 国产在视频线精品| 亚洲天堂av无毛| 精品视频人人做人人爽| 99久久综合免费| 超色免费av| 午夜影院在线不卡| 亚洲精品av麻豆狂野| 日韩一本色道免费dvd| 亚洲少妇的诱惑av| 亚洲情色 制服丝袜| 男女免费视频国产| 蜜桃久久精品国产亚洲av| 久久ye,这里只有精品| 最近中文字幕2019免费版| 久久精品久久久久久久性| 国产男人的电影天堂91| 成人免费观看视频高清| 日韩一区二区三区影片| 成人国语在线视频| 久久av网站| 国产老妇伦熟女老妇高清| 亚洲欧美成人综合另类久久久| 国产成人午夜福利电影在线观看| 春色校园在线视频观看| av福利片在线| 18+在线观看网站| 国产亚洲精品久久久com| 精品国产乱码久久久久久小说| 国产精品国产av在线观看| 久久午夜福利片| 久久99蜜桃精品久久| 黄色欧美视频在线观看| 国产成人a∨麻豆精品| 2022亚洲国产成人精品| 国产色爽女视频免费观看| 亚洲,一卡二卡三卡| 亚洲av中文av极速乱| 丝袜喷水一区| 亚洲精品久久午夜乱码| 一级,二级,三级黄色视频| 午夜福利视频在线观看免费| 亚洲欧美日韩卡通动漫| 亚洲精品久久久久久婷婷小说| 免费少妇av软件| 嫩草影院入口| 日本av免费视频播放| 精品人妻在线不人妻| 婷婷色综合www| 国产欧美日韩综合在线一区二区| 亚洲欧美成人精品一区二区| 久久国内精品自在自线图片| 成人国产麻豆网| 国产毛片在线视频| a 毛片基地| 国产精品国产三级国产av玫瑰| 看免费成人av毛片|