第二章:弹性力学基础——应力张量、应变张量、广义胡克定律
各位同学,欢迎来到弹性力学基础这一讲。说实话,很多做仿真的工程师,软件操作溜得很,但一遇到非线性收敛问题、或者结果明显不对时,就抓瞎了。为什么?因为底层的力学概念没吃透。我个人习惯是,不管用什么软件,先把弹性力学这三板斧——应力张量、应变张量、广义胡克定律——搞明白。今天咱们就掰开揉碎了聊一聊。
2.1 应力张量:不只是“力除以面积”那么简单
你想想看,中学物理里学的应力,就是力除以面积,一个数值搞定。但在三维世界里,事情没那么简单。一个微元体上,每个面都有三个方向的力分量。九个分量凑在一起,就构成了应力张量。
我习惯用矩阵形式来表示它:
σ = [ σxx σxy σxz ]
[ σyx σyy σyz ]
[ σzx σzy σzz ]
这里有个关键点:应力张量是对称的。也就是说,σxy = σyx,σxz = σzx,σyz = σzy。为什么?因为微元体不能自己转起来,力矩必须平衡。我在项目中遇到过一位同事,手动计算时把剪应力分量搞反了,结果整个结构的变形方向都错了。嗯,这个坑希望大家别踩。
核心概念:应力张量描述的是物体内部一点处的应力状态。它不依赖于你选的坐标系,但它的分量会随坐标系旋转而改变。这就是“张量”的本质——客观存在,但表达方式可以不同。
2.2 应变张量:变形到底有多“细腻”
有了应力,自然要谈应变。应变张量描述的是物体变形后的几何变化。说白了,就是看一个微元体被拉长了多少、扭曲了多少。
应变张量的定义式是这样的:
εij = (1/2) * ( ∂ui/∂xj + ∂uj/∂xi )
其中u是位移场,x是坐标。这个公式看着有点吓人,其实逻辑很简单:对角线上的εxx、εyy、εzz是正应变,表示拉伸或压缩;非对角线上的εxy、εyz、εzx是剪应变,表示角度的变化。
我记得刚做仿真那会儿,总搞不清工程剪应变γ和张量剪应变ε的区别。后来被一个老工程师点醒了:工程剪应变γ = 2ε。如果你在软件里输入材料参数时搞混了这两个,结果会差一倍。我曾经因为这个错误,把一个悬臂梁的挠度算小了一半,查了三天才找到原因。
实用技巧:在ABAQUS或ANSYS中,输出的应变分量默认是张量应变。如果你需要工程应变,记得手动乘以2。尤其是做复合材料层合板分析时,这个细节特别重要。
2.3 广义胡克定律:连接应力与应变的桥梁
胡克定律大家都不陌生,F = kx嘛。但在三维弹性力学里,它变成了一个张量方程:
σij = Cijkl * εkl
这里Cijkl是弹性刚度张量,有81个分量。别慌,对于最常见的各向同性材料,它只取决于两个独立常数:杨氏模量E和泊松比ν。
写成矩阵形式,各向同性材料的广义胡克定律长这样:
[ σxx ] [ λ+2μ λ λ 0 0 0 ] [ εxx ]
[ σyy ] [ λ λ+2μ λ 0 0 0 ] [ εyy ]
[ σzz ] = [ λ λ λ+2μ 0 0 0 ] [ εzz ]
[ σxy ] [ 0 0 0 μ 0 0 ] [ γxy ]
[ σyz ] [ 0 0 0 0 μ 0 ] [ γyz ]
[ σzx ] [ 0 0 0 0 0 μ ] [ γzx ]
其中λ和μ是拉梅常数,μ也叫剪切模量G。它们与E、ν的关系是:
| 参数 | 表达式 |
|---|---|
| λ(拉梅第一常数) | λ = Eν / [(1+ν)(1-2ν)] |
| μ = G(剪切模量) | μ = E / [2(1+ν)] |
⚠️ 特别注意:当ν接近0.5时(比如橡胶材料),λ会趋于无穷大,这就是“不可压缩”现象。在仿真中,如果你用普通单元处理ν=0.5的材料,会出现体积锁死。我建议使用杂交单元或减缩积分单元来避免这个问题。
2.4 知识体系总览
说了这么多,咱们用一张图把这三者的关系串起来。下面是我手绘的SVG结构图,展示了弹性力学基础的核心逻辑:
这张图把咱们今天讲的内容串起来了。你看,应力张量和应变张量通过广义胡克定律建立联系,再结合平衡方程和几何方程,就构成了弹性力学问题的完整求解框架。说白了,你只要搞懂这三者,大部分线弹性仿真问题都能拿捏住。
2.5 实际应用中的几点提醒
- 坐标系一致性:定义材料方向时,一定要和你的全局坐标系或局部坐标系对齐。我见过有人把正交各向异性材料的1方向和2方向搞反了,结果算出来的变形完全不对。
- 单位制:应力张量的单位是Pa,应变张量无量纲。在输入材料参数时,E和ν的单位要统一。我习惯用MPa和mm制,这样应力数值比较直观。
- 后处理解读:看应力结果时,别只看Mises应力。主应力、最大剪应力有时更能揭示失效机理。比如脆性材料,我更关注最大主应力;韧性材料,Mises应力更有参考价值。
我的一个小习惯:每次建完模型,我都会手动算一个简单工况(比如单轴拉伸)来验证材料参数和边界条件是否正确。花10分钟做这个验证,能省下后面几天的调试时间。
好了,弹性力学基础就聊到这儿。记住,应力张量、应变张量、广义胡克定律这三者,是仿真分析的基石。你越理解它们,软件在你手里就越听话。下次遇到收敛问题或者结果异常,不妨回头看看是不是这些基础概念上出了偏差。