Nonlinear vibrations and their control are critical in improving the magnetic bearings system performance and in the more widely spread use of magnetic bearings system. Multiple objective genetic algorithms (MOGAs) simultaneously optimize a vibration control law and geometrical features of a set of nonlinear magnetic bearings supporting a generic flexible, spinning shaft. The objectives include minimization of the actuator mass, minimization of the power loss, and maximization of the external static load capacity of the rotor. Levitation of the spinning rotor and the nonlinear vibration amplitude by rotor unbalance are constraint conditions according to International Organization for Standardization (ISO) specified standards for the control law search. The finite element method (FEM) was applied to determine the temperature distribution and identify the hot spot of the actuator during steady-state operation. Nonlinearities include magnetic flux saturation, and current and voltage limits of power amplifiers. Pareto frontiers were applied to identify and visualize the best-compromised solutions, which give a most compact design with minimum power loss whose vibration amplitudes satisfy ISO standards.