懸鏈線構(gòu)件的有限元方法

王瀚，解德，2，3

1華中科技大學(xué)船舶與海洋工程學(xué)院，湖北武漢430074

2船舶與海洋水動力湖北省重點(diǎn)實驗室，湖北武漢430074

3高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海200240

懸鏈線構(gòu)件的有限元方法

王瀚1，解德1，2，3

1華中科技大學(xué)船舶與海洋工程學(xué)院，湖北武漢430074

2船舶與海洋水動力湖北省重點(diǎn)實驗室，湖北武漢430074

3高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海200240

懸鏈線構(gòu)件在船舶與海洋工程中是廣泛存在的一種結(jié)構(gòu)，無論是船舶碼頭旁靠、錨泊定位中的系泊線，還是各種系泊結(jié)構(gòu)物中各種深海懸鏈線結(jié)構(gòu)（如深海柔性立管、深海錨泊系統(tǒng)、復(fù)合材料懸鏈線立管等）的系泊力響應(yīng)，對整個結(jié)構(gòu)物都顯得非常重要。其中，由于懸鏈線構(gòu)件具有很強(qiáng)的幾何非線性，其受力分析的精度對其運(yùn)動響應(yīng)的分析尤為重要。根據(jù)懸鏈線構(gòu)件靜力分析，得到懸鏈線單元的柔度矩陣，從而進(jìn)一步推導(dǎo)出懸鏈線單元的剛度矩陣和單元節(jié)點(diǎn)力向量，最后以商業(yè)有限元軟件ABAQUS的用戶自定義單元（UEL）為工具，開發(fā)出懸鏈線單元。結(jié)果顯示：運(yùn)用此單元得到的簡單懸鏈線結(jié)構(gòu)的數(shù)值解與解析解完全一致，同時也基本符合實驗結(jié)果，為解決各種復(fù)雜懸鏈線構(gòu)件問題的直接計算法提供了一種新的思路。

系泊系統(tǒng)；懸鏈線單元；幾何非線性；實驗驗證

0 引 言

21世紀(jì)是海洋的世紀(jì)，海洋中蘊(yùn)含了大量的資源，如今海洋開發(fā)正在向深海發(fā)展［1］。僅就石油資源而言，目前進(jìn)行的海洋石油開發(fā)主要模式為采用浮式或者固定式的各種適用于深海開發(fā)的新型海上石油生產(chǎn)處理裝置。其中主要包括半潛式平臺、張力腿平臺（TLP）、單柱式平臺（Spar）和浮式生產(chǎn)儲油輪（FPSO）等系泊結(jié)構(gòu)物。這些系泊結(jié)構(gòu)物中的立管系統(tǒng)、錨泊系統(tǒng)等包含有大量的懸鏈線構(gòu)件。對于立管系統(tǒng)而言，無論是傳統(tǒng)的柔性立管、頂張力立管［2］，還是適用于高溫、高壓工作環(huán)境的鋼懸鏈線立管［3-4］，亦或是結(jié)合上述各自優(yōu)點(diǎn)還處于研究階段的復(fù)合材料懸鏈線立管［5-6］，以及系泊結(jié)構(gòu)物中進(jìn)行海上定位的錨泊系統(tǒng)［7-11］，由于受風(fēng)、浪、流等惡劣海洋環(huán)境的影響，它們要抵御外界的各種載荷以保證系泊結(jié)構(gòu)物能正常工作運(yùn)轉(zhuǎn)。而對船舶而言，當(dāng)進(jìn)行碼頭旁靠或者水鼓系泊等各種停泊狀態(tài)的分析時，懸鏈線構(gòu)件都是重點(diǎn)關(guān)注的對象。因此，分析懸鏈線構(gòu)件的受力無論是對海洋平臺還是對船舶整個結(jié)構(gòu)的運(yùn)動響應(yīng)都起著至關(guān)重要的作用。懸鏈線構(gòu)件的受力主要表現(xiàn)為：在幾何上為大變形，具有很強(qiáng)的幾何非線性；懸鏈線的線長較大，其自重垂度不容忽略。針對懸鏈線的這些特點(diǎn)，必須選擇正確的有限元模型對其進(jìn)行非線性分析。Peuker等［12］將柔性管的每個組分都作為獨(dú)立元素進(jìn)行分析，并使其滿足平衡條件和連續(xù)條件，建立了理論解析模型，但該模型無法計算立管的整體剛度值。Féret等［13］提出在柔性立管的總體分析中采用集中質(zhì)量法，把柔性立管劃分為有限數(shù)目的離散構(gòu)件，將連續(xù)分布在柔性立管長度方向上的各種載荷等效集中在所建模型系統(tǒng)的節(jié)點(diǎn)處。對于錨泊系統(tǒng)而言，最常用的方法是采用包括非線性分析的通用有限元軟件（如，ANSYS軟件［14］）進(jìn)行計算，然而軟件的單元庫中并不包含索單元。經(jīng)過多年的發(fā)展，目前主要形成了3種主要數(shù)值模型［15］：用于靜力分析的懸鏈線模型和用于動力分析的集中質(zhì)量模型、細(xì)長桿模型。

