A procedure is presented for fatigue crack growth simulations of multiple cracks.Stress intensity factors are calculated using three-dimensional enriched elements.During simulations, models are re-meshed to reflect updated crack geometries.Results show excellent agreement when compared with data from the literature.The method can also handle multiple-cracks with high growth rate differences. A methodology is presented for fracture and fatigue crack growth simulations of multiple cracks. Computed stress intensity factors and a crack growth law are used to predict incremented positions of crack front nodes. The procedure is repeated until failure criterion is reached. The cases presented in this study include structures containing two or more cracks which propagate, coalesce or grow as planar and non-planar cracks. Comparisons of results with literature data shows excellent agreement. Therefore, it is concluded that the presented procedure can be used to accurately to assess fatigue crack propagation behavior of structures with multiple three-dimensional cracks.