Error Distribution of a 5-Axis Measuring Machine Based on Sensitivity Analysis of Geometric Errors

Geometric errors are inevitably introduced into any multiaxis measuring system, and the geometric error is one of the main factors that seriously affects the measurement accuracy. ,e present work investigates the error distribution of the prototype of a 5-axis measuring machine based on sensitivity analysis of geometric errors. ,e measurement error modeling of the 5axis measuring machine is first established via the homogeneous coordinate transformation, and the Sobol global sensitivity analysis method is then employed to quantify the influence of geometric errors on the measurement result with the sensitivity index. ,e result shows that most of the angular errors are the crucial geometric errors seriously affecting the measurement result. ,ese errors are supposed to be fully considered in the accuracy design and manufacturing stages. ,e error levels of the crucial geometric errors were distributed and readjusted according to the sensitivity analysis result. Some practical approaches to distribute and improve the crucial geometric errors have been given in detail. ,e error distribution method is effective to equalize the influence of the crucial geometric errors on the measurement result as possible. ,e findings of this study provide significant meanings for the optimal design and accurate manufacturing of the 5-axis measuring machine, and the proposed method can be used to improve the measurement accuracy of the 5-axis measuring machine.


Introduction
e high machining accuracy of mechanical parts becomes increasingly essential in the fields of electronic communications, biomedicine, and aerospace [1][2][3]. e evaluation of machining accuracy requires the support of measurement technology, and the processing technology is supplementary to the measurement technology. e measurement technology is a comprehensive and interdisciplinary subject involving the optics, electronics, sensors, mechanical manufacturing, and computer technology; it provides evaluation and measurement methods for the processing technology. e rapid development of measurement technology is the driving force to increasingly promote the machining accuracy and vice versa.
Unfortunately, a measuring machine is not an ideal machine that operates without any error. e measuring machine inevitably operates under the influence of various error factors such as geometric errors, thermal error, and static force error, which lead to measurement inaccuracy. Geometric errors contribute serious influence to the measurement result of a measuring machine during the measuring process. erefore, it is crucial to reveal the relationship between geometric errors and the measurement result and to further investigate the influence of geometric errors on the measurement result. e relationship between geometric errors and the measurement result can be described by the measurement error modeling based on homogeneous coordinate transformation while the influence of geometric errors on the measurement result can be obtained using the sensitivity analysis method. e homogeneous coordinate transformation method can accurately express the position and posture between two adjacent rigid bodies; it is very powerful to establish the measurement error modeling of a measuring machine.
Wei and Chen [4] built the geometric errors of X slideway in a coordinate-measuring machine to study the effect of the geometric errors on the positioning error of the probe tip. Liu et al. [5] used homogeneous coordinate transformation to obtain the geometric error modeling of a five-axis machine tool, the geometric errors were compensated according to the geometric error modeling, and the compensation method had remarkable improvement on the machining accuracy of the machined workpiece. Yao et al. [6] obtained a geometric error modeling of an ultraprecision multiaxis machine tool and successfully found out the important error sources that contribute to the more influence on the error modeling via sensitivity analysis. Chen et al. [7] established the geometric error modeling of a five-axis ultraprecision turning machine by homogeneous coordinate transformation, and the machine tool was redesigned for optimization according to the sensitivity analysis results. Li et al. [8] described the error modeling of a five-axis machine tool based on homogeneous coordinate transformation and used a sensitivity analysis method to determine the crucial geometric error items contributing to the processing accuracy. Liu et al. [9,10] investigated the geometric errors of rotary axes for a five-axis machine tool and proposed the identification and compensation method by constructing an equivalent rotary axis. A blade machining test was employed to verify the effectiveness of the proposed method.
Sensitivity analysis is one of the most effective methods to quantify the influence of each input parameter on the output result; it is classified into local sensitivity analysis and global sensitivity analysis due to the corresponding characteristic. Local sensitivity analysis can be achieved via partial differentiation and evaluates the sensitivity index of error items at a fixed position in the workspace. In contrast, global sensitivity analysis can be performed by a variety of ways including Fourier amplitude sensitivity testing (FAST) [11], Sobol method [12], and the integration method [8]. Global sensitivity analysis represents the comprehensive sensitivity over the entire workspace. Global sensitivity analysis has become the focus of attention on the research because of its excellent performance. Cheng et al. [13] employed extended Fourier amplitude sensitivity testing to investigate the interaction among geometric error items of a 3-axis milling machine and identified the dominant geometric error items by a high-dimensional nonlinear model. Zhang et al. [14] computed the influence of geometric errors on production performance by a global sensitivity analysis method. Zhang et al. [15] utilized Sobol method to identify the significant geometric errors of a vertical three-axis milling machine center. Cheng et al. [16] introduced Sobol method into a vertical three-axis milling machine center to identify the geometric errors that most strongly affect the machining accuracy and successfully determined the vital geometric errors that contributed the most to the processing inaccuracy. Several optimization measures have been then adopted, and the improvement of machining accuracy verified the effectiveness of Sobol method.
is literature review demonstrates that the sensitivity analysis is very efficient to quantify the influence of the geometric errors on the error modeling; thus, sensitivity analysis has received lots of attention in the past time. e global sensitivity analysis method is suitable to discriminate the crucial geometric errors in the error modeling of a measuring machine or a processing machine tool. However, relatively few investigations have been performed for a 5axis measuring machine, and the influence of geometric errors on the measurement result is not well understood. To overcome this limitation, the present study employs Sobol method to perform the sensitivity analysis of geometric errors of a 5-axis measuring machine and to quantify the effect of the geometric error items on the measurement result of the 5-axis measuring machine.
In general, there are two approaches to significantly decrease the measurement error for a measuring machine: the first one is to completely eliminate the geometric errors of a measuring machine while the second one is to reasonably distribute and control the error level of each geometric error according to the sensitivity analysis result in order to equalize the effect of geometric errors on the measurement result. Obviously, the first approach is more difficult to achieve because the geometric errors are inevitably introduced in the assembly and manufacturing processes. In contrast, the second approach is a more feasible way and has been widely used in practical production, which requires to moderately reduce the error levels of some crucial geometric errors with a larger sensitivity index and to finally ensure all the sensitivity indices of geometric errors equivalent as possible. e present study aims to achieve the error distribution of the prototype of a 5-axis measuring machine based on sensitivity analysis of geometric errors. e measurement error modeling is first obtained by homogeneous coordinate transformation, Sobol method is then employed to determine the crucial geometric errors that contribute strong influence to the measurement result of a 5-axis measuring machine, and the implementation solution of error distribution in the design stage is finally discussed in detail. e homogeneous coordinate transformation method is briefly reviewed in Section 2, and the measurement error modeling of a 5-axis measuring machine is obtained. Section 3 introduces the Sobol sensitivity analysis method and its working principle. Section 4 discusses the sensitivity analysis result in detail and proposes the error distribution measures. Finally, Section 5 provides significant conclusions.

