Next Article in Journal
An Object Detection Method Based on Feature Uncertainty Domain Adaptation for Autonomous Driving
Previous Article in Journal
Particle-Velocity-Based Mixed-Source Sound Field Translation for Binaural Reproduction
Previous Article in Special Issue
Effect of Robot Motion Accuracy on Surface Form during Computer-Controlled Optical Surfacing Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of Pitch Cycle-Based Iterative Learning Contour Control for Thread Milling Operations in CNC Machine Tools

1
Department of Mechanical Engineering, National Taipei University of Technology, Taipei 10608, Taiwan
2
Institute of Mechatronic Engineering, National Taipei University of Technology, Taipei 10608, Taiwan
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(11), 6447; https://doi.org/10.3390/app13116447
Submission received: 21 April 2023 / Revised: 22 May 2023 / Accepted: 23 May 2023 / Published: 25 May 2023
(This article belongs to the Special Issue Advances in Intelligent Robots and Precision Machining)

Abstract

:
The helical contour motion accuracy of feed drive axes is important for thread milling operations in computer numerical control (CNC) machine tools. However, the motion dynamics and external disturbances significantly affect the contour motion results, while the feed drive axes perform helical motions in thread milling operations. Although existing iterative learning contour control (ILCC) methods can improve contour motion accuracy, the problems of data recording and processing on memory usage and computational burden in control systems, wasted materials, and increased costs in thread manufacturing still limit the practical applications of ILCC. Therefore, considering the similar motion dynamics and external disturbances of the feed drive axes during the pitch cycle motions of a helical path, this study developed a pitch cycle-based iterative learning contour control (PCB-ILCC) method to address the control system and thread manufacturing problems caused by the use of ILCC. For PCB-ILCC, this study adopted contour error vector estimation by referring to the interpolated positions on the pitch cycle of the helical path to simplify the computational complexity and designed the ILCC using the cycle learning method to easily implement the ILCC structure. Thus, this study developed a permanent magnet synchronous motor (PMSM) driving control utilizing the robust control method to mitigate the problems of motion dynamics and external disturbances on the feed drive axes. Thread milling experiments performed on a five-axis CNC machining center demonstrated the feasibility of the PCB-ILCC and validated that it can significantly improve the helical contour motion accuracy of the feed drive axes and achieve an 80% contour error reduction rate in comparison with the proportional–proportional–integral control, which is extensively used in commercialized PMSM drivers and CNC controllers.

1. Introduction

