• Views 1,739
• Citations 6
• ePub 21
• PDF 965
`Mathematical Problems in EngineeringVolume 2014, Article ID 745250, 8 pageshttp://dx.doi.org/10.1155/2014/745250`
Research Article

## Error Modeling and Sensitivity Analysis of a Five-Axis Machine Tool

1Key Laboratory of Mechanism Theory and Equipment Design of Ministry of Education, Tianjin University, Tianjin 300072, China

2Beijing Precision Machinery & Engineering Research Co., Ltd., Beijing 101312, China

Received 20 October 2013; Revised 27 February 2014; Accepted 27 February 2014; Published 29 May 2014

Academic Editor: Miguel A. F. Sanjuán

Copyright © 2014 Wenjie Tian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Geometric error modeling and its sensitivity analysis are carried out in this paper, which is helpful for precision design of machine tools. Screw theory and rigid body kinematics are used to establish the error model of an RRTTT-type five-axis machine tool, which enables the source errors affecting the compensable and uncompensable pose accuracy of the machine tool to be explicitly separated, thereby providing designers and/or field engineers with an informative guideline for the accuracy improvement by suitable measures, that is, component tolerancing in design, manufacturing, and assembly processes, and error compensation. The sensitivity analysis method is proposed, and the sensitivities of compensable and uncompensable pose accuracies are analyzed. The analysis results will be used for the precision design of the machine tool.

#### 1. Introduction

Geometric accuracy is a crucially important performance factor for machine tools, especially under circumstances where relatively high precision is one of the basic requirements [1]. There are two ways to improve the geometric accuracy of machine tools: design and manufacture for precision and error compensation. Both require a parametric model that relates the geometric source errors to the pose accuracy of the cutting tool relative to the workpiece. Theoretically, the effects of source errors on pose accuracy of 3-, 4-, or 5-axis machine tools cannot fully be compensated by software, and only those pose errors corresponding to the permissible motion types are compensable by means of error compensation. Therefore, a comprehensive error model is essential in order to distinguish the source errors affecting compensable pose errors from those affecting uncompensable ones. Only then can suitable measures be adopted for accuracy improvement [2, 3] via component tolerance design, manufacturing, and assembly, as well as error compensation.

Of all types of source errors, geometric source errors have the largest influences on the workpiece accuracy, which should be considered first in design. Over the past decades, there has been a great deal of intensive research into error modeling of machine tools. Homogeneous transformation matrix (HTM) or Denavit-Hartenberg (D-H) transformation matrix is widely used to describe the rigid body kinematics. Kiridena and Ferreira [4] classified 5-axis machine tools into TTTRR, RTTTR, and RRTTT systems and used the D-H convention to develop kinematic models for each of these three machine types. However, their model considered only five parametric errors (one positioning error for each axis). Srivastava et al. [5] established a combined error model, which considered both geometric errors and thermal errors. For error compensation of a 5-axis machine tool, Jha and Kumar. [6] established the volumetric error model with D-H method. Lin and Shen [7] proposed a new matrix summation method; the kinematic equation was converted into six components to give each component a clear physical meaning. However, the physical meaning of geometric source error was not clear. Yang [8] dedicated one chapter of his dissertation to the formulation of a generalized 5D error synthesis model. However, he defined only 27 geometric error components, which is incomplete. Similar work has been done for the purpose of error modeling of machine tools [914]. In the past decades, most studies mainly aimed at error compensation, and little attention has been paid to the important design question of separating the source errors into two subsets associated with the compensable and uncompensable pose accuracy.

Concerning the error sensitivity analysis, Huang et al. [15] established a statistical sensitivity model and analyzed the sensitivity of uncompensable pose accuracy with respect to the geometric source errors for a 3-DOF parallel kinematic machine with parallelogram struts. Chanal et al. [16] and Fan et al. [17] did the sensitivity analysis of another two PKMs with similar method. And many similar work [18, 19] has been done. Chen et al. [20] analyzed the sensitivity of the volumetric error regarding 37 source error components for designing a five-axis machine tool. However, only three special configurations were chosen to perform the analysis, which were probably insufficient.

