# A New Strain-Softening Constitutive Model for Circular Opening considering Plastic Bearing Behavior and Its Engineering Application

Authors: Liang Chen; Xianbiao Mao; Ming Li; Wei Zhou
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

Geomaterials generally show strain-softening characteristics after peak-load. Based on the triaxial test for sandy mudstone, a simple elastopeak plastic-strain-softening-damage model (EPSDM) was proposed. Compared with the traditional strain-softening model, EPSDM shows obvious plastic bearing characteristics before strain softening. Then, the closed-formed solution of circular opening was deduced based on the newly proposed model. A plastic shear strain increment was introduced as the extension constraint condition of peak plastic zone. The solution correctness of EPSDM was also verified by comparing with other research results. In addition, the solution based on EPSDM could degenerate for a series of results obtained by elastobrittle plastic model (EBM), elasto-strain-softening model (ESM), and elasto-perfectly plastic model (EPM) under certain conditions. Hence, it could be regarded as a unified solution. Finally, the research results denoted that when the inner pressure was fully released, the maximum postpeak failure radii and surface displacement of surrounding rock indicated the characteristics of EBM>ESM>EPSDM>EPM. Therefore, the plastic bearing behavior could effectively decrease the postpeak failure zone radii and surface displacement. The dilation coefficient noticeably influenced postpeak failure range and surface displacement, particularly the damage zone radii and tunnel wall convergence. The research results can provide very important theoretical bases for evaluating the tunnel stability and support design reliability for underground engineering.


---

## Body

## 1. Introduction

The stresses and deformation of surrounding rock are an important basis for evaluating the tunnel stability and support design reliability for underground engineering. The ground reaction curve (GRC) of circular tunnel is generally required to be given to optimize support design. Originally, it was researched by regarding the rock mass as the elastic material. However, Placidi L and Altenbach H et al. pointed out that the inelastic phenomena like damage, brittle fracture, and plasticity are also extremely important for the mechanical behavior of materials [1–3]. Therefore, many studies began to predict the evolution law of surrounding rock states by using different constitutive models, such as elasto-perfectly plastic (EPM) model [4–8], elastobrittle plastic (EBM) model [9–16], elastopeak plastic-brittle plastic (EPBM) model [17, 18], and elasto-strain-softening (ESM) model [19–28] (see Figure 1). Nevertheless, each constitutive model has its own application scope. Based on the geological quality, the application scope of the above models was discussed by Hoek and Brown [29]. The results show that EBM applies only to the poor quality rock mass, while EPM is suitable for high quality rock mass and ESM for the average quality rock mass. Moreover, the EPBM is suitable for brittle rock mass with plastic bearing behavior [18]. In fact, a number of strain-softening rock masses also obviously presented plastic bearing behavior. Since the early 1970s till now, the study on strain-softening model can be divided into two categories. As shown in Figure 1(c), one is that the strain-softening zone is simplified into a continuous slash of which the slope can be obtained by the stress-strain curve [25, 26]. The other is that the strain-softening zone is considered to be composed of a series of straight lines which show the instantaneous drop characteristics (see Figure 1(d)) [19–24, 27, 28]. In addition, for the first model, when the slope of slash is equal to zero, the strain-softening models will be converted into elasto-perfectly plastic model; when the slope of slash is up to infinity, the strain-softening models will become elastobrittle plastic model. Geomaterials are easily affected by internal fissures, joints, components, and external environment. Therefore, its constitutive relation shows obvious diversity and complexity.Figure 1
Previous research model.
(a) (b) (c) (d)As shown in Figure2, the No. 4 coal mine of Saier Group in China is located between Altai mountain and Junger basin. The sandy mudstone is taken from the haulage roadway of the No. 4 coal mine, and the triaxial test result was also given. From the total stress-strain curve, it can be seen that the sandy mudstone firstly shows obvious peak plastic state (“AB”) when the initial yield occurs; then, as the strain continues to increase, the strain-softening stage will begin to develop (“BC”) with the rock mass parameter gradually decreasing. Until a certain residual value is reached, the surrounding rock starts to enter the residual flow stage. And, so, the stress-strain curve approximately experienced four phases during the triaxial test. However, above four models (EPM, EBM, and ESM) cannot be consistent with this experiment result. In other words, if the test results are simplified to the above four models, the calculation result deviation may be larger compared with the actual engineering. Therefore, the elastopeak plastic-strain-softening-damage model (EPSDM) was proposed and applied to the underground engineering. The EPSDM, including all the features of the EPM, EBM, and ESM, represents the influence of the plastic bearing behavior and strain softening on surrounding rock state and can be also transformed into the above three models under certain conditions. Hence, it can be regarded as a unified constitutive model.Figure 2
Sandy mudstone triaxial test.This current study mainly deals with a new constitutive model of strain-softening rock mass with consideration given to the plastic bearing behavior and strain-softening characteristics. Moreover, a uniform analytical solution for circular tunnel based on the EPSDM is also derived.


## 2. Problem Descriptions

### 2.1. Constitutive Model

As shown in Figure2, the completed stress-strain curve of sandy mudstone was obtained in MTS816 experimental system under different confining pressures. It can be seen that the rock mass firstly experiences elastic stage (“OA”), and then the plastic stage begins to occur when the principal stress satisfies the initial yield function. This stage represents the plastic bearing capacity of sandy mudstone. However, its extension restricted by some factors is not infinite. This paper assumes that once the plastic shear strain increment of “AB” segment reaches a certain value, the strain-softening phase (“BC”) begins to happen. As the deformation continues to increase, the parameters of rock mass will gradually decrease. Until a residual parameter is reached, the surrounding rock starts to enter the damage stage.In other words, the mechanical behavior for sandy mudstone can be represented by four phases as follows: elastic stage, peak plastic stage, strain-softening stage, and damage stage. It is obvious that the proposed EPSDM combines the plastic bearing characteristics of EPM, strain-softening feature of ESM, and residual flow characteristics of EBM as well as ESM. Hence, it can be regarded as a unified constitutive model.


### 2.2. Establishment of Mechanical Model

Figure3 shows a circular opening being excavated in an infinite, homogeneous, isotropic EPSDM rock mass subjected to a hydrostatic pressure (p0) at the infinity boundary and inner pressure (pin) at inner radii r=R0. As pin gradually reduced, the displacement occurs and the peak plastic zone firstly develops around the circular opening when pin is less than the initial yield stress. Then, the change of surrounding rock state is consistent with the material postpeak behavior of the elastopeak plastic-strain-softening-damage model. When the inner pressure is completely released, there are four different zones around the circular opening: elastic zone (“e”), peak plastic zone (“p”), strain-softening zone (“s”), and damage zone (“c”). The radii of the peak plastic zone, strain-softening zone, and damage zone are denoted as R3, R2, and R1 respectively. The mechanical model should meet the following assumption conditions:(i)
The tunnel is excavated in an infinite geological body, so the problem can be regarded as plane strain problem to be solved.(ii)
The maximum and minimum principal strains are only composed of plastic strain at the postpeak stage, respectively.Figure 3
Mechanical model of circular opening considering EPSDM.For axisymmetric plane strain problems, whenpin<p0, the tangential and radial stresses, σθ and σr, are, respectively, the maximum and minimum principal stresses. εθ and εr are the maximum and minimum principal strain, respectively. Accordingly, the linear Mohr-Coulomb (M-C) criterion as the rock mass yield condition can be expressed as(1)σθ=ξσr+σcPeakfor peak plastic zone(2)σθ=ξσr+σc¯for strain softening zone(3)σθ=ξσr+σcRfor damage zonewhere ξ=1+sin⁡φ/1-sin⁡φ; φ is the internal friction angle; σcPeak and σcR are the initial and residual uniaxial compressive strength of rock mass; σcPeak=2ccos⁡φ/1-sin⁡φ; and σcR=2crcos⁡φ/1-sin⁡φ. c and cr represent initial and residual cohesion; σc¯ represents the softening strength in strain-softening zone. It should be noted that this paper assumes that the decrease of stress-strain curve is only related to the cohesion attenuation, which is independent of the internal friction angle [25].


