Підпрограми для пружної деформації (uelastic.f90)

Надаються підпрограми для обчислення матриці пружної жорсткості та оновлення напружень у задачах пружності й гіперпружності. Щоб використати визначене користувачем пружне або гіперпружне визначальне співвідношення, спочатку задайте !ELASTIC, TYPE=USER або !HYPERELASTIC, TYPE=USER у вхідному файлі та введіть потрібні матеріальні константи, а потім реалізуйте підпрограми uElasticMatrix і uElasticUpdate.

(1) Підпрограма для обчислення матриці пружної жорсткості

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: Масив для зберігання матеріальних констант (максимум 100)
• strain: Деформація Гріна–Лагранжа
• D: Матриця пружності

(2) Підпрограма для обчислення напружень

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: Масив для зберігання матеріальних констант (максимум 100)
• strain: Деформація Гріна–Лагранжа
• stress: Напруження