An immersed peridynamics method for fluid-driven damage and failure of anisotropic materials

TL;DR

This paper extends the immersed peridynamics method to simulate fluid-driven damage and failure in anisotropic materials with complex geometries, incorporating hyperelastic and ductile failure models for realistic biomaterial behavior.

Contribution

It introduces a novel extension of the immersed peridynamics method that handles anisotropic materials, complex geometries, and realistic failure processes within a fluid-structure interaction framework.

Findings

Non-uniform discretizations improve accuracy and volume conservation.

The method accurately captures damage growth, crack propagation, and rupture.

Simulations show grid convergence and realistic failure behaviors in anisotropic materials.

Abstract

The immersed peridynamics (IPD) method is a fluid-structure interaction (FSI) model to simulate fluid-driven material damage and failure of an immersed structure, in which a peridynamic (PD) constitutive correspondence model is employed within a classical immersed boundary (IB)-type framework to describe stresses, forces, and structural deformations of a structural body, instead of classical continuum mechanics. This paper introduces an extension of the IPD method to simulate fluid-driven structural deformation, damage, and failure of anisotropic materials with complex geometries. We use quadrature rules attached to finite element (FE) meshes to generate both the PD points and their associated weights, which are used to approximate the PD integrals. We demonstrate that non-uniform discretizations improve both accuracy and volume conservation of hyperelastic materials along with accurately represented boundaries. To capture realistic biomaterial behaviors, we incorporate hyperelastic constitutive models including both isotropy and anisotropy into the proposed IPD method. In addition, a ductile failure model is adopted to simulate realistic failure processes of anisotropic materials. For non-failure cases, our numerical simulations demonstrate that the extended IPD method yields comparable accuracy with similar numbers of structural degrees of freedom for different choices of peridynamic horizon sizes. For failure tests, we investigate the effect of a fiber orientation on deformations and failure processes using realistic biomaterial models with varying fiber directions. We further demonstrate that the developed method generates grid-converged simulations of damage growth, crack formation and propagation, and rupture under large deformations, including purely fluid-driven failure processes.