Large Strain Elasto Plasticity

Juraj Šarloši, Jozef Bocko

American Journal of Mechanical Engineering

Large Strain Elasto Plasticity

Juraj Šarloši1,, Jozef Bocko1

1Department of Applied Mechanics and Mechanical Engineering, Technical University in Košice, Košice


In this paper we present a continuum theory for large strain Elasto-Plasticity based on formulations: Eulerian Formulation and Multiplicative Elasto-Plasticity. The theories includes Cauchy and Kirchhoff stress tensor as well as Truesdell rate and Jauman rate for the Cauchy stress and the multiplicative elasto-plastic decomposition. We show detailed derivative for mentioned formulations.

Cite this article:

  • Juraj Šarloši, Jozef Bocko. Large Strain Elasto Plasticity. American Journal of Mechanical Engineering. Vol. 4, No. 7, 2016, pp 344-348.
  • Šarloši, Juraj, and Jozef Bocko. "Large Strain Elasto Plasticity." American Journal of Mechanical Engineering 4.7 (2016): 344-348.
  • Šarloši, J. , & Bocko, J. (2016). Large Strain Elasto Plasticity. American Journal of Mechanical Engineering, 4(7), 344-348.
  • Šarloši, Juraj, and Jozef Bocko. "Large Strain Elasto Plasticity." American Journal of Mechanical Engineering 4, no. 7 (2016): 344-348.

Import into BibTeX Import into EndNote Import into RefMan Import into RefWorks

At a glance: Figures

1. Introduction

In recent decades improved numerical solutions of continua subjected to large inelastic strains have been obtained by changing the initial hypo-elastic-plastic algorithms to procedures based on hyperplastic formulation. Typically, a hyperelastic constitutive relation governs the elastic deformations, but the main purpose of this approach is to by-pass the integration of the stress rates in the plasticity part, rather than considering large elastic strains. Indeed, most approaches that have been presented assume that the elastic strains remain small, at least compared with the plastic In this article we will first outline the Eulerian approach, and then directly proceed with the multiplicative decomposition.

2. Eulerian Formulation

An Eulerian finite element formulation is most convenient developed starting from the virtual work expression in the current configuration :


Using Equations:


which is mathematically expression for conservation of mass of an elementary volume and other one


This identity can be also expressed in terms of Kirchhoff stress tenzor [5]:


which has advantage that the integration can now be carried out for the known volume . As with the Lagrange formulation we use the kinematic relation between the virtual strains and the virtual displacement, Equation


relation between the variation of the strain tensor and that of the continuous displacement field and with Equation


which is a continuous displacement field for field for all points within an element, to arrive at:


where for notational simplicity, the element index has been omitted. For general three-dimensional conditions we have [2]


with the inverse of the Jacobian matrix. Substitution of the discrete kinematics relation, and requiring that the resulting discrete equations hold for any virtual nodal displacement , yields the discrete equilibrium equations:

end, for the Cauchy stress tensor, the internal force vector reads:

Whereas for the Kirchhoff stress tensor we have:


The tangential stiffness matrix can be obtained in a standard manner, namely by differentiating the internal virtual work. For the Kirchhoff stress this results in [5]:


where the latter equality sign holds because of the symmetry of the Kirchhoff stress tensor. We now elaborate and first note because is fixed,


So that, using the velocity gradient as:

we get


Next, we choose the Truesdell rate of the Kirchhoff stress as the objective stress rate [1].

We can to express



- is material tangential stiffness tensor,

Note that the superscript continues to denote the transpose of a quantity. We next substitute Equations (12) into (10) to obtain:


Exploiting the symmetry of the Kirchhoff stress tensor this identity can be reworked to give:


This equation can be discretised in a straightforward manner, yielding:




the tangential stiffness matrix, and given in Equation (8), while and are given by [5]




Evidently, the choice of the Truesdell rate of the Kirchhoff stress yields a particularly simple expression for the tangential stiffness matrix. This is not so for other choice of the objective stress rate and the stress tensor. For instance, if we take the Jaumann rate in conjunction with Kirchhoff stress tensor, the resulting expression for the tangential stiffness matrix is [3]:


Where the superscript “JK” at the constitutive matrix denotes the Jaumann derivative of the Kirchhoff stress. Using of the term [4]


So that, in index notation, both constitutive matrices are related through:


where it is noted that expressing one objective stress rate into another through a modification of the constitutive matrix is always possible. Rewriting in matrix-vector format yields [5]:


which explains (12.16), with


Using a similar procedure the tangential stiffness matrices for other objective rates in conjuction with different stress measures can be derived. For instance, when using the Jaumann rate for the Cauchy stress, one arrives at




and as ,but with Cauchy stresses in lieu of Kirchhoff stresses. The matrix is derived from the relation between the Jaumann rates of the Cauchy and the Kirchhoff stress tensors, resulting in [1]:


In index notation, this equation can be reworked to give:


so that we finally obtain the matrix:


Indeed, the non-symmetry in can be balanced by a non-symmetry in , which, for instance, is the case when a hyperelastic constitutive relation is adopted [3].

By contrast, the use of the Truesdell rate of the Cauchy stress does not result in a non-symmetric tangential stiffness matrix. Substituting expression


