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

Go to the source code of this file.

Functions/Subroutines

subroutine i7cork3 (jlt, x, irect, nsv, cand_e, cand_n, 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, v, igsti, kmin, kmax, gapmax, gapmin, nin, ity, nsn)
subroutine i7corp3 (jlt, x, irect, cand_e, cand_n, stif, h1, h2, h3, h4, n1, n2, n3, ix1, ix2, ix3, ix4, nsvg, vxi, vyi, vzi, msi, dxi, dyi, dzi, nsn, nin, jlt_new, lrem)
subroutine i7corkp3 (jlt, xi, yi, zi, d, dxi, dyi, dzi, nsvg, nin, iupd)

Function/Subroutine Documentation

◆ i7cork3()

subroutine i7cork3 ( integer jlt,
x,
integer, dimension(4,*) irect,
integer, dimension(*) nsv,
integer, dimension(*) cand_e,
integer, dimension(*) cand_n,
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,
v,
integer igsti,
kmin,
kmax,
gapmax,
gapmin,
integer nin,
integer ity,
integer nsn )

Definition at line 35 of file i7cork3.F.

45C-----------------------------------------------
46C M o d u l e s
47C-----------------------------------------------
48 USE tri7box
49C-----------------------------------------------
50C I m p l i c i t T y p e s
51C-----------------------------------------------
52#include "implicit_f.inc"
53C-----------------------------------------------
54C G l o b a l P a r a m e t e r s
55C-----------------------------------------------
56#include "mvsiz_p.inc"
57C-----------------------------------------------
58C C o m m o n B l o c k s
59C-----------------------------------------------
60C-----------------------------------------------
61C D u m m y A r g u m e n t s
62C-----------------------------------------------
63 INTEGER IRECT(4,*), NSV(*), CAND_E(*), CAND_N(*),
64 . JLT,IDT, NOINT,IGAP , IGSTI,NIN,ITY,NSN
65 INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
66 . NSVG(MVSIZ)
67 my_real
68 . gap, x(3,*), stf(*), stfn(*),gap_s(*),gap_m(*),
69 . ms(*), v(3,*)
70 my_real
71 . x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz),
72 . y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
73 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
74 . xi(mvsiz), yi(mvsiz), zi(mvsiz), stif(mvsiz),
75 . gapv(mvsiz),
76 . vxi(mvsiz), vyi(mvsiz), vzi(mvsiz), msi(mvsiz),
77 . kmin, kmax, gapmax,gapmin
78C-----------------------------------------------
79C L o c a l V a r i a b l e s
80C-----------------------------------------------
81 INTEGER I ,J ,IL, L, NN, IG,JFT,NI,IX
82C-----------------------------------------------
83C
84 IF(igap==0)THEN
85 DO i=1,jlt
86 gapv(i)=gap
87 ENDDO
88 ELSE
89 DO i=1,jlt
90 IF(cand_n(i)<=nsn) THEN
91 gapv(i)=gap_s(cand_n(i))+gap_m(cand_e(i))
92 ELSE
93 gapv(i)=gapfi(nin)%P(cand_n(i)-nsn)+gap_m(cand_e(i))
94 ENDIF
95 gapv(i)=min(gapv(i),gapmax)
96 gapv(i)=max(gapmin,gapv(i))
97 ENDDO
98 ENDIF
99C
100 IF(ity==7) THEN
101 DO i=1,jlt
102 ni = cand_n(i)
103 l = cand_e(i)
104 IF(ni<=nsn)THEN
105 ig = nsv(ni)
106 nsvg(i) = ig
107C KINI(I) = KINET(IG)
108 xi(i) = x(1,ig)
109 yi(i) = x(2,ig)
110 zi(i) = x(3,ig)
111 vxi(i) = v(1,ig)
112 vyi(i) = v(2,ig)
113 vzi(i) = v(3,ig)
114 msi(i)= ms(ig)
115 ELSE
116 nn = ni - nsn
117 nsvg(i) = -nn
118C qz KINI(I) = KINFI(NIN)%P(NN)
119 xi(i) = xfi(nin)%P(1,nn)
120 yi(i) = xfi(nin)%P(2,nn)
121 zi(i) = xfi(nin)%P(3,nn)
122 vxi(i)= vfi(nin)%P(1,nn)
123 vyi(i)= vfi(nin)%P(2,nn)
124 vzi(i)= vfi(nin)%P(3,nn)
125 msi(i)= msfi(nin)%P(nn)
126C
127 END IF
128C
129 ix=irect(1,l)
130 ix1(i)=ix
131 x1(i)=x(1,ix)
132 y1(i)=x(2,ix)
133 z1(i)=x(3,ix)
134C
135 ix=irect(2,l)
136 ix2(i)=ix
137 x2(i)=x(1,ix)
138 y2(i)=x(2,ix)
139 z2(i)=x(3,ix)
140C
141 ix=irect(3,l)
142 ix3(i)=ix
143 x3(i)=x(1,ix)
144 y3(i)=x(2,ix)
145 z3(i)=x(3,ix)
146C
147 ix=irect(4,l)
148 ix4(i)=ix
149 x4(i)=x(1,ix)
150 y4(i)=x(2,ix)
151 z4(i)=x(3,ix)
152C
153 END DO
154 IF(igsti<=1)THEN
155 DO i=1,jlt
156 l = cand_e(i)
157 ni = cand_n(i)
158 IF(ni<=nsn)THEN
159 stif(i)=stf(l)*abs(stfn(ni))
160 ELSE
161 nn = ni - nsn
162 stif(i)=stf(l)*abs(stifi(nin)%P(nn))
163 END IF
164 ENDDO
165 ELSEIF(igsti==2)THEN
166 DO i=1,jlt
167 l = cand_e(i)
168 ni = cand_n(i)
169 IF(ni<=nsn)THEN
170 stif(i)=abs(stfn(ni))
171 ELSE
172 nn = ni - nsn
173 stif(i)=abs(stifi(nin)%P(nn))
174 END IF
175 stif(i)=half*(stf(l)+stif(i))
176 stif(i)=max(kmin,min(stif(i),kmax))
177 ENDDO
178 ELSEIF(igsti==3)THEN
179 DO i=1,jlt
180 l = cand_e(i)
181 ni = cand_n(i)
182 IF(ni<=nsn)THEN
183 stif(i)=abs(stfn(ni))
184 ELSE
185 nn = ni - nsn
186 stif(i)=abs(stifi(nin)%P(nn))
187 END IF
188 stif(i)=max(stf(l),stif(i))
189 stif(i)=max(kmin,min(stif(i),kmax))
190 ENDDO
191 ELSEIF(igsti==4)THEN
192 DO i=1,jlt
193 l = cand_e(i)
194 ni = cand_n(i)
195 IF(ni<=nsn)THEN
196 stif(i)=abs(stfn(ni))
197 ELSE
198 nn = ni - nsn
199 stif(i)=abs(stifi(nin)%P(nn))
200 END IF
201 stif(i)=min(stf(l),stif(i))
202 stif(i)=max(kmin,min(stif(i),kmax))
203 ENDDO
204 ELSEIF(igsti==5)THEN
205 DO i=1,jlt
206 l = cand_e(i)
207 ni = cand_n(i)
208 IF(ni<=nsn)THEN
209 stif(i)=abs(stfn(ni))
210 ELSE
211 nn = ni - nsn
212 stif(i)=abs(stifi(nin)%P(nn))
213 END IF
214 stif(i)=stf(l)*stif(i)/
215 . max(em30,(stf(l)+stif(i)))
216 stif(i)=max(kmin,min(stif(i),kmax))
217 ENDDO
218 ENDIF
219 ELSE
220C type7 KINET en -
221 DO i=1,jlt
222 ni = cand_n(i)
223 l = cand_e(i)
224 IF(ni<=nsn)THEN
225 ig = nsv(ni)
226 nsvg(i) = ig
227C KINI(I) = KINET(IG)
228 xi(i) = x(1,ig)
229 yi(i) = x(2,ig)
230 zi(i) = x(3,ig)
231 vxi(i) = v(1,ig)
232 vyi(i) = v(2,ig)
233 vzi(i) = v(3,ig)
234 msi(i)= ms(ig)
235 stif(i)=stf(l)*abs(stfn(ni))
236 ELSE
237 nn = ni - nsn
238 nsvg(i) = -nn
239C KINI(I) = KINFI(NIN)%P(NN)
240 xi(i) = xfi(nin)%P(1,nn)
241 yi(i) = xfi(nin)%P(2,nn)
242 zi(i) = xfi(nin)%P(3,nn)
243 vxi(i)= vfi(nin)%P(1,nn)
244 vyi(i)= vfi(nin)%P(2,nn)
245 vzi(i)= vfi(nin)%P(3,nn)
246 msi(i)= msfi(nin)%P(nn)
247 stif(i)=stf(l)*abs(stifi(nin)%P(nn))
248C
249 END IF
250C
251 ix=irect(1,l)
252 ix1(i)=ix
253 x1(i)=x(1,ix)
254 y1(i)=x(2,ix)
255 z1(i)=x(3,ix)
256C
257 ix=irect(2,l)
258 ix2(i)=ix
259 x2(i)=x(1,ix)
260 y2(i)=x(2,ix)
261 z2(i)=x(3,ix)
262C
263 ix=irect(3,l)
264 ix3(i)=ix
265 x3(i)=x(1,ix)
266 y3(i)=x(2,ix)
267 z3(i)=x(3,ix)
268C
269 ix=irect(4,l)
270 ix4(i)=ix
271 x4(i)=x(1,ix)
272 y4(i)=x(2,ix)
273 z4(i)=x(3,ix)
274C
275 END DO
276 END IF
277C
278 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
tri7box::vfi
type(real_pointer2), dimension(:), allocatable vfi
Definition tri7box.F:459
tri7box::stifi
type(real_pointer), dimension(:), allocatable stifi
Definition tri7box.F:449
tri7box::gapfi
type(real_pointer), dimension(:), allocatable gapfi
Definition tri7box.F:449
tri7box::msfi
type(real_pointer), dimension(:), allocatable msfi
Definition tri7box.F:449
tri7box::xfi
type(real_pointer2), dimension(:), allocatable xfi
Definition tri7box.F:459

