A numerical procedure is described for extracting mixed-mode stress intensity factors along the fronts of three-dimensional, nonplanar cracks embedded in solids. The mixed-mode stress intensity factors at points along the crack front are obtained by evaluating interaction energy integrals for three-dimensional, non-planar cracks. To assess the validity of the numerical procedure, two numerical examples are considered. First, we consider the problem of a non-planar, lens-shaped crack in an infinite solid subjected to hydrostatic tension. The numerical results are shown to be in excellent agreement with available analytical results. We then consider the case of a non-planar, warped elliptical crack surface, where to our knowledge no analytical solution exists, and the results are discussed.