晃荡-结构耦合
晃荡-结构耦合的理论基础
晃荡的概述
什么是晃荡现象？
容器内的液体因外部励振（地震、波浪、车辆加速等）而激烈摇晃的现象。在液化天然气（LNG）运输船货舱、原子能发电站乏燃料池、火箭推进剂储箱等中变得重要。
支配方程
晃荡的数学模型是什么？
在线性理论中用速度势 $\phi$ 来描述。对于矩形储罐，自由表面晃荡固有频率为：
其中 $L$ 为储罐长度，$h$ 为液面高度，$n$ 为模态阶次。
对于非线性晃荡（大振幅、破碎、伴随冲击压力），用VOF法求解Navier-Stokes方程。
储罐壁的冲击压力（晃荡撞击）局部可达数兆帕。
与结构的耦合如何处理？
当储罐壁的弹性变形影响液面行为时需要进行流固耦合分析。特别是膜式LNG储罐（Mark III、NO96）的薄波纹板结构在晃荡冲击下会发生变形，改变压力分布。
晃荡的"固有频率"——储罐的形状决定一切
晃荡中最可怕的是"共振"。储罐内液体的固有频率可由 f₁ ≈ (1/2π)√(πg/L·tanh(πh/L)) 计算（L为储罐长度，h为液深）。当液体充填约一半时，该固有频率与船舶横摇周期相同，会导致振幅爆炸性增长。1960年代初期的LNG运输船设计中，许多项目忽视了这一共振条件，导致储罐壁遭受重复冲击荷载，保温材料脱落事故频繁发生。此后，"保持满载或接近空载状态，避免通过共振范围"成为业界的操作规则。
晃荡-结构耦合的数值计算方法
数值手法
晃荡流固耦合分析采用什么方法？
MPS法是日本开发的方法吗？
是的。由东京大学越塚诚一教授开发的Moving Particle Semi-implicit（MPS）法。已由Prometech公司商业化为Particleworks。在自由表面的大变形和飞溅追踪上表现稳定，在晃荡分析中有丰富实绩。
冲击压力的评估
晃荡冲击压力如何评估？
冲击压力有两种类型：空气囊型（压缩空气缓冲）和flip-through型（直接撞击）。
- 空气囊型：液面与壁面之间的空气被压缩，产生振荡压力波形。峰值压力较低但持续时间较长
- Flip-through型：液面沿壁面上升，直接撞击。峰值压力非常高但持续时间短
考虑压缩性时，需要包括空气多项式过程的两相流分析。
SPH与VOF——晃荡分析的"两大派系"
在晃荡数值分析中，使用网格的VOF（Volume of Fluid）法和不使用网格的SPH（粒子法）长期竞争。VOF在OpenFOAM和STAR-CCM+中有丰富的实绩，气液界面的捕捉方法直观。SPH则在处理液体大变形（与储罐壁碰撞和飞溅）时更具优势。造船业中，日本的IHI和韩国三星重工在晃荡冲击压力分析中积极应用SPH。最近还出现了在GPU上并行化SPH计算的技术，可在数小时内完成百万粒子级模拟。在实务中，"哪种正确"不是关键，而是根据分析目标选择合适工具。
晃荡-结构耦合的实务应用
LNG储罐的分析步骤
LNG运输船的晃荡分析如何进行？
1. 通过耐航性分析计算船体运动（6自由度运动的时间历程）
2. 储罐形状和液位设定（部分充填率很关键。20～80%为危险区间）
3. 进行CFD（VOF/SPH/MPS）晃荡分析
4. 进行壁面压力的统计处理（短期、长期极值分布）
5. 进行结构响应分析（将冲击压力输入FE模型）
为什么中等充填率最危险？
液位太低时液体质量小，太高时液体运动受限。40～70%的中等充填率下晃荡最剧烈。国际气体运输船规则（IGC Code）对这一危险充填率区间的航行进行了限制。
冲击压力的统计处理
冲击压力的离散程度如何处理？
晃荡冲击压力在统计上呈现极大的离散性。在相同条件下，不同冲击事件的峰值压力可相差10倍以上。
采用Gumbel分布或Weibull分布进行极值统计处理，对设计寿命下的最大压力进行推定。从短期3小时分析中获得数百次冲击事件，推定极值分布的参数。
LNG运输船的No.96型储罐——晃荡对策的结晶
占全球LNG运输主流的膜式储罐中，GTT公司的No.96型通过复杂的绝热结构来吸收晃荡荷载。大型LNG运输船（容量14～17万立方米）的晃荡分析需要数百个仿真案例，组合测试波浪条件、装载率和航行速度。日本大型造船厂（三菱重工、川崎重工）拥有专用晃荡实验装置（6自由度振动台加1/50缩尺模型），用于验证CFD结果。单艘LNG运输船建造所需的晃荡分析成本可达数千万日元。
晃荡-结构耦合的软件比较
工具比较
晃荡分析可用的工具有哪些？
为什么Particleworks在LNG业界广泛应用？
MPS法能稳定追踪自由表面的大变形和破碎。GUI可轻松设置储罐摇晃条件，还具备壁面压力自动统计处理功能。三菱重工、三星重工等造船厂的采用实绩丰富。
包括流固耦合时如何处理？
LS-DYNA的SPH-FEM耦合最直接。将Particleworks或FLOW-3D的结果以压力数据形式导出，输入到FEM求解器（Nastran、Abaqus等）作为mapped load的单向耦合在实务中也很常见。
晃荡分析软件的市场竞争——专用工具vs通用工具
晃荡分析工具市场呈现有趣的竞争格局。ANSYS Fluent、STAR-CCM+等通用CFD占有较大份额，同时船舶专用工具如DNV Wasim也获得稳固支持。海洋和石油巨头（Shell、Total、BP）多采用STAR-CCM+，而LNG运输船专业造船厂倾向使用DNV-GT（GTT技术评价工具）的组合。国内方面，国土交通省"晃荡评价指南"多采用ANSYS案例，行政提交文件中ANSYS系更具优势。但GTT的专利管理严格，限制了独立设计储罐的普及。
晃荡-结构耦合的先端研究
极低温FSI
液化天然气（LNG）的温度是-162°C。需要考虑极低温的影响吗？
液化天然气的物性（密度、粘度、表面张力）随温度变化。此外，液化天然气的蒸发（BOG：boil-off gas）产生的气相存在也影响晃荡。气液两相晃荡中，气相压缩性改变空气囊型冲击压力。
绝热材料（聚氨酯泡沫等）在极低温下的脆性断裂也应与晃荡冲击组合评估。GTT的Mark III和NO96储罐中，绝热材料的冲击响应是设计的关键。
对比晃荡实验
与实验的对比如何进行？
GTT公司的大型晃荡试验装置（6自由度摇晃台+缩尺储罐）为业界标准验证数据源。此外，ISOPE（国际海洋与极地工程师学会）的晃荡基准研讨会定期进行多个代码的对比。
减振装置
有抑制晃荡的方法吗？
隔板（挡板）和TLD（调谐液体阻尼器）是代表方法。隔板形状最优化利用CFD。
其中 $A_b$ 为隔板面积，$C_D$ 为阻力系数。隔板的开口率和位置的最优化决定晃荡抑制效果。
火箭燃料储箱的晃荡——阿波罗计划的隐藏难题
晃荡不仅是海洋问题。在阿波罗计划中，液氧和液氢储箱内燃料晃荡威胁飞行稳定性。火箭姿态控制系统对燃料摇晃产生反应，反复过度修正，存在控制不稳定风险。美国宇航局（NASA）称之为"推进剂晃荡"，在储箱内设置隔板来抑制晃荡。现在的SpaceX Falcon 9在大型推进剂储箱的晃荡抑制方面仍是重要设计课题，CFD-FSI耦合模拟与姿态控制算法的开发并行进行。
晃荡-结构耦合的故障处理
冲击压力的再现性问题
用相同条件计算时，冲击压力每次都不同。
这是晃荡冲击的本质特性。初始条件的微小差异导致破碎形状改变，产生冲击压力的大幅波动。应对措施为：
- 执行多次计算并进行统计评估
- 使用空间平均压力（面积压力）
- 采用Weibull/Gumbel分布进行极值统计处理
网格密度与冲击压力的关系？
峰值压力对网格有很强的依赖性。冲击力（积分值）对网格更稳定。将结构响应作为最终评估指标，由于结构的惯性滤波作用，高频成分衰减。
液面追踪的精度
VOF法中液面扩散的问题。
| 对策 | 详细 |
|---|---|
| 界面压缩方案 | 启用OpenFOAM的interFoam的界面压缩（cAlpha=1） |
| AMR | 液面附近网格动态细分化 |
| HRIC/CICSAM | 选择高精度VOF方案 |
| Level Set + VOF | CLSVOF法改善界面锐度 |
计算耗时过长的处理方法？
晃荡分析需要长时间计算（物理时间数千秒以上）。通过GPU求解器（Particleworks支持多GPU）和AMR动态网格优化来加速。另外，等效不规则励振转化为短时间规则励振的等效规则波方法也能缩短计算时间。
晃荡分析的"压力峰值再现性问题"
在晃荡冲击压力的CFD分析中，经常出现"相对实验值，峰值压力差异大"的问题。用同一实验条件、同一求解器多次重复计算，冲击压力的峰值会产生±50%以上的变动。原因之一是混沌液面行为——初始条件的微小差异将液面形状拉向完全不同的状态，出现"蝴蝶效应"。另一个原因是计算网格大小，网格越粗，峰值压力越易被低估。实务中采取的对策是"评估大量分析案例的统计分布""输入不规则波而非规则波"等方法逐步标准化。
有用
更详细
报告