OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
weights_p4linux964_spmd_sse3.inc
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! Material & element costs
24! optimized for Sandy Bridge machine (AVX)
25! Used machine : Pajero --> Intel(R) Xeon(R) CPU E5-4650 0 @ 2.70GHz
26
27
28! temps de reference correspond au temps /elt/materiaux loi2 0 pt integ
29 TPSREF = .0000007658
30
31! -------------
32! SOLIDES
33! -------------
34
35
36! Tableau SOL1TNL
37! ---------------
38!
39! Temps par elt/cycle pour les materiaux solides
40! SOL1TNL(I,J)=
41! 1er Champ : numero du materiaux
42! 2eme Option particuliere
43!
44
45! Tableau SOL8TNL
46! ---------------
47!
48! SOL8TNL(I,J)=
49! 1er Champ : numero du materiux
50! 2eme Option particuliere
51!
52
53! Tableau SOLTELT
54! ---------------
55!
56! Temps de traitement de l'hourglas et des routines elts /elt/cycle
57!
58! SOLTELT(I)=
59! 1er Champ : option d'hourglas
60
61 DO J=1,3
62 DO I=0,MAXLAW
63 SOL1TNL(I,J) = .0000001512
64 SOL8TNL(I,J) = .0000006297
65 TET4TNL(I,J) = .0000001542
66 TET10TNL(I,J)= .0000005545
67 ENDDO
68 ENDDO
69! special option for solid
70 SOL_OPTION(1:10,1:5) = 0.
71! special option for non-local
72 NLOCAL_OPTION(1:10) = 0.
73
74! LAW 01
75 SOL1TNL(1,1) = .0000001301
76 SOL8TNL(1,1) = .0000003212
77 TET4TNL(1,1) = .0000001229
78 TET10TNL(1,1)= .0000004662
79! LAW 01 + ISMSTR = 10 or 12
80 SOL1TNL(1,2) = .0000002254
81 SOL8TNL(1,2) = .0000018500
82 TET4TNL(1,2) = .0000002254
83 TET10TNL(1,2)= .0000009016
84! LAW 02
85 SOL1TNL(2,1) = .0000001512
86 SOL8TNL(2,1) = .0000006297
87 TET4TNL(2,1) = .0000001542
88 TET10TNL(2,1)= .0000005545
89! LAW 02 ICC=1
90 SOL1TNL(2,2) = .0000002051
91 SOL8TNL(2,2) = .0000010080
92 TET4TNL(2,2) = .0000001820
93 TET10TNL(2,2)= .0000007670
94! LAW 03
95 SOL1TNL(3,1) = .0000002145
96 SOL8TNL(3,1) = .0000006086
97 TET4TNL(3,1) = .0000002292
98 TET10TNL(3,1)= .0000009112
99! LAW 04
100 SOL1TNL(4,1) = .0000001555
101 SOL8TNL(4,1) = .0000006050
102 TET4TNL(4,1) = .0000002547
103 TET10TNL(4,1)= .0000009900
104! LAW 06
105 SOL1TNL(6,1) = .0000002483
106 SOL8TNL(6,1) = .0000006004
107 TET4TNL(6,1) = .0000002412
108 TET10TNL(6,1)= .0000007858
109! LAW 10
110 SOL1TNL(10,1) = .0000002024
111 SOL8TNL(10,1) = .0000013650
112 TET4TNL(10,1) = .0000001815
113 TET10TNL(10,1)= .0000006534
114! LAW 12
115 SOL1TNL(12,1) = .0000003941
116 SOL8TNL(12,1) = .0000038700
117 TET4TNL(12,1) = .0000003670
118 TET10TNL(12,1)= .0000013560
119! LAW 21
120 SOL1TNL(21,1) = .0000001872
121 SOL8TNL(21,1) = .0000013840
122 TET4TNL(21,1) = .0000001811
123 TET10TNL(21,1)= .0000006885
124! LAW 22
125 SOL1TNL(22,1) = .0000001854
126 SOL8TNL(22,1) = .0000013990
127 TET4TNL(22,1) = .0000001663
128 TET10TNL(22,1)= .0000007349
129! LAW 23
130 SOL1TNL(23,1) = .0000001890
131 SOL8TNL(23,1) = .0000014360
132 TET4TNL(23,1) = .0000001850
133 TET10TNL(23,1)= .0000007381
134! LAW 24
135 SOL1TNL(24,1) = .0000006082
136 SOL8TNL(24,1) = .0000042810
137 TET4TNL(24,1) = .0000004260
138 TET10TNL(24,1)= .0000020050
139! LAW 25
140 SOL1TNL(25,1) = .0000005039
141 SOL8TNL(25,1) = .0000040300
142 TET4TNL(25,1) = .0000003445
143 TET10TNL(25,1)= .0000020100
144! LAW 25
145 SOL1TNL(25,1) = .0000006250
146 SOL8TNL(25,1) = .0000050290
147 TET4TNL(25,1) = .0000006250
148 TET10TNL(25,1)= .0000025400
149! LAW 25 Crasurv
150 SOL1TNL(25,2) = .0000007571
151 SOL8TNL(25,2) = .0000060500
152 TET4TNL(25,2) = .0000007571
153 TET10TNL(25,2)= .0000030280
154! LAW 28
155 SOL1TNL(28,1) = .0000007328
156 SOL8TNL(28,1) = .0000058600
157 TET4TNL(28,1) = .0000007328
158 TET10TNL(28,1)= .0000029300
159! LAW 33 0 fct
160 SOL1TNL(33,1) = .0000004453
161 SOL8TNL(33,1) = .0000021870
162 TET4TNL(33,1) = .0000004027
163 TET10TNL(33,1)= .0000015750
164! LAW 33 1 fct
165 SOL1TNL(33,2) = .0000006103
166 SOL8TNL(33,2) = .0000025320
167 TET4TNL(33,2) = .0000005659
168 TET10TNL(33,2)= .0000018270
169! LAW 34
170 SOL1TNL(34,1) = .0000002516
171 SOL8TNL(34,1) = .0000012240
172 TET4TNL(34,1) = .0000002258
173 TET10TNL(34,1)= .0000008617
174! LAW 35
175 SOL1TNL(35,1) = .0000003070
176 SOL8TNL(35,1) = .0000017020
177 TET4TNL(35,1) = .0000002691
178 TET10TNL(35,1)= .0000010870
179! LAW 36 1 fonction
180 SOL1TNL(36,1) = .0000003096
181 SOL8TNL(36,1) = .0000016040
182 TET4TNL(36,1) = .0000002877
183 TET10TNL(36,1)= .0000009584
184! LAW 36 5 fonctions
185 SOL1TNL(36,2) = .0000003358
186 SOL8TNL(36,2) = .0000016420
187 TET4TNL(36,2) = .0000002724
188 TET10TNL(36,2)= .0000010120
189! LAW 36 10 fonctions
190 SOL1TNL(36,3) = .0000003702
191 SOL8TNL(36,3) = .0000016950
192 TET4TNL(36,3) = .0000002854
193 TET10TNL(36,3)= .0000010970
194! LAW 38
195 SOL1TNL(38,1) = .0000008197
196 SOL8TNL(38,1) = .0000057300
197 TET4TNL(38,1) = .0000006935
198 TET10TNL(38,1)= .0000025160
199! LAW 38 + ISMSTR=10 or 12
200 SOL1TNL(38,2) = .0000009800
201 SOL8TNL(38,2) = .0000080100
202 TET4TNL(38,2) = .0000009800
203 TET10TNL(38,2)= .0000037800
204! LAW 40
205 SOL1TNL(40,1) = .0000004917
206 SOL8TNL(40,1) = .0000030250
207 TET4TNL(40,1) = .0000004015
208 TET10TNL(40,1)= .0000013700
209! LAW 41
210 SOL1TNL(41,1) = .0000004092
211 TET4TNL(41,1) = .0000004068
212 TET10TNL(41,1)= .0000017520
213! LAW 42
214 SOL1TNL(42,1) = .0000006159
215 SOL8TNL(42,1) = .0000041410
216 TET4TNL(42,1) = .0000006126
217 TET10TNL(42,1)= .0000023340
218C M = 1
219 SOL1TNL(42,2) = .0000006763
220 SOL8TNL(42,2) = .0000045820
221 TET4TNL(42,2) = .0000006537
222 TET10TNL(42,2)= .0000023430
223C M = 2
224 SOL1TNL(42,3) = .0000007634
225 SOL8TNL(42,3) = .0000051210
226 TET4TNL(42,3) = .0000006823
227 TET10TNL(42,3)= .0000027020
228! LAW 44
229 SOL1TNL(44,1) = .0000003088
230 SOL8TNL(44,1) = .0000016320
231 TET4TNL(44,1) = .0000002999
232 TET10TNL(44,1)= .0000010160
233! LAW 44
234 SOL1TNL(44,1) = .0000004850
235 SOL8TNL(44,1) = .0000037800
236 TET4TNL(44,1) = .0000004850
237 TET10TNL(44,1)= .0000018700
238! LAW 48
239 SOL1TNL(48,1) = .0000002994
240 SOL8TNL(48,1) = .0000015830
241 TET4TNL(48,1) = .0000002660
242 TET10TNL(48,1)= .0000009488
243! LAW 50
244 SOL1TNL(50,1) = .0000004482
245 SOL8TNL(50,1) = .0000027100
246 TET4TNL(50,1) = .0000004178
247 TET10TNL(50,1)= .0000018130
248! LAW 52
249 SOL1TNL(52,1) = .0000002826
250 SOL8TNL(52,1) = .0000015560
251 TET4TNL(52,1) = .0000002650
252 TET10TNL(52,1)= .0000010590
253! LAW 53
254 SOL1TNL(53,1) = .0000004006
255 SOL8TNL(53,1) = .0000022760
256 TET4TNL(53,1) = .0000003889
257 TET10TNL(53,1)= .0000014540
258! LAW 59 (user solid property) IGT=43
259 SOL1TNL(59,1) = .0000002207
260 SOL8TNL(59,1) = .0000002207
261 TET4TNL(59,1) = .0000002207
262 TET10TNL(59,1)= .0000002207
263! LAW 60
264 SOL1TNL(60,1) = .0000004239
265 SOL8TNL(60,1) = .0000024010
266 TET4TNL(60,1) = .0000003718
267 TET10TNL(60,1)= .0000014290
268! LAW 62
269 ! M = 0
270 SOL1TNL(62,1) = .0000009680
271 SOL8TNL(62,1) = .0000068900
272 TET4TNL(62,1) = .0000009896
273 TET10TNL(62,1)= .0000039430
274 ! M = 1
275 SOL1TNL(62,2) = .0000010850
276 SOL8TNL(62,2) = .0000078600
277 TET4TNL(62,2) = .0000010760
278 TET10TNL(62,2)= .0000043120
279 ! M = 2
280 SOL1TNL(62,3) = .0000011420
281 SOL8TNL(62,3) = .0000084600
282 TET4TNL(62,3) = .0000011150
283 TET10TNL(62,3)= .0000045435
284! LAW 65
285 SOL1TNL(65,1) = .0000004282
286 SOL8TNL(65,1) = .0000030270
287 TET4TNL(65,1) = .0000003897
288 TET10TNL(65,1)= .0000015160
289! LAW 66
290 SOL1TNL(66,1) = .0000005367
291 SOL8TNL(66,1) = .0000031200
292 TET4TNL(66,1) = .0000003860
293 TET10TNL(66,1)= .0000015470
294! LAW 68
295 SOL1TNL(68,1) = .0000006413
296! SOL8TNL(68,1) = .0000028940
297 TET4TNL(68,1) = .0000006475
298 TET10TNL(68,1)= .0000025900
299! LAW 69
300 ! M = 0 (no visc)
301 SOL1TNL(69,1) = .0000005245
302 SOL8TNL(69,1) = .0000040350
303 TET4TNL(69,1) = .0000005135
304 TET10TNL(69,1)= .0000020020
305 ! M = 1 (visc)
306 SOL1TNL(69,2) = .0000006120
307 SOL8TNL(69,2) = .0000048180
308 TET4TNL(69,2) = .0000006347
309 TET10TNL(69,2)= .0000024700
310 ! M = 2 (visc)
311 SOL1TNL(69,3) = .0000006370
312 SOL8TNL(69,3) = .0000051600
313 TET4TNL(69,3) = .0000006740
314 TET10TNL(69,3)= .0000025500
315! LAW 70
316! 1 funct / 22 points
317 SOL1TNL(70,1) = .0000003413
318 SOL8TNL(70,1) = .0000026830
319 TET4TNL(70,1) = .0000003395
320 TET10TNL(70,1)= .0000013970
321! overcost number of point
322 SOL1TNL(70,2) = .0000000020
323 SOL8TNL(70,2) = .0000000160
324 TET4TNL(70,2) = .0000000018
325 TET10TNL(70,2)= .0000000072
326! overcost number of function
327 SOL1TNL(70,3) = .0000000010
328 SOL8TNL(70,3) = .0000000078
329 TET4TNL(70,3) = .0000000009
330 TET10TNL(70,3)= .0000000036
331! LAW 71
332 SOL1TNL(71,1) = .0000006332
333 SOL8TNL(71,1) = .0000021590
334 TET4TNL(71,1) = .0000004840
335 TET10TNL(71,1)= .0000018860
336! LAW 72
337 SOL1TNL(72,1) = .0000004335
338 SOL8TNL(72,1) = .0000032810
339 TET4TNL(72,1) = .0000003947
340 TET10TNL(72,1)= .0000016170
341! LAW 74
342 SOL1TNL(74,1) = .0000003704
343 SOL8TNL(74,1) = .0000027610
344 TET4TNL(74,1) = .0000003402
345 TET10TNL(74,1)= .0000013590
346! LAW 76
347 SOL1TNL(76,1) = .0000002745
348 SOL8TNL(76,1) = .0000021280
349 TET4TNL(76,1) = .0000002739
350 TET10TNL(76,1)= .0000011500
351! LAW 77
352 SOL1TNL(77,1) = .0000004887
353 SOL8TNL(77,1) = .0000038100
354 TET4TNL(77,1) = .0000005200
355 TET10TNL(77,1)= .0000020120
356! LAW 78
357 SOL1TNL(78,1) = .0000003087
358 SOL8TNL(78,1) = .0000023090
359 TET4TNL(78,1) = .0000002902
360 TET10TNL(78,1)= .0000012520
361! LAW 79
362 SOL1TNL(79,1) = .0000003060
363 SOL8TNL(79,1) = .0000020220
364 TET4TNL(79,1) = .0000002684
365 TET10TNL(79,1)= .0000011400
366! LAW 81
367 SOL1TNL(81,1) = .0000004223
368 SOL8TNL(81,1) = .0000031430
369 TET4TNL(81,1) = .0000002750
370 TET10TNL(81,1)= .0000012470
371! LAW 82
372 ! M = 1
373 SOL1TNL(82,1) = .0000005645
374 SOL8TNL(82,1) = .0000043090
375 TET4TNL(82,1) = .0000005553
376 TET10TNL(82,1)= .0000020840
377 ! M = 2
378 SOL1TNL(82,2) = .0000006735
379 SOL8TNL(82,2) = .0000045120
380 TET4TNL(82,2) = .0000006743
381 TET10TNL(82,2)= .0000025890
382 ! M = 3
383 SOL1TNL(82,3) = .0000007826
384 SOL8TNL(82,3) = .0000052410
385 TET4TNL(82,3) = .0000008028
386 TET10TNL(82,3)= .0000030570
387! LAW 83 (user solid property) IGT=43
388 SOL1TNL(83,1) = .0000013690
389 SOL8TNL(83,1) = .0000013690
390 TET4TNL(83,1) = .0000013690
391 TET10TNL(83,1)= .0000013690
392! LAW 84
393 SOL1TNL(84,1) = .0000002771
394 SOL8TNL(84,1) = .0000018140
395 TET4TNL(84,1) = .0000002489
396 TET10TNL(84,1)= .0000007861
397! LAW 88
398 SOL1TNL(88,1) = .0000080930
399 SOL8TNL(88,1) = .0000600000
400 TET4TNL(88,1) = .0000080930
401 TET10TNL(88,1)= .0000302000
402! LAW 90
403 SOL1TNL(90,1) = .0000011200
404 SOL8TNL(90,1) = .0000089600
405 TET4TNL(90,1) = .0000011200
406 TET10TNL(90,1)= .0000044800
407! LAW 90 + ISMSTR=10 or 12
408 SOL1TNL(90,1) = .0000014000
409 SOL8TNL(90,1) = .0000096600
410 TET4TNL(90,1) = .0000014200
411 TET10TNL(90,1)= .0000047000
412! LAW 92
413 SOL1TNL(92,1) = .0000009800
414 SOL8TNL(92,1) = .0000088000
415 TET4TNL(92,1) = .0000009800
416 TET10TNL(92,1)= .0000041400
417! LAW 92 + ISMSTR=10 or 12
418 SOL1TNL(92,2) = .0000011400
419 SOL8TNL(92,2) = .0000090000
420 TET4TNL(92,2) = .0000011400
421 TET10TNL(92,2)= .0000044000
422! LAW 94
423 SOL1TNL(94,1) = .0000006230
424 SOL8TNL(94,1) = .0000051800
425 TET4TNL(94,1) = .0000006230
426 TET10TNL(94,1)= .0000024900
427! LAW 94 + ISMSTR=10 or 12
428 SOL1TNL(94,2) = .0000007010
429 SOL8TNL(94,2) = .0000056800
430 TET4TNL(94,2) = .0000007010
431 TET10TNL(94,2)= .0000027900
432! LAW 95
433 SOL1TNL(95,1) = .0000005950
434 SOL8TNL(95,1) = .0000045150
435 TET4TNL(95,1) = .0000005905
436 TET10TNL(95,1)= .0000021230
437! LAW 97
438 SOL1TNL(97,1) = .0000003053
439 SOL8TNL(97,1) = .0000022300
440 TET4TNL(97,1) = .0000002500
441 TET10TNL(97,1)= .0000009410
442! LAW 100 : material cost = base cost + viscosity cost + N * network cost
443! LAW 100 (base)
444 SOL1TNL(100,1) = .0000002280
445 SOL8TNL(100,1) = .0000002280
446 TET4TNL(100,1) = .0000002280
447 TET10TNL(100,1)= .0000002280
448
449! LAW 100 (viscosity)
450 SOL1TNL(100,2) = .0000001390
451 SOL8TNL(100,2) = .0000001390
452 TET4TNL(100,2) = .0000001390
453 TET10TNL(100,2)= .0000001390
454
455! LAW 100 (network)
456 SOL1TNL(100,3) = .0000002000
457 SOL8TNL(100,3) = .0000002000
458 TET4TNL(100,3) = .0000002000
459 TET10TNL(100,3)= .0000002000
460! LAW 101
461 SOL1TNL(101,1) = .0000090800
462 SOL8TNL(101,1) = .0000749000
463 TET4TNL(101,1) = .0000090800
464 TET10TNL(101,1)= .0000359000
465! LAW 102
466 SOL1TNL(102,1) = .0000002770
467 SOL8TNL(102,1) = .0000020330
468 TET4TNL(102,1) = .0000002560
469 TET10TNL(102,1)= .0000010010
470! LAW 103
471 SOL1TNL(103,1) = .0000003287
472 SOL8TNL(103,1) = .0000023320
473 TET4TNL(103,1) = .0000003141
474 TET10TNL(103,1)= .0000011940
475! LAW 104
476 ! Nice algo
477 SOL1TNL(104,1) = .0000002150
478 SOL8TNL(104,1) = .0000002150
479 TET4TNL(104,1) = .0000002150
480 TET10TNL(104,1) = .0000002150
481 ! Newton algo
482 SOL1TNL(104,2) = .0000003200
483 SOL8TNL(104,2) = .0000003200
484 TET4TNL(104,2) = .0000003200
485 TET10TNL(104,2) = .0000003200
486 ! iloc=1
487 SOL_OPTION(1,1) = .0000000600
488 ! iloc=2
489 SOL_OPTION(1,2) = .0000000800
490 ! iloc=3
491 SOL_OPTION(1,3) = .0000001100
492! LAW 106
493! SOL1TNL(106,1) = .0000002493
494! SOL8TNL(106,1) = .0000018090
495! TET4TNL(106,1) = .0000002246
496! TET10TNL(106,1)= .0000009194
497! LAW 187
498 SOL1TNL(187,1) = .0000003805
499 SOL8TNL(187,1) = .0000029720
500 TET4TNL(187,1) = .0000003709
501 TET10TNL(187,1)= .0000021430
502
503 IF (IPARI0.EQ.1) THEN
504! SOLTELT :
505! 1 --> ISOLID = 1
506! 2 --> ISOLID = 2
507! 3 --> ISOLID = 24
508! 4 --> ISOLID = 12
509! 5 --> ISOLID = 14 HA8 with NPT=2*2*2 (fixed part of element cost)
510! 6 --> ISOLID = 14 HA8 (overcost for NPT > 2*2*2)
511! HA8 cost : COST = HA8 (222) cost + (NPT-8) * overcost
512! 7 --> ISOLID = 17 H8C
513! 8 --> IGTYP = 43
514! 9 --> ISOLID = 18 H8C
515! 10--> IGTYP = Thickshell 20 or 21 or 22
516! NLOCAL_OPTION
517! 1 --> sol
518! 2 --> shell
519!
520! TET4ELT :
521! 1 --> TETRA 4
522! 2 --> TETRA 4 with IROT = 1
523! TETRA10ELT
524! 1 --> TETRA10
525!
526! ELM_OVER_COST :
527! 1 --> ismstr=10 or 12
528
529 SOLTELT(1) =.0000009528
530 SOLTELT(2) =.0000008540
531 SOLTELT(3) =.0000013420
532 SOLTELT(4) =.0000015180
533 SOLTELT(5) =.0000041200
534 SOLTELT(6) =.0000003812
535 SOLTELT(7) =.0000026600
536 SOLTELT(8) =.0000012750
537 SOLTELT(9) =.0000030520
538 SOLTELT(10) =.0000005769
539 TET4TELT(1) =.0000004744
540 TET4TELT(2) =.0000016460
541 TET10TELT(1)=.0000018340
542 NLOCAL_OPTION(1) = .0000002450
543
544 ELM_OVER_COST(1) = .0000001700
545 ELSE
546 SOLTELT(1) =.0000006901
547 SOLTELT(2) =.0000005650
548 SOLTELT(3) =.0000011240
549 SOLTELT(4) =.0000013490
550 SOLTELT(5) =.0000036580
551 SOLTELT(6) =.0000004100
552 SOLTELT(7) =.0000022490
553 SOLTELT(8) =.0000009907
554 SOLTELT(9) =.0000023800
555 SOLTELT(10) =.0000004480
556 TET4TELT(1) =.0000002688
557 TET4TELT(2) =.0000014960
558 TET10TELT(1)=.0000016200
559 NLOCAL_OPTION(1) = .0000002450
560
561 ELM_OVER_COST(1) = .0000001400
562 ENDIF
563
564! -------------
565! Coques
566! -------------
567
568
569
570! Tableau SHTNL
571! ---------------
572!
573! Temps par elt/cycle pour les materieux coques
574! SHTNL(I,J,K)=
575! 1er Champ : numero du materiaux
576! 2 champs : nbre de pts d integration, de taille 0:4
577! 4 sont consideres :
578! * 0 : 0pt d integ
579! * 1 : 1pt d integ
580! * 2 : 3pt d integ
581! * 3 : 5pt d integ
582! Les autres temps pour les NPT, sont calculees par interpolation avec ces 3
583! valeurs.
584!
585! 3e champ : option particuliere au besoin
586!
587! -------------
588! Triangles
589! -------------
590
591! Tableau TRITNL
592! ---------------
593!
594! Temps par elt/cycle pour les materieux coques
595! TRITNL(I,J,K)=
596! 1er Champ : numero du materiaux
597! 2e champ : nbre de pts d integration, de taille 0:4
598! 4 sont consideres :
599! * 0 : 0pt d integ
600! * 1 : 1pt d integ
601! * 2 : 3pt d integ
602! * 3 : 5pt d integ
603! Les autres temps pour les NPT, sont calculees par interpolation avec ces 3
604! valeurs.
605!
606! 3e champ : option particuliere au besoin
607!
608!
609! Tableau TRITELT
610! ---------------
611!
612! Temps par elt/cycle pour le traitement de l hourglas en fct de /IPARITH
613
614 DO I=1,3
615 DO J = 0, MAXLAW
616 SHTNL(J,0,I) = .0000001480
617 SHTNL(J,1,I) = .0000002013
618 SHTNL(J,2,I) = .0000003894
619 SHTNL(J,3,I) = .0000005561
620C
621 TRITNL(J,0,I)= .0000001373
622 TRITNL(J,1,I)= .0000002035
623 TRITNL(J,2,I)= .0000003892
624 TRITNL(J,3,I)= .0000005410
625 ENDDO
626 ENDDO
627! special option for shell or shell3n
628 SHTNL_OPTION(1:10,1:5) = 0.
629! LAW 01
630 SHTNL(1,0,1) = .0000000897
631 SHTNL(1,1,1) = .0000000811
632 SHTNL(1,2,1) = .0000001004
633 SHTNL(1,3,1) = .0000001032
634C
635 TRITNL(1,0,1)= .0000000999
636 TRITNL(1,1,1)= .0000000797
637 TRITNL(1,2,1)= .0000001024
638 TRITNL(1,3,1)= .0000000999
639! LAW 02
640 SHTNL(2,0,1) = .0000001480
641 SHTNL(2,1,1) = .0000002013
642 SHTNL(2,2,1) = .0000003894
643 SHTNL(2,3,1) = .0000005561
644C
645 TRITNL(2,0,1)= .0000001373
646 TRITNL(2,1,1)= .0000002035
647 TRITNL(2,2,1)= .0000003892
648 TRITNL(2,3,1)= .0000005410
649! LAW 02 CC
650 SHTNL(2,0,2) = .0000001606
651 SHTNL(2,1,2) = .0000002124
652 SHTNL(2,2,2) = .0000004265
653 SHTNL(2,3,2) = .0000006057
654C
655 TRITNL(2,0,2)= .0000001751
656 TRITNL(2,1,2)= .0000002033
657 TRITNL(2,2,2)= .0000004272
658 TRITNL(2,3,2)= .0000006010
659! LAW 15
660 SHTNL(15,0,1) = .0000003761
661 SHTNL(15,1,1) = .0000003761
662 SHTNL(15,2,1) = .0000008287
663 SHTNL(15,3,1) = .0000012890
664C
665 TRITNL(15,0,2)= .0000003652
666 TRITNL(15,1,2)= .0000003711
667 TRITNL(15,2,2)= .0000008201
668 TRITNL(15,3,2)= .0000011750
669
670! LAW 19
671 SHTNL(19,0,1) = .0000002260
672 SHTNL(19,1,1) = .0000002260
673 SHTNL(19,2,1) = .0000004591
674 SHTNL(19,3,1) = .0000006188
675C
676 TRITNL(19,0,1)= .0000002100
677 TRITNL(19,1,1)= .0000002100
678 TRITNL(19,2,1)= .0000004189
679 TRITNL(19,3,1)= .0000005672
680! LAW 22
681 SHTNL(22,0,1) = .0000002349
682 SHTNL(22,1,1) = .0000002947
683 SHTNL(22,2,1) = .0000005349
684 SHTNL(22,3,1) = .0000008162
685C
686 TRITNL(22,0,1)= .0000002209
687 TRITNL(22,1,1)= .0000003017
688 TRITNL(22,2,1)= .0000005351
689 TRITNL(22,3,1)= .0000008017
690! LAW 25
691 SHTNL(25,0,1) = .0000002511
692 SHTNL(25,1,1) = .0000002509
693 SHTNL(25,2,1) = .0000005094
694 SHTNL(25,3,1) = .0000008026
695C
696 TRITNL(25,0,1)= .0000002519
697 TRITNL(25,1,1)= .0000002455
698 TRITNL(25,2,1)= .0000005096
699 TRITNL(25,3,1)= .0000007655
700! LAW 25 Crasurv (/TYPE11)
701 SHTNL(25,0,2) = .0000004487
702 SHTNL(25,1,2) = .0000004487
703 SHTNL(25,2,2) = .0000009845
704 SHTNL(25,3,2) = .0000015870
705C
706 TRITNL(25,0,2)= .0000003866
707 TRITNL(25,1,2)= .0000003786
708 TRITNL(25,2,2)= .0000009420
709 TRITNL(25,3,2)= .0000013980
710! LAW 27 (0pt incompatible)
711 SHTNL(27,0,1) = .0000003490
712 SHTNL(27,1,1) = .0000003490
713 SHTNL(27,2,1) = .0000007537
714 SHTNL(27,3,1) = .0000011150
715C
716 TRITNL(27,0,1)= .0000003438
717 TRITNL(27,1,1)= .0000003396
718 TRITNL(27,2,1)= .0000007083
719 TRITNL(27,3,1)= .0000010850
720! LAW 32
721 SHTNL(32,0,1) = .0000001841
722 SHTNL(32,1,1) = .0000002502
723 SHTNL(32,2,1) = .0000004972
724 SHTNL(32,3,1) = .0000006134
725C
726 TRITNL(32,0,1)= .0000002125
727 TRITNL(32,1,1)= .0000002362
728 TRITNL(32,2,1)= .0000005149
729 TRITNL(32,3,1)= .0000007843
730! LAW 35 (0pt incompatible)
731 SHTNL(35,0,1) = .0000002696
732 SHTNL(35,1,1) = .0000002696
733 SHTNL(35,2,1) = .0000005399
734 SHTNL(35,3,1) = .0000008694
735C
736 TRITNL(35,0,1)= .0000002629
737 TRITNL(35,1,1)= .0000002629
738 TRITNL(35,2,1)= .0000005365
739 TRITNL(35,3,1)= .0000009575
740! LAW 36 - 1 fct
741 SHTNL(36,0,1) = .0000002453
742 SHTNL(36,1,1) = .0000002975
743 SHTNL(36,2,1) = .0000005370
744 SHTNL(36,3,1) = .0000007957
745C
746 TRITNL(36,0,1)= .0000002315
747 TRITNL(36,1,1)= .0000003005
748 TRITNL(36,2,1)= .0000005141
749 TRITNL(36,3,1)= .0000008021
750! LAW 36 - 5 fct
751 SHTNL(36,0,2) = .0000002832
752 SHTNL(36,1,2) = .0000003655
753 SHTNL(36,2,2) = .0000005902
754 SHTNL(36,3,2) = .0000008500
755C
756 TRITNL(36,0,2)= .0000002975
757 TRITNL(36,1,2)= .0000003608
758 TRITNL(36,2,2)= .0000005945
759 TRITNL(36,3,2)= .0000008260
760! LAW 36 - 10 fct
761 SHTNL(36,0,3) = .0000003180
762 SHTNL(36,1,3) = .0000003850
763 SHTNL(36,2,3) = .0000006500
764 SHTNL(36,3,3) = .0000009020
765C
766 TRITNL(36,0,3)= .0000003060
767 TRITNL(36,1,3)= .0000003820
768 TRITNL(36,2,3)= .0000006540
769 TRITNL(36,3,3)= .0000009000
770! LAW 42
771 SHTNL(42,0,1) = .0000005524
772 SHTNL(42,1,1) = .0000005524
773 SHTNL(42,2,1) = .0000013820
774 SHTNL(42,3,1) = .0000022760
775! LAW 42 M = 1
776 SHTNL(42,0,2) = .0000015100
777 SHTNL(42,1,2) = .0000015100
778 SHTNL(42,2,2) = .0000045280
779 SHTNL(42,3,2) = .0000074640
780! LAW 42 M = 2
781 SHTNL(42,0,3) = .0000017720
782 SHTNL(42,1,3) = .0000017720
783 SHTNL(42,2,3) = .0000053370
784 SHTNL(42,3,3) = .0000088600
785! LAW 42
786 TRITNL(42,0,1)= .0000005586
787 TRITNL(42,1,1)= .0000005586
788 TRITNL(42,2,1)= .0000015540
789 TRITNL(42,3,1)= .0000022710
790! LAW 42 M = 1
791 TRITNL(42,0,2)= .0000015000
792 TRITNL(42,1,2)= .0000015000
793 TRITNL(42,2,2)= .0000044580
794 TRITNL(42,3,2)= .0000073400
795! LAW 42 M = 2
796 TRITNL(42,0,3)= .0000017200
797 TRITNL(42,1,3)= .0000017200
798 TRITNL(42,2,3)= .0000052700
799 TRITNL(42,3,3)= .0000087000
800! LAW 43
801 SHTNL(43,0,1) = .0000002796
802 SHTNL(43,1,1) = .0000002918
803 SHTNL(43,2,1) = .0000006119
804 SHTNL(43,3,1) = .0000009136
805C
806 TRITNL(43,0,1)= .0000002856
807 TRITNL(43,1,1)= .0000002743
808 TRITNL(43,2,1)= .0000006279
809 TRITNL(43,3,1)= .0000008749
810! LAW 44
811 SHTNL(44,0,1) = .0000002640
812 SHTNL(44,1,1) = .0000002640
813 SHTNL(44,2,1) = .0000005714
814 SHTNL(44,3,1) = .0000008505
815C
816 TRITNL(44,0,1)= .0000002802
817 TRITNL(44,1,1)= .0000002802
818 TRITNL(44,2,1)= .0000005076
819 TRITNL(44,3,1)= .0000007598
820! LAW 48
821 SHTNL(48,0,1) = .0000002811
822 SHTNL(48,1,1) = .0000002811
823 SHTNL(48,2,1) = .0000005285
824 SHTNL(48,3,1) = .0000007925
825C
826 TRITNL(48,0,1)= .0000002764
827 TRITNL(48,1,1)= .0000002764
828 TRITNL(48,2,1)= .0000005447
829 TRITNL(48,3,1)= .0000007677
830! LAW 52
831 SHTNL(52,0,1) = .0000002313
832 SHTNL(52,1,1) = .0000002313
833 SHTNL(52,2,1) = .0000004775
834 SHTNL(52,3,1) = .0000007161
835C
836 TRITNL(52,0,1)= .0000002422
837 TRITNL(52,1,1)= .0000002422
838 TRITNL(52,2,1)= .0000004612
839 TRITNL(52,3,1)= .0000007237
840! LAW 57
841 SHTNL(57,0,1) = .0000005039
842 SHTNL(57,1,1) = .0000005039
843 SHTNL(57,2,1) = .0000013330
844 SHTNL(57,3,1) = .0000021130
845C
846 TRITNL(57,0,1)= .0000005568
847 TRITNL(57,1,1)= .0000005568
848 TRITNL(57,2,1)= .0000014330
849 TRITNL(57,3,1)= .0000022840
850! LAW 58
851! No visc (FLAGUL==0)
852 SHTNL(58,0,1) = .0000005435
853 SHTNL(58,1,1) = .0000006044
854 SHTNL(58,2,1) = .0000017910
855 SHTNL(58,3,1) = .0000028510
856C
857 TRITNL(58,0,1)= .0000005733
858 TRITNL(58,1,1)= .0000005733
859 TRITNL(58,2,1)= .0000015320
860 TRITNL(58,3,1)= .0000023640
861! LAW 60
862 SHTNL(60,0,1) = .0000003498
863 SHTNL(60,1,1) = .0000003992
864 SHTNL(60,2,1) = .0000008897
865 SHTNL(60,3,1) = .0000013350
866!
867 TRITNL(60,0,1)= .0000003850
868 TRITNL(60,1,1)= .0000003997
869 TRITNL(60,2,1)= .0000009103
870 TRITNL(60,3,1)= .0000014110
871! LAW 62
872 ! M = 0
873 SHTNL(62,0,1) = .0000008089
874 SHTNL(62,1,1) = .0000008089
875 SHTNL(62,2,1) = .0000021360
876 SHTNL(62,3,1) = .0000034950
877 ! M = 1
878 SHTNL(62,0,2) = .0000014640
879 SHTNL(62,1,2) = .0000014640
880 SHTNL(62,2,2) = .0000041200
881 SHTNL(62,3,2) = .0000067100
882 ! M = 2
883 SHTNL(62,0,3) = .0000017390
884 SHTNL(62,1,3) = .0000017390
885 SHTNL(62,2,3) = .0000048970
886 SHTNL(62,3,3) = .0000079110
887C
888 ! M = 0
889 TRITNL(62,0,1)= .0000008167
890 TRITNL(62,1,1)= .0000008167
891 TRITNL(62,2,1)= .0000021660
892 TRITNL(62,3,1)= .0000034760
893 ! M = 1
894 TRITNL(62,0,1)= .0000014950
895 TRITNL(62,1,1)= .0000014950
896 TRITNL(62,2,1)= .0000040800
897 TRITNL(62,3,1)= .0000067340
898 ! M = 2
899 TRITNL(62,0,1)= .0000017140
900 TRITNL(62,1,1)= .0000017140
901 TRITNL(62,2,1)= .0000048680
902 TRITNL(62,3,1)= .0000079480
903! LAW 63
904 SHTNL(63,0,1) = .0000002904
905 SHTNL(63,1,1) = .0000002904
906 SHTNL(63,2,1) = .0000006582
907 SHTNL(63,3,1) = .0000009983
908C
909 TRITNL(63,0,1)= .0000002766
910 TRITNL(63,1,1)= .0000002766
911 TRITNL(63,2,1)= .0000005871
912 TRITNL(63,3,1)= .0000008944
913! LAW 64
914 SHTNL(64,0,1) = .0000003364
915 SHTNL(64,1,1) = .0000003364
916 SHTNL(64,2,1) = .0000007866
917 SHTNL(64,3,1) = .0000011920
918C
919 TRITNL(64,0,1)= .0000003426
920 TRITNL(64,1,1)= .0000003426
921 TRITNL(64,2,1)= .0000007550
922 TRITNL(64,3,1)= .0000011970
923! LAW 65
924 SHTNL(65,0,1) = .0000003350
925 SHTNL(65,1,1) = .0000003350
926 SHTNL(65,2,1) = .0000007280
927 SHTNL(65,3,1) = .0000011240
928C
929 TRITNL(65,0,1)= .0000003324
930 TRITNL(65,1,1)= .0000003324
931 TRITNL(65,2,1)= .0000007364
932 TRITNL(65,3,1)= .0000011200
933! LAW 66
934 SHTNL(66,0,1) = .0000003310
935 SHTNL(66,1,1) = .0000003310
936 SHTNL(66,2,1) = .0000007200
937 SHTNL(66,3,1) = .0000011010
938C
939 TRITNL(66,0,1)= .0000003725
940 TRITNL(66,1,1)= .0000003666
941 TRITNL(66,2,1)= .0000007924
942 TRITNL(66,3,1)= .0000011650
943! LAW 69
944 ! M = 0 (no visc)
945 SHTNL(69,0,1) = .0000013330
946 SHTNL(69,1,1) = .0000013330
947 SHTNL(69,2,1) = .0000032970
948 SHTNL(69,3,1) = .0000054280
949 ! M = 1 (visc)
950 SHTNL(69,0,2) = .0000014210
951 SHTNL(69,1,2) = .0000014210
952 SHTNL(69,2,2) = .0000039620
953 SHTNL(69,3,2) = .0000064560
954 ! M = 2 (visc)
955 SHTNL(69,0,3) = .0000016580
956 SHTNL(69,1,3) = .0000016580
957 SHTNL(69,2,3) = .0000045270
958 SHTNL(69,3,3) = .0000074700
959C
960 ! M = 0 (no visc)
961 TRITNL(69,0,1)= .0000010400
962 TRITNL(69,1,1)= .0000010400
963 TRITNL(69,2,1)= .0000027280
964 TRITNL(69,3,1)= .0000044630
965 ! M = 1
966 TRITNL(69,0,2)= .0000012490
967 TRITNL(69,1,2)= .0000012490
968 TRITNL(69,2,2)= .0000035000
969 TRITNL(69,3,2)= .0000055580
970 ! M = 2
971 TRITNL(69,0,3)= .0000014430
972 TRITNL(69,1,3)= .0000014430
973 TRITNL(69,2,3)= .0000040510
974 TRITNL(69,3,3)= .0000065830
975! LAW 71 (incompatible npt=0)
976 SHTNL(71,0,1) = .0000005878
977 SHTNL(71,1,1) = .0000005878
978 SHTNL(71,2,1) = .0000014690
979 SHTNL(71,3,1) = .0000024250
980C
981 TRITNL(71,0,1)= .0000006319
982 TRITNL(71,1,1)= .0000006919
983 TRITNL(71,2,1)= .0000015420
984 TRITNL(71,3,1)= .0000024620
985! LAW 72
986 SHTNL(72,0,1) = .0000003671
987 SHTNL(72,1,1) = .0000003671
988 SHTNL(72,2,1) = .0000007838
989 SHTNL(72,3,1) = .0000012730
990C
991 TRITNL(72,0,1)= .0000003591
992 TRITNL(72,1,1)= .0000003591
993 TRITNL(72,2,1)= .0000007633
994 TRITNL(72,3,1)= .0000011460
995! LAW 73
996 SHTNL(73,0,1) = .0000003456
997 SHTNL(73,1,1) = .0000003689
998 SHTNL(73,2,1) = .0000007752
999 SHTNL(73,3,1) = .0000012010
1000C
1001 TRITNL(73,0,1)= .0000003618
1002 TRITNL(73,1,1)= .0000003371
1003 TRITNL(73,2,1)= .0000007637
1004 TRITNL(73,3,1)= .0000012300
1005! LAW 76
1006 SHTNL(76,0,1) = .0000007383
1007 SHTNL(76,1,1) = .0000007476
1008 SHTNL(76,2,1) = .0000018800
1009 SHTNL(76,3,1) = .0000029980
1010C
1011 TRITNL(76,0,1)= .0000007008
1012 TRITNL(76,1,1)= .0000007008
1013 TRITNL(76,2,1)= .0000018360
1014 TRITNL(76,3,1)= .0000029000
1015! LAW 78
1016 SHTNL(78,0,1) = .0000002921
1017 SHTNL(78,1,1) = .0000002921
1018 SHTNL(78,2,1) = .0000005980
1019 SHTNL(78,3,1) = .0000009682
1020C
1021 TRITNL(78,0,1)= .0000002936
1022 TRITNL(78,1,1)= .0000002939
1023 TRITNL(78,2,1)= .0000005670
1024 TRITNL(78,3,1)= .0000008430
1025! LAW 80
1026 SHTNL(80,0,1) = .0000015110
1027 SHTNL(80,1,1) = .0000014400
1028 SHTNL(80,2,1) = .0000040060
1029 SHTNL(80,3,1) = .0000065840
1030C
1031 TRITNL(80,0,1)= .0000013660
1032 TRITNL(80,1,1)= .0000013660
1033 TRITNL(80,2,1)= .0000039870
1034 TRITNL(80,3,1)= .0000063540
1035! LAW 82 (incompatible NPT=0)
1036! N = 1
1037 SHTNL(82,0,1) = .0000006991
1038 SHTNL(82,1,1) = .0000007022
1039 SHTNL(82,2,1) = .0000017930
1040 SHTNL(82,3,1) = .0000028300
1041!
1042 TRITNL(82,0,1)= .0000006886
1043 TRITNL(82,1,1)= .0000007115
1044 TRITNL(82,2,1)= .0000017680
1045 TRITNL(82,3,1)= .0000028740
1046! N = 2
1047 SHTNL(82,0,2) = .0000010470
1048 SHTNL(82,1,2) = .0000010430
1049 SHTNL(82,2,2) = .0000027830
1050 SHTNL(82,3,2) = .0000045710
1051!
1052 TRITNL(82,0,2)= .0000010690
1053 TRITNL(82,1,2)= .0000010330
1054 TRITNL(82,2,2)= .0000028430
1055 TRITNL(82,3,2)= .0000046350
1056! N = 3
1057 SHTNL(82,0,3) = .0000013880
1058 SHTNL(82,1,3) = .0000014040
1059 SHTNL(82,2,3) = .0000038590
1060 SHTNL(82,3,3) = .0000063230
1061!
1062 TRITNL(82,0,3)= .0000014340
1063 TRITNL(82,1,3)= .0000014220
1064 TRITNL(82,2,3)= .0000038320
1065 TRITNL(82,3,3)= .0000063450
1066! LAW 87
1067 SHTNL(87,0,1) = .0000003592
1068 SHTNL(87,1,1) = .0000003592
1069 SHTNL(87,2,1) = .0000007907
1070 SHTNL(87,3,1) = .0000011520
1071!
1072 TRITNL(87,0,1)= .0000003754
1073 TRITNL(87,1,1)= .0000003874
1074 TRITNL(87,2,1)= .0000008403
1075 TRITNL(87,3,1)= .0000012850
1076! LAW 93
1077 SHTNL(93,0,1) = .0000002801
1078 SHTNL(93,1,1) = .0000002801
1079 SHTNL(93,2,1) = .0000005238
1080 SHTNL(93,3,1) = .0000007912
1081!
1082 TRITNL(93,0,1)= .0000002629
1083 TRITNL(93,1,1)= .0000002629
1084 TRITNL(93,2,1)= .0000005262
1085 TRITNL(93,3,1)= .0000008110
1086! LAW 104
1087 ! Nice algo
1088 SHTNL(104,0,1) = .0000002823
1089 SHTNL(104,1,1) = .0000002823
1090 SHTNL(104,2,1) = .0000005580
1091 SHTNL(104,3,1) = .0000008341
1092 ! Newton algo
1093 SHTNL(104,0,2) = .0000004489
1094 SHTNL(104,1,2) = .0000004489
1095 SHTNL(104,2,2) = .0000008839
1096 SHTNL(104,3,2) = .0000011360
1097 ! Nice algo
1098 TRITNL(104,0,1) = .0000002823
1099 TRITNL(104,1,1) = .0000002823
1100 TRITNL(104,2,1) = .0000005580
1101 TRITNL(104,3,1) = .0000008341
1102 ! Newton algo
1103 TRITNL(104,0,2) = .0000004489
1104 TRITNL(104,1,2) = .0000004489
1105 TRITNL(104,2,2) = .0000008839
1106 TRITNL(104,3,2) = .0000011360
1107 ! iloc=1
1108 SHTNL_OPTION(1,1) = .0000000900
1109 ! iloc=2
1110 SHTNL_OPTION(1,2) = .0000001000
1111 ! iloc=3
1112 SHTNL_OPTION(1,3) = .0000002300
1113C Temps de traitement de la partie elementaire sans la loi materielle fonction de IPARITH
1114 IF (IPARI0.EQ.1) THEN
1115! SHELT :
1116! 1-->5 --> /SHELL
1117! 6-->10 --> /SH_COMP or TYPE10
1118! 11-->15 --> /SH_SANDW or TYPE11
1119! 15-->20 --> /SH_ORTH or TYPE9
1120! 20-->25 --> /SH_FABR or TYPE16
1121! IDRILL = 0
1122! 1/6/11/16/21 --> Q4
1123! 2/7/12/17/22 --> QEPH
1124! 3/8/13/18/23 --> BATOZ
1125! IDRILL = 1
1126! none --> Q4
1127! 4/9 /14/19/24 --> QEPH
1128! 5/10/15/20/25 --> BATOZ
1129! TRIELT :
1130! 1 --> tri
1131! NLOCAL_OPTION
1132! 1 --> sol
1133! 2 --> shell
1134! /SHELL
1135 SHTELT(1) = .0000005629
1136 SHTELT(2) = .0000007050
1137 SHTELT(3) = .0000013190
1138 SHTELT(4) = .0000009213
1139 SHTELT(5) = .0000016710
1140! /SH_COMP
1141 SHTELT(6) = .0000007170
1142 SHTELT(7) = .0000009040
1143 SHTELT(8) = .0000015630
1144 SHTELT(9) = .0000011660
1145 SHTELT(10)= .0000018400
1146! /SH_SANDW
1147 SHTELT(11) = .0000006329
1148 SHTELT(12) = .0000008915
1149 SHTELT(13) = .0000014100
1150 SHTELT(14) = .0000012220
1151 SHTELT(15) = .0000018450
1152! /SH_ORTH
1153 SHTELT(16) = .0000004117
1154 SHTELT(17) = .0000007741
1155 SHTELT(18) = .0000014700
1156 SHTELT(19) = .0000010510
1157 SHTELT(20) = .0000016530
1158! /SH_FABR
1159 SHTELT(21) = .0000006739
1160 SHTELT(22) = .0000007792
1161 SHTELT(23) = .0000014680
1162 SHTELT(24) = .0000011780
1163 SHTELT(25) = .0000018450
1164
1165! /TYPE51
1166 SHTELT(26) = .0000014800
1167 SHTELT(27) = .0000017300
1168 SHTELT(28) = .0000035000
1169 SHTELT(29) = .0000017700
1170 SHTELT(30) = .0000051700
1171
1172 TRITELT(1)= .0000004724
1173
1174 NLOCAL_OPTION(2) = .0000001000
1175 ELSE
1176! /SHELL
1177 SHTELT(1) = .0000004647
1178 SHTELT(2) = .0000005736
1179 SHTELT(3) = .0000012330
1180 SHTELT(4) = .0000007605
1181 SHTELT(5) = .0000013700
1182! /SH_COMP
1183 SHTELT(6) = .0000005935
1184 SHTELT(7) = .0000007462
1185 SHTELT(8) = .0000013800
1186 SHTELT(9) = .0000009625
1187 SHTELT(10)= .0000015200
1188! /SH_SANDW
1189 SHTELT(11) = .0000005224
1190 SHTELT(12) = .0000007354
1191 SHTELT(13) = .0000012600
1192 SHTELT(14) = .0000010100
1193 SHTELT(15) = .0000015200
1194! /SH_ORTH
1195 SHTELT(16) = .0000003400
1196 SHTELT(17) = .0000006300
1197 SHTELT(18) = .0000013040
1198 SHTELT(19) = .0000008676
1199 SHTELT(20) = .0000013600
1200! /SH_FABR
1201 SHTELT(21) = .0000005500
1202 SHTELT(22) = .0000006440
1203 SHTELT(23) = .0000012700
1204 SHTELT(24) = .0000009658
1205 SHTELT(25) = .0000015100
1206
1207! /TYPE51
1208 SHTELT(26) = .0000012300
1209 SHTELT(27) = .0000014400
1210 SHTELT(28) = .0000029000
1211 SHTELT(29) = .0000014700
1212 SHTELT(30) = .0000038500
1213
1214 TRITELT(1)= .0000003936
1215
1216 NLOCAL_OPTION(2) = .0000001000
1217
1218 ENDIF
1219
1220! il y a un facteur multiplicatif sur les coques de Batoz, car on rentre 4x dans cmain3.
1221
1222 BATOZMULT = 3.8
1223! /PRONY or /LPRONY overcost
1224 VISC_PRONY_COST = .0000000243
1225
1226! *******************************************************************
1227! Failure model
1228! *******************************************************************
1229! Johnson : 1
1230! Tbutcher : 2
1231! Wilkins : 3
1232! user1 : 4
1233! user2 : 5
1234! user3 : 6
1235! FLD : 7 (for shell only)
1236! Spalling : 8 (for solid only)
1237! Wierzbicki : 9
1238! Tensstrain : 10
1239! Energy : 11
1240! ??? : 12
1241! Chang : 13 (for shell only)
1242! Hashin : 14
1243! Yamada : 15
1244! Puck : 16
1245! XFEM_Tbutcher : 17
1246! Lad_dama : 18 (for solid only)
1247! Power_da : 19
1248! Connect : 20 (for solid only)
1249! XFEM_Johns : 21
1250! XFEM_FLD : 22
1251! TAB1 : 23
1252! Orthstra : 24 (for shell only)
1253! NXT : 25 (for shell only)
1254! Snconnect : 26 (for solid only)
1255! EMC : 27 (for solid only)
1256! Wind : 28 (for shell only)
1257! Sahraei1 : 29 (for solid only)
1258! Biquadratic : 30
1259! _
1260! SHELL : 1--> 0 ou 1pt |
1261! 2--> 3 pts |--- MLN < 28 or ==32
1262! 3--> 5 pts _|
1263! 4--> 0 ou 1pt |
1264! 5--> 3 pts |--- MLN >= 28 or /=32
1265! 6--> 5 pts _|
1266!
1267! SOLIDE : 1 ----> MLN <28 or ==49
1268! 2 ----> MLN >=28 or /=49
1269! --------------------------------------------------------
1270 DO I=1,MAXFAIL
1271 DO J=1,3
1272 RUPTURE_SOL(I,J) = 0.
1273 RUPTURE_SOL12(I,J) = 0.
1274 RUPTURE_TET4(I,J) = 0.
1275 RUPTURE_TET10(I,J) = 0.
1276 ENDDO
1277 DO J=1,6
1278 RUPTURE_SHELL(I,J) = 0.
1279 ENDDO
1280 ENDDO
1281! --------------------
1282! Johnson
1283 ! MLN <28 or ==32
1284 RUPTURE_SHELL(1,1) = .00000007517
1285 RUPTURE_SHELL(1,2) = .00000019090
1286 RUPTURE_SHELL(1,3) = .00000031460
1287
1288 ! MLN >=28 or /=32
1289 RUPTURE_SHELL(1,4) = .00000003021
1290 RUPTURE_SHELL(1,5) = .00000007567
1291 RUPTURE_SHELL(1,6) = .00000011620
1292
1293 ! MLN <28 or ==49
1294 RUPTURE_SOL(1,1) = .00000008839
1295 RUPTURE_SOL12(1,1) = .00000017710
1296 RUPTURE_TET4(1,1) = .00000008021
1297 RUPTURE_TET10(1,1) = .00000033420
1298 ! MLN >=28 or /=49
1299 RUPTURE_SOL(1,2) = .00000002696
1300 RUPTURE_SOL12(1,2) = .00000015790
1301 RUPTURE_TET4(1,2) = .00000002663
1302 RUPTURE_TET10(1,2) = .00000010360
1303! --------------------
1304! Tbutcher
1305 ! MLN <28 or ==49
1306 RUPTURE_SHELL(2,1) = .00000004873
1307 RUPTURE_SHELL(2,2) = .00000010840
1308 RUPTURE_SHELL(2,3) = .00000016340
1309
1310 RUPTURE_SHELL(2,4) = .00000003994
1311 RUPTURE_SHELL(2,5) = .00000010130
1312 RUPTURE_SHELL(2,6) = .00000015690
1313
1314 ! MLN <28 or ==49
1315 RUPTURE_SOL(2,1) = .00000004718
1316 RUPTURE_SOL12(2,1) = .00000031390
1317 RUPTURE_TET4(2,1) = .00000010090
1318 RUPTURE_TET10(2,1) = .00000040130
1319
1320 ! MLN >=28 or /=49
1321 RUPTURE_SOL(2,2) = .00000004949
1322 RUPTURE_SOL12(2,2) = .00000030290
1323 RUPTURE_TET4(2,2) = .00000004449
1324 RUPTURE_TET10(2,2) = .00000016800
1325! --------------------
1326! Wilkins
1327 ! MLN <28 or ==49
1328 RUPTURE_SHELL(3,1) = .00000011660
1329 RUPTURE_SHELL(3,2) = .00000032760
1330 RUPTURE_SHELL(3,3) = .00000052680
1331
1332 RUPTURE_SHELL(3,4) = .00000011410
1333 RUPTURE_SHELL(3,5) = .00000032760
1334 RUPTURE_SHELL(3,6) = .00000052910
1335
1336 ! MLN <28 or ==49
1337 RUPTURE_SOL(3,1) = .00000023210
1338 RUPTURE_SOL12(3,1) = .00000131700
1339 RUPTURE_TET4(3,1) = .00000016440
1340 RUPTURE_TET10(3,1) = .00000075060
1341
1342 ! MLN >=28 or /=49
1343 RUPTURE_SOL(3,2) = .00000016870
1344 RUPTURE_SOL12(3,2) = .00000129800
1345 RUPTURE_TET4(3,2) = .00000011060
1346 RUPTURE_TET10(3,2) = .00000051000
1347! --------------------
1348! FLD (only for shell)
1349 ! MLN <28 or ==32
1350 RUPTURE_SHELL(7,1) = .00000010900
1351 RUPTURE_SHELL(7,2) = .00000021590
1352 RUPTURE_SHELL(7,3) = .00000031890
1353
1354 ! MLN >=28 or /=32
1355 RUPTURE_SHELL(7,4) = .00000007162
1356 RUPTURE_SHELL(7,5) = .00000019550
1357 RUPTURE_SHELL(7,6) = .00000031220
1358! --------------------
1359! Spalling (only for solid)
1360
1361 ! MLN <28 or ==49
1362 RUPTURE_SOL(8,1) = .00000009428
1363 RUPTURE_SOL12(8,1) = .00000025070
1364 RUPTURE_TET4(8,1) = .00000009402
1365 RUPTURE_TET10(8,1) = .00000038110
1366
1367 ! MLN >=28 or /=49
1368 RUPTURE_SOL(8,2) = .00000003871
1369 RUPTURE_SOL12(8,2) = .00000023730
1370 RUPTURE_TET4(8,2) = .00000003716
1371 RUPTURE_TET10(8,2) = .00000014820
1372! --------------------
1373! Wierzbicki
1374 ! MLN <28 or ==49
1375 RUPTURE_SHELL(9,1) = .00000025370
1376 RUPTURE_SHELL(9,2) = .00000071910
1377 RUPTURE_SHELL(9,3) = .00000118100
1378
1379 RUPTURE_SHELL(9,4) = .00000024610
1380 RUPTURE_SHELL(9,5) = .00000071940
1381 RUPTURE_SHELL(9,6) = .00000118500
1382
1383 ! MLN <28 or ==49
1384 RUPTURE_SOL(9,1) = .00000008443
1385 RUPTURE_SOL12(9,1) = .00000011190
1386 RUPTURE_TET4(9,1) = .00000008236
1387 RUPTURE_TET10(9,1) = .00000032530
1388
1389 ! MLN >=28 or /=49
1390 RUPTURE_SOL(9,2) = .00000002591
1391 RUPTURE_SOL12(9,2) = .00000016310
1392 RUPTURE_TET4(9,2) = .00000002483
1393 RUPTURE_TET10(9,2) = .00000008821
1394! --------------------
1395! Tensstrain
1396 ! MLN <28 or ==32
1397 RUPTURE_SHELL(10,1) = .00000009068
1398 RUPTURE_SHELL(10,2) = .00000016430
1399 RUPTURE_SHELL(10,3) = .00000023120
1400
1401 ! MLN >=28 or /=32
1402 RUPTURE_SHELL(10,4) = .00000004473
1403 RUPTURE_SHELL(10,5) = .00000012650
1404 RUPTURE_SHELL(10,6) = .00000019170
1405
1406 ! MLN <28 or ==49
1407 RUPTURE_SOL(10,1) = .00000013990
1408 RUPTURE_SOL12(10,1) = .00000025600
1409 RUPTURE_TET4(10,1) = .00000010960
1410 RUPTURE_TET10(10,1) = .00000047680
1411
1412 ! MLN >=28 or /=49
1413 RUPTURE_SOL(10,2) = .00000008052
1414 RUPTURE_SOL12(10,2) = .00000056320
1415 RUPTURE_TET4(10,2) = .00000005430
1416 RUPTURE_TET10(10,2) = .00000025250
1417! --------------------
1418! Energy
1419 ! MLN <28 or ==49
1420 RUPTURE_SHELL(11,1) = .00000005295
1421 RUPTURE_SHELL(11,2) = .00000011720
1422 RUPTURE_SHELL(11,3) = .00000017760
1423
1424 RUPTURE_SHELL(11,4) = .00000003770
1425 RUPTURE_SHELL(11,5) = .00000010100
1426 RUPTURE_SHELL(11,6) = .00000015840
1427
1428 ! MLN <28 or ==49
1429 RUPTURE_SOL(11,1) = .00000009413
1430 RUPTURE_SOL12(11,1) = .00000009015
1431 RUPTURE_TET4(11,1) = .00000009164
1432 RUPTURE_TET10(11,1) = .00000037960
1433
1434 ! MLN >=28 or /=49
1435 RUPTURE_SOL(11,2) = .00000003883
1436 RUPTURE_SOL12(11,2) = .00000025910
1437 RUPTURE_TET4(11,2) = .00000003753
1438 RUPTURE_TET10(11,2) = .00000014030
1439! --------------------
1440! Chang (only for shell)
1441 ! MLN <28 or ==49
1442 RUPTURE_SHELL(13,1) = .00000004822
1443 RUPTURE_SHELL(13,2) = .00000010740
1444 RUPTURE_SHELL(13,3) = .00000016000
1445 ! MLN >=28 or /=49
1446 RUPTURE_SHELL(13,4) = .00000003770
1447 RUPTURE_SHELL(13,5) = .00000010100
1448 RUPTURE_SHELL(13,6) = .00000015840
1449! --------------------
1450! Hashin
1451 ! MLN <28 or ==49
1452 RUPTURE_SHELL(14,1) = .00000007989
1453 RUPTURE_SHELL(14,2) = .00000025240
1454 RUPTURE_SHELL(14,3) = .00000040300
1455
1456 RUPTURE_SHELL(14,4) = .00000007989
1457 RUPTURE_SHELL(14,5) = .00000025240
1458 RUPTURE_SHELL(14,6) = .00000040300
1459
1460 ! MLN <28 or ==49
1461 RUPTURE_SOL(14,1) = .00000013480
1462 RUPTURE_SOL12(14,1) = .00000107400
1463 RUPTURE_TET4(14,1) = .00000012500
1464 RUPTURE_TET10(14,1) = .00000054770
1465
1466 ! MLN >=28 or /=49
1467 RUPTURE_SOL(14,2) = .00000001289
1468 RUPTURE_SOL12(14,2) = .00000010990
1469 RUPTURE_TET4(14,2) = .00000001302
1470 RUPTURE_TET10(14,2) = .00000005823
1471! --------------------
1472! Puck (compatible MLN<28)
1473 ! MLN <28 or ==49
1474 RUPTURE_SHELL(16,1) = .00000007462
1475 RUPTURE_SHELL(16,2) = .00000018470
1476 RUPTURE_SHELL(16,3) = .00000028840
1477
1478 RUPTURE_SHELL(16,4) = .00000006606
1479 RUPTURE_SHELL(16,5) = .00000017190
1480 RUPTURE_SHELL(16,6) = .00000027600
1481
1482 ! MLN <28 or ==49
1483 RUPTURE_SOL(16,1) = .00000013060
1484 RUPTURE_SOL12(16,1) = .00000008806
1485 RUPTURE_TET4(16,1) = .00000012170
1486 RUPTURE_TET10(16,1) = .00000042400
1487
1488 RUPTURE_SOL(16,2) = .00000005419
1489 RUPTURE_SOL12(16,2) = .00000007646
1490 RUPTURE_TET4(16,2) = .00000001338
1491 RUPTURE_TET10(16,2) = .00000004460
1492! --------------------
1493! Lad_Dama (compatible MLN<28) (only for solid)
1494 ! MLN <28 or ==49
1495 RUPTURE_SOL(18,1) = .00000016890
1496 RUPTURE_SOL12(18,1) = .00000008820
1497 RUPTURE_TET4(18,1) = .00000013980
1498 RUPTURE_TET10(18,1) = .00000058630
1499
1500 RUPTURE_SOL(18,2) = .00000016890
1501 RUPTURE_SOL12(18,2) = .00000008820
1502 RUPTURE_TET4(18,2) = .00000013980
1503 RUPTURE_TET10(18,2) = .00000058630
1504! --------------------
1505! Connect (compatible MLN==59 et 83)
1506 ! MLN==59
1507 RUPTURE_SOL(20,1) = .00000054140
1508 RUPTURE_SOL12(20,1) = .00000054140
1509 RUPTURE_TET4(20,1) = .00000054140
1510 RUPTURE_TET10(20,1) = .00000054140
1511 ! MLN==83
1512 RUPTURE_SOL(20,2) = .00000054140
1513 RUPTURE_SOL12(20,2) = .00000054140
1514 RUPTURE_TET4(20,2) = .00000054140
1515 RUPTURE_TET10(20,2) = .00000054140
1516! --------------------
1517! Tab1
1518 ! MLN <28 or ==49
1519 RUPTURE_SHELL(23,1) = .00000021020
1520 RUPTURE_SHELL(23,2) = .00000060850
1521 RUPTURE_SHELL(23,3) = .00000098850
1522
1523 RUPTURE_SHELL(23,4) = .00000021590
1524 RUPTURE_SHELL(23,5) = .00000062580
1525 RUPTURE_SHELL(23,6) = .0000009949
1526
1527 ! MLN <28 or ==49
1528 RUPTURE_SOL(23,1) = .00000031120
1529 RUPTURE_SOL12(23,1) = .00000192900
1530 RUPTURE_TET4(23,1) = .00000026620
1531 RUPTURE_TET10(23,1) = .00000118100
1532
1533 ! MLN >=28 or /=49
1534 RUPTURE_SOL(23,2) = .00000025940
1535 RUPTURE_SOL12(23,2) = .00000198600
1536 RUPTURE_TET4(23,2) = .00000022330
1537 RUPTURE_TET10(23,2) = .00000096400
1538! --------------------
1539! NXT (only for shell)
1540 ! MLN <28 or ==32
1541 RUPTURE_SHELL(25,1) = .00000014880
1542 RUPTURE_SHELL(25,2) = .00000032860
1543 RUPTURE_SHELL(25,3) = .00000051100
1544
1545 ! MLN >=28 or /=32
1546 RUPTURE_SHELL(25,4) = .00000008550
1547 RUPTURE_SHELL(25,5) = .00000024030
1548 RUPTURE_SHELL(25,6) = .00000038060
1549! --------------------
1550! Snconnect (compatible MLN==59 et 83)
1551 ! MLN==59
1552 RUPTURE_SOL(26,1) = .00000065150
1553 RUPTURE_SOL12(26,1) = .00000065150
1554 RUPTURE_TET4(26,1) = .00000065150
1555 RUPTURE_TET10(26,1) = .00000065150
1556 ! MLN==83
1557 RUPTURE_SOL(26,2) = .00000065150
1558 RUPTURE_SOL12(26,2) = .00000065150
1559 RUPTURE_TET4(26,2) = .00000065150
1560 RUPTURE_TET10(26,2) = .00000065150
1561! --------------------
1562! EMC (only for solid)
1563 ! MLN <28 or ==49
1564 RUPTURE_SOL(27,1) = .00000007774
1565 RUPTURE_SOL12(27,1) = .00000013030
1566 RUPTURE_TET4(27,1) = .00000007539
1567 RUPTURE_TET10(27,1) = .00000029830
1568
1569 ! MLN >=28 or /=49
1570 RUPTURE_SOL(27,2) = .00000001913
1571 RUPTURE_SOL12(27,2) = .00000010660
1572 RUPTURE_TET4(27,2) = .00000001654
1573 RUPTURE_TET10(27,2) = .00000006008
1574! --------------------
1575! Alter (only shell/tri)
1576 ! MLN <28 or ==49
1577 RUPTURE_SHELL(28,1) = .0000001580
1578 RUPTURE_SHELL(28,2) = .0000003324
1579 RUPTURE_SHELL(28,3) = .0000006238
1580
1581 RUPTURE_SHELL(28,4) = .0000001580
1582 RUPTURE_SHELL(28,5) = .0000003324
1583 RUPTURE_SHELL(28,6) = .0000006238
1584! --------------------
1585! Biquad
1586 ! MLN <28 or ==49
1587 RUPTURE_SHELL(30,1) = .00000003069
1588 RUPTURE_SHELL(30,2) = .00000005533
1589 RUPTURE_SHELL(30,3) = .00000008252
1590
1591 RUPTURE_SHELL(30,4) = .00000003133
1592 RUPTURE_SHELL(30,5) = .00000006587
1593 RUPTURE_SHELL(30,6) = .00000010300
1594
1595 ! MLN <28 or ==49
1596 RUPTURE_SOL(30,1) = .00000008026
1597 RUPTURE_SOL12(30,1) = .00000013200
1598 RUPTURE_TET4(30,1) = .00000007523
1599 RUPTURE_TET10(30,1) = .00000030150
1600
1601 RUPTURE_SOL(30,2) = .00000001929
1602 RUPTURE_SOL12(30,2) = .00000011930
1603 RUPTURE_TET4(30,2) = .00000001793
1604 RUPTURE_TET10(30,2) = .00000006588
1605! --------------------
cmain3
subroutine cmain3(timers, elbuf_str, jft, jlt, nft, iparg, nel, mtn, ipla, ithk, group_param, pm, geo, npf, tf, bufmat, ssp, rho, viscmx, dt1c, sigy, area, exx, eyy, exy, exz, eyz, kxx, kyy, kxy, nu, off, thk0, mat, pid, mat_elem, for, mom, gstr, failwave, fwave_el, thk, eint, iofc, g, a11, a12, vol0, indxdel, ngl, zcfac, shf, gs, epsd_pg, kfts, jhbe, alpe, dir_a, dir_b, igeo, ipm, ifailure, npg, fheat, tempel, die, jthe, iexpan, tempel0, ishplyxfem, ply_exx, ply_eyy, ply_exy, ply_exz, ply_eyz, ply_f, del_ply, th_iply, sig_iply, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, ng, table, ixfem, offi, sensors, a11_iply, elcrkini, dir1_crk, dir2_crk, aldt, idt_therm, theaccfact, ismstr, ir, is, nlay, npt, ixlay, ixel, isubstack, stack, f_def, itask, drape, varnl, nloc_dmg, indx_drape, thke, sedrape, numel_drape, dt, ncycle, snpc, stf, nxlaymax, idel7nok, userl_avail, maxfunc, varnl_npttot, sbufmat, sdir_a, sdir_b, for_g)
Definition cmain3.F:87
the
end diagonal values have been computed in the(sparse) matrix id.SOL
N
#define N
Definition example_user.c:10
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
mumps_ooc_common::fct
integer, parameter fct
Definition mumps_ooc_common.F:16
y
y
Definition schur_example.m:54
sol
sol
Definition schur_example.m:54
x
x
Definition schur_example.m:53