د الاستوپلاستیک بدلون اړوند فرعي‌روټینونه (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: د دویم 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) محاسبه کوي