1. 复合材料冲击损伤仿真概述
复合材料层合板在冲击载荷下的失效行为与金属材料有本质区别。金属靠塑性变形吸收能量，损伤通常从冲击点向外以凹坑和裂纹形式扩展；而碳纤维增强树脂基复合材料（CFRP）呈各向异性、层间强度低、脆性明显，冲击能量往往在极短时间内转化为基体开裂、分层和纤维断裂等多种损伤模式，表面却可能只留下一个很浅的凹坑。这种"目视不可见、内部已严重损伤"的特征，正是复合材料冲击问题必须借助数值仿真来评估的根本原因。
冲击载荷下的典型损伤模式
理解损伤模式是选择材料模型和失效准则的前提。在低速冲击（通常指冲击速度低于 10 m/s、接触时间在毫秒量级）中，损伤按出现顺序大致可分为四类：
- 基体开裂：基体强度远低于纤维，横向拉伸或面内剪切首先使基体产生平行于纤维方向的裂纹，是分层的前驱。
- 分层：相邻铺层间因刚度不匹配产生层间剪应力和剥离应力，裂纹沿层间界面扩展。分层是低速冲击后剩余压缩强度下降的主因。
- 纤维断裂：冲击能量较高时纤维被拉断或剪断，承载能力急剧下降，常出现在冲击背面（背层受拉）。
- 纤维拔出与界面脱粘：纤维与基体界面失效，影响载荷传递效率。
高速冲击（如弹道冲击，速度数百至数千米每秒）则以纤维剪切冲塞、穿透和碎片飞散为主，应变率效应显著，需要引入应变率相关的材料参数。本课程以低速冲击损伤为主线，兼顾高速冲击的基本处理思路。
为什么选择 LS-DYNA
LS-DYNA 以显式动力学为核心，在复合材料冲击仿真中具备几项关键优势。第一，显式时间积分天然适合处理接触、大变形和短时瞬态问题，无需反复迭代求解非线性方程组。第二，内置丰富的复合材料本构，包括 *MAT_022（Chang-Chang）、*MAT_054/055（增强复合材料）、*MAT_058（LARC 失效）、*MAT_161/162（连续壳）以及基于断裂力学的 *MAT_ADD_EROSION 和 *MAT_COHESIVE 系列，可分别描述面内失效与层间分层。第三，支持壳单元、连续壳单元和实体单元混合建模，配合 *CONTACT_AUTOMATIC_SINGLE_SURFACE 与 *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE 处理层间与冲头接触。第四，可通过 *DATABASE_BINARY_D3PLOT 输出损伤变量场，结合 LS-PrePost 查看各层损伤云图。
课程学习路径
本课程按"建模—材料—接触—求解—判读"的工程链路组织，共 30 章。前段建立几何与网格基础，中段深入材料模型与失效准则，后段聚焦分层建模、接触设置、结果判读与工程验证。学习路径可概括为：
- 几何与铺层定义：建立层合板、冲头与支撑夹具，明确铺层角度与厚度。
- 单元与网格：壳单元、连续壳单元、实体单元的适用场景与网格尺寸对分层的影响。
- 材料模型：从 Chang-Chang 到 LARC，理解各参数的物理含义与取值来源。
- 层间建模：cohesive 单元、接触式分层与 VCCT 方法的选择。
- 接触与边界：冲头接触、自接触、支撑约束与初始速度施加。
- 求解控制：时间步、质量缩放、沙漏控制与能量平衡检查。
- 结果判读：损伤变量、层间应力、力-位移曲线与剩余强度评估。
工程应用场景
复合材料冲击仿真在工程中主要服务于三类问题：一是适航与安全评估，如机翼前缘、风扇机匣受鸟撞或碎片冲击后的损伤容限；二是工艺与设计优化，通过改变铺层顺序、增加层间增韧来提升抗冲击性能；三是维修决策，根据仿真预测的损伤面积判断是否需要修补或更换。这些场景对仿真的共同要求是：损伤模式判断准确、能量吸收量可信、结果可复现。
上图把本节的核心逻辑串成一条链路：冲头以给定初速度冲击层合板，能量在毫秒内转化为四类损伤模式；LS-DYNA 通过材料模型、cohesive 层间和接触设置分别描述这些模式；求解后从力-位移曲线、分层面积云图和剩余压缩强度三个角度判读结果。需要强调的是，图中速度 3 m/s、铺层 [0/90]₂ₛ、单层厚 0.125 mm 均为教学示例参数，实际工程中应根据材料数据表和工况重新确定。