Abstract

Swelling of soft rock, such as gypsum rock, is one of the major threats in tunnel engineering, causing structure damages such as floor heave and inward movement of sidewalls during construction and operation. It is of practical significance to study the swelling mechanical behavior of such rocks by tests. Swelling strain tests and swelling stress tests were performed by swelling test apparatus to study the variation of swelling strain with time and the swelling stress-strain relationship for gypsum rock samples, respectively. Three stages of the swelling strain on the time-strain curve of gypsum rock samples were noticed, which are defined as rapid swelling stage, slow swelling stage, and steady stage. And it was further found that the swelling strain caused in the slow swelling stage is of 76% of the total swelling strain. A constitutive model is proposed to describe the stress-strain relationship in swelling considering the swelling deformation and swelling pressure. The proposed model was verified using test data, which shows good agreements in describing the relationship between swelling strain and swelling stress, also in the conditions of maximum swelling strain and maximum swelling stress under lateral restraint situations.

1. Introduction

A growing number of road and railway tunnels in complex geological conditions are being constructed with the implementation of the China Western Development Strategy, which are considered as the important structure of transportation infrastructure. There are large amounts of formations containing swelling rocks, such as gypsum and anhydrite, widespread in Western China [1]. The swelling of rocks causes severe problems in tunnel engineering, such as floor heave, inward movement of sidewalls, and destruction of tunnel linings [2]. The swelling behavior of rocks was usually investigated in laboratory condition. Most laboratory experiments or theoretical analyses on the swelling characteristics of a shield tunnel in swelling rock are under certain hypotheses which have not been verified in actual situation. Currently, it is impossible to predict the swelling behavior of an actual geotechnical project due to our limited knowledge of the processes involved in the swelling of rock and of the geological, mineralogical, chemical, hydraulic, and mechanical controls of the swelling [3]. Historically, swelling rock studies have been largely concentrated on swelling mechanism [4,5] and factors affecting the swelling [69]. Several swelling tests are available on the variation of swelling deformation [10, 11] and swelling constitutive relationship [1219], which help geotechnical engineers to understand the swelling behavior of rock well.

The swelling of gypsum rocks is a major threat in tunnel engineering, causing serious damage to tunnels and producing high additional costs during tunnel construction and operation. Steiner et al. [20] proposed brittle failure of the anhydrite rock after tunnel excavation to be a major process leading to rock swelling. Liu et al. [21] studied the swelling characteristics of gypsum rock and obtained the relationships between swelling stress and water absorption also initial dry density. Gypsum rocks immersed by fresh water within different times were investigated according to softening tests by Ren et al. [22]. The variations of uniaxial compressive strength, shear strength, elastic modulus, and Poisson’s ratio of gypsum rocks with increasing immersion time were analyzed. The factors influencing the swelling characteristics of swelling rock were investigated [23, 24]; it was found that the density and the initial moisture content of rock significantly affected the swelling characteristics of rock.

Huder and Amberg [25] conducted experiments to investigate the swelling characteristics of swelling rock; a relation between normal stress and normal swelling strain was developed. Based on the experimental results of Huder and Amberg [25], Grob [26] proposed a mathematical model for relating final swelling strain to stress, the swelling strain decreases with logarithm of stress, and maximum swelling occurs when there is a minimum normal stress acting on the specimen. The one-dimensional mathematical model of Grob was expanded to three-dimensional constitutive law by Wittke-Gattemann and Wittke [27]. Einstein [28] examined the available swelling models in the literature and concluded that, although the swelling deformation due to chemical changes may be adequately predicted using available models, but time-dependent swelling models need to be developed from a fundamental knowledge of swelling mechanism. Miao et al. [29] proposed the humidity field theory which provided theories foundation for constitutive model of swelling rock. Based on the humidity stress field theory, Chen et al. [30] carried out swelling tests on anhydrite rock mass and proposed a time-dependent swelling constitutive model, and Yu et al. [31] proposed a physical equation to describe a tunnel’s mechanical responses in swelling rock. Most of the swelling constitutive models were developed based on the swelling tests and theoretical analysis, while the condition of rock tested in swelling tests is different from that in actual situation. Up to now, no common opinion is offered about the applicability of constitutive laws with which the swelling phenomena of rocks can be described. A swelling model considering the swelling behavior in engineering is needed to solve actual problems.

