33
34
35
38 USE multi_fvm_mod
39
40
41
42#include "implicit_f.inc"
43
44
45
46 TYPE (H3D_KEYWORD) H3D_KEYWORD_SOLID_SCALAR(*)
48 INTEGER IGEO(NPROPGI,NUMGEO)
49 TYPE(MULTI_FVM_STRUCT), INTENT(IN) :: MULTI_FVM
50 INTEGER, INTENT(INOUT) :: NKEY
51
52
53
54#include "com04_c.inc"
55#include "param_c.inc"
56
57
58
59 INTEGER I,II,J,NIPMAX
60
61
62
63 i = 0
64
65 i = i + 1
66 h3d_keyword_solid_scalar(i)%KEY3 = 'MASS'
67 h3d_keyword_solid_scalar(i)%TEXT1 = 'Mass'
68
69 i = i + 1
70 h3d_keyword_solid_scalar(i)%KEY3 = 'EPSP'
71 h3d_keyword_solid_scalar(i)%IS_LAYER = 1
72 h3d_keyword_solid_scalar(i)%IS_LAYER_ALL = 1
73 h3d_keyword_solid_scalar(i)%TEXT1 = 'Plastic Strain'
74
75 i = i + 1
76 h3d_keyword_solid_scalar(i)%KEY3 = 'NL_EPSP'
77 h3d_keyword_solid_scalar(i)%IS_LAYER = 1
78 h3d_keyword_solid_scalar(i)%IS_LAYER_ALL = 1
79 h3d_keyword_solid_scalar(i)%TEXT1 = 'Non-local Plastic Strain'
80
81 i = i + 1
82 h3d_keyword_solid_scalar(i)%KEY3 = 'DENS'
83 h3d_keyword_solid_scalar(i)%TEXT1 = 'Density'
84
85 i = i + 1
86 h3d_keyword_solid_scalar(i)%KEY3 = 'DT'
87 h3d_keyword_solid_scalar(i)%TEXT1 = 'Time Step'
88
89 i = i + 1
90 h3d_keyword_solid_scalar(i)%KEY3 = 'ENER'
91 h3d_keyword_solid_scalar(i)%TEXT1 = 'Internal Energy per unit mass'
92
93 i = i + 1
94 h3d_keyword_solid_scalar(i)%KEY3 = 'TEMP'
95 h3d_keyword_solid_scalar(i)%TEXT1 = 'temperature'
96
97 I = I + 1
98 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'p'
99 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'pressure'
100
101 I = I + 1
102 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vonm'
103 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'von mises'
104
105 I = I + 1
106 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'k'
107 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'turbulent energy'
108 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'specific
for turbulent energy in cfd
'
109
110 I = I + 1
111 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'tvis'
112 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'turbulent viscosity'
113 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'specific
for turbulent viscosity in cfd
'
114
115 I = I + 1
116 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vortx'
117 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'vorticity-x'
118
119 I = I + 1
120 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vorty'
121 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'vorticity-y'
122
123 I = I + 1
124 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vortz'
125 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'vorticity-z'
126
127 I = I + 1
128 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vort'
129 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'vorticity'
130
131 I = I + 1
132 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'dam1'
133 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'damage 1'
134
135 I = I + 1
136 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'dam2'
137 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'damage 2'
138
139 I = I + 1
140 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'dam3'
141 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'damage 3'
142
143 I = I + 1
144 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'sigx'
145 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'stress x'
146
147 I = I + 1
148 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'sigy'
149 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'stress y'
150
151 I = I + 1
152 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'sigz'
153 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'stress z'
154
155 I = I + 1
156 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'sigxy'
157 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'stress'
158
159 I = I + 1
160 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'sigyz'
161 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'stress yz'
162
163 I = I + 1
164 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'sigzx'
165 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'stress zx'
166
167 I = I + 1
168 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'user'
169 H3D_KEYWORD_SOLID_SCALAR(I)%IS_UVAR = 1
170 H3D_KEYWORD_SOLID_SCALAR(I)%IS_UVAR_ALL = 1
171 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'user var'
172
173 I = I + 1
174 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'hourglass'
175 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'hourglass energy per unit mass'
176
177 I = I + 1
178 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'epsd'
179 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'strain rate'
180 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'equivalent strain rate(only available in case of strain rate filtering)'
181
182 I = I + 1
183 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'nl_epsd'
184 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'non-local plastic strain rate'
185
186 I = I + 1
187 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'flay'
188 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'failed layers'
189 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'number of failed layers'
190
191 IF (MULTI_FVM%IS_USED) THEN
192
193
194
195
196 DO II = 1, MULTI_FVM%NBMAT
197
198 I = I + 1
199 WRITE(H3D_KEYWORD_SOLID_SCALAR(I)%KEY3, '(a9,i0)') 'm151vfrac', II
200 WRITE(H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1, '(a17,i0)') 'volume fraction', II
201 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'volume fractions(
for ale multi-material law151)
'
202
203 I = I + 1
204 WRITE(H3D_KEYWORD_SOLID_SCALAR(I)%KEY3, '(a8,i0)') 'm151dens', II
205 WRITE(H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1, '(a8,i0)') 'density ', II
206 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'density(
for ale multi-material law151)
'
207
208 I = I + 1
209 WRITE(H3D_KEYWORD_SOLID_SCALAR(I)%KEY3, '(a8,i0)') 'm151ener', II
210 WRITE(H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1, '(a7,i0)') 'energy ', II
211 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'energy(
for ale multi-material law151)
'
212
213 I = I + 1
214 WRITE(H3D_KEYWORD_SOLID_SCALAR(I)%KEY3, '(a8,i0)') 'm151pres', II
215 WRITE(H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1, '(a10,i0)') 'pressure ', II
216 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'pressure(
for ale multi-material law151)
'
217
218 ENDDO
219 ELSE
220
221 I = I + 1
222 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vfrac1'
223 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'volumetric fraction 1'
224 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'volumetric fractions(
for ale multi-material laws: law20, law37 and law51)
'
225
226 I = I + 1
227 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vfrac2'
228 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'volumetric fraction 2'
229 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'volumetric fractions(
for ale multi-material laws: law20, law37 and law51)
'
230
231 I = I + 1
232 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vfrac3'
233 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'volumetric fraction 3'
234 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'volumetric fractions(
for ale multi-material laws: law20, law37 and law51)
'
235
236 I = I + 1
237 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vfrac4'
238 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'volumetric fraction 4'
239 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'volumetric fractions (
for ale multi-material laws: law20, law37 and law51)
'
240
241 ENDIF ! MULTI_FVM%IS_USED
242
243 I = I + 1
244 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'orthd/psi'
245 H3D_KEYWORD_SOLID_SCALAR(I)%IS_LAYER = 1
246 H3D_KEYWORD_SOLID_SCALAR(I)%IS_LAYER_ALL = 1
247 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR = 1
248 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR_ALL = 1
249 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS = 1
250 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS_ALL = 1
251 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT = 1
252 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT_ALL = 1
253 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'psi'
254 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'euler angles psi, global reference system -> orthotropic reference system'
255
256 I = I + 1
257 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'orthd/theta'
258 H3D_KEYWORD_SOLID_SCALAR(I)%IS_LAYER = 1
259 H3D_KEYWORD_SOLID_SCALAR(I)%IS_LAYER_ALL = 1
260 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR = 1
261 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR_ALL = 1
262 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS = 1
263 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS_ALL = 1
264 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT = 1
265 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT_ALL = 1
266 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'teta'
267 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'euler angles theta, global reference system -> orthotropic reference system'
268
269 I = I + 1
270 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'orthd/phi'
271 H3D_KEYWORD_SOLID_SCALAR(I)%IS_LAYER = 1
272 H3D_KEYWORD_SOLID_SCALAR(I)%IS_LAYER_ALL = 1
273 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR = 1
274 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR_ALL = 1
275 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS = 1
276 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS_ALL = 1
277 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT = 1
278 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT_ALL = 1
279 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'phi'
280 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'euler angles phi, global reference system -> orthotropic reference system'
281
282
283
284
285
286 I = I + 1
287 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'bfrac'
288 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'burn fraction'
289 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'burn fraction(
for jwl eos: law5 and law51)
'
290
291 I = I + 1
292 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vdam1'
293 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR = 1
294 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR_ALL = 1
295 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS = 1
296 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS_ALL = 1
297 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT = 1
298 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT_ALL = 1
299 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'damage variable1'
300 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'Value of damage factor d'
301
302 I = I + 1
303 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vdam2'
304 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR = 1
305 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR_ALL = 1
306 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS = 1
307 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS_ALL = 1
308 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT = 1
309 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT_ALL = 1
310 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'damage variable2'
311 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'Value of damage surface sigmay0'
312
313 I = I + 1
314 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'vdam3'
315 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR = 1
316 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR_ALL = 1
317 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS = 1
318 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS_ALL = 1
319 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT = 1
320 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT_ALL = 1
321 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'damage variable3'
322 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'Value of rupture surface sigmayf'
323
324 I = I + 1
325 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'dama'
326 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR = 1
327 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR_ALL = 1
328 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS = 1
329 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS_ALL = 1
330 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT = 1
331 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT_ALL = 1
332 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'max damage
'
333 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'maximum of damage over time and of all failure criteria acting in one material'
334
335 I = I + 1
336 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'damg'
337 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR = 1
338 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR_ALL = 1
339 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS = 1
340 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS_ALL = 1
341 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT = 1
342 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT_ALL = 1
343 H3D_KEYWORD_SOLID_SCALAR(I)%IS_LAYER = 1
344 H3D_KEYWORD_SOLID_SCALAR(I)%IS_LAYER_ALL = 1
345 H3D_KEYWORD_SOLID_SCALAR(I)%IS_ID = 1
346 H3D_KEYWORD_SOLID_SCALAR(I)%IS_MODE = 1
347 H3D_KEYWORD_SOLID_SCALAR(I)%IS_MODE_ALL = 1
348 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'damage'
349 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'damage variable from material law'
350
351 I = I + 1
352 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'damini'
353 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR = 1
354 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR_ALL = 1
355 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS = 1
356 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS_ALL = 1
357 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT = 1
358 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT_ALL = 1
359 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'damage initiation'
360 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'damage initiation'
361
362 I = I + 1
363 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'failure'
364 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR = 1
365 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IR_ALL = 1
366 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS = 1
367 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IS_ALL = 1
368 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT = 1
369 H3D_KEYWORD_SOLID_SCALAR(I)%IS_IT_ALL = 1
370 H3D_KEYWORD_SOLID_SCALAR(I)%IS_ID = 1
371 H3D_KEYWORD_SOLID_SCALAR(I)%IS_MODE = 1
372 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'failure'
373 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'damage of current failure criterion'
374
375 I = I + 1
376 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'tsaiwu'
377 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'tsai-wu criterion'
378
379 I = I + 1
380 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'tdel'
381 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'time deletion element'
382 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'output only
for elements which are deleted, due to a /fail criterion.
'
383
384 I = I + 1
385 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'ssp'
386 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'sound speed'
387 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'sound speed. only available with ale material laws'
388
389 I = I + 1
390 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'volu'
391 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'volume'
392 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'element volume'
393
394 I = I + 1
395 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'schlieren'
396 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'schlieren'
397 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'schlieren image(optical method widely used in cfd field). ale material laws'
398
399 I = I + 1
400 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'domain'
401 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'domain'
402
403 I = I + 1
404 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'fill'
405 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'filling percentage'
406
407 I = I + 1
408 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'sigeq'
409 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'equiv stress'
410 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'equivalent stress, based on yield criteria used
for the corresponding material
'
411
412 I = I + 1
413 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'bulk'
414 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'artificial viscosity'
415
416 I = I + 1
417 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'tdet'
418 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'detonation time'
419
420 I = I + 1
421 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'momx'
422 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'momentum density x'
423
424 I = I + 1
425 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'momy'
426 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'momentum density y'
427
428 I = I + 1
429 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'momz'
430 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'momentum density z'
431
432 i = i + 1
433 h3d_keyword_solid_scalar(i)%KEY3 = 'MOMXY'
434 h3d_keyword_solid_scalar(i)%TEXT1 = 'Momentum Density XY'
435
436 i = i + 1
437 h3d_keyword_solid_scalar(i)%KEY3 = 'MOMYZ'
438 h3d_keyword_solid_scalar(i)%TEXT1 = 'Momentum Density YZ'
439
440 i = i + 1
441 h3d_keyword_solid_scalar(i)%KEY3 = 'MOMXZ'
442 h3d_keyword_solid_scalar(i)%TEXT1 = 'Momentum Density XZ'
443
444 i = i + 1
445 h3d_keyword_solid_scalar(i)%KEY3 = '|MOM|'
446 h3d_keyword_solid_scalar(i)%TEXT1 = 'Momentum Density ABS'
447
448 i = i + 1
449 h3d_keyword_solid_scalar(i)%KEY3 = 'VELX'
450 h3d_keyword_solid_scalar(i)%TEXT1 = 'Velocity X'
451
452 i = i + 1
453 h3d_keyword_solid_scalar(i)%KEY3 = 'VELY'
454 h3d_keyword_solid_scalar(i)%TEXT1 = 'Velocity Y'
455
456 i = i + 1
457 h3d_keyword_solid_scalar(i)%KEY3 = 'VELZ'
458 h3d_keyword_solid_scalar(i)%TEXT1 = 'Velocity Z'
459
460 i = i + 1
461 h3d_keyword_solid_scalar(i)%KEY3 = 'VELXY'
462 h3d_keyword_solid_scalar(i)%TEXT1 = 'Velocity XY'
463
464 i = i + 1
465 h3d_keyword_solid_scalar(i)%KEY3 = 'VELYZ'
466 h3d_keyword_solid_scalar(i)%TEXT1 = 'Velocity YZ'
467
468 i = i + 1
469 h3d_keyword_solid_scalar(i)%KEY3 = 'VELXZ'
470 h3d_keyword_solid_scalar(i)%TEXT1 = 'Velocity XZ'
471
472 i = i + 1
473 h3d_keyword_solid_scalar(i)%KEY3 = '|VEL|'
474 h3d_keyword_solid_scalar(i)%TEXT1 = 'Velocity ABS'
475
476 i = i + 1
477 h3d_keyword_solid_scalar(i)%KEY3 = 'ENTH'
478 h3d_keyword_solid_scalar(i)%TEXT1 = 'Enthalpy'
479
480 i = i + 1
481 h3d_keyword_solid_scalar(i)%KEY3 = 'ENTHM'
482 h3d_keyword_solid_scalar(i)%TEXT1 = 'Enthalpy per unit mass'
483
484 i = i + 1
485 h3d_keyword_solid_scalar(i)%KEY3 = 'ENTHV'
486 h3d_keyword_solid_scalar(i)%TEXT1 = 'Enthalpy per unit volume'
487
488 i = i + 1
489 h3d_keyword_solid_scalar(i)%KEY3 = 'EINT'
490 h3d_keyword_solid_scalar(i)%TEXT1 = 'Internal Energy'
491
492 i = i + 1
493 h3d_keyword_solid_scalar(i)%KEY3 = 'EINTM'
494 h3d_keyword_solid_scalar(i)%TEXT1 = 'Internal Energy per unit mass'
495
496 i = i + 1
497 h3d_keyword_solid_scalar(i)%KEY3 = 'EINTV'
498 h3d_keyword_solid_scalar(i)%TEXT1 = 'Internal Energy per unit volume'
499
500 i = i + 1
501 h3d_keyword_solid_scalar(i)%KEY3 = 'AMS'
502 h3d_keyword_solid_scalar(i)%TEXT1 = 'AMS selection'
503
504 i = i + 1
505 h3d_keyword_solid_scalar(i)%KEY3 = 'OFF'
506 h3d_keyword_solid_scalar(i)%TEXT1 = 'Element status'
507 h3d_keyword_solid_scalar(i)%COMMENT = '0.0:deleted >0.&<1.:under failure process 1.0:activated <0:Standby'
508
509 i = i + 1
510 h3d_keyword_solid_scalar(i)%KEY3 = 'MACH'
511 h3d_keyword_solid_scalar(i)%TEXT1 = 'Mach Number'
512 h3d_keyword_solid_scalar(i)%COMMENT = 'Mach Number. Only available with ALE material law 151'
513
514 i = i + 1
515 h3d_keyword_solid_scalar(i)%KEY3 = 'COLOR'
516 h3d_keyword_solid_scalar(i)%TEXT1 = 'Color Function'
517 h3d_keyword_solid_scalar(i)%COMMENT = 'Color function for multimaterial and multifluid ALE laws'
518
519 i = i + 1
520 h3d_keyword_solid_scalar(i)%KEY3 = 'WPLA'
521 h3d_keyword_solid_scalar(i)%TEXT1 = 'Plastic Work'
522
523 i = i + 1
524 h3d_keyword_solid_scalar(i)%KEY3 = 'FLDF/OUTER'
525 h3d_keyword_solid_scalar(i)%IS_SKIN = 1
526 h3d_keyword_solid_scalar(i)%TEXT1 = 'FLD Failure Factor OUTER'
527 h3d_keyword_solid_scalar(i)%COMMENT = 'FLD Failure Factor on the external surface of solid'
528
529 i = i + 1
530 h3d_keyword_solid_scalar(i)%KEY3 = 'FLDZ/OUTER'
531 h3d_keyword_solid_scalar(i)%IS_SKIN = 1
532 h3d_keyword_solid_scalar(i)%TEXT1 = 'FLD Zone Index OUTER'
533 h3d_keyword_solid_scalar(i)%COMMENT = '1:HighWrinkle 2:Compression 3:LooseMetal 4:SafeZone 5:MarginToFail 6:FailureZone'
534
535 i = i + 1
536 h3d_keyword_solid_scalar(i)%KEY3 = 'MDS'
537 h3d_keyword_solid_scalar(i)%TEXT1 = 'MDS'
538 h3d_keyword_solid_scalar(i)%IS_MDSVAR = 1
539 h3d_keyword_solid_scalar(i)%IS_MDSVAR_ALL = 1
540 h3d_keyword_solid_scalar(i)%IS_MDSVAR_DEF = 1
541
542 i = i + 1
543 h3d_keyword_solid_scalar(i)%KEY3 = 'PEXT'
544 h3d_keyword_solid_scalar(i)%IS_SKIN = 1
545 h3d_keyword_solid_scalar(i)%TEXT1 = 'Pressure OUTER'
546 h3d_keyword_solid_scalar(i)%COMMENT = 'External pressure on the Skin of solid'
547
548 i = i + 1
549 h3d_keyword_solid_scalar(i)%KEY3 = 'GROUP'
550 h3d_keyword_solid_scalar(i)%TEXT1 = 'Group identifier'
551
552 i = i + 1
553 h3d_keyword_solid_scalar(i)%KEY3 = 'INTERNAL.ID'
554 h3d_keyword_solid_scalar(i)%TEXT1 = 'Internal identifier'
555
556 i = i + 1
557 h3d_keyword_solid_scalar(i)%KEY3 = 'LOCAL.ID'
558 h3d_keyword_solid_scalar(i)%TEXT1 = 'Local identifier within the group'
559
560 i = i + 1
561 h3d_keyword_solid_scalar(i)%KEY3 = 'THICK'
562 h3d_keyword_solid_scalar(i)%TEXT1 = 'Thickness : direction s'
563
564 i = i + 1
565 h3d_keyword_solid_scalar(i)%KEY3 = 'THIN'
566 h3d_keyword_solid_scalar(i)%TEXT1 = 'Thinning(percentage) : direction s'
567
568 i = i + 1
569 h3d_keyword_solid_scalar(i)%KEY3 = 'ENER/TMAX'
570 h3d_keyword_solid_scalar(i)%TEXT1 = 'Max of Internal Energy per unit mass'
571 h3d_keyword_solid_scalar(i)%COMMENT = 'Max of Internal Energy per unit mass over time'
572
573 i = i + 1
574 h3d_keyword_solid_scalar(i)%KEY3 = 'vonm/tmax'
575 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'max of von mises
'
576 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'max of von mises stress over time
'
577
578 I = I + 1
579 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'sigeq/tmax'
580 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'max of equiv stress
'
581 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'max of equivalent stress over time&ipt
'
582
583 I = I + 1
584 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'dama/tmax'
585 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'max of damage
'
586 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'maximum of damage over time&ipt'
587
588 I = I + 1
589 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'tillotson'
590 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'region identifier in p,v diagram
'
591 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = '1:compressed, 2:expanded cold, 3:transition, 4:hot expansion'
592
593 I = I + 1
594 H3D_KEYWORD_SOLID_SCALAR(I)%KEY3 = 'fldf/outer_average'
595 H3D_KEYWORD_SOLID_SCALAR(I)%IS_SKIN = 1
596 H3D_KEYWORD_SOLID_SCALAR(I)%TEXT1 = 'fld failure factor outer_average'
597 H3D_KEYWORD_SOLID_SCALAR(I)%COMMENT = 'fld failure factor average of solid skin'
598
599 i = i + 1
600 h3d_keyword_solid_scalar(i)%KEY3 = 'FLDZ/OUTER_AVERAGE'
601 h3d_keyword_solid_scalar(i)%IS_SKIN = 1
602 h3d_keyword_solid_scalar(i)%TEXT1 = 'FLD Zone Index OUTER_AVERAGE'
603 h3d_keyword_solid_scalar(i)%COMMENT = '1:HighWrinkle 2:Compression 3:LooseMetal 4:SafeZone 5:MarginToFail 6:FailureZone'
604
605 i = i + 1
606 h3d_keyword_solid_scalar(i)%KEY3 = 'DIV(U)'
607 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric dilatation rate'
608 h3d_keyword_solid_scalar(i)%COMMENT = 'div(u)'
609
610 i = i + 1
611 h3d_keyword_solid_scalar(i)%KEY3 = 'ECONTROL'
612 h3d_keyword_solid_scalar(i)%TEXT1 = 'Internal Energy of distortion control'
613
614 i = i + 1
615 h3d_keyword_solid_scalar(i)%KEY3 = 'VSTRAIN'
616 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric Strain'
617 h3d_keyword_solid_scalar(i)%COMMENT = 'mu=rho/rho0-1'
618
619 i = i + 1
620 h3d_keyword_solid_scalar(i)%KEY3 = 'VSTRAIN/1'
621 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric Strain 1'
622 h3d_keyword_solid_scalar(i)%COMMENT = 'mu=rho/rho0-1'
623
624 i = i + 1
625 h3d_keyword_solid_scalar(i)%KEY3 = 'VSTRAIN/2'
626 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric Strain 2'
627 h3d_keyword_solid_scalar(i)%COMMENT = 'mu=rho/rho0-1'
628
629 i = i + 1
630 h3d_keyword_solid_scalar(i)%KEY3 = 'VSTRAIN/3'
631 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric Strain 3'
632 h3d_keyword_solid_scalar(i)%COMMENT = 'mu=rho/rho0-1'
633
634 i = i + 1
635 h3d_keyword_solid_scalar(i)%KEY3 = 'VSTRAIN/4'
636 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric Strain 4'
637 h3d_keyword_solid_scalar(i)%COMMENT = 'mu=rho/rho0-1'
638
639 i = i + 1
640 h3d_keyword_solid_scalar(i)%KEY3 = 'VSTRAIN/5'
641 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric Strain 5'
642 h3d_keyword_solid_scalar(i)%COMMENT = 'mu=rho/rho0-1'
643
644 i = i + 1
645 h3d_keyword_solid_scalar(i)%KEY3 = 'VSTRAIN/6'
646 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric Strain 6'
647 h3d_keyword_solid_scalar(i)%COMMENT = 'mu=rho/rho0-1'
648
649 i = i + 1
650 h3d_keyword_solid_scalar(i)%KEY3 = 'VSTRAIN/7'
651 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric Strain 7'
652 h3d_keyword_solid_scalar(i)%COMMENT = 'mu=rho/rho0-1'
653
654 i = i + 1
655 h3d_keyword_solid_scalar(i)%KEY3 = 'VSTRAIN/8'
656 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric Strain 8'
657 h3d_keyword_solid_scalar(i)%COMMENT = 'mu=rho/rho0-1'
658
659 i = i + 1
660 h3d_keyword_solid_scalar(i)%KEY3 = 'VSTRAIN/9'
661 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric Strain 9'
662 h3d_keyword_solid_scalar(i)%COMMENT = 'mu=rho/rho0-1'
663
664 i = i + 1
665 h3d_keyword_solid_scalar(i)%KEY3 = 'VSTRAIN/10'
666 h3d_keyword_solid_scalar(i)%TEXT1 = 'Volumetric Strain 10'
667 h3d_keyword_solid_scalar(i)%COMMENT = 'mu=rho/rho0-1'
668
669
670
671
672 nkey=i
end diagonal values have been computed in the(sparse) matrix id.SOL
subroutine jwl(iflag, nel, pm, off, eint, mu, mu2, espe, dvol, df, vnew, mat, psh, pnew, dpdm, dpde)
for(i8=*sizetab-1;i8 >=0;i8--)
subroutine tillotson(iflag, nel, pm, off, eint, mu, mu2, espe, dvol, df, vnew, mat, psh, pnew, dpdm, dpde, vareos, nvareos)