sigeps112.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "comlock.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sigeps112 (nel, ngl, nuparam, nuvar, time, timestep, uparam, uvar, jthe, off, rho0, rho, pla, dpla, epsd, soundsp, epszz, depsxx, depsyy, depszz, depsxy, depsyz, depszx, sigoxx, sigoyy, sigozz, sigoxy, sigoyz, sigozx, signxx, signyy, signzz, signxy, signyz, signzx, sigy, et, nvartmp, numtabl, vartmp, itable, table)

Function/Subroutine Documentation

sigeps112()

subroutine sigeps112	(	integer	nel,
		integer, dimension(nel), intent(in)	ngl,
		integer	nuparam,
		integer	nuvar,
			time,
			timestep,
		intent(in)	uparam,
		intent(inout)	uvar,
		integer	jthe,
		intent(inout)	off,
		intent(in)	rho0,
		intent(in)	rho,
		intent(inout)	pla,
		intent(inout)	dpla,
		intent(inout)	epsd,
		intent(out)	soundsp,
		intent(in)	epszz,
		intent(in)	depsxx,
		intent(in)	depsyy,
		intent(in)	depszz,
		intent(in)	depsxy,
		intent(in)	depsyz,
		intent(in)	depszx,
		intent(in)	sigoxx,
		intent(in)	sigoyy,
		intent(in)	sigozz,
		intent(in)	sigoxy,
		intent(in)	sigoyz,
		intent(in)	sigozx,
		intent(out)	signxx,
		intent(out)	signyy,
		intent(out)	signzz,
		intent(out)	signxy,
		intent(out)	signyz,
		intent(out)	signzx,
		intent(out)	sigy,
		intent(out)	et,
		integer	nvartmp,
		integer	numtabl,
		integer, dimension(nel,nvartmp)	vartmp,
		integer, dimension(numtabl)	itable,
		type(ttable), dimension(ntable)	table )

Definition at line 35 of file sigeps112.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE table_mod
      USE interface_table_mod
      USE elbufdef_mod     
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C O M M O N
C-----------------------------------------------
#include      "com04_c.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL,NUPARAM,NUVAR,JTHE,NUMTABL,NVARTMP,
     .        ITABLE(NUMTABL)
      INTEGER ,DIMENSION(NEL), INTENT(IN) :: NGL
      my_real 
     .   time,timestep
      INTEGER :: VARTMP(NEL,NVARTMP)
      my_real,DIMENSION(NUPARAM), INTENT(IN) :: 
     .   uparam
      my_real,DIMENSION(NEL), INTENT(IN)     :: 
     .   rho0,rho,epszz,
     .   depsxx,depsyy,depszz,depsxy,depsyz,depszx,
     .   sigoxx,sigoyy,sigozz,sigoxy,sigoyz,sigozx
      my_real ,DIMENSION(NEL), INTENT(OUT)   :: 
     .   soundsp,sigy,et,
     .   signxx,signyy,signzz,signxy,signyz,signzx
      my_real ,DIMENSION(NEL), INTENT(INOUT)       :: 
     .   dpla,off
      my_real ,DIMENSION(NEL,4), INTENT(INOUT)       :: 
     .   pla,epsd
      my_real ,DIMENSION(NEL,NUVAR), INTENT(INOUT) :: 
     .   uvar
c
      TYPE(TTABLE), DIMENSION(NTABLE) ::  TABLE
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IRES
C=======================================================================
      ires  = nint(uparam(13)) ! Plastic projection method
                               !  = 1 => Nice method
                               !  = 2 => Newton-iteration method
C=======================================================================
        ! Choice of return mapping algorithm
        SELECT CASE (ires)
          ! Nice explicit return mapping      
          CASE(1)        
            CALL mat112_xia_nice(
     1         nel     ,ngl     ,nuparam ,nuvar   ,time    ,timestep,
     2         uparam  ,uvar    ,jthe    ,off     ,rho0    ,rho     ,
     3         pla     ,dpla    ,epsd    ,soundsp ,epszz   ,
     4         depsxx  ,depsyy  ,depszz  ,depsxy  ,depsyz  ,depszx  ,
     5         sigoxx  ,sigoyy  ,sigozz  ,sigoxy  ,sigoyz  ,sigozx  ,
     6         signxx  ,signyy  ,signzz  ,signxy  ,signyz  ,signzx  ,
     7         sigy    ,et      ,
     8         nvartmp ,numtabl ,vartmp  ,itable  ,table   )     
          ! Newton cutting plane return mapping
          CASE(2)
            CALL mat112_xia_newton(
     1         nel     ,ngl     ,nuparam ,nuvar   ,time    ,timestep,
     2         uparam  ,uvar    ,jthe    ,off     ,rho0    ,rho     ,
     3         pla     ,dpla    ,epsd    ,soundsp ,epszz   ,
     4         depsxx  ,depsyy  ,depszz  ,depsxy  ,depsyz  ,depszx  ,
     5         sigoxx  ,sigoyy  ,sigozz  ,sigoxy  ,sigoyz  ,sigozx  ,
     6         signxx  ,signyy  ,signzz  ,signxy  ,signyz  ,signzx  ,
     7         sigy    ,et      ,
     8         nvartmp ,numtabl ,vartmp  ,itable  ,table   )
        END SELECT        
c-----------