Abstract
In this paper we present a new approach for solving elastoplastic problems as second order cone complementarity problems (SOCCPs). Specially, two classes of elastoplastic problems, i.e. the J 2 plasticity problems with combined linear kinematic and isotropic hardening laws and the Drucker-Prager plasticity problems with associative or non-associative flow rules, are taken as the examples to illustrate the main idea of our new approach. In the new approach, firstly, the classical elastoplastic constitutive equations are equivalently reformulated as second order cone complementarity conditions. Secondly, by employing the finite element method and treating the nodal displacements and the plasticity multiplier vectors of Gaussian integration points as the unknown variables, we obtain a standard SOCCP formulation for the elastoplasticity analysis, which enables the using of general SOCCP solvers developed in the field of mathematical programming be directly available in the field of computational plasticity. Finally, a semi-smooth Newton algorithm is suggested to solve the obtained SOCCPs. Numerical results of several classical plasticity benchmark problems confirm the effectiveness and robustness of the SOCCP approach.
Similar content being viewed by others
References
Hill R (1950) The mathematical theory of plasticity. Oxford University Press, London
Han W, Reddy BD (1999) Plasticity: mathematical theory and numerical analysis. Springer, Berlin
Simo JC, Hughes TJR (1998) Computational inelasticity. Springer, New York
Luenberger DG (1973) Introduction to linear and nonlinear programming. Addison-Wesley, Reading
Maier G (1970) A matrix structural theory of piece-wise linear plasticity with interacting yield planes. Meccanica 6: 55–66
Martin JB, Reddy BD, Griffin TB (1987) Applications of mathematical programming concepts to incremental elastic-plastic analysis. Eng Struct 9: 171–176
Feijoo RA, Zouain N (1988) Formulations in rates and increments for elastic-plastic analysis. Int J Numer Methods Eng 26: 2031–2048
Reddy BD, Martin JB (1991) Algorithms for the solution of internal variable problems in plasticity. Comput Methods Appl Mech Eng 93: 253–273
Comi C, Corigliano A, Maier G (1991) Extremum properties of finite-step solutions in elastoplasticity with nonlinear mixed hardening. Int J Solids Struct 27(8): 965–981
Romano G, Rosati L, Marottide Sciarra M (1993) Variational principles for a class of finite step elastoplastic problems with non-linear mixed hardening. Comput Methods Appl Mech Eng 109: 293–314
Simo JC, Taylor RL (1986) A return mapping algorithm for plane stress elastoplasticity. Int J Numer Methods Eng 22: 649–670
Cuomo M, Contrafatto L (2000) Stress rate formulation for elastoplastic models with internal variables based on augmented Lagrangian regularization. Int J Solids Struct 37: 3935–3964
Zhang HW, He SY, Li XS, Wriggers P (2004) A new algorithm for numerical solution of 3D elastoplastic contact problems with orthotropic friction law. Comput Mech 34: 1–14
Krabbenhøft K, Lyamin AV, Sloan SW, Wriggers P (2007) an interior point algorithm for elastoplasticity. Int J Numer Methods Eng 69: 592–626
Zhu C (1995) A finite element-mathematical programming method for elastoplastic contact problems with friction. Finite Elem Anal Des 20: 273–282
Tin-Loi F, Pang JS (1993) Elastoplastic analysis of structures with nonlinear hardening: a nonlinear complementarity approach. Comput Methods Appl Mech Eng 107: 299–312
Christensen PW (2002) A non-smooth Newton method for elastoplastic problems. Comput Methods Appl Mech Eng 191: 1189–1219
Christian Wieners (2008) SQP methods for incremental plasticity with kinematic hardening. In: IUTAM symposium on theoretical, computational and modeling aspects of inelastic media. IUTAM book series, vol 11, Part 4, pp 143–153. Springer, Berlin
Hjiaj M, Fortin J, de Saxcé G (2003) A complete stress update algorithm for the non-associated Drucker-Prager model including treatment of the apex. Int J Eng Sci 41: 1109–1143
Ben-Tal A, Nemirovski A (2001) Lectures on modern convex optimization: analysis, algorithms, and engineering applications. SIAM, Philadelphia
Bisbos CD, Pardalos PM (2007) Second-order cone and semi-definite representations of material failure criteria. J Optim Theory Appl 134: 275–301
Krabbenhøft K, Lyamin AV, Sloan SW (2007) Formulation and solution of some plasticity problems as conic programs. Int J Solids Struct 44: 1533–1549
Krabbenhøft K, Lyamin AV, Sloan SW (2008) Three-dimensional Mohr–Coulomb limit analysis using semidefinite programming. Commun Numer Methods Eng 24: 1107–1119
Makrodimopoulos A (2006) Computational formulation of shakedown analysis as a conic quadratic optimization problem. Mech Res Commun 33: 72–83
Chen JS, Pan SH (2012) A survey on SOC complementarity functions and solution methods for SOCPs and SOCCPs. Pacific J Optim 8:33–74
Hayashi S (2004) Studies on second order cone complementarity problems. PhD thesis, Kyoto University
de Souza Neto EA, Perić D, Owen DRJ (2008) Computational methods for plasticity: theory and applications. Wiley, Chichester
Wanxie Z, Hongwu Z, Chengwei W (1997) Parametric variational principle and its application in engineering. Scientific and Technical Publishers, Beijing
Reddy BD, Wieners C, Wohlmuth B (2012) Finite element analysis and algorithms for single-crystal strain-gradient plasticity. Int J Numer Methods Eng. http://www.math.kit.edu/iwrmm/seite/preprints/media/preprint%20nr.%2011-4.pdf
Pan SH, Chen JS (2010) A semi-smooth Newton method for SOCCPs based on a one-parametric class of complementarity functions. Comput Optim Appl 45: 59–88
Fukushima M, Luo ZQ, Tseng P (2002) Smoothing functions for second-order cone complementarity problems. SIAM J Optim 12: 436–460
Hayashi S, Yamashita N, Fukushima M (2005) A combined smoothing and regularization method for monotone second-order cone complementarity problems. SIAM J Optim 15: 593–615
Chen JS, Tseng P (2005) An unconstrained smooth minimization reformulation of the second-order cone complementarity problem. Math Program 104: 293–327
Hager C, Wohlmuth B (2009) Nonlinear complementarity functions for plasticity problems with frictional contact. Comput Methods Appl Mech Eng 198: 3411–3427
Zhang HW, Li JY, Pan SH (2011) New second-order cone linear complementarity formulation and semi-smooth Newton algorithm for finite element analysis of 3D frictional contact problem. Comput Methods Appl Mech Eng 200: 77–88
Nocedal J, Wright SJ (1999) Numerical optimization. Springer, Berlin
De Saxce G, Bousshine L (1998) Limit analysis theorems for implicit standard materials: application to the unilateral contact with dry friction and the non associated flow rules in soils and rocks. Int J Mech Sci 40: 387–398
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Zhang, L.L., Li, J.Y., Zhang, H.W. et al. A second order cone complementarity approach for the numerical solution of elastoplasticity problems. Comput Mech 51, 1–18 (2013). https://doi.org/10.1007/s00466-012-0698-6
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00466-012-0698-6