Mga subroutine para sa elastoplastic deformation (uyield.f90)¶

May mga subroutine para sa pagkalkula ng elastoplastic stiffness matrix at stress return mapping. Kapag gumagamit ng yield function na tinukoy ng user, itakda muna ang !PLASTIC, YIELD=USER sa input file at ilagay ang kinakailangang material constants, pagkatapos ay kailangang gawin ang function na uElastoPlasticNumStatus at ang mga subroutine na uElastoPlasticMatrix at uBackwardEuler.

(1) Function na nagbabalik ng bilang ng real state variables¶

integer function uElastoPlasticNumStatus( matl )
    real(kind=kreal),   intent(in)  :: matl(:)

• matl: array na naglalaman ng material constants (1–100: material constants na tinukoy ng system, 101–200: material constants na tinukoy ng user)

(2) Subroutine para sa pagkalkula ng elastoplastic stiffness matrix¶

  subroutine uElastoPlasticMatrix( matl, stress, istat, fstat, plstrain, D, temp, hdflag )
    real(kind=kreal),   intent(in)  :: matl(:)
    real(kind=kreal),   intent(in)  :: stress(6)
    integer(kind=kint), intent(in)  :: istat
    real(kind=kreal),   intent(in)  :: fstat(:)
    real(kind=kreal),   intent(in)  :: plstrain
    real(kind=kreal),   intent(out) :: D(:,:)
    real(kind=kreal),   intent(in)  :: temp
    integer(kind=kint), intent(in)  :: hdflag

• matl: array na naglalaman ng material constants (1–100: material constants na tinukoy ng system, 101–200: material constants na tinukoy ng user)
• stress: ika-2 Piola-Kirchhoff stress
• istat: integer state variable
• fstat: array ng real state variables
• plstrain: plastic strain sa simula ng kasalukuyang substep
• D: elastoplastic matrix
• temp: temperatura
• hdflag: kalkulahin ang lahat ng component (0), deviatoric components lamang (1), o volumetric components lamang (2)

(3) Subroutine para sa pagkalkula ng stress return mapping¶

  subroutine uBackwardEuler( matl, stress, plstrain, istat, fstat, temp, hdflag )
    real(kind=kreal),   intent(in)    :: matl(:)
    real(kind=kreal),   intent(inout) :: stress(6)
    real(kind=kreal),   intent(in)    :: plstrain
    integer(kind=kint), intent(inout) :: istat
    real(kind=kreal),   intent(inout) :: fstat(:)
    real(kind=kreal),   intent(in)    :: temp
    integer(kind=kint), intent(in)    :: hdflag

• matl: array na naglalaman ng material constants (1–100: material constants na tinukoy ng system, 101–200: material constants na tinukoy ng user)
• stress: ika-2 Piola-Kirchhoff stress na nakuha sa pag-aakalang elastic deformation para sa trial stress
• plstrain: plastic strain sa simula ng kasalukuyang substep
• istat: integer state variable
• fstat: array ng real state variables
• temp: temperatura
• hdflag: kalkulahin ang lahat ng component (0), deviatoric components lamang (1), o volumetric components lamang (2)