into (a) and using the relation between the Cauchy and the Kirchhoff stress tensor, we yields the tangential stiffness matrix [4]:

Figure 1. Multiplicative elasto-plastic decomposition: (a) initial state, (b) intermediate state,(c)final state [5]

3. Multiplicative Elasto – Plasticity

The multiplicative elasto-plastic decomposition assumes the existence of three configurations: the initial, undeformed configuration, with a line segments , which is first moved, by a purely plastic deformation, into an intermediate configuration , and subsequently, into the final configuration through a pure elastic deformation (Figure 1). For each configuration we have [4]


for the mapping from the initial state to the intermediate configuration, with the plastic part of the deformation gradient, and


the mapping from the intermediate state to the final configuration, with the elastic part of the deformation gradient, which, in standard manner, can be decomposed into a rotational part [3].

and a contribution that stems from a pure deformation, ,as follows:


Subsequently, definition (33) can be used to define the right Cauchy – Green deformation tensor referred to the intermediate, elastic reference state ,


the ”elastic” Green – Lagrange strain tensor based upon ,


and the left Cauchy – Green deformation tensor referred to the intermediate, elastic reference state:


Considering the definition of the deformation gradient tensor and combining Equations (32) and (33) give the multiplicative decomposition of the deformation gradient for elasto – plastic deformations [1]:


For the one – dimensional case, Equation (38) particularizes as:


Whit the stretch ratio, which is multiplicative decomposed into an elastic and plastic stretch ratio, and , respectively [2].

The multiplicative decomposition for elasto – plasticity is not unique. For instance, it would be equally possible to rotate the intermediate configuration by , such that:


whence, using Equation (32) and (33) we straightforwardly arrive at the following, alternative elasto – plastic multiplicative decomposition:


which is equivalent to the decomposition of Equation (38), since [4]


For crystalline materials, it is physically reasonable to consider that the plastic deformation gradient purely represents the plastic sliding between crystals lattice and its rotation. This is represented by the decomposition of Equation (38) and is shown graphically in Figure 1.

Unfortunately, the Second Piola – Kirchhoff stress tensor cannot be related generally to the Green – Lagrange strain tensor than can be constructed on the basis of the elastic deformation gradient:


Since it is not invariant with respect to a rotation of the intermediate configuration. In particular, using Equation (42), one obtains:


and the stress is not affected by the frame in which the intermediate configuration is represented. The important consequence of the decomposition of Equation (38) is that, although the elastic and plastic deformations are decomposed in a multiplicative sense, this is not so for the strains rates [2].

Using the following relations


where, in view of the definition of the velocity gradient, the additively decomposed elastic and plastic velocity gradients, and , refer to the current configuration. It is emphasized that this additive decomposition of the velocity gradient depends crucially on the definition for as given in Equation (45). for instance, when [1]


is substituted for , which would then be similar to the definition of , an additive decomposition is not obtained. Nevertheless, the symmetric port of ,


Is a measure for the plastic stretching, as the eigenvalues of the spectral decomposition,


The anti – symmetrical part of is named the plastic spin tensor [3],


And represents the instantaneous rate of plastic spin of the intermediate configuration. In principle, a constitutive equation must be postulated for the plastic spin tensor, but in this treatment the hypothesis is made that the plastic spin vanishes:


This hypothesis holds rigorously for plastic isotropy, but not necessarily for plastic anisotropy.

Using the hypothesis of Equating (50) it directly follows that:

Using the definition of Equation (46), the plastic strain rate in the currently configuration can be written as:


and, with expression can be obtained in the intermediate configuration [5]:


with as in Equation (53). Clearly, there is no unequivocal definition for the plastic strain rate, neither in the current configuration, nor in the intermediate configuration. For instance, under the assumption that the elastic strains remain small, so that , Equation (53) can be approximated as [1]:


with using relations


we get:


When the elastic strains remain small, , and using Equation (34), and .

- plastic stretching, which is transformed to the deformed configuration through the elastic rotation [4].


I would like to thank to Technical School in Kosice Department of applied mechanics and mechanical engineering. This research was supported by grants from the Slovak Grant Agency VEGA No. 1/0731/16.


[1]  Š. Benča, “Riešenie nelineárnych pevnostných úloh pomocou MKP,” 1000 2013.
In article      
[2]  Š. Benča, “Aplikovaná nelineárna mechanika kontinua”. 1000
In article      
[3]  J. Bonet, A. J. Gill F. D. Wood, “Worked Examples in Nonlinear Continuum Mechanics for Finite Element Analysis,” Cambridge University Press,New York, 2012.
In article      View Article
[4]  D. G. Fetis, “Nonlinear Structural Engineering,” Springer-Verlag Berlin Heidelberg, 2006.
In article      
[5]  R. Borst, M. A. Crisfield, J. J. C. Remmers, C. V. Verhoosel, “Non-Linear Finite Element Analysis of Solid and Structures ,” J. Wiley and Sons Ltd. 2012.
In article      View Article
  • CiteULikeCiteULike
  • MendeleyMendeley
  • StumbleUponStumbleUpon
  • Add to DeliciousDelicious
  • FacebookFacebook
  • TwitterTwitter
  • LinkedInLinkedIn