胥元剛 郝帆帆

西安石油大學石油工程學院

超高強度抽油桿P-S-N曲線擬合新方法

胥元剛 郝帆帆

西安石油大學石油工程學院

P-S-N曲線是超高強度抽油桿力學性能分析的主要依據(jù)，其擬合方法通常參照于普通抽油桿。然而實驗數(shù)據(jù)表明，超高強度抽油桿并不滿足普通抽油桿P-S-N曲線擬合模型的假設條件。為了研究適應于超高強度抽油桿的P-S-N曲線擬合方法，提出了基于Box-Cox變量變換法的三參數(shù)P-S-N曲線擬合模型，其中Box-Cox變量變換法解決了一定循環(huán)應力下，疲勞壽命不滿足正態(tài)分布的問題；三參數(shù)P-S-N曲線模型更符合典型S-N曲線的趨勢。實驗計算結果表明，新方法下的P-S-N曲線擬合度更高，相對誤差僅為1.46%～4.12%，僅為以往方法的11%～30%，新方法擬合的曲線更能體現(xiàn)超高強度抽油桿力學性能，精確度更高。

P-S-N曲線；超高強度；抽油桿；Box-Cox變量變換法

P-S-N曲線是指在給定存活概率下，疲勞實驗所得循環(huán)應力與疲勞壽命的關系曲線，簡稱概率疲勞壽命曲線［1］。抽油桿高載荷往復運動過程中受到交變應力作用，達到一定限度時抽油桿斷裂，即發(fā)生疲勞破壞。為研究抽油桿疲勞性能，進行室內疲勞實驗研究，通過擬合數(shù)據(jù)獲得的P-S-N曲線對曲線擬合程度要求很高。

針對抽油桿P-S-N曲線的擬合方法，多位學者進行了相關研究。宋開利（2003）提出利用K-S方法驗證疲勞實驗的數(shù)據(jù)分布，以及P-S-N曲線線性擬合模型［2］；林元華等（2005）基于抽油桿軸向載荷的數(shù)學模型，建立了疲勞壽命的預測模型［3］；李大建等（2015）基于疲勞實驗數(shù)據(jù)計算了H級抽油桿P-S-N曲線的方程［4］。以上研究均作出假設：同一循環(huán)應力下，疲勞壽命滿足正態(tài)分布；指定存活概率下疲勞壽命和循環(huán)應力對數(shù)線性相關。然而，大量疲勞實驗數(shù)據(jù)表明：不同循環(huán)應力下，疲勞壽命的分布不完全滿足正態(tài)分布，且方差不同；中、長壽命區(qū)材料在指定存活概率下，疲勞壽命和循環(huán)應力的對數(shù)線性相關程度低，從而導致回歸誤差大［5-8］。超高強度抽油桿的主要特點為疲勞強度高，屬于中、長壽命區(qū)材料，其擬合模型不再滿足擬合普通抽油桿的假設條件，需要提出新的擬合方法。

利用Box-Cox變量變換法解決疲勞壽命不服從正態(tài)分布的問題，能夠實現(xiàn)高擬合優(yōu)度，又能保證信息不丟失。而針對模型的不適用性，根據(jù)典型S-N曲線趨勢提出三參數(shù)擬合模型。

1 Box-Cox變量變換法

Box-Cox variable transformation method

Box和Cox于1964年提出Box-Cox變量變換法，原理為對因變量y作以下變換［9］

式中，λ為待定參數(shù)。

Box-Cox變量變換法的特點：變換因變量y前后，自身遵循規(guī)律不發(fā)生變化，信息不丟失，且回歸模型滿足正態(tài)分布。對于X=［x1，x2，…，xn］T和Y=［y1，y2，…，yn］T，變換后Y（λ）=［y1（λ），y2（λ），…，yn（λ）］T滿足方程

式中，α為矩陣Y（λ）和X的相關系數(shù)。

式中，zi（λ）為yi經過λ值轉換的求參因變量。

