Analysis of mechanical response and failure characteristics of tunnels crossing active faults considering axial force and geometrical nonlinearity

Tunnels subjected to active fault dislocation may experience significant damage. This paper establishes a novel methodology and solution procedure for analyzing the mechanical response and failure characteristics of tunnels subjected to active fault dislocation based on an elastic foundation beam model. The proposed methodology includes axial, transverse, and vertical soil-tunnel interaction terms, in addition to geometrical nonlinearity and axial force terms in the governing equation. This approach has a significantly extended application range, effectively addressing the problems encountered in tunnels crossing active faults with diverse crossing angles, dip angles, and fault types. The proposed methodology is verified by comparison to a 3D FEM model with various fault types, experimental tests, and on-site case, and the results are in excellent quantitative and qualitative agreement with the numerical, experimental, and on-site results. When fault displacement is below 0.5 m, disregarding geometric nonlinearity results in calculation errors of approximately 10% to 17% for peak axial force (Nmax), 18% to 22% for peak shear force (Vmax), and 20% to 30% for peak bending moment (Mmax). Finally, the responses caused by different factors, i.e., fault type, fault displacement, tunnel stiffness, and tunnel diameter, are investigated in detail to better understand tunnels crossing active faults. The results show that amongst the various fault types, the Vmax and the Mmax experienced by tunnels subjected to oblique slip-fault dislocation surpass those of other fault types, accompanied by the most extensive failure range. The augmentation of fault displacement, tunnel stiffness, and tunnel diameter precipitates a corresponding escalation in the Nmax, Vmax, Mmax, and failure range. Under oblique-slip fault dislocation, the tunnel undergoes an initial phase of shear failure, followed by tension-bending failure, delineated by distinct fault displacement thresholds of 0.1 m and 0.2 m, respectively. The proposed methodology provides the advantage of reliable stability analysis and design of tunnels crossing active faults. •Geometric nonlinearity and axial force are applied into the model.•The method is applicable to various fault types, crossing angles, and dip angles.