Measurement Error Modeling of the 5-Axis
Measuring Machine 2.1. e Description of the 5-Axis Measuring Machine. A 5axis measuring machine is employed to measure a workpiece designed for a special purpose. e principal feature of the workpiece is a sphere with radius R s 25 mm, and it is supported by a cylinder with radius R d 10 mm in the Z positive direction, as shown in Figure 1. Besides, five planes on the workpiece surface are parallel to the coordinate system, and the distances of these planes to the sphere center are 18 mm (plane 1 in the X positive direction), 20 mm (plane 3 in the X negative direction), 19 mm (plane 2 in the Y positive direction), 20 mm (plane 4 in the Y negative direction), and 18.5 mm (plane 5 in the Z negative direction), respectively. e required technical parameters are as follows: the surface roughness is R a ≤ 50 nm, the sphericity of the main body is ≤2 μm, the flatness of these five planes is ≤1 μm, and the perpendicularity and parallelism between planes are ≤2 μm. ese main measurement parameters are strictly required during the measuring process. Correspondingly, the proposed schematic of the 5-axis measuring machine is shown in Figure 1. e 5-axis measuring machine composes of two rotational axes (the Baxis and C-axis) and three translational axes (the X-axis, Yaxis, and Z-axis, respectively). e coordinate system is a right-hand rectangular system, and the nomenclature of each axis is defined in accordance with the rules mentioned in ISO 841 [17]. e X-axis and Z-axis are horizontally fastened to the base of the 5-axis measuring machine, the operation reference line of the X-axis is required to be perpendicular to that of the Z-axis, and the Y-axis is mounted on the top of the Z-axis and carries the C-axis in which the Y-axis moves along the vertical direction. e B-axis is installed on the Xaxis, and the axis of rotation of the B-axis is parallel to the operation reference line of the Y-axis. Analogously, the axis of rotation of the C-axis is parallel to the operation reference line of the Z-axis. A vacuum chuck is mounted to the C-axis to support the workpiece while the measurement probe holder is supported by the B-axis. Both B-axis and C-axis serve as the angular servo axis. e 5-axis measuring machine composes of eight-part rigid bodies [18], as shown in Figure 2. e base of the 5-axis measuring machine is selected as the zeroth-order body, and the remaining adjacent bodies are numbered in sequence; thus, the Z-axis, Y-axis, C-axis, workpiece, X-axis, B-axis, and the measurement probe are represented by rigid body 1 to 7 in succession, respectively. ere exist two branches in the 5-axis measuring machine: the first branch (workpiece branch) comprises the Z-axis, Y-axis, C-axis, workpiece, and the machine base while the second branch (measurement probe branch) comprises the measurement probe, B-axis, Xaxis, and the machine base; these two branches are connected by the machine base.