A thread milling operation is a machining process in which the feed drive axes of a machine tool simultaneously move a thread milling tool along a helical path to perform a cutting operation. Therefore, the helical contour motion accuracy of feed drive axes while performing thread milling operations becomes crucial. In recent decades, owing to the importance of the contour motion accuracy of computer numerical control (CNC) machine tools, many contour control designs have been proposed to improve the contour motions of feed drive axes. Liu et al. [1] proposed a velocity-based contour control method that used a B-spline-based wavelet neural network to approximate unknown dynamics and used neural network identification to design a robust control to improve the contour motion of dual-linear-motor-driven gantry stages. Chen et al. [2] transformed the desired contour path and tool orientation in the workpiece coordinate frame to a five-dimensional hypercurve in a machine coordinate frame and then adopted the equivalent error method to design a contouring controller to reduce the errors of the contour path and tool orientation. Zhao et al. [3] compensated weighted contour error components to velocity commands in the robotic task space of a six-degrees-of-freedom robot arm. They proposed a component-based contour control mechanism that can improve the contour motion accuracy of the robot arm performing machining tasks. Liu et al. [4] integrated generalized predictive control and feedback correction algorithms to develop a real-time contour error control design for five-axis machine tools. Bo et al. [5] improved the conventional cross-coupled control design to increase the contour accuracy of biaxial motion control systems. Li et al. [6] adopted a neural memory network to develop a contour error estimation and compensation method for CNC systems. Song et al. [7] designed a five-axis cross-coupling control method in joint space to simultaneously reduce task-space and joint-space contour errors. Wang et al. [8] referred to the calculations of distance and contour errors and adopted the trajectory compensation method with gain adjustment to modify the reference trajectory for improving the contour control performance of motion control systems. Zhang et al. [9] considered the input and output data of a two-dimensional linear motor system while performing reciprocating motions and then proposed a model-free adaptive control algorithm that can improve the contour motion accuracy of the linear motor system. Kuang et al. [10] proposed a discrete-time fractional-order sliding-mode contour error control design based on fractional-order sliding surface stability analysis and linear matrix inequality. Ma et al. [11] proposed an equivalent-plane cross-coupling control design to apply the developed biaxial cross-coupling control to three-axis contour error control. Shi and Lou [12] transformed a three-dimensional contour coordinate frame into a two-dimensional contour coordinate frame represented by radial and angular axes and then adopted the computed-torque control and proportional-derivative control to design a motion controller. Wang et al. [13] proposed a contour error control design based on the generalized extended state observer and model predictive control to reduce the contour errors of networked multi-axis motion systems caused by network-induced delays. Many contour control designs that can effectively improve the contour motion accuracy of the feed drive axis have been proposed; however, most of them focus on improving contour motion properties, but the adverse influence of external disturbances on the feed drive axes has yet to be carefully discussed. Particularly, the cutting force in the thread milling operations affects the contour motion accuracy of feed drive axes while performing a helical path. Disturbance estimation and friction compensation can effectively suppress the influence of external disturbances but cannot mitigate the motion perturbation induced by model uncertainty and external disturbance variation.
For repeatedly performed positional command paths, if repeated external disturbances influence the motion of a feed drive axis, iterative learning control can effectively suppress the influence of external disturbances and provide better positional motion results to the feed drive axis. Therefore, iterative learning control was applied to design iterative learning contour control (ILCC) to improve the contour motion accuracy of feed drive axes. Li et al. [14] adopted a hidden Markov model to spatially map the command path and actual trajectory, and proposed a hybrid evaluation function to estimate the contour error of the tool position and orientation. Furthermore, they used spatial iterative learning control to compensate for the contour errors and generate numerical control commands to improve the contour accuracy of five-axis CNC machine tools. Wang et al. [15] transformed the planar motion in a Cartesian coordinate frame into a polar coordinate frame, defined the polar-radius error, and designed a real-time iterative compensation method to reduce the polar-radius error-equivalent contour error for planar motion control systems. Based on a third-order contour error prediction model, Wang et al. [16] proposed an online iterative trajectory pre-compensation control design to improve the contour motion performance of multi-axis motion systems. Considering that motion control systems are usually influenced by factors such as model parameter variations, external disturbance, and non-linear friction, Yuan and Zhao [17] adopted the asymmetric Gaussian-like function to design an adaptive learning law and then suggested contour error control using the iterative learning method. Xu et al. [18] considered the properties of motion trajectory and time-varying random disturbance to design a time-varying weighting matrix and then proposed a norm-optimal cross-coupling iterative learning control using a time-varying control method. Li et al. [19] transformed the dynamic model of the motion control system from time domain to space domain and then proposed the spatial iterative learning control. Xu et al. [20] proposed a double-iterative learning and cross-coupling control design to improve the individual axis and contour tracking performance of biaxial motion systems. Hendrawan et al. [21] proposed an ILCC design that modified numerical control programs by an iterative approach to improve the contour motion performance of the feed drive axes of a CNC machine tool. Hendrawan et al. [22] referred to contour errors and proposed an embedded ILCC that can diminish the contour errors of a CNC machine tool during feed drive motions by iteratively modifying the reference trajectory of each feed drive axis. Due to the complexity of contour error models and the requirement of a large amount of collected data and processing time while using the offline iterative learning control, Dao and Chen [23] developed an online iterative learning control that referred to equivalent contour errors as the control target. Ling et al. [24] adopted the concept of position domain control and then designed the cross-coupled iterative learning control to improve the synchronization and performance of non-linear contour tracking motions. Simba et al. [25] integrated the proportional-derivative control and disturbance observer to develop an ILCC that can simultaneously improve contour motion accuracy and degrade the influence of external disturbance on the feed drive axes of a CNC machine tool. Based on the Newton extremum-seeking algorithm, Wang et al. [26] proposed an ILCC design to project contour errors to each motion axis and then adjust the axial position reference commands using an iterative learning approach. Although the ILCC can effectively reduce the contour errors of feed drive axes while performing repeated positional command paths, the motion control design of each feed drive axis must mitigate the adverse influence of model uncertainty and external disturbance variation. In addition, feed drive axes must have similar motion dynamics and external disturbance properties to match multi-axis motions and increase contour accuracy. Moreover, an ILCC that combines contour error estimation and iterative learning control must have an easy-to-implement control structure and must not excessively increase the computational burden of motion control systems to facilitate the realization of ILCC in industrial applications.
During thread milling operations, the three feed drive axes of a CNC machine tool perform synchronous motions along a helical path so that the thread milling tool can move along the helical path for thread cutting. Therefore, the cutting force in milling operations and the motion dynamics and friction disturbance of each feed drive axis affect the contour motion accuracy of feed drive axes while performing helical paths. Existing ILCC methods can improve the contour motion accuracy of feed drive axes; however, motion control systems were required to record and process a large amount of collected data, requiring longer data processing and computing time that usually increases the computational burden. In addition, before the convergence of the applied iterative learning processes, the existing ILCC could produce a few defective products due to the low contour control accuracy, which waste materials and increase manufacturing costs. A helical path consists of multiple pitch cycles, and, thus, feed drive axes perform repeated synchronous motions in pitch cycles. Moreover, during the motions of pitch cycles in a thread milling operation, feed drive axes possess the same motion dynamics and influence of the cutting force and friction disturbance. Therefore, this study developed a pitch cycle-based ILCC (PCB-ILCC) to address the control system and thread manufacturing problems caused by applying the ILCC to thread milling operations. The PCB-ILCC applied the ILCC to pitch cycles so that the contour error produced by feed drive axes while performing synchronous motions in the previous pitch cycle can be used to compute and modify the positional commands of synchronous motions in the current pitch cycle. Thus, the PCB-ILCC, pitch by pitch, gradually improves the contour motion accuracy of feed drive axes while performing helical paths. To achieve the design goal, the PCB-ILCC must integrate with the contour error vector (CEV) estimation method without excessive computational burden, the ILCC design with a simple and easy-to-implement structure, and the driving control design of the feed drive axis with robust motion properties. Therefore, this study estimated the CEV by referring to the interpolated positions on the helical cycle and the response positions on the response cycle produced by the feed drive axes. The ILCC design with the previous cycle learning method was modified and applied to generate a command cycle for the feed drive axes. This study also considered the dynamic characteristics of the feed drive axis actuated by a permanent magnet synchronous motor (PMSM) and applied the robust control with friction and dynamic compensations to design a PMSM driving control so that the feed drive axes can have robust motion performance.
To evaluate and validate the feasibility and performance of the PCB-ILCC design, this study performed a PMSM motion control experiment on a PMSM testbench and a machine tool thread milling experiment on a five-axis CNC machining center. The results of this experiment showed that in comparison to the proportional–proportional–integral control extensively used in commercialized PMSM driving control, the rate of reduction achieved by the PCB-ILCC is significant. The reduction rates of the average, root-mean-square, and maximum contour errors achieved were 96.27%, 95.92%, and 94.76%, respectively. The machine tool thread milling experiment results indicated that the PCB-ILCC could provide good helical contour accuracy even though the feed drive axes suffered the adverse influence of the friction disturbance and cutting force. The reduction rates of the average, root-mean-square, and maximum contour errors achieved 90.56%, 89.68%, and 83.35%, respectively. Therefore, the experimental results validated that the PCB-ILCC can effectively mitigate the adverse influence of the friction disturbance and cutting force to significantly improve the contour accuracy of the feed drive axes while performing the helical path in thread milling operations.
This paper is organized as follows. Section 2 describes the design of the PCB-ILCC and the difference between the PCB-ILCC and conventional ILCC designs. Section 3 presents the CEV estimation method and the ILCC design with previous cycle learning. Section 4 presents the PMSM driving control design for the feed drive axes with robust motion performance. Section 5 presents the experimental results performed on a PMSM testbench and a five-axis CNC machining center to demonstrate the feasibility and performance of the PCB-ILCC developed in this study. Section 6 concludes this paper.

2. Design of PCB-ILCC

Three paths are considered in thread milling operations, including the helical, command, and response paths, as shown in Figure 1. Here, the helical path denotes the path formed by the thread contour to be machined, the command path represents the path referenced by the feed drive axes, and the response path denotes the path produced by the feed drive axes. For thread milling operations, the feed drive axes perform synchronous motions for the given command path so that a thread milling tool can move along the response path for thread cutting. Conventionally, the command path is given as the helical path, and the feed drive axes are controlled so that the response path can accurately follow the command to make the thread milling tool perform precise helical motions. However, due to the external disturbances and motion dynamics of feed drive axes, the response path is generally different from the command path; thus, the accuracy of helical motions is limited.
Applying conventional ILCC methods, as shown in Figure 2, the command path can be iteratively modified to reduce the contour error between the helical and response paths. Here,  P k  and  R k  denote the command and response paths at the  k th iteration, respectively. The helical path  H  and response path  R k  were used to calculate the contour error  ε k  at the  k th iteration. Then, the contour error  ε k  was recorded in memories for further use in the next iteration.  P k 1 , that has been recorded in memories, denotes the command path at the  k 1 th iteration. The recorded contour error  ε k 1  and command path  P k 1  were adopted by an iterative learning contour controller to modify the command path  P k . The command path  P k  was referenced by the feed drive axes to produce the response path  R k , and was recorded in memories for further use in the next iteration. Although conventional ILCC methods can improve the accuracy of helical motions, ILCC methods usually require huge memories to record a large amount of collected data and take longer data processing and computing time. Therefore, ILCC methods are usually difficult to implement on compact controllers with a limited number of memories and computation speed. Furthermore, ILCC methods could produce a few defective products before achieving convergence, increasing manufacturing costs.
Since a helical path is composed of multiple pitch cycles, and feed drive axes generally possess similar external disturbances and motion dynamics during pitch cycle motions, the PCB-ILCC developed in this study applies the ILCC to the pitch cycles so that the contour error of the previous pitch cycle  ε j 1  can be considered with the previous command cycle  P j 1  to modify the current command cycle  P j  through the iterative learning contour controller, as shown in Figure 3. Compared with the conventional ILCC structure shown in Figure 2, the helical path  H  must be segmented to obtain the helical cycle  H j  and the response path is naturally amalgamated by the response cycle  R j  through the motions of the feed drive axes. In addition, the operation of the pitch increment is required in the segmentation process and iterative learning contour controller to make the thread milling operation move from the current pitch cycle to the next pitch cycle. The contour error of the  j th pitch cycle  ε j  was calculated by referring to the  j th helical cycle  H j  and the  j th response cycle  R j . Both the contour error  ε j  and command cycle  P j  were recorded in different memories for the modification of the next command cycle. Therefore, the PCB-ILCC can gradually improve the accuracy of helical motions pitch by pitch and mitigate the problems of applying conventional ILCC methods to thread milling operations. Nevertheless, as indicated by Figure 3, some issues must be carefully considered in this study when applying the PCB-ILCC in practice, including the estimation of the CEV (adopted in the contour error calculation to approximately calculate the contour error of each pitch cycle), the design of the ILCC (adopted in the iterative learning contour controller to gradually reduce the contour errors of the pitch cycles), and the design of the PMSM driving control (adopted in the feed drive axes to actuate the feed drive axes with a stable and improved motion performance).

