sigeps107.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	sigeps107 (nel, ngl, nuparam, nuvar, time, timestep, uparam, uvar, jthe, off, rho0, rho, pla, dpla, epsd, soundsp, 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

sigeps107()

subroutine sigeps107	(	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)	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 sigeps107.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,
     .   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,6), 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 mat107_nice(
     1         nel     ,ngl     ,nuparam ,nuvar   ,time    ,timestep,
     2         uparam  ,uvar    ,jthe    ,off     ,rho0    ,rho     ,
     3         pla     ,dpla    ,epsd    ,soundsp ,
     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 mat107_newton(
     1         nel     ,ngl     ,nuparam ,nuvar   ,time    ,timestep,
     2         uparam  ,uvar    ,jthe    ,off     ,rho0    ,rho     ,
     3         pla     ,dpla    ,epsd    ,soundsp ,
     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-----------