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,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 line
      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
      IF(.NOT. 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
      IF(.NOT. 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 ! -- VELOCITIES convected frameE (V_REPCVT)
      prop_tag(igtyp)%G_VR_REPCVT = 3  ! -- VELOCITIES convected frameE (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