Fracture mechanics analysis is nowadays adopted in several industrial fields to assess the capability of components to withstand fatigue loads. Finite Element Method (FEM) is a well-established tool for the evaluation of flaw Stress Intensity Factors (SIF) and for the survey of its propagation. Nevertheless the study of the growth of near-surface circular and elliptical cracks is still an arduous task to be faced with FEM. In fact, the interaction of the flaw with free surfaces leads the crack front to assume complex shapes, whose simulation cannot be easily accomplished. A possible answer to deal with such a problem is to use the mesh morphing technique, a nodal relocation methodology, that allows to cover different problems. In fact, with mesh morphing, it is possible to fit the baseline flaw front with the desired shape (generic shape) and to automatically simulate its evolution at a certain number of cycles. In the proposed work this approach is demonstrated exploiting ANSYS Mechanical as FEM tool and RBF Morph ACT Extension as mesh-morpher. The results of the proposed workflow are compared with those available in literature.