基于 SESAM 软件的张力腿平台疲劳强度分析

嵇春艳，郭建廷，吴　帅，陆　韵（江苏科技大学 船舶与海洋工程学院，江苏 镇江 212003）

基于 SESAM 软件的张力腿平台疲劳强度分析

嵇春艳，郭建廷，吴帅，陆韵
（江苏科技大学 船舶与海洋工程学院，江苏 镇江 212003）

以张力腿平台为研究对象，建立有限元模型，通过对结构进行总体强度分析，确定疲劳关键区域和关键点。然后建立关键区域的局部模型，对局部模型的网格划分方法进行研究，利用 SESAM 软件中的 Sestra 模块计算应力传递函数。最后基于 Stofat 模块，采用线性外推法得到疲劳关键点处的热点应力，根据 S-N 曲线，计算出结构的疲劳寿命。结果表明，张力腿平台的疲劳强度满足设计寿命要求。

张力腿平台；疲劳强度；有限元

0　引　言

随着陆地石油资源的不断减少，人们把目光转移到海洋。近年来，海洋石油的开发和利用越来越频繁，而且逐渐从浅水向深水发展［1］。张力腿平台是一种适用于深水的海洋平台，因此对张力腿平台进行探索和研究成为热点［2］。

张力腿平台与其他平台相比，具有诸多优势［3］。首先，受张力筋腱的影响，当平台处于平衡位置时，垂荡非常小，在风浪等载荷的作用下，移动相对较小，平台安全稳定；其次，张力腿平台在采油时，不需要辅助钻井船，降低了生产成本；而且张力腿平台既适用于几百米的浅水，也适用于上千米的深水甚至是超深水，既适用于小型油田，又适用于大规模的区块油田。但张力腿平台在运营过程中会因为种种原因毁坏，其中疲劳破坏是最常见的原因之一［4］。而结构的疲劳破坏会对海洋石油的开采、工作人员的安全及海洋环境都会产生影响。因此，为保证张力腿平台的安全性，并使其经济效益最大化，对张力腿平台进行疲劳寿命预报，疲劳强度分析与校核就显得尤为重要［5-7］。

目前，对于张力腿平台，研究总体强度和耦合运动响应的比较多［8］；而对于张力腿平台的疲劳问题，其研究对象大多是系泊系统［9］，对张力腿平台主体结构的疲劳强度研究较少。本文利用 SESAM 软件，对张力腿平台主体结构进行疲劳强度分析，并在分析过程中，研究关键节点有限元模型的划分方法，总结网格划分过程中的经验。

1　疲劳分析的谱分析法

谱分析法是海洋平台疲劳强度分析过程中比较常用的分析方法，中心思想是首先求解出海洋平台受到的波浪载荷，然后通过解析方法得到结构的波浪力谱，最后施加到海洋平台上，得到结构应力响应谱［10］；也可通过数值分析方法，直接获得海洋平台的应力响应谱，然后结合结构响应谱的统计特性、应力响应的分布特性及 Miner 累积损伤理论来预报海洋平台的疲劳寿命。采用谱分析法计算结构疲劳寿命的具体过程如下：

1）计算应力传递函数和应力响应谱

假设波浪是一个平稳的随机过程，经过变换得到的交变应力也是一个平稳的随机过程，由随机过程理论，上述 2 个平稳随机过程的功率谱密度之间有下列关系：

2）计算应力响应谱的 n 阶谱矩

式中：mn为应力响应谱的 n 阶矩；σ 为短期瑞利分布的均方差；TZ为波浪平均跨零周期。

3）求应力响应短期分布概率密度函数

假定各个短期海况应力响应分布若符合 Rayleigh分布，概率密度函数形式如下：

4）计算某一短期海况的疲劳损伤

波浪载荷是随机载荷，其在结构内引起的交变应力也是一随机过程，因此结构在某浪向、某短期海况作用下的疲劳累积损伤度可用下式表示：

式中：NL为在所考虑的整个设计期内应力范围的总循环次数；A和m 为 S-N 曲线的参数。

5）计算总的疲劳损伤值

考虑各个浪向的概率以及该浪向下的各短期海况，可以计算得到总的疲劳累积损伤为：

式中：n 为短期海况的个数；k 为浪向的个数；ps， i为某一海况 i 的发生概率；pd， j为某一浪向 j 的分布概率；dij为在某一海况、某一浪向下的疲劳损伤度。

2　张力腿平台疲劳寿命分析

2.1有限元模型的建立

张力腿平台主要由下浮体、立柱及上部甲板组成，采用 SESAM 软件的 GeniE 模块建立有限元模型，如图1 所示，平台主尺度信息如表1 所示。

