In this paper, the steps in the derivation of the analytical expressions, which could be used to evaluate directly stress intensity factors for cracks from the numerical results of boundary integral equation axisymmetric fracture mechanics analysis are reviewed. These expressions have been found previously to provide an efficient means by which stress intensity factors can be obtained without unnecessarily refined mesh discretisations. They include those which could be used to analyse interface cracks between dissimilar materials. Such crack problems present modelling difficulties when using conventional procedures for obtaining the fracture parameters because of the oscillatorily singular stresses in the vicinity of the interface crack tip. Three axisymmetric problems, each involving a crack at an elastic inclusion which, is embedded in an elastic matrix, are treated in this study, and stress intensity factors are obtained for a range of crack sizes. Both mechanical and thermal loading are considered, and the effects of the mismatch of the material properties of the inclusion and the matrix are also investigated. The numerical results are of interest to the study of, for example, microcracking in particulate composites.