摘要
本文采用了面向对象的程序设计方法,以C++作为编程语言,构造了一个局部自适应有限元P加密程序。研究内容涉及从Serendipity函数的构造、单元的形成、自升阶算法以及误差估计和判断。最后采用C++的一系列类及成员函数实现了这些功能,并构造了相应的输入输出模块,使其成为一个完整的自适应有限元分析程序。文中给出的计算结果表明,在P加密有限元分析程序中,误差随着阶数的升高而下降。而且,局部P加密的区域也集中在应力奇异的部位。这验证了文中的理论推导及程序实现的正确性。
In this paper, C++ langrage with Object-oriented features was used design a P-version adaptive FEM program, which contains the modules of high-order and Serendipity shape function, high-order element, adaptive course and error estimate. In the end, besides a series of class and their functions, the input and output modules were constructed, which make the program become whole adaptive FEM software. Examples show that the errors decrease with the increasing of p refinement. The local P-version a-daptive zone is upgraded concentrate in stress-singularity zone. The abouve results show that the formulate derivation and FEM program implementation are correct.
出处
《力学季刊》
CSCD
北大核心
2002年第4期463-470,共8页
Chinese Quarterly of Mechanics
基金
国家自然科学基金(59879012)