Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools

You are here: metafor.ltas.ulg.ac.be » doc » Documentation » elements » General Points » Traditional Materials
Trace:
doc:user:elements:volumes:iso_hypo_materials

Differences

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

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
doc:user:elements:volumes:iso_hypo_materials [2015/07/14 14:40] – [KelvinVoigtViscoElastHypoMaterial] papeleuxdoc:user:elements:volumes:iso_hypo_materials [2022/07/14 12:41] (current) – papeleux
Line 17: Line 17:
 | Young's modulus                                    |  ''ELASTIC_MODULUS''  |              | | Young's modulus                                    |  ''ELASTIC_MODULUS''  |              |
 | Poisson Ratio                                      |   ''POISSON_RATIO''   |              | | Poisson Ratio                                      |   ''POISSON_RATIO''   |              |
 +| Material Stiffness  \\ (STIFF_ANALYTIC - STIFF_NUMERIC) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 | Objectivity Method \\ (Jaumann = 0, GreenNaghdi = 1) |    ''OBJECTIVITY''    |      -       | | Objectivity Method \\ (Jaumann = 0, GreenNaghdi = 1) |    ''OBJECTIVITY''    |      -       |
 | Orthotropic axis                                   |    ''ORTHO_AX1_X''    |      -       | | Orthotropic axis                                   |    ''ORTHO_AX1_X''    |      -       |
Line 37: Line 38:
 | Young's modulus                                |  ''ELASTIC_MODULUS''  |             | | Young's modulus                                |  ''ELASTIC_MODULUS''  |             |
 | Poisson ratio                        |   ''POISSON_RATIO''   |             | | Poisson ratio                        |   ''POISSON_RATIO''   |             |
 +| Material Stiffness  \\ (STIFF_ANALYTIC - STIFF_NUMERIC) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 | Objectivity method \\ (Jaumann = 0, GreenNaghdi = 1)   |    ''OBJECTIVITY''    |     -       | | Objectivity method \\ (Jaumann = 0, GreenNaghdi = 1)   |    ''OBJECTIVITY''    |     -       |
 | Orthotropic axis                              |    ''ORTHO_AX1_X''    |     -       | | Orthotropic axis                              |    ''ORTHO_AX1_X''    |     -       |
Line 64: Line 66:
 \begin{cases} \begin{cases}
 p^{1}  = p^{0} + 3K {\Delta\epsilon}_{ii} \\ p^{1}  = p^{0} + 3K {\Delta\epsilon}_{ii} \\
-s^{1}_{ij}  = s^{0}_{ij} + 2G {\Delta\epsilon}_{ij} + \eta \frac{{\Delta\epsilon}_{ij}}{\Delta t}+s^{1}_{ij}  = s^{0}_{ij} + 2G {\Delta\hat{\epsilon}}_{ij} + \eta \frac{{\Delta\hat{\epsilon}}_{ij}}{\Delta t}
 \end{cases} \end{cases}
 $$ $$
