Diagnosis of two-axis cylindrical grinder machining error

This paper attempts to develop a machining error diagnosis method for classifying the error sources on the external/internal contour surfaces which are fabricated using a two-axis cylindrical grinding machine. A derivation of the contour profile to grinding wheel motion trajectory and also the grinding wheel motion trajectory to the contour profile equations has been performed. Due to the contours having a sinusoidal polynomial polar curve instead of a rose curve, the grinding wheel motion trajectory equation is then used to generate the NC code for the machining. The grinding wheel motion trajectory equation is needed for the synthesis of the profile error caused by the motion error as well as the variation of the grinding wheel diameter. The two-axis motion equation is obtained by transforming the non-circular contour equation through the coordinates of the grinding wheel and workpiece. The influence of different machining parameters such as grinding wheel outer diameter error, grinding wheel feed rate, and workpiece rotation speed on the external/internal contour profile is being studied and analyzed. The diagnosis method is based on the least square method to assemble the components of the profile error due to each individual error sources. The actual machining with probe system is used to verify the diagnosis result.


Introduction
A continuous development of modern manufacturing technology and the improvement of product quality requirements has led the mechanical manufacturing industry rapidly proceed to the direction of high precision, high speed, and high efficiency. An important and crucial performance factor being geometric accuracy for machine products, especially under circumstances where relatively high precision is one of the basic requirements. Afterwards, an error compensation plays a significant role in ensuring the overall product quality [1]. Error diagnosis and compensation has three engineering steps: kinematic modeling, error diagnosis, and error compensation which is based on certain dominating error factors. Error diagnosis plays a vital role for grinding machines, as the final accuracy of any part is directly related to the grinding machine processing accuracy. Its improvement is critical for machine tool design research. The authors in this paper tried to shed some light in analyzing the different depending components for the grinding error and also provide a definitive compensation method for each.
The kinematics of the machine have been modeled by applying the homogeneous transformation matrix. The contour profile relating to the grinding wheel motion trajectory and also the grinding wheel motion relating to the contour profile has been extensively derived initially. To improve the precision of machined products, geometrical error diagnosis is generally employed by different researchers which further paves the way for several compensation techniques. Wang et al. developed a simulation software depicting the influence of X-axis feed error, the C-axis feed error, and the wheel radius error on the grinding process [2], which ultimately provides the basis for error compensation. In order to increase grinding accuracy, Tao and Tao revealed a coupling relationship between the thermal and geometrical error [3] which laid a foundation for subsequent error compensation. An On-Line Asynchronous Compensation Method (OACM) is proposed by Shen et al. for static/quasi-static error caused by thermal deformation and machine geometry [4]. The method reduces the complexity of compensation system while an offset compensation method (OCM) adopts Bayesian network to predict error. Tian et al. presents a general and systematic approach for geometric error modeling of machine tools which provides a direct linear mapping from each source error to pose accuracy of machine tool and also has clear physical meanings [5]. However, motion or structural errors of spindle, cutting tools, and worktable are not taken into consideration. Okafor and Ertekin [6] developed mathematical model to calculate and predict the resultant error vector at the tool workpiece interface for a three-axis vertical machining center. Xiang and Altintas proposed a method [7] to measure, model, and compensate both geometrically dependent and independent volumetric errors of 5-axis, serial CNC machine tools, where a total a 41 errors of all five axes are represented by error motion twists. Several authors have proposed an integrated geometric error prediction and compensation method [8,9]. Zuo et al. presents the model [10] which is modeled by the propagation and the accumulation of errors based on Jacobian-Torsor theory. Intelligent AI techniques has also been introduced for the improvement of geometrical errors in grinding process [11,12].
This research article deals with an example of a non-circular grinding to validate its results. A non-circular grinding is very different from conventional grinding methods. Onishi et al. analyzes the elastic deformation for a long workpiece in a cylindrical grinding process [13]. A noncircular grinding suffers a disparity between the differential geometry and computational geometry. An undercut is a problem that is very related to computational geometry. A non-circular grinding geometric error is also a major cause of the undercut. Jang and Choi [14] studied the grinding force, considering the characteristics of non-circular grinding. A compensation model was thus established that can estimate the set depth of cut to further obtain the real depth of cut. A major variable considered for geometrical error also adheres to an extensive tool condition monitoring. Arun et al. provides a tool condition monitoring approach using an acoustic emission sensor [15]. Various compensation methods are also established in terms of CNC machines regarding dominant errors [16]. Liu et al. had concluded a problem in the grinding wheel which needed to be frequently dressed due to inevitable grinding wear. A novel method has thus been proposed and employed to compensate the profile error caused by grinding wheel wear and improve grinding efficiency [17]. A recent article by Fa and Ye [18] proposed a modeling method in contrast to the multi-body system, which can simplify the calculation process and reflect the influence of each error components on the tool. A comprehensive error was determined by combining 24 fitting functions, the known error, the differential matrix of each part, and the Jacobian matrix.
In this article, the authors attempt to develop a diagnosis method for classifying the machining error sources on the external/internal contour surfaces which are fabricated using a two-axis cylindrical grinding machine. Differential geometry is used for a thorough derivation between the grinding wheel position to contour and vice-versa. A disparity between the concepts of computational geometry and differential geometry which can produce an undercut and ways to reduce it has been explained. A diagnosis of geometric error, linear and non-linear motion error components is shown according to the control system for both axes in the CNC system that is required to control the workpiece rotation and movement of the grinding wheel. A compensation method has thus been established.
This paper is organized as follows. Section 2 describes the cylindrical grinding machine and its machining process. Section 3 describes the differential geometry required for the surface profile in terms of the grinding wheel and vice-versa. Section 4 studies the disparity in the concepts of differential geometry to computational geometry leading to undercut and establish the grinding wheel constraints. Section 5 provides a detailed description for the diagnosis of the machining error and its error shape function. A least square curve fitting has been established. Section 6 describes the error compensation method based on the diagnosis. Section 7 substantiates the practicality of the diagnosis through experiments involving internal and external grinding of a workpiece. Section 8 concludes the research.