確定λ值的方法：先確定λ值的取值范圍，給定初值，通過計算得出Z（λ）= ［z1（λ），z2（λ）…zn（λ）］T，誤差殘差平方和R（λ，Z（λ））=Z（λ）T［I–X（XTX）–1XT］Z（λ），λ所求值即為滿足minR（λ，Z（λ））的值。

疲勞實驗數(shù)據(jù)為某一循環(huán)應力對應的多組疲勞壽命數(shù)值，即X樣本為定值，認為minR（λ，Z（λ））=Z（λ）TZ（λ）。據(jù)此求得不同循環(huán)應力下對應的λ值，完成所有疲勞實驗數(shù)據(jù)的Box-Cox變量變換。解決了不同循環(huán)應力下，疲勞壽命的分布不完全滿足正態(tài)分布的問題。

2 三參數(shù)P-S-N曲線

Three-parameterP-S-Ncurve

以往的抽油桿P-S-N曲線模型原理：指定存活概率P下，疲勞壽命N與循環(huán)應力S的對數(shù)線性相關，為

式中，NP為存活概率（即可靠度）P下的疲勞壽命；S為循環(huán)應力，MPa；mP，CP為存活概率P下的待定系數(shù)。

令YP=lgNP，X=lgS，β1=mP，β2=lgCP。式（4）變換為線性關系式為

當疲勞壽命N趨近無限大時，循環(huán)應力趨近疲勞極限S0，S-N曲線趨勢與指數(shù)函數(shù)相近。提出適用于中、長壽命區(qū)材料S-N曲線的模型為

式中，m，C為材料常數(shù)。

給定存活率P下，變形為三參數(shù)P-S-N曲線模型方程為

式中，SP0為存活率P下的疲勞極限。令X=lg（S–SP0）、YP=lgNP、β1=mP、β2=lgCP，式（7）可以變形為式（5）的形式，與以往抽油桿擬合過程不同的是X值的計算方法以及YP經過了Box-Cox變量變換這個步驟。因此雖然最后擬合模式均為線性擬合，但參數(shù)發(fā)生變化。

3 最小二乘法線性擬合

Alignment fitting of least square method

對線性方程式（5）的擬合使用最小二乘法

據(jù)此獲得k組循環(huán)應力X=［lg（S1–SP0），lg（S2–SP0），…，lg（Sk–SP0）］T下的Yp=［yP1，yP2，…，yPk］T。最小二乘法擬合線性關系的原理為

擬合產生的誤差V表示為

最后得出擬合公式為

4 P-S-N曲線擬合思路

P-S-Ncurve fitting principle

根據(jù)上述理論，確定P-S-N曲線擬合步驟。

（1）進行疲勞實驗，測定循環(huán)應力Si（i=1，2，…，k）的疲勞壽命Nj（j=1，2，…，n）。

（2）Box-Cox變量變換，循環(huán)應力Si下，先給定λ值的取值范圍，給定初始值，利用式（3）由Nj求得Z（λ），在取值范圍內試算滿足minR（λ，Z（λ））=Z（λ）TZ（λ）的λ值。以此方法計算出其他循環(huán)應力下的λ值。用式（1）對Si下Nj進行Box-Cox變量變換，使其滿足正態(tài)分布，得出Si對應新的Nj值。

（3）給定可靠度P，根據(jù)式（8）計算YP=［yP1，yP2，…，yPk］T。

（4）假定SP0，計算X=［lg （S1–SP0），lg（S2–SP0），…，lg（Sk–SP0）］T。

（5）根據(jù)式（9）和式（10）最小二乘法擬合線性方程計算β1，β2。

（6）根據(jù)式（11）計算V值，然后返回第（4）步重新假設SP0值，V達到最小值則停止計算，記錄試算最后一步的SP0值。

（7）得出可靠度P下，超高強度抽油桿P-S-N曲線方程式（12）。

5 實驗數(shù)據(jù)

Experimental data calculation