3. CEV Estimation and ILCC Design

The CEV is the shortest distance vector between the response position (on the response path) and the helical path; thus, the precise calculation of the CEV is generally difficult and time-consuming. Although many CEV estimation methods have been proposed in recent decades, the accuracy, complexity, and calculation time of the estimated CEV must be carefully considered in the PCB-ILCC design. Therefore, this study used the interpolated positions on a helical path to estimate the CEV [27]. As shown in Figure 4 R  denotes the response position on the response path;  H i  and  H i 1  denote the interpolated positions that are the positions on the helical path and close to the response position  R H c  denotes the position on the  H i 1 H i  vector, and the  R H c  vector is perpendicular to the  H i 1 H i  vector so that the  R H c  vector can be obtained as the estimated CEV  ε  (i.e.,  ε = R H c ). Here,  ε = ε x ε y ε z T , where  ε x ε y , and  ε z  are the axial components of the estimated CEV  ε  on the X-axis, Y-axis, and Z-axis, respectively.
Let  d i 1  and  d i , respectively, denote the distance from  R  to  H i 1  and  H i , and let  D i 1  and  D i , respectively, denote the distance from  H c  to  H i 1  and  H i . Referring to the geometric relationship shown in Figure 4, the distance relations can be obtained as Equations (1) and (2):
d i 1 2 D i 1 2 = d i 2 D i 2
D i 1 + D i = H i H i 1
Here,  D i 1  and  H c  can be, respectively, calculated by Equations (3) and (4):
D i 1 = d i 1 2 d i 2 + H i H i 1 2 2 H i H i 1
H c = H i 1 + H i H i 1 H i H i 1 D i 1
Therefore, the estimated CEV  ε  can be obtained as:
ε = R H c = H c R = H i 1 R + H i H i 1 H i H i 1 D i 1
Because the helical path is adopted in thread milling operations, the estimated CEV  ε  can provide a good approximation to the actual CEV when the helical path has a large helical radius, and the thread milling tool moves with a low feed rate. While the response position  R  moves on the response path  R , the contour error  ε  can be obtained by collecting the estimated CEV  ε  with the given helical path  H , as shown in Figure 2 for the conventional ILCC methods; therefore, for the simplified expression, Equation (6) is used to represent the calculation of the contour error  ε :
ε = c t r e r r   H R
where  c t r e r r        denotes the function of the contour error calculation. Specifically, as shown in Figure 3 for the PCB-ILCC method, Equation (6) is modified as Equation (7) to represent the calculation of the contour error  ε j :
ε j = c t r e r r   H j R j
This study developed the cycle learning method shown in Equation (8) to design the PCB-ILCC shown in Figure 3:
P j = P j 1 + Z L + K s   ε j 1
where  Z L  denotes the vector of the pitch increment and  K s  denotes the transfer function matrix of the iterative learning contour controllers.  Z L  is a constant vector.  K s = d i a g K m s m = x , y , z  is a diagonal matrix, and the diagonal entries  K x s K y s , and  K z s , respectively, denote the transfer functions of the iterative learning contour controllers for the feed drive axes X, Y, and Z. Because  H j = H j 1 + Z L R j = G s   P j , the contour error  ε j  shown in Equation (7) can be rewritten as Equation (9) by substituting Equation (8) to Equation (7):
ε j = c t r e r r   H j 1 + Z L G s   P j 1 G s   Z L G s K s   ε j 1
where  G s  denotes the transfer function matrix of the feed drive axes.  G s = d i a g G m s m = x , y , z  is a diagonal matrix, and the diagonal entries  G x s G y s , and  G z s  denote the transfer functions of the feed drive axes X, Y, and Z, respectively. Because  R j 1 = G s   P j 1  and  Z L = G s   Z L  are at a steady state, Equation (9) can be further rewritten as Equations (10) and (11):
ε j = c t r e r r   H j 1 R j 1 G s K s   ε j 1
ε j = ε j 1 G s K s   ε j 1
Accordingly, from Equation (11), the iterative equation of the contour error  ε j  can be obtained as Equation (12):
ε j = I G s K s   ε j 1 ,
where  I  denotes the identity matrix. Therefore, for each feed drive axis, the relationship between the contour error components in cycles  j  and  j 1  can be represented as Equation (13):
ε j , m = 1 G m s K m s   ε j 1 , m ,   m = x , y , z ,
where  ε j , x ε j , y , and  ε j , z  are the axial components of the contour error  ε j  on the X-axis, Y-axis, and Z-axis, respectively. According to Equation (13), for the cycle learning method presented in Equation (8) to achieve convergence, the inequality shown in Equation (14) must be satisfied:
1 G m s K m s < 1 ,   m = x , y , z ,
where       denotes the infinity norm. Figure 5 shows the PCB-ILCC structure in detail using the cycle learning method; here, the operation of the PCB-ILCC illustrated in Figure 5 can be described as the following steps:
  • Step 1: compute the command cycle  P j  by combining the recorded command cycle  P j 1 , the filtered contour error  K s   ε j 1 , and the pitch increment  Z L ; here, the command cycle  P j  must also be recorded to memory.
  • Step 2: feed drive axes perform motions with the command cycle  P j  so that the response cycle  R j  can be obtained from the motion results.
  • Step 3: calculate the contour error  ε j  by using the contour error calculation function  c t r e r r     , the helical cycle  H j , and the response cycle  R j ; here, the contour error  ε j  must also be recorded to memory.
  • Step 4: if the  j th pitch cycle is not the last cycle, the PCB-ILCC increases the cycle number and then moves to Step 1.
In Figure 5, the pitch cycle number  j  is increased from  j = 1  (the first cycle); moreover, for  j = 1 , the command cycle  P 1 = H 1 , and the recorded command cycle  P 0  and the recorded contour error  ε 0  are both set to be zeros. A conventional proportional–integral controller was employed in this study as the iterative learning contour controller. The proportional and integral control gains were tuned manually according to the influence of the proportional control gain on the rate of convergence and the influence of the integral control gain on the converged contour error. Moreover, the control gains must be tuned to satisfy the inequality shown in Equation (14) to achieve convergence.

4. PMSM Driving Control Design

