大型油罐应力分析与屈曲稳定性研究

陈志平, 曾明, 余雏麟, 沈建民, 葛颂, 朱劲平

陈志平, 曾明, 余雏麟, 沈建民, 葛颂, 朱劲平. 大型油罐应力分析与屈曲稳定性研究[J]. 油气储运, 2008, 27(12): 7-12, 46. DOI: 10.6047/j.issn.1000-8241.2008.12.003
引用本文: 陈志平, 曾明, 余雏麟, 沈建民, 葛颂, 朱劲平. 大型油罐应力分析与屈曲稳定性研究[J]. 油气储运, 2008, 27(12): 7-12, 46. DOI: 10.6047/j.issn.1000-8241.2008.12.003
CHEN Zhiping, ZENG Ming, . A Research on Stress Analysis and Bulking Reliability of Large Oil Storage Tank[J]. Oil & Gas Storage and Transportation, 2008, 27(12): 7-12, 46. DOI: 10.6047/j.issn.1000-8241.2008.12.003
Citation: CHEN Zhiping, ZENG Ming, . A Research on Stress Analysis and Bulking Reliability of Large Oil Storage Tank[J]. Oil & Gas Storage and Transportation, 2008, 27(12): 7-12, 46. DOI: 10.6047/j.issn.1000-8241.2008.12.003

大型油罐应力分析与屈曲稳定性研究

基金项目: 

国家高技术研究发展计划资助项目 2007AA04Z427

详细信息
    作者简介:

    陈志平, 教授, 1965年生, 2006年博士毕业于浙江大学化工过程机械专业, 现在浙江大学化工机械研究所从事大型石油储罐结构优化与安全性分析、高参数压力容器研发, 新型节能单元设备研制等研究工作

    **310027,浙江省杭州市浙大路38号;电话:(0571)87953007。

