林坤,馬朝選,董強,黃亞,彭曉光

(中國船舶重工集團公司第七一八研究所,河北 邯鄲 056027)

0 引 言

石油測井儀器的工作環(huán)境極其惡劣。隨著測井深度的增加,溫度不斷升高,壓力也不斷加大。在井下深度至7 000 m時,溫度會達(dá)到175 ℃,壓力也會達(dá)到140 MPa,儀器外殼體承受的壓力不斷增大。

隨著鉆井的井眼直徑越來越小,測井儀器的外徑越來越小,在滿足外部壓力的條件下,外殼體的壁厚越小,內(nèi)徑就越大,外殼體內(nèi)部就有較大的空間布置電路和探測器。

在文獻[1]和[2]中,應(yīng)用傳統(tǒng)的設(shè)計方法以及力學(xué)的第三和第四強度理論進行測井儀外殼的強度設(shè)計與校核,為計算測井儀外殼體壁厚及強度校核提供了理論依據(jù),但安全系數(shù)的選擇如果不合適的話,都會造成對儀器外殼強度的影響。同時,在測井儀外殼設(shè)計好后,生產(chǎn)加工好后進行壓力試驗,如果不滿足要求,再重新設(shè)計,這樣測井儀外殼的設(shè)計生產(chǎn)周期大大延長,導(dǎo)致時間和成本上升。

本文應(yīng)用有限元分析ANSYS軟件的結(jié)構(gòu)分析模塊與可靠性分析模塊PDS(Probability Design System)相結(jié)合,采用Monte-Carlo蒙特卡羅法,以四探頭高精度自然伽馬測井儀殼體壁厚、壓力載荷及材料的彈性模量等作為隨機輸入變量,對外殼體的蒙特卡羅模擬并進行分析,得出不確定設(shè)計因素與目標(biāo)變量的靈敏度關(guān)系,影響強度的主要影響因素和殼體的可靠度,為結(jié)構(gòu)可靠性優(yōu)化設(shè)計提供有用數(shù)據(jù),減少試驗時間,提高設(shè)計效率。

1 外殼體應(yīng)力-強度模型

機械結(jié)構(gòu)的破壞概率可以表示為

Pf=P{g(X)<0}

(1)

式中,X={x1,x2,…,xn}為具有n維隨機變量的向量;g(X)為結(jié)構(gòu)的極限狀態(tài)方程。

機械結(jié)構(gòu)可靠性研究的主要對象是結(jié)構(gòu)的強度R與應(yīng)力S,通過研究兩者之間的功能函數(shù)確定機械結(jié)構(gòu)的可靠性。由此得到的結(jié)構(gòu)功能函數(shù)可以表述為應(yīng)力-強度模型

Z=R-S

(3)

概率極限狀態(tài)結(jié)構(gòu)設(shè)計中必須滿足條件(4)

Z=g(R,S)=R-S≥0

(4)

根據(jù)應(yīng)力強度干涉理論,在四探頭高精度自然伽馬測井儀外殼結(jié)構(gòu)的使用過程中,不允許出現(xiàn)應(yīng)力超過屈服強度的事件發(fā)生,所以外殼結(jié)構(gòu)失效的準(zhǔn)則為

σmax≥σs

(5)

式中,σmax為外殼結(jié)構(gòu)在工作中出現(xiàn)的最大應(yīng)力;σs為外殼材料的屈服強度。于是應(yīng)力-強度模型為

Z=g(X)=σs-σmax≥0

(6)

蒙特卡羅法又稱為隨機模擬法或統(tǒng)計試驗法[3],將隨機變量xi代入功能函數(shù)g(x),得出一個函數(shù)值g(xi)。若g(xi)<0,則在計算機程序中記入一次失效的實現(xiàn);若g(x)>0,則不記入。這樣就完成了一次計算,再產(chǎn)生下一個隨機數(shù),重復(fù)上面的計算,直至完成預(yù)定的實驗次數(shù)為止。用蒙特卡羅方法表示的失效概率可以表示為

(7)