根據(jù)SY/T 5029—2013《抽油桿》附錄L中鋼制抽油桿疲勞性能要求和實驗方法，利用PLG-300C電磁諧振高頻疲勞實驗機進行了18組超高強度抽油桿的疲勞實驗，每組7根試件（2根備用）。試件發(fā)生疲勞斷裂時，如果疲勞壽命不低于106，則為合格，認為是有效數(shù)據(jù)，否則無效。篩選出2種類型的超高強度抽油桿（HY型和HL型）在不同循環(huán)應力下的疲勞壽命數(shù)據(jù)，見表1。分別利用新方法和以往方法對實驗所得數(shù)據(jù)進行了P-S-N曲線擬合，并計算相應的誤差值，計算結果見表2。

表1 超高強度抽油桿疲勞實驗Table 1 Fatigue test on ultra-high-strength pumping rod

表2 計算結果Table 2 Calculation results

由表2可以看出，可靠度為50%和95%時，采用新方法擬合的曲線偏差更小，其相對誤差降低到以往方法的11%～30%，說明基于Box-Cox變量變換法的三參數(shù)P-S-N曲線模型與實驗數(shù)據(jù)吻合程度更高，更能代表超高強度抽油桿的力學性能。

6 結論

Conclusions

（1）針對一定循環(huán)應力下的超高強度抽油桿，其疲勞壽命不完全滿足正態(tài)分布以及以往模型不適用的問題，提出了基于Box-Cox變量變換法的三參數(shù)P-S-N曲線模型：Box-Cox變量變換法能夠保證信息不丟失，且模型滿足正態(tài)分布；三參數(shù)P-S-N曲線模型符合典型S-N曲線的趨勢，是對以往模型的優(yōu)化。

（2）實驗結果表明，新方法將擬合曲線的相對誤差降低到了1.46%～4.12%，僅為以往模型的11%～30%，提高了P-S-N曲線的擬合精度和可信度。因此，新方法更適用于超高強度抽油桿P-S-N曲線的擬合，更能代表超高強度抽油桿的力學性能。。

References:

［1］謝里陽，劉建中.樣本信息聚集原理與P-S-N曲線擬合方法［J］.機械工程學報，2013，49（15）：96-103.XIE Liyang,LIU Jianzhong.Principle of sample polymerization and method ofP-S-Ncurve fitting ［J］.Journal of Mechanical Engineering,2013,49（15）: 96-103.

［2］宋開利.超高強度抽油桿P-S-N曲線實驗研究［J］.石油礦場機械，2003，32（5）：27-29.SONG Kaili.Experimental research ofP-S-Ncurves for sucker rod with super high strength［J］.Oil Field Equipment,2003,32（5）: 27-29.

［3］林元華，張德平，駱發(fā)前，施太和.抽油桿疲勞壽命計算研究［J］.石油鉆采工藝，2005，27（6）：66-69.LIN Yuanhua,ZHANG Deping,LU Faqian,SHI Taihe.Research on fatigue life calculation method of sucker rod［J］.Oil Drilling & Production Technology,2005,27（6）:66-69.

［4］李大建，陸梅，郭靖，張建偉，劉濤，梁毅.H級抽油桿疲勞性能實驗分析［J］.石油礦場機械，2015，44（12）：42-44.LI Dajian,LU Mei,GUO Jing,ZHANG Jianwei,LIU Tao,LIANG Yi.Fatigue test analysis of H sucker rod ［J］.Oil Field Equipment,2015,44（12）: 42-44.

［5］李洪雙，呂震宙.估計疲勞壽命三參數(shù)P-S-N曲線的新方法［J］.機械強度，2007，29（2）：300-304.LI Hongshuang,LYU Zhengzhou.New method for estimating fatigue life three parameterP-S-Ncurves ［J］.Journal of Mechanical Strength,2007,29（2）: 300-304.

