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

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

1. 有限元位移解的下限性

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

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

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

$$\tilde{u}^TP \le u^TP$$

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

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

attachments-2019-04-QREb3FJH5caec3847cff5.jpeg

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

2. 应力结果的精度

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

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

  • 引入位移边界条件;
  • 求解方程组得到各结点的位移a;

$$Ka=P$$

公式1

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

$$\varepsilon=Lu=Ba^e$$

公式2

$$\sigma=D\varepsilon=DBa^e$$

公式3

例如对于平面问题,

$$\varepsilon=\lbrack\varepsilon_x \varepsilon_y \varepsilon_z\rbrack$$

$$ L=\begin{bmatrix}{ \partial}/{ \partial{x} } & 0 & { \partial }/{ \partial{y}}\\0 & { \partial }/{ \partial{x} } & { \partial}/{ \partial{y} }\end{bmatrix}^T$$

$$N=\begin{bmatrix} N_1& 0& N_ 2& 0& N_ 3&0&...\\ 0 &N_1&0&N_2&0&N_3 &... \end{bmatrix}$$

$$\large B=\large L \large N$$

公式4

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

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

  • 单元内部一般不满足平衡方程;
  • 单元与单元交界面上应力一般不连续;
  • 力的边界上一般也不满足力边界条件。

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

3. 应力解的震荡性质

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

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

attachments-2019-04-TAuQDmWf5caec57981a88.jpeg

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

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

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

4. 减缩积分和完全积分单元

完全积分是指当单元具有规则形状时,所用的高斯积分点的数目足以对单元刚度矩阵中的多项式进行精确积分。

减缩积分是指对单元刚度阵进行积分时,所用的高斯积分点数低于精确积分的要求,一般是按照下式来确定高斯积分点数,

$$n=p-m+1 $$

公式5

式中,n为高斯积分点数,p是插值函数中完全多项式的方次,m是微分算子L中的导数阶次(对于用于弹塑性力学分析的实体单元,m=1)。

按照式(5)确定积分点数的减缩积分单元比完全积分单元在每个方向上减少了一个积分点。

对于大部分商用有限元软件,二维实体单元中的线性与二次四边形单元有减缩积分形式,积分点数分别为1*1和2*2。三维实体单元中的线性和二次六面体单元有减缩积分形式,积分点数分别为1*1*1和2*2*2。三角形单元、四面体单元和楔形单元一般没有减缩积分形式。

attachments-2019-04-jEbzExRi5caec5ccbfde3.png

图3 完全积分和减缩积分单元

因为降低了刚度矩阵的积分精度,减缩积分单元看似会影响结果的准确性。但实际计算表明,采用减缩积分往往可以取得比精确积分更好的精度,主要原因如下:

  • 分享于 2019-04-11 12:11
  • 阅读 ( 447 )

[版权声明] :本文系网友分享,仅以非商业性的交流和科研为目的,并不代表本网赞同其观点和对其真实性负责。如涉及作品内容、版权和其它问题,请在30日内与本网( friends@stuch.cn )联系!我们将协调给予处理。阅读原文(请登录)..

0 条评论

请先 登录 后评论
猜猜我是谁
甜也 -研究生

9
提问
7
回答
4
文章
注册推广