一种用于浮式海洋结构物静水平衡状态计算的有限元方法技术

一种用于浮式海洋结构物静水平衡状态计算的有限元模型，属于海洋工程结构物设计方法领域。其特点在于：将浮式海洋结构物有限元模型中湿表面单元的浮力作为单元的固有属性，建立单元浮力修正矩阵，通过浮力修正矩阵修正结构刚度矩阵。模型的垂向位置由重力与浮力平衡关系确定，避免了施加垂向自由度边界约束。通过求解总体平衡方程，计算各节点的广义位移向量，得到计入浮态与变形耦合作用的静水平衡状态。本发明专利技术适用于结构相对刚度较低的浮式海洋结构物，如下水工作船、超大型浮体、及玻璃钢等复合材料船舶的浮态与结构变形耦合作用下的平衡状态分析，可用于精确计算上述浮式海洋结构物的浮态、总纵强度、结构变形和结构应力。 1

【技术实现步骤摘要】
一种用于浮式海洋结构物静水平衡状态计算的有限元方法
本专利技术涉及一种用于浮式海洋结构物静水平衡状态计算的有限元方法，其属于海洋工程结构物设计方法领域。
技术介绍
船舶、浮船坞、浮箱和超大型浮体等浮式海洋结构物的船体结构具有一定的柔性，在水中漂浮时受不均衡的浮力与重力作用，将发生一定的弹性变形。结构弹性变形量与浮体的浮态相关，同时弹性变形的发生将进一步引起浮体所受静水浮力分布的变化，引起浮态的进一步变化。浮式海洋结构物的浮态与结构弹性变形相互影响达到平衡，所以其全浮状态下静水平衡状态，本质上是一种流固耦合现象。现有的浮式海洋结构物浮态计算方法中，忽略浮态与结构变形的耦合问题，即在计算浮态时，将船体(广义船体，包括船舶和各类浮式海洋结构物)视为刚体，不考虑弹性变形引起的浮力变化。现有的船体结构变形分析方法中，将浮力作为固定载荷施加，不考虑结构变形进一步引起的浮力载荷的变化。对于常规运输船舶，船体结构刚度较大，结构变形对浮体分布的影响及浮力变化引起的船体结构变形均相对较小，所以采用上述简化原则加以简化，并对设计结果留有适当的安全余量，仍能够满足安全性、经济性等要求。但对于超大型浮体、浮箱、下水工作船、浮船坞等型长型深比相对较大的特殊浮式海洋结构物，以及采用玻璃钢等高强度、低刚度复合材料的船舶，船体的相对刚度显著低于常规船舶。同样载荷作用下，这类浮式海洋结构物的船体变形量显著大于钢制运输船舶，此时如果忽略船体结构变形与浮态的耦合作用，对设计中的浮态计算、总纵强度计算、结构局部强度分析等计算精度均会造成较大的影响。
技术实现思路
现有的浮式海洋结构物设计方法不能准确求解计入结构变形与浮态耦合作用的静水平衡状态。针对该问题，本专利技术提出一种用于浮式海洋结构物静水平衡状态计算的有限元方法。该有限元方法将浮式海洋结构物模型中湿表面单元的浮力作为单元的固有属性，建立单元的浮力修正矩阵，通过浮力修正矩阵修正单元的结构刚度矩阵。单元的结构刚度矩阵经过浮力修正矩阵修正后，单元受到的浮力载荷与其变形量满足流体压力原理，可准确模拟浮式海洋结构物在静水中漂浮时，结构湿表面所受压力随浮态与变形的变化而改变的问题。模型的垂向位置由重力与浮力平衡关系确定，不需要施加垂向自由度边界约束。通过求解总体平衡方程，得到各节点的广义位移向量，其中包括浮态变化引起的刚性位移及结构弹性变形，实现浮态与变形耦合求解，得到浮式海洋结构物自由漂浮时的静水平衡状态。一种用于浮式海洋结构物静水平衡状态计算的有限元方法，其
技术实现思路
如下。所述的用于浮式海洋结构物静水平衡状态计算的有限元方法，其单元刚度矩阵、浮力单元的浮力修正矩阵、总体刚度矩阵、载荷及边界条件、模型的求解方法如下。1)单元刚度矩阵。对于浮式海洋结构物湿表面的壳单元或者膜单元(简称浮力单元)，采用浮力修正矩阵修正单元的结构刚度矩阵其中，Ke为经浮力修正矩阵修正后的单元刚度矩阵，为单元的结构刚度矩阵，为浮力修正矩阵。对于模型中的非浮力单元，不对单元的结构刚度矩阵做修正。2)浮力单元的浮力修正矩阵。如图1所示，对于三角形、四边形浮力单元，C为单元的坐标中心点。将C与各边的中点Mi(1≤i≤NE)连接得到NE条直线段，NE为单元的边数。NE条直线段将单元的几何面划分为NE个面片Si。通过曲面积分算法，计算Si的面积Ai和法向对于4节点四边形壳单元、4节点四边形膜单元、8节点四边形壳单元及8节点四边形膜单元，所述浮力修正矩阵为其中，ρ为水密度，g为重力加速度。vx＝(1,0,0)，vy＝(0,1,0)，vz＝(0,0,1)分别为X轴、Y轴和Z轴的矢量方向，其中X方向、Y方向和Z方向分别为浮式海洋结构物长度方向、宽度方向及型深方向。对于3节点三角形壳单元、3节点三角形膜单元、6节点三角形壳单元及6节点三角形膜单元，所述浮力修正矩阵为所述的浮力修正矩阵，对于4节点四边形壳单元、3节点三角形壳单元，Rf取值为2，Rm取值为5，Re取值为3，Cm取值为3，Ce取值为3，对于4节点四边形膜单元、3节点三角形膜单元，Rf取值为2，Rm取值为2，Re取值为0，Cm取值为0，Ce取值为0。所述的浮力修正矩阵，对于8节点四边形壳单元、6节点三角形壳单元，Rf取值为2，Rm取值为5，Re取值为6NE+3，Cm取值为3，Ce取值为6NE+3；对于8节点四边形膜单元、6节点三角形膜单元，Rf取值为2，Rm取值为2，Re取值为3NE，Cm取值为0，Ce取值为3NE，其中，NE为单元的边数。3)总体刚度矩阵。所述用于浮式海洋结构物静水平衡状态计算的有限元方法，采用常规有限元法，将单元刚度矩阵组装，得到总体刚度矩阵。4)载荷与边界约束。所述用于浮式海洋结构物静水平衡状态计算的有限元方法，将浮力载荷作为单元的固有属性，体现于单元刚度矩阵中，所以载荷向量中不包含浮力分量。所述用于浮式海洋结构物静水平衡状态计算的有限元方法，模型的Z方向位置由隐含于总体平衡方程中的重力浮力平衡关系确定，所以边界约束中不包括Z方向线位移约束以及绕X轴、绕Y轴的转动约束。5)总体平衡方程求解。所述用于浮式海洋结构物静水平衡状态计算的有限元方法，通过求解总体平衡方程，计算各节点的广义位移向量，完成浮体平衡状态求解。广义位移向量包括两部分，一是因浮式海洋结构物浮态变化引起的位移向量，二是浮式海洋结构物结构弹性变形引起的位移向量。附图说明图1为单元曲面片分割算法示意图，其中，图1(a)为三角形单元示意图，图1(b)为四边形单元示意图。图2为实施例中浮体示意图。图3为实施例中有限元网格划分示意图。图4为实施例中单元划分及曲面片编号示意图。具体实施方式以实例说明所述一种用于浮式海洋结构物静水平衡状态计算的有限元方法的具体实施方式。如图2所示组合浮体由一块钢板与一块泡沫构成。钢板尺寸为2000mm×1000mm，厚度为10mm，其上方粘连的硬泡沫尺寸为2000mm×1000mm，厚度400mm。将该组合浮体置于淡水中，钢板中间位置处，左右两边各施加1000N集中力。钢的密度为7860kg/m3,杨氏模量为2.06×106Mpa,泊松比为0.3。水密度为1000kg/m3。重力加速度取9.8。忽略泡沫的质量，忽略泡沫因水压力产生的变形，求组合体在静水中的平衡状态。采用所述用于浮式海洋结构物静水平衡状态计算的有限元方法，计算上述浮体平衡状态具体实施方式如下。该算例中，泡沫仅起提供浮力作用，所以只建立钢板的有限元方法，把泡沫提供的浮力作为钢板单元的属性，将钢板单元定义为浮力单元。1)划分网格，建立单元结构刚度矩阵。如图3所示，将板划分为两个边长为1000mm的正方形4节点四边形壳单元e1和e2。e1依次由节点1,2,3,4构成，e2依次由节点2,5,6,3构成。对于壳单元e1和e2，结构刚度矩阵为24×24对称矩阵，且两个单元的结构刚度矩阵相同。采用常规有限元法，计算单元结构刚度矩阵KeS其中，2)建立单元的浮力修正矩阵。如图4所示，采用本专利技术提出的单元曲面片分割算法，取4个边的中点M1,M2,M3和M4，连接单元的坐标中点C与四个中点，得到四条直线段。四条直线段将单元划分为S1,S2,S3和S4四个曲面片。四个曲面片的面积均为250000mm2，法向均为(0,0,1)。建立单本文档来自技高网...