Line 78: Line 80:
 | Young's modulus                                   |  ''ELASTIC_MODULUS''  |    TO/TM     | | Young's modulus                                   |  ''ELASTIC_MODULUS''  |    TO/TM     |
 | Poisson Ratio                                     |   ''POISSON_RATIO''   |    TO/TM     | | Poisson Ratio                                     |   ''POISSON_RATIO''   |    TO/TM     |
 +| Material Stiffness  \\ (STIFF_ANALYTIC - STIFF_NUMERIC) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 | Viscosity coefficient                             |   ''VISCO_ETA''       |    TO/TM     | | Viscosity coefficient                             |   ''VISCO_ETA''       |    TO/TM     |
  
 +
 +===== GeneralizedMaxwellViscoElastHypoMaterial=====
 +
 +=== Description ===
 +
 +The GeneralizedMaxwell viscoelastic model result in adding of up to now, maximum 2 visco-elastic Maxwell branches to the elastic material. The similar spring/damper model is shown at the figure below
 +{{ :doc:user:elements:volumes:generalizedmaxwell.png?400 |}}
 +
 +Defining Maxwell viscous parameters from materials data (for each Maxwell branch): 
 +
 +$$
 +\begin{cases}
 +\Gamma_{i} &= \frac{\mu_{i}}  {2G}\\
 +\tau_{i}   &= \frac{\eta_{i}} {\mu_{i}}
 +\end{cases}
 +$$
 +
 +The stresses in each Maxwell branch consist in 2 effects : 
 +  - the relaxation of previous time step stresses in this Maxwell branch
 +  - the stress modification due to strain increment ${\Delta\epsilon}_{ij}$
 +
 +
 +$$
 +\begin{cases}
 +p^{1}      &= p^{0} + 3K {\Delta\epsilon}_{ii} \\
 +\underline{s}^{1}  &= \underline{s}_{E}^{1} + \underline{s}^{1}_{M1} + \underline{s}^{1}_{M2}
 +\end{cases}
 +$$
 +Stresses in each branch are computed using : 
 +$$
 +\begin{cases}
 +\underline{s}_{E}^{1}  &= \underline{s}_{E}^{0} + 2G {\Delta\hat{\underline{\epsilon}}} \\
 +\underline{s}^{1}_{M1} &= e^{(\frac{-\Delta t}{\tau_{1}})} \underline{s}^{0}_{M1} + \Gamma_{1} (1-e^{\frac{-\Delta t}{\tau_{1}}}) \frac{\tau_{1}}{\Delta t} 2G {\Delta\hat{\underline{\epsilon}}} \\
 +\underline{s}^{1}_{M2} &= e^{(\frac{-\Delta t}{\tau_{2}})} \underline{s}^{0}_{M2} + \Gamma_{2} (1-e^{\frac{-\Delta t}{\tau_{2}}}) \frac{\tau_{2}}{\Delta t} 2G {\Delta\hat{\underline{\epsilon}}} \\
 +\end{cases}
 +$$
 +
 +where $K$ and $G$ are compressibility and shear modulus, ${\Delta\underline{\epsilon}}$ is the stain increment during the time step  $\Delta t$.
 +these equations shown clearly the necessity to track history of the total stresses but also to each Maxwell branch stresses $\underline{s}^{1}_{M1}$ and $\underline{s}^{1}_{M2}$.
 +
 +
 +=== Parameters ===
 +       
 +
 +^   Name                                             ^     Metafor Code     ^  Dependency  ^
 +| Density                                           |   ''MASS_DENSITY''    |      -       |
 +| Young's modulus                                   |  ''ELASTIC_MODULUS''  |    TO/TM     |
 +| Poisson Ratio                                     |   ''POISSON_RATIO''   |    TO/TM     |
 +| Material Stiffness  \\ (STIFF_ANALYTIC - STIFF_NUMERIC) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 +| Number of Maxwell branches 1 (default) or 2       |  ''NBMAXWELLBRANCH''  |      -       |
 +| Maxwel 1 Stiffness                                |   ''VISCO_MU1''       |      -       |
 +| Maxwel 1 Viscosity                                |   ''VISCO_ETA1''      |      -       |
 +| Maxwel 2 Stiffness                                |   ''VISCO_MU2''       |      -       |
 +| Maxwel 2 Viscosity                                |   ''VISCO_ETA2''      |      -       |
 +
 +model implemented based on {{:doc:user:references:materials:1997_formulation_and_implementation_of_three-dimensional_viscoelasticity_at_small_and_finite_strains_kaliske_rothert.pdf|Kaliske M,
 +Rothert H. Formulation and implementation of three-dimensional viscoelasticity at small and finite strains. Computational Mechanics 1997;19:228-239.}}
 ===== EvpIsoHHypoMaterial ===== ===== EvpIsoHHypoMaterial =====
  
Line 102: Line 162:
 | Young's Modulus                                      |  ''ELASTIC_MODULUS''     |  ''TM''  | | Young's Modulus                                      |  ''ELASTIC_MODULUS''     |  ''TM''  |
 | Poisson Ratio                                     |   ''POISSON_RATIO''      |    -     | | Poisson Ratio                                     |   ''POISSON_RATIO''      |    -     |
 +| Material Stiffness  \\ (STIFF_ANALYTIC - STIFF_NUMERIC) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 | Thermal Expansion                                 |  ''THERM_EXPANSION''     |  ''TM''  | | Thermal Expansion                                 |  ''THERM_EXPANSION''     |  ''TM''  |
 | Number of the material law which defines the yield stress $\sigma_{yield}$    |      ''YIELD_NUM''       |    -     | | Number of the material law which defines the yield stress $\sigma_{yield}$    |      ''YIELD_NUM''       |    -     |
