hm_read_inter_type14.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_type14   ../starter/source/interfaces/int14/hm_read_inter_type14.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       -----------------------------------------------------
!||    message_mod            ../starter/share/message_module/message_mod.F
!||    submodel_mod           ../starter/share/modules1/submodel_mod.F
!||====================================================================
      SUBROUTINE hm_read_inter_type14(
     1        IPARI       ,STFAC     ,FRIGAP     ,NOINT     ,
     2        IGRNOD     ,IGRSURF    ,NPC        ,TITR      ,LSUBMODEL,
     3        UNITAB    )
C============================================================================
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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NOINT
      INTEGER IPARI(*),NPC(*)
      my_real STFAC
      my_real frigap(*)
      CHARACTER(LEN=NCHARTITLE),INTENT(IN) :: TITR
      TYPE(submodel_data),INTENT(IN)::LSUBMODEL(*)
      TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
C-----------------------------------------------
      TYPE (GROUP_)  ,TARGET, DIMENSION(NGRNOD)  :: IGRNOD
      TYPE (SURF_)   ,TARGET , DIMENSION(NSURF)   :: IGRSURF
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "units_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      CHARACTER(LEN=NCHARTITLE) :: TITR1
      INTEGER ISU1,ISU2,I,J,L, NTYP,IS1, IS2,NLO,NFRIC,NDAMP1,NDAMP2,NCURS,ISU20,INTKG
      my_real FRIC,GAP,STARTT,STOPT,VISC
      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=======================================================================
C     READING PENALTY INTERFACE /INTER/TYPE14
C=======================================================================
 
C Initializations
      is1=0
      is2=0
      nlo = 0
      nfric = 0
      ndamp1 = 0
      ndamp2 = 0
      intkg=0
C
      fric = zero
      gap = zero
      startt = zero
      stopt=ep30
      visc = zero
C
      ntyp = 14
      ipari(15)=noint
      ipari(7)=ntyp
 
      is_available = .false.
C--------------------------------------------------
C EXTRACT DATAS (INTEGER VALUES)
C--------------------------------------------------
      CALL hm_get_intv('secondaryentityids',isu1,is_available,lsubmodel)
      CALL hm_get_intv('mainentityids',isu2,is_available,lsubmodel)
      CALL hm_get_intv('Iload',nlo,is_available,lsubmodel)
      CALL hm_get_intv('IFRIC',nfric,is_available,lsubmodel)
      CALL hm_get_intv('FUN_A1',ndamp1,is_available,lsubmodel)
      CALL hm_get_intv('FUN_A2',ndamp2,is_available,lsubmodel)
C--------------------------------------------------
C EXTRACT DATAS (REAL VALUES)
C--------------------------------------------------
      CALL hm_get_floatv('STIFF1',stfac,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('FRIC',fric,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('VISC',visc,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('GAP',gap,is_available,lsubmodel,unitab)
C
C....* CHECKS *.............
 
      is1=2
      is2=4
      ingr2usr => igrnod(1:ngrnod)%ID
      IF(isu1/=0)isu1=ngr2usr(isu1,ingr2usr,ngrnod)
         isu20=isu2
         ingr2usr => igrsurf(1:nsurf)%ID
         isu2=ngr2usr(isu2,ingr2usr,nsurf)
         IF ( igrsurf(isu2)%TYPE/=100
     .    .AND.igrsurf(isu2)%TYPE/=101) THEN
          titr1 = igrsurf(isu20)%TITLE
             CALL ancmsg(msgid=111,
     .                   msgtype=msgerror,
     .                   anmode=aninfo,
     .                   i1=noint,
     .                   c1=titr,
     .                   i2=isu20,
     .                   c2=titr1)
      END IF
C-----
 
        IF (nlo==0) GOTO 11
        DO ncurs=1,nfunct
          IF (nlo==npc(nfunct+1+ncurs)) THEN
            ipari(8)=ncurs
            GOTO 11
          ENDIF
        ENDDO
        CALL ancmsg(msgid=113,
     .              msgtype=msgerror,
     .              anmode=aninfo,
     .              i1=noint,
     .              c1=titr,
     .              i2=nlo)
 11     CONTINUE
        IF (nfric==0) GOTO 12
        DO ncurs=1,nfunct
          IF (nfric==npc(nfunct+1+ncurs)) THEN
            ipari(9)=ncurs
            GOTO 12
          ENDIF
        ENDDO
        CALL ancmsg(msgid=113,
     .              msgtype=msgerror,
     .              anmode=aninfo,
     .              i1=noint,
     .              c1=titr,
     .              i2=nfric)
 12     CONTINUE
        IF (ndamp1==0) GOTO 13
        DO ncurs=1,nfunct
          IF (ndamp1==npc(nfunct+1+ncurs)) THEN
            ipari(10)=ncurs
            GOTO 13
          ENDIF
        ENDDO
        CALL ancmsg(msgid=113,
     .              msgtype=msgerror,
     .              anmode=aninfo,
     .              i1=noint,
     .              c1=titr,
     .              i2=ndamp1)
 13     CONTINUE
        IF (ndamp2==0) GOTO 14
        DO ncurs=1,nfunct
          IF (ndamp2==npc(nfunct+1+ncurs)) THEN
            ipari(11)=ncurs
            GOTO 14
          ENDIF
        ENDDO
        CALL ancmsg(msgid=113,
     .              msgtype=msgerror,
     .              anmode=aninfo,
     .              i1=noint,
     .              c1=titr,
     .              i2=ndamp2)
 14     CONTINUE
 
C.......* Storage IPARI FRIGAP *........
      ipari(45)=isu1
      ipari(46)=isu2
      ipari(13)=is1*10+is2
C
      startt=zero
      stopt =ep30
C
C.....* Storage IPARI FRIGAP *.......
      frigap(1)=fric
      frigap(2)=gap
      frigap(3)=startt
      frigap(11)=stopt
      frigap(14)=visc
 
C------------------------------------------------------------
C  General Storage IPARI FRIGAP 
C------------------------------------------------------------
      ipari(65) = intkg
C------------------------------------------------------------
C     PRINTOUT
C------------------------------------------------------------
C
      WRITE(iout,1514)
     .     stfac,nlo,fric,nfric,visc,ndamp1,ndamp2,gap,
     .     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--------------------------------------------------------------
 1000 FORMAT(/1x,'  INTERFACE NUMBER :',i10,1x,a)
C------------
      RETURN
 
 
 1514 FORMAT(//
     .    '    TYPE==14   NODE to SURFACE                ' //,
     .    '    INTERFACE STIFFNESS . . . . . . . . . . . . ',1pg20.13/,
     .    '    FUNCTION FOR ELASTIC CONTACT  . . . . . . . ',i10/,
     .    '    FRICTION COEFFICIENT  . . . . . . . . . . . ',1pg20.13/,
     .    '    FUNCTION FOR FRICTION . . . . . . . . . . . ',i10/,
     .    '    NORMAL DAMPING FACTOR . . . . . . . . . . . ',1pg20.13/,
     .    '    FUNCTION FOR DAMPING VERSUS VELOCITY  . . . ',i10/,
     .    '    FUNCTION FOR DAMPING VERSUS ELASTIC FORCE . ',i10/,
     .    '    MINIMUM GAP . . . . . . . . . . . . . . . . ',1pg20.13/,
     .    '    START TIME. . . . . . . . . . . . . . . . . ',1pg20.13/,
     .    '    STOP TIME . . . . . . . . . . . . . . . . . ',1pg20.13/)
 
      END