3. 各向异性弹性与断裂力学基础

复合材料层合板的断裂行为之所以比金属复杂，根源在于它的弹性本构是各向异性的。金属可以近似用两个常数（弹性模量 E 和泊松比 ν）描述，而单向纤维复合材料在纤维方向、垂直纤维方向和厚度方向上的刚度差异可达一个数量级。理解这种差异如何进入断裂力学的表达，是后续分析分层、纤维断裂和基体开裂的前提。

一、正交各向异性本构：从广义胡克定律到工程常数

单向复合材料在宏观上可视为正交各向异性体，具有三个相互垂直的弹性对称面。取材料主轴坐标系 1-2-3，其中 1 为纤维方向，2 为面内垂直纤维方向，3 为厚度方向。广义胡克定律写作 σ = Cε，刚度矩阵 C 为 9×9，但正交各向异性使其独立常数减少到 9 个。工程上更常用柔度形式 ε = Sσ，面内应力状态下（σ₃ = τ₂₃ = τ₁₃ = 0）简化为：

ε₁ = σ₁/E₁ − ν₁₂σ₂/E₁
ε₂ = −ν₁₂σ₁/E₁ + σ₂/E₂
γ₁₂ = τ₁₂/G₁₂

这里出现了四个面内工程常数：纵向模量 E₁、横向模量 E₂、主泊松比 ν₁₂ 和面内剪切模量 G₁₂。它们并非独立，柔度矩阵对称性要求 ν₁₂/E₁ = ν₂₁/E₂，因此只需三个独立常数即可确定面内弹性行为。以典型 T300/环氧单向板为例，E₁ ≈ 135 GPa，E₂ ≈ 9 GPa，G₁₂ ≈ 5 GPa，ν₁₂ ≈ 0.30。可以看到 E₁/E₂ ≈ 15，这种强烈的方向性直接导致裂纹尖端应力场不再具有金属那样的自相似性。

当层合板中存在离轴铺层时，需要将主轴刚度通过转轴公式变换到层合板坐标系。变换后会出现耦合项 Q₁₆、Q₂₆，它们把正应力与剪应变耦合起来。这正是 ±45° 铺层在拉伸时产生剪切变形、在剪切时产生正应力的原因，也是复合材料断裂分析中必须考虑铺层顺序的力学基础。

二、各向异性断裂力学的基本量

在线弹性断裂力学框架下，裂纹尖端附近的应力场由应力强度因子 K 控制。对正交各向异性材料，I 型裂纹尖端应力场可写为：

σ_ij = K_I / √(2πr) · f_ij(θ, 弹性常数)

与各向同性材料不同，角分布函数 f_ij 不再只依赖 θ，还依赖材料的弹性常数。这意味着同样的 K_I，在 E₁/E₂ 不同的材料中，裂纹尖端的应力分布形状不同。对于复合材料，K 的定义仍然成立，但必须明确裂纹面相对于材料主轴的取向——沿纤维裂纹、垂直纤维裂纹和纤维/基体界面裂纹，其 K 的标定和临界值都不同。

能量释放率 G 与应力强度因子的关系也随各向异性改变。对各向同性材料，平面应变下 G = K²(1−ν²)/E。对正交各向异性材料，I 型裂纹在自相似扩展条件下有：

G_I = K_I² · √( (a₁₁a₂₂/2) · [√(a₂₂/a₁₁) + (2a₁₂+a₆₆)/(2a₁₁)] )

其中 a_ij 为柔度矩阵面内分量。这个表达式说明：即使 K_I 相同，不同材料因柔度组合不同，G_I 也不同。工程中更常用的是直接测量 G 而非 K，因为层合板的 G 可以通过双悬臂梁（DCB）等试样直接获得，且与铺层和界面状态直接相关。

三、层间断裂与混合模式

复合材料最典型的断裂形式是分层，它发生在层间，属于 I/II 混合模式。I 型分量由层间正应力驱动，II 型分量由层间剪应力驱动。混合模式下的能量释放率通常用幂律准则或 B-K 准则描述：

G_I/G_Ic + G_II/G_IIc = 1        （线性）
(G_I/G_Ic)^α + (G_II/G_IIc)^β = 1 （幂律）
G_c = G_Ic + (G_IIc − G_Ic)·(G_II/(G_I+G_II))^η  （B-K）

以碳纤维/环氧层合板为例，典型值 G_Ic ≈ 0.2～0.3 kJ/m²，G_IIc ≈ 0.6～1.0 kJ/m²，η 约 2～3。这些数值说明：层间抵抗剪切分层的能力明显高于抵抗张开分层的能力，因此纯 I 型分层最危险。工程判断中，若计算得到某界面 G_I/G_Ic + G_II/G_IIc ≥ 1，即认为该界面可能起裂。

四、从本构到断裂判据的完整链条

把上述内容串起来，复合材料断裂分析的基本链条是：铺层角度与材料主轴 → 转轴得到层合板刚度 → 计算层间应力或裂纹尖端 K → 换算为 G_I、G_II → 与 G_Ic、G_IIc 比较 → 判断是否起裂及扩展方向。这条链条中任何一步的坐标系或符号约定出错，都会导致判据失效。因此在实际分析中，必须始终明确：1 方向是纤维方向，裂纹面法向与铺层界面的关系，以及 K 和 G 所对应的模式。

这张图把本节的逻辑关系集中呈现出来：左侧是材料参数，右侧是从主轴到判据的传递链。需要特别注意的是，图中 K 到 G 的换算不是简单的系数关系，而是包含柔度组合的平方根项。当铺层角度改变时，a₁₁、a₂₂、a₁₂、a₆₆ 都会变化，因此同一裂纹在不同铺层中的 G 值可能相差数倍。这也是为什么复合材料分层分析必须结合具体铺层顺序，而不能套用金属的 K-G 换算公式。

掌握以上内容后，面对一个具体的层合板分层问题，就可以按“确定铺层与主轴 → 查取或实测工程常数 → 计算层间应力或 K → 换算 G_I、G_II → 代入混合模式判据”的顺序完成定量判断。下一阶段将在此基础上引入界面断裂韧度的测试方法和数值实现细节。