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

Go to the source code of this file.

Functions/Subroutines

subroutine inelts (x, irect, ixs, nint, nel, i, area, noint, ir, surf_eltyp, surf_elem)
subroutine ineltc (nelc, neltg, is, surf_eltyp, surf_elem)
subroutine i20nelts (x, irect, ixs, nint, nel, i, area, noint, ir, surf_eltyp, surf_elem)
subroutine ineltigeo (xe, irect, ixs, nint, nelig3d, i, area, noint, ir, surf_eltyp_ige, ixig3d, kxig3d, igeo, surf_elem_ige)

Function/Subroutine Documentation

◆ i20nelts()

subroutine i20nelts ( x,
integer, dimension(4) irect,
integer, dimension(nixs,*) ixs,
integer nint,
integer nel,
integer i,
area,
integer noint,
integer ir,
integer, dimension(*) surf_eltyp,
integer, dimension(*) surf_elem )

Definition at line 165 of file inelt.F.

168C-----------------------------------------------
169C I m p l i c i t T y p e s
170C-----------------------------------------------
171#include "implicit_f.inc"
172C-----------------------------------------------
173C C o m m o n B l o c k s
174C-----------------------------------------------
175C-----------------------------------------------
176C D u m m y A r g u m e n t s
177C-----------------------------------------------
178 INTEGER NINT, NEL, I, NOINT,IR,SURF_ELTYP(*),SURF_ELEM(*)
179C REAL
180 my_real
181 . area
182 INTEGER IRECT(4), IXS(NIXS,*)
183C REAL
184 my_real
185 . x(3,*)
186C-----------------------------------------------
187C L o c a l V a r i a b l e s
188C-----------------------------------------------
189 INTEGER IY(4), N, JJ, II, K, NN, KK, IC, IAD,
190 . NUSER, NUSERM
191C REAL
192 my_real
193 . n1, n2, n3, dds
194 my_real xx1(4), xx2(4),xx3(4),xs1,ys1,zs1,xc,yc,zc
195C-----------------------------------------------
196C E x t e r n a l F u n c t i o n s
197C-----------------------------------------------
198C
199 ic =0
200 nel=0
201
202 IF (i <= 0) RETURN
203 IF (surf_eltyp(i) /= 1) RETURN
204C
205 nel=surf_elem(i)
206C-----------------------------------------------
207C VERIFICATION DE L'ORIENTATION DES SEGMENTS
208C-----------------------------------------------
209 xs1=0.
210 ys1=0.
211 zs1=0.
212 DO 100 jj=1,4
213 nn=irect(jj)
214 iy(jj)=nn
215 xx1(jj)=x(1,nn)
216 xx2(jj)=x(2,nn)
217 xx3(jj)=x(3,nn)
218 xs1=xs1+.25*x(1,nn)
219 ys1=ys1+.25*x(2,nn)
220 100 zs1=zs1+.25*x(3,nn)
221C
222 CALL norma1(n1,n2,n3,area,xx1,xx2,xx3)
223 xc=0.
224 yc=0.
225 zc=0.
226 DO 110 k=1,8
227 kk=ixs(k+1,nel)
228 xc=xc+x(1,kk)
229 yc=yc+x(2,kk)
230 zc=zc+x(3,kk)
231 110 CONTINUE
232 xc=xc*one_over_8
233 yc=yc*one_over_8
234 zc=zc*one_over_8
235 IF(ir/=0) RETURN
236 IF(ic>=2)RETURN
237 dds=n1*(xc-xs1)+n2*(yc-ys1)+n3*(zc-zs1)
238 IF(dds<0) RETURN
239 IF(iy(3)==iy(4)) THEN
240 irect(1)=iy(2)
241 irect(2)=iy(1)
242 ELSE
243 DO 120 kk=1,4
244 120 irect(kk)=iy(4-kk+1)
245 ENDIF
246 RETURN
247C
my_real
#define my_real
Definition cppsort.cpp:32
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
norma1
subroutine norma1(n1, n2, n3, area, xx1, xx2, xx3)
Definition norma1.F:38

◆ ineltc()

subroutine ineltc ( integer nelc,
integer neltg,
integer is,
integer, dimension(*) surf_eltyp,
integer, dimension(*) surf_elem )

Definition at line 131 of file inelt.F.

