Processing math: 100%

Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools

You are here: metafor.ltas.ulg.ac.be » doc » Documentation » elements » General Points » Continuous isotropic damage
Trace:
doc:user:elements:volumes:continuousdamage

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:continuousdamage [2014/06/30 12:32] – [AlvBoneRemodContinousDamage2constant] bomandoc:user:elements:volumes:continuousdamage [2021/04/09 11:35] (current) – [LinGeersContinuousDamage] tanaka
Line 1: Line 1:
-====== Endommagement continu Isotrope ====== +====== Continuous isotropic damage ======
- +
-La classe ''ContinousDamage'' gère les différentes lois d'endommagement continu. Lorsqu'on veut implémenter une nouvelle loi d'endommagement, il faut définir l'évolution de la variable d'endommagement ΔD ainsi que ses dérivées par rapport à la pression, la déformation plastique et l'endommagement +
- +
-Lois implémentées dans Metafor +
  
 +The ''ContinousDamage'' class manages all continuous damage evolution laws. When a new law is defined, the evolution of the damage variable δD must be defined, and so must be its derivatives with respect to pressure, plastic strain and damage.
  
 +Laws implemented in Metafor
  
 ===== LemaitreChabocheContinuousDamage ===== ===== LemaitreChabocheContinuousDamage =====
