2. 单层材料本构与工程常数

复合材料层合板的宏观力学行为，最终都要追溯到每一个单层的应力—应变关系。单层是层合结构的最小力学单元，只有先把单层的本构关系、工程弹性常数以及材料主方向与整体坐标之间的变换关系弄清楚，后续的层合板刚度合成、铺层应力计算才有可靠基础。本节围绕正交各向异性单层的本构方程展开，重点说明柔度矩阵与刚度矩阵的构成、工程弹性常数的物理含义，以及坐标变换在有限元建模中的实际用法。

一、正交各向异性单层的应力—应变关系

连续纤维单层在宏观上可视为正交各向异性材料：沿纤维方向（记为材料主方向 1）、垂直纤维的面内方向（记为 2）以及厚度方向（记为 3）各有一个弹性对称面。对于平面应力状态下的单层（厚度方向应力 σ₃、τ₂₃、τ₁₃ 近似为零），面内应力—应变关系可写成：

σ₁ = Q₁₁ ε₁ + Q₁₂ ε₂
σ₂ = Q₁₂ ε₁ + Q₂₂ ε₂
τ₁₂ = Q₆₆ γ₁₂

其中 Q 为平面应力刚度矩阵分量，ε₁、ε₂ 为正应变，γ₁₂ 为工程剪应变。写成矩阵形式即 {σ} = [Q]{ε}。注意工程剪应变 γ₁₂ 与张量剪应变 ε₁₂ 相差因子 2，这是有限元实现中极易出错的地方：多数商业软件内部使用张量剪应变，而材料手册给出的往往是工程剪应变，输入时必须确认约定一致。

反过来，用柔度矩阵表示应变为 {ε} = [S]{σ}，平面应力下柔度矩阵为：

S₁₁ = 1/E₁
S₂₂ = 1/E₂
S₁₂ = -ν₁₂/E₁ = -ν₂₁/E₂
S₆₆ = 1/G₁₂

柔度矩阵是对称的，因此 ν₁₂/E₁ = ν₂₁/E₂，即两个主泊松比并不独立，这一互等关系是校核材料数据是否自洽的基本依据。

二、工程弹性常数的含义与量级

单层正交各向异性材料在平面应力下由四个独立工程常数描述：纵向弹性模量 E₁、横向弹性模量 E₂、主泊松比 ν₁₂ 和面内剪切模量 G₁₂。以典型碳纤维/环氧单层为例（示例参数，仅用于教学演示）：E₁ ≈ 135 GPa，E₂ ≈ 9 GPa，G₁₂ ≈ 5 GPa，ν₁₂ ≈ 0.30。可以看到 E₁ 比 E₂ 高一个数量级以上，这种强烈的方向性正是复合材料可设计性的来源，也是层合板中铺层角度直接影响刚度的根本原因。

由这些常数可算出平面应力刚度分量：Q₁₁ = E₁/(1−ν₁₂ν₂₁)，Q₂₂ = E₂/(1−ν₁₂ν₂₁)，Q₁₂ = ν₁₂E₂/(1−ν₁₂ν₂₁)，Q₆₆ = G₁₂。代入上述示例值，ν₂₁ = ν₁₂E₂/E₁ ≈ 0.02，分母 1−ν₁₂ν₂₁ ≈ 0.994，可得 Q₁₁ ≈ 135.8 GPa、Q₂₂ ≈ 9.05 GPa、Q₁₂ ≈ 2.72 GPa、Q₆₆ = 5 GPa。这些数值将直接进入层合板 ABD 矩阵的积分计算。

三、材料主方向与坐标变换

有限元模型的整体坐标系（x, y）通常与单层的材料主方向（1, 2）不一致，铺层角 θ 定义为材料主方向 1 与整体 x 轴的夹角，逆时针为正。此时需要把材料主方向下的刚度矩阵 [Q] 变换到整体坐标下的 [Q̄]：

[Q̄] = [T]⁻ᵀ [Q] [T]⁻¹

其中 [T] 为应力变换矩阵，与 θ 的余弦、正弦有关。工程上更常用展开后的显式公式，例如：

Q̄₁₁ = Q₁₁cos⁴θ + 2(Q₁₂+2Q₆₆)sin²θcos²θ + Q₂₂sin⁴θ
Q̄₁₂ = (Q₁₁+Q₂₂−4Q₆₆)sin²θcos²θ + Q₁₂(sin⁴θ+cos⁴θ)
Q̄₆₆ = (Q₁₁+Q₂₂−2Q₁₂−2Q₆₆)sin²θcos²θ + Q₆₆(sin⁴θ+cos⁴θ)
Q̄₁₆ = (Q₁₁−Q₁₂−2Q₆₆)sinθcos³θ − (Q₂₂−Q₁₂−2Q₆₆)sin³θcosθ
Q̄₂₆ = (Q₁₁−Q₁₂−2Q₆₆)sin³θcosθ − (Q₂₂−Q₁₂−2Q₆₆)sinθcos³θ

变换后出现 Q̄₁₆、Q̄₂₆ 两个耦合项，它们把正应力与剪应变、剪应力与正应变联系起来，即拉伸—剪切耦合。当 θ = 0° 或 90° 时耦合项为零；θ = ±45° 时耦合最强。这正是 ±45° 铺层能显著提高层合板剪切刚度和改善耦合效应的原因。

四、从单层到有限元输入的实操要点

实际建模时，单层材料数据通常以工程常数形式输入，软件内部自动完成 [Q] 与 [Q̄] 的计算。需要重点检查三件事：其一，E₁、E₂、G₁₂ 的单位与模型长度单位一致；其二，ν₁₂ 与 ν₂₁ 的互等关系是否满足，若手册只给一个泊松比，应确认软件是否自动换算；其三，剪切模量 G₁₂ 是否与工程剪应变约定匹配。对于三维实体单元，还需补充 E₃、G₁₃、G₂₃、ν₁₃、ν₂₃ 等面外常数，此时材料不再是平面应力状态，柔度矩阵为 6×6 满阵。

下面用一张图概括单层本构的核心关系：从工程常数出发，经柔度矩阵与刚度矩阵，再经坐标变换得到整体坐标下的刚度，最终进入层合板分析。

理解这张关系链后，有限元建模中的材料卡片就不再是一堆孤立数字：E₁、E₂、G₁₂、ν₁₂ 决定 [Q]，铺层角 θ 决定 [Q̄]，而 [Q̄] 沿厚度积分得到层合板的 ABD 矩阵。当计算结果出现异常耦合或刚度明显偏离预期时，应优先回到本节检查工程常数单位、泊松比互等关系以及铺层角正负号约定，这三处是单层本构环节最常见的错误来源。