Abstract
Weak and hard inhomogeneous rock formations are typically encountered during tunnel excavations. The physical and mechanical properties and geological conditions of these rock formations vary significantly; thus, it is crucial to investigate the mechanical characteristics of deep bedded composite rock formations. Three-dimensional (3D) scanning and 3D printing were used to prepare composite rock specimens to simulate natural rock laminae. Triaxial compression tests were conducted to determine the influence of the bedding angle, rock composition, and confining pressure on the mechanical properties of the composite rock specimens. The anisotropic strength characteristics and the damage patterns of the composite rock specimens were analyzed under different confining pressures, and the failure mechanism during triaxial loading was revealed. The results show that the damage of the composite rock specimens with a bedding structure depends on the bedding dip angle and the rock formation. The stress-strain curves and peak strengths of the composite rock specimens have anisotropic characteristics corresponding to their failure modes. As the bedding dip angle increases, the peak strength of the three groups of specimens first decreases and then increases under different confining pressure levels. The compressive strength has a nonlinear relationship with the confining pressure, and the difference between the compressive strengths of specimens with different inclination angles decreases as the confining pressure increases. The Hoek–Brown strength criterion is a good predictor of the nonlinear increase in peak strength of the composite rock specimens under different confining pressures. The specimen with a = 60°dip angle shows the most significant increase in the strength difference with increasing confining pressure. The results can be used as a reference for testing and analyzing the anisotropic mechanical properties of bedded rock masses.
1. Introduction
Layered rocks are extensively encountered in mining, civil engineering, and geological engineering projects [1–4]. The rocks’ bedding planes are weak links that may result in geological disasters, such as landslides, collapses, and water bursts. Figure 1 shows a weak and hard inhomogeneous stratum encountered during tunneling in a practical project. Therefore, the mechanical performance of layered rocks has become a crucial issue in recent decades. An accurate determination of the deformation characteristics and failure mechanics requires an in-depth understanding of the influence of the bedding plane orientation and nonuniform deformation of layers on the strength characteristics of composite rock formations.

(a)

(b)
Numerous studies have been conducted to investigate the mechanical behavior of layered rocks. Tavallali and Vervoort [5] noted that the crack pattern of layered rocks under compression consisted of sliding along the bedding plane and subvertical cracks. Debecker and Vervoort [6] used acoustic emission (AE) data to describe the crack evolution of layered rocks. Many investigations have been conducted to predict strength changes in layered rock as a result of the layer orientation. The weakness in the in-lane directions and the associated theory on variable cohesive strength have been used to describe the influence of weak bedding on the strength of layered rocks [7, 8]. This model was further modified to accommodate changes in the friction angle within the layered rocks [9]. Nova [10] proposed a generalized failure criterion to describe the failure strength of layered rocks in three-dimensional (3D) stress states. It is difficult to obtain adequate layered rock specimens with uniform properties and an ideal layered structure because of the high variability of natural rock arising from the formation processes and geological environments. Tien and Tsao [11] proposed a technique for preparing layered rocks and conducted experiments and analyses to determine the rocks’ failure modes and specimen strengths. Chang and Deng [12] investigated the failure mechanisms of unreinforced and reinforced layered rocks with various layer dip angles using AE tests to obtain the failure modes and load-displacement responses.
These studies significantly advanced the understanding of the mechanical behavior of layered rocks. It should be noted that the macro- and microstructure characteristics of the interface between the rock layers were not considered in these studies. Therefore, our understanding of the fracture behavior of rock layers is limited, and the prediction of the mechanical performances of layered rock may be inaccurate. 3D printing (3DP) technology has developed rapidly in recent years and can be used to create structurally complex 3D objects. Several scholars [13–16] have utilized 3DP in their research on rock mechanics. Therefore, this study uses 3D scanning of natural rock laminae, mortar casting, and 3DP to produce rock specimens with bedded rock laminae to study the mechanical response of the composite rock specimens under triaxial loading conditions. The stress-strain curves, deformation, and failure modes of the composite rock specimens with different bedding dip angles, different rock formations, and different confining pressures are analyzed. The Mohr–Coulomb criterion and the Hoek–Brown criterion are used to fit the test results and analyze the effects of the bedding dip angle and confining pressure on the anisotropic characteristics of the bedded rock mass.
2. Specimen Preparation
The hardness of the rocks in the engineering rock classification standard [17] is used as the basis for sample preparation. The initial strengths of the soft and hard rocks in the composite rock formations were approximately 10 MPa, 20 MPa, and 30 MPa. The physical and mechanical parameters of the soft and hard rocks are listed in Table 1. The combinations are Group AA, Group AB, and Group AC. It should be noted that soft and hard rocks are relative concepts in this paper and are used only to distinguish between materials of different strengths.
The bedding plane of natural rock laminae was obtained by splitting a kilim rock, performing 3D scanning of the layers, and symmetrically combining the layers using 3DP to obtain a natural bedding plane with a dimension of 450 mm × 250 mm × 10 mm. A rectangular plastic mold with internal dimensions of 450 mm × 250 mm × 200 mm was used for preparing the composite specimen. During the pouring process, the mold was placed on a vibrating machine to compact the mix. The printed specimen was pressed onto the mortar to provide the morphological characteristics of a natural bedding plane. After 5 h of drying, the surface was sufficiently imprinted, and the bedding plane was removed before continuing with the next layer. The raw materials used in the tests and the pouring, drying, and processing methods were consistent to ensure uniform bedded composite rock specimens. Black ink was added to the water when the hard rock layers were fabricated to differentiate between soft and hard layers in the same specimen. The process is illustrated in Figure 2.

