ელასტიკურ დეფორმაციასთან დაკავშირებული ქვეპროგრამები (uelastic.f90)¶

უზრუნველყოფილია ქვეპროგრამები ელასტიკურ და ჰიპერელასტიკურ ამოცანებში ელასტიკური სიხისტის მატრიცის გამოსათვლელად და ძაბვის გასაახლებლად. მომხმარებლის მიერ განსაზღვრული ელასტიკური ან ჰიპერელასტიკური კონსტიტუციური დამოკიდებულების გამოსაყენებლად, ჯერ შემავალ ფაილში დააყენეთ !ELASTIC, TYPE=USER ან !HYPERELASTIC, TYPE=USER და შეიტანეთ საჭირო მასალის მუდმივები, შემდეგ კი შექმენით uElasticMatrix და uElasticUpdate ქვეპროგრამები.

(1) ელასტიკური სიხისტის მატრიცის გამოთვლის ქვეპროგრამა¶

subroutine uElasticMatrix( matl, strain, D )
    REAL(KIND=kreal), INTENT(IN) :: matl(:)
    REAL(KIND=kreal), INTENT(IN) :: strain(6)
    REAL(KIND=kreal), INTENT(OUT) :: D(6,6)

• matl: მასალის მუდმივების შემნახველი მასივი (მაქსიმუმ 100)
• strain: Green-Lagrange-ის დეფორმაცია
• D: ელასტიკურობის მატრიცა

(2) ძაბვის გამოთვლის ქვეპროგრამა¶

subroutine uElasticUpdate ( matl, strain, stress )
    REAL(KIND=kreal), INTENT(IN) :: matl(:)
    REAL(KIND=kreal), INTENT(IN) :: strain(6)
    REAL(KIND=kreal), INTENT(OUT) :: stress(6)

• matl: მასალის მუდმივების შემნახველი მასივი (მაქსიმუმ 100)
• strain: Green-Lagrange-ის დეფორმაცია
• stress: ძაბვა