孙浩涵, 袁驷. 基于EEP超收敛解的自适应有限元法特性分析[J]. 工程力学, 2019, 36(2): 17-25. DOI: 10.6052/j.issn.1000-4750.2018.10.0439
引用本文: 孙浩涵, 袁驷. 基于EEP超收敛解的自适应有限元法特性分析[J]. 工程力学, 2019, 36(2): 17-25. DOI: 10.6052/j.issn.1000-4750.2018.10.0439
SUN Hao-han, YUAN Si. PERFORMANCE OF THE ADAPTIVE FINITE ELEMENT METHOD BASED ON THE ELEMENT-ENERGY-PROJECTION TECHNIQUE[J]. Engineering Mechanics, 2019, 36(2): 17-25. DOI: 10.6052/j.issn.1000-4750.2018.10.0439
Citation: SUN Hao-han, YUAN Si. PERFORMANCE OF THE ADAPTIVE FINITE ELEMENT METHOD BASED ON THE ELEMENT-ENERGY-PROJECTION TECHNIQUE[J]. Engineering Mechanics, 2019, 36(2): 17-25. DOI: 10.6052/j.issn.1000-4750.2018.10.0439

基于EEP超收敛解的自适应有限元法特性分析

PERFORMANCE OF THE ADAPTIVE FINITE ELEMENT METHOD BASED ON THE ELEMENT-ENERGY-PROJECTION TECHNIQUE

  • 摘要: 基于EEP (单元能量投影)超收敛计算的自适应有限元法,已对一系列问题取得成功,但其自适应特性尚缺乏相关研究。该文以二阶常微分方程为模型问题,同时考察基于EEP和SPR (超收敛分片恢复)超收敛解的自适应分析方法,与有限元最优网格进行了比较分析,进而提出反映自适应有限元收敛特性的估计式,并给出了自适应收敛率β的定义。该文给出的数值试验表明:采用m次单元,对于解答光滑的问题,SPR法与EEP法均可有效用于自适应求解,其位移可按最大模获得m+1的自适应收敛率;对于奇异因子为α(<1)的奇异问题,SPR法失效,而基于EEP法的自适应求解,其位移按最大模可获得m+α的自适应收敛率,远高于α的常规有限元收敛率。

     

    Abstract: The adaptive finite element method (AFEM) based on the element energy projection (EEP) technique has succeeded in solving a wide range of problems, while few studies have been done on its adaptive characteristics. Taking the second-order ordinary differential equation as the model problem, the adaptive methods based on both the EEP and super-convergent patch recovery (SPR) solutions were studied. The meshes generated by the two adaptive methods were compared with the optimum mesh. Furthermore, an estimate formula reflecting the characteristics of AFEM was proposed with an adaptive convergence rate β. Numerical experiments show that, with the element of degree m, both SPR and EEP can be well applied for the smooth problem and the convergence rate achieved m+1. For the singular problem with a singularity factor α(<1), SPR failed while EEP-based adaptive procedure gained the convergence rate of m+α, a much higher rate than that of the conventional finite element method α.

     

/

返回文章
返回