i22cor3.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/.
!||====================================================================
!||    i22cor3          ../engine/source/interfaces/int22/i22cor3.F
!||--- called by ------------------------------------------------------
!||    i22mainf         ../engine/source/interfaces/int22/i22mainf.F
!||--- uses       -----------------------------------------------------
!||    i22bufbric_mod   ../common_source/modules/interfaces/cut-cell-search_mod.F
!||    i22edge_mod      ../common_source/modules/interfaces/cut-cell-buffer_mod.F
!||    i22tri_mod       ../common_source/modules/interfaces/cut-cell-search_mod.F
!||    tri7box          ../engine/share/modules/tri7box.F
!||====================================================================
      SUBROUTINE i22cor3(
     1                JLT     ,X         ,IRECT  ,NSV     ,CAND_E ,
     2                CAND_B  ,STF       ,STFN   ,X1      ,X2     ,
     3                X3      ,X4        ,Y1     ,Y2      ,Y3     ,
     4                Y4      ,Z1        ,Z2     ,Z3      ,Z4     ,
     5                XI      ,YI        ,ZI     ,STIF    ,IX1    ,
     6                IX2     ,IX3       ,IX4    ,NSVG    ,IGAP   ,
     7                GAP     ,GAP_S     ,GAP_M  ,GAPV    ,MS     ,
     8                VXI    ,VYI  ,
     A                VZI     ,MSI       ,NSN    ,V       ,KINET  ,
     B                KINI    ,ITY       ,NIN    ,IGSTI   ,KMIN   ,
     C                KMAX    ,GAPMAX    ,GAPMIN ,IADM    ,INDEX  , 
     D                INTTH   ,TEMP      ,CAND__E,CAND__B ,
     E                TEMPI   ,PHI       ,AREAS  ,IELEC   ,AREASI ,
     F                IELECI  ,NODNX_SMS ,NSMS   ,GAP_S_L ,GAP_M_L,
     G                BUFBRIC ,JTASK)
C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C Interface Type22 (/INTER/TYPE22) is an FSI coupling method based on cut cell method. 
C   This experimental cut cell method is not completed, abandoned, and is not an official option.
C
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE tri7box
      USE i22bufbric_mod
      USE i22tri_mod
      USE i22edge_mod
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   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "comlock.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IRECT(4,*), NSV(*), CAND_E(*), CAND_B(*),KINET(*),KINI(*),
     .        JLT,IDT, NOINT,IGAP ,NSN, ITY, NIN, IGSTI,
     .        IADM,INTTH, CAND__E(*), CAND__B(*)
      INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
     .        NSVG(MVSIZ),IELEC(*),IELECI(MVSIZ),  NSMS(MVSIZ), 
     .        NODNX_SMS(*), BUFBRIC(*), JTASK,INDEX(*)
      my_real
     .   GAP, X(3,*), STF(*), STFN(*),GAP_S(*),GAP_M(*),
     .   MS(*), V(3,*), TEMP(*),AREAS(*),PHI(*),
     .   TEMPI(*),AREASI(*),GAP_S_L(*),GAP_M_L(*)
      my_real
     .   X1(MVSIZ), X2(MVSIZ), X3(MVSIZ), X4(MVSIZ),
     .   Y1(MVSIZ), Y2(MVSIZ), Y3(MVSIZ), Y4(MVSIZ),
     .   Z1(MVSIZ), Z2(MVSIZ), Z3(MVSIZ), Z4(MVSIZ),
     .   XI(MVSIZ), YI(MVSIZ), ZI(MVSIZ), STIF(MVSIZ),
     .   gapv(mvsiz),
     .   vxi(mvsiz), vyi(mvsiz), vzi(mvsiz), msi(mvsiz),
     .   kmin, kmax, gapmax, gapmin,scut
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I ,J  ,IL, EI, NN, IG,JFT, IX,  BI, NI
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------  
C -------------------------------------------------------------
C     DEFINITION OF GLOBAL ARRAY FOR SHELLS (LOCAL + REMOTE ONES)
C -------------------------------------------------------------
C      IRECT_L( 1: 4, :)                         : Four Node ITAB( ID) : global id to manage unicity with XREM
C      IRECT_L( 5: 8, :)                         : X-coordinates
C      IRECT_L( 9:12, :)                         : Y-coordinates
C      IRECT_L(13:16, :)                         : Z-coordinates
C      IRECT_L(17:19, :)                         : Minimum X,Y,Z
C      IRECT_L(20:22, :)                         : Maximum X,Y,Z
C      IRECT_L(23   , :)                         : Stiffness
C      IRECT_L(24:26, :)                         : Lagrangian Velocity (mean)
C      IRECT_L(    :,1+4*ESHIFT:4*ESHIFT+NRTM_T) : element data on current Thread
C      IRECT_L(    :,NRTM+1:NRTM + NSNR)         : remotes elements data
C-----------------------------------------------
C
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
 
          !candidates retaines from each nodes
          !INDEX is index array for current thread
          DO i=1,jlt
            cand_e(i) = cand__e(index(i))
            cand_b(i) = cand__b(index(i))      
          ENDDO
          
 
          DO i=1,jlt
           bi = cand_b(i)    !up to i22intersect() : cand_b is related to couple I. It corresponds to brick IID = bufbric(cand_b(i))
                               !after that cand_b(i) is replaced by address in brick_list
           ei = iabs(cand_e(i))   
           ni=bi
           
           !-------------------------------!
           !RETRIEVING FACES COORDINATES   !
           !-------------------------------! 
 
           ix=irect(1,ei)
           ix1(i)=ix
           x1(i)=x(1,ix)
           y1(i)=x(2,ix)
           z1(i)=x(3,ix)
 
           ix=irect(2,ei)
           ix2(i)=ix
           x2(i)=x(1,ix)
           y2(i)=x(2,ix)
           z2(i)=x(3,ix)
 
           ix=irect(3,ei)
           ix3(i)=ix
           x3(i)=x(1,ix)
           y3(i)=x(2,ix)
           z3(i)=x(3,ix)
 
           ix=irect(4,ei)
           ix4(i)=ix
           x4(i)=x(1,ix)
           y4(i)=x(2,ix)
           z4(i)=x(3,ix)
           
           stif(i) = zero
 
          END DO
 
      RETURN
      END
C