屈博志 李玉田
(中海油天津分公司LD27-2/32-2项目组,天津 300452)
LD27-2生产储罐平台PSP为集储罐、工艺、电站和火炬的一体化平台。根据旅大平台和储罐所处的环境工况,本文采用ANSYS软件中的固液耦合的性能对该结构整体进行力学有限元计算,从而得到结构在所处的环境工况下的力学响应,为项目决策提供理论计算依据。
1.1 瞬态动力学分析
瞬态动力学方法又被成为时间历程方法,是一种确定结构在时间相关的载荷作用下的动态响应的技术方法。采用这种分析方法可以确定随时间变化的位移、应变、应力,以及结构在任何静态、动态、谐相应载荷的联合作用下的响应。
瞬态动力学分析的基本方程式为:
(M){a}+(C){v}+(K){u}={F(t)}
其中:
(M)=质量矩阵
(C)=阻尼矩阵
(K)=刚度矩阵
{a}=节点加速度矢量
{v}=节点速度矢量
{u}=节点位移矢量
{F(t)}=载荷矢量
在任何给定的时间t,这些方程可以认为是一组考虑惯性力(M){a}和阻力(C){v}的准静态平衡方程。ANSYS程序使用Newmark时间积分方法或者一种称为HHT的改进方法,在离散的时间点上来求解这些方程。在连续的时间点之间的时间增量称作积分时间步长。
1.2 液体载荷的计算依据
储罐中的原油在受到载荷作用时,会产生晃动,这种晃动会对储槽壁面和底面以及固体支撑结构的受力状况造成一定的影响。
ANSYS软件中专门有一种单元FLUID80适用于计算处理液体晃动问题。该流体单元由8个节点构成,每个节点具有3个自由度。FLUID80流体单元主要用来模拟没有净流率、盛装在容器中的流体。这种单元还特别适用于计算流体固体相互作用的静水压强问题,并且可以考虑加速度效应,特别适合于处理流体晃动问题,同时还可以考虑温度对粘度和阻尼的影响作用。
1.3 固液耦合的计算依据
对于旅大平台上的原油储罐,其内部流体所受的载荷是由平台结构传导而来:平台在各种载荷的作用下产生随时间变化的各个方向的位移,这些位移可以分解为3个自由度的分量。
在计算液体晃动力学问题的过程中,一方面,要把每个时间点上壳体单元各个节点的3自由度位移,实时地传递给储槽中液体边界处各个节点的3自由度位移,就可以计算在外载荷作用下液体的晃动位移、速度等物理参量;而另一方面,获得了液体单元在外载荷下的响应之后,也将液体的位移传递到外围的壳体节点之上,壳体单元则可以通过该位移和材料模量、泊松比等参数,得到壳体自身的应力响应。图2所示为固液耦合作用计算流程图。
LD27-2PSP平台储罐外壁为波纹板结构,有限元建模过程选用SHELL63单元模拟储罐外壁;PIPE59模拟储罐支撑构件;FLIUD80模拟储罐内储液。
3.1 粘度选择
由于储罐内的原油粘度不是一个恒定值,为了论证不同粘度水平对波高的影响,本报告选取了水平振幅40cm,垂直振幅1cm,对不同粘度的原油进行了模拟。结果如下表1所示,没有显著区别。
3.2 液体晃动冲击压强分析
通过选取了水平振幅50cm,垂直振幅1cm,对满罐工况条件下不同周期激励源激励状态下的罐内液体运动模式进行了有限元瞬态动力学分析,取得了罐壁冲击压强变化指数。
?
1)x=50cm,y=1cm 激励周期4.5s
2)x=50cm,y=1cm 激励周期9s
根据晃动冲击压强数值显示,最大值 Pmax=1477.67Pa;无自重的液体晃动所产生的对罐侧壁的液动冲击压强,占储罐整体结构应力的极小部分。
PSP平台设置了13000方的原油储罐,平台在外输油轮靠泊以及地震外力作用下会出现周期性的整体位移,当平台出现位移时,储罐内液体即发生整体波动,本文介绍了储罐内储液固液耦合的计算方法和步骤。
4.1 从流体动力学原理出发,计算结果表明储罐原油在振动工况作用下潜在2种情况,一是液体表面出现跳动性激荡;二是液体整体波动现象。
4.2 储液晃动分析可以模拟分析平台整体在一定位移条件下,储罐内整体液体晃动的情况。
[1]美国船级社,海上固定平台规划、设计和建造的推荐作法--工作应力设计法,2000
[2]胡军,李长升;影响浅海平台自振特性的敏感性因素分析,中国海洋平台,2005.06
[3]杨树耕,孟昭瑛等.有限元分析软件ANSYS在海洋工程中的应用,中国海洋平台,第15卷,第 2期,P41~P44
[4]中国船级社,《海上固定平台入级与建造规范》(1992).
[5]ANSYS,Inc.Release8.0 Documentation for ANSYS.
[6]刘涛,杨凤鹏.《精通ANSYS》清华大学出版社,2002
[7]胡军,张阳.坐底式平台地震响应谱分析[J].中国海洋平台,2008,12.