33
34
35
38 USE output_mod , ONLY : noda_surf, noda_pext
39
40
41
42#include "implicit_f.inc"
43
44
45
46#include "param_c.inc"
47#include "com04_c.inc"
48
49
50
51 TYPE (H3D_KEYWORD) H3D_KEYWORD_NODAL_SCALAR(*)
53 INTEGER, INTENT(INOUT) :: NKEY
54
55
56
57 INTEGER I,J,NIPMAX
58
59 i = 0
60
61 i = i + 1
62 h3d_keyword_nodal_scalar(i)%KEY3 = 'MASS'
63 h3d_keyword_nodal_scalar(i)%TEXT1 = 'Node Mass'
64
65 i = i + 1
66 h3d_keyword_nodal_scalar(i)%KEY3 = 'DT'
67 h3d_keyword_nodal_scalar(i)%TEXT1 = 'Nodal Time Step'
68
69 i = i + 1
70 h3d_keyword_nodal_scalar(i)%KEY3 = 'DMASS'
71 h3d_keyword_nodal_scalar(i)%TEXT1 = 'Mass Change'
72 h3d_keyword_nodal_scalar(i)%COMMENT = 'Nodal mass variation ( see time step control /DT/Eltyp/Iflag )'
73
74 i = i + 1
75 h3d_keyword_nodal_scalar(i)%KEY3 = 'P'
76 h3d_keyword_nodal_scalar(i)%TEXT1 = 'Nodal Pressure'
77 h3d_keyword_nodal_scalar(i)%COMMENT = 'Element pressure reported to node'
78
79 i = i + 1
80 h3d_keyword_nodal_scalar(i)%KEY3 = 'DENS'
81 h3d_keyword_nodal_scalar(i)%TEXT1 = 'Nodal Density'
82 h3d_keyword_nodal_scalar(i)%COMMENT = 'Element density reported to node'
83
84 i = i + 1
85 h3d_keyword_nodal_scalar(i)%KEY3 = 'ENER'
86 h3d_keyword_nodal_scalar(i)%TEXT1 = 'Nodal Specific Energy'
87 h3d_keyword_nodal_scalar(i)%COMMENT = 'Element specific energy reported to node'
88
89 i = i + 1
90 h3d_keyword_nodal_scalar(i)%KEY3 = 'temp'
91 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal temperature'
92
93 I = I + 1
94 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'nvar1'
95 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal variable 1'
96
97 I = I + 1
98 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'nvar2'
99 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal variable 2'
100
101 I = I + 1
102 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'nvar3'
103 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal variable 3'
104
105 I = I + 1
106 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'nvar4'
107 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal variable 4'
108
109 I = I + 1
110 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'nvar5'
111 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal variable 5'
112
113 I = I + 1
114 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'diner'
115 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'inertia change'
116 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'output of added inertia per nodes'
117
118 i = i + 1
119 h3d_keyword_nodal_scalar(i)%KEY3 = 'fpot'
120 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal potential'
121
122 I = I + 1
123 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'ndmass'
124 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'non diagonal mass change'
125 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'non-dagonal mass variation(see time step control /dt/ams)'
126
127 I = I + 1
128 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'dama2/normal'
129 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = '%damage(type2 interface) / normal'
130
131 I = I + 1
132 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'dama2/tangent'
133 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = '%damage(type2 interface) / tangent'
134
135 I = I + 1
136 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'schli'
137 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal schlieren'
138
139 I = I + 1
140 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'stifr'
141 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal rotational stiffness'
142
143 I = I + 1
144 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'stif'
145 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal stiffness'
146
147 I = I + 1
148 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'vfrac/1'
149 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal volumetric fraction - 1'
150 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'element volumetric fraction reported to node, used with law37 and law51'
151
152 I = I + 1
153 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'vfrac/2'
154 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal volumetric fraction - 2'
155 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'element volumetric fraction reported to node, used with law37 and law51'
156
157 I = I + 1
158 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'vfrac/3'
159 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal volumetric fraction - 3'
160 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'element volumetric fraction reported to node, used with law37 and law51'
161
162 I = I + 1
163 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'vfrac/4'
164 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal volumetric fraction - 4'
165 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'element volumetric fraction reported to node, used with law37 and law51'
166
167 I = I + 1
168 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'zvfrac/1'
169 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal volumetric fraction - 1'
170
171 I = I + 1
172 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'zvfrac/2'
173 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal volumetric fraction - 2'
174
175 I = I + 1
176 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'zvfrac/3'
177 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal volumetric fraction - 3'
178
179 I = I + 1
180 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'zvfrac/4'
181 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal volumetric fraction - 4'
182
183 I = I + 1
184 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'zvol'
185 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'centroid new volume'
186
187 I = I + 1
188 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'zvold'
189 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'centroid old volume'
190
191 I = I + 1
192 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'gps1/p'
193 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'gps1 pressure'
194 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'mean pressure of elements connected to node'
195
196 I = I + 1
197 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'gps1/vonm'
198 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'gps1 von mises'
199 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'mean von mises of elements connected to node'
200
201 I = I + 1
202 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'skid_line'
203 H3D_KEYWORD_NODAL_SCALAR(I)%IS_INTER = 1
204 H3D_KEYWORD_NODAL_SCALAR(I)%IS_INTER_ALL = 1
205 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'skid line
for Interface '
206 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'the ratio between maximum tangential force and s*pmax/sqrt3
'
207
208 I = I + 1
209 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'ssp'
210 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal sound speed'
211 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'element sound speed reported to node'
212
213 I = I + 1
214 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'internal.
id'
215 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'internal identifier'
216
217 I = I + 1
218 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'cse_fric'
219 H3D_KEYWORD_NODAL_SCALAR(I)%IS_INTER = 1
220 H3D_KEYWORD_NODAL_SCALAR(I)%IS_INTER_ALL = 1
221 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'surface frictional contact energy
for interface
'
222
223 I = I + 1
224 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'cse_fricg'
225 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'surface frictional contact energy'
226 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'node energy density(sum of all contacts) '
227
228 I = I + 1
229 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'vol'
230 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal volumes'
231 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'volume at element centroid(fvmbag)'
232
233 I = I + 1
234 H3D_KEYWORD_NODAL_SCALAR(I)%KEY3 = 'pext'
235 H3D_KEYWORD_NODAL_SCALAR(I)%TEXT1 = 'nodal External pressure'
236 H3D_KEYWORD_NODAL_SCALAR(I)%COMMENT = 'pressure applied by load case'
237
238
239 NKEY=I
end diagonal values have been computed in the(sparse) matrix id.SOL
for(i8=*sizetab-1;i8 >=0;i8--)