Crack propagation during the indentation test of a ferroelectric single crystal is simulated using a phase- eld model. This model is based on variational formulations of brittle crack propagation and domain evolution in ferroelectric materials. Due to the high compressive stresses near the indenter contact faces, a modi ed regularized formulation of the variational brittle fracture is coupled with the material model to prevent crack formation and interpenetration in the compressed regions. The simulation results show that the radial cracks perpendicular to the poling direction of the material propagate faster than the parallel ones, which is in agreement with experimental observations. This anisotropy in the crack propagation is due to interactions between the material microstructure and the radial cracks, as captured by the phase- eld simulation. Crack propagation during the indentation test of a ferroelectric single crystal is simulated using a phase-eld model. This model is based on variational formulations of brittle crack propagation and domain evolution in ferroelectric materials. Due to the high compressive stresses near the indenter contact faces, a modied regularized formulation of the variational brittle fracture is coupled with the material model to prevent crack formation and interpenetration in the compressed regions. The simulation results show that the radial cracks perpendicular to the poling direction of the material propagate faster than the parallel ones, which is in agreement with experimental observations. This anisotropy in the crack propagation is due to interactions between the material microstructure and the radial cracks, as captured by the phase-eld simulation
