During strike-slip fault dislocation, multiple fault planes are commonly observed. The resulting permanent ground deformation can lead to profound structural damage to tunnels. However, existing analytical models do not consider multiple fault planes. Instead, they concentrate the entire fault displacement onto a single fault plane for analysis, thereby giving rise to notable errors in the calculated results. To address this issue, a refined nonlinear theoretical model was established to analyze the mechanical responses of the tunnels subjected to multiple strike-slip fault dislocations. The analytical model considers the number of fault planes, nonlinear soil‒tunnel interactions, geometric nonlinearity, and fault zone width, leading to a significant improvement in its range of applicability and calculation accuracy. The results of the analytical model are in agreement, both qualitatively and quantitatively, with the model test and numerical results. Then, based on the proposed theoretical model, a sensitivity analysis of parameters was conducted, focusing on the variables such as the number of fault planes, fault plane distance (d), fault displacement ratio (η), burial depth (C), crossing angle (β), tunnel diameter (D), fault zone width (Wf), and strike-slip fault displacement (Δfs). The results show that the peak shear force (Vmax), bending moment (Mmax), and axial force (Nmax) decrease with increasing d. The Vmax of the tunnel is found at the fault plane with the largest fault displacement. C, D, and Δfs contribute to the increases in Vmax, Mmax, and Nmax. Additionally, increasing the number of fault planes reduces Vmax and Mmax, whereas the variation in Nmax remains minimal.