One of the reasons why the phase field approach has become very popular for modeling fracture processes is that therein the entire evolution of fracture follows from energetic principles without the need for different criteria for crack initiation, nucleation, kinking or branching. Furthermore, these models allow for a straightforward numerical implementation with standard finite elements, since displacement jumps and stress singularities are avoided in these models. Consequently, phase field fracture models are now also used to model and simulate fracture of heterogeneous materials such as composites, where cracks may be arrested, deflected or bifurcated at interfaces between the different components of the composite material, resulting in complicated crack patterns. The crack propagation is not only controlled by the elastic and fracture properties of the bulk phases of the material but also by the properties of the interfaces in between. Thus, it is crucial to take these interface properties into account. Therefore in this work, interfaces are equipped with a fracture energy corresponding to the phase field fracture energy of the bulk phases in order to model the interface fracture properties. In the finite element implementation, the interfaces are modeled by means of surface elements without introducing additional degrees of freedom.
Phase field modeling of interfacial fractureThe overall fracture properties of heterogeneous materials not only depend on the elastic and fracture mechanical properties of the individual constituents α, i.e. the elastic stiffness tensor C (α) and the fracture resistance G (α) c , but also on the properties of the interfaces in between the different materials. While different properties of individual bulk phases can easily be dealt with in phase field fracture models, there are several approaches on how to integrate on how to incorporate the effect of interfaces. In [1] interfaces are modeled as zones of finite width with appropriately adjusted fracture parameters. Contrary, a combination with a cohesive type of modeling interface fracture is proposed in [2]. Another approach, which makes use of a new energetic formulation mixing bulk and cohesive surface energies, is proposed in [3].In this work, the bulk phase field energy accounting for elastic and fracture energy contributions of the different phases is supplemented by an additional interface term accounting for the effect of the interfaces, such that the total energy of a volume V = V (α) with interfaces Γ = Γ (αβ) readsTherein ε and s denote the infinitesimal strain tensor and the fracture field, respectively. The size of the transition zone between sound (s = 1) and broken (s = 0) material is controlled by the regularization parameter . The small stiffness margin η ensures well posedness when s = 0. By means of the interface parameter G (αβ) Γ , the fracture resistance at an interface between phases α and β can be adjusted to the effective valueNote that neglecting the interface contribution G (αβ) Γ yields the average of ...