Dynamic analysis of a viscoelastic orthotropic cracked body using the extended finite element method

•XFEM simulations were conducted on a pulse loaded cracked body.•The material is assumed orthotropic and strain rate sensitive (viscoelastic).•SIF’s and J-integral were calculated using the in-house XFEM code.•A Material User Subroutine UMAT has been developed in ABAUQS.•Strong correlation is observed between XFEM and FEM. The extended finite element method (XFEM) is found promising in approximating solutions to locally non-smooth features such as jumps, kinks, high gradients, inclusions, voids, shocks, boundary layers or cracks in solid or fluid mechanics problems. The XFEM uses the properties of the partition of unity finite element method (PUFEM) to represent the discontinuities without the corresponding finite element mesh requirements. In the present study numerical simulations of a dynamically loaded orthotropic viscoelastic cracked body are performed using XFEM and the J-integral and stress intensity factors (SIF’s) are calculated. This is achieved by fully (reproducing elements) or partially (blending elements) enriching the elements in the vicinity of the crack tip or body. The enrichment type is restricted to extrinsic mesh-based topological local enrichment in the current work. Thus two types of enrichment functions are adopted viz. the Heaviside step function replicating a jump across the crack and the asymptotic crack tip function particular to the element containing the crack tip or its immediately adjacent ones. A constitutive model for strain-rate dependent moduli and Poisson ratios (viscoelasticity) is formulated. A symmetric double cantilever beam (DCB) of a generic orthotropic material (mixed mode fracture) is studied using the developed XFEM code. The same problem is studied using the viscoelastic constitutive material model implemented in ABAQUS through an implicit user defined material subroutine (UMAT). The results from XFEM correlate well with those of the finite element method (FEM). Three cases viz. static, dynamic and viscoelastic dynamic are studied. It is shown that there is an increase in the value of maximum J-integral when the material exhibits strain rate sensitivity.