复合材料冲击损伤分析
理论与物理
复合材料冲击损伤
老师，复合材料受到冲击会怎么样？
金属受冲击会凹陷（塑性变形），但复合材料内部会发生破坏扩展。表面可能几乎没有痕迹，但内部会发生基体开裂、层间剥离、纤维断裂等大范围损伤。
是BVID（Barely Visible Impact Damage，目视勉强可见冲击损伤）吧。
对。飞机设计中，BVID（目视难以检测的冲击损伤）是最严苛的设计条件。工具掉落（低速冲击）可能导致BVID，在此状态下承受压缩载荷时的剩余强度（CAI）决定了设计许用值。
冲击的分类
| 种类 | 速度 | 示例 | 主要损伤 |
|---|---|---|---|
| 低速冲击 | < 10 m/s | 工具掉落、冰雹 | 基体开裂、层间剥离 |
| 中速冲击 | 10〜100 m/s | 跑道碎石、鸟撞 | 穿透损伤、大范围剥离 |
| 高速冲击 | > 100 m/s | 弹道冲击 | 穿透、塞块形成 |
| 超高速冲击 | > 1000 m/s | 太空碎片 | 凹坑、完全破坏 |
低速冲击最常见，是BVID的原因呢。
飞机运营中最常发生的就是低速冲击。设计上依据ICAO/FAA规定，假设特定能量（例如：35 J/波音，50 J/空客）的冲击会导致BVID。
冲击损伤机制
低速冲击的损伤机制（时间序列）：
1. 接触开始 — 冲击体与板接触
2. 基体开裂发生 — 弯曲应力导致正交方向裂纹
3. 层间剥离扩展 — 基体裂纹到达层间界面引发剥离
4. 纤维断裂 — 能量较大时纤维也断裂
5. 反弹 — 冲击体弹回。损伤残留
剥离是从基体开裂开始的啊。
当基体裂纹到达层间界面时，如果裂纹尖端的能量超过层间的断裂韧性，就会发生剥离。剥离优先发生在纤维角不同的层间（例如：0°/90°界面）。
FEM中的冲击分析
需要在所有层间都加入CZM吗？
理想情况是这样，但计算成本巨大。仅在主要层间（纤维角急剧变化的界面）布置CZM是实际工程中的做法。
总结
我来整理一下复合材料冲击损伤的理论。
要点：
- BVID — 表面不可见的内部损伤。飞机设计最严苛的条件
- 损伤连锁 — 基体开裂→层间剥离→纤维断裂
- 低速冲击最常见 — 工具掉落、冰雹、跑道碎石
- 用显式解法（Explicit）进行仿真 — Hashin + CZM
- 纤维角急剧变化的层间优先发生剥离
鸟撞试验与CFRP冲击设计的历史
飞机CFRP结构的鸟撞试验是美国FAR 25.571强制要求的，假设1.8kg的鸟以270km/h的速度撞击。1970年代CFRP开始应用时，其冲击特性被认为远逊于金属，但通过铺层设计和冲击后压缩强度（CAI: Compression After Impact）的优化，现已证实A350机翼的冲击耐受性优于铝材。
各项的物理意义
- 惯性项（质量项）：$\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统一 |
数值解法与实现
冲击分析的FEM模型
请告诉我冲击分析的具体设置。
Abaqus/Explicit典型的低速冲击模型：
模型构成
- 冲击体: 刚体球（半径12.5 mm, 质量3 kg, 初速5 m/s → 能量37.5 J）
- 板: 实体单元（C3D8R）。每层0.125 mm。24层 = 3 mm厚
- CZM: 内聚力单元（COH3D8）布置在主要层间
- 接触: General Contact（全自动接触）
网格
- 冲击点附近: 0.5〜1 mm单元
- 远处: 2〜5 mm单元
- 板厚方向: 每层1个单元（实体）＋层间CZM
24层每层1个单元，全层加CZM…单元数真多啊。
冲击区域设为20×20 mm的精细区域，大约50万个单元。计算时间数小时〜1天（可用GPU加速）。
损伤评估
结果确认项目：
- 力-时间曲线 — 冲击体接触力的时间历程。峰值力和接触时间
- 能量-时间曲线 — 吸收能量。计算回弹系数
- 损伤面积 — 层间剥离的投影面积。与C扫描比较
- 各层损伤变量 — Hashin的$d_{ft}, d_{mt}$云图
损伤面积与C扫描比较是验证的关键呢。
对。FEM的剥离面积与超声波C扫描的剥离面积在30%以内一致即可认为预测良好。形状（椭圆还是花生形）也要比较。
总结
我来整理一下冲击分析的数值方法。
要点：
- Abaqus/Explicit + C3D8R + Hashin + CZM — 标准配置
- 每层1单元 + 层间CZM — 精细模型
- 力-时间、能量-时间曲线验证 — 与实验比较
- 损伤面积比较 — FEM与C扫描。30%以内为良好
- 计算成本高 — 50万单元。考虑使用GPU
CAI（冲击后压缩强度）评估流程
CAI（冲击后压缩强度）是CFRP设计最重要的特性之一。ASTM D7137规定使用直径16mm的冲头以特定能量（6.7J/mm）冲击后，用同一试件进行压缩以测量剩余强度。通过改变冲击能量获取CAI曲线，确认目视勉强可见损伤（BVID: Barely Visible Impact Damage）下的CAI高于许用压缩强度。
线性单元（一阶单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切自锁（可通过减缩积分或B-bar法缓解）。
二阶单元（带中间节点）
可表现曲线变形。应力精度大幅提高，但自由度约增加2〜3倍。推荐：应力评估重要时使用。
完全积分 vs 减缩积分
相关主题
なった
詳しく
報告