Subrutines per a la deformació elàstica (uelastic.f90)¶

Es proporcionen subrutines per calcular la matriu de rigidesa elàstica i actualitzar la tensió en problemes elàstics i hiperelàstics. Per utilitzar una relació constitutiva elàstica o hiperelàstica definida per l'usuari, primer establiu !ELASTIC, TYPE=USER o !HYPERELASTIC, TYPE=USER al fitxer d'entrada i introduïu les constants de material necessàries; a continuació, implementeu les subrutines uElasticMatrix i uElasticUpdate.

(1) Subrutina per calcular la matriu de rigidesa elàstica¶

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: Matriu que emmagatzema les constants del material (màxim 100)
• strain: Deformació de Green-Lagrange
• D: Matriu elàstica

(2) Subrutina per calcular la tensió¶

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: Matriu que emmagatzema les constants del material (màxim 100)
• strain: Deformació de Green-Lagrange
• stress: Tensió