ANumericalMethod for Solving Elasticity Equations with Interfaces  

作　　者：Songming Hou Zhilin Li Liqun Wang Wei Wang 

机构地区：[1]Department of Mathematics and Statistics,Louisiana Tech University,Ruston,LA,71272,USA [2]Department of Mathematics,North Carolina State University,Raleigh,NC 27695,USA [3]Nanjing Normal University,China

出　　处：《Communications in Computational Physics》2012年第7期595-612,共18页计算物理通讯（英文）

基　　金：supported by Louisiana Board of Regents RCS Grant No.LEQSF(2008-11)-RD-A-18;The second author is partially supported by theAROgrants 56349MA MA,the AFSOR grant FA9550-09-1-0520;the NSF grant DMS-0911434,the NIH grant 096195-01 and CNFS 11071123.

摘　　要：Solving elasticity equationswith interfaces is a challenging problemformost existing methods.Nonetheless,it has wide applications in engineering and science.An accurate and efficient method is desired.In this paper,an efficient non-traditional finite element method with non-body-fitting grids is proposed to solve elasticity equations with interfaces.The main idea is to choose the test function basis to be the standard finite element basis independent of the interface and to choose the solution basis to be piecewise linear satisfying the jump conditions across the interface.The resulting linear system of equations is shown to be positive definite under certain assumptions.Numerical experiments show that thismethod is second order accurate in the L¥norm for piecewise smooth solutions.More than 1.5th order accuracy is observed for solution with singularity(second derivative blows up)on the sharp-edged interface corner.

关 键 词：Elasticity equations non-body fitted mesh finite element method jump condition 

分 类 号：O17[理学—数学]