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

    鉆柱表面裂紋在拉扭作用下的彈塑性斷裂性能*

    2014-06-07 05:59:48趙廣慧汪浩瀚
    關(guān)鍵詞:鉆柱彈塑性井口

    趙廣慧,汪浩瀚,石 健,趙 莉

    西南石油大學(xué)機(jī)電工程學(xué)院,四川 成都 610500

    鉆柱表面裂紋在拉扭作用下的彈塑性斷裂性能*

    趙廣慧,汪浩瀚,石 健,趙 莉

    西南石油大學(xué)機(jī)電工程學(xué)院,四川 成都 610500

    在鉆井作業(yè)過程中,井口附近鉆桿受到的拉力和扭矩往往達(dá)到整個(gè)鉆柱的最大值,尤其是在上提鉆柱解卡、震擊解卡以及蹩鉆、卡鉆情況發(fā)生時(shí),井口附近鉆桿受力更為惡劣。針對井口附近鉆桿的斷裂問題,以某超深直井的實(shí)際工況載荷為例,將鉆桿外表面的缺陷簡化為橫向半橢圓裂紋,以30CrMo作為鉆桿材料并通過實(shí)驗(yàn)擬合得到其Ramberg–Osgood(R–O)本構(gòu)模型,研究了拉力和扭矩組合作用下裂紋的彈塑性斷裂特性。得到鉆進(jìn)工況下鉆柱外表面裂紋擴(kuò)展至不同深度的臨界形狀比以及塑性成分在J積分中所占的比例;定量分析了拉力波動和扭矩波動對井口附近鉆桿彈塑性斷裂性能的影響,發(fā)現(xiàn)扭矩波動對裂紋擴(kuò)展的貢獻(xiàn)遠(yuǎn)大于拉力波動,避免鉆柱受到扭轉(zhuǎn)沖擊對延長鉆柱壽命有重要意義。

    鉆柱;半橢圓表面裂紋;J積分;彈塑性斷裂;有限元

    趙廣慧,汪浩瀚,石 健,等.鉆柱表面裂紋在拉扭作用下的彈塑性斷裂性能[J].西南石油大學(xué)學(xué)報(bào):自然科學(xué)版,2014,36(3):157–164.

    Zhao Guanghui,Wang Haohan,Shi Jian,et al.Elastic-plastic Fracture Mechanics Analyses of Surface Cracks in Drill String Subjected to Combined Loading[J].Journal of Southwest Petroleum University:Science&Technology Edition,2014,36(3):157–164.

    引言

    隨著鉆井技術(shù)向深井、超深井以及深海油氣田開采的發(fā)展,鉆柱的可靠性也受到了越來越嚴(yán)峻的挑戰(zhàn)。在鉆進(jìn)和起下鉆過程中,井口處鉆柱的受力很復(fù)雜。起下鉆時(shí)井口處鉆柱受到最大拉力,如果猛提、猛剎,會使井口處鉆柱受到的軸向拉力大大增加;鉆進(jìn)時(shí),井口處鉆柱所受的拉力和扭矩都很大[1-2],成為鉆柱發(fā)生斷裂事故的危險(xiǎn)位置。

    對于油氣田鉆井作業(yè)中的鉆具斷裂問題,從現(xiàn)場大多數(shù)鉆具的斷口分析可以判斷為裂紋的動態(tài)擴(kuò)展所致[3]。文獻(xiàn)[4-5]的研究表明,具有任何形狀的初始表面缺陷在經(jīng)過幾個(gè)疲勞裂紋擴(kuò)展的循環(huán)后,都會發(fā)展成半橢圓裂紋的形狀,因此將鉆柱外表面的橫向裂紋取為半橢圓形能較好地反映實(shí)際問題。由于裂紋形狀的復(fù)雜性,對于實(shí)心或空心圓柱表面裂紋的研究主要是通過數(shù)值方法進(jìn)行研究,且多為線彈性計(jì)算[6-14]。其中,文獻(xiàn)[8]檢驗(yàn)了二維和三維裂紋模型的位移外推法和J積分方法,對于是否在裂尖附近采用奇異元以及是否采用正交網(wǎng)格得到的J積分結(jié)果進(jìn)行了對比研究,說明J積分不需要在裂尖周圍采用奇異網(wǎng)格,也不要求劃分正交網(wǎng)格,只需要保證回路積分的守恒性,很適合于研究具有復(fù)雜結(jié)構(gòu)的裂紋體。關(guān)于圓柱體表面半橢圓裂紋的彈塑性斷裂性能研究的文獻(xiàn)很少,且多為單向荷載作用的情況。文獻(xiàn)[15]利用J積分計(jì)算了彎矩作用下圓管的彈塑性斷裂性能,通過和參考應(yīng)力法的結(jié)果進(jìn)行對照,說明了參考應(yīng)力法的應(yīng)用局限。文獻(xiàn)[16]計(jì)算了拉力作用下圓柱的低周疲勞特性,表明塑性情況下,J積分強(qiáng)烈依賴于圓柱和裂紋的幾何參數(shù)以及材料Ramberg–Osgood(R–O)本構(gòu)模型中的硬化指數(shù)。

    鉆桿在正常情況下,都是在彈性范圍內(nèi)工作。但是當(dāng)裂紋出現(xiàn)時(shí),特別是相對裂紋深度較大、剩余連接韌帶較薄的情況下,或者是非常惡劣的受力條件下,裂紋體的非彈性變形對斷裂性能的影響就不再是小量。本文采用J積分?jǐn)?shù)值研究鉆柱外表面橫向半橢圓裂紋在拉力和扭矩組合作用下的彈塑性斷裂特性,研究鉆井工況載荷條件下裂紋擴(kuò)展至不同深度時(shí)裂紋面的形狀,評價(jià)拉力波動和扭矩波動對井口附近鉆柱斷裂性能的影響。

    1 力學(xué)模型

    井口附近的鉆桿主要受到軸向拉力T和扭矩M作用,其裂紋面與鉆柱軸線相垂直的情況最為危險(xiǎn)[17]。以某超深直井用5 in.(1 in.=25.4 mm)鉆桿為例,研究其本體外表面橫向半橢圓裂紋在鉆井載荷作用下的斷裂特性。鉆柱的外徑D=127.0 mm,內(nèi)徑d=108.6 mm,壁厚t=9.19 mm。鉆桿及其裂紋的受力和尺寸如圖1所示,其中裂紋的幾何形狀可由無因次的形狀比a/c和相對裂紋深度a/t表示(a—半橢圓裂紋的深度,m;c—半橢圓裂紋的半寬度,m)。

    圖1 拉力和扭矩組合作用下,鉆柱及其含裂紋橫截面示意圖Fig.1 Sketch of forced drill string and arbitrary crack shape

    以30CrMo作為鉆柱材料,通過拉伸試驗(yàn)測出其名義應(yīng)力–應(yīng)變曲線,再轉(zhuǎn)換成真應(yīng)力–真應(yīng)變曲線[18],以R–O本構(gòu)模型對實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合

    式中:

    ε—應(yīng)變,無因次;

    εs—屈服應(yīng)變,無因次;

    σ—應(yīng)力,MPa;

    σs—屈服應(yīng)力,MPa;

    α,n—材料本構(gòu)常數(shù),無因次。

    30CrMo的材料常數(shù):

    彈性模量E=225.3 GPa;

    泊松比μ=0.28;

    屈服應(yīng)力σs=700.1 MPa;

    屈服應(yīng)變εs=0.003 1。

    材料本構(gòu)常數(shù):

    α=3.024;

    n=12.89。

    材料本構(gòu)曲線如圖2所示。

    圖2 30CrMo材料拉伸實(shí)驗(yàn)得到的真應(yīng)力–真應(yīng)變及其擬合的R–O本構(gòu)曲線Fig.2 True stress-strain relationship and fitted R–O constitutional model of 30CrMo

    2 有限元模型及其驗(yàn)證

    由于問題非對稱,數(shù)值計(jì)算中需要模擬整個(gè)圓柱,將裂紋置于鉆柱模型的中部,其網(wǎng)格劃分如圖3所示。

    圖3 有限元網(wǎng)格模型與裂紋面截面的局部網(wǎng)格劃分Fig.3 The finite element model and enlarged drawing of crack section

    計(jì)算模型的長度取為2×150 mm[5],下端面約束其全部自由度;上端面自由,并施加拉力和扭矩。利用ABAQUS對裂紋體進(jìn)行彈塑性分析,在裂紋前緣線的垂直面內(nèi)取6個(gè)積分回路,以中間的4個(gè)回路積分的平均值作為該裂紋點(diǎn)的J積分[19],描述裂紋體拉伸和扭轉(zhuǎn)組合載荷作用下的彈塑性斷裂性能。

    由于缺少可借鑒的中空圓柱外表面橫向半橢圓裂紋的彈塑性J積分結(jié)果,本文關(guān)于J積分的數(shù)值計(jì)算有效性是通過與前人得到的彈性計(jì)算結(jié)果進(jìn)行對比驗(yàn)證的,同時(shí)驗(yàn)證回路積分的路徑無關(guān)性。以文獻(xiàn)[20]中提到的Newman J C和Wang Q Z關(guān)于該裂紋體分別在拉力和扭矩單獨(dú)作用下裂紋最深點(diǎn)A的K因子結(jié)果為參照。拉力單獨(dú)作用下,A點(diǎn)的J積分按照式(2)轉(zhuǎn)換成K因子

    式中:

    K—應(yīng)力強(qiáng)度因子,N/m3/2;

    E—楊氏彈性模量,Pa;

    J—J積分,kN/m;

    μ—泊松比,無因次。

    進(jìn)一步利用式(3)進(jìn)行無因次化

    式中:

    F—無因次K因子;

    σ—裂紋處的名義拉應(yīng)力,MPa。

    對于裂紋最深點(diǎn)A,拉力單獨(dú)作用下的有限元計(jì)算值與Newman J C所得結(jié)果之間的最大相對誤差僅為0.70%,如圖4a所示;扭矩單獨(dú)作用下的計(jì)算值與Wang Q Z所得結(jié)果的最大相對誤差為2.60%,如圖4b所示;在所進(jìn)行的8組計(jì)算中,沿所取的4個(gè)積分回路得到的J積分的最大相對誤差只有0.12%,數(shù)值結(jié)果很好地保證了J積分的守恒性,由此說明,本文數(shù)值計(jì)算模型的可靠性,并將其直接推廣應(yīng)用到該裂紋體的彈塑性計(jì)算。

    圖4 有限元結(jié)果驗(yàn)證Fig.4 Finite element model validation

    3 計(jì)算結(jié)果

    對于某超深井的井口附近鉆桿,考慮拉力T為1 050 kN、扭矩M為30 kN·m同時(shí)作用的鉆井作業(yè)載荷,鉆桿外表面橫向半橢圓裂紋的相對深度a/t在0.1~0.9變化,裂紋形狀比a/c在0.5~0.9變化。

    基于方程(1)擬合得到的彈塑性本構(gòu)關(guān)系,對含裂紋的鉆桿進(jìn)行彈塑性斷裂力學(xué)分析。鑒于無因次分析的優(yōu)點(diǎn),將J積分按式(4)進(jìn)行無因次化處理

    式中:J0—無因次J積分。

    3.1 裂紋的形狀比

    在交變應(yīng)力作用下,鉆柱外表面的初始缺陷通常會擴(kuò)展成為半橢圓裂紋。根據(jù)Paries公式,在加載方式、環(huán)境等外部條件相同的情況下,裂紋的擴(kuò)展速率只與應(yīng)力強(qiáng)度因子的變化幅值有關(guān),在彈塑性情況下,則由J積分的變化幅值表示[16]。半橢圓裂紋前緣各點(diǎn)的擴(kuò)展速率隨相對位置θ、相對裂紋深度a/t和裂紋形狀比a/c的不同而變化。對于同一相對裂紋深度的不同裂紋,存在一個(gè)臨界形狀比(a/c)cr,此時(shí)裂紋前緣最深點(diǎn)A和邊界相鄰點(diǎn)B處的J積分相等,該臨界形狀比確定的半橢圓可近似作為裂紋擴(kuò)展至此深度時(shí)的裂紋面形狀,不同裂紋相對深度條件下臨界形狀比的確定對于裂紋疲勞擴(kuò)展的研究有重要的意義。

    對于鉆井載荷(T=1 050 kN,M=30 kN·m)作用下的井口附近鉆桿,計(jì)算了相對裂紋深度a/t分別為0.2,0.4,0.6及0.8時(shí),裂紋前緣最深點(diǎn)A和邊界相鄰點(diǎn)B的無因次J積分隨形狀比a/c的變化(如圖5所示,圖中,實(shí)線對應(yīng)裂紋最深點(diǎn)A,虛線對應(yīng)邊界相鄰點(diǎn)B),兩條線的交點(diǎn)對應(yīng)的a/c即為臨界形狀比(如表1所示)??梢钥吹剑S著裂紋相對深度的增加,臨界形狀比是逐漸減小的,裂紋前緣會越來越平坦。

    圖5 無因次的J0隨裂紋形狀比的變化Fig.5 The variation of J0along with crack aspect ratios

    表1 不同裂紋相對深度對應(yīng)的臨界形狀比Tab.1 The transition aspect ratio varying with the relative crack depth

    3.2 裂紋的彈塑性斷裂性能

    鉆桿在正常工作狀態(tài)下是處于彈性范圍內(nèi)的,但是當(dāng)初始裂紋發(fā)生較大范圍的擴(kuò)展時(shí),材料的彈塑性影響可能就不再是小量了。假設(shè)裂紋在擴(kuò)展過程中,其裂紋尖端的曲線近似為具有臨界形狀比的半橢圓,研究拉力和扭矩組合作用下,鉆柱表面裂紋的彈塑性斷裂性能隨裂紋深度的變化,評價(jià)塑性變形對J積分的影響。

    圖6給出裂紋相對深度a/t分別為0.2,0.5,0.7及0.9等4種情況下,裂紋尖端J積分在裂紋前緣分布的線彈性計(jì)算結(jié)果(圖中虛線)和彈塑性計(jì)算結(jié)果(圖中實(shí)線);表2給出了不同相對裂紋深度條件下,裂紋前緣最深點(diǎn)A的塑性J積分JP在整個(gè)J積分中所占的比例。

    圖6和表2中,除了a/t=0.9的結(jié)果對應(yīng)于a/c=0.68以外,其他相對裂紋深度條件下的結(jié)果均對應(yīng)于各自的臨界形狀比(見表1)??梢钥吹剑S著裂紋的擴(kuò)展,裂紋前緣的J積分是單調(diào)增大的,且裂紋越深,其增大的速度越快;在裂紋相對深度較小時(shí)(比如a/t=0.2),裂尖附近的塑性區(qū)很小,可以近似為線彈性斷裂;隨著裂紋深度的增大,裂紋前緣附近的塑性區(qū)越來越大,其J積分沿裂紋前緣的分布曲線與線彈性計(jì)算結(jié)果的差異也越來越大,當(dāng)a/t>0.7以后塑性J積分JP在整個(gè)J積分中所占的比例迅速增大,塑性變形對裂紋斷裂性能的影響不再是小量;在裂紋相對深度達(dá)到0.9時(shí),由于裂紋最深點(diǎn)附近的韌帶很薄,發(fā)生了明顯的塑性變形,使得J積分在裂紋最深點(diǎn)附近有顯著的增大,JP/J達(dá)到了23.9%。

    圖6 T=1 050 kN,M=30 kN·m條件下,無因次J積分沿裂紋前沿的分布Fig.6 The variations of J0along the crack front under T=1 050 kN,M=30 kN·m,in which real lines are elastic-plastic results and dashed lines are elastic results

    表2 T=1 050 kN,M=30 kN·m條件下,JP/J隨裂紋相對深度a/t變化Tab.2 The variations of JP/J along with a/t under T=1 050 kN,M=30 kN·m

    3.3 載荷對斷裂性能的影響

    在需要進(jìn)行上提鉆桿解卡或使用震擊器解卡操作中,上提力或震擊力可能達(dá)到整個(gè)鉆柱重量的數(shù)倍或接近最大鉤載,而在鉆頭發(fā)生蹩鉆或卡鉆時(shí),鉆桿會承受瞬時(shí)超大扭矩的作用。對于井口附近的鉆桿,以相對裂紋深度a/t=0.2、形狀比a/c=0.79為例,研究拉力或扭矩變化對其斷裂性能的影響,評價(jià)這兩種主要載荷的波動對鉆桿安全的損害。

    為便于定量研究J積分關(guān)于拉力或扭矩的變化率,將J積分表示為

    式中:

    τ—剪應(yīng)力,Pa;

    J01—對應(yīng)于正應(yīng)力σ的無因次J積分;

    J02—對應(yīng)于剪應(yīng)力τ的無因次J積分。

    用軸向拉力T和扭矩M可以將式(5)表示為

    式中:

    S—鉆柱的橫截面積,m2;

    T—軸向拉力,kN;

    IP—極慣性矩,m4;

    D—鉆柱外徑,m;

    M—扭矩,kN·m。

    這樣在得到無因次的J01和J02沿裂紋前緣的分布曲線后(J01和J02可以通過數(shù)值計(jì)算兩組拉扭組合作用下的J積分后由式(6)換算得到),就可以方便地得到該裂紋在任意拉力和扭矩作用下裂紋前緣的J積分,同時(shí)由于式(6)將J積分表示成了關(guān)于T和M的顯式形式,可以直接得到J積分關(guān)于T和M的變化率。

    式(5)和式(6)的實(shí)質(zhì)是將拉扭組合作用下的J積分等價(jià)為拉力和扭矩分別單獨(dú)作用下的J積分之和,即利用了疊加原理。對于拉扭組合作用下含半橢圓表面裂紋的圓柱體,文獻(xiàn)[11]在彈性范圍內(nèi)研究發(fā)現(xiàn),當(dāng)扭矩引起的剪應(yīng)力小于拉力引起的正應(yīng)力時(shí),扭矩導(dǎo)致的裂紋面翹曲對疊加原理的影響非常微小,此時(shí)組合載荷作用下的K因子與各種載荷單獨(dú)作用下的K因子之和近似相等,疊加原理成立。對于本文所研究的井口附近鉆桿,考慮其拉力圍繞T=1 050 kN上下波動30%,扭矩圍繞M=30 kN·m上下波動30%,其最大扭矩對鉆桿產(chǎn)生的剪應(yīng)力仍小于最小拉力產(chǎn)生的正應(yīng)力,在彈性范圍內(nèi)符合疊加原理適用的條件。

    為驗(yàn)證疊加原理在井口附近鉆桿彈塑性斷裂研究中的適用性,任取兩組載荷:(1)T=1 050 kN和M=20 kN·m;(2)T=1 050 kN和M=40 kN·m。通過有限元計(jì)算得到其各自的J積分后,由式(6)反算得到J01與J02沿裂紋前緣的分布如圖7所示。

    圖7 無因次的J01和J02沿裂紋前緣的分布Fig.7 The variations of J01and J02along the crack front

    首先固定扭矩 M=30 kN·m,拉力 T圍繞1 050 kN上下波動30%,分別通過疊加原理式(6)和有限元計(jì)算得到的A點(diǎn)J積分隨拉力的變化如圖8所示,兩者最大相對誤差為9.1%。然后固定拉力T=1 050 kN,扭矩M圍繞30 kN·m上下波動30%,通過疊加原理和有限元計(jì)算得到A點(diǎn)J積分隨扭矩的變化如圖9所示,兩者最大相對誤差僅為2.1%。由此說明對于拉扭組合作用下的井口附近鉆桿,疊加原理也是近似成立的。

    圖8 M=30 kN·m時(shí)A點(diǎn)的J積分隨拉力T的變化Fig.8 J-integral of point A varying along tension force T under M=30 kN·m

    圖9 T=1 050 kN時(shí)A點(diǎn)的J積分隨扭矩M的變化Fig.9 J-integral of point A varying along torque M under T=1 050 kN

    J積分對拉力變化和對扭矩變化的敏感度可以通過對式(6)求偏導(dǎo)得到

    在本文所研究的拉力 T=1 050 kN、扭矩M=30 kN·m的條件下,裂紋最深點(diǎn)A的J積分隨拉力的變化率為 2.14×10?3m?1,隨扭矩的變化率為2.83×10?2m?2。保持扭矩 M=30 kN·m不變,當(dāng)軸向拉力圍繞 1 050 kN上下波動 30%時(shí),在(1.50~2.78)×10?3m?1變化;保持拉力T=1 050 kN不變,扭矩圍繞30 kN·m上下波動30%時(shí),在(1.98~3.69)×10?2m?2變化??紤]到井口附近鉆桿扭矩波動的幅值遠(yuǎn)大于拉力波動的幅值[1],扭矩波動對井口附近鉆桿斷裂的威脅遠(yuǎn)大于拉力波動。

    4 結(jié) 論

    將鉆柱外表面的缺陷簡化為橫向半橢圓裂紋,以30CrMo作為鉆柱材料并通過實(shí)驗(yàn)擬合得到R–O本構(gòu)模型,以某超深直井井口附近鉆桿的實(shí)際工況載荷為例,數(shù)值研究了井口附近鉆桿的彈塑性斷裂性能,得到以下結(jié)論:

    (1)半橢圓表面裂紋的臨界形狀比隨裂紋相對深度的增大而減小,當(dāng)裂紋相對深度由0.1增大到0.8,其臨界形狀比由0.80降低到0.68,即隨著裂紋擴(kuò)展,其裂紋前緣越來越平坦。

    (2)隨著裂紋的擴(kuò)展,裂紋前緣附近的塑性區(qū)越來越大,當(dāng)相對裂紋深度a/t>0.7以后,J積分中的塑性成分JP在整個(gè)J積分中所占的比例迅速增大,采用計(jì)入鉆柱彈塑性斷裂特性的J積分才能更好地描述鉆柱疲勞裂紋擴(kuò)展的全壽命。

    (3)對于井口附近鉆桿所受到的拉扭組合載荷,驗(yàn)證了疊加原理對其外表面半橢圓裂紋斷裂性能的適用性,通過研究J積分分別相對于拉力和扭矩的變化率發(fā)現(xiàn),扭矩波動對J積分的影響遠(yuǎn)大于拉力波動對J積分的影響,鉆柱表面裂紋在扭矩波動下更容易擴(kuò)展,對鉆柱的危害也更大。預(yù)防鉆井過程中的蹩鉆和卡鉆現(xiàn)象,避免鉆柱受到扭轉(zhuǎn)沖擊對延長鉆柱壽命有重要意義。

    [1]林鐵軍,練章華,張俊良,等.氣體鉆與泥漿鉆全井段鉆柱動力學(xué)對比研究[J].西南石油大學(xué)學(xué)報(bào):自然科學(xué)版,2011,33(1):139–143. Lin Tiejun,Lian Zhanghua,Zhang Junliang,et al.Comparative study on dynamics behaviors of the full drillstring between gas drilling and mud drilling[J].Journal of Southwest Petroleum University:Science&Technology Edition,2011,33(1):139–143.

    [2]林鐵軍,練章華,曾曉健,等.應(yīng)用XFEM模擬研究鉆桿裂紋擴(kuò)展過程[J].重慶大學(xué)學(xué)報(bào),2010,33(7):123–128. Lin Tiejun,Lian Zhanghua,Zeng Xiaojian,et al.Simulation on crack growth of drill pipe with XFEM[J].Journal of Chongqing University,2010,33(7):123–128.

    [3]董洪生.川東北地區(qū)鉆具故障現(xiàn)狀分析[R].上海:川東北工區(qū)鉆具應(yīng)用研討會,2007.

    [4]Lin X B,Smith R A.Shape growth simulation of surface cracks in tension fatigued round bars[J].International Journal of Fatigue,1997,19(6):461–469.

    [5]Carpinteri A.Shape change of surface cracks in round bars under cyclic axial loading[J].International Journal of Fatigue,1993,15(1):21–26.

    [6]Manuel da Fonte,Manuel de Freitas.Stress intensity factors for semi-elliptical surface cracks in round bars under bending and torsion[J].International Journal of Fatigue,1999,21:457–463.

    [7]Toribio J,′Alvarez N,Gonz′alez B,et al.A critical review of stress intensity factor solutions for surface cracks in round bars subjected to tension loading[J].Engineering Failure Analysis,2009,16:794–809.

    [8]Courtin S,Gardin C,Bezine G,et al.Advantages of the J-integral approach for calculating stress intensity factors when using the commercial finite element software ABAQUS[J].Engineering Fracture Mechanics,2005,72:2174–2185.

    [9]Toribio J,Matos J C,Gonzalez B,et al.Numerical modelling of crack shape evolution for surface flaws in round bars under tensile loading[J].Engineering Failure Analysis,2009,16:618–630.

    [10]Carpinteri A,Brighenti R.Part-through cracks in round bars under cyclic combined axial and bending loading[J]. International Journal of Fatigue,1996,18(1):33–39.

    [11]Ismail A E,Ariffin A K,Abdullah S,et al.Stress intensity factors for surface cracks in round bar under single and combined loadings[J].Meccanica,2012,47:1141–1156.

    [12]Shahani A R,Habibi S E.Stress intensity factors in a hollow cylinder containing a circumferential semi-elliptical crack subjected to combined loading[J].International Journal of Fatigue,2007,29:128–140.

    [13]Kou K P,Burdekin F M.Stress intensity factors for a wide range of long-deep semi-elliptical surface cracks,partly through-wall cracks and fully through-wall cracks in tubular members[J].Engineering Fracture Mechanics,2006,73:1693–1710.

    [14]Predan J,Mocilnik V,Gubeljak N.Stress intensity factors for circumferential semi-elliptical surface cracks in a hollow cylinder subjected to pure torsion[J].Engineering Fracture Mechanics,2013,105:152–168.

    [15]Mechab B,Serier B,Bouiadjra B B,et al.Linear and nonlinear analyses for semi-elliptical surface cracks in pipes under bending[J].International Journal of Pressure Vessels and Piping,2011,88:57–63.

    [16]Findley K O,Koh S W,Saxena A.J-integral expressions for semi-elliptical cracks in round bars[J].International Journal of Fatigue,2007,29:822–828.

    [17]張俊清.高速列車空心車軸表面裂紋應(yīng)力強(qiáng)度因子研究[D].北京:北京交通大學(xué),2011.

    [18]鄒吉權(quán),荊洪陽,霍立興,等.31Si2MnCrMoVE鋼薄板試樣表面裂紋斷裂韌度測試[J].機(jī)械工程學(xué)報(bào),2007,43(8):212–217. Zou Jiquan,Jing Hongyang,Huo Lixing,et al.Surfacecrack fracture toughness test of 31Si2MnCrMoVE steelon thin plate[J].Chinese Journal of Mechanical Engineering,2007,43(8):212–217.

    [19]Kim Y J,Kim J S,Park Y J,et al.Elastic-plastic fracture mechanics method for finite internal axial surface cracks in cylinders[J].Engineering Fracture Mechanics,2004,71:925–944.

    [20]中國航空研究院.應(yīng)力強(qiáng)度因子手冊[M].增訂版.北京:科學(xué)出版社,1993.

    編輯:張?jiān)圃?/p>

    編輯部網(wǎng)址:http://zk.swpuxb.com

    Elastic-plastic Fracture Mechanics Analyses of Surface Cracks in Drill String Subjected to Combined Loading

    Zhao Guanghui,Wang Haohan,Shi Jian,Zhao Li
    School of Mechanical Engineering,Southwest Petroleum University,Chengdu,Sichuan 610500,China

    The drill pipe near the surface stands the largest tension and torsion in the whole drillstring especially during drill pipe sticking or bit bouncing.As an example,a 5 inch drill pipe made of 30CrMo that was near the well head of an ultradeep straight well was analyzed here.Simplifying the crack of the drill pipe as circumferential semi-elliptical surface crack and using the fitted Ramberg-Osgood equation as constitutional relation,we simulated elastic-plastic fracture feature of the drill string under combined loading of axial force and torsion corresponding to its drilling behavior.The transition aspect ratios are determined for the given combined loading and the crack shape could be determined approximately corresponding to different crack depths.The contribution of the plastic J-integral is estimated and the effect of material nonlinearity on the crack propagation is described.Variations of J-integral as a function of axial force and torsional moment respectively are compared,and the dangers from tension and torque fluctuation are evaluated quantitatively.Compared with tension,fluctuation in torque contributes much more to crack propagation.It is an important significance for increasing drill string′s service life to avoid the drill string from severe torsion impact.

    drill string;semi-elliptical surface crack;J-integral;elastic-plastic fracture;finite element

    http://www.cnki.net/kcms/doi/10.11885/j.issn.1674-5086.2013.12.09.01.html

    趙廣慧,1971年生,女,漢族,河北石家莊人,教授,博士(后),主要從事與石油天然氣裝備相關(guān)的非線性動力學(xué)、斷裂力學(xué)和流體力學(xué)方面的研究。E-mail:wy_zgh@126.com

    汪浩瀚,1989年生,男,漢族,四川阿壩人,碩士研究生,主要從事桿管柱斷裂力學(xué)方面的研究。E-mail:286741140@qq.com

    石健,1971年生,男,回族,遼寧法庫人,副教授,博士,主要從事固體力學(xué)及斷裂實(shí)驗(yàn)方法的研究。E-mail:jianshia@swpu.edu.cn

    趙莉,1990年生,女,漢族,四川達(dá)州人,碩士研究生,主要從事桿管柱斷裂力學(xué)方面的研究。E-mail:993945753@qq.com

    10.11885/j.issn.1674-5086.2013.12.09.01

    1674-5086(2014)03-0157-08

    TE921

    A

    2013–12–09 < class="emphasis_bold">網(wǎng)絡(luò)出版時(shí)間:

    時(shí)間:2014–05–21

    國家自然科學(xué)基金重點(diǎn)項(xiàng)目“氣體鉆井技術(shù)基礎(chǔ)研究”(51134004);國家自然科學(xué)基金青年項(xiàng)目“應(yīng)力波作用下鉆具表面三維裂紋的動態(tài)擴(kuò)展研究”(11102173)。

    猜你喜歡
    鉆柱彈塑性井口
    140MPa井口壓裂四通管道沖蝕分析
    矮塔斜拉橋彈塑性地震響應(yīng)分析
    彈塑性分析在超高層結(jié)構(gòu)設(shè)計(jì)中的應(yīng)用研究
    江西建材(2018年4期)2018-04-10 12:36:52
    自適應(yīng)BPSK在井下鉆柱聲波傳輸中的應(yīng)用
    CO2驅(qū)油系統(tǒng)的井口預(yù)警系統(tǒng)設(shè)計(jì)
    水平段鉆柱失穩(wěn)后自鎖分析*
    動載荷作用下冪硬化彈塑性彎曲裂紋塑性區(qū)
    氣田開發(fā)試采井口流程探討
    井底之蛙
    結(jié)構(gòu)動力彈塑性與倒塌分析(Ⅱ)——SAP2ABAQUS接口技術(shù)、開發(fā)與驗(yàn)證
    地震研究(2014年1期)2014-02-27 09:29:55
    少妇被粗大的猛进出69影院 | 熟女人妻精品中文字幕| 少妇人妻一区二区三区视频| 寂寞人妻少妇视频99o| www.av在线官网国产| 欧美高清成人免费视频www| 国产又色又爽无遮挡免| av黄色大香蕉| 老女人水多毛片| 国语对白做爰xxxⅹ性视频网站| 免费观看a级毛片全部| 熟女电影av网| 国产精品一区二区三区四区免费观看| 精品一品国产午夜福利视频| 久久久久久久国产电影| 99热全是精品| 成人毛片a级毛片在线播放| 王馨瑶露胸无遮挡在线观看| av在线观看视频网站免费| 黄色一级大片看看| 日韩欧美一区视频在线观看 | 精品少妇内射三级| av.在线天堂| 一级黄片播放器| videos熟女内射| av在线播放精品| 日本猛色少妇xxxxx猛交久久| 国产精品三级大全| 亚洲经典国产精华液单| 成人亚洲精品一区在线观看| 99re6热这里在线精品视频| 久久精品国产a三级三级三级| 国产极品天堂在线| 成人黄色视频免费在线看| 国产淫语在线视频| 一级毛片我不卡| 免费大片18禁| 久久久久久伊人网av| 免费观看av网站的网址| 高清在线视频一区二区三区| 亚洲精品日本国产第一区| 女人精品久久久久毛片| 80岁老熟妇乱子伦牲交| 美女内射精品一级片tv| 人人妻人人澡人人看| 欧美日韩av久久| 插逼视频在线观看| 人妻少妇偷人精品九色| 日本wwww免费看| 国产淫语在线视频| 三级国产精品片| 国产在线男女| 日本欧美视频一区| 国产亚洲最大av| 精品久久久噜噜| 日韩三级伦理在线观看| 麻豆精品久久久久久蜜桃| 97超碰精品成人国产| 一区二区三区四区激情视频| 日日撸夜夜添| 熟妇人妻不卡中文字幕| 亚洲国产精品一区三区| 国产成人精品福利久久| 久久ye,这里只有精品| 极品人妻少妇av视频| 亚洲天堂av无毛| 九草在线视频观看| 国产高清不卡午夜福利| 伦理电影大哥的女人| kizo精华| 日本黄色日本黄色录像| 午夜福利,免费看| 久久久久久久亚洲中文字幕| 精品亚洲乱码少妇综合久久| 中文字幕人妻熟人妻熟丝袜美| 一本一本综合久久| 大片电影免费在线观看免费| 亚洲真实伦在线观看| 内射极品少妇av片p| 3wmmmm亚洲av在线观看| 日本黄大片高清| 蜜臀久久99精品久久宅男| 又大又黄又爽视频免费| 最黄视频免费看| av播播在线观看一区| 午夜日本视频在线| 欧美精品一区二区大全| 久久99精品国语久久久| 婷婷色麻豆天堂久久| 制服丝袜香蕉在线| 国产永久视频网站| 欧美精品一区二区大全| 亚洲av成人精品一二三区| av又黄又爽大尺度在线免费看| 久久综合国产亚洲精品| 亚洲国产精品一区二区三区在线| 国产爽快片一区二区三区| 亚洲成人手机| 亚洲欧美日韩东京热| 午夜福利在线观看免费完整高清在| 边亲边吃奶的免费视频| 少妇人妻一区二区三区视频| 777米奇影视久久| 日本wwww免费看| 国产 一区精品| 欧美bdsm另类| 亚洲av男天堂| 日韩不卡一区二区三区视频在线| 街头女战士在线观看网站| 人妻系列 视频| 国产黄色免费在线视频| 成人亚洲精品一区在线观看| 麻豆成人av视频| 性高湖久久久久久久久免费观看| av在线老鸭窝| 国产精品.久久久| 日韩欧美 国产精品| 久久国产乱子免费精品| 搡女人真爽免费视频火全软件| 久久久精品94久久精品| 国产极品粉嫩免费观看在线 | 啦啦啦中文免费视频观看日本| 综合色丁香网| 久久青草综合色| av国产久精品久网站免费入址| 精品一区二区三区视频在线| 香蕉精品网在线| 日本欧美视频一区| 97超碰精品成人国产| 中文字幕人妻丝袜制服| av播播在线观看一区| 在线观看国产h片| 久久国产亚洲av麻豆专区| 国产美女午夜福利| 熟妇人妻不卡中文字幕| a级一级毛片免费在线观看| 欧美人与善性xxx| .国产精品久久| 久久久午夜欧美精品| 色吧在线观看| 丰满迷人的少妇在线观看| 欧美 亚洲 国产 日韩一| 国产精品伦人一区二区| 免费人成在线观看视频色| 久久久久久久久久久免费av| 亚洲真实伦在线观看| 午夜福利在线观看免费完整高清在| 五月伊人婷婷丁香| 国产精品福利在线免费观看| 深夜a级毛片| 丰满人妻一区二区三区视频av| 97在线视频观看| 男男h啪啪无遮挡| 久久午夜福利片| 亚洲欧美清纯卡通| av女优亚洲男人天堂| 99热全是精品| 成年av动漫网址| 久久人人爽人人爽人人片va| 久久女婷五月综合色啪小说| 高清黄色对白视频在线免费看 | 国产亚洲午夜精品一区二区久久| 老司机亚洲免费影院| 天堂俺去俺来也www色官网| 夜夜爽夜夜爽视频| 人人妻人人看人人澡| 全区人妻精品视频| 日韩中字成人| 有码 亚洲区| 亚洲va在线va天堂va国产| 丝瓜视频免费看黄片| 久久久国产欧美日韩av| 久久亚洲国产成人精品v| 黄色一级大片看看| 伊人亚洲综合成人网| 国产精品秋霞免费鲁丝片| av免费在线看不卡| 亚洲国产精品一区三区| 国产老妇伦熟女老妇高清| 日本欧美视频一区| 免费观看av网站的网址| 日韩伦理黄色片| 高清av免费在线| 亚洲成人手机| 精品人妻熟女av久视频| 人妻少妇偷人精品九色| 在线观看www视频免费| 亚洲美女搞黄在线观看| 伊人久久精品亚洲午夜| 久久久久久久国产电影| 亚洲欧美清纯卡通| 我要看黄色一级片免费的| 狂野欧美激情性bbbbbb| 美女内射精品一级片tv| 中文乱码字字幕精品一区二区三区| 蜜桃久久精品国产亚洲av| 男人舔奶头视频| 亚洲情色 制服丝袜| 视频中文字幕在线观看| 交换朋友夫妻互换小说| 高清午夜精品一区二区三区| 日日摸夜夜添夜夜添av毛片| 丝袜脚勾引网站| 免费在线观看成人毛片| 成人黄色视频免费在线看| 精品国产一区二区久久| 有码 亚洲区| 下体分泌物呈黄色| 亚洲三级黄色毛片| 男女免费视频国产| 欧美日韩视频精品一区| 我要看日韩黄色一级片| 在线看a的网站| 亚洲欧美精品自产自拍| 99久久精品热视频| 视频区图区小说| 午夜激情久久久久久久| 黑人高潮一二区| 久久ye,这里只有精品| 三级国产精品欧美在线观看| 久久精品国产自在天天线| 夫妻午夜视频| av女优亚洲男人天堂| 国产成人一区二区在线| 国产女主播在线喷水免费视频网站| 国产亚洲91精品色在线| 国产精品99久久99久久久不卡 | 中国国产av一级| 妹子高潮喷水视频| 看免费成人av毛片| 亚洲怡红院男人天堂| 国产免费一区二区三区四区乱码| 中文字幕免费在线视频6| 国产精品.久久久| 99久久中文字幕三级久久日本| 老司机亚洲免费影院| 成人美女网站在线观看视频| 人妻夜夜爽99麻豆av| 亚洲国产av新网站| 能在线免费看毛片的网站| 国产乱人偷精品视频| av网站免费在线观看视频| 国国产精品蜜臀av免费| 欧美xxxx性猛交bbbb| 两个人免费观看高清视频 | 美女福利国产在线| 熟妇人妻不卡中文字幕| 精品酒店卫生间| av网站免费在线观看视频| 欧美精品国产亚洲| 一区二区三区免费毛片| 欧美丝袜亚洲另类| 国产精品麻豆人妻色哟哟久久| 国产成人精品福利久久| 久久精品国产自在天天线| 久久国产亚洲av麻豆专区| 亚洲国产成人一精品久久久| 欧美激情极品国产一区二区三区 | 亚洲三级黄色毛片| 成人毛片a级毛片在线播放| 精品久久久久久久久av| 熟女人妻精品中文字幕| h视频一区二区三区| 亚洲激情五月婷婷啪啪| 国产日韩一区二区三区精品不卡 | 91精品国产国语对白视频| 又爽又黄a免费视频| 伦理电影大哥的女人| 在线观看人妻少妇| 美女脱内裤让男人舔精品视频| 高清av免费在线| 国产精品久久久久久久久免| 九九在线视频观看精品| 欧美少妇被猛烈插入视频| 秋霞在线观看毛片| 国产在线免费精品| 国产精品国产三级国产av玫瑰| 欧美xxxx性猛交bbbb| 久久精品国产亚洲av涩爱| 最黄视频免费看| 大片免费播放器 马上看| videos熟女内射| 日日啪夜夜撸| 精品国产国语对白av| 免费人成在线观看视频色| 免费高清在线观看视频在线观看| 亚洲,一卡二卡三卡| tube8黄色片| 伊人亚洲综合成人网| 少妇猛男粗大的猛烈进出视频| 成人美女网站在线观看视频| av在线观看视频网站免费| 国精品久久久久久国模美| 午夜影院在线不卡| 日产精品乱码卡一卡2卡三| 久久国产亚洲av麻豆专区| 久久精品熟女亚洲av麻豆精品| 国产精品福利在线免费观看| 欧美日韩国产mv在线观看视频| 亚洲第一av免费看| 中文天堂在线官网| 亚洲欧洲日产国产| 男人狂女人下面高潮的视频| 精品国产露脸久久av麻豆| 国产精品秋霞免费鲁丝片| 秋霞在线观看毛片| 黄色日韩在线| 有码 亚洲区| 久热久热在线精品观看| 欧美日本中文国产一区发布| 日韩av在线免费看完整版不卡| 久久久国产一区二区| 国产中年淑女户外野战色| 日日啪夜夜爽| 日日撸夜夜添| 亚洲欧美中文字幕日韩二区| 建设人人有责人人尽责人人享有的| 在线观看免费视频网站a站| a级片在线免费高清观看视频| 精品卡一卡二卡四卡免费| 99久国产av精品国产电影| 日本91视频免费播放| 两个人的视频大全免费| 2018国产大陆天天弄谢| 国产成人aa在线观看| 午夜福利,免费看| 亚洲精品色激情综合| 在线观看www视频免费| 高清黄色对白视频在线免费看 | 日本黄色片子视频| 日日摸夜夜添夜夜添av毛片| 亚洲真实伦在线观看| 一区二区三区乱码不卡18| 在线观看免费视频网站a站| 国内揄拍国产精品人妻在线| 亚洲国产精品成人久久小说| 国产在线免费精品| 中文精品一卡2卡3卡4更新| 男人爽女人下面视频在线观看| 久久久国产欧美日韩av| 精品国产一区二区三区久久久樱花| 国产一区二区在线观看av| 青青草视频在线视频观看| av播播在线观看一区| 日本与韩国留学比较| 国产国拍精品亚洲av在线观看| 极品教师在线视频| 你懂的网址亚洲精品在线观看| 少妇裸体淫交视频免费看高清| 中文在线观看免费www的网站| 深夜a级毛片| 麻豆成人午夜福利视频| 麻豆成人av视频| 国产极品粉嫩免费观看在线 | 国产男女内射视频| 91精品国产九色| 插逼视频在线观看| 日韩三级伦理在线观看| 男人和女人高潮做爰伦理| 男女无遮挡免费网站观看| h日本视频在线播放| 久久久久久久大尺度免费视频| 欧美高清成人免费视频www| 少妇丰满av| 丝袜脚勾引网站| 国产熟女欧美一区二区| 日本av手机在线免费观看| 丝袜喷水一区| 免费高清在线观看视频在线观看| 国产高清有码在线观看视频| 免费人成在线观看视频色| 最黄视频免费看| 精品一区二区三区视频在线| 久久久久国产网址| 男女国产视频网站| 午夜影院在线不卡| 国产午夜精品久久久久久一区二区三区| 精品少妇久久久久久888优播| 韩国高清视频一区二区三区| 黄色日韩在线| 国产一区亚洲一区在线观看| 国产精品国产三级国产专区5o| 免费看光身美女| 一本色道久久久久久精品综合| 亚洲精品日本国产第一区| 亚洲人成网站在线观看播放| 亚洲伊人久久精品综合| 我要看黄色一级片免费的| 欧美97在线视频| 日韩中字成人| 国产日韩欧美视频二区| 波野结衣二区三区在线| 国产欧美另类精品又又久久亚洲欧美| 久久精品久久久久久久性| 欧美最新免费一区二区三区| 2018国产大陆天天弄谢| 99热这里只有是精品在线观看| 在线观看免费视频网站a站| 久久久久精品久久久久真实原创| 亚洲欧美清纯卡通| 尾随美女入室| 日本vs欧美在线观看视频 | 美女视频免费永久观看网站| 日韩成人av中文字幕在线观看| 狠狠精品人妻久久久久久综合| 人妻夜夜爽99麻豆av| 嫩草影院入口| 男人添女人高潮全过程视频| 国产精品偷伦视频观看了| 人妻人人澡人人爽人人| 国产成人freesex在线| 久久久精品94久久精品| 亚洲成人一二三区av| 久久女婷五月综合色啪小说| 男女国产视频网站| 女的被弄到高潮叫床怎么办| 看十八女毛片水多多多| 国产伦精品一区二区三区四那| 午夜福利,免费看| 亚洲国产av新网站| 极品教师在线视频| 伊人亚洲综合成人网| 成人特级av手机在线观看| 国产深夜福利视频在线观看| 一区二区三区精品91| 色婷婷久久久亚洲欧美| 亚洲中文av在线| 五月玫瑰六月丁香| 在线亚洲精品国产二区图片欧美 | 亚洲欧美日韩东京热| 久久久精品94久久精品| 纯流量卡能插随身wifi吗| 99re6热这里在线精品视频| 亚洲欧美日韩另类电影网站| tube8黄色片| 亚洲av国产av综合av卡| av国产久精品久网站免费入址| 两个人免费观看高清视频 | 国产精品久久久久成人av| 在线精品无人区一区二区三| 我的女老师完整版在线观看| 欧美成人精品欧美一级黄| 成人亚洲精品一区在线观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲内射少妇av| 亚洲精品一区蜜桃| 男女啪啪激烈高潮av片| 一区二区三区免费毛片| 99国产精品免费福利视频| 啦啦啦在线观看免费高清www| 嘟嘟电影网在线观看| 亚洲,一卡二卡三卡| 97在线人人人人妻| 国产一区二区三区av在线| 成人无遮挡网站| 成人美女网站在线观看视频| 在线天堂最新版资源| 亚洲精品色激情综合| 亚洲精品第二区| 国产成人精品久久久久久| 91精品国产国语对白视频| 欧美成人午夜免费资源| 亚洲av成人精品一二三区| 亚洲精品456在线播放app| 两个人的视频大全免费| 夜夜骑夜夜射夜夜干| 日韩在线高清观看一区二区三区| 18禁裸乳无遮挡动漫免费视频| 国产一级毛片在线| 男女啪啪激烈高潮av片| 一级毛片我不卡| 国产成人freesex在线| 狂野欧美激情性xxxx在线观看| 精品亚洲成国产av| 制服丝袜香蕉在线| 大片免费播放器 马上看| 午夜免费男女啪啪视频观看| 在线观看三级黄色| 伊人久久国产一区二区| 自拍偷自拍亚洲精品老妇| 一区二区三区乱码不卡18| 日韩熟女老妇一区二区性免费视频| 免费大片18禁| 青春草视频在线免费观看| 亚洲怡红院男人天堂| 妹子高潮喷水视频| 亚洲激情五月婷婷啪啪| 在线 av 中文字幕| 99精国产麻豆久久婷婷| 啦啦啦视频在线资源免费观看| 色5月婷婷丁香| 日韩av免费高清视频| av在线老鸭窝| 日本vs欧美在线观看视频 | 国产在线男女| 啦啦啦在线观看免费高清www| 久久久久久久亚洲中文字幕| 亚洲美女视频黄频| 晚上一个人看的免费电影| 日产精品乱码卡一卡2卡三| 亚洲国产精品成人久久小说| 日本-黄色视频高清免费观看| 啦啦啦视频在线资源免费观看| 色5月婷婷丁香| 18禁裸乳无遮挡动漫免费视频| 国产精品99久久99久久久不卡 | 成人无遮挡网站| 成人综合一区亚洲| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久久久免| 日日啪夜夜撸| 亚洲人与动物交配视频| 制服丝袜香蕉在线| 亚洲色图综合在线观看| 美女中出高潮动态图| 99久久精品热视频| 性色av一级| 日本-黄色视频高清免费观看| 久久国产亚洲av麻豆专区| 国产成人一区二区在线| 久久久精品94久久精品| 久久国产精品男人的天堂亚洲 | 中文字幕制服av| 欧美日韩在线观看h| 精品人妻熟女av久视频| 三级经典国产精品| 一级二级三级毛片免费看| 97超碰精品成人国产| av天堂久久9| 人人妻人人爽人人添夜夜欢视频 | 欧美日韩国产mv在线观看视频| 熟女电影av网| 日韩欧美精品免费久久| 国产精品.久久久| 午夜福利视频精品| 伦理电影大哥的女人| 亚洲成人手机| 夫妻性生交免费视频一级片| 久久97久久精品| 欧美人与善性xxx| 日韩一本色道免费dvd| 一区在线观看完整版| 日韩 亚洲 欧美在线| 亚洲精品自拍成人| 久久精品国产亚洲网站| 国产午夜精品久久久久久一区二区三区| 麻豆精品久久久久久蜜桃| 国产精品国产三级专区第一集| 狂野欧美白嫩少妇大欣赏| 久久久久久久久久久久大奶| 乱系列少妇在线播放| 亚洲一区二区三区欧美精品| 亚洲国产色片| av网站免费在线观看视频| 精品一区二区三区视频在线| 伊人亚洲综合成人网| 黄色视频在线播放观看不卡| 在线观看人妻少妇| 欧美日韩综合久久久久久| 欧美日韩国产mv在线观看视频| av福利片在线观看| 日韩免费高清中文字幕av| 一级毛片黄色毛片免费观看视频| 日韩欧美精品免费久久| 久久久久久久久久久免费av| 26uuu在线亚洲综合色| 国产欧美另类精品又又久久亚洲欧美| 乱人伦中国视频| 草草在线视频免费看| 中文在线观看免费www的网站| 亚洲自偷自拍三级| 人人妻人人看人人澡| 国产精品女同一区二区软件| 黄色配什么色好看| 日日爽夜夜爽网站| 色94色欧美一区二区| 高清欧美精品videossex| 日本免费在线观看一区| 免费观看性生交大片5| 国产黄片美女视频| 精品久久久精品久久久| 成人亚洲欧美一区二区av| 亚洲欧美清纯卡通| 国产精品一区二区性色av| 午夜激情久久久久久久| 91精品伊人久久大香线蕉| 欧美变态另类bdsm刘玥| 极品人妻少妇av视频| 三级国产精品欧美在线观看| 美女cb高潮喷水在线观看| 制服丝袜香蕉在线| 人妻人人澡人人爽人人| 黄色一级大片看看| 好男人视频免费观看在线| 精品酒店卫生间| 免费观看性生交大片5| 最近中文字幕2019免费版| 五月天丁香电影| 免费人妻精品一区二区三区视频| 晚上一个人看的免费电影| 国产成人免费观看mmmm| 国产日韩欧美在线精品| 国产伦精品一区二区三区视频9| 99热这里只有是精品50| 少妇的逼好多水| 97超碰精品成人国产| 熟女电影av网| 各种免费的搞黄视频| 国产色婷婷99| 久久久精品94久久精品| 欧美激情国产日韩精品一区| 国产成人免费无遮挡视频| 免费观看在线日韩| 韩国高清视频一区二区三区| 九九久久精品国产亚洲av麻豆|