Line 172: Line 172:
  
 
<span id="eq-12"></span>
 
<span id="eq-12"></span>
<span id="eq-13"></span>
 
 
{| class="formulaSCP" style="width: 100%; text-align: left;"  
 
{| class="formulaSCP" style="width: 100%; text-align: left;"  
 
|-
 
|-
Line 179: Line 178:
 
|-
 
|-
 
| style="text-align: center;" | <math>D = \left(\bar{\boldsymbol{M}}-    {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}\right):\bar{\boldsymbol{L}}_{\mathrm{p}} </math>
 
| style="text-align: center;" | <math>D = \left(\bar{\boldsymbol{M}}-    {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}\right):\bar{\boldsymbol{L}}_{\mathrm{p}} </math>
 +
|}
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (12)
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (12)
|-
+
|}
| style="text-align: left;" |
+
  
 
where the kinematic hardening stress of Mandel type (also denoted back-stress) is defined as
 
where the kinematic hardening stress of Mandel type (also denoted back-stress) is defined as
 
+
<span id="eq-13"></span>
 +
{| class="formulaSCP" style="width: 100%; text-align: left;"
 +
|-
 +
|
 +
{| style="text-align: left; margin:auto;width: 100%;"
 +
|-
 +
| style="text-align: center;" |
 
<math>    {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}= \frac{\partial }{\partial \boldsymbol{F}_{\mathrm{ke}}}\boldsymbol{F}^{\mathrm{t}}_{\mathrm{ke}} </math>
 
<math>    {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}= \frac{\partial }{\partial \boldsymbol{F}_{\mathrm{ke}}}\boldsymbol{F}^{\mathrm{t}}_{\mathrm{ke}} </math>
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (13)
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (13)

Revision as of 12:22, 13 May 2017

A Comparison of two frameworks for kinematic hardening in hyperelasto-plasticity

Knut A. MeyerFailed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): ^{*} , Magnus Ekh

Department of Applied Mechanics, Division of Material and Computational Mechanics,
Chalmers University of Technology, 412 96 Gothenburg, Sweden
e-mailFailed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): ^{*}
knut.andreas.meyer@chalmers.se, web page: http://www.chalmers.se/

Abstract

In this work we compare two frameworks for thermodynamically consistent hyperelasto-plasticity with kinematic hardening. The first was formulated by Dettmer and Reese (2004), inspired by Lion (2000), and has been used to model sheet metal forming. The second, formulated by Wallin et al. (2003), has been used to model large shear strains and cyclic ratcheting behavior of pearlitic steel (Johansson et al. 2006). In this paper we show that these frameworks can result in equivalent models for certain choices of free energies. Furthermore, it is shown that the choices of free energy found in the literature only result in minor differences. These differences are discussed theoretically and investigated numerically.

keywords Hyperelasto-plasticity, Finite strains, Kinematic hardening


1 INTRODUCTION

Large strains in metals during room temperature occur in many technical applications, often during manufacturing, such as sheet metal forming. Some components are also subjected to large strains during service, for example in the surface layer of railway rails and wheels (see e.g. [1, 2]). Experiments have shown that the Bauchinger effect, often modeled with kinematic hardening, is pronounced in many metals. Kinematic hardening can be modeled with different thermodynamically consistent hyperelasto-plastic frameworks found in the literature, and two of them are considered here. The first framework is based on rheological models with an Armstrong-Frederick (AF) type of kinematic hardening, and was proposed by Lion [3] and further developed by Dettmer and Reese [4]. The second framework, introduced by Wallin et al. [5], also features an AF type of kinematic hardening, and has been used to model the Swift effect [6] and large deformations in railway applications [7, 8]. In this paper we compare these frameworks, both theoretically and numerically.

2 DESCRIPTION OF FRAMEWORKS

In this section, the modeling frameworks are presented in terms of their assumed kinematics and thermodynamics. The common parts are presented in Subsections 2.1 and 2.2, followed by a description of how the frameworks differ in Subsections 2.3 and 2.4. Three specific models are defined in Subsections 2.5 and 2.6, which are compared in the numerical examples in Section 3. For the clarity of the presentation, only linear kinematic hardening is considered in the current section. In Section 4 we investigate nonlinear kinematic hardening, e.g. of Armstrong-Frederick type.

File:COMPLAS FIG1.pdf
Figure 1: Configurations and deformation gradients

2.1 Kinematics and notations

Figure 1 shows the different configurations used in both [4] and [5]. Dettmer and Reese [4] introduce the inelastic plastic deformation gradient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{kp}}}

connecting the fictitious kinematic configuration Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\bar{\Omega }}}
to the initial configuration Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Omega }

