Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 81
PROCEEDINGS OF THE TENTH INTERNATIONAL CONFERENCE ON CIVIL, STRUCTURAL AND ENVIRONMENTAL ENGINEERING COMPUTING Edited by: B.H.V. Topping
Paper 144
An Efficient Return Algorithm for NonAssociated MohrCoulomb Plasticity J. Clausen+, L. Damkilde+ and L. Andersen*
+Esbjerg Department of Engineering, Aalborg University Esbjerg, Esbjerg, Denmark
J. Clausen, L. Damkilde, L. Andersen, "An Efficient Return Algorithm for NonAssociated MohrCoulomb Plasticity", in B.H.V. Topping, (Editor), "Proceedings of the Tenth International Conference on Civil, Structural and Environmental Engineering Computing", CivilComp Press, Stirlingshire, UK, Paper 144, 2005. doi:10.4203/ccp.81.144
Keywords: stress update, return mapping, consistent constitutive matrix, MohrCoulomb yield criterion, nonlinear FEM.
Summary
The paper presents an improved method for updating the stresses in finite
element calculations involving MohrCoulomb plasticity. The method belongs to
the class of methods called "return mappings", or "backward Euler schemes",
which are often used for the integration of stress increments in plastic
material behaviour. The key point is a transformation from the general
sixdimensional stress space into the threedimensional principal stress space, where
the mathematical manipulations simplify those based on a geometrical interpretation.
The efficiency of the method is illustrated through numerical examples.
The returnmap algorithm starts from an elastic stress prediction based on the current stress state and a finite strain increment. The predictor stress is then returned to the yield surface through an iteration procedure in which the first and second derivatives of the yield surface with respect to the stress components are required, see e.g. [1]. In the general case with six independent stress components the calculation is somewhat cumbersome for MohrCoulomb plasticity, and a transformation into the principal stress directions simplifies the calculation as the criterion is formulated directly in the principal stresses. The key observation is that the return from the predictor stress will leave the principal stress directions unchanged during the iteration. The stress state at the yield surface and the corresponding material stiffness matrix can afterwards be transformed back to the original system by standard methods. The MohrCoulomb yield criterion consists in the general case of six independent yield surfaces, which gives intersection curves and an apex. The return of the predictor stress must identify the correct surface and this is very complicated in the general stress space whereas in the principal stress space it is fairly simple including the identification of intersections and apex. In order to have a quadratic rate of convergence of the global iteration the socalled consistent constitutive matrix is applied. The calculation is in the general stress space quite complicated, but in the principal stress space the calculations can be simplified significantly. For stress points at the intersection curves it is quite difficult in the general case to include two or more independent plastic strain directions, but again in the principal direction space it can be given a straight geometric interpretation. The method works in both associated and nonassociated plasticity, and numerical examples are presented to illustrate the efficiency. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
