hm_read_inter_type10.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_type10   ../starter/source/interfaces/int10/hm_read_inter_type10.F
!||--- called by ------------------------------------------------------
!||    hm_read_inter_struct   ../starter/source/interfaces/reader/hm_read_inter_struct.F
!||--- calls      -----------------------------------------------------
!||    ancmsg                 ../starter/source/output/message/message.F
!||    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       -----------------------------------------------------
!||    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_inter_type10(
     .        IPARI    ,STFAC    ,FRIGAP   ,IGRNOD   ,IGRSURF  ,
     .        LSUBMODEL,UNITAB   ,TITR     ,NOINT    )
C=======================================================================
C   M o d u l e s
C-----------------------------------------------
      USE unitab_mod
      USE message_mod
      USE groupdef_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   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr06_c.inc"
#include      "com01_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
      my_real  :: STFAC
      INTEGER, DIMENSION(*)  :: IPARI
      my_real, DIMENSION(*)  :: frigap
      CHARACTER(LEN=NCHARTITLE) ,INTENT(IN)             :: TITR
      TYPE (UNIT_TYPE_) ,INTENT(IN) ::  UNITAB
      TYPE (SUBMODEL_DATA), DIMENSION(*),INTENT(IN) :: LSUBMODEL
C-----------------------------------------------
      TYPE (GROUP_)  ,TARGET ,DIMENSION(NGRNOD)  :: IGRNOD
      TYPE (SURF_)   ,TARGET ,DIMENSION(NSURF)   :: IGRSURF
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: NTYP,IS1,IS2,ISU1,ISU2,IDELKEEP,ILEV,IGAP,INACTI,IDEL10,MULTIMP,ITIED
      my_real :: FRIC,GAP,STARTT,STOPT,VISC,BUMULT,GAPMAX,FPENMAX
      INTEGER, DIMENSION(:), POINTER :: INGR2USR
      LOGICAL :: IS_AVAILABLE
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER NGR2USR
C=======================================================================
c     Initializations
c--------------------------------------------------------
      ntyp     = 10
      igap     = 0
      ilev     = 0
      multimp  = 4
      idelkeep = 0
      idel10   = 0
      gapmax   = infinity
      fpenmax  = zero
      is_available = .false.
c--------------------------------------------------------
c     Read input fields
c--------------------------------------------------------
card1
      CALL hm_get_intv  ('secondaryentityids'  ,isu1   ,is_available, lsubmodel)      
      CALL hm_get_intv  ('mainentityids' ,isu2   ,is_available, lsubmodel)      
      CALL hm_get_intv  ('NodDel10'        ,idel10 ,is_available, lsubmodel)      
card2
      CALL hm_get_floatv('TYPE10_SCALE'    ,stfac  ,is_available, lsubmodel, unitab)
      CALL hm_get_floatv('GAP'             ,gap    ,is_available, lsubmodel, unitab)
      CALL hm_get_floatv('TSTART'          ,startt ,is_available, lsubmodel, unitab)
      CALL hm_get_floatv('TSTOP'           ,stopt  ,is_available, lsubmodel, unitab)
card2
      CALL hm_get_intv  ('Itied'           ,itied  ,is_available, lsubmodel)      
      CALL hm_get_intv  ('INACTIV'         ,inacti ,is_available, lsubmodel)      
      CALL hm_get_floatv('STIFF_DC'        ,visc   ,is_available, lsubmodel, unitab)
      CALL hm_get_floatv('SORT_FACT'       ,bumult ,is_available, lsubmodel, unitab)
c--------------------------------------------------------
c     Checks
c--------------------------------------------------------
        IF (inacti == 5) CALL ancmsg(msgid=1162,
     .                        msgtype=msgerror,
     .                        anmode=aninfo,
     .                        i1=noint,
     .                        c1=titr)
c
      is1 = 2   ! SECONDARY  surface input by node group
      is2 = 1   ! main surface input by surface Id
      ingr2usr => igrnod(1:ngrnod)%ID
      IF (isu1 /= 0) isu1 = ngr2usr(isu1,ingr2usr,ngrnod)
      ingr2usr => igrsurf(1:nsurf)%ID
      isu2 = ngr2usr(isu2,ingr2usr,nsurf)
      IF (idel10 < 0) THEN
         idelkeep=1
         idel10 = abs(idel10)
      END IF
      IF (idel10 > 2.OR. n2d == 1) idel10 = 0
c
      fric = itied
      IF (stfac == zero) stfac = one_fifth
      IF (visc == zero)  visc  = fiveem2
      IF (bumult == zero) bumult = bmul0
      IF (stopt == zero) stopt = infinity
c--------------------------------------------------------
c     Fill interface buffers
c--------------------------------------------------------
      frigap(1)  = fric
      frigap(2)  = gap
      frigap(3)  = startt
      frigap(4)  = bumult
      frigap(10) = zero     ! only in engine for storing number of couples candidates 
      frigap(11) = stopt
      frigap(14) = visc
      frigap(16) = gapmax
      frigap(27) = fpenmax
      frigap(15) = zero 
c
      ipari(7)  = ntyp
      ipari(11) = 0
      ipari(17) = idel10
      ipari(13) = is1*10 + is2
      ipari(20) = ilev
      ipari(21) = igap
      ipari(22) = inacti
      ipari(65) = 0
      ipari(15) = noint
      ipari(23) = multimp
      ipari(45) = isu1
      ipari(46) = isu2
      ipari(61) = idelkeep
c--------------------------------------------------------
c     Printout
c--------------------------------------------------------
       itied = nint(fric)
       WRITE(iout,1510) itied,stfac,gap,startt,stopt,
     .                  bumult,inacti,visc,multimp
c
       IF (idel10 /= 0) THEN
         WRITE(iout,'(A,A,I5/)')
     .'    DELETION FLAG ON FAILURE OF MAIN ELEMENT',
     .' (1:YES-ALL/2:YES-ANY) SET TO ',idel10
         IF (idelkeep == 1) THEN
           WRITE(iout,'(A)')
     .'    IDEL: DO NOT REMOVE NON-CONNECTED NODES FROM SECONDARY SURFACE'
         ENDIF         
       ENDIF
       WRITE(iout,'(6X,A)')'SECONDARY SURFACE INPUT BY NODES'
       WRITE(iout,'(6X,A)')'MAIN SURFACE INPUT BY SEGMENTS'
 
c--------------------------------------------------------
 1510 FORMAT(//
     .    '    TYPE 10 TIED - AUTO IMPACTING               ' //,
     .    '    ITIED . . . . . . . . . . . . . . . . . . . ',i1/,
     .    '       0: TIED DURING IMPACT - REBOUND AUTORIZED'/,
     .    '       1: TIED AFTER IMPACT NO REBOUND AUTORIZED'/,
     .    '    STIFFNESS FACTOR. . . . . . . . . . . . . ',1pg20.13/,
     .    '    MINIMUM GAP . . . . . . . . . . . . . . . ',1pg20.13/,
     .    '    START TIME. . . . . . . . . . . . . . . . ',1pg20.13/,
     .    '    STOP TIME . . . . . . . . . . . . . . . . ',1pg20.13/,
     .    '    BUCKET FACTOR . . . . . . . . . . . . . . ',1pg20.13/,
     .    '    DE-ACTIVATION OF INITIAL PENETRATIONS . . ',i10/,
     .    '    CRITICAL DAMPING FACTOR . . . . . . . . . ',1pg20.13/,
     .    '    MEAN POSSIBLE NUMBER OF IMPACT/NODE . . . ',i5/)
c-----------
      RETURN
      END