1引言

作为一名汽车结构CAE工程师，使用有限元法进行应力分析是我多年的日常工作。但是关于有限元应力结果的一些技术点一直未能吃透，查到的相关文献也不多，文献内容也存在很多不尽不实之处。最近一段时间，抽空重新读了一遍王勖成老师的《有限单元法》，也查阅了一些其它资料，简单整理出了一些关于应力结果的笔记，再加上一点自己的理解，与大家分享交流。

2有限元位移解的下限性

利用最小位能原理求得的位移近似解所对应的弹性变形能是精确解变形能的下限，即位移近似解在总体上偏小，模型偏于刚硬。利用最小余能原理得到的应力近似解所对应的弹性余能是精确解余能的上限，即应力近似解在总体上偏大，结构的计算模型偏于柔软。

我们常用的结构单元大都是位移元，以位移为未知量，基于最小位能原理建立。位移元得到的位移解具有下限性质，在给定的载荷之下，计算模型的变形比实际要小。当单元网格分割得越来越细时，位移数值解将由下方收敛于精确解，即得到真实解的下界。

需要注意的是，位移数值解并不是在每个点上都小于精确解，数值解只是总体上小于精确解，更准确的说法应该是外载荷在数值解上做的功小于在精确解上做的功，即

ũTP≤uTP（1）

上式中，ũ和u分别为位移数值解和精确解，P为外载荷。

位移解的下界性质从物理上非常好理解。连续体具有无限多个自由度，可以有无数种复杂的变形模式。划分单元后，位移场用有限个结点的位移和对应的形函数来表示，即用一些简单变形模式来逼近实际变形，如图1所示。这就意味着连续体的变形受到了约束和限制，即刚度较实际增加了。由刚度方程可知，在外力相同的情况下，所求得的位移近似解将在总体上偏小。

图1 用特定的形函数逼近准确解

3应力结果的精度

使用位移有限元法进行结构分析时，未知的场函数是结构位移。利用最小位能原理建立的求解方程是系统的平衡方程，求解方程得到的是各结点的位移，但实际工程问题往往更关注结构应力分布。

位移有限元法求解应力的基本步骤如下：

1. 引入位移边界条件；

2. 求解方程组得到各结点的位移a；

   Ka=P                 （2） 

3. 根据单元各结点位移，通过导数运算求得应变和应力。

ε=Lu=Baeσ=Dε=DBae（3）

例如对于平面问题，

应变矩阵B是对插值函数求导得到的矩阵，每求导一次，插值多项式的次数就降低一次。所以通过导数运算得到的应力解的精度较位移解降低了一阶。例如线性单元的应力是近似均布的，二次单元的应力是近似线性分布的，其精度都比位移解低（虽然应力分布函数还包含了一些高次项，但这些高次项都是非完全项，并不能提高应力精度）。

有限元应力解的近似性表现在：

1. 单元内部一般不满足平衡方程；

2. 单元与单元交界面上应力一般不连续；

3. 力的边界上一般也不满足力边界条件。

只有单元尺寸无限趋近于0时，即自由度数趋近无穷时，才能精确满足平衡方程、力边界条件和单元交界面上的应力连续性。在单元数量有限时，这些条件只能近似满足，除非实际应力变化的阶次等于或低于所用单元的应力分布函数阶次。

4应力解的震荡性质

位移有限元法从力学上解释是求位移变分所引起的应变能为极小值的问题，从数学上解释是求解应力近似解与精确解差值的加权最小二乘问题。

与位移结果不同，位移元的应力结果并没有下限性质。应力结果是精确应力���加权最小二乘意义上的近似解，应力近似解必然在精确解上下震荡；并且在某些点上，近似解恰好等于精确解，即单元内存在最佳应力点，如图2所示。

图2 有限元应力解的震荡性质

我们在有限元强度分析中，经常发现在应力集中部位的应力解低于精确解，所以有人误以为位移有限元法的应力解也有下限性质。但真正原因是单元数目有限，且应力近似解的阶次低，无**确描述应力的剧烈变化。对于高应力梯度区域，有限元解通常给出的是一个比实际平滑的结果。如果结构上某部位的应力突然上升，该部位的应力解将低于精确解，但如果某个部位应力突然下降，该处的应力解反而会高于精确解。

位移元的应力解没有下限性质从物理上也容易解释，有限元相当于在结构内部增加了约束，即增加了整体刚度。在外载荷相同的情况下，刚度增加能导致位移数值解偏小，但解得的应力应该在总体上保持不变，否则就无法满足平衡方程。

点击完整阅读全文 

关于有限元应力结果精度的几点讨论

相关阅读推荐：​​​​​​​

有限元分析工程应用实例详解