. This connection is not introduced in Wallin et al. [5], but otherwise the same configurations and remaining deformation gradients are present in both frameworks.

Tensors on the current configuration Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \omega }

are denoted with lower case letters and no bars, e.g. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{b}}

. Tensors on the intermediate Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\Omega }}

and kinematic Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\bar{\Omega }}}
configurations are denoted by one bar, e.g. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{C}}_{\mathrm{e}}}

, and two bars, e.g. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\bar{\boldsymbol{C}}}_{\mathrm{ke}}} , respectively. The following decompositions of the deformation gradients and definitions of the deformation tensors will be used:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \begin{array}{rlrlrl}\boldsymbol{F} &= \boldsymbol{F}_{\mathrm{e}}\boldsymbol{F}_{\mathrm{p}} \quad &\bar{\boldsymbol{C}}_{\mathrm{e}}&= \boldsymbol{F}^{\mathrm{t}}_{\mathrm{e}}\boldsymbol{F}_{\mathrm{e}} \quad &\boldsymbol{c}_{\mathrm{e}}&= \boldsymbol{F}^{-\mathrm{t}}_{\mathrm{e}}\boldsymbol{F}^{-1}_{\mathrm{e}}= \boldsymbol{b}^{-1}_{\mathrm{e}} \\[1mm] \boldsymbol{F}_{\mathrm{p}}&= \boldsymbol{F}_{\mathrm{ke}}\boldsymbol{F}_{\mathrm{kp}} \quad &\bar{\bar{\boldsymbol{C}}}_{\mathrm{ke}}&= \boldsymbol{F}^{\mathrm{t}}_{\mathrm{ke}}\boldsymbol{F}_{\mathrm{ke}} \quad &\bar{\boldsymbol{c}}_{\mathrm{ke}}&= \boldsymbol{F}^{-\mathrm{t}}_{\mathrm{ke}}\boldsymbol{F}^{-1}_{\mathrm{ke}}= \bar{\boldsymbol{b}}_{\mathrm{ke}}^{-1} \end{array}
(1)

The velocity gradients on the intermediate and kinematic configurations are defined as

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \bar{\boldsymbol{L}}_{\mathrm{p}}= \dot{\boldsymbol{F}}_{\mathrm{p}}\boldsymbol{F}^{-1}_{\mathrm{p}} (2)
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \bar{\boldsymbol{L}}_{\mathrm{ke}}= \dot{\boldsymbol{F}}_{\mathrm{ke}}\boldsymbol{F}^{-1}_{\mathrm{ke}} (3)
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \bar{\bar{\boldsymbol{L}}}_{\mathrm{kp}}= \dot{\boldsymbol{F}}_{\mathrm{kp}}\boldsymbol{F}^{-1}_{\mathrm{kp}} (4)

and hence the velocity gradient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{l}}

on the current configuration can be written as

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \boldsymbol{l} = \dot{\boldsymbol{F}}\boldsymbol{F}^{-1}=\dot{\boldsymbol{F}}_{\mathrm{e}}\boldsymbol{F}^{-1}_{\mathrm{e}}+ \boldsymbol{F}_{\mathrm{e}}\bar{\boldsymbol{L}}_{\mathrm{p}}\boldsymbol{F}^{-1}_{\mathrm{e}}
(5)

2.2 Thermodynamics

The free energy Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi }

is introduced with the additive split according to

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \Psi = \Psi _{\mathrm{e}}\left(\bar{\boldsymbol{C}}_{\mathrm{e}}\right)+ \Psi _{\mathrm{kin}}\left(\boldsymbol{F}_{\mathrm{ke}}\right)
(6)

whereby the dissipation inequality (see e.g. Simo (1998) [9]) becomes

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): D = \boldsymbol{\tau }:\boldsymbol{l} - \dot{\Psi } = \boldsymbol{\tau }:\boldsymbol{l} - \frac{\partial }{\partial \bar{\boldsymbol{C}}_{\mathrm{e}}}\dot{\bar{\boldsymbol{C}}}_{\mathrm{e}}- \frac{\partial }{\partial \boldsymbol{F}_{\mathrm{ke}}}:\dot{\boldsymbol{F}}_{\mathrm{ke}}\geq 0
(7)

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{\tau }}

is the Kirchhoff stress.  Using the requirement of zero dissipation during elastic loading and Equation 5, the reduced dissipation inequality becomes

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): D = \bar{\boldsymbol{M}}:\bar{\boldsymbol{L}}_{\mathrm{p}}- \left(\frac{\partial }{\partial \boldsymbol{F}_{\mathrm{ke}}}\boldsymbol{F}^{\mathrm{t}}_{\mathrm{ke}}\right):\bar{\boldsymbol{L}}_{\mathrm{ke}}
(8)

