Undirforrit fyrir teygjuplastíska formbreytingu (uyield.f90)¶

Undirforrit eru tiltæk til að reikna teygjuplastíska stífleikafylkið og vörpun spennu aftur á flotflöt. Þegar notendaskilgreint flotfall er notað skal fyrst stilla !PLASTIC, YIELD=USER í inntaksskránni og færa inn nauðsynlega efnisfasta, síðan þarf að útfæra fallið uElastoPlasticNumStatus og undirforritin uElastoPlasticMatrix og uBackwardEuler.

(1) Fall sem skilar fjölda raungildra ástandsbreyta¶

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

• matl: Fylki sem geymir efnisfasta (1-100: kerfisskilgreindir efnisfastar, 101-200: notendaskilgreindir efnisfastar)

(2) Undirforrit til að reikna teygjuplastíska stífleikafylkið¶

  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: Fylki sem geymir efnisfasta (1-100: kerfisskilgreindir efnisfastar, 101-200: notendaskilgreindir efnisfastar)
• stress: 2. Piola-Kirchhoff-spenna
• istat: Heiltöluástandsbreyta
• fstat: Fylki raungildra ástandsbreyta
• plstrain: Plastísk formbreyting í upphafi núverandi undirskrefs
• D: Teygjuplastískt fylki
• temp: Hitastig
• hdflag: Reikna alla þætti (0), aðeins fráviksþætti (1) eða aðeins rúmmálsþætti (2)

(3) Undirforrit fyrir útreikning á vörpun spennu aftur á flotflöt¶

  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: Fylki sem geymir efnisfasta (1-100: kerfisskilgreindir efnisfastar, 101-200: notendaskilgreindir efnisfastar)
• stress: 2. Piola-Kirchhoff-spenna sem fæst með því að gera ráð fyrir teygjanlegri formbreytingu fyrir prófspennuna
• plstrain: Plastísk formbreyting í upphafi núverandi undirskrefs
• istat: Heiltöluástandsbreyta
• fstat: Fylki raungildra ástandsbreyta
• temp: Hitastig
• hdflag: Reikna alla þætti (0), aðeins fráviksþætti (1) eða aðeins rúmmálsþætti (2)