OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
weights_p4linuxa964_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 ARM machine
25! Used machine : Cavium --> Cavium ThunderX2(R) CPU CN9975 v2.1 @ 2.20GHz
26
27
28! temps de reference correspond au temps /elt/materiaux loi2 0 pt integ
29 TPSREF = .0000008130 !.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) = .0000002812 !.0000001512
64 SOL8TNL(I,J) = .0000012920 !.0000006297
65 TET4TNL(I,J) = .0000002524 !.0000001542
66 TET10TNL(I,J)= .0000009849 !.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) = .0000002419
76 SOL8TNL(1,1) = .0000007397
77 TET4TNL(1,1) = .0000002074
78 TET10TNL(1,1)= .0000008310
79! LAW 01 + ISMSTR = 10 or 12
80 SOL1TNL(1,2) = .0000003500
81 SOL8TNL(1,2) = .0000027800
82 TET4TNL(1,2) = .0000003500
83 TET10TNL(1,2)= .0000014050
84! LAW 02
85 SOL1TNL(2,1) = .0000002812 !.0000001512
86 SOL8TNL(2,1) = .0000012920 !.0000006297
87 TET4TNL(2,1) = .0000002524 !.0000001542
88 TET10TNL(2,1)= .0000009849 !.0000005545
89! LAW 02 ICC=1
90 SOL1TNL(2,2) = .0000003010 !.0000002051
91 SOL8TNL(2,2) = .0000013050 !.0000010080
92 TET4TNL(2,2) = .0000002700 !.0000001820
93 TET10TNL(2,2)= .0000011980 !.0000007670
94! LAW 03
95 SOL1TNL(3,1) = .0000002932 !.0000002145
96 SOL8TNL(3,1) = .0000022800 !.0000006086
97 TET4TNL(3,1) = .0000002666 !.0000002292
98 TET10TNL(3,1)= .0000010070 !.0000009112
99! LAW 04
100 SOL1TNL(4,1) = .0000004431 !.0000001555
101 SOL8TNL(4,1) = .0000034500 !.0000006050
102 TET4TNL(4,1) = .0000003978 !.0000002547
103 TET10TNL(4,1)= .0000007207 !.0000009900
104! LAW 06
105 SOL1TNL(6,1) = .0000003679 !.0000002483
106 SOL8TNL(6,1) = .0000047100 !.0000006004
107 TET4TNL(6,1) = .0000003344 !.0000002412
108 TET10TNL(6,1)= .0000012400 !.0000007858
109! LAW 10
110 SOL1TNL(10,1) = .0000002950 !.0000002024
111 SOL8TNL(10,1) = .0000022940 !.0000013650
112 TET4TNL(10,1) = .0000002633 !.0000001815
113 TET10TNL(10,1)= .0000010430 !.0000006534
114! LAW 12
115 SOL1TNL(12,1) = .0000005186 !.0000003941
116 SOL8TNL(12,1) = .0000041180 !.0000038700
117 TET4TNL(12,1) = .0000004858 !.0000003670
118 TET10TNL(12,1)= .0000020250 !.0000013560
119! LAW 21
120 SOL1TNL(21,1) = .0000002851 !.0000001872
121 SOL8TNL(21,1) = .0000022340 !.0000013840
122 TET4TNL(21,1) = .0000002537 !.0000001811
123 TET10TNL(21,1)= .0000010050 !.0000006885
124! LAW 22
125 SOL1TNL(22,1) = .0000003647 !.0000001854
126 SOL8TNL(22,1) = .0000028350 !.0000013990
127 TET4TNL(22,1) = .0000003236 !.0000001663
128 TET10TNL(22,1)= .0000013300 !.0000007349
129! LAW 23
130 SOL1TNL(23,1) = .0000003580 !.0000001890
131 SOL8TNL(23,1) = .0000028330 !.0000014360
132 TET4TNL(23,1) = .0000003287 !.0000001850
133 TET10TNL(23,1)= .0000013210 !.0000007381
134! LAW 24
135 SOL1TNL(24,1) = .0000006084 !.0000006082
136 SOL8TNL(24,1) = .0000047450 !.0000042810
137 TET4TNL(24,1) = .0000005787 !.0000004260
138 TET10TNL(24,1)= .0000023220 !.0000020050
139! LAW 25
140 SOL1TNL(25,1) = .0000006250
141 SOL8TNL(25,1) = .0000050290
142 TET4TNL(25,1) = .0000006250
143 TET10TNL(25,1)= .0000025400
144! LAW 25 Crasurv
145 SOL1TNL(25,2) = .0000007571
146 SOL8TNL(25,2) = .0000060500
147 TET4TNL(25,2) = .0000007571
148 TET10TNL(25,2)= .0000030280
149! LAW 28
150 SOL1TNL(28,1) = .0000005104 !.0000003664
151 SOL8TNL(28,1) = .0000035740 !.0000021385
152 TET4TNL(28,1) = .0000005126 !.0000003168
153 TET10TNL(28,1)= .0000020540 !.0000013620
154! LAW 33 0 fct
155 SOL1TNL(33,1) = .0000008186 !.0000004453
156 SOL8TNL(33,1) = .0000062330 !.0000021870
157 TET4TNL(33,1) = .0000012030 !.0000004027
158 TET10TNL(33,1)= .0000043420 !.0000015750
159! LAW 33 1 fct
160 SOL1TNL(33,2) = .0000009964 !.0000006103
161 SOL8TNL(33,2) = .0000076200 !.0000025320
162 TET4TNL(33,2) = .0000013680 !.0000005659
163 TET10TNL(33,2)= .0000050210 !.0000018270
164! LAW 34
165 SOL1TNL(34,1) = .0000003399 !.0000002516
166 SOL8TNL(34,1) = .0000021420 !.0000012240
167 TET4TNL(34,1) = .0000003245 !.0000002258
168 TET10TNL(34,1)= .0000013360 !.0000008617
169! LAW 35
170 SOL1TNL(35,1) = .0000003610 !.0000003070
171 SOL8TNL(35,1) = .0000024910 !.0000017020
172 TET4TNL(35,1) = .0000003610 !.0000002691
173 TET10TNL(35,1)= .0000014590 !.0000010870
174! LAW 36 1 fonction
175 SOL1TNL(36,1) = .0000003810 !.0000003096
176 SOL8TNL(36,1) = .0000025850 !.0000016040
177 TET4TNL(36,1) = .0000003863 !.0000002877
178 TET10TNL(36,1)= .0000015590 !.0000009584
179! LAW 36 5 fonctions
180 SOL1TNL(36,2) = .0000003900 !.0000003358
181 SOL8TNL(36,2) = .0000026000 !.0000016420
182 TET4TNL(36,2) = .0000003978 !.0000002724
183 TET10TNL(36,2)= .0000015870 !.0000010120
184! LAW 36 10 fonctions
185 SOL1TNL(36,3) = .0000003950 !.0000003702
186 SOL8TNL(36,3) = .0000026250 !.0000016950
187 TET4TNL(36,3) = .0000004001 !.0000002854
188 TET10TNL(36,3)= .0000016000 !.0000010970
189! LAW 38
190 SOL1TNL(38,1) = .0000011530
191 SOL8TNL(38,1) = .0000089810
192 TET4TNL(38,1) = .0000011530
193 TET10TNL(38,1)= .0000054800
194! LAW 38 + ISMSTR=10 or 12
195 SOL1TNL(38,2) = .0000012800
196 SOL8TNL(38,2) = .0000120810
197 TET4TNL(38,2) = .0000012800
198 TET10TNL(38,2)= .0000061700
199! LAW 40
200 SOL1TNL(40,1) = .0000007278 !.0000004917
201 SOL8TNL(40,1) = .0000054950 !.0000030250
202 TET4TNL(40,1) = .0000006422 !.0000004015
203 TET10TNL(40,1)= .0000028470 !.0000013700
204! LAW 41
205 SOL1TNL(41,1) = .0000004948 !.0000004092
206 TET4TNL(41,1) = .0000005529 !.0000004068
207 TET10TNL(41,1)= .0000022130 !.0000017520
208! LAW 42
209 SOL1TNL(42,1) = .0000011990 !.0000006159
210 SOL8TNL(42,1) = .0000091810 !.0000041410
211 TET4TNL(42,1) = .0000012690 !.0000006126
212 TET10TNL(42,1)= .0000050100 !.0000023340
213C M = 1
214 SOL1TNL(42,2) = .0000012680 !.0000006763
215 SOL8TNL(42,2) = .0000098450 !.0000045820
216 TET4TNL(42,2) = .0000013590 !.0000006537
217 TET10TNL(42,2)= .0000053250 !.0000023430
218C M = 2
219 SOL1TNL(42,3) = .0000012730 !.0000007634
220 SOL8TNL(42,3) = .0000099550 !.0000051210
221 TET4TNL(42,3) = .0000013730 !.0000006823
222 TET10TNL(42,3)= .0000054540 !.0000027020
223! LAW 44
224 SOL1TNL(44,1) = .0000005450
225 SOL8TNL(44,1) = .0000043000
226 TET4TNL(44,1) = .0000005450
227 TET10TNL(44,1)= .0000025600
228! LAW 48
229 SOL1TNL(48,1) = .0000004502 !.0000002994
230 SOL8TNL(48,1) = .0000032390 !.0000015830
231 TET4TNL(48,1) = .0000004059 !.0000002660
232 TET10TNL(48,1)= .0000018150 !.0000009488
233! LAW 50
234 SOL1TNL(50,1) = .0000005967 !.0000004482
235 SOL8TNL(50,1) = .0000043780 !.0000027100
236 TET4TNL(50,1) = .0000005735 !.0000004178
237 TET10TNL(50,1)= .0000024870 !.0000018130
238! LAW 52
239 SOL1TNL(52,1) = .0000003722 !.0000002826
240 SOL8TNL(52,1) = .0000025360 !.0000015560
241 TET4TNL(52,1) = .0000003521 !.0000002650
242 TET10TNL(52,1)= .0000014420 !.0000010590
243! LAW 53
244 SOL1TNL(53,1) = .0000005613 !.0000004006
245 SOL8TNL(53,1) = .0000039870 !.0000022760
246 TET4TNL(53,1) = .0000005425 !.0000003889
247 TET10TNL(53,1)= .0000021850 !.0000014540
248! LAW 59 (user solid property) IGT=43
249 SOL1TNL(59,1) = .0000002950 !.0000002207
250 SOL8TNL(59,1) = .0000002950 !.0000002207
251 TET4TNL(59,1) = .0000002950 !.0000002207
252 TET10TNL(59,1)= .0000002950 !.0000002207
253! LAW 60
254 SOL1TNL(60,1) = .0000005551 !.0000004239
255 SOL8TNL(60,1) = .0000039320 !.0000024010
256 TET4TNL(60,1) = .0000005112 !.0000003718
257 TET10TNL(60,1)= .0000022000 !.0000014290
258! LAW 62
259 ! M = 0
260 SOL1TNL(62,1) = .0000013480 !.0000009680
261 SOL8TNL(62,1) = .0000105600 !.0000068900
262 TET4TNL(62,1) = .0000013650 !.0000009896
263 TET10TNL(62,1)= .0000055380 !.0000039430
264 ! M = 1
265 SOL1TNL(62,2) = .0000014040 !.0000010850
266 SOL8TNL(62,2) = .0000113000 !.0000078600
267 TET4TNL(62,2) = .0000014380 !.0000010760
268 TET10TNL(62,2)= .0000058710 !.0000043120
269 ! M = 2
270 SOL1TNL(62,3) = .0000014900 !.0000011420
271 SOL8TNL(62,3) = .0000161000 !.0000084600
272 TET4TNL(62,3) = .0000014790 !.0000011150
273 TET10TNL(62,3)= .0000060010 !.0000045435
274! LAW 65
275 SOL1TNL(65,1) = .0000004749 !.0000004282
276 SOL8TNL(65,1) = .0000036940 !.0000030270
277 TET4TNL(65,1) = .0000004395 !.0000003897
278 TET10TNL(65,1)= .0000018470 !.0000015160
279! LAW 66
280 SOL1TNL(66,1) = .0000004490 !.0000005367
281 SOL8TNL(66,1) = .0000035110 !.0000031200
282 TET4TNL(66,1) = .0000004426 !.0000003860
283 TET10TNL(66,1)= .0000017880 !.0000015470
284! LAW 68
285 SOL1TNL(68,1) = .0000008295 !.0000006413
286! SOL8TNL(68,1) = .0000028940
287 TET4TNL(68,1) = .0000008205 !.0000006475
288 TET10TNL(68,1)= .0000032800 !.0000025900
289! LAW 69
290 ! M = 0 (no visc)
291 SOL1TNL(69,1) = .0000010910 !.0000005245
292 SOL8TNL(69,1) = .0000089540 !.0000040350
293 TET4TNL(69,1) = .0000012170 !.0000005135
294 TET10TNL(69,1)= .0000047590 !.0000020020
295 ! M = 1 (visc)
296 SOL1TNL(69,2) = .0000011460 !.0000006120
297 SOL8TNL(69,2) = .0000093160 !.0000048180
298 TET4TNL(69,2) = .0000013480 !.0000006347
299 TET10TNL(69,2)= .0000053070 !.0000024700
300 ! M = 2 (visc)
301 SOL1TNL(69,3) = .0000011860 !.0000006370
302 SOL8TNL(69,3) = .0000095160 !.0000051600
303 TET4TNL(69,3) = .0000013850 !.0000006740
304 TET10TNL(69,3)= .0000054380 !.0000025500
305! LAW 70
306! 0 funct
307 SOL1TNL(70,1) = .0000006092 !.0000003853
308 SOL8TNL(70,1) = .0000052170 !.0000026830
309 TET4TNL(70,1) = .0000005915 !.0000003395
310 TET10TNL(70,1)= .0000023540 !.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) = .0000000006
318 SOL8TNL(70,3) = .0000000048
319 TET4TNL(70,3) = .0000000005
320 TET10TNL(70,3)= .0000000020
321! LAW 71
322 SOL1TNL(71,1) = .0000012840 !.0000006332
323 SOL8TNL(71,1) = .0000059560 !.0000021590
324 TET4TNL(71,1) = .0000012720 !.0000004840
325 TET10TNL(71,1)= .0000046600 !.0000018860
326! LAW 72
327 SOL1TNL(72,1) = .0000006166 !.0000004335
328 SOL8TNL(72,1) = .0000048530 !.0000032810
329 TET4TNL(72,1) = .0000004856 !.0000003947
330 TET10TNL(72,1)= .0000020700 !.0000016170
331! LAW 74
332 SOL1TNL(74,1) = .0000004670 !.0000003704
333 SOL8TNL(74,1) = .0000037240 !.0000027610
334 TET4TNL(74,1) = .0000004741 !.0000003402
335 TET10TNL(74,1)= .0000019010 !.0000013590
336! LAW 76
337 SOL1TNL(76,1) = .0000003535 !.0000002745
338 SOL8TNL(76,1) = .0000027540 !.0000021280
339 TET4TNL(76,1) = .0000003369 !.0000002739
340 TET10TNL(76,1)= .0000013940 !.0000011500
341! LAW 77
342 SOL1TNL(77,1) = .0000006680 !.0000004887
343 SOL8TNL(77,1) = .0000052100 !.0000038100
344 TET4TNL(77,1) = .0000006640 !.0000005200
345 TET10TNL(77,1)= .0000026230 !.0000020120
346! LAW 78
347 SOL1TNL(78,1) = .0000003481 !.0000003087
348 SOL8TNL(78,1) = .0000027990 !.0000023090
349 TET4TNL(78,1) = .0000003558 !.0000002902
350 TET10TNL(78,1)= .0000014310 !.0000012520
351! LAW 79
352 SOL1TNL(79,1) = .0000003856 !.0000003060
353 SOL8TNL(79,1) = .0000029650 !.0000020220
354 TET4TNL(79,1) = .0000003752 !.0000002684
355 TET10TNL(79,1)= .0000015180 !.0000011400
356! LAW 81
357 SOL1TNL(81,1) = .0000004321 !.0000004223
358 SOL8TNL(81,1) = .0000033820 !.0000031430
359 TET4TNL(81,1) = .0000003985 !.0000002750
360 TET10TNL(81,1)= .0000016420 !.0000012470
361! LAW 82
362 ! M = 1
363 SOL1TNL(82,1) = .0000012400
364 SOL8TNL(82,1) = .0000109000
365 TET4TNL(82,1) = .0000012300
366 TET10TNL(82,1)= .0000051000
367 ! M = 2
368 SOL1TNL(82,2) = .0000015600
369 SOL8TNL(82,2) = .0000120000
370 TET4TNL(82,2) = .0000015500
371 TET10TNL(82,2)= .0000060300
372 ! M = 3
373 SOL1TNL(82,3) = .0000019700
374 SOL8TNL(82,3) = .0000071000
375 TET4TNL(82,3) = .0000019400
376 TET10TNL(82,3)= .0000078800
377! LAW 83 (user solid property) IGT=43
378 SOL1TNL(83,1) = .0000008830 !.0000013690
379 SOL8TNL(83,1) = .0000008830 !.0000013690
380 TET4TNL(83,1) = .0000008830 !.0000013690
381 TET10TNL(83,1)= .0000008830 !.0000013690
382! LAW 84
383 SOL1TNL(84,1) = .0000002868 !.0000002771
384 SOL8TNL(84,1) = .0000022330 !.0000018140
385 TET4TNL(84,1) = .0000002866 !.0000002489
386 TET10TNL(84,1)= .0000011630 !.0000007861
387! LAW 88
388 SOL1TNL(88,1) = .0000049000
389 SOL8TNL(88,1) = .0000379000
390 TET4TNL(88,1) = .0000048190
391 TET10TNL(88,1)= .0000192760
392! LAW 90
393 SOL1TNL(90,1) = .0000011200
394 SOL8TNL(90,1) = .0000089600
395 TET4TNL(90,1) = .0000011200
396 TET10TNL(90,1)= .0000044800
397! LAW 90 + ISMSTR=10 or 12
398 SOL1TNL(90,2) = .0000014500
399 SOL8TNL(90,2) = .0000100400
400 TET4TNL(90,2) = .0000014500
401 TET10TNL(90,2)= .0000052700
402! LAW 92
403 SOL1TNL(92,1) = .0000013000
404 SOL8TNL(92,1) = .0000099000
405 TET4TNL(92,1) = .0000013000
406 TET10TNL(92,1)= .0000047000
407! LAW 92 + ISMSTR=10 or 12
408 SOL1TNL(92,2) = .0000014500
409 SOL8TNL(92,2) = .0000112000
410 TET4TNL(92,2) = .0000014500
411 TET10TNL(92,2)= .0000056400
412! LAW 94
413 SOL1TNL(94,1) = .0000009800
414 SOL8TNL(94,1) = .0000078300
415 TET4TNL(94,1) = .0000009800
416 TET10TNL(94,1)= .0000038900
417! LAW 94 + ISMSTR=10 or 12
418 SOL1TNL(94,2) = .0000010900
419 SOL8TNL(94,2) = .0000089400
420 TET4TNL(94,2) = .0000010900
421 TET10TNL(94,2)= .0000043500
422! LAW 95
423 SOL1TNL(95,1) = .0000008605 !.0000005950
424 SOL8TNL(95,1) = .0000069280 !.0000045150
425 TET4TNL(95,1) = .0000008297 !.0000005905
426 TET10TNL(95,1)= .0000033970 !.0000021230
427
428! LAW 100 : material cost = base cost + viscosity cost + N * network cost
429! LAW 100 (base)
430 SOL1TNL(100,1) = .0000002280
431 SOL8TNL(100,1) = .0000002280
432 TET4TNL(100,1) = .0000002280
433 TET10TNL(100,1)= .0000002280
434
435! LAW 100 (viscosity)
436 SOL1TNL(100,2) = .0000001390
437 SOL8TNL(100,2) = .0000001390
438 TET4TNL(100,2) = .0000001390
439 TET10TNL(100,2)= .0000001390
440
441! LAW 100 (network)
442 SOL1TNL(100,3) = .0000002000
443 SOL8TNL(100,3) = .0000002000
444 TET4TNL(100,3) = .0000002000
445 TET10TNL(100,3)= .0000002000
446! LAW 101
447 SOL1TNL(101,1) = .0000102200
448 SOL8TNL(101,1) = .0000852000
449 TET4TNL(101,1) = .0000102200
450 TET10TNL(101,1)= .0000421000
451! LAW 102
452 SOL1TNL(102,1) = .0000003368
453 SOL8TNL(102,1) = .0000026200
454 TET4TNL(102,1) = .0000003100
455 TET10TNL(102,1)= .0000011700
456! LAW 103
457 SOL1TNL(103,1) = .0000003760
458 SOL8TNL(103,1) = .0000029060
459 TET4TNL(103,1) = .0000003590
460 TET10TNL(103,1)= .0000013640
461! LAW 104
462 ! Nice algo
463 SOL1TNL(104,1) = .0000002691
464 SOL8TNL(104,1) = .0000002691
465 TET4TNL(104,1) = .0000002691
466 TET10TNL(104,1) = .0000002691
467 ! Newton algo
468 SOL1TNL(104,2) = .0000003400
469 SOL8TNL(104,2) = .0000003400
470 TET4TNL(104,2) = .0000003400
471 TET10TNL(104,2) = .0000003400
472 ! iloc=1
473 SOL_OPTION(1,1) = .0000000800
474 ! iloc=2
475 SOL_OPTION(1,2) = .0000001000
476 ! iloc=3
477 SOL_OPTION(1,3) = .0000002000
478! LAW 106
479! SOL1TNL(106,1) = .0000002493
480! SOL8TNL(106,1) = .0000018090
481! TET4TNL(106,1) = .0000002246
482! TET10TNL(106,1)= .0000009194
483! LAW 187
484 SOL1TNL(187,1) = .0000009595
485 SOL8TNL(187,1) = .0000059000
486 TET4TNL(187,1) = .0000008900
487 TET10TNL(187,1)= .0000033600
488
489 IF (IPARI0.EQ.1) THEN
490! SOLTELT :
491! 1 --> ISOLID = 1
492! 2 --> ISOLID = 2
493! 3 --> ISOLID = 24
494! 4 --> ISOLID = 12
495! 5 --> ISOLID = 14 HA8 with NPT=2*2*2 (fixed part of element cost)
496! 6 --> ISOLID = 14 HA8 (overcost for NPT > 2*2*2)
497! HA8 cost : COST = HA8 (222) cost + (NPT-8) * overcost
498! 7 --> ISOLID = 17 H8C
499! 8 --> IGTYP = 43
500! 9 --> ISOLID = 18 H8C
501! 10--> IGTYP = Thickshell 20 or 21 or 22
502! 11--> JHBE=15 : fixed part of element cost
503! 12--> JHBE=15 : overcost for N>1
504!
505! NLOCAL_OPTION
506! 1 --> sol
507! 2 --> shell
508!
509! TET4ELT :
510! 1 --> TETRA 4
511! 2 --> TETRA 4 with IROT = 1
512! TETRA10ELT
513! 1 --> TETRA10
514!
515! ELM_OVER_COST :
516! 1 --> ismstr=10 or 12
517
518 SOLTELT(1) =.0000009427 !.0000009528
519 SOLTELT(2) =.0000007901 !.0000008540
520 SOLTELT(3) =.0000013950 !.0000013420
521 SOLTELT(4) =.0000020790 !.0000015180
522 SOLTELT(5) =.0000060280 !.0000041200
523 SOLTELT(6) =.0000006515 !.0000003812
524 SOLTELT(7) =.0000034520 !.0000026600
525 SOLTELT(8) =.0000013670 !.0000012750
526 SOLTELT(9) =.0000032000
527 SOLTELT(10) =.0000005800
528 SOLTELT(11) =.0000017250
529 SOLTELT(12) =.0000001389
530 TET4TELT(1) =.0000003669 !.0000004744
531 TET4TELT(2) =.0000021890 !.0000016460
532 TET10TELT(1)=.0000021640 !.0000018340
533 NLOCAL_OPTION(1) = .0000002500
534
535 ELM_OVER_COST(1) = .0000001800
536 ELSE
537 SOLTELT(1) =.0000009361 !.0000006901
538 SOLTELT(2) =.0000007850 !.0000005650
539 SOLTELT(3) =.0000013970 !.0000011240
540 SOLTELT(4) =.0000020790 !.0000013490
541 SOLTELT(5) =.0000060280 !.0000036580
542 SOLTELT(6) =.0000006515 !.0000004100
543 SOLTELT(7) =.0000034520 !.0000022490
544 SOLTELT(8) =.0000013670 !.0000009907
545 SOLTELT(9) =.0000025000
546 SOLTELT(10) =.0000005000
547 SOLTELT(11) =.0000015000
548 SOLTELT(12) =.0000001100
549 TET4TELT(1) =.0000003220 !.0000002688
550 TET4TELT(2) =.0000021660 !.0000014960
551 TET10TELT(1)=.0000021130 !.0000016200
552 NLOCAL_OPTION(1) = .0000002500
553
554 ELM_OVER_COST(1) = .0000001500
555 ENDIF
556
557! -------------
558! Coques
559! -------------
560
561
562
563! Tableau SHTNL
564! ---------------
565!
566! Temps par elt/cycle pour les materieux coques
567! SHTNL(I,J,K)=
568! 1er Champ : numero du materiaux
569! 2 champs : nbre de pts d integration, de taille 0:4
570! 4 sont consideres :
571! * 0 : 0pt d integ
572! * 1 : 1pt d integ
573! * 2 : 3pt d integ
574! * 3 : 5pt d integ
575! Les autres temps pour les NPT, sont calculees par interpolation avec ces 3
576! valeurs.
577!
578! 3e champ : option particuliere au besoin
579!
580! -------------
581! Triangles
582! -------------
583
584! Tableau TRITNL
585! ---------------
586!
587! Temps par elt/cycle pour les materieux coques
588! TRITNL(I,J,K)=
589! 1er Champ : numero du materiaux
590! 2e champ : nbre de pts d integration, de taille 0:4
591! 4 sont consideres :
592! * 0 : 0pt d integ
593! * 1 : 1pt d integ
594! * 2 : 3pt d integ
595! * 3 : 5pt d integ
596! Les autres temps pour les NPT, sont calculees par interpolation avec ces 3
597! valeurs.
598!
599! 3e champ : option particuliere au besoin
600!
601!
602! Tableau TRITELT
603! ---------------
604!
605! Temps par elt/cycle pour le traitement de l hourglas en fct de /IPARITH
606
607 DO I=1,3
608 DO J = 0, MAXLAW
609 SHTNL(J,0,I) = .0000002210
610 SHTNL(J,1,I) = .0000002328
611 SHTNL(J,2,I) = .0000004351
612 SHTNL(J,3,I) = .0000006344
613C
614 TRITNL(J,0,I)= .0000002259
615 TRITNL(J,1,I)= .0000002230
616 TRITNL(J,2,I)= .0000004896
617 TRITNL(J,3,I)= .0000006200
618 ENDDO
619 ENDDO
620! special option for shell or shell3n
621 SHTNL_OPTION(1:10,1:5) = 0.
622
623! LAW 01
624 SHTNL(1,0,1) = .0000001083 !.0000000897
625 SHTNL(1,1,1) = .0000001065 !.0000000811
626 SHTNL(1,2,1) = .0000001154 !.0000001004
627 SHTNL(1,3,1) = .0000001465 !.0000001032
628C
629 TRITNL(1,0,1)= .0000001072 !.0000000999
630 TRITNL(1,1,1)= .0000001047 !.0000000797
631 TRITNL(1,2,1)= .0000001157 !.0000001024
632 TRITNL(1,3,1)= .0000001367 !.0000000999
633! LAW 02
634 SHTNL(2,0,1) = .0000002210
635 SHTNL(2,1,1) = .0000002328
636 SHTNL(2,2,1) = .0000004351
637 SHTNL(2,3,1) = .0000006344
638C
639 TRITNL(2,0,1)= .0000002259
640 TRITNL(2,1,1)= .0000002230
641 TRITNL(2,2,1)= .0000004896
642 TRITNL(2,3,1)= .0000006200
643! LAW 02 CC
644 SHTNL(2,0,2) = .0000003308
645 SHTNL(2,1,2) = .0000003314
646 SHTNL(2,2,2) = .0000006472
647 SHTNL(2,3,2) = .0000009295
648C
649 TRITNL(2,0,2)= .0000003358
650 TRITNL(2,1,2)= .0000003388
651 TRITNL(2,2,2)= .0000006350
652 TRITNL(2,3,2)= .0000009150
653! LAW 15
654 SHTNL(15,0,1) = .0000003704
655 SHTNL(15,1,1) = .0000003704
656 SHTNL(15,2,1) = .0000008502
657 SHTNL(15,3,1) = .0000013410
658C
659 TRITNL(15,0,2)= .0000003694
660 TRITNL(15,1,2)= .0000003694
661 TRITNL(15,2,2)= .0000008633
662 TRITNL(15,3,2)= .0000013700
663
664! LAW 19
665 SHTNL(19,0,1) = .0000002629
666 SHTNL(19,1,1) = .0000002629
667 SHTNL(19,2,1) = .0000005877
668 SHTNL(19,3,1) = .0000008430
669C
670 TRITNL(19,0,1)= .0000002600
671 TRITNL(19,1,1)= .0000002600
672 TRITNL(19,2,1)= .0000005621
673 TRITNL(19,3,1)= .0000008520
674! LAW 22
675 SHTNL(22,0,1) = .0000002741
676 SHTNL(22,1,1) = .0000003629
677 SHTNL(22,2,1) = .0000007896
678 SHTNL(22,3,1) = .0000012290
679C
680 TRITNL(22,0,1)= .0000003869
681 TRITNL(22,1,1)= .0000003869
682 TRITNL(22,2,1)= .0000008346
683 TRITNL(22,3,1)= .0000013230
684! LAW 25
685 SHTNL(25,0,1) = .0000004012 !.0000002511
686 SHTNL(25,1,1) = .0000003969 !.0000002509
687 SHTNL(25,2,1) = .0000009354 !.0000005094
688 SHTNL(25,3,1) = .0000014900 !.0000008026
689C
690 TRITNL(25,0,1)= .0000004072 !.0000002519
691 TRITNL(25,1,1)= .0000003967 !.0000002455
692 TRITNL(25,2,1)= .0000009554 !.0000005096
693 TRITNL(25,3,1)= .0000014930 !.0000007655
694! LAW 25 Crasurv (/TYPE11)
695 SHTNL(25,0,2) = .0000005295 !.0000003794
696 SHTNL(25,1,2) = .0000005275 !.0000003794
697 SHTNL(25,2,2) = .0000012420 !.0000009519
698 SHTNL(25,3,2) = .0000019600 !.0000014590
699C
700 TRITNL(25,0,2)= .0000004956 !.0000003866
701 TRITNL(25,1,2)= .0000004925 !.0000003786
702 TRITNL(25,2,2)= .0000012270 !.0000009420
703 TRITNL(25,3,2)= .0000019650 !.0000013980
704! LAW 27 (0pt incompatible)
705! IPLA = 0
706 SHTNL(27,0,1) = .0000004609
707 SHTNL(27,1,1) = .0000004609
708 SHTNL(27,2,1) = .0000011290
709 SHTNL(27,3,1) = .0000018040
710C
711 TRITNL(27,0,1)= .0000004384
712 TRITNL(27,1,1)= .0000004384
713 TRITNL(27,2,1)= .0000010750
714 TRITNL(27,3,1)= .0000016840
715
716! LAW 32
717 SHTNL(32,0,1) = .0000003927
718 SHTNL(32,1,1) = .0000003927
719 SHTNL(32,2,1) = .0000009378
720 SHTNL(32,3,1) = .0000014870
721C
722 TRITNL(32,0,1)= .0000003198
723 TRITNL(32,1,1)= .0000003198
724 TRITNL(32,2,1)= .0000007151
725 TRITNL(32,3,1)= .0000011200
726! LAW 35 (0pt incompatible)
727 SHTNL(35,0,1) = .0000003552 !.0000002696
728 SHTNL(35,1,1) = .0000003552 !.0000002696
729 SHTNL(35,2,1) = .0000007749 !.0000005399
730 SHTNL(35,3,1) = .0000013190 !.0000008694
731C
732 TRITNL(35,0,1)= .0000003552 !.0000002629
733 TRITNL(35,1,1)= .0000003552 !.0000002629
734 TRITNL(35,2,1)= .0000007749 !.0000005365
735 TRITNL(35,3,1)= .0000013190 !.0000009575
736! LAW 36 - 1 fct
737 SHTNL(36,0,1) = .0000003031 !.0000002453
738 SHTNL(36,1,1) = .0000003591 !.0000002975
739 SHTNL(36,2,1) = .0000007585 !.0000005370
740 SHTNL(36,3,1) = .0000011620 !.0000007957
741C
742 TRITNL(36,0,1)= .0000003023 !.0000002315
743 TRITNL(36,1,1)= .0000003646 !.0000003005
744 TRITNL(36,2,1)= .0000007627 !.0000005141
745 TRITNL(36,3,1)= .0000011640 !.0000008021
746! LAW 36 - 5 fct
747 SHTNL(36,0,2) = .0000003160 !.0000002832
748 SHTNL(36,1,2) = .0000003657 !.0000003655
749 SHTNL(36,2,2) = .0000007700 !.0000005902
750 SHTNL(36,3,2) = .0000011800 !.0000008500
751C
752 TRITNL(36,0,2)= .0000003052 !.0000002975
753 TRITNL(36,1,2)= .0000003628 !.0000003608
754 TRITNL(36,2,2)= .0000007831 !.0000005945
755 TRITNL(36,3,2)= .0000011790 !.0000008260
756! LAW 36 - 10 fct
757 SHTNL(36,0,3) = .0000003205 !.0000003180
758 SHTNL(36,1,3) = .0000003701 !.0000003850
759 SHTNL(36,2,3) = .0000007750 !.0000006500
760 SHTNL(36,3,3) = .0000011850 !.0000009020
761C
762 TRITNL(36,0,3)= .0000003087 !.0000003060
763 TRITNL(36,1,3)= .0000003687 !.0000003820
764 TRITNL(36,2,3)= .0000008116 !.0000006540
765 TRITNL(36,3,3)= .0000012330 !.0000009000
766! LAW 42
767 SHTNL(42,0,1) = .0000007810 !.0000005524
768 SHTNL(42,1,1) = .0000007810 !.0000005524
769 SHTNL(42,2,1) = .0000020610 !.0000013820
770 SHTNL(42,3,1) = .0000035030 !.0000022760
771! LAW 42 M = 1
772 SHTNL(42,0,2) = .0000017890 !.0000015100
773 SHTNL(42,1,2) = .0000017890 !.0000015100
774 SHTNL(42,2,2) = .0000048870 !.0000045280
775 SHTNL(42,3,2) = .0000099100 !.0000074640
776! LAW 42 M = 2
777 SHTNL(42,0,3) = .0000021680 !.0000017720
778 SHTNL(42,1,3) = .0000021680 !.0000017720
779 SHTNL(42,2,3) = .0000062360 !.0000053370
780 SHTNL(42,3,3) = .0000104900 !.0000088600
781! LAW 42
782 TRITNL(42,0,1)= .0000008042 !.0000005586
783 TRITNL(42,1,1)= .0000008042 !.0000005586
784 TRITNL(42,2,1)= .0000021190 !.0000015540
785 TRITNL(42,3,1)= .0000034490 !.0000022710
786! LAW 42 M = 1
787 TRITNL(42,0,2)= .0000017160 !.0000015000
788 TRITNL(42,1,2)= .0000017160 !.0000015000
789 TRITNL(42,2,2)= .0000047870 !.0000044580
790 TRITNL(42,3,2)= .0000099040 !.0000073400
791! LAW 42 M = 2
792 TRITNL(42,0,3)= .0000021940 !.0000017200
793 TRITNL(42,1,3)= .0000021940 !.0000017200
794 TRITNL(42,2,3)= .0000062930 !.0000052700
795 TRITNL(42,3,3)= .0000104200 !.0000087000
796! LAW 43
797 SHTNL(43,0,1) = .0000005895 !.0000002796
798 SHTNL(43,1,1) = .0000005865 !.0000002918
799 SHTNL(43,2,1) = .0000008782 !.0000006119
800 SHTNL(43,3,1) = .0000013660 !.0000009136
801C
802 TRITNL(43,0,1)= .0000003910 !.0000002856
803 TRITNL(43,1,1)= .0000003925 !.0000002743
804 TRITNL(43,2,1)= .0000008619 !.0000006279
805 TRITNL(43,3,1)= .0000013340 !.0000008749
806! LAW 44
807 SHTNL(44,0,1) = .0000003952 !.0000002640
808 SHTNL(44,1,1) = .0000003952 !.0000002640
809 SHTNL(44,2,1) = .0000008956 !.0000005714
810 SHTNL(44,3,1) = .0000012950 !.0000008505
811C
812 TRITNL(44,0,1)= .0000003687 !.0000002802
813 TRITNL(44,1,1)= .0000003687 !.0000002802
814 TRITNL(44,2,1)= .0000007813 !.0000005076
815 TRITNL(44,3,1)= .0000011960 !.0000007598
816! LAW 48
817 SHTNL(48,0,1) = .0000003756 !.0000002811
818 SHTNL(48,1,1) = .0000003756 !.0000002811
819 SHTNL(48,2,1) = .0000008432 !.0000005285
820 SHTNL(48,3,1) = .0000012950 !.0000007925
821C
822 TRITNL(48,0,1)= .0000003634 !.0000002764
823 TRITNL(48,1,1)= .0000003634 !.0000002764
824 TRITNL(48,2,1)= .0000007780 !.0000005447
825 TRITNL(48,3,1)= .0000012040 !.0000007677
826! LAW 52
827 SHTNL(52,0,1) = .0000003163 !.0000002313
828 SHTNL(52,1,1) = .0000003163 !.0000002313
829 SHTNL(52,2,1) = .0000006817 !.0000004775
830 SHTNL(52,3,1) = .0000010500 !.0000007161
831C
832 TRITNL(52,0,1)= .0000003249 !.0000002422
833 TRITNL(52,1,1)= .0000003249 !.0000002422
834 TRITNL(52,2,1)= .0000006726 !.0000004612
835 TRITNL(52,3,1)= .0000010100 !.0000007237
836! LAW 57
837 SHTNL(57,0,1) = .0000005852 !.0000005039
838 SHTNL(57,1,1) = .0000005852 !.0000005039
839 SHTNL(57,2,1) = .0000014180 !.0000013330
840 SHTNL(57,3,1) = .0000022860 !.0000021130
841C
842 TRITNL(57,0,1)= .0000005813 !.0000005568
843 TRITNL(57,1,1)= .0000005813 !.0000005568
844 TRITNL(57,2,1)= .0000014150 !.0000014330
845 TRITNL(57,3,1)= .0000022800 !.0000022840
846! LAW 58
847! No visc (FLAGUL==0)
848 SHTNL(58,0,1) = .0000007351 !.0000005435
849 SHTNL(58,1,1) = .0000007336 !.0000006044
850 SHTNL(58,2,1) = .0000021520 !.0000017910
851 SHTNL(58,3,1) = .0000035180 !.0000028510
852C
853 TRITNL(58,0,1)= .0000006849 !.0000005733
854 TRITNL(58,1,1)= .0000006859 !.0000005733
855 TRITNL(58,2,1)= .0000019760 !.0000015320
856 TRITNL(58,3,1)= .0000031160 !.0000023640
857! LAW 60
858 SHTNL(60,0,1) = .0000004346 !.0000003498
859 SHTNL(60,1,1) = .0000005408 !.0000003992
860 SHTNL(60,2,1) = .0000013020 !.0000008897
861 SHTNL(60,3,1) = .0000020480 !.0000013350
862!
863 TRITNL(60,0,1)= .0000005665 !.0000003850
864 TRITNL(60,1,1)= .0000005625 !.0000003997
865 TRITNL(60,2,1)= .0000013970 !.0000009103
866 TRITNL(60,3,1)= .0000022120 !.0000014110
867! LAW 62
868 ! M = 0
869 SHTNL(62,0,1) = .0000013230 !.0000008089
870 SHTNL(62,1,1) = .0000013230 !.0000008089
871 SHTNL(62,2,1) = .0000037110 !.0000021360
872 SHTNL(62,3,1) = .0000060800 !.0000034950
873 ! M = 1
874 SHTNL(62,0,2) = .0000024050 !.0000014640
875 SHTNL(62,1,2) = .0000024050 !.0000014640
876 SHTNL(62,2,2) = .0000068930 !.0000041200
877 SHTNL(62,3,2) = .0000114300 !.0000067100
878 ! M = 2
879 SHTNL(62,0,3) = .0000027790 !.0000017390
880 SHTNL(62,1,3) = .0000027790 !.0000017390
881 SHTNL(62,2,3) = .0000080810 !.0000048970
882 SHTNL(62,3,3) = .0000133600 !.0000079110
883C
884 ! M = 0
885 TRITNL(62,0,1)= .0000013400 !.0000008167
886 TRITNL(62,1,1)= .0000013400 !.0000008167
887 TRITNL(62,2,1)= .0000037250 !.0000021660
888 TRITNL(62,3,1)= .0000061040 !.0000034760
889 ! M = 1
890 TRITNL(62,0,1)= .0000023880 !.0000014950
891 TRITNL(62,1,1)= .0000023880 !.0000014950
892 TRITNL(62,2,1)= .0000068910 !.0000040800
893 TRITNL(62,3,1)= .0000114200 !.0000067340
894 ! M = 2
895 TRITNL(62,0,1)= .0000027850 !.0000017140
896 TRITNL(62,1,1)= .0000027850 !.0000017140
897 TRITNL(62,2,1)= .0000080740 !.0000048680
898 TRITNL(62,3,1)= .0000133400 !.0000079480
899! LAW 63
900 SHTNL(63,0,1) = .0000004878 !.0000002904
901 SHTNL(63,1,1) = .0000004878 !.0000002904
902 SHTNL(63,2,1) = .0000011750 !.0000006582
903 SHTNL(63,3,1) = .0000018530 !.0000009983
904C
905 TRITNL(63,0,1)= .0000004602 !.0000002766
906 TRITNL(63,1,1)= .0000004602 !.0000002766
907 TRITNL(63,2,1)= .0000010660 !.0000005871
908 TRITNL(63,3,1)= .0000016880 !.0000008944
909! LAW 64
910 SHTNL(64,0,1) = .0000004179 !.0000003364
911 SHTNL(64,1,1) = .0000004179 !.0000003364
912 SHTNL(64,2,1) = .0000009481 !.0000007866
913 SHTNL(64,3,1) = .0000014900 !.0000011920
914C
915 TRITNL(64,0,1)= .0000003821 !.0000003426
916 TRITNL(64,1,1)= .0000003821 !.0000003426
917 TRITNL(64,2,1)= .0000008718 !.0000007550
918 TRITNL(64,3,1)= .0000013280 !.0000011970
919! LAW 65
920 SHTNL(65,0,1) = .0000004278 !.0000003350
921 SHTNL(65,1,1) = .0000004278 !.0000003350
922 SHTNL(65,2,1) = .0000009702 !.0000007280
923 SHTNL(65,3,1) = .0000014950 !.0000011240
924C
925 TRITNL(65,0,1)= .0000004232 !.0000003324
926 TRITNL(65,1,1)= .0000004232 !.0000003324
927 TRITNL(65,2,1)= .0000009631 !.0000007364
928 TRITNL(65,3,1)= .0000014680 !.0000011200
929! LAW 66
930 SHTNL(66,0,1) = .0000004841 !.0000003310
931 SHTNL(66,1,1) = .0000004841 !.0000003310
932 SHTNL(66,2,1) = .0000010910 !.0000007200
933 SHTNL(66,3,1) = .0000017330 !.0000011010
934C
935 TRITNL(66,0,1)= .0000004288 !.0000003725
936 TRITNL(66,1,1)= .0000004288 !.0000003666
937 TRITNL(66,2,1)= .0000010020 !.0000007924
938 TRITNL(66,3,1)= .0000015510 !.0000011650
939! LAW 69
940 ! M = 0 (no visc)
941 SHTNL(69,0,1) = .0000017600 !.0000013330
942 SHTNL(69,1,1) = .0000017600 !.0000013330
943 SHTNL(69,2,1) = .0000050370 !.0000032970
944 SHTNL(69,3,1) = .0000083550 !.0000054280
945 ! M = 1 (visc)
946 SHTNL(69,0,2) = .0000018210 !.0000014210
947 SHTNL(69,1,2) = .0000018210 !.0000014210
948 SHTNL(69,2,2) = .0000051910 !.0000039620
949 SHTNL(69,3,2) = .0000086440 !.0000064560
950 ! M = 2 (visc)
951 SHTNL(69,0,3) = .0000018810 !.0000016580
952 SHTNL(69,1,3) = .0000018810 !.0000016580
953 SHTNL(69,2,3) = .0000053700 !.0000045270
954 SHTNL(69,3,3) = .0000089440 !.0000074700
955C
956 ! M = 0 (no visc)
957 TRITNL(69,0,1)= .0000017730 !.0000010400
958 TRITNL(69,1,1)= .0000017730 !.0000010400
959 TRITNL(69,2,1)= .0000050500 !.0000027280
960 TRITNL(69,3,1)= .0000083560 !.0000044630
961 ! M = 1
962 TRITNL(69,0,2)= .0000018310 !.0000012490
963 TRITNL(69,1,2)= .0000018310 !.0000012490
964 TRITNL(69,2,2)= .0000052510 !.0000035000
965 TRITNL(69,3,2)= .0000086270 !.0000055580
966 ! M = 2
967 TRITNL(69,0,3)= .0000018820 !.0000014430
968 TRITNL(69,1,3)= .0000018820 !.0000014430
969 TRITNL(69,2,3)= .0000053890 !.0000040510
970 TRITNL(69,3,3)= .0000088960 !.0000065830
971! LAW 71 (incompatible npt=0)
972 SHTNL(71,0,1) = .0000006504 !.0000005878
973 SHTNL(71,1,1) = .0000006504 !.0000005878
974 SHTNL(71,2,1) = .0000016610 !.0000014690
975 SHTNL(71,3,1) = .0000026980 !.0000024250
976C
977 TRITNL(71,0,1)= .0000006568 !.0000006319
978 TRITNL(71,1,1)= .0000006568 !.0000006919
979 TRITNL(71,2,1)= .0000016750 !.0000015420
980 TRITNL(71,3,1)= .0000027120 !.0000024620
981! LAW 72
982 SHTNL(72,0,1) = .0000005550 !.0000003671
983 SHTNL(72,1,1) = .0000005550 !.0000003671
984 SHTNL(72,2,1) = .0000013380 !.0000007838
985 SHTNL(72,3,1) = .0000020770 !.0000012730
986C
987 TRITNL(72,0,1)= .0000005225 !.0000003591
988 TRITNL(72,1,1)= .0000005225 !.0000003591
989 TRITNL(72,2,1)= .0000012400 !.0000007633
990 TRITNL(72,3,1)= .0000019560 !.0000011460
991! LAW 73
992 SHTNL(73,0,1) = .0000003868 !.0000003456
993 SHTNL(73,1,1) = .0000003868 !.0000003689
994 SHTNL(73,2,1) = .0000008484 !.0000007752
995 SHTNL(73,3,1) = .0000013300 !.0000012010
996C
997 TRITNL(73,0,1)= .0000003816 !.0000003618
998 TRITNL(73,1,1)= .0000003816 !.0000003371
999 TRITNL(73,2,1)= .0000008511 !.0000007637
1000 TRITNL(73,3,1)= .0000013100 !.0000012300
1001! LAW 76
1002 SHTNL(76,0,1) = .0000006827 !.0000007383
1003 SHTNL(76,1,1) = .0000006827 !.0000007476
1004 SHTNL(76,2,1) = .0000017400 !.0000018800
1005 SHTNL(76,3,1) = .0000027800 !.0000029980
1006C
1007 TRITNL(76,0,1)= .0000006633 !.0000007008
1008 TRITNL(76,1,1)= .0000006633 !.0000007008
1009 TRITNL(76,2,1)= .0000017240 !.0000018360
1010 TRITNL(76,3,1)= .0000027750 !.0000029000
1011! LAW 78
1012 SHTNL(78,0,1) = .0000003194 !.0000002921
1013 SHTNL(78,1,1) = .0000003194 !.0000002921
1014 SHTNL(78,2,1) = .0000006581 !.0000005980
1015 SHTNL(78,3,1) = .0000009861 !.0000009682
1016C
1017 TRITNL(78,0,1)= .0000003151 !.0000002936
1018 TRITNL(78,1,1)= .0000003151 !.0000002939
1019 TRITNL(78,2,1)= .0000006536 !.0000005670
1020 TRITNL(78,3,1)= .0000010030 !.0000008430
1021! LAW 80
1022 SHTNL(80,0,1) = .0000015110
1023 SHTNL(80,1,1) = .0000014400
1024 SHTNL(80,2,1) = .0000040060
1025 SHTNL(80,3,1) = .0000065840
1026C
1027 TRITNL(80,0,1)= .0000013660
1028 TRITNL(80,1,1)= .0000013660
1029 TRITNL(80,2,1)= .0000039870
1030 TRITNL(80,3,1)= .0000063540
1031! LAW 82 (incompatible NPT=0)
1032! N = 1
1033 SHTNL(82,0,1) = .0000010020 !.0000006991
1034 SHTNL(82,1,1) = .0000010020 !.0000007022
1035 SHTNL(82,2,1) = .0000026670 !.0000017930
1036 SHTNL(82,3,1) = .0000043910 !.0000028300
1037!
1038 TRITNL(82,0,1)= .0000010910 !.0000006886
1039 TRITNL(82,1,1)= .0000010910 !.0000007115
1040 TRITNL(82,2,1)= .0000029140 !.0000017680
1041 TRITNL(82,3,1)= .0000048540 !.0000028740
1042! N = 2
1043 SHTNL(82,0,2) = .0000015170 !.0000010470
1044 SHTNL(82,1,2) = .0000015170 !.0000010430
1045 SHTNL(82,2,2) = .0000042590 !.0000027830
1046 SHTNL(82,3,2) = .0000070430 !.0000045710
1047!
1048 TRITNL(82,0,2)= .0000016810 !.0000010690
1049 TRITNL(82,1,2)= .0000016810 !.0000010330
1050 TRITNL(82,2,2)= .0000047360 !.0000028430
1051 TRITNL(82,3,2)= .0000077400 !.0000046350
1052! N = 3
1053 SHTNL(82,0,3) = .0000021070 !.0000013880
1054 SHTNL(82,1,3) = .0000021070 !.0000014040
1055 SHTNL(82,2,3) = .0000060060 !.0000038590
1056 SHTNL(82,3,3) = .0000099310 !.0000063230
1057!
1058 TRITNL(82,0,3)= .0000023390 !.0000014340
1059 TRITNL(82,1,3)= .0000023390 !.0000014220
1060 TRITNL(82,2,3)= .0000066970 !.0000038320
1061 TRITNL(82,3,3)= .0000110100 !.0000063450
1062! LAW 87
1063 SHTNL(87,0,1) = .0000004768 !.0000003592
1064 SHTNL(87,1,1) = .0000004768 !.0000003592
1065 SHTNL(87,2,1) = .0000010600 !.0000007907
1066 SHTNL(87,3,1) = .0000016710 !.0000011520
1067!
1068 TRITNL(87,0,1)= .0000004606 !.0000003754
1069 TRITNL(87,1,1)= .0000004606 !.0000003874
1070 TRITNL(87,2,1)= .0000010630 !.0000008403
1071 TRITNL(87,3,1)= .0000016680 !.0000012850
1072! LAW 93
1073 SHTNL(93,0,1) = .0000003473 !.0000002801
1074 SHTNL(93,1,1) = .0000003773 !.0000002801
1075 SHTNL(93,2,1) = .0000007325 !.0000005238
1076 SHTNL(93,3,1) = .0000011140 !.0000007912
1077!
1078 TRITNL(93,0,1)= .0000003500 !.0000002629
1079 TRITNL(93,1,1)= .0000003500 !.0000002629
1080 TRITNL(93,2,1)= .0000007356 !.0000005262
1081 TRITNL(93,3,1)= .0000011220 !.0000008110
1082! LAW 104
1083 ! Nice algo
1084 SHTNL(104,0,1) = .0000002978
1085 SHTNL(104,1,1) = .0000002978
1086 SHTNL(104,2,1) = .0000007051
1087 SHTNL(104,3,1) = .0000011480
1088 ! Newton algo
1089 SHTNL(104,0,2) = .0000004890
1090 SHTNL(104,1,2) = .0000004890
1091 SHTNL(104,2,2) = .0000009500
1092 SHTNL(104,3,2) = .0000014800
1093 ! Nice algo
1094 TRITNL(104,0,1) = .0000002978
1095 TRITNL(104,1,1) = .0000002978
1096 TRITNL(104,2,1) = .0000007051
1097 TRITNL(104,3,1) = .0000011480
1098 ! Newton algo
1099 TRITNL(104,0,2) = .0000004890
1100 TRITNL(104,1,2) = .0000004890
1101 TRITNL(104,2,2) = .0000009500
1102 TRITNL(104,3,2) = .0000014800
1103 ! iloc=1
1104 SHTNL_OPTION(1,1) = .0000000700
1105 ! iloc=2
1106 SHTNL_OPTION(1,2) = .0000000800
1107 ! iloc=3
1108 SHTNL_OPTION(1,3) = .0000003500
1109
1110C Temps de traitement de la partie elementaire sans la loi materielle fonction de IPARITH
1111 IF (IPARI0.EQ.1) THEN
1112! SHELT :
1113! 1-->5 --> /SHELL
1114! 6-->10 --> /SH_COMP or TYPE10
1115! 11-->15 --> /SH_SANDW or TYPE11
1116! 15-->20 --> /SH_ORTH or TYPE9
1117! 20-->25 --> /SH_FABR or TYPE16
1118! IDRILL = 0
1119! 1/6/11/16/21 --> Q4
1120! 2/7/12/17/22 --> QEPH
1121! 3/8/13/18/23 --> BATOZ
1122! IDRILL = 1
1123! none --> Q4
1124! 4/9 /14/19/24 --> QEPH
1125! 5/10/15/20/25 --> BATOZ
1126! TRIELT :
1127! 1 --> tri
1128! NLOCAL_OPTION
1129! 1 --> sol
1130! 2 --> shell
1131! /SHELL
1132 SHTELT(1) = .0000006048
1133 SHTELT(2) = .0000007659
1134 SHTELT(3) = .0000019960
1135 SHTELT(4) = .0000009898
1136 SHTELT(5) = .0000018360
1137! /SH_COMP
1138 SHTELT(6) = .0000007310
1139 SHTELT(7) = .0000010510
1140 SHTELT(8) = .0000021310
1141 SHTELT(9) = .0000011870
1142 SHTELT(10) = .0000018760
1143
1144! /SH_SANDW
1145 SHTELT(11) = .0000007010
1146 SHTELT(12) = .0000008640
1147 SHTELT(13) = .0000020750
1148 SHTELT(14) = .0000013510
1149 SHTELT(15) = .0000020500
1150
1151! /SH_ORTH
1152 SHTELT(16) = .0000006410
1153 SHTELT(17) = .0000008440
1154 SHTELT(18) = .0000020310
1155 SHTELT(19) = .0000016400
1156 SHTELT(20) = .0000025800
1157
1158! /SH_FABR
1159 SHTELT(21) = .0000009400
1160 SHTELT(22) = .0000010520
1161 SHTELT(23) = .0000023400
1162 SHTELT(24) = .0000016400
1163 SHTELT(25) = .0000025800
1164
1165! /TYPE51
1166 SHTELT(26) = .0000015500
1167 SHTELT(27) = .0000018100
1168 SHTELT(28) = .0000036800
1169 SHTELT(29) = .0000018600
1170 SHTELT(30) = .0000054300
1171
1172 TRITELT(1) = .0000004504
1173
1174 NLOCAL_OPTION(2) = .0000001100
1175 ELSE
1176! /SHELL
1177 SHTELT(1) = .0000006048
1178 SHTELT(2) = .0000007659
1179 SHTELT(3) = .0000019960
1180 SHTELT(4) = .0000009898
1181 SHTELT(5) = .0000018360
1182! /SH_COMP
1183 SHTELT(6) = .0000007310
1184 SHTELT(7) = .0000010510
1185 SHTELT(8) = .0000021310
1186 SHTELT(9) = .0000011870
1187 SHTELT(10) = .0000018760
1188
1189! /SH_SANDW
1190 SHTELT(11) = .0000007010
1191 SHTELT(12) = .0000008640
1192 SHTELT(13) = .0000020750
1193 SHTELT(14) = .0000013510
1194 SHTELT(15) = .0000020500
1195
1196! /SH_ORTH
1197 SHTELT(16) = .0000006410
1198 SHTELT(17) = .0000008440
1199 SHTELT(18) = .0000020310
1200 SHTELT(19) = .0000016400
1201 SHTELT(20) = .0000025800
1202
1203! /SH_FABR
1204 SHTELT(21) = .0000009400
1205 SHTELT(22) = .0000010520
1206 SHTELT(23) = .0000023400
1207 SHTELT(24) = .0000016400
1208 SHTELT(25) = .0000025800
1209
1210! /TYPE51
1211 SHTELT(26) = .0000015500
1212 SHTELT(27) = .0000018100
1213 SHTELT(28) = .0000036800
1214 SHTELT(29) = .0000018600
1215 SHTELT(30) = .0000054300
1216
1217 TRITELT(1) = .0000004504
1218
1219 NLOCAL_OPTION(2) = .0000001100
1220 ENDIF
1221
1222! il y a un facteur multiplicatif sur les coques de Batoz, car on rentre 4x dans cmain3.
1223
1224 BATOZMULT = 3.8
1225! /PRONY or /LPRONY overcost
1226 VISC_PRONY_COST = .0000000297
1227
1228! *******************************************************************
1229! Failure model
1230! *******************************************************************
1231! Johnson : 1
1232! Tbutcher : 2
1233! Wilkins : 3
1234! user1 : 4
1235! user2 : 5
1236! user3 : 6
1237! FLD : 7 (for shell only)
1238! Spalling : 8 (for solid only)
1239! Wierzbicki : 9
1240! Tensstrain : 10
1241! Energy : 11
1242! ??? : 12
1243! Chang : 13 (for shell only)
1244! Hashin : 14
1245! Yamada : 15
1246! Puck : 16
1247! XFEM_Tbutcher : 17
1248! Lad_dama : 18 (for solid only)
1249! Power_da : 19
1250! Connect : 20 (for solid only)
1251! XFEM_Johns : 21
1252! XFEM_FLD : 22
1253! TAB1 : 23
1254! Orthstra : 24 (for shell only)
1255! NXT : 25 (for shell only)
1256! Snconnect : 26 (for solid only)
1257! EMC : 27 (for solid only)
1258! Wind : 28 (for shell only)
1259! Sahraei1 : 29 (for solid only)
1260! Biquadratic : 30
1261! _
1262! SHELL : 1--> 0 ou 1pt |
1263! 2--> 3 pts |--- MLN < 28 or ==32
1264! 3--> 5 pts _|
1265! 4--> 0 ou 1pt |
1266! 5--> 3 pts |--- MLN >= 28 or /=32
1267! 6--> 5 pts _|
1268!
1269! SOLIDE : 1 ----> MLN <28 or ==49
1270! 2 ----> MLN >=28 or /=49
1271! --------------------------------------------------------
1272 DO I=1,MAXFAIL
1273 DO J=1,3
1274 RUPTURE_SOL(I,J) = 0.
1275 RUPTURE_SOL12(I,J) = 0.
1276 RUPTURE_TET4(I,J) = 0.
1277 RUPTURE_TET10(I,J) = 0.
1278 ENDDO
1279 DO J=1,6
1280 RUPTURE_SHELL(I,J) = 0.
1281 ENDDO
1282 ENDDO
1283! --------------------
1284! Johnson
1285 ! MLN <28 or ==32
1286 RUPTURE_SHELL(1,1) = .00000009911 !.00000007517
1287 RUPTURE_SHELL(1,2) = .00000026240 !.00000019090
1288 RUPTURE_SHELL(1,3) = .00000043620 !.00000031460
1289
1290 ! MLN >=28 or /=32
1291 RUPTURE_SHELL(1,4) = .00000003700 !.00000003021
1292 RUPTURE_SHELL(1,5) = .00000009762 !.00000007567
1293 RUPTURE_SHELL(1,6) = .00000015010 !.00000011620
1294
1295 ! MLN <28 or ==49
1296 RUPTURE_SOL(1,1) = .00000010440 !.00000008839
1297 RUPTURE_SOL12(1,1) = .00000032440 !.00000017710
1298 RUPTURE_TET4(1,1) = .00000010240 !.00000008021
1299 RUPTURE_TET10(1,1) = .00000040700 !.00000033420
1300 ! MLN >=28 or /=49
1301 RUPTURE_SOL(1,2) = .00000002771 !.00000002696
1302 RUPTURE_SOL12(1,2) = .00000021210 !.00000015790
1303 RUPTURE_TET4(1,2) = .00000003000 !.00000002663
1304 RUPTURE_TET10(1,2) = .00000011590 !.00000010360
1305! --------------------
1306! Tbutcher
1307 ! MLN <28 or ==49
1308 RUPTURE_SHELL(2,1) = .00000006051 !.00000004873
1309 RUPTURE_SHELL(2,2) = .00000016470 !.00000010840
1310 RUPTURE_SHELL(2,3) = .00000025060 !.00000016340
1311
1312 RUPTURE_SHELL(2,4) = .00000004710 !.00000003994
1313 RUPTURE_SHELL(2,5) = .00000013300 !.00000010130
1314 RUPTURE_SHELL(2,6) = .00000021180 !.00000015690
1315
1316 ! MLN <28 or ==49
1317 RUPTURE_SOL(2,1) = .00000013190 !.00000004718
1318 RUPTURE_SOL12(2,1) = .00000060160 !.00000031390
1319 RUPTURE_TET4(2,1) = .00000012970 !.00000010090
1320 RUPTURE_TET10(2,1) = .00000055000 !.00000040130
1321
1322 ! MLN >=28 or /=49
1323 RUPTURE_SOL(2,2) = .00000006271 !.00000004949
1324 RUPTURE_SOL12(2,2) = .00000047430 !.00000030290
1325 RUPTURE_TET4(2,2) = .00000005780 !.00000004449
1326 RUPTURE_TET10(2,2) = .00000024800 !.00000016800
1327! --------------------
1328! Wilkins
1329 ! MLN <28 or ==49
1330 RUPTURE_SHELL(3,1) = .00000016500 !.00000011660
1331 RUPTURE_SHELL(3,2) = .00000047410 !.00000032760
1332 RUPTURE_SHELL(3,3) = .00000079230 !.00000052680
1333
1334 RUPTURE_SHELL(3,4) = .00000015530 !.00000011410
1335 RUPTURE_SHELL(3,5) = .00000044620 !.00000032760
1336 RUPTURE_SHELL(3,6) = .00000075810 !.00000052910
1337
1338 ! MLN <28 or ==49
1339 RUPTURE_SOL(3,1) = .00000034600 !.00000023210
1340 RUPTURE_SOL12(3,1) = .00000224100 !.00000131700
1341 RUPTURE_TET4(3,1) = .00000024000 !.00000016440
1342 RUPTURE_TET10(3,1) = .00000114000 !.00000075060
1343
1344 ! MLN >=28 or /=49
1345 RUPTURE_SOL(3,2) = .00000026980 !.00000016870
1346 RUPTURE_SOL12(3,2) = .00000212500 !.00000129800
1347 RUPTURE_TET4(3,2) = .00000016890 !.00000011060
1348 RUPTURE_TET10(3,2) = .00000077140 !.00000051000
1349! --------------------
1350! FLD (only for shell)
1351 ! MLN <28 or ==32
1352 RUPTURE_SHELL(7,1) = .00000008881 !.00000010900
1353 RUPTURE_SHELL(7,2) = .00000024680 !.00000021590
1354 RUPTURE_SHELL(7,3) = .00000039630 !.00000031890
1355
1356 ! MLN >=28 or /=32
1357 RUPTURE_SHELL(7,4) = .00000008301 !.00000007162
1358 RUPTURE_SHELL(7,5) = .00000022760 !.00000019550
1359 RUPTURE_SHELL(7,6) = .00000036510 !.00000031220
1360! --------------------
1361! Spalling (only for solid)
1362
1363 ! MLN <28 or ==49
1364 RUPTURE_SOL(8,1) = .00000010860 !.00000009428
1365 RUPTURE_SOL12(8,1) = .00000039150 !.00000025070
1366 RUPTURE_TET4(8,1) = .00000010170 !.00000009402
1367 RUPTURE_TET10(8,1) = .00000044560 !.00000038110
1368
1369 ! MLN >=28 or /=49
1370 RUPTURE_SOL(8,2) = .00000003551 !.00000003871
1371 RUPTURE_SOL12(8,2) = .00000026700 !.00000023730
1372 RUPTURE_TET4(8,2) = .00000003400 !.00000003716
1373 RUPTURE_TET10(8,2) = .00000013870 !.00000014820
1374! --------------------
1375! Wierzbicki
1376 ! MLN <28 or ==49
1377 RUPTURE_SHELL(9,1) = .00000036760 !.00000025370
1378 RUPTURE_SHELL(9,2) = .00000099280 !.00000071910
1379 RUPTURE_SHELL(9,3) = .00000165500 !.00000118100
1380
1381 RUPTURE_SHELL(9,4) = .00000034950 !.00000024610
1382 RUPTURE_SHELL(9,5) = .00000100800 !.00000071940
1383 RUPTURE_SHELL(9,6) = .00000177900 !.00000118500
1384
1385 ! MLN <28 or ==49
1386 RUPTURE_SOL(9,1) = .00000010240 !.00000008443
1387 RUPTURE_SOL12(9,1) = .00000032750 !.00000011190
1388 RUPTURE_TET4(9,1) = .00000010000 !.00000008236
1389 RUPTURE_TET10(9,1) = .00000040690 !.00000032530
1390
1391 ! MLN >=28 or /=49
1392 RUPTURE_SOL(9,2) = .00000002500 !.00000002591
1393 RUPTURE_SOL12(9,2) = .00000019410 !.00000016310
1394 RUPTURE_TET4(9,2) = .00000002880 !.00000002483
1395 RUPTURE_TET10(9,2) = .00000012130 !.00000008821
1396! --------------------
1397! Tensstrain
1398 ! MLN <28 or ==32
1399 RUPTURE_SHELL(10,1) = .00000007661 !.00000009068
1400 RUPTURE_SHELL(10,2) = .00000020120 !.00000016430
1401 RUPTURE_SHELL(10,3) = .00000033860 !.00000023120
1402
1403 ! MLN >=28 or /=32
1404 RUPTURE_SHELL(10,4) = .00000005930!.00000004473
1405 RUPTURE_SHELL(10,5) = .00000014790 !.00000012650
1406 RUPTURE_SHELL(10,6) = .00000025670 !.00000019170
1407
1408 ! MLN <28 or ==49
1409 RUPTURE_SOL(10,1) = .00000020430 !.00000013990
1410 RUPTURE_SOL12(10,1) = .00000024800 !.00000025600
1411 RUPTURE_TET4(10,1) = .00000014700 !.00000010960
1412 RUPTURE_TET10(10,1) = .00000070500 !.00000047680
1413
1414 ! MLN >=28 or /=49
1415 RUPTURE_SOL(10,2) = .00000012620 !.00000008052
1416 RUPTURE_SOL12(10,2) = .00000098420 !.00000056320
1417 RUPTURE_TET4(10,2) = .00000007410 !.00000005430
1418 RUPTURE_TET10(10,2) = .00000040420 !.00000025250
1419! --------------------
1420! Energy
1421 ! MLN <28 or ==49
1422 RUPTURE_SHELL(11,1) = .00000006261 !.00000005295
1423 RUPTURE_SHELL(11,2) = .00000015600 !.00000011720
1424 RUPTURE_SHELL(11,3) = .00000026420 !.00000017760
1425
1426 RUPTURE_SHELL(11,4) = .00000005000 !.00000003770
1427 RUPTURE_SHELL(11,5) = .00000013750 !.00000010100
1428 RUPTURE_SHELL(11,6) = .00000021900 !.00000015840
1429
1430 ! MLN <28 or ==49
1431 RUPTURE_SOL(11,1) = .00000011580 !.00000009413
1432 RUPTURE_SOL12(11,1) = .00000025430 !.00000009015
1433 RUPTURE_TET4(11,1) = .00000011910 !.00000009164
1434 RUPTURE_TET10(11,1) = .00000048380 !.00000037960
1435
1436 ! MLN >=28 or /=49
1437 RUPTURE_SOL(11,2) = .00000004931 !.00000003883
1438 RUPTURE_SOL12(11,2) = .00000034180 !.00000025910
1439 RUPTURE_TET4(11,2) = .00000004170 !.00000003753
1440 RUPTURE_TET10(11,2) = .00000018220 !.00000014030
1441! --------------------
1442! Chang (only for shell)
1443 ! MLN <28 or ==49
1444 RUPTURE_SHELL(13,1) = .00000006261 !.00000004822
1445 RUPTURE_SHELL(13,2) = .00000017280 !.00000010740
1446 RUPTURE_SHELL(13,3) = .0000002850 !.00000016000
1447 ! MLN >=28 or /=49
1448 RUPTURE_SHELL(13,4) = .00000005630 !.00000003770
1449 RUPTURE_SHELL(13,5) = .00000015880 !.00000010100
1450 RUPTURE_SHELL(13,6) = .00000025590 !.00000015840
1451! --------------------
1452! Hashin
1453 ! MLN <28 or ==49
1454 RUPTURE_SHELL(14,1) = .00000009000 !.00000007989
1455 RUPTURE_SHELL(14,2) = .00000025000 !.00000025240
1456 RUPTURE_SHELL(14,3) = .00000040000 !.00000040300
1457
1458 RUPTURE_SHELL(14,4) = .00000009000 !.00000007989
1459 RUPTURE_SHELL(14,5) = .00000025000 !.00000025240
1460 RUPTURE_SHELL(14,6) = .00000045000 !.00000040300
1461
1462 ! MLN <28 or ==49
1463 RUPTURE_SOL(14,1) = .00000014000 !.00000013480
1464 RUPTURE_SOL12(14,1) = .00000080000 !.00000107400
1465 RUPTURE_TET4(14,1) = .00000011000 !.00000012500
1466 RUPTURE_TET10(14,1) = .00000043000 !.00000054770
1467
1468 ! MLN >=28 or /=49
1469 RUPTURE_SOL(14,2) = .00000010000 !.00000001289
1470 RUPTURE_SOL12(14,2) = .00000080000 !.00000010990
1471 RUPTURE_TET4(14,2) = .00000010000 !.00000001302
1472 RUPTURE_TET10(14,2) = .00000040000 !.00000005823
1473! --------------------
1474! Puck (compatible MLN<28)
1475 ! MLN <28 or ==49
1476 RUPTURE_SHELL(16,1) = .00000007741 !.00000007462
1477 RUPTURE_SHELL(16,2) = .00000020750 !.00000018470
1478 RUPTURE_SHELL(16,3) = .00000033010 !.00000028840
1479
1480 RUPTURE_SHELL(16,4) = .00000006340 !.00000006606
1481 RUPTURE_SHELL(16,5) = .00000017550 !.00000017190
1482 RUPTURE_SHELL(16,6) = .00000028220 !.00000027600
1483
1484 ! MLN <28 or ==49
1485 RUPTURE_SOL(16,1) = .00000015910 !.00000013060
1486 RUPTURE_SOL12(16,1) = .00000025340 !.00000008806
1487 RUPTURE_TET4(16,1) = .00000016450 !.00000012170
1488 RUPTURE_TET10(16,1) = .00000066190 !.00000042400
1489
1490 RUPTURE_SOL(16,2) = .00000002100 !.00000005419
1491 RUPTURE_SOL12(16,2) = .00000013110 !.00000007646
1492 RUPTURE_TET4(16,2) = .00000002000 !.00000001338
1493 RUPTURE_TET10(16,2) = .00000008000 !.00000004460
1494! --------------------
1495! Lad_Dama (compatible MLN<28) (only for solid)
1496 ! MLN <28 or ==49
1497 RUPTURE_SOL(18,1) = .00000016390 !.00000016890
1498 RUPTURE_SOL12(18,1) = .00000024830 !.00000008820
1499 RUPTURE_TET4(18,1) = .00000016810 !.00000013980
1500 RUPTURE_TET10(18,1) = .00000069150 !.00000058630
1501
1502 RUPTURE_SOL(18,2) = .00000002000 !.00000016890
1503 RUPTURE_SOL12(18,2) = .00000013110 !.00000008820
1504 RUPTURE_TET4(18,2) = .00000002000 !.00000013980
1505 RUPTURE_TET10(18,2) = .00000008000 !.00000058630
1506! --------------------
1507! Connect (compatible MLN==59 et 83)
1508 ! MLN==59
1509 RUPTURE_SOL(20,1) = .00000050000 !.00000054140
1510 RUPTURE_SOL12(20,1) = .00000050000 !.00000054140
1511 RUPTURE_TET4(20,1) = .00000050000 !.00000054140
1512 RUPTURE_TET10(20,1) = .00000050000 !.00000054140
1513 ! MLN==83
1514 RUPTURE_SOL(20,2) = .00000050000 !.00000054140
1515 RUPTURE_SOL12(20,2) = .00000050000 !.00000054140
1516 RUPTURE_TET4(20,2) = .00000050000 !.00000054140
1517 RUPTURE_TET10(20,2) = .00000050000 !.00000054140
1518! --------------------
1519! Tab1
1520 ! MLN <28 or ==49
1521 RUPTURE_SHELL(23,1) = .00000026330 !.00000021020
1522 RUPTURE_SHELL(23,2) = .00000075580 !.00000060850
1523 RUPTURE_SHELL(23,3) = .00000125000 !.00000098850
1524
1525 RUPTURE_SHELL(23,4) = .00000025340 !.00000021590
1526 RUPTURE_SHELL(23,5) = .00000074690 !.00000062580
1527 RUPTURE_SHELL(23,6) = .00000123700 !.0000009949
1528
1529 ! MLN <28 or ==49
1530 RUPTURE_SOL(23,1) = .00000039050 !.00000031120
1531 RUPTURE_SOL12(23,1) = .00000262300 !.00000192900
1532 RUPTURE_TET4(23,1) = .00000031240 !.00000026620
1533 RUPTURE_TET10(23,1) = .00000138800 !.00000118100
1534
1535 ! MLN >=28 or /=49
1536 RUPTURE_SOL(23,2) = .00000031310 !.00000025940
1537 RUPTURE_SOL12(23,2) = .00000248200 !.00000198600
1538 RUPTURE_TET4(23,2) = .00000024570 !.00000022330
1539 RUPTURE_TET10(23,2) = .00000108100 !.00000096400
1540! --------------------
1541! NXT (only for shell)
1542 ! MLN <28 or ==32
1543 RUPTURE_SHELL(25,1) = .00000011680 !.00000014880
1544 RUPTURE_SHELL(25,2) = .00000031700 !.00000032860
1545 RUPTURE_SHELL(25,3) = .00000051900 !.00000051100
1546
1547 ! MLN >=28 or /=32
1548 RUPTURE_SHELL(25,4) = .00000008961 !.00000008550
1549 RUPTURE_SHELL(25,5) = .00000024460 !.00000024030
1550 RUPTURE_SHELL(25,6) = .00000041580 !.00000038060
1551! --------------------
1552! Snconnect (compatible MLN==59 et 83)
1553 ! MLN==59
1554 RUPTURE_SOL(26,1) = .00000065000 !.00000065150
1555 RUPTURE_SOL12(26,1) = .00000065000 !.00000065150
1556 RUPTURE_TET4(26,1) = .00000065000 !.00000065150
1557 RUPTURE_TET10(26,1) = .00000065000 !.00000065150
1558 ! MLN==83
1559 RUPTURE_SOL(26,2) = .00000065000 !.00000065150
1560 RUPTURE_SOL12(26,2) = .00000065000 !.00000065150
1561 RUPTURE_TET4(26,2) = .00000065000 !.00000065150
1562 RUPTURE_TET10(26,2) = .00000065000 !.00000065150
1563! --------------------
1564! EMC (only for solid)
1565 ! MLN <28 or ==49
1566 RUPTURE_SOL(27,1) = .00000009892 !.00000007774
1567 RUPTURE_SOL12(27,1) = .00000029550 !.00000013030
1568 RUPTURE_TET4(27,1) = .00000009891 !.00000007539
1569 RUPTURE_TET10(27,1) = .00000040460 !.00000029830
1570
1571 ! MLN >=28 or /=49
1572 RUPTURE_SOL(27,2) = .00000003000 !.00000001913
1573 RUPTURE_SOL12(27,2) = .00000017000 !.00000010660
1574 RUPTURE_TET4(27,2) = .00000003000 !.00000001654
1575 RUPTURE_TET10(27,2) = .00000010840 !.00000006008
1576! --------------------
1577! Alter (only shell/tri)
1578 ! MLN <28 or ==49
1579 RUPTURE_SHELL(28,1) = .00000001664
1580 RUPTURE_SHELL(28,2) = .00000003502
1581 RUPTURE_SHELL(28,3) = .00000006578
1582
1583 RUPTURE_SHELL(28,4) = .00000001664
1584 RUPTURE_SHELL(28,5) = .00000003502
1585 RUPTURE_SHELL(28,6) = .00000006578
1586! --------------------
1587! Biquad
1588 ! MLN <28 or ==49
1589 RUPTURE_SHELL(30,1) = .00000004610 !.00000003069
1590 RUPTURE_SHELL(30,2) = .00000009822 !.00000005533
1591 RUPTURE_SHELL(30,3) = .00000016600 !.00000008252
1592
1593 RUPTURE_SHELL(30,4) = .00000003880 !.00000003133
1594 RUPTURE_SHELL(30,5) = .00000008732 !.00000006587
1595 RUPTURE_SHELL(30,6) = .00000014760 !.00000010300
1596
1597 ! MLN <28 or ==49
1598 RUPTURE_SOL(30,1) = .00000009382 !.00000008026
1599 RUPTURE_SOL12(30,1) = .00000028110 !.00000013200
1600 RUPTURE_TET4(30,1) = .00000009440 !.00000007523
1601 RUPTURE_TET10(30,1) = .00000038780 !.00000030150
1602
1603 RUPTURE_SOL(30,2) = .00000003000 !.00000001929
1604 RUPTURE_SOL12(30,2) = .00000017000 !.00000011930
1605 RUPTURE_TET4(30,2) = .00000002240 !.00000001793
1606 RUPTURE_TET10(30,2) = .00000010200 !.00000006588
1607! --------------------
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