hm_read_fail_sahraei.F File Reference

#include "implicit_f.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	hm_read_fail_sahraei (fail, mat_id, fail_id, irupt, lsubmodel, unitab)

Function/Subroutine Documentation

hm_read_fail_sahraei()

subroutine hm_read_fail_sahraei	(	type(fail_param_), intent(inout)	fail,
		integer, intent(in)	mat_id,
		integer, intent(in)	fail_id,
		integer, intent(in)	irupt,
		type(submodel_data), dimension(*), intent(in)	lsubmodel,
		type(unit_type_), intent(in)	unitab )

Definition at line 37 of file hm_read_fail_sahraei.F.

C-----------------------------------------------
c    ROUTINE DESCRIPTION :
c    Sahraei failure criterion (irup = 29)
C-----------------------------------------------
C    M o d u l e s
C-----------------------------------------------
      USE fail_param_mod
      USE unitab_mod
      USE message_mod
      USE submodel_mod
      USE hm_option_read_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   B l o c k s
C-----------------------------------------------
#include      "units_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER            ,INTENT(IN) :: FAIL_ID         ! failure model ID
      INTEGER            ,INTENT(IN) :: MAT_ID          ! material law ID
      INTEGER            ,INTENT(IN) :: IRUPT           ! failure model type number
      TYPE(UNIT_TYPE_)   ,INTENT(IN) :: UNITAB          ! table of input units
      TYPE(SUBMODEL_DATA),INTENT(IN) :: LSUBMODEL(*)    ! submodel table
      TYPE(FAIL_PARAM_)  ,INTENT(INOUT) :: FAIL         ! failure model data structure
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: NUM,DEN,ORDI,COMP_DIR,IDEL,NUPARAM,NUVAR,NFUNC
      INTEGER ,PARAMETER :: NSIZE = 2
      INTEGER ,DIMENSION(NSIZE) :: IFUNC
      my_real :: vol_strain,max_comp_strain,ratio,el_ref,el_ref_unit
      LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
C-----------------------------------------------
C   S o u r c e   L i n e s
C=======================================================================
      is_encrypted = .false.
      is_available = .false.
 
      CALL hm_option_is_encrypted(is_encrypted)
c
c!#
c!# NUMERATOR    =  1  ==> Eps_xx 
c!# NUMERATOR    =  2  ==> Eps_yy 
c!# NUMERATOR    =  3  ==> Eps_zz 
c!# NUMERATOR    =  4  ==> Eps_1   1st. 3d-principal
c!# NUMERATOR    =  5  ==> Eps_2   2nd. 3d-principal
c!# NUMERATOR    =  6  ==> Eps_3   3rd. 3d-principal
c!#
c!# Denominator  =  1  ==> 2D e1 = x-z - plane 
c!# Denominator  =  2  ==> 2D e1 = x-y - plane 
c!# Denominator  =  3  ==> 2D e1 = y-z - plane 
c!# Denominator  =  4  ==> Eps_1   1st. 3d-principal
c!# Denominator  =  5  ==> Eps_2   2nd. 3d-principal 
c!# Denominator  =  6  ==> Eps_3   3rd. 3d-principal 
c!#
c!# Ordinate     =  1  ==> MAX(Eps_xx,Eps_yy,Eps_zz)
c!# Ordinate     =  2  ==> Eps_xx
c!# Ordinate     =  3  ==> Eps_yy
c!# Ordinate     =  4  ==> Eps_zz
c!# Ordinate     =  5  ==> Eps_1   1st. 3d-principal
c!# Ordinate     =  6  ==> 2D - e1 = x-z - plane
c!# Ordinate     =  7  ==> 2D - e1 = x-y - plane
c!# Ordinate     =  8  ==> 2D - e1 = y-z - plane
c!#
c!# COMP_DIR     = 0   ==> Feature is not taken into account
c!# COMP_DIR     = 1   ==> Normal to layers is in xx
c!# COMP_DIR     = 2   ==> Normal to layers is in yy
c!# COMP_DIR     = 3   ==> Normal to layers is in zz
c!#
c!# MAX_COMP_STRAIN = In-plane failure compression strain (not deleting, only damage=1)
c!# if value is negative, the element will NOT be deleted, only DAMAGE set to 1
c!# RATIO        = RATIO of the two other failure strains (default = 1.0)
c!# 
C===================================================================================
C READING CARD 1
C===================================================================================
      CALL hm_get_intv   ('Fct_Ratio'      ,ifunc(1)        ,is_available,lsubmodel) 
      CALL hm_get_intv   ('NUM'            ,num             ,is_available,lsubmodel) 
      CALL hm_get_intv   ('DEN'            ,den             ,is_available,lsubmodel)
      CALL hm_get_intv   ('ORDI'           ,ordi            ,is_available,lsubmodel)  
      IF (ordi == 0) ordi = 1
      CALL hm_get_floatv ('vol_strain'     ,VOL_STRAIN      ,IS_AVAILABLE,LSUBMODEL,UNITAB)
      CALL HM_GET_INTV   ('fct_elsize'     ,IFUNC(2)        ,IS_AVAILABLE,LSUBMODEL)
      CALL HM_GET_FLOATV ('el_ref'         ,EL_REF          ,IS_AVAILABLE,LSUBMODEL,UNITAB)
