Podprogramy pro elastickou deformaci (uelastic.f90)

Jsou k dispozici podprogramy pro výpočet matice elastické tuhosti a aktualizaci napětí v elastických a hyperelastických úlohách. Chcete-li použít uživatelsky definovaný elastický nebo hyperelastický konstitutivní vztah, nejprve v souboru vstupních dat nastavte !ELASTIC, TYPE=USER nebo !HYPERELASTIC, TYPE=USER a zadejte požadované materiálové konstanty a poté implementujte podprogramy uElasticMatrix a uElasticUpdate.

(1) Podprogram pro výpočet matice elastické tuhosti

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: Pole uchovávající materiálové konstanty (maximálně 100)
• strain: Greenovo-Lagrangeovo přetvoření
• D: Elastická matice

(2) Podprogram pro výpočet napětí

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: Pole uchovávající materiálové konstanty (maximálně 100)
• strain: Greenovo-Lagrangeovo přetvoření
• stress: Napětí