ຊັບຣູທີນທີ່ກ່ຽວກັບການປ່ຽນຮູບແບບຍືດຫຍຸ່ນ (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: ຄວາມເຄັ່ງຮູບ Green-Lagrange
• 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: ຄວາມເຄັ່ງຮູບ Green-Lagrange
• stress: ຄວາມເຄັ້ນ