where the Mandel stress Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{M}}}

is defined as

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \bar{\boldsymbol{M}} = 2\bar{\boldsymbol{C}}_{\mathrm{e}}\frac{\partial }{\partial \bar{\boldsymbol{C}}_{\mathrm{e}}}
(9)

In a standard fashion, we adopt an associative evolution of the plastic deformation gradient in this paper:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \bar{\boldsymbol{L}}_{\mathrm{p}}= \dot{\lambda }\frac{\partial }{\partial \bar{\boldsymbol{M}}}
(10)

where the functional dependence of the yield function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Phi \leq 0}

will be specified later. Up until this point the two frameworks are identical. We first describe the framework by Dettmer and Reese [4], before proceeding with the framework by Wallin et al. [5].

2.3 1st framework [4]

In the first framework proposed by Lion [3] and further developed by Dettmer and Reese [4], the plastic deformation gradient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{p}}}

is multiplicatively decomposed into an elastic part Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{ke}}}
and a plastic part Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{kp}}}

. The physical motivation is that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{ke}}}

represents local elastic deformations on the microscale caused by dislocations and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{kp}}}
represents irreversible displacements in the slip systems. The assumption is that development of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{ke}}}
results in linear kinematic hardening and the development of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{kp}}}
reduces Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{ke}}}

, hence causing saturation (dynamic recovery) of the kinematic hardening. This is illustrated using a rheological model by Lion [3]. The multiplicative split of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{p}}}

results in the following additive split of the plastic velocity gradient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{L}}_{\mathrm{p}}}


Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \bar{\boldsymbol{L}}_{\mathrm{p}}= \bar{\boldsymbol{L}}_{\mathrm{ke}}+ \boldsymbol{F}_{\mathrm{ke}}\bar{\bar{\boldsymbol{L}}}_{\mathrm{kp}}\boldsymbol{F}^{-1}_{\mathrm{ke}}
(11)

For the case of purely linear hardening (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{kp}}=\boldsymbol{I}} ), the reduced dissipation, Equation 8, can be written as

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): D = \left(\bar{\boldsymbol{M}}- {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}\right):\bar{\boldsymbol{L}}_{\mathrm{p}}
(12)

where the kinematic hardening stress of Mandel type (also denoted back-stress) is defined as

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}= \frac{\partial }{\partial \boldsymbol{F}_{\mathrm{ke}}}\boldsymbol{F}^{\mathrm{t}}_{\mathrm{ke}}

(13)

This motivates that the driving force for plastic flow is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{M}}- {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}}

and thereby a yield criterion expressed as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Phi (\bar{\boldsymbol{M}}-   {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{})}

.

2.4 2nd framework [5]

In the second framework, proposed by Wallin et al. [5], the deformation gradient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}^{-1}_{\mathrm{ke}}}

is introduced to model the deformation of the crystal lattice, due to the residual micro stresses responsible for the Bauchinger effect.  From this deformation gradient the kinematic hardening stress of Mandel type is defined as

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^2}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}= -\frac{\partial }{\partial \boldsymbol{F}_{\mathrm{ke}}}\boldsymbol{F}^{\mathrm{t}}_{\mathrm{ke}}
(14)

which yields that the reduced dissipation inequality 8 is

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): D = \bar{\boldsymbol{M}}:\bar{\boldsymbol{L}}_{\mathrm{p}} + {}\mathrm{^2}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}:\bar{\boldsymbol{L}}_{\mathrm{ke}}
(15)

Using the standard interpretation of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {}\mathrm{^2}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}}

as a back-stress that reduces the driving force for plasticity, motivates the yield function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Phi (\bar{\boldsymbol{M}}-   {}\mathrm{^2}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{})}

. This gives, by the postulate of maximum dissipation, the kinematic relation

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \bar{\boldsymbol{L}}_{\mathrm{p}}= - \bar{\boldsymbol{L}}_{\mathrm{ke}}
(16)

and the same reduced dissipation inequality as in Equation 12 is obtained:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): D = \left(\bar{\boldsymbol{M}}- {}\mathrm{^2}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}\right):\bar{\boldsymbol{L}}_{\mathrm{p}}
(17)

2.5 Specific formats for free energy

