Chương trình con liên quan đến biến dạng đàn hồi (uelastic.f90)

Cung cấp các chương trình con để tính ma trận độ cứng đàn hồi và cập nhật ứng suất cho các bài toán đàn hồi và siêu đàn hồi. Khi sử dụng quan hệ cấu thành đàn hồi hoặc siêu đàn hồi do người dùng định nghĩa, trước tiên cần thiết lập !ELASTIC, TYPE=USER hoặc !HYPERELASTIC, TYPE=USER trong tệp đầu vào và nhập các hằng số vật liệu cần thiết, sau đó tạo các chương trình con uElasticMatrix và uElasticUpdate.

(1) Chương trình con tính ma trận độ cứng đàn hồi

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: Mảng lưu các hằng số vật liệu (tối đa 100)
• strain: Biến dạng Green-Lagrange
• D: Ma trận đàn hồi

(2) Chương trình con tính ứng suất

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: Mảng lưu các hằng số vật liệu (tối đa 100)
• strain: Biến dạng Green-Lagrange
• stress: Ứng suất