有限元分析数学原理

有限元分析的数学求解原理

梁弯曲问题的近似求解：

问题的提法为，求解挠度函数$$v(x)$$，使其满足以下方程和边界条件：
$$\begin{gathered} L(v(x))+\bar{b}=0 \\ BC(u):g_u(v(x))=0 \\ BC(u):g_p(v(x))=0 \end{gathered}$$
对于平面梁的弯曲问题，$L()=-\frac{d^4}{d{x}^4},\bar{b}=\bar{p}(x)$

1.加权残值法

试函数方法基本原理：先假定一满足边界条件的试函数，然后将其带入控制方程，通过使原方程的残差值最小来确定试函数中的待定系数。

将试函数带入微分控制方程，则一定存在残差，记为，
$$R=L(\hat{v}(x))+\bar{b}\ne 0$$
对于更一般的情形，设有一组满足边界条件的试函数${\varphi }_i(x)$​​，将其线性组合为新的试函数
$$\hat{v}(x)=c_1{\varphi }_1(x)+c_2{\varphi }_2(x)+\cdot \cdot \cdot +c_n{\varphi }_n(x)$$
将其带入原始方程，必定存在残差R，存在真实的$c_1,c_2,\cdots ,c_n$使得加权残差的积分为0。
$$\begin{gathered} {\int }_{\Omega }{w}_{t1}\cdot R(c_1,c_2,\cdots ,{\varphi }_1,{\varphi }_2,\cdots ,{\varphi }_n)\mathrm{d}\Omega =0 \\ {\int }_{\Omega }{w}_{t2}\cdot R(c_1,c_2,\cdots ,{\varphi }_1,{\varphi }_2,\cdots ,{\varphi }_n)\mathrm{d}\Omega =0 \\ ⋮ \\ {\int }_{\Omega }{w}_{t3}\cdot R(c_1,c_2,\cdots ,{\varphi }_1,{\varphi }_2,\cdots ,{\varphi }_n)\mathrm{d}\Omega =0 \end{gathered}$$
其中，$w_{t1},w_{t2}\cdots ,w_{tn}$为权函数。

Galerkin 加权残值法

如果将权函数$w_{t1},w_{t2}\cdots ,w_{tn}$取${\varphi }_1,{\varphi }_2,\cdots ,{\varphi }_n$,则该方法称为Galerkin加权残值方法。

对于简支梁弯曲问题，可将试函数$\hat{v}(x)$​取为，注意：此处 试函数需满足边界条件
$$\hat{v}(x)=c_1{\varphi }_1=c_1\cdot \sin \frac{\pi x}{l}$$
残差为
$$R=EI\frac{{\mathrm{d}}^4(c_1{\varphi }_1(x))}{\mathrm{d}{x}^4}-\overline{p_0}\ne 0$$
带入式（4），求解得
$$c_1=\frac{4l^4}{\pi EI}\overline{p_0}$$
则最终的结果为：
$$\hat{v}(x)=\frac{4l^4}{\pi EI}\overline{p_0}sin\frac{\pi x}{l}$$
若提高精度，使
$${\hat{v}}_2(x)=c_1{\varphi }_1+c_2{\varphi }_2=c_1\cdot \sin \frac{\pi x}{l}+c_2\cdot \sin \frac{3\pi x}{l}$$
带入，可以解出待定系数，
$$c_1=\frac{4l^4{\bar{p}}_0}{{\pi }^{5}EI},c_2=\frac{4l^4{\bar{p}}_0}{243{\pi }^{5}EI}$$

残值最小二乘法

同样，设有一个试函数满足所有的边界条件，将其带入原始方程，则必定有残差值R，存在$c_1,c_2\cdots ,c_n$​,使得残差的平方加权积分取极小值，即
$$min[E_{rr}={\int }_{\Omega }{w}_t\cdot R^2(c_1,c_2,\cdots ,{\varphi }_1,{\varphi }_2,\cdots ,{\varphi }_n)\mathrm{d}\Omega ]$$
其中，$w_t$为权函数，一般取1，将式（16）取极值，有
$$\begin{gathered} \frac{\partial E_{rr}}{\partial c_1}=0 \\ \frac{\partial E_{rr}}{\partial c_2}=0 \\ ⋮ \\ \frac{\partial E_{rr}}{\partial c_3}=0 \end{gathered}$$
通过求解可以得出试函数的待定系数。

特点

• 试函数需要满足所有的边界条件，即位移边界条件与力边界条件；

• 积分中试函数的最高阶导数较高（对于梁的弯曲问题，导数为4阶，对于一般弹性问题，导数为2阶），因此对试函数的连续性要求较高；

• 整个方法为计算一个几何域的积分

• 通过求取积分问题的最小值==（即误差最小）==，将原方程的求解变为线性方程组的求解。

2.虚功原理、最小势能原理

加权残值法对试函数连续性的要求过高，以下方法可以在只满足位移边界条件的同时，放松连续性的要求。

虚功原理

假想: 在平衡力系上作用微小的扰动==（不影响原平衡条件）==，且外力作用的位置产生了微小的位移，该位移称为虚位移

**弹性力学中的虚功原理：**在外力作用下处于平衡状态的变形体，当给物体以微小虚位移时，外力所做的总功等于物体的总虚应变能（即应力在由虚位移产生的虚应变上所做的功）。

​ 梁弯曲求解

设满足边界条件的位移场为：
$$\hat{v}(x)=c_1\sin \frac{\pi x}{l}$$
其中，$c_1$为待定系数，则虚位移场为
$$\delta \hat{v}(x)=\delta c_1\sin \frac{\pi x}{l}$$
则，该梁的虚应变能为(注意：虚应变能没有1/2)
$$\begin{gathered} \delta U={\int }_{\Omega }{\sigma }_x\delta {\epsilon }_x\mathrm{d}\Omega \\ ={\int }_0^l{\int }_{A}E\cdot {\epsilon }_x\delta {\epsilon }_x\mathrm{d}A\cdot \mathrm{d}x \\ =\frac{EIl}{2}(\frac{\pi }{l}{)}^4\cdot c_1\cdot \delta c_1 \end{gathered}$$
外力虚功为：
$$\delta W={\int }_0^l{\bar{p}}_0\delta \hat{v}\mathrm{d}x$$
由虚功原理，有$\delta W=\delta U$

可以解出$c_1$​.

最小势能原理

在所有满足位移边界条件的许可位移场中，真实的位移场${\hat{u}}_i$使得物体的总势能取得最小值，即
$$min[\Pi ({\hat{u}}_i)=U-W]$$
其中$U=\frac{1}{2}{\int }_{\Omega }{\sigma }_x\cdot {\epsilon }_x\cdot \mathrm{d}\Omega$; $W={\int }_0^l{\bar{p}}_0\cdot \hat{v}(x)\mathrm{d}x$​,

利用偏导为零可以求得待定系数。