Geometric Errors and Characteristic Matrices.
Geometric errors of machine tools are classified into position-dependent geometric errors and position-independent geometric errors [9,10].
Position-dependent geometric errors are defined as the error motion caused by the linear movement or rotational movement of each axis. A kinematic axis introduces six geometric errors comprising linear and angular errors [19,20]. e six geometric errors of each translational axis (X-axis, Y-axis, and Z-axis) are a single positioning error, two straightness errors, and three angular errors (pitch, yaw, and roll). Similarly, the six geometric errors of each rotational axis (B-axis and C-axis) include three linear error components (a single axial error and two radial errors) and three angular errors (two tilt errors and an angular positioning error).
Position-independent geometric errors, also called location errors, are used to describe the reference lines of the axes of motion, including position errors, orientation errors, and zero position errors. For a translational axis, there exist three position-independent geometric errors, i.e., two orientation errors and the zero position error. For a rotational axis, there exist five position-independent geometric errors, i.e., two position errors, two orientation errors, and the zero angular position error (in analogy to the zero position error of a translational axis) [21][22][23][24][25].
erefore, there exist 49 geometric errors (30 positionindependent geometric errors and 19 position-independent geometric errors) in a 5-axis measuring machine in total. Parts of these 19 position-independent geometric errors can be appropriately simplified and ignored. e zero position error can be set to zero [26]. Additionally, the position errors and some orientation errors can be Contact measuring point  negligible in measurement error modeling because of the choice of the coordinate system and reference position, and there remain only seven orientation errors, which refer to the volumetric squareness errors among the five axes [7]. It means that the location errors can be ignored when the squareness errors are fully considered; therefore, it can be considered that a 5-axis measuring machine totally introduces 37 geometric error items; this simplification will not affect the analysis result and is generally accepted in the field of machine tool [7,18,27]. e influence of 37 geometric errors on the measurement result for a 5-axis measuring machine is investigated in this study, and the label from h 1 to h 37 is successively employed to stand for all 37 geometric errors for convenience, as listed in Table 1.
A 4 × 4 homogeneous transformation matrix is very powerful to express the position and posture between the two adjacent bodies in the 5-axis measuring machine. Usually, the transformation matrix between the ith-order body and the jth-order body is expressed as where j i T p represents a position transformation matrix; j i T pe represents a position error transformation matrix; j i T s represents a motion transformation matrix; and j i T se represents a motion error transformation matrix.
Correspondingly, the characteristic matrices of the adjacent bodies for the 5-axis measuring machine are given in Tables 2 and 3. p ijx , p ijy , and p ijz in position transformation matrix j i T p are the position coordinates of a rigid body i in the coordinates of its adjacent high-order rigid body j, and the corresponding position coordinates of the 5-axis measuring machine are listed in Table 4.

Measurement Error Modeling.
e measuring process can be regarded as the motion of the contact measuring point of measure probe within the workpiece coordinate system. Under the ideal situation without the influence of geometric errors, the measurement result is in accordance with the ideal trace of the workpiece. e ideal contact measuring point in the workpiece coordinate could be written as where P w � [P wx , P wy , P wz , 1] T is the position of the contact measuring point in the workpiece coordinate while P m � [P mx , P my , P mz , 1] T is the position of the contact measuring point in the measurement probe coordinate. However, the 5-axis measuring machine unavoidably operates with the effect of 37 geometric errors; these geometric errors will cause the actual trajectory between the probe and workpiece deviated from the ideal one; the actual measurement result is inaccurate relative to the nominal one, and the difference is defined as the measurement error.
us, the actual measurement result can be regarded as a superposition of the ideal measurement result of a workpiece and the measurement error caused by geometric errors. erefore, the actual contact measuring point in the workpiece coordinate can be written as e measurement error modeling caused by geometric errors can be written as where E x , E y , and E z are the components of measurement error E along the X, Y, and Z directions, respectively. Table 1, equation (4) can be given as e scale of comprehensive error modeling can be expressed as To improve the accuracy and reliability of the measurement result, the measurement error caused by geometric errors is supposed to be reduced or removed, which directly requires decreasing or eliminating the geometric error as possible. In fact, geometric errors cannot be completely eradicated; it means that the 5-axis measuring machine inevitably operates under geometric errors; thus, the measurement error always exists. e relationship between the measurement error and geometric errors is not well understood; it is of great significance to reveal the influence of geometric errors on the measurement result. Sobol method is adopted to qualify the influence using the sensitivity index, which is beneficial to achieve the error distribution of a 5axis measuring machine.

