具有最佳超收敛阶的EEP法计算格式:I算法公式

对一维C0问题的高次有限元后处理中超收敛计算的EEP(单元能量投影)法提出改进的最佳超收敛计算格式,即用m次单元对足够光滑问题的有限元解答,采用该格式计算的任意一点的位移和应力都可以达到h2m阶的最佳超收敛结果.整个工作分为3个部分,分别给出算法公式、数值算例和数学证明.该文是系列工作的第一部分,针对高次单元提出了凝聚形函数的概念,并证明了相关的逼近定理和等价定理,在此基础上给出了具体的算法公式.     

二维自由振动问题的自适应有限元分析初探

自由振动反映结构动力特性,是抗震分析和结构设计的重要基础。近年来,基于单元能量投影（EEP）法的自适应有限元分析已在一系列线弹性及非线性问题中取得成功,而有限元线法（FEMOL）自适应分析在二维自由振动问题中的应用也被证实是有效的。在此基础上,该文进一步提出二维自由振动问题的自适应有限元分析方法。通过将特征值问题线性化,合理引入二维线性问题的EEP超收敛计算和自适应求解技术,该法可得到满足精度要求的自振频率和按最大模度量满足用户给定误差限的振型。该文以弹性薄膜为例,介绍了这一进展,并给出数值算例以表明该方法的有效性和可靠性。     

一个高效的一维有限元自适应求解的新方案 第十三届全国结构工程学术大会特邀报告

基于新近提出的一维有限元后处理超收敛算法——单元能量投影(EEP)法,将有限元自适应求解问题转化为对超收敛解答的自适应分段多项式插值问题,一步便可获得最优的有限元网格划分,在该网格上再次进行有限元计算,即可获得满足用户给定的误差限的有限元解答。该法简单实用、快速高效,是一个颇具优势和潜力的自适应方法。文中以二阶常微分方程模型问题为例,对该法的形成思路和实施策略做一介绍,并给出有代表性的数值算例用以展示该法的优良性能和效果。     

二维有限元线法超收敛解答计算的EEP法

有限元线法(FEMOL)是一种优良的半解析、半离散方法,但其解答存在解析方向和离散方向的精度不相称的弱点。本文提出将二维有限元线法比拟为广义一维问题的概念,遂可将新近提出的一维有限元超收敛计算的单元能量投影(EEP)法推广到二维有限元线法分析中。经有限元线法后处理中EEP超收敛计算而获得的解答,继承和保留了一维有限元中的出色表现,不但使任意一点的位移和应力的解答在两个方向具有相当的精度,而且都具有超收敛性质。文中以二维Poisson方程问题为例,具体给出了有限元线法EEP超收敛的公式,并给出了数值算例,用以表明本法的可行性和有效性。     

具有最佳超收敛阶的EEP法计算格式:Ⅱ数值算例

对一维C0问题的高次有限元后处理中超收敛计算的EEP(单元能量投影)法提出改进的最佳超收敛计算格式,即用m次单元对足够光滑问题的有限元解答,采用该格式计算的任意一点的位移和应力都可以达到h2m阶的最佳超收敛结果。整个工作分为3个部分,分别给出算法公式、数值算例和数学证明。该文是系列工作的第二部分,给出实施算法和数值算例,用以验证理论公式的有效性和正确性。     

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

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

一维有限元后处理超收敛解答计算的EEP法

提出一维有限元法后处理中超收敛解答的一种自然合理的算法,称为单元能量投影法(EEP)。理论分析和数值算例表明,提出的方法简便易行、行之有效、效果显著;此外,还有一些颇合人意的优点,如:任一点的应力和位移的误差与结点位移的误差具有相同的收敛阶(m次单元可达mh2阶)、结点两边单元各自算出的应力自动平衡、自由端点的应力自动为精确值等。     

具有最佳超收敛阶的Galerkin有限元EEP法计算格式

对二阶非自伴问题的一维Galerkin有限元法提出其后处理超收敛计算的EEP(单元能量投影)法改进的最佳超收敛计算格式,即用m次单元对足够光滑问题的Galerkin有限元解答,采用该格式计算的任一点的位移和应力都可以达到h2m阶的最佳超收敛结果。该文首先针对高次单元提出了凝聚试探形函数和凝聚检验形函数的概念,证明了相关的逼近定理和等价定理,然后给出了具体的算法公式。最后给出了一系列典型的数值算例用以验证这种最新的EEP法改进格式确实能够使位移和导数逐点达到最佳收敛阶。     

The multi-projection method for weakly singular Fredholm integral equations of the second kind

In this paper, we propose a multi-projection method and its re-iterated algorithm for solving weakly singular Fredholm integral equations of the second kind. We apply our methods to Petrov–Galerkin versions to establish excellent superconvergence results, and we illustrate our theoretical results with a numerical example.