The dynamic characteristics of multiple combinations of angular contact ball bearings have significant effects on the bearing-rotor system. Based on the dynamics analysis and raceway control theory of the rolling bearing, the nonlinear equations were solved using Powell optimization and the Newton- Raphson algorithm, and a corresponding program was proposed to calculate the stiffness of multiple combinations of angular contact ball bearings. Taking the DBD-7016A5 bearings as an example, the results show that the axial preload value and rotational speed were nonlinear to both the stiffness and displacement of each single bearing or the combined multi-row bearing. The combined axial stiffness of the whole bearing was greater than that of a single one, but the radial stiffness was less. In order to fully realize the performance of the bearings, the total axial displacements of the multi-row bearing must be less than the clearance of the inner ring of unilateral and bilateral bearings.