Basic Principle of Sobol Method.
Sobol method has been one of the most effective global sensitivity analysis methods over the past several decades. Compared with other methods, Sobol method is of the advantages of fast calculation, high efficiency, and easy application.
Sobol method calculates the sensitivity index of the input parameter using the ratio of the variance of the input parameter to the total variance of the error model, in which the quasi-Monte Carlo sampling method is employed to improve the convergence and calculation accuracy. Sobol method can perform the sensitivity analysis of a mathematical model with up to 50 input parameters no matter whether the model was linear or nonlinear. One of the most outstanding advantages of Sobol method is that it is capable of simultaneously evaluating the first-order sensitivity index and the high-order sensitivity index, which are helpful to reveal the interaction effect between or among the input parameters.
erefore, Sobol method can be applied to quantify the influence of geometric errors on the measurement result for a 5-axis measuring machine. e basic flow of sensitivity analysis using Sobol method is shown in Figure 3 in detail.
(1) Generate two random sampling matrices A k×n and B k×n based on the Sobol sequence and the probability distribution of the geometric errors, which is required as a pre-requisite to evaluate the sensitivity index, the Sobol sequences method is freely available in MATLAB software. e matrices A k×n and B k×n both contain k rows and n columns; k is the sampling number of each geometric error and n corresponds to the number of geometric errors. Matrix A (i) B is defined that only the i-th column of matrix A k×n is replaced by the i-th column of matrix B k×n , and all the remaining columns of matrix A k×n remain unchanged [12].
(2) Construct a mathematical model that accurately represents the actual problem, and the model is the basis of the sensitivity analysis. e mathematical model in this study is the measurement error  36 , h 37 * where δ mn delegates the linear errors, θ mn delegates the angular errors, the left subscript m delegates the kinematic axis while the right subscript n delegates the direction of linear error or the rotational axis of the angular error, and β mn delegates the squareness errors of m axis to n axis, respectively. Table 2: Position and position error transformation matrices of the 5-axis measuring machine (as referred to in Figure 2), where I represents the identity matrix.

Adjacent bodies
Position transformation matrix According to the thought of Sobol method [28,29], the Sobol model suitable for a 5-axis measuring machine can be h 2 , . . . , h 37 ) represents 37 geometric error items in a 5-axis measuring machine. According to the uniqueness of analysis of variance representation [30,31], the measurement error modeling function Y � f(h) can be further expressed by the increments series where Y 0 is the expected value of the model for input variables, , the similar regularity is suitable for the high-order items. e subitems in equation (7) are required to be orthogonal to each other.  Table 4: e position coordinates between body i and body j.
Mathematical Problems in Engineering Based on the mathematical statistics, equation (7) can be expressed in the form of variance V as After the normalization of dividing by the variance V, equation (8) is transformed into n i�1 S i + n i�1 n j>i S ij + · · · + S 1,2,...,n � 1.
e variance ratio is defined as the first-order global sensitivity index [12,32]. S i evaluates the main influence of the geometric error h i on the measurement result, which augments with the increasing importance of the geometric error. S ij accounts for the interaction effect between two inputs h i and h j when i ≠ j, similarly to high-order terms. S Ti is another commonly used sensitivity index, named the total sensitivity index, and it corresponds to the sum of the firstorder sensitivity index of h i and all the interaction indices including h i [33].

e Estimation Method of S i .
As mentioned above, Sobol method employs Sobol sequences to improve the efficiency of the estimators [34]. Two random sampling matrices A k×n and B k×n and one deduced matrix are required for evaluating the sensitivity index. e commonly used estimators for S i and S Ti are [29] where (A) m and (Β) m denote the m-th column of matrices A and B, respectively.