◆ i7corkp3()

subroutine i7corkp3 ( integer jlt,
xi,
yi,
zi,
d,
dxi,
dyi,
dzi,
integer, dimension(mvsiz) nsvg,
integer nin,
integer iupd )

Definition at line 372 of file i7cork3.F.

376C-----------------------------------------------
377C M o d u l e s
378C-----------------------------------------------
379 USE imp_intm
380C-----------------------------------------------
381C I m p l i c i t T y p e s
382C-----------------------------------------------
383#include "implicit_f.inc"
384C-----------------------------------------------
385C G l o b a l P a r a m e t e r s
386C-----------------------------------------------
387#include "mvsiz_p.inc"
388C-----------------------------------------------
389C D u m m y A r g u m e n t s
390C-----------------------------------------------
391 INTEGER
392 . JLT,NIN,IUPD,NSVG(MVSIZ)
393C REAL
394C REAL
395 my_real
396 . xi(mvsiz), yi(mvsiz), zi(mvsiz),
397 . dxi(mvsiz),dyi(mvsiz),dzi(mvsiz),d(3,*)
398C-----------------------------------------------
399C L o c a l V a r i a b l e s
400C-----------------------------------------------
401 INTEGER I ,J ,NI, L, NN, NS ,NE
402C-----------------------------------------------
403C
404C
405 DO i=1,jlt
406 ni = nsvg(i)
407 IF(ni<0)THEN
408 nn = -ni
409 ns=ind_int(nin)%P(nn)
410C----------displacement
411 dxi(i)=dfi(1,ns)
412 dyi(i)=dfi(2,ns)
413 dzi(i)=dfi(3,ns)
414 ELSE
415 dxi(i) = d(1,ni)
416 dyi(i) = d(2,ni)
417 dzi(i) = d(3,ni)
418 END IF
419C----------update
420 IF(iupd>0.AND.ni<0)THEN
421 xi(i) = xi(i) + dxi(i)
422 yi(i) = yi(2) + dyi(i)
423 zi(i) = zi(3) + dzi(i)
424 END IF
425C
426 END DO
427C
428 RETURN
imp_intm
Definition imp_intm.F:118
imp_intm::ind_int
type(int_pointer2), dimension(:), allocatable ind_int
Definition imp_intm.F:133

