基于ABAQUS的鉆井平臺鎖緊系統有限元分析

王建才，楊立東，田 俊，馬樹宇

基于ABAQUS的鉆井平臺鎖緊系統有限元分析

王建才1，楊立東2，田 俊1，馬樹宇1

（1．渤海裝備遼河重工有限公司，遼寧盤錦124010；2．盤錦凱瑞能源有限公司，遼寧盤錦124010）①

鎖緊系統是自升式鉆井平臺的關鍵部件，其可靠性對于平臺自身安全至關重要。利用傳統的力學公式和計算方法很難驗算鎖緊系統各部件應力的大小及分布。為驗證鎖緊系統結構的合理性，以有限元軟件ABAQUS將鎖緊系統各部件劃分為規則的六面體網格單元，并以非線性模塊作為計算分析平臺對其進行強度分析，根據計算結果進行強度校核。結果表明，鎖緊系統各部件均能滿足強度要求，同時也為開發鎖緊系統相關系列產品提供了可靠數據。

鉆井平臺；鎖緊系統；非線性；強度校核

自升式鉆井平臺是近年來近海油氣勘探開發中使用最廣泛的海洋工程裝備之一，平臺鎖緊系統是自升式鉆井平臺的關鍵核心部件。對于鎖緊裝置而言，安全可靠性設計是極其重要的環節，研制性能良好的鎖緊裝置對于平臺的安全穩定起關鍵作用。因此，需要設計出合理的結構，選取合適的零部件材料，采取適宜的計算方法，來驗證鎖緊系統在實際作業工況中能否達到使用要求。

1　結構組成

同以前的鎖緊機構相比，新研制的鎖緊系統無論在結構和承載能力上都有了較大的提升和改進。鎖緊系統的整體鋼架采用拼接式結構，內部零件采取鍛、鑄結合的工藝，不但結構簡單，而且焊接組裝等工序較為方便，降低了裝配后的成本和風險，也使得設計人員在研發初期就能對鎖緊系統的整體性有進一步的理性認識［1］。新研制的鎖緊系統結構如圖1所示。

圖1　鎖緊系統結構

2　鎖緊系統有限元分析理論

六面體單元是擁有局部和整體兩個坐標系的三維單元，可用線性函數描述［2-3］。其基本參數如應力、應變等與彈性模量E和泊松比μ有關聯，每個六面體單元有8個節點且每個節點有3個自由度。分析計算的基本步驟如下。

2．1 結構離散化

將整體離散成有限個微分體dxdydz，任取一個來探討應力應變之間的關系。整體在外部載荷的作用下，體內任一點的應力狀態由6個應力分σx、σy、σz、τxy、τyz、τzx來表示。其中，σx、σy、σz表示各坐標軸上的正應力，τxy、τyz、τzx表示各平面內的剪應力，微觀模型如圖2所示。

圖2　微六面體單元的力學模型

2．2 單元分析

由彈性力學可知，彈性體內任取一點處的應力、應變和位移間的關系可根據以下公式描述出：

1） 位移和應變間的關系。

式中：（ε）為單元內任一點處的應變矩陣；εx、εy、εz為任一點處的正應變；γxy、γyz、γzx為任一點處的剪應變；（B）為幾何矩陣；（δe）為單元節點位移矩陣。

2） 應變和應力間的關系。

式中：（σ）為單元內任一點處的應力矩陣；（D）為彈性矩陣；E為單元所屬材料的彈性模量；μ為單元所屬材料的泊松比。

3） 單元剛度方程。

根據虛功原理建立單元剛度方程，單元節點力和節點位移間的關系即可在方程中顯現出來。

式中：（Fe）為單元節點載荷矩陣；（Ke）為單元剛度矩陣。

2．3 整體分析

在單個單元剛度方程的基礎上，通過集成得到整體剛度方程，在此基礎上進行整體結構分析。集成的原則是相鄰單元間在共用節點處的位移也相同。

式中：（F）為整體節點載荷矩陣；（K）為整體剛度矩陣；（δ）為整體節點位移矩陣。

2．4 求解節點位移及單元的應變與應力

求解方程式（4）。根據得出的節點位移，代入式（1）～（2）并計算，即可得到相應單元節點的應變和應力。

3　鎖緊系統的有限元模型建立

