Subroutine zinazohusiana na mgeuko elastiki (uelastic.f90)¶

Hutoa subroutine za kukokotoa matriki ya ugumu elastiki na kusasisha mkazo kwa matatizo ya elastiki na hyperelastic. Unapotumia constitutive equation ya elastiki au hyperelastic iliyofafanuliwa na mtumiaji, kwanza weka !ELASTIC, TYPE=USER au !HYPERELASTIC, TYPE=USER katika faili ya ingizo na uingize material constants zinazohitajika, kisha unda subroutine uElasticMatrix na uElasticUpdate.

(1) Subroutine ya kukokotoa matriki ya ugumu elastiki¶

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: Array inayohifadhi material constants (kiwango cha juu 100)
• strain: Green-Lagrange strain
• D: Matriki ya elastiki

(2) Subroutine ya kukokotoa mkazo¶

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: Array inayohifadhi material constants (kiwango cha juu 100)
• strain: Green-Lagrange strain
• stress: Mkazo