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

    某支線飛機(jī)概念設(shè)計(jì)階段機(jī)翼氣動(dòng)彈性設(shè)計(jì)

    2016-10-09 08:32:49張旭岳良明王斌
    航空工程進(jìn)展 2016年3期
    關(guān)鍵詞:氣動(dòng)彈性概念設(shè)計(jì)機(jī)翼

    張旭,岳良明,王斌

    (中國(guó)航天空氣動(dòng)力技術(shù)研究院 第二研究所,北京 100074)

    ?

    某支線飛機(jī)概念設(shè)計(jì)階段機(jī)翼氣動(dòng)彈性設(shè)計(jì)

    張旭,岳良明,王斌

    (中國(guó)航天空氣動(dòng)力技術(shù)研究院 第二研究所,北京100074)

    在飛機(jī)的總體方案概念設(shè)計(jì)階段,需要對(duì)飛機(jī)機(jī)翼進(jìn)行氣動(dòng)彈性設(shè)計(jì),以避免在后期設(shè)計(jì)中,因氣動(dòng)彈性問(wèn)題而對(duì)設(shè)計(jì)方案進(jìn)行較大更改。以某支線飛機(jī)的總體方案論證為背景,研究并歸納該飛機(jī)概念設(shè)計(jì)階段氣動(dòng)彈性設(shè)計(jì)與分析的理論基礎(chǔ)和計(jì)算方法,建立機(jī)翼的梁架式模型,初步設(shè)計(jì)模型剛度與質(zhì)量分布,并進(jìn)行機(jī)翼的靜氣動(dòng)彈性響應(yīng)與載荷分析、振動(dòng)特性及顫振特性分析。結(jié)果表明:該支線飛機(jī)機(jī)翼的氣動(dòng)彈性特性合理,符合設(shè)計(jì)要求。

    支線飛機(jī);總體方案設(shè)計(jì);氣動(dòng)彈性;靜氣動(dòng)彈性;飛行載荷;固有振動(dòng)特性;顫振

    0 引 言

    在飛機(jī)的大部分設(shè)計(jì)和研發(fā)階段,都需要考慮氣動(dòng)彈性和載荷的影響。為了避免在設(shè)計(jì)時(shí)因氣動(dòng)彈性問(wèn)題而作出較大更改,通常在一些總體參數(shù)確定之初,即飛機(jī)總體方案概念設(shè)計(jì)階段,對(duì)飛機(jī)尤其是機(jī)翼進(jìn)行氣動(dòng)彈性初步設(shè)計(jì),防止飛機(jī)氣動(dòng)彈性失穩(wěn)現(xiàn)象的發(fā)生[1]。

    國(guó)內(nèi)外對(duì)機(jī)翼的氣動(dòng)彈性設(shè)計(jì)進(jìn)行了大量研究,例如,Gareth A.Vio等[2]以氣動(dòng)彈性響應(yīng)為優(yōu)化目標(biāo),采用進(jìn)化算法對(duì)某典型復(fù)合材料機(jī)翼進(jìn)行了優(yōu)化設(shè)計(jì);Edward J.Alyanak[3]將多學(xué)科工具引入概念設(shè)計(jì)流程,在很大程度上提高了設(shè)計(jì)效率;霍應(yīng)元等[4]聯(lián)系工程實(shí)際對(duì)翼面剛度特性、機(jī)翼-發(fā)動(dòng)機(jī)短艙顫振特性等進(jìn)行了討論;蔡天星等[5]采用基于Euler方程和有限元的CFD/CSD耦合方法分析了大展弦比機(jī)翼的氣動(dòng)彈性時(shí)域響應(yīng),實(shí)現(xiàn)了靜、動(dòng)氣動(dòng)彈性一體化計(jì)算。目前,大部分設(shè)計(jì)手段均適用于初步設(shè)計(jì)及詳細(xì)設(shè)計(jì)階段,但針對(duì)概念設(shè)計(jì)階段的氣動(dòng)彈性設(shè)計(jì)并不多見(jiàn)。

    本文以某支線飛機(jī)的總體方案論證為背景,研究并歸納飛機(jī)概念設(shè)計(jì)階段氣動(dòng)彈性設(shè)計(jì)與分析的理論基礎(chǔ)和計(jì)算方法,初步設(shè)計(jì)機(jī)翼的剛度模型與質(zhì)量模型,分析并檢驗(yàn)其氣動(dòng)彈性特性。

    1 氣動(dòng)彈性模型

    1.1結(jié)構(gòu)模型

    飛機(jī)的基本數(shù)學(xué)模型必須能夠在整個(gè)關(guān)注的頻率范圍內(nèi)正確模擬其振動(dòng)特性。大型商用飛機(jī)的頻率范圍典型值為0~40 Hz,小型商用飛機(jī)則為0~60 Hz。因此,模型需要在上述頻率范圍內(nèi)產(chǎn)生固有頻率、模態(tài)質(zhì)量和正則模態(tài)形狀,且模型應(yīng)能夠充分模擬飛機(jī)的復(fù)雜性,包括操縱面、發(fā)動(dòng)機(jī)特性以及得到精確的模態(tài)形狀。

    1.1.1剛度模型

    對(duì)于具有細(xì)長(zhǎng)大展弦比機(jī)翼的飛機(jī),建立數(shù)學(xué)模型的傳統(tǒng)方法是在其擁有“梁式”結(jié)構(gòu)的基礎(chǔ)上,采用沿部件參考軸(例如彈性軸)布置的梁模型來(lái)模擬每個(gè)飛機(jī)部件(例如機(jī)翼、機(jī)身、平尾、垂尾等)。

    在細(xì)節(jié)結(jié)構(gòu)還未定義的概念設(shè)計(jì)階段,通過(guò)比例的方法沿用已有飛機(jī)的剛度值和質(zhì)量值,因此“梁式”剛度模型適用于飛機(jī)概念設(shè)計(jì)階段的氣動(dòng)彈性設(shè)計(jì)。

    1.1.2質(zhì)量模型

    為了進(jìn)行質(zhì)量分布的模擬,可以將機(jī)翼分成若干段(或片條),每段中心位于梁式模型的結(jié)構(gòu)參考點(diǎn)(節(jié)點(diǎn))。對(duì)于每一段機(jī)翼,質(zhì)量集中在參考位置上,通過(guò)剛性連接元素連接在梁軸節(jié)點(diǎn)上。這種形式的集中質(zhì)量可以表示各分段的質(zhì)量、慣性矩和質(zhì)量矩。剛性元素用來(lái)考慮各段質(zhì)心相對(duì)參考軸的偏離,從而定義與梁式剛度矩陣伴生的結(jié)構(gòu)質(zhì)量矩陣。

    1.2氣動(dòng)模型

    考慮氣動(dòng)彈性影響的場(chǎng)合一般需要采用非定常氣動(dòng)力模型來(lái)表示作用在彈性結(jié)構(gòu)上的力,氣動(dòng)力模型必須能夠計(jì)及氣動(dòng)力相對(duì)結(jié)構(gòu)運(yùn)動(dòng)的衰減和相位滯后,因此由偶極子格網(wǎng)法得到振蕩運(yùn)動(dòng)的氣動(dòng)力影響系數(shù)矩陣(AIC)是復(fù)數(shù)矩陣,且為減縮頻率的函數(shù),必須在一定范圍的減縮頻率和馬赫數(shù)下求取AIC矩陣。

    氣動(dòng)模型中包括與剛性模態(tài)、彈性模態(tài)相關(guān)的項(xiàng),即在零減縮頻率下可求得線性剛性飛機(jī)導(dǎo)數(shù)。減縮頻率下的氣動(dòng)矩陣應(yīng)用于顫振分析和突風(fēng)/湍流計(jì)算,而準(zhǔn)定常(零減縮頻率)氣動(dòng)力則應(yīng)用于靜氣彈和機(jī)動(dòng)分析。

    2 氣動(dòng)彈性分析基礎(chǔ)

    2.1振動(dòng)特性分析方程

    結(jié)構(gòu)固有振動(dòng)特性計(jì)算是在無(wú)阻尼結(jié)構(gòu)自由振動(dòng)的基礎(chǔ)上求解結(jié)構(gòu)振動(dòng)的頻率與模態(tài)。當(dāng)把結(jié)構(gòu)離散為有限元素后,其無(wú)阻尼自由振動(dòng)運(yùn)動(dòng)方程[6-7]為

    (1)

    式中:Ms為結(jié)構(gòu)整體質(zhì)量矩陣;Ks為結(jié)構(gòu)整體剛度矩陣;x為結(jié)點(diǎn)位移列陣。

    設(shè)x=feiωt,代入式(1)中,得

    (Ks-ω2Ms)f=0

    (2)

    要使式(2)中的f有非零解,則

    det(Ks-ω2Ms)=0

    (3)

    式(3)為系統(tǒng)的特征方程,求解該特征方程,得到ωi(固有振動(dòng)頻率)后,再代回式(2),即可求得特征向量fi(固有振型(模態(tài)))。通常,可根據(jù)所研究問(wèn)題的具體需要,只算出結(jié)構(gòu)的前若干階頻率與模態(tài)即可。

    2.2靜氣動(dòng)彈性響應(yīng)分析方程

    一般地,靜力學(xué)分析平衡方程可寫為

    (4)

    將氣動(dòng)力項(xiàng)引入式(4),去除單點(diǎn)約束、多點(diǎn)約束,并進(jìn)行Guyan減縮(靜態(tài)縮聚)后,可得靜氣動(dòng)彈性運(yùn)動(dòng)方程:

    (5)

    對(duì)式(5)求導(dǎo)并進(jìn)行相關(guān)計(jì)算, 可以求出彈性穩(wěn)定性、操縱性導(dǎo)數(shù)及相應(yīng)的配平參數(shù)值,此外, 還可得出變形、應(yīng)力和應(yīng)變。

    2.3顫振特性分析方程

    氣動(dòng)彈性運(yùn)動(dòng)方程[9]為

    (6)

    式中:q為廣義坐標(biāo)列陣;M為廣義質(zhì)量矩陣;K為廣義剛度矩陣;A為廣義氣動(dòng)力影響系數(shù)矩陣。

    其中,A的表達(dá)式[10]為

    (7)

    式中:FP為網(wǎng)格氣動(dòng)作用點(diǎn)處的模態(tài)矩陣;S=diag(ΔS1,…,ΔSn) 為面積加權(quán)陣,對(duì)角項(xiàng)為各氣動(dòng)網(wǎng)格的面積;D為非定常氣動(dòng)力影響系數(shù);FH為控制點(diǎn)的模態(tài)矩陣;k為減縮頻率;b為參考長(zhǎng)度。

    廣義非定常氣動(dòng)力矩陣(A)是關(guān)于馬赫數(shù)與減縮頻率k=ωb/V的復(fù)函數(shù)。顫振方程求解常用的方法包括V-g法、p-k法等,由于p-k法可以反映一定的亞臨界特性[11],本文采用該方法進(jìn)行顫振求解。

    3 機(jī)翼氣動(dòng)彈性設(shè)計(jì)

    3.1梁架式模型設(shè)計(jì)

    某支線飛機(jī)機(jī)翼展弦比為9.5,屬于大展弦比機(jī)翼,其特點(diǎn)是翼根效應(yīng)區(qū)相對(duì)較小,弦向變形可以忽略。

    綜合考慮該飛機(jī)總體方案概念設(shè)計(jì)階段的特點(diǎn),本文所提方案的氣動(dòng)彈性模型采用經(jīng)典的梁架式,根據(jù)初步確定的機(jī)翼氣動(dòng)外形、結(jié)構(gòu)布局形式,設(shè)計(jì)如下:

    (1) 主梁位置由翼肋剛心連線確定,自全機(jī)對(duì)稱面至翼梢共由23段梁元組成;

    (2) 在機(jī)翼主梁兩側(cè)沿展向布置剛度較大的無(wú)質(zhì)量桿,用來(lái)模擬翼肋和氣動(dòng)插值;

    (3) 沿機(jī)翼展向設(shè)計(jì)11個(gè)彈性桿,在彈性桿的兩端結(jié)點(diǎn)處布置設(shè)計(jì)的集中質(zhì)量;

    (4) 發(fā)動(dòng)機(jī)吊艙由剛性元模擬,掛架用彈性桿模擬,剛度可調(diào)。

    干凈機(jī)翼和帶發(fā)動(dòng)機(jī)機(jī)翼模型分別如圖1~圖2所示。

    圖1 干凈機(jī)翼梁架式模型Fig.1 Beam model of clean wing

    圖2 帶發(fā)動(dòng)機(jī)機(jī)翼梁架式模型Fig.2 Beam model of wing with engine

    3.2模型剛度設(shè)計(jì)

    根據(jù)設(shè)計(jì)經(jīng)驗(yàn)和大展弦比機(jī)翼主梁剛度沿展向呈指數(shù)分布的特點(diǎn),對(duì)某支線飛機(jī)的機(jī)翼主梁剛度進(jìn)行設(shè)計(jì)和多輪變參數(shù)分析,最終選擇氣動(dòng)彈性特性較優(yōu)的一組剛度(主梁各梁元的楊氏模量相同)。機(jī)翼主梁剛度沿展向的變化趨勢(shì)分別如圖3~圖5所示。

    圖3 機(jī)翼主梁垂直彎曲慣量沿展向變化曲線Fig.3 Vertical bending inertia of wing mainbeam vs. spanwise position

    圖4 機(jī)翼主梁水平彎曲慣量沿展向變化曲線Fig.4 Horizontal bending inertia of wing mainbeam vs. spanwise position

    圖5 機(jī)翼主梁扭轉(zhuǎn)慣量沿展向變化曲線Fig.5 Torsion inertia of wing main beam vs.spanwise position

    3.3模型質(zhì)量設(shè)計(jì)

    根據(jù)重量估算結(jié)果,機(jī)翼重量2 583kg,發(fā)動(dòng)機(jī)與掛架總重3 254kg,參考同類機(jī)型集中質(zhì)量布置,對(duì)該方案的機(jī)翼質(zhì)量進(jìn)行設(shè)計(jì)和變參數(shù)分析。

    3.4機(jī)翼靜氣動(dòng)彈性響應(yīng)與載荷分析

    根據(jù)設(shè)計(jì)要求,選取計(jì)算狀態(tài)為:飛行高度35 000ft,Ma=0.78,動(dòng)壓9 221.5Pa,攻角3°,計(jì)算結(jié)果如表 1所示。

    表1 機(jī)翼靜氣動(dòng)彈性響應(yīng)與載荷分析結(jié)果

    3.5干凈機(jī)翼振動(dòng)特性分析

    機(jī)翼根部采用固支約束,得出該飛機(jī)機(jī)翼的固有振動(dòng)特性如表2所示,由于篇幅所限,表2僅列出了前5階振動(dòng)模態(tài)的固有頻率,且未給出振型圖。

    表2 干凈機(jī)翼振動(dòng)特性

    3.6干凈機(jī)翼顫振特性分析

    在該機(jī)翼顫振計(jì)算中,非定常氣動(dòng)力計(jì)算采用亞聲速偶極子格網(wǎng)法,顫振求解采用p-k法。顫振分析V-g和V-f曲線分別如圖6~圖7所示。

    圖6 干凈機(jī)翼顫振分析V-g曲線Fig.6 V-g curve of clean wing

    計(jì)算所得機(jī)翼的顫振速度Vf=296m/s,顫振頻率ωf=11.4Hz,顫振型式為機(jī)翼三階彎曲與一階扭轉(zhuǎn)耦合。在飛行包線內(nèi),該干凈機(jī)翼不會(huì)發(fā)生顫振失穩(wěn)現(xiàn)象。

    3.7帶發(fā)動(dòng)機(jī)機(jī)翼振動(dòng)特性分析

    機(jī)翼根部采用固支約束,得出該飛機(jī)帶發(fā)動(dòng)機(jī)機(jī)翼的固有振動(dòng)特性如表3所示。

    表3 帶發(fā)動(dòng)機(jī)機(jī)翼振動(dòng)特性

    從表2可以看出:帶發(fā)動(dòng)機(jī)機(jī)翼的一階彎曲、三階彎曲頻率無(wú)明顯變化;但發(fā)動(dòng)機(jī)安裝使得機(jī)翼二階彎曲頻率降低,而一階扭轉(zhuǎn)頻率明顯提高。

    3.8帶發(fā)動(dòng)機(jī)機(jī)翼顫振特性分析

    采用3.6節(jié)所述計(jì)算方法,顫振分析V-g和V-f曲線分別如圖8~圖9所示。

    圖8 帶發(fā)動(dòng)機(jī)機(jī)翼顫振分析V-g曲線Fig.8 V-g curve of wing with engine

    圖9 帶發(fā)動(dòng)機(jī)機(jī)翼顫振分析V-f曲線Fig.9 V-f curve of wing with engine

    計(jì)算所得帶發(fā)動(dòng)機(jī)機(jī)翼的顫振速度Vf=318m/s,顫振頻率ωf=13.48Hz,顫振型式為機(jī)翼三階彎曲與一階扭轉(zhuǎn)耦合。同樣,該機(jī)翼在飛行包線內(nèi)不會(huì)發(fā)生顫振。

    4 結(jié) 論

    (1) 在概念設(shè)計(jì)方案中,本文所研究的某支線飛機(jī)巡航時(shí)翼尖最大變形為0.65m,約為半翼展的4.3%,與同級(jí)別機(jī)型相仿,該飛機(jī)的機(jī)翼剛度設(shè)計(jì)符合要求。

    (2) 機(jī)翼顫振特性分析結(jié)果表明,該設(shè)計(jì)方案中的飛機(jī)機(jī)翼在飛行包線內(nèi)不會(huì)發(fā)生顫振,進(jìn)一步說(shuō)明該支線飛機(jī)的機(jī)翼剛度設(shè)計(jì)、質(zhì)量分布等是合理的。

    (3) 發(fā)動(dòng)機(jī)的安裝在一定程度上提高了機(jī)翼的扭轉(zhuǎn)剛度,從而提高了扭轉(zhuǎn)頻率,減弱了與機(jī)翼三階彎曲的慣性耦合,提高了顫振速度。

    [1] 陳桂彬, 楊超, 鄒叢青. 氣動(dòng)彈性設(shè)計(jì)基礎(chǔ)[M]. 北京: 北京航空航天大學(xué)出版社, 2004.

    ChenGuibin,YangChao,ZouCongqing.Elementsofaeroelasticdynamics[M].Beijing:BeihangUniversityPress, 2004. (inChinese)

    [2]GarethAVio,GeorgiaGeorgiou,JonathanECooper.Designofcompositestructurestoimprovetheaeroelasticperformance[C].AIAA-2012-1458, 2012.

    [3]EdwardJAlyanak.Modelingforconceptualdesign:anaeroelasticapproach[C].AIAA-2012-1425, 2012.

    [4] 霍應(yīng)元, 譚申剛, 郭潤(rùn)江. 大型飛機(jī)氣動(dòng)彈性設(shè)計(jì)的工程考慮[C]. 昆明:第十一屆全國(guó)空氣彈性學(xué)術(shù)交流會(huì), 2009.

    HuoYingyuan,TanShengang,GuoRunjiang.Engineeringconsiderationsforaeroelasticdesignoflargeaircraft[C].Kunming:The11thNationalSymposiumonAeroelasticity, 2009.(inChinese)

    [5] 蔡天星, 徐敏, 郭敏杰, 等. 基于CFD/CSD的大展弦比機(jī)翼氣動(dòng)彈性研究[J]. 強(qiáng)度與環(huán)境, 2011, 38(3): 19-23.

    CaiTianxing,XuMin,GuoMinjie,etal.Aeroelasticanalysisofhigh-aspectratiowingbasedonCFD/CSD[J].Structure&EnvironmentEngineering, 2011, 38(3): 19-23.(inChinese)

    [6] 王勖成, 紹敏. 有限單元法基本原理和數(shù)值方法[M]. 北京: 清華大學(xué)出版社, 1997: 278-282.

    WangXucheng,ShaoMin.Basicprinciplesandnumericalmethodsoffiniteelementtechnique[M].Beijing:TsinghuaUniversityPress, 1997: 278-282.(inChinese)

    [7] 王尚文, 余旭東. 飛行器結(jié)構(gòu)動(dòng)力學(xué)[M]. 西安: 西北工業(yè)大學(xué)出版社, 1986: 97-112.

    WangShangwen,YuXudong.Structuraldynamicsofaircraft[M].Xi’an:NorthwesternPolytechnicalUniversityPress, 1986: 97-112.(inChinese)

    [8]GuyanRJ.Reductionofstiffnessandmassmatrices[J].AIAAJournal, 1965, 3(2): 380.

    [9]WrightJR,CooperJE. 飛機(jī)氣動(dòng)彈性力學(xué)及載荷導(dǎo)論[M]. 上海: 上海交通大學(xué)出版社, 2010.

    WrightJR,CooperJE.Introductiontoaircraftaeroelasticityandloads[M].Shanghai:ShanghaiJiaoTongUniversityPress, 2010.(inChinese)

    [10] 呂斌, 吳志剛, 楊超.T型尾翼顫振特性分析方法[J]. 工程力學(xué), 2008, 25(2): 230-234.

    LüBin,WuZhigang,YangChao.AnalysisofT-tailflutter[J].EngineeringMechanics, 2008, 25(2): 230-234.(inChinese)

    [11] 管德. 飛機(jī)氣動(dòng)彈性力學(xué)手冊(cè)[M]. 北京: 航空工業(yè)出版社, 1994.

    GuanDe.Manualofaircraftaeroelasticdynamics[M].Beijing:AviationIndustryPress, 1994.(inChinese)

    (編輯:馬文靜)

    Aeroelastic Design for Wing of a Regional Jet in Conceptual Design Phase

    Zhang Xu, Yue Liangming, Wang Bin

    (The Second Research Institute, China Academy of Aerospace Aerodynamics, Beijing 100074, China)

    Preventing from considerable modifications for design afterwards, aeroelastic design of wing should be considered in conceptual design of airplanes. With the background of conceptual design for a regional jet, theoretical bases and calculation methods of aeroelastic design in aircraft conceptual design are summarized, and a beam-frame model of wing is developed. Based on the model, the stiffness and mass distribution is designed. Then an analysis of static aeroelastic responses, flight loads, vibration characteristics, and flutter characteristics is carried out. The results indicate that the aeroelastic design of the regional jet is well satisfied with the requirements.

    regional jet; conceptual design; aeroelasticity; static aeroelasticity; flight loads; natural vibration characteristics; flutter

    2016-03-03;

    2016-03-28

    岳良明,250819185@qq.com

    1674-8190(2016)03-349-06

    V215.3

    A

    10.16615/j.cnki.1674-8190.2016.03.013

    張旭(1985-),男,工程師。主要研究方向:飛機(jī)總體設(shè)計(jì)、氣動(dòng)彈性分析與設(shè)計(jì)。

    岳良明(1982-),男,工程師。主要研究方向:飛機(jī)總體、氣動(dòng)設(shè)計(jì)。

    王斌(1980-),男,高級(jí)工程師。主要研究方向:飛機(jī)總體設(shè)計(jì)、氣動(dòng)優(yōu)化。

    猜你喜歡
    氣動(dòng)彈性概念設(shè)計(jì)機(jī)翼
    淺析概念設(shè)計(jì)在建筑結(jié)構(gòu)設(shè)計(jì)中的應(yīng)用
    概念設(shè)計(jì)在建筑結(jié)構(gòu)設(shè)計(jì)中的應(yīng)用論述
    變時(shí)滯間隙非線性機(jī)翼顫振主動(dòng)控制方法
    孫薈CG概念設(shè)計(jì)作品
    2016紅點(diǎn)獎(jiǎng)最佳概念設(shè)計(jì)TOP10
    飛翼無(wú)人機(jī)嗡鳴氣動(dòng)彈性響應(yīng)分析
    模態(tài)選取對(duì)靜氣動(dòng)彈性分析的影響
    機(jī)翼跨聲速抖振研究進(jìn)展
    直升機(jī)的氣動(dòng)彈性問(wèn)題
    大型風(fēng)力機(jī)整機(jī)氣動(dòng)彈性響應(yīng)計(jì)算
    好看av亚洲va欧美ⅴa在| 人人妻,人人澡人人爽秒播| 欧美色欧美亚洲另类二区| 变态另类丝袜制服| 国产精品99久久99久久久不卡| 九色国产91popny在线| 免费在线观看视频国产中文字幕亚洲| 国产亚洲欧美精品永久| 满18在线观看网站| 国产欧美日韩精品亚洲av| 日本精品一区二区三区蜜桃| 成年女人毛片免费观看观看9| 麻豆av在线久日| 听说在线观看完整版免费高清| 亚洲片人在线观看| 国产蜜桃级精品一区二区三区| 成人特级黄色片久久久久久久| 99久久99久久久精品蜜桃| av视频在线观看入口| 久久久久久久久久黄片| 精华霜和精华液先用哪个| 亚洲三区欧美一区| e午夜精品久久久久久久| 视频区欧美日本亚洲| 精品国产超薄肉色丝袜足j| 可以在线观看的亚洲视频| 视频在线观看一区二区三区| 久久久水蜜桃国产精品网| 国产男靠女视频免费网站| 国产伦在线观看视频一区| 国产亚洲欧美98| 久久亚洲真实| 国产一卡二卡三卡精品| 好男人电影高清在线观看| 欧美中文综合在线视频| 欧美黑人巨大hd| 午夜视频精品福利| 日韩国内少妇激情av| xxxwww97欧美| 久久九九热精品免费| 久久婷婷人人爽人人干人人爱| 欧美日本亚洲视频在线播放| 亚洲精品一卡2卡三卡4卡5卡| 最新在线观看一区二区三区| 香蕉国产在线看| 十八禁人妻一区二区| 国产成年人精品一区二区| 美女免费视频网站| 最好的美女福利视频网| 国产黄片美女视频| 99久久精品国产亚洲精品| 欧美乱妇无乱码| 在线天堂中文资源库| 美女午夜性视频免费| 一二三四在线观看免费中文在| 不卡av一区二区三区| 一区二区三区国产精品乱码| 国产乱人伦免费视频| 成人特级黄色片久久久久久久| 夜夜看夜夜爽夜夜摸| 一区二区三区国产精品乱码| 日韩大尺度精品在线看网址| 人妻丰满熟妇av一区二区三区| 观看免费一级毛片| 最新美女视频免费是黄的| 男女那种视频在线观看| 伦理电影免费视频| aaaaa片日本免费| 精品第一国产精品| 怎么达到女性高潮| 久久热在线av| 国产人伦9x9x在线观看| 日本一区二区免费在线视频| 国产精品美女特级片免费视频播放器 | 一区二区三区国产精品乱码| 亚洲avbb在线观看| aaaaa片日本免费| 国产成人一区二区三区免费视频网站| 午夜老司机福利片| 亚洲国产毛片av蜜桃av| 99久久99久久久精品蜜桃| 亚洲av中文字字幕乱码综合 | АⅤ资源中文在线天堂| 成人亚洲精品一区在线观看| 一边摸一边抽搐一进一小说| 国产极品粉嫩免费观看在线| 女人被狂操c到高潮| 亚洲avbb在线观看| 在线播放国产精品三级| 制服人妻中文乱码| 国产精品亚洲av一区麻豆| e午夜精品久久久久久久| 国产成人欧美在线观看| 亚洲欧美激情综合另类| 亚洲va日本ⅴa欧美va伊人久久| 一本精品99久久精品77| 国产av一区在线观看免费| 人妻久久中文字幕网| 桃色一区二区三区在线观看| 在线av久久热| e午夜精品久久久久久久| 婷婷丁香在线五月| 麻豆成人午夜福利视频| 精品熟女少妇八av免费久了| 欧美大码av| 黄色毛片三级朝国网站| 国产精品一区二区精品视频观看| 午夜两性在线视频| 亚洲自拍偷在线| 精品国产超薄肉色丝袜足j| 欧美av亚洲av综合av国产av| 免费在线观看完整版高清| 免费人成视频x8x8入口观看| 亚洲avbb在线观看| 又黄又爽又免费观看的视频| 少妇粗大呻吟视频| 免费在线观看亚洲国产| 欧美在线黄色| 国产精品综合久久久久久久免费| 久久久久精品国产欧美久久久| 国产精品一区二区免费欧美| 50天的宝宝边吃奶边哭怎么回事| 一区二区三区国产精品乱码| 中文资源天堂在线| 欧美性长视频在线观看| 动漫黄色视频在线观看| 波多野结衣高清无吗| 女人高潮潮喷娇喘18禁视频| 亚洲av成人不卡在线观看播放网| 日韩三级视频一区二区三区| 午夜亚洲福利在线播放| 十分钟在线观看高清视频www| 91国产中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 在线国产一区二区在线| av天堂在线播放| 天堂动漫精品| 1024手机看黄色片| 在线av久久热| 日日摸夜夜添夜夜添小说| 国产亚洲av高清不卡| 日韩欧美一区二区三区在线观看| 最近最新中文字幕大全免费视频| 久久性视频一级片| 成人亚洲精品av一区二区| 精品国产乱子伦一区二区三区| 99国产精品一区二区三区| 午夜福利高清视频| 久久久久久国产a免费观看| 麻豆国产av国片精品| 18美女黄网站色大片免费观看| 精品一区二区三区四区五区乱码| 亚洲第一欧美日韩一区二区三区| 国产1区2区3区精品| 不卡一级毛片| 国产高清视频在线播放一区| 久久精品影院6| 午夜免费激情av| 天堂影院成人在线观看| 久久婷婷成人综合色麻豆| 丰满的人妻完整版| 午夜福利免费观看在线| 在线观看一区二区三区| 搡老妇女老女人老熟妇| 美女扒开内裤让男人捅视频| 99国产精品一区二区三区| 亚洲国产精品999在线| 国产精品九九99| 又大又爽又粗| 日日摸夜夜添夜夜添小说| 国产午夜福利久久久久久| 午夜福利高清视频| 黄频高清免费视频| 成人亚洲精品av一区二区| av电影中文网址| 伊人久久大香线蕉亚洲五| 亚洲av美国av| 淫秽高清视频在线观看| 69av精品久久久久久| 国产av又大| 亚洲成av人片免费观看| 特大巨黑吊av在线直播 | 日韩欧美 国产精品| 黄片小视频在线播放| 午夜福利高清视频| 久久久久久九九精品二区国产 | 视频在线观看一区二区三区| av在线天堂中文字幕| 国产私拍福利视频在线观看| 久久婷婷成人综合色麻豆| 夜夜夜夜夜久久久久| 成人三级做爰电影| 黄色丝袜av网址大全| www日本黄色视频网| 欧美黑人精品巨大| 国产免费av片在线观看野外av| 国内揄拍国产精品人妻在线 | 欧美日韩黄片免| 国产野战对白在线观看| 精品福利观看| 日韩免费av在线播放| 久久婷婷人人爽人人干人人爱| 国产亚洲精品久久久久5区| 黄色 视频免费看| 国产成人精品久久二区二区免费| 叶爱在线成人免费视频播放| 淫妇啪啪啪对白视频| avwww免费| 在线观看免费视频日本深夜| 国产一卡二卡三卡精品| 国产精华一区二区三区| 亚洲精品av麻豆狂野| 亚洲欧美日韩无卡精品| 亚洲五月天丁香| 色综合欧美亚洲国产小说| 亚洲一区中文字幕在线| 一级a爱片免费观看的视频| a级毛片a级免费在线| 国产精品久久久久久亚洲av鲁大| 搡老岳熟女国产| 亚洲精品久久成人aⅴ小说| 免费一级毛片在线播放高清视频| 亚洲欧美精品综合久久99| x7x7x7水蜜桃| 成人免费观看视频高清| 免费女性裸体啪啪无遮挡网站| 99国产精品一区二区蜜桃av| 午夜免费激情av| 午夜成年电影在线免费观看| 中文字幕高清在线视频| 久久婷婷成人综合色麻豆| av欧美777| 色播亚洲综合网| 非洲黑人性xxxx精品又粗又长| 久久午夜综合久久蜜桃| 黑人操中国人逼视频| 我的亚洲天堂| 欧美乱色亚洲激情| www.自偷自拍.com| 欧美av亚洲av综合av国产av| 久久国产精品男人的天堂亚洲| 午夜视频精品福利| 一区二区三区高清视频在线| 美女大奶头视频| avwww免费| av欧美777| 老熟妇仑乱视频hdxx| 亚洲中文字幕日韩| 又黄又爽又免费观看的视频| 怎么达到女性高潮| 日日摸夜夜添夜夜添小说| 亚洲欧美激情综合另类| 极品教师在线免费播放| 精品国产美女av久久久久小说| 黄色a级毛片大全视频| 中出人妻视频一区二区| 欧美色欧美亚洲另类二区| 操出白浆在线播放| 欧美另类亚洲清纯唯美| 一级毛片精品| 欧美激情久久久久久爽电影| 亚洲精品国产区一区二| 啦啦啦免费观看视频1| 亚洲,欧美精品.| 大型av网站在线播放| 午夜成年电影在线免费观看| 在线观看www视频免费| 十八禁网站免费在线| 99在线视频只有这里精品首页| 日本一本二区三区精品| 韩国av一区二区三区四区| 热re99久久国产66热| 日本成人三级电影网站| 欧美av亚洲av综合av国产av| 色在线成人网| 色综合亚洲欧美另类图片| 国产亚洲av高清不卡| 长腿黑丝高跟| 99国产综合亚洲精品| 999久久久精品免费观看国产| 女性生殖器流出的白浆| 亚洲五月婷婷丁香| 国产激情偷乱视频一区二区| 欧美最黄视频在线播放免费| 中文亚洲av片在线观看爽| 少妇熟女aⅴ在线视频| 成人国产一区最新在线观看| 成人亚洲精品一区在线观看| 免费看十八禁软件| 波多野结衣高清无吗| 青草久久国产| 国产精品亚洲av一区麻豆| 精品久久久久久成人av| 欧美日本视频| 亚洲欧美精品综合一区二区三区| www.自偷自拍.com| 搡老岳熟女国产| 美女高潮到喷水免费观看| 欧美乱妇无乱码| 午夜成年电影在线免费观看| 国产精品亚洲一级av第二区| 久久久久久久久久黄片| 日本免费a在线| 熟女少妇亚洲综合色aaa.| 免费无遮挡裸体视频| 午夜影院日韩av| av在线天堂中文字幕| 国内久久婷婷六月综合欲色啪| 18禁观看日本| 欧美黑人巨大hd| 无限看片的www在线观看| 亚洲,欧美精品.| 国产成人一区二区三区免费视频网站| 欧美性长视频在线观看| 日韩欧美三级三区| 国产成人av教育| 欧美日韩瑟瑟在线播放| 国产一区二区在线av高清观看| 亚洲欧美精品综合久久99| 精品福利观看| 日韩视频一区二区在线观看| 久久国产精品人妻蜜桃| 老汉色∧v一级毛片| 在线永久观看黄色视频| cao死你这个sao货| 人人澡人人妻人| 精品国产亚洲在线| 亚洲aⅴ乱码一区二区在线播放 | 88av欧美| 女警被强在线播放| 老司机在亚洲福利影院| 999久久久国产精品视频| 免费观看人在逋| 成人18禁高潮啪啪吃奶动态图| 一进一出抽搐gif免费好疼| 久久中文字幕一级| 观看免费一级毛片| 精品国产超薄肉色丝袜足j| 国产av不卡久久| 亚洲国产中文字幕在线视频| 最好的美女福利视频网| 后天国语完整版免费观看| 热re99久久国产66热| 久久中文看片网| 成年免费大片在线观看| 好男人电影高清在线观看| 免费看十八禁软件| 亚洲一区二区三区不卡视频| 亚洲国产欧洲综合997久久, | 欧美绝顶高潮抽搐喷水| 日韩欧美国产一区二区入口| 国产成人av教育| 黄网站色视频无遮挡免费观看| 99国产综合亚洲精品| 国内精品久久久久精免费| 亚洲一区二区三区色噜噜| 中文字幕人妻丝袜一区二区| 成年免费大片在线观看| 夜夜看夜夜爽夜夜摸| 免费在线观看完整版高清| 亚洲一区高清亚洲精品| 久久精品人妻少妇| 麻豆av在线久日| 久久精品亚洲精品国产色婷小说| 巨乳人妻的诱惑在线观看| 在线观看免费日韩欧美大片| 757午夜福利合集在线观看| 中文字幕人成人乱码亚洲影| 欧美丝袜亚洲另类 | 国产野战对白在线观看| 欧美日韩中文字幕国产精品一区二区三区| 色av中文字幕| 日本成人三级电影网站| 国产黄色小视频在线观看| 成人一区二区视频在线观看| 国产精品香港三级国产av潘金莲| 十八禁人妻一区二区| 欧美黄色片欧美黄色片| 99热6这里只有精品| 一区福利在线观看| 9191精品国产免费久久| 亚洲 国产 在线| 亚洲精品美女久久久久99蜜臀| 亚洲色图 男人天堂 中文字幕| 色av中文字幕| 国语自产精品视频在线第100页| 国产精品综合久久久久久久免费| 日韩欧美三级三区| 又大又爽又粗| 免费在线观看日本一区| 91字幕亚洲| 精品国产乱码久久久久久男人| 亚洲人成电影免费在线| 亚洲国产中文字幕在线视频| 黄色a级毛片大全视频| 久久精品人妻少妇| 美女免费视频网站| 少妇裸体淫交视频免费看高清 | 1024视频免费在线观看| 亚洲 欧美 日韩 在线 免费| 久久天堂一区二区三区四区| 狠狠狠狠99中文字幕| 免费高清在线观看日韩| 男女床上黄色一级片免费看| 18禁国产床啪视频网站| 在线av久久热| 亚洲人成网站在线播放欧美日韩| 免费电影在线观看免费观看| 一卡2卡三卡四卡精品乱码亚洲| av天堂在线播放| 亚洲av日韩精品久久久久久密| 久久精品人妻少妇| 少妇被粗大的猛进出69影院| 好看av亚洲va欧美ⅴa在| 欧美黄色片欧美黄色片| 久久国产乱子伦精品免费另类| 男女床上黄色一级片免费看| 国产亚洲精品久久久久5区| 欧美在线一区亚洲| 欧美日韩亚洲国产一区二区在线观看| 亚洲av五月六月丁香网| 美女大奶头视频| 黄片大片在线免费观看| 国产激情久久老熟女| 中文字幕精品免费在线观看视频| 丝袜在线中文字幕| 成人免费观看视频高清| 免费在线观看亚洲国产| 女性被躁到高潮视频| 在线永久观看黄色视频| 99在线视频只有这里精品首页| 国产在线精品亚洲第一网站| 欧美三级亚洲精品| 黄色视频,在线免费观看| 在线观看日韩欧美| 欧美不卡视频在线免费观看 | 国产在线观看jvid| 久99久视频精品免费| 麻豆国产av国片精品| 我的亚洲天堂| 久久久久国产精品人妻aⅴ院| 亚洲aⅴ乱码一区二区在线播放 | 2021天堂中文幕一二区在线观 | 亚洲av美国av| 国产精品野战在线观看| 18禁美女被吸乳视频| 久热爱精品视频在线9| 久久久久九九精品影院| 韩国精品一区二区三区| 欧美一级毛片孕妇| 国产精品久久久久久精品电影 | 免费看美女性在线毛片视频| 老熟妇乱子伦视频在线观看| 18美女黄网站色大片免费观看| 亚洲人成网站在线播放欧美日韩| 久99久视频精品免费| 久久人人精品亚洲av| 亚洲avbb在线观看| 亚洲精品美女久久久久99蜜臀| 亚洲国产欧美网| 精品免费久久久久久久清纯| 变态另类成人亚洲欧美熟女| 一级作爱视频免费观看| 久久国产乱子伦精品免费另类| 国产精品久久电影中文字幕| 波多野结衣高清作品| 在线国产一区二区在线| 亚洲一码二码三码区别大吗| 亚洲电影在线观看av| 可以在线观看的亚洲视频| 色在线成人网| 国产人伦9x9x在线观看| 哪里可以看免费的av片| 一本大道久久a久久精品| 免费人成视频x8x8入口观看| 日本黄色视频三级网站网址| 免费在线观看亚洲国产| 国产av一区二区精品久久| 欧美最黄视频在线播放免费| 国内毛片毛片毛片毛片毛片| 一二三四在线观看免费中文在| 99在线视频只有这里精品首页| 国产野战对白在线观看| 国产久久久一区二区三区| 国产精品二区激情视频| 久久国产亚洲av麻豆专区| 日日干狠狠操夜夜爽| 女同久久另类99精品国产91| 非洲黑人性xxxx精品又粗又长| 午夜福利欧美成人| 欧美最黄视频在线播放免费| 亚洲精品中文字幕一二三四区| 亚洲三区欧美一区| av片东京热男人的天堂| 久热这里只有精品99| 成人一区二区视频在线观看| 日日干狠狠操夜夜爽| 动漫黄色视频在线观看| 最新在线观看一区二区三区| 亚洲成人国产一区在线观看| 99久久99久久久精品蜜桃| 日本 欧美在线| 国内少妇人妻偷人精品xxx网站 | 国产熟女午夜一区二区三区| 国产精品久久久av美女十八| 精品一区二区三区av网在线观看| 91av网站免费观看| 一边摸一边做爽爽视频免费| 精品一区二区三区四区五区乱码| 国产蜜桃级精品一区二区三区| 亚洲 欧美一区二区三区| 日本a在线网址| 国产熟女午夜一区二区三区| 麻豆成人av在线观看| 日本 欧美在线| 精品欧美一区二区三区在线| 草草在线视频免费看| 亚洲人成网站高清观看| 国产三级黄色录像| 好看av亚洲va欧美ⅴa在| 老鸭窝网址在线观看| 婷婷精品国产亚洲av| 欧美激情高清一区二区三区| 国产精品亚洲一级av第二区| 别揉我奶头~嗯~啊~动态视频| 日韩成人在线观看一区二区三区| 亚洲天堂国产精品一区在线| 欧美日韩中文字幕国产精品一区二区三区| 超碰成人久久| 久久午夜亚洲精品久久| 亚洲欧美日韩无卡精品| 两个人看的免费小视频| 极品教师在线免费播放| 国产精品98久久久久久宅男小说| av视频在线观看入口| 久久国产精品影院| 热99re8久久精品国产| 国产亚洲精品第一综合不卡| 欧美大码av| 国产精品爽爽va在线观看网站 | 久久精品国产亚洲av高清一级| 成人特级黄色片久久久久久久| 亚洲欧美激情综合另类| 日日干狠狠操夜夜爽| 一区二区三区激情视频| 夜夜夜夜夜久久久久| 香蕉av资源在线| 亚洲国产欧美网| 日韩大码丰满熟妇| av中文乱码字幕在线| 久久久久精品国产欧美久久久| 欧美日韩中文字幕国产精品一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 午夜激情福利司机影院| 黄色女人牲交| 色综合婷婷激情| 久久天躁狠狠躁夜夜2o2o| 嫩草影视91久久| 欧美国产日韩亚洲一区| 哪里可以看免费的av片| 成人国产综合亚洲| 可以在线观看的亚洲视频| 美女大奶头视频| 国产精品一区二区免费欧美| 制服丝袜大香蕉在线| 妹子高潮喷水视频| 国产v大片淫在线免费观看| 日韩欧美三级三区| 亚洲专区字幕在线| 久久午夜亚洲精品久久| 一区福利在线观看| 精品国产国语对白av| 国产成人影院久久av| 国产野战对白在线观看| 一本精品99久久精品77| 久久人妻福利社区极品人妻图片| 日本黄色视频三级网站网址| 少妇的丰满在线观看| 日本精品一区二区三区蜜桃| 天天一区二区日本电影三级| 亚洲人成电影免费在线| 日韩欧美一区视频在线观看| 波多野结衣巨乳人妻| 久久精品91无色码中文字幕| 色在线成人网| 人妻丰满熟妇av一区二区三区| 最近最新中文字幕大全免费视频| 一区二区三区高清视频在线| www.精华液| 精品久久久久久久毛片微露脸| 中文亚洲av片在线观看爽| 十分钟在线观看高清视频www| 亚洲成av片中文字幕在线观看| 白带黄色成豆腐渣| 精品无人区乱码1区二区| 1024香蕉在线观看| 亚洲自偷自拍图片 自拍| 国产熟女午夜一区二区三区| 国产蜜桃级精品一区二区三区| 中文资源天堂在线| 青草久久国产| 国产真人三级小视频在线观看| 久久 成人 亚洲| 欧美成人一区二区免费高清观看 | 在线观看www视频免费| 久久精品国产99精品国产亚洲性色| 亚洲av成人av| 久久久久国产精品人妻aⅴ院| 国产精品美女特级片免费视频播放器 | 此物有八面人人有两片| 一本精品99久久精品77| 国产高清有码在线观看视频 | 亚洲第一av免费看| 男男h啪啪无遮挡| 高清毛片免费观看视频网站|