A Reduced Crouzeix–Raviart Immersed Finite Element Method for Elasticity Problems with Interfaces

2020 ◽  
Vol 20 (3) ◽  
pp. 501-516
Author(s):  
Gwanghyun Jo ◽  
Do Young Kwak

AbstractThe purpose of this paper is to develop a reduced Crouzeix–Raviart immersed finite element method (RCRIFEM) for two-dimensional elasticity problems with interface, which is based on the Kouhia–Stenberg finite element method (Kouhia et al. 1995) and Crouzeix–Raviart IFEM (CRIFEM) (Kwak et al. 2017). We use a {P_{1}}-conforming like element for one of the components of the displacement vector, and a {P_{1}}-nonconforming like element for the other component. The number of degrees of freedom of our scheme is reduced to two thirds of CRIFEM. Furthermore, we can choose penalty parameters independent of the Poisson ratio. One of the penalty parameters depends on Lamé’s second constant μ, and the other penalty parameter is independent of both μ and λ. We prove the optimal order error estimates in piecewise {H^{1}}-norm, which is independent of the Poisson ratio. Numerical experiments show optimal order of convergence both in {L^{2}} and piecewise {H^{1}}-norms for all problems including nearly incompressible cases.

2017 ◽  
Vol 9 (2) ◽  
pp. 407-428 ◽  
Author(s):  
Daehyeon Kyeong ◽  
Do Young Kwak

AbstractIn this paper, we propose a finite element method for the elasticity problems which have displacement discontinuity along the material interface using uniform grids. We modify the immersed finite element method introduced recently for the computation of interface problems having homogeneous jumps [20, 22]. Since the interface is allowed to cut through the element, we modify the standard Crouzeix-Raviart basis functions so that along the interface, the normal stress is continuous and the jump of the displacement vector is proportional to the normal stress. We construct the broken piecewise linear basis functions which are uniquely determined by these conditions. The unknowns are only associated with the edges of element, except the intersection points. Thus our scheme has fewer degrees of freedom than most of the XFEM type of methods in the existing literature [1,8,13]. Finally, we present numerical results which show optimal orders of convergence rates.


Sign in / Sign up

Export Citation Format

Share Document