i10corp3.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/.
!||====================================================================
!||    i10corp3    ../engine/source/interfaces/int10/i10corp3.F
!||--- called by ------------------------------------------------------
!||    i10forcf3   ../engine/source/interfaces/int10/i10ke3.F
!||--- uses       -----------------------------------------------------
!||    imp_intm    ../engine/share/modules/imp_intm.F
!||    tri7box     ../engine/share/modules/tri7box.F
!||====================================================================
      SUBROUTINE i10corp3(JLT   ,X     ,IRECT ,CAND_E  ,CAND_N  ,
     1                   STIF  ,H1    ,H2    ,H3      ,H4       ,
     2                   N1    ,N2    ,N3    ,NT1     ,NT2      ,
     3                   NT3   ,IX1   ,IX2   ,IX3     ,IX4      ,
     4                   NSVG  ,VXI   ,VYI   ,VZI     ,MSI      ,
     5                   DXI   ,DYI   ,DZI   ,NSN     ,NIN      ,
     6                   JLT_NEW,LREM )
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE tri7box
      USE imp_intm
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IRECT(4,*), CAND_E(*), CAND_N(*),
     .         JLT,JLT_NEW,NIN,NSN,LREM
      INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
     .        NSVG(MVSIZ)
C     REAL
C     REAL
      my_real
     .   vxi(mvsiz), vyi(mvsiz), vzi(mvsiz), msi(mvsiz),
     .   x(3,*), stif(*),n1(mvsiz), n2(mvsiz), n3(mvsiz),
     .   nt1(mvsiz), nt2(mvsiz), nt3(mvsiz),
     .   h1(mvsiz), h2(mvsiz), h3(mvsiz), h4(mvsiz),
     .   dxi(mvsiz),dyi(mvsiz),dzi(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I ,J  ,NI, L, NN, NS ,NE
      my_real
     .   T1,T2,T3,NORM
C-----------------------------------------------
C
C
        DO i=1,jlt
         ni = cand_n(i)
         l  = cand_e(i)
         IF(ni>nsn)THEN
           nn = ni - nsn
           jlt_new = jlt_new + 1
           vxi(jlt_new)= vfi(nin)%P(1,nn)
           vyi(jlt_new)= vfi(nin)%P(2,nn)
           vzi(jlt_new)= vfi(nin)%P(3,nn)
           msi(jlt_new)= msfi(nin)%P(nn)
           ix1(jlt_new)  = irect(1,l)
           ix2(jlt_new)  = irect(2,l)
           ix3(jlt_new)  = irect(3,l)
           ix4(jlt_new)  = irect(4,l)
             ns=ind_int(nin)%P(nn)
             ne=shf_int(nin)+jlt_new+lrem
           nsvg(jlt_new) = ns
           h1(jlt_new)  = h_e(1,ne)
           h2(jlt_new)  = h_e(2,ne)
           h3(jlt_new)  = h_e(3,ne)
           h4(jlt_new)  = h_e(4,ne)
             stif(jlt_new)=stifs(ne)
             n1(jlt_new)=n_e(1,ne)
             n2(jlt_new)=n_e(2,ne)
             n3(jlt_new)=n_e(3,ne)
C----------displacement
             dxi(jlt_new)=dfi(1,ns)
             dyi(jlt_new)=dfi(2,ns)
             dzi(jlt_new)=dfi(3,ns)
           t1 = x(1,ix3(jlt_new)) - x(1,ix1(jlt_new))
           t2 = x(2,ix3(jlt_new)) - x(2,ix1(jlt_new))
           t3 = x(3,ix3(jlt_new)) - x(3,ix1(jlt_new))
           norm = one/sqrt(t1*t1+t2*t2+t3*t3)
           nt1(jlt_new) = t1*norm
           nt2(jlt_new) = t2*norm
           nt3(jlt_new) = t3*norm
C
         END IF
C
        END DO
C
      RETURN
      END