Three-dimensional enriched finite elements are used to compute mixed-mode stress intensity factors (SIFs) for three-dimensional cracks in elastic functionally graded materials (FGMs) that are subject to general mixed-mode loading. The method, which advantageously does not require special mesh configuration/modifications and post-processing of finite element results, is an enhancement of previou...