Next Article in Journal
Multiscale Modelling and Analysis for Design and Development of a High-Precision Aerostatic Bearing Slideway and Its Digital Twin
Previous Article in Journal
An Experimental Test Bench for Cable-Driven Transmission
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Efficient Calculation Method for Sensitive PDGEs of Five-Axis Reconfigurable Machine Tool

1
College of Mechanical Engineering, Yangzhou University, Yangzhou 225127, China
2
Jiangsu Key Laboratory of Digital Manufacturing for Industrial Equipment and Control Technology, Nanjing 211816, China
*
Author to whom correspondence should be addressed.
Submission received: 28 March 2021 / Revised: 16 April 2021 / Accepted: 23 April 2021 / Published: 25 April 2021
(This article belongs to the Section Advanced Manufacturing)

Abstract

:
Sensitive geometric errors of a machine tool have significant influence on machining accuracy, and it is important to identify them. Complex modeling and analysis must be carried out to identify the sensitive geometric errors of a five-axis machine tool by using the traditional method. Once the configuration structure of the machine tools is reconstructed, repetitive error modeling and analysis are required, and the identification cycle of sensitive geometric errors is long. Therefore, this paper proposes a high-efficient calculation method for sensitive position-dependent geometric error (PDGEs) identification of a five-axis reconfigurable machine tool. According to the results of sensitive geometric errors of the RTTTR-type and TTTRR-type five-axis machine tools, the mapping expressions between sensitive PDGEs and the configuration structure of machine tools was established. In this method, sensitive PDGEs can be calculated directly according to the mapping expression, which eliminates the process of error modeling and analysis. Taking a RTTTR-type five-axis machine tool as an example, the sensitive PDGEs were calculated according to the presented mapping expressions without error modeling and analysis. A series of analysis points in the machining area were selected to compare the machining errors before and after sensitive PDGE compensation. The results show that this calculation method is accurate.

1. Introduction

Five-axis computer numeric controlled (CNC) machine tools play a leading role in the complex surface manufacturing of high-end parts, and precision is one of the key indexes to evaluate its performance. The main errors that affect the machining accuracy of machine tools include geometric errors, thermal errors, and cutting force errors. In particular, geometric errors account for more than 50% of the total errors [1], which have considerable influence on the machining accuracy, and therefore, deserve more attention. It is of great significance to recognize and control the sensitive geometric errors through precision design and error compensation.
The reconfigurable machine tool was first proposed by Professor Koren in the University of Michigan [2]. By changing the structure module of the machine tool, the transformation and replacement of products and machine tools were realized. In order to adapt to the challenge of economic benefit and processing efficiency brought by the change of processing demand, customized and flexible services of reconfigurable machine tools are provided. Many scholars have carried out a lot of work on the design of reconfigurable machine tools, optimization of configuration, and reconfigurable research of the numerical control system [3,4]. The machine tool can be automatically reconfigured to a different structure by adding, deleting, or changing the machine tool modules when the processing requirements change [5].
After the reconstruction of a five-axis machine tool, the same parts are reorganized into different structures. Change in action direction of the sensitive geometric errors leads to a change in machining accuracy of a machine tool. It is necessary to analyze the sensitive geometric errors of the reconstructed machine tool for high-efficient precision design and error compensation. If the sensitive geometric errors can be calculated quickly according to the machine tool structure, design and optimization of the structure of a reconfigurable machine tool can be conducted to avoid or to reduce the influence of sensitive geometric errors. In this way, the cost for accuracy improvement of the reconstructed machine tool can be significantly reduced. In the existing literature, complex error modeling and analysis are usually implemented for each machine tool structure to identify the sensitive geometric errors, which are unable to satisfy the rapid response requirements. Therefore, an efficient algorithm is urgently needed to identify sensitive geometric errors directly according to the configuration of a machine tool.
Sensitivity error analysis of machine tools can be used to determine the most critical errors on machine accuracy. The sensitivity error analysis methods of machine tools are usually divided into local sensitivity analysis (LSA) and global sensitivity analysis (GSA) [6]. Both of them are based on spatial geometric error modeling. LSA emphasizes the impact of a single parameter and can only judge the sensitivity of one processing point at a time. Fan et al. established the sensitivity model of a machine tool based on the first-order sensitivity analysis method [7]. Li et al. used the LSA method to analyze the error sensitivity of the spatial accuracy [8]. Chen et al. established an error model of an RTTTR-type five-axis machine tool based on the MBS method. The error values were uniformly set as determinate values, and the LSA method was used to analyze the sensitivity of the error [9]. CHENG et al. established the sensitivity matrix of a four-axis precision machine tool based on the matrix differential method and used the LSA method to identify the key geometric errors of the machine tool [10]. Yao et al. analyzed the local sensitivity of the errors based on the multiple linear regression method [11].
GSA considers the coupling relationship, which can determine the influence intensity of each error term and strength of interaction among various error terms. Cheng et al. used the Sobol global sensitivity analysis method to determine the critical geometric errors of machine tools [12]. Guo et al. used the first-order sensitivity analysis method to identify the sensitivity error of the five-axis machine tool and conducted a global sensitivity analysis of the machine tool based on an extend Fourier amplitude sensitivity test (EFAST) method [13]. Cheng et al. conducted a global sensitivity analysis of the error of a five-axis machine tool based on the EFAST method [14]. Cheng et al. conducted a global sensitivity analysis of machine tool errors based on the product of exponential (POE) screw theory and the Morris sensitivity analysis method [15]. Fu et al. established the error sensitivity matrix for each axis based on POE theory and evaluated the influence degree of each axis on the errors [16]. Xia et al. used the Morris global sensitivity analysis method to quantify the error sensitivity and identified the key errors and sensitive parts of the five-axis gear grinding machine [17]. Guo et al. modeled the error of a five-axis machine tool based on the MBS method and analyzed the sensitivity of the error model [18]. Li et al. considered the position and attitude errors of the tool and improved the traditional sensitivity analysis method by simplifying the output number of error sensitivity [19]. Zhang et al. used the multiplicative dimensional reduction method (M-DRM) to analyze the sensitivity of machine tools [20]. Zou et al. used the Sobol global sensitivity analysis method to analyze the error sensitivity of the three-axis diamond lathe [21]. Fan et al. developed a quantitative interval sensitivity analysis method (QISA) to calculate the sensitivity of geometric errors [22]. Li et al. defined a new machine tool error sensitivity analysis method by expressing the tool path and path error on a tool coordinate system [23]. Li et al. simplified the expression form of the machine tool error sensitivity index and proposed a new index to express error fluctuations [24].
At present, GSA is widely used in error sensitivity analysis. Among various methods, the Morris global sensitivity analysis is mainly used to analyze the impact of input parameters on output parameters when input parameters change in the global scope. The five-axis machine tool has a large number of PDGEs, and there is a coupling relationship between the error terms. The Morris global sensitivity analysis method, based on variable discrete and random sampling, can effectively analyze the sensitivity of the PDGEs of five-axis machine tools.
Although both LSA and GSA are effective methods for sensitive geometric error identification, they need error modeling and analysis process based on HTM, which has a low response speed when structure is reconfigured. Therefore, this paper tries to find a new way to calculate the sensitive geometric errors according to the machine tool structure without error modeling and analysis. The mapping expressions between the sensitive geometric errors and the machine tool structure are established and organized.
According to the different structural configurations of the five-axis machine tools, the sensitive PDGEs were analyzed by using the traditional error modeling method. Through the systematic analysis for different structures of RTTTR and TTTRR-type five-axis machine tools, the mapping expression between the sensitive PDGEs and the configuration of machine tools was established. Additionally, a high-efficient calculation method of sensitive PDGEs was proposed according to the mapping expressions without error modeling and analysis. The proposed method can directly identify sensitive PDGEs according to the configuration structure of the machine tool without requiring theoretical modeling and sensitivity analysis, which is suitable for quick error analysis of a machine tool.
The following sections are divided into four parts to display the contents of this article. Section 2 introduces the theoretical basis of the high-efficient calculation method. In Section 3, the high-efficient calculation method of sensitive PDGEs is introduced. In Section 4, the related simulation is described. Section 5 summarizes the results of this paper.

