Development of New Measurement System of Errors in the Multiaxial Machine Tool for an Active Compensation

Error compensation techniques have been widely applied to improve multiaxis machine accuracy. However, due to the lack of reliable instrumentation for direct and overall measurements, all the compensation methods are based on offline measurements of each error component separately. The results of these measurements are static in nature and can only reflect the conditions at the moment of measurement. These results are not representative under real working conditions because of disturbances from load deformations, thermal distortions, and dynamic perturbations. This present approach involves the development of a new measurement system capable of dynamically evaluating the errors according to the six degrees of freedom. The developed system allows the generation of useful data that cover all machine states regardless of the operating conditions.The obtainedmeasurements can be used to evaluate the performance of the machine, calibration, and real time compensation of errors. This system is able to perform dynamicmeasurements reflecting the global accuracy of themachine tool without a long and expensive analysis of various error sources contribution. Finally, the system exhibits compatible metrological characteristics with high precision applications.


Introduction
The contribution of machine tools in the development of several industrial fields is demonstrated.However, during recent decades, the economic environment has forced the industry to new standards of performance in terms of quality, productivity, costs, and production time.These market constraints have evolved much faster than the performance of existing machines.The technology behind these machines is based on concepts dating back several decades and is almost to the limit of their potential.This relative stagnation is explained by the fact that technological developments in this sector are more motivated by the needs of users than by the manufacturer's initiatives.This situation highlights the need to develop strategies to respond to these new challenges by improving the performance of the machines and making them evolved to a higher level of technology [1][2][3].
The loss of precision in the machine tools is due to the geometrical imperfections of the mechanical structure of the machine or follows the modification of this structure under static, thermal, and dynamic stresses.Even if some of these imperfections can be reduced by the improvement of machine design, it is nevertheless difficult to eliminate completely their effects.The technological limitations and the costs associated with this solution led to the introduction of the concept of errors compensation [4][5][6][7][8].However, the compensation techniques such as they are practised imply the need for measuring, under typical operating and environmental conditions, all the components of the error, before combining them through mathematical models for then generating the data necessary to make the correction in real time [9,10].Such a step, although it makes possible the correction of a great proportion of the systematic error, does not make it possible to evaluate the random components of the error.Moreover, several other limitations largely compromise the potential of this approach.Among these limitations are (i) the importance of the deadlines and the costs necessary for the remote measurement of the errors under various operating conditions; (ii) the difficulty of covering all sources of errors being able to have effects on the precision of the machine; (iii) the need for carrying out the periodic calibration of the machine whose geometry undergoes continuous variations because of wear, tiredness, and other static and dynamic effects; (iv) the errors being evaluated under static conditions, with the variations they undergo under dynamic conditions not being taken into account; (v) the errors being measured individually and the mechanisms of interaction between the various sources of errors not considered.To circumvent these limitations, it becomes necessary to develop a dynamic and integrated approach of measurement able to evaluate, in real time, the individual components of the errors associated with the geometrical, thermal, and dynamic effects.This step makes it possible to provide a reliable representation of the machine state for robust and effective adaptive compensation [11][12][13][14][15][16].
The objective of this paper is to develop a measurement system in real time to identify simultaneously all components of the volumetric error in the multiaxial machine tools.Indeed, the most recent advances in instrumentation have made possible the development of techniques for measuring performance both in terms of accuracy and in terms of speed and cost.The prototype is able to monitor and quantify the components of the error generated when moving along a main axis.Based on an optoelectronics structure, the proposed approach makes it possible to evaluate simultaneously the errors according to 5 of the 6 degrees of freedom of each machine head, regardless of the operating conditions.Moreover, this approach makes it possible to simplify in a substantial way the procedures of calibration and metrological performance evaluation of the machine tools.Compared with reference measurement techniques, this approach presents some significant advantages in terms of precision and in terms of manipulation speed and cost.