Cylindrical grinding and its machining process
The concept of cylindrical grinding [19][20][21] refers to the grinding of bores and holes. Internal/External grinders are internal/external grinding tools or machine tools used for cylindrical profiles, which uses grinding wheels as cutting tools. One possible arrangement for the two-axis cylindrical grinding machine may be arrangement as shown in Fig. 1, that the workpiece is fixed to a rotary axis with rotation index and the grinding wheel is fixed to a spindle. The grinding wheel spindle is horizontally movable with the grinding wheel into the grinding position.
Despite the actual grinding machine with two axes independently controlled from the ground table. We express the grinding process in a kinematic inversion way that makes the grinding wheel to move in a plane which representing a specific cross-section of the workpiece. The planar motion of the grinding wheel is specified in the polar coordinate using the radius r and radian respectively. The grinding wheel remove the material from the workpiece to form the inner surface profile. It is also assumed that the slurry and fluid used in grinding process do not affect the profile contour.

Parametric contour of non-circular surface cross-section profile
The cavity of the workpiece is generally circular shape with variation described by two concentric circles centered at the origin O. One is the circumcircle of the cavity with diameter D a and the other is the inscribed circle of the cavity with diameter D i . A nominal circle with diameter D m which takes the mean of the two circles as follows.
A 2D parametric curve Q( ) bounded by the circumcircle and inscribed circle is used to describe the inner surface of the cavity. The Euclidean distance between Q( ) and the origin O denoted by ‖Q( )‖ is confined by the relation as follows. (1) There is at least one point on Q( ) touches the circumcircle and the other point on Q( ) touches the inscribed circle. The peak-to-peak amplitude of Q( ) is denoted by 2e = D a − D i . The properties of the 2D, C ∞ , sinusoidal parametric curve Q( ) where min ≤ ≤ max include.
1. G∞ continuity: The joint points Q( min ) and Q( max ) sat- A specific parametric curve Q( ) with min = 0 , max = 2 , and N = 3, expanded using the product-to-sum formula is proposed as follows.
Five properties are obtained from the proposed parametric curve (Fig. 2).
1. Since all terms are sinusoidal functions, the curve is C ∞ continuous. 2. Since all terms are harmonics function of integer multiples of a fundamental frequency , Q (n) (0) = Q (n) (2 ) for n = 0, 1, 2, … ∞. 3. Since   5. The extreme radii of the parametric contour from the stationary condition as follows.
To satisfy Eq. (2) for the contour design, the above equation imposes a constraint for coefficient b to avoid the extreme radii contour point is larger thana , i.e.  Figure 3 shows different parametric contours generated from different b/e ratio according to Eq. (3). The result shows that the larger the b/e ratio the sharper the corner is. Due to the a/e must be no less than 1 according to Eq. (5), it is then derived that when b/e ratio is no less than 0.414 the result in Eq. (5) is guaranteed. On the other when we set b/e ratio to be 2, then the lowest a/e is 8.

