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

    船體縱骨典型節(jié)點(diǎn)疲勞裂紋擴(kuò)展壽命評(píng)估

    2016-05-04 18:40:28何文濤劉敬喜
    船舶力學(xué) 2016年11期
    關(guān)鍵詞:縱骨尖端腹板

    何文濤,劉敬喜,2,解 德,2

    (1.華中科技大學(xué) 船舶與海洋工程學(xué)院,武漢 430074;2.船舶和海洋水動(dòng)力湖北省重點(diǎn)實(shí)驗(yàn)室,武漢 430074)

    船體縱骨典型節(jié)點(diǎn)疲勞裂紋擴(kuò)展壽命評(píng)估

    何文濤1,劉敬喜1,2,解 德1,2

    (1.華中科技大學(xué) 船舶與海洋工程學(xué)院,武漢 430074;2.船舶和海洋水動(dòng)力湖北省重點(diǎn)實(shí)驗(yàn)室,武漢 430074)

    基于有限元軟件ABAQUS,結(jié)合虛擬裂紋閉合法、裂紋擴(kuò)展判據(jù)及子結(jié)構(gòu)技術(shù),應(yīng)用腳本語言Python開發(fā)出模擬疲勞裂紋擴(kuò)展的程序(FCG-System)。對(duì)含初始裂紋的油船縱骨典型節(jié)點(diǎn)在側(cè)面壓力作用下進(jìn)行疲勞裂紋擴(kuò)展數(shù)值模擬,并探討了軟趾、背肘板及防傾肘板對(duì)疲勞裂紋擴(kuò)展路徑和壽命的影響。結(jié)果表明,增設(shè)軟趾、背肘板或防傾肘板都會(huì)使裂紋擴(kuò)展路徑曲率增大,且軟趾、防傾肘板可使裂紋擴(kuò)展壽命增大,背肘板可使裂紋擴(kuò)展壽命減小。

    裂紋擴(kuò)展路徑;裂紋擴(kuò)展壽命;虛擬裂紋閉合法;子結(jié)構(gòu)技術(shù);FCG-System;

    0 引 言

    隨著船舶大型化及高強(qiáng)度鋼在船體結(jié)構(gòu)中的廣泛應(yīng)用,船體結(jié)構(gòu)疲勞強(qiáng)度的評(píng)估顯得尤為重要。日本船級(jí)社[1]對(duì)上世紀(jì)八十年代中期采用大量高強(qiáng)度鋼建造并入該船級(jí)社的第二代VLCC進(jìn)行了資料收集。這些資料表明,76條船舶中的65條有不同程度的破損,其中縱骨與強(qiáng)框架或橫艙壁交界處的疲勞損傷情況尤為嚴(yán)重,約占全部損傷統(tǒng)計(jì)的70%。

    目前,各船級(jí)社的疲勞強(qiáng)度校核主要是基于S-N曲線的累積損傷方法[2-4],而針對(duì)船體結(jié)構(gòu)中疲勞裂紋擴(kuò)展路徑、擴(kuò)展速率以及擴(kuò)展壽命的研究較少[5-6]。Okawa[7]開發(fā)了數(shù)值模擬程序?qū)Υw縱骨焊趾處疲勞裂紋擴(kuò)展進(jìn)行研究,在模擬過程中將穿透型裂紋看作為二維平面問題。法國一個(gè)公司[8]采用線彈簧法計(jì)算應(yīng)力強(qiáng)度因子,進(jìn)而對(duì)疲勞壽命進(jìn)行評(píng)估,并開發(fā)了軟件SAPHIRS。然而,該方法需預(yù)先知道裂紋擴(kuò)展路徑。Jang[9]通過比較實(shí)驗(yàn)數(shù)據(jù)和現(xiàn)存的公式,提出了若干預(yù)測疲勞裂紋擴(kuò)展的公式。王麗麗[5]對(duì)FPSO縱骨處肘板上兩裂紋擴(kuò)展路徑及擴(kuò)展速率進(jìn)行研究,采用等效應(yīng)力強(qiáng)度因子計(jì)算裂紋擴(kuò)展壽命,然而裂紋僅分布在肘板平面內(nèi)。

    本文作者已通過面向?qū)ο蟮某绦蛘Z言Python及ABAQUS腳本接口開發(fā)了疲勞裂紋擴(kuò)展程序(FCG-System),能夠模擬復(fù)雜加載下混合模式的裂紋擴(kuò)展路徑和擴(kuò)展壽命[10]。本文主要對(duì)含初始裂紋的油船縱骨典型節(jié)點(diǎn)進(jìn)行疲勞裂紋擴(kuò)展數(shù)值模擬,并探討了軟趾、背肘板及防傾肘板對(duì)疲勞裂紋擴(kuò)展路徑和壽命的影響。

    1 基本原理

    圖1為該疲勞裂紋擴(kuò)展程序的流程圖?,F(xiàn)將該程序的主要執(zhí)行步驟總結(jié)如下:

    (1)前處理:建立模型、邊界條件、定義初始裂紋和網(wǎng)格劃分,為后面有限元計(jì)算和參數(shù)求解做準(zhǔn)備。為了模擬大型海洋結(jié)構(gòu)物疲勞裂紋擴(kuò)展,將裂紋擴(kuò)展區(qū)以外其他區(qū)域采用子結(jié)構(gòu)技術(shù)建模。

    (2)有限元分析:建立有限元模型并提交給Abaqus/standard進(jìn)行有限元分析,輸出計(jì)算結(jié)果。

    (3)應(yīng)力強(qiáng)度因子計(jì)算:從.odb文件中提取裂紋尖端的節(jié)點(diǎn)力及節(jié)點(diǎn)位移,并結(jié)合虛擬裂紋閉合法計(jì)算應(yīng)力強(qiáng)度因子;

    (4)裂紋路徑預(yù)測:通過最大周向應(yīng)力準(zhǔn)則預(yù)測疲勞裂紋擴(kuò)展;

    (5)裂紋擴(kuò)展計(jì)算:通過裂紋擴(kuò)展速率公式計(jì)算疲勞裂紋擴(kuò)展量;

    圖1 FCG-System流程圖Fig.1 Flow chart of FCG-System

    (6)下一步模擬:返回(1),更新裂紋的幾何形狀和拓?fù)浣Y(jié)構(gòu),為下一步循環(huán)中網(wǎng)格重新劃分做準(zhǔn)備。

    (7)后處理:將上述計(jì)算結(jié)果輸出到外部文件,如有必要還可以調(diào)用模塊畫出圖表。

    圖2 裂紋擴(kuò)展區(qū)網(wǎng)格劃分Fig.2 Mesh pattern of crack growth domain

    圖3 裂紋擴(kuò)展方向及裂紋尖端區(qū)域網(wǎng)格劃分Fig.3 Crack growth direction and rosette meshes of the crack tip

    1.1 自動(dòng)網(wǎng)格劃分

    在有限元建模時(shí),將結(jié)構(gòu)分為全局模型和局部模型,其中局部模型為裂紋擴(kuò)展區(qū),全局模型為除裂紋擴(kuò)展區(qū)以外的其他區(qū)域,通常采用子結(jié)構(gòu)技術(shù)建模。在擴(kuò)展過程中,僅對(duì)裂紋擴(kuò)展區(qū)域進(jìn)行網(wǎng)格重新劃分,而全局模型僅在第一個(gè)分析步中劃分網(wǎng)格。此外,裂紋擴(kuò)展區(qū)又分為裂紋尖端區(qū)域和裂紋尖端周圍區(qū)域,如圖2所示。在網(wǎng)格劃分時(shí)采用非均勻的網(wǎng)格密度,其中裂紋尖端區(qū)域采用密網(wǎng)格,而粗網(wǎng)格涵蓋其余的子域。采用這種劃分方式,可以達(dá)到計(jì)算效率和計(jì)算精度之間的平衡。

    值得注意的是,本文采用虛擬裂紋閉合法計(jì)算裂紋尖端的應(yīng)力強(qiáng)度因子,所以無需奇異的裂紋尖端單元,整個(gè)裂紋擴(kuò)展區(qū)都采用ABAQUS標(biāo)準(zhǔn)三角形單元。

    1.2 應(yīng)力強(qiáng)度因子計(jì)算

    對(duì)于船舶及海洋結(jié)構(gòu)物來說,全部使用三維實(shí)體單元代價(jià)太大,因此常用殼體單元代替塊體單元。Wang和Raju[11]提出了一些公式,用來研究加強(qiáng)筋和蒙皮的分離,但是這種公式只適用于面狀裂紋界面分離的情況,對(duì)于船體結(jié)構(gòu)中常見的線狀裂紋不適用。為此,本文采用適用于線狀裂紋的虛擬裂紋閉合法公式[10]。對(duì)于混合型加載,能量釋放率GⅠ、GⅡ和GⅢ能夠通過裂紋尖端的節(jié)點(diǎn)1的節(jié)點(diǎn)力和尖端后節(jié)點(diǎn)3、4的張開位移計(jì)算得到,如圖3所示。因此,各能量釋放率分量在局部坐標(biāo)系下可表示為:

    在混合加載模式下,為了考慮斷裂模式Ⅰ、Ⅱ和Ⅲ的共同作用需引入一個(gè)等效應(yīng)力強(qiáng)度因子Keff:

    1.3 裂紋擴(kuò)展準(zhǔn)則及裂紋擴(kuò)展角

    對(duì)于各向彈性材料,最常用的擴(kuò)展準(zhǔn)則為最大周向應(yīng)力準(zhǔn)則,其基本假設(shè)為:(1)裂紋向周向應(yīng)力最大的方向擴(kuò)展;(2)當(dāng)此方向的周向應(yīng)力大于臨界值時(shí),裂紋失穩(wěn)擴(kuò)展。

    裂紋尖端周向應(yīng)力的表達(dá)式[12]為:

    式中:r和θ為裂紋尖端的極距和角度。裂紋擴(kuò)展角θ0可由(4)式對(duì)θ一階偏導(dǎo)求得:

    因此,可得裂紋擴(kuò)展角θ0:

    1.4 裂紋擴(kuò)展率公式

    在Paris公式的基礎(chǔ)上,考慮應(yīng)力比、應(yīng)力強(qiáng)度因子門檻值的影響,得出了很多種修正公式,本文采用如下的裂紋擴(kuò)展公式[7]:

    其中:R為應(yīng)力比。因此,裂紋增長量可表示為:

    2 驗(yàn)證例題

    如圖4所示,以含裂紋的三維板結(jié)構(gòu)為例研究裂紋擴(kuò)展問題。此結(jié)構(gòu)的兩個(gè)平面各含有一個(gè)裂紋尖端,在單軸拉伸載荷作用下兩條裂紋同時(shí)擴(kuò)展。此結(jié)構(gòu)在有限元計(jì)算中分為兩部分:一部分為裂紋擴(kuò)展區(qū),此區(qū)域含有裂紋及裂紋擴(kuò)展區(qū)域;另一部分為子結(jié)構(gòu)區(qū),即除裂紋擴(kuò)展區(qū)以外的其他區(qū)域。裂紋擴(kuò)展區(qū)與子結(jié)構(gòu)區(qū)之間通過子結(jié)構(gòu)技術(shù)中的保留節(jié)點(diǎn)相連接。該結(jié)構(gòu)的有限元模型如圖5所示,整個(gè)模型均采用殼體單元。

    圖4 含裂紋的三維板結(jié)構(gòu)(單位:mm)Fig.4 A cracked three-dimensional panel structure (unit:mm)

    圖5 有限元模型:(a)子結(jié)構(gòu)區(qū);(b)裂紋擴(kuò)展區(qū)Fig.5 Finite element model:(a)Substructure domain;(b)Crack growth domain

    圖6給出了采用本文程序計(jì)算的應(yīng)力強(qiáng)度因子隨裂紋擴(kuò)展量變化的曲線。此曲線與參考文獻(xiàn)[7]中計(jì)算結(jié)果吻合得很好。

    3 典型節(jié)點(diǎn)疲勞裂紋擴(kuò)展分析

    油船船體由于受到波浪誘導(dǎo)載荷引起的船體總縱彎曲、波浪載荷引起的動(dòng)壓力以及艙內(nèi)液艙和壓載水的慣性力作用,使得縱骨與橫框架及橫艙壁的焊接處易于產(chǎn)生疲勞裂紋并沿縱骨的面板或腹板擴(kuò)展。本部分采用開發(fā)的疲勞裂紋擴(kuò)展程序?qū)v骨與橫艙壁連接節(jié)點(diǎn)在側(cè)面壓力作用下進(jìn)行疲勞裂紋擴(kuò)展及壽命分析,探討軟趾、含背肘板的扶強(qiáng)材、防傾肘板的加強(qiáng)結(jié)構(gòu)對(duì)疲勞裂紋的影響??v骨與橫艙壁的連接處的局部模型如圖7所示。此部分采用子結(jié)構(gòu)建模,橫向間距為800 mm,縱向間距為4 400 mm。模型的材料參數(shù)及裂紋擴(kuò)展模型參數(shù)列于表1。

    圖6 應(yīng)力強(qiáng)度因子比較Fig.6 Comparison of stress intensity factors

    圖7 縱骨結(jié)構(gòu)及裂紋擴(kuò)展區(qū)Fig.7 The longitudinal structure and its crack growth domain

    表1 計(jì)算模型的材料屬性(單位:應(yīng)力/MPa,長度/m)Tab.1 Material properties of the analysis model(unit:stress/MPa,length/m)

    根據(jù)縱骨與橫艙壁連接形式的不同,本文對(duì)四種常用的含扶強(qiáng)材或防傾肘板的節(jié)點(diǎn)進(jìn)行疲勞裂紋擴(kuò)展分析,研究不同的加強(qiáng)方式及軟趾設(shè)置對(duì)裂紋擴(kuò)展路徑及擴(kuò)展壽命的影響。表2對(duì)所涉及的四種結(jié)構(gòu)形式進(jìn)行了簡單描述,各種結(jié)構(gòu)形式示意圖[7]如圖8所示。

    表2 節(jié)點(diǎn)結(jié)構(gòu)形式描述Tab.2 Description of connection details

    圖8 節(jié)點(diǎn)結(jié)構(gòu)形式(單位:mm)Fig.8 Details of connection

    3.1 裂紋擴(kuò)展過程描述

    圖9給出了裂紋的擴(kuò)展過程。在整個(gè)裂紋擴(kuò)展過程中,裂紋被認(rèn)為是穿透性裂紋,模型采用殼體單元。在裂紋時(shí),一旦有裂紋擴(kuò)展至板邊致其斷裂,其余板中的裂紋尖端需重新定義,而裂紋擴(kuò)展區(qū)之外的其他區(qū)域不需要重新定義。也就是說,在整個(gè)裂紋擴(kuò)展過程中,僅對(duì)裂紋擴(kuò)展區(qū)進(jìn)行網(wǎng)格重劃分,而子結(jié)構(gòu)區(qū)僅在裂紋擴(kuò)展分析的第一步中進(jìn)行網(wǎng)格劃分。裂紋起源于面板與扶強(qiáng)材相交處,且含有三個(gè)裂紋尖端。假設(shè)初始裂紋長度為:面板A,20 mm;面板B,10 mm;腹板,10 mm(如圖9(a))。該裂紋在面板與腹板中同時(shí)擴(kuò)展,直至擴(kuò)展至縱骨與橫艙壁的交線或縱骨與內(nèi)底板的交線處。這類裂紋的擴(kuò)展過程分為三個(gè)階段:

    (1)在開始時(shí),裂紋在腹板與面板中同時(shí)擴(kuò)展,此時(shí)具有三個(gè)裂紋尖端(如圖9(b))。

    (2)面板A斷裂后,裂紋繼續(xù)在面板B與腹板中擴(kuò)展,此時(shí)裂紋具有兩個(gè)裂紋尖端(如圖9(c))。

    (3)面板B斷裂后,裂紋繼續(xù)在腹板中擴(kuò)展,直至達(dá)到縱骨與底板的交線處。此時(shí)裂紋僅含有一個(gè)裂紋尖端(如圖9(d))。

    圖9 有限元模型及裂紋擴(kuò)展過程Fig.9 The finite element model and its crack growth process

    圖10 面板裂紋擴(kuò)展路徑Fig.10 Crack growth paths in flange

    圖11 腹板裂紋擴(kuò)展路徑Fig.11 Crack growth paths in web-plate

    3.2 裂紋擴(kuò)展路徑預(yù)測

    各種形式典型節(jié)點(diǎn)的裂紋擴(kuò)展路徑如圖10和11所示,坐標(biāo)原點(diǎn)為中間縱骨、底板與橫艙壁的連接處。由此可見,無論是何種形式的節(jié)點(diǎn),面板中的裂紋都基本沿直線擴(kuò)展。而在腹板中可明顯看出不同形式典型節(jié)點(diǎn)對(duì)裂紋擴(kuò)展路徑的影響:節(jié)點(diǎn)A中裂紋轉(zhuǎn)折最小,節(jié)點(diǎn)B、C次之,節(jié)點(diǎn)D最大。結(jié)果表明,軟趾、防傾肘板及背肘板對(duì)裂紋擴(kuò)展路徑有較大的影響。增設(shè)軟趾、背肘板或設(shè)置防傾肘板都會(huì)使裂紋擴(kuò)展路徑曲率增大,且向縱骨及橫艙壁的連接處擴(kuò)展。此外,還可以看出,面板A斷裂后,裂紋路徑發(fā)生偏折,其中節(jié)點(diǎn)D最為明顯。

    3.3 應(yīng)力強(qiáng)度因子分析

    圖12、13分別給出了各種形式典型節(jié)點(diǎn)中面板和腹板的應(yīng)力強(qiáng)度因子隨裂紋擴(kuò)展長度的變化曲線??擅黠@看出,在相同裂紋長度處節(jié)點(diǎn)A中應(yīng)力強(qiáng)度因子最大,節(jié)點(diǎn)D次之,節(jié)點(diǎn)B、C最小且相差不大。結(jié)果表明,軟趾、防傾肘板及背肘板對(duì)應(yīng)力強(qiáng)度因子有較大的影響。由節(jié)點(diǎn)A和B可以看出,軟趾可使裂紋尖端應(yīng)力強(qiáng)度因子減?。挥晒?jié)點(diǎn)C和D可以看出,背肘板可使裂紋尖端的應(yīng)力強(qiáng)度因子增加;由節(jié)點(diǎn)A和C可以看出,防傾肘板可使裂紋尖端的應(yīng)力強(qiáng)度因子減小。

    圖12 面板應(yīng)力強(qiáng)度因子的變化 (a)面板A;(b)面板BFig.12 Variation of stress intensity factors in flange(a)Side-A;(b)Side-B

    圖13 腹板應(yīng)力強(qiáng)度因子的變化Fig.13 Variation of stress intensity factors in web-plate

    此外,可發(fā)現(xiàn)面板B中裂紋在擴(kuò)展至大約32~35 mm處面板A斷裂。在面板A斷裂處,其余板中的應(yīng)力強(qiáng)度因子發(fā)生跳躍式增大。這是由于面板A斷裂后,結(jié)構(gòu)中的應(yīng)力重新分布,造成面板B和腹板承受較大應(yīng)力。在模擬開始時(shí),三條裂紋分別在面板兩側(cè)及腹板中同時(shí)擴(kuò)展,其相應(yīng)的應(yīng)力強(qiáng)度因子相對(duì)較小。面板斷裂后,結(jié)構(gòu)中的應(yīng)力重新分配,腹板中裂紋尖端的應(yīng)力強(qiáng)度因子急劇增加。在裂紋擴(kuò)展至接近底板處,節(jié)點(diǎn)A中的應(yīng)力強(qiáng)度因子減小,裂紋擴(kuò)展速率減小,而其他形式的節(jié)點(diǎn)在裂紋擴(kuò)展至接近底板時(shí)出現(xiàn)速率增大的現(xiàn)象,這是由于不同形式典型節(jié)點(diǎn)的裂紋擴(kuò)展路徑不同導(dǎo)致的。

    3.4 裂紋擴(kuò)展壽命預(yù)報(bào)

    圖14給出了各種形式典型節(jié)點(diǎn)的面板和腹板裂紋擴(kuò)展壽命隨裂紋長度的變化曲線。曲線上某一點(diǎn)的斜率表示該點(diǎn)的裂紋擴(kuò)展速率,可見隨著循環(huán)次數(shù)的增加,裂紋擴(kuò)展速率逐漸加快。例如,節(jié)點(diǎn)A在循環(huán)載荷作用下初始裂紋橫穿整個(gè)腹板的循環(huán)次數(shù)為1.70×106次,此壽命可認(rèn)為是構(gòu)件破壞的總壽命。從初始裂紋至面板A斷裂的循環(huán)次數(shù)為1.10×106次,其擴(kuò)展壽命占總壽命的64.7%;至面板B斷裂的循環(huán)次數(shù)為1.19×106,其擴(kuò)展壽命占總壽命的70.0%?,F(xiàn)將各種形式典型節(jié)點(diǎn)的疲勞擴(kuò)展壽命總結(jié)如表3所示。

    可明顯看出,不同形式典型節(jié)點(diǎn)對(duì)裂紋擴(kuò)展壽命的影響:節(jié)點(diǎn)A的疲勞壽命最小,節(jié)點(diǎn)D、C次之,節(jié)點(diǎn)B最大。結(jié)果表明,軟趾、防傾肘板及背肘板對(duì)裂紋擴(kuò)展壽命有較大的影響。由節(jié)點(diǎn)A和B可以看出,軟趾使裂紋擴(kuò)展壽命增大;由節(jié)點(diǎn)C和D可以看出,背肘板可使裂紋擴(kuò)展壽命減小;由節(jié)點(diǎn)A和C可以看出,防傾肘板可使裂紋擴(kuò)展壽命增大。

    圖14 面板和腹板裂紋擴(kuò)展壽命Fig.14 Crack growth lives in flange and web-plate

    此外還可以發(fā)現(xiàn),從初始裂紋至面板A斷裂這一過程裂紋擴(kuò)展速率較慢,其擴(kuò)展壽命占總壽命的60%以上,至面板B斷裂的擴(kuò)展壽命占總壽命的70%以上。因此,面板斷裂前的擴(kuò)展壽命占總壽命的比重很大。一旦面板A斷裂,腹板中裂紋擴(kuò)展速率迅速增大,裂紋很快橫穿整個(gè)腹板,造成構(gòu)件破壞。

    表3 各節(jié)點(diǎn)裂紋擴(kuò)展壽命Tab.3 Crack growth lives of each connection detail

    位于面板及腹板上的初始裂紋在循環(huán)載荷作用下不斷擴(kuò)展,當(dāng)面板斷裂后腹板中的裂紋會(huì)迅速擴(kuò)展至船體底板或舷側(cè)板,不僅造成了結(jié)構(gòu)承載能力的下降,還會(huì)造成水/油密部件泄漏等災(zāi)難性事故。因此有關(guān)人員應(yīng)加強(qiáng)這類裂紋的管理,對(duì)容易發(fā)生疲勞裂紋的部位進(jìn)行定期檢測,并采取相應(yīng)的措施修復(fù)和制止裂紋擴(kuò)展,以確保構(gòu)件的完整性和連續(xù)性。

    4 結(jié) 論

    本文基于斷裂力學(xué)的方法對(duì)船體縱骨與橫艙壁連接處的典型節(jié)點(diǎn)進(jìn)行了疲勞裂紋擴(kuò)展及壽命分析,得出了如下結(jié)論:

    (1)軟趾、防傾肘板及背肘板對(duì)裂紋擴(kuò)展路徑有較大的影響。增設(shè)軟趾、背肘板或設(shè)置防傾肘板都會(huì)使裂紋擴(kuò)展路徑曲率增大,且向縱骨及橫艙壁的連接處擴(kuò)展。

    (2)不同形式典型節(jié)點(diǎn)對(duì)裂紋尖端應(yīng)力強(qiáng)度因子有較大的影響。軟趾、防傾肘板可使裂紋尖端應(yīng)力強(qiáng)度因子減??;背肘板可使裂紋尖端的應(yīng)力強(qiáng)度因子增加。

    (3)不同形式典型節(jié)點(diǎn)對(duì)裂紋擴(kuò)展壽命有明顯影響。軟趾、防傾肘板可使裂紋擴(kuò)展壽命增大;背肘板可使裂紋擴(kuò)展壽命減小。

    (4)有關(guān)人員應(yīng)加強(qiáng)這類裂紋的管理,對(duì)容易發(fā)生疲勞裂紋的部位進(jìn)行定期檢測,并采取相應(yīng)的措施修復(fù)和制止裂紋擴(kuò)展,以確保構(gòu)件的完整性和連續(xù)性。

    [1]NK.Casualty review-comparative damage review of 2nd generation VLCCs[R].1998.

    [2]DNV.Fatigue strength analysis of offshore steel structures(RP-203)[S].Det Norske Veritas,Hovik,2005.

    [3]ABS.Rules for building and classing steel vessels[S].American Bureau of Shipping,Houston,2002.

    [4]崔維成,蔡新剛,冷建興.船舶結(jié)構(gòu)疲勞強(qiáng)度校核研究現(xiàn)狀及我國的進(jìn)展[J].船舶力學(xué),1998,2(4):63-81. Cui Weicheng,Cai Xingang,Leng Jianxing.A state-of-the-art review for the fatigue strength assessment of ship structures[J].Journal of Ship Mechanics,1998,2(4):63-81.

    [5]王麗麗,黃小平,崔維成.復(fù)雜應(yīng)力場中裂紋疲勞擴(kuò)展壽命預(yù)報(bào)[J].船舶力學(xué),2011,15(4):383-388. Wang Lili,Huang Xiaoping,Cui Weicheng.Fatigue crack growth life prediction of cracks in complex stress field[J].Journal of Ship Mechanics,2011,15(4):383-388.

    [6]Cui Weicheng.A feasible study of fatigue life prediction for marine structures based on crack propagation analysis[J].Journal of Engineering for the Maritime Environment,2003,217(5):11-23.

    [7]Okawa T,Sumi Y,Mohri M.Simulation-based fatigue crack management of ship structural details applied to longitudinal and transverse connections[J].Marine Structures,2006,19(4):217-240.

    [8]Serror M,Marchal N.Simulation of behavior of fatigue cracks:a complete industrial process on a typical connection in a FPSO[C]//Proceedings of the ASME 2009 28th International Conference on Ocean,Offshore and Arctic Engineering OMAE2009 May 31-June 5,2009.Honolulu,Hawaii,2009.

    [9]Jang B S,Ito H,Kim K S,Suh Y S,Jeon H T,Ha Y S.A study of fatigue crack propagation at a web stiffener on a longitudinal stiffener[J].J Mar Sci Technol,2010,15:176-189.

    [10]He Wentao,Liu Jingxi,Xie De.Numerical study on fatigue crack growth at a web-stiffener of ship structural details by an objected-oriented approach in conjunction with ABAQUS[J].Marine Structures,2014,35:45-69.

    [11]解 德,錢 勤,李長安.斷裂力學(xué)中的數(shù)值計(jì)算方法及工程應(yīng)用[M].北京:科學(xué)出版社,2009.

    [12]程 靳,趙樹山.斷裂力學(xué)[M].北京:科學(xué)出版社,2006.

    Life assessment of fatigue crack growth of typical details in hull longitudinals

    HE Wen-tao1,LIU Jing-xi1,2,XIE De1,2
    (1.School of Naval Architecture and Ocean Engineering,Huazhong University of Science and Technology,Wuhan 430074, China;2.Hubei Key Laboratory of Naval Architecture&Ocean Engineering Hydrodynamics,Wuhan 430074,China)

    A program(FCG-System)is developed to simulate fatigue crack growth combined with Virtual crack closure technique,crack growth criterion and substructure technique,based on ABAQUS using python scripting language.The fatigue crack growth in common longitudinal connection of oil tanker is simulated under lateral pressure and the effect of soft toe,back bracket and tripping bracket on crack growth path and life is discussed.Results indicate that soft toe,back brackets or tripping brackets will cause the curvature of crack growth path increasing and soft toe,tripping brackets can make crack growth life increase,and the back bracket can reduce crack growth life.

    crack growth path;crack growth life;virtual crack closure technique; substructure technique;FCG-system

    U661.43

    A

    10.3969/j.issn.1007-7294.2016.07.014

    1007-7294(2016)11-1475-10

    2016-05-13

    國家自然科學(xué)基金項(xiàng)目(NO.51609089);中國博士后科學(xué)基金項(xiàng)目(No.2016M592338)

    何文濤(1986-),男,博士研究生,E-mail:hewtsd@163.com;解 德(1964-),男,教授,博士生導(dǎo)師,E-mail:dexie@hust.edu.cn。

    猜你喜歡
    縱骨尖端腹板
    腹板開口對(duì)復(fù)合材料梁腹板剪切承載性能的影響
    無頂?shù)什坌螜M艙壁附近甲板構(gòu)件布置及型式優(yōu)化
    船海工程(2021年5期)2021-10-25 08:50:34
    變截面波形鋼腹板組合箱梁的剪應(yīng)力計(jì)算分析
    科學(xué)中國人(2018年8期)2018-07-23 02:26:56
    鋼箱梁超高腹板設(shè)計(jì)理論與方法
    上海公路(2018年3期)2018-03-21 05:55:50
    超大型集裝箱船縱骨貫穿孔的形狀優(yōu)化
    江蘇船舶(2017年4期)2017-10-12 01:32:46
    承受擱置腳載荷的艦船甲板縱骨設(shè)計(jì)方法
    基于艙段模型的大開口甲板結(jié)構(gòu)穩(wěn)定性分析與設(shè)計(jì)
    鏡頭看展
    基于位移相關(guān)法的重復(fù)壓裂裂縫尖端應(yīng)力場研究
    斷塊油氣田(2014年5期)2014-03-11 15:33:49
    看免费成人av毛片| 极品教师在线视频| 国产69精品久久久久777片| 九色国产91popny在线| 久久久久久久久久成人| 亚洲av中文av极速乱 | 欧美性感艳星| 国产中年淑女户外野战色| 国产久久久一区二区三区| 国产成人一区二区在线| 老熟妇仑乱视频hdxx| 国产不卡一卡二| 波多野结衣高清作品| 日本 欧美在线| 麻豆av噜噜一区二区三区| av专区在线播放| 全区人妻精品视频| 丰满人妻一区二区三区视频av| 欧美性猛交╳xxx乱大交人| 久久久国产成人免费| 亚洲欧美精品综合久久99| 在线国产一区二区在线| 成人av一区二区三区在线看| 国产精品人妻久久久影院| 99热网站在线观看| 亚洲狠狠婷婷综合久久图片| 18禁在线播放成人免费| 最好的美女福利视频网| 男女啪啪激烈高潮av片| 精品久久久久久久久亚洲 | 91在线精品国自产拍蜜月| 18禁黄网站禁片午夜丰满| 别揉我奶头 嗯啊视频| 不卡视频在线观看欧美| 五月伊人婷婷丁香| 在线观看美女被高潮喷水网站| 最近在线观看免费完整版| 国产大屁股一区二区在线视频| 国产免费男女视频| 亚洲成人精品中文字幕电影| 一夜夜www| 国产成人av教育| 国产精品一区二区免费欧美| 99久久九九国产精品国产免费| av天堂在线播放| 国产成年人精品一区二区| 国产精品久久视频播放| 国内久久婷婷六月综合欲色啪| 成人永久免费在线观看视频| 久久欧美精品欧美久久欧美| 我要搜黄色片| 精品日产1卡2卡| 欧美激情国产日韩精品一区| 2021天堂中文幕一二区在线观| 国产精品国产三级国产av玫瑰| 91久久精品国产一区二区成人| 亚洲欧美激情综合另类| 尾随美女入室| 永久网站在线| 国产精品1区2区在线观看.| 听说在线观看完整版免费高清| 此物有八面人人有两片| 我的老师免费观看完整版| 亚洲电影在线观看av| 亚洲一区二区三区色噜噜| 欧美日本亚洲视频在线播放| or卡值多少钱| 中文字幕精品亚洲无线码一区| 国产成人a区在线观看| 国产视频一区二区在线看| eeuss影院久久| 最近中文字幕高清免费大全6 | 午夜视频国产福利| 大型黄色视频在线免费观看| 国产一区二区激情短视频| 国产亚洲精品久久久com| 亚洲天堂国产精品一区在线| 国产女主播在线喷水免费视频网站 | av天堂中文字幕网| 亚洲国产精品成人综合色| 两人在一起打扑克的视频| 免费在线观看日本一区| 22中文网久久字幕| 国产一区二区亚洲精品在线观看| 日本撒尿小便嘘嘘汇集6| 99精品在免费线老司机午夜| 91久久精品电影网| 少妇裸体淫交视频免费看高清| 亚洲人成伊人成综合网2020| 国产亚洲精品av在线| 午夜免费成人在线视频| 欧美日韩乱码在线| 欧美另类亚洲清纯唯美| 婷婷精品国产亚洲av| 久久亚洲精品不卡| 在线观看午夜福利视频| 亚洲在线自拍视频| 人妻夜夜爽99麻豆av| 长腿黑丝高跟| 亚洲人与动物交配视频| 亚洲三级黄色毛片| 日韩欧美国产一区二区入口| 一夜夜www| 久久久国产成人精品二区| 内射极品少妇av片p| eeuss影院久久| 欧美黑人巨大hd| 国内久久婷婷六月综合欲色啪| 亚洲经典国产精华液单| 久久久久久久久久黄片| 中文字幕精品亚洲无线码一区| 午夜视频国产福利| 精华霜和精华液先用哪个| 欧美zozozo另类| 日本三级黄在线观看| 日本成人三级电影网站| 男女边吃奶边做爰视频| 女的被弄到高潮叫床怎么办 | 久久午夜福利片| 中文在线观看免费www的网站| 在线观看66精品国产| 久久亚洲精品不卡| 性插视频无遮挡在线免费观看| 精品欧美国产一区二区三| 春色校园在线视频观看| 成人av一区二区三区在线看| 99久久九九国产精品国产免费| 欧美一区二区精品小视频在线| 3wmmmm亚洲av在线观看| 国产欧美日韩精品亚洲av| 欧美性感艳星| 在现免费观看毛片| 一级av片app| 亚洲人成网站高清观看| 亚洲色图av天堂| 亚洲成人久久性| 久久热精品热| 大型黄色视频在线免费观看| 日韩精品有码人妻一区| 午夜福利视频1000在线观看| 俺也久久电影网| 国模一区二区三区四区视频| 久久精品久久久久久噜噜老黄 | 欧美丝袜亚洲另类 | 亚洲电影在线观看av| 国产欧美日韩精品亚洲av| 成人性生交大片免费视频hd| 国产精品乱码一区二三区的特点| 日本 av在线| 最后的刺客免费高清国语| 乱码一卡2卡4卡精品| 老师上课跳d突然被开到最大视频| 亚洲精品456在线播放app | 久久久久免费精品人妻一区二区| 尤物成人国产欧美一区二区三区| 久久久久久久亚洲中文字幕| 亚洲,欧美,日韩| 国产精品1区2区在线观看.| 国产精品久久久久久亚洲av鲁大| 婷婷色综合大香蕉| 天堂√8在线中文| 国产精品三级大全| 99riav亚洲国产免费| 在线看三级毛片| 久久久久久久亚洲中文字幕| 日本熟妇午夜| 夜夜看夜夜爽夜夜摸| 国产亚洲精品久久久久久毛片| 精品久久久噜噜| 亚洲av二区三区四区| 国产男靠女视频免费网站| 久久精品国产亚洲av香蕉五月| 色哟哟·www| 国产高清有码在线观看视频| 一本一本综合久久| 日韩欧美在线二视频| 搡老熟女国产l中国老女人| 欧美最新免费一区二区三区| 亚洲专区中文字幕在线| 成人三级黄色视频| 色哟哟·www| 久久国内精品自在自线图片| 乱系列少妇在线播放| 亚洲四区av| av女优亚洲男人天堂| 国产激情偷乱视频一区二区| 精品国产三级普通话版| 日韩 亚洲 欧美在线| 亚洲欧美日韩东京热| 91麻豆精品激情在线观看国产| a在线观看视频网站| 国产av一区在线观看免费| 一区二区三区激情视频| 日日啪夜夜撸| 少妇丰满av| 99久久中文字幕三级久久日本| 久久人人爽人人爽人人片va| 欧美色欧美亚洲另类二区| 99久久九九国产精品国产免费| 精品一区二区免费观看| 韩国av在线不卡| 国产黄片美女视频| 91久久精品国产一区二区三区| 最新在线观看一区二区三区| 麻豆成人av在线观看| 亚洲一级一片aⅴ在线观看| 国产成人aa在线观看| 国产精品国产高清国产av| 国产亚洲91精品色在线| 九九爱精品视频在线观看| 成人av一区二区三区在线看| 男女之事视频高清在线观看| 久久久久久久久久久丰满 | or卡值多少钱| 午夜精品一区二区三区免费看| 欧美又色又爽又黄视频| 天天躁日日操中文字幕| 99久久精品一区二区三区| www日本黄色视频网| 成年人黄色毛片网站| 波多野结衣巨乳人妻| 久久精品夜夜夜夜夜久久蜜豆| 久久久久九九精品影院| 伊人久久精品亚洲午夜| 成年女人永久免费观看视频| 国产精品伦人一区二区| 日本免费一区二区三区高清不卡| 亚洲av熟女| 观看免费一级毛片| 少妇人妻一区二区三区视频| 欧美色视频一区免费| 日韩在线高清观看一区二区三区 | 狂野欧美白嫩少妇大欣赏| 国产黄片美女视频| 少妇裸体淫交视频免费看高清| 99久久成人亚洲精品观看| 神马国产精品三级电影在线观看| 成人无遮挡网站| 在线播放无遮挡| 亚洲一级一片aⅴ在线观看| 美女高潮的动态| 国产午夜福利久久久久久| 嫩草影院新地址| 男女那种视频在线观看| 最新中文字幕久久久久| 女人被狂操c到高潮| 99视频精品全部免费 在线| 国产精品精品国产色婷婷| 国产又黄又爽又无遮挡在线| 精品久久久久久,| 春色校园在线视频观看| 亚洲av二区三区四区| 久久久精品欧美日韩精品| 久久国产精品人妻蜜桃| 97人妻精品一区二区三区麻豆| 禁无遮挡网站| 免费无遮挡裸体视频| 国产av不卡久久| 欧洲精品卡2卡3卡4卡5卡区| 人人妻,人人澡人人爽秒播| 国内毛片毛片毛片毛片毛片| 69人妻影院| avwww免费| 欧美日本亚洲视频在线播放| 亚洲av.av天堂| 午夜精品久久久久久毛片777| 久久精品国产鲁丝片午夜精品 | 三级国产精品欧美在线观看| 十八禁网站免费在线| 亚洲最大成人手机在线| 麻豆精品久久久久久蜜桃| 国产伦一二天堂av在线观看| 国产不卡一卡二| 69人妻影院| 一级黄片播放器| 欧美精品啪啪一区二区三区| 国产精品av视频在线免费观看| 两个人视频免费观看高清| 亚洲精品乱码久久久v下载方式| 又黄又爽又免费观看的视频| 97人妻精品一区二区三区麻豆| 偷拍熟女少妇极品色| 成人鲁丝片一二三区免费| 成年版毛片免费区| 婷婷六月久久综合丁香| 在现免费观看毛片| 国内精品久久久久精免费| 欧美不卡视频在线免费观看| 男人舔女人下体高潮全视频| 欧美xxxx黑人xx丫x性爽| 国产高清三级在线| 日韩欧美 国产精品| 99热网站在线观看| 国产成人影院久久av| 国产亚洲欧美98| 美女被艹到高潮喷水动态| 亚洲中文字幕日韩| 亚洲精品成人久久久久久| 我要搜黄色片| 一级黄片播放器| 国产在线男女| 久久热精品热| 有码 亚洲区| 三级毛片av免费| 波多野结衣巨乳人妻| 男女啪啪激烈高潮av片| 亚洲最大成人av| 婷婷丁香在线五月| 欧美潮喷喷水| 亚洲成人免费电影在线观看| 亚洲精品粉嫩美女一区| 欧美3d第一页| 亚洲无线观看免费| 在线天堂最新版资源| 亚洲综合色惰| 18禁裸乳无遮挡免费网站照片| 尤物成人国产欧美一区二区三区| 91久久精品电影网| 久久久国产成人免费| h日本视频在线播放| 午夜亚洲福利在线播放| 国产精品不卡视频一区二区| 乱人视频在线观看| 18+在线观看网站| 亚洲内射少妇av| 亚洲av成人精品一区久久| www.色视频.com| 人妻久久中文字幕网| 91午夜精品亚洲一区二区三区 | 精品福利观看| 亚洲av一区综合| 亚洲成av人片在线播放无| 亚洲精品粉嫩美女一区| 人妻丰满熟妇av一区二区三区| 成人国产麻豆网| 国产v大片淫在线免费观看| 婷婷精品国产亚洲av在线| 简卡轻食公司| 亚洲专区中文字幕在线| 久久午夜亚洲精品久久| 成年版毛片免费区| 两个人视频免费观看高清| 国产精品99久久久久久久久| 一级毛片久久久久久久久女| 一个人看视频在线观看www免费| 欧美日韩中文字幕国产精品一区二区三区| 成人综合一区亚洲| 亚洲内射少妇av| 联通29元200g的流量卡| 亚洲国产欧美人成| 成人精品一区二区免费| 国产一区二区在线av高清观看| 亚洲天堂国产精品一区在线| 国产精品国产三级国产av玫瑰| 国产精品久久电影中文字幕| 国产免费av片在线观看野外av| 国产精品久久久久久久久免| 色噜噜av男人的天堂激情| 极品教师在线视频| 黄片wwwwww| 特大巨黑吊av在线直播| eeuss影院久久| 精品午夜福利视频在线观看一区| 亚洲自拍偷在线| 欧美丝袜亚洲另类 | 国产精品久久久久久久久免| 少妇被粗大猛烈的视频| 亚洲图色成人| 国产91精品成人一区二区三区| 精品久久久久久成人av| 色综合亚洲欧美另类图片| 久久草成人影院| 国产男靠女视频免费网站| 可以在线观看的亚洲视频| 午夜福利高清视频| 51国产日韩欧美| 久久热精品热| 亚洲成a人片在线一区二区| 男女啪啪激烈高潮av片| xxxwww97欧美| 一级a爱片免费观看的视频| aaaaa片日本免费| 亚洲专区中文字幕在线| 国产成人aa在线观看| 日韩人妻高清精品专区| 国产精品美女特级片免费视频播放器| 22中文网久久字幕| 国产日本99.免费观看| 亚洲精品在线观看二区| 国产视频内射| 国产一区二区在线av高清观看| aaaaa片日本免费| 欧美+亚洲+日韩+国产| 1024手机看黄色片| 97碰自拍视频| 成人鲁丝片一二三区免费| 日韩欧美精品v在线| 国产亚洲欧美98| 国产av在哪里看| 麻豆av噜噜一区二区三区| 精品久久久久久,| 桃色一区二区三区在线观看| 国产三级在线视频| 亚洲av电影不卡..在线观看| 日日撸夜夜添| 人妻丰满熟妇av一区二区三区| 亚洲av美国av| 亚洲人成网站高清观看| 成年女人永久免费观看视频| 久久久久免费精品人妻一区二区| 日韩欧美免费精品| 日本熟妇午夜| 97热精品久久久久久| 亚洲国产高清在线一区二区三| www.www免费av| 成人特级黄色片久久久久久久| 久久久久国内视频| 成年女人看的毛片在线观看| 午夜爱爱视频在线播放| 国产精品福利在线免费观看| 欧美成人免费av一区二区三区| 18禁黄网站禁片免费观看直播| 又爽又黄无遮挡网站| av中文乱码字幕在线| 男人的好看免费观看在线视频| 91狼人影院| 亚洲四区av| 熟妇人妻久久中文字幕3abv| 亚洲内射少妇av| 午夜久久久久精精品| 亚洲电影在线观看av| 亚洲av熟女| 老女人水多毛片| 人人妻人人澡欧美一区二区| 最后的刺客免费高清国语| 免费av观看视频| 婷婷六月久久综合丁香| 国产69精品久久久久777片| 我的老师免费观看完整版| 俄罗斯特黄特色一大片| 亚洲久久久久久中文字幕| 麻豆久久精品国产亚洲av| 欧美成人a在线观看| 不卡视频在线观看欧美| 亚洲无线在线观看| 精品久久久噜噜| 美女黄网站色视频| 在线播放国产精品三级| 2021天堂中文幕一二区在线观| 成年人黄色毛片网站| 久久久久国产精品人妻aⅴ院| 免费不卡的大黄色大毛片视频在线观看 | 日日啪夜夜撸| bbb黄色大片| 国内精品美女久久久久久| 久久国内精品自在自线图片| 日本黄色视频三级网站网址| 在线免费观看的www视频| 国产在视频线在精品| 久久精品影院6| 可以在线观看毛片的网站| 变态另类丝袜制服| 真实男女啪啪啪动态图| 国产亚洲精品综合一区在线观看| 国产精品国产高清国产av| 男女边吃奶边做爰视频| 老女人水多毛片| 婷婷六月久久综合丁香| 精品99又大又爽又粗少妇毛片 | 久久这里只有精品中国| 两人在一起打扑克的视频| 欧美日韩综合久久久久久 | 午夜激情福利司机影院| 悠悠久久av| 日韩亚洲欧美综合| 九九在线视频观看精品| 丝袜美腿在线中文| 国产成人一区二区在线| 国产一区二区三区在线臀色熟女| 亚洲在线自拍视频| 搡女人真爽免费视频火全软件 | 一本久久中文字幕| 一级黄色大片毛片| 亚洲av第一区精品v没综合| 赤兔流量卡办理| 欧美日韩中文字幕国产精品一区二区三区| 淫妇啪啪啪对白视频| 丰满的人妻完整版| 精品午夜福利视频在线观看一区| 97碰自拍视频| 可以在线观看的亚洲视频| 久久热精品热| 国产麻豆成人av免费视频| 日韩一本色道免费dvd| 国产真实乱freesex| av在线蜜桃| 在线免费十八禁| 精品人妻偷拍中文字幕| 18禁黄网站禁片午夜丰满| 国产免费一级a男人的天堂| 简卡轻食公司| 亚洲专区中文字幕在线| av女优亚洲男人天堂| 身体一侧抽搐| 成年人黄色毛片网站| 国产亚洲精品久久久久久毛片| 免费搜索国产男女视频| 国产午夜精品论理片| 狠狠狠狠99中文字幕| 欧美激情久久久久久爽电影| 欧美xxxx黑人xx丫x性爽| 亚洲第一区二区三区不卡| 亚洲国产色片| 一本精品99久久精品77| 日韩欧美精品v在线| 成人亚洲精品av一区二区| 日韩欧美一区二区三区在线观看| 美女高潮的动态| 精品人妻视频免费看| 亚洲成a人片在线一区二区| 日本与韩国留学比较| 欧美日韩亚洲国产一区二区在线观看| 深夜a级毛片| 女的被弄到高潮叫床怎么办 | 日韩亚洲欧美综合| 国产三级中文精品| 久99久视频精品免费| 午夜福利高清视频| 麻豆av噜噜一区二区三区| 欧美日韩精品成人综合77777| 日本黄色片子视频| 亚洲一区二区三区色噜噜| 丝袜美腿在线中文| 在线天堂最新版资源| 熟女人妻精品中文字幕| 波多野结衣高清无吗| 麻豆成人午夜福利视频| 亚洲精品亚洲一区二区| 极品教师在线视频| 香蕉av资源在线| 欧美最黄视频在线播放免费| 亚洲性夜色夜夜综合| 欧美日韩国产亚洲二区| 夜夜爽天天搞| 国产乱人伦免费视频| 欧美激情在线99| 白带黄色成豆腐渣| 欧美成人a在线观看| 看十八女毛片水多多多| 欧美xxxx黑人xx丫x性爽| 亚洲成人免费电影在线观看| 国产成人影院久久av| 精品不卡国产一区二区三区| 高清毛片免费观看视频网站| 亚洲一区二区三区色噜噜| 可以在线观看的亚洲视频| av天堂在线播放| 国产日本99.免费观看| 日本色播在线视频| 日本三级黄在线观看| www.色视频.com| 亚洲人成网站高清观看| 亚洲精品粉嫩美女一区| 可以在线观看毛片的网站| 中文字幕熟女人妻在线| 伦精品一区二区三区| 免费人成视频x8x8入口观看| 人妻制服诱惑在线中文字幕| 日本 欧美在线| 动漫黄色视频在线观看| 99久久精品国产国产毛片| 国产精品综合久久久久久久免费| 麻豆成人午夜福利视频| 日日啪夜夜撸| 老司机福利观看| 亚洲综合色惰| 精品久久久久久久人妻蜜臀av| 国内毛片毛片毛片毛片毛片| 久久久久久久久久成人| 亚洲欧美日韩卡通动漫| 亚洲成人久久性| 欧美国产日韩亚洲一区| 偷拍熟女少妇极品色| av在线蜜桃| 亚洲国产精品久久男人天堂| 极品教师在线免费播放| 熟女电影av网| 国产v大片淫在线免费观看| 麻豆国产97在线/欧美| 欧美最新免费一区二区三区| 最新在线观看一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 极品教师在线免费播放| 亚洲欧美日韩高清在线视频| 亚洲18禁久久av| 国产精品野战在线观看| 欧美性感艳星| 99久久九九国产精品国产免费| 亚洲黑人精品在线| 一个人看的www免费观看视频| 丝袜美腿在线中文| 在线天堂最新版资源| 国产美女午夜福利| 亚洲国产精品sss在线观看| 黄色丝袜av网址大全| 国产精品人妻久久久影院| 少妇的逼水好多| 岛国在线免费视频观看| 又紧又爽又黄一区二区| 男人的好看免费观看在线视频| av在线蜜桃| 免费看av在线观看网站| 免费av毛片视频| 小蜜桃在线观看免费完整版高清| 日韩欧美免费精品|