The finished layered rock specimens were immersed in water for 3 d to reach a certain initial strength. The specimens were drilled at 15 intervals using a 50 mm diameter drill and the drilling method shown in Figure 3; the angle between the loading direction and the normal direction of the laminae is defined as . The drilled specimens were then fabricated into standard cylindrical specimens of Ф50 mm × L100 mm (Ф is diameter of section, L is the length of the specimen), immersed in water, cured to final strength, and subjected to conventional triaxial tests. The completed specimens are shown in Figure 4.

(a)

(b)

3. Test System and Method
The conventional triaxial compression test was carried out on the RMT-150B testing machine developed by the Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, as shown in Figure 5. A total of four confining pressures (0, 2, 4, and 6 MPa) were used during the test. Seven specimens with different angles were tested under each confining pressure. The confining pressure was loaded in a force-controlled manner, starting with a hydrostatic state () and loading to a predetermined confining pressure at a loading rate of 2 MPa/min. Subsequently, the confining pressure was kept constant, and the axial force was applied in a displacement-controlled mode at a displacement rate of 0.005 mm/s. The axial force was loaded until the specimens were damaged.

4. Results
4.1. Stress-Strain Curve Analysis
Figure 6 shows the typical stress-strain curves for the group AA specimens ( = 0°–90°) under different confining pressure conditions. It is observed that the specimens undergo an elastic phase, a plastic phase, and a postpeak residual phase during the loading process. Under the effect of the confining pressure, the specimens with the 0°, 15°, 30°, 45°, 75°, and 90° angles exhibit tensile damage with a certain residual strength. After the specimens with the 60° angle reach their peak, the stress rapidly drops, indicating brittle characteristics due to the shear sliding of the specimen along the layered surface during the loading process. This behavior is also observed during the failure.

(a)

(b)

(c)

(d)

(e)

(f)

(g)
Figure 7 shows the typical stress-strain curves for the group AB specimens ( = 0°–90°) under different confining pressure conditions. The ductility characteristics of this group of rocks become more pronounced with increasing confining pressure due to the presence of the B-layer soft rock. The specimens with = 0°–45° angles only show peaks at low confining pressure (2 MPa), and at high confining pressure (4 MPa, 6 MPa), there is no peak, and the curves show ductility characteristics. The specimens with = 60° show peaks at 2 MPa and 4 MPa due to slip damage in the bedding direction and no peaks at 6 MPa.

(a)

(b)

(c)

(d)

(e)

(f)

