ელასტოპლასტიკურ დეფორმაციასთან დაკავშირებული ქვეპროგრამები (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) გამოთვლა