Contour profile measurement
An articulated coordinate measuring machine (CMM) can be used to measure the geometry of physical objects by sensing discrete points on the surface of the cavity with a probe. CMM specifies a probe's position in terms of its displacement from a reference position in a two-dimensional Cartesian coordinate system with XY axes. The sensed discrete taken from the probe written as follows are used to compare with the desired profile. The planar angle of measurement is defined as i . The corresponding parameter ( i ) is obtained from solving the desired profile such as stated in Eq. (4) iteratively as follows.
After ( i ) is obtained from a specific sensed discrete points M i , the error between the desired inner surface profile Q( ( i )) can then by compared with the specific sensed discrete points M i as follows.
The magnified plot, as shown in Fig. 4, of a specific measurement may be drawn to magnify the error at the specific polar location is performed using E i = Θ i + Q i , where Θ denotes the magnification factor. An error plot example from the CMM measurement is shown in Fig. 4.

Grinding wheel motion to contour disregarding undercut
For the convenience of analysis, the grinding process with the grinding wheel position control is easier to be shown in the kinematic equivalence scheme as shown in Fig. 5a. The point P on the grinding wheel is derived as follows. where The differential geometry equation for P being the profile point is derived as follows.
The above equation can be subsequently derived by substituting the vector derivatives. The contact point of the grinding wheel is derived as follows. where, The angle is the pressure angle with respect to the line ⃖���� ⃗ OC which connects the grinding wheel center and workpiece center. The positive sign on both the numerator and denominator of is applicable to the internal grinding and the negative sign is used in the external grinding. Using , it can be seen in the internal grinding case as shown in Fig. 5b that the contour point P of the grinding wheel determined by is the same as the grinding wheel center orientation angle at the position r∕ = 0. The angle < when r∕ is positive and vice versa. The pressure angle of the grinding machining is defined differently from the gear meshing definition. In mechanical gear pair, the pressure angle at a pitch point between the line of pressure (which is normal to the tooth surface) and the plane tangent to the pitch surface. The pressure angle in the grinding process is a complementary angle of the pressure angle defined in power transmission using gears. High pressure angles are not preferable in both definitions due to different purposes that one is power transmission which want to reduce the abrasive force and the other one is material removal which want to increase the abrasive force between two rigid bodies.
It may be given with an arbitrary trajectory as follow.
The polar coordinate for the curve is defined using r and . The derivative of r with respect to is defined using the derivative chain rule as follows. Thus The profile generated by the wheel to form the profile Q at a specific grinding wheel position C( ) ignoring the undercut is derived as follows.
We obtain ( ) from the procedures similar to in Sect. 3.1 as follows.

Contour to grinding wheel position
For a given parametric contour Q( ) for the internal surface profile, the grinding wheel position can be derived from offsetting the contour with the grinding wheel radius.
N( ) is the unit vector normal to the unit tangential vector T of contour at contour position Q( ) , which can be obtained from the first derivative of T , while the unit tangential vector can be obtained from the first derivative of Q . The first derivative of Q is derived from Eq. (3). The curvature ( ) of the parametric contour Q( ) is derived as follows.
The degeneracy of the curvature ( ) occurs when Q � ( ) becomes a null vector having zero length. By re-examining Eq. (3) for the proposed contour, we find the degeneracy will start with the positions having largest curvature locations at = (2n+1) 3 , n = 0, 1, 2.
Therefore, the design constraint for the parametric contour must follow the constraint equation to avoid the curvature degeneracy, i.e., a e > 2 + 3 b e . It is then observed that the a/e ratio must be larger than 2. The upper bound of b/e ratio in Sect. 3.1 yields a dominant bound when a/e ratio is less than 8, while the recent upper bound of b/e ratio becomes dominant when a/e ratio is higher than 8 (Fig. 6).

