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。你看,纵向刚度远大于横向刚度,这就是纤维增强的效果。
为了让你更直观地理解本章的知识体系,我画了一张图:
这张图清晰地展示了本章的逻辑链条:从最一般的各向异性弹性力学出发,利用正交各向异性的对称性简化本构方程,最后通过平面应力假设得到工程实用的单层板二维本构关系。每一步都是必要的简化,每一步都有明确的物理背景。
好了,这一章的内容就到这里。记住,力学基础是复合材料结构设计的根基。你花时间把这里搞懂了,后面铺层设计、强度分析、甚至损伤容限分析都会轻松很多。
本章要点回顾:
- 各向异性材料有21个独立弹性常数
- 正交各向异性材料简化为9个独立常数
- 平面应力假设下,单层板用缩减刚度矩阵Qᵢⱼ描述
- 工程常数(E₁, E₂, ν₁₂, G₁₂)是工程应用中最常用的参数