## 2.1. Constitutive Model

As shown in Figure2, the completed stress-strain curve of sandy mudstone was obtained in MTS816 experimental system under different confining pressures. It can be seen that the rock mass firstly experiences elastic stage (“OA”), and then the plastic stage begins to occur when the principal stress satisfies the initial yield function. This stage represents the plastic bearing capacity of sandy mudstone. However, its extension restricted by some factors is not infinite. This paper assumes that once the plastic shear strain increment of “AB” segment reaches a certain value, the strain-softening phase (“BC”) begins to happen. As the deformation continues to increase, the parameters of rock mass will gradually decrease. Until a residual parameter is reached, the surrounding rock starts to enter the damage stage.In other words, the mechanical behavior for sandy mudstone can be represented by four phases as follows: elastic stage, peak plastic stage, strain-softening stage, and damage stage. It is obvious that the proposed EPSDM combines the plastic bearing characteristics of EPM, strain-softening feature of ESM, and residual flow characteristics of EBM as well as ESM. Hence, it can be regarded as a unified constitutive model.


## 2.2. Establishment of Mechanical Model

Figure3 shows a circular opening being excavated in an infinite, homogeneous, isotropic EPSDM rock mass subjected to a hydrostatic pressure (p0) at the infinity boundary and inner pressure (pin) at inner radii r=R0. As pin gradually reduced, the displacement occurs and the peak plastic zone firstly develops around the circular opening when pin is less than the initial yield stress. Then, the change of surrounding rock state is consistent with the material postpeak behavior of the elastopeak plastic-strain-softening-damage model. When the inner pressure is completely released, there are four different zones around the circular opening: elastic zone (“e”), peak plastic zone (“p”), strain-softening zone (“s”), and damage zone (“c”). The radii of the peak plastic zone, strain-softening zone, and damage zone are denoted as R3, R2, and R1 respectively. The mechanical model should meet the following assumption conditions:(i)
The tunnel is excavated in an infinite geological body, so the problem can be regarded as plane strain problem to be solved.(ii)
The maximum and minimum principal strains are only composed of plastic strain at the postpeak stage, respectively.Figure 3
Mechanical model of circular opening considering EPSDM.For axisymmetric plane strain problems, whenpin<p0, the tangential and radial stresses, σθ and σr, are, respectively, the maximum and minimum principal stresses. εθ and εr are the maximum and minimum principal strain, respectively. Accordingly, the linear Mohr-Coulomb (M-C) criterion as the rock mass yield condition can be expressed as(1)σθ=ξσr+σcPeakfor peak plastic zone(2)σθ=ξσr+σc¯for strain softening zone(3)σθ=ξσr+σcRfor damage zonewhere ξ=1+sin⁡φ/1-sin⁡φ; φ is the internal friction angle; σcPeak and σcR are the initial and residual uniaxial compressive strength of rock mass; σcPeak=2ccos⁡φ/1-sin⁡φ; and σcR=2crcos⁡φ/1-sin⁡φ. c and cr represent initial and residual cohesion; σc¯ represents the softening strength in strain-softening zone. It should be noted that this paper assumes that the decrease of stress-strain curve is only related to the cohesion attenuation, which is independent of the internal friction angle [25].


## 3. Analytical Solution of EPSDM

### 3.1. Basic Equations

The equilibrium differential equation for axisymmetric plane strain problem in the “i” zone can be expressed as (ignoring the body force of rock mass) [4, 9, 12, 16, 22, 25](4)dσridr+σri-σθir=0where σri and σθi are radial stress and tangential stress in the “i” zone, respectively. The subscript symbol “i” represents different zone of surrounding rock, which can be replaced by the symbols “e,” “p,” “s,”and “c.”Based on the small deformation assumption, the geometric equation can be denoted as [9, 17, 24](5)εri=duridr,εθi=urirwhere εri and εθi are, respectively, radial and circumferential strain in different zones and uri represents the radial displacement.Supposing that the volume of the rock mass is changing in the postpeak failure zone, the plastic-strain relationships can be established by adopting a nonassociated linear flow rule and small strain theory as follows [11, 14, 22]:(6)εri+βiεθi=0where βi=1+sin⁡ψi/1-sin⁡ψi which is the dilation coefficient in the “i” zone. ψi is the dilation angle.


### 3.2. Elastic Zone

According to the elasticity theory, the stress functionΦ can be expressed as [30](7)Φ=Mln⁡r+Nr2where M and N are the constants. The radial and tangential stresses are then given by(8)σre=1rdΦdr=Mr2+2N(9)σθe=d2Φdr2=-Mr2+2N.Substituting the boundary conditionsσre=p0 at r→∞ and σre=σR3 at r=R3 into (8), where σR3 is the radial stress at the elastic-peak plastic zone interface, (8) and (9) can be rewritten as(10)σre=p0+σR3-p0R32r2(11)σθe=p0-σR3-σ0R32r2.Whenr=R3, (10) and (11) should satisfy (1). Thus, σR3 can be deduced as follows:(12)σR3=2p0-σcPeak1+ξ.The displacement and strains can be easily obtained as follows [25]:(13)ure=A0R32r,εre=-A0R32r2,εθe=A0R32r2where A0=(1+ν)(σ0-σR3)/E. ν and E are Poisson’s ratio and Young’s modulus. When r=R3, the radial displacement uR3 at the elastic-peak plastic zone interface is(14)uR3=A0R3.


### 3.3. Peak Plastic Zone

Substituting (1) into (4), the equilibrium differential equation in peak plastic zone can be rewritten as follows [25]:(15)dσrpdr+1-ξσrp-σcPeakr=0.Solving (15), the radial stress in peak plastic zone can be derived by considering the boundary conditions σrp=σR3 at r=R3.(16)σrp=σR3-σcPeak1-ξrR3ξ-1+σcPeak1-ξ.And then, introducing (16) into (1), the circumferential stress will be also be deduced as(17)σθp=ξσR3-σcPeak1-ξrR3ξ-1+σcPeak1-ξ.The radial contact stressesσR2 at the peak plastic-strain-softening zone interface can be easily obtained by substituting r=R2 into (16):(18)σR2=σR3-σcPeak1-ξR2R3ξ-1+σcPeak1-ξ.Then, the displacement differential equation can be also given by integrating (5) and (6) [15, 24].(19)durpdr+βpurpr=0.Solving (19), the radial displacement will be deduced in peak plastic zone by considering the boundary condition urp=uR3 at r=R3.(20)urp=A0R31+βpr-βp.According to (5), the strains in peak plastic zone can be also obtained as follows:(21)εrp=-A0βpR31+βpr-βp-1(22)εθp=A0R31+βpr-βp-1.According to (20), the radial displacement uR2 at the peak plastic-strain-softening zone interface can be achieved as follows:(23)uR2=A0R31+βpR2-βp.In addition, the maximum principal strainε1sPeak at the peak plastic-strain-softening zone interface can be also deduced by substituting r=R2 into (20) and (22).(24)ε1sPeak=A0R31+βpR2-βp-1.


