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

    含表面雙裂紋的鉆桿管體斷裂性能研究

    2024-06-19 00:00:00張智亮林瑩瑩錢光陽李文軍李安軍

    收稿日期:2022-08-12""" 修回日期:2023-02-02

    基金項(xiàng)目:國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃資助項(xiàng)目(No.2012CB215004);國家科技重大專項(xiàng)課題資助項(xiàng)目(No.2011ZX05021-001);中國博士后科學(xué)基金資助項(xiàng)目(No.2019M663865)

    通信作者:張智亮。E-mail:280238724@qq.com

    引用格式:

    張智亮,林瑩瑩,錢光陽,等. 含表面雙裂紋的鉆桿管體斷裂性能研究[J].應(yīng)用力學(xué)學(xué)報(bào),2024,41(3):516-523.

    ZHANG Zhiliang,LIN Yingying,QIAN Guangyang,et al. Study on fracture behavior of drill pipe body with surface double cracks[J].Chinese journal of applied mechanics,2024,41(3):516-523.

    文章編號:1000-4939(2024)03-0516-08

    摘" 要:鉆桿刺漏是鉆桿在服役過程中最常見的失效形式,而鉆桿表面裂紋的擴(kuò)展是鉆桿發(fā)生刺漏的直接原因。為更有效地評估鉆桿表面含有多條裂紋時(shí)的安全性能,以某超深直井G105鉆桿井口段為研究對象,建立管體外表面雙裂紋有限元模型,分析了拉扭載荷下裂紋間距、次裂紋尺寸以及走向?qū)Νh(huán)向主裂紋前緣J積分的影響。結(jié)果表明:軸向排列的環(huán)向次裂紋是相對最安全的,可大幅消減鉆桿主裂紋前緣的J積分;環(huán)向排列的環(huán)向次裂紋相對最危險(xiǎn),可大幅增大鉆桿主裂紋前緣的J積分;軸向次裂紋可增大鉆桿主裂紋前緣或部分前緣的J積分,但增幅相對環(huán)向排列的環(huán)向次裂紋來說要微弱得多。次裂紋對主裂紋前緣J積分的影響強(qiáng)度隨兩裂紋間距的增大或次裂紋尺寸的減小而減小。因此在對含表面雙裂紋的鉆桿進(jìn)行分析時(shí),在次裂紋尺寸較大且距離較近的情況下應(yīng)考慮次裂紋的影響;在次裂紋尺寸較小或距離較遠(yuǎn)時(shí)則可忽略次裂紋,將模型簡化為單裂紋鉆桿模型進(jìn)行分析。

    關(guān)鍵詞:鉆桿;雙裂紋;J積分;安全評估

    中圖分類號:TE921" 文獻(xiàn)標(biāo)志碼:A

    DOI:10.11776/j.issn.1000-4939.2024.03.004

    Study on fracture behavior of drill pipe body with surface double cracks

    ZHANG Zhiliang1,2,LIN Yingying1,QIAN Guangyang3,LI Wenjun4," LI Anjun5

    (1.School of Mechatronic Engineering,Southwest Petroleum University,610500

    Chengdu,China;

    2.Sichuan Zhongliyuan Construction Engineering Co.,Ltd.,610044 Chengdu,China;

    3.Southwestern Pipeline Co.,Ltd.,State Pipeline Network Group,610213 Chengdu,China;

    4.Petro China Urumqi Petrochemical Company,831400 Urumqi,China;

    5.Western Pipeline Co.,Ltd.,State Pipeline Network Group,831400 Urumqi,China)

    Abstract:Drilling pipe leakage is the most common failure form of drilling pipe in service,and the expansion of surface crack of drilling pipe is the direct cause of drilling pipe leakage.In order to evaluate the safety performance of drill pipe with multiple cracks on the surface,a finite element model of double cracks on the surface of drill pipe was established by taking the wellhead section of G105 drill pipe of an ultra-deep vertical well as the research object.The effects of crack spacing,crack size and strike on the J-integral at the leading edge of the circular main crack under tensile and torsional loading were analyzed.The results show that the circumferential secondary cracks arranged in the axial direction are relatively the safest,which can greatly reduce the J-integral of the leading edge of the main crack.The circumferential secondary cracks arranged in the circumferential direction are the most dangerous,which can greatly increase the J-integral of the leading edge of the main crack.The axial secondary cracks can increase the J-integral of the leading edge or part of the leading edge of the main crack in the drill pipe,but the increase is much weaker than that of the circumferential secondary cracks.The influence intensity of secondary cracks on the leading edge J-integral of the main crack decreases with the increase of the distance between two cracks and the decrease of the size of secondary cracks.Therefore,the influence of secondary cracks should be considered when the size of secondary cracks is large and the distance between them is relatively close.When the size of the secondary crack is small or the distance is long,the secondary crack can be ignored,and the model is simplified to a single-crack drill pipe model for analysis.

    Key words:drill pipe; double crack; J-integral; safety assessment

    在鉆井和起下鉆過程中,鉆桿井口處受力狀態(tài)復(fù)雜且承受最大拉力和扭矩,是鉆桿失效的危險(xiǎn)位置[1]。據(jù)統(tǒng)計(jì),目前大部分鉆桿失效事故是由鉆桿刺漏引起的[2],而鉆桿裂紋的擴(kuò)展是導(dǎo)致鉆桿刺漏的直接原因[3]。

    目前國內(nèi)外學(xué)者針對鉆桿裂紋擴(kuò)展規(guī)律進(jìn)行了大量研究。雒設(shè)計(jì)等[4]對S135鉆桿鋼的疲勞裂紋在不同應(yīng)力比下的擴(kuò)展行為進(jìn)行了研究,發(fā)現(xiàn)隨著應(yīng)力比的增加,疲勞裂紋擴(kuò)展呈速率曲線向左上方移動的趨勢。尹茂等[5]對含不同幾何參數(shù)裂紋的鉆桿模型進(jìn)行了研究,發(fā)現(xiàn)

    Ⅰ型裂紋是鉆桿裂紋在拉扭載荷作用下擴(kuò)展的主要推動力。趙廣慧等[6-8]對鉆桿本體、接頭以及加厚過渡段的表面裂紋進(jìn)行了研究,得到了不同載荷比和不同深度比條件下的裂紋臨界形狀比。LEMMA等[9]為預(yù)測組合載荷作用下鉆桿疲勞壽命,采用模糊系統(tǒng)和歐拉方法等確定了載荷大小及變化、裂紋形狀和幾何參數(shù)對裂紋擴(kuò)展規(guī)律的影響。DAO等[10]研究了空心圓柱在旋轉(zhuǎn)彎曲和拉伸作用下圓周半橢圓表面裂紋的應(yīng)力強(qiáng)度因子,利用三維有限元模型建立了3種荷載情況下適用于多種結(jié)構(gòu)的應(yīng)力強(qiáng)度因子數(shù)據(jù)庫。LIU等[11]利用有限元分析軟件模擬了無機(jī)械損傷和有初始裂紋兩種不同情況下鉆桿體的應(yīng)力分布,發(fā)現(xiàn)最大應(yīng)力隨初始裂紋深度的增加而增加,最大應(yīng)力位置由裂紋徑向根部向裂紋周向根部轉(zhuǎn)移。OJANOMARE等 [12]利用多參考狀態(tài)權(quán)函數(shù)法對拉彎載荷下的鉆桿裂紋進(jìn)行了研究,得到了鉆桿半橢圓內(nèi)裂紋前沿的應(yīng)力強(qiáng)度因子,并在MATLAB代碼中實(shí)現(xiàn)了多參考狀態(tài)權(quán)函數(shù)法。

    綜上所述,國內(nèi)外學(xué)者對于工況、載荷條件、初始裂紋尺寸以及裂紋所處位置對鉆桿裂紋擴(kuò)展規(guī)律的影響進(jìn)行了大量研究,但這些研究通?;趩瘟鸭y模型,對于含雙裂紋乃至多裂紋的鉆桿模型研究較少,而工程實(shí)際中鉆桿裂紋通常是兩條甚至多條同時(shí)存在的。本研究在前人研究的基礎(chǔ)上,以某超深直井井口附近鉆桿實(shí)際尺寸以及載荷條件為依托,對雙裂紋鉆桿中不同走向及分布的次裂紋對主裂紋擴(kuò)展規(guī)律的影響進(jìn)行研究,為鉆桿的安全性評估提供可靠依據(jù)。

    1" 理論基礎(chǔ)

    1.1" 應(yīng)力強(qiáng)度因子

    應(yīng)力強(qiáng)度因子K是表征裂紋斷裂性能的重要參數(shù)。以裂紋中心為原點(diǎn),創(chuàng)建局部三維直角坐標(biāo)系(圖1),z軸指向裂紋深度方向,x軸沿裂紋延長線方向,y軸指向裂紋張開方向。應(yīng)用位移法,裂紋前緣上某點(diǎn)的應(yīng)力強(qiáng)度因子K可以通過局部坐標(biāo)系下兩個(gè)裂紋面的節(jié)點(diǎn)位移表示為[13]

    KⅠ=2G2π1/21+κΔvr

    KⅡ=2G2π1/21+κΔur

    KⅢ=2G2π1/2Δwr(1)

    式中:KⅠ、KⅡ、KⅢ 分別為Ⅰ、Ⅱ、Ⅲ型應(yīng)力強(qiáng)度因子;G為剪切模量;Δu、Δv、Δw分別為兩裂紋面在x軸、y軸以及z軸方向上的相對節(jié)點(diǎn)位移;平面應(yīng)力條件下κ=3-μ/(1+μ),平面應(yīng)變條件下κ=3-4μ,μ為泊松比;r為裂紋前緣任意一點(diǎn)到裂紋尖端的距離。

    Fig.1 3D crack schematic diagram

    1.2" J積分

    J積分是彈塑性斷裂力學(xué)中表征裂紋尖端應(yīng)力應(yīng)變大小的斷裂參數(shù),由外載荷通過施力點(diǎn)位移對試樣所做的形變功率給出[14]。J積分定義為

    J=∫Γωdy-Txuxx+Tyuyyds(2)

    式中:Γ為圍繞裂尖的任意回路(J積分與積分路徑無關(guān)),如圖2所示;ω為應(yīng)變能密度;Tx、Ty為作用在Γ單位長度上的分力;ux、uy為位移分量;ds為積分弧長。

    Fig.2 Schematic diagram of integral loop

    組合載荷作用下復(fù)合型裂紋的等效應(yīng)力強(qiáng)度因子K與J積分的關(guān)系式為

    K=JE-(3)

    式中,平面應(yīng)力條件下E-=E,平面應(yīng)變條件下E-=E/(1-μ2),E為材料彈性模量。

    2" 含裂紋鉆桿有限元模型

    2.1" 鉆桿力學(xué)模型

    本研究選取某超深直井的G105鉆桿井口附近部分為研究對象。鉆桿內(nèi)徑d=117.81mm,壁厚t=9.19mm。根據(jù)鉆桿的實(shí)際工況條件,鉆桿上端面固定,鉆桿下端面承受扭矩N=10kN·m和拉伸載荷P=1050kN。鉆桿力學(xué)模型與裂紋截面如圖3所示。

    鉆桿材料力學(xué)特性如表1所示[15-16]。

    2.2" 有限元模型

    2.2.1" 單裂紋鉆桿有限元模型

    為消除固定端面對應(yīng)力強(qiáng)度因子的影響,模型長度取1000mm,鉆桿中部外表面有一半橢圓裂紋,裂紋深度(半短軸)為a,裂紋長度(長軸)為2c,a/t=0.2,a/c=0.8。裂紋前緣選用1/4奇異單元代替;裂尖加密區(qū)域最內(nèi)層圓周半徑0.01mm采用楔形體網(wǎng)格,其外半徑0.25mm的環(huán)形加密區(qū)采用單精度六面體網(wǎng)格,圓周局部種子數(shù)目為18;裂紋周圍加密區(qū)選用四面體自由網(wǎng)格,其他區(qū)域選用六面體自由網(wǎng)格。在模型下端面分別施加扭矩和拉伸載荷。單裂紋鉆桿有限元模型及其網(wǎng)格劃分如圖4所示。

    2.2.2" 雙裂紋鉆桿有限元模型

    除裂紋數(shù)量為兩條且模型下端面同時(shí)施加扭矩以及拉伸載荷外,雙裂紋鉆桿有限元模型的建模過程與單裂紋模型基本相同。定義位于鉆桿中部或左側(cè)的環(huán)向裂紋為主裂紋,兩裂紋中心連線與鉆桿軸線平行則稱軸向排列,兩裂紋中心在鉆桿同一周線上則稱環(huán)向排列。對于軸向排列的雙裂紋組合,裂紋間距以兩裂紋中心連線長度f表示,對于環(huán)向排列的雙裂紋組合,裂紋間距以兩裂紋所處截面圓心與兩裂紋中心連線的夾角β表示。兩裂紋的走向及分布如圖5所示。

    2.3" 有限元模型驗(yàn)證

    此處驗(yàn)證參數(shù)選用無因次應(yīng)力強(qiáng)度因子κ。

    κ=KΙσπa(4)

    σ=Pπ4D2-d2×106 (拉伸載荷下)

    σ=τ=NπD3161-dD4×106 (扭轉(zhuǎn)載荷下)(5)

    式中:KⅠ可根據(jù)J積分由式(3)計(jì)算得到;σ為裂紋所處截面上的名義應(yīng)力;D為鉆桿外徑。

    由于缺少含裂紋鉆桿斷裂參數(shù)的實(shí)驗(yàn)數(shù)據(jù),通過將有限元模型計(jì)算結(jié)果與中國航空研究院編著的《應(yīng)力強(qiáng)度因子手冊》(增訂版)[17]中的參考值進(jìn)行對比,以此來驗(yàn)證有限元模型的可靠性。采用2.2.1中所述方法構(gòu)建鉆桿單裂紋有限元模型,取該模型裂紋前緣最深點(diǎn)C以及邊界點(diǎn)A的計(jì)算結(jié)果與參考值進(jìn)行對比,結(jié)果如表2所示。

    由表2可看出,有限元模型計(jì)算值與參考值的相對誤差在5%以內(nèi),因此可認(rèn)為該有限元模型是可靠的。

    3" 雙裂紋有限元模型計(jì)算結(jié)果分析

    為更直觀地表征次裂紋對主裂紋的影響,定義F為影響因子,用以表征次裂紋對主裂紋前緣某一點(diǎn)的影響程度,其表達(dá)式為

    F(θ)=J-J0J0" (6)

    式中:J表示次裂紋存在時(shí)主裂紋前緣該點(diǎn)處的J積分;J0表示無次裂紋存在時(shí)主裂紋前緣該點(diǎn)處的J積分;θ表示主裂紋前緣某點(diǎn)相對位置的參數(shù),且0°≤θ≤180°。在某一點(diǎn)當(dāng)Fgt;0時(shí),表明次裂紋增大了裂紋前緣該點(diǎn)處的J積分;而當(dāng)Flt;0時(shí),表明次裂紋減小了主裂紋前緣該點(diǎn)處的J積分。F絕對值|F|的大小即次裂紋對裂紋前緣影響強(qiáng)度的大小。

    3.1" 次裂紋走向及分布類型的影響

    本節(jié)主要研究次裂紋走向以及分布類型對主裂紋前緣F的影響。以a1、a2分別表示主裂紋和次裂紋的深度,以2c1、2c2分別表示主裂紋和次裂紋的裂紋長度。保持主裂紋與次裂紋尺寸相同,取裂紋形狀比a/c=a1/c1=a2/c2=0.9,裂紋半長c=c1=c2=3.335mm。4個(gè)模型計(jì)算結(jié)果如圖6所示。

    由圖6(a)可見,當(dāng)雙裂紋鉆桿存在軸向排列的環(huán)向次裂紋時(shí),主裂紋裂紋前緣的Flt;0,即此時(shí)次裂紋對主裂紋前緣的J積分有減小的作用。由圖6(b)可見,當(dāng)雙裂紋鉆桿存在環(huán)向排列的環(huán)向次裂紋時(shí),主裂紋裂紋前緣的Fgt;0,即此時(shí)次裂紋對主裂紋前緣的J積分有著增大的效果。由圖6(c)可見,當(dāng)雙裂紋鉆桿含有環(huán)向排列的軸向次裂紋時(shí),主裂紋裂紋前緣的Fgt;0,即此時(shí)次裂紋對主裂紋前緣的J積分有著增大的效果。由圖6(d)可見,當(dāng)雙裂紋鉆桿存在軸向排列的軸向次裂紋時(shí),主裂紋前緣的F值隨θ近似線性變化,且當(dāng)0°≤

    θ<100°時(shí)Fgt;0,即主裂紋前緣該處的J積分被增大;θ>100°時(shí)Flt;0,即主裂紋前緣該處的J積分被減小。

    綜合圖6可以看出,對于含有雙裂紋的鉆桿模型,當(dāng)兩裂紋間距較近時(shí),在次裂紋的影響下主裂紋前緣的F都是隨θ變化的,且在θ=15°或θ=165°附近達(dá)到極值;隨著兩裂紋間距逐漸增大,主裂紋前緣各點(diǎn)的F值隨θ變化也逐漸減小。在兩裂紋間距較近時(shí),對于軸向排列的次裂紋F隨θ近似線性減小的變化趨勢;對于環(huán)向排列的環(huán)向次裂紋F隨θ近似指數(shù)形式遞減。

    由圖7可見,當(dāng)鉆桿存在環(huán)向排列的次裂紋時(shí),主裂紋前緣的Fgt;0,但環(huán)向次裂紋下的主裂紋前緣的F始終大于軸向次裂紋下主裂紋前緣的F,且兩條曲線距離較遠(yuǎn)。這說明雖環(huán)向排列的軸向次裂紋和環(huán)向次裂紋都對主裂紋前緣的J有增強(qiáng)的作用,但在裂紋間距相同時(shí)環(huán)向次裂紋比軸向次裂紋產(chǎn)生的影響明顯更強(qiáng)烈。

    3.2" 兩裂紋間距的影響

    本節(jié)主要研究兩裂紋間距對主裂紋前緣F的影響。保持主裂紋與次裂紋尺寸相同,分別計(jì)算形狀比a/c=0.6、裂紋半長c=3.335mm,以及形狀比a/c=0.9、裂紋半長c= 3.335mm的兩組模型,計(jì)算結(jié)果如圖8所示。

    從圖8可以看出,不同的裂紋間距下主裂紋前緣的F值也不同。由圖8(a)可見,隨著距離f的增大主裂紋前緣的|F|逐漸減小,表明該次裂紋對主裂紋前緣J積分的影響在減弱,在f ≥32mm后|F|值逐漸趨近于0,表明此時(shí)該次裂紋對主裂紋的影響已經(jīng)非常微弱,在這種情況下可忽略次裂紋將雙裂紋鉆桿視作單裂紋鉆桿進(jìn)行研究。由圖8(b)可見,隨著β的增大主裂紋前緣的|F|逐漸減小,表明此時(shí)次裂紋對主裂紋前緣J積分的影響在減弱,在β≥35°后|F|值逐漸趨近于0,表明此時(shí)該次裂紋對主裂紋的影響已經(jīng)非常微弱,在這種情況下可忽略次裂紋將雙裂紋鉆桿視作單裂紋鉆桿進(jìn)行研究。由圖8(c)可見,隨著β的增大主裂紋前緣的|F|逐漸減小,表明此時(shí)次裂紋對主裂紋前緣J積分的影響在減弱,在β≥15°后|F|值逐漸趨近于0,表明此時(shí)該次裂紋對主裂紋的影響已經(jīng)非常微弱,在這種情況下可忽略次裂紋將雙裂紋鉆桿視作單裂紋鉆桿進(jìn)行研究。

    由圖8(d)可見,隨著f的增大主裂紋前緣的|F|逐漸減小,表明此時(shí)次裂紋對主裂紋前緣J積分的影響在減弱,在f ≥14mm后|F|值趨近于0,表明此時(shí)該次裂紋對主裂紋的影響已經(jīng)非常微弱,在這種情況下可忽略次裂紋將雙裂紋鉆桿視作單裂紋鉆桿進(jìn)行研究。

    3.3" 次裂紋尺寸的影響

    以3.1中的計(jì)算結(jié)果為基礎(chǔ),選取各模型下|F|最大時(shí)的裂紋間距,研究裂紋間距以及主裂紋尺寸一定時(shí),不同形狀比和深度比的次裂紋對主裂紋的影響因子。主裂紋尺寸a1/c1=0.9,a1=2mm,保持次裂紋深度不變,即a2=1.5mm,分別計(jì)算a2/c2=0.5,0.6,0.7,0.8,0.9,1.0時(shí),主裂紋前緣的J積分,并將其換算成|F|,取主裂紋前緣θ=15°時(shí)的|F|值進(jìn)行作圖;保持次裂紋形狀比不變,即a2/c2=0.9,分別計(jì)算a2/t=0.10,0.15,0.20,0.25,0.30時(shí),主裂紋前緣的J積分,并將其換算成|F|,取主裂紋前緣θ=15°時(shí)的|F|值進(jìn)行作圖。

    由圖9(a)可見,雙裂紋鉆桿主裂紋前緣|F|值隨次裂紋的形狀比增大呈指數(shù)形式遞減,即當(dāng)雙裂紋鉆桿的主裂紋尺寸、兩裂紋間距以及次裂紋深度不變時(shí),隨著次裂紋形狀比增大主裂紋前緣的J積分先急劇減小而后慢慢趨于平穩(wěn);由圖9(b)可見,雙裂紋鉆桿主裂紋前緣|F|值隨次裂紋的深度比增大而增大,即當(dāng)雙裂紋鉆桿的主裂紋尺寸、兩裂紋間距以及次裂紋形狀比不變時(shí),次裂紋對主裂紋前緣J積分的影響隨次裂紋的深度比增大而增強(qiáng)。

    由圖9可見,雙裂紋鉆桿主裂紋前緣|F|值隨次裂紋尺寸增大而增大,隨次裂紋尺寸的減小而減小。當(dāng)次裂紋深度比小于0.1時(shí),可忽略次裂紋的影響,將模型視為單裂紋鉆桿模型進(jìn)行研究。結(jié)合圖9與圖6可見,4種次裂紋中環(huán)向排列的環(huán)向次裂紋對主裂紋前緣J積分的增大作用最強(qiáng),在對雙裂紋鉆桿模型進(jìn)行研究時(shí),應(yīng)重視這類次裂紋的存在。

    4" 結(jié)" 論

    通過構(gòu)建雙裂紋鉆桿模型對次裂紋的走向及分布對環(huán)向主裂紋擴(kuò)展規(guī)律的影響進(jìn)行研究,得出以下結(jié)論。

    1)當(dāng)兩裂紋間距較小時(shí),次裂紋對主裂紋的影響隨裂紋角度θ變化而變化,并在θ=15°或θ=165°附近達(dá)到最大值。隨著兩裂紋間距的增大,主裂紋前緣的F隨裂紋角度θ的變化逐漸減小。

    2) 環(huán)向排列的次裂紋對主裂紋的J積分有增強(qiáng)作用,且當(dāng)裂紋中心連線夾角相同時(shí)環(huán)向次裂紋對主裂紋的影響明顯比軸向次裂紋更劇烈,因而在對鉆桿進(jìn)行安全性能評估時(shí)應(yīng)尤為注意環(huán)向排列的環(huán)向次裂紋。

    3) 軸向排列的環(huán)向次裂紋對主裂紋前緣的J積分有減弱的作用,但軸向排列的軸向次裂紋

    0°<θ<100°時(shí)對主裂紋前緣的J積分有增強(qiáng)的作用。

    4) 次裂紋對主裂紋的影響在兩裂紋間距較小時(shí)較為強(qiáng)烈,且隨兩裂紋間距的增大而減?。淮瘟鸭y對主裂紋的影響隨次裂紋尺寸的減小而減小。因而當(dāng)次裂紋尺寸較小或兩裂紋間距較大時(shí)可忽略次裂紋,將雙裂紋鉆桿簡化為單裂紋鉆桿模型。

    參考文獻(xiàn):

    [1]" 趙萬春,劉瑩.鉆井與完井工程[M].哈爾濱:哈爾濱工業(yè)大學(xué)出版社,2012.

    [2]" 練章華,駱發(fā)前,龔建文,等.塔里木油田鉆桿刺漏原因分析[J].鉆采工藝,2003,26(6):62-65.

    LIAN Zhanghua,LUO Faqian,GONG Jianwen,et al.Analysis of drill pipe washout causes in Tarim Oilfield[J].Drilling amp; production technology,2003,26(6):62-65(in Chinese).

    [3]" 趙鵬.石油鉆桿加厚過渡區(qū)三維有限元分析[J].寶鋼技術(shù),2006(3):13-17.

    ZHAO Peng.3D finite element analysis of upset transition area for drill pipes[J].Baosteel technology,2006(3):13-17(in Chinese).

    [4]" 雒設(shè)計(jì),趙康,王榮.不同應(yīng)力比下S135鉆桿鋼疲勞裂紋的擴(kuò)展行為[J].機(jī)械工程材料,2013,37(7):72-76.

    LUO Sheji,ZHAO Kang,WANG Rong.Fatigue cracks propagation behavior of S135 drill pipe steel at different stress ratios[J].Materials for mechanical engineering,2013,37(7):72-76(in Chinese).

    [5]" 尹茂,周井玲,許波兵.含損傷裂紋的鉆桿斷裂與疲勞分析[J].煤礦機(jī)械,2020,41(12):85-87.

    YIN Mao,ZHOU Jingling,XU Bobing.Fracture and fatigue analysis of drill pipe with damage crack[J].Coal mine machinery,2020,41(12):85-87(in Chinese).

    [6]" 趙廣慧,趙莉,石健,等.鉆桿接頭表面裂紋的斷裂性能研究[J].應(yīng)用力學(xué)學(xué)報(bào),2016,33(5):838-844.

    ZHAO Guanghui,ZHAO Li,SHI Jian,et al.Fracture behaviors of surface cracks in the toll joint[J].Chinese journal of applied mechanics,2016,33(5):838-844(in Chinese).

    [7]" 趙莉.鉆桿本體及接頭表面裂紋的斷裂性能研究[D].成都:西南石油大學(xué),2016.

    [8]" 趙廣慧,趙莉,唐松,等.鉆桿內(nèi)加厚過渡段表面裂紋的斷裂性能研究[J].機(jī)械強(qiáng)度,2017,39(4):920-926.

    ZHAO Guanghui,ZHAO Li,TANG Song,et al.Fracture characters of surface crack in the inner upset transitional section of drill pipe[J].Journal of mechanical strength,2017,39(4):920-926(in Chinese).

    [9]" LEMMA T A,NANJI P,GEBREMARIAM M A,et al.Remaining useful life prediction of drill string using fuzzy systems and cumulative damage theory[J].Key engineering materials,2019,796:145-154.

    [10]DAO N H,SELLAMI H.Stress intensity factors and fatigue growth of a surface crack in a drill pipe during rotary drilling operation[J].Engineering fracture mechanics,2012,96:626-640.

    [11]LIU Y G,LI F P,YE D P,et al.Simulation on stress distribution of drill pipes surface initial crack[J].Applied mechanics and materials,2012,184-185:460-464.

    [12]OJANOMARE C,CORNETTI P,ROMAGNOLI R,et al.Fatigue crack growth analysis of drill pipes during rotary drilling operations by the multiple reference state weight function approach[J].Engineering failure analysis,2017,74:11-34.

    [13]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(5):1141-1156.

    [14]酈正能.應(yīng)用斷裂力學(xué)[M].北京:北京航空航天大學(xué)出版社,2012.

    [15]API.Specification for drill pipe:API SPEC 5DP[S].Washington,DC,USA:API,2009.

    [16]蔣林,施太和,郭西水,等.石油管Ⅰ型裂紋斷裂韌度與屈服強(qiáng)度的相關(guān)性研究[J].重慶科技學(xué)院學(xué)報(bào)(自然科學(xué)版),2015,17(6):84-88.

    JIANG Lin,SHI Taihe,GUO Xishui,et al.Interrelated law between mode-I fracture toughness and yield strength of OCTG[J].Journal of Chongqing University of Science and Technology(natural sciences edition),2015,17(6):84-88(in Chinese).

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

    (編輯" 李坤璐)

    成人毛片60女人毛片免费| 日韩成人伦理影院| 久久久久国产网址| 纵有疾风起免费观看全集完整版 | 91久久精品国产一区二区三区| 中文字幕熟女人妻在线| av国产久精品久网站免费入址| 麻豆成人午夜福利视频| 国产精品久久电影中文字幕| 91精品一卡2卡3卡4卡| 熟女电影av网| 国产精品福利在线免费观看| 亚洲av福利一区| 两个人视频免费观看高清| 简卡轻食公司| 国产av一区在线观看免费| 在线播放无遮挡| 1000部很黄的大片| 亚洲怡红院男人天堂| 久久久久久大精品| 男人的好看免费观看在线视频| 晚上一个人看的免费电影| 看片在线看免费视频| 国产老妇女一区| 国产真实伦视频高清在线观看| 变态另类丝袜制服| 亚洲精品影视一区二区三区av| 桃色一区二区三区在线观看| 天堂网av新在线| 亚洲经典国产精华液单| 久久国产乱子免费精品| 欧美一区二区国产精品久久精品| 插阴视频在线观看视频| 特大巨黑吊av在线直播| 久久久欧美国产精品| 国产精品久久久久久久电影| 男人舔女人下体高潮全视频| 亚洲av电影在线观看一区二区三区 | 国产精品嫩草影院av在线观看| 精品99又大又爽又粗少妇毛片| 中文资源天堂在线| 一区二区三区高清视频在线| 国产成人91sexporn| 最近中文字幕高清免费大全6| 亚洲综合精品二区| 一卡2卡三卡四卡精品乱码亚洲| 国产成人午夜福利电影在线观看| 我的老师免费观看完整版| 日韩制服骚丝袜av| 校园人妻丝袜中文字幕| 99在线视频只有这里精品首页| 一卡2卡三卡四卡精品乱码亚洲| 欧美日韩国产亚洲二区| 久久久亚洲精品成人影院| 黑人高潮一二区| 一区二区三区免费毛片| 国产一级毛片在线| 日韩成人伦理影院| 又粗又爽又猛毛片免费看| 国产成人免费观看mmmm| 欧美zozozo另类| 麻豆乱淫一区二区| 国产日韩欧美在线精品| av天堂中文字幕网| av免费观看日本| 噜噜噜噜噜久久久久久91| 日日干狠狠操夜夜爽| 国产人妻一区二区三区在| 老女人水多毛片| 女人久久www免费人成看片 | 午夜激情欧美在线| 亚洲在线自拍视频| av天堂中文字幕网| 看非洲黑人一级黄片| 嫩草影院精品99| 99在线视频只有这里精品首页| 别揉我奶头 嗯啊视频| 国产激情偷乱视频一区二区| 日本免费在线观看一区| 搡女人真爽免费视频火全软件| 日本免费a在线| 一级黄片播放器| 国产伦精品一区二区三区四那| 国产av码专区亚洲av| 日本与韩国留学比较| 国产亚洲5aaaaa淫片| 日韩欧美三级三区| 日韩三级伦理在线观看| 亚洲欧美日韩卡通动漫| 午夜激情福利司机影院| 中文字幕久久专区| 精品熟女少妇av免费看| 免费av毛片视频| 国产精品久久视频播放| 九九爱精品视频在线观看| 看免费成人av毛片| 亚洲精品日韩av片在线观看| 超碰av人人做人人爽久久| 亚洲国产色片| 亚洲av成人精品一区久久| 在线天堂最新版资源| 看非洲黑人一级黄片| 国产精品一区二区三区四区久久| 丝袜美腿在线中文| 午夜福利成人在线免费观看| www.av在线官网国产| 五月伊人婷婷丁香| 午夜亚洲福利在线播放| 建设人人有责人人尽责人人享有的 | 久热久热在线精品观看| 欧美激情国产日韩精品一区| 只有这里有精品99| 嫩草影院精品99| 久久精品国产亚洲网站| kizo精华| av在线观看视频网站免费| 我要看日韩黄色一级片| 三级国产精品片| 久久亚洲国产成人精品v| 亚洲精品乱码久久久久久按摩| 男女边吃奶边做爰视频| 在线a可以看的网站| 亚洲激情五月婷婷啪啪| 国产精品日韩av在线免费观看| 亚洲人与动物交配视频| 欧美高清成人免费视频www| 亚洲精品亚洲一区二区| 97在线视频观看| 中文字幕精品亚洲无线码一区| 国产激情偷乱视频一区二区| 91午夜精品亚洲一区二区三区| 欧美性猛交黑人性爽| 亚洲人与动物交配视频| 免费av观看视频| 久久久午夜欧美精品| 我要看日韩黄色一级片| 桃色一区二区三区在线观看| 熟妇人妻久久中文字幕3abv| 哪个播放器可以免费观看大片| 日本爱情动作片www.在线观看| 嘟嘟电影网在线观看| 人妻制服诱惑在线中文字幕| 五月伊人婷婷丁香| 少妇猛男粗大的猛烈进出视频 | av视频在线观看入口| 女的被弄到高潮叫床怎么办| 黄色一级大片看看| 精品久久久久久电影网 | 日本黄色视频三级网站网址| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费看光身美女| 26uuu在线亚洲综合色| 午夜福利在线观看吧| 中文字幕人妻熟人妻熟丝袜美| 国产熟女欧美一区二区| 99久久中文字幕三级久久日本| 国产精华一区二区三区| 啦啦啦观看免费观看视频高清| 欧美一区二区精品小视频在线| a级毛片免费高清观看在线播放| 精品无人区乱码1区二区| 欧美潮喷喷水| 午夜日本视频在线| 乱码一卡2卡4卡精品| 两性午夜刺激爽爽歪歪视频在线观看| 免费看a级黄色片| 岛国在线免费视频观看| 高清毛片免费看| 插逼视频在线观看| 国产精品人妻久久久影院| 精品一区二区三区视频在线| 亚洲人成网站在线观看播放| 免费看a级黄色片| 亚洲av熟女| 国产精品永久免费网站| 成年女人看的毛片在线观看| 高清av免费在线| 少妇熟女欧美另类| 两性午夜刺激爽爽歪歪视频在线观看| 26uuu在线亚洲综合色| 日韩高清综合在线| 亚洲三级黄色毛片| 国产乱人视频| 美女xxoo啪啪120秒动态图| 中文字幕亚洲精品专区| 国产精品一区二区在线观看99 | 久久精品国产自在天天线| 亚洲高清免费不卡视频| 日韩视频在线欧美| 久久久午夜欧美精品| 久久99热这里只频精品6学生 | 成人毛片a级毛片在线播放| 欧美性猛交黑人性爽| 超碰97精品在线观看| 一个人看的www免费观看视频| 最近最新中文字幕免费大全7| 在现免费观看毛片| 麻豆一二三区av精品| 国产午夜精品久久久久久一区二区三区| 国产黄片视频在线免费观看| 亚洲国产精品国产精品| 成人毛片60女人毛片免费| 成人漫画全彩无遮挡| 中文字幕免费在线视频6| 精品人妻一区二区三区麻豆| 亚洲在久久综合| 色综合亚洲欧美另类图片| 国内精品宾馆在线| 日韩欧美在线乱码| 欧美日韩精品成人综合77777| 美女内射精品一级片tv| 成人美女网站在线观看视频| 蜜臀久久99精品久久宅男| av视频在线观看入口| 老女人水多毛片| 亚洲激情五月婷婷啪啪| 嫩草影院新地址| 高清在线视频一区二区三区 | 一级毛片久久久久久久久女| 成人亚洲精品av一区二区| 插逼视频在线观看| 久久久久久久午夜电影| 国产精品久久久久久久久免| 欧美区成人在线视频| 97超视频在线观看视频| 毛片一级片免费看久久久久| 欧美另类亚洲清纯唯美| 国产毛片a区久久久久| 亚洲av免费在线观看| 久久人人爽人人片av| 精品酒店卫生间| 国产精品一区二区三区四区久久| 三级国产精品片| 白带黄色成豆腐渣| 五月玫瑰六月丁香| 国产乱来视频区| 成年版毛片免费区| 国产亚洲91精品色在线| 免费在线观看成人毛片| 色综合色国产| 丰满人妻一区二区三区视频av| 最近2019中文字幕mv第一页| 人妻制服诱惑在线中文字幕| 日韩三级伦理在线观看| 久久人人爽人人片av| 国语自产精品视频在线第100页| 五月玫瑰六月丁香| 成人鲁丝片一二三区免费| 九九在线视频观看精品| 麻豆国产97在线/欧美| 久久精品国产99精品国产亚洲性色| 99热全是精品| www.av在线官网国产| 国产欧美另类精品又又久久亚洲欧美| 三级国产精品片| 国产老妇女一区| 综合色av麻豆| 极品教师在线视频| 天堂中文最新版在线下载 | 亚洲无线观看免费| 国产私拍福利视频在线观看| 国产真实乱freesex| 搞女人的毛片| 久久久久久久久久黄片| 99热全是精品| 亚洲成人av在线免费| 国产精品日韩av在线免费观看| 成人亚洲精品av一区二区| 两个人的视频大全免费| av在线天堂中文字幕| 乱码一卡2卡4卡精品| 国产极品天堂在线| 天天躁日日操中文字幕| 亚洲丝袜综合中文字幕| 在线观看一区二区三区| 国产精品嫩草影院av在线观看| 亚洲,欧美,日韩| 免费不卡的大黄色大毛片视频在线观看 | 岛国在线免费视频观看| 联通29元200g的流量卡| 日本午夜av视频| 国产精品久久视频播放| 国产又色又爽无遮挡免| 亚洲怡红院男人天堂| 在线免费观看不下载黄p国产| 在线观看66精品国产| 99在线人妻在线中文字幕| 极品教师在线视频| 欧美精品一区二区大全| 少妇猛男粗大的猛烈进出视频 | 熟女人妻精品中文字幕| 女人久久www免费人成看片 | 国产成人午夜福利电影在线观看| 日本wwww免费看| 国产麻豆成人av免费视频| 波野结衣二区三区在线| 我要看日韩黄色一级片| 日本一本二区三区精品| 美女大奶头视频| 精品欧美国产一区二区三| 伊人久久精品亚洲午夜| 国产亚洲91精品色在线| 免费黄网站久久成人精品| 日韩av在线免费看完整版不卡| 午夜福利在线观看吧| 中国美白少妇内射xxxbb| 日韩制服骚丝袜av| 国产精品国产三级国产专区5o | 国产精品熟女久久久久浪| 不卡视频在线观看欧美| 国产成人福利小说| 别揉我奶头 嗯啊视频| 国产亚洲5aaaaa淫片| 国产精华一区二区三区| 看免费成人av毛片| 91久久精品国产一区二区三区| 亚洲天堂国产精品一区在线| 亚洲国产精品sss在线观看| 国产精华一区二区三区| 国内精品一区二区在线观看| 久久久久久久午夜电影| 精品人妻一区二区三区麻豆| 少妇人妻一区二区三区视频| 亚洲国产精品成人久久小说| 亚洲av二区三区四区| 最近视频中文字幕2019在线8| 久久久a久久爽久久v久久| 黄片wwwwww| 中文字幕精品亚洲无线码一区| 高清午夜精品一区二区三区| 日日干狠狠操夜夜爽| 成人午夜精彩视频在线观看| 高清日韩中文字幕在线| 嘟嘟电影网在线观看| 日韩人妻高清精品专区| 久久久久九九精品影院| 不卡视频在线观看欧美| 一级毛片电影观看 | 精品人妻偷拍中文字幕| 能在线免费看毛片的网站| 国产精品99久久久久久久久| 99在线人妻在线中文字幕| av线在线观看网站| 91aial.com中文字幕在线观看| 能在线免费观看的黄片| 熟女电影av网| 日本爱情动作片www.在线观看| 赤兔流量卡办理| 黑人高潮一二区| 精品99又大又爽又粗少妇毛片| 男女啪啪激烈高潮av片| 国产亚洲av嫩草精品影院| 成年版毛片免费区| 青青草视频在线视频观看| 中文字幕久久专区| 亚洲av.av天堂| 久久久久久久国产电影| 2021天堂中文幕一二区在线观| 好男人在线观看高清免费视频| 老司机影院毛片| 午夜亚洲福利在线播放| 久久久欧美国产精品| 国产国拍精品亚洲av在线观看| 成年av动漫网址| 黄片无遮挡物在线观看| 成人性生交大片免费视频hd| 日本黄色视频三级网站网址| 69av精品久久久久久| 亚洲无线观看免费| 亚洲欧洲日产国产| 99热这里只有是精品50| 婷婷六月久久综合丁香| 久热久热在线精品观看| 日产精品乱码卡一卡2卡三| 99久久精品国产国产毛片| 三级男女做爰猛烈吃奶摸视频| 最近中文字幕高清免费大全6| 在线播放无遮挡| 亚洲国产欧美人成| 国产色爽女视频免费观看| www日本黄色视频网| 午夜精品在线福利| 国产精品电影一区二区三区| 深夜a级毛片| 久久人人爽人人片av| 欧美一区二区精品小视频在线| 18禁裸乳无遮挡免费网站照片| 国产精品国产高清国产av| 深爱激情五月婷婷| 少妇的逼水好多| 国产成人福利小说| 男女那种视频在线观看| 亚洲伊人久久精品综合 | 人体艺术视频欧美日本| 免费搜索国产男女视频| 精品少妇黑人巨大在线播放 | 国产午夜精品论理片| 丰满乱子伦码专区| 免费人成在线观看视频色| 国产精品一及| 蜜桃久久精品国产亚洲av| 91av网一区二区| 中文字幕精品亚洲无线码一区| 国产精品综合久久久久久久免费| 午夜福利视频1000在线观看| 人人妻人人澡人人爽人人夜夜 | 亚洲av电影在线观看一区二区三区 | 亚洲av福利一区| 国产精品久久久久久久久免| 我的老师免费观看完整版| 免费观看的影片在线观看| 人体艺术视频欧美日本| 免费搜索国产男女视频| 国产成人freesex在线| 国产精品精品国产色婷婷| 小蜜桃在线观看免费完整版高清| 国产成人福利小说| 中文字幕久久专区| 高清日韩中文字幕在线| 别揉我奶头 嗯啊视频| 26uuu在线亚洲综合色| 免费看日本二区| 久久99精品国语久久久| 午夜激情福利司机影院| 日韩一本色道免费dvd| 国产精品一区www在线观看| 精品熟女少妇av免费看| 欧美成人一区二区免费高清观看| 网址你懂的国产日韩在线| 成人午夜高清在线视频| 日本熟妇午夜| 日本爱情动作片www.在线观看| 伦精品一区二区三区| av在线天堂中文字幕| 中文资源天堂在线| 精品久久国产蜜桃| 精品久久久久久久人妻蜜臀av| 日韩欧美精品v在线| 成年免费大片在线观看| 欧美又色又爽又黄视频| 免费黄色在线免费观看| 中文天堂在线官网| 精品一区二区三区视频在线| 两个人的视频大全免费| 国产片特级美女逼逼视频| 午夜福利网站1000一区二区三区| 亚洲av免费高清在线观看| 国内精品美女久久久久久| 日本欧美国产在线视频| 日本黄大片高清| 亚洲四区av| 精品一区二区三区人妻视频| 国产老妇伦熟女老妇高清| 1000部很黄的大片| 成人性生交大片免费视频hd| 国产女主播在线喷水免费视频网站 | 欧美一区二区亚洲| 欧美日韩国产亚洲二区| 日本免费a在线| 亚洲欧美日韩东京热| 岛国在线免费视频观看| 国产精品一区二区三区四区免费观看| 国产高清有码在线观看视频| 欧美激情久久久久久爽电影| 91久久精品国产一区二区成人| 国产精品国产高清国产av| 日日摸夜夜添夜夜爱| 国产毛片a区久久久久| 99久久无色码亚洲精品果冻| 国语自产精品视频在线第100页| 六月丁香七月| 人体艺术视频欧美日本| 男人的好看免费观看在线视频| 成年女人看的毛片在线观看| 春色校园在线视频观看| 精品国产露脸久久av麻豆 | 欧美一区二区亚洲| 91精品一卡2卡3卡4卡| a级毛片免费高清观看在线播放| 国产午夜精品一二区理论片| 变态另类丝袜制服| 日韩一区二区视频免费看| 久久精品综合一区二区三区| 国产成人午夜福利电影在线观看| 中国国产av一级| 欧美激情国产日韩精品一区| 91在线精品国自产拍蜜月| 黄色一级大片看看| 天天躁日日操中文字幕| 美女脱内裤让男人舔精品视频| 2021少妇久久久久久久久久久| 1000部很黄的大片| 国产麻豆成人av免费视频| 日日啪夜夜撸| 啦啦啦啦在线视频资源| 精华霜和精华液先用哪个| 晚上一个人看的免费电影| 午夜a级毛片| 国产av一区在线观看免费| 美女黄网站色视频| 成人美女网站在线观看视频| 麻豆国产97在线/欧美| 国产精品人妻久久久久久| 国产精品女同一区二区软件| 日韩,欧美,国产一区二区三区 | 精品国产一区二区三区久久久樱花 | 欧美zozozo另类| av天堂中文字幕网| 免费观看的影片在线观看| 欧美精品一区二区大全| 久久人人爽人人片av| 深夜a级毛片| 九九爱精品视频在线观看| 国产精品一及| 青春草视频在线免费观看| 淫秽高清视频在线观看| 在线播放无遮挡| 有码 亚洲区| 中文在线观看免费www的网站| 国产精品久久久久久精品电影| 亚洲国产精品国产精品| 国产精品人妻久久久影院| 免费一级毛片在线播放高清视频| 亚洲最大成人av| 国产午夜精品论理片| 亚洲欧美日韩卡通动漫| 最近中文字幕2019免费版| 午夜精品在线福利| 日韩av在线免费看完整版不卡| 最近的中文字幕免费完整| 久久99蜜桃精品久久| 成年av动漫网址| 亚洲av免费在线观看| 99久久精品一区二区三区| 久久久精品欧美日韩精品| 两性午夜刺激爽爽歪歪视频在线观看| 男女视频在线观看网站免费| 最近中文字幕2019免费版| 国产片特级美女逼逼视频| 日韩一本色道免费dvd| 又爽又黄无遮挡网站| 亚洲美女视频黄频| 一级爰片在线观看| 欧美性猛交黑人性爽| 亚洲,欧美,日韩| 亚洲五月天丁香| 午夜福利在线观看免费完整高清在| 亚洲国产精品成人久久小说| 日本五十路高清| 国产精品伦人一区二区| 久久99热这里只有精品18| 国产精品麻豆人妻色哟哟久久 | 日日摸夜夜添夜夜爱| 高清午夜精品一区二区三区| 日韩在线高清观看一区二区三区| 狠狠狠狠99中文字幕| 国产精品一区二区在线观看99 | 天天躁日日操中文字幕| 国产午夜福利久久久久久| 国产淫片久久久久久久久| 免费看a级黄色片| 国产精品久久电影中文字幕| 国产精品一及| 国产在线男女| 91在线精品国自产拍蜜月| 精品国产露脸久久av麻豆 | 最近手机中文字幕大全| 边亲边吃奶的免费视频| 一级黄片播放器| 最新中文字幕久久久久| 国产久久久一区二区三区| 黑人高潮一二区| 91精品一卡2卡3卡4卡| 色播亚洲综合网| 天堂中文最新版在线下载 | 日韩一本色道免费dvd| 国产单亲对白刺激| 亚洲av男天堂| 可以在线观看毛片的网站| 精品免费久久久久久久清纯| 成人亚洲欧美一区二区av| 久99久视频精品免费| 少妇人妻一区二区三区视频| 亚洲精品一区蜜桃| 久久久久久国产a免费观看| 国产伦在线观看视频一区| 97热精品久久久久久| 大话2 男鬼变身卡| 麻豆精品久久久久久蜜桃| 亚洲美女视频黄频| 啦啦啦观看免费观看视频高清| 深爱激情五月婷婷| 免费av不卡在线播放| 国产人妻一区二区三区在| 中文欧美无线码| 国产大屁股一区二区在线视频| 国产精品日韩av在线免费观看| 日韩国内少妇激情av| 性色avwww在线观看| 神马国产精品三级电影在线观看| 亚洲av不卡在线观看| 亚洲精品色激情综合| 日日摸夜夜添夜夜添av毛片| 99久久无色码亚洲精品果冻| 1000部很黄的大片| 超碰97精品在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 日韩人妻高清精品专区| 欧美成人精品欧美一级黄| 亚洲自偷自拍三级| 国产免费男女视频| 2021少妇久久久久久久久久久| 免费看a级黄色片| 欧美成人午夜免费资源| 色哟哟·www|