基于新IGC的液化气船独立液舱内部压力计算分析

郑 雷，肖 蕾，李海洲

(1. 上海交通大学，上海 200240；2. 江南造船（集团）有限责任公司，上海 201913)

基于新IGC的液化气船独立液舱内部压力计算分析

郑 雷1,2，肖 蕾2，李海洲2

(1. 上海交通大学，上海 200240；2. 江南造船（集团）有限责任公司，上海 201913)

新版IGC明确要求应采用三维加速度椭球法计算独立液舱的内部压力，相对于传统的二维计算法，该方法能够更加准确的确定液舱各位置的内部压力，但是计算也更加复杂，尤其对于几何形状不规则的液舱。本文研究了液化气船独立液舱内部压力的计算方法，给出了三维加速度及最大液柱高度的计算公式，编制了专用计算程序，并将其应用于某型实船液舱内部压力计算，得出若干结论，对于独立液舱的结构设计具有重要的参考意义。

IGC；独立液舱；内部压力

0 引 言

国际海事组织（IMO）的“国际散装运输液化气体船舶构造与设备规则”（IGC CODE）定义的独立液舱系指自身支持的液舱，不构成船体结构的一部分，对船体强度不是必需的[1]。独立液舱分为A型、B型和C型三类：A型独立液舱主要由平面结构组成，应用传统的船舶结构分析程序的公认标准（“公认标准”系指由主管机关承认的船级社制定和保持的标准）进行设计；B型独立液舱有平面结构型式和球罐型式2种，设计应采用模型试验、精确的分析手段和分析方法确定其应力水平、疲劳寿命和裂纹扩展特性；C型独立液舱一般为球型、单耳型或双耳型，应满足压力容器的设计要求[2]。

2016年1月1号生效的新版IGC明确要求应采用三维加速度椭球法计算独立液舱的内部压力，相对于传统的二维计算法，该方法能够更加准确地确定液舱各位置的内部压力，但是计算也更加复杂，尤其对于几何形状不规则的液舱。本文以采用蝶形封头的C型双耳液舱为例，研究了三维加速度的计算方法，并分别采用解析法和离散法计算最大液柱高度，进而求得液舱内部压力，具有重要的实际应用意义。

1 内部压力计算公式

根据IGC规则，内部压力Peq由设计蒸汽压力P0和内部液体压力Pgd组成，按下式计算[2]：

内部液体压力Pgd由重力和动力加速度的联合作用引起，按下式计算：

2 aβr,βp计算公式

根据新IGC规则，船舶在北大西洋中以相当于概率水平为10-8的运动而产生的加速度分量按照下式计算[2]：

式中：V为营运速度；K为13 GM/B或1，取大者；GM为静稳心高度；L0为规范船长；Cb为方形系数；B为型宽；X为液舱重心到船中的纵向距离（重心位于船中前，X为正值，船中后，X为负值）；Y为液舱重心到船中纵剖面的横向距离；Z为液舱重心到船舶实际水线的垂向距离（重心位于水线以上，Z为正值，水线以下，Z为负值）。

如图1所示，以椭球中心为原点建立空间坐标系，x轴沿船长方向（指向船首为正），y轴沿船宽方向（指向左舷为正），z轴沿型深方向（指向上为正），则加速度椭球方程为：

过空间点（0，0，1）建立一条空间直线，并假设该直线与xoz平面的夹角为横倾角βr，与yoz平面的夹角为βp，则该直线的方向向量为，则该空间直线的方程为：

省略计算过程，可得[3]：

且横倾角βr和纵倾角βp应满足以下范围[3]：

3 Zβr,βp计算公式

由于A型、B型、C型独立液舱形状各异，且因型线布置等原因，每种液舱在基本形状的基础上又有所变化，因此如何准确、快速地求得是液舱内部压力计算的关键。

3.1 解析法

采用蝶形封头的C型双耳液舱广泛应用于中小型LPG和LNG船，其边界比主要由平面构成的A型或B型液舱复杂很多，因此本文以此为例展开研究。

如图2所示，以封头与回转轴的交点为原点建立空间坐标系，x轴沿船长方向（指向船首为正），y轴沿船宽方向（指向左舷为正），z轴沿型深方向（指向上为正）。x1，x2，x3，y1，z2，z3，R，R1与 R2均为已知。

简单分析易知，考虑到液舱的对称性，仅需计算y≥0部分的后半舱；且对位于上述区域的任意计算点，产生最大液柱高度的液舱边界点总是位于前半舱的封头上，因此只需用数学表达式对组成封头的2个曲面进行数值描述即可。

该型液舱封头由圆弧回转面和球面2部分组成，分析可知，圆弧回转面的数学表达式为：

