40#include "implicit_f.inc"
44 TYPE (H3D_KEYWORD) H3D_KEYWORD_SPH_SCALAR(*)
46 INTEGER IGEO(NPROPGI,NUMGEO)
47 INTEGER,
INTENT(INOUT) :: NKEY
63 h3d_keyword_sph_scalar(i)%KEY3 =
'MASS'
64 h3d_keyword_sph_scalar(i)%TEXT1 =
'Mass'
67 h3d_keyword_sph_scalar(i)%KEY3 =
'DIAMETER'
68 h3d_keyword_sph_scalar(i)%TEXT1 =
'Diameter'
71 h3d_keyword_sph_scalar(i)%KEY3 =
'EPSP'
72 h3d_keyword_sph_scalar(i)%TEXT1 =
'Plastic Strain'
75 h3d_keyword_sph_scalar(i)%KEY3 =
'DENS'
76 h3d_keyword_sph_scalar(i)%TEXT1 =
'Density'
79 h3d_keyword_sph_scalar(i)%KEY3 =
'ENER'
80 h3d_keyword_sph_scalar(i)%TEXT1 =
'Internal Energy per unit mass'
83 h3d_keyword_sph_scalar(i)%KEY3 =
'TEMP'
84 h3d_keyword_sph_scalar(i)%TEXT1 =
'Temperature'
87 h3d_keyword_sph_scalar(i)%KEY3 =
'P'
88 h3d_keyword_sph_scalar(i)%TEXT1 =
'Pressure'
91 h3d_keyword_sph_scalar(i)%KEY3 =
'VONM'
92 h3d_keyword_sph_scalar(i)%TEXT1 =
'Von Mises'
95 h3d_keyword_sph_scalar(i)%KEY3 =
'K'
96 h3d_keyword_sph_scalar(i)%TEXT1 =
'Turbulent Energy'
99 h3d_keyword_sph_scalar(i)%KEY3 =
'TVIS'
100 h3d_keyword_sph_scalar(i)%TEXT1 =
'Turbulent Viscosity'
103 h3d_keyword_sph_scalar(i)%KEY3 =
'VORTX'
104 h3d_keyword_sph_scalar(i)%TEXT1 =
'Vorticity'
107 h3d_keyword_sph_scalar(i)%KEY3 =
'DAM1'
108 h3d_keyword_sph_scalar(i)%TEXT1 =
'Damage 1'
111 h3d_keyword_sph_scalar(i)%KEY3 =
'DAM2'
112 h3d_keyword_sph_scalar(i)%TEXT1 =
'Damage 2'
115 h3d_keyword_sph_scalar(i)%KEY3 =
'DAM3'
116 h3d_keyword_sph_scalar(i)%TEXT1 =
'Damage 3'
119 h3d_keyword_sph_scalar(i)%KEY3 =
'SIGX'
120 h3d_keyword_sph_scalar(i)%TEXT1 =
'Stress X'
123 h3d_keyword_sph_scalar(i)%KEY3 =
'SIGY'
124 h3d_keyword_sph_scalar(i)%TEXT1 =
'Stress Y'
127 h3d_keyword_sph_scalar(i)%KEY3 =
'SIGZ'
128 h3d_keyword_sph_scalar(i)%TEXT1 =
'Stress Z'
131 h3d_keyword_sph_scalar(i)%KEY3 =
'SIGXY'
132 h3d_keyword_sph_scalar(i)%TEXT1 =
'Stress XY'
135 h3d_keyword_sph_scalar(i)%KEY3 =
'SIGYZ'
136 h3d_keyword_sph_scalar(i)%TEXT1 =
'Stress YZ'
139 h3d_keyword_sph_scalar(i)%KEY3 =
'SIGZX'
140 h3d_keyword_sph_scalar(i)%TEXT1 =
'Stress ZX'
143 h3d_keyword_sph_scalar(i)%KEY3 =
'USER'
144 h3d_keyword_sph_scalar(i)%IS_UVAR = 1
145 h3d_keyword_sph_scalar(i)%IS_UVAR_ALL = 1
146 h3d_keyword_sph_scalar(i)%TEXT1 =
'User Var'
149 h3d_keyword_sph_scalar(i)%KEY3 =
'HOURGLASS'
150 h3d_keyword_sph_scalar(i)%TEXT1 =
'Hourglass Energy'
153 h3d_keyword_sph_scalar(i)%KEY3 =
'BFRAC'
154 h3d_keyword_sph_scalar(i)%TEXT1 =
'Burn Fraction'
157 h3d_keyword_sph_scalar(i)%KEY3 =
'DAMA'
158 h3d_keyword_sph_scalar(i)%TEXT1 =
'MAX DAMAGE'
159 h3d_keyword_sph_scalar(i)%COMMENT =
'maximum of damage over time and of all failure criteria acting in one material'
162 h3d_keyword_sph_scalar(i)%KEY3 =
'DOMAIN'
163 h3d_keyword_sph_scalar(i)%TEXT1 =
'Domain'
166 h3d_keyword_sph_scalar(i)%KEY3 =
'FILL'
167 h3d_keyword_sph_scalar(i)%TEXT1 =
'Filling percentage'
170 h3d_keyword_sph_scalar(i)%KEY3 =
'SIGEQ'
171 h3d_keyword_sph_scalar(i)%TEXT1 =
'Equiv stress'
174 h3d_keyword_sph_scalar(i)%KEY3 =
'TDET'
175 h3d_keyword_sph_scalar(i)%TEXT1 = 'detonation time
'
176 H3D_KEYWORD_SPH_SCALAR(I)%COMMENT = 'only available
for jwl eos or time controlled detonation
'
179 H3D_KEYWORD_SPH_SCALAR(I)%KEY3 = 'off
'
180 H3D_KEYWORD_SPH_SCALAR(I)%TEXT1 = 'element status
'
181 H3D_KEYWORD_SPH_SCALAR(I)%COMMENT = '0.0:deleted >0.&<1.:under failure process 1.0:activated <0:standby
'
184 H3D_KEYWORD_SPH_SCALAR(I)%KEY3 = 'enth
'
185 H3D_KEYWORD_SPH_SCALAR(I)%TEXT1 = 'enthalpy
'
188 H3D_KEYWORD_SPH_SCALAR(I)%KEY3 = 'enthm
'
189 H3D_KEYWORD_SPH_SCALAR(I)%TEXT1 = 'enthalpy per unit mass
'
192 H3D_KEYWORD_SPH_SCALAR(I)%KEY3 = 'enthv
'
193 H3D_KEYWORD_SPH_SCALAR(I)%TEXT1 = 'enthalpy per unit volume
'
196 H3D_KEYWORD_SPH_SCALAR(I)%KEY3 = 'eint
'
197 H3D_KEYWORD_SPH_SCALAR(I)%TEXT1 = 'internal energy
'
200 H3D_KEYWORD_SPH_SCALAR(I)%KEY3 = 'eintm
'
201 H3D_KEYWORD_SPH_SCALAR(I)%TEXT1 = 'internal energy per unit mass
'
204 H3D_KEYWORD_SPH_SCALAR(I)%KEY3 = 'eintv
'
205 H3D_KEYWORD_SPH_SCALAR(I)%TEXT1 = 'internal energy per unit volume
'
208 H3D_KEYWORD_SPH_SCALAR(I)%KEY3 = 'group
'
209 H3D_KEYWORD_SPH_SCALAR(I)%TEXT1 = 'group identifier
'
212 H3D_KEYWORD_SPH_SCALAR(I)%KEY3 = 'internal.
id'
213 H3D_KEYWORD_SPH_SCALAR(I)%TEXT1 = 'internal identifier
'
216 H3D_KEYWORD_SPH_SCALAR(I)%KEY3 = 'local.
id'
217 H3D_KEYWORD_SPH_SCALAR(I)%TEXT1 = 'local identifier within
the group
'
220 H3D_KEYWORD_SPH_SCALAR(I)%KEY3 = 'tillotson'
221 h3d_keyword_sph_scalar(i)%TEXT1 =
'Region identifier in p,v diagram'
222 h3d_keyword_sph_scalar(i)%COMMENT =
'1:compressed, 2:expanded cold, 3:transition, 4:hot expansion'
225 h3d_keyword_sph_scalar(i)%KEY3 =
'FAILURE'
226 h3d_keyword_sph_scalar(i)%IS_ID = 1
227 h3d_keyword_sph_scalar(i)%IS_MODE = 1
228 h3d_keyword_sph_scalar(i)%TEXT1 =
'Failure'
229 h3d_keyword_sph_scalar(i)%COMMENT =
'Damage of current failure criterion'
232 h3d_keyword_sph_scalar(i)%KEY3 =
'NEIGHBOURS'
233 h3d_keyword_sph_scalar(i)%TEXT1 =
'Number of Neighbours'
236 h3d_keyword_sph_scalar(i)%KEY3 =
'VSTRAIN'
237 h3d_keyword_sph_scalar(i)%TEXT1 =
'Volumetric Strain'
238 h3d_keyword_sph_scalar(i)%COMMENT =
'mu=rho/rho0-1'
subroutine jwl(iflag, nel, pm, off, eint, mu, mu2, espe, dvol, df, vnew, mat, psh, pnew, dpdm, dpde)
subroutine tillotson(iflag, nel, pm, off, eint, mu, mu2, espe, dvol, df, vnew, mat, psh, pnew, dpdm, dpde, vareos, nvareos)