hm_read_fail_tsaihill.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_tsaihill (fail, mat_id, fail_id, irupt, lsubmodel, unitab, fail_tag)

Function/Subroutine Documentation

hm_read_fail_tsaihill()

subroutine hm_read_fail_tsaihill	(	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,
		type (fail_tag_), intent(inout)	fail_tag )

Parameters

[in]	fail_id	failure model ID
[in]	mat_id	material law ID
[in]	irupt	failure model number
[in]	unitab	table of input units
[in]	lsubmodel	submodel table
[in,out]	fail	failure model data structure
[in,out]	fail_tag	failure model tag for buffer allocation

Definition at line 37 of file hm_read_fail_tsaihill.F.

C-----------------------------------------------
c    ROUTINE DESCRIPTION :
c    Read Tsai-Hill failure model parameters
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 
      USE elbuftag_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 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
      TYPE (FAIL_TAG_)    ,INTENT(INOUT) :: FAIL_TAG   !< failure model tag for buffer allocation
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: IFAIL_SH,IFAIL_SO,NUVAR
      my_real :: x11,x22,s12,tmax,fcut,pthkf
      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)
 
      ! 1st line of the card
      CALL hm_get_floatv('X11'      ,x11      ,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('X22'      ,x22      ,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('S12'      ,s12      ,is_available,lsubmodel,unitab)
      CALL hm_get_intv  ('ifail_sh' ,IFAIL_SH ,IS_AVAILABLE,LSUBMODEL)
      CALL HM_GET_INTV  ('ifail_so' ,IFAIL_SO ,IS_AVAILABLE,LSUBMODEL)
 
      ! 2nd line of the card
      CALL HM_GET_FLOATV('tau_max'  ,TMAX     ,IS_AVAILABLE,LSUBMODEL,UNITAB)
      CALL HM_GET_FLOATV('fcut'     ,FCUT     ,IS_AVAILABLE,LSUBMODEL,UNITAB)
      
      ! Check parameters and set defaults value
      IF (X11 == ZERO)  X11  = INFINITY
      IF (X22 == ZERO)  X22  = INFINITY
      IF (S12 == ZERO)  S12  = INFINITY
      IF (TMAX == ZERO) TMAX = INFINITY
      IFAIL_SH = MIN(IFAIL_SH,2)
      IFAIL_SH = MAX(0,IFAIL_SH)
      IFAIL_SO = MIN(IFAIL_SO,2)
      IFAIL_SO = MAX(0,IFAIL_SO)
      FCUT = MAX(ZERO,FCUT)  
c
      ! Shell element deletion flag
      IF (IFAIL_SH == 1) THEN
        PTHKF = EM06
      ELSEIF (IFAIL_SH == 2) THEN 
        PTHKF = ONE
      ENDIF
c---------------------------
      IF (FCUT > ZERO) THEN 
        NUVAR = 7
      ELSE
        NUVAR = 1
      ENDIF
c---------------------------
      FAIL%KEYWORD = 'tsai-hill' 
      FAIL%IRUPT   = IRUPT 
      FAIL%FAIL_ID = FAIL_ID 
      FAIL%NUPARAM = 7
      FAIL%NIPARAM = 0
      FAIL%NUVAR   = NUVAR
      FAIL%NFUNC   = 0
      FAIL%NTABLE  = 0
      FAIL%NMOD    = 1
      FAIL%PTHK    = PTHKF
c            
      ALLOCATE (FAIL%UPARAM(FAIL%NUPARAM))
      ALLOCATE (FAIL%IPARAM(FAIL%NIPARAM))
      ALLOCATE (FAIL%IFUNC (FAIL%NFUNC))
      ALLOCATE (FAIL%TABLE (FAIL%NTABLE))
c
      ! Modes of failure
      FAIL_TAG%LF_DAMMX = FAIL_TAG%LF_DAMMX + FAIL%NMOD
      ALLOCATE (FAIL%MODE(FAIL%NMOD))
      FAIL%MODE(1) = "Reserve factor"
c
      FAIL%NUPARAM   = 7
      FAIL%UPARAM(1) = X11
      FAIL%UPARAM(2) = X22
      FAIL%UPARAM(3) = S12
      FAIL%UPARAM(4) = TMAX
      FAIL%UPARAM(5) = FCUT  
      FAIL%UPARAM(6) = IFAIL_SH
      FAIL%UPARAM(7) = IFAIL_SO
c---------------------------
C
      ! Print parameter in 0.out file
      IF (IS_ENCRYPTED) THEN
        WRITE(IOUT,'(5x,a,//)')'confidential data'
      ELSE
        WRITE(IOUT,1000)
        WRITE(IOUT,1100) X11,X22,S12
.OR.        IF (IFAIL_SH > 0  IFAIL_SO > 0) THEN 
          WRITE(IOUT,1200) TMAX
        ENDIF
        IF (FCUT > ZERO) THEN 
          WRITE(IOUT,1300) FCUT 
        ENDIF 
        WRITE(IOUT,1400) IFAIL_SH
        WRITE(IOUT,1500) IFAIL_SO   
      ENDIF
c----------- 
 1000 FORMAT(
     & 5X,'  ----------------------------------------------------   ',/
     & 5X,'            failure criterion : tsai-hill                ',/,
     & 5X,'  ----------------------------------------------------   ',/)
 1100 FORMAT(
     & 5X,'critical strength in direction 1 x11  . . . . . . . . . =',1PG20.13/
     & 5X,'critical strength in direction 2 x22  . . . . . . . . . =',1PG20.13/
     & 5X,'critical shear strength s12 . . . . . . . . . . . . . . =',1PG20.13/)
 1200 FORMAT(
     & 5X,'stress softening activated:                              ',/,
     & 5X,'relaxation time tau_max . . . . . . . . . . . . . . . . =',1PG20.13/)
 1300 FORMAT(
     & 5X,'stress tensor filtering activated:                       ',/,
     & 5X,'filtering frequency fcut  . . . . . . . . . . . . . . . =',1PG20.13/)
 1400 FORMAT(
     & 5X,'shell deletion flag ifail_sh  . . . . . . . . . . . . . =',I10,/
     & 5X,'    = 0: shell never deleted and no stress softening      ',/,
     & 5X,'    = 1: shell deleted IF damage is reached for one layer ',/,
     & 5X,'    = 2: shell deleted IF damage is reached for all layers',/) 
 1500 FORMAT(
     & 5X,'solid deletion flag ifail_so  . . . . . . . . . . . . . =',I10,/
     & 5X,'    = 0: solid never deleted and no stress softening                  ',/,
     & 5X,'    = 1: solid deleted IF damage is reached for one integration point ',/,
     & 5X,'    = 2: solid deleted IF damage is reached for all integration points',/)   
c----------- 
      RETURN