Method for predicting the stress state of the lining of underground structures of quasi-rectangular and arched forms

. A method for predicting the stress-strain state of the lining of underground structures, the shape of the cross-section of which is different from the circular outline, is considered. The main task of the study is to develop a methodology for assessing the influence of the parameters of the cross-section shape of underground structures on the stress state of the lining. To solve this problem, a method for calculating the stress state of the lining for arched tunnels with a reverse arch and quasi-rectangular forms is substantiated and developed. The methodology was tested, which showed that the accuracy of the prediction of the stress state of the lining is sufficient to perform practical calculations. An algorithm for multivariate analysis of the influence of the cross-sectional shape of underground structures of arched and quasi-rectangular shapes on the stress state of the lining is proposed. Parametric calculations were performed using the developed algorithm and regularities of the formation of the stress state of the lining of underground structures for various engineering and geological conditions, as well as the initial stress state field, were obtained. A quantitative assessment of the influence of geometric parameters of tunnels on their stress-strain state was performed.

The presence of springs along the perimeter of the tunnel affects the stiffness of the considered system of rod elements. This is taken into account by changing the stiffness of each of the rod elements, consequently, the global stiffness matrix.
Methodology for calculating the loads on the lining of tunnels of various shapes. There are various hypotheses of the formation of rock pressure on the lining of underground structures, each of which has its own field of application, based on engineering and geological conditions and depending on the position of the underground structure relative to the Earth's surface. The applicability of methods for calculating loads on the lining has been studied by many scientists and is widely presented in the literature [32][33][34][35], therefore it is not the main object of this study. The method of calculating the load on the support proposed by Tertsagi [36] was adopted, which considered the formation of a descending rock pillar with a width of B1, part of the load from which is transferred to the surrounding rock mass due to friction, and the other component is held by the lining of the underground structure, which allows maintaining the equilibrium of the considered system of interaction "underground structuresoil mass".
According to the methodology proposed by Tertsagi, when the depth of the underground structure is less than two widths, the value of the vertical pressure on the lining of the underground structure is determined by the total weight of the rock pillar: where γaverage volume weight of overlying rocks; zthe distance from the Earth's surface to the design mark. For calculated cases when the depth of the underground structure is more than two widths, the calculation of the load on the lining is determined through the value of the effective power of the rocks h0, the weight of the rocks in which determines the amount of pressure: It is assumed that the descending rock pillar has a width B1, and sliding prisms are formed in the sides of the workings, the slope angle of which is assumed to be equal to (π/4 + φ/2). The value of the effective power of the descending rock pillar is determined by the formula where B1half the width of the descending rock pillar (Table 1);coupling; K0lateral pressure coefficient; Hdepth of the underground structure; P0additional load applied on the Earth's surface.
Formulas for determining the width B1 of underground structures of various cross-sectional shapes are summarized in Table 1. Features of preparation of the calculation model within the HRM method. When performing calculations of the stress state of the tunnel lining, the load transfer schemes indicated in Fig.1, a-c are adopted. The interaction between the lining of the underground structure described by the rod elements and the rock mass is carried out through springs oriented in the normal and tangential directions relative to the contour of the underground structure ( Fig.1, d). The stiffness of normal and tangential springs is defined as kn and ks, respectively. The values of kn and ks are determined through the values ηn and ηs, which characterize the ability of the host rock mass to resist deformation.
The nonlinear dependence determining the relation between the deformation of the lining δ and the reaction of the rock mass p to these deformations is presented in [13,30]: where plimthe limit value of the rock mass reaction; η0the initial value of the rock mass stiffness at a value of δ close to zero.
The indicator of the relative rock mass stiffness can be expressed as follows: The stiffness of the springs in the normal direction varies depending on the radius of the underground structure and is determined by the formula where Esmodulus of the rock mass deformation; vscoefficient of lateral deformation of the rock mass; Rithe radius of curvature of the lining at each considered point of its perimeter.
The limiting values of the reaction of the rock mass in the radial pn,lim and tangential ps,lim directions depend on the strength characteristics of the rocks (coupling and internal friction angle) of the magnitude characterizing the comprehensive compression of the rocks Δσconf : where Li -1distance between nodes i and i -1; Lidistance between nodes i and i + 1. When performing calculations, it was noted that the formation of the stress state of the lining of an underground structure is significantly influenced by its weight. Thus, in [37] it was shown that pressure reduction in the lower part of the underground structure allows for better convergence with the forecast results based on numerical modeling. When performing calculations, it is assumed that the value of the active pressure is reduced in relation to its maximum value by 10 % [38].
An algorithm for determining the shape of the cross-section of an underground structure by the force factor. The algorithm of multivariate analysis of the influence of the cross-section shape of a quasi-rectangular tunnel (Fig.2) on the stress state of the lining consists in splitting the tunnel section into separate segments, the parameters (the position of the center of the segment and the radius) of which are controlled through reference points. First of all, the cross-sectional area of the tunnel is determined in accordance with the size of the approach of buildings for a double-track subway tunnel [39,40]. The contour of the section in the light of the tunnel of a quasi-rectangular shape should overlap this zone. Taking into account this requirement, a number of mathematical dependencies are proposed to determine the geometric parameters of the quasi-rectangular shape of the tunnel section from the point of view of internal forces arising in the lining of the tunnel.
The following sequence of determining the geometric parameters of a quasi-rectangular tunnel is adopted (Fig.3).
At the initial stage, the geometric center O2 is selected and the radius R2 is calculated. The center O2 (x2, y2) of arc 2 is located on the bisector of the line segment A1A2, which is described by the linear equation The radius is determined by the formula where (d2, h1)coordinates of the point A2. (10) The location of the centers O1, O3 is determined. Arc 2 is bounded by the points 1 A and 2 A , located on a circle with the center O2 and radius R2. The angles between the straight line O2 1

