2. 力学基础(一):各向异性弹性力学基础、正交各向异性材料本构方程、单层板的宏观力学分析

各位同行,咱们今天聊点硬核的。复合材料的结构分析,说白了就是跟「方向」打交道。你想想看,一块金属,你从哪个方向拉它,它反应都差不多。但复合材料不一样——纤维方向和非纤维方向,那简直是两个世界。

我个人习惯,在讲任何复合材料结构设计之前,先把力学基础夯扎实。因为后面所有的铺层设计、强度校核、失效分析,追根溯源都回到这里。今天这一章,我们就从最底层的各向异性弹性力学开始,一步步走到单层板的宏观力学分析。

2.1 各向异性弹性力学基础

先问一个问题:为什么复合材料不能用传统的各向同性力学?

答案很简单——因为它的性能随方向变化。在经典弹性力学里,我们描述一个点的应力状态需要6个应力分量,应变状态也需要6个应变分量。它们之间通过广义胡克定律联系起来:

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

这里的Cᵢⱼ就是刚度矩阵,一共36个分量。听起来很吓人对吧?

但别慌。根据应变能密度函数的对称性,Cᵢⱼ = Cⱼᵢ,所以独立分量从36个降到了21个。这就是各向异性材料的一般情况——21个独立弹性常数。

关键概念: 各向异性材料是指材料在所有方向上力学性能都不同。对于复合材料来说,这通常是因为纤维的定向排列造成的。

我在项目中遇到过一位刚入行的同事,他拿着各向同性的公式去算复合材料板的弯曲变形,结果差了30%以上。嗯,从那以后他再也不敢乱套公式了。

2.2 正交各向异性材料本构方程

好在,我们常用的复合材料(比如碳纤维/环氧树脂)并不是完全的各向异性。它们有一个重要的对称性——正交各向异性

什么叫正交各向异性?说白了,就是材料有三个互相垂直的对称面。对于单向带(UD层),我们通常定义:

  • 1方向:纤维方向(纵向)
  • 2方向:垂直于纤维方向(横向)
  • 3方向:厚度方向

在这种坐标系下,刚度矩阵大大简化。21个独立分量降到了9个:

| σ₁₁ |   | C₁₁  C₁₂  C₁₃  0    0    0   | | ε₁₁ |
| σ₂₂ |   | C₁₂  C₂₂  C₂₃  0    0    0   | | ε₂₂ |
| σ₃₃ | = | C₁₃  C₂₃  C₃₃  0    0    0   | | ε₃₃ |
| τ₂₃ |   | 0    0    0    C₄₄  0    0   | | γ₂₃ |
| τ₁₃ |   | 0    0    0    0    C₅₅  0   | | γ₁₃ |
| τ₁₂ |   | 0    0    0    0    0    C₆₆ | | γ₁₂ |

你看,很多位置都是0。这意味着正应力与剪应变之间没有耦合,剪应力与正应变之间也没有耦合。这就是正交各向异性带来的简化。

实用技巧: 在实际工程中,我们更常用工程常数(E₁, E₂, ν₁₂, G₁₂等)而不是刚度矩阵Cᵢⱼ。因为工程常数更直观,测试方法也更成熟。两者可以通过公式互相转换,我建议你手边常备一份转换公式表。

为什么会这样?因为纤维增强复合材料在制造过程中,纤维是平行排列的,这就天然形成了三个对称面。你想想看,如果纤维乱成一团,那性能就不可控了。

2.3 单层板的宏观力学分析

好了,现在我们有了正交各向异性的本构方程。但实际结构分析中,我们很少直接处理三维问题。对于薄壁结构(比如机翼蒙皮、机身壁板),我们通常做平面应力假设

所谓平面应力,就是假设σ₃₃ = τ₂₃ = τ₁₃ = 0。这个假设在薄板分析中非常合理,因为厚度方向的应力确实很小。

于是,三维的本构方程简化为二维形式:

| σ₁₁ |   | Q₁₁  Q₁₂  0   | | ε₁₁ |
| σ₂₂ | = | Q₁₂  Q₂₂  0   | | ε₂₂ |
| τ₁₂ |   | 0    0    Q₆₆ | | γ₁₂ |

这里的Qᵢⱼ称为缩减刚度矩阵,它是由三维刚度矩阵Cᵢⱼ经过平面应力条件推导出来的。具体公式我就不列了,但你要记住:Q₁₁ = E₁/(1-ν₁₂ν₂₁),Q₂₂ = E₂/(1-ν₁₂ν₂₁),Q₁₂ = ν₁₂E₂/(1-ν₁₂ν₂₁),Q₆₆ = G₁₂。

注意: 这里的ν₂₁不是独立的,它满足ν₂₁/E₂ = ν₁₂/E₁。这个关系叫互等定理,千万别搞反了。我曾经见过有人把ν₁₂和ν₂₁当成两个独立参数输入,结果算出来的刚度矩阵不对称——这在物理上是不可能的。

单层板的宏观力学分析,核心就是建立这个二维本构关系。有了它,我们就可以计算单层板在任意载荷下的应力应变响应。

举个例子:假设一个碳纤维/环氧单向带,E₁ = 140 GPa,E₂ = 10 GPa,ν₁₂ = 0.3,G₁₂ = 5 GPa。如果沿纤维方向施加100 MPa的应力,你会得到什么结果?

嗯,算一下:ε₁₁ = σ₁₁/E₁ = 100/140000 ≈ 0.000714,也就是0.0714%的应变。而横向应变ε₂₂ = -ν₁₂ε₁₁ = -0.3 × 0.000714 ≈ -0.000214。你看,纵向刚度远大于横向刚度,这就是纤维增强的效果。

为了让你更直观地理解本章的知识体系,我画了一张图:

第二章 知识体系:各向异性弹性力学 → 正交各向异性 → 单层板分析 各向异性弹性力学基础 广义胡克定律:σᵢ = Cᵢⱼ εⱼ(i, j = 1,...,6) 21个独立弹性常数 → 对称性简化 正交各向异性材料本构方程 三个互相垂直的对称面 → 9个独立弹性常数 刚度矩阵Cᵢⱼ → 工程常数(E₁, E₂, ν₁₂, G₁₂) 正应力与剪应变解耦,剪应力与正应变解耦 单层板的宏观力学分析 平面应力假设:σ₃₃ = τ₂₃ = τ₁₃ = 0 缩减刚度矩阵Qᵢⱼ → 二维本构关系 应用:计算单层板在任意载荷下的应力应变响应 核心:从三维到二维,从复杂到实用

这张图清晰地展示了本章的逻辑链条:从最一般的各向异性弹性力学出发,利用正交各向异性的对称性简化本构方程,最后通过平面应力假设得到工程实用的单层板二维本构关系。每一步都是必要的简化,每一步都有明确的物理背景。

好了,这一章的内容就到这里。记住,力学基础是复合材料结构设计的根基。你花时间把这里搞懂了,后面铺层设计、强度分析、甚至损伤容限分析都会轻松很多。

本章要点回顾:

  • 各向异性材料有21个独立弹性常数
  • 正交各向异性材料简化为9个独立常数
  • 平面应力假设下,单层板用缩减刚度矩阵Qᵢⱼ描述
  • 工程常数(E₁, E₂, ν₁₂, G₁₂)是工程应用中最常用的参数

专注资料整理