“…The price paid is that the inclusions have to discretized into some finite elements. The domain integral equation method was also used to investigate the propagation of the crack near the inclusion (Dong et al, 2003b) and 3D inclusion-crack interaction (Dong et al, 2002).…”