Elasto-thermo-visco-plastic numerical modelling from a laboratory to geodynamic scale: implications for convergence-driven experiments (2212.12701v1)

Abstract: The development of a subduction zone, whether spontaneous or induced, encompasses a stage of strain localization and is epitomized by the growth of lithospheric-scale shear bands. Our aim in this paper, using a solid-mechanical constitutive description relevant for oceanic lithosphere, is to investigate factors that promote or inhibit localization of deformation in brittle and ductile regimes in convergence-driven numerical experiments. We used the Drucker-Prager yield criterion and a non-associative flow rule, allowing viscoplastic deformation to take directions independent of the preferred direction of yield. We present a step-by-step description of the constitutive law and the consistent algorithmic tangent modulus. The model domain contains an initial weak-zone on which localization can potentially nucleate. In solving the energy conservation problem, we incorporate a heat source term from the mechanical deformations which embodies the irreversible plastic work done. This work term couples the energy equation to the constitutive description, and hence hence the stress balance, via the evolving temperature field. On a sample-scale, we first conduct a series of isothermal benchmark tests. We then explore behavior including shear heating and volumetric work both separately and in concert. and thereby address the (in)significance of the latter, and hence assess their potential importance. We find that dilatational effects mostly enhance both shear band development and shear heating. We also observe that high temperature promotes shear band development whereas high confining pressure inhibits it, and infer that the competition between these factors is likely to be the major influence on the position within the lithosphere where shear bands nucleate.