Grinding wheel radius constraint to avoid the undercut in internal grinding
The center of curvature ( ) of the contour Q( ) , known as the fixed centrode of the grinding wheel motion, can be written as follows.
As shown in Fig. 7, the minimum distance between ( ) and the corresponding Q( ) of our propose contour is located at = (2n+1) The above equation yields an upper limit for the radius of the grinding wheel chosen to be a > min(3b + 2e, 6b − 4e) . That is if the radius of the grinding wheel is greater than the minimum radius of curvature then there is an undercut problem.
When the undercut problem occurs the profile of the inner surface undergoes an interference between the contour and the grinding wheel movement. The position control of the grinding wheel is subjected to a non-monotonic increasing function to rotate the workpiece as shown in Fig. 7. The interference suffers an invalidity of the differential geometry relation in Eq. (10) to form an actual profile.

Case study
For a particular arrangement, we have b/e ratio set to be 2, where a > min(3b + 2e, 6b − 4e) can be reduced into an expression as follows.
The above equation can then be written subsequently in terms of a/e ratio as follows.

Undercut
Enlarged view of fixed centrode locus According to the above equation, the maximum grinding wheel diameter to the internal profile nominal diameter ratio R w ∕a can be determined. Depending on the choice of a/e ratio, the maximum grinding wheel diameter can be selected to minimize the rest time when the diameter of the grinding wheel is known keeping reduced during abrasive grinding process. As shown in Fig. 8, when the a/e ratio is above 16 then the maximum grinding wheel diameter can be higher than half of the internal profile nominal diameter.

Diagnosis of machining error
Grinding machining is the process of removing metal from a workpiece in the form of tiny chips by the action of irregularly shaped abrasive particles. Grinding wheels are composed of thousands of small abrasive grains [22] held together by a bonding material. Each abrasive grain is a cutting edge. As the grain passes over the workpiece, it cuts a small chip, leaving a smooth surface. As each abrasive grain becomes dull, it breaks away from the bonding material. Diameter of the grinding wheel will keep on reducing during the machining process, therefore the NC codes must be updated to compensate the diameter change in order to maintain a higher accuracy on the workpiece dimension. Measurement of these grinding wheels often relies on tactile devices in a closed-loop process, which alternates between grinding and measuring in iterative stages to achieve diameter tolerances of ± 1.5 µm.
The servo system for both axes in the CNC system to control the workpiece rotation and the movement of the (23) grinding wheel can involve the servo error according to the servo control system as shown in Fig. 9. The servo system consists of the position error gain K p , the velocity error gain K v , the velocity feedforward gain K FF , the velocity error integral gain K v , the speed limit, and the torque limit. The coordination between two axes are assumed to be managed during the path planning process before the position command cmd is generated. There are other factors such as the anti-wind up, the torque disturbance due to G D , the back EMF from the velocity G BEMF are assumed to be managed in the electronics system such as power amplifier and current sensor feedback.

