The hp-version Legendre spectral collocation method for solving nonlinear fractional differential equations was studied. At first, the multifractional differential equation was transformed into an equivalent Volterra integral equation. Then a numerical method to approximate the original equation was constructed. Finally, the correctness of the algorithm theory and the effectiveness of the proposed numerical method were demonstrated by numerical experiments.