Coupling a discrete twin model with cohesive elements to understand twin-induced fracture

The interplay between twinning and fracture in metals under deformation is an open question. The plastic strain concentration created by twin bands can induce large stresses on the grain boundaries. We present simulations in which a continuum model describing discrete twins is coupled with a crystal plasticity finite element model and a cohesive zone model for intergranular fracture. The discrete twin model can predict twin nucleation, propagation, growth and the correct twin thickness. Therefore, the plastic strain concentration in the twin band can be modelled. The cohesive zone model is based on a bilinear traction-separation law in which the damage is caused by the normal stress on the grain boundary. An algorithm is developed to generate interface elements at the grain boundaries that satisfy the traction-separation law. The model is calibrated by comparing polycrystal simulations with the experimentally observed strain to failure and maximum stress. The dynamics of twin and crack nucleation have been investigated. First, twins nucleate and propagate in a grain, then, microcracks form near the intersection between twin tips and grain boundaries. Microcracks appear at multiple locations before merging. A propagating crack can nucleate additional twins starting from the grain boundary, a few micrometres away from the original crack nucleation site. This model can be used to understand which type of texture is more resistant against crack nucleation and propagation in cast metals in which twinning is a deformation mechanism. The code is available online at https://github.com/TarletonGroup/CrystalPlasticity