晃动-结构耦合
理论与物理
晃荡概述
晃荡是什么现象？
是指容器内的液体因外部激励（地震、波浪、车辆加速等）而发生剧烈晃动的现象。在LNG运输船的货物舱、核电站的乏燃料池、火箭推进剂储箱等场合尤为重要。
控制方程
晃荡的数学模型是怎样的？
在线性理论中，用速度势 $\phi$ 描述。对于矩形储箱，自由液面的晃荡固有频率为：
其中 $L$ 是储箱长度，$h$ 是液面高度，$n$ 是模态阶数。
对于非线性晃荡（大振幅、伴随破碎波和冲击压力），需要用VOF法求解Navier-Stokes方程。
储箱壁受到的冲击压力（晃荡冲击）在局部可达数MPa。
与结构的耦合是如何处理的？
当储箱壁的弹性变形影响液面行为时，就需要进行流固耦合分析。特别是对于薄膜型LNG储箱（Mark III, NO96），其薄波纹板结构在晃荡冲击下会发生变形，从而改变压力分布。
晃荡的“固有频率”——储箱形状决定一切
晃荡中最可怕的是“共振”。储箱内液体的固有频率，对于长方形储箱，可通过 f₁ ≈ (1/2π)√(πg/L・tanh(πh/L)) 计算（L为储箱长度，h为液深）。当液体处于半载左右，且此固有频率与船体横摇周期一致时，振幅会急剧增大。在1960年代早期的LNG运输船设计中，曾因忽视此共振条件，导致储箱壁反复承受冲击载荷、隔热材料脱落的事故频发。此后，“保持满载或接近空载状态，避免通过共振区间”便成为业界的操作准则。
各项的物理意义
- 结构-热耦合项：温度变化引起的热膨胀诱发结构变形，变形又影响温度场。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常示例】夏天铁轨因热膨胀导致间隙变小——温度升高→热膨胀→产生应力的典型例子。电子基板在焊接后翘曲，也是不同材料热膨胀系数差异所致。发动机的气缸体因高温区与低温区的温差产生热应力，最坏情况下会导致裂纹。
- 流体-结构耦合项：流体压力及剪切力使结构变形，结构变形又改变流体域，是双向的相互作用。【日常示例】强风下悬索桥缆索振动（涡激振动）——风力使结构摇晃，摇晃的结构改变风流，进而放大振动。心脏血流与血管壁的弹性变形、飞机机翼的颤振（气动弹性失稳）也是典型的FSI问题。有时单向耦合即可满足要求，但变形较大时双向耦合必不可少。
- 电磁-热耦合项：焦耳发热 $Q = J^2/\sigma$ 引起温度升高，温度变化又改变电阻，形成反馈回路。【日常示例】电暖炉的镍铬丝通电后发热（焦耳热）变红——温度升高导致电阻变化，电流分布也随之改变。IH电磁炉的涡流发热、输电线路因温升导致的垂度增加也是此类耦合的例子。
- 数据传递项：通过插值解决不同物理场间网格不匹配的问题。【日常示例】天气预报中，将“气温数据”与“风力数据”结合计算体感温度时，若观测点不同就需要插值——在CAE耦合分析中，结构网格与CFD网格通常也不一致，因此界面处的数据传递（插值）精度直接关系到结果的可靠性。
假设条件与适用范围
- 弱耦合假设（单向耦合）：当一方物理场影响另一方，而反向影响可忽略时有效
- 需要强耦合的情况：FSI中的大变形、电磁-热耦合中温度依赖性较强时
- 时间尺度分离：各物理场特征时间差异较大时，可采用子循环提高效率
- 界面条件的协调性：需确认耦合界面处的能量和动量守恒在数值上得到满足
- 不适用的情形：当三个及以上物理场同时发生强耦合时，有时需要整体式方法
量纲分析与单位制
| 变量 | SI单位 | 注意事项・换算备忘 |
|---|---|---|
| 热膨胀系数 $\alpha$ | 1/K | 钢：约12×10⁻⁶，铝：约23×10⁻⁶ |
| 耦合界面力 | N/m²（压力）或 N（集中力） | 需确认流体侧与结构侧的力平衡 |
| 数据传递误差 | 无量纲（%） | 插值精度取决于网格密度比。建议控制在5%以下 |
数值解法与实现
数值方法
晃荡FSI分析使用哪些方法？
MPS法是日本开发的方法吗？
是由越塚诚一教授（东京大学）开发的Moving Particle Semi-implicit法。已作为Prometech公司的Particleworks产品商业化。能够稳定地追踪自由液面的大变形和飞溅，在晃荡分析方面拥有众多实际业绩。
冲击压力评估
如何评估晃荡冲击压力？
冲击压力分为气囊型和直击型。
- 气囊型：液面与壁面之间困有空气并被压缩。压力波形呈振荡性。峰值压力较低但持续时间长
- 直击型：液面沿壁面上升并直接撞击。峰值压力极高但持续时间极短
考虑压缩性时，需要进行包含空气多方过程的两相流分析。
SPH与VOF——晃荡分析的“两大流派”
在晃荡数值分析领域，使用网格的VOF法和不使用网格的SPH法长期以来相互竞争。VOF在OpenFOAM和STAR-CCM+中拥有丰富的实际业绩，气液界面的捕捉方式直观。而SPH的优势在于能自然地处理液体的大变形（如撞击储箱壁、飞溅）。在造船业界，IHI和三菱重工等公司积极将SPH用于晃荡冲击压力分析。近年来，利用GPU并行化SPH计算，可在数小时内完成百万粒子规模的仿真。答案并非哪种方法“正确”，而是根据分析目的灵活选用。
整体式方法
将所有物理场作为一个联立方程组同时求解。对于强耦合问题稳定，但实现复杂，内存消耗大。
分区法（分离迭代法）
独立求解各物理场，在界面交换数据。易于实现，可利用现有求解器。适用于弱耦合。
界面数据传递
最近邻法（最简单但精度低）、投影法（具有守恒性）、RBF插值（对网格不一致适应性强）。需权衡守恒性与精度。
子迭代
在每个耦合步内进行充分迭代，确保界面条件的协调性。残差基准基于各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛因子。是一种自适应方法，可防止过松弛导致发散，并加速收敛。
稳定性条件
注意附加质量效应（在流固耦合中，当结构密度≈流体密度时）。不稳定时，可应用Robin型界面条件或IQN-ILS法。
Aitken松弛的比喻
Aitken松弛类似于“平衡跷跷板”。一方推得太用力，另一方就会弹起，反弹又导致推得更用力——为了抑制这种振荡，自动调整推力大小的就是Aitken松弛。当耦合迭代振荡不收敛时，这是一种根据前次修正量自动调整下次修正量的自适应方法。
实践指南
LNG储箱分析步骤
LNG运输船的晃荡分析如何进行？
1. 通过耐波性分析计算船体运动（6自由度运动时程）
2. 设定储箱形状与液位（部分填充率是关键。20%~80%为危险区间）
3. CFD（VOF/SPH/MPS）晃荡分析
4. 壁面压力统计处理（短期/长期极值分布）
5. 结构响应分析（将冲击压力输入有限元模型）
为什么填充率处于中间范围时危险？
液位过低则液体质量小，过高则液体运动受限。在40%~70%的中等填充率下，晃荡最为剧烈。IGC Code（国际气体运输船规则）限制在危险填充率范围内航行。
冲击压力统计处理
如何处理冲击压力的离散性？
晃荡冲击压力在概率上表现出极大的离散性。即使在相同条件下，不同冲击事件的峰值压力也可能相差10倍以上。
采用Gumbel分布或Weibull分布进行极值统计处理，估算相对于设计寿命的最大压力。通过短期3小时的分析获取数百次冲击事件，据此估算极值分布的参数。
LNG运输船的No.96型储箱——晃荡对策的结晶
在全球LNG运输主流的薄膜型储箱中，GTT公司的No.96型采用了复杂的隔热结构设计，以吸收晃荡载荷。对于容量14万至17万立方米的大型LNG运输船，晃荡分析需要进行数百个工况的仿真，全面验证波浪条件、装载率、航行速度的组合。日本的大型造船厂（三菱重工、川崎重工）拥有专用的晃荡试验装置（在6自由度振动台上放置1/50缩尺模型），用于验证CFD结果。一艘LNG运输船的建造，其晃荡分析成本可达数千万日元规模。
分析流程的比喻
你吹过气球吗？那个瞬间，其实发生了高级的流固耦合。内部气压（流体）推动橡胶壁（结构）扩张→扩张的壁改变了内部压力分布→改变的压力进一步使壁变形……在计算中逐步重复这种“传球”过程，就是FSI分析。
初学者易犯的错误
“单向耦合足够了吧？”——这个判断失误在耦合分析中最危险。如果结构变形微小，单向耦合确实足够。但对于像心脏瓣膜开闭这样变形会大幅改变流道的情况，单向耦合完全不行。判断标准是“变形量是否超过特征长度的1%”。若超过，则双向耦合必不可少。如果错误地采用单向耦合，结果可能“看似合理实则大错特错”——这是最可怕的情况。
边界条件的思考方式
耦合界面的数据交换如同“边境的出入境管理”。各国（物理场）有自己的法律（控制方程），但若不在边境（界面）准确管理人员与物资（力、温度、位移）的往来，两国的经济（能量平衡）就会崩溃。网格不一致时的插值如同“翻译”——误译（插值误差）越小，结果越好。
软件比较
工具比较
可用于晃荡分析的工具？
なった
詳しく
報告