2. 弹性力学基础：应力与应变概念、胡克定律、弹性矩阵、平面应力与平面应变问题

各位同学，欢迎来到第二章。这一章我们聊聊弹性力学基础。说实话，这部分内容看起来全是公式和概念，但它是整个有限元分析的根基。我当年刚接触有限元时，就是没把应力应变的关系吃透，结果算出来的结果跟实验对不上，折腾了好几天才发现是弹性矩阵搞错了。嗯，咱们今天就把这些基础打牢。

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

先说说应力。你想想看，一根杆子被拉长，内部到底发生了什么？说白了，应力就是单位面积上的内力。我习惯把它想象成「材料内部的抵抗感」。

应力的数学表达是：

σ = F / A

其中 F 是力，A 是截面积。但实际工程中，应力可不是这么简单的一个数。它是一个张量，有六个独立分量：三个正应力（σx, σy, σz）和三个剪应力（τxy, τyz, τzx）。

应变呢？它是描述材料变形的量。我经常跟年轻工程师说：「应力是原因，应变是结果。」应变的定义是：

ε = ΔL / L

同样，应变也有六个分量：三个正应变（εx, εy, εz）和三个剪应变（γxy, γyz, γzx）。

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

胡克定律，说白了就是「你拉我，我伸长；你压我，我缩短」。这个关系在弹性范围内是线性的。我刚开始做有限元时，总觉得胡克定律太简单了，后来才发现它的精妙之处。

一维胡克定律：

σ = E · ε

E 是弹性模量，也叫杨氏模量。不同材料的 E 值差别很大：

材料	弹性模量 E (GPa)	泊松比 ν
钢	200-210	0.27-0.30
铝	68-72	0.33-0.35
铜	110-130	0.34-0.36
钛	100-120	0.32-0.34

三维情况下，胡克定律就复杂多了。它把六个应力分量和六个应变分量联系起来。我建议你记住这个核心关系：

{σ} = [D] · {ε}

其中 [D] 就是弹性矩阵，也叫刚度矩阵。

2.3 弹性矩阵：有限元的核心密码

弹性矩阵 [D] 是连接应力和应变的桥梁。对于各向同性材料，它只依赖于两个参数：弹性模量 E 和泊松比 ν。

三维情况下的弹性矩阵形式：

[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 ]

这个矩阵看着复杂，但它的物理意义很清晰：对角线上的元素代表「拉压刚度」，非对角线上的元素代表「泊松效应」——也就是一个方向上的拉伸会引起另一个方向上的收缩。

2.4 平面应力与平面应变：简化但不简单

实际工程中，很多问题可以简化成二维。这就引出了两个经典假设：平面应力和平面应变。

平面应力问题：

适用于薄板结构，比如压力容器的壁面。假设厚度方向的应力为零（σz = 0），但应变不为零。我做过一个汽车覆盖件的分析，用的就是平面应力假设，结果跟实验吻合得很好。

平面应力的弹性矩阵简化为：

[D] = E / (1-ν²) ×
[ 1   ν   0 ]
[ ν   1   0 ]
[ 0   0  (1-ν)/2 ]

平面应变问题：

适用于长条形结构，比如水坝、隧道。假设厚度方向的应变为零（εz = 0），但应力不为零。

平面应变的弹性矩阵：

[D] = E / [(1+ν)(1-2ν)] ×
[ 1-ν   ν    0 ]
[  ν   1-ν   0 ]
[  0    0  (1-2ν)/2 ]

2.5 知识体系总览

下面这张图是我自己整理的弹性力学基础框架，帮你理清思路：

这张图把本章的核心逻辑串起来了。你从应力、应变出发，通过胡克定律和弹性矩阵，最终落到平面应力或平面应变的工程简化上。我个人觉得，理解了这个框架，后面学有限元就顺多了。

好了，这一章的内容就到这里。弹性力学基础虽然看起来理论性强，但它是有限元分析的「内功」。内功练好了，后面学什么招式都快。记住我说的：应力是原因，应变是结果，弹性矩阵是桥梁。下次见。

---

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