Dynamic crack propagation with cohesive elements: a methodology to address mesh dependency

In this paper, two brittle fracture problems are numerically simulated: the failure of a ceramic ring under centrifugal loading and crack branching in a PMMA strip. A three-dimensional finite element package in which cohesive elements are dynamically inserted has been developed. The cohesive elements' strength is chosen to follow a modified weakest link Weibull distribution. The probability of introducing a weak cohesive element is set to increase with the cohesive element size. This reflects the physically based effect according to which larger elements are more likely to contain defects. The calculations illustrate how the area dependence of the Weibull model can be used to effectively address mesh dependency. On the other hand, regular Weibull distributions have failed to reduce mesh dependency for the examples shown in this paper. The ceramic ring calculations revealed that two distinct phenomena appear depending on the magnitude of the Weibull modulus. For low Weibull modulus, the fragmentation of the ring is dominated by heterogeneities. Whereas many cracks were generated, few of them could propagate to the outer surface. Monte Carlo simulations revealed that for highly heterogeneous rings, the number of small fragments was large and that few large fragments were generated. For high Weibull modulus, signifying that the ring is close to being homogeneous, the fragmentation process was very different. Monte Carlo simulations highlighted that a larger number of large fragments are generated due to crack branching. Copyright (C) 2003 John Wiley Sons, Ltd.