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

    前緣附近微小板對(duì)S809 翼型氣動(dòng)性能的影響

    2021-03-25 01:27:08許勝焱黃典貴
    關(guān)鍵詞:吸力振幅氣動(dòng)

    許勝焱, 許 斌, 黃典貴

    (上海理工大學(xué) 能源與動(dòng)力工程學(xué)院,上海 200093)

    隨著全球化石燃料的大量消耗和環(huán)境問題的日益增加,核能、風(fēng)能和地?zé)崮艿雀鞣N新能源已成為各國(guó)關(guān)注的熱點(diǎn)[1]。風(fēng)能是一種低成本與可再生的清潔能源,風(fēng)能利用是通過風(fēng)力機(jī)將風(fēng)的動(dòng)能轉(zhuǎn)換為機(jī)械能、電能和熱能等。葉片氣動(dòng)性能的研究是風(fēng)力機(jī)設(shè)計(jì)的關(guān)鍵,葉片最基本的元素是翼型,翼型氣動(dòng)性能直接決定風(fēng)輪葉片性能,進(jìn)而影響風(fēng)力機(jī)風(fēng)能利用率。因此,翼型升阻比是決定風(fēng)力機(jī)風(fēng)能利用率的關(guān)鍵因素。但是,在大來(lái)流攻角下,翼型繞流表面曲率大,尾緣附近逆壓梯度大,翼型表面流體脫離吸力面,形成層流分離泡或發(fā)展成逐漸向尾緣脫落的分離渦,最終導(dǎo)致翼型氣動(dòng)性能降低,嚴(yán)重影響風(fēng)能利用率。翼型吸力面流場(chǎng)參數(shù)波動(dòng)大,葉片易發(fā)生結(jié)構(gòu)疲勞而導(dǎo)致壽命縮短或引發(fā)安全問題。因此,翼型繞流的分離控制一直受到關(guān)注。

    翼型在深失速后,吸力面上方存在大分離區(qū),分離點(diǎn)在翼型前緣附近位置,分離渦向尾緣脫落時(shí)逐漸變大,形成前緣失速(leading-edge stall)現(xiàn)象。對(duì)于翼型流動(dòng)分離引起失速的問題,目前主要采用流動(dòng)控制手段進(jìn)行緩解或抑制。流動(dòng)控制方法主要分為被動(dòng)控制和主動(dòng)控制。被動(dòng)控制的方法主要是通過附加機(jī)械結(jié)構(gòu),不需要能量的注入。例如,在邊界層加肋、開槽等。以Gurney 襟翼為代表的被動(dòng)型流動(dòng)控制技術(shù)在水平軸風(fēng)力機(jī)中得到廣泛應(yīng)用與發(fā)展[2-3]。主動(dòng)控制方法需要注入能量,通過控制能量輸入實(shí)現(xiàn)流動(dòng)控制。例如,合成射流(synthetic jet)是一種由激勵(lì)器交替吹吸周圍流體而產(chǎn)生的非連續(xù)零質(zhì)量射流。Maldonado等[4]通過在風(fēng)洞實(shí)驗(yàn)使用合成射流器,探討合成射流對(duì)風(fēng)力機(jī)葉片氣動(dòng)性能的影響,發(fā)現(xiàn)在合成射流器的作用下,葉片動(dòng)態(tài)失速振動(dòng)減少且翼型升力系數(shù)大大提高。

    目前在翼型上加肋、開槽和射流口都是直接開在翼型吸力面上,雖然尺寸很小,但還是會(huì)影響原始經(jīng)典翼型的幾何外形。Wang 等[5]對(duì)翼型前緣加小圓柱的流動(dòng)分離控制方法進(jìn)行了研究,結(jié)果發(fā)現(xiàn),該控制方法可以有效地運(yùn)用在風(fēng)力機(jī)上,提高風(fēng)力機(jī)的風(fēng)能利用率。徐帥等[6]將翼型前緣加裝微小圓柱的方法應(yīng)用于S809 翼型,結(jié)果表明,振動(dòng)小圓柱在合適的控制參數(shù)下能提高翼型的氣動(dòng)性能。Zhou 等[7]和Shi 等[8]對(duì)微小結(jié)構(gòu)控制流場(chǎng)進(jìn)行研究,結(jié)果表明,微小結(jié)構(gòu)可以有效控制流場(chǎng)流動(dòng)。這些研究都是在翼型前緣附近設(shè)置離體結(jié)構(gòu),從而有效控制流場(chǎng)的流動(dòng)分離,同時(shí)減小對(duì)原始翼型幾何形狀的影響,改造較為簡(jiǎn)單,經(jīng)濟(jì)實(shí)用,通用性強(qiáng)。為此,本文采用一種在葉片前緣附近設(shè)置微小板的流動(dòng)控制方法,相對(duì)于文獻(xiàn)[5-6, 8],加設(shè)微小板如同擾流板,利用微小板產(chǎn)生誘導(dǎo)渦改善翼型的流場(chǎng)特征。相對(duì)于Zhou 等[7],振動(dòng)的微小板對(duì)翼型流場(chǎng)的擾動(dòng)更劇烈,本文探討其抑制翼型表面流動(dòng)分離的可能性。同時(shí)分析微小板的位置、尺寸和振動(dòng)參數(shù)對(duì)S809 翼型氣動(dòng)性能的影響。

    1 模型驗(yàn)證

    參考Shi 等[8]的計(jì)算模型,在FLUNENT 軟件中選用κ-ωShear Stress Transfer 湍流模型,κ為湍動(dòng)能,ε為比耗散率。采用二階精度計(jì)算,速度入口,來(lái)流馬赫數(shù)為0.084 193,壓力出口邊界條件,翼型及微小板邊界設(shè)置無(wú)滑移壁面,Y+為壁面率,Y+約為1,網(wǎng)格量約為18 萬(wàn),翼型第一層網(wǎng)格高度為2.33×10-5,翼型弦長(zhǎng)c=0.5 m。圖1 為葉片與微小板的位置關(guān)系示意圖。圖2 為計(jì)算域拓?fù)浣Y(jié)構(gòu)圖,計(jì)算域A 區(qū)為結(jié)構(gòu)化網(wǎng)格,最大區(qū)域半徑為50 倍翼型弦長(zhǎng),B 區(qū)為非結(jié)構(gòu)網(wǎng)格,最大區(qū)域半徑為2.4c,微小板附近區(qū)域C 及葉片附近區(qū)域D 均為結(jié)構(gòu)化網(wǎng)格。圖3 為結(jié)構(gòu)網(wǎng)格和非結(jié)構(gòu)網(wǎng)格交界處示意圖。

    圖1 葉片及微小板位置關(guān)系示意圖Fig.1 Coordinate diagram of the relative position between the blade and tiny sheet

    為了驗(yàn)證本文數(shù)值模擬結(jié)果的可靠性,在雷諾數(shù)Re=1×106下模擬單個(gè)S809 翼型的升阻力系數(shù)隨攻角的變化情況,并與Hand 等[9]的實(shí)驗(yàn)結(jié)果進(jìn)行比較,結(jié)果如圖4 所示。CL,CD分別為升力系數(shù)和阻力系數(shù)。本文模擬的升阻力系數(shù)及變化趨勢(shì)與實(shí)驗(yàn)值較為吻合,所采用的模擬方法和模型滿足本文計(jì)算要求。

    圖4 本文的模擬結(jié)果與實(shí)驗(yàn)的對(duì)比Fig.4 Comparison between the results of simulation and experiments

    2 計(jì)算結(jié)果及分析

    2.1 靜止微小板對(duì)流場(chǎng)影響分析

    微小板板厚為翼型弦長(zhǎng)的千分之一,研究靜止微小板長(zhǎng)度L和擺放位置(X,Y)對(duì)S809 翼型氣動(dòng)性能的影響。為了方便表達(dá),L=0.005c表示微小板板長(zhǎng)L為翼型弦長(zhǎng)的0.5%;(X,Y)=(-2,4)表示微小板幾何中心位置距翼型前緣點(diǎn)的水平距離為0.02c,垂直距離為0.04c。

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

    當(dāng)微小板L=0.02c在(0,8)時(shí),設(shè)置4 套網(wǎng)格數(shù)量進(jìn)行無(wú)關(guān)性檢驗(yàn),計(jì)算設(shè)置與單個(gè)翼型數(shù)值模擬設(shè)置一樣,仿真模擬達(dá)到穩(wěn)態(tài)后,得出的結(jié)果如表1 所示。

    由表1 可知,網(wǎng)格數(shù)量對(duì)仿真模擬結(jié)果具有一定的影響,當(dāng)網(wǎng)格數(shù)量超過18 萬(wàn)時(shí),模擬結(jié)果變化差異不大,即網(wǎng)格數(shù)量增大到一定程度后計(jì)算結(jié)果趨于穩(wěn)定。綜合考慮計(jì)算時(shí)間與仿真誤差,可以近似認(rèn)為第三套35 萬(wàn)的網(wǎng)格數(shù)量可滿足網(wǎng)格數(shù)量的無(wú)關(guān)性要求,故選擇第三套網(wǎng)格數(shù)量作為模擬計(jì)算的網(wǎng)格量。

    表1 網(wǎng)格無(wú)關(guān)性驗(yàn)證結(jié)果Tab.1 Results of the grid independence verification

    2.1.2 結(jié)果分析

    圖5 是在22°攻角下,4 種不同板長(zhǎng)的靜止微小板在相同位置下翼型升阻比的變化規(guī)律。圖中橫坐標(biāo)為微小板中心與前緣點(diǎn)的垂直距離Y,圖例0.005c等表示板長(zhǎng)L。

    圖5(a)中靜止微小板的中心位置橫坐標(biāo)X=-4,當(dāng)板長(zhǎng)為0.005c或0.01c時(shí),微小板尺寸較小,距離翼型前緣位置較遠(yuǎn),靜止板對(duì)翼型吸力面的流體流動(dòng)影響小,升阻比變化不大;增大板長(zhǎng)為0.02c或0.03c時(shí),微小板對(duì)流場(chǎng)控制效果明顯,翼型升阻比隨Y的增加均先增大后減小,Y=6~11 范圍內(nèi)翼型的升阻比較大,Y=9 時(shí)獲得最佳升阻比,比原始翼型提升了1 倍;Y值繼續(xù)增大,即板在Y方向逐漸遠(yuǎn)離翼型,微小板對(duì)吸力面流場(chǎng)流動(dòng)控制能力越來(lái)越小,翼型升阻比迅速下降。圖5(b)中微小板位置在X方向增大至X=-2,板長(zhǎng)為0.005c或0.01c的微小板對(duì)翼型吸力面流體流動(dòng)的影響增大,板長(zhǎng)為0.005c的翼型升阻比隨Y的增大而減小,板長(zhǎng)為0.01c的翼型升阻比隨Y的增大而先增大后減小。增大板長(zhǎng)為0.02c或0.03c時(shí),翼型升阻比變化與微小板在X=-4 時(shí)具有相似的規(guī)律,隨Y的增加先增大后減小。

    圖5 靜止微小板不同位置(X,Y)下,翼型升阻比的變化規(guī)律Fig.5 Variation of the lift-drag ratio of airfoils with a stationary tiny sheet set at different positions (X, Y)

    由圖5(c)和5(d)可知,微小板位置在X方向增大至X=0, 2 處,板長(zhǎng)為0.005c或0.01c的微小板對(duì)流場(chǎng)控制效果更明顯,翼型的升阻比提高;當(dāng)微小板在Y=3~7 時(shí),板長(zhǎng)為0.01c的翼型升阻比大于加設(shè)其他板長(zhǎng)的翼型升阻比;增大板長(zhǎng)為0.02c或0.03c時(shí),翼型升阻比隨Y的增加先增大后減小。圖5(e)中微小板位置在X方向增大至X=4 處,板長(zhǎng)為0.005c的翼型升阻比隨Y的增大而減小,板長(zhǎng)為0.02c或0.03c時(shí),翼型升阻比隨Y的增加先增大后減小。板長(zhǎng)為0.01c的微小板在(4,6)附近時(shí)翼型升阻比小于在其他位置的翼型升阻比。

    圖6 是4 種尺寸微小板在不同安裝位置時(shí)的翼型升阻比變化云圖,紅色區(qū)域?yàn)樯璞认鄬?duì)較高區(qū),即有效控制區(qū),藍(lán)色區(qū)域?yàn)樯璞认鄬?duì)較低區(qū)。隨著板長(zhǎng)的增大,翼型升阻比較高區(qū)域面積逐漸增大且高效區(qū)逐漸遠(yuǎn)離翼型。微小板尺寸較小時(shí),如圖6(a)和6(b),板安裝位置只在翼型附近升阻比才有提高。增大微小板尺寸,如圖6(c)和6(d),在翼型附近控制效果差。此外,4 種尺寸微小板所對(duì)應(yīng)的最佳安裝位置如表2 所示,隨著板長(zhǎng)的增加,翼型升阻比提升比率越大,微小板安裝位置距離翼型也越遠(yuǎn)。

    表2 最佳靜止微小板參數(shù)表Tab.2 Data sheets of the best stationary tiny sheet

    通過分析翼型附近流場(chǎng)的渦量流線圖,可進(jìn)一步探討翼型的流動(dòng)分離機(jī)理。為了更加直觀地顯示微小板板長(zhǎng)對(duì)翼型氣動(dòng)性能的影響,圖7 和圖8 分別給出了22°攻角下,原始翼型和在(-4,9)處添加不同板長(zhǎng)微小板后的翼型渦量流線圖。圖7是原始S809 翼型的渦量流線圖,從翼型前緣位置起出現(xiàn)了一個(gè)細(xì)長(zhǎng)上翹的負(fù)渦區(qū),翼型整個(gè)吸力面被巨大的分離泡所覆蓋。當(dāng)在(-4,9)處添加L=0.005c的微小板后,如圖8(a)所示,與翼型前緣位置距離遠(yuǎn),尺寸較小時(shí),微小板對(duì)翼型周圍流場(chǎng)的影響較小,流動(dòng)還是在前緣附近發(fā)生了分離;圖8(a)中的流線顯示,其產(chǎn)生的分離泡幾乎將翼型整個(gè)吸力面覆蓋;由圖8(a)中的渦量云圖看出,上翹的負(fù)渦區(qū)依然存在,這是誘發(fā)流動(dòng)分離的關(guān)鍵因素。此時(shí),靜止微小板板長(zhǎng)過小,抑制流動(dòng)分離的能力有限,因而升阻比較低。當(dāng)增大微小板板長(zhǎng)為0.01c時(shí),如圖8(b)所示,吸力面上分離泡減小,負(fù)渦區(qū)受到微小板產(chǎn)生的渦的影響,不再是細(xì)長(zhǎng)上翹的。當(dāng)繼續(xù)增大板長(zhǎng)為0.02c時(shí),如圖8(c)所示,從微小板脫落的正負(fù)渦明顯影響了翼型的負(fù)渦區(qū)不再往上翹,分離泡顯著減小,該微小板抑制流動(dòng)分離效果顯著。當(dāng)繼續(xù)增大板長(zhǎng)為0.03c時(shí),如圖8(d)所示,微小板對(duì)翼型分離區(qū)有明顯干涉作用,微小板生成的旋渦與翼型的旋渦之間存在較強(qiáng)的相互作用,流動(dòng)分離點(diǎn)往后移,吸力面上分離泡減小,翼型脫落尾渦減小,翼型流動(dòng)分離得到控制,控制效果優(yōu)于圖8(c)。

    圖7 原始S809 翼型渦量流線圖Fig.7 Vorticity flow diagram of the original S809 airfoil

    通過分析在不同位置下同一微小板對(duì)翼型流場(chǎng)流動(dòng)的控制,可以更加直觀地顯示微小板的位置對(duì)翼型氣動(dòng)性能的影響。由圖8(d)~8(f)可知,在(-4,9)處的微小板對(duì)翼型流場(chǎng)的流動(dòng)分離控制較好,使得翼型流動(dòng)分離點(diǎn)往后移,吸力面上分離泡減?。划?dāng)微小板在翼型表面較近時(shí),在(0,3)處的微小板與翼型上表面的空隙小,對(duì)來(lái)流氣流產(chǎn)生阻滯作用,翼型吸力面附近的氣流總體動(dòng)能減小,抵抗逆壓梯度的能力減弱,從而在翼型吸力面尾緣出現(xiàn)大范圍的流動(dòng)分離,微小板產(chǎn)生的渦直接與翼型表面接觸,翼型上表面流體流動(dòng)紊亂無(wú)序,分離泡大且斜向上,此時(shí)微小板對(duì)翼型的控制效果呈負(fù)作用;當(dāng)微小板離翼型表面較遠(yuǎn)時(shí),在(-4,16)處的微小板與翼型上表面的空隙足夠大,此垂直位置添加微小板已沒有明顯的流動(dòng)控制效果,從翼型前緣位置開始流動(dòng)分離,翼型整個(gè)吸力面被巨大的分離泡所覆蓋。

    圖8 翼型渦量流線圖(1)Fig.8 Streamline-vorticity diagram of the airfoil (1)

    在(4,6)附近是翼型前緣生成負(fù)渦向后發(fā)展的位置。分離渦離開翼面,該微小板就像導(dǎo)流板一樣,將分離渦導(dǎo)離翼型,微小板和翼型形成的2 個(gè)分離渦沒有抵消,而是合成為一個(gè)較強(qiáng)的分離渦,該較強(qiáng)分離渦在來(lái)流的作用下向尾緣運(yùn)動(dòng),分離泡逐漸變大,翼型氣動(dòng)性能變差。這種現(xiàn)象在Pauley 等[10]的平板分離泡數(shù)值模擬中同樣出現(xiàn)。在(4,5)和(4,7)處的微小板對(duì)翼型有著明顯的流動(dòng)控制效果,分離泡減小。

    翼型的氣動(dòng)性能對(duì)微小板的位置和長(zhǎng)度十分敏感,適當(dāng)尺寸的微小板安裝在合適位置才能達(dá)到流動(dòng)控制效果。當(dāng)Y值較小,即微小板在離翼型表面較近時(shí),微小板作用于翼型表面邊界層的流動(dòng),小尺寸微小板形成的渦,與翼型產(chǎn)生的渦相互影響,對(duì)吸力面的流體流動(dòng)有著積極作用;大尺寸微小板與翼型表面靠近時(shí),沒有改善流動(dòng)分離,反而對(duì)吸力面流場(chǎng)流動(dòng)有堵塞作用,導(dǎo)致翼型氣動(dòng)性能變差。如圖8(a)~8(d)所示,當(dāng)Y值增大至9 時(shí),大尺寸微小板對(duì)翼型流動(dòng)分離控制效果明顯,翼型流動(dòng)分離點(diǎn)明顯后移,分離泡變??;Y值繼續(xù)增加,如圖8(f)所示,微小板遠(yuǎn)離翼型,流動(dòng)控制能力小,翼型流場(chǎng)與圖7 原始S809 翼型的流場(chǎng)相似。

    圖9 翼型渦量流線圖(2)Fig.9 Streamline-vorticity diagram of the airfoil(2)

    2.2 翼型前緣點(diǎn)附近設(shè)置振動(dòng)微小板對(duì)流場(chǎng)影響分析

    在來(lái)流攻角22°的工況下,賦予放置于(0,6)處的微小板一定的振幅和頻率的振動(dòng),進(jìn)行數(shù)值模擬分析,以此探討微小板振動(dòng)效應(yīng)對(duì)翼型氣動(dòng)性能的影響。

    定義f為微小板無(wú)量綱振動(dòng)頻率,其數(shù)值根據(jù)式(1)確定。

    式中:fs為實(shí)際振動(dòng)頻率;L為微小板板長(zhǎng);U∞為來(lái)流風(fēng)速。

    設(shè)定微小板受迫振動(dòng)的運(yùn)動(dòng)方程為

    式中,A為Y方向的幅值,即定義微小板的振動(dòng)形式為在Y軸上平動(dòng)。

    2.2.1 無(wú)關(guān)性驗(yàn)證

    針對(duì)微小板L=0.01c在(0,6)處以振幅A為0.025c振動(dòng),按時(shí)間步長(zhǎng)為0.000 1 s,設(shè)置a,b,c,d 這4 套網(wǎng)格數(shù)量進(jìn)行無(wú)關(guān)性檢驗(yàn),仿真模擬達(dá)到穩(wěn)態(tài)后,得出的結(jié)果如圖10(a)所示。圖例0 是振幅A為零,圖例a,b,c,d 分別是網(wǎng)格數(shù)量為148 233,320 389,608 632,100 908 個(gè)。

    由圖10(a)可知,網(wǎng)格數(shù)量對(duì)結(jié)果具有一定的影響,當(dāng)頻率f< 1.5 時(shí),模擬結(jié)果變化差異不大;當(dāng)頻率增大時(shí),只有網(wǎng)格數(shù)量增大到一定程度后,計(jì)算結(jié)果趨于穩(wěn)定。綜合考慮計(jì)算時(shí)間與仿真誤差,可以近似認(rèn)為網(wǎng)格數(shù)量c 可滿足網(wǎng)格數(shù)量的無(wú)關(guān)性要求。

    以網(wǎng)格數(shù)量c 為網(wǎng)格模型,針對(duì)微小板L=0.01c在(0,6)處以振幅A為0.025c振動(dòng),按時(shí)間步長(zhǎng)0.005,0.001,0.000 1,0.000 05 s,對(duì)所選網(wǎng)格數(shù)量進(jìn)行時(shí)間獨(dú)立性檢驗(yàn)。在模擬進(jìn)行到穩(wěn)態(tài)后,得出的計(jì)算結(jié)果如圖10(b)所示。在網(wǎng)格數(shù)量達(dá)到一定程度后,時(shí)間步長(zhǎng)數(shù)值的大小對(duì)微小板是小頻率振動(dòng)的模擬結(jié)果差異不大,但振動(dòng)頻率較大時(shí),模擬結(jié)果的差異明顯;當(dāng)時(shí)間步長(zhǎng)達(dá)到0.001 s和0.000 1 s 時(shí),模擬結(jié)果較為接近,綜合考慮時(shí)間步長(zhǎng)對(duì)誤差的影響,同時(shí)結(jié)合時(shí)間步長(zhǎng)對(duì)仿真計(jì)算時(shí)間的影響,選取0.000 1 s 的時(shí)間步長(zhǎng)更為適宜。

    2.2.2 結(jié)果分析

    改變微小板的振動(dòng)頻率和振幅,翼型升阻比變化規(guī)律如圖11 所示。圖11(a)中,前緣加了振動(dòng)微小板的翼型最大升阻比提高80%以上,原因是微小板尺寸小,注入的能量在微小板上的損耗小,能量能較多較快地傳遞到流場(chǎng)中,翼型升阻比顯著提高[11]。圖11(b)中,L=0.005c微小板在小振幅小頻率時(shí)類似于靜止板,故在A=0.003c,f< 1時(shí)的翼型與加設(shè)靜止板的翼型的升阻比相差不大;當(dāng)f> 1,該振幅下的翼型升阻比隨著頻率的增大而增大。

    圖10 無(wú)關(guān)性驗(yàn)證結(jié)果Fig.10 Results of the independence verification

    圖11 翼型升阻比隨微小板的振動(dòng)頻率及振幅變化的規(guī)律Fig.11 Lift-drag ratio of the airfoil varied with the vibration frequency and amplitude of the tiny sheet

    圖11 (c)中,當(dāng)L=0.01c,f=0.5,振幅A=0.003c,0.005c或0.01c時(shí),翼型的升阻比與加設(shè)靜止微小板的翼型升阻比差不多,原因是小振幅小頻率的振動(dòng)板類似于靜止板,微小板振動(dòng)與靜止對(duì)翼型的氣動(dòng)性能的影響相差不大;當(dāng)振幅增大為0.02c或0.03c時(shí),微小板顯著提高翼型的升阻比。當(dāng)L=0.01c,A=0.01c,0.015c,0.02c時(shí),翼型升阻比變化規(guī)律相似,小頻率振動(dòng)對(duì)翼型升阻比影響不大,隨著頻率的增加,較大振幅的振動(dòng)微小板控制翼型升阻比先迅速增大;頻率繼續(xù)增加,翼型的升阻比有所下降,較大振幅的振動(dòng)微小板控制翼型升阻比先減小,但仍對(duì)控制翼型氣動(dòng)性能有著積極作用;在大頻率下,翼型升阻比有所下降,原因是微小板大頻率的振動(dòng),對(duì)翼型吸力面的來(lái)流擾動(dòng)過于激烈,導(dǎo)致吸力面流體流動(dòng)紊亂,翼型氣動(dòng)性能變差,升阻比下降。當(dāng)振幅增大為0.025c或0.03c時(shí),振動(dòng)頻率是小頻率的微小板,對(duì)翼型的氣動(dòng)性能有較好的控制效果;隨著頻率增大,翼型升阻比先增大后減小。圖11(c)中,增大板長(zhǎng)為0.02c時(shí),微小板尺寸大,靜止時(shí)對(duì)流場(chǎng)的擾動(dòng)大,振動(dòng)時(shí)會(huì)在板的上下表面形成正負(fù)壓區(qū)域,影響了翼型上表面的負(fù)壓區(qū),導(dǎo)致翼型吸力面的流體流動(dòng)更紊亂,翼型氣動(dòng)性能變差,升阻比反而沒有提高。此外,4 種尺寸振動(dòng)微小板所對(duì)應(yīng)的最佳設(shè)置參數(shù)如表3 所示,板長(zhǎng)越小,翼型升阻比提升比率越大。

    通過分析翼型附近流場(chǎng)的渦量流線圖,可以進(jìn)一步探討振動(dòng)微小板對(duì)翼型流動(dòng)控制的機(jī)理。圖12 是靜止或振動(dòng)的微小板對(duì)翼型流場(chǎng)影響的對(duì)比圖,圖中微小板振動(dòng)相位與靜止微小板在Y軸上位置一致。

    表3 不同板長(zhǎng)最佳振動(dòng)微小板參數(shù)Tab.3 Data sheets of best oscillating tiny sheets with different sheet lengths

    圖12 翼型渦量流線圖Fig.12 Streamline-vorticity diagram of the S809 airfoil

    圖12(a)和圖12(b)中顯示翼型吸力面上整體呈順時(shí)針渦量,翼型分離泡大小相似,原因是圖12(b)中微小板振動(dòng)參數(shù)是A=0.003c,f=0.5,小頻率小振幅的振動(dòng)板與靜止板控制流動(dòng)效果相似;增大頻率,f=5 時(shí),振動(dòng)微小板產(chǎn)生2 個(gè)細(xì)長(zhǎng)的正負(fù)渦量,微小板的渦對(duì)翼型的旋渦影響很大,吸力面上分離泡明顯小于圖12(a)加設(shè)靜止微小板控制翼型時(shí)的分離泡,該工況下的振動(dòng)微小板能有效控制流場(chǎng)流動(dòng)。

    圖12(d)和圖12(e)中顯示翼型吸力面上分離泡大小相似,原因是圖12(e)中微小板振動(dòng)參數(shù)是A=0.005c,f=1.5,小頻率小振幅的振動(dòng)板與靜止板控制流動(dòng)效果一致;增加振幅A=0.03c,由圖12(f)中可以看出,流體的流線緊貼翼型上表面前半部分,且有著薄薄的一層較低值的渦量層,脫落的尾渦較小,該工況下的振動(dòng)微小板能有效控制流場(chǎng)流動(dòng)。

    增大板長(zhǎng)為0.02c,由圖12(g)~12(i)可知,振動(dòng)微小板對(duì)翼型的流動(dòng)有阻礙作用,翼型表面出現(xiàn)多個(gè)分離渦,在靠近尾緣處的分離渦明顯大于靜止板控制翼型時(shí)產(chǎn)生的尾渦,振動(dòng)會(huì)惡化翼型吸力面流體流動(dòng),導(dǎo)致氣動(dòng)性能變差。該工況下的微小板振動(dòng)振幅越大,翼型氣動(dòng)性能越差,微小板控制翼型流場(chǎng)流動(dòng)能力越差。

    綜上,小尺寸微小板賦予合適的振動(dòng)能控制翼型流場(chǎng)流動(dòng)分離,提高翼型的升阻比,改善翼型氣動(dòng)性能。增大板長(zhǎng)為0.02c時(shí),振動(dòng)微小板對(duì)翼型氣動(dòng)性能有消極影響,原因是微小板尺寸大,振動(dòng)會(huì)惡化翼型吸力面流體流動(dòng),導(dǎo)致翼型氣動(dòng)性能變差。

    3 結(jié) 論

    S809 翼型表面附近設(shè)置靜止微小板裝置,具有一定的流動(dòng)分離控制效果;在考察的工況范圍內(nèi),隨著微小板板長(zhǎng)的增加,流動(dòng)分離控制效果更明顯,同時(shí)也受到板位置的影響,微小板安裝在適當(dāng)?shù)奈恢貌拍苓_(dá)到較好的控制效果。

    對(duì)抑制流動(dòng)分離效果不佳的尺寸較小的靜止板施加以合適振幅和頻率的振動(dòng)后,可明顯減小流動(dòng)分離泡,提升升力系數(shù),同時(shí)可以有效地降低阻力系數(shù),使得翼型的氣動(dòng)性能得到顯著提升。

    猜你喜歡
    吸力振幅氣動(dòng)
    深水大型吸力錨測(cè)試技術(shù)
    中寰氣動(dòng)執(zhí)行機(jī)構(gòu)
    ROV在海上吸力樁安裝場(chǎng)景的應(yīng)用及安裝精度和風(fēng)險(xiǎn)控制
    化工管理(2022年11期)2022-06-03 07:08:24
    基于NACA0030的波紋狀翼型氣動(dòng)特性探索
    深水吸力樁施工技術(shù)研究
    基于反饋線性化的RLV氣動(dòng)控制一體化設(shè)計(jì)
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    滬市十大振幅
    人人妻人人爽人人添夜夜欢视频 | 大又大粗又爽又黄少妇毛片口| 日本熟妇午夜| 日本与韩国留学比较| 国产精品一及| 一级毛片我不卡| 一级av片app| 如何舔出高潮| 在线观看一区二区三区激情| 大陆偷拍与自拍| 国产精品国产三级专区第一集| 亚洲自偷自拍三级| 亚洲欧美清纯卡通| 成人鲁丝片一二三区免费| 男女无遮挡免费网站观看| 久久女婷五月综合色啪小说 | 国产亚洲精品久久久com| 国产爽快片一区二区三区| 亚洲,欧美,日韩| 日韩av在线免费看完整版不卡| 97超视频在线观看视频| 日韩欧美精品免费久久| 精品亚洲乱码少妇综合久久| 成人免费观看视频高清| 亚洲自拍偷在线| 激情 狠狠 欧美| 午夜福利视频1000在线观看| 黄色欧美视频在线观看| 熟女电影av网| 亚洲综合色惰| 毛片女人毛片| 涩涩av久久男人的天堂| 国产免费福利视频在线观看| 日韩av免费高清视频| 欧美另类一区| 国产探花极品一区二区| 久久久久国产精品人妻一区二区| 一二三四中文在线观看免费高清| 成人特级av手机在线观看| 免费观看a级毛片全部| 国产一区二区三区av在线| 亚洲精品日韩av片在线观看| 亚洲av成人精品一二三区| 精品久久久久久电影网| 五月玫瑰六月丁香| 一级毛片久久久久久久久女| 亚洲欧美中文字幕日韩二区| 成人二区视频| 亚洲国产精品国产精品| 欧美潮喷喷水| 精品少妇黑人巨大在线播放| .国产精品久久| 精品人妻偷拍中文字幕| 国产 一区精品| 日韩大片免费观看网站| 夫妻午夜视频| 卡戴珊不雅视频在线播放| 国语对白做爰xxxⅹ性视频网站| 久久精品国产鲁丝片午夜精品| 成人亚洲精品一区在线观看 | 日韩一本色道免费dvd| 午夜激情久久久久久久| 狂野欧美白嫩少妇大欣赏| 日韩欧美一区视频在线观看 | 国产成人精品福利久久| 国产成人免费观看mmmm| 精品一区二区三区视频在线| 免费av观看视频| 亚洲av中文av极速乱| 欧美另类一区| 国产精品久久久久久久久免| 亚洲av二区三区四区| 亚洲av成人精品一二三区| 六月丁香七月| 精品久久久噜噜| 制服丝袜香蕉在线| 国产亚洲最大av| 狂野欧美白嫩少妇大欣赏| 小蜜桃在线观看免费完整版高清| 亚洲自拍偷在线| 日本爱情动作片www.在线观看| 嘟嘟电影网在线观看| 国产成人freesex在线| 久久国产乱子免费精品| 精品久久国产蜜桃| 欧美日韩国产mv在线观看视频 | 一级毛片久久久久久久久女| 性色avwww在线观看| 熟女人妻精品中文字幕| 女的被弄到高潮叫床怎么办| 国产淫语在线视频| 九草在线视频观看| 在线看a的网站| 亚洲精品成人久久久久久| 国产永久视频网站| 在线观看av片永久免费下载| 免费观看a级毛片全部| 亚洲熟女精品中文字幕| 日本与韩国留学比较| 精品人妻视频免费看| 亚洲av成人精品一二三区| av国产久精品久网站免费入址| 极品少妇高潮喷水抽搐| 亚洲第一区二区三区不卡| 大码成人一级视频| 青春草国产在线视频| 在线精品无人区一区二区三 | 久久99热这里只频精品6学生| 国产v大片淫在线免费观看| 黄色一级大片看看| 日日撸夜夜添| 亚洲人与动物交配视频| 高清视频免费观看一区二区| 日韩不卡一区二区三区视频在线| 亚洲精品日本国产第一区| 国产伦精品一区二区三区四那| 亚洲精品自拍成人| 最近中文字幕高清免费大全6| 亚洲国产成人一精品久久久| 国产免费福利视频在线观看| 九九在线视频观看精品| 最近最新中文字幕免费大全7| 色播亚洲综合网| 在线免费观看不下载黄p国产| 久久久a久久爽久久v久久| 亚洲第一区二区三区不卡| 国产黄a三级三级三级人| 欧美xxxx性猛交bbbb| 97在线人人人人妻| 不卡视频在线观看欧美| 久久精品人妻少妇| 3wmmmm亚洲av在线观看| 国产男人的电影天堂91| 嫩草影院新地址| 黄色视频在线播放观看不卡| 亚洲精品乱码久久久v下载方式| 看黄色毛片网站| 你懂的网址亚洲精品在线观看| 欧美日韩视频精品一区| 一边亲一边摸免费视频| 亚洲aⅴ乱码一区二区在线播放| 亚洲人与动物交配视频| 国产亚洲av片在线观看秒播厂| 久久亚洲国产成人精品v| 国产女主播在线喷水免费视频网站| 国产男女内射视频| 成年av动漫网址| 国产av码专区亚洲av| 五月玫瑰六月丁香| 特大巨黑吊av在线直播| 色网站视频免费| 大片电影免费在线观看免费| 免费看av在线观看网站| 亚洲第一区二区三区不卡| 国国产精品蜜臀av免费| 日本三级黄在线观看| 亚洲欧美一区二区三区黑人 | 精品久久久久久久久亚洲| av国产久精品久网站免费入址| 天堂网av新在线| 日日啪夜夜撸| 亚洲高清免费不卡视频| 国产精品人妻久久久久久| av国产免费在线观看| 久久午夜福利片| 天堂网av新在线| 汤姆久久久久久久影院中文字幕| 一本久久精品| 久热这里只有精品99| 亚洲精品,欧美精品| 少妇猛男粗大的猛烈进出视频 | av一本久久久久| 国产精品蜜桃在线观看| 我的老师免费观看完整版| 在线观看国产h片| 如何舔出高潮| 国产乱人视频| 97热精品久久久久久| 成年女人在线观看亚洲视频 | 亚洲精品一区蜜桃| 国产亚洲5aaaaa淫片| av卡一久久| 精品久久久噜噜| 国产精品偷伦视频观看了| 六月丁香七月| 国产伦精品一区二区三区四那| 欧美日韩视频精品一区| 人妻少妇偷人精品九色| 国产精品一区二区在线观看99| 午夜激情久久久久久久| 日本三级黄在线观看| 欧美精品人与动牲交sv欧美| 大片免费播放器 马上看| 黄色欧美视频在线观看| 精品国产三级普通话版| 秋霞在线观看毛片| av.在线天堂| 久久久精品欧美日韩精品| 内地一区二区视频在线| 男女边摸边吃奶| 卡戴珊不雅视频在线播放| 大又大粗又爽又黄少妇毛片口| 街头女战士在线观看网站| 日韩欧美精品免费久久| 天天一区二区日本电影三级| 97在线视频观看| 午夜精品国产一区二区电影 | 女人十人毛片免费观看3o分钟| 亚洲人成网站高清观看| 免费播放大片免费观看视频在线观看| 99热6这里只有精品| 看免费成人av毛片| 五月玫瑰六月丁香| 晚上一个人看的免费电影| 制服丝袜香蕉在线| 亚洲高清免费不卡视频| 内射极品少妇av片p| 日韩三级伦理在线观看| 免费黄网站久久成人精品| 交换朋友夫妻互换小说| 欧美zozozo另类| 成人综合一区亚洲| 亚洲综合色惰| 18禁动态无遮挡网站| 天堂中文最新版在线下载 | 成年女人在线观看亚洲视频 | 免费少妇av软件| av播播在线观看一区| 欧美xxxx黑人xx丫x性爽| 久久久精品免费免费高清| 一区二区av电影网| 97精品久久久久久久久久精品| 亚洲av二区三区四区| 精品一区在线观看国产| 国产成人a区在线观看| 少妇人妻一区二区三区视频| 黄色视频在线播放观看不卡| 亚洲第一区二区三区不卡| 日韩成人av中文字幕在线观看| 丝袜喷水一区| 亚洲欧洲国产日韩| 人妻系列 视频| 亚洲国产高清在线一区二区三| 欧美一区二区亚洲| 国产探花在线观看一区二区| 亚洲av成人精品一二三区| 在线免费观看不下载黄p国产| 色网站视频免费| 18禁在线播放成人免费| tube8黄色片| 国产午夜精品久久久久久一区二区三区| 性色avwww在线观看| 99re6热这里在线精品视频| 亚洲国产欧美在线一区| 亚洲精品国产av蜜桃| 中文天堂在线官网| 97热精品久久久久久| 麻豆成人av视频| 一级片'在线观看视频| 国产91av在线免费观看| 欧美日韩视频精品一区| 久久精品国产亚洲网站| 天堂中文最新版在线下载 | 午夜免费男女啪啪视频观看| 亚洲av福利一区| 蜜桃久久精品国产亚洲av| 亚洲精品日韩在线中文字幕| 精品一区二区三区视频在线| 一级毛片黄色毛片免费观看视频| 欧美xxxx性猛交bbbb| 99久久九九国产精品国产免费| 久久久久久久午夜电影| 国产成人a区在线观看| 少妇丰满av| 秋霞在线观看毛片| 天天躁夜夜躁狠狠久久av| 王馨瑶露胸无遮挡在线观看| 国产精品一区二区三区四区免费观看| 99九九线精品视频在线观看视频| 另类亚洲欧美激情| 一级片'在线观看视频| 别揉我奶头 嗯啊视频| 欧美日韩在线观看h| 国产精品久久久久久精品电影小说 | 观看免费一级毛片| 日韩制服骚丝袜av| 久久久久久久国产电影| 日本色播在线视频| 欧美3d第一页| 亚洲激情五月婷婷啪啪| 午夜福利视频1000在线观看| 老司机影院毛片| 精品少妇久久久久久888优播| 午夜免费观看性视频| 亚洲av男天堂| 久久精品国产亚洲av天美| 亚洲成人精品中文字幕电影| 麻豆乱淫一区二区| 亚洲精品亚洲一区二区| 国产成人aa在线观看| 免费观看av网站的网址| 精品久久久精品久久久| 黄色日韩在线| 91久久精品国产一区二区成人| 久久精品熟女亚洲av麻豆精品| 99热网站在线观看| 亚洲精品一区蜜桃| 久久99热这里只有精品18| 成人二区视频| 亚洲成人精品中文字幕电影| 日本黄色片子视频| 亚洲真实伦在线观看| 国产精品一区www在线观看| 国产有黄有色有爽视频| 97热精品久久久久久| 九草在线视频观看| 亚洲色图综合在线观看| 亚洲欧美日韩东京热| a级毛片免费高清观看在线播放| tube8黄色片| 草草在线视频免费看| 国产亚洲一区二区精品| 国产精品福利在线免费观看| 最近2019中文字幕mv第一页| 99视频精品全部免费 在线| 高清av免费在线| 丝瓜视频免费看黄片| 免费观看a级毛片全部| 国产精品一二三区在线看| 国内少妇人妻偷人精品xxx网站| 亚洲av免费高清在线观看| 大陆偷拍与自拍| av卡一久久| 岛国毛片在线播放| 日日摸夜夜添夜夜爱| 国产精品不卡视频一区二区| 舔av片在线| 80岁老熟妇乱子伦牲交| 欧美少妇被猛烈插入视频| 夜夜爽夜夜爽视频| av又黄又爽大尺度在线免费看| 又爽又黄a免费视频| 精品国产三级普通话版| 亚洲av在线观看美女高潮| 日韩中字成人| 国产爽快片一区二区三区| 精品久久久久久久久亚洲| 男人狂女人下面高潮的视频| 亚洲精品日韩av片在线观看| 欧美日韩亚洲高清精品| 美女主播在线视频| 美女高潮的动态| 国产精品国产三级专区第一集| 久久久久国产精品人妻一区二区| 国产黄色免费在线视频| 青春草国产在线视频| 国内少妇人妻偷人精品xxx网站| 亚洲国产精品成人综合色| 联通29元200g的流量卡| 亚洲精品日韩av片在线观看| 亚洲综合精品二区| 精品久久久久久久人妻蜜臀av| 久久99蜜桃精品久久| 国产免费视频播放在线视频| 日韩视频在线欧美| 国产精品无大码| 最近手机中文字幕大全| 波多野结衣巨乳人妻| 成人国产麻豆网| 丝袜脚勾引网站| 久久鲁丝午夜福利片| 国产成人精品婷婷| 成人二区视频| 1000部很黄的大片| 免费观看性生交大片5| 可以在线观看毛片的网站| 在线观看一区二区三区| 少妇 在线观看| 国产老妇女一区| 18禁在线播放成人免费| 内射极品少妇av片p| 天美传媒精品一区二区| 少妇的逼水好多| 国产精品av视频在线免费观看| 国产淫片久久久久久久久| 肉色欧美久久久久久久蜜桃 | 欧美精品人与动牲交sv欧美| 免费黄色在线免费观看| 亚洲欧美精品专区久久| 久久久久久久大尺度免费视频| 国产一区二区在线观看日韩| 亚洲精品一区蜜桃| 亚州av有码| 成人美女网站在线观看视频| 有码 亚洲区| 小蜜桃在线观看免费完整版高清| 国产精品女同一区二区软件| 成人高潮视频无遮挡免费网站| 欧美激情在线99| 成年版毛片免费区| 久久久久九九精品影院| 国产免费一级a男人的天堂| 麻豆成人午夜福利视频| 日日啪夜夜撸| 欧美成人精品欧美一级黄| 少妇被粗大猛烈的视频| 国产乱人偷精品视频| 久久久久久伊人网av| 在线观看人妻少妇| 亚洲内射少妇av| 波多野结衣巨乳人妻| 久久影院123| 99re6热这里在线精品视频| 久久97久久精品| 成人亚洲精品一区在线观看 | 下体分泌物呈黄色| 亚洲精品久久久久久婷婷小说| 麻豆精品久久久久久蜜桃| av播播在线观看一区| 极品教师在线视频| 99久久精品热视频| 亚洲在线观看片| 2018国产大陆天天弄谢| 欧美日本视频| 国精品久久久久久国模美| 久久女婷五月综合色啪小说 | 岛国毛片在线播放| 18禁裸乳无遮挡免费网站照片| 免费观看av网站的网址| 最近2019中文字幕mv第一页| 国产一区二区三区av在线| 国产欧美日韩一区二区三区在线 | 丝袜喷水一区| 久久久久久国产a免费观看| 国产成人午夜福利电影在线观看| 久久精品久久精品一区二区三区| 最近中文字幕2019免费版| 热99国产精品久久久久久7| 精品人妻偷拍中文字幕| 男人狂女人下面高潮的视频| 菩萨蛮人人尽说江南好唐韦庄| 在线观看人妻少妇| 欧美国产精品一级二级三级 | 免费少妇av软件| 国产一区二区亚洲精品在线观看| 久久久久久久久久成人| 王馨瑶露胸无遮挡在线观看| 两个人的视频大全免费| 国产亚洲一区二区精品| 国产精品av视频在线免费观看| 久久久a久久爽久久v久久| 国产黄a三级三级三级人| 欧美老熟妇乱子伦牲交| 卡戴珊不雅视频在线播放| 一级a做视频免费观看| 免费看不卡的av| 国产黄片视频在线免费观看| 久久久欧美国产精品| 欧美少妇被猛烈插入视频| 人人妻人人爽人人添夜夜欢视频 | 日本色播在线视频| 涩涩av久久男人的天堂| 亚洲一级一片aⅴ在线观看| 大又大粗又爽又黄少妇毛片口| 狂野欧美激情性xxxx在线观看| 内射极品少妇av片p| 王馨瑶露胸无遮挡在线观看| 91久久精品电影网| 蜜桃亚洲精品一区二区三区| av免费在线看不卡| 人妻 亚洲 视频| 黄色怎么调成土黄色| 在线a可以看的网站| 亚洲国产精品成人久久小说| 国产午夜精品久久久久久一区二区三区| 亚洲人成网站在线播| 一本色道久久久久久精品综合| 日韩一区二区视频免费看| 欧美性感艳星| 国产精品一区二区在线观看99| 久久久色成人| 国产精品无大码| 日日啪夜夜爽| 欧美国产精品一级二级三级 | 成年女人在线观看亚洲视频 | 美女脱内裤让男人舔精品视频| 久久久久精品性色| 日本色播在线视频| 又黄又爽又刺激的免费视频.| 有码 亚洲区| 久久精品国产亚洲av涩爱| 欧美成人a在线观看| 1000部很黄的大片| 91精品国产九色| 18禁裸乳无遮挡动漫免费视频 | 国产精品女同一区二区软件| 人妻一区二区av| 少妇裸体淫交视频免费看高清| 国产精品99久久久久久久久| 免费观看无遮挡的男女| 在线看a的网站| 国产69精品久久久久777片| 久久99热这里只有精品18| 一级毛片aaaaaa免费看小| 久久久久久久精品精品| 天堂中文最新版在线下载 | 国产精品一区二区性色av| 成年免费大片在线观看| 美女国产视频在线观看| 精品久久久久久久人妻蜜臀av| 亚洲,一卡二卡三卡| 亚洲欧美一区二区三区黑人 | 国产在视频线精品| 国产亚洲5aaaaa淫片| av线在线观看网站| 久久久久久久午夜电影| 久久久久久久久久久免费av| 久久久久精品性色| 国产一区二区在线观看日韩| 大又大粗又爽又黄少妇毛片口| 综合色丁香网| 天堂中文最新版在线下载 | 国产高潮美女av| 久热久热在线精品观看| 在线观看三级黄色| 国产成人freesex在线| 99久国产av精品国产电影| 免费看a级黄色片| 伊人久久国产一区二区| 国产高清国产精品国产三级 | 国产精品久久久久久精品电影| 成年女人在线观看亚洲视频 | 亚洲欧美精品自产自拍| 干丝袜人妻中文字幕| 丰满人妻一区二区三区视频av| 久久亚洲国产成人精品v| 能在线免费看毛片的网站| 午夜福利在线在线| 国产日韩欧美亚洲二区| 在线观看国产h片| 人体艺术视频欧美日本| 天堂俺去俺来也www色官网| 我的女老师完整版在线观看| 久久亚洲国产成人精品v| 插阴视频在线观看视频| 只有这里有精品99| 真实男女啪啪啪动态图| 亚洲av一区综合| 制服丝袜香蕉在线| 91精品国产九色| 啦啦啦在线观看免费高清www| 国产一区二区三区综合在线观看 | 精品99又大又爽又粗少妇毛片| 国产国拍精品亚洲av在线观看| 免费大片黄手机在线观看| 中文乱码字字幕精品一区二区三区| 成人黄色视频免费在线看| 在线观看一区二区三区激情| 免费观看性生交大片5| 边亲边吃奶的免费视频| 国产91av在线免费观看| 丰满乱子伦码专区| 久久久久久久久久人人人人人人| 麻豆精品久久久久久蜜桃| 亚洲一级一片aⅴ在线观看| 尾随美女入室| 2021少妇久久久久久久久久久| 97人妻精品一区二区三区麻豆| 一级爰片在线观看| 免费观看性生交大片5| 国产高潮美女av| 亚洲美女搞黄在线观看| 激情 狠狠 欧美| av天堂中文字幕网| 青春草国产在线视频| 国产欧美另类精品又又久久亚洲欧美| 日本熟妇午夜| 熟女人妻精品中文字幕| 国产成人91sexporn| 中国美白少妇内射xxxbb| 禁无遮挡网站| 亚洲自拍偷在线| 免费黄频网站在线观看国产| 免费观看性生交大片5| 22中文网久久字幕| 欧美成人午夜免费资源| 中国国产av一级| 亚洲欧美精品自产自拍| 狠狠精品人妻久久久久久综合| 男女边吃奶边做爰视频| 国产男女超爽视频在线观看| 麻豆精品久久久久久蜜桃| 免费不卡的大黄色大毛片视频在线观看| 免费看光身美女| 国产黄色视频一区二区在线观看| 亚洲av男天堂| 国产精品99久久久久久久久| 国产精品嫩草影院av在线观看| 最新中文字幕久久久久| 小蜜桃在线观看免费完整版高清| 插阴视频在线观看视频| 日本午夜av视频| 国产在线男女| 日本一二三区视频观看| 精品国产乱码久久久久久小说| 午夜免费男女啪啪视频观看| 免费观看的影片在线观看| eeuss影院久久| 另类亚洲欧美激情| 一本久久精品| eeuss影院久久| 欧美高清成人免费视频www| 99久国产av精品国产电影| 看非洲黑人一级黄片| 亚洲欧美成人精品一区二区|