<p>This study presents a comprehensive analytical solution, under plane strain conditions, to evaluate the evolution of shear and normal stresses along a fault following a circular tunnel excavation. The formulation incorporates the effects of overburden pressure and deconfinement rate. Key parameters investigated include the deconfinement rate, tunnel depth, fault angle, tunnel size, fault-tunnel distance, fault pore pressure, and the initial stress state of the surrounding rock. These factors play a significant role in the redistribution of stresses and the potential for fault slip. The proposed analytical solution is first validated against a numerical finite volume model. The evolution of the fault factor of safety is analyzed for each parameter and compared with numerical results to evaluate the combined influence of all investigated parameters. Closed-form equations using dimensionless parameters are developed to estimate both the fault factor of safety and slip length, showing strong agreement with the analytical data. In addition, a case study is presented to validate the proposed approach, demonstrating its effectiveness in predicting the safety factor, slip length, and slip locations.</p>