Difference between revisions of "Steinberg-Lund Hardening"

From OSUPDOCS
Jump to navigation Jump to search
 
(10 intermediate revisions by the same user not shown)
Line 1: Line 1:
== Yield Stress ==
In the Steinberg-Lund [[Hardening Laws|hardening law]], the yield stress is given by
In the Steinberg-Lund [[Hardening Laws|hardening law]], the yield stress is given by


     
     
<math>\sigma_y = \left(Y_T(\dot\alpha,T)+\sigma_0\bigl(1 + \beta \alpha\bigr)^n\right){G(T,P)\over G_0}</math>
<math>\sigma_y = \left(\min\left\{\sigma_0\bigl(1 + \beta \alpha\bigr)^n,\sigma_y^{max}\right\}+\min\left\{Y_T(\dot\alpha,T),Y_P\right\}\right){G(T,P)\over G_0}</math>


where <math>Y_T(\dot\alpha,T)</math> is a strain rate and temperature dependent term, &sigma;<sub>0</sub> is the initial yield stress, &beta; and n are hardening law properties, &alpha; is the cumulative plastic strain, G(T,P) is the shear modulus (which may depend on temperature and pressure), and G<sub>0</sub> is the initial shear modulus. The strain rate and temperature dependent term is defined in inverse form:
where <math>Y_T(\dot\alpha,T)</math> is a strain rate and temperature dependent term, &sigma;<sub>0</sub> is the initial yield stress, &beta; and n are hardening law properties, &alpha; is the cumulative equivalent plastic strain, G(T,P) is the shear modulus (which may depend on temperature and pressure), and G<sub>0</sub> is the initial shear modulus. The strain rate and temperature dependent term is defined in inverse form:


&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<math>{d\alpha(Y_T,T)\over dt} = \left[{1\over C_1}\exp\left({2U_k\over kT}\left(1 - {Y_T\over Y_P}\right)^2\right)+{C_2\over Y_T} \right]^{-1}
<math>{d\alpha(Y_T,T)\over dt} = \left[{1\over C_1}\exp\left({2U_k\over kT}\left(1 - {Y_T\over Y_P}\right)^2\right)+{C_2\over Y_T} \right]^{-1}</math>
\qquad Y_T < Y_P</math>


where C<sub>1</sub>, C<sub>2</sub>, U<sub>k</sub>, and Y<sub>P</sub> are hardening law parameters.
where C<sub>1</sub>, C<sub>2</sub>, U<sub>k</sub>, and Y<sub>P</sub> are hardening law parameters.
The shear modulus temperature and pressure dependence are given by the same function as in the [[Steinberg-Cochran-Guinan Hardening]] law:
The shear modulus temperature and pressure dependence is given by the same function as in the [[Steinberg-Cochran-Guinan Hardening]] law:


&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<math>{G(T,P)\over G_0} = 1 + {G_P'\over G_0} P J^{1/3} + {G_T'\over G_0}(T-T_0)</math>
<math>{G(T,P)\over G_0} = 1 + {G_P'\over G_0} P J^{1/3} + {G_T'\over G_0}(T-T_0)</math>


where J is the relative volume change (V/V<sub>0</sub>), G<sub>P</sub>' and G<sub>T</sub>' are coefficients for pressure and temperature affects, T is current temperature, and T<sub>0</sub> is a reference temperature. For more details, see paper by Steinberg and Lund<ref>D. J. Steinberg and C. M. Lund, &quot;A constitutive model for strain rates from 10-4 to 106,&quot; J. Appl. Phys., 65, 1528-1533 (1989).</ref>.
where J is the relative volume change (V/V<sub>0</sub>), G<sub>P</sub>' and G<sub>T</sub>' are coefficients for pressure and temperature affects, T is current temperature, and T<sub>0</sub> is a reference temperature. For more details, see paper by Steinberg and Lund.<ref>D. J. Steinberg and C. M. Lund, &quot;A constitutive model for strain rates from 10-4 to 106,&quot; J. Appl. Phys., 65, 1528-1533 (1989).</ref>
 