(g)
Figure 8 shows the typical stress-strain curves for Group AC specimens ( = 0° to 90°) at different confining pressures. There are no peaks in the stress-strain curves for the soft rock at any confining pressures; thus, the strength is the axial compressive strength at 0.02% strain. Due to the lower strength of the C-layer soft rock, the specimens ( = 0° to 60°) exhibit ductility characteristics at different confining pressures, with no peaks. For the specimens = 75° and 90°, the more brittle hard rock layer A played a supporting role due to the boundary conditions. The soft and hard rock layers deformed almost simultaneously in the axial direction during loading, and shear failure occurred; thus, these two groups of specimens exhibited brittle failure characteristics.

(a)

(b)

(c)

(d)

(e)

(f)

(g)
4.2. Effect of Bedding Dip Angle on Strength Characteristics
Figure 9 shows the peak strengths of the three groups of specimens with different dip angles and confining pressures. As the dip angle increases, the peak strengths of the three groups of specimens at different confining pressures show an asymmetric “U” pattern with an initial decrease followed by an increase. In Group AA, the peak strengths of the specimens with = 0° are similar to those with = 90° at each confining pressure, and the peak strength of specimens with a = 60°angle is the lowest, indicating that the strengths of the layered specimens have significant anisotropic characteristics. The patterns are similar for all three groups, but due to the presence of soft rock layers, the peak strengths of the specimens with low dip angles ( = 0°–60°) in the two groups are the same under low confining pressure, and there are no anisotropic characteristics. The peak strength of the specimens with large angles ( = 75°–90°) is significantly higher than that of the specimens with low angles, especially in the AC group. The reason is that the specimen is fixed at the lower end and pressurized at the upper end in the test, and the boundary conditions of the rock specimens prevent it from completely sliding along the bedding plane to produce slip damage, which is clearly visible in the damage pattern. The specimens that exhibited splitting as the dominant damage mode had much higher peak stresses than those with complete shear failure and shear failure.

(a)

(b)

(c)
The anisotropy ratio was used to evaluate the degree of anisotropy of the bedded rock mass for the three sets of peak test results [18]. The degree of anisotropy of the laminated rock mass is described by the ratio of , which is the maximum compressive strength at a loading angle of 90° to , that is, the minimum compressive strength at a loading angle between 0° and 90°:
The ranges of the anisotropy ratios and the anisotropy classes are listed in Table 2, and the relationship between the anisotropy ratio and the confining pressure is shown in Figure 10. As can be seen from Figure 10, of the specimens in the AA, AB, and AC groups under low confining pressure (2 MPa) are 2.36, 2.2, and 2.02, respectively. All specimens had moderate anisotropy. At high confining pressures (4–6 MPa), for the AA group were 2.21 and 2.1, respectively, whereas for the AB and AC groups decreased to below 2.0 at high confining pressures, falling into the low anisotropy category. The level of anisotropy of all three groups decreased with an increase in the confining pressure. As the confining pressure increased, the specimens containing laminate structures exhibited consolidation, restoring the mechanical continuity of the material; thus, the stratification effect decreased.

4.3. Failure Mode
Figures 11(a)–11(c) show the failure characteristics of the layered composite rocks of Group AA, Group AB, and Group AC, respectively, after the compression tests under different confining pressures. Tien et al. [4] divided the failure modes of composite rock formations under uniaxial loading conditions into two major categories: slip failure in the stratification direction (Type II) and nonslip failure in the stratification direction. In this paper, the nonslip failure mode is further subdivided into five subcategories, namely, axial splitting (Type I), diagonal shear failure in the bedding plane (Type III), slip and axial splitting along the bedding plane (Type IV), Y-shaped splitting (Type V), and shearing in weak layers (Type VI). The failure modes are shown in Figure 11(d).

(a)

(b)

(c)

