The nonlinear global behavior evolution process of the non-orthogonal face gear-bearing system was studied in this paper. Firstly, a multi-degree-of-freedom coupled vibration model of the non-orthogonal face gear-bearing system was established, which took into account factors such as time-varying meshing stiffness, tooth backlash, bearing clearance, and transmission error. Then the cell mapping method was used to calculate the phase trajectory of each cell on the Poincaré section to 1000 steps, and observed the influence of the initial state of the system on the dynamic behavior of the system after a long time history. This paper used the global bifurcation diagram, the global behavior diagram obtained by the cell mapping method, the multi-initial value phase diagram and the Poincaré section to illustrate the evolution process of the global behavior of the system with the excitation frequency, meshing damping and load fluctuation. The research results show that with the change of control parameters, the system gradually changed from periodic response to chaos after a long time history, and presented complex phenomena such as the coexistence of multiple periodic responses and the coexistence of periodic responses and chaotic responses in the evolution process.