An investigation of numerical method for three-dimensional elastic-plastic flows

TY - JOUR

T1 - An investigation of numerical method for three-dimensional elastic-plastic flows

AU - Niu, Panyu

AU - Wang, Cheng

AU - Jia, Zupeng

N1 - Publisher Copyright: © 2026 Elsevier B.V. All rights are reserved, including those for text and data mining, AI training, and similar technologies.

PY - 2026/6/15

Y1 - 2026/6/15

N2 - A cell-centered Lagrangian discontinuous Galerkin (DG) method is developed for simulating elastic-plastic flows on three-dimensional unstructured hexahedral meshes. The elastic-plastic deformation is described by the hypo-elastic model with Jaumann stress rate, while plasticity is computed using the radial return algorithm. The elastic limit state is governed by the von Mises yield condition, with yield strength determined by either perfectly plastic model, linear hardening model, or Johnson-Cook model. A new Riemann solver for three-dimensional elastic-plastic problems is constructed, in which the Cauchy stress tensor is calculated at the centers of the triangular facets of the hexahedral elements and then used to compute the numerical fluxes of the semi-discrete equations for the conserved variables. Furthermore, a novel nodal solver is developed for calculating velocities of both interior nodes and boundary nodes subject to various boundary conditions. This nodal solver operates on all sub-triangular facets surrounding each node, rigorously enforcing global conservation of both total energy and momentum while satisfying entropy inequality constraints. The DG coefficients of the conserved variables and the deviatoric stress tensor are updated by using a second-order total variation diminishing Runge-Kutta time integration scheme. A least squares reconstruction algorithm with Barth-Jesperson limiter is used to eliminate the non-physical numerical oscillation of the conserved variables and the deviatoric stress tensor. Due to the high computational cost of three-dimensional simulations, a parallel algorithm on unstructured hexahedral meshes is implemented to improve efficiency. Several numerical tests and a projectile impact experiment are conducted to validate the proposed scheme. Experimental validations and numerical results confirm the accuracy and robustness of the scheme developed in this work.

AB - A cell-centered Lagrangian discontinuous Galerkin (DG) method is developed for simulating elastic-plastic flows on three-dimensional unstructured hexahedral meshes. The elastic-plastic deformation is described by the hypo-elastic model with Jaumann stress rate, while plasticity is computed using the radial return algorithm. The elastic limit state is governed by the von Mises yield condition, with yield strength determined by either perfectly plastic model, linear hardening model, or Johnson-Cook model. A new Riemann solver for three-dimensional elastic-plastic problems is constructed, in which the Cauchy stress tensor is calculated at the centers of the triangular facets of the hexahedral elements and then used to compute the numerical fluxes of the semi-discrete equations for the conserved variables. Furthermore, a novel nodal solver is developed for calculating velocities of both interior nodes and boundary nodes subject to various boundary conditions. This nodal solver operates on all sub-triangular facets surrounding each node, rigorously enforcing global conservation of both total energy and momentum while satisfying entropy inequality constraints. The DG coefficients of the conserved variables and the deviatoric stress tensor are updated by using a second-order total variation diminishing Runge-Kutta time integration scheme. A least squares reconstruction algorithm with Barth-Jesperson limiter is used to eliminate the non-physical numerical oscillation of the conserved variables and the deviatoric stress tensor. Due to the high computational cost of three-dimensional simulations, a parallel algorithm on unstructured hexahedral meshes is implemented to improve efficiency. Several numerical tests and a projectile impact experiment are conducted to validate the proposed scheme. Experimental validations and numerical results confirm the accuracy and robustness of the scheme developed in this work.

KW - Cell-centered Lagrangian scheme

KW - Discontinuous Galerkin method

KW - Nodal solver

KW - Riemann solver

KW - Three-dimensional elastic-plastic flows

UR - https://www.scopus.com/pages/publications/105035605494

U2 - 10.1016/j.cma.2026.118879

DO - 10.1016/j.cma.2026.118879

M3 - Article

AN - SCOPUS:105035605494

SN - 0045-7825

VL - 455

JO - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

M1 - 118879

ER -