Math @ Duke

Publications [#268313] of John E. Dolbow
Papers Published
 Gosz, M; Dolbow, J; Moran, B, Domain integral formulation for stress intensity factor computation along curved threedimensional interface cracks,
International Journal of Solids and Structures, vol. 35 no. 15
(1998),
pp. 17631783 [S00207683(97)001327]
(last updated on 2018/03/16)
Abstract: New domain integrals for extracting mixedmode stress intensity factors along curved, threedimensional bimaterial interface cracks are derived. In the derivation, the asymptotic auxiliary fields for the plane problem of a bimaterial interface crack are imposed along a curved crack front. The general cracktip interaction integral (a contour integral surrounding a point on the crack front which is evaluated in the limit as the contour is shrunk onto the crack tip) is then expressed in domain form which is more suitable for numerical computation. A consequence of imposing the auxiliary fields along a curved crack front is that the auxiliary stress fields do not satisfy equilibrium, and the auxiliary strain fields do not satisfy compatibility. The terms which arise due to the lack of equilibrium and compatibility are not sufficiently singular to contribute to the cracktip interaction integral or to affect its path independence; however, these terms become important when domain integral representations are introduced, because they involve fields which are not asymptotically close to the crack tip. In order to compute the pointwise stress intensity factors along the crack front, the domain integrals are evaluated as a postprocessing step in the finite element method. In the numerical results, it is demonstrated that it is crucial to incorporate the terms that arise due to lack of equilibrium and compatibility of the auxiliary fields, especially in regions where the local crack front curvature is high. In the paper, we present two numerical examples. As a benchmark, we first consider the problem of a pennyshaped interface crack embedded in a cylinder. The results for the complex stress intensity factor and phase angle are found to be in excellent agreement with the analytical solution. The problem of an elliptical crack embedded between two dissimilar isotropic materials is also considered, and the results are discussed. © 1998 Elsevier Science Ltd.
Keywords: crackedge stress field analysis;delamination;fracture mechanics;interface phenomena;stressstrain relations;


dept@math.duke.edu
ph: 919.660.2800
fax: 919.660.2821
 
Mathematics Department
Duke University, Box 90320
Durham, NC 277080320

