3、复合材料力学基础：各向异性、正交各向异性、单层板与层合板的基本力学概念
复合材料结构设计与金属结构设计的根本差异，在于材料属性的“可设计性”与“方向依赖性”。金属材料通常可视为各向同性体，其弹性模量、泊松比等力学参数在不同方向上保持一致；而复合材料（尤其是纤维增强树脂基复合材料）在纤维方向、基体方向以及厚度方向上的力学性能截然不同。这种各向异性特征决定了我们在进行结构设计时，不能简单套用传统材料力学公式，而必须建立以“单层板—层合板”为主线的分析体系。本章将系统梳理各向异性、正交各向异性、单层板与层合板的核心力学概念，为后续的强度校核、铺层设计与失效分析奠定理论基础。
3.1 核心原理：从各向异性到正交各向异性的简化逻辑
连续介质力学中，各向异性体的广义胡克定律可表示为 σij = Cijkl·εkl（i, j, k, l = 1, 2, 3）。由于应力张量和应变张量均为对称张量，弹性刚度张量 Cijkl 的独立分量从理论上可缩减为21个。这21个独立弹性常数是描述任意各向异性材料弹性行为的最一般形式，但工程上极少有材料需要如此复杂的描述。当材料内部存在三个相互正交的弹性对称面时，即构成正交各向异性体，其独立弹性常数从21个锐减至9个。对于纤维增强复合材料，我们通常假设纤维在单层板内沿某一方向（1方向）均匀排列，且纤维与基体界面结合良好，此时单层板在面内呈现正交各向异性，在厚度方向（3方向）上则表现出横观各向同性特征。
进一步简化发生在“单层板”层面。由于复合材料层合板的面内尺寸远大于厚度尺寸，且层合板在面内承受载荷时，厚度方向的正应力 σ3 通常可忽略不计，因此引入平面应力假设。在该假设下，单层板的本构关系简化为：
σ1 = Q11·ε1 + Q12·ε2
σ2 = Q12·ε1 + Q22·ε2
τ12 = Q66·γ12其中 Qij 为缩减刚度系数，可由工程弹性常数 E1（纤维方向弹性模量）、E2（横向弹性模量）、G12（面内剪切模量）、ν12（主泊松比）计算得到。值得注意的是，ν21 并非独立常数，而是由 ν21 = ν12·E2/E1 确定，这体现了正交各向异性材料的能量守恒约束。
3.2 判断方法：单层板刚度矩阵的构建与坐标变换
在实际结构中，单层板的纤维方向往往与层合板的主坐标轴（如结构的长边方向）存在一个铺层角 θ。因此，我们需要将材料主方向（1-2坐标系）下的刚度矩阵 [Q] 变换到结构坐标系（x-y坐标系）下，得到变换刚度矩阵 [Q̄]。这一变换过程是层合板理论的核心操作步骤。
- 确定材料主方向刚度矩阵 [Q]：根据单层板的工程常数 E1、E2、G12、ν12，计算 Q11 = E1/(1-ν12·ν21)，Q22 = E2/(1-ν12·ν21)，Q12 = ν12·E2/(1-ν12·ν21)，Q66 = G12。
- 计算坐标变换矩阵 [T]：对于铺层角 θ（纤维方向与x轴的夹角），变换矩阵 [T] 是 θ 的正弦、余弦函数组合。
- 求解变换刚度矩阵 [Q̄]：通过 [Q̄] = [T]⁻¹[Q][T] 计算得到。展开后，[Q̄] 的各个分量是 θ 的周期函数，其中 Q̄11、Q̄22、Q̄12、Q̄66 随 θ 变化，且会出现耦合项 Q̄16、Q̄26（拉剪耦合项）。
- 判断耦合效应：当 θ = 0° 或 90° 时，Q̄16 = Q̄26 = 0，此时单层板在结构坐标系下不产生拉剪耦合；当 θ 为其他角度时，拉剪耦合不可忽略。
| 铺层角 θ | Q̄16 / Q̄26 | 工程应用场景 | 典型特性 |
|---|---|---|---|
| 0° | 0 | 承受轴向拉伸/压缩 | 纤维方向刚度最大 |
| 90° | 0 | 承受横向载荷 | 横向刚度最大，纵向刚度低 |
| ±45° | 非零（值较大） | 承受剪切载荷 | 剪切模量贡献显著，拉剪耦合明显 |
| 任意角 θ | 非零 | 复杂应力状态 | 需通过变换矩阵精确计算 |
3.3 案例拆解：碳纤维/环氧树脂单层板的刚度计算与验证
以典型的T300碳纤维/环氧树脂单层板为例，其工程常数为：E1 = 135 GPa，E2 = 9.5 GPa，G12 = 5.2 GPa，ν12 = 0.32。我们计算其在0°、45°、90°铺层角下的变换刚度矩阵，并分析其力学响应差异。
背景与目标：某无人机机翼蒙皮采用该材料体系，设计人员需要评估不同铺层角对蒙皮面内刚度的影响，以确定最优铺层方案。
计算过程：首先计算 ν21 = ν12·E2/E1 = 0.32 × 9.5 / 135 ≈ 0.0225。然后计算 [Q] 矩阵：Q11 = 135 / (1 - 0.32×0.0225) ≈ 136.0 GPa，Q22 = 9.5 / (1 - 0.32×0.0225) ≈ 9.57 GPa，Q12 = 0.32 × 9.5 / (1 - 0.32×0.0225) ≈ 3.06 GPa，Q66 = 5.2 GPa。接着对 θ = 45° 进行坐标变换，得到 Q̄11 = 43.2 GPa，Q̄22 = 43.2 GPa，Q̄12 = 28.6 GPa，Q̄66 = 36.5 GPa，Q̄16 = 30.1 GPa，Q̄26 = 30.1 GPa。
结果对比：0°铺层时，Q̄11 = 136 GPa，刚度极高；90°铺层时，Q̄11 = 9.57 GPa，刚度很低；45°铺层时，Q̄11 介于两者之间，但出现了显著的拉剪耦合项 Q̄16 = 30.1 GPa。这意味着在45°铺层下，施加纯拉伸载荷会引起剪切变形，反之亦然。
3.4 层合板理论：从单层到整体的刚度合成
层合板由多个单层板按特定顺序和角度叠合而成。经典层合板理论（CLT）基于以下假设：层合板各层之间完美粘接、无相对滑移；层合板厚度方向的应变 ε3 忽略不计；变形前垂直于中面的直线段在变形后仍保持直线且垂直于中面（直法线假设）。在这些假设下，层合板中任意一点的应变可表示为中面应变 ε⁰ 与曲率 κ 的线性组合：ε = ε⁰ + z·κ。
通过对每一单层板在厚度方向积分，可以得到层合板的内力（合力 N 和合力矩 M）与中面应变和曲率之间的关系，即著名的 ABD 矩阵方程：
[ N ] [ A B ] [ ε⁰ ]
[ M ] = [ B D ] [ κ ]其中，A 矩阵（面内刚度矩阵）表示层合板抵抗面内变形的能力，D 矩阵（弯曲刚度矩阵）表示抵抗弯曲变形的能力，B 矩阵（耦合刚度矩阵）则表征了面内变形与弯曲变形之间的耦合关系。当铺层顺序对称于中面时（如 [0/90]s），B 矩阵为零矩阵，层合板不会产生拉弯耦合；反之，非对称铺层会产生显著的 B 矩阵，导致结构在承受拉伸时发生弯曲变形，这在工程中通常需要避免。
| 矩阵 | 物理含义 | 计算公式 | 设计关注点 |
|---|---|---|---|
| A | 面内刚度 | Aij = Σ(Q̄ij)k·(zk - zk-1) | 决定层合板拉伸/压缩刚度 |
| B | 拉弯耦合 | Bij = ½Σ(Q̄ij)k·(zk² - zk-1²) | 非对称铺层时不为零，引起翘曲 |
| D | 弯曲刚度 | Dij = ⅓Σ(Q̄ij)k·(zk³ - zk-1³) | 决定层合板抗弯能力 |
3.5 常见误区与总结复盘
误区一：将单层板的工程常数直接用于层合板整体分析。层合板的刚度不是各单层板刚度的简单平均，必须通过 ABD 矩阵进行合成。例如，[0/90]s 层合板的面内刚度 A11 并非 (Q11_0° + Q11_90°)/2，而是需要按各层厚度加权求和。
误区二：忽视 B 矩阵的影响。许多初学者在设计铺层时只关注 A 矩阵的刚度值，而忽略了非对称铺层带来的 B 矩阵。一旦 B 矩阵非零，层合板在固化冷却过程中就会因热残余应力而产生翘曲变形，导致制件尺寸超差。
误区三：混淆材料主方向与结构坐标方向。在有限元分析中，必须为每个单元指定材料方向。若材料方向定义错误，计算结果将完全失真。建议在建模时使用局部坐标系，并仔细核对铺层角度。
总结复盘：本章从各向异性的一般理论出发，逐步简化至正交各向异性单层板，再通过坐标变换和层合理论合成层合板的 ABD 刚度矩阵。这一逻辑链条是复合材料结构设计的力学基石。在实际工程中，我们应始终牢记：单层板是“材料”，层合板是“结构”。材料性能通过试验测定，结构性能则通过 CLT 计算或有限元仿真获得。下一阶段的学习将围绕强度理论（如 Tsai-Wu 准则）和失效分析展开，届时本章的刚度矩阵知识将作为核心输入。