基于改进位移模式的有限元超收敛算法研究

【摘要】：超收敛计算的研究是近年来有限元领域研究中的热点与难点之一。通过对基于改进位移模式的积分形式与基于常规位移模式伽辽金方程的比较，导出了伽辽金方程精确成立的条件。在以此条件作为位移模式分析的理论依据前提下，本文提出了改进位移模式，即将高阶有限元解的位移模式用常规有限元解的位移模式表示，用两种位移模式之和构造新的位移模式，结合多尺度方法的思想和超收敛计算的解析公式，提出了一种全新的前处理超收敛的计算方法。该方法无需任何人为的磨光，对一维有限元法的线性单元，本文结点和单元的位移、导数都达到了h4阶的超收敛精度。应力跨单元自动平衡，可精确满足自然边界条件，且实施简便，计算量增加很小。 有限元线法是一个不断发展的课题，国内外对有限元线法超收敛的计算研究得较少。本文利用有限元线法的半解析性，将一维有限元中获得全面成功的改进位移模式的超收敛算法推广至二维有限元线法，取得了良好的效果。 全文主要工作如下： 1、通过基于改进位移模式与基于常规位移模式伽辽金方程的比较，导出了伽辽金方程精确成立的条件。论证了位移模式与泡函数的相关性。提出了一次元与泡函数结合的单元方案。利用精确成立的条件，通过量级分析，保留泡函数主要的影响项，避免了求解泡函数的解析表达式。为位移模式的分析提供了理论依据。 2、以一维C0问题为模型问题，提出有限元法的基于改进位移模式的前处理超收敛算法。在原有试函数的基础上，增加了高阶试函数，使得单元内平衡方程的残差减少，从而达到提高精度的目标。对于近似单元，根据单元内部平衡条件，导出单元上任一点的位移和导数的超收敛解的计算公式。基于伽辽金方法，采用积分形式推导了单元刚度矩阵。 3、将基于改进位移模式的前处理超收敛算法成功推广到一维有限元法其它问题，包括(1)一维C1问题；(2)二阶非自伴两点边值问题Galerkin有限元问题，亦即将基于改进位移模式的前处理超收敛算法推广到了非自伴算子问题而不仅仅是自伴算子的问题。（3）一维n阶问题的超收敛算法，本部分工作为该法广泛应用于一般一维问题的有限元法的超收敛计算打下了良好的基础。 4、将基于改进位移模式的前处理超收敛算法成功推广应用到二维有限元线法的Poisson方程问题，基于线性形函数，采用变分形式推导了有限元线法求解的修正的常微分方程组。算例结果表明：结点和单元内的位移、导数的收敛精度得到了极大的提高。 由于本文得出的应力和位移是逐点超收敛的，因而有望在此基础上发展出不同于目前常规的误差估计和自适应求解方法。