Technical Brief

Efficient Fitting of the Feed Correction Polynomial for Real-Time Spline Interpolation

[+] Author and Article Information
Kaan Erkorkmaz

Associate Professor
Precision Controls Laboratory,
Department of Mechanical and Mechatronics Engineering,
University of Waterloo,
200 University Avenue West,
Waterloo, ON N2L 3G1, Canada
e-mail: kaane@uwaterloo.ca

Contributed by the Manufacturing Engineering Division of ASME for publication in the JOURNAL OF MANUFACTURING SCIENCE AND ENGINEERING. Manuscript received March 16, 2014; final manuscript received March 29, 2015; published online July 8, 2015. Assoc. Editor: Xiaoping Qian.

J. Manuf. Sci. Eng 137(4), 044501 (Aug 01, 2015) (8 pages) Paper No: MANU-14-1119; doi: 10.1115/1.4030300 History: Received March 16, 2014; Revised March 29, 2015; Online July 08, 2015

In spline toolpath interpolation, a crucial point is solving the mapping between the spline parameter (u) and actual arc length (s) accurately, so that the toolpath is traveled without undesirable fluctuations or discontinuities in the feedrate profile. To achieve this, various techniques have been proposed in literature, including Taylor series interpolation, iterative numerical methods, and approximating the mapping between u and s with a feed correction polynomial. This paper presents a new way to parameterize the seventh order feed correction polynomial, which was introduced by Erkorkmaz and Altintas (2005, “Quintic Spline Interpolation With Minimal Feed Fluctuation,” ASME J. Manuf. Sci. Eng., 127(2), pp. 339–349). The proposed technique has a closed-form solution that can be efficiently implemented in real-time, rather than having to construct and solve a linear equation system with 14 unknowns for each spline segment. In this paper, the new solution is derived step by step, and simulation case studies are presented which demonstrate that the new method accurately parameterizes the feed correction polynomial in approximately 43% less computational time, compared to applying the former solution of Erkorkmaz and Altintas (2005, “Quintic Spline Interpolation With Minimal Feed Fluctuation,” ASME J. Manuf. Sci. Eng., 127(2), pp. 339–349). This is because matrix multiplication operations and a dedicated linear equation solver, which are cumbersome to implement inside a real-time computer numerical controller (CNC), are avoided in the new solution.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Li, Z.-L., and Zhu, L.-M., 2014, “Envelope Surface Modeling and Tool Path Optimization for Five-Axis Flank Milling Considering Cutter Runout,” ASME J. Manuf. Sci. Eng., 136(4), p. 041021. [CrossRef]
Pan, Y., Zhou, C., Chen, Y., and Partanen, J., 2014, “Multitool and Multi-Axis Computer Numerically Controlled Accumulation for Fabricating Conformal Features on Curved Surfaces,” ASME J. Manuf. Sci. Eng., 136(3), p. 031007. [CrossRef]
Erkorkmaz, K., and Altintas, Y., 2005, “Quintic Spline Interpolation With Minimal Feed Fluctuation,” ASME J. Manuf. Sci. Eng., 127(2), pp. 339–349. [CrossRef]
Wang, F.-C., and Yang, D. C. H., 1993, “Nearly Arc-Length Parameterized Quintic-Spline Interpolation for Precision Machining,” Comput. Aided Des., 25(5), pp. 281–288. [CrossRef]
Wang, F.-C., Wright, P. K., Barsky, B. A., and Yang, D. C. H., 1999, “Approximately Arc-Length Parameterized C3 Quintic Interpolatory Splines,” ASME J. Mech. Des., 121(3), pp. 430–439. [CrossRef]
Farouki, R. T., and Shah, S., 1996, “Real-Time CNC Interpolators for Pythagorean-Hodograph Curves,” Comput. Aided Geom. Des., 13(7), pp. 583–600. [CrossRef]
Farouki, R. T., Al-Kandari, M., and Sakkalis, T., 2002, “Hermite Interpolation by Rotation-Invariant Spatial Pythagorean-Hodograph Curves,” Adv. Comput. Math., 17(4), pp. 369–383. [CrossRef]
Huang, J.-T., and Yang, D. C. H., 1992, “Precision Command Generation for Computer Controlled Machines,” Precision Machining: Technology and Machine Development and Improvement, Vol. PED 58, ASME, New York, pp. 89–104.
Shpitalni, M., Koren, Y., and Lo, C.-C., 1994, “Realtime Curve Interpolators,” Comput. Aided Des., 26(11), pp. 832–838. [CrossRef]
Otsuki, T., Kozai, H., and Wakinotani, Y., 1998, “Free-Form Curve Interpolation Method and Apparatus,” U.S. Patent No. 5,815,401.
Lin, R.-S., 2000, “Real-Time Surface Interpolator for 3-D Parametric Surface Machining on 3-Axis Machine Tools,” Int. J. Mach. Tools Manuf., 40(10), pp. 1513–1526. [CrossRef]
Heng, M., and Erkorkmaz, K., 2010, “Design of a NURBS Interpolator With Minimal Feed Fluctuation and Continuous Feed Modulation Capability,” Int. J. Mach. Tools Manuf., 50(3), pp. 281–293. [CrossRef]
Yuen, A., Zhang, K., and Altintas, Y., 2013, “Smooth Trajectory Generation for Five-Axis Machine Tools,” Int. J. Mach. Tools Manuf., 71, pp. 11–19. [CrossRef]
Lei, W. T., Sung, M. P., Lin, L. Y., and Huang, J. J., 2007, “Fast Real-Time NURBS Path Interpolation for CNC Machine Tools,” Int. J. Mach. Tools Manuf., 47(10), pp. 1530–1541. [CrossRef]
Cheng, C. W., Tsai, M. C., and Maciejowski, J., 2006, “Feedrate Control for Non-Uniform Rational B-Spline Motion Command Generation,” Proc. Inst. Mech. Eng., Part B, 220(11), pp. 1855–1861. [CrossRef]
Press, W. H., Flannery, B. P., Teukolsky, S. A., and Vetterling, W. T., 1992, Numerical Recipes in C: The Art of Scientific Computing, 2nd ed., Cambridge University Press, New York.
MathWorks, “MATLAB R2014b Documentation,” MathWorks, Inc., Natick, MA, http://www.mathworks.com/help/matlab/ref/mldivide.html
Altintas, Y., 2000, Manufacturing Automation: Metal Cutting Mechanics, Machine Tool Vibrations, and CNC Design, Cambridge University Press, Cambridge, UK.


