氷-構造相互作用
理论与物理
现象概述
冰与结构物的相互作用，具体是什么问题？
这是评估北极海石油平台或破冰船与海冰碰撞时的载荷问题。冰会表现出破碎、屈曲、蠕变等复杂的破坏行为，因此简单的载荷模型有时是不够的。
控制方程
冰的力学模型有哪些？
冰常被建模为粘弹塑性体。在Sinha（1978）的模型中，
将应变率分解为弹性应变、延迟弹性应变（晶界滑移）和粘性蠕变应变三个分量。在高应变速率（碰撞时）也使用Maxwell-Brittle型模型。
结构侧是通常的弹塑性FEM。作为冰-结构接触问题，使用罚函数法或增广拉格朗日法处理界面。冰的破碎用单元删除或内聚力模型来表现。
也有冰载荷的经验公式吗？
ISO 19906（北极海上结构）中，对于接触面积 $A$ 的冰压力，给出公式：
其中 $C_R$ 是冰参考强度，$h$ 是冰厚。数值模拟的验证中，与这个经验公式的比较很有用。
冰“比金属更复杂”——随温度、速率、盐分变化的冰力学特性
在结构分析中处理冰时，首先遇到的壁垒就是“该用什么材料模型来模拟冰”这个问题。冰在-2℃和-20℃时强度相差一倍以上，应变速率慢时会发生蠕变（粘性流动），快时则发生脆性破坏。此外，海冰含有盐分，比纯冰更弱，即使是同样的“冰”，特性也大不相同。工程上认为“应变速率导致破坏模式转变”的边界大约在10⁻³/s左右，船舶碰撞（高速）时以脆性破坏为主，桥墩冰压（低速）时则以蠕变为主。用一个材料模型来表现这两种破坏模式，是冰结构相互作用模拟的理论核心。
各项的物理意义
- 结构-热耦合项：温度变化引起的热膨胀诱发结构变形，变形又影响温度场。$\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%以下为目标 |
数值解法与实现
离散化方法
冰的破碎在数值上如何处理？
主要有三种方法。
| 方法 | 特点 | 适用软件 |
|---|---|---|
| FEM + 单元删除 | 达到破坏准则时删除单元 | LS-DYNA, Abaqus/Explicit |
| DEM（离散元法） | 将冰表现为粒子集合体 | PFC, YADE |
| SPH | 无网格法。易于追踪破碎 | LS-DYNA SPH |
| 近场动力学 | 非局部模型。裂纹自然产生 | Peridigm |
用DEM处理冰听起来很有趣。
DEM将冰表现为大量圆盘（2D）或球体（3D）的集合体，粒子间的键在超过破坏准则时断裂。其优点是可以自然地再现冰的破碎模式（径向裂纹、环向裂纹）。
接触算法
冰与结构物的接触如何处理？
在LS-DYNA中使用*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE定义摩擦接触。冰的摩擦系数强烈依赖于温度，在 $\mu = 0.01$〜$0.3$ 的范围内变动。
冰与结构物的接触如何处理？
在LS-DYNA中使用*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE定义摩擦接触。冰的摩擦系数强烈依赖于温度，在 $\mu = 0.01$〜$0.3$ 的范围内变动。
显式解法的时间步长遵循由最小单元尺寸和声速决定的Courant条件。
冰中的声速约为3,000 m/s，因此单元尺寸为0.01 m时，$\Delta t \approx 3 \times 10^{-6}$ s，非常短。
用DEM（离散元法）计算“冰群”——FEM无法实现的大规模破冰
破冰船实际航行时，破碎后的无数冰片会持续在船体周围流动。如果用FEM（有限元法）一个个地建模这些“碎冰片群”，单元数将达到数百万至数千万，计算将不可能。这时使用的就是DEM（离散元法）。DEM将每个冰片表现为“刚体+弹簧+阻尼器”的简单模型，高效计算冰片之间、冰片与船体之间的接触力。挪威SINTEF研究所进行了包含100万个以上冰片的模拟，使破冰船阻力的模拟结果与实测值误差在±15%以内。DEM不仅用于冰，还广泛用于砂、岩石、颗粒材料的分析，其在冰-结构问题上的发展也促进了颗粒体力学整体技术的提升。
整体式法
将所有物理场作为一个联立方程组同时求解。对于强耦合问题稳定，但实现复杂，内存消耗大。
分区法（分离迭代法）
独立求解各物理场，在界面交换数据。易于实现，可利用现有求解器。适用于弱耦合。
界面数据传递
最近邻法（最简单但精度低）、投影法（具有守恒性）、径向基函数插值（对网格不一致鲁棒性强）。守恒性与精度的平衡很重要。
子迭代
在每个耦合步内进行充分迭代，确保界面条件的协调性。残差基准基于各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛系数。防止因过度松弛导致发散，是加速收敛的自适应方法。
稳定性条件
注意附加质量效应（流体-结构耦合中结构密度≈流体密度时）。不稳定时可应用Robin型界面条件或IQN-ILS法。
Aitken松弛的比喻
Aitken松弛类似于“平衡跷跷板”。一方推得太用力，另一方就会弹起，其反作用力又导致推得更用力——为了抑制这种振荡，自动调整推力大小的就是Aitken松弛。当耦合迭代振荡不收敛时，它会根据上次的修正量自动调整下一次的修正量，是一种自适应方法。
实践指南
模型构建步骤
请告诉我开始冰-结构相互作用模拟的步骤。
1. 创建结构物的3D FE模型（钢结构时用壳单元）
2. 创建冰板模型（实体单元。若要处理破坏，需添加单元删除设置）
3. 定义接触（面面接触，设定摩擦系数）
4. 设定冰的初始速度/漂移速度
5. 定义材料模型（冰：Tsai-Wu破坏准则等，结构：弹塑性）
6. 用显式解法执行
冰的材料参数
冰的物性值如何确定？
海冰的物性强烈依赖于温度、盐分浓度和应变速率。
| 参数 | 一年冰（-10°C） | 多年冰（-10°C） |
|---|---|---|
| 杨氏模量 | 3〜9 GPa | 5〜10 GPa |
| 抗压强度 | 2〜10 MPa | 5〜15 MPa |
| 抗拉强度 | 0.5〜2 MPa | 1〜3 MPa |
| 泊松比 | 0.33 | 0.33 |
| 密度 | 900 kg/m³ | 910 kg/m³ |
物性的变化范围很大呢。
因此参数化研究是必须的。ISO 19906规定使用50年重现期的值作为特征值。有时也通过蒙特卡洛模拟进行概率论评估。
破冰船设计——通过模拟优化“破冰过程”
破冰船并非简单地“撞上去割开”，而是通过船首的倾斜角使冰发生“弯曲破坏”，从而高效地破冰。倾斜角太浅，冰不会破碎而会潜入船底；太陡，则会将冰向压缩方向推压，产生极高的载荷。最优倾斜角随冰的厚度和强度而变化，因此在实务中，会以“连续破碎厚度1.5m、抗压强度2MPa的海冰”为设计条件进行模拟，优化船首形状与推进力的组合。据报道，在俄罗斯的原子号（核动力破冰船）设计中，破冰模拟的计算结果直接导致了船首形状的修改，使破冰能力比初期设计提高了20%。
分析流程的比喻
你吹过气球吗？那个瞬间，实际上发生了高级的流体-结构耦合。内部气压（流体）将橡胶壁（结构）撑开→撑开的壁改变了内部压力分布→改变的压力进一步使壁变形……在每个计算步中重复这种“投接球”的就是FSI分析。
初学者容易陷入的误区
“单向耦合足够了吧？”——这个判断失误在耦合分析中最危险。如果结构变形微小，单向耦合确实足够，但像心脏瓣膜开闭那样变形会大幅改变流路的情况，单向耦合就完全不行了。判断标准是“变形量是否超过特征长度的1%”。超过的话，双向耦合是必须的。如果误用单向耦合，结果会“看起来合理但实际上大错特错”——这是最可怕的模式。
边界条件的思考方式
耦合界面的数据交换就像“国境的出入境管理”。各国（物理场）有自己的法律（控制方程），但如果在国境（界面）不能准确管理人和物（力、温度、位移）的往来，两国的经济（能量平衡）就会崩溃。网格不一致时的插值就像“翻译”——误译（插值误差）越小，结果越好。
相关主题
なった
詳しく
報告