基于異形梁模型的海洋柔性管纜防彎器數值模擬

席勇輝+閻軍+楊志勛+郭宏+屈衍+李陽+趙娜

作者簡介： 席勇輝(1986—)，男，湖南永州人，助理工程師，碩士，研究方向為計算力學，(Email)xyh1115@foxmail.com；

閻軍(1978—)，男，遼寧大連人，副教授，博士，研究方向為海洋工程柔性管纜結構設計分析，(Email)yanjun@dlut.edu.cn

Numerical simulation on bending stiffener of flexible marine

pipe/cable based on special shape beam model

XI Yonghui1， YAN Jun2， YANG Zhixun2， GUO Hong3，

QU Yan3， LI Yang3, ZHAO Na4

(1. CSIC Kunming Shipborne Equipment Research & Test Center, Kunming 650051, China;

2. State Key Laboratory of Structure Analysis for Industrial Equipment, Dalian University of Technology,

Dalian 116024, Liaoning, China; 3. CNOOC Research Center, Beijing 100027, China;

4. Design Co. of CNOOC Offshore Oil Engineering Co., Ltd., Tianjin 300451， China）

Abstract： To prevent the structural failure of flexible marine pipe/cable due to the too large local curvature generated by selfgravity and environment load effects, the bending stiffeners are installed at the joints which can increase the bending stiffness of pipe/cable and make the curvature distribution uniform. Through the analysis on the bending stiffeners mechanical properties and working conditions, the calculation difficulties and the feasibility and effect of several kinds of numerical modeling methods are described, and the advantages of building finite element model in beam elements are further discussed in detail. According to the characteristics of plane beam elements, the bending stiffener structure is analyzed by equivalent bending stiffness. So a simple and efficient numerical analysis method is given. A numerical example is shown to evaluate the applicability of the calculation method.

Key words： flexible marine pipe/cable; bending stiffener; special shape beam; bending stiffness; numerical simulation

1引言

海洋柔性立管、動態海纜和臍帶纜等管線為使用方便采用螺旋纏繞型結構，在滿足抗拉伸能力的同時保持一定的易彎曲性，能夠承受一定的撓度或曲率.［1］這種多材料復合、非黏結纏繞的結構形式使得管纜設計過程中的力學分析和計算難度加大.由于管纜跟浮體的連接為固接而非鉸接形式，當管纜在自重或環境載荷作用下發生彎曲時，連接浮體處管纜會受到拉彎作用，其截面角度發生急劇變化［2］，使管纜的曲率超出其承受范圍而發生結構破壞.因此，在管纜與浮體連接的部位須配置保護構件［3］，使管纜在拉彎作用下的曲率也能夠保持安全狀態.

防彎器是套在管纜上的一個由聚合物材料制造的錐形構件，具有良好的彈性和變形能力［34］，在海洋工程中作為防止管纜過度彎曲保護的主要構件.防彎器由于安裝方便、要求空間小且在動態應用時對管纜有更好的保護效果，使其在海洋工程中得到廣泛使用［5］，如連接管纜與其他剛性設施（如浮體、ROV和水下井口等）或在其他截面角度易發生劇變處.管纜結構計算難度大，且與防彎器接觸并協同變形，為防彎器在安裝條件下的分析計算帶來困難.

隨著海洋油氣生產裝備的技術發展，國外關于防彎器方面的研究已較為成熟，但主要設計和制造技術由少數幾家大型海洋工程裝備制造公司掌握.國內海洋柔性管纜研發［69］取得一定進展，但對防彎器等附件研究［1011］還相對不足，目前僅制造出防彎器試驗樣品.本文從防彎器和管纜力學性能出發，提出采用異形梁單元模擬分析其主要力學響應的方法，并建立計算模型，對某防彎器實例進行設計和分析，驗證本文設計分析方法的適用性.

2防彎器力學性能特點

防彎器由聚合物材料組成，具有良好的變形能力.其幾何形狀和尺寸可根據需要進行設計制作，受工藝限制通常有圖1所示的幾種不同形式.FPSO和半潛式平臺等浮體受海洋環境影響產生各種運動，使管纜與浮體之間有相對位移和角度變化，同時管纜自身重力使管纜始終受到拉力作用，此時防彎器段的局部受力狀態可等效為一端與剛性面固接，另一端受斜拉載荷作用，見圖2.在此受力狀態下，管纜所受彎矩大小沿軸向方向呈非線性變化，并且越靠近固定端彎矩越大.

（a）單段結構（b）兩段結構（c）三段結構圖 1不同結構防彎器形式

Fig.1Different structures of bending stiffeners

圖 2防彎器與管纜所受載荷作用情況

Fig.2Load on bending stiffener and pipe/cable

在上述彎矩下，管纜的曲率與其自身截面彎曲剛度有關.在沒有防彎器保護時，靠近固定端處管纜由于受到的彎矩最大［12］，而纜/管體的截面彎曲剛度一般不足以使其曲率保持在安全限定范圍，并且只是局部曲率過大［13］，見圖3(a).防彎器的作用可視為在危險段管纜上附加一定的截面彎曲剛度，使此段的管纜曲率不會出現局部過大，見圖3(b).

(a)無防彎器(b)有防彎器圖 3在有、無防彎器時管纜的不同形態和曲率分布

Fig.3Different shapes and curvature distributions of

pipe/cable with and without bending stiffener

管纜的變形可由管纜上的彎矩分布和彎曲剛度分布計算得到.若對某防彎器的初始特征參數進行調整，即改變附加在管纜上的彎曲剛度后，在管纜端部拉力T作用下的管纜變形會發生改變，而相同管纜端部拉力T在不同的管纜形態下對管纜產生的彎矩分布也不同.彎矩分布改變又會影響管纜，因此管纜的變形和彎矩分布為一對相互影響的耦合變量［1415］，對防彎器進行設計時就需要進行管纜變形和彎矩分布的迭代計算.防彎器本身不直接承受外載荷，而是與管纜接觸，通過限制管纜變形時受到的反作用力產生自身變形.防彎器接觸區域管纜的曲率大小直接影響防彎器變形情況，因此在進行防彎器的力學分析時需要同時考慮管纜和防彎器.

3數值計算方法分析

由以上防彎器特點分析可知，用理論計算方法無法同時解決管纜與防彎器計算時的諸多迭代和耦合問題，一般需用有限元法進行數值模擬分析.在建立管纜和防彎器的模型后在管纜端部處施加拉力載荷T，并迭代計算得出管纜的曲率分布.

防彎器是簡單的均勻材料組成的實體結構，其模型見圖4(a).防彎器載荷施加需通過管纜傳遞，而管纜為復雜的多層螺旋纏繞結構［16］，見圖4(b).

(a)防彎器結構,m

(b)管纜結構

圖 4防彎器有限元模型和管纜結構

Fig.4Finite element models of bending stiffener and

pipe/cable structure

在有限元建模時可能需要模擬的結構包括內部大量獨立單元、外層鎧裝和特殊carcass層等，見圖5.若對管纜和防彎器進行完全的結構建模計算，需消耗大量計算資源和時間解決管纜變形迭代和內部單元之間接觸摩擦問題，甚至很難得到可用的計算結果.因此，對管纜和防彎器模型進行簡化處理，使其同時能夠滿足模擬管纜和防彎器的力學特點，以得到正確的變形結果.

(a)內部結構(b)鎧裝層結構(c)carcass結構

圖 5管纜構件有限元模型

Fig.5Finite element models of pipe/cable components

彎矩和截面彎曲剛度是防彎器分析的兩個重要因素.對管纜和防彎器的力學性能模擬主要為彎曲剛度的模擬，在有限元中輸入實參數的梁單元能夠很好地模擬這一特性.梁單元節點和自由度等比實體單元少，計算效率高，因此根據防彎器特點建立以梁單元為基礎的數值模型.根據管纜軸線上各點處的剛度值建立此處的梁單元，然后對模型施加相應大小和角度的拉力，通過幾何非線性計算得到最終的變形狀態.

以圖4所示二段防彎器為例建立模型，見圖6.

圖 6防彎器基本幾何參數

Fig.6Basic geometric parameters of bending stiffener

由于防彎器和管纜的軸對稱性，當端部拉力T與管纜頂部軸線呈α角度且在任意周向角度下其變形相同，因此某一周向角度下的平面分析狀態可代替其三維情況下的力學特性.建立數值模型時分析防彎器和管纜平面變形的力學特性，并用ANSYS中的BEAM3單元對其任意微段進行模擬.

防彎器的一般材料為聚氨酯，其力學本構關系具有非線性特點.［17］本文采用一種調整防彎器參數進行快速計算的方法，找到防彎器各參數設計規律，并在某一工況條件下得到防彎器近似的合理設計參數，以減少計算時間.根據設計經驗，在防彎器保護下管纜連同防彎器的最大曲率一般小于0.1，由此產生的防彎器最大應變在0.05d1以內.聚氨酯在此應變范圍內的本構關系可近似為線性，不會對管纜參數設計規律產生顯著影響.因此，將防彎器材料本構關系作為線彈性材料處理，彈性模量為E.在距離防彎器頂端距離為lx處的Δl微段，其彎曲剛度為防彎器彎曲剛度EIb與管纜彎曲剛度EIt之和，且EIb=πd4x－d43E64（1）式中：dx為防彎器Δl微段的直徑大小，dx=lx－l1d2－d1l2+d1（2）管纜彎曲剛度EIt可由管纜的彎曲剛度試驗得到.在管纜段由于沒有防彎器的保護，EIb為0，梁單元的截面彎曲剛度僅為EIt.防彎器段BEAM3單元實參數H的賦值可根據式（2）設置，但由于實參數H僅在后處理分析起作用，當只需得到管纜應變時，可將實參數H設成管纜的直徑.對某算例進行計算，得到管纜參數和極端外部載荷，見表1.計算時選取2組防彎器幾何尺寸和材料，見表2.

表 1防彎器算例基本信息

Tab.1Basic information of bending stiffener example管纜截面

彎曲剛度/

(N?m2)管纜截面

拉伸剛度/

N管纜

直徑d3/

m最大

拉力

Fmax/kN最大

轉角

α/(°)2 7001090.120010

表 2防彎器計算參數

Tab.2Calculation parameters of bending stiffener參數ls/ml1/md1/md2/mE/MPal/m材料130.50.60.14206材料230.50.80.14306

按表1中載荷及表2中參數1信息建立數值模型，見圖7(a).計算得到管纜的彎矩沿軸向分布情況見圖7(b),可知彎矩主要分布在防彎器固定端（圖7中左端）附近，并沿載荷端（圖7中右端）方向趨近于零.在防彎器保護下管纜的變形及等效應變分布見圖7(c).

(a)數值計算模型

(b)彎矩沿軸向分布，N?m

(c)管纜等效應變

圖 7模型及部分計算結果

Fig.7Model and some calculation results

設管纜軸線上每處截面的最大應變為ε，曲率為k，則k=2εd3（3）ε與k為正比關系,因此可由管纜應變分布得到在管纜軸線方向上的曲率分布情況，表2中兩組參數防彎器的曲率計算結果見圖8.

圖 8管纜軸線方向上的曲率分布

Fig.8Curvature distribution in axial direction of pipe/cable

由圖8可知，按參數2計算比按參數1計算得到的管纜最大曲率小70%以上，表明調整防彎器的幾何尺寸和材料屬性可有效改變彎矩作用下曲率的分布規律，為進一步采用結構優化理論和方法開展防彎器結構優化設計奠定基礎.在防彎器參數設計時通常需要進行多次計算才能得到類似參數2的較優設計.采用簡潔有效的有限元模型可以盡量減少計算資源消耗，對防彎器結構的基本設計起到重要參考.

4結束語

防彎器是保護海洋柔性管纜的關鍵附屬構件，其結構設計和分析往往需要對各項參數多次調整計算，且每一組參數計算時需要解決包括幾何及材料等非線性計算的困難.為提高計算效率，滿足工程設計大量計算的需要，提出包括防彎器以及管纜結構在內的防彎器結構的簡化分析模型.采用梁單元建立防彎器等效平面模型并進行數值分析，模擬防彎器和管纜在拉、彎組合載荷作用下彎矩與曲率分布的關系，給出防彎器和管纜在位分析的高效計算方法.參考文獻：

［1］席勇輝. 海洋柔性管纜疲勞半物理仿真試驗研究［D］. 大連: 大連理工大學, 2012.

［2］de SOUZA J R, ROBERTO R J. Bending stiffeners: a parametric structural analysis［C］// Proc ASME 27th Int Conf Offshore Mech & Arctic Eng, Estoril, 2008: 165174.

［3］API Recommended Practice 17B. Recommended practice for flexible pipe［M］. Washington DC: API Publications, 2008: 1126.

［4］TONG Dongjin, LOW Yingmin. A simple mathematical formulation for bend stiffener analysis［J］. Civil Eng Res, 2010(23): 9496.

［5］VAZ M A, de LEMOS C A D. Geometrical and material nonlinear formulation for bend stiffeners［C］// Proc ASME 23rd Int Conf Offshore Mech & Arctic Eng, British Columbia, 2004: 201206.

［6］薛鴻祥, 唐文勇, 張圣坤. 非均勻來流下深海立管渦激振動響應研究［J］. 振動與沖擊, 2007, 26(12): 1316.

XUE Hongxiang, TANG Wenyong, ZHANG Shengkun. Research on shock excitation response of standpipe in deep sea under heterogeneous incoming flow［J］. J Vibration &d Shock, 2007, 26(12): 1316.

［7］姚宗, 陳剛, 楊建民, 等. 流速分層流場中細長柔性立管渦激振動試驗研究［J］. 上海交通大學學報, 2009, 43(8): 94100.

YAO Zong, CHEN Gang, YANG Jianmin, et al. Test research on shock excitation of flexibility long and thin standpipe in flow field of velocity stratification［J］. J Shanghai Jiaotong Univ, 2009, 43(8): 94100.

［8］王一飛, 潘志遠, 黃小平, 等. 深海立管渦激振動疲勞損傷影響因素分析［J］. 船舶力學, 2006, 10(5): 7683.

WANG Yifei, PAN Zhiyuan, HUANG Xiaoping, et al. Effect of factors on the VIVinduced fatigue damage in deep sea risers［J］. J Ship Mech, 2006, 10(5): 7683.

［9］孫麗萍, 周佳. Abaqus二次開發在海洋柔性立管設計分析中的應用［J］. 船舶工程, 2011, 33(6): 8891.

SUN Liping, ZHOU Jia. Application of Abaqus seconddeveloped on design of flexible risers［J］. Ship Eng, 2011, 33(6): 8891.

［10］王愛軍, 楊和振. 深海臍帶纜彎曲加強器時域疲勞敏感性分析［J］. 上海交通大學學報, 2012, 46(10): 16371641.

WANG Aijun, YANG Hezhen. Sensitivity analyses for timedomain fatigue life of umbilical bending stiffener［J］. J Shanghai Jiaotong Univ, 2012, 46(10): 16371641.

［11］王愛軍. 深海臍帶纜動力響應及彎曲加強器可靠性研究［D］. 上海: 上海交通大學, 2012.

［12］CAIRE M, VAZ M A. Viscoelastic analysis of bend stiffeners［C］// Proc ASME 2005 24th Int Conf Offshore Mech & Arctic Eng, OMAE 200567321. Halkidiki, 2005.

［13］LANE M, MCNAMARA J F. Bend stiffeners for flexible risers［C］// Proc Offshore Technol Conf 1995. Houston, 1995.

［14］de MARCO MENICONI L C, LOPES T A P. Fatigue analysis of bend stiffeners［C］// Proc ASME 20th Int Conf Offshore Mech & Arctic Eng, OMAE 2001/OFT1215. Rio de Janeiro, 2001.

［15］VAZ M A, de LEMOS C A D. Geometrical and material nonlinear formulation for bend stiffeners［C］// Proc ASME 2004 23rd Int Conf Offshore Mech & Arctic Eng, OMAE200451366. British Columbia, 2004: 201206.

［16］楊志勛. 臍帶纜結構設計與分析方法研究［D］. 大連: 大連理工大學, 2012.

［17］BOEF W J C, OUT J M M.