球面的数学表达式为：

上述两曲面上任一点（xi，yi，zi）到上述平面的距离即为对应该点的液柱高度：

分别将式（14）与式（15）代入式（16），求得各自极值点分别为：

将满足曲面方程取值范围的极值点坐标代入式（16）即可求得最大液柱高度

对于主要由平面结构构成的菱形液舱（A型、B型），仅需将所有顶点坐标代入式（16）即可求得对应的最大液柱高度

3.2 离散法

当液舱边界几何形状相对复杂时，如图3所示（后半舱采用蝶形封头，前半舱采用球形封头，中间通过两部分柱面连接过渡），很难用数学表达式描述全部曲面及各自取值范围，此时便可使用离散法[5]。

离散法的核心思想是将液舱边界曲面通过有限元建模软件（MSC.Patran，HyperMesh等）离散为一定密度的节点[6]，假设任一节点坐标为（xi，yi，zi），则其至计算点（xc，yc，zc）沿三维加速度方向的高度便可使用式（16）计算[7]。

4 内部压力编程计算

对于任意类型的独立液舱，内部压力的计算流程如下：

步骤 3 按第3节中解析法或离散法求得相应的最大液柱高度

由于内部压力计算涉及大量数据迭代，无法手动完成，因此本文编制了专门的计算程序，流程如图4所示。

5 实船应用

以某型22 000 m3LPG实船3#液舱为例，采用上述计算程序计算得到的内部压力如图5所示。

6 结 语

相对于传统的二维计算法，三维计算法更加复杂[8]。本文详细研究了内部压力的计算方法，给出了三维加速度及最大液柱高度的计算公式，编制了专用计算程序，并将其应用于某型实船液舱内部压力计算。本文还得出以下结论：

1）采用解析法计算最大液柱高度的关键和难点在于如何用数学表达式描述液舱边界曲面；

2）相对于离散法，解析法具有计算量小、精度高等优点，对于形状规则的独立液舱，应首选该方法；

3）离散法的能够适用于各型边界形状复杂的液舱，缺点是计算量大、耗时长；

4）角度迭代步长及节点密度直接影响到离散法的计算时间，应选取合适的值以在保证计算精度的前提下尽量缩短计算时间。

[1]HU Ke-yi, LI Xiao-ling. The First VLGC designed and constructed in China[C]//Proceedings of the 6th PAAMES and AMEC 2014, SR-SH-04.

[2]IMO. 国际散装运输液化气体船舶构造与设备规则(IGC Code)[S]. 2014.

[3]刘文华. 中小型 LNG 船C 型独立液舱载荷分析[J]. 船舶与海洋工程, 2012(2): 1-6.LIU Wenhua. Load analysis of Type-C cargo tank in LNG carrier [J]. Naval Architecture and Ocean Engineering, 2012(2):1-6.

[4]American Bureau of Shipping. Guide for building and Classing liquefied gas carrier with independent tanks[S]. 2010, 1.

[5]Lloyd’s Register. Rules and regulations for the constructions and classification of ships for the carriage of liquefied gases in bulk[S]. 2010, 7.

[6]孙丽萍. 船舶结构有限元分析[M]. 哈尔滨: 哈尔滨工程大学出版社, 2012.

[7]刘兵山, Patran从入门到精通[M]. 北京: 中国水利水电出版社, 2002.

[8]秦斌. 液化气船独立C型液舱直接强度分析研究[C]//上海市船舶与海洋工程学会学术年会论文, 2014.

Internal pressure analysis of independent cargo tank based on New-IGC code

ZHENG Lei1,2, XIAO Lei2, LI Hai-zhou2
(1. Shanghai Jiaotong University, Shanghai 200240, China; 2. Jiangnan Shipyard (Group) Co., Ltd, Shanghai 201913, China)

As the New-IGC required, 3-D acceleration ellipsoid method should be used to calculate the internal pressure of independent cargo tanks. Comparing to traditional 2-D method, it could get the internal pressure of any location accurately, but more complicated, especially for tanks with anomalous boundary. In this paper, calculation method of internal pressure has been investigated, and the formulae of 3-D acceleration and largest liquid height were given. On purpose of improve efficiency, a dedicated program has been developed. A target ship has been selected to execute the analysis, and some conclusions have been acquired. It could be taken as a significant reference for structure design of independent cargo tank.

IGC；independent cargo tank；internal pressure

U661

A

1672-7649(2017)11-0005-04

10.3404/j.issn.1672-7649.2017.11.002

2016-12-06；

2017-02-23

郑雷(1986-)，男，工程师，主要从事船舶结构开发、设计工作。