Mobilized Mohr-Coulomb and Hoek-Brown Strength Parameters during Failure of Granite in Alxa Area in China for High-Level Radioactive Waste Disposal

: Strength parameters of the host rock is of paramount importance for modelling the behaviors of underground disposal repository of high-level radioactive waste (HLW). Mobilization of strength parameters should be studied for a better understanding and modelling on the mechanical behaviors of the surrounding rock, considering the e ﬀ ect of temperature induced by the nuclear waste. The granite samples cored from NRG01 borehole in Alxa candidate area in China for HLW disposal are treated by di ﬀ erent temperatures (T = 20 ◦ C, 100 ◦ C and 200 ◦ C), and then are used to carry out a series of uniaxial and tri-axial compression experiments under various conﬁning pressures ( σ 3 = 0, 5, 10, 20, and 30 MPa) in this study. With the recorded axial stress—axial strain and axial stress—lateral strain curves, mobilization of both Mohr-Coulomb and Hoek-Brown strength parameters are analyzed with the increasing plastic shear strain. It has been found that NRG01 granite samples show generally similar cohesion weakening and friction strengthening behaviors, as well as the non-simultaneous mobilization of Hoek-Brown strength parameters ( m b and s ), under the e ﬀ ect of various treatment temperatures. Furthermore, the samples treated by higher temperatures show lower initial values of cohesion, but their initial friction angle and m b values are relatively higher. This should be mainly owing to the thermally induced cracks in the samples. This study should be helpful for a better modelling on the mechanical behaviors of NRG01 granite samples as the host rock of a possible HLW disposal repository.