A PMSM mainly actuates feed drive axes; thus, the PMSM driving control design can significantly affect the motion performance of each feed drive axis and further affects the operating performance of the PCB-ILCC. Considering the dynamic equation of the feed drive axis actuated by the PMSM as Equation (15):
τ m = J θ ¨ m + B θ ˙ m + τ s c + τ e ,
where  J  and  B  denote the equivalent moment of inertia and equivalent viscous friction coefficient, respectively.  θ m  denotes the rotating angular position of the PMSM,  τ m  and  τ e , respectively, denote the driving torque and external disturbance torque of the PMSM, and  τ s c  as Equation (16) denotes the compound friction torque composed of the static friction torque and Coulomb friction torque:
τ s c = τ c + τ s τ c e θ ˙ m ω s δ sgn θ ˙ m ,
where  τ s  and  τ c  denote the static friction torque and Coulomb friction torque, respectively.  ω s  denotes the Stribeck velocity, and  δ  denotes the Stribeck shape factor.
Supposing that the feed drive axis would not be influenced by the external disturbance torque (i.e.,  τ e = 0 ), this study designed the driving control as Equation (17):
τ m = J ^ a m + B ^ θ ˙ m + τ ^ s c ,
where  J ^  denotes the estimated equivalent moment of inertia,  B ^  denotes the estimated equivalent viscous friction coefficient,  τ ^ s c  denotes the estimated compound friction torque, and  a m  denotes the control term required to be further designed. Therefore, by substituting Equation (17) into Equation (15), the dynamic equation can be rewritten as Equation (18):
J θ ¨ m = J ^ a m + B ^ B θ ˙ m + τ ^ s c τ s c = J ^ a m + B ˜ θ ˙ m + τ ˜ s c ,
where  B ˜ = B ^ B  and  τ ˜ s c = τ ^ s c τ s c  denote the estimation errors of the equivalent viscous friction coefficient and compound friction torque, respectively. The estimation error of the equivalent moment of inertia can also be defined as  J ˜ = J ^ J , and Equation (18) can be rewritten as Equation (19):
θ ¨ m = J 1 J ^ a m + B ˜ θ ˙ m + τ ˜ s c = a m + J 1 J ˜ a m + B ˜ θ ˙ m + τ ˜ s c
The uncertainty  η  can be defined as Equation (20):
η = J 1 J ˜ a m + B ˜ θ ˙ m + τ ˜ s c
Then, Equation (19) can be rewritten as Equation (21):
θ ¨ m = a m + η
To overcome the influence of the uncertainty  η , this study designed the control term  a m  as shown in Equation (22):
a m = θ ¨ d K 1 θ ˙ m θ ˙ d K 0 θ m θ d + δ a m ,
where  θ ¨ d θ ˙ d , and  θ d  denote the commands for the rotating angular acceleration, angular velocity, and angular position of the PMSM, respectively.  K 1 > 0  and  K 0 > 0  denote the control gains;  δ a m  denotes the additional term and is used to attenuate the influence of the uncertainty. Based on Equations (21) and (22), this study derived the error equation as Equation (23):
θ ¨ m θ ¨ d + K 1 θ ˙ m θ ˙ d + K 0 θ m θ d = δ a m + η
Moreover, through the definition of the following error vector  e = θ m θ d θ ˙ m θ ˙ d Τ , Equation (23) can be rewritten as Equation (24):
e ˙ = A e + b δ a m + η ,
where  A  and  b  denote the error state matrix and uncertain input vector, respectively. Moreover,
A = 0 1 K 0 K 1   and   b = 0 1 .
Herein, in case the uncertainty  η  is bounded by  ρ , i.e.,
η < ρ
The additional term  δ a m  can be designed as Equation (27) so that the solution trajectory of the error equation shown in Equation (24) is bounded. The Appendix A shows the detailed derivation.
δ a m = ρ b T P e b T P e , b T P e > σ ρ σ b T P e , b T P e σ ,
where  σ  denotes a threshold constant. The matrix  A  in Equation (25) must be a stable matrix. Moreover, matrix  P = P T > 0  and matrix  Q = Q T > 0  must satisfy the Lyapunov equation as Equation (28):
A T P + P A = Q
Equations (17), (22), (27) and (28) construct the PMSM driving control designed in this study. Figure 6 shows the control structure developed in this study. The generated commands  θ d θ ˙ d , and  θ ¨ d  and the actual responses  θ m  and  θ ˙ m  are adopted to calculate the control term  a m  and the additional term  δ a m  by using Equation (22) and Equations (27) and (28), respectively. Then, Equation (17) calculates the driving torque  τ m  to actuate the feed drive axis through the PMSM. Here, the commands  θ d θ ˙ d , and  θ ¨ d  are generated from the axial component of the command cycle, and the actual response  θ m  produces the axial component of the response cycle, as shown in Figure 5.
In this study, for simplicity, the matrix  Q  is designed as  Q = q I ; where  q  is a positive and real variable. Therefore, the PMSM driving control includes five control parameters  K 1 K 0 ρ σ , and  q . The conservative design of the control parameters generally improves the stability but limits the performance of the PMSM driving control. On the other hand, to improve the performance, the design of the control parameters usually deteriorates the stability of the PMSM driving control. Therefore, due to the control parameters’ design complexity, the control parameters were tuned manually to achieve a balanced performance and stability considering their control characteristics. Referring to Equations (24) and (25), the control parameters  K 1  and  K 0  determine the dynamic response of the solution trajectory. Thus, the control parameters  K 1  and  K 0  were tuned to rapidly attenuate errors without significant oscillations. Referring to Equation (27), because the threshold constant  σ  affects the chattering phenomenon of the solution trajectory, the threshold constant  σ  was tuned to suppress chattering. Furthermore, the larger the bounded value  ρ , the more stable control is achieved; however, a large bounded value usually increases the response time and steady-state errors of the PMSM driving control. Referring to the detailed derivation in the Appendix A, the Euclidean norm of the solution trajectory at the steady state (i.e., the magnitude of the error vector at the steady state) is proportional to the product of the bounded value  ρ  and threshold constant  σ , and is inversely proportional to the control parameter  q . Therefore, to achieve stable control and improved steady-state errors, the control parameters  ρ σ , and  q  were tuned repeatedly by reducing the product of the bounded value  ρ  and threshold constant  σ  and increasing the control parameter  q .

5. Experimental Results

The PMSM motion control experiment and machine tool thread milling experiment were used to validate the PCB-ILCC approach developed in this study. An industrial computer with a steel chassis, backplane, and CPU card manufactured by Advantech Co., Ltd. was used in this study. The CPU card, which utilized a 3.4-GHz INTEL CORE i7-2600 CPU and 16-Gigabyte memory, is the main component of the industrial computer. The computer used the Linux operating system; moreover, programming language C was used to realize the control system that included a human–machine operating interface, signal measurement, feedback control, and recording and displaying of experimental data. Furthermore, the computer utilized the EtherCAT fieldbus industrial network system for sending control commands to and receiving feedback signals from the PMSM drivers at a sampling period of 0.5 ms so that the feed drive axes perform helical motions. In the experiments, the thread milling path of  M 40 × 3.5  ISO metric external threads was adopted as the helical path. Here, Equation (29) was used to calculate the helical radius  R H :
R H = 1 2 D M + D B 3 3 8 2 Ρ H ,
where  D M P H , and  D B  denote the major diameter, helical pitch of the external thread, and the blade diameter of the thread milling tool, respectively.
The following control structures were adopted in the experiments for validating the feasibility and evaluating the motion performance of the PCB-ILCC approach developed in this study.
  • PMSM-PPIC: the proportional–proportional–integral control extensively used in commercialized PMSM drivers and CNC controllers.
  • PMSM-DC: the PMSM driving control developed in this study (Section 4).
  • PCB-ILCC: the PCB-ILCC (Section 2) with the integration of the CEV estimation and ILCC design (Section 3) and the PMSM driving control (Section 4).
