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

    基于Qblade和Matlab的風(fēng)力機(jī)葉片設(shè)計(jì)與氣動性能分析*

    2016-12-25 02:04:25肖云峰高鵬遠(yuǎn)張志蓮周秀博
    化工機(jī)械 2016年4期
    關(guān)鍵詞:葉素利用系數(shù)風(fēng)力機(jī)

    肖云峰 高鵬遠(yuǎn) 張志蓮 呂 濤 周秀博

    (1. 北京石油化工學(xué)院機(jī)械工程學(xué)院; 2. 北京化工大學(xué)機(jī)電工程學(xué)院)

    基于Qblade和Matlab的風(fēng)力機(jī)葉片設(shè)計(jì)與氣動性能分析*

    肖云峰**1高鵬遠(yuǎn)2張志蓮1呂 濤1周秀博2

    (1. 北京石油化工學(xué)院機(jī)械工程學(xué)院; 2. 北京化工大學(xué)機(jī)電工程學(xué)院)

    首先使用Qblade計(jì)算翼型的氣動參數(shù),將數(shù)據(jù)進(jìn)行處理后,利用Wilson設(shè)計(jì)模型,結(jié)合Matlab編程軟件對某小型風(fēng)力機(jī)葉片進(jìn)行設(shè)計(jì)。再基于葉素動量理論,利用Qblade對所設(shè)計(jì)的葉片進(jìn)行氣動性能計(jì)算,并根據(jù)weibull風(fēng)速分布模型計(jì)算風(fēng)力機(jī)的性能。計(jì)算結(jié)果表明:編寫的程序正確,Qblade風(fēng)力機(jī)性能計(jì)算軟件能正確反映葉片氣動模型,并在保證計(jì)算結(jié)果準(zhǔn)確的同時節(jié)省風(fēng)力機(jī)葉片設(shè)計(jì)前期計(jì)算的時間。

    風(fēng)力機(jī) 葉片設(shè)計(jì) 葉素動量理論 Wilson模型 Qblade軟件 氣動性能計(jì)算

    目前在石油、煤炭等傳統(tǒng)能源短缺的環(huán)境下,世界各國均將目光投向新能源領(lǐng)域。近年來,風(fēng)能備受關(guān)注,我國已經(jīng)建立了多個大型風(fēng)電場[1]。風(fēng)力機(jī)葉片是風(fēng)力機(jī)將風(fēng)能轉(zhuǎn)化為機(jī)械能的核心部件,葉片參數(shù)直接影響風(fēng)力發(fā)電機(jī)的整體效率,因此葉片設(shè)計(jì)是風(fēng)力機(jī)的首要部分[2]。風(fēng)力機(jī)葉片設(shè)計(jì)分為氣動設(shè)計(jì)和結(jié)構(gòu)設(shè)計(jì)兩大部分[3],其中氣動設(shè)計(jì)包括氣動外形設(shè)計(jì)和氣動性能計(jì)算兩部分。目前,國內(nèi)外學(xué)者針對葉片外形設(shè)計(jì)一般采用葉片外形設(shè)計(jì)理論結(jié)合Matlab軟件的方法[4,5]。氣動性能計(jì)算為氣動設(shè)計(jì)結(jié)果提供評價和反饋,并為葉片的結(jié)構(gòu)設(shè)計(jì)提供氣動載荷等原始數(shù)據(jù)。氣動性能計(jì)算的準(zhǔn)確性關(guān)系到葉片的氣動性能和結(jié)構(gòu)安全。氣動性能計(jì)算的方法主要有:基于葉素動量理論建立氣動計(jì)算模型、基于渦流理論的渦尾跡方法和CFD方法[6]。其中葉素動量理論使用最多,葉素動量方法最主要的優(yōu)點(diǎn)在于既能保證計(jì)算結(jié)果準(zhǔn)確性,又能節(jié)省計(jì)算時間。國內(nèi)外學(xué)者在使用葉素動量方法進(jìn)行氣動計(jì)算時,大都通過編程的方式進(jìn)行,由于個人水平的不同,所編程序的準(zhǔn)確性和編程所花費(fèi)的時間也是不同的。而筆者所使用的基于葉素動量理論的包含葉片氣動性能計(jì)算、轉(zhuǎn)子葉片性能和風(fēng)力機(jī)性能計(jì)算的Qblade軟件在對風(fēng)力機(jī)葉片進(jìn)行評價的過程中能提高效率,因此筆者利用Matlab和Qblade軟件分別對葉片進(jìn)行外形設(shè)計(jì)和氣動性能計(jì)算。

    1 理論基礎(chǔ)

    1.1 葉素動量理論

    當(dāng)前國內(nèi)外計(jì)算風(fēng)力機(jī)葉片氣動性能的理論有貝茨理論、動量理論、葉素理論及葉素動量理論[7,8]等。其中,葉素動量理論結(jié)合了動量理論和葉素理論,能夠計(jì)算出風(fēng)輪掃掠面中的迭代變量軸向誘導(dǎo)因子a和周向誘導(dǎo)因子b。如圖1所示,風(fēng)輪流動模型簡化為一個理想的單元流管,并將它離散成N個高度為dr的環(huán)形單元,單元之間沒有流動。

    圖1 葉素動量理論單元流管模型

    由動量理論可得到作用在dr微段上的推力dT和扭矩dM:

    (1)

    dM=4πρV1b(1-a)r3dr

    (2)

    通過葉素理論可以得到作用在dr微段上的推力dT與扭矩dM:

    (3)

    (4)

    由于葉素理論和動量理論得出的推力和扭矩分別相等,并結(jié)合速度三角形可得:

    (5)

    (6)

    1.2 葉片設(shè)計(jì)理論

    翼型氣動特性決定了整個轉(zhuǎn)子葉片的氣動特性[9],在選取二維翼型之后,對葉片進(jìn)行優(yōu)化設(shè)計(jì),然后適當(dāng)修正設(shè)計(jì)結(jié)果,得出槳葉的弦長C和安裝角β分布。目前,常用的葉片設(shè)計(jì)模型有Glauert模型和Wilson模型,其理論模型相當(dāng)成熟。Wilson設(shè)計(jì)模型考慮了升阻比和葉尖損失對葉片最佳性能的影響,并且研究了風(fēng)輪在非設(shè)計(jì)工況下的性能。因此筆者選取Wilson方法。

    由于阻力對軸向和切向誘導(dǎo)因子影響較小,因此Wilson模型忽略了阻力的影響,但考慮了葉梢損失的影響,因此,軸向誘導(dǎo)因子a和周向誘導(dǎo)因子b的關(guān)系式分別為:

    (7)

    (8)

    葉尖損失系數(shù)F為:

    (9)

    (10)

    根據(jù)以上關(guān)系式可得考慮葉梢損失的能量方程:

    (11)

    局部風(fēng)能利用系數(shù)的計(jì)算式為:

    (12)

    在滿足能量方程的前提下,用迭代法計(jì)算誘導(dǎo)因子a、b,再根據(jù)式(12)求出dCp的最大值。將求得的a、b和相應(yīng)的F代入誘導(dǎo)因子關(guān)系式,可得弦長C和安裝角β的關(guān)系式:

    (13)

    (14)

    葉片設(shè)計(jì)流程如下:根據(jù)葉素理論,將葉片沿徑向平均分成若干等份,稱為葉素;針對每一個葉素,求解以式(9)為目標(biāo)函數(shù)、以式(10)為條件函數(shù)的最優(yōu)化問題,得出每一截面的a、b和F,進(jìn)而得出式(12)~(14)。

    1.3 風(fēng)速函數(shù)

    由于風(fēng)速變化范圍大,難以用準(zhǔn)確的數(shù)學(xué)模型加以描述,因此常用概率分布模型來表示,weibull函數(shù)分布形式簡單,并且能較好地模擬實(shí)際風(fēng)速分布,被認(rèn)為是風(fēng)能分析的有效工具。weibull分布模型的風(fēng)速概率密度分布函數(shù)為[10]:

    (15)

    式中Ch——尺寸參數(shù),m;

    K——形狀參數(shù),無量綱;

    V——來流速度,m/s。

    2 實(shí)例分析

    2.1 葉片參數(shù)設(shè)計(jì)

    筆者以額定功率為400W的小型風(fēng)力發(fā)電機(jī)為例,使用Wilson模型,利用Matlab進(jìn)行葉片外形設(shè)計(jì),隨后利用基于葉素動量理論的Qblade軟件計(jì)算葉片氣動性能和風(fēng)力機(jī)性能。具體參數(shù)如下:

    額定功率 400W

    葉尖速比 6

    葉片數(shù) 3

    額定風(fēng)速 7m/s

    風(fēng)能利用系數(shù) 0.45

    空氣密度 1.225

    風(fēng)力機(jī)發(fā)電效率 0.81

    風(fēng)輪直徑 2.5m

    2.2 翼型的選擇

    翼型的外形決定了翼型的氣動性能,綜合考慮選擇NACA4412翼型。通過Qblade軟件Xfoil部分計(jì)算雷諾數(shù)為50 000時翼型NACA4412的氣動性能,得出翼型設(shè)計(jì)點(diǎn)的氣流攻角θ=6.5°,對應(yīng)的最佳升阻比Cl/Cd=109.8173,此時升力系數(shù)Cl=1.162682,阻力系數(shù)Cd=0.010587。升阻比隨氣流攻角的變化曲線如圖2所示。

    圖2 升阻比隨氣流攻角的變化曲線

    2.3 優(yōu)化設(shè)計(jì)

    從葉尖到葉根部,將葉片平均分為9段,一共10個截面。使用Wilson模型,結(jié)合Matlab計(jì)算出每個截面的弦長C和安裝角β。具體截面參數(shù)如圖3所示。

    圖3 Qblade中葉片截面參數(shù)

    其中,1、2截面是為了滿足安裝要求設(shè)置的圓形截面,對截面進(jìn)行編輯后,再通過Qblade對葉片進(jìn)行線性優(yōu)化,以滿足加工、結(jié)構(gòu)等方面的要求。

    2.4 氣動性能計(jì)算

    氣動性能計(jì)算對于風(fēng)力發(fā)電機(jī)的初步設(shè)計(jì)具有重要的參考意義,對風(fēng)力性能具有重要表征意義的指標(biāo)分為有量綱指標(biāo)和無量綱指標(biāo)兩部分。其中,有量綱指標(biāo)主要包括:風(fēng)力機(jī)扭矩M、風(fēng)輪推力T。無量綱指標(biāo)包括:風(fēng)能利用系數(shù)Cp、扭矩系數(shù)CM和推力系數(shù)CT[11]。

    將計(jì)算出來的截面弦長和安裝角數(shù)據(jù)輸入到Qblade軟件的HAWT部分,建立葉片模型,并對所設(shè)計(jì)的葉片進(jìn)行氣動性能計(jì)算。

    2.4.1風(fēng)輪氣動性能計(jì)算

    由風(fēng)電機(jī)組葉片、風(fēng)輪的受載特性可知,葉輪氣動模型中氣動轉(zhuǎn)矩、氣動推力等參數(shù)隨風(fēng)速或葉尖速比的變化及時響應(yīng)。仿真結(jié)果須與理論分析結(jié)果相一致,氣動模型才能正確反映出葉片的氣動性能。

    風(fēng)能利用系數(shù)Cp決定了風(fēng)力機(jī)風(fēng)輪所能獲取能量的總量,即反映風(fēng)力發(fā)電機(jī)從自然風(fēng)中捕獲風(fēng)能程度的系數(shù)。從圖4中可以看出,在葉尖速比達(dá)到設(shè)計(jì)點(diǎn)前,風(fēng)能利用系數(shù)隨葉尖速比呈上升趨勢,在一定范圍內(nèi)風(fēng)能利用系數(shù)保持在較高位置,獲能效率高,當(dāng)葉尖速比再增大時,葉片的風(fēng)能利用系數(shù)呈下降趨勢,直至為0。此變化過程符合葉片實(shí)際運(yùn)行特征,且當(dāng)葉尖速比達(dá)到6時,風(fēng)能利用系數(shù)達(dá)到最大值0.485。實(shí)際生產(chǎn)中的風(fēng)機(jī)當(dāng)風(fēng)能利用系數(shù)達(dá)到0.460時即可加工[12]。

    圖4 風(fēng)能利用系數(shù)隨葉尖速比的變化曲線

    推力系數(shù)CT很大程度上影響了塔架的設(shè)計(jì),利用推力系數(shù)結(jié)合塔筒計(jì)算可以對風(fēng)力機(jī)氣動載荷進(jìn)行計(jì)算。扭矩系數(shù)CM決定了齒輪箱的尺寸和發(fā)電機(jī)的選型。推力系數(shù)和扭矩系數(shù)隨葉尖速比的變化曲線如圖5、6所示。

    圖5 推力系數(shù)隨葉尖速比的變化曲線

    圖6 扭矩系數(shù)隨葉尖速比的變化曲線

    定義風(fēng)速范圍為1~18m/s,旋轉(zhuǎn)速度范圍為200~400r/min。選取扭矩隨風(fēng)速的變化曲線,得出在額定風(fēng)速和額定轉(zhuǎn)速下扭矩的值。由圖7可知,在1~18m/s內(nèi)風(fēng)力機(jī)的扭矩隨風(fēng)速增大而增大,扭矩的計(jì)算結(jié)果可以為傳動鏈設(shè)計(jì)提供依據(jù)。圖8表明,在風(fēng)速1~18m/s范圍內(nèi),推力隨風(fēng)速的增大而增大,在達(dá)到12m/s時,增幅下降。

    圖7 扭矩隨風(fēng)速的變化曲線

    圖8 推力隨風(fēng)速的變化曲線

    2.4.2風(fēng)力機(jī)性能計(jì)算

    整機(jī)的推力隨風(fēng)速的變化趨勢與圖8近似,扭矩隨風(fēng)速的變化趨勢與圖7一致,具體值相差不大,沒有出現(xiàn)突變,表明風(fēng)力機(jī)葉片設(shè)計(jì)合理,具有實(shí)際工程應(yīng)用價值。從風(fēng)能利用系數(shù)隨風(fēng)速的變化曲線(圖9)可知,該風(fēng)力機(jī)能有效捕捉來流的風(fēng)能,與最初設(shè)計(jì)目標(biāo)吻合。圖10表明,在3~18m/s范圍內(nèi)風(fēng)力機(jī)推力系數(shù)隨風(fēng)速的增大而減小,與風(fēng)輪氣動性能計(jì)算結(jié)果相一致。從扭矩系數(shù)隨風(fēng)速的變化曲線(圖11)可知扭矩系數(shù)逐漸增大直到風(fēng)速為10m/s時開始降低,因此來流風(fēng)速10m/s附近是最佳捕風(fēng)范圍。

    圖9 風(fēng)能利用系數(shù)隨風(fēng)速的變化曲線

    圖10 推力系數(shù)隨風(fēng)速的變化曲線

    圖11 扭矩系數(shù)隨風(fēng)速的變化曲線

    3 結(jié)論

    3.1基于Wilson模型,400W小型風(fēng)力機(jī)采用NACA4412翼型時,葉片設(shè)計(jì)合理,并為之后使用Qblade軟件進(jìn)行氣動性能分析提供了基礎(chǔ)。

    3.2利用Qblade軟件對所設(shè)計(jì)的葉片進(jìn)行氣動性能計(jì)算,避免了氣動計(jì)算復(fù)雜的編程過程和反復(fù)迭代循環(huán)計(jì)算,計(jì)算結(jié)果能準(zhǔn)確表現(xiàn)葉片氣動模型,氣動性能一維計(jì)算結(jié)果符合氣動規(guī)律,驗(yàn)證了設(shè)計(jì)葉片的合理性。筆者的設(shè)計(jì)和計(jì)算體現(xiàn)了Qblade軟件具有較好的工程實(shí)用性。

    3.3利用Qblade軟件對所設(shè)計(jì)的葉片進(jìn)行了氣動性能計(jì)算,結(jié)果表明,Qblade軟件在計(jì)算時間和精度上有很好的平衡,工程應(yīng)用價值強(qiáng),并為之后的結(jié)構(gòu)設(shè)計(jì)和模擬分析奠定了基礎(chǔ)。

    3.4通過風(fēng)輪氣動計(jì)算和風(fēng)力機(jī)性能仿真結(jié)果的比較,本次400W風(fēng)力機(jī)氣動設(shè)計(jì)合理,該葉型可直接投放生產(chǎn)。

    [1] 巴麗合亞,陳華.風(fēng)光互補(bǔ)發(fā)電技術(shù)在新疆的應(yīng)用及展望[J].化工自動化及儀表,2015:42(1):1~2.

    [2] 張禮達(dá),任臘春,陳榮盛,等.風(fēng)力機(jī)葉片外形設(shè)計(jì)及三維實(shí)體建模研究[J].太陽能學(xué)報(bào),2008,29(9):1177~1180.

    [3] 陳云程.風(fēng)力機(jī)設(shè)計(jì)與應(yīng)用[M].上海:上海科學(xué)技術(shù)出版社,1998.

    [4] 王軍,周丙超.基于Matlab的小型風(fēng)力機(jī)葉片設(shè)計(jì)[J].水電能源科學(xué),2007,25(5):142~144.

    [5] 宋文龍,劉湘晨,周曉.基于Matlab和Solidworks的水輪機(jī)葉片設(shè)計(jì)和三維建模[J].化工機(jī)械,2015,42(1):93~96.

    [6] 張仲柱.水平軸風(fēng)力機(jī)葉片氣動性能計(jì)算模型研究[D].北京:中國科學(xué)院研究生院(工程熱物理研究所),2007.

    [7] Noda M,F(xiàn)lay R G J.A Simulation Model for Wind Turbine Blade Fatigue Loads[J].Journal of Wind Engineering and Industrial Aerodynamics,1999,83(1/2/3):527~540.

    [8] Maalawi K Y,Badawy M T S.A Direct Method for Evaluating Performance of Horizontal Axis Wind Turbines[J]. Renewable & Sustainable Energy Reviews,2001,5(2):175~190.

    [9] 李春,葉舟,高偉,等.現(xiàn)代大型風(fēng)力機(jī)設(shè)計(jì)原理[M].上海:上海科學(xué)技術(shù)出版社,2013.

    [10] Weisser D.A Wind Energy Analysis of Grenada:an Estimation Using the ‘Weibull’Bensity Function [J].Renewable Energy,2003,28(11):1803~1812.

    [11] 姚興佳.風(fēng)力發(fā)電機(jī)組理論與設(shè)計(jì)[M].北京:機(jī)械工業(yè)出版社,2012.

    [12] 張湘東,聶國華.大型水平軸風(fēng)力機(jī)葉片氣動性能優(yōu)化[J].計(jì)算機(jī)輔助工程,2009,18(1):47~50.

    [13] 楊新彥.小型水平軸風(fēng)力發(fā)電機(jī)葉片設(shè)計(jì)研究[D].南寧:廣西大學(xué),2006.

    WindTurbineBladeDesignandAerodynamicPerformanceAnalysisBasedonQbladeandMatlab

    XIAO Yun-feng1, GAO Peng-yuan2, ZHANG Zhi-lian1,LV Tao1,ZHOU Xiu-bo2

    (1.CollegeofMechanicalEngineering,BeijingInstituteofPetrochemicalTechnology,Beijing102617,China;
    2.CollegeofMechanicalEngineering,BeijingUniversityofChemicalTechnology,Beijing100029,China)

    Qblade was used to calculate airfoil’s aerodynamic parameters; having Wilson design model adopted and Matlab programming software considered to design small wind turbine blades after processing the data concerned. Basing on the blade element momentum theory, Qblade was employed to calculate aerodynamic performance of the blades designed, including having Weibull distribution based to calculate wind turbine performance. The calculation results prove correctness of the program compiled and Qblade calculation software’s performance in correctly reflecting blade aerodynamic model; the better balance between the time and accuracy of the calculation can ensure right calculation results and save the time in pre-design phase.

    wind turbine, blade design, blade element momentum theory, Wilson model, Qblade software, aerodynamic performance calculation

    *北京石油化工學(xué)院優(yōu)秀青年教師培育計(jì)劃項(xiàng)目(08031862008/038)。

    **肖云峰,男,1976年10月生,副教授。北京市,102617。

    TQ050.2

    A

    0254-6094(2016)04-0467-05

    2015-06-15,

    2016-07-19)

    猜你喜歡
    葉素利用系數(shù)風(fēng)力機(jī)
    燒結(jié)漏風(fēng)率對利用系數(shù)影響關(guān)系的探討
    昆鋼科技(2020年4期)2020-10-23 09:31:58
    基于UIOs的風(fēng)力機(jī)傳動系統(tǒng)多故障診斷
    梅鋼提高三號燒結(jié)利用系數(shù)的實(shí)踐
    四川冶金(2018年1期)2018-09-25 02:39:28
    灌區(qū)渠系水利用系數(shù)推求方法優(yōu)選研究
    不能說的秘密
    精心設(shè)計(jì)老婆出軌釀滅門慘案
    新傳奇(2016年29期)2016-09-18 08:30:12
    那個不能說的秘密
    中外文摘(2016年12期)2016-08-29 08:14:40
    瘋狂的秘密:離婚5天后滅門岳父家
    大型風(fēng)力機(jī)整機(jī)氣動彈性響應(yīng)計(jì)算
    小型風(fēng)力機(jī)葉片快速建模方法
    太陽能(2015年6期)2015-02-28 17:09:35
    亚洲久久久久久中文字幕| 12—13女人毛片做爰片一| 狂野欧美激情性xxxx| 夜夜爽天天搞| 午夜a级毛片| av中文乱码字幕在线| 亚洲精品一卡2卡三卡4卡5卡| 熟女少妇亚洲综合色aaa.| 日日干狠狠操夜夜爽| 久久中文看片网| 国产午夜福利久久久久久| 国产色爽女视频免费观看| 国产欧美日韩一区二区精品| 免费高清视频大片| 亚洲电影在线观看av| 99久久成人亚洲精品观看| 久久久国产成人免费| 国产亚洲欧美在线一区二区| 亚洲精品亚洲一区二区| 一个人看视频在线观看www免费 | 又黄又爽又免费观看的视频| 欧美大码av| 精品一区二区三区视频在线 | 夜夜夜夜夜久久久久| 欧美一区二区国产精品久久精品| 亚洲真实伦在线观看| 少妇高潮的动态图| 久久国产精品人妻蜜桃| 国产成人av教育| 俺也久久电影网| 99国产精品一区二区蜜桃av| 床上黄色一级片| 精华霜和精华液先用哪个| 老汉色∧v一级毛片| 熟妇人妻久久中文字幕3abv| 天堂影院成人在线观看| 少妇人妻精品综合一区二区 | 香蕉丝袜av| 国产亚洲精品久久久久久毛片| 两人在一起打扑克的视频| 宅男免费午夜| 一个人看视频在线观看www免费 | 国产精品野战在线观看| 丝袜美腿在线中文| 国产精品久久久久久亚洲av鲁大| 中文字幕人妻熟人妻熟丝袜美 | 99久国产av精品| 亚洲av一区综合| 国产三级中文精品| 人妻丰满熟妇av一区二区三区| 国模一区二区三区四区视频| 欧美黑人欧美精品刺激| 国产91精品成人一区二区三区| 一区福利在线观看| 美女免费视频网站| 人人妻,人人澡人人爽秒播| 三级国产精品欧美在线观看| 黄色日韩在线| 国产毛片a区久久久久| 国产欧美日韩精品一区二区| or卡值多少钱| 国产精品av视频在线免费观看| avwww免费| 国产精品,欧美在线| 欧美不卡视频在线免费观看| 精品久久久久久久人妻蜜臀av| 久久国产精品人妻蜜桃| 精品人妻一区二区三区麻豆 | 国产伦精品一区二区三区视频9 | av天堂在线播放| 国产一区二区三区视频了| 国产精品野战在线观看| 亚洲人成网站在线播| 99久国产av精品| 两人在一起打扑克的视频| 俺也久久电影网| 琪琪午夜伦伦电影理论片6080| 国产av一区在线观看免费| 久久久久九九精品影院| 99久久无色码亚洲精品果冻| 国产欧美日韩精品亚洲av| 色综合亚洲欧美另类图片| 久久欧美精品欧美久久欧美| 中文亚洲av片在线观看爽| 亚洲成人免费电影在线观看| 久久精品亚洲精品国产色婷小说| 亚洲av中文字字幕乱码综合| 精品午夜福利视频在线观看一区| 午夜免费激情av| 国产精品久久视频播放| 亚洲国产精品成人综合色| 性色av乱码一区二区三区2| 久久久久久久午夜电影| 男女床上黄色一级片免费看| 91av网一区二区| 又粗又爽又猛毛片免费看| 欧美乱色亚洲激情| 麻豆一二三区av精品| 久久人妻av系列| 国产精品 欧美亚洲| 观看免费一级毛片| 国内久久婷婷六月综合欲色啪| 又爽又黄无遮挡网站| 亚洲av免费高清在线观看| 在线观看午夜福利视频| 色精品久久人妻99蜜桃| 国产精品女同一区二区软件 | 亚洲av免费高清在线观看| 精品久久久久久久末码| 亚洲aⅴ乱码一区二区在线播放| 国产黄色小视频在线观看| 国产伦人伦偷精品视频| 一个人免费在线观看的高清视频| 男女床上黄色一级片免费看| 热99在线观看视频| 国产美女午夜福利| 在线观看日韩欧美| 亚洲成人久久爱视频| 久久久国产成人免费| 久久久色成人| 99在线视频只有这里精品首页| 亚洲最大成人手机在线| 毛片女人毛片| tocl精华| 午夜精品久久久久久毛片777| 最近在线观看免费完整版| 真人一进一出gif抽搐免费| 久99久视频精品免费| 亚洲人成电影免费在线| 精品熟女少妇八av免费久了| 国产美女午夜福利| 在线免费观看不下载黄p国产 | 成年版毛片免费区| 亚洲av五月六月丁香网| h日本视频在线播放| 又黄又爽又免费观看的视频| 老司机深夜福利视频在线观看| 亚洲一区二区三区不卡视频| 亚洲最大成人中文| 免费看光身美女| 老司机在亚洲福利影院| 国产视频一区二区在线看| 在线天堂最新版资源| h日本视频在线播放| 欧美日韩综合久久久久久 | 国产熟女xx| 欧美高清成人免费视频www| 91九色精品人成在线观看| 人人妻人人看人人澡| av专区在线播放| 国产熟女xx| 老熟妇乱子伦视频在线观看| 婷婷精品国产亚洲av| 九九在线视频观看精品| 亚洲一区二区三区色噜噜| 精品国产亚洲在线| 好男人电影高清在线观看| 不卡一级毛片| 亚洲成人免费电影在线观看| 变态另类成人亚洲欧美熟女| 国产一区二区三区视频了| 又粗又爽又猛毛片免费看| 国产乱人视频| av中文乱码字幕在线| 国产一区在线观看成人免费| 国产蜜桃级精品一区二区三区| 午夜免费观看网址| 哪里可以看免费的av片| 深爱激情五月婷婷| 午夜福利在线观看吧| 日本免费a在线| 国产精品女同一区二区软件 | 中文字幕久久专区| 人人妻人人看人人澡| 婷婷亚洲欧美| 麻豆久久精品国产亚洲av| 波多野结衣高清无吗| 婷婷六月久久综合丁香| 十八禁网站免费在线| 免费观看人在逋| 欧美zozozo另类| 高清毛片免费观看视频网站| 可以在线观看的亚洲视频| 欧美色欧美亚洲另类二区| 少妇的丰满在线观看| 国产色爽女视频免费观看| 一级毛片高清免费大全| 国产黄片美女视频| 日本免费一区二区三区高清不卡| 在线视频色国产色| 国产私拍福利视频在线观看| АⅤ资源中文在线天堂| 国产精品自产拍在线观看55亚洲| 69人妻影院| 黄片大片在线免费观看| 91麻豆av在线| 精品久久久久久久久久免费视频| 一级毛片女人18水好多| 久久久久久久久中文| 两性午夜刺激爽爽歪歪视频在线观看| 波野结衣二区三区在线 | 亚洲专区中文字幕在线| 日本五十路高清| 草草在线视频免费看| 美女黄网站色视频| 亚洲狠狠婷婷综合久久图片| 中国美女看黄片| 免费在线观看成人毛片| 国产久久久一区二区三区| 啦啦啦韩国在线观看视频| 一个人看视频在线观看www免费 | 国产成人福利小说| 亚洲精品久久国产高清桃花| 男人的好看免费观看在线视频| 国产真实乱freesex| 久久久国产成人精品二区| 亚洲在线观看片| 成人高潮视频无遮挡免费网站| 亚洲国产高清在线一区二区三| 国产三级中文精品| 久久这里只有精品中国| 亚洲18禁久久av| 18美女黄网站色大片免费观看| 欧美日韩一级在线毛片| 性欧美人与动物交配| 国产 一区 欧美 日韩| 欧美乱色亚洲激情| 此物有八面人人有两片| 日本精品一区二区三区蜜桃| 成人午夜高清在线视频| 中出人妻视频一区二区| 18+在线观看网站| 麻豆成人av在线观看| 午夜影院日韩av| 变态另类丝袜制服| 最近最新免费中文字幕在线| 很黄的视频免费| 欧美国产日韩亚洲一区| 久久久久久久久中文| 久久久久久九九精品二区国产| 久久这里只有精品中国| 欧美一区二区国产精品久久精品| 免费看十八禁软件| 精品无人区乱码1区二区| 精品久久久久久久久久久久久| 亚洲成av人片在线播放无| 亚洲中文字幕日韩| 精品国产美女av久久久久小说| АⅤ资源中文在线天堂| 色精品久久人妻99蜜桃| 美女 人体艺术 gogo| 亚洲天堂国产精品一区在线| 国产熟女xx| 人人妻,人人澡人人爽秒播| 免费高清视频大片| 女人十人毛片免费观看3o分钟| 中文字幕精品亚洲无线码一区| 国产又黄又爽又无遮挡在线| 又黄又粗又硬又大视频| 国产免费一级a男人的天堂| 日韩 欧美 亚洲 中文字幕| 亚洲精品国产精品久久久不卡| 亚洲五月天丁香| 午夜激情福利司机影院| 精品国内亚洲2022精品成人| 日本成人三级电影网站| 日本与韩国留学比较| 国产精品日韩av在线免费观看| 天天躁日日操中文字幕| 国产单亲对白刺激| 欧美高清成人免费视频www| 国内少妇人妻偷人精品xxx网站| 88av欧美| 日本一二三区视频观看| 精品久久久久久久久久免费视频| 国产精品久久久久久久电影 | bbb黄色大片| 精品国内亚洲2022精品成人| 两个人的视频大全免费| 精品久久久久久久毛片微露脸| 18禁美女被吸乳视频| 免费在线观看成人毛片| 国产伦精品一区二区三区四那| 欧美日韩黄片免| 亚洲国产欧洲综合997久久,| 18禁裸乳无遮挡免费网站照片| 午夜福利成人在线免费观看| 亚洲一区二区三区不卡视频| 国产精品久久久久久亚洲av鲁大| 亚洲精品美女久久久久99蜜臀| 俺也久久电影网| 最近最新中文字幕大全电影3| 在线天堂最新版资源| 免费av毛片视频| 欧美成人a在线观看| 成人18禁在线播放| 免费av不卡在线播放| 男女下面进入的视频免费午夜| а√天堂www在线а√下载| 国产精品久久久久久久久免 | e午夜精品久久久久久久| 欧美日韩国产亚洲二区| 97碰自拍视频| 成人特级av手机在线观看| av片东京热男人的天堂| 18+在线观看网站| 午夜影院日韩av| 性色av乱码一区二区三区2| 少妇的丰满在线观看| 国产精品国产高清国产av| 啦啦啦免费观看视频1| 搡老妇女老女人老熟妇| 久久久国产成人精品二区| 国产免费男女视频| 精华霜和精华液先用哪个| 欧美日韩乱码在线| 亚洲一区高清亚洲精品| 久久久久免费精品人妻一区二区| 人人妻人人澡欧美一区二区| 18禁在线播放成人免费| 美女高潮的动态| 国产探花极品一区二区| 国产男靠女视频免费网站| 舔av片在线| 久久久国产成人精品二区| 久久精品国产亚洲av涩爱 | 亚洲第一欧美日韩一区二区三区| 欧美绝顶高潮抽搐喷水| av天堂在线播放| 舔av片在线| 桃色一区二区三区在线观看| 日日摸夜夜添夜夜添小说| 久久久久久九九精品二区国产| 亚洲国产欧美人成| 99久久精品国产亚洲精品| 国产成+人综合+亚洲专区| 美女免费视频网站| 亚洲aⅴ乱码一区二区在线播放| 亚洲avbb在线观看| h日本视频在线播放| 国产伦精品一区二区三区视频9 | 亚洲avbb在线观看| 欧美最新免费一区二区三区 | 国产精品一区二区三区四区免费观看 | 老司机福利观看| 亚洲无线观看免费| www日本黄色视频网| 亚洲国产精品成人综合色| 国产久久久一区二区三区| 熟女电影av网| www.999成人在线观看| av女优亚洲男人天堂| 中文字幕av在线有码专区| 少妇人妻一区二区三区视频| 欧美黑人巨大hd| 午夜久久久久精精品| 亚洲中文字幕一区二区三区有码在线看| 黄片小视频在线播放| 一个人观看的视频www高清免费观看| 免费搜索国产男女视频| 免费大片18禁| 日本 av在线| 此物有八面人人有两片| 中文字幕av在线有码专区| 国产综合懂色| 国产亚洲精品一区二区www| 女同久久另类99精品国产91| 少妇人妻一区二区三区视频| 国产精品精品国产色婷婷| 亚洲人与动物交配视频| 看片在线看免费视频| 人妻夜夜爽99麻豆av| 五月玫瑰六月丁香| 久99久视频精品免费| 国产欧美日韩精品一区二区| 亚洲中文字幕日韩| 欧美黑人欧美精品刺激| 免费观看的影片在线观看| 精品人妻一区二区三区麻豆 | 男女床上黄色一级片免费看| a在线观看视频网站| 午夜久久久久精精品| 黄色成人免费大全| 久久精品国产亚洲av香蕉五月| 99久久无色码亚洲精品果冻| 特大巨黑吊av在线直播| 亚洲最大成人中文| 在线观看一区二区三区| 91久久精品国产一区二区成人 | 日韩欧美国产在线观看| 成人特级黄色片久久久久久久| 天堂动漫精品| 99在线视频只有这里精品首页| a级一级毛片免费在线观看| 亚洲成人精品中文字幕电影| 亚洲欧美一区二区三区黑人| 免费看日本二区| 欧美不卡视频在线免费观看| 国产欧美日韩一区二区三| 久9热在线精品视频| 久久久久免费精品人妻一区二区| eeuss影院久久| 俺也久久电影网| 久久久久久久久中文| 中文字幕人成人乱码亚洲影| 性色avwww在线观看| 免费一级毛片在线播放高清视频| www日本黄色视频网| 欧美成狂野欧美在线观看| 99久久99久久久精品蜜桃| 日韩成人在线观看一区二区三区| 91久久精品国产一区二区成人 | 国产精品98久久久久久宅男小说| 99精品欧美一区二区三区四区| 亚洲欧美激情综合另类| 一个人看的www免费观看视频| 欧美黑人欧美精品刺激| 欧美色欧美亚洲另类二区| 亚洲精品一卡2卡三卡4卡5卡| 天堂动漫精品| 99久国产av精品| 色哟哟哟哟哟哟| 免费大片18禁| 国产在视频线在精品| 成人av一区二区三区在线看| 少妇裸体淫交视频免费看高清| 欧美国产日韩亚洲一区| 一进一出抽搐gif免费好疼| 极品教师在线免费播放| 级片在线观看| 亚洲精品亚洲一区二区| 国产高清有码在线观看视频| 亚洲人成网站高清观看| 国产男靠女视频免费网站| 最后的刺客免费高清国语| 又紧又爽又黄一区二区| 国产午夜福利久久久久久| 精品久久久久久久毛片微露脸| 久久伊人香网站| 伊人久久大香线蕉亚洲五| 成人高潮视频无遮挡免费网站| 久久99热这里只有精品18| 免费观看的影片在线观看| 国产免费男女视频| 一二三四社区在线视频社区8| 婷婷丁香在线五月| 国产一级毛片七仙女欲春2| 女人十人毛片免费观看3o分钟| 免费观看精品视频网站| 3wmmmm亚洲av在线观看| 99精品久久久久人妻精品| 内射极品少妇av片p| 99精品欧美一区二区三区四区| 少妇的丰满在线观看| 亚洲 国产 在线| 色播亚洲综合网| 男女那种视频在线观看| 无遮挡黄片免费观看| 嫁个100分男人电影在线观看| 国产精品亚洲美女久久久| 深夜精品福利| 一卡2卡三卡四卡精品乱码亚洲| ponron亚洲| 嫩草影院入口| 国产亚洲精品综合一区在线观看| 免费电影在线观看免费观看| 美女高潮的动态| 悠悠久久av| 亚洲欧美日韩卡通动漫| 成人精品一区二区免费| 久久久成人免费电影| 两人在一起打扑克的视频| av在线蜜桃| 天天躁日日操中文字幕| 一本综合久久免费| 一进一出好大好爽视频| e午夜精品久久久久久久| 国产视频一区二区在线看| 乱人视频在线观看| 亚洲欧美日韩高清在线视频| av欧美777| 桃色一区二区三区在线观看| 国产一区二区三区视频了| а√天堂www在线а√下载| 看片在线看免费视频| 99久久精品热视频| 99久久精品国产亚洲精品| 老汉色∧v一级毛片| 亚洲精品影视一区二区三区av| 国产亚洲精品av在线| 99在线人妻在线中文字幕| 中文字幕人妻熟人妻熟丝袜美 | 精品一区二区三区视频在线 | svipshipincom国产片| 亚洲av五月六月丁香网| 国产精品爽爽va在线观看网站| 国产成人啪精品午夜网站| 神马国产精品三级电影在线观看| 毛片女人毛片| 黄色成人免费大全| 偷拍熟女少妇极品色| 精品一区二区三区av网在线观看| 男人舔女人下体高潮全视频| 精品无人区乱码1区二区| 国产伦精品一区二区三区四那| 久久久久亚洲av毛片大全| 母亲3免费完整高清在线观看| 中文字幕久久专区| 久久婷婷人人爽人人干人人爱| 亚洲一区高清亚洲精品| 亚洲成人久久性| 欧美激情在线99| 欧美乱码精品一区二区三区| 99精品欧美一区二区三区四区| 亚洲不卡免费看| 亚洲真实伦在线观看| 成年人黄色毛片网站| 国产精品98久久久久久宅男小说| 琪琪午夜伦伦电影理论片6080| 国产在视频线在精品| 又粗又爽又猛毛片免费看| 高清毛片免费观看视频网站| 亚洲av五月六月丁香网| 欧美一区二区亚洲| 真人做人爱边吃奶动态| 91麻豆精品激情在线观看国产| 校园春色视频在线观看| 国产精品香港三级国产av潘金莲| 狂野欧美激情性xxxx| 国产野战对白在线观看| 精品熟女少妇八av免费久了| 亚洲人成网站在线播放欧美日韩| 亚洲狠狠婷婷综合久久图片| 免费人成在线观看视频色| 久久久久久久久大av| 欧美一级毛片孕妇| 午夜福利18| 精品无人区乱码1区二区| 成人av一区二区三区在线看| 精品久久久久久久末码| 欧美不卡视频在线免费观看| 亚洲av第一区精品v没综合| 看黄色毛片网站| 国产av不卡久久| 亚洲天堂国产精品一区在线| 亚洲美女视频黄频| 高清在线国产一区| eeuss影院久久| 久久草成人影院| 日本五十路高清| 少妇人妻精品综合一区二区 | 欧美一区二区精品小视频在线| 成年人黄色毛片网站| 亚洲电影在线观看av| 免费搜索国产男女视频| 亚洲欧美日韩东京热| 性色av乱码一区二区三区2| 欧美高清成人免费视频www| 五月玫瑰六月丁香| 日韩大尺度精品在线看网址| 99精品欧美一区二区三区四区| 国产在视频线在精品| 国产不卡一卡二| 蜜桃亚洲精品一区二区三区| 无遮挡黄片免费观看| 午夜激情欧美在线| 人人妻,人人澡人人爽秒播| 日本撒尿小便嘘嘘汇集6| 又粗又爽又猛毛片免费看| 男人的好看免费观看在线视频| 全区人妻精品视频| 午夜福利在线观看吧| 久久久久国产精品人妻aⅴ院| а√天堂www在线а√下载| 国产伦在线观看视频一区| 久久精品综合一区二区三区| 高潮久久久久久久久久久不卡| 宅男免费午夜| 精品一区二区三区av网在线观看| 88av欧美| 久久久精品大字幕| 中文字幕人妻丝袜一区二区| 国产成人a区在线观看| 国产精品一区二区三区四区久久| 国产精品一区二区免费欧美| 久久久久久久久大av| 天堂影院成人在线观看| 国产精品久久久人人做人人爽| 丰满人妻熟妇乱又伦精品不卡| 久久99热这里只有精品18| 精品日产1卡2卡| 日本 av在线| 精品人妻一区二区三区麻豆 | 欧美日韩精品网址| 色播亚洲综合网| 免费av不卡在线播放| 日韩欧美精品免费久久 | 成人av一区二区三区在线看| 国产真实乱freesex| 日本与韩国留学比较| 蜜桃久久精品国产亚洲av| 丰满人妻一区二区三区视频av | 国产三级中文精品| 国产精品久久久久久久久免 | 亚洲真实伦在线观看| 亚洲久久久久久中文字幕| 天堂动漫精品| 亚洲一区二区三区不卡视频| 欧美成人a在线观看| 成人av一区二区三区在线看| 国产97色在线日韩免费| 国产高清激情床上av| 国产色婷婷99| 亚洲av一区综合| 变态另类成人亚洲欧美熟女|