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

    葉片尾緣加彎板垂直軸風(fēng)力機(jī)轉(zhuǎn)矩特性數(shù)值模擬

    2017-04-01 06:10:26李巖張婷婷田川公太朗馮放鄭玉芳公維佳
    關(guān)鍵詞:垂直軸尾緣風(fēng)輪

    李巖,張婷婷,田川公太朗,馮放,鄭玉芳,公維佳

    (1.東北農(nóng)業(yè)大學(xué)工程學(xué)院,哈爾濱 150030;2.寒地農(nóng)業(yè)可再生資源利用技術(shù)與裝備黑龍江省重點(diǎn)實(shí)驗(yàn)室,哈爾濱 150030;3.日本鳥取大學(xué)地域?qū)W部,日本鳥取 680855 2;4.東北農(nóng)業(yè)大學(xué)理學(xué)院,哈爾濱 150030)

    葉片尾緣加彎板垂直軸風(fēng)力機(jī)轉(zhuǎn)矩特性數(shù)值模擬

    李巖1,2,張婷婷1,田川公太朗3,馮放4,鄭玉芳1,公維佳1,2

    (1.東北農(nóng)業(yè)大學(xué)工程學(xué)院,哈爾濱 150030;2.寒地農(nóng)業(yè)可再生資源利用技術(shù)與裝備黑龍江省重點(diǎn)實(shí)驗(yàn)室,哈爾濱 150030;3.日本鳥取大學(xué)地域?qū)W部,日本鳥取 680855 2;4.東北農(nóng)業(yè)大學(xué)理學(xué)院,哈爾濱 150030)

    文章提出一種在葉片尾緣加裝彎板改善直線翼垂直軸風(fēng)力機(jī)葉片周圍流場,提高風(fēng)力機(jī)轉(zhuǎn)矩特性的方法。選取NACA0024和NACA0012兩種直線翼垂直軸風(fēng)力機(jī)常用翼型,利用數(shù)值模擬計(jì)算包括不加裝彎板在內(nèi)具有7種彎板長度2葉片直線翼垂直軸風(fēng)力機(jī)輸出轉(zhuǎn)矩特性和靜態(tài)起動轉(zhuǎn)矩特性。結(jié)果表明,在NACA系列對稱翼型尾緣加彎板葉片可在一定程度上改善直線翼垂直軸風(fēng)力機(jī)輸出轉(zhuǎn)矩特性和靜態(tài)起動特性。葉片尾緣加裝彎板后可在某些旋轉(zhuǎn)角下改變?nèi)~片上下表面壓力分布,減少渦旋和流動分離產(chǎn)生,提高葉片氣動特性,改善風(fēng)力機(jī)轉(zhuǎn)矩特性。研究發(fā)現(xiàn),彎板長度為弦長40%風(fēng)力機(jī)轉(zhuǎn)矩特性改善效果最好,最大輸出轉(zhuǎn)矩系數(shù)較無彎板風(fēng)力機(jī)提高17%,靜態(tài)平均起動轉(zhuǎn)矩提高26.4%。

    垂直軸風(fēng)力機(jī);葉片尾緣;彎板;轉(zhuǎn)矩特性;數(shù)值模擬

    與水平軸風(fēng)力機(jī)相比,垂直軸風(fēng)力機(jī)最大特征為無需迎風(fēng)轉(zhuǎn)向裝置,應(yīng)用于中小型風(fēng)能利用系統(tǒng)、分布式發(fā)電和微網(wǎng)[1]。近年來,作為垂直軸風(fēng)力機(jī)代表直線翼垂直軸風(fēng)力機(jī)以形式簡單、容易加工、結(jié)構(gòu)新穎等特點(diǎn)受到關(guān)注。然而,該風(fēng)力機(jī)旋轉(zhuǎn)形成風(fēng)輪空間內(nèi)部流動復(fù)雜,為典型大分離非定常流動[2],因此其氣動特性分析與設(shè)計(jì)相對難于中小型水平軸風(fēng)力機(jī)。隨著風(fēng)力機(jī)轉(zhuǎn)速提高,風(fēng)輪內(nèi)外流動愈加復(fù)雜,伴隨葉片邊界層分離和大量渦旋形成,風(fēng)力機(jī)在某些轉(zhuǎn)速下動態(tài)失速,影響風(fēng)力機(jī)葉片氣動力作用在風(fēng)輪中心轉(zhuǎn)矩特性,風(fēng)力機(jī)整體輸出特性降低,為目前該種風(fēng)力機(jī)發(fā)展制約因素。因此,要有效控制分離與渦旋產(chǎn)生,研究適合直線翼垂直軸風(fēng)力機(jī)葉片與翼型尤為重要[3]。為此,研究者從不同角度開展相關(guān)研究。廖書學(xué)等分析翼型厚度對垂直軸風(fēng)力機(jī)氣動性能影響,結(jié)果表明翼型厚度對風(fēng)力機(jī)各氣動參數(shù)影響較大[4]。趙萬里等對增設(shè)Gurney襟翼風(fēng)力機(jī)專用翼型作數(shù)值研究,分析Gurney襟翼對流場影響和揭示繞流機(jī)理[5]。韓中合等模擬NACA4412翼型流場數(shù)值,分析鈍尾緣翼型風(fēng)力機(jī)性能優(yōu)于原翼型風(fēng)力機(jī)[6]。楊科等研究上風(fēng)面開縫、下風(fēng)面開縫及由下向上開縫三種方案,評價(jià)不同開縫方案對風(fēng)力機(jī)靜態(tài)失速特性影響[7]。李巖等對風(fēng)輪上方加聚風(fēng)罩及在葉片后加小翼開展研究,分析風(fēng)機(jī)氣動性能提高原因[8-9]。王建明等對尾緣噴氣垂直軸風(fēng)力機(jī)氣動性能研究,發(fā)現(xiàn)尾緣噴氣可改變?nèi)~片表面壓力系數(shù)[10]。楊從新等對低風(fēng)速下直線翼垂直軸風(fēng)力機(jī)氣動性能作大量研究[11]。趙振宙等基于擾流技術(shù)模擬研究直線翼垂直軸風(fēng)輪葉片性能,提出采用擾流方法解決因局部方位角葉片攻角極小導(dǎo)致的整體風(fēng)輪性能較低問題[12]。張立茹等研究風(fēng)力機(jī)葉尖有無小翼三維流場特性數(shù)值,雖針對水平軸風(fēng)力機(jī),但可對垂直軸風(fēng)力機(jī)提供借鑒作用[13]。

    本研究采用現(xiàn)有直線翼垂直軸風(fēng)力機(jī)典型翼型,即NACA系列對稱翼型,保留其最主要主體氣動外形,僅改造尾緣部分,加裝具有一定形狀彎板,改善葉片尾緣附近流動狀況,提高葉片轉(zhuǎn)矩特性。選取兩種對稱翼型NACA0024和NACA0012,用于對比彎板參數(shù)對厚翼和薄翼影響規(guī)律。采用數(shù)值模擬方法,計(jì)算包括不加裝彎板在內(nèi)具有7種彎板長度2葉片直線翼垂直軸風(fēng)力機(jī)轉(zhuǎn)矩特性,包括動態(tài)和靜態(tài),分析加裝彎板后葉片周圍流場變化,探討尾緣加彎板對風(fēng)輪轉(zhuǎn)矩特性影響機(jī)理。

    1 尾緣加彎板葉片設(shè)計(jì)

    整體看,應(yīng)用NACA四位數(shù)系列對稱翼型可使直線翼垂直軸風(fēng)力機(jī)獲得較好輸出特性[4]。但在高轉(zhuǎn)速下,產(chǎn)生分離流與大量渦旋影響葉片周圍流動,降低氣動特性。為此,本研究選擇兩種NACA系列對稱翼型為基礎(chǔ)翼型,即NACA0024和NACA0012,二者外形輪廓一致,最大厚度相差一倍,對比翼型厚度影響。在兩種翼型尾緣后部加入彎板,改善尾緣部分流動情況。原始翼型與加入彎板翼型見圖1、2。圖中V-來流風(fēng)速,C-翼型弦長,W-彎板末端在來流方向投影長度,F(xiàn)-彎板末端在垂直來流方向投影長度。本研究設(shè)計(jì)直線翼垂直軸風(fēng)力機(jī)尺寸以東北農(nóng)業(yè)大學(xué)低速風(fēng)洞試驗(yàn)段尺寸為依據(jù)(1m×1m)。風(fēng)輪直徑(D)為0.8m,風(fēng)輪高度(H)為0.5m,葉片數(shù)為2個(gè)?;A(chǔ)葉片翼型弦長(C)為0.125m,彎板安裝在尾緣部分。為簡化設(shè)計(jì),彎板參數(shù)僅考慮:彎板長度與彎度,兼顧考慮實(shí)際葉片與彎板強(qiáng)度問題。因此,彎板彎度曲率不宜過大,設(shè)定F固定為10%C。彎板始端與翼型結(jié)合寬度為0.003m并作一定圓角過渡,且在結(jié)合處彎板與來流方向相切。本研究定義彎板長度與弦長比為彎弦比(ε),見式(1)。

    圖1 NACA0024翼型與尾緣加彎板翼型Fig.1 NACA0024 air foiland theair foilwith curved plateat trailing edge

    圖2 NACA0012翼型與尾緣加彎板翼型Fig.1 NACA0012 air foiland theair foilw ith curved p lateat trailing edge

    因彎弦比過小無法發(fā)揮作用,過大影響原始翼型氣動特性,為探討不同彎弦比影響,本文定義包括原始翼型在內(nèi)的7種彎弦比,分別為0、10%、20%、30%、40%、50%和60%。并給出各種彎弦比對應(yīng)負(fù)荷,結(jié)果見表1。

    表1 彎弦比Tab le1 Ratio of curved plate length and airfoil chord

    2 數(shù)值模擬方法

    2.1 計(jì)算模型

    為充分模擬葉片周圍流場特性,葉片周圍流場取長方形,風(fēng)輪中心前方區(qū)域長度為5D,后方區(qū)域長度為10D,風(fēng)輪中心上下寬度各為5D。整個(gè)計(jì)算域分為兩部分,內(nèi)計(jì)算域?yàn)?.5倍風(fēng)輪直徑,外計(jì)算域?yàn)殚L方形除去內(nèi)計(jì)算域剩余部分。兩個(gè)計(jì)算域采用非結(jié)構(gòu)網(wǎng)格,葉片表面采用網(wǎng)格加密。葉片起始位置即為0°,逆時(shí)針旋轉(zhuǎn)。計(jì)算域如圖3所示,葉片表面網(wǎng)格加密如圖4所示。

    入口采用速度入口邊界條件,風(fēng)速為10m·s-1,出口為自由出流壓力出口邊界條件,上下邊界為壁面,葉片及彎板為無滑移邊界條件。計(jì)算采用RNG k-ε湍流模型,擴(kuò)散項(xiàng)采用中心差分格式,對流項(xiàng)采用二階迎風(fēng)格式,壓力速度耦合采用SIMPLE算法。

    圖3 計(jì)算域Fig.3 Com putationaldomain diagram

    圖4 葉片附近網(wǎng)格Fig.4 Air foil Localgrid

    2.2 網(wǎng)格無關(guān)性驗(yàn)證

    為確定合適網(wǎng)格數(shù)目,保證計(jì)算精度,作網(wǎng)格無關(guān)性驗(yàn)證。以NACA0024為例,分別采用153988、211095、290358和370081計(jì)算網(wǎng)格數(shù)量,計(jì)算結(jié)果見圖5。當(dāng)網(wǎng)格數(shù)>211095時(shí),其轉(zhuǎn)矩系數(shù)趨于穩(wěn)定,力矩系數(shù)波動<3%,基本達(dá)到本研究計(jì)算要求。綜合考慮計(jì)算時(shí)間和周期影響,本研究采用網(wǎng)格數(shù)量為212000。

    圖5 網(wǎng)格數(shù)與轉(zhuǎn)矩系數(shù)Fig.5 Torque coefficient under d ifferent numbersof Grid

    3 結(jié)果與分析

    3.1 尾緣加彎板對風(fēng)力機(jī)輸出轉(zhuǎn)矩特性影響

    風(fēng)速為10m·s-1時(shí),具有不同彎弦比風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)(CM)與葉尖速比(λ)間關(guān)系見圖6。CM與λ定義為:

    式中,M-風(fēng)力機(jī)轉(zhuǎn)矩(N·m);U-風(fēng)速(m·s-1);R-風(fēng)輪半徑(m);A-風(fēng)輪掃略面積(m2);ρ-空氣密度(kg·m-3);ω-風(fēng)輪角速度(rad·s-1)。

    從整體上看,NACA0024翼型風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)高于相同條件下采用NACA0012翼型風(fēng)力機(jī)。在安裝彎板后,從風(fēng)力機(jī)轉(zhuǎn)矩特性改善幅度上看,NACA0024翼型風(fēng)力機(jī)效果明顯。由此可知,采用厚翼型NACA系列對稱翼型風(fēng)力機(jī)轉(zhuǎn)矩特性更好,彎板對改善轉(zhuǎn)矩系數(shù)作用效果明顯。從旋轉(zhuǎn)特性看,采用NACA0024翼型風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)最大值出現(xiàn)在尖速比1.6附近,而采用NACA0012翼型風(fēng)力機(jī)則出現(xiàn)在尖速比2附近,說明薄翼型風(fēng)力機(jī)可獲得較高轉(zhuǎn)速,而厚翼型可在較低轉(zhuǎn)速下獲得相對較高轉(zhuǎn)矩系數(shù)。由此可知,如需單純轉(zhuǎn)矩特性較高風(fēng)力機(jī),如風(fēng)力提水機(jī),可采用厚翼型;如需較高輸出功率,如風(fēng)力發(fā)電機(jī),可考慮稍薄翼型。

    采用NACA0012翼型風(fēng)力機(jī),加裝彎板后,在低尖速比下,如λ=0.4時(shí),風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)呈下降趨勢。但隨著尖速比增大,彎板作用顯現(xiàn),直到λ=1.6前,帶有彎板風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)均高于無彎板風(fēng)力機(jī)。但當(dāng)尖速比繼續(xù)增大時(shí),無彎板風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)逐漸高于有彎板風(fēng)力機(jī),在λ=2.4時(shí),無彎板風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)達(dá)到最大且高于有彎板轉(zhuǎn)矩系數(shù)。但有彎板風(fēng)力機(jī)最大轉(zhuǎn)矩系數(shù)出現(xiàn)在λ=2時(shí),彎弦比在40%以下風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)高于無彎板風(fēng)力機(jī),其中,彎弦比為30%風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)最大。

    圖6 不同彎弦比風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)Fig.6 Torque coefficien t ofw ind turbinesw ith differen tε

    對于采用NACA0024翼型風(fēng)力機(jī),除彎弦比為10%風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)整體低于無彎板風(fēng)力機(jī),彎弦比為20%風(fēng)力機(jī)基本與無彎板風(fēng)力機(jī)轉(zhuǎn)矩特性相近外,其他彎弦比風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)均較無彎板風(fēng)力機(jī)有不同程度改善。

    λ>0.4后,有彎板風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)大幅高于無彎板風(fēng)力機(jī),說明加彎板后葉片在低尖速比下性能較好,可有效改善低轉(zhuǎn)速下性能不穩(wěn)定缺點(diǎn)。在λ=1.6時(shí),風(fēng)力機(jī)獲得最大轉(zhuǎn)矩系數(shù),彎弦比為40%風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)最大,較無彎板風(fēng)力機(jī)提高17%。

    3.2 尾緣加彎板對葉片周圍流場影響

    研究一個(gè)旋轉(zhuǎn)周期內(nèi)不同風(fēng)輪旋轉(zhuǎn)角θ下風(fēng)力機(jī)轉(zhuǎn)矩特性。由于所選風(fēng)力機(jī)為兩葉片,風(fēng)力機(jī)從起始位置0°(見圖3)逆時(shí)針旋轉(zhuǎn)180°即為一個(gè)旋轉(zhuǎn)周期。由圖6可知,采用NACA0024翼型風(fēng)力機(jī)轉(zhuǎn)矩特性較NACA0012翼型風(fēng)力機(jī)性能好,所以采用NACA0024翼型風(fēng)力機(jī)獲得最大轉(zhuǎn)矩系數(shù)尖速比λ=1.6為對象,選取三種典型彎弦比:10%c、40% c、60%c與基礎(chǔ)翼型,計(jì)算風(fēng)力機(jī)在該尖速比下,一個(gè)旋轉(zhuǎn)周期內(nèi)轉(zhuǎn)矩系數(shù)隨旋轉(zhuǎn)角變化情況,并對比分析。計(jì)算間隔為10°,計(jì)算結(jié)果見圖7。

    圖7 不同旋轉(zhuǎn)角處風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)(λ=1.6)Fig.7 Torque Coefficient at different azim uth angles(λ=1.6)

    由圖7可知,與整體轉(zhuǎn)矩特性一致,采用NACA0024翼型風(fēng)力機(jī)在一個(gè)旋轉(zhuǎn)周期內(nèi)轉(zhuǎn)矩系數(shù)整體高于NACA0012翼型風(fēng)力機(jī)。兩種風(fēng)力機(jī)獲得最大轉(zhuǎn)矩系數(shù)旋轉(zhuǎn)角均在90°附近,而加入彎板后,獲得最大轉(zhuǎn)矩系數(shù)角度略滯后,最大轉(zhuǎn)矩系數(shù)有不同程度提高。在獲得最大轉(zhuǎn)矩系數(shù)之前旋轉(zhuǎn)角處,無彎板風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)均高于有彎板風(fēng)力機(jī),而在此旋轉(zhuǎn)角后,有彎板風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)逐漸高于無彎板風(fēng)力機(jī),說明彎板發(fā)揮改善轉(zhuǎn)矩特性作用實(shí)際旋轉(zhuǎn)角為90°~180°。采用NACA0024翼型風(fēng)力機(jī),彎弦比為10%風(fēng)力機(jī)在一個(gè)旋轉(zhuǎn)周期內(nèi)轉(zhuǎn)矩特性較無彎板風(fēng)力機(jī)低,說明彎板過短不利改善各個(gè)旋轉(zhuǎn)角處轉(zhuǎn)矩特性,使轉(zhuǎn)矩系數(shù)降低。而對于彎弦比為60%情況,雖然在旋轉(zhuǎn)角90°以后轉(zhuǎn)矩特系數(shù)明顯提高,但在0°~90°,轉(zhuǎn)矩系數(shù)大幅減少,整體不理想。而彎弦比為40%風(fēng)力機(jī),在0°~90°轉(zhuǎn)矩系數(shù)略下降,后轉(zhuǎn)矩系數(shù)大幅提升,在各種彎弦比中轉(zhuǎn)矩改善效果最理想。采用NACA0012翼型風(fēng)力機(jī)變化趨勢與NACA0024翼型風(fēng)力機(jī)基本一致。

    以NACA0024翼型風(fēng)力機(jī)為例,選取兩個(gè)典型旋轉(zhuǎn)角,即30°和100°,計(jì)算在λ=1.6時(shí)具有4種彎弦比風(fēng)力機(jī)葉片周圍速度和壓力分布,分析彎板對翼型周圍流場的影響(見圖8、9)。

    由圖8可知,對無彎板風(fēng)力機(jī),上葉片腹面存在一片負(fù)壓區(qū),在該葉片上產(chǎn)生氣動力合力較大,為提供風(fēng)力機(jī)逆時(shí)針旋轉(zhuǎn)主要作用力。而下葉片腹面與背面壓力差較小,無法提供使葉片旋轉(zhuǎn)氣動力??梢娫谖簿壸髠?cè)存在一個(gè)負(fù)壓區(qū),使風(fēng)作用在該葉片上產(chǎn)生氣動力合力指向這一區(qū)域,該氣動力將驅(qū)使風(fēng)力機(jī)沿順時(shí)針旋轉(zhuǎn),抵消一部分上葉片逆時(shí)針轉(zhuǎn)矩。但由于上葉片提供轉(zhuǎn)矩較大,在該角度下風(fēng)輪轉(zhuǎn)矩仍為正值。而對有彎板風(fēng)力機(jī),上下葉片周圍流動情況發(fā)生變化較大。當(dāng)彎弦比為10%時(shí),上葉片腹面負(fù)壓區(qū)減小,該葉片提供逆時(shí)針轉(zhuǎn)矩減少;在彎板兩側(cè)出現(xiàn)回流現(xiàn)象,形成渦旋,產(chǎn)生能量損失,加劇上葉片氣動力降低。而下葉片整體流動情況無較大改變,僅尾緣左側(cè)負(fù)壓區(qū)減少。因此,綜合作用結(jié)果使該角度下風(fēng)力機(jī)轉(zhuǎn)矩幾乎為零,甚至出現(xiàn)負(fù)值,即有反轉(zhuǎn)趨勢。彎弦比為40%時(shí),雖然上葉片腹面負(fù)壓區(qū)較彎弦比為10%時(shí)增大,但仍低于無彎板風(fēng)力機(jī)。而下葉片尾緣由于彎板長度增加,負(fù)壓區(qū)明顯增大,該葉片提供負(fù)轉(zhuǎn)矩也增大。因此,總體來看,在該角度下風(fēng)力機(jī)轉(zhuǎn)矩仍很小,且為負(fù)值。當(dāng)彎弦比增大到60%時(shí),上葉片負(fù)壓區(qū)基本與無彎板相同,但在下葉片尾緣左側(cè)負(fù)壓區(qū)顯著增大,且負(fù)壓值較大,導(dǎo)致該葉片提供順時(shí)針轉(zhuǎn)矩比上葉片提供逆時(shí)針轉(zhuǎn)矩大,在該角度下風(fēng)力機(jī)轉(zhuǎn)矩出現(xiàn)較大負(fù)值。由此可見,葉片尾緣加裝彎板后,主要影響上葉片腹面負(fù)壓區(qū)和下葉片尾緣負(fù)壓區(qū)域變化,改變?nèi)~片氣動特性和風(fēng)輪轉(zhuǎn)矩特性。

    圖8 不同彎弦比風(fēng)力機(jī)葉片周圍流場(NACA0024,λ=1.6,θ=30°)Fig.8 Flow fieldsaround bladeairfoilwith differentε(NACA0024,λ=1.6,θ=30°)

    由圖9可知,在該角度下下葉片提供風(fēng)力機(jī)旋轉(zhuǎn)轉(zhuǎn)矩,在該葉片背面有一個(gè)正壓區(qū),腹面有一個(gè)較大負(fù)壓區(qū),形成較大氣動合力,且該壓力區(qū)集中在前緣,為該葉片提供逆時(shí)針轉(zhuǎn)矩。而對于上葉片,壓力變化集中在前緣附近,作用在該葉片氣動力為風(fēng)力機(jī)提供逆時(shí)針轉(zhuǎn)矩。彎弦比為10%,上葉片周圍流動變化較小,但下葉片腹面負(fù)壓區(qū)變小,因此風(fēng)輪轉(zhuǎn)矩較無彎板情況略低。彎弦比為40%,同樣上葉片周圍流動變化不大,但下葉片腹面壓力區(qū)域顯著增大,且存在一定環(huán)流,因此該葉片提供轉(zhuǎn)矩增大,風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)增大。而當(dāng)彎弦比增大到60%時(shí),彎板長度增加,該趨勢愈加明顯,雖然上葉片迎風(fēng)面出現(xiàn)較大渦旋影響該葉片氣動力,但在總體上風(fēng)輪轉(zhuǎn)矩仍繼續(xù)增大。由此可見,對于圖7所示在旋轉(zhuǎn)角90°以后有彎板風(fēng)力機(jī)轉(zhuǎn)矩系數(shù)提高主要原因,是彎板存在增大下葉片腹面負(fù)壓區(qū)。

    3.3 尾緣加彎板對風(fēng)力機(jī)靜態(tài)起動轉(zhuǎn)矩特性影響

    直線翼垂直軸風(fēng)力機(jī)在靜止?fàn)顟B(tài)下起動轉(zhuǎn)矩較低是制約風(fēng)力機(jī)表態(tài)起動特性主要因素。在此,以采用NACA0024翼型風(fēng)力機(jī)為例,選取兩種彎弦比,即動態(tài)輸出轉(zhuǎn)矩特性較好彎弦比40%與基礎(chǔ)翼型,計(jì)算風(fēng)力機(jī)在靜止?fàn)顟B(tài)(λ=0)時(shí)一個(gè)轉(zhuǎn)動周期內(nèi)靜態(tài)轉(zhuǎn)矩系數(shù)隨旋轉(zhuǎn)角變化情況,計(jì)算角度間隔為10°。

    靜態(tài)轉(zhuǎn)矩計(jì)算方法按照公式(2),計(jì)算結(jié)果如圖10所示。

    圖10 不同彎弦比風(fēng)力機(jī)靜態(tài)起動轉(zhuǎn)矩Fig.10 Static torque coefficien t of wind turbinesw ith differentε

    圖11 不同彎弦比風(fēng)力機(jī)葉片周圍流場(NACA0024,λ=0,θ=0°)Fig.11 Flow fieldsaround bladeair foilwith differentε(NACA0024,λ=0,θ=0°)

    由圖10可見,兩條曲線變化趨勢基本一致,當(dāng)旋轉(zhuǎn)角為0~120°時(shí),加彎板風(fēng)力機(jī)靜態(tài)轉(zhuǎn)矩系數(shù)與無彎板風(fēng)力機(jī)基本一致,僅在某些角度下略有提高。但當(dāng)旋轉(zhuǎn)角為120~180°時(shí),靜態(tài)轉(zhuǎn)矩系數(shù)提高顯著,尤其旋轉(zhuǎn)角為150~180°時(shí),無彎板風(fēng)力機(jī)靜態(tài)轉(zhuǎn)矩較低,甚至有負(fù)值(反向轉(zhuǎn)矩),但加裝彎板后,不但轉(zhuǎn)矩系數(shù)大幅提高,且消除反向轉(zhuǎn)矩。進(jìn)一步計(jì)算可得到兩種風(fēng)機(jī)一個(gè)旋轉(zhuǎn)周期內(nèi)平均靜態(tài)轉(zhuǎn)矩系數(shù),采用原始翼型風(fēng)力機(jī)平均靜態(tài)轉(zhuǎn)矩系數(shù)為0.0250,而采用具有彎弦比40%變形葉片風(fēng)力機(jī)平均靜態(tài)轉(zhuǎn)矩系數(shù)為0.0316,提升26.4%,說明葉片尾緣加彎板對風(fēng)力機(jī)整體起動性能提高效果良好。

    選取兩個(gè)典型旋轉(zhuǎn)角,即0°和160°,計(jì)算這兩種風(fēng)力機(jī)在靜止?fàn)顟B(tài)下葉片周圍速度和壓力分布,分析彎板對翼型周圍流場影響,如圖11~12所示。

    圖12 不同彎弦比風(fēng)力機(jī)葉片周圍流場(NACA0024,λ=0,θ=160°)Fig.12 Flow fieldsaround b ladeairfoilw ith differentε(NACA0024,λ=0,θ=160°)

    由圖11可知,在旋轉(zhuǎn)角為0°時(shí),上、下葉片周圍流場為對稱分布,不產(chǎn)生氣動力,在此角度下靜態(tài)轉(zhuǎn)矩為0。當(dāng)葉片尾緣安裝彎板后,改變翼型對稱性,在上葉片背面出現(xiàn)負(fù)壓區(qū),在下葉片尾緣彎板處與前緣前部均出現(xiàn)環(huán)流,改變?nèi)~片上下表面壓力分布,產(chǎn)生氣動力,形成轉(zhuǎn)矩,有效改善該角度下風(fēng)力機(jī)起動性能。

    由圖12可知,在該角度下,無彎板風(fēng)力機(jī)上葉片負(fù)壓區(qū)僅出現(xiàn)在葉片背面前緣附近,背面存在較大渦旋區(qū),流動分離現(xiàn)象嚴(yán)重。而下葉片腹面同樣存在兩個(gè)較大渦旋和流動分離,即該角度下葉片氣動特性較差,風(fēng)輪轉(zhuǎn)矩較低。對于有彎板風(fēng)力機(jī),由于在尾緣加入彎板,接受來流面積增加,使上葉片背面負(fù)壓區(qū)域增大,腹面正壓區(qū)增大,渦旋移動到葉片后方,減少流動分離;同時(shí)改善下葉片流動分離,葉片氣動特性充分發(fā)揮,風(fēng)輪轉(zhuǎn)矩大幅提升。

    4 結(jié)論

    本研究提出一種在葉片尾緣加裝彎板提高直線翼垂直軸風(fēng)力機(jī)轉(zhuǎn)矩特性方法,利用數(shù)值模擬計(jì)算具有不同彎板長度風(fēng)力機(jī)動態(tài)輸出轉(zhuǎn)矩和靜態(tài)起動特性,結(jié)論如下:

    a.在NACA系列對稱翼型尾緣加彎板葉片,可在一定程度上有效改善直線翼垂直軸風(fēng)力機(jī)輸出轉(zhuǎn)矩特性和靜態(tài)起動特性,可作為提升直線翼垂直軸風(fēng)力機(jī)氣動特性方法,為相關(guān)研究提供借鑒。

    b.彎板加在厚翼型葉片改善風(fēng)力機(jī)輸出轉(zhuǎn)矩特性效果顯著。彎板長度與翼型弦長比例要適當(dāng),彎板長度為弦長40%風(fēng)力機(jī)轉(zhuǎn)矩特性改善效果最佳,最大輸出轉(zhuǎn)矩系數(shù)較無彎板風(fēng)力機(jī)提高17%,一個(gè)旋轉(zhuǎn)周期內(nèi)靜態(tài)平均起動轉(zhuǎn)矩提高26.4%,可消除反向轉(zhuǎn)矩,風(fēng)機(jī)整體起動性能提高。

    c.葉片尾緣加裝彎板后可在某些旋轉(zhuǎn)角下改變?nèi)~片上下表面壓力分布,減少渦旋和流動分離,提高該旋轉(zhuǎn)角下氣動特性,改善風(fēng)力機(jī)轉(zhuǎn)矩特性。

    [1]李巖.垂直軸風(fēng)力機(jī)應(yīng)用及其發(fā)展前景[J].可再生能源,2009 (6):118-120.

    [2]高偉,李春,葉舟,等.垂直軸風(fēng)力機(jī)動態(tài)特性及氣動性能[J].排灌機(jī)械工程學(xué)報(bào),2014,32(2):146-150.

    [3]舒超,蔡新,潘盼,等.垂直軸風(fēng)力機(jī)改進(jìn)翼型氣動及功率特性研究[J].水電能源科學(xué),2015(1):162-165.

    [4]廖書學(xué),李春,聶佳斌,等.不同翼型對垂直軸風(fēng)力機(jī)性能的影響[J].機(jī)械設(shè)計(jì)與研究,2011,27(3):108-111.

    [5]趙萬里,劉沛清,朱建勇,等.Gurney襟翼對風(fēng)力機(jī)流動控制數(shù)值研究[J].電網(wǎng)與清潔能源,2011,27(9):85-92.

    [6]韓中合,焦紅瑞等.加裝鈍尾緣改善風(fēng)力機(jī)槳葉氣動性能的研究[J].動力工程學(xué)報(bào),2009,29(11):1073-1077.

    [7]楊科,王會社,徐建中,等.開縫式風(fēng)力機(jī)靜態(tài)失速特性的研究[J].工程熱物理學(xué)報(bào),2008,29(1):32-35.

    [8]李巖,唐靜,田川公太朗,等.圓臺型聚風(fēng)罩對垂直軸風(fēng)力機(jī)起動性影響[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2016,47(4):95-101.

    [9]李巖,鄭玉芳,唐靜,等.葉片后加小翼垂直軸風(fēng)力機(jī)氣動特性數(shù)值模擬[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2016,47(7):76-81.

    [10]王建明,譚永志,溫學(xué)兵,等.葉片尾緣噴氣垂直軸風(fēng)力機(jī)氣動性能研究[J].渤海大學(xué)學(xué)報(bào):自然科學(xué)版,2013,34(4):386-392.

    [11]楊從新,李壽圖,王秀勇.低風(fēng)速下H型垂直軸風(fēng)力機(jī)氣動性能[J].排灌機(jī)械工程學(xué)報(bào),2014,32(10):888-893.

    [12]趙振宙,陳潘浩,王同光,等.基于擾流技術(shù)直葉片升力型垂直軸風(fēng)輪性能改善[J].機(jī)械工程學(xué)報(bào),2016,52:146-152.

    [13]張立茹,汪建文,曲立群,等.風(fēng)力機(jī)葉尖有無小翼三維流場的數(shù)值研究[J].太陽能學(xué)報(bào),2011(11):1599-1604.

    Numerical simulation on torque characteristics of vertical axis windturbine with curved plate adding at blade trailing edge

    LI Yan1,2,ZHANGTingt ing1,Tagawa Kotaro3,FENG Fang4,ZHENG Yufang1,GONG Weij ia1,2
    (1. School of Engineering,Northeast Agricultural University, Harbin 150030, China; 2. Heilongjiang Provincial Key Laboratory ofTechnology and Equipment for Utilization of Agricultural Renewable Resources in Cold Region, Harbin150030, China; 3. Faculty of Regional Sciences, Tottori University, Tottori 6808552, Japan; 4. School ofSciences, Northeast Agricultural University, Harbin 150030, China)

    In this paper, a kind of blade with curved plate adding at its trailing edge was proposed forstraight-bladed vertical axis wind turbine (SB-VAWT) in order to improve the flow field around blade and tothe torque characteristics. Two typical blade airfoils, NACA0024 and NACA0012, often used for SB-VAWTwere selected for numerical simulation. The torque characteristics and static torque performance werecalculated for a SB-VAWT with 7 kinds of designed curved plate length including without curved plate. Thesimulation results showed that the curved plate adding at the trailing edge of the NACA series symmetricalairfoil could improve the torque characteristics and static starting performance to a certain degree. Thepressure distribution can be changed and the vortex and separation could be restrained due to the presenceof a curved plate, which improve the aerodynamics performance of blade and the torque characteristics ofwind turbine. Among all the curved plates researched, the maximum torque coefficient comes for the turbinewith the curved plate which its ratio of curved plate length to chord equal 40%, is 17% larger than the turbinewithout curved plate. The statics torque coefficient was also 26.4% lager.

    vertical axis wind turbine (VAWT); trailing edge of blade; curved plate; torqueperformance; numerical simulation

    TM315

    A

    1005-9369(2017)03-0072-08

    時(shí)間2017-3-21 14:04:00[URL]http://kns.cnki.net/kcms/detail/23.1391.S.20170321.1404.018.htm l

    李巖,張婷婷,田川公太朗,等.葉片尾緣加彎板垂直軸風(fēng)力機(jī)轉(zhuǎn)矩特性數(shù)值模擬[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2017,48(3):72-79.

    Li Yan, Zhang Tingting, Tagawa Kotaro, et al. Numerical simulation on torque characteristics of vertical axis wind turbinewith curved plate adding at blade trailing edge[J]. Journal of Northeast Agricultural University, 2017, 48(3): 72-79. (in Chinesewith English abstract)

    2017-02-12

    黑龍江省教育廳科學(xué)技術(shù)研究項(xiàng)目(12541012);國家自然科學(xué)基金項(xiàng)目(10702015)

    李巖(1972-),男,教授,研究方向?yàn)轱L(fēng)能及可再生能源綜合利用。E-mail:liyanneau@163.com

    猜你喜歡
    垂直軸尾緣風(fēng)輪
    基于強(qiáng)化換熱的偏斜尾緣設(shè)計(jì)
    能源工程(2021年1期)2021-04-13 02:05:46
    垂直軸風(fēng)力機(jī)主軸直徑對氣動性能的影響
    葉片數(shù)目對風(fēng)輪位移和應(yīng)力的影響
    太陽能(2019年10期)2019-10-29 07:25:08
    從五臟相關(guān)理論淺析祛風(fēng)退翳法在風(fēng)輪疾病的應(yīng)用
    翼型湍流尾緣噪聲半經(jīng)驗(yàn)預(yù)測公式改進(jìn)
    具有尾緣襟翼的風(fēng)力機(jī)動力學(xué)建模與恒功率控制
    伸縮葉片式垂直軸風(fēng)機(jī)葉輪的數(shù)值仿真
    并網(wǎng)垂直軸風(fēng)力發(fā)電機(jī)組的大型化技術(shù)研究
    鈍尾緣葉片三維建模方法的研究
    新型雙風(fēng)輪風(fēng)力機(jī)氣動特性的三維流場數(shù)值模擬
    一级av片app| 久久精品国产鲁丝片午夜精品 | 亚洲av美国av| 99久久无色码亚洲精品果冻| 91精品国产九色| 极品教师在线视频| 亚洲中文日韩欧美视频| 黄色配什么色好看| 悠悠久久av| 精品乱码久久久久久99久播| 亚洲四区av| 精品人妻一区二区三区麻豆 | 一区二区三区激情视频| 久久久久久伊人网av| 久久6这里有精品| 欧美中文日本在线观看视频| 国产伦在线观看视频一区| 欧美+亚洲+日韩+国产| 国产欧美日韩精品亚洲av| 高清毛片免费观看视频网站| 老女人水多毛片| 人妻制服诱惑在线中文字幕| 少妇丰满av| 校园人妻丝袜中文字幕| 国产单亲对白刺激| 国产精品一区二区三区四区免费观看 | 真人做人爱边吃奶动态| 两个人的视频大全免费| 国产精品三级大全| 一边摸一边抽搐一进一小说| 国产精品久久电影中文字幕| 亚洲经典国产精华液单| 日韩欧美一区二区三区在线观看| 18禁裸乳无遮挡免费网站照片| 如何舔出高潮| 成人综合一区亚洲| 国产单亲对白刺激| 日本黄大片高清| 亚洲,欧美,日韩| netflix在线观看网站| 亚洲精华国产精华液的使用体验 | 国产亚洲91精品色在线| 欧美黑人欧美精品刺激| 乱人视频在线观看| 精品免费久久久久久久清纯| 色精品久久人妻99蜜桃| 国产欧美日韩一区二区精品| 精品欧美国产一区二区三| 午夜激情福利司机影院| 色视频www国产| 亚洲久久久久久中文字幕| 老师上课跳d突然被开到最大视频| 亚洲av成人精品一区久久| 欧美成人免费av一区二区三区| 亚洲av中文字字幕乱码综合| 婷婷精品国产亚洲av| 精品人妻1区二区| a级毛片免费高清观看在线播放| 国产麻豆成人av免费视频| 国产激情偷乱视频一区二区| 自拍偷自拍亚洲精品老妇| 中文字幕熟女人妻在线| 色av中文字幕| 亚洲中文字幕日韩| 有码 亚洲区| 国产真实乱freesex| 免费观看精品视频网站| 日本五十路高清| or卡值多少钱| 老司机福利观看| 日本熟妇午夜| 国产伦精品一区二区三区视频9| 啪啪无遮挡十八禁网站| 中文字幕人妻熟人妻熟丝袜美| 波野结衣二区三区在线| АⅤ资源中文在线天堂| 99热这里只有是精品在线观看| 91在线观看av| 精品久久久久久久久亚洲 | 日韩精品中文字幕看吧| 一进一出好大好爽视频| 亚洲国产色片| 午夜日韩欧美国产| 欧美日韩精品成人综合77777| 嫩草影院新地址| 亚洲欧美精品综合久久99| 色综合色国产| 男女下面进入的视频免费午夜| 久久久久久久精品吃奶| 男女那种视频在线观看| 国内少妇人妻偷人精品xxx网站| 国产一区二区激情短视频| 日韩 亚洲 欧美在线| 波多野结衣巨乳人妻| 97碰自拍视频| 亚洲欧美日韩卡通动漫| 亚洲av成人精品一区久久| 老熟妇仑乱视频hdxx| 成年女人永久免费观看视频| 国产单亲对白刺激| 成年女人看的毛片在线观看| 日韩欧美三级三区| 免费在线观看成人毛片| 国产毛片a区久久久久| 搡女人真爽免费视频火全软件 | 国产亚洲精品久久久com| 亚洲在线自拍视频| 免费看a级黄色片| 最近中文字幕高清免费大全6 | 午夜福利18| 免费在线观看日本一区| 亚洲乱码一区二区免费版| 国产黄色小视频在线观看| 成年免费大片在线观看| 国产av麻豆久久久久久久| 在线观看免费视频日本深夜| 亚洲精品成人久久久久久| 亚洲av.av天堂| 免费观看人在逋| 中文字幕精品亚洲无线码一区| 色尼玛亚洲综合影院| 欧美绝顶高潮抽搐喷水| 国产91精品成人一区二区三区| 久久九九热精品免费| 啦啦啦韩国在线观看视频| 国产乱人视频| а√天堂www在线а√下载| 亚洲久久久久久中文字幕| 国产精品久久久久久久电影| 欧美一区二区国产精品久久精品| 九九爱精品视频在线观看| 国产熟女欧美一区二区| 欧美性猛交╳xxx乱大交人| 一级av片app| 国语自产精品视频在线第100页| 一进一出好大好爽视频| 美女 人体艺术 gogo| 91狼人影院| 欧美bdsm另类| 久久久国产成人免费| 中文资源天堂在线| 美女高潮喷水抽搐中文字幕| 两人在一起打扑克的视频| 老熟妇乱子伦视频在线观看| 亚洲熟妇熟女久久| 老司机福利观看| 国产精品久久久久久久电影| 美女xxoo啪啪120秒动态图| 人妻久久中文字幕网| 欧美日韩综合久久久久久 | 蜜桃久久精品国产亚洲av| 99视频精品全部免费 在线| 波多野结衣高清作品| 一级黄片播放器| 欧美不卡视频在线免费观看| 亚洲精品粉嫩美女一区| 老司机午夜福利在线观看视频| 变态另类成人亚洲欧美熟女| 国产不卡一卡二| 亚洲av二区三区四区| 亚洲图色成人| 国产精品久久电影中文字幕| 真实男女啪啪啪动态图| 18+在线观看网站| 老熟妇乱子伦视频在线观看| 日本-黄色视频高清免费观看| h日本视频在线播放| 成人国产一区最新在线观看| 少妇丰满av| 99久久精品热视频| 亚洲最大成人av| 美女cb高潮喷水在线观看| 精品午夜福利视频在线观看一区| 在线国产一区二区在线| 色综合亚洲欧美另类图片| 国产真实乱freesex| 亚洲va日本ⅴa欧美va伊人久久| 久久精品国产亚洲av香蕉五月| 国产淫片久久久久久久久| 尾随美女入室| 一区福利在线观看| 国产麻豆成人av免费视频| 少妇熟女aⅴ在线视频| 色噜噜av男人的天堂激情| 最后的刺客免费高清国语| 成人一区二区视频在线观看| 亚洲美女视频黄频| 最新在线观看一区二区三区| 日本成人三级电影网站| 国产精品一区二区三区四区久久| 极品教师在线视频| 悠悠久久av| 亚洲av第一区精品v没综合| 精品人妻一区二区三区麻豆 | 女同久久另类99精品国产91| 日本爱情动作片www.在线观看 | 动漫黄色视频在线观看| 男女之事视频高清在线观看| 真人做人爱边吃奶动态| 久久中文看片网| 亚洲人与动物交配视频| 全区人妻精品视频| 欧美日韩乱码在线| 联通29元200g的流量卡| 99热网站在线观看| 久久精品国产99精品国产亚洲性色| 国产精品亚洲一级av第二区| 天美传媒精品一区二区| 久久久久久久久中文| 男女做爰动态图高潮gif福利片| 亚洲国产色片| 国产亚洲精品久久久久久毛片| 91久久精品国产一区二区三区| 亚洲图色成人| 午夜福利18| 日本欧美国产在线视频| 亚洲色图av天堂| 床上黄色一级片| 精品日产1卡2卡| 国内精品美女久久久久久| 69人妻影院| 99热精品在线国产| 美女黄网站色视频| 91麻豆av在线| 此物有八面人人有两片| 一本一本综合久久| 婷婷精品国产亚洲av| 狂野欧美激情性xxxx在线观看| 尤物成人国产欧美一区二区三区| 自拍偷自拍亚洲精品老妇| 午夜福利视频1000在线观看| 婷婷精品国产亚洲av| 亚洲狠狠婷婷综合久久图片| 中文字幕久久专区| 99久久中文字幕三级久久日本| 亚洲一区高清亚洲精品| 91久久精品电影网| a在线观看视频网站| 在线播放无遮挡| 亚洲综合色惰| 女人被狂操c到高潮| 两个人的视频大全免费| 午夜福利视频1000在线观看| .国产精品久久| 俄罗斯特黄特色一大片| 午夜视频国产福利| 国产av不卡久久| 欧美日韩亚洲国产一区二区在线观看| 婷婷精品国产亚洲av在线| 久久久久国内视频| 国语自产精品视频在线第100页| 露出奶头的视频| 最近最新中文字幕大全电影3| 国产精品无大码| 在现免费观看毛片| .国产精品久久| 最近视频中文字幕2019在线8| 亚洲第一电影网av| av国产免费在线观看| av女优亚洲男人天堂| 亚洲色图av天堂| av在线天堂中文字幕| 国产又黄又爽又无遮挡在线| 亚洲一区二区三区色噜噜| 成人国产麻豆网| 国产精品99久久久久久久久| 99国产极品粉嫩在线观看| 亚洲av美国av| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩高清在线视频| 熟女电影av网| 我的女老师完整版在线观看| 日本欧美国产在线视频| 精品人妻一区二区三区麻豆 | 国产 一区 欧美 日韩| 免费在线观看影片大全网站| 校园人妻丝袜中文字幕| 欧美日韩亚洲国产一区二区在线观看| 人妻夜夜爽99麻豆av| 69人妻影院| 国产单亲对白刺激| 色综合色国产| 免费高清视频大片| 国产麻豆成人av免费视频| 成年版毛片免费区| 亚洲国产精品sss在线观看| 99国产精品一区二区蜜桃av| 久久人妻av系列| 老熟妇仑乱视频hdxx| 亚洲最大成人av| 精品日产1卡2卡| 在线观看av片永久免费下载| 我要看日韩黄色一级片| 亚洲精品乱码久久久v下载方式| 国产一区二区三区视频了| 少妇的逼好多水| 成人二区视频| 在线观看66精品国产| 亚洲精品色激情综合| 国产欧美日韩一区二区精品| 成人精品一区二区免费| 国产乱人视频| 舔av片在线| 别揉我奶头 嗯啊视频| 精品久久久久久久久久久久久| 日本黄色视频三级网站网址| 亚洲无线在线观看| 久久久国产成人精品二区| 成人高潮视频无遮挡免费网站| 精品午夜福利在线看| 国产v大片淫在线免费观看| 女同久久另类99精品国产91| 五月玫瑰六月丁香| 久久精品国产亚洲网站| 国产成人一区二区在线| 久久精品国产亚洲av香蕉五月| a在线观看视频网站| 成人一区二区视频在线观看| 午夜免费成人在线视频| 久久精品影院6| 国内精品久久久久久久电影| 欧美激情国产日韩精品一区| 搡女人真爽免费视频火全软件 | 久久精品国产清高在天天线| 亚洲最大成人手机在线| 国产伦精品一区二区三区四那| 国内少妇人妻偷人精品xxx网站| 97人妻精品一区二区三区麻豆| 欧美日韩精品成人综合77777| 亚洲精品国产成人久久av| 亚洲精品久久国产高清桃花| 欧美高清性xxxxhd video| 亚洲人成网站在线播放欧美日韩| 国内揄拍国产精品人妻在线| 男女之事视频高清在线观看| 免费黄网站久久成人精品| 日本成人三级电影网站| 99在线人妻在线中文字幕| 桃色一区二区三区在线观看| 欧美精品啪啪一区二区三区| 亚洲av美国av| 欧美xxxx性猛交bbbb| 亚洲最大成人中文| 久久久久久伊人网av| 久久午夜亚洲精品久久| 午夜免费男女啪啪视频观看 | 亚洲国产高清在线一区二区三| 国产欧美日韩精品亚洲av| 亚洲人成伊人成综合网2020| 国内精品一区二区在线观看| 国产精品久久久久久精品电影| 男人和女人高潮做爰伦理| 人人妻人人看人人澡| 国产不卡一卡二| 日韩亚洲欧美综合| 欧美绝顶高潮抽搐喷水| 美女黄网站色视频| 久久精品国产亚洲网站| 国产精品98久久久久久宅男小说| 美女黄网站色视频| 很黄的视频免费| 欧美不卡视频在线免费观看| 88av欧美| bbb黄色大片| 国产乱人伦免费视频| 精品午夜福利在线看| x7x7x7水蜜桃| 国产精品乱码一区二三区的特点| 国产不卡一卡二| 国产精华一区二区三区| 国产精品不卡视频一区二区| 日韩中字成人| 好男人在线观看高清免费视频| 欧美人与善性xxx| 高清日韩中文字幕在线| 真人一进一出gif抽搐免费| 91在线观看av| 日韩欧美一区二区三区在线观看| 欧美色视频一区免费| 伦精品一区二区三区| 国产精品一及| 精品无人区乱码1区二区| 伊人久久精品亚洲午夜| 精品乱码久久久久久99久播| aaaaa片日本免费| 国语自产精品视频在线第100页| a级毛片免费高清观看在线播放| 一级av片app| 久久香蕉精品热| 亚洲狠狠婷婷综合久久图片| 国产单亲对白刺激| 国产午夜福利久久久久久| 啪啪无遮挡十八禁网站| 真人做人爱边吃奶动态| 在线免费观看的www视频| 成年免费大片在线观看| 久久久久久九九精品二区国产| 99久久无色码亚洲精品果冻| 91狼人影院| 久久久久久久午夜电影| 日本一本二区三区精品| 亚洲美女视频黄频| 99热这里只有是精品在线观看| 麻豆久久精品国产亚洲av| 欧美黑人巨大hd| 欧美日韩国产亚洲二区| 日本精品一区二区三区蜜桃| 久久精品国产亚洲av天美| 神马国产精品三级电影在线观看| 免费观看的影片在线观看| 久久久久国产精品人妻aⅴ院| 91久久精品国产一区二区三区| 在线观看舔阴道视频| 我要搜黄色片| 网址你懂的国产日韩在线| 欧美成人一区二区免费高清观看| 不卡一级毛片| 午夜激情福利司机影院| 国产午夜精品久久久久久一区二区三区 | 美女xxoo啪啪120秒动态图| 一个人看视频在线观看www免费| 日本a在线网址| 国产白丝娇喘喷水9色精品| 中国美女看黄片| 永久网站在线| 99热6这里只有精品| 国产精品无大码| 精品久久久久久成人av| 国产成年人精品一区二区| 国产精品野战在线观看| 免费大片18禁| 成年版毛片免费区| 久久午夜亚洲精品久久| av天堂在线播放| 女的被弄到高潮叫床怎么办 | 欧美最黄视频在线播放免费| 极品教师在线免费播放| 网址你懂的国产日韩在线| 亚洲欧美清纯卡通| 亚洲 国产 在线| 成人二区视频| 国产精品一及| 少妇的逼水好多| 观看美女的网站| 精品福利观看| 国产淫片久久久久久久久| av天堂在线播放| 色av中文字幕| 亚洲18禁久久av| 日韩欧美在线乱码| 少妇猛男粗大的猛烈进出视频 | 国产精品99久久久久久久久| 一个人观看的视频www高清免费观看| h日本视频在线播放| 亚洲最大成人手机在线| av在线观看视频网站免费| 成年女人毛片免费观看观看9| 3wmmmm亚洲av在线观看| 国产又黄又爽又无遮挡在线| 国国产精品蜜臀av免费| 国产精品久久久久久亚洲av鲁大| 亚洲精品影视一区二区三区av| 午夜福利在线观看免费完整高清在 | 久久中文看片网| 免费人成在线观看视频色| av在线天堂中文字幕| 伦理电影大哥的女人| 给我免费播放毛片高清在线观看| 午夜免费成人在线视频| 亚洲美女搞黄在线观看 | 最后的刺客免费高清国语| 悠悠久久av| 亚洲一区高清亚洲精品| 亚洲综合色惰| 小蜜桃在线观看免费完整版高清| 免费看美女性在线毛片视频| 国产精品1区2区在线观看.| 日本欧美国产在线视频| 国产三级在线视频| 午夜久久久久精精品| 亚洲国产欧美人成| 婷婷亚洲欧美| 麻豆国产av国片精品| 99国产极品粉嫩在线观看| 精品99又大又爽又粗少妇毛片 | 特级一级黄色大片| 俺也久久电影网| 国产三级中文精品| 欧美高清性xxxxhd video| 国产精品一区二区三区四区久久| 中文资源天堂在线| 欧美成人一区二区免费高清观看| 91午夜精品亚洲一区二区三区 | 日本精品一区二区三区蜜桃| 搡老熟女国产l中国老女人| 精品久久久久久久末码| 久久热精品热| 婷婷丁香在线五月| 亚洲国产欧洲综合997久久,| 熟妇人妻久久中文字幕3abv| 春色校园在线视频观看| 日本五十路高清| 日本在线视频免费播放| 色综合站精品国产| 成年版毛片免费区| 乱系列少妇在线播放| 欧美日韩精品成人综合77777| 十八禁国产超污无遮挡网站| 日韩强制内射视频| 亚洲经典国产精华液单| 狂野欧美白嫩少妇大欣赏| 黄片wwwwww| 99热这里只有是精品50| 听说在线观看完整版免费高清| 欧美最新免费一区二区三区| 少妇人妻精品综合一区二区 | 欧美区成人在线视频| bbb黄色大片| 成人国产麻豆网| av专区在线播放| а√天堂www在线а√下载| 999久久久精品免费观看国产| 久久久久性生活片| 日韩欧美一区二区三区在线观看| 亚洲中文日韩欧美视频| 国模一区二区三区四区视频| 99国产极品粉嫩在线观看| 高清在线国产一区| 国产精品自产拍在线观看55亚洲| .国产精品久久| 三级毛片av免费| 国产精品av视频在线免费观看| 天堂av国产一区二区熟女人妻| 丝袜美腿在线中文| 身体一侧抽搐| 老司机午夜福利在线观看视频| 国产高清视频在线观看网站| 不卡视频在线观看欧美| 麻豆成人午夜福利视频| 日韩精品有码人妻一区| av专区在线播放| 99久久无色码亚洲精品果冻| 久久精品人妻少妇| 日本三级黄在线观看| 国产精品爽爽va在线观看网站| 免费av不卡在线播放| 韩国av在线不卡| 我的女老师完整版在线观看| 亚洲性久久影院| 亚洲四区av| 国产精品无大码| 天堂网av新在线| av视频在线观看入口| 国产免费一级a男人的天堂| 成熟少妇高潮喷水视频| 琪琪午夜伦伦电影理论片6080| 亚洲经典国产精华液单| 很黄的视频免费| 国产精品免费一区二区三区在线| 久久久久久久久久成人| 在现免费观看毛片| 国产不卡一卡二| 国产在线男女| 18禁裸乳无遮挡免费网站照片| 亚洲乱码一区二区免费版| 精品久久久噜噜| 精品免费久久久久久久清纯| 老熟妇仑乱视频hdxx| 国产久久久一区二区三区| 亚洲最大成人av| 免费观看在线日韩| 特大巨黑吊av在线直播| 久久精品综合一区二区三区| 又紧又爽又黄一区二区| 免费看av在线观看网站| АⅤ资源中文在线天堂| 欧美不卡视频在线免费观看| 日韩 亚洲 欧美在线| 亚洲av.av天堂| 日日摸夜夜添夜夜添av毛片 | 大型黄色视频在线免费观看| 亚洲精品乱码久久久v下载方式| 色综合站精品国产| 中文亚洲av片在线观看爽| 日本撒尿小便嘘嘘汇集6| 深夜精品福利| 熟妇人妻久久中文字幕3abv| 99精品久久久久人妻精品| 久久精品夜夜夜夜夜久久蜜豆| 免费av观看视频| www.www免费av| 美女免费视频网站| 日韩大尺度精品在线看网址| 国产精品久久久久久av不卡| 日韩欧美国产一区二区入口| 国产在线精品亚洲第一网站| 日本爱情动作片www.在线观看 | 国产精品精品国产色婷婷| 在线a可以看的网站| 色视频www国产| 国产精品一区二区三区四区久久| 亚洲色图av天堂| 少妇人妻精品综合一区二区 | 成人毛片a级毛片在线播放| 国内久久婷婷六月综合欲色啪| 噜噜噜噜噜久久久久久91| 亚洲中文字幕一区二区三区有码在线看| 亚洲欧美精品综合久久99| 露出奶头的视频| 亚洲中文字幕一区二区三区有码在线看| 精品99又大又爽又粗少妇毛片 | 国产老妇女一区| 老司机午夜福利在线观看视频| 久久精品综合一区二区三区| 久久精品91蜜桃|