sigeps22c.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
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)

Function/Subroutine Documentation

sigeps22c()

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

Definition at line 30 of file sigeps22c.F.

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