Two boundary elements have been proposed for simulation of variable order s
ingularities in two dimensions. The first can model the variable order stra
in singularity near a crack tip and the second can simulate both the strain
and traction. These elements can be easily incorporated in a standard boun
dary element computer programme. The elements are useful for computation of
stress intensity factors (SIFs) in fracture mechanics. To improve the accu
racy of such computations further, a modified crack closure integral (MCCI)
based method for mechanical loading are presented. Examples of mode I and
mixed mode crack problems are examined to illustrate the performance of the
proposed elements and the MCCI based calculations. The effects of order of
Gauss quadrature associated with such elements on the accuracy of the SIFs
are also reported. (C) 2000 Elsevier Science Ltd. All rights reserved.