Swelling characteristics of gypsum rock varies as diagenesis environment and occurrence varies. Swelling strain tests were conducted to investigate the variation of swelling strain with time for gypsum rock. Based on the results of swelling strain test, swelling stress tests under conditions of controlling the swelling strain were carried out to study the swelling stress-strain relationship for gypsum rock. A formula describing the swelling potential of rock was proposed by introducing virtual swelling stress and virtual swelling strain. A swelling constitutive model was developed and verified using three sets of experimental data from the present research and literature.

2. Swelling Test

2.1. Materials

The gypsum rock blocks, taken from the surrounding rock of Dugongling tunnel in Shanxi Province, China, were used to prepare the samples of swelling tests. The gypsum rock blocks were transported to the laboratory and processed into the samples with 50 mm diameter and 36 mm height. The powder of natural gypsum rock passing sieve number 500 was used to test the mineralogical composition of gypsum rock by D8-ADVANCE type X-ray diffractometer. A typical XRD pattern is shown in Figure 1. Results of XRD test in Table 1 show that the gypsum rock is mainly composed of dihydrate gypsum, anhydrite, and quartz. The density and natural moisture content were tested according to the Chinese Soil Test Standard SL237. The measured density of gypsum rock ranged from 2240 kg·m−3 to 2370 kg·m−3. The natural moisture content of gypsum rock was 0.12%.

2.2. Test Procedure

The samples of gypsum rock were dried under 220°C condition for 48 hours using a drying oven. The dried samples were put in a sealed bag for sealing. Lateral restraint swelling equipment SCY-2 shown in Figure 2 was used to conduct the swelling strain tests and swelling stress tests. The primary components of SCY-2 are a water container, a rigid snug-fitting restraint ring, dial gage, and a loading device.

The swelling strain is the change in height of the sample divided by the initial height of the sample immersed in water and expressed in percentage. The procedure of the swelling strain test is as follows:(1)Apply some Vaseline lotion around the inside of a rigid snug-fitting ring and place the rock sample in the ring. The ring containing the sample is placed in the water container.(2)Place the porous disks at the top and bottom of the sample.(3)Assemble the water container and dial gage. Check the dial gage and make sure it is precise. The top of the sample is not subjected to any axial load.(4)Add enough water into the water container so that swelling can start. The water level should be higher than the surface of upper porous disk.(5)Record the axial swelling movement of the sample as a function of time using the dial gage. The swelling test is finished when the axial swelling movement reaches an equilibrium value, that strain is the maximum swelling strain ε0.

Swelling stress is the ratio of the force which can inhibit the swelling strain caused by the water seeped into the sample to the sectional area of the sample. Swelling stress tests were performed to measure the swelling stress using balanced pressure method according to the Chinese Soil Test Standard SL237. Four sets of swelling stress tests were conducted to measure the swelling stress by controlling the axial strain equal to 0, 20, 40, and 60 percent of the maximum swelling strain measured in the swelling strain test. The procedure of the swelling stress tests is as follows:(1)Apply some Vaseline lotion around the inside of a rigid snug-fitting ring and place the rock sample in the ring. The ring containing the sample is placed in the water container.(2)Place the porous disks at the top and bottom of the sample. The pressure sensor is placed at the top of the upper porous disk.(3)Assemble the water container, dial gage, and loading device. Check the dial gage and make sure it is precise.(4)Add enough water into the water container so that swelling can start. The water level should be higher than the surface of upper porous disk.(5)Record the swelling movement by the dial gage. Apply axial load to maintain the swelling movement constant when the swelling strain reaches the required value (0%, 20%, 40%, and 60% maximum swelling strain). The swelling test is finished when the axial swelling stress reaches an equilibrium value, that stress is the swelling stress under the induced swelling strain condition.

