sigeps115.F File Reference

#include "implicit_f.inc"
#include "scr17_c.inc"
#include "com08_c.inc"
#include "units_c.inc"
#include "comlock.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sigeps115 (nel, ngl, nuparam, nuvar, grho, time, timestep, uparam, uvar, off, sigy, rho0, pla, dpla, soundsp, et, seq, depsxx, depsyy, depszz, depsxy, depsyz, depszx, sigoxx, sigoyy, sigozz, sigoxy, sigoyz, sigozx, signxx, signyy, signzz, signxy, signyz, signzx)

Function/Subroutine Documentation

sigeps115()

subroutine sigeps115	(	integer	nel,
		integer, dimension(nel), intent(in)	ngl,
		integer	nuparam,
		integer	nuvar,
		intent(in)	grho,
			time,
			timestep,
		intent(in)	uparam,
		intent(inout)	uvar,
		intent(inout)	off,
		intent(out)	sigy,
		intent(in)	rho0,
		intent(inout)	pla,
		intent(inout)	dpla,
		intent(out)	soundsp,
		intent(out)	et,
		intent(inout)	seq,
		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 )

Definition at line 31 of file sigeps115.F.

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      "scr17_c.inc"
#include      "com08_c.inc"
#include      "units_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
      INTEGER ,DIMENSION(NEL), INTENT(IN) :: NGL
      my_real 
     .   time,timestep
      my_real,DIMENSION(NUPARAM), INTENT(IN) :: 
     .   uparam
      my_real,DIMENSION(2*NEL), INTENT(IN)     :: 
     .   grho
      my_real,DIMENSION(NEL), INTENT(IN)     :: 
     .   rho0,
     .   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)       :: 
     .   pla,dpla,off,seq
      my_real ,DIMENSION(NEL,NUVAR), INTENT(INOUT) :: 
     .   uvar
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IRES,NINDX,I,J,INDX(NEL)
      my_real
     .   cfail,pfail
C=======================================================================
c
      ires  = nint(uparam(11)) ! Plastic projection method
                               !  = 1 => Nice method
                               !  = 2 => Newton
      cfail = uparam(14)       ! Tensile volumetric strain at failure
      pfail = uparam(15)       ! Maximum principal stress at failure
c
c--------------------------                              
      SELECT CASE (ires)
c      
        CASE(1) 
          CALL mat115_nice(
     1         nel     ,ngl     ,nuparam ,nuvar   ,grho    , 
     2         time    ,timestep,uparam  ,uvar    ,off     ,sigy    ,
     3         rho0    ,pla     ,dpla    ,soundsp ,et      ,seq     ,
     4         depsxx  ,depsyy  ,depszz  ,depsxy  ,depsyz  ,depszx  ,
     5         sigoxx  ,sigoyy  ,sigozz  ,sigoxy  ,sigoyz  ,sigozx  ,
     6         signxx  ,signyy  ,signzz  ,signxy  ,signyz  ,signzx  )
        CASE(2) 
          CALL mat115_newton(
     1         nel     ,ngl     ,nuparam ,nuvar   ,grho    , 
     2         time    ,timestep,uparam  ,uvar    ,off     ,sigy    ,
     3         rho0    ,pla     ,dpla    ,soundsp ,et      ,seq     ,
     4         depsxx  ,depsyy  ,depszz  ,depsxy  ,depsyz  ,depszx  ,
     5         sigoxx  ,sigoyy  ,sigozz  ,sigoxy  ,sigoyz  ,sigozx  ,
     6         signxx  ,signyy  ,signzz  ,signxy  ,signyz  ,signzx  )
c
      END SELECT         
c  
c-----------
C
      ! Counting broken elements
      IF ((cfail > zero).OR.(pfail > zero)) THEN
        nindx = 0
        DO i = 1,nel
          IF (off(i) == four_over_5) THEN
            nindx       = nindx+1
            indx(nindx) = i
          ENDIF
        ENDDO
        ! Printout element deletion
        IF (nindx > 0) THEN
          DO j=1,nindx
#include "lockon.inc"
            WRITE(iout, 1000) ngl(indx(j))
            WRITE(istdo,1100) ngl(indx(j)),tt
#include "lockoff.inc"
          ENDDO
        ENDIF
      ENDIF
!---
 1000 FORMAT(1x,'RUPTURE OF SOLID ELEMENT NUMBER ',i10)
 1100 FORMAT(1x,'RUPTURE OF SOLID ELEMENT NUMBER ',i10,
     .          ' AT TIME :',g11.4)                          
!---