📚 课程目录第2章复合材料层合板基础

2. 复合材料层合板基础

分层是层合板最典型的失效形式之一，而分层的驱动力几乎都来自层间应力。要理解层间应力从何而来，必须先把四件事串起来：铺层坐标系如何定义、单层刚度如何表达、层合板如何把各层“组装”成整体、以及层间应力在哪个环节被计算出来。本节按这条逻辑链展开，为后续分层分析建立统一的坐标与力学基础。

一、铺层坐标系：先约定，再计算

层合板分析涉及三套坐标系，必须严格区分并保持一致：

材料主坐标系 1-2-3：1 为纤维方向，2 为面内垂直纤维方向，3 为单层厚度方向。单层刚度在 1-2-3 下定义。
层合板坐标系 x-y-z：x、y 为层合板面内参考方向，z 为厚度方向，通常取中面为 z=0，向上为正。
铺层角 θ：从 x 轴逆时针转到 1 轴为正，即 θ 是材料主方向相对层合板参考方向的角度。

例如 [0/90/±45]s 中，0° 层纤维沿 x，90° 层纤维沿 y，+45° 层从 x 逆时针转 45°，−45° 层顺时针转 45°，s 表示对称铺层。铺层顺序从底面到顶面书写，z 坐标随之确定。所有正负号、角度、厚度方向约定一旦确定，全文与图示必须一致，否则层间剪应力方向会整体反号。

二、单层刚度：从工程常数到 Q 矩阵

单层在材料主坐标系下是正交各向异性材料，平面应力状态下由四个独立工程常数描述：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̄。转换的核心是方向余弦 m=cosθ、n=sinθ，Q̄₁₁、Q̄₂₂、Q̄₁₂、Q̄₁₆、Q̄₂₆、Q̄₆₆ 中会出现 m²n² 与 m³n、mn³ 项。正是这些交叉项，使 ±45° 层在面内剪切与正应变之间产生耦合，也是层间剪应力出现的重要来源。

教学示例（合理假设，非实测）：取典型碳纤维/环氧单层 E₁=135 GPa、E₂=9 GPa、ν₁₂=0.30、G₁₂=5 GPa、单层厚度 0.125 mm。可算得 Q₁₁≈136.6 GPa、Q₂₂≈9.11 GPa、Q₁₂≈2.73 GPa、Q₆₆=5 GPa。这些数值用于后续示例计算，不代表任何具体产品实测值。

三、层合板经典理论：ABD 矩阵把各层组装起来

经典层合板理论（CLT）建立在直法线假设与平面应力假设之上，核心是把每层的刚度沿厚度积分，得到面内力 N、面内力矩 M 与中面应变 ε⁰、曲率 κ 的关系：

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

A_ij = Σ Q̄_ij^(k) · (z_k - z_{k-1})
B_ij = ½ Σ Q̄_ij^(k) · (z_k² - z_{k-1}²)
D_ij = ⅓ Σ Q̄_ij^(k) · (z_k³ - z_{k-1}³)

A 是拉伸刚度，D 是弯曲刚度，B 是拉弯耦合刚度。对称铺层时 B=0，层合板不产生拉弯耦合；非对称铺层 B≠0，受拉即产生弯曲，层间应力显著增大。求解时先由外载得到 N、M，再解 6×6 方程组得到 ε⁰ 与 κ，最后由 ε(z)=ε⁰+zκ 得到各层应变，乘 Q̄ 得到各层面内应力 σ_x、σ_y、τ_xy。

四、层间应力计算入口：从面内应力到三维应力

CLT 只给出面内应力，层间应力 σ_z、τ_xz、τ_yz 并不在 CLT 的直接输出中。层间应力的计算入口通常有三条：

层间应力恢复：利用三维平衡方程 ∂σ_x/∂x+∂τ_xy/∂y+∂τ_xz/∂z=0 等，对已知面内应力沿厚度积分，逐层恢复 τ_xz、τ_yz，再由 ∂τ_xz/∂x+∂τ_yz/∂y+∂σ_z/∂z=0 恢复 σ_z。
自由边效应：在自由边附近，面内应力沿厚度突变，层间应力在边界层内急剧升高，是分层起裂的高发位置。
分层驱动力：层间正应力 σ_z 与层间剪应力 τ_xz、τ_yz 共同构成分层驱动力，配合界面强度或断裂韧性判据判断是否起裂与扩展。

下面用一张信号链式示意图，把“铺层角 → 单层刚度 → ABD 组装 → 面内应力 → 层间应力恢复 → 分层驱动力”这条主线串起来，便于快速定位每一步的输入与输出。

五、示例：对称层合板的层间剪应力估算

以 [0/90]s 层合板受单轴拉伸 N_x 为例，说明层间应力恢复的基本做法。取单层厚度 t=0.125 mm，总厚 0.5 mm，中面为 z=0，各层 z 坐标为：−0.25～−0.125（0°）、−0.125～0（90°）、0～0.125（90°）、0.125～0.25（0°）。

由 CLT 解出中面应变后，各层面内 σ_x 在 0° 层与 90° 层中数值不同，沿厚度呈阶梯分布。对 τ_xz 用平衡方程 ∂σ_x/∂x+∂τ_xz/∂z=0 沿厚度积分：从自由表面 σ_z=0、τ_xz=0 出发，逐层累加，得到 τ_xz 沿厚度的分布。由于 0° 与 90° 层 σ_x 差异大，界面处 τ_xz 出现峰值，这正是分层容易在 0/90 界面萌生的力学原因。

注意：CLT 假设直法线，无法直接给出层间应力；层间应力恢复是后处理步骤，其精度依赖面内应力梯度 ∂σ_x/∂x 的准确程度。自由边附近梯度极大，恢复结果对网格与边界条件敏感，工程上常配合三维有限元或分层单元校核。

六、本节要点回顾

铺层角 θ 从 x 轴逆时针到 1 轴为正，铺层顺序从底面到顶面书写，z 以中面为零、向上为正。
单层刚度由 E₁、E₂、ν₁₂、G₁₂ 确定 Q，再经 θ 转换得 Q̄，交叉项带来拉剪耦合。
ABD 矩阵把各层沿厚度积分组装，对称铺层 B=0，非对称铺层 B≠0。
CLT 给出面内应力；层间应力 σ_z、τ_xz、τ_yz 需经三维平衡方程恢复，自由边是分层起裂高发区。
层间正应力与剪应力共同构成分层驱动力，是后续分层判据与扩展分析的输入。