Line 147: Line 208:
 | Young's Modulus                                                              |  ''ELASTIC_MODULUS''  |  ''TO/TM''  | | Young's Modulus                                                              |  ''ELASTIC_MODULUS''  |  ''TO/TM''  |
 | Poisson Ratio                                                          |   ''POISSON_RATIO''   |  ''TO/TM''  | | Poisson Ratio                                                          |   ''POISSON_RATIO''   |  ''TO/TM''  |
 +| Material Stiffness  \\ (STIFF_ANALYTIC - STIFF_NUMERIC) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 | Thermal expansion                                                         |  ''THERM_EXPANSION''  |  ''TO/TM''  | | Thermal expansion                                                         |  ''THERM_EXPANSION''  |  ''TO/TM''  |
 | Number of the material law which defines the yield stress $\sigma_{yield}$    |      ''YIELD_NUM''    |    -        | | Number of the material law which defines the yield stress $\sigma_{yield}$    |      ''YIELD_NUM''    |    -        |
Line 180: Line 242:
 | Young's Modulus                                      |  ''ELASTIC_MODULUS''     |  ''TM''  | | Young's Modulus                                      |  ''ELASTIC_MODULUS''     |  ''TM''  |
 | Poisson Ratio                                  |   ''POISSON_RATIO''      |    -     | | Poisson Ratio                                  |   ''POISSON_RATIO''      |    -     |
 +| Material Stiffness  \\ (STIFF_ANALYTIC - STIFF_NUMERIC) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 | Thermal expansion                                 |  ''THERM_EXPANSION''     |  ''TM''  | | Thermal expansion                                 |  ''THERM_EXPANSION''     |  ''TM''  |
 | Number of the material law which defines the yield stress $\sigma_{yield}$  |      ''YIELD_NUM''       |    -     | | Number of the material law which defines the yield stress $\sigma_{yield}$  |      ''YIELD_NUM''       |    -     |
Line 213: Line 276:
 | Young's Modulus                                      |  ''ELASTIC_MODULUS''  |  ''TO/TM''  | | Young's Modulus                                      |  ''ELASTIC_MODULUS''  |  ''TO/TM''  |
 | Poisson Ratio                                  |   ''POISSON_RATIO''   |  ''TO/TM''  | | Poisson Ratio                                  |   ''POISSON_RATIO''   |  ''TO/TM''  |
 +| Material Stiffness  \\ (STIFF_ANALYTIC - STIFF_NUMERIC) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 | Thermal Expansion                                |  ''THERM_EXPANSION''  |  ''TO/TM''  | | Thermal Expansion                                |  ''THERM_EXPANSION''  |  ''TO/TM''  |
 | Number of the material law which defines the yield stress $\sigma_{yield}$     |    ''YIELD_NUM''  |  -   | | Number of the material law which defines the yield stress $\sigma_{yield}$     |    ''YIELD_NUM''  |  -   |
Line 252: Line 316:
 | Young's Modulus                           |   ''ELASTIC_MODULUS''  | | | Young's Modulus                           |   ''ELASTIC_MODULUS''  | |
 | Poisson Ratio                    |    ''POISSON_RATIO''   | | | Poisson Ratio                    |    ''POISSON_RATIO''   | |
 +| Material Stiffness  \\ (0 : Ana - 1 : Num) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 | Number of the material law which defines the yield stress $\sigma_{yield}$ |      ''YIELD_NUM''       |    -     | | Number of the material law which defines the yield stress $\sigma_{yield}$ |      ''YIELD_NUM''       |    -     |
 | Number of the damage evolution law              |     ''DAMAGE_NUM''     | | | Number of the damage evolution law              |     ''DAMAGE_NUM''     | |
Line 294: Line 359:
  
 $$ $$
-f=\dfrac{\overline{\sigma}}{1-D}-\sigma_{yield}=0+f=\dfrac{\overline{\sigma}}{1-h\cdot D}-\sigma_{yield}=0
 $$ $$
  
