平恩順，王 林，徐慶祥，李 楠，黃 其，汪 強(qiáng)

（中國石油集團(tuán)渤海鉆探工程有限公司井下技術(shù)服務(wù)分公司，天津 300283）

油氣工程

基于ANSYS Workbench的可變徑膨脹錐有限元分析

平恩順，王 林，徐慶祥，李 楠，黃 其，汪 強(qiáng)

（中國石油集團(tuán)渤海鉆探工程有限公司井下技術(shù)服務(wù)分公司，天津 300283）

根據(jù)彈塑性有限元理論，應(yīng)用非線性有限元分析軟件ANSYS Workbench模塊建立了可變徑膨脹錐實(shí)體模型，動態(tài)模擬了可變徑膨脹錐的膨脹過程，得到了可變徑膨脹錐的力學(xué)狀態(tài)圖，可為膨脹錐的結(jié)構(gòu)設(shè)計(jì)以及膨脹方案的優(yōu)化提供技術(shù)依據(jù)。

有限元分析；ANSYS Workbench；可變徑膨脹錐；膨脹力

隨著陸地鉆井逐漸向深層、超深層以及復(fù)雜地層發(fā)展，套管補(bǔ)貼、側(cè)鉆井完井和膨脹管鉆井技術(shù)的勘探開發(fā)日趨重要。膨脹管技術(shù)是將管柱下到井內(nèi)，以機(jī)械或液壓的方式使生產(chǎn)管柱內(nèi)徑膨脹使其發(fā)生永久性塑性變形，從而達(dá)到套管內(nèi)徑膨脹的工藝目標(biāo)?？勺儚脚蛎涘F作為膨脹管技術(shù)中的一個(gè)重要組成部分，可變徑膨脹錐的研發(fā)對套管補(bǔ)貼、側(cè)鉆井完井和膨脹管鉆井技術(shù)有著極其深刻的影響。膨脹管受力緩慢膨脹是一個(gè)復(fù)雜的大塑性變形問題，材料強(qiáng)化現(xiàn)象明顯，采用理想彈塑性模型不能反映材料在實(shí)際膨脹過程中對膨脹力的影響，由于膨脹過程中大塑性變形的復(fù)雜性，很難進(jìn)行準(zhǔn)確的數(shù)學(xué)推導(dǎo)，對可變徑膨脹錐的數(shù)值計(jì)算相對復(fù)雜，但可以借助有限元分析軟件對其進(jìn)行準(zhǔn)確的數(shù)值模擬，分析整個(gè)膨脹過程的受力情況。秦國明等[1]利用ANSYS軟件的LS-DYNA模塊動態(tài)模擬了膨脹套管的全過程，得到了膨脹力隨時(shí)間變化的曲線，并分析了在不同的工藝參數(shù)下的膨脹力變化情況。徐丙貴等[2]利用ANSYS軟件建立了膨脹套管在套損井修復(fù)應(yīng)用中的有限元模型，優(yōu)選模型的邊界條件，進(jìn)行模型的數(shù)值模擬分析。于洋等[3]根據(jù)彈塑性有限元理論，利用有限元數(shù)值模擬研究了Φ139.7 mm實(shí)體膨脹管的膨脹特性，探討了膨脹率、屈服強(qiáng)度、摩擦系數(shù)和膨脹錐錐角對膨脹力的影響規(guī)律。郭慧娟等[4]采用非線性有限元分析軟件ABAQUS對實(shí)體膨脹管不同錐角時(shí)的膨脹過程進(jìn)行了詳細(xì)的數(shù)值模擬，研究了膨脹錐錐角與膨脹壓力等的關(guān)系。唐興波等[5]利用有限元軟件對變徑膨脹工具的錐角角度，長度和厚度進(jìn)行了結(jié)構(gòu)優(yōu)化。以上文獻(xiàn)主要集中于對實(shí)體膨脹套管的有限元分析，而對可變徑膨脹錐的有限元分析較少。因此，本文應(yīng)用ANSYS Workbench模塊對可變徑膨脹錐實(shí)體模型進(jìn)行數(shù)值模擬分析，動態(tài)模擬了膨脹過程，得到了可變徑膨脹錐力學(xué)狀態(tài)圖，可為膨脹錐的結(jié)構(gòu)設(shè)計(jì)以及膨脹方案的優(yōu)化提供技術(shù)依據(jù)。

1 有限元方程

可膨脹管的塑性大變形膨脹過程，是一種既含有材料非線性，又含有幾何非線性的高度非線性行為。因此，在進(jìn)行有限元分析時(shí)，不但要采用非線性的應(yīng)變和位移之間的本構(gòu)關(guān)系,而且平衡方程也必須建立于變形后的狀態(tài)，以考慮變形對平衡的影響[6]。

