Deformazio elastoplastikorako azpirrutinak (uyield.f90)¶

Zurruntasun-matrize elastoplastikoa eta tentsioaren itzulera-mapatzea kalkulatzeko azpirrutinak eskaintzen dira. Erabiltzaileak definitutako isurpen-funtzio bat erabiltzean, lehenik !PLASTIC, YIELD=USER ezarri sarrera-fitxategian eta beharrezko material-konstanteak sartu; ondoren, uElastoPlasticNumStatus funtzioa eta uElastoPlasticMatrix eta uBackwardEuler azpirrutinak sortu behar dira.

(1) Egoera-aldagai errealen kopurua itzultzen duen funtzioa¶

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

• matl: material-konstanteak gordetzen dituen arraya (1–100: sistemak definitutako material-konstanteak, 101–200: erabiltzaileak definitutako material-konstanteak)

(2) Zurruntasun-matrize elastoplastikoa kalkulatzeko azpirrutina¶

  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: material-konstanteak gordetzen dituen arraya (1–100: sistemak definitutako material-konstanteak, 101–200: erabiltzaileak definitutako material-konstanteak)
• stress: 2. Piola-Kirchhoff tentsioa
• istat: zenbaki osoko egoera-aldagaia
• fstat: egoera-aldagai errealen arraya
• plstrain: deformazio plastikoa uneko azpiurratsaren hasieran
• D: matrize elastoplastikoa
• temp: tenperatura
• hdflag: osagai guztiak (0), osagai desbideratzaileak bakarrik (1) edo osagai bolumetrikoak bakarrik (2) kalkulatzen ditu

(3) Tentsioaren itzulera-mapatzea kalkulatzeko azpirrutina¶

  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: material-konstanteak gordetzen dituen arraya (1–100: sistemak definitutako material-konstanteak, 101–200: erabiltzaileak definitutako material-konstanteak)
• stress: probako tentsioarentzat deformazio elastikoa suposatuz lortutako 2. Piola-Kirchhoff tentsioa
• plstrain: deformazio plastikoa uneko azpiurratsaren hasieran
• istat: zenbaki osoko egoera-aldagaia
• fstat: egoera-aldagai errealen arraya
• temp: tenperatura
• hdflag: osagai guztiak (0), osagai desbideratzaileak bakarrik (1) edo osagai bolumetrikoak bakarrik (2) kalkulatzen ditu