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

    TBCC發(fā)動(dòng)機(jī)進(jìn)氣機(jī)匣流道評(píng)估與優(yōu)化

    2015-01-06 02:47:36王俊余華蔚蔡留成楊偉
    燃?xì)鉁u輪試驗(yàn)與研究 2015年1期
    關(guān)鍵詞:支板恢復(fù)系數(shù)流道

    王俊,余華蔚,蔡留成,楊偉

    (中國(guó)燃?xì)鉁u輪研究院,四川成都610500)

    TBCC發(fā)動(dòng)機(jī)進(jìn)氣機(jī)匣流道評(píng)估與優(yōu)化

    王俊,余華蔚,蔡留成,楊偉

    (中國(guó)燃?xì)鉁u輪研究院,四川成都610500)

    對(duì)渦輪基組合循環(huán)(TBCC)發(fā)動(dòng)機(jī)進(jìn)氣機(jī)匣進(jìn)行全三維流場(chǎng)計(jì)算,評(píng)估進(jìn)氣機(jī)匣流道的流場(chǎng)和性能,發(fā)現(xiàn)進(jìn)氣機(jī)匣流道中氣流損失過(guò)大,總壓恢復(fù)系數(shù)不滿足指標(biāo)要求,需進(jìn)行優(yōu)化。然后根據(jù)流場(chǎng)分析結(jié)果和流道截面積變化,使用三次樣條插值方法對(duì)進(jìn)氣機(jī)匣流道進(jìn)行優(yōu)化,使流道截面積變化更為平緩。對(duì)修改設(shè)計(jì)后的流道再進(jìn)行全三維流場(chǎng)計(jì)算分析,反復(fù)迭代使得最終的進(jìn)氣機(jī)匣流道和支板的總壓恢復(fù)系數(shù)大幅提高,并滿足指標(biāo)要求。

    渦輪基組合循環(huán)發(fā)動(dòng)機(jī);進(jìn)氣機(jī)匣;三次樣條;優(yōu)化設(shè)計(jì);三維計(jì)算

    1 引言

    渦輪基組合循環(huán)(TBCC)發(fā)動(dòng)機(jī)作為一種復(fù)雜度相當(dāng)高的新型吸氣式推進(jìn)系統(tǒng),進(jìn)氣機(jī)匣是其重要組成部分[1]。設(shè)計(jì)進(jìn)氣機(jī)匣時(shí),流過(guò)其中的氣流應(yīng)滿足其后壓氣機(jī)的使用需要,并使氣流以最小的總壓損失和氣動(dòng)阻力流過(guò)進(jìn)氣機(jī)匣的彎曲流道[2]。另外,進(jìn)排氣機(jī)匣設(shè)計(jì)對(duì)壓氣機(jī)整機(jī)效率的影響不能忽略。進(jìn)排氣機(jī)匣設(shè)計(jì)不合理,有時(shí)會(huì)造成總壓損失過(guò)大,需增加級(jí)數(shù)才能保證壓氣機(jī)的整機(jī)性能。而合理的壓氣機(jī)進(jìn)氣機(jī)匣設(shè)計(jì),能保證葉片入口處的流場(chǎng)均勻[3]。為使TBCC發(fā)動(dòng)機(jī)在寬廣的飛行馬赫數(shù)范圍都能有效工作并具有高的氣動(dòng)性能,必須保證飛行器進(jìn)氣系統(tǒng)與發(fā)動(dòng)機(jī)氣動(dòng)性能和結(jié)構(gòu)布局匹配[4]。

    國(guó)內(nèi)外已對(duì)各類飛行器進(jìn)氣機(jī)匣流道的優(yōu)化進(jìn)行了很多研究[1-5],但公開(kāi)文獻(xiàn)中還沒(méi)有采用三次樣條插值方法對(duì)流道進(jìn)行優(yōu)化的先例。本文以TBCC組合發(fā)動(dòng)機(jī)進(jìn)氣機(jī)匣的設(shè)計(jì)工況為研究對(duì)象,對(duì)進(jìn)氣機(jī)匣中的流場(chǎng)進(jìn)行三維計(jì)算與特性分析,并根據(jù)計(jì)算分析結(jié)果,采用三次樣條插值方法對(duì)進(jìn)氣機(jī)匣進(jìn)行優(yōu)化設(shè)計(jì),改善其流場(chǎng),降低總壓損失。

    2 優(yōu)化思路及方法

    2.1 優(yōu)化思路

    根據(jù)進(jìn)氣機(jī)匣流道截面積變化來(lái)優(yōu)化流道形狀。這是由于流道中的損失主要由壓力梯度和曲率共同決定[5],根據(jù)流量連續(xù)方程(式(1))可知:當(dāng)流量不變、流體不可壓時(shí),流體速度v與流道截面積A成反比,若流道截面積連續(xù)變化,則流體速度也將連續(xù)、光滑變化,所導(dǎo)致的壓力損失也最小。當(dāng)流過(guò)進(jìn)氣機(jī)匣的流量為定值且流速較低時(shí),進(jìn)口氣流馬赫數(shù)較小,因此可認(rèn)為機(jī)匣內(nèi)氣流為不可壓流,密度基本不變。這時(shí),若希望氣流流過(guò)流道的壓力損失最小,除要求流道壁面光滑轉(zhuǎn)接、保證流道幾何曲率連續(xù)變化外,還應(yīng)要求流道截面積平滑變化,以保證流體速度變化連續(xù)、順暢,從而使壓力梯度光滑變化。

    2.2 優(yōu)化方法

    流道優(yōu)化過(guò)程中,為保證流道曲率光滑,采用三次樣條插值函數(shù)[6]對(duì)流道數(shù)據(jù)點(diǎn)進(jìn)行插值。

    三次樣條插值函數(shù)的構(gòu)造方法[7]是:假設(shè)給定n+1個(gè)點(diǎn)

    三次樣條插值函數(shù)S(x)在每個(gè)子區(qū)間[xk,xk+1]上是三次多項(xiàng)式,故S″(x)是線性函數(shù)。令mk=S″(xk) (k=0,1,2,…,n),設(shè)x∈[xk,xk+1],則過(guò)兩點(diǎn)(xk,mk)與(xk+1,mk+1)的線性函數(shù)可表示為:

    式中:hk=xk+1-xk。

    對(duì)上式兩端在區(qū)間[xk,xk+1]上連續(xù)兩次求積分,即可得到:

    因此,確定了mk,就可確定各子區(qū)間上的S(x)。由于S(x)在[a,b]上的一階導(dǎo)數(shù)連續(xù),即S′(x-0)=S′(x+0),則由上式可推出:

    這樣共得到n-1個(gè)方程,但由于有n+1個(gè)變量mk,故還需兩個(gè)方程才能求得唯一解。為此,本文選擇自然邊界條件作為附加條件,即樣條函數(shù)S(x)在插值區(qū)間[a,b]兩端點(diǎn)處的二階導(dǎo)數(shù)為零,產(chǎn)生S″(x0)=0、S″(xn)=0兩個(gè)方程。將以上n+1個(gè)方程聯(lián)立求得mk之后,即可確定S(x)。

    綜上所述,三次樣條插值計(jì)算步驟為:

    (1)確定邊界條件;

    (2)根據(jù)所確定的邊界條件計(jì)算各值,形成方程組;

    (3)解方程組,求得mk(k=0,1,2,…,n);

    (4)將所求mk值代入S(x)的表達(dá)式,求得區(qū)間[a,b]內(nèi)任意一點(diǎn)的近似值。

    3 進(jìn)氣機(jī)匣流場(chǎng)計(jì)算分析

    采用全三維計(jì)算軟件NUMECA的Spalart-All?maras(S-A)模型,進(jìn)行全三維粘性流場(chǎng)分析。為便于分析,在計(jì)算中設(shè)置了幾個(gè)關(guān)鍵分析截面,如圖1所示。進(jìn)氣機(jī)匣段(0-2段)流道的總壓恢復(fù)系數(shù)指標(biāo)要求不低于0.985 0。

    圖1 進(jìn)氣機(jī)匣各關(guān)鍵分析截面示意圖Fig.1 Sketch of each analytical section

    圖2~圖4分別為流道各分析截面總壓恢復(fù)系數(shù)、氣流角、相對(duì)馬赫數(shù)沿徑向的分布曲線,5%、50%和95%葉高擬S1流面馬赫數(shù)等值線分布見(jiàn)圖5,近葉面擬S2流面馬赫數(shù)等值線分布見(jiàn)圖6。

    圖2 各段流道總壓恢復(fù)系數(shù)的徑向分布Fig.2 The radial distribution of the pressure recovery coefficient

    從圖2中喇叭口到支板前(0-1段)的總壓恢復(fù)系數(shù)徑向分布可以看出,氣流從入口到支板前,根部和尖部能量損失較大。氣流在流過(guò)支板的過(guò)程中,也是根部和尖部的損失較中部截面大。其原因可能是,越靠近流道壁面,氣流受壁面的影響越大。圖3中,氣流從喇叭口到流過(guò)支板,氣流角都沒(méi)有變化。圖4中,在支板前(1-1截面)流道中間位置,氣流馬赫數(shù)接近0.7,而根部和尖部約為0.4,這就使得支板前流場(chǎng)不均勻,加大了氣流經(jīng)過(guò)支板時(shí)的壓力損失,從而使得整個(gè)進(jìn)氣段流道(0-2段)的氣流損失偏大。若對(duì)0-2段流道進(jìn)行優(yōu)化,降低氣流馬赫數(shù),使流道中的氣流更均勻,就能減少損失,提高總壓恢復(fù)系數(shù),這樣支板前進(jìn)氣馬赫數(shù)降低也更均勻,氣流經(jīng)過(guò)支板的壓力損失就會(huì)降低。

    圖3 各截面氣流角的徑向分布Fig.3 The radial distribution of relative flow angle

    圖4 各截面馬赫數(shù)的徑向分布Fig.4 The radial distribution of relative Mach number

    從圖4、圖5可以看出,支板表面中部截面位置的馬赫數(shù)最高,根部和尖部的馬赫數(shù)較低。從圖6可以看出,當(dāng)氣流進(jìn)入流道流過(guò)進(jìn)氣罩時(shí),流道截面積突然變小,造成氣流速度迅速上升,氣流馬赫數(shù)最高點(diǎn)出現(xiàn)在截面積最小位置;氣流經(jīng)喉道到達(dá)支板前,又由于流道截面積變大而流速迅速降低,損失增加;而進(jìn)氣罩正前方氣流由于進(jìn)氣罩的阻擋,速度迅速降低。氣流流過(guò)支板時(shí),也因流道截面減小而加速。

    流道優(yōu)化前的三維計(jì)算結(jié)果如表1所示,可見(jiàn)支板(1-2段)的總壓恢復(fù)系數(shù)為0.981 5,沒(méi)能達(dá)到指標(biāo)要求。所以應(yīng)對(duì)進(jìn)氣帽罩段流道進(jìn)行優(yōu)化設(shè)計(jì),提高進(jìn)氣機(jī)匣流道和支板的總壓恢復(fù)系數(shù)。

    表1 進(jìn)氣機(jī)匣流道優(yōu)化前后三維計(jì)算結(jié)果對(duì)比Table 1 The contrast of the 3D calculation results

    圖5 流道優(yōu)化前擬S1流面馬赫數(shù)分布Fag.5 Relative Mach number distribution at S1 stream surface before optimization

    圖6 流道優(yōu)化前近葉面擬S2流面馬赫數(shù)分布Fag.6 Relative Mach number distribution at S2 stream surface before optimization

    4 進(jìn)氣機(jī)匣流道優(yōu)化設(shè)計(jì)

    4.1 流道優(yōu)化

    基于2.1節(jié)的優(yōu)化思路,通過(guò)調(diào)節(jié)帽罩前端流道的控制點(diǎn)對(duì)流道進(jìn)行優(yōu)化,使帽罩段的流道截面積擴(kuò)大并平緩變化,同時(shí)擴(kuò)大支板最大厚度處的流道面積;然后再進(jìn)行三維計(jì)算,根據(jù)計(jì)算結(jié)果對(duì)流道進(jìn)行修改。經(jīng)過(guò)近20輪優(yōu)化迭代后,進(jìn)氣機(jī)匣帽罩部分流道優(yōu)化前后對(duì)比如圖7所示(圖中沒(méi)有顯示的部分未進(jìn)行修改)??梢?jiàn),流道內(nèi)表面只在帽罩前端有微小調(diào)整,流道外表面部分外擴(kuò)。圖8示出了優(yōu)化前后流道截面積變化曲線對(duì)比,可見(jiàn)優(yōu)化后流道截面積在喉道部位明顯擴(kuò)大,且變化也更加緩和。

    圖7 進(jìn)氣機(jī)匣流道優(yōu)化前后對(duì)比Fig.7 The contrast of flow passage before and after optimization

    圖8 進(jìn)氣機(jī)匣流道優(yōu)化前后截面積對(duì)比Fig.8 The contrast of flow passage area before and after optimization

    4.2 全三維流場(chǎng)分析

    流道優(yōu)化后各截面總壓恢復(fù)系數(shù)、氣流角及馬赫數(shù)的徑向分布如圖2~圖4所示,優(yōu)化后流道中氣流馬赫數(shù)分布如圖9、圖10所示。

    從圖2中看,進(jìn)氣機(jī)匣流道優(yōu)化后,中部主流區(qū)的總壓恢復(fù)系數(shù)基本沒(méi)變,但根部和尖部明顯提高。從圖3中看,流道優(yōu)化前后,各截面氣流角徑向分布并無(wú)變化,說(shuō)明對(duì)流道的修改并未影響流道中的氣流方向。從圖4中看,支板前(1-1截面)主流區(qū)氣流馬赫數(shù)明顯降低,且與根部和尖部的馬赫數(shù)差值范圍縮小,說(shuō)明流場(chǎng)分布比優(yōu)化前均勻;支板后(2-2截面)的馬赫數(shù)徑向分布也比優(yōu)化前均勻。

    對(duì)比圖5和圖9發(fā)現(xiàn),流道優(yōu)化后支板葉片5%、50%、95%葉高附近的氣流馬赫數(shù)都明顯降低。對(duì)比圖6和圖10可以看到,流道優(yōu)化后從進(jìn)氣帽罩到支板后的整個(gè)進(jìn)氣機(jī)匣流道中,氣流速度都明顯降低。氣流在流過(guò)帽罩經(jīng)過(guò)喉道達(dá)到支板前都比較平穩(wěn),氣流速度變化較小,沒(méi)有明顯的局部高速氣流;氣流流過(guò)喉道后也未明顯減速,因此可降低損失。

    從表1中也可看出,與優(yōu)化前相比,優(yōu)化后進(jìn)氣機(jī)匣和支板的總壓恢復(fù)系數(shù)都有較大提高,流道后方壓氣機(jī)的工作效率也有一定程度的提高。

    圖9 流道優(yōu)化后擬S1流面馬赫數(shù)分布Fag.9 Relative Mach number distribution at S1 stream surface after optimization

    圖10 流道優(yōu)化后擬S2流面馬赫數(shù)分布Fag.10 Relative Mach number distribution at S2 stream surface after optimization

    5 結(jié)論

    (1)原進(jìn)氣機(jī)匣流道的總壓損失較大,總壓恢復(fù)系數(shù)不滿足指標(biāo)要求,需進(jìn)行優(yōu)化設(shè)計(jì)。

    (2)對(duì)進(jìn)氣機(jī)匣流道多輪優(yōu)化設(shè)計(jì)后,進(jìn)氣機(jī)匣流道和支板的總壓恢復(fù)系數(shù)都有明顯提高,且滿足指標(biāo)要求。

    (3)優(yōu)化后進(jìn)氣機(jī)匣流道中的流場(chǎng)得到很大改善,性能得以提高,說(shuō)明采用三次樣條插值方法能對(duì)流道進(jìn)行有效優(yōu)化。

    [1]Seddon J,Goldsmith E L.Intake aerodynamics[M].2nd ed.Washington:American Institute of Aeronautics and As?tronautics,Inc.,1999.

    [2]Mahoney J J.Inlets for supersonic missiles[M].Washing?ton:American Institute of Aeronautics and Astronautics, Inc.,1991.

    [3]陳瑩,王世安,鄒積國(guó),等.某型壓氣機(jī)高壓進(jìn)氣機(jī)匣的改進(jìn)設(shè)計(jì)[J].熱能動(dòng)力工程,2006,21(3):314—316.

    [4]蔡元虎,張建東,王占學(xué).TBCC發(fā)動(dòng)機(jī)用進(jìn)氣道設(shè)計(jì)及沿飛行軌跡斜板角度優(yōu)化分析[J].西北工業(yè)大學(xué)學(xué)報(bào),2007,25(5):615—619.

    [5]Baaley D W.Performance assessment of an annular S-shaped duct[J].Journal of Turbomachanery,1997,119:149—156.

    [6]李鵬,顧宏斌,高振興.三次樣條插值法在氣動(dòng)導(dǎo)數(shù)計(jì)算中的應(yīng)用[J].飛行力學(xué),2008,26(2):74—76.

    [7]張?zhí)m新,王長(zhǎng)富,劉占全,等.C語(yǔ)言程序設(shè)計(jì)與計(jì)算方法基礎(chǔ)[M].北京:北京航空航天大學(xué)出版社,2000.

    Evaluation and improvement of a TBCC engine inlet-case flow passage

    WANG Jun,YU Hua-wei,CAI Liu-cheng,YANG Wei
    (China Gas Turbine Establishment,Chengdu 610500,China)

    Firstly,the original flow-passage was computed by three-dimensional CFD software,and the flow field was analyzed.It could find that the airflow loss was too much in the original inlet-case flow pas?sage,and the total pressure recovery coefficient didn’t meet the design requirements,so the flow-passage should be improved.Then,the flow passage was optimized by the cubic spline method to make the cross sec?tion flatter.The optimized flow passage was computed by three-dimensional CFD software,and the flow field was analyzed.Finally,the optimized inlet-case was gained after several rounds of iteration,the total pressure recovery coefficient of the new flow passage and strut was higher and could meet the targets.

    TBCC engine;inlet case;cubic spline;optimization design;3D calculation

    V231.3

    A

    1672-2620(2015)01-0049-05

    2014-05-06;

    2014-07-30

    王俊(1983-),男,四川樂(lè)山人,工程師,碩士,主要從事壓氣機(jī)設(shè)計(jì)工作。

    猜你喜歡
    支板恢復(fù)系數(shù)流道
    剛體彈性碰撞中恢復(fù)系數(shù)的探討
    考慮河道影響的泵站進(jìn)水前池水力優(yōu)化設(shè)計(jì)研究
    利用恢復(fù)系數(shù)巧解碰撞問(wèn)題
    泵站底板與流道層側(cè)墻臺(tái)階式施工縫研究與應(yīng)用
    基于逆向氣體噴注的支板熱防護(hù)研究
    壁面噴射當(dāng)量比對(duì)支板凹腔耦合燃燒的影響
    一種可便捷更換外殼的計(jì)算機(jī)機(jī)箱
    主副流道旋轉(zhuǎn)微噴頭的研制和性能測(cè)試
    基于有限元的流道出口膠料速率均勻性分析
    落石碰撞法向恢復(fù)系數(shù)的模型試驗(yàn)研究
    97热精品久久久久久| av在线蜜桃| 午夜福利在线在线| 久久精品综合一区二区三区| 国产亚洲精品久久久com| 偷拍熟女少妇极品色| 18禁在线无遮挡免费观看视频 | 欧美人与善性xxx| 国产高潮美女av| 91麻豆精品激情在线观看国产| 一个人免费在线观看电影| 国产淫片久久久久久久久| 国产单亲对白刺激| 国产亚洲精品综合一区在线观看| 久久久午夜欧美精品| 十八禁网站免费在线| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美成人精品一区二区| 99国产极品粉嫩在线观看| 天堂√8在线中文| 一级毛片我不卡| 99热这里只有精品一区| 老司机福利观看| av卡一久久| 国内精品美女久久久久久| 最近最新中文字幕大全电影3| 精品无人区乱码1区二区| 国产真实乱freesex| 亚洲性夜色夜夜综合| 日本免费一区二区三区高清不卡| 免费无遮挡裸体视频| 成人一区二区视频在线观看| 精品午夜福利视频在线观看一区| 国产麻豆成人av免费视频| 精品免费久久久久久久清纯| 日韩三级伦理在线观看| 最新中文字幕久久久久| 午夜福利高清视频| 欧美不卡视频在线免费观看| 免费av观看视频| 麻豆久久精品国产亚洲av| av免费在线看不卡| 久久久国产成人免费| 国产又黄又爽又无遮挡在线| av.在线天堂| 亚洲综合色惰| 俺也久久电影网| 国产乱人偷精品视频| 成年版毛片免费区| 国产精品一及| 亚洲欧美日韩高清专用| 日韩一区二区视频免费看| 九九在线视频观看精品| 精品午夜福利在线看| 亚洲精品国产av成人精品 | 18禁在线播放成人免费| 国产片特级美女逼逼视频| 天堂影院成人在线观看| 人人妻人人澡人人爽人人夜夜 | 久久久色成人| 国产伦精品一区二区三区四那| 三级国产精品欧美在线观看| 日韩欧美三级三区| 国产一区二区在线av高清观看| 亚州av有码| 99国产极品粉嫩在线观看| 少妇被粗大猛烈的视频| 久久久久国内视频| 久久久精品大字幕| 男人舔奶头视频| 男女边吃奶边做爰视频| 麻豆av噜噜一区二区三区| 在线播放无遮挡| 亚洲精品亚洲一区二区| 国产黄色小视频在线观看| 啦啦啦啦在线视频资源| 亚洲色图av天堂| 久久精品夜夜夜夜夜久久蜜豆| 69av精品久久久久久| 成熟少妇高潮喷水视频| 高清午夜精品一区二区三区 | 高清毛片免费观看视频网站| 国内精品久久久久精免费| 不卡一级毛片| 中文字幕熟女人妻在线| 人妻久久中文字幕网| 久久久久国产精品人妻aⅴ院| 亚洲乱码一区二区免费版| 我要看日韩黄色一级片| 国产精品久久久久久久久免| 在线a可以看的网站| 日韩制服骚丝袜av| 精品欧美国产一区二区三| 国产精品一区二区免费欧美| 日韩av不卡免费在线播放| 男女边吃奶边做爰视频| 国产精品一二三区在线看| 午夜视频国产福利| 色在线成人网| 日日干狠狠操夜夜爽| 亚洲av电影不卡..在线观看| 中文亚洲av片在线观看爽| 国产精品免费一区二区三区在线| 久99久视频精品免费| 青春草视频在线免费观看| 免费在线观看影片大全网站| 18禁在线播放成人免费| 久久久久精品国产欧美久久久| 亚洲第一区二区三区不卡| 亚洲熟妇熟女久久| 国产精品伦人一区二区| 欧美日韩一区二区视频在线观看视频在线 | 亚洲性久久影院| 国产 一区精品| 天美传媒精品一区二区| 真实男女啪啪啪动态图| 国产男靠女视频免费网站| 我要看日韩黄色一级片| 国内久久婷婷六月综合欲色啪| 夜夜看夜夜爽夜夜摸| 国产av不卡久久| 午夜激情福利司机影院| h日本视频在线播放| 日本 av在线| 舔av片在线| 男人和女人高潮做爰伦理| 国产在视频线在精品| 国产精品久久久久久久电影| 国产一区二区三区在线臀色熟女| 精品99又大又爽又粗少妇毛片| 亚洲专区国产一区二区| 国产精品三级大全| 欧美另类亚洲清纯唯美| 12—13女人毛片做爰片一| 最近视频中文字幕2019在线8| 一级av片app| 三级经典国产精品| a级毛片免费高清观看在线播放| 最新在线观看一区二区三区| 你懂的网址亚洲精品在线观看 | 亚洲av.av天堂| 亚洲四区av| 联通29元200g的流量卡| 男人舔女人下体高潮全视频| 成人性生交大片免费视频hd| 伦精品一区二区三区| 久久这里只有精品中国| 国产精品美女特级片免费视频播放器| 特级一级黄色大片| 国产高清视频在线观看网站| 成熟少妇高潮喷水视频| 国内少妇人妻偷人精品xxx网站| 亚洲性夜色夜夜综合| 亚洲久久久久久中文字幕| 99视频精品全部免费 在线| 色哟哟·www| 成人鲁丝片一二三区免费| 国语自产精品视频在线第100页| 深夜精品福利| 欧美成人免费av一区二区三区| 夜夜看夜夜爽夜夜摸| 国产精品爽爽va在线观看网站| 晚上一个人看的免费电影| 国产视频一区二区在线看| 国国产精品蜜臀av免费| 一卡2卡三卡四卡精品乱码亚洲| 国产精品亚洲一级av第二区| 欧美一区二区亚洲| 精品久久久久久久久久免费视频| 美女大奶头视频| 亚洲精品久久国产高清桃花| 伦精品一区二区三区| 免费av毛片视频| 夜夜爽天天搞| 99热这里只有是精品在线观看| 99久久精品一区二区三区| 久久国产乱子免费精品| 能在线免费观看的黄片| 黄色日韩在线| 欧美一区二区亚洲| 老师上课跳d突然被开到最大视频| 22中文网久久字幕| 99国产精品一区二区蜜桃av| 久久久久久久久中文| 在线观看一区二区三区| 少妇熟女欧美另类| 女同久久另类99精品国产91| 日本免费a在线| 国产aⅴ精品一区二区三区波| 大型黄色视频在线免费观看| 草草在线视频免费看| 插阴视频在线观看视频| 国产亚洲精品综合一区在线观看| 亚洲精品色激情综合| 麻豆av噜噜一区二区三区| 全区人妻精品视频| 少妇裸体淫交视频免费看高清| 精品人妻一区二区三区麻豆 | 老司机影院成人| 国产精品亚洲美女久久久| 国产精品免费一区二区三区在线| 国产69精品久久久久777片| 99视频精品全部免费 在线| 中国国产av一级| 亚洲第一电影网av| 国产 一区精品| 午夜精品在线福利| 少妇熟女aⅴ在线视频| 亚洲熟妇中文字幕五十中出| 成人鲁丝片一二三区免费| 亚洲第一区二区三区不卡| 在线观看美女被高潮喷水网站| 久久人妻av系列| 亚洲av中文av极速乱| 麻豆国产97在线/欧美| 性插视频无遮挡在线免费观看| 亚洲电影在线观看av| 亚洲熟妇中文字幕五十中出| 夜夜夜夜夜久久久久| 久久久久免费精品人妻一区二区| 91久久精品电影网| 色综合站精品国产| 欧美成人精品欧美一级黄| 熟妇人妻久久中文字幕3abv| 又粗又爽又猛毛片免费看| 久久人人爽人人爽人人片va| 精品久久久久久久末码| 欧美3d第一页| 久久久久久久久久成人| 看非洲黑人一级黄片| 一边摸一边抽搐一进一小说| 国产又黄又爽又无遮挡在线| 熟妇人妻久久中文字幕3abv| 日韩欧美一区二区三区在线观看| 国产一区二区三区在线臀色熟女| 蜜桃亚洲精品一区二区三区| 黄色一级大片看看| 一a级毛片在线观看| 日韩中字成人| 免费看美女性在线毛片视频| 国产淫片久久久久久久久| 人妻丰满熟妇av一区二区三区| 一个人看视频在线观看www免费| 亚洲美女黄片视频| 有码 亚洲区| 在线观看午夜福利视频| 亚洲av不卡在线观看| 日韩国内少妇激情av| 精品人妻偷拍中文字幕| 天堂√8在线中文| 亚洲五月天丁香| 中出人妻视频一区二区| 国产亚洲精品综合一区在线观看| 一区二区三区高清视频在线| 深夜精品福利| 国内精品宾馆在线| 亚洲欧美日韩高清在线视频| 国产精品福利在线免费观看| 日韩av在线大香蕉| 久久草成人影院| 久久久久国产精品人妻aⅴ院| 男人舔女人下体高潮全视频| 免费av毛片视频| 麻豆av噜噜一区二区三区| 国产伦在线观看视频一区| 亚洲精品亚洲一区二区| 我要搜黄色片| 99热这里只有精品一区| 日本-黄色视频高清免费观看| 香蕉av资源在线| 国产精品女同一区二区软件| 亚洲人与动物交配视频| 亚洲人成网站在线观看播放| 国内少妇人妻偷人精品xxx网站| 亚洲一级一片aⅴ在线观看| 色播亚洲综合网| 亚洲图色成人| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品日韩av片在线观看| 色综合站精品国产| 久99久视频精品免费| 亚洲av成人精品一区久久| 在线观看av片永久免费下载| 九九热线精品视视频播放| 久久久a久久爽久久v久久| 国产一区二区在线观看日韩| 欧美成人精品欧美一级黄| 亚洲精品在线观看二区| 中文字幕精品亚洲无线码一区| 亚洲中文字幕日韩| 亚洲欧美日韩卡通动漫| 亚洲国产精品合色在线| 国产黄色小视频在线观看| 久久久欧美国产精品| 长腿黑丝高跟| 日韩大尺度精品在线看网址| 国产高清视频在线观看网站| 亚洲在线观看片| 日韩欧美一区二区三区在线观看| 特大巨黑吊av在线直播| 国产午夜精品久久久久久一区二区三区 | 69人妻影院| 午夜视频国产福利| 色哟哟·www| 精品久久久久久久久久免费视频| 久久精品国产清高在天天线| 国产av一区在线观看免费| 黄片wwwwww| 欧美3d第一页| 深夜a级毛片| 在线看三级毛片| 国产高潮美女av| 国产亚洲精品久久久com| 变态另类成人亚洲欧美熟女| 成人午夜高清在线视频| 少妇熟女aⅴ在线视频| 99九九线精品视频在线观看视频| 色5月婷婷丁香| 久久久精品94久久精品| 久久久欧美国产精品| 九色成人免费人妻av| 亚洲av成人精品一区久久| 香蕉av资源在线| 高清毛片免费看| 精品久久久久久久久亚洲| 美女高潮的动态| 波多野结衣高清无吗| 亚洲欧美中文字幕日韩二区| 老熟妇仑乱视频hdxx| 精品欧美国产一区二区三| 夜夜看夜夜爽夜夜摸| 一边摸一边抽搐一进一小说| 黄色日韩在线| 在线观看午夜福利视频| 欧美日韩一区二区视频在线观看视频在线 | 免费高清视频大片| 久久久久精品国产欧美久久久| 欧美日韩精品成人综合77777| 久久久国产成人精品二区| 神马国产精品三级电影在线观看| 国产视频内射| 村上凉子中文字幕在线| 性欧美人与动物交配| 一进一出抽搐gif免费好疼| 欧美色视频一区免费| 精品不卡国产一区二区三区| 午夜免费激情av| 亚洲人成网站在线播放欧美日韩| 国产亚洲av嫩草精品影院| 高清午夜精品一区二区三区 | 国产色爽女视频免费观看| 国产精品久久久久久精品电影| 欧美另类亚洲清纯唯美| 国产蜜桃级精品一区二区三区| 亚洲人成网站在线播| 国产精品无大码| 日本黄色片子视频| 亚洲精品成人久久久久久| 免费看av在线观看网站| 亚洲精品成人久久久久久| 91精品国产九色| 国产一区二区三区av在线 | 久久久久国产网址| 欧美在线一区亚洲| 变态另类成人亚洲欧美熟女| 春色校园在线视频观看| 悠悠久久av| 别揉我奶头~嗯~啊~动态视频| 99久久久亚洲精品蜜臀av| 麻豆国产av国片精品| 国产女主播在线喷水免费视频网站 | 97超碰精品成人国产| 淫秽高清视频在线观看| 91狼人影院| 国产精品女同一区二区软件| 成人高潮视频无遮挡免费网站| 亚洲七黄色美女视频| 日韩,欧美,国产一区二区三区 | 深爱激情五月婷婷| 久久鲁丝午夜福利片| 午夜亚洲福利在线播放| 少妇高潮的动态图| 成人无遮挡网站| 精品不卡国产一区二区三区| 国产精品嫩草影院av在线观看| 国产午夜福利久久久久久| 99久国产av精品国产电影| 色哟哟哟哟哟哟| 国产人妻一区二区三区在| 大型黄色视频在线免费观看| 美女高潮的动态| 亚洲精品日韩av片在线观看| 亚洲久久久久久中文字幕| 草草在线视频免费看| 级片在线观看| 亚洲精品色激情综合| 青春草视频在线免费观看| 男女边吃奶边做爰视频| 欧美日韩国产亚洲二区| 日日摸夜夜添夜夜爱| 国产精品1区2区在线观看.| 女人十人毛片免费观看3o分钟| 真人做人爱边吃奶动态| 久久婷婷人人爽人人干人人爱| 欧美中文日本在线观看视频| 你懂的网址亚洲精品在线观看 | 在线观看av片永久免费下载| 久久久久国内视频| 欧美人与善性xxx| 国产精品人妻久久久影院| 一级av片app| 我的老师免费观看完整版| 国产一区二区亚洲精品在线观看| 亚洲欧美精品综合久久99| 美女大奶头视频| 免费黄网站久久成人精品| 免费av毛片视频| 精品久久久久久久久久免费视频| 欧美性猛交黑人性爽| 日日摸夜夜添夜夜添av毛片| 女人十人毛片免费观看3o分钟| 国产精品1区2区在线观看.| 男插女下体视频免费在线播放| av在线蜜桃| 成人性生交大片免费视频hd| 午夜福利18| 如何舔出高潮| 日日撸夜夜添| 国产激情偷乱视频一区二区| 偷拍熟女少妇极品色| 黄色一级大片看看| 国产亚洲精品久久久com| 国产爱豆传媒在线观看| 亚洲一区高清亚洲精品| 99久久精品热视频| 免费人成在线观看视频色| or卡值多少钱| 老熟妇乱子伦视频在线观看| 午夜精品在线福利| 日本五十路高清| 欧美成人免费av一区二区三区| 国产蜜桃级精品一区二区三区| 99热这里只有是精品50| 舔av片在线| 黄色配什么色好看| 免费观看的影片在线观看| 22中文网久久字幕| 看黄色毛片网站| 人人妻人人澡欧美一区二区| 精品99又大又爽又粗少妇毛片| av免费在线看不卡| 小蜜桃在线观看免费完整版高清| 狠狠狠狠99中文字幕| 三级男女做爰猛烈吃奶摸视频| 99久国产av精品国产电影| 嫩草影视91久久| 亚洲精品456在线播放app| 久久婷婷人人爽人人干人人爱| 成人美女网站在线观看视频| 亚洲性夜色夜夜综合| 伦理电影大哥的女人| 禁无遮挡网站| 亚洲精品在线观看二区| 九九爱精品视频在线观看| 亚洲18禁久久av| 国产精品人妻久久久久久| 亚洲七黄色美女视频| 亚洲人成网站在线播放欧美日韩| 久久精品夜色国产| 国产不卡一卡二| 国产高清三级在线| 亚洲性夜色夜夜综合| 国产白丝娇喘喷水9色精品| 国产亚洲av嫩草精品影院| 亚洲av中文av极速乱| 人妻久久中文字幕网| 久久久久久久久久黄片| 久久亚洲国产成人精品v| 99久久无色码亚洲精品果冻| 我的老师免费观看完整版| 最新在线观看一区二区三区| 香蕉av资源在线| 看片在线看免费视频| 在线观看美女被高潮喷水网站| 国产高清不卡午夜福利| 亚洲人与动物交配视频| 午夜a级毛片| 天堂网av新在线| 欧美性猛交╳xxx乱大交人| 人人妻人人澡人人爽人人夜夜 | 在线观看66精品国产| 在线免费十八禁| 色在线成人网| av在线老鸭窝| 超碰av人人做人人爽久久| 国产老妇女一区| 美女 人体艺术 gogo| avwww免费| 啦啦啦观看免费观看视频高清| 一级毛片久久久久久久久女| 高清午夜精品一区二区三区 | 伦精品一区二区三区| 日日啪夜夜撸| 蜜桃久久精品国产亚洲av| 久久精品影院6| 成人永久免费在线观看视频| 在线观看午夜福利视频| 此物有八面人人有两片| 亚洲久久久久久中文字幕| 国产欧美日韩精品亚洲av| 国产精品人妻久久久久久| 嫩草影院入口| 日韩精品青青久久久久久| 蜜桃亚洲精品一区二区三区| 天堂动漫精品| 久久鲁丝午夜福利片| 久久久色成人| 精品久久久久久久久亚洲| 啦啦啦观看免费观看视频高清| 校园人妻丝袜中文字幕| 精品一区二区三区人妻视频| 国产v大片淫在线免费观看| 99热网站在线观看| 国产激情偷乱视频一区二区| 舔av片在线| 亚洲性久久影院| 亚洲精品日韩在线中文字幕 | 国产成人福利小说| 国产av麻豆久久久久久久| 国产精品爽爽va在线观看网站| 欧美人与善性xxx| 毛片一级片免费看久久久久| 看十八女毛片水多多多| 夜夜爽天天搞| 免费无遮挡裸体视频| 日本免费a在线| 国产探花极品一区二区| 久久久久久大精品| 国产色婷婷99| 美女黄网站色视频| 亚洲在线自拍视频| 变态另类丝袜制服| 91在线观看av| 亚洲av第一区精品v没综合| 少妇高潮的动态图| 久久久久九九精品影院| 三级国产精品欧美在线观看| 麻豆成人午夜福利视频| 天天躁夜夜躁狠狠久久av| 啦啦啦啦在线视频资源| 国产不卡一卡二| 尤物成人国产欧美一区二区三区| 亚洲一级一片aⅴ在线观看| 国产色婷婷99| 97在线视频观看| 国产在线精品亚洲第一网站| 亚洲欧美成人综合另类久久久 | 男人舔女人下体高潮全视频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一及| 亚洲丝袜综合中文字幕| 亚洲最大成人av| 国产伦在线观看视频一区| 嫩草影院新地址| av福利片在线观看| 亚洲精品国产av成人精品 | 国产精品美女特级片免费视频播放器| 人人妻,人人澡人人爽秒播| 亚洲精品456在线播放app| 精品久久久久久久久亚洲| 午夜精品国产一区二区电影 | 亚洲图色成人| 高清日韩中文字幕在线| 99久久成人亚洲精品观看| 国产一区亚洲一区在线观看| 日本一本二区三区精品| 亚洲18禁久久av| av专区在线播放| 国产精品精品国产色婷婷| 亚洲av二区三区四区| 成人av在线播放网站| 看非洲黑人一级黄片| 午夜精品一区二区三区免费看| 欧美性猛交黑人性爽| 三级经典国产精品| 99在线视频只有这里精品首页| 村上凉子中文字幕在线| 国产久久久一区二区三区| 亚州av有码| 亚洲丝袜综合中文字幕| 黄色视频,在线免费观看| 婷婷精品国产亚洲av| a级毛色黄片| 免费不卡的大黄色大毛片视频在线观看 | 午夜福利高清视频| 久久久久久久久久久丰满| 一级毛片久久久久久久久女| 国产高潮美女av| 精品人妻一区二区三区麻豆 | 麻豆一二三区av精品| 一本久久中文字幕| 色噜噜av男人的天堂激情| 免费av观看视频| 久久国内精品自在自线图片| www日本黄色视频网| 国产精品人妻久久久久久| 亚洲av免费高清在线观看| 久久亚洲精品不卡| 成人永久免费在线观看视频| 能在线免费观看的黄片| 精品国产三级普通话版| 99精品在免费线老司机午夜| 我要搜黄色片| 在线a可以看的网站| 一本一本综合久久|