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

    前緣和轉(zhuǎn)軸影響翼?yè)u滾特性的數(shù)值模擬*

    2021-06-07 08:33:16楊小亮黃世林柴振霞孫喜萬(wàn)
    關(guān)鍵詞:三角翼翼面轉(zhuǎn)軸

    楊小亮,黃世林,楊 起,柴振霞,孫喜萬(wàn),劉 偉

    (1. 國(guó)防科技大學(xué) 空天科學(xué)學(xué)院, 湖南 長(zhǎng)沙 410073; 2. 中國(guó)人民解放軍93108部隊(duì), 黑龍江 齊齊哈爾 161007)

    三角翼因結(jié)構(gòu)強(qiáng)度高、可提供非線性渦升力等優(yōu)點(diǎn),是現(xiàn)代戰(zhàn)斗機(jī)主要采用的機(jī)翼構(gòu)型,其氣動(dòng)特性受到了學(xué)界的普遍關(guān)注,尤其是其動(dòng)態(tài)特性密切關(guān)系到現(xiàn)代戰(zhàn)斗機(jī)大攻角機(jī)動(dòng)飛行的安全[1]。大攻角條件下,三角翼會(huì)產(chǎn)生大幅度自維持的翼?yè)u滾現(xiàn)象。雖然早在20世紀(jì)中葉,人們就觀察到了翼?yè)u滾現(xiàn)象,但由于其復(fù)雜的非線性氣動(dòng)特性,直到20世紀(jì)80年代才陸續(xù)出現(xiàn)對(duì)三角翼?yè)u滾現(xiàn)象較為系統(tǒng)的試驗(yàn)研究[2-5],而耦合求解N-S方程和Euler剛體動(dòng)力學(xué)方程[6-10]完成三角翼?yè)u滾現(xiàn)象的數(shù)值模擬則是二十多年來(lái)才開(kāi)始的事情。

    準(zhǔn)確捕捉位于三角翼背風(fēng)面的非線性集中渦是數(shù)值模擬三角翼繞流時(shí)必須解決的問(wèn)題。近年來(lái),大量文獻(xiàn)在模擬、分析及控制三角翼上渦流特性方面不遺余力[11-15],進(jìn)展卓有成效[16],而數(shù)值模擬渦流主導(dǎo)的三角翼?yè)u滾現(xiàn)象則主要集中在捕捉搖滾現(xiàn)象的定性研究。三角翼上集中渦形成和發(fā)展具有強(qiáng)非線性,故對(duì)三角翼構(gòu)型及影響因素極為敏感[17],造成了數(shù)值模擬與實(shí)驗(yàn)之間、實(shí)驗(yàn)與實(shí)驗(yàn)之間往往存在明顯的分歧[5]。提升數(shù)值模擬和風(fēng)洞實(shí)驗(yàn)反映翼?yè)u滾現(xiàn)象物理本質(zhì)的水平,減小數(shù)值模擬和實(shí)驗(yàn)差距,必須研究翼?yè)u滾特性的影響因素。前緣構(gòu)型和滾轉(zhuǎn)軸安裝位置是實(shí)驗(yàn)研究必須陳述的關(guān)鍵細(xì)節(jié),且文獻(xiàn)[4-5]表明不同的處理可導(dǎo)致?lián)u滾特性顯著不同。閻超等[18]研究了前緣形狀對(duì)三角翼上渦流特性的影響,但對(duì)搖滾特性影響的研究則相對(duì)較少。本文耦合求解N-S方程和剛體動(dòng)力學(xué)方程,針對(duì)80°后掠三角翼研究前緣構(gòu)型和轉(zhuǎn)軸安裝位置影響三角翼?yè)u滾特性的程度和規(guī)律。

    1 模型及計(jì)算方法

    1.1 三角翼模型及網(wǎng)格

    研究采用80°后掠尖前緣平板三角翼模型,前緣分別采用下削尖、雙面削尖、上削尖三種方式,削尖角均為45°,尾緣不做削尖處理,垂直翼面,如圖1所示。為研究滾轉(zhuǎn)軸對(duì)細(xì)長(zhǎng)三角翼?yè)u滾特性的影響,基于下削尖前緣三角翼分別研究了滾轉(zhuǎn)軸安裝在四種不同位置時(shí)的搖滾特性,如圖2所示,滾轉(zhuǎn)軸位置分別安裝在下翼面下方5 mm處(記為Axel_DD)、下翼面上(記為Axel_D)、上翼面上(記為Axel_U)、上翼面上方5 mm處(記為Axel_UU),四種安裝位置的滾轉(zhuǎn)軸均位于三角翼體對(duì)稱面上。網(wǎng)格生成采用OH型拓?fù)浣Y(jié)構(gòu),周向O型、其余方向H型,以平衡描述精度和網(wǎng)格規(guī)模[19]。計(jì)算網(wǎng)格如圖3所示,網(wǎng)格上游邊界距離三角翼頂點(diǎn)2.5倍根弦長(zhǎng),周向邊界距離體軸2.5倍根弦長(zhǎng),下游邊界距離三角翼尾緣5倍根弦長(zhǎng)。網(wǎng)格數(shù)為112×226×50(流向×周向×法向),翼尖、尾緣及翼前緣附近流動(dòng)較為復(fù)雜,局部適當(dāng)加密。

    (a) 下削尖前緣(a) Down beveled leading edge

    (b) 雙面削尖前緣(b) Double beveled leading edge

    (c) 上削尖前緣(c) Up beveled leading edge圖1 三種不同前緣形狀的80°后掠三角翼模型Fig.1 Three 80° swept delta wing model with different leading edge configuration

    (a) 滾轉(zhuǎn)軸位于下翼面下方(Axel_DD)(a) Roll-axis under the lower wing surface(Axel_DD)

    (b) 滾轉(zhuǎn)軸位于下翼面(Axel_D)(b) Roll-axis on the lower wing surface(Axel_D)

    (c) 滾轉(zhuǎn)軸位于上翼面(Axel_U)(c) Roll-axis on the upper wing surface(Axel_U)

    (d)滾轉(zhuǎn)軸位于上翼面上方(Axel_UU)(d) Roll-axis above the upper wing surface(Axel_UU)圖2 四種不同滾轉(zhuǎn)軸位置示意圖Fig.2 Schema of four roll-axis installation of delta wing

    圖3 80°后掠三角翼網(wǎng)格分布Fig.3 Grid distribution of the 80° swept delta wing

    1.2 數(shù)值方法

    采用非定常N-S方程組和滾轉(zhuǎn)動(dòng)力學(xué)方程描述三角翼單自由度搖滾運(yùn)動(dòng)[20],貼體坐標(biāo)系(t,ξ,η,ζ)下的方程寫(xiě)為:

    (1)

    基于結(jié)構(gòu)網(wǎng)格的有限體積法,采用Roe格式和含雙時(shí)間步的LU-SGS方法離散非定常N-S方程組。滾轉(zhuǎn)動(dòng)力學(xué)方程采用時(shí)間二階精度的單邊差分離散,得到時(shí)間二階精度的滾轉(zhuǎn)差分方程:

    (2)

    非定常N-S方程組和滾轉(zhuǎn)動(dòng)力學(xué)方程的耦合求解采用松耦合模式,即由滾轉(zhuǎn)動(dòng)力學(xué)方程提供三角翼n+1時(shí)刻的位置,再由非定常N-S方程組計(jì)算n+1時(shí)刻的滾轉(zhuǎn)力矩系數(shù)。

    遠(yuǎn)場(chǎng)采用無(wú)反射邊界條件,壁面邊界取絕熱壁,壓力條件為?p/?n|wall=-ρa(bǔ)wall·nwall,速度條件為無(wú)滑移條件,即V=Vwall。

    1.3 驗(yàn)證與確認(rèn)

    本文的計(jì)算基于國(guó)防科技大學(xué)劉偉教授團(tuán)隊(duì)開(kāi)發(fā)的飛行器動(dòng)態(tài)特性研究程序(Aircraft Dynamic Characteristics Research Program, ADCRP)。該程序針對(duì)多塊對(duì)接結(jié)構(gòu)網(wǎng)格,求解三維非定常RANS方程,經(jīng)過(guò)二十余年的發(fā)展,具備復(fù)雜飛行器流動(dòng)/運(yùn)動(dòng)耦合的多自由度非定過(guò)程的數(shù)值模擬能力。該程序采用系列標(biāo)模外形進(jìn)行了系統(tǒng)的驗(yàn)證與確認(rèn)計(jì)算,適用于三角翼等定、動(dòng)姿態(tài)非定常問(wèn)題的分析[20-21],在多個(gè)型號(hào)項(xiàng)目中得到了應(yīng)用。該程序具備模擬細(xì)長(zhǎng)三角翼?yè)u滾現(xiàn)象的研究能力。

    2 前緣構(gòu)型對(duì)搖滾特性的影響

    考察前緣構(gòu)型對(duì)細(xì)長(zhǎng)三角翼?yè)u滾特性的影響。計(jì)算時(shí),給定來(lái)流馬赫數(shù)Ma=0.2,基于根弦長(zhǎng)度的雷諾數(shù)Re=0.4×106,無(wú)量綱滾轉(zhuǎn)轉(zhuǎn)動(dòng)慣量Ixx=0.1,滾轉(zhuǎn)軸與x軸重合且位于上翼面,不計(jì)軸承阻尼。統(tǒng)一給定10°的初始滾轉(zhuǎn)角,研究不同前緣三角翼的滾轉(zhuǎn)響應(yīng),分別對(duì)分岔攻角和搖滾振蕩的振幅進(jìn)行討論。

    2.1 前緣形狀對(duì)分岔攻角的影響

    文獻(xiàn)[22]指出,三角翼攻角α變化會(huì)導(dǎo)致其橫向穩(wěn)定性發(fā)生Hopf分岔,存在攻角的臨界值,攻角小于該值,三角翼橫向動(dòng)態(tài)穩(wěn)定,大于該值,三角翼發(fā)生橫向失穩(wěn),形成大幅度自維持的極限環(huán)振蕩,該臨界值稱為分岔攻角。

    對(duì)分岔攻角的數(shù)值模擬結(jié)果如圖4所示,顯示了下削尖、雙面削尖、上削尖三種前緣構(gòu)型的三角翼不同攻角下的滾轉(zhuǎn)角時(shí)間歷程。按照文獻(xiàn)描述,易于判斷下削尖前緣三角翼分岔攻角約為22.4°,雙面削尖前緣三角翼分岔攻角約為16.9°,上削尖前緣三角翼分岔攻角約為15.5°。其中:下削尖前緣三角翼能夠在較大的攻角范圍保持滾轉(zhuǎn)穩(wěn)定性,分岔攻角最大;上削尖前緣三角翼分岔攻角最小,即最容易發(fā)生橫向失穩(wěn);雙面削尖前緣的三角翼橫向穩(wěn)定性介于二者之間。

    (a) 下削尖前緣(a) Down beveled leading edge

    (b) 雙面削尖前緣(b) Double beveled leading edge

    (c) 上削尖前緣(c) Up beveled leading edge圖4 不同前緣形狀三角翼的滾轉(zhuǎn)角時(shí)間歷程Fig.4 Time history of roll angle of delta wing with different leading edge

    為進(jìn)一步理解前緣構(gòu)型對(duì)三角翼動(dòng)態(tài)特性的影響,基于強(qiáng)迫簡(jiǎn)諧分析法計(jì)算三角翼滾轉(zhuǎn)穩(wěn)定性參數(shù),如表1和表2所示。穩(wěn)定性參數(shù)為負(fù)則表示具有穩(wěn)定性,為正則表示不穩(wěn)定。由表1可知:不同前緣構(gòu)型的三角翼,在所研究的攻角范圍內(nèi)都是靜穩(wěn)定的,整體而言,靜穩(wěn)定性相當(dāng)。由表2可知:隨著攻角增大,滾轉(zhuǎn)阻尼導(dǎo)數(shù)由負(fù)變正,表明從滾轉(zhuǎn)動(dòng)穩(wěn)定變?yōu)闈L轉(zhuǎn)動(dòng)不穩(wěn)定,從動(dòng)態(tài)穩(wěn)定性的角度說(shuō)明了三角翼滾轉(zhuǎn)方向發(fā)生Hopf分岔的這一現(xiàn)象。其中,下削尖三角翼在22.4°、雙面削尖三角翼在16.9°、上削尖三角翼在15.5°時(shí),滾轉(zhuǎn)阻尼導(dǎo)數(shù)的絕對(duì)值顯著小于其他角度的絕對(duì)值,這表明此時(shí)三角翼幾乎處于無(wú)阻尼狀態(tài),正好對(duì)應(yīng)三角翼在跨過(guò)該攻角時(shí)的滾轉(zhuǎn)穩(wěn)定性性態(tài)的變化,分岔攻角在以上角度附近。

    表1 三角翼滾轉(zhuǎn)靜穩(wěn)定性導(dǎo)數(shù)

    表2 三角翼滾轉(zhuǎn)阻尼導(dǎo)數(shù)

    2.2 前緣形狀對(duì)振幅特性的影響

    上一節(jié)的研究表明,攻角大于分岔攻角以后,均會(huì)發(fā)生橫向失穩(wěn),形成大幅度自維持極限環(huán)形式的搖滾現(xiàn)象。為進(jìn)一步考察不同前緣構(gòu)型的三角翼橫向失穩(wěn)以后的動(dòng)態(tài)特性,選取25°和30°攻角,給定0°初始滾轉(zhuǎn)角,研究前緣形狀對(duì)三角翼振幅特性的影響,結(jié)果如圖5所示。

    25°攻角時(shí),如圖5(a)所示,三種前緣構(gòu)型的三角翼均形成了極限環(huán)形式的等幅振蕩,下削尖前緣三角翼的振幅約為19.2°,雙面削尖前緣三角翼的振幅約為65.2°,上削尖前緣三角翼的振幅約為75.6°。其中,下削尖三角翼振幅最小,雙面削尖三角翼次之,上削尖三角翼的振幅最大,對(duì)比表2數(shù)據(jù),三種前緣構(gòu)型的三角翼的滾轉(zhuǎn)阻尼導(dǎo)數(shù)均為正,從數(shù)值上而言,下削尖三角翼最小,雙面削尖三角翼次之,上削尖三角翼最大,三角翼自激滾轉(zhuǎn)振蕩的表現(xiàn)與動(dòng)態(tài)穩(wěn)定性參數(shù)預(yù)示的規(guī)律一致。

    30°攻角時(shí),如圖5(b)所示,三種機(jī)翼均最終形成了極限環(huán)形式的等幅振蕩,但雙面削尖和上削尖前緣的三角翼在發(fā)生翻轉(zhuǎn)后繞滾轉(zhuǎn)軸進(jìn)行滾轉(zhuǎn)振蕩。其中,雙面削尖前緣三角翼?yè)u滾振幅約為49.8°,上削尖前緣三角翼?yè)u滾振幅約為31.4°,下削尖前緣三角翼未發(fā)生翻轉(zhuǎn),其振幅約為37.1°。

    理論上講,雙面削尖前緣三角翼,上下對(duì)稱,翻轉(zhuǎn)后與自身同構(gòu),上削尖前緣三角翼與下削尖前緣三角翼在結(jié)構(gòu)上相同,僅與來(lái)流在空間的相對(duì)位置不同。上削尖前緣三角翼翻轉(zhuǎn)后,空間位置與下削尖前緣三角翼相同。但由圖5(b)觀察到的結(jié)果,上削尖前緣三角翼翻轉(zhuǎn)后的振幅(31.4°)與下削尖前緣三角翼的振幅(37.1°)并不相同。經(jīng)分析,主要由于翻轉(zhuǎn)后滾轉(zhuǎn)軸位置不相同,對(duì)于雙面削尖的三角翼同樣存在該問(wèn)題,需要進(jìn)一步分析。

    (a) α=25°

    (b) α=30° 圖5 不同前緣構(gòu)型三角翼在大攻角下的搖滾特性Fig.5 Wing rock characteristics of delta wing with different leading edge in high angle of attack

    3 滾轉(zhuǎn)軸位置對(duì)搖滾特性的影響

    基于對(duì)三角翼?yè)u滾現(xiàn)象的理解,進(jìn)一步考察滾轉(zhuǎn)軸位置對(duì)細(xì)長(zhǎng)三角翼?yè)u滾特性的影響。計(jì)算時(shí),給定來(lái)流馬赫數(shù)Ma=0.2,基于根弦長(zhǎng)度的雷諾數(shù)Re=0.4×106,無(wú)量綱滾轉(zhuǎn)轉(zhuǎn)動(dòng)慣量Ixx=0.1,統(tǒng)一指定10°的初始滾轉(zhuǎn)角,不計(jì)軸承阻尼。選取下削尖三角翼,指定不同滾轉(zhuǎn)軸位置,如圖2所示,研究滾轉(zhuǎn)軸位置對(duì)搖滾特性的影響,分別對(duì)分岔攻角和搖滾振幅進(jìn)行了討論。

    3.1 滾轉(zhuǎn)軸位置對(duì)分岔攻角的影響

    如圖6(a)所示,滾轉(zhuǎn)軸位于下翼面下方5 mm處時(shí),三角翼從10°初始滾轉(zhuǎn)角釋放后,在25°攻角下?lián)u滾振幅呈衰減狀態(tài),在30°攻角下?lián)u滾振幅呈增大狀態(tài),而在27°攻角時(shí)搖滾振幅基本保持10°不變,說(shuō)明滾轉(zhuǎn)軸位于下翼面下方5 mm的三角翼的分岔攻角約為27°。由圖6(b)~(d)易知,滾轉(zhuǎn)軸位于下翼面的三角翼分岔攻角為26.5°,滾轉(zhuǎn)軸位于上翼面的三角翼分岔攻角為22.4°,滾轉(zhuǎn)軸位于上翼面上方5 mm處的三角翼的分岔攻角為17.4°。由此可見(jiàn),隨著滾轉(zhuǎn)軸從下翼面下方5 mm處移動(dòng)到上翼面上方5 mm處,三角翼的分岔攻角逐漸減小??梢缘贸鼋Y(jié)論,在所研究的滾轉(zhuǎn)軸位置范圍內(nèi),滾轉(zhuǎn)軸位置越高,三角翼?yè)u滾的分岔攻角越小,越易于橫向失穩(wěn),反之亦然。

    (a) 滾轉(zhuǎn)軸位于下翼面以下(Axel_DD)(a) Roll-axis under the lower wing surface(Axel_DD)

    (b) 滾轉(zhuǎn)軸位于下翼面(Axel_D)(b) Roll-axis on the lower wing surface(Axel_D)

    (c) 滾轉(zhuǎn)軸位于上翼面(Axel_U)(c) Roll-axis on the upper wing surface(Axel_U)

    (d) 滾轉(zhuǎn)軸位于上翼面以上(Axel_UU)(d) Roll-axis above the upper wing surface(Axel_UU)圖6 不同滾轉(zhuǎn)軸位置的三角翼滾轉(zhuǎn)角時(shí)間歷程曲線Fig.6 Roll angle time historical solution of delta wing installed at different roll-axes

    3.2 滾轉(zhuǎn)軸位置對(duì)振幅特性的影響

    進(jìn)一步研究滾轉(zhuǎn)軸位置對(duì)三角翼振幅特性的影響,給定30°、40°、50°和60°攻角,模擬四種滾轉(zhuǎn)軸位置的三角翼?yè)u滾振幅特性,計(jì)算結(jié)果如圖7所示。

    從圖7中可以看出,四種滾轉(zhuǎn)軸位置下?lián)u滾振幅隨攻角變化的規(guī)律基本一致,30°攻角時(shí),四種滾轉(zhuǎn)軸位置的三角翼均形成了極限環(huán)自激搖滾,滾轉(zhuǎn)軸位置越靠上振幅越大;40°攻角左右建立最大振幅,隨著攻角進(jìn)一步增大,由于受到渦破裂的影響,搖滾振幅隨攻角增大而減??;當(dāng)攻角增加到60°以后,四種滾轉(zhuǎn)軸位置的三角翼均難以形成大幅度的搖滾振幅。

    圖7 搖滾振幅隨攻角變化曲線Fig.7 Amplitude of roll varies with attack angle

    4 結(jié)論

    耦合求解三維非定常N-S方程和歐拉剛體動(dòng)力學(xué)方程組,研究了前緣構(gòu)型與滾轉(zhuǎn)軸位置對(duì)80°后掠平板三角翼?yè)u滾特性的影響,得到以下結(jié)論:

    1)前緣削尖方式不改變?nèi)且頇M向動(dòng)態(tài)穩(wěn)定性的性態(tài),但影響三角翼?yè)u滾的分岔攻角和搖滾振幅;不同前緣構(gòu)型三角翼的自激搖滾特性與動(dòng)態(tài)穩(wěn)定性參數(shù)表征的穩(wěn)定性一致。

    2)對(duì)于轉(zhuǎn)軸位于上表面的三角翼,攻角25°時(shí),前緣下削尖時(shí)搖滾振幅最小,上削尖時(shí)搖滾振幅最大;攻角30°時(shí),雙面削尖前緣和上削尖前緣的三角翼自激搖滾時(shí)發(fā)生了翻轉(zhuǎn),翻轉(zhuǎn)后繞滾轉(zhuǎn)軸等幅振蕩,雙面削尖前緣與上削尖前緣三角翼先翻轉(zhuǎn)后再形成極限環(huán)形式的等幅振蕩。

    3)滾轉(zhuǎn)軸的安裝位置影響三角翼的搖滾特性,滾轉(zhuǎn)軸位置越高,三角翼?yè)u滾的分岔攻角越小,極限環(huán)滾轉(zhuǎn)振蕩的振幅越大,反之亦然。

    三角翼?yè)u滾現(xiàn)象研究多年依然是熱點(diǎn)問(wèn)題,對(duì)其產(chǎn)生和維持機(jī)理已有較多深刻的認(rèn)識(shí),但對(duì)敏感因素的分析還不足,旋渦等流動(dòng)結(jié)構(gòu)與飛行器之間作用的過(guò)程和機(jī)理等問(wèn)題仍值得進(jìn)一步的分析。

    猜你喜歡
    三角翼翼面轉(zhuǎn)軸
    基于拔銷器鎖定的飛行器氣動(dòng)控制面解鎖控制方法
    三角翼機(jī)翼?yè)u滾主動(dòng)控制多學(xué)科耦合數(shù)值模擬
    大型汽輪發(fā)電機(jī)轉(zhuǎn)軸接地方式及軸電流分析
    大擴(kuò)散角明渠中導(dǎo)流三角翼體形優(yōu)化試驗(yàn)
    人民黃河(2020年11期)2020-12-14 04:32:27
    軋機(jī)轉(zhuǎn)軸無(wú)損檢測(cè)及修復(fù)技術(shù)
    山東冶金(2019年2期)2019-05-11 09:12:24
    固定翼二維彈道修正引信升力翼面位置的影響
    基于Aerobook平臺(tái)的復(fù)合材料翼面結(jié)構(gòu)設(shè)計(jì)流程
    CY—06三角翼無(wú)人機(jī)
    航空模型(2016年10期)2017-05-09 06:22:13
    小細(xì)節(jié)大功效 淺談筆記本屏幕轉(zhuǎn)軸設(shè)計(jì)
    不同后掠角三角翼的靜態(tài)地面效應(yīng)數(shù)值模擬
    精品国产露脸久久av麻豆| 成人影院久久| 热99国产精品久久久久久7| 国产91av在线免费观看| 国产一区二区三区综合在线观看 | 性色av一级| 国产色爽女视频免费观看| 女性生殖器流出的白浆| 午夜福利在线在线| 国产精品三级大全| 国产免费又黄又爽又色| 夫妻性生交免费视频一级片| 成人国产麻豆网| 日韩成人av中文字幕在线观看| 免费在线观看成人毛片| 国产av一区二区精品久久 | 国产av国产精品国产| 人妻系列 视频| 在线观看免费视频网站a站| 在线天堂最新版资源| 3wmmmm亚洲av在线观看| 亚洲色图综合在线观看| 热re99久久精品国产66热6| 国产在线男女| 亚洲精品自拍成人| 中国国产av一级| 夜夜骑夜夜射夜夜干| 午夜精品国产一区二区电影| 国产精品一区二区三区四区免费观看| 日韩不卡一区二区三区视频在线| 国产高清有码在线观看视频| 欧美高清性xxxxhd video| 久久精品国产亚洲av涩爱| 91久久精品电影网| 亚洲国产欧美在线一区| 嫩草影院入口| 国产日韩欧美亚洲二区| 日本-黄色视频高清免费观看| 99久久精品热视频| 97热精品久久久久久| 狠狠精品人妻久久久久久综合| 国产av一区二区精品久久 | 欧美xxxx黑人xx丫x性爽| 久久久久久久大尺度免费视频| 国产毛片在线视频| 午夜福利在线在线| 99九九线精品视频在线观看视频| 嘟嘟电影网在线观看| 色5月婷婷丁香| 亚洲av中文av极速乱| 男女边吃奶边做爰视频| av女优亚洲男人天堂| 国产成人精品一,二区| 亚洲国产高清在线一区二区三| 国产成人免费观看mmmm| 国产精品嫩草影院av在线观看| 国产在线一区二区三区精| 伦理电影大哥的女人| 观看av在线不卡| 一区二区三区四区激情视频| 国产乱人偷精品视频| 五月开心婷婷网| 亚洲成人中文字幕在线播放| 王馨瑶露胸无遮挡在线观看| 高清毛片免费看| .国产精品久久| 国产伦精品一区二区三区视频9| 91精品国产国语对白视频| 内地一区二区视频在线| 波野结衣二区三区在线| 大码成人一级视频| 狂野欧美激情性bbbbbb| 亚洲国产精品999| 久久99精品国语久久久| 啦啦啦视频在线资源免费观看| 国产精品蜜桃在线观看| 国产男女超爽视频在线观看| 亚州av有码| 亚洲av免费高清在线观看| 亚洲av成人精品一区久久| 九九在线视频观看精品| 亚洲美女视频黄频| 99久国产av精品国产电影| 国产色爽女视频免费观看| 亚洲在久久综合| 亚洲精品aⅴ在线观看| 在线观看一区二区三区| 噜噜噜噜噜久久久久久91| 欧美xxⅹ黑人| 在线观看一区二区三区| 青春草亚洲视频在线观看| 亚洲成人手机| 99久久精品国产国产毛片| 美女高潮的动态| 亚洲性久久影院| 哪个播放器可以免费观看大片| av国产精品久久久久影院| 中文字幕免费在线视频6| 看免费成人av毛片| 亚洲av.av天堂| 一本一本综合久久| 蜜桃在线观看..| 蜜臀久久99精品久久宅男| 亚洲国产精品一区三区| 国产无遮挡羞羞视频在线观看| 在线精品无人区一区二区三 | 国产久久久一区二区三区| 女性生殖器流出的白浆| 亚洲av福利一区| 日本欧美国产在线视频| 在线播放无遮挡| 亚洲美女黄色视频免费看| 久久97久久精品| 亚洲精品乱码久久久v下载方式| .国产精品久久| 黑人高潮一二区| 亚洲欧美中文字幕日韩二区| 不卡视频在线观看欧美| 亚州av有码| 国产精品99久久99久久久不卡 | av在线老鸭窝| 国产亚洲5aaaaa淫片| 黄色怎么调成土黄色| 亚洲美女搞黄在线观看| av.在线天堂| 精品久久久久久久久av| 国产精品蜜桃在线观看| 亚洲av日韩在线播放| 多毛熟女@视频| 啦啦啦啦在线视频资源| 国产av精品麻豆| 大片电影免费在线观看免费| 国产亚洲一区二区精品| 国产熟女欧美一区二区| 国产精品欧美亚洲77777| 99热全是精品| 久久国产精品大桥未久av | 日本欧美国产在线视频| 久久女婷五月综合色啪小说| 十分钟在线观看高清视频www | 三级国产精品片| 不卡视频在线观看欧美| 色吧在线观看| 国产无遮挡羞羞视频在线观看| h日本视频在线播放| 色吧在线观看| 干丝袜人妻中文字幕| 少妇高潮的动态图| 十分钟在线观看高清视频www | 午夜精品国产一区二区电影| 男的添女的下面高潮视频| 在线观看国产h片| 91精品国产九色| av天堂中文字幕网| 亚洲精品国产色婷婷电影| 特大巨黑吊av在线直播| 国产精品一区二区性色av| 国产精品一区二区在线不卡| 中文字幕亚洲精品专区| 欧美精品人与动牲交sv欧美| 少妇猛男粗大的猛烈进出视频| 观看免费一级毛片| 日本av免费视频播放| 欧美日韩国产mv在线观看视频 | 超碰av人人做人人爽久久| 精品99又大又爽又粗少妇毛片| 国产熟女欧美一区二区| 天美传媒精品一区二区| 精品一品国产午夜福利视频| 久久精品久久久久久噜噜老黄| 毛片女人毛片| 国产成人午夜福利电影在线观看| 亚洲怡红院男人天堂| 国产精品嫩草影院av在线观看| 欧美日本视频| 色吧在线观看| 欧美 日韩 精品 国产| 久久久久久久久久人人人人人人| 免费av中文字幕在线| 免费看日本二区| 肉色欧美久久久久久久蜜桃| 亚洲色图综合在线观看| 免费人成在线观看视频色| 国产真实伦视频高清在线观看| 成人无遮挡网站| 国产在线男女| 中国三级夫妇交换| 啦啦啦啦在线视频资源| 欧美精品人与动牲交sv欧美| 成人无遮挡网站| 免费看不卡的av| 久久久久精品久久久久真实原创| 97精品久久久久久久久久精品| 精品酒店卫生间| 高清av免费在线| 一级av片app| 成年美女黄网站色视频大全免费 | 欧美精品国产亚洲| 亚洲人成网站在线播| 在线看a的网站| 视频中文字幕在线观看| 韩国av在线不卡| 国产成人a区在线观看| 欧美三级亚洲精品| 精品午夜福利在线看| 一区二区三区免费毛片| 熟妇人妻不卡中文字幕| 欧美老熟妇乱子伦牲交| 成人国产麻豆网| 国产精品一区www在线观看| 欧美+日韩+精品| av免费在线看不卡| 蜜桃亚洲精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 国产成人免费无遮挡视频| 久久国产精品男人的天堂亚洲 | 蜜桃久久精品国产亚洲av| 边亲边吃奶的免费视频| 欧美xxxx黑人xx丫x性爽| kizo精华| 日本猛色少妇xxxxx猛交久久| 嘟嘟电影网在线观看| av在线蜜桃| 嫩草影院新地址| 久久久久久久亚洲中文字幕| 日韩伦理黄色片| 丝袜脚勾引网站| www.av在线官网国产| 少妇的逼好多水| 亚洲美女搞黄在线观看| 成人黄色视频免费在线看| 久久国内精品自在自线图片| 少妇高潮的动态图| 在线观看国产h片| 国产色爽女视频免费观看| 国产一区二区三区av在线| 久久久色成人| 看免费成人av毛片| 一级毛片久久久久久久久女| 国产午夜精品一二区理论片| 日本vs欧美在线观看视频 | 国产精品不卡视频一区二区| 免费久久久久久久精品成人欧美视频 | 夜夜骑夜夜射夜夜干| 亚洲久久久国产精品| 日韩三级伦理在线观看| 久久久久久久久久久丰满| 国产亚洲5aaaaa淫片| 99久久精品国产国产毛片| 亚洲最大成人中文| 我要看日韩黄色一级片| 男人添女人高潮全过程视频| 美女国产视频在线观看| 三级国产精品片| 亚洲欧美一区二区三区黑人 | 亚洲av不卡在线观看| 国产免费一级a男人的天堂| 国产av国产精品国产| 赤兔流量卡办理| 亚洲av男天堂| 国产一级毛片在线| av播播在线观看一区| 久久久久网色| 亚洲欧美精品自产自拍| 国模一区二区三区四区视频| 网址你懂的国产日韩在线| 激情 狠狠 欧美| 一个人免费看片子| 久久久成人免费电影| 日本黄色片子视频| 日产精品乱码卡一卡2卡三| 亚洲av免费高清在线观看| 性高湖久久久久久久久免费观看| 男女国产视频网站| 秋霞伦理黄片| 色婷婷av一区二区三区视频| 国产深夜福利视频在线观看| 亚洲欧美日韩无卡精品| 秋霞在线观看毛片| 亚洲激情五月婷婷啪啪| 精品人妻熟女av久视频| 亚洲欧美一区二区三区黑人 | 国产在线一区二区三区精| 国产午夜精品久久久久久一区二区三区| 高清黄色对白视频在线免费看 | 九九爱精品视频在线观看| 亚洲真实伦在线观看| av女优亚洲男人天堂| 人体艺术视频欧美日本| 国模一区二区三区四区视频| 久久女婷五月综合色啪小说| 午夜福利视频精品| 女性生殖器流出的白浆| 日日摸夜夜添夜夜爱| 免费看不卡的av| 又粗又硬又长又爽又黄的视频| 国产久久久一区二区三区| 国产真实伦视频高清在线观看| 狂野欧美激情性bbbbbb| 日韩人妻高清精品专区| 国产精品99久久99久久久不卡 | 黄色日韩在线| 高清不卡的av网站| 国产在线视频一区二区| 成人漫画全彩无遮挡| 夜夜爽夜夜爽视频| 啦啦啦在线观看免费高清www| 高清毛片免费看| 成人黄色视频免费在线看| 韩国av在线不卡| 国产高潮美女av| 日本一二三区视频观看| 美女脱内裤让男人舔精品视频| 91精品国产国语对白视频| 在线精品无人区一区二区三 | 看十八女毛片水多多多| 亚洲,一卡二卡三卡| 欧美日本视频| 久久影院123| 成人一区二区视频在线观看| 又粗又硬又长又爽又黄的视频| 老师上课跳d突然被开到最大视频| 五月伊人婷婷丁香| 中文字幕亚洲精品专区| 日韩精品有码人妻一区| 蜜桃久久精品国产亚洲av| 亚洲精品日韩av片在线观看| 如何舔出高潮| 26uuu在线亚洲综合色| 美女视频免费永久观看网站| 日本色播在线视频| 18禁动态无遮挡网站| 91aial.com中文字幕在线观看| 99热6这里只有精品| 22中文网久久字幕| 乱系列少妇在线播放| 我的老师免费观看完整版| 国产免费一区二区三区四区乱码| 免费av不卡在线播放| 国产av国产精品国产| 内射极品少妇av片p| 深爱激情五月婷婷| 国产成人一区二区在线| 精品少妇久久久久久888优播| 日日摸夜夜添夜夜添av毛片| 久久影院123| 国产成人午夜福利电影在线观看| 免费av中文字幕在线| 久久婷婷青草| 久久久精品94久久精品| 91精品一卡2卡3卡4卡| 国国产精品蜜臀av免费| 蜜桃在线观看..| 国产亚洲最大av| 亚洲av中文字字幕乱码综合| 国产高清国产精品国产三级 | 免费观看无遮挡的男女| 国产精品一区二区三区四区免费观看| av卡一久久| 国产白丝娇喘喷水9色精品| 亚洲色图综合在线观看| 成人亚洲精品一区在线观看 | 久久av网站| 国产精品爽爽va在线观看网站| 赤兔流量卡办理| 欧美另类一区| 欧美高清成人免费视频www| 成人特级av手机在线观看| 国产欧美日韩精品一区二区| 夫妻午夜视频| 久久久欧美国产精品| 日本猛色少妇xxxxx猛交久久| 国产淫语在线视频| 欧美另类一区| 国产伦在线观看视频一区| 欧美性感艳星| 国产欧美亚洲国产| 精品久久久久久久末码| 肉色欧美久久久久久久蜜桃| 国产淫片久久久久久久久| 永久免费av网站大全| 在线 av 中文字幕| 22中文网久久字幕| 免费黄色在线免费观看| 久久99精品国语久久久| 晚上一个人看的免费电影| 久久精品国产自在天天线| 狠狠精品人妻久久久久久综合| 国内揄拍国产精品人妻在线| 久久久久久久久久久丰满| 亚洲av.av天堂| 丝袜喷水一区| 麻豆国产97在线/欧美| 国产av一区二区精品久久 | 国产免费又黄又爽又色| 多毛熟女@视频| 国产精品人妻久久久久久| 99九九线精品视频在线观看视频| 久久鲁丝午夜福利片| 国产免费视频播放在线视频| 免费av不卡在线播放| 18禁裸乳无遮挡免费网站照片| 99视频精品全部免费 在线| 18禁裸乳无遮挡免费网站照片| 日本免费在线观看一区| 久久久色成人| 国产黄色免费在线视频| 久久精品久久久久久久性| 成人黄色视频免费在线看| 欧美极品一区二区三区四区| 精品久久久久久电影网| 国产一区有黄有色的免费视频| 亚洲精品色激情综合| 日韩电影二区| 成人影院久久| 免费看光身美女| 一二三四中文在线观看免费高清| 欧美最新免费一区二区三区| 夜夜骑夜夜射夜夜干| 国产精品一区二区性色av| 国产 精品1| 亚洲精品乱码久久久久久按摩| 国产乱来视频区| 免费观看a级毛片全部| 亚洲人成网站在线播| 国产成人91sexporn| 精品亚洲成国产av| 国国产精品蜜臀av免费| 99re6热这里在线精品视频| 国产av国产精品国产| 夫妻性生交免费视频一级片| 女性被躁到高潮视频| 男女边摸边吃奶| 国产69精品久久久久777片| 一区在线观看完整版| 国产深夜福利视频在线观看| 777米奇影视久久| 午夜视频国产福利| 日韩成人av中文字幕在线观看| 看十八女毛片水多多多| 在现免费观看毛片| 亚洲av国产av综合av卡| 美女视频免费永久观看网站| 久久久久久久久大av| 亚洲国产精品国产精品| 亚洲内射少妇av| 日日摸夜夜添夜夜添av毛片| 国产免费视频播放在线视频| 国产免费福利视频在线观看| 国产精品一区www在线观看| 自拍偷自拍亚洲精品老妇| 午夜精品国产一区二区电影| 国产一区亚洲一区在线观看| 久热这里只有精品99| 天堂8中文在线网| 精品久久国产蜜桃| videos熟女内射| 免费久久久久久久精品成人欧美视频 | 成人毛片60女人毛片免费| 男人爽女人下面视频在线观看| 国产黄色免费在线视频| 中文字幕免费在线视频6| 18禁裸乳无遮挡动漫免费视频| 久久这里有精品视频免费| 日韩精品有码人妻一区| 久久久a久久爽久久v久久| 国国产精品蜜臀av免费| 看免费成人av毛片| 国内精品宾馆在线| 亚洲精品,欧美精品| 嘟嘟电影网在线观看| 亚洲av.av天堂| 午夜免费观看性视频| 天美传媒精品一区二区| 一级毛片黄色毛片免费观看视频| 国产 精品1| 久久国产乱子免费精品| 精品人妻一区二区三区麻豆| 欧美少妇被猛烈插入视频| 亚洲久久久国产精品| 最后的刺客免费高清国语| 成人毛片60女人毛片免费| 韩国高清视频一区二区三区| 久久久久精品久久久久真实原创| 黄色欧美视频在线观看| 视频中文字幕在线观看| 欧美日韩一区二区视频在线观看视频在线| 99热6这里只有精品| 在线观看一区二区三区| 大又大粗又爽又黄少妇毛片口| 国产精品精品国产色婷婷| 亚洲经典国产精华液单| 简卡轻食公司| 欧美bdsm另类| 国产淫语在线视频| av在线播放精品| av.在线天堂| 麻豆精品久久久久久蜜桃| 搡女人真爽免费视频火全软件| 久久久久久久亚洲中文字幕| 偷拍熟女少妇极品色| 亚洲经典国产精华液单| 国产成人aa在线观看| 观看av在线不卡| 99热这里只有是精品50| 一级毛片久久久久久久久女| 国产精品久久久久久久电影| 99热6这里只有精品| 亚洲av成人精品一二三区| 久久婷婷青草| 一二三四中文在线观看免费高清| 两个人的视频大全免费| 欧美日韩亚洲高清精品| 亚洲精品色激情综合| 国产毛片在线视频| 最近手机中文字幕大全| 国产精品成人在线| av在线老鸭窝| 综合色丁香网| 永久网站在线| 免费看av在线观看网站| 各种免费的搞黄视频| 91狼人影院| 亚洲精品国产av成人精品| 蜜臀久久99精品久久宅男| 国产乱人偷精品视频| 狂野欧美激情性xxxx在线观看| 女人久久www免费人成看片| 色婷婷久久久亚洲欧美| 美女主播在线视频| 久久精品熟女亚洲av麻豆精品| 亚洲色图综合在线观看| 边亲边吃奶的免费视频| 欧美极品一区二区三区四区| 日本猛色少妇xxxxx猛交久久| 黄色一级大片看看| 乱系列少妇在线播放| 视频区图区小说| 日韩欧美一区视频在线观看 | 高清毛片免费看| 欧美少妇被猛烈插入视频| 国产黄片美女视频| 亚洲欧美日韩另类电影网站 | 干丝袜人妻中文字幕| 免费不卡的大黄色大毛片视频在线观看| kizo精华| 春色校园在线视频观看| 国产黄频视频在线观看| 国产成人精品一,二区| 嫩草影院新地址| 午夜福利网站1000一区二区三区| 纯流量卡能插随身wifi吗| 国产国拍精品亚洲av在线观看| 黄色配什么色好看| 麻豆成人午夜福利视频| 日韩不卡一区二区三区视频在线| 国产精品国产av在线观看| 丰满少妇做爰视频| 国产亚洲5aaaaa淫片| 观看免费一级毛片| 国产黄片视频在线免费观看| 91久久精品电影网| 久久精品国产a三级三级三级| 日韩强制内射视频| 大陆偷拍与自拍| 99热全是精品| 亚洲精品乱久久久久久| 伊人久久精品亚洲午夜| 简卡轻食公司| 国产伦精品一区二区三区视频9| 日韩在线高清观看一区二区三区| 婷婷色综合www| av免费观看日本| 久久这里有精品视频免费| 国产免费一区二区三区四区乱码| 亚洲精品,欧美精品| 日本欧美视频一区| 卡戴珊不雅视频在线播放| 这个男人来自地球电影免费观看 | 如何舔出高潮| 国国产精品蜜臀av免费| 欧美区成人在线视频| 妹子高潮喷水视频| 日韩av免费高清视频| 欧美最新免费一区二区三区| 色视频www国产| 国产国拍精品亚洲av在线观看| 网址你懂的国产日韩在线| 国产视频内射| 啦啦啦在线观看免费高清www| 久久久久久久久久久免费av| 国产精品国产三级国产专区5o| 国产欧美另类精品又又久久亚洲欧美| 九草在线视频观看| 国产中年淑女户外野战色| 免费不卡的大黄色大毛片视频在线观看| 国产一区二区三区综合在线观看 | 黄片wwwwww| 免费看不卡的av| 亚洲av免费高清在线观看| 国内揄拍国产精品人妻在线| 中国美白少妇内射xxxbb| 五月天丁香电影| 菩萨蛮人人尽说江南好唐韦庄| 伦精品一区二区三区| 久久久久精品性色| 七月丁香在线播放| 日本黄色日本黄色录像| 少妇人妻精品综合一区二区| 99久久中文字幕三级久久日本| 亚洲国产日韩一区二区| 亚洲欧美日韩另类电影网站 | 直男gayav资源| 校园人妻丝袜中文字幕| 男女边摸边吃奶|