📚 课程目录第2章有限元理论基础

2、有限元理论基础：弹性力学基本方程、虚功原理、最小势能原理

各位工程师朋友，咱们今天聊聊有限元背后的“硬核”理论。说实话，我刚入行那会儿，也觉得这些公式离实际工程很远。但后来我发现，不懂这些，你连结果误差出在哪都说不清楚。

有限元分析，说白了就是用数学去模拟物理。你想想看，一个零件受力变形，它内部到底发生了什么？这就要用到弹性力学的基本方程。而虚功原理和最小势能原理，则是把物理问题转化成我们能求解的代数方程的关键桥梁。

这部分是基础中的基础。我习惯把弹性力学问题看成是“找三个量”：应力、应变、位移。它们之间由三类方程联系着。

物体内部任意一个微元体，受到的力必须平衡。这就像你站在地上，重力向下，地面支持力向上，你才能不动。在三维问题中，平衡方程是这样的：

∂σx/∂x + ∂τxy/∂y + ∂τxz/∂z + fx = 0
∂τyx/∂x + ∂σy/∂y + ∂τyz/∂z + fy = 0
∂τzx/∂x + ∂τzy/∂y + ∂σz/∂z + fz = 0

其中，σ 是正应力，τ 是剪应力，f 是体积力（比如重力）。

我的经验： 检查平衡方程是否满足，是判断网格是否足够细的一个快速方法。我曾经遇到一个悬臂梁模型，粗网格下支反力总和与外力差了好几个百分点，加密网格后立马就对上了。

物体变形了，我们得知道它“变”了多少。应变就是描述这种变形的量。小变形情况下，几何方程是线性的：

εx = ∂u/∂x
εy = ∂v/∂y
εz = ∂w/∂z
γxy = ∂u/∂y + ∂v/∂x
γyz = ∂v/∂z + ∂w/∂y
γzx = ∂w/∂x + ∂u/∂z

这里 u, v, w 是三个方向的位移，ε 是正应变，γ 是剪应变。

注意： 这个方程只适用于小变形。如果你做的是大变形分析（比如橡胶密封圈压缩），几何方程就复杂多了，得用格林应变这类非线性描述。

这是材料的“脾气”。对于最常见的线弹性材料，就是胡克定律：

{σ} = [D] {ε}

其中 [D] 是弹性矩阵，由材料的弹性模量 E 和泊松比 ν 决定。对于各向同性材料，它长这样：

这三类方程加起来，一共15个方程（3个平衡，6个几何，6个物理），对应15个未知量（6个应力，6个应变，3个位移）。理论上可以求解，但实际边界条件复杂，直接解偏微分方程组几乎不可能。所以，我们需要更聪明的办法——虚功原理。

虚功原理，听起来很玄乎，其实核心思想很简单：如果一个物体处于平衡状态，那么对于任意微小的、满足边界条件的虚位移，外力做的虚功等于物体内部产生的虚应变能。

用公式表达就是：

δW_ext = δU_int

其中，δW_ext 是外力虚功，δU_int 是虚应变能。

为什么说它重要？因为它把复杂的微分方程变成了一个积分方程。你想想看，我们不需要去管每个点上的力是否平衡，只需要保证整个系统在“能量”上是平衡的。这大大简化了问题。

  核心理解： 虚功原理是有限元法的“灵魂”。它让我们从“力”的视角，切换到了“能量”的视角。我个人觉得，理解了这个切换，你就抓住了有限元的一半。

最小势能原理是虚功原理的一个特例，适用于弹性体。它说：在所有可能的位移中，真实的位移会使系统的总势能取最小值。

总势能 Π 等于应变能 U 减去外力势能 V：

Π = U - V

这个原理有什么用？它给了我们一个“搜索”真实解的方法。我们假设一个位移场（用形函数和节点位移表示），然后让总势能对这个位移场的参数（节点位移）求导，并令其等于零。这样，我们就得到了有限元的核心方程：

[K] {u} = {F}

你看，[K] 是刚度矩阵，{u} 是节点位移向量，{F} 是节点力向量。整个有限元求解，本质上就是在解这个线性方程组。

避坑指南： 我曾经犯过一个错误，在分析一个过盈配合问题时，没有正确施加接触边界条件，导致模型的总势能没有达到最小值，结果算出来的应力分布完全不对。后来检查才发现，是接触刚度设置得太低了。所以，边界条件和接触设置，直接决定了你的模型是否满足最小势能原理的前提。

为了让你更直观地理解这三者之间的关系，我画了一张图。你可以把它当作本章的“地图”。

从这张图你可以看到，弹性力学方程是“原材料”，描述了物理本质。虚功原理是“加工方法”，把微分方程转化成积分形式。最小势能原理是“优化准则”，帮我们找到那个最真实的解。这三者层层递进，最终导向了有限元的核心方程。

嗯，理论部分就讲到这里。别怕，这些公式你不需要死记硬背，关键是理解它们背后的物理意义和逻辑关系。下次当你看到软件报错或者结果异常时，能想起今天聊的这些，知道从哪个方向去排查，那这章就没白学。