2. 弹性力学基础：应力与应变概念、胡克定律、平衡方程与几何方程、物理方程

各位同学，咱们今天聊聊弹性力学基础。说实话，这章是整个有限元分析的根基。我当年刚接触有限元时，总觉得这些概念太抽象，后来做项目摔过几次跟头，才明白——嗯，基础不牢，地动山摇。

2.1 应力与应变：从直觉到数学

先说说应力。你想想看，一根杆子拉长了，内部到底发生了什么？

应力，说白了就是单位面积上的内力。我习惯把它想象成「材料内部的抵抗感」。你拉它，它不想被拉长，这种「不情愿」就是应力。

但实际工程中，应力没那么简单。一个微元体上有九个应力分量：

• 正应力：σ_xx、σ_yy、σ_zz——垂直于截面的应力
• 剪应力：τ_xy、τ_yz、τ_zx——平行于截面的应力

我在项目中遇到过一件事：一个同事算出来的应力值明明没超限，结果构件还是断了。后来一查，是剪应力被忽略了。所以啊，千万别只看正应力。

应变就好理解了——变形量与原始尺寸的比值。同样有正应变和剪应变：

• 正应变 ε = ΔL / L
• 剪应变 γ = 角度变化量（弧度）

2.2 胡克定律：材料最朴素的脾气

胡克定律，说白了就是「你拉我多少，我反抗多少」。线性关系：

σ = E · ε

这个 E 就是弹性模量，也叫杨氏模量。它反映了材料抵抗变形的能力。钢的 E 大约是 210 GPa，铝是 70 GPa——所以钢比铝硬三倍。

但注意了，胡克定律只适用于线弹性范围。我曾经见过一个新手，把塑性区的数据也用胡克定律算，结果偏差了30%以上。

对于三维问题，胡克定律要写成矩阵形式：

σ = D · ε

其中 D 是弹性矩阵，对于各向同性材料：
D = E/(1+ν)(1-2ν) × 
[1-ν  ν   ν   0   0   0
 ν  1-ν  ν   0   0   0
 ν   ν  1-ν  0   0   0
 0   0   0   (1-2ν)/2  0   0
 0   0   0   0   (1-2ν)/2  0
 0   0   0   0   0   (1-2ν)/2]

看着复杂，其实核心就两个参数：E 和 ν（泊松比）。

2.3 平衡方程：力不能白来

平衡方程讲的是——物体内部任意一点，受力必须平衡。说白了就是：

力的总和 = 0，力矩的总和 = 0

在三维情况下，平衡方程是三个偏微分方程：

∂σ_xx/∂x + ∂τ_xy/∂y + ∂τ_xz/∂z + f_x = 0
∂τ_yx/∂x + ∂σ_yy/∂y + ∂τ_yz/∂z + f_y = 0
∂τ_zx/∂x + ∂τ_zy/∂y + ∂σ_zz/∂z + f_z = 0

其中 f_x、f_y、f_z 是体力分量（比如重力）。

我刚开始学的时候，觉得这方程太完美了——现实中哪有这么理想？后来做有限元才明白，平衡方程是「强形式」，有限元用的是「弱形式」（虚功原理）。但强形式是基础，不理解它，你就不懂有限元在近似什么。

2.4 几何方程：变形与位移的关系

几何方程，说白了就是「位移怎么变成应变」。假设物体内一点位移为 u、v、w，那么：

ε_xx = ∂u/∂x
ε_yy = ∂v/∂y
ε_zz = ∂w/∂z
γ_xy = ∂u/∂y + ∂v/∂x
γ_yz = ∂v/∂z + ∂w/∂y
γ_zx = ∂w/∂x + ∂u/∂z

这里有个关键点：这些公式假设变形很小。如果变形很大，就要用格林应变张量了。我在做橡胶密封件分析时就吃过这个亏——小变形假设下算出来密封性很好，实际一测，漏了。

2.5 物理方程：把前面串起来

物理方程就是胡克定律的广义形式。它把应力和应变联系起来：

σ = D · ε

到这里，弹性力学的三大方程就齐了：

方程类型	作用	未知量
平衡方程	力的平衡	应力 σ
几何方程	位移→应变	应变 ε、位移 u
物理方程	应变→应力	应力 σ、应变 ε

你看，三个方程，六个未知量（三个应力、三个应变、三个位移），刚好能解。这就是弹性力学的基本框架。

2.6 知识体系总览

下面这张图，是我自己总结的弹性力学核心逻辑。每次做有限元分析前，我都会在脑子里过一遍：

这张图我建议你存下来。每次做有限元分析前，对着它捋一遍思路，能避免很多低级错误。

好了，弹性力学基础就讲到这里。记住：平衡方程管力，几何方程管变形，物理方程管材料。三者缺一不可。下一章咱们聊聊有限元法的基本思想——怎么把这些连续方程变成计算机能算的离散形式。

---

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