在建立整体有限元模型的过程中，采用 2 m 的网格密度，由板单元、梁单元和质量块组成，共有 34 588个节点和 48 058 个单元。

图1　有限元模型Fig. 1　Finite element model

表1　主尺度信息Tab. 1　Information of main scale

2.2疲劳关键节点的选取

研究表明，疲劳破坏往往发生在结构的不连续处和应力集中的焊接节点处，对张力腿平台所有可能发生疲劳破坏的区域进行疲劳强度分析的工作量大，可行性差。比较合理的做法是：根据总体强度分析结果来筛选关键区域，在此过程中，要重点关注对总体强度影响较大的不连续处的焊接节点［11］。

本文通过 SESAM 的 Wadam 模块对张力腿平台进行水动力分析，得到一阶波浪力、附加质量和附加阻尼等水动力基本参数，然后通过 SESAM 的 Sestra 模块对张力腿平台进行总体强度分析，总体强度分析结果如图2 所示，图3 和图4 分别为张力腿平台所受的静水压力示意图和波浪载荷示意图。

从图2 可看出，此张力腿平台整体应力分布均匀合理，在立柱与下浮体连接处和立柱与上甲板连接处的应力比较大，图5 和图6 分别为这两处的局部放大图，而且这两处恰恰是结构不连续处的焊接节点。因此疲劳关键区域选为立柱与下浮体连接处和立柱与上甲板连接处。确定疲劳关键区域后，疲劳关键节点的坐标也随之确定，如表2 所示。

图2　总体强度分析Fig. 2　Analysis of overall strength

图3　静水压力示意图Fig. 3　Diagram of hydrostatic pressure

图4　波浪载荷示意图Fig. 4　Diagram of wave load

图5　立柱与下浮体连接处Fig. 5　Connection of column and pontoon

图6　立柱与甲板连接处Fig. 6　Connection of column and deck

2.3关键节点有限元模型的划分方法和疲劳载荷计算

如图2 所示，把立柱与下浮体连接处定为关键节点 1，立柱与上甲板连接处定为关键节点 2，通过GeniE 模块建立关键节点的有限元模型，网格划分后的模型如图7 和图9 所示。

根据局部模型边界节点的坐标，在整体模型中取出相应位置节点的位移，作为局部模型的边界条件，利用 SESAM 的 Submod 模块，将总体强度模型中节点的状态传递给子模型。因此，在局部有限元模型网格划分时，特别注意边界处网格单元的大小和位置应与整体模型完全一致。

利用 SESAM 软件对网格进行划分的过程中，如果处理的不恰当，在不同网格密度的过渡区域会出现畸形网格，这样将无法进行疲劳动载荷的计算，也就无法对热点进行疲劳分析，因此要花费很多时间来调整网格密度的大小和不同网格密度板的大小。有时对于个别出现畸形网格的地方，通过插入 Feature edge 来调整局部网格。

表2　疲劳关键节点坐标Tab. 2　Coordinate of fatigue critical node

在网格划分过程中，网格要实现平稳过渡，且要求疲劳热点部位的网格大小为 t×t，其中 t 为疲劳关键节点处的板厚。在实际划分网格的过程中，用的网格密度为 0.03 m，0.1 m，0.5 m，1 m 和 2 m，从里向外逐渐增大，热点区域的网格密度约为板厚 0.03 m，这样算出来的热点应力更为准确，图8 为立柱与下浮体连接处密度为 0.03 m 的有限元模型，图10 为立柱与甲板连接处密度为 0.03 的有限元模型。边界区域的网格密度为 2 m，这是因为总体模型中的网格密度为 2 m，在获取子模型的边界条件时，边界上的节点可以吻合。

在 Wadam 中进行水动力分析时，环境参数为：浪向角为 180°～270°，间隔 10°，等概率分布；周期为 2～32 s，间隔 2 s；海浪谱采用 Pierson Moskowitz（P-M）谱，根据中国南海海浪长期资料，可得到此海域波浪参数的统计值如表3 所示。

2.4疲劳寿命预报

通过 Wadam 进行水动力分析，然后 Sestra 计算关键节点的应力传递函数，最后导入 Stofat 进行疲劳寿命预报。

Stofat 在计算热点应力时，根据 ABS 规范［12］采用外推法。S-N 曲线的选取，同样根据 ABS 规范。对于非管节点结构，在计算结构的疲劳寿命时，采用 E 曲线；关键节点 1（立柱与下浮体连接处）在海水中，所以采用 ABS-E-CP 曲线；关键节点 2（立柱与上甲板连接处）暴露在空气中，所以采用 ABS-E-A 曲线。它们的参数如表4 所示。

本文张力腿平台的疲劳设计寿命为 20 年，根据ABS 规范［12］，疲劳安全因子应分别取为 3，最后通过Stofat 计算出关键节点处的疲劳寿命如表5 所示。

