Difference between revisions of "Steinberg-Lund Hardening"
(18 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( | | ||
<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, σ<sub>0</sub> is the initial yield stress, β and n are hardening law properties, α 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, σ<sub>0</sub> 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 G<sub>0</sub> is the initial shear modulus. The strain rate and temperature dependent term is defined in inverse form: | ||
<math>{d\alpha(Y_T,T)\over dt} = \left[ | | ||
<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> | |||
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 | The shear modulus temperature and pressure dependence is given by the same function as in the [[Steinberg-Cochran-Guinan Hardening]] law: | ||
| |||
<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, "A constitutive model for strain rates from 10-4 to 106," 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, "A constitutive model for strain rates from 10-4 to 106," 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 22: | Line 31: | ||
! Property !! Description | ! Property !! Description | ||
|- | |- | ||
| yield || Initial yield stress (σ<sub>0</sub> at zero pressure and the reference temperature). Enter in units | | yield || Initial yield stress (σ<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 β. It is dimensionless. | | betahard || Yield stress hardening term β. 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 | | 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 | | 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 | | 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 | | 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 | | 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]]. | ||
|- | |- | ||
| | | 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 | | 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 | 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 49: | Line 58: | ||
This [[Hardening Laws|hardening law]] defines three history variables, wihich are the cumulative equivalent plastic strain (absolute) defined as: | This [[Hardening Laws|hardening law]] defines three history variables, wihich are the cumulative equivalent plastic strain (absolute) defined as: | ||
<math>\alpha = \sqrt{2\over3}\ ||d\varepsilon_p||</math> | | ||
<math>\alpha = \sum \sqrt{2\over3}\ ||d\varepsilon_p||</math> | |||
the current rate- and temperature-dependent yield stress (Y<sub>T</sub>) in | where dε<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 [[ConsistentUnits Command#Legacy and Consistent Units|pressure units]], and the current plastic strain rate (dα/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
- ↑ 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).