=== Plasticity Theory ===
 
The implementation of this hardening law is formally non-associative J<sub>2</sub> plasticity theory as explained [[Steinberg-Cochran-Guinan Hardening#Plasticity Theory|here]].


== Hardening Law Properties ==
== Hardening Law Properties ==
Line 26: Line 31:
! Property !! Description  
! Property !! Description  
|-
|-
| yield ||  Initial yield stress (&sigma;<sub>0</sub> at zero pressure and the reference temperature). Enter in units of MPa.
| yield ||  Initial yield stress (&sigma;<sub>0</sub> at zero pressure and the reference temperature). Enter in [[ConsistentUnits Command#Legacy and Consistent Units|pressure units]].
|-
|-
| betahard || Yield stress hardening term &beta;. It is dimensionless.
| betahard || Yield stress hardening term &beta;. It is dimensionless.
|-
|-
| nhard || Exponent on cumulative plastic strain in hardening term. It is dimensionless.
| nhard || Exponent on cumulative equivalent plastic strain in hardening term. It is dimensionless.
|-
|-
| GPpG0 || The (G<sub>p</sub>'/G<sub>0</sub>) ratio term for pressure dependence of shear modulus. Enter in units MPa<sup>-1</sup>. Enter 0 to omit pressure dependence in shear modulus.
| GPpG0 || The (G<sub>p</sub>'/G<sub>0</sub>) ratio term for pressure dependence of shear modulus. Enter in [[ConsistentUnits Command#Legacy and Consistent Units|pressure<sup>-1</sup> units]]. Enter 0 to omit pressure dependence in shear modulus.
|-
|-
| GTpG0 || The (G<sub>T</sub>'/G<sub>0</sub>) ratio term for temperature dependence of shear modulus. Enter in units MPa<sup>-1</sup>. Enter 0 to omit temperature dependence in shear modulus.
| GTpG0 || The (G<sub>T</sub>'/G<sub>0</sub>) ratio term for temperature dependence of shear modulus. Enter in units K<sup>-1</sup>. Enter 0 to omit temperature dependence in shear modulus.
|-
|-
| C1SL || The C<sub>1</sub> constant in this law entered in units of 1/sec.
| C1SL || The C<sub>1</sub> constant in this law entered in [[ConsistentUnits Command#Legacy and Consistent Units|1/time units]].
|-
|-
| C2SL || The C<sub>2</sub> constant in this law entered in units of MPa-sec.
| C2SL || The C<sub>2</sub> constant in this law entered in [[ConsistentUnits Command#Legacy and Consistent Units|pressure-time units]].
|-
|-
| YP || The Peierls stress (Y<sub>P</sub>) and also the maximum rate-dependent yield stress. Enter in units of MPa.
| YP || The Peierls stress (Y<sub>P</sub>) and also the maximum rate-dependent yield stress. Enter in [[ConsistentUnits Command#Legacy and Consistent Units|pressure units]].
|-
|-
| Uk || An energy associated with forming kinks (U<sub>k</sub>). It has units eV.
| UkOverk || An energy associated with forming kinks (entered as U<sub>k</sub>/k). It has units degrees K.
|-
|-
| yieldMax || Maximum yield stress. Enter in units of MPa.
| yieldMax || Maximum yield stress. Enter in [[ConsistentUnits Command#Legacy and Consistent Units|pressure units]].
|}
|}


The reference temperature, T<sub>0</sub>, is set using the simulations stress free temperature and not in the [[Hardening Laws|hardening law]] properties.
The reference temperature, T<sub>0</sub>, is set using the simulation's [[Thermal Calculations#Stress Free Temperature|stress free temperature]] and not set using [[Hardening Laws|hardening law]] properties.


== History Data ==
== History Data ==
Line 57: Line 62:


where d&epsilon;<sub>p</sub> is the incremental plastic strain tensor in one time step,
where d&epsilon;<sub>p</sub> is the incremental plastic strain tensor in one time step,
the current rate- and temperature-dependent yield stress (Y<sub>T</sub>) in MPa, and the current equivalent plastic strain rate (d&alpha;/dt in 1/sec). These variables are stored as as history variables #1, #2, and #3.
the current rate- and temperature-dependent yield stress (Y<sub>T</sub>) in [[ConsistentUnits Command#Legacy and Consistent Units|pressure units]], and the current plastic strain rate (d&alpha;/dt in [[ConsistentUnits Command#Legacy and Consistent Units|1/time units]]). These variables are stored as history variables #1, #2, and #3.


== References ==


<references/>
<references/>

Latest revision as of 22:57, 29 October 2019

Yield Stress

In the Steinberg-Lund hardening law, the yield stress is given by

      [math]\displaystyle{ \sigma_y = \left(\min\left\{\sigma_0\bigl(1 + \beta \alpha\bigr)^n,\sigma_y^{max}\right\}+\min\left\{Y_T(\dot\alpha,T),Y_P\right\}\right){G(T,P)\over G_0} }[/math]

where [math]\displaystyle{ Y_T(\dot\alpha,T) }[/math] is a strain rate and temperature dependent term, σ0 is the initial yield stress, β and n are hardening law properties, α is the cumulative equivalent plastic strain, G(T,P) is the shear modulus (which may depend on temperature and pressure), and G0 is the initial shear modulus. The strain rate and temperature dependent term is defined in inverse form:

      [math]\displaystyle{ {d\alpha(Y_T,T)\over dt} = \left[{1\over C_1}\exp\left({2U_k\over kT}\left(1 - {Y_T\over Y_P}\right)^2\right)+{C_2\over Y_T} \right]^{-1} }[/math]

where C1, C2, Uk, and YP are hardening law parameters. The shear modulus temperature and pressure dependence is given by the same function as in the Steinberg-Cochran-Guinan Hardening law:

      [math]\displaystyle{ {G(T,P)\over G_0} = 1 + {G_P'\over G_0} P J^{1/3} + {G_T'\over G_0}(T-T_0) }[/math]

where J is the relative volume change (V/V0), GP' and GT' are coefficients for pressure and temperature affects, T is current temperature, and T0 is a reference temperature. For more details, see paper by Steinberg and Lund.[1]

Plasticity Theory

The implementation of this hardening law is formally non-associative J2 plasticity theory as explained here.

Hardening Law Properties

This hardening law can set the following properties:

Property Description
yield Initial yield stress (σ0 at zero pressure and the reference temperature). Enter in pressure units.
betahard Yield stress hardening term β. It is dimensionless.
nhard Exponent on cumulative equivalent plastic strain in hardening term. It is dimensionless.
GPpG0 The (Gp'/G0) ratio term for pressure dependence of shear modulus. Enter in pressure-1 units. Enter 0 to omit pressure dependence in shear modulus.
GTpG0 The (GT'/G0) ratio term for temperature dependence of shear modulus. Enter in units K-1. Enter 0 to omit temperature dependence in shear modulus.
C1SL The C1 constant in this law entered in 1/time units.
C2SL The C2 constant in this law entered in pressure-time units.
YP The Peierls stress (YP) and also the maximum rate-dependent yield stress. Enter in pressure units.
UkOverk An energy associated with forming kinks (entered as Uk/k). It has units degrees K.
yieldMax Maximum yield stress. Enter in pressure units.

The reference temperature, T0, is set using the simulation's stress free temperature and not set using hardening law properties.

History Data

This hardening law defines three history variables, wihich are the cumulative equivalent plastic strain (absolute) defined as:

      [math]\displaystyle{ \alpha = \sum \sqrt{2\over3}\ ||d\varepsilon_p|| }[/math]

where dεp is the incremental plastic strain tensor in one time step, the current rate- and temperature-dependent yield stress (YT) in pressure units, and the current plastic strain rate (dα/dt in 1/time units). These variables are stored as history variables #1, #2, and #3.

References

  1. D. J. Steinberg and C. M. Lund, "A constitutive model for strain rates from 10-4 to 106," J. Appl. Phys., 65, 1528-1533 (1989).