◆ i7corp3()

subroutine i7corp3 ( integer jlt,
x,
integer, dimension(4,*) irect,
integer, dimension(*) cand_e,
integer, dimension(*) cand_n,
stif,
h1,
h2,
h3,
h4,
n1,
n2,
n3,
integer, dimension(mvsiz) ix1,
integer, dimension(mvsiz) ix2,
integer, dimension(mvsiz) ix3,
integer, dimension(mvsiz) ix4,
integer, dimension(mvsiz) nsvg,
vxi,
vyi,
vzi,
msi,
dxi,
dyi,
dzi,
integer nsn,
integer nin,
integer jlt_new,
integer lrem )

Definition at line 288 of file i7cork3.F.

294C-----------------------------------------------
295C M o d u l e s
296C-----------------------------------------------
297 USE tri7box
298 USE imp_intm
299C-----------------------------------------------
300C I m p l i c i t T y p e s
301C-----------------------------------------------
302#include "implicit_f.inc"
303C-----------------------------------------------
304C G l o b a l P a r a m e t e r s
305C-----------------------------------------------
306#include "mvsiz_p.inc"
307C-----------------------------------------------
308C D u m m y A r g u m e n t s
309C-----------------------------------------------
310 INTEGER IRECT(4,*), CAND_E(*), CAND_N(*),
311 . JLT,JLT_NEW,NIN,NSN,LREM
312 INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
313 . NSVG(MVSIZ)
314C REAL
315C REAL
316 my_real
317 . vxi(mvsiz), vyi(mvsiz), vzi(mvsiz), msi(mvsiz),
318 . x(3,*), stif(*),n1(mvsiz), n2(mvsiz), n3(mvsiz),
319 . h1(mvsiz), h2(mvsiz), h3(mvsiz), h4(mvsiz),
320 . dxi(mvsiz),dyi(mvsiz),dzi(mvsiz)
321C-----------------------------------------------
322C L o c a l V a r i a b l e s
323C-----------------------------------------------
324 INTEGER I ,J ,NI, L, NN, NS ,NE
325C-----------------------------------------------
326C
327C
328 DO i=1,jlt
329 ni = cand_n(i)
330 l = cand_e(i)
331 IF(ni>nsn)THEN
332 nn = ni - nsn
333 jlt_new = jlt_new + 1
334 vxi(jlt_new)= vfi(nin)%P(1,nn)
335 vyi(jlt_new)= vfi(nin)%P(2,nn)
336 vzi(jlt_new)= vfi(nin)%P(3,nn)
337 msi(jlt_new)= msfi(nin)%P(nn)
338 ix1(jlt_new) = irect(1,l)
339 ix2(jlt_new) = irect(2,l)
340 ix3(jlt_new) = irect(3,l)
341 ix4(jlt_new) = irect(4,l)
342 ns=ind_int(nin)%P(nn)
343 ne=shf_int(nin)+jlt_new+lrem
344 nsvg(jlt_new) = ns
345 h1(jlt_new) = h_e(1,ne)
346 h2(jlt_new) = h_e(2,ne)
347 h3(jlt_new) = h_e(3,ne)
348 h4(jlt_new) = h_e(4,ne)
349 stif(jlt_new)=stifs(ne)
350 n1(jlt_new)=n_e(1,ne)
351 n2(jlt_new)=n_e(2,ne)
352 n3(jlt_new)=n_e(3,ne)
353C----------displacement
354 dxi(jlt_new)=dfi(1,ns)
355 dyi(jlt_new)=dfi(2,ns)
356 dzi(jlt_new)=dfi(3,ns)
357C
358 END IF
359C
360 END DO
361C
362 RETURN
imp_intm::shf_int
integer, dimension(:), allocatable shf_int
Definition imp_intm.F:136