工程力学

• 基本方法 • 上一篇    下一篇

弹性力学问题自适应有限元及其局部多重网格法

刘春梅1, 肖映雄2, 舒适1, 钟柳强3   

  1. 1. 湘潭大学数学与计算科学学院,湖南,湘潭 411105;
    2. 湘潭大学土木工程与力学学院,湖南,湘潭 411105;
    3. 华南师范大学数学科学学院,广东,广州 510631
  • 收稿日期:2010-12-17 修回日期:2011-07-21 出版日期:2012-09-25 发布日期:2012-09-25
  • 通讯作者: 肖映雄(1970―),男,湖南城步人,教授,博士,从事计算固体力学研究(E-mail: xyx610xyx@yahoo.com.cn). E-mail:xyx610xyx@yahoo.com.cn
  • 作者简介:刘春梅(1981―),女,山西五台人,博士生,从事偏微分方程数值解研究(E-mail: liuchunmei8080@yahoo.com.cn);舒 适(1962―),男,湖南双峰人,教授,博士,博导,从事偏微分方程数值解及多重网格法研究(E-mail: shushi@xtu.edu.cn);钟柳强(1980―),男,广东兴宁人,博士,从事偏微分方程数值解研究(E-mail: zlq1980@gmail.com).
  • 基金资助:

    国家自然科学基金重大研究计划重点支持项目(91130002);国家自然科学基金项目(10972191);长江学者和创新团队发展计划项目(IRT1179);湖南省研究生创新基金项目(CX2010B245,CX2010B246)

ADAPTIVE FINITE ELEMENT METHOD AND LOCAL MULTIGRID METHOD FOR ELASTICITY PROBLEMS

LIU Chun-mei1, XIAO Ying-xiong2, SHU Shi1, ZHONG Liu-qiang3   

  1. 1. School of Mathematics and Computational Science, Xiangtan University, Xiangtan, Hunan 411105, China;
    2. Civil Engineering and Mechanics College, Xiangtan University, Xiangtan, Hunan 411105, China;
    3. School of Mathematical Sciences, South China Normal University, Guangzhou, Guangdong 510631, China
  • Received:2010-12-17 Revised:2011-07-21 Online:2012-09-25 Published:2012-09-25
  • Contact: xiao ying-xiong E-mail:xyx610xyx@yahoo.com.cn

摘要:

针对平面弹性力学问题,利用最新顶点二分法,设计了一种不需要标记振荡项和加密单元不需要满足“内节点”性质的自适应有限元法;利用自适应加密过程中每层网格上只有局部单元需要加密这一特性,设计了一种基于局部松弛的多重网格法.数值实验结果表明:该文设计的自适应有限元法具有一致收敛性和拟最优计算复杂度,基于局部松弛的多重网格法对求解弹性力学问题自适应网格下的有限元方程具有很好的计算效率和鲁棒性.

关键词: 弹性力学问题, 自适应有限元, 拟最优复杂性, 局部松弛, 多重网格法

Abstract:

In this paper, an adaptive finite element (AFEM) method is designed by using the newest vertex bisection for linear elasticity problems in two dimensions. This method marks exclusively according to the error estimator without special treatment of oscillation and performs a minimal element refinement without the interior node property. Furthermore, a type of multigrid method based on the local relaxation is applied to the AFEM discrete systems by using the special properties during refinement. The results of various numerical experiments are shown that the proposed AFEM method is uniformly convergent and has quasi-optimal numerical complexity. The resulting multigrid method is much more robust and efficient in CPU times than the usual multigrid methods.

Key words: elasticity problems, adaptive finite element method, quasi-optimal complexity, local relaxation, multigrid method

中图分类号: 

  • O343.3
[1] Senturia S, Aluru N, White J. Simulating the behavior ofMEMS devices: Computational methods and needs [J].IEEE Computational Science and Engineering, 1997, 4:30-43.  



