4. 层合结构有限元建模总体流程

层合结构有限元建模与各向同性金属件建模的最大区别在于：材料方向随铺层角变化、厚度方向存在逐层刚度差异、层间界面可能成为失效起点。因此，建模流程不能简单套用“画几何—赋材料—划网格—加约束”的通用顺序，而要把铺层信息、单元类型和坐标系约定作为贯穿始终的主线。本节按几何简化、材料定义、铺层设置、单元选择、边界条件、载荷施加与结果后处理七个环节，说明各环节的关键输入、操作要点和相互制约关系。

几何简化：先决定“用什么几何承载铺层”

层合板在有限元中通常简化为壳、连续壳或实体三类几何。壳几何只需中面，厚度和铺层由截面属性承载；连续壳用一层实体单元代表一层或若干层，厚度方向至少一个单元；实体几何则显式建出每层厚度。选择依据是关注量：若只关心面内应力与整体变形，壳足够；若关心层间剪应力、自由边效应或分层，必须用连续壳或实体。

几何简化时要去掉对刚度贡献很小的倒角、小孔和工艺凸台，但孔边应力集中区、连接区和自由边不能随意简化。一个常见做法是：整体模型用壳，局部危险区用实体子模型，通过位移边界传递。简化后要检查中面位置是否与铺层基准面一致，否则后续铺层偏移会引入附加弯矩。

材料定义：单层工程常数与强度参数

层合结构的基本材料是单层（lamina），其面内工程常数包括纵向模量 E1、横向模量 E2、主泊松比 ν12、面内剪切模量 G12，以及层间模量 E3、G13、G23。以典型碳纤维/环氧单层为例（示例假设值）：E1=135 GPa，E2=9.5 GPa，E3=9.0 GPa，ν12=0.30，G12=5.0 GPa，G13=4.8 GPa，G23=3.2 GPa。这些值必须与材料坐标系绑定，1 方向为纤维方向。

强度参数至少需要 Xt、Xc、Yt、Yc、S12，若做层间分析还需 S13、S23 和层间断裂韧性。材料定义阶段要明确单位制并全文统一，例如统一用 mm、N、MPa，避免模量用 GPa 而几何用 mm 导致量级错误。

铺层设置：角度、顺序、厚度与坐标系

铺层设置是把单层材料“堆”成层合板的过程，需要给出每层材料、厚度、铺层角和参考方向。以 [0/45/-45/90]s 为例，共 8 层，若单层厚 0.125 mm，总厚 1.0 mm。铺层角定义为纤维方向与参考轴（通常为局部 1 轴）的夹角，逆时针为正。参考轴必须与几何边或全局坐标明确关联，否则同一串角度在不同参考下结果完全不同。

在壳单元中，铺层通过截面属性赋予，需指定积分点数量。若每层至少 3 个积分点，8 层至少 24 个，才能较好捕捉层内弯曲应力梯度。连续壳和实体则直接按层划分单元。铺层设置完成后应输出一张铺层表，核对总厚度、对称性和角度顺序。

单元选择：壳、连续壳与实体的分工

壳单元计算效率高，适合大尺寸层合板整体分析，但厚度方向应力需后处理近似。连续壳单元在厚度方向布置一层或多层实体单元，能直接给出层间应力，代价是网格量增加。实体单元最通用，但层合板厚度方向尺寸远小于面内尺寸，若每层一个单元，长厚比可能达到数百，容易导致剪切自锁，需要控制面内尺寸或使用减缩积分。

选择原则可归纳为：整体刚度与面内应力用壳；层间应力与自由边用连续壳或实体；冲击、接触、分层扩展用实体。同一模型可混合使用，但壳与实体连接处要处理转动自由度传递，常用刚性连接或壳到实体耦合。

边界条件与载荷施加

边界条件要反映真实约束，而不是“随便固定几个点”。例如四点弯曲示例中，下支撑为简支，上加载为线载荷；若把支撑点全固定，会人为限制转动，导致跨中弯矩偏大。约束不足会产生刚体位移，约束过多会改变内力路径，两者都要避免。

载荷施加需注意层合结构的耦合效应：拉伸载荷可能引起弯扭耦合（非对称铺层），因此载荷作用点应尽量与中面一致，或明确偏心量。分布载荷要按实际作用面施加，集中力要扩散到足够节点，避免单点应力奇异。多工况分析时，每个工况单独定义载荷步，后处理再组合。

结果后处理：逐层看，而不是只看整体

层合结构后处理的核心是“逐层”。壳单元要输出每层积分点的应力，检查 σ1、σ2、τ12 是否超过对应强度；连续壳和实体还要看 σ3、τ13、τ23。以示例 [0/45/-45/90]s 板在面内拉伸下，0°层 σ1 最高，90°层 σ2 较高，45°层 τ12 较大，各层失效模式不同，不能只看最大位移。

判断依据可用失效因子：例如最大应力准则下，失效因子 = max(σ1/Xt, σ2/Yt, |τ12|/S12)。若某层失效因子大于 1，说明该层已失效，但结构未必整体破坏，需结合刚度退化继续分析。后处理还应检查层间应力在自由边是否集中，必要时细化网格验证收敛性。

建模检查清单（可直接对照使用）

把上述七个环节串起来，层合结构有限元建模的实质是：用几何承载铺层，用材料定义单层，用铺层设置决定方向与顺序，用单元选择决定厚度方向分辨率，用边界与载荷还原工况，最后用逐层后处理判断失效。每一步的输入都来自上一步，任何一步的近似都会传递到结果，因此建模过程中保留完整的铺层表、材料卡和工况记录，比事后补记更可靠。