e Workspace and Geometric Error Intervals.
As discussed, the prototype of a 5-axis measuring machine is designed as shown in Figure 1; the stroke of X-axis, Y-axis, and Z-axis is, respectively, ±75 mm, ±50 mm, and ±75 mm, and the rotation range of C-axis is 360°while that of B-axis is from − 120°to +120°. erefore, there exist four body diagonals in the whole workspace of 150 mm × 100 mm × 150 mm, as shown in Figure 4; each diagonal of the workspace is evenly divided into six segments before conducting the sensitivity analysis. e five points along each diagonal near the workspace center are selected as test points because these points are located in the most frequently used coordinate range of the 5axis measuring machine, and four body diagonals generate 17 test points in total in the workspace, performing the sensitivity analysis at these points will provide more meaningful results. It is noted that B-axis is parallel to Y-axis; therefore, the rotation angle of B-axis is arranged synchronously with the location coordinate of Y-axis in the process of sensitivity analysis, similarly C-axis to Z-axis. Correspondingly, the coordinates of each test point are listed in Table 5.
Based on the previous experience in developing an ultraprecision machine tool, the error levels of geometric error  Figure 3: e flowchart of sensitivity analysis using Sobol method.

Mathematical Problems in Engineering
for the prototype of a 5-axis measuring machine are summarized and listed in Table 6, all the positioning errors of translational axes are ±1 μm/100 mm, all the linear errors except the positioning errors are ±0.5 μm/100 mm, and all the angular errors and squareness errors are ±5″. A Monte Carlo simulation of sensitivity analysis could involve thousands or tens of thousands of iteration calculations. With the increasing sampling number of each geometric error, the simulation accuracy and time consuming could be simultaneously increased. For the objective of utilizing the lower sampling number to attain the same simulation accuracy in a shorter time, it is necessary to balance the relationship between the sampling number and time consuming. 10000 samplings are sufficient for each geometric error under the premise of ensuring the simulation accuracy [16]; thus, two independent 10000 × 37 sampling matrixes A and B are constructed by Sobol sequences for the simulations, and the sensitivity index of geometric errors at each test point can be obtained by performing Sobol method; the sensitivity index is finally calculated using the average method. e sensitivity analysis result will be discussed in Section 4.

Results and Discussion
It is commonly accepted that the first-order sensitivity index S i represents the main influence of the geometric error h i on the measurement result. It is proved that S i is capable of adequately representing S Ti when the error levels of geometric errors are relatively small [35], which means that the interaction between any two arbitrary geometric errors can be ignored in the sensitivity analysis process, and this simplification can still provide valuable guidance without affecting the analysis result. erefore, only S i is henceforth discussed in the study. e sensitivity analysis was conducted at each test point which is listed in Table 5 to investigate the influence of geometric errors on the comprehensive error modeling and its component along the X, Y, and Z directions. e corresponding results of the sensitivity analysis are subsequently discussed in detail.
As mentioned, the sum of all the first-order sensitivity indices is regarded as 1 when the higher-order sensitivity indices are neglected. ere are 37 geometric errors in the 5axis measuring machine; thus, the sensitivity index 0.05 (approximately twice as the mean sensitivity index) is selected as the critical threshold to determine whether the geometric error provides a significant influence on the measurement result. If the sensitivity index S i of a geometric error h i is greater than 0.05, the geometric error h i can be regarded as a crucial sensitive error, and the geometric errors whose sensitivity index is below 0.05 are hereby nonsensitive error items.

