لچک-پلاسٹک تبدیلیِ شکل سے متعلق سب روٹینز (uyield.f90)¶
لچک-پلاسٹک سختی میٹرکس اور تناؤ کی ریٹرن میپنگ کا حساب کرنے کے لیے سب روٹینز فراہم کی گئی ہیں۔ صارف کی متعین کردہ تسلیم تفاعل استعمال کرنے کے لیے پہلے ان پٹ فائل میں !PLASTIC, YIELD=USER متعین کر کے درکار مادی مستقلات درج کریں، پھر تفاعل uElastoPlasticNumStatus اور سب روٹینز uElastoPlasticMatrix اور uBackwardEuler تیار کرنا ضروری ہے۔
(1) حقیقی عددی حالت کے متغیرات کی تعداد واپس کرنے والا تفاعل¶
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: دوسرا 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: لچک دار تبدیلیِ شکل فرض کر کے حاصل کیا گیا آزمائشی تناؤ، یعنی دوسرا Piola-Kirchhoff تناؤplstrain: موجودہ ذیلی مرحلے کے آغاز پر پلاسٹک کرنشistat: صحیح عددی حالت کا متغیرfstat: حقیقی عددی حالت کے متغیرات کی صفtemp: درجۂ حرارتhdflag: تمام اجزا (0)، صرف انحرافی اجزا (1)، یا صرف حجمی اجزا (2) کا حساب