3. 复合材料层合板几何建模与网格划分

复合材料层合板的冲击仿真，成败往往在建模阶段就已决定大半。与各向同性金属不同，层合板具有方向性：面内纤维方向刚度高、横向与层间靠基体承载，厚度方向几乎没有纤维增强。因此“用什么单元、铺层方向怎么定、网格划多密”这三件事不是彼此独立的操作步骤，而是一组相互制约的建模决策。本节围绕壳单元与实体单元的选择、铺层方向与网格取向的对应关系、网格密度对冲击响应的影响与收敛性验证展开。

一、壳单元与实体单元的选择策略

壳单元（如 LS-DYNA 中的 *SECTION_SHELL 配合 *PART_COMPOSITE）把层合板沿厚度方向压缩为参考面，通过积分点定义各铺层，计算量小、适合大尺寸板件与低速冲击的全局响应。实体单元（*SECTION_SOLID，常用 ELFORM=1 或 -1）则显式建立每一层的几何厚度，能捕捉厚度方向应力梯度、层间剪切与分层前的应力集中，代价是单元数量与计算时间成倍增加。

选择的核心判据是：冲击问题中真正需要被解析的物理量是否沿厚度方向变化。若关注的是冲击力时程、板中心挠度、整体吸能，且冲头直径远大于板厚，壳单元通常足够；若关注分层起始位置、层间应力、基体开裂或需要与 cohesive 单元配合模拟层间失效，则必须使用实体单元，或采用“壳+实体”混合建模：冲击区局部用实体层合，外围用壳单元，通过 *CONSTRAINED_SHELL_TO_SOLID 或共节点过渡。

一个实用的经验规则是：当冲头直径与板厚之比大于约 10、且不关心层间失效时优先壳单元；当该比值小于 5、或明确要研究分层时优先实体单元。壳单元建模时，厚度方向积分点数量（NIP）应至少等于铺层数，否则无法区分各层应力；对多层合板，建议每个铺层至少 1 个积分点，并启用 *CONTROL_SHELL 中的 ISTUPD 等选项控制厚度更新。

二、铺层方向与网格取向的对应关系

这是复合材料建模中最容易被忽视、后果却最严重的一环。铺层角（如 0°、90°、±45°）是相对于单元坐标系定义的，而单元坐标系又由网格的节点编号顺序决定。若网格取向与设计意图不一致，0° 铺层可能实际变成了 90°，整个刚度矩阵随之错误。

在 LS-DYNA 中，壳单元的局部坐标系由单元节点顺序确定：节点 1→2 方向为局部 x 轴，法向为局部 z 轴。因此，必须保证网格的 1→2 方向与铺层定义中的 0° 参考方向一致。对于规则矩形板，通常让网格 1→2 方向沿板长方向，并在 *MAT 或 *PART_COMPOSITE 中按此基准定义各层角度。对于曲面或复杂外形，应通过 *DEFINE_COORDINATE_SYSTEM 或 *ELEMENT_SHELL 的节点顺序统一控制，必要时用 *INITIAL_STRESS_SHELL 或坐标系统卡片校正。

实体单元没有“单元坐标系角度”这一直接概念，铺层方向需通过材料坐标系（*MAT_ORTHOTROPIC 配合 *DEFINE_COORDINATE_SYSTEM）或 *PART_COMPOSITE 的铺层方向向量来指定。此时网格取向影响的是材料主轴的插值精度，规则六面体网格配合明确的材料坐标系最为可靠。

三、网格密度对冲击响应的影响与收敛性验证

冲击仿真中，网格密度直接决定能否捕捉接触区的高应力梯度与损伤演化。冲头与板接触区存在剧烈的应力集中，若网格过粗，接触力峰值被平滑、损伤起始被延迟，甚至出现“穿模”或能量异常。一般建议接触区最小单元尺寸不超过冲头半径的 1/5～1/10，且沿厚度方向实体单元至少每层 1 个，关键层可加密到 2～3 个。

但网格并非越密越好。显式动力学中，最小单元尺寸决定临界时间步长 Δt ≈ L_min/c（c 为材料波速），网格加密会使时间步长成比例减小，计算时间急剧上升。因此需要在精度与成本之间寻找收敛点。

收敛性验证的标准做法是：固定其他条件，逐级加密网格（如接触区单元尺寸取 2.0、1.5、1.0、0.75 mm），提取冲击力峰值、最大挠度、吸能量、分层面积等关键响应量，观察其随网格尺寸的变化。当相邻两级网格的结果差异小于 5%（工程上常取 3%～5%）时，可认为已收敛。下表给出一个示例性的收敛验证记录（数值为教学假设，非实测）：

从示例数据可见，1.0 mm 到 0.75 mm 时冲击力峰值仅变化约 2%，可认为该工况下 1.0 mm 已接近收敛。实际项目中还应检查能量比（总能量/初始动能）、沙漏能占比（一般应小于 5%）以及接触能是否合理，作为收敛判据的补充。

综上，壳与实体的选择取决于是否需要解析厚度方向物理量，铺层方向必须与网格取向严格绑定，网格密度则通过分级加密与关键响应量对比来确定收敛点。三者共同构成层合板冲击建模的质量基础，任何一环出错都会使后续的损伤分析与实验对标失去意义。