A, O2
2 A and the horizontal axis must satisfy the condition: 1 = (0, 11) и 2 = (22, 90). The center O1 of arc 1 is the intersection of a line passing through point O2 and having a slope tg(1), and the horizontal x axis: O1(х1,0), A , respectively. Adaptation of the HRM method for a tunnel of non-circular cross-section shape was performed in the MatLab software package. The algorithm for calculating the parameters of the tunnel crosssection is also an element of the general program code for calculating the stress state of the lining. The output data are integral indicators of the stress state of the lining (shearing force, bending moment, axial force) and the position of the nodal points of the system after optimization.
A similar algorithm for assessing the influence of the cross-section shape of tunnels on the stress state of the lining has also been developed for tunnels of arched cross-section with a reverse arch (see Fig.1, b), but it is not considered in this paper.
Discussion of the results. Assessment of the reliability of the stress state prediction of the lining in the framework of the HRM method. The formation of the stress state of the lining of underground structures of circular, quasi-rectangular and arched cross-section forms, the parameters of which are shown in Fig.1 and in Table 2, is considered. The calculation of the stress state of the lining is carried out on the basis of the above methodology. The results obtained by the HRM method were compared with the results obtained on the basis of numerical simulation of the interaction of the lining with the rock mass by the finite element method. A detailed statement of the problem within the framework of the finite element method is presented in [38]. When performing calculations, the thickness of the lining was assumed to be equal to 0.5 m. The problem was solved in a plane-deformation formulation, the influence of the face zone on the formation of the load on the lining was not taken into account. The depth of the tunnel (from the ground surface to the crown of the arch of the underground structure) is 10 m, the lateral pressure coefficient for all design cases is 0.6. The rock mass was considered as a nonlinear deformable medium with a limit on the magnitude of maximum tangential stresses. Physical and mechanical characteristics of the rock mass, showing the mechanical behavior of rocks in the area under consideration, and the lining material : γ -18/25 kN/m 3 ; Е -3.6/35000 MPa; v -0.495/0.15; с -25.6/-kPa; φ -16.5/-respectively.
The results of the calculation of the stress state of the lining by the HRM method and its comparison with the results obtained based on the solution of the problem by the finite element method are presented in the form of diagrams of bending moments and axial forces along the perimeter of the (13) tunnel lining (Fig.4). The results show good convergence between different schemes of interaction of the rock mass and the lining of the underground structure. The greatest difference in the results between the two methods in terms of the bending moment is observed when considering a tunnel of circular cross-section and is 8.3 %, while for tunnels of quasi-rectangular and arched cross-sections the difference is 1.2 and 2.0 %, respectively. In terms of the magnitude of the axial forces, the difference between the two methods is 5.8, 5.4 and 6.0 %, respectively. From a qualitative point of view, there is a good convergence. The calculations performed with other parameters of the rock mass, support and stress state [38] also show good convergence of the two calculation methods.
It can be concluded that the presented method makes it possible to estimate the stress state of the lining with sufficient accuracy for practical calculations. Taking into account the fact that the speed of calculations within the HRM method is several orders of magnitude higher compared to numerical modeling methods, it can be effectively used to solve optimization problems, where the number of computational iterations can range from several hundred to tens of thousands.
The results of the evaluation of the influence of a curved tunnel cross-section on the stress state of the lining. Let's consider an example of the practical implementation of the proposed algorithm for assessing the influence of the tunnel cross-section shape on the stress state of the lining. Engineering and geological conditions for the construction of a quasi-rectangular cross-section tunnel were adopted for Hanoi, Vietnam. The size of the approach of the tunnel structures for the conditions under consideration is shown in Fig.5. The depth of the tunnel is 10 m, the thickness of the lining is 0.5 m. The geometric parameters of the cross-section are calculated according to the algorithm described earlier.
To reduce the calculation time, each calculation case is estimated using a step of 1 for angles 1 and 2 from bounded angles 11 and 22 respectively; 0.1 m to change the parameter x2. When calculating the geometric parameters of the tunnel cross-section, 12479 calculations of the stress state of the lining of the quasi-rectangular cross-section were performed, the estimated time for each case was about two seconds.
The dependence of internal forces and radial displacements of the lining on the radius of the reverse arch R3 of the tunnel is shown in Fig.6. The influence of the radius R3 on the value of the maximum bending moment in the cross section of the lining is significant. An increase in the radius R3 leads to an increase in the maximum bending moment in the lining. The maximum normal force changes slightly when the radius R3 changes.
Also, an increase in the radius R3 leads to the fact that the arch and the reverse arch of the tunnel become flatter. Then the normal vertical load acting on these segments of the lining increases, and the vertical lateral load decreases. The vertical load transmitted to the side wall of the tunnel is also Processing of the calculation results allowed to establish that the ratio between the radius R3 and the maximum bending moment in the lining of the tunnel of a quasi-rectangular cross-section shape can be described by the equation The relation between the maximum bending moment in the tunnel lining and the radii R2, R3, R1 with sufficient accuracy for practical calculations can be described by the equation  (14) and (15) that the maximum bending moment in the lining is mainly affected by the radius of the reverse arch and slightly by the radius of the side wall of the tunnel.
A detailed analysis of the calculation results allowed to establish that the value of the maximum bending moment in the section of the tunnel lining for the conditions under consideration varies from 1.05 (the best shape of the cross-section according to the force factor) to 1.85 MN·m/m (the worst shape), i.e. the difference in the values of the bending moment is 75 %. At the same time, the difference in area between the best and worst forms of cross-section was 21 %. It can be concluded that the analysis of the data obtained on the basis of a multivariate calculation allowed to significantly reduce the stress state of the lining while maintaining the minimum size of the cross-section area of the tunnel.
To assess the influence of various parameters and geometry of the tunnel on the stress state of the lining, a multivariate calculation was performed. In a certain range, the modulus of deformation of the rock mass, the depth of the tunnel and the lateral pressure coefficient changed. The modulus of soil deformation is one of the main parameters determining the stress state of the tunnel lining. Its value during the calculations varied from 5 to 200 MPa, while the remaining physical and mechanical characteristics of the rocks remained the same. A total of 26031 calculations were performed. The results are presented in the form of dependences of the change in the maximum value of the bending moment and the axial force in the lining of the tunnel (Fig.7, a). As can be seen from the presented dependencies, the following relation is true for all values of the deformation modules of the rock mass: the higher the value of R3, the higher the absolute value of the bending moment; at the same time, the magnitude of the axial force varies slightly. The obtained results correlate well with the data presented in [41].
With a change in R3 from 10 to 15 m, a sharp increase in the bending moment values is observed, while the values of the maximum normal forces decrease rapidly. A further increase in the R3 index does not have a serious effect on the stressed state of the lining.
The obtained results can be generalized in the form of a dependence that characterizes the relation between the maximum value of the bending moment, geometric parameters and the modulus of deformation of the rock mass: The assessment of the influence of the initial stress state of the rock mass, expressed in terms of the value of the lateral pressure coefficient K0 and the depth of the tunnel, is presented in the form of plots of changes in the integral indicators of the stress state of the lining (Fig.7, b). The value of the lateral pressure coefficient K0 varied from 0.4 to 1.0. The number of calculations performed was 22778. Based on the analysis of the calculation results, it can be concluded that an increase in K0 leads to a decrease in the maximum value of the bending moment, but the value of the axial force as a whole increases. This tendency is typical for any values of R3. When R3 reaches a value of 10 m, the influence of the coefficient K0 on the stress state of the lining increases. With an increase in the value of R3, the compression effect decreases, therefore, the influence of the R3 indicator on the stressed state of the lining also decreases. According to the results of calculations, a dependence is obtained for determining the maximum value of the bending moment through the value R3 and the lateral pressure coefficient K0: The lateral pressure coefficient K0 also influences the formation of normal forces in the lining section, however, the degree of this influence is significantly lower than the influence on the bending moment. This result was obtained at different values of R3.
To assess the effect of the radius R3 on the stress state of the lining at different depths of the laying tunnel, additional calculations were performed in the amount of 19524 (Fig.7, с). The depth of laying H varied from 5 to 50 m. Based on the performed studies, it was found that the influence of the parameter R3 on the bending moment decreases as the depth of the tunnel increases. This can be explained by the fact that with increasing depth, the active loads on the lining also increase, which leads to an increase in the values of the integral indicators of the stress state of the lining. An increase in the R3 indicator leads to the fact that the reverse arch becomes flatter, and the ratio between normal and tangential loads increases, which reduces the vertical load in the tunnel walls. To calculate the value of the maximum bending moment in the lining, depending on the depth of the tunnel and the geometric parameters of the tunnel, the formula is used Below is a formula that takes into account the influence of a combination of the factors considered: soil deformation modulus Es, tunnel depth H, soil pressure coefficient on the maximum bending moment in the lining. The equation is obtained based on the analysis of 956675 cases: Mmax = 1.10608 -0.26759lnEs -1.06703K + 0.03098H -0.01048lnR1 + + 0.00395expR2 + 0.50104lnR3 MN·m/m.
As can be seen from the presented results, the proposed algorithm makes it possible to effectively assess the influence of the controlling parameters of the tunnel section on the stress state of the lining of various cross-section shapes. At the same time, an automated method for optimizing the shape of the tunnel according to the specified criteria (material strength, structural stability, cross-section area, etc.) has not been proposed, but will be studied in further studies.
Conclusion. As the practice of construction of underground structures shows, their cross-section shapes can be circular, arched with a reverse arch, quasi-rectangular, etc. The HRM method can be used to calculate the stress state of the lining of non-circular shapes, which allows solving the problem with sufficient accuracy for practical calculations. Taking into account the fact that the calculation speed is several orders of magnitude higher than in the finite element method, where the rock mass is considered as a continuous medium described by the corresponding model of its deformation, the HRM method can be used to solve optimization problemsfinding the optimal shape of the crosssection of an underground structure under given mining and geological conditions. Based on the studies performed: • the HRM method has been developed for calculating underground structures of non-circular cross-section shape taking into account the nonlinear behavior of the rock mass; • an algorithm for the multivariate calculation of the stress state of the lining of tunnels of arched and quasi-rectangular cross-section forms has been developed; • the proposed method has been tested for various tunnel construction conditions, and good convergence with the results of numerical modeling obtained in the PLAXIS 2D software package has been obtained; • a set of parametric calculations has been performed and dependences have been obtained linking the integral parameters of the stress state of the lining (bending moment and axial forces) with the modulus of deformation of the rock, the stress state of the rock mass and the geometric parameters of the cross-section of tunnels of quasi-rectangular and arched shapes; • analytical expressions are obtained for calculating the maximum values of bending moments in the lining cross-section; • the direction of further research on the development of an algorithm for finding the optimal shape of the cross-section of tunnels according to the specified criteria is indicated.