OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
weights_p4linux964_spmd_avx512.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 Skylake machine (AVX512)
25! Used machine : Pico --> Intel(R) Xeon(R) Gold 6248 CPU @ 2.50GHz
26
27
28! temps de reference correspond au temps /elt/materiaux loi2 0 pt integ
29 TPSREF = .0000002492 !.0000001927 !.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) = .0000001047 !.0000000700 !.0000001512
64 SOL8TNL(I,J) = .0000004822 !.0000003250 !.0000006297
65 TET4TNL(I,J) = .0000000937 !.0000000640 !.0000001542
66 TET10TNL(I,J)= .0000003522 !.0000002339 !.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! LAW 01
74 SOL1TNL(1,1) = .0000000858
75 SOL8TNL(1,1) = .0000001996
76 TET4TNL(1,1) = .0000000762
77 TET10TNL(1,1)= .0000002815
78! LAW 01 + ISMSTR = 10 or 12
79 SOL1TNL(1,2) = .0000001487
80 SOL8TNL(1,2) = .0000012190
81 TET4TNL(1,2) = .0000001487
82 TET10TNL(1,2)= .0000006840
83! LAW 02
84 SOL1TNL(2,1) = .0000001047 !.0000000700 !.0000001512
85 SOL8TNL(2,1) = .0000004822 !.0000003250 !.0000006297
86 TET4TNL(2,1) = .0000000937 !.0000000640 !.0000001542
87 TET10TNL(2,1)= .0000003522 !.0000002339 !.0000005545
88! LAW 02 ICC=1
89 SOL1TNL(2,2) = .0000001302 !.0000000951 !.0000002051
90 SOL8TNL(2,2) = .0000006944 !.0000004699 !.0000010080
91 TET4TNL(2,2) = .0000001198 !.0000000890 !.0000001820
92 TET10TNL(2,2)= .0000004568 !.0000003051 !.0000007670
93! LAW 03
94 SOL1TNL(3,1) = .0000001450 !.0000001010 !.0000002145
95 SOL8TNL(3,1) = .0000011050 !.0000007280 !.0000006086
96 TET4TNL(3,1) = .0000001360 !.0000001006 !.0000002292
97 TET10TNL(3,1)= .0000005400 !.0000003600 !.0000009112
98! LAW 04
99 SOL1TNL(4,1) = .0000001810 !.0000001150 !.0000001555
100 SOL8TNL(4,1) = .0000014110 !.0000009200 !.0000006050
101 TET4TNL(4,1) = .0000001699 !.0000001134 !.0000002547
102 TET10TNL(4,1)= .0000006460 !.0000004300 !.0000009900
103! LAW 06
104 SOL1TNL(6,1) = .0000001300 !.0000000864 !.0000002483
105 SOL8TNL(6,1) = .0000009100 !.0000006910 !.0000006004
106 TET4TNL(6,1) = .0000001200 !.0000000820 !.0000002412
107 TET10TNL(6,1)= .0000004700 !.0000003280 !.0000007858
108! LAW 10
109 SOL1TNL(10,1) = .0000001128 !.0000000732 !.0000002024
110 SOL8TNL(10,1) = .0000008037 !.0000005992 !.0000013650
111 TET4TNL(10,1) = .0000001001 !.0000000660 !.0000001815
112 TET10TNL(10,1)= .0000003798 !.0000002570 !.0000006534
113! LAW 12
114 SOL1TNL(12,1) = .0000002394 !.0000001630 !.0000003941
115 SOL8TNL(12,1) = .0000018070 !.0000013930 !.0000038700
116 TET4TNL(12,1) = .0000002244 !.0000001545 !.0000003670
117 TET10TNL(12,1)= .0000008822 !.0000005915 !.0000013560
118! LAW 21
119 SOL1TNL(21,1) = .0000001126 !.0000000700 !.0000001872
120 SOL8TNL(21,1) = .0000008074 !.0000005600 !.0000013840
121 TET4TNL(21,1) = .0000000996 !.0000000906 !.0000001811
122 TET10TNL(21,1)= .0000003798 !.0000002529 !.0000006885
123! LAW 22
124 SOL1TNL(22,1) = .0000001180 !.0000000740 !.0000001854
125 SOL8TNL(22,1) = .0000008551 !.0000006025 !.0000013990
126 TET4TNL(22,1) = .0000001050 !.0000000696 !.0000001663
127 TET10TNL(22,1)= .0000004120 !.0000002694 !.0000007349
128! LAW 23
129 SOL1TNL(23,1) = .0000001200 !.0000000780 !.0000001890
130 SOL8TNL(23,1) = .0000008699 !.0000006269 !.0000014360
131 TET4TNL(23,1) = .0000001073 !.0000000708 !.0000001850
132 TET10TNL(23,1)= .0000004130 !.0000002767 !.0000007381
133! LAW 24
134 SOL1TNL(24,1) = .0000003024 !.0000002029 !.0000006082
135 SOL8TNL(24,1) = .0000022760 !.0000015870 !.0000042810
136 TET4TNL(24,1) = .0000002885 !.0000001965 !.0000004260
137 TET10TNL(24,1)= .0000011390 !.0000008094 !.0000020050
138! LAW 25
139 SOL1TNL(25,1) = .0000002549
140 SOL8TNL(25,1) = .0000019950
141 TET4TNL(25,1) = .0000002605
142 TET10TNL(25,1)= .0000010080
143! LAW 25 Crasurv
144 SOL1TNL(25,2) = .0000003029
145 SOL8TNL(25,2) = .0000024380
146 TET4TNL(25,2) = .0000003029
147 TET10TNL(25,2)= .0000011930
148! LAW 28
149 SOL1TNL(28,1) = .0000002598 !.0000002796 !.0000003664
150 SOL8TNL(28,1) = .0000017400 !.0000019590 !.0000021385
151 TET4TNL(28,1) = .0000002400 !.0000002706 !.0000003168
152 TET10TNL(28,1)= .0000010550 !.0000012226 !.0000013620
153! LAW 33 0 fct
154 SOL1TNL(33,1) = .0000002037 !.0000001307 !.0000004453
155 SOL8TNL(33,1) = .0000012670 !.0000008472 !.0000021870
156 TET4TNL(33,1) = .0000001778 !.0000001152 !.0000004027
157 TET10TNL(33,1)= .0000006775 !.0000004627 !.0000015750
158! LAW 33 1 fct
159 SOL1TNL(33,2) = .0000003314 !.0000002168 !.0000006103
160 SOL8TNL(33,2) = .0000023110 !.0000015580 !.0000025320
161 TET4TNL(33,2) = .0000003050 !.0000001975 !.0000005659
162 TET10TNL(33,2)= .0000012100 !.0000007945 !.0000018270
163! LAW 34
164 SOL1TNL(34,1) = .0000001404 !.0000000885 !.0000002516
165 SOL8TNL(34,1) = .0000007640 !.0000005305 !.0000012240
166 TET4TNL(34,1) = .0000001307 !.0000000864 !.0000002258
167 TET10TNL(34,1)= .0000005092 !.0000003467 !.0000008617
168! LAW 35
169 SOL1TNL(35,1) = .0000001512 !.0000000964 !.0000003070
170 SOL8TNL(35,1) = .0000008688 !.0000005900 !.0000017020
171 TET4TNL(35,1) = .0000001460 !.0000000961 !.0000002691
172 TET10TNL(35,1)= .0000005651 !.0000003747 !.0000010870
173! LAW 36 1 fonction
174 SOL1TNL(36,1) = .0000001721 !.0000001103 !.0000003096
175 SOL8TNL(36,1) = .0000009919 !.0000006723 !.0000016040
176 TET4TNL(36,1) = .0000001632 !.0000001070 !.0000002877
177 TET10TNL(36,1)= .0000006291 !.0000004051 !.0000009584
178! LAW 36 5 fonctions
179 SOL1TNL(36,2) = .0000001877 !.0000001200 !.0000003358
180 SOL8TNL(36,2) = .0000014650 !.0000006900 !.0000016420
181 TET4TNL(36,2) = .0000001750 !.0000001100 !.0000002724
182 TET10TNL(36,2)= .0000006650 !.0000004100 !.0000010120
183! LAW 36 10 fonctions
184 SOL1TNL(36,3) = .0000002064 !.0000001300 !.0000003702
185 SOL8TNL(36,3) = .0000016090 !.0000007000 !.0000016950
186 TET4TNL(36,3) = .0000002000 !.0000001150 !.0000002854
187 TET10TNL(36,3)= .0000007800 !.0000004200 !.0000010970
188! LAW 38
189 SOL1TNL(38,1) = .0000005560
190 SOL8TNL(38,1) = .0000041000
191 TET4TNL(38,1) = .0000005560
192 TET10TNL(38,1)= .0000022240
193! LAW 38 + ISMSTR=10 or 12
194 SOL1TNL(38,2) = .0000006012
195 SOL8TNL(38,2) = .0000047800
196 TET4TNL(38,2) = .0000006012
197 TET10TNL(38,2)= .0000023900
198! LAW 40
199 SOL1TNL(40,1) = .0000002040 !.0000001303 !.0000004917
200 SOL8TNL(40,1) = .0000013660 !.0000008978 !.0000030250
201 TET4TNL(40,1) = .0000001898 !.0000001282 !.0000004015
202 TET10TNL(40,1)= .0000007526 !.0000004790 !.0000013700
203! LAW 41
204 SOL1TNL(41,1) = .0000002549 !.0000001671 !.0000004092
205 SOL8TNL(41,1) = .0000019880 !.0000020390 !.0000030250
206 TET4TNL(41,1) = .0000002874 !.0000001923 !.0000004068
207 TET10TNL(41,1)= .0000011320 !.0000007760 !.0000017520
208! LAW 42
209 SOL1TNL(42,1) = .0000002754 !.0000001783 !.0000006159
210 SOL8TNL(42,1) = .0000018460 !.0000012270 !.0000041410
211 TET4TNL(42,1) = .0000002556 !.0000001650 !.0000006126
212 TET10TNL(42,1)= .0000009823 !.0000006596 !.0000023340
213C M = 1
214 SOL1TNL(42,2) = .0000003066 !.0000002003 !.0000006763
215 SOL8TNL(42,2) = .0000020470 !.0000013980 !.0000045820
216 TET4TNL(42,2) = .0000002822 !.0000001855 !.0000006537
217 TET10TNL(42,2)= .0000010900 !.0000007557 !.0000023430
218C M = 2
219 SOL1TNL(42,3) = .0000003275 !.0000002193 !.0000007634
220 SOL8TNL(42,3) = .0000021700 !.0000014950 !.0000051210
221 TET4TNL(42,3) = .0000002901 !.0000001899 !.0000006823
222 TET10TNL(42,3)= .0000011330 !.0000007831 !.0000027020
223! LAW 44
224 SOL1TNL(44,1) = .0000002750
225 SOL8TNL(44,1) = .0000022000
226 TET4TNL(44,1) = .0000002775
227 TET10TNL(44,1)= .0000011000
228! LAW 48
229 SOL1TNL(48,1) = .0000001872 !.0000001223 !.0000002994
230 SOL8TNL(48,1) = .0000012240 !.0000008252 !.0000015830
231 TET4TNL(48,1) = .0000001781 !.0000001190 !.0000002660
232 TET10TNL(48,1)= .0000007016 !.0000004441 !.0000009488
233! LAW 50
234 SOL1TNL(50,1) = .0000003108 !.0000002051 !.0000004482
235 SOL8TNL(50,1) = .0000021920 !.0000014810 !.0000027100
236 TET4TNL(50,1) = .0000002728 !.0000001858 !.0000004178
237 TET10TNL(50,1)= .0000012690 !.0000008257 !.0000018130
238! LAW 52
239 SOL1TNL(52,1) = .0000001674 !.0000001121 !.0000002826
240 SOL8TNL(52,1) = .0000010590 !.0000007335 !.0000015560
241 TET4TNL(52,1) = .0000001549 !.0000001024 !.0000002650
242 TET10TNL(52,1)= .0000006149 !.0000003935 !.0000010590
243! LAW 53
244 SOL1TNL(53,1) = .0000002179 !.0000001408 !.0000004006
245 SOL8TNL(53,1) = .0000013770 !.0000009253 !.0000022760
246 TET4TNL(53,1) = .0000002065 !.0000001358 !.0000003889
247 TET10TNL(53,1)= .0000008063 !.0000005273 !.0000014540
248! LAW 59 (user solid property) IGT=43
249 SOL1TNL(59,1) = .0000001496 !.0000001028 !.0000002207
250 SOL8TNL(59,1) = .0000001496 !.0000001028 !.0000002207
251 TET4TNL(59,1) = .0000001496 !.0000001028 !.0000002207
252 TET10TNL(59,1)= .0000001496 !.0000001028 !.0000002207
253! LAW 60
254 SOL1TNL(60,1) = .0000002415 !.0000001584 !.0000004239
255 SOL8TNL(60,1) = .0000015600 !.0000010710 !.0000024010
256 TET4TNL(60,1) = .0000002174 !.0000001485 !.0000003718
257 TET10TNL(60,1)= .0000008970 !.0000006158 !.0000014290
258! LAW 62
259 ! M = 0
260 SOL1TNL(62,1) = .0000004174 !.0000002929 !.0000009680
261 SOL8TNL(62,1) = .0000032200 !.0000023400 !.0000068900
262 TET4TNL(62,1) = .0000003950 !.0000002800 !.0000009896
263 TET10TNL(62,1)= .0000015010 !.0000011200 !.0000039430
264 ! M = 1
265 SOL1TNL(62,2) = .0000004554 !.0000003280 !.0000010850
266 SOL8TNL(62,2) = .0000035520 !.0000026700 !.0000078600
267 TET4TNL(62,2) = .0000004300 !.0000003200 !.0000010760
268 TET10TNL(62,2)= .0000016300 !.0000013200 !.0000043120
269 ! M = 2
270 SOL1TNL(62,3) = .0000004789 !.0000003460 !.0000011420
271 SOL8TNL(62,3) = .0000037800 !.0000028700 !.0000084600
272 TET4TNL(62,3) = .0000004700 !.0000003400 !.0000011150
273 TET10TNL(62,3)= .0000017860 !.0000013600 !.0000045435
274! LAW 65
275 SOL1TNL(65,1) = .0000002518 !.0000001685 !.0000004282
276 SOL8TNL(65,1) = .0000018560 !.0000014530 !.0000030270
277 TET4TNL(65,1) = .0000002325 !.0000001563 !.0000003897
278 TET10TNL(65,1)= .0000009361 !.0000006260 !.0000015160
279! LAW 66
280 SOL1TNL(66,1) = .0000001995 !.0000001312 !.0000005367
281 SOL8TNL(66,1) = .0000014480 !.0000009792 !.0000031200
282 TET4TNL(66,1) = .0000001903 !.0000001303 !.0000003860
283 TET10TNL(66,1)= .0000007459 !.0000004765 !.0000015470
284! LAW 68
285 SOL1TNL(68,1) = .0000004163 !.0000002812 !.0000006413
286 SOL8TNL(68,1) = .0000032500 !.0000022450 !.0000028940
287 TET4TNL(68,1) = .0000003950 !.0000002800 !.0000006475
288 TET10TNL(68,1)= .0000015100 !.0000011200 !.0000025900
289! LAW 69
290 ! M = 0 (no visc)
291 SOL1TNL(69,1) = .0000002307 !.0000001474 !.0000005245
292 SOL8TNL(69,1) = .0000016880 !.0000012990 !.0000040350
293 TET4TNL(69,1) = .0000002146 !.0000001417 !.0000005135
294 TET10TNL(69,1)= .0000008195 !.0000005402 !.0000020020
295 ! M = 1 (visc)
296 SOL1TNL(69,2) = .0000003525 !.0000002314 !.0000006120
297 SOL8TNL(69,2) = .0000026350 !.0000016670 !.0000048180
298 TET4TNL(69,2) = .0000003335 !.0000002211 !.0000006347
299 TET10TNL(69,2)= .0000012980 !.0000008437 !.0000024700
300 ! M = 2 (visc)
301 SOL1TNL(69,3) = .0000003577 !.0000002422 !.0000006370
302 SOL8TNL(69,3) = .0000026840 !.0000016960 !.0000051600
303 TET4TNL(69,3) = .0000003379 !.0000002260 !.0000006740
304 TET10TNL(69,3)= .0000013290 !.0000008731 !.0000025500
305! LAW 70
306! 1 funct / 22 points
307 SOL1TNL(70,1) = .0000001832 !.0000002354 !.0000003853
308 SOL8TNL(70,1) = .0000015060 !.0000015510 !.0000026830
309 TET4TNL(70,1) = .0000001800 !.0000002312 !.0000003395
310 TET10TNL(70,1)= .0000007280 !.0000008120 !.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) = .0000003106
323 SOL8TNL(71,1) = .0000022870
324 TET4TNL(71,1) = .0000001963
325 TET10TNL(71,1)= .0000007662
326! LAW 72
327 SOL1TNL(72,1) = .0000001604 !.0000001024 !.0000004335
328 SOL8TNL(72,1) = .0000011490 !.0000007460 !.0000032810
329 TET4TNL(72,1) = .0000001526 !.0000000999 !.0000003947
330 TET10TNL(72,1)= .0000005878 !.0000003764 !.0000016170
331! LAW 74
332 SOL1TNL(74,1) = .0000002353 !.0000001546 !.0000003704
333 SOL8TNL(74,1) = .0000017440 !.0000013200 !.0000027610
334 TET4TNL(74,1) = .0000002244 !.0000001514 !.0000003402
335 TET10TNL(74,1)= .0000008874 !.0000005876 !.0000013590
336! LAW 76
337 SOL1TNL(76,1) = .0000001691 !.0000001097 !.0000002745
338 SOL8TNL(76,1) = .0000012200 !.0000007951 !.0000021280
339 TET4TNL(76,1) = .0000001618 !.0000001073 !.0000002739
340 TET10TNL(76,1)= .0000006237 !.0000004039 !.0000011500
341! LAW 77
342 SOL1TNL(77,1) = .0000003195 !.0000002322 !.0000004887
343 SOL8TNL(77,1) = .0000024900 !.0000018570 !.0000038100
344 TET4TNL(77,1) = .0000003000 !.0000002280 !.0000005200
345 TET10TNL(77,1)= .0000012000 !.0000009120 !.0000020120
346! LAW 78
347 SOL1TNL(78,1) = .0000001681 !.0000001091 !.0000003087
348 SOL8TNL(78,1) = .0000011910 !.0000007906 !.0000023090
349 TET4TNL(78,1) = .0000001589 !.0000001072 !.0000002902
350 TET10TNL(78,1)= .0000006161 !.0000003944 !.0000012520
351! LAW 79
352 SOL1TNL(79,1) = .0000001536 !.0000001013 !.0000003060
353 SOL8TNL(79,1) = .0000010960 !.0000007296 !.0000020220
354 TET4TNL(79,1) = .0000001460 !.0000000969 !.0000002684
355 TET10TNL(79,1)= .0000005627 !.0000003681 !.0000011400
356! LAW 81
357 SOL1TNL(81,1) = .0000003597 !.0000002316 !.0000004223
358 SOL8TNL(81,1) = .0000027200 !.0000016720 !.0000031430
359 TET4TNL(81,1) = .0000003380 !.0000002159 !.0000002750
360 TET10TNL(81,1)= .0000012840 !.0000008423 !.0000012470
361! LAW 82
362 ! M = 1
363 SOL1TNL(82,1) = .0000004500
364 SOL8TNL(82,1) = .0000036000
365 TET4TNL(82,1) = .0000004400
366 TET10TNL(82,1)= .0000017600
367 ! M = 2
368 SOL1TNL(82,2) = .0000005500
369 SOL8TNL(82,2) = .0000045000
370 TET4TNL(82,2) = .0000005400
371 TET10TNL(82,2)= .0000020600
372 ! M = 3
373 SOL1TNL(82,3) = .0000007300
374 SOL8TNL(82,3) = .0000058400
375 TET4TNL(82,3) = .0000007180
376 TET10TNL(82,3)= .0000028700
377! LAW 83 (user solid property) IGT=43
378 SOL1TNL(83,1) = .0000011010 !.0000087740 !.0000013690
379 SOL8TNL(83,1) = .0000011010 !.0000087740 !.0000013690
380 TET4TNL(83,1) = .0000011010 !.0000087740 !.0000013690
381 TET10TNL(83,1)= .0000011010 !.0000087740 !.0000013690
382! LAW 84
383 SOL1TNL(84,1) = .0000001577 !.0000001032 !.0000002771
384 SOL8TNL(84,1) = .0000011260 !.0000007414 !.0000018140
385 TET4TNL(84,1) = .0000001490 !.0000001001 !.0000002489
386 TET10TNL(84,1)= .0000005753 !.0000003790 !.0000007861
387! LAW 88
388 SOL1TNL(88,1) = .0000039930
389 SOL8TNL(88,1) = .0000301000
390 TET4TNL(88,1) = .0000039930
391 TET10TNL(88,1)= .0000156000
392! LAW 90
393 SOL1TNL(90,1) = .0000005812
394 SOL8TNL(90,1) = .0000040500
395 TET4TNL(90,1) = .0000005812
396 TET10TNL(90,1)= .0000022080
397! LAW 90 + ISMSTR=10 or 12
398 SOL1TNL(90,2) = .0000007010
399 SOL8TNL(90,2) = .0000051090
400 TET4TNL(90,2) = .0000007010
401 TET10TNL(90,2)= .0000026700
402! LAW 92
403 SOL1TNL(92,1) = .0000005000
404 SOL8TNL(92,1) = .0000038000
405 TET4TNL(92,1) = .0000005000
406 TET10TNL(92,1)= .0000019000
407! LAW 92 + ISMSTR=10 or 12
408 SOL1TNL(92,2) = .0000005300
409 SOL8TNL(92,2) = .0000041000
410 TET4TNL(92,2) = .0000005300
411 TET10TNL(92,2)= .0000021000
412! LAW 94
413 SOL1TNL(94,1) = .0000002842
414 SOL8TNL(94,1) = .0000022700
415 TET4TNL(94,1) = .0000002842
416 TET10TNL(94,1)= .0000011040
417! LAW 94 + ISMSTR=10 or 12
418 SOL1TNL(94,2) = .0000003093
419 SOL8TNL(94,2) = .0000025000
420 TET4TNL(94,2) = .0000003093
421 TET10TNL(94,2)= .0000012300
422! LAW 95
423 SOL1TNL(95,1) = .0000002493
424 SOL8TNL(95,1) = .0000018090
425 TET4TNL(95,1) = .0000002246
426 TET10TNL(95,1)= .0000009194
427! LAW 97
428 SOL1TNL(97,1) = .0000001729
429 SOL8TNL(97,1) = .0000012030
430 TET4TNL(97,1) = .0000001560
431 TET10TNL(97,1)= .0000006190
432
433! LAW 100 : material cost = base cost + viscosity cost + N * network cost
434! LAW 100 (base)
435 SOL1TNL(100,1) = .0000002280
436 SOL8TNL(100,1) = .0000002280
437 TET4TNL(100,1) = .0000002280
438 TET10TNL(100,1)= .0000002280
439
440! LAW 100 (viscosity)
441 SOL1TNL(100,2) = .0000001390
442 SOL8TNL(100,2) = .0000001390
443 TET4TNL(100,2) = .0000001390
444 TET10TNL(100,2)= .0000001390
445
446! LAW 100 (network)
447 SOL1TNL(100,3) = .0000002000
448 SOL8TNL(100,3) = .0000002000
449 TET4TNL(100,3) = .0000002000
450 TET10TNL(100,3)= .0000002000
451! LAW 101
452 SOL1TNL(101,1) = .0000045120
453 SOL8TNL(101,1) = .0000356000
454 TET4TNL(101,1) = .0000045120
455 TET10TNL(101,1)= .0000179000
456! LAW 102
457 SOL1TNL(102,1) = .0000001649
458 SOL8TNL(102,1) = .0000011640
459 TET4TNL(102,1) = .0000001553
460 TET10TNL(102,1)= .0000006200
461! LAW 103
462 SOL1TNL(103,1) = .0000001758
463 SOL8TNL(103,1) = .0000012850
464 TET4TNL(103,1) = .0000001667
465 TET10TNL(103,1)= .0000006700
466! LAW 104
467 ! Nice algo
468 SOL1TNL(104,1) = .0000001464
469 SOL8TNL(104,1) = .0000001464
470 TET4TNL(104,1) = .0000001464
471 TET10TNL(104,1) = .0000001464
472 ! Newton algo
473 SOL1TNL(104,2) = .0000001843
474 SOL8TNL(104,2) = .0000001843
475 TET4TNL(104,2) = .0000001843
476 TET10TNL(104,2) = .0000001843
477 ! iloc=1
478 SOL_OPTION(1,1) = .0000000400
479 ! iloc=2
480 SOL_OPTION(1,2) = .0000000500
481 ! iloc=3
482 SOL_OPTION(1,3) = .0000000900
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) = .0000007818
490 SOL8TNL(187,1) = .0000058570
491 TET4TNL(187,1) = .0000002434
492 TET10TNL(187,1)= .0000012160
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) =.0000003163
523 SOLTELT(2) =.0000002903
524 SOLTELT(3) =.0000004425
525 SOLTELT(4) =.0000006316
526 SOLTELT(5) =.0000014720
527 SOLTELT(6) =.0000001380
528 SOLTELT(7) =.0000011885
529 SOLTELT(8) =.0000005360
530 SOLTELT(9) =.0000014650
531 SOLTELT(10) =.0000001923
532 SOLTELT(11) =.0000005789
533 SOLTELT(12) =.0000000466
534
535 NLOCAL_OPTION(1) = .0000001100
536 TET4TELT(1) =.0000001979
537 TET4TELT(2) =.0000007503
538 TET10TELT(1)=.0000008108
539
540 ELM_OVER_COST(1) = .0000000600
541 ELSE
542 SOLTELT(1) =.0000002635
543 SOLTELT(2) =.0000002419
544 SOLTELT(3) =.0000003687
545 SOLTELT(4) =.0000005263
546 SOLTELT(5) =.0000012780
547 SOLTELT(6) =.0000001219
548 SOLTELT(7) =.0000009904
549 SOLTELT(8) =.0000004466
550 SOLTELT(9) =.0000012200
551 SOLTELT(10) =.0000001600
552 SOLTELT(11) =.0000004847
553 SOLTELT(12) =.0000000388
554 NLOCAL_OPTION(1) = .0000001100
555 TET4TELT(1) =.0000001646
556 TET4TELT(2) =.0000006252
557 TET10TELT(1)=.0000006756
558
559 ELM_OVER_COST(1) = .0000000500
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) = .0000000585
615 SHTNL(J,1,I) = .0000000903
616 SHTNL(J,2,I) = .0000001578
617 SHTNL(J,3,I) = .0000002326
618C
619 TRITNL(J,0,I)= .0000000620
620 TRITNL(J,1,I)= .0000000686
621 TRITNL(J,2,I)= .0000001132
622 TRITNL(J,3,I)= .0000001630
623 ENDDO
624 ENDDO
625! special option for shell or shell3n
626 SHTNL_OPTION(1:10,1:5) = 0.
627! --------------------
628! LAW 01
629 SHTNL(1,0,1) = .0000000440
630 SHTNL(1,1,1) = .0000000440
631 SHTNL(1,2,1) = .0000000440
632 SHTNL(1,3,1) = .0000000453
633C
634 TRITNL(1,0,1)= .0000000360
635 TRITNL(1,1,1)= .0000000343
636 TRITNL(1,2,1)= .0000000361
637 TRITNL(1,3,1)= .0000000368
638! LAW 02
639 SHTNL(2,0,1) = .0000000585 !.0000000633
640 SHTNL(2,1,1) = .0000000903 !.0000000804
641 SHTNL(2,2,1) = .0000001578 !.0000001340
642 SHTNL(2,3,1) = .0000002326 !.0000001894
643C
644 TRITNL(2,0,1)= .0000000620
645 TRITNL(2,1,1)= .0000000686 !.0000000750
646 TRITNL(2,2,1)= .0000001132 !.0000001250
647 TRITNL(2,3,1)= .0000001630 !.0000001800
648! LAW 02 CC
649 SHTNL(2,0,2) = .0000000684 !.0000000730
650 SHTNL(2,1,2) = .0000000974 !.0000000877
651 SHTNL(2,2,2) = .0000001785 !.0000001517
652 SHTNL(2,3,2) = .0000002663 !.0000002208
653C
654 TRITNL(2,0,2)= .0000000720
655 TRITNL(2,1,2)= .0000000750 !.0000000850
656 TRITNL(2,2,2)= .0000001355 !.0000001500
657 TRITNL(2,3,2)= .0000001979 !.0000002150
658! LAW 15
659 SHTNL(15,0,1) = .0000001666 !.0000001550
660 SHTNL(15,1,1) = .0000001670 !.0000001550 !.0000001547
661 SHTNL(15,2,1) = .0000003758 !.0000003610
662 SHTNL(15,3,1) = .0000005853 !.0000005627
663C
664 TRITNL(15,0,2)= .0000001480
665 TRITNL(15,1,2)= .0000001475
666 TRITNL(15,2,2)= .0000003550
667 TRITNL(15,3,2)= .0000005580
668
669! LAW 19
670 SHTNL(19,0,1) = .0000000707 !.0000000734
671 SHTNL(19,1,1) = .0000000707 !.0000000738
672 SHTNL(19,2,1) = .0000001154 !.0000001203
673 SHTNL(19,3,1) = .0000001609 !.0000001628
674C
675 TRITNL(19,0,1)= .0000000650
676 TRITNL(19,1,1)= .0000000650
677 TRITNL(19,2,1)= .0000001150
678 TRITNL(19,3,1)= .0000001550
679! LAW 22
680 SHTNL(22,0,1) = .0000000908 !.0000000684
681 SHTNL(22,1,1) = .0000001131 !.0000001121
682 SHTNL(22,2,1) = .0000002179 !.0000001954
683 SHTNL(22,3,1) = .0000003294 !.0000003015
684C
685 TRITNL(22,0,1)= .0000000908
686 TRITNL(22,1,1)= .0000000963 !.0000000984
687 TRITNL(22,2,1)= .0000001935 !.0000001865
688 TRITNL(22,3,1)= .0000002945 !.0000002989
689! LAW 25
690 SHTNL(25,0,1) = .0000001202
691 SHTNL(25,1,1) = .0000001208
692 SHTNL(25,2,1) = .0000002481
693 SHTNL(25,3,1) = .0000003783
694C
695 TRITNL(25,0,1)= .0000001100
696 TRITNL(25,1,1)= .0000001120
697 TRITNL(25,2,1)= .0000002400
698 TRITNL(25,3,1)= .0000003800
699! LAW 25 Crasurv (/TYPE11)
700 SHTNL(25,0,2) = .0000002129
701 SHTNL(25,1,2) = .0000002120
702 SHTNL(25,2,2) = .0000004848
703 SHTNL(25,3,2) = .0000007574
704C
705 TRITNL(25,0,2)= .0000001750
706 TRITNL(25,1,2)= .0000001750
707 TRITNL(25,2,2)= .0000004600
708 TRITNL(25,3,2)= .0000006600
709! LAW 27 (0pt incompatible)
710 SHTNL(27,0,1) = .0000001270 !.0000001400
711 SHTNL(27,1,1) = .0000001270 !.0000001416
712 SHTNL(27,2,1) = .0000002660 !.0000003213
713 SHTNL(27,3,1) = .0000004185 !.0000005054
714C
715 TRITNL(27,0,1)= .0000001053
716 TRITNL(27,1,1)= .0000001119 !.0000001053
717 TRITNL(27,2,1)= .0000002537 !.0000002544
718 TRITNL(27,3,1)= .0000004037 !.0000004107
719! LAW 32
720 SHTNL(32,0,1) = .0000000792 !.0000000665
721 SHTNL(32,1,1) = .0000000792 !.0000000665
722 SHTNL(32,2,1) = .0000001343 !.0000001161
723 SHTNL(32,3,1) = .0000001922 !.0000001711
724C
725 TRITNL(32,0,1)= .0000000713 !.0000000601
726 TRITNL(32,1,1)= .0000000713 !.0000000601
727 TRITNL(32,2,1)= .0000001360 !.0000001219
728 TRITNL(32,3,1)= .0000002002 !.0000001836
729! LAW 35 (0pt incompatible)
730 SHTNL(35,0,1) = .0000001391
731 SHTNL(35,1,1) = .0000001391
732 SHTNL(35,2,1) = .0000002713
733 SHTNL(35,3,1) = .0000004261
734C
735 TRITNL(35,0,1)= .0000000980
736 TRITNL(35,1,1)= .0000000980
737 TRITNL(35,2,1)= .0000001960
738 TRITNL(35,3,1)= .0000003919
739! LAW 36 - 1 fct
740 SHTNL(36,0,1) = .0000001000 !.0000000964 !.0000002453
741 SHTNL(36,1,1) = .0000001281 !.0000001215 !.0000002975
742 SHTNL(36,2,1) = .0000002456 !.0000002050 !.0000005370
743 SHTNL(36,3,1) = .0000003851 !.0000002888 !.0000007957
744! 5 fct
745 SHTNL(36,0,2) = .0000001050 !.0000001076 !.0000002832
746 SHTNL(36,1,2) = .0000001438 !.0000001417 !.0000003655
747 SHTNL(36,2,2) = .0000002724 !.0000002480 !.0000005902
748 SHTNL(36,3,2) = .0000004310 !.0000003154 !.0000008500
749! 10 fct
750 SHTNL(36,0,3) = .0000001150 !.0000001301 !.0000003180
751 SHTNL(36,1,3) = .0000001620 !.0000001550 !.0000003850
752 SHTNL(36,2,3) = .0000003183 !.0000002624 !.0000006500
753 SHTNL(36,3,3) = .0000005000 !.0000003550 !.0000009020
754
755! LAW 36 - 1 fct
756 TRITNL(36,0,1)= .0000000824 !.0000000949 !.0000002315
757 TRITNL(36,1,1)= .0000001153 !.0000001208 !.0000003005
758 TRITNL(36,2,1)= .0000002406 !.0000002062 !.0000005141
759 TRITNL(36,3,1)= .0000003828 !.0000002898 !.0000008021
760! 5 fct
761 TRITNL(36,0,2)= .0000000999 !.0000001062 !.0000002975
762 TRITNL(36,1,2)= .0000001420 !.0000001418 !.0000003608
763 TRITNL(36,2,2)= .0000002780 !.0000002480 !.0000005945
764 TRITNL(36,3,2)= .0000004353 !.0000003102 !.0000008260
765! 10 fct
766 TRITNL(36,0,3)= .0000001100 !.0000001300 !.0000003060
767 TRITNL(36,1,3)= .0000001600 !.0000001541 !.0000003820
768 TRITNL(36,2,3)= .0000003101 !.0000002622 !.0000006540
769 TRITNL(36,3,3)= .0000004980 !.0000003552 !.0000009000
770! LAW 42
771 SHTNL(42,0,1) = .0000002912 !.0000001981 !.0000005524
772 SHTNL(42,1,1) = .0000002912 !.0000001981 !.0000005524
773 SHTNL(42,2,1) = .0000007362 !.0000004928 !.0000013820
774 SHTNL(42,3,1) = .0000011840 !.0000008058 !.0000022760
775! LAW 42 M = 1
776 SHTNL(42,0,2) = .0000004000 !.0000002666 !.0000015100
777 SHTNL(42,1,2) = .0000004000 !.0000002666 !.0000015100
778 SHTNL(42,2,2) = .0000010475 !.0000006981 !.0000045280
779 SHTNL(42,3,2) = .0000017115 !.0000011410 !.0000074640
780! LAW 42 M = 2
781 SHTNL(42,0,3) = .0000004210 !.0000002796 !.0000017720
782 SHTNL(42,1,3) = .0000004210 !.0000002796 !.0000017720
783 SHTNL(42,2,3) = .0000011280 !.0000007499 !.0000053370
784 SHTNL(42,3,3) = .0000018560 !.0000012640 !.0000088600
785! LAW 42
786 TRITNL(42,0,1)= .0000002710 !.0000001907 !.0000005586
787 TRITNL(42,1,1)= .0000002710 !.0000001907 !.0000005586
788 TRITNL(42,2,1)= .0000007088 !.0000004940 !.0000015540
789 TRITNL(42,3,1)= .0000011500 !.0000007806 !.0000022710
790! LAW 42 M = 1
791 TRITNL(42,0,2)= .0000003832 !.0000002720 !.0000015000
792 TRITNL(42,1,2)= .0000003832 !.0000002720 !.0000015000
793 TRITNL(42,2,2)= .0000010340 !.0000007114 !.0000044580
794 TRITNL(42,3,2)= .0000016890 !.0000011630 !.0000073400
795! LAW 42 M = 2
796 TRITNL(42,0,3)= .0000004018 !.0000002772 !.0000017200
797 TRITNL(42,1,3)= .0000004018 !.0000002772 !.0000017200
798 TRITNL(42,2,3)= .0000011170 !.0000007621 !.0000052700
799 TRITNL(42,3,3)= .0000018640 !.0000012440 !.0000087000
800! LAW 43
801 SHTNL(43,0,1) = .0000001440 !.0000000958 !.0000002796
802 SHTNL(43,1,1) = .0000001449 !.0000000970 !.0000002918
803 SHTNL(43,2,1) = .0000002896 !.0000001942 !.0000006119
804 SHTNL(43,3,1) = .0000004214 !.0000002907 !.0000009136
805C
806 TRITNL(43,0,1)= .0000001350 !.0000000954 !.0000002856
807 TRITNL(43,1,1)= .0000001350 !.0000000945 !.0000002743
808 TRITNL(43,2,1)= .0000002700 !.0000001911 !.0000006279
809 TRITNL(43,3,1)= .0000004050 !.0000002852 !.0000008749
810! LAW 44
811 SHTNL(44,0,1) = .0000001390 !.0000000940 !.0000002640
812 SHTNL(44,1,1) = .0000001390 !.0000000940 !.0000002640
813 SHTNL(44,2,1) = .0000002716 !.0000001841 !.0000005714
814 SHTNL(44,3,1) = .0000004043 !.0000002761 !.0000008505
815C
816 TRITNL(44,0,1)= .0000001128 !.0000000891 !.0000002802
817 TRITNL(44,1,1)= .0000001128 !.0000000891 !.0000002802
818 TRITNL(44,2,1)= .0000002328 !.0000001766 !.0000005076
819 TRITNL(44,3,1)= .0000003504 !.0000002572 !.0000007598
820! LAW 48
821 SHTNL(48,0,1) = .0000001529 !.0000001037 !.0000002811
822 SHTNL(48,1,1) = .0000001529 !.0000001037 !.0000002811
823 SHTNL(48,2,1) = .0000003086 !.0000002113 !.0000005285
824 SHTNL(48,3,1) = .0000004696 !.0000003253 !.0000007925
825C
826 TRITNL(48,0,1)= .0000001258 !.0000000985 !.0000002764
827 TRITNL(48,1,1)= .0000001258 !.0000000985 !.0000002764
828 TRITNL(48,2,1)= .0000002718 !.0000002039 !.0000005447
829 TRITNL(48,3,1)= .0000004167 !.0000003064 !.0000007677
830! LAW 52
831 SHTNL(52,0,1) = .0000001398 !.0000000930 !.0000002313
832 SHTNL(52,1,1) = .0000001398 !.0000000930 !.0000002313
833 SHTNL(52,2,1) = .0000002699 !.0000001841 !.0000004775
834 SHTNL(52,3,1) = .0000003992 !.0000002788 !.0000007161
835C
836 TRITNL(52,0,1)= .0000001154 !.0000000910 !.0000002422
837 TRITNL(52,1,1)= .0000001154 !.0000000910 !.0000002422
838 TRITNL(52,2,1)= .0000002415 !.0000001797 !.0000004612
839 TRITNL(52,3,1)= .0000003696 !.0000002650 !.0000007237
840! LAW 57
841 SHTNL(57,0,1) = .0000001965 !.0000001303 !.0000005039
842 SHTNL(57,1,1) = .0000001965 !.0000001303 !.0000005039
843 SHTNL(57,2,1) = .0000004390 !.0000002981 !.0000013330
844 SHTNL(57,3,1) = .0000006959 !.0000004671 !.0000021130
845C
846 TRITNL(57,0,1)= .0000001867 !.0000001377 !.0000005568
847 TRITNL(57,1,1)= .0000001867 !.0000001377 !.0000005568
848 TRITNL(57,2,1)= .0000004672 !.0000003261 !.0000014330
849 TRITNL(57,3,1)= .0000007421 !.0000005104 !.0000022840
850! LAW 58
851! No visc (FLAGUL==0)
852 SHTNL(58,0,1) = .0000003516 !.0000002025 !.0000005435
853 SHTNL(58,1,1) = .0000003516 !.0000003516 !.0000006044
854 SHTNL(58,2,1) = .0000010550 !.0000007109 !.0000017910
855 SHTNL(58,3,1) = .0000017070 !.0000011530 !.0000028510
856C
857 TRITNL(58,0,1)= .0000003400 !.0000002113 !.0000005733
858 TRITNL(58,1,1)= .0000003400 !.0000002117 !.0000005733
859 TRITNL(58,2,1)= .0000009800 !.0000006316 !.0000015320
860 TRITNL(58,3,1)= .0000016500 !.0000009908 !.0000023640
861! LAW 60
862 SHTNL(60,0,1) = .0000001525 !.0000001044 !.0000003498
863 SHTNL(60,1,1) = .0000001744 !.0000001173 !.0000003992
864 SHTNL(60,2,1) = .0000003723 !.0000002531 !.0000008897
865 SHTNL(60,3,1) = .0000005732 !.0000003902 !.0000013350
866!
867 TRITNL(60,0,1)= .0000001598 !.0000001209 !.0000003850
868 TRITNL(60,1,1)= .0000001660 !.0000001273 !.0000003997
869 TRITNL(60,2,1)= .0000003935 !.0000002837 !.0000009103
870 TRITNL(60,3,1)= .0000006257 !.0000004377 !.0000014110
871! LAW 62
872 ! M = 0
873 SHTNL(62,0,1) = .0000002349 !.0000001553 !.0000008089
874 SHTNL(62,1,1) = .0000002349 !.0000001553 !.0000008089
875 SHTNL(62,2,1) = .0000005647 !.0000003720 !.0000021360
876 SHTNL(62,3,1) = .0000009094 !.0000005925 !.0000034950
877 ! M = 1
878 SHTNL(62,0,2) = .0000005137 !.0000003347 !.0000014640
879 SHTNL(62,1,2) = .0000005137 !.0000003347 !.0000014640
880 SHTNL(62,2,2) = .0000014080 !.0000009365 !.0000041200
881 SHTNL(62,3,2) = .0000023170 !.0000015170 !.0000067100
882 ! M = 2
883 SHTNL(62,0,3) = .0000007152 !.0000004757 !.0000017390
884 SHTNL(62,1,3) = .0000007152 !.0000004757 !.0000017390
885 SHTNL(62,2,3) = .0000020190 !.0000013210 !.0000048970
886 SHTNL(62,3,3) = .0000033370 !.0000022390 !.0000079110
887C
888 ! M = 0
889 TRITNL(62,0,1)= .0000002153 !.0000001570 !.0000008167
890 TRITNL(62,1,1)= .0000002153 !.0000001570 !.0000008167
891 TRITNL(62,2,1)= .0000005425 !.0000003764 !.0000021660
892 TRITNL(62,3,1)= .0000008716 !.0000005973 !.0000034760
893 ! M = 1
894 TRITNL(62,0,2)= .0000004882 !.0000003355 !.0000014950
895 TRITNL(62,1,2)= .0000004882 !.0000003355 !.0000014950
896 TRITNL(62,2,2)= .0000013850 !.0000009285 !.0000040800
897 TRITNL(62,3,2)= .0000022810 !.0000015420 !.0000067340
898 ! M = 2
899 TRITNL(62,0,3)= .0000006940 !.0000004649 !.0000017140
900 TRITNL(62,1,3)= .0000006940 !.0000004649 !.0000017140
901 TRITNL(62,2,3)= .0000020020 !.0000013230 !.0000048680
902 TRITNL(62,3,3)= .0000033140 !.0000022210 !.0000079480
903! LAW 63
904 SHTNL(63,0,1) = .0000001307 !.0000000872 !.0000002904
905 SHTNL(63,1,1) = .0000001307 !.0000000872 !.0000002904
906 SHTNL(63,2,1) = .0000002469 !.0000001682 !.0000006582
907 SHTNL(63,3,1) = .0000003666 !.0000002467 !.0000009983
908C
909 TRITNL(63,0,1)= .0000001053 !.0000000826 !.0000002766
910 TRITNL(63,1,1)= .0000001053 !.0000000826 !.0000002766
911 TRITNL(63,2,1)= .0000002133 !.0000001578 !.0000005871
912 TRITNL(63,3,1)= .0000003204 !.0000002313 !.0000008944
913! LAW 64
914 SHTNL(64,0,1) = .0000001496 !.0000001006 !.0000003364
915 SHTNL(64,1,1) = .0000001496 !.0000001006 !.0000003364
916 SHTNL(64,2,1) = .0000003055 !.0000002075 !.0000007866
917 SHTNL(64,3,1) = .0000004623 !.0000003197 !.0000011920
918C
919 TRITNL(64,0,1)= .0000001231 !.0000000966 !.0000003426
920 TRITNL(64,1,1)= .0000001231 !.0000000966 !.0000003426
921 TRITNL(64,2,1)= .0000002664 !.0000001997 !.0000007550
922 TRITNL(64,3,1)= .0000004126 !.0000003000 !.0000011970
923! LAW 65
924 SHTNL(65,0,1) = .0000001840 !.0000001242 !.0000003350
925 SHTNL(65,1,1) = .0000001840 !.0000001242 !.0000003350
926 SHTNL(65,2,1) = .0000003948 !.0000002677 !.0000007280
927 SHTNL(65,3,1) = .0000006070 !.0000004125 !.0000011240
928C
929 TRITNL(65,0,1)= .0000001563 !.0000001220 !.0000003324
930 TRITNL(65,1,1)= .0000001563 !.0000001220 !.0000003324
931 TRITNL(65,2,1)= .0000003650 !.0000002652 !.0000007364
932 TRITNL(65,3,1)= .0000005787 !.0000004097 !.0000011200
933! LAW 66
934 SHTNL(66,0,1) = .0000001628 !.0000001094 !.0000003310
935 SHTNL(66,1,1) = .0000001628 !.0000001094 !.0000003310
936 SHTNL(66,2,1) = .0000003334 !.0000002274 !.0000007200
937 SHTNL(66,3,1) = .0000005119 !.0000003468 !.0000011010
938C
939 TRITNL(66,0,1)= .0000001333 !.0000001057 !.0000003725
940 TRITNL(66,1,1)= .0000001333 !.0000001057 !.0000003666
941 TRITNL(66,2,1)= .0000002982 !.0000002192 !.0000007924
942 TRITNL(66,3,1)= .0000004634 !.0000003287 !.0000011650
943! LAW 69
944 ! M = 0 (no visc)
945 SHTNL(69,0,1) = .0000003985 !.0000002644 !.0000013330
946 SHTNL(69,1,1) = .0000003985 !.0000002644 !.0000013330
947 SHTNL(69,2,1) = .0000010580 !.0000007107 !.0000032970
948 SHTNL(69,3,1) = .0000017370 !.0000011590 !.0000054280
949 ! M = 1 (visc)
950 SHTNL(69,0,2) = .0000005288 !.0000003489 !.0000014210
951 SHTNL(69,1,2) = .0000005288 !.0000003489 !.0000014210
952 SHTNL(69,2,2) = .0000014430 !.0000009592 !.0000039620
953 SHTNL(69,3,2) = .0000023900 !.0000015630 !.0000064560
954 ! M = 2 (visc)
955 SHTNL(69,0,3) = .0000006553 !.0000004397 !.0000016580
956 SHTNL(69,1,3) = .0000006553 !.0000004397 !.0000016580
957 SHTNL(69,2,3) = .0000018340 !.0000012050 !.0000045270
958 SHTNL(69,3,3) = .0000030470 !.0000019880 !.0000074700
959C
960 ! M = 0 (no visc)
961 TRITNL(69,0,1)= .0000003792 !.0000002656 !.0000010400
962 TRITNL(69,1,1)= .0000003792 !.0000002656 !.0000010400
963 TRITNL(69,2,1)= .0000010490 !.0000007163 !.0000027280
964 TRITNL(69,3,1)= .0000017070 !.0000011470 !.0000044630
965 ! M = 1
966 TRITNL(69,0,2)= .0000005049 !.0000003447 !.0000012490
967 TRITNL(69,1,2)= .0000005049 !.0000003447 !.0000012490
968 TRITNL(69,2,2)= .0000014300 !.0000009625 !.0000035000
969 TRITNL(69,3,2)= .0000023740 !.0000015500 !.0000055580
970 ! M = 2
971 TRITNL(69,0,3)= .0000006319 !.0000004448 !.0000014430
972 TRITNL(69,1,3)= .0000006319 !.0000004448 !.0000014430
973 TRITNL(69,2,3)= .0000018180 !.0000012030 !.0000040510
974 TRITNL(69,3,3)= .0000030050 !.0000020200 !.0000065830
975! LAW 71 (incompatible npt=0)
976 SHTNL(71,0,1) = .0000002779 !.0000001823 !.0000005878
977 SHTNL(71,1,1) = .0000002779 !.0000001823 !.0000005878
978 SHTNL(71,2,1) = .0000006921 !.0000004629 !.0000014690
979 SHTNL(71,3,1) = .0000011280 !.0000007471 !.0000024250
980C
981 TRITNL(71,0,1)= .0000002571 !.0000001860 !.0000006319
982 TRITNL(71,1,1)= .0000002571 !.0000001860 !.0000006919
983 TRITNL(71,2,1)= .0000006719 !.0000004627 !.0000015420
984 TRITNL(71,3,1)= .0000010950 !.0000007391 !.0000024620
985! LAW 72
986 SHTNL(72,0,1) = .0000002082 !.0000001390 !.0000003671
987 SHTNL(72,1,1) = .0000002082 !.0000001390 !.0000003671
988 SHTNL(72,2,1) = .0000004876 !.0000003235 !.0000007838
989 SHTNL(72,3,1) = .0000007582 !.0000005036 !.0000012730
990C
991 TRITNL(72,0,1)= .0000001787 !.0000001334 !.00000s03591
992 TRITNL(72,1,1)= .0000001787 !.0000001334 !.0000003591
993 TRITNL(72,2,1)= .0000004340 !.0000003038 !.0000007633
994 TRITNL(72,3,1)= .0000006868 !.0000004765 !.0000011460
995! LAW 73
996 SHTNL(73,0,1) = .0000001921 !.0000001248 !.0000003456
997 SHTNL(73,1,1) = .0000001921 !.0000001248 !.0000003689
998 SHTNL(73,2,1) = .0000004181 !.0000002821 !.0000007752
999 SHTNL(73,3,1) = .0000006513 !.0000004450 !.0000012010
1000C
1001 TRITNL(73,0,1)= .0000001655 !.0000001245 !.0000003618
1002 TRITNL(73,1,1)= .0000001655 !.0000001245 !.0000003371
1003 TRITNL(73,2,1)= .0000003938 !.0000002802 !.0000007637
1004 TRITNL(73,3,1)= .0000006235 !.0000004278 !.0000012300
1005! LAW 76
1006 SHTNL(76,0,1) = .0000003014 !.0000002022 !.0000007383
1007 SHTNL(76,1,1) = .0000003014 !.0000002022 !.0000007476
1008 SHTNL(76,2,1) = .0000007581 !.0000005089 !.0000018800
1009 SHTNL(76,3,1) = .0000012190 !.0000008119 !.0000029980
1010C
1011 TRITNL(76,0,1)= .0000002751 !.0000001994 !.0000007008
1012 TRITNL(76,1,1)= .0000002751 !.0000001994 !.0000007008
1013 TRITNL(76,2,1)= .0000007339 !.0000005036 !.0000018360
1014 TRITNL(76,3,1)= .0000011940 !.0000008260 !.0000029000
1015! LAW 78
1016 SHTNL(78,0,1) = .0000001539 !.0000001045 !.0000002921
1017 SHTNL(78,1,1) = .0000001539 !.0000001045 !.0000002921
1018 SHTNL(78,2,1) = .0000003119 !.0000002158 !.0000005980
1019 SHTNL(78,3,1) = .0000004750 !.0000003308 !.0000009682
1020C
1021 TRITNL(78,0,1)= .0000001307 !.0000001017 !.0000002936
1022 TRITNL(78,1,1)= .0000001307 !.0000001017 !.0000002939
1023 TRITNL(78,2,1)= .0000002885 !.0000002104 !.0000005670
1024 TRITNL(78,3,1)= .0000004450 !.0000003165 !.0000008430
1025! LAW 80
1026 SHTNL(80,0,1) = .0000010640 !.0000007081 !.0000015110
1027 SHTNL(80,1,1) = .0000010640 !.0000007081 !.0000014400
1028 SHTNL(80,2,1) = .0000030400 !.0000020610 !.0000040060
1029 SHTNL(80,3,1) = .0000045130 !.0000029980 !.0000065840
1030C
1031 TRITNL(80,0,1)= .0000001015 !.0000006949 !.0000013660
1032 TRITNL(80,1,1)= .0000001015 !.0000006949 !.0000013660
1033 TRITNL(80,2,1)= .0000002963 !.0000019830 !.0000039870
1034 TRITNL(80,3,1)= .0000004913 !.0000032500 !.0000063540
1035! LAW 82 (incompatible NPT=0)
1036! N = 1
1037 SHTNL(82,0,1) = .0000002316 !.0000001536 !.0000006991
1038 SHTNL(82,1,1) = .0000002316 !.0000001536 !.0000007022
1039 SHTNL(82,2,1) = .0000005547 !.0000003693 !.0000017930
1040 SHTNL(82,3,1) = .0000008899 !.0000005914 !.0000028300
1041!
1042 TRITNL(82,0,1)= .0000002102 !.0000001558 !.0000006886
1043 TRITNL(82,1,1)= .0000002108 !.0000001558 !.0000007115
1044 TRITNL(82,2,1)= .0000005342 !.0000003788 !.0000017680
1045 TRITNL(82,3,1)= .0000008588 !.0000005881 !.0000028740
1046! N = 2
1047 SHTNL(82,0,2) = .0000003484 !.0000002314 !.0000010470
1048 SHTNL(82,1,2) = .0000003484 !.0000002314 !.0000010430
1049 SHTNL(82,2,2) = .0000009013 !.0000006029 !.0000027830
1050 SHTNL(82,3,2) = .0000014760 !.0000009733 !.0000045710
1051!
1052 TRITNL(82,0,2)= .0000003256 !.0000002278 !.0000010690
1053 TRITNL(82,1,2)= .0000003256 !.0000002278 !.0000010330
1054 TRITNL(82,2,2)= .0000008872 !.0000006062 !.0000028430
1055 TRITNL(82,3,2)= .0000014420 !.0000009569 !.0000046350
1056! N = 3
1057 SHTNL(82,0,3) = .0000004674 !.0000003158 !.0000013880
1058 SHTNL(82,1,3) = .0000004674 !.0000003158 !.0000014040
1059 SHTNL(82,2,3) = .0000012710 !.0000008427 !.0000038590
1060 SHTNL(82,3,3) = .0000020930 !.0000013890 !.0000063230
1061!
1062 TRITNL(82,0,3)= .0000004538 !.0000003084 !.0000014340
1063 TRITNL(82,1,3)= .0000004538 !.0000003084 !.0000014220
1064 TRITNL(82,2,3)= .0000012530 !.0000008498 !.0000038320
1065 TRITNL(82,3,3)= .0000020610 !.0000013620 !.0000063450
1066! LAW 87
1067 SHTNL(87,0,1) = .0000002036 !.0000001386 !.0000003592
1068 SHTNL(87,1,1) = .0000002036 !.0000001386 !.0000003592
1069 SHTNL(87,2,1) = .0000004685 !.0000003151 !.0000007907
1070 SHTNL(87,3,1) = .0000007389 !.0000004885 !.0000011520
1071!
1072 TRITNL(87,0,1)= .0000001826 !.0000001347 !.0000003754
1073 TRITNL(87,1,1)= .0000001826 !.0000001347 !.0000003874
1074 TRITNL(87,2,1)= .0000004412 !.0000003124 !.0000008403
1075 TRITNL(87,3,1)= .0000007025 !.0000004843 !.0000012850
1076! LAW 93
1077 SHTNL(93,0,1) = .0000001371 !.0000000928 !.0000002801
1078 SHTNL(93,1,1) = .0000001371 !.0000000928 !.0000002801
1079 SHTNL(93,2,1) = .0000002701 !.0000001849 !.0000005238
1080 SHTNL(93,3,1) = .0000004036 !.0000002761 !.0000007912
1081!
1082 TRITNL(93,0,1)= .0000001175 !.0000000922 !.0000002629
1083 TRITNL(93,1,1)= .0000001175 !.0000000922 !.0000002629
1084 TRITNL(93,2,1)= .0000002418 !.0000001835 !.0000005262
1085 TRITNL(93,3,1)= .0000003755 !.0000002740 !.0000008110
1086
1087! LAW 104
1088 ! Nice algo
1089 SHTNL(104,0,1) = .0000001034
1090 SHTNL(104,1,1) = .0000001034
1091 SHTNL(104,2,1) = .0000002044
1092 SHTNL(104,3,1) = .0000003111
1093 ! Newton algo
1094 SHTNL(104,0,2) = .0000001400
1095 SHTNL(104,1,2) = .0000001400
1096 SHTNL(104,2,2) = .0000002944
1097 SHTNL(104,3,2) = .0000004683
1098 ! Nice algo
1099 TRITNL(104,0,1) = .0000001034
1100 TRITNL(104,1,1) = .0000001034
1101 TRITNL(104,2,1) = .0000002044
1102 TRITNL(104,3,1) = .0000003111
1103 ! Newton algo
1104 TRITNL(104,0,2) = .0000001400
1105 TRITNL(104,1,2) = .0000001400
1106 TRITNL(104,2,2) = .0000002944
1107 TRITNL(104,3,2) = .0000004683
1108 ! iloc=1
1109 SHTNL_OPTION(1,1) = .0000000400
1110 ! iloc=2
1111 SHTNL_OPTION(1,2) = .0000000500
1112 ! iloc=3
1113 SHTNL_OPTION(1,3) = .0000000900
1114
1115C Temps de traitement de la partie elementaire sans la loi materielle fonction de IPARITH
1116 IF (IPARI0.EQ.1) THEN
1117! SHELT :
1118! 1-->5 --> /SHELL
1119! 6-->10 --> /SH_COMP or TYPE10
1120! 11-->15 --> /SH_SANDW or TYPE11
1121! 15-->20 --> /SH_ORTH or TYPE9
1122! 20-->25 --> /SH_FABR or TYPE16
1123! IDRILL = 0
1124! 1/6/11/16/21 --> Q4
1125! 2/7/12/17/22 --> QEPH
1126! 3/8/13/18/23 --> BATOZ
1127! IDRILL = 1
1128! none --> Q4
1129! 4/9 /14/19/24 --> QEPH
1130! 5/10/15/20/25 --> BATOZ
1131! TRIELT :
1132! 1 --> tri
1133! NLOCAL_OPTION
1134! 1 --> sol
1135! 2 --> shell
1136! /SHELL
1137 SHTELT(1) = .0000002670
1138 SHTELT(2) = .0000003380
1139 SHTELT(3) = .0000008070
1140 SHTELT(4) = .0000004730
1141 SHTELT(5) = .0000010390
1142! /SH_COMP
1143 SHTELT(6) = .0000002870
1144 SHTELT(7) = .0000003550
1145 SHTELT(8) = .0000008350
1146 SHTELT(9) = .0000004560
1147 SHTELT(10) = .0000010620
1148! /SH_SANDW
1149 SHTELT(11) = .0000003040
1150 SHTELT(12) = .0000004510
1151 SHTELT(13) = .0000009120
1152 SHTELT(14) = .0000005710
1153 SHTELT(15) = .0000003040
1154! /SH_ORTH
1155 SHTELT(16) = .0000002820
1156 SHTELT(17) = .0000003820
1157 SHTELT(18) = .0000008800
1158 SHTELT(19) = .0000004960
1159 SHTELT(20) = .0000003040
1160! /SH_FABR
1161 SHTELT(21) = .0000003800
1162 SHTELT(22) = .0000004410
1163 SHTELT(23) = .0000009780
1164 SHTELT(24) = .0000004600
1165 SHTELT(25) = .0000016180
1166
1167! /TYPE51
1168 SHTELT(26) = .0000008400
1169 SHTELT(27) = .0000009840
1170 SHTELT(28) = .0000021800
1171 SHTELT(29) = .0000010040
1172 SHTELT(30) = .0000029300
1173
1174 TRITELT(1)= .0000002410
1175 NLOCAL_OPTION(2) = .0000000600
1176 ELSE
1177! /SHELL
1178 SHTELT(1) = .0000001695 !.0000004647
1179 SHTELT(2) = .0000002289 !.0000005736
1180 SHTELT(3) = .0000006200 !.0000012330
1181 SHTELT(4) = .0000003115
1182 SHTELT(5) = .0000008131
1183! /SH_COMP
1184 SHTELT(6) = .0000001868 !.0000005935
1185 SHTELT(7) = .0000002431 !.0000007462
1186 SHTELT(8) = .0000006430 !.0000013800
1187 SHTELT(9) = .0000003276
1188 SHTELT(10)= .0000008321
1189! /SH_SANDW
1190 SHTELT(11) = .0000002008 !.0000005224
1191 SHTELT(12) = .0000003229 !.0000007354
1192 SHTELT(13) = .0000007075 !.0000012600
1193 SHTELT(14) = .0000004230
1194 SHTELT(15) = .0000008941
1195! /SH_ORTH
1196 SHTELT(16) = .0000001827 !.0000003400
1197 SHTELT(17) = .0000002658 !.0000006300
1198 SHTELT(18) = .0000006808 !.0000013040
1199 SHTELT(19) = .0000003611
1200 SHTELT(20) = .0000008800
1201! /SH_FABR
1202 SHTELT(21) = .0000002641 !.0000005500
1203 SHTELT(22) = .0000003150 !.0000006440
1204 SHTELT(23) = .0000007625 !.0000012700
1205 SHTELT(24) = .0000003333
1206 SHTELT(25) = .0000012900
1207
1208! /TYPE51
1209 SHTELT(26) = .0000005800
1210 SHTELT(27) = .0000006900
1211 SHTELT(28) = .0000015570
1212 SHTELT(29) = .0000007000
1213 SHTELT(30) = .0000020710
1214
1215 TRITELT(1)= .0000001391 !.0000003721
1216 NLOCAL_OPTION(2) = .0000000600
1217
1218 ENDIF
1219
1220! il y a un facteur multiplicatif sur les coques de Batoz, car on rentre 4x dans cmain3.
1221
1222 BATOZMULT = 3.8
1223! /PRONY or /LPRONY overcost
1224 VISC_PRONY_COST = .0000000088
1225
1226! *******************************************************************
1227! Failure model
1228! *******************************************************************
1229! Johnson : 1
1230! Tbutcher : 2
1231! Wilkins : 3
1232! user1 : 4
1233! user2 : 5
1234! user3 : 6
1235! FLD : 7 (for shell only)
1236! Spalling : 8 (for solid only)
1237! Wierzbicki : 9
1238! Tensstrain : 10
1239! Energy : 11
1240! ??? : 12
1241! Chang : 13 (for shell only)
1242! Hashin : 14
1243! Yamada : 15
1244! Puck : 16
1245! XFEM_Tbutcher : 17
1246! Lad_dama : 18 (for solid only)
1247! Power_da : 19
1248! Connect : 20 (for solid only)
1249! XFEM_Johns : 21
1250! XFEM_FLD : 22
1251! TAB1 : 23
1252! Orthstra : 24 (for shell only)
1253! NXT : 25 (for shell only)
1254! Snconnect : 26 (for solid only)
1255! EMC : 27 (for solid only)
1256! Wind : 28 (for shell only)
1257! Sahraei1 : 29 (for solid only)
1258! Biquadratic : 30
1259! _
1260! SHELL : 1--> 0 ou 1pt |
1261! 2--> 3 pts |--- MLN < 28 or ==32
1262! 3--> 5 pts _|
1263! 4--> 0 ou 1pt |
1264! 5--> 3 pts |--- MLN >= 28 or /=32
1265! 6--> 5 pts _|
1266!
1267! SOLIDE : 1 ----> MLN <28 or ==49
1268! 2 ----> MLN >=28 or /=49
1269! --------------------------------------------------------
1270 DO I=1,MAXFAIL
1271 DO J=1,3
1272 RUPTURE_SOL(I,J) = 0.
1273 RUPTURE_SOL12(I,J) = 0.
1274 RUPTURE_TET4(I,J) = 0.
1275 RUPTURE_TET10(I,J) = 0.
1276 ENDDO
1277 DO J=1,6
1278 RUPTURE_SHELL(I,J) = 0.
1279 ENDDO
1280 ENDDO
1281! --------------------
1282! Johnson
1283 ! MLN <28 or ==32
1284 RUPTURE_SHELL(1,1) = .00000002342 !.00000007517
1285 RUPTURE_SHELL(1,2) = .00000005640 !.00000019090
1286 RUPTURE_SHELL(1,3) = .00000008993 !.00000031460
1287
1288 ! MLN >=28 or /=32
1289 RUPTURE_SHELL(1,4) = .00000001638 !.00000003021
1290 RUPTURE_SHELL(1,5) = .00000003725 !.00000007567
1291 RUPTURE_SHELL(1,6) = .00000005660 !.00000011620
1292
1293 ! MLN <28 or ==49
1294 RUPTURE_SOL(1,1) = .00000004261 !.00000008839
1295 RUPTURE_SOL12(1,1) = .00000012740 !.00000017710
1296 RUPTURE_TET4(1,1) = .00000004250 !.00000008021
1297 RUPTURE_TET10(1,1) = .00000015910 !.00000033420
1298 ! MLN >=28 or /=49
1299 RUPTURE_SOL(1,2) = .00000002573 !.00000002696
1300 RUPTURE_SOL12(1,2) = .00000015790 !.00000015790
1301 RUPTURE_TET4(1,2) = .00000002443 !.00000002663
1302 RUPTURE_TET10(1,2) = .00000009469 !.00000010360
1303! --------------------
1304! Tbutcher
1305 ! MLN <28 or ==49
1306 RUPTURE_SHELL(2,1) = .00000001895 !.00000004873
1307 RUPTURE_SHELL(2,2) = .00000004270 !.00000010840
1308 RUPTURE_SHELL(2,3) = .00000006646 !.00000016340
1309
1310 RUPTURE_SHELL(2,4) = .00000001827 !.00000003994
1311 RUPTURE_SHELL(2,5) = .00000004190 !.00000010130
1312 RUPTURE_SHELL(2,6) = .00000006501 !.00000015690
1313
1314 ! MLN <28 or ==49
1315 RUPTURE_SOL(2,1) = .00000004250 !.00000004718
1316 RUPTURE_SOL12(2,1) = .00000013750 !.00000031390
1317 RUPTURE_TET4(2,1) = .00000004034 !.00000010090
1318 RUPTURE_TET10(2,1) = .00000016070 !.00000040130
1319
1320 ! MLN >=28 or /=49
1321 RUPTURE_SOL(2,2) = .00000002613 !.00000004949
1322 RUPTURE_SOL12(2,2) = .00000017280 !.00000030290
1323 RUPTURE_TET4(2,2) = .00000002267 !.00000004449
1324 RUPTURE_TET10(2,2) = .00000009430 !.00000016800
1325! --------------------
1326! Wilkins
1327 ! MLN <28 or ==49
1328 RUPTURE_SHELL(3,1) = .00000006864 !.00000011660
1329 RUPTURE_SHELL(3,2) = .00000019790 !.00000032760
1330 RUPTURE_SHELL(3,3) = .00000033320 !.00000052680
1331
1332 RUPTURE_SHELL(3,4) = .00000003115 !.00000011410
1333 RUPTURE_SHELL(3,5) = .00000008054 !.00000032760
1334 RUPTURE_SHELL(3,6) = .00000013990 !.00000052910
1335
1336 ! MLN <28 or ==49
1337 RUPTURE_SOL(3,1) = .00000013710 !.00000023210
1338 RUPTURE_SOL12(3,1) = .00000091220 !.00000131700
1339 RUPTURE_TET4(3,1) = .00000008845 !.00000016440
1340 RUPTURE_TET10(3,1) = .00000043350 !.00000075060
1341
1342 ! MLN >=28 or /=49
1343 RUPTURE_SOL(3,2) = .00000011900 !.00000016870
1344 RUPTURE_SOL12(3,2) = .00000100800 !.00000129800
1345 RUPTURE_TET4(3,2) = .00000007117 !.00000011060
1346 RUPTURE_TET10(3,2) = .00000034650 !.00000051000
1347! --------------------
1348! FLD (only for shell)
1349 ! MLN <28 or ==32
1350 RUPTURE_SHELL(7,1) = .00000003036 !.00000010900
1351 RUPTURE_SHELL(7,2) = .00000007421 !.00000021590
1352 RUPTURE_SHELL(7,3) = .00000011910 !.00000031890
1353
1354 ! MLN >=28 or /=32
1355 RUPTURE_SHELL(7,4) = .00000003115 !.00000007162
1356 RUPTURE_SHELL(7,5) = .00000008054 !.00000019550
1357 RUPTURE_SHELL(7,6) = .00000013990 !.00000031220
1358! --------------------
1359! Spalling (only for solid)
1360
1361 ! MLN <28 or ==49
1362 RUPTURE_SOL(8,1) = .00000003544 !.00000009428
1363 RUPTURE_SOL12(8,1) = .00000008080 !.00000025070
1364 RUPTURE_TET4(8,1) = .00000003630 !.00000009402
1365 RUPTURE_TET10(8,1) = .00000013920 !.00000038110
1366
1367 ! MLN >=28 or /=49
1368 RUPTURE_SOL(8,2) = .00000002000 !.00000003871
1369 RUPTURE_SOL12(8,2) = .00000012370 !.00000023730
1370 RUPTURE_TET4(8,2) = .00000001807 !.00000003716
1371 RUPTURE_TET10(8,2) = .00000007336 !.00000014820
1372! --------------------
1373! Wierzbicki
1374 ! MLN <28 or ==49
1375 RUPTURE_SHELL(9,1) = .00000016750 !.00000025370
1376 RUPTURE_SHELL(9,2) = .00000049190 !.00000071910
1377 RUPTURE_SHELL(9,3) = .00000082040 !.00000118100
1378
1379 RUPTURE_SHELL(9,4) = .00000016670 !.00000024610
1380 RUPTURE_SHELL(9,5) = .00000049650 !.00000071940
1381 RUPTURE_SHELL(9,6) = .00000080950 !.00000118500
1382
1383 ! MLN <28 or ==49
1384 RUPTURE_SOL(9,1) = .00000003236 !.00000008443
1385 RUPTURE_SOL12(9,1) = .00000005825 !.00000011190
1386 RUPTURE_TET4(9,1) = .00000003239 !.00000008236
1387 RUPTURE_TET10(9,1) = .00000012800 !.00000032530
1388
1389 ! MLN >=28 or /=49
1390 RUPTURE_SOL(9,2) = .00000001600 !.00000002591
1391 RUPTURE_SOL12(9,2) = .00000009595 !.00000016310
1392 RUPTURE_TET4(9,2) = .00000001494 !.00000002483
1393 RUPTURE_TET10(9,2) = .00000005972 !.00000008821
1394! --------------------
1395! Tensstrain
1396 ! MLN <28 or ==32
1397 RUPTURE_SHELL(10,1) = .00000002166 !.00000009068
1398 RUPTURE_SHELL(10,2) = .00000005092 !.00000016430
1399 RUPTURE_SHELL(10,3) = .00000008145 !.00000023120
1400
1401 ! MLN >=28 or /=32
1402 RUPTURE_SHELL(10,4) = .00000002328 !.00000004473
1403 RUPTURE_SHELL(10,5) = .00000005666 !.00000012650
1404 RUPTURE_SHELL(10,6) = .00000009416 !.00000019170
1405
1406 ! MLN <28 or ==49
1407 RUPTURE_SOL(10,1) = .00000006278 !.00000013990
1408 RUPTURE_SOL12(10,1) = .00000012500 !.00000025600
1409 RUPTURE_TET4(10,1) = .00000005218 !.00000010960
1410 RUPTURE_TET10(10,1) = .00000022200 !.00000047680
1411
1412 ! MLN >=28 or /=49
1413 RUPTURE_SOL(10,2) = .00000004650 !.00000008052
1414 RUPTURE_SOL12(10,2) = .00000035060 !.00000056320
1415 RUPTURE_TET4(10,2) = .00000002950 !.00000005430
1416 RUPTURE_TET10(10,2) = .00000015020 !.00000025250
1417! --------------------
1418! Energy
1419 ! MLN <28 or ==49
1420 RUPTURE_SHELL(11,1) = .00000001914 !.00000005295
1421 RUPTURE_SHELL(11,2) = .00000004461 !.00000011720
1422 RUPTURE_SHELL(11,3) = .00000006982 !.00000017760
1423
1424 RUPTURE_SHELL(11,4) = .00000002371 !.00000003770
1425 RUPTURE_SHELL(11,5) = .00000005454 !.00000010100
1426 RUPTURE_SHELL(11,6) = .00000008598 !.00000015840
1427
1428 ! MLN <28 or ==49
1429 RUPTURE_SOL(11,1) = .00000003780 !.00000009413
1430 RUPTURE_SOL12(11,1) = .00000003700 !.00000009015
1431 RUPTURE_TET4(11,1) = .00000004242 !.00000009164
1432 RUPTURE_TET10(11,1) = .00000015180 !.00000037960
1433
1434 ! MLN >=28 or /=49
1435 RUPTURE_SOL(11,2) = .00000002200 !.00000003883
1436 RUPTURE_SOL12(11,2) = .00000014500 !.00000025910
1437 RUPTURE_TET4(11,2) = .00000002011 !.00000003753
1438 RUPTURE_TET10(11,2) = .00000008094 !.00000014030
1439! --------------------
1440! Chang (only for shell)
1441 ! MLN <28 or ==49
1442 RUPTURE_SHELL(13,1) = .00000001598 !.00000004822
1443 RUPTURE_SHELL(13,2) = .00000003435 !.00000010740
1444 RUPTURE_SHELL(13,3) = .00000005443 !.00000016000
1445 ! MLN >=28 or /=49
1446 RUPTURE_SHELL(13,4) = .00000001948 !.00000003770
1447 RUPTURE_SHELL(13,5) = .00000004518 !.00000010100
1448 RUPTURE_SHELL(13,6) = .00000006896 !.00000015840
1449! --------------------
1450! Hashin
1451 ! MLN <28 or ==49
1452 RUPTURE_SHELL(14,1) = .00000002500 !.00000007989
1453 RUPTURE_SHELL(14,2) = .00000007500 !.00000025240
1454 RUPTURE_SHELL(14,3) = .00000013000 !.00000040300
1455
1456 RUPTURE_SHELL(14,4) = .00000002500 !.00000007989
1457 RUPTURE_SHELL(14,5) = .00000007500 !.00000025240
1458 RUPTURE_SHELL(14,6) = .00000013000 !.00000040300
1459
1460 ! MLN <28 or ==49
1461 RUPTURE_SOL(14,1) = .00000007500 !.00000013480
1462 RUPTURE_SOL12(14,1) = .00000005400 !.00000107400
1463 RUPTURE_TET4(14,1) = .00000006250 !.00000012500
1464 RUPTURE_TET10(14,1) = .00000023010 !.00000054770
1465
1466 ! MLN >=28 or /=49
1467 RUPTURE_SOL(14,2) = .00000002000 !.00000001289
1468 RUPTURE_SOL12(14,2) = .00000010000 !.00000010990
1469 RUPTURE_TET4(14,2) = .00000001500 !.00000001302
1470 RUPTURE_TET10(14,2) = .00000006000 !.00000005823
1471! --------------------
1472! Puck (compatible MLN<28)
1473 ! MLN <28 or ==49
1474 RUPTURE_SHELL(16,1) = .00000002028 !.00000007462
1475 RUPTURE_SHELL(16,2) = .00000004716 !.00000018470
1476 RUPTURE_SHELL(16,3) = .00000007455 !.00000028840
1477
1478 RUPTURE_SHELL(16,4) = .00000002229 !.00000006606
1479 RUPTURE_SHELL(16,5) = .00000005406 !.00000017190
1480 RUPTURE_SHELL(16,6) = .00000008620 !.00000027600
1481
1482 ! MLN <28 or ==49
1483 RUPTURE_SOL(16,1) = .00000004857 !.00000013060
1484 RUPTURE_SOL12(16,1) = .00000003453 !.00000008806
1485 RUPTURE_TET4(16,1) = .00000004501 !.00000012170
1486 RUPTURE_TET10(16,1) = .00000018820 !.00000042400
1487
1488 RUPTURE_SOL(16,2) = .00000002500 !.00000005419
1489 RUPTURE_SOL12(16,2) = .00000011000 !.00000007646
1490 RUPTURE_TET4(16,2) = .00000001240 !.00000001338
1491 RUPTURE_TET10(16,2) = .00000004832 !.00000004460
1492! --------------------
1493! Lad_Dama (compatible MLN<28) (only for solid)
1494 ! MLN <28 or ==49
1495 RUPTURE_SOL(18,1) = .00000005564 !.00000016890
1496 RUPTURE_SOL12(18,1) = .00000003508 !.00000008820
1497 RUPTURE_TET4(18,1) = .00000005428 !.00000013980
1498 RUPTURE_TET10(18,1) = .00000021920 !.00000058630
1499
1500 RUPTURE_SOL(18,2) = .00000070000 !.00000016890
1501 RUPTURE_SOL12(18,2) = .00000250000 !.00000008820
1502 RUPTURE_TET4(18,2) = .00000070000 !.00000013980
1503 RUPTURE_TET10(18,2) = .00000280000 !.00000058630
1504! --------------------
1505! Connect (compatible MLN==59 et 83)
1506 ! MLN==59
1507 RUPTURE_SOL(20,1) = .00000025000 !.00000054140
1508 RUPTURE_SOL12(20,1) = .00000025000 !.00000054140
1509 RUPTURE_TET4(20,1) = .00000025000 !.00000054140
1510 RUPTURE_TET10(20,1) = .00000025000 !.00000054140
1511 ! MLN==83
1512 RUPTURE_SOL(20,2) = .00000025000 !.00000054140
1513 RUPTURE_SOL12(20,2) = .00000025000 !.00000054140
1514 RUPTURE_TET4(20,2) = .00000025000 !.00000054140
1515 RUPTURE_TET10(20,2) = .00000025000 !.00000054140
1516! --------------------
1517! Tab1
1518 ! MLN <28 or ==49
1519 RUPTURE_SHELL(23,1) = .00000010690 !.00000021020
1520 RUPTURE_SHELL(23,2) = .00000031560 !.00000060850
1521 RUPTURE_SHELL(23,3) = .00000052190 !.00000098850
1522
1523 RUPTURE_SHELL(23,4) = .00000011350 !.00000021590
1524 RUPTURE_SHELL(23,5) = .00000031430 !.00000062580
1525 RUPTURE_SHELL(23,6) = .00000046800 !.0000009949
1526
1527 ! MLN <28 or ==49
1528 RUPTURE_SOL(23,1) = .00000017800 !.00000031120
1529 RUPTURE_SOL12(23,1) = .00000123900 !.00000192900
1530 RUPTURE_TET4(23,1) = .00000016280 !.00000026620
1531 RUPTURE_TET10(23,1) = .00000068920 !.00000118100
1532
1533 ! MLN >=28 or /=49
1534 RUPTURE_SOL(23,2) = .00000016170 !.00000025940
1535 RUPTURE_SOL12(23,2) = .00000142600 !.00000198600
1536 RUPTURE_TET4(23,2) = .00000014240 !.00000022330
1537 RUPTURE_TET10(23,2) = .00000061180 !.00000096400
1538! --------------------
1539! NXT (only for shell)
1540 ! MLN <28 or ==32
1541 RUPTURE_SHELL(25,1) = .00000003909 !.00000014880
1542 RUPTURE_SHELL(25,2) = .00000010420 !.00000032860
1543 RUPTURE_SHELL(25,3) = .00000017070 !.00000051100
1544
1545 ! MLN >=28 or /=32
1546 RUPTURE_SHELL(25,4) = .00000003648 !.00000008550
1547 RUPTURE_SHELL(25,5) = .00000009550 !.00000024030
1548 RUPTURE_SHELL(25,6) = .00000015980 !.00000038060
1549! --------------------
1550! Snconnect (compatible MLN==59 et 83)
1551 ! MLN==59
1552 RUPTURE_SOL(26,1) = .00000030000 !.00000065150
1553 RUPTURE_SOL12(26,1) = .00000030000 !.00000065150
1554 RUPTURE_TET4(26,1) = .00000030000 !.00000065150
1555 RUPTURE_TET10(26,1) = .00000030000 !.00000065150
1556 ! MLN==83
1557 RUPTURE_SOL(26,2) = .00000030000 !.00000065150
1558 RUPTURE_SOL12(26,2) = .00000030000 !.00000065150
1559 RUPTURE_TET4(26,2) = .00000030000 !.00000065150
1560 RUPTURE_TET10(26,2) = .00000030000 !.00000065150
1561! --------------------
1562! EMC (only for solid)
1563 ! MLN <28 or ==49
1564 RUPTURE_SOL(27,1) = .00000003106 !.00000007774
1565 RUPTURE_SOL12(27,1) = .00000004772 !.00000013030
1566 RUPTURE_TET4(27,1) = .00000003162 !.00000007539
1567 RUPTURE_TET10(27,1) = .00000012270 !.00000029830
1568
1569 ! MLN >=28 or /=49
1570 RUPTURE_SOL(27,2) = .00000005000 !.00000001913
1571 RUPTURE_SOL12(27,2) = .00000020000 !.00000010660
1572 RUPTURE_TET4(27,2) = .00000002000 !.00000001654
1573 RUPTURE_TET10(27,2) = .00000008000 !.00000006008
1574! --------------------
1575! Alter (only shell/tri)
1576 ! MLN <28 or ==49
1577 RUPTURE_SHELL(28,1) = .0000000647
1578 RUPTURE_SHELL(28,2) = .0000001670
1579 RUPTURE_SHELL(28,3) = .0000002749
1580
1581 RUPTURE_SHELL(28,4) = .0000000647
1582 RUPTURE_SHELL(28,5) = .0000001670
1583 RUPTURE_SHELL(28,6) = .0000002749
1584! --------------------
1585! Biquad
1586 ! MLN <28 or ==49
1587 RUPTURE_SHELL(30,1) = .00000001321 !.00000003069
1588 RUPTURE_SHELL(30,2) = .00000002521 !.00000005533
1589 RUPTURE_SHELL(30,3) = .00000003750 !.00000008252
1590
1591 RUPTURE_SHELL(30,4) = .00000001680 !.00000003133
1592 RUPTURE_SHELL(30,5) = .00000003637 !.00000006587
1593 RUPTURE_SHELL(30,6) = .00000005447 !.00000010300
1594
1595 ! MLN <28 or ==49
1596 RUPTURE_SOL(30,1) = .00000003124 !.00000008026
1597 RUPTURE_SOL12(30,1) = .00000004905 !.00000013200
1598 RUPTURE_TET4(30,1) = .00000003162 !.00000007523
1599 RUPTURE_TET10(30,1) = .00000012270 !.00000030150
1600
1601 RUPTURE_SOL(30,2) = .00000002500 !.00000001929
1602 RUPTURE_SOL12(30,2) = .00000008000 !.00000011930
1603 RUPTURE_TET4(30,2) = .00000002500 !.00000001793
1604 RUPTURE_TET10(30,2) = .00000007500 !.00000006588
1605! --------------------
cmain3
subroutine cmain3(timers, elbuf_str, jft, jlt, nft, iparg, nel, mtn, ipla, ithk, group_param, pm, geo, npf, tf, bufmat, ssp, rho, viscmx, dt1c, sigy, area, exx, eyy, exy, exz, eyz, kxx, kyy, kxy, nu, off, thk0, mat, pid, mat_elem, for, mom, gstr, failwave, fwave_el, thk, eint, iofc, g, a11, a12, vol0, indxdel, ngl, zcfac, shf, gs, epsd_pg, kfts, jhbe, alpe, dir_a, dir_b, igeo, ipm, ifailure, npg, fheat, tempel, die, jthe, iexpan, tempel0, ishplyxfem, ply_exx, ply_eyy, ply_exy, ply_exz, ply_eyz, ply_f, del_ply, th_iply, sig_iply, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, ng, table, ixfem, offi, sensors, a11_iply, elcrkini, dir1_crk, dir2_crk, aldt, idt_therm, theaccfact, ismstr, ir, is, nlay, npt, ixlay, ixel, isubstack, stack, f_def, itask, drape, varnl, nloc_dmg, indx_drape, thke, sedrape, numel_drape, dt, ncycle, snpc, stf, nxlaymax, idel7nok, userl_avail, maxfunc, varnl_npttot, sbufmat, sdir_a, sdir_b, for_g)
Definition cmain3.F:87
the
end diagonal values have been computed in the(sparse) matrix id.SOL
N
#define N
Definition example_user.c:10
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
mumps_ooc_common::fct
integer, parameter fct
Definition mumps_ooc_common.F:16
y
y
Definition schur_example.m:54
sol
sol
Definition schur_example.m:54
x
x
Definition schur_example.m:53