Sensitivity Analysis Result of Geometric Errors.
e Xdirectional result of sensitivity index S i is displayed in Figure 5. It is obvious that h 5 and h 28 are the two outstanding geometric errors contributing enormously to the X-directional measurement error. Additionally, it is observed that the values of S i for h 6 , h 11 , h 27 , and h 29 are greater than 0.05; thus these four geometric errors are important geometric errors. In addition, h 27 and h 29 have an equivalent degree of effect on the X-directional measurement error. e remaining geometric error items have little or no influence on the measurement result and are therefore not considered as crucial geometric errors. From these results, it is concluded that h 5 , h 6 , h 11 , h 27 , h 28 , and h 29 are the crucial geometric errors that seriously impress the X-directional measurement error.
According to the labels listed in Table 1, these above six error items are θ xy , θ yy , θ zy , θ xz , θ zz , and β zx , and it is noted that the crucial geometric errors influencing strongly the Xdirectional measurement error are angular and squareness errors, i.e., three angular errors (θ xy , θ yy , and θ zy ) along the Y direction, respectively, attributed to the motion of X-axis, Yaxis, and Z-axis, two angular errors (θ xz and θ zz ) along the Z y x z Figure 4: e test points in the workspace.   Mathematical Problems in Engineering direction, respectively, attributed to the motion of X-axis and Z-axis, and one squareness error (β zx ) of X-axis to Zaxis.
As was conducted previous, the sensitivity analysis result along the Y direction is shown in Figure 6, and h 4 is the most important factor affecting the Y-directional measurement error and contributes over 40% of the total sum of S i . Besides, the sensitivity index S i of h 10 , h 13 , h 26 , h 28 , h 33, and h 36 is greater than 0.05. Based on the above information, it can be reasonably concluded that h 4 , h 10 , h 13 , h 26 , h 28 , h 33 , and h 36 are the crucial geometric errors influencing the Y-directional measurement error, and h 4 plays the most important role in the Y-directional measurement error. Accordingly, θ xx , θ yx , θ zx , θ bx , β yz , and β bx are the geometric errors that have a serious impact on the measurement result along the Y direction, i.e., four angular errors (θ xx , θ yx , θ zx , and θ bx ) along the X direction, respectively, attributed to the motion of X-axis, Y-axis, Z-axis, and B-axis, one angular error (θ xz ) along the Z direction attributed to the motion of X-axis, and two squareness errors (β yz and β bx ) of Y-axis to Z-axis and B-axis to X-axis. It is noted that the angular and squareness errors remained dominant comparing to the linear errors. Figure 7 shows the sensitivity analysis result along the Z direction, and h 26 is the most important geometric error with the sensitivity index 0.39. In addition, h 4 is another critical geometric error with the sensitivity index greater than 0.20. h 27 , h 29 , h 34 , and h 37 are considered as the critical geometric errors, and these geometric errors provide an equivalent degree of influence on the Z-directional measurement error. erefore, θ xx , θ zx , θ xy , θ by , β zx , and β bz are the crucial geometric errors contributing serious effect on the measurement result along the Z direction, i.e., two angular errors (θ xx and θ zx ) along the X direction, respectively, attributed to the motion of X-axis and Z-axis, two angular errors (θ xy and θ by ) along the Y direction, respectively, attributed to the motion of X-axis and B-axis, and two squareness errors (β zx and β bz ) of Z-axis to X-axis and B-axis to Z-axis. e sensitivity indices of geometric errors with respect to comprehensive error modeling are investigated, and the corresponding result is given in Figure 8. It is obvious that each sensitivity index of h 4 , h 5 , h 6 , h 11 , h 26 , h 27 , h 28 , and h 29 is greater than the critical threshold 0.05, which leads to the conclusion that the comprehensive error modeling is sensitive to these geometric errors including θ xx , θ xy , θ xz , θ yy , θ zx , θ zy , θ zz , and β zx . θ xx , θ xy , and θ xz represent the angular errors, respectively, along the X, Y, and Z direction caused by the motion of X-axis, analogously for θ zx , θ zy , and θ zz . θ yy is the angular error along the Y direction attributed to the motion of Y-axis, and β zx is the squareness error of Z-axis to X-axis.
From these results, it is noted that most of the crucial geometric errors in Figure 8 are consistent with those in Figures 5-7; this indicates that these crucial geometric errors, whose sensitivity index is greater than the critical threshold 0.05, have a serious effect on both the comprehensive error modeling and its components along the X, Y, and Z direction. erefore, the crucial geometric errors identified by sensitivity analysis are supposed to be fully considered in the design and manufacturing stages of the 5-axis measuring machine, which is helpful to reduce the influence of geometric error on the measurement result. erefore, the error levels of these crucial geometric errors should be improved and distributed as possible in order to equalize their influence on the measurement result. e geometric errors, as mentioned before, are difficult to completely eliminate because of the assembly and Mathematical Problems in Engineering manufacturing processes; it is a more feasible and reasonable approach of moderately distributing the error levels of the crucial geometric errors with larger sensitivity index in order to decrease the influence of some crucial geometric errors on measurement error. is method is employed in the study to improve the performance of the 5axis measuring machine, and the error levels of some crucial geometric errors are distributed and readjusted according to the sensitivity analysis results, as listed in Table 7; the corresponding reduction range is from 20% to 70%, and the error levels of the remaining geometric errors remain constant.
After the distribution of geometric errors, the sensitivity analysis is re-executed to investigate the efficiency of the proposed error distribution method, and the result is available in the following subsection.

