导管架非节点有限元简化疲劳分析方法
包清华,冯 珩
(中国石油集团海洋工程有限公司工程设计院,北京100176)
摘要:文章以一座已投入使用的导管架平台裙桩套筒的设计为例 详细介绍了复杂非管节点裙装套筒结构的简化疲劳分析方法:首先采用SACS程序对导管架平台在波浪荷载作用下进行结构整体分析,选取简化有限元疲劳分析中所需的荷载工况组、荷载和边界条件,然后应用ALGOR程序建立裙桩套筒的有限元模型, 施加荷载和边界条件,并对应力分析输出结果进行后处理得到疲劳分析结果。这种简化疲劳分析方法及其分析结果得到了挪威船级社(DNV)的认可,对导管架平台非管节点,特别是对适度水深和深水导管架的裙桩套筒的疲劳设计具有参考意义。
关键词:海上导管架;裙桩套筒;有限元;疲劳评估
文章编号:l671-9433(201l)03-032l-04
1 引言
海上导管架在大量波浪荷载往复作用下,即使所产生的名义应力较小,也会引起疲劳损伤。对于导管架典型节点,如T、Y、K、TK,双面T和X型等节点,相关规范给出了用于计算节点焊缝的应力集中系数的计算公式,应用海洋工程结构专业分析软件如SACS对这些典型节点的应力集中系数进行设定,从而可以通过软件计算出典型管节点的应力,进而可以直接采用SACS等海洋工程结构专业软件对其进行疲劳评估。因此可以看出,对于典型管节点疲劳特性的评估,通常采用软件直接进行计算分析;而对于复杂非管节点的热点应力计算 这些应力集中系数计算公式不再适用,SACS软件也无法对其进行直接的疲劳特性评估。鉴于此,本文对简化疲劳有限元计算方法进行了探讨,介绍了一种用于海上导管架平台复杂非管节点结构在波浪荷载作用下的有限元简化疲劳分析方法。并以一座已投入使用的导管架平台裙桩套筒的设计为例,详细介绍了有限元简化疲劳分析的方法和过程。这种有限元简化疲劳分析方法及其分析结果得到了挪威船级社(DNV)的认可,对海上导管架平台非管节点的疲劳设计具有参考意义。
2 分析过程
有限元简化疲劳分析方法包括如下分析过程:首先,根据导管架设计图建立如图1所示的导管架SACS整体分析模型,对其进行整体分析,从而得到用于简化疲劳分析的危险工况和位移边界条件;然后采用ALGOR有限元分析软件建立导管架最复杂部分裙桩套筒的有限元分析模型,如图2所示,并对裙桩套筒有限元分析模型施加荷载和位移边界条件,之后进行有限元分析;最后根据DNV疲劳规范进行导管架裙桩套筒的疲劳特性评估。
2.1模型
导管架裙桩套筒作为导管结构的重要组成部分,起着将导管架及上部组块荷载传递到桩上的作用,模型包括撑杆、导管架腿、YOKE板、剪切板和带有灌浆的桩和套(见图2)。在模型中,采用了4节点壳单元对撑杆、导管架腿、YOKE板、剪切板进行了模拟。管节点处有限元网格的尺寸接近于杆件的厚度,并且经过多次网格划分,直到每次网格划分后的计算分析结果只有略微变化为止。
图1 导管架的SACS整体分析模型 图2 导管架裙桩套筒有限元模型
裙桩套筒的灌浆连接部分采用8节点实体单元进行模拟,桩和套筒采用4节点壳单元进行模拟,灌浆部分的实体单元与桩、套筒之间的边界条件采取耦合模拟,这表明灌浆实体单元与桩套筒壳单元的有限元网格连接节点将会有一致的位移边界条件;导管架撑杆、导管架腿、桩的端部中心点被模拟成一个刚性节点,端部荷载及位移边界条件被施加在这个点上【4】。
2.2 材料属性
裙桩套筒各个构件组成部分及灌浆连接部分分别为高强度钢材和高标号水泥灌浆,其材料属性如表1所示。
表1 材料属性
项目 |
钢材 |
水泥浆 |
弹性模量/MPa |
2.1times;105 |
2.4times;104 |
泊松比 |
0.3 |
0.167 |
密度(kg/msup3;) |
7.85times;104 |
3.45times;103 |
屈服强度/MPa |
345 |
3 疲劳强度分析
3.1 加载
导管架裙桩套筒整体分析模型如图3所示。采用SACS软件对整体分析模型进行多个方向的波浪荷载分析,莫里森方程被用来计算由波浪产生的拖曳力和惯性力。根据整体分析结果选出波浪条件下的危险疲劳工况组,如表2所示,这个疲劳工况组将用来计算裙桩套筒节点处的应力幅。疲劳工况组是根据API规范进行选择的,为具有一百年重现期的波浪荷 载工况,其产生最大/最小基底剪力的相位与产生最大/最小倾覆力矩的波浪荷载工况相同, 因此具有最大/最小基底剪力的荷载疲劳工况组被用于进行疲劳分析。
由节点205、201、207确定的平面节点编号 由节点101、103、125确定的平面节点编号
由节点105、101、205确定的平面节点编号 由节点107、103、207确定的平面节点编号
表2 疲劳工况
SACS分析工况 |
工况描述 |
最大 |
产生最大基底剪力的波浪荷载工况 |
最小 |
产生最小基底剪力的波浪荷载工况 |
从SACS整体分析结果输出的导管架上每个节点的线位移和转角位移分别如表3和表4所示。表中的节点编号所对应的节点如图1和图2所示。在有限元局部疲劳分析中,所有的这些 节点线位移和转角位移均作为边界条件施加到模型相应节点端部。
表3 从SACS分析输出的工况MAX下的节点线位移和转角位移
节点号 |
工况 |
X方向位移/cm |
Y方向位移/cm |
Z方向位移/cm |
绕X方向转角/rad |
绕Y方向转角/rad |
绕Z方向转角/rad |
37 |
MAX |
-2.8631668 |
2.7766593 |
-1.9302737 |
-0.0007072 |
-0.0007468 |
0.0000472 |
124 |
MAX |
-3.6661503 |
2.8477535 |
-1.4024048 |
0.0000899 |
-0.0005899 |
-0.0009466 |
125 |
MAX |
-2.9236839 |
3.6827505 |
-1.4064275 |
-0.0006188 |
0.0001194 |
0.0009817 |
203 |
MAX |
-3.7072046 |
3.7753172 |
-0.9421117 |
-0.0002081 |
-0.0006414 |
-0.0018477 |
表4 从SACS分析输出的工况MIN下的节点线位移和转角位移
节点号 |
工况 |
X方向位移/cm |
Y方向位移/cm |
Z方向位移/cm |
绕X方向转角/rad |
绕Y方向转角/rad |
绕Z方向转角/rad |
37 |
MIN |
1.5454831 |
-1.5635995 |
-0.4188206 |
0.0005342 |
0.0005278 |
-0.0000953 |
124 |
MIN |
2.0688691 |
-1.7463771 |
-0.7680346 |
-0.0002331 |
0.0005006 |
0.0006445 |
125 |
MIN |
1.6793743 |
-2.1677744 |
-0.7651194 |
0.0005241 |
-0.0002314 |
-0.0005687 |
203 |
MIN |
2.0657177 |
-2.1292496 |
-0.9039879 |
0.0004154 |
0.0000269 |
0.0011738 |
3.2疲劳评价方法
根据DNV规范,上述有限元简化疲劳分析方法被用来进行疲劳特性评估。基于疲劳试验的S-N曲线长期应力幅的分布是用来进行疲劳分析的一个基本要求,作为一种替代,威布尔分布被用来表示应力幅的长期分布,采用威布尔应力分布和S-N双曲线,疲劳损伤表达式如下所示:
节点疲劳设计的设计寿命为二倍25年的服役寿命再加上2年的运输疲劳特性,为52年的设计疲劳特性;平均波周期为9.312秒。 根据eta;=1.0,不同节点类型的允许应力幅如表5所示。
表5 允许应力幅
节点类型 |
描述 |
允许应力幅/MPa |
T |
管节点 |
231 |
F |
管与剪切板或YOKE板连接节点 |
182 |
当管节点板材的厚度大于规范规定的参考厚度时,允许应力幅值可以按下式进行计算得到。
3.3 应力幅计算
与单元局部坐标系相应的每个工况的六个应力分量可以用来计算单元的正应力和剪应力, 并根据如下公式可计算出单元的主应力,从而得到节点处单元的最大应力幅。
3.4 疲劳分析结果
根据3.3节所述,可以计算出导管架裙桩套筒在疲劳工况组条件下节点处的最大应力幅, 如表6
全文共4760字,剩余内容已隐藏,支付完成后下载完整资料
英语原文共 4 页,剩余内容已隐藏,支付完成后下载完整资料
资料编号:[13502],资料为PDF文档或Word文档,PDF文档可免费转换为Word
以上是毕业论文外文翻译,课题毕业论文、任务书、文献综述、开题报告、程序设计、图纸设计等资料可联系客服协助查找。