sigeps122c.F File Reference

#include "implicit_f.inc"
#include "tabsiz_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sigeps122c (nel, nuparam, nuvar, uparam, uvar, epsxx, epsyy, rho, pla, dpla, depsxx, depsyy, depsxy, depsyz, depszx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, thk, thkly, off, sigy, etse, dmg, seq, shf, soundsp, asrate, epsd_pg, nfunc, ifunc, npf, tf, nvartmp, vartmp, ioff_duct)

Function/Subroutine Documentation

sigeps122c()

subroutine sigeps122c	(	integer, intent(in)	nel,
		integer, intent(in)	nuparam,
		integer, intent(in)	nuvar,
		dimension(nuparam), intent(in)	uparam,
		intent(inout)	uvar,
		intent(in)	epsxx,
		intent(in)	epsyy,
		intent(in)	rho,
		intent(inout)	pla,
		intent(inout)	dpla,
		intent(in)	depsxx,
		intent(in)	depsyy,
		intent(in)	depsxy,
		intent(in)	depsyz,
		intent(in)	depszx,
		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)	thk,
		intent(in)	thkly,
		intent(inout)	off,
		intent(out)	sigy,
		intent(out)	etse,
		intent(inout)	dmg,
		intent(inout)	seq,
		intent(in)	shf,
		intent(out)	soundsp,
		intent(in)	asrate,
		intent(in)	epsd_pg,
		integer, intent(in)	nfunc,
		integer, dimension(nfunc), intent(in)	ifunc,
		integer, dimension(snpc), intent(in)	npf,
		dimension(stf), intent(in)	tf,
		integer, intent(in)	nvartmp,
		integer, dimension(nel,nvartmp), intent(inout)	vartmp,
		integer, dimension(nel), intent(inout)	ioff_duct )

Definition at line 33 of file sigeps122c.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      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      "tabsiz_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL,NUPARAM,NUVAR,
     .   NFUNC,IFUNC(NFUNC),NPF(SNPC),NVARTMP
      my_real, INTENT(IN) :: asrate
      my_real, INTENT(IN) ::
     .   uparam(nuparam),tf(stf)
      my_real,DIMENSION(NEL), INTENT(IN) :: 
     .   rho,epsxx,epsyy,
     .   depsxx,depsyy,depsxy,depsyz,depszx,
     .   sigoxx,sigoyy,sigoxy,sigoyz,sigozx,
     .   thkly,shf
      my_real ,DIMENSION(NEL), INTENT(OUT)   :: 
     .   soundsp,signxx,signyy,signxy,signyz,signzx,
     .   sigy,etse
      my_real ,DIMENSION(NEL), INTENT(IN)    :: epsd_pg
      my_real ,DIMENSION(NEL), INTENT(INOUT) :: pla,thk,off,seq,dpla
      my_real ,DIMENSION(NEL,6), INTENT(INOUT) :: dmg
      my_real ,DIMENSION(NEL,NUVAR), INTENT(INOUT) :: 
     .   uvar
      INTEGER ,DIMENSION(NEL,NVARTMP), INTENT(INOUT) :: VARTMP
      INTEGER, DIMENSION(NEL), INTENT(INOUT) :: IOFF_DUCT
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IRES,I
      my_real ,DIMENSION(NEL) :: epsd
C=======================================================================
c
      ires = nint(uparam(17)) ! Plastic projection method
                              !  = 1 => Nice method
                              !  = 2 => Cutting plane
                              
      ! strain rate filtering, using global element strain rate
      do i=1,nel
        epsd(i) = asrate*epsd_pg(i) + (one-asrate)*uvar(i,18)
        uvar(i,18) = epsd(i)
      end do
c--------------------------                        
      SELECT CASE (ires)
c      
        CASE(1)
c     
          CALL mat122c_nice(
     1         nel     ,nuparam ,nuvar   ,uparam  ,uvar    ,
     2         epsxx   ,epsyy   ,rho     ,pla     ,dpla    ,
     3         depsxx  ,depsyy  ,depsxy  ,depsyz  ,depszx  ,
     4         sigoxx  ,sigoyy  ,sigoxy  ,sigoyz  ,sigozx  ,
     5         signxx  ,signyy  ,signxy  ,signyz  ,signzx  ,
     6         thk     ,thkly   ,off     ,sigy    ,etse    ,
     7         dmg     ,seq     ,epsd    ,shf     ,soundsp ,
     8         nfunc   ,ifunc   ,npf     ,tf      ,nvartmp ,
     9         vartmp  )
c      
        CASE(2)
c
          CALL mat122c_newton(
     1         nel     ,nuparam ,nuvar   ,uparam  ,uvar    ,        
     2         epsxx   ,epsyy   ,rho     ,pla     ,dpla    ,
     3         depsxx  ,depsyy  ,depsxy  ,depsyz  ,depszx  ,
     4         sigoxx  ,sigoyy  ,sigoxy  ,sigoyz  ,sigozx  ,
     5         signxx  ,signyy  ,signxy  ,signyz  ,signzx  ,
     6         thk     ,thkly   ,off     ,sigy    ,etse    ,
     7         dmg     ,seq     ,epsd    ,shf     ,soundsp ,
     8         nfunc   ,ifunc   ,npf     ,tf      ,nvartmp ,
     9         vartmp  )
    
      END SELECT     
c
      ! Ductile failure activation
      ioff_duct(1:nel) = 1  
c  
c-----------