In this paper, after Section 1 has briefly addressed current challenges in error modeling and sensitivity analysis of machine tools, Section 2 introduces two kinds of geometric source errors within the machine tools. In Section 3, a systematic model is established, which formulates the linear map between the pose error twist and geometric source errors within a machine tool using the homogeneous transformation matrix method. Section 4 gives an error model that fully clarifies which source errors affect the compensable and uncompensable pose accuracy. Finally, Section 5 proposes a new method for sensitivity analysis of multiaxis machine tools and takes a typical five-axis machine tool as an example to illustrate the generality and effectiveness of this approach before conclusions are drawn in Section 6.

#### 2. Geometric Source Errors of a Five-Axis Machine Tool

##### 2.1. Machine Tool Structure

The schematic diagram of a five-axis machine tool (RRTTT-type) is shown in Figure 1. The machine tool has three translational axes and two rotational axes, that is, -axis, -axis, -axis, -axis, and -axis, respectively. The tool system is supported by the -axis and the workpiece is supported by the -axis. The -axis is supported by the -, - and -axis sequentially.

Figure 1: Schematic diagram of a 5-axis machine tool.
##### 2.2. Geometric Source Errors

According to the definitions proposed in [20], the geometric error of the machine tool refers to the error of individual axes and those between axes, such as linear/angular positioning errors, straightness errors, orientational errors, offset errors, and squareness errors. Every translational/rotational axis has six position-dependent geometric source errors. All position-dependent errors of this machine tool are shown in Table 1. Position-independent errors of a machine tool refer to the location errors of an axis from its nominal pose [20]. In this paper, we use to express the squareness error between -axis and -axis. For the machine tool shown in Figure 1, with the method presented in [20], the real motion direction of -axis is selected as the reference direction; so -axis has no location errors. The plane through -axis and -axis is selected as the reference plane; so -axis has one squareness error , and -axis has two squareness errors and . Similarly, -axis has two squareness errors and , and -axis has two squareness errors and . , , and are the linear offset errors of -axis with respect to -axis in , , and directions, respectively. is the offset error of the -axis with respect to -axis in direction, since the two rotary axes may not intersect perfectly.

Table 1: Geometric source errors of 5-axis machine tools.

#### 3. Error Modeling

Error modeling aims to establish a map from geometric source errors to the pose errors of machine tool, which is the common premise of precision design and error compensation. All kinds of machine tools have two structural chains which contain several components connected in series by prismatic and rotational joints. One is from machine bed to workpiece, and the other is from machine bed to cutting tool, which are represented by workpiece kinematic chain and cutting tool kinematic chain, respectively. Thus, establishing the error model of kinematic chain is the foundation of the whole machine tool error model.

##### 3.1. Error Modeling for Workpiece Kinematic Chain

Figure 2(a) shows the workpiece kinematic chain of the five-axis machine tool. This chain has four joints, two translational joints and two rotational joints, and, here, we use to present the -, -, -, and -axis sequentially. Modeling the chain error uses two global reference frames and : is located at point on the machine bed, while is located at point on the workpiece and remains parallel to . In order to describe the source errors within the chain and relate them to pose accuracy of workpiece, body fixed frame is placed at point , and body fixed frames are placed on the th joint (). As previously mentioned, the source errors can be divided into two categories: position-dependent errors, that is, the location errors of relative to , and position-independent errors, that is, the joint motion errors along/about axes measured in .

Figure 2: Two kinematic chains of the five-axis machine tool: (a) workpiece kinematic chain; (b) cutting tool kinematic chain.

At the nominal configuration, the pose (position and orientation) of the workpiece measured in can be expressed by successive homogeneous transformations. Consider where is the homogeneous transformation matrix of , with respect to , is the position vector of the origin of , and is its orientation matrix.

Consider initially that the only location errors are of relative to and the only joint motion errors are of the th joint. The actual workpiece pose measured in is then approximately expressed by where is the fourth-order identity matrix, denotes the pose error matrix of the workpiece measured in , and is the location error matrix of relative to expressed by Here, , , and are the location error components of relative to along (about) three orthogonal axes of , is the skew-symmetric matrix of vector relating to its vector product, and is the joint motion error matrix measured in and has the form Here, and are the movement error vectors along and about the axes of and is the linear or angular displacement of the th joint. The elements of and have been listed in Table 1.

