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

    新型水輪機(jī)聚能導(dǎo)流結(jié)構(gòu)的設(shè)計(jì)與參數(shù)優(yōu)化

    2017-08-01 00:02:51李龍敬周慎杰
    船舶力學(xué) 2017年7期
    關(guān)鍵詞:凸面凹面彎度

    李龍敬,周慎杰

    (山東大學(xué) 機(jī)械工程學(xué)院,山東 濟(jì)南 250061)

    新型水輪機(jī)聚能導(dǎo)流結(jié)構(gòu)的設(shè)計(jì)與參數(shù)優(yōu)化

    李龍敬,周慎杰

    (山東大學(xué) 機(jī)械工程學(xué)院,山東 濟(jì)南 250061)

    文章以提高葉片姿勢(shì)可變的水輪機(jī)獲能能力為目的,結(jié)合水輪機(jī)的運(yùn)動(dòng)規(guī)律,提出了一種新的三段式聚能導(dǎo)流結(jié)構(gòu),借助數(shù)值模擬的方法,對(duì)不同參數(shù)的NACA翼型的增速效果進(jìn)行了研究,確定了聚能導(dǎo)流結(jié)構(gòu)的具體結(jié)構(gòu)形式。通過改變影響聚能導(dǎo)流結(jié)構(gòu)水動(dòng)力學(xué)性能的幾個(gè)重要參數(shù),綜合實(shí)際應(yīng)用得到了水動(dòng)力學(xué)性能較優(yōu)的聚能導(dǎo)流結(jié)構(gòu),確定了導(dǎo)流結(jié)構(gòu)的安裝參數(shù)和結(jié)構(gòu)參數(shù)。

    翼型;聚能導(dǎo)流結(jié)構(gòu);增速比

    0 引 言

    通常潮流能水輪機(jī)的獲能效率在25%到30%之間[1],相對(duì)風(fēng)力發(fā)電,整體獲能效率較低,為了最大限度地獲取潮流的動(dòng)能,在潮流能獲能機(jī)構(gòu)外部增加聚能導(dǎo)流結(jié)構(gòu)是提高獲能效率的有效方法之一,同時(shí)導(dǎo)流結(jié)構(gòu)還能起到支撐體的作用,將水輪機(jī)固定在工作海域內(nèi)。因此,國內(nèi)外不少學(xué)者和科研機(jī)構(gòu)都在致力于聚能導(dǎo)流結(jié)構(gòu)的研究[2-3]。

    陳晗[4]在對(duì)彈簧控角豎軸直葉片水輪機(jī)的實(shí)驗(yàn)研究中,提出了凹形聚能導(dǎo)流結(jié)構(gòu)和凸形聚能導(dǎo)流結(jié)構(gòu)兩種結(jié)構(gòu),通過實(shí)驗(yàn)研究發(fā)現(xiàn)凹形聚能導(dǎo)流結(jié)構(gòu)水輪機(jī)的獲能效率要高于凸形聚能導(dǎo)流結(jié)構(gòu)的水輪機(jī)。孫科[5,8-9]借鑒風(fēng)力機(jī)擴(kuò)散器的理論,設(shè)計(jì)了一種對(duì)稱式的導(dǎo)流結(jié)構(gòu),提出了S型、W型等多種線型結(jié)構(gòu),借助數(shù)值模擬的方法對(duì)不同的線型結(jié)構(gòu)的性能進(jìn)行了分析;提出了設(shè)計(jì)聚能導(dǎo)流結(jié)構(gòu)擴(kuò)張門方案,通過研究擴(kuò)張門不僅可以提高罩內(nèi)流速,而且還有利于提高流場(chǎng)內(nèi)的穩(wěn)定性。劉斌[6]通過非定常渦面元法對(duì)加裝聚能導(dǎo)流結(jié)構(gòu)水輪機(jī)的水動(dòng)力性能進(jìn)行理論計(jì)算,分析了水動(dòng)力機(jī)理,得出聚能導(dǎo)流結(jié)構(gòu)幾何參數(shù)和安裝位置對(duì)于輪機(jī)水動(dòng)力性能的影響規(guī)律。王樹杰[7]針對(duì)柔性葉片的水輪機(jī)設(shè)計(jì)了一種中心對(duì)稱式的導(dǎo)流結(jié)構(gòu),通過改變影響聚能導(dǎo)流結(jié)構(gòu)性能的幾個(gè)重要參數(shù),對(duì)聚能導(dǎo)流結(jié)構(gòu)的水動(dòng)力學(xué)性能進(jìn)行了數(shù)值模擬,得出了水動(dòng)力學(xué)性能較優(yōu)的導(dǎo)流結(jié)構(gòu)。

    加拿大的Faure和Pratte[10]提出在河道內(nèi)垂直軸水輪機(jī)上安裝聚能導(dǎo)流結(jié)構(gòu)的想法,通過聚能導(dǎo)流結(jié)構(gòu)的引流作用改變水流沖擊水輪機(jī)的角度和速度,解決低流速時(shí)水輪機(jī)獲能效率偏低的問題。他們研究了聚能導(dǎo)流結(jié)構(gòu)入流角度、葉輪密實(shí)度、葉片與聚能導(dǎo)流結(jié)構(gòu)之間的縫隙等因素對(duì)水輪機(jī)獲能效率的影響。

    我們?cè)趯?duì)新型水輪機(jī)的前期研究中,設(shè)計(jì)了葉片姿勢(shì)可變的潮流能獲能水輪機(jī),并對(duì)水輪機(jī)的水動(dòng)力學(xué)性能進(jìn)行了研究,已經(jīng)取得了階段性的成果[11]。為了進(jìn)一步提高水輪機(jī)的獲能效率,根據(jù)水輪機(jī)自身的運(yùn)動(dòng)特點(diǎn),提出了一種新型的三段式導(dǎo)流結(jié)構(gòu),并對(duì)不同翼型導(dǎo)流結(jié)構(gòu)的增速效果進(jìn)行了研究,確定了增速段的NACA翼型結(jié)構(gòu),并針對(duì)不同的安裝參數(shù)和結(jié)構(gòu)參數(shù)進(jìn)行了數(shù)值計(jì)算,確定了聚能導(dǎo)流結(jié)構(gòu)的具體參數(shù)。

    1 聚能導(dǎo)流結(jié)構(gòu)的設(shè)計(jì)

    流體在流動(dòng)過程中沿流動(dòng)方向作用于固體壁面的總力稱為曳力,而固體壁面在流動(dòng)方向?qū)α黧w的反作用力稱為流體阻力。對(duì)于繞流流動(dòng),流體作用與固體壁面的力一般可分為正壓力和切應(yīng)力兩部分。相應(yīng)地,固體壁面的反作用力也由形狀阻力和摩擦阻力兩部分構(gòu)成。其中,形狀阻力是固體壁面上正壓力分布不均所產(chǎn)生的,又稱壓差阻力;而摩擦阻力則是固體壁面上切應(yīng)力分布不均所產(chǎn)生的。因此在導(dǎo)流結(jié)構(gòu)的設(shè)計(jì)當(dāng)中要盡可能減小流體阻力,從而最大限度地獲取潮流的動(dòng)能。

    1.1 聚能導(dǎo)流結(jié)構(gòu)的確定

    如圖1所示,根據(jù)水輪機(jī)自身運(yùn)動(dòng)規(guī)律,聚能導(dǎo)流結(jié)構(gòu)的設(shè)計(jì)應(yīng)該保證兩側(cè)來流時(shí)都有明顯的增速效果,且左側(cè)來流時(shí),葉片從正上方打開,水輪機(jī)的主增速區(qū)域應(yīng)處于轉(zhuǎn)子正上方,根據(jù)流體力學(xué)的理論,當(dāng)流體流過曲面物體時(shí),沿曲面方向的速度是改變的,所以考慮采用凹面的結(jié)構(gòu)會(huì)使流體沿切線方向向上運(yùn)動(dòng),使轉(zhuǎn)子上方的增速效果明顯。當(dāng)從右側(cè)來流時(shí),葉片從正下方打開,水輪機(jī)的主增速區(qū)域應(yīng)處于轉(zhuǎn)子正下方,同樣基于流體力學(xué)的理論,采用凸面結(jié)構(gòu)使流體沿著曲面的切線方向流動(dòng),從而達(dá)到了轉(zhuǎn)子下放流速增大的目的。

    圖1 聚能導(dǎo)流結(jié)構(gòu)和水輪機(jī)示意圖Fig.1 The diagrammatic sketch of diffuser and turbine

    基于上述考慮,擬采用單側(cè)三段式結(jié)構(gòu),根據(jù)每段的特點(diǎn),分別命名為凹面增速段、過渡段和凸面增速段。

    人們?cè)谝郧暗难芯恐?,凹面增速段大多采用圓弧結(jié)構(gòu)。由于圓弧的變化形式單一,導(dǎo)流后對(duì)水流的增速作用不夠強(qiáng),不便于最大限度地獲取潮流的動(dòng)能。為進(jìn)一步提高導(dǎo)流結(jié)構(gòu)的增速能力,我們借鑒NACA翼型結(jié)構(gòu),提出了在凹面增速段采用NACA翼型導(dǎo)流結(jié)構(gòu),并與之前的圓弧導(dǎo)流結(jié)構(gòu)進(jìn)行對(duì)比分析。

    由于凹面增速段與凸面增速段無法直接相連,故考慮采用圓弧的過渡結(jié)構(gòu)連接兩部分,即過渡段。

    在凸面增速段的結(jié)構(gòu)設(shè)計(jì)上,結(jié)合流體力學(xué)知識(shí),當(dāng)流體繞彎曲壁面流動(dòng)時(shí),為了盡可能地減小流動(dòng)阻力,需保證過渡段與凸面增速段相切,故凸面增速段優(yōu)先考慮圓弧結(jié)構(gòu)。

    NACA翼型是美國國家航空咨詢委員會(huì)(NACA)開發(fā)的一系列翼型。每個(gè)翼型的代號(hào)由“NACA”這四個(gè)字母與一串?dāng)?shù)字組成,將這串?dāng)?shù)字所描述的幾何參數(shù)代入特定方程中即可得到翼型的精確形狀。NACA四位數(shù)字翼型是美國NACA最早建立的一個(gè)低速翼型系列,與早期的其他翼型相比,有較高的最大升力系數(shù)和較低的阻力系數(shù)。在每個(gè)翼型代號(hào)NACA XYZZ中,后面四位數(shù)字的含義是:X—相對(duì)彎度;Y—最大彎度位置;ZZ—相對(duì)厚度。例如,NACA2412表示翼型的相對(duì)彎度為2%,最大彎度位置在弦長的0.4,相對(duì)厚度為12%。主要的翼型結(jié)構(gòu)如表1所示。

    表1 不同的翼型結(jié)構(gòu)Tab.1 Different kinds of airfoil profile

    本文分別對(duì)不同相對(duì)彎度、不同的最大彎度位置和不同的相對(duì)厚度進(jìn)行研究,然后取最優(yōu)結(jié)果與圓弧型迎流面進(jìn)行比較。

    1.2 主要參數(shù)

    聚能導(dǎo)流結(jié)構(gòu)的參數(shù)主要分為安裝參數(shù)和結(jié)構(gòu)參數(shù)兩部分。安裝參數(shù)主要包括聚能導(dǎo)流結(jié)構(gòu)在高度方向上與水輪機(jī)轉(zhuǎn)子中心的相對(duì)位置h和水平方向上與水輪機(jī)轉(zhuǎn)子中心的相對(duì)位置e;結(jié)構(gòu)參數(shù)根據(jù)實(shí)際情況主要考慮凹面增速段高度H和凸面增速段圓弧半徑R。

    凹面增速段高度H:由于增速段采用翼型結(jié)構(gòu),翼型的高度與長度是成對(duì)應(yīng)關(guān)系的,翼型高度確定后,其長度也是確定的,所以對(duì)于增速段只需要考慮翼型的相對(duì)厚度即可確定增速段的結(jié)構(gòu)。

    凸面增速段圓弧半徑R:對(duì)于凸面增速段結(jié)構(gòu),首先從結(jié)構(gòu)設(shè)計(jì)角度凸面增速段的圓弧要與過渡段的圓弧滿足相切的關(guān)系,所以凸面增速段的圓心相對(duì)位置是確定的,由于增速段高度已經(jīng)確定,所以凸面增速段結(jié)構(gòu)的高度也確定,通過改變不同的圓弧半徑,即可得到不同的凸面增速段的長度,可通過研究不同的圓弧半徑下聚能導(dǎo)流結(jié)構(gòu)的增速效果來確定凸面增速段圓弧半徑。

    2 水動(dòng)力學(xué)特性的數(shù)值模擬

    從優(yōu)化導(dǎo)流結(jié)構(gòu)來說,水輪機(jī)的轉(zhuǎn)動(dòng)并不會(huì)對(duì)增速效果有很大的影響,反而是水輪機(jī)轉(zhuǎn)動(dòng)的網(wǎng)格處理與計(jì)算量加大了模擬的復(fù)雜程度與難度,所以在模擬過程中考慮放置一靜止?fàn)顟B(tài)的水輪機(jī)對(duì)不同形式的導(dǎo)流結(jié)構(gòu)的增速效果進(jìn)行對(duì)比,選出最優(yōu)的導(dǎo)流結(jié)構(gòu)。

    2.1 數(shù)值模擬方法

    基于Fluent軟件穩(wěn)態(tài)計(jì)算方法,進(jìn)行聚能導(dǎo)流結(jié)構(gòu)在敞水狀態(tài)下的水動(dòng)力性能模擬。如圖2所示。計(jì)算區(qū)域:流場(chǎng)計(jì)算域以轉(zhuǎn)子中心為基準(zhǔn),轉(zhuǎn)子直徑2R為100mm,設(shè)置前方區(qū)域10倍于轉(zhuǎn)子直徑,后方10倍于轉(zhuǎn)子直徑,上取6倍于轉(zhuǎn)子直徑的矩形形區(qū)域,下到聚能導(dǎo)流結(jié)構(gòu)的下側(cè)邊緣。

    網(wǎng)格劃分:整個(gè)流場(chǎng)區(qū)域采用結(jié)構(gòu)化網(wǎng)格結(jié)構(gòu),在靠近聚能導(dǎo)流結(jié)構(gòu)的壁面區(qū)域增加網(wǎng)格密度,而在聚能導(dǎo)流結(jié)構(gòu)外圍則適當(dāng)放寬網(wǎng)格尺寸,以減小數(shù)值模擬的運(yùn)算量。整個(gè)計(jì)算區(qū)域的網(wǎng)格數(shù)量為20 989。

    邊界條件:定義為左側(cè)為流場(chǎng)入口,入口速度為1 m/s,邊界類型為速度入口(velocity-inlet);右側(cè)為流場(chǎng)出口,邊界類型為壓力出口(pressure-outlet);流場(chǎng)上下邊界、聚能導(dǎo)流結(jié)構(gòu)和轉(zhuǎn)子的壁面均采用無滑移壁面(wall)。

    針對(duì)聚能導(dǎo)流結(jié)構(gòu)流線型較好的特點(diǎn),數(shù)值仿真中選用k-ω湍流模型;時(shí)間步長為0.01s,迭代步數(shù)為1 500。

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

    在模擬過程中,網(wǎng)格的疏密程度會(huì)直接影響計(jì)算結(jié)果的準(zhǔn)確性。分別對(duì)導(dǎo)流罩內(nèi)流場(chǎng)進(jìn)行不同疏密程度的網(wǎng)格劃分,得到的網(wǎng)格數(shù)分別為10 388、20 989、30 660和43 311。

    結(jié)果如表2所示,將網(wǎng)格數(shù)為1萬與2萬的相比其改變量大于3%,而2萬與3萬和4萬相比,改變量均小于3%??紤]到網(wǎng)格數(shù)越多,計(jì)算時(shí)間越長,為了提高計(jì)算效率,采用了2萬網(wǎng)格的加密方式對(duì)流場(chǎng)區(qū)域進(jìn)行網(wǎng)格劃分。

    圖2 模擬流場(chǎng)模型Fig.2 Fluid domain model

    表2 不同單元尺寸下的增速比Tab.2 Speed increasing ratio under different element size

    2.3 模擬有效性的驗(yàn)證

    利用數(shù)值模擬的計(jì)算方法對(duì)流體結(jié)構(gòu)進(jìn)行設(shè)計(jì)與參數(shù)優(yōu)化是目前國內(nèi)外流體力學(xué)領(lǐng)域普遍采取的一種基本方法,計(jì)算結(jié)果也大都在后續(xù)的實(shí)驗(yàn)中得到了有效的驗(yàn)證。在前面引言中已經(jīng)提到的陳晗、孫科、王樹杰等均采用數(shù)值模擬的方法對(duì)不同類型的導(dǎo)流結(jié)構(gòu)進(jìn)行了研究,取得了一定的成果。本文也采用相類似的數(shù)值模擬方法,對(duì)新提出的三段式聚能導(dǎo)流結(jié)構(gòu)進(jìn)行了參數(shù)優(yōu)化。

    孫科[5]在論文中結(jié)合對(duì)水輪機(jī)的實(shí)驗(yàn)數(shù)據(jù)對(duì)幾種導(dǎo)流結(jié)構(gòu)模型進(jìn)行了數(shù)值模擬,計(jì)算結(jié)果已得到業(yè)內(nèi)專家的廣泛認(rèn)可。采用與本文數(shù)值研究相同的模擬方法,對(duì)孫科提出的有關(guān)導(dǎo)流結(jié)構(gòu)的3種模型進(jìn)行了數(shù)值計(jì)算,并與孫科論文中的原模擬結(jié)果進(jìn)行了對(duì)比,本次模擬結(jié)果與孫科論文中的模擬結(jié)果相吻合,誤差率均在3%以內(nèi)。由此可見,本文所采用的模擬方法是有效的。詳見表3。

    表3 導(dǎo)流罩模型內(nèi)的最大流速(m/s)Tab.3 The maximum velocity of dome models(m/s)

    3 數(shù)值結(jié)果分析

    考慮通過引入增速比來評(píng)定聚能導(dǎo)流結(jié)構(gòu)的增速效果。

    由于水流從兩側(cè)來流時(shí)主要的增速區(qū)域有一定的區(qū)別,所以以轉(zhuǎn)子中心為原點(diǎn),左側(cè)來流時(shí),取轉(zhuǎn)子上方55 mm到95 mm的范圍得到該區(qū)間的瞬時(shí)速度,然后對(duì)速度取平均值,計(jì)算出轉(zhuǎn)子上方的瞬時(shí)平均流速,然后再除以來流速度,得到水輪機(jī)上方的增速比。取轉(zhuǎn)子下方55 mm到75 mm的范圍得到該區(qū)間的瞬時(shí)速度,然后對(duì)速度取平均值,計(jì)算出轉(zhuǎn)子下方的瞬時(shí)平均流速,然后再除以來流速度,得到水輪機(jī)下方的增速比。右側(cè)來流時(shí),取轉(zhuǎn)子上方55 mm到75 mm的范圍得到該區(qū)間的瞬時(shí)速度,然后對(duì)速度取平均值,計(jì)算出轉(zhuǎn)子上方的瞬時(shí)平均流速,然后再除以來流速度,得到水輪機(jī)上方的增速比。取轉(zhuǎn)子下方55 mm到95 mm的范圍得到該區(qū)間的瞬時(shí)速度,然后對(duì)速度取平均值,計(jì)算出轉(zhuǎn)子下方的瞬時(shí)平均流速,然后再除以來流速度,得到水輪機(jī)下方的增速比。

    本文通過對(duì)不同翼型的聚能導(dǎo)流結(jié)構(gòu)進(jìn)行數(shù)值模擬,結(jié)合數(shù)值模擬結(jié)果分析確定最優(yōu)的翼型結(jié)構(gòu);在此基礎(chǔ)上再對(duì)聚能導(dǎo)流結(jié)構(gòu)進(jìn)行參數(shù)優(yōu)化,確定最優(yōu)的安裝參數(shù)和結(jié)構(gòu)參數(shù)。

    3.1 凹面增速段結(jié)構(gòu)的確定

    3.1.1 不同參數(shù)NACA翼型結(jié)構(gòu)的數(shù)值模擬結(jié)果

    在模擬中固定增速段高度為2.2R,圖3和圖4是聚能導(dǎo)流結(jié)構(gòu)在不同相對(duì)彎度(NACA2215、NACA4215、NACA5215和NACA6215)下的增速效果。隨著相對(duì)彎度的增大,充分利用翼型曲面對(duì)流場(chǎng)的干擾,改變了流體的運(yùn)動(dòng)軌跡,起到了更優(yōu)的引流作用,使得轉(zhuǎn)子正上方的流速進(jìn)一步增大,當(dāng)增加到一定程度后,彎度過大會(huì)導(dǎo)致聚能導(dǎo)流結(jié)構(gòu)的迎流面對(duì)流體的運(yùn)動(dòng)形成阻礙作用,所以聚能導(dǎo)流結(jié)構(gòu)的增速比呈現(xiàn)先大后小的過程,故NACA4215的翼型結(jié)構(gòu)增速效果更佳。根據(jù)翼型的代號(hào)規(guī)則確定翼型相對(duì)彎度為4。

    圖3 不同相對(duì)彎度下左側(cè)來流時(shí)聚能導(dǎo)流結(jié)構(gòu)上方增速比Fig.3 Speed increasing ratio above diffuser under different relative curvature in case of incoming flow on the left

    圖4不同相對(duì)彎度下左側(cè)來流時(shí)聚能導(dǎo)流結(jié)構(gòu)下方增速比Fig.4 Speed increasing ratio below diffuser under different relative curvature in case of incoming flow on the left

    圖5 和圖6是不同相對(duì)厚度(NACA4206、NACA4212、NACA4215和NACA4221)下聚能導(dǎo)流結(jié)構(gòu)的增速效果,從圖中可以看出,隨著相對(duì)厚度的增加,針對(duì)固定高度的導(dǎo)流結(jié)構(gòu),相當(dāng)于增速段的長度由長變短,當(dāng)長度過長時(shí),阻力形式主要為摩擦阻力,但過于平緩增速效果不佳,隨著長度變短,阻力形式以壓差阻力為主,對(duì)流體的阻礙作用慢慢增大,聚能導(dǎo)流結(jié)構(gòu)的增速效果呈現(xiàn)了迅速增加到快速下降然后趨于平緩的過程,對(duì)比可以看出NACA4212翼型結(jié)構(gòu)的增速效果更佳。根據(jù)翼型的代號(hào)規(guī)則確定翼型相對(duì)厚度為12。

    圖5 不同翼型厚度下左側(cè)來流時(shí)聚能導(dǎo)流結(jié)構(gòu)上方增速比Fig.5 Speed increasing ratio above diffuser under different airfoilthickness in case ofincoming flow on the left

    圖6 不同翼型厚度下左側(cè)來流時(shí)聚能導(dǎo)流結(jié)構(gòu)下方增速比Fig.6 Speed increasing ratio below diffuser under differentairfoilthickness in case ofincoming flow on the left

    圖7 不同最大彎度下左側(cè)來流時(shí)聚能導(dǎo)流結(jié)構(gòu)上方增速比Fig.7 Speed increasing ratio above diffuser under different circle radius in case of incoming flow on the left

    圖8不同最大彎度下左側(cè)來流時(shí)聚能導(dǎo)流結(jié)構(gòu)下方增速比Fig.8 Speed increasing ratio below diffuser under different circle radius in case of incoming flow on the left

    圖7 和圖8是不同最大彎度(NACA4212、NACA4312、NACA4412和NACA4512)下聚能導(dǎo)流結(jié)構(gòu)的增速效果。從圖中可以看出隨著最大彎度的變大,相當(dāng)于改變了翼型曲面的彎曲程度,會(huì)改變曲面對(duì)流場(chǎng)的干擾作用,引起轉(zhuǎn)子上方流速的變化,增速比呈現(xiàn)先下降后逐漸增大的趨勢(shì),聚能導(dǎo)流結(jié)構(gòu)下方的增速比呈現(xiàn)先增大后逐漸下降的趨勢(shì),在這個(gè)變化過程中聚能導(dǎo)流結(jié)構(gòu)的主作用增速區(qū)域發(fā)生了變換,綜合聚能導(dǎo)流結(jié)構(gòu)上下的增速比考慮采用NACA4212的翼型結(jié)構(gòu)。根據(jù)翼型的代號(hào)規(guī)則確定翼型最大彎度位置為2。

    3.1.2 圓弧結(jié)構(gòu)的數(shù)值模擬結(jié)果

    給定相同的凹面增速段高度與長度參數(shù),對(duì)采用不同圓弧半徑的導(dǎo)流結(jié)構(gòu)進(jìn)行數(shù)值模擬,根據(jù)圖9和圖10,可以得到在半徑為300時(shí),聚能導(dǎo)流結(jié)構(gòu)有更佳的增速效果,導(dǎo)流結(jié)構(gòu)上方增速比為1.531,下方增速比1.518。

    圖9 不同圓弧半徑下左側(cè)來流時(shí)聚能導(dǎo)流結(jié)構(gòu)上方增速比Fig.9 Speed increasing ratio above diffuser under different airfoil thickness in case of incoming flow on the left

    圖10 不同圓弧半徑下左側(cè)來流時(shí)聚能導(dǎo)流結(jié)構(gòu)下方增速比Fig.10 Speed increasing ratio below diffuser under different airfoil thickness in case of incoming flow on the left

    3.1.3 翼型凹面增速段與圓弧型凹面增速段的增速分析

    由表4可以看出,NACA4212翼型結(jié)構(gòu)的增速效果優(yōu)于圓弧結(jié)構(gòu),故增速段結(jié)構(gòu)應(yīng)采用NACA4212的翼型結(jié)構(gòu)。

    表4 不同線型的增速比Tab.4 Speed increasing ratio under different linetypes

    3.2 安裝參數(shù)—長度方向上不同的安裝距離e

    安裝參數(shù)主要考慮聚能導(dǎo)流結(jié)構(gòu)在高度方向上和長度方向上與水輪機(jī)的相對(duì)位置??紤]讓聚能導(dǎo)流結(jié)構(gòu)的內(nèi)壁與水輪機(jī)葉片完全張開時(shí)的軌跡相切,故聚能導(dǎo)流結(jié)構(gòu)在高度方向上與水輪機(jī)的相對(duì)位置可以直接確定,在水平方向上的距離e要通過模擬得到。圖11為聚能導(dǎo)流結(jié)構(gòu)安裝位置參數(shù)示意圖。

    取R=50,結(jié)合實(shí)際考慮取安裝距離分別為0.6R、0.8R、1R、1.2R、1.4R、1.6R和1.8R七種工況,采用控制變量法進(jìn)行研究。數(shù)值模擬結(jié)果如圖12-15所示。隨著聚能導(dǎo)流結(jié)構(gòu)在長度方向上的安裝距離e從大到小不斷變化,聚能導(dǎo)流結(jié)構(gòu)上下方的平均增速比呈現(xiàn)先增大后減小的過程。當(dāng)聚能導(dǎo)流結(jié)構(gòu)的頂點(diǎn)位置距離水輪機(jī)轉(zhuǎn)子中心位置一個(gè)R的距離時(shí),平均增速比達(dá)到最大值。隨著安裝位置的移動(dòng),聚能導(dǎo)流結(jié)構(gòu)下放的增速比變化趨勢(shì)更為明顯,當(dāng)聚能導(dǎo)流結(jié)構(gòu)逐漸靠近水輪機(jī)轉(zhuǎn)子后,轉(zhuǎn)子上下均有較好的增速效果。

    圖11 聚能導(dǎo)流結(jié)構(gòu)的安裝位置參數(shù)Fig.11 The installation parameters ofdiffuser

    圖12 不同安裝位置e下聚能導(dǎo)流結(jié)構(gòu)上方的增速比Fig.12 Speed increasing ratio under different installation sites

    圖15 安裝位置為1.8R時(shí)的速度矢量圖Fig.15 The velocity vector when installed at location of 1.8R

    3.3 結(jié)構(gòu)參數(shù)

    結(jié)構(gòu)參數(shù)如圖16所示,根據(jù)翼型結(jié)構(gòu)的特點(diǎn)以及設(shè)計(jì)的原則,主要參數(shù)包括凹面增速段高度H和凸面增速段圓弧半徑r。

    3.3.1 不同高度H下聚能導(dǎo)流結(jié)構(gòu)的增速效果分析

    考慮研究不同高度下的NACA4212的結(jié)構(gòu),探究增速段高度對(duì)導(dǎo)流增速效果的影響。取2.2R、2.6R、3R和3.4R四種情況進(jìn)行分析。模擬結(jié)果如圖17和圖18所示。隨著聚能導(dǎo)流結(jié)構(gòu)高度的增加,聚能導(dǎo)流結(jié)構(gòu)對(duì)流場(chǎng)的干擾作用發(fā)生了改變,對(duì)水流的阻礙作用由摩擦阻力變成壓差阻力,聚能導(dǎo)流結(jié)構(gòu)上方增速比會(huì)逐漸增大然后慢慢減小并趨于平緩;而聚能導(dǎo)流結(jié)構(gòu)下放增速比一直處于逐漸增大的過程,轉(zhuǎn)子下方流速明顯增大,但轉(zhuǎn)子上方的流速相對(duì)減小,這是聚能導(dǎo)流結(jié)構(gòu)整體結(jié)構(gòu)的增大與轉(zhuǎn)子相對(duì)大小綜合作用的結(jié)果。綜合考慮增速效果與成本,考慮選用高為2.6R的導(dǎo)流結(jié)構(gòu)。

    3.3.2 不同半徑r聚能導(dǎo)流結(jié)構(gòu)的增速效果分析

    考慮水流從右側(cè)流過時(shí),研究圓弧半徑對(duì)聚能導(dǎo)流結(jié)構(gòu)增速效果的影響。分別取4.5R、5R、5.5R、6R、8R和10R六種不同的圓弧半徑,模擬結(jié)果如圖21所示,隨著半徑的不斷增大,聚能導(dǎo)流結(jié)構(gòu)本身對(duì)于水流的阻礙作用會(huì)相對(duì)減小,增速比會(huì)逐漸增大,但是當(dāng)聚能導(dǎo)流結(jié)構(gòu)凸面增速段圓弧半徑增大到6R后,隨著半徑的再度增大,增速效果會(huì)逐漸減小然后趨于平緩,這是因?yàn)殡S著半徑的進(jìn)一步增大,聚能導(dǎo)流結(jié)構(gòu)的凸面增速段越接近于直線結(jié)構(gòu),圓弧對(duì)流場(chǎng)幾乎起不到干擾作用,增速效果也逐漸趨近于直線結(jié)構(gòu)。

    圖16 聚能導(dǎo)流結(jié)構(gòu)的結(jié)構(gòu)參數(shù)Fig.16 The structure parameters of diffuser

    圖17 不同高度H下聚能導(dǎo)流結(jié)構(gòu)上方的增速比Fig.17 Speed increasing ratio above diffuser under different H

    圖18 不同高度H下聚能導(dǎo)流結(jié)構(gòu)下方的增速比Fig.18 Speed increasing ratio below diffuser under different H

    圖19 高度為2.6R時(shí)速度矢量圖Fig.19 Velocity vector atthe heightof2.6R

    圖20 高度為3.4R時(shí)速度矢量圖Fig.20 Velocity vector at the heightof3.4R

    4 結(jié) 論

    本文根據(jù)新型水輪機(jī)運(yùn)動(dòng)特點(diǎn),以新型水輪機(jī)聚能導(dǎo)流結(jié)構(gòu)為研究對(duì)象,提出了新型三段式聚能導(dǎo)流結(jié)構(gòu),通過數(shù)值模擬的方法確定了增速段NACA翼型結(jié)構(gòu)選型,對(duì)影響聚能導(dǎo)流結(jié)構(gòu)的主要參數(shù)進(jìn)行了優(yōu)化研究,得出以下結(jié)論:

    (1)本文首次對(duì)NACA翼型結(jié)構(gòu)的增速效果進(jìn)行了研究,通過對(duì)不同參數(shù)的翼型結(jié)構(gòu)進(jìn)行數(shù)值計(jì)算,確定選用NACA4212的翼型作為凹面增速段結(jié)構(gòu)。根據(jù)模擬結(jié)果,NACA4212的翼型結(jié)構(gòu)聚能導(dǎo)流結(jié)構(gòu)上方增速比為1.591,下方增速比為1.581,比同樣條件下采用圓弧型聚能導(dǎo)流結(jié)構(gòu)的增速比分別提高了3.92%和4.15%。

    (2)在增速段選用翼型結(jié)構(gòu)的基礎(chǔ)上進(jìn)行結(jié)構(gòu)參數(shù)優(yōu)化,根據(jù)翼型結(jié)構(gòu)的特點(diǎn),增速段的高度與長度成比例關(guān)系,故選取不同的增速段高度作為優(yōu)化參數(shù)進(jìn)行優(yōu)化,最終得到聚能導(dǎo)流結(jié)構(gòu)的凹面增速段的高度為2.6R;考慮到水流從右側(cè)來流的增速效果,針對(duì)右側(cè)凸面增速段圓弧確定的結(jié)構(gòu)特點(diǎn),選取了圓弧半徑作為優(yōu)化參數(shù),最終確定凸面增速段的圓弧半徑為6R。

    (3)根據(jù)水輪機(jī)的運(yùn)動(dòng)特點(diǎn),考慮讓聚能導(dǎo)流結(jié)構(gòu)的內(nèi)壁與水輪機(jī)葉片完全張開時(shí)的軌跡盡量相切,故聚能導(dǎo)流結(jié)構(gòu)在高度方向上與水輪機(jī)的相對(duì)位置可以直接確定,因此選取了豎直方向上的距離e作為優(yōu)化參數(shù)進(jìn)行優(yōu)化,確定了聚能導(dǎo)流結(jié)構(gòu)凹面增速段與過渡段連接處與葉片姿勢(shì)可變水輪機(jī)轉(zhuǎn)子之間水平方向的距離為R。

    圖21 右側(cè)來流時(shí)不同半徑下聚能導(dǎo)流結(jié)構(gòu)下方的增速比Fig.21 Speed increasing ratio below diffuser under different radius in case of incoming flow on the right

    參考文獻(xiàn):

    [1]余 寅,唐宏德,郭家寶.中國可再生能源發(fā)展前景分析[J].華東電力,2009,37(8):1306-1038.Yu Yin,Tang Hongde,Guo Jiabao.Developing prospect analysis of renewable energy resources in China[J].East China Electric Power,2009,37(8):1306-1038.

    [2]Cresswell N W,Ingram G L,Dominy R G.The impact of diffuser augmentation on a tidal stream turbine[J].Ocean Engineering,2015,108(12):155-163.

    [3]Daniele E,Coiro D P.Optimization of diffuser geometry for an horizontal axis shrouded hydro turbine[C].International Conference on Clean Electrical Power.IEEE,2013:240-247.

    [4]陳 晗.彈簧控角豎軸直葉水輪機(jī)水動(dòng)力性能研究[D].哈爾濱:哈爾濱工程大學(xué),2006.Chen Han.On study of hydrodynamic performance of the spring-controlling of vertical-axis variable-pitch turbine[D]. Harbin:Harbin Engineering University,2006.

    [5]孫 科.豎軸H型葉輪及聚能導(dǎo)流結(jié)構(gòu)流體動(dòng)力性能參數(shù)模擬[D].哈爾濱:哈爾濱工程大學(xué),2008.Sun Ke.Numerical simulation on fluid dynamic performance of h-shaped vertical axis turbine and duct[D].Harbin:Harbin Engineering University,2008.

    [6]劉 斌.帶聚能導(dǎo)流結(jié)構(gòu)的豎軸潮流水輪機(jī)水動(dòng)力性能研究[D].哈爾濱:哈爾濱工程大學(xué),2007.Liu Bin.Study on the hydrodynamic performance of vertical-axis turbine with demos[D].Harbin:Harbin Engineering U-niversity,2007.

    [7]鹿蘭帥.潮流能驅(qū)動(dòng)的柔性葉片水輪機(jī)水動(dòng)力性能試驗(yàn)研究[D].青島:中國海洋大學(xué),2009.Lu Lanshuai.Experimental study on hydrodynamic performances of flexible blade turbine driven by tidal energy[D]. Qingdao:Ocean University of China,2009.

    [8]張 亮,孫 科,羅慶杰.潮流水輪機(jī)聚能導(dǎo)流結(jié)構(gòu)的水動(dòng)力設(shè)計(jì)[J].哈爾濱工程大學(xué)學(xué)報(bào),2007,28(7):734-737.Zhang Liang,Sun Ke,Luo Qingjie.Hydrodynamic design of diversion cover for a tidal-stream hydro turbine[J].Journal of Harbin Engineering University,2007,28(7):734-737.

    [9]荊豐梅,張 亮,張鵬遠(yuǎn),等.潮流能發(fā)電增速聚能導(dǎo)流結(jié)構(gòu)研究[J].哈爾濱工程大學(xué)學(xué)報(bào),2012,33(4):409-413.Jing Fengmei,Zhang Liang,Zhang Pengyuan.Study on an accelerated diffuser of tidal current power generation[J].Journal of Harbin Engineering University,2007,28(7):734-737.

    [10]Batten W M J,Bahaj A S,Molland A F,et al.Experimentally validated numerical method for the hydrodynamic design of horizontal axis tidal turbines[J].Ocean Engineering,2007,34(7):1013-1020.

    [11]Li L J,Zhou S J.Numerical simulation of hydrodynamic performance of blade position-variable hydraulic turbine[J].Journal of Hydrodynamics Ser B,2017,29(2):314-321.

    Design and parameter optimization of novel energy-gathered and diversion mechanism

    LI Long-jing,ZHOU Shen-jie
    (Shandong University School of Mechanical Engineering,Jinan 250061,China)

    In order to improve the efficiency of the Blade Position-Variable Hydraulic Turbine,a threepiece energy-gathered and diversion mechanism is proposed,combined with motion law of hydraulic turbine. By means of the numerical simulation method,the accelerating effect of different parameters of NACA airfoilis studied,then the specific structure of energy-gathered and diversion mechanism is determined.A kind of energy-gathered and diversion mechanism with excellent hydraulic performance is obtained through changing severalimportantparameters.

    airfoil;energy-gathered and diversion mechanism;speed increasing ratio

    TB126 TK79

    A

    10.3969/j.issn.1007-7294.2017.07.005

    1007-7294(2017)07-0832-10

    2017-03-08

    山東省科技開發(fā)項(xiàng)目(2014GGX103028)

    李龍敬(1988-),男,博士研究生,E-mail:kglilongjing@163.com;

    周慎杰(1958-),男,教授。

    猜你喜歡
    凸面凹面彎度
    活動(dòng)期油源斷裂輸導(dǎo)油氣有利凸面脊分布預(yù)測(cè)方法
    連續(xù)變彎度翼型動(dòng)態(tài)氣動(dòng)特性數(shù)值模擬
    全息離子束擺動(dòng)刻蝕凸面閃耀光柵制備技術(shù)
    翼型彎度變化自適應(yīng)規(guī)律研究
    搜集凹面錐體
    機(jī)翼后緣連續(xù)變彎度對(duì)客機(jī)氣動(dòng)特性影響
    減壓小心機(jī)
    決策探索(2016年21期)2016-11-28 09:30:57
    彩寶藝術(shù)大師 Bulgari
    鳳凰生活(2016年8期)2016-08-10 12:38:53
    生物型人工硬腦膜預(yù)防大腦凸面腦膜瘤切除術(shù)后癲癇的療效觀察
    紫外全息凹面光柵離子束刻蝕技術(shù)
    爱豆传媒免费全集在线观看| 免费看美女性在线毛片视频| 亚洲激情五月婷婷啪啪| 99热精品在线国产| 国产精品伦人一区二区| 色播亚洲综合网| 青春草亚洲视频在线观看| 亚洲熟妇中文字幕五十中出| 欧美日韩精品成人综合77777| 小说图片视频综合网站| 国产伦理片在线播放av一区| 又黄又爽又刺激的免费视频.| 最近2019中文字幕mv第一页| 免费观看性生交大片5| 99在线人妻在线中文字幕| 免费av毛片视频| 午夜a级毛片| 成人三级黄色视频| 亚洲精品,欧美精品| 一区二区三区高清视频在线| 在线a可以看的网站| 久久鲁丝午夜福利片| 赤兔流量卡办理| 天堂影院成人在线观看| 亚洲精品乱码久久久v下载方式| 国产人妻一区二区三区在| 国产精品久久视频播放| av女优亚洲男人天堂| 99国产精品一区二区蜜桃av| 中国国产av一级| 亚洲精品色激情综合| 看非洲黑人一级黄片| 久久99蜜桃精品久久| 日韩精品青青久久久久久| 精品国产一区二区三区久久久樱花 | 尾随美女入室| 大又大粗又爽又黄少妇毛片口| 高清在线视频一区二区三区 | 麻豆一二三区av精品| 久久久久久久亚洲中文字幕| 国产久久久一区二区三区| 婷婷色av中文字幕| a级毛色黄片| 在线免费观看不下载黄p国产| 99久国产av精品国产电影| 卡戴珊不雅视频在线播放| 欧美激情在线99| 亚洲欧美日韩卡通动漫| 狠狠狠狠99中文字幕| 一级二级三级毛片免费看| 国产精品一及| 男的添女的下面高潮视频| 看免费成人av毛片| 久久精品熟女亚洲av麻豆精品 | 国产色爽女视频免费观看| 一二三四中文在线观看免费高清| 欧美又色又爽又黄视频| 精品久久久久久久人妻蜜臀av| 国产精品熟女久久久久浪| 亚洲美女视频黄频| 日本熟妇午夜| 国产精品久久久久久久久免| 日本免费a在线| 97超视频在线观看视频| 日韩制服骚丝袜av| 日日撸夜夜添| 我的女老师完整版在线观看| 国产不卡一卡二| 久久久久精品久久久久真实原创| 黄片无遮挡物在线观看| 只有这里有精品99| 汤姆久久久久久久影院中文字幕 | 国产色婷婷99| videos熟女内射| 国产日韩欧美在线精品| 国产色婷婷99| 人人妻人人看人人澡| 国产精品国产高清国产av| 欧美日本视频| 国产伦精品一区二区三区四那| videos熟女内射| 又黄又爽又刺激的免费视频.| 亚洲美女搞黄在线观看| 蜜桃亚洲精品一区二区三区| 99热这里只有精品一区| 嘟嘟电影网在线观看| 欧美成人免费av一区二区三区| 国产精品蜜桃在线观看| 亚洲国产精品国产精品| 中文字幕人妻熟人妻熟丝袜美| 乱码一卡2卡4卡精品| 乱码一卡2卡4卡精品| 国产成人a区在线观看| 美女被艹到高潮喷水动态| 国产精品不卡视频一区二区| 伦精品一区二区三区| 国产精品国产高清国产av| av福利片在线观看| 国产一级毛片七仙女欲春2| 久久精品人妻少妇| 国产 一区精品| 哪个播放器可以免费观看大片| www.色视频.com| 日日撸夜夜添| 91久久精品国产一区二区三区| 成人午夜高清在线视频| 最近中文字幕高清免费大全6| 美女内射精品一级片tv| 午夜福利视频1000在线观看| 美女cb高潮喷水在线观看| 国内精品宾馆在线| 久久精品久久久久久久性| 亚洲av电影不卡..在线观看| 日本三级黄在线观看| 国产成人精品一,二区| 亚洲人成网站在线观看播放| 特级一级黄色大片| 日本免费一区二区三区高清不卡| 九九在线视频观看精品| 人体艺术视频欧美日本| 精品人妻偷拍中文字幕| 国产精品无大码| www.色视频.com| 亚洲av成人av| 有码 亚洲区| 一个人免费在线观看电影| 日韩高清综合在线| 国产成人aa在线观看| 亚洲一区高清亚洲精品| 久久精品国产鲁丝片午夜精品| 久久人妻av系列| 午夜视频国产福利| 国产探花极品一区二区| 夫妻性生交免费视频一级片| 成人毛片60女人毛片免费| 日本与韩国留学比较| 亚洲怡红院男人天堂| 久久久久久久久久黄片| 中文精品一卡2卡3卡4更新| 国产黄片视频在线免费观看| 国产精华一区二区三区| 少妇熟女欧美另类| 国产精品久久电影中文字幕| 久久人妻av系列| 欧美+日韩+精品| 一级黄色大片毛片| 纵有疾风起免费观看全集完整版 | 久热久热在线精品观看| 亚洲av电影不卡..在线观看| 又爽又黄a免费视频| 欧美变态另类bdsm刘玥| 亚洲性久久影院| 亚洲av日韩在线播放| 日本与韩国留学比较| 99热网站在线观看| 丰满人妻一区二区三区视频av| 麻豆国产97在线/欧美| 亚洲av一区综合| 精品免费久久久久久久清纯| 日本一二三区视频观看| 国产亚洲一区二区精品| 日韩欧美国产在线观看| 精品久久久久久久人妻蜜臀av| 久久精品国产亚洲av涩爱| 欧美最新免费一区二区三区| 日韩人妻高清精品专区| av免费观看日本| 久久这里只有精品中国| 亚洲一级一片aⅴ在线观看| 亚洲内射少妇av| 日本爱情动作片www.在线观看| 免费搜索国产男女视频| 国产淫片久久久久久久久| 久久久国产成人精品二区| 高清毛片免费看| 成人欧美大片| 成人高潮视频无遮挡免费网站| 精品不卡国产一区二区三区| 午夜日本视频在线| 久久婷婷人人爽人人干人人爱| 久久精品综合一区二区三区| 精品久久久久久久末码| 麻豆成人午夜福利视频| 中文字幕人妻熟人妻熟丝袜美| 一级毛片aaaaaa免费看小| 日韩国内少妇激情av| 日本-黄色视频高清免费观看| 国产成人精品久久久久久| 午夜福利在线在线| 99热网站在线观看| 国产亚洲最大av| 色播亚洲综合网| 亚洲成人久久爱视频| 欧美成人免费av一区二区三区| 看黄色毛片网站| 国产老妇女一区| 美女大奶头视频| 天堂√8在线中文| 春色校园在线视频观看| 桃色一区二区三区在线观看| 在线播放国产精品三级| 亚洲无线观看免费| 1000部很黄的大片| 亚洲综合色惰| 国产精品日韩av在线免费观看| 欧美成人a在线观看| 青春草国产在线视频| 免费观看精品视频网站| 三级经典国产精品| 亚洲真实伦在线观看| 少妇的逼好多水| 一级毛片我不卡| 国产91av在线免费观看| 中国美白少妇内射xxxbb| 精品熟女少妇av免费看| 免费看a级黄色片| 麻豆国产97在线/欧美| 日韩欧美精品免费久久| 老司机福利观看| 日韩欧美在线乱码| 欧美区成人在线视频| 午夜a级毛片| 日本免费a在线| 日韩亚洲欧美综合| 成年av动漫网址| 18+在线观看网站| 欧美三级亚洲精品| 日本wwww免费看| 搡老妇女老女人老熟妇| 中文天堂在线官网| 久久99热这里只有精品18| 最近中文字幕高清免费大全6| 国产麻豆成人av免费视频| 国产精品麻豆人妻色哟哟久久 | 成人美女网站在线观看视频| 国产精品一区二区在线观看99 | 国产精品一区二区三区四区免费观看| 男人舔奶头视频| 国产成人福利小说| 国产在视频线精品| 日本黄色片子视频| 久久亚洲国产成人精品v| 成人高潮视频无遮挡免费网站| a级一级毛片免费在线观看| 久久精品国产亚洲网站| 免费大片18禁| 免费播放大片免费观看视频在线观看 | 国产一区二区在线av高清观看| 日本与韩国留学比较| 中文精品一卡2卡3卡4更新| 亚洲欧美清纯卡通| 亚洲图色成人| 亚洲欧美日韩东京热| 自拍偷自拍亚洲精品老妇| 久久精品综合一区二区三区| 成人性生交大片免费视频hd| 亚洲在线自拍视频| 看十八女毛片水多多多| 91av网一区二区| 狂野欧美激情性xxxx在线观看| 简卡轻食公司| 亚洲三级黄色毛片| 少妇熟女aⅴ在线视频| 成年版毛片免费区| 亚洲乱码一区二区免费版| 91久久精品国产一区二区三区| 午夜亚洲福利在线播放| 性色avwww在线观看| 少妇熟女欧美另类| 韩国高清视频一区二区三区| 亚洲三级黄色毛片| 亚洲中文字幕一区二区三区有码在线看| 久久久久久久久大av| 卡戴珊不雅视频在线播放| 久久精品国产99精品国产亚洲性色| 啦啦啦韩国在线观看视频| 精品国产一区二区三区久久久樱花 | 亚洲国产色片| 高清在线视频一区二区三区 | 男人舔奶头视频| 亚洲内射少妇av| 国产精品久久久久久精品电影小说 | 午夜老司机福利剧场| 99久国产av精品| 亚洲在久久综合| av免费观看日本| 亚洲欧美日韩卡通动漫| 久久精品久久久久久久性| 91狼人影院| 在线观看一区二区三区| 久久久久久大精品| 春色校园在线视频观看| 国产精品永久免费网站| 岛国毛片在线播放| 婷婷六月久久综合丁香| 一个人观看的视频www高清免费观看| 久热久热在线精品观看| 国产伦一二天堂av在线观看| 亚洲精品自拍成人| 欧美性感艳星| 成人鲁丝片一二三区免费| 91狼人影院| 91午夜精品亚洲一区二区三区| 亚洲伊人久久精品综合 | 久久精品久久久久久久性| 尾随美女入室| 欧美日韩在线观看h| av在线天堂中文字幕| 免费在线观看成人毛片| 看黄色毛片网站| 久久99蜜桃精品久久| 综合色av麻豆| 乱系列少妇在线播放| 十八禁国产超污无遮挡网站| 国产熟女欧美一区二区| 26uuu在线亚洲综合色| 91精品伊人久久大香线蕉| 欧美高清性xxxxhd video| 内射极品少妇av片p| 精品人妻视频免费看| 久久久久久大精品| 非洲黑人性xxxx精品又粗又长| 国产成人午夜福利电影在线观看| 亚洲激情五月婷婷啪啪| 亚洲精品国产成人久久av| 日韩亚洲欧美综合| 久久这里有精品视频免费| 成人毛片a级毛片在线播放| 亚洲熟妇中文字幕五十中出| 99久国产av精品国产电影| 天堂av国产一区二区熟女人妻| 亚洲精品自拍成人| av在线天堂中文字幕| 欧美成人免费av一区二区三区| 免费无遮挡裸体视频| 久久久久久久久久黄片| 美女cb高潮喷水在线观看| 欧美一级a爱片免费观看看| 欧美三级亚洲精品| 丰满乱子伦码专区| 亚洲精品乱码久久久久久按摩| 亚洲熟妇中文字幕五十中出| 欧美xxxx性猛交bbbb| 高清视频免费观看一区二区 | 久久久亚洲精品成人影院| 一级毛片久久久久久久久女| 成人二区视频| 亚洲欧美日韩东京热| 久久久成人免费电影| av国产久精品久网站免费入址| 欧美另类亚洲清纯唯美| 波野结衣二区三区在线| 久久久久久久久久黄片| 亚洲欧美中文字幕日韩二区| 18+在线观看网站| 欧美日韩精品成人综合77777| 免费搜索国产男女视频| 亚洲人成网站在线播| 亚洲成人av在线免费| 亚洲熟妇中文字幕五十中出| 在线免费十八禁| 免费人成在线观看视频色| 身体一侧抽搐| 久久精品综合一区二区三区| 在线免费观看不下载黄p国产| 91精品一卡2卡3卡4卡| 欧美日韩在线观看h| 男女视频在线观看网站免费| 十八禁国产超污无遮挡网站| ponron亚洲| 舔av片在线| 亚洲精品乱码久久久v下载方式| 亚洲经典国产精华液单| 国产私拍福利视频在线观看| 大香蕉97超碰在线| 亚洲丝袜综合中文字幕| 精品人妻视频免费看| 好男人视频免费观看在线| 亚洲天堂国产精品一区在线| 亚洲中文字幕一区二区三区有码在线看| 一边亲一边摸免费视频| 日韩一区二区三区影片| 精品酒店卫生间| 日韩欧美三级三区| 在线免费十八禁| 国产激情偷乱视频一区二区| 日产精品乱码卡一卡2卡三| 日韩成人av中文字幕在线观看| 极品教师在线视频| av天堂中文字幕网| 久久久久久国产a免费观看| 女人被狂操c到高潮| 三级经典国产精品| 婷婷六月久久综合丁香| 麻豆成人午夜福利视频| 亚洲精华国产精华液的使用体验| 黄片无遮挡物在线观看| 嫩草影院入口| 亚洲人成网站高清观看| 边亲边吃奶的免费视频| 最近视频中文字幕2019在线8| 99久国产av精品| 九九在线视频观看精品| 长腿黑丝高跟| 国产中年淑女户外野战色| 国产探花极品一区二区| 久久这里只有精品中国| 国产精品嫩草影院av在线观看| 欧美变态另类bdsm刘玥| 亚洲欧美成人综合另类久久久 | 一级黄色大片毛片| 久久久久久久午夜电影| 国产精品蜜桃在线观看| 亚洲电影在线观看av| 狠狠狠狠99中文字幕| 真实男女啪啪啪动态图| 国产亚洲5aaaaa淫片| 欧美成人一区二区免费高清观看| 99久久精品一区二区三区| 欧美bdsm另类| 国产人妻一区二区三区在| 黄色一级大片看看| 国产国拍精品亚洲av在线观看| 18禁动态无遮挡网站| 最近中文字幕2019免费版| 久久久久久久午夜电影| 国产精品永久免费网站| 九草在线视频观看| 国产精品,欧美在线| 大香蕉97超碰在线| 国产精品不卡视频一区二区| 日本wwww免费看| 人妻系列 视频| av.在线天堂| 99九九线精品视频在线观看视频| 国产乱人偷精品视频| 国产三级在线视频| 美女cb高潮喷水在线观看| 国产探花极品一区二区| 又粗又硬又长又爽又黄的视频| 18禁动态无遮挡网站| 亚洲精品亚洲一区二区| 成人无遮挡网站| 久久久久久久久中文| 国产精品人妻久久久影院| 婷婷色av中文字幕| ponron亚洲| 精品久久久久久久末码| 看片在线看免费视频| 久久99热这里只有精品18| 国产精品一二三区在线看| 亚洲国产高清在线一区二区三| 视频中文字幕在线观看| 联通29元200g的流量卡| 国产一区二区在线观看日韩| 在线免费十八禁| 久久久久久久久久久免费av| 国产精品日韩av在线免费观看| 天堂影院成人在线观看| 婷婷色综合大香蕉| 国产精品一区www在线观看| 亚洲国产色片| 色综合站精品国产| 人妻系列 视频| 日本免费在线观看一区| 亚洲av一区综合| 午夜精品一区二区三区免费看| 久久亚洲精品不卡| 欧美精品一区二区大全| 国产精品女同一区二区软件| 久久精品熟女亚洲av麻豆精品 | 国产精品综合久久久久久久免费| 一级毛片久久久久久久久女| 亚洲在线自拍视频| 欧美日韩精品成人综合77777| 久久99热这里只有精品18| 精品熟女少妇av免费看| 黄片无遮挡物在线观看| 69人妻影院| 日日干狠狠操夜夜爽| 麻豆av噜噜一区二区三区| 免费av观看视频| 汤姆久久久久久久影院中文字幕 | 中文亚洲av片在线观看爽| 国产乱人视频| 亚洲国产精品专区欧美| 国产一区亚洲一区在线观看| av又黄又爽大尺度在线免费看 | 国产av不卡久久| av卡一久久| 最近的中文字幕免费完整| 精品久久久久久成人av| 欧美日韩国产亚洲二区| 国产精品一区二区性色av| 草草在线视频免费看| 18禁在线播放成人免费| 亚洲国产色片| 国产大屁股一区二区在线视频| 尾随美女入室| 一区二区三区高清视频在线| 美女国产视频在线观看| 狂野欧美激情性xxxx在线观看| 亚洲精品国产av成人精品| 国产精品久久久久久久久免| 久久久久久久久中文| 人妻夜夜爽99麻豆av| 日韩欧美三级三区| 99热6这里只有精品| 精品久久久噜噜| 最近中文字幕2019免费版| 91av网一区二区| 有码 亚洲区| 久久久色成人| 国产色爽女视频免费观看| 亚洲av二区三区四区| 亚洲成人中文字幕在线播放| 淫秽高清视频在线观看| 国产亚洲av嫩草精品影院| 国产午夜精品一二区理论片| 国产爱豆传媒在线观看| 五月玫瑰六月丁香| 99久国产av精品国产电影| 午夜久久久久精精品| 国语自产精品视频在线第100页| 深夜a级毛片| 国产精品久久久久久精品电影| 中文字幕熟女人妻在线| 亚洲五月天丁香| 熟女人妻精品中文字幕| 大又大粗又爽又黄少妇毛片口| 日本黄色视频三级网站网址| 免费av毛片视频| 九九在线视频观看精品| 欧美性猛交╳xxx乱大交人| 亚洲在线自拍视频| 少妇丰满av| 国产精品国产三级国产专区5o | 男女视频在线观看网站免费| 国产成人午夜福利电影在线观看| 久久久久久久久大av| 亚洲av不卡在线观看| 日韩精品青青久久久久久| 最近手机中文字幕大全| 国产精品乱码一区二三区的特点| 精品久久国产蜜桃| 亚洲精品一区蜜桃| 久久国产乱子免费精品| 国产美女午夜福利| 高清午夜精品一区二区三区| 岛国在线免费视频观看| 国产成人午夜福利电影在线观看| 久久人人爽人人爽人人片va| 在线播放国产精品三级| 国产 一区精品| 非洲黑人性xxxx精品又粗又长| 免费观看在线日韩| 毛片一级片免费看久久久久| h日本视频在线播放| av卡一久久| 国产免费男女视频| 国产精品爽爽va在线观看网站| 成人美女网站在线观看视频| 久久婷婷人人爽人人干人人爱| 久久精品国产99精品国产亚洲性色| 国内少妇人妻偷人精品xxx网站| 亚洲av免费高清在线观看| .国产精品久久| 日韩av在线大香蕉| 熟妇人妻久久中文字幕3abv| 亚洲精品国产成人久久av| 免费观看精品视频网站| 我的女老师完整版在线观看| av线在线观看网站| 中文欧美无线码| 最近的中文字幕免费完整| ponron亚洲| 精品欧美国产一区二区三| 精品国产一区二区三区久久久樱花 | 能在线免费观看的黄片| 国产一级毛片七仙女欲春2| 97热精品久久久久久| 国产精品.久久久| 搡老妇女老女人老熟妇| 国产精品不卡视频一区二区| 色哟哟·www| 免费观看精品视频网站| 国产亚洲91精品色在线| 日本一本二区三区精品| 美女国产视频在线观看| 日韩欧美国产在线观看| 国产成人免费观看mmmm| 久久精品国产亚洲av天美| eeuss影院久久| 国产av一区在线观看免费| 久久这里只有精品中国| 中文字幕亚洲精品专区| 少妇被粗大猛烈的视频| 久久草成人影院| 欧美日韩在线观看h| 国产亚洲一区二区精品| 国产极品天堂在线| 少妇熟女欧美另类| 日日啪夜夜撸| 少妇熟女欧美另类| 2021天堂中文幕一二区在线观| 久久久久性生活片| 国产极品天堂在线| 美女高潮的动态| 国产亚洲5aaaaa淫片| 少妇人妻精品综合一区二区| 永久免费av网站大全| 久久久国产成人免费| 我的老师免费观看完整版| 欧美色视频一区免费| 在线观看av片永久免费下载|