Traditional interpolation algorithms often fail to meet the precision requirements of ultra-precision machining when applied to super-precision polishing machines. Moreover, the processing of fused silica glass optical components is frequently threatened by mechanical impacts due to their fragility. In order to address this issue, this paper proposes a high-precision motion interpolation method based on fourth-order Runge-Kutta and velocity-flexible planning. This method is designed for open-architecture small multi-axis optical polishing machines to polish quartz glass. The algorithm initially employs composite Simpson's rule to calculate the lengths of sub-paths within the polishing trajectory. Based on these length values, flexible velocity planning is executed to ensure the smooth continuity of velocity, acceleration, and jerk during motion interpolation. This reduces the risk of mechanical impacts that could damage the components during the machining process. The introduction of the adaptive fourth-order Runge-Kutta method significantly enhances the parameter point calculation accuracy of NURBS curves. The incorporation of adaptive principles also maintains a higher processing speed, thereby greatly improving processing efficiency. This method comprehensively addresses both the precision of curve interpolation and execution efficiency. Finally, experimental validation is conducted on an open-architecture small multi-axis optical polishing machine. The proposed method based on FRKVF not only mitigates mechanical impacts resulting from discontinuous acceleration, thereby ensuring the machining quality of optical components, but also satisfies the high-precision requirements for processing fused silica optical elements.