### 3.4. Strain-Softening Zone

Based on the calculation method of peak plastic zone, the radial displacement in strain-softening zoneurs is easily acquired as follows [24]:(25)urs=A0R31+βpR2βs-βpr-βs.Then, the radial and tangential strain in strain-softening zone are also deduced by integrating (5) and (25).(26)εrs=-A0βsR31+βpR2βs-βpr-βs-1(27)εθs=A0R31+βpR2βs-βpr-βs-1.Meanwhile, according to (25), the radial displacement uR1 at the strain-softening-damage zone interface can be obtained as well.(28)uR1=A0R31+βpR2βs-βpR1-βs.In strain-softening zone, it is assumed that the strength attenuation of rock mass is only related to cohesion (c). Then, the compressive strength at any point can be expressed as(29)σc¯=σcPeak-σcPeak-σcRεcR-ε1sPeakεθs-ε1sPeak=σcPeak-λεθs-ε1sPeakwhere εcR is the maximum principal strain at the point “C” in Figure 3. λ may be called strain-softening modulus and can be determined by the slope of “BC” segment in Figure 3. Introducing (24) and (27) into (29), it can be rewritten as(30)σc¯=σcPeak-α1+νp0-σR3R31+βpR2-βp-1R2r1+βs-1where α=λ/E, which can be defined as a strain-softening coefficient.Substituting (2) and (30) into (4), the equilibrium differential equation in strain-softening zone can be expressed as(31)dσrsdr+1-ξσrsr-σcPeak-α1+νp0-σR3R31+βpR2-βp-1R2/r1+βs-1r=0.Combined with the boundary conditionσrs=σR2 at r=R2, the radial stress in strain-softening zone can be derived by solving (31).(32)σrs=σR3-σcPeak1-ξR3r1-ξ+σcPeak1-ξ+α1+νp0-σR3βs+ξR3R21+βpR2r1+βs-R2r1-ξ+α1+νp0-σR31-ξR3R21+βp1-R2r1-ξ.The circumferential stress is also obtained by substituting (30) and (32) into (2).(33)σθs=ξσR3-σcPeak1-ξR3r1-ξ+σcPeak1-ξ+α1+νp0-σR3βs+ξR3R21+βpR2r1+βs-R2r1-ξ+α1+νp0-σR31-ξR3R21+βp1-R2r1-ξ+σcPeak-α1+νp0-σR3R3R21+βpR2r1+βs-1.


### 3.5. Damage Zone

Substituting (3) into (4), the equilibrium differential equation in damage zone can be rewritten as follows:(34)dσrcdr+1-ξσrc-σcRr=0.Solving (34) and considering the boundary conditions σrc=pin at r=R0, the radial stress in damage zone can be derived as(35)σrc=pin-σcR1-ξrR0ξ-1+σcR1-ξ.Then, introducing (35) into (3), the circumferential stress in damage zone will be also obtained as(36)σθc=ξpin-σcR1-ξrR0ξ-1+σcR1-ξ.As shown in (25), the radial displacement of damage zone can be deduced by considering the boundary condition (urc)r=R1=uR1.(37)urc=A0R31+βpR2βs-βpR1βc-βsr-βc.The radial and circumferential strains in damage zone can be also obtained by substituting (37) into (5).(38)εrc=-βcA0R31+βpR2βs-βpR1βc-βsr-βc-1(39)εθc=A0R31+βpR2βs-βpR1βc-βsr-βc-1.According to (37), when r=R0, the surface displacement of surrounding rock u0 is as follows:(40)u0=A0R31+βpR2βs-βpR1βc-βsR0-βc.


### 3.6. Radius of Postpeak Failure Zone

In order to obtain the closed-form solution of the stresses and deformation of EPSDM, the postpeak failure zone radius (R1,R2,R3) should be firstly determined. According to the stress contact condition (σrs)r=R1=(σrc)r=R1 at r=R1, one of the relationship among R1, R2, and R3 can be established by combining with (32) and (35).(41)σR3-σcPeak1-ξR3R11-ξ+σcPeak1-ξ+α1+νp0-σR3βs+ξR3R21+βpR2R11+βs-R2R11-ξ+α1+νp0-σR31-ξR3R21+βp1-R2R11-ξ-pin-σcR1-ξR1R0ξ-1-σcR1-ξ=0.According to (30), when r=R1, the residual strength can be expressed as follows:(42)σcR=σcPeak-α1+νp0-σR3R31+βpR2-βp-1R2R11+βs-1.In addition to (41) and (42), a new condition is still needed to make the EPSDM form a closed solution. Zhang and Jiang et al. assumed that the surrounding rock begins to enter the next stage when the plastic shear strain or the equivalent plastic shear strain in peak plastic zone reaches a critical value [17, 18]. This paper assumes that when the plastic shear stress increment Δγ∗ between the yield points “A” and “B” (see Figure 3) satisfies certain conditions, the surrounding rock will enter the strain-softening state. Thus, the other relationship among R3, R2, and R1 can be expressed as(43)Δγ∗=γB∗-γA∗=ε1sPeak-ε3sPeak-ε1pPeak-ε3pPeak=εθp-εrpr=R2-εθp-εrpr=R3where γB∗ and γA∗ are the plastic shear strain at points “B” and “A,” respectively. ε1pPeak is the maximum principal strain at the point “A.” ε3sPeak and ε3pPeak are the minimum principal strain at the points “B” and “A.” As shown in Figure 3, the maximum and minimum principal strain of rock mass can be easily obtained in the triaxial test. Substituting (21) and (22) into (43), it can be rewritten as(44)R3=Δγ∗A01+βp+11/1+βpR2=TR2.Introducing (44) into (42), the following expression can be obtained as follows:(45)R2=σcPeak-σcRα1+νp0-σR3T1+βp+11/1+βsR1=tR1.Then, the damage zone radius can be deduced by substituting (44) and (45) into (41).(46)R1=R0σR3-σcPeak1-ξT1-ξt1-ξ-α1+νp0-σR3βs+ξ+α1+νp0-σR31-ξT1+βpt1-ξ+α1+νp0-σR3βs+ξt1+βs+α1+νp0-σR31-ξT1+βp+σcPeak-σcR1-ξ×pin-σcR1-ξ-11/ξ-1.Subsequently, the radiusesR3 and R2 can be also derived by substituting (46) into (45) and (44), respectively.


### 3.7. Critical Inner Pressure under EPSDM

