ABSTRACT. Shape optimization of components subjected to multiaxial fatigue is studied. The application of static optimization algorithms, generally based on Von-Mises stress, is no more applicable in these cases. In this paper an optimizing routine for multiaxial fatigue is developed on the basis of the CAO technique proposed by Mattheck. According to a criterion for fatigue strength estimation of notched specimens made of ductile materials and subjected to mutiaxial fatigue: Liu-Zenner. Abaqus 6.8-1is used as the commercial software to develop finite element simulation and Python 2.4 and Matlab2007 are used as the subroutine programs