This study presents a developed finite element code written by Visual Fortran to computationally model fatigue crack growth (FCG) in arbitrary 2D structures with constant amplitude loading, using the linear elastic fracture mechanics (LEFM) concept. Accordingly, optimizing an FCG analysis, it is necessary to describe all the characteristics of the 2D model of the cracked component, including loads, support conditions, and material characteristics. The advancing front method has been used to generate the finite element mesh. The equivalent stress intensity factor was used as the onset criteria of crack propagation, since it is the main significant parameter that must be precisely predicted. As such, a criterion premised on direction (maximum circumferential stress theory) was implemented. After pre-processing, the analysis continues with incremental analysis of the crack growth, which is discretized into short straight segments. The adaptive mesh finite element method was used to perform the stress analysis for each increment. The displacement extrapolation technique was employed at each crack extension increment to compute the SIFs, which are then assessed by the maximum circumferential stress theory to determine the direction of the crack growth and predict the fatigue life as a function of crack length using a modified form of Paris’ law. The application examples demonstrate the developed program’s capability and performance.