第二章 设计基础与力学原理

各位工程师朋友,欢迎来到复合材料结构设计的核心地带。说实话,这一章的内容是整门课程的基石。你想想看,连材料怎么受力、怎么变形都没搞明白,后面谈铺层设计、谈强度校核,那都是空中楼阁。

我刚开始接触复合材料时,也犯过这样的错误——拿着各向同性材料的思维去套。结果呢?算出来的变形方向完全不对,被老工程师狠狠批了一顿。从那以后,我老老实实把各向异性弹性力学和层合板理论啃了下来。今天,我就把这些核心内容掰开了讲给你听。

2.1 各向异性弹性力学基础

复合材料跟金属最大的区别在哪?说白了,就是它的性能跟方向有关。金属你往哪个方向拉,弹性模量基本都一样。但复合材料不一样——沿着纤维方向硬得像钢,垂直纤维方向软得像塑料。

为什么会这样?因为复合材料是由纤维和基体两种材料组成的。纤维承担主要载荷,基体负责传递和固定。这个结构决定了它的力学行为天生就是各向异性的。

描述这种各向异性行为,我们需要用到广义胡克定律:

σᵢ = Cᵢⱼ · εⱼ   (i, j = 1, 2, ..., 6)

这里的Cᵢⱼ是刚度矩阵,一共有36个分量。不过别怕,实际工程中没这么复杂。对于正交各向异性材料(大多数复合材料层合板都属于这一类),独立的分量只有9个。

我给大家整理了一个常用材料的刚度矩阵形式:

材料类型 独立弹性常数个数 典型常数
各向同性(金属) 2 E, ν
横观各向同性 5 E₁, E₂, ν₁₂, G₁₂, ν₂₃
正交各向异性 9 E₁, E₂, E₃, ν₁₂, ν₁₃, ν₂₃, G₁₂, G₁₃, G₂₃
个人经验: 我在做风电叶片设计时,经常用横观各向同性来近似单向带。因为厚度方向的性能跟面内横向差不多,少测几个参数能省不少试验费。但如果你做的是高精度分析,还是老老实实用正交各向异性吧。

2.2 经典层合板理论(CLT)

CLT,全称Classical Laminate Theory。这个名字听起来挺唬人,其实核心思想就一句话:把多层板当成一个整体来算

它的基本假设有三个:

  • 层间粘结完美,没有滑移
  • 变形前垂直于中面的直线,变形后仍然是直线且垂直于中面
  • 厚度方向的应变忽略不计

嗯,这里要注意。第三个假设在薄板分析时没问题,但如果你做的是厚板或者局部冲击分析,这个假设就不成立了。我曾经在一个船用夹芯板项目里吃过这个亏——用CLT算出来的变形偏小,后来换成一阶剪切变形理论才把问题搞定。

CLT的核心公式是:

⎧ N ⎫   ⎡ A  B ⎤ ⎧ ε⁰ ⎫
⎨   ⎬ = ⎢      ⎥ ⎨     ⎬
⎩ M ⎭   ⎣ B  D ⎦ ⎩ κ  ⎭

其中:

  • A矩阵:面内刚度,控制拉伸和剪切
  • D矩阵:弯曲刚度,控制弯曲和扭转
  • B矩阵:耦合刚度,拉伸和弯曲会互相影响
避坑指南: 我曾经遇到一个新手,他设计的对称铺层B矩阵算出来是零,但实际加工时因为铺层角度偏差,B矩阵非零了。结果一加载,板子自己翘起来了。所以我的建议是:即使理论上是零,也要留一点余量。

2.3 应力-应变关系

搞清楚了层合板的整体行为,我们还得回到单层板层面。毕竟,破坏往往是从某一层开始的。

单层板的应力-应变关系,在材料主方向坐标系下是这样的:

⎧ σ₁ ⎫   ⎡ Q₁₁  Q₁₂   0  ⎤ ⎧ ε₁ ⎫
⎪    ⎪   ⎢               ⎥ ⎪    ⎪
⎨ σ₂ ⎬ = ⎢ Q₁₂  Q₂₂   0  ⎥ ⎨ ε₂ ⎬
⎪    ⎪   ⎢               ⎥ ⎪    ⎪
⎩ τ₁₂⎭   ⎣  0    0   Q₆₆ ⎦ ⎩ γ₁₂⎭