According to (44),(45), and (46), when the relationship between postpeak failure zone radius (R1,R2,R3) and inner radius (R0) satisfies certain condition, the rock mass will be in the critical state.(i)
WhenR0=R1=R2=R3, the rock mass is in the critical state that the peak plastic is just not arisen yet. Then, the critical inner pressure pcrmaxE-P at this state will be deduced according to (12)(47)pcrmaxE-P=σR3=2p0-σcPeak1+ξ.(ii)
WhenR0=R1=R2<R3, the rock mass is in the critical state that the peak plastic zone has reached the maximum and the stain softening zone has not yet appeared. Consequently, the critical radii of peak plastic zone can be obtained by (44).(48)R3=Δγ∗A01+βp+11/1+βpR0.Then, introducing (48) into (16), the critical inner pressure pcrmaxP-S at this state will be deduced as follows:(49)pcrmaxP-S=σR3-σcPeak1-ξΔγ∗A01+βp+1-ξ-1/1+βp+σcPeak1-ξ.(iii)
WhenR0=R1<R2<R3, the rock mass is in the critical state that the damage zone begins to form. According to (45), the strain-softening zone radius can be expressed as(50)R2=σcPeak-σcRα1+νp0-σR3T1+βp+11/1+βsR1=tR0.The peak plastic zone radius can be calculated by (44) at this state. Then, when r=R0, the critical inner pressure pcrmaxS-D can be also obtained by substituting (44) and (50) into (32).(51)pcrmaxS-D=σR3-σcPeak1-ξT1-ξt1-ξ+σcPeak1-ξ+α1+νp0-σR3βs+ξt1+βs+α1+νp0-σR31-ξT1+βp-α1+νp0-σR3βs+ξ+α1+νp0-σR31-ξT1+βpt1-ξ.


### 3.8. The Special Cases

#### 3.8.1. Li’s Formula

WhenΔγ∗=0, T=limΔγ∗=0⁡T=1, and βs=1, the EPSDM degenerates for ESM. By solving (45) and (46), the radii of damage zone and strain-softening zone will be, respectively, rewritten as(52a)R1=R0σR3-σcPeak/1-ξ-α1+νp0-σR3/1+ξ-α1+νp0-σR3/1-ξt1-ξ+α1+νp0-σR3/1+ξt2+α1+νp0-σR3/1-ξ+σcPeak-σcR/1-ξpin-σcR/1-ξ1/ξ-1(52b)R2=σcPeak-σcRα1+νp0-σR3+11/2R1.Equations (52a) and (52b) are Li’s formula [25].Then, the critical inner pressurepcrmaxS-D between strain-softening and damage zone can be also rewritten as follows:(53)pcrmaxS-D=σR3-σcPeak1-ξt1-ξ+σcPeak1-ξ+α1+νp0-σR31+ξt2+α1+νp0-σR31-ξ-α1+νp0-σR31+ξ+α1+νp0-σR31-ξt1-ξ.


#### 3.8.2. Kastner’s Formula

Whenθ=0∘, α=0, and σcPeak=σcR, the EPSDM degenerates for EPM. Combining with (16) and the stress boundary condition (σrp)r=R0=pin at r=R0, the postpeak failure zone radii R3 can be determined by the following equation:(54)R3=R02p0-σcPeak/1+ξ-σcPeak/1-ξpin-σcPeak/1-ξ1/ξ-1.Equation (54) is Kastner’s formula [8]


#### 3.8.3. Wilson’s Formula

WhenΔγ∗=0, α→∞, T=limΔγ∗=0⁡T=1, and t=limα→∞⁡t=1, the EPSDM degenerates for EBM. According to (46), the damage zone radius R1 will be rewritten as(55)R1=R02p0-σcPeak/1+ξ-σcR/1-ξpin-σcR/1-ξ1/ξ-1.Equation (55) is the formula of Wilson [13].


## 3.1. Basic Equations

The equilibrium differential equation for axisymmetric plane strain problem in the “i” zone can be expressed as (ignoring the body force of rock mass) [4, 9, 12, 16, 22, 25](4)dσridr+σri-σθir=0where σri and σθi are radial stress and tangential stress in the “i” zone, respectively. The subscript symbol “i” represents different zone of surrounding rock, which can be replaced by the symbols “e,” “p,” “s,”and “c.”Based on the small deformation assumption, the geometric equation can be denoted as [9, 17, 24](5)εri=duridr,εθi=urirwhere εri and εθi are, respectively, radial and circumferential strain in different zones and uri represents the radial displacement.Supposing that the volume of the rock mass is changing in the postpeak failure zone, the plastic-strain relationships can be established by adopting a nonassociated linear flow rule and small strain theory as follows [11, 14, 22]:(6)εri+βiεθi=0where βi=1+sin⁡ψi/1-sin⁡ψi which is the dilation coefficient in the “i” zone. ψi is the dilation angle.


## 3.2. Elastic Zone

According to the elasticity theory, the stress functionΦ can be expressed as [30](7)Φ=Mln⁡r+Nr2where M and N are the constants. The radial and tangential stresses are then given by(8)σre=1rdΦdr=Mr2+2N(9)σθe=d2Φdr2=-Mr2+2N.Substituting the boundary conditionsσre=p0 at r→∞ and σre=σR3 at r=R3 into (8), where σR3 is the radial stress at the elastic-peak plastic zone interface, (8) and (9) can be rewritten as(10)σre=p0+σR3-p0R32r2(11)σθe=p0-σR3-σ0R32r2.Whenr=R3, (10) and (11) should satisfy (1). Thus, σR3 can be deduced as follows:(12)σR3=2p0-σcPeak1+ξ.The displacement and strains can be easily obtained as follows [25]:(13)ure=A0R32r,εre=-A0R32r2,εθe=A0R32r2where A0=(1+ν)(σ0-σR3)/E. ν and E are Poisson’s ratio and Young’s modulus. When r=R3, the radial displacement uR3 at the elastic-peak plastic zone interface is(14)uR3=A0R3.


## 3.3. Peak Plastic Zone

Substituting (1) into (4), the equilibrium differential equation in peak plastic zone can be rewritten as follows [25]:(15)dσrpdr+1-ξσrp-σcPeakr=0.Solving (15), the radial stress in peak plastic zone can be derived by considering the boundary conditions σrp=σR3 at r=R3.(16)σrp=σR3-σcPeak1-ξrR3ξ-1+σcPeak1-ξ.And then, introducing (16) into (1), the circumferential stress will be also be deduced as(17)σθp=ξσR3-σcPeak1-ξrR3ξ-1+σcPeak1-ξ.The radial contact stressesσR2 at the peak plastic-strain-softening zone interface can be easily obtained by substituting r=R2 into (16):(18)σR2=σR3-σcPeak1-ξR2R3ξ-1+σcPeak1-ξ.Then, the displacement differential equation can be also given by integrating (5) and (6) [15, 24].(19)durpdr+βpurpr=0.Solving (19), the radial displacement will be deduced in peak plastic zone by considering the boundary condition urp=uR3 at r=R3.(20)urp=A0R31+βpr-βp.According to (5), the strains in peak plastic zone can be also obtained as follows:(21)εrp=-A0βpR31+βpr-βp-1(22)εθp=A0R31+βpr-βp-1.According to (20), the radial displacement uR2 at the peak plastic-strain-softening zone interface can be achieved as follows:(23)uR2=A0R31+βpR2-βp.In addition, the maximum principal strainε1sPeak at the peak plastic-strain-softening zone interface can be also deduced by substituting r=R2 into (20) and (22).(24)ε1sPeak=A0R31+βpR2-βp-1.


## 3.4. Strain-Softening Zone