本文將基于靜力分析的懸鏈線模型推導(dǎo)出懸鏈線單元［16］，其可以根據(jù)一根索的任何一個參數(shù)值（如，索端張拉力或下垂高度）為已知條件來確定索的狀態(tài)。并將進(jìn)行相對應(yīng)的實驗［17］，用結(jié)果證明在商業(yè)有限元軟件ABAQUS中運(yùn)用該單元計算可以獲得與解析解完全一致的結(jié)果，并且能與實驗結(jié)果相吻合，以證明所推導(dǎo)出的單元的準(zhǔn)確性。

1 懸鏈線的靜力解析解

圖1所示為無彈性懸鏈線段上任意一小段的靜力圖，其中W為系泊線單元長度的重量，dl為系泊線單位長度，T為系泊線的張力，dT為張力在dl上的增量，θ為張力T的方向與水平線的夾角，dθ為θ在dl上的增量。沿水平和垂直方向建立平衡方程，可以得到：

圖2所示為無彈性懸鏈線段靜力圖。其中：Ta，θa與Tb，θb為懸鏈線段ab兩端的張力及其傾角；T0為Ta和Tb的水平分量；ab長為l；a端和b端的垂直距離為y，水平距離為x。T0與Ta，Tb的關(guān)系如下：

圖1 無彈性懸鏈線靜力圖Fig.1 Statics graph of inelastic catenary

圖2 無彈性懸鏈線單元靜力圖Fig.2 Statics graph of inelastic catenary element

由式（1）～式（5）可推導(dǎo)出懸鏈線段拉力的水平分量T0與其坐標(biāo)的關(guān)系：

2 懸鏈線單元的剛度矩陣和節(jié)點(diǎn)力向量

將圖2所示的懸鏈線段視為一個無彈性懸鏈線單元，如圖3所示，由力的平衡條件，得

因此，式（6）和式（7）可改寫為：

對式（8）和式（9）進(jìn)行微分，得

所以，可以得到

則剛度矩陣

式中，detF=f11f22-f12f21。

所以，單元的切線剛度矩陣為

由于ABAQUS是一個通用的非線性有限元程序，因此它的求解過程是迭代法。因此結(jié)構(gòu)的平衡方程為

式中：ΔU為節(jié)點(diǎn)增量位移矩陣；P為各種外力形成的等效節(jié)點(diǎn)力矩陣。并以此為基礎(chǔ)編寫了二維問題的懸鏈線ABAQUS-UEL子程序。

3 懸鏈線單元的數(shù)值驗證

以1根懸鏈線為例，鏈重W=100 kg/m，總長l=480 m。圖4與圖5分別為水深100 m和水深145 m時解析解與ABAQUS數(shù)值解的對比圖。

圖4 水深100 m時結(jié)果對比Fig.4 The results contrast when depth is 100 m

圖5 水深145 m時結(jié)果對比Fig.5 The results contrast when depth is 145 m

可以看到，解析解與ABAQUS所求的數(shù)值解完全一致。

4 懸鏈線單元的實驗驗證

為進(jìn)一步驗證該單元的正確性，實驗驗證顯然必不可少。鐵鏈參數(shù)如下：

鏈重W/（N·mm-1） 0.002 844

總長l/mm 1 740

分別對500，1 000和1 500 mm這3個不同跨距下的懸鏈線進(jìn)行了實驗驗證，與實驗結(jié)果對比后的結(jié)果如圖6～圖8所示，圖中帶圓圈的線條為利用MATLAB得到的解析解。

圖6 跨距500 mm形狀對比Fig.6 Shape comparison when span is 500 mm

圖7 跨距1 000 mm形狀對比Fig.7 Shape comparison when span is 1 000 mm

圖8 跨距1 500 mm結(jié)果對比Fig.8 Shape comparison when span is 1 500 mm

利用測力計測得的最大拉力與解析解和數(shù)值解基本一致，具體結(jié)果如表1～表3所示。

表1 跨距500 mm結(jié)果對比Tab.1 The results contrast when span is 500 mm

表2 跨距1 000 mm結(jié)果對比Tab.2 The results contrast when span is 1 000 mm

