Subroutine zinazohusiana na mgeuko elastoplastiki (uyield.f90)¶

Hutoa subroutine za kukokotoa matriki ya ugumu elastoplastiki na return mapping ya mkazo. Unapotumia yield function iliyofafanuliwa na mtumiaji, kwanza weka !PLASTIC, YIELD=USER kwenye faili ya ingizo na uingize material constants zinazohitajika, kisha unda function uElastoPlasticNumStatus na subroutine uElastoPlasticMatrix pamoja na uBackwardEuler.

(1) Function inayorudisha idadi ya state variables za real¶

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

• matl: Array inayohifadhi material constants (1-100: material constants zilizofafanuliwa na mfumo, 101-200: material constants zilizofafanuliwa na mtumiaji)

(2) Subroutine ya kukokotoa matriki ya ugumu elastoplastiki¶

  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 inayohifadhi material constants (1-100: material constants zilizofafanuliwa na mfumo, 101-200: material constants zilizofafanuliwa na mtumiaji)
• stress: Mkazo wa 2nd Piola-Kirchhoff
• istat: State variable ya nambari kamili
• fstat: Array ya state variables za real
• plstrain: Mgeuko plastiki mwanzoni mwa substep ya sasa
• D: Matriki ya elastoplastiki
• temp: Halijoto
• hdflag: Hukokotoa vipengele vyote (0), vipengele vya deviatoric pekee (1), au vipengele vya volumetric pekee (2)

(3) Subroutine ya kukokotoa Return mapping ya mkazo¶

  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 inayohifadhi material constants (1-100: material constants zilizofafanuliwa na mfumo, 101-200: material constants zilizofafanuliwa na mtumiaji)
• stress: Mkazo wa 2nd Piola-Kirchhoff uliopatikana kama trial stress kwa kudhani mgeuko elastiki
• plstrain: Mgeuko plastiki mwanzoni mwa substep ya sasa
• istat: State variable ya nambari kamili
• fstat: Array ya state variables za real
• temp: Halijoto
• hdflag: Hukokotoa vipengele vyote (0), vipengele vya deviatoric pekee (1), au vipengele vya volumetric pekee (2)