La classe Damage gère les différentes lois d'endommagement. Lorsqu'on veut implémenter une nouvelle loi d'endommagement, il faut définir plusieurs choses :

• La contrainte liée à l'endommagement à prendre en compte dans le critère de plasticité $ \sigma_{damage} $
• L'évolution de la variable d'endommagement $ D $
• La fonction d'adoucissement de la limite élastique $ \omega $

Lois implémentées dans Metafor

• La contrainte liée à l'endommagement $ \sigma_{damage} $ s'écrit:

$$ \sigma_{damage}=K\left[\dfrac{\dot{\bar{\varepsilon}}^{vp}\bar{\sigma}\left[1-D\right]}{\bar{\sigma}+p\dfrac{\partial f}{\partial p}}\right]^{m}\left[\bar{\varepsilon}^{vp}\right]^{n} $$

• La fonction d’adoucissement de la limite élastique $ \omega $ s'écrit:

$$ \omega \left(D,p,\sigma_{yield}\right) = \sqrt{1-2\dfrac{D}{D_{ult}}\cosh\left(\dfrac{3\alpha p}{\sigma_{yield}}\right)+\left[\dfrac{D}{D_{ult}}\right]^2} $$

• La loi d'évolution de la variable d'endommagement $ D $ s'écrit:

$$ \begin{align} \dot{D}&=\dfrac{D_{N}}{s_N\sqrt{2\pi}}\exp\left[-\dfrac{1}{2}\left[\dfrac{\bar{\varepsilon}^{vp}-\varepsilon_N}{s_N}\right]^2\right]\dot{\bar{\varepsilon}}^{vp}+\left[1-D\right]\text{tr}\left(\mathbf{D^{irr}}\right) &\text{ si }D<D_{crit} \notag \\ &=\dfrac{D_{ult}-D_{crit}}{\Delta \varepsilon}\dot{\bar{\varepsilon}}^{vp}& \text{ si } D>D_{crit} \notag \end{align} $$ === Paramètres === ^ Nom ^ Codes Metafor ^ Type de dépendance ^ | Viscosité ($ K $) | GURSON_K | - |

• La contrainte liée à l'endommagement $ \sigma_{damage} $ s'écrit:

$$ \sigma_{damage}=K\left(1-\sqrt{D}\right)\left[\dfrac{\dot{\bar{\varepsilon}}^{vp}\bar{\sigma}\left[1-D\right]}{\bar{\sigma}+p\dfrac{\partial f}{\partial p}}\right]^{m}\left[\bar{\varepsilon}^{vp}\right]^{n} $$

• La fonction d’adoucissement de la limite élastique $ \omega $ s'écrit:

$$ \omega \left(D,p,\sigma_{yield}\right) = \left(1-\sqrt{D}\right) - \dfrac{\alpha_\omega 3p}{\sigma_{yield}} $$

• La loi d'évolution de la variable d'endommagement $ D $ s'écrit:

$$ \begin{align} \dot{D}&=B\sigma_{vp}^v \left(\bar{\varepsilon}^{vp}\right)^b \dot{\bar{\varepsilon}}^{vp}+\left[1-D\right]E_v\eta\left(p\right)\text{tr}\left(\mathbf{D^{irr}}\right)&\text{ si }D<D_{crit} \notag \\ &=B\sigma_{vp}^v \left(\bar{\varepsilon}^{vp}\right)^b \dot{\bar{\varepsilon}}^{vp}+\left[1-D\right]F E_v\eta\left(p\right)\text{tr}\left(\mathbf{D^{irr}}\right)&\text{ si }D>D_{crit} \notag \end{align} $$ où $\eta\left(p\right)$ est donné par:

$$ \eta = \dfrac{3}{2} \dfrac{m+1}{m} \sinh\left(2\dfrac{2-m}{2+m}\dfrac{p}{\bar{\sigma}}\right) $$

• La contrainte liée à l'endommagement $ \sigma_{damage} $ s'écrit:

$$ \sigma_{damage}=K\left(1-\sqrt{D}\right)\left[\dfrac{\dot{\bar{\varepsilon}}^{vp}\bar{\sigma}\left[1-D\right]}{\bar{\sigma}+p\dfrac{\partial f}{\partial p}}\right]^{m}\left[\bar{\varepsilon}^{vp}\right]^{n} $$

• La fonction d’adoucissement de la limite élastique $ \omega $ s'écrit:

$$ \begin{align*} \omega \left(D,p,\sigma_{yield}\right) &= 1-\sqrt{D}\left(1+\dfrac{\alpha_\omega 3|p|}{\sigma_{yield}}\right) &\text{ si } |p| > \dfrac{p_{lim}}{PLIM} \\ &= \sqrt{\dfrac{3}{2}}\dfrac{\zeta+\sqrt{\beta^2-p^2}}{\sigma_{yield}} &\text{ si } |p| < \dfrac{p_{lim}}{PLIM} \end{align*} $$ où $$ \begin{eqnarray*} &p_{lim} &= \dfrac{1-\sqrt{D}}{\sqrt{D}} \dfrac{\sigma_{yield}}{3\alpha_\omega} \\ &\zeta &= \sqrt{\dfrac{2}{3}} \left(1 - \left(1+\dfrac{3\alpha p_{lim}}{\sigma_{yield}PLIM}\right) \sqrt{D}\right) \sigma_{yield} - \sqrt{\dfrac{3}{2}} \dfrac{p_{lim}}{3\alpha\sqrt{D}PLIM} \\ &\beta &= \sqrt{ \left(\dfrac{p_{lim}}{PLIM}\right)^2 + \dfrac{3}{2}\left(\dfrac{p_{lim}}{3\alpha\sqrt{D}PLIM}\right)^2} \end{eqnarray*} $$ * La loi d'évolution de la variable d'endommagement $ D $ s'écrit:

$$\begin{align} \dot{D}&=B\sigma_{vp}^v \left(\bar{\varepsilon}^{vp}\right)^b \dot{\bar{\varepsilon}}^{vp}+\left[1-D\right]E_v\eta\left(p\right)\text{tr}\left(\mathbf{D^{irr}}\right)&\text{ si }D<D_{crit} \notag \\ &=B\sigma_{vp}^v \left(\bar{\varepsilon}^{vp}\right)^b \dot{\bar{\varepsilon}}^{vp}+\left[1-D\right]F E_v\eta\left(p\right)\text{tr}\left(\mathbf{D^{irr}}\right)&\text{ si }D>D_{crit} \notag \end{align} $$ où $\eta\left(p\right)$ est donné par:

$$ \eta = \dfrac{3}{2} \dfrac{m+1}{m} \sinh\left(2\dfrac{2-m}{2+m}\dfrac{|p|}{\alpha_\eta\sigma_{yield}}\right) $$