ซับรูทีนที่เกี่ยวข้องกับการเปลี่ยนรูปแบบยืดหยุ่น (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: ความเค้น