2. LS-DYNA复合材料建模基础
在LS-DYNA中做复合材料冲击仿真，材料模型选错，后面网格再细、接触再准，结果也很难可信。复合材料与金属最大的不同在于：它是有方向性的层合结构，面内纤维方向与垂直纤维方向的模量、强度可以相差一个数量级，层间又靠基体和界面传递载荷。因此，LS-DYNA的复合材料建模必须同时回答三个问题：用什么本构描述单层板、材料坐标系如何与铺层方向对齐、每一层的厚度和角度怎样输入。这三者构成一个整体，缺一不可。
一、常用材料模型的选择依据
LS-DYNA中复合材料单层板常用模型主要有MAT_54/55、MAT_58、MAT_59，以及更复杂的MAT_261、MAT_262等。选择时先看两个维度：是否需要渐进损伤，以及是否需要考虑层间失效。
MAT_54/55（*MAT_ENHANCED_COMPOSITE_DAMAGE）是较早的增强复合材料损伤模型。MAT_54用于壳单元，MAT_55用于实体单元。它基于Chang-Chang准则，能判断纤维拉伸、纤维压缩、基体拉伸、基体压缩四种面内失效模式，并通过单元删除模拟损伤扩展。优点是参数少、计算稳定、适合快速评估；缺点是层间剪切和分层无法直接描述，损伤演化较粗糙。若只关心面内冲击损伤范围、对分层不敏感，MAT_54/55是性价比很高的选择。
MAT_58（*MAT_LAMINATED_COMPOSITE_FABRIC）是LS-DYNA中应用最广的连续壳/壳单元复合材料模型。它采用Hashin准则，支持纤维拉伸/压缩、基体拉伸/压缩以及层间剪切等多种失效模式，并可设置损伤软化参数，使损伤后刚度逐步退化而不是瞬间删除单元。MAT_58还支持应变率效应，这对冲击问题很重要。它适合模拟碳纤维层合板在低速冲击、高速冲击下的面内损伤与刚度退化，是大多数层合板冲击仿真的首选。
MAT_59（*MAT_COMPOSITE_FAILURE_SOLID_MODEL）用于实体单元，采用Tsai-Wu或最大应力准则，适合厚截面、需要三维应力状态的复合材料结构。当层合板较厚、冲击过程中厚度方向应力不可忽略，或需要模拟挤压、穿透时，可考虑MAT_59配合实体单元。
选择逻辑可以概括为：薄板面内冲击优先MAT_58；快速评估或参数不足时用MAT_54/55；厚截面三维应力显著时用MAT_59；需要显式分层时再叠加内聚单元或MAT_261/262。本节后续以MAT_58为主线讲解，因为它的坐标系和铺层输入方式最具代表性。
二、材料坐标系定义
复合材料单层板的刚度矩阵是在材料坐标系下定义的，通常记为1-2-3：1轴沿纤维方向，2轴垂直纤维方向（面内横向），3轴为单层板厚度方向。LS-DYNA中壳单元的局部坐标系与材料坐标系是两套概念，必须区分清楚。
对于壳单元，MAT_58通过积分点定义每一层的材料方向。关键字*MAT_LAMINATED_COMPOSITE_FABRIC中的BETA角用于定义参考方向，而每一层的铺层角通过*INTEGRATION_SHELL中的角度参数输入。实际纤维方向由单元局部坐标系、BETA角和铺层角共同决定。若使用*PART_COMPOSITE，则直接在部件层面定义每层厚度和角度，更为直观。
对于实体单元，材料坐标系通过*MAT_59中的AOPT和V1、V2向量定义。AOPT=2时，用V1指定1轴方向，V2指定2轴方向；AOPT=3时用角度定义。无论哪种方式，都要保证1轴与真实纤维方向一致，否则模量会张冠李戴。
三、单层板属性设置
以MAT_58为例，单层板需要输入的核心参数包括：密度RO、弹性模量EA（纤维方向）、EB（横向）、EC（厚度方向）、泊松比PRBA、PRCA、PRCB、剪切模量GAB、GBC、GCA，以及各方向的拉伸/压缩强度XT、XC、YT、YC、SC。这些参数通常来自材料供应商数据表或标准试验（如ASTM D3039、D3518）。
教学示例（合理假设，非实测）：某T700/环氧单层板，厚度0.125 mm，RO=1.6e-9 t/mm³，EA=135000 MPa，EB=9000 MPa，EC=9000 MPa，PRBA=0.30，GAB=5000 MPa，GBC=3000 MPa，GCA=5000 MPa，XT=2200 MPa，XC=1200 MPa，YT=50 MPa，YC=200 MPa，SC=80 MPa。这些数值仅用于演示输入格式，实际项目必须用真实材料数据替换。
在*PART_COMPOSITE中，每一层需要给出厚度、材料ID和角度。例如四层板[0/90/90/0]，总厚0.5 mm，可写成：
*PART_COMPOSITE
$ PID ELFORM SHRF NLOC MAREA HGID ADPOS PSETID
100 2 1.0 4
$ MID THICK BETA TFAIL SOF ...
1 0.125 0.0
1 0.125 90.0
1 0.125 90.0
1 0.125 0.0其中BETA为铺层角，单位度，绕材料3轴旋转，正负号遵循右手定则。若使用*INTEGRATION_SHELL，则需在壳单元积分点定义中逐层给出角度，并确保与*MAT_58的BETA角配合。
四、铺层角度输入方法与验证
铺层角输入有三种常见方式：直接在*PART_COMPOSITE中逐层写角度、通过*INTEGRATION_SHELL定义积分点角度、用*DEFINE_CURVE或参数化脚本批量生成。对于规则层合板，第一种最直观；对于变角度铺层或需要参数扫描，第三种更高效。
角度定义必须与材料坐标系一致。通常约定：0°层纤维沿参考方向，+45°层绕3轴逆时针旋转45°，-45°层顺时针旋转45°。若模型中参考方向与全局X轴不一致，需通过BETA角修正。建议在LS-PrePost中打开“Material Orientation”显示，检查每层纤维方向是否与设计铺层一致。
验证方法：先做一个单层板拉伸或弯曲的简单算例，对比解析解或试验值，确认模量和强度方向正确；再逐步增加层数和角度，观察损伤起始位置是否符合预期。若0°层在横向载荷下先失效，说明材料方向可能反了。
最后强调一个工程习惯：每次修改铺层或材料方向后，先跑一个极短时间的显式分析，检查能量曲线和材料方向云图，确认没有异常后再进行完整冲击仿真。复合材料建模的可靠性，往往就取决于这些看似琐碎的坐标系和角度核对。