# Online Prediction of Milling Inner Hole Roundness Error Based on Accurate SSEM Value Extraction

Authors: Zhenbang Hu; Gedong Jiang; Xuesong Mei; Xialun Yun; Yun Zhang
Journal: Shock and Vibration
 (2019)

---

## Abstract

To improve the machining accuracy and production efficiency of precision components with deep hole structures, an online prediction method of the inner hole roundness error, which cannot be directly measured in real time during the machining process, is proposed in this paper. For online prediction of the workpiece roundness error (WRE) during machining, a predictive model based on correlation analysis and a proportional method is proposed according to the spindle synchronous error motion (SSEM) by three-probe method testing. To improve the prediction accuracy of the WRE, a particle swarm optimization (PSO) algorithm is introduced for optimizing a probe mounting angle of a three-probe method, and a harmonic wavelet method for SSEM feature extraction is proposed. Using the PSO algorithm, the optimal probe mounting angle of the three-probe method is obtained, the influence of spindle surface roundness on SSEM is eliminated, and the higher-order harmonic suppression of the three-probe method is avoided effectively. By the harmonic wavelet method, the accurate SSEM extraction is enhanced and the WRE prediction accuracy is promoted. The experiments show that the inner hole roundness error online prediction method proposed in this paper has high prediction accuracy.


---

## Body

## 1. Introduction

