OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
h3d_list_shell_scalar.F
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!||====================================================================
24!|| h3d_list_shell_scalar ../engine/source/output/h3d/input_list/h3d_list_shell_scalar.F
25!||--- called by ------------------------------------------------------
26!|| lech3d ../engine/source/output/h3d/h3d_build_fortran/lech3d.f
27!||--- uses -----------------------------------------------------
28!|| h3d_mod ../engine/share/modules/h3d_mod.F
29!|| message_mod ../engine/share/message_module/message_mod.F
30!|| multi_fvm_mod ../common_source/modules/ale/multi_fvm_mod.F90
31!|| stack_mod ../engine/share/modules/stack_mod.F
32!||====================================================================
33 SUBROUTINE h3d_list_shell_scalar(H3D_KEYWORD_SHELL_SCALAR,GEO,IGEO,MULTI_FVM,NKEY)
34C-----------------------------------------------
35C M o d u l e s
36C-----------------------------------------------
37 USE message_mod
38 USE h3d_mod
39 USE stack_mod
40 USE multi_fvm_mod
41C-----------------------------------------------
42C I m p l i c i t T y p e s
43C-----------------------------------------------
44#include "implicit_f.inc"
45C-----------------------------------------------
46C C o m m o n B l o c k s
47C-----------------------------------------------
48#include "param_c.inc"
49#include "com04_c.inc"
50C-----------------------------------------------
51C D u m m y A r g u m e n t s
52C-----------------------------------------------
53 TYPE (H3D_KEYWORD) H3D_KEYWORD_SHELL_SCALAR(*)
54 my_real geo(npropg,numgeo)
55 INTEGER IGEO(NPROPGI,NUMGEO)
56 TYPE(multi_fvm_struct), INTENT(IN) :: MULTI_FVM
57 INTEGER, INTENT(INOUT) :: NKEY
58C-----------------------------------------------
59C L o c a l V a r i a b l e s
60C-----------------------------------------------
61 INTEGER I,II
62C=========================================================================
63 i = 0
64c-----------------------------------------------
65 i = i + 1
66 h3d_keyword_shell_scalar(i)%KEY3 = 'MASS'
67 h3d_keyword_shell_scalar(i)%TEXT1 = 'Mass'
68c-----------------------------------------------
69 i = i + 1
70 h3d_keyword_shell_scalar(i)%KEY3 = 'DENS'
71 h3d_keyword_shell_scalar(i)%TEXT1 = 'Density'
72c-----------------------------------------------
73 i = i + 1
74 h3d_keyword_shell_scalar(i)%KEY3 = 'ENER'
75 h3d_keyword_shell_scalar(i)%TEXT1 = 'Internal Energy per unit mass'
76c-----------------------------------------------
77 i = i + 1
78 h3d_keyword_shell_scalar(i)%KEY3 = 'EINTM'
79 h3d_keyword_shell_scalar(i)%TEXT1 = 'Internal Energy per unit mass'
80c-----------------------------------------------
81 i = i + 1
82 h3d_keyword_shell_scalar(i)%KEY3 = 'EINTV'
83 h3d_keyword_shell_scalar(i)%TEXT1 = 'Internal Energy per unit volume'
84c-----------------------------------------------
85 i = i + 1
86 h3d_keyword_shell_scalar(i)%KEY3 = 'EINT'
87 h3d_keyword_shell_scalar(i)%TEXT1 = 'Internal Energy'
88c-----------------------------------------------
89 i = i + 1
90 h3d_keyword_shell_scalar(i)%KEY3 = 'TEMP'
91 h3d_keyword_shell_scalar(i)%TEXT1 = 'Temperature'
92c-----------------------------------------------
93 i = i + 1
94 h3d_keyword_shell_scalar(i)%KEY3 = 'THICK'
95 h3d_keyword_shell_scalar(i)%TEXT1 = 'Thickness'
96c-----------------------------------------------
97 i = i + 1
98 h3d_keyword_shell_scalar(i)%KEY3 = 'P'
99 h3d_keyword_shell_scalar(i)%TEXT1 = 'Pressure'
100c-----------------------------------------------
101 i = i + 1
102 h3d_keyword_shell_scalar(i)%KEY3 = 'VONM'
103 h3d_keyword_shell_scalar(i)%TEXT1 = 'Von Mises'
104c-----------------------------------------------
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'
110c-----------------------------------------------
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'
115c-----------------------------------------------
116 i = i + 1
117 h3d_keyword_shell_scalar(i)%KEY3 = 'VORTX'
118 h3d_keyword_shell_scalar(i)%TEXT1 = 'Vorticity'
119c-----------------------------------------------
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'
130c-----------------------------------------------
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'
141c-----------------------------------------------
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'
152c-----------------------------------------------
153 i = i + 1
154 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGX'
155 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress X'
156c-----------------------------------------------
157 i = i + 1
158 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGY'
159 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress Y'
160c-----------------------------------------------
161 i = i + 1
162 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGZ'
163 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress Z'
164c-----------------------------------------------
165 i = i + 1
166 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGXY'
167 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress XY'
168c-----------------------------------------------
169 i = i + 1
170 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGYZ'
171 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress YZ'
172c-----------------------------------------------
173 i = i + 1
174 h3d_keyword_shell_scalar(i)%KEY3 = 'SIGZX'
175 h3d_keyword_shell_scalar(i)%TEXT1 = 'Stress ZX'
176c-----------------------------------------------
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'
188c-----------------------------------------------
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'
198c-----------------------------------------------
199 i = i + 1
200 h3d_keyword_shell_scalar(i)%KEY3 = 'HOURGLASS'
201 h3d_keyword_shell_scalar(i)%TEXT1 = 'Hourglass Energy per unit mass'
202c-----------------------------------------------
203 IF (multi_fvm%IS_USED) THEN
204c
205c Number of present material is dynamic for law151, hence create full list
206c of output keywords at runtime by looping over all present materials
207c
208 DO ii = 1, multi_fvm%NBMAT
209c-----------------------------------------------
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)'
214c-----------------------------------------------
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)'
219c-----------------------------------------------
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)'
224c-----------------------------------------------
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)'
229c-----------------------------------------------
230 ENDDO
231c-----------------------------------------------
232 i = i + 1
233 h3d_keyword_shell_scalar(i)%KEY3 = 'ENTH'
234 h3d_keyword_shell_scalar(i)%TEXT1 = 'Enthalpy'
235c-----------------------------------------------
236 i = i + 1
237 h3d_keyword_shell_scalar(i)%KEY3 = 'ENTHM'
238 h3d_keyword_shell_scalar(i)%TEXT1 = 'Enthalpy per unit mass'
239c-----------------------------------------------
240 i = i + 1
241 h3d_keyword_shell_scalar(i)%KEY3 = 'ENTHV'
242 h3d_keyword_shell_scalar(i)%TEXT1 = 'Enthalpy per unit volume'
243c-----------------------------------------------
244 ENDIF ! MULTI_FVM%IS_USED
245c-----------------------------------------------
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)'
250c-----------------------------------------------
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 '
260c-----------------------------------------------
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'
266c-----------------------------------------------
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 '
272c-----------------------------------------------
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)'
280c-----------------------------------------------
281 i = i + 1
282 h3d_keyword_shell_scalar(i)%KEY3 = 'CRMAS'
283 h3d_keyword_shell_scalar(i)%TEXT1 = 'Airbag crossing mass'
284c-----------------------------------------------
285 i = i + 1
286 h3d_keyword_shell_scalar(i)%KEY3 = 'CRVEL'
287 h3d_keyword_shell_scalar(i)%TEXT1 = 'airbag crossing velocity'
288c-----------------------------------------------
289 I = I + 1
290 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'aleba/mass'
291 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fvmbag - mass'
292c-----------------------------------------------
293 I = I + 1
294 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'aleba/pres'
295 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fvmbag - pressure'
296c-----------------------------------------------
297 I = I + 1
298 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'aleba/momx'
299 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fvmbag - fluid velocity x'
300c-----------------------------------------------
301 I = I + 1
302 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'aleba/momy'
303 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fvmbag - fluid velocity y'
304c-----------------------------------------------
305 I = I + 1
306 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'aleba/momz'
307 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fvmbag - fluid velocity z'
308c-----------------------------------------------
309 I = I + 1
310 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'aleba/dens'
311 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fvmbag - density'
312c-----------------------------------------------
313 I = I + 1
314 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'aleba/ener'
315 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fvmbag - energy'
316c-----------------------------------------------
317 I = I + 1
318 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'aleba/ssp'
319 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fvmbag - sound speed'
320c-----------------------------------------------
321 I = I + 1
322 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'aleba/gama'
323 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fvmbag - gama'
324c-----------------------------------------------
325 I = I + 1
326 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'aleba/visu'
327 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'fvmbag - visu polyhedra'
328c-----------------------------------------------
329 I = I + 1
330 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'thin'
331 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = '%Thinning'
332c-----------------------------------------------
333 I = I + 1
334 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'error/thick'
335 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'estimated error on thickness'
336c-----------------------------------------------
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'
348c-----------------------------------------------
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'
359c-----------------------------------------------
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'
372c-----------------------------------------------
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'
383c-----------------------------------------------
384 I = I + 1
385 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'dama/memb'
386 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'max damage membrane'
387c-----------------------------------------------
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'
401c-----------------------------------------------
402 I = I + 1
403 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'damg/memb'
404 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'damage(membrane value)'
405c-----------------------------------------------
406 I = I + 1
407 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = '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'
413c-----------------------------------------------
414 I = I + 1
415 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'nxtf/memb'
416 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'nxt failure factor membrane'
417c-----------------------------------------------
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.'
422c-----------------------------------------------
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'
427c-----------------------------------------------
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'
432c-----------------------------------------------
433 I = I + 1
434 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'phi/memb'
435 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'phi membrane'
436c-----------------------------------------------
437 I = I + 1
438 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'domain'
439 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'domain'
440c-----------------------------------------------
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'
445c-----------------------------------------------
446 I = I + 1
447 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'bulk'
448 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'artificial viscosity'
449c-----------------------------------------------
450 I = I + 1
451 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'tdet'
452 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'detonation time'
453c-----------------------------------------------
454 I = I + 1
455 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'bfrac'
456 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'burn fraction'
457c-----------------------------------------------
458 I = I + 1
459 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'dt'
460 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'time step'
461c-----------------------------------------------
462 I = I + 1
463 H3D_KEYWORD_SHELL_SCALAR(I)%KEY3 = 'ams'
464 H3D_KEYWORD_SHELL_SCALAR(I)%TEXT1 = 'ams selection'
465c-----------------------------------------------
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 - change of angle between two anisotropy directions'
472c-----------------------------------------------
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'
477c-----------------------------------------------
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'
483c-----------------------------------------------
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'
489c-----------------------------------------------
490 i = i + 1
491 h3d_keyword_shell_scalar(i)%KEY3 = 'FLDF/MEMB'
492 h3d_keyword_shell_scalar(i)%TEXT1 = 'FLD Failure Factor Membrane'
493c-----------------------------------------------
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
501c-----------------------------------------------
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'
506c-----------------------------------------------
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'
515c-----------------------------------------------
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'
519c-----------------------------------------------
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
529c-----------------------------------------------
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 and below HC 3: Above HC and DSSE'
534c-----------------------------------------------
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'
545c-----------------------------------------------
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'
555c-----------------------------------------------
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
568c-----------------------------------------------
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'
574c-----------------------------------------------
575 i = i + 1
576 h3d_keyword_shell_scalar(i)%KEY3 = 'GROUP'
577 h3d_keyword_shell_scalar(i)%TEXT1 = 'Group identifier'
578c-----------------------------------------------
579 i = i + 1
580 h3d_keyword_shell_scalar(i)%KEY3 = 'INTERNAL.ID'
581 h3d_keyword_shell_scalar(i)%TEXT1 = 'Internal identifier'
582c-----------------------------------------------
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'
586c-----------------------------------------------
587 i = i + 1
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'
591c-----------------------------------------------
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'
596c-----------------------------------------------
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 = 1
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'
607c-----------------------------------------------
608 i = i + 1
609 h3d_keyword_shell_scalar(i)%KEY3 = 'DAMA/TMAX'
610 h3d_keyword_shell_scalar(i)%IS_LAYER = 1
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'
618c-----------------------------------------------
619 !2d solids (triangles)
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'
624c-----------------------------------------------
625c-----------------------------------------------
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'
630c-----------------------------------------------
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'
635c-----------------------------------------------
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'
640c-----------------------------------------------
641 i = i + 1
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'
645c-----------------------------------------------
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'
650c-----------------------------------------------
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'
655c-----------------------------------------------
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'
660c-----------------------------------------------
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'
665c-----------------------------------------------
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'
670c-----------------------------------------------
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'
675c-----------------------------------------------
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'
680c-----------------------------------------------
681
682 NKEY=I
683 END
my_real
#define my_real
Definition cppsort.cpp:32
the
end diagonal values have been computed in the(sparse) matrix id.SOL
alpha
#define alpha
Definition eval.h:35
h3d_list_shell_scalar
subroutine h3d_list_shell_scalar(h3d_keyword_shell_scalar, geo, igeo, multi_fvm, nkey)
Definition h3d_list_shell_scalar.F:34
lech3d
subroutine lech3d(output, geo, igeo, ipm, ipart, h3d_data, multi_fvm, ipari, iparg, tag_skins6, mds_label, mds_output_table, mds_nmat, max_depvar, mds_ndepsvar, elbuf_str, stack, ibcl, iloadp, lloadp, loads, mat_param, pblast, igrpart, npc, tf, snpc, stf)
Definition lech3d.F:97
max
#define max(a, b)
Definition macros.h:21
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
h3d_mod
Definition h3d_mod.F:311
message_mod
Definition message_mod.F:1257
stack_mod
Definition stack_mod.F:75
velocity
subroutine velocity(a, ar, v, vr, fzero, itab, nale)
Definition velocity.F:29