Validation of the Effectiveness of Error Distribution.
e sensitivity analysis is re-executed after the error levels are distributed, and Figures 9-12, respectively, show the sensitivity analysis result along the X, Y, and Z directions and to comprehensive error modeling after error distribution.
Compared to Figure 5, it is obviously observed in Figure 9 that the sensitivity index S 5 rapidly decreases, the   maximum sensitivity index along the X direction is within 0.15 after the error levels are distributed, and the sensitivity indices of the remaining crucial geometric errors show a slight fluctuation and tend to be equivalent. As shown in Figure 10, the sensitivity index of h 4 dramatically drops to 0.168 after the error levels are distributed, the maximum sensitivity index along the Y direction is within 0.2, and the sensitivity indices of the crucial geometric errors are nearly equivalent. e sensitivity indices along the Z direction and with respect to comprehensive error modeling are, respectively, shown in Figures 11 and 12. It is noted that all the sensitivity indices rapidly reduced to less than 0.2 after error distribution; it means that the error distribution method is effective to reduce the influence of the crucial geometric errors on the measurement result, which leads to the similar regularity obtained from Figures 9 and 10.
From the sensitivity analysis result after error distribution, it is reasonably concluded that error distribution is effective to reduce the influence of the crucial geometric errors with large sensitivity index on measurement error. All the sensitivity indices rapidly decrease to near or below 0.2, the crucial geometric errors having the most important influence along the X, Y, and Z direction disappear, and the sensitivity indices of the crucial geometric errors tend to be nearly equivalent. e sensitivity indices of the crucial geometric errors tend to be at the approximately equivalent level by the reasonable geometric error distribution and adjustment. e above results of sensitivity analysis are beneficial for the error distribution of the prototype of a 5-axis measuring machine, and the results provide significant meaning to the optimal design and accurate manufacturing of the 5-axis measuring machine, which is helpful to equalize the influence of the crucial geometric errors as possible.

Error Distribution and Adjustment
Methods. Based on the detailed discussion of the forgoing sensitivity analysis result, the crucial geometric errors for the prototype of a 5axis measuring machine are three angular errors (θ xx , θ xy , and θ xz ) caused by the motion of X-axis, two angular errors (θ yx and θ yy ) caused by the motion of Y-axis, three angular errors (θ zx , θ zy , and θ zz ) caused by the motion of Z-axis, and an angular error θ by caused by the motion of B-axis, as well as two squareness errors (β zx and β yz ) of Z-axis to X-axis and Y-axis to Z-axis. It is summarized that the crucial geometric errors are composed of the angular errors caused by the translational axes, the angular error caused by the motion of B-axis, and the squareness errors.
Due to the characteristics of the crucial geometric errors, some practical measures can be taken to distribute and adjust the levels of the crucial geometric errors.

e Angular Errors of the Translational and Rotational Axes Can Be Reduced by Improving the Geometric Imperfections of the Manufacturing and Assembly of Machine
Components.
e translational axes (X-axis, Y-axis, and Zaxis) and B-axis in the 5-axis measuring machine are in the form of orifice-compensated hydrostatic bearing, as shown in Figure 13.
Improving the machining accuracy (dimensional accuracy and form accuracy) of the matching components is an effective approach to guarantee the matching clearance and oil film thickness. For a translational axis, the poor surface flatnesses for slideway carriages and slideway base are the main factors causing angular errors. For a rotational axis, the roundness errors and cylindricity errors of the shaft and shaft sleeve lead to the angular errors. Dimensional accuracy of the matching surfaces affects the oil film thickness for both a translational axis and a rotational axis, and the excessive oil film thickness causes the increasing angular errors. erefore, these above form errors and dimensional errors need to be strictly monitored and corrected in the precision machining stage. Although the assembly accuracy is strongly dependent of the machining accuracy of the matching components, the assembly error induced in the assembly process can also cause the angular error, and the assembly inaccuracy should be fully considered to ensure the closure of the dimensional accuracy chain, which is beneficial to improve the angular errors of the translational and rotational axes.
Additionally, the increasing oil supply pressure is beneficial to reduce the angular errors because of the error averaging effect of the oil film. e oil film acts as a lubricant when a translational axis or a rotational axis is running, and the increasing oil supply pressure can directly improve the oil film stiffness and the error averaging effect, which is effective for decreasing the angular errors caused by the imperfections of the matching components [36].

