一维变分不等式问题的自适应有限元分析新探  

A NEW APPROACH TO SELF-ADAPTIVE FEM FOR ONEDIMENSIONAL VARIATIONAL INEQUALITY PROBLEMS

作　　者：袁驷[1] 刘泽洲[1] 邢沁妍[1] 

机构地区：[1]清华大学土木工程系土木工程安全与耐久教育部重点实验室,北京100084

出　　处：《工程力学》2015年第7期11-16,共6页Engineering Mechanics

基　　金：国家自然科学基金项目(51378293;51078199;50678093;50278046);长江学者和创新团队发展计划项目(IRT00736)

摘　　要：结构工程中的弹性薄膜接触和杆件弹塑性扭转等问题是典型的变分不等式问题,对其高效精确求解,特别是满足给定精度要求下的自适应求解,是挑战性课题。该文作者新近成功实现了一维变分不等式问题的自适应有限元分析,该文对此进展作一报道。对于变分不等式的有限元求解,该文提出区域二分法和C检验技术,极大提升了松弛迭代的收敛速度,一般4次~5次线性解即可得到收敛的有限元解答,进而采用作者提出的EEP（单元能量投影）超收敛公式计算超收敛解答,用其检验误差并指导网格细分,逐步得到堪称为数值精确解的解答,亦即得到按照最大模度量逐点满足精度要求的解答。该文给出的数值算例表明所提出的算法具有高效、可靠、精确的优良特性。Elastic membrane contact problems and elasto-plastic torsion problems are typical variational inequality problems, and the accurate and efficient analyses of these problems, especially in a self-adaptive manner to satisfy the user-preset error tolerance, are still of great challenges. This paper reports a recently developed approach to the self-adaptive finite element （FE） analysis of one-dimensional （1D） variational inequality problems. To obtain the FE solution rapidly based on variational inequalities, the present paper proposed two techniques, i.e. bisection bounding and C-check technique, which greatly accelerate the convergence rate of the conventional relaxation iteration, and in average, it suffices to implement only 4 to 5 steps of linear FE analysis. With the converged FE solution obtained, the super-convergent solution by EEP （Element Energy Projection） method, also developed by the first author＇s research team, is calculated to estimate the errors and then to guide a mesh refinement with final results being virtually ＇numerically exact＇, i.e. satisfying the error tolerance by the maximum norm. The numerical examples presented show that the proposed approach is efficient, accurate and reliable.

关 键 词：变分不等式 有限元 单元能量投影 超收敛 自适应 

分 类 号：TU311.4[建筑科学—结构工程]