laws:vmvp
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
laws:vmvp [2024/11/15 14:06] carlos [Numerical model] |
laws:vmvp [2024/11/29 15:53] (current) carlos [VMVP: formulation generalities] |
||
|---|---|---|---|
| Line 4: | Line 4: | ||
| 3D isotropic viscoplastic damage law that enables the modeling of non-classical creep responses via Graham-Walles and a modified activation function-Norton viscosity function. | 3D isotropic viscoplastic damage law that enables the modeling of non-classical creep responses via Graham-Walles and a modified activation function-Norton viscosity function. | ||
| - | ==== Numerical model ==== | + | === Overview === |
| This law was implemented in the context of C.Rojas-Ulloa's PhD. project (01/2021-12/2025) on the modeling of the non-classical long-term creep response of Incoloy 800H (see [[https://doi.org/10.1016/j.camwa.2023.12.002|(C.Rojas-Ulloa et al., 2024)]] for more details). This colaw is based in the work of Hélène Morch ([[laws:chab|CHAB]], a von-Mises yield function combined with a Norton-type viscosity function, Kachanov-Lemaitre creep-fatigue damage and high flexibility for introducing parameters as $f(T)$). **VMVP** incorporates two new viscoplastic functions (Graham-Walles & AFN), a simple //IfW// creep-fatiogue damage formulation, and new parameter interpolation methods available for material parameters. | This law was implemented in the context of C.Rojas-Ulloa's PhD. project (01/2021-12/2025) on the modeling of the non-classical long-term creep response of Incoloy 800H (see [[https://doi.org/10.1016/j.camwa.2023.12.002|(C.Rojas-Ulloa et al., 2024)]] for more details). This colaw is based in the work of Hélène Morch ([[laws:chab|CHAB]], a von-Mises yield function combined with a Norton-type viscosity function, Kachanov-Lemaitre creep-fatigue damage and high flexibility for introducing parameters as $f(T)$). **VMVP** incorporates two new viscoplastic functions (Graham-Walles & AFN), a simple //IfW// creep-fatiogue damage formulation, and new parameter interpolation methods available for material parameters. | ||
| - | === VMVP: formulation generalities=== | + | ==== VMVP: formulation generalities==== |
| The yield surface is defined by the von-Mises yield criterion: an isotropic $J_{2}$-type function of the form:\\ | The yield surface is defined by the von-Mises yield criterion: an isotropic $J_{2}$-type function of the form:\\ | ||
| Line 23: | Line 23: | ||
| $\Phi = J_{2}(\underline{\tilde{\sigma}}-\underline{\mathbb{X}})-\sigma_{y}-R \leq 0$\\ | $\Phi = J_{2}(\underline{\tilde{\sigma}}-\underline{\mathbb{X}})-\sigma_{y}-R \leq 0$\\ | ||
| - | where:\\ | + | >where: |
| - | * $\sigma_{y}$ is the yield stress of the material.\\ | + | >* $\sigma_{y}$ is the yield stress of the material. |
| - | * $R$ is the isotorpic hardening, calculated as the Voce iso. hardening law:\\ | + | >* $R$ is the isotorpic hardening, calculated as the Voce iso. hardening law: |
| - | $\dot{R} = -b\cdot (Q-R)\cdot \dot{p} \hspace{1cm} -> \hspace{1cm} R(p)=Q\cdot (1-\exp(-b\cdot p))$\\ where $p$ is the equivalent plastic strain $(-)$, and $\dot{p}$ is the equivalent plastic strain rate $(s^{-1})$\\ | + | >$\hspace{1cm} \dot{R} = -b\cdot (Q-R)\cdot \dot{p} \hspace{1cm} -> \hspace{1cm} R(p)=Q\cdot (1-\exp(-b\cdot p))$\\ where $p$ is the equivalent plastic strain $(-)$, and $\dot{p}$ is the equivalent plastic strain rate $(s^{-1})$\\ |
| Inheriting the classical Chaboche-type formulation, the backstress tensor $\underline{\mathbb{X}}$ is calculated as the sum of a total of $nAF$ user-defined Armstrong-Frederick (AF) backstress equations. \\ | Inheriting the classical Chaboche-type formulation, the backstress tensor $\underline{\mathbb{X}}$ is calculated as the sum of a total of $nAF$ user-defined Armstrong-Frederick (AF) backstress equations. \\ | ||
| Line 46: | Line 46: | ||
| >> This term is important in presence of thermal gradients $\dot{T}$, and if we introduce the parameter $C_i$ as function of temperature.\\ | >> This term is important in presence of thermal gradients $\dot{T}$, and if we introduce the parameter $C_i$ as function of temperature.\\ | ||
| - | === VMVP: viscoplasticity === | + | ==== VMVP: viscoplasticity ==== |
| As a viscoplastiv law, **VMVP** includes a total of 3 viscoplastic functions. | As a viscoplastiv law, **VMVP** includes a total of 3 viscoplastic functions. | ||
| > **option ivp=1: Norton law** | > **option ivp=1: Norton law** | ||
| Line 61: | Line 61: | ||
| > | > | ||
| > **option ivp=4: AFN** (Activation Function $\times$ Norton) | > **option ivp=4: AFN** (Activation Function $\times$ Norton) | ||
| - | >> This equation is new, and was formulated after a thorough analysis on the creep response of Incoloy 800H under lows-tress and high-temperature loadings. It is is intended to follow the classical Norton behavior, where the material reaches a steady-state creep rate $p_{\text{ss}}$. However, the phenomena such as solid-solution hardening or particle-strengthening may induce an initial creep hardening effect where the material can reach a creep rate lower that that of the steady state. To address this initial hardening, a phenomenological temperature and stress-dependent activation function has been added. The form of AFN is:\\ $\exp \Bigl[ -a(1.05 + p - 0.001b)^{-c} \Bigr] \times \langle \frac{\sigma_v}{K} \rangle ^{n}$ \\ Let us now walk you through the Activation Function particulars:\\ | + | >> This equation is new, and was formulated after a thorough analysis on the creep response of Incoloy 800H under lows-tress and high-temperature loadings. It is is intended to follow the classical Norton behavior, where the material reaches a steady-state creep rate $p_{\text{ss}}$. However, the phenomena such as solid-solution hardening or particle-strengthening may induce an initial creep hardening effect where the material can reach a creep rate lower that that of the steady state. To address this initial hardening, a phenomenological temperature and stress-dependent activation function has been added. The form of AFN is:\\ $\dot{p}= \exp \Bigl[ -a(1.05 + p - 0.001b)^{-c} \Bigr] \times \langle \frac{\sigma_v}{K} \rangle ^{n}$ \\ Let us now walk you through the Activation Function particulars:\\ |
| >>> __Activation Function__. The Activation Function is of the form: | >>> __Activation Function__. The Activation Function is of the form: | ||
| >>> $\exp \Bigl[ -a(1.05 + p - 0.001b)^{-c} \Bigr]$ | >>> $\exp \Bigl[ -a(1.05 + p - 0.001b)^{-c} \Bigr]$ | ||
| Line 72: | Line 72: | ||
| >>>> where $p_{21}, p_{22}, p_{23}$ are user-defined material parameters. | >>>> where $p_{21}, p_{22}, p_{23}$ are user-defined material parameters. | ||
| - | === VMVP: damage formulations === | + | ==== VMVP: damage formulations ==== |
| Damage computation is added to enable the prediction of the effects of microcracks and its evolution towards the final fracture of the material. To use this formulation, the user must first define the damage methodology (dam_type) to be used in the simulation. A brief explanation of the avalable possibilities is given below:\\ | Damage computation is added to enable the prediction of the effects of microcracks and its evolution towards the final fracture of the material. To use this formulation, the user must first define the damage methodology (dam_type) to be used in the simulation. A brief explanation of the avalable possibilities is given below:\\ | ||
| > **dam_type = 0** No damage | > **dam_type = 0** No damage | ||
| Line 86: | Line 86: | ||
| >> · $\underline{\tilde{\sigma}}_{k+1}=\underline{\sigma}_{k+1} \cdot (1-D_{k})^{-1}$ | >> · $\underline{\tilde{\sigma}}_{k+1}=\underline{\sigma}_{k+1} \cdot (1-D_{k})^{-1}$ | ||
| >> · $D_{k+1}=f \Bigl( \underline{\sigma}_{k+1}, D_{k},... \Bigr)$\\ | >> · $D_{k+1}=f \Bigl( \underline{\sigma}_{k+1}, D_{k},... \Bigr)$\\ | ||
| - | > **dam_type = 3** Fully-coupled damage. | + | > **dam_type = 2** Fully-coupled damage. |
| >In this final case, the damage equations are included within the big Jacobian unknown set of equations, and are solved together with all the system in an implicit manner. It is a time-consuming and computationally heavy approach, where: | >In this final case, the damage equations are included within the big Jacobian unknown set of equations, and are solved together with all the system in an implicit manner. It is a time-consuming and computationally heavy approach, where: | ||
| >> · $\underline{\tilde{\sigma}}_{k+1}=\underline{\sigma}_{k+1} \cdot (1-D_{k+1})^{-1}$ | >> · $\underline{\tilde{\sigma}}_{k+1}=\underline{\sigma}_{k+1} \cdot (1-D_{k+1})^{-1}$ | ||
| Line 114: | Line 114: | ||
| >> Similar to the previous law, Kachanov formulation is introduced as: | >> Similar to the previous law, Kachanov formulation is introduced as: | ||
| >> $\dot{D}_{c}=k_3 \Bigl[\frac{\mathcal{Y(k_{4}\cdot \underline{\sigma})}}{S_{c}} \Bigr]^{S_{ce}} \frac{1}{(1-D)^{k_k}} $ | >> $\dot{D}_{c}=k_3 \Bigl[\frac{\mathcal{Y(k_{4}\cdot \underline{\sigma})}}{S_{c}} \Bigr]^{S_{ce}} \frac{1}{(1-D)^{k_k}} $ | ||
| - | >> where $S_{c}$, $S_{ce}$, $k_3$, $k_4$, $k_k$ are material constants. In particular, $S_c$ is conceived as $f(T)$. | + | >> where $S_{c}$, $S_{ce}$, $k_3$, $k_4$, $k_k$ are material constants. In particular, $S_c$ is conceived as $f(T)$.\\ |
| + | > **idam=2**: IfW Creep-Fatigue damage formulation | ||
| + | >This creep-fatigue damage law was developed by IfW institute in Aachen University. It came to our knowledge after being used by N. K. Karthik in the context of his PhD. project on the life prediction of metallic components subjected to creep-fatigue loadings. It is inspired in the Graham-Walles formulation, and is of the form: | ||
| + | > $\dot{D}=K_{\text{D}} J_{2}\Bigl( \underline{\tilde{\sigma}} - \underline{\mathbb{X}} \Bigr) + |\dot{T}|K_{\text{Td}}p^{m_{\text{Td}}}$ | ||
| + | > where $K_{\text{Td}}$, $K_{\text{Td}}$ and $m_{\text{Td}}$ are material parameters. | ||
| ===== Input file ===== | ===== Input file ===== | ||
| ==== Parameters defining the type of constitutive law ==== | ==== Parameters defining the type of constitutive law ==== | ||
| Line 120: | Line 124: | ||
| ^Line 01; format (2I5, 60A1)^^ | ^Line 01; format (2I5, 60A1)^^ | ||
| |IL|Law number| | |IL|Law number| | ||
| - | |ITYPE| 271| | + | |ITYPE| 282| |
| |COMMENT| Any comment (up to 60 characters) that will be reproduced on the output listing| | |COMMENT| Any comment (up to 60 characters) that will be reproduced on the output listing| | ||
laws/vmvp.1731675963.txt.gz · Last modified: 2024/11/15 14:06 by carlos
Page Tools
