In this paper, we have investigated the problem of the convergence rate of the multiple integral
{∫T0···∫T0 f(t1, ··· , tn)dBH, πt1 ···dBH, πtn},
where f ∈ Cn+1([0, T]n) is a given function, π is a partition of the interval [0, T] and {BH, πti } is a family of interpolation approximation of fractional Brownian motion BHt with Hurst parameter H < 1/2. The limit process is the multiple Stratonovich integral of the function f. In view of known results, the convergence rate is different for different multiplicity n. Under some mild conditions, we obtain that the uniform convergence rate is Δ2H in the mean square sense, where Δ is the norm of the partition generating the approximations.