Skip to content

ຊັບຣູທີນສຳລັບວັດສະດຸທີ່ຜູ້ໃຊ້ກຳນົດ (umat.f)

ສະໜອງອິນເຕີເຟດສຳລັບການວິເຄາະການປ່ຽນຮູບຂອງວັດສະດຸທົ່ວໄປ ໂດຍບໍ່ຈຳກັດວ່າເປັນວັດສະດຸ ຍືດຫຍຸ່ນ, ໄຮເປີອີລາສຕິກ ຫຼື ອີລາສໂຕພລາສຕິກ. ເມື່ອໃຊ້ວັດສະດຸທີ່ຜູ້ໃຊ້ກຳນົດ, ກ່ອນອື່ນໃຫ້ກຳນົດ !USER_MATERIAL ໃນໄຟລ໌ຂໍ້ມູນເຂົ້າ ແລະໃສ່ຄ່າຄົງທີ່ວັດສະດຸທີ່ຈຳເປັນ, ຈາກນັ້ນສ້າງຊັບຣູທີນ uMatlMatrix ແລະ uUpdate.

(1) ຊັບຣູທີນສຳລັບຄຳນວນເມທຣິກຄວາມແຂງ

subroutine uMatlMatrix( mname, matl, strain, stress, fstat, D, dtime, ttime, temperature )
        character(len=*), intent(in)  :: mname
        real(kind=kreal), intent(in)  :: matl(:)
        real(kind=kreal), intent(in)  :: strain(6)
        real(kind=kreal), intent(in)  :: stress(6)
        real(kind=kreal), intent(in)  :: fstat(:)
        real(kind=kreal), intent(out) :: D(:,:)
        real(kind=kreal), intent(in)  :: dtime
        real(kind=kreal), intent(in)  :: ttime
        real(kind=kreal), optional    :: temperature
  • mname: ຊື່ວັດສະດຸ
  • matl: ອາເຣທີ່ເກັບຄ່າຄົງທີ່ວັດສະດຸ (ສູງສຸດ 100)
  • strain: ຄວາມເຄັ່ງຮູບ Green-Lagrange
  • stress: ຄວາມເຄັ້ນ 2nd Piola-Kirchhoff
  • fstat: ຕົວແປສະຖານະ
  • D: ຄວາມສຳພັນ ຄວາມສຳພັນວັດສະດຸ
  • dtime: ຂັ້ນເພີ່ມເວລາ
  • ttime: ເວລາລວມໃນຈຸດເລີ່ມຂອງ ຂັ້ນເພີ່ມເວລາ ປັດຈຸບັນ
  • temperature: ອຸນຫະພູມ

(2) ຊັບຣູທີນສຳລັບຄຳນວນອັບເດດຄວາມເຄັ່ງຮູບ ແລະຄວາມເຄັ້ນ

subroutine uUpdate(  mname, matl, strain, stress, fstat, dtime, ttime, temperature )
        character(len=\*), intent(in)    :: mname
        real(kind=kreal), intent(in)    :: matl(:)
        real(kind=kreal), intent(in)    :: strain(6)
        real(kind=kreal), intent(inout) :: stress(6)
        real(kind=kreal), intent(inout) :: fstat(:)
        real(kind=kreal), intent(in)    :: dtime
        real(kind=kreal), intent(in)    :: ttime
        real(kind=kreal), optional      :: temperature
  • mname: ຊື່ວັດສະດຸ
  • matl: ອາເຣທີ່ເກັບຄ່າຄົງທີ່ວັດສະດຸ (ສູງສຸດ 100)
  • strain: ຄວາມເຄັ່ງຮູບ
  • stress: ຄວາມເຄັ້ນ 2nd Piola-Kirchhoff
  • fstat: ຕົວແປສະຖານະ
  • dtime: ຂັ້ນເພີ່ມເວລາ
  • ttime: ເວລາລວມໃນຈຸດເລີ່ມຂອງ ຂັ້ນເພີ່ມເວລາ ປັດຈຸບັນ
  • temperature: ອຸນຫະພູມ