Nënrutinat për deformimin elastik (uelastic.f90)¶

Ofrohen nënrutina për llogaritjen e matricës së ngurtësisë elastike dhe përditësimin e sforcimit në probleme elastike dhe hiperelastike. Kur përdoret një ligj konstitutiv elastik ose hiperelastik i përcaktuar nga përdoruesi, fillimisht duhet të vendoset !ELASTIC, TYPE=USER ose !HYPERELASTIC, TYPE=USER në skedarin e hyrjes dhe të jepen konstantet e nevojshme të materialit; më pas duhet të krijohen nënrutinat uElasticMatrix dhe uElasticUpdate.

(1) Nënrutina për llogaritjen e matricës së ngurtësisë elastike¶

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: vargu që ruan konstantet e materialit (maksimumi 100)
• strain: deformimi Green-Lagrange
• D: matrica elastike

(2) Nënrutina për llogaritjen e sforcimit¶

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: vargu që ruan konstantet e materialit (maksimumi 100)
• strain: deformimi Green-Lagrange
• stress: sforcimi