自然边界元与有限元求解平面弹性问题的耦合法.docVIP

自然边界元与有限元求解平面弹性问题的 耦合法 臧彤,赵慧明,杨敏 (中国矿业大学力学与建筑工程学院,江苏 徐州 221008 摘要:为了更充分地利用有限元与自然边界元各自的优点,并尽量减少由于方法的局限性造成的在计算量及计算精度上的不足。本文引入了有限元与自然边界元耦合的方法,在文中简单介绍了自然边界元法及其与有限元的耦合的原理,通过设置含有重叠区域的圆形人工边界,实现自然边界元法与有限元法的耦合,并把该方法应用到无界区域上的实际算例中,从计算结果的比较中可以看出自然边界元与有限元耦合算法的收敛速度更快,充分体现了耦合法在解决无界区域问题上的优越性。 关键词:自然边界元法;有限元法;无界区域;耦合法;D-N迭代 中图分类号:TB112 A coupling method of natural boundary element and finite element for the planar elastic problem Zang Tong, Zhao Huiming, Yang Min (School of Mechanices and Civil Engineering, China University of Mining and Technology, JiangSu XuZhou 221008 Abstract: In order to make better use of finite element and boundary element’s the merits, and to minimize the deficiencies in calculating the amount and accuracy on the deficiencies caused by method limitations, this article introduced the natural boundary element method and its coupling with the finite element process, and described simply the natural boundary element method and its coupling with the finite element process in the text. This paper gave a briefing on the natural boundary element iterative realization of D-N by setting a circular artificial boundary containing the overlap region, and applied the method to the real unbounded example. By the comparing of results from the different calculation methods, the precision and algorithm converges of results by using coupling method were more good. The analysis fully reflected the superiority of coupling method to solve unbounded regional problems. Keywords:NBEM; FEM; unbounded regional; direct coupling method; D-N iteration 0引言 自然边界元法与有限元法自创立以来,在众多领域都取得了瞩目的成就。但它们也都存在着一定的不足,有限元方法在处理无界问题时有着本质的困难,而自然边界元法则由于一般区域的格林函数难以求得,使得它的应用范围较窄,但二者基于相同的变分原理,这为两种方法的耦合奠定了坚实的理论基础。自然边界元方法,是由Green函数和Green公式出发,将微分方程边值问题归化为边界上强奇异积分方程(或称为超奇异积分方程的一种数值计算方法,在一定程度上它保持了原边值问题的一些有用的性质,保证了解的唯一性和稳定性【1】。另外从数值计算的角度,与一般边界元方法相比,自然边界元方法使刚度矩阵的计算量也大为减少。耦合法既吸取了有限元能适应较任意的有限区域的优点,又克服了自然边界归化对区域一些要求的限制,其总体刚度矩阵正是自然边界元法及有限元法得到的刚度矩阵之和【2】。 作者简介:臧彤,(1985-,女,硕士研究生,主要研究方向:计算方法. E-mail: cumt_zt@163.com 这与其他类型的边界