A phase-field model for the brittle fracture of Euler–Bernoulli beams coupling stretching and bending

Damage gradient models seek to simulate fracture mechanics by modulating the material stiffness. Within this framework, a singular scalar field representing damage is commonly utilized to globally decrease elastic energy. However, when considering structural models like beams and plates, this approach often fails to adequately capture important aspects due to the interaction between stretching and bending contributions. In this work, we propose a model for planar Euler–Bernoulli beams that incorporates the following features: firstly, we utilize two phase-field damage parameters to describe material damaging, specifically addressing the ‘erosion’ occurring above and below the original, undamaged beam; secondly, we assume a simple linear dependence of the axial stress field on the thickness coordinate, along with linear dependence on the axial force and bending moment. By appropriately identifying the constitutive response, our model effectively considers the coupling between stretching and bending induced during through-the-thickness damage, demonstrating good agreement with 2D observations.