3. LS-DYNA复合材料建模流程
复合材料冲击仿真的建模流程与金属结构有本质差异：铺层方向、层间界面、失效准则和单元厚度积分都会直接影响计算结果。本节按一条可复现的建模链路展开，从几何与网格准备，到材料、单元、边界与求解控制，给出每一步的关键输入与判断依据。示例参数均为教学假设值，用于说明配置方式，不代表特定材料的实测数据。
一、前处理工具选择
LS-DYNA本身是求解器，几何清理、网格划分和关键字生成通常由前处理器完成。常见组合有三类：
- LS-PrePost：LS-DYNA自带，关键字支持最完整，适合检查
*PART_COMPOSITE、*SECTION_SHELL、*CONTACT等卡片，也便于后处理云图与层间结果提取。适合作为最终检查与提交工具。 - ANSYS Workbench + LS-DYNA：适合已有ANSYS环境的用户，可用ACP模块做铺层定义，再导出关键字。优点是铺层可视化好，缺点是部分关键字需手工补写。
- HyperMesh / Abaqus + 转换：HyperMesh的复合材料铺层工具成熟，Abaqus模型可通过转换工具导入，但需核对单元类型与材料卡片的对应关系。
无论选哪种工具，最终都要回到关键字文件核对三件事：铺层角度序列、单元法向与厚度积分点、接触与边界卡片是否完整。建议以LS-PrePost作为提交前的统一检查入口。
二、材料模型分类
复合材料冲击问题中，材料模型按“是否考虑损伤演化”和“是否考虑层间界面”两个维度分类：
| 类型 | 典型关键字 | 适用场景 |
|---|---|---|
| 线弹性正交各向异性 | *MAT_002（OPTION=ORTHO） | 只关心刚度与应力波传递，不关心失效 |
| 渐进损伤（面内） | *MAT_054/055（ENHANCED_COMPOSITE_DAMAGE） | 面内纤维/基体失效、冲击后剩余强度 |
| 连续损伤力学 | *MAT_058（LAMINATED_COMPOSITE_FABRIC） | 需要平滑刚度退化的织物或层压板 |
| 层间界面 | *MAT_138（COHESIVE_MIXED_MODE） | 分层起始与扩展，配合*CONTACT_..._TIED |
| 短纤维/各向同性近似 | *MAT_024、*MAT_003 | 快速估算或填充区域 |
选择原则：若研究目标是“是否分层”，必须引入*MAT_138内聚单元或*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_TIEBREAK；若只关注面内损伤，用*MAT_054即可。注意*MAT_054需要输入SOFT、SOFTG等软化参数，取值过大会导致过早删除单元，建议先用默认值试算再调整。
三、单元类型选择
复合材料层合板常用壳单元或实体单元：
- 壳单元：
*SECTION_SHELL配合ELFORM=2（Belytschko-Tsay）或ELFORM=16（全积分）。通过NIP设置厚度方向积分点，每个积分点对应一个铺层。例如4层板可设NIP=4，再用*INTEGRATION_SHELL指定各层角度与厚度。壳单元计算效率高，适合大尺寸板冲击。 - 实体单元：每层单独建一层实体，用
*SECTION_SOLID（ELFORM=1或-1）。能更真实地捕捉厚度方向应力与层间剪切，但网格量大，冲击接触区域需要细化。
教学示例：300 mm × 300 mm × 2 mm 的4层碳纤维层合板，铺层[0/90/90/0]，采用壳单元，NIP=4，面内网格尺寸2 mm，冲击区域局部加密到0.5 mm。若改用实体单元，厚度方向每层至少1个单元，总单元数约增加一个量级。
四、边界条件与初始速度设置
典型落锤冲击或弹丸冲击的边界设置如下：
- 夹具约束：用
*BOUNDARY_SPC_SET对板四周节点约束全部平动自由度，模拟环形夹具。若夹具直径为200 mm，则约束半径100 mm以外的节点。 - 对称边界：若几何与载荷对称，可用
*BOUNDARY_SPC_SET施加对称约束，减少计算量。 - 初始速度：用
*INITIAL_VELOCITY对冲击体节点集施加初速度。例如弹丸质量5 kg、速度3 m/s，方向沿板法向（-Z）。注意速度方向与坐标系一致，若模型Z轴向上，则VZ=-3。 - 接触：用
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE定义冲击体与板的接触，设置FS、FD摩擦系数（示例取0.2）。
五、求解控制参数
显式求解的核心控制参数包括：
- 终止时间：
*CONTROL_TERMINATION的ENDTIM。冲击问题通常取冲击体接触后2~5倍接触时间，示例取5 ms。 - 时间步：
*CONTROL_TIMESTEP的DTINIT可留空由程序计算，TSSFAC默认0.9，若出现负体积可降到0.6~0.8。 - 质量缩放：
*CONTROL_TIMESTEP中DT2MS为负值时启用质量缩放，可显著减少小单元导致的时间步过小，但需监控added mass不超过总质量的1%~5%。 - 沙漏控制：壳单元用
*CONTROL_HOURGLASS，IHQ=4（类型4）配合QM=0.05，抑制沙漏能。 - 能量控制：
*CONTROL_ENERGY中HGEN=2、RWEN=2，输出沙漏能与刚性墙能量，便于判断结果可靠性。 - 输出：
*DATABASE_BINARY_D3PLOT设置DT=0.05 ms，*DATABASE_ASCII输出接触力与能量曲线。
下面用一张信号链与参数关系图，把上述建模流程中的关键对象、输入与判断关系串联起来，便于对照检查。
六、建模检查清单
提交计算前，建议按以下顺序核对关键字文件：
- 铺层角度序列与厚度是否与设计一致，
*INTEGRATION_SHELL的积分点数量是否等于层数。 - 材料卡片中
SOFT、SOFTG、FS等失效参数是否合理，避免默认值导致过早删除。 - 接触卡片是否覆盖冲击体与板、板与夹具之间的所有可能接触面。
- 初始速度方向与坐标系是否一致，速度值是否与冲击能量匹配。
- 时间步与质量缩放是否导致
added mass超过5%，若超过需减小DT2MS或细化网格。 - 输出频率是否足够捕捉接触力峰值，
D3PLOT间隔建议不大于0.05 ms。
NIP与铺层数不一致，导致角度错位；或*MAT_054未设置SOFT，单元在低应力下被删除。建议先用线弹性材料试算，确认接触力与能量曲线正常后再替换为损伤模型。完成上述配置后，即可提交LS-DYNA求解。下一阶段将围绕冲击响应曲线、损伤云图和分层面积展开结果分析，本节建立的建模链路是后续所有判读工作的基础。