Error sources modeling
We employed several error sources that cause the grinding machining error into the analysis. These error sources are encoded into either the grinding wheel radius or the motion error of C( ) . The error contours are enumerated based on Eq. (14), which are then compared with the desired contour to obtain the error shape functions. It is first required to obtain a template for the desired contour Q( ) from Eq. (3). The discrete contour points in terms of ( ) are obtained from Eq. (2). The error shape function is the individual error function due to a very small amount of error source, which is obtained as follows.
( ) is obtained from comparing the contour point angle ( ) obtained from with the contour point on Q (̃ ) to reach the same contour point angle ( ) . The error shape function g( ) is then in terms of the parameter (Fig. 10). According to Eqs. (11) and (13), as the angle is also a function of C( ) , the machining error due to the grinding wheel outer diameter error is not a perfect circle. Due to Q( ) is changing with respect to ( ) , the corresponding error shape function shown as in Fig. 12b approaches the original shape of Q( ).
r-axis Servo Velocity Gain Adjustment ΔK v,r : Variation of the servo velocity gain of the r-axis control affects the position control response of grinding wheel radial motion in the simplified expression as follows.
The term d d |C( )| can also affect the profile contour derived in Eq. (13). The dot term in the numerator of Eq. (13) does not change much due to this small deviation but the cross-product term on the denominator does. As a result, the r∕ term in Eq. (13) affects to the angle in Eq. (11). Hence for those r∕ = 0 regions, the servo velocity gain cannot affect the contour error much when it affect the contour accuracies on the significant rate of grinding wheel radial position change r∕ zones. The corresponding error shape function shown as in Fig. 12c is compared to a zero value circle since it has both positive and negative signs due to the sign change of r∕ in different zones.
r-axis Speed Limit ΔK v_limit,r : The speed limit of the raxis control affects the position control of grinding wheel radial motion in the simplified expression as follows.
The power z is used to represent the nonlinearity of the velocity, which is 3 in our study. The corresponding error shape function is calculated as shown in Fig. 12d.
r-axis Servo Acceleration Gain Adjustment ΔK a,r : Variation of the servo torque gain of the r-axis control affects the position control response of grinding wheel radial motion in the simplified expression as follows. The corresponding error shape function is calculated as shown in Fig. 12e.
r-axis Torque Limit ΔK a_limit,r : The torque limit of the r-axis control affects the position control of grinding wheel radial motion in the simplified expression as follows. The corresponding error shape function is calculated as shown in Fig. 12f.
θ-axis Servo Velocity Gain Adjustment ΔK v, : Variation of the servo velocity gain of the -axis control affects the rotation control response of workpiece rotational motion in the simplified expression as follows. The corresponding error shape function is calculated as shown in Fig. 12g.
θ-axis Speed Limit ΔK v_limit, : The speed limit of theaxis control affects the rotation control response of workpiece rotational motion in the simplified expression as follows. The corresponding error shape function is calculated as shown in Fig. 12h.
θ-axis Servo Acceleration Gain Adjustment ΔK a, : Variation of the servo torque gain of the -axis control affects the rotation control response of workpiece rotational motion in the simplified expression as follows. The corresponding error shape function is calculated as shown in Fig. 12i.
θ-axis Torque Limit ΔK a_limit, : The torque limit of the -axis control affects the rotation control response of workpiece rotational motion in the simplified expression as follows. The corresponding error shape function is calculated as shown in Fig. 12j.

Least square fitting
The error shape functions g i ( ) , I = 1, 2, …, n, can be considered as the Fourier expansion [23] to the first order of the error source i effect as follow.
The matrix spanning the error space for m different contour point polar angles ′ j s can be presented in form as follows.
The small error quantities used to construct the error shape functions is expressed in a diagonal matrix as follows.
The actual shape function matrix is J E . The contour error vector e can be formed from the CMM measurements which are then interpolated into the discrete points e( j ) , using (J E) = e , where is the scaled error quantity. The least square method is applicable to obtain the estimated quantity E of each individual error sources using the equation as follows. where The total error is then calculated as a linear combination of all errors as follows.
The above equation is actually the linear combination of all error shapes which can cause linear approximation error for the nonlinear system.

Compensation
Three categories of machining error sources have been included in this paper, which are the geometric error, the linear servo motion error, and the nonlinear servo motion error. The geometric error category includes the grinding wheel radius error. The linear motion error category includes the position gain error, the velocity gain errors of two axes, and the torque gain errors of two axes. The nonlinear motion error category includes the speed limit and torque limit of two axes. The compensation of each individual machine error sources are.
1. Grinding wheel radius error: The measurement of grinding wheel radius must be performed constantly to maintain the accuracy of the machining result. The grinding wheel radius error detected in the diagnosis process can be used to calculate the grinding wheel path C( ) according to Eqs. (16) and (17). 2. The position gain of grinding wheel position control error: The position gain may be either adjusted the servo control system or compensated in the pre-filtering compensation which offsets the contour Q( ) inward or outward toward the center of the workpiece according to the position gain compensation required. 3. The velocity gain of grinding wheel control error: To reduce the machining error of this kind may be achieved by either adjusting the servo control system or reducing the speed of the grinding motion commands. 4. The torque gain of grinding wheel control error: To reduce the machining error of this kind may be achieved by either adjusting the servo control system or reducing the speed override as well as the feed rate of the grinding motion commands. 5. The nonlinear motion error: Both the speed and torque limits affect the machining accuracy in a more unpredictable way. The closed loop control system with higher gain can result the nonlinearities happening. However, reducing the control gains may suffer the reduction of trajectory following accuracy simultaneously. Hence the high level motion planning including the NURB interpolation of the motion path with better speed and acceleration planning is necessary.
The motion planning to minimize the velocity variation as well as the acceleration is achieved by performing the NC code re-interpolation. There are two types of velocity minimization which can be done; minimize the -axis (may also referred to as the C-axis on the CNC machine) velocity while the other is to minimize the contact point velocity. The original was interpolated at the beginning using a linear function. Let v( (u)) represent the velocity of concern, which is obtained from (u) interpolation. The velocity minimization algorithm can be stated in the following four iteration steps.
Step 1. Set (u) = u;u = [0, 2 ] Step 2. Calculate v( (u)).  Fig. 13b that, both the velocity of C and X (or r-axis) axes are minimized due to minimizing C-axis velocity based on the algorithm stated above, which shows a 25% reduction in both maximum velocities and accelerations. When minimization of the grinding contact point velocity is concerned, the result shown in Fig. 13c demonstrates a smoother C-axis acceleration (or -axis). The smoother acceleration implies a smaller jerk which is obtained from the corresponding optimization. The differences of (u) for different interpolations is shown in Fig. 13d.