Hydraulic cylinders are the key actuators in hydraulic systems and are widely used in the aerospace, automobile, and medical device industries. With further industrial developments, there is a tremendous `demand for hydraulic cylinders in the manufacturing industry. This brings with it the challenge of high-accuracy machining and high-efficiency manufacturing of hydraulic cylinder components. Typically, an inner reciprocating cylinder moves along an outer cylinder, which is essentially a deep inner hole structure; therefore, high-precision machining is required to manufacture the outer cylinder. In particular, the surface roundness error of the inner hole is important due to its direct impact on the working performance and fault rates of hydraulic cylinders.Traditional technology used for manufacturing inner holes is limited to drilling, boring, reaming, or broaching processes, which can be affected by the environment; moreover, precision is difficult to guarantee. Some advantages can be achieved by combining different technologies; however, there are still shortcomings, such as complicated processing methods and low efficiency. For high-precision machining of inner holes, various approaches have been presented and focus on improving the processing technology. Special tools have been integrated into various processing technologies designed as possible solutions. Zhang et al. studied on a customized gun drill regrinding system and the effect of apex offset variation on the hole straightness deviation though gun drilling experiments, and the research results show that the hole straightness deviation is the critical performance indicator in deep hole gun drilling process, and it is very difficult to achieve precise control for the gun drill tool geometry, so the customized drilling tool and the gun drill measurement system are hard to be widely popularized and applied [1]. Chin et al. investigated the effect of internally boring and trepanning association (BTA) drilling round profile, and they used the ratio of frequencies for modeling of deep hole geometry and the harmonic model for predicting the hole profile, but the method needs to use a special tool for BTA drilling, and the prediction model is unquantifiable processing [2]. Zhang et al. proposed a boring and trepanning association deep hole drilling mechanism, which uses a specialized staggered-edge drilling head with an asymmetric structure for machining holes with a large length-to-diameter ratio, and the influence of different factors on the functional behavior of deep hole machining and machining quality was investigated [3]. Aziz et al. developed a microcompound tool consisting of a drilling part and microdiamond-electroplated-grinding part to minimize burr formation and improve the surface roughness of the drilled holes [4]. However, developing special tools for different hole sizes results in excessive costs, and the methods often cannot be applied to other applications. High-speed milling (HSM) is an advanced processing technology that offers excellent efficiency and good quality and has gradually replaced surface grinding. HSM is a promising method as a half-finishing and finishing process for inner holes and can be used to machine inner holes of different sizes under different processing environments, with the potential to greatly improve hole machining efficiency. Nonetheless, surface roundness error of the inner hole remains a bottleneck to improving the precision of inner holes. Measuring and monitoring the surface roundness error during the milling of inner holes is therefore significant to improving the precision of HSM inner holes.At present, surface roundness error measurements of inner holes are preformed offline, and milling efficiency is very low. Thus, studies of online methods to predict the inner hole roundness error are required. In milling, the spindle rotates relative to the average axis, and the radial error motion is defined as the relative displacement of the spindle from the average axis in the perpendicular direction on a certain axial plan. In general, the radial error motion is measured in a certain axial plane and includes both synchronous and asynchronous error motions. Synchronous error motion is an integer multiple of the rotational frequency, reflects the mean contour of the total error motion, and has good repeatability under certain conditions. During milling, the cutter rotates synchronously and repetitively with the spindle; therefore, the surface roundness error of the inner hole can be mapped by the spindle synchronous error motion. Ryua et al. used the character to flat end milling for generating the plane surface, and the surface texture characteristics are evaluated and predicted through a series of vector transformations considering tool run out and tool setting error [5]. Marsh et al. found that the workpiece diameter was affected by frequency content and phase of spindle, and they used simple grinding force models to predict best roundness of workpiece; the experimental results show the rate of roundness improvement matches theoretical predictions [6]. So the spindle error motion has a direct impact on WRE. With the spindle error motion measurement accuracy improvement, the predictive capability of workpiece roundness error will be enhanced.The three-probe method is commonly used for spindle radial error motion measurements and is more appropriate for in situ measurements than conventional methods such as the Donaldson reversal method [7] and multiorientation method [8]. To accurately measure the error motion of the spindle using the three-probe method, the angles between the three measurement probes must be optimized and play a decisive role. The mounting angle of the probes determines whether the method is suitable for the required scenario, and a suitable angle combination could prevent the occurrence of higher-order harmonic suppression. Therefore, the influence of the probe mounting angles has widely been investigated. Obi et al. systematically demonstrated a sequential three-probe method and investigated the relevance of the measurement angle and coefficient of roundness, in addition to analyzing the characteristics and measurement errors of various angle combinations [9]. Similarly, Okuyama et al. demonstrated high-accuracy radial motion measurements using the three-probe method, along with a weighting equation and inverse filtering [10, 11]. Shi et al. proposed a hybrid method to correctly estimate roundness error and considered the effects of stochastic disturbances using different measurement angles (0°, 42.8°, and 126.4°) [12]. Moreover, Grejda implemented a unique adaptation of the Donaldson reversal technique to advance the state of spindle metrology using different angle combinations (0°, 122.5°, and 202.5°) [13]. The methods are only applicable to the experimental environment, and the steps in these methods are cumbersome; in addition, the measurement accuracy needs to be further improved. In this paper, an optimization is performed to acquire the mounting angle of the probes, as well as the separation accuracy, and the performance of the novel angle combinations are compared to existing angle combinations.The extraction of synchronous features of the measurement signal can have an important influence on predicting the inner hole roundness error. Synchronous feature extraction is centered on the characteristics associated with the integer multiples of the rotational frequency. Moreover, the synchronous feature extraction process can eliminate interference caused by noise in the signal. Many studies have focused on rotating machinery fault diagnosis using single signal analysis based on synchronous feature extraction. Siegel et al. used a tachometer-less synchronously averaged envelope method that uses signal processing and feature extraction techniques to diagnose the health and degradation of a rolling element bearings based on the synchronous characteristic signal [14]. Wang et al. performed signal feature extraction of rolling bearings with early weak faults: first, using the EEMD method to decompose the original signal into several IMFs; then, selecting the IMF with biggest kurtosis index value; and finally, using the tunable Q-factor wavelet transform for the fault signal feature extraction [15]. Synchronous feature extraction can also be used for signal purification. In this way, An et al. employed a mathematical morphology filter and combined the morphology filter with a simple algorithm, facile calculation, and effective filter to eliminate noise interference and purify the signal; nonetheless, the method was still affected by the characteristic of structural elements [16]. Zhang et al. combined synchronous feature selection and parameter optimization using an ant colony algorithm and support vector machine, to differentiate between different fault conditions, and demonstrated good final diagnostic accuracy for the intelligent fault diagnosis of rotating machinery [17]. Numerous methods have attempted feature extraction from characteristic signals. Chen et al. used the daubechies wavelet transform [18], and Lei et al. used the EEMD method [19]. However, most of these methods require complex algorithms; moreover, the harsh working conditions of many processing environments make it difficult to adopt these methods on-site for measuring rotating machinery error motion. Therefore, the harmonic wavelet method has been introduced to extract the synchronous features of measurement signals since it is suitable for in situ use and capable of accurate signal feature extraction, and importantly, can retain the basic characteristics of the measurement signal.The online prediction method of the inner hole roundness error is needed for manufacturing inner holes as hydraulic cylinders. The method proposed in this paper has advantages of cost-effectiveness, high adaptability, and high accuracy. The work to improve prediction property mainly lies in the high-accuracy SSEM measuring and extracting method. The accurate predictive model is also essential. The rest of the paper is organized as follows: in Section2, a novel method for the online prediction of the inner hole roundness error is introduced. In Section 3, the SSEM measurement and feature extraction method are presented, with a focus on improving measurement accuracy. Exactness of the SSEM values were calculated. In Section 4, a predictive model based on a correlation analysis and experimental data is presented. In Section 5, results of experimental validation are presented, and in Section 6, discussions on the various angle combinations are provided. Finally, the conclusions are summarized in Section 7.


## 2. Method for WRE Online Prediction

The high-speed spindle and cutting tool rotate together during the inner hole milling process; therefore, if the spindle and the cutting tool are ideally connected, their motions are synchronized and consistent. The internal surface of the workpiece is affected by cutting tool error motion; in particular, the workpiece roundness profile is highly influenced by the cutting tool error motion during hole milling. It follows that the error motion radius will lead to a corresponding roundness profile radius.In the general case, the roundness profile of the machined internal surface is definite, owing to the repeatability of spindle radial error motion. The SSEM is extracted from the spindle radial error motion and plays a crucial role in this repeatability. This paper proposes a novel method for WRE online prediction using the accurate SSEM value along with the mathematical model. The procedure can be outlined as follows:(1)
Construct the experimental system for accurate spindle radial error motion measurement, which must satisfy the requirements for online measurements under certain machining conditions. Select a measurement method that can eliminate the interference factor and effectively improve measurement accuracy.(2)
Extract the spindle radial error motion signals to obtain the SSEM. Synchronous components of the signal have a much larger effect on WRE, and extracting the SSEM can effectively improve the signal-to-noise ratio (SNR) and prediction accuracy.(3)
Establish the predictive model using the extracted SSEM and workpiece roundness profile and determine whether there is a correlation between them.(4)
Predict the WRE according to the predictive models and experimentally validate the SSEM under the same machining conditions.


## 3. Accurate Measurement of SSEM

The proposed method aims to solve critical challenges that have an important influence on the final predicted result. One particularly important factor is the SSEM value. The first step in determining the SSEM is to acquire the spindle radial error motion. The key elements of this step are selecting a suitable measurement angle combination and extracting an accurate SSEM value.


### 3.1. Three-Probe Method for Measuring Spindle Radial Error Motion

The three-probe method is often used to calibrate the rotational accuracy of precision spindles and to acquire the surface roundness profile of an object; however, the method should not be limited to these applications. In this study, the three-probe method is used to measure spindle radial error motion and eliminate surface roundness profile. Consequently, the spindle radial error motion can be obtained with relatively high accuracy.To perform this method, three probes are arranged at certain angles (θ1,θ2,θ3) within the same cross section of a particular axial location. The intersection of the different directions is identified as the artifact origin O; moreover, rL is the radius of the least square circle for measuring rotor surface.The displacement probes measure the relative displacement in the sensitive direction between the probe and surface of the specified axial location of the rotor. The surface roundness profile can be expressed asrθθ∈0,2π, and the displacement from the probe to the artifact origin can be defined as Ri (i = 1, 2, 3). Supposing the radial error motion of the spindle is comprised of the axis components exθ and eyθ and considering the effects of the radial error motion in the sensitive direction, the signal measured by the displacement probe can be described by the following relational expression:(1)siθ=Ri−rL−rθ+θi−exθcosθi−eyθsinθi,i=1,2,3;θ∈0,2π.In general, one of the probe directions will overlap with thex-axis when θ1=0°. The relationship between the radial error motion, surface roundness profile, and measurement signal can be expressed as(2)exθ=s1θ−rθ,(3)eyθ=s2θ−rθ+θ2−s1θcosθ2−rθcosθ2sinθ2.From equations (2) and (3), the critical calculation of radial error motion can be obtained from the expression of rθ, and in addition, the weighted combination process adopts the weighted coefficient ai (i = 1, 2, 3), which transforms the measurement signal into a weighted combination function using a1=1, a2=sinθ3/sinθ2−θ3, and a3=sinθ2/sinθ3−θ2.At this point, the weighted combination functionsθ only contains information about the surface roundness profile characteristics and demonstrates that sθ has periodicity; therefore, Sλ=∫−∞∞sθe−jλθdθ is the frequency-domain expression of the weighted combination function of the measurement signal, λ is the harmonic number, and the weighted function Fλ can be given as(4)Fλ=a1ejλ0+a2ejλθ2+a3ejλθ3.According to equation (4), rθ can be derived using transposition and the inverse Fourier transform as(5)rθ=12π∫−∞∞SλFλejλθdλ.Thus, once the surface roundness profile of a specified axial location has been acquired from equation (5), by substituting equation (5) into equations (2) and (3), the radial error motion components exθ and eyθ can be determined.In addition, the three-probe method can be used to separate surface roundness profile and radial error motion from the measurement signals since it is also an error separation method; however, the measurement angle is the critical component in obtaining accurate separation results using the above analysis. Improper probe mounting angles will cause higher-order harmonic suppression to occur, and it will not be possible to obtain predetermined separation results from the measurements; thus, it is vital to select an appropriate set of measurement angles.In the separation process, the weighted functionFλ and weighted combination function Sλ are crucial to obtaining the surface roundness profile rθ, and the expression of the critical Fλ will directly influence future results.


### 3.2. Enhancing the Measurement Precision

#### 3.2.1. Optimization of Defined Target

In the weighted functionFλ, the measurement angle is the denominator, as shown in equation (5), such that a zero value of Fλ causes the method to lose significance; therefore, each value of the weighted function must be kept far from zero. First, it is necessary to ensure that each individual weighted value is kept far away from zero. If the minimum value of Fλ is far from zero, the other values must also be far from zero. Thus, if an arbitrary angle combination is given, as λ changes, there must be a corresponding sequence of values, Fλ. Find the minimum value, expressed as Fλ−∞. The notation ⋅−∞ represents the minimum value of the numeric sequence. Since the angle combinations change, there will be many sequences Fλ and many minimum values from which the farthest from zero should be chosen. In other words, finding the maximum value can be expressed as Fλ−∞∞, where the notation ⋅∞ represents the maximum value of the numeric sequence. Here, Aλ=Fλ−∞∞ is used to express the maximum value. Furthermore, the entire weighted function, i.e., the summation of the individual weighted values, should also be kept far away from zero and can be expressed as Bλ=Fλ1∞, where the notation ⋅1 represents the summation of all values.Hence, the weighted function optimization is a multiobjective optimization, which incorporates the locality and global optimization goals, which are not independent and have mutual effects. To solve the multiobjective optimization problem, both the locality and global optimization goals must be considered. Then, the linear weighted combination method can be used to construct the combined optimization goals given as(6)Cλ=αAλ+βBλ=αFλ−∞∞+βFλ1∞,where α and β are the weighting factors and must satisfy α+β=1. The values of α and β are determined based on the importance of the optimization goals. Generally, the locality optimization goals have a much greater influence on the separation accuracy and results for larger values of α.


#### 3.2.2. Optimized Angle Based on PSO

The PSO method is a multivariate and global optimization method with a number of advantages, including quick optimization, high efficiency, and reasonable accuracy, and thus has significance in many engineering applications.The solution to the optimization problem corresponds to a particle in the PSO. The optimal searching process is used to first initialize a group of random particles, often referred to as the initial particle swarm, wherein each particle has a randomly generated initial value for its location and velocity. The particle swarm invariably follows the current optimal particle and iteratively updates the optimal particle to search for the optimal solution. During the iterative searching process, each particle is constantly updating its location and velocity within the search scope, and the updates during each iterative process are based on both the current optimal solution of the particle and the particle swarm. The current optimal solution of the particle is often referred to as the personal bestpb, and analogously, the current optimal solution of the particle swarm is referred to as the global best gb. The location and velocity of the particle are updated according to(7)vijk+1=σvijk+c1r1pbijk−lijk+c2r2gbijk−lijk,lijk+1=lijk+vijk+1,where i is the particle number; j is the location or velocity dimension; k is the number of iterations; σ is the inertia weight; lijk and vijk are the location and velocity of the particle, respectively; c1 and c2 are the learning factors, usually c1=c2=2; and r1 and r2 are the random parameters in the range [0, 1]. The method for obtaining the optimal solution of the angle combination based on the PSO search process is illustrated in Figure 1. The optimization procedure is as follows:(1)
Initialize the particle parameters and set the initial angle combinations and angular increments in advance(2)
Calculate the solution for the target function by using particle parameters and record them(3)
Obtain the local and global optimization targets based on the current particle parameters(4)
Update the particle parameters and identify the optimal solution(5)
Compare the current optimal value with the previous value and select the larger one(6)
Repeat steps (2) to (5) until the termination conditions are metFigure 1
Process diagram for finding the optimal solution of the angle combination based on PSO.As shown in Figure1, the measurement angles θ2 and θ3 are the search parameters, and the search range is (0°, 360°). In applying the PSO method, the angle combination of (0°, θ2, θ3) refers to the particle’s location lijk, and the number of particles are randomly generated as i such that a series of angle combinations (0°, θ2i, θ3i) can be acquired. In other words, the original particle swarm and angular increment (0°, dθ2i, dθ3i) are equivalent to the particle velocity vijk.In the optimal solution searching process, every particle locationlijk corresponds to a serialized set of weighted functions Fλ, and the corresponding value of Cλ can be obtained by continuously searching for the optimal solution, while updating the location and velocity of the particle according to equation (7). pbijk is the particle location that fits the requirements of the locality optimized solution Aλ, and gbijk is the particle location that meets the requirements of the globally optimized solution for the combination of optimization goals Cλ.During the optimization process, as the current angle combination gets closer and closer to the optimal solution, the angular increments become smaller and smaller; therefore, the stopping criterion of the PSO is when the deviation of the angular increment is smaller than a predefined value or a preset number of iterations is reached. Finally, using the PSO approach, the optimal measurement angle was acquired as (0°, 45°, 96.5°).


### 3.3. Feature Extraction for Accurately Predicting SSEM

The rotor of a high-speed spindle and cutting tool rotates together during the milling process; therefore, if error motion of the rotor and the cutting tool are ideally connected, they are synchronized and consistent. The milling workpiece internal surface is influenced by the cutting tool error motion, and in particular, the workpiece roundness profile is highly relevant to cutting tool error motion during hole milling. It follows that a larger error motion radius will lead to a larger roundness profile radius and vice versa.In the general case, the roundness profile of the machined internal surface is repeatable due to the repetitive radial error motion. Therefore, the synchronous error motion of the cutting tool plays a particularly crucial role in determining the internal surface roundness profile. Hence, the synchronous error motion of the rotor could be used to predict the machined surface roundness, based on the feature extraction method, which extracts the synchronous error motion from the total error motion.Thus, the accurate SSEM extraction problem can be transformed into a feature signal extraction problem. The harmonic wavelet was introduced to investigate the feature signal extraction of integer multiples of the rotational frequency since the characteristics of the extracted target frequency signal include strong background noise.


#### 3.3.1. Accurate SSEM Extraction

The harmonic wavelet is a complex wavelet, first proposed by Newland [20], comprising a complex transformation of the time-domain and frequency-domain functions. The harmonic wavelet is suitable for extracting the signal of interest even when the signal is weak and difficult to acquire.The harmonic wavelet transform can easily be realized in the frequency domain for practical engineering applications. First, Fourier transformation is applied to the signal. Then, it is transformed into the harmonic wavelet in the frequency domain. Finally, the inverse Fourier transform is used to reconstruct the signal in the time domain.By investigating the feature signal extraction method, the ultimate aim is to acquire repeatable synchronous error motionesrθ, which is comprised of the integer-harmonics signals, esxθ and esyθ, extracted from the total error motions exθ and eyθ. Using the horizontal direction as an illustration, exθ is transformed into the frequency domain using Fourier transformation, and Heximi,ni,λ is obtained from the harmonic wavelet, where mi,ni∈R+, and mi<ni. Hence, one of the integer-harmonic signals can be acquired by performing an inverse Fourier transformation on Heximi,ni,λ, resulting in(8)heximi,ni,θ=12π∫−∞∞Heximi,ni,λejλθdλ.According to the above approach, the sequential extraction of the integer-harmonic signals can be realized according to the number of harmonicsi changed. The horizontal direction synchronous error motion esxθ is the superposition of all the integer-harmonic signals given as(9)esxθ=∑i=1Nheximi,ni,θ.Noise is extremely common in industrial engineering, and in particular, measurement signals acquired during the milling process usually contain significant noise components. Unfortunately, it is difficult to extract feature signals for a low SNR condition using conventional methods. In such cases, the harmonic wavelet is useful for performing feature signal extraction in the presence of strong background noise; hence, extracting both the exact amplitude and phase angle of integer-harmonics signals and fraction-harmonics signals is performed.


#### 3.3.2. Calculation of SSEM Value

To verify whether the method is effective for feature extraction, spectrum analysis was performed on the measurement signal and extracted signal. The experiments were used to validate the accurate SSEM feature extraction and performed under idle conditions. The measurement signals acquired at a rotational speed of 4800 r/min were separated into the surface roundness profile and spindle radial error motionerθ using the three-probe method. Subsequently, SSEM esrθ was extracted from spindle radial error motion erθ, and a frequency spectrum analysis was carried out. The spindle radial error motion and SSEM are illustrated in Figures 2 and 3, respectively.Figure 2
Spindle radial error motion separated from the measurement signal using the three-probe method.Figure 3
SSEM extracted from the spindle radial error motion using the harmonic wavelet.The SSEM demonstrates better repeatability in the polar plot and the integer-harmonics are more distinct in the frequency domain; moreover, the extracted signals are smoother, and the influence of noise is negligible in the time domain, as well as the frequency domain.Hence, the harmonic wavelet meets the extraction requirements for spindle radial error motion, and the harmonic wavelet is suitable for feature signal extraction in the presence of strong background noise. In addition, the repeatability and periodicity of the total error motion can be determined.The SSEM value is calculated as the difference between the radii of the maximum external circlermec and minimum inner circle rmic. The spindle radial error motion value was calculated at 4.3 μm, as shown in Figure 2, and the corresponding SSEM value was calculated at 3.58 μm, as shown in Figure 3. The signals contain noise which will lead to lower accuracy since noisy signals expand the outer boundary and reduce the inner boundary, resulting in an inner minimum circle with a smaller radius and a maximum external circle with a larger radius. The difficulty degree for feature signal extraction in complex algorithms and the harsh working conditions is significantly reduced compared to previous approaches [16–19]. Thus, feature signal extraction can effectively improve the accuracy of the SSEM value, thereby improving the overall accuracy of the predictive model.


## 3.1. Three-Probe Method for Measuring Spindle Radial Error Motion

The three-probe method is often used to calibrate the rotational accuracy of precision spindles and to acquire the surface roundness profile of an object; however, the method should not be limited to these applications. In this study, the three-probe method is used to measure spindle radial error motion and eliminate surface roundness profile. Consequently, the spindle radial error motion can be obtained with relatively high accuracy.To perform this method, three probes are arranged at certain angles (θ1,θ2,θ3) within the same cross section of a particular axial location. The intersection of the different directions is identified as the artifact origin O; moreover, rL is the radius of the least square circle for measuring rotor surface.The displacement probes measure the relative displacement in the sensitive direction between the probe and surface of the specified axial location of the rotor. The surface roundness profile can be expressed asrθθ∈0,2π, and the displacement from the probe to the artifact origin can be defined as Ri (i = 1, 2, 3). Supposing the radial error motion of the spindle is comprised of the axis components exθ and eyθ and considering the effects of the radial error motion in the sensitive direction, the signal measured by the displacement probe can be described by the following relational expression:(1)siθ=Ri−rL−rθ+θi−exθcosθi−eyθsinθi,i=1,2,3;θ∈0,2π.In general, one of the probe directions will overlap with thex-axis when θ1=0°. The relationship between the radial error motion, surface roundness profile, and measurement signal can be expressed as(2)exθ=s1θ−rθ,(3)eyθ=s2θ−rθ+θ2−s1θcosθ2−rθcosθ2sinθ2.From equations (2) and (3), the critical calculation of radial error motion can be obtained from the expression of rθ, and in addition, the weighted combination process adopts the weighted coefficient ai (i = 1, 2, 3), which transforms the measurement signal into a weighted combination function using a1=1, a2=sinθ3/sinθ2−θ3, and a3=sinθ2/sinθ3−θ2.At this point, the weighted combination functionsθ only contains information about the surface roundness profile characteristics and demonstrates that sθ has periodicity; therefore, Sλ=∫−∞∞sθe−jλθdθ is the frequency-domain expression of the weighted combination function of the measurement signal, λ is the harmonic number, and the weighted function Fλ can be given as(4)Fλ=a1ejλ0+a2ejλθ2+a3ejλθ3.According to equation (4), rθ can be derived using transposition and the inverse Fourier transform as(5)rθ=12π∫−∞∞SλFλejλθdλ.Thus, once the surface roundness profile of a specified axial location has been acquired from equation (5), by substituting equation (5) into equations (2) and (3), the radial error motion components exθ and eyθ can be determined.In addition, the three-probe method can be used to separate surface roundness profile and radial error motion from the measurement signals since it is also an error separation method; however, the measurement angle is the critical component in obtaining accurate separation results using the above analysis. Improper probe mounting angles will cause higher-order harmonic suppression to occur, and it will not be possible to obtain predetermined separation results from the measurements; thus, it is vital to select an appropriate set of measurement angles.In the separation process, the weighted functionFλ and weighted combination function Sλ are crucial to obtaining the surface roundness profile rθ, and the expression of the critical Fλ will directly influence future results.


## 3.2. Enhancing the Measurement Precision

### 3.2.1. Optimization of Defined Target

In the weighted functionFλ, the measurement angle is the denominator, as shown in equation (5), such that a zero value of Fλ causes the method to lose significance; therefore, each value of the weighted function must be kept far from zero. First, it is necessary to ensure that each individual weighted value is kept far away from zero. If the minimum value of Fλ is far from zero, the other values must also be far from zero. Thus, if an arbitrary angle combination is given, as λ changes, there must be a corresponding sequence of values, Fλ. Find the minimum value, expressed as Fλ−∞. The notation ⋅−∞ represents the minimum value of the numeric sequence. Since the angle combinations change, there will be many sequences Fλ and many minimum values from which the farthest from zero should be chosen. In other words, finding the maximum value can be expressed as Fλ−∞∞, where the notation ⋅∞ represents the maximum value of the numeric sequence. Here, Aλ=Fλ−∞∞ is used to express the maximum value. Furthermore, the entire weighted function, i.e., the summation of the individual weighted values, should also be kept far away from zero and can be expressed as Bλ=Fλ1∞, where the notation ⋅1 represents the summation of all values.Hence, the weighted function optimization is a multiobjective optimization, which incorporates the locality and global optimization goals, which are not independent and have mutual effects. To solve the multiobjective optimization problem, both the locality and global optimization goals must be considered. Then, the linear weighted combination method can be used to construct the combined optimization goals given as(6)Cλ=αAλ+βBλ=αFλ−∞∞+βFλ1∞,where α and β are the weighting factors and must satisfy α+β=1. The values of α and β are determined based on the importance of the optimization goals. Generally, the locality optimization goals have a much greater influence on the separation accuracy and results for larger values of α.


### 3.2.2. Optimized Angle Based on PSO

The PSO method is a multivariate and global optimization method with a number of advantages, including quick optimization, high efficiency, and reasonable accuracy, and thus has significance in many engineering applications.The solution to the optimization problem corresponds to a particle in the PSO. The optimal searching process is used to first initialize a group of random particles, often referred to as the initial particle swarm, wherein each particle has a randomly generated initial value for its location and velocity. The particle swarm invariably follows the current optimal particle and iteratively updates the optimal particle to search for the optimal solution. During the iterative searching process, each particle is constantly updating its location and velocity within the search scope, and the updates during each iterative process are based on both the current optimal solution of the particle and the particle swarm. The current optimal solution of the particle is often referred to as the personal bestpb, and analogously, the current optimal solution of the particle swarm is referred to as the global best gb. The location and velocity of the particle are updated according to(7)vijk+1=σvijk+c1r1pbijk−lijk+c2r2gbijk−lijk,lijk+1=lijk+vijk+1,where i is the particle number; j is the location or velocity dimension; k is the number of iterations; σ is the inertia weight; lijk and vijk are the location and velocity of the particle, respectively; c1 and c2 are the learning factors, usually c1=c2=2; and r1 and r2 are the random parameters in the range [0, 1]. The method for obtaining the optimal solution of the angle combination based on the PSO search process is illustrated in Figure 1. The optimization procedure is as follows:(1)
Initialize the particle parameters and set the initial angle combinations and angular increments in advance(2)
Calculate the solution for the target function by using particle parameters and record them(3)
Obtain the local and global optimization targets based on the current particle parameters(4)
Update the particle parameters and identify the optimal solution(5)
Compare the current optimal value with the previous value and select the larger one(6)
Repeat steps (2) to (5) until the termination conditions are metFigure 1
Process diagram for finding the optimal solution of the angle combination based on PSO.As shown in Figure1, the measurement angles θ2 and θ3 are the search parameters, and the search range is (0°, 360°). In applying the PSO method, the angle combination of (0°, θ2, θ3) refers to the particle’s location lijk, and the number of particles are randomly generated as i such that a series of angle combinations (0°, θ2i, θ3i) can be acquired. In other words, the original particle swarm and angular increment (0°, dθ2i, dθ3i) are equivalent to the particle velocity vijk.In the optimal solution searching process, every particle locationlijk corresponds to a serialized set of weighted functions Fλ, and the corresponding value of Cλ can be obtained by continuously searching for the optimal solution, while updating the location and velocity of the particle according to equation (7). pbijk is the particle location that fits the requirements of the locality optimized solution Aλ, and gbijk is the particle location that meets the requirements of the globally optimized solution for the combination of optimization goals Cλ.During the optimization process, as the current angle combination gets closer and closer to the optimal solution, the angular increments become smaller and smaller; therefore, the stopping criterion of the PSO is when the deviation of the angular increment is smaller than a predefined value or a preset number of iterations is reached. Finally, using the PSO approach, the optimal measurement angle was acquired as (0°, 45°, 96.5°).


## 3.2.1. Optimization of Defined Target

In the weighted functionFλ, the measurement angle is the denominator, as shown in equation (5), such that a zero value of Fλ causes the method to lose significance; therefore, each value of the weighted function must be kept far from zero. First, it is necessary to ensure that each individual weighted value is kept far away from zero. If the minimum value of Fλ is far from zero, the other values must also be far from zero. Thus, if an arbitrary angle combination is given, as λ changes, there must be a corresponding sequence of values, Fλ. Find the minimum value, expressed as Fλ−∞. The notation ⋅−∞ represents the minimum value of the numeric sequence. Since the angle combinations change, there will be many sequences Fλ and many minimum values from which the farthest from zero should be chosen. In other words, finding the maximum value can be expressed as Fλ−∞∞, where the notation ⋅∞ represents the maximum value of the numeric sequence. Here, Aλ=Fλ−∞∞ is used to express the maximum value. Furthermore, the entire weighted function, i.e., the summation of the individual weighted values, should also be kept far away from zero and can be expressed as Bλ=Fλ1∞, where the notation ⋅1 represents the summation of all values.Hence, the weighted function optimization is a multiobjective optimization, which incorporates the locality and global optimization goals, which are not independent and have mutual effects. To solve the multiobjective optimization problem, both the locality and global optimization goals must be considered. Then, the linear weighted combination method can be used to construct the combined optimization goals given as(6)Cλ=αAλ+βBλ=αFλ−∞∞+βFλ1∞,where α and β are the weighting factors and must satisfy α+β=1. The values of α and β are determined based on the importance of the optimization goals. Generally, the locality optimization goals have a much greater influence on the separation accuracy and results for larger values of α.


## 3.2.2. Optimized Angle Based on PSO

The PSO method is a multivariate and global optimization method with a number of advantages, including quick optimization, high efficiency, and reasonable accuracy, and thus has significance in many engineering applications.The solution to the optimization problem corresponds to a particle in the PSO. The optimal searching process is used to first initialize a group of random particles, often referred to as the initial particle swarm, wherein each particle has a randomly generated initial value for its location and velocity. The particle swarm invariably follows the current optimal particle and iteratively updates the optimal particle to search for the optimal solution. During the iterative searching process, each particle is constantly updating its location and velocity within the search scope, and the updates during each iterative process are based on both the current optimal solution of the particle and the particle swarm. The current optimal solution of the particle is often referred to as the personal bestpb, and analogously, the current optimal solution of the particle swarm is referred to as the global best gb. The location and velocity of the particle are updated according to(7)vijk+1=σvijk+c1r1pbijk−lijk+c2r2gbijk−lijk,lijk+1=lijk+vijk+1,where i is the particle number; j is the location or velocity dimension; k is the number of iterations; σ is the inertia weight; lijk and vijk are the location and velocity of the particle, respectively; c1 and c2 are the learning factors, usually c1=c2=2; and r1 and r2 are the random parameters in the range [0, 1]. The method for obtaining the optimal solution of the angle combination based on the PSO search process is illustrated in Figure 1. The optimization procedure is as follows:(1)
Initialize the particle parameters and set the initial angle combinations and angular increments in advance(2)
Calculate the solution for the target function by using particle parameters and record them(3)
Obtain the local and global optimization targets based on the current particle parameters(4)
Update the particle parameters and identify the optimal solution(5)
Compare the current optimal value with the previous value and select the larger one(6)
Repeat steps (2) to (5) until the termination conditions are metFigure 1
Process diagram for finding the optimal solution of the angle combination based on PSO.As shown in Figure1, the measurement angles θ2 and θ3 are the search parameters, and the search range is (0°, 360°). In applying the PSO method, the angle combination of (0°, θ2, θ3) refers to the particle’s location lijk, and the number of particles are randomly generated as i such that a series of angle combinations (0°, θ2i, θ3i) can be acquired. In other words, the original particle swarm and angular increment (0°, dθ2i, dθ3i) are equivalent to the particle velocity vijk.In the optimal solution searching process, every particle locationlijk corresponds to a serialized set of weighted functions Fλ, and the corresponding value of Cλ can be obtained by continuously searching for the optimal solution, while updating the location and velocity of the particle according to equation (7). pbijk is the particle location that fits the requirements of the locality optimized solution Aλ, and gbijk is the particle location that meets the requirements of the globally optimized solution for the combination of optimization goals Cλ.During the optimization process, as the current angle combination gets closer and closer to the optimal solution, the angular increments become smaller and smaller; therefore, the stopping criterion of the PSO is when the deviation of the angular increment is smaller than a predefined value or a preset number of iterations is reached. Finally, using the PSO approach, the optimal measurement angle was acquired as (0°, 45°, 96.5°).


## 3.3. Feature Extraction for Accurately Predicting SSEM

The rotor of a high-speed spindle and cutting tool rotates together during the milling process; therefore, if error motion of the rotor and the cutting tool are ideally connected, they are synchronized and consistent. The milling workpiece internal surface is influenced by the cutting tool error motion, and in particular, the workpiece roundness profile is highly relevant to cutting tool error motion during hole milling. It follows that a larger error motion radius will lead to a larger roundness profile radius and vice versa.In the general case, the roundness profile of the machined internal surface is repeatable due to the repetitive radial error motion. Therefore, the synchronous error motion of the cutting tool plays a particularly crucial role in determining the internal surface roundness profile. Hence, the synchronous error motion of the rotor could be used to predict the machined surface roundness, based on the feature extraction method, which extracts the synchronous error motion from the total error motion.Thus, the accurate SSEM extraction problem can be transformed into a feature signal extraction problem. The harmonic wavelet was introduced to investigate the feature signal extraction of integer multiples of the rotational frequency since the characteristics of the extracted target frequency signal include strong background noise.


### 3.3.1. Accurate SSEM Extraction

The harmonic wavelet is a complex wavelet, first proposed by Newland [20], comprising a complex transformation of the time-domain and frequency-domain functions. The harmonic wavelet is suitable for extracting the signal of interest even when the signal is weak and difficult to acquire.The harmonic wavelet transform can easily be realized in the frequency domain for practical engineering applications. First, Fourier transformation is applied to the signal. Then, it is transformed into the harmonic wavelet in the frequency domain. Finally, the inverse Fourier transform is used to reconstruct the signal in the time domain.By investigating the feature signal extraction method, the ultimate aim is to acquire repeatable synchronous error motionesrθ, which is comprised of the integer-harmonics signals, esxθ and esyθ, extracted from the total error motions exθ and eyθ. Using the horizontal direction as an illustration, exθ is transformed into the frequency domain using Fourier transformation, and Heximi,ni,λ is obtained from the harmonic wavelet, where mi,ni∈R+, and mi<ni. Hence, one of the integer-harmonic signals can be acquired by performing an inverse Fourier transformation on Heximi,ni,λ, resulting in(8)heximi,ni,θ=12π∫−∞∞Heximi,ni,λejλθdλ.According to the above approach, the sequential extraction of the integer-harmonic signals can be realized according to the number of harmonicsi changed. The horizontal direction synchronous error motion esxθ is the superposition of all the integer-harmonic signals given as(9)esxθ=∑i=1Nheximi,ni,θ.Noise is extremely common in industrial engineering, and in particular, measurement signals acquired during the milling process usually contain significant noise components. Unfortunately, it is difficult to extract feature signals for a low SNR condition using conventional methods. In such cases, the harmonic wavelet is useful for performing feature signal extraction in the presence of strong background noise; hence, extracting both the exact amplitude and phase angle of integer-harmonics signals and fraction-harmonics signals is performed.


### 3.3.2. Calculation of SSEM Value

To verify whether the method is effective for feature extraction, spectrum analysis was performed on the measurement signal and extracted signal. The experiments were used to validate the accurate SSEM feature extraction and performed under idle conditions. The measurement signals acquired at a rotational speed of 4800 r/min were separated into the surface roundness profile and spindle radial error motionerθ using the three-probe method. Subsequently, SSEM esrθ was extracted from spindle radial error motion erθ, and a frequency spectrum analysis was carried out. The spindle radial error motion and SSEM are illustrated in Figures 2 and 3, respectively.Figure 2
Spindle radial error motion separated from the measurement signal using the three-probe method.Figure 3
SSEM extracted from the spindle radial error motion using the harmonic wavelet.The SSEM demonstrates better repeatability in the polar plot and the integer-harmonics are more distinct in the frequency domain; moreover, the extracted signals are smoother, and the influence of noise is negligible in the time domain, as well as the frequency domain.Hence, the harmonic wavelet meets the extraction requirements for spindle radial error motion, and the harmonic wavelet is suitable for feature signal extraction in the presence of strong background noise. In addition, the repeatability and periodicity of the total error motion can be determined.The SSEM value is calculated as the difference between the radii of the maximum external circlermec and minimum inner circle rmic. The spindle radial error motion value was calculated at 4.3 μm, as shown in Figure 2, and the corresponding SSEM value was calculated at 3.58 μm, as shown in Figure 3. The signals contain noise which will lead to lower accuracy since noisy signals expand the outer boundary and reduce the inner boundary, resulting in an inner minimum circle with a smaller radius and a maximum external circle with a larger radius. The difficulty degree for feature signal extraction in complex algorithms and the harsh working conditions is significantly reduced compared to previous approaches [16–19]. Thus, feature signal extraction can effectively improve the accuracy of the SSEM value, thereby improving the overall accuracy of the predictive model.


## 3.3.1. Accurate SSEM Extraction

The harmonic wavelet is a complex wavelet, first proposed by Newland [20], comprising a complex transformation of the time-domain and frequency-domain functions. The harmonic wavelet is suitable for extracting the signal of interest even when the signal is weak and difficult to acquire.The harmonic wavelet transform can easily be realized in the frequency domain for practical engineering applications. First, Fourier transformation is applied to the signal. Then, it is transformed into the harmonic wavelet in the frequency domain. Finally, the inverse Fourier transform is used to reconstruct the signal in the time domain.By investigating the feature signal extraction method, the ultimate aim is to acquire repeatable synchronous error motionesrθ, which is comprised of the integer-harmonics signals, esxθ and esyθ, extracted from the total error motions exθ and eyθ. Using the horizontal direction as an illustration, exθ is transformed into the frequency domain using Fourier transformation, and Heximi,ni,λ is obtained from the harmonic wavelet, where mi,ni∈R+, and mi<ni. Hence, one of the integer-harmonic signals can be acquired by performing an inverse Fourier transformation on Heximi,ni,λ, resulting in(8)heximi,ni,θ=12π∫−∞∞Heximi,ni,λejλθdλ.According to the above approach, the sequential extraction of the integer-harmonic signals can be realized according to the number of harmonicsi changed. The horizontal direction synchronous error motion esxθ is the superposition of all the integer-harmonic signals given as(9)esxθ=∑i=1Nheximi,ni,θ.Noise is extremely common in industrial engineering, and in particular, measurement signals acquired during the milling process usually contain significant noise components. Unfortunately, it is difficult to extract feature signals for a low SNR condition using conventional methods. In such cases, the harmonic wavelet is useful for performing feature signal extraction in the presence of strong background noise; hence, extracting both the exact amplitude and phase angle of integer-harmonics signals and fraction-harmonics signals is performed.


## 3.3.2. Calculation of SSEM Value

To verify whether the method is effective for feature extraction, spectrum analysis was performed on the measurement signal and extracted signal. The experiments were used to validate the accurate SSEM feature extraction and performed under idle conditions. The measurement signals acquired at a rotational speed of 4800 r/min were separated into the surface roundness profile and spindle radial error motionerθ using the three-probe method. Subsequently, SSEM esrθ was extracted from spindle radial error motion erθ, and a frequency spectrum analysis was carried out. The spindle radial error motion and SSEM are illustrated in Figures 2 and 3, respectively.Figure 2
Spindle radial error motion separated from the measurement signal using the three-probe method.Figure 3
SSEM extracted from the spindle radial error motion using the harmonic wavelet.The SSEM demonstrates better repeatability in the polar plot and the integer-harmonics are more distinct in the frequency domain; moreover, the extracted signals are smoother, and the influence of noise is negligible in the time domain, as well as the frequency domain.Hence, the harmonic wavelet meets the extraction requirements for spindle radial error motion, and the harmonic wavelet is suitable for feature signal extraction in the presence of strong background noise. In addition, the repeatability and periodicity of the total error motion can be determined.The SSEM value is calculated as the difference between the radii of the maximum external circlermec and minimum inner circle rmic. The spindle radial error motion value was calculated at 4.3 μm, as shown in Figure 2, and the corresponding SSEM value was calculated at 3.58 μm, as shown in Figure 3. The signals contain noise which will lead to lower accuracy since noisy signals expand the outer boundary and reduce the inner boundary, resulting in an inner minimum circle with a smaller radius and a maximum external circle with a larger radius. The difficulty degree for feature signal extraction in complex algorithms and the harsh working conditions is significantly reduced compared to previous approaches [16–19]. Thus, feature signal extraction can effectively improve the accuracy of the SSEM value, thereby improving the overall accuracy of the predictive model.


## 4. Modeling with Correlation Analysis of SSEM and WRE

### 4.1. Experimental Setup

The proposed method was validated experimentally and the relationship between WRE and SSEM values during the milling process were examined. The experiment was used to characterize the direct influence of the error motion and surface roundness error and to derive a mathematical relationship between them. Inner shallow holes were machined at different rotational speeds, and the spindle radial error motion under the milling process was acquired.


#### 4.1.1. SSEM Measurement

Experiments were conducted using a vertical-type machine tool (BCH850, SFMTP). The experimental setup is shown in Figure4. A multichannel synchronization signal acquisition system, developed by our group, was also used. The core component (USB-4432, NI corporation) has a 24 bit resolution and was used to acquire analog signals and convert them into digital signals. The system used three displacement sensing systems (NCDT3010-M-S1) and S1 eddy current probes (Micro-Epsilon) with a range of 1 mm, and a resolution of 0.05 μm was used for obtaining measurements. A signal conditioner was used for signal conditioning and signal amplification, and a stand was used to hold the probes at the optimal angles (0°, 45°, 96.5°). The program for controlling the system was developed in LabVIEW, using the signal acquire, process, plot, and store functions.Figure 4
Experimental setup.
(a)
(b)The sampling frequency of the recorded displacement signals was 10 kHz, and the signal acquisition card ensured the signals were based on the different channel measurements from each probe, with synchronous sampling at evenly spaced angular increments. In other words, the number of sampled points per revolution remained consistent even as the rotational speed changed.In the angular domain, every sampling point corresponded to an angle of uniform distribution; therefore, the spindle radial error motion, or more specifically, the accurate SSEM, was displayed as the relative position in the angular domain. Thus, the position of the accurate SSEM corresponds to the position of the cutting tool under the milling process in the angular domain, which means the accurate SSEM corresponds to the workpiece surface roundness profile in the angular domain.


#### 4.1.2. Roundness Measurement to Determine WRE

The measurement signals of the probes distributed by different angles were acquired by the system at a stable rotational spindle speed during the milling process and separated into the spindle radial error motion using the three-probe method presented in Section3.1. Then, feature extraction was performed to obtain the spindle radial error motion and SSEM.The SSEM was obtained, and the corresponding workpiece surface roundness profile was measured using a roundness measuring instrument (YD200, LHCT company). In the measuring process, the machining workpiece first needed to be placed in the appropriate position. The inner hole roundness of machining workpiece measurement is shown in Figure5.Figure 5
Roundness measurements using YD200.


### 4.2. Establishment of Predictive Model

The predictive model is based on actual measurements. Signals were acquired from a spindle working to cut the workpiece at a stable rotational speed of 4200 r/min, and accurate SSEM values is shown in Figure6(a). Measurements of the inner hole roundness of the machined workpiece is presented in Figure 6(b).Figure 6
Measurements obtained during the machining process. (a) SSEM. (b) Workpiece roundness profile.
(a)
(b)It can be seen from Figure6 that the accurate SSEM has a similar shape to the workpiece roundness profile. To verify this, a similarity analysis of the image histograms was performed and to determine the degree of similarity. The histogram, usually referred to as a fingerprint image, depicts the intrinsic properties of the image and is not affected by size, position, or rotational angle.


#### 4.2.1. Correlation of Measurement Results

To investigate morphological similarities between accurate SSEM and workpiece roundness profile, the plots were extracted, converted to grayscale images (190 × 200 pixels), and histograms of the grayscale images were generated. As an example, histograms for a spindle rotational speed of 4200 r/min are shown in Figure7.Figure 7
Initial images and histograms. Grayscale image of (a) accurate SSEM and (b) workpiece roundness profile. Histogram of (c) accurate SSEM and (d) workpiece roundness profile.
(a)
(b)
(c)
(d)From Figure7, it can be seen that the distributions of pixels amongst the grayscale values are close to each other in both histograms. The degree of similarity between histograms can be judged in many ways, including distance judgment methods commonly selected for comparing image morphologies. The histogram intersection distance and the distance proposed by Bhattacharyya et al. can indicate similarity between histograms and only differ in their expression.The histogram intersection distance reflects the number of repeated pixels in the histogram, and the result can be normalized to be within the range [0, 1] by dividing the histogram intersection distance by the number of pixels in one of the histograms. Similarity between histograms is higher as the normalized result gets closer to 1. The normalized histogram intersection distanceDnhi can be calculated as(10)Dnhi=∑i=1NpminHGSi,HGRi∑i=1NpHGSi,where HGSi is the histogram of the accurate SSEM, HGRi is the histogram of workpiece roundness profile, and Np is the total number of pixels. The normalized histogram intersection distance Dnhi is 0.9989.Analogously, the Bhattacharyya distance reflects the probability distribution of the pixels of different histograms, and its normalized functionDnB can be given as(11)DnB=1−∑i=1NpHGSi⋅HGRi∑i=1NpHGSi⋅∑i=1NpHGRi.In this case, the normalized histogram intersection distanceDnB is 0.0047.According to the distance judgment methods, the similarity between the accurate SSEM and workpiece roundness profile shows that the normalized histogram intersection distance is very close to 1 and the normalized Bhattacharyya distance is close to 0. It follows that, under a steady rotational speed, both similarity judgment methods indicate that the histograms have a high degree of similarity in relation to image morphology.


#### 4.2.2. Predictive Model

To establish the predictive model via numerical analysis, the relationship between the SSEM value and the WRE must be determined, as well as the relationship between data sequences of SSEM and workpiece roundness profile. The SSEM value and the WRE can be calculated from the SSEM and workpiece roundness profile, respectively. The SSEM value, the WRE, and the ratio between them are presented in Table1. A correlation analysis can be introduced to obtain the relationship between data sequences of the accurate SSEM and workpiece roundness profile. The correlation coefficient for the two correlated sequences can be defined as(12)ρesrθ,rtwθ=∑θθendesrθ−esrθ¯⋅rtwθ−rtwθ¯∑θθendesrθ−esrθ¯2⋅∑θθendrtwθ−rtwθ¯2.Table 1
Results of proportional analysis and correlation analysis.
Rotational speed (r/min)
SSEM valueResr (μm)
WRERtw (μm)
RatioKRa=Resr/Rtw
Correlation coefficientρesr,rtw
4200
5.04
5.75
0.88
0.91The correlation coefficient for the SSEM value and WRE was acquired and is presented in Table1.As seen in Table1, the SSEM value and WRE are directly proportional and the proportionality coefficient KRa is 0.88. Assuming the SSEM value is known, the proportionality coefficient Resr can be used to establish the predictive model, and WRE can be expressed as Rtw=cpResr=1/KRa⋅Resr=1.14Resr. The proportional relationship between the SSEM value and WRE is a comprehensive technical indicator but only reflects the numerical relationship in a single dimension. This oversimplified model for online WRE prediction is not accurate enough.Furthermore, the correlation coefficientρesrθ,rtwθ is 0.91, which is very large and suggests the SSEM and workpiece roundness profile are highly correlated, according to the numerical analysis. The correlation coefficient consistently reflects trends in the signal sequences of the SSEM and the workpiece roundness profile more detailed, equivalent to extending the dimensions of the analysis, and can therefore be used as a correction factor in the predictive model along with proportionality coefficient for WRE, which comprehensively reflects the performance of the tool. Based on the proportionality coefficient and correction factor, the predictive model can be established as(13)Rtw=cpcResr=1KRa⋅ρesrθ,rtwθ⋅Resr=1.25Resr.The necessary correlation between the SSEM value and WRE was found to exist, and the relationship can be defined using a precise mathematical formula, which can then be incorporated into the predictive model. At the same time, based on proportionality and the correlation analyses, the predictive model can be obtained. Both the proportionality coefficient and the correlation coefficient play an important role in the predictive model.According to the description above, when the inner surface quality of the machined workpiece is difficult to determine by direct measurement during the machining process, the WRE is the most concerned parameter in deep hole machining and can be indirectly predicted using the online prediction method. For SSEM and workpiece roundness profile measurements exhibiting high morphological similarity, a correlation analysis was performed. Based on the results, the predictive model was established and used for online prediction of the inner hole roundness error. Thus, the method can improve upon existing online monitoring techniques and enhance machining accuracy.


## 4.1. Experimental Setup

The proposed method was validated experimentally and the relationship between WRE and SSEM values during the milling process were examined. The experiment was used to characterize the direct influence of the error motion and surface roundness error and to derive a mathematical relationship between them. Inner shallow holes were machined at different rotational speeds, and the spindle radial error motion under the milling process was acquired.


### 4.1.1. SSEM Measurement

Experiments were conducted using a vertical-type machine tool (BCH850, SFMTP). The experimental setup is shown in Figure4. A multichannel synchronization signal acquisition system, developed by our group, was also used. The core component (USB-4432, NI corporation) has a 24 bit resolution and was used to acquire analog signals and convert them into digital signals. The system used three displacement sensing systems (NCDT3010-M-S1) and S1 eddy current probes (Micro-Epsilon) with a range of 1 mm, and a resolution of 0.05 μm was used for obtaining measurements. A signal conditioner was used for signal conditioning and signal amplification, and a stand was used to hold the probes at the optimal angles (0°, 45°, 96.5°). The program for controlling the system was developed in LabVIEW, using the signal acquire, process, plot, and store functions.Figure 4
Experimental setup.
(a)
(b)The sampling frequency of the recorded displacement signals was 10 kHz, and the signal acquisition card ensured the signals were based on the different channel measurements from each probe, with synchronous sampling at evenly spaced angular increments. In other words, the number of sampled points per revolution remained consistent even as the rotational speed changed.In the angular domain, every sampling point corresponded to an angle of uniform distribution; therefore, the spindle radial error motion, or more specifically, the accurate SSEM, was displayed as the relative position in the angular domain. Thus, the position of the accurate SSEM corresponds to the position of the cutting tool under the milling process in the angular domain, which means the accurate SSEM corresponds to the workpiece surface roundness profile in the angular domain.


### 4.1.2. Roundness Measurement to Determine WRE

The measurement signals of the probes distributed by different angles were acquired by the system at a stable rotational spindle speed during the milling process and separated into the spindle radial error motion using the three-probe method presented in Section3.1. Then, feature extraction was performed to obtain the spindle radial error motion and SSEM.The SSEM was obtained, and the corresponding workpiece surface roundness profile was measured using a roundness measuring instrument (YD200, LHCT company). In the measuring process, the machining workpiece first needed to be placed in the appropriate position. The inner hole roundness of machining workpiece measurement is shown in Figure5.Figure 5
Roundness measurements using YD200.


## 4.1.1. SSEM Measurement

Experiments were conducted using a vertical-type machine tool (BCH850, SFMTP). The experimental setup is shown in Figure4. A multichannel synchronization signal acquisition system, developed by our group, was also used. The core component (USB-4432, NI corporation) has a 24 bit resolution and was used to acquire analog signals and convert them into digital signals. The system used three displacement sensing systems (NCDT3010-M-S1) and S1 eddy current probes (Micro-Epsilon) with a range of 1 mm, and a resolution of 0.05 μm was used for obtaining measurements. A signal conditioner was used for signal conditioning and signal amplification, and a stand was used to hold the probes at the optimal angles (0°, 45°, 96.5°). The program for controlling the system was developed in LabVIEW, using the signal acquire, process, plot, and store functions.Figure 4
Experimental setup.
(a)
(b)The sampling frequency of the recorded displacement signals was 10 kHz, and the signal acquisition card ensured the signals were based on the different channel measurements from each probe, with synchronous sampling at evenly spaced angular increments. In other words, the number of sampled points per revolution remained consistent even as the rotational speed changed.In the angular domain, every sampling point corresponded to an angle of uniform distribution; therefore, the spindle radial error motion, or more specifically, the accurate SSEM, was displayed as the relative position in the angular domain. Thus, the position of the accurate SSEM corresponds to the position of the cutting tool under the milling process in the angular domain, which means the accurate SSEM corresponds to the workpiece surface roundness profile in the angular domain.


## 4.1.2. Roundness Measurement to Determine WRE

The measurement signals of the probes distributed by different angles were acquired by the system at a stable rotational spindle speed during the milling process and separated into the spindle radial error motion using the three-probe method presented in Section3.1. Then, feature extraction was performed to obtain the spindle radial error motion and SSEM.The SSEM was obtained, and the corresponding workpiece surface roundness profile was measured using a roundness measuring instrument (YD200, LHCT company). In the measuring process, the machining workpiece first needed to be placed in the appropriate position. The inner hole roundness of machining workpiece measurement is shown in Figure5.Figure 5
Roundness measurements using YD200.


## 4.2. Establishment of Predictive Model

The predictive model is based on actual measurements. Signals were acquired from a spindle working to cut the workpiece at a stable rotational speed of 4200 r/min, and accurate SSEM values is shown in Figure6(a). Measurements of the inner hole roundness of the machined workpiece is presented in Figure 6(b).Figure 6
Measurements obtained during the machining process. (a) SSEM. (b) Workpiece roundness profile.
(a)
(b)It can be seen from Figure6 that the accurate SSEM has a similar shape to the workpiece roundness profile. To verify this, a similarity analysis of the image histograms was performed and to determine the degree of similarity. The histogram, usually referred to as a fingerprint image, depicts the intrinsic properties of the image and is not affected by size, position, or rotational angle.


### 4.2.1. Correlation of Measurement Results

To investigate morphological similarities between accurate SSEM and workpiece roundness profile, the plots were extracted, converted to grayscale images (190 × 200 pixels), and histograms of the grayscale images were generated. As an example, histograms for a spindle rotational speed of 4200 r/min are shown in Figure7.Figure 7
Initial images and histograms. Grayscale image of (a) accurate SSEM and (b) workpiece roundness profile. Histogram of (c) accurate SSEM and (d) workpiece roundness profile.
(a)
(b)
(c)
(d)From Figure7, it can be seen that the distributions of pixels amongst the grayscale values are close to each other in both histograms. The degree of similarity between histograms can be judged in many ways, including distance judgment methods commonly selected for comparing image morphologies. The histogram intersection distance and the distance proposed by Bhattacharyya et al. can indicate similarity between histograms and only differ in their expression.The histogram intersection distance reflects the number of repeated pixels in the histogram, and the result can be normalized to be within the range [0, 1] by dividing the histogram intersection distance by the number of pixels in one of the histograms. Similarity between histograms is higher as the normalized result gets closer to 1. The normalized histogram intersection distanceDnhi can be calculated as(10)Dnhi=∑i=1NpminHGSi,HGRi∑i=1NpHGSi,where HGSi is the histogram of the accurate SSEM, HGRi is the histogram of workpiece roundness profile, and Np is the total number of pixels. The normalized histogram intersection distance Dnhi is 0.9989.Analogously, the Bhattacharyya distance reflects the probability distribution of the pixels of different histograms, and its normalized functionDnB can be given as(11)DnB=1−∑i=1NpHGSi⋅HGRi∑i=1NpHGSi⋅∑i=1NpHGRi.In this case, the normalized histogram intersection distanceDnB is 0.0047.According to the distance judgment methods, the similarity between the accurate SSEM and workpiece roundness profile shows that the normalized histogram intersection distance is very close to 1 and the normalized Bhattacharyya distance is close to 0. It follows that, under a steady rotational speed, both similarity judgment methods indicate that the histograms have a high degree of similarity in relation to image morphology.


### 4.2.2. Predictive Model

To establish the predictive model via numerical analysis, the relationship between the SSEM value and the WRE must be determined, as well as the relationship between data sequences of SSEM and workpiece roundness profile. The SSEM value and the WRE can be calculated from the SSEM and workpiece roundness profile, respectively. The SSEM value, the WRE, and the ratio between them are presented in Table1. A correlation analysis can be introduced to obtain the relationship between data sequences of the accurate SSEM and workpiece roundness profile. The correlation coefficient for the two correlated sequences can be defined as(12)ρesrθ,rtwθ=∑θθendesrθ−esrθ¯⋅rtwθ−rtwθ¯∑θθendesrθ−esrθ¯2⋅∑θθendrtwθ−rtwθ¯2.Table 1
Results of proportional analysis and correlation analysis.
Rotational speed (r/min)
SSEM valueResr (μm)
WRERtw (μm)
RatioKRa=Resr/Rtw
Correlation coefficientρesr,rtw
4200
5.04
5.75
0.88
0.91The correlation coefficient for the SSEM value and WRE was acquired and is presented in Table1.As seen in Table1, the SSEM value and WRE are directly proportional and the proportionality coefficient KRa is 0.88. Assuming the SSEM value is known, the proportionality coefficient Resr can be used to establish the predictive model, and WRE can be expressed as Rtw=cpResr=1/KRa⋅Resr=1.14Resr. The proportional relationship between the SSEM value and WRE is a comprehensive technical indicator but only reflects the numerical relationship in a single dimension. This oversimplified model for online WRE prediction is not accurate enough.Furthermore, the correlation coefficientρesrθ,rtwθ is 0.91, which is very large and suggests the SSEM and workpiece roundness profile are highly correlated, according to the numerical analysis. The correlation coefficient consistently reflects trends in the signal sequences of the SSEM and the workpiece roundness profile more detailed, equivalent to extending the dimensions of the analysis, and can therefore be used as a correction factor in the predictive model along with proportionality coefficient for WRE, which comprehensively reflects the performance of the tool. Based on the proportionality coefficient and correction factor, the predictive model can be established as(13)Rtw=cpcResr=1KRa⋅ρesrθ,rtwθ⋅Resr=1.25Resr.The necessary correlation between the SSEM value and WRE was found to exist, and the relationship can be defined using a precise mathematical formula, which can then be incorporated into the predictive model. At the same time, based on proportionality and the correlation analyses, the predictive model can be obtained. Both the proportionality coefficient and the correlation coefficient play an important role in the predictive model.According to the description above, when the inner surface quality of the machined workpiece is difficult to determine by direct measurement during the machining process, the WRE is the most concerned parameter in deep hole machining and can be indirectly predicted using the online prediction method. For SSEM and workpiece roundness profile measurements exhibiting high morphological similarity, a correlation analysis was performed. Based on the results, the predictive model was established and used for online prediction of the inner hole roundness error. Thus, the method can improve upon existing online monitoring techniques and enhance machining accuracy.


## 4.2.1. Correlation of Measurement Results

To investigate morphological similarities between accurate SSEM and workpiece roundness profile, the plots were extracted, converted to grayscale images (190 × 200 pixels), and histograms of the grayscale images were generated. As an example, histograms for a spindle rotational speed of 4200 r/min are shown in Figure7.Figure 7
Initial images and histograms. Grayscale image of (a) accurate SSEM and (b) workpiece roundness profile. Histogram of (c) accurate SSEM and (d) workpiece roundness profile.
(a)
(b)
(c)
(d)From Figure7, it can be seen that the distributions of pixels amongst the grayscale values are close to each other in both histograms. The degree of similarity between histograms can be judged in many ways, including distance judgment methods commonly selected for comparing image morphologies. The histogram intersection distance and the distance proposed by Bhattacharyya et al. can indicate similarity between histograms and only differ in their expression.The histogram intersection distance reflects the number of repeated pixels in the histogram, and the result can be normalized to be within the range [0, 1] by dividing the histogram intersection distance by the number of pixels in one of the histograms. Similarity between histograms is higher as the normalized result gets closer to 1. The normalized histogram intersection distanceDnhi can be calculated as(10)Dnhi=∑i=1NpminHGSi,HGRi∑i=1NpHGSi,where HGSi is the histogram of the accurate SSEM, HGRi is the histogram of workpiece roundness profile, and Np is the total number of pixels. The normalized histogram intersection distance Dnhi is 0.9989.Analogously, the Bhattacharyya distance reflects the probability distribution of the pixels of different histograms, and its normalized functionDnB can be given as(11)DnB=1−∑i=1NpHGSi⋅HGRi∑i=1NpHGSi⋅∑i=1NpHGRi.In this case, the normalized histogram intersection distanceDnB is 0.0047.According to the distance judgment methods, the similarity between the accurate SSEM and workpiece roundness profile shows that the normalized histogram intersection distance is very close to 1 and the normalized Bhattacharyya distance is close to 0. It follows that, under a steady rotational speed, both similarity judgment methods indicate that the histograms have a high degree of similarity in relation to image morphology.


## 4.2.2. Predictive Model

To establish the predictive model via numerical analysis, the relationship between the SSEM value and the WRE must be determined, as well as the relationship between data sequences of SSEM and workpiece roundness profile. The SSEM value and the WRE can be calculated from the SSEM and workpiece roundness profile, respectively. The SSEM value, the WRE, and the ratio between them are presented in Table1. A correlation analysis can be introduced to obtain the relationship between data sequences of the accurate SSEM and workpiece roundness profile. The correlation coefficient for the two correlated sequences can be defined as(12)ρesrθ,rtwθ=∑θθendesrθ−esrθ¯⋅rtwθ−rtwθ¯∑θθendesrθ−esrθ¯2⋅∑θθendrtwθ−rtwθ¯2.Table 1
Results of proportional analysis and correlation analysis.
Rotational speed (r/min)
SSEM valueResr (μm)
WRERtw (μm)
RatioKRa=Resr/Rtw
Correlation coefficientρesr,rtw
4200
5.04
5.75
0.88
0.91The correlation coefficient for the SSEM value and WRE was acquired and is presented in Table1.As seen in Table1, the SSEM value and WRE are directly proportional and the proportionality coefficient KRa is 0.88. Assuming the SSEM value is known, the proportionality coefficient Resr can be used to establish the predictive model, and WRE can be expressed as Rtw=cpResr=1/KRa⋅Resr=1.14Resr. The proportional relationship between the SSEM value and WRE is a comprehensive technical indicator but only reflects the numerical relationship in a single dimension. This oversimplified model for online WRE prediction is not accurate enough.Furthermore, the correlation coefficientρesrθ,rtwθ is 0.91, which is very large and suggests the SSEM and workpiece roundness profile are highly correlated, according to the numerical analysis. The correlation coefficient consistently reflects trends in the signal sequences of the SSEM and the workpiece roundness profile more detailed, equivalent to extending the dimensions of the analysis, and can therefore be used as a correction factor in the predictive model along with proportionality coefficient for WRE, which comprehensively reflects the performance of the tool. Based on the proportionality coefficient and correction factor, the predictive model can be established as(13)Rtw=cpcResr=1KRa⋅ρesrθ,rtwθ⋅Resr=1.25Resr.The necessary correlation between the SSEM value and WRE was found to exist, and the relationship can be defined using a precise mathematical formula, which can then be incorporated into the predictive model. At the same time, based on proportionality and the correlation analyses, the predictive model can be obtained. Both the proportionality coefficient and the correlation coefficient play an important role in the predictive model.According to the description above, when the inner surface quality of the machined workpiece is difficult to determine by direct measurement during the machining process, the WRE is the most concerned parameter in deep hole machining and can be indirectly predicted using the online prediction method. For SSEM and workpiece roundness profile measurements exhibiting high morphological similarity, a correlation analysis was performed. Based on the results, the predictive model was established and used for online prediction of the inner hole roundness error. Thus, the method can improve upon existing online monitoring techniques and enhance machining accuracy.


## 5. Results and Discussion

### 5.1. Experimental Validation of Predictive Model

To validate the inner hole roundness error online prediction method, inner holes were machined by a continuous milling process using the parameters of the model. The SSEM was measured online at the rotational speed of 4200 r/min by the experimental setup mentioned in Section4.1, and the values were obtained using the method described in Section 3.3 and are presented in Figure 8.Figure 8
Online measurements of SSEM for experimental validation of the predictive model.
(a)
(b)
(c)
(d)
(e)From Figure8, it can clearly be observed that SSEM is irregular in morphology. The calculated SSEM values and predicted results of the different models are presented in Table 2, including the model using only the proportionality coefficient as well as the corrected model.Table 2
Predictive performance analysis.
SSEM value (μm)
WRE (μm)
Naive modelRtw=cpResr
Composite modelRtw=cpcResr
Prediction (μm)
Deviation (μm)
Difference (%)
Prediction (μm)
Deviation (μm)
Difference (%)
(a)
6.40
7.70
7.30
−0.40
−5.25
8
0.30
3.90
(b)
8.93
11.45
10.18
−1.27
−11.09
11.16
−0.29
−2.51
(c)
5.14
6.52
5.86
−0.66
−10.13
6.46
−0.09
−1.46
(d)
6.15
7.19
7.01
−0.18
−2.49
7.69
0.50
6.92
(e)
7.40
9.12
8.44
−0.68
−7.50
9.25
0.13
1.43To verify exactness of the predicted results, the workpiece surface roundness profiles were tested by the roundness measuring instrument mentioned in Section4.1 and shown in Figure 9. Numerical characteristics of workpiece surface roundness profiles are listed in Table 2, and a comparative analysis of the results is also provided.Figure 9
Workpiece surface roundness profile for verifying exactness of the predicted solutions.
(a)
(b)
(c)
(d)
(e)Figure9 shows that similarity also exists between the SSEM and workpiece surface roundness profile. The comparative analysis of the numerical characteristics provides accurate and intuitive data, as shown in Table 2. Advantages of the different models were compared by extracting the absolute value of the deviation and percentage differences. In the composite model, deviations from the predicted results can be controlled within 0.5 μm (a difference of less than 7%) and the prediction accuracy is significantly improved compared to results of the naive model, with deviations reaching 1.27 μm (a difference of 11.09%). Comparing with research results of Ryua et al. [5] and Marsh et al. [6], the predicted results here not only give the qualitative analysis results but also provide the quantitative analysis results, and it has high accuracy. Thus, it can be concluded that the method proposed in this paper could be used to accurately predict WRE online.


## 5.1. Experimental Validation of Predictive Model

To validate the inner hole roundness error online prediction method, inner holes were machined by a continuous milling process using the parameters of the model. The SSEM was measured online at the rotational speed of 4200 r/min by the experimental setup mentioned in Section4.1, and the values were obtained using the method described in Section 3.3 and are presented in Figure 8.Figure 8
Online measurements of SSEM for experimental validation of the predictive model.
(a)
(b)
(c)
(d)
(e)From Figure8, it can clearly be observed that SSEM is irregular in morphology. The calculated SSEM values and predicted results of the different models are presented in Table 2, including the model using only the proportionality coefficient as well as the corrected model.Table 2
Predictive performance analysis.
SSEM value (μm)
WRE (μm)
Naive modelRtw=cpResr
Composite modelRtw=cpcResr
Prediction (μm)
Deviation (μm)
Difference (%)
Prediction (μm)
Deviation (μm)
Difference (%)
(a)
6.40
7.70
7.30
−0.40
−5.25
8
0.30
3.90
(b)
8.93
11.45
10.18
−1.27
−11.09
11.16
−0.29
−2.51
(c)
5.14
6.52
5.86
−0.66
−10.13
6.46
−0.09
−1.46
(d)
6.15
7.19
7.01
−0.18
−2.49
7.69
0.50
6.92
(e)
7.40
9.12
8.44
−0.68
−7.50
9.25
0.13
1.43To verify exactness of the predicted results, the workpiece surface roundness profiles were tested by the roundness measuring instrument mentioned in Section4.1 and shown in Figure 9. Numerical characteristics of workpiece surface roundness profiles are listed in Table 2, and a comparative analysis of the results is also provided.Figure 9
Workpiece surface roundness profile for verifying exactness of the predicted solutions.
(a)
(b)
(c)
(d)
(e)Figure9 shows that similarity also exists between the SSEM and workpiece surface roundness profile. The comparative analysis of the numerical characteristics provides accurate and intuitive data, as shown in Table 2. Advantages of the different models were compared by extracting the absolute value of the deviation and percentage differences. In the composite model, deviations from the predicted results can be controlled within 0.5 μm (a difference of less than 7%) and the prediction accuracy is significantly improved compared to results of the naive model, with deviations reaching 1.27 μm (a difference of 11.09%). Comparing with research results of Ryua et al. [5] and Marsh et al. [6], the predicted results here not only give the qualitative analysis results but also provide the quantitative analysis results, and it has high accuracy. Thus, it can be concluded that the method proposed in this paper could be used to accurately predict WRE online.


## 6. Discussion

The combination of probe angles for measuring SSEM has a large effect on the weighted functionFλ. Herein, we discuss how certain angle combinations affect the performance of Fλ as well as characteristics of the optimized angle combination proposed in this paper.The fundamental harmonic of the measurement signal (λ=1) is related to eccentric mounting (rotor eccentricity), and the effects of rotor eccentricity can be eliminated by using eccentricity compensation. The 2nd to 50th harmonics of the measurement signal (2≤λ≤50) make up the surface roundness profile. Higher frequency harmonics (λ>50) are usually presented as surface waviness or roughness. Therefore, in analyzing the surface roundness profile, a performance assessment of the weighted function Fλ needs to only consider the range 2≤λ≤50.The optimal angle actually corresponds to the value ofFλ, according to every ω kept away from zero. The measurement angle using the conventional angle combination without optimization corresponding to Fλ is shown graphically in Figure 10.Figure 10
Impact of the measurement angle combination on the weighted function ∣F(λ)∣.
(a)
(b)As seen in Figure10, the conventional measurement angle combination without optimization shows many values of Fλ satisfying the usage requirements; however, many values of the harmonic order λ lead values of Fλ quite close to zero, and zero values of Fλ lead to harmonic suppression. Due to the harmonic suppression, the results for the surface roundness profile which are separated from the measurement signal using the three‐probe method will therefore be incorrect; moreover, an accurate measurement of the spindle radial error motion cannot be obtained.For further verification, the angle combination acquired using the PSO was compared with the angle combinations proposed by Shi et al. (0°, 42.8°, 126.4°) [12] and Grejda (0°, 122.5°, 202.5°) [13]. Then, Fλ was determined by substituting the measurement angle combinations into equation (4), as illustrated in Figure 11.Figure 11
Comparison of the performance of different optimization measurement angle combinations.
(a)
(b)
(c)The results suggest that by using the optimized angle combination (0°, 45°, 96.5°), the minimum value of the weighted functionFλ is furthest from zero, as shown in Figure 11. Moreover, the entire weighted function appears to be further away from zero compared to functions associated with other angle combinations. Additional analyses of the results are necessary to verify this. The performance of different angle combinations were assessed based on several indexes including the minimum value Fmin=minFλ, mean value μ, difference between the mean value and standard deviation δ=μ−sd, and the number of values close to zero ncz (Fλ<0.2), as shown in Table 3.Table 3
Influence of different measurement angles on various indexes.
Measurement angle
F
min
μ
δ
=
μ
−
s
d
n
cz
(0°, 30°, 90°)
0
1.39
0.67
8
(0°, 30°, 120°)
0
1.24
0.76
8
(0°, 30°, 150°)
0
1.14
0.77
8
(0°, 60°, 90°)
0
2.41
0.26
8
(0°, 60°, 120°)
0
1.35
0.12
16
(0°, 60°, 150°)
0
1.23
0.68
8
(0°, 42.8°, 126.4°)
0.17
1.34
0.93
1
(0°, 122.5°, 202.5°)
0.27
1.27
0.96
0
(0°, 45°, 96.5°)
0.74
1.71
1.21
0The conventional measurement angle combinations tend to be close to zero at the minimum valueFmin, which results in harmonic suppression. In comparison, the optimized angle combinations push Fmin further away from zero, and as the minimum value becomes larger, the harmonics become unsuppressed. Analogously, the number of values close to zero should be as small as possible and is an important index to invalidate the error separation method. As an initial criterion, each individual index can be used to screen for the most suitable angle combinations.In contrast to the individual index, for situations in which the minimum value is not close to zero, however, the individual weighted value is close to the minimum, i.e., the entire weighted function is close to zero, so there is a higher probability of harmonic suppression. Therefore, it is necessary to introduce global indexes to identify the angle combinations with superior performance. The mean valueμ and the difference between the mean value and standard deviation δ are global indexes that can be used to ensure the entire weighted function is kept far from zero.Of particular importance is that there is some interaction between individual and global indexes, which means that for largerμ values, Fmin may be equal to zero and ncz will be large. Examples of conventional measurement angle combinations are shown in Figure 10; however, these angle combinations are not suitable for high-accuracy measurement of the spindle radial error motion. This highlights the importance of the global index δ since the standard deviation represents volatility of weighted function and δ can be used to determine whether the entire weighted function with eliminated volatility is kept far away from zero. Taking the angle combination of (0°, 60°, 90°) as an example, μ is larger compared to the values of other combinations, but δ is smaller. This suggests the value of Fλ has a large amount of volatility, the individual index Fmin is equal zero, and ncz is equal to 8. Thus, the angle combination is not suitable. Both the individual and global indexes should comprehensively be considered when screening for the most appropriate angle combination.Based on the results, the measurement angle combination of (0°, 45°, 96.5°), which was optimized using PSO, offers superior performance, as shown in Table3 and Figure 11. Therefore, the optimal angle of (0°, 45°, 96.5°) should be used to achieve high accuracy during the milling process and to investigate the spindle radial error motion of a high-speed spindle.


## 7. Conclusion

This paper proposed a novel method for online prediction of WRE based on the SSEM value, and mainly in better performance by comparison with previous studies. Results of the analysis show the method is feasible and practical. The main findings can be summarized as follows:(1)
The manufacturing accuracy of deep inner hole structure as the hydraulic cylinder is restricted to online measurement of machining quality; to solve the problem, the method based on accurate SSEM value extraction is proposed to predicted WRE during machining process, the concrete implementation processes, and technical essentials of this method are introduced in this paper, and this method has important academic meaning and value of practical application on high-accuracy manufacturing of deep inner hole.(2)
During high-speed hole milling, the accurate spindle radial error motion can be measured in real time using the three-probe method. The influence of surface roundness error was eliminated thereby improving measurement accuracy. The PSO algorithm was used to optimize the mounting angles of the probes according to the set of optimization goals. The optimized angle combination (0°, 45°, 96.5°) prevented higher-order harmonic suppression and enhanced measurement accuracy.(3)
The SSEM can be extracted from the spindle radial error motion with high accuracy using the harmonic wavelet even in the presence of strong background noise. Feature signal extraction of integer multiples of the rotational frequency provides the necessary analytical foundation for subsequent accurate calculations of the SSEM value.(4)
The predictive model established using the proportionality coefficient and correlation coefficient has high prediction accuracy, and deviations from true values can be controlled to be within an extremely small range.Hence, the novel method could be used for online monitoring of machining quality during high-precision hole processing. Furthermore, the machining state of the machine tool can be adjusted online based on the predicted results, the control decision reference for dynamic balance can be provided online, and processing accidents can be avoided. Thus, the method offers the potential to drive high processing efficiency, advanced processing technologies, and good processing quality.



---

*Source: 6049316-2019-08-05.xml*