132C-----------------------------------------------
133C I m p l i c i t T y p e s
134C-----------------------------------------------
135#include "implicit_f.inc"
136C-----------------------------------------------
137C D u m m y A r g u m e n t s
138C-----------------------------------------------
139 INTEGER NELC, IS, NELTG,SURF_ELTYP(*),SURF_ELEM(*)
140C-----------------------------------------------
141C L o c a l V a r i a b l e s
142C-----------------------------------------------
143 INTEGER ELTYP
144C-----------------------------------------------
145 nelc=0
146 neltg=0
147
148 IF (is <= 0) RETURN
149 eltyp = surf_eltyp(is)
150 IF (eltyp == 3) THEN
151 nelc = surf_elem(is)
152 ELSEIF (eltyp == 7) THEN
153 neltg = surf_elem(is)
154 ENDIF
155!---
156 RETURN

◆ ineltigeo()

subroutine ineltigeo ( xe,
integer, dimension(4,*) irect,
integer, dimension(nixs,*) ixs,
integer nint,
integer nelig3d,
integer i,
area,
integer noint,
integer ir,
integer, dimension(*) surf_eltyp_ige,
integer, dimension(*) ixig3d,
integer, dimension(nixig3d,*) kxig3d,
integer, dimension(npropgi,*) igeo,
integer, dimension(*) surf_elem_ige )

Definition at line 256 of file inelt.F.

259C-----------------------------------------------
260C I m p l i c i t T y p e s
261C-----------------------------------------------
262#include "implicit_f.inc"
263C-----------------------------------------------
264C C o m m o n B l o c k s
265C-----------------------------------------------
266#include "param_c.inc"
267C-----------------------------------------------
268C D u m m y A r g u m e n t s
269C-----------------------------------------------
270 INTEGER NINT, NELIG3D, I, NOINT,IR
271C REAL
272 my_real
273 . area
274 INTEGER IRECT(4,*), IXS(NIXS,*),SURF_ELTYP_IGE(*),
275 . SURF_ELEM_IGE(*),KXIG3D(NIXIG3D,*),IGEO(NPROPGI,*),
276 . IXIG3D(*)
277C REAL
278 my_real
279 . xe(3,*)
280C-----------------------------------------------
281C L o c a l V a r i a b l e s
282C-----------------------------------------------
283 INTEGER IY(4), N, JJ, II, K, NN, KK, IC, IAD,
284 . NUSER, NUSERM, COIN_IGE(8), PX, PY, PZ, IPID
285C REAL
286 my_real
287 . n1, n2, n3, dds
288 my_real :: xx1(4), xx2(4),xx3(4),xc,yc,zc,xs1,ys1,zs1
289C-----------------------------------------------
290C E x t e r n a l F u n c t i o n s
291C-----------------------------------------------
292C
293 ic =0
294 nelig3d=0
295 IF (surf_eltyp_ige(i) /= 101) RETURN
296C
297 nelig3d=surf_elem_ige(i)
298
299C-----------------------------------------------
300C VERIFICATION DE L'ORIENTATION DES SEGMENTS
301C-----------------------------------------------
302 xs1=0.
303 ys1=0.
304 zs1=0.
305 DO jj=1,4
306 nn=irect(jj,i)
307 iy(jj)=nn
308 xx1(jj)=xe(1,nn)
309 xx2(jj)=xe(2,nn)
310 xx3(jj)=xe(3,nn)
311 xs1=xs1+.25*xe(1,nn)
312 ys1=ys1+.25*xe(2,nn)
313 zs1=zs1+.25*xe(3,nn)
314 ENDDO
315C
316 CALL norma1(n1,n2,n3,area,xx1,xx2,xx3)
317 xc=0.
318 yc=0.
319 zc=0.
320
321 ipid = kxig3d(2,nelig3d)
322 px = igeo(41,ipid)-1
323 py = igeo(42,ipid)-1
324 pz = igeo(43,ipid)-1
325 coin_ige(1) = (px+1)*py+1
326 coin_ige(2) = (px+1)*(py+1)
327 coin_ige(3) = px+1
328 coin_ige(4) = 1
329 coin_ige(5) = (px+1)*(py+1)*pz+(px+1)*py+1
330 coin_ige(6) = (px+1)*(py+1)*(pz+1)
331 coin_ige(7) = (px+1)*(py+1)*pz+px+1
332 coin_ige(8) = (px+1)*(py+1)*pz+1
333
334 DO k=1,8
335 xc=xc+xe(1,ixig3d(kxig3d(4,nelig3d)+coin_ige(k)-1))
336 yc=yc+xe(2,ixig3d(kxig3d(4,nelig3d)+coin_ige(k)-1))
337 zc=zc+xe(3,ixig3d(kxig3d(4,nelig3d)+coin_ige(k)-1))
338 ENDDO
339 xc=xc*one_over_8
340 yc=yc*one_over_8
341 zc=zc*one_over_8
342 IF(ir/=0) RETURN
343 IF(ic>=2)RETURN
344 dds=n1*(xc-xs1)+n2*(yc-ys1)+n3*(zc-zs1)
345 IF(dds<0) RETURN
346 DO kk=1,4
347 irect(kk,i)=iy(4-kk+1)
348 ENDDO
349 RETURN
350C

