This paper addresses the numerical simulation of fatigue crack growth in arbitrary 2D geometries under constant amplitude loading by the using a new finite element software. The purpose of this software is on the determination of 2D crack paths and surfaces as well as on the evaluation of components Lifetimes as a part of the damage tolerant assessment. Throughout the simulation of fatigue crack propagation an automatic adaptive mesh is carried out in the vicinity of the crack front nodes and in the elements which represent the higher stresses distribution. The fatigue crack direction and the corresponding stress-intensity factors are estimated at each small crack increment by employing the displacement extrapolation technique under facilitation of singular crack tip elements. The propagation is modeled by successive linear extensions, which are determined by the stress intensity factors under linear elastic fracture mechanics (LEFM) assumption. The stress intensity factors range history must be recorded along the small crack increments. Upon completion of the stress intensity factors range history recording, fatigue crack propagation life of the examined specimen is predicted. A consistent transfer algorithm and a crack relaxation method are proposed and implemented for this purpose. Verification of the predicted fatigue life is validated with relevant experimental data and numerical results obtained by other researchers. The comparisons show that the program is capable of demonstrating the fatigue life prediction results as well as the fatigue crack path satisfactorily.