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

    深海頂張式立管組合參激共振的非線性振動(dòng)分析

    2015-06-05 15:30:53唐友剛朱龍歡李楊青
    關(guān)鍵詞:渦激浮體立管

    唐友剛,朱龍歡,李楊青

    (1. 天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300072;2. 天津大學(xué)建筑工程學(xué)院,天津 300072)

    深海頂張式立管組合參激共振的非線性振動(dòng)分析

    唐友剛1,2,朱龍歡1,2,李楊青1,2

    (1. 天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300072;2. 天津大學(xué)建筑工程學(xué)院,天津 300072)

    研究了深海頂張式立管參數(shù)激勵(lì)和渦激共同作用下的非線性振動(dòng)特性. 考慮平臺(tái)升沉運(yùn)動(dòng)激勵(lì)和渦激力建立立管振動(dòng)方程,采用多尺度方法求解立管振動(dòng)方程的近似解析解. 考慮和型組合參激共振Ω≈ω1+ω2情況研究立管的振動(dòng)特性,計(jì)算得到了立管的幅頻響應(yīng)曲線,分析了平臺(tái)升沉運(yùn)動(dòng)對(duì)深海立管非線性振動(dòng)的影響. 結(jié)果表明:當(dāng)參激頻率滿足和型組合參激共振條件時(shí),立管振動(dòng)響應(yīng)中頻率為1/2參激頻率的亞諧波成分明顯;且由于內(nèi)共振關(guān)系的存在,立管1階模態(tài)被激發(fā),其幅值遠(yuǎn)大于2階模態(tài)幅值;隨著平臺(tái)升沉運(yùn)動(dòng)幅值的增大,立管橫向振動(dòng)幅值顯著增大,這表明平臺(tái)運(yùn)動(dòng)對(duì)于立管彎曲振動(dòng)有重要影響.

    頂張式立管(TTRs);組合參激共振;多尺度法;渦激振動(dòng)

    立管是海洋油氣開(kāi)發(fā)中連接海底井口和上部浮體的結(jié)構(gòu)[1].當(dāng)水深達(dá)到1,500,m以上時(shí),海洋環(huán)境尤為惡劣.立管頂部承受平臺(tái)升沉作用,內(nèi)部有高壓的油氣通過(guò),外部有復(fù)雜的浪、流等環(huán)境荷載作用,此環(huán)境下立管會(huì)發(fā)生碰撞、波頻振動(dòng)、渦激振動(dòng)、參激振動(dòng),這些都能引起立管的疲勞損傷和斷裂失效,不僅給工程本身帶來(lái)巨大經(jīng)濟(jì)損失,而且還會(huì)對(duì)自然環(huán)境造成嚴(yán)重的次生災(zāi)害[2].

    關(guān)于深海細(xì)長(zhǎng)結(jié)構(gòu)的參激振動(dòng)響應(yīng)問(wèn)題,國(guó)內(nèi)外學(xué)者進(jìn)行了大量理論和實(shí)驗(yàn)研究.Patel等[3-4]研究了低軸力條件下張力腿的不穩(wěn)定區(qū)域,用半解析的方法討論了張力腿在復(fù)合激勵(lì)下的動(dòng)力響應(yīng)問(wèn)題;Park等[5]用有限元方法研究了海洋細(xì)長(zhǎng)結(jié)構(gòu)物在參數(shù)激勵(lì)和波流激勵(lì)聯(lián)合作用下的動(dòng)力響應(yīng),并計(jì)算得到了馬蒂厄不穩(wěn)定參數(shù)區(qū)域;Chatjigeorgiou[6]針對(duì)參數(shù)激勵(lì)頻率接近立管非線性共振范圍的情況,研究參數(shù)激勵(lì)對(duì)立管振動(dòng)的影響,還討論了阻尼對(duì)于立管參數(shù)激勵(lì)振動(dòng)穩(wěn)定性的影響;Chatjigeorgiou等[7]在考慮內(nèi)流和參數(shù)激勵(lì)的情況下,針對(duì)參數(shù)激勵(lì)頻率接近兩倍1階固有頻率的特殊情況,得到立管振動(dòng)響應(yīng)的近似解析解.國(guó)內(nèi)學(xué)者朱克強(qiáng)[8]最早研究了考慮船舶的升沉和橫蕩運(yùn)動(dòng)情況時(shí)立管的振動(dòng)響應(yīng),用有限差分法求解了振動(dòng)方程,得到了立管在參數(shù)激勵(lì)和莫里森波浪力共同作用下的時(shí)間歷程響應(yīng);楊和振等[9]忽略波流力,考慮平臺(tái)浮體簡(jiǎn)諧升沉運(yùn)動(dòng)作為參數(shù)激勵(lì),研究深海立管參數(shù)激勵(lì)振動(dòng)響應(yīng);邵衛(wèi)東等[10]考慮立管頂端動(dòng)邊界條件及立管軸力隨水深的變化,給出了考慮浮體升沉及張緊環(huán)運(yùn)動(dòng)時(shí)頂張力管的固有特征.

    本文將立管上部浮體的升沉運(yùn)動(dòng)簡(jiǎn)化為參數(shù)激勵(lì),考慮立管大長(zhǎng)徑比引起的幾何非線性,綜合渦激力的影響,將立管振動(dòng)處理為參激-強(qiáng)迫激勵(lì)振動(dòng)系統(tǒng),建立深海立管振動(dòng)方程,采用多尺度方法求立管振動(dòng)的解析解,得到立管升沉運(yùn)動(dòng)激勵(lì)和渦激共同作用下非線性振動(dòng)的響應(yīng)特點(diǎn),并研究浮體升沉幅值對(duì)立管振動(dòng)的影響.

    1 立管分析模型

    立管為大長(zhǎng)細(xì)比結(jié)構(gòu),內(nèi)有不可壓縮無(wú)旋流體經(jīng)過(guò),管外承受的載荷包括順流向拖曳力、垂直于流速的渦激力以及上部浮體升沉運(yùn)動(dòng)誘發(fā)的軸向力.考慮底部約束為鉸支,立管坐標(biāo)系原點(diǎn)位于底端,軸豎直向上,建立參數(shù)激勵(lì)下深海立管模型,如圖1所示.

    圖1 參數(shù)激勵(lì)下深海立管模型示意Fig.1 Schematic diagram of deep-water riser under parametric excitation

    圖1 中,L為立管總長(zhǎng)度,v為均勻流速,上部浮體運(yùn)動(dòng)為簡(jiǎn)諧運(yùn)動(dòng),其位移為ua( t)=u0cos(Ωt),u0為位移幅值,上部浮體的運(yùn)動(dòng)會(huì)引起立管上端部的簡(jiǎn)諧升沉運(yùn)動(dòng).

    單位長(zhǎng)度流體和立管的受力情況如圖2和圖3所示.

    圖2 單位長(zhǎng)度流體受力Fig.2 Forces acting on the unit length of fluid

    圖3 單位長(zhǎng)度立管受力Fig.3 Forces acting on the unit length of riser

    考慮單位長(zhǎng)度流體和立管法向上受力的平衡條件[11-12],可得

    式中:m、ma和ml分別為單位長(zhǎng)度的立管質(zhì)量、附連水質(zhì)量和流體質(zhì)量;vl為立管管內(nèi)流體的速度;w為立管法向的位移;m0=(m+ma+m1) g ;S為剪切力;ε為局部應(yīng)變;χ為彎曲曲率;φ為立管切線方向與水平方向的夾角;EI、EA分別為立管彎曲剛度和軸向剛度;F為作用在立管單位長(zhǎng)度上的水動(dòng)力.F包括由漩渦泄放產(chǎn)生的渦激升力fL和立管在渦激力方向運(yùn)動(dòng)產(chǎn)生的拖曳力fD,即

    式中:wρ為周?chē)K拿芏?;D為立管外徑;sω為渦激頻率;LC、DC分別為升力系數(shù)和流體阻尼系數(shù).

    僅考慮立管在垂直來(lái)流平面的運(yùn)動(dòng),立管局部應(yīng)變?chǔ)欧譃殪o應(yīng)變?chǔ)?和動(dòng)應(yīng)變?chǔ)?兩部分,且任意點(diǎn)的切向位移u( z, t)為法向位移的高階小量O( w2),即

    此外,曲率、剪力、應(yīng)變及轉(zhuǎn)角之間有如下關(guān)系[7]:

    同時(shí),結(jié)構(gòu)所受的軸向張力eT分為靜頂張力e0T和動(dòng)張力e1T兩部分,其中靜頂張力e0T為定常值,且張力與應(yīng)變呈線性關(guān)系,即

    將式(2)~(11)代入式(1)中并進(jìn)行化簡(jiǎn),得到立管法向的運(yùn)動(dòng)微分方程

    應(yīng)用伽遼金法和振型疊加原理對(duì)式(12)進(jìn)行處理,設(shè)解的形式為

    式中:αu(z )、αw(z)滿足立管邊界條件,且僅與位置坐標(biāo)z有關(guān);gi( z)、fi( z)分別為立管軸向振動(dòng)和彎曲振動(dòng)的振型函數(shù);ua( t)、wa( t)分別為施加于立管頂端的軸向和橫向位移邊界條件,即ua( t)= u0cos(Ωt),wa( t)=0.將式(13)和(14)代入式(12),考慮結(jié)構(gòu)阻尼的作用,在等式左邊添加2μiq˙項(xiàng),其中μi為結(jié)構(gòu)阻尼系數(shù),結(jié)合振型的正交性進(jìn)行化簡(jiǎn),得到關(guān)于qi( t)的微分方程為

    式中iω為立管的第i階固有頻率.本文用解析方法求解式(15),僅考慮上部浮體的升沉運(yùn)動(dòng)和立管的橫向振動(dòng),因此可忽略式中含a()w t的乘積項(xiàng)和jkp q項(xiàng).由于阻尼項(xiàng)的存在,高階模態(tài)迅速收斂,因此僅保留前兩階模態(tài)進(jìn)行分析,由此整理式(15)得

    其中

    2 多尺度法求解模態(tài)振動(dòng)微分方程

    考慮整個(gè)振動(dòng)系統(tǒng)為弱非線性系統(tǒng),可將方程(16)和(17)中各非線性項(xiàng)系數(shù)看作小量,運(yùn)用多尺度法求解模態(tài)振動(dòng)方程.其中渦激升力項(xiàng)為周期變化的外激勵(lì),本文考慮渦激升力頻率遠(yuǎn)離立管固有振動(dòng)頻率,即外激勵(lì)為非共振的情形.

    式中e為無(wú)量綱小參數(shù).設(shè)方程(16)和(17)的解的形式為

    將式(27)代入式(16)和(17),并考慮式(26),令e同冪次項(xiàng)系數(shù)相等,并統(tǒng)一省略各變量中的“^”,得到式(28)~(31),即

    式中D0、D1分別是對(duì)T0、T1求偏導(dǎo)數(shù)的運(yùn)算符號(hào). 系數(shù)表達(dá)式分別為

    設(shè)式(28)和(29)解的形式為

    式中:cc代表前面表達(dá)式的共軛;參數(shù)iΛ的表達(dá)式為

    將式(35)和(36)代入式(30)和(31),消除永年項(xiàng),得到關(guān)于11()A T和22()A T的方程組.此振動(dòng)系統(tǒng)中存在幾種參激共振的情形,包括12Ωω≈,22Ωω≈,.這里Ω為參激頻率,1ω和 2ω分別為橫向彎曲振動(dòng)1階和2階固有頻率.本文考慮組合參激共振情形,對(duì)12+Ωωω≈的情況進(jìn)行詳細(xì)討論.同時(shí),該立管模型兩端為鉸支邊界條件,前兩階固有頻率之間存在2倍關(guān)系,即前兩階模態(tài)間存在內(nèi)共振關(guān)系,由此可以得到可解性條件為

    代入式(38)和(39)中,分離虛部與實(shí)部,為求振動(dòng)的穩(wěn)態(tài)響應(yīng),令10a′=,即可得到前兩階近似解的表達(dá)形式.

    3 算例及結(jié)果

    本文算例以深海中1,500,m長(zhǎng)的生產(chǎn)立管為例,計(jì)算和分析立管的振動(dòng)特性.對(duì)上文討論的解析近似解進(jìn)行編程,得到立管參數(shù)激勵(lì)下非線性振動(dòng)的幅頻響應(yīng)曲線,運(yùn)用數(shù)值解法求解響應(yīng)的時(shí)間歷程曲線,并做了對(duì)比分析.立管參數(shù)如表1所示.

    表1 深海1,500,m立管參數(shù)Tab.1 Parameters of 1,500,m deep-water riser

    將表1中參數(shù)代入式(38)和(39)中求解,得到組合參數(shù)激勵(lì)下立管振動(dòng)幅值隨解諧參數(shù)σ的變化情況,即幅頻響應(yīng)曲線,如圖4和圖5所示.

    組合參激共振情況下,由圖4和圖5分析可得,在σ趨近零的范圍內(nèi),兩階振幅均隨頻率增大而減小,且1階振幅1a明顯大于2階振幅2a,頻響曲線中的實(shí)線和虛線部分分別表示穩(wěn)定解和不穩(wěn)定解,這是由解的形式中兩個(gè)不同分支決定的.

    圖4 立管的1階模態(tài)幅頻響應(yīng)曲線Fig.4 Amplitude-frequency response curve of the 1st modal riser

    圖5 立管的2階模態(tài)幅頻響應(yīng)曲線Fig.5 Amplitude-frequency response curve of the 2nd modal riser

    對(duì)式(38)、(39)分離虛、實(shí)部的結(jié)果進(jìn)行直接數(shù)值求解,可以得到1階模態(tài)坐標(biāo)1q的時(shí)間歷程曲線,見(jiàn)圖6.

    圖6 1階模態(tài)坐標(biāo)1q的時(shí)間歷程曲線Fig.6 Time history curve of the 1st modal cordinate1q

    再對(duì)圖6所示振動(dòng)歷程進(jìn)行傅里葉變換,得到圖7所示振動(dòng)響應(yīng)譜.

    由圖7頻譜分析可得,響應(yīng)成分主要包括3個(gè):其一為參激頻率Ω的1/2的亞諧振動(dòng),即按照頻率0.296,9,rad/s進(jìn)行振動(dòng)的成分,此時(shí)立管產(chǎn)生明顯的亞諧參激振動(dòng);其二為渦激振動(dòng)部分,即按照渦激升力頻率1.609,0,rad/s進(jìn)行振動(dòng)的成分;其三為按照立管1階固有頻率0.199,2,rad/s進(jìn)行振動(dòng)的成分,產(chǎn)生該成分的主要原因是由于本文立管模型為兩端鉸支的梁結(jié)構(gòu),其橫向振動(dòng)的前兩階固有頻率之間存在內(nèi)共振關(guān)系,即212ωω≈,使2階固有頻率相對(duì)1階固有頻率而言恰好滿足參激共振條件,激發(fā)了立管1階模態(tài)的振動(dòng),由此說(shuō)明此時(shí)立管1階模態(tài)仍為主要激發(fā)模態(tài),且在參激共振和內(nèi)共振共同作用下,亞諧振動(dòng)諧波占主要部分,這表明了組合參數(shù)激勵(lì)對(duì)于立管振動(dòng)的重要影響..從圖4和圖5可得,0σ=時(shí),

    圖7 1階模態(tài)振動(dòng)響應(yīng)譜Fig.7 Spectrum of the 1st modal vibration

    改變頂部升沉幅值,得到1,maxa,如表2所示.由表2中數(shù)據(jù)可知,在參激頻率不變的情況下,立管1階模態(tài)幅值最大值隨浮體升沉幅值的增加而增加,當(dāng)浮體升沉幅值達(dá)到3,m時(shí),幅值響應(yīng)發(fā)生很明顯的增大.立管處于參數(shù)共振區(qū)時(shí)易發(fā)生失穩(wěn),需控制立管振動(dòng)響應(yīng).

    表2 1階模態(tài)幅值的最大值1,maxaTab.2 Maximum response value of the 1st modal1,maxa

    4 結(jié) 論

    (1) 理論分析結(jié)果與數(shù)值模擬分析結(jié)果基本吻合,說(shuō)明該理論分析結(jié)果是合理的,分析方法是可靠的.

    (2) 當(dāng)立管頂部參激頻率接近立管1階固有頻率與2階固有頻率之和時(shí),1階響應(yīng)明顯大于2階響應(yīng),為主要響應(yīng)成分.

    (3) 參數(shù)激勵(lì)使響應(yīng)中產(chǎn)生明顯的亞諧參激振動(dòng)成分,并且結(jié)構(gòu)本身的內(nèi)共振關(guān)系使2階固有頻率相對(duì)1階固有頻率滿足參激共振條件,在立管振動(dòng)系統(tǒng)中產(chǎn)生重要作用.

    (4) 在參激頻率不變的情況下,立管1階模態(tài)幅值最大值隨浮體升沉幅值的增加而增加,并且當(dāng)升沉幅值較大時(shí),立管橫向振動(dòng)幅值顯著增大.因此在實(shí)際設(shè)計(jì)中,應(yīng)該避免立管頂部發(fā)生較大的升沉幅值.

    (5) 本文選取渦泄頻率遠(yuǎn)離立管橫向振動(dòng)固有頻率,而當(dāng)外激頻率與某階固有頻率接近時(shí),渦激共振與參激共振同時(shí)發(fā)生,可能導(dǎo)致立管振動(dòng)大幅加劇,使立管處于危險(xiǎn)狀態(tài),因此在實(shí)際工程中,還需要關(guān)注兩種諧振同時(shí)發(fā)生的情形,確保立管安全工作.

    [1] 王一飛. 深海立管渦激振動(dòng)疲勞損傷預(yù)報(bào)方法研究[D]. 上海:上海交通大學(xué)船舶海洋與建筑工程學(xué)院,2008. Wang Yifei. Prediction Methods for the VIV-induced Fatigue Damage in Deep Sea Risers [D]. Shanghai:School of Naval Architecture,Ocean and Civil Engineering,Shanghai Jiao Tong University,2008(in Chinese).

    [2] 陳云水,王德禹. 深水立管的若干結(jié)構(gòu)力學(xué)研究進(jìn)展[J]. 中國(guó)海洋平臺(tái),2007,22(5):1-5. Chen Yunshui,Wang Deyu. Some advance of research on structural mechanics of deep water risers[J]. China Offshore Platform,2007,22(5):1-5(in Chinese).

    [3] Patel M H,Park H I. Dynamics of tension leg platform tethers at low tension (Part I):Mathieu stability at large parameters[J]. Marine Structures,1991,4(3):257-273.

    [4] Patel M H,Park H I. Combined axial and lateral responses of tensioned buoyant platform tethers[J]. Engineering Structures,1995,17(10):687-695.

    [5] Park H I,Jung D H. A finite element method for dynamic analysis of long slender marine structures under combined parametric and forcing excitations[J]. Ocean Engineering,2002,29(11):1313-1325.

    [6] Chatjigeorgiou I K. On the parametric excitation of vertical elastic slender structures and the effect of damping in marine applications[J]. Applied Ocean Research,2004,26(1/2):23-33.

    [7] Chatjigeorgiou I K,Mavrakos S A. Nonlinear resonances of parametrically excited risers numerical and analytical investigation for12Ωω?[J]. Computers and Structures,2005,83(8/9):560-573.

    [8] 朱克強(qiáng). 駁船升沉和橫蕩運(yùn)動(dòng)下海洋立管動(dòng)力響應(yīng)[J]. 中外船舶科技,2002(3):1-4. Zhu Keqiang. Dynamic response of submarine riser connected to a vessel with heave and lateral motion[J]. Shipbuilding Science and Technology,2002(3):1-4(in Chinese).

    [9] 楊和振,李華軍. 參數(shù)激勵(lì)下深海立管動(dòng)力特性研究[J]. 振動(dòng)與沖擊,2009,28(9):65-69. Yang Hezhen,Li Huajun. Vibration analysis of deep-sea risers under parametric excitations[J]. Journal of Vibration and Shock,2009,28(9):65-69(in Chinese).

    [10] 邵衛(wèi)東,唐友剛,樊娟娟,等. 考慮浮體升沉及張緊環(huán)運(yùn)動(dòng)深海立管固有振動(dòng)特征研究[J]. 海洋工程,2012,30(2):8-13. Shao Weidong,Tang Yougang,F(xiàn)an Juanjuan,et al. Study of natural vibration characteristics of deep-water riser considering heave motion of platform and tensionring’s motion[J]. The Ocean Engineering,2012,30(2):8-13(in Chinese).

    [11] Paidoussis M P,Li G X. Pipes conveying fluid:A model dynamical problem[J]. Journal of Fluids and Structures,1993,7(2):137-204.

    [12] 李楊青. 深海頂張力立管非線性振動(dòng)特性研究[D]. 天津:天津大學(xué)建筑工程學(xué)院,2013. Li Yangqing. Study on the Nonlinear Vibration Behaviors of Top-Tensioned Risers in Deep Water[D]. Tianjin:School of Civil Engineering,Tianjin University,2013(in Chinese).

    (責(zé)任編輯:樊素英)

    Nonlinear Vibration Analysis of Combination Parametric Resonance for TTRs in Deep Water

    Tang Yougang1,2,Zhu Longhuan1,2,Li Yangqing1,2
    (1. State Key Laboratory of Hydraulic Engineering Simulation and Safety,Tianjin University,Tianjin 300072,China;2. School of Civil Engineering,Tianjin University,Tianjin 300072,China)

    The nonlinear vibration characteristics induced by parameters and vortex were investigated for toptensioned risers (TTRs) in deep water. The vibration equations were established with consideration of both the parametric excitation caused by the heave motion of platform and vortex induced force. Approximated analytical solution to the equations was calculated using the method of multiple scales. Considering the case of combination parametric resonance,Ω≈ω1+ω2, the frequency response curves of the riser were obtained and the influence of the heave motion of platform on the nonlinear vibration of the riser was analyzed. The results show that when the parametric excitation frequency satisfies the combination parametric resonance conditions, the sub-harmonic resonances (response frequency equal to 1/2 parametric frequency) is obvious. Because of the existence of internal resonance, the vibration response of the first mode is excited and the amplitude is much larger than that of the second mode. As the heave motion amplitude increases, the amplitude of the lateral vibration response of TTRs increases obviously, which indicates that the heave motion of platform has important influence on the lateral vibration of TTRs。

    top-tensioned risers(TTRs);combination parametric resonance;method of multiple scales;vortex induced vibration

    P756.2

    A

    0493-2137(2015)09-0811-06

    10.11784/tdxbz201402022

    2014-02-18;

    2014-09-02.

    國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)資助項(xiàng)目(2014CB0468);國(guó)家自然科學(xué)基金重點(diǎn)資助項(xiàng)目(51239008).

    唐友剛(1952— ),博士,教授.

    唐友剛,tangyougang_td@163.com.

    猜你喜歡
    渦激浮體立管
    人民長(zhǎng)江(2023年6期)2023-07-25 12:24:14
    不同間距比下串聯(lián)圓柱渦激振動(dòng)數(shù)值模擬研究
    物探船硬浮體陣列自擴(kuò)變量分析與應(yīng)用
    超大型浮體結(jié)構(gòu)碰撞損傷研究
    渦激振動(dòng)發(fā)電裝置及其關(guān)鍵技術(shù)
    常見(jiàn)高層建筑物室內(nèi)給水立管材質(zhì)解析
    有限流動(dòng)水域浮體受力及側(cè)傾研究
    盤(pán)球立管結(jié)構(gòu)抑制渦激振動(dòng)的數(shù)值分析方法研究
    電子制作(2018年14期)2018-08-21 01:38:42
    深水鋼懸鏈立管J型鋪設(shè)研究
    柔性圓管在渦激振動(dòng)下的模態(tài)響應(yīng)分析
    国产精品偷伦视频观看了| 悠悠久久av| 日韩制服丝袜自拍偷拍| 国产不卡av网站在线观看| videosex国产| 国产成人精品无人区| 欧美精品一区二区大全| 免费在线观看黄色视频的| 亚洲精品国产一区二区精华液| 少妇的丰满在线观看| 大片电影免费在线观看免费| 一区二区三区四区激情视频| 狂野欧美激情性bbbbbb| 精品亚洲成国产av| 亚洲精品久久成人aⅴ小说| videosex国产| 亚洲国产av影院在线观看| 伦理电影免费视频| 久久精品亚洲熟妇少妇任你| 亚洲精品中文字幕在线视频| 国产伦理片在线播放av一区| 亚洲国产欧美一区二区综合| 一区二区av电影网| 2018国产大陆天天弄谢| netflix在线观看网站| 国产精品99久久99久久久不卡| 97精品久久久久久久久久精品| 亚洲专区国产一区二区| 日本五十路高清| 久久毛片免费看一区二区三区| av视频免费观看在线观看| videosex国产| 最新的欧美精品一区二区| 操美女的视频在线观看| 久久99热这里只频精品6学生| 天堂8中文在线网| 丝袜美足系列| 黄网站色视频无遮挡免费观看| 国产精品久久久久久精品电影小说| 国产精品久久久久久精品电影小说| 久久久久国产精品人妻一区二区| 欧美成人午夜精品| 啦啦啦视频在线资源免费观看| 久久精品久久久久久噜噜老黄| 国产片特级美女逼逼视频| 搡老乐熟女国产| 天天添夜夜摸| 一区福利在线观看| 国产成人系列免费观看| 久久久久久久国产电影| 亚洲色图综合在线观看| 亚洲伊人久久精品综合| 99热网站在线观看| 中文精品一卡2卡3卡4更新| 少妇的丰满在线观看| 国产成人免费观看mmmm| 可以免费在线观看a视频的电影网站| 嫁个100分男人电影在线观看 | a级毛片在线看网站| 亚洲av日韩在线播放| 免费女性裸体啪啪无遮挡网站| 国产精品国产av在线观看| 悠悠久久av| 一本久久精品| 在线 av 中文字幕| 女性被躁到高潮视频| 亚洲专区国产一区二区| 国产精品国产三级国产专区5o| 91精品伊人久久大香线蕉| 日日夜夜操网爽| 少妇猛男粗大的猛烈进出视频| 高清av免费在线| 欧美少妇被猛烈插入视频| 狂野欧美激情性xxxx| 男女边摸边吃奶| 制服诱惑二区| 操出白浆在线播放| 免费在线观看黄色视频的| 国产精品一区二区在线不卡| 免费在线观看视频国产中文字幕亚洲 | 精品少妇内射三级| 手机成人av网站| 秋霞在线观看毛片| 精品亚洲成a人片在线观看| 国产欧美日韩一区二区三 | 晚上一个人看的免费电影| 久久久久久久大尺度免费视频| 女人高潮潮喷娇喘18禁视频| 亚洲国产日韩一区二区| 超色免费av| av片东京热男人的天堂| 美女扒开内裤让男人捅视频| netflix在线观看网站| 亚洲熟女毛片儿| 午夜免费成人在线视频| 性高湖久久久久久久久免费观看| 天天躁狠狠躁夜夜躁狠狠躁| a级片在线免费高清观看视频| 日本黄色日本黄色录像| 国产黄色免费在线视频| 精品第一国产精品| 九色亚洲精品在线播放| 无遮挡黄片免费观看| 久久女婷五月综合色啪小说| 日韩伦理黄色片| 国产成人啪精品午夜网站| 9191精品国产免费久久| 国产日韩欧美亚洲二区| 成人免费观看视频高清| 日本午夜av视频| 一级毛片电影观看| 欧美国产精品va在线观看不卡| 不卡av一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 成人午夜精彩视频在线观看| 少妇裸体淫交视频免费看高清 | 不卡av一区二区三区| 黄网站色视频无遮挡免费观看| 精品少妇一区二区三区视频日本电影| 99久久人妻综合| kizo精华| 视频在线观看一区二区三区| 可以免费在线观看a视频的电影网站| 精品欧美一区二区三区在线| 真人做人爱边吃奶动态| 欧美亚洲日本最大视频资源| 1024香蕉在线观看| 欧美日韩福利视频一区二区| 制服人妻中文乱码| 久久人人爽人人片av| 亚洲成人手机| 欧美另类一区| 美女大奶头黄色视频| 欧美人与性动交α欧美精品济南到| a级片在线免费高清观看视频| 久久精品熟女亚洲av麻豆精品| 热99久久久久精品小说推荐| 久久精品国产亚洲av高清一级| 一级,二级,三级黄色视频| 色婷婷av一区二区三区视频| 欧美中文综合在线视频| 操美女的视频在线观看| 色综合欧美亚洲国产小说| 欧美国产精品一级二级三级| 亚洲,一卡二卡三卡| 久久精品久久精品一区二区三区| 日日爽夜夜爽网站| 欧美日韩av久久| 91成人精品电影| 国产av一区二区精品久久| 亚洲国产中文字幕在线视频| 久久久久国产一级毛片高清牌| 婷婷丁香在线五月| 久久人人爽av亚洲精品天堂| 亚洲国产欧美网| 一边摸一边做爽爽视频免费| 亚洲精品国产av蜜桃| 午夜免费观看性视频| 国产不卡av网站在线观看| 亚洲国产精品国产精品| 男女床上黄色一级片免费看| 人体艺术视频欧美日本| 中文字幕另类日韩欧美亚洲嫩草| 午夜福利,免费看| 久久久久精品人妻al黑| 国产亚洲欧美在线一区二区| 老司机影院毛片| 天天添夜夜摸| 国产成人影院久久av| 欧美日韩成人在线一区二区| 香蕉国产在线看| 可以免费在线观看a视频的电影网站| 天天躁夜夜躁狠狠久久av| 欧美日韩一级在线毛片| 操美女的视频在线观看| 午夜久久久在线观看| 久久久久精品国产欧美久久久 | 午夜免费成人在线视频| 亚洲精品一二三| 久久人妻福利社区极品人妻图片 | 日日爽夜夜爽网站| 亚洲五月色婷婷综合| 嫩草影视91久久| 国产精品二区激情视频| 国产成人精品在线电影| 欧美日韩成人在线一区二区| 色综合欧美亚洲国产小说| 国产高清国产精品国产三级| 午夜免费鲁丝| 亚洲国产最新在线播放| 欧美人与性动交α欧美软件| 日韩免费高清中文字幕av| 久久人人爽人人片av| 中文精品一卡2卡3卡4更新| 国产伦人伦偷精品视频| 日韩,欧美,国产一区二区三区| 久久精品aⅴ一区二区三区四区| 日韩制服丝袜自拍偷拍| 在线 av 中文字幕| 777久久人妻少妇嫩草av网站| 国产精品久久久久成人av| av不卡在线播放| 每晚都被弄得嗷嗷叫到高潮| 久久精品aⅴ一区二区三区四区| 天天影视国产精品| 亚洲精品美女久久av网站| 一级毛片黄色毛片免费观看视频| 色94色欧美一区二区| 久久鲁丝午夜福利片| 国产主播在线观看一区二区 | 国产成人精品无人区| 国产精品国产三级专区第一集| 亚洲黑人精品在线| 中文字幕人妻熟女乱码| 18禁裸乳无遮挡动漫免费视频| 欧美黄色淫秽网站| 一级毛片电影观看| 人人妻人人添人人爽欧美一区卜| 两个人看的免费小视频| 日韩 亚洲 欧美在线| 久热爱精品视频在线9| 日韩制服丝袜自拍偷拍| 91字幕亚洲| 999久久久国产精品视频| 午夜福利乱码中文字幕| 交换朋友夫妻互换小说| 最新的欧美精品一区二区| 亚洲欧洲国产日韩| av又黄又爽大尺度在线免费看| 午夜激情久久久久久久| 各种免费的搞黄视频| 欧美黄色淫秽网站| av天堂久久9| 国产免费一区二区三区四区乱码| 国产精品久久久久成人av| 黄色 视频免费看| 成年动漫av网址| 大香蕉久久网| 亚洲av国产av综合av卡| 在线av久久热| 欧美日韩国产mv在线观看视频| 我要看黄色一级片免费的| 男人舔女人的私密视频| 国产欧美亚洲国产| 操美女的视频在线观看| 一级黄片播放器| 亚洲男人天堂网一区| 欧美在线一区亚洲| 午夜福利影视在线免费观看| avwww免费| 精品国产乱码久久久久久小说| 爱豆传媒免费全集在线观看| 成年人午夜在线观看视频| 99精品久久久久人妻精品| 50天的宝宝边吃奶边哭怎么回事| 成年av动漫网址| 欧美激情高清一区二区三区| 亚洲欧美日韩高清在线视频 | 成人国语在线视频| 9191精品国产免费久久| 国产无遮挡羞羞视频在线观看| 岛国毛片在线播放| 国产精品人妻久久久影院| 美女国产高潮福利片在线看| 国产精品久久久av美女十八| 欧美亚洲 丝袜 人妻 在线| 一本大道久久a久久精品| 久久久欧美国产精品| 欧美亚洲 丝袜 人妻 在线| av在线播放精品| 成年av动漫网址| 一级毛片黄色毛片免费观看视频| 国产精品欧美亚洲77777| 久久精品久久久久久久性| 亚洲av日韩精品久久久久久密 | 校园人妻丝袜中文字幕| 乱人伦中国视频| 汤姆久久久久久久影院中文字幕| 国产高清国产精品国产三级| 韩国高清视频一区二区三区| 男女下面插进去视频免费观看| 欧美黄色淫秽网站| 日韩一区二区三区影片| 少妇裸体淫交视频免费看高清 | 尾随美女入室| 大香蕉久久成人网| 亚洲成人手机| 久久久久精品人妻al黑| www.精华液| 久久久久国产精品人妻一区二区| 亚洲成人国产一区在线观看 | 亚洲国产欧美日韩在线播放| 99精国产麻豆久久婷婷| 亚洲情色 制服丝袜| 欧美精品亚洲一区二区| 国产亚洲欧美在线一区二区| 国产男人的电影天堂91| 在线亚洲精品国产二区图片欧美| 久久中文字幕一级| 亚洲精品久久午夜乱码| 男人添女人高潮全过程视频| e午夜精品久久久久久久| 国产日韩欧美视频二区| 久久天堂一区二区三区四区| 成人亚洲欧美一区二区av| 成在线人永久免费视频| 亚洲,欧美,日韩| 叶爱在线成人免费视频播放| 只有这里有精品99| 亚洲天堂av无毛| 看十八女毛片水多多多| 美女大奶头黄色视频| 老司机午夜十八禁免费视频| 免费一级毛片在线播放高清视频 | 真人做人爱边吃奶动态| 国产精品久久久av美女十八| 日韩av在线免费看完整版不卡| 成年女人毛片免费观看观看9 | 亚洲人成电影观看| 丰满人妻熟妇乱又伦精品不卡| 97人妻天天添夜夜摸| 中文字幕人妻丝袜制服| 国产成人免费无遮挡视频| 久久亚洲国产成人精品v| 18禁裸乳无遮挡动漫免费视频| 9热在线视频观看99| 国产亚洲精品久久久久5区| √禁漫天堂资源中文www| 亚洲伊人色综图| 只有这里有精品99| 精品福利观看| 亚洲精品自拍成人| 黄色毛片三级朝国网站| 国产91精品成人一区二区三区 | 天天躁夜夜躁狠狠躁躁| 深夜精品福利| av网站在线播放免费| 国产片特级美女逼逼视频| 国产在线免费精品| 99久久精品国产亚洲精品| 中文字幕亚洲精品专区| av福利片在线| 国产精品免费视频内射| 国产福利在线免费观看视频| 亚洲中文日韩欧美视频| 亚洲第一青青草原| 一级毛片电影观看| 18禁国产床啪视频网站| 美女扒开内裤让男人捅视频| 亚洲精品在线美女| 天堂中文最新版在线下载| 日日夜夜操网爽| 在线精品无人区一区二区三| 国产视频首页在线观看| 久久毛片免费看一区二区三区| 一区在线观看完整版| 中文字幕高清在线视频| 超碰97精品在线观看| 国产精品二区激情视频| 亚洲色图综合在线观看| www.999成人在线观看| 色婷婷久久久亚洲欧美| 51午夜福利影视在线观看| 精品少妇内射三级| 国产精品九九99| 亚洲国产精品一区三区| 国产不卡av网站在线观看| 狂野欧美激情性xxxx| 丝袜喷水一区| 99热网站在线观看| 少妇粗大呻吟视频| 国产精品成人在线| 日韩一区二区三区影片| 香蕉国产在线看| xxx大片免费视频| 国产精品国产三级国产专区5o| 国产免费福利视频在线观看| 国产精品一区二区免费欧美 | 美女中出高潮动态图| 精品国产乱码久久久久久小说| 在线观看免费午夜福利视频| 啦啦啦在线观看免费高清www| 麻豆乱淫一区二区| 狠狠婷婷综合久久久久久88av| 十分钟在线观看高清视频www| 亚洲成人免费电影在线观看 | 黄片播放在线免费| 亚洲成人免费av在线播放| 老司机午夜十八禁免费视频| 亚洲 国产 在线| 欧美人与性动交α欧美精品济南到| 最近手机中文字幕大全| 人人妻人人爽人人添夜夜欢视频| 2018国产大陆天天弄谢| 国产亚洲av片在线观看秒播厂| av福利片在线| 在线亚洲精品国产二区图片欧美| 最黄视频免费看| 亚洲中文字幕日韩| 国产成人精品久久久久久| 欧美乱码精品一区二区三区| 久久久国产精品麻豆| 国产成人免费无遮挡视频| 久热这里只有精品99| 一本—道久久a久久精品蜜桃钙片| 99国产精品99久久久久| av在线app专区| 久久影院123| 又紧又爽又黄一区二区| 国产av一区二区精品久久| 亚洲人成电影免费在线| 免费在线观看视频国产中文字幕亚洲 | 香蕉国产在线看| 欧美乱码精品一区二区三区| 在线观看免费视频网站a站| 国产一区二区三区av在线| 色婷婷久久久亚洲欧美| 男女床上黄色一级片免费看| 欧美精品av麻豆av| 国产伦理片在线播放av一区| 中文字幕最新亚洲高清| 精品人妻一区二区三区麻豆| 大话2 男鬼变身卡| 丰满人妻熟妇乱又伦精品不卡| 国产黄色视频一区二区在线观看| 国产精品久久久人人做人人爽| 一区二区三区四区激情视频| 亚洲精品第二区| 女性生殖器流出的白浆| 美女扒开内裤让男人捅视频| 国语对白做爰xxxⅹ性视频网站| 97精品久久久久久久久久精品| 一区福利在线观看| 亚洲av成人精品一二三区| 亚洲三区欧美一区| 黄色一级大片看看| xxx大片免费视频| 热re99久久精品国产66热6| 男女之事视频高清在线观看 | 日韩 亚洲 欧美在线| 亚洲人成电影观看| 麻豆乱淫一区二区| 婷婷色综合大香蕉| 国产高清视频在线播放一区 | 99国产综合亚洲精品| 色视频在线一区二区三区| 丝袜喷水一区| 成年人黄色毛片网站| 国产爽快片一区二区三区| 极品人妻少妇av视频| av天堂在线播放| 国产精品秋霞免费鲁丝片| 国产成人系列免费观看| 午夜激情久久久久久久| 欧美老熟妇乱子伦牲交| 两个人看的免费小视频| 亚洲人成77777在线视频| 91精品国产国语对白视频| 精品国产一区二区三区久久久樱花| av国产精品久久久久影院| 老汉色av国产亚洲站长工具| 亚洲国产最新在线播放| 精品国产乱码久久久久久男人| 又黄又粗又硬又大视频| 大话2 男鬼变身卡| 一级黄片播放器| 1024视频免费在线观看| 欧美成人午夜精品| 波多野结衣一区麻豆| 亚洲少妇的诱惑av| 少妇的丰满在线观看| 国产又色又爽无遮挡免| 国产视频首页在线观看| 久久久欧美国产精品| 亚洲色图综合在线观看| 免费观看人在逋| 欧美人与性动交α欧美软件| 中文字幕人妻丝袜制服| 日韩一区二区三区影片| 日本av免费视频播放| 老熟女久久久| www.精华液| 欧美激情 高清一区二区三区| 久久久国产一区二区| 午夜免费男女啪啪视频观看| 久久国产精品人妻蜜桃| 午夜福利在线免费观看网站| 又黄又粗又硬又大视频| 久久人妻熟女aⅴ| 久久国产精品影院| 日本色播在线视频| 狠狠精品人妻久久久久久综合| 在线观看免费日韩欧美大片| 十八禁高潮呻吟视频| 色婷婷av一区二区三区视频| 国产午夜精品一二区理论片| 日日夜夜操网爽| 久久精品人人爽人人爽视色| 无限看片的www在线观看| 亚洲国产精品一区二区三区在线| av国产久精品久网站免费入址| 99国产精品一区二区蜜桃av | 美女国产高潮福利片在线看| 叶爱在线成人免费视频播放| 日韩一区二区三区影片| 国产精品一国产av| 91字幕亚洲| 国产一区二区在线观看av| 亚洲国产毛片av蜜桃av| 成年女人毛片免费观看观看9 | 性色av一级| 又大又爽又粗| 国产在线一区二区三区精| 18禁裸乳无遮挡动漫免费视频| 男人添女人高潮全过程视频| √禁漫天堂资源中文www| 精品人妻熟女毛片av久久网站| 国产高清videossex| 美女午夜性视频免费| 欧美大码av| 欧美人与善性xxx| 热re99久久国产66热| 叶爱在线成人免费视频播放| 老司机在亚洲福利影院| 在线观看免费高清a一片| 两个人免费观看高清视频| 在线观看国产h片| 欧美激情极品国产一区二区三区| 精品一区在线观看国产| av视频免费观看在线观看| av在线播放精品| 国产免费现黄频在线看| 少妇 在线观看| 91国产中文字幕| 波多野结衣av一区二区av| 欧美日韩亚洲国产一区二区在线观看 | 9色porny在线观看| 赤兔流量卡办理| 久久综合国产亚洲精品| 国产一区二区三区综合在线观看| 国产成人av激情在线播放| 国产在线免费精品| 曰老女人黄片| 在线观看免费日韩欧美大片| 悠悠久久av| av又黄又爽大尺度在线免费看| 久热这里只有精品99| 青青草视频在线视频观看| 亚洲九九香蕉| 狂野欧美激情性bbbbbb| 国产一区二区在线观看av| 国产xxxxx性猛交| 免费在线观看日本一区| 国产精品av久久久久免费| 国产视频一区二区在线看| 久9热在线精品视频| 男女无遮挡免费网站观看| 亚洲国产精品一区三区| 99精国产麻豆久久婷婷| 欧美精品亚洲一区二区| 一区二区av电影网| 精品国产一区二区三区四区第35| 一级毛片电影观看| 老汉色av国产亚洲站长工具| 下体分泌物呈黄色| 如日韩欧美国产精品一区二区三区| 亚洲精品一区蜜桃| 首页视频小说图片口味搜索 | 日韩制服骚丝袜av| 两个人免费观看高清视频| 在线观看国产h片| 黄频高清免费视频| 亚洲精品成人av观看孕妇| 七月丁香在线播放| 久久精品aⅴ一区二区三区四区| 青青草视频在线视频观看| 国产一区二区激情短视频 | 中文字幕高清在线视频| 一二三四社区在线视频社区8| 精品久久蜜臀av无| 97精品久久久久久久久久精品| 国产一区二区激情短视频 | 欧美成狂野欧美在线观看| 午夜视频精品福利| 欧美亚洲 丝袜 人妻 在线| a 毛片基地| 黄频高清免费视频| 无限看片的www在线观看| 国产成人精品无人区| 亚洲欧美清纯卡通| 国产精品二区激情视频| 中文字幕高清在线视频| 岛国毛片在线播放| 丁香六月欧美| 在线观看免费视频网站a站| 国产高清国产精品国产三级| 美女中出高潮动态图| 大片电影免费在线观看免费| 国产不卡av网站在线观看| www.熟女人妻精品国产| www.999成人在线观看| 黄色视频在线播放观看不卡| 观看av在线不卡| 赤兔流量卡办理| 婷婷色综合www| 嫁个100分男人电影在线观看 | 欧美人与性动交α欧美精品济南到| 国产日韩欧美亚洲二区| 国产视频首页在线观看| 中文字幕人妻丝袜制服| 亚洲精品在线美女| 国产男女内射视频| 久久精品国产亚洲av涩爱| 水蜜桃什么品种好| 亚洲人成电影免费在线|