read_material_models.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/.
!||====================================================================
!||    read_material_models       ../starter/source/materials/read_material_models.F
!||--- called by ------------------------------------------------------
!||    lectur                     ../starter/source/starter/lectur.F
!||--- calls      -----------------------------------------------------
!||    fill_buffer_51_0           ../starter/source/materials/mat/mat051/fill_buffer_51_0.F
!||    hm_read_eos                ../starter/source/materials/eos/hm_read_eos.F
!||    hm_read_fail               ../starter/source/materials/fail/hm_read_fail.F
!||    hm_read_leak               ../starter/source/airbag/hm_read_leak.F
!||    hm_read_mat                ../starter/source/materials/mat/hm_read_mat.F90
!||    hm_read_nonlocal           ../starter/source/materials/nonlocal/hm_read_nonlocal.F
!||    hm_read_therm              ../starter/source/materials/therm/hm_read_therm.F
!||    hm_read_therm_stress       ../starter/source/materials/therm/hm_read_therm_stress.F90
!||    hm_read_visc               ../starter/source/materials/visc/hm_read_visc.F
!||    ini_eos_vars               ../starter/source/elements/elbuf_init/ini_eos_vars.F
!||    ini_mat_elem               ../starter/source/materials/mat/ini_mat_elem.F
!||    ini_mlaw_vars              ../starter/source/elements/elbuf_init/ini_mlaw_vars.F
!||    multi_check_eos            ../starter/source/multifluid/multi_check_eos.F
!||    multi_check_psh            ../starter/source/multifluid/multi_check_psh.F
!||    read_ale_mat               ../starter/source/materials/ale/read_ale_mat.F
!||    read_euler_mat             ../starter/source/materials/ale/read_euler_mat.f
!||--- uses       -----------------------------------------------------
!||    elbuftag_mod               ../starter/share/modules1/elbuftag_mod.F
!||    hm_option_read_mod         ../starter/share/modules1/hm_option_read_mod.F
!||    hm_read_mat_mod            ../starter/source/materials/mat/hm_read_mat.F90
!||    hm_read_therm_stress_mod   ../starter/source/materials/therm/hm_read_therm_stress.F90
!||    message_mod                ../starter/share/message_module/message_mod.F
!||    submodel_mod               ../starter/share/modules1/submodel_mod.F
!||    table_mod                  ../starter/share/modules1/table_mod.F
!||====================================================================
      SUBROUTINE read_material_models(
     .             MAT_ELEM    ,MLAW_TAG ,FAIL_TAG ,EOS_TAG  ,
     .             BUFMAT      ,SBUFMAT  ,IPM      ,PM       ,UNITAB   ,
     .             MULTI_FVM   ,FAILWAVE ,NLOC_DMG ,LSUBMODEL,TABLE    ,
     .             LTITR       ,USERL_AVAIL,MAT_NUMBER,
     .             NPC         ,TF       ,SNPC     ,NPTS    ,BUFLEN)
C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C This subroutine is managing all readers related to material model :
C    - Constitutive Laws  : /MAT
C    - Equations of State : /EOS
C    - Failure models     : /FAIL
C    - viscosity          : /VISC
C    - thermal parameters : /HEAT
C    - etc ...
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE unitab_mod
      USE elbuftag_mod            
      USE message_mod
      USE multi_fvm_mod
      USE failwave_mod
      USE nlocal_reg_mod
      USE submodel_mod
      USE hm_option_read_mod
      USE table_mod
      USE mat_elem_mod
      USE hm_read_mat_mod
      use hm_read_therm_stress_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      "param_c.inc"
#include      "tablen_c.inc"
#include      "units_c.inc"
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "scr19_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ,INTENT(IN)    :: LTITR
      INTEGER ,INTENT(IN)    :: USERL_AVAIL
      INTEGER ,INTENT(IN)    :: SBUFMAT
      INTEGER ,INTENT(INOUT) :: MAT_NUMBER
      INTEGER ,INTENT(INOUT) :: BUFLEN
      INTEGER ,DIMENSION(NPROPMI,NUMMAT), INTENT(INOUT) :: IPM
      my_real ,DIMENSION(NPROPM ,NUMMAT), INTENT(INOUT) :: pm
      my_real ,DIMENSION(SBUFMAT)       , INTENT(INOUT) :: bufmat
      TYPE (UNIT_TYPE_),INTENT(IN) :: UNITAB 
      TYPE(mlaw_tag_) , DIMENSION(NUMMAT)  , INTENT(INOUT) :: MLAW_TAG
      TYPE(eos_tag_)  ,DIMENSION(0:MAXEOS) , INTENT(INOUT) :: EOS_TAG
      TYPE(fail_tag_) ,DIMENSION(0:MAXFAIL), INTENT(INOUT) :: FAIL_TAG
      TYPE(multi_fvm_struct) :: MULTI_FVM
      TYPE(failwave_str_) :: FAILWAVE 
      TYPE(nlocal_str_) :: NLOC_DMG 
      TYPE(submodel_data) ,INTENT(IN)  :: LSUBMODEL(NSUBMOD)
      TYPE(ttable)  ,DIMENSION(NTABLE) :: TABLE 
      TYPE(mat_elem_) ,INTENT(INOUT)   :: MAT_ELEM
      INTEGER,INTENT(IN) :: SNPC, NPTS
      INTEGER,INTENT(IN) :: NPC(SNPC)
      my_real,INTENT(IN) :: tf(npts)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: IADBUF,ILAW
      CHARACTER MESS*40
