A numerical scheme, based on the concept of domain integral method, is
developed for evaluation of J(2) in two-dimensional anisotropic linea
r elastic solids. The scheme is designed to be in conjunction with fin
ite element solutions. Although the method is verified to be patch-ind
ependent, in a modified sense, the near tip region is always included
in the calculation due to the nature of J(2). Special attention is add
ressed on appropriate modelling of the singular behavior around the cr
ack tip area. The numerical examples show that J(2) can be calculated
with good accuracy without using;any particular singular elements. The
order of singularity of the strain energy density jump along the crac
k surfaces in the near tip region is also properly estimated. Copyrigh
t (C) 1996 Elsevier Science Ltd.