◆ inelts()

subroutine inelts ( x,
integer, dimension(4,*) irect,
integer, dimension(nixs,*) ixs,
integer nint,
integer nel,
integer i,
area,
integer noint,
integer ir,
integer, dimension(*) surf_eltyp,
integer, dimension(*) surf_elem )

Definition at line 36 of file inelt.F.

39C-----------------------------------------------
40C I m p l i c i t T y p e s
41C-----------------------------------------------
42#include "implicit_f.inc"
43C-----------------------------------------------
44C C o m m o n B l o c k s
45C-----------------------------------------------
46C-----------------------------------------------
47C D u m m y A r g u m e n t s
48C-----------------------------------------------
49 INTEGER NINT, NEL, I, NOINT,IR,SURF_ELTYP(*),SURF_ELEM(*)
50C REAL
51 my_real
52 . area
53 INTEGER IRECT(4,*), IXS(NIXS,*)
54C REAL
55 my_real
56 . x(3,*)
57C-----------------------------------------------
58C L o c a l V a r i a b l e s
59C-----------------------------------------------
60 INTEGER IY(4), N, JJ, II, K, NN, KK, IC, IAD,
61 . NUSER, NUSERM
62C REAL
63 my_real
64 . n1, n2, n3, dds
65 my_real :: xx1(4), xx2(4),xx3(4),xs1,ys1,zs1,xc,yc,zc
66C-----------------------------------------------
67C E x t e r n a l F u n c t i o n s
68C-----------------------------------------------
69C
70 ic =0
71 nel=0
72 IF (i <= 0) RETURN
73 IF (surf_eltyp(i) /= 1) RETURN
74C
75 nel=surf_elem(i)
76C-----------------------------------------------
77C VERIFICATION DE L'ORIENTATION DES SEGMENTS
78C-----------------------------------------------
79 xs1=0.
80 ys1=0.
81 zs1=0.
82 DO 100 jj=1,4
83 nn=irect(jj,i)
84 iy(jj)=nn
85 xx1(jj)=x(1,nn)
86 xx2(jj)=x(2,nn)
87 xx3(jj)=x(3,nn)
88 xs1=xs1+.25*x(1,nn)
89 ys1=ys1+.25*x(2,nn)
90 100 zs1=zs1+.25*x(3,nn)
91C
92 CALL norma1(n1,n2,n3,area,xx1,xx2,xx3)
93 xc=0.
94 yc=0.
95 zc=0.
96 DO 110 k=1,8
97 kk=ixs(k+1,nel)
98 xc=xc+x(1,kk)
99 yc=yc+x(2,kk)
100 zc=zc+x(3,kk)
101 110 CONTINUE
102 xc=xc*one_over_8
103 yc=yc*one_over_8
104 zc=zc*one_over_8
105 IF(ir/=0) RETURN
106 IF(ic>=2)RETURN
107 dds=n1*(xc-xs1)+n2*(yc-ys1)+n3*(zc-zs1)
108 IF(dds<0) RETURN
109 IF(iy(3)==iy(4)) THEN
110 irect(1,i)=iy(2)
111 irect(2,i)=iy(1)
112 ELSE
113 DO 120 kk=1,4
114 120 irect(kk,i)=iy(4-kk+1)
115 ENDIF
116 RETURN
117C