第二章 各向异性弹性力学基础
各位工程师朋友,大家好。欢迎来到复合材料结构强度计算的第一道坎——各向异性弹性力学。说实话,我当年刚接触这个领域时,也被一堆张量公式搞得头晕。但后来我发现,只要抓住几个核心概念,这东西其实没那么可怕。
2.1 应力与应变张量
先说说应力。我们熟悉的金属材料,通常用三个正应力和三个剪应力就能描述一个点的受力状态。但复合材料不一样,它内部结构复杂,应力分布也更微妙。
应力的数学描述
在三维空间中,一点的应力状态需要用二阶张量来表示。我习惯用矩阵形式写出来:
σ = [σ11 σ12 σ13]
[σ21 σ22 σ23]
[σ31 σ32 σ33]
这里σ11、σ22、σ33是正应力,其余是剪应力。注意,根据剪应力互等定理,σ12 = σ21,σ13 = σ31,σ23 = σ32。所以实际上只有6个独立分量。
应变张量
应变张量的形式和应力类似:
ε = [ε11 ε12 ε13]
[ε21 ε22 ε23]
[ε31 ε32 ε33]
工程上常用工程剪应变γ,它与张量剪应变的关系是:γ12 = 2ε12。这个转换我刚开始也搞混过,后来在编写材料子程序时吃了亏才记住。
2.2 广义胡克定律
对于一般各向异性材料,应力与应变的关系可以写成:
σij = Cijkl · εkl
其中Cijkl是刚度张量,有81个分量。但别被这个数字吓到,实际上由于对称性,独立分量只有21个。
写成矩阵形式更直观:
[σ11] [C11 C12 C13 C14 C15 C16] [ε11]
[σ22] [C21 C22 C23 C24 C25 C26] [ε22]
[σ33] = [C31 C32 C33 C34 C35 C36] [ε33]
[σ23] [C41 C42 C43 C44 C45 C46] [ε23]
[σ13] [C51 C52 C53 C54 C55 C56] [ε13]
[σ12] [C61 C62 C63 C64 C65 C66] [ε12]
2.3 正交各向异性材料
复合材料中最常见的就是正交各向异性材料。说白了,就是材料有三个互相垂直的对称面。比如单向纤维增强复合材料,沿着纤维方向、垂直于纤维方向、厚度方向,性能都不一样。
对于正交各向异性材料,刚度矩阵简化为:
[C11 C12 C13 0 0 0 ]
[C12 C22 C23 0 0 0 ]
[C13 C23 C33 0 0 0 ]
[ 0 0 0 C44 0 0 ]
[ 0 0 0 0 C55 0 ]
[ 0 0 0 0 0 C66]
独立分量从21个减少到9个。嗯,这里要注意,虽然矩阵简化了,但每个分量背后都有明确的物理意义。
2.4 单层板的工程常数
在实际工程中,我们很少直接用刚度矩阵Cij。我更习惯用工程常数——弹性模量、泊松比、剪切模量。这些参数更直观,也更容易从试验中获得。
正交各向异性单层板的工程常数
| 符号 | 含义 | 典型值(碳纤维/环氧) |
|---|---|---|
| E1 | 纵向弹性模量 | 140 GPa |
| E2 | 横向弹性模量 | 10 GPa |
| E3 | 厚度方向弹性模量 | 10 GPa |
| ν12 | 纵向-横向泊松比 | 0.3 |
| ν13 | 纵向-厚度方向泊松比 | 0.3 |
| ν23 | 横向-厚度方向泊松比 | 0.4 |
| G12 | 纵向-横向剪切模量 | 5 GPa |
| G13 | 纵向-厚度方向剪切模量 | 5 GPa |
| G23 | 横向-厚度方向剪切模量 | 3.5 GPa |
这些工程常数与刚度矩阵的转换关系是:
C11 = (1 - ν23·ν32) / (E2·E3·Δ)
C12 = (ν21 + ν31·ν23) / (E1·E3·Δ)
...
其中 Δ = (1 - ν12·ν21 - ν23·ν32 - ν31·ν13 - 2·ν21·ν32·ν13) / (E1·E2·E3)
知识体系总览
下面这张图是我自己总结的,把本章的核心逻辑串起来了:
这张图把本章的知识脉络理清了。从应力应变张量出发,通过广义胡克定律建立关系,再引入正交各向异性的简化,最后落到工程常数上。每一步都有明确的物理意义和工程背景。
好了,这一章的内容就到这里。工程常数与刚度矩阵的转换公式,我建议你亲手推导一遍,这样印象更深。下一章我们会讨论单层板的强度理论,到时候这些基础知识都会用上。
公众号:蓝海资料掘金营,微信deep3321