The following performance indices were also adopted:
  • AVG: average value.
  • SD: standard deviation.
  • RMS: root-mean-square value.
  • MAX: maximum value.

5.1. PMSM Motion Control Experiment

This study used 1.5 kW commercialized PMSM packs for the motion control experiment, as shown in Figure 7. The PMSM has a rated speed and torque of 2000 rpm and 7.16 Nm, respectively. In addition, the PMSM is coupled to a 23-bit absolute optical encoder that can measure the rotating angular position of the PMSM. This study identified the moment of inertia and the friction parameters of the PMSM through a constant speed experiment and estimation experiment [28]. Table 1 shows the identification results.
Figure 8 shows the experimental results of 11 cycles of the helical path. The vertical axis is the RMS of the contour error, and the horizontal axis is the cycle number of the helical path. Since the PMSM does not significantly suffer from external disturbances, for the results of the PMSM-PPIC, the AVG of the RMS variation is  0.129  rad with an SD of  6.1811 × 10 5  rad; hence, the ratio between the SD and AVG is  0.48 % . Accordingly, the RMS contour error of the PMSM-PPIC is almost unchanged in each cycle of the helical path. Similarly, the RMS contour error is almost unchanged because the PMSM-DC does not perform cycle learning. The AVG of the RMS variation is  0.0012  rad with an SD of  8.4230 × 10 6  rad so that the ratio of the SD and AVG is  0.69 % . Nevertheless, compared with the PMSM-PPIC, the contour errors of the PMSM-DC were significantly reduced, and the AVG of the RMS variation achieves a rate of reduction of  90.48 % . The PCB-ILCC can further reduce the RMS contour errors through cycle learning. Figure 8 shows that, after using the PCB-ILCC, the RMS contour error starts to converge at the fourth cycle; the RMS contour error is reduced from  0.0012  rad (the first cycle) to  5.2681 × 10 4  rad (the eleventh cycle), achieving a rate of reduction of  57.61 % .
Figure 9 shows the magnitude of the CEVs in the last cycle (the eleventh cycle) of Figure 8. Figure 10 compares the different control structures in the AVG, RMS, and MAX errors. The experimental results show that, compared with the PMSM-PPIC, the PMSM-DC significantly reduces contour errors, and the PCB-ILCC can further reduce contour errors. The AVG, RMS, and MAX errors for the PMSM-PPIC are  0.0117 0.0129 , and  0.0314  rad, respectively. As the PMSM-PPIC is susceptible to the influence of PMSM motion dynamics and friction disturbances, the contour errors using the PMSM-PPIC vary considerably. However, because the PMSM-DC and PCB-ILCC use the PMSM driving control developed in this study, which can not only compensate for the PMSM motion dynamics and friction disturbances but can also provide robust performance to PMSM motions, the PMSM-DC and PCB-ILCC have robust control results and small contour errors compared with the PMSM-PPIC. For the PMSM-DC, the AVG error is reduced to  0.0011  rad at a rate of reduction of  90.74 % , the RMS error is reduced to  0.0012  rad at a rate of reduction of  90.50 % , and the MAX error is reduced to  0.0031  rad at a rate of reduction of  90.18 % . Although the PMSM-DC can mitigate the adverse influence of external disturbances, few contour errors remain. Since external disturbances can repeatedly act on the PMSM while performing cycle motions along the helical path, the PCB-ILCC helps reduce contour errors. Compared with the PMSM-PPIC, the rate of reduction of the AVG, RMS, and MAX errors are  96.27 % 95.92 % , and  94.76 % , respectively. Moreover, compared with the PMSM-DC, the rate of reduction of the AVG, RMS, and MAX errors are  59.72 % 57.04 % , and  46.64 % , respectively. Therefore, the PMSM motion control experiment can validate that, for the helical path, the PCB-ILCC can effectively reduce contour errors and exhibit performance with robust control and disturbance suppression.

5.2. Machine Tool Thread Milling Experiment

This study utilized the five-axis CNC machining center, as shown in Figure 11, for the thread milling experiment to validate and evaluate the feasibility and performance of the PCB-ILCC approach developed in this study under the cutting force disturbance. The machining center mainly comprises a vertical axis (Z feed drive axis) and two planar axes (X feed drive axis and Y feed drive axis). The vertical axis and planar axes are actuated by PMSMs so that the thread milling tool can perform linear and circular motions. Each feed drive axis is equipped with a ballscrew lead of 10 mm and adopts a 1.5 kW PMSM pack. Table 2 shows the PMSM-equivalent parameters of the feed drive axes, respectively, identified by the constant speed experiment and estimation experiment [28]. The thread milling experiment used a 90 mm long thread milling tool with a tool holder diameter of 12 mm. The material of the disk-type thread milling blade was tungsten steel. The blade diameter is 20 mm with a blade angle of 60 degrees. The workpiece material is AL6061.
Figure 12 shows the experimental results while performing thread milling operations on the five-axis CNC machining center; here, because the PMSM-DC has been validated to provide superior control results with small contour errors, the PCB-ILCC was compared with the PMSM-PPIC to clarify the control performance between the proposed and conventional control methods under the adverse influence of external disturbances. The vertical axis denotes the RMS contour error of each cycle, and the horizontal axis represents the cycle number of the helical path. Figure 12 illustrates that the RMS contour error of the PMSM-PPIC is almost unchanged and remains at a large value. The AVG of the RMS variation is  10.6930  μm, the SD is  0.0286  μm, and the ratio between the SD and AVG is  0.27 % . Compared with the PMSM-PPIC, although the cutting force makes the RMS contour error fluctuate, the PCB-ILCC still significantly reduces the RMS contour error. The PCB-ILCC achieves that the RMS contour error starts to converge at the fifth cycle and is reduced from  3.0734  μm (the first cycle) to  1.1055  μm (the tenth cycle) with a rate of reduction of  64.03 % .
Figure 13 shows the magnitude of the CEVs in the last cycle (the tenth cycle) of Figure 12. Figure 14 compares the control performance of the PMSM-PPIC and PCB-ILCC. Because of the influence of friction disturbance on each feed drive axis, the contour error of the PMSM-PPIC has been significantly increased at quadrants of circular motion, with the MAX error of  21.7032  μm. In addition, because of the influence of the cutting force and motion dynamics of each feed drive axis, the suppression of the contour error is limited, and there were larger contour errors while using the PMSM-PPIC; the AVG error is  10.4015  μm and the RMS error is  10.7086  μm. Compared with the PMSM-PPIC, the PCB-ILCC can effectively suppress the adverse influence of the cutting force and friction disturbance, and improves the motion dynamics of each feed drive axis. Therefore, the PCB-ILCC can significantly reduce the contour error; the MAX error is reduced to    3.6132  μm at a rate of  83.35 % , the RMS error is reduced to  1.1055  μm with a rate of  89.68 % , and the AVG error is reduced to  0.9818  μm with a rate of  90.56 % . Therefore, the machine tool thread milling experiment indicated that, compared with the PMSM-PPIC, the PCB-ILCC could significantly improve the contour accuracy of the feed drive axes while performing the helical path in thread milling operations.