C-----------------------------------------------
      DATA IADBUF /1/
      DATA MESS/'MATERIAL DEFINITION                     '/
C-----------------------------------------------
C   S o u r c e   L i n e s
C=======================================================================
      CALL ini_mat_elem (mat_elem)
    
      !Set default parameters used in material laws
      CALL ini_mlaw_vars(mlaw_tag ,nummat)   
      CALL ini_eos_vars (eos_tag)
c----------      
      ilaw    = 0
      buflen  = 0
      ilag    = 0
      iale    = 0
      ieuler  = 0
      iturb   = 0
 
      pm(1:100,nummat)=zero
      multi_fvm%NBMAT   = 0 !Initialize number of material in law151  
      failwave%WAVE_MOD = 0 !Initialization of failwave option
      failwave%NDDL = 0
c----------       
      !PRECONDITION - RETURN
      IF (nummat == 1) THEN
        ilag=1
        RETURN
      ENDIF
c---------- 
    
      WRITE(iout,1000)nummat-1
 
c-------------------------------------------------------------------
c     /MAT : Material laws (Constitutive law : sigma,epsilon)
c-------------------------------------------------------------------      
      CALL hm_read_mat(
     .     mat_elem%MAT_PARAM       ,mlaw_tag    ,eos_tag     ,bufmat      ,
     .     buflen      ,iadbuf      ,ipm         ,pm          ,
     .     multi_fvm   ,unitab      ,lsubmodel   ,table       ,
     .     sbufmat     ,npropmi     ,npropm      ,trimat      ,
     .     ialelag     ,ntable      ,nummat      ,hm_nummat   ,
     .     ltitr       ,userl_avail,mat_number   )
 
c-------------------------------------------------------------------
c     /EOS : Equations of State P=P(µ,E)
c-------------------------------------------------------------------
      CALL hm_read_eos(mat_elem%MAT_PARAM ,ipm      ,pm       ,bufmat   ,nummat, nfunct,
     .                 buflen   ,iadbuf   ,eos_tag  ,unitab   ,lsubmodel,
     .                 mlaw_tag ,npc      ,tf       ,snpc     ,npts     ,sbufmat,
     .                 ntable   ,table)
 
c-------------------------------------------------------------------
c     /FAIL : Failure Models
c-------------------------------------------------------------------     
      CALL hm_read_fail(mat_elem%MAT_PARAM,nummat,maxfail  ,fail_tag,
     .                  ntable  ,table    ,failwave ,nloc_dmg,
     .                  unitab  ,lsubmodel)
     
c-------------------------------------------------------------------
c     /VISC : Visco elastic  Models
c-------------------------------------------------------------------    
      CALL hm_read_visc(mat_elem%MAT_PARAM ,ipm     ,bufmat, 
     .                  unitab  ,lsubmodel ,table   )
 
c-------------------------------------------------------------------
c     /LEAK : Fabric Leakage Models
c-------------------------------------------------------------------
      CALL hm_read_leak(ipm, pm, unitab, lsubmodel)
 
c-------------------------------------------------------------------
C     /ALE/MAT
c-------------------------------------------------------------------
        CALL read_ale_mat(lsubmodel, unitab, ipm, pm)
 
c-------------------------------------------------------------------
C     /EULER/MAT
c-------------------------------------------------------------------
        CALL read_euler_mat(lsubmodel, unitab, ipm, pm)
        
c-------------------------------------------------------------------
c     /MAT/LAW51 - FILL BUFFER (modern input)
c-------------------------------------------------------------------
      CALL fill_buffer_51_0(mat_elem%MAT_PARAM, ipm, pm, bufmat, mat_elem%MAT_PARAM, mlaw_tag )
 
c-------------------------------------------------------------------
c     /HEAT/MAT : thermal parameters
c-------------------------------------------------------------------
      CALL hm_read_therm(mat_elem%MAT_PARAM,bufmat, buflen, ipm, pm, unitab, lsubmodel)
 
c-------------------------------------------------------------------
c     /THERM/STRESS : thermal expansion
c-------------------------------------------------------------------
      CALL hm_read_therm_stress(nummat   ,mat_elem%MAT_PARAM,mlaw_tag, unitab, lsubmodel, 
     .                          iout     ,npropm   ,npropmi   ,ipm      ,pm     )
 
c-------------------------------------------------------------------
c     /NONLOCAL : Non-local regularization
c------------------------------------------------------------------
      CALL hm_read_nonlocal(mat_elem%MAT_PARAM,nloc_dmg ,mlaw_tag ,ipm, unitab, lsubmodel)
 
c-------------------------------------------------------------------
c      ALE EULER SPECIFIC TREATMENTS 
c-------------------------------------------------------------------
      
      !SPECIFIC TREATMENTS FOR LAW151 MATERIAL LAWS RELATED TO PSH PARAMETERS & EOS
      IF (multi_fvm%IS_USED) THEN
         !Check pressure shift consistency between submaterials
         !They must be equal one to another, otherwise, throw an error
         CALL multi_check_psh(multi_fvm, nummat, npropmi, npropm, ipm, pm)
         CALL multi_check_eos(multi_fvm, nummat, npropmi, npropm, ipm, pm)
      ENDIF   
 
C------------------------------
      RETURN
c------------------------------
 1000 FORMAT(//
     .    '    MATERIAL LAWS'/
     .    '    -------------'/,
     .    '    NUMBER OF MATERIALS. . . . . . . . . =',i10//)
c------------------------------
      END