A Research on Stress Analysis and Bulking Reliability of Large Oil Storage Tank

  • 摘要: 初步分析了罐壁应力分布的基本特点,提出了一种组合圆柱壳理论的应力计算解析方法。比较了已有罐底板应力计算中罐一土耦合有限元分析的建模方法,建立了以基础沉降量为边界条件的新模型。采用数值模拟优化分析方法,研究了不同结构参数对油罐安全性的影响。利用非线性有限元方法,进行了准静态弹塑性屈曲分析。研究结果表明,新的计算方法和计算模型可大大提高罐体应力计算结果的准确度,部分研究成果可为制订我国大型立式圆筒形钢制焊接油罐标准提供技术数据。
    Abstract: After analyzing the basic features of stress distribution over tank shell, a stress calculation method is proposed based on the theory of combined cylindrical shells. Compared with the existing modeling methods for stress calculation of bottom by using finite element analysis of tank-soil coupling, a new model which takes foundation settlement as the boundary condition is established. The influence o safe reliability from different structural parameters is studied by using analytical methods of numerical simulated optimization. And elastic-plastic buckling research under quasi-static is carried out by using nonlinear finite element analysis.The results show that new methods and models can significantly improve the accuracy of stress calculation, and part of the achieved results can provide technical data for making standards, which are concerned with large stand cylindrical steel welded oil storage tank.
  • 近年来对冷热油交替顺序输送过程的热力问题进行了相关领域的研究, 建立了有关冷热油交替顺序输送过程的数学模型, 其解法中选择的关键参数都是从以往管道投产及正常生产中经过多年验证的数值。在对文献〔1〕中不妥的边界条件和相关公式进行修正后, 得出了初步计算结果, 同时结合温度场蓄热量计算和东营—滨州管道投产测温结果加以进一步论证, 指出了与文献〔1〕存在的差异, 以供商榷。

    文献〔1〕中采用二维热传导偏微分方程, 描述的土壤导热方程为:

    (1)

    x=0时, 且0≤yhC-D0/2或x=xd, hC+D0/2≤yyd时, 边界条件为:

    (2)

    r=D0时, 且Tw=Ts, 边界条件为:

    (3)

    式中  Ts——土壤温度, ℃;

    τ——时间, s;

    αs——土壤导温系数, m2/s;

    hC——管轴埋深, m;

    xd——直角坐标系中x方向管道热力影响区域的距离, m;

    D0——管道的最外层直径, m;

    yd——直角坐标系中y方向管道热力影响区域的距离, m;

    λs——土壤的导热系数, W/m·℃);

    r——管子半径, m;

    αh——当量换热系数;

    Ty——管内油流的平均温度, ℃;

    Tw——管道最外壁面的温度, ℃。

    从根据选择文献〔1〕中的边界条件所画出的求解区域示意图(见图 1)可以看出, 线段OE和BC为当x=0, 且0≤yhC-D0/2或x=xd, hC+D0/2≤yyd两个限定边界表示的线段, 而线段FG和AB则不包括在内, 这显然是不妥的。

    图  1  土壤求解区域示意图

    事实上, 文献〔1〕中式(2)的真正限定边界应为:

    (4)

    边界条件为式(2), 这在形式上就对了。但是, 当x=xd仍存在如下问题。

    (1) 在式(4)限定边界条件下, 给出的式(2)意味着在边界x=xd处直线AC上是绝热边界, 即热量传到xd处就不再向外传播, 如沿着边界线AC存在一块大的绝热材料。此问题引起的后果是, 在稳定运行状态下, 总能使管道与AC线之间达到热平衡, 从管道放出的热量只能向管道上方或下方传播, 而不能向左右传播, 提高了稳定状况下温度场的蓄热量。使得进行冷热油交替输送热力计算前, 计算出的稳定温度场本身就存在较大的失真。

    (2) 文献〔1〕中没有考虑自然地温(T0)随土壤深度的变化。从文献〔2〕中可以看出, 采用二维热传导方程应考虑自然地温随深度的变化, 否则, 计算出的输油管道稳定温度场也不能反映温度场的实际情况。因此, 实际上边界条件式(2)应改为当x=0且0≤yhC-D0/2或hC+D0/2≤yyd时, 边界条件为; 当x=xd时, 边界条件为Ts=T0(y, t)。后者属于第一类边界条件, 即把直线AC以外部分看成无限大的具有温度为T0(y, t)的热源, 而且应明确给出T0(y, t)函数关系。

    在文献〔1〕的边界条件式(3)中, 当量换热系数(αh)除包括原油至管内壁的对流放热系数(α1)外, 还包括管壁、绝热层和保温层等在内的各层热阻。其表达式为:

    (5)

    式中  d——钢管的内径, m;

    λi——第i环层导热系数, W/(m·℃);

    Di——第i环层的外径, m。

    对于稳定运行情况而言, 式(5)为多层合在一起的当量换热系数。在稳定运行中, 虽然各层存在热流, 但各层温度是恒定不变的, 且热流量是连续的, 可以用式(5)总的换热系数来描述热交换过程。但在冷热油交替输送时, 各层间有热量传播且各处温度也在频繁变化, 处于极不稳定的传热状态, 各层间热流量并不完全相等。一个时刻管内油流温度(Ty)不能对应于该时刻管道最外壁的温度(Tw), 因此边界条件式(3)不成立, 特别是对于含有保温层的情况更是如此。一般的处理方法是另外列出各层的热传导微分方程, 进行联立求解3, 4, 除非经过多次实际管道的运行考验, 证明这种近似不带来大的误差, 否则式(3)边界条件是不妥的。

    东营-滨州输油管道改线段56 km, 管道直径为529×8 mm, 防腐层厚度为6 mm。1990年6月15日投产5时管道埋深处自然地温为17.5℃。该管道投产采用热水一次性预热后输油, 即准冷投投产方法。在总水量为8 100 m3, 输水排量为800 m3/h, 热水出站温度为65℃时, 经过现场测量, 滨州进站的水温与时间的关系见表 1

    表  1  滨州站进站水温实测值与时间的关系
    下载: 导出CSV 
    | 显示表格

    表 1可知, 热水头的进站温度只有19.2℃, 经过4 h后, 进站水温升高到34.4℃。

    在文献〔1〕的计算结果中, 俄罗斯原油出站温度为16℃, 输送93.6 h后的油尾进站温度为18.4℃。转换为输送大庆油时, 其油头所感受到的管道内壁温度在16~18.4℃之间, 平均温度为17.2℃, 这个平均温度与东营—滨州管道投产时的地温17.5℃基本相同。如果说铁秦输油管道内壁温度较低(17.2℃)并不代表管外土壤温度, 由于原温度场的存在, 管外土壤温度会稍微高一些。但由于连续近4天输送温度只有16℃的冷油, 管外温度场中靠近管道处土壤的热量要向管道内部不断地传播, 这些热量被冷油吸收。而较远地方土壤的热量还要继续向更远处传播, 土壤温度场, 特别是靠近管道的部分将有明显的耗散。因此, 从实际情况分析, 不可能使大庆油头的进站温度仍然维持在与原来稳态运行时很接近的数值, 即从俄罗斯原油的油尾温度18.4℃快速升到36.2℃(见文献〔1〕中的表 3)。

    从能量比较角度分析, 如果管输原油携带的热能多, 则进站温度会快速升高。只有当管输的液流温度高于环境温度时, 才能使得液流的热量向周围散失, 而单位时间进入管道的热量(P)的数值为:

    (6)

    式中  P——单位时间进入管道的热量, J/h;

    ρ——流体的密度, kg/m3;

    Q——输量, m3/h;

    c——热容量, J/(kg·℃);

    ΔT——出站温度与环境温度之差, ℃。

    文献〔1〕中, 大庆油的输量(Q1)为2 000 m3/h, 油品密度(ρ1)为862 kg/m3, 热容量(c1)为2 094 J/kg·℃), 出站温度与环境温度之差(ΔT1)为25.8℃。而东营-滨州管道输水预热的输量(Q2)为800 m3/h, 密度(ρ2)为1 000kg/m3, 热容量(c2)为4 187 J/(kg·℃), 出站温度与环境温度之差(ΔT2)为47.5℃。

    由以上参数可得, 铁岭—新民站间输送大庆原油的热流量(P1)为93.14 kMJ/h, 东营至滨州间输水时热流量(P2)为159 kMJ/h, 且经过一个站间管段, 前者需18.4 h, 后者只需14.5 h。由于热流量P2远大于P1, 且在管道中行走时间较短, 两者外部环境情况相类似, 因此滨州进站水温及其增长速度比新民站进站油温增长更快才最为合理。

    然而, 测量得到滨州进站水温经过了4 h才上升到34.4℃, 而用文献〔1〕计算的新民进站油温(大庆原油)在油头部分可达到36.2℃, 这一计算结果在现实中是难以实现的。

    修正文献〔1〕中的边界条件后, 加入了描述防腐层(或保温层)的二维热传导微分方程, 与土壤热传导方程和油流方程联立求解, 得到如下初步结果。

    以铁岭—新民站间管段为例计算。该管段直径为720 mm, 壁厚为8 mm, 站间距为92 km, 热油(大庆原油)出站温度为43℃, 稳态运行时大庆油进站温度36℃, 冷油(俄罗斯原油)出站温度为16℃。按输量2 000 m3/h, 拟定7天为一个输送循环周期(冷油4天, 热油3天), 在每个循环中, 先输冷油(俄罗斯原油)4天, 再输热油(大庆原油)3天。自然地温为16℃, 输油方向为从铁岭输往新民, 铁岭站为出站端, 新民站为进站端, 计算得到第一个循环周期的比较结果(见图 2图 3)。

    图  2  第一个循环中的冷油进站温度随时间变化关系
    图  3  第一个循环中的热油进站温度随时间变化关系

    图 2给出了第一个循环中冷油进站温度随时间的变化关系。可以看出, 由于刚经过热油稳态运行后, 温度场温度较高, 因此, 冷油在92km管道中长时间行走后, 受到原先温度场的反加热, 使得油头温度上升较高, 达到30.2℃, 这个温度比稳态运行时热油进站温度(36℃)约低6℃左右, 这是符合客观规律的。随着冷油进站时间的延长, 冷油进站温度逐渐下降, 经过12h下降到27.5℃, 4天(96h)后, 冷油油尾进站温度下降到23.5℃。

    图 3给出了第一个循环中, 冷油输完以后, 接着输热油时, 热油进站温度随时间变化关系。可以看出, 由于冷油已输送了4天, 温度场温度消耗了许多, 因此, 热油油头的进站温度仅为29.7℃, 比稳态运行时热油进站温度低了约6℃, 比冷油油尾进站温度高约6℃, 而且随着运行时间的推移, 热油进站温度逐步上升。因此, 要使热油的进站温度上升到目前输送的安全进站温度33℃, 需要约24 h。经计算得知, 输送72 h后, 热油尾进站温度为35.2℃。

    以上给出的仅是第一个循环周期的结果。不难预测, 随着循环次数的增加, 冷油头及冷油尾的进站温度与第一个循环相比将还要不断下降, 而且热油头的进站温度也同样不断降低, 要使热油的进站温度达到安全运行数值33℃, 进站时间将拖得更长。

    冷油头进站温度30.2℃, 冷油尾进站温度为23.5℃, 由于冷油在铁岭的出站温度为16℃, 输往新民站后, 冷油温度平均升高为10.85℃。周围温度场有没有那么多热量来加热冷油?为直观和明了起见, 对能影响到冷油反加热的那部分土壤温度场所包含的热量作了近似的计算。

    在热油稳态运行中, 出站温度43℃, 进站温度36℃, 据现场测温得知, 对于这种非保温管道, 管外壁温度低于油温大约2℃, 铁岭出站端管外壁温度为41℃, 新民进站端管外壁温度为34℃。由于冷油在管道中仅行走4天, 热影响半径定为距管道中心1.5 m是可行的, 并且取距管中心1.5 m处的土壤温度等于相应管外壁温度的1/2, 即靠近出站端为20.5℃, 靠近新民进站端为17℃。用线性回归方法描述土壤各点温度随轴向(z)和径向(r)的变化, 不难导出。

    (7)

    从而, 土壤蓄热量的微分形式为:

    (8)

    式中  ws——土壤蓄热量, J;

    Cs——土壤比热, J/(kg·℃);

    ρs——土壤密度, kg/m3

    z方向从0~92 km, r方向从管外壁0.36~1.5 m, 做二重积分, 得到土壤蓄热量(ws)为22 148 kMJ, 俄罗斯原油(冷油)密度为840 kg/m3, 比热为2 094 J/(kg·℃), 则冷油输送到新民站后, 总吸收的热量(we)为3 664 kMJ。所以得到土壤蓄热量(ws)与冷油吸收的热量(we)的比值为6.045。这就是说, 冷油吸收的热量(we)只有所能感受到的土壤蓄热量(ws)的1/6。因此, 冷油经过稳态温度场反加热后, 温度平均升高10.85℃是可能的。比较图 3和文献〔1〕的计算结果可以看出, 文献〔1〕的热油进站温度算得过高, 在实际中难以实现。

    (1) 认为文献〔1〕中给出的微分计算边界条件及计算公式存在不妥之处, 其计算结果与以往相关测试有很大差距, 因此, 在对管道实际运行起指导作用前, 还需要经过一定的修正。

    (2) 类似于大庆原油和俄罗斯原油这样流动性相差较大的两种原油, 进行冷热油交替输送时, 其一要进行认真对待, 因为它相当于每循环一次都进行了一次准冷投; 其二是在进行冷热油交替输送时, 若循环周期较短(如7天以内)且冷油总量与热油总量相近, 必须及早大幅提升冷油的出站温度, 保持温度场, 若循环周期长, 且冷油的批次量又大, 则在低温冷油输送一定时间后, 对后续的冷油应大幅度提升出站温度, 重新建立温度场, 确保输送大庆油时的进站温度能够达到安全温度。

  • 图  1   底板与壁板力学分析模型

    图  2   15×104 m3石油储罐壁板外表面环向应力分布

    图  3   充液后罐底板的实际变形示意图

    图  4   基于沉降量法的大型油罐基本模型

    图  5   “象足”屈曲后的非锚固罐有限元网格图

    表  1   不同设计标准对15×104 m3油罐参数和应力程度的影响[15, 16]

    下载: 导出CSV

    表  2   环向应力值一定时系数C的取值

    下载: 导出CSV
  • [1] 武铜柱: 大型立式油罐发展综述, 石油化工设备技术, 2004, 25(3)。 https://www.cnki.com.cn/Article/CJFDTOTAL-SHSJ200403015.htm
    [2] 李国琛: 大型变截面圆柱罐壁和罐底的应力分析, 力学与实践, 1979, 1(4)。 https://www.cnki.com.cn/Article/CJFDTOTAL-LXYS197904011.htm
    [3] 傅强 陈志平等: 弹性基础上大型石油储罐的应力分析, 化工机械, 2002, 29(4)。 https://www.cnki.com.cn/Article/CJFDTOTAL-HGJX200204005.htm
    [4] 陈志平 沈建民等: 基于组合圆柱壳理论的大型油罐应力分析﹐浙江大学学报(工学版), 2006, 40(9)。 https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC200609034.htm
    [5]

    Chen Z P, Duan Y Y, Shen J M, Jiang J LIA simplified method for calculating the stress of large storage tank wall, Proc. IMechE Part E: J. Process Mechanical Engineering, 2007, 221(3).

    [6] 陈志平: 大型非锚固储油罐应力分析与抗震研究(博上论文), 浙江大学, 2006。
    [7] 陈志平 王飞等: 大型非锚固原油储罐应力分析, 机械工程学报, 2006, 42(11)。 https://www.cnki.com.cn/Article/CJFDTOTAL-JXXB200611035.htm
    [8]

    WU Tianyun: More accurate method devised for tank—bottom annular plate design, Oil & Gas, 1996, 94(21).

    [9]

    WU Tianyun, LIU Guorong: Comparison of design methods for tank-bottom annular plate and concrete ring wall, Internation Journal of Pressure Vessels and Piping, 2000, 77(3).

    [10] 李多民 宋虎堂等: 12.5×104m3原油储罐应力测试及分析, 压力容器. 2005, 22(3)。 https://www.cnki.com.cn/Article/CJFDTOTAL-YLRQ200503004.htm
    [11] 沈建民: 大型油罐的静强度及动力响应分析(硕上论文), 浙江大学, 2006。
    [12] 陈志平 葛颂 沈建民等: 大型原油储罐有限元分析建模的新方法, 浙江大学学报(工学版), 2006, 40(6)。 https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC200606011.htm
    [13] 陈志平 程香: 在役人型原油储罐安全评估的新方法, 第十届全国高等学校过程装备与控制工程专业校际交流会论文集, 化学工业出版社(北京), 2006。
    [14] 程香 陈志平等: 基于基础沉降监测的在役大型油罐寿命预测技术, 压力容器, 2007, 24(1)。 https://www.cnki.com.cn/Article/CJFDTOTAL-YLRQ200701007.htm
    [15] 王飞: 大型非锚固原油储罐数值模拟及结构优化(硕士论文). 浙江大学, 2006。
    [16] 段嫒嫒 刘义君等: 15×104 m3超人型石油储罐结构优化分析, 压力容器, 2006, 23(12)。 https://www.cnki.com.cn/Article/CJFDTOTAL-YLRQ200612004.htm
    [17] 陈志平 蒋伟华 沈建民等: 大型油罐大角焊缝处峰值应力分析, 压力容器, 2005, 22(5)。 https://www.cnki.com.cn/Article/CJFDTOTAL-YLRQ200505003.htm
    [18] 庄茁 刘焕忠等: 储液罐动力响应中的能量分析, 第六届全国压力容器学术会议——压力容器先进技术精选集, 机械工业出版社(北京), 2005。 https://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGVE200509002038.htm
    [19]

    Hamdan F H: Seismic behavior of cylindrical steel liquid storage tanks, Journal of Constructional Steel Research, 2000(53).

    [20]

    Haroun, Medhat A: Implications of recent nonlinear analyses on seismic standards of liquid storage tanks, Technical Council on Lifeline Earthquake Engineering Monograph Seattle, WA, USA, American Society of Civil Engineers, Reston, VA, USA, 1999(16).

    [21]

    Priestley M J N et al: Seismic Design of Storage Tanks, Recommendations of a Study Group of the New Zealand National Society for Earth quake Engineering, 1986.

    [22]

    American Water Works Association: Welded Steel Tanks for Water Storage. AWWA Standard D100.1984.

    [23] JIS B8501—1985钢制焊接油罐结构. 日本工业标准委员会。
    [24] 葛颂 陈志平等: 立式储液罐象足屈曲的准静态数值模拟, 压力容器, 2006, 23(1)。 https://www.cnki.com.cn/Article/CJFDTOTAL-YLRQ200601001.htm
    [25] 葛松: 大型立式储液罐抗震分析的数值模拟研究(硕士论文), 浙江大学, 2006。
    [26] 中国石油化工总公司抗震办公室(译): 洛马普里埃塔地震考察, 地震出版社(北京), 1991。
    [27]

    Uwe Hornung, Helmut SaM: Buckling loads of tank shells with imperfection, Non-linear Mechanics, 2002(37).

    [28]

    Htibner A J, Teng G. SaM H: Buckling behaviour of large steel cylinders with patterned welds, International Journal of Pressure Vessels and Piping, 2006(83).

    [29]

    Aghajari S, Abedi K. Showkati H: Buckling and post-buckling behavior of thin-walled cylindrical steel shells with varying thickness subjected to uniform external pressure, Thin-Walled Structures, 2006(44).

    [30] 陈宏湛 沈成康: 组合圆柱壳静态屈曲的几个影响因素分析, 力学季刊, 2001, 22(2)。 https://www.cnki.com.cn/Article/CJFDTOTAL-SHLX200102010.htm
    [31] 朱劲平: 大型立式圆柱形油罐"象足"屈曲研究(硕士论文), 浙江大学, 2007。
图(5)  /  表(2)
计量
  • 文章访问数: 
  • HTML全文浏览量:  0
  • PDF下载量: 
  • 被引次数: 0
出版历程
  • 收稿日期:  2007-11-01
  • 网络出版日期:  2023-08-21
  • 刊出日期:  2008-12-24

目录

/

返回文章
返回