Now, by taking similar account of all location errors and joint motion errors within every joint of the kinematic chain, the actual workpiece pose measured in is then expressed by

Expanding (6), ignoring higher-order terms in the errors, and applying appropriate adjoint transformations lead to a vector representation of the linearized pose error evaluated in . Consider where is the adjoint transformation matrix of relative to and , , and , obtained via identification mapping [21], are the vector representations of , , and , respectively.

Finally, premultiplying both sides of (7) with results in the pose error twist of the workpiece evaluated in : where is the adjoint transformation matrix of relative to .

##### 3.2. Error Modeling for Cutting Tool Kinematic Chain

Figure 2(b) shows the cutting tool kinematic chain of the five-axis machine tool, and this chain has only one joint. Similarly, global reference frame is located at point on the tool tip. The pose error twist of the tool evaluated in can be expressed as where and denotes the -axis.

##### 3.3. Error Modeling for Machine Tool

For the two kinematic chains, since the reference point of the tool and workpiece is the same one and the pose error twists of the tool and workpiece are both evaluated in , the pose error twist of machine tool can be written as

In (10), the pose error twist of the machine tool is represented as the linear combination of the two kinds of geometric source errors, joint motion errors, and location errors. Therefore, the mapping from each source error to the pose accuracy can directly be obtained with explicit physical meaning. can also be written in matrix form: where

In (11), all geometric source errors are written in the vector form as follows:

Equation (11) can also be written in the following form: where and are the positional and orientational error vectors measured in , respectively.

#### 4. Error Separation

Theoretically, the effects of source errors on pose accuracy of 5-axis machine tools cannot fully be compensated by software, and only those pose errors corresponding to the permissible motion types are compensable by means of error compensation. Therefore, it is essential to distinguish the source errors affecting uncompensable pose errors from those affecting compensable ones. Only then can suitable measures be adopted for accuracy improvement via component tolerance design, manufacturing, and assembly, as well as error compensation.

In order to identify the uncompensable source errors from , we should first ascertain the directions of uncompensable pose error of machine tool at any configuration. From the view of motion and constrain, uncompensable directions are the constrained directions of machine tool system. Since they are not parallel to any actuated joint’s axis, in other words, the actuating forces or torques of machine tool cannot do work in these uncompensable directions, the uncompensable pose errors cannot be compensated by modifying the ideal output of actuated joints.

Generally, the position errors of cutting tool with respect to workpiece can be compensated by the three mutually perpendicular translational axes. Therefore, only the orientational errors need to be analyzed. Figure 3 shows the direction of uncompensable orientational error of the 5-axis machine tool (only the tilting rotary table is drawn). and denote the unit vectors of rotary axes and , and . Obviously, the orientational error about direction is uncompensable.

Figure 3: Direction of uncompensable orientational error of 5-axis machine tool.

According to (14), we have

Taking dot products on both sides of (16) with uncompensable direction vector leads to where is the th component in , is the coefficient of , is the number of geometric source errors in , is the th component in , and is the element of row and column in . The geometric source error whose coefficient will affect the uncompensable orientational error about direction .

Let . , , and are the orientational errors about , , and , respectively. With the method mentioned above, all the uncompensable source errors can be identified. Thus, (14) can be rewritten in terms of the compensable and uncompensable degrees of freedom such that where denotes the 5-dimensional compensable pose error vector, denotes the 1-dimensional uncompensable pose error vector, and denotes the compensable (uncompensable) geometric source error vector. For the 5-axis machine tool analyzed in this paper, , and , which should be eliminated or at least minimized during the manufacturing and assembly processes.

#### 5. Sensitivity Analysis

After the error model is established and the compensable/uncompensable source errors are separated, a probability global sensitivity model is proposed. And, then, the sensitivities of compensable/uncompensable pose accuracies of the machine tool with respect to geometric source errors are analyzed, respectively.

##### 5.1. Probability Model

In order to investigate the influences of the geometric errors on the volumetric accuracy in a statistical sense, a probability model is needed. Taking the norm on both sides of (14) yields where is the th element in , is the element of row and column in , and is the number of geometric errors. Assume that all elements in are independent and of zero mean. This allows the mean value of to be zero; that is, , where is the statistical expectation. As a result, the variance of can be expressed by

Thus, the standard deviation of \$ can be represented by where is defined as the volumetric sensitivity of with respect to , representing the standard deviation of caused by the unit standard deviation of . Note that varies with the configuration of the machine tool; therefore, its mean value can be employed as an index to evaluate the influence of on the pose accuracy in the overall workspace. This index is defined as where is the volume of the prescribed workspace.

According to the sensitivity analysis model proposed above, the sensitivities of the elements of with respect to the elements of can be expressed as

##### 5.2. Sensitivity Analysis

Utilizing the indices given by (23) and (24), the sensitivity analysis is carried out for the machine tool, whose prescribed workspace is listed in Table 2. The structure parameters are shown in Table 3. It should be noted that the parameter vector , the position vector of cutting point measured in frame , changes with different machine configuration. Therefore, it should be calculated at every given configuration. Table 4 shows the numbers for corresponding error components [20].

Table 2: Workspace of the machine tool.
Table 3: Structure parameters of the machine tool.
Table 4: Number of error components.

In the processes of sensitivity analysis, all error components are considered as the same value, that is, μm for position errors and μm/m for angular errors [20]. The sensitivity analysis results of (23) and (24) are shown in Figures 4 and 5, respectively. It can be found that linear source errors only have influences on position errors of the tool relative to the workpiece, while angular source errors having influences on both position and orientation errors. In Figure 4, we can easily find the key error components for each compensable direction. For example, the error components 5, 6, 11, 12, 24, 30, 32, 33, 34, and 35 have more influences on the orientational error than other error components. Figure 5 shows the sensitivity of uncompensable orientational error with respect to every error component in . The sensitivities of and are 1, since the direction of and always coincides with the uncompensable direction . The other components in have similar value of 0.6~0.7. The above sensitivity analysis results can be used for the accuracy design of the machine tool, and the key error components should be paid more attention to in manufacture and assembly processes.

Figure 4: Sensitivities of with respect to .
Figure 5: Sensitivities of with respect to .

#### 6. Conclusions

The proposed error modeling methodology in this paper allows the geometric source errors affecting the compensable and uncompensable pose accuracy of the five-axis machine tool to be identified in an explicit manner. The merit of this approach is that it can provide designers and/or field engineers with an informative guideline to adopt suitable measures for improving accuracy. The model has clear physical meanings and shows higher computational efficiency in numerical simulations without considering the nonlinear effects.

The sensitivities of compensable and uncompensable pose accuracies have been analyzed, which will be helpful in the accuracy design of a machine tool. It will enhance the understanding of which error components are more important or influence the machining target, and these error components should be paid more consideration to in design processes.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

This work was supported by a Grant from the National High Technology Research and Development Program of China (no. 2012AA040701).

#### References

