Skip to content

ຊັບຣູທີນທີ່ກ່ຽວກັບການປ່ຽນຮູບອີລາສໂຕພລາສຕິກ (uyield.f90)

ມີຊັບຣູທີນສຳລັບຄຳນວນເມທຣິກຄວາມແຂງອີລາສໂຕພລາສຕິກ ແລະ ການແມັບກັບຄືນ ຂອງຄວາມເຄັ້ນ. ເມື່ອໃຊ້ຟັງຊັນຄ່າຄອບຄົວທີ່ຜູ້ໃຊ້ກຳນົດ, ກ່ອນອື່ນໃຫ້ຕັ້ງ !PLASTIC, YIELD=USER ໃນໄຟລ໌ຂາເຂົ້າ ແລະໃສ່ຄ່າຄົງທີ່ວັດສະດຸທີ່ຈຳເປັນ, ຈາກນັ້ນຈຳເປັນຕ້ອງສ້າງຟັງຊັນ uElastoPlasticNumStatus ແລະຊັບຣູທີນ uElastoPlasticMatrix ແລະ uBackwardEuler .

(1) ຟັງຊັນທີ່ສົ່ງຄືນຈຳນວນຕົວແປສະຖານະຈຳນວນຈິງ

integer function uElastoPlasticNumStatus( matl )
    real(kind=kreal),   intent(in)  :: matl(:)
  • matl: ອາເຣທີ່ເກັບຄ່າຄົງທີ່ວັດສະດຸ (1-100: ຄ່າຄົງທີ່ວັດສະດຸທີ່ລະບົບກຳນົດ, 101-200: ຄ່າຄົງທີ່ວັດສະດຸທີ່ຜູ້ໃຊ້ກຳນົດ)

(2) ຊັບຣູທີນສຳລັບຄຳນວນເມທຣິກຄວາມແຂງອີລາສໂຕພລາສຕິກ

  subroutine uElastoPlasticMatrix( matl, stress, istat, fstat, plstrain, D, temp, hdflag )
    real(kind=kreal),   intent(in)  :: matl(:)
    real(kind=kreal),   intent(in)  :: stress(6)
    integer(kind=kint), intent(in)  :: istat
    real(kind=kreal),   intent(in)  :: fstat(:)
    real(kind=kreal),   intent(in)  :: plstrain
    real(kind=kreal),   intent(out) :: D(:,:)
    real(kind=kreal),   intent(in)  :: temp
    integer(kind=kint), intent(in)  :: hdflag
  • matl: ອາເຣທີ່ເກັບຄ່າຄົງທີ່ວັດສະດຸ (1-100: ຄ່າຄົງທີ່ວັດສະດຸທີ່ລະບົບກຳນົດ, 101-200: ຄ່າຄົງທີ່ວັດສະດຸທີ່ຜູ້ໃຊ້ກຳນົດ)
  • stress: ຄວາມເຄັ້ນ 2nd Piola-Kirchhoff
  • istat: ຕົວແປສະຖານະຈຳນວນເຕັມ
  • fstat: ອາເຣຂອງຕົວແປສະຖານະຈຳນວນຈິງ
  • plstrain: ຄວາມເຄັ່ງພລາສຕິກທີ່ຕົ້ນສັບສະເຕັບປັດຈຸບັນ
  • D: ເມທຣິກອີລາສໂຕພລາສຕິກ
  • temp: ອຸນຫະພູມ
  • hdflag: ຄຳນວນອົງປະກອບທັງໝົດ (0), ສະເພາະອົງປະກອບ ເບື່ອງເບນ (1), ຫຼືສະເພາະອົງປະກອບປະລິມາດ (2)

(3) ຊັບຣູທີນສຳລັບຄຳນວນ ການແມັບກັບຄືນ ຂອງຄວາມເຄັ້ນ

  subroutine uBackwardEuler( matl, stress, plstrain, istat, fstat, temp, hdflag )
    real(kind=kreal),   intent(in)    :: matl(:)
    real(kind=kreal),   intent(inout) :: stress(6)
    real(kind=kreal),   intent(in)    :: plstrain
    integer(kind=kint), intent(inout) :: istat
    real(kind=kreal),   intent(inout) :: fstat(:)
    real(kind=kreal),   intent(in)    :: temp
    integer(kind=kint), intent(in)    :: hdflag
  • matl: ອາເຣທີ່ເກັບຄ່າຄົງທີ່ວັດສະດຸ (1-100: ຄ່າຄົງທີ່ວັດສະດຸທີ່ລະບົບກຳນົດ, 101-200: ຄ່າຄົງທີ່ວັດສະດຸທີ່ຜູ້ໃຊ້ກຳນົດ)
  • stress: ຄວາມເຄັ້ນ 2nd Piola-Kirchhoff ທີ່ໄດ້ຮັບໂດຍສົມມຸດການປ່ຽນຮູບຍືດຫຍຸ່ນສຳລັບ ຄວາມເຄັ້ນທົດລອງ
  • plstrain: ຄວາມເຄັ່ງພລາສຕິກທີ່ຕົ້ນສັບສະເຕັບປັດຈຸບັນ
  • istat: ຕົວແປສະຖານະຈຳນວນເຕັມ
  • fstat: ອາເຣຂອງຕົວແປສະຖານະຈຳນວນຈິງ
  • temp: ອຸນຫະພູມ
  • hdflag: ຄຳນວນອົງປະກອບທັງໝົດ (0), ສະເພາະອົງປະກອບ ເບື່ອງເບນ (1), ຫຼືສະເພາະອົງປະກອບປະລິມາດ (2)