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

    適用于獨(dú)立變槳控制方式的風(fēng)機(jī)氣動(dòng)模型和風(fēng)模型的研究

    2011-04-28 07:49:38王維慶王海云
    水力發(fā)電 2011年4期
    關(guān)鍵詞:距角變槳主軸

    張 迪,王維慶,王海云,李 強(qiáng)

    (1.新疆大學(xué)電氣工程學(xué)院,新疆 烏魯木齊 830047;2.北京金風(fēng)科創(chuàng)風(fēng)電設(shè)備有限公司,北京 100176)

    0 引言

    風(fēng)機(jī)變槳控制分為集體變槳和獨(dú)立變槳,隨著風(fēng)機(jī)的容量,葉輪和塔架的尺寸的增加,如何降低風(fēng)機(jī)載荷的問(wèn)題日益凸顯出來(lái),獨(dú)立變槳同集體變槳相比的優(yōu)勢(shì)在于降低風(fēng)機(jī)載荷同時(shí)保證風(fēng)力發(fā)電機(jī)的輸出額定功率。獨(dú)立變槳是變槳控制發(fā)展的趨勢(shì)。少數(shù)國(guó)外先進(jìn)的風(fēng)機(jī)制造商已經(jīng)將獨(dú)立變槳應(yīng)用于兆瓦級(jí)風(fēng)機(jī)上,國(guó)內(nèi)先進(jìn)的風(fēng)機(jī)制造商也加大了對(duì)獨(dú)立變槳控制技術(shù)的研究。

    隨著風(fēng)機(jī)變槳控制方式從集體變槳過(guò)渡到獨(dú)立變槳,由于在傳統(tǒng)的應(yīng)用于3葉片集體變槳控制策略的風(fēng)機(jī)氣動(dòng)模型中,3只槳葉的槳距角對(duì)風(fēng)機(jī)氣動(dòng)轉(zhuǎn)矩的影響是耦合在一起的, 無(wú)法利用模型來(lái)分析單只葉片的槳距角對(duì)單只葉片的拍擊力和升力的影響,所以需要回到葉素理論中來(lái)研究適用于獨(dú)立變槳風(fēng)機(jī)氣動(dòng)模型的建立方法,因此需要一種新的建模方式。

    1 風(fēng)機(jī)氣動(dòng)模型的建立

    1.1 應(yīng)用于集體變槳的風(fēng)機(jī)氣動(dòng)模型

    [1]中的風(fēng)機(jī)風(fēng)能利用系數(shù)Cp的經(jīng)驗(yàn)計(jì)算公式為:

    式中,λ為葉尖速比;λi為中間變量;β為槳距角;C1=0.5176, C2=116, C3=0.4, C4=5, C5=21, C6=0.006 8。

    風(fēng)輪吸收的功率表達(dá)式為:

    式中,ρ為空氣密度;R為葉片半徑;V為風(fēng)速。風(fēng)輪產(chǎn)生的轉(zhuǎn)矩為:

    式中,ω為風(fēng)輪轉(zhuǎn)速。

    葉尖速比λ的計(jì)算公式:

    將式(4)和(5)代入式(3)中得出風(fēng)機(jī)氣動(dòng)轉(zhuǎn)矩 T 的表達(dá)式為:

    基于式 (1)、(2)、(6) 在 simulink 中搭建模型結(jié)構(gòu),仿真條件為槳距角β取0°,空氣密度為1.29 kg/m3,葉片半徑40 m;風(fēng)輪轉(zhuǎn)速ω從0 r/min變化到35 r/min, 風(fēng)速?gòu)? m/s到14 m/s。

    仿真結(jié)果:①風(fēng)機(jī)的輸出轉(zhuǎn)矩T變化情況如圖1中曲線(xiàn)所示;②根據(jù)文獻(xiàn)[1],Cp最大值對(duì)應(yīng)的減速比λopt為 8.1,帶入式(5)可以計(jì)算出最優(yōu)轉(zhuǎn)速ωopt和與其對(duì)應(yīng)的最優(yōu)轉(zhuǎn)矩Topt。用光滑的曲線(xiàn)連接形成圖1中虛線(xiàn)。虛線(xiàn)中的低風(fēng)速下段可以作為低風(fēng)速下最大功率跟蹤的T-ω參考曲線(xiàn)。

    圖1 槳距角為0°,不同風(fēng)速下葉輪轉(zhuǎn)矩隨轉(zhuǎn)速變化曲線(xiàn)

    當(dāng)風(fēng)速高于額定風(fēng)速時(shí),通過(guò)變槳電機(jī)驅(qū)動(dòng)槳葉增大槳距角,降低了Cp,從而降低風(fēng)機(jī)氣動(dòng)轉(zhuǎn)矩,維持風(fēng)輪轉(zhuǎn)速額定。圖2繪制了風(fēng)速為12 m/s、槳距角為0°,風(fēng)速為13 m/s、槳距角為1°,風(fēng)速14 m/s、槳距角5°轉(zhuǎn)矩同轉(zhuǎn)速變化的對(duì)比曲線(xiàn),可以看出增大槳距角對(duì)降低氣動(dòng)轉(zhuǎn)矩的效果是非常明顯的。

    圖2 變槳對(duì)氣動(dòng)轉(zhuǎn)矩的影響

    這種建模的方式可以用于對(duì)集體變槳策略的仿真驗(yàn)證,但是對(duì)于獨(dú)立變槳,由于模型沒(méi)有獨(dú)立的槳距角輸入,而只能輸入3個(gè)槳葉集體的槳距角,所以無(wú)法應(yīng)用于獨(dú)立變槳策略的驗(yàn)證。

    1.2 應(yīng)用于獨(dú)立變槳的風(fēng)機(jī)氣動(dòng)模型

    根據(jù)文獻(xiàn)[2]的葉素理論如圖3所示。若以葉素為參考系,風(fēng)將會(huì)疊加一個(gè)同葉素轉(zhuǎn)動(dòng)方向相反,大小相等的線(xiàn)速度-u,風(fēng)速V和-u矢量合為葉素的相對(duì)風(fēng)速w。w為葉素真正感受到的風(fēng)速。由于相對(duì)風(fēng)速w的作用,在風(fēng)機(jī)葉素產(chǎn)生了與w同一方向的阻力Fdrag與垂直于阻力的升力Flift。

    圖3 葉素受力分析

    式中,Cd為阻力系數(shù),Cl為升力系數(shù),同Cp一樣沒(méi)有量綱,二者的取值變化決定于攻角i,影響Cd、Cl的因素有彎曲度,葉片厚度,表面粗糙度和雷諾數(shù),Cd、Cl數(shù)據(jù)參考翼型為NACA0015;L為弦長(zhǎng);dr可以理解為一段很短的葉展。攻角變化范圍為0-2π。

    為了計(jì)算出作用在風(fēng)機(jī)主軸上的氣動(dòng)轉(zhuǎn)矩Tair和作用在葉片上的拍擊力Fflap,首先需要將葉素上的阻力和推力經(jīng)行軸向和切向分解與合成 (如圖4)。

    每一小段葉素上產(chǎn)生的升力dFu和拍擊力dFa由公式(10)和公式(11)計(jì)算出來(lái)

    圖4 Fdrag和Flift的分解與合成

    式中,I為入流角。

    為了實(shí)現(xiàn)仿真,將微分放大為變化量。將每只槳葉被分為均勻的10段,每段的弦長(zhǎng)為此段兩端弦長(zhǎng)的平均值,分別建立每段葉片的拍擊力△Fa和升力△Fu的模型。

    將一只葉片上每段的拍擊力△Fai與對(duì)輪轂處的轉(zhuǎn)矩△Fuiri求和,最終計(jì)算出每只葉片上的拍擊力F(flap)j與 3只葉片總體在主軸上的氣動(dòng)轉(zhuǎn)矩 Tair。

    式中,i為每段葉片的標(biāo)號(hào);ri為第i段葉片中心距離輪轂的距離;j為葉片數(shù)。

    基于式(12)~(17)在 simulink 中建立仿真模型,仿真條件為:空氣密度1.25 kg/m3, 葉片長(zhǎng)度40 m,風(fēng)速?gòu)? m/s至13 m/s,風(fēng)輪轉(zhuǎn)速?gòu)?到35 r/min,槳距角保持-1°。

    仿真結(jié)果如圖5,主軸氣動(dòng)轉(zhuǎn)矩隨著風(fēng)速,風(fēng)輪轉(zhuǎn)速變化曲線(xiàn)。由于槳距角為-1°,所以整體的主軸氣動(dòng)轉(zhuǎn)矩大于圖1中的主軸氣動(dòng)轉(zhuǎn)矩。

    圖5中橫坐標(biāo)為主軸轉(zhuǎn)速,縱坐標(biāo)為氣動(dòng)轉(zhuǎn)矩,曲線(xiàn)下方的面積為葉輪吸收的風(fēng)能,根據(jù)最大面積的方法同樣可以得到最佳轉(zhuǎn)速值以及對(duì)應(yīng)的最佳轉(zhuǎn)矩值作為低風(fēng)速下的最大功率跟蹤T-ω參考曲線(xiàn)。

    1.3 兩種建模方式的比較

    (1)第一種風(fēng)機(jī)氣動(dòng)模型是否可以反映風(fēng)機(jī)的氣動(dòng)轉(zhuǎn)矩關(guān)鍵在于對(duì)Cp(λ,β)曲線(xiàn)的擬合準(zhǔn)確度,對(duì)于不同類(lèi)型的葉片需要大量的數(shù)據(jù)進(jìn)行擬合,而第二種氣動(dòng)模型的建立只需要知道葉片的升力系數(shù)Cl和阻力系數(shù)Cd,這兩個(gè)參數(shù)可以在對(duì)應(yīng)的葉片參數(shù)中找到,能夠直接放入模型中,而且由于單只葉片模型由一段段葉片拼接而成,由此就能實(shí)現(xiàn)對(duì)單獨(dú)一只葉片的槳距角和風(fēng)速進(jìn)行獨(dú)立修改,可應(yīng)用于獨(dú)立變槳中。

    (2)第一種風(fēng)機(jī)氣動(dòng)模型無(wú)法計(jì)算單只葉片的拍擊力而第二種模型可以計(jì)算。

    (3)第二種模型的分段數(shù)可以增加,將會(huì)更加準(zhǔn)確模擬葉片,具有發(fā)展?jié)摿Α?/p>

    圖5 新模型不同風(fēng)速下葉輪轉(zhuǎn)矩隨轉(zhuǎn)速變化曲線(xiàn)

    2 風(fēng)剪塔影對(duì)風(fēng)機(jī)氣動(dòng)模型輸出的影響

    2.1 風(fēng)模型同風(fēng)機(jī)氣動(dòng)模型的影響與配合

    風(fēng)剪塔影最終對(duì)于風(fēng)機(jī)氣動(dòng)模型模型的影響包括:①在單只葉片上會(huì)產(chǎn)生1P的拍擊力波動(dòng),波動(dòng)頻率為葉輪旋轉(zhuǎn)周期的倒數(shù)。②在風(fēng)機(jī)主軸轉(zhuǎn)矩中產(chǎn)生3P的波動(dòng),波動(dòng)頻率為1P波動(dòng)的頻率的三倍。塔影引起的波動(dòng)要比風(fēng)剪更大。

    由于風(fēng)機(jī)的每只葉片被分為10段,風(fēng)速模型輸出的風(fēng)速也與這10段每一段中心的位置對(duì)應(yīng)。即風(fēng)速最終都應(yīng)化為葉片位置角度θ與葉片位置點(diǎn)到輪轂中心距離ri的函數(shù)從而風(fēng)速模型輸出30個(gè)風(fēng)速變量與3只葉片中的30個(gè)葉片段位置對(duì)應(yīng)。

    2.2 風(fēng)剪切以及與模型配合

    風(fēng)剪切是指當(dāng)隨著高度的不同,風(fēng)速隨之變化從而造成轉(zhuǎn)矩和拍擊力的周期性波動(dòng),參考文獻(xiàn)[3]中的風(fēng)剪切公式:

    式中,V(h)為某一高度處風(fēng)速;V(h0)為輪轂高度處風(fēng)速;α為風(fēng)剪切經(jīng)驗(yàn)指數(shù),取1/7。

    將式(18)轉(zhuǎn)化為每段葉片位置角度的函數(shù)如式(19):

    式中,ri為第i段葉片距離輪轂的半徑長(zhǎng)度;ω為風(fēng)輪轉(zhuǎn)速;根據(jù)式(20)和式(21),作用在單只葉片某段上的風(fēng)速,會(huì)疊加一個(gè)頻率為 (ω/60)Hz的分量。仿真設(shè)定條件為:葉片半徑40 m,輪轂高度70 m。轉(zhuǎn)速為19 r/min,槳距角為4°,輪轂處的風(fēng)速為13 m/s,觀(guān)察與一段葉片對(duì)應(yīng)的風(fēng)剪模型的輸出風(fēng)速如圖6。

    圖6 風(fēng)剪切效應(yīng)下的風(fēng)速波動(dòng)

    對(duì)風(fēng)速模型輸出的數(shù)據(jù)經(jīng)行快速傅里葉變化得到頻譜圖如圖7。

    圖7 風(fēng)速頻譜圖

    由圖7可以看出,風(fēng)速中除了直流分量外,還有明顯疊加了0.2 Hz諧波含量。

    風(fēng)剪切造成葉片拍擊力的1 P波動(dòng)如圖8,頻率為0.2 Hz,即等于葉輪旋轉(zhuǎn)周期5 s的倒數(shù),而仿真設(shè)定轉(zhuǎn)速為12 r/min,也就是0.2 r/s, 即5 s/r,即葉輪旋轉(zhuǎn)周期周期為5 s仿真結(jié)果驗(yàn)證了風(fēng)剪對(duì)單葉片拍擊力和氣動(dòng)轉(zhuǎn)矩影響的結(jié)論。

    圖8 風(fēng)剪切造成葉片拍擊力的1P波動(dòng)

    風(fēng)剪切造成主軸轉(zhuǎn)矩3P波動(dòng),頻率為0.6Hz如圖9。

    2.3 塔影效應(yīng)及與模型配合

    塔影效應(yīng)指塔架正面的風(fēng)速由于塔架的作用而減小,所以當(dāng)葉片經(jīng)過(guò)塔架的豎直位置時(shí),葉片上的拍擊力會(huì)減小,單只葉片旋轉(zhuǎn)一周后由風(fēng)速引起的拍擊力會(huì)產(chǎn)生1P波動(dòng)。而對(duì)于3葉片風(fēng)機(jī),葉輪旋轉(zhuǎn)一周后,每只葉片都經(jīng)過(guò)一次塔架位置,所以在主軸轉(zhuǎn)矩中會(huì)產(chǎn)生3P波動(dòng)。結(jié)合文獻(xiàn)[3]的塔影效應(yīng)計(jì)算公式得出每段葉片對(duì)應(yīng)位置的風(fēng)速表達(dá)式:

    圖9 風(fēng)剪切造成主軸轉(zhuǎn)矩3P波動(dòng)

    式中,a為塔架的半徑 ;R為葉片半徑;H為塔架高度;x為葉片與塔架的距離;塔影效應(yīng)只在位置角度θ為π到π內(nèi)起作用。仿真的其他條件不變考慮塔影效應(yīng)的情況下觀(guān)察到單只葉片的拍擊中產(chǎn)生了1P波動(dòng),主軸轉(zhuǎn)矩中產(chǎn)生3P波動(dòng),塔影效應(yīng)造成葉片拍擊力的1P波動(dòng),頻率為0.2 Hz,如圖10。

    圖10 塔影效應(yīng)造成葉片拍擊力的1P波動(dòng)

    塔影效應(yīng)造成主軸轉(zhuǎn)矩3P波動(dòng),頻率為0.6 Hz如圖11。

    圖11 塔影效應(yīng)造成主軸轉(zhuǎn)矩3P波動(dòng)

    2.4 風(fēng)剪塔影對(duì)風(fēng)機(jī)的影響

    風(fēng)剪塔影都會(huì)造成單只葉片的拍擊力中的1P波動(dòng),主軸轉(zhuǎn)矩中的3P波動(dòng),但是由于塔影效應(yīng)對(duì)風(fēng)速的影響比風(fēng)剪切更加嚴(yán)重。所以塔影效應(yīng)產(chǎn)生的波動(dòng)更加嚴(yán)重。

    3 獨(dú)立變槳對(duì)葉片拍擊力的影響

    為了減小風(fēng)對(duì)葉片拍擊力的影響,需要對(duì)每只葉片的槳距角進(jìn)行獨(dú)立的控制,風(fēng)速為13 m/s,轉(zhuǎn)速為12 r/min的條件下,3只槳葉的槳距角分別為2°,4°,6°,從仿真結(jié)果可以看出槳距角的增大對(duì)葉片的拍擊力的影響如圖12。

    圖12 獨(dú)立變槳對(duì)葉片拍擊力的影響

    4 結(jié)論

    (1)綜合比較了基于Cp經(jīng)驗(yàn)公式與基于葉素理論的風(fēng)機(jī)氣動(dòng)模型,通過(guò)分析兩種模型構(gòu)架思想和仿真結(jié)果,說(shuō)明了基于葉素理論的風(fēng)機(jī)氣動(dòng)模型應(yīng)用于獨(dú)立變槳的正確性。

    (2)在基于葉素理論的風(fēng)機(jī)氣動(dòng)模型與基于風(fēng)剪切塔影效應(yīng)理論的風(fēng)模型基礎(chǔ)上,通過(guò)分析二者聯(lián)合仿真的結(jié)果,驗(yàn)證了風(fēng)剪塔影導(dǎo)致單只葉片拍擊力的1P波動(dòng)和3只葉片產(chǎn)生主軸轉(zhuǎn)矩中的3P波動(dòng)。

    參考文獻(xiàn):

    [1] HeierS.GridIntegrationofWindEnergyConversionSystems[M].1998.

    [2] 葉杭冶.風(fēng)力發(fā)電機(jī)組的控制技術(shù)[M].機(jī)械工業(yè)出版社2002:53-33.

    [3] 李少林,張興,謝震,等.雙饋風(fēng)力發(fā)電系統(tǒng)3次功率脈動(dòng)的研究[J].電網(wǎng)技術(shù) 2010, 34(4):38-39.

    [4] Song S H,Jeong B C,Oh J H,et al.Simulation of output characteristics of rotor blades using a hardware-in-loop wind turbine simulator[J].IEEE Trans on Energy Conversion,2005,3(6):1791-1796.

    [5] Anderson C G,Richon J-B,Campbell T J.An Aerodynamic Moment-Controlled Surface for Gust Load Alleviation on Wind Turbine Rotors[J].IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY,1998,6(5):578-580.

    [6] Dolan D S L,Lehn P W.Simulation model of wind turbine 3p torque oscillations due to wind shear and tower shadow[J].IEEE Trans on Energy Conversion, 2006, 21(3):717-724.

    [7] Geyler M,Caselitz P.Individual Blade Pitch Control Design for Load Reduction on Large Wind Turbines[C]//Proceedings of European Wind Energy Conference,2007:1-4.

    [8] 董萍,吳捷,楊金明,等.風(fēng)力發(fā)電機(jī)組建模研究現(xiàn)狀[J].太陽(yáng)能學(xué)報(bào) 2004, 25(5):613-615.

    [9] Van Baars G E,Bongers P M M.Flexible wind turbine model validation[J].Wind Engineering,1992,16(4):247-256.

    [10] Muljadi E,Butterfield C P.Pitch-Controlled Variable-Speed Wind Turbine Generation[J].IEEE Trans on Industry Applications,2001,37(1):240-246.

    猜你喜歡
    距角變槳主軸
    基于數(shù)據(jù)驅(qū)動(dòng)的風(fēng)電機(jī)組最優(yōu)槳距角 辨識(shí)方法
    考慮風(fēng)速差異的風(fēng)電場(chǎng)減載方案與一次調(diào)頻策略
    風(fēng)力發(fā)電機(jī)組最優(yōu)槳距角自尋優(yōu)控制算法設(shè)計(jì)
    兆瓦級(jí)風(fēng)電機(jī)組變槳距系統(tǒng)設(shè)計(jì)
    雙主軸雙排刀復(fù)合機(jī)床的研制
    基于FANUC-31i外部一轉(zhuǎn)信號(hào)在三檔主軸定向中的應(yīng)用
    變速風(fēng)力發(fā)電機(jī)組變槳距建模與仿真
    基于最小化多變量的獨(dú)立變槳距控制研究
    應(yīng)對(duì)最大360mm×360mm的加工物研發(fā)了雙主軸·半自動(dòng)切割機(jī)※1「DAD3660」
    變槳距風(fēng)力機(jī)在全風(fēng)速段內(nèi)的槳距角控制策略仿真研究
    一二三四在线观看免费中文在| 亚洲成人国产一区在线观看| 十八禁高潮呻吟视频| 天天躁日日躁夜夜躁夜夜| 国产97色在线日韩免费| 伦理电影免费视频| 精品国产国语对白av| 丰满饥渴人妻一区二区三| 涩涩av久久男人的天堂| a级片在线免费高清观看视频| 午夜福利在线免费观看网站| 大型av网站在线播放| 男人操女人黄网站| 一区二区av电影网| 男女高潮啪啪啪动态图| 在线十欧美十亚洲十日本专区| 亚洲熟女毛片儿| avwww免费| 国产免费现黄频在线看| 男女下面插进去视频免费观看| 两性夫妻黄色片| av又黄又爽大尺度在线免费看| 日韩制服骚丝袜av| av片东京热男人的天堂| 777米奇影视久久| 首页视频小说图片口味搜索| 亚洲av日韩在线播放| 美女视频免费永久观看网站| 欧美在线黄色| 91麻豆精品激情在线观看国产 | 国产成人欧美| 久久国产精品男人的天堂亚洲| 国产精品九九99| 老司机影院毛片| 青青草视频在线视频观看| 他把我摸到了高潮在线观看 | 国产免费av片在线观看野外av| 日本精品一区二区三区蜜桃| 香蕉丝袜av| 久久久久国内视频| 国产国语露脸激情在线看| 亚洲国产中文字幕在线视频| 热re99久久国产66热| 97精品久久久久久久久久精品| 亚洲精品国产av蜜桃| 97在线人人人人妻| 成人av一区二区三区在线看 | 三上悠亚av全集在线观看| av在线app专区| 国产在线一区二区三区精| 欧美在线黄色| 亚洲久久久国产精品| 99精国产麻豆久久婷婷| 亚洲激情五月婷婷啪啪| 两性夫妻黄色片| 成人国产av品久久久| 亚洲国产av新网站| 夫妻午夜视频| 日日摸夜夜添夜夜添小说| av电影中文网址| 手机成人av网站| 亚洲精品国产av蜜桃| 极品少妇高潮喷水抽搐| 精品高清国产在线一区| 一本久久精品| 国产精品一区二区精品视频观看| 国产精品一二三区在线看| 黑人巨大精品欧美一区二区蜜桃| 妹子高潮喷水视频| 伊人久久大香线蕉亚洲五| 日韩视频在线欧美| 亚洲,欧美精品.| 黑丝袜美女国产一区| 日韩大码丰满熟妇| 免费观看a级毛片全部| 一本色道久久久久久精品综合| 国产一区二区 视频在线| 午夜91福利影院| 侵犯人妻中文字幕一二三四区| 新久久久久国产一级毛片| 精品人妻在线不人妻| 欧美精品人与动牲交sv欧美| 国产成人av激情在线播放| 蜜桃国产av成人99| 最近中文字幕2019免费版| 日韩大片免费观看网站| av不卡在线播放| 欧美97在线视频| 男女免费视频国产| 青春草亚洲视频在线观看| 丰满饥渴人妻一区二区三| 亚洲精品美女久久av网站| 免费在线观看黄色视频的| 91精品伊人久久大香线蕉| av天堂在线播放| 99热全是精品| 性高湖久久久久久久久免费观看| 99久久精品国产亚洲精品| 中文字幕精品免费在线观看视频| a级片在线免费高清观看视频| 欧美日韩视频精品一区| 亚洲精品国产av成人精品| 搡老岳熟女国产| 纵有疾风起免费观看全集完整版| 久久久久国产一级毛片高清牌| 欧美激情久久久久久爽电影 | 亚洲av日韩精品久久久久久密| 一区二区av电影网| 下体分泌物呈黄色| 亚洲男人天堂网一区| 国产成人影院久久av| 国产成人av激情在线播放| 狠狠精品人妻久久久久久综合| 国产一区二区三区在线臀色熟女 | 在线永久观看黄色视频| 后天国语完整版免费观看| 黄色毛片三级朝国网站| 国产在线观看jvid| 两性夫妻黄色片| 国产成人a∨麻豆精品| 亚洲av成人不卡在线观看播放网 | 亚洲成人免费电影在线观看| 女人久久www免费人成看片| 99热国产这里只有精品6| 国产免费一区二区三区四区乱码| 午夜免费观看性视频| 亚洲精品日韩在线中文字幕| 窝窝影院91人妻| 精品人妻1区二区| 亚洲精品中文字幕一二三四区 | 亚洲五月色婷婷综合| 免费在线观看完整版高清| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美 亚洲 国产 日韩一| 精品一品国产午夜福利视频| 99香蕉大伊视频| 精品久久久精品久久久| 欧美激情 高清一区二区三区| 精品国产乱码久久久久久小说| 黑人操中国人逼视频| 搡老岳熟女国产| 亚洲av国产av综合av卡| 亚洲av片天天在线观看| 另类亚洲欧美激情| 菩萨蛮人人尽说江南好唐韦庄| 在线十欧美十亚洲十日本专区| 久久天堂一区二区三区四区| 女人精品久久久久毛片| 亚洲精品在线美女| 大陆偷拍与自拍| 亚洲avbb在线观看| 久久久久久人人人人人| 亚洲九九香蕉| 国产片内射在线| 成人av一区二区三区在线看 | 亚洲成人免费电影在线观看| 精品国产乱码久久久久久男人| 午夜成年电影在线免费观看| 成人国产一区最新在线观看| 亚洲精品日韩在线中文字幕| 夜夜骑夜夜射夜夜干| 日韩欧美一区视频在线观看| 三级毛片av免费| 在线十欧美十亚洲十日本专区| 亚洲第一青青草原| 少妇人妻久久综合中文| 天堂俺去俺来也www色官网| 国产精品.久久久| av又黄又爽大尺度在线免费看| 国产精品国产三级国产专区5o| 大香蕉久久网| 丰满人妻熟妇乱又伦精品不卡| 蜜桃国产av成人99| 老鸭窝网址在线观看| 久久青草综合色| 日本91视频免费播放| 国产精品 国内视频| 啦啦啦中文免费视频观看日本| 丝瓜视频免费看黄片| 在线看a的网站| 国产精品一区二区免费欧美 | 精品一区在线观看国产| 伊人久久大香线蕉亚洲五| 大型av网站在线播放| xxxhd国产人妻xxx| 午夜福利乱码中文字幕| 日本一区二区免费在线视频| 丝袜脚勾引网站| 啦啦啦 在线观看视频| 日韩中文字幕视频在线看片| 亚洲国产成人一精品久久久| 亚洲一区二区三区欧美精品| 啦啦啦中文免费视频观看日本| 不卡一级毛片| 宅男免费午夜| 黑丝袜美女国产一区| 国产一区二区三区av在线| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩精品网址| 波多野结衣av一区二区av| 精品少妇久久久久久888优播| 秋霞在线观看毛片| 日韩精品免费视频一区二区三区| 大香蕉久久网| 亚洲人成电影观看| 制服人妻中文乱码| 日韩 亚洲 欧美在线| 91精品三级在线观看| 日本wwww免费看| 青春草亚洲视频在线观看| 在线观看免费日韩欧美大片| 首页视频小说图片口味搜索| 黑人操中国人逼视频| 女人高潮潮喷娇喘18禁视频| 欧美精品人与动牲交sv欧美| 美国免费a级毛片| 久久av网站| 国产区一区二久久| 国产精品一二三区在线看| 91av网站免费观看| 国精品久久久久久国模美| 成年动漫av网址| 国产欧美日韩综合在线一区二区| 麻豆乱淫一区二区| 丝袜美腿诱惑在线| 久久人人爽av亚洲精品天堂| 国产主播在线观看一区二区| 欧美黄色淫秽网站| 999久久久国产精品视频| 电影成人av| 色播在线永久视频| 免费在线观看完整版高清| 成年女人毛片免费观看观看9 | 国产91精品成人一区二区三区 | 亚洲精品av麻豆狂野| 在线观看一区二区三区激情| 亚洲欧美成人综合另类久久久| 777久久人妻少妇嫩草av网站| 一个人免费在线观看的高清视频 | 男人添女人高潮全过程视频| 这个男人来自地球电影免费观看| 黄片播放在线免费| 99久久99久久久精品蜜桃| 久久热在线av| 天天躁夜夜躁狠狠躁躁| 欧美黑人精品巨大| 国产一区二区激情短视频 | 成人av一区二区三区在线看 | 91九色精品人成在线观看| 丝袜喷水一区| 一级片免费观看大全| 久久人人爽av亚洲精品天堂| 不卡一级毛片| 欧美日韩中文字幕国产精品一区二区三区 | 午夜91福利影院| 一个人免费看片子| 中文字幕人妻丝袜一区二区| 欧美日韩一级在线毛片| 美女福利国产在线| 亚洲专区中文字幕在线| 日本精品一区二区三区蜜桃| av福利片在线| 天天躁日日躁夜夜躁夜夜| 国产在线观看jvid| 母亲3免费完整高清在线观看| www.熟女人妻精品国产| 又黄又粗又硬又大视频| 国产野战对白在线观看| 精品免费久久久久久久清纯 | av有码第一页| 国产极品粉嫩免费观看在线| 色播在线永久视频| 无遮挡黄片免费观看| 在线看a的网站| 国产人伦9x9x在线观看| 另类精品久久| 国产精品九九99| 国产区一区二久久| 五月天丁香电影| 久久精品aⅴ一区二区三区四区| 欧美亚洲 丝袜 人妻 在线| 啦啦啦在线免费观看视频4| 三上悠亚av全集在线观看| 亚洲精华国产精华精| 国产99久久九九免费精品| 国产深夜福利视频在线观看| 国产欧美日韩一区二区精品| 精品国产超薄肉色丝袜足j| 啪啪无遮挡十八禁网站| 免费在线观看视频国产中文字幕亚洲 | 爱豆传媒免费全集在线观看| 丁香六月天网| 国产97色在线日韩免费| 新久久久久国产一级毛片| 亚洲人成电影免费在线| 久久久久久人人人人人| 超色免费av| 曰老女人黄片| 亚洲av欧美aⅴ国产| 欧美人与性动交α欧美软件| 90打野战视频偷拍视频| 国产成人精品久久二区二区免费| 国产免费现黄频在线看| 热re99久久国产66热| 日韩,欧美,国产一区二区三区| 国产精品一区二区在线不卡| 欧美97在线视频| 桃花免费在线播放| 午夜免费成人在线视频| 国产成人av教育| 飞空精品影院首页| 免费高清在线观看视频在线观看| 侵犯人妻中文字幕一二三四区| 最近最新中文字幕大全免费视频| av不卡在线播放| 99国产综合亚洲精品| 久久 成人 亚洲| 欧美黄色片欧美黄色片| 国产精品一区二区在线不卡| 久久亚洲国产成人精品v| 日本猛色少妇xxxxx猛交久久| 亚洲国产毛片av蜜桃av| 久久久久久亚洲精品国产蜜桃av| 丁香六月天网| 国产成人欧美在线观看 | 三上悠亚av全集在线观看| 中文欧美无线码| 一二三四在线观看免费中文在| 亚洲av日韩在线播放| 高清欧美精品videossex| 久久国产精品男人的天堂亚洲| 男女床上黄色一级片免费看| 老汉色∧v一级毛片| 欧美av亚洲av综合av国产av| 久久人妻福利社区极品人妻图片| 欧美中文综合在线视频| 婷婷色av中文字幕| 亚洲av日韩精品久久久久久密| 久久久国产一区二区| 国产一级毛片在线| 日本a在线网址| 制服人妻中文乱码| 精品亚洲成国产av| 又大又爽又粗| 丰满少妇做爰视频| 国产91精品成人一区二区三区 | 欧美日韩一级在线毛片| a级毛片黄视频| 国产亚洲欧美在线一区二区| 日本猛色少妇xxxxx猛交久久| 亚洲熟女精品中文字幕| 淫妇啪啪啪对白视频 | 宅男免费午夜| 香蕉丝袜av| 9色porny在线观看| 欧美av亚洲av综合av国产av| 91国产中文字幕| 黑人操中国人逼视频| 国产免费av片在线观看野外av| bbb黄色大片| 精品福利观看| 少妇被粗大的猛进出69影院| 欧美日韩黄片免| 免费观看人在逋| 久久午夜综合久久蜜桃| 极品少妇高潮喷水抽搐| 国产精品成人在线| 国产精品久久久久久精品电影小说| av免费在线观看网站| 久久久水蜜桃国产精品网| 精品久久蜜臀av无| 妹子高潮喷水视频| 满18在线观看网站| av在线播放精品| 久久毛片免费看一区二区三区| 国产成人免费观看mmmm| 悠悠久久av| 人成视频在线观看免费观看| 99久久99久久久精品蜜桃| 成年动漫av网址| 国产国语露脸激情在线看| 欧美国产精品va在线观看不卡| 一本一本久久a久久精品综合妖精| 操美女的视频在线观看| 一区二区三区激情视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av电影在线观看一区二区三区| 人人妻,人人澡人人爽秒播| 免费人妻精品一区二区三区视频| 美国免费a级毛片| 精品少妇黑人巨大在线播放| 亚洲第一欧美日韩一区二区三区 | 国产男女内射视频| 国产精品影院久久| 亚洲美女黄色视频免费看| 欧美人与性动交α欧美软件| 精品人妻熟女毛片av久久网站| 日本猛色少妇xxxxx猛交久久| h视频一区二区三区| 法律面前人人平等表现在哪些方面 | 日本一区二区免费在线视频| 涩涩av久久男人的天堂| 亚洲国产欧美在线一区| 久久亚洲精品不卡| 久久人人爽av亚洲精品天堂| 女性生殖器流出的白浆| 亚洲精品国产av成人精品| 51午夜福利影视在线观看| 手机成人av网站| 免费看十八禁软件| 国产无遮挡羞羞视频在线观看| 制服诱惑二区| 亚洲av成人不卡在线观看播放网 | 亚洲视频免费观看视频| 久久影院123| 日本黄色日本黄色录像| 精品国产国语对白av| 国产亚洲精品久久久久5区| 超碰成人久久| 亚洲精品日韩在线中文字幕| 黄色 视频免费看| 热99re8久久精品国产| 在线观看人妻少妇| 国产精品欧美亚洲77777| 国产不卡av网站在线观看| 色视频在线一区二区三区| 国产精品九九99| 咕卡用的链子| 伊人亚洲综合成人网| 麻豆乱淫一区二区| 国产极品粉嫩免费观看在线| 欧美av亚洲av综合av国产av| www.熟女人妻精品国产| 久久久水蜜桃国产精品网| 王馨瑶露胸无遮挡在线观看| 久久精品国产亚洲av高清一级| 亚洲欧美精品综合一区二区三区| 波多野结衣一区麻豆| www.自偷自拍.com| 自拍欧美九色日韩亚洲蝌蚪91| 国产免费av片在线观看野外av| 黄片播放在线免费| 免费高清在线观看日韩| 国产精品久久久久久精品电影小说| 又大又爽又粗| 好男人电影高清在线观看| 18禁裸乳无遮挡动漫免费视频| 亚洲精品久久成人aⅴ小说| 两个人看的免费小视频| 精品久久久久久久毛片微露脸 | 黑丝袜美女国产一区| 欧美精品啪啪一区二区三区 | 亚洲午夜精品一区,二区,三区| 丝袜人妻中文字幕| 考比视频在线观看| 满18在线观看网站| 久久精品国产综合久久久| 日韩三级视频一区二区三区| 国产高清国产精品国产三级| 午夜免费观看性视频| 91麻豆av在线| 十八禁高潮呻吟视频| 午夜老司机福利片| 黄色 视频免费看| 青草久久国产| 久久人人爽人人片av| 精品福利永久在线观看| 欧美人与性动交α欧美软件| 久久国产精品大桥未久av| 亚洲 欧美一区二区三区| 亚洲精品一二三| 国产片内射在线| 欧美变态另类bdsm刘玥| 国产精品免费大片| 亚洲第一av免费看| 狂野欧美激情性xxxx| 亚洲国产毛片av蜜桃av| 国产97色在线日韩免费| 久久精品国产亚洲av香蕉五月 | 高清av免费在线| 国产日韩欧美在线精品| 亚洲成国产人片在线观看| 亚洲第一欧美日韩一区二区三区 | 亚洲avbb在线观看| 99久久人妻综合| 午夜老司机福利片| 91字幕亚洲| 欧美日韩一级在线毛片| 每晚都被弄得嗷嗷叫到高潮| 国产成人影院久久av| 秋霞在线观看毛片| 亚洲av国产av综合av卡| 国产视频一区二区在线看| 中国美女看黄片| 久热爱精品视频在线9| 最新的欧美精品一区二区| 亚洲欧美激情在线| 国产黄色免费在线视频| 精品少妇一区二区三区视频日本电影| 亚洲一区二区三区欧美精品| 黑人巨大精品欧美一区二区mp4| 久久久久国内视频| 亚洲三区欧美一区| 母亲3免费完整高清在线观看| 老司机在亚洲福利影院| 老鸭窝网址在线观看| 亚洲专区国产一区二区| 18在线观看网站| 一区二区三区激情视频| 亚洲美女黄色视频免费看| 丰满饥渴人妻一区二区三| 国产黄频视频在线观看| 午夜两性在线视频| 亚洲美女黄色视频免费看| 成人av一区二区三区在线看 | 乱人伦中国视频| 亚洲免费av在线视频| 亚洲欧美色中文字幕在线| av又黄又爽大尺度在线免费看| 日韩 欧美 亚洲 中文字幕| 一级毛片电影观看| 捣出白浆h1v1| 久久精品国产综合久久久| 天堂中文最新版在线下载| www日本在线高清视频| 1024视频免费在线观看| 老汉色∧v一级毛片| 免费在线观看黄色视频的| 日韩欧美国产一区二区入口| 高清av免费在线| 亚洲精品国产精品久久久不卡| www日本在线高清视频| 国产熟女午夜一区二区三区| 男人添女人高潮全过程视频| 国产av国产精品国产| 男女免费视频国产| 国产欧美日韩一区二区三区在线| 法律面前人人平等表现在哪些方面 | 国产男人的电影天堂91| 汤姆久久久久久久影院中文字幕| 99国产极品粉嫩在线观看| 国产人伦9x9x在线观看| 叶爱在线成人免费视频播放| 丰满饥渴人妻一区二区三| 一本色道久久久久久精品综合| 日本精品一区二区三区蜜桃| 一区二区三区乱码不卡18| 免费高清在线观看视频在线观看| 午夜福利视频在线观看免费| 久久人人爽av亚洲精品天堂| 性色av乱码一区二区三区2| 欧美成人午夜精品| 少妇的丰满在线观看| 日本黄色日本黄色录像| 精品福利观看| 日韩一卡2卡3卡4卡2021年| 久久中文看片网| 国产精品久久久久成人av| 叶爱在线成人免费视频播放| 曰老女人黄片| 久久九九热精品免费| 久久精品亚洲熟妇少妇任你| 极品少妇高潮喷水抽搐| tocl精华| 飞空精品影院首页| 视频区图区小说| 精品亚洲成国产av| 国产精品 欧美亚洲| 19禁男女啪啪无遮挡网站| 少妇猛男粗大的猛烈进出视频| 精品免费久久久久久久清纯 | av欧美777| 99久久国产精品久久久| 97在线人人人人妻| 成人国产av品久久久| 纯流量卡能插随身wifi吗| 国产精品亚洲av一区麻豆| 久热这里只有精品99| 欧美精品av麻豆av| 久久久水蜜桃国产精品网| 一本久久精品| 老鸭窝网址在线观看| 啦啦啦啦在线视频资源| 欧美日韩一级在线毛片| 天天躁夜夜躁狠狠躁躁| 99久久精品国产亚洲精品| 国产xxxxx性猛交| 国产亚洲午夜精品一区二区久久| 亚洲成人国产一区在线观看| 国精品久久久久久国模美| 91av网站免费观看| 国产精品免费视频内射| 波多野结衣一区麻豆| 精品人妻一区二区三区麻豆| 欧美日韩亚洲国产一区二区在线观看 | 亚洲国产精品一区二区三区在线| 一边摸一边做爽爽视频免费| 日本一区二区免费在线视频| tube8黄色片| 伦理电影免费视频| 国产亚洲欧美在线一区二区| 美女国产高潮福利片在线看| 一区二区三区精品91| 午夜免费成人在线视频| 亚洲欧洲日产国产| 蜜桃在线观看..| 中文字幕高清在线视频| 最近最新中文字幕大全免费视频| 亚洲中文字幕日韩| 视频区图区小说| 久久狼人影院| 精品国产乱子伦一区二区三区 | 国产精品久久久久成人av| 免费高清在线观看视频在线观看| 嫩草影视91久久|