**4. Conclusions**

The results of the developed program simulation were compared with experimental and numerical data for the two internal non-colinear cracks and the three-point bending beam with three holes with two different configurations. The developed program combines the adaptive mesh refinement with increasing mesh density in the required area only in order to reduce the computational time while increasing the solution accuracy. The norm stress error is taken as a posterior estimator for the h-type adaptive refinement. With this series of simulations, the capability of the developed program was demonstrated to accurately predict the crack path trajectory, stress intensity factors, and fatigue life under constant amplitude loading. In these simulation sequences, holes act as a crack stopper and attract a crack trajectory to growth. Such findings support that the algorithm can be used to identify crack-stopping holes used in damage tolerance designs.

**Author Contributions:** Conceptualization, A.M.A.; methodology, A.M.A.; software, A.M.A.; validation, A.M.A., and Y.A.F.; formal analysis, A.M.A.; investigation, A.M.A. and Y.A.F.; resources, A.M.A.; data curation, A.M.A. and Y.A.F.; writing—original draft preparation, A.M.A.; writing—review and editing, A.M.A. and Y.A.F.; visualization, A.M.A. and Y.A.F.; supervision, A.M.A.; project administration, A.M.A. and Y.A.F.; funding acquisition, A.M.A. and Y.A.F. Both authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Data Availability Statement:** The data presented in this study are available in the article.

**Conflicts of Interest:** The authors declare no conflict of interest.
