3. 有限元方法在复合材料中的实现要点
复合材料层合结构在有限元中的实现，与各向同性金属件有本质差别：材料刚度随铺层方向变化，面内与面外响应尺度差异大，层间应力往往决定失效。因此，从位移场插值到收敛判断，每一步都需要围绕“方向性”和“分层性”做专门处理。本节按建模链条梳理六个关键实现要点，并给出可操作的参数与判断依据。
位移场插值：先确定单元阶次与积分策略
位移场插值决定了单元能表达什么样的变形模式。对于以弯曲为主的层合板壳，若采用线性插值的四节点壳单元，在厚度方向只有一层积分点时，面外弯曲刚度容易被高估，出现“剪切自锁”和“沙漏”现象。工程上常用两种做法：一是采用二次插值的八节点壳单元（如 S8R 类），配合减缩积分；二是对线性单元使用增强假设应变（EAS）或选择性减缩积分，使弯曲与剪切项分别积分。
以一块 200 mm × 100 mm、厚 2 mm 的碳纤维层合板为例，若用四节点壳单元，面内网格取 2 mm 时约 5000 个单元；改用八节点单元后，同样精度下网格可放宽到 4 mm，单元数降至约 1250，计算量明显下降。判断插值是否足够，可观察厚度方向弯曲应力沿板厚的分布是否光滑；若出现锯齿状振荡，通常说明积分点不足或单元阶次偏低。
壳单元与实体单元的选择
壳单元把三维问题降为二维，计算效率高，适合厚度远小于面内尺寸（一般厚度/边长 < 1/10）的薄壁层合板。实体单元保留厚度方向的正应力与层间剪应力，适合厚截面、自由边附近、连接区以及需要提取层间应力的部位。实际建模中常采用混合策略：主体用壳单元，螺栓孔、胶接端部等局部区域用实体单元，两者通过壳-实体耦合或绑定约束连接。
选择时可用一个简单判据：若关注的是面内应力与整体刚度，壳单元足够；若关注分层起始、层间剪切或厚度方向压缩，必须用实体单元，且厚度方向至少布置 2～4 个单元。对于 2 mm 厚层合板，若用实体单元，厚度方向单元尺寸约 0.5 mm，单层 0.125 mm 的铺层无法逐层建模，此时应改用分层壳或等效实体。
层合结构建模策略
层合板建模有三种常见策略。第一种是“分层壳”：在壳截面中定义每一铺层的材料、厚度和方向角，单元仍为壳，但沿厚度积分时逐层计算应力。第二种是“等效单层”：用经典层合理论算出等效 ABD 刚度矩阵，再赋给壳单元，适合只关心整体变形、不关心层间应力的场景。第三种是“逐层实体”：每层单独建一层实体单元，精度最高但计算量最大。
以 [0/90/±45]s 共 8 层、单层 0.125 mm 的层合板为例，分层壳只需在截面定义中写入 8 行铺层信息，单元数不变；逐层实体则厚度方向至少 8 层单元，单元数增加约 8 倍。若只做整体刚度校核，等效单层即可；若要判断 90° 层是否先出现横向裂纹，必须用分层壳或逐层实体。
材料方向定义
复合材料刚度矩阵必须建立在材料主方向上，再通过方向角转换到整体坐标。壳单元中，材料 1 方向通常沿纤维方向，2 方向为面内垂直纤维方向，3 方向为厚度方向。定义时需明确：铺层角是相对于整体坐标 x 轴还是参考轴，旋转正负号按右手定则。若定义混乱，0° 层与 90° 层会互换，计算结果完全失真。
实操中建议在每层属性中显式写出方向角，并在后处理中检查 0° 层的纤维方向应力是否明显大于 90° 层。例如 [0/90]s 板受 x 向拉伸时，0° 层 σ11 应远大于 90° 层 σ11；若两者接近，说明方向定义有误。对于实体单元，还需为每个单元指定局部材料坐标系，常用方法是通过单元方向或参考坐标系赋值。
边界条件与载荷施加
边界条件要尽量模拟真实约束，避免过约束或欠约束。简支板可在四边施加法向位移约束，同时释放面内转动；固支板则约束全部平动与转动。载荷施加时，集中力应加在节点或耦合参考点上，分布压力加在单元面上，螺栓载荷可通过耦合到孔壁节点实现。
以三点弯曲为例，支座处约束竖向位移与面外位移，压头处施加位移控制。若采用位移加载，可避免收敛困难；若采用力加载，需注意初始增量步不宜过大。对于层合板，还应注意载荷方向与材料方向的夹角，必要时在加载点建立局部坐标系。
收敛性判断
收敛性判断包括网格收敛与非线性迭代收敛两方面。网格收敛可通过逐步加密网格，观察关键应力或位移的变化；当两次网格结果差异小于 5% 时，可认为网格收敛。非线性问题中，需检查残差是否低于容差、增量步是否完成、是否出现负特征值。
若出现不收敛，常见原因有：材料软化导致刚度矩阵奇异、接触状态突变、载荷步过大。可尝试减小初始增量、开启自动稳定、改用位移控制或调整接触刚度。对于分层损伤分析，还应检查损伤变量是否突变，避免因单元删除导致能量不守恒。
把上述六点串起来，实际建模流程可以概括为：先根据厚度与关注量选择壳或实体，再确定插值阶次与积分方案；随后在截面或单元属性中逐层定义材料方向与铺层角；接着施加符合物理实际的边界与载荷；最后通过网格加密和迭代残差判断结果是否可信。每一步都直接影响后续，方向定义错误会让应力分布完全颠倒，单元选择不当则可能掩盖层间应力。建议在提交计算前，先用一个简单铺层（如 [0/90]s）做小模型验证，确认方向、约束和收敛行为符合预期，再扩展到完整结构。