hm_read_inter_hertz_type17.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_inter_hertz_type17   ../starter/source/interfaces/int17/hm_read_inter_hertz_type17.F
!||--- called by ------------------------------------------------------
!||    hm_read_inter_struct         ../starter/source/interfaces/reader/hm_read_inter_struct.F
!||--- calls      -----------------------------------------------------
!||    hm_get_floatv                ../starter/source/devtools/hm_reader/hm_get_floatv.F
!||    hm_get_intv                  ../starter/source/devtools/hm_reader/hm_get_intv.F
!||    ngr2usr                      ../starter/source/system/nintrr.F
!||--- uses       -----------------------------------------------------
!||    message_mod                  ../starter/share/message_module/message_mod.F
!||    submodel_mod                 ../starter/share/modules1/submodel_mod.F
!||====================================================================
      SUBROUTINE hm_read_inter_hertz_type17(
     1        IPARI       ,FRIGAP       ,NOINT    ,NTYP     ,
     2        IGRBRIC     ,UNITAB       ,LSUBMODEL,TITR     )
C============================================================================
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE groupdef_mod
      USE submodel_mod
      USE unitab_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   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr06_c.inc"
#include      "com04_c.inc"
#include      "units_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NOINT,NI,NTYP
      INTEGER IPARI(*)
      my_real
     .   frigap(*)
      CHARACTER(LEN=NCHARTITLE)::TITR
      TYPE(SUBMODEL_DATA), DIMENSION(NSUBMOD), INTENT(IN) :: LSUBMODEL
      TYPE (UNIT_TYPE_), INTENT(IN) :: UNITAB 
      TYPE (GROUP_)  ,TARGET, DIMENSION(NGRBRIC) :: IGRBRIC
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER 
     .   IS1,IS2,ISU1,ISU2,ISU10,ISU20,MULTIMP,ITIED,UID
      my_real
     .   startt,bumult,stopt,fric,gap
      LOGICAL IS_AVAILABLE
      INTEGER, DIMENSION(:), POINTER :: INGR2USR
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER NGR2USR
C=======================================================================
C     READING INTERFACE /INTER/HERTZ/TYPE17
C=======================================================================
C Initializations of reals and integers
      is1     = 0
      is2     = 0
      startt  = zero
      stopt   = infinity
      fric    = zero
      gap     = zero
      bumult  = bmul0
      multimp = 24
      itied   = 0
      is_available = .false.
C------------------------------------------------------------
!  Printing out titles and interface ID
      WRITE(iout,1000) noint,trim(titr)
C------------------------------------------------------------
!  Reading cards
C------------------------------------------------------------
C  Card1
C------------------------------------------------------------
      CALL hm_get_intv('secondaryentityids'  ,isu10 ,is_available, lsubmodel)  
      CALL hm_get_intv('mainentityids' ,isu20 ,is_available, lsubmodel)  
C------------------------------------------------------------
C  Card2
C------------------------------------------------------------
      CALL hm_get_floatv('FRIC',fric ,is_available, lsubmodel, unitab)
C------------------------------------------------------------
!     Checking values
      is1 = 5
      is2 = 5
      ingr2usr => igrbric(1:ngrbric)%ID
      isu1=ngr2usr(isu10,ingr2usr,ngrbric)
      isu2=ngr2usr(isu20,ingr2usr,ngrbric)
C------------------------------------------------------------
! Storage in IPARI and FRIGAP
C------------------------------------------------------------
      ipari(7)  = ntyp
      ipari(13) = is1*10+is2
      ipari(15) = noint
      ipari(23) = multimp
      ipari(30) = itied
      ipari(45) = isu1
      ipari(46) = isu2
C
      frigap(1)  = fric
      frigap(2)  = gap
      frigap(3)  = startt
      frigap(4)  = bumult
      frigap(5)  = one
      frigap(11) = stopt
C------------------------------------------------------------
C     PRINTOUT
C------------------------------------------------------------
C
       WRITE(iout,1717) isu10,isu20,fric,startt,stopt
C
      IF(is1==0)THEN
        WRITE(iout,'(6X,A)')'NO SECONDARY SURFACE INPUT'
      ELSEIF(is1==1)THEN
        WRITE(iout,'(6X,A)')'SECONDARY SURFACE INPUT BY SEGMENTS'
      ELSEIF(is1==2)THEN
        WRITE(iout,'(6X,A)')'SECONDARY SURFACE INPUT BY NODES'
      ELSEIF(is1==3)THEN
        WRITE(iout,'(6X,A)')'SECONDARY SURFACE INPUT BY SEGMENTS'
      ELSEIF(is1==4 )THEN
        WRITE(iout,'(6X,A)')'SECONDARY SIDE INPUT BY BRICKS'
      ELSEIF(is1==5 )THEN
        WRITE(iout,'(6X,A)')'SECONDARY SIDE INPUT BY SOLID ELEMENTS'        
      ENDIF
      IF(is2==0)THEN
        WRITE(iout,'(6X,A)')'NO MAIN SURFACE INPUT'
      ELSEIF(is2==1)THEN
        WRITE(iout,'(6X,A)')'MAIN SURFACE INPUT BY SEGMENTS'
      ELSEIF(is2==2)THEN
        WRITE(iout,'(6X,A)')'MAIN SURFACE INPUT BY NODES'
      ELSEIF(is2==3)THEN
        WRITE(iout,'(6X,A)')'MAIN SURFACE INPUT BY SEGMENTS'
      ELSEIF(is2==4)THEN
        WRITE(iout,'(6X,A)')'MAIN SURFACE REFERS ',
     .                      'TO HYPER-ELLIPSOIDAL SURFACE'
      ENDIF
C  
C------------
      RETURN
C--------------------------------------------------------------
 1000 FORMAT(/1x,'  INTERFACE NUMBER :',i10,1x,a)
C--------------------------------------------------------------
 
 1717 FORMAT(//
     .    '    TYPE==17 SURFACE/SURFACE 16 NODES THICK SHELL' /,
     .    '              HERTZ PRESSURE FORMULATION' //,
     .    '    FIRST SOLID ELEMENT GROUP. . . . . . .  . . ',i10/,
     .    '    SECOND SOLID ELEMENT GROUP . . . . . .  . . ',i10/,
     .    '    FRICTION FACTOR . . . . . . . . . . . . . . ',1pg20.13/,
     .    '    START TIME. . . . . . . . . . . . . . . . . ',1pg20.13/,
     .    '    STOP TIME . . . . . . . . . . . . . . . . . ',1pg20.13/)
      END