Rugalmas deformációhoz tartozó alprogramok (uelastic.f90)

A rugalmas merevségi mátrix kiszámítására és a feszültség frissítésére szolgáló alprogramok állnak rendelkezésre rugalmas és hiperelasztikus feladatokhoz. Felhasználó által definiált rugalmas vagy hiperelasztikus konstitutív összefüggés használatához először állítsa be a bemeneti fájlban az !ELASTIC, TYPE=USER vagy az !HYPERELASTIC, TYPE=USER értéket, adja meg a szükséges anyagállandókat, majd valósítsa meg az uElasticMatrix és uElasticUpdate alprogramokat.

(1) A rugalmas merevségi mátrix kiszámítására szolgáló alprogram

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: Anyagállandókat tároló tömb (legfeljebb 100)
• strain: Green–Lagrange-alakváltozás
• D: Rugalmassági mátrix

(2) A feszültség kiszámítására szolgáló alprogram

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: Anyagállandókat tároló tömb (legfeljebb 100)
• strain: Green–Lagrange-alakváltozás
• stress: Feszültség