The swelling strain is maintained as 0% maximum swelling strain during the swelling stress tests; when the axial stress reaches an equilibrium value, that stress is the maximum swelling stress σ0.

3. Results of the Swelling Test

3.1. Swelling Strain

Figure 3 presents a plot showing the variation of swelling strain with time. It can be seen from Figure 3(a) that the procedure of swelling for gypsum rock lasted for about 300 minutes. The maximum swelling strain was equal to 4.397%. Regression works were conducted based on maximizing the coefficient of determination R2. The regression result shows that there is a positive exponential relationship between swelling strain and time:where ε is the swelling stain and t is the swelling time.

The slope of the dash curve (Figure 3(a)) is the derivative of equation (1), which increases first and then decreases with increasing the time as shown in Figure 4. A point of the dash curve represents the maximum derivative of equation (1), which corresponds to the swelling time of 18.1 minutes. The swelling strain of the point corresponding to the swelling time of 18.1 minutes was about 0.56%; this point on the swelling strain versus time curve was named as A. Once the swelling time exceeds 18.1 minutes, the slope of dash curve shows a decreasing trend as the time increases. The data point at which the swelling strain is 90% of the maximum swelling strain was selected as point B shown in Figure 3(a). The data point on the swelling strain versus time curve corresponding to the final time was named as point C. According to the three points, A, B, and C, the swelling strain versus time curve shown in Figure 3(a) can be divided into three parts: rapid swelling stage (path 0A), slow swelling stage (path AB), and stable stage (path BC).

The rapid swelling stage lasted for about 18 minutes. For rapid swelling stage, the swelling deformation was too small to be measured at the first three minutes. The voids in the gypsum rock were filled with the swelling deformations of rock that occurred at the first three minutes. As the filling of swelling deformations in the rock voids, the increment of swelling strain increases rapidly at the rapid swelling stage (Figure 3(b)). The slow swelling stage lasted for about 115 minutes, and the swelling strain (ε) was 3.96%. The swelling strain generated in the slow swelling stage is 76% of the maximum swelling strain. The path BC in Figure 3(a) shows that the increment of swelling deformation gradually tends to zero with an increase in time, and the swelling strain ends up reaching an equilibrium value called the maximum swelling strain ε0.

3.2. Swelling Stress

Results of swelling stress under conditions of controlling the swelling strain of 0%, 20%, 40%, and 60% maximum swelling strain are shown in Table 2. The measured maximum swelling stress was 3340 kPa for the gypsum rock tested in the present research.

4. Swelling Constitutive Model

4.1. Development of the Swelling Constitutive Model

For a certain type of swelling rock, the swelling characteristics are closely related to the intrinsic rock property without considering the environmental factors. The swelling behavior of rock is externally expressed in the form of swelling pressure and swelling deformation which are parameterized by swelling stress and swelling strain, respectively [32]. Different from the rock stress-strain relationship measured under external load, there is a negative correlation between swelling stress and swelling strain when considering the swelling behavior of rock. The swelling stress shows a decreasing trend with an increase in swelling strain. For the condition of maximum swelling stress, the swelling stress is the largest and the swelling strain is 0. For the condition of maximum swelling strain, the swelling stress is 0 and the swelling strain is the largest.

Assumptions were made to propose a swelling constitutive model considering the swelling behavior of rock as follows:(1)The swelling potential Qs of a certain swelling rock only depends on its intrinsic property, which can be expressed by swelling stress or swelling strain.(2)Swelling pressure and swelling deformation are the external form of swelling behavior. The swelling potential can be decomposed into swelling potential expressed by swelling pressure Qσ and swelling potential expressed by swelling deformation Qε, which can be parameterized by swelling stress or swelling strain. The expression of swelling potential of rock is as follows:(3)The virtual swelling moduli and were introduced into the expressions of virtual swelling stress (σε, kPa) and virtual swelling strain (εσ), respectively:where ε is the actual swelling strain, and σ (kPa) is the actual swelling stress.

In this research, the swelling potential was parameterized by swelling stress. The swelling stress-strain relationship was analyzed using the experimental data from the swelling tests under lateral restraint condition. It is assumed that there are linear relationships between swelling stress and swelling strain, so that equations (3) and (4) become

