sigeps22c_8F_source.html

Copyright>        OpenRadioss

Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.

Copyright>

Copyright>        This program is free software: you can redistribute it and/or modify

Copyright>        it under the terms of the GNU Affero General Public License as published by

Copyright>        the Free Software Foundation, either version 3 of the License, or

Copyright>        (at your option) any later version.

Copyright>

Copyright>        This program is distributed in the hope that it will be useful,

Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of

Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

Copyright>        GNU Affero General Public License for more details.

Copyright>

Copyright>        You should have received a copy of the GNU Affero General Public License

Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.

Copyright>

Copyright>

Copyright>        Commercial Alternative: Altair Radioss Software

Copyright>

Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss

Copyright>        software under a commercial license.  Contact Altair to discuss further if the

Copyright>        commercial version may interest you: https://www.altair.com/radioss/.

!||====================================================================

!||    sigeps22c   ../engine/source/materials/mat/mat022/sigeps22c.F

!||--- called by ------------------------------------------------------

!||    mulawc      ../engine/source/materials/mat_share/mulawc.F90

!||--- calls      -----------------------------------------------------

!||    m22cplr     ../engine/source/materials/mat/mat022/m22cplr.F

!||====================================================================


      SUBROUTINE sigeps22c(

     1                   JFT      ,JLT    ,PM       ,THK      ,OFF     ,

     2                   SIGY     ,DT1C   ,IPLA     ,NEL      ,OFF_OLD ,

     3                   GS       ,DPLA   ,IOFF_DUCT,NPTT     ,IPT     ,

     4                   EPCHK    ,ALPE   ,THKLYL   ,IMAT     ,PLA     ,

     5                   DEPSXX   ,DEPSYY ,DEPSXY   ,DEPSYZ   ,DEPSZX  ,

     6                   SIGOXX   ,SIGOYY ,SIGOXY   ,SIGOYZ   ,SIGOZX  ,

     7                   SIGNXX   ,SIGNYY ,SIGNXY   ,SIGNYZ   ,SIGNZX  ,

     8                   INLOC    ,DPLANL ,LOFF     )

C-----------------------------------------------

C   I m p l i c i t   T y p e s

C-----------------------------------------------

#include      "implicit_f.inc"

C-----------------------------------------------

C   G l o b a l   P a r a m e t e r s

C-----------------------------------------------

#include      "mvsiz_p.inc"

C-----------------------------------------------

C   C o m m o n   B l o c k s

C-----------------------------------------------

#include      "param_c.inc"

C-----------------------------------------------

C   D u m m y   A r g u m e n t s

C-----------------------------------------------

      INTEGER JFT,JLT,IPLA,NEL,NPT,NPTT,IMAT,IPT

      INTEGER IOFF_DUCT(MVSIZ),INLOC

C     REAL

      my_real

     .   PM(NPROPM,*),THK(*),OFF(*),OFF_OLD(*),SIGY(*),DT1C(*),GS(*)

      my_real

     .   DPLA(*),PLA(NEL),

     .   SIGOXX(NEL),SIGOYY(NEL),SIGOXY(NEL),SIGOYZ(NEL),

     .   sigozx(nel),signxx(nel),signyy(nel),signxy(nel),

     .   signyz(nel),signzx(nel),depsxx(mvsiz),depsyy(mvsiz),

     .   depsxy(mvsiz),depsyz(mvsiz),depszx(mvsiz),epchk(mvsiz),

     .   alpe(mvsiz),thklyl(nel),dplanl(nel)

      my_real, DIMENSION(NEL), INTENT(IN) :: loff

C-----------------------------------------------

C   L o c a l   V a r i a b l e s

C-----------------------------------------------

      INTEGER ICC(MVSIZ),I,ILAYER,ICC_1,J

C     REAL

      my_real

     .   A1(MVSIZ),YLD(MVSIZ),

     .   ca(mvsiz),cb(mvsiz),cn(mvsiz),ymax(mvsiz),

     .   young(mvsiz),cc(mvsiz),epdr(mvsiz),nu(mvsiz),

     .   epsl(mvsiz),hl(mvsiz),ym(mvsiz),yldl(mvsiz),ezz(mvsiz),

     .   ym_1, a1_1, ca_1, cb_1, cn_1,

     .   epmx, ymax_1, cc_1, epsl_1,

     .   hl_1, yldl_1,nu_1

C-----------------------------------------------

