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

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. [https://doi.org/10.1061/9780784413449.040GongC.J.DingW.Q.JinY.L.GuoX.H.TuoY.F.

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.

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.

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.

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.

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

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

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.

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.

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.

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.

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.

Material parameters of rubber _{10}, _{01}. 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%.

The center coordinates of all the holes in the cross section _{i}_{i}

Radius of each round hole _{i}

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 _{r}_{01}, _{10} are:
_{10}=0.700MPa, _{01}=0.035MPa, _{01} and _{10} of the equations (1) and (2) can be obtained as follows:

The center coordinates of the holes are _{1}=_{3}=2.75mm, _{2}=3.5mm, _{4}=9.5mm, _{1}=3.5mm, _{2}=8.3mm, _{3}=13mm, _{4}=11.5mm, respectively. The apertures are _{1}=_{3}=1.75mm, _{2}=_{4}=2mm. The _{01}, _{10}, center coordinates _{i}_{i}_{i}

Statistical characteristics of random input variables

Variable name | Distribution pattern | Mean value | Standard deviation | Min limit value | Max limit value |
---|---|---|---|---|---|

_{01}/MPa | truncated normal | 0.035 | 0.05_{01} | 0.0279 | 0.044 |

_{10}/MPa | truncated normal | 0.700 | 0.05_{10} | 0.5578 | 0.88 |

_{i} | truncated normal | _{i} | 0.15 | _{i} | _{i} |

_{i} | truncated normal | _{i} | 0.15 | _{i} | _{i} |

_{i} | truncated normal | _{i} | 0.03 | _{i} | _{i} |

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.

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).

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 [https://doi.org/10.1016/S0022-1694(01)00594-7YueS.PilonP.CavadiasG.

The sequence of random structure parameters of the sealing gasket is _{10}, _{01}, _{i}_{i}_{i}

Let
_{j}_{1}, _{j}_{2}, . . . , _{jn}

For each data pair, the Spearman rank correlation coefficient _{sj}

Let _{k j}_{j}_{1}, _{2}, . . . ,_{14}; _{k j}

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].

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

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

The sensitivity values of the input variables

Var. | _{01} | _{10} | _{1} | _{2} | _{3} | _{4} | _{1} | _{2} | _{3} | _{4} | _{1} | _{2} | _{3} | _{4} |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|

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 | |

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 _{10} is the most effective factor for the average contact stress of the bottom surface, followed by _{3}, _{1}, _{4}, _{2}, and the influence of _{10} is greater than that of the other four variables. The _{10} is also the most significant factor for the average contact stress of the upper surface, followed by _{4}, _{2}, and the influence of _{10} 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 _{10}. Besides, it will make the average contact stress of the bottom surface increase by increasing _{1}. The average contact stresses on both the upper and bottom surfaces will decrease with the increase of _{2}. The average contact stress on the upper surface will decrease with the increase of _{4}. The average contact stress on the bottom surface will decrease with the increase of _{3} or _{4}.

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 _{01} and _{10} 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.

