董 问，袁 奕，张正艺,2,3，解 德,2,3

(1. 华中科技大学 船舶与海洋工程学院，湖北 武汉 430074；2. 高新船舶与深海开发装备协同创新中心，上海 200240；3. 船舶和海洋水动力湖北省重点实验室，湖北 武汉 430074)

0 引 言

随着各国对清洁能源的需求不断增加，全球范围内天然气运输量也在快速增长。与气态形式相比，LNG 能量密度高、易于储存，非常适合陆上或海上长距离运输。LNG 船液货舱通常可分为2 种类型：薄膜型和独立型[1]。独立型液货舱又可以分为A 型、B 型和C 型[2]。B 型液货舱的货物维护系统主要包括液货舱和支撑结构，其特点在于液货舱具有部分次屏蔽且不作为压力容器。液货舱和船体通过弹性支撑结构连接。船体运动时，支撑可将作用力传递给液货舱，引起液货舱运动，进而激发舱内液体的晃荡。此时，晃荡载荷不仅会威胁液货舱的结构安全，也会使支撑结构受到的作用力增大。因此，在设计阶段，必须考虑晃荡载荷对支撑结构作用力的影响。

随着计算机性能的提高和相关理论的发展，越来越多的数值方法和求解技术被应用于液舱晃荡的研究[3]，如有限元法（Finite Element，FE）[4]，流体体积法（Volume of Fluid，VOF）[5]，任意拉格朗日-欧拉法（Arbitrary Lagrangian-Eulerian，ALE）[6]等。

本文以1 艘LNG 船B 型独立液货舱为研究对象，建立二维模型，基于ALE 方法和体积模量缩减技术，运用有限元软件Abaqus 对液舱晃荡进行分析，得到横摇状态下支撑结构作用力的时间响应。讨论了船体横摇幅值、货物状态以及液货舱挡板布置对支撑结构作用力的影响。

1 研究对象

B 型液货舱按形状可分为棱柱形（SPB 型）和球形（MOSS 型），如图1 所示[7]。本文的研究对象是棱柱形液货舱，主要由耐低温的平面板材制成。

图 1 B 型独立液货舱Fig. 1 Independent type B tank

支撑结构用于连接船体和液货舱。根据其安装位置和功能的不同，支撑结构可分为垂向支撑、防横摇支撑、防纵摇支撑以及止浮支撑等，如图2 所示。垂向支撑和止浮支撑用于支持液货舱，防止其与船体发生碰撞，可提供y 方向的作用力，且只受压不受拉。防横摇支撑和防纵摇支撑用于减轻由船体运动引起的液货舱晃动，可提供x 或z 方向的作用力。

图 2 B 型独立液货舱支撑结构Fig. 2 Support structure of independent type B tank

2 数值仿真

2.1 数值方法

本文利用有限元软件Abaqus 的显式动态分析模块（Abaqus/Explicit，Dynamic），采用任意拉格朗日-欧拉自适应网格技术（Arbitrary Lagrangian-Eulerian adaptive meshing technique）和体积模量缩减技术（Bulk modulus reduction technique），对B 型独立液货舱的晃荡问题进行研究。

使用有限元法对液舱晃荡进行分析时，若流体单元变形过大则极易导致计算不收敛。ALE 自适应网格技术可以有效解决这一问题，它结合了拉格朗日方法（网格点与物质点重合，网格之间没有物质的转移）和欧拉法（网格点固定，物质可以在网格间自由转移）的优点，能够在不改变网格原有拓扑结构的情况下保证网格质量，增大计算结果的收敛性和准确性。ALE 方法的物质导数描述为[8]：

该方法的质量守恒、动量守恒和能量守恒控制方程分别为：

本文将LNG 假定为无粘性不可压缩流体，在Abaqus 中使用Mie-Grüneisen 状态方程对其建模。一种常见的Mie-Grüneisen 状态方程表达为[9-10]：

数值分析中采用体积模量缩减技术能够在保证结果精度的情况下缩短计算时间。有限元分析时，稳态时间增量步长约等于通过网格中任意单元的膨胀波的最小传播时间，即

考虑到体积模量的表达式为：

将式(12)代入式(11)，得到：

2.2 方法验证

为验证上述数值方法的准确性，选取参考文献[15]中的二维矩形容器晃荡算例进行分析，如图3 所示。

该矩形容器在水平方向外力作用下发生受迫振动，其水平加速度可表示为：

其中：加速度幅值X0=0.002 m；角速度ω=5.5 rad/s；重力加速度g=9.8 m/s2。容器内液体为水，密度9 8 3.2 k g/m3，真实体积模量t=0 时液体完全静止。计算得到自由液面最右端A 点的液面高度随时间的变化。为考察体积模量缩减对晃荡分析的影响，选取4 个声速不同的模型进行计算，将结果与文献数据进行对比，如表1 和图4 所示。

可以看出，随着体积模量的缩减，计算时间也在减少。将体积模量缩减为真实值的1/‰，即得到的结果依然与文献数据吻合较好。因此，采用ALE 自适应网格技术和体积模量缩减技术能够对液舱晃荡问题进行有效分析，并在保证结果精度的同时缩短计算时间。

表 1 不同模型的计算时间Tab. 1 Computation time of different models

图 4 自由液面上A 点的波高变化Fig. 4 Variation of wave height at point A

2.3 有限元模型

