Aller au contenu

Sous-routines relatives à la déformation élastique (uelastic.f90)

Des sous-routines sont fournies pour calculer la matrice de rigidité élastique et mettre à jour les contraintes dans les problèmes élastiques et hyperélastiques. Pour utiliser une loi constitutive élastique ou hyperélastique définie par l'utilisateur, il faut d'abord définir !ELASTIC, TYPE=USER ou !HYPERELASTIC, TYPE=USER dans le fichier d'entrée et saisir les constantes de matériau nécessaires, puis créer les sous-routines uElasticMatrix et uElasticUpdate.

(1) Sous-routine de calcul de la matrice de rigidité élastique

subroutine uElasticMatrix( matl, strain, D )
    REAL(KIND=kreal), INTENT(IN) :: matl(:)
    REAL(KIND=kreal), INTENT(IN) :: strain(6)
    REAL(KIND=kreal), INTENT(OUT) :: D(6,6)
  • matl: Tableau stockant les constantes de matériau (maximum 100)
  • strain: Déformation de Green-Lagrange
  • D: Matrice élastique

(2) Sous-routine de calcul des contraintes

subroutine uElasticUpdate ( matl, strain, stress )
    REAL(KIND=kreal), INTENT(IN) :: matl(:)
    REAL(KIND=kreal), INTENT(IN) :: strain(6)
    REAL(KIND=kreal), INTENT(OUT) :: stress(6)
  • matl: Tableau stockant les constantes de matériau (maximum 100)
  • strain: Déformation de Green-Lagrange
  • stress: Contrainte