Premium
Trust‐region based return mapping algorithm for implicit integration of elastic‐plastic constitutive models
Author(s) -
Lester B. T.,
Scherzinger W. M.
Publication year - 2017
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.5515
Subject(s) - robustness (evolution) , constitutive equation , quadratic equation , line search , trust region , mathematical optimization , algorithm , computer science , numerical integration , projection method , mathematics , finite element method , dykstra's projection algorithm , engineering , geometry , mathematical analysis , structural engineering , biochemistry , chemistry , computer security , radius , gene
Summary A new method for the solution of the non‐linear equations forming the core of constitutive model integration is proposed. Specifically, the trust‐region method that has been developed in the numerical optimization community is successfully modified for use in implicit integration of elastic‐plastic models. Although attention here is restricted to these rate‐independent formulations, the proposed approach holds substantial promise for adoption with models incorporating complex physics, multiple inelastic mechanisms, and/or multiphysics. As a first step, the non‐quadratic Hosford yield surface is used as a representative case to investigate computationally challenging constitutive models. The theory and implementation are presented, discussed, and compared with other common integration schemes. Multiple boundary value problems are studied and used to verify the proposed algorithm and demonstrate the capabilities of this approach over more common methodologies. Robustness and speed are then investigated and compared with existing algorithms. Through these efforts, it is shown that the utilization of a trust‐region approach leads to superior performance versus a traditional closest‐point projection Newton–Raphson method and comparable speed and robustness to a line search augmented scheme. Copyright © 2017 John Wiley & Sons, Ltd.