表3　南海海况统计资料Tab. 3　Statistics of the South China Sea

表4　S-N 曲线参数Tab. 4　Parameter of S-N curve

表5　关键节点疲劳寿命结果Tab. 5　Fatigue life of critical node

3　结　语

海洋平台结构疲劳强度评估是平台结构安全评估的重要内容之一。本文应用 SESAM 软件，对一张力腿平台关键部位进行了疲劳强度研究，其主要结论如下：

1）通过对张力腿平台的总体强度分析，获得了平台容易产生疲劳破坏的 2 个关键节点，分别为立柱与下浮体连接处和立柱与甲板连接处。

2）利用 SESAM 软件对结构的疲劳强度进行分析的过程中，特别要注意关键节点有限元模型的网格划分，建议采用多梯度平稳过渡的方法，这样可以避免畸形网格的产生，从而使疲劳关键节点的应力更加准确。应该在局部模型的边界处设置与整体模型中相同的网格密度，并且局部模型的边界应在舱壁、T型材和加强筋处，这样与整体模型中的网格能够吻合，从而使局部模型获得更加精确的边界条件。

3）对关键节点进行疲劳寿命分析结果表明，立柱与下浮体连接处的疲劳寿命为 62.7 年；立柱与甲板连接处的疲劳寿命为 148 年。本平台的设计寿命是 20年，安全因子为 3，疲劳寿命预报结果表明关键节点满足设计要求。相比较而言，立柱与下浮体连接处的疲劳损伤度大，疲劳寿命短，而且该位置处于海水中，如果发生疲劳破坏不易被发现，因此应进行定期的检查和维护。

［1］李润培， 谢永和， 舒志. 深海平台技术的研究现状与发展趋势［J］. 中国海洋平台， 2003， 18（3）： 1-5.

［7］杨鹏， 顾学康. 半潜平台结构疲劳寿命评估方法比较［J］. 舰船科学技术， 2012， 34（8）： 112-118. YANG Peng， GU Xue-kang. Comparative researches of structural fatigue life assessment procedures for a semi-submersible platform［J］. Ship Science and Technology， 2012， 34（8）：112-118.

［8］谷家扬， 吕海宁， 杨建民. 张力腿平台在随机波浪中的耦合运动响应研究［J］. 船舶力学， 2013， 17（8）： 888-900. GU Jia-yang， LV Hai-ning， YANG Jian-min. Studies on coupling dynamic response of TLP in stochastic waves［J］. Journal of Ship Mechanics， 2013， 17（8）： 888-900.

［9］余建星， 张中华. 波、流联合作用下张力腿平台系索系统疲劳可靠性分析方法研究［J］. 海洋技术， 2009， 28（2）： 69-72. YU Jian-xing， ZHANG Zhong-hua. Research on fatigue reliability of tension leg platform tendon system under wave and current［J］. Ocean Technology， 2009， 28（2）： 69-72.

［10］梁瑜. 基于谱分析法的半潜平台疲劳强度研究［J］. 船海工程，2012， 41（1）： 77-80. LIANG Yu. Fatigue strength assessment for semi-submerged platform based on the spectral-based analysis［J］. Ship & ocean Engineering， 2012， 41（1）： 77-80.

［11］白艳彬. 深水半潜式钻井平台总体强度分析及疲劳强度评估［D］. 上海： 上海交通大学， 2010.

［12］ABS. Guide for the fatigue assessment of offshore structures［S］. Houston， TX： American Bureau of Shipping，2003.

Fatigue strength analysis of tension leg platform based on SESAM

JI Chun-yan， GUO Jian-ting， WU Shuai， LU Yun
（School of Naval Architecture and Offshore Engineering， Jiangsu University of Science and Technology， Zhenjiang 212003， China）

The finite element model was made taking the tension leg platform as the research object. The key area of the fatigue and the key points' coordinates of the fatigue were determined through overall strength analysis. Then the local model of the key area was made. The approach of dividing grids was researched about the local model. The stress transfer function was calculated through using Sestra module in SESAM. In the end， based on the Stofat module， using the linear extrapolation， the stress of the key fatigue points was obtained. According to the S-N curve， the fatigue life of tension leg platform was calculated. The results show that the fatigue strength of tension leg platform meets the requirements of the design life.

tension leg platform；fatigue strength；finite element

P751

A

1672-7619（2016）05-0039-05

10.3404/j.issn.1672-7619.2016.05.009

2015-10-27；

2015-11-30

国家自然科学青年基金资助项目（51409129）

嵇春艳（1976-），女，教授，研究方向为船舶与海洋结构物的力学性能、动力性能及减振方法。