OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s8zderi3.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| s8zderi3 ../starter/source/elements/solid/solide8z/s8zderi3.F
25!||--- called by ------------------------------------------------------
26!|| s8cinit3 ../starter/source/elements/thickshell/solide8c/s8cinit3.F
27!|| s8zinit3 ../starter/source/elements/solid/solide8z/s8zinit3.F
28!||--- calls -----------------------------------------------------
29!|| ancmsg ../starter/source/output/message/message.F
30!||--- uses -----------------------------------------------------
31!|| message_mod ../starter/share/message_module/message_mod.F
32!||====================================================================
33 SUBROUTINE s8zderi3(VOL,VEUL,GEO,
34 . KSI,ETA,ZETA,WI,
35 . HX, HY, HZ,
36 . CJ1,CJ2,CJ3,
37 . CJ4,CJ5,CJ6,
38 . CJ7,CJ8,CJ9,
39 . SMAX, DELTAX, NGL,VOLDP)
40C-----------------------------------------------
41C M o d u l e s
42C-----------------------------------------------
43 USE message_mod
44C-----------------------------------------------
45C I m p l i c i t T y p e s
46C-----------------------------------------------
47#include "implicit_f.inc"
48C-----------------------------------------------
49C G l o b a l P a r a m e t e r s
50C-----------------------------------------------
51#include "mvsiz_p.inc"
52C-----------------------------------------------
53C C o m m o n B l o c k s
54C-----------------------------------------------
55#include "vect01_c.inc"
56#include "param_c.inc"
57C-----------------------------------------------
58C D u m m y A r g u m e n t s
59C-----------------------------------------------
60 my_real
61 . vol(*), veul(lveul,*),geo(npropg,*)
62 my_real
63 . ksi,eta,zeta,wi,
64 . hx(4,*), hy(4,*), hz(4,*),
65 . cj1(*),cj2(*),cj3(*),
66 . cj4(*),cj5(*),cj6(*),
67 . cj7(*),cj8(*),cj9(*),
68 . smax(*), deltax(*)
69 INTEGER NGL(*)
70 DOUBLE PRECISION
71 . VOLDP(*),DETDP
72C-----------------------------------------------
73C L o c a l V a r i a b l e s
74C-----------------------------------------------
75 INTEGER I,J
76C 12
77 my_real
78 . jac1(mvsiz),jac2(mvsiz),jac3(mvsiz),
79 . jac4(mvsiz),jac5(mvsiz),jac6(mvsiz),
80 . jac7(mvsiz),jac8(mvsiz),jac9(mvsiz),
81 . jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz), det(mvsiz)
82C-----------------------------------------------
83C Jacobian matrix
84 DO i=lft,llt
85 jac1(i)=cj1(i)+hx(3,i)*eta+(hx(2,i)+hx(4,i)*eta)*zeta
86 jac2(i)=cj2(i)+hy(3,i)*eta+(hy(2,i)+hy(4,i)*eta)*zeta
87 jac3(i)=cj3(i)+hz(3,i)*eta+(hz(2,i)+hz(4,i)*eta)*zeta
88C
89 jac4(i)=cj4(i)+hx(1,i)*zeta+(hx(3,i)+hx(4,i)*zeta)*ksi
90 jac5(i)=cj5(i)+hy(1,i)*zeta+(hy(3,i)+hy(4,i)*zeta)*ksi
91 jac6(i)=cj6(i)+hz(1,i)*zeta+(hz(3,i)+hz(4,i)*zeta)*ksi
92C
93 jac7(i)=cj7(i)+hx(2,i)*ksi+(hx(1,i)+hx(4,i)*ksi)*eta
94 jac8(i)=cj8(i)+hy(2,i)*ksi+(hy(1,i)+hy(4,i)*ksi)*eta
95 jac9(i)=cj9(i)+hz(2,i)*ksi+(hz(1,i)+hz(4,i)*ksi)*eta
96 ENDDO
97C
98 DO i=lft,llt
99 jac_59_68(i)=jac5(i)*jac9(i)-jac6(i)*jac8(i)
100 jac_67_49(i)=jac6(i)*jac7(i)-jac4(i)*jac9(i)
101 jac_48_57(i)=jac4(i)*jac8(i)-jac5(i)*jac7(i)
102 ENDDO
103C
104 DO i=lft,llt
105 detdp=one_over_512*(jac1(i)*jac_59_68(i)+jac2(i)*jac_67_49(i)+jac3(i)*jac_48_57(i))
106 det(i)=detdp
107 voldp(i)=wi*detdp
108 vol(i)= voldp(i)
109 ENDDO
110 DO i=lft,llt
111 IF(det(i)>0.0) cycle
112 CALL ancmsg(msgid=245,
113 . msgtype=msgerror,
114 . anmode=aninfo,
115 . i1=ngl(i))
116 END DO
117C
118C
119 DO i=lft,llt
120 deltax(i)=min(deltax(i),hundred28*det(i)*smax(i))
121 ENDDO
122 RETURN
123 END
124!||====================================================================
125!|| s8zderic3 ../starter/source/elements/solid/solide8z/s8zderi3.F
126!||--- called by ------------------------------------------------------
127!|| s8cinit3 ../starter/source/elements/thickshell/solide8c/s8cinit3.F
128!|| s8zinit3 ../starter/source/elements/solid/solide8z/s8zinit3.F
129!||--- calls -----------------------------------------------------
130!|| ancmsg ../starter/source/output/message/message.F
131!||--- uses -----------------------------------------------------
132!|| message_mod ../starter/share/message_module/message_mod.F
133!||====================================================================
134 SUBROUTINE s8zderic3(VOL,
135 . HX, HY, HZ,
136 . JAC1,JAC2,JAC3,
137 . JAC4,JAC5,JAC6,
138 . JAC7,JAC8,JAC9,SMAX, DET, NGL,
139 . XD1, XD2, XD3, XD4, XD5, XD6, XD7, XD8,
140 . YD1, YD2, YD3, YD4, YD5, YD6, YD7, YD8,
141 . ZD1, ZD2, ZD3, ZD4, ZD5, ZD6, ZD7, ZD8)
142C-----------------------------------------------
143C M o d u l e s
144C-----------------------------------------------
145 USE message_mod
146C-----------------------------------------------
147C I m p l i c i t T y p e s
148C-----------------------------------------------
149#include "implicit_f.inc"
150C-----------------------------------------------
151C G l o b a l P a r a m e t e r s
152C-----------------------------------------------
153#include "mvsiz_p.inc"
154C-----------------------------------------------
155C C o m m o n B l o c k s
156C-----------------------------------------------
157#include "vect01_c.inc"
158C-----------------------------------------------
159C D u m m y A r g u m e n t s
160C-----------------------------------------------
161 my_real
162 . vol(*),hx(4,*), hy(4,*), hz(4,*),
163 . jac1(*),jac2(*),jac3(*),
164 . jac4(*),jac5(*),jac6(*),
165 . jac7(*),jac8(*),jac9(*),smax(*), det(*)
166 double precision
167 . xd1(mvsiz), xd2(mvsiz), xd3(mvsiz), xd4(mvsiz),
168 . xd5(mvsiz), xd6(mvsiz), xd7(mvsiz), xd8(mvsiz),
169 . yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
170 . yd5(mvsiz), yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
171 . zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),
172 . zd5(mvsiz), zd6(mvsiz), zd7(mvsiz), zd8(mvsiz)
173 INTEGER NGL(*)
174C-----------------------------------------------
175C L o c a l V a r i a b l e s
176C-----------------------------------------------
177 INTEGER I, J
178
179C 12
180 my_real
181 . x17(mvsiz) , x28(mvsiz) , x35(mvsiz) , x46(mvsiz),
182 . y17(mvsiz) , y28(mvsiz) , y35(mvsiz) , y46(mvsiz),
183 . z17(mvsiz) , z28(mvsiz) , z35(mvsiz) , z46(mvsiz),
184 . jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz),
185 . jac_38_29(mvsiz), jac_19_37(mvsiz), jac_27_18(mvsiz),
186 . jac_26_35(mvsiz), jac_34_16(mvsiz), jac_15_24(mvsiz),
187 . x_17_46(mvsiz) , x_28_35(mvsiz) ,
188 . y_17_46(mvsiz) , y_28_35(mvsiz) ,
189 . z_17_46(mvsiz) , z_28_35(mvsiz)
190C-----------------------------------------------
191 DO i=lft,llt
192 x17(i)=xd7(i)-xd1(i)
193 x28(i)=xd8(i)-xd2(i)
194 x35(i)=xd5(i)-xd3(i)
195 x46(i)=xd6(i)-xd4(i)
196 y17(i)=yd7(i)-yd1(i)
197 y28(i)=yd8(i)-yd2(i)
198 y35(i)=yd5(i)-yd3(i)
199 y46(i)=yd6(i)-yd4(i)
200 z17(i)=zd7(i)-zd1(i)
201 z28(i)=zd8(i)-zd2(i)
202 z35(i)=zd5(i)-zd3(i)
203 z46(i)=zd6(i)-zd4(i)
204 END DO
205C
206 DO i=lft,llt
207 jac4(i)=x17(i)+x28(i)-x35(i)-x46(i)
208 jac5(i)=y17(i)+y28(i)-y35(i)-y46(i)
209 jac6(i)=z17(i)+z28(i)-z35(i)-z46(i)
210 x_17_46(i)=x17(i)+x46(i)
211 x_28_35(i)=x28(i)+x35(i)
212 y_17_46(i)=y17(i)+y46(i)
213 y_28_35(i)=y28(i)+y35(i)
214 z_17_46(i)=z17(i)+z46(i)
215 z_28_35(i)=z28(i)+z35(i)
216 END DO
217C
218C Jacobian matrix
219 DO i=lft,llt
220 jac7(i)=x_17_46(i)+x_28_35(i)
221 jac8(i)=y_17_46(i)+y_28_35(i)
222 jac9(i)=z_17_46(i)+z_28_35(i)
223 jac1(i)=x_17_46(i)-x_28_35(i)
224 jac2(i)=y_17_46(i)-y_28_35(i)
225 jac3(i)=z_17_46(i)-z_28_35(i)
226 END DO
227C
228 DO i=lft,llt
229 jac_59_68(i)=jac5(i)*jac9(i)-jac6(i)*jac8(i)
230 jac_67_49(i)=jac6(i)*jac7(i)-jac4(i)*jac9(i)
231 jac_38_29(i)=(-jac2(i)*jac9(i)+jac3(i)*jac8(i))
232 jac_19_37(i)=( jac1(i)*jac9(i)-jac3(i)*jac7(i))
233 jac_27_18(i)=(-jac1(i)*jac8(i)+jac2(i)*jac7(i))
234 jac_26_35(i)=( jac2(i)*jac6(i)-jac3(i)*jac5(i))
235 jac_34_16(i)=(-jac1(i)*jac6(i)+jac3(i)*jac4(i))
236 jac_15_24(i)=( jac1(i)*jac5(i)-jac2(i)*jac4(i))
237 jac_48_57(i)=jac4(i)*jac8(i)-jac5(i)*jac7(i)
238 END DO
239C
240 DO i=lft,llt
241 det(i)=one_over_64*(jac1(i)*jac_59_68(i)+jac2(i)*jac_67_49(i)+jac3(i)*jac_48_57(i))
242 vol(i)=det(i)
243 END DO
244C
245 DO i=lft,llt
246 IF(det(i)>0.0) cycle
247 CALL ancmsg(msgid=245,
248 . msgtype=msgerror,
249 . anmode=aninfo,
250 . i1=ngl(i))
251 END DO
252C
253C
254C mode 1
255C 1 1 -1 -1 -1 -1 1 1
256 DO i=lft,llt
257 hx(1,i)=(xd1(i)+xd2(i)-xd3(i)-xd4(i)-xd5(i)-xd6(i)+xd7(i)+xd8(i))
258 hy(1,i)=(yd1(i)+yd2(i)-yd3(i)-yd4(i)-yd5(i)-yd6(i)+yd7(i)+yd8(i))
259 hz(1,i)=(zd1(i)+zd2(i)-zd3(i)-zd4(i)-zd5(i)-zd6(i)+zd7(i)+zd8(i))
260 ENDDO
261C mode 2
262C 1 -1 -1 1 -1 1 1 -1
263 DO i=lft,llt
264 hx(2,i)=(xd1(i)-xd2(i)-xd3(i)+xd4(i)-xd5(i)+xd6(i)+xd7(i)-xd8(i))
265 hy(2,i)=(yd1(i)-yd2(i)-yd3(i)+yd4(i)-yd5(i)+yd6(i)+yd7(i)-yd8(i))
266 hz(2,i)=(zd1(i)-zd2(i)-zd3(i)+zd4(i)-zd5(i)+zd6(i)+zd7(i)-zd8(i))
267 ENDDO
268C mode 3
269C 1 -1 1 -1 1 -1 1 -1
270 DO i=lft,llt
271 hx(3,i)=(xd1(i)-xd2(i)+xd3(i)-xd4(i)+xd5(i)-xd6(i)+xd7(i)-xd8(i))
272 hy(3,i)=(yd1(i)-yd2(i)+yd3(i)-yd4(i)+yd5(i)-yd6(i)+yd7(i)-yd8(i))
273 hz(3,i)=(zd1(i)-zd2(i)+zd3(i)-zd4(i)+zd5(i)-zd6(i)+zd7(i)-zd8(i))
274 ENDDO
275C mode 4
276C -1 1 -1 1 1 -1 1 -1
277 DO i=lft,llt
278 hx(4,i)=(-xd1(i)+xd2(i)-xd3(i)+xd4(i)+xd5(i)-xd6(i)+xd7(i)-xd8(i))
279 hy(4,i)=(-yd1(i)+yd2(i)-yd3(i)+yd4(i)+yd5(i)-yd6(i)+yd7(i)-yd8(i))
280 hz(4,i)=(-zd1(i)+zd2(i)-zd3(i)+zd4(i)+zd5(i)-zd6(i)+zd7(i)-zd8(i))
281 ENDDO
282C----surface max mediane-- *16
283 DO i=lft,llt
284 smax(i)= jac_59_68(i)*jac_59_68(i)+jac_67_49(i)*jac_67_49(i)
285 . +jac_48_57(i)*jac_48_57(i)
286 smax(i)= max(smax(i),jac_38_29(i)*jac_38_29(i)+jac_19_37(i)*jac_19_37(i)
287 . +jac_27_18(i)*jac_27_18(i))
288 smax(i)= max(smax(i),jac_26_35(i)*jac_26_35(i)+jac_34_16(i)*jac_34_16(i)
289 . +jac_15_24(i)*jac_15_24(i))
290 ENDDO
291 DO i=lft,llt
292 IF(smax(i)<=zero)THEN
293 CALL ancmsg(msgid=245,
294 . msgtype=msgerror,
295 . anmode=aninfo,
296 . i1=ngl(i))
297 ENDIF
298 smax(i)= one/sqrt(smax(i))
299 ENDDO
300 RETURN
301C
302 END
303!||====================================================================
304!|| s8zjac_ic ../starter/source/elements/solid/solide8z/s8zderi3.F
305!||--- called by ------------------------------------------------------
306!|| s8zinit3 ../starter/source/elements/solid/solide8z/s8zinit3.F
307!||--- uses -----------------------------------------------------
308!|| message_mod ../starter/share/message_module/message_mod.F
309!||====================================================================
310 SUBROUTINE s8zjac_ic(
311 . XD1 ,XD2 ,XD3 ,XD4 ,XD5 ,XD6 ,XD7 ,XD8 ,
312 . YD1 ,YD2 ,YD3 ,YD4 ,YD5 ,YD6 ,YD7 ,YD8 ,
313 . ZD1 ,ZD2 ,ZD3 ,ZD4 ,ZD5 ,ZD6 ,ZD7 ,ZD8 ,
314 . JAC1 ,JAC2 ,JAC3 ,
315 . JAC4 ,JAC5 ,JAC6 ,
316 . JAC7 ,JAC8 ,JAC9 ,
317 . HX, HY, HZ,
318 . JAC_I)
319C-----------------------------------------------
320C M o d u l e s
321C-----------------------------------------------
322 USE message_mod
323C-----------------------------------------------
324C I m p l i c i t T y p e s
325C-----------------------------------------------
326#include "implicit_f.inc"
327C-----------------------------------------------
328C G l o b a l P a r a m e t e r s
329C-----------------------------------------------
330#include "mvsiz_p.inc"
331C-----------------------------------------------
332C C o m m o n B l o c k s
333C-----------------------------------------------
334#include "vect01_c.inc"
335C-----------------------------------------------
336C D u m m y A r g u m e n t s
337C-----------------------------------------------
338 my_real
339 . hx(4,*), hy(4,*), hz(4,*),
340 . jac1(*),jac2(*),jac3(*),
341 . jac4(*),jac5(*),jac6(*),
342 . jac7(*),jac8(*),jac9(*),jac_i(10,*)
343 double precision
344 . xd1(mvsiz), xd2(mvsiz), xd3(mvsiz), xd4(mvsiz),
345 . xd5(mvsiz), xd6(mvsiz), xd7(mvsiz), xd8(mvsiz),
346 . yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
347 . yd5(mvsiz), yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
348 . zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),
349 . zd5(mvsiz), zd6(mvsiz), zd7(mvsiz), zd8(mvsiz)
350C-----------------------------------------------
351C L o c a l V a r i a b l e s
352C-----------------------------------------------
353 INTEGER I, J
354
355 my_real
356 . DET(MVSIZ) ,DETT(MVSIZ)
357C 12
358 my_real
359 . x17(mvsiz) , x28(mvsiz) , x35(mvsiz) , x46(mvsiz),
360 . y17(mvsiz) , y28(mvsiz) , y35(mvsiz) , y46(mvsiz),
361 . z17(mvsiz) , z28(mvsiz) , z35(mvsiz) , z46(mvsiz),
362 . jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz),
363 . jac_38_29(mvsiz), jac_19_37(mvsiz), jac_27_18(mvsiz),
364 . jac_26_35(mvsiz), jac_34_16(mvsiz), jac_15_24(mvsiz),
365 . x_17_46(mvsiz) , x_28_35(mvsiz) ,
366 . y_17_46(mvsiz) , y_28_35(mvsiz) ,
367 . z_17_46(mvsiz) , z_28_35(mvsiz)
368
369C-----------------------------------------------
370 DO i=lft,llt
371 x17(i)=xd7(i)-xd1(i)
372 x28(i)=xd8(i)-xd2(i)
373 x35(i)=xd5(i)-xd3(i)
374 x46(i)=xd6(i)-xd4(i)
375 y17(i)=yd7(i)-yd1(i)
376 y28(i)=yd8(i)-yd2(i)
377 y35(i)=yd5(i)-yd3(i)
378 y46(i)=yd6(i)-yd4(i)
379 z17(i)=zd7(i)-zd1(i)
380 z28(i)=zd8(i)-zd2(i)
381 z35(i)=zd5(i)-zd3(i)
382 z46(i)=zd6(i)-zd4(i)
383 END DO
384C
385 DO i=lft,llt
386 jac4(i)=x17(i)+x28(i)-x35(i)-x46(i)
387 jac5(i)=y17(i)+y28(i)-y35(i)-y46(i)
388 jac6(i)=z17(i)+z28(i)-z35(i)-z46(i)
389 x_17_46(i)=x17(i)+x46(i)
390 x_28_35(i)=x28(i)+x35(i)
391 y_17_46(i)=y17(i)+y46(i)
392 y_28_35(i)=y28(i)+y35(i)
393 z_17_46(i)=z17(i)+z46(i)
394 z_28_35(i)=z28(i)+z35(i)
395 END DO
396C
397 DO i=lft,llt
398 jac7(i)=x_17_46(i)+x_28_35(i)
399 jac8(i)=y_17_46(i)+y_28_35(i)
400 jac9(i)=z_17_46(i)+z_28_35(i)
401 jac1(i)=x_17_46(i)-x_28_35(i)
402 jac2(i)=y_17_46(i)-y_28_35(i)
403 jac3(i)=z_17_46(i)-z_28_35(i)
404 END DO
405C
406 DO i=lft,llt
407 jac_59_68(i)=jac5(i)*jac9(i)-jac6(i)*jac8(i)
408 jac_67_49(i)=jac6(i)*jac7(i)-jac4(i)*jac9(i)
409 jac_38_29(i)=(-jac2(i)*jac9(i)+jac3(i)*jac8(i))
410 jac_19_37(i)=( jac1(i)*jac9(i)-jac3(i)*jac7(i))
411 jac_27_18(i)=(-jac1(i)*jac8(i)+jac2(i)*jac7(i))
412 jac_26_35(i)=( jac2(i)*jac6(i)-jac3(i)*jac5(i))
413 jac_34_16(i)=(-jac1(i)*jac6(i)+jac3(i)*jac4(i))
414 jac_15_24(i)=( jac1(i)*jac5(i)-jac2(i)*jac4(i))
415 jac_48_57(i)=jac4(i)*jac8(i)-jac5(i)*jac7(i)
416 END DO
417C
418 DO i=lft,llt
419 det(i)=one_over_64*(jac1(i)*jac_59_68(i)+jac2(i)*jac_67_49(i)+jac3(i)*jac_48_57(i))
420 dett(i)=one_over_64/det(i)
421 END DO
422C
423 DO i=lft,llt
424 jac_i(1,i)=dett(i)*jac_59_68(i)
425 jac_i(4,i)=dett(i)*jac_67_49(i)
426 jac_i(7,i)=dett(i)*jac_48_57(i)
427 jac_i(2,i)=dett(i)*jac_38_29(i)
428 jac_i(5,i)=dett(i)*jac_19_37(i)
429 jac_i(8,i)=dett(i)*jac_27_18(i)
430 jac_i(3,i)=dett(i)*jac_26_35(i)
431 jac_i(6,i)=dett(i)*jac_34_16(i)
432 jac_i(9,i)=dett(i)*jac_15_24(i)
433 jac_i(10,i)=det(i)
434 ENDDO
435C mode 1
436C 1 1 -1 -1 -1 -1 1 1
437 DO i=lft,llt
438 hx(1,i)=(xd1(i)+xd2(i)-xd3(i)-xd4(i)-xd5(i)-xd6(i)+xd7(i)+xd8(i))
439 hy(1,i)=(yd1(i)+yd2(i)-yd3(i)-yd4(i)-yd5(i)-yd6(i)+yd7(i)+yd8(i))
440 hz(1,i)=(zd1(i)+zd2(i)-zd3(i)-zd4(i)-zd5(i)-zd6(i)+zd7(i)+zd8(i))
441 ENDDO
442C mode 2
443C 1 -1 -1 1 -1 1 1 -1
444 DO i=lft,llt
445 hx(2,i)=(xd1(i)-xd2(i)-xd3(i)+xd4(i)-xd5(i)+xd6(i)+xd7(i)-xd8(i))
446 hy(2,i)=(yd1(i)-yd2(i)-yd3(i)+yd4(i)-yd5(i)+yd6(i)+yd7(i)-yd8(i))
447 hz(2,i)=(zd1(i)-zd2(i)-zd3(i)+zd4(i)-zd5(i)+zd6(i)+zd7(i)-zd8(i))
448 ENDDO
449C mode 3
450C 1 -1 1 -1 1 -1 1 -1
451 DO i=lft,llt
452 hx(3,i)=(xd1(i)-xd2(i)+xd3(i)-xd4(i)+xd5(i)-xd6(i)+xd7(i)-xd8(i))
453 hy(3,i)=(yd1(i)-yd2(i)+yd3(i)-yd4(i)+yd5(i)-yd6(i)+yd7(i)-yd8(i))
454 hz(3,i)=(zd1(i)-zd2(i)+zd3(i)-zd4(i)+zd5(i)-zd6(i)+zd7(i)-zd8(i))
455 ENDDO
456C mode 4
457C -1 1 -1 1 1 -1 1 -1
458 DO i=lft,llt
459 hx(4,i)=(-xd1(i)+xd2(i)-xd3(i)+xd4(i)+xd5(i)-xd6(i)+xd7(i)-xd8(i))
460 hy(4,i)=(-yd1(i)+yd2(i)-yd3(i)+yd4(i)+yd5(i)-yd6(i)+yd7(i)-yd8(i))
461 hz(4,i)=(-zd1(i)+zd2(i)-zd3(i)+zd4(i)+zd5(i)-zd6(i)+zd7(i)-zd8(i))
462 ENDDO
463C
464 RETURN
465C
466 END
467!||====================================================================
468!|| s8zpij_ic ../starter/source/elements/solid/solide8z/s8zderi3.F
469!||--- called by ------------------------------------------------------
470!|| s8zinit3 ../starter/source/elements/solid/solide8z/s8zinit3.F
471!||--- uses -----------------------------------------------------
472!|| message_mod ../starter/share/message_module/message_mod.f
473!||====================================================================
474 SUBROUTINE s8zpij_ic(
475 . XD1 ,XD2 ,XD3 ,XD4 ,XD5 ,XD6 ,XD7 ,XD8 ,
476 . YD1 ,YD2 ,YD3 ,YD4 ,YD5 ,YD6 ,YD7 ,YD8 ,
477 . ZD1 ,ZD2 ,ZD3 ,ZD4 ,ZD5 ,ZD6 ,ZD7 ,ZD8 ,
478 . JAC1 ,JAC2 ,JAC3 ,
479 . JAC4 ,JAC5 ,JAC6 ,
480 . JAC7 ,JAC8 ,JAC9 ,
481 . HX, HY, HZ,
482 . PXC1, PXC2, PXC3, PXC4,
483 . PYC1, PYC2, PYC3, PYC4,
484 . PZC1, PZC2, PZC3, PZC4)
485C-----------------------------------------------
486C M o d u l e s
487C-----------------------------------------------
488 USE message_mod
489C-----------------------------------------------
490C I m p l i c i t T y p e s
491C-----------------------------------------------
492#include "implicit_f.inc"
493C-----------------------------------------------
494C G l o b a l P a r a m e t e r s
495C-----------------------------------------------
496#include "mvsiz_p.inc"
497C-----------------------------------------------
498C C o m m o n B l o c k s
499C-----------------------------------------------
500#include "vect01_c.inc"
501C-----------------------------------------------
502C D u m m y A r g u m e n t s
503C-----------------------------------------------
504 my_real
505 . hx(4,*), hy(4,*), hz(4,*),
506 . pxc1(*), pxc2(*), pxc3(*), pxc4(*),
507 . pyc1(*), pyc2(*), pyc3(*), pyc4(*),
508 . pzc1(*), pzc2(*), pzc3(*), pzc4(*),
509 . jac1(*),jac2(*),jac3(*),
510 . jac4(*),jac5(*),jac6(*),
511 . jac7(*),jac8(*),jac9(*)
512 double precision
513 . xd1(mvsiz), xd2(mvsiz), xd3(mvsiz), xd4(mvsiz),
514 . xd5(mvsiz), xd6(mvsiz), xd7(mvsiz), xd8(mvsiz),
515 . yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
516 . yd5(mvsiz), yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
517 . zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),
518 . zd5(mvsiz), zd6(mvsiz), zd7(mvsiz), zd8(mvsiz)
519C-----------------------------------------------
520C L o c a l V a r i a b l e s
521C-----------------------------------------------
522 INTEGER I, J
523
524 my_real
525 . DET(MVSIZ) ,DETT(MVSIZ)
526C 12
527 my_real
528 . X17(MVSIZ) , X28(MVSIZ) , X35(MVSIZ) , X46(MVSIZ),
529 . Y17(MVSIZ) , Y28(MVSIZ) , Y35(MVSIZ) , Y46(MVSIZ),
530 . Z17(MVSIZ) , Z28(MVSIZ) , Z35(MVSIZ) , Z46(MVSIZ),
531 . JAC_59_68(MVSIZ), JAC_67_49(MVSIZ), JAC_48_57(MVSIZ),
532 . JAC_38_29(MVSIZ), JAC_19_37(MVSIZ), JAC_27_18(MVSIZ),
533 . JAC_26_35(MVSIZ), JAC_34_16(MVSIZ), JAC_15_24(MVSIZ),
534 . JACI1, JACI2, JACI3,JACI4, JACI5, JACI6,JACI7, JACI8, JACI9,
535 . jaci12, jaci45, jaci78,jaci12p, jaci45p, jaci78p,
536 . x_17_46(mvsiz) , x_28_35(mvsiz) ,
537 . y_17_46(mvsiz) , y_28_35(mvsiz) ,
538 . z_17_46(mvsiz) , z_28_35(mvsiz)
539
540C-----------------------------------------------
541 DO i=lft,llt
542 x17(i)=xd7(i)-xd1(i)
543 x28(i)=xd8(i)-xd2(i)
544 x35(i)=xd5(i)-xd3(i)
545 x46(i)=xd6(i)-xd4(i)
546 y17(i)=yd7(i)-yd1(i)
547 y28(i)=yd8(i)-yd2(i)
548 y35(i)=yd5(i)-yd3(i)
549 y46(i)=yd6(i)-yd4(i)
550 z17(i)=zd7(i)-zd1(i)
551 z28(i)=zd8(i)-zd2(i)
552 z35(i)=zd5(i)-zd3(i)
553 z46(i)=zd6(i)-zd4(i)
554 END DO
555C
556 DO i=lft,llt
557 jac4(i)=x17(i)+x28(i)-x35(i)-x46(i)
558 jac5(i)=y17(i)+y28(i)-y35(i)-y46(i)
559 jac6(i)=z17(i)+z28(i)-z35(i)-z46(i)
560 x_17_46(i)=x17(i)+x46(i)
561 x_28_35(i)=x28(i)+x35(i)
562 y_17_46(i)=y17(i)+y46(i)
563 y_28_35(i)=y28(i)+y35(i)
564 z_17_46(i)=z17(i)+z46(i)
565 z_28_35(i)=z28(i)+z35(i)
566 END DO
567C
568 DO i=lft,llt
569 jac7(i)=x_17_46(i)+x_28_35(i)
570 jac8(i)=y_17_46(i)+y_28_35(i)
571 jac9(i)=z_17_46(i)+z_28_35(i)
572 jac1(i)=x_17_46(i)-x_28_35(i)
573 jac2(i)=y_17_46(i)-y_28_35(i)
574 jac3(i)=z_17_46(i)-z_28_35(i)
575 END DO
576C
577 DO i=lft,llt
578 jac_59_68(i)=jac5(i)*jac9(i)-jac6(i)*jac8(i)
579 jac_67_49(i)=jac6(i)*jac7(i)-jac4(i)*jac9(i)
580 jac_38_29(i)=(-jac2(i)*jac9(i)+jac3(i)*jac8(i))
581 jac_19_37(i)=( jac1(i)*jac9(i)-jac3(i)*jac7(i))
582 jac_27_18(i)=(-jac1(i)*jac8(i)+jac2(i)*jac7(i))
583 jac_26_35(i)=( jac2(i)*jac6(i)-jac3(i)*jac5(i))
584 jac_34_16(i)=(-jac1(i)*jac6(i)+jac3(i)*jac4(i))
585 jac_15_24(i)=( jac1(i)*jac5(i)-jac2(i)*jac4(i))
586 jac_48_57(i)=jac4(i)*jac8(i)-jac5(i)*jac7(i)
587 END DO
588C
589 DO i=lft,llt
590 det(i)=one_over_64*(jac1(i)*jac_59_68(i)+jac2(i)*jac_67_49(i)+jac3(i)*jac_48_57(i))
591 dett(i)=one_over_64/det(i)
592 END DO
593C
594 DO i=lft,llt
595 jaci1=dett(i)*jac_59_68(i)
596 jaci4=dett(i)*jac_67_49(i)
597 jaci7=dett(i)*jac_48_57(i)
598 jaci2=dett(i)*jac_38_29(i)
599 jaci5=dett(i)*jac_19_37(i)
600 jaci8=dett(i)*jac_27_18(i)
601 jaci3=dett(i)*jac_26_35(i)
602 jaci6=dett(i)*jac_34_16(i)
603 jaci9=dett(i)*jac_15_24(i)
604 jaci12=jaci1-jaci2
605 jaci45=jaci4-jaci5
606 jaci78=jaci7-jaci8
607 pxc2(i)= jaci12-jaci3
608 pyc2(i)= jaci45-jaci6
609 pzc2(i)= jaci78-jaci9
610 pxc4(i)=-jaci12-jaci3
611 pyc4(i)=-jaci45-jaci6
612 pzc4(i)=-jaci78-jaci9
613 jaci12p=jaci1+jaci2
614 jaci45p=jaci4+jaci5
615 jaci78p=jaci7+jaci8
616 pxc1(i)=-jaci12p-jaci3
617 pyc1(i)=-jaci45p-jaci6
618 pzc1(i)=-jaci78p-jaci9
619 pxc3(i)= jaci12p-jaci3
620 pyc3(i)= jaci45p-jaci6
621 pzc3(i)= jaci78p-jaci9
622 ENDDO
623C mode 1
624C 1 1 -1 -1 -1 -1 1 1
625 DO i=lft,llt
626 hx(1,i)=(xd1(i)+xd2(i)-xd3(i)-xd4(i)-xd5(i)-xd6(i)+xd7(i)+xd8(i))
627 hy(1,i)=(yd1(i)+yd2(i)-yd3(i)-yd4(i)-yd5(i)-yd6(i)+yd7(i)+yd8(i))
628 hz(1,i)=(zd1(i)+zd2(i)-zd3(i)-zd4(i)-zd5(i)-zd6(i)+zd7(i)+zd8(i))
629 ENDDO
630C mode 2
631C 1 -1 -1 1 -1 1 1 -1
632 DO i=lft,llt
633 hx(2,i)=(xd1(i)-xd2(i)-xd3(i)+xd4(i)-xd5(i)+xd6(i)+xd7(i)-xd8(i))
634 hy(2,i)=(yd1(i)-yd2(i)-yd3(i)+yd4(i)-yd5(i)+yd6(i)+yd7(i)-yd8(i))
635 hz(2,i)=(zd1(i)-zd2(i)-zd3(i)+zd4(i)-zd5(i)+zd6(i)+zd7(i)-zd8(i))
636 ENDDO
637C mode 3
638C 1 -1 1 -1 1 -1 1 -1
639 DO i=lft,llt
640 hx(3,i)=(xd1(i)-xd2(i)+xd3(i)-xd4(i)+xd5(i)-xd6(i)+xd7(i)-xd8(i))
641 hy(3,i)=(yd1(i)-yd2(i)+yd3(i)-yd4(i)+yd5(i)-yd6(i)+yd7(i)-yd8(i))
642 hz(3,i)=(zd1(i)-zd2(i)+zd3(i)-zd4(i)+zd5(i)-zd6(i)+zd7(i)-zd8(i))
643 ENDDO
644C mode 4
645C -1 1 -1 1 1 -1 1 -1
646 DO i=lft,llt
647 hx(4,i)=(-xd1(i)+xd2(i)-xd3(i)+xd4(i)+xd5(i)-xd6(i)+xd7(i)-xd8(i))
648 hy(4,i)=(-yd1(i)+yd2(i)-yd3(i)+yd4(i)+yd5(i)-yd6(i)+yd7(i)-yd8(i))
649 hz(4,i)=(-zd1(i)+zd2(i)-zd3(i)+zd4(i)+zd5(i)-zd6(i)+zd7(i)-zd8(i))
650 ENDDO
651C
652 RETURN
653C
654 END
655!||====================================================================
656!|| s8zjac_i3 ../starter/source/elements/solid/solide8z/s8zderi3.F
657!||--- called by ------------------------------------------------------
658!|| s8zinit3 ../starter/source/elements/solid/solide8z/s8zinit3.F
659!||--- uses -----------------------------------------------------
660!|| message_mod ../starter/share/message_module/message_mod.F
661!||====================================================================
662 SUBROUTINE s8zjac_i3(
663 . KSI,ETA,ZETA,WI,
664 . HX, HY, HZ,
665 . CJ1,CJ2,CJ3,
666 . CJ4,CJ5,CJ6,
667 . CJ7,CJ8,CJ9,
668 . JAC_I,L_PIJ,PIJ,NEL)
669C-----------------------------------------------
670 USE message_mod
671C-----------------------------------------------
672C I m p l i c i t T y p e s
673C-----------------------------------------------
674#include "implicit_f.inc"
675C-----------------------------------------------
676C G l o b a l P a r a m e t e r s
677C-----------------------------------------------
678#include "mvsiz_p.inc"
679C-----------------------------------------------
680C C o m m o n B l o c k s
681C-----------------------------------------------
682#include "vect01_c.inc"
683C-----------------------------------------------
684C D u m m y A r g u m e n t s
685C-----------------------------------------------
686 INTEGER L_PIJ,NEL
687 my_real
688 . ksi,eta,zeta,wi,
689 . hx(4,*), hy(4,*), hz(4,*),
690 . cj1(*),cj2(*),cj3(*),
691 . cj4(*),cj5(*),cj6(*),
692 . cj7(*),cj8(*),cj9(*),
693 . jac_i(10,*),pij(nel,*)
694C-----------------------------------------------
695C L o c a l V a r i a b l e s
696C-----------------------------------------------
697 INTEGER I, J
698 my_real
699 . det(mvsiz) , dett(mvsiz) ,
700 . jac1(mvsiz),jac2(mvsiz),jac3(mvsiz),
701 . jac4(mvsiz),jac5(mvsiz),jac6(mvsiz),
702 . jac7(mvsiz),jac8(mvsiz),jac9(mvsiz),
703 . jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz),
704 . jac_38_29(mvsiz), jac_19_37(mvsiz), jac_27_18(mvsiz),
705 . jac_26_35(mvsiz), jac_34_16(mvsiz), jac_15_24(mvsiz)
706 my_real
707 . jaci1(mvsiz), jaci2(mvsiz), jaci3(mvsiz),
708 . jaci4(mvsiz), jaci5(mvsiz), jaci6(mvsiz),
709 . jaci7(mvsiz), jaci8(mvsiz), jaci9(mvsiz),
710 . a1pr1(mvsiz),a1pr3(mvsiz),a1pr5(mvsiz),a1pr7(mvsiz),
711 . a4pr1(mvsiz),a4pr3(mvsiz),a4pr5(mvsiz),a4pr7(mvsiz),
712 . a7pr1(mvsiz),a7pr3(mvsiz),a7pr5(mvsiz),a7pr7(mvsiz),
713 . a2ps1(mvsiz),a2ps2(mvsiz),a2ps5(mvsiz),a2ps6(mvsiz),
714 . a5ps1(mvsiz),a5ps2(mvsiz),a5ps5(mvsiz),a5ps6(mvsiz),
715 . a8ps1(mvsiz),a8ps2(mvsiz),a8ps5(mvsiz),a8ps6(mvsiz),
716 . a3pt1(mvsiz),a3pt2(mvsiz),a3pt3(mvsiz),a3pt4(mvsiz),
717 . a6pt1(mvsiz),a6pt2(mvsiz),a6pt3(mvsiz),a6pt4(mvsiz),
718 . a9pt1(mvsiz),a9pt2(mvsiz),a9pt3(mvsiz),a9pt4(mvsiz),
719 . pr(8),ps(8),pt(8),rp, sp, tp, rm, sm, tm,
720 . px1, px2, px3, px4,px5, px6, px7, px8,
721 . py1, py2, py3, py4,py5, py6, py7, py8,
722 . pz1, pz2, pz3, pz4,pz5, pz6, pz7, pz8
723C-----------------------------------------------
724 DO i=lft,llt
725 jac1(i)=cj1(i)+hx(3,i)*eta+(hx(2,i)+hx(4,i)*eta)*zeta
726 jac2(i)=cj2(i)+hy(3,i)*eta+(hy(2,i)+hy(4,i)*eta)*zeta
727 jac3(i)=cj3(i)+hz(3,i)*eta+(hz(2,i)+hz(4,i)*eta)*zeta
728C
729 jac4(i)=cj4(i)+hx(1,i)*zeta+(hx(3,i)+hx(4,i)*zeta)*ksi
730 jac5(i)=cj5(i)+hy(1,i)*zeta+(hy(3,i)+hy(4,i)*zeta)*ksi
731 jac6(i)=cj6(i)+hz(1,i)*zeta+(hz(3,i)+hz(4,i)*zeta)*ksi
732C
733 jac7(i)=cj7(i)+hx(2,i)*ksi+(hx(1,i)+hx(4,i)*ksi)*eta
734 jac8(i)=cj8(i)+hy(2,i)*ksi+(hy(1,i)+hy(4,i)*ksi)*eta
735 jac9(i)=cj9(i)+hz(2,i)*ksi+(hz(1,i)+hz(4,i)*ksi)*eta
736 ENDDO
737C
738C JACOBIAN
739C
740 DO i=lft,llt
741 jac_59_68(i)=jac5(i)*jac9(i)-jac6(i)*jac8(i)
742 jac_67_49(i)=jac6(i)*jac7(i)-jac4(i)*jac9(i)
743 jac_38_29(i)=(-jac2(i)*jac9(i)+jac3(i)*jac8(i))
744 jac_19_37(i)=( jac1(i)*jac9(i)-jac3(i)*jac7(i))
745 jac_27_18(i)=(-jac1(i)*jac8(i)+jac2(i)*jac7(i))
746 jac_26_35(i)=( jac2(i)*jac6(i)-jac3(i)*jac5(i))
747 jac_34_16(i)=(-jac1(i)*jac6(i)+jac3(i)*jac4(i))
748 jac_15_24(i)=( jac1(i)*jac5(i)-jac2(i)*jac4(i))
749 jac_48_57(i)=jac4(i)*jac8(i)-jac5(i)*jac7(i)
750 ENDDO
751C
752 DO i=lft,llt
753 det(i)=one_over_512*(jac1(i)*jac_59_68(i)+jac2(i)*jac_67_49(i)+jac3(i)*jac_48_57(i))
754 dett(i)=one_over_512/det(i)
755 ENDDO
756C
757 IF (l_pij ==0) THEN
758 DO i=lft,llt
759 jac_i(1,i)=dett(i)*jac_59_68(i)
760 jac_i(4,i)=dett(i)*jac_67_49(i)
761 jac_i(7,i)=dett(i)*jac_48_57(i)
762 jac_i(2,i)=dett(i)*jac_38_29(i)
763 jac_i(5,i)=dett(i)*jac_19_37(i)
764 jac_i(8,i)=dett(i)*jac_27_18(i)
765 jac_i(3,i)=dett(i)*jac_26_35(i)
766 jac_i(6,i)=dett(i)*jac_34_16(i)
767 jac_i(9,i)=dett(i)*jac_15_24(i)
768 jac_i(10,i)=det(i)
769 ENDDO
770 ELSE
771C--------calculate Pij----
772 DO i=lft,llt
773 jaci1(i)=dett(i)*jac_59_68(i)
774 jaci4(i)=dett(i)*jac_67_49(i)
775 jaci7(i)=dett(i)*jac_48_57(i)
776 jaci2(i)=dett(i)*jac_38_29(i)
777 jaci5(i)=dett(i)*jac_19_37(i)
778 jaci8(i)=dett(i)*jac_27_18(i)
779 jaci3(i)=dett(i)*jac_26_35(i)
780 jaci6(i)=dett(i)*jac_34_16(i)
781 jaci9(i)=dett(i)*jac_15_24(i)
782 ENDDO
783C------NI,ksi-------
784 rp = one + ksi
785 sp = one + eta
786 tp = one + zeta
787 rm = one - ksi
788 sm = one - eta
789 tm = one - zeta
790 pr(1)=-sm*tm
791 pr(2)=-pr(1)
792 pr(3)= sp*tm
793 pr(4)=-pr(3)
794 pr(5)=-sm*tp
795 pr(6)=-pr(5)
796 pr(7)= sp*tp
797 pr(8)=-pr(7)
798 ps(1)=-rm*tm
799 ps(2)=-rp*tm
800 ps(3)=-ps(2)
801 ps(4)=-ps(1)
802 ps(5)=-rm*tp
803 ps(6)=-rp*tp
804 ps(7)=-ps(6)
805 ps(8)=-ps(5)
806 pt(1)=-rm*sm
807 pt(2)=-rp*sm
808 pt(3)=-rp*sp
809 pt(4)=-rm*sp
810 pt(5)=-pt(1)
811 pt(6)=-pt(2)
812 pt(7)=-pt(3)
813 pt(8)=-pt(4)
814cc
815 DO i=lft,llt
816 a1pr1(i)=jaci1(i)*pr(1)
817 a1pr3(i)=jaci1(i)*pr(3)
818 a1pr5(i)=jaci1(i)*pr(5)
819 a1pr7(i)=jaci1(i)*pr(7)
820 a4pr1(i)=jaci4(i)*pr(1)
821 a4pr3(i)=jaci4(i)*pr(3)
822 a4pr5(i)=jaci4(i)*pr(5)
823 a4pr7(i)=jaci4(i)*pr(7)
824 a7pr1(i)=jaci7(i)*pr(1)
825 a7pr3(i)=jaci7(i)*pr(3)
826 a7pr5(i)=jaci7(i)*pr(5)
827 a7pr7(i)=jaci7(i)*pr(7)
828 ENDDO
829 DO i=lft,llt
830 a2ps1(i)=jaci2(i)*ps(1)
831 a2ps2(i)=jaci2(i)*ps(2)
832 a2ps5(i)=jaci2(i)*ps(5)
833 a2ps6(i)=jaci2(i)*ps(6)
834 a5ps1(i)=jaci5(i)*ps(1)
835 a5ps2(i)=jaci5(i)*ps(2)
836 a5ps5(i)=jaci5(i)*ps(5)
837 a5ps6(i)=jaci5(i)*ps(6)
838 a8ps1(i)=jaci8(i)*ps(1)
839 a8ps2(i)=jaci8(i)*ps(2)
840 a8ps5(i)=jaci8(i)*ps(5)
841 a8ps6(i)=jaci8(i)*ps(6)
842 ENDDO
843C
844 DO i=lft,llt
845 a3pt1(i)=jaci3(i)*pt(1)
846 a3pt2(i)=jaci3(i)*pt(2)
847 a3pt3(i)=jaci3(i)*pt(3)
848 a3pt4(i)=jaci3(i)*pt(4)
849 a6pt1(i)=jaci6(i)*pt(1)
850 a6pt2(i)=jaci6(i)*pt(2)
851 a6pt3(i)=jaci6(i)*pt(3)
852 a6pt4(i)=jaci6(i)*pt(4)
853 a9pt1(i)=jaci9(i)*pt(1)
854 a9pt2(i)=jaci9(i)*pt(2)
855 a9pt3(i)=jaci9(i)*pt(3)
856 a9pt4(i)=jaci9(i)*pt(4)
857 ENDDO
858C
859 DO i=lft,llt
860 px1= a1pr1(i)+a2ps1(i)+a3pt1(i)
861 px2=-a1pr1(i)+a2ps2(i)+a3pt2(i)
862 px3= a1pr3(i)-a2ps2(i)+a3pt3(i)
863 px4=-a1pr3(i)-a2ps1(i)+a3pt4(i)
864 px5= a1pr5(i)+a2ps5(i)-a3pt1(i)
865 px6=-a1pr5(i)+a2ps6(i)-a3pt2(i)
866 px7= a1pr7(i)-a2ps6(i)-a3pt3(i)
867 px8=-a1pr7(i)-a2ps5(i)-a3pt4(i)
868C
869 py1= a4pr1(i)+a5ps1(i)+a6pt1(i)
870 py2=-a4pr1(i)+a5ps2(i)+a6pt2(i)
871 py3= a4pr3(i)-a5ps2(i)+a6pt3(i)
872 py4=-a4pr3(i)-a5ps1(i)+a6pt4(i)
873 py5= a4pr5(i)+a5ps5(i)-a6pt1(i)
874 py6=-a4pr5(i)+a5ps6(i)-a6pt2(i)
875 py7= a4pr7(i)-a5ps6(i)-a6pt3(i)
876 py8=-a4pr7(i)-a5ps5(i)-a6pt4(i)
877C
878 pz1= a7pr1(i)+a8ps1(i)+a9pt1(i)
879 pz2=-a7pr1(i)+a8ps2(i)+a9pt2(i)
880 pz3= a7pr3(i)-a8ps2(i)+a9pt3(i)
881 pz4=-a7pr3(i)-a8ps1(i)+a9pt4(i)
882 pz5= a7pr5(i)+a8ps5(i)-a9pt1(i)
883 pz6=-a7pr5(i)+a8ps6(i)-a9pt2(i)
884 pz7= a7pr7(i)-a8ps6(i)-a9pt3(i)
885 pz8=-a7pr7(i)-a8ps5(i)-a9pt4(i)
886C
887 pij(i,1) = px1
888 pij(i,2) = py1
889 pij(i,3) = pz1
890 pij(i,4) = px2
891 pij(i,5) = py2
892 pij(i,6) = pz2
893 pij(i,7) = px3
894 pij(i,8) = py3
895 pij(i,9) = pz3
896 pij(i,10) = px4
897 pij(i,11) = py4
898 pij(i,12) = pz4
899 pij(i,13) = px5
900 pij(i,14) = py5
901 pij(i,15) = pz5
902 pij(i,16) = px6
903 pij(i,17) = py6
904 pij(i,18) = pz6
905 pij(i,19) = px7
906 pij(i,20) = py7
907 pij(i,21) = pz7
908 pij(i,22) = px8
909 pij(i,23) = py8
910 pij(i,24) = pz8
911 ENDDO
912 END IF !(L_JAC >0) THEN
913C
914 RETURN
915 END
916!||====================================================================
917!|| s8e_pij ../starter/source/elements/solid/solide8z/s8zderi3.F
918!||--- called by ------------------------------------------------------
919!|| s8zinit3 ../starter/source/elements/solid/solide8z/s8zinit3.F
920!||--- calls -----------------------------------------------------
921!|| s8ederi_bij ../starter/source/elements/solid/solide8z/s8zderi3.F
922!|| s8ederil2g ../starter/source/elements/solid/solide8z/s8zderi3.F
923!|| s8ederish2 ../starter/source/elements/solid/solide8z/s8zderi3.F
924!|| s8egetpij ../starter/source/elements/solid/solide8z/s8zderi3.F
925!|| s8eselecsh ../starter/source/elements/solid/solide8z/s8zderi3.F
926!|| setmatrans ../starter/source/elements/solid/solide8z/s8zderi3.F
927!||--- uses -----------------------------------------------------
928!||====================================================================
929 SUBROUTINE s8e_pij(NPTR,NPTS,NPTT,NNPT,NEL,
930 . PXC1, PXC2, PXC3, PXC4,
931 . PYC1, PYC2, PYC3, PYC4,
932 . PZC1, PZC2, PZC3, PZC4,
933 . E1X ,E1Y ,E1Z ,E2X ,E2Y ,E2Z ,E3X ,E3Y ,E3Z ,
934 . NU ,ELBUF_STR)
935C-----------------------------------------------
936C M o d u l e s
937C-----------------------------------------------
938 USE elbufdef_mod
939C-----------------------------------------------
940C I m p l i c i t T y p e s
941C-----------------------------------------------
942#include "implicit_f.inc"
943C-----------------------------------------------
944C G l o b a l P a r a m e t e r s
945C-----------------------------------------------
946#include "mvsiz_p.inc"
947C-----------------------------------------------
948C C o m m o n B l o c k s
949C-----------------------------------------------
950#include "vect01_c.inc"
951C-----------------------------------------------
952C D u m m y A r g u m e n t s
953C-----------------------------------------------
954 INTEGER NPTR,NPTS,NPTT,NNPT,NEL
955 my_real
956 . pxc1(*), pxc2(*), pxc3(*), pxc4(*),
957 . pyc1(*), pyc2(*), pyc3(*), pyc4(*),
958 . pzc1(*), pzc2(*), pzc3(*), pzc4(*),
959 . e1x(mvsiz),e1y(mvsiz),e1z(mvsiz),e2x(mvsiz),
960 . e2y(mvsiz),e2z(mvsiz),e3x(mvsiz),e3y(mvsiz),e3z(mvsiz),nu(*)
961 TYPE(ELBUF_STRUCT_), TARGET :: ELBUF_STR
962C-----------------------------------------------
963C L o c a l V a r i a b l e s
964C-----------------------------------------------
965 INTEGER IP,IR,IS,IT,IL,i
966C-----------------------------------------------
967 TYPE(L_BUFEL_) ,POINTER :: LBUF
968 TYPE(G_BUFEL_) ,POINTER :: GBUF
969c-----------------------------------------------------
970 my_real
971 . bxy1(mvsiz),bxy2(mvsiz),bxy3(mvsiz),bxy4(mvsiz),
972 . bxy5(mvsiz),bxy6(mvsiz),bxy7(mvsiz),bxy8(mvsiz),
973 . byx1(mvsiz),byx2(mvsiz),byx3(mvsiz),byx4(mvsiz),
974 . byx5(mvsiz),byx6(mvsiz),byx7(mvsiz),byx8(mvsiz),
975 . bxz1(mvsiz),bxz2(mvsiz),bxz3(mvsiz),bxz4(mvsiz),
976 . bxz5(mvsiz),bxz6(mvsiz),bxz7(mvsiz),bxz8(mvsiz),
977 . bzx1(mvsiz),bzx2(mvsiz),bzx3(mvsiz),bzx4(mvsiz),
978 . bzx5(mvsiz),bzx6(mvsiz),bzx7(mvsiz),bzx8(mvsiz),
979 . byz1(mvsiz),byz2(mvsiz),byz3(mvsiz),byz4(mvsiz),
980 . byz5(mvsiz),byz6(mvsiz),byz7(mvsiz),byz8(mvsiz),
981 . bzy1(mvsiz),bzy2(mvsiz),bzy3(mvsiz),bzy4(mvsiz),
982 . bzy5(mvsiz),bzy6(mvsiz),bzy7(mvsiz),bzy8(mvsiz)
983 my_real
984 . pxy1(mvsiz),pxy2(mvsiz),pxy3(mvsiz),pxy4(mvsiz),
985 . pxy5(mvsiz),pxy6(mvsiz),pxy7(mvsiz),pxy8(mvsiz),
986 . pyx1(mvsiz),pyx2(mvsiz),pyx3(mvsiz),pyx4(mvsiz),
987 . pyx5(mvsiz),pyx6(mvsiz),pyx7(mvsiz),pyx8(mvsiz),
988 . pxz1(mvsiz),pxz2(mvsiz),pxz3(mvsiz),pxz4(mvsiz),
989 . pxz5(mvsiz),pxz6(mvsiz),pxz7(mvsiz),pxz8(mvsiz),
990 . pzx1(mvsiz),pzx2(mvsiz),pzx3(mvsiz),pzx4(mvsiz),
991 . pzx5(mvsiz),pzx6(mvsiz),pzx7(mvsiz),pzx8(mvsiz),
992 . pyz1(mvsiz),pyz2(mvsiz),pyz3(mvsiz),pyz4(mvsiz),
993 . pyz5(mvsiz),pyz6(mvsiz),pyz7(mvsiz),pyz8(mvsiz),
994 . pzy1(mvsiz),pzy2(mvsiz),pzy3(mvsiz),pzy4(mvsiz),
995 . pzy5(mvsiz),pzy6(mvsiz),pzy7(mvsiz),pzy8(mvsiz),
996 . px1(mvsiz,nnpt),px2(mvsiz,nnpt),px3(mvsiz,nnpt),px4(mvsiz,nnpt),
997 . px5(mvsiz,nnpt),px6(mvsiz,nnpt),px7(mvsiz,nnpt),px8(mvsiz,nnpt),
998 . py1(mvsiz,nnpt),py2(mvsiz,nnpt),py3(mvsiz,nnpt),py4(mvsiz,nnpt),
999 . py5(mvsiz,nnpt),py6(mvsiz,nnpt),py7(mvsiz,nnpt),py8(mvsiz,nnpt),
1000 . pz1(mvsiz,nnpt),pz2(mvsiz,nnpt),pz3(mvsiz,nnpt),pz4(mvsiz,nnpt),
1001 . pz5(mvsiz,nnpt),pz6(mvsiz,nnpt),pz7(mvsiz,nnpt),pz8(mvsiz,nnpt),
1002 . p0xy1(mvsiz,2),p0xy2(mvsiz,2),p0xy3(mvsiz,2),p0xy4(mvsiz,2),
1003 . p0xy5(mvsiz,2),p0xy6(mvsiz,2),p0xy7(mvsiz,2),p0xy8(mvsiz,2),
1004 . p0yx1(mvsiz,2),p0yx2(mvsiz,2),p0yx3(mvsiz,2),p0yx4(mvsiz,2),
1005 . p0yx5(mvsiz,2),p0yx6(mvsiz,2),p0yx7(mvsiz,2),p0yx8(mvsiz,2),
1006 . p0xz1(mvsiz,2),p0xz2(mvsiz,2),p0xz3(mvsiz,2),p0xz4(mvsiz,2),
1007 . p0xz5(mvsiz,2),p0xz6(mvsiz,2),p0xz7(mvsiz,2),p0xz8(mvsiz,2),
1008 . p0zx1(mvsiz,2),p0zx2(mvsiz,2),p0zx3(mvsiz,2),p0zx4(mvsiz,2),
1009 . p0zx5(mvsiz,2),p0zx6(mvsiz,2),p0zx7(mvsiz,2),p0zx8(mvsiz,2),
1010 . p0yz1(mvsiz,2),p0yz2(mvsiz,2),p0yz3(mvsiz,2),p0yz4(mvsiz,2),
1011 . p0yz5(mvsiz,2),p0yz6(mvsiz,2),p0yz7(mvsiz,2),p0yz8(mvsiz,2),
1012 . p0zy1(mvsiz,2),p0zy2(mvsiz,2),p0zy3(mvsiz,2),p0zy4(mvsiz,2),
1013 . p0zy5(mvsiz,2),p0zy6(mvsiz,2),p0zy7(mvsiz,2),p0zy8(mvsiz,2)
1014 my_real
1015 . qc(mvsiz,3,3),qgc(mvsiz,3,3),qt(mvsiz,3,3)
1016C-----------------------------------------------
1017 il = 1
1018 gbuf => elbuf_str%GBUF
1019 lbuf => elbuf_str%BUFLY(il)%LBUF(1,1,1)
1020 DO ir=1,nptr
1021 DO is=1,npts
1022 DO it=1,nptt
1023C-----------
1024 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
1025 ip = ir + ( (is-1) + (it-1)*npts )*nptr
1026C
1027 CALL s8egetpij(
1028 . px1(1,ip), px2(1,ip), px3(1,ip), px4(1,ip),
1029 . px5(1,ip), px6(1,ip), px7(1,ip), px8(1,ip),
1030 . py1(1,ip), py2(1,ip), py3(1,ip), py4(1,ip),
1031 . py5(1,ip), py6(1,ip), py7(1,ip), py8(1,ip),
1032 . pz1(1,ip), pz2(1,ip), pz3(1,ip), pz4(1,ip),
1033 . pz5(1,ip), pz6(1,ip), pz7(1,ip), pz8(1,ip),
1034 . lbuf%PIJ,llt)
1035c
1036 ENDDO
1037 ENDDO
1038 ENDDO
1039C
1040 CALL s8eselecsh(nptr,npts,nptt,nnpt,
1041 . px1, px2, px3, px4, px5, px6, px7, px8,
1042 . py1, py2, py3, py4, py5, py6, py7, py8,
1043 . pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8,
1044 . p0xy1,p0xy2,p0xy3,p0xy4,p0xy5,p0xy6,p0xy7,p0xy8,
1045 . p0yx1,p0yx2,p0yx3,p0yx4,p0yx5,p0yx6,p0yx7,p0yx8,
1046 . p0xz1,p0xz2,p0xz3,p0xz4,p0xz5,p0xz6,p0xz7,p0xz8,
1047 . p0zx1,p0zx2,p0zx3,p0zx4,p0zx5,p0zx6,p0zx7,p0zx8,
1048 . p0yz1,p0yz2,p0yz3,p0yz4,p0yz5,p0yz6,p0yz7,p0yz8,
1049 . p0zy1,p0zy2,p0zy3,p0zy4,p0zy5,p0zy6,p0zy7,p0zy8)
1050C--------conctruction trandfert matrix QQ,Q
1051 CALL setmatrans(
1052 . e1x ,e1y ,e1z ,e2x ,e2y ,e2z ,e3x ,e3y ,e3z ,
1053 . qt ,qc ,qgc )
1054 DO ir=1,nptr
1055 DO is=1,npts
1056 DO it=1,nptt
1057 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
1058 ip = ir + ( (is-1) + (it-1)*npts )*nptr
1059 CALL s8ederish2(
1060 . pxy1,pxy2,pxy3,pxy4,pxy5,pxy6,pxy7,pxy8,
1061 . pyx1,pyx2,pyx3,pyx4,pyx5,pyx6,pyx7,pyx8,
1062 . pxz1,pxz2,pxz3,pxz4,pxz5,pxz6,pxz7,pxz8,
1063 . pzx1,pzx2,pzx3,pzx4,pzx5,pzx6,pzx7,pzx8,
1064 . pyz1,pyz2,pyz3,pyz4,pyz5,pyz6,pyz7,pyz8,
1065 . pzy1,pzy2,pzy3,pzy4,pzy5,pzy6,pzy7,pzy8,
1066 . p0xy1(1,ir),p0xy2(1,ir),p0xy3(1,ir),p0xy4(1,ir),
1067 . p0xy5(1,ir),p0xy6(1,ir),p0xy7(1,ir),p0xy8(1,ir),
1068 . p0yx1(1,ir),p0yx2(1,ir),p0yx3(1,ir),p0yx4(1,ir),
1069 . p0yx5(1,ir),p0yx6(1,ir),p0yx7(1,ir),p0yx8(1,ir),
1070 . p0xz1(1,it),p0xz2(1,it),p0xz3(1,it),p0xz4(1,it),
1071 . p0xz5(1,it),p0xz6(1,it),p0xz7(1,it),p0xz8(1,it),
1072 . p0zx1(1,it),p0zx2(1,it),p0zx3(1,it),p0zx4(1,it),
1073 . p0zx5(1,it),p0zx6(1,it),p0zx7(1,it),p0zx8(1,it),
1074 . p0yz1(1,is),p0yz2(1,is),p0yz3(1,is),p0yz4(1,is),
1075 . p0yz5(1,is),p0yz6(1,is),p0yz7(1,is),p0yz8(1,is),
1076 . p0zy1(1,is),p0zy2(1,is),p0zy3(1,is),p0zy4(1,is),
1077 . p0zy5(1,is),p0zy6(1,is),p0zy7(1,is),p0zy8(1,is))
1078 CALL s8ederi_bij(
1079 . px1(1,ip), px2(1,ip), px3(1,ip), px4(1,ip),
1080 . px5(1,ip), px6(1,ip), px7(1,ip), px8(1,ip),
1081 . py1(1,ip), py2(1,ip), py3(1,ip), py4(1,ip),
1082 . py5(1,ip), py6(1,ip), py7(1,ip), py8(1,ip),
1083 . pz1(1,ip), pz2(1,ip), pz3(1,ip), pz4(1,ip),
1084 . pz5(1,ip), pz6(1,ip), pz7(1,ip), pz8(1,ip),
1085 . pxc1, pxc2, pxc3, pxc4,
1086 . pyc1, pyc2, pyc3, pyc4,
1087 . pzc1, pzc2, pzc3, pzc4,
1088 . pxy1,pxy2,pxy3,pxy4,pxy5,pxy6,pxy7,pxy8,
1089 . pyx1,pyx2,pyx3,pyx4,pyx5,pyx6,pyx7,pyx8,
1090 . pxz1,pxz2,pxz3,pxz4,pxz5,pxz6,pxz7,pxz8,
1091 . pzx1,pzx2,pzx3,pzx4,pzx5,pzx6,pzx7,pzx8,
1092 . pyz1,pyz2,pyz3,pyz4,pyz5,pyz6,pyz7,pyz8,
1093 . pzy1,pzy2,pzy3,pzy4,pzy5,pzy6,pzy7,pzy8,
1094 . bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8,
1095 . byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8,
1096 . bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8,
1097 . bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8,
1098 . byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8,
1099 . bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8,
1100 . nu )
1101C ------------return to global system for PIJ and store in PIJ(1:24->PX,PY,PZ)
1102 CALL s8ederil2g(
1103 . px1(1,ip), px2(1,ip), px3(1,ip), px4(1,ip),
1104 . px5(1,ip), px6(1,ip), px7(1,ip), px8(1,ip),
1105 . py1(1,ip), py2(1,ip), py3(1,ip), py4(1,ip),
1106 . py5(1,ip), py6(1,ip), py7(1,ip), py8(1,ip),
1107 . pz1(1,ip), pz2(1,ip), pz3(1,ip), pz4(1,ip),
1108 . pz5(1,ip), pz6(1,ip), pz7(1,ip), pz8(1,ip),
1109 . pxy1,pxy2,pxy3,pxy4,pxy5,pxy6,pxy7,pxy8,
1110 . pyx1,pyx2,pyx3,pyx4,pyx5,pyx6,pyx7,pyx8,
1111 . pxz1,pxz2,pxz3,pxz4,pxz5,pxz6,pxz7,pxz8,
1112 . pzx1,pzx2,pzx3,pzx4,pzx5,pzx6,pzx7,pzx8,
1113 . pyz1,pyz2,pyz3,pyz4,pyz5,pyz6,pyz7,pyz8,
1114 . pzy1,pzy2,pzy3,pzy4,pzy5,pzy6,pzy7,pzy8,
1115 . bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8,
1116 . byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8,
1117 . bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8,
1118 . bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8,
1119 . byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8,
1120 . bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8,
1121 . qt , qc , qgc , lbuf%PIJ,llt)
1122C
1123 ENDDO
1124 ENDDO
1125 ENDDO
1126C
1127 RETURN
1128 END
1129!||====================================================================
1130!|| s8egetpij ../starter/source/elements/solid/solide8z/s8zderi3.F
1131!||--- called by ------------------------------------------------------
1132!|| s8e_pij ../starter/source/elements/solid/solide8z/s8zderi3.F
1133!||====================================================================
1134 SUBROUTINE s8egetpij(
1135 . PX1, PX2, PX3, PX4,
1136 . PX5, PX6, PX7, PX8,
1137 . PY1, PY2, PY3, PY4,
1138 . PY5, PY6, PY7, PY8,
1139 . PZ1, PZ2, PZ3, PZ4,
1140 . PZ5, PZ6, PZ7, PZ8,
1141 . PIJ,NEL)
1142C-----------------------------------------------
1143C I m p l i c i t T y p e s
1144C-----------------------------------------------
1145#include "implicit_f.inc"
1146C-----------------------------------------------
1147C G l o b a l P a r a m e t e r s
1148C-----------------------------------------------
1149#include "mvsiz_p.inc"
1150C-----------------------------------------------
1151C D u m m y A r g u m e n t s
1152C-----------------------------------------------
1153 INTEGER NEL
1154C REAL
1155 my_real
1156 . PX1(MVSIZ),PX2(MVSIZ),PX3(MVSIZ),PX4(MVSIZ),
1157 . PX5(MVSIZ),PX6(MVSIZ),PX7(MVSIZ),PX8(MVSIZ),
1158 . PY1(MVSIZ),PY2(MVSIZ),PY3(MVSIZ),PY4(MVSIZ),
1159 . PY5(MVSIZ),PY6(MVSIZ),PY7(MVSIZ),PY8(MVSIZ),
1160 . PZ1(MVSIZ),PZ2(MVSIZ),PZ3(MVSIZ),PZ4(MVSIZ),
1161 . PZ5(MVSIZ),PZ6(MVSIZ),PZ7(MVSIZ),PZ8(MVSIZ),
1162 . PIJ(NEL,*)
1163C-----------------------------------------------
1164C L o c a l V a r i a b l e s
1165C-----------------------------------------------
1166 INTEGER I,J,K
1167
1168C-----------------------------------------------
1169 DO I=1,nel
1170 px1(i) = pij(i,1)
1171 py1(i) = pij(i,2)
1172 pz1(i) = pij(i,3)
1173 px2(i) = pij(i,4)
1174 py2(i) = pij(i,5)
1175 pz2(i) = pij(i,6)
1176 px3(i) = pij(i,7)
1177 py3(i) = pij(i,8)
1178 pz3(i) = pij(i,9)
1179 px4(i) = pij(i,10)
1180 py4(i) = pij(i,11)
1181 pz4(i) = pij(i,12)
1182 px5(i) = pij(i,13)
1183 py5(i) = pij(i,14)
1184 pz5(i) = pij(i,15)
1185 px6(i) = pij(i,16)
1186 py6(i) = pij(i,17)
1187 pz6(i) = pij(i,18)
1188 px7(i) = pij(i,19)
1189 py7(i) = pij(i,20)
1190 pz7(i) = pij(i,21)
1191 px8(i) = pij(i,22)
1192 py8(i) = pij(i,23)
1193 pz8(i) = pij(i,24)
1194 END DO
1195C-----------
1196 RETURN
1197 END
1198!||====================================================================
1199!|| s8eselecsh ../starter/source/elements/solid/solide8z/s8zderi3.F
1200!||--- called by ------------------------------------------------------
1201!|| s8e_pij ../starter/source/elements/solid/solide8z/s8zderi3.F
1202!||====================================================================
1203 SUBROUTINE s8eselecsh(NPTR,NPTS,NPTT,NNPT,
1204 . PX1, PX2, PX3, PX4, PX5, PX6, PX7, PX8,
1205 . PY1, PY2, PY3, PY4, PY5, PY6, PY7, PY8,
1206 . PZ1, PZ2, PZ3, PZ4, PZ5, PZ6, PZ7, PZ8,
1207 . P0XY1,P0XY2,P0XY3,P0XY4,P0XY5,P0XY6,P0XY7,P0XY8,
1208 . P0YX1,P0YX2,P0YX3,P0YX4,P0YX5,P0YX6,P0YX7,P0YX8,
1209 . P0XZ1,P0XZ2,P0XZ3,P0XZ4,P0XZ5,P0XZ6,P0XZ7,P0XZ8,
1210 . P0ZX1,P0ZX2,P0ZX3,P0ZX4,P0ZX5,P0ZX6,P0ZX7,P0ZX8,
1211 . P0YZ1,P0YZ2,P0YZ3,P0YZ4,P0YZ5,P0YZ6,P0YZ7,P0YZ8,
1212 . P0ZY1,P0ZY2,P0ZY3,P0ZY4,P0ZY5,P0ZY6,P0ZY7,P0ZY8)
1213C-----------------------------------------------
1214C I m p l i c i t T y p e s
1215C-----------------------------------------------
1216#include "implicit_f.inc"
1217C-----------------------------------------------
1218C G l o b a l P a r a m e t e r s
1219C-----------------------------------------------
1220#include "mvsiz_p.inc"
1221C-----------------------------------------------
1222C C o m m o n B l o c k s
1223C-----------------------------------------------
1224#include "vect01_c.inc"
1225C-----------------------------------------------
1226C D u m m y A r g u m e n t s
1227C-----------------------------------------------
1228 INTEGER NPTR,NPTS,NPTT,NNPT
1229C REAL
1230 my_real
1231 . px1(mvsiz,nnpt),px2(mvsiz,nnpt),px3(mvsiz,nnpt),px4(mvsiz,nnpt),
1232 . px5(mvsiz,nnpt),px6(mvsiz,nnpt),px7(mvsiz,nnpt),px8(mvsiz,nnpt),
1233 . py1(mvsiz,nnpt),py2(mvsiz,nnpt),py3(mvsiz,nnpt),py4(mvsiz,nnpt),
1234 . py5(mvsiz,nnpt),py6(mvsiz,nnpt),py7(mvsiz,nnpt),py8(mvsiz,nnpt),
1235 . pz1(mvsiz,nnpt),pz2(mvsiz,nnpt),pz3(mvsiz,nnpt),pz4(mvsiz,nnpt),
1236 . pz5(mvsiz,nnpt),pz6(mvsiz,nnpt),pz7(mvsiz,nnpt),pz8(mvsiz,nnpt),
1237 . p0xy1(mvsiz,2),p0xy2(mvsiz,2),p0xy3(mvsiz,2),p0xy4(mvsiz,2),
1238 . p0xy5(mvsiz,2),p0xy6(mvsiz,2),p0xy7(mvsiz,2),p0xy8(mvsiz,2),
1239 . p0yx1(mvsiz,2),p0yx2(mvsiz,2),p0yx3(mvsiz,2),p0yx4(mvsiz,2),
1240 . p0yx5(mvsiz,2),p0yx6(mvsiz,2),p0yx7(mvsiz,2),p0yx8(mvsiz,2),
1241 . p0xz1(mvsiz,2),p0xz2(mvsiz,2),p0xz3(mvsiz,2),p0xz4(mvsiz,2),
1242 . p0xz5(mvsiz,2),p0xz6(mvsiz,2),p0xz7(mvsiz,2),p0xz8(mvsiz,2),
1243 . p0zx1(mvsiz,2),p0zx2(mvsiz,2),p0zx3(mvsiz,2),p0zx4(mvsiz,2),
1244 . p0zx5(mvsiz,2),p0zx6(mvsiz,2),p0zx7(mvsiz,2),p0zx8(mvsiz,2),
1245 . p0yz1(mvsiz,2),p0yz2(mvsiz,2),p0yz3(mvsiz,2),p0yz4(mvsiz,2),
1246 . p0yz5(mvsiz,2),p0yz6(mvsiz,2),p0yz7(mvsiz,2),p0yz8(mvsiz,2),
1247 . p0zy1(mvsiz,2),p0zy2(mvsiz,2),p0zy3(mvsiz,2),p0zy4(mvsiz,2),
1248 . p0zy5(mvsiz,2),p0zy6(mvsiz,2),p0zy7(mvsiz,2),p0zy8(mvsiz,2)
1249C-----------------------------------------------
1250C L o c a l V a r i a b l e s
1251C-----------------------------------------------
1252 INTEGER I, J ,IR,IS,IT,IP,I1,I2,I3,I4
1253 INTEGER IXY(4,2),IXZ(4,2),IYZ(4,2)
1254 DATA ixy/1,5,3,7,2,6,4,8/
1255 DATA ixz/1,3,2,4,5,7,6,8/
1256 DATA iyz/1,5,2,6,3,7,4,8/
1257C------ moyen in 3 directions
1258 DO ir=1,nptr
1259 i1 = ixy(1,ir)
1260 i2 = ixy(2,ir)
1261 i3 = ixy(3,ir)
1262 i4 = ixy(4,ir)
1263 DO i=lft,llt
1264 p0yx1(i,ir)=fourth*(px1(i,i1)+px1(i,i2)+px1(i,i3)+px1(i,i4))
1265 p0yx2(i,ir)=fourth*(px2(i,i1)+px2(i,i2)+px2(i,i3)+px2(i,i4))
1266 p0yx3(i,ir)=fourth*(px3(i,i1)+px3(i,i2)+px3(i,i3)+px3(i,i4))
1267 p0yx4(i,ir)=fourth*(px4(i,i1)+px4(i,i2)+px4(i,i3)+px4(i,i4))
1268 p0yx5(i,ir)=fourth*(px5(i,i1)+px5(i,i2)+px5(i,i3)+px5(i,i4))
1269 p0yx6(i,ir)=fourth*(px6(i,i1)+px6(i,i2)+px6(i,i3)+px6(i,i4))
1270 p0yx7(i,ir)=fourth*(px7(i,i1)+px7(i,i2)+px7(i,i3)+px7(i,i4))
1271 p0yx8(i,ir)=fourth*(px8(i,i1)+px8(i,i2)+px8(i,i3)+px8(i,i4))
1272 ENDDO
1273 DO i=lft,llt
1274 p0xy1(i,ir)=fourth*(py1(i,i1)+py1(i,i2)+py1(i,i3)+py1(i,i4))
1275 p0xy2(i,ir)=fourth*(py2(i,i1)+py2(i,i2)+py2(i,i3)+py2(i,i4))
1276 p0xy3(i,ir)=fourth*(py3(i,i1)+py3(i,i2)+py3(i,i3)+py3(i,i4))
1277 p0xy4(i,ir)=fourth*(py4(i,i1)+py4(i,i2)+py4(i,i3)+py4(i,i4))
1278 p0xy5(i,ir)=fourth*(py5(i,i1)+py5(i,i2)+py5(i,i3)+py5(i,i4))
1279 p0xy6(i,ir)=fourth*(py6(i,i1)+py6(i,i2)+py6(i,i3)+py6(i,i4))
1280 p0xy7(i,ir)=fourth*(py7(i,i1)+py7(i,i2)+py7(i,i3)+py7(i,i4))
1281 p0xy8(i,ir)=fourth*(py8(i,i1)+py8(i,i2)+py8(i,i3)+py8(i,i4))
1282 ENDDO
1283 ENDDO
1284C
1285 DO it=1,nptt
1286 i1 = ixz(1,it)
1287 i2 = ixz(2,it)
1288 i3 = ixz(3,it)
1289 i4 = ixz(4,it)
1290 DO i=lft,llt
1291 p0zx1(i,it)=fourth*(px1(i,i1)+px1(i,i2)+px1(i,i3)+px1(i,i4))
1292 p0zx2(i,it)=fourth*(px2(i,i1)+px2(i,i2)+px2(i,i3)+px2(i,i4))
1293 p0zx3(i,it)=fourth*(px3(i,i1)+px3(i,i2)+px3(i,i3)+px3(i,i4))
1294 p0zx4(i,it)=fourth*(px4(i,i1)+px4(i,i2)+px4(i,i3)+px4(i,i4))
1295 p0zx5(i,it)=fourth*(px5(i,i1)+px5(i,i2)+px5(i,i3)+px5(i,i4))
1296 p0zx6(i,it)=fourth*(px6(i,i1)+px6(i,i2)+px6(i,i3)+px6(i,i4))
1297 p0zx7(i,it)=fourth*(px7(i,i1)+px7(i,i2)+px7(i,i3)+px7(i,i4))
1298 p0zx8(i,it)=fourth*(px8(i,i1)+px8(i,i2)+px8(i,i3)+px8(i,i4))
1299 ENDDO
1300 DO i=lft,llt
1301 p0xz1(i,it)=fourth*(pz1(i,i1)+pz1(i,i2)+pz1(i,i3)+pz1(i,i4))
1302 p0xz2(i,it)=fourth*(pz2(i,i1)+pz2(i,i2)+pz2(i,i3)+pz2(i,i4))
1303 p0xz3(i,it)=fourth*(pz3(i,i1)+pz3(i,i2)+pz3(i,i3)+pz3(i,i4))
1304 p0xz4(i,it)=fourth*(pz4(i,i1)+pz4(i,i2)+pz4(i,i3)+pz4(i,i4))
1305 p0xz5(i,it)=fourth*(pz5(i,i1)+pz5(i,i2)+pz5(i,i3)+pz5(i,i4))
1306 p0xz6(i,it)=fourth*(pz6(i,i1)+pz6(i,i2)+pz6(i,i3)+pz6(i,i4))
1307 p0xz7(i,it)=fourth*(pz7(i,i1)+pz7(i,i2)+pz7(i,i3)+pz7(i,i4))
1308 p0xz8(i,it)=fourth*(pz8(i,i1)+pz8(i,i2)+pz8(i,i3)+pz8(i,i4))
1309 ENDDO
1310 ENDDO
1311C
1312 DO is=1,npts
1313 i1 = iyz(1,is)
1314 i2 = iyz(2,is)
1315 i3 = iyz(3,is)
1316 i4 = iyz(4,is)
1317 DO i=lft,llt
1318 p0zy1(i,is)=fourth*(py1(i,i1)+py1(i,i2)+py1(i,i3)+py1(i,i4))
1319 p0zy2(i,is)=fourth*(py2(i,i1)+py2(i,i2)+py2(i,i3)+py2(i,i4))
1320 p0zy3(i,is)=fourth*(py3(i,i1)+py3(i,i2)+py3(i,i3)+py3(i,i4))
1321 p0zy4(i,is)=fourth*(py4(i,i1)+py4(i,i2)+py4(i,i3)+py4(i,i4))
1322 p0zy5(i,is)=fourth*(py5(i,i1)+py5(i,i2)+py5(i,i3)+py5(i,i4))
1323 p0zy6(i,is)=fourth*(py6(i,i1)+py6(i,i2)+py6(i,i3)+py6(i,i4))
1324 p0zy7(i,is)=fourth*(py7(i,i1)+py7(i,i2)+py7(i,i3)+py7(i,i4))
1325 p0zy8(i,is)=fourth*(py8(i,i1)+py8(i,i2)+py8(i,i3)+py8(i,i4))
1326 ENDDO
1327 DO i=lft,llt
1328 p0yz1(i,is)=fourth*(pz1(i,i1)+pz1(i,i2)+pz1(i,i3)+pz1(i,i4))
1329 p0yz2(i,is)=fourth*(pz2(i,i1)+pz2(i,i2)+pz2(i,i3)+pz2(i,i4))
1330 p0yz3(i,is)=fourth*(pz3(i,i1)+pz3(i,i2)+pz3(i,i3)+pz3(i,i4))
1331 p0yz4(i,is)=fourth*(pz4(i,i1)+pz4(i,i2)+pz4(i,i3)+pz4(i,i4))
1332 p0yz5(i,is)=fourth*(pz5(i,i1)+pz5(i,i2)+pz5(i,i3)+pz5(i,i4))
1333 p0yz6(i,is)=fourth*(pz6(i,i1)+pz6(i,i2)+pz6(i,i3)+pz6(i,i4))
1334 p0yz7(i,is)=fourth*(pz7(i,i1)+pz7(i,i2)+pz7(i,i3)+pz7(i,i4))
1335 p0yz8(i,is)=fourth*(pz8(i,i1)+pz8(i,i2)+pz8(i,i3)+pz8(i,i4))
1336 ENDDO
1337 ENDDO
1338C
1339 RETURN
1340 END
1341!||====================================================================
1342!|| s8ederish2 ../starter/source/elements/solid/solide8z/s8zderi3.F
1343!||--- called by ------------------------------------------------------
1344!|| s8e_pij ../starter/source/elements/solid/solide8z/s8zderi3.F
1345!||--- calls -----------------------------------------------------
1346!|| s8ea2bp8 ../starter/source/elements/solid/solide8z/s8zderi3.F
1347!||====================================================================
1348 SUBROUTINE s8ederish2(
1349 . PXY1,PXY2,PXY3,PXY4,PXY5,PXY6,PXY7,PXY8,
1350 . PYX1,PYX2,PYX3,PYX4,PYX5,PYX6,PYX7,PYX8,
1351 . PXZ1,PXZ2,PXZ3,PXZ4,PXZ5,PXZ6,PXZ7,PXZ8,
1352 . PZX1,PZX2,PZX3,PZX4,PZX5,PZX6,PZX7,PZX8,
1353 . PYZ1,PYZ2,PYZ3,PYZ4,PYZ5,PYZ6,PYZ7,PYZ8,
1354 . PZY1,PZY2,PZY3,PZY4,PZY5,PZY6,PZY7,PZY8,
1355 . P0XY1,P0XY2,P0XY3,P0XY4,P0XY5,P0XY6,P0XY7,P0XY8,
1356 . P0YX1,P0YX2,P0YX3,P0YX4,P0YX5,P0YX6,P0YX7,P0YX8,
1357 . P0XZ1,P0XZ2,P0XZ3,P0XZ4,P0XZ5,P0XZ6,P0XZ7,P0XZ8,
1358 . P0ZX1,P0ZX2,P0ZX3,P0ZX4,P0ZX5,P0ZX6,P0ZX7,P0ZX8,
1359 . P0YZ1,P0YZ2,P0YZ3,P0YZ4,P0YZ5,P0YZ6,P0YZ7,P0YZ8,
1360 . P0ZY1,P0ZY2,P0ZY3,P0ZY4,P0ZY5,P0ZY6,P0ZY7,P0ZY8)
1361C-----------------------------------------------
1362C I m p l i c i t T y p e s
1363C-----------------------------------------------
1364#include "implicit_f.inc"
1365C-----------------------------------------------
1366C D u m m y A r g u m e n t s
1367C-----------------------------------------------
1368C REAL
1369 my_real
1370 . pxy1(*),pxy2(*),pxy3(*),pxy4(*),
1371 . pxy5(*),pxy6(*),pxy7(*),pxy8(*),
1372 . pyx1(*),pyx2(*),pyx3(*),pyx4(*),
1373 . pyx5(*),pyx6(*),pyx7(*),pyx8(*),
1374 . pxz1(*),pxz2(*),pxz3(*),pxz4(*),
1375 . pxz5(*),pxz6(*),pxz7(*),pxz8(*),
1376 . pzx1(*),pzx2(*),pzx3(*),pzx4(*),
1377 . pzx5(*),pzx6(*),pzx7(*),pzx8(*),
1378 . pyz1(*),pyz2(*),pyz3(*),pyz4(*),
1379 . pyz5(*),pyz6(*),pyz7(*),pyz8(*),
1380 . pzy1(*),pzy2(*),pzy3(*),pzy4(*),
1381 . pzy5(*),pzy6(*),pzy7(*),pzy8(*),
1382 . p0xy1(*),p0xy2(*),p0xy3(*),p0xy4(*),
1383 . p0xy5(*),p0xy6(*),p0xy7(*),p0xy8(*),
1384 . p0yx1(*),p0yx2(*),p0yx3(*),p0yx4(*),
1385 . p0yx5(*),p0yx6(*),p0yx7(*),p0yx8(*),
1386 . p0xz1(*),p0xz2(*),p0xz3(*),p0xz4(*),
1387 . p0xz5(*),p0xz6(*),p0xz7(*),p0xz8(*),
1388 . p0zx1(*),p0zx2(*),p0zx3(*),p0zx4(*),
1389 . p0zx5(*),p0zx6(*),p0zx7(*),p0zx8(*),
1390 . p0yz1(*),p0yz2(*),p0yz3(*),p0yz4(*),
1391 . p0yz5(*),p0yz6(*),p0yz7(*),p0yz8(*),
1392 . p0zy1(*),p0zy2(*),p0zy3(*),p0zy4(*),
1393 . p0zy5(*),p0zy6(*),p0zy7(*),p0zy8(*)
1394C-----------------------------------------------
1395C L o c a l V a r i a b l e s
1396C-----------------------------------------------
1397 INTEGER I, J
1398C 12
1399 CALL s8ea2bp8(
1400 . pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8,
1401 . p0xy1, p0xy2, p0xy3, p0xy4,
1402 . p0xy5, p0xy6, p0xy7, p0xy8)
1403 CALL s8ea2bp8(
1404 . pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8,
1405 . p0yx1,p0yx2,p0yx3,p0yx4,
1406 . p0yx5,p0yx6,p0yx7,p0yx8)
1407 CALL s8ea2bp8(
1408 . pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8,
1409 . p0xz1,p0xz2,p0xz3,p0xz4,
1410 . p0xz5,p0xz6,p0xz7,p0xz8)
1411 CALL s8ea2bp8(
1412 . pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8,
1413 . p0zx1,p0zx2,p0zx3,p0zx4,
1414 . p0zx5,p0zx6,p0zx7,p0zx8)
1415 CALL s8ea2bp8(
1416 . pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8,
1417 . p0yz1,p0yz2,p0yz3,p0yz4,
1418 . p0yz5,p0yz6,p0yz7,p0yz8)
1419 CALL s8ea2bp8(
1420 . pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8,
1421 . p0zy1,p0zy2,p0zy3,p0zy4,
1422 . p0zy5,p0zy6,p0zy7,p0zy8)
1423C
1424 RETURN
1425 END
1426!||====================================================================
1427!|| s8ea2bp8 ../starter/source/elements/solid/solide8z/s8zderi3.F
1428!||--- called by ------------------------------------------------------
1429!|| s8ederish2 ../starter/source/elements/solid/solide8z/s8zderi3.F
1430!||====================================================================
1431 SUBROUTINE s8ea2bp8(
1432 . PX1, PX2, PX3, PX4, PX5, PX6, PX7, PX8,
1433 . P0X1, P0X2, P0X3, P0X4, P0X5, P0X6, P0X7, P0X8)
1434C-----------------------------------------------
1435C I m p l i c i t T y p e s
1436C-----------------------------------------------
1437#include "implicit_f.inc"
1438C-----------------------------------------------
1439C C o m m o n B l o c k s
1440C-----------------------------------------------
1441#include "vect01_c.inc"
1442C-----------------------------------------------
1443C D u m m y A r g u m e n t s
1444C-----------------------------------------------
1445C REAL
1446 my_real
1447 . px1(*), px2(*), px3(*), px4(*),
1448 . px5(*), px6(*), px7(*), px8(*),
1449 . p0x1(*), p0x2(*), p0x3(*), p0x4(*),
1450 . p0x5(*), p0x6(*), p0x7(*), p0x8(*)
1451C-----------------------------------------------
1452C L o c a l V a r i a b l e s
1453C-----------------------------------------------
1454 INTEGER I
1455C
1456 PX1(LFT:LLT) = p0x1(lft:llt)
1457 px2(lft:llt) = p0x2(lft:llt)
1458 px3(lft:llt) = p0x3(lft:llt)
1459 px4(lft:llt) = p0x4(lft:llt)
1460 px5(lft:llt) = p0x5(lft:llt)
1461 px6(lft:llt) = p0x6(lft:llt)
1462 px7(lft:llt) = p0x7(lft:llt)
1463 px8(lft:llt) = p0x8(lft:llt)
1464C
1465 RETURN
1466 END
1467!||====================================================================
1468!|| s8ederi_bij ../starter/source/elements/solid/solide8z/s8zderi3.F
1469!||--- called by ------------------------------------------------------
1470!|| s8e_pij ../starter/source/elements/solid/solide8z/s8zderi3.F
1471!||====================================================================
1472 SUBROUTINE s8ederi_bij(
1473 . PX1, PX2, PX3, PX4, PX5, PX6, PX7, PX8,
1474 . PY1, PY2, PY3, PY4, PY5, PY6, PY7, PY8,
1475 . PZ1, PZ2, PZ3, PZ4, PZ5, PZ6, PZ7, PZ8,
1476 . PXC1, PXC2, PXC3, PXC4,
1477 . PYC1, PYC2, PYC3, PYC4,
1478 . PZC1, PZC2, PZC3, PZC4,
1479 . PXY1,PXY2,PXY3,PXY4,PXY5,PXY6,PXY7,PXY8,
1480 . PYX1,PYX2,PYX3,PYX4,PYX5,PYX6,PYX7,PYX8,
1481 . PXZ1,PXZ2,PXZ3,PXZ4,PXZ5,PXZ6,PXZ7,PXZ8,
1482 . PZX1,PZX2,PZX3,PZX4,PZX5,PZX6,PZX7,PZX8,
1483 . PYZ1,PYZ2,PYZ3,PYZ4,PYZ5,PYZ6,PYZ7,PYZ8,
1484 . PZY1,PZY2,PZY3,PZY4,PZY5,PZY6,PZY7,PZY8,
1485 . BXY1, BXY2, BXY3, BXY4, BXY5, BXY6, BXY7, BXY8,
1486 . BYX1, BYX2, BYX3, BYX4, BYX5, BYX6, BYX7, BYX8,
1487 . BXZ1, BXZ2, BXZ3, BXZ4, BXZ5, BXZ6, BXZ7, BXZ8,
1488 . BZX1, BZX2, BZX3, BZX4, BZX5, BZX6, BZX7, BZX8,
1489 . BYZ1, BYZ2, BYZ3, BYZ4, BYZ5, BYZ6, BYZ7, BYZ8,
1490 . BZY1, BZY2, BZY3, BZY4, BZY5, BZY6, BZY7, BZY8,
1491 . NU )
1492C-----------------------------------------------
1493C I m p l i c i t T y p e s
1494C-----------------------------------------------
1495#include "implicit_f.inc"
1496C-----------------------------------------------
1497C G l o b a l P a r a m e t e r s
1498C-----------------------------------------------
1499#include "mvsiz_p.inc"
1500C-----------------------------------------------
1501C C o m m o n B l o c k s
1502C-----------------------------------------------
1503#include "vect01_c.inc"
1504C-----------------------------------------------
1505C D u m m y A r g u m e n t s
1506C-----------------------------------------------
1507C REAL
1508 my_real
1509 . px1(*), px2(*), px3(*), px4(*),
1510 . px5(*), px6(*), px7(*), px8(*),
1511 . py1(*), py2(*), py3(*), py4(*),
1512 . py5(*), py6(*), py7(*), py8(*),
1513 . pz1(*), pz2(*), pz3(*), pz4(*),
1514 . pz5(*), pz6(*), pz7(*), pz8(*),
1515 . pxy1(*),pxy2(*),pxy3(*),pxy4(*),
1516 . pxy5(*),pxy6(*),pxy7(*),pxy8(*),
1517 . pyx1(*),pyx2(*),pyx3(*),pyx4(*),
1518 . pyx5(*),pyx6(*),pyx7(*),pyx8(*),
1519 . pxz1(*),pxz2(*),pxz3(*),pxz4(*),
1520 . pxz5(*),pxz6(*),pxz7(*),pxz8(*),
1521 . pzx1(*),pzx2(*),pzx3(*),pzx4(*),
1522 . pzx5(*),pzx6(*),pzx7(*),pzx8(*),
1523 . pyz1(*),pyz2(*),pyz3(*),pyz4(*),
1524 . pyz5(*),pyz6(*),pyz7(*),pyz8(*),
1525 . pzy1(*),pzy2(*),pzy3(*),pzy4(*),
1526 . pzy5(*),pzy6(*),pzy7(*),pzy8(*),
1527 . pxc1(*), pxc2(*), pxc3(*), pxc4(*),
1528 . pyc1(*), pyc2(*), pyc3(*), pyc4(*),
1529 . pzc1(*), pzc2(*), pzc3(*), pzc4(*),
1530 . bxy1(*),bxy2(*),bxy3(*),bxy4(*),
1531 . bxy5(*),bxy6(*),bxy7(*),bxy8(*),
1532 . byx1(*),byx2(*),byx3(*),byx4(*),
1533 . byx5(*),byx6(*),byx7(*),byx8(*),
1534 . bxz1(*),bxz2(*),bxz3(*),bxz4(*),
1535 . bxz5(*),bxz6(*),bxz7(*),bxz8(*),
1536 . bzx1(*),bzx2(*),bzx3(*),bzx4(*),
1537 . bzx5(*),bzx6(*),bzx7(*),bzx8(*),
1538 . byz1(*),byz2(*),byz3(*),byz4(*),
1539 . byz5(*),byz6(*),byz7(*),byz8(*),
1540 . bzy1(*),bzy2(*),bzy3(*),bzy4(*),
1541 . bzy5(*),bzy6(*),bzy7(*),bzy8(*),
1542 . nu(*)
1543C-----------------------------------------------
1544C L o c a l V a r i a b l e s
1545C-----------------------------------------------
1546 INTEGER I, J ,ICOR
1547C REAL
1548C 12
1549 my_real
1550 . nu1(mvsiz),x1,x3,y1,y2,z2,z3,x24,y34,z14
1551C-----------------------------------------------
1552C---PX_I=X1234_I+PXC_I; PYX_I=X23_I+PXC_I; PZX_I=X12_I+PXC_I
1553C-----X2_I=0->X1_I = X12_I ; X3_I = X23_I X24_I =X1234_I-X1_I-X3_I
1554C-----verifie valeur de X4_I otherwise to be calculated
1555C---PY_I=Y1234_I+PYC_I; PXY_I=Y23_I+PYC_I; PZY_I=Y13_I+PYC_I
1556C-----Y3_I=0->Y1_I = Y13_I ; Y2_I = Y23_I Y34_I =Y1234_I-Y1_I-Y2_I
1557C---PZ_I=Z1234_I+PZC_I; PXZ_I=Z12_I+PZC_I; PYZ_I=Z13_I+PZC_I
1558C-----Z1_I=0->Z2_I = Z12_I ; Z3_I = Z13_I Z14_I =Z1234_I-Z2_I-Z3_I
1559C-------|0 BYX BZX |
1560C-------|BXY 0 BZY | same than H8A
1561C-------|BXZ BYZ 0 |
1562C-BXY_I=-nu1*X1_I-nu*X24_I
1563C-BXZ_I=-nu1*X3_I-nu*X24_I
1564C-BYX_I=-nu1*Y1_I-nu*Y34_I
1565C-BYZ_I=-nu1*Y2_I-nu*Y34_I
1566C-BZX_I=-nu1*Z2_I-nu*Z14_I
1567C-BZY_I=-nu1*Z3_I-nu*Z14_I
1568 DO i=lft,llt
1569 nu1(i) = nu(i)/(one - nu(i))
1570 ENDDO
1571 DO i=lft,llt
1572 x1 = pzx1(i)-pxc1(i)
1573 x3 = pyx1(i)-pxc1(i)
1574 y1 = pzy1(i)-pyc1(i)
1575 y2 = pxy1(i)-pyc1(i)
1576 z2 = pxz1(i)-pzc1(i)
1577 z3 = pyz1(i)-pzc1(i)
1578 x24= nu(i)*(px1(i)-pxc1(i)-x1-x3)
1579 y34= nu(i)*(py1(i)-pyc1(i)-y1-y2)
1580 z14= nu(i)*(pz1(i)-pzc1(i)-z2-z3)
1581 bxy1(i)=-nu1(i)*x1-x24
1582 bxz1(i)=-nu1(i)*x3-x24
1583 byx1(i)=-nu1(i)*y1-y34
1584 byz1(i)=-nu1(i)*y2-y34
1585 bzx1(i)=-nu1(i)*z3-z14
1586 bzy1(i)=-nu1(i)*z2-z14
1587 ENDDO
1588C
1589 DO i=lft,llt
1590 x1 = pzx2(i)-pxc2(i)
1591 x3 = pyx2(i)-pxc2(i)
1592 y1 = pzy2(i)-pyc2(i)
1593 y2 = pxy2(i)-pyc2(i)
1594 z2 = pxz2(i)-pzc2(i)
1595 z3 = pyz2(i)-pzc2(i)
1596 x24= nu(i)*(px2(i)-pxc2(i)-x1-x3)
1597 y34= nu(i)*(py2(i)-pyc2(i)-y1-y2)
1598 z14= nu(i)*(pz2(i)-pzc2(i)-z2-z3)
1599 bxy2(i)=-nu1(i)*x1-x24
1600 bxz2(i)=-nu1(i)*x3-x24
1601 byx2(i)=-nu1(i)*y1-y34
1602 byz2(i)=-nu1(i)*y2-y34
1603 bzx2(i)=-nu1(i)*z3-z14
1604 bzy2(i)=-nu1(i)*z2-z14
1605 ENDDO
1606C
1607 DO i=lft,llt
1608 x1 = pzx3(i)-pxc3(i)
1609 x3 = pyx3(i)-pxc3(i)
1610 y1 = pzy3(i)-pyc3(i)
1611 y2 = pxy3(i)-pyc3(i)
1612 z2 = pxz3(i)-pzc3(i)
1613 z3 = pyz3(i)-pzc3(i)
1614 x24= nu(i)*(px3(i)-pxc3(i)-x1-x3)
1615 y34= nu(i)*(py3(i)-pyc3(i)-y1-y2)
1616 z14= nu(i)*(pz3(i)-pzc3(i)-z2-z3)
1617 bxy3(i)=-nu1(i)*x1-x24
1618 bxz3(i)=-nu1(i)*x3-x24
1619 byx3(i)=-nu1(i)*y1-y34
1620 byz3(i)=-nu1(i)*y2-y34
1621 bzx3(i)=-nu1(i)*z3-z14
1622 bzy3(i)=-nu1(i)*z2-z14
1623 ENDDO
1624C
1625 DO i=lft,llt
1626 x1 = pzx4(i)-pxc4(i)
1627 x3 = pyx4(i)-pxc4(i)
1628 y1 = pzy4(i)-pyc4(i)
1629 y2 = pxy4(i)-pyc4(i)
1630 z2 = pxz4(i)-pzc4(i)
1631 z3 = pyz4(i)-pzc4(i)
1632 x24= nu(i)*(px4(i)-pxc4(i)-x1-x3)
1633 y34= nu(i)*(py4(i)-pyc4(i)-y1-y2)
1634 z14= nu(i)*(pz4(i)-pzc4(i)-z2-z3)
1635 bxy4(i)=-nu1(i)*x1-x24
1636 bxz4(i)=-nu1(i)*x3-x24
1637 byx4(i)=-nu1(i)*y1-y34
1638 byz4(i)=-nu1(i)*y2-y34
1639 bzx4(i)=-nu1(i)*z3-z14
1640 bzy4(i)=-nu1(i)*z2-z14
1641 ENDDO
1642C
1643 DO i=lft,llt
1644 x1 = pzx5(i)+pxc3(i)
1645 x3 = pyx5(i)+pxc3(i)
1646 y1 = pzy5(i)+pyc3(i)
1647 y2 = pxy5(i)+pyc3(i)
1648 z2 = pxz5(i)+pzc3(i)
1649 z3 = pyz5(i)+pzc3(i)
1650 x24= nu(i)*(px5(i)+pxc3(i)-x1-x3)
1651 y34= nu(i)*(py5(i)+pyc3(i)-y1-y2)
1652 z14= nu(i)*(pz5(i)+pzc3(i)-z2-z3)
1653 bxy5(i)=-nu1(i)*x1-x24
1654 bxz5(i)=-nu1(i)*x3-x24
1655 byx5(i)=-nu1(i)*y1-y34
1656 byz5(i)=-nu1(i)*y2-y34
1657 bzx5(i)=-nu1(i)*z3-z14
1658 bzy5(i)=-nu1(i)*z2-z14
1659 ENDDO
1660C
1661 DO i=lft,llt
1662 x1 = pzx6(i)+pxc4(i)
1663 x3 = pyx6(i)+pxc4(i)
1664 y1 = pzy6(i)+pyc4(i)
1665 y2 = pxy6(i)+pyc4(i)
1666 z2 = pxz6(i)+pzc4(i)
1667 z3 = pyz6(i)+pzc4(i)
1668 x24= nu(i)*(px6(i)+pxc4(i)-x1-x3)
1669 y34= nu(i)*(py6(i)+pyc4(i)-y1-y2)
1670 z14= nu(i)*(pz6(i)+pzc4(i)-z2-z3)
1671 bxy6(i)=-nu1(i)*x1-x24
1672 bxz6(i)=-nu1(i)*x3-x24
1673 byx6(i)=-nu1(i)*y1-y34
1674 byz6(i)=-nu1(i)*y2-y34
1675 bzx6(i)=-nu1(i)*z3-z14
1676 bzy6(i)=-nu1(i)*z2-z14
1677 ENDDO
1678C
1679 DO i=lft,llt
1680 x1 = pzx7(i)+pxc1(i)
1681 x3 = pyx7(i)+pxc1(i)
1682 y1 = pzy7(i)+pyc1(i)
1683 y2 = pxy7(i)+pyc1(i)
1684 z2 = pxz7(i)+pzc1(i)
1685 z3 = pyz7(i)+pzc1(i)
1686 x24= nu(i)*(px7(i)+pxc1(i)-x1-x3)
1687 y34= nu(i)*(py7(i)+pyc1(i)-y1-y2)
1688 z14= nu(i)*(pz7(i)+pzc1(i)-z2-z3)
1689 bxy7(i)=-nu1(i)*x1-x24
1690 bxz7(i)=-nu1(i)*x3-x24
1691 byx7(i)=-nu1(i)*y1-y34
1692 byz7(i)=-nu1(i)*y2-y34
1693 bzx7(i)=-nu1(i)*z3-z14
1694 bzy7(i)=-nu1(i)*z2-z14
1695 ENDDO
1696C
1697 DO i=lft,llt
1698 x1 = pzx8(i)+pxc2(i)
1699 x3 = pyx8(i)+pxc2(i)
1700 y1 = pzy8(i)+pyc2(i)
1701 y2 = pxy8(i)+pyc2(i)
1702 z2 = pxz8(i)+pzc2(i)
1703 z3 = pyz8(i)+pzc2(i)
1704 x24= nu(i)*(px8(i)+pxc2(i)-x1-x3)
1705 y34= nu(i)*(py8(i)+pyc2(i)-y1-y2)
1706 z14= nu(i)*(pz8(i)+pzc2(i)-z2-z3)
1707 bxy8(i)=-nu1(i)*x1-x24
1708 bxz8(i)=-nu1(i)*x3-x24
1709 byx8(i)=-nu1(i)*y1-y34
1710 byz8(i)=-nu1(i)*y2-y34
1711 bzx8(i)=-nu1(i)*z3-z14
1712 bzy8(i)=-nu1(i)*z2-z14
1713 ENDDO
1714C
1715 RETURN
1716 END
1717!||====================================================================
1718!|| setmatrans ../starter/source/elements/solid/solide8z/s8zderi3.f
1719!||--- called by ------------------------------------------------------
1720!|| s8e_pij ../starter/source/elements/solid/solide8z/s8zderi3.f
1721!||====================================================================
1722 SUBROUTINE setmatrans(
1723 . E1X ,E1Y ,E1Z ,E2X ,E2Y ,E2Z ,E3X ,E3Y ,E3Z ,
1724 . QT ,QC ,QGC )
1725C-----------------------------------------------
1726C I m p l i c i t T y p e s
1727C-----------------------------------------------
1728#include "implicit_f.inc"
1729C-----------------------------------------------
1730C G l o b a l P a r a m e t e r s
1731C-----------------------------------------------
1732#include "mvsiz_p.inc"
1733C-----------------------------------------------
1734C C o m m o n B l o c k s
1735C-----------------------------------------------
1736#include "vect01_c.inc"
1737C-----------------------------------------------
1738C D u m m y A r g u m e n t s
1739C-----------------------------------------------
1740C REAL
1741 my_real
1742 . e1x(mvsiz),e1y(mvsiz),e1z(mvsiz),e2x(mvsiz),
1743 . e2y(mvsiz),e2z(mvsiz),e3x(mvsiz),e3y(mvsiz),e3z(mvsiz),
1744 . qt(mvsiz,3,3),qc(mvsiz,3,3),qgc(mvsiz,3,3)
1745C-----------------------------------------------
1746C L o c a l V a r i a b l e s
1747C-----------------------------------------------
1748 INTEGER I,J,K
1749C REAL
1750 my_real
1751 . e1(mvsiz,3),e2(mvsiz,3),e3(mvsiz,3)
1752C-----------------------------------------------
1753 DO i=lft,llt
1754 qt(i,1,1)=e1x(i)
1755 qt(i,2,1)=e2x(i)
1756 qt(i,3,1)=e3x(i)
1757 qt(i,1,2)=e1y(i)
1758 qt(i,2,2)=e2y(i)
1759 qt(i,3,2)=e3y(i)
1760 qt(i,1,3)=e1z(i)
1761 qt(i,2,3)=e2z(i)
1762 qt(i,3,3)=e3z(i)
1763 ENDDO
1764 DO j= 1,3
1765 DO i=lft,llt
1766 e1(i,j)=qt(i,1,j)
1767 e2(i,j)=qt(i,2,j)
1768 e3(i,j)=qt(i,3,j)
1769 ENDDO
1770 ENDDO
1771C
1772 DO j= 1,3
1773 DO i=lft,llt
1774 qc(i,1,j)=e1(i,j)*e1(i,j)
1775 qc(i,2,j)=e2(i,j)*e2(i,j)
1776 qc(i,3,j)=e3(i,j)*e3(i,j)
1777C
1778 qgc(i,1,j)=e1(i,j)*e2(i,j)
1779 qgc(i,2,j)=e2(i,j)*e3(i,j)
1780 qgc(i,3,j)=e1(i,j)*e3(i,j)
1781 ENDDO
1782 ENDDO
1783C
1784 RETURN
1785 END
1786!||====================================================================
1787!|| bmtrl2g ../starter/source/elements/solid/solide8z/s8zderi3.F
1788!||--- called by ------------------------------------------------------
1789!|| s8el2gnj ../starter/source/elements/solid/solide8z/s8zderi3.F
1790!||--- calls -----------------------------------------------------
1791!|| ckktran3v ../starter/source/elements/solid/solide8z/s8zderi3.F
1792!||====================================================================
1793 SUBROUTINE bmtrl2g(BM ,BC ,QC ,QGC ,QT )
1794C-----------------------------------------------
1795C I m p l i c i t T y p e s
1796C-----------------------------------------------
1797#include "implicit_f.inc"
1798C-----------------------------------------------
1799C G l o b a l P a r a m e t e r s
1800C-----------------------------------------------
1801#include "mvsiz_p.inc"
1802C-----------------------------------------------
1803C C o m m o n B l o c k s
1804C-----------------------------------------------
1805#include "vect01_c.inc"
1806C-----------------------------------------------
1807C D u m m y A r g u m e n t s
1808C-----------------------------------------------
1809C REAL
1810 my_real
1811 . qc(mvsiz,3,3),qgc(mvsiz,3,3),bm(mvsiz,3,3),
1812 . bc(mvsiz,3,3),qt(mvsiz,3,3)
1813C----------[BM]g=Qc^t[BM]Qt+QGc^t[BC]Qt
1814C L o c a l V a r i a b l e s
1815C-----------------------------------------------
1816 INTEGER I,J,EP
1817C REAL
1818C----------Qc^t[BM]Qt; QGc^t[BC]Qt
1819 CALL CKKTRAN3V(LFT ,LLT ,QC ,BM ,QT )
1820 CALL CKKTRAN3V(LFT ,LLT ,QGC ,BC ,QT )
1821 DO I=1,3
1822 do j=1,3
1823 DO ep=lft,llt
1824 bm(ep,i,j)= bm(ep,i,j)+bc(ep,i,j)
1825 ENDDO
1826 ENDDO
1827 ENDDO
1828C
1829 RETURN
1830 END
1831!||====================================================================
1832!|| s8ederil2g ../starter/source/elements/solid/solide8z/s8zderi3.F
1833!||--- called by ------------------------------------------------------
1834!|| s8e_pij ../starter/source/elements/solid/solide8z/s8zderi3.F
1835!||--- calls -----------------------------------------------------
1836!|| s8el2gnj ../starter/source/elements/solid/solide8z/s8zderi3.F
1837!||====================================================================
1838 SUBROUTINE s8ederil2g(
1839 . PX1, PX2, PX3, PX4, PX5, PX6, PX7, PX8,
1840 . PY1, PY2, PY3, PY4, PY5, PY6, PY7, PY8,
1841 . PZ1, PZ2, PZ3, PZ4, PZ5, PZ6, PZ7, PZ8,
1842 . PXY1,PXY2,PXY3,PXY4,PXY5,PXY6,PXY7,PXY8,
1843 . PYX1,PYX2,PYX3,PYX4,PYX5,PYX6,PYX7,PYX8,
1844 . PXZ1,PXZ2,PXZ3,PXZ4,PXZ5,PXZ6,PXZ7,PXZ8,
1845 . PZX1,PZX2,PZX3,PZX4,PZX5,PZX6,PZX7,PZX8,
1846 . PYZ1,PYZ2,PYZ3,PYZ4,PYZ5,PYZ6,PYZ7,PYZ8,
1847 . PZY1,PZY2,PZY3,PZY4,PZY5,PZY6,PZY7,PZY8,
1848 . BXY1, BXY2, BXY3, BXY4, BXY5, BXY6, BXY7, BXY8,
1849 . BYX1, BYX2, BYX3, BYX4, BYX5, BYX6, BYX7, BYX8,
1850 . BXZ1, BXZ2, BXZ3, BXZ4, BXZ5, BXZ6, BXZ7, BXZ8,
1851 . BZX1, BZX2, BZX3, BZX4, BZX5, BZX6, BZX7, BZX8,
1852 . BYZ1, BYZ2, BYZ3, BYZ4, BYZ5, BYZ6, BYZ7, BYZ8,
1853 . BZY1, BZY2, BZY3, BZY4, BZY5, BZY6, BZY7, BZY8,
1854 . QT , QC , QGC , PIJ , NEL )
1855C-----------------------------------------------
1856C I m p l i c i t T y p e s
1857C-----------------------------------------------
1858#include "implicit_f.inc"
1859C-----------------------------------------------
1860C G l o b a l P a r a m e t e r s
1861C-----------------------------------------------
1862#include "mvsiz_p.inc"
1863C-----------------------------------------------
1864C C o m m o n B l o c k s
1865C-----------------------------------------------
1866#include "vect01_c.inc"
1867C-----------------------------------------------
1868C D u m m y A r g u m e n t s
1869C-----------------------------------------------
1870 INTEGER NEL
1871C REAL
1872 my_real
1873 . px1(*), px2(*), px3(*), px4(*),
1874 . px5(*), px6(*), px7(*), px8(*),
1875 . py1(*), py2(*), py3(*), py4(*),
1876 . py5(*), py6(*), py7(*), py8(*),
1877 . pz1(*), pz2(*), pz3(*), pz4(*),
1878 . pz5(*), pz6(*), pz7(*), pz8(*),
1879 . pxy1(*),pxy2(*),pxy3(*),pxy4(*),
1880 . pxy5(*),pxy6(*),pxy7(*),pxy8(*),
1881 . pyx1(*),pyx2(*),pyx3(*),pyx4(*),
1882 . pyx5(*),pyx6(*),pyx7(*),pyx8(*),
1883 . pxz1(*),pxz2(*),pxz3(*),pxz4(*),
1884 . pxz5(*),pxz6(*),pxz7(*),pxz8(*),
1885 . pzx1(*),pzx2(*),pzx3(*),pzx4(*),
1886 . pzx5(*),pzx6(*),pzx7(*),pzx8(*),
1887 . pyz1(*),pyz2(*),pyz3(*),pyz4(*),
1888 . pyz5(*),pyz6(*),pyz7(*),pyz8(*),
1889 . pzy1(*),pzy2(*),pzy3(*),pzy4(*),
1890 . pzy5(*),pzy6(*),pzy7(*),pzy8(*),
1891 . bxy1(*),bxy2(*),bxy3(*),bxy4(*),
1892 . bxy5(*),bxy6(*),bxy7(*),bxy8(*),
1893 . byx1(*),byx2(*),byx3(*),byx4(*),
1894 . byx5(*),byx6(*),byx7(*),byx8(*),
1895 . bxz1(*),bxz2(*),bxz3(*),bxz4(*),
1896 . bxz5(*),bxz6(*),bxz7(*),bxz8(*),
1897 . bzx1(*),bzx2(*),bzx3(*),bzx4(*),
1898 . bzx5(*),bzx6(*),bzx7(*),bzx8(*),
1899 . byz1(*),byz2(*),byz3(*),byz4(*),
1900 . byz5(*),byz6(*),byz7(*),byz8(*),
1901 . bzy1(*),bzy2(*),bzy3(*),bzy4(*),
1902 . bzy5(*),bzy6(*),bzy7(*),bzy8(*),
1903 . qt(mvsiz,9) , qc(mvsiz,9) , qgc(mvsiz,9) , pij(nel,*)
1904C-----------------------------------------------
1905C L o c a l V a r i a b l e s
1906C-----------------------------------------------
1907 INTEGER I, J ,ICOR
1908C REAL
1909C 12
1910 my_real
1911 . XA
1912C-----------------------------------------------
1913 CALL S8EL2GNJ(
1914 . PX1, PY1, PZ1,
1915 . PXY1, PYX1, PXZ1, PZX1, PYZ1, PZY1,
1916 . BXY1, BYX1, BXZ1, BZX1, BYZ1, BZY1,
1917 . QT , QC , QGC )
1918 CALL S8EL2GNJ(
1919 . px2, py2, pz2,
1920 . pxy2, pyx2, pxz2, pzx2, pyz2, pzy2,
1921 . bxy2, byx2, bxz2, bzx2, byz2, bzy2,
1922 . qt , qc , qgc )
1923 CALL s8el2gnj(
1924 . px3, py3, pz3,
1925 . pxy3, pyx3, pxz3, pzx3, pyz3, pzy3,
1926 . bxy3, byx3, bxz3, bzx3, byz3, bzy3,
1927 . qt , qc , qgc )
1928 CALL s8el2gnj(
1929 . px4, py4, pz4,
1930 . pxy4, pyx4, pxz4, pzx4, pyz4, pzy4,
1931 . bxy4, byx4, bxz4, bzx4, byz4, bzy4,
1932 . qt , qc , qgc )
1933 CALL s8el2gnj(
1934 . px5, py5, pz5,
1935 . pxy5, pyx5, pxz5, pzx5, pyz5, pzy5,
1936 . bxy5, byx5, bxz5, bzx5, byz5, bzy5,
1937 . qt , qc , qgc )
1938 CALL s8el2gnj(
1939 . px6, py6, pz6,
1940 . pxy6, pyx6, pxz6, pzx6, pyz6, pzy6,
1941 . bxy6, byx6, bxz6, bzx6, byz6, bzy6,
1942 . qt , qc , qgc )
1943 CALL s8el2gnj(
1944 . px7, py7, pz7,
1945 . pxy7, pyx7, pxz7, pzx7, pyz7, pzy7,
1946 . bxy7, byx7, bxz7, bzx7, byz7, bzy7,
1947 . qt , qc , qgc )
1948 CALL s8el2gnj(
1949 . px8, py8, pz8,
1950 . pxy8, pyx8, pxz8, pzx8, pyz8, pzy8,
1951 . bxy8, byx8, bxz8, bzx8, byz8, bzy8,
1952 . qt , qc , qgc )
1953C
1954 DO i=lft,llt
1955 pij(i,1) = px1(i)
1956 pij(i,2) = py1(i)
1957 pij(i,3) = pz1(i)
1958 pij(i,4) = px2(i)
1959 pij(i,5) = py2(i)
1960 pij(i,6) = pz2(i)
1961 pij(i,7) = px3(i)
1962 pij(i,8) = py3(i)
1963 pij(i,9) = pz3(i)
1964 pij(i,10) = px4(i)
1965 pij(i,11) = py4(i)
1966 pij(i,12) = pz4(i)
1967 pij(i,13) = px5(i)
1968 pij(i,14) = py5(i)
1969 pij(i,15) = pz5(i)
1970 pij(i,16) = px6(i)
1971 pij(i,17) = py6(i)
1972 pij(i,18) = pz6(i)
1973 pij(i,19) = px7(i)
1974 pij(i,20) = py7(i)
1975 pij(i,21) = pz7(i)
1976 pij(i,22) = px8(i)
1977 pij(i,23) = py8(i)
1978 pij(i,24) = pz8(i)
1979C
1980 pij(i,25) = bxy1(i)
1981 pij(i,26) = byx1(i)
1982 pij(i,27) = bxz1(i)
1983 pij(i,28) = bzx1(i)
1984 pij(i,29) = byz1(i)
1985 pij(i,30) = bzy1(i)
1986C
1987 pij(i,31) = bxy2(i)
1988 pij(i,32) = byx2(i)
1989 pij(i,33) = bxz2(i)
1990 pij(i,34) = bzx2(i)
1991 pij(i,35) = byz2(i)
1992 pij(i,36) = bzy2(i)
1993C
1994 pij(i,37) = bxy3(i)
1995 pij(i,38) = byx3(i)
1996 pij(i,39) = bxz3(i)
1997 pij(i,40) = bzx3(i)
1998 pij(i,41) = byz3(i)
1999 pij(i,42) = bzy3(i)
2000C
2001 pij(i,43) = bxy4(i)
2002 pij(i,44) = byx4(i)
2003 pij(i,45) = bxz4(i)
2004 pij(i,46) = bzx4(i)
2005 pij(i,47) = byz4(i)
2006 pij(i,48) = bzy4(i)
2007C
2008 pij(i,49) = bxy5(i)
2009 pij(i,50) = byx5(i)
2010 pij(i,51) = bxz5(i)
2011 pij(i,52) = bzx5(i)
2012 pij(i,53) = byz5(i)
2013 pij(i,54) = bzy5(i)
2014C
2015 pij(i,55) = bxy6(i)
2016 pij(i,56) = byx6(i)
2017 pij(i,57) = bxz6(i)
2018 pij(i,58) = bzx6(i)
2019 pij(i,59) = byz6(i)
2020 pij(i,60) = bzy6(i)
2021C
2022 pij(i,61) = bxy7(i)
2023 pij(i,62) = byx7(i)
2024 pij(i,63) = bxz7(i)
2025 pij(i,64) = bzx7(i)
2026 pij(i,65) = byz7(i)
2027 pij(i,66) = bzy7(i)
2028C
2029 pij(i,67) = bxy8(i)
2030 pij(i,68) = byx8(i)
2031 pij(i,69) = bxz8(i)
2032 pij(i,70) = bzx8(i)
2033 pij(i,71) = byz8(i)
2034 pij(i,72) = bzy8(i)
2035 ENDDO
2036C
2037 RETURN
2038 END
2039!||====================================================================
2040!|| s8el2gnj ../starter/source/elements/solid/solide8z/s8zderi3.F
2041!||--- called by ------------------------------------------------------
2042!|| s8ederil2g ../starter/source/elements/solid/solide8z/s8zderi3.F
2043!||--- calls -----------------------------------------------------
2044!|| bmtrl2g ../starter/source/elements/solid/solide8z/s8zderi3.F
2045!||====================================================================
2046 SUBROUTINE s8el2gnj(
2047 . PX1, PY1, PZ1,
2048 . PXY1, PYX1, PXZ1, PZX1, PYZ1, PZY1,
2049 . BXY1, BYX1, BXZ1, BZX1, BYZ1, BZY1,
2050 . QT , QC , QGC )
2051C-----------------------------------------------
2052C I m p l i c i t T y p e s
2053C-----------------------------------------------
2054#include "implicit_f.inc"
2055C-----------------------------------------------
2056C G l o b a l P a r a m e t e r s
2057C-----------------------------------------------
2058#include "mvsiz_p.inc"
2059C-----------------------------------------------
2060C C o m m o n B l o c k s
2061C-----------------------------------------------
2062#include "vect01_c.inc"
2063C-----------------------------------------------
2064C D u m m y A r g u m e n t s
2065C-----------------------------------------------
2066C REAL
2067 my_real
2068 . px1(*), py1(*), pz1(*),
2069 . pxy1(*), pyx1(*), pxz1(*), pzx1(*), pyz1(*), pzy1(*),
2070 . bxy1(*), byx1(*), bxz1(*), bzx1(*), byz1(*), bzy1(*),
2071 . qt(mvsiz,9) , qc(mvsiz,9) , qgc(mvsiz,9)
2072C-----------------------------------------------
2073C L o c a l V a r i a b l e s
2074C-----------------------------------------------
2075 INTEGER I, J ,ICOR
2076C REAL
2077C 12
2078 my_real
2079 . bm(mvsiz,3,3) , bc(mvsiz,3,3)
2080C-----------------------------------------------
2081C-------|PX BYX BZX | |PXY PYX 0 |
2082C-------|BXY PY BZY | :BM |PXZ 0 PZX| :BC
2083C-------|BXZ BYZ PZ | | 0 PYZ PZY|
2084C --- n1
2085 DO i=lft,llt
2086 bm(i,1,1) = px1(i)
2087 bm(i,2,2) = py1(i)
2088 bm(i,3,3) = pz1(i)
2089 bm(i,1,2) = byx1(i)
2090 bm(i,1,3) = bzx1(i)
2091 bm(i,2,1) = bxy1(i)
2092 bm(i,2,3) = bzy1(i)
2093 bm(i,3,1) = bxz1(i)
2094 bm(i,3,2) = byz1(i)
2095 ENDDO
2096 DO i=lft,llt
2097 bc(i,1,1) = py1(i)
2098 bc(i,1,2) = px1(i)
2099 bc(i,1,3) = zero
2100 bc(i,3,1) = pz1(i)
2101 bc(i,3,2) = zero
2102 bc(i,3,3) = px1(i)
2103 bc(i,2,1) = zero
2104 bc(i,2,2) = pz1(i)
2105 bc(i,2,3) = py1(i)
2106 ENDDO
2107 CALL bmtrl2g(bm ,bc ,qc ,qgc ,qt )
2108 DO i=lft,llt
2109 px1(i) = bm(i,1,1)
2110 py1(i) = bm(i,2,2)
2111 pz1(i) = bm(i,3,3)
2112 byx1(i)= bm(i,1,2)
2113 bzx1(i)= bm(i,1,3)
2114 bxy1(i)= bm(i,2,1)
2115 bzy1(i)= bm(i,2,3)
2116 bxz1(i)= bm(i,3,1)
2117 byz1(i)= bm(i,3,2)
2118 ENDDO
2119C
2120 RETURN
2121 END
2122!||====================================================================
2123!|| ckktran3v ../starter/source/elements/solid/solide8z/s8zderi3.F
2124!||--- called by ------------------------------------------------------
2125!|| bmtrl2g ../starter/source/elements/solid/solide8z/s8zderi3.F
2126!||====================================================================
2127 SUBROUTINE ckktran3v(JFT ,JLT ,VQI ,KK,VQJ)
2128C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
2129#include "implicit_f.inc"
2130#include "mvsiz_p.inc"
2131C-----------------------------------------------
2132C D U M M Y A R G U M E N T S
2133C-----------------------------------------------
2134 INTEGER JFT,JLT
2135 my_real
2136 . vqi(mvsiz,3,3), vqj(mvsiz,3,3),kk(mvsiz,3,3)
2137C-----------------------------------------------
2138C LOCAL A R G U M E N T S
2139C-----------------------------------------------
2140 INTEGER I,J,EP
2141 my_real
2142 . k(mvsiz,3,3)
2143C-----------------------------------------------
2144 DO i=1,3
2145 DO j=1,3
2146 DO ep=jft,jlt
2147 k(ep,i,j)=vqi(ep,1,i)*(kk(ep,1,1)*vqj(ep,1,j)+
2148 1 kk(ep,1,2)*vqj(ep,2,j)+kk(ep,1,3)*vqj(ep,3,j))+
2149 2 vqi(ep,2,i)*(kk(ep,2,1)*vqj(ep,1,j)+
2150 3 kk(ep,2,2)*vqj(ep,2,j)+kk(ep,2,3)*vqj(ep,3,j))+
2151 4 vqi(ep,3,i)*(kk(ep,3,1)*vqj(ep,1,j)+
2152 5 kk(ep,3,2)*vqj(ep,2,j)+kk(ep,3,3)*vqj(ep,3,j))
2153 ENDDO
2154 ENDDO
2155 ENDDO
2156C
2157 DO i=1,3
2158 DO j=1,3
2159 DO ep=jft,jlt
2160 kk(ep,i,j)= k(ep,i,j)
2161 ENDDO
2162 ENDDO
2163 ENDDO
2164C
2165 RETURN
2166 END
2167!||====================================================================
2168!|| s8ejacip3 ../starter/source/elements/solid/solide8z/s8zderi3.F
2169!||--- called by ------------------------------------------------------
2170!|| s8zinit3 ../starter/source/elements/solid/solide8z/s8zinit3.F
2171!||====================================================================
2172 SUBROUTINE s8ejacip3(
2173 . HX, HY, HZ,
2174 . CJ1,CJ2,CJ3,
2175 . CJ4,CJ5,CJ6,
2176 . CJ7,CJ8,CJ9,
2177 . JAC1,JAC2,JAC3,
2178 . JAC4,JAC5,JAC6,
2179 . JAC7,JAC8,JAC9)
2180C-----------------------------------------------
2181C I m p l i c i t T y p e s
2182C-----------------------------------------------
2183#include "implicit_f.inc"
2184C-----------------------------------------------
2185C G l o b a l P a r a m e t e r s
2186C-----------------------------------------------
2187#include "mvsiz_p.inc"
2188C-----------------------------------------------
2189C C o m m o n B l o c k s
2190C-----------------------------------------------
2191#include "vect01_c.inc"
2192C-----------------------------------------------
2193C D u m m y A r g u m e n t s
2194C-----------------------------------------------
2195C REAL
2196 my_real
2197 . hx(4,mvsiz), hy(4,mvsiz), hz(4,mvsiz),
2198 . cj1(*),cj2(*),cj3(*),
2199 . cj4(*),cj5(*),cj6(*),
2200 . cj7(*),cj8(*),cj9(*),
2201 . jac1(mvsiz,8),jac2(mvsiz,8),jac3(mvsiz,8),
2202 . jac4(mvsiz,8),jac5(mvsiz,8),jac6(mvsiz,8),
2203 . jac7(mvsiz,8),jac8(mvsiz,8),jac9(mvsiz,8)
2204C-----------------------------------------------
2205C L o c a l V a r i a b l e s
2206C-----------------------------------------------
2207 INTEGER I, J ,IP
2208C REAL
2209C 12
2210 my_real
2211 . hx1pg(mvsiz), hx2pg(mvsiz), hx3pg(mvsiz),hx4pg2(mvsiz),
2212 . hy1pg(mvsiz), hy2pg(mvsiz), hy3pg(mvsiz),hy4pg2(mvsiz),
2213 . hz1pg(mvsiz), hz2pg(mvsiz), hz3pg(mvsiz),hz4pg2(mvsiz),
2214 . pg2,pg
2215 parameter(pg=.577350269189625d0)
2216C-----------------------------------------------
2217 pg2=pg*pg
2218C
2219 DO i=lft,llt
2220 hx1pg(i)=hx(1,i)*pg
2221 hx2pg(i)=hx(2,i)*pg
2222 hx3pg(i)=hx(3,i)*pg
2223 hx4pg2(i)=hx(4,i)*pg2
2224 hy1pg(i)=hy(1,i)*pg
2225 hy2pg(i)=hy(2,i)*pg
2226 hy3pg(i)=hy(3,i)*pg
2227 hy4pg2(i)=hy(4,i)*pg2
2228 hz1pg(i)=hz(1,i)*pg
2229 hz2pg(i)=hz(2,i)*pg
2230 hz3pg(i)=hz(3,i)*pg
2231 hz4pg2(i)=hz(4,i)*pg2
2232 ENDDO
2233 ip=1
2234C - - -
2235 DO i=lft,llt
2236 jac1(i,ip)=cj1(i)-hx3pg(i)-hx2pg(i)+hx4pg2(i)
2237 jac2(i,ip)=cj2(i)-hy3pg(i)-hy2pg(i)+hy4pg2(i)
2238 jac3(i,ip)=cj3(i)-hz3pg(i)-hz2pg(i)+hz4pg2(i)
2239C
2240 jac4(i,ip)=cj4(i)-hx1pg(i)-hx3pg(i)+hx4pg2(i)
2241 jac5(i,ip)=cj5(i)-hy1pg(i)-hy3pg(i)+hy4pg2(i)
2242 jac6(i,ip)=cj6(i)-hz1pg(i)-hz3pg(i)+hz4pg2(i)
2243C
2244 jac7(i,ip)=cj7(i)-hx2pg(i)-hx1pg(i)+hx4pg2(i)
2245 jac8(i,ip)=cj8(i)-hy2pg(i)-hy1pg(i)+hy4pg2(i)
2246 jac9(i,ip)=cj9(i)-hz2pg(i)-hz1pg(i)+hz4pg2(i)
2247 ENDDO
2248 ip=2
2249C + - -
2250 DO i=lft,llt
2251 jac1(i,ip)=cj1(i)-hx3pg(i)-hx2pg(i)+hx4pg2(i)
2252 jac2(i,ip)=cj2(i)-hy3pg(i)-hy2pg(i)+hy4pg2(i)
2253 jac3(i,ip)=cj3(i)-hz3pg(i)-hz2pg(i)+hz4pg2(i)
2254C
2255 jac4(i,ip)=cj4(i)-hx1pg(i)+hx3pg(i)-hx4pg2(i)
2256 jac5(i,ip)=cj5(i)-hy1pg(i)+hy3pg(i)-hy4pg2(i)
2257 jac6(i,ip)=cj6(i)-hz1pg(i)+hz3pg(i)-hz4pg2(i)
2258C
2259 jac7(i,ip)=cj7(i)+hx2pg(i)-hx1pg(i)-hx4pg2(i)
2260 jac8(i,ip)=cj8(i)+hy2pg(i)-hy1pg(i)-hy4pg2(i)
2261 jac9(i,ip)=cj9(i)+hz2pg(i)-hz1pg(i)-hz4pg2(i)
2262 ENDDO
2263 ip=3
2264C - + -
2265 DO i=lft,llt
2266 jac1(i,ip)=cj1(i)+hx3pg(i)-hx2pg(i)-hx4pg2(i)
2267 jac2(i,ip)=cj2(i)+hy3pg(i)-hy2pg(i)-hy4pg2(i)
2268 jac3(i,ip)=cj3(i)+hz3pg(i)-hz2pg(i)-hz4pg2(i)
2269C
2270 jac4(i,ip)=cj4(i)-hx1pg(i)-hx3pg(i)+hx4pg2(i)
2271 jac5(i,ip)=cj5(i)-hy1pg(i)-hy3pg(i)+hy4pg2(i)
2272 jac6(i,ip)=cj6(i)-hz1pg(i)-hz3pg(i)+hz4pg2(i)
2273C
2274 jac7(i,ip)=cj7(i)-hx2pg(i)+hx1pg(i)-hx4pg2(i)
2275 jac8(i,ip)=cj8(i)-hy2pg(i)+hy1pg(i)-hy4pg2(i)
2276 jac9(i,ip)=cj9(i)-hz2pg(i)+hz1pg(i)-hz4pg2(i)
2277 ENDDO
2278 ip=4
2279C + + -
2280 DO i=lft,llt
2281 jac1(i,ip)=cj1(i)+hx3pg(i)-hx2pg(i)-hx4pg2(i)
2282 jac2(i,ip)=cj2(i)+hy3pg(i)-hy2pg(i)-hy4pg2(i)
2283 jac3(i,ip)=cj3(i)+hz3pg(i)-hz2pg(i)-hz4pg2(i)
2284C
2285 jac4(i,ip)=cj4(i)-hx1pg(i)+hx3pg(i)-hx4pg2(i)
2286 jac5(i,ip)=cj5(i)-hy1pg(i)+hy3pg(i)-hy4pg2(i)
2287 jac6(i,ip)=cj6(i)-hz1pg(i)+hz3pg(i)-hz4pg2(i)
2288C
2289 jac7(i,ip)=cj7(i)+hx2pg(i)+hx1pg(i)+hx4pg2(i)
2290 jac8(i,ip)=cj8(i)+hy2pg(i)+hy1pg(i)+hy4pg2(i)
2291 jac9(i,ip)=cj9(i)+hz2pg(i)+hz1pg(i)+hz4pg2(i)
2292 ENDDO
2293 ip=5
2294C - - +
2295 DO i=lft,llt
2296 jac1(i,ip)=cj1(i)-hx3pg(i)+hx2pg(i)-hx4pg2(i)
2297 jac2(i,ip)=cj2(i)-hy3pg(i)+hy2pg(i)-hy4pg2(i)
2298 jac3(i,ip)=cj3(i)-hz3pg(i)+hz2pg(i)-hz4pg2(i)
2299C
2300 jac4(i,ip)=cj4(i)+hx1pg(i)-hx3pg(i)-hx4pg2(i)
2301 jac5(i,ip)=cj5(i)+hy1pg(i)-hy3pg(i)-hy4pg2(i)
2302 jac6(i,ip)=cj6(i)+hz1pg(i)-hz3pg(i)-hz4pg2(i)
2303C
2304 jac7(i,ip)=cj7(i)-hx2pg(i)-hx1pg(i)+hx4pg2(i)
2305 jac8(i,ip)=cj8(i)-hy2pg(i)-hy1pg(i)+hy4pg2(i)
2306 jac9(i,ip)=cj9(i)-hz2pg(i)-hz1pg(i)+hz4pg2(i)
2307 ENDDO
2308 ip=6
2309C + - +
2310 DO i=lft,llt
2311 jac1(i,ip)=cj1(i)-hx3pg(i)+hx2pg(i)-hx4pg2(i)
2312 jac2(i,ip)=cj2(i)-hy3pg(i)+hy2pg(i)-hy4pg2(i)
2313 jac3(i,ip)=cj3(i)-hz3pg(i)+hz2pg(i)-hz4pg2(i)
2314C
2315 jac4(i,ip)=cj4(i)+hx1pg(i)+hx3pg(i)+hx4pg2(i)
2316 jac5(i,ip)=cj5(i)+hy1pg(i)+hy3pg(i)+hy4pg2(i)
2317 jac6(i,ip)=cj6(i)+hz1pg(i)+hz3pg(i)+hz4pg2(i)
2318C
2319 jac7(i,ip)=cj7(i)+hx2pg(i)-hx1pg(i)-hx4pg2(i)
2320 jac8(i,ip)=cj8(i)+hy2pg(i)-hy1pg(i)-hy4pg2(i)
2321 jac9(i,ip)=cj9(i)+hz2pg(i)-hz1pg(i)-hz4pg2(i)
2322 ENDDO
2323 ip=7
2324C - + +
2325 DO i=lft,llt
2326 jac1(i,ip)=cj1(i)+hx3pg(i)+hx2pg(i)+hx4pg2(i)
2327 jac2(i,ip)=cj2(i)+hy3pg(i)+hy2pg(i)+hy4pg2(i)
2328 jac3(i,ip)=cj3(i)+hz3pg(i)+hz2pg(i)+hz4pg2(i)
2329C
2330 jac4(i,ip)=cj4(i)+hx1pg(i)-hx3pg(i)-hx4pg2(i)
2331 jac5(i,ip)=cj5(i)+hy1pg(i)-hy3pg(i)-hy4pg2(i)
2332 jac6(i,ip)=cj6(i)+hz1pg(i)-hz3pg(i)-hz4pg2(i)
2333C
2334 jac7(i,ip)=cj7(i)-hx2pg(i)+hx1pg(i)-hx4pg2(i)
2335 jac8(i,ip)=cj8(i)-hy2pg(i)+hy1pg(i)-hy4pg2(i)
2336 jac9(i,ip)=cj9(i)-hz2pg(i)+hz1pg(i)-hz4pg2(i)
2337 ENDDO
2338 ip=8
2339C + + +
2340 DO i=lft,llt
2341 jac1(i,ip)=cj1(i)+hx3pg(i)+hx2pg(i)+hx4pg2(i)
2342 jac2(i,ip)=cj2(i)+hy3pg(i)+hy2pg(i)+hy4pg2(i)
2343 jac3(i,ip)=cj3(i)+hz3pg(i)+hz2pg(i)+hz4pg2(i)
2344C
2345 jac4(i,ip)=cj4(i)+hx1pg(i)+hx3pg(i)+hx4pg2(i)
2346 jac5(i,ip)=cj5(i)+hy1pg(i)+hy3pg(i)+hy4pg2(i)
2347 jac6(i,ip)=cj6(i)+hz1pg(i)+hz3pg(i)+hz4pg2(i)
2348C
2349 jac7(i,ip)=cj7(i)+hx2pg(i)+hx1pg(i)+hx4pg2(i)
2350 jac8(i,ip)=cj8(i)+hy2pg(i)+hy1pg(i)+hy4pg2(i)
2351 jac9(i,ip)=cj9(i)+hz2pg(i)+hz1pg(i)+hz4pg2(i)
2352 ENDDO
2353C
2354 RETURN
2355 END
2356!||====================================================================
2357!|| s8ederi3 ../starter/source/elements/solid/solide8z/s8zderi3.F
2358!||--- called by ------------------------------------------------------
2359!|| s8zinit3 ../starter/source/elements/solid/solide8z/s8zinit3.F
2360!||--- calls -----------------------------------------------------
2361!|| ancmsg ../starter/source/output/message/message.f
2362!||--- uses -----------------------------------------------------
2363!|| message_mod ../starter/share/message_module/message_mod.F
2364!||====================================================================
2365 SUBROUTINE s8ederi3(VOL,VEUL,GEO,WI,
2366 . JACP1,JACP2,JACP3,
2367 . JACP4,JACP5,JACP6,
2368 . JACP7,JACP8,JACP9,
2369 . SMAX, DELTAX, NGL,VOLDP)
2370 USE message_mod
2371C-----------------------------------------------
2372C I m p l i c i t T y p e s
2373C-----------------------------------------------
2374#include "implicit_f.inc"
2375C-----------------------------------------------
2376C G l o b a l P a r a m e t e r s
2377C-----------------------------------------------
2378#include "mvsiz_p.inc"
2379C-----------------------------------------------
2380C C o m m o n B l o c k s
2381C-----------------------------------------------
2382#include "vect01_c.inc"
2383#include "param_c.inc"
2384C-----------------------------------------------
2385C D u m m y A r g u m e n t s
2386C-----------------------------------------------
2387 my_real
2388 . vol(*), veul(lveul,*),geo(npropg,*),wi
2389 my_real
2390 . jacp1(*),jacp2(*),jacp3(*),
2391 . jacp4(*),jacp5(*),jacp6(*),
2392 . jacp7(*),jacp8(*),jacp9(*),
2393 . smax(*), deltax(*)
2394 INTEGER NGL(*)
2395 double precision
2396 . voldp(*), detdp
2397C-----------------------------------------------
2398C L o c a l V a r i a b l e s
2399C-----------------------------------------------
2400 INTEGER I,J
2401C 12
2402 my_real
2403 . jac1(mvsiz),jac2(mvsiz),jac3(mvsiz),
2404 . jac4(mvsiz),jac5(mvsiz),jac6(mvsiz),
2405 . jac7(mvsiz),jac8(mvsiz),jac9(mvsiz),
2406 . jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz), det(mvsiz)
2407C-----------------------------------------------
2408 jac1(lft:llt)=jacp1(lft:llt)
2409 jac2(lft:llt)=jacp2(lft:llt)
2410 jac3(lft:llt)=jacp3(lft:llt)
2411 jac4(lft:llt)=jacp4(lft:llt)
2412 jac5(lft:llt)=jacp5(lft:llt)
2413 jac6(lft:llt)=jacp6(lft:llt)
2414 jac7(lft:llt)=jacp7(lft:llt)
2415 jac8(lft:llt)=jacp8(lft:llt)
2416 jac9(lft:llt)=jacp9(lft:llt)
2417C
2418C JACOBIAN
2419C
2420 DO i=lft,llt
2421 jac_59_68(i)=jac5(i)*jac9(i)-jac6(i)*jac8(i)
2422 jac_67_49(i)=jac6(i)*jac7(i)-jac4(i)*jac9(i)
2423 jac_48_57(i)=jac4(i)*jac8(i)-jac5(i)*jac7(i)
2424 ENDDO
2425C
2426 DO i=lft,llt
2427 detdp=one_over_512*(jac1(i)*jac_59_68(i)+jac2(i)*jac_67_49(i)+jac3(i)*jac_48_57(i))
2428 det(i)=detdp
2429 voldp(i)= wi*detdp
2430 vol(i)= voldp(i)
2431 ENDDO
2432 DO 130 i=lft,llt
2433 IF(det(i)>0.0) GO TO 130
2434 CALL ancmsg(msgid=245,
2435 . msgtype=msgerror,
2436 . anmode=aninfo,
2437 . i1=ngl(i))
2438 130 CONTINUE
2439C
2440C
2441 DO i=lft,llt
2442 deltax(i)=min(deltax(i),hundred28*det(i)*smax(i))
2443 ENDDO
2444 RETURN
2445 END
2446
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
message_mod
Definition message_mod.F:1249
s8zpij_ic
subroutine s8zpij_ic(xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9, hx, hy, hz, pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4)
Definition s8zderi3.F:485
ckktran3v
subroutine ckktran3v(jft, jlt, vqi, kk, vqj)
Definition s8zderi3.F:2128
s8egetpij
subroutine s8egetpij(px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pij, nel)
Definition s8zderi3.F:1142
s8ea2bp8
subroutine s8ea2bp8(px1, px2, px3, px4, px5, px6, px7, px8, p0x1, p0x2, p0x3, p0x4, p0x5, p0x6, p0x7, p0x8)
Definition s8zderi3.F:1434
s8el2gnj
subroutine s8el2gnj(px1, py1, pz1, pxy1, pyx1, pxz1, pzx1, pyz1, pzy1, bxy1, byx1, bxz1, bzx1, byz1, bzy1, qt, qc, qgc)
Definition s8zderi3.F:2051
s8eselecsh
subroutine s8eselecsh(nptr, npts, nptt, nnpt, px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, p0xy1, p0xy2, p0xy3, p0xy4, p0xy5, p0xy6, p0xy7, p0xy8, p0yx1, p0yx2, p0yx3, p0yx4, p0yx5, p0yx6, p0yx7, p0yx8, p0xz1, p0xz2, p0xz3, p0xz4, p0xz5, p0xz6, p0xz7, p0xz8, p0zx1, p0zx2, p0zx3, p0zx4, p0zx5, p0zx6, p0zx7, p0zx8, p0yz1, p0yz2, p0yz3, p0yz4, p0yz5, p0yz6, p0yz7, p0yz8, p0zy1, p0zy2, p0zy3, p0zy4, p0zy5, p0zy6, p0zy7, p0zy8)
Definition s8zderi3.F:1213
s8ederi3
subroutine s8ederi3(vol, veul, geo, wi, jacp1, jacp2, jacp3, jacp4, jacp5, jacp6, jacp7, jacp8, jacp9, smax, deltax, ngl, voldp)
Definition s8zderi3.F:2370
s8zderi3
subroutine s8zderi3(vol, veul, geo, ksi, eta, zeta, wi, hx, hy, hz, cj1, cj2, cj3, cj4, cj5, cj6, cj7, cj8, cj9, smax, deltax, ngl, voldp)
Definition s8zderi3.F:40
setmatrans
subroutine setmatrans(e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, qt, qc, qgc)
Definition s8zderi3.F:1725
s8ederish2
subroutine s8ederish2(pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8, pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8, pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8, pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8, pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8, pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8, p0xy1, p0xy2, p0xy3, p0xy4, p0xy5, p0xy6, p0xy7, p0xy8, p0yx1, p0yx2, p0yx3, p0yx4, p0yx5, p0yx6, p0yx7, p0yx8, p0xz1, p0xz2, p0xz3, p0xz4, p0xz5, p0xz6, p0xz7, p0xz8, p0zx1, p0zx2, p0zx3, p0zx4, p0zx5, p0zx6, p0zx7, p0zx8, p0yz1, p0yz2, p0yz3, p0yz4, p0yz5, p0yz6, p0yz7, p0yz8, p0zy1, p0zy2, p0zy3, p0zy4, p0zy5, p0zy6, p0zy7, p0zy8)
Definition s8zderi3.F:1361
bmtrl2g
subroutine bmtrl2g(bm, bc, qc, qgc, qt)
Definition s8zderi3.F:1794
s8ejacip3
subroutine s8ejacip3(hx, hy, hz, cj1, cj2, cj3, cj4, cj5, cj6, cj7, cj8, cj9, jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9)
Definition s8zderi3.F:2180
s8ederil2g
subroutine s8ederil2g(px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8, pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8, pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8, pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8, pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8, pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, qt, qc, qgc, pij, nel)
Definition s8zderi3.F:1855
s8zderic3
subroutine s8zderic3(vol, hx, hy, hz, jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9, smax, det, ngl, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8)
Definition s8zderi3.F:142
s8ederi_bij
subroutine s8ederi_bij(px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4, pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8, pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8, pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8, pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8, pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8, pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, nu)
Definition s8zderi3.F:1492
s8zjac_ic
subroutine s8zjac_ic(xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9, hx, hy, hz, jac_i)
Definition s8zderi3.F:319
s8zjac_i3
subroutine s8zjac_i3(ksi, eta, zeta, wi, hx, hy, hz, cj1, cj2, cj3, cj4, cj5, cj6, cj7, cj8, cj9, jac_i, l_pij, pij, nel)
Definition s8zderi3.F:669
s8e_pij
subroutine s8e_pij(nptr, npts, nptt, nnpt, nel, pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, nu, elbuf_str)
Definition s8zderi3.F:935
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
starter
program starter
Definition starter.F:39