Shafts are among the most common machine elements. The typical shape used to reduce stress concentrations is circulararches due to the simplicity. A shaft is typically loaded by axial, bending, and torsional loads in different combinations.The stress concentration factors are found in tables and charts. The circular design is not optimal from a strength pointof view, and the strength can be increased using shape optimization. It is in this article shown how the maximum stressfrom the combined loads can be minimized, when the shape is parameterized using the simple super ellipse. This makesthe resulting optimized designs easily transferable to practical design. The stress evaluation is numerically performedusing the finite element method using harmonic elements that facilitates an axisymmetric model although the loading isunsymmetric. The stresses are reduced by up to a factor of 15%.