Numerical Simulation of the Dynamic Responses and Cumulative Damage of Underground Caverns under Multiple Explosions

Table 2

Parameters for Riedel–Hiermaier–Thoma (RHT) model of the rock mass.

Parameter	Value	Parameter	Value
Mass density RO (kg/m^3)	1700	Porosity exponent NP	0.0
Initial porosity ALPHA	0.0	Reference compressive strain rate EOC	3.0E−11
Crush pressure PEL (GPa)	0.0	Reference tensile strain rate EOT	3.0E−12
Compaction pressure PCO (GPa)	0.36	Break compressive strain rate EC	3.0E + 19
Hugoniot polynomial coefficient A1 (GPa)	2.1162	Break tensile strain rate ET	3.0E + 19
Hugoniot polynomial coefficient A2 (GPa)	2.3748	Compressive strain rate dependence exponent BETAC	0.032
Hugoniot polynomial coefficient A3 (GPa)	0.5424	Tensile strain rate dependence exponent BETAT	0.036
Parameter for polynomial EOS B0	1.22	Volumetric plastic strain fraction in tension PTF	0.001
Parameter for polynomial EOS B1	1.22	Compressive yield surface parameter GC	0.53
Parameter for polynomial EOS T1 (GPa)	2.1162	Tensile yield surface parameter GT	0.7
Parameter for polynomial EOS T2 (GPa)	0.0	Erosion plastic strain EPSF	2.0
Elastic shear modulus SHEAR (GPa)	1.002	Shear modulus reduction factor XI	0.5
Compressive strength FC (MPa)	2.1	Damage parameter D1	0.04
Relative tensile strength FT	0.1	Damage parameter D2	1.0
Relative shear strength FS	0.18	Minimum damaged residual strain EPM	0.01
Failure surface parameter A	1.6	Residual surface parameter AF	1.6
Failure surface parameter N	0.61	Residual surface parameter NF	0.61
Lode angle dependence factor Q0	0.6805	Gruneisen gamma GAMMA	0.0
Lode angle dependence factor B	0.0105	—