[2] Senturia S, Harris R, Johnson B, et al. A computer-aideddesign system for micro-electromechanical systems [J].Journal of Micro-electromechanical Systems, 1992, 1:3-13.



[3] Brenner C, Li Yengsung. Linear finite element methodsfor planar linear elasticity [J]. Mathematics ofComputation, 1992, 59: 321-338.  



[4] 陈竹昌, 王建华, 王卫中. 自适应多层网格有限元求解应力集中问题[J]. 同济大学学报, 1994, 22(3): 203-208.Chen Zhuchang, Wang Jianhua, Wang Weizhong.Adaptive multigrid FEM for stress concentration [J].Journal of Tongji University, 1994, 22(3): 203-208. (inChinese)



[5] 梁力, 林韵梅. 有限元网格修正的自适应分析及其应用[J]. 工程力学, 1995, 12(2): 109-118.Liang Li, Lin Yunmei. Adaptive mesh refinement offinite element method and its application [J]. EngineeringMechanics, 1995, 12(2): 109-118. (in Chinese)



[6] 王建华. 线弹性有限元的自适应加密与多层网格法求解[J]. 河海大学学报, 1994, 22(3): 16-22.Wang Jianhua. Adaptive refinement and multigridsolution for linear finite element method [J]. Journal ofHohai University, 1994, 22(3): 16-22. (in Chinese)



[7] 王建华, 殷宗泽, 赵维炳. 自适应多层网格有限元网格生成器研制[J]. 计算结构力学与其应用, 1995, 12(1):86-92.Wang Jianhua, Yin Zongze, Zhao Weibing.Implementation of the mesh generator for adaptivemultigrid finite element method [J]. ComputationalStructural Mechanics and Applications, 1995, 12(1):86-92. (in Chinese)



[8] Cai Z Q, Korsawe J, Starke G. An adaptive least squaresmixed finite element method for the stress displacementformulation of linear elasticity [J]. Numerical Methodsfor Partial Differential Equations, 2005, 21(1): 132-148.  



[9] Whiler T P. Locking free adaptive discontinuous galerkinFEM for linear elasticity problem [J]. Mathematics ofComputation, 2006, 75(255): 1087-1102.  



[10] Chen L, Zhang C S. A coarsening algorithm on adaptivegrids by newest vertex bisection and its applications [J].Journal of Computational Mathematics, 2010, 28(6):767-789.



[11] Verfürth R. A review of a posteriori error estimationtechniques for elasticity problem [J]. Computer Methodsin Applie Mechanics and Engineering, 1999, 176: 419-440.  



[12] Carstensen C, Dolzmann G, Funken S A, Helm D S.Locking-free adaptive mixed finite element methods inlinear elasticity [J]. Computer Methods in ApplieMechanics and Engineering, 2000, 190(13): 1701-1718.  



[13] Lonsing M, Verfürth R. A posteriori error estimators formixed finite element methods in linear elasticity [J].Numerische Mathematik, 2004, 97: 757-778.  



[14] Mitchell W F. Optimal multilevel iterative methods foradaptive grids [J]. SIAM Journal on Scientific andStatistical Computing, 1992, 13: 146-167.  



[15] Wu H J, Chen Z M. Uniform convergence of multigridV-cycle on adaptively refined finite element meshes forsecond order elliptic problems [J]. Science in ChinaSeries A: Mathematics, 2006, 49(10): 1405-1429.  



[16] Carstemsen C. Convergence of adaptive finite elementmethods in computations mechanics [J]. AppliedMathematics and Computation, 2009, 59(9): 2119-2130.



[17] Cascon J, Kreuzer C, Nochetto R, Siebert K.Quasi-optimal convergence rate for an adaptive finiteelement method [J]. SIAM Journal on NumericalAnalysis, 2008, 46: 2524-2550.



