OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
weights_p4linux964_spmd.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 Broadwell machine (AVX2)
25! Used machine : Presto --> Intel(R) Xeon(R) CPU E7-8890 v4 @ 2.20GHz
26
27
28! temps de reference correspond au temps /elt/materiaux loi2 0 pt integ
29 TPSREF = .0000004784 !.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) = .0000001235 !.0000001512
64 SOL8TNL(I,J) = .0000004478 !.0000006297
65 TET4TNL(I,J) = .0000001078 !.0000001542
66 TET10TNL(I,J)= .0000004025 !.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) = .0000001139
76 SOL8TNL(1,1) = .0000002589
77 TET4TNL(1,1) = .0000000916
78 TET10TNL(1,1)= .0000003281
79! LAW 01 + ISMSTR = 10 or 12
80 SOL1TNL(1,2) = .0000002010
81 SOL8TNL(1,2) = .0000016080
82 TET4TNL(1,2) = .0000002010
83 TET10TNL(1,2)= .0000008040
84! LAW 02
85 SOL1TNL(2,1) = .0000001235 !.0000001512
86 SOL8TNL(2,1) = .0000004478 !.0000006297
87 TET4TNL(2,1) = .0000001078 !.0000001542
88 TET10TNL(2,1)= .0000004025 !.0000005545
89! LAW 02 ICC=1
90 SOL1TNL(2,2) = .0000001811 !.0000002051
91 SOL8TNL(2,2) = .0000006586 !.0000010080
92 TET4TNL(2,2) = .0000001496 !.0000001820
93 TET10TNL(2,2)= .0000005725 !.0000007670
94! LAW 03
95 SOL1TNL(3,1) = .0000001620 !.0000002145
96 SOL8TNL(3,1) = .0000004681 !.0000006086
97 TET4TNL(3,1) = .0000001564 !.0000002292
98 TET10TNL(3,1)= .0000006055 !.0000009112
99! LAW 04
100 SOL1TNL(4,1) = .0000001800 !.0000001555
101 SOL8TNL(4,1) = .0000012600 !.0000006050
102 TET4TNL(4,1) = .0000001727 !.0000002547
103 TET10TNL(4,1)= .0000006580 !.0000009900
104! LAW 06
105 SOL1TNL(6,1) = .0000001680 !.0000002483
106 SOL8TNL(6,1) = .0000011900 !.0000006004
107 TET4TNL(6,1) = .0000001600 !.0000002412
108 TET10TNL(6,1)= .0000004800 !.0000007858
109! LAW 10
110 SOL1TNL(10,1) = .0000001807 ! .0000002024
111 SOL8TNL(10,1) = .0000011540 ! .0000013650
112 TET4TNL(10,1) = .0000001395 ! .0000001815
113 TET10TNL(10,1)= .0000005183 ! .0000006534
114! LAW 12
115 SOL1TNL(12,1) = .0000003025 ! .0000003941
116 SOL8TNL(12,1) = .0000022960 ! .0000038700
117 TET4TNL(12,1) = .0000002812 ! .0000003670
118 TET10TNL(12,1)= .0000011590 ! .0000013560
119! LAW 21
120 SOL1TNL(21,1) = .0000001211 ! .0000001872
121 SOL8TNL(21,1) = .0000011280 ! .0000013840
122 TET4TNL(21,1) = .0000001338 ! .0000001811
123 TET10TNL(21,1)= .0000005061 ! .0000006885
124! LAW 22
125 SOL1TNL(22,1) = .0000001341 !.0000001854
126 SOL8TNL(22,1) = .0000011850 !.0000013990
127 TET4TNL(22,1) = .0000001409 !.0000001663
128 TET10TNL(22,1)= .0000005420 !.0000007349
129! LAW 23
130 SOL1TNL(23,1) = .0000001537 !.0000001890
131 SOL8TNL(23,1) = .0000012070 !.0000014360
132 TET4TNL(23,1) = .0000001438 !.0000001850
133 TET10TNL(23,1)= .0000005498 !.0000007381
134! LAW 24
135 SOL1TNL(24,1) = .0000003916 !.0000006082
136 SOL8TNL(24,1) = .0000029130 !.0000042810
137 TET4TNL(24,1) = .0000003826 !.0000004260
138 TET10TNL(24,1)= .0000014510 !.0000020050
139! LAW 25
140 SOL1TNL(25,1) = .0000003084
141 SOL8TNL(25,1) = .0000024600
142 TET4TNL(25,1) = .0000003084
143 TET10TNL(25,1)= .0000012330
144! LAW 25 Crasurv
145 SOL1TNL(25,2) = .0000003634
146 SOL8TNL(25,2) = .0000029070
147 TET4TNL(25,2) = .0000003634
148 TET10TNL(25,2)= .0000014570
149! LAW 28
150 SOL1TNL(28,1) = .0000002758 !.0000003664
151 SOL8TNL(28,1) = .0000016610 !.0000021385
152 TET4TNL(28,1) = .0000002680 !.0000003168
153 TET10TNL(28,1)= .0000011490 !.0000013620
154! LAW 33 0 fct
155 SOL1TNL(33,1) = .0000002991 !.0000004453
156 SOL8TNL(33,1) = .0000019190 !.0000021870
157 TET4TNL(33,1) = .0000002858 !.0000004027
158 TET10TNL(33,1)= .0000011570 !.0000015750
159! LAW 33 1 fct
160 SOL1TNL(33,2) = .0000004217 !.0000006103
161 SOL8TNL(33,2) = .0000027120 !.0000025320
162 TET4TNL(33,2) = .0000004088 !.0000005659
163 TET10TNL(33,2)= .0000015910 !.0000018270
164! LAW 34
165 SOL1TNL(34,1) = .0000001847 !.0000002516
166 SOL8TNL(34,1) = .0000008083 !.0000012240
167 TET4TNL(34,1) = .0000001696 !.0000002258
168 TET10TNL(34,1)= .0000006402 !.0000008617
169! LAW 35
170 SOL1TNL(35,1) = .0000001947 !.0000003070
171 SOL8TNL(35,1) = .0000010040 !.0000017020
172 TET4TNL(35,1) = .0000002021 !.0000002691
173 TET10TNL(35,1)= .0000007814 !.0000010870
174! LAW 36 1 fonction
175 SOL1TNL(36,1) = .0000001962 !.0000003096
176 SOL8TNL(36,1) = .0000011260 !.0000016040
177 TET4TNL(36,1) = .0000001950 !.0000002877
178 TET10TNL(36,1)= .0000008260 !.0000009584
179! LAW 36 5 fonctions
180 SOL1TNL(36,2) = .0000002010 !.0000003358
181 SOL8TNL(36,2) = .0000011400 !.0000016420
182 TET4TNL(36,2) = .0000002000 !.0000002724
183 TET10TNL(36,2)= .0000008650 !.0000010120
184! LAW 36 10 fonctions
185 SOL1TNL(36,3) = .0000002100 !.0000003702
186 SOL8TNL(36,3) = .0000011800 !.0000016950
187 TET4TNL(36,3) = .0000002050 !.0000002854
188 TET10TNL(36,3)= .0000008800 !.0000010970
189! LAW 38
190 SOL1TNL(38,1) = .0000007220
191 SOL8TNL(38,1) = .0000053000
192 TET4TNL(38,1) = .0000007220
193 TET10TNL(38,1)= .0000028910
194! LAW 38 + ISMSTR=10 or 12
195 SOL1TNL(38,2) = .0000007815
196 SOL8TNL(38,2) = .0000061000
197 TET4TNL(38,2) = .0000007815
198 TET10TNL(38,2)= .0000030900
199! LAW 40
200 SOL1TNL(40,1) = .0000002842 !.0000004917
201 SOL8TNL(40,1) = .0000017380 !.0000030250
202 TET4TNL(40,1) = .0000002744 !.0000004015
203 TET10TNL(40,1)= .0000011310 !.0000013700
204! LAW 41
205 SOL1TNL(41,1) = .0000002894 !.0000004092
206 TET4TNL(41,1) = .0000003329 !.0000004068
207 TET10TNL(41,1)= .0000013090 !.0000017520
208! LAW 42
209 SOL1TNL(42,1) = .0000004594 !.0000006159
210 SOL8TNL(42,1) = .0000030370 !.0000041410
211 TET4TNL(42,1) = .0000004061 !.0000006126
212 TET10TNL(42,1)= .0000017340 !.0000023340
213C M = 1
214 SOL1TNL(42,2) = .0000005017 !.0000006763
215 SOL8TNL(42,2) = .0000034450 !.0000045820
216 TET4TNL(42,2) = .0000005178 !.0000006537
217 TET10TNL(42,2)= .0000019600 !.0000023430
218C M = 2
219 SOL1TNL(42,3) = .0000005585 !.0000007634
220 SOL8TNL(42,3) = .0000036910 !.0000051210
221 TET4TNL(42,3) = .0000005400 !.0000006823
222 TET10TNL(42,3)= .0000022000 !.0000027020
223! LAW 44
224 SOL1TNL(44,1) = .0000003162
225 SOL8TNL(44,1) = .0000025400
226 TET4TNL(44,1) = .0000003162
227 TET10TNL(44,1)= .0000012600
228! LAW 48
229 SOL1TNL(48,1) = .0000002001 !.0000002994
230 SOL8TNL(48,1) = .0000011410 !.0000015830
231 TET4TNL(48,1) = .0000001934 !.0000002660
232 TET10TNL(48,1)= .0000008028 !.0000009488
233! LAW 50
234 SOL1TNL(50,1) = .0000003576 !.0000004482
235 SOL8TNL(50,1) = .0000022190 !.0000027100
236 TET4TNL(50,1) = .0000003318 !.0000004178
237 TET10TNL(50,1)= .0000013960 !.0000018130
238! LAW 52
239 SOL1TNL(52,1) = .0000001956 !.0000002826
240 SOL8TNL(52,1) = .0000010740 !.0000015560
241 TET4TNL(52,1) = .0000001982 !.0000002650
242 TET10TNL(52,1)= .0000007707 !.0000010590
243! LAW 53
244 SOL1TNL(53,1) = .0000002705 !.0000004006
245 SOL8TNL(53,1) = .0000014390 !.0000022760
246 TET4TNL(53,1) = .0000002479 !.0000003889
247 TET10TNL(53,1)= .0000009377 !.0000014540
248! LAW 59 (user solid property) IGT=43
249 SOL1TNL(59,1) = .0000001809 !.0000002207
250 SOL8TNL(59,1) = .0000001809 !.0000002207
251 TET4TNL(59,1) = .0000001809 !.0000002207
252 TET10TNL(59,1)= .0000001809 !.0000002207
253! LAW 60
254 SOL1TNL(60,1) = .0000002901 !.0000004239
255 SOL8TNL(60,1) = .0000017500 !.0000024010
256 TET4TNL(60,1) = .0000002759 !.0000003718
257 TET10TNL(60,1)= .0000011720 !.0000014290
258! LAW 62
259 ! M = 0
260 SOL1TNL(62,1) = .0000006549 !.0000009680
261 SOL8TNL(62,1) = .0000052300 !.0000068900
262 TET4TNL(62,1) = .0000006193 !.0000009896
263 TET10TNL(62,1)= .0000026920 !.0000039430
264 ! M = 1
265 SOL1TNL(62,2) = .0000007000 !.0000010850
266 SOL8TNL(62,2) = .0000056000 !.0000078600
267 TET4TNL(62,2) = .0000006892 !.0000010760
268 TET10TNL(62,2)= .000002692 !.0000043120
269 ! M = 2
270 SOL1TNL(62,3) = .0000007570 !.0000011420
271 SOL8TNL(62,3) = .0000059760 !.0000084600
272 TET4TNL(62,3) = .0000007288 !.0000011150
273 TET10TNL(62,3)= .0000027780 !.0000045435
274! LAW 65
275 SOL1TNL(65,1) = .0000002856 !.0000004282
276 SOL8TNL(65,1) = .0000021460 !.0000030270
277 TET4TNL(65,1) = .0000002801 !.0000003897
278 TET10TNL(65,1)= .0000011220 !.0000015160
279! LAW 66
280 SOL1TNL(66,1) = .0000002659 !.0000005367
281 SOL8TNL(66,1) = .0000018670 !.0000031200
282 TET4TNL(66,1) = .0000002465 !.0000003860
283 TET10TNL(66,1)= .0000009512 !.0000015470
284! LAW 68
285 SOL1TNL(68,1) = .0000004811 !.0000006413
286! SOL8TNL(68,1) = .0000028940
287 TET4TNL(68,1) = .0000004870 !.0000006475
288 TET10TNL(68,1)= .0000019480 !.0000025900
289! LAW 69
290 ! M = 0 (no visc)
291 SOL1TNL(69,1) = .0000004161 !.0000005245
292 SOL8TNL(69,1) = .0000028460 !.0000040350
293 TET4TNL(69,1) = .0000003590 !.0000005135
294 TET10TNL(69,1)= .0000014360 !.0000020020
295 ! M = 1 (visc)
296 SOL1TNL(69,2) = .0000004501 !.0000006120
297 SOL8TNL(69,2) = .0000035560 !.0000048180
298 TET4TNL(69,2) = .0000004622 !.0000006347
299 TET10TNL(69,2)= .0000018040 !.0000024700
300 ! M = 2 (visc)
301 SOL1TNL(69,3) = .0000005034 !.0000006370
302 SOL8TNL(69,3) = .0000036270 !.0000051600
303 TET4TNL(69,3) = .0000004980 !.0000006740
304 TET10TNL(69,3)= .0000019950 !.0000025500
305! LAW 70
306! 1 funct / 22 points
307 SOL1TNL(70,1) = .0000002656 !.0000003853
308 SOL8TNL(70,1) = .0000018990 !.0000026830
309 TET4TNL(70,1) = .0000002435 !.0000003395
310 TET10TNL(70,1)= .0000009684 !.0000013970
311! overcost number of point
312 SOL1TNL(70,2) = .0000000014
313 SOL8TNL(70,2) = .0000000112
314 TET4TNL(70,2) = .0000000013
315 TET10TNL(70,2)= .0000000052
316! overcost number of function
317 SOL1TNL(70,3) = .0000000008
318 SOL8TNL(70,3) = .0000000064
319 TET4TNL(70,3) = .0000000008
320 TET10TNL(70,3)= .0000000032
321! LAW 71
322 SOL1TNL(71,1) = .0000004733 !.0000006332
323 SOL8TNL(71,1) = .0000016480 !.0000021590
324 TET4TNL(71,1) = .0000003537 !.0000004840
325 TET10TNL(71,1)= .0000014450 !.0000018860
326! LAW 72
327 SOL1TNL(72,1) = .0000003084 !.0000004335
328 SOL8TNL(72,1) = .0000022690 !.0000032810
329 TET4TNL(72,1) = .0000002603 !.0000003947
330 TET10TNL(72,1)= .0000011850 !.0000016170
331! LAW 74
332 SOL1TNL(74,1) = .0000002712 !.0000003704
333 SOL8TNL(74,1) = .0000019300 !.0000027610
334 TET4TNL(74,1) = .0000002488 !.0000003402
335 TET10TNL(74,1)= .0000009575 !.0000013590
336! LAW 76
337 SOL1TNL(76,1) = .0000001949 !.0000002745
338 SOL8TNL(76,1) = .0000015130 !.0000021280
339 TET4TNL(76,1) = .0000002020 !.0000002739
340 TET10TNL(76,1)= .0000007802 !.0000011500
341! LAW 77
342 SOL1TNL(77,1) = .0000003605 !.0000004887
343 SOL8TNL(77,1) = .0000028100 !.0000038100
344 TET4TNL(77,1) = .0000003500 !.0000005200
345 TET10TNL(77,1)= .0000014000 !.0000020120
346! LAW 78
347 SOL1TNL(78,1) = .0000002057 !.0000003087
348 SOL8TNL(78,1) = .0000015950 !.0000023090
349 TET4TNL(78,1) = .0000002297 !.0000002902
350 TET10TNL(78,1)= .0000008553 !.0000012520
351! LAW 79
352 SOL1TNL(79,1) = .0000002027 !.0000003060
353 SOL8TNL(79,1) = .0000014940 !.0000020220
354 TET4TNL(79,1) = .0000002031 !.0000002684
355 TET10TNL(79,1)= .0000007830 !.0000011400
356! LAW 81
357 SOL1TNL(81,1) = .0000003878 !.0000004223
358 SOL8TNL(81,1) = .0000028400 !.0000031430
359 TET4TNL(81,1) = .0000003539 !.0000002750
360 TET10TNL(81,1)= .0000014370 !.0000012470
361! LAW 82
362 ! M = 1
363 SOL1TNL(82,1) = .0000005600
364 SOL8TNL(82,1) = .0000045000
365 TET4TNL(82,1) = .0000005500
366 TET10TNL(82,1)= .0000022000
367 ! M = 2
368 SOL1TNL(82,2) = .0000006800
369 SOL8TNL(82,2) = .0000054400
370 TET4TNL(82,2) = .0000006700
371 TET10TNL(82,2)= .0000026900
372 ! M = 3
373 SOL1TNL(82,3) = .0000008800
374 SOL8TNL(82,3) = .0000071000
375 TET4TNL(82,3) = .0000008800
376 TET10TNL(82,3)= .0000036000
377! LAW 83 (user solid property) IGT=43
378 SOL1TNL(83,1) = .0000010200 !.0000013690
379 SOL8TNL(83,1) = .0000010200 !.0000013690
380 TET4TNL(83,1) = .0000010200 !.0000013690
381 TET10TNL(83,1)= .0000010200 !.0000013690
382! LAW 84
383 SOL1TNL(84,1) = .0000001938 !.0000002771
384 SOL8TNL(84,1) = .0000014240 !.0000018140
385 TET4TNL(84,1) = .0000001955 !.0000002489
386 TET10TNL(84,1)= .0000007306 !.0000007861
387! LAW 88
388 SOL1TNL(88,1) = .0000051800
389 SOL8TNL(88,1) = .0000421000
390 TET4TNL(88,1) = .0000051800
391 TET10TNL(88,1)= .0000210000
392! LAW 90
393 SOL1TNL(90,1) = .0000007640
394 SOL8TNL(90,1) = .0000053480
395 TET4TNL(90,1) = .0000007640
396 TET10TNL(90,1)= .0000029032
397! LAW 90 + ISMSTR=10 or 12
398 SOL1TNL(90,2) = .0000009110
399 SOL8TNL(90,2) = .0000065000
400 TET4TNL(90,2) = .0000009110
401 TET10TNL(90,2)= .0000034900
402! LAW 92
403 SOL1TNL(92,1) = .0000006540
404 SOL8TNL(92,1) = .0000052800
405 TET4TNL(92,1) = .0000006540
406 TET10TNL(92,1)= .0000026100
407! LAW 92 + ISMSTR=10 or 12
408 SOL1TNL(92,2) = .0000006900
409 SOL8TNL(92,2) = .0000054000
410 TET4TNL(92,2) = .0000006900
411 TET10TNL(92,2)= .0000028500
412! LAW 94
413 SOL1TNL(94,1) = .0000003650
414 SOL8TNL(94,1) = .0000029200
415 TET4TNL(94,1) = .0000003650
416 TET10TNL(94,1)= .0000014600
417! LAW 94 + ISMSTR=10 or 12
418 SOL1TNL(94,2) = .0000003900
419 SOL8TNL(94,2) = .0000033900
420 TET4TNL(94,2) = .0000003900
421 TET10TNL(94,2)= .0000016900
422
423! LAW 95
424 SOL1TNL(95,1) = .0000004170 !.0000005950
425 SOL8TNL(95,1) = .0000031200 !.0000045150
426 TET4TNL(95,1) = .0000004068 !.0000005905
427 TET10TNL(95,1)= .0000015740 !.0000021230
428! LAW 97
429 SOL1TNL(97,1) = .0000002348
430 SOL8TNL(97,1) = .0000018790
431 TET4TNL(97,1) = .0000002200
432 TET10TNL(97,1)= .0000008921
433! LAW 100 : material cost = base cost + viscosity cost + N * network cost
434! LAW 100 (base)
435 SOL1TNL(100,1) = .0000001753 !.0000002280
436 SOL8TNL(100,1) = .0000001753 !.0000002280
437 TET4TNL(100,1) = .0000001753 !.0000002280
438 TET10TNL(100,1)= .0000001753 !.0000002280
439
440! LAW 100 (viscosity)
441 SOL1TNL(100,2) = .0000001080 !.0000001390
442 SOL8TNL(100,2) = .0000001080 !.0000001390
443 TET4TNL(100,2) = .0000001080 !.0000001390
444 TET10TNL(100,2)= .0000001080 !.0000001390
445
446! LAW 100 (network)
447 SOL1TNL(100,3) = .0000001550 !.0000002000
448 SOL8TNL(100,3) = .0000001550 !.0000002000
449 TET4TNL(100,3) = .0000001550 !.0000002000
450 TET10TNL(100,3)= .0000001550 !.0000002000
451! LAW 101
452 SOL1TNL(101,1) = .0000058500
453 SOL8TNL(101,1) = .0000468000
454 TET4TNL(101,1) = .0000058500
455 TET10TNL(101,1)= .0000234900
456! LAW 102
457 SOL1TNL(102,1) = .0000002074
458 SOL8TNL(102,1) = .0000014790
459 TET4TNL(102,1) = .0000001970
460 TET10TNL(102,1)= .0000007800
461! LAW 103
462 SOL1TNL(103,1) = .0000002258
463 SOL8TNL(103,1) = .0000016170
464 TET4TNL(103,1) = .0000002150
465 TET10TNL(103,1)= .0000008450
466! LAW 104
467 ! Nice algo
468 SOL1TNL(104,1) = .0000002095
469 SOL8TNL(104,1) = .0000002095
470 TET4TNL(104,1) = .0000002095
471 TET10TNL(104,1) = .0000002095
472 ! Newton algo
473 SOL1TNL(104,2) = .0000003078
474 SOL8TNL(104,2) = .0000003078
475 TET4TNL(104,2) = .0000003078
476 TET10TNL(104,2) = .0000003078
477 ! iloc=1
478 SOL_OPTION(1,1) = .0000000600
479 ! iloc=2
480 SOL_OPTION(1,2) = .0000000700
481 ! iloc=3
482 SOL_OPTION(1,3) = .0000001200
483! LAW 106
484! SOL1TNL(106,1) = .0000002493
485! SOL8TNL(106,1) = .0000018090
486! TET4TNL(106,1) = .0000002246
487! TET10TNL(106,1)= .0000009194
488! LAW 187
489 SOL1TNL(187,1) = .0000002973
490 SOL8TNL(187,1) = .0000021260
491 TET4TNL(187,1) = .0000010430
492 TET10TNL(187,1)= .0000018088
493
494 IF (IPARI0.EQ.1) THEN
495! SOLTELT :
496! 1 --> ISOLID = 1
497! 2 --> ISOLID = 2
498! 3 --> ISOLID = 24
499! 4 --> ISOLID = 12
500! 5 --> ISOLID = 14 HA8 with NPT=2*2*2 (fixed part of element cost)
501! 6 --> ISOLID = 14 HA8 (overcost for NPT > 2*2*2)
502! HA8 cost : COST = HA8 (222) cost + (NPT-8) * overcost
503! 7 --> ISOLID = 17 H8C
504! 8 --> IGTYP = 43
505! 9 --> ISOLID = 18 H8C
506! 10--> IGTYP = Thickshell 20 or 21 or 22
507! 11--> JHBE=15 : fixed part of element cost
508! 12--> JHBE=15 : overcost for N>1
509! NLOCAL_OPTION
510! 1 --> sol
511! 2 --> shell
512!
513! TET4ELT :
514! 1 --> TETRA 4
515! 2 --> TETRA 4 with IROT = 1
516! TETRA10ELT
517! 1 --> TETRA10
518!
519! ELM_OVER_COST :
520! 1 --> ismstr=10 or 12
521
522 SOLTELT(1) =.0000005799 !.0000009528
523 SOLTELT(2) =.0000005446 !.0000008540
524 SOLTELT(3) =.0000006983 !.0000013420
525 SOLTELT(4) =.0000011707 !.0000015180
526 SOLTELT(5) =.0000030407 !.0000041200
527 SOLTELT(6) =.0000003003 !.0000003812
528 SOLTELT(7) =.0000018808 !.0000026600
529 SOLTELT(8) =.0000008152 !.0000012750
530 SOLTELT(9) =.0000023960
531 SOLTELT(10) =.0000003461
532 SOLTELT(11) =.0000009980
533 SOLTELT(12) =.0000000854
534
535 TET4TELT(1) =.0000002920 !.0000004744
536 TET4TELT(2) =.0000013046 !.0000016460
537 TET10TELT(1)=.0000014747 !.0000018340
538 NLOCAL_OPTION(1) = .0000001500
539
540 ELM_OVER_COST(1) = .0000000780
541 ELSE
542 SOLTELT(1) =.0000004832
543 SOLTELT(2) =.0000004538
544 SOLTELT(3) =.0000005818
545 SOLTELT(4) =.0000009755
546 SOLTELT(5) =.0000025400
547 SOLTELT(6) =.0000002500
548 SOLTELT(7) =.0000015700
549 SOLTELT(8) =.0000006793
550 SOLTELT(9) =.0000019970
551 SOLTELT(10) =.0000002914
552 SOLTELT(11) =.0000008400
553 SOLTELT(12) =.0000000702
554 TET4TELT(1) =.0000002433
555 TET4TELT(2) =.0000010870
556 TET10TELT(1)=.0000012200
557 NLOCAL_OPTION(1) = .0000001500
558
559 ELM_OVER_COST(1) = .0000000650
560 ENDIF
561
562! -------------
563! Coques
564! -------------
565
566
567
568! Tableau SHTNL
569! ---------------
570!
571! Temps par elt/cycle pour les materieux coques
572! SHTNL(I,J,K)=
573! 1er Champ : numero du materiaux
574! 2 champs : nbre de pts d integration, de taille 0:4
575! 4 sont consideres :
576! * 0 : 0pt d integ
577! * 1 : 1pt d integ
578! * 2 : 3pt d integ
579! * 3 : 5pt d integ
580! Les autres temps pour les NPT, sont calculees par interpolation avec ces 3
581! valeurs.
582!
583! 3e champ : option particuliere au besoin
584!
585! -------------
586! Triangles
587! -------------
588
589! Tableau TRITNL
590! ---------------
591!
592! Temps par elt/cycle pour les materieux coques
593! TRITNL(I,J,K)=
594! 1er Champ : numero du materiaux
595! 2e champ : nbre de pts d integration, de taille 0:4
596! 4 sont consideres :
597! * 0 : 0pt d integ
598! * 1 : 1pt d integ
599! * 2 : 3pt d integ
600! * 3 : 5pt d integ
601! Les autres temps pour les NPT, sont calculees par interpolation avec ces 3
602! valeurs.
603!
604! 3e champ : option particuliere au besoin
605!
606!
607! Tableau TRITELT
608! ---------------
609!
610! Temps par elt/cycle pour le traitement de l hourglas en fct de /IPARITH
611
612 DO I=1,3
613 DO J = 0, MAXLAW
614 SHTNL(J,0,I) = .0000001059 !.0000001480
615 SHTNL(J,1,I) = .0000001537
616 SHTNL(J,2,I) = .0000002827
617 SHTNL(J,3,I) = .0000003898
618C
619 TRITNL(J,0,I)= .0000000959 !.0000001373
620 TRITNL(J,1,I)= .0000001358
621 TRITNL(J,2,I)= .0000002297
622 TRITNL(J,3,I)= .0000003223
623 ENDDO
624 ENDDO
625! special option for shell or shell3n
626 SHTNL_OPTION(1:10,1:5) = 0.
627! LAW 01
628 SHTNL(1,0,1) = .0000000606 !.0000000897
629 SHTNL(1,1,1) = .0000000660 !.0000000811
630 SHTNL(1,2,1) = .0000000680 !.0000001004
631 SHTNL(1,3,1) = .0000000690 !.0000001032
632C
633 TRITNL(1,0,1)= .0000000521 !.0000000999
634 TRITNL(1,1,1)= .0000000474 !.0000000797
635 TRITNL(1,2,1)= .0000000508 !.0000001024
636 TRITNL(1,3,1)= .0000000550 !.0000000999
637! LAW 02
638 SHTNL(2,0,1) = .0000001059 !.0000001480
639 SHTNL(2,1,1) = .0000001537
640 SHTNL(2,2,1) = .0000002827
641 SHTNL(2,3,1) = .0000003898
642C
643 TRITNL(2,0,1)= .0000000959 !.0000001373
644 TRITNL(2,1,1)= .0000001358
645 TRITNL(2,2,1)= .0000002297
646 TRITNL(2,3,1)= .0000003223
647! LAW 02 CC
648 SHTNL(2,0,2) = .0000001450 !.0000001606
649 SHTNL(2,1,2) = .0000001600
650 SHTNL(2,2,2) = .0000003066
651 SHTNL(2,3,2) = .0000004153
652C
653 TRITNL(2,0,2)= .0000001400
654 TRITNL(2,1,2)= .0000001520
655 TRITNL(2,2,2)= .0000002700
656 TRITNL(2,3,2)= .0000003900
657! LAW 15
658 SHTNL(15,0,1) = .0000002391
659 SHTNL(15,1,1) = .0000002391
660 SHTNL(15,2,1) = .0000005294
661 SHTNL(15,3,1) = .0000008390
662C
663 TRITNL(15,0,2)= .0000002414
664 TRITNL(15,1,2)= .0000002429
665 TRITNL(15,2,2)= .0000005681
666 TRITNL(15,3,2)= .0000008166
667
668! LAW 19
669 SHTNL(19,0,1) = .0000001566
670 SHTNL(19,1,1) = .0000001566
671 SHTNL(19,2,1) = .0000003025
672 SHTNL(19,3,1) = .0000004075
673C
674 TRITNL(19,0,1)= .0000001540
675 TRITNL(19,1,1)= .0000001586
676 TRITNL(19,2,1)= .0000002802
677 TRITNL(19,3,1)= .0000004140
678! LAW 22
679 SHTNL(22,0,1) = .0000001741
680 SHTNL(22,1,1) = .0000002228
681 SHTNL(22,2,1) = .0000004186
682 SHTNL(22,3,1) = .0000006076
683C
684 TRITNL(22,0,1)= .0000001559 !.0000002209
685 TRITNL(22,1,1)= .0000001452 !.0000001903
686 TRITNL(22,2,1)= .0000002874 !.0000004236
687 TRITNL(22,3,1)= .0000004223 !.0000005610
688! LAW 25
689 SHTNL(25,0,1) = .0000002534 !.0000002511
690 SHTNL(25,1,1) = .0000002534 !.0000002509
691 SHTNL(25,2,1) = .0000005014 !.0000005094
692 SHTNL(25,3,1) = .0000007885 !.0000008026
693C
694 TRITNL(25,0,1)= .0000002519 !.0000002519
695 TRITNL(25,1,1)= .0000002455 !.0000002455
696 TRITNL(25,2,1)= .0000005096 !.0000005096
697 TRITNL(25,3,1)= .0000007655 !.0000007655
698! LAW 25 Crasurv (/TYPE11)
699 SHTNL(25,0,2) = .0000004094 !.0000004487
700 SHTNL(25,1,2) = .0000004094 !.0000004487
701 SHTNL(25,2,2) = .0000008112 !.0000009845
702 SHTNL(25,3,2) = .0000013858 !.0000015870
703C
704 TRITNL(25,0,2)= .0000004000 !.0000003866
705 TRITNL(25,1,2)= .0000004000 !.0000003786
706 TRITNL(25,2,2)= .0000007900 !.0000009420
707 TRITNL(25,3,2)= .0000013500 !.0000013980
708! LAW 27 (0pt incompatible)
709 SHTNL(27,0,1) = .0000002399
710 SHTNL(27,1,1) = .0000002399
711 SHTNL(27,2,1) = .0000005425
712 SHTNL(27,3,1) = .0000008176
713C
714 TRITNL(27,0,1)= .0000002148
715 TRITNL(27,1,1)= .0000002148
716 TRITNL(27,2,1)= .0000004802
717 TRITNL(27,3,1)= .0000007094
718! LAW 32
719 SHTNL(32,0,1) = .0000001526
720 SHTNL(32,1,1) = .0000001526
721 SHTNL(32,2,1) = .0000003010
722 SHTNL(32,3,1) = .0000004443
723C
724 TRITNL(32,0,1)= .0000001743
725 TRITNL(32,1,1)= .0000001743
726 TRITNL(32,2,1)= .0000003877
727 TRITNL(32,3,1)= .0000005818
728! LAW 35 (0pt incompatible)
729 SHTNL(35,0,1) = .0000002142 !.0000002696
730 SHTNL(35,1,1) = .0000002142 !.0000002696
731 SHTNL(35,2,1) = .0000003850 !.0000005399
732 SHTNL(35,3,1) = .0000006158 !.0000008694
733C
734 TRITNL(35,0,1)= .0000002010 !.0000002629
735 TRITNL(35,1,1)= .0000002010 !.0000002629
736 TRITNL(35,2,1)= .0000003750 !.0000005365
737 TRITNL(35,3,1)= .0000006010 !.0000009575
738! LAW 36 - 1 fct
739 SHTNL(36,0,1) = .0000001685 !.0000002453
740 SHTNL(36,1,1) = .0000002119 !.0000002975
741 SHTNL(36,2,1) = .0000003913 !.0000005370
742 SHTNL(36,3,1) = .0000005846 !.0000007957
743C
744 TRITNL(36,0,1)= .0000001572 !.0000002315
745 TRITNL(36,1,1)= .0000002020 !.0000003005
746 TRITNL(36,2,1)= .0000003910 !.0000005141
747 TRITNL(36,3,1)= .0000005784 !.0000008021
748! LAW 36 - 5 fct
749 SHTNL(36,0,2) = .0000001700 !.0000002832
750 SHTNL(36,1,2) = .0000002210 !.0000003655
751 SHTNL(36,2,2) = .0000003980 !.0000005902
752 SHTNL(36,3,2) = .0000005900 !.0000008500
753C
754 TRITNL(36,0,2)= .0000001650 !.0000002975
755 TRITNL(36,1,2)= .0000002060 !.0000003608
756 TRITNL(36,2,2)= .0000003960 !.0000005945
757 TRITNL(36,3,2)= .0000005860 !.0000008260
758! LAW 36 - 10 fct
759 SHTNL(36,0,3) = .0000001750 !.0000003180
760 SHTNL(36,1,3) = .0000002250 !.0000003850
761 SHTNL(36,2,3) = .0000004010 !.0000006500
762 SHTNL(36,3,3) = .0000005950 !.0000009020
763C
764 TRITNL(36,0,3)= .0000001710 !.0000003060
765 TRITNL(36,1,3)= .0000002100 !.0000003820
766 TRITNL(36,2,3)= .0000004020 !.0000006540
767 TRITNL(36,3,3)= .0000005950 !.0000009000
768! LAW 42
769 SHTNL(42,0,1) = .0000004390 !.0000005524
770 SHTNL(42,1,1) = .0000004390 !.0000005524
771 SHTNL(42,2,1) = .0000010540 !.0000013820
772 SHTNL(42,3,1) = .0000016890 !.0000022760
773! LAW 42 M = 1
774 SHTNL(42,0,2) = .0000015080 !.0000015100
775 SHTNL(42,1,2) = .0000015080 !.0000015100
776 SHTNL(42,2,2) = .0000041230 !.0000045280
777 SHTNL(42,3,2) = .0000060650 !.0000074640
778! LAW 42 M = 2
779 SHTNL(42,0,3) = .0000016050 !.0000017720
780 SHTNL(42,1,3) = .0000016050 !.0000017720
781 SHTNL(42,2,3) = .0000042500 !.0000053370
782 SHTNL(42,3,3) = .0000062320 !.0000088600
783! LAW 42
784 TRITNL(42,0,1)= .0000003891 !.0000005586
785 TRITNL(42,1,1)= .0000003891 !.0000005586
786 TRITNL(42,2,1)= .0000010570 !.0000015540
787 TRITNL(42,3,1)= .0000016780 !.0000022710
788! LAW 42 M = 1
789 TRITNL(42,0,2)= .0000015100 !.0000015000
790 TRITNL(42,1,2)= .0000015100 !.0000015000
791 TRITNL(42,2,2)= .0000041690 !.0000044580
792 TRITNL(42,3,2)= .0000069380 !.0000073400
793! LAW 42 M = 2
794 TRITNL(42,0,3)= .0000016000 !.0000017200
795 TRITNL(42,1,3)= .0000016000 !.0000017200
796 TRITNL(42,2,3)= .0000043000 !.0000052700
797 TRITNL(42,3,3)= .0000063020 !.0000087000
798! LAW 43
799 SHTNL(43,0,1) = .0000001555 !.0000002796
800 SHTNL(43,1,1) = .0000001555 !.0000002918
801 SHTNL(43,2,1) = .0000004576 !.0000006119
802 SHTNL(43,3,1) = .0000007020 !.0000009136
803C
804 TRITNL(43,0,1)= .0000001500 !.0000002856
805 TRITNL(43,1,1)= .0000001500 !.0000002743
806 TRITNL(43,2,1)= .0000004476 !.0000006279
807 TRITNL(43,3,1)= .0000006920 !.0000008749
808! LAW 44
809 SHTNL(44,0,1) = .0000002156 !.0000002640
810 SHTNL(44,1,1) = .0000002156 !.0000002640
811 SHTNL(44,2,1) = .0000004318 !.0000005714
812 SHTNL(44,3,1) = .0000006232 !.0000008505
813C
814 TRITNL(44,0,1)= .0000001982 !.0000002802
815 TRITNL(44,1,1)= .0000001982 !.0000002802
816 TRITNL(44,2,1)= .0000004025 !.0000005076
817 TRITNL(44,3,1)= .0000006002 !.0000007598
818! LAW 48
819 SHTNL(48,0,1) = .0000002079 !.0000002811
820 SHTNL(48,1,1) = .0000002079 !.0000002811
821 SHTNL(48,2,1) = .0000003971 !.0000005285
822 SHTNL(48,3,1) = .0000006038 !.0000007925
823C
824 TRITNL(48,0,1)= .0000001915 !.0000002764
825 TRITNL(48,1,1)= .0000001915 !.0000002764
826 TRITNL(48,2,1)= .0000003940 !.0000005447
827 TRITNL(48,3,1)= .0000005780 !.0000007677
828! LAW 52
829 SHTNL(52,0,1) = .0000001984 !.0000002313
830 SHTNL(52,1,1) = .0000001984 !.0000002313
831 SHTNL(52,2,1) = .0000003811 !.0000004775
832 SHTNL(52,3,1) = .0000005483 !.0000007161
833C
834 TRITNL(52,0,1)= .0000001848 !.0000002422
835 TRITNL(52,1,1)= .0000001848 !.0000002422
836 TRITNL(52,2,1)= .0000003654 !.0000004612
837 TRITNL(52,3,1)= .0000005430 !.0000007237
838! LAW 57
839 SHTNL(57,0,1) = .0000003247 !.0000005039
840 SHTNL(57,1,1) = .0000003247 !.0000005039
841 SHTNL(57,2,1) = .0000007464 !.0000013330
842 SHTNL(57,3,1) = .0000011440 !.0000021130
843C
844 TRITNL(57,0,1)= .0000003169 !.0000005568
845 TRITNL(57,1,1)= .0000003169 !.0000005568
846 TRITNL(57,2,1)= .0000007755 !.0000014330
847 TRITNL(57,3,1)= .0000011380 !.0000022840
848! LAW 58
849! No visc (FLAGUL==0)
850 SHTNL(58,0,1) = .0000003467 !.0000005435
851 SHTNL(58,1,1) = .0000003467 !.0000006044
852 SHTNL(58,2,1) = .0000010060 !.0000017910
853 SHTNL(58,3,1) = .0000015890 !.0000028510
854C
855 TRITNL(58,0,1)= .0000003367 !.0000005733
856 TRITNL(58,1,1)= .0000003367 !.0000005733
857 TRITNL(58,2,1)= .0000009850 !.0000015320
858 TRITNL(58,3,1)= .0000014800 !.0000023640
859! LAW 60
860 SHTNL(60,0,1) = .0000002189 !.0000003498
861 SHTNL(60,1,1) = .0000002731 !.0000003992
862 SHTNL(60,2,1) = .0000005789 !.0000008897
863 SHTNL(60,3,1) = .0000008774 !.0000013350
864!
865 TRITNL(60,0,1)= .0000002565 !.0000003850
866 TRITNL(60,1,1)= .0000002737 !.0000003997
867 TRITNL(60,2,1)= .0000006308 !.0000009103
868 TRITNL(60,3,1)= .0000009732 !.0000014110
869! LAW 62
870 ! M = 0
871 SHTNL(62,0,1) = .0000005627 !.0000008089
872 SHTNL(62,1,1) = .0000005627 !.0000008089
873 SHTNL(62,2,1) = .0000014310 !.0000021360
874 SHTNL(62,3,1) = .0000023190 !.0000034950
875 ! M = 1
876 SHTNL(62,0,2) = .0000010350 !.0000014640
877 SHTNL(62,1,2) = .0000010350 !.0000014640
878 SHTNL(62,2,2) = .0000028030 !.0000041200
879 SHTNL(62,3,2) = .0000046550 !.0000067100
880 ! M = 2
881 SHTNL(62,0,3) = .0000012120 !.0000017390
882 SHTNL(62,1,3) = .0000012120 !.0000017390
883 SHTNL(62,2,3) = .0000033150 !.0000048970
884 SHTNL(62,3,3) = .0000055170 !.0000079110
885C
886 ! M = 0
887 TRITNL(62,0,1)= .0000005550 !.0000008167
888 TRITNL(62,1,1)= .0000005550 !.0000008167
889 TRITNL(62,2,1)= .0000014490 !.0000021660
890 TRITNL(62,3,1)= .0000022850 !.0000034760
891 ! M = 1
892 TRITNL(62,0,1)= .0000010070 !.0000014950
893 TRITNL(62,1,1)= .0000010070 !.0000014950
894 TRITNL(62,2,1)= .0000028340 !.0000040800
895 TRITNL(62,3,1)= .0000046390 !.0000067340
896 ! M = 2
897 TRITNL(62,0,1)= .0000012110 !.0000017140
898 TRITNL(62,1,1)= .0000012110 !.0000017140
899 TRITNL(62,2,1)= .0000033420 !.0000048680
900 TRITNL(62,3,1)= .0000054980 !.0000079480
901! LAW 63
902 SHTNL(63,0,1) = .0000002200 !.0000002904
903 SHTNL(63,1,1) = .0000002200 !.0000002904
904 SHTNL(63,2,1) = .0000004518 !.0000006582
905 SHTNL(63,3,1) = .0000006700 !.0000009983
906C
907 TRITNL(63,0,1)= .0000001996 !.0000002766
908 TRITNL(63,1,1)= .0000001996 !.0000002766
909 TRITNL(63,2,1)= .0000004432 !.0000005871
910 TRITNL(63,3,1)= .0000006509 !.0000008944
911! LAW 64
912 SHTNL(64,0,1) = .0000002615 !.0000003364
913 SHTNL(64,1,1) = .0000002615 !.0000003364
914 SHTNL(64,2,1) = .0000005904 !.0000007866
915 SHTNL(64,3,1) = .0000009000 !.0000011920
916C
917 TRITNL(64,0,1)= .0000002433 !.0000003426
918 TRITNL(64,1,1)= .0000002433 !.0000003426
919 TRITNL(64,2,1)= .0000005692 !.0000007550
920 TRITNL(64,3,1)= .0000008979 !.0000011970
921! LAW 65
922 SHTNL(65,0,1) = .0000002530 !.0000003350
923 SHTNL(65,1,1) = .0000002530 !.0000003350
924 SHTNL(65,2,1) = .0000005146 !.0000007280
925 SHTNL(65,3,1) = .0000007834 !.0000011240
926C
927 TRITNL(65,0,1)= .0000002375 !.0000003324
928 TRITNL(65,1,1)= .0000002375 !.0000003324
929 TRITNL(65,2,1)= .0000005089 !.0000007364
930 TRITNL(65,3,1)= .0000007783 !.0000011200
931! LAW 66
932 SHTNL(66,0,1) = .0000002521 !.0000003310
933 SHTNL(66,1,1) = .0000002521 !.0000003310
934 SHTNL(66,2,1) = .0000004955 !.0000007200
935 SHTNL(66,3,1) = .0000007601 !.0000011010
936C
937 TRITNL(66,0,1)= .0000002317 !.0000003725
938 TRITNL(66,1,1)= .0000002317 !.0000003666
939 TRITNL(66,2,1)= .0000004875 !.0000007924
940 TRITNL(66,3,1)= .0000007270 !.0000011650
941! LAW 69
942 ! M = 0 (no visc)
943 SHTNL(69,0,1) = .0000008395 !.0000013330
944 SHTNL(69,1,1) = .0000008395 !.0000013330
945 SHTNL(69,2,1) = .0000022420 !.0000032970
946 SHTNL(69,3,1) = .0000036290 !.0000054280
947 ! M = 1 (visc)
948 SHTNL(69,0,2) = .0000009622 !.0000014210
949 SHTNL(69,1,2) = .0000009622 !.0000014210
950 SHTNL(69,2,2) = .0000026220 !.0000039620
951 SHTNL(69,3,2) = .0000041070 !.0000064560
952 ! M = 2 (visc)
953 SHTNL(69,0,3) = .0000010930 !.0000016580
954 SHTNL(69,1,3) = .0000010930 !.0000016580
955 SHTNL(69,2,3) = .0000029100 !.0000045270
956 SHTNL(69,3,3) = .0000049120 !.0000074700
957C
958 ! M = 0 (no visc)
959 TRITNL(69,0,1)= .0000007899 !.0000010400
960 TRITNL(69,1,1)= .0000007899 !.0000010400
961 TRITNL(69,2,1)= .0000022350 !.0000027280
962 TRITNL(69,3,1)= .0000035990 !.0000044630
963 ! M = 1
964 TRITNL(69,0,2)= .0000009697 !.0000012490
965 TRITNL(69,1,2)= .0000009697 !.0000012490
966 TRITNL(69,2,2)= .000002603 !.0000035000
967 TRITNL(69,3,2)= .0000042470 !.0000055580
968 ! M = 2
969 TRITNL(69,0,3)= .0000010700 !.0000014430
970 TRITNL(69,1,3)= .0000010700 !.0000014430
971 TRITNL(69,2,3)= .0000029360 !.0000040510
972 TRITNL(69,3,3)= .0000049240 !.0000065830
973! LAW 71 (incompatible npt=0)
974 SHTNL(71,0,1) = .0000003612 !.0000005878
975 SHTNL(71,1,1) = .0000003612 !.0000005878
976 SHTNL(71,2,1) = .0000011340 !.0000014690
977 SHTNL(71,3,1) = .0000017410 !.0000024250
978C
979 TRITNL(71,0,1)= .0000004469 !.0000006319
980 TRITNL(71,1,1)= .0000004469 !.0000006919
981 TRITNL(71,2,1)= .0000011820 !.0000015420
982 TRITNL(71,3,1)= .0000018630 !.0000024620
983! LAW 72
984 SHTNL(72,0,1) = .0000002688 !.0000003671
985 SHTNL(72,1,1) = .0000002688 !.0000003671
986 SHTNL(72,2,1) = .0000005694 !.0000007838
987 SHTNL(72,3,1) = .0000008790 !.0000012730
988C
989 TRITNL(72,0,1)= .0000002433 !.0000003591
990 TRITNL(72,1,1)= .0000002433 !.0000003591
991 TRITNL(72,2,1)= .0000005342 !.0000007633
992 TRITNL(72,3,1)= .0000008420 !.0000011460
993! LAW 73
994 SHTNL(73,0,1) = .0000002180 !.0000003456
995 SHTNL(73,1,1) = .0000002180 !.0000003689
996 SHTNL(73,2,1) = .0000004810 !.0000007752
997 SHTNL(73,3,1) = .0000007301 !.0000012010
998C
999 TRITNL(73,0,1)= .0000002209 !.0000003618
1000 TRITNL(73,1,1)= .0000002209 !.0000003371
1001 TRITNL(73,2,1)= .0000004888 !.0000007637
1002 TRITNL(73,3,1)= .0000007160 !.0000012300
1003! LAW 76
1004 SHTNL(76,0,1) = .0000004953 !.0000007383
1005 SHTNL(76,1,1) = .0000004953 !.0000007476
1006 SHTNL(76,2,1) = .0000012140 !.0000018800
1007 SHTNL(76,3,1) = .0000019470 !.0000029980
1008C
1009 TRITNL(76,0,1)= .0000004596 !.0000007008
1010 TRITNL(76,1,1)= .0000004596 !.0000007008
1011 TRITNL(76,2,1)= .0000011830 !.0000018360
1012 TRITNL(76,3,1)= .0000019440 !.0000029000
1013! LAW 78
1014 SHTNL(78,0,1) = .0000002184 !.0000002921
1015 SHTNL(78,1,1) = .0000002184 !.0000002921
1016 SHTNL(78,2,1) = .0000004568 !.0000005980
1017 SHTNL(78,3,1) = .0000006627 !.0000009682
1018C
1019 TRITNL(78,0,1)= .0000002040 !.0000002936
1020 TRITNL(78,1,1)= .0000002040 !.0000002939
1021 TRITNL(78,2,1)= .0000004389 !.0000005670
1022 TRITNL(78,3,1)= .0000006650 !.0000008430
1023! LAW 80
1024 SHTNL(80,0,1) = .0000010260 !.0000015110
1025 SHTNL(80,1,1) = .0000010260 !.0000014400
1026 SHTNL(80,2,1) = .0000027800 !.0000040060
1027 SHTNL(80,3,1) = .0000042040 !.0000065840
1028C
1029 TRITNL(80,0,1)= .0000009896 !.0000013660
1030 TRITNL(80,1,1)= .0000009896 !.0000013660
1031 TRITNL(80,2,1)= .0000027450 !.0000039870
1032 TRITNL(80,3,1)= .0000044780 !.0000063540
1033! LAW 82 (incompatible NPT=0)
1034! N = 1
1035 SHTNL(82,0,1) = .0000005243 !.0000006991
1036 SHTNL(82,1,1) = .0000005243 !.0000007022
1037 SHTNL(82,2,1) = .0000013350 !.0000017930
1038 SHTNL(82,3,1) = .0000021650 !.0000028300
1039!
1040 TRITNL(82,0,1)= .0000005115 !.0000006886
1041 TRITNL(82,1,1)= .0000005115 !.0000007115
1042 TRITNL(82,2,1)= .0000013590 !.0000017680
1043 TRITNL(82,3,1)= .0000021880 !.0000028740
1044! N = 2
1045 SHTNL(82,0,2) = .0000007921 !.0000010470
1046 SHTNL(82,1,2) = .0000007921 !.0000010430
1047 SHTNL(82,2,2) = .0000020750 !.0000027830
1048 SHTNL(82,3,2) = .0000034520 !.0000045710
1049!
1050 TRITNL(82,0,2)= .0000007918 !.0000010690
1051 TRITNL(82,1,2)= .0000007918 !.0000010330
1052 TRITNL(82,2,2)= .0000021550 !.0000028430
1053 TRITNL(82,3,2)= .0000035270 !.0000046350
1054! N = 3
1055 SHTNL(82,0,3) = .0000010730 !.0000013880
1056 SHTNL(82,1,3) = .0000010730 !.0000014040
1057 SHTNL(82,2,3) = .0000029110 !.0000038590
1058 SHTNL(82,3,3) = .0000048550 !.0000063230
1059!
1060 TRITNL(82,0,3)= .0000010500 !.0000014340
1061 TRITNL(82,1,3)= .0000010500 !.0000014220
1062 TRITNL(82,2,3)= .0000029440 !.0000038320
1063 TRITNL(82,3,3)= .0000048030 !.0000063450
1064! LAW 87
1065 SHTNL(87,0,1) = .0000002700 !.0000003592
1066 SHTNL(87,1,1) = .0000002700 !.0000003592
1067 SHTNL(87,2,1) = .0000005776 !.0000007907
1068 SHTNL(87,3,1) = .0000008880 !.0000011520
1069!
1070 TRITNL(87,0,1)= .0000002542 !.0000003754
1071 TRITNL(87,1,1)= .0000002542 !.0000003874
1072 TRITNL(87,2,1)= .0000005625 !.0000008403
1073 TRITNL(87,3,1)= .0000008637 !.0000012850
1074! LAW 93
1075 SHTNL(93,0,1) = .0000001974 !.0000002801
1076 SHTNL(93,1,1) = .0000001974 !.0000002801
1077 SHTNL(93,2,1) = .0000003624 !.0000005238
1078 SHTNL(93,3,1) = .0000005093 !.0000007912
1079!
1080 TRITNL(93,0,1)= .0000001775 !.0000002629
1081 TRITNL(93,1,1)= .0000001775 !.0000002629
1082 TRITNL(93,2,1)= .0000003555 !.0000005262
1083 TRITNL(93,3,1)= .0000005344 !.0000008110
1084! LAW 104
1085 ! Nice algo
1086 SHTNL(104,0,1) = .0000001691
1087 SHTNL(104,1,1) = .0000001691
1088 SHTNL(104,2,1) = .0000003493
1089 SHTNL(104,3,1) = .0000005263
1090 ! Newton algo
1091 SHTNL(104,0,2) = .0000002163
1092 SHTNL(104,1,2) = .0000002163
1093 SHTNL(104,2,2) = .0000004636
1094 SHTNL(104,3,2) = .0000007109
1095 ! Nice algo
1096 TRITNL(104,0,1) = .0000001691
1097 TRITNL(104,1,1) = .0000001691
1098 TRITNL(104,2,1) = .0000003493
1099 TRITNL(104,3,1) = .0000005263
1100 ! Newton algo
1101 TRITNL(104,0,2) = .0000002163
1102 TRITNL(104,1,2) = .0000002163
1103 TRITNL(104,2,2) = .0000004636
1104 TRITNL(104,3,2) = .0000007109
1105 ! iloc=1
1106 SHTNL_OPTION(1,1) = .0000000600
1107 ! iloc=2
1108 SHTNL_OPTION(1,2) = .0000000700
1109 ! iloc=3
1110 SHTNL_OPTION(1,3) = .0000001500
1111
1112C Temps de traitement de la partie elementaire sans la loi materielle fonction de IPARITH
1113 IF (IPARI0.EQ.1) THEN
1114! SHELT :
1115! 1-->5 --> /SHELL
1116! 6-->10 --> /SH_COMP or TYPE10
1117! 11-->15 --> /SH_SANDW or TYPE11
1118! 15-->20 --> /SH_ORTH or TYPE9
1119! 20-->25 --> /SH_FABR or TYPE16
1120! IDRILL = 0
1121! 1/6/11/16/21 --> Q4
1122! 2/7/12/17/22 --> QEPH
1123! 3/8/13/18/23 --> BATOZ
1124! IDRILL = 1
1125! none --> Q4
1126! 4/9 /14/19/24 --> QEPH
1127! 5/10/15/20/25 --> BATOZ
1128! TRIELT :
1129! 1 --> tri
1130! NLOCAL_OPTION
1131! 1 --> sol
1132! 2 --> shell
1133! /SHELL
1134 SHTELT(1) = .0000003760 !.0000005629
1135 SHTELT(2) = .0000005334 !.0000007050
1136 SHTELT(3) = .0000009042 !.0000013190
1137 SHTELT(4) = .0000007847
1138 SHTELT(5) = .0000014430
1139! /SH_COMP
1140 SHTELT(6) = .0000004436 !.0000007170
1141 SHTELT(7) = .0000006470 !.0000009040
1142 SHTELT(8) = .0000011140 !.0000015630
1143 SHTELT(9) = .0000009621
1144 SHTELT(10)= .0000015060
1145! /SH_SANDW
1146 SHTELT(11) = .0000004598 !.0000006329
1147 SHTELT(12) = .0000006824 !.0000008915
1148 SHTELT(13) = .0000012190 !.0000014100
1149 SHTELT(14) = .0000009744
1150 SHTELT(15) = .0000015670
1151! /SH_ORTH
1152 SHTELT(16) = .0000004210 !.0000004117
1153 SHTELT(17) = .0000005792 !.0000007741
1154 SHTELT(18) = .0000010670 !.0000014700
1155 SHTELT(19) = .0000008670
1156 SHTELT(20) = .0000015420
1157! /SH_FABR
1158 SHTELT(21) = .0000004556 !.0000006739
1159 SHTELT(22) = .0000005507 !.0000007792
1160 SHTELT(23) = .0000011250 !.0000014680
1161 SHTELT(24) = .0000009110
1162 SHTELT(25) = .0000015050
1163
1164! /TYPE51
1165 SHTELT(26) = .0000010400
1166 SHTELT(27) = .0000012200
1167 SHTELT(28) = .0000027200
1168 SHTELT(29) = .0000012400
1169 SHTELT(30) = .0000036500
1170
1171 NLOCAL_OPTION(2) = .0000000900
1172
1173 TRITELT(1)= .0000002951 !.0000004724
1174 ELSE
1175! /SHELL
1176 SHTELT(1) = .0000003133
1177 SHTELT(2) = .0000004445
1178 SHTELT(3) = .0000007535
1179 SHTELT(4) = .0000006540
1180 SHTELT(5) = .0000012025
1181! /SH_COMP
1182 SHTELT(6) = .0000003696
1183 SHTELT(7) = .0000005391
1184 SHTELT(8) = .0000009283
1185 SHTELT(9) = .0000008017
1186 SHTELT(10)= .0000012550
1187! /SH_SANDW
1188 SHTELT(11) = .0000003831
1189 SHTELT(12) = .0000005686
1190 SHTELT(13) = .0000010158
1191 SHTELT(14) = .0000008120
1192 SHTELT(15) = .0000013060
1193! /SH_ORTH
1194 SHTELT(16) = .0000003508
1195 SHTELT(17) = .0000004826
1196 SHTELT(18) = .0000008891
1197 SHTELT(19) = .0000007225
1198 SHTELT(20) = .0000012850
1199! /SH_FABR
1200 SHTELT(21) = .0000003796
1201 SHTELT(22) = .0000004586
1202 SHTELT(23) = .0000009375
1203 SHTELT(24) = .0000007591
1204 SHTELT(25) = .0000012541
1205
1206! /TYPE51
1207 SHTELT(26) = .0000008400
1208 SHTELT(27) = .0000009840
1209 SHTELT(28) = .0000021800
1210 SHTELT(29) = .0000010040
1211 SHTELT(30) = .0000029300
1212
1213 NLOCAL_OPTION(2) = .0000000900
1214
1215 TRITELT(1)= .0000002450
1216
1217 ENDIF
1218
1219! il y a un facteur multiplicatif sur les coques de Batoz, car on rentre 4x dans cmain3.
1220
1221 BATOZMULT = 3.8
1222! /PRONY or /LPRONY overcost
1223 VISC_PRONY_COST = .0000000133
1224
1225! *******************************************************************
1226! Failure model
1227! *******************************************************************
1228! Johnson : 1
1229! Tbutcher : 2
1230! Wilkins : 3
1231! user1 : 4
1232! user2 : 5
1233! user3 : 6
1234! FLD : 7 (for shell only)
1235! Spalling : 8 (for solid only)
1236! Wierzbicki : 9
1237! Tensstrain : 10
1238! Energy : 11
1239! ??? : 12
1240! Chang : 13 (for shell only)
1241! Hashin : 14
1242! Yamada : 15
1243! Puck : 16
1244! XFEM_Tbutcher : 17
1245! Lad_dama : 18 (for solid only)
1246! Power_da : 19
1247! Connect : 20 (for solid only)
1248! XFEM_Johns : 21
1249! XFEM_FLD : 22
1250! TAB1 : 23
1251! Orthstra : 24 (for shell only)
1252! NXT : 25 (for shell only)
1253! Snconnect : 26 (for solid only)
1254! EMC : 27 (for solid only)
1255! Wind : 28 (for shell only)
1256! Sahraei1 : 29 (for solid only)
1257! Biquadratic : 30
1258! _
1259! SHELL : 1--> 0 ou 1pt |
1260! 2--> 3 pts |--- MLN < 28 or ==32
1261! 3--> 5 pts _|
1262! 4--> 0 ou 1pt |
1263! 5--> 3 pts |--- MLN >= 28 or /=32
1264! 6--> 5 pts _|
1265!
1266! SOLIDE : 1 ----> MLN <28 or ==49
1267! 2 ----> MLN >=28 or /=49
1268! --------------------------------------------------------
1269 DO I=1,MAXFAIL
1270 DO J=1,3
1271 RUPTURE_SOL(I,J) = 0.
1272 RUPTURE_SOL12(I,J) = 0.
1273 RUPTURE_TET4(I,J) = 0.
1274 RUPTURE_TET10(I,J) = 0.
1275 ENDDO
1276 DO J=1,6
1277 RUPTURE_SHELL(I,J) = 0.
1278 ENDDO
1279 ENDDO
1280! --------------------
1281! Johnson
1282 ! MLN <28 or ==32
1283 RUPTURE_SHELL(1,1) = .00000007517
1284 RUPTURE_SHELL(1,2) = .00000019090
1285 RUPTURE_SHELL(1,3) = .00000031460
1286
1287 ! MLN >=28 or /=32
1288 RUPTURE_SHELL(1,4) = .00000003021
1289 RUPTURE_SHELL(1,5) = .00000007567
1290 RUPTURE_SHELL(1,6) = .00000011620
1291
1292 ! MLN <28 or ==49
1293 RUPTURE_SOL(1,1) = .00000008839
1294 RUPTURE_SOL12(1,1) = .00000017710
1295 RUPTURE_TET4(1,1) = .00000008021
1296 RUPTURE_TET10(1,1) = .00000033420
1297 ! MLN >=28 or /=49
1298 RUPTURE_SOL(1,2) = .00000002696
1299 RUPTURE_SOL12(1,2) = .00000015790
1300 RUPTURE_TET4(1,2) = .00000002663
1301 RUPTURE_TET10(1,2) = .00000010360
1302! --------------------
1303! Tbutcher
1304 ! MLN <28 or ==49
1305 RUPTURE_SHELL(2,1) = .00000004873
1306 RUPTURE_SHELL(2,2) = .00000010840
1307 RUPTURE_SHELL(2,3) = .00000016340
1308
1309 RUPTURE_SHELL(2,4) = .00000003994
1310 RUPTURE_SHELL(2,5) = .00000010130
1311 RUPTURE_SHELL(2,6) = .00000015690
1312
1313 ! MLN <28 or ==49
1314 RUPTURE_SOL(2,1) = .00000004718
1315 RUPTURE_SOL12(2,1) = .00000031390
1316 RUPTURE_TET4(2,1) = .00000010090
1317 RUPTURE_TET10(2,1) = .00000040130
1318
1319 ! MLN >=28 or /=49
1320 RUPTURE_SOL(2,2) = .00000004949
1321 RUPTURE_SOL12(2,2) = .00000030290
1322 RUPTURE_TET4(2,2) = .00000004449
1323 RUPTURE_TET10(2,2) = .00000016800
1324! --------------------
1325! Wilkins
1326 ! MLN <28 or ==49
1327 RUPTURE_SHELL(3,1) = .00000011660
1328 RUPTURE_SHELL(3,2) = .00000032760
1329 RUPTURE_SHELL(3,3) = .00000052680
1330
1331 RUPTURE_SHELL(3,4) = .00000011410
1332 RUPTURE_SHELL(3,5) = .00000032760
1333 RUPTURE_SHELL(3,6) = .00000052910
1334
1335 ! MLN <28 or ==49
1336 RUPTURE_SOL(3,1) = .00000023210
1337 RUPTURE_SOL12(3,1) = .00000131700
1338 RUPTURE_TET4(3,1) = .00000016440
1339 RUPTURE_TET10(3,1) = .00000075060
1340
1341 ! MLN >=28 or /=49
1342 RUPTURE_SOL(3,2) = .00000016870
1343 RUPTURE_SOL12(3,2) = .00000129800
1344 RUPTURE_TET4(3,2) = .00000011060
1345 RUPTURE_TET10(3,2) = .00000051000
1346! --------------------
1347! FLD (only for shell)
1348 ! MLN <28 or ==32
1349 RUPTURE_SHELL(7,1) = .00000010900
1350 RUPTURE_SHELL(7,2) = .00000021590
1351 RUPTURE_SHELL(7,3) = .00000031890
1352
1353 ! MLN >=28 or /=32
1354 RUPTURE_SHELL(7,4) = .00000007162
1355 RUPTURE_SHELL(7,5) = .00000019550
1356 RUPTURE_SHELL(7,6) = .00000031220
1357! --------------------
1358! Spalling (only for solid)
1359
1360 ! MLN <28 or ==49
1361 RUPTURE_SOL(8,1) = .00000009428
1362 RUPTURE_SOL12(8,1) = .00000025070
1363 RUPTURE_TET4(8,1) = .00000009402
1364 RUPTURE_TET10(8,1) = .00000038110
1365
1366 ! MLN >=28 or /=49
1367 RUPTURE_SOL(8,2) = .00000003871
1368 RUPTURE_SOL12(8,2) = .00000023730
1369 RUPTURE_TET4(8,2) = .00000003716
1370 RUPTURE_TET10(8,2) = .00000014820
1371! --------------------
1372! Wierzbicki
1373 ! MLN <28 or ==49
1374 RUPTURE_SHELL(9,1) = .00000025370
1375 RUPTURE_SHELL(9,2) = .00000071910
1376 RUPTURE_SHELL(9,3) = .00000118100
1377
1378 RUPTURE_SHELL(9,4) = .00000024610
1379 RUPTURE_SHELL(9,5) = .00000071940
1380 RUPTURE_SHELL(9,6) = .00000118500
1381
1382 ! MLN <28 or ==49
1383 RUPTURE_SOL(9,1) = .00000008443
1384 RUPTURE_SOL12(9,1) = .00000011190
1385 RUPTURE_TET4(9,1) = .00000008236
1386 RUPTURE_TET10(9,1) = .00000032530
1387
1388 ! MLN >=28 or /=49
1389 RUPTURE_SOL(9,2) = .00000002591
1390 RUPTURE_SOL12(9,2) = .00000016310
1391 RUPTURE_TET4(9,2) = .00000002483
1392 RUPTURE_TET10(9,2) = .00000008821
1393! --------------------
1394! Tensstrain
1395 ! MLN <28 or ==32
1396 RUPTURE_SHELL(10,1) = .00000009068
1397 RUPTURE_SHELL(10,2) = .00000016430
1398 RUPTURE_SHELL(10,3) = .00000023120
1399
1400 ! MLN >=28 or /=32
1401 RUPTURE_SHELL(10,4) = .00000004473
1402 RUPTURE_SHELL(10,5) = .00000012650
1403 RUPTURE_SHELL(10,6) = .00000019170
1404
1405 ! MLN <28 or ==49
1406 RUPTURE_SOL(10,1) = .00000013990
1407 RUPTURE_SOL12(10,1) = .00000025600
1408 RUPTURE_TET4(10,1) = .00000010960
1409 RUPTURE_TET10(10,1) = .00000047680
1410
1411 ! MLN >=28 or /=49
1412 RUPTURE_SOL(10,2) = .00000008052
1413 RUPTURE_SOL12(10,2) = .00000056320
1414 RUPTURE_TET4(10,2) = .00000005430
1415 RUPTURE_TET10(10,2) = .00000025250
1416! --------------------
1417! Energy
1418 ! MLN <28 or ==49
1419 RUPTURE_SHELL(11,1) = .00000005295
1420 RUPTURE_SHELL(11,2) = .00000011720
1421 RUPTURE_SHELL(11,3) = .00000017760
1422
1423 RUPTURE_SHELL(11,4) = .00000003770
1424 RUPTURE_SHELL(11,5) = .00000010100
1425 RUPTURE_SHELL(11,6) = .00000015840
1426
1427 ! MLN <28 or ==49
1428 RUPTURE_SOL(11,1) = .00000009413
1429 RUPTURE_SOL12(11,1) = .00000009015
1430 RUPTURE_TET4(11,1) = .00000009164
1431 RUPTURE_TET10(11,1) = .00000037960
1432
1433 ! MLN >=28 or /=49
1434 RUPTURE_SOL(11,2) = .00000003883
1435 RUPTURE_SOL12(11,2) = .00000025910
1436 RUPTURE_TET4(11,2) = .00000003753
1437 RUPTURE_TET10(11,2) = .00000014030
1438! --------------------
1439! Chang (only for shell)
1440 ! MLN <28 or ==49
1441 RUPTURE_SHELL(13,1) = .00000004822
1442 RUPTURE_SHELL(13,2) = .00000010740
1443 RUPTURE_SHELL(13,3) = .00000016000
1444 ! MLN >=28 or /=49
1445 RUPTURE_SHELL(13,4) = .00000003770
1446 RUPTURE_SHELL(13,5) = .00000010100
1447 RUPTURE_SHELL(13,6) = .00000015840
1448! --------------------
1449! Hashin
1450 ! MLN <28 or ==49
1451 RUPTURE_SHELL(14,1) = .00000007989
1452 RUPTURE_SHELL(14,2) = .00000025240
1453 RUPTURE_SHELL(14,3) = .00000040300
1454
1455 RUPTURE_SHELL(14,4) = .00000007989
1456 RUPTURE_SHELL(14,5) = .00000025240
1457 RUPTURE_SHELL(14,6) = .00000040300
1458
1459 ! MLN <28 or ==49
1460 RUPTURE_SOL(14,1) = .00000013480
1461 RUPTURE_SOL12(14,1) = .00000107400
1462 RUPTURE_TET4(14,1) = .00000012500
1463 RUPTURE_TET10(14,1) = .00000054770
1464
1465 ! MLN >=28 or /=49
1466 RUPTURE_SOL(14,2) = .00000001289
1467 RUPTURE_SOL12(14,2) = .00000010990
1468 RUPTURE_TET4(14,2) = .00000001302
1469 RUPTURE_TET10(14,2) = .00000005823
1470! --------------------
1471! Puck (compatible MLN<28)
1472 ! MLN <28 or ==49
1473 RUPTURE_SHELL(16,1) = .00000007462
1474 RUPTURE_SHELL(16,2) = .00000018470
1475 RUPTURE_SHELL(16,3) = .00000028840
1476
1477 RUPTURE_SHELL(16,4) = .00000006606
1478 RUPTURE_SHELL(16,5) = .00000017190
1479 RUPTURE_SHELL(16,6) = .00000027600
1480
1481 ! MLN <28 or ==49
1482 RUPTURE_SOL(16,1) = .00000013060
1483 RUPTURE_SOL12(16,1) = .00000008806
1484 RUPTURE_TET4(16,1) = .00000012170
1485 RUPTURE_TET10(16,1) = .00000042400
1486
1487 RUPTURE_SOL(16,2) = .00000005419
1488 RUPTURE_SOL12(16,2) = .00000007646
1489 RUPTURE_TET4(16,2) = .00000001338
1490 RUPTURE_TET10(16,2) = .00000004460
1491! --------------------
1492! Lad_Dama (compatible MLN<28) (only for solid)
1493 ! MLN <28 or ==49
1494 RUPTURE_SOL(18,1) = .00000016890
1495 RUPTURE_SOL12(18,1) = .00000008820
1496 RUPTURE_TET4(18,1) = .00000013980
1497 RUPTURE_TET10(18,1) = .00000058630
1498
1499 RUPTURE_SOL(18,2) = .00000016890
1500 RUPTURE_SOL12(18,2) = .00000008820
1501 RUPTURE_TET4(18,2) = .00000013980
1502 RUPTURE_TET10(18,2) = .00000058630
1503! --------------------
1504! Connect (compatible MLN==59 et 83)
1505 ! MLN==59
1506 RUPTURE_SOL(20,1) = .00000054140
1507 RUPTURE_SOL12(20,1) = .00000054140
1508 RUPTURE_TET4(20,1) = .00000054140
1509 RUPTURE_TET10(20,1) = .00000054140
1510 ! MLN==83
1511 RUPTURE_SOL(20,2) = .00000054140
1512 RUPTURE_SOL12(20,2) = .00000054140
1513 RUPTURE_TET4(20,2) = .00000054140
1514 RUPTURE_TET10(20,2) = .00000054140
1515! --------------------
1516! Tab1
1517 ! MLN <28 or ==49
1518 RUPTURE_SHELL(23,1) = .00000021020
1519 RUPTURE_SHELL(23,2) = .00000060850
1520 RUPTURE_SHELL(23,3) = .00000098850
1521
1522 RUPTURE_SHELL(23,4) = .00000021590
1523 RUPTURE_SHELL(23,5) = .00000062580
1524 RUPTURE_SHELL(23,6) = .0000009949
1525
1526 ! MLN <28 or ==49
1527 RUPTURE_SOL(23,1) = .00000031120
1528 RUPTURE_SOL12(23,1) = .00000192900
1529 RUPTURE_TET4(23,1) = .00000026620
1530 RUPTURE_TET10(23,1) = .00000118100
1531
1532 ! MLN >=28 or /=49
1533 RUPTURE_SOL(23,2) = .00000025940
1534 RUPTURE_SOL12(23,2) = .00000198600
1535 RUPTURE_TET4(23,2) = .00000022330
1536 RUPTURE_TET10(23,2) = .00000096400
1537! --------------------
1538! NXT (only for shell)
1539 ! MLN <28 or ==32
1540 RUPTURE_SHELL(25,1) = .00000014880
1541 RUPTURE_SHELL(25,2) = .00000032860
1542 RUPTURE_SHELL(25,3) = .00000051100
1543
1544 ! MLN >=28 or /=32
1545 RUPTURE_SHELL(25,4) = .00000008550
1546 RUPTURE_SHELL(25,5) = .00000024030
1547 RUPTURE_SHELL(25,6) = .00000038060
1548! --------------------
1549! Snconnect (compatible MLN==59 et 83)
1550 ! MLN==59
1551 RUPTURE_SOL(26,1) = .00000065150
1552 RUPTURE_SOL12(26,1) = .00000065150
1553 RUPTURE_TET4(26,1) = .00000065150
1554 RUPTURE_TET10(26,1) = .00000065150
1555 ! MLN==83
1556 RUPTURE_SOL(26,2) = .00000065150
1557 RUPTURE_SOL12(26,2) = .00000065150
1558 RUPTURE_TET4(26,2) = .00000065150
1559 RUPTURE_TET10(26,2) = .00000065150
1560! --------------------
1561! EMC (only for solid)
1562 ! MLN <28 or ==49
1563 RUPTURE_SOL(27,1) = .00000007774
1564 RUPTURE_SOL12(27,1) = .00000013030
1565 RUPTURE_TET4(27,1) = .00000007539
1566 RUPTURE_TET10(27,1) = .00000029830
1567
1568 ! MLN >=28 or /=49
1569 RUPTURE_SOL(27,2) = .00000001913
1570 RUPTURE_SOL12(27,2) = .00000010660
1571 RUPTURE_TET4(27,2) = .00000001654
1572 RUPTURE_TET10(27,2) = .00000006008
1573! --------------------
1574! Alter (only shell/tri)
1575 ! MLN <28 or ==49
1576 RUPTURE_SHELL(28,1) = .0000000979
1577 RUPTURE_SHELL(28,2) = .0000002061
1578 RUPTURE_SHELL(28,3) = .0000003868
1579
1580 RUPTURE_SHELL(28,4) = .0000000979
1581 RUPTURE_SHELL(28,5) = .0000002061
1582 RUPTURE_SHELL(28,6) = .0000003868
1583! --------------------
1584! Biquad
1585 ! MLN <28 or ==49
1586 RUPTURE_SHELL(30,1) = .00000003069
1587 RUPTURE_SHELL(30,2) = .00000005533
1588 RUPTURE_SHELL(30,3) = .00000008252
1589
1590 RUPTURE_SHELL(30,4) = .00000003133
1591 RUPTURE_SHELL(30,5) = .00000006587
1592 RUPTURE_SHELL(30,6) = .00000010300
1593
1594 ! MLN <28 or ==49
1595 RUPTURE_SOL(30,1) = .00000008026
1596 RUPTURE_SOL12(30,1) = .00000013200
1597 RUPTURE_TET4(30,1) = .00000007523
1598 RUPTURE_TET10(30,1) = .00000030150
1599
1600 RUPTURE_SOL(30,2) = .00000001929
1601 RUPTURE_SOL12(30,2) = .00000011930
1602 RUPTURE_TET4(30,2) = .00000001793
1603 RUPTURE_TET10(30,2) = .00000006588
1604! --------------------
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