# An Improved Numerical Model of Shield Tunnel with Double Lining and Its Applications

Authors: Qi-xiang Yan; Chao-fan Yao; Wen-bo Yang; Chuan He; Ping Geng
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

Based on the existing numerical models of shield tunnel with double lining, an improved numerical model is developed and its rationality is verified by a similar model test. In the improved numerical model, lining, joint, and junction surface, respectively, are simulated by beam, spring, and a combination of compression bar and spring. Through the comparison of the numerical analysis results of the improved model and existing models, it turns out that the defects or problems in the existing numerical models are resolved; tension appearance on the contact surface and junction surface and the abrupt change of bending moment in linings are solved in the improved model because the compression bar element and the coupling technology of node displacement in the junction surface is applied. Therefore, the improved numerical model could be applied to analyze double lining with waterproof on the junction surface and separation of the junction surface under an unfavorable load. In this paper, the parameter formulas of element stiffness, applicable to the junction surface and contact surface of double lining, are given definitely, and the influence of the element amount of junction surface on the analysis results is discussed. Based on the improved numerical model, the mechanical behavior of the double lining of the Huangpu River Tunnel in China is studied, and some conclusions are obtained as follows. (1) The thickness increase of the double lining will clearly increase its bending moment, but it has little influence on its axial force. (2) The allocation proportion of the bending moment between the segment lining and secondary lining has no linear relationship with the ratio of the lining thickness.


---

## Body

## 1. Introduction

Double lining is a kind of tunnel structure composed of an outer segment lining and an inner cast-in-place concrete secondary lining formed after tunnel boring machine installation of a segment lining, as shown in Figure1. As the double lining structure has advantages with respect to structure strengthening, anticollision, water proof strengthening, tunnel inner surface smoothing, tunnel S-shape correction, vibration absorption, and settlement resistance, many shield tunnels with double lining have appeared worldwide with the spread of tunnel boring machine technology. For example, the primary lining of a Dradenau sewer syphon tunnel in Hamburg is segment lining, and its secondary lining is cast-in-place concrete [1]. A tunnel whose external and internal diameters are 13.9 meters and 11.9 meters was constructed beneath a street carrying heavy traffic in the vicinity of Nishi-Shinjuku: its segment lining thickness is 0.7 meters and a secondary lining thickness is 0.3 meters [2]. The cracking in the extension shield tunnel of the blue metro line in Lisbon, Portugal, eventually led to corrosion of the rebar due to the marine environment, and the capacity of the tunnel would decrease in the long term, leading to a substantial risk of instability. For this reason, a secondary lining was constructed inside the existing tunnel as a durable remedial measure [3]. Shiziyang undersea tunnel of the Guangzhou-Shenzhen-Hong Kong express rail link in China has a segment lining of 0.5 meters and a secondary lining of 0.2 meters [4]. The Trans-Tokyo Bay highway tunnel, which is subjected to extremely soft strata under the sea, extremely high water pressure, and active seismic conditions in the Tokyo Bay area, has a large external diameter of the primary lining (approximately 14 m). To ensure the safety of the Trans-Tokyo Bay highway tunnel, a secondary lining inside with a thickness of 35 cm and segment lining with a thickness of 65 cm were applied [5]. As there are more and more underwater railway tunnels crossing rivers or seas in China, the application of double lining will become more frequent. The Huangpu River Tunnel in Shanghai-Nantong Railway, the Qiantang River Tunnel in Hangzhou-Changsha Railway, and the Qiongzhou Strait Tunnel from Zhanjiang to Haikou all will adopt the shield tunnel with double lining [6].Figure 1
Sketch of shield tunnel with double lining.The researches on shield tunnel mainly focus on single lining tunnel. Do et al. [7] applied a numerical study, which has been performed under seismic circumstance, to investigate the factors that affect segmental tunnel lining behavior. Gomes [8] examined the effect of the stress disturbance induced by tunnel construction on the completed tunnel’s seismic response. El Naggar et al. [9] developed an analytical procedure for evaluating in-plane moments and thrusts in composite and jointed tunnel linings during earthquakes.Some researchers also studied the behaviour of double lining structures. Murakami and Koizumi [10] investigated the behavior of shield tunnels reinforced by secondary lining, based on experiment results. Test models were composed of notched steel rings with mortar cast-in-place to their interiors. Based on the results, analytic models for double lining structures were proposed. Takamatsu et al. [11] used tests and analyses to investigate the behavior in the longitudinal direction of shield tunnels provided with secondary linings. Through studying a structural model for the evaluation of behavior in the longitudinal direction, confirmation of the validity of the structural model and a rational design method for shield tunnels with secondary linings were proposed. It was clarified that the member forces acting on the linings could be reduced by creating joints at appropriate positions on the secondary linings. Zhang et al. [12] suggested three joint interaction models corresponding to different forms of joint between primary and secondary linings and it was recommended that the coupling between primary and secondary linings must be reinforced, forming an integrated structure to bear the external pressure in the design of engineering.The existing double lining models have some defects: no analytic model for the above-mentioned design is so scientific that it is widely recognized and adopted. ITA’s (International Tunneling Association) Guidelines [13] for the design of shield tunnel lining (draft) do not involve the problems of double lining. Within China, there is usually one rough guiding principle for design: when there is enough shear resisting strength between the layers of double lining, it can be calculated as a whole structure and otherwise calculated as a composites plate only transmitting radial pressure. But no specific analysis model and calculating method are recommended for the design of double lining [14]. Therefore, on the basis of some existing double lining models, it is necessary to propose an improved double lining analysis model which can overcome some existing defects as far as possible. At present, the existing numerical models [10–12] of double lining mainly include Model I, Model II, and Model III, which are based on the two-dimensional beam-spring model and can be seen in Figure 2. The common characteristics of the three models are that the segment lining and secondary lining all adopt a beam element to simulate, the segment joint uses a spring element, and the contact surface between the segment lining and ground (hereinafter referred to as contact surface) employs radial spring element and tangential spring element to simulate. As all of these three models employ a radial spring element in the contact surface, radial tension will inevitably appear at the contact surface.Figure 2
Present numerical models of shield tunnel with double linings.The difference between these three models lies in the junction surface between the segment lining and the secondary lining (hereinafter referred to as the junction surface). In Model I, a radial beam element is used to simulate the transmission of the stress between the segment lining and secondary lining. Two nodes of the beam element are fixed in the beam element of the segment lining and secondary lining, respectively, which will lead to an abrupt change in the bending moment at the fixing point of the segment lining or secondary lining. Meanwhile, radial tension that does not conform to the fact will appear in this radial beam. In Model II, the radial spring was adopted in the junction surface, which could not avoid the appearance of radial tension in the junction surface as well. In Model III, a contact friction element was employed to simulate the force transfer on the junction surface. Its feature is that the tangential force value on junction surface is in direct proportion to the radial force and friction coefficient, which could not reflect the relationship between the tangential force and tangential relative displacement. Moreover, as there is no experimental data, the friction coefficient value of the junction surface is difficult to determine.In conclusion, the defects of the three models mentioned above are mainly as follows: the appearance of radial tension on the junction surface and the contact surface, an abrupt change in the bending moment, discordance of the radial force, and radial displacement of junction surface, in addition to evidence deficiency of friction coefficient determination. These defects would greatly limit the application of the models mentioned above in specific engineering. Aiming at these aforementioned defects, some improvements of the numerical model are conducted. A model of lining, joint, and junction surface, respectively, simulated by beam, spring, and combination of compression bar and spring is put forward to solve problems. Furthermore, the number of junction surface elements and their stiffness are discussed. Secondly, model tests are conducted to verify the improved model, and the comparison between the improved model and existing numerical models reveals the advantages of the improved model. Finally, the mechanical behavior of the Qiantang River Shield Tunnel with double lining in the Hangzhou-Changsha Railway in China is studied in terms of the improved model.


## 2. The Improved Numerical Model and Method

### 2.1. The Improved Numerical Model for Shield Tunnel with Double Lining

To overcome the defects of the aforementioned models, an improved model, whose lining, joint, and junction surface are, respectively, simulated by beam, spring, and a combination of compression bar and spring, is put forward, as shown in Figure3. The bar element bearing compression, which moves perpendicular with junction surface, is used only to simulate the radial force transmission on the contact surface and the junction surface. When the compression bar element is in tension, its axial force becomes zero naturally, so that the problem that radial tension appears on contact surface and junction surface can be solved. The coupled tangential spring, which runs parallel to the junction surface, could be used to extract the displacement of the segment lining node and its corresponding secondary lining node in the tangential direction to the two ends of the tangential spring. With extracted tangential displacement and calculated stiffness, the coupled tangential spring could simulate the tangential force transmission on the junction surface, so as to solve the problem that the tangential force and tangential displacement are unrelated to each other in Model III: see Figure 3(a). As the junction surface is simulated by the hinged radial compression bar and the tangential spring, both of which are unidirectional stress elements, it avoids the appearance of the abrupt change phenomenon of bending moment caused by the fixed beam element in Model I. The tension and compression spring, shear spring, and bend spring are set in the position of the segmental joint to simulate the passing of forces in every direction: see Figure 3(b). The comparison of the improved model and the existing three models can be seen in Table 1.Table 1
Comparison of elements between improved model and existing models.
Items
Segment lining and secondary lining
Segment lining joint
Junction surface
Contact surface
Model I
Beam element
Spring elements with functions of tension and compression, shear, and bending resistance
Fixed radial beam element
Radial and tangential spring elements with functions of tension and compression resistance
Model II
Beam element
Spring elements with functions of tension and compression, shear, and bending resistance
Radial and tangential spring elements with function of tension and compression resistance
Radial and tangential spring elements with functions of tension and compression resistance
Model III
Beam element
Spring elements with functions of tension and compression, shear, and bending resistance
Friction elements
Radial and tangential spring elements with functions of tension and compression resistance
Improved model
Beam element
Spring elements of node displacement coupling with functions of tension and compression, shear, and bending resistance
Radial compression bar with functions of compression resistance and tangential spring with coupling functions of elements node displacement
Radial compression bar with functions of compression resistance and tangential  spring elementFigure 3
Improved model.
(a)
Junction surface and contact surface
(b)
Segmental joint


### 2.2. Stiffness of Junction Surface and Contact Surface in Improved Model

