hm__read__mat94_8F_source.html

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_mat94            ../starter/source/materials/mat/mat094/hm_read_mat94.F

!||--- called by ------------------------------------------------------

!||    hm_read_mat              ../starter/source/materials/mat/hm_read_mat.F90

!||--- calls      -----------------------------------------------------

!||    hm_get_floatv            ../starter/source/devtools/hm_reader/hm_get_floatv.F

!||    hm_option_is_encrypted   ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F

!||    init_mat_keyword         ../starter/source/materials/mat/init_mat_keyword.F

!||--- uses       -----------------------------------------------------

!||    elbuftag_mod             ../starter/share/modules1/elbuftag_mod.F

!||    message_mod              ../starter/share/message_module/message_mod.f

!||    submodel_mod             ../starter/share/modules1/submodel_mod.F

!||====================================================================


      SUBROUTINE hm_read_mat94(UPARAM ,MAXUPARAM,NUPARAM  ,ISRATE   , IMATVIS  ,

     .                         NUVAR  ,IFUNC    ,MAXFUNC  ,NFUNC    , PARMAT   ,

     .                         UNITAB ,MAT_ID   ,TITR     ,MTAG     , LSUBMODEL,

     .                         PM     ,IPM      ,MATPARAM )

C-----------------------------------------------

C   D e s c r i p t i o n

C-----------------------------------------------

C   READ MAT LAW92 WITH HM READER ( TO BE COMPLETED )

C

C   DUMMY ARGUMENTS DESCRIPTION:

C   ===================

C

C     NAME            DESCRIPTION

C

C     PM              MATERIAL ARRAY(REAL)

C     UNITAB          UNITS ARRAY

C     ID              MATERIAL ID(INTEGER)

C     TITR            MATERIAL TITLE

C     LSUBMODEL       SUBMODEL STRUCTURE

C

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 names_and_titles_mod , ONLY : nchartitle, ncharfield

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

      my_real, INTENT(INOUT)                :: PM(NPROPM),PARMAT(100),UPARAM(MAXUPARAM)

      INTEGER, INTENT(INOUT)                :: IPM(NPROPMI),ISRATE,IFUNC(MAXFUNC),NFUNC,MAXFUNC,MAXUPARAM,NUPARAM, NUVAR,IMATVIS

      TYPE(mlaw_tag_),INTENT(INOUT)         :: MTAG

      INTEGER,INTENT(IN)                    :: MAT_ID

      CHARACTER(LEN=NCHARTITLE) ,INTENT(IN)       :: TITR

      TYPE(submodel_data),INTENT(IN)        :: LSUBMODEL(*)

      TYPE(matparam_struct_) ,INTENT(INOUT) :: MATPARAM

C-----------------------------------------------

C   L o c a l   V a r i a b l e s

C-----------------------------------------------

      INTEGER ID,IDAM,ILAW

      my_real E,NU,G,MU,RBULK,C10,C20,C30,D1,D2,D3

      CHARACTER(LEN=NCHARFIELD) ::  STRING

      my_real :: rho0, rhor

      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.

      idam = 0

      israte = 0

      imatvis = 1

      nuvar = 0

      ilaw = 94


      CALL hm_option_is_encrypted(is_encrypted)

      !line+1

      CALL hm_get_floatv('MAT_RHO'      ,rho0     ,is_available, lsubmodel, unitab)

      !! CALL HM_GET_FLOATV('MAT_RHO'      ,RHOR     ,IS_AVAILABLE, LSUBMODEL, UNITAB)

       !line-1 Is not active (Blank)

      !CALL HM_GET_INTV  ('Itype'        ,ITEST    ,IS_AVAILABLE, LSUBMODEL)

      !CALL HM_GET_INTV  ('MAT_FCT_IDI'  ,IFUNC(1) ,IS_AVAILABLE, LSUBMODEL)

      !CALL HM_GET_FLOATV('MAT_NU'       ,NU       ,IS_AVAILABLE, LSUBMODEL, UNITAB)

      !line-2

      CALL hm_get_floatv('LAW94_C01'     ,c10      ,is_available, lsubmodel, unitab)

      CALL hm_get_floatv('LAW94_C02'     ,c20      ,is_available, lsubmodel, unitab)

      CALL hm_get_floatv('LAW94_C03'     ,c30      ,is_available, lsubmodel, unitab)

      !line-3

      CALL hm_get_floatv('LAW94_D1'     ,d1      ,is_available, lsubmodel, unitab)

      CALL hm_get_floatv('LAW94_D2'     ,d2      ,is_available, lsubmodel, unitab)

      CALL hm_get_floatv('LAW94_D3'     ,d3      ,is_available, lsubmodel, unitab)

      !========== DEFAULT VALUES=============!

      rhor=zero

      !!IF(RHOR==ZERO)RHOR=RHO0

      pm(1) =rhor

      pm(89)=rho0


C

C Uniaxial - itest = 1

C equibiaxial - Itest = 2

c Planar - Itest = 3

C

!the proposed one      READ(LINE,FMT=FMT_A_I_F)STRING,IFUNC(1),NU

CSTRING = UNIAXIAL

CSTRING = EQUIBIAXIAL

CSTRING = PLANAR


C

      g = zero

      rbulk = zero

      e  = zero

C I

      IF(d2 /=  zero ) d2 = one/d2

      IF(d3 /= zero )  d3 = one/d3

C

      g = two*c10

      IF(d1 == zero) THEN

        d2 = zero

        d3 = zero

        nu = 0.495

        rbulk = two_third*g*(one + nu)/(one-two*nu)

        d1 = rbulk / two ! 1/D1

        e = two*g*(one + nu)

      ELSE

        d1 = one/d1

        rbulk= two*d1

        nu = (three*rbulk -two*g)/(three*rbulk + g)/two

        e = nine*rbulk*g/(three*rbulk + g)

      ENDIF