Qσ and Qε are expressed using the actual swelling stress and virtual swelling stress, respectively, which can be written as follows:

Substituting equations (6) and (7) into (2) leads to

Equation (8) is a general formula, which can represent the swelling behavior of a certain type of swelling rock under different conditions. The swelling rock under the condition of maximum swelling stress and maximum swelling strain are the two special conditions in equation (8).

The maximum swelling stress σ0 can be measured using the swelling stress test under lateral restraint condition. The axial swelling strain of a sample immersed in water is kept 0 by applying axial load; the maximum loading stress is the maximum swelling stress. For the condition of maximum swelling stress, equation (8) can be expressed as follows:

The maximum swelling strain ε0 is the final swelling strain measured from the swelling strain test under lateral restraint condition. The top of the rock sample is not subjected to any axial load. For the condition of maximum swelling strain, equation (8) can be written as

For a certain swelling rock, the swelling potential Qs is constant under different swelling conditions. Equations (9) and (10) can be rewritten as

Substituting equations (9) and (11) into (8), one can derive

Equation (12) can be rewritten as

Equation (13) is the swelling constitutive model developed in present research.

4.2. Test of the Proposed Model

The proposed constitutive model (equation (13)) was tested using the experimental data presented in this paper and the collected data from Chang et al. [14] and Liu et al. [11]. In the present research, the measured maximum swelling stress (σ0) and maximum swelling strain (ε0) were 3340 kPa and 4.397%, respectively. Figure 5 compares the measured swelling stress and those calculated using equation (13) for gypsum rock. Equation (13) has a R2 of 0.993 for the data in Table 2, which indicates that the proposed swelling constitutive model reproduces well the swelling stress-strain relationship for gypsum rock tested in the present research.

The experiments of Chang et al. [14] were conducted on red sandstone taken from Zhuzhou in China to measure the swelling strain by loading the specimen under one certain axial pressure in the condition of lateral restraint. The pressure acted on the rock sample can be considered as the swelling stress. The measured swelling strain is shown in Table 3. The maximum swelling stress (σ0) and maximum swelling strain (ε0) were not measured in the experiments of Chang et al. [14]. Linear interpolation method was used to determine σ0 and ε0 according to two sets of data points on line No. 1 and on another line shown in Table 3. For example, the σ0 on line No. 2 was calculated using the data of σ and ε on line No. 1 and on line 2, the same as ε0. The average of the calculated σ0 was considered as the parameter of the proposed formula, the same as ε0. The maximum swelling stress (σ0) and maximum swelling strain (ε0) of red sandstone were estimated as 163.10 kPa and 0.923%, respectively.

Figure 6 compares the predictions and measurements for the red sandstone tested by Chang et al. [14]. The R2 between the measured swelling stress and the calculated swelling stress is 0.975, which reveals that the proposed swelling model can be used to describe the swelling stress-strain relationship for red sandstone.

Liu et al. [11] studied the swelling properties of remolded rock in Middle Route Project of South-to-North Water Diversion using swelling strain tests and swelling stress tests. The swelling stress of rock was measured by remaining the swelling strain constant at 0%, 1%, 2%, and 3% for swelling rock with different initial water content , 14%, 17%, and 20% (Table 4). The maximum swelling stress of rock with initial water content of 14%, 17%, and 20% was 1044, 1014, and 924 kPa, respectively. The maximum swelling strain was not measured by Liu et al. [11], which was calculated based on the experimental data shown in Table 4 using linear interpolation method mentioned in previous paragraph. The calculated maximum swelling strain for rock with initial water content of 14%, 17%, and 20% was 6.93%, 6.72%, and 6.67%, respectively. The variation of calculated maximum swelling strain with initial water content indicates that an increase in the initial water content leads to a decrease in swelling strain in the same external load state. This is consistent with the finding of Wang et al. [23].

Figure 7 shows the calculated and measured swelling strain for swelling rock tested by Liu et al. [11]. The R2 between the measured and calculated swelling strain for the three sets of tests is larger than 0.95, which reveals that the proposed swelling constitutive model performs well in describing the swelling stress-strain relationship for remolded swelling rock.

