Energy release rate along a three-dimensional crack front in a thermally stressed body

Based on a line-integral expression for the energy release rate in terms of crack tip fields, which is valid for general material response, a (area/volume) domain integral expression for the energetic force in a thermally stressed body is derived. The general three-dimensional finite domain integral expression and the two-dimensional and axisymmetric specializations for the energy release rate are given. The domain expression is naturally compatible with the finite element formulation of the field equations. As such it is ideally suited for efficient and accurate calculation of the pointwise values of the energy release rate along a three-dimensional crack front. The finite element implementation of the domain integral corresponds to the virtual crack extension technique. Procedures for calculating the energy release rate using the numerically determined field solutions are discussed. For illustrative purposes several numerical examples are presented.