椭圆型偏微分方程的三角形单元有限元的数值解法  

Triangle finite element numerical solution of elliptic partial differential equations

作　　者：俞辉[1] 杨臻豪 奉继明 瞿少成[2] 

机构地区：[1]三峡大学理学院,湖北宜昌443000 [2]华中师范大学电子信息工程系,武汉430079

出　　处：《华中师范大学学报（自然科学版）》2016年第4期489-495,共7页Journal of Central China Normal University：Natural Sciences

基　　金：国家自然科学基金项目(61273183;61374028;61304162)

摘　　要：本文以椭圆型偏微分方程问题为背景研究三角形单元的有限元方法.随着计算机图形学的发展,三角形单元划分取得了巨大的成就,可以得到质量非常好的三角形单元,进而提高偏微分方程的有限元方法的数值解的精度.本文采用节点增量算法,对问题区域进行三角形单元划分,得到的三角形单元满足Delaunay条件,再对三角形单元的所有节点采用自适应编号,最后运用三角形单元的有限元方法得到椭圆型偏微分方程的数值解.通过数值实验,得出相比传统的三角形单元的有限元方法,本文的三角形单元的有限元方法减小了舍入误差,提高了计算精度.In the present work,the finite element method for triangular element is studied based on the elliptic partial differential equation.With the development of computer graphics,the triangle element division has made great achievements.High-quality triangular elements are able to be generated,which improved the accuracy of numerical solutions for partial differential equations.Here,the triangular element satisfying the Delaunay conditions is calculated through unit division of the region with the node incremental algorithm.All nodes of the triangular element are adopted with adaptive numbers and the numerical solution of elliptic partial differential equations is obtained utilizing the finite element method of triangular elements.Finally,through numerical experiments,the finite element method is compared with the traditional one.Results indicate that the finite element method of triangular element reduces the rounding error and elevates the calculation precision.

关 键 词：Delaunay三角形单元 自适应编号 舍入误差 

分 类 号：O175.25[理学—数学]