1. 复合材料冲击仿真概述
复合材料在冲击载荷下的响应与金属截然不同。金属靠塑性变形吸收能量，而碳纤维增强树脂基复合材料（CFRP）几乎没有屈服阶段，能量吸收主要依赖基体开裂、层间分层、纤维断裂等损伤机制。理解这些失效模式，是建立可信有限元模型的前提。
一、冲击载荷下的典型失效模式
低速冲击（通常低于10 m/s，如工具跌落、跑道碎石撞击）与高速冲击（如鸟撞、弹道冲击）在损伤形貌上有本质差异。低速冲击下，结构整体弯曲变形主导，损伤以层间分层和基体裂纹为主，表面往往只留下一个浅凹坑，内部却已大面积分层——这就是著名的"目视不可见冲击损伤（BVID）"问题。高速冲击则表现为局部穿透、纤维剪切断裂和冲塞破坏。
从损伤演化的时间尺度看，一次典型的低速冲击事件持续约1~5 ms，其中分层扩展主要发生在冲击后0.5~2 ms的卸载阶段。这一阶段层间剪应力达到峰值，是分层萌生的关键窗口。工程上常用冲击后压缩强度（CAI）来量化损伤对承载能力的影响：2 mm厚CFRP层合板在15 J冲击能量下，CAI强度可能下降40%~60%。
具体到失效模式，可归纳为五类：基体拉伸/压缩开裂、层间分层、纤维拉伸断裂、纤维压缩屈曲（kink band）、以及冲剪破坏。它们并非独立发生，而是相互耦合——基体裂纹会触发分层，分层又会导致纤维失去侧向支撑而提前屈曲。
二、LS-DYNA在冲击仿真中的优势
LS-DYNA的核心竞争力在于其显式时间积分算法。冲击问题本质上是瞬态动力学问题，惯性效应和应力波传播不可忽略。隐式方法需要求解大型刚度矩阵并迭代收敛，而显式方法通过中心差分法直接推进时间步，单步计算量小，天然适合处理接触、大变形和材料失效等强非线性问题。
具体到复合材料仿真，LS-DYNA提供了几类关键能力。第一，丰富的复合材料本构模型，包括MAT_54/55（增强型复合材料损伤模型）、MAT_58（连续损伤力学模型）、MAT_261（基于Hashin准则的渐进损伤模型）等，可分别描述面内损伤和层间分层。第二，内聚单元（Cohesive Zone Model）和基于接触的分层模拟方法，可在层间预置零厚度界面单元，通过双线性牵引-分离法则模拟分层萌生与扩展。第三，高效的接触算法，如*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE和*CONTACT_AUTOMATIC_SINGLE_SURFACE，可处理冲击头与层合板、层合板自身各层之间的复杂接触状态。
此外，LS-DYNA支持多尺度建模策略：在损伤区域采用精细壳单元或实体单元，在远场采用计算量更小的壳单元，通过*CONSTRAINED_ADAPTIVE或子模型技术实现计算效率与精度的平衡。对于工程中常见的15 J低速冲击案例，一个合理规模的模型（约5万~20万单元）在32核工作站上通常可在2~6小时内完成求解。
三、课程学习路径与工程案例导入
本课程围绕"建模—求解—判读—验证"的完整闭环展开，共30章，分为五个阶段。第一阶段（第1~6章）建立基础：复合材料力学基础、LS-DYNA显式算法原理、单元类型选择、材料模型参数获取。第二阶段（第7~14章）聚焦建模实操：几何清理、铺层定义、接触设置、边界条件与载荷施加、网格收敛性验证。第三阶段（第15~22章）深入损伤模拟：Hashin准则实现、内聚单元建模、分层模拟、损伤演化参数标定。第四阶段（第23~27章）面向工程案例：低速冲击、高速穿透、鸟撞、冲击后压缩（CAI）仿真。第五阶段（第28~30章）综合提升：实验对标、模型验证与不确定度分析。
贯穿全课程的主线案例是一块300 mm × 200 mm × 2.4 mm的CFRP层合板，铺层为[0/45/-45/90]₂ₛ，承受直径16 mm、质量5.5 kg的半球形冲击头以3.5 m/s速度的垂直冲击。这个案例来自典型的航空结构低速冲击适航验证场景，后续各章将逐步搭建其完整仿真模型，并对比实验测得的冲击力-时间曲线和超声C扫描分层形貌。
需要特别强调的是，复合材料冲击仿真的可信度高度依赖材料参数的准确性。文献中同一牌号CFRP的GIC值可能相差30%以上，直接套用他人参数往往导致分层面积预测偏差超过50%。因此，课程中会专门讲解如何通过双悬臂梁（DCB）和端部缺口弯曲（ENF）试验反推内聚单元参数，以及如何用冲击力-位移曲线的初始斜率校验面内刚度。
最后，冲击仿真不是"跑完就算"的一次性任务。一个合格的仿真报告至少应包含：网格收敛性证据、能量平衡检查（沙漏能占比应低于5%）、与至少一组实验数据的对比、以及关键参数的不确定度分析。这些要求将贯穿全部30章，逐步内化为你的建模习惯。