The elastic and kinematic free energies (with the third invariant Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I_{\mathrm{3\bullet }} = \det (\bullet )} ) proposed by Vladimirov et al [10] are

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^A}\Psi _{\mathrm{e}}\mathrm{} = \frac{1}{2}G\left(\mathrm{tr}(\bar{\boldsymbol{C}}_{\mathrm{e}}) - 3 - \ln \left(I_{\mathrm{3C}_{\mathrm{e}}}\right)\right) + \frac{\Lambda }{4}\left(I_{\mathrm{3C}_{\mathrm{e}}} - 1 - \ln \left(I_{\mathrm{3C}_{\mathrm{e}}}\right)\right) (18)
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^A}\Psi _{\mathrm{kin}}\mathrm{} = \frac{1}{2}H_{\mathrm{kin}}\left(\mathrm{tr}(\bar{\bar{\boldsymbol{C}}}_{\mathrm{ke}}) - 3 - \ln \left(I_{\mathrm{3C}_{\mathrm{ke}}}\right)\right) (19)

The part of the elastic free energy, corresponding to Lamé's second parameter Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Lambda } , is thus not included in the kinematic free energy.

A similar split is introduced in [5], but the free energy in that work is decomposed into an isochoric and a volumetric part. The formulation for the volumetric part is not the same in [5] and [7, 8], and here we use the formulation from [7, 8]. This difference only affects the bulk elastic response and the influence on the numerical results studied in this paper is therefore negligible.

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^B}\Psi _{\mathrm{e}}\mathrm{} = \frac{1}{2}G\left(\mathrm{tr}\left(I_{\mathrm{3C}_{\mathrm{e}}}^{-1/3} \bar{\boldsymbol{C}}_{\mathrm{e}}\right)- 3\right)+ \frac{1}{2}K\left(I_{\mathrm{3C}_{\mathrm{e}}}^{1/2} - 1 \right)^2 (20)
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^B}\Psi _{\mathrm{kin}}\mathrm{} = \frac{1}{2}H_{\mathrm{kin}}\left(\mathrm{tr}\left(I_{\mathrm{3c}_{\mathrm{ke}}}^{-1/3} \bar{\boldsymbol{c}}_{\mathrm{ke}}\right)- 3\right) (21)

From the discussion so far, there seem to be several differences between the frameworks: (1) the definition of the Mandel back-stress (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}}

or Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle    {}\mathrm{^2}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}}

), (2) the variable of which the kinematic free energy depends on (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\bar{\boldsymbol{C}}}_{\mathrm{ke}}}

or Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{c}}_{\mathrm{ke}}}

) and (3) what part of the elastic free energy formulation that is used to formulate the kinematic free energy. The third of these can be investigated by taking the format of free energy from the second framework, but using the definitions and variables from the first framework to obtain model Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{C}}

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^C}\Psi _{\mathrm{e}}\mathrm{} = {}\mathrm{^B}\Psi _{\mathrm{e}}\mathrm{} (22)
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^C}\Psi _{\mathrm{kin}}\mathrm{} = \frac{1}{2}H_{\mathrm{kin}}\left(\mathrm{tr}\left(I_{\mathrm{3C}_{\mathrm{ke}}}^{-1/3} \bar{\bar{\boldsymbol{C}}}_{\mathrm{ke}}\right)- 3\right) (23)

2.6 Stresses for each model

We have now described both frameworks, the first by Dettmer and Reese [4] and the second by Wallin et al. [5]. By letting Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _{\mathrm{kin}}}

depend on Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\bar{\boldsymbol{C}}}_{\mathrm{ke}}}
or Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{c}}_{\mathrm{ke}}}
we can use 13 and 14, respectively, to obtain the Mandel back-stresses for the two frameworks:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}= \quad \left(\frac{\partial }{\partial \bar{\bar{\boldsymbol{C}}}_{\mathrm{ke}}}:\frac{\partial }{\partial \boldsymbol{F}_{\mathrm{ke}}}\right)\boldsymbol{F}^{\mathrm{t}}_{\mathrm{ke}} = 2\boldsymbol{F}_{\mathrm{ke}}\left(\frac{\partial }{\partial \bar{\bar{\boldsymbol{C}}}_{\mathrm{ke}}} \right)\boldsymbol{F}^{\mathrm{t}}_{\mathrm{ke}} (24)
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^2}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}= -\left(\frac{\partial }{\partial \bar{\boldsymbol{c}}_{\mathrm{ke}}}:\frac{\partial }{\partial \boldsymbol{F}_{\mathrm{ke}}}\right)\boldsymbol{F}^{\mathrm{t}}_{\mathrm{ke}} = 2\bar{\boldsymbol{c}}_{\mathrm{ke}}\frac{\partial }{\partial \bar{\boldsymbol{c}}_{\mathrm{ke}}} (25)

The Mandel stresses for model Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{A}} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {}\mathrm{^A}\bar{\boldsymbol{M}}\mathrm{}}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle    {}\mathrm{^A}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}}