本文采用更新的拉格朗日格式（Updated Lagrange Formulation）來建立平衡方程。更新的拉格朗日格式在求t+Δt時(shí)刻物體的平衡解時(shí)，選擇的參考構(gòu)形是最后一個(gè)已知的平衡狀態(tài)，即將t+Δt時(shí)刻作為參考構(gòu)形，去求t時(shí)刻物體的各個(gè)未知變量。這時(shí)的虛功方程為[7]：

式中：Sij-Kirchhoff應(yīng)力張量；δEij-Green 應(yīng)變張量；W-外力虛功；σij-Cauchy應(yīng)力張量；ui-位移分量；xi-坐標(biāo)分量。式中各變量的左上標(biāo)表示所處的時(shí)刻，左下標(biāo)表示參照的時(shí)刻。

增量形式的本構(gòu)方程為[8]：

將式（2）、（3）、（4）、（5）代入虛功方程式（1）后再線性化，最終可得到更新的拉格朗日格式有限元平衡方程[9]：

式中：［KL］、［KNL］-線性和非線性剛度矩陣；｛F｝-初應(yīng)力載荷列陣。

式中：t［Dep］-小變形情況下的彈塑性矩陣；t［σd］-由Cauchy應(yīng)力分量構(gòu)成的矩陣。

2 有限元模型建立

2.1 實(shí)體模型

利用非線性有限元分析軟件ANSYS Workbench模塊對膨脹過程進(jìn)行數(shù)值模擬，建立可變徑膨脹錐實(shí)體模型（見圖1）。在有限元模型中，膨脹錐為主面（剛體），膨脹管內(nèi)壁為從面（可變形體），計(jì)算屬于剛體-柔體的接觸問題，且存在摩擦和大變形。因此在建模及模擬過程中采用如下簡化假設(shè)[4]：

規(guī)?；酿B(yǎng)羊場中羊群數(shù)量龐大，一經(jīng)發(fā)現(xiàn)死亡病例時(shí)，應(yīng)及時(shí)檢查并明確死因，確保羊群的健康安全。因此，大腸桿菌病的檢測與診斷，可以通過采取羊的心頭血和肝組織位置，制成涂片放在顯微鏡下觀察，如果很容易觀察到大腸桿菌，便可以確診羊的死亡原因，并做好防治措施。之所以要對羊的死亡診斷透徹，是因?yàn)榇竽c桿菌病與羊快疫和猝狙的外觀病癥十分相似，而大腸桿菌病的死亡容易引起羊群的高發(fā)感染，因此要明確診斷羊的死因。并實(shí)施明確的治療方案。

（1）由于膨脹速度較慢，忽略膨脹過程中的熱效應(yīng)；

（2）采用V.Mises初始屈服條件；

（3）應(yīng)用運(yùn)動硬化法則；

（4）采用雙線性理想彈塑性應(yīng)力流動模型；

（5）采用庫倫摩擦模型，且摩擦只發(fā)生在膨脹錐與膨脹管的接觸界面，摩擦因數(shù)在膨脹過程中保持不變；

（6）假設(shè)膨脹錐為剛體，不發(fā)生任何變形，膨脹管為變形體，即接觸為剛?cè)峤佑|問題。

圖1 可變徑膨脹錐實(shí)體模型

表1 TWIP鋼力學(xué)性能表

2.2 實(shí)體模型材料性能參數(shù)

通過對TWIP鋼的力學(xué)拉伸性能等相關(guān)試驗(yàn)的確定[10]，得出其材料性能參數(shù)（見表1）。將得出的數(shù)據(jù)通過Workbench數(shù)據(jù)輸入到材料的力學(xué)性能中，做好力學(xué)分析的準(zhǔn)備工作。

2.3 有限元模型

膨脹管的膨脹過程是一種高度的非線性行為，將其作為三維剛?cè)峤佑|問題來研究，即假設(shè)膨脹管為柔性體，而膨脹錐為剛體。膨脹錐與膨脹管之間的摩擦采用庫侖摩擦模型。為了提高求解精度及單元對大變形的適應(yīng)能力，對作為彈塑性體的膨脹管采用三維高階20節(jié)點(diǎn)六面體單元（SOLID95）進(jìn)行掃掠網(wǎng)格劃分，然后采用三維接觸單元（CONTA174）對接觸表面劃分接觸單元。由于將膨脹錐視為剛體，因此直接應(yīng)用三維目標(biāo)接觸單元（TARGE170）對其進(jìn)行網(wǎng)絡(luò)劃分[8]。

邊界條件為固定基管后膨脹的一端，使軸向位移為零，環(huán)向、徑向自由，膨脹錐上表面有垂直方向的位移，其他兩個(gè)方向自由。研究可變徑膨脹在膨脹管為TWIP鋼的情況下的結(jié)構(gòu)薄弱環(huán)節(jié)和受力情況。

3 有限元數(shù)值模擬分析

圖2 可變徑膨脹錐力學(xué)狀態(tài)圖