The problem of the parameter value to choose relates to the lack of theory evidence in Model III. The improved model learns from the existing research results [15, 16] and gives the following calculation methods of the radial and tangential stiffness of the junction surface and the contact surface.(1) Stiffness of junction surface: after the segment lining is assembled, the secondary lining is manufactured by cast-in-place concrete directly. Double lining can be taken as a superimposition lining structure, and its radial compression can be simplified as the mechanical model shown as in Figure4. In the figure, B is the width of single segment lining in the longitudinal direction. Supposing that the compression stiffness between the layers is equal to the compression stiffness of the concrete column in the compression section, which has half thickness of the interior lining and external lining, the compression stiffness in the compression section of the interior lining and the external lining is as follows:(1a)
K
Y
1
=
E
1
×
A
t
1
′
=
2
E
1
×
A
t
1
,
(1b)
K
Y
2
=
E
2
×
A
t
2
′
=
2
E
2
×
A
t
2
.In the formulas, K
Y
1 and K
Y
2 represent the compression stiffness of the segment lining and secondary lining, respectively, kN/m; E
1 and E
2 represent the elasticity modules of concrete in the segment lining and secondary lining, respectively, kN/m2; t
1 and t
2 represent the thickness of the segment lining and secondary lining, respectively, m; t
1
′ and t
2
′ represent the thickness of the compression zone in the segment lining and secondary lining, respectively, m; A is the corresponding area of single radial compression bar element in the junction surface: A
=
B
×
θ
×
R
H, m2; B is the width of the single segment lining in the longitudinal direction, m; θ is the corresponding calculating angle of the single compression bar element in the junction surface; R
H is the radius of the junction surface, m.Figure 4
Compression stiffness calculation schematic diagram of junction surface.
(a)
Cross section
(b)
A-A sectionCombine Formulas (1a) and (1b), according to the calculating formula of spring in series stiffness coefficient: that is,(2)
K
Y
=
K
Y
1
×
K
Y
2
K
Y
1
+
K
Y
2
=
2
E
1
E
2
A
E
1
t
2
+
E
2
t
1
.In the formula, K
Y is the radial compression stiffness of a single compression bar in the junction surface, kN/m.The junction surface tangential relationship is similar to the radial relationship, and the tangential stiffness of the double lining is as follows:(3a)
K
T
1
=
G
1
×
A
t
1
′
=
2
G
1
×
A
t
1
,
(3b)
K
T
2
=
G
2
×
A
t
2
′
=
2
G
2
×
A
t
2
,
(4)
K
T
=
K
T
1
×
K
T
2
K
T
1
+
K
T
2
=
2
G
1
G
2
A
G
1
t
2
+
G
2
t
1
.In the formula, K
T is the tangential stiffness of a single tangential spring in the junction surface, kN/m; K
T
1 and K
T
2 represent the tangential shear stiffness of the segment lining and secondary lining, respectively, kN/m; G
1 and G
2 represent the shear modulus of concrete in the segment lining and secondary lining, respectively, kN/m2.(2) Stiffness of contact surface: the calculating model of the contact surface is shown as in Figure5. As the segment lining stiffness is far larger than that of stratum, it can be considered that the segment lining stiffness has no influence on the stiffness of the contact surface approximately. The radial stiffness of a single compression bar at the contact surface can be presented as(5)
K
R
=
k
r
×
A
′
.In the formula, K
R is the radial stiffness of single compression bar, kN/m; k
r is the radial resistance coefficient of stratum, kN/m3; A
′ is the area corresponding to a single radial compression bar element at the contact surface: A
′
=
B
×
θ
′
×
R
w, m2; θ
′ is the calculating angle of a single radial compression bar element at the contact surface; R
w is the external radius of the segment lining.Figure 5
Schematic diagram of calculation model of contact face between segment and strata.
(a)
Mechanical model of contact surface
(b)
Calculation areaAccording to the experience of Japan, the contact surface tangential stiffness is approximately half the radial stiffness [5].


### 2.3. Comparison between Numerical Simulation and Experimental Test

To verify the rationality of the improved model and compare the advantages and defects of the existing numerical models, the shield tunnel-ground complex simulation equipment is employed to conduct the model experiment on the Shiziyang Tunnel in China with double lining. Without water pressure, the test results with the numerical analysis results of the aforementioned Model I and Model II in the same load and same structure parameters are compared.


#### 2.3.1. Experimental Subject

The double lining section of the Shiziyang Tunnel in Guangzhou-Shenzhen railway, which is buried at the depth of 35.35 m, is taken as the prototype structure. The external and internal radiuses of the tunnel structure are 10.8 m and 9.4 m, and the thicknesses of the segment lining and the secondary lining are 0.5 m and 0.2 m. The width of the single segment lining in the longitudinal direction is 2.0 m. The segment lining is assembled with a staggered joint, and “7
+
1” block divided mode is employed, which can be seen in Figure 6. The central angle of the key block is 16.37°, and the central angle of the 2 adjacent blocks and 5 standard blocks is 49.09°. The calculation section of this tunnel mainly crosses argillaceous siltstone stratum.Figure 6
Double lining of Shiziyang Tunnel.


#### 2.3.2. Experimental Device

The tunnel-ground complex simulation equipment is a plane strain loading device of a shield tunnel, which mainly includes three parts: a hydraulic loading unit, stratum and structural body, and signal collection transaction and analysis apparatus. The whole test platform’s length, height, and width are 6.0 m × 6.0 m × 2.55 m, and the tunnel model and model soil around it are put inside the device. The available filling range of the model soil is 3.0 m × 3.0 m × 0.2 m, as shown in Figure6. The loading panel is set in the vertical direction of the tunnel crossing section (direction Z), and the vertical pressure increased to the loading panel through the jack to guarantee the tunnel at the plane strain state. Horizon directions X and Y are set appositional to the jack, respectively, and make some impact on the removable baffle. Through the force transfer device the needed load is added to direction X (simulated horizon soil pressure) and direction Y (simulated vertical soil pressure). The schematic diagram of the experiment device and physical device is shown in Figures 7 and 8.Figure 7
Sketch of shield tunnel-ground complex simulation equipment (unit: mm).
(a)
Plan
(b)
ElevationFigure 8
Tunnel-ground complex simulation equipment.


#### 2.3.3. Method

The test takes the 1 : 20 geometric similarity ratio and 1 : 1 bulk density similarity ratio as the basic similarity ratios. According to the elastic similarity theory, Poisson’s ratio and strain and the internal friction angle similarity ratio can be deduced asC
μ
=
C
ε
=
C
ϕ
=
1, as well as the strength, stress, cohesion, and elastic modulus similarity ratio C
R
=
C
σ
=
C
c
=
C
E
=
20. The model soil is filled around the double lining model and passes horizontal jack output to the double lining model through the model soil. The similarity model satisfies the similarity relationship mainly in the following aspects: the ground cohesion, compression modulus, internal friction angle and bulk density, the uniaxial compressive strength of the segment lining, double lining and elastic modulus, and the equivalent tensile (compression) strength of main reinforcement.The model stratum employs a mixture made of barite powder, fly ash, river sand, coarse quartz, fine quartz, rosin, and oil with the specific ratio, and the mass ratio is 1 : 0.341 : 0.676 : 0.334 : 0.334 : 0.08 : 0.14. The mechanical properties of the argillaceous siltstone stratum prototype and its model material can be seen in Table2.Table 2
Parameters of argillaceous siltstone strata in model and prototype.
Name
Cohesionc/kPa
Internal friction angleϕ/(°)
Model elastic quantityE/MPa
Bulk densityγ/(kN·m−3)
The prototype stratum
22.0
23.5°
28
22.3
The model stratum
1.1
23.5°
1.4
22.3
Similarity ratio
1 : 20
1 : 1
1 : 20
1 : 1The material of the segment lining is C50 reinforced concrete, and the concrete model adopts the composite material made of water, gypsum, and diatomite with the ratio of 1 : 1.40 : 0.1, and the ring direction main reinforcement employs the fine iron with the diameter ofφ1.2 mm to model. The material of the secondary lining is C30 plain concrete, which uses the composite material made of water and gypsum with the ratio of 1 : 1.40, and is added through the cast-in-place concrete. Slits are cut with specific depth in the joint location to weaken the flexural stiffness, which is a way to simulate the transverse joints, and the slits’ depth is calculated on the basis of the equivalent principal of the prototype joints’ flexural ability, which can be seen in Table 3. As there is tiny dislocation in the longitudinal joint parts of the segment in the practical engineering, steel wire is used with specific diameter in the section of the longitudinal joint linking to simulate the longitudinal joint. Steel wire with the diameter of 4 mm and length of 40 mm is employed in this test. The segment lining model and double lining after casting the secondary lining can be seen in Figure 9.Table 3
Depths of slits in test models.
Bending stiffness/(N·m·rad−1)
Equivalent slit depth in prototype/mm
Equivalent slit depth in prototype of model/mm
2.0 × 107 (positive bending)
240
12 (tunnel crown/tunnel bottom)
2.6 × 107 (negative bending)
200
10 (left/right side lumbers)Figure 9
Test model of double lining.
(a)
Segment lining
(b)
Segment lining and secondary liningThe electric resistance strain gauges are laid out in a transverse direction inside and outside the segment lining and secondary lining, and the test employs the electric resistance strain indicator to test the values of the gauges. The structural internal force of the segment lining and secondary lining could be calculated based on the values of gauges.Based on the calculating results of the section, the test takes the vertical soil pressure of the section (the loading directionY of the test) as the uniform soil pressure, and the value is 285.3 kN/m; the horizontal earth pressure (the loading direction X of the test) is also taken as the uniform soil pressure, and the value is 112.5 kN/m. The external load of Model I, Model II, and the improved model is in accord with the test, and the lining beam element amount of these two kinds of numerical models is 400. According to Formulas (1a) and (1b) to (5), the radial compression bar and tangential spring element amount of the junction surface are both 400, and their stiffness is 7.28 × 106 kN/m and 3.03 × 106 kN/m, respectively. The radial compression bar and tangential spring stiffness of the ground’s contact surface are 4.41 × 105 kN/m/n and 2.21 × 105 kN/m/n, respectively, and n is the tangential compression bar amount of the contact surface.


## 2.1. The Improved Numerical Model for Shield Tunnel with Double Lining

To overcome the defects of the aforementioned models, an improved model, whose lining, joint, and junction surface are, respectively, simulated by beam, spring, and a combination of compression bar and spring, is put forward, as shown in Figure3. The bar element bearing compression, which moves perpendicular with junction surface, is used only to simulate the radial force transmission on the contact surface and the junction surface. When the compression bar element is in tension, its axial force becomes zero naturally, so that the problem that radial tension appears on contact surface and junction surface can be solved. The coupled tangential spring, which runs parallel to the junction surface, could be used to extract the displacement of the segment lining node and its corresponding secondary lining node in the tangential direction to the two ends of the tangential spring. With extracted tangential displacement and calculated stiffness, the coupled tangential spring could simulate the tangential force transmission on the junction surface, so as to solve the problem that the tangential force and tangential displacement are unrelated to each other in Model III: see Figure 3(a). As the junction surface is simulated by the hinged radial compression bar and the tangential spring, both of which are unidirectional stress elements, it avoids the appearance of the abrupt change phenomenon of bending moment caused by the fixed beam element in Model I. The tension and compression spring, shear spring, and bend spring are set in the position of the segmental joint to simulate the passing of forces in every direction: see Figure 3(b). The comparison of the improved model and the existing three models can be seen in Table 1.Table 1
Comparison of elements between improved model and existing models.
Items
Segment lining and secondary lining
Segment lining joint
Junction surface
Contact surface
Model I
Beam element
Spring elements with functions of tension and compression, shear, and bending resistance
Fixed radial beam element
Radial and tangential spring elements with functions of tension and compression resistance
Model II
Beam element
Spring elements with functions of tension and compression, shear, and bending resistance
Radial and tangential spring elements with function of tension and compression resistance
Radial and tangential spring elements with functions of tension and compression resistance
Model III
Beam element
Spring elements with functions of tension and compression, shear, and bending resistance
Friction elements
Radial and tangential spring elements with functions of tension and compression resistance
Improved model
Beam element
Spring elements of node displacement coupling with functions of tension and compression, shear, and bending resistance
Radial compression bar with functions of compression resistance and tangential spring with coupling functions of elements node displacement
Radial compression bar with functions of compression resistance and tangential  spring elementFigure 3
Improved model.
(a)
Junction surface and contact surface
(b)
Segmental joint


