Elastiseen muodonmuutokseen liittyvät aliohjelmat (uelastic.f90)

Tarjolla on aliohjelmat elastisen jäykkyysmatriisin laskemiseen ja jännityksen päivittämiseen elastisissa ja hyperelastisissa ongelmissa. Käyttäjän määrittämää elastista tai hyperelastista konstitutiivista relaatiota käytettäessä syötetiedostoon asetetaan ensin !ELASTIC, TYPE=USER tai !HYPERELASTIC, TYPE=USER ja annetaan tarvittavat materiaalivakiot, minkä jälkeen toteutetaan aliohjelmat uElasticMatrix ja uElasticUpdate.

(1) Elastisen jäykkyysmatriisin laskeva aliohjelma

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: Materiaalivakiot sisältävä taulukko (enintään 100)
• strain: Green–Lagrangen venymä
• D: Elastinen matriisi

(2) Jännityksen laskeva aliohjelma

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: Materiaalivakiot sisältävä taulukko (enintään 100)
• strain: Green–Lagrangen venymä
• stress: Jännitys