e Squareness Errors (β zx and β yz ) Can Be Improved by Adjusting the Corresponding Spatial Relationship among Axes and Software Correction Method.
e squareness error of axis of motion with respect to each other is generally affected by the orientation of the reference line of these axes [26], and the squareness errors can be calibrated by the laser interferometer and double ball bar in the assembly process of machine components; therefore, adjusting the alignment and orientation of the reference line of these axes is helpful to enhance the squareness errors in the assembly process of the 5-axis measuring machine.
To overcome the accuracy limitation in the mechanical adjustment, the software correction method can be used to further improve the squareness errors. e software correction method uses UMAC (Universal Motion and Automation Controller) to redefine the coordinate system and correct the motor parameters in order to orthogonally compensate for the offset caused by the squareness errors.
e UMAC provides the function of creating a coordinate system, defining an axis system, and configuring the motor parameters. is function is employed in this study to perform orthogonal correction of the squareness errors β zx and β yz . e schematic of the squareness errors β zx and β yz is shown in Figure 14; the squareness error β zx causes the xdirectional offset relative to the nominal position when Z-axis is moving, and the squareness error β yz causes the zdirectional offset relative to the nominal position when Yaxis is moving. e offset can be approximately represented by the reference line.
As the pulse signal is preferred by UMAC, the position of each axis should be first transformed into the corresponding number of pulse counts that can be recognized by UMAC. e linear encoder with the linear resolution of 5 nm is used as the position feedback unit in the 5-axis measuring machine, so the displacement of 1 mm corresponds to 200,000 pulse counts. e squareness errors β zx and β yz can be corrected by loading the following codes into UMAC: Undefine All; % Clear all axis-definitions in the coordinate system. &1; % Use the symbol "&" to define a coordinate system and label it as No. 1. #1->200000Y; % Assign motor 1 to Y-axis, and convert the Y-axis position into the number of pulse counts. #2->200000Z + 200000 * Y * tan(β yz ); % Assign motor 2 to Z-axis, convert the Z-axis position into the number of pulse counts, correct the z-directional offset caused by the squareness errors β yz . #3->200000X + 200000 * Z * tan(β zx ); % Assign motor 3 to X-axis, convert the X-axis position into the number of pulse counts, correct the x-directional offset caused by the squareness errors β zx .
After the correction code has been stored in UMAC, the offset caused by the squareness errors can be compensated via a software correction method to improve the squareness errors β zx and β yz . e 5-axis measuring machine could be optimal designed and accurate manufactured for improving the efficiency and performance, and an experimental prototype is fabricated, as shown in Figure 15. e future work will mainly focus on the research of software development, measurement method, and data processing method.

Conclusions
e present work focuses on the error distribution of a 5-axis measuring machine based on sensitivity analysis of geometric errors. e objective of this paper is to equalize the influence of the crucial geometric errors on the measurement result by distributing reasonably the error levels of the crucial geometric errors. is ensures that the measurement accuracy of a 5-axis measuring machine is no longer strictly limited to some geometric errors with larger sensitivity index, which provides significant meaning to the optimal design and accurate manufacturing of the 5-axis measuring machine. e homogeneous coordinate transformation method was first employed to establish the measurement error modeling of the 5-axis measuring machine. e Sobol global sensitivity analysis method was then used to investigate the sensitivity indices of geometric errors and to  further quantify the influence of geometric errors on the measurement result. Finally, the error levels of the crucial geometric errors were distributed and readjusted according to the sensitivity analysis result, and the sensitivity analysis result after error distribution demonstrates that the error distribution is effective to equalize the influence of the crucial geometric errors on the measurement result. Some practical approaches to distribute and improve the crucial geometric errors were given in detail.
A summary of the significant conclusions are presented as follows: (1) e sensitivity analysis can investigate the influence of geometric errors on the measurement result in a 5axis measuring machine, which is beneficial for the accuracy design and manufacture of a 5-axis measuring machine. e sensitivity analysis result shows that most of the crucial geometric errors seriously affecting the measurement error are the angular errors caused by the translational axis and B-axis, as well as the squareness errors β zx and β yz . e crucial geometric errors are supposed to pay more attention in the design and manufacturing stage. (2) e error distribution method is proposed to equalize the influence of the crucial geometric errors on the measurement result, and the effectiveness of the proposed error distribution method is validated by the re-executed sensitivity analysis; some practical measures to distribute and adjust the levels of the crucial geometric errors are provided in detail. (3) e error distribution method is effective to reduce the influence of the crucial geometric errors on the measurement result, and this finding is significant for the optimal design and accurate manufacturing of the 5-axis measuring machine.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.