hm_read_inter_lagmul_type16.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_lagmul_type16   ../starter/source/interfaces/int16/hm_read_inter_lagmul_type16.F
!||--- called by ------------------------------------------------------
!||    hm_read_inter_lagmul          ../starter/source/interfaces/reader/hm_read_inter_lagmul.F
!||--- calls      -----------------------------------------------------
!||    ancmsg                        ../starter/source/output/message/message.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_lagmul_type16(
     .           IPARI      ,STFAC      ,FRIGAP     ,NOINT     ,
     .           IGRNOD     ,IGRBRIC    ,LSUBMODEL )
C============================================================================
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE groupdef_mod
      USE submodel_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ISU1,ISU2,NOINT
      INTEGER IPARI(*)
      my_real frigap(*) , stfac
      TYPE (GROUP_)  ,TARGET ,DIMENSION(NGRNOD)  :: IGRNOD
      TYPE (GROUP_)  ,TARGET, DIMENSION(NGRBRIC) :: IGRBRIC
      TYPE(submodel_data) LSUBMODEL(*)
C----------------------s-------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "scr06_c.inc"
#include      "units_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: ISU10,ISU20,IDUM,MULTIMP,ITIED,NTYP,IS1,IS2
      my_real :: STARTT,BUMULT,STOPT,GAP
      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     reading /INTER/LAGMUL/TYPE16
C=======================================================================
 
C Initializations
      ntyp = 16
      startt = zero
      stopt  = ep20
      stfac  = one_fifth
      gap    = zero
c------------------------------------------------------------
c  Card1
      CALL hm_get_intv('secondaryentityids ',isu10,is_available,lsubmodel)
      CALL hm_get_intv('mainentityids',isu20,is_available,lsubmodel)
c
c  Card2
      CALL hm_get_intv('Itied' ,itied,is_available,lsubmodel)
c------------------------------------------------------------
c     CHECKS
c------------------------------------------------------------
      is1 = 0
      is2 = 5
      ingr2usr => igrnod(1:ngrnod)%ID
      isu1 = ngr2usr(isu10,ingr2usr,ngrnod)
      ingr2usr => igrbric(1:ngrbric)%ID
      isu2 = ngr2usr(isu20,ingr2usr,ngrbric)
c
      bumult = bmul0
      multimp=12
c
      IF (nspmd > 1) THEN
          CALL ancmsg(msgid=755,msgtype=msgerror,anmode=aninfo,c1='TYPE 16 LAGRANGE INTERFACE')
      END IF
c------------------------------------------------------------
c     STORAGE
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(2) = gap
      frigap(3) = startt
      frigap(4) = bumult
      frigap(5) = one
      frigap(11)= stopt
C------------------------------------------------------------
C     PRINTOUT
C------------------------------------------------------------
      WRITE(iout,1000)isu10,isu20,itied,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
      WRITE(iout,'(A)')' '
c-----------
      RETURN
c--------------------------------------------------------------
 1000 FORMAT(//
     .    '    TYPE 16   NODES TO 16 NODES THICK SHELL  ' //,
     .    '    NODE GROUP . . . . . . .  . . . . . . . . . ',i10/,
     .    '    SOLID ELEMENT GROUP. . . . . . .  . . . . . ',i10/,
     .    '    ITIED . . . . . . . . . . . . . . . . . . . ',i10/,
     .    '       0: SLIDING '/,
     .    '       1: TIED (during contact)'/,
     .    '       2: TIED (no rebound)'/,
     .    '    START TIME. . . . . . . . . . . . . . . . . ',1pg20.13/,
     .    '    STOP TIME . . . . . . . . . . . . . . . . . ',1pg20.13/)
c--------------------------------------------------------------
      END