(d)
The failure modes of the specimens in Group AA with a low dip angle ( = 0°–30°) were similar (Figure 11(a)). At lower confining pressures ( = 0–2 MPa), axial splitting damage occurred (Type I), and the specimens showed longitudinal cracks along the cylindrical surface with significant damage. As the confining pressure increased ( = 4, 6 MPa), the failure mode changed to shear failure with an oblique intersection in the bedding plane (Type III). The = 45° specimens subjected to low confining pressure conditions ( = 0–2 MPa) exhibited slip damage along the bedding plane (Type II). When the confining pressure was increased ( = 4, 6 MPa), the anisotropy of the specimens weakened with increasing confining pressure, and coupled slip and axial splitting failures occurred along the bedding plane (Type IV). The reason is that the presence of the confining pressure restricted the bedding slip in the stratification direction, and splitting failure occurred in the rock mass. The = 60°–75° specimens had a larger dip angle, and the failure mode under different confining pressures was slip failure in the stratification direction. For the specimens with = 75°, the damage was accompanied by a few splitting cracks along the slip of the bedding plane due to the influence of the boundary conditions. The = 90° specimens exhibited axial splitting failure (Type I) at lower confining pressures ( = 0–2 MPa), whereas at high confining pressures ( = 4, 6 MPa), the rupture mode was diagonal shear failure in the bedding plane direction (Type III).
Figure 11(b) shows the failure mode of the composite rock mass of Group AB under different confining pressures. The specimens with lower bedding dip angles ( = 0°–30°) showed Y-shaped splitting failure (Type V) under different confining pressures. The strengths of the two rock layers A and B in the specimens were different, resulting in cracks mainly in the lower strength layer B. The specimens with = 5° exhibited slip failure in the stratification direction (Type II) under no confining pressure ( = 0 MPa). When confining pressure was applied ( = 2–6 MPa), the specimens exhibited coupled slip and axial splitting along the bedding plane (Type IV), and cracks only appeared in the weaker B formations. The specimens with = 60° to 75° exhibited slip damage along the bedding plane under different confining pressures. The = 90° specimens all showed axial splitting failure (Type I) at lower confining pressures ( = 0–2 MPa), whereas at confining pressures of = 4–6 MPa, the primary fracture mode was diagonal shear failure in the bedding plane (Type III).
Figure 11(c) shows the failure mode of the composite rock masses of the AC group under different confining pressures. The specimens with a lower bedding angle ( = 0°–30°) showed Y-shaped splitting failure (Type V) under no confining pressure. After the confining pressure was applied ( = 2–6 MPa), unlike the first two groups, the damage to the specimens only occurred in the C formation, whereas no damage was observed in the A formation, indicating shear damage in the weak layer (Type VI). The reason is that the strength of the C formation is much less than that of the A formation; thus, the former is more likely to be damaged in the weak layer under confining pressure. The = 45°–75° specimens under different confining pressures exhibited slip failure in the stratification direction (Type II). The = 90° specimens showed axial splitting failure (Type I) under uniaxial loading, whereas under confining pressure ( = 2–6 MPa), the dominant rupture mode was diagonal shear failure in the bedding plane (Type III).
The experimental results show that the bedding angle and lithological composition are the dominant factors influencing the compressional fracture behavior of the layered rock masses under triaxial compression.
5. Discussion of the Strength Characteristics
The linear Mohr–Coulomb strength criterion [19] is the most widely used strength criterion in rock mechanics. It is defined aswhere is the shear strength of the rock, is the positive stress in the shear plane, and and are the cohesion and angle of the internal friction of the rock material, respectively.
The Mohr–Coulomb criterion can be rewritten in the form of principal stresses as
Equation (3) can be used to describe the strength of different groups of specimens under triaxial stress conditions as follows:where and are the maximum and minimum principal stresses at the failure of the rock mass, respectively, is the uniaxial compressive strength of the rock specimen, and is a coefficient.
The characteristic strength values of the three groups of layered specimens in Figure 8 at different confining pressures were fitted using the Mohr–Coulomb strength criterion. The results are shown in Figure 12. The calculated coefficients of the Mohr–Coulomb criterion for each group of rocks are listed in Table 3. The correlation coefficients ranged from 0.85 to 0.98.

(a)

(b)

