层间剥离解析
层间剥离的理论基础
层间剥离是什么
老师，"层间剥离（Delamination）"是复合材料最危险的破坏模式吗？
没错。层间剥离是积层板的层与层之间剥离的破坏，是复合材料最常见和最危险的损伤模式。这是表面看不见的"内部损伤"，检测困难增加了危险性。
看不见！那太可怕了。
低速冲击（坠落物、工具坠落等）表面几乎无痕迹，但内部可能广泛发生层间剥离。这称为BVID（极难察觉的冲击损伤），是航空机复合材料设计中最重要的设计条件。
层间剥离的机制
层间剥离发生在层间剪切应力（$\tau_{xz}, \tau_{yz}$）和剥离应力（$\sigma_z$）超过界面强度时。
从断裂力学的角度有3种模式：
| 模式 | 应力 | 变形 | 代表试验 |
|---|---|---|---|
| Mode I（开口） | $\sigma_z$（拉伸） | 层间张开 | DCB |
| Mode II（面内剪切） | $\tau_{xz}$ | 层相对滑动 | ENF / 4ENF |
| Mode III（面外剪切） | $\tau_{yz}$ | 层扭转 | ECT |
Mode I最危险吗？
Mode I的临界能量释放率 $G_{Ic}$ 最低（$\approx 0.1 \sim 0.3$ kJ/m²）。Mode II的 $G_{IIc}$ 是其2〜4倍。因此Mode I开口通常先行发生。
能量释放率
层间剥离的扩展条件由能量释放率（Energy Release Rate, ERR）描述：
$G$ 是当前的能量释放率，$G_c$ 是临界值（材料特性）。
对于混合模式（Mode I + Mode II 同时作用），使用混合模式准则：
$\alpha = \beta = 1$ 是Benzeggagh-Kenane（BK）准则的特殊情况。
FEM中的层间剥离建模
怎样在FEM中建模层间剥离？
主要有2种方法：
1. VCCT（Virtual Crack Closure Technique）
从裂纹尖端的节点力和开口位移计算能量释放率。这是一种追踪现有裂纹的扩展的方法。
2. CZM（Cohesive Zone Model）
在界面放置内聚单元，用应力-开口变位的构成关系（Traction-Separation法则）表达剥离。可处理裂纹核生成和扩展两个过程。
CZM更通用吗？
CZM不需要事先假设裂纹位置（只需在所有界面放置内聚单元）。VCCT仅用于现有裂纹的扩展。CZM是当前主流。
总结
让我整理一下层间剥离的理论。
要点：
- 层间剥离是复合材料最危险的损伤 — 表面看不见（BVID）
- 3种破坏模式 — Mode I（开口）、II（剪切）、III（扭转）
- 能量释放率 $G \geq G_c$ 时扩展 — 混合模式准则
- CZM（内聚区模型）是主流 — 可处理核生成＋扩展
- VCCT仅用于现有裂纹的扩展 — 比CZM限制更多，但计算更轻
BVID（无法看见的冲击损伤）由来已久。这就像与一个"看不见的敌人"作战的设计…
完全正确。航空机复合材料设计是按照"最坏的BVID存在的前提"进行的。因此CAI强度成为设计许用值。
层间剥离的发现与航空航天的影响
CFRP积层间剥离（Delamination）问题在1970年代航空机CFRP导入时成为一个严重课题。在F-14战斗机CFRP水平尾翼上发现的积层间剥离导致设计强度下降40%，美国海军在1975年对全机进行了NDT检查。这个事件成为了CFRP设计中把"Delamination"作为主要破损模式处理的契机。
层间剥离的数值计算方法
CZM（内聚区模型）的实现
请告诉我CZM的具体实现。
CZM用Traction-Separation法则描述界面行为。
双线性的Traction-Separation法则：
1. 线性弹性阶段 — 初始刚度 $K$ 应力增长
2. 损伤开始 — 应力达到界面强度 $t^0$
3. 软化阶段 — 应力下降，开口位移增加
4. 完全分离 — 能量释放量达到 $G_c$ 破坏
需要什么参数？
Abaqus中的设置
```
*COHESIVE SECTION, RESPONSE=TRACTION SEPARATION
1.0,
*SURFACE INTERACTION, NAME=cohesive_prop
*COHESIVE BEHAVIOR
1e6, 1e6, 1e6
*DAMAGE INITIATION, CRITERION=QUADS
60., 90., 90.
*DAMAGE EVOLUTION, TYPE=ENERGY, MIXED MODE BEHAVIOR=BK, POWER=1.5
0.28, 0.79, 0.79
```
初始刚度是 $10^6$？相当大啊。
初始刚度是表现"粘合状态下无变形"的惩罚参数。太大会导致条件数恶化，收敛困难。太小会导致分离前变形。目标是 $K \approx E / t_{ply}$（层弹性率/层厚）的10〜100倍。
VCCT vs. CZM
| 特性 | VCCT | CZM |
|---|---|---|
| 裂纹核生成 | × | ○ |
| 现有裂纹扩展 | ○ | ○ |
| 网格依赖性 | 有 | 少（由 $G_c$ 正则化） |
| 参数 | 仅 $G_{Ic}, G_{IIc}$ | 强度＋ $G_c$ ＋刚度 |
| 计算成本 | 低 | 高 |
| 稳定性 | 稍微不稳定 | 更稳定 |
CZM更通用，更稳定，但参数多，成本高。
正是。VCCT适合DCB/ENF试验等简单层间剥离扩展的模拟。CZM适合冲击后损伤等复杂问题。
网格要求
CZM的网格要求是什么？
内聚区（处理区）内至少需要3〜5个单元。处理区长度为：
对于CFRP，$l_{cz} \approx 0.5 \sim 2$ mm。即单元大小0.1〜0.5 mm。
0.1 mm的单元…网格非常细致。
CZM最大的成本是网格密度。在全界面配置0.1 mm单元会导致计算膨胀。仅在预期发生剥离的界面配置CZM单元，其他界面保持结合，这是实务中的方法。
总结
让我整理层间剥离的数值方法。
要点：
- CZM是主流 — Traction-Separation法则。可处理核生成＋扩展
- 3个参数组 — 界面强度、临界能量释放率、初始刚度
- BK准则处理混合模式 — 幂律参数 $\eta$
- 网格要求严格 — 处理区内3〜5个单元（单元大小0.1〜0.5 mm）
- 仅在预期剥离界面配置CZM — 全界面计算成本太大
通过DCB和ENF试验测定GIc·GIIc
作为Delamination解析输入的积层间破坏韧性GIc由DCB试验（双悬臂梁）测定，GIIc由ENF试验（端部凹口弯曲）测定。ASTM D5528（DCB）和ASTM D7905（ENF）是标准试验规格。CFRP T800/3900-2的GIc代表值约500J/m²（Mode I），GIIc约1200J/m²（Mode II），这些值输入到CZM（内聚力模型）中。
层间剥离的实务应用
层间剥离解析的实务应用
层间剥离解析在实务中怎样应用？
主要在航空航天中，用于认证试验削减和损伤容限设计。
组件法
航空机复合材料结构的认证基于组件法：
| 级别 | 试件 | FEM的作用 |
|---|---|---|
| 样品 | 材料试验（DCB, ENF等） | 获取材料参数 |
| 单元 | 孔板、接头 | PDA（Hashin + CZM） |
| 副部件 | 面板、加强筋 | 屈曲 + 层间剥离 |
| 部件 | 机翼盒、机身截面 | 整体强度·稳定性 |
用FEM削减下位级的试验。
样品级试验是必须的（获取材料数据），但在单元〜副部件级利用FEM削减试验数。传统需要数百次试验，FEM可削减30〜50%。
CAI解析的实施步骤
航空航天最重要的层间剥离相关解析是CAI（冲击后压缩）。
1. 冲击分析（显式法） — 落锤冲击模拟。用Hashin + CZM预测损伤
2. 与无损检测比较 — 超声Cscans的剥离面积与FEM比较
3. 压缩分析（隐式/Riks） — 含损伤状态下的面内压缩
4. 预测残余压缩强度 — 用FEM与试验验证
与无损检测结果比较FEM，很有说服力。
用超声Cscans测定实际剥离面积，与FEM预测比较。一致的话就能认证"FEM能正确预测剥离"。
实务检查清单
请给我层间剥离解析的检查清单。
DCB/ENF试验的FEM重现是"入场券"。这个不合就无法继续。
完全正确。无样品级验证，进行结构级层间剥离解析没有意义。
直升机旋翼叶片的层间剥离检查
直升机旋翼叶片采用CFRP，层间剥离会改变共振特性，影响飞行安全。AgustaWestland每500小时进行一次超声探伤（A-scan·C-scan）检查层间剥离有无，FEM模拟的容许层间剥离大小（最大直径50mm、面积100cm²）为管理基准。
层间剥离的软件比较
层间剥离解析的工具
层间剥离解析有什么工具？
Abaqus和LS-DYNA是主流？
航空航天用Abaqus的CZM，汽车冲撞用LS-DYNA的*TIEBREAK。Ansys也支持，但层间剥离专题的研究论文绝大多数用Abaqus。
选择指南
层间剥离是Abaqus的"论文标准工具"。
是的。在CZM的实现品质和验证实绩上，Abaqus遥遥领先。
GENOA Progressive Failure解析工具
GENOA（Alpha STAR Corporation）是复合材渐进失效分析（含层间剥离）的专用商用工具。内部集成FEM求解器，自动将发生破坏的单元材料劣化，逐步加载。Boeing用GENOA模拟X-45无人机CFRP机翼梁的破坏，最终破坏载荷预测精度在试验值±5%以内。
层间剥离的前沿研究
层间剥离的前沿研究
层间剥离解析的前沿是什么？
有3个很活跃的方向。
面内损伤与层间剥离的耦合
实际的复合材破坏中，基体开裂（面内损伤）会诱发层间剥离。裂纹达到界面时，剥离开始。正确表现这种耦合需要完全耦合的Hashin损伤（面内）＋CZM（层间）模型。
Abaqus能同时用两者吗？
能的。在壳单元加Hashin损伤，层间加CZM单元。但计算成本膨胀，实务上仅在关注的层间配置CZM。
疲劳层间剥离
不仅是静态 $G_c$，疲劳荷载下的剥离扩展也是重要研究课题。复合材的Paris法则：
和金属疲劳裂纹扩展是一样的形式。
Abaqus实现了疲劳CZM（*DAMAGE EVOLUTION, CYCLIC）。可直接模拟循环次数对剥离的扩展。但计算成本较大。
高速冲击导致的剥离
鸟撞或弹道冲击等高速冲击，冲击波传播与剥离发生耦合。应力波在层间界面反射、透射，引张波导致剥离（spalling现象）。
高速冲击的层间剥离用LS-DYNA显式法解析。也使用SPH（光滑粒子流体动力学）耦合。
总结
层间剥离前沿研究，整理一下。
层间剥离是影响复合材可靠性的最关键课题，处于研究最前沿。
Z-pinning技术用于抑制层间剥离
Z-pinning（Z针钉）是在CFRP积层板厚度方向打入细碳纤维针（直径0.3〜0.5mm），将层间破坏韧性提高1.5〜3倍的技术。2000年代由Baggage Air Vehicles Corp（BAI）开发，Airbus A340机身补强材的接合部加强采用。针钉打入导致的面内强度降低（10〜15%）需要设计最优化权衡。
层间剥离的故障处理
层间剥离解析的故障
层间剥离解析常见故障是什么？
CZM是非线性解析中最难收敛的。
CZM不收敛
剥离开始扩展时不收敛。
应对方案（逐步进行）：
1. 粘度正则化 — *DAMAGE STABILIZATION 设 $\eta = 10^{-5}$ 左右
2. 缩小增分 — 最小增分 $10^{-10}$
3. 稳定化法 — *STATIC, STABILIZE
4. 切换显式法 — Abaqus/Explicit
经常切换到显式法？
冲撞解析原本就是显式法，没问题。但准静态问题（DCB试验等）有时也用显式法。用质量缩放大化稳定时间增分，确认运动能量/全能量 < 5%。
剥离面积与实验不符
FEM的剥离面积与实验Cscan不一致。
检查项目：
- 界面强度正确吗 — 太高则剥离迟缓，太低则快速扩展
- $G_c$ 正确吗 — DCB/ENF试验模拟是否已验证
- CZM网格够细吗 — 处理区内3个单元以上
- 层间位置正确吗 — CZM是否放在实验发生剥离的层间
"哪个层间会剥离"的预测也很难。
全层间加CZM就能自动选出最弱界面剥离，但计算成本巨大。实务上优先在纤维角急变的层间（如0°/90°界面）配置CZM。
初始刚度问题
CZM初始刚度不合适会怎样？
初始刚度目标：$K = \alpha E_3 / t_{ply}$，$\alpha = 10 \sim 100$。Abaqus默认PENALTY STIFFNESS自动计算。手动设置则在这个范围试验。
总结
层间剥离故障对应，整理一下。
"DCB/ENF不合则无法继续"。层间剥离解析的铁律。
无样品级验证，信任结构级预测是愚蠢的。这是全部FEM解析的通则。
CZM层间剥离解析中界面刚度惩罚过大的情况
CZM的惩罚刚度（cohesive stiffness K=E/t₀）过大时FEM条件数恶化，不收敛。K=10⁶〜10⁷ N/mm³（CFRP积层板界面厚t₀=0.01mm情形）为实用范围，超过则数值不稳定。不收敛时，K降1位数，单元大小细分至相干区长的1/3以下再重试。
相关主题
细节
错误