Research Article  Open Access
Modification of Vertical Earth Pressure Formulas for High Fill CutandCover Tunnels Using Experimental and Numerical Methods
Abstract
The highfilled cutandcover tunnel (HFCCT) is a solution to reclaim more useable lands due to the unique landforms of Loess Plateau in northwestern China. Because of the ultrahigh backfill, the estimation of vertical earth pressure will significantly affect the design and safety of the cutandcover tunnel (CCT). The current methods for estimating the vertical earth pressure are either to overestimate or underestimate the vertical earth pressure on the top of HFCCT. To more precisely estimate the vertical earth pressure distribution, the vertical earth pressure based on the soil column pressure, γh (h: the height of backfill above the CCT), needs to be properly modified. Considering different influential factors, four corresponding coefficients are proposed: : crosssectional shape of CCT effect, : stiffness of backfill effect, : width of CCT effect, and : coupling effect of slope angle, θ, and the ratio of B/D. It is found that has little influence; the and reduce and amplifies the γh. The corresponding general forms for these coefficients are determined based on finite element analysis results. A general equation for estimating vertical earth pressure for the HFCCT including these four coefficients is proposed. Meanwhile, this general form is verified by the numerical analysis results and experimental results for different cases. Therefore, this proposed equation is applicable to estimate the vertical earth pressure for existing or newly designed HFCCT. Furthermore, this proposed method can significantly reduce the computational work in engineering analysis.
1. Introduction
With the rapid development of economy and the massive constructions of highspeed railway, the demand of useable land has been increased tremendously in northwestern China. The Northwest Loess Plateau is mountainous and deep valley. In order to create more useable lands in this area, reclamation above cutandcover tunnel (CCT) is a good strategy. Due to these unique landforms in northwestern China, the required backfill above the CCT is usually massive and high. The highfilled (HF) materials will significantly increase the overburden pressure on this CCT. Therefore, to more accurately estimate the magnitude of the resulting increase in vertical earth pressure is necessary for the safety and operation of these existing CCTs. Additionally, the vertical pressure should be properly considered in the design when new highfilled cutandcover tunnels (HFCCT) are being built in this area.
For a newly designed HFCCT, undoubtedly, the CCT can avoid the curved design and ensure the safety of highspeed railway in operation. Moreover, compared to the conventional CCT, the required thickness of the lining structure and bearing capacity of foundation is greater. In addition, when the lining is made thicker, the concrete cracks due to shrinkage after hardening will also significantly affect the durability of the tunnel structure.
Currently, the available theories for estimating vertical earth pressure cannot be properly applied to these necessary highfilled backfills. The underestimation of the vertical earth pressure for this HFCCT may result in nonconservative designs in the bearing capacity of lining structure and foundation beneath the HFCCT. This would result in the failure of the lining structure and the excessive differential settlements of the foundation. These factors will significantly affect the safety and durability of the HFCCT. On the contrary, if the vertical earth pressure is overestimated, a relatively conservative design will be made and significantly increase the cost of the project. Accordingly, a better understanding of earth pressure distribution mechanisms enables the engineers to precisely estimate the vertical earth pressure around these underground structures and make the proper designs for HFCCT.
2. Methods of Vertical Earth Pressure Estimation
The analytical and experimental solutions for earth pressure on culvert were first proposed by Marston [1, 2]. Later, it was determined that the loads on the underground structure are highly dependent on the relative movement between the central soil prism over the structure and the adjacent soil columns. Total vertical earth loads on the bottom slab of box culverts are sum of earth load on the top slab, frictional force on the sidewalls, and dead load of the structure [3, 4]. Therefore, the vertical earth pressure, based on embankment installation (EI) culvert, is usually larger than γh (γ: unit weight of backfill (kN/m^{3}), h: the height of backfill above the top of culvert (m)) [5–7]. The solution also shows that a greater difference between the actual vertical earth pressure and γh resulted when increases [8–14]. Conversely, another type of installation for culverts called trench installation (TI) results in a lower vertical earth pressure when compared to γh[7, 14–20]. Marston [2] found imperfect trench installation (ITI) based on this principles of mechanics, and engineers have developed this ITI method by means of adapting several types of compressible material over a buried pipe, such as leaves [21], baled straw [22], sawdust, woodchips [23, 24], or EPS ([25], Vaslestad and Johansen, 1993, [16, 17], Gu et al., 2005).
Currently, these analytical and experimental methods are widely used in practical engineering settings for estimating vertical earth pressure on the culverts with backfills and shallow pipe installations. Several specifications regarding vertical earth pressure estimation such as AASHTO [26], General Specifications for Design of Highway Bridges and Culverts [27], and Code for Design on Railway Bridges and Culvert [28] are based on MarstonSpangler (MS) theory and empirical results. However, these methods mainly consider the vertical earth pressure on culverts or pipes having crosssectional areas relatively small compared with HFCCTs. Furthermore, these methods lack proper considerations for the influences of slope angle of trench and the modulus of elasticity for the backfill materials. The estimation of the vertical earth pressure on the top of the CCT using the methods in the specifications of Code for Design of Road Tunnel [29] and Code for Design of Railway Tunnel [30] is based on the soil column pressure. These methods are applicable to these cases with the backfill heights which are smaller (usually between 24 meters). However, in the loading process of the HFCCTs, the computed magnitude of vertical earth pressure tends to be overestimated and underestimated with increasing backfill height. The loading mechanism and backfillCCT interactions were not properly considered in these specifications.
In a recent HFCCT project on Yuli railway, which is between Chongqing and Lichuan in China, the vertical earth pressure was computed using a method which considers the effect of slopes adjacent to the backfill. The vertical earth pressure was found to be reduced and agreed with the soilheld function proposed by Luo and Zhang [31]. However, the width of bottom in the valley between two slopes in which the CCT is sitting and the characteristics of HF soils are not considered in the equation.
Practically, the slope angle of the valley, the width of the bottom of valleys where the HFCCT is sitting, the size of CCT, and the modulus of elasticity for backfill are found to have significant influence on the vertical earth pressure distribution on the top of HFCCT [32]. Accordingly, these factors need to be well taken account of when estimating the vertical earth pressure formula for HFCCTs. The objective of this paper is to further investigate these influential factors such as shape of CCT’s crosssection, modulus of elasticity for backfill soil, the size of CCT’s crosssection and the coupling of slope angle and the width of bottom of valleys where the HFCCT is sitting. Through the parametric study using finite element analysis, it has been determined that the soil column pressure, γh, needs to be modified. Several coefficients corresponding to influential factors are proposed. A general formula for vertical earth pressure distribution on the top of HFCCT is recommended to satisfy the need for the HFCCT design.
3. HighFilled CutandCover Tunnel
The HFCCT includes the CCTs and massive backfills lying on top of those structures. Figure 1 shows a site of CCT called Changshoushan No. 2 in Gansu Province in Northwest China. This CCT was built in valleys, and new land will be reclaimed for the use of Auto Expo garden. The crosssection of the CCT and the needed backfills are shown in Figure 2. A total of 49 meters of backfill is needed from the top of the CCT to the level of the new land (see Figure 2). The backfills are mainly composed of the borrowed loess from surrounding areas in Loess Plateau in western China. The soil is compacted mechanically for each 0.5m. In the process of soil reclamation, the presupporting arch is used to lower the overburden pressure on the lining structure of tunnel to ensure the safety. The physical properties of the backfill soil are summarized in Table 1.

