📚 课程目录第3章复合材料层合板理论回顾

3. 复合材料层合板理论回顾

在 OptiStruct 中做复合材料优化，无论后续是自由尺寸优化、铺层顺序优化还是尺寸优化，求解器最终都要回到同一套层合板力学模型上：用经典层合板理论（CLT）把单层材料性能、铺层角度和铺层厚度，折算成层合板的等效刚度，再据此计算每层的应力和应变，最后用失效准则判断该层是否安全。理解这条链路，才能看懂 OptiStruct 输出的 PCOMP/PCOMPG 卡片、MAT8 材料卡以及优化响应（如失效指标、强度比）之间的对应关系。

一、经典层合板理论的三个基本假设

CLT 把层合板简化为一块薄板，并作如下假设：

层合板各层之间完全粘接，无相对滑移，变形连续；
板厚度方向的应变可忽略，即 ε_z ≈ 0，各层处于平面应力状态；
直法线假设：变形前垂直于中面的直线，变形后仍为直线且保持垂直，因此横向剪应变 γ_xz、γ_yz 为零。

由此，板内任意一点的位移可由中面位移 u₀、v₀、w₀ 表示，面内应变沿厚度方向呈线性分布：

εx = εx0 + z·κx
εy = εy0 + z·κy
γxy = γxy0 + z·κxy

其中 ε⁰ 为中面应变，κ 为中面曲率，z 为该点距中面的坐标。这一线性关系是后面 ABD 矩阵能够成立的根本原因。

二、单层材料在材料主轴下的应力应变关系

单层板（lamina）在材料主方向 1-2 下是正交各向异性材料，平面应力状态下有：

σ1   = Q11·ε1 + Q12·ε2
σ2   = Q12·ε1 + Q22·ε2
τ12 = Q66·γ12

其中折算刚度 Q 由工程常数给出：Q11 = E1/(1−ν12ν21)，Q22 = E2/(1−ν12ν21)，Q12 = ν12E2/(1−ν12ν21)，Q66 = G12。注意 ν21 = ν12·E2/E1，这是 OptiStruct 中 MAT8 卡片必须同时给出 E1、E2、ν12、G12 的原因。

当铺层角度为 θ 时，需要把材料主方向的 Q 转换到层合板整体坐标系（x-y）下，得到 Q̄。转换涉及 m = cosθ、n = sinθ 的四次项，是铺层角度影响刚度的直接来源。θ = 0° 时 Q̄ = Q；θ = 90° 时 1、2 方向互换；θ = ±45° 时会出现 Q̄16、Q̄26 耦合项，这正是 ±45° 铺层能提供剪切刚度的原因。

三、ABD 矩阵：层合板的等效刚度

把各层的 Q̄ 沿厚度积分，就得到层合板的内力—应变关系：

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

其中 A 为拉伸刚度矩阵，B 为拉伸—弯曲耦合矩阵，D 为弯曲刚度矩阵。对第 k 层，其贡献为：

Aij = Σ Q̄ij(k)·(zk − zk−1)
Bij = ½·Σ Q̄ij(k)·(zk² − zk−1²)
Dij = ⅓·Σ Q̄ij(k)·(zk³ − zk−1³)

这里 z_k 是第 k 层上表面到中面的距离。可以看出：A 只与各层厚度有关，D 与厚度三次方有关，B 则取决于铺层关于中面的对称性。对称铺层时 B = 0，层合板不产生拉弯耦合，这也是工程上优先采用对称铺层的原因。OptiStruct 的 PCOMP 卡片中，通过 LAYER 行的 T（厚度）和 THETA（角度）以及 SOUT 选项，实际上就是在定义这些 z_k 和 Q̄_ij。

_allow / σ_applied R > 1 → 安全 R < 1 → 失效层内应力应变关系（逐层计算） ε⁰, κ → ε(z) = ε⁰ + z·κ → σ(z) = Q̄·ε(z) → 转换到材料主方向 σ₁, σ₂, τ₁₂ OptiStruct 输出：每层应力、应变、失效指标、强度比优化响应可选用 FAILURE 或 STRESS 作为约束

四、逐层应力计算与失效准则

得到 ABD 后，对给定载荷 {N, M} 求解线性方程组，得到中面应变 ε⁰ 和曲率 κ。然后对每一层，在其上下表面 z 处计算：

εx(z) = εx0 + z·κx
εy(z) = εy0 + z·κy
γxy(z) = γxy0 + z·κxy

再通过坐标转换得到材料主方向应变 ε₁、ε₂、γ₁₂，乘以 Q 得到 σ₁、σ₂、τ₁₂。由于每层角度不同，同一块板内各层的应力状态差异很大，这正是复合材料分析必须逐层进行的原因。