Based on the calculation method of peak plastic zone, the radial displacement in strain-softening zoneurs is easily acquired as follows [24]:(25)urs=A0R31+βpR2βs-βpr-βs.Then, the radial and tangential strain in strain-softening zone are also deduced by integrating (5) and (25).(26)εrs=-A0βsR31+βpR2βs-βpr-βs-1(27)εθs=A0R31+βpR2βs-βpr-βs-1.Meanwhile, according to (25), the radial displacement uR1 at the strain-softening-damage zone interface can be obtained as well.(28)uR1=A0R31+βpR2βs-βpR1-βs.In strain-softening zone, it is assumed that the strength attenuation of rock mass is only related to cohesion (c). Then, the compressive strength at any point can be expressed as(29)σc¯=σcPeak-σcPeak-σcRεcR-ε1sPeakεθs-ε1sPeak=σcPeak-λεθs-ε1sPeakwhere εcR is the maximum principal strain at the point “C” in Figure 3. λ may be called strain-softening modulus and can be determined by the slope of “BC” segment in Figure 3. Introducing (24) and (27) into (29), it can be rewritten as(30)σc¯=σcPeak-α1+νp0-σR3R31+βpR2-βp-1R2r1+βs-1where α=λ/E, which can be defined as a strain-softening coefficient.Substituting (2) and (30) into (4), the equilibrium differential equation in strain-softening zone can be expressed as(31)dσrsdr+1-ξσrsr-σcPeak-α1+νp0-σR3R31+βpR2-βp-1R2/r1+βs-1r=0.Combined with the boundary conditionσrs=σR2 at r=R2, the radial stress in strain-softening zone can be derived by solving (31).(32)σrs=σR3-σcPeak1-ξR3r1-ξ+σcPeak1-ξ+α1+νp0-σR3βs+ξR3R21+βpR2r1+βs-R2r1-ξ+α1+νp0-σR31-ξR3R21+βp1-R2r1-ξ.The circumferential stress is also obtained by substituting (30) and (32) into (2).(33)σθs=ξσR3-σcPeak1-ξR3r1-ξ+σcPeak1-ξ+α1+νp0-σR3βs+ξR3R21+βpR2r1+βs-R2r1-ξ+α1+νp0-σR31-ξR3R21+βp1-R2r1-ξ+σcPeak-α1+νp0-σR3R3R21+βpR2r1+βs-1.


## 3.5. Damage Zone

Substituting (3) into (4), the equilibrium differential equation in damage zone can be rewritten as follows:(34)dσrcdr+1-ξσrc-σcRr=0.Solving (34) and considering the boundary conditions σrc=pin at r=R0, the radial stress in damage zone can be derived as(35)σrc=pin-σcR1-ξrR0ξ-1+σcR1-ξ.Then, introducing (35) into (3), the circumferential stress in damage zone will be also obtained as(36)σθc=ξpin-σcR1-ξrR0ξ-1+σcR1-ξ.As shown in (25), the radial displacement of damage zone can be deduced by considering the boundary condition (urc)r=R1=uR1.(37)urc=A0R31+βpR2βs-βpR1βc-βsr-βc.The radial and circumferential strains in damage zone can be also obtained by substituting (37) into (5).(38)εrc=-βcA0R31+βpR2βs-βpR1βc-βsr-βc-1(39)εθc=A0R31+βpR2βs-βpR1βc-βsr-βc-1.According to (37), when r=R0, the surface displacement of surrounding rock u0 is as follows:(40)u0=A0R31+βpR2βs-βpR1βc-βsR0-βc.


## 3.6. Radius of Postpeak Failure Zone

In order to obtain the closed-form solution of the stresses and deformation of EPSDM, the postpeak failure zone radius (R1,R2,R3) should be firstly determined. According to the stress contact condition (σrs)r=R1=(σrc)r=R1 at r=R1, one of the relationship among R1, R2, and R3 can be established by combining with (32) and (35).(41)σR3-σcPeak1-ξR3R11-ξ+σcPeak1-ξ+α1+νp0-σR3βs+ξR3R21+βpR2R11+βs-R2R11-ξ+α1+νp0-σR31-ξR3R21+βp1-R2R11-ξ-pin-σcR1-ξR1R0ξ-1-σcR1-ξ=0.According to (30), when r=R1, the residual strength can be expressed as follows:(42)σcR=σcPeak-α1+νp0-σR3R31+βpR2-βp-1R2R11+βs-1.In addition to (41) and (42), a new condition is still needed to make the EPSDM form a closed solution. Zhang and Jiang et al. assumed that the surrounding rock begins to enter the next stage when the plastic shear strain or the equivalent plastic shear strain in peak plastic zone reaches a critical value [17, 18]. This paper assumes that when the plastic shear stress increment Δγ∗ between the yield points “A” and “B” (see Figure 3) satisfies certain conditions, the surrounding rock will enter the strain-softening state. Thus, the other relationship among R3, R2, and R1 can be expressed as(43)Δγ∗=γB∗-γA∗=ε1sPeak-ε3sPeak-ε1pPeak-ε3pPeak=εθp-εrpr=R2-εθp-εrpr=R3where γB∗ and γA∗ are the plastic shear strain at points “B” and “A,” respectively. ε1pPeak is the maximum principal strain at the point “A.” ε3sPeak and ε3pPeak are the minimum principal strain at the points “B” and “A.” As shown in Figure 3, the maximum and minimum principal strain of rock mass can be easily obtained in the triaxial test. Substituting (21) and (22) into (43), it can be rewritten as(44)R3=Δγ∗A01+βp+11/1+βpR2=TR2.Introducing (44) into (42), the following expression can be obtained as follows:(45)R2=σcPeak-σcRα1+νp0-σR3T1+βp+11/1+βsR1=tR1.Then, the damage zone radius can be deduced by substituting (44) and (45) into (41).(46)R1=R0σR3-σcPeak1-ξT1-ξt1-ξ-α1+νp0-σR3βs+ξ+α1+νp0-σR31-ξT1+βpt1-ξ+α1+νp0-σR3βs+ξt1+βs+α1+νp0-σR31-ξT1+βp+σcPeak-σcR1-ξ×pin-σcR1-ξ-11/ξ-1.Subsequently, the radiusesR3 and R2 can be also derived by substituting (46) into (45) and (44), respectively.


## 3.7. Critical Inner Pressure under EPSDM

According to (44),(45), and (46), when the relationship between postpeak failure zone radius (R1,R2,R3) and inner radius (R0) satisfies certain condition, the rock mass will be in the critical state.(i)
WhenR0=R1=R2=R3, the rock mass is in the critical state that the peak plastic is just not arisen yet. Then, the critical inner pressure pcrmaxE-P at this state will be deduced according to (12)(47)pcrmaxE-P=σR3=2p0-σcPeak1+ξ.(ii)
WhenR0=R1=R2<R3, the rock mass is in the critical state that the peak plastic zone has reached the maximum and the stain softening zone has not yet appeared. Consequently, the critical radii of peak plastic zone can be obtained by (44).(48)R3=Δγ∗A01+βp+11/1+βpR0.Then, introducing (48) into (16), the critical inner pressure pcrmaxP-S at this state will be deduced as follows:(49)pcrmaxP-S=σR3-σcPeak1-ξΔγ∗A01+βp+1-ξ-1/1+βp+σcPeak1-ξ.(iii)
WhenR0=R1<R2<R3, the rock mass is in the critical state that the damage zone begins to form. According to (45), the strain-softening zone radius can be expressed as(50)R2=σcPeak-σcRα1+νp0-σR3T1+βp+11/1+βsR1=tR0.The peak plastic zone radius can be calculated by (44) at this state. Then, when r=R0, the critical inner pressure pcrmaxS-D can be also obtained by substituting (44) and (50) into (32).(51)pcrmaxS-D=σR3-σcPeak1-ξT1-ξt1-ξ+σcPeak1-ξ+α1+νp0-σR3βs+ξt1+βs+α1+νp0-σR31-ξT1+βp-α1+νp0-σR3βs+ξ+α1+νp0-σR31-ξT1+βpt1-ξ.


