Towards optimization of crack resistance of composite materials by adjustment of fiber shapes

We investigate the evolution and propagation of cracks in 2-d elastic domains, which are subjected to quasi-static loading scenarios. In addition to the classical variational formulation, where the standard potential energy is minimized over the cracked domain under physical conditions characterizing the behavior of the material close to the crack (e.g. non-penetration conditions), we include a ‘cohesive traction term’ in the energy expression. In this way we obtain a mathematically concise set of partial differential equations with non-linear boundary conditions at the crack interfaces. We perform a finite element discretization using a combination of standard continuous finite elements and so-called cohesive elements. During the simulation process cohesive elements are adaptively inserted at positions where a certain stress bound is exceeded. In our numerical studies we consider domains consisting of a matrix material with fiber inclusions. Beyond pure crack path simulation, our ultimate goal is to determine an optimal shape of the fibers resulting in a crack path that releases for a given load scenario as much energy as possible without destroying the specimen completely. We develop a corresponding optimization model and propose a solution algorithm for the same. The article is concluded by numerical results.