• <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ù)
    全区人妻精品视频| 久久人人爽人人爽人人片va | 国产精品不卡视频一区二区 | 日韩有码中文字幕| 波多野结衣高清作品| 白带黄色成豆腐渣| 一卡2卡三卡四卡精品乱码亚洲| 嫩草影院入口| 国产精品女同一区二区软件 | 国产成+人综合+亚洲专区| 国产精品伦人一区二区| 久久午夜福利片| 亚洲第一区二区三区不卡| 99热这里只有是精品在线观看 | 高潮久久久久久久久久久不卡| 国产大屁股一区二区在线视频| 国产亚洲av嫩草精品影院| 久久久久久国产a免费观看| 国产午夜精品久久久久久一区二区三区 | 亚洲aⅴ乱码一区二区在线播放| 国产精品,欧美在线| 成人特级黄色片久久久久久久| 午夜福利在线观看吧| 国产日本99.免费观看| 国产91精品成人一区二区三区| 人妻久久中文字幕网| 黄色女人牲交| 国产成+人综合+亚洲专区| av在线老鸭窝| 亚洲欧美日韩高清专用| 99久久成人亚洲精品观看| 国产一区二区在线av高清观看| 黄色日韩在线| 欧美xxxx性猛交bbbb| 变态另类丝袜制服| 国产极品精品免费视频能看的| 欧美最黄视频在线播放免费| 搡女人真爽免费视频火全软件 | 日韩欧美 国产精品| 国产真实乱freesex| 免费观看人在逋| 少妇丰满av| 午夜视频国产福利| 高清毛片免费观看视频网站| 亚洲美女视频黄频| 午夜日韩欧美国产| 简卡轻食公司| 俄罗斯特黄特色一大片| 成人特级黄色片久久久久久久| avwww免费| 3wmmmm亚洲av在线观看| 特大巨黑吊av在线直播| 成人性生交大片免费视频hd| 国产精品女同一区二区软件 | 欧美最新免费一区二区三区 | 国产一级毛片七仙女欲春2| 成人永久免费在线观看视频| 男女那种视频在线观看| 国产精品亚洲av一区麻豆| 国产亚洲av嫩草精品影院| 午夜亚洲福利在线播放| 欧美在线黄色| 欧美三级亚洲精品| 国产免费男女视频| 国产伦精品一区二区三区四那| 亚洲不卡免费看| 真人一进一出gif抽搐免费| 成人特级av手机在线观看| 我的女老师完整版在线观看| 亚洲无线观看免费| 真人一进一出gif抽搐免费| 国产黄色小视频在线观看| 一进一出好大好爽视频| 看黄色毛片网站| 亚洲精品影视一区二区三区av| 12—13女人毛片做爰片一| 免费无遮挡裸体视频| 欧美另类亚洲清纯唯美| 日韩av在线大香蕉| 中文资源天堂在线| 69人妻影院| 18禁黄网站禁片免费观看直播| 搞女人的毛片| 亚洲第一电影网av| 一区二区三区免费毛片| 欧美性猛交黑人性爽| 免费看a级黄色片| 麻豆成人av在线观看| 国产精品亚洲av一区麻豆| 亚洲五月婷婷丁香| 久久人人爽人人爽人人片va | 亚洲无线观看免费| 国产亚洲欧美98| 淫秽高清视频在线观看| 深夜a级毛片| 91午夜精品亚洲一区二区三区 | aaaaa片日本免费| 最近在线观看免费完整版| 九九热线精品视视频播放| 久久精品久久久久久噜噜老黄 | 精品午夜福利在线看| 日韩精品青青久久久久久| 国产精品98久久久久久宅男小说| 国产色爽女视频免费观看| 色av中文字幕| 在现免费观看毛片| 51午夜福利影视在线观看| 国产精品免费一区二区三区在线| 日韩精品中文字幕看吧| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 观看美女的网站| 在线观看免费视频日本深夜| 亚洲av中文字字幕乱码综合| 窝窝影院91人妻| 国产aⅴ精品一区二区三区波| 最近在线观看免费完整版| 午夜亚洲福利在线播放| 深夜精品福利| 精品久久久久久久久亚洲 | 久久人妻av系列| 99热6这里只有精品| 一区二区三区免费毛片| 国产伦精品一区二区三区视频9| 亚洲中文字幕日韩| 欧美最新免费一区二区三区 | 18+在线观看网站| 国产精品永久免费网站| 国产美女午夜福利| 嫁个100分男人电影在线观看| 久久久成人免费电影| 内射极品少妇av片p| 麻豆av噜噜一区二区三区| 国产淫片久久久久久久久 | 色吧在线观看| 99国产精品一区二区三区| 欧美高清性xxxxhd video| 天美传媒精品一区二区| 日本精品一区二区三区蜜桃| 亚洲精品在线美女| 国产在线精品亚洲第一网站| 51午夜福利影视在线观看| 久久精品影院6| 热99在线观看视频| 国产一区二区三区视频了| 欧美日本视频| www.熟女人妻精品国产| 又爽又黄无遮挡网站| 久久亚洲真实| 丰满人妻一区二区三区视频av| 国产高潮美女av| 亚洲av五月六月丁香网| 老熟妇乱子伦视频在线观看| 亚洲国产精品合色在线| 免费在线观看日本一区| 日本黄色片子视频| 熟女人妻精品中文字幕| 俄罗斯特黄特色一大片| 亚洲精品日韩av片在线观看| 麻豆久久精品国产亚洲av| 欧美不卡视频在线免费观看| 黄色日韩在线| 三级国产精品欧美在线观看| 日本黄大片高清| 国内精品一区二区在线观看| 欧美乱色亚洲激情| 少妇的逼水好多| 看黄色毛片网站| 可以在线观看毛片的网站| 日韩欧美在线二视频| 精品人妻熟女av久视频| 国产一区二区亚洲精品在线观看| 欧美午夜高清在线| 精品国产三级普通话版| 内射极品少妇av片p| 天堂av国产一区二区熟女人妻| 亚洲人成网站在线播| 久久久久性生活片| 美女cb高潮喷水在线观看| 日韩av在线大香蕉| 99在线人妻在线中文字幕| 男女床上黄色一级片免费看| 91麻豆精品激情在线观看国产| av欧美777| 亚洲专区中文字幕在线| 在线观看免费视频日本深夜| 亚洲人成网站在线播放欧美日韩| 亚洲电影在线观看av| 变态另类成人亚洲欧美熟女| 在线观看免费视频日本深夜| 久久久久国产精品人妻aⅴ院| 亚洲va日本ⅴa欧美va伊人久久| 欧美+日韩+精品| 真人做人爱边吃奶动态| 在线观看舔阴道视频| 特级一级黄色大片| 动漫黄色视频在线观看| 一本综合久久免费| 中文字幕av在线有码专区| 久久人妻av系列| 级片在线观看| 精品午夜福利视频在线观看一区| 两人在一起打扑克的视频| 国产精品嫩草影院av在线观看 | 欧美另类亚洲清纯唯美| 亚洲欧美激情综合另类| 成人午夜高清在线视频| 国产一区二区三区在线臀色熟女| 久久精品综合一区二区三区| 国内久久婷婷六月综合欲色啪| 99久国产av精品| 人妻久久中文字幕网| 精品日产1卡2卡| 性欧美人与动物交配| 国产亚洲欧美98| 麻豆av噜噜一区二区三区| 夜夜夜夜夜久久久久| 美女 人体艺术 gogo| 亚洲av第一区精品v没综合| 极品教师在线免费播放| 亚洲av熟女| 亚洲av熟女| 99国产精品一区二区蜜桃av| 一夜夜www| 俄罗斯特黄特色一大片| 亚洲内射少妇av| 久久亚洲精品不卡| 我的老师免费观看完整版| 精品不卡国产一区二区三区| 三级男女做爰猛烈吃奶摸视频| 激情在线观看视频在线高清| 99热这里只有是精品50| 欧美又色又爽又黄视频| 日韩欧美 国产精品| 日韩欧美国产一区二区入口| 婷婷精品国产亚洲av在线| 一夜夜www| 久久久久久国产a免费观看| 超碰av人人做人人爽久久| 亚洲一区二区三区色噜噜| 欧美乱色亚洲激情| 久久久久免费精品人妻一区二区| 非洲黑人性xxxx精品又粗又长| 国产色婷婷99| 久久久精品大字幕| 国产精品人妻久久久久久| 国产午夜精品论理片| 免费电影在线观看免费观看| 国产高清激情床上av| 精华霜和精华液先用哪个| 欧美bdsm另类| 久久精品夜夜夜夜夜久久蜜豆| 色噜噜av男人的天堂激情| 搡女人真爽免费视频火全软件 | 精品久久久久久久久久免费视频| 国产大屁股一区二区在线视频| 精品一区二区三区av网在线观看| 免费人成在线观看视频色| 赤兔流量卡办理| 国产午夜精品论理片| 搡老熟女国产l中国老女人| 91九色精品人成在线观看| 9191精品国产免费久久| 精品久久久久久久久久免费视频| 高清在线国产一区| 伦理电影大哥的女人| 日韩国内少妇激情av| 又黄又爽又免费观看的视频| 免费av不卡在线播放| ponron亚洲| 草草在线视频免费看| 久久婷婷人人爽人人干人人爱| 高清日韩中文字幕在线| 国产一区二区亚洲精品在线观看| a级毛片免费高清观看在线播放| av视频在线观看入口| 欧美成人一区二区免费高清观看| 国模一区二区三区四区视频| 青草久久国产| 蜜桃亚洲精品一区二区三区| 久久人妻av系列| 99久久成人亚洲精品观看| 亚洲人成网站在线播| 丰满人妻熟妇乱又伦精品不卡| 乱码一卡2卡4卡精品| 高潮久久久久久久久久久不卡| 天堂av国产一区二区熟女人妻| 国产精品不卡视频一区二区 | 欧美+亚洲+日韩+国产| 又粗又爽又猛毛片免费看| 最好的美女福利视频网| av国产免费在线观看| 男女做爰动态图高潮gif福利片| 深夜精品福利| 精品久久国产蜜桃| 欧美性感艳星| 最近最新免费中文字幕在线| 亚洲成av人片免费观看| 免费在线观看亚洲国产| 精品午夜福利视频在线观看一区| 欧美3d第一页| 中文字幕免费在线视频6| 精品无人区乱码1区二区| 一级黄片播放器| 一本一本综合久久| av女优亚洲男人天堂| 成人性生交大片免费视频hd| xxxwww97欧美| 别揉我奶头 嗯啊视频| 99热这里只有精品一区| 亚洲国产精品久久男人天堂| 亚洲av成人精品一区久久| 午夜精品一区二区三区免费看| 好看av亚洲va欧美ⅴa在| 成年版毛片免费区| 亚洲av二区三区四区| 身体一侧抽搐| 91午夜精品亚洲一区二区三区 | 亚洲av不卡在线观看| 国产大屁股一区二区在线视频| 国产精品电影一区二区三区| 成人高潮视频无遮挡免费网站| 1000部很黄的大片| 无遮挡黄片免费观看| 在线免费观看的www视频| 免费大片18禁| 日韩欧美精品免费久久 | 在线a可以看的网站| 精品一区二区三区人妻视频| 亚洲成人精品中文字幕电影| 舔av片在线| 国产中年淑女户外野战色| 亚洲成av人片在线播放无| 精品不卡国产一区二区三区| 精品久久久久久,| 国产成人福利小说| 久久久久性生活片| 啦啦啦韩国在线观看视频| 天天一区二区日本电影三级| 999久久久精品免费观看国产| 搡女人真爽免费视频火全软件 | 国产精品影院久久| 精品国产三级普通话版| 18禁在线播放成人免费| 国产欧美日韩一区二区精品| 亚洲av电影不卡..在线观看| 国产成+人综合+亚洲专区| 欧美日韩黄片免| 黄色女人牲交| 亚洲五月婷婷丁香| 国产精品亚洲av一区麻豆| 国产av麻豆久久久久久久| 亚洲精品一卡2卡三卡4卡5卡| 大型黄色视频在线免费观看| 亚洲国产精品999在线| 国产一区二区在线av高清观看| а√天堂www在线а√下载| 99久久精品国产亚洲精品| 18+在线观看网站| 国产精品98久久久久久宅男小说| 亚洲专区中文字幕在线| 天堂√8在线中文| 午夜精品在线福利| 久久国产精品人妻蜜桃| 日日摸夜夜添夜夜添av毛片 | 99在线视频只有这里精品首页| 亚洲人成网站在线播| 欧美日本视频| 国产一级毛片七仙女欲春2| 婷婷色综合大香蕉| 久久久久久大精品| 一夜夜www| 国产麻豆成人av免费视频| 欧美高清性xxxxhd video| 国产熟女xx| 亚洲中文字幕一区二区三区有码在线看| 国产精品影院久久| 精品乱码久久久久久99久播| 少妇裸体淫交视频免费看高清| 99精品久久久久人妻精品| 国产 一区 欧美 日韩| 国产 一区 欧美 日韩| 中文字幕av在线有码专区| 久久香蕉精品热| 亚洲片人在线观看| 午夜亚洲福利在线播放| bbb黄色大片| av黄色大香蕉| 国产探花极品一区二区| 麻豆成人av在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲av美国av| 成人特级av手机在线观看| 国产精品久久久久久人妻精品电影| 欧美黄色淫秽网站| 国产精品野战在线观看| 12—13女人毛片做爰片一| 黄色女人牲交| a级毛片a级免费在线| 日韩欧美精品免费久久 | 九色国产91popny在线| 看黄色毛片网站| 网址你懂的国产日韩在线| 日韩中文字幕欧美一区二区| 丰满的人妻完整版| 啪啪无遮挡十八禁网站| 99在线视频只有这里精品首页| 97超级碰碰碰精品色视频在线观看| 亚洲在线自拍视频| 亚洲经典国产精华液单 | 国产伦人伦偷精品视频| 一卡2卡三卡四卡精品乱码亚洲| 婷婷精品国产亚洲av| 他把我摸到了高潮在线观看| 夜夜爽天天搞| 午夜激情欧美在线| 在线播放无遮挡| 国产亚洲精品久久久久久毛片| 国产免费av片在线观看野外av| 一区二区三区免费毛片| 看片在线看免费视频| 中文字幕人成人乱码亚洲影| 精品欧美国产一区二区三| 18+在线观看网站| 免费看a级黄色片| 色播亚洲综合网| 观看免费一级毛片| 永久网站在线| 美女xxoo啪啪120秒动态图 | 国产白丝娇喘喷水9色精品| 日韩中字成人| 少妇裸体淫交视频免费看高清| 国产成人欧美在线观看| 国产免费男女视频| 欧美在线黄色| 一个人观看的视频www高清免费观看| 特大巨黑吊av在线直播| 小说图片视频综合网站| 国产主播在线观看一区二区| 久久6这里有精品| 99热这里只有精品一区| 99久久九九国产精品国产免费| 国产精品av视频在线免费观看| 亚洲精华国产精华精| 亚洲美女搞黄在线观看 | 成人国产综合亚洲| 搞女人的毛片| 少妇高潮的动态图| 久久这里只有精品中国| 色播亚洲综合网| 日本三级黄在线观看| 看十八女毛片水多多多| 国产蜜桃级精品一区二区三区| 最近视频中文字幕2019在线8| 亚洲成a人片在线一区二区| 国产精品人妻久久久久久| 亚洲精华国产精华精| 日本黄色片子视频| 两性午夜刺激爽爽歪歪视频在线观看| av国产免费在线观看| 91久久精品国产一区二区成人| 欧美日韩亚洲国产一区二区在线观看| 国产免费男女视频| 97人妻精品一区二区三区麻豆| 免费看美女性在线毛片视频| 亚洲成a人片在线一区二区| 亚洲一区二区三区不卡视频| 精品人妻1区二区| 一级毛片久久久久久久久女| 一级av片app| 老司机深夜福利视频在线观看| 日本五十路高清| av在线天堂中文字幕| 午夜福利成人在线免费观看| 亚洲电影在线观看av| 久久精品久久久久久噜噜老黄 | 精品一区二区三区视频在线| 在线播放国产精品三级| 国产精品一区二区三区四区免费观看 | 性色avwww在线观看| 色尼玛亚洲综合影院| h日本视频在线播放| 丰满人妻一区二区三区视频av| x7x7x7水蜜桃| 小说图片视频综合网站| 亚洲狠狠婷婷综合久久图片| 亚洲七黄色美女视频| 中文字幕久久专区| 看免费av毛片| 免费在线观看影片大全网站| 亚洲欧美日韩无卡精品| 欧美性猛交黑人性爽| 国产美女午夜福利| 综合色av麻豆| 最新在线观看一区二区三区| 欧美高清成人免费视频www| 日本黄大片高清| 亚洲,欧美精品.| 神马国产精品三级电影在线观看| 国产91精品成人一区二区三区| www.www免费av| 国产在视频线在精品| 看十八女毛片水多多多| 宅男免费午夜| 欧美色视频一区免费| 国内久久婷婷六月综合欲色啪| 欧美日韩综合久久久久久 | 舔av片在线| 亚洲人与动物交配视频| 看片在线看免费视频| 成人毛片a级毛片在线播放| 精品午夜福利在线看| 久久午夜亚洲精品久久| 99热这里只有精品一区| 国产精品乱码一区二三区的特点| 成人三级黄色视频| 毛片一级片免费看久久久久 | 日本与韩国留学比较| a级毛片免费高清观看在线播放| 国产视频内射| 9191精品国产免费久久| 免费看a级黄色片| 91久久精品国产一区二区成人| 欧美国产日韩亚洲一区| 女同久久另类99精品国产91| 免费电影在线观看免费观看| 1024手机看黄色片| 午夜视频国产福利| 午夜福利18| 夜夜爽天天搞| 成人一区二区视频在线观看| 97碰自拍视频| 天美传媒精品一区二区| 欧美国产日韩亚洲一区| 三级国产精品欧美在线观看| 久久久久久国产a免费观看| 亚洲人成伊人成综合网2020| 国产伦在线观看视频一区| 久久国产乱子免费精品| 啪啪无遮挡十八禁网站| 亚洲天堂国产精品一区在线| 中文亚洲av片在线观看爽| 亚洲精品粉嫩美女一区| 国产视频内射| 久久国产精品影院| 亚洲精品色激情综合| 很黄的视频免费| 亚洲av五月六月丁香网| 哪里可以看免费的av片| 99热这里只有是精品在线观看 | 久久这里只有精品中国| 18禁裸乳无遮挡免费网站照片| 久久人妻av系列| av天堂在线播放| 丝袜美腿在线中文| 免费人成在线观看视频色| 午夜福利18| 毛片女人毛片| 在线播放无遮挡| 国产一区二区在线观看日韩| 毛片女人毛片| 国产免费男女视频| 国产91精品成人一区二区三区| 亚洲欧美日韩高清专用| 国产午夜精品论理片| 久久精品国产亚洲av香蕉五月| 免费在线观看亚洲国产| 日本成人三级电影网站| 欧美日韩中文字幕国产精品一区二区三区| 久久精品国产99精品国产亚洲性色| 久久国产精品影院| 成年女人永久免费观看视频| 首页视频小说图片口味搜索| 757午夜福利合集在线观看| 精品乱码久久久久久99久播| x7x7x7水蜜桃| 国产精品影院久久| 国产单亲对白刺激| www日本黄色视频网| 国产精品伦人一区二区| 成年女人看的毛片在线观看| 90打野战视频偷拍视频| 国产久久久一区二区三区| 2021天堂中文幕一二区在线观| 久久午夜福利片| 波多野结衣巨乳人妻| 色综合亚洲欧美另类图片| 亚洲人成电影免费在线| 国产黄片美女视频| 18禁在线播放成人免费| 十八禁国产超污无遮挡网站| 18禁黄网站禁片免费观看直播| 欧美成狂野欧美在线观看| 亚洲国产精品合色在线| 国产一区二区激情短视频| 久久伊人香网站| 自拍偷自拍亚洲精品老妇| 一本一本综合久久| 网址你懂的国产日韩在线| 深爱激情五月婷婷| 乱码一卡2卡4卡精品| 此物有八面人人有两片| 18禁裸乳无遮挡免费网站照片| 久久精品国产亚洲av涩爱 | 久久香蕉精品热| 国模一区二区三区四区视频| 变态另类成人亚洲欧美熟女| 日韩高清综合在线| 午夜两性在线视频| 18美女黄网站色大片免费观看| 国产精品1区2区在线观看.| 99国产精品一区二区蜜桃av| 久久久久国产精品人妻aⅴ院| 日日干狠狠操夜夜爽| 欧美日韩瑟瑟在线播放| 女人十人毛片免费观看3o分钟|