(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.

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 _{i}_{i}_{i}_{i}_{i}_{1} and _{4} deviations, the less influence on the average contact stress on the upper surface. The stricter the control of the _{2} and _{3} 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.

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.

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 _{10}. 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.

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.

#### Statistical characteristics of random input variables

Variable name | Distribution pattern | Mean value | Standard deviation | Min limit value | Max limit value |
---|---|---|---|---|---|

_{01}/MPa | truncated normal | 0.035 | 0.05_{01} | 0.0279 | 0.044 |

_{10}/MPa | truncated normal | 0.700 | 0.05_{10} | 0.5578 | 0.88 |

_{i} | truncated normal | _{i} | 0.15 | _{i} | _{i} |

_{i} | truncated normal | _{i} | 0.15 | _{i} | _{i} |

_{i} | truncated normal | _{i} | 0.03 | _{i} | _{i} |

#### The sensitivity values of the input variables

Var. | _{01} | _{10} | _{1} | _{2} | _{3} | _{4} | _{1} | _{2} | _{3} | _{4} | _{1} | _{2} | _{3} | _{4} |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|

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 | |

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 |

Regarding new wave distributions of the non-linear integro-partial Ito differential and fifth-order integrable equations Nonlinear Mathematical Modelling of Bone Damage and Remodelling Behaviour in Human Femur Value Creation of Real Estate Company Spin-off Property Service Company Listing Entrepreneur's Passion and Entrepreneurial Opportunity Identification: A Moderated Mediation Effect Model Applications of the extended rational sine-cosine and sinh-cosh techniques to some nonlinear complex models arising in mathematical physics Study on the Classification of Forestry Infrastructure from the Perspective of Supply Based on the Classical Quartering Method A Modified Iterative Method for Solving Nonlinear Functional Equation A comprehensive evaluation of county economies in the Beijing-Tianjin-Hebei Region based on entropy TOPSIS analysis New Principles of Non-Linear Integral Inequalities on Time Scales Has the belt and road initiative boosted the resident consumption in cities along the domestic route? – evidence from credit card consumption Analysis of the agglomeration of Chinese manufacturing industries and its effect on economic growth in different regions after entering the new normal Study on the social impact Assessment of Primary Land Development: Empirical Analysis of Public Opinion Survey on New Town Development in Pinggu District of Beijing Possible Relations between Brightest Central Galaxies and Their Host Galaxies Clusters and Groups Attitude control for the rigid spacecraft with the improved extended state observer An empirical investigation of physical literacy-based adolescent health promotion MHD 3-dimensional nanofluid flow induced by a power-law stretching sheet with thermal radiation, heat and mass fluxes The research of power allocation algorithm with lower computational complexity for non-orthogonal multiple access The art design of industrialised manufacturing furniture products based on the simulation of mathematical curves Research on the normalisation method of logging curves: taking XJ Oilfield as an example A Method of Directly Defining the inverse Mapping for a HIV infection of CD4+ T-cells model On the interaction of species capable of explosive growth Research on Evaluation of Intercultural Competence of Civil Aviation College Students Based on Language Operator Combustion stability control of gasoline compression ignition (GCI) under low-load conditions: A review Research on the Psychological Distribution Delay of Artificial Neural Network Based on the Analysis of Differential Equation by Inequality Expansion and Contraction Method The Comprehensive Diagnostic Method Combining Rough Sets and Evidence Theory Study on Establishment and Improvement Strategy of Aviation Equipment Design of software-defined network experimental teaching scheme based on virtualised Environment Research on Financial Risk Early Warning of Listed Companies Based on Stochastic Effect Mode System dynamics model of output of ball mill The Model of Sugar Metabolism and Exercise Energy Expenditure Based on Fractional Linear Regression Equation Constructing Artistic Surface Modeling Design Based on Nonlinear Over-limit Interpolation Equation Statistical analysis of typical elevator accidents in China from 2002 to 2019 Optimal allocation of microgrid using a differential multi-agent multi-objective evolution algorithm About one method of calculation in the arbitrary curvilinear basis of the Laplace operator and curl from the vector function Numerical Simulation Analysis Mathematics of Fluid Mechanics for Semiconductor Circuit Breaker Dynamics of infectious diseases: A review of the main biological aspects and their mathematical translation Cartesian space robot manipulator clamping movement in ROS simulation and experiment Effects of internal/external EGR and combustion phase on gasoline compression ignition at low-load condition Research of urban waterfront space planning and design based on children-friendly idea Characteristics of Mathematical Statistics Model of Student Emotion in College Physical Education Human Body Movement Coupling Model in Physical Education Class in the Educational Mathematical Equation of Reasonable Exercise Course Sensitivity Analysis of the Waterproof Performance of Elastic Rubber Gasket in Shield Tunnel Impact of Web Page House Listing Cues on Internet Rental Research on management and control strategy of E-bikes based on attribute reduction method A study of aerial courtyard of super high-rise building based on optimisation of space structure Exact solutions of (2 + 1)-Ablowitz-Kaup-Newell-Segur equation