常用失效准则包括：

最大应力/最大应变准则：各分量独立与许用值比较，简单但忽略耦合；
Tsai-Hill 准则：二次交互式，适用于拉压强度相近的材料；
Tsai-Wu 准则：引入一次项，能区分拉压强度不同，是工程上最常用的准则之一；
Hashin 准则：区分纤维失效与基体失效，适合需要判断失效模式的场景。

以 Tsai-Wu 为例，其失效判据为：

F1·σ1 + F2·σ2 + F11·σ1² + F22·σ2²
+ 2·F12·σ1·σ2 + F66·τ12² = 1

其中 F1 = 1/Xt − 1/Xc，F11 = 1/(Xt·Xc)，F2、F22 同理，F66 = 1/S²。当左边小于 1 时该层安全，等于 1 时达到失效边界。

五、强度比 R 的计算与工程判读

强度比（Strength Ratio）R 定义为：把当前应力状态按比例放大，直到刚好满足失效判据时的放大倍数。对 Tsai-Wu 这类二次判据，令 σ → R·σ 代入，得到关于 R 的二次方程：

a·R² + b·R − 1 = 0
R = [−b + √(b² + 4a)] / (2a)

其中 a 为二次项系数，b 为一次项系数。R > 1 表示该层在当前载荷下安全，R 越大裕度越高；R < 1 表示已经失效；R = 1 表示刚好达到临界。OptiStruct 在复合材料分析中会输出每层的强度比或失效指标，优化时通常把最小强度比作为约束，要求其不小于 1（或按设计规范取 1.5 等安全系数）。

在 OptiStruct 中，MAT8 卡片定义单层工程常数和强度值（XT、XC、YT、YC、S），PCOMP/PCOMPG 定义铺层角度、厚度和失效准则（通过 FT 字段选择）。优化响应中可用 FAILURE 提取失效指标，用 STRESS 提取层应力。若使用 Tsai-Wu，需注意 F12 的取值，OptiStruct 默认采用 F12 = −0.5√(F11·F22)，也可通过 MAT8 的 F12 字段显式指定。

强度比 R 与安全系数不是同一个概念。R 是基于当前应力状态和所选失效准则计算出的裕度指标，换用不同准则（如从 Tsai-Hill 改为 Tsai-Wu）会得到不同的 R 值。因此在优化设置中必须固定失效准则，并在结果对比时保持同一准则，否则约束值不可比。

六、一个可手算的示例

取单层材料：E₁ = 135 GPa，E₂ = 9 GPa，G₁₂ = 5 GPa，ν₁₂ = 0.3，Xt = 1500 MPa，Xc = 1200 MPa，Yt = 50 MPa，Yc = 200 MPa，S = 70 MPa。假设某 0° 层计算得到 σ₁ = 600 MPa，σ₂ = 20 MPa，τ₁₂ = 30 MPa。用 Tsai-Wu 计算：

F1 = 1/1500 − 1/1200 = −1.667e-4
F11 = 1/(1500×1200) = 5.556e-7
F2 = 1/50 − 1/200 = 0.015
F22 = 1/(50×200) = 1.0e-4
F66 = 1/70² = 2.041e-4
F12 = −0.5×√(5.556e-7 × 1.0e-4) = −3.727e-6

a = F11σ1² + F22σ2² + 2F12σ1σ2 + F66τ12²
  = 5.556e-7×360000 + 1.0e-4×400
    + 2×(−3.727e-6)×12000 + 2.041e-4×900
  = 0.200 + 0.040 − 0.0894 + 0.1837 = 0.3343
b = F1σ1 + F2σ2 = −1.667e-4×600 + 0.015×20
  = −0.100 + 0.300 = 0.200

R = [−0.200 + √(0.04 + 4×0.3343)] / (2×0.3343)
  = [−0.200 + √1.3772] / 0.6686
  = (−0.200 + 1.1735) / 0.6686 ≈ 1.456

R ≈ 1.46 > 1，该层安全。若把 σ₁ 提高到 1000 MPa，重新计算可得 R 明显下降，接近 1 时即进入临界状态。这个手算过程与 OptiStruct 内部逐层计算逻辑一致，只是软件会对所有层、所有单元批量执行。

掌握这条从铺层定义 → Q̄ 转换 → ABD 集成 → 逐层应力 → 失效判据 → 强度比的完整链路，后续在 OptiStruct 中设置复合材料优化时，就能清楚每个参数和响应分别作用在链路的哪一环，也能判断优化结果是否合理。