Simulation and experiment
To characterize the whole experimental process, the quantitative parameters are hereby mentioned: The workpiece material is a SUJ2 with tempered hardness 58/62 (hardness Rockwell C). The grinding wheel was an 80-grit aluminum oxide grinding wheel. For the internal grinding with the workpiece radius 44 mm with a grinding wheel radius around 20 mm, the workpiece spindle rotated on around 30 rpm, the grinding wheel speed is around 3000 rpm, and the axial depth cut is usually 0.3 mm/min during fine grinding. The dressing tool is a single-point diamond dresser. The cooling lubricant is water-soluble cutting fluid. Table 1 shows the design parameters for both external and internal grinding. The software for the machining error diagnosis is implemented using VC + + . In the simulation session, we provided error to the servo system as stated in Eqs. (26) to (35) as well as the grinding wheel radius error. In this study, we ignored the asymmetric machining errors and only include the symmetric ones including the grinding wheel radius error, the servo position gains, the servo acceleration gains, and the torque limits in the analysis due to the fact that the asymmetric error sources are insignificant in our experiments. In the verification stage, it was found that the position gain error of grinding wheel position control may be linearly dependent with the diameter error of the grinding wheel when the position gain error is very small, therefore the position gain effect shall also be ignored from the diagnosis. Table 2 shows the assigned value to the individual error sources and corresponding result obtained from the diagnosis. The diagnosis result shows an overall within 10% error. Figure 14 depicts the simulation result. The red line contour is the desired contour. The purple line contour is the magnified contour error between the contour subjected to the error sources and the desired contour. The contour error is in the order of 100um error. The blue line contour is the least square fitting result. The least square fitting error in blue line closely matches the purple line contour. The remaining difference between the least square fitting error and the contour error is again magnified using green line, which is referred to as sum of the unclassified error and the linear fitting error. The unclassified error shall be zero since all error sources are predetermined from the assigned values of error source. Therefore, the remaining difference is simply the linear interpolation error due to the error shape function which is derived from the linearization stated in Eq. (36). The linear interpolation error in 1-µm scale shows a maximum of 5 µm which is about 3% of the maximum contour error which is 150 µm. The linear approximation error shows five-times frequency compared to the frequency of the design contour, which is an important feature to distinguish the linear approximation error from the unclassified error, if any found in the experiment.
The experiment uses a Studer® cylindrical grinding machine fabricated by Fritz Studer AG. The probe system is VAST® Stylus system equipment from Carl Zeiss Co. The internal grinding uses the grinding wheel in blue color as shown in Fig. 15a. The measurement of grinding wheel is done by caliper. We feed the diameter data of grinding wheel 17 mm taken from the catalog to the controller. The tool holder is fabricated by feeding the profile data according to Table 1 and contour Eq. (3) by having a = D M ∕2 . The internal profile of the machining part was measured using the probe system as shown in Fig. 15b. The CMM data was fed into the diagnosis software and the diagnosis result is shown in Fig. 16. A grinding wheel being nearly 100 µm smaller than the new grinding wheel was purposely chosen. The diagnosis result shows that not only large servo acceleration gains but also the speed and torque limits were presented. The caliper was used before the internal grinding process and the measurement shows the diameter was 100 µm shorter than the specification which is inconsistent with the diagnosis result. The remaining error of diagnosis is 20 µm, which is once again in the form of linear approximation error. The error shape function of the speed limit and torque limit was formed from ΔK a_limit,r = 1E-11 and ΔK a_limit, = 1E-12 respectively. The results are actually 50 to 100 times of the source of the error shape function. The error sources were converted into the actual resulting machining error in um scale as shown in Table 3.
In experiment 2, we changed to a new grinding wheel and tuned the motion trajectory into the version with a  smoother workpiece speed. The second tool holder was fabricated using the new grinding wheel. The profile was again measured using the probe system as shown in Fig. 15b. The CMM data was fed into the diagnosis software and the diagnosis result is shown in Fig. 17. The diameter of the grinding wheel is new. The diagnosis result shows that the servo acceleration gain was smaller due to the said motion planning. The remaining error of diagnosis is 10 µm which is better than experiment 1; however, not acceptable as a final product whose specification is 2 µm (Table 4). A large amount of error was still present (6.56 and 8.84 µm) due to the motion planning shown in experiment 2. We thus apply the (u) mapping from the velocity minimization on the contact point as shown in Fig. 13 c and d to the following experiment. In the third experiment, the contour of the fabricated tool holder was measured using the same probe system as stated before. The CMM data was fed into the diagnosis software and the diagnosis result is shown in Fig. 18. The diameter of the grinding wheel is 1.48 µm larger than the nominal value of the catalog. The remaining error of diagnosis is 1 µm which is in the form of white noise. Therefore, it concludes that the new grinding wheel is of a larger diameter. From Fig. 18, we observed that there is 1 µm stick motion error on r-axis axis which results into six peaks marked in circular red. The backlash or other mechanical errors including stick motion, axis alignment and many others [24] may need to be examined for higher precision ( Table 5).
As a result, the tool holder is made, which satisfies the specification 2um as the final product. Following the similar line of reasoning the tool axis is also made to satisfy the same specification in order to mate with the tool holder. The two individual test parts are shown in Fig. 19.

Conclusion
The inaccuracy of the grinding machining can cause the contour of the male and female workpieces to not match each other. The precision measurement of the diameter of the grinding wheel is the key factor for the precision machining. However, the grinding wheel is abrasive and the diameter of which is constantly reducing. The in-situ measurement of the workpiece and the diagnosis software together can help the precision machining during the light-duty grinding. The diagnosis software is used to guide the operator to adjust the control input such as the diameter of grinding wheel information and speed override as well. The compensation means are machining tool and precision dependent. Due to a high cost of the grinding wheel, the diagnosis software can help the manufacturer to regenerate the new motion NC codes by monitoring radius change of the grinding wheel. The errors due to the servo gain and the motor torque limits can also be detected through the same diagnosis process, which can yield certain advices to the machine operators to make high quality products when the production throughput is not reduced. The asymmetric error sources are important when the precision requirement is higher.
The extension work of this study including the motion error of the machine and the compensation of the profile grinding error using the pre-filtering techniques will be addressed in the future.