.AND.      IF ((EL_REF == ZERO)(IFUNC(2) > 0)) THEN 
        CALL HM_GET_FLOATV_DIM('el_ref' ,EL_REF_UNIT,IS_AVAILABLE, LSUBMODEL, UNITAB)
        EL_REF = ONE*EL_REF_UNIT
      ENDIF 
C===================================================================================
C READING CARD 2 
C===================================================================================    
      CALL HM_GET_INTV   ('comp_dir'       ,COMP_DIR        ,IS_AVAILABLE,LSUBMODEL) 
      CALL HM_GET_INTV   ('idel'           ,IDEL            ,IS_AVAILABLE,LSUBMODEL) 
      CALL HM_GET_FLOATV ('max_comp_strain',MAX_COMP_STRAIN ,IS_AVAILABLE,LSUBMODEL,UNITAB)
      IF (MAX_COMP_STRAIN == ZERO) MAX_COMP_STRAIN = INFINITY
      CALL HM_GET_FLOATV ('ratio'          ,RATIO           ,IS_AVAILABLE,LSUBMODEL,UNITAB)
      IF (RATIO == ZERO) RATIO = ONE
C===================================================================================
C NUMBER OF PARAMETERS,FUNCTIONS, AND USER VARIABLE
C===================================================================================
      NUPARAM = 9
      IF (IFUNC(2) > 0) THEN
        NFUNC = 2
        NUVAR = 1
      ELSE
        NFUNC = 1
        NUVAR = 0
      ENDIF
C===================================================================================
C STORING PARAMETERS
C===================================================================================
      FAIL%KEYWORD = 'sahraei' 
      FAIL%IRUPT   = IRUPT 
      FAIL%FAIL_ID = FAIL_ID 
      FAIL%NUPARAM = NUPARAM
      FAIL%NIPARAM = 0
      FAIL%NUVAR   = NUVAR
      FAIL%NFUNC   = NFUNC
      FAIL%NTABLE  = 0
      FAIL%NMOD    = 0
c            
      ALLOCATE (FAIL%UPARAM(FAIL%NUPARAM))
      ALLOCATE (FAIL%IPARAM(FAIL%NIPARAM))
      ALLOCATE (FAIL%IFUNC (FAIL%NFUNC))
      ALLOCATE (FAIL%TABLE (FAIL%NTABLE))
      
      FAIL%IFUNC(1:NFUNC) = IFUNC(1:NFUNC)
c
      FAIL%UPARAM(1) = VOL_STRAIN
      FAIL%UPARAM(2) = NUM
      FAIL%UPARAM(3) = DEN
      FAIL%UPARAM(4) = ORDI
      IF (COMP_DIR == ZERO) THEN
        FAIL%UPARAM(5) = ZERO
      ELSE
        FAIL%UPARAM(5) = COMP_DIR
      ENDIF
      FAIL%UPARAM(6) = MAX_COMP_STRAIN
      FAIL%UPARAM(7) = RATIO
      FAIL%UPARAM(8) = IDEL
      FAIL%UPARAM(9) = EL_REF
C===================================================================================
C PRINTING OUT DATA
C===================================================================================
      IF (IS_ENCRYPTED) THEN
        WRITE (IOUT,'(5x,a,//)') 'confidential data'
      ELSE 
        WRITE (IOUT, 1000) IFUNC(1),VOL_STRAIN
        WRITE (IOUT, 2000) NUM,DEN,ORDI
        IF (IFUNC(2) > 0) THEN 
          WRITE (IOUT, 3000) IFUNC(2),EL_REF
        ENDIF
        IF (COMP_DIR /= ZERO) THEN
          WRITE (IOUT, 4000) COMP_DIR, IDEL, MAX_COMP_STRAIN, RATIO
        ENDIF
      ENDIF
C===================================================================================
C PRINTOUT FORMAT
C===================================================================================
 1000 FORMAT(
     & 5X,'    sahraei local electric battery failure      ',/,
     & 5X,'    --------------------------------------      ',/,
     & 5X,'strain ratio FUNCTION id . . . . . . . . . . . . .=',I10/,
     & 5X,'volumetric strain limit. . . . . . . . . . . . . .=',1PG20.13/)
 2000 FORMAT(
     & 5X,'numerator flag for strain ratio. . . . . . . . . .=',I10/,
     & 5X,'denominator flag for strain ratio. . . . . . . . .=',I10/,
     & 5X,'failure ordinate for failure limit . . . . . . . .=',I10/)
 3000 FORMAT(
     & 5X,'element size regularization function id. . . . . .=',I10/,
     & 5X,'element reference size . . . . . . . . . . . . . .=',1PG20.13/)
 4000 FORMAT(
     & 5X,'normal in-plane compression direction. . . . . . .=',I10/,
     & 5X,'flag for element deletion in compression . . . . .=',I10/,
     & 5X,'in-plane compression failure strain  . . . . . . .=',1PG20.13/,
     & 5X,'ratio for in-plane 2nd direction failure strain . =',1PG20.13/)
C-----------
      RETURN