Differences

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

--- doc:user:elements:volumes:hyper_dev_potential [2025/11/14 14:56] – [HolzapfelGasserOgdenHyperPotential] vanhulle
+++ doc:user:elements:volumes:hyper_dev_potential [2025/11/14 15:33] (current) – [Anisotropic Elastic Potentials] vanhulle
@@ Line 78: / Line 78: @@
 W_{dev} = \sum^n_{i=1} W_{dev}^{(i)} \left(\bar{I}_1, \bar{I}_2, \bar{I}_3, \bar{I}_4^{(i)}, \bar{I}_5^{(i)} \right) = \sum^n_{i=1} W_{dev}^{(i)} \left(\bar{I}_1, \bar{I}_2, J, \bar{I}_4^{(i)}, \bar{I}_5^{(i)} \right)
 $$
+At the moment, a maximum of 3 separate principal directions can be given to the material law.
 === Reminders ===
@@ Line 122: / Line 123: @@
 | Holzapfel-Gasser-Ogden coefficient ($k_2$)  |  ''HYPER_HGO_K2''  |  ''TO/TM''  |
 | Fiber dispersion fraction ($d$)  |  ''HYPER_HGO_DISP''  |  ''TO/TM''  |
+| Direction of the first principal (fiber) direction ($a^1_x$) |  ''HYPER_FIB1_X''  |  -  |
+| Direction of the first principal (fiber) direction ($a^1_y$) |  ''HYPER_FIB1_Y''  |  -  |
+| Direction of the first principal (fiber) direction ($a^1_z$) |  ''HYPER_FIB1_Z''  |  -  |
+| Direction of the second principal (fiber) direction ($a^2_x$) |  ''HYPER_FIB2_X''  |  -  |
+| Direction of the second principal (fiber) direction ($a^2_y$) |  ''HYPER_FIB2_Y''  |  -  |
+| Direction of the second principal (fiber) direction ($a^2_z$) |  ''HYPER_FIB2_Z''  |  -  |
+| Direction of the third principal (fiber) direction ($a^3_x$) |  ''HYPER_FIB3_X''  |  -  |
+| Direction of the third principal (fiber) direction ($a^3_y$) |  ''HYPER_FIB3_Y''  |  -  |
+| Direction of the third principal (fiber) direction ($a^3_z$) |  ''HYPER_FIB3_Z''  |  -  |
@@ Line 130: / Line 140: @@
 W_{\text{BB},~dev}^{(i)}\left(\bar{I}_1, \bar{I}_4^{(i)}, \bar{I}_5^{(i)} \right) = \left[\alpha + \beta \left( \bar{I}_1-3 \right) + \gamma \left( \bar{I}^{(i)}_4 -1\right)\right]\left(\bar{I}^{(i)}_4 - 1\right) - \frac{1}{2}\alpha \left(\bar{I}^{(i)}_5 -1\right)
 $$
