OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i22cor3.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "comlock.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine i22cor3 (jlt, x, irect, nsv, cand_e, cand_b, stf, stfn, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, stif, ix1, ix2, ix3, ix4, nsvg, igap, gap, gap_s, gap_m, gapv, ms, vxi, vyi, vzi, msi, nsn, v, kinet, kini, ity, nin, igsti, kmin, kmax, gapmax, gapmin, iadm, index, intth, temp, cand__e, cand__b, tempi, phi, areas, ielec, areasi, ieleci, nodnx_sms, nsms, gap_s_l, gap_m_l, bufbric, jtask)

Function/Subroutine Documentation

◆ i22cor3()

subroutine i22cor3 ( integer jlt,
x,
integer, dimension(4,*) irect,
integer, dimension(*) nsv,
integer, dimension(*) cand_e,
integer, dimension(*) cand_b,
stf,
stfn,
x1,
x2,
x3,
x4,
y1,
y2,
y3,
y4,
z1,
z2,
z3,
z4,
xi,
yi,
zi,
stif,
integer, dimension(mvsiz) ix1,
integer, dimension(mvsiz) ix2,
integer, dimension(mvsiz) ix3,
integer, dimension(mvsiz) ix4,
integer, dimension(mvsiz) nsvg,
integer igap,
gap,
gap_s,
gap_m,
gapv,
ms,
vxi,
vyi,
vzi,
msi,
integer nsn,
v,
integer, dimension(*) kinet,
integer, dimension(*) kini,
integer ity,
integer nin,
integer igsti,
kmin,
kmax,
gapmax,
gapmin,
integer iadm,
integer, dimension(*) index,
integer intth,
temp,
integer, dimension(*) cand__e,
integer, dimension(*) cand__b,
tempi,
phi,
areas,
integer, dimension(*) ielec,
areasi,
integer, dimension(mvsiz) ieleci,
integer, dimension(*) nodnx_sms,
integer, dimension(mvsiz) nsms,
gap_s_l,
gap_m_l,
integer, dimension(*) bufbric,
integer jtask )

Definition at line 33 of file i22cor3.F.

49C-----------------------------------------------
50C D e s c r i p t i o n
51C-----------------------------------------------
52C Interface Type22 (/INTER/TYPE22) is an FSI coupling method based on cut cell method.
53C This experimental cut cell method is not completed, abandoned, and is not an official option.
54C
55C-----------------------------------------------
56C M o d u l e s
57C-----------------------------------------------
58 USE tri7box
59 USE i22bufbric_mod
60 USE i22tri_mod
61 USE i22edge_mod
62C-----------------------------------------------
63C I m p l i c i t T y p e s
64C-----------------------------------------------
65#include "implicit_f.inc"
66C-----------------------------------------------
67C G l o b a l P a r a m e t e r s
68C-----------------------------------------------
69#include "mvsiz_p.inc"
70C-----------------------------------------------
71C C o m m o n B l o c k s
72C-----------------------------------------------
73#include "comlock.inc"
74C-----------------------------------------------
75C D u m m y A r g u m e n t s
76C-----------------------------------------------
77 INTEGER IRECT(4,*), NSV(*), CAND_E(*), CAND_B(*),KINET(*),KINI(*),
78 . JLT,IDT, NOINT,IGAP ,NSN, ITY, NIN, IGSTI,
79 . IADM,INTTH, CAND__E(*), CAND__B(*)
80 INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
81 . NSVG(MVSIZ),IELEC(*),IELECI(MVSIZ), NSMS(MVSIZ),
82 . NODNX_SMS(*), BUFBRIC(*), JTASK,INDEX(*)
83 my_real
84 . gap, x(3,*), stf(*), stfn(*),gap_s(*),gap_m(*),
85 . ms(*), v(3,*), temp(*),areas(*),phi(*),
86 . tempi(*),areasi(*),gap_s_l(*),gap_m_l(*)
87 my_real
88 . x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz),
89 . y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
90 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
91 . xi(mvsiz), yi(mvsiz), zi(mvsiz), stif(mvsiz),
92 . gapv(mvsiz),
93 . vxi(mvsiz), vyi(mvsiz), vzi(mvsiz), msi(mvsiz),
94 . kmin, kmax, gapmax, gapmin,scut
95C-----------------------------------------------
96C L o c a l V a r i a b l e s
97C-----------------------------------------------
98 INTEGER I ,J ,IL, EI, NN, IG,JFT, IX, BI, NI
99C-----------------------------------------------
100C S o u r c e L i n e s
101C-----------------------------------------------
102C -------------------------------------------------------------
103C DEFINITION OF GLOBAL ARRAY FOR SHELLS (LOCAL + REMOTE ONES)
104C -------------------------------------------------------------
105C IRECT_L( 1: 4, :) : Four Node ITAB( ID) : global id to manage unicity with XREM
106C IRECT_L( 5: 8, :) : X-coordinates
107C IRECT_L( 9:12, :) : Y-coordinates
108C IRECT_L(13:16, :) : Z-coordinates
109C IRECT_L(17:19, :) : Minimum X,Y,Z
110C IRECT_L(20:22, :) : Maximum X,Y,Z
111C IRECT_L(23 , :) : Stiffness
112C IRECT_L(24:26, :) : Lagrangian Velocity (mean)
113C IRECT_L( :,1+4*ESHIFT:4*ESHIFT+NRTM_T) : element data on current Thread
114C IRECT_L( :,NRTM+1:NRTM + NSNR) : remotes elements data
115C-----------------------------------------------
116C
117C-----------------------------------------------
118C S o u r c e L i n e s
119C-----------------------------------------------
120
121 !candidates retaines from each nodes
122 !INDEX is index array for current thread
123 DO i=1,jlt
124 cand_e(i) = cand__e(index(i))
125 cand_b(i) = cand__b(index(i))
126 ENDDO
127
128
129 DO i=1,jlt
130 bi = cand_b(i) !up to i22intersect() : cand_b is related to couple I. It corresponds to brick IID = bufbric(cand_b(i))
131 !after that cand_b(i) is replaced by address in brick_list
132 ei = iabs(cand_e(i))
133 ni=bi
134
135 !-------------------------------!
136 !RETRIEVING FACES COORDINATES !
137 !-------------------------------!
138
139 ix=irect(1,ei)
140 ix1(i)=ix
141 x1(i)=x(1,ix)
142 y1(i)=x(2,ix)
143 z1(i)=x(3,ix)
144
145 ix=irect(2,ei)
146 ix2(i)=ix
147 x2(i)=x(1,ix)
148 y2(i)=x(2,ix)
149 z2(i)=x(3,ix)
150
151 ix=irect(3,ei)
152 ix3(i)=ix
153 x3(i)=x(1,ix)
154 y3(i)=x(2,ix)
155 z3(i)=x(3,ix)
156
157 ix=irect(4,ei)
158 ix4(i)=ix
159 x4(i)=x(1,ix)
160 y4(i)=x(2,ix)
161 z4(i)=x(3,ix)
162
163 stif(i) = zero
164
165 END DO
166
167 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
i22bufbric_mod
Definition cut-cell-search_mod.F:95
i22edge_mod
Definition cut-cell-buffer_mod.F:48
i22tri_mod
Definition cut-cell-search_mod.F:212