Research Article  Open Access
Elastoplastic Analysis of TwoLayered Circular Lining Based on the Unified Strength Theory
Abstract
Based on the Unified Strength Theory (UST), elastoplastic analysis of the twolayered circular lining is carried out. The stresses, displacements, and the elastic and plastic zones in both layers are discussed under different values of Young’s moduli of the inner and outer layers. The results reveal that, compared to the singlelayered lining, the tangential stress distributions in the twolayered linings are more reasonable along the radial direction, which is beneficial to enhance the overall elastic and plastic ultimate bearing capacities. When considering the intermediate stress (i.e., the axial load), the elastic ultimate bearing capacity will be higher. However, the plastic ultimate bearing capacity remains unchanged. Moreover, a comparison between the Unified Strength Theory and Tresca Criterion is analyzed as well.
1. Introduction
Recently, along with the construction of underground caverns in depth, the support technique needs to develop accordingly, to assess the stability of the structures and surrounding rock [1]. Thickwalled hollow cylinder is one of the widely used support structures. To optimize the structural design and make the best usage of construction materials, pressurized hollow cylinders are common cases in deepburied circular caverns.
For the elastic plane problem on circular structures, a lot of literature can be found [2]. The wellknown Lame solution presents the stress and displacement fields for a singlelayered hollow cylinder subjected to axisymmetric loads. Wu [3] solved the solutions of the multilayered case. For the nonaxisymmetric case, Lu et al. [4], by using the semiinverse method, presented the elastic plane stress solution of a lined vertical shaft in isotropic ground. Wu and Lu [5] solved the stress fields for the plane problem of a thickwalled cylinder, by utilizing the complex variable function method, and they found that there exists a tensile tangential stress.
For the elastoplastic analysis on singlelayered thickwalled hollow cylinder, there also have been many solutions on stress, displacement distributions, and elastic, plastic zones based on different strength theories [6–11]. Durban and Kube [12] investigated the problem of a thickwalled cylindrical tube subjected to internal pressure using a finite strain elastoplastic flow theory. Detournay and Fairhurst [13] describe a twodimensional elastoplastic model of a long, cylindrical cavity in an infinite rock mass subject to nonhydrostatic farfield stress loading. To discuss the plastic ultimate bearing capacity, Gao et al. [14] considered the size effect of strainhardening level. And in the fields of elasticbrittleplastic analysis, Leu [15] applied sequential limit analysis, dealing with the quasistatic problem of hardening and weakening behaviors. Lee and Pietruszczak [16] calculated the stresses and displacements fields of a circular tunnel excavated in strainsoftening MohrCoulomb or generalized HoekBrown rock mass. Wang et al. [17, 18] and Zhang et al. [19, 20] presented a solution of a spherical cavity in brittle–plastic medium under hydrostatic in situ stress, in which the linear MohrCoulomb and nonlinear HoekBrown yield criteria were applied.
For nonplane problems, Lu et al. [21] presented the elastoplastic solutions for a circular tunnel with considering the axial stress, showing that the distributions of plastic zones are dependent on the axial and horizontal in situ stresses. Wang et al. [22] suggested an elasticbrittleplastic constitutive model with a nonassociated flow rule.
From the above mentioned literature, there are few results on elastoplastic analysis for the multilayered thickwalled hollow cylinder. We herein will focus on the elastoplastic analysis of a multilayered thickwalled hollow cylinder (see Figure 1). Moreover, the ultimate bearing capacities are to be investigated as well. In Figure 1, E_{i}, μ_{i}, and (i=1 and 2) denote the elastic modulus, Poisson’s ratio, and compressive strength of the inner and outer layers, respectively. The axisymmetric external load on the outer boundary is p. Here, we assume that the tunnel is deepburied and materials of both layers are homogeneous and isotropic and satisfy the ideal elastoplastic model. Moreover, the interface boundary of the two layers is a pure bond, no slide happening along the tangential or longitudinal direction.
(a)
(b)
2. Basic Equations and Stresses and Displacements in Elastic Zones
The cylindrical coordinate system is employed. Moreover, for simplicity, the characteristic of axial symmetry is considered. The straindisplacement relationships for an axisymmetric problem arewhere and are the radial and tangential strains, respectively. Moreover, u is the radial displacement, and r is the radial distance from the center of opening.
The differential equation of equilibrium for an axisymmetric problem iswhere and are the tangential and radial stresses, respectively.
The compatibility equation of strain can be given as
The generalized Hooke’s Law is given as follows:where is the intermediate principal stress.
The expressions of stress components in elastic state for an axisymmetric problem can be given according to the Lame solutions as
The subscript i=1 and 2, denoting the inner and outer layers, respectively (similarly hereinafter). Substituting Eq. (8) into Eq. (6) and considering Eq. (2), the radial displacement in elastic zones is obtained as
Note that the expressions of the stress and displacement components, (8) and (10), are suitable for both inner and outer layers, just with different values of the coefficients A_{i} and B_{i} (i=1 and 2).
The Unified Strength Theory is applied herein. This strength theory can take the intermediate principal stress into consideration. For the axisymmetric plane strain problem, there exist three nonzero principal stresses (i.e., the radial, circumferential, and axial stresses, , , and ). Compression is taken as positive. For the plane strain problem, the axial stress can be given as follows [23]:where m is the parameter reflecting the state of the intermediate principal stress. In elastic zones, . As the inplane plastic strains increase, approaches the average of the minor and major principal stresses (i.e., m_{i} tends to 1). For simplicity, in plastic zones, m is assumed to be 1 [24]. Then, (11) can be expressed as
The parameter m takes 0 and denoting the plane stress and strain problems, respectively. The intermediate principal stress in the elastic zone can be given by substituting (8) and (9) into (11) as
And the radial displacement in elastic zones can be obtained by substituting (14) into (10) as
3. Elastoplastic Analysis
3.1. Elastic Analysis
As the axisymmetric external load p is relatively small, two layers are in the elastic state. Based on the stress boundary conditions on the inner and outer boundaries, and the stress and displacement continuity conditions on the interface between two layers, we have
Based on (16)–(19), the undetermined coefficients A_{i}, B_{i} (i=1, 2) can be solved, which are presented in the Appendix.
3.2. Expressions of Stress and Displacement in the Plastic Zone
The Unified Strength Theory was proposed based on an orthogonal octahedron of a twin shear element model [23, 25]. It is a twinshear element model [26]. It can be expressed using the principal stresses aswhere σ_{1}> σ_{2}> σ_{3}, and α denotes the ratio of the uniaxial tensile to the compressive strengths (i.e., α=f_{t}/f_{c}) and 0<α⩽1 for most engineering materials. In this work, we take α=(1sinφ)/(1+sinφ). Moreover, f_{t} and f_{c} are the uniaxial tensile and compressive strengths, respectively; b reflects the effect of the intermediate principal stress. Theoretically, b can take any value. Actually, the Unified Strength Theory can represent different criteria as α and b take specific values [23]. For example, if b=0 and α=1, it is the Tresca Criterion. If b=1/2 or α=1 and , it is the linear approximations of the von Mises criterion. The MohrCoulomb criterion can be obtained as b=0 and 0<α<1. If b=1 and 0<α<1, it becomes the generalized twin shear stress criterion. Specifically, when α=1 and b=1, it becomes the twin shear stress criterion. Thus, it has been widely used due to the generalization since it was proposed [7, 9, 23, 27, 28].
According to (11), the axial stress σ_{z} should be the intermediate principal stress. The three principal stresses satisfy the relationship, σ_{θ}> σ_{z}> σ_{r}; thus, we have
And we also have
According to (22) and (23), (20) can be satisfied. Then,
Equation (24) can be expressed as the form of yield function as follows:where and .
Combining (3), (11), and (25), we havewhere , .
Combining Eq. (11), substituting Eq. (26) into Eq. (24), we havewhere and .
The flow rule related with the Unified Strength Theory is employed:
According to (31) and (32), we have
When the material enters the plastic state, the elastic and plastic strains constitute the strain, which iswhere and are the radial and tangential stresses in the plastic zone (see (26) and (27)), respectively; σ_{z} is the axial stress (i.e., (11)).
Substituting (11) into (38) and combining (26) and (27), the tangential strain in the plastic zone is derived as follows:where
Substituting (39) into (4), we have
Combining (2) and (41), the radial displacement in the plastic zone can be obtained as
To obtain the stresses and displacements in plastic zones, the coefficients A_{i}, B_{i}, C_{i} and D_{i} (i=1, 2) need to be calculated.
3.3. Elastoplastic Analysis
When a multilayered thickwalled hollow cylinder is subjected to hydrostatic pressure on the outer boundary, the plastic yielding will always initiate from the inner boundary of each layer. Moreover, the extending path of plastic yielding is dependent on the elastic modulus E_{i}, material strength σ_{ci}, and dimensions R_{i}. Specifically, for a twolayered hollow lining, plastic yielding may initiate in three different ways (see Figure 2; and are the dimensions of plastic zones). Hereby, we only discuss the solution of the most complicated situation (Figure 2(c)). The other two simpler cases can be discussed in a similar way.
(a) Plastic zone in inner layer (R0≤ rp≤ R1)
(b) Plastic zone in outer layer (R1≤ rp≤ R2)
(c) Plastic zones in both layers (R0≤ rp1≤ R1, R1≤ rp2≤ R2)
The equations for the undetermined coefficients can be obtained according to the boundary and continuity conditions, which include the stress boundary conditions at the inner boundary of the inner layer r=R_{0}, ; the radial stress and displacement continuity conditions at three interfaces (i.e., at the elasticplastic interface of the inner layer r=r_{p1}, and ; at the interboundary r=R_{1}, and ; at the elasticplastic interface of the outer layer r=r_{p2}, and ); and the stress boundary conditions at the outer boundary of the outer layer r=R_{2}, . According to the conditions, we have
There are eight undetermined variables A_{i}, B_{i}, C_{i}, D_{i} (i=1 and 2) and eight equations (i.e., (43)–(50)). To obtain the radii of the plastic zones (i.e., r_{p1} and r_{p2}), two more equations are required. From the extending path of plastic yielding, the stresses on the inner boundaries of the two elastic zones should also satisfy the Unified Strength Theory. Thus, we have
Substituting (8) and (14) into (51) and (52), we have
With the increase of the load p, the extending path of plastic yielding can be obtained from (53) and (54).
4. Results and Discussions
4.1. Material Properties
For comparing, the geometry and material properties (see Table 1) are taken the same as in [29]. Let K denote the ratio of Young’s moduli; that is, K=E_{1}/E_{2}.