式中,N為抽樣模擬總數(shù);當(dāng)g(xi)<0,I[g(xi)]=1,反之,I[g(xi)]=0。

2 基于ANSYS/PDS外殼體強度蒙特卡羅模擬

四探頭高精度自然伽馬測井儀外殼體初步設(shè)計參數(shù):承受外部壓強140 MPa;外殼體長度3 000 mm;殼體外徑45 mm;殼體內(nèi)徑35 mm。外殼體的材料采用一種超高強度鈦合金TC11,屈服點為σs=1 030 MPa,并采用一定的熱處理工藝,達(dá)到高強度。

應(yīng)用ANSYS軟件中PDS(Probability Design System)概率設(shè)計模塊對該測井儀外殼體結(jié)構(gòu)進行蒙特卡羅模擬與可靠性分析。

2.1 生成分析文件

該階段包括對輸入?yún)?shù)進行參數(shù)化設(shè)定,模型的建立、網(wǎng)格劃分,加載求解和后處理[4-6]。

(1) 定義輸入?yún)?shù)。輸入?yún)?shù)主要是指影響結(jié)構(gòu)行為的不確定因素如材質(zhì)的不確定性,結(jié)構(gòu)幾何尺寸的不確定性等,輸入的相關(guān)參數(shù)見表1。

表1 輸入?yún)?shù)表

(2) 建立模型。本文采用圖形用戶接口GUI方式,用自下而上的方法建立外殼體模型,通過創(chuàng)建關(guān)鍵點、線,建立立柱截面,然后通過拉伸、鏡像、布爾運算等命令,生成了外殼體的參數(shù)化有限元模型。

(3) 網(wǎng)格的劃分。有限元模型采用SOLID92單元,并設(shè)定軸對稱選項,建立1/2軸對稱分析模型。用SizeControl對網(wǎng)格單元的大小進行了控制,用volume sweep(體掃描)的方式對立柱進行了網(wǎng)格劃分。

(4) 加載并求解。殼體外壁施加均勻壓力面載荷,殼體XOZ平面施加固定約束。加載和約束后的有限元模型見圖1。

圖1 施加約束后的外殼體有限元模型

(5) 后處理。求解后,提取外殼體的Y向最大變形和最大等效應(yīng)力作為后續(xù)結(jié)構(gòu)可靠性分析的輸出參數(shù)。

2.2 外殼體蒙特卡羅模擬

采用蒙特卡羅法[7-8],循環(huán)次數(shù)初選loop=500,輸入?yún)?shù)的分布類型如表2所示,進行蒙特卡羅模擬并輸出結(jié)果。

表2 輸入變量概率分布類型及參數(shù)

3 外殼體蒙特卡羅模擬結(jié)果分析

3.1 統(tǒng)計結(jié)果分析

選擇了工作壓力p、壁厚H這2個具有代表性的變量,觀察其函數(shù)取值的分布,Y向最大變形和最大等效應(yīng)力的均值變化曲線分別如圖2和圖3所示。中間的1條曲線表示函數(shù)的平均值,而上下2條曲線則表示了抽樣過程中向平均值趨近的過程,曲線越接近,表明平均值越趨向于收斂。從圖2和圖3中可看出,Y向最大變形和最大等效應(yīng)力的平均值均趨向收斂,說明本文采用500次循環(huán)模擬次數(shù)是足夠的。

圖2 UDMAX的均值變化曲線

從圖2和圖3中可得到Y(jié)向最大變形的平均值約為0.06 mm,標(biāo)準(zhǔn)差為0.008 mm,最大等效應(yīng)力的平均值為845.6 MPa,標(biāo)準(zhǔn)差為38.25 MPa。

圖3 最大等效應(yīng)力SMAX的均值變化曲線

3.2 輸出變量敏感性分析

輸出變量敏感性分析結(jié)果見圖4。其中,輸入?yún)?shù)對輸出參數(shù)影響水平在2.5%以下的歸為相對影響不大的因素;影響水平在2.5%以上的歸為有顯著影響的因素。

