ABSTRACT
The elastic-brittle-damage constitutive model for three-dimensional stress state is incorporated into the commercial finite element code ABAQUS. To take into account the heterogeneity effect, mesoscale elements in finite element analysis are assigned by different strength and stiffness properties according to the Weibull distribution function. A FORTRAN code is developed and employed within the static equilibrium equation to analyse three-dimensional progressive failure of rock masses. The numerical model is extended to the coupled hydro-mechanical problem by adding the fluid mass conservation equation to the governing equation system and relating hydraulic conductivity of the elements to their damage and hydrostatic stress states. The results for two typical problems with the pure dry and coupled solid/fluid conditions, respectively, are presented. Good agreement between the numerical results and previously findings in the literature verify such implementation and demonstrates its significance.