, are found using Equations 9 and 24 with the free energies in Equations 18 and 19:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^A}\bar{\boldsymbol{M}}\mathrm{} = G (\bar{\boldsymbol{C}}_{\mathrm{e}}- \boldsymbol{I}) + \frac{\Lambda }{2} (I_\mathrm{3C_e} - 1)\boldsymbol{I}, \quad {}\mathrm{^A}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{} = H_{\mathrm{kin}}\left(\bar{\boldsymbol{b}}_{\mathrm{ke}}-\boldsymbol{I}\right) (26)

The stresses for model Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \mathrm{B} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^B}\bar{\boldsymbol{M}}\mathrm{}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":):    {}\mathrm{^B}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}

, are given by using Equations 9 and 25 with the free energies in Equations 20 and 21:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^B}\bar{\boldsymbol{M}}\mathrm{} = G I_\mathrm{3C_e}^{-1/3} \bar{\boldsymbol{C}}_{\mathrm{e}}^{\mathrm{dev}}+ K(I_\mathrm{3C_e} - I_\mathrm{3C_e}^{1/2})\boldsymbol{I}, \quad {}\mathrm{^B}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{} = H_{\mathrm{kin}}I_{\mathrm{c}_{\mathrm{ke}}}^{-1/3} \bar{\boldsymbol{c}}_{\mathrm{ke}}^{\mathrm{dev}}

(27)

For model Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \mathrm{C} , we use the first framework, i.e. the stresses Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^C}\bar{\boldsymbol{M}}\mathrm{}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":):    {}\mathrm{^C}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}
are given by using Equations 9 and 24, but with the free energies in Equations 22 and 23. We further note that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): I_{3\mathrm{C}_{\mathrm{ke}}}=I_{3\mathrm{b}_{\mathrm{ke}}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \mathrm{tr}(\bar{\bar{\boldsymbol{C}}}_{\mathrm{ke}})=\mathrm{tr}(\bar{\boldsymbol{b}}_{\mathrm{ke}})

, which leads to:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^C}\bar{\boldsymbol{M}}\mathrm{} = {}\mathrm{^B}\bar{\boldsymbol{M}}\mathrm{} \quad {}\mathrm{^C}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{} = H_{\mathrm{kin}}I_{\mathrm{b}_{\mathrm{ke}}}^{-1/3} \bar{\boldsymbol{b}}_{\mathrm{ke}}^{\mathrm{dev}}

(28)

If Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {}\mathrm{^B}\bar{\boldsymbol{c}}_{\mathrm{ke}}\mathrm{} = {}\mathrm{^C}\bar{\boldsymbol{b}}_{\mathrm{ke}}\mathrm{}}

then clearly model Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{B}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{C}}
are equivalent. Assuming that is the case for some point in time, we also have Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle    {}\mathrm{^B}\bar{\boldsymbol{L}}_{\mathrm{p}}\mathrm{}=   {}\mathrm{^C}\bar{\boldsymbol{L}}_{\mathrm{p}}\mathrm{}=\bar{\boldsymbol{L}}_{\mathrm{p}}}

. Since model Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{B}}

is using the second framework, and model Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{C}}
is using the first, we also have Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle -   {}\mathrm{^B}\bar{\boldsymbol{L}}_{\mathrm{ke}}\mathrm{}=   {}\mathrm{^C}\bar{\boldsymbol{L}}_{\mathrm{ke}}\mathrm{}=\bar{\boldsymbol{L}}_{\mathrm{p}}}

, hence

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^B}\dot{\bar{\boldsymbol{c}}}_{\mathrm{ke}}\mathrm{} = -\left( {}\mathrm{^B}\bar{\boldsymbol{L}}\mathrm{_{\mathrm{ke}}^{\mathrm{t}}} {}\mathrm{^B}\bar{\boldsymbol{c}}_{\mathrm{ke}}\mathrm{} + {}\mathrm{^B}\bar{\boldsymbol{c}}_{\mathrm{ke}}\mathrm{} {}\mathrm{^B}\bar{\boldsymbol{L}}_{\mathrm{ke}}\mathrm{}\right) = \bar{\boldsymbol{L}}_{\mathrm{p}}^{\mathrm{t}} {}\mathrm{^B}\bar{\boldsymbol{c}}_{\mathrm{ke}}\mathrm{} + {}\mathrm{^B}\bar{\boldsymbol{c}}_{\mathrm{ke}}\mathrm{}\bar{\boldsymbol{L}}_{\mathrm{p}} (29)
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^C}\dot{\bar{\boldsymbol{b}}}_{\mathrm{ke}}\mathrm{} = \quad \left( {}\mathrm{^C}\bar{\boldsymbol{L}}_{\mathrm{ke}}\mathrm{} {}\mathrm{^C}\bar{\boldsymbol{b}}_{\mathrm{ke}}\mathrm{} + {}\mathrm{^C}\bar{\boldsymbol{b}}_{\mathrm{ke}}\mathrm{} {}\mathrm{^C}\bar{\boldsymbol{L}}\mathrm{_{\mathrm{ke}}^{\mathrm{t}}}\right)= \bar{\boldsymbol{L}}_{\mathrm{p}} {}\mathrm{^C}\bar{\boldsymbol{b}}_{\mathrm{ke}}\mathrm{} + {}\mathrm{^C}\bar{\boldsymbol{b}}_{\mathrm{ke}}\mathrm{}\bar{\boldsymbol{L}}_{\mathrm{p}}^{\mathrm{t}} (30)

As Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {}\mathrm{^B}\bar{\boldsymbol{c}}_{\mathrm{ke}}\mathrm{} = {}\mathrm{^C}\bar{\boldsymbol{b}}_{\mathrm{ke}}\mathrm{} = \boldsymbol{I}}

initially, the statement Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle    {}\mathrm{^B}\bar{\boldsymbol{c}}_{\mathrm{ke}}\mathrm{} =    {}\mathrm{^C}\bar{\boldsymbol{b}}_{\mathrm{ke}}\mathrm{}}
is true for all points in time under the assumption that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{L}}_{\mathrm{p}}}
is symmetric. If the free energy is isotropic, the Mandel stresses Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{M}}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{M}}_{\mathrm{k}}}
are symmetric, and hence Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{L}}_{\mathrm{p}}}
becomes symmetric for the associative choice of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{L}}_{\mathrm{p}}}
in Equation 10. This leads to the conclusion that model B and C are equivalent, which is verified numerically later. Furthermore, this proof leads to the interesting conclusion that the frameworks can give exactly the same model with proper choices of free energy. The possibility to formulate an unsymmetric Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{L}}_{\mathrm{p}}}
for isotropy is discussed in e.g. Wallin et al. [5] and Wallin and Ristinmaa [6], but is not investigated in this paper.

3 NUMERICAL RESULTS

In this section we evaluate the response of the material models for uniaxial loading and simple shear loading. The models are implemented using a standard Backward Euler integration scheme, for which Vladimirov et al. [10] noted that the accuracy suffers at large time steps. To avoid these accuracy problems, approximately Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 4\cdot{10}^4}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 5\cdot 10^4}
load steps are used for the uniaxial and simple shear loading, respectively.

The von Mises effective stress is used to define the yield function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Phi }

according to
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \Phi = \sqrt{\frac{3}{2}}\sqrt{\mathrm{dev}\left(\bar{\boldsymbol{M}}^\mathrm{t}-\bar{\boldsymbol{M}}_{\mathrm{k}}^\mathrm{t}\right):\mathrm{dev}\left(\bar{\boldsymbol{M}}-\bar{\boldsymbol{M}}_{\mathrm{k}}\right)} - Y_0 = f\left(\bar{\boldsymbol{M}}-\bar{\boldsymbol{M}}_{\mathrm{k}}\right)- Y_0 \leq 0
(31)

whereby the evolution of the plastic deformation gradient in 10 becomes

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \bar{\boldsymbol{L}}_{\mathrm{p}}= \dot{\lambda }\frac{3}{2} \frac{\mathrm{dev}\left(\bar{\boldsymbol{M}}^\mathrm{t}-\bar{\boldsymbol{M}}_{\mathrm{k}}^\mathrm{t}\right)}{f\left(\bar{\boldsymbol{M}}-\bar{\boldsymbol{M}}_{\mathrm{k}}\right)}
(32)

From this it follows that the plastic deformation is isochoric: By time differentiation of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \det \left(\boldsymbol{F}_{\mathrm{p}}\right)}

and using 2
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \frac{\partial }{\partial t} \det \left(\boldsymbol{F}_{\mathrm{p}}\right) = \det \left(\boldsymbol{F}_{\mathrm{p}}\right)\boldsymbol{F}^{-\mathrm{t}}_{\mathrm{p}}:\left(\bar{\boldsymbol{L}}_{\mathrm{p}}\boldsymbol{F}_{\mathrm{p}}\right) = \det \left(\boldsymbol{F}_{\mathrm{p}}\right)\mathrm{tr}\left(\bar{\boldsymbol{L}}_{\mathrm{p}}\right)= 0
(33)

hence, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \det \left(\boldsymbol{F}_{\mathrm{p}}\right)= 1} .

