摘要
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.
基金
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.