Timoshenko梁单元超收敛结点应力的EEP法计算

将新近提出的单元能量投影(ElementEnergyProjection,简称EEP)法应用于Timoshenko梁单元的超收敛结点应力计算· 根据单元投影定理具体推导了一般单元的计算公式,并对两个有代表性的单元给出了数值算例· 分析和算例表明,EEP法对于解答是向量函数(即常微分方程组)的问题具有同样优良的表现,不仅能给出与结点位移精度同阶、同量级的超收敛结点应力。

基于应力超收敛恢复技术的广义特征值问题后验误差估计

根据有限元解的超收敛特性提出了一种基于应力超收敛恢复技术的广义特征值问题后验误差估计。通过对单元内的应力超收敛点以及相邻单元的应力超收敛点进行插值或外推处理,得到单元内其它点处更高精度的应力解。通过高精度的应力值可以得到结构处理后改进的势能。将改进的势能代入瑞利商,最终得到比原始有限元解更高精度的特征解。将后处理特征解作为“准精确解”代替误差估计因子中未知的精确解,实现后验误差估计过程。数值计算结果表明,所提出的后验误差估计是渐进精确的,因此可作为结构广义特征值问题自适应有限元方法的误差估计因子。

一维有限元后处理的EEP法的数学分析

利用一维投影型插值与有限元超收敛基本估计,对一类两点边值问题,严格证明了袁驷等人由单元能量投影(EEP)法获得的节点恢复导数,当有限元空间的次数不超过4时,具有最佳阶超收敛.理论分析圆满地解释了已有的数值结果.