C

      uparam(1) = g

      uparam(2) = rbulk

      uparam(3) = nu

      uparam(4) = c10

      uparam(5) = c20

      uparam(6) = c30

      uparam(7) = d1

      uparam(8) = d2

      uparam(9) = d3

!! not used     UPARAM(10) = ITEST

C

      nuparam = 9

      IF(nuparam>maxuparam)THEN

        WRITE(iout,*)' ** ERROR : NUPARAM GT MAXUPARAM'

        WRITE(iout,*)'      NUPARAM =',nuparam,

     .                ' MAXUPARAM =',maxuparam

      ENDIF

C

      parmat(1) = rbulk

      parmat(2) = e

      parmat(3) = nu

C

C     Formulation for solid elements time step computation.

      parmat(16) = 2

      parmat(17) = two*g/(rbulk + four_over_3*g)

C

       imatvis = 1

       pm(100) = rbulk

c-----------------

      CALL init_mat_keyword(matparam,"INCOMPRESSIBLE")

      CALL init_mat_keyword(matparam,"TOTAL")

      CALL init_mat_keyword(matparam,"HOOK")

      ! Properties compatibility

      CALL init_mat_keyword(matparam,"SOLID_ISOTROPIC")

c-----------------

      !========== OUTPUT =============!


      WRITE(iout,1010) trim(titr),mat_id,94

      WRITE(iout,1000)

      IF(is_encrypted)THEN

        WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'

      ELSE

        WRITE(iout,1020)rho0

! not used

!!        IF (IFUNC(1) > 0) THEN

!!           IF(ITEST == 1) THEN

!!              WRITE(IOUT,1300)IFUNC(1), NU

!!           ELSEIF(ITEST == 2) THEN

!!              WRITE(IOUT,1310)IFUNC(1), NU

!!           ELSEIF(ITEST == 3) THEN

!!              WRITE(IOUT,1320)IFUNC(1), NU

!!           ENDIF

!!        ELSE

          WRITE(iout,1100)c10,c20,c30,d1,d2,d3

          WRITE(iout,1200)g,rbulk, nu

!!        ENDIF

      ENDIF

C

 1000 FORMAT(

     & 5x,'  YEOH LAW ',/,

     & 5x,'  ---------- ',//)

 1010 FORMAT(/

     & 5x,a,/,

     & 5x,'MATERIAL NUMBER. . . . . . . . . . . . . =',i10/,

     & 5x,'MATERIAL LAW . . . . . . . . . . . . . . =',i10/)

 1020 FORMAT(

     & 5x,'INITIAL DENSITY . . . . . . . . . . . . .=',1pg20.13/)

 1100 FORMAT(

     & 5x,'C10 . . . . . . . . . . . . . . . . . . .=',1pg20.13/

     & 5x,'C20 . . . . . . . . . . . . . . . . . . .=',1pg20.13/

     & 5x,'C30 . . . . . . . . . . . . . . . . . . .=',1pg20.13/

     & 5x,'1/D1  . . . . . . . . . . . . . . . . . .=',1pg20.13/

     & 5x,'1/D2  . . . . . . . . . . . . . . . . . .=',1pg20.13/

     & 5x,'1/D3  . . . . . . . . . . . . . . . . . .=',1pg20.13)

 1200 FORMAT(

     & 5x,'INITIAL SHEAR MODULUS . . . . . . . . . .=',1pg20.13/

     & 5x,'BULK MODULUS. . . . . . . . . . . . . . .=',1pg20.13/

     & 5x,'POISSON RATIO . . . . . . . . . . . . . .=',1pg20.13/)

C1300 FORMAT(

C    & 5X,'UNIAXIAL DATA TEST CURVE. . . . . . . . .=',I10/

C    & 5X,'POISSON RATIO . . . . . . . . . . . . . .=',1PG20.13///)

C1310 FORMAT(

C    & 5X,'EQUIBIAXIAL DATA TEST CURVE . . . . . . .=',I10/

C    & 5X,'POISSON RATIO . . . . . . . . . . . . . .=',1PG20.13///)

C 1320 FORMAT(

C     & 5X,'PLANAR DATA TEST CURVE. . . . . . . . . .=',I10/

C     & 5X,'POISSON RATIO . . . . . . . . . . . . . .=',1PG20.13///)

C

      RETURN


      END

C-----------------------------------------------

hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:418

hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:288

hm_read_mat94
subroutine hm_read_mat94(uparam, maxuparam, nuparam, israte, imatvis, nuvar, ifunc, maxfunc, nfunc, parmat, unitab, mat_id, titr, mtag, lsubmodel, pm, ipm, matparam)
Definition hm_read_mat94.F:40

init_mat_keyword
subroutine init_mat_keyword(matparam, keyword)
Definition init_mat_keyword.F:153

elbuftag_mod
Definition elbuftag_mod.F:230

message_mod
Definition message_mod.F:1257

names_and_titles_mod
Definition names_and_titles_mod.F:1002

names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1008

names_and_titles_mod::ncharfield
integer, parameter ncharfield
Definition names_and_titles_mod.F:1010

submodel_mod
Definition submodel_mod.F:717

unitab_mod
Definition unitab_mod.F:517

starter
program starter
Definition starter.F:39

elbuftag_mod::mlaw_tag_
Definition elbuftag_mod.F:235

submodel_mod::submodel_data
Definition submodel_mod.F:728