本文以1 艘典型LNG 船B 型独立液货舱为例，建立横剖面二维模型。液货舱宽11.0 m，高6.0 m，位于中线处的纵向水密舱壁将其分为2 个体积相等的舱室。液货舱主体结构，如外板、纵向舱壁、桁材等，均由刚体单元（R2D2）建模。模型包含5 个支撑结构，中纵舱壁底部有1 个防横摇支撑，液货舱底部关于中纵舱壁对称分布有4 个垂向支撑。为方便计算，将这些支撑简化为弹性连接，由连接单元（CONN2D2）建模。防横摇支撑具有x 方向刚度，垂向支撑具有y 方向受压刚度，刚度系数取值见表2。有限元模型如图5 所示。其中，VS1～VS4 表示垂向支撑，AR1 表示防横摇支撑。

表 2 支撑结构刚度系数Tab. 2 Stiffness coefficient of support structure

图 5 B 型独立液货舱有限元模型Fig. 5 FE model of the independent type B tank

液货舱装载率为40%，LNG 由平面应变单元（CPE4R）建模。对自由表面附近的网格进行加密处理。不考虑液体和液货舱之间的摩擦力。LNG 物理属性如表3 所示[16]。

2.4 计算工况

从运动状态、货物状态和液货舱内有无挡板3 个方面考虑，选择7 种计算工况，如表4 所示。

表 3 LNG 物理属性Tab. 3 Material properties of LNG

表 4 计算工况Tab. 4 Computation case

液货舱按如下正弦规律做横摇运动：其中，横摇幅值为φ0，角速度ω=0.748 rad/s，相应的运动周期T=8.4 s，为获得稳定结果，计算时长选为5 个周期。

为验证连接单元的准确性与可靠性，计算LC1 工况，比较支撑结构作用力与液货舱自身重力。为分析不同横摇幅值对支撑结构作用力的影响，分别对φ0=5°和φ0=8°两种情况下的液舱晃荡进行计算。为分析液体和固体货物对支撑结构作用力的影响，计算LC4 和LC7工况，货物状态为液体表示装载LNG，货物状态为固体表示装载与LNG 具有相同密度和体积的固体，计算时当作刚体处理。为分析液货舱内挡板对支撑结构作用力的影响，分别计算有挡板和无挡板2 种情况。

3 计算结果

对于LC1 工况，液货舱静止时，仅有4 个垂向支撑受到y 方向的作用力。将计算得到的连接单元作用力与液货舱重力进行对比，结果如表5 所示。

表 5 LC1 工况下垂向支撑的作用力Tab. 5 Force on the vertical support structure in LC1

4 个垂向支撑的总支持力与液货舱重力的误差小于1%，说明连接单元能够较好的模拟支撑结构，传递作用力。

图6 和图7 分别对比了无挡板和有挡板的液货舱模型在不同横摇幅值下支撑结构的受力情况。可以看出，支撑结构的作用力随着船体横摇呈现出周期性变化。对于垂向支撑VS，液货舱静止时平均作用力横摇运动时作用力以为平衡位置，在其附近波动，周期与横摇运动相同。远离液

图 6 LC2 和LC5 工况支撑结构作用力对比Fig. 6 Comparison of force on supports in LC2 and LC5

图 7 LC3 和LC6 工况支撑结构作用力对比Fig. 7 Comparison of forces on supports in LC3 and LC6

货舱中心线的VS1 的作用力变化范围明显大于靠近中心线的VS2。对于防横摇支撑AR，液货舱静止时平均作用力，横摇运动时作用力以为平衡位置，波动规律与垂向支撑类似。对于所有支撑，横摇幅值越大，支撑结构作用力的变化范围越大，作用力绝对值的最大值越大。

图8 和图9 分别对比了横摇幅值不同的液货舱模型在无挡板、有挡板以及装载固体货物时支撑结构的受力情况。可以看出，无论液货舱内有挡板还是没有挡板，装载液体还是固体，支撑结构作用力均具有类似的周期性变化规律。对于所有支撑，液货舱无挡板时支撑结构作用力变化范围最大，有挡板时次之，装载固体货物时最小。这说明，由于晃荡效应的存在，相对于固体货物，液体货物使支撑结构承受了更大的作用力。同时，在液货舱设置挡板能够有效降低作用力变化范围。

图 8 LC2，LC3 和LC4 工况支撑结构作用力对比Fig. 8 Comparison of forces on supports in LC2，LC3 and LC4

图 9 LC5，LC6 和LC7 工况支撑结构作用力对比Fig. 9 Comparison of forces on supports in LC5，LC6 and LC7

4 结 语

本文利用有限元法，采用任意拉格朗日-欧拉自适应网格技术和体积模量缩减技术，对LNG 船B 型独立液货舱进行晃荡分析，并探讨了晃荡对支撑结构作用力的影响。通过与文献数据对比，验证了数值方法的可靠性；建立有二维限元模型，根据横摇幅值、液货舱布置以及货物形态选取7 个工况进行计算；比较不同工况下支撑结构作用力随时间的变化，得到以下结论：

1）船体做周期性横摇运动时，支撑结构作用力也呈现出周期性变化，且在平衡位置附近波动；

2）船体横摇幅值越大，支撑结构作用力的变化范围越大；

3）由于液舱晃荡，装载液体货物时支撑结构作用力的变化范围大于装载相同密度和体积的固体货物时的变化范围；

4）液货舱内设置挡板能够降低支撑结构作用力的变化范围。