有限元軟件自身的建模功能與專業的三維繪圖軟件對比而言，其功能相對單一，但均與各繪圖軟件有著良好的接口程序，因此鎖緊系統實體模型采用U G進行建模，之后導入ABAQUS軟件中將模型進行規則的六面體網格劃分，以提高分析計算結果的精度。確保模型在分析計算過程中能夠達到所需精度的前提下，將實體結構中的部分細小倒角，螺栓孔等對整體強度影響較小的結構進行簡化處理，以縮短不必要的計算時間［4］。以ABAQUS軟件的非線性分析模塊作為分析平臺，對鎖緊系統加載極限載荷的工況下進行極限強度分析，并對計算結果進行必要的后處理。

在對鎖緊系統進行強度分析的同時，將其內部結構如鎖緊齒條、壓塊等鎖緊裝置中的關鍵部件進行建模，盡量做到與實際情況相一致，從而減少因外部條件改變造成的計算誤差。因此在分析過程中采取的分析類型為邊界非線性。此類問題在物體相互接觸邊界位置以及接觸面上的應力大小和分布事先是未知的，將整個問題求解完成后才能確定各接觸面的受力情況以及位移變化。

前處理在整個有限元分析過程中所花費的時間和精力是最多的，占整個過程的70%左右。而前處理過程中，有限元模型網格質量對于求解的精度和準確度有著決定性的影響。前處理的效率高低以及質量好壞對整個有限元分析重要性很高。在接觸分析中，同等條件下，六面體網格要比四面體網格仿真的結果更為準確，在接觸起作用的同時也能使分析更容易收斂，這也是本次鎖緊系統采用六面體單元為主分析的原因。

鎖緊系統實體模型選用8節點六面體單元對其進行規則的網格劃分，壓塊加載端面采用主控制點與MPC剛性梁單元結合的形式進行連接，以模仿極端載荷的加載形式。樁腿齒條施加對稱的固定約束方式，鎖緊系統內部各接觸件如導軌架與壓塊，壓塊與鎖緊齒條，鎖緊齒條與樁腿件均建立接觸對，以模仿實際情況下兩兩物體的接觸狀態。

因鎖緊裝置為左右對稱結構，因此在有限元建模過程中只建立1／2結構，這樣不但節省建模周期，也能簡化計算過程，節省前處理的時間。有限元模型如圖3所示。

圖3　鎖緊系統有限元模型

本次分析選取極限載荷工況進行應力分析，極端載荷為7．2×104k N。因為模型成對稱性模型，所以計算采用半模型進行分析計算，對稱面上施加對稱邊界約束，加載取3．6×104k N進行分析計算。

鎖緊齒條和楔塊的材質均為40CrNiMo，其綜合力學性能良好，有高的強度、韌度和良好的淬透性和抗過熱的穩定性，使其廣泛用于制造強度高、塑性好的重要零部件。樁腿齒條采用A517Q，導軌架材料選為E550，其基本屬性如表1。

表1　材料基本屬性

本次有限元分析，在對模型附著材料屬性時的單位制是mmNMPa，其原因是有限元軟件中沒有固定的單位制，因此在對模型附著材料屬性時，其數值本身所含有的特定單位將被軟件作為默認值［5］。

4　鎖緊裝置的強度分析

有限元位移法的分析方程可歸納為：

式中：K為結構剛度矩陣；U為位移列向量；F為外載荷列向量。

在分析過程中用到了多點約束即MPC。MPC通過多點約束的形式將不同的節點進行耦合連接，MPC內部合力為零，其不會增加結構的剛度，并且運用梁單元將各不同部件連接，其計算過程仍可看作是線性的，此舉節省了大量的前后處理時間。MPC單元即多點約束剛性梁有如下特性：

式中：ui為節點i處的平動矢量；ei為節點i處的轉動矢量；φ為當前坐標系下的旋轉角。

由式（6）可得到3個附加方程來處理實體單元與梁單元間不同自由度數之間耦合的問題［6］。

根據材料力學方面的理論可知，在一定厚度的范圍內，板材厚度對材料的許用應力是沒有影響的。但隨著板材厚度的增加以及制造工藝的改變可能使板材內部存在缺陷等，就會造成材料的力學性能隨著板材厚度的增加而下降現象的普遍發生。為此，在鎖緊系統的設計過程中，這方面的因素也在考慮范圍之內，從而加強了設計的合理性和可行性［7-8］。根據《海上移動平臺入級規范》，取材料的一般局部安全系數γm＝1．1。鎖緊系統關鍵部件鎖緊齒條、楔塊、導軌架的有限元分析應力云圖分別如圖4～6所示。

