渐进损伤分析
理论与物理
什么是渐进损伤
老师，什么是“渐进损伤分析（PDA）”？
复合材料与金属不同，即使发生局部破坏，整体结构也可能不会崩溃。即使基体出现裂纹，纤维仍能继续承担载荷。PDA就是模拟这种损伤的阶段性扩展和载荷的重新分配。
也就是说，仅靠破坏判据（Tsai-Wu, Hashin）是不够的吗？
Tsai-Wu和Hashin是预测“初始破坏（First Ply Failure）”的。但复合材料结构的最终破坏载荷（Last Ply Failure）有时可达初始破坏载荷的2~3倍。PDA追踪的是从“初始破坏到最终崩溃”的全过程。
PDA的三要素
PDA由三个要素构成：
1. 损伤起始准则（Damage Initiation）
“损伤何时开始”。Hashin准则、Puck准则、LaRC准则等。判定各破坏模式的起始。
2. 损伤演化法则（Damage Evolution）
“损伤如何发展”。如何描述损伤起始后的刚度退化。
- 瞬时退化 — 损伤起始时刚度骤降至零（突然退化模型）。简单但网格依赖性大
- 渐进退化 — 基于断裂能量逐步降低刚度。网格依赖性小
3. 单元删除（Element Deletion）
将损伤变量达到1.0（完全损伤）的单元从分析中移除。表示材料完全破坏的区域。
渐进退化在物理上更准确吗？
是的。使用断裂能量 $G_c$ 的渐进退化，能得出与网格尺寸无关（正则化的）结果。Abaqus的Hashin Damage Evolution就是这种方法。
损伤变量
用损伤变量 $d$（0~1）表示损伤程度：
- $d = 0$: 完好（无损伤）
- $0 < d < 1$: 部分损伤
- $d = 1$: 完全破坏
Hashin的4种模式各有损伤变量吗？
是的。定义了纤维拉伸 $d_{ft}$、纤维压缩 $d_{fc}$、基体拉伸 $d_{mt}$、基体压缩 $d_{mc}$ 四个独立的损伤变量。
损伤后的刚度矩阵：
各模式的损伤独立发展，刚度矩阵的相应分量随之退化。
这就是基于CDM（连续损伤力学）的PDA。在连续体力学的框架内处理损伤，因此可以自然地融入FEM的标准框架。
总结
我来整理一下PDA的理论。
要点：
- 模拟从初始破坏到最终崩溃的全过程 — 从FPF到LPF
- 三要素 — 损伤起始准则 + 损伤演化法则 + 单元删除
- 损伤变量 $d$（0~1） — 追踪各破坏模式的损伤程度
- 基于CDM（连续损伤力学） — 可自然地融入FEM
- 基于断裂能量的正则化 — 消除网格依赖性
PDA是一种“将复合材料强度利用到极限”的分析方法呢。
能够评估复合材料后破坏强度（post-failure strength）的PDA，是提升轻量化极限不可或缺的技术。
断裂力学与渐进损伤的统一
复合材料的渐进损伤（Progressive Damage）是指随着载荷增加，局部材料损伤依次扩展的过程。1990年代，Chaboche和Ladevez等人将连续损伤力学（CDM）应用于CFRP，并公式化了损伤变量d（0=未损伤，1=完全破坏）的发展方程。基于CDM的渐进分析使得在FEM中预测最终破坏载荷成为可能，并用于量化设计的安全裕度。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出去的经历吗？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是假设“因为缓慢施加载荷，加速度可以忽略”。但在冲击载荷或振动问题中绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想要恢复原状的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——用相同的力拉铁棒和橡皮筋，哪个伸得更长？当然是橡皮筋。这种“不易伸长”的性质就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高=强度高”是不对的。刚度是“不易变形的程度”，强度是“不易破坏的程度”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（重力等）和表面力 $f_s$（压力、接触力等）。可以这样想——桥上卡车的重量是“作用于整个内部物体的力”（体积力），轮胎压路面的力是“仅作用于表面的力”（表面力）。风压、水压、螺栓预紧力…这些都是外力。这里容易犯的错误：弄错载荷方向。本想施加“拉伸”却变成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。弹一下吉他弦试试。声音会一直响吗？不，会逐渐变小对吧。这是因为振动能量通过空气阻力和弦的内部摩擦转化成了热能。汽车的减震器也是同样原理——特意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，所以设定适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入mm时，载荷·弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢：约210 GPa，铝：约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm制时为tonne/mm³（钢约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm制用N，m制也用N统一 |
数值解法与实现
PDA的实现
请告诉我PDA的具体实现方法。
两种主要的实现方法：
1. Abaqus内置的Hashin Damage
Abaqus的DAMAGE INITIATION (HASHIN) + DAMAGE EVOLUTION。使用前一页说明的设置，损伤起始→渐进退化→单元删除的全过程会自动运行。
2. 用户子程序（UMAT/VUMAT）
如果想使用更高级的损伤模型（Puck、LaRC05、自定义CDM），需要自行编写用户子程序。对于Abaqus：
VUMAT更容易实现吗？
VUMAT不需要推导切线刚度矩阵，只需根据给定的应变增量更新应力。对于PDA的首次实现，推荐使用VUMAT + Explicit。隐式解法的UMAT切线刚度推导困难，还会影响收敛性。
收敛性问题
PDA很难收敛对吧。
损伤引起的刚度退化会引发局部软化，可能导致载荷-位移路径出现回弹。对策：
| 方法 | 特点 |
|---|---|
| 粘性正则化 | 用微小粘性“平滑”软化。$\eta \approx 10^{-5}$ |
| 显式解法 | 无收敛问题（无迭代）。计算成本高 |
| 弧长法（Riks） | 可追踪回弹。设置复杂 |
| 稳定化法 | *STATIC, STABILIZE。用能量比验证 |
显式解法最稳定吗？
在不用担心收敛的意义上最稳定。但要用显式解法求解准静态问题需要质量缩放，需注意惯性效应。在实际应用中，显式解法+质量缩放正成为PDA的主流。
LS-DYNA中的PDA
在LS-DYNA中，MAT54（Chang-Chang准则）和MAT58（连续损伤力学）是PDA的标准：
| 模型 | 特点 |
|---|---|
| MAT54 | 突然退化模型。简单但网格依赖性大 |
| MAT58 | 基于CDM。渐进退化。比MAT54更准确 |
汽车碰撞分析中广泛使用MAT54/58对吧。
CFRP防撞盒和保险杠的碰撞分析，LS-DYNA MAT54/58是事实上的标准。用于预测碰撞时的能量吸收量。
总结
我来整理一下PDA的数值方法。
要点：
- Abaqus内置的Hashin Damage最便捷 — 仅需设置即可运行
- VUMAT（显式解法）是自定义PDA的推荐方法 — 无需切线刚度
- 显式解法+质量缩放是PDA的主流 — 收敛
相关主题
なった
詳しく
報告