To cite this version:Ioannis Stefanou, Hadrien Rattez, Jean Sulem. The importance of Thermo-Hydro-Mechanical couplings and microstructure to strain localization in 3D continua with application to seismic faults. Part I: Theory and Linear Stability Analysis.Abstract A Thermo-Hydro-Mechanical (THM) model for Cosserat continua is developed to explore the influence of frictional heating and thermal pore fluid pressurization on the strain localization phenomenon. A general framework is presented to conduct a bifurcation analysis for elasto-plastic Cosserat continua with THM couplings and predict the onset of instability. Furthermore, the presence of an internal length in Cosserat continua enables to estimate the thickness of the localization zone. This is done by performing a linear stability analysis of the system and looking for the selected wavelength corresponding to the instability mode with fastest finite growth coefficient. These concepts are applied to the study of fault zones under fast shearing. For doing so, we consider a model of a sheared saturated infinite granular layer. The influence of THM couplings on the bifurcation state and the shear band width is investigated. Taking representative parameters for a centroidal fault gouge, the evolution of the thickness of the localized zone under continuous shear is studied. Furthermore, the effect of grain crushing inside the shear band is explored by varying the internal length of the constitutive law. 45 the fact that it has been recognized to play an important role on fault gouge behavior (Anthony and Marone, 2005, Cashman et al., 2007, Phillips and White, 2017).Cosserat continuum (Cosserat and Cosserat, 1909) is a special case of micromorphic continua (Germain, 1973, Godio et al., 2016, also called generalized 50 or higher order continua. In addition to the translational degrees of freedom of the Cauchy continuum, Cosserat theory considers rotational degrees of freedom at the material point that allow for a better representation of the physics and the mechanical behavior of heterogeneous solids with non-negligible microstructure. Cosserat continuum theory naturally incorporates one or several 55 material lengths related to the microstructure in the constitutive equations of the material (see Appendix B).Cosserat continuum has been previously used for studying the behavior of fault gouges and strain localization , Veveakis et al., 2013. In 3 these works the conditions for the onset of localization were investigated under 60 THM couplings with a microstructure of given size. In the present papers (Part I and II) we extend the aforementioned works by studying (a) the evolution of the localization zone thickness and its dependency on various parameters such the size of the microstructure, (b) the full stress-strain response of the fault gouge, which is related to the transition from seismic to aseismic slip (Scholz, 65 2002, Tse andRice, 1986) and (c) the apparent rate dependency of the system due to THM couplings even under rate-independent constitutive law...