hm_read_prop46.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_prop46           ../starter/source/properties/spring/hm_read_prop46.F
!||--- called by ------------------------------------------------------
!||    hm_read_prop_generic     ../starter/source/properties/hm_read_prop_generic.F
!||--- calls      -----------------------------------------------------
!||    hm_get_floatv            ../starter/source/devtools/hm_reader/hm_get_floatv.F
!||    hm_get_floatv_dim        ../starter/source/devtools/hm_reader/hm_get_floatv_dim.F
!||    hm_get_intv              ../starter/source/devtools/hm_reader/hm_get_intv.F
!||    hm_option_is_encrypted   ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
!||    set_u_geo                ../starter/source/user_interface/uaccess.F
!||    set_u_pnu                ../starter/source/user_interface/uaccess.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_prop46(IOUT,NUVAR ,PARGEO, UNITAB,IGTYP,
     .                          ID,PROP_TAG,TITR,LSUBMODEL,IUNIT)
C-----------------------------------------------
      USE unitab_mod
      USE message_mod
      USE elbuftag_mod
      USE submodel_mod
      USE names_and_titles_mod , ONLY : nchartitle
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C----------+---------+---+---+--------------------------------------------
C VAR      | SIZE    |TYP| RW| DEFINITION
C----------+---------+---+---+--------------------------------------------
C IIN      |  1      | I | R | INPUT FILE UNIT (D00 file) 
C IOUT     |  1      | I | R | OUTPUT FILE UNIT (L00 file)
C NUVAR    |  1      | I | W | NUMBER OF USER ELEMENT VARIABLES
C----------+---------+---+---+--------------------------------------------
C PARGEO   |  *      | F | W | 1)SKEW NUMBER
C          |         |   |   | 2)STIFNESS FOR INTERFACE
C          |         |   |   | 3)FRONT WAVE OPTION
C          |         |   |   | 4)... not yet used
C----------+---------+---+---+--------------------------------------------
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "tablen_c.inc"
C----------+---------+---+---+--------------------------------------------
C
C     This subroutine read the user geometry parameters.
C
C     The geometry datas has to bee stored in radioss storage 
C     with the function SET_U_GEO(value_index,value). 
C
C     If some standard radioss functions (time function or 
C     x,y function) are used, this function IDs has to 
C     bee stored with the function SET_U_PNU(func_index,func_id,KFUNC).
C
C     If this property refers to a user material, this 
C     material IDs has to bee stored with the function 
C     SET_U_PNU(mat_index,mat_id,KMAT).
C
C     If this property refers to a user property, this 
C     sub-property IDs has to bee stored with the function 
C     SET_U_PNU(sub_prop_index,sub_prop_id,KMAT).
C
C     SET_U_GEO and SET_U_PNU return 0 if no error 
C     SET_U_GEO and SET_U_PNU return the maximum allowed index 
C     if index is larger than this maximum
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB 
      INTEGER IOUT,NUVAR,IGTYP,ID,IUNIT
      my_real pargeo(*)
      CHARACTER(LEN=NCHARTITLE) :: TITR
      TYPE(prop_tag_) , DIMENSION(0:MAXPROP) :: PROP_TAG
      TYPE(submodel_data),INTENT(IN)::LSUBMODEL(*)
C=======================================================================
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IFUNC1,IFUNC2,IFUNC3,IFUNC4,IFUNC5,IERROR,KFUNC,
     .   epsi,idens
      my_real
     .  amas,elastif,xlim1,xlim2,xk,damp,fac_m, fac_l, fac_t,fscale,
     .  scalet,scalex,scalev,scalef,r_epsi,scalex_unit,scalev_unit
      LOGICAL IS_AVAILABLE, IS_ENCRYPTED
C-----------------------------------------------
 
      INTEGER SET_U_PNU,SET_U_GEO
      EXTERNAL set_u_pnu,set_u_geo
      parameter(kfunc=29)
      !PARAMETER (KMAT=31)
      !PARAMETER (KPROP=30)
C             
C=======================================================================
C
      is_encrypted = .false.
      is_available = .false.
C
C--------------------------------------------------
C EXTRACT DATA (IS OPTION CRYPTED)
C--------------------------------------------------
      CALL hm_option_is_encrypted(is_encrypted)
C--------------------------------------------------
C EXTRACT DATAS (INTEGER VALUES)
C--------------------------------------------------
      CALL hm_get_intv('FUN_A1',ifunc1,is_available,lsubmodel)
      CALL hm_get_intv('FUN_B1',ifunc2,is_available,lsubmodel)
      CALL hm_get_intv('FUN_C1',ifunc3,is_available,lsubmodel)
      CALL hm_get_intv('FUN_D1',ifunc4,is_available,lsubmodel)
      CALL hm_get_intv('MAT_IMASS',idens,is_available,lsubmodel)
      CALL hm_get_intv('EPSI',epsi,is_available,lsubmodel)
