大型外浮顶立式储液罐地震反应分析.docxVIP

大型外浮顶立式储液罐地震反应分析 石油和天然气公司的主要储油罐大多采用外浮屋顶的形状，安全合理地设计大型垂直储液罐的抗弯设计。这一直是国内外科学家关注的问题。到目前为止，国内外对大型外浮屋顶储油罐的地震资料进行了深入研究，并对其进行了充分的理论研究、计算和分析，特别是由于液固耦合振动的困难和计算手段的缺乏，大部分研究都是对开瓶液罐进行的，忽视了对开瓶液体的影响。事实上，外浮屋顶可以减少液体的流动。随着计算机和金元软件的发展，利用大型通用开元分析软件对大型储液罐的地震反应进行了分析。考虑到外部浮屋顶的影响，本文首先建立了一种适用于该模型的垂直储液罐元音模型。该模型的可靠性得到了改进。在分析不同形状的浮顶（平浮顶、正锥形和反锥形）下的大型储存罐和开放式气罐的地震反应时，不同形式的浮顶对存储罐的地震反应的影响，并确定最合适的浮顶形式。 1 罐体模型的建立 根据薄壁罐的特点,选择具有四个节点,每个节点有六个自由度的SHELL181单元建立罐壁模型,在动力荷载下,允许单元的横向剪切变形,允许单元伸展,因而厚度设置为固定的单元参数;罐内液体采用基于Housner模型简化的三维液体单元FLUID80,此单元适合于容器液体、无净流速的液体模型分析,以及有关静水压力计算、固-液耦合计算和加速度效应等;基础采用SOLID45单元,该单元具有塑性、蠕变、膨胀、应力刚化、大变形、大应变的特性;浮顶由于其厚度相对较小,也采用SOLID45实体单元模拟;灌顶抗风圈采用BEAM4单元,该单元可用于计算应力硬化及大变形问题.此模型建立的关键是应用接触单元和目标单元来模拟移动边界的非线性,以此实现了液体的大幅晃动、罐底提离和浮顶与液面随动,采用的单元是与以上罐体单元相匹配的CONTA173单元和TARGE170单元.状态非线性分析的关键在于罐体与液体接触面之间设置接触对,罐底与基础之间设置接触对,且在液体上表面与浮顶下表面设置接触对. 所建的模型以容积为10000 m3外浮顶储液罐为原型,ANSYS模型罐尺寸参考实际罐规格.罐体直径为30.4 m,罐高14.27 m,罐内液体高度12 m.罐壁采取变厚度,这是由于储罐的罐壁承受储液的静压,此静压是按三角形分布,由上至下逐渐增大,故罐壁厚度由上至下逐渐增厚.0～1.5 m为18 mm,1.5～3.0 m为17 mm,3.0 ～4.5 m为15 mm,4.5 ～6.0 m为13 mm,6.0～7.5 m为11 mm,7.5～9.0m为10 mm,9.0 ～10.5 m为8 mm,10.5 m-14.3 m为7 mm,罐底厚度为8 mm.浮顶浮船外径比储罐内径小500 mm,即为29.9 m,平浮顶厚度取600 mm. 罐壁材料采用490 MPa级的高强度钢材,考虑材料的非线性,采用双线性随动强化模型,弹性模量为2.06×1011N/m2,屈服强度为4.9×108N/m2,切线模量为2.06×109N/m2.浮顶弹性模量为2.06×1011N/m2,密度近似采用实际重量与体积之比为1×102kg/m3.基础高度为1 m,直径为60 m,考虑基础的压缩性,采用弹性模量为1.8×109N/m2.罐内液体采用原油,密度为(0.7～1)×103kg/m3,取1×103kg/m3,弹性模量3×109N/m2,粘滞系数为0.00113 Ns/m. 整个模型网格划分采用映射网格划分,在网格划分时力求罐壁网格与罐内液体网格划分结点一致,罐底结点与基础结点一致,浮顶结点与罐内液体结点一致.外覆顶立式储液罐有限元模型见图1,接触单元模型见图2. 2 垂直储液罐地震反应分析 地震动采用水平单向EL CENTRO波,加速度最大峰值为0.1g,出现在2.14 s,持续时间12 s,时间步长t为0.02 s. 2.1 罐底液体支撑应力场 由储液罐震害可知,储液罐震害大部分是罐壁底部损坏;一部分储液罐发生罐顶破坏,而且明显地看出储液溢出,这都说明了地震时储液晃动.计算所得与罐壁毗邻左右结点的液体高度变化时程曲线见图3, 从图中可以看出,液面晃动最大高度为17.3 cm,液面晃动是以长周期波动为主.计算所得罐底侧壁左右结点的竖向压应力时程曲线见图4,罐底左右结点竖向压应力峰值为11.8 MPa.计算所得罐底左右结点的竖向位移时程曲线见图5,罐底左右结点竖向位移峰值为5.6 mm.显然罐底产生了提离. 2.2 浮顶对储液摆动的影响 常规采用的平浮顶是指浮顶底面是平面的浮顶,计算所得平浮顶储液罐与罐壁毗邻左右结点的液体高度变化时程曲线见图6,从图中可以看出,液面晃动最大高度为16.1cm,与敞口罐相比液面晃动幅度减小了1.2 cm,显然,平浮顶起到了减小液面晃动的作用.液面晃动仍是以长周期波动为主.计算所得罐底左右结点的竖向压应力时程曲线见图7,可以