The phase field fracture model is attracting significant interest. To model fracture in hyperelastic material and hydrogel, we have implemented a robust two- and three-dimensional phase field method in the commercial finite element code ABAQUS/Standard. The method is based on the rate-independent variational principle of diffuse fracture and also exploits the analogy between the phase field evolution law and the heat transfer equation, enabling the use of Abaqus’ in-built features and sparing the need for defining user elements. The framework is shown to accommodate both staggered and monolithic solution schemes. The approach can properly simulate the fracture for both hyperelastic material and hydrogel under different boundary conditions. Several examples are provided to demonstrate the robustness of the method. The provided source codes and the tutorials make it easy for practicing engineers and scientists to model crack propagation in hyperelastic and gel materials.
Polymer network is a crucial component of hydrogels, and network imperfection is a prominent feature of polymer networks, significantly influencing the performance of hydrogels. Two essential features of network imperfection are unequal chain lengths and dangling chains, both of which have a significant impact on the mechanical properties of single-network (SN) hydrogels. However, a theoretical framework considering network imperfection in SN hydrogels is still lacking. Here, we propose a theoretical model for SN hydrogels with network imperfection to study the damage behavior during deformation, in which we adopt different chain length distributions to accurately depict the real physical characteristics of the polymer network and incorporate the normalized critical chain force for a more precise measurement of network damage. To verify our theory, we discuss the effects of model parameters on the stress-stretch responses of SN hydrogels and predict the results of uniaxial loading-unloading tests of SN hydrogels, which agree well with experimentally measured stress-stretch behaviors. Finally, we implement the constitutive model into ABAQUS as a user subroutine to study the inhomogeneous deformation of hydrogels.