Subroutines voor elastische vervorming (uelastic.f90)

Er zijn subroutines beschikbaar voor het berekenen van de elastische stijfheidsmatrix en het bijwerken van spanning in elastische en hyperelastische problemen. Om een door de gebruiker gedefinieerde elastische of hyperelastische constitutieve relatie te gebruiken, stelt u eerst !ELASTIC, TYPE=USER of !HYPERELASTIC, TYPE=USER in het invoerbestand in en voert u de vereiste materiaalconstanten in; vervolgens implementeert u de subroutines uElasticMatrix en uElasticUpdate.

(1) Subroutine voor het berekenen van de elastische stijfheidsmatrix

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: Array waarin materiaalconstanten worden opgeslagen (maximaal 100)
• strain: Green-Lagrange-rek
• D: Elastische matrix

(2) Subroutine voor het berekenen van spanning

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: Array waarin materiaalconstanten worden opgeslagen (maximaal 100)
• strain: Green-Lagrange-rek
• stress: Spanning