Grahic Jump Location
Fig. 2

Three commonly used spline interpolation methods [3]: (a) natural interpolation, (b) interpolation using Taylor series expansion, and (c) interpolation with feedrate correction polynomial

Grahic Jump Location
Fig. 3

Parameterization of the feed correction polynomial using data generated from the numerical integration of the arc length (s) over the spline parameter (u) [3]

Grahic Jump Location
Fig. 4

Implementation comparing the new parameterization method for the feed correction polynomial with the earlier results reported in Ref. [3]. (a) 88-Segment toolpath, (b) natural interpolation, (c) first order Taylor, (d) feed correction polynomial (earlier solution, using matrix calculations [3]), and (e) feed correction polynomial (computed with the proposed solution).

Grahic Jump Location
Fig. 5

Comparison between seventh order feed correction polynomial and fifth order variations (i.e., justification of using seventh order): (a) seventh order feed correction polynomial parameterized considering u–s data, as well as zeroth, first, and second order boundary conditions (proposed), (b) fifth order feed correction polynomial parameterized considering u–s data and only zeroth and first order boundary conditions, and (c) fifth order feed correction polynomial parameterized considering only zeroth, first, and second order boundary conditions

Grahic Jump Location
Fig. 6

(a) Computational load and (b) average computational time comparison between the earlier solution [3] and proposed new method for computing the feed correction polynomial for one segment. The reported results are based on the 88-segment toolpath in Fig. 4(a), which was split, on average, to 440 divisions per segment (= Mk) during the computation of the segment arc-lengths and feed correction polynomials.

Grahic Jump Location
Fig. 7

Implementation of the proposed feed correction method on an airfoil profile. (a) 48-Segment toolpath, (b) natural interpolation, and (c) feed correction polynomial (computed with the proposed solution).



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In