4.2. Effect of K on Stresses, Displacements, and Plastic Zones
Under the parameters listed in Table 1, the radial and tangential stresses are discussed. To compare the stresses and displacements distributions in different situations, we take the value of support pressure p=16 and 20 MPa as examples. In this section, the intermediate principal stress is considered, taking b_{1}=b_{2}=0.5. The ratio K=25/41 and 41/25, denoting two kinds of multilayered structures, and K=41/41 representing the traditional homogenous singlelayered lining. The radial and tangential stresses in the radial direction under different confining pressures are illustrated in Figures 3 and 4. The distributions of the tangential stress differ significantly for the three conditions, but the difference in the radial stress is relatively small.
The tangential stress concentration is the main reason of the plastic yielding. Usually, the largest tangential stress concentration always occurs on the inner boundary of each layer. Besides, the stresses there are twodimensional, so the inner boundary is the most yieldingprone zone. As p=16 MPa (see Figure 3), the plastic yielding only can be found in the area of the inner boundary of the inner layer as K=41/25. For this case, the greatest tangential stress concentration factor σ_{θ}/p=3.538, while, for the other two cases K=41/41 and 25/41, at the same position, the tangential stress concentration factor σ_{θ}/p is 3.125 and 2.486, decreasing by 11.7% and 29.7% compared to the case of K=41/25, respectively. Moreover, for the case K=25/41, the greatest tangential stress concentration, σ_{θ}/p=3.078, occurs at the inner boundary of the outer layer where the radial stress is nonzero, so the stress there is threedimensional. As known to all, the material under threedimensional stress has higher strength than that under twodimensional stress. Thus, for the three cases, the most reasonable stress distribution is in the lining of K=25/41, and the worst is that of K=41/25. As p=20 MPa (see Figure 4), the plastic yielding can be found in the cases of K=41/25 and 41/41, with the thicknesses of 0.097 m and 0.028 m, respectively. However, the whole lining of K=25/41 is still in the elastic state. Thus, according to the stress distributions and the plastic area, we can come to the fact that the structure of K=25/41 is the optimal one in the three cases.
Figure 5 illustrates the extending path of the plastic yielding for K=25/41, 41/41, and 41/25. To illustrate the elastoplastic stages, the radial displacement at the inner boundary of the linings is plotted in Figure 6. With the increase of p, for the case of K=41/25, the plastic yielding first happens in the inner layer, initiating from the inner boundary as p=14.639 MPa (i.e., the elastic ultimate bearing capacity, point H in Figures 5 and 6) and finishing at the outer boundary as p=20.080 MPa (point I). After that, the structure turns into the elastic state (stage IJ) before plastic yielding happens in the outer layer as p=24.425 MPa (point J), and finally the outer layer totally turns into the plastic yielding state as p=26.399 MPa (i.e., the plastic ultimate bearing capacity point K).
For the conventional singlelayered structure (K=41/41), the plastic yielding is found as p=17.360 MPa (point E, the elastic ultimate bearing capacity) at the inner boundary and extends to the outer boundary of the whole lining as p=26.399 MPa (point G, the plastic ultimate bearing capacity). In contrast with the other twolayered structures, the singlelayered lining is considered as a twolayered lining as well, just with the same material parameters in both layers. In this way, point F corresponds to the moment that the plastic yielding extends from the inner layer to the outer layer, of which the confining pressure p=24.425 MPa. Note that the plastic zone extending stage FG is exactly the same as the stage JK of the twolayered lining (K=41/25), but with different radial displacements (Figure 6).
By contrast, for the other twolayered lining (K=25/41), the plastic yielding happens as p=20.560 MPa (point A) from the inner boundary of the outer layer, not from the inner layer. The plastic yielding cannot be found in the inner layer until the confining pressure p increases to 21.769 MPa, when the dimension of the plastic zone in the outer layer is 0.014 m, about 14% of the layer thickness. With the increase of p, the plastic yielding extends in both layers simultaneously (stages B′C and BC′), then the outer layer first turns into the plastic state totally as p=25.925 MPa (point C), and finally the inner layer as p=26.399 MPa (point D, the plastic ultimate bearing capacity).
The case of K=25/41 has the greatest elastic ultimate bearing capacity 20.560 MPa, about 40.4% and 18.4% higher than the cases of K=41/25 and 41/41, respectively. Moreover, the three cases have the same plastic ultimate bearing capacity 26.399 MPa. Thus, from the aspect of the elastic bearing capacity, the lining of K=25/41 performs the best among the three structures. Actually, for a multilayered lining, the most expected condition is like this: the plastic yielding initiates and finishes in every layer simultaneously; in this way, the material strength can be most utilized. This expected condition can be realized by optimizing the ratio K, the radius R_{0}, and the thicknesses of each layer [29]. According to the inverse analysis, the general rule of the optimized result is that Young’s modulus needs to be increased monotonously from the inner to the outer layer [29].
4.3. Comparison with the Tresca Criterion
As mentioned in Section 3.2, parameter b can reflect the effect of the intermediate principal stress on the plastic yielding. Usually, the intermediate principal stress is helpful in increasing the material strength. To analyze this issue, the variations of the elastic bearing capacity with b are plotted in Figure 7. The analytical results support the above point: with the increase in value b, the elastic ultimate bearing capacity increases monotonously (see Figure 7). Specifically, as the value b increases from 0.0 to 1.0, the elastic ultimate bearing capacities of the three cases of K=25/41, 41/41, and 41/25 increase by 6.5%, 13.3%, and 13.3%, respectively. In other words, when the intermediate stress is considered, the elastic ultimate bearing capacity will increase, which is in line with the expectations.
To compare with the Tresca Criterion [29], taking b=0 and α=1, UST becomes the Tresca Criterion. The extending paths of the plastic zone and the radial displacement of the inner boundary based on the two criteria are plotted in Figures 8 and 9, respectively. Under a given confining pressure, the Trescabased plastic zone is larger than the USTbased plastic zone. The USTbased elastic and plastic bearing capacities are greater than the corresponding Trescabased bearing capacities. Specifically, for the singlelayered lining (i.e., K=41/41), the USTbased elastic bearing capacity is 17.360 MPa, 8.5% higher than the Trescabased value. For the twolayered lining (i.e., K=25/41), the corresponding increment is 4.1%, from 19.741 to 20.560 MPa. For a certain yielding criterion, the plastic bearing capacity keeps unchanged, which depends on the material strength. The USTbased plastic bearing capacities of the singlelayered and twolayered structures are the same, 26.399 MPa, 3.4% higher than the Trescabased plastic bearing capacity, 25.541 MPa. The two figures also show that, for any of the two criteria, the twolayered lining (K=25/41) has greater elastic bearing capacity than the singlelayered lining. The test data of the two cases K=25/41 and 41/41 are plotted in Figure 6 as well [29], which verifies the conclusion. The radial displacements by tests are smaller than the theoretical results, which is related to the methods of the tests.
5. Conclusion
The multilayered circular structures have been proved to have some advantages compared to the traditional singlelayered structure. An elastoplastic analysis details such structure; a twolayered circular lining is conducted based on the Unified Strength Theory, which can take the intermediate principal stress into consideration. The purpose of the work is to discuss the bearing capacities of different kinds of structures. The results show that the distributions of stresses in the twolayered lining are more reasonable compared to the traditional singlelayered lining. The twolayered Lining has greater elastic bearing capacity, but the plastic bearing capacity keeps unchanged for the same criterion. Considering the intermediate principal stress, the elastic bearing capacity increases correspondingly. Compared to the Trescabased results, the USTbased elastic bearing capacities are greater, and the plastic zones are smaller, which is due to the intermediate principal stress. The analytical conclusions are also verified by the results of the model tests.
Appendix
where
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research work was supported by the Natural Science Foundation of China (Grant No. 51704117, 11572126), the Fundamental Research Funds for the Central Universities (Grant No. 2017MS025), the Project of Shandong Province Higher Educational Science and Technology Program (Grant No. J15LG01) and the Natural Science Foundation of Shandong Province (Grant No. ZR2018MEE051), State key laboratory for Mine disaster prevention and control, and cultivation base cobuilt by province and ministry of Shandong university of science and technology (Grant No. MDPC201707). Their financial support is gratefully acknowledged.
References
 T. Zhao, W. Guo, Y. Tan, Y. Yin, L. Cai, and J. Pan, “Case studies of rock bursts under complicated geological conditions during multiseam mining at a depth of 800 m,” Rock Mechanics and Rock Engineering, vol. 51, no. 5, pp. 1539–1564, 2018. View at: Publisher Site  Google Scholar
 I. D. Peila and I. P. P. Oreste, “Axisymmetric analysis of ground reinforcing in tunnelling design,” Computers & Geosciences, vol. 17, no. 2, pp. 253–274, 1995. View at: Publisher Site  Google Scholar
 Q. Wu, Stress Analysis on ThickWalled Cylinder of Multilayer Concrete Material with Functionally Graded Characteristics, North China Electric Power University, 2009.
 A.Z. Lu, N. Zhang, Y.Q. Xu, and P. Cui, “Stressdisplacement solution for a lined vertical borehole due to nonaxisymmetric in situ stresses,” International Journal of Rock Mechanics and Mining Sciences, vol. 57, pp. 64–74, 2013. View at: Publisher Site  Google Scholar
 Q. Wu and A. Lu, “Stress analytical solution for plane problem of a thickwalled cylinder subjected to a type of nonuniform distributed pressures,” Engineering Mechanics, vol. 28, no. 6, pp. 6–10, 2011. View at: Google Scholar
 J. Feng, Plastic Limit Loading Analysis of thickwalled tube based on Twin Shear Unified Strength Theory, Xiangtan University, 2002.
 S. Xu, M. Yu, and X. Hu, “The stability analysis of circular tunnel based the twin shear unified strength theory,” Journal of China Coal Society, vol. 28, no. 5, pp. 522–526, 2003. View at: Google Scholar
 J. Zhao, Y. Zhang, H. Liao, and Z. Yin, “Unified limit solutions of thick wall cylinder and thick wall spherical shell with unified strength theory,” Chinese Journal of Applied Mechanics, vol. 17, no. 1, pp. 157–161, 2000. View at: Google Scholar
 M. Guowei, S. Iwasaki, Y. Miyamoto, and H. Deto, “Plastic limit analyses of circular plates with respect to unified yield criterion,” International Journal of Mechanical Sciences, vol. 40, no. 10, pp. 963–976, 1998. View at: Publisher Site  Google Scholar
 C. CarranzaTorres and C. Fairhurst, “The elastoplastic response of underground excavations in rock masses that satisfy the HoekBrown failure criterion,” International Journal of Rock Mechanics and Mining Sciences, vol. 36, no. 6, pp. 777–809, 1999. View at: Publisher Site  Google Scholar
 C. D. Martin, P. K. Kaiser, and D. R. McCreath, “HoekBrown parameters for predicting the depth of brittle failure around tunnels,” Canadian Geotechnical Journal, vol. 36, no. 1, pp. 136–151, 1999. View at: Publisher Site  Google Scholar
 D. Durban and M. Kubi, “General solution for the pressurized elastoplastic tube,” Journal of Applied Mechanics, vol. 59, no. 1, pp. 20–26, 1992. View at: Publisher Site  Google Scholar
 E. Detournay and C. Fairhurst, “Twodimensional elastoplastic analysis of a long, cylindrical cavity under nonhydrostatic loading,” International Journal of Rock Mechanics and Mining Sciences, vol. 24, no. 4, pp. 197–211, 1987. View at: Google Scholar
 H.W. Gao, L.S. He, and Y.Q. Zhang, “Unified plastic limit solution for thickwalled cylinder with size effect,” Journal of Zhejiang University (Engineering Science), vol. 45, no. 4, pp. 684–687, 2011. View at: Google Scholar
 S.Y. Leu, “Analytical and numerical investigation of strainhardening viscoplastic thickwalled cylinders under internal pressure by using sequential limit analysis,” Computer Methods Applied Mechanics and Engineering, vol. 196, no. 2528, pp. 2713–2722, 2007. View at: Publisher Site  Google Scholar
 Y.K. Lee and S. Pietruszczak, “A new numerical procedure for elastoplastic analysis of a circular opening excavated in a strainsoftening rock mass,” Tunnelling and Underground Space Technology, vol. 23, no. 5, pp. 588–599, 2008. View at: Publisher Site  Google Scholar
 S. Wang and S. Yin, “A closedform solution for a spherical cavity in the elasticbrittleplastic medium,” Tunnelling and Underground Space Technology, vol. 26, no. 1, pp. 236–241, 2011. View at: Publisher Site  Google Scholar
 S. Wang, X. Yin, H. Tang, and X. Ge, “A new approach for analyzing circular tunnel in strainsoftening rock masses,” International Journal of Rock Mechanics and Mining Sciences, vol. 47, no. 1, pp. 170–178, 2010. View at: Publisher Site  Google Scholar
 Q. Zhang, B.S. Jiang, S.L. Wang, X.R. Ge, and H.Q. Zhang, “Elastoplastic analysis of a circular opening in strainsoftening rock mass,” International Journal of Rock Mechanics and Mining Sciences, vol. 50, pp. 38–46, 2012. View at: Publisher Site  Google Scholar
 Q. Zhang, B. S. Jiang, X. S. Wu, H. Q. Zhang, and L. J. Han, “Elastoplastic coupling analysis of circular openings in elastobrittleplastic rock mass,” Theoretical and Applied Fracture Mechanics, vol. 60, no. 1, pp. 60–67, 2012. View at: Publisher Site  Google Scholar
 A. Z. Lu, G. S. Xu, F. Sun, and W. Q. Sun, “Elastoplastic analysis of a circular tunnel including the effect of the axial in situ stress,” International Journal of Rock Mechanics and Mining Sciences, vol. 47, no. 1, pp. 50–59, 2010. View at: Publisher Site  Google Scholar
 S. Wang, Z. Wu, M. Guo, and X. Ge, “Theoretical solutions of a circular tunnel with the influence of axial in situ stress in elasticbrittleplastic rock,” Tunnelling and Underground Space Technology, vol. 30, pp. 155–168, 2012. View at: Publisher Site  Google Scholar
 J. Li, G. Ma, and M. Yu, “Penetration analysis for geomaterial based on unified strength criterion,” International Journal of Impact Engineering, vol. 35, no. 10, pp. 1154–1163, 2008. View at: Publisher Site  Google Scholar
 C. Zhang, J. Zhao, Q. Zhang, and X. Hu, “A new closedform solution for circular openings modeled by the unified strength theory and radiusdependent Young’s modulus,” Computers and Geosciences, vol. 42, pp. 118–128, 2012. View at: Publisher Site  Google Scholar
 M.H. Yu, Y.W. Zan, J. Zhao, and M. Yoshimine, “A Unified Strength criterion for rock material,” International Journal of Rock Mechanics and Mining Sciences, vol. 39, no. 8, pp. 975–989, 2002. View at: Publisher Site  Google Scholar
 M. Yu, V. A. Kolupaev, Y. Li, and J. Li, “Advances in unified strength theory and its generalization,” Procedia Engineering, vol. 10, pp. 2508–2513, 2011. View at: Google Scholar
 M. Yu, S. Yang, S. C. Fan, and G. Ma, “Unified elastoplastic associated and nonassociated constitutive model and its engineering applications,” Computers and Structures, vol. 71, no. 6, pp. 627–636, 1999. View at: Google Scholar
 J. Feng, J. Zhang, P. Zhang, and L. Han, “Plastic limit load analyses of thick walled tube based on twin shear unified strength theory,” Acta Mechnica Solida Sinica, vol. 25, no. 2, pp. 208–212, 2004. View at: Google Scholar
 N. Zhang, A. Lu, C. C. Li et al., “Support performance of functionally graded concrete lining,” Construction and Building Materials, vol. 147, pp. 35–47, 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Ning Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.