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

Go to the source code of this file.

Functions/Subroutines

subroutine s16deri3 (ngl, off, r, s, t, w, dnidr, dnids, dnidt, dxdr, dydr, dzdr, dxds, dyds, dzds, dxdt, dydt, dzdt, xx, yy, zz, px, py, pz, vol, deltax, kxx, ni, volg, ul, ir, is, it, voldp, nel)

Function/Subroutine Documentation

◆ s16deri3()

subroutine s16deri3 ( integer, dimension(*) ngl,
off,
r,
s,
t,
w,
dnidr,
dnids,
dnidt,
dxdr,
dydr,
dzdr,
dxds,
dyds,
dzds,
dxdt,
dydt,
dzdt,
xx,
yy,
zz,
px,
py,
pz,
vol,
deltax,
kxx,
ni,
volg,
ul,
integer ir,
integer is,
integer it,
double precision, dimension(*) voldp,
integer, intent(in) nel )

Definition at line 33 of file s16deri3.F.

43C-----------------------------------------------
44C M o d u l e s
45C-----------------------------------------------
46 USE message_mod
47C-----------------------------------------------
48C I m p l i c i t T y p e s
49C-----------------------------------------------
50#include "implicit_f.inc"
51#include "comlock.inc"
52C-----------------------------------------------
53C G l o b a l P a r a m e t e r s
54C-----------------------------------------------
55#include "mvsiz_p.inc"
56C-----------------------------------------------
57C D u m m y A r g u m e n t s
58C-----------------------------------------------
59 INTEGER, INTENT(IN) :: NEL
60 INTEGER NGL(*),IR,IS,IT
61C REAL
62 my_real
63 . r,s,t,w,
64 . dnidr(16),dnids(16),dnidt(16),dxdr(*),dydr(*),dzdr(*),
65 . dxds(*),dyds(*),dzds(*),dxdt(*),dydt(*),dzdt(*),
66 . drdx(mvsiz), dsdx(mvsiz), dtdx(mvsiz),
67 . drdy(mvsiz), dsdy(mvsiz), dtdy(mvsiz),
68 . drdz(mvsiz), dsdz(mvsiz), dtdz(mvsiz),
69 . xx(mvsiz,16),yy(mvsiz,16),zz(mvsiz,16),
70 . px(mvsiz,16),py(mvsiz,16),pz(mvsiz,16),
71 . vol(*),off(*),deltax(*),volg(mvsiz) ,
72 . kxx(mvsiz,16),ni(16),ul(mvsiz,16)
73 DOUBLE PRECISION
74 . VOLDP(*)
75C-----------------------------------------------
76C C o m m o n B l o c k s
77C-----------------------------------------------
78C-----------------------------------------------
79C L o c a l V a r i a b l e s
80C-----------------------------------------------
81 INTEGER I,N,ICOR
82 my_real
83 . d, aa, bb, det(mvsiz),r9 ,r13 ,s9 ,s10 ,s11 ,s12 ,t10 ,t14
84 DOUBLE PRECISION
85 . DETDP
86C-----------------------------------------------
87C
88C
89C
90C ^ S _ T
91C | /|
92C | /
93C 7 | 14 /
94C O-----------|-----O-----------------O 6
95C /. | / /|
96C / . | / / |
97C / . | / |
98C / . | / / |
99C / . | / / |
100C 15 O . * + * * O |
101C / . | / / 13 |
102C / . / / |
103C / | + / |
104C / .# # / # / |
105C / . * 16 | * / * |
106C 8 O-----------------O-----------------O 5 |
107C | . | / | |
108C | @ . @ @ | |
109C R <------|- - -+ - - -#- - - - - # - - - - -#|- - -+ |
110C | . * / * | * |
111C | . 3 | 10 | |
112C | O......../.. .....O.....|...........O 2
113C | @ ' @ | @ | /
114C | ' # / # #| /
115C | ' + | | /
116C | ' | /
117C | 11 ' | | /
118C | O @ @ + @ | O
119C | ' | / 9
120C | ' | /
121C | ' | /
122C | ' | /
123C |' |/
124C O-----------------O-----------------O
125C 4 12 1
126C
127C
128C
129C
130C-----------------------------------------------
131C
132C
133C ^ S _ T
134C | /|
135C | /
136C | /
137C ( 7)==========|===(14)===============( 6)
138C //| | / //|
139C // | | / //||
140C // | | // ||
141C // | | / // ||
142C // | | / // ||
143C (15) | *- - - + - * - - - - -* (13) ||
144C // | /| | /| / /| // ||
145C // | / // ||
146C // | / | |/ | + / // ||
147C // |#- - - - - # - -/- - -# // ||
148C // | * - - /|- -*- - -/ -//* ||
149C ( 8)===============(16)==============( 5) ||
150C || / || / | / | / | / ||| | ||
151C || @- | - - - @ - - - - -@ || ||
152C R <-----||- - -+ -|- -# - - -| - # - - -|- -#|| - -+ ||
153C || | * - - /| - *- - -/-|| * ||
154C || | || / | / | |||/ ||
155C || ( 3)-------/--|---(10)----||---------( 2)
156C || @ /- / - - @ - -/ - - @ || //
157C || |/ #- - -/| - # - - -|- -#|| //
158C || / / + /| / || //
159C || /| | | || //
160C || / / / | / || //
161C || (11) @- - - - - @ - + - - -@ || ( 9)
162C || / || //
163C || / || //
164C || / || //
165C || / ||//
166C ||/ ||/
167C ( 4)===============(12)==============( 1)
168C
169C
170C
171C
172C
173C*/
174C-----------------------------------------------
175C
176C _
177C \
178C f = /_ (fi * Ni)
179C _
180C \
181C df/dx = /_ (fi * dNi/dx)
182C
183C dNi/dx = dNi/dr dr/dx + dNi/ds ds/dx + dNi/dt dt/dx
184C
185C-----------------------------------------------
186C _
187C \
188C x(r,s,t) = /_ (xi * Ni(r,s,t))
189C _
190C \
191C y(r,s,t) = /_ (yi * Ni(r,s,t))
192C _
193C \
194C z(r,s,t) = /_ zi * Ni(r,s,t))
195C
196C _
197C \
198C dx/dr = /_ (xi * dNi/dr)
199C ...
200C
201C [dx/dr dy/dr dz/dr]
202C [J] = |dx/ds dy/ds dz/ds|
203C [dx/dt dy/dt dz/dt]
204C
205C |dNi/dx| -1 |dNi/dr|
206C {dNi/dy} = [J] {dNi/ds}
207C |dNi/dz| |dNi/dt|
208C
209C
210C [dNi/dx 0 0 ]
211C | 0 dNi/dy 0 |
212C [Bi] = | 0 0 dNi/dz|
213C |dNi/dy dNi/dx 0 |
214C | 0 dNi/dz dNi/dy|
215C [dNi/dz 0 dNi/dx]
216C
217C
218C [dNi/dx 0 0 dNi/dy 0 dNi/dz]
219C [Bi]t = | 0 dNi/dy 0 dNi/dx dNi/dz 0 |
220C [ 0 0 dNi/dz 0 dNi/dy dNi/dx]
221C
222C [D11 D12 D12 0 0 0 ]
223C |D12 D11 D12 0 0 0 |
224C [D] = |D12 D12 D11 0 0 0 |
225C | 0 0 0 G 0 0 |
226C | 0 0 0 0 G 0 |
227C [ 0 0 0 0 0 G ]
228C
229C [D11*dNi/dx D12*dNi/dx D12*dNi/dx G*dNi/dy 0 G*dNi/dz]
230C [Bi]t [D] = |D12*dNi/dy D11*dNi/dy D12*dNi/dy G*dNi/dx G*dNi/dz 0 |
231C [D12*dNi/dz D12*dNi/dz D11*dNi/dz 0 G*dNi/dy G*dNi/dx]
232C _
233C / t
234C eps = [B] u sig = [D] eps F = _/ [B] sig dvol
235C
236C _
237C / t
238C [Kij] = _/ [Bi] [D] [Bj] dvol
239C
240C
241C
242C [ D11*dNi/dx*dNj/dx D12*dNi/dx*dNj/dy D12*dNi/dx*dNj/dz ]
243C | + G*dNi/dy*dNj/dy + G*dNi/dy*dNj/dx + G*dNi/dz*dNj/dx |
244C | + G*dNi/dz*dNj/dz |
245C [Kij] | D11*dNi/dy*dNj/dy D12*dNi/dy*dNj/dz |
246C ---- = | D12*dNi/dy*dNj/dx + G*dNi/dx*dNj/dx + G*dNi/dz*dNj/dy |
247C Volp | + G*dNi/dx*dNj/dy + G*dNi/dz*dNj/dz |
248C | D11*dNi/dz*dNj/dz |
249C | D12*dNi/dz*dNj/dx D12*dNi/dz*dNj/dy + G*dNi/dy*dNj/dy |
250C [ + G*dNi/dx*dNj/dz + G*dNi/dy*dNj/dz + G*dNi/dx*dNj/dx ]
251C
252C-----------------------------------------------
253C-----------------------------------------------------------------------
254C dx/dr; dx/ds; dx/dt
255C-----------------------------------------------------------------------
256 DO i=1,nel
257C
258 dxdr(i) = dnidr(1)*xx(i,1) + dnidr(2)*xx(i,2) + dnidr(3)*xx(i,3)
259 + + dnidr(4)*xx(i,4) + dnidr(5)*xx(i,5) + dnidr(6)*xx(i,6)
260 + + dnidr(7)*xx(i,7) + dnidr(8)*xx(i,8)
261 + + dnidr(9)*(xx(i,9) - xx(i,11)) + dnidr(10)*xx(i,10)
262 + + dnidr(12)*xx(i,12) + dnidr(13)*(xx(i,13) - xx(i,15))
263 + + dnidr(14)*xx(i,14) + dnidr(16)*xx(i,16)
264C
265 dxds(i) = dnids(1)*xx(i,1) + dnids(2)*xx(i,2) + dnids(3)*xx(i,3)
266 + + dnids(4)*xx(i,4) + dnids(5)*xx(i,5) + dnids(6)*xx(i,6)
267 + + dnids(7)*xx(i,7) + dnids(8)*xx(i,8)
268 + + dnids(9)* (xx(i,9) - xx(i,13))
269 + + dnids(10)*(xx(i,10) - xx(i,14))
270 + + dnids(11)*(xx(i,11) - xx(i,15))
271 + + dnids(12)*(xx(i,12) - xx(i,16))
272C
273 dxdt(i) = dnidt(1)*xx(i,1) + dnidt(2)*xx(i,2) + dnidt(3)*xx(i,3)
274 + + dnidt(4)*xx(i,4) + dnidt(5)*xx(i,5) + dnidt(6)*xx(i,6)
275 + + dnidt(7)*xx(i,7) + dnidt(8)*xx(i,8)
276 + + dnidt(9)*xx(i,9) + dnidt(10)*(xx(i,10) - xx(i,12))
277 + + dnidt(11)*xx(i,11) + dnidt(13)*xx(i,13)
278 + + dnidt(14)*(xx(i,14) - xx(i,16)) + dnidt(15)*xx(i,15)
279C-----------------------------------------------------------------------
280C dy/dr; dy/ds; dy/dt
281C-----------------------------------------------------------------------
282 dydr(i) = dnidr(1)*yy(i,1) + dnidr(2)*yy(i,2) + dnidr(3)*yy(i,3)
283 + + dnidr(4)*yy(i,4) + dnidr(5)*yy(i,5) + dnidr(6)*yy(i,6)
284 + + dnidr(7)*yy(i,7) + dnidr(8)*yy(i,8)
285 + + dnidr(9)*(yy(i,9) - yy(i,11)) + dnidr(10)*yy(i,10)
286 + + dnidr(12)*yy(i,12) + dnidr(13)*(yy(i,13) - yy(i,15))
287 + + dnidr(14)*yy(i,14) + dnidr(16)*yy(i,16)
288C
289 dyds(i) = dnids(1)*yy(i,1) + dnids(2)*yy(i,2) + dnids(3)*yy(i,3)
290 + + dnids(4)*yy(i,4) + dnids(5)*yy(i,5) + dnids(6)*yy(i,6)
291 + + dnids(7)*yy(i,7) + dnids(8)*yy(i,8)
292 + + dnids(9)* (yy(i,9) - yy(i,13))
293 + + dnids(10)*(yy(i,10) - yy(i,14))
294 + + dnids(11)*(yy(i,11) - yy(i,15))
295 + + dnids(12)*(yy(i,12) - yy(i,16))
296C
297 dydt(i) = dnidt(1)*yy(i,1) + dnidt(2)*yy(i,2) + dnidt(3)*yy(i,3)
298 + + dnidt(4)*yy(i,4) + dnidt(5)*yy(i,5) + dnidt(6)*yy(i,6)
299 + + dnidt(7)*yy(i,7) + dnidt(8)*yy(i,8)
300 + + dnidt(9)*yy(i,9) + dnidt(10)*(yy(i,10) - yy(i,12))
301 + + dnidt(11)*yy(i,11) + dnidt(13)*yy(i,13)
302 + + dnidt(14)*(yy(i,14) - yy(i,16)) + dnidt(15)*yy(i,15)
303C-----------------------------------------------------------------------
304C dz/dr; dz/ds; dz/dt
305C-----------------------------------------------------------------------
306 dzdr(i) = dnidr(1)*zz(i,1) + dnidr(2)*zz(i,2) + dnidr(3)*zz(i,3)
307 + + dnidr(4)*zz(i,4) + dnidr(5)*zz(i,5) + dnidr(6)*zz(i,6)
308 + + dnidr(7)*zz(i,7) + dnidr(8)*zz(i,8)
309 + + dnidr(9)*(zz(i,9) - zz(i,11)) + dnidr(10)*zz(i,10)
310 + + dnidr(12)*zz(i,12) + dnidr(13)*(zz(i,13) - zz(i,15))
311 + + dnidr(14)*zz(i,14) + dnidr(16)*zz(i,16)
312C
313 dzds(i) = dnids(1)*zz(i,1) + dnids(2)*zz(i,2) + dnids(3)*zz(i,3)
314 + + dnids(4)*zz(i,4) + dnids(5)*zz(i,5) + dnids(6)*zz(i,6)
315 + + dnids(7)*zz(i,7) + dnids(8)*zz(i,8)
316 + + dnids(9)* (zz(i,9) - zz(i,13))
317 + + dnids(10)*(zz(i,10) - zz(i,14))
318 + + dnids(11)*(zz(i,11) - zz(i,15))
319 + + dnids(12)*(zz(i,12) - zz(i,16))
320C
321 dzdt(i) = dnidt(1)*zz(i,1) + dnidt(2)*zz(i,2) + dnidt(3)*zz(i,3)
322 + + dnidt(4)*zz(i,4) + dnidt(5)*zz(i,5) + dnidt(6)*zz(i,6)
323 + + dnidt(7)*zz(i,7) + dnidt(8)*zz(i,8)
324 + + dnidt(9)*zz(i,9) + dnidt(10)*(zz(i,10) - zz(i,12))
325 + + dnidt(11)*zz(i,11) + dnidt(13)*zz(i,13)
326 + + dnidt(14)*(zz(i,14) - zz(i,16)) + dnidt(15)*zz(i,15)
327c print *,'dx y z/dr'
328c print *,DXDR(1),DYDR(1),DZDR(1)
329c print *,'dx y z/ds'
330c print *,DXDS(1),DYDS(1),DZDS(1)
331c print *,'dx y z/dt='
332c print *,DXDT(1),DYDT(1),DZDT(1)
333C-----------------------------------------------------------------------
334C -1
335C [J] Inversion du jacobien
336C-----------------------------------------------------------------------
337 drdx(i)=dyds(i)*dzdt(i)-dzds(i)*dydt(i)
338 drdy(i)=dzds(i)*dxdt(i)-dxds(i)*dzdt(i)
339 drdz(i)=dxds(i)*dydt(i)-dyds(i)*dxdt(i)
340C
341 dsdz(i)=dxdt(i)*dydr(i)-dydt(i)*dxdr(i)
342 dsdy(i)=dzdt(i)*dxdr(i)-dxdt(i)*dzdr(i)
343 dsdx(i)=dydt(i)*dzdr(i)-dzdt(i)*dydr(i)
344C
345 dtdx(i)=dydr(i)*dzds(i)-dzdr(i)*dyds(i)
346 dtdy(i)=dzdr(i)*dxds(i)-dxdr(i)*dzds(i)
347 dtdz(i)=dxdr(i)*dyds(i)-dydr(i)*dxds(i)
348C
349 detdp = dxdr(i) * drdx(i)
350 . + dydr(i) * drdy(i)
351 . + dzdr(i) * drdz(i)
352 det(i) = detdp
353 voldp(i) = w * detdp
354 vol(i) = voldp(i)
355C
356 ENDDO
357C
358C IF(IDTMIN(1)==1)THEN
359C ...
360C ENDIF
361C
362C
363 icor=0
364 DO i=1,nel
365 IF(off(i)==zero)THEN
366 vol(i)=1.
367 voldp(i) = one
368 ELSEIF(vol(i)<=zero)THEN
369 icor=1
370 ENDIF
371 ENDDO
372 IF(icor/=0)THEN
373 DO i=1,nel
374 IF(vol(i)<=zero)THEN
375 CALL ancmsg(msgid=170,anmode=aninfo,
376 . i1=ngl(i),i2=ir,i3=is,i4=it)
377 CALL arret(2)
378 ENDIF
379 ENDDO
380 ENDIF
381C
382 DO i=1,nel
383C-----------------------------------------------------------------------
384C -1 Inversion du jacobien suite
385C [J] et repere local r,s,t
386C-----------------------------------------------------------------------
387 d = one/det(i)
388 drdx(i)=d*drdx(i)
389 dsdx(i)=d*dsdx(i)
390 dtdx(i)=d*dtdx(i)
391C
392 drdy(i)=d*drdy(i)
393 dsdy(i)=d*dsdy(i)
394 dtdy(i)=d*dtdy(i)
395C
396 drdz(i)=d*drdz(i)
397 dsdz(i)=d*dsdz(i)
398 dtdz(i)=d*dtdz(i)
399C-----------------------------------------------------------------------
400C |dNi/dx| -1 |dNi/dr|
401C {dNi/dy} = [J] {dNi/ds}
402C |dNi/dz| |dNi/dt|
403C-----------------------------------------------------------------------
404 px(i,1) = dnidr(1)*drdx(i) + dnids(1)*dsdx(i) + dnidt(1)*dtdx(i)
405 px(i,2) = dnidr(2)*drdx(i) + dnids(2)*dsdx(i) + dnidt(2)*dtdx(i)
406 px(i,3) = dnidr(3)*drdx(i) + dnids(3)*dsdx(i) + dnidt(3)*dtdx(i)
407 px(i,4) = dnidr(4)*drdx(i) + dnids(4)*dsdx(i) + dnidt(4)*dtdx(i)
408 px(i,5) = dnidr(5)*drdx(i) + dnids(5)*dsdx(i) + dnidt(5)*dtdx(i)
409 px(i,6) = dnidr(6)*drdx(i) + dnids(6)*dsdx(i) + dnidt(6)*dtdx(i)
410 px(i,7) = dnidr(7)*drdx(i) + dnids(7)*dsdx(i) + dnidt(7)*dtdx(i)
411 px(i,8) = dnidr(8)*drdx(i) + dnids(8)*dsdx(i) + dnidt(8)*dtdx(i)
412 r9 = dnidr(9) *drdx(i)
413 r13 = dnidr(13)*drdx(i)
414 s9 = dnids(9) *dsdx(i)
415 s10 = dnids(10)*dsdx(i)
416 s11 = dnids(11)*dsdx(i)
417 s12 = dnids(12)*dsdx(i)
418 t10 = dnidt(10)*dtdx(i)
419 t14 = dnidt(14)*dtdx(i)
420 px(i,9) = r9 + s9 + dnidt(9)*dtdx(i)
421 px(i,10)= dnidr(10)*drdx(i)+ s10 + t10
422 px(i,11)=-r9 + s11 + dnidt(11)*dtdx(i)
423 px(i,12)= dnidr(12)*drdx(i)+ s12 - t10
424 px(i,13)= r13 - s9 + dnidt(13)*dtdx(i)
425 px(i,14)= dnidr(14)*drdx(i)- s10 + t14
426 px(i,15)=-r13 - s11 + dnidt(15)*dtdx(i)
427 px(i,16)= dnidr(16)*drdx(i)- s12 - t14
428C
429 py(i,1) = dnidr(1)*drdy(i) + dnids(1)*dsdy(i) + dnidt(1)*dtdy(i)
430 pz(i,1) = dnidr(1)*drdz(i) + dnids(1)*dsdz(i) + dnidt(1)*dtdz(i)
431 py(i,2) = dnidr(2)*drdy(i) + dnids(2)*dsdy(i) + dnidt(2)*dtdy(i)
432 py(i,3) = dnidr(3)*drdy(i) + dnids(3)*dsdy(i) + dnidt(3)*dtdy(i)
433 py(i,4) = dnidr(4)*drdy(i) + dnids(4)*dsdy(i) + dnidt(4)*dtdy(i)
434 py(i,5) = dnidr(5)*drdy(i) + dnids(5)*dsdy(i) + dnidt(5)*dtdy(i)
435 py(i,6) = dnidr(6)*drdy(i) + dnids(6)*dsdy(i) + dnidt(6)*dtdy(i)
436 py(i,7) = dnidr(7)*drdy(i) + dnids(7)*dsdy(i) + dnidt(7)*dtdy(i)
437 py(i,8) = dnidr(8)*drdy(i) + dnids(8)*dsdy(i) + dnidt(8)*dtdy(i)
438 r9 = dnidr(9) *drdy(i)
439 r13 = dnidr(13)*drdy(i)
440 s9 = dnids(9) *dsdy(i)
441 s10 = dnids(10)*dsdy(i)
442 s11 = dnids(11)*dsdy(i)
443 s12 = dnids(12)*dsdy(i)
444 t10 = dnidt(10)*dtdy(i)
445 t14 = dnidt(14)*dtdy(i)
446 py(i,9) = r9 + s9 + dnidt(9)*dtdy(i)
447 py(i,10)= dnidr(10)*drdy(i)+ s10 + t10
448 py(i,11)=-r9 + s11 + dnidt(11)*dtdy(i)
449 py(i,12)= dnidr(12)*drdy(i)+ s12 - t10
450 py(i,13)= r13 - s9 + dnidt(13)*dtdy(i)
451 py(i,14)= dnidr(14)*drdy(i)- s10 + t14
452 py(i,15)=-r13 - s11 + dnidt(15)*dtdy(i)
453 py(i,16)= dnidr(16)*drdy(i)- s12 - t14
454C
455 pz(i,1) = dnidr(1)*drdz(i) + dnids(1)*dsdz(i) + dnidt(1)*dtdz(i)
456 pz(i,2) = dnidr(2)*drdz(i) + dnids(2)*dsdz(i) + dnidt(2)*dtdz(i)
457 pz(i,3) = dnidr(3)*drdz(i) + dnids(3)*dsdz(i) + dnidt(3)*dtdz(i)
458 pz(i,4) = dnidr(4)*drdz(i) + dnids(4)*dsdz(i) + dnidt(4)*dtdz(i)
459 pz(i,5) = dnidr(5)*drdz(i) + dnids(5)*dsdz(i) + dnidt(5)*dtdz(i)
460 pz(i,6) = dnidr(6)*drdz(i) + dnids(6)*dsdz(i) + dnidt(6)*dtdz(i)
461 pz(i,7) = dnidr(7)*drdz(i) + dnids(7)*dsdz(i) + dnidt(7)*dtdz(i)
462 pz(i,8) = dnidr(8)*drdz(i) + dnids(8)*dsdz(i) + dnidt(8)*dtdz(i)
463 r9 = dnidr(9) *drdz(i)
464 r13 = dnidr(13)*drdz(i)
465 s9 = dnids(9) *dsdz(i)
466 s10 = dnids(10)*dsdz(i)
467 s11 = dnids(11)*dsdz(i)
468 s12 = dnids(12)*dsdz(i)
469 t10 = dnidt(10)*dtdz(i)
470 t14 = dnidt(14)*dtdz(i)
471 pz(i,9) = r9 + s9 + dnidt(9)*dtdz(i)
472 pz(i,10)= dnidr(10)*drdz(i)+ s10 + t10
473 pz(i,11)=-r9 + s11 + dnidt(11)*dtdz(i)
474 pz(i,12)= dnidr(12)*drdz(i)+ s12 - t10
475 pz(i,13)= r13 - s9 + dnidt(13)*dtdz(i)
476 pz(i,14)= dnidr(14)*drdz(i)- s10 + t14
477 pz(i,15)=-r13 - s11 + dnidt(15)*dtdz(i)
478 pz(i,16)= dnidr(16)*drdz(i)- s12 - t14
479C
480 ENDDO
481C
482C
483C
484 DO n=1,16
485 DO i=1,nel
486 kxx(i,n) = vol(i)*
487 . (px(i,n)*px(i,n) + py(i,n)*py(i,n) + pz(i,n)*pz(i,n))
488 ul(i,n) = ul(i,n) + kxx(i,n)
489 ENDDO
490 ENDDO
491 DO i=1,nel
492 volg(i) =volg(i) + vol(i)
493 ENDDO
494C
495c AA = 0
496c DO N=1,16
497c AA = AA + NI(N)*NI(N)
498c ENDDO
499c DO I=1,NEL
500c BB = 1.E20
501c DO N=1,16
502c BB = MIN(BB,NI(N)*NI(N)/UL(I,N))
503c ENDDO
504c DELTAX(I) = SQRT(2.*BB/AA)
505c ENDDO
506C
507C
508 2000 FORMAT(/' ZERO OR NEGATIVE VOLUME : DELETE 3D-ELEMENT NB',i10/)
509C
510 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
message_mod
Definition message_mod.F:1249
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