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

    張力腿平臺局部系泊失效模式下動力響應分析

    2019-12-21 09:03:36余建星高曉東徐立新
    關(guān)鍵詞:筋腱系泊瞬態(tài)

    余建星,高曉東,余?楊,郝?帥,徐立新

    ?

    張力腿平臺局部系泊失效模式下動力響應分析

    余建星1, 2,高曉東1, 2,余?楊1, 2,郝?帥1, 2,徐立新1, 2

    (1. 天津大學水利工程仿真與安全國家重點實驗室,天津 300350;2. 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海 200240)

    針對浮體-張力筋腱非線性耦合系統(tǒng),采用時域分析方法對局部系泊失效下的張力腿平臺進行動力響應分析,重點研究兩根筋腱失效時筋腱失效位置、波浪載荷參數(shù)與平臺各自由度動力響應的關(guān)系.結(jié)果表明:局部系泊失效對平臺橫蕩、縱蕩、艏搖自由度的動力響應影響遠小于橫搖、縱搖、垂蕩自由度;失效位置與平臺動力響應關(guān)系密切,失效筋腱位置越集中,平臺響應越大,反之亦然;筋腱失效后平臺在不同波浪載荷參數(shù)下的運動特性及頂端張力表現(xiàn)各異.垂蕩瞬態(tài)階段響應對波高和波周期敏感性較強,垂蕩穩(wěn)態(tài)階段響應對浪向、波高和波周期敏感性較強;縱搖(橫搖)瞬態(tài)和穩(wěn)態(tài)階段響應對浪向、波高和波周期敏感性均較強;波浪載荷參數(shù)對筋腱頂端張力的影響較?。?/p>

    浮體-張力筋腱耦合動動力學分析;局部系泊失效;失效位置;波浪載荷參數(shù);動力響應;時域分析

    張力腿平臺(tension leg platform,TLP)是一種垂直系泊的半順應性半剛性平臺,是由剛性系統(tǒng)(平臺上部結(jié)構(gòu)、立柱與浮箱)和彈性系統(tǒng)(張力筋腱與立管)組成的復雜非線性動力耦合系統(tǒng)[1].隨著TLP作業(yè)水深的不斷增加,筋腱載荷顯著提高,結(jié)構(gòu)系統(tǒng)運動的不穩(wěn)定性增大,極端環(huán)境載荷易導致平臺大幅升沉運動,甚至產(chǎn)生共振,這些都容易使筋腱張力超過結(jié)構(gòu)承受極限,發(fā)生瞬態(tài)失效.同時,由于累積疲勞、腐蝕缺陷、意外碰撞等產(chǎn)生的局部結(jié)構(gòu)損傷也可能會造成張力筋腱的失效,進而通過平臺主體與張力筋腱的耦合作用影響整個系統(tǒng)的性能,最終導致平臺傾覆、人員傷亡、環(huán)境污染等特大災難性事故的發(fā)生.

    為了準確預測TLP與其系泊系統(tǒng)的耦合動力響應,國內(nèi)外學者開展了大量研究[2-5],多數(shù)采用時域方法.Tahar等[6]根據(jù)非線性梁理論對系纜進行簡化,開發(fā)出一種可以模擬浮體-系纜-立管系統(tǒng)非線性動力響應的時域耦合分析方法,并進行了截斷模型實驗驗證.Ahmad等[7]和Siddiqui等[8]考慮了不同浸水、浮力和垂直方向波浪力引起的張力變化等各種非線性效應,對筋腱張力變化引起的次諧波共振以及浸水效果引起的次諧波共振進行了分析.Yang等[9-10]運用高階邊界元方法處理波浪與結(jié)構(gòu)的關(guān)系,基于細長桿理論對系纜、立管進行簡化,建立了深水浮式結(jié)構(gòu)-系泊系統(tǒng)的耦合分析方法,并研究了各結(jié)構(gòu)間的動力耦合效應.谷家揚[11]對多工況下TLP的運動耦合響應、張力筋腱和立管的張力特性進行了數(shù)值求解,并進行了模型實驗驗證(1∶40).Jayalekshmi等[12]考慮了艙體-筋腱的耦合動態(tài)響應,針對一類單柱TLP開展實驗研究,研究表明艙體-筋腱的耦合會產(chǎn)生顯著的動力放大效應,從而進一步說明耦合分析在深水工程中的重要性.

    深海作業(yè)環(huán)境異常復雜,TLP系統(tǒng)往往存在著發(fā)生各種失效的風險,學者們對于筋腱失效的瞬態(tài)響應也逐漸展開研究,Mansour等[13]分析了筋腱松弛、筋腱突然破裂以及移去單根筋腱后TLP在惡劣環(huán)境條件下的自存能力.Kim等[14]和Yang等[15]的研究主要是針對極端颶風載荷導致張力筋腱連接失效之后的平臺系統(tǒng)瞬態(tài)響應特性,并初步分析了平臺的生存性與完整性.Kim研究了筋腱在初始時刻失效和時域分析過程中失效對TLP瞬態(tài)響應的影響,并證明筋腱端部失效較底部脫落更加危險[15].然而,對于筋腱失效位置、波浪載荷參數(shù)與筋腱失效后平臺動力響應的規(guī)律性探究很少被提及.

    針對傳統(tǒng)TLP平臺,在兩根筋腱失效后剩余筋腱連續(xù)破壞的可能性激增,平臺處于危險狀態(tài),且工業(yè)界要求兩根筋腱失效后平臺仍能保證自存,因此對兩根筋腱失效后平臺的動力響應研究有重要的現(xiàn)實意義.基于此,本文建立了浮體-筋腱非線性耦合模型,計及6個自由度之間的耦合作用,在規(guī)則波下模擬兩根筋腱失效的時域響應全過程,對筋腱失效位置、組合形式與平臺動力響應的關(guān)系進行了規(guī)律性探究;并分析了平臺各階段響應的運動特性(如響應幅值、平衡位置等)及筋腱頂端張力與入射波角度、波高和波周期參數(shù)的敏感性關(guān)系.

    1?數(shù)值模型

    1.1?幾何模型

    本研究選取一個典型的四立柱式張力腿平臺“ISSC TLP”作為計算模型[16],模型由浮體結(jié)構(gòu)(立柱以及旁通)、張力腿筋腱以及海底基礎(chǔ)裝置構(gòu)成,如圖1所示.

    圖1?張力腿平臺浮體-筋腱布置

    張力腿平臺上部浮體結(jié)構(gòu)以及張力筋腱的具體參數(shù)如表1所示.

    表1?ISSCTLP主要參數(shù)

    Tab.1?Parameters of ISSC TLP

    1.2?網(wǎng)格劃分

    利用ANSYS?AQWATM軟件進行有限元計算,對TLP浮體結(jié)構(gòu)進行網(wǎng)格劃分時,一個波浪長度至少覆蓋5~8個最大尺寸單元[17],劃分網(wǎng)格時要根據(jù)需要計算的最大波浪頻率設定網(wǎng)格的控制尺寸,網(wǎng)格越細,可計算的波浪載荷頻率越大,同時對應著計算耗時的增加.本研究在保證精度的前提下選擇效率最高的劃分方式:整個浮體結(jié)構(gòu)被劃分為6,029個六面體單元,單元尺寸最大為5,m,最小為3,m.張力筋腱采用長徑比較大的柔性圓柱形管模型進行模擬,使用tether單元,每根筋腱被劃分為25個單元.

    1.3?數(shù)值分析過程

    數(shù)值分析過程如圖2所示,首先,在幾何模塊中建立浮體和系泊系統(tǒng)的幾何模型,浮體結(jié)構(gòu)由面元模型和莫里森模型組成,筋腱用tether單元進行模擬;然后,通過接觸模塊設置結(jié)構(gòu)之間的系泊屬性,在水動力輻射模塊進行頻域計算,得到附加質(zhì)量系數(shù)、附加阻尼系數(shù)、一階波浪力以及傳遞函數(shù)等,為時域分析提供基礎(chǔ)水動力數(shù)據(jù);最后,在水動力響應模塊中進行波浪載荷參數(shù)設置并完成局部筋腱失效后的時域計算.

    圖2?數(shù)值分析流程

    2?理論方法

    2.1?三維勢流理論

    基于三維勢流理論進行TLP水動力分析,假設TLP附近的流場理想,即無黏、無旋、有勢,波浪場速度勢是空間位置及時間的函數(shù),其應滿足如下邊界?條件.

    控制方程

    線性化自由面條件

    水底條件

    物面條件

    在直角坐標系中,沿波傳播方向,沿靜止水面法向;是水深參數(shù);是物體的運動速度;是物面上的單位外法向矢量.

    求出速度勢之后,根據(jù)Bernoulli方程計算平臺濕表面上的壓力分布.

    2.2?時域運動方程

    在時域分析中,張力腿平臺的運動方程為

    2.3?筋腱模型

    本研究采用tether模型模擬實際作業(yè)中的TLP筋腱,每個單元受力包括3部分.

    3?模型驗證

    本研究使用的張力腿平臺模型來源于Gie和Boom的張力腿平臺縮尺比實驗[16].包括英國海事技術(shù)有限公司、布朗根(英國)有限公司、英國倫敦大學學院等在內(nèi)的17家組織機構(gòu)對該TLP做過頻域分析,計算得到的水動力參數(shù)與本文數(shù)值模擬得到的結(jié)果幾乎相同.頻域計算得到的水動力數(shù)據(jù)是時域計算的基礎(chǔ),因此充分證明了本次數(shù)值模擬的準確性.此前17家組織計算的水動力參數(shù)包括0°、22.5°、45°浪向下各自由度的RAO(幅值響應算子)、附加質(zhì)量、輻射阻尼以及平均漂移力等[16].因篇幅有限,選取部分對比數(shù)據(jù)如圖3和圖4所示.

    圖3?22.5°浪向下垂蕩響應模量

    圖4?0°浪向下縱蕩平均漂移力

    4?計算結(jié)果及討論

    4.1?不同數(shù)量筋腱失效后平臺動力響應分析

    平臺筋腱位置編號及波浪方向如圖5所示,定義波浪首先接觸的位置為迎浪位置,反之為背浪位置.針對本典型張力腿平臺,分別建立1根、2根、3根筋腱失效模型,在規(guī)則波下進行時域計算分析,取規(guī)則波波高6,m、周期10,s、浪向0°,設置筋腱在297.2,s同時失效(筋腱張力達到最大時),以垂蕩響應為例,對3種失效工況進行對比分析,如圖6所示.

    圖5?0°浪向下筋腱失效位置示意

    圖6?不同數(shù)量筋腱失效后平臺垂蕩時歷曲線

    當1根筋腱失效后,平臺瞬態(tài)響應最大值為3.22,m,較初始響應幅值僅增大了5.57%,.在整個時域響應過程中,剩余筋腱張力值均遠小于極限張力值,故不會發(fā)生連續(xù)失效,平臺仍處于安全狀態(tài).當2根筋腱失效后,系統(tǒng)垂蕩自由度剛度嚴重缺失,平臺產(chǎn)生了較大的瞬態(tài)響應,響應最大值較初始幅值增大了30.2%,;瞬態(tài)響應之后的穩(wěn)態(tài)平衡位置和響應幅值都有明顯提高,靠近失效位置的筋腱張力激增,接近其極限破斷張力,平臺處于臨界危險狀態(tài).當3根筋腱失效后,平臺垂蕩自由度產(chǎn)生劇烈響應,瞬態(tài)響應幅值較初始幅值增大了56.4%,,且緊鄰失效位置的#4筋腱頂端張力在筋腱失效瞬間激增,遠超過極限張力值,將會引發(fā)連續(xù)失效,平臺終將失去自存能力.綜上所述,1根筋腱失效后平臺處于較安全狀態(tài),2根筋腱失效后平臺處于臨界危險狀態(tài),3根筋腱失效后平臺將連續(xù)破壞直至傾覆.通常工業(yè)界要求平臺設計時需具備2根筋腱失效之后TLP仍能自存的能力,因此本文重點討論臨近危險狀態(tài)即2根筋腱失效后平臺的動力響應特征.

    4.2?筋腱失效位置與平臺動力響應的關(guān)系

    由于長期作業(yè),實際在位的TLP某些筋腱處于疲勞損傷狀態(tài),儲備張力值遠低于設計張力.在極端環(huán)境載荷下,筋腱失效會導致剩余筋腱張力水平在短時間內(nèi)迅速增加,誘發(fā)連續(xù)失效.當2根筋腱連續(xù)失效的過程極其短暫時,可以認為失效同時發(fā)生,這具有很大的危險性.本節(jié)將不同位置的2根筋腱組合,對TLP進行局部失效后的時域分析,研究筋腱失效位置與張力腿平臺各自由度動力響應的關(guān)系.6組不同位置的失效組合具體工況分類見表2.

    表2?筋腱失效組合位置

    Tab.2?Classification of tendon failure position

    4.2.1?規(guī)則波下局部筋腱失效后平臺動力響應分析

    張力腿平臺在平面內(nèi)的運動響應(縱蕩、橫蕩和艏搖)具有柔性特點,即順應性.由于TLP結(jié)構(gòu)的對稱性,局部筋腱失效后,平臺橫蕩和縱蕩瞬態(tài)響應規(guī)律一致,只需考慮其一即可.

    計算選取規(guī)則波波高為8,m、周期為12,s,經(jīng)試算選擇時域時長1,000,s,計算步長0.1,s,筋腱失效時刻為295,s(筋腱張力達到最大時).圖7是部分筋腱組合失效后TLP橫蕩自由度的時歷響應曲線,將各種工況下的最大響應值變化量統(tǒng)計如表3所示.

    圖7?局部系泊失效過程TLP橫蕩時歷曲線

    由表3可知,筋腱失效后TLP橫蕩響應最大值較初始響應幅值的最大偏移量為0.5,m(工況1),僅占作業(yè)水深的0.11%,.結(jié)合橫蕩自由度順應性特征,該瞬態(tài)效應引起的位移變化可以忽略不計.值得注意的是,#3#4筋腱失效后,TLP橫蕩最大響應出現(xiàn)在303,s,距離失效時刻僅8,s;但#1#8筋腱失效后,TLP橫蕩最大響應并沒有在失效瞬間發(fā)生,而是出現(xiàn)在距離失效時刻約50,s的時間節(jié)點處.因此對于橫蕩自由度,TLP在局部失效后的瞬態(tài)響應不明顯,且部分失效組合引起的瞬態(tài)效果在時間上具有一定滯后性.

    表3?橫蕩和艏搖響應統(tǒng)計值

    Tab.3?Statistics of sway and yaw responses

    圖8是部分筋腱組合失效后TLP艏搖時歷響應曲線,各種失效工況下的艏搖最大響應值見表3.對比各失效工況下的艏搖響應,發(fā)現(xiàn)#3#4、#2#3筋腱組合失效后TLP艏搖響應發(fā)生激增,響應值遠大于其他失效組合.原因為當上述筋腱組合失效后,在筋腱失效位置、波浪作用方向等一系列非對稱因素作用下,平臺在橫蕩和縱蕩自由度產(chǎn)生了非對稱響應,而艏搖與縱蕩、橫蕩有明顯的耦合效應,間接導致艏搖響應產(chǎn)生跨量級的激增.與橫蕩自由度類似,#3#4筋腱失效后,艏搖響應最大值出現(xiàn)在350,s左右,說明部分失效組合引起的瞬態(tài)效果具有一定滯后性.

    圖8?局部系泊失效過程TLP艏搖時歷曲線

    張力腿平臺的筋腱系統(tǒng)往往在預張力作用下時刻處于張緊狀態(tài),巨大的張力使平臺豎直面上運動響應較小,幾乎近似于剛性.因為筋腱為剛性自由度提供了直接約束,所以筋腱失效會直接影響橫搖、縱搖和垂蕩的動力響應,這將是本文的研究重點.

    圖9、圖10是在0°浪向規(guī)則波下不同筋腱組合失效后TLP垂蕩、縱搖響應時歷曲線.

    現(xiàn)將時域響應過程進行劃分,以#3#4筋腱失效后垂蕩響應為例:0~295,s為筋腱系統(tǒng)完好的初始響應階段;295~315,s為筋腱失效后的瞬態(tài)響應階段,在此期間響應發(fā)生較大波動;315,s以后為穩(wěn)態(tài)響應階段,響應趨于穩(wěn)定平衡.通過研究對比不同響應階段的運動特性(如響應幅值、平衡位置等)來更全面地分析響應規(guī)律.

    圖9?局部系泊失效過程TLP垂蕩時歷曲線

    圖10?局部系泊失效過程TLP縱搖時歷曲線

    為了清晰全面地表現(xiàn)出局部筋腱失效后平臺動力響應大小,設置如下無量綱參數(shù):

    就垂蕩和縱搖自由度而言,同一立柱下的2根筋腱失效響應參數(shù)最大,同一旁通上的2根筋腱失效響應參數(shù)次之,對角線立柱下的筋腱失效響應參數(shù)最?。驗楫敳煌铍旖M合失效后平臺垂蕩剛度相近,但縱搖剛度差別迥異,其中同一立柱下筋腱組合失效使平臺一角約束缺失,縱搖剛度產(chǎn)生嚴重折減,因此在波浪載荷作用下平臺產(chǎn)生較強的縱搖響應.垂蕩與縱搖有強耦合作用,所以也表現(xiàn)出明顯的動力響應,使TLP處于最危險狀態(tài).

    表4?垂蕩和縱搖響應參數(shù)統(tǒng)計

    Tab.4?Parameters of heave and pitch responses

    #1#8失效工況的垂蕩響應參數(shù)均大于#2#7失效工況,說明同一旁通的外側(cè)筋腱較內(nèi)側(cè)筋腱對TLP垂蕩響應的影響更大.原因為#1#8筋腱失效之后,TLP會繞著次轉(zhuǎn)軸(次轉(zhuǎn)軸平行且更靠近迎浪旁通)轉(zhuǎn)動,#2#7筋腱失效亦然,但#2#7筋腱到次轉(zhuǎn)軸的距離較#1#8筋腱更近,所以局部失效之后,#2#7筋腱張力水平高于#1#8筋腱.基于小變形前提,考慮筋腱的軸向剛度相同,#2#7筋腱的伸長量大于#1#8筋腱,也就是說#1#8筋腱組合失效引起的垂蕩響應幅值大于#2#7筋腱組合失效.

    #1#8失效后的響應參數(shù)和#4#5失效后的響應參數(shù)相差均在0.03%,之內(nèi),說明失效筋腱所處迎浪位置和背浪位置對TLP垂蕩響應影響不大.

    由于垂蕩和縱搖、橫搖的耦合關(guān)系,筋腱失效位置和TLP縱搖響應參數(shù)變化規(guī)律與垂蕩相似,但存在特殊情況,從表4可知#2#3失效后縱搖響應參數(shù)很小,因為#2#3筋腱失效后暫存筋腱組成的系泊系統(tǒng)仍關(guān)于軸對稱,結(jié)構(gòu)對稱性未被破壞,整個系統(tǒng)縱搖剛度較大,所以TLP縱搖響應不明顯.因此,同一旁通下筋腱失效,若筋腱失效位置關(guān)于縱(橫)蕩自由度對稱,則TLP縱(橫)搖響應很小.

    4.2.2?不規(guī)則波下局部筋腱失效后平臺的動力響應分析

    考慮中國南海百年一遇的海況(表5),設置系泊系統(tǒng)兩根筋腱在為清晰量化平臺響應大小,綜合對比不同階段運動特性參數(shù):瞬態(tài)響應最大值,穩(wěn)態(tài)響應幅值以及穩(wěn)態(tài)平衡位置.由圖11(a)可知,#3#4筋腱(相同立柱)失效后,垂蕩響應參數(shù)最大,其中瞬態(tài)響應最大值較初始響應幅值增大了27%,,TLP處于最危險狀態(tài);#4#5筋腱(相同旁通)失效后,垂蕩響應居中;#1#6筋腱(對角線立柱)失效后,垂蕩響應最?。蓤D11(b)可知,在不規(guī)則波下,#1#8筋腱失效后的垂蕩響應參數(shù)均大于#2#7筋腱失效,可知在不規(guī)則波作用下外側(cè)筋腱較內(nèi)側(cè)筋腱對TLP垂蕩影響更大.

    表5?中國南海百年一遇海況參數(shù)

    Tab.5 Parameters of 100-year sea conditions in South China Sea

    圖11不規(guī)則波下局部筋腱失效后TLP垂蕩時歷曲線

    由于垂蕩和縱搖的耦合作用,局部失效后平臺縱搖響應規(guī)律與垂蕩類似.由圖12(a)可知,對于縱搖,#3#4筋腱失效后,TLP縱搖響應最大;#4#5筋腱失效后,縱搖響應居中;#1#6筋腱失效后,TLP縱搖響應與完整筋腱系統(tǒng)響應相差無幾,影響較?。蓤D12(b)可知,#1#8筋腱失效后,縱搖響應大于#2#7筋腱失效;#2#3筋腱失效后,平臺縱搖響應很小,因為失效后平臺縱搖剛度較大,響應不明顯.

    圖12?不規(guī)則波下局部筋腱失效后TLP縱搖時歷曲線

    綜上對比可知,規(guī)則波下得到的筋腱失效位置與TLP平臺垂蕩和縱搖(橫搖)的瞬態(tài)響應規(guī)律在不規(guī)則波作用下同樣適用.

    4.3?波浪載荷參數(shù)與筋腱失效后平臺動力響應關(guān)系

    通過研究筋腱失效位置與平臺動力響應的關(guān)系,發(fā)現(xiàn)筋腱失效后平臺橫蕩、縱蕩、艏搖自由度的瞬態(tài)響應效果遠小于垂蕩、橫搖、縱搖自由度,因此,在研究波浪載荷參數(shù)與筋腱失效后平臺動力響應的關(guān)系時,重點考慮垂蕩、橫搖和縱搖自由度,并將其作為主控自由度.

    計算工況由不同波高(=6,m,8,m,10,m)、不同波周期(=10,s,12,s,14,s)及不同入射波角度?(=0°,10°,20°)組合而成,如A0H6T10表示入射波角度0°、波高6,m、波周期10,s的波浪載荷;筋腱失效組合選擇最危險的失效模式,相同立柱下的兩根筋腱同時失效.本文對不同入射波角度、波高和波周期展開討論.

    4.3.1?不同浪向?qū)铍焓Ш笃脚_動力響應的影響

    選取規(guī)則波的波高為6,m,波周期為10,s,入射角度分別取0°、10°、20°;設置#1#2筋腱在297.2,s失效.

    因為張力腿平臺結(jié)構(gòu)的對稱性,平臺橫搖和縱搖響應有相似的變化規(guī)律,所以只需分析固定浪向下主控自由度的響應規(guī)律即可.在不同浪向下,對局部失效后平臺進行時域響應求解,并通過快速傅里葉變換(FFT)對穩(wěn)態(tài)響應階段進行頻譜分析,平臺垂蕩和縱搖的動力響應結(jié)果如圖13和圖14所示.

    根據(jù)局部系泊失效后平臺時域響應過程的劃分,重點對比不同響應階段的運動參數(shù):瞬態(tài)響應階段的最大值、穩(wěn)態(tài)響應階段的幅值以及平衡位置.

    圖13的計算結(jié)果表明:在筋腱完好的初始階段,不同浪向下平臺的垂蕩響應幅值和平衡位置幾乎相同.局部筋腱失效后,在瞬態(tài)響應階段,平臺運動最大值相同,均為3.92,m;在穩(wěn)態(tài)響應階段,隨著浪向的增大,平臺的響應幅值逐漸增大,最大相差約75.1%,;平臺在穩(wěn)態(tài)響應階段的平衡位置保持不變;通過快速傅里葉變換對穩(wěn)態(tài)響應階段進行頻譜分析,發(fā)現(xiàn)存在較強的2倍頻成分以及微弱的3倍、4倍頻成分,局部系泊失效后TLP運動穩(wěn)定性下降.隨著入射波角度的增加,2倍頻分量幅值增長明顯.綜上可知,浪向?qū)植拷铍焓Ш笃脚_穩(wěn)態(tài)階段響應有影響,對瞬態(tài)階段響應影響較小.

    圖13?不同浪向下筋腱失效過程TLP垂蕩動力響應

    從圖14可知,筋腱未失效時,入射波角對縱搖的影響甚微.局部筋腱失效后,隨著入射波角度的增大,瞬態(tài)響應階段最大響應值逐漸增大;在穩(wěn)態(tài)響應階段中外激勵頻率分量隨著入射角度增加而不斷增加,幅值最大相差87.7%,.雖出現(xiàn)一定量的二倍頻成分,但其與入射波角度相關(guān)性較小.因此,局部失效后平臺縱搖各階段響應對浪向參數(shù)的敏感度較強,且與垂蕩響應變化規(guī)律又不盡相同.

    圖14?不同浪向下筋腱失效過程TLP縱搖動力響應

    因為TLP系泊系統(tǒng)為垂蕩、縱(橫)搖自由度提供較強的剛性約束,所以當局部筋腱失效后,垂蕩和縱(橫)搖自由度剛度較初始剛度均大幅減小.當浪向改變時,波浪力沿垂向分量變化不大,所以垂蕩響應對浪向參數(shù)敏感性較弱;而波浪力的水平分量隨著浪向的改變而產(chǎn)生明顯改變,所以縱搖響應對浪向表現(xiàn)出較強的敏感性.

    4.3.2?不同波高對筋腱失效后平臺動力響應的影響

    選取規(guī)則波的浪向為0°,波周期為10,s,波高分別取6,m、8,m、10,m.在0°浪向下,垂蕩與縱搖是主控搖動自由度.不同波高下,平臺在局部系泊失效后垂蕩和縱搖動力響應結(jié)果如圖15和圖16所示.

    從圖15可知,在筋腱完好的初始階段,隨著波高的增大,垂蕩幅值逐漸增大,垂蕩平衡位置逐漸減?。植拷铍焓Ш?,在瞬態(tài)響應階段,隨著波高的增大,運動最大值逐漸增大.在穩(wěn)態(tài)響應階段出現(xiàn)的2倍頻成分較明顯,外激勵頻率分量與2倍頻分量隨著波高提升,均不斷增長,穩(wěn)態(tài)響應幅值最大相差93.3%,.穩(wěn)態(tài)平衡位置保持不變,且存在一定的相位差,這是由于不同波高產(chǎn)生不同垂向分力造成的.波高越大,穩(wěn)態(tài)階段響應幅值較初始階段幅值的增加量越大,這是由于系泊系統(tǒng)破壞使得垂向約束減小所致.因此,局部筋腱失效后垂蕩各階段響應對波高敏感度均較強.

    圖15?不同波高下筋腱失效過程TLP垂蕩動力響應

    如圖16所示,在初始階段,不同波高下的縱搖平衡位置幾乎不變,但響應幅值隨著波高增大而增大.在局部筋腱失效之后,隨著波高的增大,瞬態(tài)響應最大值增大.穩(wěn)態(tài)響應幅值增大,最多相差72.7%,,但平衡位置幾乎不變.與垂蕩響應類似,穩(wěn)態(tài)響應也出現(xiàn)較強的2倍頻成分,且外載荷頻率分量與2倍頻分量均與波高有一定的相關(guān)性.由于縱搖和垂蕩自由度的強耦合性,波高對局部失效后平臺各階段縱搖響應的影響與垂蕩相似,但敏感性比垂蕩稍微弱一些.

    出現(xiàn)上述規(guī)律原因為局部系泊失效會使平臺各自由度的剛度有不同程度的折損.針對本次研究的典型TLP模型,當#1#2筋腱失效后,平臺垂蕩自由度的剛度較初始剛度減小比重大于縱搖自由度,即平臺垂蕩自由度的回復能力損失較縱搖更嚴重,導致平臺垂蕩響應更劇烈,因此垂蕩對波高參數(shù)的敏感性略強于縱搖.

    圖16?不同波高下筋腱失效過程TLP縱搖動力響應

    4.3.3?不同波浪周期對筋腱失效后平臺動力響應的影響

    選取規(guī)則波的浪向為0°,波高為6,m,波浪周期分別取10,s、12,s、14,s.不同波浪周期下,平臺在局部筋腱失效后垂蕩和縱搖動力響應結(jié)果如圖17和圖18所示.

    不同波浪周期下,在初始階段平臺垂蕩與縱搖的變化極為相似:時歷曲線都存在相位差;平衡位置都基本不變.如圖17所示,局部筋腱失效后,隨著波浪周期的增大,瞬態(tài)階段垂蕩響應最大值越來越?。M入穩(wěn)態(tài)階段,垂蕩平衡位置幾乎保持不變,穩(wěn)態(tài)幅值最大相差114.3%,,且與波浪周期并不是線性遞變,3種工況中周期為12,s的波浪載荷對應的響應幅值最大,平臺響應出現(xiàn)了較強的2倍頻成分,周期為10,s載荷的影響次之,周期為14,s載荷的影響最?。c完整平臺的在不同波周期下的運動幅值變化(如圖3)相比,周期為10,s的波浪載荷下運動幅值上升明顯.原因為:完整筋腱系統(tǒng)下平臺的垂蕩固有周期約為3.4,s,兩根筋腱失效后,垂蕩自由度的剛度減小,對與其周期接近的波浪載荷會產(chǎn)生一定的動力放大效果.綜上所述,局部筋腱失效后平臺垂蕩各階段響應對波浪周期的敏感度較強,且接近失效系統(tǒng)垂蕩固有周期的波浪載荷在穩(wěn)態(tài)階段更容易產(chǎn)生較大的響應幅值.

    就縱搖自由度而言,筋腱失效后,隨著波浪周期增大,縱搖瞬態(tài)響應最大值減小,最大相差20.3%,;不同波浪周期的穩(wěn)態(tài)響應幅值不同,最大相差79.3%,,平衡位置相差無幾;在頻譜分析中,不同波浪周期載荷參數(shù)下平臺縱搖運動均出現(xiàn)了不同程度的二倍頻成分,其規(guī)律與垂蕩基本一致.因此,波浪周期對局部失效后縱搖各階段響應也有較大的影響.

    圖17?不同波周期下筋腱失效過程TLP垂蕩動力響應

    圖18?不同波周期下筋腱失效過程TLP縱搖動力響應

    4.3.4?不同波浪參數(shù)對筋腱頂端張力的影響

    相對于平臺的運動響應,筋腱張力的變化更加直接反映結(jié)構(gòu)的危險狀態(tài).本次研究的筋腱采用X65鋼管(參數(shù)見表1),預張力值是1.7×107,N,在一根或多根筋腱失效后,與其緊鄰的筋腱張力會激增,產(chǎn)生連續(xù)失效破壞,計算設置#1#2筋腱失效,研究不同波浪參數(shù)與#3、#8筋腱(緊鄰失效筋腱)頂端張力最大值的關(guān)系.

    圖19是不同波浪載荷參數(shù)下#3和#8筋腱在整個時域分析過程中頂端張力最大值統(tǒng)計圖,結(jié)果顯示:局部筋腱失效后,不同入射角度下#3和#8筋腱頂端最大張力值幾乎沒有差別,也就是說,筋腱頂端張力對浪向敏感性很??;隨著波高的增大,#3和#8筋腱頂端張力值逐漸增大;隨著波周期的增大,#3和#8筋腱頂端張力值也逐漸遞減.但是總體而言,不同波高下筋腱張力值最大僅相差7.9%,,不同波周期下筋腱張力值最大相差9.2%,,說明不同波浪載荷參數(shù)對筋腱頂端張力的影響較?。?/p>

    圖19?不同波浪載荷參數(shù)下筋腱頂端張力統(tǒng)計值

    圖20是不同波浪載荷參數(shù)下平臺垂蕩瞬態(tài)響應最大值統(tǒng)計,可以發(fā)現(xiàn),在不同浪向、波高和波周期參數(shù)下,垂蕩瞬態(tài)響應最大值與剩余筋腱頂端最大張力值的變化規(guī)律完全一致,充分體現(xiàn)了垂蕩與頂端張力的強相關(guān)性.

    圖20?不同波浪載荷參數(shù)下平臺垂蕩響應統(tǒng)計值

    5?結(jié)?論

    針對擁有8根筋腱的經(jīng)典張力腿平臺結(jié)構(gòu),在規(guī)則波下采用時域分析方法,對局部系泊失效后張力腿平臺動力響應和系泊特性開展了研究,主要結(jié)論如下.

    (1) 由于TLP在平面內(nèi)的運動響應具有順應性特征,局部系泊失效后,橫蕩、縱蕩和艏搖自由度的瞬態(tài)響應并不明顯,部分筋腱組合失效引起的瞬態(tài)效果在時間上具有一定的滯后性.

    (2) 系泊系統(tǒng)為垂蕩、橫搖和縱搖自由度提供主要剛性約束,所以局部系泊失效后,這3個自由度的瞬態(tài)效果明顯;失效筋腱位置越集中,TLP響應越大,失效筋腱位置越分散,TLP響應越?。?/p>

    (3) 系泊系統(tǒng)完好時,TLP艏搖響應遠小于橫搖、縱搖響應;部分筋腱組合失效后,由于筋腱失效位置及環(huán)境載荷作用方向等不對稱因素,TLP艏搖響應最大值激增,數(shù)值接近橫搖和縱搖響應的數(shù)量級.

    (4) 局部系泊失效后,平臺垂蕩瞬態(tài)階段響應對波高和波周期敏感性較強,穩(wěn)態(tài)階段響應對浪向、波高和波周期敏感性較強;縱搖(橫搖)瞬態(tài)和穩(wěn)態(tài)階段響應對浪向、波高和波周期敏感性均較強.

    (5) 波高越大,波浪周期越小,在時域分析中剩余筋腱頂端張力最大值越大,但總體而言,波浪參數(shù)對筋腱頂端張力的影響較?。唤铍祉敹藦埩εc平臺垂蕩響應有強相關(guān)性.

    [1] API(American Petroleum Institute). Draft Recommended Practice 2,T,Rev. 3. Planning,Designing,and Constructing Tension Leg Platforms[S]. Washington,DC:API,2007.

    [2] Jameel M,Oyejobi D O,Siddiqui N A,et al. Nonlinear dynamic response of tension leg platform under environmental loads[J]. KSCE Journal of Civil Engineering,2017,21(3):1022-1030.

    [3] Oyejobi D O,Jameel M,Sulong N H R. Nonlinear response of tension leg platform subjected to wave,current and wind forces[J]. International Journal of Civil Engineering,2016,14(8,A):521-533.

    [4] 谷家揚,楊建民,呂海寧. 考慮風浪流的TLP張力腿張力響應規(guī)律分析[J]. 船舶力學,2012,16(9):1041-1055. Gu Jiayang,Yang Jianmin,Lü Haining. Comprehensive study on tendon tension of TLP considering wind,wave and current[J]. Journal of Ship Mechanics,2012,16(9):1041-1055(in Chinese).

    [5] Gu Jiayang,Yang Jianmin,Lu Haining. Studies of TLP dynamic response under wind,waves and current[J]. China Ocean Engineering,2012,26(3):363-378.

    [6] Tahar A,Kim M H. Hull/mooring/riser coupled dynamic analysis and sensitivity study of a tanker-based FPSO [J]. Applied Ocean Research,2003,25(6):367-382.

    [7] Ahmad S. Stochastic TLP response under long crested random sea[J]. Computers & Structure,1996,61(6):975-993.

    [8] Siddiqui N A,Ahmad S. Fatigue and fracture reliability of TLP tethers under random loading[J]. Marine Structures,2001,14(3):331-352.

    [9] Yang Mindong,Teng Bin,Xiao Longfei. Full time-domain nonlinear coupled dynamic analysis of a truss spar and its mooring/riser system in irregular wave[J]. Science China Physics Mechanics & Astronomy,2014,57(1):152-165.

    [10] Yang Mindong,Teng Bin,Ning Ding,et al. Coupled dynamic analysis for wave interaction with a truss spar and its mooring line/riser system in time domain[J]. Ocean Engineering,2012,39:72-87.

    [11] 谷家揚. 張力腿平臺復雜動力響應及渦激特性研究[D]. 上海:上海交通大學船舶海洋與建筑工程學院,2013. Gu Jiayang. Study on the Complex Dynamic Response and Vortex-Induced Motion Characteristics of Tension Leg Platform[D]. Shanghai:School of Naval Architecture,Ocean and Civil Engineering,Shanghai Jiao Tong University,2013(in Chinese).

    [12] Jayalekshmi R,Idichandy V G,Sundaravadivelu R. Physical simulation of the hull-tether coupled dynamics of a deepwater TLP[C]// ASME 26th International Conference on Offshore Mechanics and Arctic Engineering 2007. Michigan:American Society of Mechanical Engineers,2007:149-157.

    [13] Mansour A M,Huang E W,Phadke A C,et al. Tension leg platform survivability analysis[C]//ASME 25th International Conference on Offshore Mechanics and Arctic Engineering. Michigan:American Society of Mechanical Engineers,2006:477-483.

    [14] Kim M H,Zhang Zhi. Transient effects of tendon disconnection of a TLP in moderate-strength hurricane conditions[J]. International Journal of Naval Architecture and Ocean Engineering,2009,1(1):13-19.

    [15] Yang Chan K,Kim M H. Transient effects of tendon disconnection of a TLP by hull-tendon-riser coupled dynamic analysis[J]. Ocean Engineering,2010,37(8/9):667-677.

    [16] Eatock Taylor R,Jefferys E R. Variability of hydrodynamic load predictions for a tension leg platform[J]. Ocean Engineering,1986,13(5):449-490.

    [17] Liu Weimin. Research on time-domain dynamic response of tension leg platform in regular wave[J]. Applied Mechanics and Materials,2014,670/671:801-804,

    [18] Faltinsen Odd M. Sea Loads on Ships and Offshore Structures[M]. United Kingdom:Cambridge University Press,1993.

    [19] Chakrabarti S K. Hydrodynamics of Offshore Structures [M]. United Kingdom:WIT Press,1987.

    [20] API(American Petroleum Institute). API Recommended Practice 2,SK,Design and Analysis of Station Keeping Systems for Floating Structures[S]. 3rd Ed. Washington,DC:API,2005.

    Dynamic Response Analysis of a Tension Leg Platform with Multiple Broken Tendons

    Yu Jianxing1, 2,Gao Xiaodong1, 2,Yu Yang1, 2,Hao Shuai1, 2,Xu Lixin1, 2

    (1.State Kay Laboratory of Hydraulic Engineering Simulation and Safety,Tianjin University,Tianjin 300350,China;2.Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration,Shanghai 200240,China)

    In terms of the nonlinear coupling effect of a tension leg platform(TLP) hull and tendons,the dynamic response of a TLP with multiple broken tendons was analyzed using a time-domain method.The relationship among tendon failure positions,wave force parameters,and dynamic responses of six degrees of freedom was then investigated.It was found that with multiple broken tendons,the dynamic responses of sway,surge and yaw were far smaller than those of roll,pitch and heave.The tendon failure positions had close correlations with the TLP dynamic responses.The more centralized the tendon failure positions were,the higher the TLP responses were,and vice versa.As wave parameters varied,the TLP performance varied across different motion characteristics and at peak tension.The transient heave response was sensitive to wave height and wave period,while its steady-state response was sensitive to wave direction,height,period.However,the transient and steady-state responses in pitch and roll were also sensitive to wave direction,height,and period.The wave force parameters had small impact at top tension of tendon.

    hull-tendon coupled dynamic analysis;multiple broken tendons;failure position;wave force parameter;dynamic response;time-domain analysis

    P751

    A

    0493-2137(2019)01-0040-12

    2018-04-11;

    2018-05-16.

    余建星(1958—??),男,教授,yjx2000@tju.edu.cn.

    余?楊,yang.yu@tju.edu.cn.

    國家重點基礎(chǔ)研究發(fā)展計劃(973計劃)資助項目(2014CB046803);國家自然科學基金青年基金資助項目(51609169);國家自然科學基金資助項目(51779173);廣西科技重大專項項目(桂科AA17292007).

    the National Basic Research Program of China(No.,2014CB046803),the National Natural Science Foundation of China for Youth(No.,51609169),the National Natural Science Foundation of China(No.,51779173),the Guangxi Science and Technology Major Project(No.,Guike AA17292007).

    10.11784/tdxbz201804032

    (責任編輯:王新英)

    猜你喜歡
    筋腱系泊瞬態(tài)
    基于CFD方法的系泊狀態(tài)FSO舵載荷計算
    高韌性海洋用三級系泊鏈鋼的研制
    山東冶金(2022年2期)2022-08-08 01:50:50
    筋腱斷裂后張力腿平臺-筋腱-立管耦合系統(tǒng)響應分析
    高壓感應電動機斷電重啟時的瞬態(tài)仿真
    防爆電機(2020年3期)2020-11-06 09:07:36
    張力腿平臺筋腱安裝與臨時浮筒設計
    張力腿筋腱樁基水平安裝誤差影響分析
    船海工程(2017年5期)2017-11-01 12:46:37
    基于懸鏈線方程的系泊系統(tǒng)分析
    十億像素瞬態(tài)成像系統(tǒng)實時圖像拼接
    中國光學(2015年5期)2015-12-09 09:00:39
    基于瞬態(tài)流場計算的滑動軸承靜平衡位置求解
    DC/DC變換器中的瞬態(tài)特性分析
    亚洲第一av免费看| 美女cb高潮喷水在线观看| 国产毛片在线视频| 乱人伦中国视频| 黄色视频在线播放观看不卡| 18禁在线无遮挡免费观看视频| 久久久久久久久大av| 狠狠婷婷综合久久久久久88av| 99国产精品免费福利视频| 精品人妻熟女毛片av久久网站| 18禁观看日本| 国产亚洲av片在线观看秒播厂| 中文字幕免费在线视频6| 亚洲精品成人av观看孕妇| 最后的刺客免费高清国语| 久久久久久久精品精品| av天堂久久9| 一区二区三区四区激情视频| 丰满乱子伦码专区| av福利片在线| 哪个播放器可以免费观看大片| 国产成人91sexporn| 18禁裸乳无遮挡动漫免费视频| 中国三级夫妇交换| 欧美另类一区| 亚洲欧美日韩卡通动漫| 一级a做视频免费观看| 18禁裸乳无遮挡动漫免费视频| 在线观看免费日韩欧美大片 | 久久久久精品久久久久真实原创| 久久久久久久久久成人| 最近的中文字幕免费完整| 一边摸一边做爽爽视频免费| 一区二区三区免费毛片| av福利片在线| 亚洲美女黄色视频免费看| 国产免费福利视频在线观看| 高清午夜精品一区二区三区| 五月开心婷婷网| 91精品一卡2卡3卡4卡| 一区二区日韩欧美中文字幕 | 欧美日韩在线观看h| 欧美精品人与动牲交sv欧美| 亚洲精品亚洲一区二区| 黄色配什么色好看| 国产欧美亚洲国产| 亚洲人成77777在线视频| 亚洲欧美精品自产自拍| 亚洲精品视频女| 国产黄色免费在线视频| 亚洲天堂av无毛| 日韩大片免费观看网站| 男女高潮啪啪啪动态图| 水蜜桃什么品种好| 国内精品宾馆在线| 91精品国产国语对白视频| 性色avwww在线观看| 美女内射精品一级片tv| 精品久久久精品久久久| 免费观看无遮挡的男女| 一二三四中文在线观看免费高清| 日韩制服骚丝袜av| 精品久久久久久电影网| 亚洲国产精品一区二区三区在线| 日韩一区二区视频免费看| 下体分泌物呈黄色| 一区二区三区免费毛片| 国产高清有码在线观看视频| 日韩av不卡免费在线播放| 国产精品成人在线| 欧美精品一区二区大全| 九九在线视频观看精品| 日韩一区二区视频免费看| 国产欧美日韩一区二区三区在线 | 欧美老熟妇乱子伦牲交| 久久毛片免费看一区二区三区| 国产日韩欧美在线精品| 日本欧美视频一区| 秋霞伦理黄片| 99re6热这里在线精品视频| 欧美+日韩+精品| 午夜影院在线不卡| 青春草视频在线免费观看| 18+在线观看网站| 亚洲精品一二三| 国产毛片在线视频| 亚洲在久久综合| 考比视频在线观看| 亚洲性久久影院| 两个人免费观看高清视频| 尾随美女入室| 中文字幕久久专区| 伊人久久精品亚洲午夜| 五月开心婷婷网| 久久国产精品大桥未久av| 成年美女黄网站色视频大全免费 | 免费久久久久久久精品成人欧美视频 | 国产精品人妻久久久影院| 日韩大片免费观看网站| 国产一区二区在线观看日韩| 美女国产视频在线观看| 丰满乱子伦码专区| 国产片内射在线| 性色avwww在线观看| 亚洲av免费高清在线观看| 下体分泌物呈黄色| 纯流量卡能插随身wifi吗| 老女人水多毛片| 国产日韩欧美亚洲二区| 日本91视频免费播放| 少妇精品久久久久久久| 中文字幕久久专区| 免费黄频网站在线观看国产| 蜜桃在线观看..| 久久久久久久久久人人人人人人| 国产熟女欧美一区二区| 免费人妻精品一区二区三区视频| 国产成人精品久久久久久| 只有这里有精品99| 久久人妻熟女aⅴ| 亚洲av欧美aⅴ国产| 视频区图区小说| 国产亚洲最大av| av国产久精品久网站免费入址| 一二三四中文在线观看免费高清| av在线播放精品| 99久久精品国产国产毛片| 99久久精品国产国产毛片| 男人爽女人下面视频在线观看| 亚洲国产欧美日韩在线播放| 午夜免费男女啪啪视频观看| 人妻一区二区av| 国产永久视频网站| 各种免费的搞黄视频| 搡老乐熟女国产| 一区二区三区乱码不卡18| 精品人妻熟女av久视频| 亚洲国产欧美日韩在线播放| 秋霞伦理黄片| 你懂的网址亚洲精品在线观看| 18在线观看网站| 男女边吃奶边做爰视频| 免费观看性生交大片5| av女优亚洲男人天堂| 伊人久久国产一区二区| 成人亚洲精品一区在线观看| 亚洲美女黄色视频免费看| 国产精品三级大全| 日韩成人av中文字幕在线观看| 男男h啪啪无遮挡| 在现免费观看毛片| 久久人人爽人人片av| 黄色欧美视频在线观看| www.色视频.com| 99re6热这里在线精品视频| 亚洲第一av免费看| 国产精品不卡视频一区二区| 成年人免费黄色播放视频| 赤兔流量卡办理| 大香蕉久久网| 国产精品国产三级专区第一集| 草草在线视频免费看| 色哟哟·www| 久久97久久精品| 久久99热6这里只有精品| 成人综合一区亚洲| 女人久久www免费人成看片| 男人操女人黄网站| 亚洲综合色惰| 成人18禁高潮啪啪吃奶动态图 | 久久久久精品性色| 国产欧美日韩一区二区三区在线 | 这个男人来自地球电影免费观看 | 69精品国产乱码久久久| 国产男人的电影天堂91| 综合色丁香网| 中文乱码字字幕精品一区二区三区| 美女视频免费永久观看网站| 久久久精品免费免费高清| 国产精品久久久久久久电影| 日韩人妻高清精品专区| 国产极品天堂在线| 国产伦精品一区二区三区视频9| 18禁在线播放成人免费| 99久久综合免费| 久久久精品区二区三区| 97在线人人人人妻| 欧美日本中文国产一区发布| 婷婷色麻豆天堂久久| 午夜激情久久久久久久| 午夜福利影视在线免费观看| 纯流量卡能插随身wifi吗| 欧美人与善性xxx| 日韩成人伦理影院| 中文天堂在线官网| 少妇的逼好多水| 免费观看无遮挡的男女| 啦啦啦啦在线视频资源| 精品久久久噜噜| 国产探花极品一区二区| 欧美日韩一区二区视频在线观看视频在线| 成人毛片a级毛片在线播放| 天堂中文最新版在线下载| 各种免费的搞黄视频| 亚洲欧美清纯卡通| 国产不卡av网站在线观看| 一区在线观看完整版| 在线观看免费高清a一片| 内地一区二区视频在线| 久久狼人影院| 国产成人精品久久久久久| 69精品国产乱码久久久| 午夜免费鲁丝| 最近的中文字幕免费完整| 免费人成在线观看视频色| 国产一区二区三区综合在线观看 | 精品熟女少妇av免费看| 久久久久视频综合| 国产精品熟女久久久久浪| 欧美一级a爱片免费观看看| 少妇的逼水好多| 国产精品 国内视频| 欧美少妇被猛烈插入视频| 国产精品一区www在线观看| 欧美变态另类bdsm刘玥| 只有这里有精品99| 久久精品久久精品一区二区三区| 啦啦啦啦在线视频资源| 欧美激情 高清一区二区三区| 亚洲精品自拍成人| a级毛片免费高清观看在线播放| 国产精品嫩草影院av在线观看| 蜜桃久久精品国产亚洲av| 国产av码专区亚洲av| 18禁裸乳无遮挡动漫免费视频| 卡戴珊不雅视频在线播放| 婷婷成人精品国产| 一个人看视频在线观看www免费| 22中文网久久字幕| 一级毛片 在线播放| 亚洲精品国产av成人精品| 黑人巨大精品欧美一区二区蜜桃 | 人人妻人人爽人人添夜夜欢视频| 亚洲精品aⅴ在线观看| 美女国产视频在线观看| 精品少妇久久久久久888优播| 亚洲精品日本国产第一区| 丝瓜视频免费看黄片| 夫妻性生交免费视频一级片| 欧美日韩精品成人综合77777| 伦理电影免费视频| 亚洲国产精品成人久久小说| 美女中出高潮动态图| 在线精品无人区一区二区三| 成人午夜精彩视频在线观看| 中文欧美无线码| 毛片一级片免费看久久久久| 丰满饥渴人妻一区二区三| 伦精品一区二区三区| 免费观看在线日韩| 亚洲综合色惰| 国产成人免费无遮挡视频| 久久久久久久大尺度免费视频| 精品亚洲成国产av| 伊人亚洲综合成人网| 国产午夜精品久久久久久一区二区三区| 日韩制服骚丝袜av| 日韩欧美精品免费久久| 超碰97精品在线观看| 亚洲精华国产精华液的使用体验| 日韩成人伦理影院| 亚洲精品456在线播放app| 精品亚洲成国产av| 人妻一区二区av| 亚洲经典国产精华液单| 午夜av观看不卡| 黄片播放在线免费| 中文字幕制服av| 亚洲国产精品999| 3wmmmm亚洲av在线观看| 亚洲欧美日韩卡通动漫| 亚洲美女搞黄在线观看| 色5月婷婷丁香| h视频一区二区三区| 高清黄色对白视频在线免费看| 成人国语在线视频| 中文乱码字字幕精品一区二区三区| 日日啪夜夜爽| 欧美日韩一区二区视频在线观看视频在线| 国产精品 国内视频| 极品人妻少妇av视频| 人妻少妇偷人精品九色| 免费人成在线观看视频色| 欧美日韩综合久久久久久| 久久久久久伊人网av| 午夜老司机福利剧场| av线在线观看网站| 九草在线视频观看| 一本—道久久a久久精品蜜桃钙片| 天美传媒精品一区二区| 日韩 亚洲 欧美在线| 黄色欧美视频在线观看| 免费观看在线日韩| 一区二区av电影网| xxxhd国产人妻xxx| 如何舔出高潮| 亚洲经典国产精华液单| 九九久久精品国产亚洲av麻豆| 国产成人a∨麻豆精品| 亚洲精品一区蜜桃| 国产高清不卡午夜福利| 五月伊人婷婷丁香| 美女内射精品一级片tv| 91成人精品电影| 高清黄色对白视频在线免费看| 精品久久蜜臀av无| 欧美精品亚洲一区二区| 国产午夜精品一二区理论片| 在线天堂最新版资源| 极品人妻少妇av视频| av电影中文网址| 成人综合一区亚洲| 日韩av在线免费看完整版不卡| 我的女老师完整版在线观看| 性色avwww在线观看| 国产乱人偷精品视频| 精品国产一区二区三区久久久樱花| 桃花免费在线播放| 人妻人人澡人人爽人人| 乱码一卡2卡4卡精品| 国产精品.久久久| 成人国语在线视频| 亚洲欧美日韩另类电影网站| 亚洲人成77777在线视频| av又黄又爽大尺度在线免费看| 插阴视频在线观看视频| 日本欧美视频一区| 亚洲精品aⅴ在线观看| 色视频在线一区二区三区| 丁香六月天网| 交换朋友夫妻互换小说| 激情五月婷婷亚洲| 91精品国产九色| 九草在线视频观看| av免费在线看不卡| 亚洲成人手机| 青春草视频在线免费观看| 国产精品成人在线| 热re99久久国产66热| 久久久精品免费免费高清| 91国产中文字幕| 国产精品久久久久久av不卡| 丝袜喷水一区| 免费人妻精品一区二区三区视频| 久久影院123| 日日啪夜夜爽| 男的添女的下面高潮视频| 美女xxoo啪啪120秒动态图| 如何舔出高潮| 精品国产一区二区三区久久久樱花| 日本与韩国留学比较| 大又大粗又爽又黄少妇毛片口| 黄色配什么色好看| 少妇人妻久久综合中文| 欧美日韩av久久| 一边摸一边做爽爽视频免费| 2021少妇久久久久久久久久久| 天堂俺去俺来也www色官网| 国产欧美亚洲国产| 亚洲情色 制服丝袜| 国产精品久久久久久精品古装| 亚洲精品美女久久av网站| 精品国产一区二区久久| a级毛片黄视频| 国产在线免费精品| 极品人妻少妇av视频| 建设人人有责人人尽责人人享有的| 久久久国产精品麻豆| 国产片内射在线| 国产深夜福利视频在线观看| 精品熟女少妇av免费看| 中国国产av一级| 午夜影院在线不卡| 大香蕉久久成人网| 国产免费一区二区三区四区乱码| 97在线人人人人妻| 日韩人妻高清精品专区| 高清av免费在线| 国产精品一二三区在线看| 国产精品一区二区在线观看99| 日韩制服骚丝袜av| 免费日韩欧美在线观看| 人人妻人人爽人人添夜夜欢视频| 国产一区二区在线观看日韩| 午夜福利视频精品| kizo精华| 久久久久久久久大av| av女优亚洲男人天堂| 亚洲人成网站在线观看播放| av福利片在线| 91精品伊人久久大香线蕉| 久久久久久久大尺度免费视频| 国产又色又爽无遮挡免| 91午夜精品亚洲一区二区三区| 国产伦精品一区二区三区视频9| 亚洲国产色片| 赤兔流量卡办理| 少妇丰满av| av电影中文网址| 久久亚洲国产成人精品v| 桃花免费在线播放| 亚洲精品456在线播放app| 久久久午夜欧美精品| av电影中文网址| 97超碰精品成人国产| 又粗又硬又长又爽又黄的视频| 国产精品一二三区在线看| 九九在线视频观看精品| 中国国产av一级| 久久国产精品男人的天堂亚洲 | kizo精华| 99久久中文字幕三级久久日本| 亚洲精品日本国产第一区| 亚洲av二区三区四区| 久久久国产一区二区| 自线自在国产av| 亚洲精品日韩在线中文字幕| 成年人午夜在线观看视频| 亚洲国产精品999| www.色视频.com| 婷婷色av中文字幕| 国产有黄有色有爽视频| 久久综合国产亚洲精品| 亚洲国产精品国产精品| 91午夜精品亚洲一区二区三区| 啦啦啦视频在线资源免费观看| 狠狠精品人妻久久久久久综合| 成人毛片60女人毛片免费| 一二三四中文在线观看免费高清| 亚洲精品第二区| 五月天丁香电影| 男人添女人高潮全过程视频| 亚洲第一区二区三区不卡| 搡女人真爽免费视频火全软件| 在现免费观看毛片| 成人综合一区亚洲| 婷婷色综合大香蕉| 在线 av 中文字幕| 99re6热这里在线精品视频| 精品午夜福利在线看| 亚洲精品日本国产第一区| 青青草视频在线视频观看| 91精品一卡2卡3卡4卡| 久久青草综合色| 夫妻性生交免费视频一级片| 国产片特级美女逼逼视频| 亚洲av欧美aⅴ国产| 国产精品国产三级国产专区5o| freevideosex欧美| 午夜精品国产一区二区电影| 亚洲精品乱久久久久久| 久久毛片免费看一区二区三区| 搡老乐熟女国产| 日韩,欧美,国产一区二区三区| 日日摸夜夜添夜夜爱| 欧美最新免费一区二区三区| 我要看黄色一级片免费的| 国产精品 国内视频| 另类精品久久| 亚洲伊人久久精品综合| 久久精品熟女亚洲av麻豆精品| 精品久久久久久久久av| 欧美变态另类bdsm刘玥| 亚洲欧美日韩卡通动漫| 99精国产麻豆久久婷婷| 精品人妻熟女av久视频| 日本与韩国留学比较| 亚洲国产欧美日韩在线播放| 亚洲av.av天堂| 久久国产精品男人的天堂亚洲 | 国产综合精华液| 久久久国产精品麻豆| 99久国产av精品国产电影| freevideosex欧美| 九九爱精品视频在线观看| 国产成人免费观看mmmm| 亚洲人成网站在线观看播放| 亚洲精品av麻豆狂野| 全区人妻精品视频| 狂野欧美白嫩少妇大欣赏| 哪个播放器可以免费观看大片| 自线自在国产av| 国产男人的电影天堂91| 考比视频在线观看| 日日摸夜夜添夜夜爱| 日本-黄色视频高清免费观看| 国产片特级美女逼逼视频| 日本与韩国留学比较| 在线亚洲精品国产二区图片欧美 | 亚洲精品久久成人aⅴ小说 | 国产色爽女视频免费观看| 亚洲精品色激情综合| 亚洲欧美成人精品一区二区| 丝袜美足系列| 日韩人妻高清精品专区| 999精品在线视频| 日韩大片免费观看网站| 91精品一卡2卡3卡4卡| 免费观看性生交大片5| 免费看光身美女| 久久久国产一区二区| 国产亚洲精品第一综合不卡 | 日韩一本色道免费dvd| 丝袜美足系列| www.av在线官网国产| 国产成人av激情在线播放 | 国产免费福利视频在线观看| 丰满少妇做爰视频| 精品国产露脸久久av麻豆| 亚洲欧美一区二区三区国产| 69精品国产乱码久久久| 久久国产精品男人的天堂亚洲 | 97精品久久久久久久久久精品| 午夜日本视频在线| 久久精品国产亚洲av天美| 中国三级夫妇交换| 国产伦理片在线播放av一区| 一级毛片我不卡| 2022亚洲国产成人精品| 日日啪夜夜爽| 国产精品一区二区在线不卡| 九草在线视频观看| 少妇猛男粗大的猛烈进出视频| 中文字幕久久专区| 欧美激情国产日韩精品一区| 亚洲怡红院男人天堂| 老熟女久久久| 国产欧美日韩综合在线一区二区| 在线免费观看不下载黄p国产| 国产免费一级a男人的天堂| 久久99热这里只频精品6学生| 尾随美女入室| 黄色配什么色好看| 亚洲经典国产精华液单| 国产亚洲av片在线观看秒播厂| 国产永久视频网站| 天堂俺去俺来也www色官网| 国产 一区精品| 亚洲欧美日韩另类电影网站| 亚洲精品乱码久久久久久按摩| 99久国产av精品国产电影| 欧美xxⅹ黑人| 丝袜脚勾引网站| 亚洲av中文av极速乱| 国产视频内射| 国产精品一国产av| 纯流量卡能插随身wifi吗| 搡女人真爽免费视频火全软件| 久久久久久久久久久免费av| 久久韩国三级中文字幕| 国产熟女欧美一区二区| 99热6这里只有精品| 99热网站在线观看| 亚洲人成网站在线观看播放| 成人无遮挡网站| 欧美精品一区二区大全| 99热这里只有精品一区| 国产男人的电影天堂91| 99国产精品免费福利视频| 在现免费观看毛片| 午夜av观看不卡| 中文字幕久久专区| 亚洲色图综合在线观看| 天天躁夜夜躁狠狠久久av| 免费黄网站久久成人精品| av又黄又爽大尺度在线免费看| 国产毛片在线视频| 国产精品熟女久久久久浪| 3wmmmm亚洲av在线观看| 男人添女人高潮全过程视频| 黄色毛片三级朝国网站| 国产永久视频网站| 十八禁高潮呻吟视频| 制服人妻中文乱码| 天堂中文最新版在线下载| 特大巨黑吊av在线直播| 欧美另类一区| 国产高清不卡午夜福利| 最近2019中文字幕mv第一页| 草草在线视频免费看| 中国美白少妇内射xxxbb| 最后的刺客免费高清国语| 亚洲av国产av综合av卡| 又黄又爽又刺激的免费视频.| 亚洲国产精品一区三区| 制服丝袜香蕉在线| 亚洲国产欧美日韩在线播放| 夫妻午夜视频| 国产av码专区亚洲av| 91精品国产九色| 国产成人精品婷婷| 亚洲成色77777| 天堂俺去俺来也www色官网| 国产成人精品婷婷| 久久女婷五月综合色啪小说| 国产日韩欧美视频二区| 在线观看三级黄色| 亚洲经典国产精华液单| 在线免费观看不下载黄p国产| 菩萨蛮人人尽说江南好唐韦庄| 精品人妻在线不人妻| 国产在视频线精品| 亚洲天堂av无毛| 亚洲人成网站在线观看播放| 春色校园在线视频观看|