3次元弾性体解析
理论与物理
三维弹性理论的基础
老师，FEM的结构分析说到底是在求解“三维弹性理论”对吧？
控制方程
请告诉我三维弹性理论的控制方程。
有三个基本方程。
1. 平衡方程（力的平衡）
2. 应变-位移关系（协调条件）
3. 本构关系（胡克定律）
6个应力分量（$\sigma_x, \sigma_y, \sigma_z, \tau_{xy}, \tau_{yz}, \tau_{xz}$）、6个应变分量、3个位移分量。总共15个未知数对应15个方程，对吧。
理解得非常完美。将这15个方程以位移为唯一未知数汇总起来，就是纳维方程（拉梅-纳维方程）：
这里 $\lambda, \mu$ 是拉梅常数。$\mu = G$（剪切弹性模量），$\lambda = \frac{E\nu}{(1+\nu)(1-2\nu)}$。
各向同性弹性体的D矩阵
FEM中使用的 $[D]$ 矩阵是怎样的？
各向同性弹性体的三维本构关系（Voigt记法）：
这个6×6的矩阵就是所有的基础啊。平面应力和平面应变的 $[D]$ 就是它的缩减版。
是的。从三维的 $[D]$ 出发，平面应力是在 $\sigma_z = 0$ 条件下缩减，平面应变是在 $\varepsilon_z = 0$ 条件下缩减。壳和梁也是基于各自的假设从三维的 $[D]$ 推导出来的。
各向异性材料
对于非各向同性材料呢？
一般各向异性弹性体的 $[D]$ 拥有21个独立常数（6×6对称矩阵）。但实用中以下特殊情况较多：
| 材料对称性 | 独立常数 | 示例 |
|---|---|---|
| 各向同性 | 2 ($E, \nu$) | 钢、铝 |
| 横观各向同性 | 5 | 单向增强CFRP单层、沉积土 |
| 正交各向异性 | 9 | 木材、织物CFRP |
| 一般各向异性 | 21 | 晶体 |
CFRP（碳纤维复合材料）是横观各向同性的吗？
单向材料（UD材）的单层是横观各向同性的。纤维方向及其正交面对称。多层叠合后，整体上会呈现正交各向异性或更复杂的对称性。
三维分析不可或缺的场景
哪些情况下无法使用二维近似，必须进行三维分析？
“不用三维就解不了”的问题意外地多呢。
随着计算机性能提升，三维分析已变得普遍。但用三维求解并不一定意味着正确。如果网格质量、边界条件、材料模型不正确，三维分析也会产生垃圾结果。
总结
我来整理一下三维弹性体分析的理论。
要点：
- 三个基本方程 — 平衡、协调、本构关系。总计15未知数15方程
- $[D]$ 矩阵是所有的基础 — 二维单元是其缩减版
- 各向同性为2常数，各向异性最多21常数
- 三维分析不可或缺的场景 — 3D应力集中、厚度方向应力、接触、复杂形状
- 三维 ≠ 一定准确 — 输入质量决定结果
FEM的所有单元理论都归结于三维弹性理论。理解了这一点，似乎就能自然地判断该使用哪种单元了。
正是如此。三维弹性理论是“主干”，各单元类型是“枝叶”。理解了主干，枝叶的选择自然就明白了。
3D弹性体的控制方程
三维弹性体的平衡方程由纳维（Navier）和柯西（Cauchy）在1820年代独立推导。连接9分量应力张量σij和6分量应变张量εij的广义胡克定律（σ=Cε，C为四阶弹性张量）最多拥有C11〜C66共21个独立分量。对于各向同性材料，可归结为拉梅常数λ和μ两个参数，并可转换为杨氏模量E=μ(3λ+2μ)/(λ+μ)。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出的经历吗？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是“缓慢施力故加速度可忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想要恢复的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长”的特性就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高=强度高”是不对的。刚度是“不易变形”，强度是“不易破坏”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（重力等）和表面力 $f_s$（压力、接触力等）。可以这样想——桥上卡车的重量是“作用在整个内容物上的力”（体积力），轮胎压路面的力是“只作用在表面的力”（表面力）。风压、水压、螺栓紧固力…全都是外力。这里容易犯的错误：弄错载荷方向。本想“拉伸”却成了“压缩”——听起来像笑话，但在3D空间中坐标系旋转时确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他的弦。声音会一直持续吗？不，会逐渐变小。这是因为振动能量通过空气阻力和弦的内部摩擦变成了热。汽车的减震器也是同样原理——故意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，所以设置适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入mm时，载荷·弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意区分工程应变与对数应变（大变形时） |
| 弹性模量 $E$ | Pa | 钢: 约210 GPa，铝: 约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
三维实体单元
三维实体单元有哪些种类？
基本形状有三种：六面体（hex）、四面体（tet）、五面体（wedge/prism）。
| 单元 | 1次 | 2次 | 精度 | 网格生成 |
|---|---|---|---|---|
| 四面体 | TET4（4节点） | TET10（10节点） | TET4: 低 / TET10: 高 | 自动网格容易 |
| 六面体 | HEX8（8节点） | HEX20（20节点） | HEX8: 中 / HEX20: 非常高 | 自动网格困难 |
| 五面体 | WEDGE6 | WEDGE15 | 中〜高 | 用于hex与tet的过渡 |
TET4精度低吗？
TET4是常应变单元（CST的三维版）。单元内应变为常数，因此无法表现应力梯度。TET4在实际工作中不应使用。使用TET10（二次四面体）则精度足够。
但是自动网格生成的话，TET4更容易生成吧。
确实如此，但TET4要得到准确的
なった
詳しく
報告