6. Conclusions

During thread milling operations, the helical contour motion accuracy of feed drive axes is critical because it can significantly affect the accuracy of manufactured threads. The contour control design can improve the contour motion accuracy of feed drive axes; however, the motion dynamics and friction disturbance of each feed drive axis and the cutting force in thread milling operations significantly affect the helical contour motion accuracy of feed drive axes. Although the ILCC can reduce contour errors while feed drive axes perform repeated motions, recording and processing large amounts of data increase the memory usage and the computational burden on the motion control system. Moreover, during the ILCC convergence, producing multiple defective products can lead to wasted materials while increasing manufacturing costs.
Because the helical path is composed of multiple pitch cycles, and feed drive axes are generally influenced similarly by motion dynamics and external disturbances such as the friction disturbance and cutting force during the motions of pitch cycles, the PCB-ILCC is developed to apply the contour error in the previous pitch cycle motion to modify the position command of the current pitch cycle motion so as to not only improve the contour motion accuracy of the feed drive axes while performing a helical path but also deal with the control system problem and thread manufacturing problem induced by applying the ILCC to thread milling operations. To realize the PCB-ILCC, the CEV estimation method must not excessively increase the computational burden, and the ILCC must be easy to implement; in addition, the PMSM driving control must mitigate the adverse influence of motion dynamics and external disturbances acting on the feed drive axes. Therefore, the PCB-ILCC estimated the CEV by referring to the interpolated positions on the helical path, designed the ILCC using a cycle learning method, and designed the PMSM driving control using the robust control method that integrates friction compensation and dynamics compensation.
The motion control experiments performed on a testbench validated that the robust PMSM driving control method developed in this study can significantly reduce contour errors compared with the extensively used proportional–proportional–integral control. As proportional–proportional–integral control is susceptible to the influence of PMSM motion dynamics and external disturbances, the contour errors vary considerably. Compared with proportional–proportional–integral control, the PCB-ILCC integrated with the PMSM driving control has robust control results and small contour errors owing to the compensation of PMSM motion dynamics and friction disturbances and the robust performance of PMSM motions. Compared with the extensively used proportional–proportional–integral control, the PCB-ILCC achieves a rate of reduction in the contour errors of over 90%. Moreover, compared with the robust PMSM driving control, the PCB-ILCC can further reduce the contour errors with a rate of reduction of over 40%. Thread milling experiments performed on a five-axis CNC machining center showed that even though the feed drive axes were adversely influenced by the motion dynamics and external disturbances, including the friction disturbance and cutting force, the PCB-ILCC achieved a significant contour error reduction because of the cycle learning method developed in this study. The reduction rate of contour errors when applying the PCB-ILCC reached 80% in thread milling operations compared to proportional–proportional–integral control. Thus, the experimental results validated that the PCB-ILCC can be applied to thread milling operations, effectively improving the helical contour motion accuracy of the feed drive axes.
The PCB-ILCC presented in this paper exhibited improved contour motion accuracy for helical contour motions. Therefore, in the future, the developed PCB-ILCC will be applied to the motion control system of a robotic manufacturing system for manufacturing mechanical parts with a non-circular helix using a robot manipulator. In future studies, some developed approaches in the presented PCB-ILCC must be carefully redesigned. For instance, for the design of the contour error vector estimation, we must consider the curvature variation in the helical path with a non-circular contour. Furthermore, for the design of the iterative learning contour control, we must consider the dynamics of the robot manipulator, and the design of the PMSM driving control must be extended to the design of robotic motion control.

Author Contributions

Conceptualization, S.-S.Y. and W.-J.J.; data curation, S.-S.Y. and W.-J.J.; formal analysis, S.-S.Y. and W.-J.J.; funding acquisition, S.-S.Y.; investigation, S.-S.Y. and W.-J.J.; methodology, S.-S.Y. and W.-J.J.; project administration, S.-S.Y.; resources, S.-S.Y.; software, W.-J.J.; supervision, S.-S.Y.; validation, S.-S.Y. and W.-J.J.; visualization, S.-S.Y. and W.-J.J.; writing—original draft, S.-S.Y. and W.-J.J.; writing—review and editing, S.-S.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the Ministry of Science and Technology, Taiwan (Grant number MOST111-2221-E-027-141).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Chen-Chou Hsieh (Shihlin Electric & Engineering Corp., Taiwan) for his valuable comments on the design of PMSM servo drives.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

A error state matrix
a m control term
B equivalent viscous friction coefficient
B ^ estimated equivalent viscous friction coefficient
B ˜ estimation error of the equivalent viscous friction coefficient
b uncertain input vector
c t r e r r       contour error calculation function
D B blade diameter
D i distance from  H c  to  H i  
D i 1 distance from  H c  to  H i 1  
D M major diameter
d i distance from  R  to  H i  
d i 1 distance from  R  to  H i 1  
δ Stribeck shape factor
δ a m additional term
e following error vector
ε contour error
ε j contour error of the  j th pitch cycle
ε j , x axial component of the contour error  ε j  on X-axis
ε j , y axial component of the contour error  ε j  on Y-axis
ε j , z axial component of the contour error  ε j  on Z-axis
ε k contour error at the  k th iteration
ε estimated CEV
ε x axial component of the estimated CEV  ε  on X-axis
ε y axial component of the estimated CEV  ε  on Y-axis
ε z axial component of the estimated CEV  ε  on Z-axis
η uncertainty
G s transfer function matrix of feed drive axes
G x s transfer function of the feed drive X-axis
G y s transfer function of the feed drive Y-axis
G z s transfer function of the feed drive Z-axis
H helical path
H j helical cycle of the  j th pitch cycle
H c position on  H i 1 H i  vector
H i i th interpolated position on a helical path
I identity matrix
J equivalent moment of inertia
J ^ estimated equivalent moment of inertia
J ˜ estimation error of the equivalent moment of inertia
j pitch cycle number
K 0 control gain with  K 0 > 0  
K 1 control gain with  K 1 > 0  
K s transfer function matrix of the iterative learning contour controllers
K x s transfer function of the iterative learning contour controller for the feed drive X-axis
K y s transfer function of the iterative learning contour controller for the feed drive Y-axis
K z s transfer function of the iterative learning contour controller for the feed drive Z-axis
k iteration number
λ min Q minimum eigenvalue of the matrix  Q  
ω s Stribeck velocity
P j command cycle of the  j th pitch cycle
P k command path at the  k th iteration
P matrix with  P = P T > 0  
P H helical pitch
Q matrix with  Q = Q T > 0  
q positive and real variable
R response path
R j response cycle of the  j th pitch cycle
R k response path at the  k th iteration
R response position on response path
R H helical radius
ρ bounded value of uncertainty
S A set in following error space
S B set in following error space
σ threshold constant
τ m driving torque
τ c Coulomb friction torque
τ e external disturbance torque
τ s static friction torque
τ s c compound friction torque
τ ^ s c estimated compound friction torque
τ ˜ s c estimation error of the compound friction torque
θ d rotating angular position command
θ m rotating angular position
V Lyapunov function
Z L vector of pitch increment

Appendix A

