Subdifferential‐based implicit return‐mapping operators in computational plasticity. Issue 11 (20th June 2016)

• Record Type: Journal Article
• Authors: Sysala, Stanislav; Cermak, Martin; Koudelka, Tomáš; Kruis, Jaroslav; Zeman, Jan; Blaheta, Radim
• Is Part Of: Zeitschrift für angewandte Mathematik und Mechanik. Volume 96:Issue 11(2016)
• Page Start: 1318
• Page End: 1338
• Publication Date: 2016-06-20
• Subjects: Elastoplasticity -- nonsmooth yield surface -- multivalued flow direction -- implicit return‐mapping scheme -- semismooth Newton method -- limit analysis -- 35Q74 -- 74C05 -- 74S05 -- 90C25; Mathematics -- Periodicals; Mechanics, Applied -- Periodicals; Engineering -- Periodicals; 519
• DOI: 10.1002/zamm.201500305
• Languages: English
• ISSNs: 0044-2267

Abstract

Abstract : The authors explore a numerical solution to elastoplastic constitutive initial value problems. An improved form of the implicit return‐mapping scheme for nonsmooth yield surfaces is proposed that systematically builds upon a subdifferential formulation of the flow rule. The main advantage of this approach is that the treatment of singular points – apices or edges at which the flow direction is multivalued – only involves a uniquely defined set of non‐linear equations, similarly to smooth yield surfaces. This paper focuses on isotropic models containing: a ) yield surfaces with one or two apices (singular points) on the hydrostatic axis, b ) plastic pseudo‐potentials that are independent of the Lode angle, and c ) possibly nonlinear isotropic hardening. Abstract : In this paper we explore a numerical solution to elastoplastic constitutive initial value problems. An improved form of the implicit return‐mapping scheme for nonsmooth yield surfaces is proposed that systematically builds upon a subdifferential formulation of the flow rule. The main advantage of this approach is that the treatment of singular points – apices or edges at which the flow direction is multivalued – only involves a uniquely defined set of non‐linear equations, similarly to smooth yield surfaces. This paper focuses on isotropic models containing: a ) yield surfaces with one or two apices (singular points) on the hydrostatic axis, b ) plastic pseudo‐potentials that are independent of the Lode angle, and c ) possibly nonlinear isotropic hardening. We show that for some models the improved integration scheme also enables us to a priori decide about a type of the return and to investigate the existence, uniqueness, and semismoothness of discretized constitutive operators. The semismooth Newton method is also introduced for solving the incremental boundary‐value problems. The paper contains numerical examples related to slope stability with publicly available Matlab implementations.