• <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ì)算
    一级毛片 在线播放| 亚洲,一卡二卡三卡| 丝袜脚勾引网站| 下体分泌物呈黄色| 亚洲三级黄色毛片| 午夜福利在线观看免费完整高清在| 五月伊人婷婷丁香| 有码 亚洲区| 啦啦啦中文免费视频观看日本| 91久久精品国产一区二区成人| 久久 成人 亚洲| 在线免费观看不下载黄p国产| 成年av动漫网址| videos熟女内射| 日韩欧美一区视频在线观看| 91精品三级在线观看| 18在线观看网站| 国产精品一区二区在线观看99| 精品国产一区二区三区久久久樱花| 成人毛片a级毛片在线播放| 狠狠婷婷综合久久久久久88av| 人人妻人人澡人人看| 777米奇影视久久| 免费大片18禁| 亚洲综合精品二区| 乱码一卡2卡4卡精品| 午夜免费男女啪啪视频观看| 伦理电影大哥的女人| 中文字幕人妻丝袜制服| 亚洲精品乱码久久久久久按摩| 国产成人91sexporn| 免费av中文字幕在线| 黄色毛片三级朝国网站| 麻豆乱淫一区二区| 久久精品国产亚洲av天美| 日韩欧美一区视频在线观看| 久久人妻熟女aⅴ| 2022亚洲国产成人精品| 18禁裸乳无遮挡动漫免费视频| 久久久久人妻精品一区果冻| 在现免费观看毛片| 亚洲人成网站在线播| 欧美成人午夜免费资源| 人体艺术视频欧美日本| 91aial.com中文字幕在线观看| av一本久久久久| 国内精品宾馆在线| 久久99热6这里只有精品| 日本91视频免费播放| 欧美激情 高清一区二区三区| 蜜桃在线观看..| 久久免费观看电影| 99久久中文字幕三级久久日本| 伦理电影免费视频| 青春草国产在线视频| av黄色大香蕉| www.av在线官网国产| 国产 一区精品| 国产成人精品婷婷| 能在线免费看毛片的网站| 亚洲精品成人av观看孕妇| 精品一品国产午夜福利视频| 成人亚洲欧美一区二区av| 汤姆久久久久久久影院中文字幕| 啦啦啦在线观看免费高清www| 久久午夜福利片| 建设人人有责人人尽责人人享有的| 精品亚洲成a人片在线观看| 国产精品国产av在线观看| 免费av不卡在线播放| freevideosex欧美| 黑人猛操日本美女一级片| 能在线免费看毛片的网站| 免费看av在线观看网站| 另类亚洲欧美激情| 亚洲欧美一区二区三区黑人 | 一区二区三区免费毛片| 另类精品久久| 97超碰精品成人国产| 女性生殖器流出的白浆| 九草在线视频观看| 色94色欧美一区二区| 精品国产一区二区三区久久久樱花| 伦理电影大哥的女人| 午夜福利视频在线观看免费| 亚洲av不卡在线观看| 美女xxoo啪啪120秒动态图| 久久精品久久久久久噜噜老黄| 亚洲av中文av极速乱| 成人国产麻豆网| 日韩人妻高清精品专区| 日本-黄色视频高清免费观看| 国产探花极品一区二区| 久久av网站| 欧美精品高潮呻吟av久久| 人妻 亚洲 视频| 亚洲经典国产精华液单| 亚洲精品,欧美精品| 亚洲国产欧美在线一区| 男男h啪啪无遮挡| 九色亚洲精品在线播放| 久久久久久人妻| 亚洲精品一二三| 国产精品99久久久久久久久| 免费高清在线观看日韩| 91久久精品电影网| 最近最新中文字幕免费大全7| 99九九线精品视频在线观看视频| 秋霞在线观看毛片| 日韩精品免费视频一区二区三区 | av免费观看日本| 欧美最新免费一区二区三区| 91精品一卡2卡3卡4卡| 国产亚洲精品第一综合不卡 | 亚洲伊人久久精品综合| 亚洲精品自拍成人| 亚洲av日韩在线播放| 亚洲人成网站在线观看播放| 老女人水多毛片| 国产日韩欧美亚洲二区| 伦精品一区二区三区| 亚洲第一区二区三区不卡| 免费播放大片免费观看视频在线观看| 成人18禁高潮啪啪吃奶动态图 | av在线老鸭窝| 美女脱内裤让男人舔精品视频| 欧美亚洲 丝袜 人妻 在线| 国产男人的电影天堂91| 高清欧美精品videossex| 18禁动态无遮挡网站| 蜜桃在线观看..| 我的老师免费观看完整版| 亚洲伊人久久精品综合| 精品卡一卡二卡四卡免费| 久久精品久久久久久久性| 免费人成在线观看视频色| 高清午夜精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 国产高清有码在线观看视频| 人妻系列 视频| 精品久久久噜噜| 亚洲中文av在线| 欧美亚洲 丝袜 人妻 在线| 一级片'在线观看视频| 精品卡一卡二卡四卡免费| 婷婷成人精品国产| 日本黄色日本黄色录像| 2022亚洲国产成人精品| 国产综合精华液| 男人爽女人下面视频在线观看| 久久精品国产鲁丝片午夜精品| 久久鲁丝午夜福利片| 亚洲av中文av极速乱| 中文精品一卡2卡3卡4更新| 99九九在线精品视频| 亚洲精品456在线播放app| 考比视频在线观看| 亚洲av欧美aⅴ国产| 桃花免费在线播放| 欧美 亚洲 国产 日韩一| 久久精品国产a三级三级三级| 人妻一区二区av| 这个男人来自地球电影免费观看 | av电影中文网址| 三上悠亚av全集在线观看| 欧美激情 高清一区二区三区| 波野结衣二区三区在线| 欧美日韩视频精品一区| 免费大片18禁| 久久人人爽人人片av| 秋霞伦理黄片| 成人毛片a级毛片在线播放| 国产精品成人在线| a级毛片免费高清观看在线播放| 女性被躁到高潮视频| av在线老鸭窝| 亚洲第一区二区三区不卡| 97超视频在线观看视频| 国产男人的电影天堂91| 亚洲综合精品二区| 精品亚洲乱码少妇综合久久| 美女脱内裤让男人舔精品视频| 国产成人91sexporn| 成人影院久久| 啦啦啦啦在线视频资源| 精品少妇黑人巨大在线播放| 久久精品久久久久久噜噜老黄| 观看美女的网站| 秋霞在线观看毛片| 亚洲人成网站在线观看播放| 成人黄色视频免费在线看| 青青草视频在线视频观看| 麻豆精品久久久久久蜜桃| 热99久久久久精品小说推荐| 人成视频在线观看免费观看| 高清av免费在线| 曰老女人黄片| 少妇的逼水好多| 天天操日日干夜夜撸| 日韩中文字幕视频在线看片| 18禁在线播放成人免费| 亚洲精品456在线播放app| 婷婷成人精品国产| 国产精品人妻久久久久久| 日韩伦理黄色片| 春色校园在线视频观看| 久久久a久久爽久久v久久| 熟女电影av网| 97超视频在线观看视频| 水蜜桃什么品种好| 日韩一区二区视频免费看| av.在线天堂| 99久久中文字幕三级久久日本| 一区二区三区乱码不卡18| 涩涩av久久男人的天堂| 久久久国产精品麻豆| 菩萨蛮人人尽说江南好唐韦庄| 99久久综合免费| 免费黄网站久久成人精品| 国产白丝娇喘喷水9色精品| 美女主播在线视频| 99久久综合免费| 精品亚洲成a人片在线观看| 狠狠婷婷综合久久久久久88av| 亚洲一级一片aⅴ在线观看| 在线观看www视频免费| 欧美日韩一区二区视频在线观看视频在线| 夫妻午夜视频| 在线观看www视频免费| 国产成人精品一,二区| 国产成人精品在线电影| 亚洲国产精品专区欧美| 久久国内精品自在自线图片| 国产色婷婷99| 精品人妻偷拍中文字幕| 黄片播放在线免费| 亚洲精品国产av蜜桃| 亚洲国产日韩一区二区| av免费在线看不卡| 免费观看的影片在线观看| 久久狼人影院| 亚洲熟女精品中文字幕| 精品久久久精品久久久| 欧美成人精品欧美一级黄| 中国美白少妇内射xxxbb| 交换朋友夫妻互换小说| 99九九线精品视频在线观看视频| 嘟嘟电影网在线观看| 国产成人精品福利久久| 亚洲av.av天堂| 久久精品国产鲁丝片午夜精品| 亚洲成色77777| 一级毛片 在线播放| 99久国产av精品国产电影| 青春草视频在线免费观看| 纵有疾风起免费观看全集完整版| 男女高潮啪啪啪动态图| 999精品在线视频| 精品酒店卫生间| 色吧在线观看| 午夜激情av网站| 色5月婷婷丁香| 欧美日韩精品成人综合77777| 少妇猛男粗大的猛烈进出视频| 两个人免费观看高清视频| 免费高清在线观看视频在线观看| 最后的刺客免费高清国语| 国产日韩欧美在线精品| av在线观看视频网站免费| 一级二级三级毛片免费看| 国产69精品久久久久777片| 看免费成人av毛片| 狂野欧美激情性xxxx在线观看| 另类精品久久| 国产精品一区二区在线观看99| 日韩中文字幕视频在线看片| 国产精品一二三区在线看| 久久鲁丝午夜福利片| 国产成人免费观看mmmm| 最近的中文字幕免费完整| 国产成人一区二区在线| 极品人妻少妇av视频| 亚洲内射少妇av| a级毛片黄视频| 国产毛片在线视频| 草草在线视频免费看| 毛片一级片免费看久久久久| 夜夜爽夜夜爽视频| 久久久午夜欧美精品| 亚洲欧美日韩卡通动漫| 国产高清三级在线| 在线观看免费日韩欧美大片 | av电影中文网址| 一区二区日韩欧美中文字幕 | 三上悠亚av全集在线观看| 人体艺术视频欧美日本| 多毛熟女@视频| 香蕉精品网在线| 一级毛片黄色毛片免费观看视频| 狂野欧美激情性xxxx在线观看| 狠狠精品人妻久久久久久综合| 亚洲成人av在线免费| 色吧在线观看| 国产无遮挡羞羞视频在线观看| videosex国产| 久久久国产一区二区| 超碰97精品在线观看| av视频免费观看在线观看| 熟女电影av网| 午夜免费鲁丝| 国产色爽女视频免费观看| a级毛片在线看网站| 国产av国产精品国产| 国产日韩欧美视频二区| 色94色欧美一区二区| 少妇的逼好多水| 久久久久久久大尺度免费视频| 精品国产露脸久久av麻豆| 国产亚洲最大av| 国产无遮挡羞羞视频在线观看| 国产av码专区亚洲av| 少妇的逼好多水| 丰满饥渴人妻一区二区三| 免费看不卡的av| 日韩精品有码人妻一区| 夜夜骑夜夜射夜夜干| 精品一品国产午夜福利视频| 日本wwww免费看| 久久久久网色| 日本免费在线观看一区| 大香蕉97超碰在线| 亚洲国产精品一区二区三区在线| 中文乱码字字幕精品一区二区三区| 99精国产麻豆久久婷婷| 午夜福利影视在线免费观看| 狠狠婷婷综合久久久久久88av| 亚洲国产最新在线播放| 女性生殖器流出的白浆| 亚洲成人手机| .国产精品久久| 午夜福利视频精品| 免费高清在线观看日韩| av国产久精品久网站免费入址| 亚洲人成网站在线观看播放| 亚洲欧洲国产日韩| 亚洲av不卡在线观看| 22中文网久久字幕| 好男人视频免费观看在线| 午夜免费观看性视频| 国产精品无大码| 性色avwww在线观看| 国产精品不卡视频一区二区| 99热网站在线观看| videossex国产| 国产午夜精品一二区理论片| 国产伦精品一区二区三区视频9| 欧美日韩av久久| 亚洲精品美女久久av网站| 18在线观看网站| av国产久精品久网站免费入址| 两个人免费观看高清视频| 婷婷色综合大香蕉| 亚洲av免费高清在线观看| 亚洲图色成人| 在线看a的网站| 人人澡人人妻人| 十八禁高潮呻吟视频| .国产精品久久| 欧美少妇被猛烈插入视频| 中文字幕亚洲精品专区| 久久人妻熟女aⅴ| 男女国产视频网站| 妹子高潮喷水视频| 我的老师免费观看完整版| 黄色毛片三级朝国网站| 人妻 亚洲 视频| 人人妻人人澡人人看| 精品一区二区三区视频在线| 久久精品人人爽人人爽视色| 国产精品三级大全| 丝瓜视频免费看黄片| 看免费成人av毛片| 久久午夜综合久久蜜桃| 青春草视频在线免费观看| 国产一区二区在线观看日韩| 久久精品夜色国产| 国产av精品麻豆| 色婷婷av一区二区三区视频| 午夜激情福利司机影院| 99精国产麻豆久久婷婷| av卡一久久| 成年av动漫网址| 18在线观看网站| 国产在视频线精品| 在线精品无人区一区二区三| 亚洲色图 男人天堂 中文字幕 | 日日撸夜夜添| 久久人人爽人人片av| 免费av不卡在线播放| 免费大片18禁| 日韩视频在线欧美| 日韩不卡一区二区三区视频在线| a级毛色黄片| 亚洲第一区二区三区不卡| 高清不卡的av网站| 狂野欧美激情性bbbbbb| 成年美女黄网站色视频大全免费 | 国产淫语在线视频| 男人操女人黄网站| 麻豆乱淫一区二区| 久久热精品热| 久久国产亚洲av麻豆专区| 国产成人一区二区在线| 欧美国产精品一级二级三级| 嘟嘟电影网在线观看| 国产免费现黄频在线看| xxxhd国产人妻xxx| 午夜免费男女啪啪视频观看| 只有这里有精品99| 丰满饥渴人妻一区二区三| 精品人妻偷拍中文字幕| 一本大道久久a久久精品| 国产成人精品福利久久| 久久人人爽av亚洲精品天堂| 91成人精品电影| 亚洲色图 男人天堂 中文字幕 | 欧美日韩国产mv在线观看视频| 午夜福利,免费看| 自拍欧美九色日韩亚洲蝌蚪91| 自线自在国产av| 搡女人真爽免费视频火全软件| 啦啦啦啦在线视频资源| 22中文网久久字幕| 美女cb高潮喷水在线观看| 亚洲,欧美,日韩| 久久人人爽人人爽人人片va| 黑人巨大精品欧美一区二区蜜桃 | videosex国产| 9色porny在线观看| 亚洲国产精品国产精品| av播播在线观看一区| 亚洲欧美日韩另类电影网站| 男女边摸边吃奶| 简卡轻食公司| 男的添女的下面高潮视频| 超碰97精品在线观看| 精品人妻熟女毛片av久久网站| 久久午夜综合久久蜜桃| 少妇被粗大的猛进出69影院 | 久久久久精品性色| 夜夜看夜夜爽夜夜摸| 卡戴珊不雅视频在线播放| 人妻 亚洲 视频| 九色成人免费人妻av| 国产色爽女视频免费观看| 亚洲美女视频黄频| 热re99久久精品国产66热6| 肉色欧美久久久久久久蜜桃| 91久久精品国产一区二区成人| 欧美激情国产日韩精品一区| 中文字幕最新亚洲高清| 一个人免费看片子| 久久精品国产自在天天线| 一边摸一边做爽爽视频免费| 日本与韩国留学比较| 国产精品99久久99久久久不卡 | av在线老鸭窝| 欧美最新免费一区二区三区| 国产免费视频播放在线视频| 中国美白少妇内射xxxbb| 各种免费的搞黄视频| 久久久精品区二区三区| 欧美亚洲 丝袜 人妻 在线| 久久人人爽人人片av| 22中文网久久字幕| 简卡轻食公司| 欧美精品人与动牲交sv欧美| av卡一久久| 91精品三级在线观看| 免费大片18禁| 国产毛片在线视频| 99久久中文字幕三级久久日本| 大又大粗又爽又黄少妇毛片口| 久久ye,这里只有精品| 伊人久久精品亚洲午夜| 午夜福利视频精品| 黄色欧美视频在线观看| 最近最新中文字幕免费大全7| 在线 av 中文字幕| 人人妻人人爽人人添夜夜欢视频| 国产精品人妻久久久影院| 少妇 在线观看| 国产精品 国内视频| 老司机影院成人| 成年美女黄网站色视频大全免费 | 少妇丰满av| 久久精品国产亚洲av天美| 九色成人免费人妻av| 久久精品国产亚洲av涩爱| 亚洲成人手机| 天天躁夜夜躁狠狠久久av| 国产精品蜜桃在线观看| 欧美精品国产亚洲| 久久精品国产亚洲av涩爱| 婷婷色综合www| 亚洲激情五月婷婷啪啪| 亚洲精品乱码久久久v下载方式| 五月玫瑰六月丁香| 亚洲在久久综合| 人妻人人澡人人爽人人| 欧美日韩成人在线一区二区| 久久精品国产亚洲av天美| 国产精品三级大全| 伦理电影大哥的女人| 美女国产高潮福利片在线看| 午夜av观看不卡| videos熟女内射| 搡老乐熟女国产| 免费av不卡在线播放| 欧美xxxx性猛交bbbb| 亚洲欧洲国产日韩| 80岁老熟妇乱子伦牲交| 丝袜美足系列| 久久久午夜欧美精品| 亚洲精品一区蜜桃| 日韩强制内射视频| 99精国产麻豆久久婷婷| 三级国产精品片| 久久久久久伊人网av| 国产不卡av网站在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲av国产av综合av卡| 极品少妇高潮喷水抽搐| 成人18禁高潮啪啪吃奶动态图 | 国产av国产精品国产| 满18在线观看网站| 嘟嘟电影网在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 一个人看视频在线观看www免费| 2021少妇久久久久久久久久久| 欧美亚洲 丝袜 人妻 在线| 国产日韩欧美亚洲二区| 黑人欧美特级aaaaaa片| 国产熟女欧美一区二区| 国产视频内射| 女人精品久久久久毛片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 如日韩欧美国产精品一区二区三区 | 中国美白少妇内射xxxbb| 午夜福利影视在线免费观看| 中文欧美无线码| 免费少妇av软件| 久久久久久久精品精品| 大话2 男鬼变身卡| 中文字幕制服av| 久久人人爽人人爽人人片va| 国产精品.久久久| 人妻少妇偷人精品九色| 国产欧美日韩综合在线一区二区| 久久精品熟女亚洲av麻豆精品| 日日摸夜夜添夜夜爱| 亚洲精品aⅴ在线观看| 久久久久久人妻| 亚洲精品国产色婷婷电影| 久久女婷五月综合色啪小说| 97精品久久久久久久久久精品| 亚洲美女视频黄频| 欧美老熟妇乱子伦牲交| 亚洲综合精品二区| 成人黄色视频免费在线看| 久久狼人影院| 嫩草影院入口| 国产精品偷伦视频观看了| 亚洲国产毛片av蜜桃av| 国产av国产精品国产| 一级,二级,三级黄色视频| 国产日韩欧美亚洲二区| 色视频在线一区二区三区| 国产女主播在线喷水免费视频网站| 女的被弄到高潮叫床怎么办| 欧美少妇被猛烈插入视频| 人体艺术视频欧美日本| 蜜桃国产av成人99| 日本猛色少妇xxxxx猛交久久| 国产精品一区二区三区四区免费观看| 51国产日韩欧美| 美女国产视频在线观看| 亚洲久久久国产精品| 欧美日韩国产mv在线观看视频| 久久久午夜欧美精品| 国产一区有黄有色的免费视频| 久久97久久精品| 亚洲激情五月婷婷啪啪| 一个人免费看片子| 高清不卡的av网站| 亚洲精品久久久久久婷婷小说| 在线观看一区二区三区激情| av在线播放精品| 亚洲一区二区三区欧美精品| 少妇精品久久久久久久| 亚洲精品久久午夜乱码| 我的老师免费观看完整版| 乱人伦中国视频| 午夜免费鲁丝| 男女啪啪激烈高潮av片| 少妇精品久久久久久久| 卡戴珊不雅视频在线播放| 只有这里有精品99| 免费人成在线观看视频色| 欧美日韩国产mv在线观看视频| 一区在线观看完整版| 久久毛片免费看一区二区三区| 美女cb高潮喷水在线观看| 最新中文字幕久久久久| 久久av网站|