A finite deformation phase field model for electromechanical fracture of flexible piezoelectric materials

Fracture failure is a major concern in mechanical engineering, particularly for piezoelectric materials. In contrast to other numerical methods, phase field method has significant advantages in addressing fracture progress. It can automatically track crack surfaces through ordered parameter evolution, which is versatile for modeling complex fracture behaviors. However, previous studies on phase field fracture in piezoelectric solids have primarily focused on brittle ceramics with small deformation. In recent years, flexible piezoelectrics with high stretchability have already been achieved in industrial production. These materials exhibit obvious nonlinear characteristics during deformation, which renders the traditional assumption of small deformations inadequate for predicting their fracture behaviors. In this work, we propose a finite deformation phase field fracture model for flexible piezoelectric materials, building upon the established nonlinear electromechanical material model. The numerical framework is carried out in the commercial software ABAQUS via a user element subroutine. Both single–pass staggered algorithm (SPSA) and residual-controlled staggered algorithm with even-odd iteration split (RCSA-EO) are employed to solve coupled electro-mechanical-phase field governing equations. The proposed model is validated through comparisons with analytical solutions and existing literature. Moreover, the developed numerical framework effectively explains the nonlinear fracture behavior observed in experiments conducted on Polyvinylidene fluoride (PVDF), a flexible piezoelectric material with a large failure strain. Numerical simulations are also performed to demonstrate the influence of the applied electric field on electromechanical fracture behavior. The results highlight that the specific impact of electric fields depends on material parameters, geometric parameters, and boundary conditions. The developed model is capable of making accurate and realistic predictions of fracture in flexible piezoelectric materials. This is particularly important for evaluating the reliability and safety of flexible piezoelectric devices.