2. Theoretical Basis of the High-Efficient Calculation Method

This section introduces the theoretical basis of the high-efficient calculation method. Based on the traditional analysis method in this section, the RTTTR and TTTRR-type five-axis machine tools with different structures were analyzed. Through the systematic analysis for different structures of machine tools, eventually the mapping expressions between the sensitive PDGEs and the configuration of machine tools was established.

2.1. Morris Global Sensitivity Analysis Method

There have been a total of 36 PDGEs of the five-axis reconfigurable machine tool. The motion axes of a five-axis machine tool consist of X, Y, Z, A, B or X, Y, Z, A, C and so on. The names and geometric meanings of error are shown in Table 1.
This section takes the RTTTR-type five-axis machine tool as an example, as seen in Figure 1. Sub-coordinate system on each axis of the machine tool was established. Based on MBS theory [25,26] and HTM, the geometric error model of the five-axis machine tool was established. The Morris global sensitivity analysis method [27] was used to analyze sensitive PDGEs of the machine tool.
The cutting point coordinate in the tool coordinate system is set as [ 0 0 0 1 ] T . The machine tool is divided into two chains according to the topological relation. The workpiece chain is bed-X-C-workpiece, and the tool chain is bed-Y-Z-A-tool. Supposing that the ideal position transformation matrix between adjacent body i and j is T i j p , among them, i ( F , X , Y , Z , A , C ) and j ( X , Y , Z , A , C ) , the F, X, Y, Z, A and C represent bed, X axis, Y axis, Z axis, A axis, and C axis, respectively. The position error transformation matrix between adjacent body i and j is Δ T i j p . The ideal motion transformation matrix of adjacent body is T i j s . The motion error transformation matrix of adjacent body is Δ T i j s . Then, the ideal forming matrix of the machine tool is shown in Equation (1), and the actual forming matrix of the machine tool is shown in Equation (2). The comprehensive error matrix of the machine tool can be given as Equation (3).
P i d e a l = ( T F X p T F X s T X C p T X C s ) 1 T F Y p T F Y s T Y Z p T Y Z s T Z A p T Z A s [ 0 0 0 1 ] T
P a c t u a l = ( T F X p Δ T F X p T F X x s Δ T F X s T X C p Δ T X C p T X C s Δ T X C s ) 1 T F Y p Δ T F Y p T F Y s Δ T F Y s T Y Z p Δ T Y Z p T Y Z s Δ T Y Z s T Z A p Δ T Z A p T Z A s Δ T Z A s [ 0 0 0 1 ] T
E = P a c t u a l P i d e a l = [ E x E y E z 0 ] T
The Morris global sensitivity analysis method can effectively sort and evaluate the influence of model parameters; it has good applicability for models with many analysis parameters. Based on the experimental design of the one-time change method, only one model parameter is changed for each sampling, and the influence of each model parameter is calculated in turn so as to realize the evaluation of the sensitivity of the model parameters. The specific form is as follows.
E E i = y ( x 1 , , x i 1 , , x n ) y ( x 1 , , x i 1 + Δ , , x n ) Δ
In Equation (4), E E i is the influence effect of the ith input parameter, y is the output function, x is the input parameter, and Δ is the variation of the input parameter Δ = 1 / ( q 1 ) , q is the number of samples for each error term. The number of cyclic sampling is set to SN, and the influence effect of each input parameter is calculated in turn. In this way, the SN effects of each input parameter can be obtained.
μ i = 1 S N j = 1 S N E E i j σ i = 1 S N j = 1 S N ( E E i j μ i ) 2
As shown in Equation (5), the mean μ i and standard deviation σ i of each effect were calculated, respectively. Where μ i reflects the influence of each input parameter on the output function; the greater the value, the greater the influence. σ i reflects the coupling strength between each input parameter and other parameters; the larger the value, the higher the coupling strength.
The 30 PDGEs of the five-axis machine tool shown in Figure 1 were used as input parameters, and the comprehensive error model in Equation (3) was taken as output function by using the Morris method. The range of linear errors in PDGEs was set as [0, 15] μ m and for angular errors was set as [0, 0.015] deg. One-hundred and twenty (SN) cycles of sampling were conducted in the working area to obtain the mean value μ i and standard deviation σ i , so as to identify the sensitive PDGEs of the machine tool.
The sensitivity results are shown in Figure 2a–c. Serial numbers 1–6 represent the PDGEs of X axis ( Δ x ( X ) , Δ y ( X ) , Δ z ( X ) , Δ α ( X ) , Δ β ( X ) , and Δ γ ( X ) ), serial numbers 7–12 represent the PDGEs of Y axis( Δ x ( Y ) , Δ y ( Y ) , Δ z ( Y ) , Δ α ( Y ) , Δ β ( Y ) , and Δ γ ( Y ) ), serial numbers 13–18 represent the PDGEs of Z axis ( Δ x ( Z ) , Δ y ( Z ) , Δ z ( Z ) , Δ α ( Z ) , Δ β ( Z ) , and Δ γ ( Z ) ), serial numbers 19–24 represent the PDGEs of A axis ( Δ x ( A ) , Δ y ( A ) , Δ z ( A ) , Δ α ( A ) , Δ β ( A ) , and Δ γ ( A ) ) and serial numbers 25–30 represent the PDGEs of C axis ( Δ x ( C ) , Δ y ( C ) , Δ z ( C ) , Δ α ( C ) , Δ β ( C ) , and Δ γ ( C ) ), respectively. The red point represents the mean value μ i of the influence effect of each error term and reflects the sensitivity of each error term. The blue point σ i represents the coupling strength between each error term and other error terms. In this way, the sensitive PDGEs of the machine tool in three coordinate directions can be obtained.

2.2. Discovery of Mapping Relationships

Based on the error modeling and sensitive PDGEs analysis method described in Section 2.1, the sensitive PDGEs analyses of RTTTR and TTTRR-type five-axis reconfigurable machine tools with different structures were carried out, respectively. After comparing the analysis results of sensitive PDGEs with the structure of the machine tool, it is found that there is a mapping relationship between the sensitive PDGEs of the machine tool and the configuration of each axis of the machine tool. Some of these relationships are described below.
Case 1: The linear positioning errors of the translational axis corresponding to the identification direction are always sensitive PDGEs [14,24]. Examples are described as follows: in x direction, sensitive linear PDGEs of translational axes are Δ x ( X ) , Δ x ( Y ) , and Δ x ( Z ) . In y direction, sensitive linear PDGEs of translational axes are Δ y ( X ) , Δ y ( Y ) , and Δ y ( Z ) . In z direction, sensitive linear PDGEs of translational axes are Δ z ( X ) , Δ z ( Y ) , and Δ z ( Z ) . The definition of each error term is shown in Table 1.
Case 2: If the translational axis belongs to the workpiece chain, the axis has two sensitive angular PDGEs in each identification direction. The examples are as follows: assuming that Y axis belongs to workpiece chain, in x direction, the sensitive angular PDGEs of Y axis are Δ β ( Y ) and Δ γ ( Y ) ; in y direction, the sensitive angular PDGEs of Y axis are Δ α ( Y ) and Δ γ ( Y ) ; in z direction, the sensitive angular PDGEs of Y axis are Δ α ( Y ) and Δ β ( Y ) .
There are other relationships besides those mentioned above. In order to facilitate the identification, this paper presents these mapping relations in the form of mathematical expressions in Section 3. Combined with these mapping expressions, the sensitive PDGEs of the machine tool can be directly obtained through simple calculations.

3. Mapping Expressions

This section describes the mapping expression between the sensitive PDGEs of the five-axis machine tool and the machine tool structure. Based on this mapping expression, the sensitive PDGEs of the five-axis reconfigurable machine tool can be quickly identified, eliminating the process of error modeling and analysis.
The traditional process of sensitive geometric error analysis is shown in Figure 3a, and the proposed method is shown in Figure 3b. It can be found that the identification process of the sensitive geometric error of a single-structure five-axis machine tool is inherently complicated. Once the structure of the machine tool is reorganized, complex error modeling and analysis need to be performed again. In contrast with traditional methods, the proposed method can not only quickly identify the sensitive PDGEs of a specific five-axis machine tool but also quickly respond to the reorganized five-axis machine tool.
Before introducing the proposed method, the relevant symbols are defined and explained. This article uses lowercase italics when referring to directions, such as: x, y, and z and when referring to axes, uses uppercase italics, such as X, Y, Z, A, B, and C. Other related expressions are described in Section 3.1,Section 3.2 and Section 3.3, which describe the high-efficient calculation methods for sensitive PDGEs identification of RTTTR and TTTRR-type five-axis reconfigurable machine tools, respectively.

3.1. The Defined Symbols and Expressions

The new definition of mapping expressions of sensitive PDGEs identification is shown as Equation (6).
T y p ( i ) A x i s D i r = Δ C a l ( p a r m ) ( A x i s )
Where superscript Dir denotes the identification direction. Generally, the identification directions of spatial positional errors of a machine tool are x, y, and z directions along X-axis, Y-axis, and Z-axis, respectively. Axis denotes the motion axis of a machine tool, A x i s ( X , Y , Z , A , B , C ) . i = 1, 2, 3, or 4, which is used to distinguish different identification axes and different error types. Cal is the mapping algorithm; the type and direction of error can be obtained by Cal calculating parm. The detailed meanings of the defined symbols are described as follows:
(1)
T y p ( i ) A x i s D i r represents the sensitive PDGEs of the motion axis (Axis) in the identification direction (Dir). T y p ( i ) is the type of sensitive geometric error to be calculated and is divided into the following four categories: T y p ( 1 ) represents the sensitive linear PDGEs of the translational axis. T y p ( 2 ) represents the sensitive linear PDGEs of the rotational axis. T y p ( 3 ) represents the sensitive angular PDGEs of the translational axis. T y p ( 4 ) represents the sensitive angular PDGEs of rotational axis.
The sensitive PDGEs of each axis in each direction are identified in x, y, and z directions; D i r is the direction of error identification, D i r ( x , y , z ) .
A x i s indicates the motion axis that needs to be identified. Once T y p ( i ) is defined, the selection range of the corresponding axis is also determined. If i = 1 or 3, A x i s ( X , Y , Z ) . If i = 2 or 4, A x i s ( A , B , C ) .
For example, T y p ( 1 ) Y x represents the sensitive linear PDGEs of the Y axis in x direction. T y p ( 2 ) A z represents the sensitive linear PDGEs of the A axis in z direction. T y p ( 3 ) Y x represents the sensitive angular PDGEs of the Y axis in x direction. T y p ( 4 ) A z represents the sensitive angular PDGEs of the A axis in z direction.
(2)
The defined mapping expressions Δ C a l ( p a r m ) ( A x i s ) is used to calculate the error items in Table 1, taking Δ x ( X ) for an example, where C a l ( p a r m ) = x and A x i s = X.
In Δ C a l ( p a r m ) ( A x i s ) , Axis∈(X, Y, Z, A, B, C); parm∈(Dir, Axis, Toolvector), Toolvector is the tool axis vector (usually in z direction); when calculating sensitive linear PDGEs, C a l ( p a r m ) (x, y, z); when calculating the sensitive angular PDGEs, C a l ( p a r m ) ∈(α, β, γ). Taking X axis as an example, the calculation results of Δ C a l ( p a r m ) ( X ) are shown in Appendix A. The calculation results of other axes (Y, Z, A, B, or C axis) can also be obtained in this form.
In Appendix A, in order to facilitate understanding, the result of the C a l ( p a r m ) in C a l ( p a r m ) ( X ) is marked with two underscores. Among them, C a l ( p a r m ) represents the axial direction of parm; C a l ( p a r m ¯ ) represents the axial direction of the other two directions excepting the direction of parm; C a l ( p a r m 1 + p a r m 2 ) represents the axial direction of the parm1 and parm2; C a l ( p a r m 1 + p a r m 2 ¯ ) indicates the axial direction excepting parm1 and parm2.
(3)
Combined with the relevant definitions in (1) and (2), the relevant calculation forms of Equation (6) are illustrated with examples as follows:
  • If Equation (6) is T y p ( 1 ) X z = Δ C a l ( z ) ( X ) , then T y p ( 1 ) X z = Δ z ( X ) ;
  • If Equation (6) is T y p ( 3 ) X z = Δ C a l ( z ) ( X ) , then T y p ( 3 ) X z = Δ γ ( X ) ;
  • If Equation (6) is T y p ( 2 ) A y = Δ C a l ( y ¯ ) ( A ) , then T y p ( 2 ) A y = { Δ x ( A ) , Δ z ( A ) };
  • If Equation (6) is T y p ( 4 ) A y = Δ C a l ( y ¯ ) ( A ) , then T y p ( 4 ) A y = { Δ α ( A ) , Δ γ ( A ) }.

3.2. Calculation Method of RTTTR-Type Five-Axis Machine Tools

This section introduces the mapping expressions for sensitive PDGE identification of the RTTTR-type five-axis reconfigurable machine tool. According to the configuration axis information of the machine tool, combined with the basic calculation formula described in Section 3.1, the sensitive PDGEs can be quickly obtained by a simple calculation.
This section takes the machine tool shown in Figure 1 as an example to illustrate the proposed method. In Figure 1, the workpiece chain is: bed-X-C-workpiece, the tool chain is: bed-Y-Z-A-tool. The rotational axis in the workpiece chain is marked as R1, the rotational axis in the tool chain is marked as R2, the translational axis closest to the bed in the tool chain is marked as TF, and the axis vector of tool is z direction.

3.2.1. Sensitive Linear PDGEs for Translational Axes

If the sensitive linear PDGEs of translational axes are identified, i = 1, A x i s ( X , Y , Z ) can be obtained according to Section 3.1. Where D i r ( x , y , z ) , then the mapping expression is given as:
T y p ( 1 ) A x i s D i r = Δ C a l ( D i r ) ( A x i s )
According to Equation (7), take the X axis as an example:
If sensitivity PDGE in x direction is identified, then T y p ( 1 ) X x = Δ C a l ( x ) ( X ) = Δ x ( X ) can be obtained;
If sensitivity PDGE in y direction is identified, then T y p ( 1 ) X y = Δ C a l ( y ) ( X ) = Δ y ( X ) can be obtained;
If sensitivity PDGE in z direction is identified, then T y p ( 1 ) X z = Δ C a l ( z ) ( X ) = Δ z ( X ) can be obtained. In the same way, the sensitive linear PDGEs of the Y axis and Z axis can be calculated.

3.2.2. Sensitive Linear PDGEs for Rotational Axes

If the sensitive linear PDGEs of rotational axes are identified, i = 2, A x i s ( A , B , C ) can be obtained according to Section 3.1. Where D i r ( x , y , z ) , then the mapping expression is given as:
{ if : C a l ( D i r ) = C a l ( A x i s ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( A x i s ) ( A x i s ) if : C a l ( D i r ) C a l ( A x i s ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( A x i s ¯ ) ( A x i s )
According to Equation (8), take the A axis as an example:
If sensitivity PDGE in x direction is identified, C a l ( x ) = C a l ( A ) = ( C a l ( D i r ) = C a l ( A x i s ) ) , then T y p ( 2 ) A x = Δ C a l ( A ) ( A ) = Δ x ( A ) can be obtained;
If sensitivity PDGE in y direction is identified, C a l ( y ) C a l ( A ) ( C a l ( D i r ) C a l ( A x i s ) ), then T y p ( 2 ) A y = Δ C a l ( A ¯ ) ( A ) = { Δ y ( A ) , Δ z ( A ) } can be obtained;
If sensitivity PDGE in z direction is identified, C a l ( z ) C a l ( A ) ( C a l ( D i r ) C a l ( A x i s ) ), then T y p ( 2 ) A z = Δ C a l ( A ¯ ) ( A ) = { Δ y ( A ) , Δ z ( A ) } , can be obtained.
Therefore, the sensitive linear PDGEs of A axis in the x direction is Δ x ( A ) , in the y direction are Δ y ( A ) and Δ z ( A ) , and in the z direction are Δ y ( A ) and Δ z ( A ) . In the same way, the sensitive linear PDGEs of the other rotational axis (C axis) in the x, y, and z directions can be calculated.

3.2.3. Sensitive Angular PDGEs of Translational Axis

If the sensitive linear PDGEs of translational axes are identified, i = 3, A x i s ( X , Y , Z ) can be obtained according to Section 3.1. Where D i r ( x , y , z ) , then the mapping expression is given as:
(1)
When A x i s belongs to the workpiece chain
T y p ( 3 ) A x i s D i r = Δ C a l ( D i r ¯ ) ( A x i s )
(2)
When A x i s belongs to the tool chain
{ if : C a l ( D i r ) = C a l ( R 2 ) , then : T y p ( 3 ) A x i s D i r = Δ C a l ( R 2 ¯ ) ( A x i s ) if : C a l ( D i r ) C a l ( R 2 ) , then : T y p ( 3 ) A x i s D i r = Δ C a l ( R 2 ) ( A x i s ) if :   tool   chain   has   3   translational   axes ,   and C a l ( T F ) C a l ( R 1 + R 2 ) , then : T y p ( 3 ) A x i s D i r = Δ C a l ( D i r ¯ ) ( T F )
According to Equation (9), take the X axis as an example (X axis belongs to the workpiece chain):
If sensitivity PDGE in x direction is identified, then T y p ( 3 ) X x = Δ C a l ( x ¯ ) ( X ) = { Δ β ( X ) , Δ γ ( X ) } can be obtained;
If sensitivity PDGE in y direction is identified, then T y p ( 3 ) X y = Δ C a l ( y ¯ ) ( X ) = { Δ α ( X ) , Δ γ ( X ) } can be obtained;
If sensitivity PDGE in z direction is identified, then T y p ( 3 ) X z = Δ C a l ( z ¯ ) ( X ) = { Δ α ( X ) , Δ β ( X ) } can be obtained.
Therefore, the sensitive angular PDGEs of X axis in the x direction are Δ β ( X ) and Δ γ ( X ) ; in the y direction are Δ α ( X ) and Δ γ ( X ) ; and in the z direction are Δ α ( X ) and Δ β ( X ) .
According to Equation (10), take the Y axis as an example (Y axis belongs to the tool chain):
If sensitivity PDGE in x direction is identified, C a l ( x ) = C a l ( A ) ( C a l ( D i r ) = C a l ( R 2 ) ), then T y p ( 3 ) Y x = Δ C a l ( A ¯ ) ( Y ) = { Δ β ( Y ) , Δ γ ( Y ) } can be obtained;
If sensitivity PDGE in y direction is identified, C a l ( y ) C a l ( A ) ( C a l ( D i r ) C a l ( R 2 ) ), then T y p ( 3 ) Y y = Δ C a l ( A ) ( Y ) = Δ α ( Y ) can be obtained;
If sensitivity PDGE in z direction is identified, C a l ( z ) C a l ( A ) ( C a l ( D i r ) C a l ( R 2 ) ), then T y p ( 3 ) Y z = Δ C a l ( A ) ( Y ) = Δ α ( Y ) can be obtained.
Therefore, the sensitive angular PDGEs of Y axis in the x direction are Δ β ( Y ) and Δ γ ( Y ) ; in the y direction is Δ α ( Y ) ; and in the z direction is Δ α ( Y ) . In the same way, the sensitive angular PDGEs for the Z axis of the tool chain in the x, y, and z directions can be calculated.

3.2.4. Sensitive Angular PDGEs of Rotational Axis

If the sensitive linear PDGEs of rotational axes are identified, i = 4, A x i s ( A , B , C ) can be obtained according to Section 3.1. Where D i r ( x , y , z ) , then the mapping expression is given as:
(1)
When A x i s = R1
{ if : C a l ( D i r ) = C a l ( R 1 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( R 1 ¯ ) ( R 1 ) if : C a l ( D i r ) C a l ( R 1 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( X + Y + Z ) ( R 1 )
(2)
When A x i s = R2
{ if : C a l ( D i r ) = C a l ( T o o l v e c t o r ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( R 2 ) ( R 2 ) if : C a l ( D i r ) C a l ( T o o l v e c t o r ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( D i r + T o o l v e c t o r ¯ ) ( R 2 )
According to Equation (11), take the C axis as an example (C axis is marked as R1 axis):
  • If sensitivity PDGE in x direction is identified, C a l ( x ) C a l ( C ) ( C a l ( D i r ) C a l ( R 1 ) ), then T y p ( 4 ) C x = Δ C a l ( X + Y + Z ) ( C ) = { Δ α ( C ) , Δ β ( C ) , Δ γ ( C ) } can be obtained;
  • If sensitivity PDGE in y direction is identified, C a l ( y ) C a l ( C ) ( C a l ( D i r ) C a l ( R 1 ) ), then T y p ( 4 ) C y = Δ C a l ( X + Y + Z ) ( C ) = { Δ α ( C ) , Δ β ( C ) , Δ γ ( C ) } can be obtained;
  • If sensitivity PDGE in z direction is identified, C a l ( z ) = C a l ( C ) ( C a l ( D i r ) = C a l ( R 1 ) ), then T y p ( 4 ) C z = Δ C a l ( C ¯ ) ( C ) = { Δ α ( C ) , Δ β ( C ) } can be obtained.
Therefore, the sensitive angular PDGEs of C axis in the x direction are Δ α ( C ) , Δ β ( C ) , and Δ γ ( C ) ; in the y direction are Δ α ( C ) , Δ β ( C ) , and Δ γ ( C ) ; and in the z direction are Δ α ( C ) and Δ β ( C ) .
According to Equation (12), take the A axis as an example (A axis is marked as R2 axis):
  • If sensitivity PDGE in x direction is identified, C a l ( x ) C a l ( T o o l v e c t o r ) ( C a l ( D i r ) C a l ( T o o l v e c t o r ) , then T y p ( 4 ) A x = Δ C a l ( x + z ¯ ) ( A ) = Δ β ( A ) can be obtained;
  • If sensitivity PDGE in y direction is identified, C a l ( y ) C a l ( T o o l v e c t o r ) ( C a l ( D i r ) C a l ( T o o l v e c t o r ) ), then T y p ( 4 ) A y = Δ C a l ( y + z ¯ ) ( A ) = Δ α ( A ) can be obtained;
  • If sensitivity PDGE in z direction is identified, C a l ( z ) = C a l ( T o o l v e c t o r ) ( C a l ( D i r ) = C a l ( T o o l v e c t o r ) ), then T y p ( 4 ) A z = Δ C a l ( A ) ( A ) = Δ α ( A ) can be obtained.
Therefore, the sensitive angular PDGEs of A axis in the x direction is Δ β ( A ) ; in the y direction is Δ α ( A ) ; and in the z direction is Δ α ( A ) .
In summary, all the sensitive geometric errors in the x, y, and z directions are shown in Table 2.

3.3. Calculation Method of TTTRR-Type five-Axis Machine Tools

This section introduces the mapping expression for sensitive PDGEs of the TTTRR-type five-axis reconfigurable machine tool. According to this mapping expression, sensitive PDGEs can be quickly calculated without error modeling and analysis.
The rotational axis closest to the tool is marked as R1, the rotational axis farther from the tool is marked as R2, and the axis vector of tool is z direction. The meaning of the symbols in the mapping expression is the same as in Section 3.2; the specific mapping expression is as follows.

3.3.1. Sensitive Linear PDGEs of Translational Axis

If the sensitive linear PDGEs of translational axes are identified, i = 1, A x i s ( X , Y , Z ) can be obtained according to Section 3.1. Where D i r ( x , y , z ) , then the mapping expression is shown as Equation (7).

3.3.2. Sensitive Linear PDGEs of Rotational Axes

If the sensitive linear PDGEs of rotational axes are identified, i = 2, A x i s ( A , B , C ) can be obtained according to Section 3.1. Where D i r ( x , y , z ) , then the mapping expression is given as:
(1)
When A x i s = R 1
{ if : C a l ( D i r ) = C a l ( R 2 ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( R 1 ¯ ) ( R 1 ) if : C a l ( D i r ) C a l ( R 2 ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( X + Y + Z ) ( R 1 )
(2)
When A x i s = R 2
{ if : C a l ( D i r ) = C a l ( R 2 ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( R 2 ) ( R 2 ) if : C a l ( D i r ) C a l ( R 2 ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( R 2 ¯ ) ( R 2 )

3.3.3. Sensitive Angular PDGEs of Translational Axis

If the sensitive linear PDGEs of rotational axes are identified, i = 3, A x i s ( X , Y , Z ) can be obtained according to Section 3.1. Where D i r ( x , y , z ) , then the mapping expression is given as:
T y p ( 3 ) A x i s D i r = Δ C a l ( D i r ¯ ) ( A x i s )

3.3.4. Sensitive Angular PDGEs of Rotational Axis

If the sensitive linear PDGEs of rotational axes are identified, i = 4, A x i s ( A , B , C ) can be obtained according to Section 3.1. Where D i r ( x , y , z ) , then the mapping expression is given as:
(1)
When A x i s = R 1
{ if : C a l ( D i r ) = C a l ( R 2 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( R 1 ) ( R 1 ) if : C a l ( D i r ) C a l ( R 2 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( T o o l v e c t o r ¯ ) ( R 1 )
(2)
When A x i s = R 2
{ if : C a l ( D i r ) = C a l ( R 2 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( R 1 ) ( R 2 ) if : C a l ( D i r ) C a l ( R 2 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( X + Y + Z ) ( R 2 )

3.4. Summary of Mapping Expressions

Using the proposed method to compare with the existing literature [19,24], the results of the literature analysis are consistent with the sensitive PDGEs identified by this method. Table 3 and Table 4 present a summary of the calculation methods for RTTTR and TTTRR-type five-axis machine tool sensitive PDGEs. Based on these method, sensitive PDGES can be quickly calculated, and it provides guidance for accurate design and error traceability of machine tools.

4. Simulation Analysis

The impeller is a typical part in five-axis machining. In order to verify the accuracy and effectiveness of this calculation method, an impeller was simulated and analyzed. The value of PDGEs were preset, the linear PDGEs were set as 0.01 mm × sin(0.5 × position), and the angular PDGEs were set as 0.01 deg × sin(0.5 × angle). The value range of PDGEs was borrowed from literature [12,14,15,28].
VERICUT software was used to simulate the machining of the impeller, and the machining errors before and after compensation were compared. Firstly, 100 analysis points were selected on the surface of one side of the impeller model, as shown in Figure 4a. Then, based on the sensitive geometric error analysis results in Table 2, the actual inverse kinematics [28] were used to calculate the NC codes after compensation in each direction.
Taking the analysis results in x direction as an example, Figure 4b shows the impeller processing diagram after compensating for the sensitive PDGEs. Figure 5 shows some NC codes before and after compensating for the sensitive PDGEs in x direction. Figure 6a shows the theoretical error and the machining error before and after compensating for the sensitive PDGEs in x direction; the values in Figure 6 are the maximum machining error before and after compensation. The machining error before and after compensation in the y and z direction sensitive PDGEs are shown in Figure 6b,c, respectively.
It can be found that the machining error after compensating for these sensitive PDGEs is almost consistent with the theoretical error, which indicates that the machining accuracy is mainly affected by these sensitive PDGEs. At the same time, the machining error after compensation is significantly lower than that before processing, which indicates the accuracy of the proposed method.
According to the proposed method, Matlab software was used to make related software. As shown in Figure 7, the sensitivity index of each PDGE was calculated through the function embedded in the software by inputting the topology of the machine tool, and the sensitive PDGEs were displayed.

5. Conclusions

A new high-efficient calculation method for sensitive PDGE identification in the five-axis reconfigurable machine tool was proposed and verified in this paper. According to machine tool configuration axis information, sensitive PDGEs of RTTTR and TTTRR-type five-axis reconfigurable machine tools can be quickly identified according to the proposed method without traditional HTM modeling and sensitivity analysis.
The simulation results show that the machining accuracy of the impeller is significantly improved after compensating for sensitive PDGEs. According to the simulation results, the proposed method is accurate. It can provide a reference for precision design and error compensation of five-axis reconfigurable machine tools.
Based on the proposed calculation method, the barrier-free output of sensitive PDGEs can be realized through the programming of related software. It can quickly respond to the sensitive PDGEs identification in the reconstructed five-axis machine tool. Rapid identification of sensitive PDGEs for five-axis reconfigurable machine tools is studied in this paper. In future work, the calculation methods of other error sources of five-axis reconfigurable machine tools are also worth further study.

Author Contributions

Z.S. contributed to the conception of the study. S.D. helped perform the analysis with constructive discussions. Z.C., and Z.L. performed the simulation. Z.S., and Z.W. contributed significantly to manuscript preparation. All authors have read and approved the manuscript being submitted and agree to its submission to this journal.

Funding

National Natural Science Foundation of China (No. 51905470), China Postdoctoral Science Foundation (No. 2020M671617), Natural Science Foundation of the Jiangsu Higher Education Institutions of China (No. 19KJB460029), and Jiangsu Postgraduate Research Innovation Project of China (No. SJCX20_1373).

Data Availability Statement

The datasets used or analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

Thanks for the support of Mechanical Engineering College of Yangzhou University.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Linear PDGEsAngular PDGEs
Δ C a l ( p a r m ) ( X ) Δ C a l ( p a r m ¯ ) ( X ) Δ C a l ( p a r m ) ( X ) Δ C a l ( p a r m ¯ ) ( X )
Δ C a l ( x ) _ _ ( X ) = Δ x _ _ ( X ) Δ C a l ( x ¯ ) _ _ ( X ) = { Δ y _ _ ( X ) , Δ z _ _ ( X ) } Δ C a l ( x ) _ _ ( X ) = Δ α _ _ ( X ) Δ C a l ( x ¯ ) _ _ ( X ) = { Δ β _ _ ( X ) , Δ γ _ _ ( X ) }
Δ C a l ( y ) _ _ ( X ) = Δ y _ _ ( X ) Δ C a l ( y ¯ ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ z _ _ ( X ) } Δ C a l ( y ) _ _ ( X ) = Δ β _ _ ( X ) Δ C a l ( y ¯ ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ γ _ _ ( X ) }
Δ C a l ( z ) _ _ ( X ) = Δ z _ _ ( X ) Δ C a l ( z ¯ ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ y _ _ ( X ) } Δ C a l ( z ) _ _ ( X ) = Δ γ _ _ ( X ) Δ C a l ( z ¯ ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ β _ _ ( X ) }
Δ C a l ( X ) _ _ ( X ) = Δ x _ _ ( X ) Δ C a l ( X ¯ ) _ _ ( X ) = { Δ y _ _ ( X ) , Δ z _ _ ( X ) } Δ C a l ( X ) _ _ ( X ) = Δ α _ _ ( X ) Δ C a l ( X ¯ ) _ _ ( X ) = { Δ β _ _ ( X ) , Δ γ _ _ ( X ) }
Δ C a l ( Y ) _ _ ( X ) = Δ y _ _ ( X ) Δ C a l ( Y ¯ ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ z _ _ ( X ) } Δ C a l ( Y ) _ _ ( X ) = Δ β _ _ ( X ) Δ C a l ( Y ¯ ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ γ _ _ ( X ) }
Δ C a l ( Z ) _ _ ( X ) = Δ z _ _ ( X ) Δ C a l ( Z ¯ ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ y _ _ ( X ) } Δ C a l ( Z ) _ _ ( X ) = Δ γ _ _ ( X ) Δ C a l ( Z ¯ ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ β _ _ ( X ) }
Δ C a l ( x + y ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ y _ _ ( X ) } Δ C a l ( x + y ¯ ) _ _ ( X ) = Δ z _ _ ( X ) Δ C a l ( x + y ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ β _ _ ( X ) } Δ C a l ( x + y ¯ ) _ _ ( X ) = Δ γ _ _ ( X )
Δ C a l ( x + z ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ z _ _ ( X ) } Δ C a l ( x + z ¯ ) _ _ ( X ) = Δ y _ _ ( X ) Δ C a l ( x + z ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ γ _ _ ( X ) } Δ C a l ( x + z ¯ ) _ _ ( X ) = Δ β _ _ ( X )
Δ C a l ( y + z ) _ _ ( X ) = { Δ y _ _ ( X ) , Δ z _ _ ( X ) } Δ C a l ( y + z ¯ ) _ _ ( X ) = Δ x _ _ ( X ) Δ C a l ( y + z ) _ _ ( X ) = { Δ β _ _ ( X ) , Δ γ _ _ ( X ) } Δ C a l ( y + z ¯ ) _ _ ( X ) = Δ α _ _ ( X )
Δ C a l ( x + y + z ) _ _ ( X )   =
{ Δ x _ _ ( X ) , Δ y _ _ ( X ) , Δ z _ _ ( X ) }
Δ C a l ( X + Y ¯ ) _ _ ( X ) = Δ z _ _ ( X ) Δ C a l ( x + y + z ) _ _ ( X )   =
{ Δ α _ _ ( X ) , Δ β _ _ ( X ) , Δ γ _ _ ( X ) }
Δ C a l ( X + Y ¯ ) _ _ ( X ) = Δ γ _ _ ( X )
Δ C a l ( X + Y ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ y _ _ ( X ) } Δ C a l ( X + Z ¯ ) _ _ ( X ) = Δ y _ _ ( X ) Δ C a l ( X + Y ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ β _ _ ( X ) } Δ C a l ( X + Z ¯ ) _ _ ( X ) = Δ β _ _ ( X )
Δ C a l ( X + Z ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ z _ _ ( X ) } Δ C a l ( Y + Z ¯ ) _ _ ( X ) = Δ x _ _ ( X ) Δ C a l ( X + Z ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ γ _ _ ( X ) } Δ C a l ( Y + Z ¯ ) _ _ ( X ) = Δ α _ _ ( X )
Δ C a l ( Y + Z ) _ _ ( X ) = { Δ y _ _ ( X ) , Δ z _ _ ( X ) } Δ C a l ( A ¯ ) _ _ ( X ) = { Δ y _ _ ( X ) Δ z _ _ ( X ) } Δ C a l ( Y + Z ) _ _ ( X ) = { Δ β _ _ ( X ) , Δ γ _ _ ( X ) } Δ C a l ( A ¯ ) _ _ ( X ) = { Δ β _ _ ( X ) , Δ γ _ _ ( X ) }
Δ C a l ( X + Y + Z ) _ _ ( X )   =
Δ x _ _ ( X ) , Δ y _ _ ( X ) , Δ z _ _ ( X ) }
Δ C a l ( B ¯ ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ z _ _ ( X ) } Δ C a l ( X + Y + Z ) _ _ ( X )   =
{ Δ α _ _ ( X ) , Δ β _ _ ( X ) , Δ γ _ _ ( X ) }
Δ C a l ( B ¯ ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ γ _ _ ( X ) }
Δ C a l ( A ) _ _ ( X ) = Δ x _ _ ( X ) Δ C a l ( C ¯ ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ y _ _ ( X ) } Δ C a l ( A ) _ _ ( X ) = Δ α _ _ ( X ) Δ C a l ( C ¯ ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ β _ _ ( X ) }
Δ C a l ( B ) _ _ ( X ) = Δ y _ _ ( X ) Δ C a l ( B ) _ _ ( X ) = Δ β _ _ ( X )
Δ C a l ( C ) _ _ ( X ) = Δ z _ _ ( X ) Δ C a l ( C ) _ _ ( X ) = Δ γ _ _ ( X )
Δ C a l ( A + B ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ y _ _ ( X ) } Δ C a l ( A + B ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ β _ _ ( X ) }
Δ C a l ( A + C ) _ _ ( X ) = { Δ x _ _ ( X ) , Δ z _ _ ( X ) } Δ C a l ( A + C ) _ _ ( X ) = { Δ α _ _ ( X ) , Δ γ _ _ ( X ) }
Δ C a l ( B + C ) _ _ ( X ) = { Δ y _ _ ( X ) , Δ z _ _ ( X ) } Δ C a l ( B + C ) _ _ ( X ) = { Δ β _ _ ( X ) , Δ γ _ _ ( X ) }

References

  1. Ni, J. CNC Machine Accuracy Enhancement Through Real-Time Error Compensation. J. Manuf. Sci. Eng. 1997, 119, 717–725. [Google Scholar] [CrossRef]
  2. Koren, Y.; Heisel, U.; Jovane, F.; Moriwaki, T.; Pritschow, G.; Ulsoy, G.; Van Brussel, H. Reconfigurable Manufacturing System. Ann. CIRP 1999, 48, 527–540. [Google Scholar] [CrossRef]
  3. Yao, Y.; Liu, M.; Du, J.; Zhou, L. Design of a machine tool control system for function reconfiguration and reuse in network environment. Robot. Comput. Manuf. 2019, 56, 117–126. [Google Scholar] [CrossRef]
  4. Gadalla, M.; Xue, D.Y. Recent advances in research on reconfigurable machine tools: A literature review. Int. J. Prod. Res. 2017, 55, 1440–1454. [Google Scholar] [CrossRef]
  5. Álvarez, Á.; Calleja, A.; Ortega, N.; López de Lacalle, L.N. Five-axis milling of large spiral bevel gears: Tool path definition, finishing, and shape errors. Metals 2018, 8, 353. [Google Scholar] [CrossRef] [Green Version]
  6. McCulloch, A. Sensitivity analysis in practice: A guide to assessing scientific models. J. R. Stat. Soc. Ser. Stat. Soc. 2005, 168, 466. [Google Scholar] [CrossRef]
  7. Fan, K.-C.; Wang, H.; Zhao, J.-W.; Chang, T.-H. Sensitivity analysis of the 3-PRS parallel kinematic spindle platform of a serial-parallel machine tool. Int. J. Mach. Tools Manuf. 2003, 43, 1561–1569. [Google Scholar] [CrossRef]
  8. Li, D.; Feng, P.; Zhang, J.; Yu, D.; Wu, Z. An identification method for key geometric errors of machine tool based on matrix differential and experimental test. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2014, 228, 3141–3155. [Google Scholar] [CrossRef]
  9. Chen, G.; Liang, Y.; Sun, Y.; Chen, W.; Wang, B. Volumetric error modeling and sensitivity analysis for designing a five-axis ultra-precision machine tool. Int. J. Adv. Manuf. Technol. 2013, 68, 2525–2534. [Google Scholar] [CrossRef]
  10. Cheng, Q.; Liu, G.; Liu, Z.; Xuan, D.; Chang, W. An identification approach for key geometric error sources of machine tool based on sensitivity analysis. J. Mech. Eng. 2012, 48, 171–179. (In Chinese) [Google Scholar] [CrossRef]
  11. Yao, H.; Li, Z.; Zhao, X.; Sun, T.; Dobrovolskyi, G.; Li, G. Modeling of kinematics errors and alignment method of a swing arm ultra-precision diamond turning machine. Int. J. Adv. Manuf. Technol. 2016, 87, 165–176. [Google Scholar] [CrossRef]
  12. Cheng, Q.; Zhao, H.; Zhang, G.; Gu, P.; Cai, L. An analytical approach for crucial geometric errors identification of multi-axis machine tool based on global sensitivity analysis. Int. J. Adv. Manuf. Technol. 2014, 75, 107–121. [Google Scholar] [CrossRef]
  13. Guo, S.; Zhang, D.; Xi, Y. Global Quantitative Sensitivity Analysis and Compensation of Geometric Errors of CNC Machine Tool. Math. Probl. Eng. 2016, 2016, 2834718. [Google Scholar] [CrossRef] [Green Version]
  14. Cheng, Q.; Sun, B.; Liu, Z.; Li, J.; Dong, X.; Gu, P. Key geometric error extraction of machine tool based on extended Fourier amplitude sensitivity test method. Int. J. Adv. Manuf. Technol. 2017, 90, 3369–3385. [Google Scholar] [CrossRef]
  15. Cheng, Q.; Feng, Q.; Liu, Z.; Gu, P.; Zhang, G. Sensitivity analysis of machining accuracy of multi-axis machine tool based on POE screw theory and Morris method. Int. J. Adv. Manuf. Technol. 2016, 84, 2301–2318. [Google Scholar] [CrossRef]
  16. Fu, G.; Gong, H.; Fu, J.; Gao, H.; Deng, X. Geometric error contribution modeling and sensitivity evaluating for each axis of five-axis machine tools based on POE theory and transforming differential changes between coordinate frames. Int. J. Mach. Tools Manuf. 2019, 147, 103455. [Google Scholar] [CrossRef]
  17. Changjiu, X.; Shilong, W.; Shouli, S.; Xiaochuan, L.; Xiaodiao, H. Geometric error to tooth surface error model and identification of crucial errors in five-axis CNC gear profile grinding machines. Comput. Integr. Manuf. Syst. 2020, 26, 1191–1201. (In Chinese) [Google Scholar]
  18. Guo, S.; Jiang, G.; Mei, X. Investigation of sensitivity analysis and compensation parameter optimization of geometric error for five-axis machine tool. Int. J. Adv. Manuf. Technol. 2017, 93, 3229–3243. [Google Scholar] [CrossRef]
  19. Li, Q.; Wang, W.; Jiang, Y.; Li, H.; Zhang, J.; Jiang, Z. A sensitivity method to analyze the volumetric error of five-axis machine tool. Int. J. Adv. Manuf. Technol. 2018, 98, 1791–1805. [Google Scholar] [CrossRef]
  20. Zhang, X.; Zhang, Y.; Pandey, M.D. Global sensitivity analysis of a CNC machine tool: Application of MDRM. Int. J. Adv. Manuf. Technol. 2015, 81, 159–169. [Google Scholar] [CrossRef]
  21. Zou, X.; Zhao, X.; Li, G.; Li, Z.; Sun, T. Sensitivity analysis using a variance-based method for a three-axis diamond turning machine. Int. J. Adv. Manuf. Technol. 2017, 92, 4429–4443. [Google Scholar] [CrossRef]
  22. Fan, J.; Tao, H.; Pan, R.; Chen, D. An approach for accuracy enhancement of fifive-axis machine tools based on quantitative interval sensitivity analysis. Mech. Mach. Theory 2020, 148, 103806. [Google Scholar] [CrossRef]
  23. Li, Z.; Sato, R.; Shirase, K.; Ihara, Y.; Milutinovic, D.S. Sensitivity analysis of relationship between error motions and machined shape errors in five-axis machining center-Peripheral milling using squareend mill as test case. Precis. Eng. 2019, 60, 28–41. [Google Scholar] [CrossRef]
  24. Li, J.; Xie, F.; Liu, X.-J. Geometric error modeling and sensitivity analysis of a five-axis machine tool. Int. J. Adv. Manuf. Technol. 2016, 82, 2037–2051. [Google Scholar] [CrossRef]
  25. Okafor, A.C.; Ertekin, Y.M. Derivation of machine tool error models and error compensation procedure for three axes vertical machining center using rigid body kinematics. Int. J. Mach. Tools Manuf. 2000, 40, 1199–1213. [Google Scholar] [CrossRef]
  26. Fan, J.W.; Guan, J.L.; Wang, W.C.; Luo, Q.; Zhang, X.; Wang, L. A universal modeling method for enhancement the volumetric accuracy of CNC machine tools. J. Mater. Process. Technol. 2002, 129, 624–628. [Google Scholar] [CrossRef]
  27. Crosetto, M.; Tarantola, S. Uncertainty and sensitivity analysis: Tools for GIS-based model implementation. Int. J. Geogr. Inf. Sci. 2001, 15, 415–437. [Google Scholar] [CrossRef]
  28. Ding, S.; Huang, X.; Yu, C.; Wang, W. Actual inverse kinematics for position-independent and position-dependent geometric error compensation of five-axis machine tools. Int. J. Mach. Tools Manuf. 2016, 111, 55–62. [Google Scholar] [CrossRef]
Figure 1. RTTTR-type five-axis machine tool structure diagram.
Figure 1. RTTTR-type five-axis machine tool structure diagram.
Machines 09 00084 g001
Figure 2. The sensitivity index of the RTTTR-type machine tool.
Figure 2. The sensitivity index of the RTTTR-type machine tool.
Machines 09 00084 g002aMachines 09 00084 g002b
Figure 3. The process of identifying the sensitive geometric errors of the five-axis reconfigurable machine tool.
Figure 3. The process of identifying the sensitive geometric errors of the five-axis reconfigurable machine tool.
Machines 09 00084 g003
Figure 4. Impeller model diagram.
Figure 4. Impeller model diagram.
Machines 09 00084 g004
Figure 5. Part of the NC codes.
Figure 5. Part of the NC codes.
Machines 09 00084 g005
Figure 6. Comparison of results before and after compensation in each direction.
Figure 6. Comparison of results before and after compensation in each direction.
Machines 09 00084 g006
Figure 7. Sensitive PDGEs calculation software.
Figure 7. Sensitive PDGEs calculation software.
Machines 09 00084 g007
Table 1. Definition of PDGEs.
Table 1. Definition of PDGEs.
The Name of the Linear Error TermThe Geometric Meaning of the Error TermThe Name of the Angular Error TermThe Geometric Meaning of the Error Term
Δ x ( X ) Positioning error of the X axis in the x direction Δ α ( X ) Rolling error of X axis in x direction
Δ y ( X ) Straightness error of X axis in y direction Δ β ( X ) X axis angular error in y direction
Δ z ( X ) Straightness error of X axis in z direction Δ γ ( X ) X axis angular error in z direction
Δ x ( Y ) Straightness error of Y axis in x direction Δ α ( Y ) Y axis angular error in x direction
Δ y ( Y ) Positioning error of Y axis in y direction Δ β ( Y ) Rolling error of Y axis in y direction
Δ z ( Y ) Straightness error of Y axis in z direction Δ γ ( Y ) Y axis angular error in z direction
Δ x ( Z ) Straightness error of Z axis in x direction Δ α ( Z ) Z axis angular error in x direction
Δ y ( Z ) Straightness error of Z axis in y direction Δ β ( Z ) Z axis angular error in y direction
Δ z ( Z ) Positioning error of Z axis in z direction Δ γ ( Z ) Rolling error of Z axis in z direction
Δ x ( A ) Positioning error of the A axis in the x direction Δ α ( A ) Rolling error of A axis in x direction
Δ y ( A ) Straightness error of A axis in y direction Δ β ( A ) A axis angular error in y direction
Δ z ( A ) Straightness error of A axis in z direction Δ γ ( A ) A axis angular error in z direction
Δ x ( B ) Straightness error of B axis in x direction Δ α ( B ) B axis angular error in x direction
Δ y ( B ) Positioning error of B axis in y direction Δ β ( B ) Rolling error of B axis in y direction
Δ z ( B ) Straightness error of B axis in z direction Δ γ ( B ) B axis angular error in z direction
Δ x ( C ) Straightness error of C axis in x direction Δ α ( C ) C axis angular error in x direction
Δ y ( C ) Straightness error of C axis in y direction Δ β ( C ) C axis angular error in y direction
Δ z ( C ) Positioning error of C axis in z direction Δ γ ( C ) Rolling error of C axis in z direction
Table 2. Sensitive PDGEs in all directions.
Table 2. Sensitive PDGEs in all directions.
Error DirectionSensitive PDGEs
x Δ x ( X ) , Δ x ( Y ) , Δ x ( Z ) , Δ x ( A ) , Δ x ( C ) , Δ y ( C ) , Δ β ( X ) , Δ γ ( X ) , Δ β ( Y ) , Δ γ ( Y ) , Δ β ( Z ) , Δ γ ( Z ) , Δ α ( C ) , Δ β ( C ) , Δ γ ( C ) , Δ β ( A )
y Δ y ( X ) ,   Δ y ( Y ) ,   Δ y ( Z ) ,   Δ y ( A ) ,   Δ z ( A ) , Δ x ( C ) , Δ y ( C ) , Δ α ( X ) , Δ γ ( X ) , Δ α ( Y ) , Δ α ( Z ) ,
Δ α ( C ) , Δ β ( C ) , Δ γ ( C ) , Δ α ( A )
z Δ z ( X ) , Δ z ( Y ) , Δ z ( Z ) , Δ y ( A ) , Δ z ( A ) , Δ z ( C ) , Δ α ( X ) , Δ β ( X ) , Δ α ( Y ) , Δ α ( Z ) , Δ α ( C ) , Δ β ( C ) , Δ α ( A )
Table 3. Mapping expressions of RTTTR-type five-axis machine tool.
Table 3. Mapping expressions of RTTTR-type five-axis machine tool.
Sensitive Linear PDGEsSensitive Angular PDGEs
Translational axis in the workpiece chain T y p ( 1 ) A x i s D i r = Δ C a l ( D i r ) ( A x i s ) T y p ( 3 ) A x i s D i r = Δ C a l ( D i r ¯ ) ( A x i s )
Translational axis in the tool chain T y p ( 1 ) A x i s D i r = Δ C a l ( D i r ) ( A x i s ) { if : C a l ( D i r ) = C a l ( R 2 ) , then : T y p ( 3 ) A x i s D i r = Δ C a l ( R 2 ¯ ) ( A x i s ) if : C a l ( D i r ) C a l ( R 2 ) , then : T y p ( 3 ) A x i s D i r = Δ C a l ( R 2 ) ( A x i s ) if :   tool   chain   has   3   translational   axes , and C a l ( T F ) C a l ( R 1 + R 2 ) , then : T y p ( 3 ) A x i s D i r = Δ C a l ( D i r ¯ ) ( T F )
Rotational axis in the workpiece chain { if : C a l ( D i r ) = C a l ( A x i s ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( A x i s ) ( A x i s ) if : C a l ( D i r ) C a l ( A x i s ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( A x i s ¯ ) ( A x i s ) { if : C a l ( D i r ) = C a l ( R 1 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( R 1 ¯ ) ( R 1 ) if : C a l ( D i r ) C a l ( R 1 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( X + Y + Z ) ( R 1 )
Rotational axis in the tool chain { if : C a l ( D i r ) = C a l ( A x i s ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( A x i s ) ( A x i s ) if : C a l ( D i r ) C a l ( A x i s ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( A x i s ¯ ) ( A x i s ) { if : C a l ( D i r ) = C a l ( T o o l v e c t o r ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( R 2 ) ( R 2 ) if : C a l ( D i r ) C a l ( T o o l v e c t o r ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( D i r + T o o l v e c t o r ¯ ) ( R 2 )
Table 4. Mapping expressions of TTTRR-type five-axis machine tool.
Table 4. Mapping expressions of TTTRR-type five-axis machine tool.
Sensitive Linear PDGEsSensitive Angular PDGEs
Translational axis in the workpiece chain T y p ( 1 ) A x i s D i r = Δ C a l ( D i r ) ( A x i s ) T y p ( 3 ) A x i s D i r = Δ C a l ( D i r ¯ ) ( A x i s )
Translational axis in the tool chain T y p ( 1 ) A x i s D i r = Δ C a l ( D i r ) ( A x i s ) T y p ( 3 ) A x i s D i r = Δ C a l ( D i r ¯ ) ( A x i s )
Rotational axis in the tool chain When   A x i s = R 1 :
{ if : C a l ( D i r ) = C a l ( R 2 ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( R 1 ¯ ) ( R 1 ) if : C a l ( D i r ) C a l ( R 2 ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( X + Y + Z ) ( R 1 )
When   A x i s = R 2 :
{ if : C a l ( D i r ) = C a l ( R 2 ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( R 2 ) ( R 2 ) if : C a l ( D i r ) C a l ( R 2 ) , then : T y p ( 2 ) A x i s D i r = Δ C a l ( R 2 ¯ ) ( R 2 )
When   A x i s = R 1 :
{ if : C a l ( D i r ) = C a l ( R 2 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( R 1 ) ( R 1 ) if : C a l ( D i r ) C a l ( R 2 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( T o o l v e c t o r ¯ ) ( R 1 )
When   A x i s = R 2 :
{ if : C a l ( D i r ) = C a l ( R 2 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( R 1 ) ( R 2 ) if : C a l ( D i r ) C a l ( R 2 ) , then : T y p ( 4 ) A x i s D i r = Δ C a l ( X + Y + Z ) ( R 2 )
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Song, Z.; Ding, S.; Chen, Z.; Lu, Z.; Wang, Z. High-Efficient Calculation Method for Sensitive PDGEs of Five-Axis Reconfigurable Machine Tool. Machines 2021, 9, 84. https://0-doi-org.brum.beds.ac.uk/10.3390/machines9050084

AMA Style

Song Z, Ding S, Chen Z, Lu Z, Wang Z. High-Efficient Calculation Method for Sensitive PDGEs of Five-Axis Reconfigurable Machine Tool. Machines. 2021; 9(5):84. https://0-doi-org.brum.beds.ac.uk/10.3390/machines9050084

Chicago/Turabian Style

Song, Zhanqun, Shuang Ding, Zhiwei Chen, Zhongwang Lu, and Zhouzhou Wang. 2021. "High-Efficient Calculation Method for Sensitive PDGEs of Five-Axis Reconfigurable Machine Tool" Machines 9, no. 5: 84. https://0-doi-org.brum.beds.ac.uk/10.3390/machines9050084

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

Article Metrics

Back to TopTop