［6］陳家權，陳國軍，陳非凡，張云泉，溫杰明.Box-Cox 變量變換擬合法估計三參數(shù)P-S-N曲線研究［J］.機械強度，2011，33（3）：358-362.CHEN Jiaquan,CHEN Guojun,CHEN Feifan,ZHANG Yunquan,WEN Jieming.Reserch on Box-Cox variable transform fitting method for estimating fatigue life three parameterP-S-Ncurves ［J］.Journal of Mechanical Strength,2011,33（3）: 358-362.

［7］傅惠民，高鎮(zhèn)同P-S-N曲線擬合法［J］.航空學報，1988，9（7）：338-341.FU Huimin,GAO Zhentong.Amethod for fittingP-S-Ncurves ［J］.Astronautica Sinica,1988,9（7）: 338-341.

［8］謝金標，姚衛(wèi)星.疲勞S-N曲線擬合的雙加權最小二乘法［J］.宇航學報，2010，31（6）：1661-1665.XIE Jinbiao,YAO Weixing.Double weighted least square method for fatigueS-Ncurve fitting ［J］.Journal of Astronautic,2010,31（6）: 1661-1665.

［9］左文鼎，夏祥謙，王力賓.Box-Cox 變換在MATLAB上的實現(xiàn)［J］.云南財經大學學報，2007，23（6）：116-119.ZUO Wending,XIA Xiangqian,WANG Libin.Box-Cox variable transform and its implementation by MATLAB［J］.Journal of Astronautic,2007,23（6）: 116-119.

（修改稿收到日期 2017-03-12）

〔編輯 李春燕〕

A new fitting method forP-S-Ncurve of ultra-high-strength sucker rod

XU Yuangang,HAO Fanfan
College of Petroleum Engineering,Xi’an Shiyou University,Xi’an710065,Shaanxi,China

P-S-Ncurve is one important basis for analyzing the mechanical performance of ultra-high-strength sucker rod,and its fitting method is generally referred to that of conventional sucker rods.Based on experimental data,however,ultra-high-strength sucker rod can’t satisfy the assumptions ofP-S-Ncurve fitting model which is used for conventional sucker rods.A three-parameterP-S-Ncurve fitting model based on Box-Cox variable transformation method was developed to search for a fitting method suitable for theP-S-Ncurve of ultra-high-strength sucker rod.By means of Box-Cox variable transformation method,fatigue life can satisfy normal distribution in a certain cyclic stress.The three-parameterP-S-Ncurve model is more accordant with the trend of typical S-N curve.It is indicated from experimental calculation results that by virtue of this new method,the fitting degree ofP-S-Ncurve is higher with relative error only 1.46%-4.12%,which is just 11%-30% of the conventional value.The fitting based on this new method can embody the mechanical performance of ultra-high-strength sucker rod more accurately.

P-S-Ncurve; ultra high strength; sucker rod; Box-Cox variable transformation method

胥元剛，郝帆帆.超高強度抽油桿P-S-N曲線擬合新方法［J］.石油鉆采工藝，2017，39（3）：334-337.

TE928

：A

1000–7393（2017 ）03–0334–04DOI:10.13639/j.odpt.2017.03.015

: XU Yuangang,HAO Fanfan.A new fitting method forP-S-Ncurve of ultra-high-strength sucker rod［J］.Oil Drilling& Production Technology,2017,39(3): 334-337.

西安石油大學研究生創(chuàng)新基金項目（編號：2015cx140104）；西安石油大學優(yōu)秀碩士學位論文培育項目（編號：2015yp1401 04）。

胥元剛（1963-），1998年畢業(yè)于中國石油大學（華東）油氣田開發(fā)工程專業(yè)，博士研究生，現(xiàn)從事采油工程理論與技術方向的研究和教學，教授。通訊地址：（710065）陜西省西安市雁塔區(qū)電子二路18號西安石油大學。E-mail：ygxu@xsyu.edu.cn

郝帆帆（1990-），西安石油大學在讀碩士研究生現(xiàn)，主要從事采油工程理論與技術方面的研究工作。通訊地址：（710065）陜西省西安市雁塔區(qū)電子二路18號西安石油大學。E-mail：fan_fanhao@126.com