sigeps119c.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sigeps119c (nel, ipt, npt, nuparam, nuvar, numtabl, itable, table, uvar, uparam, thk, thklyl, shf, ssp, off, iqeph, zcfac, depsxx, depsyy, depsxy, epsxx, epsyy, epsxy, epsyz, epszx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, wmc, thly, thk0)

Function/Subroutine Documentation

sigeps119c()

subroutine sigeps119c	(	integer, intent(in)	nel,
		integer, intent(in)	ipt,
		integer, intent(in)	npt,
		integer, intent(in)	nuparam,
		integer, intent(in)	nuvar,
		integer, intent(in)	numtabl,
		integer, dimension(numtabl), intent(in)	itable,
		type(ttable), dimension(ntable)	table,
		intent(inout)	uvar,
		intent(in)	uparam,
		intent(inout)	thk,
		intent(inout)	thklyl,
		intent(in)	shf,
		intent(out)	ssp,
		intent(in)	off,
		logical, intent(in)	iqeph,
		intent(inout)	zcfac,
		intent(in)	depsxx,
		intent(in)	depsyy,
		intent(in)	depsxy,
		intent(in)	epsxx,
		intent(in)	epsyy,
		intent(in)	epsxy,
		intent(in)	epsyz,
		intent(in)	epszx,
		intent(in)	sigoxx,
		intent(in)	sigoyy,
		intent(in)	sigoxy,
		intent(in)	sigoyz,
		intent(in)	sigozx,
		intent(out)	signxx,
		intent(out)	signyy,
		intent(out)	signxy,
		intent(out)	signyz,
		intent(out)	signzx,
		intent(inout)	wmc,
		intent(inout)	thly,
		intent(in)	thk0 )

Definition at line 33 of file sigeps119c.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE table_mod
      USE interface_table_mod
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      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      LOGICAL ,INTENT(IN) :: IQEPH
      INTEGER ,INTENT(IN) :: NEL,IPT,NPT,NUPARAM,NUVAR,NUMTABL
      INTEGER ,DIMENSION(NUMTABL) ,INTENT(IN)  :: ITABLE
      my_real ,DIMENSION(NUPARAM) ,INTENT(IN)  :: uparam
      my_real ,DIMENSION(NEL)     ,INTENT(IN)  :: shf,off,
     .   depsxx,depsyy,depsxy,epsxx,epsyy,epsxy,epsyz,epszx,
     .   sigoxx,sigoyy,sigoxy,sigoyz,sigozx,thk0
      my_real ,DIMENSION(MVSIZ,2),INTENT(INOUT) :: zcfac
      my_real ,DIMENSION(MVSIZ)  ,INTENT(INOUT) :: wmc
      my_real, DIMENSION(NEL)    ,INTENT(INOUT) :: thk,thly,thklyl
      my_real ,DIMENSION(NEL,NUVAR) ,INTENT(INOUT) :: uvar
      my_real ,DIMENSION(NEL) ,INTENT(OUT) :: ssp,
     .   signxx,signyy,signxy,signyz,signzx
      TYPE(TTABLE), DIMENSION(NTABLE) ::  TABLE
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I
      my_real :: nu,g12,gs,gc,a1c,a2c,tcoat,dezz
      my_real ,DIMENSION(NEL) :: et,thfac
C=======================================================================
      nu      = uparam(5)
      g12     = uparam(6)
      a1c     = uparam(14)
      a2c     = uparam(15)
      tcoat   = uparam(16)
      ssp     = uparam(17)
      gc      = uparam(20)
      gs      = g12*shf(1) 
c-------------------------------------------------------------------------
      IF (npt == 1) THEN ! membrane only
 
        CALL law119_membrane(
     .       nel    ,nuparam,nuvar  ,uparam ,uvar   ,
     .       gs     ,et     ,depsxx ,depsyy ,depsxy ,
     .       epsxx  ,epsyy  ,epsxy  ,epsyz  ,epszx  ,
     .       sigoxx ,sigoyy ,sigoxy ,sigoyz ,sigozx ,
     .       signxx ,signyy ,signxy ,signyz ,signzx ,
     .       numtabl,itable ,table  )
 
      ELSE  ! membrane + coating
 
        IF (ipt == 2) THEN
          CALL law119_membrane(
     .       nel    ,nuparam,nuvar  ,uparam ,uvar   ,
     .       gs     ,et     ,depsxx ,depsyy ,depsxy ,
     .       epsxx  ,epsyy  ,epsxy  ,epsyz  ,epszx  ,
     .       sigoxx ,sigoyy ,sigoxy ,sigoyz ,sigozx ,
     .       signxx ,signyy ,signxy ,signyz ,signzx ,
     .       numtabl,itable ,table  )
          thklyl(1:nel) = thk0(1:nel) - two*tcoat
          thly(1:nel)= one - two * tcoat / thk0(1:nel)
          wmc(1:nel) = zero
        ELSE    ! IPT = 1,3 => coating
          DO i=1,nel
            signxx(i) = a1c*epsxx(i) + a2c*epsyy(i)
            signyy(i) = a2c*epsxx(i) + a1c*epsyy(i)
            signxy(i) = gc*epsxy(i)
            signyz(i) = gc*epsyz(i)*shf(1)
            signzx(i) = gc*epszx(i)*shf(1)
          END DO
          thklyl(1:nel) = tcoat
          DO i=1,nel
            thly(i) = tcoat*(thk0(i)-tcoat) / (thk0(i)**2)
          END DO
          IF (ipt == 1) THEN
            wmc(1:nel) =-half*thly(1:nel)
          ELSE IF(ipt == 3) THEN
            wmc(1:nel) = half*thly(1:nel)
          END IF
        END IF
      END IF
c------------------------------------            
c       thickness update:
c------------------------------------            
      DO i=1,nel
        dezz = -(depsxx(i)+depsyy(i))*nu/(one-nu)
        thk(i) = thk(i) + dezz * thklyl(i)*off(i)
      ENDDO
C
C------for QEPH-------       
      IF (iqeph) THEN
        DO i=1,nel
          zcfac(i,1) = zcfac(i,1) + et(i) / npt
          zcfac(i,2) = min(et(i),zcfac(i,2))
        ENDDO
      ENDIF
c-----------
      RETURN