As extensions to our Lagrangian finite-strain plasticity framework based on the approximate exponential integrator, we introduce two new algorithms: (i) a fixed-radius trust region root finder for the nonlinear system involving the elastic right Cauchy Green tensor and plastic multiplier (ii) a partitioned approach for the hardening variables, which are determined in a staggered form using the strongly-coupled concept typically adopted for multiphysics problems. This allows the use of intricate hyperelastic laws combined with recent yield functions, which otherwise would involve a laborious treatment, and the use of corresponding work-hardening. Work-hardening would introduce significant nonlinearities in the constitutive system if used in a fully-coupled form. For the partitioned approach, a dynamic relaxation algorithm is adopted. This allows the efficient solution of the two nonlinear equations without significant drifting. Results show that robustness and efficiency are significantly improved. Herein, algorithms are described in detail. Significant testing is performed with imposed strains up to 100 for a carbon steel. This contributes to the robustness of equilibrium iterations. Drifting is also assessed as a function of number of steps in the dynamic relaxation algorithm. Numerical experimentation is performed in the 3D tension test for an anisotropic yield function.
AreiasPRosaPARRabczukT. Fully anisotropic hyperelasto-plasticity with exponential approximation by power series and scaling/squaring. Comput Mech2021; 68: 391–404.
2.
AreiasPRabczukTVazMF, et al. A consistent algorithm for finite-strain visco-hyperelasticity and visco-plasticity of amorphous polymers. Comp Method Appl M2022; 400: 115528.
3.
AreiasPDias-da CostaDPiresEB, et al. A new semi-implicit formulation for multiple-surface flow rules in multiplicative plasticity. Comput Mech2012; 49: 545–564.
4.
BonetJWoodRD. Nonlinear continuum mechanics for finite element analysis. 2nd ed. Cambridge: Cambridge University Press, 2008.
5.
PowellMJD. Chapter A hybrid method for nonlinear equations. In: Numerical methods for nonlinear algebraic equations. London, U.K: Gordon and Breach Science, 1970, pp.87–144.
6.
DennisJESchnabelRB. Number 16 in classicals in applied mathematics. In: Numerical methods for unconstrained optimization and nonlinear equations. Philadelphia, PA: SIAM, 1996.
7.
ConnARGouldNIMTointPL. MPS-SIAM series on optimization. In: Trust-region methods. Philadelphia, USA: SIAM - Society for Industrial and Applied Mathematics, 2000.
ErbtsPDüsterA. Accelerated staggered coupling schemes for problems of thermoelasticity at finite strains. Comput Math Appl2012; 64: 2408–2430.
11.
ZhouXLuDZhangY, et al. An open-source unconstrained stress updating algorithm for the modified Cam-clay model. Comp Method Appl M2022; 390: 114356.
12.
BennettKCRegueiroRALuscherDJ. Anisotropic finite hyper-elastoplasticity of geomaterials with Drucker-Prager/Cap type constitutive model formulation. Int J Plast2019; 123: 224–250.
13.
DegrooteJHaeltermanRAnnerelS, et al. Performance of partitioned procedures in fluid-structure interaction. Comput Struct2010; 88: 446–457.
14.
EidelBGruttmannF. Elastoplastic orthotropy at finite strains: Multiplicative formulation and numerical implementation. Comput Mater Sci2003; 28: 732–742.
15.
MandelJ. Equations constitutives et directeurs dans les milieux plastiques et viscoplastiques. Int J Solids Struct1973; 9: 725–740.
16.
SimoJC. Algorithms for static and dynamic multiplicative plasticity that preserve the classical return mapping schemes of the infinitesimal theory. Comp Method Appl M1992; 99: 61–112.
17.
CuitiñoAOrtizM. A material-independent method for extending stress update algorithms from small-strain plasticity to finite plasticity with multiplicative kinematics. Eng Computation1992; 9: 437–451.
18.
de Souza NetoEAPerićDOwenDRJ. Computational methods for plasticity: Theory and applications. West Sussex, England: John Wiley and Sons, Ltd., 2008.
19.
ShutovAV. Efficient time stepping for the multiplicative Maxwell fluid including the Mooney-Rivlin hyperelasticity. Int J Numer Meth Eng2018; 113: 1851–1869.
20.
MolerCVan LoanC. Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later. SIAM Rev2003; 45: 1–46.
21.
WardRC. Numerical computation of the matrix exponential with accuracy estimate. SIAM J Numer Anal1977; 14: 600–610.
22.
HighamNJ. The scaling and squaring method for the matrix exponential revisited. SIAM Rev2009; 51: 747–764.
23.
FanJ. Convergence rate of the trust region method for nonlinear equations under local error bound condition. Comput Optim Appl2005; 34: 215–227.
24.
ByrdRHKhalfanHFSchnabelRB. Analysis of a symmetric rank-one trust region method. SIAM J Opt1996; 6: 1025–1039.
25.
DeuflhardP. Newton methods for nonlinear problems. Affine invariance and adaptive algorithms. Heidelberg: Springer, 2004.
LublinerJ. Plasticity theory. New York, USA: Macmillan, 1990.
31.
SimoJCHughesTJR. Computational inelasticity. 2nd ed. New York: Springer, 2000.
32.
BessonJCailletaudGChabocheJ-L, et al. Non-linear mechanics of materials, volume 167of Solid Mechanics and its Applications. Springer, 2010.
33.
WriggersP. Nonlinear finite element methods. Berlin, Heidelberg: Springer, 2008.
34.
BelytschkoTLiuWKMoranB. Nonlinear finite elements for continua and structures. New York: John Wiley & Sons, 2000.
35.
OppermannPDenzerRMenzelA. A thermo-viscoplasticity model for metals over wide temperature ranges-application to case hardening steel. Comput Mech2022; 69: 541–563.
36.
SimoJCArmeroF. Geometrically non-linear enhanced strain mixed methods and the method of incompatible modes. Int J Numer Meth Eng1992; 33: 1413–1449.
37.
SimoJCArmeroFTaylorRL. Improved versions of assumed strain tri-linear elements for 3D finite deformation problems. Comp Method Appl M1993; 110: 359–386.
38.
HillR. A theory of yielding and plastic flow of anisotropic metals. Proc Roy Soc London1948; 193: 281–297.
39.
SimoJCTaylorRLPisterKS. Variational and projection methods for the volume constraint in finite deformation elasto-plasticity. Comp Method Appl M1985; 51: 177–208.
LesterBTScherzingerWM. Trust-region based return mapping algorithm for implicit integration of elastic-plastic constitutive models. Int J Numer Meth Eng2017; 112: 257–282.
42.
BickartTA. Matrix exponential: Approximation by truncated power series. Proc IEEE1968; 56: 372–373.