4. Analytical Model
According to previous studies, the vertical pressure was known to be either overestimated or underestimated. Both overestimation and underestimation have adverse effects on the design of CCT. It was also found that the holding effect due to the adjacent soil columns with slope can reduce the vertical earth pressure. The backfillstructuresurrounding soil interactions mechanism can be illustrated using two analytical models, shown in Figure 3. Figure 3(a) shows that the CCT is sitting on a valley bottom. This case shows steeper slopes adjacent to the CCT and the distance of valley bottom, B (B=D+2b, D: width of CCT (m), b: bottom distance between wall of the CCT to slope (m)), between two slopes is relatively narrow. Figure 3(b) shows another example which has less steep slopes adjacent to the CCT and a relatively wider bottom between two slopes (B is greater).
(a)
(b) : unit weight of backfill; : the average of earth pressure at a certain depth
For the former case, the backfill soil tends to slide along the slopes due to the gravitational force. As a result, the interaction force N and τ acting on the soil provided by the slopes contribute to the uplift forces resisting the relative movement between slopes and backfills subjected to selfweight. Therefore, the soil arching is formed to reduce the vertical earth pressure on the top of CCT (see Figure 3(a)). With the decrease of the slope angle, θ (°), the rotation of the normal stress, N, will reduce the arching effect formed in the backfill. If the slope angle, θ, further decreases and becomes less than the slip surface developed in the backfill due to Rankine active earth pressure (45°+φ′/2), the contribution of normal force, N, from the adjacent slopes will be less significant. The normal forces, , and the gravitational force will result in the relative settlements for the soils between the slip surfaces of Rankine active earth pressure and the slopes. Meanwhile, due to the settlement, the shear stress, , developed along the slip surfaces of Rankine active earth pressure will drag the backfill on the top of CCT downwards. Thus, with the increase of distance of b, the soil mass between the slip surface and the adjacent slope also increases and results in a greater drag down forces on the backfill on the top of CCT (see Figure 3(b)). However, with an increase of backfill height, the differential settlement between different soil columns (backfill column above the CCT and two backfill columns between slopes and slip surfaces) will be reduced. Although the vertical earth pressure increases, the rate of the increase is reduced. Therefore, for a HFCCT, the phenomenon showing the vertical earth pressure influenced by slope angle and B/D is evident.
Based on the analyses, the vertical earth pressure is relevant to several parameters such as D, B, θ, and γh. Meanwhile the relative settlement between different soil columns which are analyzed is related to stiffness of backfill soils, called the modulus of elasticity, E (MPa). In addition, the patterns of soil stress distribution on the top of soil are different between varying shapes of crosssection [33].
Therefore, the vertical earth pressure, q (kPa), on top of the HFCCT can be expressed as a function of the parameters mentioned above, using the following form: where represents the influence of crosssectional of the CCT.
5. Numerical Analysis
In order to verify the modified formulas for the vertical earth pressure on the top of HFCCT, finite element analysis was utilized. The parametric studies were conducted to correlate the results to each influential factor for earth pressure.
5.1. Finite Element Analysis
The ANSYS finite element code was used to investigate the earth pressure distribution on the HFCCT subjected to different types of reclamations. The twodimensional (2D) finite element analysis was used to conduct the parametric studies. Figure 4(a) shows the meshed model for EI and Figure 4(b) presents the meshed model for TI and the slope angle of the trench is 70°. The element used to model the HFCCT is planestrain. For the boundary conditions, the two vertical boundaries were restrained horizontally using rollers which only allow vertical deformation; the bottom boundary is fixed. The CCT and the slope were modeled using linearly elastic materials while the backfill for HF and the foundation soil beneath the CCT were modeled using MohrCoulomb elastoplastic materials. All the required parameters in the numerical analysis including the material properties for CCT, the engineering properties for backfills, and the foundation soils and the slope characteristics are summarized in Table 2. The soil properties were obtained from laboratory tests, including direct shear tests and the unconfined compression tests (see Table 2). In the finite element analysis, the interface elements between the soil and structure were simulated using contact element CONTA172 (backfill) and the target element TARGE169 (structure) which are defaulted in ANSYS.

(a)
(b) TI
5.2. Parametric Studies
In order to investigate the influence of each factor which was mentioned previously, the four corresponding coefficients, , , and , are proposed to modify the soil column pressure, γh, for backfills. The parametric studies are based on different shapes of CCT sections, the modulus of elasticity for backfill, width of CCT, slope angle, and the ratio of B/D. The present study considers two different shapes of CCT crosssections, rectangular and arch, with the following parameters: width of tunnel, D = 6.85m, 9 m and 12 m; modulus of elasticity for backfill, E = 10 MPa, 20 MPa, 40 MPa, 80 MPa and 120 MPa; slope angle, θ = 40°, 50°, 60° and 70°; the ratio of B/D = 1, 1.5 and 2. A maximum backfill height = 100 m was assumed for finite element analysis.
5.3. Analysis Results
Figure 5 shows the results of finite element analysis for the cases of D = 12 m and the height of backfill = 50 m. Figure 5 shows the contours of vertical stress distribution around the HFCCT for arch (see Figure 5(a)) and rectangular (see Figure 5(b)) crosssections, respectively. It can be seen that the stress distribution patterns are different around the CCT. The comparisons of the vertical earth pressure distribution on the top of HFCCT with different crosssectional shape are shown in Figure 6(a). For the rectangular crosssection, the peak value of earth pressure appears on both edges of the CCT, and the lowest value is shown in the middle of the CCT (U distribution). The maximum value appears symmetrically at nearly 0.5D away from the symmetric axis of CCT. For the arch crosssectional CCT, the maximum earth pressure appears in the middle (on the symmetric axis of CCT) and the lowest value at the edges of the top of CCT (inversed U distribution). Using the load equivalent method, the earth pressure distribution on the top of HFCCT can be presented using an equivalently uniformly distributed pressure shown in Figures 6(b) and 6(c). The comparisons between the maximum pressure and averaged equivalent pressure are shown in Figure 6(b) for arch crosssection and in Figure 6(c) for rectangular crosssection.
(a)
(b)
(a)
(b)
(c)
5.3.1. Influence of CrossSectional Shape
is the coefficient used to modify the overburden pressure based on different crosssectional shape of the HFCCT. The dimensionless ratio, , is defined as the equivalent pressure distribution for rectangular crosssection CCT (see Figure 6(b)) to the equivalent pressure for arch crosssection CCT (see Figure 6(c)) on the HFCCT with similar dimension, D. The coefficient can be correlated to another dimensionless ratio, h/D, which is the backfill height, h to the diameter of CCT, D. Three different dimensions of CCT for two different crosssectional shapes were analyzed: D=6.85m, 9 m and 12 m, respectively. The results showing the relationship between and h/D are presented in Figure 7.
As it can be seen in Figure 7, the effect of is significant when the h/D ratio is less than 2. If the h/D ratio is greater than 2, the coefficient of the three widths studied is close to 1.0 regardless of the h/D ratio. These results mean that the coefficient, , is independent on shape effect for a higher h/D. Thus, for the HFCCT, the = 1.0, regardless of the shape of CCT’s crosssections (rectangle or arch). That would also mean the effect of shape on the vertical earth pressure can be ignored for the HFCCT, but should be considered if h/D less than 2.
5.3.2. Influence of Stiffness of Backfill
According to previous studies, the influences of internal friction angle, ϕ, and the cohesion, c, of the backfill on vertical earth pressure at the top of HFCCT are less significant [32]. However, the modulus of backfill soil, E, was found to have significant influence on vertical earth pressure. In the numerical analysis, the strength parameters, c and ϕ of the backfill soil, were assumed to be invariable and homogeneous, and different values (10MPa, 20 MPa, 40 MPa, 80 Mpa, and 120 MPa) for backfill modulus of elasticity, E, were studied, respectively. Figure 8 shows the contours of vertical earth pressure distribution for the cases using E=10 MPa and 40 MPa. The equivalent earth pressures on the top of HFCCT for these two examples are 1.13MPa and 1.07MPa, respectively. Similarly, the analysis for other modulus values was also conducted. For E=20 MPa, 80 MPa, and 120 MPa, the equivalent pressures on the top of HFCCT are 1.10MPa, 1.04MPa, and 1.02MPa. is a coefficient which is corresponding to the influence of the stiffness of the backfill. The dimensionless ratio, , is defined as “the equivalent pressure distribution for different modulus of elasticity, E, to the equivalent pressure for modulus of elasticity, E=5.4 MPa on the HFCCT.” The relationship between and h/D is shown in Figure 9.
(a)
(b)
(a)
(b)
In Figure 9, the coefficient decreases when the ratio, h/D, increases. Meanwhile, the reduction rate increases if the modulus of elasticity, E, is greater (see Figure 8(a)). Based on these findings, each curve for the relationships between and h/D, considering the findings of different modulus of elasticity studies, can be mathematically expressed using the functions of h/D. The individual equation for subjected to different modulus of elasticity, E, is shown in Figure 9(a). To further demonstrate a possible correlation between the influence of the modulus, E, to , assuming = ·ln(h/D)+1, and =f(E), the relationship between and is plotted in Figure 9(b). The equation provided by simple logarithmic regression with R^{2}= 0.9721 is expressed in
Therefore, the general form of the can be written in a mathematical form as
Thus, the modulus of elasticity, E, has a positive influence on earth pressure reduction for a higher backfill. With the increase of fill height, the reduction is more notable. Thus, the modification is necessary for the HFCCT.
5.3.3. Influence of Dimension
The dimension of CCT specifically refers to the width of CCT. Figure 10 shows the results of finite element analysis for the vertical earth pressure distribution for the two different widths of CCT: 6.85m (Figure 10(a)) and 12 m (Figure 10(b)). It can be seen from the stress contours that the equivalent earth pressures are 1.30MPa and 1.13MPa, respectively. The results show that the tunnel width is relevant to the induced equivalent earth pressure, and a CCT with a smaller section will result in a higher equivalent earth pressure. A coefficient, , is introduced to modify the soil column pressure, γh, taking into account the effect of dimension.
(a)
(b)
Three different widths of CCT were studied: D=6.85 m, 9 m and 12 m. The coefficient for modifying the size effect of CCT was defined to be , having the equivalent pressure distribution for differing widths of tunnel to weight of soil column, γh, on the HFCCT. The results of versus h/D are plotted and shown in Figure 11. As it can be seen from Figure 10, the influence of the size of HFCCT nonlinearly increases with h/D. Meanwhile, the smaller size of CCT will result in a higher value of . The three curves of regression are found to be highly relevant to h/D and the corresponding equations for subjected to different widths of CCT are also shown in Figure 11. These corresponding to different width, D, are expressed in mathematical forms based on the function of h/D. The can be written in the form of = ·(h/D)^{0.1}, such that values of (1.2787, 1.2131 and 1.1604) are corresponding to D = 6.85m, 9 m, and 12 m, respectively. To further correlate the equation to the width of CCT, D, the relationship between and is plotted in Figure 11(b). Assuming = f(D), the equation for this linear regression curve with R^{2} = 0.9753 on Figure 11(b) can be expressed in
(a)
(b)
Accordingly, the general form for can be mathematically expressed using (5). The is a value which is greater than 1.0, and the value increases when the width of CCT decreases. As such, the effect of will amplify the overburden pressure, γh, and the amplification increases when the size of CCT decreases.
5.3.4. Influence of Slope Angle and B/D
Figure 12 shows that the results of finite element analysis for the vertical earth pressure distribution considering two different slop angles, θ, and the ratio, B/D. Figure 12(a) shows the example for θ = 40° and B/D = 2, while Figure 12(b) shows the case for θ = 70 ° and B/D = 1. It can be seen from the stress contours that the equivalent earth pressures are 1.02 MPa and 0.52 MPa, respectively. The results show that the coupling of θ and B/D are highly relevant to the induced equivalent vertical earth pressure. A CCT with a smaller θ and a greater B/D was found to result in a higher equivalent vertical earth pressure.
(a)
(b)
A coefficient, , is defined to modify the soil column pressure, γh, based on the slope angle and the ratio of B/D. It has been found that the width of bottom of valleys and the slope angle have some constraints on the high backfill. Therefore, the slope angle and the ratio of B/D should not be considered separately. The coupling influence of θ and B/D are presented in the mathematical form of .
In the finite element analysis, three different B/D (1.0, 1.5, and 2.0) and four slope angles, θ, (40°, 50°, 60° and 70°) were considered, respectively. For different B/D values, the results are as follows.
For B/D = 1.0 with θ = 40°, 50°, 60° and 70°, the analysis results are shown in Figure 13. In Figure 13(a), based on the relationships of versus h/D, four equations for corresponding to different θ are presented. To further include the effects of θ in the these forms, it is assumed = ·ln(h/D)+, while = f(θ) and = f(tan θ). Figure 13(b) shows the relationship between and θ. Through the regression using a straight line with R^{2} = 0.9985, the equation for is shown in Figure 13(b). Figure 13(c) shows the relationship between and tan θ. The equation of is given. Therefore, the can be mathematically expressed using
(a)
(b)
(c)
Similarly, for B/D = 1.5 and 2.0 and θ = 40°, 50°, 60° and 70°, the methods to obtain the general forms for and are shown in Figures 14 and 15, respectively. The mathematical forms of and are shown as
(a)
(b)
(c)
(a)
(b)
(c)
In order to obtain the general forms for , it is assumed =(k_{31}θ+) ln(h/D)k_{33} tan θ+1 (, and are the functions of B/D). The relationships for versus B/D, versus B/D and versus B/D are plotted in Figure 16, respectively. Based on the results shown in Figures 16(a), 16(b), and 16(c), the general form for is established. The mathematical for of can be expressed using
(a)
(b)
(c)
is a coefficient which is usually less than 1.0 and can be significantly reduced up to 70%, depending on the coupling effect of slope angle, θ, and the ratio, B/D. This reduction comes from the contribution of the normal holding stress of the valley slopes adjacent to the backfills. With the increase of the fill height, the reduction is more obvious. In practical engineering, if the influence for coefficient of can be considered carefully, the significant reduction on vertical earth pressure will enhance the safety of existing CCT and avoid an overconservative design by considering the vertical earth pressure for the HFCCT.
With these four coefficients determined and demonstrated as influential factors, the vertical earth pressure on the top of the HFCCT should be modified as follows:where , , and are the influence coefficients of crosssection shape, elasticity modulus, and the width of CCT, respectively. is the coupling effect of slope angle, θ, and the ratio, B/D.
6. Numerical Verification of Proposed Method
In order to verify the feasibility of this proposed equation which includes four coefficients of modification: , , , and , an example shown below is a HFCCT that has a dimension in width, D = 13.76. The modulus of elasticity of soil, E, is 20 MPa, and the other properties are consistent with those parameters shown in Table 1. Two cases (Case I: B/D= 1.5, θ= 60°; and Case II: B/D = 2.0, θ= 40°) were considered, and the computed results using this proposed method are compared to results using soil column theory (γh) and finite element analysis using ANSYS. The comparisons are shown in Figure 17. The average relative error is 1.14% for the results using the proposed method and the numerical analysis. The proposed method is effective in avoiding the complicated numerical analysis considering all influential factors.
Selecting h=50m, the detailed steps to obtain each coefficient to modify the soil column pressure, γh, and use the proposed method are shown as follows.
(1) Calculation of . Since the sectional shape of a CCT was found to have little influence on the vertical earth pressure for the HFCCT, thus, assume = 1 is reasonable.
(2) Calculation of . Calculating using equation (3), =0.959 is obtained.
(3) Calculation of . According to equation (5), the =1.273 is computed as follows:
(4) Calculation of . Calculating according to equation (9), the detailed calculations and = 0.669 are shown as follows:After determining , , , and , the final vertical earth pressure can be modified using equation (10) for the HFCCT (see below). The vertical earth pressure is 722.80 kPa on the top of HFCCT for the backfill height, h=50m:The results for other backfill heights can be obtained using the same mathematical procedures and substituting the appropriate values. Close examination of Figure 17 leads to several findings that can be addressed herein.
7. Experimental Verification of Proposed Method
In order to verify the feasibility of this proposed equation which includes four coefficients of modification: , , , and , we designed the experimental test. A rigid steel tank with dimensions of 180 cm ×120 cm × 150 cm (length×width×height) for the tests and the dimensions of the model CCT with either an arched or rectangular crosssectional shape is 120 cm in length, 28 cm in height, and 56 cm in width. Figure 18(a) presents a crosssectional representation of the HFCCT apparatus that includes both the rectangular and arched CCT’s design. In the model, the two slopes were formed using concrete, which is assumed to be rigid. The slope angle of the two slopes in the tank was defaulted to be 70°. In order to record the distribution of the earth pressure on top of the CCT, nine earth pressure cells (No. 1 to No. 9) were installed on top of the firstlayer soil to record the earth pressure distribution on top of the CCT, as shown in Figure 18(b).
(a)
(b)
Table 3 summarizes the conditions for the two tests. S1 and S2 had rectangle and arch crosssectional shape, respectively.

The physical and engineering properties of the backfill soil from the laboratory tests are summarized in Table 4.

In order to model the actual backfill construction over an HFCCT, the CCT was installed between the slope firstly; secondly, fill and compact the first layer of soil above the CCT, and install 9 earth pressure cells on top of the first layer; thirdly, repeat the fill and compaction process for the subsequent layer(s) until the maximum backfill height is reached. Meanwhile, measure the earth pressure while filling the backfill at each stage.
Figure 19 compares the computational results for experiments, numerical analyses, and the modified equation for the backfill height of 90 cm. The results summarized in Figure 19 are based on a vertical earth pressure coefficient, λ = vertical earth pressure on top of CCT/γh. This results show that the estimations of vertical earth pressure on the top of HFCCT using the proposed method are nearly consistent with the results using finite element analysis; the experimental results are slightly smaller than the results using proposed method. The average error between experimental results and proposed method results is about 4.2%. This error occurs due to the soil compaction process for each layer of backfill.
8. Discussions
The results show that the modification for γh is necessary for the HFCCT; the deviation increases as the backfill height increases. On the contrary, for the shallow backfill, the results have little difference regardless of estimation methods used.
The vertical earth pressure may be greater or smaller than γh, which is as shown as Figure 17. The greater deviation between the computed value and γh is highly relevant to the coupling effect of slope angel and D/B ratio. To avoid overestimating and underestimating the vertical earth pressure for the HFCCT in design, the is the most significant parameter to modify and reduce the vertical earth pressure.
9. Conclusion
The conventional methods for estimating vertical earth pressure are not applicable to the HFCCT. Several influential factors such as shape of CCT’s section, stiffness of backfill soil, dimension of CCT, slope angle, and the D/B ratio should be carefully considered in order to more precisely estimate the vertical earth pressure on the top of HFCCT. Based on these parametric studies using finite element analysis, the coefficients for modifying the soil column pressure, γh, are proposed. The results were verified using both numerical analyses and experiments. Several conclusions can be drawn.(1)The theory of soil column method is not applicable to the HFCCT. The deviation increases with the height of backfill increasing. Therefore, the estimation of vertical earth pressure must be modified.(2)Four coefficients are proposed to modify the soil column pressure, γh, using finite element analysis. Four general forms corresponding to these coefficients are , crosssectional shape of CCT effect; , stiffness of backfill effect; , width of CCT effect; and , coupling effect of slope angle, θ, and the ratio of B/D effect.(3)For existing or newly designed HFCCT, the coefficient of is 1.0, regardless of crosssectional shapes; is less than 1.0 and can be as low as about 0.8, depending on the height of backfills; is a coefficient which is always greater than 1.0. The increases with the height of backfills and decreases with the dimension of CCT; is equal to or less than 1.0 depending on the coupling effect of slope angle of the valley adjacent to the backfills and the ratio, B/D.(4)For the overall impacts of these coefficients on earth pressure for the HFCCT, the coefficients of and reduce the vertical earth pressure, the will amplify the vertical earth pressure, and has little influence on earth pressure. Therefore, to ensure the safety and reduce the vertical earth pressure on the top of HFCCT, proper consideration of the coefficients of and is necessary before and during the construction process.
Data Availability
The summarized resulted data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was supported by the National Science Foundation of China (51668036, 51868041), General Projects of Scientific Research of High School in Gansu (2017A111), the Program of Changjiang Scholars and Innovative Research Team in University (IRT_15R29), Support Program Funding of Collaborative Innovation Technology Team in Gansu Provincial (2017C08) and Projects of Excellent Research Team in Lanzhou Jiaotong University (201606), China, Energy Geomechanics Laboratory at the University of North Dakota, and Mr. Daniel Brunson at the University of North Dakota.
References
 A. Marston and A. O. Anderson, The Theory of Loads on Pipes in Ditch and Tests of Cement and Clay Drain Tile and Sever Pipe, Iowa Engineering Experiment Station, Ames, Iowa, USA, 1913.
 A. Marston, The Theory of External Loads on Closed Conduits in the Light of the Latest Experiments, Iowa Engineering Experiment Station, Ames, IA, USA, 1930.
 J. Kang, F. Parker, and C. H. Yoo, “Soilstructure interaction and imperfect trench installations for deeply buried concrete pipes,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 133, no. 3, pp. 277–285, 2007. View at: Publisher Site  Google Scholar
 J. Kang, F. Parker, Y. J. Kang, and C. H. Yoo, “Effects of frictional forces acting on sidewalls of buried box culverts,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 32, no. 3, pp. 289–306, 2008. View at: Publisher Site  Google Scholar
 W. V. Binger, “Discussion to underground conduitsan appraisal of modern research,” Proceedings of the American Society of Civil Engineers, vol. 73, pp. 1543–1545, 1947. View at: Google Scholar
 M. Z. Yang, Evaluation of factors affecting earth pressures on buried box culverts [Ph.D. thesis], University of Tennessee, Knoxville, Tenn, USA, 2000.
 B. G. Chen and L. Sun, “Performance of a reinforced concrete box culvert installed in trapezoidal trenches,” Journal of Bridge Engineering, vol. 19, no. 1, pp. 120–130, 2014. View at: Publisher Site  Google Scholar
 M. G. Spangler, “Underground conduits—an appraisal of modern research,” Transactions of the American Society of Civil Engineers, vol. 113, no. 1, pp. 316–345, 1948. View at: Google Scholar
 D. H. Trollope, M. G. Speedie, and I. K. Lee, “Pressure measurements on tullaroop dam culvert,” in Proceedings of the 4th Australia–New Zealand Conf. on Soil Mechanics and Foundation Engineering, In stitution of Engineers, pp. 81–92, Barton, ACT, Australia, 1963. View at: Google Scholar
 A. D. M. Penman, J. A. Charles, J. K. T. L. Nash, and J. D. Humphreys, “Performance of culvert under winscardam,” Geotechnique, vol. 25, no. 4, pp. 713–730, 1975. View at: Google Scholar
 A. Q. Gu, “Investigation of the vertical earth pressure on projecting conduit and underground chamber under a high embankment,” Chinese Journal of Geotechnical Engineering, vol. 3, no. 1, pp. 3–15, 1981 (Chinese). View at: Google Scholar
 A. Dasgupta and B. Sengupta, “Largescale model test on square box culvert backfilled with sand,” Journal of Geotechnical Engineering, vol. 117, no. 1, pp. 156–161, 1991. View at: Publisher Site  Google Scholar
 R. M. Bennett, S. M. Wood, E. C. Drumm, and N. R. Rainwater, “Vertical loads on concrete box culverts under high embankments,” Journal of Bridge Engineering, vol. 10, no. 6, pp. 643–649, 2005. View at: Publisher Site  Google Scholar
 K. Kim and C. H. Yoo, “Design loading on deeply buried box culverts,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 131, no. 1, pp. 20–27, 2005. View at: Publisher Site  Google Scholar
 J. Vaslestad, T. H. Johansen, and W. Holm, “Load reduction on rigid culverts beneath high fills: longterm behavior,” Transportation Research Record, Transportation Research Board, no. 1415, pp. 58–68, 1993. View at: Google Scholar
 R. P. McAffee and A. J. Valsangkar, “Field performance, centrifuge testing, and numerical modelling of an induced trench installation,” Canadian Geotechnical Journal, vol. 45, no. 1, pp. 85–101, 2008. View at: Publisher Site  Google Scholar
 B. L. McGuigan and A. J. Valsangkar, “Centrifuge testing and numerical analysis of box culverts installed in induced trenches,” Canadian Geotechnical Journal, vol. 47, no. 2, pp. 147–163, 2010. View at: Publisher Site  Google Scholar
 X.W. Yang and Y.X. Zhang, “Study on arch action and earth pressure theory for culverts under high embankment,” Journal of Rock Mechanics and Engineering, vol. 24, no. 21, pp. 3887–3893, 2005 (Chinese). View at: Google Scholar
 J.J. Zheng, B.G. Chen, and S.B. Zhang, “Experimental investigation and numerical simulation of nonlinear earth pressure for trenchburied culverts,” Chinese Journal of Geotechnical Engineering, vol. 30, no. 12, pp. 1771–1777, 2008 (Chinese). View at: Google Scholar
 Y.G. Li and S.Y. Zhang, “Experimental and theoretical study on earth pressure on top of rectangular trenchburied culverts,” Yantu Lixue/Rock and Soil Mechanics, vol. 29, no. 4, pp. 1081–1086, 2008 (Chinese). View at: Google Scholar
 M. G. Spangler, “A Practical application of the imperfect ditch method of construction,” in Proceedings of the ThirtySeventh Annual Meeting of the Highway Research Board, Highway Research Board, vol. 37, pp. 271–277, Washington, DC, USA, 1958. View at: Google Scholar
 N. G. Larsen and J. G. Hendrickson, “A practical method for constructing rigid conduits under high fills,” in Proceedings of the 41st Annual Meeting of the Highway Research Board, Highway Research Board, vol. 41, pp. 273–280, Washington, DC, USA, 1962. View at: Google Scholar
 R. P. McAffee and A. J. Valsangkar, “Geotechnical properties of compressible materials used for induced trench construction,” Journal of Testing and Evaluation, vol. 32, no. 2, pp. 143–152, 2004. View at: Publisher Site  Google Scholar
 R. P. McAffee and A. J. Valsangkar, “Performance of an induced trench installation,” Transportation Research Record, Transportation Research Board, vol. 1936, no. 1, pp. 230–237, 2005. View at: Google Scholar
 J. A. Sladen and J. M. Oswell, “The induced trench method—a critical review and case history,” Canadian Geotechnical Journal, vol. 25, no. 3, pp. 541–549, 1988. View at: Publisher Site  Google Scholar
 AASHTO, AASHTO Standard Specifications for Highway Bridges, Washington, DC, USA, 2010.
 China Communications Construction Company Highway Consultants CO., General Specifications for Design of Highway Bridges and Culverts, Beijing, China, 2015.
 China Railway Design Corporation, Code for Design on Railway Bridges and Culvert, Beijing, China, 2017.
 China Merchants Chongqing Communications Technology Research & Design Institute CO., Code for Design of Road Tunnel, Beijing, China, 2004.
 China Railway Eryuan Engineering Group CO., Code for Design of Railway Tunnel, Beijing, China, 2017.
 L. S. Luo and P. Zhang, “Preliminary study on load of high fill open cut tunnel,” High Speed Railway Technology, vol. 5, no. 3, pp. 90–93, 2014 (Chinese). View at: Google Scholar
 S. Li, L. Ma, Q. C. Wang, J. X. Li, W. W. Li, and Y. J. Zhang, “Calculation method and influencing factors for earth pressure of high fill open cut tunnel,” China Railway Science, vol. 37, no. 5, pp. 41–49, 2016 (Chinese). View at: Google Scholar
 A. N. Dancygier, Y. S. Karinski, and A. Chacha, “A model to assess the response of an arched roof of a lined tunnel,” Tunnelling and Underground Space Technology, vol. 56, pp. 211–225, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Sheng Li 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.