面向变质量船舶结构的创新数学-数值模型与有限元联合评估方法及其在耙吸挖泥船上的应用

《Journal of Ocean Engineering and Science》：An innovative mathematical and numerical model combined with FEA for structural assessment of marine vessels with time-varying mass characteristics

本文针对耙吸挖泥船(TSHD)作业过程中因泥沙装载/卸载导致的快速时变质量特性所带来的复杂动态响应与结构疲劳风险，提出了一种结合改进Euler-Bernoulli梁理论、修正Kane动力学方程与多点多约束(MPC)技术的流固耦合分析方法，实现了从水动力参数(SESAM预计算)到时域动态响应(FORTRAN程序TVM_HullGirder)，再到精细化应力应变评估(ANSYS FEA)的全流程模拟，揭示了波浪方向(0°/45°)与时变质量对结构振动模态、全局运动及局部应力的显著影响，为特种船舶的结构优化与安全设计提供了高精度预测工具。

　　

在海洋工程领域，耙吸式挖泥船(Trailing Suction Hopper Dredgers, TSHD)扮演着维持航道水深、实施海岸防护与土地复垦的关键角色。然而，这类特种船舶在“空载航行-疏浚-满载航行-卸载”的循环作业中，船体质量与吃水深度会经历剧烈变化——例如卸载过程可能在几分钟内让数万吨泥沙迅速排空，导致船体刚度分布改变、重心偏移，进而引发复杂的动态响应。更棘手的是，TSHD中部为容纳泥舱而设计的结构不连续区域（如大型开口），本就削弱了整体刚度，叠加波浪、螺旋桨与柴油机激励，使得传统静态或恒定质量假设下的结构分析方法难以捕捉真实工况下的应力集中与疲劳风险。这些挑战不仅威胁船舶的结构可靠性，还可能因未预见的动态效应导致安全隐患。为此，Yerong Zhang、Yongchang Pu与Zhiqiang Hu在《Journal of Ocean Engineering and Science》发表研究，旨在开发一套创新的数学-数值模型，将时变质量特性、流体-结构相互作用与精细有限元分析相结合，为变质量船舶的结构评估提供可预测、可优化的工程框架。

研究团队采用“理论建模-数值求解-实验验证（基于工业软件数据）”的思路，核心方法包括：① 从全尺度三维有限元模型(FEM)提取等效变截面梁属性，构建考虑时变质量与附加质量的改进Euler-Bernoulli梁模型，用于描述船体梁的弯曲振动；② 基于改进的Kane动力学方程推导刚柔耦合系统的控制方程，整合六自由度(6 DOF)刚体运动与弹性变形；③ 利用边界元法(BEM)软件SESAM/HydroD预计算10种载荷工况下的频域水动力参数（附加质量A^a、阻尼系数、响应幅值算子RAO），通过插值实现时变吃水条件下的实时载荷更新；④ 开发FORTRAN程序TVM_HullGirder，采用四阶龙格-库塔(RK4)方法进行时域积分，求解刚性-柔性耦合动态响应；⑤ 借助多点多约束(MPC)技术将梁模型的位移/转角结果映射到ANSYS三维FEM，完成局部应力、应变与变形的精细化评估。研究以某138米长TSHD为案例，设定360秒卸载周期（60-300秒为质量变化阶段），对比了恒定质量/时变质量、0°迎浪/45°斜浪等场景，并通过时间步收敛性分析验证了数值稳定性。

3.1. 建模与仿真设置

船体梁被离散为20个梁单元，其中07-14号单元定义为时变质量段（模拟泥舱区域），其余为恒定质量段。结构参数源自实际设计：杨氏模量2.10×10¹¹Pa，选取前2阶模态，单元长度6.75米，通过ANSYS截面计算获取等效面积(A)、线密度(ρA)及惯性矩(I_yy、I_zz)。水动力输入基于SESAM计算的10种载荷工况（LC01满载至LC10轻载），覆盖吃水1.56米至-1.56米的变化范围；环境载荷采用JONSWAP不规则波谱（有效波高0.8米，峰值周期20秒），包含12Hz柴油机激励与螺旋桨脉动推力。仿真定义了多种对比组：动态分析类型（5/6为刚体，7/8为刚柔耦合）、水动力更新策略（0/1/2对应静态/分段/插值）、干湿模态开关（0/1对应忽略/考虑附加质量）。

3.2. 预水动力分析结果

附加质量与辐射阻尼表现出显著的吃水依赖性：纵荡、横荡与艏摇的附加质量随吃水减小而降低，而垂荡、横摇与纵摇则在轻载时更大——这与浅浸没下的自由表面效应增强有关。响应幅值算子(RAO)显示，纵摇在0°迎浪时呈现双峰特征：低频峰对应共振频率（随质量减小而降低），高频峰约0.80 rad/s（由波长与船长比λ/L_pp≈0.75时的波浪激励主导，与吃水无关）。45°斜浪条件下，船首球鼻在吃水变化时暴露引起的非线性效应更为明显，尤其在0.5-1.5 rad/s频段出现强烈的RAO波动。插值误差分析表明，LC03-LC04区间（球鼻接近水面）的非线性最强，但整体误差在工程可接受范围内。

3.3. 时变质量下的振动特性

考虑湿模态（附加质量参与）使一、二阶自然频率较干模态显著降低。随着卸载进行，质量减少导致结构柔性增加，自然频率呈非线性上升；同时，节点位置偏移、模态振幅放大，高阶模态在船首附近出现明显畸变。时变质量还引起频率变化率(dω/dt)的突变，尤其在水动力参数分段更新（方法811）时更为剧烈，而插值策略（方法821）则产生平滑过渡，凸显了连续水动力更新的重要性。

3.4. 水动力-结构动态响应

时间步收敛性测试（Δt=0.00625/0.0125/0.025秒）证实，中等步长足以捕获动态特征，谱分析显示低频频能量集中且无明显偏移。全局运动方面，刚柔耦合对垂荡与纵摇的影响微弱，但时变质量显著改变了响应轨迹——卸载后期振幅增大，且45°斜浪下垂荡更突出，0°迎浪下纵摇更敏感。结构响应上，斜浪工况的平均（Typical）与极限（ULS）位移/转角均高于迎浪，Y向位移增幅达20%-30%；时变质量条件（Case 821）进一步放大了响应，因其改变了惯性分布与刚度平衡，导致局部变形加剧。

3.5. 有限元分析与MPC映射

通过MPC技术在6个中间截面建立约束关系，将梁模型的动态结果传递至ANSYS三维壳单元模型。结果显示：45°斜浪下的von Mises等效应力比0°迎浪高20%-30%，最大总变形增加15%-25%；时变质量场景（821）的应力、弹性应变与变形均高于恒定质量场景（701），且应力集中出现在船底与舷侧连接处（如07号与12号截面）。这些峰值区域正是疲劳裂纹易萌生的高风险部位。

本研究成功构建了一套从宏观水动力到微观结构响应的完整分析链条，证明时变质量与波浪方向的耦合效应是TSHD结构评估中的关键变量。斜浪工况引发的复合载荷（横向+扭转）比迎浪更危险，而卸载过程中的质量再分配会放大动态响应，使局部应力达到不可忽视的水平。所提出的“改进梁理论+Kane方程+MPC/FEA”框架，不仅填补了变质量船舶精细化动态分析的空白，还为其他类似作业船舶（如矿砂船、浮式生产储油船）的结构设计与安全评估提供了通用方法论。未来若扩展Timoshenko梁模型以纳入剪切变形，并增加模态截断收敛性测试，将进一步增强该方法对高频与复杂耦合响应的预测能力。