[18] Dörfler W. A convergent adaptive algorithm forpoisson’s equation [J]. SIAM Journal on NumericalAnalysis, 1996, 33(3): 1106-1124.
[1] 蔡勇, 李光耀, 王琥. 基于多重网格法和GPU并行计算的大规模壳结构快速计算方法[J]. 工程力学, 2014, 31(5): 20-26.
[2] 张 肖;张建海;周 钟;饶宏玲. 锦屏一级高拱坝大垫座稳定性研究[J]. 工程力学, 2010, 27(增刊I): 232-235,.
[3] 袁 驷;叶康生;王 珂. 平面曲梁面内自由振动分析的自适应有限元法[J]. 工程力学, 2009, 26(增刊Ⅱ): 126-132.
[4] 杨强;吴浩;周维垣. 大坝有限元分析应力取值的研究[J]. 工程力学, 2006, 23(1): 69-73.
[5] 谢学斌;肖映雄;潘长良;舒适. 代数多重网格法在岩体力学有限元分析中的应用[J]. 工程力学, 2005, 22(5): 165-170.
[6] 林晓辉;吴京荣. 求解线接触弹流润滑完全数值解的多重网格复合直接迭代法[J]. 工程力学, 1998, 15(3): 126-132,.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 原 园;徐颖强;吕国志;朱贤飞. 齿轮啮合过程中安定状态残余应力的数值方法研究[J]. 工程力学, 2008, 25(10): 0 -211, .
[2] 邢德进;李忠献. 应用SMA智能阻尼器的结构模糊控制[J]. 工程力学, 2008, 25(10): 0 -228, .
[3] 周小平;杨海清;张永兴. 有限宽偏心裂纹板在裂纹面受两对集中拉力作用时裂纹线的弹塑性解析解[J]. 工程力学, 2008, 25(1): 0 -027 .
[4] 龚耀清;包世华. 超高层建筑空间巨型框架自由振动计算的新方法[J]. 工程力学, 2008, 25(10): 0 -140 .
[5] 刘金兴;邓守春;张 晶;梁乃刚. 颗粒复合材料断裂的梁网格模型[J]. 工程力学, 2008, 25(10): 0 -037 .
[6] 郎风超;邢永明;朱 静. 应用纳米压痕技术研究表面纳米化后316L 不锈钢力学性能[J]. 工程力学, 2008, 25(10): 0 -071 .
[7] 郭小刚;刘人怀;曾 娜;金 星. 子结构位移迭代法修正软管空间形态[J]. 工程力学, 2008, 25(10): 0 -032 .
[8] 邢静忠;柳春图. 线弹性土壤中埋设悬跨管道的屈曲分析[J]. 工程力学, 2008, 25(10): 0 -075 .
[9] 刘祥庆;刘晶波. 基于纤维模型的拱形断面地铁车站结构弹塑性地震反应时程分析[J]. 工程力学, 2008, 25(10): 0 -157 .
[10] 郝庆多;王言磊;侯吉林;欧进萍;. GFRP带肋筋粘结性能试验研究[J]. 工程力学, 2008, 25(10): 0 -165, .
X

近日,本刊多次接到来电,称有不法网站冒充《工程力学》杂志官网,并向投稿人收取高额费用。在此,我们郑重申明:

1.《工程力学》官方网站是本刊唯一的投稿渠道(原网站已停用),《工程力学》所有刊载论文必须经本刊官方网站的在线投稿审稿系统完成评审。我们不接受邮件投稿,也不通过任何中介或编辑收费组稿。

2.《工程力学》在稿件符合投稿条件并接收后会发出接收通知,请作者在接到版面费或审稿费通知时,仔细检查收款人是否为“《工程力学》杂志社”,千万不要汇款给任何的个人账号。请广大读者、作者相互转告,广为宣传!如有疑问,请来电咨询:010-62788648。

感谢大家多年来对《工程力学》的支持与厚爱,欢迎继续关注我们!

《工程力学》杂志社

2018年11月15日