跳轉至

彈性變形相關副程式（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：應力