5. Discussion

A constitutive model (equation (14)) proposed based on the experimental results of Huder and Amberg [25] was widely used in describing the swelling stress-strain relationship. The difference between the proposed model in present research and the Huder–Amberg model was discussed. The Huder–Amberg model is written as follows:where K is the axial swelling strain when the swelling stress σ is 0.1 MPa.

The proposed constitutive model was compared with the Huder–Amberg model; it can be found that the expressions of the two models are similar. There is a linear relationship between swelling stress and swelling strain in equation (13), and the swelling strain shows a linear relationship with logarithm swelling stress in the Huder–Amberg model.

Some disadvantages of the Huder–Amberg model were found based on the analysis of equation (14). Variation of axial swelling strain varies with axial swelling stress when the maximum swelling stress σ0 is larger than 10 MPa. For example, when the axial swelling stress σ is larger than 1 MPa, lgσ is larger than 0, the value of 1 − lgσ/lgσ0 is less than 1. While, in the case of σ less than 1 MPa, the value of 1 − lgσ/lgσ0 is larger than 1.

In the case of 0 < σ < 0.1 MPa and σ0 > 1 MPa, the calculated axial swelling strain is larger than the value of K, which is not in accordance with the definition of K. It is difficult to describe the stress-strain relationship in the case of σ ≈ 0. The calculated swelling strain will be infinite if swelling stress is approximately 0 MPa, which is not consistent with the results of swelling test. The Huder–Amberg model is shown to be an inappropriate constitutive model to describe the swelling stress-strain relationship in the case of 0 < σ < 0.1 MPa. The analyzed disadvantages of the Huder–Amberg model did not exist in the proposed model (equation (13)).

The constitutive model (equation (13)) was proposed depending on the relationship between virtual swelling stress and swelling strain shown in equation (3) and the relationship between virtual swelling strain and swelling stress shown in equation (4). A linear assumption was introduced in equation (3) and equation (4), which leads to a linear relation between swelling strain and swelling stress. Different types of swelling rock show different swelling stress-strain relationships; a linear stress-strain relationship will not be suitable for all types of swelling rock. In order to optimize the proposed constitutive model, additional swelling tests and research need to be conducted.

6. Conclusions

Swelling strain and swelling stress are the main test results to evaluate the characteristics of swelling rock. In order to study the properties of swelling rock, the swelling strain test and swelling stress tests were performed on gypsum rock samples drilled from a tunnel jobsite in Shanxi province, China. The study on the swelling constitutive model is an important issue for swelling rock. Based on the experimental data of swelling stress tests, a constitutive law describing the relationship between swelling strain and swelling stress was proposed by introducing the virtual swelling strain and virtual swelling stress. As a result of this study, the following conclusions were derived:(1)The measured maximum swelling strain of gypsum rock sample in the swelling strain test was 4.397%. Three stages of the swelling strain of gypsum rock were defined by the tests results, namely, rapid swelling stage, slow swelling stage, and steady stage. The slow swelling stage lasted for about 115 minutes in which the swelling strain accounted for 76% of the total swelling strain.(2)The measured maximum swelling stress of gypsum rock was 3340 kPa. There was a negative linear relationship between swelling stress and swelling strain for the gypsum rock tested in the present research.(3)The constitutive model was verified using the experimental data from the present research and the collected data from other literatures. The proposed constitutive model performs well in describing the variation of swelling strain with swelling stress for different types of swelling rock. The developed model has good agreements with test data in the condition of maximum swelling stress and maximum swelling strain.

In order to optimize the proposed constitutive model, additional swelling tests and research needs would be further conducted.

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 regarding the publication of this paper.

Acknowledgments

This work was financially supported by the Special Science and Technology Innovation fund of Research Institute of Highway Ministry of Transport (Grant no. 2019-C505), the Central Public-Interest Scientific Institution Basal Research Fund (Grant no. 2020-9032), and the Science and Technology Project of Department of Transportation of Shanxi Province (Grant no. 2017-1-6).