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/)

 1300 FORMAT(

     & 5x,'UNIAXIAL DATA TEST CURVE. . . . . . . . .=',i10/

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

 1310 FORMAT(

     & 5x,'EQUIBIAXIAL DATA TEST CURVE . . . . . . .=',i10/

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

 1320 FORMAT(

     & 5x,'PLANAR DATA TEST CURVE. . . . . . . . . .=',i10/

     & 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:415

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

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:150

elbuftag_mod
Definition elbuftag_mod.F:227

message_mod
Definition message_mod.F:1249

names_and_titles_mod
Definition names_and_titles_mod.F:997

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

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

submodel_mod
Definition submodel_mod.F:712

unitab_mod
Definition unitab_mod.F:513

starter
program starter
Definition starter.F:39

elbuftag_mod::mlaw_tag_
Definition elbuftag_mod.F:232

submodel_mod::submodel_data
Definition submodel_mod.F:723