OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s20coork.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!|| s20coork ../engine/source/elements/solid/solide20/s20coork.F
25!||--- called by ------------------------------------------------------
26!|| s20ke3 ../engine/source/elements/solid/solide20/s20ke3.F
27!||--- uses -----------------------------------------------------
28!|| element_mod ../common_source/modules/elements/element_mod.F90
29!||====================================================================
30 SUBROUTINE s20coork(
31 1 X, IXS, IXS20, NPE,
32 2 XX, YY, ZZ, OFFG,
33 3 OFF, SAV, NC, NGL,
34 4 MXT, NGEO, NEL, K11,
35 5 K12, K13, K14, K15,
36 6 K16, K17, K18, K19,
37 7 K110, K111, K112, K113,
38 8 K114, K115, K116, K117,
39 9 K118, K119, K120, K22,
40 A K23, K24, K25, K26,
41 B K27, K28, K29, K210,
42 C K211, K212, K213, K214,
43 D K215, K216, K217, K218,
44 E K219, K220, K33, K34,
45 F K35, K36, K37, K38,
46 G K39, K310, K311, K312,
47 H K313, K314, K315, K316,
48 I K317, K318, K319, K320,
49 J K44, K45, K46, K47,
50 K K48, K49, K410, K411,
51 L K412, K413, K414, K415,
52 M K416, K417, K418, K419,
53 N K420, K55, K56, K57,
54 O K58, K59, K510, K511,
55 P K512, K513, K514, K515,
56 Q K516, K517, K518, K519,
57 R K520, K66, K67, K68,
58 S K69, K610, K611, K612,
59 T K613, K614, K615, K616,
60 U K617, K618, K619, K620,
61 V K77, K78, K79, K710,
62 W K711, K712, K713, K714,
63 X K715, K716, K717, K718,
64 Y K719, K720, K88, K89,
65 Z K810, K811, K812, K813,
66 1 K814, K815, K816, K817,
67 2 K818, K819, K820, K99,
68 3 K910, K911, K912, K913,
69 4 K914, K915, K916, K917,
70 5 K918, K919, K920, K1010,
71 6 K1011, K1012, K1013, K1014,
72 7 K1015, K1016, K1017, K1018,
73 8 K1019, K1020, K1111, K1112,
74 9 K1113, K1114, K1115, K1116,
75 A K1117, K1118, K1119, K1120,
76 B K1212, K1213, K1214, K1215,
77 C K1216, K1217, K1218, K1219,
78 D K1220, K1313, K1314, K1315,
79 E K1316, K1317, K1318, K1319,
80 F K1320, K1414, K1415, K1416,
81 G K1417, K1418, K1419, K1420,
82 H K1515, K1516, K1517, K1518,
83 I K1519, K1520, K1616, K1617,
84 J K1618, K1619, K1620, K1717,
85 K K1718, K1719, K1720, K1818,
86 L K1819, K1820, K1919, K1920,
87 M K2020, ISMSTR)
88 use element_mod , only : nixs
89C-----------------------------------------------
90C I m p l i c i t T y p e s
91C-----------------------------------------------
92#include "implicit_f.inc"
93C-----------------------------------------------
94C G l o b a l P a r a m e t e r s
95C-----------------------------------------------
96#include "mvsiz_p.inc"
97C-----------------------------------------------
98C C o m m o n B l o c k s
99C-----------------------------------------------
100#include "scr18_c.inc"
101C-----------------------------------------------
102C D u m m y A r g u m e n t s
103C-----------------------------------------------
104 INTEGER, INTENT(IN) :: ISMSTR
105 INTEGER NEL
106C REAL
107 my_real
108 . X(3,*),
109 . XX(MVSIZ,20), YY(MVSIZ,20), ZZ(MVSIZ,20),
110 . OFFG(*), OFF(*), SAV(NEL,30)
111 my_real
112 . K11(9,*) ,K12(9,*) ,K13(9,*) ,K14(9,*) ,K15(9,*) ,
113 . K16(9,*) ,K17(9,*) ,K18(9,*) ,K19(9,*) ,K110(9,*) ,
114 . K111(9,*) ,K112(9,*) ,K113(9,*) ,K114(9,*) ,K115(9,*) ,
115 . K116(9,*) ,K117(9,*) ,K118(9,*) ,K119(9,*) ,K120(9,*) ,
116 . K22(9,*) ,K23(9,*) ,K24(9,*) ,K25(9,*) ,K26(9,*) ,
117 . K27(9,*) ,K28(9,*) ,K29(9,*) ,K210(9,*) ,K211(9,*) ,
118 . K212(9,*) ,K213(9,*) ,K214(9,*) ,K215(9,*) ,K216(9,*) ,
119 . K217(9,*) ,K218(9,*) ,K219(9,*) ,K220(9,*) ,K33(9,*) ,
120 . K34(9,*) ,K35(9,*) ,K36(9,*) ,K37(9,*) ,K38(9,*) ,
121 . K39(9,*) ,K310(9,*) ,K311(9,*) ,K312(9,*) ,K313(9,*) ,
122 . K314(9,*) ,K315(9,*) ,K316(9,*) ,K317(9,*) ,K318(9,*) ,
123 . K319(9,*) ,K320(9,*) ,K44(9,*) ,K45(9,*) ,K46(9,*) ,
124 . K47(9,*) ,K48(9,*) ,K49(9,*) ,K410(9,*) ,K411(9,*) ,
125 . K412(9,*) ,K413(9,*) ,K414(9,*) ,K415(9,*) ,K416(9,*) ,
126 . K417(9,*) ,K418(9,*) ,K419(9,*) ,K420(9,*) ,K55(9,*) ,
127 . K56(9,*) ,K57(9,*) ,K58(9,*) ,K59(9,*) ,K510(9,*) ,
128 . K511(9,*) ,K512(9,*) ,K513(9,*) ,K514(9,*) ,K515(9,*) ,
129 . K516(9,*) ,K517(9,*) ,K518(9,*) ,K519(9,*) ,K520(9,*) ,
130 . K66(9,*) ,K67(9,*) ,K68(9,*) ,K69(9,*) ,K610(9,*) ,
131 . K611(9,*) ,K612(9,*) ,K613(9,*) ,K614(9,*) ,K615(9,*) ,
132 . K616(9,*) ,K617(9,*) ,K618(9,*) ,K619(9,*) ,K620(9,*) ,
133 . K77(9,*) ,K78(9,*) ,K79(9,*) ,K710(9,*) ,K711(9,*) ,
134 . K712(9,*) ,K713(9,*) ,K714(9,*) ,K715(9,*) ,K716(9,*) ,
135 . K717(9,*) ,K718(9,*) ,K719(9,*) ,K720(9,*) ,K88(9,*) ,
136 . K89(9,*) ,K810(9,*) ,K811(9,*) ,K812(9,*) ,K813(9,*) ,
137 . K814(9,*) ,K815(9,*) ,K816(9,*) ,K817(9,*) ,K818(9,*) ,
138 . K819(9,*) ,K820(9,*) ,K99(9,*) ,K910(9,*) ,K911(9,*) ,
139 . K912(9,*) ,K913(9,*) ,K914(9,*) ,K915(9,*) ,K916(9,*) ,
140 . K917(9,*) ,K918(9,*) ,K919(9,*) ,K920(9,*) ,K1010(9,*) ,
141 . K1011(9,*) ,K1012(9,*) ,K1013(9,*) ,K1014(9,*) ,K1015(9,*) ,
142 . K1016(9,*) ,K1017(9,*) ,K1018(9,*) ,K1019(9,*) ,K1020(9,*) ,
143 . K1111(9,*) ,K1112(9,*) ,K1113(9,*) ,K1114(9,*) ,K1115(9,*) ,
144 . K1116(9,*) ,K1117(9,*) ,K1118(9,*) ,K1119(9,*) ,K1120(9,*) ,
145 . k1212(9,*) ,k1213(9,*) ,k1214(9,*) ,k1215(9,*) ,k1216(9,*) ,
146 . k1217(9,*) ,k1218(9,*) ,k1219(9,*) ,k1220(9,*) ,k1313(9,*) ,
147 . k1314(9,*) ,k1315(9,*) ,k1316(9,*) ,k1317(9,*) ,k1318(9,*) ,
148 . k1319(9,*) ,k1320(9,*) ,k1414(9,*) ,k1415(9,*) ,k1416(9,*) ,
149 . k1417(9,*) ,k1418(9,*) ,k1419(9,*) ,k1420(9,*) ,k1515(9,*) ,
150 . k1516(9,*) ,k1517(9,*) ,k1518(9,*) ,k1519(9,*) ,k1520(9,*) ,
151 . k1616(9,*) ,k1617(9,*) ,k1618(9,*) ,k1619(9,*) ,k1620(9,*) ,
152 . k1717(9,*) ,k1718(9,*) ,k1719(9,*) ,k1720(9,*) ,k1818(9,*) ,
153 . k1819(9,*) ,k1820(9,*) ,k1919(9,*) ,k1920(9,*) ,k2020(9,*)
154 INTEGER NC(MVSIZ,20),MXT(*), NGL(*),NGEO(*)
155 INTEGER IXS(NIXS,*), IXS20(12,*),NPE
156C-----------------------------------------------
157C L o c a l V a r i a b l e s
158C-----------------------------------------------
159 INTEGER I, IPERM1(20),IPERM2(20),N,N1,N2,NN,IUN,J,MXT_1
160 DATA IPERM1/0,0,0,0,0,0,0,0,1,2,3,4,1,2,3,4,5,6,7,8/
161 DATA iperm2/0,0,0,0,0,0,0,0,2,3,4,1,5,6,7,8,6,7,8,5/
162C-----------------------------------------------
163 iun=1
164C
165 mxt_1 = ixs(1,1)
166 DO i=1,nel
167 ngeo(i) =ixs(10,i)
168 ngl(i) =ixs(11,i)
169 mxt(i) =mxt_1
170 nc(i,1) =ixs(2,i)
171 nc(i,2) =ixs(3,i)
172 nc(i,3) =ixs(4,i)
173 nc(i,4) =ixs(5,i)
174 nc(i,5) =ixs(6,i)
175 nc(i,6) =ixs(7,i)
176 nc(i,7) =ixs(8,i)
177 nc(i,8) =ixs(9,i)
178 nc(i,9) =ixs20(1,i)
179 nc(i,10) =ixs20(2,i)
180 nc(i,11) =ixs20(3,i)
181 nc(i,12) =ixs20(4,i)
182 nc(i,13) =ixs20(5,i)
183 nc(i,14) =ixs20(6,i)
184 nc(i,15) =ixs20(7,i)
185 nc(i,16) =ixs20(8,i)
186 nc(i,17) =ixs20(9,i)
187 nc(i,18) =ixs20(10,i)
188 nc(i,19) =ixs20(11,i)
189 nc(i,20) =ixs20(12,i)
190 ENDDO
191
192C----------------------------
193C NODAL COORDINATES |
194C----------------------------
195 DO n=1,npe
196
197 IF(ismstr==1.OR.(ismstr==2.AND.idtmin(1)==3))THEN
198C
199 DO i=1,nel
200 IF(offg(i)>one)THEN
201 xx(i,n)=sav(i,n)
202 yy(i,n)=sav(i,n+npe)
203 zz(i,n)=sav(i,n+2*npe)
204 off(i) = offg(i)-one
205 ELSE
206 nn = max(iun,nc(i,n))
207 xx(i,n)=x(1,nn)
208 yy(i,n)=x(2,nn)
209 zz(i,n)=x(3,nn)
210 sav(i,n)=xx(i,n)
211 sav(i,n+npe)=yy(i,n)
212 sav(i,n+npe*2)=zz(i,n)
213 off(i) = offg(i)
214 ENDIF
215 ENDDO
216C
217 ELSE
218C
219 DO i=1,nel
220 nn = max(iun,nc(i,n))
221 xx(i,n)=x(1,nn)
222 yy(i,n)=x(2,nn)
223 zz(i,n)=x(3,nn)
224 off(i) = min(one,abs(offg(i)))
225 ENDDO
226 ENDIF
227 ENDDO
228C
229
230 DO n=9,npe
231 n1=iperm1(n)
232 n2=iperm2(n)
233 DO i=1,nel
234 IF(nc(i,n)==0)THEN
235 xx(i,n) = half*(xx(i,n1)+xx(i,n2))
236 yy(i,n) = half*(yy(i,n1)+yy(i,n2))
237 zz(i,n) = half*(zz(i,n1)+zz(i,n2))
238 ENDIF
239 ENDDO
240 ENDDO
241
242
243C-----------
244 DO i=1,nel
245 DO j=1,9
246 k11(j,i) =zero
247 k12(j,i) =zero
248 k13(j,i) =zero
249 k14(j,i) =zero
250 k15(j,i) =zero
251 k16(j,i) =zero
252 k17(j,i) =zero
253 k18(j,i) =zero
254 k19(j,i) =zero
255 k110(j,i) =zero
256 k111(j,i) =zero
257 k112(j,i) =zero
258 k113(j,i) =zero
259 k114(j,i) =zero
260 k115(j,i) =zero
261 k116(j,i) =zero
262 k117(j,i) =zero
263 k118(j,i) =zero
264 k119(j,i) =zero
265 k120(j,i) =zero
266 k22(j,i) =zero
267 k23(j,i) =zero
268 k24(j,i) =zero
269 k25(j,i) =zero
270 k26(j,i) =zero
271 k27(j,i) =zero
272 k28(j,i) =zero
273 k29(j,i) =zero
274 k210(j,i) =zero
275 k211(j,i) =zero
276 k212(j,i) =zero
277 k213(j,i) =zero
278 k214(j,i) =zero
279 k215(j,i) =zero
280 k216(j,i) =zero
281 k217(j,i) =zero
282 k218(j,i) =zero
283 k219(j,i) =zero
284 k220(j,i) =zero
285 k33(j,i) =zero
286 k34(j,i) =zero
287 k35(j,i) =zero
288 k36(j,i) =zero
289 k37(j,i) =zero
290 k38(j,i) =zero
291 k39(j,i) =zero
292 k310(j,i) =zero
293 k311(j,i) =zero
294 k312(j,i) =zero
295 k313(j,i) =zero
296 k314(j,i) =zero
297 k315(j,i) =zero
298 k316(j,i) =zero
299 k317(j,i) =zero
300 k318(j,i) =zero
301 k319(j,i) =zero
302 k320(j,i) =zero
303 k44(j,i) =zero
304 k45(j,i) =zero
305 k46(j,i) =zero
306 k47(j,i) =zero
307 k48(j,i) =zero
308 k49(j,i) =zero
309 k410(j,i) =zero
310 k411(j,i) =zero
311 k412(j,i) =zero
312 k413(j,i) =zero
313 k414(j,i) =zero
314 k415(j,i) =zero
315 k416(j,i) =zero
316 k417(j,i) =zero
317 k418(j,i) =zero
318 k419(j,i) =zero
319 k420(j,i) =zero
320 k55(j,i) =zero
321 k56(j,i) =zero
322 k57(j,i) =zero
323 k58(j,i) =zero
324 k59(j,i) =zero
325 k510(j,i) =zero
326 k511(j,i) =zero
327 k512(j,i) =zero
328 k513(j,i) =zero
329 k514(j,i) =zero
330 k515(j,i) =zero
331 k516(j,i) =zero
332 k517(j,i) =zero
333 k518(j,i) =zero
334 k519(j,i) =zero
335 k520(j,i) =zero
336 k66(j,i) =zero
337 k67(j,i) =zero
338 k68(j,i) =zero
339 k69(j,i) =zero
340 k610(j,i) =zero
341 k611(j,i) =zero
342 k612(j,i) =zero
343 k613(j,i) =zero
344 k614(j,i) =zero
345 k615(j,i) =zero
346 k616(j,i) =zero
347 k617(j,i) =zero
348 k618(j,i) =zero
349 k619(j,i) =zero
350 k620(j,i) =zero
351 k77(j,i) =zero
352 k78(j,i) =zero
353 k79(j,i) =zero
354 k710(j,i) =zero
355 k711(j,i) =zero
356 k712(j,i) =zero
357 k713(j,i) =zero
358 k714(j,i) =zero
359 k715(j,i) =zero
360 k716(j,i) =zero
361 k717(j,i) =zero
362 k718(j,i) =zero
363 k719(j,i) =zero
364 k720(j,i) =zero
365 k88(j,i) =zero
366 k89(j,i) =zero
367 k810(j,i) =zero
368 k811(j,i) =zero
369 k812(j,i) =zero
370 k813(j,i) =zero
371 k814(j,i) =zero
372 k815(j,i) =zero
373 k816(j,i) =zero
374 k817(j,i) =zero
375 k818(j,i) =zero
376 k819(j,i) =zero
377 k820(j,i) =zero
378 k99(j,i) =zero
379 k910(j,i) =zero
380 k911(j,i) =zero
381 k912(j,i) =zero
382 k913(j,i) =zero
383 k914(j,i) =zero
384 k915(j,i) =zero
385 k916(j,i) =zero
386 k917(j,i) =zero
387 k918(j,i) =zero
388 k919(j,i) =zero
389 k920(j,i) =zero
390 k1010(j,i)=zero
391 k1011(j,i)=zero
392 k1012(j,i)=zero
393 k1013(j,i)=zero
394 k1014(j,i)=zero
395 k1015(j,i)=zero
396 k1016(j,i)=zero
397 k1017(j,i)=zero
398 k1018(j,i)=zero
399 k1019(j,i)=zero
400 k1020(j,i)=zero
401 k1111(j,i)=zero
402 k1112(j,i)=zero
403 k1113(j,i)=zero
404 k1114(j,i)=zero
405 k1115(j,i)=zero
406 k1116(j,i)=zero
407 k1117(j,i)=zero
408 k1118(j,i)=zero
409 k1119(j,i)=zero
410 k1120(j,i)=zero
411 k1212(j,i)=zero
412 k1213(j,i)=zero
413 k1214(j,i)=zero
414 k1215(j,i)=zero
415 k1216(j,i)=zero
416 k1217(j,i)=zero
417 k1218(j,i)=zero
418 k1219(j,i)=zero
419 k1220(j,i)=zero
420 k1313(j,i)=zero
421 k1314(j,i)=zero
422 k1315(j,i)=zero
423 k1316(j,i)=zero
424 k1317(j,i)=zero
425 k1318(j,i)=zero
426 k1319(j,i)=zero
427 k1320(j,i)=zero
428 k1414(j,i)=zero
429 k1415(j,i)=zero
430 k1416(j,i)=zero
431 k1417(j,i)=zero
432 k1418(j,i)=zero
433 k1419(j,i)=zero
434 k1420(j,i)=zero
435 k1515(j,i)=zero
436 k1516(j,i)=zero
437 k1517(j,i)=zero
438 k1518(j,i)=zero
439 k1519(j,i)=zero
440 k1520(j,i)=zero
441 k1616(j,i)=zero
442 k1617(j,i)=zero
443 k1618(j,i)=zero
444 k1619(j,i)=zero
445 k1620(j,i)=zero
446 k1717(j,i)=zero
447 k1718(j,i)=zero
448 k1719(j,i)=zero
449 k1720(j,i)=zero
450 k1818(j,i)=zero
451 k1819(j,i)=zero
452 k1820(j,i)=zero
453 k1919(j,i)=zero
454 k1920(j,i)=zero
455 k2020(j,i)=zero
456 ENDDO
457 ENDDO
458C-----------
459 RETURN
460 END
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
s20coork
subroutine s20coork(x, ixs, ixs20, npe, xx, yy, zz, offg, off, sav, nc, ngl, mxt, ngeo, nel, k11, k12, k13, k14, k15, k16, k17, k18, k19, k110, k111, k112, k113, k114, k115, k116, k117, k118, k119, k120, k22, k23, k24, k25, k26, k27, k28, k29, k210, k211, k212, k213, k214, k215, k216, k217, k218, k219, k220, k33, k34, k35, k36, k37, k38, k39, k310, k311, k312, k313, k314, k315, k316, k317, k318, k319, k320, k44, k45, k46, k47, k48, k49, k410, k411, k412, k413, k414, k415, k416, k417, k418, k419, k420, k55, k56, k57, k58, k59, k510, k511, k512, k513, k514, k515, k516, k517, k518, k519, k520, k66, k67, k68, k69, k610, k611, k612, k613, k614, k615, k616, k617, k618, k619, k620, k77, k78, k79, k710, k711, k712, k713, k714, k715, k716, k717, k718, k719, k720, k88, k89, k810, k811, k812, k813, k814, k815, k816, k817, k818, k819, k820, k99, k910, k911, k912, k913, k914, k915, k916, k917, k918, k919, k920, k1010, k1011, k1012, k1013, k1014, k1015, k1016, k1017, k1018, k1019, k1020, k1111, k1112, k1113, k1114, k1115, k1116, k1117, k1118, k1119, k1120, k1212, k1213, k1214, k1215, k1216, k1217, k1218, k1219, k1220, k1313, k1314, k1315, k1316, k1317, k1318, k1319, k1320, k1414, k1415, k1416, k1417, k1418, k1419, k1420, k1515, k1516, k1517, k1518, k1519, k1520, k1616, k1617, k1618, k1619, k1620, k1717, k1718, k1719, k1720, k1818, k1819, k1820, k1919, k1920, k2020, ismstr)
Definition s20coork.F:88