Mga subroutine para sa elastikong deformasyon (uelastic.f90)¶

Nagbibigay ng mga subroutine para kalkulahin ang elastikong matris ng katigasan at i-update ang stress sa mga problemang elastiko at hiperelastiko. Upang gumamit ng constitutive relation na elastiko o hiperelastiko na tinukoy ng user, itakda muna ang !ELASTIC, TYPE=USER o !HYPERELASTIC, TYPE=USER sa input file at ilagay ang kinakailangang mga konstante ng materyal, at pagkatapos ay buuin ang mga subroutine na uElasticMatrix at uElasticUpdate.

(1) Subroutine para sa pagkalkula ng elastikong matris ng katigasan¶

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: array na nag-iimbak ng mga konstante ng materyal (maximum 100)
• strain: Green-Lagrange na deformasyon
• D: elastikong matris

(2) Subroutine para sa pagkalkula ng stress¶

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: array na nag-iimbak ng mga konstante ng materyal (maximum 100)
• strain: Green-Lagrange na deformasyon
• stress: stress