Cohesive crack propagation in a random elastic medium

The issue of generating non-Gaussian, multivariate and correlated random fields, while preserving the internal auto-correlation structure of each single-parameter field, is discussed with reference to the problem of cohesive crack propagation. Three different fields are introduced to model the spatial variability of the Young modulus, the tensile strength of the material, and the fracture energy, respectively. Within a finite-element context, the crack-propagation phenomenon is analyzed by coupling a Monte Carlo simulation scheme with an iterative solution algorithm based on a truly-mixed variational formulation which is derived from the Hellinger–Reissner principle. The selected approach presents the advantage of exploiting the finite-element technology without the need to introduce additional modes to model the displacement discontinuity along the crack boundaries. Furthermore, the accuracy of the stress estimate pursued by the truly-mixed approach is highly desirable, the direction of crack propagation being determined on the basis of the principal-stress criterion. The numerical example of a plain concrete beam with initial crack under a three-point bending test is considered. The statistics of the response is analyzed in terms of peak load and load–mid-deflection curves, in order to investigate the effects of the uncertainties on both the carrying capacity and the post-peak behaviour. A sensitivity analysis is preliminarily performed and its results emphasize the negative effects of not accounting for the auto-correlation structure of each random field. A probabilistic method is then applied to enforce the auto-correlation without significantly altering the target marginal distributions. The novelty of the proposed approach with respect to other methods found in the literature consists of not requiring the a priori knowledge of the global correlation structure of the multivariate random field.