Skip to content

Subrutini zinazohusiana na mgeuko elastiki (uelastic.f90)

Hutoa subrutini za kukokotoa matriki ya ugumu elastiki na kusasisha mkazo kwa matatizo ya elastiki na hiperelastiki. Unapotumia mlinganyo wa konstitutivu wa elastiki au hiperelastiki iliyofafanuliwa na mtumiaji, kwanza weka !ELASTIC, TYPE=USER au !HYPERELASTIC, TYPE=USER katika faili ya ingizo na uingize konstanti za nyenzo zinazohitajika, kisha unda subrutini uElasticMatrix na uElasticUpdate.

(1) Subrutini 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: Safu inayohifadhi konstanti za nyenzo (kiwango cha juu 100)
  • strain: mnyumbuko wa Green-Lagrange
  • D: Matriki ya elastiki

(2) Subrutini 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: Safu inayohifadhi konstanti za nyenzo (kiwango cha juu 100)
  • strain: mnyumbuko wa Green-Lagrange
  • stress: Mkazo