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