## 3.8. The Special Cases

### 3.8.1. Li’s Formula

WhenΔγ∗=0, T=limΔγ∗=0⁡T=1, and βs=1, the EPSDM degenerates for ESM. By solving (45) and (46), the radii of damage zone and strain-softening zone will be, respectively, rewritten as(52a)R1=R0σR3-σcPeak/1-ξ-α1+νp0-σR3/1+ξ-α1+νp0-σR3/1-ξt1-ξ+α1+νp0-σR3/1+ξt2+α1+νp0-σR3/1-ξ+σcPeak-σcR/1-ξpin-σcR/1-ξ1/ξ-1(52b)R2=σcPeak-σcRα1+νp0-σR3+11/2R1.Equations (52a) and (52b) are Li’s formula [25].Then, the critical inner pressurepcrmaxS-D between strain-softening and damage zone can be also rewritten as follows:(53)pcrmaxS-D=σR3-σcPeak1-ξt1-ξ+σcPeak1-ξ+α1+νp0-σR31+ξt2+α1+νp0-σR31-ξ-α1+νp0-σR31+ξ+α1+νp0-σR31-ξt1-ξ.


### 3.8.2. Kastner’s Formula

Whenθ=0∘, α=0, and σcPeak=σcR, the EPSDM degenerates for EPM. Combining with (16) and the stress boundary condition (σrp)r=R0=pin at r=R0, the postpeak failure zone radii R3 can be determined by the following equation:(54)R3=R02p0-σcPeak/1+ξ-σcPeak/1-ξpin-σcPeak/1-ξ1/ξ-1.Equation (54) is Kastner’s formula [8]


### 3.8.3. Wilson’s Formula

WhenΔγ∗=0, α→∞, T=limΔγ∗=0⁡T=1, and t=limα→∞⁡t=1, the EPSDM degenerates for EBM. According to (46), the damage zone radius R1 will be rewritten as(55)R1=R02p0-σcPeak/1+ξ-σcR/1-ξpin-σcR/1-ξ1/ξ-1.Equation (55) is the formula of Wilson [13].


## 3.8.1. Li’s Formula

WhenΔγ∗=0, T=limΔγ∗=0⁡T=1, and βs=1, the EPSDM degenerates for ESM. By solving (45) and (46), the radii of damage zone and strain-softening zone will be, respectively, rewritten as(52a)R1=R0σR3-σcPeak/1-ξ-α1+νp0-σR3/1+ξ-α1+νp0-σR3/1-ξt1-ξ+α1+νp0-σR3/1+ξt2+α1+νp0-σR3/1-ξ+σcPeak-σcR/1-ξpin-σcR/1-ξ1/ξ-1(52b)R2=σcPeak-σcRα1+νp0-σR3+11/2R1.Equations (52a) and (52b) are Li’s formula [25].Then, the critical inner pressurepcrmaxS-D between strain-softening and damage zone can be also rewritten as follows:(53)pcrmaxS-D=σR3-σcPeak1-ξt1-ξ+σcPeak1-ξ+α1+νp0-σR31+ξt2+α1+νp0-σR31-ξ-α1+νp0-σR31+ξ+α1+νp0-σR31-ξt1-ξ.


## 3.8.2. Kastner’s Formula

Whenθ=0∘, α=0, and σcPeak=σcR, the EPSDM degenerates for EPM. Combining with (16) and the stress boundary condition (σrp)r=R0=pin at r=R0, the postpeak failure zone radii R3 can be determined by the following equation:(54)R3=R02p0-σcPeak/1+ξ-σcPeak/1-ξpin-σcPeak/1-ξ1/ξ-1.Equation (54) is Kastner’s formula [8]


## 3.8.3. Wilson’s Formula

WhenΔγ∗=0, α→∞, T=limΔγ∗=0⁡T=1, and t=limα→∞⁡t=1, the EPSDM degenerates for EBM. According to (46), the damage zone radius R1 will be rewritten as(55)R1=R02p0-σcPeak/1+ξ-σcR/1-ξpin-σcR/1-ξ1/ξ-1.Equation (55) is the formula of Wilson [13].


## 4. Example Study

### 4.1. Case I: Verification for EPSDM

In this paper, the analytical solution for the circular opening based on the EPSDM can be regarded as an extension of strain-softening model. For the contraction problem of the circular hole excavated in the ESM rock mass, Li et al. [25] also carried out relevant studies. However, the influence of plastic bearing behavior and dilation coefficient on rock mass was neglected. The geometrical and physical parameters of circular opening are shown in Table 1. For guaranteeing the same condition as Li’s research, the plastic shear strain increment is assumed as 0.Table 1
Geometrical and physical parameters of circular opening.
Parameters Value Radius of opening,R0 (m) 0.1 In-situ stress,p0 (MPa) 15 Internal pressure,pin (MPa) 0 Young’s modulus,E (MPa) 1500 Dilatancy coefficient,βi 1.0 Poisson’s ratio,ν 0.3 Strain softening coefficient,α 0.5 internal friction angle,φ (°) 30 Peak compressive strength,σcpeak (MPa) 8 Residual compressive strength,σcR (MPa) 1The stress and displacement of surrounding rock along the radius direction are shown in Figure4 and Figure 5. In addition, Table 2 also lists the calculated results under EPSDM, in accordance with Li’s closed solution. Therefore, the solution by Li et al. [25] is a special case of this paper. In other words, the EPSDM has more extensive practicality.Table 2
Calculation results of postpeak failure zone radii, critical inner pressure, and displacement. (Note that Li’s calculations dates are in brackets.)
Parameter R 2 / R 0 or R1/R0 p c r  m a x E - P / p 0 or pcrmaxS-D/p0 u 0 / a Δ γ ∗ = 0 2.834 (2.834) 0.367 (0.367) 0.066 (0.066) 1.940 (1.940) 0.092 (0.092)Figure 4
Stress distribution along the radius direction.Figure 5
Displacement distribution along the radius direction.


### 4.2. Case II: Compared with the Other Constitutive Models