The following material parameters: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle G=81\,\mathrm{GPa}} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle K=174\,\mathrm{GPa}} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Lambda = K - 2G/3} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle Y_0=100\,\mathrm{MPa}}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle H_{\mathrm{kin}}=1000\,\mathrm{MPa}}
, are used in the numerical examples in Figure 2.
Uniaxial response Simple shear response
(a) Uniaxial response (b) Simple shear response
Figure 2: Numerical results

We first consider uniaxial stress in Figure 2a, by letting the normal deformation gradient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle F_{11}}

increase from 1 to 10 while keeping the Cauchy stresses Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \sigma _{22}=\sigma _{33}=0}
(Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{\sigma }=\boldsymbol{\tau }/I_{3\mathrm{F}}}

). As previously shown theoretically, model Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{B}}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{C}}
give the same response. Model Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{A}}
gives a somewhat stiffer response at large deformations, due to the different choice of free energy.

In the case of simple shear loading, Figure 2b shows that the response of all the models coincide. While this is expected for model Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{B}}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{C}}

, the fact that model Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{A}}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathrm{C}}
coincide is explained with simple shear being an isochoric process (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I_{3\mathrm{F}}=I_{3\mathrm{F}_{\mathrm{p}}}=1 \Rightarrow I_{3\mathrm{F}_{\mathrm{e}}}=1}

)

The results in Figure 2 show negligible differences between the different formulations of free energy. From a theoretical point of view, one could argue that the deviatoric dependence of the back-stress is more correct, based on the experimental evidence of volume preservation for metal plasticity.

4 ARMSTRONG-FREDERICK SATURATION

Linear kinematic hardening was considered in Section 2, for which Dettmer and Reese [4] set Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{kp}}=\boldsymbol{I}} , yielding the reduced dissipation inequality in Equation 12. If the general case with an evolving Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{F}_{\mathrm{kp}}}

is considered, the reduced dissipation inequality, using Equation 11, becomes
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): D = \left(\bar{\boldsymbol{M}}- {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}\right):\bar{\boldsymbol{L}}_{\mathrm{p}}+ {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}:\left(\boldsymbol{F}_{\mathrm{ke}}\bar{\bar{\boldsymbol{L}}}_{\mathrm{kp}}\boldsymbol{F}^{-1}_{\mathrm{ke}}\right)
(34)

From this the kinematic stress of Mandel type on the kinematic configuration Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {}\mathrm{^1}\bar{\bar{\boldsymbol{M}}}_{\mathrm{k}}\mathrm{}}

is defined as

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {}\mathrm{^1}\bar{\bar{\boldsymbol{M}}}_{\mathrm{k}}\mathrm{} = \boldsymbol{F}^{\mathrm{t}}_{\mathrm{ke}} {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}\boldsymbol{F}^{-\mathrm{t}}_{\mathrm{ke}}
(35)

Saturation is motivated by the rheological model, setting the evolution on Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\bar{\boldsymbol{L}}}_{\mathrm{kp}}^{\mathrm{sym}}} . As the Mandel stresses in [4] are symmetric, this can be written as

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \bar{\bar{\boldsymbol{L}}}_{\mathrm{kp}}= \dot{\lambda }\frac{ {}\mathrm{^1}\bar{\bar{\boldsymbol{M}}}^\mathrm{t}_{\mathrm{k}}\mathrm{}}{b_\infty } (36)

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): b_\infty

is a material parameter controlling the kinematic saturation. Equation 11 and 35 then yield

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \bar{\boldsymbol{L}}_{\mathrm{ke}}= \bar{\boldsymbol{L}}_{\mathrm{p}}- \dot{\lambda }\frac{ {}\mathrm{^1}\bar{\boldsymbol{M}}_{\mathrm{k}}^\mathrm{t}\mathrm{}}{b_\infty }


This Equation can be compared with Wallin et al. [5], who use a modified potential Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \Phi ^*_{\mathrm{kin}}

to obtain

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \bar{\boldsymbol{L}}_{\mathrm{ke}}= \dot{\lambda }\frac{\partial }{\partial \left( {}\mathrm{^2}\bar{\boldsymbol{M}}_{\mathrm{k}}\mathrm{}\right)} = -\bar{\boldsymbol{L}}_{\mathrm{p}}+ \dot{\lambda }\frac{ {}\mathrm{^2}\bar{\boldsymbol{M}}_{\mathrm{k}}^\mathrm{t}\mathrm{}}{b_\infty }

(37)

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{L}}_{\mathrm{ke}}}