Introduction
Granite is considered as one of the most important types of host rock for geological disposal of high-level radioactive waste (HLW) [1][2][3][4][5][6][7][8][9]. An appropriate modelling on the mechanical behavior of granite is of great importance for site selection and design of the repository [9][10][11][12][13]. Specially, it should be noted that the heat induced by the nuclear waste may have considerable influences on the mechanical behavior of the host rock, so the thermal effect cannot be ignored [8,12,[14][15][16][17].
There have been extensive studies on the mechanical behaviors of the host rock for HLW disposal [1,9,[18][19][20][21][22]. In many studies, simultaneously mobilized Mohr-Coulomb strength parameters (cohesion and friction angle) were used in the modelling of the underground excavation [23][24][25][26][27]. In order to describe the plastic strain softening behavior of the rock, they assumed that both cohesion and friction angle degrade from the initial value to the residual value with the increasing plastic strain, and piecewise linear models were usually adopted [23][24][25][26][27]. However, based on a series of  [46]. TMS01, TMS02, NRG01 and NRG02 are four boreholes drilled in TMS and NRG sub-areas.

Samples and Experimental Methods
NRG01 granite samples treated by different temperatures are used to conduct a series of uniaxial and tri-axial compression tests with various confining pressures. The obtained stress-strain data will be used for analyzing the mobilization of strength components during the brittle failure of granite considering the thermal effect.

Samples
The granite samples are cored from NRG01 borehole, which is one of the four 600 m-deep boreholes in Alxa area. According to the field investigations on the corresponding outcrops, RQD analyses on the drilling cores, as well as the mechanical experiments on the cored specimens in laboratory, NRG01 samples show the best structural and strength quality and thus are selected to be used for further studies [1].
The pink samples are cored from the depth of 500-600 m. The typical samples are presented in Figure 2a,b. It can be found that the samples are heterogeneous and have coarse particles. According to the observation on thin sections under polarized microscopy, the mineral contents and the grain sizes of NRG01 granite samples are analyzed and listed in Table 1 [47]. Based on the mineral components, the samples should be named as biotite syenogranite. Nevertheless, they are still called as granite samples in this paper for simplicity. Figure 2c shows a comparison on the strength values of different granite samples under various confining pressures. Apparently, NRG01 samples have higher strength than TMS01 granite samples cored from TMS01 borehole in TMS sub-area (shown in Figure 1) of Alxa candidate area. Compared with BS06 granite samples cored from Beishan candidate area in Gansu Province [48], NRG01 samples show a little lower strength under lower confinements  [46]. TMS01, TMS02, NRG01 and NRG02 are four boreholes drilled in TMS and NRG sub-areas.

Samples and Experimental Methods
NRG01 granite samples treated by different temperatures are used to conduct a series of uniaxial and tri-axial compression tests with various confining pressures. The obtained stress-strain data will be used for analyzing the mobilization of strength components during the brittle failure of granite considering the thermal effect.

Samples
The granite samples are cored from NRG01 borehole, which is one of the four 600 m-deep boreholes in Alxa area. According to the field investigations on the corresponding outcrops, RQD analyses on the drilling cores, as well as the mechanical experiments on the cored specimens in laboratory, NRG01 samples show the best structural and strength quality and thus are selected to be used for further studies [1].
The pink samples are cored from the depth of 500-600 m. The typical samples are presented in Figure 2a,b. It can be found that the samples are heterogeneous and have coarse particles. According to the observation on thin sections under polarized microscopy, the mineral contents and the grain sizes of NRG01 granite samples are analyzed and listed in Table 1 [47]. Based on the mineral components, the samples should be named as biotite syenogranite. Nevertheless, they are still called as granite samples in this paper for simplicity. Figure 2c shows a comparison on the strength values of different granite samples under various confining pressures. Apparently, NRG01 samples have higher strength than TMS01 granite samples cored from TMS01 borehole in TMS sub-area (shown in Figure 1) of Alxa candidate area. Compared with BS06 granite samples cored from Beishan candidate Hoek-Brown criterion is also used to analyze the data, and the non-linear Hoek-Brown fitting curve is shown in Figure 2c. Hoek-Brown criterion is also used to analyze the data, and the non-linear Hoek-Brown fitting curve is shown in Figure 2c.

Experimental Methods
A series of cylindrical NRG01 granite samples are well prepared (listed in Table 2) for uniaxial and tri-axial compression experiments under various confining pressures (σ3 = 0, 5, 10, 20, and 30 MPa). Concerning the heat produced by the high-level radioactive waste during the long-term disposal period, the effect of temperature should also be considered in this study. According to an extensive review on the conceptual design of repositories [7,[49][50][51][52], the temperature applied on the host rock will be no higher than 100 °C-120 °C. Consequently, this study focuses on the range from 20 °C (room temperature) to 200 °C.

Experimental Methods
A series of cylindrical NRG01 granite samples are well prepared (listed in Table 2) for uniaxial and tri-axial compression experiments under various confining pressures (σ 3 = 0, 5, 10, 20, and 30 MPa). Concerning the heat produced by the high-level radioactive waste during the long-term disposal period, the effect of temperature should also be considered in this study. According to an extensive review on the conceptual design of repositories [7,[49][50][51][52], the temperature applied on the host rock will be no higher than 100 • C-120 • C. Consequently, this study focuses on the range from 20 • C (room temperature) to 200 • C.
The specimens are firstly heated in a heating cabinet to the designed temperatures as shown in Table 2. The heating rate is set as 2 • C/min. When the target temperatures are reached, the heat treated samples are used for a series of uniaxial and tri-axial compression experiments with the TAW2000 servo-control tri-axial compression test system in Key Laboratory of Shale Gas and Geoengineeirng, Chinese Academy of Sciences. It should be noted that the rock specimens cannot remain their treatment temperatures as there is not a heating system during the compression tests. The confining pressures are applied to the target values as presented in Table 2, followed by the axial loading at a constant strain rate of 1.0 × 10 −5 ·s −1 . During each test, the axial and lateral strain are both measured with a set of extensometers, and the axial stress is obtained according to the axial load monitored by a force sensor. Consequently, the axial stress -axial strain curve and axial stress-lateral strain curve can be recorded for each test, and the failure characteristics of the specimens will also be observed after the experiments are completed.

Experimental Results
The differential stress-axial strain curves and differential stress-lateral strain curves for all the tests are presented in Figure 3a-e. It should be noted that differential stress (σ 1 − σ 3 ) is used in these curves in order for a more consistent observation. The peak strength values and the strength envelopes fitted with Hoek-Brown criterion are shown in Figure 3f. Based on these test results, some features can be observed as follows: (1) For the NRG01 granite samples treated by different temperatures (T = 20 • C, 100 • C and 200 • C), the differential stress -axial strain curves show the similar brittle -ductile transition behaviors with the increasing confining pressures (σ 3 = 0-30 MPa); (2) According to the experimental results, the heat treatment by temperatures no higher than 200 • C does not have very significant influence on the stress -strain curves of the NRG01 granite samples under various confining pressures (σ 3 = 0-30 MPa). However, if we make a more careful observation, it can be found that the samples treated by higher temperature show relatively more ductile behavior during the post-peak stage; The specimens are firstly heated in a heating cabinet to the designed temperatures as shown in Table 2. The heating rate is set as 2 °C/min. When the target temperatures are reached, the heat treated samples are used for a series of uniaxial and tri-axial compression experiments with the TAW2000 servo-control tri-axial compression test system in Key Laboratory of Shale Gas and Geoengineeirng, Chinese Academy of Sciences. It should be noted that the rock specimens cannot remain their treatment temperatures as there is not a heating system during the compression tests. The confining pressures are applied to the target values as presented in Table 2, followed by the axial loading at a constant strain rate of 1.0 × 10 −5 •s −1 . During each test, the axial and lateral strain are both measured with a set of extensometers, and the axial stress is obtained according to the axial load monitored by a force sensor. Consequently, the axial stress -axial strain curve and axial stress-lateral strain curve can be recorded for each test, and the failure characteristics of the specimens will also be observed after the experiments are completed.

Experimental Results
The differential stress-axial strain curves and differential stress-lateral strain curves for all the tests are presented in Figure 3a-e. It should be noted that differential stress (σ1 − σ3) is used in these curves in order for a more consistent observation. The peak strength values and the strength envelopes fitted with Hoek-Brown criterion are shown in Figure 3f. Based on these test results, some features can be observed as follows:

Analytical Method
According to the previous studies [26,28,31,33,53], the Mohr-Coulomb strength parameters (cohesion c and inner frictional angle φ) of rock should be mobilized dependent on rock damage or plastic parameters of geo-materials. The most widely accepted plastic parameter is the plastic shear strain p , which can be obtained as the difference between the maximum and minimum principal plastic strains ( p and p , respectively) [27,30,33,54]: There are usually two methods to obtain the plastic strain values. One method is to differentiate the recoverable and irrecoverable strain by taking cyclic loading-unloading experiments. The plastic strain can be obtained from the irrecoverable strain in each cycle of the tests directly, however, it is quite complicated to control this type of experiment, and the data is limited by the numbers of cycles [27,33,54,55]. Therefore, another method is developed based on the assumption that the unloading curve in each cycle has the same modulus as the initial deformation modulus. In this way, a series of

Analytical Method
According to the previous studies [26,28,31,33,53], the Mohr-Coulomb strength parameters (cohesion c and inner frictional angle ϕ) of rock should be mobilized dependent on rock damage or plastic parameters of geo-materials. The most widely accepted plastic parameter is the plastic shear strain γ p , which can be obtained as the difference between the maximum and minimum principal plastic strains (ε p 1 and ε p 3 , respectively) [27,30,33,54]: There are usually two methods to obtain the plastic strain values. One method is to differentiate the recoverable and irrecoverable strain by taking cyclic loading-unloading experiments. The plastic strain can be obtained from the irrecoverable strain in each cycle of the tests directly, however, it is quite complicated to control this type of experiment, and the data is limited by the numbers of cycles [27,33,54,55]. Therefore, another method is developed based on the assumption that the unloading curve in each cycle has the same modulus as the initial deformation modulus. In this way, a series of plastic strains can be obtained with a series of assumed loading-unloading cycles by just carrying out conventional uniaxial and tri-axial compression experiments. This method has been widely accepted and used in many studies [27,33,54] and is also employed here in this research. Figure 4 gives a sketch to illustrate this method for determining the plastic axial and lateral strains, as well as the corresponding axial stress values. A series of lines parallel with the tangent lines at the linear elastic stage of the σ 1 -ε 1 curves are drawn to determine the plastic axial strain ε p 1,i and plastic lateral strain ε p 3,i , respectively. The symbol i here is a series of positive integers, showing that a series of plastic strain values can be collected with this method. It should be noted that there is a gap |OA| owing to the crack closure stage of σ 1 -1 curves, so this gap should be removed for determining the plastic axial strain ε p 1,i : ε Energies 2019, 12, x FOR PEER REVIEW 7 of 15 plastic strains can be obtained with a series of assumed loading-unloading cycles by just carrying out conventional uniaxial and tri-axial compression experiments. This method has been widely accepted and used in many studies [27,33,54] and is also employed here in this research. Figure 4 gives a sketch to illustrate this method for determining the plastic axial and lateral strains, as well as the corresponding axial stress values. A series of lines parallel with the tangent lines at the linear elastic stage of the σ1-ε1 curves are drawn to determine the plastic axial strain , p and plastic lateral strain , p , respectively. The symbol i here is a series of positive integers, showing that a series of plastic strain values can be collected with this method. It should be noted that there is a gap | | owing to the crack closure stage of σ1-1 curves, so this gap should be removed for determining the plastic axial strain , p : Then the plastic shear strain can be obtained as: Figure 4. Sketch of the method for determining plastic strains and the corresponding stress values.
For each plastic shear strain, the corresponding maximum principal stress σ1 is collected under different confining pressures σ3. Thereafter, cohesion c and internal friction angle can be calculated by drawing Mohr circles or by linear fitting of σ1 and σ3 with the following equation: 2 cos 1 sin 1 sin 1 sin The values of c and can then be plotted with the increasing plastic shear strain. As the temperature induced by the high-level radioactive waste may affect the mechanical behavior of the host rock, the evolutionary characteristics of c and are also studied for NRG01 granite under different heat treatment (20 °C-200 °C).

Data Analyses and Discussion
Based on the stress-strain curves of NRG01 granite under different heat treatment and confinements presented in Figure 3, as well as the methodology demonstrated in Sub-section 4.1.1, a series of axial stress at different plastic shear strains can be plotted in Figure 5. For each plastic shear strain, a set of axial stress values under various confining pressures can be obtained to calculate the cohesion and friction angle. These values are shown in Figure 5. Then the plastic shear strain can be obtained as: For each plastic shear strain, the corresponding maximum principal stress σ 1 is collected under different confining pressures σ 3 . Thereafter, cohesion c and internal friction angle ϕ can be calculated by drawing Mohr circles or by linear fitting of σ 1 and σ 3 with the following equation: The values of c and ϕ can then be plotted with the increasing plastic shear strain. As the temperature induced by the high-level radioactive waste may affect the mechanical behavior of the host rock, the evolutionary characteristics of c and ϕ are also studied for NRG01 granite under different heat treatment (20 • C-200 • C).

Data Analyses and Discussion
Based on the stress-strain curves of NRG01 granite under different heat treatment and confinements presented in Figure 3, as well as the methodology demonstrated in Section 4.1.1, a series of axial stress at different plastic shear strains can be plotted in Figure 5. For each plastic shear strain, a set of axial  The results obtained in Figure 6 show that NRG01 granite samples have the generally cohesion weakening and friction strengthening (CWFS) behaviors for various treatment temperatures (room temperature to 200 °C). According to the published references [26,29,31,36], the cohesion component should be weakened to the residual value before the full mobilization of friction angle, however, it is not true for NRG01 granite treated by different temperatures. It is shown that cohesion is weakened in a gradual manner with increasing plastic shear strain, nevertheless, the friction angle increases to the peak value more immediately. As an actual fact, the test results similar to this study can also be found in references [30,31,33]. This difference has also been discussed in [29,36], and it is believed that the plastic strain limit at which the cohesion reaches the residual value or the friction angle is fully mobilized is dependent on many factors such as the rock type, grain size, heterogeneity, as well as the hoop effect owing to the cylindrical shape of the specimens, etc. More systematic studies should be carried out to learn more clearly about the exact influencing factors and the mechanism. The results obtained in Figure 6 show that NRG01 granite samples have the generally cohesion weakening and friction strengthening (CWFS) behaviors for various treatment temperatures (room temperature to 200 • C). According to the published references [26,29,31,36], the cohesion component should be weakened to the residual value before the full mobilization of friction angle, however, it is not true for NRG01 granite treated by different temperatures. It is shown that cohesion is weakened in a gradual manner with increasing plastic shear strain, nevertheless, the friction angle increases to the peak value more immediately. As an actual fact, the test results similar to this study can also be found in references [30,31,33]. This difference has also been discussed in [29,36], and it is believed that the plastic strain limit at which the cohesion reaches the residual value or the friction angle is fully mobilized is dependent on many factors such as the rock type, grain size, heterogeneity, as well as the hoop effect owing to the cylindrical shape of the specimens, etc. More systematic studies should be carried out to learn more clearly about the exact influencing factors and the mechanism. Based on the characteristics of the mobilized cohesion and friction angle presented in Figure 6, the generally used linear CWFS model [30,[37][38][39] may not be suitable for NRG01 granite. For a better description of the rock behaviors, a non-linear model should be used. With the fitting Equations (5) and (6) proposed in reference [31], the mobilized cohesion and friction angle values can be well fitted as shown in Figure 6. The fitted coefficients are listed in Table 3.
where, i and r are the initial and residual values of cohesion, respectively. ,r p is the plastic shear strain when the cohesion is close to the residual value.
where, i is the initial value of friction angle, while max is the maximum value, and r is the residual value. ,max p is the plastic shear strain when the friction angle is close to its peak value, while ,r p is the plastic shear strain when the friction angle is close to its residual value. According to Figure 6, it can also be observed that the different treated temperatures may lead to a few different evolutionary behaviors of cohesion and friction angle for NRG01 granite samples.   Based on the characteristics of the mobilized cohesion and friction angle presented in Figure 6, the generally used linear CWFS model [30,[37][38][39] may not be suitable for NRG01 granite. For a better description of the rock behaviors, a non-linear model should be used. With the fitting Equations (5) and (6) proposed in reference [31], the mobilized cohesion and friction angle values can be well fitted as shown in Figure 6. The fitted coefficients are listed in Table 3.
where, c i and c r are the initial and residual values of cohesion, respectively. γ p c,r is the plastic shear strain when the cohesion is close to the residual value.
where, ϕ i is the initial value of friction angle, while ϕ max is the maximum value, and ϕ r is the residual value. γ p ϕ,max is the plastic shear strain when the friction angle is close to its peak value, while γ p ϕ,r is the plastic shear strain when the friction angle is close to its residual value. According to Figure 6, it can also be observed that the different treated temperatures may lead to a few different evolutionary behaviors of cohesion and friction angle for NRG01 granite samples. The more obvious influences are shown for the behaviors of cohesion component, i.e., the sample under room temperature (T = 20 • C) shows an apparently higher initial cohesion value (c = 31.99 MPa), and a more obvious decrease with the increasing plastic shear strain, compared with the samples treated by higher temperatures (c = 24.67 MPa and 25.88 MPa for T = 100 • C and 200 • C, respectively). This should be explained by the more thermally induced cracks in the heat treated samples, which decreased the initial cohesive strength of the specimens. According to Figure 6b, the thermally treated samples by higher temperatures present higher initial friction angles (ϕ = 38.87 • for T = 200 • C, compared with ϕ = 32.54 • and 33.74 • for T = 20 • C and 100 • C, respectively), which should also be resulted from the more thermally induced crack surfaces treated by higher temperatures. Based on these observations, a general trend can be concluded that higher treatment temperatures may lead to relatively lower initial values of cohesion (c) and higher initial values of friction angles (ϕ) for NRG01 granite samples. This phenomenon should be owing to the different amounts of thermally induced cracks inside the rock under the effects of different temperatures.
This section demonstrates the characteristics of mobilized cohesion and friction angle during the failure process of NRG01 granite samples treated by different temperatures. The non-simultaneous mobilization should be considered in the constitutive models when analyzing the stability of the host rock for site selection or design of a HLW disposal repository.

Mobilized Hoek-Brown Strength Parameters During Failure of NRG01 Granite
As discussed above, the non-simultaneous mobilization of strength parameters mainly focused on the linear Mohr-Coulomb criterion. For the widely used non-linear Hoek-Brown criterion, only simultaneous mobilization of strength parameters (m and s) can be found to be considered in the published studies [23][24][25][26]. It is quite necessary to research the mobilization behaviors of Hoek-Brown strength parameters during the failure process of rock. This section will present such a study based on the laboratory experiments on NRG01 granite samples treated by different temperatures.

Analytical Method
The similar method as illustrated in Figure 4 is used here for determining a series of plastic strain and axial stress values. So the same results presented in Figure 5 can be used in this part of analyses. For each certain plastic shear strain value, the set of stress values under different confining pressures are used to fit the Hoek-Brown criterion [26,42]: where, σ 1 and σ 3 are the maximum and minimum principal stresses, respectively; σ ci is the uniaxial compression strength of the intact rock; m b , s and a are the constants for the damaged rock specimens. It should be noted that the parameter s is related to the fracturing degree of the samples. When σ 3 = 0, the uniaxial compression strength of the damaged samples can be obtained as: As the constant a is defined as: where, GSI is the Geological Strength Index indicating the rock quality with the number ranging from 0 (totally fractured) to 100 (intact) [41,42]. It can be seen that for 25 < GSI < 100, the value of a is very close to 0.5. Therefore, a is reasonable enough to be set as a constant number 0.5 for simplicity in this study. Thereafter, for the damaged samples at each plastic shear strain, s can be identified based on the corresponding uniaxial compression strength values presented in Figure 5 by the following equation: With the determined s values, a series of m b values can be obtained by fitting the data in Figure 5 with Equation (7). As a result, the mobilized m b and s values varying with plastic shear strain for NRG01 granite samples treated by different temperatures are presented in Figure 7a,b, respectively.

Data Analyses and Discussion
According to Figure 7, for NRG01 granite samples treated by different temperatures (T = 20 °C, 100 °C and 200 °C), several features can be observed as follows: (1) Similar to the mobilization of cohesion and friction angle, the mobilization of Hoek-Brown strength parameters ( b and ) is also non-simultaneous during the failure process of NRG01 granite treated by different temperatures no higher than 200 °C; (2) With increasing plastic shear strain, b increases significantly to a maximum value and then decreases until a residual value; (3) decreases gradually with the increasing plastic shear strain. This is related to the damage and fracturing process during the tests.
The mobilization of b can be fitted with Equation (11) where, bi and br are the initial and residual value of b , respectively; bmax is the maximum value of b ; b,max p is the plastic shear strain when b is close to its peak value, while b,r p is the plastic shear strain when b is close to its residual value.
The mobilization of can be fitted with Equation (12): where, i and r are the initial and residual value of , respectively. 0 p is the transitional plastic shear strain when the value turns to decrease in a gradual manner. is a constant determining the shape of the curve.

Data Analyses and Discussion
According to Figure 7, for NRG01 granite samples treated by different temperatures (T = 20 • C, 100 • C and 200 • C), several features can be observed as follows: (1) Similar to the mobilization of cohesion and friction angle, the mobilization of Hoek-Brown strength parameters (m b and s) is also non-simultaneous during the failure process of NRG01 granite treated by different temperatures no higher than 200 • C; (2) With increasing plastic shear strain, m b increases significantly to a maximum value and then decreases until a residual value; (3) s decreases gradually with the increasing plastic shear strain. This is related to the damage and fracturing process during the tests.
The mobilization of m b can be fitted with Equation (11): where, m bi and m br are the initial and residual value of m b , respectively; m bmax is the maximum value of m b ; γ p mb,max is the plastic shear strain when m b is close to its peak value, while γ p mb,r is the plastic shear strain when m b is close to its residual value.
The mobilization of s can be fitted with Equation (12): where, s i and s r are the initial and residual value of s, respectively. γ p 0 is the transitional plastic shear strain when the s value turns to decrease in a gradual manner. n is a constant determining the shape of the curve.
The fitted curves are presented in Figure 7, and the fitted coefficients are shown in Table 4. According to Figure 7a, it can be observed that the mobilization of m b value has very similar characteristics with the increasing plastic shear strain, for the NRG01 granite samples treated by different temperatures (T = 20 • C, 100 • C and 200 • C). Based on a more detailed observation on Figure 7a, we can find that for the NRG01 granite samples treated by the temperature T = 200 • C, the initial value of m b = 7.24. This is apparently higher than the value for the cases of lower temperatures (m b = 3.54 and 2.59, for T = 20 • C and 100 • C, respectively). It is always believed that m b value is more related to the frictional strength in Mohr-Coulomb criterion [29,56]. For the case of T = 200 • C, there should be more crack surfaces induced by the heat in the granite samples, and this should be the reason why the initial values of friction angle and m b are both higher. According to Figure 7b, the mobilization of s value are also quite similar with the increasing plastic shear strain, for the NRG01 granite samples treated by different temperatures (T = 20 • C, 100 • C and 200 • C). There are not enough evidence to prove how heat treatment influence the mobilization of s value based on this study. More systematic experimental studies should be carried out in order to make clear the characteristics of mobilized s during the failure process of granite treated by different temperatures.

Conclusions
NRG01 granite samples cored from Alxa candidate area for HLW disposal were treated by different temperatures, and then were used to carry out a series of uniaxial and tri-axial compression experiments under different confining pressures. Complete axial stress-axial strain curves and axial stress-lateral strain curves were recorded. These data were collected to study the mobilization of both Mohr-Coulomb and Hoek-Brown strength parameters during the damage and failure of NRG01 granites samples considering the effect of heat induced by the nuclear waste. According to the analyses in this study, several conclusions can be drawn as follows: (1) Cohesion weakening and friction angle strengthening occurs during the damage and failure process of NRG01 granite samples treated by different temperatures. However, compared with the findings in the previous studies, cohesion decreases in a more gradual manner for NRG01 granite samples, and the friction angle increases immediately to its maximum value before the cohesion approaching to the residual value. This may be owing to the grain size, heterogeneity, or even the hoop effect induced by the cylindrical shape of the samples. More systematic studies are required to make clear the exact influencing factors, as well as the mechanism. (2) The temperatures of no higher than 200 • C do not have significant influence on the characteristics of mobilized cohesion or friction angle during the damage and failure process of NRG01 granite samples. However, the samples under room temperature (20 • C) have higher initial cohesion than the samples treated by higher temperatures (T=100 • C and 200 • C). In addition, the samples treated by temperature of 200 • C have higher friction angle than the samples treated by lower temperatures. This should be caused by the cracks induced by the heat treatment.
(3) The Hoek-Brown strength parameters m b and s are also observed to show non-simultaneous mobilization behaviors during the failure process of NRG01 granite samples treated by different temperatures. It is found that m b increases significantly to a maximum value and then decreases until a residual value, and s decreases gradually with the increasing plastic shear strain. The general characteristics of the mobilized m b and s are similar for NRG01 granite samples treated by different temperatures, and the fitted equations for modelling the mobilization of both parameters are proposed. The samples treated by temperature of 200 • C have higher initial m b value, this should also be caused by the cracks induced by the heat treatment.
These findings on the mobilization of strength parameters provide a better understanding on the strength properties of NRG01 granite samples, and can be used for building a plastic constitutive model in the next step. This study should also be helpful for guiding the selection and design of HLW disposal repository in Alxa area in China. This study put forward the research on non-simultaneous mobilization of strength parameters to Hoek-Brown strength criterion, and more experimental studies are required to consolidate the results. The methods used in this paper can also be used for this kind of analyses in the other candidate areas for HLW disposal.