“…The enrichment can also be intrinsic, based on the recent work by Fries and Belytschko [8]. In this paper, we focus on the extrinsic partition of unity enrichment and in general, the field variables are approximated by [1,4,9,10,6,11,7,12]: u h (x) = I∈N fem N I (x)q I + enrichment functions (1) where N I (x) are standard finite element shape functions, q I are nodal variables associated with node I. XFEM, one of the aforementioned partition of unity methods, was successfully applied for crack propagation and other fields in computational physics [13,14,15,16,17,18,19,20,21,22] and recently open source XFEM codes were released to help the development of the method [23] and numerical implementation and efficiency aspects were studied [24]. XFEM is quite a robust and popular method which is now used for industrial problems [25] and under implementation by leading computational software companies.…”