hm_read_mat25.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	hm_read_mat25 (uparam, maxuparam, nuparam, unitab, lsubmodel, mtag, id, titr, pm, israte, parmat, matparam)

Function/Subroutine Documentation

hm_read_mat25()

subroutine hm_read_mat25	(	dimension(maxuparam)	uparam,
		integer	maxuparam,
		integer	nuparam,
		type (unit_type_), intent(in)	unitab,
		type(submodel_data), dimension(*), intent(in)	lsubmodel,
		type(mlaw_tag_), intent(inout)	mtag,
		integer, intent(in)	id,
		character(len=nchartitle), intent(in)	titr,
		intent(inout)	pm,
		integer, intent(inout)	israte,
		intent(inout)	parmat,
		type(matparam_struct_), intent(inout)	matparam )

Definition at line 39 of file hm_read_mat25.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE unitab_mod
      USE elbuftag_mod            
      USE message_mod
      USE submodel_mod
      USE matparam_def_mod
      USE read_mat25_crasurv_mod
      USE read_mat25_tsaiwu_mod
      USE names_and_titles_mod , ONLY : nchartitle
C-----------------------------------------------
C   ROUTINE DESCRIPTION :
C   ===================
C   READ MAT LAW25 WITH HM READER
C-----------------------------------------------
C   DUMMY ARGUMENTS DESCRIPTION:
C   ===================
C     UNITAB          UNITS ARRAY
C     MAT_ID              MATERIAL MAT_ID(INTEGER)
C     TITR            MATERIAL TITLE
C     LSUBMODEL       SUBMODEL STRUCTURE    
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"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB 
      INTEGER, INTENT(IN)                          :: ID
      TYPE(SUBMODEL_DATA), DIMENSION(*),INTENT(IN) :: LSUBMODEL
      CHARACTER(LEN=NCHARTITLE) ,INTENT(IN)             :: TITR
c
      INTEGER, INTENT(INOUT)                         :: ISRATE
      my_real, DIMENSION(100)       ,INTENT(INOUT)   :: parmat     
      my_real, DIMENSION(NPROPM)    ,INTENT(INOUT)   :: pm     
      TYPE(MLAW_TAG_)               ,INTENT(INOUT) :: MTAG
      TYPE(MATPARAM_STRUCT_) ,INTENT(INOUT) :: MATPARAM
c
      INTEGER :: MAXUPARAM,NUPARAM
      my_real, DIMENSION(MAXUPARAM)    :: uparam
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: ILAW,IFLAG
      LOGICAL :: IS_AVAILABLE
C=======================================================================
      is_available = .false.
      ilaw = 25
c--------------------------------------------------------
c     Formulation type
c--------------------------------------------------------
c     
      CALL hm_get_intv  ('MAT_Iflag' ,iflag    ,is_available,lsubmodel)
      
      IF (iflag == 0) then   ! Tsai-Wu formulation
        CALL read_mat25_tsaiwu(
     .           matparam ,parmat   ,unitab   ,lsubmodel,
     .           id       ,titr     ,pm       ,israte   ,
     .           iout     ,npropm   )
        matparam%NMOD = 3
      ELSE     
       CALL read_mat25_crasurv(
     .           matparam ,parmat   ,unitab   ,lsubmodel,
     .           id       ,titr     ,pm       ,israte   ,
     .           iout     ,npropm   )
        matparam%NMOD = 6
      END IF
c-------------------------------------------------------------
c---- Element buffer variable activation (for automatic allocation)
c-------------------------------------------------------------
c
      mtag%G_IERR   = 1  ! tag element for projection error message
      mtag%G_GAMA   = 6  
      mtag%G_PLA    = 1    
      mtag%G_EPSD   = 1    
      mtag%G_DAMDL  = 1  
      mtag%G_TSAIWU = 1 
      mtag%G_DMG    = 1 + matparam%NMOD 
c
      mtag%L_GAMA   = 6  
      mtag%L_PLA    = 1
      mtag%L_EPSD   = 1
      mtag%L_CRAK   = 2
      mtag%L_STRA   = 6
      mtag%L_TSAIWU = 1
      mtag%L_DMG    = 1 + matparam%NMOD
c
      ALLOCATE(matparam%MODE(matparam%NMOD))
      matparam%MODE(1) = "Tensile damage in direction 1"
      matparam%MODE(2) = "Tensile damage in direction 2"
      matparam%MODE(3) = "Global maximum plastic work"
      IF (iflag > 0) THEN 
        matparam%MODE(4) = "Failure index in direction 1"
        matparam%MODE(5) = "Failure index in direction 2"
        matparam%MODE(6) = "Failure index in shear plane 12"
      ENDIF
c
c-----------
      CALL init_mat_keyword(matparam,"HOOK")
      CALL init_mat_keyword(matparam,"ELASTO_PLASTIC")
      CALL init_mat_keyword(matparam,"COMPRESSIBLE")
      CALL init_mat_keyword(matparam,"ORTHOTROPIC")
c-----------
      ! Properties compatibility
      CALL init_mat_keyword(matparam,"SOLID_ORTHOTROPIC")
      IF (iflag == 0) THEN 
        CALL init_mat_keyword(matparam,"SHELL_ORTHOTROPIC")
      ELSE 
        CALL init_mat_keyword(matparam,"SHELL_COMPOSITE")    
      ENDIF
c--------
      RETURN