## 1 Introduction
In the last decade, the JH-series constitutive models have begun to be used in the research of rock materials involving the excavation of mining, tunnels, and quarries. [PERSON] and [PERSON] [26] simulated shock-induced damage in granite targets impacted by projectiles at different velocities using the JH-2 model. Additionally, [PERSON] and [PERSON] [27], [PERSON] et al. [28], [PERSON] and [PERSON] [29; 30], and [PERSON] and [PERSON] [31] simulated stress-wave-induced fracture in granite due to explosive action using the JH-2 model, and then compared the crack patterns and pressures with lab-scale blast experiments and empirical formulae. However, the simulation methodologies above still remained in two dimensions, which may limit the total presentation for the transmission of stress waves and affect the final blasting results. [PERSON] et al. [32] thus reproduced the rock dynamic fracture process of tunnel boring machine (TBM) tunneling in three dimensions. Moreover, the rock constitutive model was improved based on the JH-2 model by combing a Rankine tensile cracking softening model, however the validity of the improved model still needs to be deeply verified if employed in blasting rock tunnels. [PERSON] et al. [33; 34] incorporated the JHB model into the particle-based numerical manifold method (PNMM) to study the dynamic behavior of rock blasting by dual-level discretization system. However, the computational method is generally used to simulate single borehole blasting under dynamic loading conditions; if employed in large-scale and three-dimensional modeling, it would be hard to complete because of the limitation of the PNMM algorithm. [PERSON] et al. [35] proposed a coupled finite element method (FEM) and smoothed particle hydrodynamics (SPH) approach capable of reproducing the blast response in rock based on the JH-2 model. Different to the general FEM, the explosive was modelled explicitly using SPH and set as a node to surface contact for the interaction between rock and explosive materials. Although the FEM-SPH can overcome the excessive deformation of explosive, the uncertainty of the artificial contact would affect the results of damage. In general, most parametersof the simulations mentioned above lack complete descriptions about their resources or cited directly from other materials. Substantial static and dynamic experiments need to be conducted to completely describe one dynamic constitutive model, which would increase time and economic costs. If some simple and reliable computational or extrapolated methods can be supplied to obtain these parameters based on some fundamental experimental data, it would be of great significance and widely accepted. In this study, a parameter determination method for the JH-2 constitutive model for Barre granite is proposed and verified. The verified method together with the JH-2 model is applied in the simulation of mechanical behaviors and damage evolution induced by blasting for tunnel excavation. **2. JH-2 Constitutive Model for Rock Materials**  
The JH-2 constitutive model was initially utilized to simulate the behavior of brittle materials, especially ceramics. Based on its original concept presented in [2], the JH-2 model adds softening characteristics and contains pressure-dependent strength, damage, and fracture; significant strength after fracture; bulking; and strain rate effects. A general overview of the model is provided in Figures 1 and 2.