Subrutine pentru deformație elastică (uelastic.f90)

Sunt furnizate subrutine pentru calculul matricei de rigiditate elastică și actualizarea tensiunii în probleme elastice și hiperelastice. Pentru a utiliza o relație constitutivă elastică sau hiperelastică definită de utilizator, setați mai întâi !ELASTIC, TYPE=USER sau !HYPERELASTIC, TYPE=USER în fișierul de intrare și introduceți constantele de material necesare, apoi implementați subrutinele uElasticMatrix și uElasticUpdate.

(1) Subrutină pentru calculul matricei de rigiditate elastică

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: Tabou care stochează constantele de material (maximum 100)
• strain: Deformație Green-Lagrange
• D: Matrice elastică

(2) Subrutină pentru calculul tensiunii

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: Tabou care stochează constantele de material (maximum 100)
• strain: Deformație Green-Lagrange
• stress: Tensiune