圖4 敏感性分析結(jié)果示意圖

從圖4可以看出,壓力載荷p是影響失效的最主要因素;其次是殼體厚度H。因此設(shè)計和制造過程中應(yīng)嚴(yán)格控制這些參數(shù)。殼體的可靠度對壓力載荷p的靈敏度為負(fù)值,表示壓力載荷p增大,結(jié)構(gòu)的等效應(yīng)力最大值增加,其強度的可靠性降低。

在設(shè)計和制造時,應(yīng)嚴(yán)格控制這些主要敏感參數(shù)。在經(jīng)濟成本允許的情況下,可考慮適當(dāng)增加外殼體的厚度H,以改善外殼體的強度和剛度,使其獲得更好的力學(xué)性能,保證外殼組件的承壓強度。

3.3 結(jié)構(gòu)可靠性分析

外殼體承受的最大壓力為140 MPa,屈服強度為1 030 MPa,最大等效應(yīng)力的平均值為845.6 MPa,標(biāo)準(zhǔn)差為38.25 MPa。因此,運用“pdprob”命令計算外殼體的最大等效應(yīng)力SMAX小于1 030 MPa的概率,在置信度為95%的情形下,最大等效應(yīng)力的可靠度均為100%。說明外殼體采用初步設(shè)計參數(shù)的結(jié)構(gòu)可靠,安全度高。

對四探頭自然伽馬測井儀外殼體做外壓密封試驗,在140 MPa壓力下,保持15 min,要求無塑性變形及泄露。壓力試驗結(jié)果表明,外殼體加壓、保壓、卸壓正常,驗證了初步設(shè)計的正確性。

4 結(jié) 論

(1) 采用采用蒙特卡羅模擬法,根據(jù)壓力外部條件初步設(shè)計的外殼體結(jié)構(gòu)的幾何尺寸、應(yīng)力等各參數(shù)視為隨機變量,直觀地顯示了設(shè)計參數(shù)對結(jié)構(gòu)應(yīng)力的影響程度,對認(rèn)識結(jié)構(gòu)的可靠性有重要作用,為結(jié)構(gòu)的參數(shù)設(shè)計及優(yōu)化分析提供了有效理論依據(jù)。

(2) 該方法對設(shè)計其他型號的測井儀外殼體具有適應(yīng)性,可在初步設(shè)計的前提下進行外殼體蒙特卡洛模擬,校驗外殼體的強度,若不合適,重新調(diào)整設(shè)計參數(shù),使設(shè)計參數(shù)達(dá)到最優(yōu)。

參考文獻:

[1] 池偉,劉東友. ZGM70型自然伽馬測井儀殼體強度設(shè)計與校核 [J]. 艦船防化,2003(3): 35-37.

[2] 嚴(yán)慧娟,岳愛忠,趙均,等. 測井儀耐壓鋼筒的應(yīng)力分析與強度校核 [J]. 機械工程師,2010(5): 36-37.

[3] 中子—γ密度測井及地層元素測井的蒙特卡羅模擬研究 [D]. 蘭州: 蘭州大學(xué)碩士學(xué)位論文,2013.

[4] 浦廣益. ANSYS Workbench 12基礎(chǔ)教程與實例詳解 [M]. 北京: 中國水利水電出版社,2010.

[5] 李兵,何正嘉,陳雪峰. ANSYS Workbench設(shè)計、仿真與優(yōu)化 [M]. 北京: 清華大學(xué)出版社,2008.

[6] 許艷玲,張保,徐曉東. 高壓球形容器強度的可靠性分析 [J]. 化學(xué)工程與裝備,2010(10): 98-100.

[7] 張建國,蘇多,劉英衛(wèi). 機械產(chǎn)品可靠性分析與優(yōu)化 [M]. 北京: 電子工業(yè)出版社,2008.

[8] 嚴(yán)慧娟,岳愛忠,趙均,等. 地層元素測井儀器結(jié)構(gòu)參數(shù)的蒙特卡羅數(shù)值模擬 [J]. 測井技術(shù),2012,36(3): 282-285.