水弾性問題
理论与物理
水弹性物理背景
水弹性分析，在哪些领域会用到呢？
当船舶和海洋结构物承受波浪载荷时，结构物的弹性变形会影响流体力，而流体力又会进一步改变变形，这种相互作用就是水弹性分析处理的内容。这对于评估VLFS（超大型浮式结构物）和集装箱船的弹振与颤振响应至关重要。
控制方程
具体是哪些方程的组合呢？
流体侧以势流理论为基础。求解速度势 $\phi$ 的拉普拉斯方程。
线性化自由表面条件后，得到：
结构侧使用模态叠加法。将位移表示为固有模态的线性组合。
流体压力 $p = -\rho \partial\phi/\partial t$ 成为作用在结构上的外力，模态坐标的运动方程为：
其中 $a_r$ 是附加质量系数，$b_r$ 是兴波阻尼系数。这些系数通过边界元法（BEM）计算。
船舶的“砰击”——船底拍击水面的冲击理论
船舶在波涛汹涌的海面上航行时，若发生剧烈的垂荡运动，船首底部会猛烈拍击水面，这种现象称为“砰击”。在这一瞬间，船底会受到数百至数千吨的冲击力，作用时间极短，仅数毫秒。从水弹性理论的角度看，砰击是“局部水的附加质量高速变化的非定常FSI现象”。冲击力的峰值与船速的平方成正比，因此对于波高4米、航速20节的集装箱船，砰击载荷可达设计静波浪载荷的2至5倍。这种冲击直接关系到船底的疲劳裂纹，因此水弹性分析不仅是振动问题，更是疲劳寿命设计的核心。
各项的物理含义
- 结构-热耦合项：温度变化引起的热膨胀诱发结构变形，变形又影响温度场。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常示例】夏天铁轨因热膨胀导致间隙变小——温度升高→热膨胀→产生应力的典型例子。电子电路板在焊接后翘曲也是由于不同材料热膨胀系数差异所致。发动机的气缸体因高温区与低温区的温差产生热应力，严重时会导致裂纹。
- 流体-结构耦合（FSI）项：流体压力及剪切力使结构变形，结构变形又改变流体域，形成双向相互作用。【日常示例】强风下悬索桥缆索振动（涡激振动）——风力使结构摇晃，摇晃的结构改变风流，进而放大振动。心脏血流与血管壁的弹性变形、飞机机翼的颤振（气动弹性不稳定性）也是典型的FSI问题。有时单向耦合即可满足要求，但变形较大时双向耦合必不可少。
- 电磁-热耦合项：焦耳发热 $Q = J^2/\sigma$ 引起温度升高，温度变化又改变电阻，形成反馈回路。【日常示例】电炉的镍铬丝通电后发热（焦耳热）变红——温度升高导致电阻变化，电流分布也随之改变。IH电磁炉的涡流发热、输电线路因温度升高导致的垂度增加也是此类耦合的例子。
- 数据传递项：通过插值解决不同物理场间网格不匹配的问题。【日常示例】天气预报中，将“气温数据”和“风速数据”结合计算体感温度时，若观测点不同就需要插值——在CAE耦合分析中，结构网格与CFD网格通常也不一致，因此界面处的数据传递（插值）精度直接关系到结果的可靠性。
假设条件与适用范围
- 弱耦合假设（单向耦合）：当一方物理场影响另一方，而反向影响可忽略时有效
- 需要强耦合的情况：FSI中的大变形、电磁-热耦合中温度依赖性较强时
- 时间尺度分离：各物理场特征时间差异较大时，可采用子循环提高效率
- 界面条件一致性：需确认耦合界面处的能量和动量守恒在数值上得到满足
- 不适用情况：三个及以上物理场同时强耦合时，可能需要采用整体式方法
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 热膨胀系数 $\alpha$ | 1/K | 钢：约12×10⁻⁶，铝：约23×10⁻⁶ |
| 耦合界面力 | N/m²（压力）或 N（集中力） | 需确认流体侧与结构侧的力平衡 |
| 数据传递误差 | 无量纲（%） | 插值精度取决于网格密度比。建议控制在5%以下 |
数值解法与实现
基于BEM的流体分析
为什么用BEM而不用CFD呢？
在势流假设成立的情况下，BEM只需离散化物体表面即可，因此计算成本大幅降低。若使用波浪格林函数，则连自由表面的离散化也无需进行。
WAMIT、AQWA、Hydrostar是采用此方法的代表性求解器。
模态传递步骤
如何将FEM的固有模态传递给BEM呢？
在Nastran或Abaqus中执行干模态分析，然后将物体表面的节点位移映射到BEM网格上。
| 步骤 | 工具示例 | 输出 |
|---|---|---|
| FE模型创建 | MSC Patran, HyperMesh | .bdf, .inp |
| 特征值分析 | MSC Nastran SOL 103 | 固有模态 |
| 模态传递 | MpCCI, 自定义脚本 | BEM输入格式 |
| 水弹性BEM分析 | WAMIT, AQWA | 附加质量、阻尼 |
| 响应分析 | HOMER, WASIM | 模态坐标时程 |
需要考虑多少阶模态呢？
通常考虑刚体6个模态加上弹性模态10到20阶左右。对于集装箱船的弹振，垂直2节点振动模态占主导，但对于颤振，高阶模态的贡献也不可忽视。
面元法 vs. CFD——船舶水弹性分析该用哪个？
船舶水弹性计算主要有两种方法。一种是“面元法（边界元法）”——仅用面元覆盖船体表面，基于势流理论计算波浪力，计算速度快，适合设计初期阶段。另一种是求解Navier-Stokes方程的“CFD（RANS法）”，可以考虑粘性效应和大波浪的破碎，但计算成本是前者的100到1000倍。选择哪种方法取决于“你想知道什么”——对于线性波响应（设计波条件），面元法通常足够；若想评估集装箱船的甲板上浪或浮体的倾覆极限，则必须使用CFD。实际工作中，“面元法设计→CFD验证最恶劣工况”的分阶段方法是主流，全CFD分析仅限于计算资源充裕的研究项目。
整体式方法
将所有物理场作为一个联立方程组同时求解。对于强耦合问题稳定，但实现复杂，内存消耗大。
分区法（分离迭代法）
独立求解各物理场，在界面交换数据。易于实现，可利用现有求解器。适用于弱耦合。
界面数据传递
最近邻法（最简单但精度低）、投影法（具有守恒性）、RBF插值（对网格不一致鲁棒性强）。需要在守恒性和精度之间取得平衡。
子迭代
在每个耦合步内进行充分的迭代，确保界面条件的协调性。残差基准基于各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛因子。防止因过度松弛导致发散，是一种加速收敛的自适应方法。
稳定性条件
注意附加质量效应（在流体-结构耦合中，当结构密度≈流体密度时）。若不稳定，可应用Robin型界面条件或IQN-ILS法。
Aitken松弛的比喻
Aitken松弛类似于“玩跷跷板时的平衡调整”。一方推得太用力，另一方就会翘起，反弹回来又推得太用力——为了抑制这种振荡，自动调整推力大小的就是Aitken松弛。当耦合迭代振荡不收敛时，它会根据前一次的修正量自动调整下一次的修正量，是一种自适应方法。
实践指南
分析步骤概述
请告诉我从头开始进行船舶水弹性分析的步骤。
基本流程如下。
1. 创建结构FE模型（船体梁模型或3D FE模型）
2. 干模态分析（通过SOL 103获取固有频率和模态振型）
3. 创建BEM模型（湿表面面元模型）
4. 水弹性频率响应分析（计算各模态的RAO）
5. 短期/长期响应统计（结合海况数据评估疲劳和极值响应）
面元密度指南
BEM的面元尺寸如何确定？
一般准则是面元尺寸 $l_p$ 应小于波长 $\lambda$ 的七分之一。
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 面元尺寸 | < $\lambda_{min}/7$ | 取决于最短目标波长 |
| 面元数量（单舷） | 300〜3000 | 取决于船型复杂程度 |
| 与FE网格的协调 | 必需 | 影响模态传递精度 |
如何验证结果？
第一步是将刚体模态（垂荡、纵摇）的RAO与实验值进行比较。对于弹性模态，可将2节点振动的固有频率与实测值比较。DNV的基准问题（如S175集装箱船等）对验证也很有用。
超大型集装箱船会弯曲——水弹性实务中的中拱与中垂
全长超过400米的超大型集装箱船（ULCC）巨大到可以同时跨越波峰和波谷。当波峰位于船首和船尾，船中部位于波谷时发生“中拱”，反之发生“中垂”，这两种状态交替出现，船体像弓一样向上或向下弯曲。这种弯曲变形（船体梁的挠曲）会引发频率为波浪周期两倍的高谐波振动，即“弹振（springing）”，这是水弹性特有的现象。如果在设计中忽略它，疲劳寿命可能降至计算值的三分之一以下。日本的造船厂已将水弹性FEM分析作为400米级船舶总体强度评估的强制性要求，每天都要进行涉及数千种波浪条件和装载工况组合的大规模批处理计算。
分析流程的比喻
你吹过气球吗？那个瞬间，实际上发生了高度复杂的流体-结构耦合。内部气压（流体）推动橡胶壁（结构）扩张→扩张的壁改变了内部压力分布→改变的压力进一步使壁变形……在计算中，每个时间步重复这种“传球”过程的就是FSI分析。
初学者容易陷入的误区
“单向耦合应该够了吧？”——这个判断失误在耦合分析中最危险。如果结构变形微小，单向耦合确实足够。但对于像心脏瓣膜开闭这样变形会大幅改变流路的情况，单向耦合就完全不行了。一个经验法则是“变形量是否超过特征长度的1%”。如果超过，双向耦合是必须的。如果错误地使用了单向耦合，结果可能会“看起来合理但实际上大错特错”——这是最可怕的情况。
边界条件的思考方式
耦合界面的数据交换就像“边境的出入境管理”。各国（物理场）有自己的法律（控制方程），但如果在边境（界面）不能准确管理人员和物资（力、温度、位移）的往来，两国的经济（能量平衡）就会崩溃。当网格不匹配时，插值就像“翻译”——误译（插值误差）越小，结果越好。
软件比较
工具比较
请告诉我支持水弹性分析的软件。
我们来整理一下主要工具。
| 工具 | 开发商 | 方法 | 水弹性支持 |
|---|---|---|---|
| WAMIT | MIT / WAMIT Inc. | 3D面元法 | 支持广义模态 |
| AQWA | Ansys Inc. | 3D面元法 | 与Ansys Mechanical协同 |
| Hydrostar | Bureau Veritas | 3D面元法 | 与HOMER时域模块协同 |
| WASIM | DNV | Rankine面元法 | SESAM环境下的非线性时域分析 |
| OrcaFlex | Orcina | Morison/BEM | 擅长线状结构物 |
OpenFOAM + 相关主题この記事の評価 ご回答ありがとうございます！ 参考に なった もっと 詳しく 誤りを 報告 |