平面弹性问题自适应有限元方法的收敛性分析-openrepositoryof.pdfVIP

　 ， ３５ ９ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１４年９月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３５，Ｎｏ．９，Ｓｅｐ．１５，２０１４ 文章编号：１０００⁃０８８７（２０１４）０９⁃０９６９⁃１０ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 平面弹性问题自适应有限元 方法的收敛性分析∗ １ ２ ３ ４ 刘春梅 ，　 钟柳强 ，　 舒　 适 ，　 肖映雄 （１． 湖南科技学院数学与计算科学系计算数学研究所，湖南 永州４２５１９９； ２． 华南师范大学 数学科学学院，广州 ５１０６３１； ３． 湘潭大学 数学与计算科学学院，湖南 湘潭４１１１０５； ４． 湘潭大学 土木工程与力学学院，湖南 湘潭４１１１０５） 摘要：　 针对平面弹性问题，首先采用基于最新顶点二分法的网格加密方法，给出一种不需要标记 振荡项和加密单元、不需要满足“内节点”性质的自适应有限元方法．其次，通过对各层网格上解函 数和误差指示子的分析，利用相邻网格层上解函数的正交性、解函数和真解函数的能量误差的上 界估计、相邻网格层上误差指示子的近似压缩性等结果，从理论上严格证明了该自适应有限元方 法是收敛的．最后数值实验验证了该自适应有限元方法是收敛的和鲁棒的． 关　 键　 词：　 平面弹性问题；　 自适应有限元方法；　 收敛性 中图分类号：　 Ｏ２４１．８；Ｏ２４２　 　 　 文献标志码：　 Ａ ｄｏｉ：１０．３８７９／ ｊ．ｉｓｓｎ．１０００⁃０８８７．２０１４．０９．００３ 引　 　 言 弹性问题的求解在科学计算和工程技术上都是很重要的问题，如对微电子机械系统中的 ［１］ 多体问题进行数值模拟时，常常需要进行弹性分析 ．但由于问题本身及求解区域的复杂性， ［２］ 一般很难得到它们的解析解，有限元法是目前进行弹性分析时最常用的数值方法之一 ．在实 际应用中，有许多因素使得弹性问题产生强奇性解，如物理区域是一个非光滑的几何区域（如 凹角或裂缝等）、材料系数不连续、应力集中等．此时在进行数值求解，若采用一致加密网格，则 将会导致自由度的指数增长，从而引起计算代价的剧烈增加．自适应有限元方法可以依据解的 性态来进行网格点的自动分布，由此有效克服由网格一致加密所导致的自由度过度增长的求 解困难，从而实现利用最小的计算量获得最大的计算精度，因而在现代工程计算中得到广泛的 ［３⁃６］ 应用 ． ∗ 收稿日期：　 ２０１４⁃０１⁃２０；修订日期：　 ２０１４⁃０６⁃０５ 基金项目：　 湖南省自然科学基金（１４ＪＪ３１３５）； 国家自然科学基金；全国博士学位论文作 者专项资金（２０１２１２）；广东省高等学校优秀青年教师培养计划（Ｙｑ２０１３０５４）；广州市珠 江科技新星项目（２０１３Ｊ２２０００６３） 作者简介：　 刘春梅（１９８１—），女，山西人，讲师，博士（Ｅ⁃ｍａｉｌ：ｌｉｕｃｈｕｎｍｅｉ０６２９＠１６３．ｃｏｍ）； 舒适（１９６２—），男，湖南人，教授，博士，博士生导师（通讯作者． Ｅ⁃ｍａｉｌ：ｓｈｕｓｈｉ＠ｘｔｕ． ｅｄｕ．ｃｎ）． ９６９ ９７０ 刘 春 梅　 　 　 钟 柳 强　 　 　 舒　 　 适　 　 　 肖映雄