Referring to the error equation as Equation (24), this study defines the Lyapunov function as Equation (A1):
V = e Τ P e ,
where the matrix  P = P T > 0 . By differentiating Equation (A1), the derivative of the Lyapunov function can be obtained as Equation (A2):
V ˙ = e ˙ Τ P e + e Τ P e ˙
By substituting Equation (24), Equation (A2) can be further rearranged as Equation (A3):
V ˙ = e Τ A Τ P e + e Τ P A e + 2 e Τ P b δ a m + η
Then, referring to Equation (28), Equation (A3) can be rewritten as Equation (A4):
V ˙ = e Τ Q e + 2 b Τ P e Τ δ a m + η
Because the uncertainty  η  is bounded by  ρ  as Equation (26), Equation (A5) can be deduced from Equation (A4):
V ˙ < e Τ Q e + 2 b Τ P e Τ δ a m + ρ b Τ P e b Τ P e
Case (A): in the case that the additional term  δ a m  is designed as Equation (27) with  b T P e > σ , then  δ a m = ρ b T P e b T P e .
Since  2 b T P e T δ a m + ρ b T P e b T P e = 0  and  Q = Q T > 0 V ˙ < 0 .
Case (B): If the additional term  δ a m  is designed as Equation (27) with  b T P e σ , then  δ a m = ρ σ b T P e , and thus  2 b T P e T δ a m + ρ b T P e b T P e = 2 ρ σ b T P e 2 + 2 ρ b T P e . Therefore, Equation (A5) can be rewritten as Equation (A6):
V ˙ < e Τ Q e + 2 ρ σ b T P e 2 + 2 ρ b T P e
Since the minimum eigenvalue of the matrix  Q λ min Q , satisfies the inequality  λ min Q e 2 e T Q e , to achieve  V ˙ < 0 , Equation (A7) must be satisfied:
2 ρ σ b T P e 2 2 ρ b T P e + λ min Q e 2 > 0
Since Equation (A7) is a quadratic inequality with a minimal value of  ρ σ 2 + λ min Q e 2  when  b T P e = σ 2 , there exists a matrix  Q  such that Equation (A8) is satisfied:
ρ σ 2 + λ min Q e 2 > 0
Then,  V ˙ < 0  and  e > ρ σ 2 λ min Q .
The description of Case (A) denotes that there exists a set  S A  so that the solution trajectory of Equation (24) will reach and enter the set  S A  since  V ˙ < 0  is outside  S A . The description of Case (B) denotes that there exists a set  S B  with  e ρ σ 2 λ min Q  and  S B S A . Therefore, the descriptions of Case (A) and Case (B) denote that the solution trajectory of Equation (24) will be bounded within the region between set  S A  and set  S B .

References

  1. Liu, Y.; Sun, W.; Buccella, C.; Cecati, C. Robust control of dual-linear-motor-driven gantry stage for coordinated contouring tasks based on feed velocity. IEEE Trans. Ind. Electron. 2023, 70, 6229–6238. [Google Scholar] [CrossRef]
  2. Chen, S.L.; Khong, M.H.; Hsieh, S.M. Contouring control of a five-axis machine tool with equivalent errors. Electronics 2022, 11, 2521. [Google Scholar] [CrossRef]
  3. Zhao, H.; Li, X.; Ge, K.; Ding, H. A contour error definition, estimation approach and control structure for six-dimensional robotic machining tasks. Robot. Comput.-Integr. Manuf. 2022, 73, 102235. [Google Scholar] [CrossRef]
  4. Liu, Y.; Wan, M.; Xiao, Q.; Qin, X. Combined predictive and feedback contour error control with dynamic contour error estimation for industrial five-axis machine tools. IEEE Trans. Ind. Electron. 2021, 69, 6668–6677. [Google Scholar] [CrossRef]
  5. Bo, L.; Taiyong, W.; Peng, W. Cross-coupled control based on real-time double circle contour error estimation for biaxial motion system. Meas. Control 2021, 54, 324–335. [Google Scholar] [CrossRef]
  6. Li, J.; Qi, C.; Li, Y.; Wu, Z. Prediction and compensation of contour error of CNC systems based on LSTM neural-network. IEEE/ASME Trans. Mechatron. 2021, 27, 572–581. [Google Scholar] [CrossRef]
  7. Song, D.N.; Ma, J.W.; Zhong, Y.G.; Yao, J.J. Definition and estimation of joint-space contour error based on generalized curve for five-axis contour following control. Precis. Eng. 2020, 65, 32–43. [Google Scholar] [CrossRef]
  8. Wang, Z.; Hu, C.; Zhu, Y.; Zhang, M.; Zhang, C. The modular design of trajectory compensation based on ATCF for precision motion control. Mech. Syst. Signal Process. 2020, 135, 106393. [Google Scholar] [CrossRef]
  9. Zhang, B.; Cao, R.; Hou, Z. The model-free adaptive cross-coupled control for two-dimensional linear motor. Trans. Inst. Meas. Control 2020, 42, 1059–1069. [Google Scholar] [CrossRef]
  10. Kuang, Z.; Sun, G.; Gao, H. Simplified Newton-based CEE and discrete-time fractional-order sliding-mode CEC. IEEE/ASME Trans. Mechatron. 2019, 24, 175–185. [Google Scholar] [CrossRef]
  11. Ma, J.W.; Song, D.N.; Jia, Z.Y.; Jiang, W.W.; Wang, F.J.; Liu, W. An equivalent-plane cross-coupling position control for contour-accuracy improvement in three-axis free-form contour following tasks. J. Dyn. Syst. Meas. Control Trans. ASME 2019, 141, 041007. [Google Scholar] [CrossRef]
  12. Shi, R.; Lou, Y. Three-dimensional contouring control: A task polar coordinate frame approach. IEEE Access 2019, 7, 63626–63637. [Google Scholar] [CrossRef]
  13. Wang, Y.W.; Zhang, W.A.; Dong, H.; Yu, L. A GESO based MPC approach to contour error control of networked motion control system. Int. J. Syst. Sci. 2019, 50, 2216–2225. [Google Scholar] [CrossRef]
  14. Li, J.; Yue, R.; Fei, Y. Five-axis contour error control based on numerical control data. Machines 2023, 11, 85. [Google Scholar] [CrossRef]
  15. Wang, Z.; Zhou, R.; Hu, C.; Zhu, Y. Real-time iterative compensation based contouring control method for polar coordinate motion systems. IEEE/ASME Trans. Mechatron. 2022, 27, 3517–3526. [Google Scholar] [CrossRef]
  16. Wang, Z.; Hu, C.; Zhu, Y.; Zhu, L. Prediction-model-based contouring error iterative precompensation scheme for precision multiaxis motion systems. IEEE/ASME Trans. Mechatron. 2021, 26, 2274–2284. [Google Scholar] [CrossRef]
  17. Yuan, H.; Zhao, X. Advanced contouring compensation approach via Newton-ILC and adaptive jerk control for biaxial motion system. IEEE Trans. Ind. Electron. 2021, 69, 5081–5090. [Google Scholar] [CrossRef]
  18. Xu, W.; Hou, J.; Li, J.; Yuan, C.; Simeone, A.; Xu, W.; Li, J. Multi-axis motion control based on time-varying norm optimal cross-coupled iterative learning. IEEE Access 2020, 8, 124802–124811. [Google Scholar] [CrossRef]
  19. Li, J.; Wang, Y.; Li, Y.; Luo, W. Reference trajectory modification based on spatial iterative learning for contour control of two-axis NC systems. IEEE/ASME Trans. Mechatron. 2020, 25, 1266–1275. [Google Scholar] [CrossRef]
  20. Xu, W.; Hou, J.; Yang, W.; Wang, C. A double-iterative learning and cross-coupling control design for high-precision motion control. Arch. Electr. Eng. 2019, 68, 427–442. [Google Scholar] [CrossRef]
  21. Hendrawan, Y.M.; Farrage, A.; Uchiyama, N. Iterative NC program modification and energy saving for a CNC machine tool feed drive system with linear motors. Int. J. Adv. Manuf. Technol. 2019, 102, 3543–3562. [Google Scholar] [CrossRef]
  22. Hendrawan, Y.M.; Simba, K.R.; Uchiyama, N. Iterative learning based trajectory generation for machine tool feed drive systems. Robot. Comput.-Integr. Manuf. 2018, 51, 230–237. [Google Scholar] [CrossRef]
  23. Dao, V.T.; Chen, S.L. Helical contouring control with online iterative learning control. In Proceedings of the 2018 IEEE/ASME International Conference on Advanced Intelligent Mechatronics, Auckland, New Zealand, 9–12 July 2018; pp. 599–603. [Google Scholar] [CrossRef]
  24. Ling, J.; Feng, Z.; Yao, D.; Xiao, X. Non-linear contour tracking using feedback PID and feedforward position domain cross-coupled iterative learning control. Trans. Inst. Meas. Control 2018, 40, 1970–1982. [Google Scholar] [CrossRef]
  25. Simba, K.R.; Bui, B.D.; Msukwa, M.R.; Uchiyama, N. Robust iterative learning contouring controller with disturbance observer for machine tool feed drives. ISA Trans. 2018, 75, 207–215. [Google Scholar] [CrossRef]
  26. Wang, Z.; Hu, C.; Zhu, Y.; He, S.; Zhang, M.; Mu, H. Newton-ILC contouring error estimation and coordinated motion control for precision multiaxis systems with comparative experiments. IEEE Trans. Ind. Electron. 2018, 65, 1470–1480. [Google Scholar] [CrossRef]
  27. Hendrawan, Y.M.; Simba, K.R.; Uchiyama, N. Iterative learning contouring controller based on trajectory generation with linearly interpolated contour error estimation and Bézier reposition trajectory for computerized numerical control machine tool feed drive systems. Adv. Mech. Eng. 2019, 11, 1687814019868108. [Google Scholar] [CrossRef]
  28. Yeh, S.S.; Su, H.C. Development of friction identification methods for feed drives of CNC machine tools. Int. J. Adv. Manuf. Technol. 2011, 52, 263–278. [Google Scholar] [CrossRef]
