Analytical Solutions for the Mechanical Responses of Shallow Double-Arched Tunnel Subjected to Symmetric Loads

Double-arched tunnel is a special kind of tunnel structure which is widely adopted for shallow rock tunnels. +e internal forces within a double-arched tunnel are difficult to be determined due to its special geometry and complex interactions with surrounding rock masses. +is paper presents a set of analytical solutions for determining the internal forces within the lining structure of shallow double-arched tunnel subjected to symmetric external loads. +e double-arched tunnel is decomposed into three main parts, i.e., the arch rings, the side walls, and the middle wall. +e force method and the elastic foundation beam model are employed for the arch rings and the side walls, respectively, while the middle wall is treated as a cantilever beam fixed at the bottom. Analytical solutions for the internal forces of the three parts are derived separately, which are continuous at the conjunctions of different parts. +e derived analytical solutions are verified by comparing with the FEM simulation results. Finally, parametric studies are performed to investigate the influences of burial depth, elastic resistance coefficient, opening angle, and tunnel span on the internal forces based on which some recommendations are provided for the construction and design of the double-arched tunnels. +e derived analytical solutions provide fast estimations for the internal forces and deformation of the double-arched lining structures, which will be a useful tool for design optimization.


Introduction
e double-arched tunnel is a special kind of tunnel whose lining structures between two adjacent tunnels are connected and supported by the middle wall [1].
Unlike single tunnels with a circular or straight wall arch cross section which have been profoundly investigated [2][3][4][5][6][7][8][9], due to its special geometry, the distribution of internal forces of double-arched tunnels has not yet been well understood. Yan et al. set up a comprehensive monitoring program at critical sections of a double-arched tunnel to measure the internal forces in both steel arches and secondary linings [10]. e monitoring and analysis results show that the forces of lining structures are directly related to multiple factors, including the construction procedures, geological conditions, and locations in the double-arched tunnel. ey further showed that the arch vault is the most critical section, which was also observed by Yuan et al. [11]. rough in-site monitoring, Lai et al. found that the critical areas of double-arched tunnels are mainly located on the middle wall, side wall, and vault [12]. Zhang et al. carried out a model test for an unequally spanned double-arched tunnel to study the distribution of surrounding rock pressure and mechanical characteristics and found that the tunnel was in a small eccentric bending state, and thrust stress in the side walls increases with increasing span [13]. He et al. conducted model tests to investigate the distribution of internal forces of nonsymmetric double-arched tunnels in sand-cobble ground [14]. ey showed that the internal force of the structure increased linearly with the increase in external load until the first crack appeared, after which the internal force redistributed but the structure still had a certain bearing capacity. Li et al. derived the mechanical responses of supporting structure for a shallow-buried double-arched tunnel under the China Great Wall based on reduced-scale model tests [15].
Numerical simulations have also been conducted to analyse the stability and mechanical behaviours of doublearched lining structure. Yu and Yang analysed the internal force of the lining of the double-arched tunnel by using the loading structure method, considering the situations of both single-side and double-side bearings [16]. Zhang and He performed numerical simulations on the mechanical behaviour of middle wall of double-arched tunnel and found that the main factors affecting the loads of middle wall include the burial depth, the properties of surrounding rock, the construction method, and the eccentric loading condition [17]. Based on the coupled fluid-solid theorem, Li et al. investigated the stability of surrounding rock and internal forces of double-arched tunnel under seepage condition [18]. Mao et al. simulated the construction process of a double-arched tunnel in loess area considering the lateral recharge and fluid-solid coupling effect, which showed that the water disaster susceptible areas mainly are concentrated in the middle wall, arch springs of both tunnels, and tunnel face [19]. Liu et al. studied the acceleration response of a shallow-buried biased double-arched tunnel under the action of Wenchuan wave [20]. Based on both physical model testing and numerical simulation results, Min et al. reported that the presence of voids on the top of the middle wall may induce adverse effects on the lining, which are characterized by significant change of internal forces and increase in lining deformation [21]. Similarly, Zhang et al. showed that the internal forces change significantly at the area in the close vicinity of the void [22]. e aforementioned studies are either costly or timeconsuming. In contrast, analytical approach is a convenient, fast, and effective method for analysing the mechanical behaviour of both surrounding ground and tunnel lining [23][24][25]. Based on Airy's stress function, Nagger and Hinchberger presented a closed-form solution for bending moment and thrust force of tunnel linings that can be idealized as an inner jointed segmental lining and an outer thick-walled cylinder embedded in a homogeneous infinite elastic soil or rock, accounting for the interaction between the lining and the surrounding ground [26]. Lee et al. adopted the force method to predict internal forces and displacements of a jointed circular segmental tunnel lining by treating the longitudinal joints as a set of linear rotational springs [27]. Wang et al. presented an analytical solution for circular segmental tunnel lining based on the state space method [28]. e curved Euler beam theory and the principle of minimum potential energy were adopted during the derivation of governing equations. Wang et al. further derived a unified analytical solution for the circular segmental tunnel lining irrespective of the layout of segmental assemblies between adjacent rings by modelling the interaction between the segments as a set of springs with threedirectional resistances [29]. Talmon and Bezuijen presented an analytical model to analyse the longitudinal bending moment in the tunnel lining during the construction phase of bored tunnels, taking into account the influence of forces from the TBM, the loading of the tunnel lining by the grout in the liquid phase, and linear elastic properties of the tunnel lining and soil [30]. Zhao et al. employed the commonly accepted analytical solutions to calculate the lining forces and deformations for both shallow and deep tunnels [31]. Park et al. developed an analytical solution by considering the relationships between displacements and interaction forces at the soil-lining interface and including the slippage effect at the soil-lining interface using a spring-type flexibility coefficient [32]. Bobet presented closed-form solutions for the liner stresses of a deep circular tunnel excavated in transversely anisotropic rock [33]. Kargar et al. provided a semianalytical elastic plane-strain solution for stress field around a lined noncircular tunnel at great depth subjected to uniform ground load [34]. e analytical solution for an arbitrarily shaped tunnel excavated in anisotropic rock mass supported by an isotropic lining was presented by Wang et al. [35]. e equations and solutions of stresses and displacements in lining and surrounding rock mass were derived based on the complex variable method.
Most of the research studies about analytical approach mentioned above mainly focus on the internal forces or stresses of single tunnel lining, especially for the circular tunnel due to the simplification and convenience in the derivation. Hu et al. presented an analytical solution for the determination of internal forces of a jointed segmental precast DOT (double-O-tube) lining by combining the unbalanced force transition method with the force method [23].
is is the only analytical work available for the double tunnel lining. However, it should be noted that the shape of the DOT lining is still relatively regular when comparing with the double-arched tunnel which is normally composed of arch rings, straight side walls, and middle wall. So far, no analytical solutions for the internal forces of double-arched tunnel have been reported in the literature.
In this paper, analytical solutions for determining the internal forces of lining structure in shallow double-arched tunnels subjected to symmetric loading are presented. e double-arched tunnel is decomposed into three parts, including the arch rings, the side walls, and the middle walls. e force method is applied to determine the internal forces of arch rings.
e elastic foundation beam model is employed for representing the side walls. e middle wall is treated as a cantilever beam fixed at the bottom. e internal forces are continuous at the conjunctions of the three different parts. e internal forces of a background project calculated using the proposed analytical solutions are compared with the values obtained by the numerical simulations. Finally, based on the obtained analytical solutions, parametric analyses are conducted to investigate the influences of some key influential factors on the internal forces of the double-arched tunnel lining.

Model Formulation
(1) e double-arched tunnel is shallowly buried (2) e middle wall of the double-arched tunnel is straight and can be taken as a cantilever beam while the side wall is assumed as an elastic foundation beam (3) e bottom of the side wall is fixed in both horizontal and vertical directions but can rotate under the elastic resistance, while the bottom of the middle wall is fixed (4) e inverted arch is not considered in the derivation of the analytical solution (5) e analytical solutions for the internal forces along the tunnel arch, the middle wall, and the side walls are derived separately (6) e force method and elastic foundation beam method are used for the derivation of internal forces along the tunnel arch and side wall, respectively (7) e thicknesses of the tunnel lining and the middle wall are not considered during the derivation and the external dimensions of the lining are employed (8) e tunnel lining is considered as linear elastic Based on the above assumptions, the original mechanical model presented in Figure 1(a) could be simplified into Figure 1(b). e definitions of symbols presented in Figure 1 are given as follows: H is the burial depth at the crown of the double-arched tunnel; H 1 is the height from the top of the middle wall to the ground surface; H 2 is the height of the double-arched tunnel; B is the total width of the double-arched tunnel; R is the outer radius of the arch ring of left tunnel; α 1 , α 2 , and α 3 can be considered as the growth angle of lateral earth pressure along depth of the corresponding three kind of pressures, which are given by tan α 1 � (e 2 − e 1 )/H 2 and tan α 2 � (e 2 ′ − e 1 ′ )/f 2 , α 3 � q ′ /l 2 ; a denotes the bottom point of the side wall in the left tunnel; h is the foot of left arch of left tunnel on the top of the side wall; d is the foot of right arch of left tunnel on the top of middle wall; l is the height of the side wall of left tunnel; l 1 and l 2 are the widths of left and right half arch rings of left tunnel, respectively; f 1 and f 2 are the height of left and right half arch rings of left tunnel; and h a denotes the bottom width of the left side wall. A prime is added to the same set of symbols to denote the corresponding quantities of the right tunnel. e external pressure from the surrounding rock mass is determined under the assumption that the doublearched tunnel is shallowly buried [36,37]. e cut-off altitude (H p ) separating deep tunnel from shallow tunnel can be considered as follows: where h p is the load equivalent height, which can be calculated according to the approach proposed by [36,37]. As shown in Figure 1(a), based on H p and h p , the pressure of shallow double-arched tunnel can be determined according to two different conditions as shown in Table 1 [36].
In Table 1, c is the unit weight of the surrounding rock mass, φ is the calculated friction angle of surrounding rock masses, λ is the lateral pressure coefficient, β is the angle between the failure surfaces and the horizontal for both sides of the tunnel, and θ is the friction angle of sliding surface, which can be determined according to Table 2.

Model Derivation
In this section, the analytical solutions for the internal forces of arch ring, side wall, and middle wall are derived separately based on the simplified model given in Figure 1(b). Since the structure configuration and pressure distribution are both symmetric, only the left part of the tunnel structure is employed in the following derivation. e internal forces of the middle wall can be obtained by the superposition of corresponding quantities from both arch rings.

Analytical Solutions for the Arch Rings.
e simplified models of the left arch ring for the force method are presented in Figure 2. All the pressures exerted on the arch rings shown in Figure 2 are considered as active pressures, which are all positive. e redundant internal forces of the crown section of the arch ring induced by the active pressures are X 1p , X 2p , and X 3p , respectively. e angular displacement of the arch foot turning outwards, the outward horizontal displacement, the upward vertical displacement of the right arch foot, and the downward vertical displacement of the left arch foot are positive.
β pl , u pl , and v pl represent the angular, horizontal, and vertical displacements of the left arch foot, respectively, while β pr , u pr , and v pr denote the same quantities of the right arch foot.
According to the displacement coordination condition of the force method at the crown of the arch ring, the displacement equations are as follows: where δ ij � (i, j � 1, 2, 3) are the flexibility coefficients of the arch ring and ∆ ip is the displacement induced by the active pressure in the direction of X ip at the cross section of the tunnel crown. By using the superposition principle, the displacements of the left and right arch feet are obtained as follows: where β il , u il , and v il are the angular, horizontal, and vertical displacements of the left arch foot under the unit loads Advances in Civil Engineering Substituting the displacements in equation (3) into equation (2), the following displacement equations containing the unknown X 1p , X 2p , and X 3p are obtained: a 11 X 1p + a 12 X 2p + a 13 X 3p + a 10 � 0 where the coefficients a ij are determined by the following: By solving equation (5), X 1p , X 2p , and X 3p can be determined as follows: where Δ � e unknown coefficients for the arch ring can be solved based on the knowledge of structural mechanics. e basic structure of the arch ring under the unit loads X i � 1(i � 1, 2, 3) is shown in Figure 3. e internal forces of the arch ring corresponding to the unit loads X 1 � 1, X 2 � 1, and X 3 � 1 are presented in the following equations, respectively: left half arch ring: right half arch ring: left half arch ring: right half arch ring: Since the effect of axial and shear forces on displacement is relatively small, usually only the effect of the bending moments is considered [27]. However, in order to fully consider the internal forces, the axial and shear forces are also involved. According to the structure mechanics, the flexibility coefficients δ ij for the arch ring can be calculated as follows: where M i (M j ), N i (N j ), and Q i (Q j ) are the bending moments, axial force, and shear force of the arch ring under the unit loads X i (X j ) � 1(i, j � 1, 2, 3) at the crown of the arch ring of the basic structure; E′ is the elastic modulus of the arch ring; G is the shear modulus of the arch ring; I is the polar moment of inertia of the arch ring's cross section with respect to the longitudinal direction; A is the area of the cross section; k is the nonuniform coefficient of the shear stress, which is taken as 1.2 for the rectangular cross section in this paper; and s is the arch length.
Substituting equations (8)-(10) into equation (11) for i, j � 1,2,3, respectively, the flexibility coefficients of the left tunnel can be obtained as follows: Advances in Civil Engineering

Advances in Civil Engineering
Note that φ i ranges from 0 to φ h for the left half arch ring and 0 to φ d for the right half arch ring. e flexibility coefficients of the right tunnel can be obtained by replacing φ h and φ d by φ d ′ and φ h ′ in the same formulae of the left tunnel. e internal forces caused by active pressures can be determined by accumulating the pressures as shown in Figure 4, namely, the vertical uniform pressure q, the vertical triangular pressure q′, and the lateral trapezoidal pressure e + Δe at both sides which can be considered as uniform pressures e � e 1 for the left half arch ring and e � e 1 ′ for the right half arch ring with the triangular pressure Δe � f 1 tan α 1 for the left half arch ring and Δe � f 2 tan α 2 for the right half ring.
where M o ip , N o ip , and Q o ip denote the bending moment, axial force, and shear force under the active pressures, respectively.
For the vertical uniform pressure q, left half arch ring : For the vertical triangular pressure q ′ , For the lateral uniform pressure e, left half arch ring: For the triangular pressure ∆e, left half arch ring : e displacement coefficients ∆ ip for the arch ring can be calculated as follows: Substituting equation (14) and equations (8)-(10) into equation (18) yields the following displacement coefficients of the arch ring under the vertical uniform pressure q: Using the same method, the corresponding displacement coefficients for other active pressures can be obtained as follows:

Advances in Civil Engineering
Advances in Civil Engineering e displacement coefficients under the active pressures are the superposition of the obtained coefficients, which can be expressed by the following: Accordingly, the displacement coefficients for the right tunnel can be obtained by utilizing the same method.
After the coefficients for the arch ring have been obtained, the remaining unknown coefficients at the top of the side wall and middle wall could be calculated by the elastic beam model and the force analysis on the top of the middle wall, which are shown in Appendix A and B, respectively.
Based on the obtained coefficients, the redundant forces X 1p , X 2p , and X 3p can be obtained by equation (7). Based on equation (7), the internal forces along the arch ring can be determined according to the static equilibrium condition, which are given by equations (22) and (23) for the left and right half arch rings, respectively. e bending moment M ia is positive when the intrados of the arch ring is in tension and the axial force N ia is positive when the arch ring is under compression, while the shear force Q ia leading to the clockwise rotation of the arch section is taken as positive.
where M o ip , N o ip , and Q o ip are the internal forces for arbitrary positions of the arch ring subject to the active pressures and φ i is the angle between the analysed cross section and crown of the arch ring, which is in the range of 0∼φ h for the left half arch ring and 0∼φ d for the right half arch ring.
Because of symmetry and by using the same methods, the analytical solutions of internal forces for the right tunnel arch ring can be obtained accordingly.

Analytical Solutions for the Side Wall.
e internal forces of the arch rings are transferred continuously to the side wall through the arch foot, as shown in Figure 5(a). For the side wall, the bending moment M is is positive when the intrados of the side wall is in tension and the axial force N is is positive when the side wall is under compression, while the shear force Q is leading to the clockwise rotation of the side wall is positive.
In general, the angle of the joint of the side wall and arch foot at the crown of the arch ring, namely, φ h for the left tunnel and φ h ′ for the right one, are usually right angle and it is unnecessary to decompose the internal forces. Nevertheless, in order to embrace more general situations, φ h and φ h ′ are considered to be arbitrary angles in the current derivation.
e corresponding decompositions of the internal forces at the top of the side wall are shown in Figure 5(b).
In the derivation of the analytical solutions, the side wall is assumed to be an elastic foundation beam. e internal forces induced by the unit bending moment, unit horizontal force, and vertical force exerted on the top of the side wall are first derived separately, which can refer to Appendix A.
e real internal forces are obtained by applying the principle of superposition. For the left tunnel, the final solutions for the internal forces of the side wall are obtained as follows: where M M�1 and Q M�1 , M H�1 and Q H�1 , and M V�1 and Q V�1 are the bending moments and shear forces of the side wall induced by the unit bending moment exerted on the top of the side wall, induced by the unit horizontal force exerted on the top of the sidewall, and induced by the unit vertical force exerted on the top of the sidewall, respectively. M e and Q e ; M ∆e and Q ∆e are the bending moments and shear forces induced by the uniform pressure exerted on the side wall along the vertical direction and induced by the triangular pressure exerted on the side wall along the vertical direction, respectively. ey are all related with the elastic resistance according to the elastic foundation beam method, with the derivation given in Appendix A. e analytical solutions of internal forces for the right tunnel side wall can be obtained in a similar way.

Analytical Solutions for the Middle Wall.
e middle wall is assumed as a cantilever beam. e internal forces of the arch rings are transferred continuously to the middle wall through the arch feet connected with the middle wall. e bending moment M im causing the counterclockwise rotation of the middle wall is positive, the shear force Q im leading the middle wall to rotate clockwise is positive, and the compressive axial force N im is positive. e forces at the top of the middle wall and their decomposition are shown in Figure 6(a).
Based on Figure 6, the analytical solutions of internal forces for the middle wall can be obtained as follows: Equations (22) to (25) constitute the complete analytical solutions of the internal forces of shallow double-arched tunnel for the left tunnel and middle wall. Detailed information on how to determine the unknown coefficients is presented in Appendix A and B. Figure 7 summarizes the stepwise procedure of deducting the internal forces of double-arched tunnel using the analytical solutions derived in the current study:

Main Steps to Determine the Internal Forces of Double-Arched Tunnel.
(1) determination of the outer dimension and lining characteristics of a double-arched tunnel. e outer dimension of the tunnel should be determined according to the geometry design, including the outer radius of arch ring, the height of side and middle walls, the total width of tunnel, as well as the opening angle. e lining characteristics including the lining thickness and properties should also be determined.
(2) determination of the pressures on a double-arched tunnel. Based on the burial depth and geological conditions, the pressures on the double-arched tunnel can be derived according to the information listed in Figure 1 and Table 1. (3) determination of the flexibility coefficients and displacement coefficients. Based on the aforementioned methods and equations, the flexibility coefficients (equation (12)) and displacement coefficients (equations (19) to (21)) of arch ring can be determined, respectively. e displacement coefficients of the side wall and middle wall can be determined according to the method in Appendix A and B. (4) determination of the internal forces of the arch ring from equations (22) and (23). (5) determination of the internal forces of the side wall and middle wall from the equations (24) and (25)     Advances in Civil Engineering e derived analytical solutions in this section mainly focus on the double-arched tunnel subjected to symmetric external loads, while the nonsymmetric scenarios are more commonly encountered in engineering practice. e main differences between the two scenarios are the pressures exerted on the tunnel. For the nonsymmetric ones, the vertical pressures q in Figure 1 will be trapezoidal instead of rectangle; thus, another growth angle of vertical earth pressure α 4 should be added [38], similar to the lateral earth pressure outside of double-arched tunnel. erefore, the derivation process of displacement coefficients of the arch ring (equations (14) and (19)) should consider the effect of α 4 , which is the main difference in the derivation process between the two scenarios. In addition, for deep tunnels, such as the right tunnel, a larger lateral pressure outside the tunnel will be present, and another growth angle α 5 of the lateral pressure may not equal to the present α1 in Figure 1. Nonetheless, the deduced process of the right-side tunnel is still similar to the left one. e remaining part of the deduction of the analytical solution for the nonsymmetric scenarios is the same with the aforementioned process. Overall, it is convenient to extend the present analytical solution to the more common nonsymmetric scenarios.

thickness and properties) Preparation
Detuction process Determine the surrounding pressures exerted on the tunnel according to the information listed in Figure 1 and Table 1 Calculate flexibility coefficients and displacement coefficients of arch ring (equationes (12), (19)- (21)) Calculate displacement coefficients of the side wall (Appendix A) and middle wall (Appendix B) Calculate the final internal forces along the arch ring (equationes (22) and (23) Figure 8(a), the adopted double-arched tunnel model referred to the project introduced by [10]. Since the dimension of the side wall was not clearly defined in [10], for simplicity, the height of the side wall is roughly set to be 3m according to the length scale as presented in [10]. e thickness of the sprayed concrete and the secondary lining is 25 cm and 50 cm, respectively. e quality of rock mass mainly reflects the pressure exerted on the tunnel lining, which is usually evaluated by the rock mass classification schemes. Over the past decades, a number of rock mass classification schemes have been proposed [39], in which the rock masses are classified into different grades according to their compactness and mechanical properties. For example, the RMR method divides the rock mass into five grades, i.e., grades I ∼IV, of which the grade I denotes the rock mass with the highest quality, while the grade V represents the rock mass with the lowest quality. Based on the data available in [10] and according to the design code of China [38], the K89 + 729 section is investigated herein and surrounding rock mass is classified as grade V by the RMR method. e main physical and mechanical properties of the grade V surrounding rock masses are listed in Table 3.

Input Parameters in the Analytical Solutions.
Although the exactly burial depth of the K89 + 729 section was not given in [10], the contact pressures (q) between the surrounding rock and the primary lining at the crown of the arch ring were presented, which equaled to 0.3 MPa. erefore, the equivalent burial depth (H) of the tunnel section can be back calculated based on Protodyakonov's theory, i.e., H � q/c � 9.5 m, which is smaller than the load equivalent height h q � 14.24 m. en, the detailed information about the ground condition and the secondary lining input in the analytical solutions can be derived as shown in Table 4.

FEM Numerical Modelling.
In order to accurately capture the key characteristics of the elastic foundation beam and the analytical solutions, the load-structure method is adopted in the 2D FEM model. e external pressures exerted on the double-arched tunnel can be calculated using the corresponding values in Table 4 based on the method introduced in Table 1. e determined external pressures will be imported into the load-structure method-based FEM simulations.
All the FEM analyses are carried out under the planestrain conditions. e numerical model is presented in Figure 8(b), in which the primary lining structure is not considered and the secondary lining is simulated by the beam elements according to the dimensions in Figure 8(a). e extrados and bottom of both side walls are restricted by unidirectional compression springs with a stiffness K � 200 MPa/m, which models the elastic resistances from the surrounding rock masses. e bottoms of the side and middle walls are also fixed in both horizontal and vertical directions. e rotation of the middle wall is also inhibited. It should be noted that the FEM model is specified to have the same loading configuration and assumptions of material properties as those adopted for analytical derivations. erefore, the dynamic ground-structure interaction and plastic yielding of lining structure are not considered. Figure 9 compares the internal forces derived from analytical solutions and numerical simulations. e variational trend and the distributive pattern of the internal forces of both methods are very close to each other. e discrepancies between the two are overall small and acceptable, except that the bending moment predicted by the analytical solutions deviates increasingly more from the numerical simulation data when approaching the middle wall. is may be because in the analytical solution the elastic resistance coefficient k d between the top of the middle wall and the arch ring is set to the elastic modulus E ′ of the lining structure in the analytical solutions. e simplification of beam element of the numerical model at the top of the middle wall may also deviate from the real situation, which may also account for the discrepancies. Further investigations may be required to determine the real contact characteristics at the top of the middle wall for double-arched tunnel.

Advances in Civil Engineering 13
Despite the discrepancies, the analytical solutions and the numerical solutions are overall in good agreement. It shows that the shear forces are smaller than the axial forces with the maximum shear force located at the conjunction point between the arch ring and the side wall, which is also an inflection point of the shear force distribution. Such an obvious inflection point may be because we only apply the elastic resistance to the side wall. e middle wall only bears the axial forces without presence of bending moment and shear forces in accordance with the symmetry condition of the lining structure and surrounding pressure during the derivation of the analytical solutions. e distribution law of internal forces can be derived from Figure 9. Specifically, the bending moment in the side wall is positive and increases from the bottom upwards but then decreases when approaching the arch foot. e maximum bending moment of the lining which is in the upper side along the side wall may be caused by the constraint to deformation provided by the external elastic resistance. is is confirmed by Figure 10, which shows that the upper side of the side wall experiences a larger outside deformation than the lower side. e maximum positive bending moment of the arch ring lies near the crown of the arch, and its magnitude is larger than that of the maximum negative bending moment, which is located at the point with an angle about 60°measured from arch crown to the half ring near both side walls. is may also be verified by the deformation shown in Figure 10, of which the larger deformation is present in the crown and outside haunch of both tunnels.
e results indicate that both the intrados of the crown and extrados of the haunch of the arch ring are in larger tension. ese characteristics should be noted during the design and construction of the tunnel lining structure. Generally, the axial force shows a relatively uniform distribution along the arch ring, while the axial force near the side wall is slightly larger than that near the crown of the arch ring. us, the side wall is in larger compression. e smaller axial force near the arch crown may be caused by the larger vertical pressure at the crown, which helps to reduce the axial force in the quasi-horizontal direction. e maximum values of shear force of arch rings in Figure 9(b) mainly appeared near the zero point of bending moment with the nonshear point located close to the place where the maximum bending moment is present, consistent with the general distribution law of forces within beam elements. Furthermore, the large shear force at the bottom of the side wall should be taken into consideration during the design and construction.

Parametric Analysis
During the derivation of the analytical solutions, most parameters can be deduced according to the mechanical properties and geometrical characteristics of the lining Advances in Civil Engineering structure. However, there are still some factors that cannot be derived directly, such as the real equivalent depth H and the elastic resistance coefficient K of the elastic foundation beam, which may affect the internal force distributions. Furthermore, the opening angle φ d of the double-arched tunnel is also a main geometrical parameter of the lining structure that may affect the internal force. Besides, to accommodate the increasing transportation demand, the mountain tunnels tend to have large dimensions, while the difference in tunnel span may also change the loading states of the lining structure. erefore, it is necessary to explore how these factors affect the internal forces and provide useful guidance for the design optimization. In this section, the four factors are varied systematically in the benchmark case presented in the section of Model Validation to investigate their influences on the internal forces. Due to the symmetry of the absolute values, only the internal forces at the left half of the double-arched tunnel are analysed in the following subsections. Note that the shear force on the right half is opposite in sign to that of the left half, while the bending moment and axial force are identical at both sides.

e Effect of the Burial Depth.
e equivalent depth H is varied between 5m and 10m at an increment of 1m.  16 Advances in Civil Engineering As shown in Figure 11, the largest positive and negative bending moment are in the side wall about l � 0.9 m and the arch ring at location about β � 125°, respectively. e maximum axial force is in the middle wall as shown in Figure 12. e largest shear force occurs at the bottom and top of the side wall as shown in Figure 13. e variations of the largest magnitudes of the internal forces with increasing H are shown in Figure 14, which all show linear correlations. Figure 14 indicates that the internal forces could be significantly affected by the burial depth as expected.

e Effect of Elastic Resistance Coefficient.
e elastic resistance coefficient K � 200 MPa/m was adopted in the analytical analysis and FEM analysis presented above for the grade V rock masses. In order to embrace the elastic resistance characteristics of surrounding rock masses with different grades, a wide range of K (K � 100, 200, 500, 1200, 1800, and 2800 MPa/m) is adopted in this section which may represent the grade I∼V rock masses. Figures 15-17 show the internal forces along the lining structure under different elastic resistance coefficients. e magnitude of bending moment along the arch ring decreases with the increasing K except the section near the top of the side wall which becomes larger as K increases. e bending moment along the side wall increases as K increases from 100 to 200 MPa/m. When K is greater than 500 MPa/m, the bending moment initially increases with the increasing K   decreases as K increases for β in the range of 0°to about 120°, while as β becomes greater than 120°, the magnitude of shear force increases with increasing K. e magnitude of shear force along the side wall decreases in the top three tenths of the side wall and increases for the remaining section for which K increases in the range of 100-200 MPa/m. When K is greater than 500 MPa/m, the magnitude of shear force decreases with increasing K for the top one tenth and bottom four tenths and increases for the remaining section of side wall. Overall, the elastic resistance coefficient has more influence on the internal forces along the side wall. is may be because the elastic resistance is only applied to the side wall but is not present over the arch ring. e situation when elastic resistance is present over the entire lining ring needs further exploration in the future.

e Effect of Opening
Angle. e opening angle φ d is one of the main geometrical parameters determining the dimension of the double-arched tunnel, which is measured from the arch crown to the top of the middle wall as shown in Figure 2. e analysis of the opening angle serves as an important reference for designing and optimizing the cross section of the double-arched tunnel. e arch ring radius R is assumed to be a constant. In such a case, the tunnel width B and height of middle wall h d would change in accordance with the change of φ d . e adopted values in the analysis are shown in Table 5. e results of the internal forces obtained for different cases are shown in  It can be found that for both side walls and middle wall, the bending moment and axial force increase as φ d increases.    Advances in Civil Engineering 23 conjunction region between the arch ring and the top of the middle wall. Since the conjunction region between the arch ring and the top of middle wall is a special location closely related to φ d where faults may exist and seepage may usually take place [40], it is of key concern for construction and design. erefore, the influence of φ d on the internal force at this location is explored in more details. Figure 21 shows the variation of the internal forces at the conjunction region

24
Advances in Civil Engineering between the arch ring and the top of the middle wall under different φ d . As can be seen from Figure 21, in the range of opening angles considered, the bending moment increases initially to about 630 kN m and then decreases as φ d further increases. In contrast, the shear force decreases initially to zero and reverses direction as φ d further increases. e axial force increases almost linearly with increasing φ d . Besides, Figure 19(a) shows that when φ d exceeds 55°, the axial force at the conjunction region will be the largest along the arch ring. Considering that the conjunction region is the most difficult part for construction and the quality of lining installation may be the lowest in comparison to other parts of the lining structure, it is practically safer to maintain the axial force at the conjunction region at a low level. erefore, it is suggested that φ d should not exceed 55°.

e effect of tunnel span.
e tunnel span of doublearched tunnel is represented by the total width B and is mainly decided by the radius of arch ring R and thickness of middle wall. In order to analyse the effect of tunnel span on the internal forces, the height of the side wall and middle wall and the opening angle in Figure 8(a) are not changed, and the burial depth remains 9.5 m, while the outer radius of arch ring increases from 5 m to 9 m at an increment of 1 m, and the outer radius 6.15 m in Figure 8(a) is still maintained. Neglecting the thickness of the middle wall in the analysis, the total width of the double-arched tunnel increases from 18.83 m to 33.89 m accordingly. e results of the internal forces obtained for different cases are shown in Figures 22-24.
A larger tunnel span means a larger area of tunnel lining subjected to the vertical pressure on the top of the tunnel. Figures 22-24 show that the variation trends of the internal forces with tunnel span are very similar to those with burial depth as shown in the Figures 11-13. As R increases, the absolute magnitude of axial force increases, while most of the absolute magnitude of the bending moment and shear force also increase except that near the inflection point. e point of the largest positive bending moment and the zero point of shear force of side wall move towards the top of side wall from about 1.2 m to 0.6 m as R increases. e largest negative bending moment is present in the arch ring at locations of β � 125°to β � 116°as R increases. ese positions are different from those in the analysis of burial depth H. e largest positive of bending moment is present in the side wall about I � 0.6 m to I � 1.2 m except that at R � 5 m, the largest positive of bending moment is located in the arch ring at location about β � 56°. e maximum axial force is in the middle wall as shown in Figure 23. e largest shear force occurs at the bottom and top of the side wall as shown in Figure 24. As shown in Figure 25, nonlinear correlations between the largest internal forces and R are observed, which are different from the linear correlations as between the largest internal forces and burial depth. is indicates that the tunnel span is more influential to the internal forces than the burial depth ( Figure 25).

Conclusions
is paper derived the analytical solutions for the internal forces of shallow double-arched tunnel based on the force method and elastic foundation beam model. e internal forces predicted by the proposed analytical solutions were in good agreement with the results of numerical simulations. e influences of different factors on the internal forces were investigated using the newly derived solutions. Based on the calculated and numerical simulation data, some conclusions are drawn as follows: (1) e variation trend of the internal forces in different burial depths is similar, but the absolute values increase as the burial depth increases. (2) e elastic resistance coefficient of the surrounding rock has a notable effect on the internal force of the side wall. In general, the internal forces decrease with increasing elastic resistance coefficient, which indicates a better performance of lining structure surrounded by rock masses with higher quality. Advances in Civil Engineering 25 (5) e results of the internal forces along the lining structure mainly show that the intrados of the crown and extrados of the haunch of the arch ring experience the largest tensile forces, while the side and middle walls experience larger compressive forces than other regions. ese zones require special attention during the construction and design.
e derived analytical solutions provide fast estimations for the internal forces and deformation of double-arched lining structures, which does not rely on either costly comprehensive field monitoring or time-consuming sophisticated numerical modelling. However, it should be noted that both the analytical solutions and numerical simulations conducted in this paper are based on the load-structure method but do not consider the dynamic interaction between tunnel structure and surrounding rock masses. Furthermore, it is likely that the double-arched tunnels will lie in the sloping region which induces unsymmetrical pressures. Both deserve further investigations in the future.

A. The Coefficients on the Top of Side Wall
e displacement coefficients at the conjunction region between the arch ring and the side wall can be obtained by solving the top coefficients of the side wall using the elastic foundation method. During the derivation of the displacement coefficients, the directions of the force and displacement are in accordance with the elastic foundation method, which are shown in Figure 26. After the displacement coefficients have been obtained, the direction of the displacement is changed to the specified direction in the derivation of the analytical solutions, in which the angular displacement is positive when the turning direction is outwards and the horizontal displacement is positive with outward movement. e equation set for the elastic foundation beam in Figure 26 is as follows: where y 0 , θ 0 , M 0 , and Q 0 are the initial parameters at the initial cross section of the elastic foundation beam; a � ����������� ((Kb)/(E ′ I))