Subrutiner för elastisk deformation (uelastic.f90)

Subrutiner tillhandahålls för beräkning av den elastiska styvhetsmatrisen och uppdatering av spänningen i elastiska och hyperelastiska problem. För att använda en användardefinierad elastisk eller hyperelastisk konstitutiv relation anger du först !ELASTIC, TYPE=USER eller !HYPERELASTIC, TYPE=USER i indatafilen och matar in de materialkonstanter som krävs, och implementerar sedan subrutinerna uElasticMatrix och uElasticUpdate.

(1) Subrutin för beräkning av den elastiska styvhetsmatrisen

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 som lagrar materialkonstanter (högst 100)
• strain: Green–Lagrange-töjning
• D: Elastisk matris

(2) Subrutin för beräkning av spänning

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 som lagrar materialkonstanter (högst 100)
• strain: Green–Lagrange-töjning
• stress: Spänning