Figure 1. Relationship between helical path, command path, and response path.
Figure 1. Relationship between helical path, command path, and response path.
Applsci 13 06447 g001
Figure 2. Control structure of ILCC methods applying to thread milling operations.
Figure 2. Control structure of ILCC methods applying to thread milling operations.
Applsci 13 06447 g002
Figure 3. Control structure of the PCB-ILCC method developed in this study.
Figure 3. Control structure of the PCB-ILCC method developed in this study.
Applsci 13 06447 g003
Figure 4. CEV estimation based on the interpolated positions on the helical path.
Figure 4. CEV estimation based on the interpolated positions on the helical path.
Applsci 13 06447 g004
Figure 5. Detailed PCB-ILCC structure developed in this study.
Figure 5. Detailed PCB-ILCC structure developed in this study.
Applsci 13 06447 g005
Figure 6. PMSM driving control structure.
Figure 6. PMSM driving control structure.
Applsci 13 06447 g006
Figure 7. Testbench with PMSM packs used in the motion control experiment.
Figure 7. Testbench with PMSM packs used in the motion control experiment.
Applsci 13 06447 g007
Figure 8. Variations in RMS of contour errors (PMSM motion control experiment).
Figure 8. Variations in RMS of contour errors (PMSM motion control experiment).
Applsci 13 06447 g008
Figure 9. Contour error variations in PMSM motion control experiment.
Figure 9. Contour error variations in PMSM motion control experiment.
Applsci 13 06447 g009
Figure 10. Control performance comparison in PMSM motion control experiment.
Figure 10. Control performance comparison in PMSM motion control experiment.
Applsci 13 06447 g010
Figure 11. Five-axis CNC machining center used in thread milling experiment.
Figure 11. Five-axis CNC machining center used in thread milling experiment.
Applsci 13 06447 g011
Figure 12. Variations in RMS of contour errors (machine tool thread milling experiment).
Figure 12. Variations in RMS of contour errors (machine tool thread milling experiment).
Applsci 13 06447 g012
Figure 13. Contour error variations in machine tool thread milling experiment.
Figure 13. Contour error variations in machine tool thread milling experiment.
Applsci 13 06447 g013
Figure 14. Control performance comparison in machine tool thread milling experiment.
Figure 14. Control performance comparison in machine tool thread milling experiment.
Applsci 13 06447 g014
Table 1. PMSM parameters of the testbench.
Table 1. PMSM parameters of the testbench.
ParametersX-AxisY-AxisZ-Axis
Moment of inertia (kg∙m2)   0.0016   0.0015   0.0015
Viscous friction coefficient (Nm/rad/s)   0.0074   0.0071   0.0064
Static friction torque (Nm)   0.1037   0.1043   0.1079
Coulomb friction torque (Nm)   0.0890   0.0916   0.0929
Stribeck velocity (rad/s)   0.2588   0.3146   0.3406
Stribeck shape factor   1.5302   1.7817   1.5339
Table 2. PMSM-equivalent parameters of feed drive axes.
Table 2. PMSM-equivalent parameters of feed drive axes.
ParametersX-AxisY-AxisZ-Axis
Moment of inertia (kg∙m2)   0.0030   0.0026   0.0024
Viscous friction coefficient (Nm/rad/s)   0.0177   0.0514   0.0232
Static friction torque (Nm)   0.6977   0.6228   1.9821
Coulomb friction torque (Nm)   0.6061   0.5680   1.1565
Stribeck velocity (rad/s)   0.5926   2.0243   1.2274
Stribeck shape factor   1.7386   2.1397   1.0293
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yeh, S.-S.; Jiang, W.-J. Development of Pitch Cycle-Based Iterative Learning Contour Control for Thread Milling Operations in CNC Machine Tools. Appl. Sci. 2023, 13, 6447. https://doi.org/10.3390/app13116447

AMA Style

Yeh S-S, Jiang W-J. Development of Pitch Cycle-Based Iterative Learning Contour Control for Thread Milling Operations in CNC Machine Tools. Applied Sciences. 2023; 13(11):6447. https://doi.org/10.3390/app13116447

Chicago/Turabian Style

Yeh, Syh-Shiuh, and Wei-Jia Jiang. 2023. "Development of Pitch Cycle-Based Iterative Learning Contour Control for Thread Milling Operations in CNC Machine Tools" Applied Sciences 13, no. 11: 6447. https://doi.org/10.3390/app13116447

APA Style

Yeh, S. -S., & Jiang, W. -J. (2023). Development of Pitch Cycle-Based Iterative Learning Contour Control for Thread Milling Operations in CNC Machine Tools. Applied Sciences, 13(11), 6447. https://doi.org/10.3390/app13116447

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop