OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i22assembly.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "comlock.inc"
#include "parit_c.inc"
#include "com01_c.inc"
#include "inter22.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine i22ass0 (jlt, ix1, ix2, ix3, ix4, nsvg, h1, h2, h3, h4, stif, fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, fxi, fyi, fzi, a, stifn, nin, intth, phi, fthe, phi1, phi2, phi3, phi4)
subroutine i22ass2 (jlt, ix1, ix2, ix3, ix4, itab, nsvg, h1, h2, h3, h4, stif, fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, fxi, fyi, fzi, fskyi, isky, niskyfi, nin, noint, intth, phi, ftheskyi, phi1, phi2, phi3, phi4, cb_loc, ce_loc, irect, ixs)

Function/Subroutine Documentation

◆ i22ass0()

subroutine i22ass0 ( integer jlt,
integer, dimension(mvsiz) ix1,
integer, dimension(mvsiz) ix2,
integer, dimension(mvsiz) ix3,
integer, dimension(mvsiz) ix4,
integer, dimension(mvsiz) nsvg,
h1,
h2,
h3,
h4,
stif,
fx1,
fy1,
fz1,
fx2,
fy2,
fz2,
fx3,
fy3,
fz3,
fx4,
fy4,
fz4,
fxi,
fyi,
fzi,
a,
stifn,
integer nin,
integer intth,
phi,
fthe,
phi1,
phi2,
phi3,
phi4 )

Definition at line 30 of file i22assembly.F.

37C-----------------------------------------------
38C D e s c r i p t i o n
39C-----------------------------------------------
40C Interface Type22 (/INTER/TYPE22) is an FSI coupling method based on cut cell method.
41C This experimental cut cell method is not completed, abandoned, and is not an official option.
42C-----------------------------------------------
43C M o d u l e s
44C-----------------------------------------------
45 USE tri7box
46C-----------------------------------------------
47C I m p l i c i t T y p e s
48C-----------------------------------------------
49#include "implicit_f.inc"
50C-----------------------------------------------
51C G l o b a l P a r a m e t e r s
52C-----------------------------------------------
53#include "mvsiz_p.inc"
54C-----------------------------------------------
55C D u m m y A r g u m e n t s
56C-----------------------------------------------
57 INTEGER JLT, NIN,INTTH,
58 . IX1(MVSIZ),IX2(MVSIZ),IX3(MVSIZ),IX4(MVSIZ),NSVG(MVSIZ)
59 my_real
60 . h1(mvsiz),h2(mvsiz),h3(mvsiz),h4(mvsiz),stif(mvsiz),
61 . fx1(mvsiz),fy1(mvsiz),fz1(mvsiz),
62 . fx2(mvsiz),fy2(mvsiz),fz2(mvsiz),
63 . fx3(mvsiz),fy3(mvsiz),fz3(mvsiz),
64 . fx4(mvsiz),fy4(mvsiz),fz4(mvsiz),
65 . fxi(mvsiz),fyi(mvsiz),fzi(mvsiz),
66 . a(3,*), stifn(*),phi(*), fthe(*),
67 . phi1(*), phi2(*), phi3(*), phi4(*)
68C-----------------------------------------------
69C L o c a l V a r i a b l e s
70C-----------------------------------------------
71 INTEGER I, J1, IG
72C
73C-----------------------------------------------
74C S o u r c e L i n e s
75C-----------------------------------------------
76 DO i=1,jlt
77 j1=ix1(i)
78 a(1,j1)=a(1,j1)+fx1(i)
79 a(2,j1)=a(2,j1)+fy1(i)
80 a(3,j1)=a(3,j1)+fz1(i)
81! STIFN(J1) = STIFN(J1) + STIF(I)*ABS(H1(I))
82C
83 j1=ix2(i)
84 a(1,j1)=a(1,j1)+fx2(i)
85 a(2,j1)=a(2,j1)+fy2(i)
86 a(3,j1)=a(3,j1)+fz2(i)
87! STIFN(J1) = STIFN(J1) + STIF(I)*ABS(H2(I))
88C
89 j1=ix3(i)
90 a(1,j1)=a(1,j1)+fx3(i)
91 a(2,j1)=a(2,j1)+fy3(i)
92 a(3,j1)=a(3,j1)+fz3(i)
93! STIFN(J1) = STIFN(J1) + STIF(I)*ABS(H3(I))
94C
95 j1=ix4(i)
96 a(1,j1)=a(1,j1)+fx4(i)
97 a(2,j1)=a(2,j1)+fy4(i)
98 a(3,j1)=a(3,j1)+fz4(i)
99! STIFN(J1) = STIFN(J1) + STIF(I)*ABS(H4(I))
100 ENDDO
101
102C
103 RETURN
my_real
#define my_real
Definition cppsort.cpp:32

◆ i22ass2()

subroutine i22ass2 ( integer jlt,
integer, dimension(mvsiz) ix1,
integer, dimension(mvsiz) ix2,
integer, dimension(mvsiz) ix3,
integer, dimension(mvsiz) ix4,
integer, dimension(*) itab,
integer, dimension(mvsiz) nsvg,
h1,
h2,
h3,
h4,
stif,
fx1,
fy1,
fz1,
fx2,
fy2,
fz2,
fx3,
fy3,
fz3,
fx4,
fy4,
fz4,
fxi,
fyi,
fzi,
fskyi,
integer, dimension(*) isky,
integer niskyfi,
integer nin,
integer noint,
integer intth,
phi,
ftheskyi,
phi1,
phi2,
phi3,
phi4,
integer, dimension(mvsiz) cb_loc,
integer, dimension(mvsiz) ce_loc,
integer, dimension(4,*) irect,
integer, dimension(nixs,*) ixs )

Definition at line 120 of file i22assembly.F.

128C-----------------------------------------------
129C M o d u l e s
130C-----------------------------------------------
131 USE tri7box
132 USE i22tri_mod
133 USE i22bufbric_mod
134 USE message_mod
135C-----------------------------------------------
136C I m p l i c i t T y p e s
137C-----------------------------------------------
138#include "implicit_f.inc"
139#include "comlock.inc"
140C-----------------------------------------------
141C G l o b a l P a r a m e t e r s
142C-----------------------------------------------
143#include "mvsiz_p.inc"
144C-----------------------------------------------
145C C o m m o n B l o c k s
146C-----------------------------------------------
147#include "parit_c.inc"
148#include "com01_c.inc"
149#include "inter22.inc"
150C-----------------------------------------------
151C D u m m y A r g u m e n t s
152C-----------------------------------------------
153 INTEGER JLT,NISKYFI,NIN,NOINT,INTTH,
154 . ISKY(*),ITAB(*),CB_LOC(MVSIZ),CE_LOC(MVSIZ),
155 . IX1(MVSIZ),IX2(MVSIZ),IX3(MVSIZ),IX4(MVSIZ),NSVG(MVSIZ),
156 . IRECT(4,*),IXS(NIXS,*)
157 my_real
158 . h1(mvsiz),h2(mvsiz),h3(mvsiz),h4(mvsiz),stif(mvsiz),
159 . fx1(mvsiz),fy1(mvsiz),fz1(mvsiz),
160 . fx2(mvsiz),fy2(mvsiz),fz2(mvsiz),
161 . fx3(mvsiz),fy3(mvsiz),fz3(mvsiz),
162 . fx4(mvsiz),fy4(mvsiz),fz4(mvsiz),
163 . fxi(mvsiz),fyi(mvsiz),fzi(mvsiz),
164 . fskyi(lskyi,nfskyi),ftheskyi(lskyi),phi(mvsiz),
165 . phi1(*),phi2(*) ,phi3(*) ,phi4(*)
166C-----------------------------------------------
167C L o c a l V a r i a b l e s
168C-----------------------------------------------
169 INTEGER I, J1, IG, NISKYL1, NISKYL,IGP,IGM,IDR,NISKYFIL, K
170C
171 niskyl1 = 0
172 niskyl1 = niskyl1 + 4*jlt
173C-----------------------------------------------
174C P r e - C o n d i t i o n
175C-----------------------------------------------
176 IF(int22==0)RETURN
177C-----------------------------------------------
178C S o u r c e L i n e s
179C-----------------------------------------------
180#include "lockon.inc"
181 niskyl = nisky ;
182 nisky = nisky + niskyl1
183 niskyfil = niskyfi
184 niskyfi = niskyfi
185#include "lockoff.inc"
186C
187 IF (niskyl+niskyl1 > lskyi) THEN
188 !CALL ANCMSG(MSGID=26,ANMODE=ANINFO)
189 print *, "** ERROR INTER22 : force assembly issue with PARITH/ON option"
190 CALL arret(2)
191 ENDIF
192 IF (niskyfil > nlskyfi(nin)) THEN
193 CALL ancmsg(msgid=26,anmode=aninfo)
194 CALL arret(2)
195 ENDIF
196C
197 IF(intth == 0 ) THEN
198 DO i=1,jlt
199 niskyl = niskyl + 1
200 fskyi(niskyl,1) = fx1(i)
201 fskyi(niskyl,2) = fy1(i)
202 fskyi(niskyl,3) = fz1(i)
203 fskyi(niskyl,4) = zero
204 isky(niskyl) = ix1(i)
205 if(ibug22_fcontass==-1)then
206 k = abs(ce_loc(i))
207 print *, " assembly iparit2 - node_id : ", itab(ix1(i))
208 print *, " +--BRIQUE : ", ixs(11,brick_list(nin,cb_loc(i))%ID)
209 print *, " +--FACETTE: ", itab(irect(1:4,k))
210 print *, " +--NCYCLE: ", ncycle
211 write(*,fmt='(A,3E30.16)') " +--Fxyz1: ", fx1(i),fy1(i),fz1(i)
212 endif
213 ENDDO
214 DO i=1,jlt
215 niskyl = niskyl + 1
216 fskyi(niskyl,1) = fx2(i)
217 fskyi(niskyl,2) = fy2(i)
218 fskyi(niskyl,3) = fz2(i)
219 fskyi(niskyl,4) = zero
220 isky(niskyl) = ix2(i)
221 if(ibug22_fcontass==-1)then
222 print *, " assembly iparit2 - node_id : ", itab(ix2(i))
223 write(*,fmt='(A,3E30.16)') " +--Fxyz2: ", fx2(i),fy2(i),fz2(i)
224 endif
225 ENDDO
226 DO i=1,jlt
227 niskyl = niskyl + 1
228 fskyi(niskyl,1) = fx3(i)
229 fskyi(niskyl,2) = fy3(i)
230 fskyi(niskyl,3) = fz3(i)
231 fskyi(niskyl,4) = zero
232 isky(niskyl) = ix3(i)
233 if(ibug22_fcontass==-1)then
234 print *, " assembly iparit2 - node_id : ", itab(ix3(i))
235 write(*,fmt='(A,3E30.16)') " +--Fxyz3: ", fx3(i),fy3(i),fz3(i)
236 endif
237 ENDDO
238 DO i=1,jlt
239 niskyl = niskyl + 1
240 fskyi(niskyl,1) = fx4(i)
241 fskyi(niskyl,2) = fy4(i)
242 fskyi(niskyl,3) = fz4(i)
243 fskyi(niskyl,4) = zero
244 isky(niskyl) = ix4(i)
245 if(ibug22_fcontass==-1)then
246 print *, " assembly iparit2 - node_id : ", itab(ix4(i))
247 write(*,fmt='(A,3E30.16)') " +--Fxyz4: ", fx4(i),fy4(i),fz4(i)
248 endif
249 ENDDO
250
251 !### Thermique ###!
252 ELSE
253 !not planned yet
254 endif!(INTTH == 0 )
255C
256 RETURN
i22bufbric_mod
Definition cut-cell-search_mod.F:95
i22bufbric_mod::brick_list
type(brick_entity), dimension(:,:), allocatable, target brick_list
Definition cut-cell-search_mod.F:105
i22tri_mod
Definition cut-cell-search_mod.F:212
i22tri_mod::ibug22_fcontass
integer ibug22_fcontass
Definition cut-cell-search_mod.F:337
message_mod
Definition message_mod.F:1249
tri7box::nlskyfi
integer, dimension(:), allocatable nlskyfi
Definition tri7box.F:512
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889
arret
subroutine arret(nn)
Definition arret.F:87