hm_read_mat25.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/.
!||====================================================================
!||    hm_read_mat25            ../starter/source/materials/mat/mat025/hm_read_mat25.F
!||--- called by ------------------------------------------------------
!||    hm_read_mat              ../starter/source/materials/mat/hm_read_mat.F90
!||--- calls      -----------------------------------------------------
!||    hm_get_intv              ../starter/source/devtools/hm_reader/hm_get_intv.F
!||    init_mat_keyword         ../starter/source/materials/mat/init_mat_keyword.F
!||    read_mat25_crasurv       ../starter/source/materials/mat/mat025/read_mat25_crasurv.F90
!||    read_mat25_tsaiwu        ../starter/source/materials/mat/mat025/read_mat25_tsaiwu.F90
!||--- uses       -----------------------------------------------------
!||    elbuftag_mod             ../starter/share/modules1/elbuftag_mod.F
!||    message_mod              ../starter/share/message_module/message_mod.F
!||    read_mat25_crasurv_mod   ../starter/source/materials/mat/mat025/read_mat25_crasurv.F90
!||    read_mat25_tsaiwu_mod    ../starter/source/materials/mat/mat025/read_mat25_tsaiwu.F90
!||    submodel_mod             ../starter/share/modules1/submodel_mod.F
!||====================================================================
      SUBROUTINE hm_read_mat25(
     .           UPARAM   ,MAXUPARAM,NUPARAM  ,UNITAB   ,LSUBMODEL,
     .           MTAG     ,ID       ,TITR     ,PM       ,ISRATE   ,
     .           PARMAT   ,MATPARAM )
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
      END