Separation of crack extension modes in orthotropic delamination models

In modeling a crack along a distinct interface between dissimilar elastic materials, the ratio of mode I to mode II stress intensity factors or energy release rates is typically not unique, due to oscillatory behavior of near-tip stresses and displacements. Although methods have been developed for comparing mode mixes for isotropic interfacial fracture problems, this behavior currently limits the applicability of interfacial fracture mechanics in predicting delamination in layered materials without isotropic symmetry. The virtual crack closure technique (VCCT) is a method used to extract mode I and mode II energy release rate components from numerical fracture solutions. Energy release rate components extracted from an oscillatory solution using the VCCT are not unique due to their dependence on the virtual crack extension length, Delta. In this work, a method is presented for using the VCCT to extract Delta-independent energy release rate quantities for the case of an interface crack between two in-plane orthotropic materials. The method does not involve altering the analysis to eliminate its oscillatory behavior and it is similar to existing methods for extracting a mode mix from isotropic interfacial fracture models. Knowledge of near-tip fields is used to determine the explicit Delta dependence of energy release rate parameters. Energy release rates are then defined that are separated from the oscillatory dependence on Delta. A modified VCCT using these energy release rate definitions is applied to results from finite element analyses, showing that Delta-independent energy release rate quantities result. The modified technique has potential as a consistent method for extracting a mode mix from numerical solutions. The Delta-independent energy release rate quantities extracted using this technique can also aid numerical modelers, serving as guides for testing the convergence of finite element models. Direct applications of this work include the analysis of planar composite delamination problems, where plies or debonded laminates are modeled as in-plane orthotropic materials.