Modeling and simulation of trajectory smoothing and feedrate scheduling for vibration-damping CNC machining
Introduction
The mathematical models for tool path smoothing and feedrate scheduling are critical for high-precision and low-vibration machining. The tool path generated by computer-aided manufacturing (CAM) system is the piecewise linear trajectory with only G° continuity, leading to large machinery vibration at the sharp corners. To decline the impact of machine tool during the process and raise the surface quality, Gk (k ≥ 2) trajectory smoothing is indispensable [1]. Recently, the scholars [2,3] found the tool path with higher geometric continuity can raise the feedrate at the sharp corner, declining the fluctuation of feedrate and acceleration, raising the machining quality. On the other hand, the feedrate profile with high smoothness is a necessary condition to guarantee the continuity of high-order axial kinematic variables (such as acceleration, jerk, snap), increasing the smoothness of feedrate trajectory and reduces the impact caused by axis actuators [4], [5], [6]. Hence it is significant to develop a mathematical model of tool-path smoothing and feedrate scheduling for damping the machinery vibration and raising the surface quality.
The trajectory smoothing methods includes approximation schemes [7], [8], [9], [10], [11], [12], [13], [14], transition schemes [1,[15], [16], [17], [18], [19], [20]], and interpolation schemes [21], [22], [23], [24], [25], [26], [27], [28], [29], [30], [31], [32], [33]. These works are valuable to improve the machining quality and efficiency. However, there still exist some shortcomings when these methods were employed in CNC machining. The comparison of Gk (k ≥ 2) tool-path smoothing schemes in recent studies is provided in Table 1. Besides the deficiencies in controlling the chord error, interpolating G01 points, the form of curvature extremum, and real-time performance, the existing tool-path smoothing schemes performs insufficiently in reducing the curvature extrema and raising the smoothness of tool path.
To generate Gk (k ≥ 2) tool path and compress the number of the data points, the approximation schemes using B-spline or Non-uniform Rational B-spline (NURBS) curve have been widely studied in the literatures [7], [8], [9], [10], [11], [12], [13], [14], including the Least Square Fitting (LSF) methods [7], [8], [9], [10] and the progressive and iterative approximation (PIA) methods [11], [12], [13], [14]. Park et al. [7,8] and Li et al. [9] adopted the conventional LSF methods to generate the smooth trajectory with the confined error between data points and the smoothed trajectory. Nevertheless, the chord error was unconstrained, and both computational efficiency and numerical stability were low in these works. Yeh et al. [10] developed a NURBS fitting method with LSF to generate G2 trajectory in real time. The error between G01 points and the smoothed path was constrained, but the chord error was not constrained. Lin et al. [11] proposed an extended PIA method to reduce the number of control points. Deng et al. [12] further presented a progressive and iterative approximation method based on least squares fitting (LSPIA). Nonetheless, the chord error constraint was not considered in the works [11,12]. He et al. [13] incorporated the energy objective function into the LSPIA method, and proposed ELSPIA method to generate the smooth trajectory with confined chord error. Nonetheless, the algorithm was time-consuming due to the iterative least square approximation, and the chord-error conformance possibly weakens the overall smoothness of tool path. Bi et al. [14] proposed a fast B-spline fitting scheme with PIA method to generate G2 tool path under confined chord error. For the approximation schemes [7], [8], [9], [10], [11], [12], [13], [14], the general deficiency is that the smoothed path has undefined analytical curvature extrema, leading to the heavy computational load in the stage of feedrate scheduling, hampering the real-time performance of tool path processing.
To conveniently control the shape of tool path, Gk (k ≥ 2) transition schemes [1,[15], [16], [17], [18], [19], [20]] are performed using the transitional curve, such as B-spline curve, Bezier curve. The authors [1,[15], [16], [17], [18]] performed G2 transition scheme at the sharp corner. Nevertheless, the chord error of G2 trajectory was not constrained in Ref. [15,16], and the analytical curvature extremum of G2 trajectory was unavailable in Ref. [16]. Zhao [1] and Farouki [17] realized G2 corner smoothing under the confined chord error in real time, and the tool path generated had the analytical curvature extrema. To generate smoother tool path and decline the jerk during the process, G3 transitional trajectory has been recently designed by the works [18], [19], [20]. Piazzi [18] and Fang [20] constructed G3 trajectory between the splines, but the trajectory had only G1 continuity at the joints. Besides, the approximation error was unconstrained in the works [18,19]. Although the transition schemes generated the shape-preserving trajectories, the trajectory did not interpolate G01 points, leading to the loss of the featured position of original contour.
To preserve the information of featured positions of original contour, Gk (k ≥ 2) interpolation methods have been extensively studied in the works [2,[21], [22], [23], [24], [25], [26], [27], [28], [29], [30], [31], [32], [33]]. Compare with approximation schemes and transition schemes, interpolation methods can guarantee the smoothed trajectory passes through data points, and avoid the error between points and the smoothed trajectory. In those works [21], [22], [23], [24], [25], [26], [27], [28], [29], [30], [31], [32], [33], G2 interpolative trajectories have been generated. Nevertheless, in these works [21], [22], [23], [24], [25], [26], [27], [28], the chord error was not constrained rigorously. The methods [25], [26], [27], [28], [29] were only applicable to the planar data or monotone data. Fan et al. [30,31] adopted Bezier curve to generate interpolative trajectory under confined chord error but the trajectories generated were not rigidly shape-preserving. G3 interpolation scheme has been recently studied by the works [2,32,33]. Bajaj [32] and Renka [33] generated a G3 shape-preserving interpolative trajectory for spatial point sequence but the chord error was not constrained and analytical curvature extremum was unavailable. Besides, in this work, a nonlinear optimization problem needed to be solved, and thereby the algorithm was only performed offline. Fan et al. [2] designed a G3 interpolative path for G01 trajectory. However, there still exists the space of declining curvature extrema and raising the smoothness at the sharp corners.
After generating the smooth trajectory, the geometric information of trajectory is transferred to the feedrate scheduling module in CNC system. The feedrate profile under confined jerk decreases greatly the excitation of the resonance frequencies of machine tool, declining machinery vibrations. In the works [34], [35], [36], [37], [38], jerk-limited feedrate scheduling schemes have been developed. Nonetheless the methods are only performed offline. Tajima et al. [39] generated jerk-limited trajectory under confined approximation error for corner smoothing in real time but the trajectory was not interpolative. Lai [40] and Annoni [41] performed real-time feedrate scheduling under confined jerk and fine interpolation error for NURBS trajectory. The kinematic smoothness helps to decline the flexible impact created by the axis actuator. To raise the kinematic smoothness, Fan [6], Liu [42], and Bharathi [43] realized time-optimal and jerk-continuous feedrate scheduling but these methods were performed offline. Jin et al. [44] proposed a real-time and jerk-continuous feedrate scheduling scheme but 7-phase ACC/DEC profile adopted was not time-optimal. In addition to the deficiencies in real-time performance and time-optimality, the kinematic smoothness in the existing feedrate scheduling schemes is insufficient, and thus there still exists the improving space for damping the impact caused by axis actuator.
High-smoothness tool path and feedrate profile are indispensable to generate smooth kinematic profiles, and helpful to damp the impact caused by axis actuator. Aiming to reduce the vibration and guarantee high machining efficiency, this work considers the various requirements of CNC machining to develop a real-time methodology, including G4-continuous tool-path interpolation scheme and time-optimal feedrate scheduling scheme with smooth jerk. The remainder of the work is arranged as follows. Section 2 proposed G4 trajectory model and provides the generation algorithm of interpolative trajectory. Section 3 proposes jerk-smooth feedrate mode and the look-ahead algorithm of feedrate planning. In Section 4, the advantages of G4 tool path and jerk-smooth feedrate mode are shown. Section 5 provides the simulation and comparison. The experiment and analysis are done in Section 6. The conclusion is presented in Section 7.
Section snippets
Geometric model of G4 continuous transitional curve
The geometric model of G4 continuous transitional curve connecting adjacent line segments is constructed in this section. As illustrated in Fig. 1, the transitional curve is composed by a symmetric nine-degree Bezier curve with the following form
The transitional curve has G4 continuity at the joint points A0 and A9, and has G∞ continuity within the curve. Thus the whole transitional curve has G4 continuity. Assume , ,, ,
Jerk-smooth feedrate scheduling for G4 tool path
After generating G4 interpolative path in the geometric module of CNC system, the geometric information of the smoothed tool path is transferred into the feedrate planning module. In this section, jerk-smooth feedrate model is proposed, and time-optimal feedrate profiles under different given conditions are computed. Finally, jerk-smooth feedrate scheduling scheme with look-ahead function is performed.
The role of G4 tool path in raising smoothness and reducing jerk
In this section, the performance of G4 tool path in raising smoothness and reducing jerk are analyzed and compared with the methods [1,2] , as illustrated in Fig. 6.
Given the same approximation error δ, the maximal curvatures of transitional curves for G2 transition [1] and G3 transition schemes [2] are respectively
The curvature of G4 transitional curve reaches the maximal value at the midpoint Q, expressed with
Simulation
In this section, the effectiveness and performance of the algorithm proposed are illustrated by several examples. Compared with the previous methods, the advantages of the proposed methods are shown in this section.
Experiment
In this section, the experiments are performed to compare the vibration intensity and machining quality among G2 transition scheme with jerk-limited feedrate scheduling [1], G3 interpolation scheme with jerk-continuous feedrate scheduling [2] and the proposed. Additionally, the experimental analysis about maximal crackle and snap is performed, and the rule of determining appropriate values of maximal crackle and snap are provided in the industrial application.
Conclusion
In recent studies on trajectory smoothing and feedrate scheduling schemes, there are still some requirements that cannot be satisfied simultaneously, including confined chord error, G01 points interpolated, analytical curvature extremum, real-time performance, kinematic time-optimality, and the smoothness of tool-path and feedrate profile. To reduce the vibration and guarantee high machining efficiency, this work develops G4 tool-path interpolation scheme under various geometric constraints,
Acknowledgements
This research is supported by Science Challenge Project (TZ2016006-0502-03) and National Natural Science Foundation of China(No.61605182). The authors gratefully acknowledge Dr. Xingzhan Li and technician Ke Wang from Institute of Mechanical Manufacturing Technology of China Academy Engineering Physics, who provided the assistance for the processing validation and vibration test.
References (49)
- et al.
A real-time look-ahead interpolation methodology with curvature-continuous B-spline transition scheme for CNC machining of short line segments
Int. J. Mach. Tools Manuf
(2013) - et al.
A realtime curvature-smooth interpolation scheme and motion planning for CNC machining of short line segments
Int. J. Mach. Tools Manuf.
(2015) Dissociated jerk-limited trajectory applied to time-varying vibration reduction
Robot. Comput. Integr. Manuf.
(2013)- et al.
New damped-jerk trajectory for vibration reduction
Control Eng. Pract.
(2014) An error-bounded approximate method for representing planar curves in B-splines
Comput. Aided Geom. Des.
(2004)- et al.
B-spline curve fitting based on adaptive curve refinement using dominant points
Comput.-Aided Des.
(2007) - et al.
Adaptive knot placement in B-spline curve approximation
Comput.-Aided Des.
(2005) - et al.
An extended iterative format for the progressive-iteration approximation
Comput. Graph.
(2011) - et al.
Progressive and iterative approximation for least squares B-spline curve and surface fitting
Comput.-Aided Des.
(2014) - et al.
A chord error conforming tool path B-spline fitting method for NC machining based on energy minimization and LSPIA
J. Comput. Des. Eng.
(2015)
G3 approximation of conic sections by quantic polynomial curves
Comput. Aided Geom. Des.
Osculatory interpolation
Comput. Aided Geom. Des.
Bezier curves and C2 interpolation in Riemannian manifolds
J. Approx. Theory
Geometric Hermite interpolation with circular precision
Comput. Aided Geom. Des.
G2 Hermite interpolation with circular precision
Comput. Aided Des.
C-shaped G2 Hermite interpolation with circular precision based on cubic PH curve interpolation
Comput.-Aided Des.
Monotonicity-preserving C2 rational cubic spline for monotone data
Appl. Math. Comput.
A-splines: local interpolation and approximation using Gk-continuous piecewise real algebraic curves
Comput. Aided Geom. Des.
Shape-preserving interpolation by fair discrete G3 space curves
Comput. Aided Geom. Des.
Feedrate optimization with jerk constraints for generating minimum-time trajectories
Int. J. Mach. Tools Manuf.
A smooth curve evolution approach to the feedrate planning on five-axis tool path with geometric and kinematic constraints
Int. J. Mach. Tools Manuf.
Kinematic corner smoothing for high speed machine tools
Int. J. Mach. Tools Manuf.
A real-time configurable NURBS interpolator with bounded acceleration, jerk and chord error
Comput.-Aided Des.
Time-optimal and jerk-continuous trajectory planning for robot manipulators with kinematic constraints
Robot. Comput.-Integr. Manuf.
Cited by (22)
Near-time optimal feedrate planning for the NURBS curve considering interpolation error constraints
2024, Robotics and Computer-Integrated ManufacturingAn integrated trajectory smoothing method for lines and arcs mixed toolpath based on motion overlapping strategy
2023, Journal of Manufacturing Processes2D NURBS Curve discretization method with double constraints of area and chord error and contour accuracy prediction modeling
2023, Simulation Modelling Practice and TheoryCitation Excerpt :Dong et al. [9] obtained a smooth feed speed based on the curvature change. Fan et al. [10] considered both trajectory smoothing and velocity scheduling problems, smoothed the corners of the G code trajectory with nine times Bezier curve, and performed velocity planning by Jerk-smooth feedrate method, and his work obtained a smooth machined surface. However, there is already an error in the short segment trajectory when born from the original curve.
An analytical method for corner smoothing of five-axis linear paths using the conformal geometric algebra
2022, CAD Computer Aided DesignGlobal smoothing of short line segment toolpaths by control-point-assigning-based geometric smoothing and FIR filtering-based motion smoothing
2021, Mechanical Systems and Signal ProcessingCitation Excerpt :Fan et al. [5] accomplished G3 continuity by approximating one corner with two quartic Bezier splines whose curvature energy was minimalized, and a 15 phases acceleration/deceleration profile was employed to realize the jerk smooth interpolation. To further improve the smoothness and reduce vibration, Fan et al. [6] and Zhang et al. [7] approximated the sharp corners with a nine-degree Bezier curve and a quantic B-spline curve, respectively, thus realizing the G4 continuous. In order to combine the geometric smoothing and motion smoothing into one step, Tajima et al. [8] and Sencer et al. [9,10] proposed a kinematic cornering strategy based on FIR (Finite Impulse Response) filtering.