整個(gè)過程中尾端的受力最大（見圖2），而最小的受力點(diǎn)來自膨脹塊的反面膨脹本體導(dǎo)向槽，在整個(gè)膨脹過程中的時(shí)間/空間狀態(tài)下膨脹塊邊緣的受力很大，膨脹塊兩側(cè)飛邊導(dǎo)致了曲線的不平滑，這是由于膨脹塊在到達(dá)膨脹軌道的最高點(diǎn)的時(shí)候膨脹塊兩側(cè)的飛邊引起的問題。

4 結(jié)論及建議

（1）建立了可變徑膨脹錐的實(shí)體模型，動態(tài)模擬了膨脹錐的膨脹過程,得到了可變徑膨脹錐力學(xué)狀態(tài)圖。

（2）整個(gè)過程中尾端的受力最大，而最小的受力點(diǎn)來自膨脹塊的反面膨脹本體導(dǎo)向槽，在整個(gè)膨脹過程中的時(shí)間/空間狀態(tài)下膨脹塊邊緣的受力很大，膨脹塊兩側(cè)飛邊導(dǎo)致了曲線的不平滑。

（3）可變徑膨脹錐膨脹過程的有限元數(shù)值模擬分析，可為后續(xù)的膨脹套管修復(fù)技術(shù)提供技術(shù)依據(jù)，建議加大該方面理論研究及先導(dǎo)性試驗(yàn)，并為套管修復(fù)技術(shù)的現(xiàn)場應(yīng)用提供可行性分析，確?，F(xiàn)場施工作業(yè)的安全可靠性。

［1］ 秦國明，何東升，張麗萍，等.基于ANSYS/LS-DYNA的實(shí)體膨脹管膨脹力分析［J］.石油礦場機(jī)械，2009，38（8）：9-11.

［2］ 徐丙貴，張燕萍，王輝，等.數(shù)值模擬法在膨脹套管修復(fù)套損井技術(shù)中的應(yīng)用［J］. 石油勘探與開發(fā)，2009，36（5）：651-657.

［3］ 于洋，周偉，劉曉民，等.實(shí)體膨脹管的膨脹力有限元數(shù)值模擬及其應(yīng)用［J］.石油鉆探技術(shù)，2013，（5）：107-110.

［4］ 郭慧娟，楊慶榜，徐丙貴，等.實(shí)體膨脹管數(shù)值模擬及膨脹錐錐角優(yōu)化設(shè)計(jì)［J］.石油機(jī)械，2010，38（7）：30-32.

［5］ 唐興波，李黔，劉永剛，等.膨脹管變徑膨脹工具結(jié)構(gòu)優(yōu)化設(shè)計(jì)［J］.石油礦場機(jī)械，2008，37（11）：23-25.

［6］ 王勖成，邵敏.有限單元法基本原理和數(shù)值方法［M］.北京：清華大學(xué)出版社，2002：483-564.

［7］ 楊玉英，徐偉力，金朝海，等.幾種典型件成形過程的彈塑性有限元數(shù)值模擬［J］.塑性工程學(xué)報(bào)，1999，6（4）：22-25.

［8］ 秦永和，付勝利，高德利.實(shí)體膨脹管膨脹后軸向位移有限元數(shù)值模擬［J］.石油鉆采工藝，2006，28（5）：1-3.

［9］ 付勝利，高德利，等.可膨脹管膨脹過程三維有限元數(shù)值模擬［J］.西安石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2006，21（1）：54-57.

［10］ 梁金豹.技術(shù)機(jī)會驅(qū)動的可變徑膨脹錐的創(chuàng)新設(shè)計(jì)［D］.天津∶河北工業(yè)大學(xué)機(jī)械工程學(xué)院，2012：46-52.

Finite element analysis of the variable diameter expansion cone based on ANSYS Workbench

PING Enshun，WANG Lin，XU Qingxiang，LI Nan，HUANG Qi，WANG Qiang
（Downhole Technology Service Company，CNPC Bohai Drilling Engineering Company Limited，Tianjin 300283，China）

According to the theory of elastoplasticity finite element analysis,the solid model of the variable diameter expansion cone was established based on nonlinear finite element analysis software ANSYS Workbench module.The expansion process of the variable diameter expansion cone was used to simulate dynamically.The mechanical state figure of the variable diameter expansion cone was achieved.Thus technical basis was laid for structural design ofexpansion cone and optimization of expansion programs.

finite element analysis；ANSYS Workbench；variable diameter expansion cone；expansive force

TE931

A

1673-5285（2017）07-0011-04

10.3969/j.issn.1673-5285.2017.07.003

2017－06-07

中國石油渤海鉆探分公司項(xiàng)目“可降解材料在封隔器元件中的應(yīng)用研究”，項(xiàng)目編號：2016JXJF-07；“電動液壓橋塞座封工具的研制”，項(xiàng)目編號：2017JXJF-06。

平恩順（1986-），工程師，2015年博士畢業(yè)于河北工業(yè)大學(xué)機(jī)械制造及其自動化專業(yè)，工學(xué)博士學(xué)位，現(xiàn)主要從事油氣田儲層增產(chǎn)措施改造方面的研究工作，郵箱：pingenshun@163.com。