the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
A dilatant visco-elasto-viscoplasticity model with globally continuous tensile cap: stable two-field mixed formulation
Abstract. Rocks break if shear stresses exceed their strength. It is therefore important for typical geoscientific applications to take shear failure mechanism and the subsequent development of mode-II shear bands or faults into account. Many existing codes incorporate non-associated Drucker-Prager or Mohr-Coulomb plasticity models to simulate this behavior. Yet, when effective mean stress becomes extensional, for example when fluid pressure becomes large, the dominant failure mode changes to a mode-I (opening) mode, which initiates plastic volumetric deformation. It is rather difficult to represent both failure modes in numerical models in a self-consistent manner, while also accounting for the nonlinear visco-elastic host rock rheology, which varies from being nearly incompressible in the mantle to being compressible in surface-near regions. Here, we present a simple plasticity model that is designed to overcome these difficulties. We employ a combination of a linearized Drucker-Prager shear failure envelope with a circular tensile cap function in way that ensures continuity and smoothness of both yield surface and flow potential in the entire stress space. A Perzyna-type viscoplastic regularization ensures that the resulting localization zones are mesh-insensitive. To deal with the near incompressibility condition, a mixed two-field finite element formulation is employed. The local nonlinear iterations at the integration-point level are used to determine the stress increments. The global Newton-Raphson iterations are applied to solve the discretized momentum and continuity residual equations. The presented plasticity model is implemented in an open-source 2D unstructured finite element code GeoTech2D. The results of several typical test cases that range from crustal scale deformation to the propagation of fluid-induced tensile failure zones demonstrate rapid convergence. The robustness of the solution scheme is enhanced by the adaptive time stepping algorithm.
Competing interests: At least one of the (co-)authors is a member of the editorial board of Geoscientific Model Development.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.- Preprint
(5589 KB) - Metadata XML
- BibTeX
- EndNote
Status: closed
-
RC1: 'Comment on egusphere-2025-2469', Anonymous Referee #1, 26 Jun 2025
This is a very solid and interesting manuscript, which addresses the proper modelling of the transition between brittle and ductile failure. Publication is recommended, but the authors may wish to pay attention to the following issues in a revised version:
* While viscosity (mostly) remedies mesh dependence, it should be pointed out that it is a weak regularisation technique in quasi-static loading conditions, as considered here. Only under dynamic loading conditions full regularisation of the governing equations can be proven, see De Borst and Duretz (2020). This should be discussed briefly in the paper in order to avoid misconceptions in the community.
* Work on combining mode-I and mode-II plasticity has been pursued before in the context of modelling of concrete, albeit for plane stress rather than for plane strain or 3D conditions as is the focus of the current paper, see e.g. Feenstra and de Borst. Int. J. Solids Structures (1996) 33, 707 - 730. It would be advisable to put the present contribution also in that context.
* The discussion on volumetric locking is confusing and perhaps misleading. It is suggested that volumetric locking occurs for isochoric plastic deformations and it is not (explicitly) pointed out that this phenomenon also occurs for dilatant or contractant plastic flow. Indeed, in all cases a kinematic constraint is imposed at failure, i.e. when the elastic deformations vanish, see De Borst and Groen, Int. J. Num. Meth. Engng (1995) 38, 2887 - 2906.Â
Citation: https://doi.org/10.5194/egusphere-2025-2469-RC1 - AC1: 'Comment on egusphere-2025-2469', Anton Popov, 15 Aug 2025
-
RC2: 'Comment on egusphere-2025-2469', Anonymous Referee #2, 14 Jul 2025
This is a robust manuscript looking at the addition of a tensile cap for modelling rock failure in a geodynamic context. The manuscript is of excellent quality and should eventually be published to benefit the community. There is however a point where the manuscript could be improved, regarding viscoplasticity and strain regularization.
- The manuscript presents Perzyna's viscoplasticity as a way to regularize the ill-posed problem of strain localization. While this is valid way to consider viscoplasticity, there is no mention in the manuscript that viscoplasticity could also capture a physical behaviour of rate-dependent plastic deformation (which is not uncommon).
- The authors claim that introducing viscoplasticity in a quasi-static approximation is sufficient to regularize the strain localization problem. I am not aware of any studies proving this statement. A number of studies analyzed the role of viscoplasticity for strain regularization and, to my knowledge, they all included the dynamic terms.
- Furthermore, the authors should consider discussing other studies about viscoplasticity and strain regularization, some of which are more critical about the role of viscoplasticity than the cited studies in the manuscript (Jacquey et al. (2021) and Stathas and Stefanou (2022) as examples).
References :
Stathas, A. and Stefanou, I. The role of viscous regularization in dynamical problems, strain localization and mesh dependency. Computer Methods in applied Mechanics and Engineering, 388, 113185, 2022.
Jacquey, A. B., Rattez, H. and Veveakis, M. Strain localization regularization and patterns formation in rate-dependent plastic materials with multiphysics coupling. Journal of the Mechanics and Physics of Solids, 152, 104422, 2021.Citation: https://doi.org/10.5194/egusphere-2025-2469-RC2 - AC1: 'Comment on egusphere-2025-2469', Anton Popov, 15 Aug 2025
- AC1: 'Comment on egusphere-2025-2469', Anton Popov, 15 Aug 2025
Status: closed
-
RC1: 'Comment on egusphere-2025-2469', Anonymous Referee #1, 26 Jun 2025
This is a very solid and interesting manuscript, which addresses the proper modelling of the transition between brittle and ductile failure. Publication is recommended, but the authors may wish to pay attention to the following issues in a revised version:
* While viscosity (mostly) remedies mesh dependence, it should be pointed out that it is a weak regularisation technique in quasi-static loading conditions, as considered here. Only under dynamic loading conditions full regularisation of the governing equations can be proven, see De Borst and Duretz (2020). This should be discussed briefly in the paper in order to avoid misconceptions in the community.
* Work on combining mode-I and mode-II plasticity has been pursued before in the context of modelling of concrete, albeit for plane stress rather than for plane strain or 3D conditions as is the focus of the current paper, see e.g. Feenstra and de Borst. Int. J. Solids Structures (1996) 33, 707 - 730. It would be advisable to put the present contribution also in that context.
* The discussion on volumetric locking is confusing and perhaps misleading. It is suggested that volumetric locking occurs for isochoric plastic deformations and it is not (explicitly) pointed out that this phenomenon also occurs for dilatant or contractant plastic flow. Indeed, in all cases a kinematic constraint is imposed at failure, i.e. when the elastic deformations vanish, see De Borst and Groen, Int. J. Num. Meth. Engng (1995) 38, 2887 - 2906.Â
Citation: https://doi.org/10.5194/egusphere-2025-2469-RC1 - AC1: 'Comment on egusphere-2025-2469', Anton Popov, 15 Aug 2025
-
RC2: 'Comment on egusphere-2025-2469', Anonymous Referee #2, 14 Jul 2025
This is a robust manuscript looking at the addition of a tensile cap for modelling rock failure in a geodynamic context. The manuscript is of excellent quality and should eventually be published to benefit the community. There is however a point where the manuscript could be improved, regarding viscoplasticity and strain regularization.
- The manuscript presents Perzyna's viscoplasticity as a way to regularize the ill-posed problem of strain localization. While this is valid way to consider viscoplasticity, there is no mention in the manuscript that viscoplasticity could also capture a physical behaviour of rate-dependent plastic deformation (which is not uncommon).
- The authors claim that introducing viscoplasticity in a quasi-static approximation is sufficient to regularize the strain localization problem. I am not aware of any studies proving this statement. A number of studies analyzed the role of viscoplasticity for strain regularization and, to my knowledge, they all included the dynamic terms.
- Furthermore, the authors should consider discussing other studies about viscoplasticity and strain regularization, some of which are more critical about the role of viscoplasticity than the cited studies in the manuscript (Jacquey et al. (2021) and Stathas and Stefanou (2022) as examples).
References :
Stathas, A. and Stefanou, I. The role of viscous regularization in dynamical problems, strain localization and mesh dependency. Computer Methods in applied Mechanics and Engineering, 388, 113185, 2022.
Jacquey, A. B., Rattez, H. and Veveakis, M. Strain localization regularization and patterns formation in rate-dependent plastic materials with multiphysics coupling. Journal of the Mechanics and Physics of Solids, 152, 104422, 2021.Citation: https://doi.org/10.5194/egusphere-2025-2469-RC2 - AC1: 'Comment on egusphere-2025-2469', Anton Popov, 15 Aug 2025
- AC1: 'Comment on egusphere-2025-2469', Anton Popov, 15 Aug 2025
Model code and software
Unstructured FEM code GeoTech2D Anton A. Popov and Boris J. P. Kaus https://doi.org/10.5281/zenodo.15496842
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
545 | 79 | 12 | 636 | 17 | 27 |
- HTML: 545
- PDF: 79
- XML: 12
- Total: 636
- BibTeX: 17
- EndNote: 27
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1