We develop a program for predicting bearing performance, including the maximum contact stress, bearing power loss, minimum lubricating film thickness, static safety factor, and dynamic life, and an optimal design for high-speed angular contact ball bearings in aircraft gearboxes using optimization algorithms. A program is developed to predict bearing performance by calculating the loads and displacements of individual balls based on the bearing’s kinematic and static equilibrium equations. Next, an optimization program is developed with non-dominated sorting genetic algorithm III (NSGA-III), wherein calculated bearing performance indicators serve as constraints or objective functions. The performance of individuals within the final Pareto front is evaluated using the inter-criteria correlation (CRITIC) and weighted product methods. The optimization performances of NSGA-III and NSGA-II are compared based on their hypervolume indicators. A solution of the kinematic and static equilibrium equations under the load cases and duty cycles of real aircraft gearboxes is obtained, enhancing the prediction accuracy. Furthermore, optimization, including bearing performance metrics rarely considered in previous studies, is performed. Optimal specifications superior to reference bearings employed in aircraft gearboxes are obtained, enhancing all three objective functions or specific objective functions. The optimization performance of NSGA-III is confirmed to surpass that of conventional NSGA-II.