Lewati ke isi

Subrutin untuk Deformasi Elastik (uelastic.f90)

Disediakan subrutin untuk menghitung matriks kekakuan elastik dan memperbarui tegangan pada masalah elastik dan hiperelastik. Untuk menggunakan hubungan konstitutif elastik atau hiperelastik yang didefinisikan pengguna, pertama-tama atur !ELASTIC, TYPE=USER atau !HYPERELASTIC, TYPE=USER dalam file input dan masukkan konstanta material yang diperlukan, kemudian implementasikan subrutin uElasticMatrix dan uElasticUpdate.

(1) Subrutin untuk Menghitung Matriks Kekakuan Elastik

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 yang menyimpan konstanta material (maksimum 100)
  • strain: Regangan Green-Lagrange
  • D: Matriks elastik

(2) Subrutin untuk Menghitung Tegangan

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 yang menyimpan konstanta material (maksimum 100)
  • strain: Regangan Green-Lagrange
  • stress: Tegangan