圖4　鎖緊齒條應力云圖

圖5　楔塊應力云圖

根據圖4鎖緊齒條應力云圖所示，在非線性計算分析中，由于楔塊與鎖緊齒條應力接觸造成鎖緊齒條后端棱角出現應力集中問題，對零件真實使用效果無直接影響，所以鎖緊齒條最大的合理計算應力值可取為720mPa，材料許用應力為800mPa，安全系數為1．11，符合規范要求。

根據圖5楔塊應力云圖所示，楔塊由于邊界條件中MPC的存在，造成了附近應力集中問題，對零件的真實使用效果無失效影響，楔塊的合理計算應力為720mPa，材料許用應力為800mPa，安全系數為1．11，符合規范要求。

圖6　導軌架應力云圖

根據圖6導軌架應力云圖所示，導軌架模型中存在90°尖角的緣故，造成了其附近應力集中現象的存在，此現象的存在對零件的真實使用效果無失效影響，導軌架計算應力為495mPa，材料許用應力為550mPa，安全系數為1．11，符合規范要求。

5　結論

1） 本鎖緊系統達到了結構強度要求，并通過了船級社的相關認證，實現了預期的設計目的，同時為相關零部件的研發設計提供了真實可信的參考數據。

2） 鎖緊系統在施加極限載荷狀態下依然可以滿足承載需求，可達到鎖緊系統最初設計的各項性能指標，而且結構簡單、工藝新穎、生產成本低，并且其獨特的結構還能消除裝置在加載時所出現的局部應力集中現象。

3） 六面體單元網格所建立起來的鎖緊系統有限元模型，在計算過程中節省了大量的計算分析時間，且應力和變形計算結果的準確性也得到了相應的提高，不但獲得了各部位均勻過渡的有限元云圖，還準確地計算出鎖緊系統各關鍵部件的應力閾值，為驗證鎖緊系統的合理性提供了可靠依據。

［1］ 劉茜，南傳立，萬淑敏．基于有限元分析的壓實機壓實效果研究［J］．天津大學學報，2004，37（6）：501-504．

［2］ 夏盛來，何景武．基于工程應用的有限元網格劃分研究［J］．飛機設計，2008，28（4）：10-13．

［3］ 董梅．有限元法在結構件強度分析中的應用［J］．指揮控制與仿真，2007，29（4）：114-117．

［4］ 平安，甘娥忠，于亞婷．有限元法—原理、建模及應用［M］．北京：國防工業出版社，2004．

［5］ 馬愛軍，周傳月，王旭．Patran和Nastran有限元分析專業教程［M］．北京：清華大學出版社．2005．

［6］ 申瑞民，許彥青，張同珍，等．基于多代理的智能型遠程教學環境研究［J］．計算機工程與應用，2002（4）：253-256．

［7］ 姚興佳，張晨晨，單光坤，等．基于ANSYS的風力發電機機艙底盤的強度分析［J］．沈陽工業大學學報，2008，30（1）：3841．

［8］ 朱愛民，其木格，鄭鑫．基于UG的輪轂結構可靠性分析［J］．機械工程與自動化，2010（1）：74-76．

Finite Element Analysis of Drilling Platform Locking System Based on A B A Q US

WANG Jiancai1，YANG Lidong2，TIAN Jun1，MA Shuyu1
（1．Bohai Equipment Liaohe Heauy Industry Com pany，Panjin124010，China；2．K airui Energy Co．，Ltd．，Panjin124010，China）

The locking system is a key component of offshore drilling platform，the reliability of the locking device is vital for the safety of platform．It is difficult to sim ulate its stress value and distribution accurately by using classicalmechanics．In order to verify the rationality of the locking structure，each component of the locking system ismeshed as hexahedral elements based on the A B A Q U S software，strength analysis of the locking device is finished based on the nonlinear analysis platform and checking the strength according to the calculation results．The results showed that，its strengthmeets the design requirements，at the same time it also provides the reliable data for the development of the locking system related products．

drilling platform；locking system；nonlinear；strength chec

TE951

A

10．3969／j．issn．10013842．2015．03．011

10013482（2015）03004804

①2014-10-12

王建才（1968-），男，甘肅武威人，高級工程師，主要從事石油鉆機和海洋工程裝備的研究，Email：wjclsnoi＠126．com。