Elastinės deformacijos paprogramės (uelastic.f90)¶

Pateikiamos elastinio standumo matricai apskaičiuoti ir įtempiams atnaujinti elastiniuose bei hiperelastiniuose uždaviniuose skirtos paprogramės. Norėdami naudoti naudotojo apibrėžtą elastinį ar hiperelastinį konstitucinį ryšį, pirmiausia įvesties faile nustatykite !ELASTIC, TYPE=USER arba !HYPERELASTIC, TYPE=USER ir įveskite reikiamas medžiagos konstantas, tada realizuokite paprogrames uElasticMatrix ir uElasticUpdate.

(1) Elastinio standumo matricos skaičiavimo paprogramė¶

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: Masyvas, kuriame saugomos medžiagos konstantos (ne daugiau kaip 100)
• strain: Greeno–Lagranžo deformacija
• D: Tamprumo matrica

(2) Įtempių skaičiavimo paprogramė¶

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: Masyvas, kuriame saugomos medžiagos konstantos (ne daugiau kaip 100)
• strain: Greeno–Lagranžo deformacija
• stress: Įtempiai