-or alternatively,
+where $\alpha$, $\beta$ and $\gamma$ are material parameters which are related to the engineering material constants from the fibers and matrix (see [[doc:user:elements:volumes:hyper_functionbased|Bonet-Burton material example]]). This model is actually directly derived from small-strain orthotropic (transversely isotropic) elasticity.
+=== Remarks ===
+Alternatively, another implementation of this material law is available where the hyperlastic law writes
+$$
+W_{\text{BB}}^{(i)}\left(J, \bar{I}_4^{(i)}, \bar{I}_5^{(i)} \right) = \left[\alpha + \beta~\text{ln}J + \gamma \left( \bar{I}^{(i)}_4 -1\right)\right]\left(\bar{I}^{(i)}_4 - 1\right) - \frac{1}{2}\alpha \left(\bar{I}^{(i)}_5 -1\right)
+$$
+by using the parameter ''HYPER_BB_USE_LNJ=true''.
+Note that in this case, $W_{\text{BB}}^{(i)}$ is not purely deviatoric since there is a coupling between $J$ and $\bar{I}_4^{(i)}$. Therefore, this formulation also contributes to the volumetric part of the deformation gradient.
+Mathematical derivations, such as the analytical tangent stiffness matrix, and information regarding the second form of the hyperelastic law can be found in {{ :doc:user:references:materials:vanhulle_251106_slides_function_based_hyper_v1.pdf |this presentation}}.
+=== Parameters ===
+^   Name                                                  ^  Metafor Code  ^ Dependency ^
+| Bonet-Burton coefficient ($\alpha$)  |  ''HYPER_BB_ALPHA''  |  ''TO/TM''  |
+| Bonet-Burton coefficient ($\beta$)  |  ''HYPER_BB_BETA''  |  ''TO/TM''  |
+| Bonet-Burton coefficient ($\gamma$)  |  ''HYPER_BB_GAMMA''  |  ''TO/TM''  |
+| Use the alternative Bonet-Burton law with $\beta~\text{ln}J$ \\ boolean: ''true'' (default) |  ''HYPER_BB_USE_LNJ''  |  ''TO/TM''  |
+| Direction of the first principal (fiber) direction ($a^1_x$) |  ''HYPER_FIB1_X''  |  -  |
+| Direction of the first principal (fiber) direction ($a^1_y$) |  ''HYPER_FIB1_Y''  |  -  |
+| Direction of the first principal (fiber) direction ($a^1_z$) |  ''HYPER_FIB1_Z''  |  -  |
+| Direction of the second principal (fiber) direction ($a^2_x$) |  ''HYPER_FIB2_X''  |  -  |
+| Direction of the second principal (fiber) direction ($a^2_y$) |  ''HYPER_FIB2_Y''  |  -  |
+| Direction of the second principal (fiber) direction ($a^2_z$) |  ''HYPER_FIB2_Z''  |  -  |
+| Direction of the third principal (fiber) direction ($a^3_x$) |  ''HYPER_FIB3_X''  |  -  |
+| Direction of the third principal (fiber) direction ($a^3_y$) |  ''HYPER_FIB3_Y''  |  -  |
+| Direction of the third principal (fiber) direction ($a^3_z$) |  ''HYPER_FIB3_Z''  |  -  |
+====== Rheological Laws ======
+{{:doc:user:ico-worker.png?25|Under construction}} Other laws will follow with the addition of visco-elasticity.
+===== CombinedElasticPotential =====
+=== Description ===
+The ''CombinedElasticPotential'' material law allows to combine two deviatoric hyperelastic potentials together as
 $$
-W_{\text{BB},~dev}^{(i)}\left(J, \bar{I}_4^{(i)}, \bar{I}_5^{(i)} \right) = \left[\alpha + \beta \text{ln}J + \gamma \left( \bar{I}^{(i)}_4 -1\right)\right]\left(\bar{I}^{(i)}_4 - 1\right) - \frac{1}{2}\alpha \left(\bar{I}^{(i)}_5 -1\right)
+\boldsymbol{\sigma}^e = \boldsymbol{\sigma}^e_1 + \boldsymbol{\sigma}^e_2
 $$
-where $\alpha$, $\beta$ and $\gamma$ are material parameters which are related to the engineering material constants from the fibers and matrix (see [[doc:user:elements:volumes:hyper_functionbased|Bonet-Burton material example]]).
+This can be illustrated using the following analogous rheological element
+{{ :doc:user:references:materials:rheoelast.png?300 |}}
-Mathematical derivations, such as the analytical tangent stiffness matrix, can be found in {{ :doc:user:references:materials:vanhulle_251106_slides_function_based_hyper_v1.pdf |this presentation}}.
+The main purpose of this element is to create anisotropic hyperelastic materials, as they are often composed of an isotropic (generally a Neo-Hookean) matrix component and an anisotropic fibrous component (see [[doc:user:elements:volumes:hyper_functionbased|anisotropic material examples]]). Nonetheless, this material law can also be used to add two or more deviatoric potentials, since ''CombinedElasticPotential'' can combine with itself.