!

      nu_1   = pm(21,imat)

      ym_1   = pm(20,imat)

      a1_1   = pm(24,imat)

      ca_1   = pm(38,imat)

      cb_1   = pm(39,imat)

      cn_1   = pm(40,imat)

      epmx   = pm(41,imat)

      ymax_1 = pm(42,imat)

      cc_1   = pm(43,imat)

      epsl_1 = pm(45,imat)

      hl_1   = pm(46,imat)

      yldl_1 = pm(47,imat)

      icc_1  = nint(pm(49,imat))

!

      DO i=jft,jlt

        ym(i)   = ym_1

        a1(i)   = a1_1

        ca(i)   = ca_1

        cb(i)   = cb_1

        cn(i)   = cn_1

        ymax(i) = ymax_1

        cc(i)   = cc_1

        epdr(i) = max(em20,pm(44,imat)*dt1c(i))

        epsl(i) = epsl_1

        hl(i)   = hl_1

        yldl(i) = yldl_1

        icc(i)  = icc_1

        nu(i) = nu_1

        ezz(i) = zero

      ENDDO

C-----------------------------------------------------------

C     PLASTICITE AUX POINTS D'INTEGRATION

C-----------------------------------------------------------

C     CONTRAINTES PLASTIQUEMENT ADMISSIBLES

C-----------------------------------------------------------

        CALL m22cplr(jft    ,jlt    ,ezz    ,off    ,pla   ,

     2               ipla   ,gs     ,yld    ,a1     ,ca    ,

     3               cb     ,cn     ,ymax   ,nu     ,dpla  ,

     4               epchk  ,young  ,cc     ,epdr   ,epsl  ,

     5               hl     ,ym     ,yldl   ,alpe   ,icc   ,

     6               depsxx ,depsyy ,depsxy ,depsyz ,depszx,

     7               sigoxx ,sigoyy ,sigoxy ,sigoyz ,sigozx,

     8               signxx ,signyy ,signxy ,signyz ,signzx,

     9               nel    )

!

        IF (ipt == 1+nptt/2) THEN

          DO i=jft,jlt

            sigy(i)=yld(i)

          ENDDO

        ENDIF

C----------------------------

C     TEST DE RUPTURE DUCTILE

C----------------------------

        DO i=jft,jlt

          IF (off(i) == off_old(i) .and. off(i) > zero) THEN

            IF (off(i) == one .and. epchk(i) >= epmx) THEN

              off(i)= four_over_5

              ioff_duct(i) = 1    ! start de la rupture

            ELSE IF (off(i) < one ) THEN

              off(i) = off(i)*four_over_5

            ENDIF

          ENDIF

        ENDDO

!---------------------------

!     thickness update:

!---------------------------

      DO i=jft,jlt

        IF (inloc > 0) THEN

          IF (loff(i) == one) THEN

            ezz(i) = -nu(i)*(signxx(i)-sigoxx(i)+signyy(i)-sigoyy(i))/young(i)

            ezz(i) = ezz(i) - max(dplanl(i),zero)*half*(signxx(i)+signyy(i))/yld(i)

          ENDIF

        ELSE

          ezz(i) = -(depsxx(i)+depsyy(i))*nu(i)-(one - two*nu(i))*ezz(i)

          ezz(i) = ezz(i)/(one-nu(i))

        ENDIF

        thk(i) = thk(i) + ezz(i) * thklyl(i)*off(i)

      ENDDO

c-----------

      RETURN


      END

ymax
subroutine ymax(idn, fac, npc, pld, stiffmin, stiffmax, stiffini, stiffavg)
Definition law100_upd.F:272

m22cplr
subroutine m22cplr(jft, jlt, ezz, off, epseq, ipla, gs, yld, a1, ca, cb, cn, ymax, nu, dpla, epchk, young, cc, epdr, epsl, hl, ym, yldl, alpe, icc, depsxx, depsyy, depsxy, depsyz, depszx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, nel)
Definition m22cplr.F:37

max
#define max(a, b)
Definition macros.h:21

sigeps22c
subroutine sigeps22c(jft, jlt, pm, thk, off, sigy, dt1c, ipla, nel, off_old, gs, dpla, ioff_duct, nptt, ipt, epchk, alpe, thklyl, imat, pla, depsxx, depsyy, depsxy, depsyz, depszx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, inloc, dplanl, loff)
Definition sigeps22c.F:39