第二章 有限元基础理论回顾：弹性力学基本方程、虚功原理、最小势能原理、有限元求解流程

各位工程师朋友，大家好。欢迎来到《光机结构有限元分析实战演练》的第二讲。

这一章，咱们聊聊理论。别急着皱眉，我知道很多人一听「理论」就犯困。但说实话，有限元这东西，你如果只懂操作软件，不懂背后的力学原理，那就像开车只会踩油门，不懂刹车和方向盘——迟早要出事。

我个人习惯是：先花半小时把核心理论捋一遍，后面建模和调试时能省下三天时间。咱们今天就把弹性力学、虚功原理、最小势能原理和求解流程这四块硬骨头，用工程师的语言啃下来。

2.1 弹性力学基本方程

弹性力学，说白了就是研究「物体在外力作用下怎么变形、内部应力怎么分布」的学问。在光机结构里，镜片受压、支架受热、螺钉受拉，全得靠它来算。

核心方程就三组，我管它叫「三大金刚」：

2.1.1 平衡方程

物体内部任意一个微元体，受力必须平衡。用数学写出来就是：

∂σ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是体力（比如重力）。

我在项目中遇到过一个案例：某光学平台在振动测试时局部应力超标，一查发现是平衡方程没满足——网格太粗，微元体上的力没传过去。后来加密网格，问题就解决了。嗯，这里要注意：平衡方程是有限元求解的底线，不满足它，结果就是错的。

2.1.2 几何方程

描述应变和位移的关系。小变形情况下：

ε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是三个方向的位移，ε是正应变，γ是剪应变。

你想想看，如果镜片支架的位移算错了，那光路偏移量就全错了。所以几何方程是连接「位移场」和「应变场」的桥梁。

2.1.3 物理方程（本构关系）

也就是大家熟悉的胡克定律。对于各向同性材料：

εx = (1/E)[σx - ν(σy + σz)]
εy = (1/E)[σy - ν(σx + σz)]
εz = (1/E)[σz - ν(σx + σy)]
γxy = τxy/G
γyz = τyz/G
γzx = τzx/G

E是弹性模量，ν是泊松比，G是剪切模量。这三个参数，做光机结构的人必须烂熟于心。

2.2 虚功原理

虚功原理，名字听着玄乎，其实道理很简单：如果一个物体处于平衡状态，那么对于任意微小的虚位移，外力做的虚功等于内力做的虚功。

用公式表达：

∫V σij δεij dV = ∫V fi δui dV + ∫S ti δui dS

左边是内力虚功（应力在虚应变上做功），右边是外力虚功（体力和面力在虚位移上做功）。

我曾经用虚功原理推导过一个复杂镜组的刚度矩阵。当时软件算出来的结果总是不收敛，我手动用虚功原理验算了一遍，发现是边界条件加错了。虚功原理就像一面照妖镜，能帮你发现模型里的妖魔鬼怪。

为什么会这样？因为虚功原理不依赖于材料的本构关系，它只要求平衡。所以它特别适合用来推导有限元方程——你不需要知道材料具体怎么变形的，只要知道它平衡就行。

2.3 最小势能原理

最小势能原理是虚功原理的「升级版」。它说：在所有可能的位移场中，真实位移场使系统的总势能取最小值。

总势能Π = 应变能U - 外力势能W

用数学表达：

Π = (1/2)∫V σij εij dV - ∫V fi ui dV - ∫S ti ui dS

对Π求变分，令δΠ=0，就得到了平衡方程。

我个人觉得，最小势能原理是有限元法最优雅的地方。它把复杂的微分方程问题，转化成了一个求极值的代数问题。你想想看，本来要解偏微分方程组，现在只需要找一组位移让势能最小——这不就是计算机最擅长的事吗？

2.4 有限元求解流程

好了，理论铺垫完毕，咱们来看看有限元到底是怎么求解的。我把流程总结成六个步骤：

1. 结构离散化：把连续体划分成有限个单元，单元之间通过节点连接。
2. 单元分析：对每个单元，建立形函数，推导单元刚度矩阵。
3. 整体组装：把所有单元的刚度矩阵组装成整体刚度矩阵。
4. 施加边界条件：固定约束、载荷、温度等。
5. 求解方程组：解 [K]{u} = {F}，得到节点位移。
6. 后处理：由位移计算应变、应力，进行结果评估。

下面这张图，是我自己画的一个流程图，把整个求解过程串起来了：

这个流程，说白了就是：把连续问题离散化，把微分方程变成代数方程，然后让计算机去解。

我刚开始学有限元时，总觉得这六个步骤是割裂的。后来做项目多了才发现，每一步都环环相扣。比如网格划分得不好，单元分析就会出错；边界条件加错了，求解出来的位移就是天方夜谭。

2.5 小结

这一章咱们把有限元的基础理论过了一遍。弹性力学三大方程是根基，虚功原理和最小势能原理是桥梁，求解流程是实操指南。

你可能会问：这些理论在实际光机结构分析中到底怎么用？别急，下一章咱们就开始动手建模了。到时候你会发现，今天讲的这些理论，每一个都会在软件操作中体现出来。

记住一句话：理论是地图，软件是交通工具。没有地图，你开再好的车也会迷路。

---

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