Modern power systems are high-dimensional, strongly coupled nonlinear systems with complex and diverse dynamic characteristics. The polynomial model of the power system is a key focus in stability research. Therefore, this paper presents a study on the approximate transient stability solution targeting the fault process in power systems. Firstly, based on the inherent sinusoidal coupling characteristics of the power system swing equations, a generalized polynomial matrix description in perturbation form is presented using the Taylor expansion formula. Secondly, considering the staged characteristics of the stability process in power systems, the approximate solutions of the polynomial model during and after the fault are provided using coordinate transformation and regular perturbation techniques. Then, the structural characteristics of the approximate solutions are analyzed, revealing the mathematical basis of the stable motion patterns of the power grid, and a monotonicity rule of the system’s power angle oscillation amplitude is discovered. Finally, the effectiveness of the proposed methods and analyses is further validated by numerical examples of the IEEE 3-machine 9-bus system and IEEE 10-machine 39-bus system.