这里的Q矩阵是缩减刚度矩阵,跟工程常数之间的关系是:

  • Q₁₁ = E₁ / (1 - ν₁₂ν₂₁)
  • Q₂₂ = E₂ / (1 - ν₁₂ν₂₁)
  • Q₁₂ = ν₁₂E₂ / (1 - ν₁₂ν₂₁)
  • Q₆₆ = G₁₂

你可能会问:那如果载荷方向跟材料主方向不一致怎么办?好问题!这时候就需要做坐标变换。说白了,就是把应力应变从材料坐标系转到全局坐标系。这个变换矩阵我就不列了,你查一下教材都有。我个人的习惯是写一个小脚本自动算,省得每次手算出错。

2.4 强度理论

强度理论,说白了就是判断材料什么时候坏。复合材料跟金属不一样,它的破坏模式多得很——纤维断裂、基体开裂、分层、纤维-基体脱粘……每种模式对应不同的判据。

2.4.1 Tsai-Wu 强度理论

Tsai-Wu是一个张量多项式形式的强度理论。它的好处是统一了各种破坏模式,用一个公式搞定:

F₁σ₁ + F₂σ₂ + F₁₁σ₁² + F₂₂σ₂² + F₆₆τ₁₂² + 2F₁₂σ₁σ₂ = 1

这个公式里的F系数,都是通过试验测出来的。我个人的经验是:Tsai-Wu在预测拉伸主导的破坏时比较准,但在压缩工况下容易偏保守。所以如果你做的是压缩设计,建议再搭配一个Hashin准则做交叉验证。

注意: Tsai-Wu中的F₁₂项很难测准。我见过不少设计报告,直接取F₁₂ = 0。这样做虽然省事,但会引入误差。如果你做的是高可靠性设计(比如航空结构),建议通过双轴试验来标定。

2.4.2 Hashin 强度理论

Hashin准则跟Tsai-Wu不一样,它把破坏模式分开了:

  • 纤维拉伸破坏:σ₁ > 0时,用(σ₁/Xᵀ)² + (τ₁₂/S)² = 1
  • 纤维压缩破坏:σ₁ < 0时,用(σ₁/Xᶜ)² = 1
  • 基体拉伸破坏:σ₂ > 0时,用(σ₂/Yᵀ)² + (τ₁₂/S)² = 1
  • 基体压缩破坏:σ₂ < 0时,用(σ₂/(2Sᵀ))² + [(Yᶜ/(2Sᵀ))² - 1]σ₂/Yᶜ + (τ₁₂/Sᵀ)² = 1

我个人更偏爱Hashin准则,原因很简单——它能告诉你哪一层、哪种模式先破坏。这对后续的铺层优化太重要了。比如你发现基体先坏了,那就调整铺层角度,让载荷更多地由纤维承担。

我记得有一次做无人机机翼设计,用Tsai-Wu算出来安全裕度0.85,感觉还行。但用Hashin一查,发现有一层的基体已经进入塑性了。后来调整了铺层顺序,把0°层从中间移到表面,问题就解决了。你看,这就是区分破坏模式的价值。

2.5 本章知识体系

为了帮你理清思路,我画了一张图,把这一章的核心逻辑串起来:

复合材料结构设计基础 · 知识体系 设计基础与力学原理 各向异性弹性力学 经典层合板理论 强度理论 广义胡克定律 刚度矩阵 · 工程常数 坐标变换 ABD矩阵 面内刚度 · 弯曲刚度 耦合刚度 Tsai-Wu · Hashin 破坏模式区分 安全裕度计算 工程应用:铺层设计 · 强度校核 · 失效分析 从材料属性 → 层合板响应 → 强度判断 → 设计迭代

这张图把本章的逻辑关系讲得很清楚。从左到右,从材料属性到层合板响应,再到强度判断,最后落到工程应用。你每次做设计时,都可以沿着这个路径走一遍,保证不会漏掉关键环节。


好了,这一章的内容就到这里。记住,理论是工具,不是目的。我见过太多人把CLT背得滚瓜烂熟,但一到实际项目就不知道怎么用。我的建议是:多动手算几个例子,把ABD矩阵自己推导一遍,把Tsai-Wu和Hashin的判据写个小程序跑一跑。只有亲手做过,这些公式才能真正变成你自己的东西。

公众号:蓝海资料掘金营,微信deep3321