As previously mentioned, the constitutive model of EPSDM reflects all the features of the EPM, EBM, EPBM, and ESM. Therefore, the closed solution of EPSDM is different from the above four models. In order to study the evolution law of surrounding rock states under different constitutive models, the input data for soft rock mass is shown in Table3.Table 3
Geometry and material mechanics parameters.
Parameters Value Radius of opening,R0 (m) 1 In situ stress,p0 (MPa) 1 Young’s modulus,E (MPa) 5000 Dilatancy angle,ψi (°) 0 or 30° Poisson’s ratio,ν 0.2 Strain-softening coefficient,α 2.5 internal friction angle,φ (°) 35 Initial cohesion,c (MPa) 0.276 Residual cohesion,cr (MPa) 0.055 Shear stain increment,Δγ∗ 2.5e-5The influences of constitutive model on the postpeak failure zone radii, surface displacement, and stresses distribution are presented in Figures6, 7, and 8, respectively. According to (47), (49), and (51), there is only elastic zone around the circular opening when pcrmaxE-P≤pin<p0. Then, the surrounding rock will consist of peak plastic zone and elastic zone when pcrmaxP-S≤pin<pcrmaxE-P. Next, the strain-softening zone will be gradually formed when the inner pressure satisfies pcrmaxS-D≤pin<pcrmaxP-S. Finally, the surrounding rock will be composed of damage zone, strain-softening zone, peak plastic zone, and elastic zone with the continuous decrease of inner pressures. From Figures 6, 7, and 8, it can be seen that the closed solution of EPSDM is in accordance with EPM for βi=1.0 and βi=3.0 when pcrmaxP-S≤pin<pcrmaxE-P. In addition, once the inner pressure is fully released, the surface displacement of surrounding rock will show the characteristics of EBM>EPBM>ESM>EPSDM>EPM (see Figure 7). Therefore, the small plastic bearing capacity can effectively decrease the postpeak failure zone radii as well as surface displacement.Figure 6
The influence of inner pressure on the radius of postpeak failure zone.
(a)
The radii evolution law (ψi=0°; βi=1.0) (b)
The radii evolution law (ψi=30°; βi=3.0)Figure 7
The influence of inner pressure on the surface displacement.
(a)
The displacement evolution law (ψi=0°; βi=1.0) (b)
The displacement evolution law (ψi=30°; βi=3.0)Figure 8
Stresses distribution under different constitutive models.
(a)
The stresses evolution law (ψi=0°; βi=1.0) (b)
The stresses evolution law (ψi=30°; βi=3.0)The support parameter design is closely related to the rock mass properties in deep underground engineering. The rock mass with certain plastic bearing capacity can effectively resist the influence of external load on the tunnel deformation. Therefore, the support strength designs for the above different characteristics tunnel should satisfy the relationship with EBM>EPBM>ESM>EPSDM>EPM.In general, most jointed and fractured rock masses show a certain plastic bearing capacity before strain softening. Therefore, the support parameters design by using EPSDM may be more reasonable.Whenpin=0, the postpeak failure zone radius and surface displacement of EPSDM are shown in Table 4. It can be seen that the postpeak failure zone radius and surface displacement increase with an increasing dilation coefficient, respectively. For instance, when βi is changed from 1.0 to 3.0, the dimensionless radius (R3/R0, R2/R0, and R1/R0) and surface displacement (2u0G/[R0(p0-σR3)]), respectively, increase 12.65%, 15.29%, 23.62%, and187.57%. Therefore, the dilation coefficient noticeably influences postpeak failure zone radius and surface displacement, particularly the damage zone radii and tunnel wall convergence.Table 4
The postpeak failure zone radius and surface displacement of EPSDM whenpin=0.
R 3 / R 0 R 2 / R 0 R 1 / R 0 2 u 0 G / [ R 0 ( p 0 - σ R 3 ) ] β i = 1.0 1.336 1.295 1.122 1.786 β i = 3.0 1.505 1.493 1.387 5.136


### 4.3. Case III: Application in No. 4 Coal Mine of Saier Group

The haulage roadway by the No. 4 coal mine of Saier Group in China was excavated in the EPSDM rock mass. It is buried at about 320m underground, the maximum vertical stress is 15.16MPa, and two horizontal principal stresses are, respectively, 18.67MPa and 22.43MPa. Thus, the average value of the in situ stress (p0) is 18.75MPa. The equivalent excavation radius (R0) is 3.54m, and Poisson’s ratio (ν) is 0.23. According to Figure 2, the relationship between maximum and minimum principal stress for the sandy mudstone can be obtained (see Figure 9). Meanwhile, the parameters variation (E and λ) with confining pressure is also shown in Figure 10. From the above two figures, It can be seen that the peak plastic strength (σcPeak) and residual strength (σcR) are, respectively, 20.68MPa and 5.15MPa, the parameter (ξ) is about 3.38, Young’s modulus (E) is 1.99GPa, and strain-softening coefficient (α) is 3.56. In addition, it is assumed that the plastic shear strain increment (Δγ∗) is 0.0001, and the dilation coefficient is 1.0.Figure 9
Fitting curves of M-C strength criterion for sandy mudstone.Figure 10
Fitting curves between E,λ and confining pressure.Figure11 shows the failure region tested by JL-IDOI (A) intelligent optical borehole imager at the haulage roadway. In the test, nine boreholes were arranged around the tunnel. According to the test result, the average failure range and surface displacement are 1.48m and 6.28cm, respectively. In addition, the calculated peak plastic zone, strain-softening zone, and damage zone radius by EPSDM are R3=5.04m, R2=5.03m, and R1=4.52, respectively, with equivalent excavation radii R0=3.54. Meanwhile, the calculated surface displacement is u0=6.62cm. It should be noted that the postpeak failure region should be determined by the strain-softening and damage zones radius. Hence, the calculated maximum failure region is 1.49m and is in good accordance with the field test results.Figure 11
Failure region tested by JL-IDOI (A) intelligent optical borehole imager.


## 4.1. Case I: Verification for EPSDM

In this paper, the analytical solution for the circular opening based on the EPSDM can be regarded as an extension of strain-softening model. For the contraction problem of the circular hole excavated in the ESM rock mass, Li et al. [25] also carried out relevant studies. However, the influence of plastic bearing behavior and dilation coefficient on rock mass was neglected. The geometrical and physical parameters of circular opening are shown in Table 1. For guaranteeing the same condition as Li’s research, the plastic shear strain increment is assumed as 0.Table 1
Geometrical and physical parameters of circular opening.
Parameters Value Radius of opening,R0 (m) 0.1 In-situ stress,p0 (MPa) 15 Internal pressure,pin (MPa) 0 Young’s modulus,E (MPa) 1500 Dilatancy coefficient,βi 1.0 Poisson’s ratio,ν 0.3 Strain softening coefficient,α 0.5 internal friction angle,φ (°) 30 Peak compressive strength,σcpeak (MPa) 8 Residual compressive strength,σcR (MPa) 1The stress and displacement of surrounding rock along the radius direction are shown in Figure4 and Figure 5. In addition, Table 2 also lists the calculated results under EPSDM, in accordance with Li’s closed solution. Therefore, the solution by Li et al. [25] is a special case of this paper. In other words, the EPSDM has more extensive practicality.Table 2
Calculation results of postpeak failure zone radii, critical inner pressure, and displacement. (Note that Li’s calculations dates are in brackets.)
Parameter R 2 / R 0 or R1/R0 p c r  m a x E - P / p 0 or pcrmaxS-D/p0 u 0 / a Δ γ ∗ = 0 2.834 (2.834) 0.367 (0.367) 0.066 (0.066) 1.940 (1.940) 0.092 (0.092)Figure 4
Stress distribution along the radius direction.Figure 5
Displacement distribution along the radius direction.


## 4.2. Case II: Compared with the Other Constitutive Models

