撰文/ 青岛远洋船员职业学院 张俊 贠亚杰
船舶液舱舱容计量的方法很多，包括B 样条法、tank 计算软件等，但目前使用的方法中均没有考虑液舱板架结构的变形这一因素。本文主要研究液舱板架结构在不同载况下的变形对舱容的影响： 首先使用有限元分析软件Patran 计算液舱变形，然后将变形后的节点数据导入三维建模软件CATIA 并生产新的三维模型， 由CATIA 软件计算变形后的舱容。
国际石油贸易中石油主要采用水运方式运输，船舶运输的同时也作为测量石油等液货体积的计量器具。近年来，液货舱计量的精确性越来越受到重视。由于船舶本身形状及内部构造复杂，加上船舶浮态和载况多样，影响舱容测定的因素较多，测量难度大。当船舶航行时，受到不同载况的影响，会使液货舱发生不同的变形导致舱容改变。
一、CATIA 软件运用宏命令计算舱容
对于船舶来说，装载之后的液货舱虽然变形较小，显示出来的仍旧是规则舱形状，实际上液舱各个舱壁已变成了有微小凸起或者凹进的曲面，因此我们在研究中可以将液货舱看做一个曲面。CATIA 软件曲面造型功能很强，我们可以利用该功能来完成曲面的创建。
（1） 对于节点数据较少的模型， 可以利用CATIA 中的GSD_PointSpline Loft From Excel.xls 手动操作建模，先输入点坐标，再由点连成线、由线形成面。但本文中液舱节点数量大、曲线多，手工输入方法不可行。
（2）对于大量数据曲面拟合我们也可以采用编译宏命令来实现自动读入数据。
首先导出变形后的液货舱节点数据，然后利用宏命令将节点数据导入CATIA 中，随即自动完成曲面造型。以某液舱为例，如图1 所示。完成节点数据的输入后，我们需要使节点按照一定的顺序形成样条曲线，如图2 所示。将曲线连成曲面，曲面拟合便可以得到变形后的液舱体积。
通过该方法的研究可以发现，该方法进行液货舱曲面拟合的前提是知道液货舱节点坐标，Patran 有限元计算后，虽然能够得到变形后节点坐标，但是不能直接导出变形后的节点坐标，而是分别导出变形前节点坐标以及变形后节点位移。因此利用该方法计算时必须手工将变形前的节点坐标以及变形后的位移导出，在进行节点数据的整理后才能得到坐标数据，并且整理时需要区分出内壳部分的节点坐标，操作比较繁琐，我们下面来介绍另一种更加简便的方法。
二、基于Patran 计算结果的容积计算方法
上述方法的主要问题集中在液货舱节点坐标导出这一环节上，现在介绍一种方法，避免了节点坐标导出这一环节，而是直接将变形之后的液舱模型导出，最终达到完成体积计算的目的。
首先在有限元软件Patran 中进行计算，得到强度及变形结果，然后利用Patran 软件将原模型的网格移动到变形后网格位置，再将变形之后的网格生成几何并输出到CAD软件中，就完成了变形后模型的导出。
以其中一种结构形式为例，图3 为该结构装载后的变形云图（液舱内壳云图）。
（1）以向量表示变形结果。由于目的是计算变形后液舱舱容，因此对于双壳船，只需要知道变形后内壳容积即可，模型其他部分隐藏。网格节点变形向量如图4 所示。
（2）创建FEM 场及移动网格。为将液舱节点移动到变形之后的位置，在菜单中创建FEM 场。
在工具中找到FEM-Nodes/Node modify By Field工具，并选择第二步中创建的FEM 场，就完成了网格移动到变形之后的位置这一步。变形后网格如图5 所示。
（3）将壳单元生成几何面。
（4）最后将几何体导出，由于我们使用CATIA 软件来进行后续的计算，因此导出为IGES 格式，导出时注意Patran 导出几何体的默认单位为英尺，因此需要将单位修改为毫米。
（5）导出文件并在CATIA 中打开。变形后的模型如图6 所示。
几何体导入后，进行简单的结合及封闭操作就得出该封闭几何体的体积，如图7 所示。该方法操作比较简单，没有大量的数据整理过程，可以减少错误出现的概率，可行性比较强，是一种比较简单有效的舱容计算方法，本论文中多采用该方法计算变形后的舱容。
三、以11800 吨油船为例计算液舱舱容
1.11800 吨油船模型简介
实船主要参数为： 总长为134.85m； 垂线间长为126.00m；型宽为22.00m；型深为10.60m；双层底高为1.50m；设计吃水为7.50m；结构吃水为7.80m；载重量为11800t。油船结构图如图8 所示，典型横剖面图如图9所示。
首先建立该油船的典型舱段模型即第三油舱，将该舱段作为舱容变形研究的典型舱段，并且按照规范，建模时采用1/2+1+1/2 形式建立模型，为#83 肋位至#131 肋位之间结构。建模时务必使模型更贴近真实结构以保证研究的准确度。
根据该油船的结构图资料，在Patran 软件内进行有限元建模。
2. 载荷计算与加载
根据规范规定的边界条件计算船体在满载工况下的弯矩与剪力，得出弯矩剪力值后，直接在Patran 中完成加载。其中载荷值可以采用Maxsurf 软件计算以下是计算出的#83 肋位与#131 肋位的弯矩剪力值。表1 为载荷表，表2 为边界条件表。
3. 结构及参数对舱容变形影响的计算
由于本文主要研究板厚和构件强度对于舱容变形的影响，为了充分研究，因此需要合理更改船舶各个参数，使研究结果更加全面。下面就分别更改板厚、构件强度来进行舱容变形计算。板厚为8mm 变形计算结果如图11 所示。
根据有限元计算结果，利用以上方法，将原模型网格移动到变形之后的网格位置后，导入CATIA 中计算舱容结果，如图12 所示。
由于油船装载时油舱内满载是指98% 装载，因此在计算体积时应该忽略甲板变形。模型导出时是将变形后的液舱模型导出，这其中也包含甲板变形。因此需要计算甲板变形而导致的舱容变化量，将该值与变形后舱容值加和算出舱容变形后的实际体积。甲板导致变形舱容量计算如图13所示。
依次计算9mm、10mm、11mm 和12mm 板厚时的舱容以及不同板厚时对应不同型号构件的舱容。舱容计算结果如表3 ～ 8 所示。
4. 板厚及构件强度变化时的舱容曲线
根据计算结果，我们可以绘制出相应的曲线，如图14所示。为了对比明显，将舱容变化率作为曲线的纵坐标，将构件横截面积作为横坐标。5 条曲线分别代表不同板厚，每一板厚对应四种强度不同的构件，其横截面积为15mm2 ～ 50mm2。
由图14 可以看到，板厚越薄，舱容变化曲线的曲率变化越大，构件横截面面积对舱容的影响也就越大。并且，随着板厚增加，不同板厚、相同构件之间相比较，板厚厚度越大则板厚差相同时舱容变化量曲线越密集，板厚是影响舱容变化的一个主要因素，构件强度这一因素次之。
四、结语
本文主要应用有限元分析软件Patran 和三维建模软件CATIA 计算了船舶受力变形后对船舶液货舱舱容的影响，可以更精确地测量采用船舶运输的液货体积。本文中的计算方法可以得到受力变形后的几何模型，对其他类似问题也有参考价值。
