4. Abaqus复合材料单元类型选择
复合材料层合板在厚度方向上的力学行为具有明显的分层特征：面内纤维主导的刚度远高于基体主导的层间刚度，层间剪切与正应力又直接决定分层、基体开裂等失效模式。因此，在 Abaqus 中选择单元类型，本质上是在回答一个问题——你希望求解器如何离散厚度方向，以及你愿意为这种离散付出多少计算代价。常规壳、连续壳与实体单元，正是三种不同取舍的答案。
三类单元的力学假设与适用边界
常规壳单元（S4R、S8R、SC8R 之外的 S 系列）基于 Kirchhoff–Love 或 Mindlin–Reissner 板壳理论，把三维问题降为二维中面问题。它通过截面定义（Section）赋予厚度，厚度方向的应力由本构关系与平面应力假设推导。其优点是自由度少、计算快，适合大尺寸、薄壁、以面内应力为主的层合板整体刚度与屈曲分析。但它默认厚度方向应力 σ33 与层间剪切由后处理近似给出，无法直接捕捉自由边层间应力集中，也难以模拟分层扩展。
连续壳单元（SC6R、SC8R）是 Abaqus 中专门为复合材料设计的“降维实体”。它保留三维本构与厚度方向积分点，但几何上仍以中面描述，节点自由度与壳类似。相比常规壳，它能更准确地给出厚度方向正应力与层间剪切，且支持多层截面定义，是层合板强度校核与首层失效判据的常用选择。代价是计算量高于常规壳，且对网格质量与厚度方向积分点数量更敏感。
实体单元（C3D8R、C3D8、C3D20R 等）完全放弃板壳假设，直接三维离散。它能捕捉自由边效应、层间应力、分层与压溃，是精细失效分析、冲击与接触问题的首选。但厚度方向必须划分足够多的单元，层数多时网格规模急剧膨胀，计算成本最高。
沙漏控制：减缩积分单元的“副作用”
S4R、SC8R、C3D8R 中的“R”代表减缩积分（Reduced Integration）。减缩积分在每个方向少用一个积分点，显著降低计算量并缓解剪切自锁，但会引入沙漏模式（Hourglassing）——一种零能量变形模式，单元可以发生锯齿状变形而不产生应变能，表现为网格出现棋盘状、无物理意义的振荡。
判断沙漏是否严重，最直接的指标是伪应变能（ALLAE）与内能（ALLIE）之比。工程经验上，该比值应控制在 5% 以内；若超过 10%，结果通常不可信。控制手段包括：
- 增强沙漏控制：在单元控制中选用 Enhanced 或 Relax Stiffness，适用于大变形与接触问题；
- 细化网格：沙漏模式波长与单元尺寸相关，网格越细越不易激发；
- 改用完全积分或非协调单元：如 C3D8I，但需注意剪切自锁与计算量上升；
- 施加合理边界与载荷：避免单点集中载荷直接作用于减缩积分单元。
厚度方向积分点：复合材料精度的关键旋钮
对连续壳与实体单元，厚度方向积分点数量直接决定层间应力能否被解析。Abaqus 中通过截面定义设置：常规壳与连续壳在 Shell Section 中指定 Integration Points，实体单元则通过 Section Points 或分层材料指派控制。
基本规则是：每个铺层至少布置 1 个积分点，层间应力梯度大的区域建议 2～3 个。若一个截面有 8 层，积分点取 8 可保证每层一个；若关注自由边层间剪切，可在靠近自由边的局部区域加密至 16 或 24。积分点过少会导致层间应力被平均化，首层失效判据失真；过多则计算量线性上升，且可能引入数值噪声。
*Shell Section, elset=Laminate, composite
0.125, 3, Ply-0, 0.0
0.125, 3, Ply-90, 90.0
0.125, 3, Ply-45, 45.0
0.125, 3, Ply-45, -45.0
0.125, 3, Ply-90, 90.0
0.125, 3, Ply-0, 0.0
0.125, 3, Ply-0, 0.0
0.125, 3, Ply-90, 90.0
上述示例中每层设置 3 个积分点，共 24 个厚度积分点，适用于需要输出层间应力的连续壳模型。若仅做整体屈曲分析，可降为每层 1 点以节省时间。
实操建议与常见误区
在具体建模时，建议先做一次“单元敏感性试算”：同一层合板分别用 S4R、SC8R 和 C3D8R 建模，比较中心挠度与首层失效载荷。若三者差异在 5% 以内，说明整体响应由面内刚度控制，可放心使用常规壳；若连续壳与实体差异显著，则说明层间效应不可忽略，应优先连续壳或实体。
常见误区有三：其一，用常规壳做自由边分层分析，结果必然偏乐观；其二，连续壳厚度积分点只设 1 点，层间应力被平均，失效判据失效；其三，看到沙漏振荡就盲目加密网格，却未检查 ALLAE 比值，可能浪费大量计算资源。正确的顺序是：先看能量比，再决定是否加密或更换沙漏控制方式。
Section Points 中显式指定，并在后处理中检查每层应力是否连续。最后需要强调：单元类型与积分点设置不是孤立的，它们与网格密度、沙漏控制、材料本构共同决定结果可信度。一个可靠的复合材料分析，应当在提交计算前就明确“厚度方向如何离散、沙漏如何抑制、失效如何判读”这三个问题，而不是等结果异常后再回头猜测。