User Tools

Site Tools


laws:cazacutn

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
laws:cazacutn [2022/03/17 10:19]
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 102: Line 111:
 |$C_{X}$ | Saturation value [MPa], //param(31, ilaw)// ​          | |$C_{X}$ | Saturation value [MPa], //param(31, ilaw)// ​          |
 ^Line 9 (4G10.0) Standard initial damage control parameters ​    ^^ ^Line 9 (4G10.0) Standard initial damage control parameters ​    ^^
-|$f_{0}$| Initial porosity ratio, // VARIN(32,ilaw) //           | +|$f_{0}$| Initial porosity ratio, // VARIN(5,ilaw) //           | 
-|$q_{1}$| Tvergaard&​Needleman parameter, // param(33,ilaw) //    | +|$q_{1}$| Tvergaard&​Needleman parameter, // param(32,ilaw) //    | 
-|$q_{2}$| Tvergaard&​Needleman parameter, // param(34,ilaw) //    | +|$q_{2}$| Tvergaard&​Needleman parameter, // param(33,ilaw) //    | 
-|$q_{3}$| Tvergaard&​Needleman parameter, // param(35,ilaw) //    |+|$q_{3}$| Tvergaard&​Needleman parameter, // param(34,ilaw) //    |
  
  
 |**SELECT CASE (IDAMAGE)** ​                                             |||| |**SELECT CASE (IDAMAGE)** ​                                             ||||
 ^    |**CASE (0)**: //No damage increment is calculated// ​               ||| ^    |**CASE (0)**: //No damage increment is calculated// ​               |||
-       ​|^Line 8 (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 9 (2G10.0) Kinematic hardening parameters ​                ^^ +
-|:::     ​||$S_{X}$ | Saturation rate  [-], //param(33, ilaw)// ​            | +
-|:::     ​||$C_{X}$ | Saturation value [MPa], //param(34, ilaw)// ​          | +
-|:::                                                                    ​||||+
 ^    |**CASE (1)**: //Growth is the only active damage mechanism// ​      ||| ^    |**CASE (1)**: //Growth is the only active damage mechanism// ​      |||
-|        |^Line 8 (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 9 (2G10.0) Kinematic hardening parameters ​                ^^ 
-|:::     ​||$S_{X}$ | Saturation rate  [-], //param(33, ilaw)// ​            | 
-|:::     ​||$C_{X}$ | Saturation value [MPa], //param(34, ilaw)// ​          | 
 |                                                                       |||| |                                                                       ||||
-^    |**CASE (2)**: //Growth and nucleation of voids are considered//    ||| +^    |**CASE (2)**: //Growth and nucleation of voids are active//        ||| 
-|        |^Line ​(3G10.0) Nucleation model parameters ​            ​^^ +|        |^Line ​10 (3G10.0) Nucleation model parameters ​                  ​^^ 
-|:::     ​||$F_{N}$ | Total nucleated porosity ratio, //param(40, ilaw)// ​  | +|:::     ​||$F_{N}$ | Total nucleated porosity ratio, //param(36, ilaw)// ​  | 
-|:::     ​||$S_{N}$ | Standard deviation, //param(40, ilaw)// ​              | +|:::     ​||$S_{N}$ | Standard deviation, //param(37, ilaw)// ​              | 
-|:::     ​||$\epsilon_{N}$ | Standard mean, //param(40, ilaw)// ​            +|:::     ​||$\epsilon_{N}$ | Standard mean, //param(38, ilaw)// ​            |
-|:::     ​|^Line 9 (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 10 (2G10.0) Kinematic hardening parameters ​               ^^ +
-|:::     ​||$S_{X}$ | Saturation rate  [-], //param(33, ilaw)// ​            | +
-|:::     ​||$C_{X}$ | Saturation value [MPa], //param(34, ilaw)// ​          |+
 |                                                                       ​|||| ​ |                                                                       ​|||| ​
 ^    |**CASE (3)**: //Growth, nucleation and coalescence are active// ​   ||| ^    |**CASE (3)**: //Growth, nucleation and coalescence are active// ​   |||
-|        |^Line ​(3G10.0) Nucleation model parameters ​                   ^^  +|        |^Line ​10 (3G10.0) Nucleation model parameters ​                  ​^^  
-|:::     ​||$F_{N}$ | Total nucleated porosity ratio, //param(40, ilaw)// ​  | +|:::     ​||$F_{N}$ | Total nucleated porosity ratio, //param(36, ilaw)// ​  | 
-|:::     ​||$S_{N}$ | Standard deviation, //param(41, ilaw)// ​              | +|:::     ​||$S_{N}$ | Standard deviation, //param(37, ilaw)// ​              | 
-|:::     ​||$\epsilon_{N}$ | Standard mean, //param(42, ilaw)// ​            | +|:::     ​||$\epsilon_{N}$ | Standard mean, //param(38, ilaw)// ​            | 
-|:::     ​|^Line ​(3G10.0) Coalescence model parameters ​                  ​^^  +|:::     ​|^Line ​11 (3G10.0) Coalescence model parameters ​                 ^^  
-|:::     ​||$f_{U}$ ​  | Ultimate porosity ratio, //param(43, ilaw)// ​       | +|:::     ​||$f_{U}$ ​  | Ultimate porosity ratio, //param(39, ilaw)// ​       | 
-|:::     ​||$f_{F}$ ​  | Fracture porosity ratio, //param(44, ilaw)// ​       | +|:::     ​||$f_{F}$ ​  | Fracture porosity ratio, //param(40, ilaw)// ​       | 
-|:::     ​||$f_{cr}$ | Critical porosity ratio for onset of coalescence, //param(45, ilaw)//   | +|:::     ​||$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)// ​  |
  
laws/cazacutn.1647508750.txt.gz · Last modified: 2022/03/17 10:19 by carlos