Line 12: Line 9:
 === Description === === Description ===
  
 +Lemaitre & Chaboche damage model [[doc:user:elements:volumes:continuousdamage#References|[1,2]]].
 $$ $$
-\dot D = \left(\dfrac{\bar \sigma^2 R_\nu}{2ES\left(1-D\right)^2}\right)^s \dot \varepsilon^{pl} \mbox{ si } \varepsilon^{pl} > \varepsilon^{pl}_D+\dot D = \left(\dfrac{\bar \sigma^2 R_\nu}{2ES\left(1-D\right)^2}\right)^s \dot \varepsilon^{pl} \mbox{, if } \varepsilon^{pl} > \varepsilon^{pl}_D \mbox{, and } \eta > \eta_D
 $$ $$
  
-où la fonction de triaxialité est définie par:+where the triaxiliaty function is defined as:
  
 $$ $$
Line 22: Line 20:
 $$ $$
  
-avec $ E $ le module de Young du matériau, ν le coefficient du matériau, p la pression et ˉσ est la contrainte équivalente de Von-Mises.+where $ p $ is the pressure, ˉσ is Von Mises stress and η is the stress triaxiality ratio.
  
-=== Paramètres === +=== Parameters ===
  
-^          Nom                                                    ^        Codes Metafor         ^ Type de dépendance ^  +^   Name                                            ^     Metafor Code     ^  Dependency ^  
-| Module de Young E                              |      ''LEMAITRE_YOUNG''      |     ''TO/TM''      | +| Young Modulus E                                |      ''LEMAITRE_YOUNG''      |     ''TO/TM''      | 
-| Coefficient de Poisson ν                       |       ''LEMAITRE_NU''        |     ''TO/TM''      | +| Poisson Ratio ν                                |       ''LEMAITRE_NU''        |     ''TO/TM''      | 
-| Exposant s                                     |     ''LEMAITRE_SMALL_S''     |     ''TO/TM''      |+| Exponent s                                     |     ''LEMAITRE_SMALL_S''     |     ''TO/TM''      |
 | Coefficient S                                  |      ''LEMAITRE_BIG_S''      |     ''TO/TM''      | | Coefficient S                                  |      ''LEMAITRE_BIG_S''      |     ''TO/TM''      |
-| Déformation plastique seuil εplD |  ''LEMAITRE_EPL_THRESHOLD''  |     ''TO/TM''      |+| Plastic strain threshold εplD    |  ''LEMAITRE_EPL_THRESHOLD''  |     ''TO/TM''      | 
 +| Triaxiality threshold ηD    |  ''LEMAITRE_TRIAX_THRESHOLD''  |     ''TO/TM''      | 
 +===== BoneRemodContinousDamage =====
  
 +This law is designed for bone remodeling (extracted from Doblaré's law, which he uses only in elasticity). Damage evolution depends mostly on damage, surface available for remodeling and a "remodelling rate" function, which depends on stress state.
  
- 
- 
- 
- 
- 
- 
- 
- 
- 
- 
- 
- 
- 
- 
- 
- 
-===== BoneRemodContinousDamage ===== 
- 
-Il s'agit d'une loi d'endommagement pour le remodelage osseux (tirée de la loi d'endo. de Doblaré qu'il utilise uniquement en élasticité). La variation d'endo dépend principalement de l'endo, de la surface disponible pour le remodelage et d'une fonction "taux de remodelage", fonction elle même de l'état de contrainte. 
  
 === Description === === Description ===
Line 59: Line 41:
 $$ $$
  
-où+where
  
-Sv(d) est la surface disponible par unité de volume pour du remodelage (polynôme 5ème puissance de l'endo)+Sv(d) is the surface per unit volume available for remodeling (polynomial of degree 5 in d)
  
-et où+and where
  
 $$ $$
 \begin{align*} \begin{align*}
-\dot r &= \ cf_1(d, \rho_0)g_f&\text{ si }g_f>0 \\ +\dot r &= \ cf_1(d, \rho_0)g_f&\text{ if }g_f>0 \\ 
-\dot r &= -cf_1(d, \rho_0)g_r&\text{ si }g_r>0 +\dot r &= -cf_1(d, \rho_0)g_r&\text{ if }g_r>0 
 \end{align*} \end{align*}
 $$ $$
  
-avec +with
  
 $$ $$
Line 81: Line 63:
 $$ $$
  
-f,f1 sont fonctions de l'endo, +f,f1 are functions in the damage variable, u is a measure of the elastic strain energy.
-u est une mesure de l'énergie de défo élastique+
  
-=== Paramètres === +=== Parameters ===
  
-^          Nom                                            ^  Codes Metafor                     ^+^   Name                                            ^     Metafor Code     ^
 | Coefficient N                                       |        ''BONE_REMOD_N''            | | Coefficient N                                       |        ''BONE_REMOD_N''            |
-| Pourcentage de surface disponible k                 |       ''BONE_REMOD_K''             | +| Percentage of available surface k                   |       ''BONE_REMOD_K''             | 
-| Energie de déformation élastique de référence ψ  |     ''BONE_REMOD_PSI''             | +| Reference elastic strain energy ψ                |     ''BONE_REMOD_PSI''             | 
-| Demi-largeur de la zone morte   ω              |      ''BONE_REMOD_OMEGA''          |  +| Half width of the dead zone ω                  |      ''BONE_REMOD_OMEGA''          |  
-| Vitesse de remodelage c                             |  ''BONE_REMOD_C''                  | +| Remodeling speed c                                  |  ''BONE_REMOD_C''                  | 
-| Densité du matériau non endommagé ρ0[kg/m3]   |  ''BONE_REMOD_MASS_DENSITY''       |+| Density of undamaged materialρ0[kg/m3]        |  ''BONE_REMOD_MASS_DENSITY''       | 
 + 
 ==== AlvBoneRemodContinousDamage ==== ==== AlvBoneRemodContinousDamage ====
  
-Il s'agit d'une loi d'endommagement pour le remodelage de l'os alvéolaire. La variation d'endo dépend en plus de la pression.+This law is defined for the remodeling of the alveolar bone. Damage evolution also depends on pressure. 
  
 === Description === === Description ===
Line 102: Line 86:
 $$ $$
  
-où+where
  
-Sv(d) est la surface disponible par unité de volume pour du remodelage (polynôme 5ème puissance de l'endo)+Sv(d) is the surface per unit volume available for remodeling (polynomial of degree 5 in d)
  
 $$ $$
 \begin{align*} \begin{align*}
-\dot r &= cf_1(d, \rho_0)g_f &\text{ si }g_f>0 \text{ et } p>0 \\ +\dot r &= cf_1(d, \rho_0)g_f &\text{ if }g_f>0 \text{ and } p>0 \\ 
-\dot r &= -cf_1(d, \rho_0)g_f &\text{ si }g_f>0 \text{ et } p<0 \\ +\dot r &= -cf_1(d, \rho_0)g_f &\text{ if }g_f>0 \text{ and } p<0 \\ 
-\dot r &= -cf_1(d, \rho_0)g_r &\text{ si }g_r>0 +\dot r &= -cf_1(d, \rho_0)g_r &\text{ if }g_r>0 
 \end{align*} \end{align*}
 $$ $$
  
-avec +with
  
 $$ $$
Line 123: Line 107:
 $$ $$
  
 +f,f1 are functions in the damage variable, u is a measure of the elastic strain energy.
  
-$ f, f_1 $ sont fonctions de l'endo, $ u $ est une mesure de l'énergie de défo élastique.+=== Parameters === 
 +^   Name                                            ^     Metafor Code     ^  
 +| Coefficient $ N $                                       |        ''BONE_REMOD_N''            |    
 +| Percentage of available surface $ k $                   |       ''BONE_REMOD_K''            |    
 +| Reference elastic strain energy ψ                |     ''BONE_REMOD_PSI''          |   
 +| Remodeling speed c                                  |  ''BONE_REMOD_C''            | 
 +| Density of undamaged material ρ0[kg/m3]  |  ''BONE_REMOD_MASS_DENSITY'' |
  
-=== Paramètres ===  
- 
-^          Nom                                            ^        Codes Metafor                     ^    
-| Coefficient N                                       |        ''BONE_REMOD_N''            |    
-| Pourcentage de surface disponible k                 |       ''BONE_REMOD_K''            |    
-| Energie de déformation élastique de référence ψ  |     ''BONE_REMOD_PSI''          |   
-| Vitesse de remodelage c                             |  ''BONE_REMOD_C''            | 
-| Densité du matériau non endommagé ρ0[kg/m3]  |  ''BONE_REMOD_MASS_DENSITY'' | 
  
 ==== AlvBoneRemodContinousDamage2constant ==== ==== AlvBoneRemodContinousDamage2constant ====
  
-idem que précédente mais constantes de remodelage différentes en formation et résorption.+Same law than the previous one, except that remodeling constants are different in formation and resorption.
  
 === Description === === Description ===
Line 143: Line 126:
 $$ $$
 \begin{align*} \begin{align*}
-\dot r &= c_ff_1(d, \rho_0)g_f\;\;&\text{ si }g_f>0\;\;\text{ et } p>0 \\ +\dot r &= c_ff_1(d, \rho_0)g_f\;\;&\text{ if }g_f>0\;\;\text{ and } p>0 \\ 
-\dot r &= -c_rf_1(d, \rho_0)g_f\;\;&\text{ si }g_f>0\;\;\text{ et } p<0 \\ +\dot r &= -c_rf_1(d, \rho_0)g_f\;\;&\text{ if }g_f>0\;\;\text{ and } p<0 \\ 
-\dot r &= -c_rf_1(d, \rho_0)g_r\;\;&\text{ si }g_r>0 +\dot r &= -c_rf_1(d, \rho_0)g_r\;\;&\text{ if  }g_r>0 
 \end{align*} \end{align*}
 $$ $$
  
-avec+with
  
-<note>Manque certaines définitions</note>+<note> 
 +some definitions are lacking 
 +</note> 
 +...
  
-=== Paramètres ===  +=== Parameters === 
- +^   Name                                            ^     Metafor Code     ^ 
-^          Nom                                            ^        Codes Metafor              ^   +
 | Coefficient N                                       |        ''BONE_REMOD_N''            |    | Coefficient N                                       |        ''BONE_REMOD_N''            |   
-| Pourcentage de surface disponible k                 |       ''BONE_REMOD_K''            |    +| Percentage of available surface k                   |       ''BONE_REMOD_K''            |    
-| Energie de déformation élastique de référence ψ  |     ''BONE_REMOD_PSI''          |   +| Reference elastic strain energy ψ                |     ''BONE_REMOD_PSI''          |    
-| Vitesse de remodelage cf                             |  ''BONE_REMOD_CF''            | +| Remodeling speed cf                                |  ''BONE_REMOD_CF''            | 
-| Vitesse de remodelage cr                             |  ''BONE_REMOD_CR''            | +| Remodeling speed cr                                |  ''BONE_REMOD_CR''            | 
-| Densité du matériau non endommagé ρ0[kg/m3]            |  ''BONE_REMOD_MASS_DENSITY'' | +| Density of undamaged material $\rho_0 [\mbox{kg}/\mbox{m}^3] $  |  ''BONE_REMOD_MASS_DENSITY'' | 
- + 
 ===== LangsethContinousDamage ===== ===== LangsethContinousDamage =====
  
Line 169: Line 153:
  
 $$ $$
-\dot D = D_C\dfrac{\dot \varepsilon^{pl}}{\varepsilon^{pl}_f-\varepsilon^{pl}_D} \mbox{ si } \varepsilon^{pl} > \varepsilon^{pl}_D+\dot D = D_C\dfrac{\dot \varepsilon^{pl}}{\varepsilon^{pl}_f-\varepsilon^{pl}_D} \mbox{ if } \varepsilon^{pl} > \varepsilon^{pl}_D
 $$ $$
  
-où la déformation plastique de rupture est définie par:+where the plastic strain at failure is defined as:
  
 $$ $$
Line 178: Line 162:
 $$ $$
  
-où p est la pression et ˉσ est la contrainte équivalente de Von-Mises.+where p is the pressure and ˉσ the Von Mises stress.
  
-=== Paramètres === +=== Parameters ===
  
-^          Nom                                     ^        Codes Metafor         ^ Type de dépendance ^ +^   Name                                            ^     Metafor Code     ^  Dependency ^ 
 | D1                             |        ''LANGSETH_D1''       |     ''TO/TM''      | | D1                             |        ''LANGSETH_D1''       |     ''TO/TM''      |
 | D2                             |       ''LANGSETH_D2''        |     ''TO/TM''      | | D2                             |       ''LANGSETH_D2''        |     ''TO/TM''      |
Line 188: Line 172:
 | D4                             |      ''LANGSETH_D4''         |     ''TO/TM''      | | D4                             |      ''LANGSETH_D4''         |     ''TO/TM''      |
 | D5                             |      ''LANGSETH_D5''         |     ''TO/TM''      | | D5                             |      ''LANGSETH_D5''         |     ''TO/TM''      |
-| Endommagement DC               |      ''LANGSETH_DC''         |     ''TO/TM''      |+| Damage DC                      |      ''LANGSETH_DC''         |     ''TO/TM''      |
 | ˙εpl0         |      ''LANGSETH_EPSP0''      |     ''TO/TM''      | | ˙εpl0         |      ''LANGSETH_EPSP0''      |     ''TO/TM''      |
-| Température ambiante Troom   |      ''LANGSETH_ROOM''       |         -          | +| Room temperature Troom       |      ''LANGSETH_ROOM''       |         -          | 
-| Température de fusion Tmelt  |      ''LANGSETH_TMELT''      |         -          | +| Melting temperature Tmelt    |      ''LANGSETH_TMELT''      |         -          | 
-| Defo plastique seuil εplD |  ''LANGSETH_EPL_THRESHOLD''  |    -    | +| Plastic strain threshold εplD |  ''LANGSETH_EPL_THRESHOLD''  |    -    |
- +
- +
- +
  
  
 ===== GeersContinuousDamage ===== ===== GeersContinuousDamage =====
  
-Loi d'endommagement continu de Geers. On trouve plusieurs variantes de cette loi, toutes du même auteur, d'où le regroupement au sein d'une même classe mère. Il faut noter que, dans le modèle complet de Geers, l'endommagement est intégré de manière globale à toute la structure et non de manière locale (à chaque point d'intégration). [[mail:ppjeunechamps@ulg.ac.be|Je]] veux bien fournir des références aux personnes intéressées. Les modèles sont basés sur une grandeur caractéristique appelée κ. +Damage evolution law following Geers's models. Several laws actually exist, all of the same author, which is why they are gathered in a same class. If the full Geers's model, damage is integrated globally on the structure, and not locally at each integration point. [[mail:ppjeunechamps@ulg.ac.be|I]] can give references if needed. All models are based on a characteristic quantity, κ.
- +
-=== Paramètres communs à tous les modèles === +
- +
-^          Nom                                       ^    Codes Metafor    ^ Type de dépendance ^  +
-| Valeur de déclenchement κi  |  ''GEERS_KAPPA_I''  |     ''TO/TM''      | +
-| Valeur critique κc          |  ''GEERS_KAPPA_C''  |     ''TO/TM''      | +
- +
- +
- +
  
 +=== Parameters common to all models ===
  
 +^   Name                                 ^     Metafor Code     ^  Dependency ^ 
 +| Initiation value κi          |  ''GEERS_KAPPA_I''  |     ''TO/TM''      |
 +| Critical value   κc          |  ''GEERS_KAPPA_C''  |     ''TO/TM''      |
  
  
 ==== PowGeersContinuousDamage ==== ==== PowGeersContinuousDamage ====
  
-Loi puissance. La grandeur caractéristique κ est la déformation plastique équivalente εpl:+Power law. κ is the equivalent plastic strain εpl:
  
 $$ $$
-D = 1 - \left(\dfrac{\kappa_i}{\kappa}\right)^{n_1} \left(\dfrac{\kappa-\kappa_i}{\kappa_c-\kappa_i}\right)^{n_2} \mbox{ si } \kappa_i\leq\kappa\leq\kappa_c+D = 1 - \left(\dfrac{\kappa_i}{\kappa}\right)^{n_1} \left(\dfrac{\kappa-\kappa_i}{\kappa_c-\kappa_i}\right)^{n_2} \mbox{ if } \kappa_i\leq\kappa\leq\kappa_c
 $$ $$
  
-^          Nom          ^  Codes Metafor  ^ Type de dépendance ^ +^   Name   ^     Metafor Code     ^  Dependency ^ 
 | n1  |  ''GEERS_N1''   |     ''TO/TM''      | | n1  |  ''GEERS_N1''   |     ''TO/TM''      |
 | n2  |  ''GEERS_N2''   |     ''TO/TM''      | | n2  |  ''GEERS_N2''   |     ''TO/TM''      |
- 
- 
- 
- 
  
  
 ==== ExpGeersContinuousDamage ==== ==== ExpGeersContinuousDamage ====
  
-Loi exponentielle. La grandeur caractéristique κ est la déformation plastique équivalente ˉεpl+Exponential law. κ is the equivalent plastic strain ˉεpl
  
 $$ $$
Line 241: Line 211:
 $$ $$
  
-^          Nom            ^  Codes Metafor   ^ Type de dépendance ^ +^   Name     ^     Metafor Code     ^  Dependency ^ 
 | β  |  ''GEERS_BETA''  |     ''TO/TM''      | | β  |  ''GEERS_BETA''  |     ''TO/TM''      |
- 
  
  
Line 249: Line 218:
 ==== TanhGeersContinuousDamage ==== ==== TanhGeersContinuousDamage ====
  
-Loi en tangente hyperbolique. La grandeur caractéristique est la déformation plastique équivalente εpl+Hyperbolic tangent. κ is the equivalent plastic strain εpl
  
 $$ $$
Line 255: Line 224:
 $$ $$
  
-^          Nom                                       ^    Codes Metafor    ^ Type de dépendance ^  +^   Name                                            ^     Metafor Code     ^  Dependency ^  
-| Valeur de déclenchement κi  |  ''GEERS_KAPPA_I''  |     ''TO/TM''      | +| Initiation value κi  |  ''GEERS_KAPPA_I''  |     ''TO/TM''      | 
-| Valeur critique κc          |  ''GEERS_KAPPA_C''  |     ''TO/TM''      |+| Critical value κc          |  ''GEERS_KAPPA_C''  |     ''TO/TM''      |
  
  
Line 265: Line 234:
 ==== LinGeersContinuousDamage ==== ==== LinGeersContinuousDamage ====
  
-Loi linéaire: La grandeur caractéristique est une fonction de la triaxialité des contraintes et de la déformation plastique équivalente εpl+Law linear. κ is a function of the stress triaxiality and the equivalent plastic strain εpl
  
 $$ $$
 \dot{\kappa} = C\left<1+A\dfrac{p}{\bar\sigma}\right> \left(\varepsilon^{pl}\right)^B \dot\varepsilon^{pl} \dot{\kappa} = C\left<1+A\dfrac{p}{\bar\sigma}\right> \left(\varepsilon^{pl}\right)^B \dot\varepsilon^{pl}
 $$ $$
-où p est la pression et ¯σ est la contrainte équivalente de Von Mises. ⟨.⟩ sont les symboles de MacCaulay ( ⟨α⟩=α si α≥0 et 0 sinon)+where p is the pressure, and ¯σ the Von Mises stress. ⟨.⟩ are Macaulay symbols( ⟨α⟩=α if α≥0 and 0 otherwise)
  
 $$ $$
Line 276: Line 245:
 $$ $$
  
-^          Nom        ^  Codes Metafor  ^ Type de dépendance ^ +^   Name                                            ^     Metafor Code     ^  Dependency ^ 
 | A  |   ''GEERS_A''   |     ''TO/TM''      | | A  |   ''GEERS_A''   |     ''TO/TM''      |
 | B  |   ''GEERS_B''   |     ''TO/TM''      | | B  |   ''GEERS_B''   |     ''TO/TM''      |
 | C  |   ''GEERS_C''   |     ''TO/TM''      | | C  |   ''GEERS_C''   |     ''TO/TM''      |
  
 +===== References =====
 +
 +[1] Lemaitre J. A continuous damage mechanics model for ductile fracture. Journal of Engineering Materials and
 +Technology 1985;107:9–83.
 +
 +[2] Chaboche JL. Description thermodynamique et phénoménologique de la viscoélasticité cyclique avec endommagement.
 +PhD Thesis, Université Pierre et Marie Curie, Paris VI, 1978.
 +
 +[3]
 +
 +[4]
doc/user/elements/volumes/continuousdamage.1404124333.txt.gz · Last modified: 2016/03/30 15:22 (external edit)
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki