A robust returnmap algorithm for general multisurface plasticity
Three new contributions to the field of multisurface plasticity are presented for general situations with an arbitrary number of nonlinear yield surfaces with hardening or softening. A method for handling linearly dependent flow directions is described. A residual that can be used in a line search is defined. An algorithm that has been implemented and comprehensively tested is discussed in detail. Examples are presented to illustrate the computational cost of various components of the algorithm. The overall result is that a single NewtonRaphson iteration of the algorithm costs between 1.5 and 2 times that of an elastic calculation. Examples also illustrate the successful convergence of the algorithm in complicated situations. For example, without using the new contributions presented here, the algorithm fails to converge for approximately 50% of the trial stresses for a common geomechanical model of sedementary rocks, while the current algorithm results in complete success. Since it involves no approximations, the algorithm is used to quantify the accuracy of an efficient, pragmatic, but approximate, algorithm used for sedimentaryrock plasticity in a commercial software package. Furthermore, the main weakness of the algorithm is identified as the difficulty of correctly choosing the set of initially active constraints in the generalmore »
 Authors:

^{[1]};
^{[1]};
^{[2]};
^{[1]}
 CSIRO Energy Flagship, Queensland Centre for Advanced Technologies, Kenmore (Australia)
 Idaho National Lab. (INL), Idaho Falls, ID (United States)
 Publication Date:
 Report Number(s):
 INL/JOU1534725
Journal ID: ISSN 00295981
 Grant/Contract Number:
 AC0705ID14517
 Type:
 Accepted Manuscript
 Journal Name:
 International Journal for Numerical Methods in Engineering
 Additional Journal Information:
 Journal Name: International Journal for Numerical Methods in Engineering; Journal ID: ISSN 00295981
 Publisher:
 Wiley
 Research Org:
 Idaho National Lab. (INL), Idaho Falls, ID (United States)
 Sponsoring Org:
 USDOE
 Country of Publication:
 United States
 Language:
 English
 Subject:
 15 GEOTHERMAL ENERGY; multisurface plasticity; returnmap algorithm; KuhnTucker conditions; geomechanics
 OSTI Identifier:
 1259311
Adhikary, Deepak P., Jayasundara, Chandana T., Podgorney, Robert K., and Wilkins, Andy H.. A robust returnmap algorithm for general multisurface plasticity. United States: N. p.,
Web. doi:10.1002/nme.5284.
Adhikary, Deepak P., Jayasundara, Chandana T., Podgorney, Robert K., & Wilkins, Andy H.. A robust returnmap algorithm for general multisurface plasticity. United States. doi:10.1002/nme.5284.
Adhikary, Deepak P., Jayasundara, Chandana T., Podgorney, Robert K., and Wilkins, Andy H.. 2016.
"A robust returnmap algorithm for general multisurface plasticity". United States.
doi:10.1002/nme.5284. https://www.osti.gov/servlets/purl/1259311.
@article{osti_1259311,
title = {A robust returnmap algorithm for general multisurface plasticity},
author = {Adhikary, Deepak P. and Jayasundara, Chandana T. and Podgorney, Robert K. and Wilkins, Andy H.},
abstractNote = {Three new contributions to the field of multisurface plasticity are presented for general situations with an arbitrary number of nonlinear yield surfaces with hardening or softening. A method for handling linearly dependent flow directions is described. A residual that can be used in a line search is defined. An algorithm that has been implemented and comprehensively tested is discussed in detail. Examples are presented to illustrate the computational cost of various components of the algorithm. The overall result is that a single NewtonRaphson iteration of the algorithm costs between 1.5 and 2 times that of an elastic calculation. Examples also illustrate the successful convergence of the algorithm in complicated situations. For example, without using the new contributions presented here, the algorithm fails to converge for approximately 50% of the trial stresses for a common geomechanical model of sedementary rocks, while the current algorithm results in complete success. Since it involves no approximations, the algorithm is used to quantify the accuracy of an efficient, pragmatic, but approximate, algorithm used for sedimentaryrock plasticity in a commercial software package. Furthermore, the main weakness of the algorithm is identified as the difficulty of correctly choosing the set of initially active constraints in the general setting.},
doi = {10.1002/nme.5284},
journal = {International Journal for Numerical Methods in Engineering},
number = ,
volume = ,
place = {United States},
year = {2016},
month = {6}
}