sigeps115.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    sigeps115       ../engine/source/materials/mat/mat115/sigeps115.F
!||--- called by ------------------------------------------------------
!||    mulaw           ../engine/source/materials/mat_share/mulaw.F90
!||--- calls      -----------------------------------------------------
!||    mat115_newton   ../engine/source/materials/mat/mat115/mat115_newton.F
!||    mat115_nice     ../engine/source/materials/mat/mat115/mat115_nice.F
!||====================================================================
      SUBROUTINE sigeps115(
     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-----------------------------------------------
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)                          
!---
      END