কনটেন্টে যান

স্থিতিস্থাপক-প্লাস্টিক বিকৃতি-সম্পর্কিত সাবরুটিন (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) গণনা করে