📚 课程目录第4章复合材料本构模型详解

4. 复合材料本构模型详解

在LS-DYNA中做复合材料冲击仿真，选错材料模型往往比网格画错更致命。冲击问题同时涉及面内纤维主导的弹性响应、层间基体开裂、厚度方向压缩与分层，单一本构很难覆盖全部机制。本节围绕四类最常用的复合材料本构——MAT_022、MAT_054/055、MAT_058、MAT_161/162——讲清它们的理论基础、适用场景和关键参数含义，帮助你在建模时做出有依据的选择。

一、MAT_022：正交各向异性弹性，只做“弹性基线”

MAT_022（*MAT_ORTHOTROPIC_ELASTIC）是最简单的正交各向异性线弹性模型。它假设材料在三个材料主方向上各自具有独立的弹性模量、剪切模量和泊松比，应力—应变关系始终线性，不产生任何损伤或失效。其本构关系可写成：

σ = C · ε
C 由 E11、E22、E33、ν12、ν13、ν23、G12、G13、G23 共 9 个独立常数确定
（LS-DYNA 中通过 ν21、ν31、ν32 或直接输入工程常数换算）

关键参数含义：E11 为纤维方向模量，通常远大于 E22、E33；G12、G13 为面内与面外剪切模量；ν12 为纤维方向拉伸时横向收缩比。MAT_022 的典型用途是：先建立无损伤的弹性基线，验证铺层方向、边界条件和载荷传递是否正确；或用于只关心刚度、不关心失效的工况。它不能预测分层、基体开裂或纤维断裂，因此冲击损伤仿真中一般只作为对照或前处理验证，而非最终损伤模型。

二、MAT_054/055：增强型复合材料损伤模型

MAT_054（*MAT_ENHANCED_COMPOSITE_DAMAGE）和 MAT_055（*MAT_ENHANCED_COMPOSITE_DAMAGE_2D）是LS-DYNA中面向壳单元的经典复合材料损伤模型。它们基于Chang-Chang类失效准则，在面内四个方向（纤维拉伸、纤维压缩、基体拉伸、基体压缩）分别判断失效，并通过刚度退化模拟损伤演化。MAT_054 支持三维应力状态，MAT_055 为二维简化版本，计算量更小。

核心参数包括：

E11、E22、G12、ν12：面内弹性常数，决定初始刚度。
XT、XC：纤维方向拉伸与压缩强度。
YT、YC：横向拉伸与压缩强度。
SC：面内剪切强度。
DFAILT、DFAILC、DFAILS：失效后刚度退化因子，控制损伤后承载能力下降幅度。

该模型的优点是参数获取相对容易、计算稳定，适合模拟层合板在冲击下的面内损伤扩展。但它对分层和厚度方向失效的刻画较弱，通常需要配合层间界面单元或接触定义使用。

三、MAT_058：连续损伤力学（CDM）模型

MAT_058（*MAT_LAMINATED_COMPOSITE_FABRIC）采用连续损伤力学框架，通过损伤变量将材料刚度随应变演化逐步退化。它比 MAT_054/055 更平滑地描述损伤累积过程，避免了突然失效带来的数值震荡。其本构中，损伤变量 d 随等效应变增长，刚度按 (1-d) 折减，并区分纤维和基体两套损伤演化。

关键参数：E11、E22、G12、ν12 为弹性常数；XT、XC、YT、YC、SC 为强度；SLIMT1、SLIMC1、SLIMT2、SLIMC2、SLIMS 等为各失效模式的极限应变或退化参数。MAT_058 适合需要观察损伤渐进扩展、且对数值稳定性要求较高的冲击仿真，例如低速冲击后压缩（CAI）分析。

四、MAT_161/162：多尺度/微观力学模型

MAT_161（*MAT_COMPOSITE_MSC）和 MAT_162（*MAT_COMPOSITE_MSC_DMG）是LS-DYNA中较新的多尺度复合材料模型。它们从纤维和基体的微观力学出发，通过均匀化方法得到宏观响应，能够更自然地描述纤维体积分数、纤维方向、基体性能对整体行为的影响。MAT_162 在 MAT_161 基础上引入损伤演化，可模拟纤维断裂、基体开裂和界面脱粘。

典型参数包括：纤维与基体的弹性模量、强度、体积分数，以及界面参数和损伤演化阈值。由于参数较多且部分需要微观实验或反演标定，MAT_161/162 更适合研究型或高精度分析场景，而非快速工程估算。

五、模型选择与参数对照

下面用一张关系图梳理四类模型在“弹性—损伤—分层—多尺度”维度上的定位，以及它们与冲击仿真关键需求的对应关系。

实际建模时，建议先用 MAT_022 跑通弹性响应，确认铺层角度和边界条件无误；再切换到 MAT_054/055 或 MAT_058 观察损伤起始与扩展；若需要研究纤维体积分数或界面性能的影响，再考虑 MAT_161/162。参数标定方面，E11、E22、G12 可由标准拉伸和剪切试验获得，XT、XC、YT、YC、SC 由相应强度试验给出，损伤演化参数则常需结合文献或反演优化确定。冲击仿真中还应关注沙漏控制、单元尺寸与损伤参数的匹配，避免因网格过粗导致损伤被“跳过”。