2. 单层复合材料力学基础
要理解层合板的整体响应，必须先弄清一块单层板在材料主轴方向上的本构关系。单层板由纤维与基体复合而成，宏观上可视为均匀的正交各向异性材料：沿纤维方向（记为 1 轴）刚度最高，垂直纤维方向（2 轴）刚度最低，面内剪切（1-2 平面）由基体主导。三个材料主方向互相垂直，因此独立的弹性常数只有 4 个：纵向模量 E₁、横向模量 E₂、主泊松比 ν₁₂ 和面内剪切模量 G₁₂。
一、工程常数与柔度矩阵
在材料主轴坐标系（1-2）下，平面应力状态（σ₃ = τ₂₃ = τ₃₁ = 0）的应力—应变关系可写成柔度形式：
ε₁ = S₁₁σ₁ + S₁₂σ₂
ε₂ = S₁₂σ₁ + S₂₂σ₂
γ₁₂ = S₆₆τ₁₂其中柔度分量与工程常数的换算关系为：
S₁₁ = 1/E₁
S₂₂ = 1/E₂
S₁₂ = -ν₁₂/E₁ = -ν₂₁/E₂
S₆₆ = 1/G₁₂注意 S₁₂ 的对称性要求 ν₁₂/E₁ = ν₂₁/E₂，即两个主泊松比并不独立，知道 ν₁₂ 后可由 ν₂₁ = ν₁₂·E₂/E₁ 推出。柔度矩阵 S 是对称矩阵，物理含义是“单位应力产生多少应变”，因此 E₁ 越大，S₁₁ 越小，材料越难被拉伸。
以典型碳纤维/环氧单层板为例（示例参数，仅用于教学演算）：E₁ = 135 GPa，E₂ = 9.0 GPa，ν₁₂ = 0.30，G₁₂ = 5.0 GPa。代入可得：
S₁₁ = 1/135 = 7.407×10⁻³ GPa⁻¹
S₂₂ = 1/9.0 = 1.111×10⁻¹ GPa⁻¹
S₁₂ = -0.30/135 = -2.222×10⁻³ GPa⁻¹
S₆₆ = 1/5.0 = 2.000×10⁻¹ GPa⁻¹可以看到 S₂₂ 比 S₁₁ 大约 15 倍，说明横向柔顺得多；S₁₂ 为负，表示纵向拉伸会引起横向收缩。
二、刚度矩阵与应力—应变关系
工程计算中更常用刚度形式，即用应变求应力：
σ₁ = Q₁₁ε₁ + Q₁₂ε₂
σ₂ = Q₁₂ε₁ + Q₂₂ε₂
τ₁₂ = Q₆₆γ₁₂刚度分量 Q 与工程常数的关系为：
Q₁₁ = E₁/(1 - ν₁₂ν₂₁)
Q₂₂ = E₂/(1 - ν₁₂ν₂₁)
Q₁₂ = ν₁₂E₂/(1 - ν₁₂ν₂₁)
Q₆₆ = G₁₂分母 Δ = 1 - ν₁₂ν₂₁ 称为正交各向异性耦合因子。对上述示例，ν₂₁ = 0.30×9.0/135 = 0.020，Δ = 1 - 0.30×0.020 = 0.994，接近 1，说明面内泊松耦合较弱。计算得 Q₁₁ ≈ 135.8 GPa，Q₂₂ ≈ 9.05 GPa，Q₁₂ ≈ 2.72 GPa，Q₆₆ = 5.0 GPa。刚度矩阵 Q 与柔度矩阵 S 互为逆矩阵，即 Q = S⁻¹，这是校验计算是否正确的重要依据。
三、主轴方向与偏轴方向的区分
上述 S、Q 矩阵只在材料主轴（1-2）下成立。当单层板铺放角 θ ≠ 0° 时，应力应变关系需通过坐标变换得到偏轴柔度 S̄ 与偏轴刚度 Q̄，此时会出现 Q̄₁₆、Q̄₂₆ 等耦合项，表示正应力与剪应变之间产生耦合（拉剪耦合）。这正是层合板设计中通过改变铺层角来“调刚度”的物理根源。本节只讨论主轴下的关系，偏轴变换将在后续章节展开。
判断单层板主轴本构是否用对，可做三项检查：一是 S 与 Q 是否互为逆矩阵；二是 ν₁₂/E₁ 是否等于 ν₂₁/E₂；三是给定应变后算出的应力是否满足 σ₁ 远大于 σ₂（当 ε₁ 与 ε₂ 同量级时）。若三项都通过，说明矩阵构建正确。
四、从矩阵到工程判断
掌握 S 与 Q 之后，就能快速估算单层板在给定载荷下的变形。例如对上述单层板施加 σ₁ = 100 MPa、σ₂ = 0、τ₁₂ = 0，则 ε₁ = S₁₁σ₁ = 7.41×10⁻³ × 0.1 = 7.41×10⁻⁴，ε₂ = S₁₂σ₁ = −2.22×10⁻³ × 0.1 = −2.22×10⁻⁴，横向收缩约为纵向伸长的 30%，与 ν₁₂ 一致。若改用刚度形式，先由 ε 求 σ，则 σ₁ = Q₁₁ε₁ + Q₁₂ε₂，代入同样数值可还原 100 MPa，验证两套矩阵自洽。
在层合板分析中，单层板的 Q 矩阵是后续计算层合板面内刚度矩阵 A、耦合矩阵 B、弯曲刚度矩阵 D 的基本输入。每一层的 Q 先经偏轴变换得到 Q̄，再按铺层角与厚度叠加。因此本节的主轴本构关系是整门课程的计算起点，务必做到能由 E₁、E₂、ν₁₂、G₁₂ 手算 S 与 Q，并能用 Q = S⁻¹ 自查。
小结：单层板在材料主轴下由 4 个独立工程常数完全描述；柔度矩阵 S 描述“应力→应变”，刚度矩阵 Q 描述“应变→应力”，两者互逆；ν₁₂/E₁ = ν₂₁/E₂ 是必须满足的对称条件。掌握这三条，就为后续偏轴变换与层合板理论打好了基础。