As previously mentioned, the constitutive model of EPSDM reflects all the features of the EPM, EBM, EPBM, and ESM. Therefore, the closed solution of EPSDM is different from the above four models. In order to study the evolution law of surrounding rock states under different constitutive models, the input data for soft rock mass is shown in Table3.Table 3
Geometry and material mechanics parameters.
Parameters Value Radius of opening,R0 (m) 1 In situ stress,p0 (MPa) 1 Young’s modulus,E (MPa) 5000 Dilatancy angle,ψi (°) 0 or 30° Poisson’s ratio,ν 0.2 Strain-softening coefficient,α 2.5 internal friction angle,φ (°) 35 Initial cohesion,c (MPa) 0.276 Residual cohesion,cr (MPa) 0.055 Shear stain increment,Δγ∗ 2.5e-5The influences of constitutive model on the postpeak failure zone radii, surface displacement, and stresses distribution are presented in Figures6, 7, and 8, respectively. According to (47), (49), and (51), there is only elastic zone around the circular opening when pcrmaxE-P≤pin<p0. Then, the surrounding rock will consist of peak plastic zone and elastic zone when pcrmaxP-S≤pin<pcrmaxE-P. Next, the strain-softening zone will be gradually formed when the inner pressure satisfies pcrmaxS-D≤pin<pcrmaxP-S. Finally, the surrounding rock will be composed of damage zone, strain-softening zone, peak plastic zone, and elastic zone with the continuous decrease of inner pressures. From Figures 6, 7, and 8, it can be seen that the closed solution of EPSDM is in accordance with EPM for βi=1.0 and βi=3.0 when pcrmaxP-S≤pin<pcrmaxE-P. In addition, once the inner pressure is fully released, the surface displacement of surrounding rock will show the characteristics of EBM>EPBM>ESM>EPSDM>EPM (see Figure 7). Therefore, the small plastic bearing capacity can effectively decrease the postpeak failure zone radii as well as surface displacement.Figure 6
The influence of inner pressure on the radius of postpeak failure zone.
(a)
The radii evolution law (ψi=0°; βi=1.0) (b)
The radii evolution law (ψi=30°; βi=3.0)Figure 7
The influence of inner pressure on the surface displacement.
(a)
The displacement evolution law (ψi=0°; βi=1.0) (b)
The displacement evolution law (ψi=30°; βi=3.0)Figure 8
Stresses distribution under different constitutive models.
(a)
The stresses evolution law (ψi=0°; βi=1.0) (b)
The stresses evolution law (ψi=30°; βi=3.0)The support parameter design is closely related to the rock mass properties in deep underground engineering. The rock mass with certain plastic bearing capacity can effectively resist the influence of external load on the tunnel deformation. Therefore, the support strength designs for the above different characteristics tunnel should satisfy the relationship with EBM>EPBM>ESM>EPSDM>EPM.In general, most jointed and fractured rock masses show a certain plastic bearing capacity before strain softening. Therefore, the support parameters design by using EPSDM may be more reasonable.Whenpin=0, the postpeak failure zone radius and surface displacement of EPSDM are shown in Table 4. It can be seen that the postpeak failure zone radius and surface displacement increase with an increasing dilation coefficient, respectively. For instance, when βi is changed from 1.0 to 3.0, the dimensionless radius (R3/R0, R2/R0, and R1/R0) and surface displacement (2u0G/[R0(p0-σR3)]), respectively, increase 12.65%, 15.29%, 23.62%, and187.57%. Therefore, the dilation coefficient noticeably influences postpeak failure zone radius and surface displacement, particularly the damage zone radii and tunnel wall convergence.Table 4
The postpeak failure zone radius and surface displacement of EPSDM whenpin=0.
R 3 / R 0 R 2 / R 0 R 1 / R 0 2 u 0 G / [ R 0 ( p 0 - σ R 3 ) ] β i = 1.0 1.336 1.295 1.122 1.786 β i = 3.0 1.505 1.493 1.387 5.136


## 4.3. Case III: Application in No. 4 Coal Mine of Saier Group

The haulage roadway by the No. 4 coal mine of Saier Group in China was excavated in the EPSDM rock mass. It is buried at about 320m underground, the maximum vertical stress is 15.16MPa, and two horizontal principal stresses are, respectively, 18.67MPa and 22.43MPa. Thus, the average value of the in situ stress (p0) is 18.75MPa. The equivalent excavation radius (R0) is 3.54m, and Poisson’s ratio (ν) is 0.23. According to Figure 2, the relationship between maximum and minimum principal stress for the sandy mudstone can be obtained (see Figure 9). Meanwhile, the parameters variation (E and λ) with confining pressure is also shown in Figure 10. From the above two figures, It can be seen that the peak plastic strength (σcPeak) and residual strength (σcR) are, respectively, 20.68MPa and 5.15MPa, the parameter (ξ) is about 3.38, Young’s modulus (E) is 1.99GPa, and strain-softening coefficient (α) is 3.56. In addition, it is assumed that the plastic shear strain increment (Δγ∗) is 0.0001, and the dilation coefficient is 1.0.Figure 9
Fitting curves of M-C strength criterion for sandy mudstone.Figure 10
Fitting curves between E,λ and confining pressure.Figure11 shows the failure region tested by JL-IDOI (A) intelligent optical borehole imager at the haulage roadway. In the test, nine boreholes were arranged around the tunnel. According to the test result, the average failure range and surface displacement are 1.48m and 6.28cm, respectively. In addition, the calculated peak plastic zone, strain-softening zone, and damage zone radius by EPSDM are R3=5.04m, R2=5.03m, and R1=4.52, respectively, with equivalent excavation radii R0=3.54. Meanwhile, the calculated surface displacement is u0=6.62cm. It should be noted that the postpeak failure region should be determined by the strain-softening and damage zones radius. Hence, the calculated maximum failure region is 1.49m and is in good accordance with the field test results.Figure 11
Failure region tested by JL-IDOI (A) intelligent optical borehole imager.


## 5. Conclusion

Considering the influence of plastic bearing capacity on surrounding rock state, a new strain-softening constitutive model was proposed. Then, a closed-form solution for circular opening based on the newly proposed constitutive model was deduced and the validity of the solution is also verified. Moreover, the postpeak failure zone radius, surface displacement, and stresses distribution of circular opening were analyzed under different constitutive models. Finally, an engineering case shows that this paper’s solution is consistent with the field test results. The primary conclusions can be summarized as follows:(1)
The peak plastic zone represents the plastic bearing behavior of strain-softening rock mass, and once plastic shear strain increment reaches the critical failure condition, the strain-softening zone begins to occur. Compared with other solutions, only when the plastic shear strain increment was zero, the closed solution of the EPSDM can be converted into ESM’s solution. Similarly, when the plastic shear strain increment was zero and strain-softening coefficient was large enough, the EPSDM’s solution will transform into EBM’s solution; only when strain-softening coefficient was zero, the calculated results of the EPSDM was found to be in accordance with the closed-form solution of the EPM.(2)
When the inner pressure is fully released, the maximum postpeak failure radii and surface displacement of surrounding rock satisfy the relationship with EBM>ESM>EPSDM>EPM. It is shown that the plastic bearing behavior of strain-softening rock mass can effectively decrease the postpeak failure zone radii as well as surface displacement. Therefore, the design of support parameters by using this paper’s EPSDM may be more reasonable for peak plastic-strain-softening rock mass.(3)
The rock mass with certain plastic bearing capacity can effectively resist the influence of external load on the tunnel deformation. Therefore, the support strength designs should satisfy the relationship with EBM tunnel>EPBM tunnel>ESM tunnel>EPSDM tunnel>EPM tunnel.(4)
The dilation coefficient noticeably influences postpeak failure zone radius and surface displacement, particularly the damage zone radii and tunnel wall convergence. Therefore, the support parameters design should take into account the influence of the dilation coefficient.



---

*Source: 4521492-2018-07-02.xml*