This paper describes a fully implicit return mapping algorithm for the numerical integration of a three-invariant non-associative plasticity model for hard rocks. Bifurcation analysis is performed to detect the onset of localization during the deformation process. A nonlinear Matsuoka-Nakai yield criterion along with an isotropic cohesive-softening law is employed to capture the inelastic material response. The integration algorithm guarantees a numerically stable and first-order accurate solution at the Gauss point level. The formulation of the model includes a consistent linearization of the weak form of the linear momentum balance to enable optimal convergence for Newton-Raphson iterations. The onset of strain localization is detected us-ing classical bifurcation theory by checking for the singularity of the elastoplastic constitutive operator at each step of the calculation process. The accuracy and convergence properties of the proposed algorithm are assessed through numerical simulations of single element tests and through a plane strain boundary value problem.

A three-invariant non-associative plasticity model for hard rocks: implicit integration and strain localization analysis

BOLDINI, DANIELA;
2007

Abstract

This paper describes a fully implicit return mapping algorithm for the numerical integration of a three-invariant non-associative plasticity model for hard rocks. Bifurcation analysis is performed to detect the onset of localization during the deformation process. A nonlinear Matsuoka-Nakai yield criterion along with an isotropic cohesive-softening law is employed to capture the inelastic material response. The integration algorithm guarantees a numerically stable and first-order accurate solution at the Gauss point level. The formulation of the model includes a consistent linearization of the weak form of the linear momentum balance to enable optimal convergence for Newton-Raphson iterations. The onset of strain localization is detected us-ing classical bifurcation theory by checking for the singularity of the elastoplastic constitutive operator at each step of the calculation process. The accuracy and convergence properties of the proposed algorithm are assessed through numerical simulations of single element tests and through a plane strain boundary value problem.
Numerical Models in Geomechanics - NUMOG X
205
210
D. Boldini; R.I. Borja; P.F. Sanz
File in questo prodotto:
Eventuali allegati, non sono esposti

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: http://hdl.handle.net/11585/45748
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 0
  • ???jsp.display-item.citation.isi??? 0
social impact