当前位置：首页 > news >正文

有限元技巧核心原理与学习路径：从一维基础到多维拓展（七步流程）

news 2026/5/12 17:04:46

有限元学习七步流程（系统优化版）

本文先从一维杆的核心力学关系入手，求解其应力、应变；随后介绍位移法的求解思路；在此基础上，引入试函数法，通过虚功原理与最小势能原理建立问题的弱形式；最后，对有限元法的基本计算流程进行总结。

如若希望直接了解整体流程，可跳过前面的PPT内容，从 “① 理论入门：掌握物理本质” 开始观看。

已知: 通过给定的材料属性（如杨氏模量 $E$ 和截面面积 $A$ ）和外力 $F$ ，计算结构在受力下的位移 $u$ 、应变 $\epsilon$ 和应力 $\sigma$ 。

1、传统材料力学求解方法

当结构较为复杂时，采用传统的材料力学方法求解将变得十分繁琐，因此引入了位移法来进行分析。

2、位移解法

接下来将介绍有限元分析的基本解题流程。

基础流程总结

有限元特点

1️⃣ 力学法的局限

最早的结构分析方法是力学法（位移解法）。
：就是它的基本思路
通过建立平衡方程 + 几何方程 + 本构关系，求解结构的位移、应变和应力分布。

这种方法理论严谨、结果精确，但困难在于——
一旦结构稍微复杂，便会出现高阶偏微分方程。
微分方程不仅难解，而且容易出错，手算几乎不可能。
因此，虽然力学法最精确，却不具备工程可操作性。

2️⃣ 引入能量原理：摆脱微分方程

为了避免直接求解微分方程，人们引入了虚功原理和最小势能原理。
这两个原理是等价的能量法，它们允许我们绕开微分方程，借助积分求解。

核心思想是：

先假设位移函数，再用能量平衡确定未知参数。就是不再直接解微分方程，而

例如：

你能够猜位移是线性的： $u(x) = a_0 + a_1 x$

或者猜是正弦的： $u(x) = b_0 + b_1 \sin x$

只要假设满足（或接近）边界条件，就能通过能量原理得到结构响应。
整个过程只涉及积分，不再需要解微分方程！

这就是所谓的：

“从强形式（微分方程）到弱形式（积分方程）”的转变。
具体推导过程可参考链接

虚功原理

最小势能原理

例子
一维

二维

三维

① 理论入门：掌握物理本质

建立力学三大基本关系（以一维杆为例）：

平衡方程：
$\frac{d\sigma}{dx} + f = 0$
几何关系：
$\varepsilon = \frac{du}{dx}$
本构关系：
$\sigma = E\,\varepsilon$

材料假设：线性弹性、均质、各向同性。
目标：明确应力–应变–位移之间的物理因果关系。

② 数学建模：位移与应变的假设

假设线性位移场：
$c\,x$
边界条件：
$u (0) = 0$
根据力平衡和边界条件：
$\frac{P}{E A}$
因此：
$\varepsilon = \frac{du}{dx} = \frac{P}{E A}, \quad \sigma = E\varepsilon = \frac{P}{A}$

③ 能量法与虚功原理：建立变分关系

采用势能极值原理：
$\Pi(u) = U - W$

平衡条件由：
$\delta \Pi = 0 \quad \Rightarrow \quad \delta U = \delta W$

对于杆件问题：
$\frac{P}{E A}x, \quad \varepsilon = \frac{P}{E A}, \quad \sigma = \frac{P}{A}$

⚙️ 虚功原理推导过程

(1) 虚功方程

$\delta U = \delta W$

(2) 应变能变化

$\delta W = P(x)\,\delta u(x)$

$\delta U = \delta \left( \frac{1}{2} \int_{\Omega} \sigma(x)\,\varepsilon(x)\, dV \right)$

(3) 试函数

$c\,x \quad \Rightarrow \quad \delta u(x) = \delta c \cdot x$

$\sigma(x) = E\,\varepsilon(x) \quad \Rightarrow \quad \delta\sigma(x) = E \cdot \delta c$

$\varepsilon(x) = \frac{du}{dx} \quad \Rightarrow \quad \delta\varepsilon(x) = \delta c$

(5) 平衡结果

$\frac{P}{E A}, \quad \varepsilon(x) = \frac{P}{E A}, \quad \sigma(x) = \frac{P}{A}, \quad u(x) = \frac{P}{E A}x$

⚙️ 最小势能原理推导

(1) 最小势能原理

$\min_{u(x) \in BC(u)} \big[\, \Pi(u) = U - W \,\big]$

(2) 环境应变能

$\frac{1}{2} \int_{\Omega} \sigma(x)\,\varepsilon(x)\, dV$

(3) 系统外力势能

$\cdot u(x = l)$

(4) 势能函数

$\Pi(u) = \frac{1}{2} E\, c^2\, A\, l - P\, c\, l$

(5) 取极值条件

$\frac{\partial \Pi}{\partial c} = 0$

(6) 试函数与结果

$c\,x, \quad \varepsilon(x) = c, \quad c = \frac{P}{E A}$

同样得到：
$\varepsilon(x) = \frac{P}{E A}, \quad \sigma(x) = \frac{P}{A}, \quad u(x) = \frac{P}{E A}\, x$

④ 弱形式与有限元基础

由强形式（微分方程）通过分部积分得到弱形式：
$\int_0^L EA \frac{du}{dx}\frac{dv}{dx}dx = \int_0^L f\,v\,dx + [\bar{T}v]_L$

边界条件：

本质边界条件：直接作用于位移，如 $u (0) = 0$ ；
自然边界条件：由弱式自动包含，如 $\sigma(L) = \bar{T}$ 。

⑤ 有限元离散与矩阵方程

单元内位移近似：
$u(\xi) = N_1 u_1 + N_2 u_2, \quad N_1 = 1 - \xi, \; N_2 = \xi$
单元刚度矩阵：
$K_e = \int_0^L B^T E A B\,dx = \frac{E A}{L} \begin{bmatrix} 1 & -1 \\ -1 & 1 \end{bmatrix}$
全局组装：按节点连接关系叠加 $K_e$ 形成整体刚度矩阵 $K$ 。

⑥ 数值计算与验证

通过求解线性方程：
$K U = F$
得到节点位移。
计算应变与应力：
$\varepsilon = B U_e, \quad \sigma = E \varepsilon$
与解析解对比，验证有限元解的正确性。

⑦ 知识融合与高维扩展

（1）从一维到二维的推广

二维问题常见于平面应力与平面应变状态：

$\quad v = v(x, y)$

几何方程（应变–位移关系）：
$\begin{aligned} \varepsilon_x &= \frac{\partial u}{\partial x}, \\ \varepsilon_y &= \frac{\partial v}{\partial y}, \\ \gamma_{xy} &= \frac{\partial u}{\partial y} + \frac{\partial v}{\partial x} \end{aligned}$

矩阵形式：

在这里插入图片描述

平面应力本构关系：
$\mathbf{D}_{\text{plane-stress}} = \frac{E}{1-\nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \frac{1-\nu}{2} \end{bmatrix}$

平面应变本构关系：
$\mathbf{D}_{\text{plane-strain}} = \frac{E(1-\nu)}{(1+\nu)(1-2\nu)} \begin{bmatrix} 1 & \frac{\nu}{1-\nu} & 0 \\ \frac{\nu}{1-\nu} & 1 & 0 \\ 0 & 0 & \frac{1-2\nu}{2(1-\nu)} \end{bmatrix}$