Numerical modeling of the viscoplastic damage behaviour of rocks and application to underground storage facilities

The long-term behavior of large, underground works of a civil engineering nature carried out in a rock mass is currently the subject of numerous studies. The object is to attain a better understanding of complex phenomena, such as the convergence of excavated cavities or the outbreak and development of damaged zones in the rock mass neighboring the works, in order to foresee them. This Ph.D. thesis is devoted to the analysis of viscoplastic strain in rocks and to the degradation of their mechanical properties with time, often referred to as deferred damage. A bibliographical record presents the current depth of understanding as regards underlying microstructural phenomena and summarizes the main theories upon which the modeling of these phenomena at the macroscopic scale is based. The formulations enabling a coupling between the viscous effects and the deferred damage are revisited and discussed in detail. One phenomenological model in particular, Lemaitre's viscoplastic constitutive damage law is retained for the numerical modeling. The calculations were performed with the help of a finite element code (CAST3M). Designs of nuclear waste disposal structures at great depth make up the subject of different case studies. The Lemaitre model, originally designed for metallic materials, is next the subject of a theoretical development of which the aim is to better adapt it to the description of the long-term mechanical behavior of rocks. The modifications focus on several points; notably that the hypotheses of anelastic strain at constant volume and of isotropy of damage are rejected. The main characteristics of time-dependent strain in rocks; in particular the phenomena of viscoplastic dilation and contraction as well as the anisotropy induced by damage to the rock matrix are reproduced by the proposed model. A parametric study is then undertaken, using the experimental results obtained on different types of rock, in order to demonstrate the model's capabilities. Finally, the model is applied to the simple case of an axisymmetric cavity. (author)