This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
laws:cazacutn [2022/03/17 10:24] carlos |
laws:cazacutn [2022/10/14 15:52] (current) carlos |
||
---|---|---|---|
Line 31: | Line 31: | ||
Where: | Where: | ||
* $\bar{\Sigma}_{CPB06}$ is the CPB06 yield stress: | * $\bar{\Sigma}_{CPB06}$ is the CPB06 yield stress: | ||
- | \[\bar{\Sigma}_{CPB06}=\tilde{m}\Big[ \overset{3}{\underset{J=1}{\Sigma}} \big(|\Sigma_{J}| - k\Sigma_{J}\big)^a \Big] ^{\frac{1}{a}}\] | + | \[\bar{\Sigma}_{CPB06}=\tilde{m}\Big[ \overset{3}{\underset{i=1}{\Sigma}} \big(|\Sigma_{i}| - k\Sigma_{i}\big)^a \Big] ^{\frac{1}{a}}\] |
* $\sigma_{y}$ is the current yield stress of the material: | * $\sigma_{y}$ is the current yield stress of the material: | ||
\[\sigma_{y} = \sigma_{0} + S_{R}\big[1-exp\big(-C_{R}\bar{\epsilon}^{p}\big)\big]\] | \[\sigma_{y} = \sigma_{0} + S_{R}\big[1-exp\big(-C_{R}\bar{\epsilon}^{p}\big)\big]\] | ||
* STF is the stress transformation function, containing all the damage-related variables: | * STF is the stress transformation function, containing all the damage-related variables: | ||
- | \[STF=2fq_{1}cosh\Big[\frac{3q_{2}\big(\sigma_{m} - X_{m}\big)}{h\sigma_{y}}\Big] - q_{3}f^{2} - 1\] | + | \[STF= 1 - 2fq_{1}cosh\Big[\frac{3q_{2}\big(\sigma_{m} - X_{m}\big)}{h\sigma_{y}}\Big] - q_{3}f^{2}\] |
+ | The corrected stress ($\hat{\sigma}$) and backstress ($\hat{X}$) are respectively calculated as: | ||
+ | \[\hat{\sigma} = L_{ijmn}T_{mnkl}\sigma_{kl}\] | ||
+ | \[\hat{X}= L_{ijmn}T_{mnkl}X_{kl}\] | ||
+ | |||
+ | |||
+ | The backstress tensor is calculated using the Armstrong-Frederick model: | ||
+ | |||
+ | \[dX = S_{x}[C_{x}d{\epsilon^p}-Xd\bar{\epsilon}^p]\] | ||
==== Files ==== | ==== Files ==== | ||
Line 46: | Line 54: | ||
==== Subroutines ==== | ==== Subroutines ==== | ||
- | ^ Files ^ Contained subroutines ^ Description ^ | + | ^ Files ^ Contained subroutines ^ Description ^ |
- | | LCAZACUTN.F | LCAZACUTN |Main Prepro LCAZACUTN subroutine | | + | | LCAZACUTN.F | LCAZACUTN |Main Prepro LCAZACUTN subroutine | |
- | | CAZACUTN.F | CAZACUTN |Main Lagamine CAZACUTN subroutine | | + | | CAZACUTN.F | CAZACUTN |Main Lagamine CAZACUTN subroutine | |
- | | ::: | CAZACUTNFUN |Calculation of CAZACUTN yield locus | | + | | ::: | CAZACUTNFUN |Calculation of CAZACUTN yield locus | |
- | | THZCOAL.F | THOMASON_ZHANG |Calculation of coalescence criterion| | + | | COALCRITERIA.F | THOMASON_ZHANG |Calculation of [[https://www.sciencedirect.com/science/article/pii/S0013794400000552|ThZ coalescence criterion]]| |
+ | | ::: | LSP_2006 |Calculation of [[https://www.sciencedirect.com/science/article/pii/S0013794405003048|LSP06 coalescence criterion]]| | ||
===== Availability ===== | ===== Availability ===== | ||
Line 120: | Line 129: | ||
| |||| | | |||| | ||
^ |**CASE (3)**: //Growth, nucleation and coalescence are active// ||| | ^ |**CASE (3)**: //Growth, nucleation and coalescence are active// ||| | ||
- | | |^Line 8 (3G10.0) Nucleation model parameters ^^ | + | | |^Line 10 (3G10.0) Nucleation model parameters ^^ |
|::: ||$F_{N}$ | Total nucleated porosity ratio, //param(36, ilaw)// | | |::: ||$F_{N}$ | Total nucleated porosity ratio, //param(36, ilaw)// | | ||
|::: ||$S_{N}$ | Standard deviation, //param(37, ilaw)// | | |::: ||$S_{N}$ | Standard deviation, //param(37, ilaw)// | | ||
|::: ||$\epsilon_{N}$ | Standard mean, //param(38, ilaw)// | | |::: ||$\epsilon_{N}$ | Standard mean, //param(38, ilaw)// | | ||
- | |::: |^Line 9 (3G10.0) Coalescence model parameters ^^ | + | |::: |^Line 11 (3G10.0) Coalescence model parameters ^^ |
|::: ||$f_{U}$ | Ultimate porosity ratio, //param(39, ilaw)// | | |::: ||$f_{U}$ | Ultimate porosity ratio, //param(39, ilaw)// | | ||
|::: ||$f_{F}$ | Fracture porosity ratio, //param(40, ilaw)// | | |::: ||$f_{F}$ | Fracture porosity ratio, //param(40, ilaw)// | | ||
- | |::: ||$f_{cr}$ | Critical porosity ratio for onset of coalescence, //VARIN(19,ilaw), If 0, Thomason criterion is applied// | | + | |::: ||$f_{cr}$ | Critical porosity ratio for coalescence onset, //VARIN(19,ilaw) [If 0, Thomason criterion is applied]// | |
- | |::: |^Line 10 (3G10.0) Isotropic hardening law parameters ^^ | + | |
- | |::: ||$\sigma_{0}$ | Initial Yield stress [MPa], //param(30, ilaw)// | | + | |
- | |::: ||$S_{R}$ | Saturation rate [MPa], //param(31, ilaw)// | | + | |
- | |::: ||$C_{R}$ | Saturation value [-], //param(32, ilaw)// | | + | |
- | |::: |^Line 11 (2G10.0) Kinematic hardening parameters ^^ | + | |
- | |::: ||$C_{X}$ | Saturation rate [-], //param(33, ilaw)// | | + | |
- | |::: ||$S_{X}$ | Saturation value [MPa], //param(34, ilaw)// | | + | |
- | | |||| | + | |
+ | ... If the previous was the $n^{th}$ line... | ||
+ | ^ |IF **IDELEM = 1** THEN ||| | ||
+ | | |^Line $(n+1)^{th}$ (2G10.0) DELEM control parameters ^^ | ||
+ | | ||FDELEM | Porosity ratio at which element starts being deleted, //param(42, ilaw)// | | ||
+ | |::: ||TDELEM | Time for deleting the element (linear interpolation), //param(43, ilaw)// | | ||