Measurement System Configuration
Recent advances in the instrumentation field made possible the development of new, powerful measurement techniques.The measurement approach, suggested during this research, is based partly on these innovations, combining the optical and electronic components to carry out a high precision system of measurement.A diagram of the device of measurement is presented in Figure 1.The basic configuration is made up of two parts: a fixed plate constituting the measurement reference and a mobile plate directly attached to the moved mechanism.The fixed part supports a He-Ne laser source, three semireflective mirrors, and three photodetectors.The moving part comprises a retroreflector and 3 flat mirrors.The reference beam generated by the Ne-He laser source is divided into 3 distinct beams.The first one is projected on flat mirror FM 1 and then reflected to photodetector PD 1 .The second beam is projected on retroreflector RR and then directed to photodetector PD 2 .The third beam is reflected to photodetector PD 3 using 2 flat mirrors FM 2 and FM 3 .Any deviation of the mobile plate compared to the reference would cause variations of beam position on the photodetectors.These variations are obtained starting from the output signals of the photodetectors connected to an electronic circuit for Y X Z

Reference plate
Laser source Mobile plate conditioning and preprocessing.A computer is used for acquisition, treatment, and exploitation of various signals.
The He-Ne laser presents low-cost solution with good repeatability.The separators of the nonpolarized laser beam prevent polarization of the incident beam.Their low absorption allows these dielectric spacers to separate the beams by 50%.They are made from a substrate in high precision, having a thermal expansion factor of the order of 7.10 −6∘ C −1 , to minimize the distortion of the wavelength.Moreover, the inclination of the substrate avoids any possible back reflection to the laser.The position sensors are made of silicon photodiodes in a single crystal with a quadruple electrode.Compatible with low-power laser beams and other infrared beams, they allow obtaining accurate information on the linear positioning of the beam in both horizontal and vertical directions.
It is important to note the degrees of freedom principle.In a linear constraint, displacement along a known axis leads to the elimination of 5 of the 6 degrees of freedom of the moving part compared to the guidance support.The movement in the direction of the remaining degree of freedom is controlled with precision using a suitable control device.A linear displacement of the moving part along a specific axis involves deviations according to three directions: three translations (positioning and straightness errors) and three rotations (yaw, pitch, and roll).Consequently, the three beams obtained from the principal beam must be aligned in the direction of the moving axis so that the considered beams are received by the photodetectors.The principle and the method of calculation of the position of each beam on a photodetector are presented in Figure 2. The beam position is calculated according to the photocurrents from the voltages   ,   ,   , and   .The positions in the horizontal (  ) and  vertical (  ) planes are calculated starting from these exits according to preestablished relations ( The photodetector PD 1 receives the beam reflected by flat mirror FM 1 , with its positions ( 1 ,  1 ) sensitive to yaw and pitch, respectively.These readings are not affected by the linear position of the mirror.The positions ( 2 ,  2 ) of the photodetector PD 2 , which receives the beam reflected by the retroreflector RR, are sensitive to the straightness errors in the horizontal and vertical plans.The beam position is reflected by 2 flat mirrors FM 2 and FM 3 (oriented at 45 ∘ ) and represented by the positions ( 3 ,  3 ) on the photodetector PD 3 .This constitutes the cumulative effects of the angular errors (pitch and roll) as well as the straightness in the horizontal plane.Based on the dimensional and geometrical configuration of the measurement system and interactions between data from the photodetectors, some arithmetic relations were established to calculate the various errors.Thus, for each relative movement between the two plates (mobile and references plates), six coordinates are available to evaluate the five components of error.Redundant measurement can be used as checking criterion.In the case of the machine tool, measurements obtained in this manner reflect the real position of the tool related to the part and make it possible to evaluate the deviations that it undergoes compared to its programmed trajectory.These deviations can be caused by geometrical, static, thermal, and dynamic sources.Once the   and   positions are calculated for PD 2 , the deviation in the horizontal plan is regarded as a horizontal straightness error while any deviation in the vertical plan is considered a vertical straightness error The two angular components   and   depend on the rotation around the two horizontal and vertical axes.Any rotation of the target with respect to the x-axis and the zaxis will be translated into yaw and pitch using trigonometric relationships ) . ( The roll error   depends on the rotation about the y-and z-axes.The angle of rotation about the y-axis, which is considered the axis of movement, is given by the trigonometric relationship

Experimental Setup and Preliminary Tests
With the aim of analyzing the responses of the new measurement system in order to determine its metrological characteristics, a coordinate measurement machine (CMM) is used to carry out evaluation tests.The diagrammatic representation of the experimental setup is illustrated in Figure 3.The overall system includes both the mobile and reference plates supporting the optical accessories.The fixed part of the measurement device is placed on the machine marble, whereas the moving part is fixed on the CMM head.Preliminary tests were carried out in order to optimize the utilisation conditions of the measurement system.In fact, it seems important to maximize the precision and the stability before starting the measurement validation.The obtained results show clearly that environmental conditions (light, vibrations, draught, and moisture) have a negative influence on the quality of the system responses.To evaluate the impact of these phenomena, an analysis of sensitivity was performed.This study made it possible to determine their actions and the mechanisms of their interactions and to release the conditions allowing the reduction of their effects.Thus, the two parts of the measurement system were covered so as to reduce the effects of the external sources of light.The use of the covers made it possible to improve substantially the stability and the luminosity of the beams.Figure 4 illustrates the improvements obtained and shows, that with this arrangement, the variation of the signals decreases from 1.50 m to 0.15 m relative to the horizontal plan and from 1.75 m to 0.25 m relative to the vertical plan.The use of filters also made it possible to reduce the noise and the adverse effects due to the variations of the other surrounding conditions.Moreover, with an aim of improving the repeatability of the system and of taking account of the effects heating, a 100-minute preheating period was observed before each test series (Figure 5).
To calibrate the system and to reduce the effect of environmental conditions, artificial neural networks present an interesting solution.In fact, models based artificial neural networks (ANN) are commonly used to compensate for the effects of temperature and relative humidity.The integration of these models in the compensation is easy and simple to implement.However, the model cannot be generalized for all surfaces of the position detector, because of nonlinearity due to the geometrical shape of the position detector.To overcome these limitations, it is necessary to calibrate the position sensor by developing an algorithm to estimate the horizontal and vertical positions while correcting for nonlinearity due to the concave shape of the surface position sensors.This algorithm must take into account the position of the beam on the projection surface from the position detector by combining all of the environmental factors.To get a rough estimation of the resolution readings, calibration tests were performed.During these tests, the mobile plate is mounted on the displacement member of the CMM related to the fixed plate so as to bring the laser beam to the center of each photodetector position.This position is then the reference position for position sensor calibration.A scan of the entire surface of the position detector is made to generate the necessary data for modeling.The pattern of data acquisition is shown in Figure 6.
The positions   and   of the laser beam on the photodetector are estimated in accordance with the voltages   ,   ,   , and   .The model consists of a representation of . . .

Output layers
Hidden layers the positions behavior based on voltages and the influence of ambient temperature (  ) and relative humidity (  ).The modelling process is performed by the ANN technique multilayer perceptron (MLP) known by its good forecast ability, simplicity, and flexibility of prediction.The MLP consists of three layers; the first and the third layers have the same number of neurons compared to the inputs and the outputs, respectively.The hidden layer is developed with taking into consideration the input and the output neurons and the complexity of the estimated parameters to avoid the overfitting.Several training attempts are important to determinate the 10 neurons of the hidden layer that met the best accuracy.An artificial neural controller capable of making appropriate corrections is developed to compensate for the effects of environmental conditions.Note that, to validate the effectiveness of the model (Figure 7), it is necessary to verify the model using new data that were not used during the training process.Tests on the coordinate measuring machine have led to the generation of useful data for training and for validation.
Figures 8 and 9 show the residual error of   and   positions during training and validations stages.During the training phase, the mean square errors (MSE) are, respectively, of the order of 76.15 m and about 92.81 m for the   and   positions.The variance is in the range of 0.601 m 2 for   position and about 0.650 m 2 for the   position.During the validation phase, the two mean square errors (MSE) are, respectively, of the order of 277.71 m and 240.27 m for the   and   positions.These errors are 3 times higher compared to training errors.The variance is of the order of 1.797 m 2 for the   position and about 1.795 m 2 for the   position.The residual errors have a relatively large variation with respect to that of the training phase.The results show that the data taken in the peripheral region of the photodetector detector have a time uncertainty larger than that measured, which corresponds exactly to the specifications given by the manufacturer of the optical element.To reduce the uncertainty, it is important to operate as close as possible to the central area by improving the total alignment of the measurement system (Table 1).
The analysis model was used to evaluate the influence of temperature on the position.On average, a temperature variation of about ±0.5 ∘ C causes a change in the beams  position of 1 m in the horizontal direction and 1.5 m in the vertical direction.However, the variation due to temperature is not linear, because it depends on the laser beam position on the photodetector surface and the uncertainty associated with the calibration of temperature sensor integrated into the experimental setup.Table 2 shows the influence of ambient temperature on the laser beam position.Also, it is possible to evaluate the effect of the relative humidity (RH) on the measurements.Consequently, each average RH variation of 0.5% HR position causes an average variation of 0.15 m in the horizontal direction and 0.20 m in the vertical direction of the beam position on the photodetector (Table 3).Following the impact assessment of the environmental conditions on the laser beam, the sensitivity analysis consists of quantifying the sensitivity of the signal under these disturbing factors.The resulting models help to estimate the influence of different sources of error.A neural compensator was necessary to reduce variations due to environmental conditions, taking into account the whole area of the detector.The stability has been improved due to the quick action of the compensator.The effort has controlled variations resulting from various sources, both systematic and random in nature.It remains to achieve metrological characterization of the new measurement system in order to determine its accuracy and uncertainty related to the operation and validate its performance using a multiaxial machine tool.

Experimental Results
The last phase of this project is to evaluate the metrological characteristics of the new measurement system (NMS) and validate their performances using a multiaxial machine tool by comparison with a precision interferometer laser (Renishaw-ML10).A protocol for the referencing and installation of the measuring system on the machine is built up to facilitate the alignment procedures and reduce setup time.
To get a good estimation of the 5 error components along the main displacement axis, a laser interferometer (ML10) with an accuracy of about 0.1 m and thermal stability on the order of ±0.0011 ppm was used to evaluate the measurements obtained by the NMS.Since the ML10 does not allow the measurement of more than one component at the same time, the adopted measurement procedure is to compare the NMS to the ML10 for each error component.In each test, a series displacement (front-back) was conducted to validate the linear (horizontal and vertical straightness) and angular deviations (yaw and pitch).With regard to the roll, there is no direct laser ML10 configuration to validate it.The calibration of the NMS generates algorithms to provide a reliable and effective measurement by reducing errors due to the total alignment of the NMS and systematic errors due to the measurement system itself.This calibration must be repeated after every system installation on a new machine.
A reset reading is required to define an absolute NMS to return to its initial position.The developed models are used to estimate the errors according to the beam positions on the photodetectors and the distance traveled by the mobile plate.
The validation tests are conducted on the CMM.The reference plate of the NMS is placed on the machine frame, while the mobile plate is fixed on the machine displacement head.The NMS is aligned according to an alignment protocol.Since the interferometer laser uses an elementary measurement approach, it cannot, however, validate more than one error component at each test.It is necessary to mount in parallel with the NMS while respecting the proposed configurations for each measure.
Figure 10 shows the result of the horizontal straightness measuring using both approaches.The NMS calibration is performed using a linear model which transforms data related to the horizontal position from the photodetector for horizontal position straightness.This calibration is to eliminate systematic errors and cosine errors that can affect the entire measurement system.Both curves (front and back) obtained by the ML10 laser have the same tendency, and they exhibit an offset in the form of mechanical hysteresis behavior, primarily reflecting geometrical errors and the guiding bridge CMM relative to the y-axis.The maximum straightness error is about 3.8 m.The two curves (return) obtained by NMS have the same tendency as reference curves.The total straightness error is on the order of 2.6 m.Regarding the back curve, the encoder receives an offset on the order of 0.25 m which is due to the mechanical hysteresis and the hysteresis CMM measuring system itself.
Figure 11 shows the result of vertical straightness measuring using both approaches.The NMS is calibrated using a linear model which transforms data of deviations in the vertical position.In this case, both curves (front and back) obtained using the ML10 laser have the same tendency and they exhibit an offset in the mechanical hysteresis.The maximum vertical straightness error is about 9.64 m.The curves obtained by NMS have the same tendency as the reference curves.The maximum vertical straightness error is on the order of 9.01 m.Regarding the back curve, the encoder receives an offset on the order of 0.15 m which is due to the mechanical hysteresis of the CMM and the hysteresis of the system measurement itself.The statistical analysis performed was used to compare the correlation between the two curves (front and back).A comparative study of the performance curves according to the displacement cycle is conducted.Regarding the sequence "front" the residual error between the curves obtained with the two techniques is on the order of 11.130 m for the horizontal straightness about 13.665 m for the vertical straightness.As for the sequence "back" the residual error is two times less than the cycle "front" for the horizontal straightness and about 9.371 m for vertical straightness.Analysis of the variance of the difference between the two curves (front and back) reveals a variation nearly equal to the vertical and horizontal straightness.This analysis was used to predict a rough estimation of the overall accuracy of the new measurement system based on the variance.The measurements obtained by the NMS are introduced through a calibration model to overcome the cosine and systematic errors as well as the correction of the basic equation used to calculate the yaw and pitch errors.This calibration takes into account the actual position in both horizontal and vertical directions as well as the distance traveled by the mobile plate.
Figure 12 illustrates the extent of the yaw using both approaches in a specific front and back cycle.The error evolution is nearly linear depending on the y-axis.The maximum error is about 5.9 arcsec.The results obtained by NMS are consistent with those taken by the ML10.The two curves (front and back) have the same tendency and they vary in a linear style as a function of displacement.However, there is an offset of about 0.45 arcsec caused by the hysteresis  between the NMS and the CMM.The maximum yaw error measured by the NMS, at about 5.9 arcsec, is comparable to that obtained by the laser ML10. Figure 13 illustrates the pitch measurement by using the two approaches in a front and back cycle.The curves obtained have the same trend as for the ML10 laser measurement.Unlike for the straightness measurements, the hysteresis error is minimal.The maximum error is bigger than the error of lace and is of the order of 9.456 arcsec.The results obtained by NMS are consistent with the measures taken by the ML10.However, there is offset of about 0.25 arcsec caused by the hysteresis between the NMS and the CMM.The maximum error is comparable to that obtained by the ML10 which is of the order of 9.474 arcsec.
It is also interesting to compare the two approaches by calculating the mean square error (MSE) and variance.Concerning the yaw, the variance is more pronounced between

Figure 1 :
Figure 1: Disposition of measurement system components.
and V D Positions P H and P V Deviations  ij and  ij

Figure 3 :
Figure 3: Schematic representation of the experimental setup.

Figure 7 :
Figure 7: ANN model architecture of positions prediction.

Figure 8 :
Figure 8: Residual error of training stage.

Figure 9 :
Figure 9: Residual error of validation stage.

Table 1 :
Comparison of training and validation performances.

Table 2 :
Influence of temperature on the positions.

Table 3 :
Effect of relative humidity (RH) on the positions.