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

    雙面焊U肋內(nèi)側(cè)焊縫初始裂紋特性研究

    2021-11-25 11:12:28郭殊倫雷俊卿黃祖慰
    關(guān)鍵詞:裂紋

    郭殊倫,雷俊卿,黃祖慰

    (1.北京交通大學(xué)土木建筑工程學(xué)院,北京,100044;2.交通運(yùn)輸部路網(wǎng)監(jiān)測(cè)與應(yīng)急處置中心,北京,100029)

    正交異性鋼橋面板在底部設(shè)置了縱、橫向正交垂直的加勁肋,可滿足2個(gè)方向上不同的剛度需求,具有強(qiáng)度高、質(zhì)量小、跨越能力強(qiáng)、制造容易和施工便利等優(yōu)點(diǎn),在工程中得到了普遍應(yīng)用。但在車輛荷載的反復(fù)作用下,其疲勞開裂問題突出[1-6],對(duì)于具有閉口縱肋的構(gòu)造來說,連接縱向U肋與頂板的焊縫是正交異性鋼橋面板中最易發(fā)生疲勞破壞的部位之一[7-12]。

    焊根位置未熔透是導(dǎo)致傳統(tǒng)單面焊構(gòu)造細(xì)節(jié)疲勞開裂的主因,其主導(dǎo)疲勞失效模式為頂板焊根開裂并沿頂板厚度方向擴(kuò)展[13]。針對(duì)該問題,國(guó)內(nèi)外研究者研發(fā)了正交異性鋼橋面板自動(dòng)化U 肋內(nèi)焊系統(tǒng),推出U 肋與頂板新型雙面焊構(gòu)造以提升該細(xì)節(jié)的疲勞性能,并成功運(yùn)用于實(shí)際工程中[14-16]。然而,目前對(duì)于雙面焊構(gòu)造的疲勞開裂問題尚缺少深入、系統(tǒng)分析,與單面焊相比,雙面焊構(gòu)造疲勞破壞模式上最大的不同在于:當(dāng)熔透率滿足一定要求時(shí)可基本消除焊根開裂,但存在U肋內(nèi)側(cè)焊趾開裂并沿腹板厚度發(fā)展的疲勞破壞模式[17],這是單面焊構(gòu)造所不具備的。這種裂紋位于橋面板下方閉口肋內(nèi)部,常規(guī)巡檢難以檢出,檢出時(shí)裂紋一般穿透鋼板,可能引發(fā)一系列的橋梁耐久性問題,進(jìn)而影響運(yùn)營(yíng)質(zhì)量[18]。

    斷裂力學(xué)的發(fā)展為疲勞問題的研究提供了理論基礎(chǔ)和有效工具,國(guó)內(nèi)外學(xué)者基于斷裂力學(xué)理論對(duì)正交異性鋼橋面板U 肋與頂板焊接構(gòu)造的疲勞裂紋擴(kuò)展特性開展了大量研究,但其中鮮有對(duì)于初始裂紋角度問題的分析。由于內(nèi)焊過程的復(fù)雜性和隨機(jī)性,焊接裂紋等初始缺陷難以避免,U肋內(nèi)側(cè)焊趾的初始裂紋也可能具有較強(qiáng)的幾何不確定性,且目前對(duì)于新型雙面焊構(gòu)造下的U 肋內(nèi)焊趾疲勞開裂問題認(rèn)識(shí)不足,亟須對(duì)該細(xì)節(jié)的疲勞裂紋擴(kuò)展特性進(jìn)行研究。

    針對(duì)上述問題,本文作者基于三維斷裂力學(xué)理論,結(jié)合正交異性鋼橋面板足尺節(jié)段疲勞試驗(yàn),研究U 肋與頂板雙面焊構(gòu)造的U 肋內(nèi)側(cè)焊趾疲勞裂紋擴(kuò)展特性,重點(diǎn)研究初始裂紋角度因素的影響。研究成果有助于研究者們更加全面地認(rèn)識(shí)該細(xì)節(jié)處的疲勞裂紋擴(kuò)展問題,并為相關(guān)研究提供依據(jù)和參考。

    1 疲勞裂紋擴(kuò)展理論及方法

    1.1 表面裂紋特征化

    圖1所示為典型的半橢圓表面裂紋。國(guó)際焊接學(xué)會(huì)(IIW)指出,當(dāng)不明確面狀表面缺陷具體形狀尺寸時(shí),可采用半橢圓形對(duì)其進(jìn)行模擬[19],裂紋深度為a,表面長(zhǎng)度為2c。

    圖1 半橢圓形表面裂紋Fig.1 Semi-elliptical surface crack

    1.2 基于斷裂力學(xué)的疲勞裂紋擴(kuò)展理論

    斷裂力學(xué)中,KNOWLES 等[20]提出的M 積分比J積分能更有效地反映裂紋擴(kuò)展規(guī)律,本文選取M 積分計(jì)算應(yīng)力強(qiáng)度因子(stress intensity factor,SIF)。在此基礎(chǔ)上采用BS7910 所推薦的有效應(yīng)力強(qiáng)度因子幅值計(jì)算公式來體現(xiàn)3種應(yīng)力強(qiáng)度因子的綜合影響[21],并采用應(yīng)用較廣泛的最大周向拉應(yīng)力理論計(jì)算裂紋面擴(kuò)展角度,表達(dá)式如下:

    式中:ν為泊松比,此處取0.3;ΔKI,ΔKII和ΔKIII分別為I,II 和III 型裂紋對(duì)應(yīng)的應(yīng)力強(qiáng)度因子幅值。

    采用NASGRO擴(kuò)展模型計(jì)算裂紋的擴(kuò)展速率,NASGRO基于FORMAN模型進(jìn)行了改進(jìn),考慮裂紋張開效應(yīng)、應(yīng)力比、擴(kuò)展門檻值及斷裂韌性的影響,因而比PARIS 和WALKER 模型更加全面,其表達(dá)式如下[22]:

    式中:R為應(yīng)力比;ΔK為應(yīng)力強(qiáng)度因子幅值;C,n,p和q為與材料相關(guān)的參數(shù);ΔKth為應(yīng)力強(qiáng)度因子幅值門檻值;Kc為材料斷裂韌性;Kmax為最大應(yīng)力強(qiáng)度因子;f為裂紋張開函數(shù),具體形式可見文獻(xiàn)[23]。以上參數(shù)可根據(jù)相關(guān)數(shù)據(jù)庫(kù)進(jìn)行選取和計(jì)算,本文選取美國(guó)材料與試驗(yàn)協(xié)會(huì)(ASTM)數(shù)據(jù)庫(kù)結(jié)果。

    裂紋擴(kuò)展采用指定應(yīng)力強(qiáng)度因子中值節(jié)點(diǎn)擴(kuò)展步長(zhǎng)的方式[23],其原理如圖2所示。指定中值節(jié)點(diǎn)i處的擴(kuò)展可由下式計(jì)算:

    圖2 指定中值點(diǎn)擴(kuò)展Fig.2 Specified median extension

    式中:Δa為擴(kuò)展步長(zhǎng);da/dN為擴(kuò)展速率;下標(biāo)i和median分別表示裂紋前緣第i點(diǎn)及應(yīng)力強(qiáng)度因子中值點(diǎn)。

    1.3 FRANC3D-ANSYS交互技術(shù)

    裂紋前緣位置更新及網(wǎng)格劃分可運(yùn)用ANSYS或ABAQUS 等有限元軟件實(shí)現(xiàn),但過程較為煩瑣且難以控制網(wǎng)格質(zhì)量。FRANC3D軟件可簡(jiǎn)單、高效地進(jìn)行裂紋網(wǎng)格劃分,并已在航空、機(jī)械等領(lǐng)域中得到廣泛應(yīng)用。本文運(yùn)用FRANC3D-ANSYS交互技術(shù)對(duì)疲勞裂紋的擴(kuò)展過程進(jìn)行模擬,交互流程如圖3所示。FRANC3D 能夠方便地對(duì)原模型進(jìn)行切分并導(dǎo)入ANSYS進(jìn)行合并計(jì)算,從而顯著減小工作量,可在裂紋前緣生成高質(zhì)量的對(duì)稱單元環(huán)形式網(wǎng)格(template)以盡量消除計(jì)算誤差,并對(duì)結(jié)果進(jìn)行光順化處理,顯著提升了計(jì)算精度。

    圖3 FRANC3D-ANSYS交互流程Fig.3 Interaction process of FRANC3D-ANSYS

    2 疲勞試驗(yàn)及數(shù)值模擬方法驗(yàn)證

    2.1 試驗(yàn)概況

    試件共3 個(gè)(編號(hào)為T1,T2 和T3),采用Q345qE 鋼材,其為某大橋正交異性鋼橋面板中的雙面焊U 肋節(jié)段,截取頂板寬度為700 mm,長(zhǎng)度為300 mm,U 肋與頂板采用雙面焊全熔透焊接方式以徹底消除焊根開裂的影響,細(xì)部尺寸、約束加載裝置、焊縫形式及試驗(yàn)現(xiàn)場(chǎng)布置如圖4所示。加載方式及設(shè)備參照曾志斌等[24]的試驗(yàn)進(jìn)行安排,作動(dòng)器荷載通過分載梁傳遞至U 肋外側(cè)頂板(圖4(b))而使U肋內(nèi)側(cè)焊趾附近受拉彎應(yīng)力,以考察其疲勞開裂性能。根據(jù)靜載試驗(yàn)結(jié)果取65 kN作為疲勞荷載最大值,應(yīng)力比為0.031,頻率為4 Hz。

    圖4 試件及試驗(yàn)布置Fig.4 Specimen and arrangement of test

    2.2 疲勞破壞現(xiàn)象

    裂紋位置及擴(kuò)展情況如圖5(a)所示。試件T1,T2 和T3 均于U 肋內(nèi)側(cè)焊趾外推0~2 mm 范圍內(nèi)起裂,裂紋沿腹板厚度方向擴(kuò)展,最終導(dǎo)致U 肋與頂板的連接失效。文獻(xiàn)[25]對(duì)9 個(gè)雙面焊U 肋試件進(jìn)行疲勞試驗(yàn)也得到了與本文相似的結(jié)果。試件整體斷面與腹板表面基本垂直,但初始裂紋面區(qū)域與腹板表面的夾角具有不確定性;初始裂紋的形狀近似為半橢圓或半圓。通過實(shí)測(cè)并計(jì)算初始裂紋區(qū)與腹板表面法向的夾角,規(guī)定垂直腹板表面為0°,偏向頂板方向?yàn)檎?,反之為?fù),測(cè)得初始裂紋區(qū)與腹板表面法向夾角的最大值約為+16°,最小值約為-10°,分別如圖5(b)和5(c)所示。

    圖5 試件的疲勞破壞現(xiàn)象Fig.5 Fatigue failures of specimens

    2.3 疲勞裂紋模擬方法驗(yàn)證

    2.3.1 標(biāo)準(zhǔn)試樣驗(yàn)證

    試驗(yàn)得到的裂紋均為表面裂紋,以梁長(zhǎng)錦等[26]選用的含表面裂紋的標(biāo)準(zhǔn)拉伸鋼板試樣作為驗(yàn)證算例。對(duì)文獻(xiàn)[26]中2 組試樣的表面裂紋應(yīng)力強(qiáng)度因子進(jìn)行計(jì)算,結(jié)果如圖6所示。從圖6可以看出運(yùn)用FRANC3D-ANSYS 交互技術(shù)可以合理計(jì)算表面裂紋的應(yīng)力強(qiáng)度因子。

    圖6 應(yīng)力強(qiáng)度因子對(duì)比Fig.6 Comparison of SIF

    2.3.2 本文試驗(yàn)驗(yàn)證

    裂紋的擴(kuò)展路徑及形狀變化對(duì)比結(jié)果如圖7和圖8所示。從圖7和圖8可以看出FRANC3DANSYS 交互技術(shù)能夠合理模擬裂紋擴(kuò)展路徑及形貌變化特征。

    圖7 裂紋擴(kuò)展路徑對(duì)比Fig.7 Comparison of crack growth path

    圖8 裂紋擴(kuò)展形狀對(duì)比Fig.8 Comparison of crack shape in growth

    綜上,可以借助FRANC3D-ANSYS 交互技術(shù)對(duì)裂紋擴(kuò)展特性進(jìn)行分析。

    3 裂紋擴(kuò)展數(shù)值模擬

    3.1 計(jì)算模型

    運(yùn)用ANSYS 按照試件實(shí)際尺寸建模,統(tǒng)一采用六面體網(wǎng)格劃分,單元類型為Solid185。圖9所示為在FRANC3D中切分子模型,于U肋內(nèi)焊趾中點(diǎn)垂直外推1 mm處插入初始裂紋,采用三層對(duì)稱單元環(huán)的形式劃分裂紋前緣網(wǎng)格以保證計(jì)算精度,單元環(huán)最內(nèi)圈采用15 節(jié)點(diǎn)楔形體單元,中間層采用20 節(jié)點(diǎn)六面體單元,裂尖與全局模型過渡區(qū)采用13 節(jié)點(diǎn)金字塔單元,單元環(huán)外部全局模型采用10節(jié)點(diǎn)四面體單元,初始狀態(tài)下單元數(shù)為166 532個(gè)。裂紋插入位置及前緣網(wǎng)格劃分如圖10所示。

    圖9 子模型切分Fig.9 Local model dividing

    圖10 裂紋位置及網(wǎng)格Fig.10 Position and mesh of crack

    3.2 初始裂紋參數(shù)

    根據(jù)試驗(yàn)結(jié)果及文獻(xiàn)[27],半橢圓表面初始裂紋深度一般不超過0.5 mm,可取該值作為初始裂紋深度以反映最不利情況。初始裂紋幾何參數(shù)如表1所示,其中,ω為初始裂紋形狀比,ω=a0/(2c0),α0為初始角度,其相關(guān)定義如圖11所示。

    圖11 初始裂紋角度定義Fig.11 Definition of initial crack angle

    表1 初始裂紋參數(shù)Table 1 Parameters of initial cracks

    3.3 裂紋擴(kuò)展

    計(jì)算模型如圖12所示。采用與試驗(yàn)相同的方式對(duì)U 肋下半部進(jìn)行固定約束,并在兩側(cè)頂板對(duì)稱垂直加載。最大拉應(yīng)力位于U 肋內(nèi)焊趾處,當(dāng)荷載為65 kN 時(shí),最大表面拉應(yīng)力模擬值為340.52 MPa,小于且接近屈服強(qiáng)度345 MPa,選取2~65 kN作為彈性階段最不利荷載幅進(jìn)行加載,應(yīng)力比R=0.031。初始裂紋前緣單元環(huán)網(wǎng)格半徑?。?.02~0.1)×min(a0,c0);裂紋擴(kuò)展時(shí),單元環(huán)網(wǎng)格半徑?。?.1~0.6)×min(Δai,Δci),擴(kuò)展步長(zhǎng)取裂紋最小特征尺寸的15%。

    圖12 試件的ANSYS有限元模型Fig.12 ANSYS FEM model of specimens

    4 裂紋擴(kuò)展特性分析

    4.1 裂紋前緣計(jì)算點(diǎn)選取

    裂紋前緣計(jì)算點(diǎn)如圖13所示,取裂紋前緣A,B,C,D,E和F點(diǎn)作為計(jì)算點(diǎn),分別對(duì)應(yīng)θ為0.5π,0.4π,0.3π,0.2π,0.1π 和0。其中,表面點(diǎn)及其附近較小區(qū)域不符合平面應(yīng)變假設(shè),計(jì)算結(jié)果存在一定的誤差,根據(jù)應(yīng)力強(qiáng)度因子手冊(cè)[28]的建議,可根據(jù)偏于內(nèi)部的計(jì)算值進(jìn)行小范圍外推得到較為準(zhǔn)確的結(jié)果。

    圖13 裂紋前緣計(jì)算點(diǎn)Fig.13 Calculation points on crack front

    4.2 基于初始角度的裂紋擴(kuò)展

    4.2.1 試驗(yàn)現(xiàn)象分析

    當(dāng)試件剩余斷面不足以承受外荷載(剩余斷面應(yīng)力達(dá)到或超過了材料的靜強(qiáng)度)時(shí),裂紋將出現(xiàn)失穩(wěn)擴(kuò)展,從而引發(fā)截面瞬斷[29]。試件斷口如圖14所示。從圖14可知:斷口存在明顯的瞬斷區(qū),其形貌較為粗糙,與疲勞條帶部分區(qū)別顯著,后者平均寬度占斷口總寬度的0.5~0.6。

    圖14 試件斷口Fig.14 Fracture section of specimen

    以半圓形初始裂紋C3 為代表,在不同的初始角度下對(duì)裂紋擴(kuò)展過程進(jìn)行計(jì)算。用歸一化裂紋深度D=at/t衡量裂紋的擴(kuò)展程度,其中at為裂紋前緣最深點(diǎn)在腹板橫截面上的投影深度,t為腹板厚度。為盡量避免裂紋擴(kuò)展過程中表面點(diǎn)計(jì)算誤差的干擾,選取裂紋前緣A點(diǎn)和E點(diǎn)進(jìn)行分析。有效應(yīng)力強(qiáng)度因子幅值ΔKeff的變化情況如圖15所示。從圖15可見:不同初始角度的裂紋應(yīng)力強(qiáng)度因子整體發(fā)展趨勢(shì)相近,當(dāng)D在0.3~0.4 之間時(shí)出現(xiàn)一定程度的離散,結(jié)合試件疲勞條帶形狀可認(rèn)為此時(shí)裂紋進(jìn)入失穩(wěn)擴(kuò)展;當(dāng)D在0.4~0.5之間時(shí),應(yīng)力強(qiáng)度因子的離散程度再次顯著增大,在D=0.5左右擴(kuò)展模擬終止。模擬結(jié)果與試驗(yàn)中試件的斷口特征基本一致。不同初始角度的裂紋達(dá)到失穩(wěn)時(shí)的歸一化裂深具有一致性,與初始角度無關(guān),進(jìn)一步表明裂紋失穩(wěn)及截面瞬斷取決于剩余截面抵抗外荷載的能力,且通過應(yīng)力強(qiáng)度因子的整體變化趨勢(shì)判斷,裂紋沿長(zhǎng)度方向的擴(kuò)展是后期導(dǎo)致截面剩余面積減小的主因。

    圖15 有效應(yīng)力強(qiáng)度因子幅值發(fā)展曲線Fig.15 Development curve of equivalent SIF amplitude

    4.2.2 裂紋擴(kuò)展階段的疲勞壽命分布

    這里僅以最深點(diǎn)A的擴(kuò)展速率為例進(jìn)行展示,如圖16所示。按速率變化情況可將裂紋最深點(diǎn)擴(kuò)展過程分為3 個(gè)階段:擴(kuò)展速率快速增長(zhǎng)階段(第一階段)、平穩(wěn)下降階段(第二階段,緩沖段)和失穩(wěn)階段(第三階段)。根據(jù)計(jì)算結(jié)果,統(tǒng)一將D=0.35 作為本文算例中裂紋達(dá)到失穩(wěn)時(shí)的歸一化裂深,此時(shí)所對(duì)應(yīng)的疲勞壽命記為擴(kuò)展階段的全壽命N;用Nu表示擴(kuò)展速率快速增長(zhǎng)階段的壽命,Ns表示速率平穩(wěn)下降階段壽命,則有N=Nu+Ns。裂紋最深點(diǎn)擴(kuò)展速率隨荷載循環(huán)次數(shù)的變化如表2所示。

    圖16 裂紋最深點(diǎn)擴(kuò)展速率Fig.16 Development curve of equivalent SIF amplitude

    由圖16和表2可知:在裂紋擴(kuò)展速率方面,α0對(duì)第一階段影響顯著,當(dāng)|α0|較大時(shí),裂紋初期擴(kuò)展速率較小,但速率增幅較大,該階段內(nèi)裂紋與外界的能量交換速度快,周邊應(yīng)力狀態(tài)不穩(wěn)定;第二階段裂紋擴(kuò)展速率較為平穩(wěn),幾乎不受α0影響。在疲勞壽命方面,|α0|越大,Nu/N就越大,相應(yīng)的Ns/N就越小,意味著|α0|較大的裂紋在整個(gè)擴(kuò)展過程中的不穩(wěn)定性更強(qiáng),緩沖段占比更小,裂紋經(jīng)歷速率快速增長(zhǎng)期后會(huì)相對(duì)更快地達(dá)到失穩(wěn),進(jìn)而導(dǎo)致截面瞬間破壞??梢婋p面焊U 肋內(nèi)側(cè)焊趾的初始裂紋角度主要影響構(gòu)造的疲勞壽命分布,初始角度較大的裂紋會(huì)對(duì)橋梁檢測(cè)周期的制定、維修加固等工作帶來不利影響。研究中應(yīng)全面考慮裂紋的初始角度因素,在進(jìn)行預(yù)制裂紋擴(kuò)展試驗(yàn)或數(shù)值模擬等之前,應(yīng)合理設(shè)置初始裂紋角度。

    表2 裂紋擴(kuò)展速率變化Table 2 Change of crack growth rate

    4.3 初始裂紋角度

    以半圓裂紋C3 為例,用|KII/KI|和|KIII/KI|反映3種應(yīng)力強(qiáng)度因子的大小關(guān)系,得到裂紋前緣在不同初始角度下的應(yīng)力強(qiáng)度因子變化規(guī)律如圖17所示。分析結(jié)果可知:

    圖17 應(yīng)力強(qiáng)度因子比值曲面及等高線圖Fig.17 SIF ratio curve surface and contour map

    1)裂紋整體上屬于I型主導(dǎo)的I-II-III型復(fù)合裂紋,在某一初始角度下可視為純I 型,該算例中α0為+ 5°左右,記該特征角度為α0-p;

    2)裂紋的|α0|越大,II和III型特征越明顯。當(dāng)|α0|較大時(shí),KII在表面和最深點(diǎn)均具有較大占比,此時(shí),裂紋前緣不同位置所處的應(yīng)力差異較大,但總體上仍為I 型主導(dǎo)的裂紋;KIII自裂紋最深點(diǎn)向表面逐漸增大。具備一定初始角度的深處和表面附近裂紋在擴(kuò)展過程中,分別在KII和KIII的主導(dǎo)作用下發(fā)生偏轉(zhuǎn),如圖18所示。

    圖18 裂紋偏轉(zhuǎn)Fig.18 Crack deflection

    裂紋的最初擴(kuò)展?fàn)顟B(tài)是指初始裂紋經(jīng)歷最初1次荷載循環(huán)前后的變化狀態(tài),根據(jù)裂紋擴(kuò)展門檻值ΔKth的定義,可通過對(duì)該狀態(tài)進(jìn)行計(jì)算來判斷裂紋在給定荷載幅下是否具備擴(kuò)展能力。根據(jù)ASTM數(shù)據(jù)庫(kù)結(jié)果,得到Q345的同性能鋼材在R=0.031 時(shí)對(duì)應(yīng)的ΔKth=168.9 MPa·mm0.5,在最不利荷載幅下,對(duì)不同α0裂紋的最初擴(kuò)展?fàn)顟B(tài)應(yīng)力強(qiáng)度因子幅值進(jìn)行計(jì)算,結(jié)果如圖19所示。從圖19可知:

    圖19 裂紋最初擴(kuò)展?fàn)顟B(tài)的應(yīng)力強(qiáng)度因子幅值Fig.19 SIF amplitude of crack at beginning of propagation

    1)以α0-p為基準(zhǔn),裂紋前緣ΔKI隨|α0|的增大而減小,ΔKII和ΔKIII隨|α0|的增大而先增大后減小,且始終小于ΔKth,可知U肋內(nèi)焊趾初始裂紋的擴(kuò)展由ΔKI>ΔKth引起,ΔKII和ΔKIII在裂紋最初擴(kuò)展時(shí)僅對(duì)擴(kuò)展角度產(chǎn)生影響;當(dāng)|α0|超過一定值時(shí),ΔKI,ΔKII和ΔKIII均小于ΔKth,此時(shí),初始裂紋不擴(kuò)展。

    2)ΔKeff與ΔKI的變化趨勢(shì)相似,說明ΔKeff可以合理反映裂紋擴(kuò)展過程中應(yīng)力強(qiáng)度因子的整體變化規(guī)律;但在|α0|較大的裂紋中(包括焊接裂紋或其他外因?qū)е碌牧鸭y),當(dāng)ΔKI,ΔKII和ΔKIII均小于ΔKth時(shí),仍可能有ΔKeff>ΔKth,故有效應(yīng)力強(qiáng)度因子幅值不能作為評(píng)估裂紋擴(kuò)展與否的依據(jù),在相關(guān)研究中應(yīng)分別通過3種應(yīng)力強(qiáng)度因子幅值與裂紋擴(kuò)展門檻值來判斷裂紋是否擴(kuò)展。

    4.4 初始裂紋擴(kuò)展的最大有效角度區(qū)間

    從前文分析可知,初始裂紋僅在某個(gè)有效角度區(qū)間內(nèi)擴(kuò)展。裂紋的幾何參數(shù)、裂紋處應(yīng)力狀態(tài)及鋼材類型均有可能對(duì)該區(qū)間產(chǎn)生影響。本文以形狀比為參數(shù)對(duì)有效角度區(qū)間進(jìn)行分析。

    4.4.1 形狀比參數(shù)分析

    取0.5 mm 作為工程中初始裂紋最大深度[27]以反映最不利情況。根據(jù)前文分析,運(yùn)用ΔKI和ΔKth的大小關(guān)系去判定雙面焊U 肋內(nèi)焊趾處初始裂紋的擴(kuò)展情況,以裂紋C1~C5 為計(jì)算對(duì)象,有效角度區(qū)與裂紋形狀比的關(guān)系分別如圖20和圖21所示。

    由圖20和圖21可得出:

    圖21 初始裂紋擴(kuò)展的有效角度區(qū)間與形狀比的關(guān)系Fig.21 Relationship between effective angle range of initial crack growth and shape ratio

    1)以|Δα0-e|=|α0-max-α0-min|定義U肋雙面焊內(nèi)焊趾初始裂紋擴(kuò)展的有效角度區(qū)間,α0-max和α0-min分別為區(qū)間的初始角度上下限。初始裂紋僅在|Δα0-e|內(nèi)具備擴(kuò)展能力,角度位于區(qū)間外的初始裂紋在結(jié)構(gòu)彈性階段內(nèi)不擴(kuò)展;該區(qū)間由裂紋前緣各點(diǎn)的ΔKI與母材ΔKth的大小關(guān)系決定;初始裂紋的α0-p不隨形狀比改變;在同一應(yīng)力比下,|Δα0-e|隨ω的增大而減小。在各形狀比下,U肋內(nèi)焊趾初始裂紋的|Δα0-e|計(jì)算結(jié)果如表3所示。

    表3 不同形狀比的初始裂紋擴(kuò)展有效角度區(qū)間Table 3 Effective angle range of initial crack growth with different shape ratios

    2)扁長(zhǎng)形裂紋(ω<0.5)在一定的初始角度范圍內(nèi)沿表面長(zhǎng)度方向不擴(kuò)展,僅在深度方向上發(fā)生局部擴(kuò)展,或前者的擴(kuò)展程度比后者的小,窄深形裂紋(ω>0.5)則與之相反。裂紋局部擴(kuò)展行為如圖22所示,該過程中裂紋形狀比發(fā)生變化,逐漸趨于整體擴(kuò)展。

    圖22 裂紋初始角度及形狀比對(duì)擴(kuò)展行為的影響Fig.22 Effect of initial angle and shape ratio of crack on propagation behavior

    4.4.2 有效角度區(qū)間的數(shù)學(xué)模型推導(dǎo)

    在推導(dǎo)|Δα0-e|~ω的數(shù)學(xué)模型前需對(duì)函數(shù)類型進(jìn)行確定,以保證模型的真實(shí)合理性。推導(dǎo)基于3個(gè)近似假設(shè)。

    1)選點(diǎn)假設(shè)。從圖20可知:①隨著|α0|的增大,初始裂紋表面點(diǎn)的應(yīng)力強(qiáng)度因子減小速度比其他點(diǎn)的大(ω=0.1 的裂紋除外,但此時(shí)表面點(diǎn)的應(yīng)力強(qiáng)度因子始終比其他點(diǎn)的?。?,該現(xiàn)象使得|Δα0-e|與表面點(diǎn)的應(yīng)力強(qiáng)度因子幅值無關(guān);②當(dāng)ω<0.5 時(shí),|Δα0-e|取決于裂紋最深點(diǎn)ΔKI與母材的ΔKth;當(dāng)ω>0.5 時(shí),|Δα0-e|取決于近表面區(qū)域(不含表面點(diǎn))的裂紋前緣最大ΔKI與母材的ΔKth;當(dāng)ω=0.5 時(shí),|Δα0-e|由以上兩者共同決定。因此,近似假設(shè)|Δα0-e|由最深點(diǎn)A和近表面點(diǎn)E(θ/π=0.1)的ΔKI與母材ΔKth的大小關(guān)系決定。

    圖20 初始裂紋擴(kuò)展的有效角度區(qū)間Fig.20 Effective angle range of initial crack growth

    2)線性假設(shè)。①不同初始角度裂紋的ΔKI-A和ΔKI-E分別在ω<0.5 和ω>0.5 時(shí)隨ω近似呈線性變化(如圖23所示),當(dāng)荷載工況確定時(shí),可假設(shè)其斜率kA和kE均為定值;②可近似假設(shè)ΔKI在接近ΔKth時(shí)隨α0呈線性變化(如圖24所示)。

    圖23 應(yīng)力強(qiáng)度因子幅值隨裂紋形狀比的變化(線性假設(shè)①)Fig.23 Change of SIF amplitude with crack shape ratio(linear assumption ①)

    圖24 應(yīng)力強(qiáng)度因子幅值隨裂紋初始角度的變化(線性假設(shè)②)Fig.24 Change of SIF amplitude with initial crack angle(linear assumption ②)

    3)對(duì)稱假設(shè):根據(jù)圖20,除表面點(diǎn)外,可近似假設(shè)裂紋前緣任一點(diǎn)處的ΔKI在任意ω時(shí)關(guān)于α0-p對(duì)稱分布,且α0-p不受ω影響。

    根據(jù)近似假設(shè)可得到如圖25所示的幾何模型,圖25中ωi(i=1,2,…)表示不同的形狀比,且ω1<ω2。

    根據(jù)假設(shè)1),可知|Δα0-e|~ω應(yīng)為分段函數(shù)形式。首先將ω=0.1 和ω=1 作為討論的形狀比上下限,應(yīng)以ω=0.5作為分段函數(shù)分界點(diǎn)。再根據(jù)假設(shè)2),結(jié)合圖25中的幾何關(guān)系可推出:

    圖25 初始裂紋擴(kuò)展的有效角度區(qū)間幾何推導(dǎo)模型Fig.25 Geometric derivation model of effective angle range of initial crack

    根據(jù)假設(shè)2),結(jié)合圖23,可得到

    將式(6)代入式(5),可得到

    根據(jù)假設(shè)2)及圖23可知kA和kE在荷載工況確定時(shí)為定值且小于0,故|Δα0-e|與ω以ω=0.5 為分界點(diǎn)呈分段線性負(fù)相關(guān),根據(jù)式(7)可得擬合公式的一般形式:

    式中:k1,b1,k2和b2為與初始裂紋ΔKI及材料門檻值ΔKth有關(guān)的常數(shù)。Q345 鋼材正交異性橋面板U肋內(nèi)側(cè)焊趾初始裂紋擴(kuò)展的最大有效角度區(qū)間擬合結(jié)果如圖26所示,在最不利應(yīng)力工況下,有k1=-66.286,b1=142.41,k2=-31.4,b2=123.48。

    圖26 初始裂紋擴(kuò)展的有效角度區(qū)間擬合結(jié)果Fig.26 Fitting result of initial crack growth effective angle range

    式(8)基于材料彈性階段內(nèi)的最大應(yīng)力幅及工程上一般認(rèn)為的最大初始裂紋深度得出,計(jì)算結(jié)果反映了該條件下|Δα0-e|的最大值,也可采用相同方法給出特定材料及應(yīng)力水平下的|Δα0-e|。該擬合計(jì)算式有助于研究者在裂紋擴(kuò)展計(jì)算中合理設(shè)置初始角度,全面考慮初始角度因素,有助于加深對(duì)疲勞裂紋擴(kuò)展特性的理解,也可為裂紋檢測(cè)周期的安排、維修加固對(duì)策的制定提供合理參考依據(jù)。

    5 結(jié)論

    1)U肋與頂板雙面焊構(gòu)造存在U肋內(nèi)焊趾開裂并沿板厚方向擴(kuò)展的疲勞破壞模式,最終裂紋的失穩(wěn)擴(kuò)展引發(fā)腹板截面的瞬斷導(dǎo)致連接失效;初始裂紋的角度具有不確定性。

    2)U 肋內(nèi)焊趾裂紋為I 型主導(dǎo)的I-II-III 型復(fù)合裂紋,裂紋的初始角度越大,II-III 型復(fù)合開裂特征越強(qiáng),KII和KIII分別影響裂紋在腹板厚度和結(jié)構(gòu)表面方向上的擴(kuò)展角度。

    3)U肋內(nèi)焊趾裂紋最深點(diǎn)的擴(kuò)展速率在擴(kuò)展過程中表現(xiàn)出“快速增長(zhǎng)+平穩(wěn)下降(緩沖)+失穩(wěn)”的三階段特性。在擴(kuò)展速率方面,裂紋初始角度僅對(duì)第一階段產(chǎn)生顯著影響;在疲勞壽命方面,裂紋初始角度主要影響結(jié)構(gòu)的疲勞壽命分布,初始角度越大,緩沖段在總疲勞壽命中的占比就越小,裂紋經(jīng)歷第一階段后會(huì)更快地達(dá)到失穩(wěn)擴(kuò)展,進(jìn)而引發(fā)截面瞬斷,這對(duì)于裂紋檢測(cè)周期的制定、維修加固安排等工作都更為不利。

    4)U 肋內(nèi)側(cè)焊趾初始裂紋是否擴(kuò)展由ΔKI與ΔKth的大小關(guān)系決定,有效應(yīng)力強(qiáng)度因子幅值可在一定程度上反映裂紋的擴(kuò)展規(guī)律,但不能作為裂紋是否擴(kuò)展的判據(jù)。

    5)U肋內(nèi)側(cè)焊趾初始裂紋僅可在某一有效角度區(qū)間內(nèi)擴(kuò)展。以初始裂紋深寬形狀比0.5 為分界,有效角度區(qū)間大小與形狀比近似呈分段線性負(fù)相關(guān),據(jù)此給出的擬合計(jì)算公式可為裂紋擴(kuò)展研究提供一定的基礎(chǔ)和依據(jù),也可用于判斷焊接裂紋或其他相似面狀缺陷的擴(kuò)展行為。本文在Q345鋼材前提下主要針對(duì)區(qū)間的最大值進(jìn)行了研究,不同鋼材種類及應(yīng)力水平下的具體公式有待于進(jìn)一步探索。

    猜你喜歡
    裂紋
    基于擴(kuò)展有限元的疲勞裂紋擴(kuò)展分析
    裂紋長(zhǎng)度對(duì)焊接接頭裂紋擴(kuò)展驅(qū)動(dòng)力的影響
    裂紋圓管彎曲承載能力研究
    裂紋敏感性鋼鑄坯表面質(zhì)量控制
    山東冶金(2019年6期)2020-01-06 07:45:58
    Epidermal growth factor receptor rs17337023 polymorphism in hypertensive gestational diabetic women: A pilot study
    42CrMo托輥裂紋的堆焊修復(fù)
    山東冶金(2019年3期)2019-07-10 00:54:06
    心生裂紋
    Overcoming scarring in the urethra:Challenges for tissue engineering
    微裂紋區(qū)對(duì)主裂紋擴(kuò)展的影響
    A7NO1鋁合金退火處理后焊接接頭疲勞裂紋擴(kuò)展特性
    焊接(2015年2期)2015-07-18 11:02:38
    亚洲va日本ⅴa欧美va伊人久久| 久久精品亚洲精品国产色婷小说| 桃色一区二区三区在线观看| 精品欧美国产一区二区三| 黄片大片在线免费观看| 黄色片一级片一级黄色片| 亚洲第一欧美日韩一区二区三区| 岛国在线观看网站| 在线观看一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 亚洲专区字幕在线| 国产一区在线观看成人免费| av视频在线观看入口| 桃色一区二区三区在线观看| 少妇熟女aⅴ在线视频| 午夜免费成人在线视频| 国产伦人伦偷精品视频| www国产在线视频色| 91在线观看av| 村上凉子中文字幕在线| 欧美日韩亚洲综合一区二区三区_| 免费在线观看黄色视频的| 黑人欧美特级aaaaaa片| 黄色 视频免费看| 亚洲久久久国产精品| 日韩大码丰满熟妇| 最近最新中文字幕大全电影3 | 丝袜美腿诱惑在线| 国产视频内射| 日本熟妇午夜| 久久亚洲精品不卡| 久久久久九九精品影院| 亚洲天堂国产精品一区在线| 国产亚洲精品久久久久久毛片| 久久精品aⅴ一区二区三区四区| 亚洲七黄色美女视频| 国产单亲对白刺激| 老司机在亚洲福利影院| 18禁国产床啪视频网站| 99在线视频只有这里精品首页| 黑人欧美特级aaaaaa片| 国产爱豆传媒在线观看 | 99久久久亚洲精品蜜臀av| 国产精品久久久久久精品电影 | 12—13女人毛片做爰片一| www.999成人在线观看| 在线观看免费日韩欧美大片| 高潮久久久久久久久久久不卡| 亚洲精品在线美女| 欧美激情极品国产一区二区三区| 很黄的视频免费| 午夜日韩欧美国产| 亚洲国产毛片av蜜桃av| 欧美性猛交╳xxx乱大交人| 久久中文字幕人妻熟女| 岛国视频午夜一区免费看| 亚洲精品一区av在线观看| 身体一侧抽搐| 精品久久久久久,| 在线免费观看的www视频| 国产精品av久久久久免费| 这个男人来自地球电影免费观看| 成年免费大片在线观看| 亚洲一区高清亚洲精品| www日本在线高清视频| 亚洲熟女毛片儿| 国产精品香港三级国产av潘金莲| 午夜福利免费观看在线| av超薄肉色丝袜交足视频| 久久久精品国产亚洲av高清涩受| 精品第一国产精品| 亚洲色图av天堂| 欧美av亚洲av综合av国产av| 久久久久久久午夜电影| 亚洲精华国产精华精| 欧美成人免费av一区二区三区| 在线播放国产精品三级| 欧美国产日韩亚洲一区| 国产片内射在线| 看免费av毛片| 麻豆一二三区av精品| 亚洲国产精品合色在线| 精品电影一区二区在线| 日韩 欧美 亚洲 中文字幕| 亚洲国产欧洲综合997久久, | 又黄又粗又硬又大视频| 国产av一区二区精品久久| 18禁裸乳无遮挡免费网站照片 | 一区福利在线观看| 午夜福利视频1000在线观看| 看片在线看免费视频| 丰满的人妻完整版| 欧美激情 高清一区二区三区| 男男h啪啪无遮挡| 岛国在线观看网站| 人人妻人人澡欧美一区二区| 亚洲男人天堂网一区| 国产成人精品久久二区二区91| 免费在线观看成人毛片| 精品国产乱子伦一区二区三区| 一二三四在线观看免费中文在| 波多野结衣高清无吗| 国产免费av片在线观看野外av| 丝袜美腿诱惑在线| 欧美日韩福利视频一区二区| 两个人看的免费小视频| 亚洲一区中文字幕在线| 此物有八面人人有两片| 欧美激情极品国产一区二区三区| 俺也久久电影网| 老司机午夜福利在线观看视频| 久久久久久久久免费视频了| 一本大道久久a久久精品| 一进一出抽搐gif免费好疼| 午夜福利一区二区在线看| 香蕉丝袜av| 女警被强在线播放| www.熟女人妻精品国产| 亚洲自偷自拍图片 自拍| 欧美一区二区精品小视频在线| 免费av毛片视频| 国产伦人伦偷精品视频| 中文字幕精品亚洲无线码一区 | netflix在线观看网站| 在线观看www视频免费| 国产又色又爽无遮挡免费看| 国产黄a三级三级三级人| 免费看十八禁软件| 久久久久久亚洲精品国产蜜桃av| 两个人看的免费小视频| 色老头精品视频在线观看| www.自偷自拍.com| 欧美黄色片欧美黄色片| 日日干狠狠操夜夜爽| 午夜精品久久久久久毛片777| 欧洲精品卡2卡3卡4卡5卡区| 在线国产一区二区在线| 成人特级黄色片久久久久久久| 久久国产乱子伦精品免费另类| 亚洲精品粉嫩美女一区| 在线观看午夜福利视频| 成年女人毛片免费观看观看9| 成人一区二区视频在线观看| 手机成人av网站| 久久久久亚洲av毛片大全| 国产精品久久久人人做人人爽| 欧美日韩亚洲综合一区二区三区_| 男人操女人黄网站| 欧美另类亚洲清纯唯美| 中文字幕高清在线视频| 欧美成人午夜精品| 免费看a级黄色片| 哪里可以看免费的av片| 在线观看午夜福利视频| 午夜福利成人在线免费观看| 窝窝影院91人妻| av超薄肉色丝袜交足视频| 亚洲第一欧美日韩一区二区三区| 精品久久久久久久人妻蜜臀av| 日韩精品青青久久久久久| 男女床上黄色一级片免费看| 大型av网站在线播放| 麻豆成人午夜福利视频| 国产高清videossex| 熟妇人妻久久中文字幕3abv| 亚洲av五月六月丁香网| 18禁观看日本| 日本黄色视频三级网站网址| 欧美日韩黄片免| 在线观看66精品国产| 久久人人精品亚洲av| 免费高清视频大片| 99国产精品一区二区蜜桃av| 免费在线观看成人毛片| 精品一区二区三区av网在线观看| 又黄又粗又硬又大视频| 每晚都被弄得嗷嗷叫到高潮| 男人舔女人下体高潮全视频| 国产av不卡久久| 一个人观看的视频www高清免费观看 | 男女之事视频高清在线观看| 老司机午夜十八禁免费视频| 亚洲人成网站高清观看| 亚洲va日本ⅴa欧美va伊人久久| 国产爱豆传媒在线观看 | 日韩欧美三级三区| 一夜夜www| 三级毛片av免费| 91大片在线观看| 欧美性长视频在线观看| 少妇粗大呻吟视频| 国产私拍福利视频在线观看| 免费女性裸体啪啪无遮挡网站| 免费在线观看成人毛片| 亚洲欧美日韩高清在线视频| 侵犯人妻中文字幕一二三四区| av中文乱码字幕在线| 亚洲成av片中文字幕在线观看| 色精品久久人妻99蜜桃| 校园春色视频在线观看| 男男h啪啪无遮挡| 亚洲一区二区三区色噜噜| 国内揄拍国产精品人妻在线 | 国产麻豆成人av免费视频| 亚洲人成网站高清观看| 欧美日韩精品网址| 国产野战对白在线观看| 久久久久久大精品| 亚洲色图 男人天堂 中文字幕| 亚洲精品国产一区二区精华液| 国产亚洲精品一区二区www| www.自偷自拍.com| 日韩三级视频一区二区三区| 男人的好看免费观看在线视频 | 精品久久久久久成人av| 色播在线永久视频| 日韩欧美 国产精品| 免费av毛片视频| 亚洲av熟女| 亚洲一区高清亚洲精品| 男男h啪啪无遮挡| 欧美又色又爽又黄视频| 黄片小视频在线播放| 亚洲黑人精品在线| 午夜福利成人在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| www日本黄色视频网| 亚洲av电影在线进入| 最近在线观看免费完整版| 黄色丝袜av网址大全| 一级a爱视频在线免费观看| 美女大奶头视频| 久久久国产精品麻豆| tocl精华| 久久久久久久久免费视频了| 一级黄色大片毛片| 巨乳人妻的诱惑在线观看| 久久人妻福利社区极品人妻图片| 非洲黑人性xxxx精品又粗又长| 高潮久久久久久久久久久不卡| 日本成人三级电影网站| 国产在线观看jvid| 久久中文字幕人妻熟女| 人人妻人人澡人人看| 国产精品野战在线观看| 在线观看www视频免费| 一进一出好大好爽视频| 欧美在线一区亚洲| 免费看日本二区| 麻豆一二三区av精品| 久久婷婷成人综合色麻豆| 亚洲成人国产一区在线观看| 久久久久九九精品影院| 正在播放国产对白刺激| 超碰成人久久| 神马国产精品三级电影在线观看 | 久久久国产精品麻豆| 国产熟女午夜一区二区三区| 黄片播放在线免费| 国产一级毛片七仙女欲春2 | 精品久久久久久久久久久久久 | 一区二区日韩欧美中文字幕| 一夜夜www| 久久精品国产99精品国产亚洲性色| 亚洲av中文字字幕乱码综合 | 给我免费播放毛片高清在线观看| 日本一本二区三区精品| 欧美日韩黄片免| 91在线观看av| 又黄又爽又免费观看的视频| 久久精品夜夜夜夜夜久久蜜豆 | 黄片小视频在线播放| 欧美成狂野欧美在线观看| 我的亚洲天堂| 老汉色∧v一级毛片| 19禁男女啪啪无遮挡网站| 麻豆成人午夜福利视频| 久久久久久久久免费视频了| 悠悠久久av| 欧美激情高清一区二区三区| 三级毛片av免费| 国产99久久九九免费精品| 村上凉子中文字幕在线| 免费在线观看完整版高清| 99riav亚洲国产免费| 人人澡人人妻人| 伊人久久大香线蕉亚洲五| 国产亚洲精品久久久久5区| 久久午夜综合久久蜜桃| 国产精品影院久久| 亚洲精华国产精华精| 国产激情欧美一区二区| 国产av在哪里看| 欧美在线一区亚洲| 黑人欧美特级aaaaaa片| 国内少妇人妻偷人精品xxx网站 | 国产1区2区3区精品| 日韩欧美 国产精品| 亚洲久久久国产精品| 男人操女人黄网站| 亚洲中文日韩欧美视频| 欧美又色又爽又黄视频| 欧美性猛交╳xxx乱大交人| 国产精品亚洲av一区麻豆| 国产精品国产高清国产av| 草草在线视频免费看| 亚洲五月天丁香| 久久伊人香网站| 亚洲三区欧美一区| 国产成年人精品一区二区| 久9热在线精品视频| 波多野结衣巨乳人妻| 人人妻人人澡欧美一区二区| 欧美激情高清一区二区三区| 国产v大片淫在线免费观看| 激情在线观看视频在线高清| 国产精品久久久久久精品电影 | 黑人巨大精品欧美一区二区mp4| 91麻豆精品激情在线观看国产| 两性夫妻黄色片| 成人欧美大片| 免费看a级黄色片| 99久久无色码亚洲精品果冻| 久久人妻av系列| 亚洲中文字幕日韩| 国产午夜精品久久久久久| 成年免费大片在线观看| 亚洲最大成人中文| 国产蜜桃级精品一区二区三区| 中亚洲国语对白在线视频| 午夜免费鲁丝| 人人妻,人人澡人人爽秒播| 欧美日韩亚洲综合一区二区三区_| 校园春色视频在线观看| 国产亚洲av高清不卡| 两个人免费观看高清视频| 日本精品一区二区三区蜜桃| 性色av乱码一区二区三区2| 欧美色视频一区免费| 精品欧美一区二区三区在线| 脱女人内裤的视频| 少妇熟女aⅴ在线视频| 国产精品综合久久久久久久免费| 久久久久久人人人人人| 日韩精品免费视频一区二区三区| 在线天堂中文资源库| 99久久综合精品五月天人人| 欧美精品啪啪一区二区三区| 88av欧美| 国产精品久久电影中文字幕| 亚洲成人免费电影在线观看| 成人国语在线视频| 国产精品一区二区免费欧美| 男人舔女人的私密视频| 国产精品日韩av在线免费观看| 最新在线观看一区二区三区| 白带黄色成豆腐渣| 亚洲久久久国产精品| 中文字幕精品亚洲无线码一区 | 久久久久亚洲av毛片大全| 一夜夜www| 精品第一国产精品| 日韩欧美国产一区二区入口| 国产成人精品久久二区二区91| 91大片在线观看| 好男人电影高清在线观看| 精品免费久久久久久久清纯| 99精品在免费线老司机午夜| 久久国产乱子伦精品免费另类| 久久久久九九精品影院| 99精品在免费线老司机午夜| 波多野结衣av一区二区av| 又大又爽又粗| 国产亚洲精品av在线| 男女下面进入的视频免费午夜 | 精品第一国产精品| 欧美又色又爽又黄视频| 99国产精品一区二区蜜桃av| 亚洲天堂国产精品一区在线| 此物有八面人人有两片| 亚洲av熟女| 欧美最黄视频在线播放免费| 美女免费视频网站| 十八禁人妻一区二区| 窝窝影院91人妻| 一a级毛片在线观看| 真人一进一出gif抽搐免费| 国产伦一二天堂av在线观看| 欧美午夜高清在线| 精品无人区乱码1区二区| 搡老熟女国产l中国老女人| 两个人免费观看高清视频| 久久久国产成人免费| 99国产精品一区二区三区| 美女高潮喷水抽搐中文字幕| 国产日本99.免费观看| 黄色成人免费大全| 国产精品自产拍在线观看55亚洲| 午夜福利成人在线免费观看| 欧美三级亚洲精品| 久久久久国产精品人妻aⅴ院| 亚洲国产看品久久| av福利片在线| 不卡av一区二区三区| 一进一出好大好爽视频| 亚洲免费av在线视频| 变态另类成人亚洲欧美熟女| 女生性感内裤真人,穿戴方法视频| 精品久久蜜臀av无| 变态另类成人亚洲欧美熟女| 国产精品自产拍在线观看55亚洲| 午夜成年电影在线免费观看| 97碰自拍视频| 最近最新免费中文字幕在线| 国产激情偷乱视频一区二区| 日韩成人在线观看一区二区三区| 亚洲五月天丁香| 久久精品国产综合久久久| 1024香蕉在线观看| 少妇的丰满在线观看| 国产午夜福利久久久久久| 狠狠狠狠99中文字幕| 免费在线观看日本一区| avwww免费| 国产成人av教育| 一级黄色大片毛片| 免费看美女性在线毛片视频| 国产在线精品亚洲第一网站| 成人三级黄色视频| 久9热在线精品视频| 老司机在亚洲福利影院| 免费高清在线观看日韩| 亚洲国产精品久久男人天堂| 男人舔奶头视频| 国产成人影院久久av| 一级毛片高清免费大全| 亚洲激情在线av| e午夜精品久久久久久久| 精品国产超薄肉色丝袜足j| а√天堂www在线а√下载| 亚洲国产精品sss在线观看| 亚洲久久久国产精品| 熟女少妇亚洲综合色aaa.| 久久久久久免费高清国产稀缺| 日本 欧美在线| 99久久99久久久精品蜜桃| 国产欧美日韩一区二区三| 国产成人欧美在线观看| 日韩免费av在线播放| 俺也久久电影网| 天天一区二区日本电影三级| 美女高潮喷水抽搐中文字幕| 亚洲avbb在线观看| 国产精品久久久久久亚洲av鲁大| 桃红色精品国产亚洲av| 在线十欧美十亚洲十日本专区| 女性生殖器流出的白浆| 久久国产精品人妻蜜桃| 精品国产一区二区三区四区第35| 真人一进一出gif抽搐免费| 国语自产精品视频在线第100页| 国产精品av久久久久免费| 高潮久久久久久久久久久不卡| 欧美黑人精品巨大| 亚洲免费av在线视频| 日韩欧美 国产精品| 色精品久久人妻99蜜桃| 又黄又粗又硬又大视频| 国产精品日韩av在线免费观看| 国产又黄又爽又无遮挡在线| 国产又爽黄色视频| 中文字幕精品亚洲无线码一区 | 日韩精品免费视频一区二区三区| 亚洲人成伊人成综合网2020| 久久婷婷成人综合色麻豆| 亚洲精品久久国产高清桃花| 久久精品亚洲精品国产色婷小说| 亚洲精品粉嫩美女一区| 久久午夜综合久久蜜桃| avwww免费| 少妇的丰满在线观看| 妹子高潮喷水视频| 一个人观看的视频www高清免费观看 | 国产亚洲精品综合一区在线观看 | 国产精品久久久久久精品电影 | 国产在线观看jvid| 热99re8久久精品国产| 日本黄色视频三级网站网址| 久久亚洲真实| 欧美日韩瑟瑟在线播放| 亚洲九九香蕉| 人人妻人人看人人澡| 亚洲精品粉嫩美女一区| 两个人看的免费小视频| 一区二区日韩欧美中文字幕| 亚洲免费av在线视频| 亚洲 国产 在线| 后天国语完整版免费观看| 侵犯人妻中文字幕一二三四区| 色播亚洲综合网| 超碰成人久久| 国产精品精品国产色婷婷| 国产1区2区3区精品| 国产精品九九99| 夜夜躁狠狠躁天天躁| 久久婷婷人人爽人人干人人爱| 午夜福利成人在线免费观看| 无遮挡黄片免费观看| 国产亚洲欧美98| av中文乱码字幕在线| 热99re8久久精品国产| 一进一出抽搐gif免费好疼| 99久久无色码亚洲精品果冻| 欧美性长视频在线观看| 香蕉丝袜av| 超碰成人久久| 51午夜福利影视在线观看| 搡老熟女国产l中国老女人| 久久久精品欧美日韩精品| 99精品在免费线老司机午夜| 变态另类丝袜制服| 怎么达到女性高潮| 国产91精品成人一区二区三区| 丰满人妻熟妇乱又伦精品不卡| av在线播放免费不卡| 午夜久久久在线观看| av在线播放免费不卡| 人妻久久中文字幕网| 美女午夜性视频免费| 亚洲第一欧美日韩一区二区三区| 欧美av亚洲av综合av国产av| 国产单亲对白刺激| 搡老岳熟女国产| 成熟少妇高潮喷水视频| 啦啦啦免费观看视频1| 精品午夜福利视频在线观看一区| 午夜福利18| 狂野欧美激情性xxxx| 亚洲国产毛片av蜜桃av| 免费看日本二区| 成人三级黄色视频| 欧美一级毛片孕妇| 巨乳人妻的诱惑在线观看| 老司机午夜福利在线观看视频| 亚洲国产毛片av蜜桃av| 欧美一级毛片孕妇| 午夜a级毛片| 成人18禁在线播放| 亚洲激情在线av| 99国产精品99久久久久| 日韩三级视频一区二区三区| 高潮久久久久久久久久久不卡| 亚洲成av片中文字幕在线观看| 高潮久久久久久久久久久不卡| 亚洲精品国产精品久久久不卡| 国产1区2区3区精品| 久久久久久国产a免费观看| 久久国产亚洲av麻豆专区| 女性生殖器流出的白浆| 日日夜夜操网爽| 90打野战视频偷拍视频| 91在线观看av| 国产伦人伦偷精品视频| 婷婷亚洲欧美| 色精品久久人妻99蜜桃| 午夜两性在线视频| 午夜免费激情av| 精品国产超薄肉色丝袜足j| 欧美成狂野欧美在线观看| 99久久久亚洲精品蜜臀av| 日本一区二区免费在线视频| 国产精品影院久久| 国产区一区二久久| 亚洲第一av免费看| 99热这里只有精品一区 | 无人区码免费观看不卡| av天堂在线播放| 91麻豆av在线| 国产主播在线观看一区二区| 天堂影院成人在线观看| 天堂动漫精品| 看黄色毛片网站| 免费观看人在逋| 国产精品国产高清国产av| 一夜夜www| 亚洲av五月六月丁香网| 视频在线观看一区二区三区| 国产午夜福利久久久久久| 国产精品美女特级片免费视频播放器 | 精品国产乱码久久久久久男人| 亚洲最大成人中文| 夜夜爽天天搞| 国产亚洲精品久久久久久毛片| 国产av又大| 亚洲精品在线观看二区| 日本撒尿小便嘘嘘汇集6| 99久久久亚洲精品蜜臀av| 母亲3免费完整高清在线观看| 国产亚洲精品综合一区在线观看 | 免费在线观看影片大全网站| 视频在线观看一区二区三区| 欧美绝顶高潮抽搐喷水| 国产精品日韩av在线免费观看| 久久人妻av系列| 亚洲片人在线观看| 可以免费在线观看a视频的电影网站| 日韩欧美 国产精品| 久久久久国内视频| 成人永久免费在线观看视频| 一二三四社区在线视频社区8| www.熟女人妻精品国产| 亚洲国产毛片av蜜桃av| x7x7x7水蜜桃| 精品电影一区二区在线| 在线观看免费日韩欧美大片|