第3章复合材料力学基础

3. 复合材料力学基础：各向异性弹性力学、经典层合板理论、失效准则

各位工程师朋友，咱们今天聊点硬核的。做防弹复合材料，你绕不开力学基础。说白了，你得知道材料怎么变形、怎么受力、什么时候会坏。我刚开始接触这行的时候，也啃过不少理论书，但真正让我理解的，还是项目里一次次试错和复盘。今天我就把压箱底的经验掏出来，跟你们捋一捋。

金属材料，你拉它一下，它各个方向反应差不多。这叫各向同性。但复合材料不一样。你想想，碳纤维沿着纤维方向硬得像钢筋，垂直方向却软得像塑料。这就是各向异性。

描述这种“方向性”的力学行为，我们需要广义胡克定律。公式长这样：

σ_ij = C_ijkl * ε_kl

这里 C_ijkl 是刚度矩阵，有 81 个分量。别怕，实际用起来，因为对称性，会简化很多。对于正交各向异性材料（比如单向带），独立常数只有 9 个：E1、E2、E3、ν12、ν13、ν23、G12、G13、G23。

关键点： 在防弹分析中，我们最关心面内性能。通常假设材料处于平面应力状态（σ3=0）。这时候，刚度矩阵简化为 3x3。

我遇到过不少新手，拿着金属那套 E 和 ν 就往复合材料上套，结果算出来的变形量差了好几倍。嗯，这里要注意：复合材料的面内剪切模量 G12 和泊松比 ν12，跟金属完全不是一个概念。你想想看，纤维和基体之间的剪切耦合，复杂得很。

防弹板很少是单层，都是好几层不同角度叠起来的。经典层合板理论，就是用来分析这种“叠层子”的力学行为的。我个人习惯，把它拆成三步走：

核心公式是 ABD 矩阵：

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

其中：

避坑指南： 我曾经设计过一块 [0/90]s 的层合板，以为对称铺层就没有耦合了。结果一算，B 矩阵确实为零，但层间剪切应力在边缘处大得惊人。所以，CLT 算的是“板内”响应，边缘效应得另想办法。

为什么 CLT 在防弹分析里这么重要？因为弹击是局部冲击，会产生剧烈的弯曲变形。这时候，D 矩阵决定了板子抵抗变形的能力。你想想看，如果 D 矩阵太小，弹头一撞，板子背面就鼓个大包，背面层先裂了。

材料再好，也有扛不住的时候。失效准则就是判断“什么时候坏”的数学判据。防弹领域，我主要用两个：Tsai-Wu 和 Hashin。

Tsai-Wu 准则是一个二次张量多项式。它用一个光滑的曲面来包络所有应力状态。公式如下：

F1*σ1 + F2*σ2 + F11*σ1^2 + F22*σ2^2 + F66*τ12^2 + 2*F12*σ1*σ2 = 1

其中 F1、F2 等系数，由材料的强度参数（Xt、Xc、Yt、Yc、S）决定。

这个准则的好处是：简单、统一。不管你是拉、压、剪，还是组合应力，一个公式搞定。但缺点也很明显：它不区分失效模式。你只知道它坏了，但不知道是纤维断了还是基体裂了。

注意： Tsai-Wu 准则中的 F12 项，需要双轴试验数据来确定。如果没有，通常取 F12 = -0.5 * sqrt(F11*F22)。我建议你保守一点，取 F12 = 0，这样预测的失效包络线更安全。

Hashin 准则就聪明多了。它把失效分成四种模式：

每种模式有各自的判据。比如纤维拉伸失效：

(σ1 / Xt)^2 + (τ12 / S)^2 = 1

基体拉伸失效：

(σ2 / Yt)^2 + (τ12 / S)^2 = 1

你看，Hashin 准则能告诉你：哦，是基体先裂了，还是纤维先断了。这对防弹设计太重要了。因为纤维断裂意味着板子彻底失去承载能力，而基体开裂可能还能扛几发。

我做过一个对比：用 Tsai-Wu 算出来失效载荷是 1000N，用 Hashin 算出来基体失效在 800N，纤维失效在 1200N。实际打靶测试，板子在 850N 左右出现背面裂纹，1100N 时穿透。你看，Hashin 更贴近物理过程。

下面这张图，是我自己总结的本章知识逻辑。你一看就明白：

理论归理论，上了项目你就知道，真实情况比书本复杂得多。比如：

所以，我的习惯是：先用 Tsai-Wu 做初步筛选，再用 Hashin 做详细分析，最后一定要用试验验证。你想想看，理论算得再漂亮，打一枪就露馅了，那多尴尬。

一个小技巧： 做层合板铺层设计时，尽量让 [B] 矩阵为零（对称铺层），这样可以避免拉弯耦合导致的翘曲变形。我在做防弹头盔面层设计时，就吃过这个亏——不对称铺层导致成型后头盔变形，贴合度很差。

好了，这一章的内容就到这里。力学基础是后面所有分析的地基，你花时间把它啃透了，后面做冲击仿真、优化设计，就会顺手很多。