【技术保护点】
一种用于浮式海洋结构物静水平衡状态计算的有限元模型，其特征在于：1)将有限元模型中浮体湿表面的壳单元或者膜单元定义为浮力单元，将浮力作为浮力单元的固有属性，通过浮力修正矩阵修正单元刚度矩阵其中，Ke为浮力单元的单元刚度矩阵，为单元的结构刚度矩阵，为单元的浮力修正矩阵；2)对于4节点四边形壳单元、4节点四边形膜单元、8节点四边形壳单元及8节点四边形膜单元，所述浮力修正矩阵为其中，ρ为水密度，g为重力加速度。根据单元分割算法，将四边形单元划分为4片，每一片包含1个非中间节点；Ai(1≤i≤4)为第i个曲面片的面积，为第i个曲面片的平均法向；vx＝(1,0,0)，vy＝(0,1,0)，vz＝(0,0,1)分别为X轴、Y轴和Z轴的方向矢量，其中X轴、Y轴和Z轴分别为浮式海洋结构物长度方向、宽度方向及型深方向；3)对于3节点三角形壳单元、3节点三角形膜单元、6节点三角形壳单元 及6节点三角形膜单元，所述浮力修正矩阵为其中，各符号含义同上。

【技术特征摘要】
1.一种用于浮式海洋结构物静水平衡状态计算的有限元方法，其特征在于：1)将有限元模型中浮体湿表面的壳单元或者膜单元定义为浮力单元，将浮力作为浮力单元的固有属性，通过浮力修正矩阵修正单元刚度矩阵其中，Ke为浮力单元的单元刚度矩阵，为单元的结构刚度矩阵，为单元的浮力修正矩阵；2)对于4节点四边形壳单元、4节点四边形膜单元、8节点四边形壳单元及8节点四边形膜单元，所述浮力修正矩阵为其中，ρ为水密度，g为重力加速度；根据单元分割算法，将四边形单元划分为4片，每一片包含1个非中间节点；Ai，1≤i≤4，为第i个曲面片的面积；1≤i≤4，为第i个曲面片的平均法向量；vx＝(1,0,0)，vy＝(0,1,0)，vz＝(0,0,1)分别为X轴、Y轴和Z轴的方向矢量，其中X轴、Y轴和Z轴分别为浮式海洋结构物长度方向、宽度方向及型深方向；Rf为浮力修正矩阵前端空列数，Rm为浮力修正矩阵中间空列数，Re为浮力修正矩阵后端空列数，Cm为浮力修正矩阵中间空行数，Ce为浮力修正矩阵底部空行数；3)对于3节点三角形壳单元、3节点三角形膜单元、6节点三角形壳单元及6节点三角形膜单元，所述浮力修正矩阵为其中，ρ为水密度，g为重力加速度；根据单元分割算法，将三角形单元划分为3片，每一片包含1个非中间节点；Ai，1≤i≤3，为第i个曲面片的面积；1≤i≤3，为第i个曲面片的平均法向量；vx＝(1,0,0)，vy＝(0,1,0)，vz＝(0,0,1)分别为X轴、Y轴和Z轴的方向矢量，其中X轴、Y轴和Z轴分别为浮式海洋结构物长度方向、宽度方向及型深方向；实施所述有限元模型的具体过程如下：1)划分网格，建立单元结构刚度矩阵将浮式海洋结构物结构的几何模型划分为壳单元，对于4节点壳单元，结构刚度矩阵为24×24对称矩阵，采用常规有限元法，计算单元结构刚度矩阵2)建立单元的浮力修正矩阵采用单元分割算法，取4个边的中点M1,M2,M3和M4，连接单元的...

【专利技术属性】
技术研发人员：于雁云，林焰，陈明，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁;21