1. R. Ramesh, M. A. Mannan, and A. N. Poo, “Error compensation in machine tools—a review. Part I. Geometric, cutting-force induced and fixture-dependent errors,” International Journal of Machine Tools and Manufacture, vol. 40, no. 9, pp. 1235–1256, 2000.
2. H. Schwenke, W. Knapp, H. Haitjema, A. Weckenmann, R. Schmitt, and F. Delbressine, “Geometric error measurement and compensation of machines-An update,” CIRP Annals—Manufacturing Technology, vol. 57, no. 2, pp. 660–675, 2008.
3. W. Tian, W. Gao, D. Zhang et al., “A general approach for error modeling of machine tools,” International Journal of Machine Tools and Manufacture, vol. 79, pp. 17–23, 2014.
4. V. Kiridena and P. M. Ferreira, “Mapping the effects of positioning errors on the volumetric accuracy of five-axis CNC machine tools,” International Journal of Machine Tools and Manufacture, vol. 33, no. 3, pp. 417–437, 1993.
5. A. K. Srivastava, S. C. Veldhuis, and M. A. Elbestawit, “Modelling geometric and thermal errors in a five-axis cnc machine tool,” International Journal of Machine Tools and Manufacture, vol. 35, no. 9, pp. 1321–1337, 1995.
6. B. K. Jha and A. Kumar, “Analysis of geometric errors associated with five-axis machining centre in improving the quality of cam profile,” International Journal of Machine Tools and Manufacture, vol. 43, no. 6, pp. 629–636, 2003.
7. Y. Lin and Y. Shen, “Modelling of five-axis machine tool metrology models using the matrix summation approach,” International Journal of Advanced Manufacturing Technology, vol. 21, no. 4, pp. 243–248, 2003.
8. S. H. Yang, Real-time compensation for geometric, thermal, and cutting force induced errors in machining tools [Ph.D. thesis], University of Michigan, Ann Arbor, Mich, USA, 1996.
9. Y. Y. Hsu and S. S. Wang, “A new compensation method for geometry errors of five-axis machine tools,” International Journal of Machine Tools and Manufacture, vol. 47, no. 2, pp. 352–360, 2007.
10. J. H. Lee, Y. Liu, and S. H. Yang, “Accuracy improvement of miniaturized machine tool: Geometric error modeling and compensation,” International Journal of Machine Tools and Manufacture, vol. 46, no. 12-13, pp. 1508–1516, 2006.
11. K. G. Ahn and D. W. Cho, “Analysis of the volumetric error uncertainty of a three-axis machine tool by beta distribution,” International Journal of Machine Tools and Manufacture, vol. 40, no. 15, pp. 2235–2248, 2000.
12. M. Rahman, J. Heikkala, and K. Lappalainen, “Modeling, measurement and error compensation of multi-axis machine tools. Part I. Theory,” International Journal of Machine Tools and Manufacture, vol. 40, no. 10, pp. 1535–1546, 2000.
13. A. C. Okafor and Y. M. Ertekin, “Derivation of machine tool error models and error compensation procedure for three axes vertical machining center using rigid body kinematics,” International Journal of Machine Tools and Manufacture, vol. 40, no. 8, pp. 1199–1213, 2000.
14. S. Zhu, G. Ding, S. Qin, J. Lei, L. Zhuang, and K. Yan, “Integrated geometric error modeling, identification and compensation of CNC machine tools,” International Journal of Machine Tools and Manufacture, vol. 52, no. 1, pp. 24–29, 2012.
15. T. Huang, D. J. Whitehouse, and D. G. Chetwynd, “A unified error model for tolerance design, assembly and error compensation of 3-DOF parallel kinematic machines with parallelogram struts,” CIRP Annals—Manufacturing Technology, vol. 51, no. 1, pp. 297–301, 2002.
16. H. Chanal, F. Paccot, and E. Duc, “Sensitivity analysis of an overconstrained parallel structure machine tool, the Tripteor X7,” Applied Mechanics and Materials, vol. 162, pp. 394–402, 2012.
17. K. C. Fan, H. Wang, J. W. Zhao, and T. H. Chang, “Sensitivity analysis of the 3-PRS parallel kinematic spindle platform of a serial-parallel machine tool,” International Journal of Machine Tools and Manufacture, vol. 43, no. 15, pp. 1561–1569, 2003.
18. S. M. Wang and J. J. Lin, “On-machine volumetric-error measurement and compensation methods for micro machine tools,” International Journal of Precision Engineering and Manufacturing, vol. 14, no. 6, pp. 989–994, 2013.
19. A. J. Patel and K. F. Ehmann, “Volumetric error analysis of a Stewart platform-based machine tool,” CIRP Annals—Manufacturing Technology, vol. 46, no. 1, pp. 287–290, 1997.
20. G. Chen, Y. Liang, Y. Sun et al., “Volumetric error modeling and sensitivity analysis for designing a five-axis ultra-precision machine tool,” The International Journal of Advanced Manufacturing Technology, vol. 68, no. 9–12, pp. 2525–2534, 2013.
21. L. W. Tsai, Robot Analysis: The Mechanics of Serial and Parallel Manipulators, John Wiley & Sons, New York, NY, USA, 1999.