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

    考慮腐蝕磨損的全生命期內(nèi)錨鏈極限強(qiáng)度和疲勞演變分析

    2023-09-22 01:48:50張科棟王文華王巍巍
    船舶力學(xué) 2023年9期
    關(guān)鍵詞:鏈環(huán)錨鏈系泊

    張科棟,王文華,楊 光,王巍巍,李 瑤,劉 剛,黃 一

    (1.大連理工大學(xué)船舶工程學(xué)院,遼寧 大連 116024;2.海洋石油工程股份有限公司,天津 300461)

    0 引 言

    系泊錨鏈長(zhǎng)期服役于海水中,海洋環(huán)境條件以及循環(huán)交變載荷的作用會(huì)使錨鏈結(jié)構(gòu)發(fā)生磨損、老化、銹蝕等現(xiàn)象,從而引起系泊構(gòu)件強(qiáng)度降低,導(dǎo)致錨鏈?zhǔn)茢嗟氖鹿蕰r(shí)有發(fā)生[1]。其中,疲勞斷裂和極限破斷是錨鏈常見的失效模式。在海洋環(huán)境條件影響下平臺(tái)產(chǎn)生往復(fù)運(yùn)動(dòng),錨鏈?zhǔn)艿窖h(huán)交變載荷作用,在關(guān)鍵位置存在應(yīng)力集中,容易誘發(fā)疲勞破壞。同時(shí),在長(zhǎng)時(shí)間的服役過程中,遭遇極端海況時(shí),錨鏈會(huì)受較大的張力載荷,極易達(dá)到強(qiáng)度極限從而發(fā)生破斷。而在系泊錨鏈服役期間腐蝕磨損的影響會(huì)降低錨鏈的極限強(qiáng)度、增加疲勞損傷,因此準(zhǔn)確分析錨鏈在全生命期內(nèi)疲勞損失和極限強(qiáng)度對(duì)于保障系泊系統(tǒng)可靠性是至關(guān)重要的。

    針對(duì)上述問題,劉金沅等[2]對(duì)系泊錨鏈的腐蝕進(jìn)行了實(shí)驗(yàn)研究,并與數(shù)值模擬方法一起對(duì)錨鏈腐蝕損傷進(jìn)行了對(duì)比分析;徐偉等[3]基于接觸算法和相對(duì)強(qiáng)度法,提出了浮式平臺(tái)系泊系統(tǒng)損傷錨鏈強(qiáng)度評(píng)估方法,結(jié)果表明隨著腐蝕程度增加錨鏈強(qiáng)度明顯減弱;Steenkiste 等[4]對(duì)錨鏈的磨損機(jī)理和磨損形式進(jìn)行了研究,將磨損分為磨粒磨損、雙體磨損、三體磨損以及腐蝕磨損等;喬?hào)|生等[5]研究了損傷錨鏈在不同平面外彎曲角度下結(jié)構(gòu)應(yīng)力變化與接觸特性;Lassen 等[6]提出了采用熱點(diǎn)應(yīng)力法,結(jié)合SN曲線計(jì)算錨鏈面外彎曲(out-of-plane bending,OPB)疲勞損傷;Kim等[7]對(duì)FPSO系泊系統(tǒng)的整體和局部進(jìn)行了疲勞分析,考慮OPB現(xiàn)象對(duì)系泊系統(tǒng)的影響。到目前為止,同時(shí)考慮錨鏈腐蝕和磨損狀態(tài)的研究相對(duì)較少,此外有關(guān)在全生命期內(nèi)錨鏈損傷狀態(tài)下極限強(qiáng)度和疲勞損傷演變的研究極少。

    據(jù)此,本文基于Archard 磨損理論和均勻腐蝕模型,采用有限元數(shù)值仿真方法計(jì)算錨鏈的累積損傷,建立系泊錨鏈的全生命期損傷演變模型?;阱^鏈損傷模型,提出系泊錨鏈在全生命期內(nèi)考慮腐蝕磨損的極限強(qiáng)度演變和疲勞損傷演變?cè)u(píng)估方法,其中采用逐步加載方式,以材料第四強(qiáng)度理論為判斷標(biāo)準(zhǔn),對(duì)系泊錨鏈極限強(qiáng)度演變進(jìn)行分析;基于熱點(diǎn)應(yīng)力集中系數(shù),利用雨流計(jì)數(shù)法、S-N曲線法和Miner疲勞累積損傷理論對(duì)錨鏈在長(zhǎng)期海況下的疲勞損傷演變進(jìn)行分析。

    1 算法理論簡(jiǎn)介

    1.1 系泊錨鏈腐蝕和磨損損傷模型

    由于系泊錨鏈長(zhǎng)期服役于海水中,常常面臨著嚴(yán)重的腐蝕作用,導(dǎo)致錨鏈的直徑變細(xì);同時(shí)由于平臺(tái)運(yùn)動(dòng)和海流的影響,系泊錨鏈長(zhǎng)期受到張力載荷作用,鏈環(huán)接觸面之間有較大的接觸壓力,而且接觸面間存在著微小的滑動(dòng),進(jìn)而導(dǎo)致磨損的存在,使鏈環(huán)的接觸面形貌發(fā)生改變。

    關(guān)于錨鏈腐蝕速率的確定,DNV-OS-E301[8]在設(shè)計(jì)系泊錨鏈時(shí),在錨鏈直徑方向上考慮了0.4 mm/y的腐蝕余量;根據(jù)Gao等[9]的研究結(jié)果,不同時(shí)間的平均腐蝕速率如表1所示。

    表1 錨鏈不同階段的腐蝕速率Tab.1 Corrosion rate of anchor chain at different stages

    但是在實(shí)際的工程應(yīng)用當(dāng)中,錨鏈局部通常會(huì)施加一些陰極保護(hù)措施,減緩腐蝕對(duì)錨鏈性能的影響。因此錨鏈的腐蝕速率會(huì)因?yàn)樗幍暮S驐l件和采取的防腐措施的不同,產(chǎn)生巨大的差異,會(huì)對(duì)錨鏈的損傷演變分析產(chǎn)生巨大的影響。因此本文將采取腐蝕速率為0~0.4 mm/y 區(qū)間定義腐蝕損傷,來反映各種條件下的錨鏈腐蝕狀態(tài),并對(duì)錨鏈損傷演變進(jìn)行分析。

    關(guān)于錨鏈磨損量的計(jì)算是基于Archard 磨損方程[10],如式(1)所示。根據(jù)有限元計(jì)算得到的接觸面法向壓力,可推導(dǎo)得到接觸區(qū)域的磨損體積。

    式中,V為磨損體積,K為無量綱的磨損系數(shù),P為接觸面的法向壓力,L為滑移距離,H為材料接觸表面的硬度。K表示磨損過程中粗糙表面相互作用形成磨損顆粒的概率,可以用k來代替K/H,表示有量綱的磨損系數(shù),根據(jù)實(shí)驗(yàn)結(jié)果[11]取k=4×10-10MPa-1;L表示鏈環(huán)間滑移距離的局部增量,可以由鏈環(huán)之間轉(zhuǎn)動(dòng)的位移關(guān)系得到。

    對(duì)于鏈環(huán)之間的磨損,本文采用Archard 磨損方程,根據(jù)ANSYS 計(jì)算得到接觸面的壓應(yīng)力和錨鏈運(yùn)動(dòng)的轉(zhuǎn)動(dòng)位移,求得鏈環(huán)之間的磨損體積。在鏈環(huán)模型進(jìn)行腐蝕修正之后,以鏈環(huán)的磨損體積來確定幾何模型布爾減運(yùn)算的切割深度,最終得到鏈環(huán)的損傷模型。因此在每次循環(huán)過程中,都要在計(jì)算完成后提取接觸面的壓應(yīng)力,計(jì)算磨損體積,并根據(jù)腐蝕后的鏈環(huán)直徑與切削體積的函數(shù)關(guān)系,得到累加的磨損深度,在下次循環(huán)開始時(shí)對(duì)模型進(jìn)行修正。本文的系泊錨鏈損傷分析流程,如圖1所示。

    圖1 系泊錨鏈損傷分析流程圖Fig.1 Damage analysis flow chart of mooring anchor chain

    1.2 錨鏈模型極限強(qiáng)度演變的分析方法

    在船舶與海洋工程領(lǐng)域,極限強(qiáng)度研究方法主要包括有限元法、經(jīng)驗(yàn)公式和解析法、模型試驗(yàn)法等[12]。經(jīng)驗(yàn)公式和解析法使用簡(jiǎn)便,但是使用對(duì)象和場(chǎng)景有限,由于錨鏈存在接觸關(guān)系和結(jié)構(gòu)的特殊性,難以找到合適的極限強(qiáng)度經(jīng)驗(yàn)公式;模型試驗(yàn)法最接近實(shí)際情況,但是難以對(duì)全生命期內(nèi)錨鏈累積損傷下的極限強(qiáng)度進(jìn)行實(shí)驗(yàn);而有限元法可以通過建立腐蝕磨損的損傷模型,對(duì)錨鏈進(jìn)行極限強(qiáng)度分析,因此本文采用非線性有限元法進(jìn)行系泊錨鏈的極限強(qiáng)度計(jì)算。

    本文基于ANSYS 有限元軟件計(jì)算錨鏈的極限張力,根據(jù)系泊錨鏈在真實(shí)海洋環(huán)境中的邊界條件,考慮彈塑性材料本構(gòu)關(guān)系以及接觸非線性,在錨鏈腐蝕磨損結(jié)構(gòu)損傷模型上采用逐步加載的方式,當(dāng)最大等效應(yīng)力值達(dá)到極限強(qiáng)度時(shí),即得到對(duì)應(yīng)的極限張力載荷。根據(jù)不同的錨鏈損傷模型,經(jīng)過計(jì)算最終得到不同腐蝕磨損損傷下的系泊錨鏈極限張力演變規(guī)律。

    1.3 錨鏈模型疲勞損傷演變的分析方法

    海洋平臺(tái)在風(fēng)浪流作用下,平臺(tái)運(yùn)動(dòng)呈現(xiàn)周期性的波動(dòng),導(dǎo)致系泊錨鏈承受著循環(huán)交變載荷的作用,極易引起錨鏈局部的疲勞破壞。針對(duì)這一問題,本文對(duì)全生命期內(nèi),考慮腐蝕磨損損傷的錨鏈進(jìn)行疲勞損傷評(píng)估,對(duì)錨鏈的熱點(diǎn)應(yīng)力進(jìn)行計(jì)算,提出考慮腐蝕磨損損傷的系泊錨鏈OPB疲勞損傷評(píng)估方法,其分析流程如圖2所示。

    圖2 系泊錨鏈疲勞損傷分析流程圖Fig.2 Flow chart of fatigue damage analysis of mooring anchor chain

    為了計(jì)算得到錨鏈在實(shí)際海況下受到的載荷,需要對(duì)浮體進(jìn)行動(dòng)力響應(yīng)分析。首先對(duì)海洋結(jié)構(gòu)物所處海域的波浪散布圖進(jìn)行簡(jiǎn)化處理[12],然后采用AWAQ軟件計(jì)算錨鏈動(dòng)力響應(yīng),得到系泊錨鏈關(guān)鍵位置的張力時(shí)程曲線和轉(zhuǎn)角時(shí)程曲線,將其作為錨鏈有限元計(jì)算施加的載荷。其次,采用ANSYS軟件進(jìn)行全生命期內(nèi)錨鏈考慮腐蝕磨損的累積損傷模擬。最后,根據(jù)不同損傷模型,采用時(shí)域分析法對(duì)錨鏈OPB疲勞損傷進(jìn)行計(jì)算,然后對(duì)錨鏈全生命期內(nèi)的疲勞損傷進(jìn)行評(píng)估。

    根據(jù)相關(guān)研究,系泊錨鏈最上端的鏈環(huán)最易發(fā)生疲勞破壞,鏈環(huán)之間受到張力鎖緊,接觸面存在著較大的摩擦力,當(dāng)鏈環(huán)間發(fā)生角度偏轉(zhuǎn)時(shí),便發(fā)生OPB 現(xiàn)象[13]。OPB 疲勞熱點(diǎn)位于鏈環(huán)彎曲段接近于接觸區(qū)域的表面,規(guī)范提供了兩個(gè)典型疲勞熱點(diǎn)位置,見圖3中A、B兩點(diǎn),A點(diǎn)為:αA=17°,βA=28°;B點(diǎn)為:αB=36°,βB=36°。錨鏈在純拉伸張力載荷下,疲勞熱點(diǎn)位于鏈環(huán)直線段到彎曲段內(nèi)側(cè)過渡處的表面,如圖3中C點(diǎn)所示。

    圖3 錨鏈疲勞熱點(diǎn)位置Fig.3 Hotspot locations of mooring chains

    錨鏈OPB 疲勞是由于純張力和OPB 彎矩聯(lián)合作用的影響,BV 規(guī)范[14]指出將熱點(diǎn)應(yīng)力σH分解為由張力作用得到的應(yīng)力分量σT-T和由OPB 彎矩作用得到的應(yīng)力分量σOPB:

    式中,σH為錨鏈的熱點(diǎn)應(yīng)力,σn,T-T為錨鏈張力作用下的名義應(yīng)力,σn,OPB為錨鏈的OPB彎矩作用下的名義應(yīng)力,T為錨鏈?zhǔn)艿降膹埩χ担琈為相鄰鏈環(huán)接觸區(qū)域的彎矩,d為錨鏈直徑,ST-T為張力作用的熱點(diǎn)應(yīng)力集中系數(shù),SOPB為彎矩作用的熱點(diǎn)應(yīng)力集中系數(shù)。

    對(duì)于疲勞損傷的計(jì)算,參考BV規(guī)范[14],取S-N曲線為

    式中,nc(S)為應(yīng)力范圍循環(huán)次數(shù);S為第一主應(yīng)力范圍,單位為MPa;lg(aD)為S-N曲線截距,取值12.436;m為S-N曲線的斜率,取值3.0。

    2 系泊錨鏈腐蝕和磨損損傷模型

    2.1 錨鏈模型參數(shù)

    本文選取的系泊錨鏈為無檔形式,如圖4所示,材質(zhì)為R3S級(jí)錨鏈鋼,無檔錨鏈直徑為158 mm,根據(jù)ISO 1704規(guī)范確定錨鏈的幾何尺寸,鏈環(huán)的尺寸及材料等主要參數(shù)見表2。

    圖4 無檔錨鏈Fig.4 Studless chain

    表2 無檔錨鏈尺寸和材料參數(shù)Tab.2 Studless chain standard sizes and material parameters

    本文的幾何模型為由三個(gè)鏈環(huán)組成的錨鏈段,以中間完整鏈環(huán)為研究對(duì)象,在兩端鏈環(huán)橫截面分別施加載荷和約束,并對(duì)錨鏈的接觸區(qū)域進(jìn)行細(xì)化切分,通過設(shè)置接觸對(duì)實(shí)現(xiàn)相鄰鏈環(huán)的接觸關(guān)系。模型采用SOLID95單元,接觸對(duì)采用CONTA174和TARGE170單元,在橫截面處建立MPC184單元,在其中心節(jié)點(diǎn)處施加載荷,其有限元模型如圖5所示。

    圖5 錨鏈有限元模型Fig.5 Finite element model of mooring chain

    2.2 錨鏈張力和直徑對(duì)錨鏈損傷的影響

    在建立全生命周期內(nèi)錨鏈損傷模型時(shí),最重要的問題是磨損與腐蝕二者對(duì)于錨鏈的影響。首先腐蝕會(huì)由于是否施加陰極保護(hù)等防腐措施以及工作海域等因素,導(dǎo)致腐蝕速率變化不一,而磨損則是一個(gè)隨服役年限的增加而不斷疊加的過程,而且二者之間的相互影響尚不明確。因此在整個(gè)生命周期內(nèi),建立一個(gè)有效準(zhǔn)確的損傷模型,有利于對(duì)錨鏈性能的準(zhǔn)確評(píng)估。

    把某海況下頂端張力平均值1718.56 kN、鏈間轉(zhuǎn)角幅值0.5°、假定的循環(huán)載荷周期24 s,作為有限元的施加載荷進(jìn)行計(jì)算。對(duì)于錨鏈損傷模型的建立,首先考慮載荷對(duì)磨損損傷的影響進(jìn)行了分析,根據(jù)Archard 磨損方程可以明顯地看到磨損體積與滑移距離成正比關(guān)系,因此在初始158 mm 的模型中,滑移距離相同的情況下,將頂端張力上下變化20%之后,計(jì)算得到錨鏈磨損體積數(shù)值,如表3 所示,由此可以得到,磨損體積隨張力近似呈正比關(guān)系。這也說明了隨著施加張力載荷的變化,接觸面的法向應(yīng)力與載荷近似成正比,根據(jù)Archard磨損方程,磨損體積也呈正相關(guān)的關(guān)系。

    表3 張力對(duì)磨損的影響Tab.3 Effect of tension on wear

    考慮到錨鏈的直徑因?yàn)楦g發(fā)生變化對(duì)磨損的影響,采取在相同的載荷條件下,計(jì)算錨鏈初始直徑為158 mm 和經(jīng)過若干年腐蝕后的直徑為142 mm 時(shí),磨損體積的變化情況。根據(jù)計(jì)算結(jié)果可以看到,錨鏈的直徑對(duì)于初始磨損的影響很小,二者磨損體積的誤差很小,如表4所示。

    表4 錨鏈直徑對(duì)初始磨損的影響Tab.4 Effect of anchor chain diameter on initial wear

    在進(jìn)行錨鏈損傷計(jì)算時(shí),要考慮累積磨損錨鏈的磨損體積變化情況,因此計(jì)算了錨鏈直徑為158 mm 和142 mm 時(shí),在相同磨損體積的前提下,施加相同的載荷,進(jìn)行下一次的磨損體積計(jì)算,其結(jié)果如表5 所示。根據(jù)計(jì)算結(jié)果可以得到,在相同磨損體積和相同載荷的條件下,不同的腐蝕直徑對(duì)于磨損的影響很小。

    表5 錨鏈直徑對(duì)累加磨損的影響Tab.5 Effect of anchor chain diameter on cumulative wear

    進(jìn)而根據(jù)之前的分析,可以看出腐蝕直徑對(duì)于磨損的影響很小,磨損只會(huì)隨著服役年限的增加而不斷累積。產(chǎn)生這樣現(xiàn)象的主要原因是,由于施加的載荷相同,在直徑較大時(shí),鏈環(huán)之間的接觸面積較大,而接觸應(yīng)力相對(duì)較小;在直徑因腐蝕變小時(shí),鏈環(huán)之間的接觸面積較小,接觸應(yīng)力值變大,但是二者的面積與應(yīng)力乘積的垂向載荷值變化不大,因此根據(jù)Archard磨損方程計(jì)算得到的磨損體積值變化不大。

    根據(jù)之前的分析,可以看出載荷對(duì)于磨損的影響幾乎呈正比的關(guān)系,因此可以根據(jù)南海實(shí)際海浪環(huán)境,首先利用簡(jiǎn)化波浪散布圖[12],然后采用AQWA軟件進(jìn)行FPSO運(yùn)動(dòng)響應(yīng)計(jì)算,可以得到短期海況下所有系泊鏈頂端張力和轉(zhuǎn)角時(shí)程,圖6為海況1的張力和轉(zhuǎn)角時(shí)程。

    圖6 海況1條件下錨鏈張力和轉(zhuǎn)角時(shí)程Fig.6 Anchor chain tension and angle time course under Sea State Condition 1

    考慮計(jì)算成本,假定系泊錨鏈在服役期間內(nèi)面臨的海況浪向的來臨是均勻分布的,即按照0°、60°、120°、180°、240°、300°進(jìn)行計(jì)算,張力取平均值,轉(zhuǎn)角產(chǎn)生的位移為轉(zhuǎn)角時(shí)程的波動(dòng)之和,然后按照海況的概率進(jìn)行求和,最終得到系泊鏈頂端錨鏈的載荷,如表6所示。

    表6 錨鏈張力和轉(zhuǎn)角Tab.6 Anchor chain tension and angle

    為了減小計(jì)算量,按照錨鏈生命期腐蝕速率進(jìn)行分段,分別計(jì)算得到0、4、8、12、16、20、25、30 年的損傷模型。從圖7 中可以看出,鏈環(huán)的磨損體積在整個(gè)服役周期內(nèi)變化較小,因此本文采取周期內(nèi)所有階段磨損體積的平均值3828.28 mm3,作為年磨損體積進(jìn)行修正損傷模型。因此,錨鏈的典型腐蝕和磨損損傷幾何模型如圖8所示。

    圖7 錨鏈磨損體積的變化趨勢(shì)Fig.7 Trend of anchor chain wear volume

    圖8 錨鏈損傷幾何模型Fig.8 Damage geometry model of anchor chain

    根據(jù)以上的分析,錨鏈直徑對(duì)磨損體積的影響不大,磨損體積可以簡(jiǎn)化處理為只與時(shí)間相關(guān)的物理量;錨鏈均勻腐蝕為與直徑相關(guān)的物理量,二者之間的相互影響很小,幾乎可以忽略不計(jì),由此建立錨鏈完整的生命周期損傷模型,并用于后面的極限強(qiáng)度演變和疲勞損傷演變分析。

    3 錨鏈極限強(qiáng)度和疲勞損傷演變特性

    3.1 錨鏈極限強(qiáng)度演變

    根據(jù)錨鏈損傷模型,采用非線性有限元分析方法,以第四強(qiáng)度理論為判斷標(biāo)準(zhǔn)計(jì)算錨鏈的極限破斷張力。在張力載荷不斷增大時(shí),錨鏈在載荷作用下嚴(yán)重變形,從彈性階段過渡到塑性階段,兩端被拉長(zhǎng),接觸面積變大,直線段中間向內(nèi)側(cè)彎曲。從圖9 的應(yīng)力分布圖中可以看到,直線段和彎曲段連接位置和冠部位置的應(yīng)力水平最高,是易發(fā)生斷裂失效的區(qū)域,因此當(dāng)直線段和彎曲段連接位置的Von Mises 等效應(yīng)力大于材料給定的極限強(qiáng)度時(shí),認(rèn)為錨鏈發(fā)生斷裂破壞,此時(shí)對(duì)應(yīng)的張力載荷即為錨鏈的極限破斷張力。

    圖9 錨鏈磨損損傷等效應(yīng)力云圖Fig.9 Von Mises stress distribution of wear anchor chain

    在施加相同的張力載荷16 000 kN 時(shí),只考慮磨損損傷的條件下,計(jì)算了相同直徑158 mm 錨鏈,生命周期時(shí)間分別為0 年和16 年?duì)顟B(tài)下的應(yīng)力分布。從圖9 的Von Mises 等效應(yīng)力云圖中可以看出,磨損后應(yīng)力集中位置同樣出現(xiàn)在彎曲段與接觸面之間的區(qū)域,但是隨著磨損程度的增加,最大應(yīng)力點(diǎn)的位置逐漸向磨損區(qū)域的邊緣靠近,接觸面附近的應(yīng)力分布有所改變,二者的應(yīng)力值變化較小。說明磨損損傷的存在會(huì)逐漸改變局部的應(yīng)力分布,是錨鏈性能分析中值得關(guān)注的問題。

    在考慮腐蝕的影響時(shí),計(jì)算了相同磨損狀態(tài)即生命周期為30年時(shí),錨鏈直徑為158 mm和144 mm時(shí)的應(yīng)力分布狀態(tài),如圖10所示。根據(jù)計(jì)算結(jié)果可知,由于腐蝕導(dǎo)致錨鏈直徑變小,在應(yīng)力云圖中可以看出腐蝕模型的應(yīng)力數(shù)值明顯增大,但是二者的應(yīng)力分布近似,都在鏈環(huán)直線段與彎曲段過渡位置存在應(yīng)力集中的現(xiàn)象。

    圖10 錨鏈腐蝕損傷等效應(yīng)力云圖Fig.10 Von Mises stress distribution of corrosion anchor chain

    根據(jù)錨鏈損傷模型,采用逐步加載的方式計(jì)算極限張力,當(dāng)最大等效應(yīng)力值達(dá)到極限強(qiáng)度時(shí)得到極限張力,采用衰減系數(shù)表示生命周期內(nèi)錨鏈極限破斷張力的演變,其結(jié)果如表7所示。

    表7 錨鏈極限破斷張力衰減系數(shù)Tab.7 Ultimate breaking tension decay factor of anchor chain

    以最大腐蝕速率為例,在腐蝕磨損的影響下,錨鏈從服役初期可以承受16 119 kN 的張力載荷減小到服役30 年后只能承受不到11 081 kN 的張力,衰減率達(dá)到了31.3%,表明錨鏈腐蝕磨損損傷對(duì)于極限強(qiáng)度的影響十分關(guān)鍵。

    從表7 中的計(jì)算結(jié)果可知,磨損對(duì)于極限破斷張力的影響較小,但是也不容忽視,生命周期內(nèi)大概有接近4%的衰減,但是腐蝕速率的快慢則會(huì)對(duì)錨鏈的極限破斷張力數(shù)值產(chǎn)生巨大的影響。

    圖11(a)為在無腐蝕損失狀態(tài)下即錨鏈直徑為158 mm時(shí),錨鏈?zhǔn)艿侥p的影響,極限強(qiáng)度隨磨損損傷的演變;圖11(b)為錨鏈磨損為30年時(shí),腐蝕速率為0.4 mm/y的條件下,即錨鏈直徑初始時(shí)為158 mm,30 年時(shí)為134 mm,極限強(qiáng)度隨腐蝕損傷的演變。從曲線圖中可以看出,腐蝕對(duì)錨鏈的極限破斷張力影響遠(yuǎn)遠(yuǎn)大于磨損,而且錨鏈腐蝕后極限破斷張力與直徑呈近似線性波動(dòng)快速下降的趨勢(shì),在整個(gè)生命周期內(nèi)下降了大約三分之一左右,因此腐蝕磨損損傷對(duì)于錨鏈極限強(qiáng)度的影響是不容忽視的。

    圖11 錨鏈損傷對(duì)極限強(qiáng)度演變的影響Fig.11 Effect of anchor chain damage on the evolution of ultimate strength

    3.2 疲勞損傷演變?cè)u(píng)估

    根據(jù)錨鏈損傷模型,采用時(shí)域分析法進(jìn)行OPB 疲勞損傷評(píng)估,其中,OPB 的應(yīng)力值由轉(zhuǎn)角時(shí)程計(jì)算得到,而張力作用的熱點(diǎn)應(yīng)力時(shí)程經(jīng)過有限元計(jì)算得到熱點(diǎn)應(yīng)力集中系數(shù)而求得。然后采用雨流計(jì)數(shù)法和S-N曲線計(jì)算得到錨鏈的疲勞壽命,然后采用Miner 線性累積損傷準(zhǔn)則計(jì)算錨鏈的疲勞損傷,最終完成對(duì)錨鏈在全生命期內(nèi)的疲勞損傷評(píng)估。

    為了探究錨鏈磨損對(duì)疲勞損傷的影響,在只考慮磨損損傷的條件下,計(jì)算了相同直徑(158 mm)錨鏈的熱點(diǎn)應(yīng)力隨服役時(shí)間的變化,圖12顯示的是生命周期時(shí)間分別為0年和16年?duì)顟B(tài)下的應(yīng)力分布。

    圖12 錨鏈磨損損傷第一主應(yīng)力云圖Fig.12 First principal stress distribution of wear anchor chain

    從圖12的應(yīng)力云圖中可看到,隨著磨損程度的增加,磨損接觸面不斷加深,導(dǎo)致接觸面積不斷增大,最大應(yīng)力點(diǎn)位置逐漸向外側(cè)移動(dòng),與接觸面的距離更近,說明接觸面磨損形貌的改變影響了錨鏈局部的應(yīng)力大小和分布,其原因是外力只影響作用區(qū)域的應(yīng)力分布,在離外力足夠遠(yuǎn)的區(qū)域,由于外力變化而產(chǎn)生的應(yīng)力變化將忽略不計(jì)。在本模型中,接觸面即可視為施加載荷位置,當(dāng)接觸面發(fā)生磨損后,接觸面形貌發(fā)生變化,而熱點(diǎn)區(qū)域應(yīng)力值隨之發(fā)生改變,說明其仍處于載荷作用區(qū)域內(nèi),所以磨損對(duì)于錨鏈疲勞熱點(diǎn)處的應(yīng)力值會(huì)產(chǎn)生影響,進(jìn)行錨鏈疲勞分析時(shí)對(duì)于磨損損傷不可忽視。同時(shí)從圖13 中可以看到,隨著磨損損傷的增加,熱點(diǎn)處的應(yīng)力值呈現(xiàn)波動(dòng)的現(xiàn)象,這是因?yàn)殒湱h(huán)間磨損的存在,改變了錨鏈局部的應(yīng)力分布,導(dǎo)致熱點(diǎn)處應(yīng)力的波動(dòng)。

    圖13 熱點(diǎn)應(yīng)力隨磨損損傷的變化Fig.13 Variation of hot spot stress with wear damage

    為了探究錨鏈腐蝕對(duì)疲勞損傷的影響,在只考慮腐蝕損傷的條件下,計(jì)算了生命周期為30 年時(shí)相同磨損狀態(tài)下,不同錨鏈直徑的熱點(diǎn)應(yīng)力值,圖14顯示的是錨鏈直徑為158 mm和144 mm時(shí)的應(yīng)力分布狀態(tài)。

    圖14 錨鏈腐蝕損傷第一主應(yīng)力云圖Fig.14 First principal stress distribution of corrosion anchor chain

    從圖14中可以看出,在同一磨損條件下,腐蝕造成的錨鏈直徑的變化,對(duì)于應(yīng)力的分布沒有太大的影響,只是由于直徑導(dǎo)致錨鏈橫截面積減小,錨鏈整體的應(yīng)力值有所上升,從而導(dǎo)致了熱點(diǎn)處應(yīng)力逐漸增加的趨勢(shì),如圖15所示。根據(jù)名義應(yīng)力計(jì)算公式可知,名義應(yīng)力與直徑呈二次反比的關(guān)系,隨著直徑的減小,名義應(yīng)力會(huì)迅速增加。但是有限元計(jì)算的熱點(diǎn)應(yīng)力值會(huì)受到實(shí)際損傷模型的影響,名義應(yīng)力與熱點(diǎn)應(yīng)力增長(zhǎng)速度不一定一致,因此熱點(diǎn)處的應(yīng)力集中系數(shù)存在減小的可能。錨鏈結(jié)構(gòu)熱點(diǎn)應(yīng)力集中系數(shù)如表8所示。

    圖15 熱點(diǎn)應(yīng)力隨腐蝕直徑的變化Fig.15 Variation of hot spot stress with corrosion diameter

    表8 生命周期內(nèi)錨鏈結(jié)構(gòu)熱點(diǎn)應(yīng)力集中系數(shù)Tab.8 SCF of anchor chain structure during life cycle

    根據(jù)公式(2)~(4),首先進(jìn)行錨鏈OPB彎矩計(jì)算,之后利用計(jì)算得到熱點(diǎn)的應(yīng)力集中系數(shù),根據(jù)上述計(jì)算張力時(shí)程曲線和轉(zhuǎn)角時(shí)程曲線,得到熱點(diǎn)應(yīng)力時(shí)程曲線。根據(jù)得到的熱點(diǎn)處應(yīng)力時(shí)程曲線,按照?qǐng)D2的方法,計(jì)算得到實(shí)際海況下系泊錨鏈結(jié)構(gòu)熱點(diǎn)位置的疲勞損傷,然后根據(jù)海況分布的概率進(jìn)行平均和累計(jì),最終得到錨鏈在該損傷模型下的年疲勞損傷值,這里以無腐蝕磨損損傷錨鏈為例,具體數(shù)據(jù)見表9。

    表9 無腐蝕磨損損傷錨鏈的年疲勞損傷Tab.9 Annual fatigue damage of anchor chain without corrosion and wear damage

    根據(jù)表8 數(shù)據(jù)計(jì)算得到錨鏈疲勞損傷,分析在全生命期內(nèi)磨損和腐蝕對(duì)錨鏈疲勞損傷的影響,如圖16 所示。在圖16(a)中,在無腐蝕損失狀態(tài)下錨鏈直徑為158 mm,隨著服役時(shí)間的增加,錨鏈磨損損傷加劇,但是錨鏈結(jié)構(gòu)熱點(diǎn)處的疲勞損傷波動(dòng)變化較小,只有10%左右;在圖16(b)中,在磨損狀態(tài)為30 年時(shí),錨鏈直徑以0.4 mm/y 的腐蝕速率從158 mm 到134 mm 變化,其疲勞損傷呈現(xiàn)非線性增長(zhǎng)的趨勢(shì),增長(zhǎng)速度逐漸加快。這表明腐蝕磨損導(dǎo)致的鏈環(huán)形貌發(fā)生損傷,進(jìn)而使得在全壽命期內(nèi)134 mm 錨鏈相比158 mm 錨鏈的疲勞損傷增加了175%。因此在進(jìn)行錨鏈疲勞分析時(shí)不能忽視腐蝕磨損的損傷對(duì)錨鏈帶來的影響。

    圖16 錨鏈損傷對(duì)疲勞演變的影響Fig.16 Effect of anchor chain damage on fatigue evolution

    4 結(jié) 論

    基于錨鏈損傷演變數(shù)值模擬方法,本文提出了考慮腐蝕磨損條件下錨鏈在全生命期內(nèi)極限強(qiáng)度演變和疲勞損傷演變的評(píng)估方法,具體結(jié)論如下:

    (1)針對(duì)南海典型海況,根據(jù)錨鏈損傷演變數(shù)值模擬方法,定量研究了典型工況錨鏈的腐蝕量和磨損量。結(jié)果表明,錨鏈直徑對(duì)磨損體積的影響不大,磨損體積可以簡(jiǎn)化處理為只與時(shí)間相關(guān)的物理量,而錨鏈均勻腐蝕表示為與直徑相關(guān)的物理量。

    (2)基于全生命期內(nèi)的錨鏈損傷模型進(jìn)行極限破斷張力演變和疲勞損傷演變分析。結(jié)果表明,隨著磨損的產(chǎn)生,接觸面形貌發(fā)生變化,磨損影響了錨鏈局部的應(yīng)力大小和分布;在腐蝕的影響下,錨鏈的應(yīng)力分布沒有太大的改變,但是應(yīng)力的數(shù)值會(huì)產(chǎn)生變化,導(dǎo)致極限破斷張力和熱點(diǎn)應(yīng)力數(shù)值會(huì)隨之改變。

    (3)在極限破斷張力演變分析中,考慮腐蝕磨損影響下錨鏈的極限破斷張力逐漸減小。在磨損的影響下,全生命周期內(nèi)錨鏈極限破斷張力減小4%左右;在腐蝕的影響下,當(dāng)腐蝕后直徑減小到一定程度時(shí),破斷張力會(huì)迅速減小,生命周期內(nèi)錨鏈極限破斷張力減小了31.3%。腐蝕和磨損對(duì)錨鏈極限強(qiáng)度影響顯著,在錨鏈極限強(qiáng)度演變分析中不可忽略。

    (4)在疲勞損傷演變分析中,考慮腐蝕磨損影響下錨鏈的疲勞損傷會(huì)逐漸增大。在無腐蝕狀態(tài)下,磨損損傷對(duì)鏈環(huán)的疲勞損傷有一定影響,最值之差約10%左右;在磨損30 年條件下,腐蝕損傷的增加導(dǎo)致鏈環(huán)的疲勞損傷呈現(xiàn)非線性增長(zhǎng)趨勢(shì),其增長(zhǎng)速度逐漸加快,相比初始狀態(tài),134 mm 錨鏈的疲勞損傷增加了175%。因此,全生命期內(nèi)腐蝕和磨損因素均會(huì)對(duì)錨鏈疲勞強(qiáng)度產(chǎn)生影響(尤其是腐蝕),需要在實(shí)際工程中著重考慮。

    猜你喜歡
    鏈環(huán)錨鏈系泊
    考慮錨鏈腐蝕因素的錨鏈艙設(shè)計(jì)
    江蘇船舶(2023年2期)2023-06-14 11:07:44
    基于CFD方法的系泊狀態(tài)FSO舵載荷計(jì)算
    簡(jiǎn)單拓?fù)鋱D及幾乎交錯(cuò)鏈環(huán)補(bǔ)中的閉曲面
    高韌性海洋用三級(jí)系泊鏈鋼的研制
    山東冶金(2022年2期)2022-08-08 01:50:50
    錨鏈和錨鏈輪剛?cè)狁詈蟿?dòng)力學(xué)建模及嚙合過程力學(xué)分析
    船海工程(2021年6期)2021-12-17 03:17:44
    船用錨鏈發(fā)展及標(biāo)準(zhǔn)化現(xiàn)狀
    圈-雙交叉多面體鏈環(huán)的Kauffman括號(hào)多項(xiàng)式和束多項(xiàng)式
    基于懸鏈線方程的系泊系統(tǒng)分析
    關(guān)于33·42格圖的鏈環(huán)分支數(shù)計(jì)數(shù)的幾個(gè)結(jié)論
    系泊布置設(shè)計(jì)中幾個(gè)問題的討論
    船海工程(2013年6期)2013-03-11 18:57:28
    国产成+人综合+亚洲专区| 嫩草影院精品99| 日本一本二区三区精品| 搡老岳熟女国产| 很黄的视频免费| 国内少妇人妻偷人精品xxx网站| 亚洲最大成人手机在线| 757午夜福利合集在线观看| 噜噜噜噜噜久久久久久91| 小蜜桃在线观看免费完整版高清| 亚洲av成人av| 高潮久久久久久久久久久不卡| 久久婷婷人人爽人人干人人爱| 蜜桃久久精品国产亚洲av| 一a级毛片在线观看| 国产亚洲欧美98| 欧美日本亚洲视频在线播放| 夜夜夜夜夜久久久久| 国产av在哪里看| 美女 人体艺术 gogo| 嫩草影视91久久| 午夜激情欧美在线| 男人舔奶头视频| 中文字幕精品亚洲无线码一区| 成人国产一区最新在线观看| 在线看三级毛片| 欧美日本亚洲视频在线播放| 欧美黄色淫秽网站| 日本与韩国留学比较| 久久久国产成人精品二区| 丁香欧美五月| 亚洲人成网站在线播| 九九久久精品国产亚洲av麻豆| 尤物成人国产欧美一区二区三区| 高清日韩中文字幕在线| 精品久久久久久久久久免费视频| 精品午夜福利在线看| 亚洲午夜理论影院| 免费看光身美女| 身体一侧抽搐| 身体一侧抽搐| 大型黄色视频在线免费观看| 日本免费一区二区三区高清不卡| 麻豆成人午夜福利视频| 亚洲 国产 在线| 国产亚洲欧美98| 精品国产三级普通话版| 国产大屁股一区二区在线视频| 亚洲自偷自拍三级| 免费一级毛片在线播放高清视频| 国产精品自产拍在线观看55亚洲| 91九色精品人成在线观看| 国产伦人伦偷精品视频| 欧美不卡视频在线免费观看| 日韩欧美精品免费久久 | 五月伊人婷婷丁香| 人人妻人人看人人澡| 国产一区二区亚洲精品在线观看| 在线播放无遮挡| 久久久久久久精品吃奶| 99视频精品全部免费 在线| 日韩亚洲欧美综合| 禁无遮挡网站| 日韩欧美国产一区二区入口| 国产亚洲欧美98| 午夜激情欧美在线| 欧美激情在线99| 天堂动漫精品| 婷婷精品国产亚洲av在线| 国产白丝娇喘喷水9色精品| 国产 一区 欧美 日韩| 中文字幕熟女人妻在线| 国产白丝娇喘喷水9色精品| 国产人妻一区二区三区在| 久久国产乱子伦精品免费另类| 日韩欧美国产在线观看| 一个人免费在线观看的高清视频| 伦理电影大哥的女人| 99精品久久久久人妻精品| 亚洲欧美日韩无卡精品| 中文亚洲av片在线观看爽| 好看av亚洲va欧美ⅴa在| 精品久久久久久久久久久久久| 三级毛片av免费| 最近最新免费中文字幕在线| 嫩草影院精品99| 麻豆成人av在线观看| 免费高清视频大片| 午夜激情欧美在线| 日本熟妇午夜| 精品一区二区免费观看| 男插女下体视频免费在线播放| 成人美女网站在线观看视频| 国产高清视频在线播放一区| 欧美日韩中文字幕国产精品一区二区三区| 在线免费观看不下载黄p国产 | 极品教师在线免费播放| 午夜视频国产福利| 性色avwww在线观看| 国产真实伦视频高清在线观看 | 精品熟女少妇八av免费久了| 亚洲片人在线观看| 国产亚洲精品av在线| 亚洲av免费高清在线观看| 日本 av在线| 亚洲午夜理论影院| 日本一二三区视频观看| 国内揄拍国产精品人妻在线| 日韩欧美在线乱码| av女优亚洲男人天堂| 国产高清激情床上av| 在线观看舔阴道视频| 国产欧美日韩一区二区精品| 久久久久久久久大av| 日日摸夜夜添夜夜添av毛片 | www.999成人在线观看| 特级一级黄色大片| 精品无人区乱码1区二区| 日韩欧美三级三区| 欧美不卡视频在线免费观看| 97超级碰碰碰精品色视频在线观看| 宅男免费午夜| 搞女人的毛片| 精品无人区乱码1区二区| 3wmmmm亚洲av在线观看| 99精品在免费线老司机午夜| 色播亚洲综合网| 久9热在线精品视频| 99精品在免费线老司机午夜| 亚洲国产精品999在线| 波多野结衣高清无吗| 成人欧美大片| 亚洲狠狠婷婷综合久久图片| 国产精品一及| av欧美777| 2021天堂中文幕一二区在线观| 精品一区二区三区视频在线观看免费| 精品福利观看| 免费无遮挡裸体视频| 国产免费男女视频| 午夜视频国产福利| 午夜精品一区二区三区免费看| 久久久久久大精品| 欧美一区二区亚洲| 国产精品综合久久久久久久免费| 亚洲美女搞黄在线观看 | 国产亚洲av嫩草精品影院| 草草在线视频免费看| 欧美绝顶高潮抽搐喷水| 少妇人妻精品综合一区二区 | 精品久久久久久久久久久久久| 91在线精品国自产拍蜜月| 99国产极品粉嫩在线观看| 天堂网av新在线| 欧美成人一区二区免费高清观看| 18禁黄网站禁片免费观看直播| 日韩大尺度精品在线看网址| 国产精品久久视频播放| 亚洲国产精品sss在线观看| 小蜜桃在线观看免费完整版高清| 国产欧美日韩一区二区精品| 我的女老师完整版在线观看| 老熟妇乱子伦视频在线观看| 老司机午夜福利在线观看视频| 国产一区二区三区在线臀色熟女| 国产高潮美女av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 美女xxoo啪啪120秒动态图 | 久久亚洲真实| 国产一区二区在线观看日韩| 搞女人的毛片| 欧美中文日本在线观看视频| 舔av片在线| av中文乱码字幕在线| 免费看a级黄色片| 一进一出好大好爽视频| 九九久久精品国产亚洲av麻豆| 男女下面进入的视频免费午夜| 村上凉子中文字幕在线| 91午夜精品亚洲一区二区三区 | 久久国产精品人妻蜜桃| 毛片女人毛片| 床上黄色一级片| 国产三级中文精品| 午夜免费男女啪啪视频观看 | 亚洲国产精品999在线| 很黄的视频免费| 亚洲 国产 在线| 能在线免费观看的黄片| 黄色视频,在线免费观看| 色在线成人网| 成人国产综合亚洲| 内地一区二区视频在线| 91在线精品国自产拍蜜月| 3wmmmm亚洲av在线观看| 变态另类成人亚洲欧美熟女| 很黄的视频免费| 久久久久久久精品吃奶| 一进一出好大好爽视频| 色哟哟哟哟哟哟| 悠悠久久av| 午夜视频国产福利| 老司机午夜福利在线观看视频| av视频在线观看入口| www.色视频.com| 不卡一级毛片| 18+在线观看网站| 如何舔出高潮| 日韩高清综合在线| 国产黄a三级三级三级人| av在线观看视频网站免费| 757午夜福利合集在线观看| 精品久久久久久久久久久久久| 日韩亚洲欧美综合| 黄色视频,在线免费观看| 日本 av在线| 国产免费一级a男人的天堂| 欧美高清成人免费视频www| 国产单亲对白刺激| 日日干狠狠操夜夜爽| 久久精品国产清高在天天线| 亚洲成人久久性| 搡老熟女国产l中国老女人| 老司机福利观看| 久久久久九九精品影院| 亚洲国产欧美人成| 国产午夜精品久久久久久一区二区三区 | 亚洲av电影在线进入| 一级a爱片免费观看的视频| 黄色日韩在线| 欧美最黄视频在线播放免费| 亚洲人成伊人成综合网2020| 69人妻影院| 久久久久久久久久黄片| 国产免费一级a男人的天堂| 国产午夜精品论理片| 国产精品久久电影中文字幕| 午夜精品在线福利| 国产精品精品国产色婷婷| 日本 av在线| 亚洲精品亚洲一区二区| av福利片在线观看| 在线国产一区二区在线| 校园春色视频在线观看| 99精品久久久久人妻精品| 日本a在线网址| 国产高清视频在线观看网站| 午夜福利18| 欧美最新免费一区二区三区 | 欧美激情国产日韩精品一区| 又黄又爽又刺激的免费视频.| 国产成人av教育| 性插视频无遮挡在线免费观看| 天堂√8在线中文| 美女免费视频网站| 免费在线观看影片大全网站| 老司机深夜福利视频在线观看| 国产在视频线在精品| 狂野欧美白嫩少妇大欣赏| 国产精品久久视频播放| 亚洲国产欧美人成| 婷婷六月久久综合丁香| 嫩草影视91久久| 欧美区成人在线视频| 亚洲无线观看免费| 精品久久久久久久久久免费视频| 国产色婷婷99| 成人av一区二区三区在线看| 亚洲精品粉嫩美女一区| 亚洲精品成人久久久久久| 国产成年人精品一区二区| 中亚洲国语对白在线视频| 成人鲁丝片一二三区免费| 国产av麻豆久久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 俄罗斯特黄特色一大片| 久久精品国产清高在天天线| 免费看a级黄色片| 日本黄大片高清| 老熟妇乱子伦视频在线观看| av专区在线播放| 免费黄网站久久成人精品 | 乱码一卡2卡4卡精品| 欧美成人免费av一区二区三区| 51国产日韩欧美| a级一级毛片免费在线观看| 欧美日韩亚洲国产一区二区在线观看| 别揉我奶头~嗯~啊~动态视频| 国产精品久久视频播放| 91在线精品国自产拍蜜月| 国产三级黄色录像| 香蕉av资源在线| 欧美在线黄色| 午夜免费成人在线视频| 全区人妻精品视频| 99视频精品全部免费 在线| 黄色女人牲交| 性色av乱码一区二区三区2| 亚洲精品久久国产高清桃花| 一本久久中文字幕| 国产精品亚洲av一区麻豆| 真实男女啪啪啪动态图| 欧美绝顶高潮抽搐喷水| 91狼人影院| 国产伦在线观看视频一区| 欧美日本视频| 人人妻人人看人人澡| 免费一级毛片在线播放高清视频| 日韩大尺度精品在线看网址| 国产色爽女视频免费观看| 天堂影院成人在线观看| 色尼玛亚洲综合影院| 校园春色视频在线观看| 欧美成人免费av一区二区三区| 国产精品伦人一区二区| 欧美色欧美亚洲另类二区| 老女人水多毛片| 成人高潮视频无遮挡免费网站| 精品不卡国产一区二区三区| 村上凉子中文字幕在线| 日韩欧美三级三区| 欧美成人性av电影在线观看| 在线看三级毛片| 蜜桃亚洲精品一区二区三区| 性欧美人与动物交配| 国产精品亚洲av一区麻豆| 女人被狂操c到高潮| 国产精品亚洲av一区麻豆| 男人狂女人下面高潮的视频| 国产精品亚洲av一区麻豆| 国产精品伦人一区二区| 如何舔出高潮| 校园春色视频在线观看| 内地一区二区视频在线| 国产在视频线在精品| 欧美一级a爱片免费观看看| 窝窝影院91人妻| 午夜日韩欧美国产| 久久精品国产亚洲av涩爱 | 国产单亲对白刺激| 特大巨黑吊av在线直播| 中文字幕免费在线视频6| 国产单亲对白刺激| 亚洲avbb在线观看| 欧美国产日韩亚洲一区| 高清在线国产一区| 日本一本二区三区精品| 亚洲av成人精品一区久久| 精品国内亚洲2022精品成人| 日日摸夜夜添夜夜添小说| 国产欧美日韩精品一区二区| 看黄色毛片网站| 国产熟女xx| 日韩精品中文字幕看吧| 亚洲经典国产精华液单 | 国内毛片毛片毛片毛片毛片| 少妇裸体淫交视频免费看高清| 性色av乱码一区二区三区2| 午夜老司机福利剧场| 在线国产一区二区在线| 3wmmmm亚洲av在线观看| 一级作爱视频免费观看| 搞女人的毛片| 99国产精品一区二区三区| 国产精品一区二区三区四区免费观看 | 国产欧美日韩精品一区二区| 中文字幕av在线有码专区| 国产伦一二天堂av在线观看| 亚洲国产精品合色在线| 国产老妇女一区| 国产免费一级a男人的天堂| 最好的美女福利视频网| 男人舔女人下体高潮全视频| 国产精品久久久久久精品电影| 我的女老师完整版在线观看| 亚洲熟妇中文字幕五十中出| 欧美性猛交黑人性爽| 国产高清有码在线观看视频| 最近最新免费中文字幕在线| 一本一本综合久久| 免费在线观看亚洲国产| 老司机福利观看| 观看免费一级毛片| 国产探花在线观看一区二区| 欧美性猛交黑人性爽| 乱人视频在线观看| 性插视频无遮挡在线免费观看| 午夜福利视频1000在线观看| 很黄的视频免费| 免费黄网站久久成人精品 | 亚洲国产精品999在线| 亚洲av.av天堂| 亚洲av电影在线进入| 中亚洲国语对白在线视频| 亚洲人成电影免费在线| 成年女人毛片免费观看观看9| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 99精品在免费线老司机午夜| 老司机深夜福利视频在线观看| 熟女电影av网| 午夜免费成人在线视频| 欧美午夜高清在线| 久久久精品欧美日韩精品| 最近最新免费中文字幕在线| 日韩有码中文字幕| 亚洲久久久久久中文字幕| 国产欧美日韩精品一区二区| 白带黄色成豆腐渣| 国产乱人视频| 久久中文看片网| 51午夜福利影视在线观看| 性色avwww在线观看| 在线看三级毛片| 高清在线国产一区| 窝窝影院91人妻| 久久久久国产精品人妻aⅴ院| 国产精品一区二区性色av| 97碰自拍视频| 搡女人真爽免费视频火全软件 | 国产欧美日韩一区二区三| 国产av不卡久久| 日本a在线网址| 能在线免费观看的黄片| 亚洲久久久久久中文字幕| 日日摸夜夜添夜夜添小说| 亚洲国产精品合色在线| 人妻夜夜爽99麻豆av| 午夜福利欧美成人| 亚洲激情在线av| 性色av乱码一区二区三区2| 久久九九热精品免费| 亚洲欧美精品综合久久99| 看免费av毛片| .国产精品久久| 色在线成人网| 18禁在线播放成人免费| 尤物成人国产欧美一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 长腿黑丝高跟| 一级黄片播放器| 国产成+人综合+亚洲专区| av欧美777| 久久性视频一级片| 久久国产精品人妻蜜桃| 一个人看的www免费观看视频| 免费看a级黄色片| 国产亚洲欧美在线一区二区| 嫁个100分男人电影在线观看| 亚洲av不卡在线观看| 亚洲午夜理论影院| 老熟妇仑乱视频hdxx| 99热只有精品国产| 嫩草影院入口| 蜜桃亚洲精品一区二区三区| 91在线观看av| 91午夜精品亚洲一区二区三区 | 亚洲精品一区av在线观看| 国产国拍精品亚洲av在线观看| 午夜免费激情av| 欧美成人免费av一区二区三区| 一区二区三区高清视频在线| www.www免费av| 少妇的逼水好多| 亚洲中文字幕日韩| 日韩欧美免费精品| 免费一级毛片在线播放高清视频| 色噜噜av男人的天堂激情| 国产淫片久久久久久久久 | 99热只有精品国产| 中文在线观看免费www的网站| 久久精品久久久久久噜噜老黄 | 一级a爱片免费观看的视频| 亚洲人与动物交配视频| 在线天堂最新版资源| 嫩草影院入口| 99视频精品全部免费 在线| 免费观看的影片在线观看| 亚洲国产高清在线一区二区三| 国产高清有码在线观看视频| 欧美不卡视频在线免费观看| 日日干狠狠操夜夜爽| 757午夜福利合集在线观看| 日日干狠狠操夜夜爽| 国产精品乱码一区二三区的特点| 深夜a级毛片| 亚洲欧美精品综合久久99| 久久亚洲精品不卡| 亚洲不卡免费看| 久久久久国内视频| 日本熟妇午夜| а√天堂www在线а√下载| 91狼人影院| 欧美一区二区国产精品久久精品| 亚洲人成网站高清观看| 国产主播在线观看一区二区| 又爽又黄无遮挡网站| 中文字幕av成人在线电影| 午夜精品在线福利| 成人高潮视频无遮挡免费网站| 日韩欧美三级三区| 天堂网av新在线| 欧美黄色片欧美黄色片| 国产成人欧美在线观看| 男人舔奶头视频| 非洲黑人性xxxx精品又粗又长| 色哟哟哟哟哟哟| 一本一本综合久久| 白带黄色成豆腐渣| 九色成人免费人妻av| 午夜久久久久精精品| 成人三级黄色视频| 亚洲美女视频黄频| 国产v大片淫在线免费观看| 男女下面进入的视频免费午夜| 一卡2卡三卡四卡精品乱码亚洲| 国产精品综合久久久久久久免费| 99riav亚洲国产免费| 国产伦一二天堂av在线观看| 俄罗斯特黄特色一大片| 精品国内亚洲2022精品成人| 亚洲熟妇中文字幕五十中出| 国内精品久久久久精免费| 精品人妻视频免费看| 久久久国产成人精品二区| 日韩免费av在线播放| 久9热在线精品视频| 麻豆成人av在线观看| 俺也久久电影网| 天天一区二区日本电影三级| 69人妻影院| 在线国产一区二区在线| 琪琪午夜伦伦电影理论片6080| 观看免费一级毛片| 亚洲18禁久久av| 国内久久婷婷六月综合欲色啪| 亚洲 欧美 日韩 在线 免费| 亚洲国产精品成人综合色| 日韩欧美在线乱码| 少妇被粗大猛烈的视频| 午夜激情福利司机影院| 可以在线观看毛片的网站| 91av网一区二区| 无人区码免费观看不卡| 欧美区成人在线视频| 91在线精品国自产拍蜜月| 搡老熟女国产l中国老女人| 亚洲av电影在线进入| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 天天一区二区日本电影三级| 亚洲黑人精品在线| 九九在线视频观看精品| 97热精品久久久久久| 欧美日韩瑟瑟在线播放| 久久精品人妻少妇| 午夜免费激情av| 亚洲色图av天堂| 久久天躁狠狠躁夜夜2o2o| 一个人看视频在线观看www免费| 琪琪午夜伦伦电影理论片6080| 亚洲中文日韩欧美视频| 国产精品一区二区三区四区久久| 又黄又爽又刺激的免费视频.| 欧美日韩乱码在线| av福利片在线观看| 美女高潮的动态| 3wmmmm亚洲av在线观看| 波多野结衣巨乳人妻| 国产精品女同一区二区软件 | 18禁黄网站禁片免费观看直播| 色av中文字幕| 天堂√8在线中文| 亚洲国产日韩欧美精品在线观看| 男插女下体视频免费在线播放| 亚洲成人久久性| 欧美另类亚洲清纯唯美| 国产黄a三级三级三级人| 久久久久久久久久成人| 1024手机看黄色片| 男人舔女人下体高潮全视频| 欧美高清成人免费视频www| 禁无遮挡网站| 一a级毛片在线观看| 赤兔流量卡办理| 一本久久中文字幕| 婷婷色综合大香蕉| 国产激情偷乱视频一区二区| 久久午夜福利片| 国产三级在线视频| 亚洲三级黄色毛片| 久久午夜福利片| 国产真实伦视频高清在线观看 | 亚洲国产高清在线一区二区三| 国产v大片淫在线免费观看| 亚洲,欧美精品.| 国产精品一区二区三区四区免费观看 | 精品一区二区三区av网在线观看| 丰满的人妻完整版| 真人做人爱边吃奶动态| 一夜夜www| 久久久久久九九精品二区国产| 精品久久久久久成人av| 91麻豆精品激情在线观看国产| 9191精品国产免费久久| 午夜免费激情av| 国产乱人视频| av黄色大香蕉| 一夜夜www| 日韩欧美国产一区二区入口| 99久国产av精品| a级毛片a级免费在线| 不卡一级毛片| 老司机午夜十八禁免费视频| 久久久久久久久久黄片| 成年版毛片免费区| 男插女下体视频免费在线播放|