C--------------------------------------------------
C EXTRACT DATAS (REAL VALUES)
C--------------------------------------------------
      CALL hm_get_floatv('PROP_MASS',amas,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('STIFF0',elastif,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('VEL_X',xlim1,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('NFORCE',xlim2,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('STIFF1',xk,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('DAMP1',damp,is_available,lsubmodel,unitab)
C-- Optional ligne
      CALL hm_get_floatv('FScale11',scalet,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('FScale22',scalex,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('FScale21',scalev,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('fscale12',SCALEF,IS_AVAILABLE,LSUBMODEL,UNITAB)
C
      CALL HM_GET_FLOATV_DIM('fscale22',SCALEX_UNIT,IS_AVAILABLE,LSUBMODEL,UNITAB)
      CALL HM_GET_FLOATV_DIM('fscale21',SCALEV_UNIT,IS_AVAILABLE,LSUBMODEL,UNITAB)
C----------------------
C
.NOT.      IF( IS_ENCRYPTED)THEN
        WRITE(IOUT,1400) ID
      ELSE
        WRITE(IOUT,1500) ID
      ENDIF
C
      FAC_M  = UNITAB%FAC_M(IUNIT)                        
      FAC_L  = UNITAB%FAC_L(IUNIT)
      FAC_T  = UNITAB%FAC_T(IUNIT)
      FSCALE = FAC_M * FAC_L / (FAC_T **TWO)
c        WRITE(IOUT,*) ' My User routine : Debug Start USER'
c        WRITE(IOUT,*) '        *** MAT MUSCLE ***'
      NUVAR     = 4
C
      !default contact stiff
      IF (XK == ZERO)THEN
         XK = ELASTIF
      ENDIF
      PARGEO(1) = ZERO
      PARGEO(2) = XK
C     front wave = 1
      PARGEO(3) = ONE
C
      IF (XLIM1 == ZERO) XLIM1 = EP30*SCALEV_UNIT
      IF (SCALET == ZERO) SCALET = ONE*FAC_T
      IF (SCALEX == ZERO) SCALEX = ONE*SCALEX_UNIT
      IF (SCALEV == ZERO) SCALEV = ONE*SCALEV_UNIT
      IF (SCALEF == ZERO) SCALEF = ONE*FSCALE
C
      IERROR = SET_U_GEO(1,AMAS)
      IERROR = SET_U_GEO(2,ELASTIF)
      IERROR = SET_U_GEO(3,XLIM1)
      IERROR = SET_U_GEO(4,XLIM2)
      IERROR = SET_U_GEO(5,XK)
      IERROR = SET_U_GEO(6,DAMP)
      R_EPSI = EPSI
      IERROR = SET_U_GEO(7,R_EPSI)
      IERROR = SET_U_GEO(8,SCALET)
      IERROR = SET_U_GEO(9,SCALEX)
      IERROR = SET_U_GEO(10,SCALEV)
      IERROR = SET_U_GEO(11,SCALEF)
      IERROR = SET_U_GEO(12,(REAL(IDENS)+EM05))
C
      IFUNC5 = 0
      IERROR = SET_U_PNU(1,IFUNC1,KFUNC)
      IERROR = SET_U_PNU(2,IFUNC2,KFUNC)
      IERROR = SET_U_PNU(3,IFUNC3,KFUNC)
      IERROR = SET_U_PNU(4,IFUNC4,KFUNC)
      IERROR = SET_U_PNU(5,IFUNC5,KFUNC)
C
.NOT.      IF( IS_ENCRYPTED)THEN
        WRITE(IOUT,1000)
     .     AMAS,ELASTIF,XLIM1,XLIM2,XK,DAMP,EPSI,
     .     IFUNC1,IFUNC2,IFUNC3,
     .     IFUNC4,IDENS,SCALEF,SCALET,SCALEX,SCALEV
      ENDIF
C
C-----------------------------
C       PROPERTY BUFFER 
C-----------------------------
 
      PROP_TAG(IGTYP)%G_FOR  = 3
      PROP_TAG(IGTYP)%G_MOM = 5
      PROP_TAG(IGTYP)%G_SKEW = 6
      PROP_TAG(IGTYP)%G_MASS = 1
      PROP_TAG(IGTYP)%G_V_REPCVT  = 3 ! -- VITESSES REPERE CONVECTEE (V_REPCVT)
      PROP_TAG(IGTYP)%G_VR_REPCVT = 3  ! -- VITESSES REPERE CONVECTEE (VR_REPCVT)
      PROP_TAG(IGTYP)%G_NUVAR = NUVAR
C
      RETURN
C
 1000 FORMAT(
     & 5X,'mass per element . . . . . . . . . . . . .=',E12.4/,
     & 5X,'stiffness per unit length. . . . . . . . .=',E12.4/,
     & 5X,'maximum strain rate. . . . . . . . . . . .=',E12.4/,
     & 5X,'maximum force. . . . . . . . . . . . . . .=',E12.4/,
     & 5X,'stiffness for interface. . . . . . . . . .=',E12.4/,
     & 5X,'damping VALUE. . . . . . . . . . . . . . .=',E12.4/,
     & 5X,'elongation calculation flag:1=l;0=l/l_o-1 =',I10/,
C     & 5X,' (1: elongation=L , 0(def): =(L / L_0) -1)',E12.4/,
     & 5X,'force vs. time active FUNCTION id. . . . .=',I10/,
     & 5X,'force vs. deflection active function id. .=',I10/,
     & 5X,'force vs. velocity active function id. . .=',I10/,
     & 5X,'force vs. deflection passive function id .=',I10/,
C     & 5X,'DEBUGGING FUNCTION ID. . . . . . . . . . .=',I10/,
     & 5X,'flag for mass input. . . . . . . . . . . .=',I10/,
     & 5X,'force scale factor . . . . . . . . . . . .=',E12.4/,
     & 5X,'time scale factor . . . . . . . . . . . .=',E12.4/,
     & 5X,'elongation scale factor . . . . . . . . .=',E12.4/,
     & 5X,'velocity scale factor . . . . . . . . . .=',E12.4//)
C
 1400 FORMAT(
     & 5X,'user property set'/,
     & 5X,'property set number . . . . . . . . . .=',I10)
C
 1500 FORMAT(
     & 5X,'user property set'/,
     & 5X,'property set number . . . . . . . . . .=',I10,
     & 5X,'confidential data'//)
C
      END