• <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视频在线免费观看| 夜夜爽天天搞| 天天躁日日操中文字幕| 91在线精品国自产拍蜜月| 最后的刺客免费高清国语| xxxwww97欧美| 国产成人av教育| 91麻豆av在线| 精品久久久久久久人妻蜜臀av| 看免费成人av毛片| 日韩国内少妇激情av| 国产毛片a区久久久久| 亚洲中文字幕日韩| 日韩人妻高清精品专区| 国产av在哪里看| 成人三级黄色视频| 国产私拍福利视频在线观看| 免费高清视频大片| 中文字幕熟女人妻在线| 欧美性感艳星| 毛片一级片免费看久久久久 | 在线免费十八禁| xxxwww97欧美| 极品教师在线视频| 免费看av在线观看网站| 黄色丝袜av网址大全| 又黄又爽又免费观看的视频| 国产免费av片在线观看野外av| 国产高清有码在线观看视频| 日韩欧美在线二视频| 欧美日韩瑟瑟在线播放| 午夜a级毛片| 两个人视频免费观看高清| 男人狂女人下面高潮的视频| 看十八女毛片水多多多| 久久精品国产清高在天天线| 久久久色成人| 波多野结衣巨乳人妻| 22中文网久久字幕| 噜噜噜噜噜久久久久久91| 一进一出抽搐动态| 欧美日本视频| www日本黄色视频网| 一个人免费在线观看电影| 亚洲成a人片在线一区二区| 国产精品人妻久久久久久| 又粗又爽又猛毛片免费看| 亚洲精品粉嫩美女一区| aaaaa片日本免费| 午夜日韩欧美国产| 久久国产乱子免费精品| 毛片一级片免费看久久久久 | 亚洲专区中文字幕在线| 日韩欧美免费精品| 欧美高清性xxxxhd video| 国产精品1区2区在线观看.| 国产 一区精品| 在线观看美女被高潮喷水网站| 两人在一起打扑克的视频| 亚洲国产日韩欧美精品在线观看| 又粗又爽又猛毛片免费看| 综合色av麻豆| 久久精品国产亚洲网站| 中文字幕高清在线视频| 一个人看视频在线观看www免费| 在线观看美女被高潮喷水网站| or卡值多少钱| 天天躁日日操中文字幕| 日本一二三区视频观看| 中亚洲国语对白在线视频| 禁无遮挡网站| 91精品国产九色| 国产色爽女视频免费观看| 国产午夜福利久久久久久| 久久久久久久久久黄片| 三级毛片av免费| 人人妻人人看人人澡| 日韩精品青青久久久久久| 成年女人永久免费观看视频| 女人被狂操c到高潮| 欧美绝顶高潮抽搐喷水| 99热网站在线观看| 三级毛片av免费| 国产精品不卡视频一区二区| 久久精品夜夜夜夜夜久久蜜豆| 久久久久久久精品吃奶| 久久人妻av系列| 成人鲁丝片一二三区免费| 国产午夜精品论理片| 国产午夜精品久久久久久一区二区三区 | 日韩欧美在线二视频| 人妻制服诱惑在线中文字幕| 久久久久久久久久黄片| 欧美色欧美亚洲另类二区| 尾随美女入室| 黄色女人牲交| 午夜爱爱视频在线播放| 免费av观看视频| 日韩欧美国产一区二区入口| 一夜夜www| 99热这里只有精品一区| 999久久久精品免费观看国产| 在线免费观看的www视频| 在现免费观看毛片| 日日摸夜夜添夜夜添av毛片 | 日韩欧美在线二视频| 九九在线视频观看精品| 女生性感内裤真人,穿戴方法视频| 国产麻豆成人av免费视频| 国产精品不卡视频一区二区| 色综合亚洲欧美另类图片| 亚洲av免费高清在线观看| 日韩一区二区视频免费看| 精品人妻1区二区| 婷婷色综合大香蕉| 亚洲av一区综合| 极品教师在线视频| 国产精品一区二区性色av| 久久久久久国产a免费观看| aaaaa片日本免费| 一个人观看的视频www高清免费观看| 国产精品三级大全| 蜜桃亚洲精品一区二区三区| 国产精品亚洲一级av第二区| 国产精品爽爽va在线观看网站| 亚洲成a人片在线一区二区| 亚洲国产精品合色在线| 久久久久久久久久黄片| 国产视频内射| 亚洲中文字幕日韩| 一级毛片久久久久久久久女| 乱码一卡2卡4卡精品| 人妻夜夜爽99麻豆av| 99热网站在线观看| 亚洲aⅴ乱码一区二区在线播放| 免费在线观看成人毛片| 亚洲av二区三区四区| 亚洲av.av天堂| 国产欧美日韩一区二区精品| 久久久久国产精品人妻aⅴ院| 久久久久久伊人网av| 国产真实乱freesex| 尾随美女入室| 亚州av有码| 热99re8久久精品国产| 色5月婷婷丁香| 九九久久精品国产亚洲av麻豆| 给我免费播放毛片高清在线观看| 韩国av一区二区三区四区| 午夜福利欧美成人| 午夜精品久久久久久毛片777| 国产亚洲精品av在线| eeuss影院久久| 男人的好看免费观看在线视频| 免费观看精品视频网站| avwww免费| 成年人黄色毛片网站| 午夜免费男女啪啪视频观看 | 在线观看美女被高潮喷水网站| 日韩欧美在线二视频| 亚洲一区二区三区色噜噜| 最近视频中文字幕2019在线8| 国产精品久久久久久久久免| 伊人久久精品亚洲午夜| av天堂在线播放| 亚洲av五月六月丁香网| 波多野结衣高清无吗| 美女xxoo啪啪120秒动态图| 乱人视频在线观看| 夜夜夜夜夜久久久久| 亚洲精华国产精华液的使用体验 | 免费看美女性在线毛片视频| 久久精品综合一区二区三区| 色哟哟·www| 免费人成视频x8x8入口观看| 99热这里只有精品一区| 婷婷亚洲欧美| 无遮挡黄片免费观看| 观看美女的网站| 国产视频内射| 欧美日韩瑟瑟在线播放| 久久久久久久久大av| 国产精品久久视频播放| 欧美一级a爱片免费观看看| 精品一区二区三区人妻视频| 国产一区二区亚洲精品在线观看| 精品一区二区三区视频在线| av专区在线播放| 亚洲专区中文字幕在线| 我要看日韩黄色一级片| 亚洲精品久久国产高清桃花| 黄色欧美视频在线观看| 亚洲人成网站在线播放欧美日韩| 在线免费观看的www视频| 国产一区二区三区在线臀色熟女| 51国产日韩欧美| 女的被弄到高潮叫床怎么办 | 小说图片视频综合网站| 搡老岳熟女国产| 99国产精品一区二区蜜桃av| bbb黄色大片| 亚洲人成网站高清观看| 麻豆久久精品国产亚洲av| 亚洲精品色激情综合| 给我免费播放毛片高清在线观看| 淫妇啪啪啪对白视频| 免费在线观看影片大全网站| 变态另类丝袜制服| 日韩亚洲欧美综合| 久久精品国产亚洲av天美| 国产精品乱码一区二三区的特点| 成人毛片a级毛片在线播放| 亚洲内射少妇av| 国产精品一区www在线观看 | 久久婷婷人人爽人人干人人爱| 国产黄色小视频在线观看| 欧美3d第一页| 在线观看舔阴道视频| 亚洲专区中文字幕在线| 深夜精品福利| 亚洲av免费高清在线观看| 国产av麻豆久久久久久久| 久久人人精品亚洲av| 国产成人影院久久av| 色5月婷婷丁香| 深夜精品福利| 国产在视频线在精品| 午夜免费成人在线视频| 真人一进一出gif抽搐免费| 又黄又爽又免费观看的视频| 美女xxoo啪啪120秒动态图| 精品久久久久久久久久免费视频| avwww免费| 三级毛片av免费| 国产高清不卡午夜福利| 国产高清有码在线观看视频| 久久久久免费精品人妻一区二区| АⅤ资源中文在线天堂| 亚洲精品乱码久久久v下载方式| 午夜福利高清视频| 禁无遮挡网站| 国产老妇女一区| 在线观看美女被高潮喷水网站| 午夜精品一区二区三区免费看| 亚洲国产精品sss在线观看| 真人一进一出gif抽搐免费| 伦理电影大哥的女人| 国产一区二区三区在线臀色熟女| 亚洲在线观看片| 男女下面进入的视频免费午夜| 无人区码免费观看不卡| 熟女电影av网| 久久亚洲真实| 国产午夜精品论理片| 久久人人爽人人爽人人片va| 国产 一区 欧美 日韩| 国产一区二区三区av在线 | 最近最新免费中文字幕在线| 一个人看视频在线观看www免费| 哪里可以看免费的av片| 观看免费一级毛片| 精品一区二区免费观看| 欧美3d第一页| 2021天堂中文幕一二区在线观| 桃红色精品国产亚洲av| 日韩精品中文字幕看吧| 日本免费a在线| 黄色丝袜av网址大全| 在线天堂最新版资源| 国产精品一区二区三区四区久久| 欧美色视频一区免费| 免费观看精品视频网站| 久久久久久九九精品二区国产| 床上黄色一级片| 国产美女午夜福利| 人人妻人人看人人澡| 亚洲乱码一区二区免费版| 制服丝袜大香蕉在线| 成人鲁丝片一二三区免费| 成人国产综合亚洲| 亚洲综合色惰| 成年人黄色毛片网站| av专区在线播放| 免费看a级黄色片| 免费大片18禁| 一本精品99久久精品77| 热99在线观看视频| 欧美高清性xxxxhd video| 日韩欧美精品v在线| 欧美xxxx性猛交bbbb| 夜夜看夜夜爽夜夜摸| 少妇猛男粗大的猛烈进出视频 | 亚洲三级黄色毛片| 99热网站在线观看| 99久久成人亚洲精品观看| www日本黄色视频网| xxxwww97欧美| 亚洲精品一区av在线观看| 麻豆国产av国片精品| 麻豆成人av在线观看| 91久久精品国产一区二区成人| 免费av观看视频| 床上黄色一级片| av在线老鸭窝| 精品乱码久久久久久99久播| 久久人妻av系列| 国产免费av片在线观看野外av| 精品一区二区三区av网在线观看| 午夜福利在线在线| www.www免费av| 看片在线看免费视频| 日韩精品有码人妻一区| 精品一区二区三区视频在线| 亚洲欧美日韩东京热| 国产高清视频在线观看网站| 亚洲男人的天堂狠狠| 永久网站在线| 嫩草影院新地址| 天堂网av新在线| 亚洲成人中文字幕在线播放| 国产精华一区二区三区| 国产成人一区二区在线| 色综合亚洲欧美另类图片| 精品午夜福利视频在线观看一区| 国产午夜精品久久久久久一区二区三区 | 色在线成人网| 非洲黑人性xxxx精品又粗又长| 国产精品一区二区性色av| 午夜福利视频1000在线观看| 欧美日韩黄片免| 亚洲精品久久国产高清桃花| 全区人妻精品视频| 亚洲四区av| 国产精品98久久久久久宅男小说| 日韩av在线大香蕉| 欧美不卡视频在线免费观看| 亚洲精品影视一区二区三区av| 日韩精品中文字幕看吧| 色播亚洲综合网| 国产精品乱码一区二三区的特点| 两个人视频免费观看高清| 老司机深夜福利视频在线观看| 91久久精品国产一区二区三区| 欧美最黄视频在线播放免费| 亚洲自偷自拍三级| 国内精品久久久久精免费| 亚洲av免费在线观看| 国产精品1区2区在线观看.| 99九九线精品视频在线观看视频| 国产亚洲精品av在线| 国产精品爽爽va在线观看网站| 国产激情偷乱视频一区二区| or卡值多少钱| 亚洲,欧美,日韩| 国产精品女同一区二区软件 | 日韩av在线大香蕉| av福利片在线观看| 日本一二三区视频观看| 亚洲电影在线观看av| 午夜福利视频1000在线观看| 我的老师免费观看完整版| 毛片女人毛片| 久久久久久九九精品二区国产| 午夜精品久久久久久毛片777| 免费观看在线日韩| 国产精品嫩草影院av在线观看 | 99九九线精品视频在线观看视频| 亚洲,欧美,日韩| 免费看美女性在线毛片视频| 亚洲,欧美,日韩| 在线观看美女被高潮喷水网站| 欧美日韩国产亚洲二区| 亚洲七黄色美女视频| 久久精品国产亚洲av涩爱 | 亚洲国产欧美人成| 国产精品不卡视频一区二区| 舔av片在线| 91在线观看av| 欧美日韩乱码在线| 国产精品一区二区三区四区免费观看 | 久久久久久大精品| 中国美白少妇内射xxxbb| 舔av片在线| 伦精品一区二区三区| 超碰av人人做人人爽久久| 国产成人一区二区在线| 亚洲aⅴ乱码一区二区在线播放| 日本熟妇午夜| 午夜福利在线观看吧| 麻豆成人av在线观看| 美女cb高潮喷水在线观看| 国语自产精品视频在线第100页| 欧美成人性av电影在线观看| .国产精品久久| 国模一区二区三区四区视频| 国产一区二区三区视频了| 午夜精品久久久久久毛片777| 亚洲成人久久性| 国产麻豆成人av免费视频| 久久久午夜欧美精品| 在线看三级毛片| 亚洲欧美日韩无卡精品| 中文字幕免费在线视频6| 一区福利在线观看| 日本一本二区三区精品| 国产成人福利小说| 91狼人影院| 一区二区三区四区激情视频 | 中文资源天堂在线| 两个人视频免费观看高清| 国产乱人伦免费视频| 色精品久久人妻99蜜桃| 人妻夜夜爽99麻豆av| 狂野欧美白嫩少妇大欣赏| 97碰自拍视频| 免费看日本二区| 久久精品夜夜夜夜夜久久蜜豆| 在线播放国产精品三级| 国产精品人妻久久久影院| 成人综合一区亚洲| 中文字幕熟女人妻在线| 国产精品一区二区三区四区免费观看 | av视频在线观看入口| 免费无遮挡裸体视频| 久久久久国内视频| 国产高潮美女av| 很黄的视频免费| 国产真实伦视频高清在线观看 | 美女黄网站色视频| 真人做人爱边吃奶动态| 日韩欧美免费精品| 欧美最新免费一区二区三区| 午夜精品在线福利| 国产av在哪里看| 午夜激情福利司机影院| 成人特级黄色片久久久久久久| 黄色视频,在线免费观看| 成人一区二区视频在线观看| 国产精品一及| 日本黄色片子视频| 久久国产乱子免费精品| 天堂影院成人在线观看| 免费高清视频大片| 色5月婷婷丁香| 亚洲黑人精品在线| 成年女人看的毛片在线观看| 18禁黄网站禁片午夜丰满| 午夜福利18| 久久精品国产99精品国产亚洲性色| 大型黄色视频在线免费观看| 中国美女看黄片| 国产精品99久久久久久久久| 国产一区二区激情短视频| 99久久成人亚洲精品观看| 国产精品福利在线免费观看| 99精品久久久久人妻精品| 国产一区二区三区视频了| 午夜免费激情av| 波野结衣二区三区在线| 亚洲欧美日韩卡通动漫| 成人高潮视频无遮挡免费网站| 亚洲无线在线观看| 国产探花在线观看一区二区| 97超视频在线观看视频| 在线播放国产精品三级| 国产高清激情床上av| 91麻豆av在线| 亚洲精品日韩av片在线观看| 午夜影院日韩av| 亚洲国产色片| 深夜a级毛片| 色尼玛亚洲综合影院| 精品午夜福利在线看| 日韩欧美在线乱码| 性插视频无遮挡在线免费观看| 日韩中文字幕欧美一区二区| 国产成人aa在线观看| 亚洲成人中文字幕在线播放| 国产精品一区二区免费欧美| 99九九线精品视频在线观看视频| 内射极品少妇av片p| 亚洲狠狠婷婷综合久久图片| 18+在线观看网站| 国产一区二区三区在线臀色熟女| 伊人久久精品亚洲午夜| 久久久久国内视频| 久久欧美精品欧美久久欧美| 亚洲国产精品久久男人天堂| 女生性感内裤真人,穿戴方法视频| 欧美激情国产日韩精品一区| 久久久久九九精品影院| 神马国产精品三级电影在线观看| 蜜桃久久精品国产亚洲av| 国产私拍福利视频在线观看| 亚洲中文日韩欧美视频| 两个人视频免费观看高清| 动漫黄色视频在线观看| 日本黄色片子视频| 丰满人妻一区二区三区视频av| 少妇的逼水好多| 99久久无色码亚洲精品果冻| 1000部很黄的大片| 久久久国产成人精品二区| 国产精品,欧美在线| 免费av观看视频| 亚洲五月天丁香| 亚洲专区中文字幕在线| 国产精品久久电影中文字幕| 国产一区二区亚洲精品在线观看| 99久久九九国产精品国产免费| 日本欧美国产在线视频| 国产 一区精品| 久久久久久久亚洲中文字幕| 观看美女的网站| 欧美三级亚洲精品| 桃色一区二区三区在线观看| 美女黄网站色视频| 成人永久免费在线观看视频| 99riav亚洲国产免费| 又粗又爽又猛毛片免费看| 亚洲美女搞黄在线观看 | 最后的刺客免费高清国语| 99国产极品粉嫩在线观看| 久久精品国产鲁丝片午夜精品 | 亚洲va日本ⅴa欧美va伊人久久| 高清毛片免费观看视频网站| 国产精品福利在线免费观看| 日本一本二区三区精品| 校园人妻丝袜中文字幕| 1000部很黄的大片| 亚洲人成网站在线播| 女生性感内裤真人,穿戴方法视频| 日韩精品中文字幕看吧| 国产精品一区www在线观看 | 国产在线精品亚洲第一网站| 在线免费十八禁| eeuss影院久久| 国产探花在线观看一区二区| 丝袜美腿在线中文| 成人毛片a级毛片在线播放| 欧美极品一区二区三区四区| 亚洲国产欧美人成| 99久久九九国产精品国产免费| 国产真实乱freesex| 在线观看一区二区三区| a级一级毛片免费在线观看| 丝袜美腿在线中文| 男女啪啪激烈高潮av片| 精品人妻一区二区三区麻豆 | 国产伦人伦偷精品视频| 国产极品精品免费视频能看的| 国产单亲对白刺激| 哪里可以看免费的av片| 国产乱人伦免费视频| 非洲黑人性xxxx精品又粗又长| 久久99热这里只有精品18| 亚洲午夜理论影院| 免费观看的影片在线观看| 色在线成人网| 大又大粗又爽又黄少妇毛片口| 女人十人毛片免费观看3o分钟| 国产免费一级a男人的天堂| 国产亚洲欧美98| 99热只有精品国产| 欧美绝顶高潮抽搐喷水| 成人av在线播放网站| 亚洲美女视频黄频| 色综合婷婷激情| 男人的好看免费观看在线视频| 2021天堂中文幕一二区在线观| 十八禁国产超污无遮挡网站| 午夜亚洲福利在线播放| 成人毛片a级毛片在线播放| 亚洲图色成人| 听说在线观看完整版免费高清| 日韩中文字幕欧美一区二区| 99久国产av精品| 琪琪午夜伦伦电影理论片6080| 三级男女做爰猛烈吃奶摸视频| 亚洲人成伊人成综合网2020| 88av欧美| 成年版毛片免费区| 热99re8久久精品国产| 国产精品免费一区二区三区在线| 成人永久免费在线观看视频| 国产淫片久久久久久久久| 国产精品av视频在线免费观看| netflix在线观看网站| 美女cb高潮喷水在线观看| www日本黄色视频网| 少妇的逼水好多| 99久久精品国产国产毛片| 色噜噜av男人的天堂激情| 波多野结衣巨乳人妻| 国产精华一区二区三区| 美女免费视频网站| 亚洲中文日韩欧美视频| 久久国内精品自在自线图片| 男女下面进入的视频免费午夜| 18+在线观看网站| 午夜精品在线福利| 国产精品av视频在线免费观看| 尾随美女入室| 国产真实乱freesex| 色综合亚洲欧美另类图片| 一个人看视频在线观看www免费| 久99久视频精品免费| 国产精品,欧美在线| 日韩国内少妇激情av| 精品人妻熟女av久视频| 精品乱码久久久久久99久播| 99久久九九国产精品国产免费| aaaaa片日本免费|