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

    基于雙向流固耦合的流線型箱梁斷面顫振分析

    2019-09-10 07:22:44黃新章
    西部交通科技 2019年12期

    摘要:文章以大型有限元分析軟件Ansys - Fluent為計(jì)算平臺(tái),基于計(jì)算流體動(dòng)力學(xué)方法,采用大渦模擬湍流模型,求解不可壓縮流體N-S方程,對(duì)在不同風(fēng)速下蘇通大橋的流線型箱梁斷面的振動(dòng)響應(yīng)進(jìn)行了CFD數(shù)值模擬研究。研究結(jié)果表明:當(dāng)來流風(fēng)速處于145 m/s附近時(shí),發(fā)生了“軟顫振”現(xiàn)象;當(dāng)風(fēng)速超過147 m/s時(shí),發(fā)生了明顯的“顫振”現(xiàn)象,說明該流固耦合計(jì)算方法所得出的數(shù)值模擬結(jié)果與已有文獻(xiàn)試驗(yàn)結(jié)果吻合較好。

    關(guān)鍵詞:箱梁斷面;顫振;動(dòng)網(wǎng)格;UDF;雙向流固耦合

    中圖分類號(hào):U448. 21+3文獻(xiàn)標(biāo)識(shí)碼:A

    DOI: 10.1 3282/j. cnki. wccst. 201 9. 12. 026

    文章編號(hào):1673 - 4874(2019)12 - 0095 - 05

    0引言

    顫振是一種典型的橫風(fēng)向發(fā)散性自激振動(dòng)現(xiàn)象,當(dāng)氣流流經(jīng)流線型斷面時(shí),氣流和結(jié)構(gòu)相互耦合作用。當(dāng)來流風(fēng)速超過某一定值時(shí),振動(dòng)的橋梁結(jié)構(gòu)會(huì)在垂直于氣流的方向發(fā)生非常顯著的大幅度振動(dòng),其特點(diǎn)是振動(dòng)的橋梁結(jié)構(gòu)源源不斷地從氣流中吸收能量,此時(shí)結(jié)構(gòu)吸收的能量大于結(jié)構(gòu)阻尼所消耗的能量。隨著結(jié)構(gòu)設(shè)計(jì)和施工水平的提高,現(xiàn)代橋粱結(jié)構(gòu)向柔性、輕質(zhì)的方向發(fā)展,這使得橋梁對(duì)風(fēng)荷載的敏感性大大增強(qiáng),特別是在施工狀態(tài),風(fēng)致振動(dòng)問題表現(xiàn)得更為激烈。因此,施工狀態(tài)下的風(fēng)致振動(dòng)問題是橋梁風(fēng)致作用研究的主要方向。起初,橋梁抗風(fēng)動(dòng)力特性研究主要是通過縮尺模型風(fēng)洞試驗(yàn)進(jìn)行分析,實(shí)驗(yàn)往往受到人力、物力、時(shí)間等因素的影響,使得橋梁結(jié)構(gòu)風(fēng)洞試驗(yàn)研究舉步維艱。隨著計(jì)算機(jī)硬件的提升以及計(jì)算流體動(dòng)力學(xué)理論的完善,CFD數(shù)值模擬逐步取代了風(fēng)洞試驗(yàn)。

    目前,大跨度橋梁的抗風(fēng)研究主要集中在渦振和抖振,對(duì)流線型箱梁斷面的顫振研究相對(duì)較少。Bleich[1]首次運(yùn)用理想平板氣動(dòng)自激力理論研究了橋梁扭轉(zhuǎn)顫振的可能性。Tomko[2]對(duì)主梁斷面顫振導(dǎo)數(shù)進(jìn)行了深入研究,奠定了橋梁顫振分析理論的基礎(chǔ)。Frandsen[3]基于任意拉格朗日一歐拉法對(duì)彈性支撐于流場(chǎng)中的橋梁斷面進(jìn)行了流固耦合數(shù)值模擬研究,并測(cè)出了顫振臨界風(fēng)速。Slevam[4-6]等通過有限元方法求解流體瞬態(tài)N-S方程以及通過ALE方法實(shí)現(xiàn)了橋梁結(jié)構(gòu)流固耦合數(shù)值計(jì)算模型,最后對(duì)其顫振穩(wěn)定性計(jì)算進(jìn)行了流固耦合數(shù)值模擬研究。劉志文[7]通過在Fluent嵌入自定義程序代碼UDF,基于Newmark-p法求解結(jié)構(gòu)振動(dòng)方程,采用兩方程模型求解不可壓縮Navier -Stokes方程,利用彈簧光順法及網(wǎng)格重構(gòu)法搭配使用的動(dòng)網(wǎng)格技術(shù)建立了流固耦合數(shù)值計(jì)算模型,對(duì)橋梁斷面進(jìn)行了風(fēng)致振動(dòng)響應(yīng)計(jì)算。李永樂[8]基于CSD和CFD耦合振動(dòng)分析方法,采用微分方程的數(shù)值解法和動(dòng)網(wǎng)格技術(shù)對(duì)方柱和平板的顫振進(jìn)行了數(shù)值模擬研究。張瑞琴[9]利用CFX軟件對(duì)典型翼型斷面進(jìn)行了流固耦合作用下的顫振研究。綜上所述,當(dāng)前的文獻(xiàn)基于CFD方法對(duì)橋梁斷面的顫振研究依然較少。

    本文以蘇通大橋的流線型箱梁斷面為研究對(duì)象,采用“剛生運(yùn)動(dòng)區(qū)域+動(dòng)網(wǎng)格區(qū)域+靜網(wǎng)格區(qū)域”的動(dòng)網(wǎng)格策略建立網(wǎng)格,基于Newmark -β法,對(duì)計(jì)算流體動(dòng)力學(xué)軟件Fluent進(jìn)行二次開發(fā),將求解的結(jié)構(gòu)振動(dòng)響應(yīng)編入U(xiǎn)DF中,對(duì)橋斷面風(fēng)致振動(dòng)響應(yīng)進(jìn)行了流固耦合數(shù)值模擬研究。本文數(shù)值模擬得到了考慮流固耦合作用時(shí)的流場(chǎng)特性,并分析了顫振響應(yīng)機(jī)理。本文數(shù)值模擬結(jié)果與已有文獻(xiàn)的結(jié)果吻合較好[10]。

    1 箱梁斷面流固耦合數(shù)值模擬方法

    1.1 流體控制方程求解

    CFD流體計(jì)算為了準(zhǔn)確地對(duì)湍流流動(dòng)進(jìn)行數(shù)值模擬,一方面要求計(jì)算域的尺寸應(yīng)該包含湍流中可能出現(xiàn)的最大渦,另外一方面要求網(wǎng)格最小尺寸足以分辨最小渦的運(yùn)動(dòng)。然而,就目前的計(jì)算機(jī)硬件能力而言,能夠采用的最小的計(jì)算網(wǎng)格的尺度仍然比最小渦的尺度大許多。權(quán)衡之下,F(xiàn)luent求解器中推出了一種比較流行的湍流模型——大渦模擬LES湍流模型,可以通過建立亞格子尺度模型來模擬小尺度渦對(duì)大尺度渦運(yùn)動(dòng)的影響。

    本文CFD數(shù)值模擬分析流體計(jì)算采用二維大渦模擬湍流模型進(jìn)行求解不可壓縮的Navier - Stokes方程,其控制方程為:

    1.2 結(jié)構(gòu)振動(dòng)控制方程求解

    需要注意的是,Ansys - Fluent僅僅是CFD流體計(jì)算分析軟件,并不能直接對(duì)結(jié)構(gòu)的變形進(jìn)行求解。因此,當(dāng)我們需要求解結(jié)構(gòu)的振動(dòng)響應(yīng)時(shí),就需要編寫相應(yīng)的振動(dòng)方程的程序。本文通過編寫自定義程序代碼UDF對(duì)Fluent進(jìn)行二次開發(fā),得到了一種箱梁斷面風(fēng)致振動(dòng)流固耦合的數(shù)值分析方法。此程序還實(shí)現(xiàn)了Fluent并行計(jì)算,提高了計(jì)算效率,并且可以對(duì)類似的二維結(jié)構(gòu)斷面進(jìn)行渦振、顫振分析,對(duì)其結(jié)構(gòu)安全性具有重要的評(píng)估意義。

    本文考慮了結(jié)構(gòu)的豎向振動(dòng)及扭轉(zhuǎn)振動(dòng),即對(duì)蘇通大橋流線型箱梁斷面進(jìn)行了兩自由度的流固耦合數(shù)值模擬研究。箱梁斷面振動(dòng)控制方程表達(dá)式為:

    本文是將Newmark -口算法嵌入U(xiǎn)DF中求解結(jié)構(gòu)的振動(dòng)方程,通過Define—CG—Motion指定動(dòng)網(wǎng)格區(qū)域的網(wǎng)格運(yùn)動(dòng)位移,進(jìn)而實(shí)現(xiàn)氣流與蘇通大橋箱梁斷面的流固耦合數(shù)值模擬方法。

    2 工程背景及計(jì)算參數(shù)

    2.1 工程背景

    為了檢驗(yàn)本文流固耦合數(shù)值模擬方法的計(jì)算精度,進(jìn)行了寬高比為10:1的簡(jiǎn)單橋斷面顫振臨界風(fēng)速計(jì)算。以蘇通大橋?yàn)橛?jì)算背景,該橋?yàn)殡p塔雙索面斜拉橋,主橋跨徑采用100 m+100 m+ 300 m+1 080m+300 m+100 m+ 100 m,主梁采用的是流線型鋼箱梁斷面型式,其橫截面橋?qū)?1 m,梁高4m。

    2.2 計(jì)算分析參數(shù)

    2.2.1 蘇通大橋流線型箱梁斷面動(dòng)力特性參數(shù)

    本文采用大型有限元分析軟件Ansys對(duì)蘇通大橋進(jìn)行了動(dòng)力特性分析,得到了豎彎頻率及扭轉(zhuǎn)頻率,分別為0. 196 Hz、0.565 Hz,如表1所示。

    2.2.2邊界條件及湍流模型

    入口采用速度入口邊界條件,切向速度為零,湍流強(qiáng)度設(shè)為0.5%,湍流黏性比設(shè)為2;出口為壓力出口邊界條件,參考?jí)毫橐粋€(gè)大氣壓。橋梁主梁斷面表面為無滑移壁面,整個(gè)流域阻塞率為0.3%,滿足阻塞率<5%的要求。計(jì)算域分塊及邊界條件如圖2所示。

    本文CFD數(shù)值模擬計(jì)算結(jié)合“剛性運(yùn)動(dòng)區(qū)域+動(dòng)網(wǎng)格運(yùn)動(dòng)區(qū)域+靜止網(wǎng)格運(yùn)動(dòng)區(qū)域”的動(dòng)網(wǎng)格劃分策略,箱梁斷面壁面處首層網(wǎng)格高度取值為1/4 0008,其中B=41 m。剛性運(yùn)動(dòng)區(qū)域采用高質(zhì)量全結(jié)構(gòu)化網(wǎng)格,將箱梁結(jié)構(gòu)壁面隨剛生層運(yùn)動(dòng)區(qū)域一起運(yùn)動(dòng),從而解決了箱梁壁面邊界層網(wǎng)格因?yàn)樽冃味霈F(xiàn)的網(wǎng)格畸變和負(fù)體積等問題,確保了邊界層關(guān)鍵部位的網(wǎng)格質(zhì)量,提高了計(jì)算精度。動(dòng)網(wǎng)格區(qū)域采用三角形網(wǎng)格進(jìn)行劃分,在結(jié)構(gòu)振動(dòng)過程中,這部分網(wǎng)格會(huì)進(jìn)行重新劃分,更新網(wǎng)格。為了提高計(jì)算效率,最外側(cè)網(wǎng)格采用靜止的結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分。

    本文湍流模型為大渦模擬LES三維湍流模型,先進(jìn)行2 000步穩(wěn)態(tài)計(jì)算,將穩(wěn)態(tài)結(jié)果作為瞬態(tài)計(jì)算的初場(chǎng)以提高計(jì)算效率。為了滿足計(jì)算結(jié)果的精度,壁面網(wǎng)格無量綱高度Y- plus值在1附近。經(jīng)過網(wǎng)格無關(guān)性和時(shí)間無關(guān)性測(cè)試,最終時(shí)間步長(zhǎng)確定為5×10 s,流體分析二維總網(wǎng)格數(shù)量為21萬,網(wǎng)格劃分及箱梁壁面Y- plus值如圖3所示。計(jì)算在國(guó)家超算高性能計(jì)算機(jī)群上進(jìn)行,每個(gè)工況申請(qǐng)6個(gè)CPU進(jìn)行并行計(jì)算,單工況耗時(shí)約為50 CPUs。

    3 數(shù)值模擬結(jié)果及分析

    文獻(xiàn)M的作者給出了蘇通大橋的顫振臨界風(fēng)速為145.5 m/s,為此,本文采用CFD流固耦合數(shù)值模擬方法對(duì)3個(gè)來流風(fēng)速進(jìn)行了流固耦合顫振響應(yīng)研究,分別為144 m/s、146 m/s、147 m/s。Fluent計(jì)算過程中,采用了SST k-w、RNG k-e和LES大渦模擬共計(jì)三種湍流模型分別對(duì)蘇通大橋計(jì)算了其箱梁斷面的三分力系數(shù),選用的時(shí)間步長(zhǎng)有0.000 1 s、0.000 5 s、0.000 05 s。最后比較得出:LES大渦模擬湍流模型及0.000 5 s的時(shí)間步長(zhǎng)模擬出的結(jié)果和參考文獻(xiàn)的結(jié)果最為接近。因此,本文的流固耦合數(shù)值模擬部分最終采用的是LES大渦模擬湍流模型,數(shù)值模擬結(jié)果與參考文獻(xiàn)的結(jié)果相對(duì)吻合。

    3.1 振動(dòng)響應(yīng)分析

    下頁圖4~6為蘇通大橋鋼箱梁斷面的風(fēng)致振動(dòng)時(shí)程位移隨來流風(fēng)速的變化曲線。當(dāng)來流風(fēng)速≤146 m/s時(shí),豎彎振動(dòng)和扭轉(zhuǎn)振動(dòng)幅值均很小,豎彎運(yùn)動(dòng)和扭轉(zhuǎn)運(yùn)動(dòng)均呈現(xiàn)出衰減運(yùn)動(dòng)的趨勢(shì),最終穩(wěn)定下來,上下波動(dòng)特別緩和;當(dāng)來流風(fēng)速達(dá)到146 m/s時(shí),此時(shí)的蘇通大橋發(fā)生了非常明顯的大幅度振動(dòng),扭轉(zhuǎn)位移發(fā)散,出現(xiàn)了“顫振”現(xiàn)象。經(jīng)過參考大量的關(guān)于顫振的文獻(xiàn)發(fā)現(xiàn),當(dāng)橋梁進(jìn)入顫振臨界狀態(tài)后,不一定一直呈現(xiàn)出發(fā)散性的運(yùn)動(dòng),而是呈現(xiàn)出穩(wěn)定的大振幅“扭轉(zhuǎn)渦激共振”現(xiàn)象,即“軟顫振”現(xiàn)象。比如當(dāng)來流風(fēng)速為147 m/s時(shí),出現(xiàn)了“軟顫振”現(xiàn)象,即“扭轉(zhuǎn)渦激共振”現(xiàn)象。

    3.2 振動(dòng)機(jī)理分析

    考慮流固耦合效應(yīng)下的風(fēng)場(chǎng)不同時(shí)刻處的風(fēng)速云圖(如圖7所示)。從風(fēng)速云圖可以明顯看出:不同時(shí)刻的漩渦脫落位置均不一樣,漩渦表現(xiàn)得非常激烈,同時(shí)夾雜了許多小渦。從這里也側(cè)面反映出蘇通大橋在高風(fēng)速下,發(fā)生了大幅度振動(dòng)。

    4 結(jié)語

    本文基于CFD方法,以蘇通大橋的流線型箱梁斷面為研究對(duì)象,成功模擬了高風(fēng)速下蘇通大橋顫振臨界風(fēng)速,結(jié)論如下:

    (1)本文在Fluent軟件的基礎(chǔ)上進(jìn)行二次開發(fā),結(jié)合“剛、生運(yùn)動(dòng)區(qū)域+動(dòng)網(wǎng)格區(qū)域+靜網(wǎng)格區(qū)域”的網(wǎng)格劃分策略,有效解決了蘇通大橋顫振響應(yīng)計(jì)算過程中由于大變形可能產(chǎn)生的網(wǎng)格畸變及負(fù)體積等問題。經(jīng)數(shù)值模擬結(jié)果與相關(guān)文獻(xiàn)進(jìn)行比較,證明了該方法的正確性,此方法可應(yīng)用于類似結(jié)構(gòu)的流固耦合分析。

    (2)本文通過CFD數(shù)值模擬分析,得到了蘇通大橋的顫振臨界風(fēng)速為146 m/s。

    (3)當(dāng)風(fēng)速增大到147 m/s時(shí),出現(xiàn)了大幅值穩(wěn)定的扭轉(zhuǎn)運(yùn)動(dòng),即發(fā)生了“軟顫振”現(xiàn)象。

    參考文獻(xiàn)

    [1] Bleich F_ Dynamic instability of truss-stiffened suspensionbridges urder wind action[J].Transocticns of the Americanscciety of civil engineers,1949( 114) : 1 177-1 222.

    [2]Scanlan R H, Tomko J J.Airtoil and Bridge Deck FlutterDerivatives[J]. Journalof the Engineering MechanicsDivision,1971 , 97( 6) :1 717 -1 733.

    [3]Frandsen J B.Numerical bridge deck studies using finiteelements Part l: flutter [J]. Journal of Fluids &Structures, 2004,19( 2) : 1 71 -191.

    [4] Bosch H, Selvam R P. Govindawamy S. Aeroelasticanalysis of bridges using FEM and moving grids[J].Wind & Structures An International Journal, 2002. 5:257 - 266.

    [5] Dettmer W, Peri D. A computational framewo rk for f|uid -rigid body interaction: Finite element formulation andapplications [ J ]. Computer Methods in Applied Mechanics and Engineering, 2006 ,195(13 - 1 6) :1 633 -1 666.

    [6]Braun A L, Awruch A M. Finite element simulation of thewind action over bridge sectional modelS: Applicationto the Guama River Bridge( Par State, Brazil)[J]. FiniteElements in Analysis&Design, 2008, 44(3) :105 - 122.

    [7]劉志文,周帥,陳政清.寬高比為4的矩形斷面渦激振動(dòng)響應(yīng)數(shù)值模擬[J].振動(dòng)與j中擊,2011,30(11):153 - 156.

    [8]李永樂.基于CFD的橋梁渦激振和顫振氣彈模擬研究[D].成都:西南交通大學(xué),2013.

    [9]張瑞琴,翁建生.基于流固耦合的葉片顫振分析[J].計(jì)算機(jī)仿真,2011(3):48- 51.

    [10] Rui Zhou, Yaojun Ge, Yongxin Yang, et al. Quasi-3Dfinite element modelling of nonlinear aerodynamic be-haviours of super long-span cable-supported bridges[J].Transactions of the American society of civil engi-neers, 2018(52):1 011 -1 040.

    [11]王林凱,夏飛龍,劉志文,等.基于流固耦合的鋼箱梁斷面顫振數(shù)值模擬[J].公路,2019, 64(7):108 - 114.

    作者簡(jiǎn)介:黃新章(1976-),工程師,從事國(guó)省道日常養(yǎng)護(hù)及項(xiàng)目大中修管理工作。

    日本黄大片高清| 神马国产精品三级电影在线观看| avwww免费| 有码 亚洲区| 制服丝袜大香蕉在线| 成人性生交大片免费视频hd| 精品国内亚洲2022精品成人| 两个人的视频大全免费| 精品乱码久久久久久99久播| 国产精华一区二区三区| 国产综合懂色| 成人国产综合亚洲| 久久精品国产自在天天线| 国产一区二区激情短视频| 中国美女看黄片| 蜜桃亚洲精品一区二区三区| 久久久久免费精品人妻一区二区| 亚洲av成人不卡在线观看播放网| 欧美性猛交黑人性爽| 免费观看的影片在线观看| 婷婷丁香在线五月| 黄色配什么色好看| 亚洲精品日韩av片在线观看| 一区二区三区四区激情视频 | av福利片在线观看| 久久这里只有精品中国| 毛片一级片免费看久久久久 | 久久久国产成人精品二区| 麻豆一二三区av精品| 午夜福利在线观看免费完整高清在 | 深爱激情五月婷婷| 99久久精品国产亚洲精品| 天堂网av新在线| 日韩中文字幕欧美一区二区| 国产精品一区二区三区四区久久| 1024手机看黄色片| 日日摸夜夜添夜夜添小说| 日本黄色视频三级网站网址| 久久久久久久久久黄片| 一本一本综合久久| 欧美最黄视频在线播放免费| 日韩免费av在线播放| 亚洲精品色激情综合| 国模一区二区三区四区视频| 在线播放国产精品三级| 亚洲精品一卡2卡三卡4卡5卡| 成人三级黄色视频| 色播亚洲综合网| 久久精品国产99精品国产亚洲性色| 日韩欧美在线乱码| 亚洲成人精品中文字幕电影| 成人午夜高清在线视频| 欧美高清性xxxxhd video| 日本熟妇午夜| 婷婷六月久久综合丁香| 性色avwww在线观看| 国产精品永久免费网站| 一区二区三区激情视频| 欧美精品国产亚洲| 国产精品美女特级片免费视频播放器| av专区在线播放| 三级国产精品欧美在线观看| 欧美性猛交黑人性爽| 简卡轻食公司| 在线观看舔阴道视频| 日韩欧美国产一区二区入口| 日韩欧美免费精品| 欧美成人免费av一区二区三区| 亚洲国产日韩欧美精品在线观看| 一级a爱片免费观看的视频| 欧美日韩乱码在线| 欧美国产日韩亚洲一区| 午夜福利在线观看免费完整高清在 | 成人鲁丝片一二三区免费| 国产精品不卡视频一区二区 | 欧美另类亚洲清纯唯美| 麻豆成人av在线观看| eeuss影院久久| 91久久精品国产一区二区成人| 在线免费观看不下载黄p国产 | 此物有八面人人有两片| 免费大片18禁| 日韩欧美在线二视频| .国产精品久久| 国产伦精品一区二区三区视频9| 亚洲国产精品成人综合色| 特大巨黑吊av在线直播| 午夜两性在线视频| 久久精品国产清高在天天线| 99在线人妻在线中文字幕| 亚洲av美国av| 亚洲成人精品中文字幕电影| 国产精品影院久久| 深夜a级毛片| 国产主播在线观看一区二区| 宅男免费午夜| 性色av乱码一区二区三区2| 亚洲五月天丁香| 可以在线观看的亚洲视频| 国产精品久久久久久久久免 | 不卡一级毛片| 五月伊人婷婷丁香| 亚洲自偷自拍三级| 观看美女的网站| 日韩有码中文字幕| 国产成人欧美在线观看| 国产精品一区二区三区四区免费观看 | 亚洲午夜理论影院| 国产一区二区在线观看日韩| 青草久久国产| 国内揄拍国产精品人妻在线| www.www免费av| 最近最新中文字幕大全电影3| 国产激情偷乱视频一区二区| 国产精品久久久久久亚洲av鲁大| 欧美激情在线99| 亚洲第一电影网av| 日韩欧美 国产精品| 亚洲熟妇熟女久久| 夜夜夜夜夜久久久久| 熟妇人妻久久中文字幕3abv| 亚洲激情在线av| 性色avwww在线观看| 午夜精品一区二区三区免费看| 日韩免费av在线播放| 国产精品精品国产色婷婷| 成年女人毛片免费观看观看9| 人妻制服诱惑在线中文字幕| 欧美黄色片欧美黄色片| 桃色一区二区三区在线观看| 国产老妇女一区| 黄片小视频在线播放| 久9热在线精品视频| a级一级毛片免费在线观看| 国产精品影院久久| 欧美xxxx黑人xx丫x性爽| 观看免费一级毛片| 久久欧美精品欧美久久欧美| 啪啪无遮挡十八禁网站| 成人欧美大片| www日本黄色视频网| 老鸭窝网址在线观看| 国产精品一区二区免费欧美| 亚洲精品一区av在线观看| 国产伦人伦偷精品视频| 九色国产91popny在线| 一区二区三区高清视频在线| 99国产精品一区二区蜜桃av| 国产欧美日韩精品一区二区| 自拍偷自拍亚洲精品老妇| 狠狠狠狠99中文字幕| 亚洲第一电影网av| 国产av在哪里看| 两个人视频免费观看高清| 国产av一区在线观看免费| 欧美日韩黄片免| 淫秽高清视频在线观看| 黄色丝袜av网址大全| av在线蜜桃| 国产一区二区三区在线臀色熟女| 丰满人妻熟妇乱又伦精品不卡| 亚洲成人久久性| 嫩草影院精品99| 亚洲av不卡在线观看| 国产精品亚洲一级av第二区| 一级毛片久久久久久久久女| 亚洲欧美日韩高清在线视频| 亚洲成人久久性| 免费观看人在逋| 国内久久婷婷六月综合欲色啪| 三级国产精品欧美在线观看| 国产精品久久久久久人妻精品电影| 亚洲人成网站高清观看| 久久精品国产自在天天线| 给我免费播放毛片高清在线观看| 赤兔流量卡办理| www日本黄色视频网| 美女cb高潮喷水在线观看| 18+在线观看网站| 97人妻精品一区二区三区麻豆| 久久久久国内视频| 一个人看视频在线观看www免费| 国产精品一及| 无人区码免费观看不卡| 免费av不卡在线播放| 九色成人免费人妻av| 男人狂女人下面高潮的视频| 12—13女人毛片做爰片一| 免费大片18禁| 亚洲av成人精品一区久久| 亚洲乱码一区二区免费版| 性插视频无遮挡在线免费观看| av专区在线播放| eeuss影院久久| 1024手机看黄色片| 国产精品1区2区在线观看.| 国产av不卡久久| 成年女人毛片免费观看观看9| 日韩中文字幕欧美一区二区| 日韩av在线大香蕉| 日韩大尺度精品在线看网址| 亚洲成人中文字幕在线播放| 美女cb高潮喷水在线观看| 国产熟女xx| 91麻豆精品激情在线观看国产| 亚洲人与动物交配视频| 国产精品亚洲一级av第二区| 日本a在线网址| 国内精品久久久久精免费| 国产男靠女视频免费网站| 欧美在线一区亚洲| 免费av观看视频| 婷婷六月久久综合丁香| 欧美中文日本在线观看视频| 亚洲专区中文字幕在线| 国产精品久久视频播放| 国产精品乱码一区二三区的特点| www.999成人在线观看| 最新在线观看一区二区三区| 夜夜看夜夜爽夜夜摸| 欧美性猛交╳xxx乱大交人| 乱码一卡2卡4卡精品| 久久人人精品亚洲av| 禁无遮挡网站| 伦理电影大哥的女人| 变态另类丝袜制服| av黄色大香蕉| 亚洲成人中文字幕在线播放| 久久久久性生活片| 在线观看免费视频日本深夜| 岛国在线免费视频观看| 中文字幕免费在线视频6| 国产伦精品一区二区三区四那| 国产探花在线观看一区二区| 一个人看视频在线观看www免费| 日本五十路高清| 极品教师在线免费播放| 色视频www国产| 国产精品不卡视频一区二区 | 亚洲人成伊人成综合网2020| 日日摸夜夜添夜夜添av毛片 | 欧美午夜高清在线| 嫁个100分男人电影在线观看| 九色国产91popny在线| 三级国产精品欧美在线观看| 90打野战视频偷拍视频| 欧美中文日本在线观看视频| 国产白丝娇喘喷水9色精品| 亚洲自拍偷在线| 真人一进一出gif抽搐免费| 黄色配什么色好看| 日韩欧美免费精品| 日韩大尺度精品在线看网址| 黄片小视频在线播放| 91麻豆av在线| 日韩欧美精品v在线| 1024手机看黄色片| 欧美三级亚洲精品| 最近最新免费中文字幕在线| 久久久久久久亚洲中文字幕 | 又爽又黄a免费视频| 久久久久亚洲av毛片大全| 青草久久国产| 三级男女做爰猛烈吃奶摸视频| 俄罗斯特黄特色一大片| 在线国产一区二区在线| 精品久久久久久久久av| 淫妇啪啪啪对白视频| 男人舔女人下体高潮全视频| 蜜桃亚洲精品一区二区三区| 日韩精品青青久久久久久| 最后的刺客免费高清国语| 两性午夜刺激爽爽歪歪视频在线观看| 丰满乱子伦码专区| 成人一区二区视频在线观看| 久久久久久久久久黄片| 中文字幕人妻熟人妻熟丝袜美| 亚洲乱码一区二区免费版| 午夜精品久久久久久毛片777| 国产乱人伦免费视频| 国产精品精品国产色婷婷| 亚洲三级黄色毛片| 深夜精品福利| 九色国产91popny在线| 久久精品久久久久久噜噜老黄 | 老司机深夜福利视频在线观看| 亚州av有码| 欧洲精品卡2卡3卡4卡5卡区| 搡老岳熟女国产| 18禁黄网站禁片午夜丰满| 一区二区三区激情视频| a在线观看视频网站| 国产精品一区二区三区四区免费观看 | 国产精品综合久久久久久久免费| 一本一本综合久久| 国产一级毛片七仙女欲春2| 亚洲五月天丁香| 天天一区二区日本电影三级| 制服丝袜大香蕉在线| 成人精品一区二区免费| 免费高清视频大片| 高潮久久久久久久久久久不卡| 国产主播在线观看一区二区| 亚洲美女搞黄在线观看 | 中文字幕av成人在线电影| 给我免费播放毛片高清在线观看| 成人特级黄色片久久久久久久| 岛国在线免费视频观看| 欧美日韩瑟瑟在线播放| 69人妻影院| 亚洲无线观看免费| 欧美性猛交╳xxx乱大交人| 国产av麻豆久久久久久久| 国产91精品成人一区二区三区| 色综合欧美亚洲国产小说| 亚洲精品亚洲一区二区| 精品无人区乱码1区二区| 99久久久亚洲精品蜜臀av| 在线观看午夜福利视频| 国产精品98久久久久久宅男小说| 日本 欧美在线| 亚洲va日本ⅴa欧美va伊人久久| 综合色av麻豆| 中文字幕高清在线视频| 琪琪午夜伦伦电影理论片6080| 精品久久久久久久久久免费视频| 美女高潮喷水抽搐中文字幕| 真实男女啪啪啪动态图| 精品乱码久久久久久99久播| 久久久久久久亚洲中文字幕 | 少妇的逼水好多| 精品午夜福利在线看| 老女人水多毛片| 色综合亚洲欧美另类图片| 亚洲av不卡在线观看| 深夜a级毛片| 夜夜夜夜夜久久久久| 久久久久亚洲av毛片大全| 国产成人aa在线观看| 亚洲一区二区三区色噜噜| 国产精华一区二区三区| 999久久久精品免费观看国产| 日韩欧美精品免费久久 | 亚洲av免费高清在线观看| 亚洲成人免费电影在线观看| 俺也久久电影网| 国产黄片美女视频| 午夜福利在线观看吧| 亚洲经典国产精华液单 | 永久网站在线| 男人和女人高潮做爰伦理| 国产毛片a区久久久久| 日韩欧美三级三区| 99久久成人亚洲精品观看| 欧美bdsm另类| 少妇人妻精品综合一区二区 | 制服丝袜大香蕉在线| 搡老熟女国产l中国老女人| 亚洲狠狠婷婷综合久久图片| 欧美日韩国产亚洲二区| 日本a在线网址| 日本熟妇午夜| 久久亚洲真实| 午夜两性在线视频| 真人一进一出gif抽搐免费| 色哟哟哟哟哟哟| 桃红色精品国产亚洲av| 国产免费一级a男人的天堂| 毛片女人毛片| 成人av一区二区三区在线看| 少妇人妻精品综合一区二区 | 精品一区二区三区视频在线| 日韩人妻高清精品专区| 免费在线观看亚洲国产| 午夜福利18| 国产单亲对白刺激| 国内揄拍国产精品人妻在线| 日韩欧美国产在线观看| 97超视频在线观看视频| 免费人成在线观看视频色| 亚洲精品粉嫩美女一区| 国产精品人妻久久久久久| 久久久久久久午夜电影| 国产伦在线观看视频一区| 99国产综合亚洲精品| 国产伦在线观看视频一区| 日本免费a在线| 亚洲,欧美,日韩| 亚洲av第一区精品v没综合| 99国产精品一区二区三区| 69人妻影院| 少妇高潮的动态图| 久久午夜福利片| 三级男女做爰猛烈吃奶摸视频| 最后的刺客免费高清国语| 一个人看的www免费观看视频| 18禁在线播放成人免费| 一级av片app| 中文字幕av在线有码专区| 日韩欧美国产一区二区入口| 亚洲av熟女| a级一级毛片免费在线观看| 国产欧美日韩一区二区三| 国产视频内射| 国产免费一级a男人的天堂| 性插视频无遮挡在线免费观看| 国产精品99久久久久久久久| 国产欧美日韩精品亚洲av| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美日韩无卡精品| 成年人黄色毛片网站| 看十八女毛片水多多多| 一个人免费在线观看的高清视频| 色尼玛亚洲综合影院| 精品久久久久久,| 日韩 亚洲 欧美在线| 国产午夜福利久久久久久| 日韩av在线大香蕉| 国产午夜精品论理片| 国产一区二区激情短视频| av在线老鸭窝| 亚洲片人在线观看| 最新在线观看一区二区三区| avwww免费| 一本精品99久久精品77| 麻豆国产97在线/欧美| 波多野结衣巨乳人妻| av在线天堂中文字幕| 757午夜福利合集在线观看| 老鸭窝网址在线观看| 国产亚洲精品久久久久久毛片| 青草久久国产| 国产精品影院久久| 亚洲男人的天堂狠狠| 少妇丰满av| 久久久久精品国产欧美久久久| 88av欧美| 亚洲午夜理论影院| 欧美在线一区亚洲| 999久久久精品免费观看国产| 男人狂女人下面高潮的视频| www日本黄色视频网| www.熟女人妻精品国产| 久久人人爽人人爽人人片va | 国产91精品成人一区二区三区| 国产探花在线观看一区二区| 淫秽高清视频在线观看| 国产精品国产高清国产av| 日韩高清综合在线| 又粗又爽又猛毛片免费看| 男人狂女人下面高潮的视频| 一进一出好大好爽视频| 久久久久久久亚洲中文字幕 | 欧美日韩瑟瑟在线播放| 男人的好看免费观看在线视频| 精品一区二区三区人妻视频| 亚洲成人久久爱视频| 无遮挡黄片免费观看| 亚洲久久久久久中文字幕| 美女黄网站色视频| 午夜视频国产福利| av视频在线观看入口| 午夜两性在线视频| 非洲黑人性xxxx精品又粗又长| 亚洲真实伦在线观看| 精品99又大又爽又粗少妇毛片 | 欧美色视频一区免费| av天堂在线播放| 中文字幕免费在线视频6| 丰满人妻一区二区三区视频av| av专区在线播放| 观看美女的网站| 三级男女做爰猛烈吃奶摸视频| 成年人黄色毛片网站| 国产精品影院久久| 亚洲国产精品久久男人天堂| 床上黄色一级片| 亚洲一区二区三区不卡视频| 免费搜索国产男女视频| 动漫黄色视频在线观看| 国产日本99.免费观看| 在线观看一区二区三区| 夜夜躁狠狠躁天天躁| 我要看日韩黄色一级片| 欧美最新免费一区二区三区 | av天堂中文字幕网| 亚洲国产色片| 日韩成人在线观看一区二区三区| 国产又黄又爽又无遮挡在线| 久久久精品欧美日韩精品| 欧美激情久久久久久爽电影| 免费观看的影片在线观看| 人人妻,人人澡人人爽秒播| 亚洲avbb在线观看| 床上黄色一级片| 一个人看视频在线观看www免费| www.www免费av| 男人舔奶头视频| 国产精品久久电影中文字幕| 亚洲熟妇熟女久久| 欧美日韩瑟瑟在线播放| 久久久色成人| 一级a爱片免费观看的视频| 男女那种视频在线观看| 偷拍熟女少妇极品色| 女同久久另类99精品国产91| 少妇人妻精品综合一区二区 | 亚洲精品456在线播放app | 欧美色视频一区免费| 亚洲欧美日韩无卡精品| 99热这里只有是精品50| 免费看美女性在线毛片视频| 美女免费视频网站| 欧美极品一区二区三区四区| 男女床上黄色一级片免费看| 好看av亚洲va欧美ⅴa在| 日韩欧美精品免费久久 | 精品久久久久久成人av| 精品久久国产蜜桃| 欧美不卡视频在线免费观看| 能在线免费观看的黄片| av欧美777| 日日摸夜夜添夜夜添小说| 久久久久久久精品吃奶| 日韩免费av在线播放| 午夜福利成人在线免费观看| 天堂动漫精品| 亚洲,欧美,日韩| 免费观看的影片在线观看| 国产野战对白在线观看| 日本一本二区三区精品| 偷拍熟女少妇极品色| 亚洲中文字幕日韩| 欧美中文日本在线观看视频| 日本精品一区二区三区蜜桃| 欧美+日韩+精品| 给我免费播放毛片高清在线观看| 亚洲国产精品999在线| 亚洲在线观看片| 51午夜福利影视在线观看| 身体一侧抽搐| 亚洲av电影不卡..在线观看| 九色成人免费人妻av| 婷婷亚洲欧美| 国产麻豆成人av免费视频| 男女做爰动态图高潮gif福利片| 美女被艹到高潮喷水动态| 亚洲黑人精品在线| 一进一出抽搐gif免费好疼| 国产v大片淫在线免费观看| 长腿黑丝高跟| 黄色一级大片看看| 韩国av一区二区三区四区| 99riav亚洲国产免费| 色在线成人网| 亚洲三级黄色毛片| 免费人成视频x8x8入口观看| 三级国产精品欧美在线观看| 中文字幕免费在线视频6| 国产一区二区亚洲精品在线观看| 白带黄色成豆腐渣| 无人区码免费观看不卡| 我的女老师完整版在线观看| 亚洲中文日韩欧美视频| 97热精品久久久久久| 久久欧美精品欧美久久欧美| 乱码一卡2卡4卡精品| 欧美性猛交╳xxx乱大交人| 一区二区三区激情视频| 如何舔出高潮| 久久久久久久午夜电影| 女同久久另类99精品国产91| 性欧美人与动物交配| 亚洲熟妇中文字幕五十中出| 亚洲av五月六月丁香网| 在线看三级毛片| 国产精品免费一区二区三区在线| 欧美乱色亚洲激情| 好男人在线观看高清免费视频| 欧美bdsm另类| 国产精品人妻久久久久久| 亚洲欧美清纯卡通| www日本黄色视频网| 成人精品一区二区免费| 嫁个100分男人电影在线观看| 97人妻精品一区二区三区麻豆| 国产欧美日韩精品一区二区| 成人无遮挡网站| 日本与韩国留学比较| 成人美女网站在线观看视频| 欧美在线一区亚洲| 国产av在哪里看| 国内精品久久久久精免费| 老鸭窝网址在线观看| 国产三级在线视频| 亚洲在线观看片| 啪啪无遮挡十八禁网站| 2021天堂中文幕一二区在线观| 亚洲欧美日韩高清在线视频| 欧美成人a在线观看| 午夜免费成人在线视频| 怎么达到女性高潮| 亚洲精品影视一区二区三区av| 免费看a级黄色片| 午夜精品一区二区三区免费看| 深爱激情五月婷婷| 国产精品98久久久久久宅男小说| 亚洲欧美日韩高清专用| 亚洲精品在线观看二区| 老司机午夜福利在线观看视频| 久久久久久大精品| 国产毛片a区久久久久| 免费人成视频x8x8入口观看| 成人特级黄色片久久久久久久| 美女 人体艺术 gogo|