2. 正交各向异性材料本构与工程常数
复合材料单层板在宏观尺度上通常被理想化为正交各向异性材料：存在三个相互垂直的弹性对称面，沿这三个材料主轴方向的力学行为彼此独立。理解本构关系的关键，是把刚度矩阵、柔度矩阵与工程常数三者之间的对应关系建立起来，并能在 Abaqus 中正确输入与验证。
一、三维正交各向异性本构关系
在材料主轴坐标系（1、2、3）下，应力与应变通过广义胡克定律联系。用柔度矩阵 S 表示更直观，因为柔度分量直接由工程常数构成：
ε = S · σ
S = [ 1/E1 -ν12/E1 -ν13/E1 0 0 0 ]
[ -ν12/E1 1/E2 -ν23/E2 0 0 0 ]
[ -ν13/E1 -ν23/E2 1/E3 0 0 0 ]
[ 0 0 0 1/G12 0 0 ]
[ 0 0 0 0 1/G23 0 ]
[ 0 0 0 0 0 1/G13 ]其中 E1、E2、E3 为三个主轴方向的杨氏模量，G12、G23、G13 为三个主平面内的剪切模量，ν12、ν13、ν23 为泊松比。注意柔度矩阵是对称的，因此泊松比满足互等关系：ν12/E1 = ν21/E2，即 ν21 = ν12·E2/E1。这一点在输入 Abaqus 时极易出错——Abaqus 要求的是 ν12、ν13、ν23（即“主方向对次方向”的泊松比），而不是反过来。
刚度矩阵 C 是柔度矩阵的逆：σ = C · ε，C = S⁻¹。工程上更常用工程常数，因为 E、G、ν 可直接由单层板拉伸、剪切试验测得，而 C 的分量是这些常数的复杂组合。对于平面应力状态（薄板，σ3 = τ23 = τ13 = 0），本构退化为：
σ1 1 ν21 0 ε1
σ2 = ────── · ν12 1 0 · ε2
τ12 1-ν12ν21 0 0 (1-ν12ν21) γ12
· G12这就是经典层合板理论中单层板的平面应力本构，也是 Abaqus 壳单元与连续壳单元中最常用的形式。
二、工程常数的物理含义与典型量级
以常见碳纤维/环氧单层板（示例参数，非实测值）为例：E1 ≈ 135 GPa，E2 ≈ 9 GPa，E3 ≈ 9 GPa，G12 ≈ 5 GPa，G23 ≈ 3 GPa，G13 ≈ 5 GPa，ν12 ≈ 0.30，ν23 ≈ 0.40，ν13 ≈ 0.30。可以看到 E1 远大于 E2，体现纤维方向承载能力远高于横向；G12 与 G23 差异说明面内剪切与层间剪切刚度不同。这些数值关系决定了单层板强烈的方向依赖性。
三、Abaqus 中 Engineering Constants 的定义
在 Abaqus/CAE 中，材料模块选择 Mechanical → Elasticity → Elastic → Type: Engineering Constants，可输入三维正交各向异性常数。若使用壳单元，则选择 Type: Lamina，只需输入 E1、E2、ν12、G12、G13、G23（平面应力假设下 E3、ν13、ν23 不参与面内计算，但 G13、G23 影响横向剪切）。
输入时需注意：材料方向必须与单元局部坐标系一致。壳单元默认材料 1 方向沿单元局部 1 轴，需通过 Assign Material Orientation 指定铺层角度；实体单元则需定义局部坐标系或使用 Orientation 赋予材料方向。若方向定义错误，即使常数正确，结果也会完全失真。
四、验证方法
最直接的验证是单单元试验：建立一个 1×1×1 的 C3D8 单元，沿 1 方向施加单位应变，检查反力是否等于 E1；沿 12 面施加单位剪应变，检查剪应力是否等于 G12。对壳单元，可建立单层板拉伸模型，将 Abaqus 计算的面内模量与经典层合板理论解析值对比。若偏差超过 1%，应检查材料方向、单元类型与边界条件。
验证时还应检查柔度矩阵的正定性。可在 Abaqus 中提交一个仅含材料定义的小模型，若材料常数不满足正定条件，求解器会在 .dat 文件中给出警告或错误。对于层合板分析，建议先用单层板验证，再逐步增加铺层，避免方向与常数错误被层间耦合掩盖。
掌握工程常数与本构矩阵的对应关系后，后续铺层设计、失效判据与层间应力分析才有可靠基础。建议在正式分析前，始终用单单元或单层板模型做一次快速验证，确认材料方向、常数输入与单元类型三者一致。