(c)
Although the Mohr–Coulomb criterion is straightforward and has a clear physical meaning, it has limited application because it is based on the assumption that the cohesion and internal friction angle of the rock are constant during the test. Due to the anisotropic nature of the composite formation specimens and the nonlinear behavior of the confining pressure and peak strength, the Hoek–Brown strength criterion [20] is used in a nonlinear regression of the triaxial peak strengths for the three sets of specimens. This criterion was first proposed by Hoek and Brown in 1980 as an empirical strength criterion for intact rocks based on the results of triaxial tests on hundreds of sets of rocks. A generalized Hoek–Brown strength criterion was later proposed based on a large number of field rock tests [21]. It is defined inwhere and are the maximum and minimum principal stresses in a rock sample at failure, is the uniaxial compressive strengths of the rock sample, and , , , and are constants related to the properties of the rock mass. Equation (5) can also be expressed as
By setting and , we obtain
The squared correlation coefficient is
When the rock sample is an intact rock, we set = , = 0.5, and = 1. The expression is
Based on this equation, the strength of the specimens with different inclination angles can be described aswhere , , , and are the uniaxial compressive strength, maximum principal stress, minimum principal stress, and parameters from the Hoek–Brown criterion for specimens with different inclination angles, respectively. The corresponding and of the criterion can be derived for each angle of the specimens. The calculated values are listed in Table 4, and the fitting results are shown in Figure 13. ranges from 0.96 to 0.99, indicating higher accuracy than that obtained from the Mohr–Coulomb strength criterion (0.85 to 0.98). A comparison of the values of the fitting results based on the two criteria is shown in Figure 14. The results indicate that the Hoek–Brown strength criterion has higher accuracy for predicting the failure strength characteristics of the composite rock specimens under different confining pressures.

(a)

(b)

(c)

Both sides of (10) are divided by to obtain a dimensionless parameter independent of the uniaxial compressive strength:
The Hoek–Brown strength curve is standardized by substituting the parameters of the specimen with different inclination angles (Table 2) into (11), as shown in Figure 15.

(a)

(b)

(c)
With the increasing of confining pressure, the closer the curve is, the lower the strength anisotropy of the specimen is, and conversely, the far the curve is, the more pronounced the anisotropy is. The specimens with = 60°in all three groups show the most significant increase in the strength difference with an increase in the confining pressure. It can be inferred that the nonlinear variation of the peak strength with the confining pressure becomes more pronounced as the confining pressure increases.
6. Conclusion
We used 3D scanning and 3DP to prepare composite specimens with a bedding structure and three different strengths and conducted triaxial compression tests to investigate the strength anisotropy and deformation of the composite rocks. We analyzed the failure mode and failure mechanism of the composite specimens. The conclusions are as follows.(1)The bedding angle and composition substantially influenced the failure mode of the rock body. When the lithology of the two layers was similar, the rock body exhibited axial splitting (Type I), diagonal shear failure in the bedding plane (Type III), and slip failure in the stratification direction (Type II) for bedding angles of [0°–90°]. When the lithology of the two layers was dissimilar, slip failure and axial splitting occurred along the bedding plane (Type IV), and Y-shaped splitting (Type V) and shearing in the weak layers (Type VI) were observed due to the weaker rock mass. The failure mode and mechanical properties depended on the rock mass, the rock formation, and the stress state.(2)The stress-strain curves and peak strengths of the composite specimens with layered structures showed anisotropic characteristics corresponding to their failure modes. As the bedding angle increased, the peak strength of the three groups of specimens showed an asymmetric “U”-shaped curve that first decreased and then increased under different confining pressure levels. The peak strength of the specimens was linearly related to the confining pressure, and the difference in the compressive strength of the specimens with different bedding angles decreased as the confining pressure level increased.(3)The anisotropy of the rock strength gradually decreased with increasing confining pressure, evidenced by (a) a decrease in the anisotropy coefficient with increasing confining pressure and (b) a change in failure mode from slips along the joints to block damage with increasing confining pressure, as in the case of the specimens with a bedding angle of = 45°.(4)The Hoek–Brown strength criterion provided better performance than the Mohr–Coulomb strength criterion for predicting the nonlinear failure strength characteristics of the composite specimens at different confining pressures. The specimens with = 60° in all three groups showed the most significant increase in strength with increasing confining pressure.
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 they have no conflicts of interest.
Acknowledgments
The work described in this paper was fully supported by National Natural Science Foundation of China (51674100).