单元应力vs节点应力
单元的应力由该单元全部节点变形与坐标计算得到,由于相邻单元仅在相邻的节点与边界上变形连续,而非相邻的边界上并不一定相等,导致两相邻单元在相邻节点上计算的应力不连续,一般采用绕节点平均或应力磨平等方法进行处理。1、节点的应力值有很多处理方法,包括绕节点平均与应力磨平等等,不能简单认为周围各单元在同一节点上应力的均值即为该点的应力值。2、减小网格尺寸可以使不同单元在同一节点处的应力值接近,但是这个方法太笨了,对复杂模型的具体应用上存在问题,会很大程度上增加机时,不建议采取这种方法,最好能利用应力磨平等方法来求得近似解答,各种有限元书籍上都有相关的介绍。
ABAQUS中应力存储在高斯积分点而不是节点的原因!!(看来还要好好挖掘一下有限元理论啊)By Bravo(姜琛 Jesse)可能很多人不知道,ABAQUS计算后的应力其实是放在积分点的,为什么这跟我们所想的不一样,应力应该跟位移一样放在节点上呢?
其实这是由于数值计算造成的,因为积分是可以提高精度的,微分是会降低精度(泰勒展开一下,微分一下,会将多项式的次数降低一次,也就降低一阶的精度)。所以为什么说有限元算出的应力一般精度较低,如果在软件 中不开应力磨平,后处理时可以看见相邻单元的应力有一个突变,即不连续的。(当然针对二次单元不适用,因为位移求导得到应力是一阶线性的,而非一次的常数,所以应力也是连续的)
但是根据有限元的应力近似解性质(见王勖成,有限单元法,P168),应力和应变近似解一定是在精确解上下震荡的,但是在某些点上这些解刚好和精确解相等,也即最佳应力点。根据以应力为自变量的最小位能原理,其中是有积分的,一般都采用高斯数值积分。由高斯数值积分的性质,必然在积分点处是精确值。
至此你懂了ABAQUS为什么要这么存储了。然后通过积分点的应力值,采用单元形函数可以插值得到单元内任意点的应力,当然这比节点插值的精度要高!
所以说有的时候,学软件的时候,是可以再review一下理论的