| Abstract: | This paper implements a domain integral energy method for modelling crack growth in composite material shell structures using the finite element method. Volume integral expressions to evaluate the dynamic energy release rate in a through‐thickness three‐dimensional crack are derived. Using the domain integral, the energy release rate computation is implemented in the DYNA3D explicit non‐linear dynamic finite element analysis program wherein crack propagation is modelled by releasing the constraints between initially constrained node pairs. The implementation enables the program to either determine the energy resistance response for the material (provided experimental data is available) or predict the rate of crack propagation in shell structures. The numerical implementation was verified by simulating mode I and mode III slow crack growth problems in semi‐infinite transversely isotropic media, for which analytic solutions are available. Oscillations of energy following the release of nodal constraints as the crack propagates in discrete increments were suppressed using light mass proportional damping and a moving averaging scheme. Copyright © 2004 John Wiley & Sons, Ltd. |
