第二章 各向异性弹性力学基础

各位工程师朋友,大家好。欢迎来到复合材料结构强度计算的第一道坎——各向异性弹性力学。说实话,我当年刚接触这个领域时,也被一堆张量公式搞得头晕。但后来我发现,只要抓住几个核心概念,这东西其实没那么可怕。

2.1 应力与应变张量

先说说应力。我们熟悉的金属材料,通常用三个正应力和三个剪应力就能描述一个点的受力状态。但复合材料不一样,它内部结构复杂,应力分布也更微妙。

应力的数学描述

在三维空间中,一点的应力状态需要用二阶张量来表示。我习惯用矩阵形式写出来:

σ = [σ11  σ12  σ13]
    [σ21  σ22  σ23]
    [σ31  σ32  σ33]

这里σ11、σ22、σ33是正应力,其余是剪应力。注意,根据剪应力互等定理,σ12 = σ21,σ13 = σ31,σ23 = σ32。所以实际上只有6个独立分量。

我的小经验: 做有限元分析时,我习惯把应力分量按 [σ11, σ22, σ33, σ23, σ13, σ12] 的顺序排列。这个顺序在ABAQUS和ANSYS里都是默认的,能减少很多麻烦。

应变张量

应变张量的形式和应力类似:

ε = [ε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]
注意: 这个6×6的刚度矩阵是对称的,即Cij = Cji。我曾经见过有人把非对称矩阵输入到有限元软件里,结果算出来的变形方向都是错的。

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)
关键点: 工程常数必须满足热力学稳定性条件。比如E1、E2、E3必须为正,泊松比必须满足一定的约束。我在做材料参数拟合时,经常遇到试验数据不满足这些条件的情况,这时候就需要重新审视试验方法了。

知识体系总览

下面这张图是我自己总结的,把本章的核心逻辑串起来了:

各向异性弹性力学知识体系 应力与应变张量 广义胡克定律 正交各向异性 应力张量 σij = σji(对称) 6个独立分量 应变张量 εij 刚度矩阵 Cij σ = C · ε 21个独立分量(一般) 9个独立分量(正交) 三个对称面 1-2-3 材料主方向 刚度矩阵简化 无拉剪耦合 工程常数:E1, E2, E3, νij, Gij 刚度矩阵 ↔ 工程常数 相互转换 核心逻辑:从张量描述 → 本构关系 → 材料对称性 → 工程应用参数

这张图把本章的知识脉络理清了。从应力应变张量出发,通过广义胡克定律建立关系,再引入正交各向异性的简化,最后落到工程常数上。每一步都有明确的物理意义和工程背景。

避坑指南: 我曾经在做一个风电叶片项目时,直接用各向同性的公式去算复合材料层合板的刚度,结果偏差了30%以上。后来才意识到,正交各向异性的耦合效应不能忽略。所以,拿到材料参数后,第一件事就是确认材料的主方向。

好了,这一章的内容就到这里。工程常数与刚度矩阵的转换公式,我建议你亲手推导一遍,这样印象更深。下一章我们会讨论单层板的强度理论,到时候这些基础知识都会用上。


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