跳转至

与弹性变形相关的子程序 (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: 应力