船舶砰击分析
理论与物理
砰击的物理
船舶的砰击是什么现象？
是船舶在波浪中航行时，船底撞击水面瞬间产生极大冲击压力的现象。可分为船首部的舷外飘击、船底的底部砰击和船尾的尾部砰击。冲击压力可达MPa量级，会引起局部结构损伤和鞭状振动。
控制方程
冲击压力有理论解吗？
Wagner（1932）的理论是经典。针对二维楔形的水面撞击，
其中 $V$ 是撞击速度，$\beta$ 是入水角，$c(t)$ 是湿宽。当 $\beta \to 0$ 时压力会发散至无穷大，因此需要考虑气垫效应的修正模型。
流体侧求解不可压缩（撞击速度 $V \ll c_{water}$ 时）或可压缩的Navier-Stokes方程。自由表面用VOF法追踪。
结构侧是弹塑性FEM。根据冲击载荷的时间尺度（毫秒量级）与结构固有周期的关系，决定是动态响应还是准静态响应。
瓦格纳楔形入水理论——"水墙"为何坚硬
船底拍击海面的砰击瞬间，局部压力可达静水压的数十至数百倍。最早将其理论化的是1932年的H.瓦格纳，他推导出楔形截面以速度V入水时的最大压力关系为 p_max ≈ ½ρ(πV/2tanβ)²（β为半顶角）。关键在于"入水角越小，压力呈爆炸性增长"。例如β=5°、入水速度5m/s时，峰值压力可达数MPa。FRP制小型高速艇若持续发生砰击，船底板会发出"砰！"的声响而凹陷，原因就在于此。该理论至今仍作为砰击设计的起点被广泛使用。
各项的物理意义
- 结构-热耦合项：温度变化引起的热膨胀诱发结构变形，变形又影响温度场。$\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%以下 |
数值解法与实现
数值方法
砰击的CFD-FSI使用什么方法？
由于涉及冲击，通常采用显式解法。
| 方法 | 流体 | 结构 | 特点 |
|---|---|---|---|
| CFD-VOF + FEM显式解法 | OpenFOAM/Fluent | LS-DYNA/Abaqus Explicit | 通用。FSI多为弱耦合 |
| SPH + FEM | SPH | FEM | 无网格。擅长处理飞溅·碎波 |
| ALE (LS-DYNA) | ALE fluid | FEM | *CONSTRAINED_LAGRANGE_IN_SOLID |
| BEM + FEM | 面元法 | FEM | 高效但无法表现飞沫 |
SPH法适合砰击分析吗？
SPH无需网格，能自然追踪自由表面的大变形和飞溅。但压力场易出现数值振荡（噪声），因此常使用黎曼SPH或δ-SPH等改进版本。LS-DYNA内置了SPH求解器，可直接与FEM结构耦合。
时空分辨率
要准确捕捉冲击压力需要多高的分辨率？
砰击压力峰值在0.1~1 ms左右消失，因此时间步长需在0.01 ms以下，空间网格在撞击面附近需达到1~5 mm以下。
| 参数 | 推荐值 |
|---|---|
| 撞击面网格尺寸 | 1~5 mm |
| 时间步长 | 0.01 ms以下 |
| VOF界面分辨率 | 至少5个单元/水膜厚度 |
| 压力采样 | 0.001 ms以下 |
水池试验与CFD的"千倍壁垒"——砰击分析的现实
砰击的数值分析持续与时间步长问题作斗争。砰击峰值压力在数毫秒内产生，因此时间步长Δt必须低于0.01ms才能解析。另一方面，实际船舶越过一个波浪需要数秒。这种时间尺度的差异形成了"千倍壁垒"。实务中，主流做法是将"砰击单体的CFD"与"船体整体的长时间FEA"分离，采用将砰击压力映射的单向耦合。与水池试验比较时，实务中一般认为峰值压力误差在±20%以内就算可以接受。CFD与实验差异的主要原因往往是"气泡混入效应"，这可能使压力降低10~30%。
整体式方法
将所有物理场作为一个联立方程组同时求解。对强耦合问题稳定，但实现复杂且内存消耗大。
分区法（分离迭代法）
各物理场独立求解，在界面交换数据。易于实现，可复用现有求解器。适用于弱耦合。
界面数据传递
最近邻法（最简单但精度低）、投影法（具有守恒性）、RBF插值（对网格不一致鲁棒性强）。守恒性与精度的平衡很重要。
子迭代
在每个耦合步内进行充分迭代，确保界面条件的相容性。残差基准基于各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛系数。防止过度松弛导致发散，是一种加速收敛的自适应方法。
稳定性条件
注意附加质量效应（流体-结构耦合中结构密度≈流体密度时）。不稳定时可应用罗宾型界面条件或IQN-ILS法。
Aitken松弛的比喻
Aitken松弛类似于"平衡跷跷板"。一方推得太用力，另一方就会弹起，反作用力又导致推得更用力——为了抑制这种振荡，自动调整推力大小的就是Aitken松弛。当耦合迭代振荡不收敛时，根据前一次的修正量自动调整下一次修正量的自适应方法。
实践指南
分析步骤
实务中推进砰击分析的步骤是？
一般采用两阶段方法。
第一阶段：全船耐波性分析
- 用势流BEM或CFD计算船体运动
- 确定砰击发生条件（相对速度、相对位移）
第二阶段：局部砰击分析
- 从第一阶段结果中提取撞击速度·角度
- 用局部CFD-FSI或SPH-FEM计算冲击压力·结构响应
为什么要分两个阶段？
因为若在全船尺度上确保砰击的时空分辨率，计算成本将变得极其高昂。第一阶段确定"何时、何地、以多大速度"发生撞击，第二阶段仅对该局部事件进行高精度分析。
验证数据
有哪些可用于验证的基准问题？
楔形水面撞击问题是标准。
| 基准 | 形状 | 验证对象 |
|---|---|---|
| Wagner理论 | 2D楔形 | 冲击压力分布 |
| Zhao & Faltinsen实验 | 2D楔形(β=30°) | 压力时程 |
| Aarsnes弹性楔形 | 2D弹性板 | FSI响应 |
| Luo实验 | 3D船首剖面 | 3D压力分布 |
集装箱船的"鞭状振动"——砰击后船体整体摇晃
砰击发生后，船体整体像拨动吉他弦一样产生纵向振动的现象称为"鞭状振动"。大型集装箱船的固有周期约为2~4秒，砰击后船体中部会出现数百MPa的弯曲应力尖峰，并重复5~10次。2013年MOL Comfort沉没事故（集装箱船断成两截）中，就指出这种鞭状振动应力可能超过了设计预期。受此事故影响，国际船级社协会（IACS）修订了规定，明确要求将鞭状振动应力纳入集装箱船设计标准（UR S11A）。这成为砰击-鞭状振动耦合分析在设计上不可或缺的转折点。
分析流程的比喻
你吹过气球吗？那个瞬间，其实发生了高级的流体-结构耦合。内部气压（流体）推挤橡胶壁（结构）使其扩张→扩张的壁改变了内部压力分布→改变的压力进一步使壁变形…这种"传球"在每个计算步中反复进行，就是FSI分析。
初学者易犯的错误
"单向耦合足够了吧？"——这个判断失误在耦合分析中最危险。结构变形微小的话，单向耦合确实足够。但像心脏瓣膜开闭那样变形会大幅改变流路的情况，单向耦合就完全不行了。判断标准是"变形量是否超过特征长度的1%"。若超过，则双向耦合必不可少。若错误地采用单向耦合，结果会"看似合理实则大错特错"——这是最可怕的情况。
边界条件的思考方式
耦合界面的数据交换如同"国境的出入境管理"。各国（物理场）有各自的法律（控制方程），但若不在国境（界面）准确管理人员与物资（力·温度·位移）的往来，两国的经济（能量平衡）就会崩溃。网格不一致时的插值如同"翻译"——误译（插值误差）越小，结果越好。
软件比较
工具比较
适合砰击分析的软件是？