2. 复合材料力学基础：各向异性材料本构关系、经典层合板理论、强度准则简介
各位同事，大家好。今天我们来聊聊复合材料力学的几个核心基础。说实话，搞叶片制造应力场分析，如果不懂材料怎么“受力变形”，后面全是瞎忙活。我个人习惯，每次拿到一个新铺层设计，第一件事不是急着算，而是先把材料本构和层合板理论在脑子里过一遍。
这一节，我打算讲三个东西：各向异性本构关系、经典层合板理论，还有强度准则。嗯，听起来有点枯燥，但我会尽量结合项目里的坑来讲。
2.1 各向异性材料本构关系
先问大家一个问题：为什么复合材料不能像金属那样用一个弹性模量E和泊松比ν就搞定？
答案很简单——因为它是各向异性的。说白了，沿着纤维方向和垂直于纤维方向，刚度差了一个数量级。你想想看，碳纤维沿轴向的模量能到230GPa，横向可能只有10GPa左右。这要是还用各向同性公式去算，结果肯定离谱。
对于正交各向异性材料，应力-应变关系在材料主方向（1-纤维方向，2-面内垂直方向，3-厚度方向）上可以写成：
| ε₁ | [ S₁₁ S₁₂ S₁₃ 0 0 0 ] | σ₁ |
| ε₂ | | S₁₂ S₂₂ S₂₃ 0 0 0 | | σ₂ |
| ε₃ | = | S₁₃ S₂₃ S₃₃ 0 0 0 | | σ₃ |
| γ₂₃| | 0 0 0 S₄₄ 0 0 | | τ₂₃|
| γ₁₃| | 0 0 0 0 S₅₅ 0 | | τ₁₃|
| γ₁₂| | 0 0 0 0 0 S₆₆ | | τ₁₂|这里的S矩阵就是柔度矩阵。工程上我们更常用刚度矩阵Q，其实就是S的逆矩阵。我记得刚入行那会儿，总搞不清S和Q的转换关系，后来发现只要记住“柔度是应变对应力的响应，刚度是应力对应变的抵抗”，就顺了。
关键点：对于平面应力问题（σ₃=0, τ₂₃=0, τ₁₃=0），本构关系可以简化为：
| σ₁ | [ Q₁₁ Q₁₂ 0 ] | ε₁ |
| σ₂ | = | Q₁₂ Q₂₂ 0 | | ε₂ |
| τ₁₂| [ 0 0 Q₆₆] | γ₁₂|其中Q₁₁ = E₁/(1-ν₁₂ν₂₁)，Q₂₂ = E₂/(1-ν₁₂ν₂₁)，Q₁₂ = ν₁₂E₂/(1-ν₁₂ν₂₁)，Q₆₆ = G₁₂。
我的小技巧：实际计算时，我习惯先把单层板的工程常数（E₁, E₂, G₁₂, ν₁₂）列个表，再代入公式算Q矩阵。这样不容易出错，也方便后续做参数敏感性分析。
2.2 经典层合板理论
单层板搞清楚了，接下来就是怎么把多层板叠在一起。经典层合板理论（CLT）是分析层合板应力应变的基础工具。它的核心假设是：变形前垂直于中面的直线，变形后仍然保持直线且垂直于中面。嗯，其实就是忽略横向剪切变形。
基于这个假设，层合板内任意一点的应变可以写成：
εₓ = εₓ⁰ + z·κₓ
εᵧ = εᵧ⁰ + z·κᵧ
γₓᵧ = γₓᵧ⁰ + z·κₓᵧ其中ε⁰是中面应变，κ是中面曲率，z是距离中面的坐标。
然后，通过积分各单层板的应力，可以得到层合板的合力与合力矩：
| N | [ A B ] | ε⁰ |
| M | = [ B D ] | κ |这里的A矩阵是拉伸刚度，D矩阵是弯曲刚度，B矩阵是耦合刚度。B矩阵不为零时，拉伸和弯曲会耦合在一起——你拉它一下，它不光伸长，还会弯。我在做叶片根部铺层设计时就遇到过这个问题，一开始没注意B矩阵，结果算出来的变形跟实测差了30%。后来把铺层顺序调成对称的，B矩阵归零，问题就解决了。
避坑指南：我曾经在分析一个非对称铺层的叶片时，直接用CLT算应力，结果在厚度方向上的应力分布完全不对。后来发现，对于厚叶片或者局部曲率大的区域，横向剪切变形不能忽略，得用一阶剪切变形理论（FSDT）才行。CLT只适用于薄板，这个边界条件一定要记牢。
2.3 强度准则简介
有了应力场，接下来就是判断材料会不会失效。复合材料的强度准则比金属复杂得多，因为不同方向的强度不一样，而且拉压强度也不同。
常用的准则有这几个：
- 最大应力准则：各方向的应力分别不超过对应的强度值。简单粗暴，但忽略了应力之间的相互作用。
- 最大应变准则：类似，但用应变来判断。我个人觉得这个比最大应力准则更准一些，尤其是对于纤维控制的情况。
- Tsai-Hill准则：考虑了应力之间的耦合，形式是：
(σ₁/X)² - (σ₁σ₂/X²) + (σ₂/Y)² + (τ₁₂/S)² = 1其中X、Y、S分别是纵向、横向和剪切强度。这个准则在工程上用得很多，我一般用它做初步筛选。
- Tsai-Wu准则：更通用，能处理拉压强度不同的情况。形式是：
F₁σ₁ + F₂σ₂ + F₁₁σ₁² + F₂₂σ₂² + F₆₆τ₁₂² + 2F₁₂σ₁σ₂ = 1这个准则的参数标定比较麻烦，尤其是F₁₂项。我记得有一次做叶片疲劳分析，Tsai-Wu准则算出来的失效指数总是偏大，后来发现是F₁₂取错了。建议各位在做实验标定时，多测几个应力组合，别偷懒。
我的建议：实际项目中，我通常先用最大应力准则做快速判断，再用Tsai-Hill或Tsai-Wu做详细校核。如果时间紧，Tsai-Hill就够了；如果涉及拉压不对称或复杂应力状态，还是得上Tsai-Wu。
知识体系总览
下面这张图是我自己整理的，把这一节的核心逻辑串起来了。你从单层板的本构出发，经过坐标变换和层合板理论，最后落到强度校核上。每一步都有坑，每一步也都有技巧。
好了，这一节的内容就这些。记住，理论是基础，但真正用的时候一定要结合工艺和实验数据。下一节我们会讲制造过程中的热-力耦合问题，到时候再聊。
公众号:蓝海资料掘金营，微信deep3321