“…Of course, more number of replications will include more possible events (provides more informative data) that result in a more accurate expected value and variance. For the quasi-brittle materials, the Monte Carlo finite element method (MC-FEM) was used to model the dependence of the computed crack probabilities on the type of perturbation in [55,56,57], and the polynomial chaos expansion in functionally graded materials with random material properties is used to model the phase-field fracture, see [58]. In computational mechanics, stochastic discretization techniques have been employed for variational theory for nonlinear problems with stochastic coefficients [59,60,61], inelastic media under uncertainty [62], elastic-plastic material with uncertain parameters [63], fatigue crack propagation due to the inherent uncertainties according to the material properties [64], nonlinear fracture mechanics of concrete [65], and stochastic fracture response and crack growth analysis of laminated composites [66].…”