2. 层合板理论快速回顾与工程简化

在ANSYS复合材料结构仿真中，真正决定分析结果是否可信的，往往不是网格有多密、求解器设置有多精细，而是对层合板力学行为的理解是否到位。经典层合板理论（Classical Laminate Theory, CLT）是连接单层材料属性与整体层合板响应的桥梁，而ABD矩阵则是这座桥梁的数学核心。本节将快速回顾CLT的基本假设与推导逻辑，重点厘清ABD矩阵各子块的物理意义，并说明在ANSYS中如何将单层正轴/偏轴刚度、工程常数输入与失效判据选择串联成一条可执行的工程链路。

一、经典层合板理论的核心假设与位移场

CLT建立在克希霍夫（Kirchhoff）薄板假设之上，其核心可归纳为三点：第一，直法线假设——变形前垂直于中面的直线，变形后仍保持直线且垂直于变形后的中面；第二，厚度方向应变忽略，即ε_z≈0；第三，各单层之间理想粘接，无滑移与脱层。由此，层合板内任意一点的位移可由中面位移u₀、v₀、w₀及转角表示，最终得到应变沿厚度方向的线性分布关系：

ε(z) = ε⁰ + z·κ

其中ε⁰为中面应变，κ为中面曲率。这一关系是后续刚度积分的起点，也是理解ABD矩阵中“A管面内、D管弯曲、B管耦合”的物理基础。

二、单层正轴刚度与偏轴转换

单层材料在正轴（1-2坐标系，1为纤维方向）下的刚度由四个独立工程常数决定：E₁、E₂、ν₁₂、G₁₂。由此可计算正轴刚度矩阵Q的各分量，例如：

Q11 = E1 / (1 - ν12·ν21)
Q22 = E2 / (1 - ν12·ν21)
Q12 = ν12·E2 / (1 - ν12·ν21)
Q66 = G12
ν21 = ν12·E2 / E1

当单层铺放角为θ时，需通过坐标变换得到偏轴刚度矩阵Q̄。变换中会出现Q̄16、Q̄26等耦合项，它们正是导致拉伸-剪切耦合的根源。在ANSYS的ACP模块中，用户只需输入E₁、E₂、ν₁₂、G₁₂及铺层角，软件会自动完成Q到Q̄的转换，但理解这一过程有助于判断输入参数是否合理。

三、ABD矩阵的物理意义与工程简化

将各单层的偏轴刚度沿厚度积分，得到层合板的本构关系：

[N]   [A  B] [ε⁰]
[M] = [B  D] [κ ]

其中A为面内刚度矩阵，D为弯曲刚度矩阵，B为耦合刚度矩阵。A_ij反映层合板抵抗面内拉压与剪切的能力；D_ij反映抵抗弯曲与扭转的能力；B_ij则描述面内变形与弯曲变形之间的耦合。对于对称铺层，B=0；对于正交铺设，A₁₆=A₂₆=0；对于准各向同性铺层，A矩阵呈现各向同性特征。工程中常通过对称、均衡、正交等铺层设计来消除或利用这些耦合。

四、ANSYS中的工程常数输入与失效判据选择

在ANSYS Composite PrepPost（ACP）中，材料定义通常采用Engineering Data中的Orthotropic Elasticity，输入E₁、E₂、E₃、ν₁₂、ν₁₃、ν₂₃、G₁₂、G₁₃、G₂₃。对于二维壳单元，E₃、ν₁₃、ν₂₃、G₁₃、G₂₃对面内响应影响较小，但仍需合理赋值以避免刚度矩阵奇异。铺层定义时，每层需指定材料、厚度和角度，ACP会自动计算ABD矩阵并生成等效截面属性。

失效判据的选择直接决定安全裕度的计算方式。常用判据包括最大应力/最大应变判据、Tsai-Hill、Tsai-Wu、Hashin及Puck等。最大应力判据简单直观，适用于初步筛选；Tsai-Wu为二次交互判据，适合纤维主导的拉压组合；Hashin区分纤维与基体失效模式，更适合工程失效分析。在ANSYS中，可通过FC（Failure Criteria）命令或ACP后处理中的Failure工具指定判据，并输出各层的失效因子。选择时应结合材料类型、载荷工况与设计阶段：概念设计可用Tsai-Wu快速评估，详细设计则推荐Hashin或Puck以区分失效模式。

五、从ABD矩阵到仿真结果的判断逻辑

完成ABD矩阵构建后，可通过几个快速检查判断模型是否合理：若B矩阵非零且量级与A、D相当，说明铺层不对称，需关注固化变形或拉弯耦合是否被正确考虑；若A₁₆、A₂₆显著非零，说明存在拉剪耦合，面内载荷下可能出现 unexpected 的剪切变形；若D矩阵中D₁₆、D₂₆较大，则弯曲-扭转耦合明显，需检查边界条件是否约束了扭转自由度。这些判断无需运行完整求解，仅通过ABD矩阵即可快速评估铺层设计的力学合理性。

综上，CLT与ABD矩阵是复合材料仿真的“底层语言”。理解单层正轴/偏轴刚度的转换逻辑，掌握ABD矩阵各子块的物理含义，并在ANSYS中正确输入工程常数、合理选择失效判据，才能让后续的应力分析、失效评估与优化设计建立在可靠的基础之上。