## 2.2. Stiffness of Junction Surface and Contact Surface in Improved Model

The problem of the parameter value to choose relates to the lack of theory evidence in Model III. The improved model learns from the existing research results [15, 16] and gives the following calculation methods of the radial and tangential stiffness of the junction surface and the contact surface.(1) Stiffness of junction surface: after the segment lining is assembled, the secondary lining is manufactured by cast-in-place concrete directly. Double lining can be taken as a superimposition lining structure, and its radial compression can be simplified as the mechanical model shown as in Figure4. In the figure, B is the width of single segment lining in the longitudinal direction. Supposing that the compression stiffness between the layers is equal to the compression stiffness of the concrete column in the compression section, which has half thickness of the interior lining and external lining, the compression stiffness in the compression section of the interior lining and the external lining is as follows:(1a)
K
Y
1
=
E
1
×
A
t
1
′
=
2
E
1
×
A
t
1
,
(1b)
K
Y
2
=
E
2
×
A
t
2
′
=
2
E
2
×
A
t
2
.In the formulas, K
Y
1 and K
Y
2 represent the compression stiffness of the segment lining and secondary lining, respectively, kN/m; E
1 and E
2 represent the elasticity modules of concrete in the segment lining and secondary lining, respectively, kN/m2; t
1 and t
2 represent the thickness of the segment lining and secondary lining, respectively, m; t
1
′ and t
2
′ represent the thickness of the compression zone in the segment lining and secondary lining, respectively, m; A is the corresponding area of single radial compression bar element in the junction surface: A
=
B
×
θ
×
R
H, m2; B is the width of the single segment lining in the longitudinal direction, m; θ is the corresponding calculating angle of the single compression bar element in the junction surface; R
H is the radius of the junction surface, m.Figure 4
Compression stiffness calculation schematic diagram of junction surface.
(a)
Cross section
(b)
A-A sectionCombine Formulas (1a) and (1b), according to the calculating formula of spring in series stiffness coefficient: that is,(2)
K
Y
=
K
Y
1
×
K
Y
2
K
Y
1
+
K
Y
2
=
2
E
1
E
2
A
E
1
t
2
+
E
2
t
1
.In the formula, K
Y is the radial compression stiffness of a single compression bar in the junction surface, kN/m.The junction surface tangential relationship is similar to the radial relationship, and the tangential stiffness of the double lining is as follows:(3a)
K
T
1
=
G
1
×
A
t
1
′
=
2
G
1
×
A
t
1
,
(3b)
K
T
2
=
G
2
×
A
t
2
′
=
2
G
2
×
A
t
2
,
(4)
K
T
=
K
T
1
×
K
T
2
K
T
1
+
K
T
2
=
2
G
1
G
2
A
G
1
t
2
+
G
2
t
1
.In the formula, K
T is the tangential stiffness of a single tangential spring in the junction surface, kN/m; K
T
1 and K
T
2 represent the tangential shear stiffness of the segment lining and secondary lining, respectively, kN/m; G
1 and G
2 represent the shear modulus of concrete in the segment lining and secondary lining, respectively, kN/m2.(2) Stiffness of contact surface: the calculating model of the contact surface is shown as in Figure5. As the segment lining stiffness is far larger than that of stratum, it can be considered that the segment lining stiffness has no influence on the stiffness of the contact surface approximately. The radial stiffness of a single compression bar at the contact surface can be presented as(5)
K
R
=
k
r
×
A
′
.In the formula, K
R is the radial stiffness of single compression bar, kN/m; k
r is the radial resistance coefficient of stratum, kN/m3; A
′ is the area corresponding to a single radial compression bar element at the contact surface: A
′
=
B
×
θ
′
×
R
w, m2; θ
′ is the calculating angle of a single radial compression bar element at the contact surface; R
w is the external radius of the segment lining.Figure 5
Schematic diagram of calculation model of contact face between segment and strata.
(a)
Mechanical model of contact surface
(b)
Calculation areaAccording to the experience of Japan, the contact surface tangential stiffness is approximately half the radial stiffness [5].


## 2.3. Comparison between Numerical Simulation and Experimental Test

To verify the rationality of the improved model and compare the advantages and defects of the existing numerical models, the shield tunnel-ground complex simulation equipment is employed to conduct the model experiment on the Shiziyang Tunnel in China with double lining. Without water pressure, the test results with the numerical analysis results of the aforementioned Model I and Model II in the same load and same structure parameters are compared.


### 2.3.1. Experimental Subject

The double lining section of the Shiziyang Tunnel in Guangzhou-Shenzhen railway, which is buried at the depth of 35.35 m, is taken as the prototype structure. The external and internal radiuses of the tunnel structure are 10.8 m and 9.4 m, and the thicknesses of the segment lining and the secondary lining are 0.5 m and 0.2 m. The width of the single segment lining in the longitudinal direction is 2.0 m. The segment lining is assembled with a staggered joint, and “7
+
1” block divided mode is employed, which can be seen in Figure 6. The central angle of the key block is 16.37°, and the central angle of the 2 adjacent blocks and 5 standard blocks is 49.09°. The calculation section of this tunnel mainly crosses argillaceous siltstone stratum.Figure 6
Double lining of Shiziyang Tunnel.


### 2.3.2. Experimental Device

The tunnel-ground complex simulation equipment is a plane strain loading device of a shield tunnel, which mainly includes three parts: a hydraulic loading unit, stratum and structural body, and signal collection transaction and analysis apparatus. The whole test platform’s length, height, and width are 6.0 m × 6.0 m × 2.55 m, and the tunnel model and model soil around it are put inside the device. The available filling range of the model soil is 3.0 m × 3.0 m × 0.2 m, as shown in Figure6. The loading panel is set in the vertical direction of the tunnel crossing section (direction Z), and the vertical pressure increased to the loading panel through the jack to guarantee the tunnel at the plane strain state. Horizon directions X and Y are set appositional to the jack, respectively, and make some impact on the removable baffle. Through the force transfer device the needed load is added to direction X (simulated horizon soil pressure) and direction Y (simulated vertical soil pressure). The schematic diagram of the experiment device and physical device is shown in Figures 7 and 8.Figure 7
Sketch of shield tunnel-ground complex simulation equipment (unit: mm).
(a)
Plan
(b)
ElevationFigure 8
Tunnel-ground complex simulation equipment.


### 2.3.3. Method

The test takes the 1 : 20 geometric similarity ratio and 1 : 1 bulk density similarity ratio as the basic similarity ratios. According to the elastic similarity theory, Poisson’s ratio and strain and the internal friction angle similarity ratio can be deduced asC
μ
=
C
ε
=
C
ϕ
=
1, as well as the strength, stress, cohesion, and elastic modulus similarity ratio C
R
=
C
σ
=
C
c
=
C
E
=
20. The model soil is filled around the double lining model and passes horizontal jack output to the double lining model through the model soil. The similarity model satisfies the similarity relationship mainly in the following aspects: the ground cohesion, compression modulus, internal friction angle and bulk density, the uniaxial compressive strength of the segment lining, double lining and elastic modulus, and the equivalent tensile (compression) strength of main reinforcement.The model stratum employs a mixture made of barite powder, fly ash, river sand, coarse quartz, fine quartz, rosin, and oil with the specific ratio, and the mass ratio is 1 : 0.341 : 0.676 : 0.334 : 0.334 : 0.08 : 0.14. The mechanical properties of the argillaceous siltstone stratum prototype and its model material can be seen in Table2.Table 2
Parameters of argillaceous siltstone strata in model and prototype.
Name
Cohesionc/kPa
Internal friction angleϕ/(°)
Model elastic quantityE/MPa
Bulk densityγ/(kN·m−3)
The prototype stratum
22.0
23.5°
28
22.3
The model stratum
1.1
23.5°
1.4
22.3
Similarity ratio
1 : 20
1 : 1
1 : 20
1 : 1The material of the segment lining is C50 reinforced concrete, and the concrete model adopts the composite material made of water, gypsum, and diatomite with the ratio of 1 : 1.40 : 0.1, and the ring direction main reinforcement employs the fine iron with the diameter ofφ1.2 mm to model. The material of the secondary lining is C30 plain concrete, which uses the composite material made of water and gypsum with the ratio of 1 : 1.40, and is added through the cast-in-place concrete. Slits are cut with specific depth in the joint location to weaken the flexural stiffness, which is a way to simulate the transverse joints, and the slits’ depth is calculated on the basis of the equivalent principal of the prototype joints’ flexural ability, which can be seen in Table 3. As there is tiny dislocation in the longitudinal joint parts of the segment in the practical engineering, steel wire is used with specific diameter in the section of the longitudinal joint linking to simulate the longitudinal joint. Steel wire with the diameter of 4 mm and length of 40 mm is employed in this test. The segment lining model and double lining after casting the secondary lining can be seen in Figure 9.Table 3
Depths of slits in test models.
Bending stiffness/(N·m·rad−1)
Equivalent slit depth in prototype/mm
Equivalent slit depth in prototype of model/mm
2.0 × 107 (positive bending)
240
12 (tunnel crown/tunnel bottom)
2.6 × 107 (negative bending)
200
10 (left/right side lumbers)Figure 9
Test model of double lining.
(a)
Segment lining
(b)
Segment lining and secondary liningThe electric resistance strain gauges are laid out in a transverse direction inside and outside the segment lining and secondary lining, and the test employs the electric resistance strain indicator to test the values of the gauges. The structural internal force of the segment lining and secondary lining could be calculated based on the values of gauges.Based on the calculating results of the section, the test takes the vertical soil pressure of the section (the loading directionY of the test) as the uniform soil pressure, and the value is 285.3 kN/m; the horizontal earth pressure (the loading direction X of the test) is also taken as the uniform soil pressure, and the value is 112.5 kN/m. The external load of Model I, Model II, and the improved model is in accord with the test, and the lining beam element amount of these two kinds of numerical models is 400. According to Formulas (1a) and (1b) to (5), the radial compression bar and tangential spring element amount of the junction surface are both 400, and their stiffness is 7.28 × 106 kN/m and 3.03 × 106 kN/m, respectively. The radial compression bar and tangential spring stiffness of the ground’s contact surface are 4.41 × 105 kN/m/n and 2.21 × 105 kN/m/n, respectively, and n is the tangential compression bar amount of the contact surface.


## 2.3.1. Experimental Subject

The double lining section of the Shiziyang Tunnel in Guangzhou-Shenzhen railway, which is buried at the depth of 35.35 m, is taken as the prototype structure. The external and internal radiuses of the tunnel structure are 10.8 m and 9.4 m, and the thicknesses of the segment lining and the secondary lining are 0.5 m and 0.2 m. The width of the single segment lining in the longitudinal direction is 2.0 m. The segment lining is assembled with a staggered joint, and “7
+
1” block divided mode is employed, which can be seen in Figure 6. The central angle of the key block is 16.37°, and the central angle of the 2 adjacent blocks and 5 standard blocks is 49.09°. The calculation section of this tunnel mainly crosses argillaceous siltstone stratum.Figure 6
Double lining of Shiziyang Tunnel.


## 2.3.2. Experimental Device

The tunnel-ground complex simulation equipment is a plane strain loading device of a shield tunnel, which mainly includes three parts: a hydraulic loading unit, stratum and structural body, and signal collection transaction and analysis apparatus. The whole test platform’s length, height, and width are 6.0 m × 6.0 m × 2.55 m, and the tunnel model and model soil around it are put inside the device. The available filling range of the model soil is 3.0 m × 3.0 m × 0.2 m, as shown in Figure6. The loading panel is set in the vertical direction of the tunnel crossing section (direction Z), and the vertical pressure increased to the loading panel through the jack to guarantee the tunnel at the plane strain state. Horizon directions X and Y are set appositional to the jack, respectively, and make some impact on the removable baffle. Through the force transfer device the needed load is added to direction X (simulated horizon soil pressure) and direction Y (simulated vertical soil pressure). The schematic diagram of the experiment device and physical device is shown in Figures 7 and 8.Figure 7
Sketch of shield tunnel-ground complex simulation equipment (unit: mm).
(a)
Plan
(b)
ElevationFigure 8
Tunnel-ground complex simulation equipment.


## 2.3.3. Method

The test takes the 1 : 20 geometric similarity ratio and 1 : 1 bulk density similarity ratio as the basic similarity ratios. According to the elastic similarity theory, Poisson’s ratio and strain and the internal friction angle similarity ratio can be deduced asC
μ
=
C
ε
=
C
ϕ
=
1, as well as the strength, stress, cohesion, and elastic modulus similarity ratio C
R
=
C
σ
=
C
c
=
C
E
=
20. The model soil is filled around the double lining model and passes horizontal jack output to the double lining model through the model soil. The similarity model satisfies the similarity relationship mainly in the following aspects: the ground cohesion, compression modulus, internal friction angle and bulk density, the uniaxial compressive strength of the segment lining, double lining and elastic modulus, and the equivalent tensile (compression) strength of main reinforcement.The model stratum employs a mixture made of barite powder, fly ash, river sand, coarse quartz, fine quartz, rosin, and oil with the specific ratio, and the mass ratio is 1 : 0.341 : 0.676 : 0.334 : 0.334 : 0.08 : 0.14. The mechanical properties of the argillaceous siltstone stratum prototype and its model material can be seen in Table2.Table 2
Parameters of argillaceous siltstone strata in model and prototype.
Name
Cohesionc/kPa
Internal friction angleϕ/(°)
Model elastic quantityE/MPa
Bulk densityγ/(kN·m−3)
The prototype stratum
22.0
23.5°
28
22.3
The model stratum
1.1
23.5°
1.4
22.3
Similarity ratio
1 : 20
1 : 1
1 : 20
1 : 1The material of the segment lining is C50 reinforced concrete, and the concrete model adopts the composite material made of water, gypsum, and diatomite with the ratio of 1 : 1.40 : 0.1, and the ring direction main reinforcement employs the fine iron with the diameter ofφ1.2 mm to model. The material of the secondary lining is C30 plain concrete, which uses the composite material made of water and gypsum with the ratio of 1 : 1.40, and is added through the cast-in-place concrete. Slits are cut with specific depth in the joint location to weaken the flexural stiffness, which is a way to simulate the transverse joints, and the slits’ depth is calculated on the basis of the equivalent principal of the prototype joints’ flexural ability, which can be seen in Table 3. As there is tiny dislocation in the longitudinal joint parts of the segment in the practical engineering, steel wire is used with specific diameter in the section of the longitudinal joint linking to simulate the longitudinal joint. Steel wire with the diameter of 4 mm and length of 40 mm is employed in this test. The segment lining model and double lining after casting the secondary lining can be seen in Figure 9.Table 3
Depths of slits in test models.
Bending stiffness/(N·m·rad−1)
Equivalent slit depth in prototype/mm
Equivalent slit depth in prototype of model/mm
2.0 × 107 (positive bending)
240
12 (tunnel crown/tunnel bottom)
2.6 × 107 (negative bending)
200
10 (left/right side lumbers)Figure 9
Test model of double lining.
(a)
Segment lining
(b)
Segment lining and secondary liningThe electric resistance strain gauges are laid out in a transverse direction inside and outside the segment lining and secondary lining, and the test employs the electric resistance strain indicator to test the values of the gauges. The structural internal force of the segment lining and secondary lining could be calculated based on the values of gauges.Based on the calculating results of the section, the test takes the vertical soil pressure of the section (the loading directionY of the test) as the uniform soil pressure, and the value is 285.3 kN/m; the horizontal earth pressure (the loading direction X of the test) is also taken as the uniform soil pressure, and the value is 112.5 kN/m. The external load of Model I, Model II, and the improved model is in accord with the test, and the lining beam element amount of these two kinds of numerical models is 400. According to Formulas (1a) and (1b) to (5), the radial compression bar and tangential spring element amount of the junction surface are both 400, and their stiffness is 7.28 × 106 kN/m and 3.03 × 106 kN/m, respectively. The radial compression bar and tangential spring stiffness of the ground’s contact surface are 4.41 × 105 kN/m/n and 2.21 × 105 kN/m/n, respectively, and n is the tangential compression bar amount of the contact surface.


## 3. Analysis of Results and Discussion

### 3.1. Results after Comparison

As the friction coefficient value of the junction surface of Model III is difficult to determine, it will only compare the result of test with the results of Model I, Model II, and the improved model here. The segment lining joint stiffness of Model I, Model II, and the improved model is in accord with that of test prototype. The measured internal force of the model test and the calculating internal force of the numerical models are shown in Figure10, and the extreme internal force values can be seen in Table 4. The bending moment takes the external tensile as the positive value and the axial force takes the compression as the positive value.Table 4
Comparison on extreme internal forces between numerical and test results.
Subject
Model I
Mode II
The improved model
Similarity model test
Maximum value
Minimum value
Maximum value
Minimum value
Maximum value
Minimum value
Maximum value
Minimum value
Bending moment (kN·m)
Segment lining
585.65
−639.55
603.44
−636.72
627.35
−641.52
663.74
−635.67
Secondary lining
163.47
−144.32
195.42
−192.54
178.48
−162.88
188.54
−148.87
Axial force (MN)
Segment lining
−10.98
−4.65
−12.03
−5.97
−11.67
−4.75
−10.28
−5.24
Secondary lining
−2.89
−1.45
−2.32
−0.92
−2.78
−1.22
−2.44
−1.07Figure 10
Diagrams for tested values and calculated values of internal forces in lining.
(a)
Bending moment of segmental lining (unit: kN·m)
(b)
Bending moment of secondary lining (unit: kN·m)
(c)
Axial force of segmental lining (unit: MN)
(d)
Axial force of secondary lining (unit: MN)From Figure10 it can be seen that the bending moment distribution and axial force distribution of the segment lining and secondary lining in the model test and improved numerical model are similar, and the extreme values of bending moment and axial force on the segment lining and secondary lining are relatively close as well. For instance, the segment lining bending moment’s maximum values from the similarity model test and improved numerical model are 663.74 kN·m and 627.35 kN·m, respectively, and the minimum values are −635.67 kN·m and −641.52 kN·m, respectively. The secondary lining maximum values from the model test and the improved numerical model are 188.54 kN·m and 178.48 kN·m, respectively, and the minimum values are −148.87 kN·m and −162.88 kN·m, respectively. Corresponding to the similarity model test, the relative error of the bending moment’s extreme value is within ±10.0%. According to the same analysis, the relative error of the axial force’s extreme value is within ±12.5%. It indicates that using the improved numerical model to analyze the mechanical behavior of double lining is feasible. The bending moment and axial force from the improved numerical model and similarity model test have some difference as well. The reason for this difference may be due to the fact that the improved numerical model is a model of load application to structure directly, but the model test is a model of load application to structure through stratum. Even though the imposed external load between the two is the same, the actual bearing load of the structure in the improved model has some difference with that of the structure in the model test. Meanwhile, adopting the imposed constraint from the radial and tangential ground spring will make some errors to the imposed constraint from similar soil material to the structure in the test.According to Figure10, it also can be seen that the internal force distribution of Model I and Model II is close to the improved model too. The reason is that Model I and Model II employ the same stiffness calculating methods and parameters, which match the improved model’s junction surface and the contact surface. But in Figure 10(a), from the bending moment’s enlarged drawing in the parts of the crown and haunch of the segment lining, it can be established that there still exist differences of internal force distribution among Model I, Model II, and the improved model. The bending moment of numerical Model I is presented with the zigzag shape in detailed parts of the crown, and the bending moment of numerical Model II has some relatively large fluctuations in the crown and haunch parts as well. And there is precisely the abrupt change of bending moment caused by the radial beam element of the junction surface in Model I and the radial spring transmitting tension in Model II, although the phenomenon should not really exist. When laying the waterproof layer between the segment lining and the secondary lining, which cannot pass tension completely, carrying on the mechanical behavior analysis will cause more obvious error when using Model I and Model II.When there is a relatively large difference between the vertical and horizontal loads, the deformation of segment lining and secondary lining will no longer be synchronous in some parts, where the junction surface will separate inevitably. And laying the waterproof layer between the segment lining and the secondary lining, the junction surface could not transfer tension. Under these conditions, the values from Model I and Model II are removed from reality. At this time, using the improved model will avoid this phenomenon effectively and acquire the internal force distribution much closer to the model test. Therefore, the improved model should be selected as the numerical analysis model of the double lining because of its greater rationality.


### 3.2. Application and Discussion of Improved Model

To know more about the mechanical behavior of double lining, taking the Huangpu River Tunnel in the Shanghai-Nantong Railway in China as the example here, the improved numerical model is adopted to simulate and analyse. The total length of the Huangpu River Tunnel is 8.3 km, and the internal and external diameters are 8.74 m and 10.3 m. The segment lining employs a reinforced concrete structure of the “5
+
2
+
1” segment block form and universal segment lining with a thickness of 0.48 m. The secondary lining is made up of a cast-in-place reinforced concrete structure with a thickness of 0.3 m. The basic material parameters of the double lining can be seen in Figure 5.The double lining structure schematic is shown as in Figure11(a), and the load distribution is shown in Figure 11(b). The earth pressure adopts the Terzaghi Theory to calculate: p
1 and p
2 are the pressures of the covering soil and reaction force of the bottom ground, respectively, q
1 and q
2 are the lateral earth pressures of tunnel crown and tunnel bottom, respectively, and w
1 and w
2 are the water pressures of the tunnel crown and tunnel bottom, respectively. A section with a buried depth of 20.51 m and the water head of 19.75 m is chosen to analyse. At this time, p
1, p
2, q
1, q
2, w
1, and w
2 are 880.7 kN/m, 638.1 kN/m, 442.5 kN/m, 1162.3 kN/m, 2051.3 kN/m, and 3081.1 kN/m, respectively.Figure 11
Schematic diagram of double linings and its load case.
(a)
Sketch of shield tunnel with double lining
(b)
Load conditionThe thickness of the segment lining and the secondary lining in the models is calculated according to its actual thickness. The amount of the junction surface compression bar, the amount of the junction surface tangential spring, the stratum compression bar amount, and the stratum tangential spring amount all are 400. Putting the parameters in Table5 into Formula (2), it could be calculated that the radial compression bar stiffness of the junction surface is 6.78 × 106 kN/m and the tangential spring stiffness is 2.82 × 106 kN/m. According to the geological data, the ground radial resistance coefficient is 7 MPa/m3. According to (5), the radial compression bar stiffness of the contact surface can be calculated: 2.26 × 105 kN/m/n, where n is the amount of the ground contact surface radial compression bar, taking half of the radial stiffness as the stiffness of tangential spring.Table 5
Material parameters of double linings.
Subject
Material
Elastic modulus (kPa)
Density (kg/m3)
Poisson’s ratio
Segment lining
C50 concrete
3.45 × 107
2450
0.2
Secondary lining
C30 concrete
3.00 × 107
2450
0.2According to the relevant research results [17, 18], the stiffness of the tension and compression spring of the segment joint may be taken as 4.14 × 109 kN/m and the shear spring stiffness as 1.34 × 109 kN/m. The transversal flexural stiffness of the segment joint may be taken as 4 × 105 kN·m when traction is in the internal surface of the lining and 2.4 × 105 kN·m when traction is in the outside surface.


#### 3.2.1. The Discussion of Junction Surface Element Amount

Theoretically, the radial compression bar and tangential spring amount of the junction surface (i.e., the element amount of the junction surface, the same below) influences the stiffness of the single radial compression bar and tangential spring. Obviously, the larger the amount, the smaller the area corresponding to the single radial compression bar and tangential spring, and the corresponding stiffness is smaller. Different junction surface quantities must be set in order to optimize and analyse, so as to determine the reasonable amount of tangential spring and radial compression bar. Setting the radial compression bar and tangential spring amount between the segment lining and the secondary lining as 20, 40, 100, 200, 400, 1000, and 2000, respectively, and also comparing the maximum deformation of the segment lining and the secondary lining, the extreme value of the bending moment (the bending moment is positive when traction is in the outside surface, the opposite is negative, the maximum is the positive peak, and the minimum is the negative peak, the same below) and the extreme value of the axial force (the maximum value of the axial force is the maximum pressure value, and the minimum value is the minimum pressure value, the same below) can be noted. The calculated results can be seen in Figure12.Figure 12
The influence of element number at junction surface on inner forces of double lining.
(a)
The max deformation
(b)
The extreme bending moment
(c)
The extreme thrust forceFigure12 shows that, with the constant increase of the junction surface element amount, the maximum deformation value of the secondary lining and segment lining will eventually trend towards convergence. And it will be basically steady when the element amount of the junction surface is more than 400; when the junction surface amount is less than 100, the lining bending moment extreme value will be more influenced by the element amount. But when the element amount is more than 400, the bending moment extreme values of the secondary lining or of the segment lining will both achieve basic stability. Meanwhile, when the junction surface element amount is less than 400, the axial force extreme value of the segment lining will increase with the element amount, and the axial force extreme value of secondary lining will decrease gradually with the element amount. But when the element amount is more than 400, the axial force extreme value of these two will both be basically steady. The aforementioned phenomenon shows that the junction surface element amount could not be chosen optionally because the amount will have a large influence on the analysis result. And among the different engineering objects, their junction surface element amount has a reasonable range.As far as the double lining of the Huangpu River Shield Tunnel is concerned, comprehensively considering all the factors including the maximum deformation value, the bending moment extreme value, and axial force extreme value of the segment lining and the double lining, a conclusion can be drawn that it is appropriate to take 400 as the junction surface amount. The discussion below will be based on the condition of 400 radial compression bars and 400 tangential springs.


#### 3.2.2. Mechanical Behavior of Double Lining

Based on the radial compression bar and tangential spring amount of the junction surface above, and with a focus on the Huangpu River Tunnel, the mechanical behavior analysis of the different segment lining thickness and different double lining thickness within the same external load is carried out.(1) Different thickness of segment lining: the thickness of the secondary lining is 0.30 m, and the thicknesses of the segment lining are 0.30 m, 0.40 m, 0.50 m, 0.60 m, and 0.70 m, respectively, and the thicknesses of the corresponding double lining are 0.60 m, 0.70 m, 0.80 m, 0.90 m, and 1.00 m, respectively. The bending moment extreme values of the double lining (the sum of the bending moment of the segment lining and the secondary lining in the same section within a single ring breadth with 2 m, the same below) and the corresponding bending moment of the segment lining can be seen in Figure13(a) and Table 6. The maximum and minimum axial force of the double lining and segment lining (the sum of axial force of the segment lining and secondary lining in the same section within a single ring breadth with 2 m, the same below) and the corresponding axial force of the segment lining can be seen in Figure 13(b) and Table 7.Table 6
The bending moments extreme values under different thicknesses of segment lining.
Thickness of segment lining (cm)
30
40
50
60
70
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum bending moment of the double lining
Bending moment of double lining (kN·m)
131.047
142.037
148.896
151.945
152.454
Bending moment of segment lining at the corresponding section (kN·m)
63.178
97.604
120.650
133.701
140.300
The allocation proportion of bending moment of segment lining
48%
69%
81%
88%
92%
The increase ratio of bending moment of double lining
100%
108%
114%
116%
116%
The section of minimum bending moment of the double lining
Bending moment of double lining (kN·m)
−158.954
−171.705
−179.542
−182.704
−183.285
Bending moment of segment lining at the corresponding section (kN·m)
−77.292
−118.510
−145.771
−161.090
−168.780
The allocation proportion of bending moment of segment lining
49%
69%
81%
88%
92%
The increase ratio of bending moment of double lining
100%
108%
113%
115%
115%Table 7
The axial force extreme values under different thicknesses of segment lining.
Thickness of segment lining (cm)
30
40
50
60
70
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum axial force of the double lining
Axial force of double lining (MN)
−1.144
−1.142
−1.135
−1.124
−1.120
Axial force of segment lining at the corresponding section (MN)
−0.744
−0.789
−0.819
−0.845
−0.875
The allocation proportion of axial force of segment lining
65%
69%
72%
75%
78%
The increase ratio of axial force of double lining
100%
99%
99%
98%
98%
The section of minimum axial force of the double lining
Axial force of double lining (MN)
−1.030
−1.031
−1.022
−1.021
−1.012
Axial force of segment lining at the corresponding section (MN)
−0.669
−0.715
−0.744
−0.778
−0.797
The allocation proportion of axial force of segment lining
65%
69%
73%
76%
79%
The increase ratio of axial force of double lining
100%
99%
99%
99%
98%Figure 13
Numerical results of double lining under different thicknesses of segment lining.
(a)
The bending moment
(b)
The axial forceFigure13(a) shows that, on condition that the secondary lining has a constant or specific thickness, with the segment lining thickness increasing constantly, the absolute values of the maximum and minimum bending moment of the double lining both tend to increase, and the absolute values of bending moments for the segment lining at each corresponding section increase as well. From Figure 13(b), it can be seen that the maximum and minimum axial force of the double lining has little change with the increase of the segment lining thickness, but each axial force of the segment lining at a corresponding section obviously increases.From Tables6 and 7, it can be seen that, with the constant increase of the segment lining thickness, the increasing ratio of the bending moment and axial force for the double lining has little relation to the increasing ratio of the double lining thickness. The allocation proportion of the bending moment and axial force between the segment lining and secondary lining has no linear relationship with the ratio of the segment lining and the secondary lining thickness. With the constant increase in the segment lining thickness, the ratio of bending moment and axial force burdened by the segment lining goes up gradually. The allocation proportion of the bending moment of the segment lining at each section of the maximum and the minimum bending moment of the double lining takes a basically consistent ratio, and the allocation proportion of the axial force at each section of the maximum and minimum axial force shows a similar law as well.(2) Different thickness of secondary lining: the thickness of segment lining is 0.50 m, and the thicknesses of the secondary lining are 0.10 m, 0.20 m, 0.30 m, 0.40 m, and 0.50 m, respectively, and the thicknesses of the corresponding double lining are 0.60 m, 0.70 m, 0.80 m, 0.90 m, and 1.00 m, respectively. Through computing, the maximum and minimum bending moment of double lining and the bending moment of the segment lining at the corresponding section can be seen in Figure14(a) and Table 8, and the maximum and minimum axial force of the double lining and the axial force of segment lining at the corresponding section can be seen in Figure 14(b) and Table 9.Table 8
The bending moments extreme values under different thicknesses of secondary lining.
Thickness of secondary lining (cm)
10
20
30
40
50
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum bending moment of the double lining
Bending moment of double lining (kN·m)
146.689
147.398
148.896
150.660
151.945
Bending moment of segment lining at the corresponding section (kN·m)
145.480
138.001
120.650
94.474
72.271
The allocation proportion of bending moment of segment lining
99%
94%
81%
62%
48%
The increase ratio of bending moment of double lining
100%
100%
102%
103%
104%
The section of minimum bending moment of the double lining
Bending moment of double lining (kN·m)
−177.101
−177.886
−179.542
−181.469
−182.814
Bending moment of segment lining at the corresponding section (kN·m)
−175.650
−166.640
−145.771
−116.662
−87.462
The allocation proportion of bending moment of segment lining
99%
94%
81%
64%
48%
The increase ratio of bending moment of double lining
100%
100%
101%
102%
103%Table 9
The axial force extreme values under different thicknesses of secondary lining.
Thickness of secondary lining (cm)
10
20
30
40
50
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum axial force of the double lining
Axial force of double lining (MN)
−1.121
−1.121
−1.121
−1.120
−1.120
Axial force of segment lining at the corresponding section (MN)
−1.009
−0.897
−0.818
−0.762
−0.728
The allocation proportion of axial force of segment lining
90%
80%
73%
68%
65%
The increase ratio of axial force of double lining
100%
100%
100%
100%
100%
The section of minimum axial force of the double lining
Axial force of double lining (MN)
−1.058
−1.057
−1.057
−1.056
−1.056
Axial force of segment lining at the corresponding section (MN)
−0.963
−0.856
−0.772
−0.718
−0.686
The allocation proportion of axial force of segment lining
91%
81%
73%
68%
65%
The increase ratio of axial force of double lining
100%
100%
100%
100%
100%Figure 14
Numerical results of double lining under different thicknesses of secondary lining.
(a)
The bending moment
(b)
The axial forceFrom Figure14(a) it can be seen that, on condition that the secondary lining has a specific thickness, with the secondary lining thickness increasing constantly, the absolute values of the maximum and minimum bending moment of the double lining both tend to increase, but the rate of increase is small. And the absolute value of the maximum and minimum bending moment of the segment lining at the corresponding section tends to decrease. From Figure 14(b), it can be seen that the maximum and minimum axial force of the double lining barely changes with the increase in the secondary lining thickness, but each axial force extreme value of the segment lining at the corresponding section tends to decrease.From Tables8 and 9 it can be seen that, with the constant increase of the thickness of the secondary lining, the increasing ratio of bending moment and axial force of the double lining has little relationship with the increasing ratio of the double lining thickness, and the allocation proportion of the bending moment and the axial force between the segment lining and the secondary lining has no linear relationship with the ratio of the segment lining and the secondary lining thickness. With the constant increase in the secondary thickness, the ratio of the bending moment and the axial force that is burdened by the secondary lining decreases quickly. The allocation proportion of the bending moment of the segment lining at each section of the maximum bending moment and the minimum bending moment of the double lining takes a basically consistent ratio, and the same law applies to the axial force as well.


## 3.1. Results after Comparison

As the friction coefficient value of the junction surface of Model III is difficult to determine, it will only compare the result of test with the results of Model I, Model II, and the improved model here. The segment lining joint stiffness of Model I, Model II, and the improved model is in accord with that of test prototype. The measured internal force of the model test and the calculating internal force of the numerical models are shown in Figure10, and the extreme internal force values can be seen in Table 4. The bending moment takes the external tensile as the positive value and the axial force takes the compression as the positive value.Table 4
Comparison on extreme internal forces between numerical and test results.
Subject
Model I
Mode II
The improved model
Similarity model test
Maximum value
Minimum value
Maximum value
Minimum value
Maximum value
Minimum value
Maximum value
Minimum value
Bending moment (kN·m)
Segment lining
585.65
−639.55
603.44
−636.72
627.35
−641.52
663.74
−635.67
Secondary lining
163.47
−144.32
195.42
−192.54
178.48
−162.88
188.54
−148.87
Axial force (MN)
Segment lining
−10.98
−4.65
−12.03
−5.97
−11.67
−4.75
−10.28
−5.24
Secondary lining
−2.89
−1.45
−2.32
−0.92
−2.78
−1.22
−2.44
−1.07Figure 10
Diagrams for tested values and calculated values of internal forces in lining.
(a)
Bending moment of segmental lining (unit: kN·m)
(b)
Bending moment of secondary lining (unit: kN·m)
(c)
Axial force of segmental lining (unit: MN)
(d)
Axial force of secondary lining (unit: MN)From Figure10 it can be seen that the bending moment distribution and axial force distribution of the segment lining and secondary lining in the model test and improved numerical model are similar, and the extreme values of bending moment and axial force on the segment lining and secondary lining are relatively close as well. For instance, the segment lining bending moment’s maximum values from the similarity model test and improved numerical model are 663.74 kN·m and 627.35 kN·m, respectively, and the minimum values are −635.67 kN·m and −641.52 kN·m, respectively. The secondary lining maximum values from the model test and the improved numerical model are 188.54 kN·m and 178.48 kN·m, respectively, and the minimum values are −148.87 kN·m and −162.88 kN·m, respectively. Corresponding to the similarity model test, the relative error of the bending moment’s extreme value is within ±10.0%. According to the same analysis, the relative error of the axial force’s extreme value is within ±12.5%. It indicates that using the improved numerical model to analyze the mechanical behavior of double lining is feasible. The bending moment and axial force from the improved numerical model and similarity model test have some difference as well. The reason for this difference may be due to the fact that the improved numerical model is a model of load application to structure directly, but the model test is a model of load application to structure through stratum. Even though the imposed external load between the two is the same, the actual bearing load of the structure in the improved model has some difference with that of the structure in the model test. Meanwhile, adopting the imposed constraint from the radial and tangential ground spring will make some errors to the imposed constraint from similar soil material to the structure in the test.According to Figure10, it also can be seen that the internal force distribution of Model I and Model II is close to the improved model too. The reason is that Model I and Model II employ the same stiffness calculating methods and parameters, which match the improved model’s junction surface and the contact surface. But in Figure 10(a), from the bending moment’s enlarged drawing in the parts of the crown and haunch of the segment lining, it can be established that there still exist differences of internal force distribution among Model I, Model II, and the improved model. The bending moment of numerical Model I is presented with the zigzag shape in detailed parts of the crown, and the bending moment of numerical Model II has some relatively large fluctuations in the crown and haunch parts as well. And there is precisely the abrupt change of bending moment caused by the radial beam element of the junction surface in Model I and the radial spring transmitting tension in Model II, although the phenomenon should not really exist. When laying the waterproof layer between the segment lining and the secondary lining, which cannot pass tension completely, carrying on the mechanical behavior analysis will cause more obvious error when using Model I and Model II.When there is a relatively large difference between the vertical and horizontal loads, the deformation of segment lining and secondary lining will no longer be synchronous in some parts, where the junction surface will separate inevitably. And laying the waterproof layer between the segment lining and the secondary lining, the junction surface could not transfer tension. Under these conditions, the values from Model I and Model II are removed from reality. At this time, using the improved model will avoid this phenomenon effectively and acquire the internal force distribution much closer to the model test. Therefore, the improved model should be selected as the numerical analysis model of the double lining because of its greater rationality.


## 3.2. Application and Discussion of Improved Model

To know more about the mechanical behavior of double lining, taking the Huangpu River Tunnel in the Shanghai-Nantong Railway in China as the example here, the improved numerical model is adopted to simulate and analyse. The total length of the Huangpu River Tunnel is 8.3 km, and the internal and external diameters are 8.74 m and 10.3 m. The segment lining employs a reinforced concrete structure of the “5
+
2
+
1” segment block form and universal segment lining with a thickness of 0.48 m. The secondary lining is made up of a cast-in-place reinforced concrete structure with a thickness of 0.3 m. The basic material parameters of the double lining can be seen in Figure 5.The double lining structure schematic is shown as in Figure11(a), and the load distribution is shown in Figure 11(b). The earth pressure adopts the Terzaghi Theory to calculate: p
1 and p
2 are the pressures of the covering soil and reaction force of the bottom ground, respectively, q
1 and q
2 are the lateral earth pressures of tunnel crown and tunnel bottom, respectively, and w
1 and w
2 are the water pressures of the tunnel crown and tunnel bottom, respectively. A section with a buried depth of 20.51 m and the water head of 19.75 m is chosen to analyse. At this time, p
1, p
2, q
1, q
2, w
1, and w
2 are 880.7 kN/m, 638.1 kN/m, 442.5 kN/m, 1162.3 kN/m, 2051.3 kN/m, and 3081.1 kN/m, respectively.Figure 11
Schematic diagram of double linings and its load case.
(a)
Sketch of shield tunnel with double lining
(b)
Load conditionThe thickness of the segment lining and the secondary lining in the models is calculated according to its actual thickness. The amount of the junction surface compression bar, the amount of the junction surface tangential spring, the stratum compression bar amount, and the stratum tangential spring amount all are 400. Putting the parameters in Table5 into Formula (2), it could be calculated that the radial compression bar stiffness of the junction surface is 6.78 × 106 kN/m and the tangential spring stiffness is 2.82 × 106 kN/m. According to the geological data, the ground radial resistance coefficient is 7 MPa/m3. According to (5), the radial compression bar stiffness of the contact surface can be calculated: 2.26 × 105 kN/m/n, where n is the amount of the ground contact surface radial compression bar, taking half of the radial stiffness as the stiffness of tangential spring.Table 5
Material parameters of double linings.
Subject
Material
Elastic modulus (kPa)
Density (kg/m3)
Poisson’s ratio
Segment lining
C50 concrete
3.45 × 107
2450
0.2
Secondary lining
C30 concrete
3.00 × 107
2450
0.2According to the relevant research results [17, 18], the stiffness of the tension and compression spring of the segment joint may be taken as 4.14 × 109 kN/m and the shear spring stiffness as 1.34 × 109 kN/m. The transversal flexural stiffness of the segment joint may be taken as 4 × 105 kN·m when traction is in the internal surface of the lining and 2.4 × 105 kN·m when traction is in the outside surface.


### 3.2.1. The Discussion of Junction Surface Element Amount

Theoretically, the radial compression bar and tangential spring amount of the junction surface (i.e., the element amount of the junction surface, the same below) influences the stiffness of the single radial compression bar and tangential spring. Obviously, the larger the amount, the smaller the area corresponding to the single radial compression bar and tangential spring, and the corresponding stiffness is smaller. Different junction surface quantities must be set in order to optimize and analyse, so as to determine the reasonable amount of tangential spring and radial compression bar. Setting the radial compression bar and tangential spring amount between the segment lining and the secondary lining as 20, 40, 100, 200, 400, 1000, and 2000, respectively, and also comparing the maximum deformation of the segment lining and the secondary lining, the extreme value of the bending moment (the bending moment is positive when traction is in the outside surface, the opposite is negative, the maximum is the positive peak, and the minimum is the negative peak, the same below) and the extreme value of the axial force (the maximum value of the axial force is the maximum pressure value, and the minimum value is the minimum pressure value, the same below) can be noted. The calculated results can be seen in Figure12.Figure 12
The influence of element number at junction surface on inner forces of double lining.
(a)
The max deformation
(b)
The extreme bending moment
(c)
The extreme thrust forceFigure12 shows that, with the constant increase of the junction surface element amount, the maximum deformation value of the secondary lining and segment lining will eventually trend towards convergence. And it will be basically steady when the element amount of the junction surface is more than 400; when the junction surface amount is less than 100, the lining bending moment extreme value will be more influenced by the element amount. But when the element amount is more than 400, the bending moment extreme values of the secondary lining or of the segment lining will both achieve basic stability. Meanwhile, when the junction surface element amount is less than 400, the axial force extreme value of the segment lining will increase with the element amount, and the axial force extreme value of secondary lining will decrease gradually with the element amount. But when the element amount is more than 400, the axial force extreme value of these two will both be basically steady. The aforementioned phenomenon shows that the junction surface element amount could not be chosen optionally because the amount will have a large influence on the analysis result. And among the different engineering objects, their junction surface element amount has a reasonable range.As far as the double lining of the Huangpu River Shield Tunnel is concerned, comprehensively considering all the factors including the maximum deformation value, the bending moment extreme value, and axial force extreme value of the segment lining and the double lining, a conclusion can be drawn that it is appropriate to take 400 as the junction surface amount. The discussion below will be based on the condition of 400 radial compression bars and 400 tangential springs.


### 3.2.2. Mechanical Behavior of Double Lining

Based on the radial compression bar and tangential spring amount of the junction surface above, and with a focus on the Huangpu River Tunnel, the mechanical behavior analysis of the different segment lining thickness and different double lining thickness within the same external load is carried out.(1) Different thickness of segment lining: the thickness of the secondary lining is 0.30 m, and the thicknesses of the segment lining are 0.30 m, 0.40 m, 0.50 m, 0.60 m, and 0.70 m, respectively, and the thicknesses of the corresponding double lining are 0.60 m, 0.70 m, 0.80 m, 0.90 m, and 1.00 m, respectively. The bending moment extreme values of the double lining (the sum of the bending moment of the segment lining and the secondary lining in the same section within a single ring breadth with 2 m, the same below) and the corresponding bending moment of the segment lining can be seen in Figure13(a) and Table 6. The maximum and minimum axial force of the double lining and segment lining (the sum of axial force of the segment lining and secondary lining in the same section within a single ring breadth with 2 m, the same below) and the corresponding axial force of the segment lining can be seen in Figure 13(b) and Table 7.Table 6
The bending moments extreme values under different thicknesses of segment lining.
Thickness of segment lining (cm)
30
40
50
60
70
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum bending moment of the double lining
Bending moment of double lining (kN·m)
131.047
142.037
148.896
151.945
152.454
Bending moment of segment lining at the corresponding section (kN·m)
63.178
97.604
120.650
133.701
140.300
The allocation proportion of bending moment of segment lining
48%
69%
81%
88%
92%
The increase ratio of bending moment of double lining
100%
108%
114%
116%
116%
The section of minimum bending moment of the double lining
Bending moment of double lining (kN·m)
−158.954
−171.705
−179.542
−182.704
−183.285
Bending moment of segment lining at the corresponding section (kN·m)
−77.292
−118.510
−145.771
−161.090
−168.780
The allocation proportion of bending moment of segment lining
49%
69%
81%
88%
92%
The increase ratio of bending moment of double lining
100%
108%
113%
115%
115%Table 7
The axial force extreme values under different thicknesses of segment lining.
Thickness of segment lining (cm)
30
40
50
60
70
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum axial force of the double lining
Axial force of double lining (MN)
−1.144
−1.142
−1.135
−1.124
−1.120
Axial force of segment lining at the corresponding section (MN)
−0.744
−0.789
−0.819
−0.845
−0.875
The allocation proportion of axial force of segment lining
65%
69%
72%
75%
78%
The increase ratio of axial force of double lining
100%
99%
99%
98%
98%
The section of minimum axial force of the double lining
Axial force of double lining (MN)
−1.030
−1.031
−1.022
−1.021
−1.012
Axial force of segment lining at the corresponding section (MN)
−0.669
−0.715
−0.744
−0.778
−0.797
The allocation proportion of axial force of segment lining
65%
69%
73%
76%
79%
The increase ratio of axial force of double lining
100%
99%
99%
99%
98%Figure 13
Numerical results of double lining under different thicknesses of segment lining.
(a)
The bending moment
(b)
The axial forceFigure13(a) shows that, on condition that the secondary lining has a constant or specific thickness, with the segment lining thickness increasing constantly, the absolute values of the maximum and minimum bending moment of the double lining both tend to increase, and the absolute values of bending moments for the segment lining at each corresponding section increase as well. From Figure 13(b), it can be seen that the maximum and minimum axial force of the double lining has little change with the increase of the segment lining thickness, but each axial force of the segment lining at a corresponding section obviously increases.From Tables6 and 7, it can be seen that, with the constant increase of the segment lining thickness, the increasing ratio of the bending moment and axial force for the double lining has little relation to the increasing ratio of the double lining thickness. The allocation proportion of the bending moment and axial force between the segment lining and secondary lining has no linear relationship with the ratio of the segment lining and the secondary lining thickness. With the constant increase in the segment lining thickness, the ratio of bending moment and axial force burdened by the segment lining goes up gradually. The allocation proportion of the bending moment of the segment lining at each section of the maximum and the minimum bending moment of the double lining takes a basically consistent ratio, and the allocation proportion of the axial force at each section of the maximum and minimum axial force shows a similar law as well.(2) Different thickness of secondary lining: the thickness of segment lining is 0.50 m, and the thicknesses of the secondary lining are 0.10 m, 0.20 m, 0.30 m, 0.40 m, and 0.50 m, respectively, and the thicknesses of the corresponding double lining are 0.60 m, 0.70 m, 0.80 m, 0.90 m, and 1.00 m, respectively. Through computing, the maximum and minimum bending moment of double lining and the bending moment of the segment lining at the corresponding section can be seen in Figure14(a) and Table 8, and the maximum and minimum axial force of the double lining and the axial force of segment lining at the corresponding section can be seen in Figure 14(b) and Table 9.Table 8
The bending moments extreme values under different thicknesses of secondary lining.
Thickness of secondary lining (cm)
10
20
30
40
50
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum bending moment of the double lining
Bending moment of double lining (kN·m)
146.689
147.398
148.896
150.660
151.945
Bending moment of segment lining at the corresponding section (kN·m)
145.480
138.001
120.650
94.474
72.271
The allocation proportion of bending moment of segment lining
99%
94%
81%
62%
48%
The increase ratio of bending moment of double lining
100%
100%
102%
103%
104%
The section of minimum bending moment of the double lining
Bending moment of double lining (kN·m)
−177.101
−177.886
−179.542
−181.469
−182.814
Bending moment of segment lining at the corresponding section (kN·m)
−175.650
−166.640
−145.771
−116.662
−87.462
The allocation proportion of bending moment of segment lining
99%
94%
81%
64%
48%
The increase ratio of bending moment of double lining
100%
100%
101%
102%
103%Table 9
The axial force extreme values under different thicknesses of secondary lining.
Thickness of secondary lining (cm)
10
20
30
40
50
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum axial force of the double lining
Axial force of double lining (MN)
−1.121
−1.121
−1.121
−1.120
−1.120
Axial force of segment lining at the corresponding section (MN)
−1.009
−0.897
−0.818
−0.762
−0.728
The allocation proportion of axial force of segment lining
90%
80%
73%
68%
65%
The increase ratio of axial force of double lining
100%
100%
100%
100%
100%
The section of minimum axial force of the double lining
Axial force of double lining (MN)
−1.058
−1.057
−1.057
−1.056
−1.056
Axial force of segment lining at the corresponding section (MN)
−0.963
−0.856
−0.772
−0.718
−0.686
The allocation proportion of axial force of segment lining
91%
81%
73%
68%
65%
The increase ratio of axial force of double lining
100%
100%
100%
100%
100%Figure 14
Numerical results of double lining under different thicknesses of secondary lining.
(a)
The bending moment
(b)
The axial forceFrom Figure14(a) it can be seen that, on condition that the secondary lining has a specific thickness, with the secondary lining thickness increasing constantly, the absolute values of the maximum and minimum bending moment of the double lining both tend to increase, but the rate of increase is small. And the absolute value of the maximum and minimum bending moment of the segment lining at the corresponding section tends to decrease. From Figure 14(b), it can be seen that the maximum and minimum axial force of the double lining barely changes with the increase in the secondary lining thickness, but each axial force extreme value of the segment lining at the corresponding section tends to decrease.From Tables8 and 9 it can be seen that, with the constant increase of the thickness of the secondary lining, the increasing ratio of bending moment and axial force of the double lining has little relationship with the increasing ratio of the double lining thickness, and the allocation proportion of the bending moment and the axial force between the segment lining and the secondary lining has no linear relationship with the ratio of the segment lining and the secondary lining thickness. With the constant increase in the secondary thickness, the ratio of the bending moment and the axial force that is burdened by the secondary lining decreases quickly. The allocation proportion of the bending moment of the segment lining at each section of the maximum bending moment and the minimum bending moment of the double lining takes a basically consistent ratio, and the same law applies to the axial force as well.


## 3.2.1. The Discussion of Junction Surface Element Amount

Theoretically, the radial compression bar and tangential spring amount of the junction surface (i.e., the element amount of the junction surface, the same below) influences the stiffness of the single radial compression bar and tangential spring. Obviously, the larger the amount, the smaller the area corresponding to the single radial compression bar and tangential spring, and the corresponding stiffness is smaller. Different junction surface quantities must be set in order to optimize and analyse, so as to determine the reasonable amount of tangential spring and radial compression bar. Setting the radial compression bar and tangential spring amount between the segment lining and the secondary lining as 20, 40, 100, 200, 400, 1000, and 2000, respectively, and also comparing the maximum deformation of the segment lining and the secondary lining, the extreme value of the bending moment (the bending moment is positive when traction is in the outside surface, the opposite is negative, the maximum is the positive peak, and the minimum is the negative peak, the same below) and the extreme value of the axial force (the maximum value of the axial force is the maximum pressure value, and the minimum value is the minimum pressure value, the same below) can be noted. The calculated results can be seen in Figure12.Figure 12
The influence of element number at junction surface on inner forces of double lining.
(a)
The max deformation
(b)
The extreme bending moment
(c)
The extreme thrust forceFigure12 shows that, with the constant increase of the junction surface element amount, the maximum deformation value of the secondary lining and segment lining will eventually trend towards convergence. And it will be basically steady when the element amount of the junction surface is more than 400; when the junction surface amount is less than 100, the lining bending moment extreme value will be more influenced by the element amount. But when the element amount is more than 400, the bending moment extreme values of the secondary lining or of the segment lining will both achieve basic stability. Meanwhile, when the junction surface element amount is less than 400, the axial force extreme value of the segment lining will increase with the element amount, and the axial force extreme value of secondary lining will decrease gradually with the element amount. But when the element amount is more than 400, the axial force extreme value of these two will both be basically steady. The aforementioned phenomenon shows that the junction surface element amount could not be chosen optionally because the amount will have a large influence on the analysis result. And among the different engineering objects, their junction surface element amount has a reasonable range.As far as the double lining of the Huangpu River Shield Tunnel is concerned, comprehensively considering all the factors including the maximum deformation value, the bending moment extreme value, and axial force extreme value of the segment lining and the double lining, a conclusion can be drawn that it is appropriate to take 400 as the junction surface amount. The discussion below will be based on the condition of 400 radial compression bars and 400 tangential springs.


## 3.2.2. Mechanical Behavior of Double Lining

Based on the radial compression bar and tangential spring amount of the junction surface above, and with a focus on the Huangpu River Tunnel, the mechanical behavior analysis of the different segment lining thickness and different double lining thickness within the same external load is carried out.(1) Different thickness of segment lining: the thickness of the secondary lining is 0.30 m, and the thicknesses of the segment lining are 0.30 m, 0.40 m, 0.50 m, 0.60 m, and 0.70 m, respectively, and the thicknesses of the corresponding double lining are 0.60 m, 0.70 m, 0.80 m, 0.90 m, and 1.00 m, respectively. The bending moment extreme values of the double lining (the sum of the bending moment of the segment lining and the secondary lining in the same section within a single ring breadth with 2 m, the same below) and the corresponding bending moment of the segment lining can be seen in Figure13(a) and Table 6. The maximum and minimum axial force of the double lining and segment lining (the sum of axial force of the segment lining and secondary lining in the same section within a single ring breadth with 2 m, the same below) and the corresponding axial force of the segment lining can be seen in Figure 13(b) and Table 7.Table 6
The bending moments extreme values under different thicknesses of segment lining.
Thickness of segment lining (cm)
30
40
50
60
70
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum bending moment of the double lining
Bending moment of double lining (kN·m)
131.047
142.037
148.896
151.945
152.454
Bending moment of segment lining at the corresponding section (kN·m)
63.178
97.604
120.650
133.701
140.300
The allocation proportion of bending moment of segment lining
48%
69%
81%
88%
92%
The increase ratio of bending moment of double lining
100%
108%
114%
116%
116%
The section of minimum bending moment of the double lining
Bending moment of double lining (kN·m)
−158.954
−171.705
−179.542
−182.704
−183.285
Bending moment of segment lining at the corresponding section (kN·m)
−77.292
−118.510
−145.771
−161.090
−168.780
The allocation proportion of bending moment of segment lining
49%
69%
81%
88%
92%
The increase ratio of bending moment of double lining
100%
108%
113%
115%
115%Table 7
The axial force extreme values under different thicknesses of segment lining.
Thickness of segment lining (cm)
30
40
50
60
70
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum axial force of the double lining
Axial force of double lining (MN)
−1.144
−1.142
−1.135
−1.124
−1.120
Axial force of segment lining at the corresponding section (MN)
−0.744
−0.789
−0.819
−0.845
−0.875
The allocation proportion of axial force of segment lining
65%
69%
72%
75%
78%
The increase ratio of axial force of double lining
100%
99%
99%
98%
98%
The section of minimum axial force of the double lining
Axial force of double lining (MN)
−1.030
−1.031
−1.022
−1.021
−1.012
Axial force of segment lining at the corresponding section (MN)
−0.669
−0.715
−0.744
−0.778
−0.797
The allocation proportion of axial force of segment lining
65%
69%
73%
76%
79%
The increase ratio of axial force of double lining
100%
99%
99%
99%
98%Figure 13
Numerical results of double lining under different thicknesses of segment lining.
(a)
The bending moment
(b)
The axial forceFigure13(a) shows that, on condition that the secondary lining has a constant or specific thickness, with the segment lining thickness increasing constantly, the absolute values of the maximum and minimum bending moment of the double lining both tend to increase, and the absolute values of bending moments for the segment lining at each corresponding section increase as well. From Figure 13(b), it can be seen that the maximum and minimum axial force of the double lining has little change with the increase of the segment lining thickness, but each axial force of the segment lining at a corresponding section obviously increases.From Tables6 and 7, it can be seen that, with the constant increase of the segment lining thickness, the increasing ratio of the bending moment and axial force for the double lining has little relation to the increasing ratio of the double lining thickness. The allocation proportion of the bending moment and axial force between the segment lining and secondary lining has no linear relationship with the ratio of the segment lining and the secondary lining thickness. With the constant increase in the segment lining thickness, the ratio of bending moment and axial force burdened by the segment lining goes up gradually. The allocation proportion of the bending moment of the segment lining at each section of the maximum and the minimum bending moment of the double lining takes a basically consistent ratio, and the allocation proportion of the axial force at each section of the maximum and minimum axial force shows a similar law as well.(2) Different thickness of secondary lining: the thickness of segment lining is 0.50 m, and the thicknesses of the secondary lining are 0.10 m, 0.20 m, 0.30 m, 0.40 m, and 0.50 m, respectively, and the thicknesses of the corresponding double lining are 0.60 m, 0.70 m, 0.80 m, 0.90 m, and 1.00 m, respectively. Through computing, the maximum and minimum bending moment of double lining and the bending moment of the segment lining at the corresponding section can be seen in Figure14(a) and Table 8, and the maximum and minimum axial force of the double lining and the axial force of segment lining at the corresponding section can be seen in Figure 14(b) and Table 9.Table 8
The bending moments extreme values under different thicknesses of secondary lining.
Thickness of secondary lining (cm)
10
20
30
40
50
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum bending moment of the double lining
Bending moment of double lining (kN·m)
146.689
147.398
148.896
150.660
151.945
Bending moment of segment lining at the corresponding section (kN·m)
145.480
138.001
120.650
94.474
72.271
The allocation proportion of bending moment of segment lining
99%
94%
81%
62%
48%
The increase ratio of bending moment of double lining
100%
100%
102%
103%
104%
The section of minimum bending moment of the double lining
Bending moment of double lining (kN·m)
−177.101
−177.886
−179.542
−181.469
−182.814
Bending moment of segment lining at the corresponding section (kN·m)
−175.650
−166.640
−145.771
−116.662
−87.462
The allocation proportion of bending moment of segment lining
99%
94%
81%
64%
48%
The increase ratio of bending moment of double lining
100%
100%
101%
102%
103%Table 9
The axial force extreme values under different thicknesses of secondary lining.
Thickness of secondary lining (cm)
10
20
30
40
50
Thickness of double lining (cm)
60
70
80
90
100
The increase ratio of double lining thickness
100%
117%
133%
150%
167%
The section of maximum axial force of the double lining
Axial force of double lining (MN)
−1.121
−1.121
−1.121
−1.120
−1.120
Axial force of segment lining at the corresponding section (MN)
−1.009
−0.897
−0.818
−0.762
−0.728
The allocation proportion of axial force of segment lining
90%
80%
73%
68%
65%
The increase ratio of axial force of double lining
100%
100%
100%
100%
100%
The section of minimum axial force of the double lining
Axial force of double lining (MN)
−1.058
−1.057
−1.057
−1.056
−1.056
Axial force of segment lining at the corresponding section (MN)
−0.963
−0.856
−0.772
−0.718
−0.686
The allocation proportion of axial force of segment lining
91%
81%
73%
68%
65%
The increase ratio of axial force of double lining
100%
100%
100%
100%
100%Figure 14
Numerical results of double lining under different thicknesses of secondary lining.
(a)
The bending moment
(b)
The axial forceFrom Figure14(a) it can be seen that, on condition that the secondary lining has a specific thickness, with the secondary lining thickness increasing constantly, the absolute values of the maximum and minimum bending moment of the double lining both tend to increase, but the rate of increase is small. And the absolute value of the maximum and minimum bending moment of the segment lining at the corresponding section tends to decrease. From Figure 14(b), it can be seen that the maximum and minimum axial force of the double lining barely changes with the increase in the secondary lining thickness, but each axial force extreme value of the segment lining at the corresponding section tends to decrease.From Tables8 and 9 it can be seen that, with the constant increase of the thickness of the secondary lining, the increasing ratio of bending moment and axial force of the double lining has little relationship with the increasing ratio of the double lining thickness, and the allocation proportion of the bending moment and the axial force between the segment lining and the secondary lining has no linear relationship with the ratio of the segment lining and the secondary lining thickness. With the constant increase in the secondary thickness, the ratio of the bending moment and the axial force that is burdened by the secondary lining decreases quickly. The allocation proportion of the bending moment of the segment lining at each section of the maximum bending moment and the minimum bending moment of the double lining takes a basically consistent ratio, and the same law applies to the axial force as well.


## 4. Conclusions

Double lining of an underwater shield tunnel has wide application in the future. Based on the comparative analysis of the existing numerical models of double lining, a new improved numerical model has been presented in this paper. The lining, joint, and junction surface are simulated by beam, spring, and a combination of compression bar and spring. On the basis of employing a similar model test to contrast the existing and the improved models, the advantages of the improved model are compared to the existing model. Eventually, adopting the improved numerical model, the reasonable element amounts of junction surface needed for implementing the accurate analysis of the double lining of the Huangpu River Tunnel are discussed and the mechanic behaviors are studied when changing the thickness of the segment lining and the secondary lining, respectively, at the same load. The main conclusions are listed as follows:(1)
The existing double lining models have the defects of an abrupt change of bending moment, the appearance of radial tensile stress on the junction surface and contact surface, and so forth, and they are inadequate for the task of setting the waterproof layer and the separation of the junction surface.
(2)
Through the radial compression bar element, the improved numerical model simulates the radial interaction of the junction surface or the contact surface, and the link between the bar element and beam element is hinged. Therefore, it can solve the problems of bending moment abrupt change in the junction surface and contact surface and the appearance of radial tensile stress. Through the coupling of node displacement, it demonstrates that the shear force of the junction surface is directly relevant to the shear displacement. Meanwhile, there are definite calculating methods of the element stiffness parameters.
(3)
Results of the improved numerical model have better consistency with results from the similar model test, and the element amount of the junction surface and contact surface should be discussed in a specific project, which is the precondition of using the improved model to carry out the double lining mechanical analysis.
(4)
With other conditions unchanged, the increase in thickness of the segment lining or the secondary lining will cause the increase of the bending moment of double lining, but this has little influence on the axial force of the double lining. The allocation proportion of the bending moment also has no linear relationship with the ratio of the thickness between the segment lining and the secondary lining.
(5)
With the condition that the thickness of the secondary lining remains unchanged, increasing the thickness of the segment lining, the bending moment and its proportion ratio of segment lining increases more at the corresponding section of the maximum and minimum bending moment of the double lining. It shows that the increased part of the bending moment of the double lining mainly depends on the thickened segment lining. And the axial force of segment lining and its proportion have only a small increase.
(6)
With the condition that the thickness of the segment lining is unchanged, increasing the thickness of the secondary lining, the bending moment and its proportion ratio of segment lining decreases sharply at the corresponding section of the maximum and minimum bending moment of the double lining. It shows that the increased part of the bending moment of the double lining mainly transfers to the thickened secondary lining. The axial force of the segment lining and its proportion have a small decrease.



---

*Source: 430879-2015-09-27.xml*