will be symmetric if a modified yield potential Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Phi ^*_{\mathrm{kin}}}
exists and we have, as before, symmetric Mandel stresses. This is the case for the considered model with Armstrong-Frederick type of non-linear kinematic saturation. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{L}}_{\mathrm{ke}}}
is also symmetric in the work by Zhu et al. [11], where the first framework was extended to include nonlinear kinematic hardening of Ohno-Wang type. When Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{L}}_{\mathrm{ke}}}
is symmetric, the same arguments as before relating to Equations 29 and 30 hold true. Hence, for appropriate free energies the two frameworks give the same model also for nonlinear kinematic hardening.

5 CONCLUDING REMARKS

We have shown that the two different frameworks, introduced by [4] and [5], can give equivalent models for isotropic free energies. The major difference between the models used for the different frameworks is the kinematic free energy. To obtain the same model, the same structure of the kinematic free energy must be used, but with a different variable (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\bar{\boldsymbol{C}}}_{\mathrm{ke}}}

or Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar{\boldsymbol{c}}_{\mathrm{ke}}}

). The numerical results confirm these theoretical findings. They further show that the difference between the formulations of free energy has a negligible effect on the material response up to a stretch of 5 for uniaxial loading, and no effect during simple shear.


6 ACKNOWLEDGMENTS

This work has been partly financed within the European Horizon 2020 Joint Technology Initiative Shift2Rail through contract no. 730841. The use of AceGen [12] has been very effective in speeding up the implementation of the material models.

REFERENCES

[1] F. A. M Alwahdi, A. Kapoor, and F. J. Franklin. "Subsurface microstructural analysis and mechanical properties of pearlitic rail steels in service". Wear 302.1-2 (2013), pp. 1453-1460. doi: 10.1016/j.wear.2012.12.058.
[2] K. Cvetkovski and J. Ahlström. "Characterisation of plastic deformation and thermal softening of the surface layer of railway passenger wheel treads". Wear 300.1-2 (2013), pp. 200-204. doi: 10.1016/j.wear.2013.01.094.
[3] A. Lion. "Constitutive modelling in finite thermoviscoplasticity: a physical approachbased on nonlinear rheological models". International Journal of Plasticity 16.5 (2000), pp. 469-494. doi: 10.1016/S0749-6419(99)00038-8.
[4] W. Dettmer and S. Reese. "On the theoretical and numerical modelling of Armstrong-Frederick kinematic hardening in the finite strain regime". In: Computer Methods in Applied Mechanics and Engineering 193.1 (2004), pp. 87-116. doi: 10.1016/j.cma.2003.09.005.
[5] M. Wallin, M. Ristinmaa, and N. S. Ottosen. "Kinematic hardening in large strain plasticity". In: European Journal of Mechanics - A/Solids 22.3 (2003), pp. 341-356. doi: 10.1016/S0997-7538(03)00026-3.
[6] M. Wallin and M. Ristinmaa. "Deformation gradient based kinematic hardening model". In: International Journal of Plasticity 21.10 (2005), pp. 2025-2050. doi:10.1016/j.ijplas.2005.01.007.
[7] G. Johansson, J. Ahlström, and M. Ekh. "Parameter identification and modeling of large ratcheting strains in carbon steel". In: Computers and Structures 84.15-16 (2006), pp. 1002-1011. doi: 10.1016/j.compstruc.2006.02.016.
[8] N. Larijani, G. Johansson, and M. Ekh. "Hybrid micro-macromechanical modeling of anisotropy evolution in pearlitic steel". In: European Journal of Mechanics - A/Solids 38 (Mar. 2013), pp. 38-47. doi: 10.1016/j.euromechsol.2012.09.011.
[9] J. C. Simo. "A Framework For Finite Strain Based on Maximum plastic dissipation and the multiplicative decomposition: Part I. Continuum formulation". In: Computer Methods in Applied Mechanics and Engineering 66 (1988), pp. 199-219.
[10] A. Lion. I. N. Vladimirov, M. P. Pietryga, and S. Reese. "On the modelling of non-linear kinematic hardening at finite strains with application to springback - Comparison of time integration algorithms". In: International Journal for Numerical Methods in Engineering 75.1 (July 2008), pp. 1-28. doi: 10.1002/nme.2234.
[11] Y. Zhu et al. "An extended cyclic plasticity model at finite deformations describing the Bauschinger effect and ratchetting behavior". In: 13th International Conference on Fracture 2013, ICF 2013 5 (2013), pp. 1-11.
[12] J. Korelc. "Multi-language and Multi-environment Generation of Nonlinear Finite Element Codes". In: Engineering with Computers 18.4 (Nov. 2002), pp. 312-327. doi: 10.1007/s003660200028.
Back to Top

Document information

Published on 20/05/17
Submitted on 12/05/17

Licence: Other

Document Score

0

Views 113
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?