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

    后退式微型后緣裝置對(duì)翼型氣動(dòng)特性影響的實(shí)驗(yàn)研究

    2015-06-21 15:08:47駿,欽,2,*
    實(shí)驗(yàn)流體力學(xué) 2015年5期
    關(guān)鍵詞:翼面襟翼后緣

    夏 駿, 代 欽,2,*

    (1. 上海大學(xué) 上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所, 上海 200072; 2.上海市力學(xué)在能源工程中的應(yīng)用重點(diǎn)實(shí)驗(yàn)室, 上海 200072)

    后退式微型后緣裝置對(duì)翼型氣動(dòng)特性影響的實(shí)驗(yàn)研究

    夏 駿1, 代 欽1,2,*

    (1. 上海大學(xué) 上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所, 上海 200072; 2.上海市力學(xué)在能源工程中的應(yīng)用重點(diǎn)實(shí)驗(yàn)室, 上海 200072)

    介紹了裝有后退式微型后緣裝置(Rearward Mini-TED)的NACA23012翼型在低雷諾數(shù)條件下的表面壓力分布、氣動(dòng)力和PIV速度場(chǎng)的風(fēng)洞實(shí)驗(yàn)結(jié)果,并與NACA23012原型翼的對(duì)應(yīng)測(cè)量結(jié)果進(jìn)行了對(duì)比分析, 以探討Mini-TED裝置對(duì)翼型流場(chǎng)、氣動(dòng)特性產(chǎn)生的影響。本實(shí)驗(yàn)風(fēng)速為15m/s,以弦長(zhǎng)為特征量的雷諾數(shù)為Re≈1.3×105,翼型表面壓力分布采用測(cè)壓孔和壓力傳感器測(cè)量,通過(guò)積分獲得翼型升力和壓差阻力,并利用尾耙測(cè)量翼型受到的總阻力。結(jié)果表明,后退式 Mini-TED翼型改變了翼型周?chē)牧鲌?chǎng)速度分布和尾流流動(dòng)結(jié)構(gòu),導(dǎo)致上翼面吸力和下翼面的壓力升高,使翼型升力增加,但壓差阻力也增加。同時(shí)發(fā)現(xiàn)后退式 Mini-TED翼型使前駐點(diǎn)位置后移,加快了上翼面的流動(dòng)速度,后緣分離受到抑制。

    后退式Mini-TED;翼型氣動(dòng)特性;低雷諾數(shù);表面壓力分布;PIV測(cè)量

    0 引 言

    微型后緣裝置(Mini-TED)是一種后緣流動(dòng)控制裝置,其幾何形態(tài)和安裝參數(shù)的變化能夠有效地改變翼型表面的壓力分布,增大翼型的升力和升阻比。由于安裝使用簡(jiǎn)易,增升效果明顯,因而受到研究人員的關(guān)注。Mini-TED種類(lèi)繁多,有Gurney襟翼、開(kāi)裂式襟翼、發(fā)散后緣和楔形襟翼等(如圖1所示,其中圖(a)~(e)引自文獻(xiàn)[1])。文獻(xiàn)[2]將Mini-TED看作是一類(lèi)具有不同幾何形態(tài)的后緣控制裝置的統(tǒng)稱(chēng),尺寸比傳統(tǒng)后緣流動(dòng)控制裝置小,通常限制在弦長(zhǎng)的2%以下。

    圖1 幾種不同形態(tài)的Mini-TED

    K.Richter和H.Roseman[2-3]對(duì)裝有Mini-TED的超臨界翼型在跨聲速流動(dòng)中的氣動(dòng)特性進(jìn)行總覽性的研究,從實(shí)驗(yàn)研究和數(shù)值模擬2方面,對(duì)Gurney襟翼、開(kāi)裂式襟翼和發(fā)散式襟翼3種形態(tài)的Mini-TED及其尺寸對(duì)翼面壓力分布、氣動(dòng)力特性的影響進(jìn)行了對(duì)比分析。研究認(rèn)為,Mini-TED使翼型后緣彎度提高,翼型后段上下翼面壓力差增大,從而產(chǎn)生增升效果。Gurney襟翼的長(zhǎng)度變化與開(kāi)裂襟翼的安裝角變化對(duì)翼型氣動(dòng)特性的影響類(lèi)似,但相同長(zhǎng)度條件下的開(kāi)裂襟翼增升效果較小,而發(fā)散式后緣與開(kāi)裂式襟翼對(duì)翼型具有同樣的氣動(dòng)影響。該文獻(xiàn)根據(jù)3種Mini-TED氣動(dòng)效率比較分析的結(jié)果,認(rèn)為Gurney襟翼和開(kāi)裂式襟翼適于發(fā)展為自適應(yīng)流動(dòng)控制裝置(Adaptive Flow Control Device)和阻力優(yōu)化設(shè)計(jì)。由于該文偏重于超臨界翼型流動(dòng)控制的工程應(yīng)用,因此,如能夠進(jìn)一步開(kāi)展對(duì)流場(chǎng)結(jié)構(gòu)的測(cè)量及其與氣動(dòng)力之間關(guān)系的討論,對(duì)于深入理解Mini-TED增升效果的機(jī)理將十分有幫助。

    Liebeck[4]在關(guān)于翼型增升的研究中,將Gurney襟翼作為高升力翼型的一種附加裝置進(jìn)行實(shí)驗(yàn)研究。在Newman翼型安裝了高度為1.25%c(c為弦長(zhǎng))的Gurney襟翼后,發(fā)現(xiàn)升力及最大升力系數(shù)比原翼型有較大的提升, 并伴隨零升迎角和同等升力系數(shù)下阻力的減小。

    王晉軍等對(duì)Gurney襟翼進(jìn)行了綜述性的回顧[5],并詳細(xì)研究了Gurney襟翼的安裝參數(shù)對(duì)低速翼型的氣動(dòng)影響和增升機(jī)理[6]。研究中發(fā)現(xiàn),隨著襟翼安裝角的增加,翼型的升力和阻力也有所增加;對(duì)于縮進(jìn)式Gurney襟翼,升力的增量隨著縮進(jìn)量的增加而減小。其中關(guān)于安裝角小于90°的縮進(jìn)式Gurney襟翼的討論可以作為開(kāi)裂式襟翼研究的參考。

    歐洲空中客車(chē)公司的AWIATOR計(jì)劃中將微型后緣裝置(Mini-TED)作為“自適應(yīng)翼型”的一種形式設(shè)計(jì)進(jìn)行研究[7-9]。其中A.D. Gardner[8]等人認(rèn)為小型開(kāi)裂式Mini-TED既可以在巡航狀態(tài)時(shí)作為發(fā)散后緣使用,也可以在低速狀態(tài)下起到與Gurney襟翼類(lèi)似的增升作用。

    諸胡冰[10]等人利用數(shù)值計(jì)算的方法研究了開(kāi)裂式Mini-TED安裝參數(shù)對(duì)翼型氣動(dòng)特性的影響以及后緣渦結(jié)構(gòu)的特性。該研究認(rèn)為,開(kāi)裂式Mini-TED打開(kāi)后,后緣形成反向駐渦,誘導(dǎo)上翼面氣流加速下偏,擴(kuò)展尾跡,加長(zhǎng)“氣動(dòng)弦長(zhǎng)”,使升力線整體上移,最大升力和升阻比提高,由于襟翼長(zhǎng)度不同,造成襟翼埋入邊界層的厚度不同,從而對(duì)升阻特性產(chǎn)生不同的影響。A.W.Bloy[1]等人對(duì)5種不同形態(tài)的Mini-TED在相同襟翼長(zhǎng)度的條件下進(jìn)行實(shí)驗(yàn)研究,包括Gurney襟翼、有安裝角的Gurney襟翼、45°楔形襟翼、發(fā)散后緣以及方截面襟翼(分別對(duì)應(yīng)了圖1中(a)、(b)、(c)、(d)和(e) 5種襟翼)。研究表明,所有形態(tài)的Mini-TED都能起到增升的作用,但是45°楔形襟翼能夠在最大升力系數(shù)提升的同時(shí)比Gurney襟翼受到更小的阻力,在中低升力狀態(tài)下達(dá)到更高的升阻比。

    周華[11]使用計(jì)算方法研究了加裝開(kāi)裂式Mini-TED后NACA0012翼型在Ma=0.8時(shí)的跨聲速氣動(dòng)特性和流場(chǎng)結(jié)構(gòu)。研究認(rèn)為,開(kāi)裂式Mini-TED對(duì)翼型的增升機(jī)理不能簡(jiǎn)單地歸因于加裝襟翼造成后緣彎度的增加,而是Mini-TED在后緣后產(chǎn)生的駐渦結(jié)構(gòu)改變了后緣附近的庫(kù)塔條件,并導(dǎo)致激波位置的大幅度后移。在進(jìn)一步的研究[12]中,周華認(rèn)為,跨聲速流動(dòng)中,開(kāi)裂式Mini-TED后產(chǎn)生的三渦結(jié)構(gòu)穩(wěn)定性不足,因而提出后退式Mini-TED (Rearward Mini-TED)設(shè)計(jì)概念,即將開(kāi)裂式Mini-TED在打開(kāi)的同時(shí)向后緣移動(dòng),伸出后緣形成后向臺(tái)階結(jié)構(gòu),以獲得更加穩(wěn)定的后緣駐渦。與在原設(shè)計(jì)條件下的翼型相比,安裝后退式Mini-TED的NACA0012翼型在亞跨聲速時(shí)與開(kāi)裂式Mini-TED升力接近而阻力較??;在跨聲速條件下,安裝后退式Mini-TED的NACA0012翼型在同等條件下激波位置后移,致使升力曲線以及升阻比曲線相比開(kāi)裂式Mini-TED翼型更加規(guī)則,并具有更高的升力及升阻比。

    由于對(duì)于后退式Mini-TED機(jī)翼的研究目前較少,僅限于以對(duì)稱(chēng)翼型在大雷諾數(shù)和亞跨聲速條件下的數(shù)值模擬結(jié)果,尚無(wú)實(shí)驗(yàn)研究的報(bào)道。同時(shí)考慮到該Mini-TED裝置安裝簡(jiǎn)便,十分適用于在較小雷諾數(shù)下微型飛行器的增升作用。因此本實(shí)驗(yàn)針對(duì)加裝后退式Mini-TED的NACA23012翼型在低雷諾數(shù)條件下進(jìn)行表面壓力分布以及尾流壓力分布測(cè)量,從而得到翼型表面的壓力分布和升阻特性,并與原型翼型的實(shí)驗(yàn)結(jié)果相對(duì)比,以研究低速條件下后退式 Mini-TED對(duì)翼型的氣動(dòng)特性產(chǎn)生的影響,同時(shí)配合以PIV速度場(chǎng)測(cè)量,對(duì)壓力分布和氣動(dòng)力變化的原因進(jìn)行探討。

    1 實(shí)驗(yàn)裝置和條件

    表面壓力分布測(cè)量和PIV速度測(cè)量在上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所風(fēng)洞實(shí)驗(yàn)室中進(jìn)行,所使用風(fēng)洞為SIAMM回流式低速低湍度風(fēng)洞,實(shí)驗(yàn)段尺寸500mm×500mm×1 500mm,湍流度0.1%。實(shí)驗(yàn)布局如圖2(a)所示,翼型模型一端豎直固定在風(fēng)洞底部迎角調(diào)節(jié)機(jī)構(gòu)的電動(dòng)轉(zhuǎn)盤(pán)上,另一端加裝端板以降低三維流動(dòng)的影響。

    圖2 翼型測(cè)壓實(shí)驗(yàn)布置和后退式Mini-TED安裝示意

    Fig.2 Schematic diagram of the pressure measurement experimental setup and Rearward Mini-TED installation

    測(cè)壓實(shí)驗(yàn)基本翼型為NACA23012,表面光滑,弦長(zhǎng)c=120mm,在翼型模型表面共設(shè)置48個(gè)測(cè)壓孔,用以測(cè)量翼型表面的壓力分布;特別地,在翼型后緣布置有一個(gè)測(cè)壓孔,以確保上下翼面壓力數(shù)據(jù)在后緣點(diǎn)保持連續(xù),翼型升力通過(guò)對(duì)表面壓力分布進(jìn)行積分后獲得。后退式Mini-TED使用薄鋁片制作,厚度0.2mm,長(zhǎng)度5%c,安裝位置為下翼面距后緣點(diǎn)1.5%c處,安裝角與弦線成30°(沿弦線順時(shí)針?lè)较颍鐖D2(b)所示)。實(shí)驗(yàn)來(lái)流速度U∞=15m/s,以翼型弦長(zhǎng)為特征長(zhǎng)度的雷諾數(shù)Re=1.3×105,經(jīng)測(cè)試該鋁制Mini-TED有較好的強(qiáng)度,在本次實(shí)驗(yàn)氣流作用下可保持無(wú)變形和無(wú)顫振。實(shí)驗(yàn)測(cè)量對(duì)象為NACA23012基本翼型和加裝后退式Mini-TED的NACA23012翼型2種模型,以對(duì)比2者壓力分布的差異, 測(cè)壓實(shí)驗(yàn)工況迎角范圍從0°到18°,間隔1°。

    翼型阻力采用動(dòng)量法測(cè)量。在翼型后緣下游安裝尾耙,垂直于模型展向并位于翼型模型高度中央。尾耙測(cè)壓管布置如圖3所示,共設(shè)有51個(gè)測(cè)壓通道,靜壓管長(zhǎng)度50mm,總壓管長(zhǎng)度41mm,相鄰測(cè)壓管間距4mm。尾耙的測(cè)壓管連接到64通道壓力傳感器上進(jìn)行實(shí)時(shí)壓力采集,而后利用動(dòng)量定理計(jì)算翼型的阻力系數(shù)。

    圖3 尾耙測(cè)壓管布置

    PIV實(shí)驗(yàn)對(duì)前述2種翼型的上、下翼面和后緣附近Mini-TED周?chē)奈擦魉俣葓?chǎng)進(jìn)行了測(cè)量,以對(duì)比2翼型流動(dòng)結(jié)構(gòu)的差異,同時(shí)分析與壓力分布所對(duì)應(yīng)的流動(dòng)結(jié)構(gòu)之間的關(guān)聯(lián)性。測(cè)量截面與自由來(lái)流方向平行,并位于展長(zhǎng)中段,以保證足夠好的二維流動(dòng)特性。用于PIV測(cè)量的CCD相機(jī)分辨率為1600pixel×1200pixel,雙曝光工作模式;雙腔Nd:YAG激光器,單脈沖能量120mJ,重復(fù)頻率為15Hz,相機(jī)與激光器的同步由信號(hào)發(fā)生器控制。實(shí)驗(yàn)在每種迎角下連續(xù)采集500組粒子圖像后經(jīng)統(tǒng)計(jì)運(yùn)算得到時(shí)均速度場(chǎng)信息。

    2 實(shí)驗(yàn)結(jié)果分析

    2.1 壓力分布實(shí)驗(yàn)結(jié)果分析

    圖4為迎角0°、5°、10°和15°條件下,后退式 Mini-TED翼型與NACA23012原型翼型表面的壓力系數(shù)分布的對(duì)比,圖中Mini-TED代表帶后退式襟翼的翼型,CLEAN代表原型翼型。

    迎角0°時(shí)如圖4(a)所示,原型翼型上、下翼面前段70%弦長(zhǎng)均受到吸力作用,而下翼面后30%弦長(zhǎng)段壓力分布接近0值。后退式Mini-TED翼型上翼面負(fù)壓力值明顯低于原型翼型,下翼面前段的壓力分布與原型翼型幾乎重合;但從45%弦長(zhǎng)位置起,壓力值逐漸從負(fù)壓增加為正壓,下翼面后段受到正壓力作用;在90%弦長(zhǎng)處正壓達(dá)最大值。與原型翼型相比,Mini-TED翼型的壓力分布曲線包圍了較大的面積,因此翼型受到更大的升力作用。Mini-TED翼型上翼面壓力曲線從10%弦長(zhǎng)處向下游均勻下降至90%弦長(zhǎng)處,此后曲線略有抬升,使得翼面后段10%弦長(zhǎng)形成順壓梯度,即后緣點(diǎn)處受到吸力作用,明顯低于原型翼型后緣的壓力值,該吸力是由于機(jī)翼后緣形成的雙渦結(jié)構(gòu)的低壓區(qū)造成(參看后文尾流PIV速度場(chǎng)測(cè)量結(jié)果圖10,迎角0°與5°的雙渦結(jié)構(gòu)形態(tài)類(lèi)似,關(guān)于流場(chǎng)PIV測(cè)量分析的詳細(xì)內(nèi)容,作者已另外撰文討論[18])。同時(shí)在下翼面80%c~90%c段弦長(zhǎng),由于后退式Mini-TED對(duì)氣流產(chǎn)生的阻擋作用,造成局部的靜壓集中,使壓力分布在后緣前迅速升高,受到正壓力作用。后緣負(fù)壓區(qū)和Mini-TED下方的正壓區(qū)引起的壓強(qiáng)差,給翼型后緣附近帶來(lái)額外的升力增量,使翼型的低頭力矩增加。

    圖4 后退式 Mini-TED和原型翼型表面壓力系數(shù)

    Fig.4 Pressure coefficient distribution of rearward Mini-TED and NACA23012 airfoil

    迎角為5°時(shí)(見(jiàn)圖4(b)),Mini-TED翼型的上翼面全部受到吸力的作用,最大負(fù)壓系數(shù)已超過(guò)-1.0,幾乎為0°迎角時(shí)的1倍,吸力值遠(yuǎn)大于原型翼型,約為原型翼型吸力值的3倍,后緣點(diǎn)仍受到負(fù)壓作用;而原型翼型在后緣處壓力達(dá)到略大于0的正值。在下翼面,原型翼型只在翼面前段20%弦長(zhǎng)承受到正壓力,之后的壓力分布保持在0值附近;后退式Mini-TED翼型下翼面氣流受到襟翼阻擋,流速減緩靜壓增加,整個(gè)下翼面均受到正壓力作用,并且從35%弦長(zhǎng)起,壓力向下游近似呈線性增長(zhǎng),沒(méi)有出現(xiàn)迎角0°時(shí)后緣附近的壓力突增。上述2迎角時(shí),前駐點(diǎn)的位置均位于前緣點(diǎn)。

    10°迎角時(shí),如圖5(c)所示,2種翼型下翼面壓力均比5°迎角時(shí)增加。比較壓力最大值及其坐標(biāo)后發(fā)現(xiàn),Mini-TED翼型的前駐點(diǎn)比原型翼型稍遠(yuǎn)離前緣,因此氣流以更大的速度繞過(guò)前緣流向上翼面,造成上翼面的吸力峰增高,且分離泡尺度較小。Mini-TED翼型后緣點(diǎn)仍為負(fù)壓,而原型翼后緣點(diǎn)保持較低的正壓,Mini-TED翼型壓力分布曲線比原型翼型壓力分布曲線包圍更大的面積。當(dāng)迎角增大到15°時(shí)(見(jiàn)圖4(d)),原型翼型上翼面已經(jīng)發(fā)生了全翼面的流動(dòng)分離,整個(gè)上翼面壓力均勻分布;而Mini-TED翼型前駐點(diǎn)比10°迎角時(shí)更加遠(yuǎn)離前緣,上翼面最大流速進(jìn)一步增高,分離泡尺度減小,其內(nèi)部具有更大的負(fù)壓值,但仍未發(fā)生全翼面的流動(dòng)分離。最顯著的變化是原型翼的后緣壓力降為負(fù)值,且低于Mini-TED翼型的后緣壓力值,這是由于原型翼型上翼面發(fā)生了全翼面分離導(dǎo)致氣流不再減速,分離后的壓強(qiáng)接近分離點(diǎn)處壓強(qiáng);而Mini-TED翼型氣流前緣分離并再附后,翼型后段的流速持續(xù)下降,壓力得到恢復(fù),并從70%弦長(zhǎng)后超過(guò)原型翼型的壓力值。原型翼型下翼面壓力分布隨迎角變化不大,但Mini-TED翼型下翼面中段和后段的壓力隨著迎角的增大而不斷上升,壓力分布遠(yuǎn)高于5°和10°迎角時(shí)的情形,前駐點(diǎn)位置繼續(xù)向下游產(chǎn)生移動(dòng)。

    對(duì)比迎角從5°升高到15°的壓強(qiáng)變化可以看到,原型翼型上翼面負(fù)壓系數(shù)峰值的增量為0.15、0.2和-0.25,在迎角從10°升高到15°后,由于分離區(qū)覆蓋了全部上翼面,負(fù)壓值不增反降;而下翼面的壓力分布只在前30%范圍內(nèi)得到提升,后70%弦長(zhǎng)段壓強(qiáng)維持在0值附近,因此原型翼型升力的變化主要依靠吸力面的影響。同時(shí),由于下翼面的幾何輪廓發(fā)生了改變,流動(dòng)受阻增強(qiáng),從而整個(gè)下翼面的壓力隨著迎角增大而明顯升高。因此后退式Mini-TED翼型升力的增加來(lái)自下翼面壓力和上翼面吸力的共同貢獻(xiàn)。

    2.2 升/阻力結(jié)果分析

    圖5給出了由原型翼型和Mini-TED翼型表面壓力分布積分后得到的升力系數(shù)隨迎角變化的曲線,其中CLEAN代表原型翼型。與原型翼型相比,后退式Mini-TED翼型的CL-α曲線整體增長(zhǎng)顯著,升力系數(shù)提高,最大升力系數(shù)CLmax提升已達(dá)1.3,而原型翼型的CLmax僅為0.9;Mini-TED翼型升力系數(shù)曲線斜率也大于原型翼,即迎角越大則Mini-TED對(duì)升力增長(zhǎng)的貢獻(xiàn)越大。2翼型臨界迎角區(qū)別不大,均在17°迎角附近失速。

    圖6為2種翼型由表面壓力分布積分得到的壓差阻力系數(shù)以及利用尾耙測(cè)量得到的型阻隨迎角變化的曲線。0°迎角時(shí),Mini-TED翼型壓差阻力比原型翼型有微弱的提升,而型阻提升略多;隨著迎角增大,Mini-TED翼型比原型翼型迎風(fēng)面積增加,導(dǎo)致2種翼型阻力的差距不斷擴(kuò)大,Mini-TED翼型的阻力增長(zhǎng)速度更快,因此迎角為16°時(shí),翼型失速前Mini-TED翼型相對(duì)于原型翼總阻力的增量已超過(guò)壓差阻力的增量,即摩阻對(duì)阻力增量的貢獻(xiàn)隨著迎角的增加而增加。原型翼的壓差阻力與型阻之差(摩阻)幾乎不隨迎角發(fā)生變化,而Mini-TED 2壓力之差(摩阻)在11°迎角之前基本保持不變;在迎角超過(guò)11°后,摩阻逐漸增加,壓差阻力在型阻中所占比例有所下降。原因一方面在于Mini-TED使上翼面流速增加,邊界層內(nèi)有較高的法向速度梯度,具有更強(qiáng)的剪切運(yùn)動(dòng),從而摩擦力得到提高,這一推斷與較大迎角時(shí),Mini-TED翼型上面靜壓分布較低的實(shí)驗(yàn)結(jié)果相吻合。另一方面,前緣分離泡再附后,翼面邊界層通常已發(fā)展為湍流,而湍流邊界層的摩阻也較層流邊界層高,隨著迎角增加,分離泡尺度減小而湍流邊界層的區(qū)域擴(kuò)大。

    圖5 2種機(jī)翼升力系數(shù)隨迎角變化曲線

    圖6 2種機(jī)翼阻力系數(shù)隨迎角變化曲線

    Fig.6 Variation of the drag coefficientCDwith angles of attack(AOA)

    圖7(a)為2種翼型的升阻比隨升力系數(shù)變化的曲線。小升力系數(shù)時(shí)Mini-TED翼型升阻比小于原型翼型,當(dāng)升力系數(shù)高于0.35后,Mini-TED翼型的升阻比始終高于原型翼型,2翼型升阻比的差別也隨著升力系數(shù)的增大而不斷擴(kuò)大,在升力系數(shù)為0.9時(shí),Mini-TED翼型升阻比仍高于原型翼型,而2翼型的升阻比的極大值相差11%,并且Mini-TED翼型在升力系數(shù)0.6時(shí)具有最大升阻比,原型翼型僅在升力系數(shù)為0.3時(shí)具有最大升阻比,因此,后退式Mini-TED翼型具有更優(yōu)的升阻性能。兩種翼型的極曲線顯示(見(jiàn)圖7(b)),升力系數(shù)小于0.3時(shí),產(chǎn)生相同升力時(shí)原型翼型所受阻力較小,但與Mini-TED翼型相差不大。而當(dāng)升力系數(shù)大于0.3后,相同升力的Mini-TED翼型所受阻力低于原型翼型,并且隨著升力的增加,該阻力的差距逐漸擴(kuò)大。表明在中高升力條件下后退式Mini-TED翼型具有更好的阻力特性。

    圖7(a) 2種翼型升阻比隨升力變化曲線

    圖7(b) 2種翼型極曲線

    圖8為2種翼型對(duì)前緣的俯仰力矩系數(shù)隨迎角的變化曲線,抬頭為正力矩。2種翼型始終受到低頭力矩的作用,后退式Mini-TED翼型比原型翼型具有更大的低頭力矩,且力矩隨著迎角的增大而不斷增加。這表明,在襟翼上下方的壓力差造成的額外升力增大了翼型的低頭力矩;另外,由于翼型有效弦長(zhǎng)的增加,使翼型氣動(dòng)中心向后緣移動(dòng),氣動(dòng)中心到前緣的力臂加長(zhǎng),也是導(dǎo)致低頭力矩的增大的原因。2種翼型低頭力矩隨著升力增大而增加,因此2種翼型都具有較好的俯仰穩(wěn)定性。對(duì)比2條曲線可以看到,原型翼型在小于2°迎角時(shí)的低頭力矩隨迎角增加較快,而Mini-TED翼型的低頭力矩在該迎角范圍內(nèi)基本保持不變。大于2°迎角后,Mini-TED翼型低頭力矩曲線的斜率則高于原型翼型,即Mini-TED對(duì)翼型低頭力矩隨著迎角增加得比原型翼型快??傮w上,2種翼型的俯仰力矩系數(shù)曲線在絕大部分迎角情況下都表現(xiàn)出較好的線性,各迎角下Mini-TED翼型的低頭力矩均大于原型翼型。

    圖8 2種翼型俯仰力矩系數(shù)隨迎角變化曲線

    Fig.8 Variation of the pitch-moment coefficientCmwith angles of attack(AOA)

    2.3 速度分布測(cè)量結(jié)果

    壓力分布與流動(dòng)結(jié)構(gòu)密切相關(guān),因此課題組還進(jìn)行了不同迎角時(shí),翼型周?chē)乃俣确植糚IV測(cè)量。由于PIV測(cè)量數(shù)據(jù)量較大,關(guān)于各迎角PIV速度場(chǎng)、渦量場(chǎng)的詳細(xì)討論,請(qǐng)參閱文獻(xiàn)[18]。本文僅給出5°迎角下的速度場(chǎng)數(shù)據(jù),與圖4(b)進(jìn)行對(duì)比,討論Mini-TED翼型翼面速度分布的對(duì)壓力分布的影響。圖9為原型翼型與Mini-TED翼型翼面無(wú)量綱平均速度分布云圖和流線。Mini-TED翼型從前緣至后緣流速下降較為平緩,說(shuō)明翼面逆壓梯度小,流體有足夠的動(dòng)量克服逆壓造成的阻力,使上翼面流速普遍高于原型翼型。而原型翼型流動(dòng)承受較大的逆壓,從最大流速起向下游迅速衰減,尤其在后緣逆壓梯度的作用下,貼近翼面的區(qū)域流動(dòng)已接近停滯,流速遠(yuǎn)低于Mini-TED翼型。在原型翼型下翼面,流體保持貼體運(yùn)動(dòng),靜壓較低,而Mini-TED翼型流動(dòng)在后緣附近受到襟翼的阻礙,速度明顯低于原型翼型,整個(gè)翼面靜壓增加,并且在Mini-TED襟翼前部角區(qū)因速度降低為0而形成局部靜壓集中。因此,翼面的速度分布與圖4中壓力結(jié)果基本對(duì)應(yīng)。

    為探討Mini-TED翼型上翼面具有較高流速的原因,測(cè)量了后緣周?chē)牧鲃?dòng)結(jié)構(gòu)。圖10為5°迎角時(shí),NACA23012翼型后緣和Mini-TED襟翼周?chē)臅r(shí)均速度矢量,背景為無(wú)量綱化速度分布云圖,2圖采用相同的色系。圖中可以看出原型翼型后緣周?chē)魉倬陀谧杂蓙?lái)流速度,且上下翼面流動(dòng)在十分接近后緣的下游平滑匯合,沒(méi)有形成尾渦結(jié)構(gòu),且速度虧損區(qū)的范圍較小,在速度虧損區(qū)內(nèi)沒(méi)有觀察到倒流現(xiàn)象;而Mini-TED翼型后緣產(chǎn)生強(qiáng)烈的流動(dòng)分離,形成了旋轉(zhuǎn)方向相反的“對(duì)渦”結(jié)構(gòu),2渦之間有明顯的逆流存在,速度虧損區(qū)范圍顯著擴(kuò)大。由于渦旋結(jié)構(gòu)使靜壓降低,在“對(duì)渦”的誘導(dǎo)下流體加速運(yùn)動(dòng),可看出來(lái)自上翼面的流速已經(jīng)高于自由來(lái)流速度,因此導(dǎo)致Mini-TED上翼面壓力分布低于原型翼型。同時(shí)渦結(jié)構(gòu)消耗了流動(dòng)的動(dòng)能,引起阻力的增加;另外,還可觀察到Mini-TED翼型后緣上方邊界層厚度小于原型翼型,而邊界層外流速卻遠(yuǎn)高于原型翼型,因此受到更大的剪切應(yīng)力作用;上述2個(gè)原因使得Mini-TED翼型比原型翼型具有更大型阻。Mini-TED襟翼上方的低壓區(qū)還使翼型后緣附近上下翼面壓差增加,后緣獲得額外的升力,且力臂較長(zhǎng),翼型受到低頭力矩的作用,有利于翼型的俯仰穩(wěn)定性。5°迎角時(shí)2種翼型尾流結(jié)構(gòu)的差異普遍存在于翼型從0°迎角~失速迎角的范圍內(nèi),區(qū)別在于渦旋尺度和強(qiáng)度的改變。

    圖9 2種翼型周?chē)俣确植荚屏骶€圖

    圖10 2種翼型后緣流場(chǎng)結(jié)構(gòu)和速度分布云圖(迎角5°)

    Fig.10 Flow structures and contour plots of the velocity distribution (AOA 5°)

    3 結(jié) 論

    本文通過(guò)實(shí)驗(yàn)研究,討論了后退式Mini-TED對(duì)翼型和NACA23012翼型表面壓力分布、氣動(dòng)力以及速度分布的差異,并分析Mini-TED對(duì)翼型氣動(dòng)力產(chǎn)生的影響及原因。實(shí)驗(yàn)來(lái)流速度U∞=15m/s,以翼型弦長(zhǎng)為特征長(zhǎng)度的雷諾數(shù)Re=1.3×105。后退式Mini-TED上表面的吸力和下表面的壓力在各迎角時(shí),均比NACA23012翼型明顯增加,因此翼型升力的提高來(lái)自上下翼面壓力增長(zhǎng)的共同貢獻(xiàn);而原型翼型下翼面壓力分布幾乎不隨迎角增加發(fā)生變化,僅上翼面吸力增加引起升力的改變。Mini-TED翼型上翼面流速普遍高于原型翼型,可以延遲后緣流動(dòng)分離的發(fā)生,其翼面流速增加的主要原因一方面在于前駐點(diǎn)后移,導(dǎo)致上翼面前緣附近最大速度增加,另一方面由于后緣形成的渦結(jié)構(gòu)使翼面逆壓梯度降低。安裝Mini-TED后,翼型的壓差阻力也大于原型翼型;在中等迎角條件下(11°~16°),摩阻也有所增加,原因在于湍流邊界層覆蓋的翼面積隨迎角的增加而擴(kuò)大。然而,在相同升力系數(shù)下,后退式Mini-TED翼型仍具有更大的升阻比,且最大升阻比明顯高于原型翼型的最大升阻比。Mini-TED翼型的低頭力矩也有所增大,但是與原型翼型相比,翼型的俯仰穩(wěn)定性有所改善。

    本實(shí)驗(yàn)僅對(duì)特定形態(tài)的后退式Mini-TED翼型進(jìn)行一般氣動(dòng)特性和流場(chǎng)結(jié)構(gòu)的研究,未能對(duì)不同安裝狀態(tài)下,如襟翼長(zhǎng)度、安裝角等對(duì)翼型氣動(dòng)特性的影響進(jìn)行研究,因此有必要對(duì)Mini-TED形態(tài)的優(yōu)化特性開(kāi)展系列研究。

    [1] Richter K, Rosemann H. Steady aerodynamics of miniature trailing-edge devices in transonic flows[J]. Journal of Aircraft, 2012, 49(3): 898-910.

    [2] Richter K, Rosemann H. Experimental investigation of trailing-edge devices at transonic speeds[C]//Proceedings of Aerodynamics Conference 2001, Royal Aeronautical Society 2001, London, Paper No. 2621

    [3] Liebeck R H. Design of subsonic airfoils for high lift[J]. Journal of Aircraft, 1978, 15(9): 547-561.

    [4] Wang J J, Li Y C, Choi K S. Gurney flap-Lift enhancement, mechanisms and applications[J]. Progress in Aerospace Sciences, 2008, 44(1): 22-47.

    [5] Li Y C, Wang J J, Zhang P F. Influences of mounting angles and locations on the effects of Gurney flaps[J]. Journal of Aircraft, 2003, 40(3): 494-498.

    [6] Richter K, Rosemann H. Numerical nvestigation on the aerodynamic effect of Mini-TEDS on the AWIATOR aircraft at cruise conditions[C]//25th International Congress of the Aeronautical Sciences, 2006.

    [7] Gardner A D, Nitzsche J. Adaptive load redistribution using Mini-TEDS[C]//25th International Congress of the Aeronautical Sciences, 2006.

    [8] Hanse H. Application of adaptive trailing edge elements the AWIATOR project[R]. CEAS/KATnet Conference on Key Aerodynamic Technologies, Bremen, 2005.

    [9] Bechert D W, Meyer R, Hage W. Drag reduction of airfoils with miniflaps. can we learn from dragonflies?[J]. AIAA Paper, 2000, 2315: 2000.

    [10] 周華. Mini-TED改變翼型跨聲速性能的數(shù)值分析[J]. 航空學(xué)報(bào), 2009, 30(8): 1367-1373.

    Zhou Hua. Numerical analysis of Mini-TED on aerodynamic characteristics of airfoils in transonic flow[J]. Acta Aeronautica et Astronautica Sinica, 2009, 30(8): 1367-1373.

    [11] 褚胡冰, 張彬乾, 陳迎春, 等. 微型后緣裝置增升效率及幾何參數(shù)影響研究[J]. 航空學(xué)報(bào), 2012, 33(3): 381-389.

    Chu Hubin, Zhang Binqian, Chen Yingchun, et al. Investigation on Mini-TED efficiency and impact of its geometrical parameters[J]. Acta Aeronautica et Astronautica Sinica, 2012, 33(3): 381-389.

    [12] Zhou H. Rearward Mini-TED: its concept design and numerical validation[C]. Asia-Pacific International Symposium on Aerospace Technology, Jeju, Korea, 2012.

    [13] Bloy A W, Tsioumanis N, Mellor N T. Enhanced aerofoil performance using small trailing-edge flaps[J]. Journal of Aircraft, 1997, 34(4): 569-571.

    [14] Balaji R, Bramkamp F, Hesse M, et al. Effect of flap and slat riggings on 2-d high-lift aerodynamics[R]. AIAA-2006-19391, 2006.

    [15] 李亞臣, 王晉軍. NACA0012翼型Gurney襟翼增升特性及其機(jī)理實(shí)驗(yàn)研究[C]//第二屆全國(guó)低跨超聲速空氣動(dòng)力學(xué)交流會(huì)論文集, 2003, 2: 149-163.

    Li Yachen, Wang Jinjun. Experimental investigation on the lift enhancement and mechanism of NACA0012 airfoil with Gurney flap[C]//Proceeding of the 2nd National Low-speed, Transonic and Supersonic Aerodynamics Symposium, 2003, 2: 149-163.

    [16] Thompson B E, Lotz R D. Divergent trailing edge airfoil flow[R]. AIAA-1996-47040, 1996.

    [17] Anderson J D Jr. Fundamentals of aerodynamics[M]. 4th Edition. McGraw Hill.

    [18] 程巨擘, 代欽. 低雷諾數(shù)下Mini-TED對(duì)翼型流動(dòng)分離特性的影響PIV實(shí)驗(yàn)研究[J]. 實(shí)驗(yàn)流體力學(xué), 2014, 28(6): 20-26.

    Chen Jubei, Daichin. Experimental investigation on the influence of a Mini-TED on the flow separation of an airfoil at low Reynolds number using PIV[J]. Journal of Experiments in Fluid Mechanics, 2014, 28(6): 20-26.

    (編輯:楊 娟)

    The experimental investigation on the effect of rearward Mini-TED to the aerodynamic characteristics of an airfoil

    Xia Jun1, Daichin1,2,*

    (1. Key Laboratory of Mechanics in Energy Engineering, Shanghai Institute of Applied Mathematics and Mechanics, Shanghai University, Shanghai 200072, China; 2. Shanghai Key Laboratory of Mechanics in Energy Engineering, Shanghai 200072, China)

    The pressure distribution, aerodynamics and PIV velocity fields of a NACA23012 airfoil mounted with a rearward Mini-TED captured in a low speed wind tunnel tests are introduced and compared with experimental results of a NACA23012 prototype airfoil in this paper, in order to investigate the influence of the Mini-TED to the flow field and aerodynamics loads of the airfoil.The Reynolds number is 1.3×105based upon the chord length of the airfoil. The pressure distribution on the airfoil is measured using the pressure tubes and pressure transducers, and then the results are integrated to obtain the lift force and pressure drag acting on the airfoil. The total drag is measured using the comb of stagnation pressure based on the momentum theorem. The velocity fields around the airfoil are captured by a PIV system. The rearward Mini-TED induces changes to the flow velocity around the airfoil and the wake flow structure, which leads to the variation of the pressure distribution on the suction and pressure surfaces and the increase of both the lift and pressure drag. The position of the front stagnation point of the airfoil with a Mini-TED shifts downstream slightly compared with that of the NACA23012 prototype airfoil, which consequently leads to the increase of the flow velocity on the suction surface and the suppression of the flow separation near the trailing edge.

    rearward Mini-TED;aerodynamics of airfoil;low Reynolds number;surface pressure distribution;PIV measurement

    1672-9897(2015)05-0001-08

    10.11729/syltlx20140148

    2014-12-18;

    2015-01-18

    國(guó)家自然科學(xué)基金項(xiàng)目(11472169和11072142)

    XiaJ,Daichin.TheexperimentalinvestigationontheeffectofrearwardMini-TEDtotheaerodynamiccharacteristicsofanairfoil.JournalofExperimentsinFluidMechanics, 2015, 29(5): 1-7,25. 夏 駿, 代 欽. 后退式微型后緣裝置對(duì)翼型氣動(dòng)特性影響的實(shí)驗(yàn)研究. 實(shí)驗(yàn)流體力學(xué), 2015, 29(5): 1-7,25.

    V211.41

    A

    夏 駿(1990-),男,上海人,碩士研究生。研究方向:實(shí)驗(yàn)流體力學(xué)。通信地址:上海大學(xué)上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所(200072)。E-mail:dult1998615@sina.com

    *通信作者 E-mail: daichin@staff.shu.edu.cn

    猜你喜歡
    翼面襟翼后緣
    基于拔銷(xiāo)器鎖定的飛行器氣動(dòng)控制面解鎖控制方法
    民用飛機(jī)襟翼交聯(lián)機(jī)構(gòu)吸能仿真技術(shù)研究
    固定翼二維彈道修正引信升力翼面位置的影響
    基于Aerobook平臺(tái)的復(fù)合材料翼面結(jié)構(gòu)設(shè)計(jì)流程
    機(jī)翼后緣連續(xù)變彎度對(duì)客機(jī)氣動(dòng)特性影響
    某型機(jī)襟翼系統(tǒng)通電檢查故障分析
    柔性后緣可變形機(jī)翼氣動(dòng)特性分析
    737NG飛機(jī)的后緣襟翼指示故障
    科技尚品(2016年6期)2016-07-06 08:54:13
    升力式再入飛行器體襟翼姿態(tài)控制方法
    TNF-α和PGP9.5在椎體后緣離斷癥軟骨終板的表達(dá)及意義
    国产三级黄色录像| 丝袜在线中文字幕| 韩国av一区二区三区四区| 日韩大尺度精品在线看网址| 好男人在线观看高清免费视频 | 午夜成年电影在线免费观看| 久久久久久亚洲精品国产蜜桃av| av片东京热男人的天堂| 女警被强在线播放| 一本一本综合久久| 啪啪无遮挡十八禁网站| 欧美zozozo另类| 一级片免费观看大全| 大型av网站在线播放| 啪啪无遮挡十八禁网站| 亚洲av熟女| 亚洲人成77777在线视频| 亚洲av日韩精品久久久久久密| 亚洲av成人av| 精品乱码久久久久久99久播| 国产欧美日韩一区二区三| 精品人妻1区二区| √禁漫天堂资源中文www| 色综合欧美亚洲国产小说| 特大巨黑吊av在线直播 | 亚洲成人久久性| 久久精品91蜜桃| 精品欧美国产一区二区三| 丝袜美腿诱惑在线| 国内毛片毛片毛片毛片毛片| 国产精品自产拍在线观看55亚洲| 国产三级黄色录像| 亚洲无线在线观看| 国产精品精品国产色婷婷| 精品国产乱码久久久久久男人| 亚洲狠狠婷婷综合久久图片| 嫩草影院精品99| or卡值多少钱| 免费在线观看完整版高清| 午夜福利一区二区在线看| 日本 av在线| 草草在线视频免费看| 国产精品久久久久久精品电影 | 国产国语露脸激情在线看| 一进一出抽搐动态| x7x7x7水蜜桃| 狂野欧美激情性xxxx| 人妻久久中文字幕网| 国产亚洲欧美精品永久| 国产成人啪精品午夜网站| 一卡2卡三卡四卡精品乱码亚洲| 国产激情久久老熟女| 老司机福利观看| 熟女少妇亚洲综合色aaa.| 午夜久久久久精精品| 最近最新中文字幕大全电影3 | 两性午夜刺激爽爽歪歪视频在线观看 | 中文字幕另类日韩欧美亚洲嫩草| 黑丝袜美女国产一区| 免费女性裸体啪啪无遮挡网站| av视频在线观看入口| 欧美黑人巨大hd| 亚洲精品久久国产高清桃花| a在线观看视频网站| www.www免费av| 又大又爽又粗| 男人的好看免费观看在线视频 | 99在线视频只有这里精品首页| 亚洲成人久久爱视频| e午夜精品久久久久久久| 一进一出抽搐动态| 亚洲人成伊人成综合网2020| e午夜精品久久久久久久| 午夜福利成人在线免费观看| 中文字幕最新亚洲高清| 亚洲 欧美一区二区三区| 高潮久久久久久久久久久不卡| 真人做人爱边吃奶动态| 国产亚洲精品综合一区在线观看 | 黄色毛片三级朝国网站| av电影中文网址| 亚洲人成网站在线播放欧美日韩| 亚洲精品中文字幕在线视频| 少妇粗大呻吟视频| 精品人妻1区二区| 草草在线视频免费看| АⅤ资源中文在线天堂| 欧美日韩中文字幕国产精品一区二区三区| 99国产综合亚洲精品| av超薄肉色丝袜交足视频| 精品久久久久久,| 国产精品久久久久久精品电影 | 欧美久久黑人一区二区| 日日爽夜夜爽网站| 国产成人精品久久二区二区免费| 不卡av一区二区三区| 中文字幕精品免费在线观看视频| 一级a爱视频在线免费观看| 妹子高潮喷水视频| 午夜视频精品福利| 亚洲国产毛片av蜜桃av| 国产一卡二卡三卡精品| 欧美日韩中文字幕国产精品一区二区三区| 亚洲性夜色夜夜综合| 亚洲精华国产精华精| 久99久视频精品免费| 校园春色视频在线观看| 久久人妻av系列| 中文字幕高清在线视频| 国产精品98久久久久久宅男小说| 99久久久亚洲精品蜜臀av| 好男人电影高清在线观看| 久久狼人影院| 亚洲国产欧洲综合997久久, | 精品久久久久久久久久免费视频| 一级毛片高清免费大全| 国产伦在线观看视频一区| 欧美成人午夜精品| 国产精品影院久久| 久热这里只有精品99| 日韩三级视频一区二区三区| 18禁美女被吸乳视频| 亚洲精华国产精华精| a在线观看视频网站| 操出白浆在线播放| 99久久精品国产亚洲精品| 欧美黄色淫秽网站| 久久久久精品国产欧美久久久| 亚洲精品av麻豆狂野| 久久精品亚洲精品国产色婷小说| 成人国产综合亚洲| 日本黄色视频三级网站网址| 亚洲av成人不卡在线观看播放网| www日本黄色视频网| 中文字幕人成人乱码亚洲影| 亚洲精品色激情综合| 免费看美女性在线毛片视频| 桃红色精品国产亚洲av| 在线观看日韩欧美| 女人被狂操c到高潮| 在线播放国产精品三级| 免费观看精品视频网站| 不卡av一区二区三区| 成人特级黄色片久久久久久久| 成人国产一区最新在线观看| 国产精品久久久久久精品电影 | www.精华液| 亚洲中文日韩欧美视频| 熟女电影av网| 国产亚洲精品第一综合不卡| 久久婷婷人人爽人人干人人爱| 欧美黑人欧美精品刺激| 久久中文字幕人妻熟女| 久久精品91无色码中文字幕| 国产极品粉嫩免费观看在线| 欧美日韩亚洲国产一区二区在线观看| 人人澡人人妻人| 中文字幕精品免费在线观看视频| 亚洲 欧美 日韩 在线 免费| 国产高清视频在线播放一区| 18禁黄网站禁片免费观看直播| 亚洲av成人av| 两个人视频免费观看高清| 国产麻豆成人av免费视频| 中出人妻视频一区二区| 九色国产91popny在线| 好看av亚洲va欧美ⅴa在| 手机成人av网站| 麻豆久久精品国产亚洲av| 亚洲第一av免费看| 欧美激情极品国产一区二区三区| 狠狠狠狠99中文字幕| tocl精华| 国产精品亚洲av一区麻豆| 50天的宝宝边吃奶边哭怎么回事| 亚洲自偷自拍图片 自拍| 欧美亚洲日本最大视频资源| 日韩欧美一区视频在线观看| 亚洲精品美女久久av网站| 男女视频在线观看网站免费 | 这个男人来自地球电影免费观看| 久久久久久亚洲精品国产蜜桃av| 级片在线观看| 香蕉久久夜色| 热re99久久国产66热| 色综合婷婷激情| 国内精品久久久久精免费| 一本大道久久a久久精品| 成人av一区二区三区在线看| 在线永久观看黄色视频| 嫩草影院精品99| 真人一进一出gif抽搐免费| 黄色片一级片一级黄色片| 两人在一起打扑克的视频| 久久久久久久久免费视频了| 757午夜福利合集在线观看| 97碰自拍视频| 中文字幕最新亚洲高清| 在线看三级毛片| 精品一区二区三区av网在线观看| 欧美不卡视频在线免费观看 | 午夜视频精品福利| 18禁国产床啪视频网站| av视频在线观看入口| 很黄的视频免费| 成人18禁在线播放| 亚洲va日本ⅴa欧美va伊人久久| 少妇粗大呻吟视频| 国产精品免费一区二区三区在线| 久久国产精品影院| 久久久久国产精品人妻aⅴ院| 女人高潮潮喷娇喘18禁视频| 亚洲国产欧洲综合997久久, | 亚洲一区中文字幕在线| 熟女少妇亚洲综合色aaa.| 黄色片一级片一级黄色片| 啦啦啦韩国在线观看视频| 国语自产精品视频在线第100页| 日韩成人在线观看一区二区三区| 女性生殖器流出的白浆| 国产精品美女特级片免费视频播放器 | 亚洲va日本ⅴa欧美va伊人久久| 2021天堂中文幕一二区在线观 | 91老司机精品| www日本在线高清视频| www.自偷自拍.com| 亚洲久久久国产精品| 一个人观看的视频www高清免费观看 | 在线观看免费午夜福利视频| 黄色丝袜av网址大全| 欧美乱色亚洲激情| 每晚都被弄得嗷嗷叫到高潮| 亚洲av美国av| 麻豆国产av国片精品| 日日爽夜夜爽网站| 亚洲中文字幕日韩| 在线永久观看黄色视频| 久久久久久久精品吃奶| 首页视频小说图片口味搜索| 一级毛片精品| 身体一侧抽搐| 国产免费av片在线观看野外av| 亚洲黑人精品在线| 亚洲第一欧美日韩一区二区三区| 免费一级毛片在线播放高清视频| 身体一侧抽搐| 免费在线观看完整版高清| 久久性视频一级片| 国产久久久一区二区三区| 国产成+人综合+亚洲专区| 此物有八面人人有两片| 在线天堂中文资源库| 成人手机av| 久久热在线av| 侵犯人妻中文字幕一二三四区| 日韩 欧美 亚洲 中文字幕| 亚洲狠狠婷婷综合久久图片| 18禁黄网站禁片午夜丰满| 中文字幕最新亚洲高清| 国产精品99久久99久久久不卡| 国产不卡一卡二| 午夜亚洲福利在线播放| 国产精品亚洲av一区麻豆| 黑人操中国人逼视频| 一区二区三区高清视频在线| 久久中文看片网| 俺也久久电影网| 国产aⅴ精品一区二区三区波| 一级a爱视频在线免费观看| 久久婷婷人人爽人人干人人爱| 日本 欧美在线| 国产v大片淫在线免费观看| 国产1区2区3区精品| 日韩精品青青久久久久久| 91成年电影在线观看| 美女 人体艺术 gogo| 黄色毛片三级朝国网站| 午夜两性在线视频| 老熟妇仑乱视频hdxx| 中文亚洲av片在线观看爽| netflix在线观看网站| 国产不卡一卡二| 一级毛片精品| 国产成+人综合+亚洲专区| 18禁观看日本| 在线天堂中文资源库| 美女 人体艺术 gogo| 超碰成人久久| 在线观看一区二区三区| 1024视频免费在线观看| 久久国产精品人妻蜜桃| 老司机靠b影院| 久久久久久国产a免费观看| 中文字幕精品亚洲无线码一区 | 亚洲 国产 在线| 亚洲狠狠婷婷综合久久图片| 少妇的丰满在线观看| 女性被躁到高潮视频| 亚洲片人在线观看| www.www免费av| 大型av网站在线播放| 国产亚洲欧美98| 成人永久免费在线观看视频| 日韩成人在线观看一区二区三区| 国产黄色小视频在线观看| 女人被狂操c到高潮| 国产成+人综合+亚洲专区| 97人妻精品一区二区三区麻豆 | 这个男人来自地球电影免费观看| 观看免费一级毛片| 少妇粗大呻吟视频| 亚洲无线在线观看| 亚洲一区高清亚洲精品| 欧美性长视频在线观看| 久久精品aⅴ一区二区三区四区| 欧美黑人巨大hd| 我的亚洲天堂| 国产精品久久久久久亚洲av鲁大| 极品教师在线免费播放| 最近最新中文字幕大全电影3 | 久久国产精品男人的天堂亚洲| 国产成人精品久久二区二区91| 法律面前人人平等表现在哪些方面| 精品少妇一区二区三区视频日本电影| 精品乱码久久久久久99久播| 一级作爱视频免费观看| 一级毛片高清免费大全| 成人国产一区最新在线观看| 久久中文字幕人妻熟女| 可以免费在线观看a视频的电影网站| 欧美国产精品va在线观看不卡| 国产主播在线观看一区二区| 国产精品 国内视频| 成人18禁高潮啪啪吃奶动态图| 成人av一区二区三区在线看| 精品久久久久久久毛片微露脸| 国产成人欧美在线观看| 777久久人妻少妇嫩草av网站| 黄片小视频在线播放| 操出白浆在线播放| netflix在线观看网站| 男人舔女人下体高潮全视频| 久久久久国内视频| 脱女人内裤的视频| 午夜福利在线在线| 欧美日韩乱码在线| 久久久久精品国产欧美久久久| 欧美日韩黄片免| 亚洲av电影在线进入| 精品久久蜜臀av无| 亚洲成人免费电影在线观看| 欧美日韩中文字幕国产精品一区二区三区| 免费在线观看亚洲国产| 亚洲最大成人中文| 少妇被粗大的猛进出69影院| 久久久国产成人免费| 99国产精品一区二区三区| 国产免费av片在线观看野外av| 制服诱惑二区| netflix在线观看网站| 一区福利在线观看| 亚洲成av人片免费观看| 亚洲av成人一区二区三| 国产在线观看jvid| 国产成人影院久久av| 丁香六月欧美| 午夜福利18| 久久久久国内视频| 亚洲 欧美一区二区三区| 亚洲国产精品sss在线观看| 久久精品成人免费网站| 看免费av毛片| 午夜影院日韩av| a在线观看视频网站| 亚洲国产日韩欧美精品在线观看 | 哪里可以看免费的av片| 久久精品国产99精品国产亚洲性色| 一级黄色大片毛片| 女性生殖器流出的白浆| 国产亚洲欧美精品永久| 999久久久精品免费观看国产| 国内毛片毛片毛片毛片毛片| 成人18禁在线播放| 亚洲成人久久性| а√天堂www在线а√下载| 正在播放国产对白刺激| 亚洲专区中文字幕在线| 国产精品av久久久久免费| www.精华液| 国产爱豆传媒在线观看 | 久久天躁狠狠躁夜夜2o2o| 精品高清国产在线一区| 日本 av在线| 97超级碰碰碰精品色视频在线观看| 国产精品98久久久久久宅男小说| 啦啦啦 在线观看视频| 在线av久久热| 色哟哟哟哟哟哟| 国产黄色小视频在线观看| 精品国产乱子伦一区二区三区| 精华霜和精华液先用哪个| 99精品久久久久人妻精品| 激情在线观看视频在线高清| 欧美激情极品国产一区二区三区| 色综合欧美亚洲国产小说| 日韩欧美 国产精品| 日本 欧美在线| 国产伦在线观看视频一区| 亚洲久久久国产精品| 成熟少妇高潮喷水视频| 欧美一区二区精品小视频在线| 51午夜福利影视在线观看| 亚洲在线自拍视频| 99国产极品粉嫩在线观看| 久久久久九九精品影院| 久久久久久久久中文| 99热只有精品国产| 中出人妻视频一区二区| 一区二区三区激情视频| 午夜亚洲福利在线播放| 夜夜躁狠狠躁天天躁| 草草在线视频免费看| 午夜老司机福利片| 免费无遮挡裸体视频| 亚洲精品粉嫩美女一区| 亚洲国产精品久久男人天堂| 日本五十路高清| 在线观看午夜福利视频| 美国免费a级毛片| АⅤ资源中文在线天堂| 国产久久久一区二区三区| 成人国产一区最新在线观看| 日本 欧美在线| 亚洲熟妇中文字幕五十中出| 国内精品久久久久久久电影| 一级黄色大片毛片| 日韩视频一区二区在线观看| 亚洲欧洲精品一区二区精品久久久| 久久精品亚洲精品国产色婷小说| 国产精品一区二区三区四区久久 | 久久久久久久久中文| 欧美+亚洲+日韩+国产| 91成年电影在线观看| 制服人妻中文乱码| 巨乳人妻的诱惑在线观看| bbb黄色大片| 亚洲精品一卡2卡三卡4卡5卡| 日韩视频一区二区在线观看| 精品久久久久久久久久免费视频| 美女免费视频网站| 亚洲人成电影免费在线| 亚洲欧美精品综合一区二区三区| 黑丝袜美女国产一区| 男女那种视频在线观看| 久久狼人影院| e午夜精品久久久久久久| 日韩欧美 国产精品| 欧美成人一区二区免费高清观看 | 男女视频在线观看网站免费 | 女性被躁到高潮视频| 国产精品爽爽va在线观看网站 | 老汉色∧v一级毛片| 母亲3免费完整高清在线观看| 精品卡一卡二卡四卡免费| 精品乱码久久久久久99久播| 特大巨黑吊av在线直播 | xxxwww97欧美| 一区二区三区高清视频在线| 在线观看午夜福利视频| 欧美性长视频在线观看| 亚洲狠狠婷婷综合久久图片| 黄片播放在线免费| 丰满的人妻完整版| 精品午夜福利视频在线观看一区| 欧美又色又爽又黄视频| 国产色视频综合| 色老头精品视频在线观看| 丝袜美腿诱惑在线| 国产1区2区3区精品| 亚洲成av片中文字幕在线观看| 亚洲avbb在线观看| 国产真人三级小视频在线观看| 中文字幕人妻熟女乱码| 日本成人三级电影网站| 日韩欧美免费精品| 精品卡一卡二卡四卡免费| 国产一区二区激情短视频| 最新美女视频免费是黄的| 男人舔女人下体高潮全视频| 国产97色在线日韩免费| 老司机在亚洲福利影院| 亚洲精品色激情综合| 欧美国产日韩亚洲一区| 妹子高潮喷水视频| 99久久无色码亚洲精品果冻| 日韩国内少妇激情av| 欧美性猛交╳xxx乱大交人| 久久 成人 亚洲| 国产黄a三级三级三级人| 淫妇啪啪啪对白视频| 国产精品久久久久久精品电影 | 亚洲成人久久性| 久久午夜亚洲精品久久| 中文亚洲av片在线观看爽| a级毛片在线看网站| 婷婷亚洲欧美| 满18在线观看网站| netflix在线观看网站| www.999成人在线观看| 久久国产精品人妻蜜桃| 久久狼人影院| 欧美一级毛片孕妇| 村上凉子中文字幕在线| 国产日本99.免费观看| 久久欧美精品欧美久久欧美| 最近在线观看免费完整版| 亚洲国产高清在线一区二区三 | 亚洲国产精品久久男人天堂| 日韩精品免费视频一区二区三区| 成熟少妇高潮喷水视频| 精华霜和精华液先用哪个| 桃色一区二区三区在线观看| 黄色 视频免费看| 啪啪无遮挡十八禁网站| 一区二区三区国产精品乱码| 日韩三级视频一区二区三区| 午夜久久久久精精品| 久久中文看片网| 亚洲免费av在线视频| 精品电影一区二区在线| 啪啪无遮挡十八禁网站| 99国产精品一区二区三区| 欧美性猛交╳xxx乱大交人| 色播在线永久视频| 国产精品爽爽va在线观看网站 | 熟妇人妻久久中文字幕3abv| 亚洲av日韩精品久久久久久密| 久久国产亚洲av麻豆专区| 在线天堂中文资源库| 黄色毛片三级朝国网站| 久久久国产精品麻豆| 99精品在免费线老司机午夜| 19禁男女啪啪无遮挡网站| 日日摸夜夜添夜夜添小说| 国产精品九九99| 日日爽夜夜爽网站| 欧美激情极品国产一区二区三区| 怎么达到女性高潮| 在线天堂中文资源库| 国产精品电影一区二区三区| 国产精品九九99| 看片在线看免费视频| 国产欧美日韩精品亚洲av| 精品国产乱子伦一区二区三区| 亚洲国产欧洲综合997久久, | 国产成年人精品一区二区| 一级毛片女人18水好多| avwww免费| 最新在线观看一区二区三区| 精品久久久久久,| 欧美黄色淫秽网站| 成人国产一区最新在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲中文字幕日韩| 欧美色视频一区免费| 亚洲专区中文字幕在线| 精品无人区乱码1区二区| 久久人人精品亚洲av| 99热6这里只有精品| 亚洲成人精品中文字幕电影| 久久草成人影院| 久久精品aⅴ一区二区三区四区| 欧美成人免费av一区二区三区| avwww免费| 国产成人啪精品午夜网站| 国产又色又爽无遮挡免费看| 亚洲精品国产一区二区精华液| 一区二区日韩欧美中文字幕| 丰满人妻熟妇乱又伦精品不卡| 亚洲激情在线av| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品美女久久av网站| 热re99久久国产66热| 日韩视频一区二区在线观看| 不卡av一区二区三区| 最近最新中文字幕大全免费视频| 国产精品1区2区在线观看.| 麻豆av在线久日| 久久久久久免费高清国产稀缺| 日韩有码中文字幕| 最好的美女福利视频网| 欧美一级毛片孕妇| 黄色女人牲交| 成人18禁高潮啪啪吃奶动态图| 免费高清视频大片| 国产亚洲精品综合一区在线观看 | 香蕉丝袜av| 欧美精品啪啪一区二区三区| 又大又爽又粗| 12—13女人毛片做爰片一| 非洲黑人性xxxx精品又粗又长| 午夜免费成人在线视频| 精品久久蜜臀av无| 亚洲精品久久成人aⅴ小说| 成熟少妇高潮喷水视频| 国产精品久久久人人做人人爽| 精品国产乱子伦一区二区三区| 午夜福利免费观看在线| 搞女人的毛片| 麻豆av在线久日| 老汉色∧v一级毛片| 97超级碰碰碰精品色视频在线观看| www日本黄色视频网| 最近最新中文字幕大全电影3 | 别揉我奶头~嗯~啊~动态视频|