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

    雷諾數(shù)對(duì)渦槳飛機(jī)增升裝置氣動(dòng)特性影響的計(jì)算研究

    2016-12-19 02:55:22許瑞飛錢瑞戰(zhàn)雷武濤
    航空工程進(jìn)展 2016年4期
    關(guān)鍵詞:附面層渦槳迎角

    許瑞飛,錢瑞戰(zhàn),雷武濤

    (中國(guó)航空工業(yè)集團(tuán)公司 第一飛機(jī)設(shè)計(jì)研究院,西安 710089)

    ?

    雷諾數(shù)對(duì)渦槳飛機(jī)增升裝置氣動(dòng)特性影響的計(jì)算研究

    許瑞飛,錢瑞戰(zhàn),雷武濤

    (中國(guó)航空工業(yè)集團(tuán)公司 第一飛機(jī)設(shè)計(jì)研究院,西安 710089)

    雷諾數(shù)是表征流體粘性對(duì)流動(dòng)影響的相似參數(shù),對(duì)飛機(jī)部件的氣動(dòng)性能具有重要的影響。利用商業(yè)軟件CFX對(duì)某民用渦槳飛機(jī)增升裝置構(gòu)型進(jìn)行變雷諾數(shù)計(jì)算研究,詳細(xì)分析雷諾數(shù)對(duì)升力系數(shù)、失速特性以及附面層速度特性的影響。結(jié)果表明:雷諾數(shù)對(duì)最大升力系數(shù)和失速迎角都有顯著的影響;不同雷諾數(shù)下機(jī)翼分離形態(tài)變化顯著,大雷諾數(shù)下機(jī)翼分離區(qū)域較??;不同雷諾數(shù)下機(jī)翼附面層狀態(tài)不一致,雷諾數(shù)增大使得附面層速度型更飽滿,機(jī)翼的抗分離能力增強(qiáng)。

    渦槳飛機(jī);增升裝置;雷諾數(shù);氣動(dòng)特性

    0 引 言

    隨著飛機(jī)和發(fā)動(dòng)機(jī)技術(shù)的不斷發(fā)展,渦槳飛機(jī)的巡航速度逐漸提高,在400 nmile的航段上運(yùn)營(yíng)的經(jīng)濟(jì)性明顯高于噴氣式飛機(jī)。這一技術(shù)優(yōu)勢(shì)有助于航空公司開辟新市場(chǎng)。采用主動(dòng)噪聲和振動(dòng)抑制系統(tǒng),渦槳飛機(jī)的客艙內(nèi)部噪聲大幅降低,乘客舒適性顯著提高。因此在高頻率、短航程及低成本市場(chǎng)上渦槳飛機(jī)有很大的發(fā)展空間。

    渦槳飛機(jī)的設(shè)計(jì)過程需要掌握一些關(guān)鍵技術(shù)。其中,復(fù)雜增升裝置的設(shè)計(jì)技術(shù)[1-3]是氣動(dòng)設(shè)計(jì)的關(guān)鍵技術(shù)之一,直接影響著飛機(jī)的起飛距離、著陸距離、起飛重量等[4]。在增升裝置氣動(dòng)性能的影響因素中,雷諾數(shù)是一個(gè)十分重要的參數(shù),它表示為慣性力與粘性力之比,是一個(gè)表征流體粘性對(duì)流動(dòng)影響的相似參數(shù),對(duì)增升裝置氣動(dòng)性能有顯著影響。

    國(guó)外型號(hào)設(shè)計(jì)單位如波音、空客等都對(duì)雷諾數(shù)影響規(guī)律進(jìn)行了大量的研究。F.M.Payne等[5]通過NTF風(fēng)洞對(duì)B777-200飛機(jī)的起飛和著陸構(gòu)型進(jìn)行了變雷諾數(shù)試驗(yàn),并利用MDOE方法對(duì)外襟翼縫道參數(shù)進(jìn)行了變雷諾數(shù)研究。M.M.Curtin等[6]對(duì)B777飛機(jī)進(jìn)行了不同雷諾數(shù)的測(cè)力和測(cè)壓試驗(yàn)研究,并將飛行雷諾數(shù)得到的試驗(yàn)數(shù)據(jù)與試飛數(shù)據(jù)比較,試圖建立風(fēng)洞試驗(yàn)數(shù)據(jù)至試飛數(shù)據(jù)的修正。R.W.Clark等[7]對(duì)MD-11飛機(jī)巡航狀態(tài)進(jìn)行了變雷諾數(shù)研究,并詳細(xì)分析了不同雷諾數(shù)下升力、阻力和力矩系數(shù)的變化。J.Quest等[8]利用ETW風(fēng)洞對(duì)空客某高升力構(gòu)型進(jìn)行了高雷諾數(shù)風(fēng)洞試驗(yàn),以探索雷諾數(shù)影響規(guī)律。為了研究雷諾數(shù)對(duì)飛機(jī)的影響規(guī)律,歐洲EUROLIFT計(jì)劃專門針對(duì)某高升力構(gòu)型進(jìn)行雷諾數(shù)影響研究[9-10]。鄭隆乾等[11]對(duì)某民用飛機(jī)增升裝置著陸構(gòu)型進(jìn)行了半模高雷諾數(shù)風(fēng)洞試驗(yàn),獲得了不同馬赫數(shù)下的雷諾數(shù)影響。

    雷諾數(shù)對(duì)增升裝置氣動(dòng)特性有顯著影響。升力系數(shù)、阻力系數(shù)、升阻比等氣動(dòng)參數(shù)在飛機(jī)設(shè)計(jì)中極為重要,直接關(guān)系到型號(hào)設(shè)計(jì)的成敗。因此,研究雷諾數(shù)對(duì)增升裝置的影響具有現(xiàn)實(shí)意義,研究成果可直接應(yīng)用于型號(hào)設(shè)計(jì)。

    本文利用CFD技術(shù)對(duì)某渦槳飛機(jī)著陸翼身組合體構(gòu)型進(jìn)行數(shù)值模擬。首先將計(jì)算結(jié)果與試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比,然后研究雷諾數(shù)對(duì)某渦槳飛機(jī)升力系數(shù)、失速特性以及附面層速度特性的影響,詳細(xì)分析雷諾數(shù)對(duì)飛機(jī)氣動(dòng)特性的影響規(guī)律,以期為渦槳飛機(jī)氣動(dòng)力設(shè)計(jì)和風(fēng)洞試驗(yàn)數(shù)據(jù)的相關(guān)性修正提供參考。

    1 計(jì)算方法

    國(guó)內(nèi)外已經(jīng)開展了有關(guān)增升裝置的數(shù)值計(jì)算方法研究[12-15]。本文采用商業(yè)軟件CFX進(jìn)行飛機(jī)增升裝置的數(shù)值模擬計(jì)算。CFX軟件是基于有限元的有限體積法,用于求解雷諾平均N-S方程。基于有限元的有限體積法兼具有限體積法的守恒特性和有限元法的數(shù)值精確性。求解控制方程為

    2 計(jì)算模型與網(wǎng)格

    計(jì)算模型采用某民用渦槳飛機(jī)著陸翼身組合體構(gòu)型。該飛機(jī)增升裝置只有后緣襟翼,無前緣縫翼,如圖1所示。

    圖1 某民用渦槳飛機(jī)著陸翼身組合體模型

    計(jì)算網(wǎng)格采用H-O拓?fù)湫问?,在物面附近生成O型網(wǎng)格,遠(yuǎn)場(chǎng)采用H型網(wǎng)格。整個(gè)流場(chǎng)域網(wǎng)格變化率合適,疏密分布合理,過渡均勻,滿足計(jì)算要求,表面網(wǎng)格如圖2所示,網(wǎng)格量為2 000萬左右。

    圖2 某民用渦槳飛機(jī)著陸翼身組合體構(gòu)型表面網(wǎng)格

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

    3.1 計(jì)算與試驗(yàn)驗(yàn)證

    測(cè)試算例采用某民用渦槳飛機(jī)風(fēng)洞試驗(yàn)?zāi)P?。試?yàn)?zāi)P椭皇窃趫D1模型(即本文計(jì)算模型)的基礎(chǔ)上對(duì)短艙進(jìn)行了下移,其他部件數(shù)模無變化。計(jì)算網(wǎng)格生成與第2節(jié)類似。計(jì)算馬赫數(shù)Ma=0.2,雷諾數(shù)Re=7.4×106,計(jì)算與試驗(yàn)結(jié)果的比較如圖3所示。

    (a) 升力曲線

    (b) 阻力曲線

    從圖3可以看出:計(jì)算得到的升力系數(shù)曲線與試驗(yàn)在小迎角吻合較好,失速迎角和最大升力系數(shù)值略有不同,失速發(fā)展趨勢(shì)基本一致,表明本文計(jì)算方法能夠很好地模擬增升裝置的氣動(dòng)特性和流動(dòng)特性。

    3.2 不同雷諾數(shù)下力系數(shù)和壓力分布比較

    為了研究雷諾數(shù)對(duì)某民用渦槳飛機(jī)增升裝置的影響,基于不同雷諾數(shù)(1.9×106,3.0×106,4.0×106,6.0×106,8.0×106,1.17×107)對(duì)著陸翼身組合體構(gòu)型進(jìn)行計(jì)算研究。計(jì)算馬赫數(shù)Ma=0.2,不同雷諾數(shù)下某民用渦槳飛機(jī)著陸翼身組合體構(gòu)型的升力系數(shù)和升阻極曲線如圖4所示。

    (a) 升力曲線

    (b) 升阻極曲線

    從圖4可以看出:①隨著雷諾數(shù)的增加,失速迎角和最大升力系數(shù)均顯著增加;在小雷諾數(shù)下,升力系數(shù)的增長(zhǎng)幅度比較大;隨著雷諾數(shù)增加,粘性效應(yīng)降低,最大升力系數(shù)的增長(zhǎng)放緩,在小迎角下,升力系數(shù)獲得提升,這主要是由于雷諾數(shù)的增加使得附面層厚度降低,飛機(jī)上的流動(dòng)分離減輕,從而使得升力系數(shù)獲得一定程度的增加,飛機(jī)失速迎角也推遲;②從圖4(b)可以看出,在相同升力條件下,阻力降低,可以提高飛機(jī)的升阻比。

    通過計(jì)算得到的某民用渦槳飛機(jī)最大升力系數(shù)和失速迎角隨雷諾數(shù)的變化曲線如圖5所示。

    (a) 最大升力系數(shù)隨雷諾數(shù)變化曲線

    (b) 失速迎角隨雷諾數(shù)變化曲線

    從圖5可以看出:隨著雷諾數(shù)的增加,最大升力系數(shù)和失速迎角逐漸增加;在雷諾數(shù)為1.17×107時(shí),最大升力系數(shù)達(dá)到最大值。

    為了深入研究雷諾數(shù)對(duì)機(jī)翼氣動(dòng)特性的影響,截取機(jī)翼展向18.0%、43.0%、64.5%、82.0%站位的壓力分布進(jìn)行比較。迎角4°下不同雷諾數(shù)機(jī)翼展向站位壓力分布比較如圖6所示。

    (a) 18.0%站位

    (b) 43.0%站位

    (c) 64.5%站位

    (d) 82.0%站位

    從圖6可以看出:隨著雷諾數(shù)的增加,主翼和后緣襟翼的吸力峰值逐漸增大;主翼后緣壓力恢復(fù)提升,襟翼上的流動(dòng)分離區(qū)域逐漸減小,由此可對(duì)大雷諾數(shù)下升力系數(shù)較大的主要原因作出解釋。

    3.3 不同雷諾數(shù)下失速特性比較

    迎角16°時(shí)著陸翼身組合體構(gòu)型不同雷諾數(shù)下的表面壓力云圖和流譜如圖7所示。

    (a) Re=1.9×106

    (b) Re=3.0×106

    (c) Re=6.0×106

    (d) Re=1.17×107

    從圖7可以看出:在小雷諾數(shù)1.9×106下機(jī)翼的中內(nèi)側(cè)出現(xiàn)了大面積的流動(dòng)分離;當(dāng)雷諾數(shù)增大到3.0×106,機(jī)翼上表面出現(xiàn)較小的流動(dòng)分離,同時(shí)在短艙后面存在向機(jī)翼兩側(cè)擴(kuò)展的橫向流動(dòng);當(dāng)雷諾數(shù)達(dá)到6.0×106以后,機(jī)翼上表面未出現(xiàn)明顯的嚴(yán)重分離。

    為了進(jìn)一步理清機(jī)理,選擇機(jī)翼后緣某點(diǎn),截取當(dāng)?shù)厮俣刃瓦M(jìn)行分析。迎角16°下不同雷諾數(shù)機(jī)翼后緣某點(diǎn)附面層速度型分布的比較如圖8所示。

    (a) Re=1.9×106(b) Re=3.0×106

    (c) Re=6.0×106(d) Re=1.17×107

    圖8 迎角16°下不同雷諾數(shù)機(jī)翼后緣某點(diǎn)附面層速度型比較

    Fig.8ThecomparisonofvelocitydistributionfortrailingedgeofwingunderdifferentReynoldsnumberatangleofattack16°

    從圖8可以看出:當(dāng)雷諾數(shù)為1.9×106時(shí),此時(shí)當(dāng)?shù)馗矫鎸颖容^厚,接近壁面的速度值比較低;這種速度形態(tài)使得當(dāng)?shù)睾缶壓茈y抵抗較強(qiáng)的逆壓梯度,容易造成流動(dòng)分離;隨著雷諾數(shù)的增加,當(dāng)?shù)厮俣刃驮絹碓斤枬M,接近壁面的速度也逐漸提高,使得壁面的抗分離能力增強(qiáng);在大雷諾數(shù)下機(jī)翼能夠抵抗較大的逆壓梯度,分離較小,這也是在大雷諾數(shù)下機(jī)翼分離區(qū)域比較小的主要原因。

    4 結(jié) 論

    (1) 隨著雷諾數(shù)的增加,失速迎角和最大升力系數(shù)均顯著增加;飛機(jī)阻力降低,升阻比提高,氣動(dòng)性能顯著增強(qiáng)。在大雷諾數(shù)下,升力系數(shù)的增長(zhǎng)幅度明顯放緩。

    (2) 隨著雷諾數(shù)的增加,主翼和后緣襟翼的吸力峰值逐漸增大;主翼后緣壓力恢復(fù)提升,襟翼上的流動(dòng)分離區(qū)域逐漸減小。負(fù)壓峰值的提升和流動(dòng)分離的減小是大雷諾數(shù)下升力系數(shù)較大的主要原因。

    (3) 低雷諾數(shù)下,由于機(jī)翼后緣附面層較厚,附面層速度型不飽滿,造成機(jī)翼后緣抗分離能力較弱,易于形成流動(dòng)分離;隨著雷諾數(shù)的增加,速度型更加飽滿,抵抗逆壓梯度能力增強(qiáng),從而降低了流動(dòng)分離區(qū)域。

    (4) 不同雷諾數(shù)下某民用渦槳飛機(jī)分離形態(tài)不一致;大雷諾數(shù)下分離較小,但是總體來說,各雷諾數(shù)下分離位置基本一致,均在機(jī)翼短艙內(nèi)側(cè)開始分離。

    綜上所述,雷諾數(shù)對(duì)民用渦槳飛機(jī)增升裝置氣動(dòng)特性有著顯著影響。常規(guī)風(fēng)洞的試驗(yàn)雷諾數(shù)很難達(dá)到飛行雷諾數(shù),給氣動(dòng)數(shù)據(jù)的準(zhǔn)確確定帶來了一定的困難。因此建議在型號(hào)設(shè)計(jì)中,盡可能開展高雷諾數(shù)風(fēng)洞試驗(yàn),提高氣動(dòng)力數(shù)據(jù)的準(zhǔn)確性。

    [1]DamvanCP.Theaerodynamicdesignofmulti-elementhigh-liftsystemsfortransportairplanes[J].ProgressinAerospaceSciences, 2002, 38(2): 101-144.

    [2] Wild J. Validation of numerical optimization of high-lift multi-element airfoils based on navier-stokes-equations[C]. AIAA-2002-2939, 2002.

    [3] Kim S, Alonso J J, Jameson A. Multi-element high-lift configuration design optimization using viscous continuous adjoint method[J]. Journal of Aircraft, 2004, 41(5): 1082-1097.

    [4] Smith A M O. High lift aerodynamics[J]. Journal of Aircraft, 1975, 12(6): 511-530.

    [5] Payne F M , Wyatt G W, Bogue D R, et al. High Reynolds number studies of a Boeing 777-200 high lift configuration in the NASA ARC 12’ pressure tunnel and NASA LaRC National Transonic Facility[C]. AIAA-2000-4220, 2000.

    [6] Curtin M M, Bogue D R, Om D, et al. Investigation of transonic Reynolds number scaling on a twin-engine transport[C]. AIAA-2002-0420, 2002.

    [7] Clark R W, Pelkmant R A. High Reynolds number testing of advanced transport aircraft wings in the National Transonic Facility[C]. AIAA-2001-0910, 2001.

    [8] Quest J, Hansen H, Mesuro G G. First measurements on an airbus high lift configuration at ETW up to flight Reynolds number[C]. AIAA-2002-0423, 2002.

    [9] Rudnik R, Germain E. Re-No. scaling sffects on the EUROLIFT high lift configurations[C]. AIAA-2007-752, 2007.

    [10] Rudnik R. Stall behaviour of the EUROLIFT high lift configurations[R]. AIAA-2008-0836, 2008.

    [11] 鄭隆乾, 陳迎春, 李亞林, 等. 某民機(jī)著陸構(gòu)型雷諾數(shù)效應(yīng)風(fēng)洞試驗(yàn)研究[J]. 航空計(jì)算技術(shù), 2012, 42(6): 70-71. Zheng Longqian, Chen Yingchun, Li Yalin, et al. Wind tunnel research about Reynolds number effects on landing configuration of a civil aircraft[J]. Aeronautical Computing Technique, 2012, 42(6): 70-71.(in Chinese)

    [12] 朱自強(qiáng). 應(yīng)用計(jì)算流體力學(xué)[M]. 北京: 北京航空航天大學(xué)出版社, 1998. Zhu Ziqiang. Application computational fluid dynamics[M]. Beijing: Beihang University Press, 1998.(in Chinese)

    [13] 朱自強(qiáng), 陳迎春, 吳宗成, 等. 高升力系統(tǒng)外形的數(shù)值模擬計(jì)算[J]. 航空學(xué)報(bào), 2005, 26(3): 257-262. Zhu Ziqiang, Chen Yingchun, Wu Zongcheng, et al. Numerical simulation of high lift system configuration[J]. Acta Aeronautica et Astronautica Sinica, 2005, 26(3): 257-262.(in Chinese)

    [14] 左志成, 段卓毅. 運(yùn)輸類飛機(jī)高升力構(gòu)型數(shù)值模擬研究[J]. 飛機(jī)設(shè)計(jì), 2007, 27(1): 16-21. Zuo Zhicheng, Duan Zhouyi. Numerical simulation of high lift system for the transport aircraft[J]. Aircraft Design, 2007, 27(1): 16-21.(in Chinese)

    [15] Dodbele S D. Three dimensional aerodynamic analysis of a high-lift transport configuration[R]. AIAA-93-3536, 1993.

    (編輯:趙毓梅)

    Numerical Investigation of Reynolds Number Effects on Aerodynamics Characteristics of High Lift System for Turboprop Aircraft

    Xu Ruifei, Qian Ruizhan, Lei Wutao

    (The First Aircraft Design and Research Institute, Aviation Industry Corporation of China, Xi’an 710089, China)

    Reynolds number is one of the most useful similarity parameters in fluid dynamics which denotes the effect of the viscosity on the fluid flow. It has an important effect on the aerodynamic characteristic of the aircraft. The variable Reynolds number research of the high lift for one of the civil turboprop aircraft is completed by the commercial software CFX. The effect of Reynolds number on the lift coefficient, stall characteristic, and the velocity characteristic of the boundary layer are also analyzed in detail. The research indicates that the Reynolds number has the notable effect on the maximum lift coefficient and the stall angle. The state of the flow separation on the wing is changed obviously under the different Reynolds number. The area of the separation becomes smaller when the Reynolds number is increased. The state of the boundary layer is different under the different Reynolds number. When the Reynolds number increases, the velocity profile is fuller, and the trend of separation is limited.

    turboprop aircraft; high lift; Reynolds number; aerodynamic characteristic

    2016-07-11;

    2016-10-09

    航空科學(xué)基金(20143203002)

    許瑞飛,xrf803@163.com

    1674-8190(2016)04-420-06

    V211

    A

    10.16615/j.cnki.1674-8190.2016.04.004

    許瑞飛(1985-),男,碩士,高級(jí)工程師。主要研究方向:飛機(jī)氣動(dòng)力設(shè)計(jì)與計(jì)算。

    錢瑞戰(zhàn)(1966-),男,博士,研究員。主要研究方向:飛機(jī)氣動(dòng)力設(shè)計(jì)。

    雷武濤(1975-),男,碩士,研究員。主要研究方向:飛機(jī)總體氣動(dòng)設(shè)計(jì)。

    猜你喜歡
    附面層渦槳迎角
    渦槳發(fā)動(dòng)機(jī)安裝技術(shù)研究
    基于網(wǎng)格框架的非結(jié)構(gòu)附面層網(wǎng)格生成技術(shù)
    連續(xù)變迎角試驗(yàn)數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    基于數(shù)值模擬的流場(chǎng)附面層邊緣識(shí)別方法
    民用飛機(jī)設(shè)計(jì)參考機(jī)種之一 ATR 72雙發(fā)渦槳支線運(yùn)輸機(jī)
    一例渦槳飛機(jī)飛行事故中的人因分析與研究
    探究渦槳飛機(jī)的客戶服務(wù)
    超聲壓氣機(jī)葉柵流場(chǎng)的數(shù)值模擬與試驗(yàn)驗(yàn)證
    失速保護(hù)系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    發(fā)動(dòng)機(jī)進(jìn)口附面層測(cè)量試驗(yàn)與數(shù)值模擬
    啦啦啦免费观看视频1| 999久久久国产精品视频| 亚洲免费av在线视频| 久久久久性生活片| 老司机午夜十八禁免费视频| 少妇裸体淫交视频免费看高清 | 国产精品九九99| 天天躁狠狠躁夜夜躁狠狠躁| 女同久久另类99精品国产91| 成人一区二区视频在线观看| 午夜老司机福利片| 免费人成视频x8x8入口观看| 男女床上黄色一级片免费看| 欧美av亚洲av综合av国产av| 国产午夜福利久久久久久| 人妻丰满熟妇av一区二区三区| 可以在线观看毛片的网站| 国产精品一区二区三区四区免费观看 | 99久久无色码亚洲精品果冻| 国产欧美日韩一区二区三| 男人舔女人的私密视频| 国产乱人伦免费视频| 精品久久久久久久人妻蜜臀av| 欧美日韩亚洲国产一区二区在线观看| 国产欧美日韩一区二区精品| 别揉我奶头~嗯~啊~动态视频| 成年人黄色毛片网站| 亚洲人与动物交配视频| 成年免费大片在线观看| bbb黄色大片| 久久久久国内视频| a级毛片在线看网站| 九九热线精品视视频播放| 精品国产超薄肉色丝袜足j| 色综合婷婷激情| 亚洲欧美日韩高清在线视频| 午夜激情福利司机影院| 国产高清视频在线观看网站| 午夜免费成人在线视频| 国产精品99久久99久久久不卡| 免费在线观看日本一区| 国产区一区二久久| 日韩大码丰满熟妇| 91成年电影在线观看| 久久亚洲精品不卡| 身体一侧抽搐| 99国产精品一区二区三区| 黑人欧美特级aaaaaa片| 久久99热这里只有精品18| 12—13女人毛片做爰片一| 超碰成人久久| 亚洲在线自拍视频| 国产成人欧美在线观看| 在线观看免费日韩欧美大片| 亚洲人成77777在线视频| 久久精品aⅴ一区二区三区四区| 12—13女人毛片做爰片一| www日本在线高清视频| 怎么达到女性高潮| 精品久久蜜臀av无| 最近在线观看免费完整版| 国产黄色小视频在线观看| 亚洲精品色激情综合| 日韩欧美精品v在线| 精品不卡国产一区二区三区| 亚洲av成人一区二区三| 校园春色视频在线观看| 在线观看免费午夜福利视频| 人成视频在线观看免费观看| 麻豆国产97在线/欧美 | 女生性感内裤真人,穿戴方法视频| 18禁黄网站禁片免费观看直播| 老汉色∧v一级毛片| 一本大道久久a久久精品| 黄色a级毛片大全视频| 狂野欧美激情性xxxx| 桃红色精品国产亚洲av| 免费观看人在逋| 一二三四社区在线视频社区8| 啦啦啦免费观看视频1| 麻豆久久精品国产亚洲av| 国产精品一区二区三区四区免费观看 | 免费人成视频x8x8入口观看| 久久久水蜜桃国产精品网| 日本 av在线| 亚洲精品在线观看二区| 一区二区三区激情视频| 全区人妻精品视频| 女人高潮潮喷娇喘18禁视频| 五月玫瑰六月丁香| 一二三四社区在线视频社区8| 国产蜜桃级精品一区二区三区| 男女床上黄色一级片免费看| 叶爱在线成人免费视频播放| 久久久水蜜桃国产精品网| 91麻豆精品激情在线观看国产| 精华霜和精华液先用哪个| 此物有八面人人有两片| 两个人的视频大全免费| 欧美精品亚洲一区二区| 无限看片的www在线观看| 久久久久久久久中文| 级片在线观看| 成熟少妇高潮喷水视频| 久久久精品国产亚洲av高清涩受| 好男人在线观看高清免费视频| 丝袜人妻中文字幕| 99久久无色码亚洲精品果冻| 国产亚洲精品久久久久5区| 日韩欧美免费精品| 国产一区二区三区在线臀色熟女| 床上黄色一级片| 亚洲熟妇中文字幕五十中出| 久久中文看片网| 他把我摸到了高潮在线观看| av欧美777| 国产精品久久久av美女十八| www日本黄色视频网| 午夜福利高清视频| 日本黄大片高清| 亚洲精品在线美女| 成人特级黄色片久久久久久久| 一级黄色大片毛片| av福利片在线观看| 午夜激情av网站| 精品福利观看| 亚洲国产日韩欧美精品在线观看 | 色综合亚洲欧美另类图片| 丰满人妻一区二区三区视频av | 日日干狠狠操夜夜爽| 久久中文字幕一级| videosex国产| 19禁男女啪啪无遮挡网站| 亚洲,欧美精品.| 一本一本综合久久| 91成年电影在线观看| 久久久久久大精品| 久久久精品欧美日韩精品| 久久性视频一级片| 精品日产1卡2卡| 国内精品久久久久精免费| 久久精品国产综合久久久| 深夜精品福利| 欧美日本亚洲视频在线播放| 午夜视频精品福利| 大型黄色视频在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 99久久久亚洲精品蜜臀av| 亚洲人成网站在线播放欧美日韩| 亚洲人成网站在线播放欧美日韩| 成人一区二区视频在线观看| 国产1区2区3区精品| 一本大道久久a久久精品| 国产精品99久久99久久久不卡| 麻豆成人av在线观看| 久久久国产欧美日韩av| 午夜亚洲福利在线播放| 国产单亲对白刺激| 日韩国内少妇激情av| 一区二区三区高清视频在线| 宅男免费午夜| 亚洲狠狠婷婷综合久久图片| cao死你这个sao货| 又大又爽又粗| 国模一区二区三区四区视频 | 久久亚洲精品不卡| 午夜a级毛片| 国产一区二区三区在线臀色熟女| 亚洲人与动物交配视频| 在线观看日韩欧美| 午夜激情福利司机影院| 国产97色在线日韩免费| 91麻豆av在线| 国内久久婷婷六月综合欲色啪| 91大片在线观看| 亚洲在线自拍视频| 国产激情偷乱视频一区二区| 在线观看舔阴道视频| 国产精品99久久99久久久不卡| 日韩 欧美 亚洲 中文字幕| 嫩草影院精品99| 成人av在线播放网站| 99精品在免费线老司机午夜| 精品少妇一区二区三区视频日本电影| 一区二区三区高清视频在线| 国产av不卡久久| 国产精品美女特级片免费视频播放器 | 国产亚洲精品久久久久5区| а√天堂www在线а√下载| 后天国语完整版免费观看| 成年女人毛片免费观看观看9| 一级毛片女人18水好多| 国产精品永久免费网站| 国产av一区在线观看免费| 可以在线观看毛片的网站| 亚洲va日本ⅴa欧美va伊人久久| 日本熟妇午夜| 精品不卡国产一区二区三区| 欧美在线黄色| 精品国产美女av久久久久小说| 欧美zozozo另类| 日韩av在线大香蕉| 午夜福利视频1000在线观看| 久久草成人影院| 欧美另类亚洲清纯唯美| 91在线观看av| 18禁黄网站禁片免费观看直播| 国产精品av久久久久免费| 国产亚洲av高清不卡| 757午夜福利合集在线观看| 成年版毛片免费区| 久久精品91蜜桃| 久久国产乱子伦精品免费另类| 亚洲熟妇中文字幕五十中出| 午夜免费激情av| 中文字幕精品亚洲无线码一区| 黑人欧美特级aaaaaa片| 不卡av一区二区三区| av在线播放免费不卡| 叶爱在线成人免费视频播放| 欧美黑人巨大hd| 三级毛片av免费| 欧美乱色亚洲激情| av国产免费在线观看| 一a级毛片在线观看| 国产精品久久久人人做人人爽| 国产精品久久久久久精品电影| 国内毛片毛片毛片毛片毛片| 丝袜人妻中文字幕| 一级作爱视频免费观看| 最近最新中文字幕大全电影3| 看黄色毛片网站| 在线视频色国产色| 18禁美女被吸乳视频| 精品久久蜜臀av无| 婷婷亚洲欧美| 国产av一区二区精品久久| 一个人观看的视频www高清免费观看 | 欧美日本亚洲视频在线播放| 叶爱在线成人免费视频播放| bbb黄色大片| 午夜日韩欧美国产| 美女大奶头视频| 日韩av在线大香蕉| 国内久久婷婷六月综合欲色啪| 精品不卡国产一区二区三区| 欧美黑人欧美精品刺激| 久99久视频精品免费| 黑人操中国人逼视频| 亚洲欧美日韩无卡精品| cao死你这个sao货| 正在播放国产对白刺激| 亚洲精品中文字幕在线视频| 黄色视频,在线免费观看| 老司机午夜福利在线观看视频| 婷婷精品国产亚洲av| www.999成人在线观看| 狂野欧美白嫩少妇大欣赏| 搡老熟女国产l中国老女人| 99国产极品粉嫩在线观看| 五月玫瑰六月丁香| 国产熟女午夜一区二区三区| 久久久久国产一级毛片高清牌| 可以在线观看毛片的网站| 中文亚洲av片在线观看爽| 国产精品一区二区三区四区久久| 成人手机av| 日日摸夜夜添夜夜添小说| 亚洲真实伦在线观看| 国产成人精品久久二区二区91| 亚洲欧洲精品一区二区精品久久久| 欧美日韩黄片免| 国产精品日韩av在线免费观看| 桃色一区二区三区在线观看| 久久久久国产精品人妻aⅴ院| 久久精品国产清高在天天线| 国产欧美日韩一区二区精品| 丁香六月欧美| 最好的美女福利视频网| 人成视频在线观看免费观看| 黄色a级毛片大全视频| 成人国产一区最新在线观看| 免费在线观看完整版高清| 首页视频小说图片口味搜索| 18禁观看日本| 真人一进一出gif抽搐免费| 国产精品98久久久久久宅男小说| 岛国在线免费视频观看| 99国产精品一区二区三区| 九九热线精品视视频播放| 50天的宝宝边吃奶边哭怎么回事| 日韩欧美三级三区| 久久人妻av系列| 国内揄拍国产精品人妻在线| 精品日产1卡2卡| 人人妻人人澡欧美一区二区| aaaaa片日本免费| 少妇粗大呻吟视频| 香蕉国产在线看| 亚洲专区国产一区二区| 岛国在线免费视频观看| 少妇熟女aⅴ在线视频| 国产精品精品国产色婷婷| 国产av麻豆久久久久久久| 在线免费观看的www视频| 99久久综合精品五月天人人| 一级a爱片免费观看的视频| 国内精品久久久久久久电影| 国产精品综合久久久久久久免费| 别揉我奶头~嗯~啊~动态视频| 亚洲中文字幕日韩| 日本一二三区视频观看| 小说图片视频综合网站| 久久久久国内视频| 久久香蕉精品热| 亚洲五月天丁香| 久久久久国产精品人妻aⅴ院| 亚洲九九香蕉| 特级一级黄色大片| 午夜亚洲福利在线播放| 欧美乱码精品一区二区三区| 欧美日韩亚洲综合一区二区三区_| 99国产精品99久久久久| 亚洲无线在线观看| 99久久99久久久精品蜜桃| 女生性感内裤真人,穿戴方法视频| 亚洲人成网站高清观看| 中文字幕人妻丝袜一区二区| 亚洲人成网站高清观看| 免费在线观看完整版高清| 欧美成人免费av一区二区三区| 欧美一区二区精品小视频在线| 极品教师在线免费播放| 在线观看免费午夜福利视频| 日韩欧美国产在线观看| 欧美国产日韩亚洲一区| 一本综合久久免费| 国产一级毛片七仙女欲春2| 久久人人精品亚洲av| 在线观看一区二区三区| 狠狠狠狠99中文字幕| 可以在线观看的亚洲视频| 欧美激情久久久久久爽电影| 美女午夜性视频免费| 午夜福利视频1000在线观看| 午夜久久久久精精品| 麻豆久久精品国产亚洲av| 最近视频中文字幕2019在线8| 最新在线观看一区二区三区| netflix在线观看网站| 黑人操中国人逼视频| 亚洲中文字幕一区二区三区有码在线看 | 亚洲一区二区三区不卡视频| 免费无遮挡裸体视频| 日韩免费av在线播放| 国产av不卡久久| 97超级碰碰碰精品色视频在线观看| 久久婷婷成人综合色麻豆| 狠狠狠狠99中文字幕| 精品国内亚洲2022精品成人| 亚洲国产看品久久| 久久久久九九精品影院| 婷婷丁香在线五月| 精品国产亚洲在线| 免费搜索国产男女视频| 亚洲av成人一区二区三| 国产精品一及| 久久久久久亚洲精品国产蜜桃av| 亚洲精品在线美女| 亚洲欧美一区二区三区黑人| 欧洲精品卡2卡3卡4卡5卡区| 国产成人精品久久二区二区91| 久久国产精品人妻蜜桃| 国产一级毛片七仙女欲春2| 国产精品爽爽va在线观看网站| 亚洲精品国产一区二区精华液| 每晚都被弄得嗷嗷叫到高潮| 久久天躁狠狠躁夜夜2o2o| 老司机午夜福利在线观看视频| 黄色 视频免费看| 免费电影在线观看免费观看| 欧美中文综合在线视频| 18禁黄网站禁片免费观看直播| 操出白浆在线播放| 在线a可以看的网站| 制服丝袜大香蕉在线| 亚洲av中文字字幕乱码综合| 岛国视频午夜一区免费看| 老司机午夜福利在线观看视频| ponron亚洲| 这个男人来自地球电影免费观看| 久久精品亚洲精品国产色婷小说| 天堂√8在线中文| 午夜成年电影在线免费观看| 最好的美女福利视频网| 亚洲欧美日韩高清专用| 黄色成人免费大全| 午夜亚洲福利在线播放| 免费一级毛片在线播放高清视频| 精品国产美女av久久久久小说| xxxwww97欧美| 两个人的视频大全免费| 久久精品91无色码中文字幕| 国产成+人综合+亚洲专区| 黄色视频不卡| 最新美女视频免费是黄的| 欧美在线黄色| 在线看三级毛片| 每晚都被弄得嗷嗷叫到高潮| 免费看十八禁软件| 欧美精品啪啪一区二区三区| 国产精品爽爽va在线观看网站| 91老司机精品| 黄色女人牲交| 成人三级做爰电影| 美女大奶头视频| 久久这里只有精品19| 妹子高潮喷水视频| 国产精品 国内视频| 欧美乱色亚洲激情| 亚洲无线在线观看| 亚洲欧美日韩高清在线视频| 国产免费av片在线观看野外av| 免费在线观看成人毛片| 每晚都被弄得嗷嗷叫到高潮| 99久久无色码亚洲精品果冻| 黑人操中国人逼视频| 国产精品 欧美亚洲| 在线永久观看黄色视频| 日韩成人在线观看一区二区三区| 1024手机看黄色片| 欧美激情久久久久久爽电影| 90打野战视频偷拍视频| √禁漫天堂资源中文www| 欧美av亚洲av综合av国产av| 黄色片一级片一级黄色片| 亚洲色图av天堂| 亚洲性夜色夜夜综合| 国产人伦9x9x在线观看| 日韩精品青青久久久久久| 男人的好看免费观看在线视频 | 男人舔女人的私密视频| 两个人视频免费观看高清| 女人被狂操c到高潮| 一边摸一边做爽爽视频免费| 男女视频在线观看网站免费 | www日本黄色视频网| 在线观看美女被高潮喷水网站 | 在线观看免费午夜福利视频| 国产激情欧美一区二区| 亚洲国产精品999在线| 久久午夜亚洲精品久久| 变态另类成人亚洲欧美熟女| 99热只有精品国产| bbb黄色大片| 不卡av一区二区三区| 亚洲国产精品999在线| 成人午夜高清在线视频| 国产片内射在线| 国产一级毛片七仙女欲春2| 午夜亚洲福利在线播放| 99热这里只有精品一区 | 99热这里只有精品一区 | √禁漫天堂资源中文www| 亚洲精品中文字幕在线视频| av有码第一页| 可以在线观看毛片的网站| 哪里可以看免费的av片| 九色国产91popny在线| 国产精品久久久人人做人人爽| 成人三级做爰电影| 一区二区三区高清视频在线| 最近最新中文字幕大全免费视频| 欧美高清成人免费视频www| 天堂动漫精品| 亚洲国产中文字幕在线视频| 亚洲 国产 在线| 桃红色精品国产亚洲av| 欧美日本视频| 国产一区二区在线观看日韩 | 成人av一区二区三区在线看| 91成年电影在线观看| 亚洲男人的天堂狠狠| 男人舔女人的私密视频| 亚洲七黄色美女视频| 日韩欧美在线乱码| 一二三四社区在线视频社区8| 欧美日韩亚洲国产一区二区在线观看| 天天躁夜夜躁狠狠躁躁| 18禁黄网站禁片午夜丰满| 成人国产一区最新在线观看| 美女 人体艺术 gogo| 亚洲精品在线观看二区| e午夜精品久久久久久久| 久久久国产欧美日韩av| 国产成人欧美在线观看| 97人妻精品一区二区三区麻豆| 久久久久国产一级毛片高清牌| 欧美在线黄色| 午夜成年电影在线免费观看| 国产免费av片在线观看野外av| 国产片内射在线| 久久中文看片网| 夜夜躁狠狠躁天天躁| 哪里可以看免费的av片| 真人做人爱边吃奶动态| 久久中文字幕人妻熟女| 亚洲一区二区三区不卡视频| 久久久久久久久久黄片| 久久伊人香网站| 丰满的人妻完整版| 丰满人妻熟妇乱又伦精品不卡| 久久精品aⅴ一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 国产一区二区在线观看日韩 | 亚洲美女黄片视频| 性欧美人与动物交配| 国产精品久久久av美女十八| 老司机福利观看| 老司机在亚洲福利影院| 宅男免费午夜| 亚洲第一欧美日韩一区二区三区| 2021天堂中文幕一二区在线观| www国产在线视频色| 2021天堂中文幕一二区在线观| 又爽又黄无遮挡网站| www日本黄色视频网| 国产激情久久老熟女| 国产av一区二区精品久久| 级片在线观看| 国产精品一区二区三区四区久久| 成人18禁高潮啪啪吃奶动态图| 国产精品一区二区三区四区久久| 亚洲 欧美一区二区三区| 女生性感内裤真人,穿戴方法视频| 人妻丰满熟妇av一区二区三区| 美女高潮喷水抽搐中文字幕| 欧美久久黑人一区二区| 天堂动漫精品| 国产激情偷乱视频一区二区| 嫩草影院精品99| 无限看片的www在线观看| 成人午夜高清在线视频| 国产精品99久久99久久久不卡| 久久天躁狠狠躁夜夜2o2o| 成人特级黄色片久久久久久久| 午夜精品在线福利| 日本免费一区二区三区高清不卡| 精品一区二区三区四区五区乱码| 国产野战对白在线观看| 欧美中文日本在线观看视频| 美女大奶头视频| 法律面前人人平等表现在哪些方面| 欧美黑人精品巨大| 国产片内射在线| 岛国视频午夜一区免费看| 精品国产乱子伦一区二区三区| 国产欧美日韩精品亚洲av| 日韩大尺度精品在线看网址| 黄色视频不卡| 精品久久久久久久人妻蜜臀av| 91麻豆精品激情在线观看国产| 夜夜夜夜夜久久久久| 久久婷婷人人爽人人干人人爱| 国内精品久久久久精免费| 国产黄a三级三级三级人| 国产精品野战在线观看| 免费人成视频x8x8入口观看| 大型黄色视频在线免费观看| 欧美国产日韩亚洲一区| 久久中文字幕人妻熟女| 麻豆久久精品国产亚洲av| 亚洲国产中文字幕在线视频| 九九热线精品视视频播放| 三级男女做爰猛烈吃奶摸视频| 亚洲免费av在线视频| 美女高潮喷水抽搐中文字幕| 不卡av一区二区三区| 日韩有码中文字幕| 久久99热这里只有精品18| 国产亚洲精品综合一区在线观看 | 天天一区二区日本电影三级| 久久久国产成人免费| 男女下面进入的视频免费午夜| 国产激情偷乱视频一区二区| 日本 欧美在线| 99国产极品粉嫩在线观看| 国产aⅴ精品一区二区三区波| 精品国产乱码久久久久久男人| 12—13女人毛片做爰片一| 亚洲国产中文字幕在线视频| 欧美性猛交╳xxx乱大交人| 亚洲 欧美 日韩 在线 免费| www国产在线视频色| 国产三级黄色录像| 在线观看66精品国产| 久久精品国产99精品国产亚洲性色| 天堂动漫精品| 欧美日韩精品网址| 在线观看免费日韩欧美大片| ponron亚洲| 欧美乱妇无乱码| 两个人看的免费小视频| 十八禁网站免费在线| 午夜日韩欧美国产| 午夜亚洲福利在线播放| 婷婷精品国产亚洲av| 一个人观看的视频www高清免费观看 | 久久久国产精品麻豆| 国产午夜精品论理片| 性色av乱码一区二区三区2| 午夜激情av网站| 少妇人妻一区二区三区视频| 一本一本综合久久| 男人舔女人的私密视频|