表3 跨距1 500 mm結(jié)果對比Tab.3 The results contrast when span is 1 500 mm

可以看到，利用MATLAB得到的解析解與實驗中的鐵鏈完全吻合，而通過測力計所得到的最大拉力以及鏈條的垂向距離也與解析解以及通過ABAQUS得到的數(shù)值解基本一致，從而進(jìn)一步驗證了此單元的精確度。

5 結(jié) 論

本文對懸鏈線構(gòu)件的有限元分析進(jìn)行了基礎(chǔ)性的研究，并在簡化模型中分析了其受力情況，主要結(jié)論如下：

1）計算出了懸鏈線的剛度矩陣，利用ABAQUS開發(fā)了一種懸鏈線單元，并且利用ABAQUS得到的數(shù)值解與解析解完全一致。

2）進(jìn)行了簡單的實驗，并與解析解和數(shù)值解進(jìn)行對比，結(jié)果基本吻合，證明了此懸鏈線單元的準(zhǔn)確性。

3）本單元為解決懸鏈線構(gòu)件運(yùn)動等復(fù)雜問題的直接計算法提供了一種新的思路。但目前為止的計算及實驗都是簡化模型，還需對單元進(jìn)行進(jìn)一步的開發(fā)以使其能夠適應(yīng)多種載荷下以及在三維情況下的運(yùn)動響應(yīng)。

［1］ BAARHOLM R，PALAZZO F G.Hybrid verification of a DICAS moored FPSO［C］//Proceedings of the 14th International Offshore and Polar Engineering Confer?ence.Toulon，F(xiàn)rance：ISOPE，2004：307-314.

［2］ 姜豪，楊和振，劉昊.深海非粘結(jié)柔性立管簡化模型數(shù)值分析及實驗研究［J］.中國艦船研究，2013，8（1）：64-72. JIANG Hao，YANG Hezhen，LIU Hao.Experimental and numerical analysis of a new simplified model for the deepwater unbonded flexible risers［J］.Chinese Journal of Ship Research，2013，8（1）：64-72.

［3］ MARTINGS C A，HIGASHI E.A parametric analysis of steel catenary risers：fatigue behavior near the top［C］//Proceedings of the 10th International Offshore and Polar Engineering Conference.Seattle，USA：ISPOE，2000：54-59.

［4］ 黃維平，白興蘭，李華軍.國外深水鋼懸鏈線立管研究發(fā)展現(xiàn)狀［J］.中國海洋大學(xué)學(xué)報，2009，39（2）：290-294. HUANG Weiping，BAI Xinglan，LI Huajun.State of the art of research and development of overseas deep?water steel catenary risers［J］.Periodical of Ocean Uni?versity of China，2009，39（2）：290-294.

［5］ OCHOA O O，SALAMA M M.Offshore composites：transition barriers to an enabling technology［J］.Com?posites Science and Technology，2005，65（15）：2588-2596.

［6］ 沈欽雄，楊和振，朱云.深海復(fù)合材料懸鏈線立管基于可靠度的優(yōu)化設(shè)計［J］.中國艦船研究，2014，9（5）：77-84. SHEN Qinxiong，YANG Hezhen，ZHU Yun.Reliabili?ty-based design optimization of deep-water composite catenary risers［J］.Chinese Journal of Ship Research，2014，9（5）：77-84.

［7］ 陳徐均，崔維成，沈慶.對稱式布置錨鏈系統(tǒng)的線性化處理［J］.海洋工程，2002，20（1）：75-79. CHEN Xujun，CUI Weicheng，SHEN Qing.Lineariza?tion for symmetrically collocated mooring system［J］. The Ocean Engineering，2002，20（1）：75-79.

［8］ 潘斌，高捷，陳小紅，等.浮標(biāo)系泊系統(tǒng)靜力計算［J］.重慶交通學(xué)院學(xué)報，1997，16（1）：68-73. PAN Bin，GAO Jie，CHEN Xiaohong，et al.Static cal?culation of buoy mooring fast［J］.Journal of Chongqing Jiaotong Institute，1997，16（1）：68-73.

［9］ 胡靈斌，唐軍.懸鏈線方程的求解及其應(yīng)用［J］.船舶，2004（1）：17-20. HU Lingbin，TANG Jun.Solemnest and application of catenary's curve equation［J］.Ship and Boat，2004（1）：17-20.

［10］ 張火明，范菊，楊建民.復(fù)雜系泊系統(tǒng)靜力特性快速計算方法研究［J］.艦船科學(xué)技術(shù)，2007，29（3）：113-118. ZHANG Huoming，F(xiàn)AN Ju，YANG Jianmin.Investi?gation on quick computation method of the static char?acteristics of complex mooring system［J］.Ship Sci?ence and Technology，2007，29（3）：113-118.

