Nënprograme për deformimin elastoplastik (uyield.f90)¶

Ofrohen nënprograme për llogaritjen e matricës së ngurtësisë elastoplastike dhe të return mapping të sforcimit. Kur përdoret një funksion rrjedhjeje i përcaktuar nga përdoruesi, fillimisht vendosni !PLASTIC, YIELD=USER në skedarin hyrës dhe jepni konstantat e nevojshme të materialit, pastaj duhet të krijohen funksioni uElastoPlasticNumStatus dhe nënprogramet uElastoPlasticMatrix dhe uBackwardEuler .

(1) Funksioni që kthen numrin e ndryshoreve reale të gjendjes¶

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

• matl: Varg që ruan konstantat e materialit (1-100: konstanta materiali të përcaktuara nga sistemi, 101-200: konstanta materiali të përcaktuara nga përdoruesi)

(2) Nënprogrami për llogaritjen e matricës së ngurtësisë elastoplastike¶

  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: Varg që ruan konstantat e materialit (1-100: konstanta materiali të përcaktuara nga sistemi, 101-200: konstanta materiali të përcaktuara nga përdoruesi)
• stress: Sforcimi i dytë Piola-Kirchhoff
• istat: Ndryshore e plotë e gjendjes
• fstat: Varg i ndryshoreve reale të gjendjes
• plstrain: Deformimi plastik në fillim të nën-hapit aktual
• D: Matrica elastoplastike
• temp: Temperatura
• hdflag: Llogarit të gjithë komponentët (0), vetëm komponentët deviatorikë (1), ose vetëm komponentët vëllimorë (2)

(3) Nënprogrami për llogaritjen e return mapping të sforcimit¶

  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: Varg që ruan konstantat e materialit (1-100: konstanta materiali të përcaktuara nga sistemi, 101-200: konstanta materiali të përcaktuara nga përdoruesi)
• stress: Sforcimi i dytë Piola-Kirchhoff i marrë duke supozuar deformim elastik për sforcimin provë
• plstrain: Deformimi plastik në fillim të nën-hapit aktual
• istat: Ndryshore e plotë e gjendjes
• fstat: Varg i ndryshoreve reale të gjendjes
• temp: Temperatura
• hdflag: Llogarit të gjithë komponentët (0), vetëm komponentët deviatorikë (1), ose vetëm komponentët vëllimorë (2)