Abstract: In the northwest, a lot of multi-stage high slopes with loess have been encountered in the construction of basic engineering. Slope engineering belongs to permanent project. So it is unavoidable to meet the conditions of rainfall. Under rainfall infiltration, the degree of saturation and mass of soil slope increase, which can further cause the change of slope matric suction and effective stress. As a result, the stability of the slope can be greatly influenced. In this paper, the coupling calculation model of seepage field and stress field is established. It considers the changes of the matrix suction as the rain water infiltration at the different depths of the slope body. Then, the stability of the multi-stage high slope with loess under the condition of rainfall infiltration can be analyzed. Then, according to an actual engineering and using the PLAXIS 3D geotechnical finite element software, the calculation model of slope stability is established. By setting the variation function of rainfall with the time and boundary condition, the stability calculation of multi-stage high slope with loess is carried out under the condition of rainfall infiltration. And then, The numerical calculation results are compared with the theoretical calculation results. According to the numerical calculation results, the changes of the deformation, matrix suction, effective stress, potential slip surface and safety factor of the multi-stage high slope with loess under the condition of rainfall infiltration are obtained. And the stability of the multi-stage high slope with loess under the condition of rainfall infiltration are analyzed. The study results can provide some guidance for the design of the multi-stage high slope with loess under the rainfall infiltration.