［11］ 郝春玲，滕斌.不均勻可拉伸單錨鏈系統(tǒng)的靜力分析［J］.中國海洋平臺，2003，18（4）：18-21，33. HAO Chunling，TENG Bin.Static analysis for a non-uniform flexible mooring cable system［J］.China Offshore Platform，2003，18（4）：18-21，33.

［12］ PEUKER M，KOKKINOWRACHOS K，GIESE K. Development of flexible riser for floating offshore pro?duction［C］//Offshore Technology Conference.Hous?ton，Texas：OTC，1987.

［13］ FéRET J J，BOURNAZEL C L.Calculation of stress?es and slip in structural layers of unbounded flexible pipes［C］//Proceeding of First Offshore Mechanics and Arctic Engineering Specialty Symposium.New Orleans，Louisiana，1986：311-319.

［14］ 劉元丹，熊治國，劉敬喜，等.基于AQWA的旁靠油輪水動力相互作用研究［J］.艦船科學(xué)技術(shù)，2012，34（5）：13-17，22. LIU Yuandan，XIONG Zhiguo，LIU Jingxi，et al. Analysis on hydrodynamic interaction between side by side mooring tankers［J］.Ship Science and Tech?nology，2012，34（5）：13-17，22.

［15］ 王冰.基于細(xì)長桿理論的系泊纜索靜力及動力分析方法研究［D］.哈爾濱：哈爾濱工程大學(xué)，2013.

［16］ 楊佑發(fā)，白文軒，郜建人.懸鏈線解答在斜拉鎖數(shù)值分析中的應(yīng)用［J］.重慶建筑大學(xué)學(xué)報，2007，29（6）：31-34. YANG Youfa，BAI Wenxuan，GAO Jianren.Applica?tion of catenary solution in numerical analysis of stayed cables［J］.Journal of Chongqing Jianzhu Uni?versity，2007，29（6）：31-34.

［17］ KIM M H，KOO B J，MERCIER R M，et al.Vessel/ mooring/riser coupled dynamic analysis of a tur?ret-moored FPSO compared with OTRC experiment［J］.Ocean Engineering， 2005， 32（14/15）：1780-1802.

［責(zé)任編輯：喻菁］

Analysis and experimetal verification of the catenary component with the finite element method

WANG Han1，XIE De1，2，3

1 School of Naval Architecture and Ocean Engineering，Huazhong University of Science and Technology，Wuhan 430074，China

2 Hubei Key Laboratory of Naval Architecture and Ocean Engineering Hydrodynamics，Wuhan 430074，China

3 Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration，Shanghai 200240，China

The catenary component is widely applied on ships and ocean engineering structures.Whether it is for the mooring line of dockside ships,the anchor moored positioning,or various mooring structures with different deep-sea catenary structures(such as deep-sea flexible risers,deepwater mooring systems, deep-sea steel catenary risers,etc.),as well as the mooring force response of the entire structure are of par?ticular significance.To be more specific,the catenary component in the mooring system is obviously a non?linear structure geometrically,while the accuracy of the stress analysis of the catenary is particularly impor?tant.The stiffness matrix and node force vector of the catenary element can be gained based on the flexibili?ty matrix obtained from static analytical solutions of the catenary.Finally,the catenary element is estab?lished by using the UEL(User Defined Element)from commercial finite element software ABAQUS.The re?sult is seen to be fully consistent with the analytical solution and the experimental results.In brief,it induc?es creativity to the direct calculation method of solving various complex catenary component problems.

mooring system；catenary element；geometrical nonlinearity；experimental verification

U663

A

10.3969/j.issn.1673-3185.2015.06.006

http://www.cnki.net/kcms/detail/42.1755.TJ.20151110.1025.012.html期刊網(wǎng)址：www.ship-research.com

王瀚，解德.懸鏈線構(gòu)件的有限元方法［J］.中國艦船研究，2015，10（6）：34-38. WANG Han，XIEDe.Analysisandexperimetalverificationofthecatenarycomponentwiththefiniteelementmethod［J］. Chinese Journal of Ship Research，2015，10（6）：.34-38

2015-03-18 < class="emphasis_bold"> 網(wǎng)絡(luò)出版時間：

時間：2015-11-10 10:25

王瀚，男，1989年生，碩士生。研究方向：系泊結(jié)構(gòu)物運(yùn)動及響應(yīng)。E-mail：327684188@qq.com解德（通信作者），男，1964年生，博士，教授，博士生導(dǎo)師。研究方向：船舶與海洋結(jié)構(gòu)物設(shè)計制造