Engineering, Technology & Applied Science Research Vol. 7, No. 2, 2017, 1504-1514 1504 www.etasr.com Rastbood et al.: Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel… Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel Lining Armin Rastbood School of Mining Engineering College of Engineering, University of Tehran Tehran, Iran arminrastbood88@gmail.com Yaghoob Gholipour School of Civil Engineering College of Engineering, University of Tehran Tehran, Iran ygholipour@ut.ac.ir Abbas Majdi School of Mining Engineering College of Engineering, University of Tehran Tehran, Iran amajdi@ut.ac.ir Abstract—The main objective of this paper is to optimize the geometrical and engineering characteristics of concrete segments of tunnel lining using Finite Element (FE) based Response Surface Methodology (RSM). Input data for RSM statistical analysis were obtained using FEM. In RSM analysis, thickness (t) and elasticity modulus of concrete segments (E), tunnel height (H), horizontal to vertical stress ratio (K) and position of key segment in tunnel lining ring (θ) were considered as input independent variables. Maximum values of Mises and Tresca stresses and tunnel ring displacement (UMAX) were set as responses. Analysis of variance (ANOVA) was carried out to investigate the influence of each input variable on the responses. Second-order polynomial equations in terms of influencing input variables were obtained for each response. It was found that elasticity modulus and key segment position variables were not included in yield stresses and ring displacement equations, and only tunnel height and stress ratio variables were included in ring displacement equation. Finally optimization analysis of tunnel lining ring was performed. Due to absence of elasticity modulus and key segment position variables in equations, their values were kept to average level and other variables were floated in related ranges. Response parameters were set to minimum. It was concluded that to obtain optimum values for responses, ring thickness and tunnel height must be near to their maximum and minimum values, respectively and ground state must be similar to hydrostatic conditions. Keywords-tunnel; segment; lining; response surface methodology; analysis of variance I. INTRODUCTION Nowadays, most number of the tunnels are excavated using shield tunneling method in soft ground conditions, because this type of tunneling is rapid, cost effectiveness and has minimum effect on surrounding environment. In general, this type of tunneling uses concrete segments as a tunnel support system. Assembling concrete segments to each other in one excavation step forms a structure named ring [1, 2]. Usually, the number of erected segments in a ring is more than four. Successive rings compose a final support system of the tunnel, called lining. One of the segments in a ring is designed usually smaller than the other segments and is called the key segment. This causes the facility both in ring erection and excavation operation in curved alignments. Key segment will be installed at the end of ring assembling operation. In Figures 1 and 2 some nomenclatures of tunnel lining support system are illustrated. Fig. 1. Assembled ring in segment manufacturing factory, Tabriz metro- Line 1 Fig. 2. Tunnel lining structure Concrete segments have a complicated process from the design and construction points of view [3]. The final design of concrete segments must provide optimized ring bearing capacity against different types of loads both in Ultimate Limit State (ULS) and Serviceability Limit State (SLS) [4, 5]. In [6], an experimental equation to calculate the reduced effective moment of inertia of a ring with joints was proposed. The proposed equation includes segment number and joint stiffness parameter: 4nI,I.I n 4 II e 2 je        (1) Where Ie is effective moment of inertia of jointed ring, Ij is moment of inertia in the joints, I is the moment of inertia of a continuous ring without joints and n is the number of segments. Engineering, Technology & Applied Science Research Vol. 7, No. 2, 2017, 1504-1514 1505 www.etasr.com Rastbood et al.: Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel… In [7], authors discussed about the main differences in the assumptions for three different models: the continuum ring model, the Muir Wood design model and the bedded beam model without bedding at the tunnel lining crown region. In [8], authors presented a comparative study between spatial (3D) computing static models and theoretical results of commonly used plain (2D) computing models for segmental tunnel lining. Comprehensive methods for determination of structural forces in tunnel lining and segment design were presented in [9]. In [10], authors categorized all tunnel lining design methods into four major types: (a) empirical design methods; (b) experimental and laboratory modelling; (c) circular ring in elastic foundation model; and (d) continuum mechanics models including analytical methods and numerical methods. Four most important models of tunnel lining structure by considering ring and segment joints are illustrated in [11]. In [12], authors discussed about the various load types imposed on shield machine and tunnel lining. In [13], authors performed a full scale test on three-rings of segmental lining structure of the Shanghai Changjiang tunnel with a diameter of 15m. In [14], author conducted a series of tests on small laboratory segmental tunnel lining were made from PVC to obtain their load bearing capacity. In this study authors focused on simple beam support tests to estimate joint stiffness. A comprehensive numerical model using 2D finite difference element model was proposed in [15]. In this paper many factors that influence tunnel lining behavior together with three different types of joint springs, i.e. rotational, axial and radial springs were modeled. Some analytical solution are available for segmental tunnel lining design, as those proposed in [7, 16-19]. In [10], authors introduced an analytical equation for moment reduction factor based on the maximum horizontal displacement of a uniform ring. An analytical solution for analyzing tunnel lining in longitudinal direction was proposed in [20]. In [21], authors studied the elastic analysis of a circular lined tunnel by considering the delayed installation of the tunnel support. Recently, many publications about numerical modelling of segmental tunnel lining both in 2D and 3D, considering more detailed aspects of shield tunneling were published [15, 22-25]. Despite most comprehensive available attempts based on different approaches to analyze segmental tunnel lining in various point of view, there are no considerable publications in this field using statistical methods. In this paper, FEM based Response Surface Methodology (RSM) approach was organized to establish statistical equations between tunnel lining characteristics as input variables and generated stresses and displacements in tunnel lining structure as responses parameters using analysis of variance (ANOVA). Finally optimization analysis was performed for each response parameter in terms of input variables. II. NUMERICAL MODELLING In this study, input data for statistical analysis, called experiments from now on, were obtained using the finite element method (ABAQUS [26]). In all FEM models, one ring of tunnel lining support system composed from 5+1 segments. Table I shows the engineering and geometrical characteristics of the concrete segments. The behavior of the tunnel lining is assumed to be elastic. In a ring, the shape of the five segments (A2-A6) are similar to each other except for the key segment (A1) which is smaller. In all numerical models, it is assumed that the ring is positioned far away from the tunnel face and is not influenced by shield machine loads and tunneling operation. Beam-spring model introduced in [27] was used in numerical modeling. This type of segmental tunnel lining is proposed in several tunneling guidelines [9, 19]. Interaction of ground on outer surface of tunnel lining was simulated using tangential and normal springs. Stiffness of ground normal springs is evaluated using (2), [28, 29]:   ν1R.A.EK n  (2) Where Kn is the stiffness of radial spring, E and ν are the elasticity modulus and poison’s ratio of soil respectively, R is the tunnel radius, and A represents the effective area of tunnel lining that is subjected to implied force from the soil, and is calculated using (3): .bRA  (3) Where θ is radial angle between two successive radial springs applied on lining surface, and b represents effective area of each spring in tunnel longitudinal direction. The stiffness of tangential springs (Kt) is assumed to be one-third of the normal spring stiffness (Figure 3) [34]. The ring under the impact of interaction springs and surrounding ground load is shown in Figures 4a and 4b respectively. 3D solid-stress elements with linear geometric order were used to model concrete segments (Figure 5a). Plane strain condition was considered in the models. Also according to the literature, in transverse direction the segmental lining structure is usually designed in plane strain condition [9, 30-33]. In this study, it is assumed that origin of angle in transverse section of the model is positioned at tunnel crown (Figure 5b). Longitudinal joints of assembled segments in a ring and key segment position at θ=90° are shown in Figure 5c. Hard contact was assumed for interaction of six concrete to concrete contact surfaces in segment joints with frictional penalty coefficient of 0.4. Fig. 3. Constitutive relationship of the springs representing the ground reaction in the bedded-spring model. Pn and Pt represent the normal and tangential load, δn and δt are the normal and tangential displacement; Kn and Kt are the normal and tangential stiffness A. Input data preparation for statistical RSM analysis The numerical model first was solved for t (concrete segment thickness) =30 cm, tunnel overburden H= 5m, K=0.5 (horizontal to vertical stress ratio), E(lining) =20 GPa and θ=0° Engineering, Technology & Applied Science Research Vol. 7, No. 2, 2017, 1504-1514 1506 www.etasr.com Rastbood et al.: Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel… (key position at crown). Then t, H, K and E values were kept constant and θ value changed 30°, 60°, 90°, 120°, 150° and 180° respectively. To prepare enough input data for statistical RSM, for each value of input variables, values of other input variables were changed according to the rule presented in Table II. These variables considered as input independent variables in statistical analysis. This table shows times of changes for each input variable and its value. Due to the axisymmetric shape of the ring, key positions at 210°, 240°, 270°, 300° and 330° at lining periphery were neglected. Finally 252 numerical models were analyzed. TABLE I. ENGINEERING AND GEOMETRICAL PROPERTIES OF CONCRETE SEGMENTS Engineering properties Geometrical properties Segment No. E* ν** ρ *** t **** Central angle(°) A1(key segment) 30 A2-A6 20 0.15 2350 30 66 *elasticity modulus (GPa), **Poisson ratio, *** density (kg/m3), **** thickness(cm) Fig. 4. Tunnel lining under ground springs and load: (a) ground radial and tangential springs, (b) surrounding ground load imposed on tunnel lining Commonly, Tresca and Von Mises yield criteria are used as failure criteria for materials. According to Tresca yield criterion, material begins to yield when maximum absolute value of shear stress reaches to a critical value, and based on Von Mises yield criterion, material begins to yield when the second deviatoric stress invariant approaches to a critical value. So, due to importance of Tresca and Mises stresses, these both type of stresses together with ring displacement are considered for prediction analysis. Figures 6a, 6b and 6c show extreme values resulted for Von Mises and Tresca stresses, and ring displacements for t=30cm, H=15m, K=0.5, E=20 GPa and θ=0°. These output parameters considered as responses in statistical analysis. Fig. 5. Assembled ring of concrete segments: (a) Meshed Model, (b) Origin of θ angle, (c) Longitudinal joints and key position TABLE II. DIFFERENT VALUES FOR INPUT VARIABLES Different values of 5 Input Variables t (cm) H(m) K E(GPa) Key Position(°) 30 40 5 15 25 0.5 1.0 1.5 20 35 0 30 60 90 120 150 180 Fig. 6. Extreme values and distribution of (a) Von Mises stresses (N/m2), (b) Tresca stresses (N/m2) and (c) ring displacements (m) for t=30cm, H=15m, K=0.5, E=20 GPa and θ=0°. III. RESPONSE SURFACE METHODOLOGY Response surface methodology (RSM) is a collection of statistical and mathematical techniques useful for developing, improving, and optimizing processes [35]. In most RSM problems, the form of the relationship between the response and the independent variables is unknown. Thus, the first step in RSM is to find a suitable approximation for the true functional relationship between y and the set of independent variables [36]. Linear or square polynomial functions are employed to explain the considered problem. If there is (a) (b) (a) (b) (c) Engineering, Technology & Applied Science Research Vol. 7, No. 2, 2017, 1504-1514 1507 www.etasr.com Rastbood et al.: Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel… curvature in the system, then a polynomial of higher degree must be used, such as the second-order model in (4):    εxxβxβxββy jiij k 1i 2 iiii k 1i i0 (4) Where y is the response, k is the number of variables, β0 is the constant term, βi represents the coefficients of the linear parameters, βii represents the coefficients of the quadratic parameter, βij represents the coefficients of the interaction parameters, xi represents the variables, and ε is the residual associated to the experiments [36]. RSM based statistical methodology is capable of optimizing the responses in terms of input variables. A flow chart of the RSM general procedure is shown in Figure 7. Quadratic regression equation was applied in RSM to establish equations between five input variables and three responses. In RSM, a series of analysis were conducted by Historical Data Design (HDD) using Design Expert 10.0.1 (State Ease, Inc., Minneapolis, MN, USA). Historical Data design enables the researcher to use data from a previous experiment or investigation. Any set of numeric or categoric data can be analyzed via Historical Data Design. As explained previously, 252 experiments conducted to provide data for statistical analysis. Some sample values of five factors and three responses are shown in Table III. Both for input variables and responses, their range, average value, ratio of maximum to minimum value, and standard deviation were as mentioned Tables IV and V, respectively. TABLE III. PREPARED EXPERIMENTS FROM FINITE ELEMENT MODELING Factor 1 Factor 2 Factor 3 Factor 4 Factor 5 Response 1 Response 2 Response 3 Run A:Thickness (cm) B:Height (m) C:K Ratio (---) D:E Modulous (Gpa) E:Theta (Degree) MisesMax (N/m 2 ) TrescaMax (N/m 2 ) UMax (m) 1 30.00 5.00 0.50 20.00 0.00 3.003E+006 3.269E+006 0.005132 2 30.00 5.00 0.50 20.00 30.00 2.768E+006 3.019E+006 0.005211 3 30.00 5.00 0.50 20.00 60.00 3.099E+006 3.381E+006 0.005242 … … … … … … … … … 252 40.00 25.00 1.50 35.00 180.00 3.222E+007 3.465E+007 0.02388 TABLE IV. INPUT INDEPENDENT VARIABLES Name Units Minimum Maximum Mean Standard deviation. Segment Thickness (A) cm 30.00 40.00 35.00 5.01 Height of tunnel (B) m 5.00 25.00 15.00 8.18 K (stress ratio) (C) --- 0.50 1.50 1.00 0.41 E (Elasticity Modulous) (D) GPa 20.00 35.00 27.50 7.51 θ (Key position) (E) Degree 0.00 180.00 90.00 60.12 TABLE V. THREE RESPONSES OF RSM ANALYSIS Response Name Units Analysis Minimum Maximum Mean Standard Deviation. Ratio Mises Max stress value N/m 2 Polynomial 1.977E+006 6.05E+007 1.72815E+007 1.37813E+007 30.6019 Tresca Max stress value N/m 2 Polynomial 2.115E+006 6.896E+007 1.91188E+007 1.55427E+007 32.6052 Displacement Max value m Polynomial 0.004799 0.02973 0.0142558 0.00811969 6.19504 Fig. 7. FEM based RSM flow chart Engineering, Technology & Applied Science Research Vol. 7, No. 2, 2017, 1504-1514 1508 www.etasr.com Rastbood et al.: Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel… A. Analysis of Variance (ANOVA) The fitted mathematical equation to the data using RSM may not adequately explain the considered experiment. Analysis of variance (ANOVA) is a method to assess the reliability of fitted function. The significance of regression can be evaluated by the ratio between the mean of the square of regression (MSreg) and the mean of the square of residuals (MSres), considering their related degrees of freedom. The higher value of this ratio (F-value) represents that the statistical model is fitted properly to the experimental data. In this study, analysis of variance was used to evaluate the influence of each independent variable and their interactions on responses. Statistical parameters of three proposed models for maximum Misses stress, maximum Tresca stress and maximum ring displacement are illustrated in Table VI. The higher the F- value, the better the significance of the model. “Adequate Precision” measures the signal to noise ratio. A ratio greater than 4 is desirable. For three responses, actual data against predicted ones are shown in Figure 8. The evaluation of the fitted models is performed using the regression coefficient. Presented models approximately show satisfactory regression values. Also the predicted R-squared is in reasonable agreement with the adjusted R-squared. Normal probability plot of the studied residuals for triple responses are illustrated in Figure 9. These plots are to check the normality of residuals. TABLE VI. STATISTICAL PARAMETERS OF PRESENTED MODELS FOR THREE RESPONSES Statistical parameter Maximum Misses stress Maximum Tresca stress Maximum ring displacement Description F-value* 178.15 173.01 397.62 Significant model Adequate Precision 46.953 46.546 49.470 Adequate signal R-Squared 0.8135 0.8090 0.8656 Predicted R-Squared 0.8021 0.7974 0.8597 Adjusted R-squared 0.8090 0.8044 0.8634 The Predicted R-Squared is in reasonable agreement with the Adjusted R-Squared * There is only a 0.01% chance that an F-value this large could occur due to noise. Fig. 8. Actual and predicted values for maximum Mises Stress, Maximum Tresca stress and Maximum ring displacement: (a) Maximum Mises stress (b) Maximum Tresca stress (c) Maximum ring displacement Fig. 9. Plot of Normal probability versus studentized residuals : (a) Maximum Mises stress (b) Maximum Tresca stress (c) Maximum ring displacement (a) (b) (c) Engineering, Technology & Applied Science Research Vol. 7, No. 2, 2017, 1504-1514 1509 www.etasr.com Rastbood et al.: Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel… If residual points follow a straight line, the normal probability plot indicates that the residuals are normally distributed. The more linear the shape of the normal probability plot, the better the quality of the model. B. Second-order polynomial models In analysis of variance, values of "Prob > F" less than 0.0500 indicate model terms are significant. For maximum Mises stress case, A, B, C, AB, B2, C2 are significant model terms, Table VII. Values greater than 0.1000 indicate the model terms are not significant. So the included variables in the equation are tunnel lining thickness, tunnel height and K ratio. Final equation in terms of actual factors for Mises maximum stress is as follows (5): Sqrt (MisesMax) =7052-37.4*(Thickness)+280.2*(Height)- 11977.4*(K Ratio)-2.16*(Thickness * Height)-3.03*(Height)2 +6766*(K Ratio)2 (5) For maximum Tresca stress case, similar to maximum Mises stress, A, B, C, AB, B2, C2 are significant model terms, Table VIII. Final equation in terms of actual factors for Tresca maximum stress is as follows (6): Sqrt (TrescaMax) =7666.685-43.982*(Thickness)+ 300.087*(Height)-12816.724*(K Ratio)-2.464*(Thickness* Height)-3.170*(Height)2+7232.985*(K Ratio)2 (6) For maximum ring displacement case, B, C, BC, C2 are significant model terms, Table IX. So the included variables are tunnel height and K ratio. Final equation in terms of actual factors for maximum ring displacement is: Sqrt (UMax) =0.189+4.215E-003*(Height)-0.34*(K Ratio)- 2.53E-00*(Height * K Ratio)+0.206*(K Ratio)2 (7) C. Response surface analysis of maximum Mises stress, Tresca stress and tunnel ring displacement The individual and interaction influence of included variables in equations (5) to (7) on responses are shown in Figures 10 to 23. It must be noticed that in the analysis of the influence of each variable on the response, other input parameters are kept in their average values. TABLE VII. ANOVA FOR MAXIMUM MISES STRESS Analysis of variance table Source Sum of Squares df Mean Square F Value p-value Prob > F Model 5.157E+008 6 8.596E+007 178.15 < 0.0001 si g n if ic an t A- Thickness 3.070E+007 1 3.070E+007 63.64 < 0.0001 B-Height 2.162E+008 1 2.162E+008 448.10 < 0.0001 C-K Ratio 1.015E+008 1 1.015E+008 210.33 < 0.0001 AB 1.962E+006 1 1.962E+006 4.07 0.0448 B2 5.169E+006 1 5.169E+006 10.71 0.0012 C2 1.602E+008 1 1.602E+008 332.08 < 0.0001 TABLE VIII. ANOVA FOR MAXIMUM TRESCA STRESS Analysis of variance Source Sum of Squares df Mean Square F Value p-value Prob > F Model 5.836E+008 6 9.727E+007 173.01 < 0.0001 si g n if ic an t A-Thickness 4.128E+007 1 4.128E+007 73.42 < 0.0001 B-Height 2.368E+008 1 2.368E+008 421.19 < 0.0001 C-K Ratio 1.142E+008 1 1.142E+008 203.20 < 0.0001 AB 2.550E+006 1 2.550E+006 4.54 0.0342 B2 5.630E+006 1 5.630E+006 10.01 0.0018 C2 1.831E+008 1 1.831E+008 325.69 < 0.0001 TABLE IX. ANOVA FOR MAXIMUM LINING RING DISPLACEMENT Analysis of variance Source Sum of Squares df Mean Square F Value p-value Prob > F Model 0.26 4 0.065 397.62 < 0.0001 significant B-Height 0.048 1 0.048 289.93 < 0.0001 C-K Ratio 0.048 1 0.048 291.23 < 0.0001 BC 0.018 1 0.018 108.81 < 0.0001 C2 0.15 1 0.15 900.50 < 0.0001 According to Figure 10, as tunnel lining thickness increases in its range, the maximum Mises stress decreases slightly. Maximum Mises stress increases when tunnel height value increases from 5.0 m up to 25 m, Figure 11. By increasing K ratio from 0.5 to unity (hydrostatic condition), maximum Mises stress decreases, and consequently increases when K ratio increases from unity up to 1.5, Figure 12. As can be seen from Figure 13, simultaneously increasing both in lining thickness and tunnel height variables cause reduction in maximum Mises stress, but tunnel height variable has more influence on maximum Mises stress reduction than lining thickness variable. In a constant value of lining thickness, as K stress ratio increases up to unity, the maximum Mises stress decreases, and then increase mutually with K ratio. This trend can be seen in all constant values of lining thickness. This effect is shown in Figure 14 three dimensionally. Analysis of variance showed that regression equations of maximum Mises stress and maximum Tresca stress, both include the same terms. As illustrated in Figures 15-18, lining thickness, tunnel height, K ratio and thickness-height interaction have the same effect on maximum Tresca stress as those effect on maximum Mises stress. In a constant value of tunnel height, as lining thickness increases from 30 cm up to 40 cm, maximum Tresca stress decreases slightly, but the rate of reduction in maximum Tresca stress is greater in higher values of tunnel height (Figure 19). For maximum displacement of a tunnel lining ring (Umax), as tunnel height increases, maximum tunnel lining displacement increases (Figure 20). By increasing K ratio from 0.5 to unity (hydrostatic condition), maximum ring displacement decreases and consequently, increase mutually when K ratio increases from unity to 1.5 (Figure 21). Engineering, Technology & Applied Science Research Vol. 7, No. 2, 2017, 1504-1514 1510 www.etasr.com Rastbood et al.: Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel… Fig. 10. The effect of lining thickness on maximum Mises thickness Fig. 11. The effect of tunnel height on maximum Mises stress Fig. 12. The effect of K ratio on maximum Mises stress Fig. 13. The effect of thickness and height interaction in 2D plane on maximum Mises stress Fig. 14. The effect of thickness and height interaction in 3D space on maximum Mises stress Fig. 15. The effect of lining thickness on maximum Tresca stress Engineering, Technology & Applied Science Research Vol. 7, No. 2, 2017, 1504-1514 1511 www.etasr.com Rastbood et al.: Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel… Fig. 16. The effect of tunnel height on maximum Tresca stress Fig. 17. The effect of K ratio on maximum Tresca stress Fig. 18. The effect of thickness and height interaction in 2D plane on maximum Tresca stress Fig. 19. The effect of thickness and height interaction in 3D space on maximum Tresca stress Fig. 20. The effect of tunnel height on maximum lining ring displacement Fig. 21. The effect of tunnel height on maximum lining ring displacement Engineering, Technology & Applied Science Research Vol. 7, No. 2, 2017, 1504-1514 1512 www.etasr.com Rastbood et al.: Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel… Fig. 22. The effect of tunnel height and k ratio interaction in 2D plane on maximum lining ring displacement As it can be seen from Figure 22, increasing simultaneously both in tunnel height and K ratio variables increases the maximum ring displacement, but tunnel height variable has more influence to increase the maximum ring displacement than K ratio variable. In all constant values of K ratio, as tunnel height increases, the maximum ring displacement increases slightly. This effect is shown dimensionally in Figure 23. Fig. 23. The effect of tunnel height and k ratio interaction in 3D space on maximum lining ring displacement D. Optimization Structural optimization based on FEM requires special computational tools. This attempt has become a major part in this field of study, where optimization seems to be a mandatory process in design and economical points of view. Optimization using numerical methods needs considerable computational cost and most likely is time consuming based on the complexity of the problem. Application of FE based response surface methodology is an alternative solution to reduce computational cost and time of analysis considerably without deficiency in the models. Therefore, to obtain an optimum structure for segmental tunnel lining both from geometrical and engineering point of view, at first constraints, goals and importance order of each input and response parameter were set according to Table X. Since the analysis of variance showed that elasticity modulus and key segment position variables are not involved in response equations, so their values were kept constant in average level during optimization process. Domain of other input variables, i.e. tunnel lining thickness, tunnel height and K ratio were set into their ranges. To optimize the tunnel lining structure, response parameters, i.e. maximum Mises and Tresca stresses and maximum lining ring displacement were set to their minimum values. After optimization process, obtained results are presented in Table XI. TABLE X. GOALS AND CONSTRAINTS OF VARIABLES AND RESPONSES Constraints Name Goal Lower Limit Upper Limit Importance A:Thickness is in range 30 40 3 B:Height is in range 5 25 3 C:K Ratio is in range 0.5 1.5 3 D:E Modulus average 20 35 3 E:Theta average 0 180 3 MisesMax minimize 1.977E+006 6.05E+007 3 TrescaMax minimize 2.115E+006 6.896E+007 3 UMax minimize 0.004799 0.02973 5 One hundred solutions were proposed by the statistical analysis. Although, desirability index was equal to unity in all proposed solutions, but the best solution is underlined in first row of Table XI. Plots of ramp display for best solution is shown in Figure 24. Optimum values of variables and responses are shown. TABLE XI. PROPOSED SOLUTIONS TO OPTIMIZE TUNNEL LINING CHARACTERISTICS Solutions: No. Thickness (cm) Height (m) K Ratio (--) E Modulous* (Gpa) Theta* (°) MisesMax (N/m2) TrescaMax (N/m2) UMax (m) Desirability Index 1 38.595 5. 0.864 27.5 90 1972897.422 2104034.290 0.004 1 2 39.9 5.3 0.88 27.5 90 1960092.90 2069289.145 0.0037 1 3 39. 5.07 0.853 27.5 90 1958046.11 2082296.28 0.0036 1 4 39.7 5.08 0.94 27.5 90 1911454.36 2020534.00 0.0038 1 5 39.7 5.02 0.80 27.5 90 1951520.38 2068656.82 0.0037 1 6 39.5 5.10 0.84 27.5 90 1917118.97 2031757.06 0.0036 1 … … … … … … … … … … 100 34.4 5.0 0.852 27.5 90 2502924.69 2740456.29 0.0036 1 Engineering, Technology & Applied Science Research Vol. 7, No. 2, 2017, 1504-1514 1513 www.etasr.com Rastbood et al.: Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel… Desirability=1.0 Fig. 24. Optimized values: red and blue circles represent input variables and responses, respectively. From the optimization process it is concluded that to obtain minimum values for each three response parameters, ring thickness and tunnel height must be near to their maximum and minimum values, respectively and the tunnel lining must be embedded nearly in hydrostatic condition of the ground, i.e. a K ratio of 1. On the other hand, the higher values for tunnel lining thickness together with the lower values for tunnel height simultaneously, will induce both the minimum values for Mises and Tresca stress values and the displacement of tunnel lining ring in nearly hydrostatic ground conditions. IV. CONCLUSIONS FEM based response surface methodology was applied to prepare prediction models for segmental tunnel lining. Quadratic regression equation was applied in RSM to establish equations between five input variables and three responses. Obtained models were used to determine the maximum Mises and Tresca stresses and lining ring displacement generated in one ring of tunnel lining. Analysis of variance was used to discuss about the influence of each independent input variable and their interactions on response parameters. For the three responses, the evaluation of the fitted models are performed by the regression coefficient. Presented models approximately had acceptable values of regression. The terms in maximum Mises and Tresca stresse models were the same: concrete segment thickness (t), tunnel overburden (H) and stress ratio (K). The terms in maximum ring displacement model were tunnel overburden (H) and stress ration (K), i.e. concrete segment thickness and elasticity modulus of concrete segments and key position had no considerable effect on maximum ring displacement. It is found that tunnel height variable influenced the responses more than both segment thickness and K ratio variables in stress and displacement models, respectively. In the optimization process, it is concluded that to obtain minimum values for responses, ring thickness and tunnel height must be near to their maximum and minimum values, respectively and the tunnel must be embedded nearly in hydrostatic condition of the ground, i.e. a K ratio of 1. REFERENCES [1] V. Guglielmetti, P. Grasso, A. Mahtab, S. Xu, Mechanized tunnelling in urban areas: design methodology and construction control, CRC Press, 2008 [2] J. B. O'Carroll, A Guide to Planning, Constructing, and Supervising Earth Pressure Balance TBM Tunneling, Parsons Brinckerhoff, 2005 [3] R. G. A. de Waal, Steel fibre reinforced tunnel segments-for the application in shield driven tunnel linings, TU Delft, Delft University of Technology, 2000 [4] C. B. M. Blom, Design philosophy of concrete linings for tunnels in soft soils, TU Delft, Delft University of Technology, 2002 [5] A. Majdi, H. Ajamzadeh, K. Emad, “Numerical evaluation of the influence of bolted-and non-bolted joints on segmental tunnel lining behavior in line 4 Tehran subways”, 6th Asian Rock Mechanics Symposium - Advances in Rock Engineering, New Delhi, India, October 23-27, 2010 [6] A. M. Wood, “The circular tunnel in elastic ground”, Geotechnique, Vol. 25, No. 1, pp. 115-127, 1975 [7] H. Duddeck, J. Erdmann, Structural design models for tunnels in soft soil, Underground Space, United States, 1985 [8] I. Hudoba, “Contribution to static analysis of load-bearing concrete tunnel lining built by shield-driven technology”, Tunnelling and Underground Space Technology, Vol. 12, No. 1, pp. 55-58, 1997 [9] I. W. Group, “Guidelines for the design of shield tunnel lining”, Tunnelling and Underground Space Technology, Vol. 15, No. 3, pp. 303-331, 2000 [10] K. M. Lee, X. Y. Hou, X. W. Ge, Y. Tang, “An analytical solution for a jointed shield-driven tunnel lining”, International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 25, No. 4, pp. 365-390, 2001 [11] Y. Koyama, “Present status and technology of shield tunneling method in Japan”, Tunnelling and Underground Space Technology, Vol. 18, No. 2, pp. 145-159, 2003 [12] F. Grübl, M. Thewes, “Recommendations of the German working group on design and calculation of shield machines”, Tunnelling and Underground Space Technology, Vol. 21, No. 3, pp. 261-261, 2006 Engineering, Technology & Applied Science Research Vol. 7, No. 2, 2017, 1504-1514 1514 www.etasr.com Rastbood et al.: Finite Element Based Response Surface Methodology to Optimize Segmental Tunnel… [13] L. Lu, X. Lu, P. Fan, “Full-ring experimental study of the lining structure of Shanghai Changjiang tunnel”, 4th International Conference on Earthquake Engineering, Taipei, Taiwan, 2006 [14] S. Teachavorasinskun, T. Chub-Uppakarn, “Experimental verification of joint effects on segmental tunnel lining”, Electronic Journal of Geotechnical Engineering, Vol. 14, pp. 1-8, 2008 [15] N. Do, D. Dias, P. Oreste, I. Djeran-Maigre, “2D numerical investigation of segmental tunnel lining behavior”, Tunnelling and Underground Space Technology, Vol. 37, pp. 115-127, 2013 [16] H. Morgan, “A contribution to the analysis of stress in a circular tunnel”, Geotechnique, Vol. 11, No. 1, pp. 37-46, 1961 [17] R. E. Ranken, J. Ghaboussi, A. Hendron Jr, Analysis of ground-liner interaction for tunnels, Department of Transportation, Office of the Secretary , USA, 1978 [18] C. W. Schwartz, H. H. Einstein, “Simplified Analysis For Ground- structurinet Eraction In Tunneling”, 21st US Symposium on Rock Mechanics (USRMS), American Rock Mechanics Association, 1980 [19] Japan Society of Civil Engineers, Japanese Standard for Shield Tunnelling, Tokyo,Japan, JSCE, 1996 [20] S. Liao, F. Peng, S. Shen, “Analysis of shearing effect on tunnel induced by load transfer along longitudinal direction”, Tunnelling and Underground Space Technology, Vol. 23, No. 4, pp. 421-430, 2008 [21] C. Carranza-Torres, B. Rysdahl, M. Kasim, “On the elastic analysis of a circular lined tunnel considering the delayed installation of the support”, International Journal of Rock Mechanics and Mining Sciences, Vol. 61, pp. 57-85, 2013 [22] N. Do, D. Dias, P. Oreste, “3D numerical investigation of mechanized twin tunnels in soft ground–Influence of lagging distance between two tunnel faces”, Engineering Structures, Vol.109, pp. 117-125, 2016 [23] N. Do, D. Dias, P. Oreste, I. Djeran-Maigre, “Behaviour of segmental tunnel linings under seismic loads studied with the hyperstatic reaction method”, Soil Dynamics and Earthquake Engineering, Vol. 79, Part A, pp. 108-117, 2015 [24] N. Do, D. Dias, P. Oreste, “3D numerical investigation on the interaction between mechanized twin tunnels in soft ground”, Environmental Earth Sciences, Vol. 73, No. 5, pp. 2101-2113, 2015 [25] N. Do, D. Dias, P. Oreste, I. Djeran-Maigre, “2D numerical investigation of segmental tunnel lining under seismic loading”, Soil Dynamics and Earthquake Engineering, Vol. 72, pp. 66-76, 2015 [26] ABAQUS/CAE User's Manual, Hibbitt, Karlsson & Sorensen, Inc. 2014 [27] H. Schulze, H. Duddeck, “Spannungen in schildvorgetriebenen Tunneln”, Beton-und Stahlbetonbau, Vol. 59, No. 8, pp. 169-175, 1964 [28] U. A. C. O. Engineering, Tunnels and Shafts in Rock Engineer Manual, EM 1110-2-2901. Dept. of the Army, Corps of Engineers, Washington, DC, 1997 [29] German Committee for Underground Construction, Recommendations for static analysis of shield tunnelling machines. Tunnel. 7, pp. 44-58, 2005 [30] K. Bakker, “Structural design of linings for bored tunnels in soft ground”, Heron, Vol. 48, No. 1, pp. 33-63, 2003 [31] T. Telford, Tunnel lining design guide, British Tunneling Society and the Institution of Civil Engineers. 2004 [32] X. Huang, H. Huang, J. Zhang, “Flattening of jointed shield-driven tunnel induced by longitudinal differential settlements”, Tunnelling and Underground Space Technology, Vol. 31, pp. 20-32, 2012 [33] F. Ye, C. Gou, H. Sun, Y. Liu, Y. Xia, Z. Zhou, “Model test study on effective ratio of segment transverse bending rigidity of shield tunnel”, Tunnelling and Underground Space Technology, Vol. 41, pp. 193-205, 2014 [34] F. Barpi, M. Barbero, D. Peila, “Numerical modelling of ground-tunnel support interaction using bedded-beam-spring model with fuzzy parameters”, Gospodarka Surowcami Mineralnymi, Vol. 27, pp. 71-87, 2011 [35] R. H. Myers, D. C. Montgomery, C. M. Anderson-Cook, Response surface methodology: process and product optimization using designed experiments, John Wiley & Sons, 2016 [36] D. C. Montgomery, Design and analysis of experiments. John Wiley & Sons, 2008