-where $ \overline{\sigma} $ is the equivalent stress computed as a function of [[doc:user:elements:volumes:plasticity_criterion#vonmisesplasticcriterion|Von Mises plasticy criterion]], $ \sigma_{yield} $ is the yield stress, $ D $ is the damage variable updated as a function of the [[doc:user:elements:volumes:continuousdamage|damage evolution law]].+where $ \overline{\sigma} $ is the equivalent stress computed as a function of [[doc:user:elements:volumes:plasticity_criterion#vonmisesplasticcriterion|Von Mises plasticy criterion]], $ \sigma_{yield} $ is the yield stress, $ D $ is the damage variable updated as a function of the [[doc:user:elements:volumes:continuousdamage|damage evolution law]]. Moreover, $h$ is the Micro-Crack Closure Effect parameter that makes the distinction of the weakening effect of damage under compressive and tensile stress states, which is defined as: 
 +$$ 
 + h = \left\{ 
 + \begin{array}{ll} 
 +  \text{DAMAGE_MCCE} &\mbox{, if } \dfrac{p}{J_2}< 0.0\\ 
 +  1.0 &\mbox{, if } \dfrac{p}{J_2}\geq 0.0\\ 
 + \end{array} 
 + \right. 
 +$$
  
 The evolution law coupled with plasticity can be integrated three ways depending on the parameter ''TYPE_INTEG'': The evolution law coupled with plasticity can be integrated three ways depending on the parameter ''TYPE_INTEG'':
Line 311: Line 384:
 | Young's Modulus                               |   ''ELASTIC_MODULUS''   |  ''TM''  | | Young's Modulus                               |   ''ELASTIC_MODULUS''   |  ''TM''  |
 | Poisson Ratio                           |    ''POISSON_RATIO''    |    -     | | Poisson Ratio                           |    ''POISSON_RATIO''    |    -     |
 +| Material Stiffness  \\ (0 : Ana - 1 : Num) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 | Number of the material law which defines the yield stress $\sigma_{yield}$ |      ''YIELD_NUM''       |    -     | | Number of the material law which defines the yield stress $\sigma_{yield}$ |      ''YIELD_NUM''       |    -     |
 | Number of the damage evolution law      |     ''DAMAGE_NUM''      |    -     | | Number of the damage evolution law      |     ''DAMAGE_NUM''      |    -     |
 | Initial damage                         |     ''DAMAGE_INIT''     |    -     | | Initial damage                         |     ''DAMAGE_INIT''     |    -     |
 | Integration method                         |     ''TYPE_INTEG''      |    -     | | Integration method                         |     ''TYPE_INTEG''      |    -     |
 +| Micro-Crack Closure Effect parameter (=1.0 by default)   |    ''DAMAGE_MCCE''     | |
  
 ===== ContinuousAnisoDamageEvpIsoHHypoMaterial ===== ===== ContinuousAnisoDamageEvpIsoHHypoMaterial =====
Line 345: Line 420:
 | Young's Modulus                               |   ''ELASTIC_MODULUS''   |  ''TM''  | | Young's Modulus                               |   ''ELASTIC_MODULUS''   |  ''TM''  |
 | Poisson Ratio                           |    ''POISSON_RATIO''    |    -     | | Poisson Ratio                           |    ''POISSON_RATIO''    |    -     |
 +| Material Stiffness  \\ (0 : Ana - 1 : Num) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 | Number of the material law which defines the yield stress $\sigma_{yield}$   |      ''YIELD_NUM''       |    -     | | Number of the material law which defines the yield stress $\sigma_{yield}$   |      ''YIELD_NUM''       |    -     |
 | Number of the damage evolution law      |     ''DAMAGE_NUM''      |    -     | | Number of the damage evolution law      |     ''DAMAGE_NUM''      |    -     |
Line 371: Line 447:
 | Young's Modulus                                                   |  ''ELASTIC_MODULUS''  |  ''TO/TM''  | | Young's Modulus                                                   |  ''ELASTIC_MODULUS''  |  ''TO/TM''  |
 | Poisson Ratio                                               |    ''POISSON_RATIO''  |    -       | | Poisson Ratio                                               |    ''POISSON_RATIO''  |    -       |
 +| Material Stiffness  \\ (0 : Ana - 1 : Num) \\ only if element Stiffness == STIFF_ANALYTIC | ''MATERIALSTIFFMETHOD''  |      -       |
 | Number of the material law which defines the yield stress $\sigma_{yield}$   |      ''YIELD_NUM''       |    -     | | Number of the material law which defines the yield stress $\sigma_{yield}$   |      ''YIELD_NUM''       |    -     |
 | Number of the damage evolution law                          |     ''DAMAGE_NUM''    |    -       | | Number of the damage evolution law                          |     ''DAMAGE_NUM''    |    -       |
doc/user/elements/volumes/iso_hypo_materials.txt · Last modified: 2022/07/14 12:41 by papeleux
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki