Sensitivity Analysis of the Waterproof Performance of Elastic Rubber Gasket in Shield Tunnel

Zhiqiang Wang; Zhenyu Lei

Journal & Issue Details

Journal Details

Format: Journal
eISSN: 2444-8656
First Published: 01 Jan 2016
Publication timeframe: 2 times per year
Languages: English

In order to study the waterproof performance of elastic rubber gasket in shield tunnel lining joints, an innovative sensitivity analysis method is proposed by combining the Monte Carlo method with the stochastic finite element method (FEM) in this paper. The sensitivity values of the waterproof performance respecting to elastic rubber gaskets are obtained via the ANSYS Probabilistic Design System (PDS) module, in which the parameters of material hardness, coordinates of the hole center, apertures are selected as random input variables. Meantime, the extent of the tolerance effect of the random parameters on the waterproof performance is explored.

Keywords

elastic rubber gasket
sensitivity analysis
waterproof performance
Monte Carlo method
stochastic FEM

Article

Introduction

Waterproofing is one of the cruxes for underground projections, which should be paid especial attention to in shield tunnels. Segment joints are the main leakage points and rubber gaskets are generally used for joints waterproofing. Therefore, waterproof properties of the rubber gasket should be strictly controlled.

There are many researches focusing on the impermeability performance of waterproof materials. Shalabi et al. [1,2] investigated the waterproof performance of gasket-in-groove under different contact stress states and studied the differences of waterproof ability between circular and longitudinal seams under the action of cyclic loading. Based on the variation law of mechanical properties of materials, Shi et al. [3] proposed a time-dependent constitutive model of rubber gaskets, and the relevant parameters were obtained by the material ageing test. Gong et al. [4] studied the water resistance characteristics of rubber gaskets for river crossing tunnels by using the new high-water pressure and automatic loading device.

In addition to the impermeability of waterproof materials, the optimization design method of rubber gaskets has also attracted considerable attention. Xiang and Shi [5] analyzed the deforming characteristics of elastic gaskets under different levels of compression, the distribution of contact stress and the compressive force needed during installation through Finite Element Method (FEM). Then the section that satisfied all design requirements was found based on the FEM analysis results by repeated adjustments and trials. Novak et al. [6] developed a probabilistic technique to evaluate gasket designs under different assembly conditions. Jones [7] outlined the procedures about optimizing gasket geometries and presented a comparison which showed the benefits of using such techniques. Krishnan [8] proposed a Computer Aided Engineering (CAE) approach to design, evaluate and optimize automotive powertrain gaskets, seals and gasket/seal assemblies.

Due to the elasticity of rubber materials, the gasket deforms easily when it is squeezed. Some sections of elastic gaskets, provided by different manufacturers, were measured by researchers which were used in the shield tunnel of a city subway, as shown in Figure 1. It can be seen from Figure 1, all of the geometry sizes, locations and shapes of the internal holes are wrong. Based on the optimization design experience, those incorrect parts have significant influences on the waterproof performance of rubber gaskets.

Comparison of the product and design drawings of two manufacturers (Note: dotted line refers to the design drawing): (a) No.1 product of manufacturers, (b) No.2 product of manufacturers.

In addition, mechanical parameters, such as hardness of synthetic rubber materials, are largely influenced by places of production, compositions and temperatures. Generally, it is required that the hardness should change in a range of ±5 degrees by considering the randomness of mechanical properties.

Although many studies have been conducted to analysis the waterproof performance of EPDM (Ethylene Propylene Diene Monomer) rubber materials and the gasket optimization design method, the effect on the waterproof performance of rubber gaskets due to the continuous change of the structural parameter still remains unknown. In addition, the effects of the uncertainties of geometry sizes, machining errors and material parameters on the waterproof performance of rubber gaskets are not considered in previous optimization processes.

Based on the aforementioned, it is known that traditional design methods for the elastic rubber gasket have many defects. Therefore, aiming at dealing with the problems involved in the traditional methods, an innovative sensitivity analysis method is proposed, where the material hardness, coordinates of the hole center, apertures are regarded as random input variables instead of deterministic ones. Monte Carlo method and FEM are combined to acquire the sensitivity values of the waterproof performance with regard to the uncertain structural parameters.

The rest of the paper is organized as follows. After the introduction, we briefly explain the background theory of stochastic FEM based on the Monte Carlo method and ANSYS PDS module in section 2. Section 3 discusses a new method for the sensitivity analysis of elastic rubber gaskets in shield tunnel lining joints. Finally, we close the paper with some conclusions and remarks.

Stochastic FEM based on the Monte Carlo method and ANSYS PDS module

Since the geometry sizes, machining errors, materials, actual working loads are all random for the same batch of products in practical engineering, we need to use stochastic FEM to carry out the detailed analysis. The stochastic FEM [9] used in this paper is a combination of the Monte Carlo method and the FEM.

2.1

The principle of Monte Carlo stochastic method

Monte Carlo method [10], also known as statistical test method or random simulation method, is based on the theory of statistical sampling. The structural stochastic analysis process carried out by using Monte Carlo method can be summarized as following three basic steps.

(1)

Sampling random variables. Random sampling based on the known probability distribution of the basic random variables.

(2)

Solving structural response. For each sample, the finite element method is used to calculate the response of the structure.

(3)

Statistical analysis of the response. According to the calculation results of all samples, the mean value, variance and even probability distribution function of the response are obtained.

2.2

Sampling methods for random variables

A major drawback of the Monte Carlo method is the need of large sampling space. In order to overcome this shortcoming, this paper uses the Latin Hypercube Sampling (LHS) method [11] which reduces the sampling number by improving the convergence rate. LHS [12] method does not change the mean and variance of sampling space, and can greatly improve the computing efficiency. Therefore, it is more suitable for solving the stochastic response of engineering structures.

2.3

Probabilistic Design System (PDS) of ANSYS

The PDS module of ANSYS combines the finite element method with the probability design technology. Based on the probability design of FEM, it can be used to study the uncertain effect on waterproof properties of rubber gaskets. The parametric modeling in ANSYS-PDS is often employed to carry out probability design analysis and the APDL (ANSYS Parametric Design Language) command flow is generally used to build the parametric model for probability design analysis. Based on the parametric model, the new model can be built by just changing some parameters, so it is more efficient and less costly to repeatedly analyze the effects of parameters on the waterproof performance, such as hardness, size, etc. The probability design analysis process based on ANSYS-PDS is shown in Figure 2.

Sensitivity analysis method

In order to tackle the problems produced in the traditional design methods, a new sensitivity analysis method of the waterproof performance with regard to the uncertain structural parameters has been presented in this section.

3.1

Constructing parameterized model

The geometry sizes of the rubber gasket and segment groove studied in this paper are shown in Figure 3. The selected random parameters are as follows.

(1)

Material parameters of rubber C₁₀, C₀₁. The constitutive model used in this paper is the Mooney-Rivlin two-parameter model [13]. The model can describe the deformation characteristics of rubber materials in 150%.

(2)

The center coordinates of all the holes in the cross section x_i, y_i, i = 1 ∼ 4. The effect of the position of the hole on the waterproof property of the sealing gasket is analyzed. The section of the sealing gasket is symmetrical, and there are eight holes totally in the cross section and the number is 1 ∼ 4.

(3)

Radius of each round hole R_i, i = 1 ∼ 4.

3.2

Statistical characteristics of random parameters

The elastic rubber gasket in shield tunnel lining joints is EPDM synthetic rubber. In general, the hardness of rubber materials is clearly defined and the hardness of IRHD (International Rubber Hardness Standard) is often adopted, and the hardness tolerance required is ±5 degrees [14]. The relationship between hardness H_r, elastic modulus E and mechanical characteristic constants C₀₁, C₁₀ are: (1) $lgE = 0.0198 H_{r} - 0.5432$ lgE = 0.0198{H_r} - 0.5432(2) $E = 6 C_{10} (1 + \frac{C_{01}}{C_{10}})$ E = 6{C_{10}}(1 + {{{C_{01}}} \over {{C_{10}}}}) where $\frac{C_{01}}{C_{10}} = 0.05$ {{{C_{01}}} \over {{C_{10}}}} = 0.05 , C₁₀=0.700MPa, C₀₁=0.035MPa, E=4.41MPa, Poisson's ratio μ=0.499. The maximum and minimum limit values of C₀₁ and C₁₀ of the equations (1) and (2) can be obtained as follows: (3) $C_{10} = {\begin{array}{l} 0.88 MPa \\ 0.5578 MPa \end{array} C_{01} = {\begin{array}{l} 0.044 MPa \\ 0.0279 MPa \end{array}$ {C_{10}} = \left\{{\matrix{{0.88\rm{MPa}} \hfill \cr {0.5578\rm{MPa}} \hfill \cr}} \right.{C_{01}} = \left\{{\matrix{{0.044\rm{MPa}} \hfill \cr {0.0279\rm{MPa}} \hfill \cr}} \right.

The center coordinates of the holes are x₁=x₃=2.75mm, x₂=3.5mm, x₄=9.5mm, y₁=3.5mm, y₂=8.3mm, y₃=13mm, y₄=11.5mm, respectively. The apertures are R₁=R₃=1.75mm, R₂=R₄=2mm. The C₀₁, C₁₀, center coordinates x_i and y_i of all holes and radius of each hole R_i are selected as random input variables. The correlation between the variables is ignored. The random distributions and corresponding characteristic values of each random variable are as shown in Table 1. The main characteristics of each variable are determined according to the reference [15].

Table 1

Statistical characteristics of random input variables

Variable name	Distribution pattern	Mean value	Standard deviation	Min limit value	Max limit value
C₀₁/MPa	truncated normal	0.035	0.05C₀₁	0.0279	0.044
C₁₀/MPa	truncated normal	0.700	0.05C₁₀	0.5578	0.88
x_i,i = 1,2,3,4/mm	truncated normal	x_i	0.15	x_i−0.15	x_i+0.15
y_i,i = 1,2,3,4/mm	truncated normal	y_i	0.15	y_i−0.15	y_i+0.15
R_i,i = 1,2,3,4/mm	truncated normal	R_i	0.03	R_i−0.05	R_i+0.01

3.3

Finite element model

In this paper, the unilateral compression model is employed to simulate the compression process of the rubber gasket, as in Figure 4. The rubber gasket can be considered as the flexible body and the concrete groove can be regarded as the rigid body. There are rigid-flexible contacts between the groove and the rubber and flexible-flexible contacts between the rubber and the rubber.

Unilateral compress FE model of the rubber gasket.

Due to the hyper elasticity characteristic of rubber materials, the model is established with the PLANE182 element. The contact elements use TARGE169 and CONTA172 elements with friction coefficient of 0.25 for rubber materials.

According to the maximum allowable opening amount in engineering, the unilateral compression can be set as 5 mm. All the degrees of the nodes of the groove are constrained and the X-direction degrees as well as the rotational degrees of the nodes of the upper pressure plate are constrained. All nodes of the upper pressure plate are specified to have a Y-direction displacement of 5 mm (compression direction).

3.4

Sensitivity calculation and analysis

3.4.1

Sensitivity analysis based on Spearman rank correlation coefficient

Sensitivity analysis (SA) [16,17,18,19,20] is an important step in the process of model establishment and result discussion. In this section, the sensitivity analysis method based on the Spearman rank correlation coefficient is adopted [21].

The sequence of random structure parameters of the sealing gasket is x={C₁₀, C₀₁, x_i, y_i, R_i}, (i = 1 ∼ 4)and the 14 random variables are all subject to the truncated normal distribution. Assuming that the n times random simulation is carried out, the structure response values of n times are obtained. Since the sealing property of non-expansion elastic gaskets is mainly subjected to the elastic pressure produced by the surface contact stress [5], the average contact stress $\bar{σ_{u}}$ \overline {{\sigma _u}} on the upper surface and the average contact stress $\bar{σ_{d}}$ \overline {{\sigma _d}} on the bottom surface are set as random output variables, recorded as $\bar{σ_{i 1}}, \bar{σ_{i 2}}, \dots, \bar{σ_{in}}$ \overline {{\sigma _{i1}}} ,\overline {{\sigma _{i2}}} , \ldots ,\overline {{\sigma _{in}}} , i = u,d, respectively.

Let ${(x_{jn}, \bar{σ_{in}})}_{j = 1}^{14}$ \{({x_{jn}},\overline {{\sigma _{in}}})\} _{j = 1}^{14} denote 14 independent and identically distributed data pairs which are composed of the sample values of the n random simulation operations x_j₁, x_j₂, . . . , x_jn and the structure response values of the n random simulation operations $\bar{σ_{i 1}}, \bar{σ_{i 2}}, \dots, \bar{σ_{in}}$ \overline {{\sigma _{i1}}} ,\overline {{\sigma _{i2}}} , \ldots ,\overline {{\sigma _{in}}} .

For each data pair, the Spearman rank correlation coefficient r_sj [22] is calculated by the equations (4) and (5). (4) $r_{sj} = - \frac{6 \sum_{k = 1}^{n} {(P_{kj} - Q_{kj})}^{2}}{n (n^{2} - 1)} (j = 1 \sim 14)$ {r_{sj}} = - {{6\sum\nolimits_{k = 1}^n {{({P_{kj}} - {Q_{kj}})}^2}} \over {n({n^2} - 1)}}(j = 1 \sim 14)(5) $- 1 \leq r_{sj} \leq 1$ - 1 \le {{\rm{r}}_{sj}} \le 1

Let P_{k j} be the rank of x_j among x₁, x₂, . . . ,x₁₄; Q_{k j} be the rank of the corresponding structural response $\bar{σ_{ij}}$ \overline {{\sigma _{ij}}} ; is the total sampling times.

The Spearman rank correlation coefficient is an important parameter used to test the correlation between variables in nonparametric statistics. If the coefficient is close to 1 or −1, it indicates that the input variables significantly affect the output variables. If the coefficient is close to 0, it illustrates that the effect is weak. Besides, if the coefficient is positive, the output variables increase with the increase of the input variables. If the coefficient is negative, the output variables decrease with the increase of the input variables [23].

3.4.2

Sensitivity calculation results

In the ANSYS-PDS module, the Monte Carlo method and Latin Hypercube Sampling method are chosen as the probability analysis methods and the cycle times is set to be 500 times for performing probability design analysis. Then, the trends of $\bar{σ_{u}}$ \overline {{\sigma _u}} and $\bar{σ_{d}}$ \overline {{\sigma _d}} can be obtained and are shown in Figure 5. It can be seen that with the increase of sampling frequency, the width of the confidence interval becomes smaller. Meantime, the trends of $\bar{σ_{u}}$ \overline {{\sigma _u}} and $\bar{σ_{d}}$ \overline {{\sigma _d}} gradually tend to horizontal, which indicates that 500 sampling times is enough.

Trends of average contact stresses:(a) Average contact stress on the bottom surface, (b) Average contact stress on the upper surface.

The specific sensitivity values of the random input variables based on the Spearman rank correlation coefficient are listed in Table 2.

Table 2

The sensitivity values of the input variables

Var.	C₀₁	C₁₀	x₁	x₂	x₃	x₄	y₁	y₂	y₃	y₄	R₁	R₂	R₃	R₄
$\bar{σ_{d}}$ \overline {{\sigma _d}}	0.049	0.971	0.126	−0.102	0.094	0.011	0.010	−0.019	0.010	0.032	−0.072	−0.040	−0.138	−0.102
$\bar{σ_{u}}$ \overline {{\sigma _u}}	0.036	0.918	0.080	−0.107	0.074	−0.277	−0.027	0.031	−0.002	0.035	−0.069	0.031	−0.062	−0.065

As can be seen from Table 2, the C₁₀ is the most effective factor for the average contact stress of the bottom surface, followed by R₃, x₁, R₄, x₂, and the influence of C₁₀ is greater than that of the other four variables. The C₁₀ is also the most significant factor for the average contact stress of the upper surface, followed by x₄, x₂, and the influence of C₁₀ is greater than that of the other two variables. Furthermore, in the random input variable space, it will make the average contact stresses of both the upper and bottom surfaces increase by increasing the rubber hardness C₁₀. Besides, it will make the average contact stress of the bottom surface increase by increasing x₁. The average contact stresses on both the upper and bottom surfaces will decrease with the increase of x₂. The average contact stress on the upper surface will decrease with the increase of x₄. The average contact stress on the bottom surface will decrease with the increase of R₃ or R₄.

3.4.3

Sensitivity influence of tolerance

The maximum and minimum limit values of random input variables are calculated by the mean values and tolerances. In this section, the tolerances of random input variables are changed in order to analyze their influence on the sensitivity.

(1) Rubber hardness tolerance

The rubber hardness tolerances are chosen as ±5 degrees, ±4 degrees, ±3 degrees, ±2 degrees, ±1 degrees, respectively. The maximum and minimum limit values of C₀₁ and C₁₀ can be calculated using the equations (2), (4) and (5). The statistical characteristics of the other input variables are the same as those in Table 1. The same analysis method and sampling times as the above are used to get the sensitivity values of average contact stresses on upper and bottom surfaces, as shown in Figure 6. From Figure 6, it can be found that the rubber hardness tolerance has little effect on the contact stress sensitivity when the rubber hardness tolerance is in the ranges of ±5∼±3 degrees. When the rubber hardness tolerance is restricted within the range of ±2 degrees by taking stricter production standards, the effect of the rubber hardness on the contact stress decreases significantly and the effects of other variables on the contact stress increase significantly, such as the hole center location and aperture.

The sensitivity influence of rubber hardness tolerance (Only list the top five variables with great influence):(a) Average contact stress on the bottom surface, (b) Average contact stress on the upper surface.

(2) Coordinate deviation

The abscissa deviations are selected as ±0.2 mm, ±0.17 mm, ±0.15 mm and ±0.13 mm, respectively. The maximum and minimum limits of abscissas of hole centers can be obtained by adding these deviations with the mean values. The statistical characteristics of the other input variables are the same as those in Table 1. The same analysis method and sampling times as the above are used to get the sensitivity values of average contact stresses on the upper and bottom surfaces, as shown in Figure 7. From Figure 7, it can be obtained that the influence of the change of abscissas on the average contact stress on the bottom surface decreases if the abscissa deviation is strictly controlled. Because the variation range of sensitivity values is small, it may be due to the influence of rubber hardness on the contact stress is more significant than other variables when the tolerance of rubber hardness is in the range of ±5 degrees in this condition. If the abscissa position deviation is strictly controlled, the influence of the abscissa change of the 4th hole center on the average contact stress on the upper surface decreases and the change value of the sensitivity is 0.113. The change of other abscissas has little effect on the average surface contact stress.

The sensitivity influence of horizontal position deviation:(a) Average contact stress on the bottom surface, (b) Average contact stress on the upper surface.

From Table 2, it can be found that the ordinates of hole centers have little effect on the contact stress. Therefore, the influence of the change of the vertical position deviation on the contact stress sensitivity is not discussed in this paper.

(3) Aperture deviation

The maximum limit of aperture is kept as R_i+0.01, and the minimum limits of aperture are separately set as R_i−0.1, R_i−0.07, R_i−0.05 and R_i−0.03. The statistical characteristics of the other input variables are the same as those in Table 1. The same analysis method and sampling times as the above are used to get the sensitivity values of average contact stresses on the upper and bottom surfaces, as shown in Figure 8. From Figure 8, it can be seen that the stricter the control of the aperture deviation, the less influence on the average contact stress on the bottom surface. The largest variation value of the sensitivity is 0.048. The stricter the control of the R₁ and R₄ deviations, the less influence on the average contact stress on the upper surface. The stricter the control of the R₂ and R₃ deviations, the more influence on the average contact stress on the upper surface. But the change extents are small because the rubber hardness has a greater effect on the contact stress than the aperture.

The sensitivity influence of aperture deviation:(a) Average contact stress on the bottom surface, (b) Average contact stress on the upper surface.

Conclusions and remarks

(1)

This paper focuses on the sensitivity analysis of elastic rubber gaskets in shield tunnel lining segments. In order to obtain the correlation between the parameters and the average contact stress of sealing gaskets, the ANSYS-PDS module is employed by combining the Monte Carlo method with the FEM, and the parameters of material hardness, hole center coordinates and apertures are selected as the random input variables. The analysis method and the results of this paper have a certain reference value on optimal design of elastic rubber gaskets in the future.

(2)

By controlling the tolerance of the parameters in turn, the influence of various parameters on the average contact stress of the sealing gasket is obtained. The mechanical parameters of rubber materials are the main factors that influence the water pressure resistance, especially C₁₀. The influence of rubber hardness on the contact stress is the most significant when the tolerance of rubber hardness is in the range of ±5∼±3 degrees. When the rubber hardness tolerance is restricted in the range of ±2 degrees under strict production standards, the influence of the rubber hardness on the contact stress decreases significantly, while the influence of other parameters on the contact stress increases remarkably.

(3)

As for geometric parameters, the greatest impact on the contact stress is the position of bottom hole 1. The effect of the size of the two round holes in the middle part of the section on the surface contact stress is larger than that of the other round holes.

Figures & Tables

Statistical characteristics of random input variables

Variable name	Distribution pattern	Mean value	Standard deviation	Min limit value	Max limit value
C₀₁/MPa	truncated normal	0.035	0.05C₀₁	0.0279	0.044
C₁₀/MPa	truncated normal	0.700	0.05C₁₀	0.5578	0.88
x_i,i = 1,2,3,4/mm	truncated normal	x_i	0.15	x_i−0.15	x_i+0.15
y_i,i = 1,2,3,4/mm	truncated normal	y_i	0.15	y_i−0.15	y_i+0.15
R_i,i = 1,2,3,4/mm	truncated normal	R_i	0.03	R_i−0.05	R_i+0.01

The sensitivity values of the input variables

Var.	C₀₁	C₁₀	x₁	x₂	x₃	x₄	y₁	y₂	y₃	y₄	R₁	R₂	R₃	R₄
σd¯\overline {{\sigma _d}}	0.049	0.971	0.126	−0.102	0.094	0.011	0.010	−0.019	0.010	0.032	−0.072	−0.040	−0.138	−0.102
σu¯\overline {{\sigma _u}}	0.036	0.918	0.080	−0.107	0.074	−0.277	−0.027	0.031	−0.002	0.035	−0.069	0.031	−0.062	−0.065

References

[1]Shalabi, F.I., Cording, E.J., Paul, S.L., 2016. Sealant behavior of gasketed segmental tunnel lining-Conceptual model. Geomechanics and Tunnelling. 9(4), 345–355. https://doi.org/10.1002/geot.201500030ShalabiF.I.CordingE.J.PaulS.L.2016Sealant behavior of gasketed segmental tunnel lining-Conceptual modelGeomechanics and Tunnelling94345355https://doi.org/10.1002/geot.20150003010.1002/geot.201500030Search in Google Scholar

[2]Shalabi, F.I., Cording, E.J., Paul, S.L., 2012. Concrete segment tunnel lining sealant performance under earthquake loading. Tunnelling and Underground Space Technology. 31(9), 51–60. https://doi.org/10.1016/j.tust.2012.04.006ShalabiF.I.CordingE.J.PaulS.L.2012Concrete segment tunnel lining sealant performance under earthquake loadingTunnelling and Underground Space Technology3195160https://doi.org/10.1016/j.tust.2012.04.00610.1016/j.tust.2012.04.006Search in Google Scholar

[3]Shi, C.H., Cao, C.Y., Lei, M.F., Peng, L.M., Shen, J.J., 2015. Time-dependent performance and constitutive model of EPDM rubber gasket used for tunnel segment joints. Tunneling and Underground Space Technology. 50(8), 490–498. https://doi.org/10.1016/j.tust.2015.09.004ShiC.H.CaoC.Y.LeiM.F.PengL.M.ShenJ.J.2015Time-dependent performance and constitutive model of EPDM rubber gasket used for tunnel segment jointsTunneling and Underground Space Technology508490498https://doi.org/10.1016/j.tust.2015.09.00410.1016/j.tust.2015.09.004Search in Google Scholar

[4]Gong, C.J., Ding, W.Q., Jin, Y.L., Guo, X.H., Tuo, Y.F., 2015. Waterproofing performance of shield-driven tunnel's segment joint under ultrahigh water pressure. Tunneling and Underground Construction. 42(1), 410–418. https://doi.org/10.1061/9780784413449.040GongC.J.DingW.Q.JinY.L.GuoX.H.TuoY.F.2015Waterproofing performance of shield-driven tunnel's segment joint under ultrahigh water pressureTunneling and Underground Construction421410418https://doi.org/10.1061/9780784413449.04010.1061/9780784413449.040Search in Google Scholar

[5]Xiang, K., Shi, X.W., 2008. Design and optimization of elastic gasket section of shield tunnel lining. Chinese Journal of Underground Space and Engineering. 4(2), 361–364. https://doi.org/CNKI:SUN:BASE.0.2008-02-031XiangK.ShiX.W.2008Design and optimization of elastic gasket section of shield tunnel liningChinese Journal of Underground Space and Engineering42361364https://doi.org/CNKI:SUN:BASE.0.2008-02-031Search in Google Scholar

[6]Novak, G., Sadowski, M., Hu, Z.Q., 1995. A probabilistic gasket design method. SAE Transactions. 1467–1474. https://doi.org/10.4271/950765NovakG.SadowskiM.HuZ.Q.1995A probabilistic gasket design methodSAE Transactions14671474https://doi.org/10.4271/95076510.4271/950765Search in Google Scholar

[7]Jones, P.A., 1995. Gasket design the way forward. SAE Technical Paper.JonesP.A.1995Gasket design the way forwardSAE Technical Paper10.4271/950319Search in Google Scholar

[8]Krishnan, M.R., Sural, V., Doshi, V., 1994. Design optimization of automotive powertrain gaskets and seals using cae techniques. SAE Transactions. 245–251. https://doi.org/10.4271/940288KrishnanM.R.SuralV.DoshiV.1994Design optimization of automotive powertrain gaskets and seals using cae techniquesSAE Transactions245251https://doi.org/10.4271/94028810.4271/940288Search in Google Scholar

[9]Peng, X.Q., Liu, G., Wu, L.Y., Liu, G.R., Lam, K.Y., 1998. A stochastic finite element method for fatigue reliability analysis of gear teeth subjected to bending. Computational Mechanics. 21(3), 253–261. https://doi.org/10.1007/s004660050300PengX.Q.LiuG.WuL.Y.LiuG.R.LamK.Y.1998A stochastic finite element method for fatigue reliability analysis of gear teeth subjected to bendingComputational Mechanics213253261https://doi.org/10.1007/s00466005030010.1007/s004660050300Search in Google Scholar

[10]Shinozuka, M., Wen, Y.K., 1972. Monte Carlo solution of nonlinear vibrations. AIAA Journal. 10(1), 37–40. https://doi.org/10.2514/3.50064ShinozukaM.WenY.K.1972Monte Carlo solution of nonlinear vibrationsAIAA Journal1013740https://doi.org/10.2514/3.5006410.2514/6.1971-213Search in Google Scholar

[11]Mckay, M.D., Conover, R.J.B.J., 1979. A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics. 21(2), 239–245. https://doi.org/10.2307/1271432MckayM.D.ConoverR.J.B.J.1979A comparison of three methods for selecting values of input variables in the analysis of output from a computer codeTechnometrics212239245https://doi.org/10.2307/127143210.2307/1271432Search in Google Scholar

[12]Shields, M.D., Zhang, J., 2016. The generalization of Latin hypercube sampling. Reliability Engineering and System Safety. 148(8), 96–108. https://doi.org/10.1016/j.ress.2015.12.002ShieldsM.D.ZhangJ.2016The generalization of Latin hypercube samplingReliability Engineering and System Safety148896108https://doi.org/10.1016/j.ress.2015.12.00210.1016/j.ress.2015.12.002Search in Google Scholar

[13]Pereira, C.E.L., Bittencourt, M.L., 2010. Topological sensitivity analysis for a two-parameter Mooney-Rivlin hyperelastic constitutive model. Latin American Journal of Solids and Structures. 7(4), 391–411. https://doi.org/10.1590/s1679-78252010000400002PereiraC.E.L.BittencourtM.L.2010Topological sensitivity analysis for a two-parameter Mooney-Rivlin hyperelastic constitutive modelLatin American Journal of Solids and Structures74391411https://doi.org/10.1590/s1679-7825201000040000210.1590/S1679-78252010000400002Search in Google Scholar

[14]GB 18173.4-2010., 2010. Polymer water-proof materials part 4: Rubber gasket for shield-driven tunnel. Standards Press of China, Beijing.GB 18173.4-20102010Polymer water-proof materials part 4: Rubber gasket for shield-driven tunnelStandards Press of ChinaBeijingSearch in Google Scholar

[15]ISO 3302-1:2014(E)., 2014. Rubber-Tolerances for products-Part 1: Dimensional tolerances.ISO 3302-1:2014(E)2014Rubber-Tolerances for products-Part 1: Dimensional tolerancesSearch in Google Scholar

[16]Emanuele, B., Elmar, P., 2016. Sensitivity analysis: A review of recent advances. European Journal of Operational Research 248(3), 869–887. https://doi.org/10.1016/j.ejor.2015.06.032EmanueleB.ElmarP.2016Sensitivity analysis: A review of recent advancesEuropean Journal of Operational Research2483869887https://doi.org/10.1016/j.ejor.2015.06.03210.1016/j.ejor.2015.06.032Search in Google Scholar

[17]Gan, Y.J., Duan, Q.Y., Gong, W., Tong, C., Sun, Y.W., Chu, W., Ye, A.Z., Miao, C.Y., Di, Z.H., 2014. A comprehensive evaluation of various sensitivity analysis methods: A case study with a hydrological model. Environmental Modelling and Software 51, 269–285. https://doi.org/10.1016/j.envsoft.2013.09.031GanY.J.DuanQ.Y.GongW.TongC.SunY.W.ChuW.YeA.Z.MiaoC.Y.DiZ.H.2014A comprehensive evaluation of various sensitivity analysis methods: A case study with a hydrological modelEnvironmental Modelling and Software51269285https://doi.org/10.1016/j.envsoft.2013.09.03110.1016/j.envsoft.2013.09.031Search in Google Scholar

[18]Hornberger, G.M., Spear, R.C., 1981. An approach to the preliminary analysis of environmental systems. Journal of Environmental Management. 12 (1), 7–18.HornbergerG.M.SpearR.C.1981An approach to the preliminary analysis of environmental systemsJournal of Environmental Management121718Search in Google Scholar

[19]Abdelal, G.F., Cooper, J.E., Robotham, A.J., 2013. Reliability assessment of 3D space frame structures applying stochastic finite element analysis. International Journal of Mechanics and Materials in Design. 9(1): 1–9. https://doi.org/10.1007/s10999-011-9168-0AbdelalG.F.CooperJ.E.RobothamA.J.2013Reliability assessment of 3D space frame structures applying stochastic finite element analysisInternational Journal of Mechanics and Materials in Design9119https://doi.org/10.1007/s10999-011-9168-010.1007/s10999-011-9168-0Search in Google Scholar

[20]Jung, D., Gea, H.C., 2004. Compliant mechanism design with non-linear materials using topology optimization. International Journal of Mechanics and Materials in Design. 1(2), 157–171. https://doi.org/10.1007/s10999-004-1494-zJungD.GeaH.C.2004Compliant mechanism design with non-linear materials using topology optimizationInternational Journal of Mechanics and Materials in Design12157171https://doi.org/10.1007/s10999-004-1494-z10.1115/DETC2002/DAC-34150Search in Google Scholar

[21]Yue, S., Pilon, P., Cavadias, G., 2002. Power of the Mann-Kendall and Spearman's rho tests for detecting monotonic trends in hydrological series. Journal of Hydrology. 259(1–4), 254–271. https://doi.org/10.1016/S0022-1694(01)00594-7YueS.PilonP.CavadiasG.2002Power of the Mann-Kendall and Spearman's rho tests for detecting monotonic trends in hydrological seriesJournal of Hydrology2591–4254271https://doi.org/10.1016/S0022-1694(01)00594-710.1016/S0022-1694(01)00594-7Search in Google Scholar

[22]Xu, W., Hou, Y., Hung, Y.S., Zou, Y.X., 2013. A comparative analysis of Spearman's rho and Kendall's tau in normal and contaminated normal models. Signal Processing. 93(1), 261–276. https://doi.org/10.1016/j.sigpro.2012.08.005XuW.HouY.HungY.S.ZouY.X.2013A comparative analysis of Spearman's rho and Kendall's tau in normal and contaminated normal modelsSignal Processing931261276https://doi.org/10.1016/j.sigpro.2012.08.00510.1016/j.sigpro.2012.08.005Search in Google Scholar

[23]Lei, Z.Y., Wang, Z.Q., 2020. Analysis of stress relaxation characteristics of rubber sealing gaskets under the influence of random parameters. Iranian Journal of Science and Technology-Transactions of Mechanical Engineering. 2, 1–8. https://doi.org/10.1007/S40997-020-00353-wLeiZ.Y.WangZ.Q.2020Analysis of stress relaxation characteristics of rubber sealing gaskets under the influence of random parametersIranian Journal of Science and Technology-Transactions of Mechanical Engineering218https://doi.org/10.1007/S40997-020-00353-w10.1007/s40997-020-00353-wSearch in Google Scholar

Applied Mathematics and Nonlinear Sciences

Sensitivity Analysis of the Waterproof Performance of Elastic Rubber Gasket in Shield Tunnel

Zhiqiang Wang andZhiqiang WangZhenyu LeiZhenyu Lei

Published Online: 22 Mar 2021

Page range: -

Received: 01 Dec 2020

Accepted: 31 Jan 2021

Keywords

Fig. 1

Fig. 2

Fig. 3

Fig. 4

Fig. 5

Fig. 6

Fig. 7

Fig. 8

Fig. 1

Fig. 2

Fig. 3

Fig. 4

Fig. 5

Fig. 6

Fig. 7

Fig. 8

Statistical characteristics of random input variables

The sensitivity values of the input variables

Zhiqiang Wang and
Zhiqiang Wang
Zhenyu Lei
Zhenyu Lei