34
35
36
40 USE multi_fvm_mod
41
42
43
44#include "implicit_f.inc"
45
46
47
48#include "param_c.inc"
49#include "com04_c.inc"
50
51
52
53 TYPE (H3D_KEYWORD) H3D_KEYWORD_SHELL_SCALAR(*)
55 INTEGER IGEO(NPROPGI,NUMGEO)
56 TYPE(MULTI_FVM_STRUCT), INTENT(IN) :: MULTI_FVM
57 INTEGER, INTENT(INOUT) :: NKEY
58
59
60
61 INTEGER I,II,J,NIPMAX,IPLYMAX,CPTPLY,IPTMAX
62
63 i = 0
64
65 i = i + 1
66 h3d_keyword_shell_scalar(i)%KEY3 = 'MASS'
67 h3d_keyword_shell_scalar(i)%TEXT1 = 'Mass'
68
69 i = i + 1
70 h3d_keyword_shell_scalar(i)%KEY3 = 'DENS'
71 h3d_keyword_shell_scalar(i)%TEXT1 = 'Density'
72
73 i = i + 1
74 h3d_keyword_shell_scalar(i)%KEY3 = 'ENER'
75 h3d_keyword_shell_scalar(i)%TEXT1 = 'Internal Energy per unit mass'
76
77 i = i + 1
78 h3d_keyword_shell_scalar(i)%KEY3 = 'EINTM'
79 h3d_keyword_shell_scalar(i)%TEXT1 = 'Internal Energy per unit mass'
80
81 i = i + 1
82 h3d_keyword_shell_scalar(i)%KEY3 = 'EINTV'
83 h3d_keyword_shell_scalar(i)%TEXT1 = 'Internal Energy per unit volume'
84
85 i = i + 1
86 h3d_keyword_shell_scalar(i)%KEY3 = 'EINT'
87 h3d_keyword_shell_scalar(i)%TEXT1 = 'Internal Energy'
88
89 i = i + 1
90 h3d_keyword_shell_scalar(i)%KEY3 = 'TEMP'
91 h3d_keyword_shell_scalar(i)%TEXT1 = 'Temperature'
92
93 i = i + 1
94 h3d_keyword_shell_scalar(i)%KEY3 = 'THICK'
95 h3d_keyword_shell_scalar(i)%TEXT1 = 'Thickness'
96
97 i = i + 1
98 h3d_keyword_shell_scalar(i)%KEY3 = 'P'
99 h3d_keyword_shell_scalar(i)%TEXT1 = 'Pressure'
100
101 i = i + 1
102 h3d_keyword_shell_scalar(i)%KEY3 = 'VONM'
103 h3d_keyword_shell_scalar(i)%TEXT1 = 'Von Mises'
104
105 i = i + 1
106 h3d_keyword_shell_scalar(i)%KEY3 = 'K'
107 h3d_keyword_shell_scalar(i)%TEXT1 = 'Turbulent Energy'
108 h3d_keyword_shell_scalar(i)%COMMENT =
109 .'Specific for turbulent energy in CFD'
110
111 i = i + 1
112 h3d_keyword_shell_scalar(i)%KEY3 = 'TVIS'
113 h3d_keyword_shell_scalar(i)%TEXT1 = 'Turbulent Viscosity'
114 h3d_keyword_shell_scalar(i)%COMMENT = 'Specific for turbulent Viscosity in CFD'
115
116 i = i + 1
117 h3d_keyword_shell_scalar(i)%KEY3 = 'VORTX'
118 h3d_keyword_shell_scalar(i)%TEXT1 = 'Vorticity'
119
120 i = i + 1
121 h3d_keyword_shell_scalar(i)%KEY3 = 'DAM1'
122 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
123 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
124 h3d_keyword_shell_scalar(i)%IS_IPT = 1
125 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
126 h3d_keyword_shell_scalar(i)%IS_PLY = 1
127 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
128 h3d_keyword_shell_scalar(i)%TEXT1 = 'Damage 1'
129 h3d_keyword_shell_scalar(i)%COMMENT = 'Damage in material direction 1'
130
131 i = i + 1
132 h3d_keyword_shell_scalar(i)%KEY3 = 'DAM2'
133 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
134 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
135 h3d_keyword_shell_scalar(i)%IS_IPT = 1
136 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
137 h3d_keyword_shell_scalar(i)%IS_PLY = 1
138 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
139 h3d_keyword_shell_scalar(i)%TEXT1 = 'Damage 2'
140 h3d_keyword_shell_scalar(i)%COMMENT = 'Damage in material direction 2'
141
142 i = i + 1
143 h3d_keyword_shell_scalar(i)%KEY3 = 'DAM3'
144 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
145 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
146 h3d_keyword_shell_scalar(i)%IS_IPT = 1
147 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
148 h3d_keyword_shell_scalar(i)%IS_PLY = 1
149 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
150 h3d_keyword_shell_scalar(i)%TEXT1 = 'Damage 3'
151 h3d_keyword_shell_scalar(i)%COMMENT = 'Damage in material direction 3'
152
153 i = i + 1
154 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGX'
155 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress X'
156
157 i = i + 1
158 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGY'
159 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress Y'
160
161 i = i + 1
162 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGZ'
163 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress Z'
164
165 i = i + 1
166 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGXY'
167 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress XY'
168
169 i = i + 1
170 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGYZ'
171 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress YZ'
172
173 i = i + 1
174 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGZX'
175 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress ZX'
176
177 i = i + 1
178 h3d_keyword_shell_scalar(i)%KEY3 = 'USER'
179 h3d_keyword_shell_scalar(i)%IS_UVAR = 1
180 h3d_keyword_shell_scalar(i)%IS_UVAR_ALL = 1
181 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
182 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
183 h3d_keyword_shell_scalar(i)%IS_IPT = 1
184 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
185 h3d_keyword_shell_scalar(i)%IS_PLY = 1
186 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
187 h3d_keyword_shell_scalar(i)%TEXT1 = 'User Var'
188
189 i = i + 1
190 h3d_keyword_shell_scalar(i)%KEY3 = 'TSAIWU'
191 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
192 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
193 h3d_keyword_shell_scalar(i)%IS_IPT = 1
194 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
195 h3d_keyword_shell_scalar(i)%IS_PLY = 1
196 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
197 h3d_keyword_shell_scalar(i)%TEXT1 = 'Tsai-Wu Criterion'
198
199 i = i + 1
200 h3d_keyword_shell_scalar(i)%KEY3 = 'HOURGLASS'
201 h3d_keyword_shell_scalar(i)%TEXT1 = 'Hourglass Energy per unit mass'
202
203 IF (multi_fvm%IS_USED) THEN
204
205
206
207
208 DO ii = 1, multi_fvm%NBMAT
209
210 i = i + 1
211 WRITE(h3d_keyword_shell_scalar(i)%KEY3, '(A9,I0)') 'M151VFRAC', ii
212 WRITE(h3d_keyword_shell_scalar(i)%TEXT1, '(A17,I0)') 'Volume fraction ', ii
213 h3d_keyword_shell_scalar(i)%COMMENT = 'Volume fractions (for ALE multi-material Law151)'
214
215 i = i + 1
216 WRITE(h3d_keyword_shell_scalar(i)%KEY3, '(A8,I0)') 'M151DENS', ii
217 WRITE(h3d_keyword_shell_scalar(i)%TEXT1, '(A8,I0)') 'Density ', ii
218 h3d_keyword_shell_scalar(i)%COMMENT = 'Density (for ALE multi-material Law151)'
219
220 i = i + 1
221 WRITE(h3d_keyword_shell_scalar(i)%KEY3, '(A8,I0)') 'M151ENER', ii
222 WRITE(h3d_keyword_shell_scalar(i)%TEXT1, '(A7,I0)') 'Energy ', ii
223 h3d_keyword_shell_scalar(i)%COMMENT = 'Energy (for ALE multi-material Law151)'
224
225 i = i + 1
226 WRITE(h3d_keyword_shell_scalar(i)%KEY3, '(A8,I0)') 'M151PRES', ii
227 WRITE(h3d_keyword_shell_scalar(i)%TEXT1, '(A10,I0)') 'Pressure ', ii
228 h3d_keyword_shell_scalar(i)%COMMENT = 'Pressure (for ALE multi-material Law151)'
229
230 ENDDO
231
232 i = i + 1
233 h3d_keyword_shell_scalar(i)%KEY3 = 'ENTH'
234 h3d_keyword_shell_scalar(i)%TEXT1 = 'Enthalpy'
235
236 i = i + 1
237 h3d_keyword_shell_scalar(i)%KEY3 = 'ENTHM'
238 h3d_keyword_shell_scalar(i)%TEXT1 = 'Enthalpy per unit mass'
239
240 i = i + 1
241 h3d_keyword_shell_scalar(i)%KEY3 = 'ENTHV'
242 h3d_keyword_shell_scalar(i)%TEXT1 = 'Enthalpy per unit volume'
243
244 ENDIF
245
246 i = i + 1
247 h3d_keyword_shell_scalar(i)%KEY3 = 'EPSD'
248 h3d_keyword_shell_scalar(i)%TEXT1 = 'Strain Rate'
249 h3d_keyword_shell_scalar(i)%COMMENT = 'Equivalent strain rate (only available in case of strain rate filtering)'
250
251 i = i + 1
252 h3d_keyword_shell_scalar(i)%KEY3 = 'EPSP'
253 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
254 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
255 h3d_keyword_shell_scalar(i)%IS_IPT = 1
256 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
257 h3d_keyword_shell_scalar(i)%IS_PLY = 1
258 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
259 h3d_keyword_shell_scalar(i)%TEXT1 = 'Plastic Strain '
260
261 i = i + 1
262 h3d_keyword_shell_scalar(i)%KEY3 = 'NL_EPSD'
263 h3d_keyword_shell_scalar(i)%IS_IPT = 1
264 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
265 h3d_keyword_shell_scalar(i)%TEXT1 = 'Non-local Plastic Strain Rate'
266
267 i = i + 1
268 h3d_keyword_shell_scalar(i)%KEY3 = 'NL_EPSP'
269 h3d_keyword_shell_scalar(i)%IS_IPT = 1
270 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
271 h3d_keyword_shell_scalar(i)%TEXT1 = 'Non-local Plastic Strain '
272
273 i = i + 1
274 h3d_keyword_shell_scalar(i)%KEY3 = 'FAIL'
275 h3d_keyword_shell_scalar(i)%IS_PLY = 1
276 h3d_keyword_shell_scalar(i)%IS_IPT = 1
277 h3d_keyword_shell_scalar(i)%TEXT1 = 'Nb of Failed layers'
278 h3d_keyword_shell_scalar(i)%COMMENT =
279 .'/PROP/TYPE10, TYPE11, TYPE17, TYPE51, PCOMPP, /MAT/LAW15 & LAW25 (COMPSH)'
280
281 i = i + 1
282 h3d_keyword_shell_scalar(i)%KEY3 = 'CRMAS'
283 h3d_keyword_shell_scalar(i)%TEXT1 = 'Airbag crossing mass'
284
285 i = i + 1
286 h3d_keyword_shell_scalar(i)%KEY3 = 'CRVEL'
287 h3d_keyword_shell_scalar(i)%TEXT1 = 'Airbag crossing velocity'
288
289 i = i + 1
290 h3d_keyword_shell_scalar(i)%KEY3 = 'ALEBA/MASS'
291 h3d_keyword_shell_scalar(i)%TEXT1 = 'FVMBAG - Mass'
292
293 i = i + 1
294 h3d_keyword_shell_scalar(i)%KEY3 = 'ALEBA/PRES'
295 h3d_keyword_shell_scalar(i)%TEXT1 = 'FVMBAG - Pressure'
296
297 i = i + 1
298 h3d_keyword_shell_scalar(i)%KEY3 = 'ALEBA/MOMX'
299 h3d_keyword_shell_scalar(i)%TEXT1 = 'FVMBAG - Fluid velocity X'
300
301 i = i + 1
302 h3d_keyword_shell_scalar(i)%KEY3 = 'ALEBA/MOMY'
303 h3d_keyword_shell_scalar(i)%TEXT1 = 'FVMBAG - Fluid velocity Y'
304
305 i = i + 1
306 h3d_keyword_shell_scalar(i)%KEY3 = 'ALEBA/MOMZ'
307 h3d_keyword_shell_scalar(i)%TEXT1 = 'FVMBAG - Fluid velocity Z'
308
309 i = i + 1
310 h3d_keyword_shell_scalar(i)%KEY3 = 'ALEBA/DENS'
311 h3d_keyword_shell_scalar(i)%TEXT1 = 'FVMBAG - Density'
312
313 i = i + 1
314 h3d_keyword_shell_scalar(i)%KEY3 = 'ALEBA/ENER'
315 h3d_keyword_shell_scalar(i)%TEXT1 = 'FVMBAG - Energy'
316
317 i = i + 1
318 h3d_keyword_shell_scalar(i)%KEY3 = 'ALEBA/SSP'
319 h3d_keyword_shell_scalar(i)%TEXT1 = 'FVMBAG - Sound Speed'
320
321 i = i + 1
322 h3d_keyword_shell_scalar(i)%KEY3 = 'ALEBA/GAMA'
323 h3d_keyword_shell_scalar(i)%TEXT1 = 'FVMBAG - Gama'
324
325 i = i + 1
326 h3d_keyword_shell_scalar(i)%KEY3 = 'ALEBA/VISU'
327 h3d_keyword_shell_scalar(i)%TEXT1 = 'FVMBAG - Visu Polyhedra'
328
329 i = i + 1
330 h3d_keyword_shell_scalar(i)%KEY3 = 'THIN'
331 h3d_keyword_shell_scalar(i)%TEXT1 = '%Thinning'
332
333 i = i + 1
334 h3d_keyword_shell_scalar(i)%KEY3 = 'ERROR/THICK'
335 h3d_keyword_shell_scalar(i)%TEXT1 = 'Estimated Error on Thickness'
336
337 i = i + 1
338 h3d_keyword_shell_scalar(i)%KEY3 = 'PHI'
339 h3d_keyword_shell_scalar(i)%IS_PLY = 1
340 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
341 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
342 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
343 h3d_keyword_shell_scalar(i)%IS_IPT = 1
344 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
345 h3d_keyword_shell_scalar(i)%TEXT1 = 'PHI'
346 h3d_keyword_shell_scalar(i)%COMMENT =
347 .'angle between element skew and direction 1 orthotropy'
348
349 i = i + 1
350 h3d_keyword_shell_scalar(i)%KEY3 = 'DAMA'
351 h3d_keyword_shell_scalar(i)%IS_PLY = 1
352 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
353 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
354 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
355 h3d_keyword_shell_scalar(i)%IS_IPT = 1
356 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
357 h3d_keyword_shell_scalar(i)%TEXT1 = 'MAX DAMAGE'
358 h3d_keyword_shell_scalar(i)%COMMENT = 'maximum of damage over time and of all failure criteria acting in one material'
359
360 i = i + 1
361 h3d_keyword_shell_scalar(i)%KEY3 = 'FAILURE'
362 h3d_keyword_shell_scalar(i)%IS_PLY = 1
363 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
364 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
365 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
366 h3d_keyword_shell_scalar(i)%IS_IPT = 1
367 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
368 h3d_keyword_shell_scalar(i)%IS_ID = 1
369 h3d_keyword_shell_scalar(i)%IS_MODE = 1
370 h3d_keyword_shell_scalar(i)%TEXT1 = 'Failure'
371 h3d_keyword_shell_scalar(i)%COMMENT = 'Damage of current failure criterion'
372
373 i = i + 1
374 h3d_keyword_shell_scalar(i)%KEY3 = 'DAMINI'
375 h3d_keyword_shell_scalar(i)%IS_PLY = 1
376 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
377 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
378 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
379 h3d_keyword_shell_scalar(i)%IS_IPT = 1
380 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
381 h3d_keyword_shell_scalar(i)%TEXT1 = 'Damage initiation'
382 h3d_keyword_shell_scalar(i)%COMMENT = 'Damage initiation variable'
383
384 i = i + 1
385 h3d_keyword_shell_scalar(i)%KEY3 = 'DAMA/MEMB'
386 h3d_keyword_shell_scalar(i)%TEXT1 = 'MAX DAMAGE MEMBRANE'
387
388 i = i + 1
389 h3d_keyword_shell_scalar(i)%KEY3 = 'DAMG'
390 h3d_keyword_shell_scalar(i)%IS_IPT = 1
391 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
392 h3d_keyword_shell_scalar(i)%IS_PLY = 1
393 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
394 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
395 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
396 h3d_keyword_shell_scalar(i)%IS_ID = 1
397 h3d_keyword_shell_scalar(i)%IS_MODE = 1
398 h3d_keyword_shell_scalar(i)%IS_MODE_ALL = 1
399 h3d_keyword_shell_scalar(i)%TEXT1 = 'Damage'
400 h3d_keyword_shell_scalar(i)%COMMENT = 'Damage variable from material law'
401
402 i = i + 1
403 h3d_keyword_shell_scalar(i)%KEY3 = 'DAMG/MEMB'
404 h3d_keyword_shell_scalar(i)%TEXT1 = 'Damage (Membrane value)'
405
406 i = i + 1
407 h3d_keyword_shell_scalar(i'NXTF'
408 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
409 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
410 h3d_keyword_shell_scalar(i)%IS_IPT = 1
411 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
412 h3d_keyword_shell_scalar(i)%TEXT1 = 'NXT FAILURE FACTOR'
413
414 i = i + 1
415 h3d_keyword_shell_scalar(i)%KEY3 = 'NXTF/MEMB'
416 h3d_keyword_shell_scalar(i)%TEXT1 = 'NXT FAILURE FACTOR MEMBRANE'
417
418 i = i + 1
419 h3d_keyword_shell_scalar(i)%KEY3 = 'TDEL'
420 h3d_keyword_shell_scalar(i)%TEXT1 = 'TIME DELETION ELEMENT'
421 h3d_keyword_shell_scalar(i)%COMMENT = 'output only for elements which are deleted, due to a /FAIL criterion.'
422
423 i = i + 1
424 h3d_keyword_shell_scalar(i)%KEY3 = 'SSP'
425 h3d_keyword_shell_scalar(i)%TEXT1 = 'Sound Speed'
426 h3d_keyword_shell_scalar(i)%COMMENT = 'Sound speed. Only available with ALE material laws'
427
428 i = i + 1
429 h3d_keyword_shell_scalar(i)%KEY3 = 'SCHLIEREN'
430 h3d_keyword_shell_scalar(i)%TEXT1 = 'Schlieren'
431 h3d_keyword_shell_scalar(i)%COMMENT = 'Schlieren image (optical method widely used in CFD field). ALE material laws'
432
433 i = i + 1
434 h3d_keyword_shell_scalar(i)%KEY3 = 'PHI/MEMB'
435 h3d_keyword_shell_scalar(i)%TEXT1 = 'phi MEMBRANE'
436
437 i = i + 1
438 h3d_keyword_shell_scalar(i)%KEY3 = 'DOMAIN'
439 h3d_keyword_shell_scalar(i)%TEXT1 = 'domain'
440
441 I = I + 1
442 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'sigeq'
443 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'equiv stress'
444 H3D_KEYWORD_SHELL_SCALAR(I)%COMMENT = 'equivalent stress, based on yield criteria used
for the corresponding material
'
445
446 I = I + 1
447 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'bulk'
448 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'artificial viscosity'
449
450 I = I + 1
451 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'tdet'
452 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'detonation time'
453
454 I = I + 1
455 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'bfrac'
456 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'burn fraction'
457
458 I = I + 1
459 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'dt'
460 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'time step'
461
462 I = I + 1
463 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'ams'
464 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'ams selection'
465
466 I = I + 1
467 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'alpha'
468 H3D_KEYWORD_SHELL_SCALAR(I)%IS_PLY = 1
469 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER = 1
470 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'alpha'
471 H3D_KEYWORD_SHELL_SCALAR(I)%COMMENT = 'shear angle'
472
473 I = I + 1
474 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'off'
475 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'element status'
476 H3D_KEYWORD_SHELL_SCALAR(I)%COMMENT = '0.0:deleted >0.&<1.:under failure process 1.0:activated <0:standby'
477
478!to be moved in h3d_tria_scalar when implemented.
479 I = I + 1
480 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'mach'
481 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'mach number'
482 H3D_KEYWORD_SHELL_SCALAR(I)%COMMENT = 'mach number. only available with ale material law 151'
483
484!to be moved in h3d_tria_scalar when implemented.
485 I = I + 1
486 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'color'
487 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'color function'
488 H3D_KEYWORD_SHELL_SCALAR(I)%COMMENT = 'color
function for multimaterial and multifluid ale laws
'
489
490 I = I + 1
491 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'fldf/memb'
492 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fld failure factor membrane'
493
494 I = I + 1
495 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'fldf'
496 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fld failure factor'
497 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT = 1
498 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT_ALL = 1
499 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER = 1
500 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER_ALL = 1
501
502 I = I + 1
503 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'fldz/memb'
504 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fld zone index membrane'
505 H3D_KEYWORD_SHELL_SCALAR(I)%COMMENT = '1:highwrinkle 2:compression 3:loosemetal 4:safezone 5:margintofail 6:failurezone'
506
507 I = I + 1
508 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'fldz'
509 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fld zone index'
510 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT = 1
511 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT_ALL = 1
512 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER = 1
513 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER_ALL = 1
514 H3D_KEYWORD_SHELL_SCALAR(I)%COMMENT = '1:highwrinkle 2:compression 3:loosemetal 4:safezone 5:margintofail 6:failurezone'
515
516 I = I + 1
517 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'hc_dsse_f/memb'
518 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'hc_dsse failure factor membrane'
519
520 I = I + 1
521 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'hc_dsse_f'
522 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'hc_dsse failure factor'
523 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT = 1
524 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT_ALL = 1
525 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER = 1
526 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER_ALL = 1
527 H3D_KEYWORD_SHELL_SCALAR(I)%IS_PLY = 1
528 H3D_KEYWORD_SHELL_SCALAR(I)%IS_PLY_ALL = 1
529
530 I = I + 1
531 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'hc_dsse_z/memb'
532 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'hc_dsse zone index membrane'
533 H3D_KEYWORD_SHELL_SCALAR(I)%COMMENT = '1:below hc and dsse 2: above dsse'
534
535 I = I + 1
536 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'hc_dsse_z'
537 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'hc_dsse zone index'
538 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT = 1
539 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT_ALL = 1
540 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER = 1
541 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER_ALL = 1
542 H3D_KEYWORD_SHELL_SCALAR(I)%IS_PLY = 1
543 H3D_KEYWORD_SHELL_SCALAR(I)%IS_PLY_ALL = 1
544 H3D_KEYWORD_SHELL_SCALAR(I)%COMMENT = '1:below hc and dsse 2: above dsse and below hc 3: above hc and dsse'
545
546 I = I + 1
547 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'wpla'
548 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER = 1
549 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER_ALL = 1
550 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT = 1
551 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT_ALL = 1
552 H3D_KEYWORD_SHELL_SCALAR(I)%IS_PLY = 1
553 H3D_KEYWORD_SHELL_SCALAR(I)%IS_PLY_ALL = 1
554 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'plastic work'
555
556 I = I + 1
557 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'mds'
558 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'mds'
559 H3D_KEYWORD_SHELL_SCALAR(I)%IS_MDSVAR = 1
560 H3D_KEYWORD_SHELL_SCALAR(I)%IS_MDSVAR_ALL = 1
561 H3D_KEYWORD_SHELL_SCALAR(I)%IS_MDSVAR_DEF = 1
562 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER = 1
563 H3D_KEYWORD_SHELL_SCALAR(I)%IS_LAYER_ALL = 1
564 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT = 1
565 H3D_KEYWORD_SHELL_SCALAR(I)%IS_IPT_ALL = 1
566 H3D_KEYWORD_SHELL_SCALAR(I)%IS_PLY = 1
567 H3D_KEYWORD_SHELL_SCALAR(I)%IS_PLY_ALL = 1
568
569 I = I + 1
570 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'pext'
571 h3d_keyword_shell_scalar(i)%IS_SKIN = 1
572 h3d_keyword_shell_scalar(i)%TEXT1 = 'Pressure OUTER'
573 h3d_keyword_shell_scalar(i)%COMMENT = 'External pressure on the Skin of shell'
574
575 i = i + 1
576 h3d_keyword_shell_scalar(i)%KEY3 = 'GROUP'
577 h3d_keyword_shell_scalar(i)%TEXT1 = 'Group identifier'
578
579 i = i + 1
580 h3d_keyword_shell_scalar(i)%KEY3 = 'INTERNAL.ID'
581 h3d_keyword_shell_scalar(i)%TEXT1 = 'Internal identifier'
582
583 i = i + 1
584 h3d_keyword_shell_scalar(i)%KEY3 = 'LOCAL.ID'
585 h3d_keyword_shell_scalar(i)%TEXT1 = 'Local identifier within the group'
586
587 i
588 h3d_keyword_shell_scalar(i)%KEY3 = 'ENER/TMAX'
589 h3d_keyword_shell_scalar(i)%TEXT1 = 'Max of Internal Energy per unit mass'
590 h3d_keyword_shell_scalar(i)%COMMENT = 'Max of Internal Energy per unit mass over time'
591
592 i = i + 1
593 h3d_keyword_shell_scalar(i)%KEY3 = 'VONM/TMAX'
594 h3d_keyword_shell_scalar(i)%TEXT1 = 'Max of Von Mises'
595 h3d_keyword_shell_scalar(i)%COMMENT = 'Max of Von Mises stress over time&ipt'
596
597 i = i + 1
598 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGEQ/TMAX'
599 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
600 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
601 h3d_keyword_shell_scalar(i)%IS_IPT = 1
602 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
603 h3d_keyword_shell_scalar(i)%IS_PLY
604 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
605 h3d_keyword_shell_scalar(i)%TEXT1 = 'Max of Equiv stress'
606 h3d_keyword_shell_scalar(i)%COMMENT = 'Max of Equivalent stress over time&ipt'
607
608
609 h3d_keyword_shell_scalar(i)%KEY3 = 'DAMA/TMAX'
610
611 h3d_keyword_shell_scalar(i)%IS_LAYER_ALL = 1
612 h3d_keyword_shell_scalar(i)%IS_IPT = 1
613 h3d_keyword_shell_scalar(i)%IS_IPT_ALL = 1
614 h3d_keyword_shell_scalar(i)%IS_PLY = 1
615 h3d_keyword_shell_scalar(i)%IS_PLY_ALL = 1
616 h3d_keyword_shell_scalar(i)%TEXT1 = 'MAX of DAMAGE'
617 h3d_keyword_shell_scalar(i)%COMMENT = 'maximum of damage over time&ipt'
618
619
620 i = i + 1
621 h3d_keyword_shell_scalar(i)%KEY3 = 'DIV(U)'
622 h3d_keyword_shell_scalar(i)%TEXT1 = 'div(u)'
623 h3d_keyword_shell_scalar(i)%COMMENT = 'Volumetric dilatation rate'
624
625
626 i = i + 1
627 h3d_keyword_shell_scalar(i)%KEY3 = 'VSTRAIN'
628 h3d_keyword_shell_scalar(i)%TEXT1 = 'Volumetric Strain'
629 h3d_keyword_shell_scalar(i)%COMMENT = 'mu=rho/rho0-1'
630
631 i = i + 1
632 h3d_keyword_shell_scalar(i)%KEY3 = 'VSTRAIN/1'
633 h3d_keyword_shell_scalar(i)%TEXT1 = 'Volumetric Strain 1'
634 h3d_keyword_shell_scalar(i)%COMMENT = 'mu=rho/rho0-1'
635
636 i = i + 1
637 h3d_keyword_shell_scalar(i)%KEY3 = 'VSTRAIN/2'
638 h3d_keyword_shell_scalar(i)%TEXT1 = 'Volumetric Strain 2'
639 h3d_keyword_shell_scalar(i)%COMMENT = 'mu=rho/rho0-1'
640
641
642 h3d_keyword_shell_scalar(i)%KEY3 = 'VSTRAIN/3'
643 h3d_keyword_shell_scalar(i)%TEXT1 = 'Volumetric Strain 3'
644 h3d_keyword_shell_scalar(i)%COMMENT = 'mu=rho/rho0-1'
645
646 i = i + 1
647 h3d_keyword_shell_scalar(i)%KEY3 = 'VSTRAIN/4'
648 h3d_keyword_shell_scalar(i)%TEXT1 = 'Volumetric Strain 4'
649 h3d_keyword_shell_scalar(i)%COMMENT = 'mu=rho/rho0-1'
650
651 i = i + 1
652 h3d_keyword_shell_scalar(i)%KEY3 = 'VSTRAIN/5'
653 h3d_keyword_shell_scalar(i)%TEXT1 = 'Volumetric Strain 5'
654 h3d_keyword_shell_scalar(i)%COMMENT = 'mu=rho/rho0-1'
655
656 i = i + 1
657 h3d_keyword_shell_scalar(i)%KEY3 = 'VSTRAIN/6'
658 h3d_keyword_shell_scalar(i)%TEXT1 = 'Volumetric Strain 6'
659 h3d_keyword_shell_scalar(i)%COMMENT = 'mu=rho/rho0-1'
660
661 i = i + 1
662 h3d_keyword_shell_scalar(i)%KEY3 = 'VSTRAIN/7'
663 h3d_keyword_shell_scalar(i)%TEXT1 = 'Volumetric Strain 7'
664 h3d_keyword_shell_scalar(i)%COMMENT = 'mu=rho/rho0-1'
665
666 i = i + 1
667 h3d_keyword_shell_scalar(i)%KEY3 = 'VSTRAIN/8'
668 h3d_keyword_shell_scalar(i)%TEXT1 = 'Volumetric Strain 8'
669 h3d_keyword_shell_scalar(i)%COMMENT = 'mu=rho/rho0-1'
670
671 i = i + 1
672 h3d_keyword_shell_scalar(i)%KEY3 = 'VSTRAIN/9'
673 h3d_keyword_shell_scalar(i)%TEXT1 = 'Volumetric Strain 9'
674 h3d_keyword_shell_scalar(i)%COMMENT = 'mu=rho/rho0-1'
675
676 i = i + 1
677 h3d_keyword_shell_scalar(i)%KEY3 = 'VSTRAIN/10'
678 h3d_keyword_shell_scalar(i)%TEXT1 = 'Volumetric Strain 10'
679 h3d_keyword_shell_scalar(i)%COMMENT = 'mu=rho/rho0-1'
680
681
682 nkey=i
end diagonal values have been computed in the(sparse) matrix id.SOL
for(i8=*sizetab-1;i8 >=0;i8--)