Issue 
Int. J. Metrol. Qual. Eng.
Volume 14, 2023



Article Number  4  
Number of page(s)  7  
DOI  https://doi.org/10.1051/ijmqe/2023003  
Published online  24 May 2023 
Research article
Estimation of parallelism measurement uncertainty according to the Geometrical Product Specifications standard using coordinate measuring machine
^{1}
Laboratoire de recherche PCMT, Mohammed V University of Rabat, Avenue de l’Armée Royale, 10100 Rabat, Maroc
^{2}
Measurement Control Center MCC, 269, Zone Industrielle SUDOUEST, 28810 Mohammedia, Maroc
^{*} Corresponding author: nabil_habibi@um5.ac.ma
Received:
8
November
2022
Accepted:
29
March
2023
The dimensional and geometrical specifications control of a mechanical part is commonly done using a coordinate measuring machine. Collecting and processing measurement data through a probe, it allows measuring the conformity of the part according to specific tolerances. An accurate estimation of the uncertainty of measurement is critical; it is within this framework that our study is oriented; we aim to estimate the measurement uncertainty of parallelism error using the Guide to the Expression of Measurement Uncertainty (GUM), then we will proceed to a Monte Carlo simulation to compare the obtained uncertainty, we will then proceed to an interlaboratory comparison to validate our model. Our contribution is based on a more detailed and precise estimation of the uncertainties of the measurement process taking into accounts the calibration of the machine and the propagation of uncertainties.
Key words: Coordinate measuring machine (CMM) / parallelism error / guide to the expression of uncertainty in measurement (GUM) / Monte Carlo simulation (MCS) / interlaboratory comparison (ILC) / parallelism uncertainty
© N. Habibi et al., Published by EDP Sciences, 2023
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Considering the increasingly demanding quality that the modern industry is experiencing, accurate and precise inspection tools are required to go along the technological progress of manufactured parts, thus making the coordinate measuring machine (CMM) one of the widely used solution to control dimensional and geometrical specifications. Collecting inspected points through a probe, the software then processes data and creates fitted surfaces using the least square method, then measures the conformity of the part according to a specific criterion, which is a parallelism error in our case. This succession of steps is subject to a propagation of uncertainties, which if not estimated right, can lead to a high risk on the declaration of conformity. An accurate estimation of the uncertainty of measurement is critical and has a direct impact on decisionmaking; it is within this framework that our study is oriented.
Several studies have been made in order to estimate the form error measurement uncertainty in a CMM. Rosenda et al. [1] proposed a simplified model to estimate the measurement uncertainty deviations in CMM when determining circularity and cylindricity based on identifying the variables that influence the determination of geometry deviation. Ilyas Khan et al. [2] coded a genetic algorithm based on the minimum zone criterion. The algorithm estimates the parallelism error considering both datum and measured planes flatness default. Zha et al. [3] propose a strategy to estimate and minimize parallelism default using CMM. Wladyslaw and Wojciech developed in reference [4] a CMM measurement uncertainty evaluation software fully consistent with the GPS norm. Wojciech [5] then proposed different models of geometrical deviations associated uncertainties, including parallelism.
Referring to the GUM approach and the Monte Carlo method is usually done to estimate the measurement uncertainty. Egidi et al. [6] developed a model using the Monte Carlo method to estimate the uncertainty during the measurement of roundness on a coordinate measuring machine. Balasubramanian et al. [7] estimated uncertainty in angle measurement using the GUM method taking into consideration the geometrical errors, temperature, vibrations, and probing strategy. Using a comparison between these two methods (GUM and the MCM) to validate an uncertainty measurement model, has been proven to give consistent results, it's within this framework that Jalid proposed a comparison of these two methods on the estimation of flatness uncertainty [8] based on the orthogonal distance regression (ODR) that provided the parameters of the substitute plane and their uncertainties, then studied the influence of sample size on the flatness estimation and uncertainty [9].
In this paper, we aim to introduce a parallelism error model and estimate the associated measurement uncertainty using the guide to the expression of measurement uncertainty [10] and compare it to the results obtained with the Monte Carlo method [11], and then we will proceed to an interlaboratory comparison to validate our model. Our model combines the experimental and the analytical methods to estimate the parallelism error and the associated uncertainty taking into consideration the calibration of the machine and the propagation of uncertainties through the measurement process and following the normative guidelines.
2 Materials and methods
In order to estimate the parallelism associated uncertainty, we applied the following approach:
Setting a parallelism error equation according to the GPS standard.
Applying the GUM uncertainty propagation model to the parallelism equation.
Coding an algorithm on Matlab to estimate the datum plane parameters and associated variancecovariance matrix.
Estimating the variancecovariance matrix of the measured points p_{max} and p_{min.}
Calculation the parallelism uncertainty using the GUM
Validating the GUM results using a Monte Carlo simulation.
Confirmation of our model through an interlaboratory comparison using normalized error.
2.1 Parallelism error mathematical model
Based on ISO 1101 [12], Parallelism is an orientation tolerance, defined in as the minimum distance between two theoretical parallel planes P1 and P2, both parallel to the datum plane P0, within which all measured points lie inside (Fig. 1).
Let p_{max} and p_{min} be the two most distant measured points such as:
p_{max} (x_{max}, y_{max}, z_{max}) ∈ P1 and p_{min} (x_{min}, y_{min}, z_{min}) ∈ P2
The parallelism can be expressed in the following way:
where (n_{xo}, n_{yo}, n_{zo}) represents the datum plane's normal vector.
Fig. 1 Parallelism error. 
Fig. 2 CMM used. 
Fig. 3 Distribution law of MCM. 
2.2 GUM uncertainty propagation model
The law of propagation of uncertainty GUM (Guide to the expression of uncertainty in measurement) is an analytic method based on a first order Taylor expansion of a function Y = f (X_{1}, X_{2}, X_{3}, … , X_{n}). The final uncertainty is obtained by propagating the elementary components through a linear approximation. To estimate the parallelism error uncertainty modeled in the equation (1), we relied on the variance propagation law given by the GUM guide applied to :
where J is the Jacobian matrix given bellow:
M is the variance covariance matrix of the variables given as follows:
where the terms of are the followings:
And the terms of [p_{i}] are the followings:
2.3 Datum plane parameters (, A_{0}) and associated variancecovariance matrix
Surface fitting is critical step; in our case, it is about associating an ideal plane from the measured points that represents the best the probed plane. To do so, several criteria are available and can be considered according to the nature of the need, the least squares method remains the most frequently used, it is a standard regression approach to approximate the solution of over determined systems:
The best fit in the leastsquares sense minimizes the sum of squared residuals
where A and are the substitute plane parameters respectively representing a point belonging to the substitute plane and its normal vector and P_{i} are the probed points. This leads us to the equation to minimize as follows:
The starting parameters are critical and can lead to local solution, hence our choice of:
the barycenter of probed points.
an initial normal vector based on the most distant points apart B, C and D.
To solve this equation, we have developed an algorithm on Matlab using the function “nlinfit” which allows us to have the parameters of the fitted plane and its associated variancecovariance matrix.
2.4 Estimation of the variancecovariance matrix of a measured point [p_{i}]
CMMs are reliable and accurate. However, several factors may be sources of measurement uncertainty including geometric errors, measuring probe errors, thermally induced errors, vibrations, etc. The measurement uncertainties are not constant over the entire measurement space and varies according to the position of the axis, these same errors are dependent.
Several estimations of the measurement uncertainty have been made, Bahassou et al. [13,14] proposed the following matrix based on the uncertainty propagation law of the CMM calibration model according to the standard ISO 10360 [15]:
The uncertainties defined in the matrix represent the standard uncertainty, which will be divided by the coverage factor k; we will take k = 2 for the rest of the calculations.
2.5 Validation of the GUM model using Monte Carlo method
The estimation of measurement uncertainties by the Monte Carlo method is based on the propagation of distributions of random input quantities through a mathematical model of the measurement model. This is a great alternative when the law of uncertainty propagation (GUM propagation of variances) presents some difficulties, such as when the probability density function for the output variable deviates significantly from a Gaussian distribution giving unrealistic confidence intervals, or when the propagation based on the Taylor expansion to the first order presents an inadequate linearization of the model.
The Monte Carlo approach can also be used as a validation method of the GUM propagation of variances results by following this procedure:
– Applying the GUM uncertainty propagation to calculate and such as:
where U (dp) represents the uncertainty associated to “dp” obtained by applying GUM
– Applying the Monte Carlo method and deduce and as they represent the limits for a 95.45% confidence interval (dp±2σ) of the generated distribution with known mean value and deviation.
– Setting the tolerance associated with the uncertainty:
– Calculating d_{low} and d_{high} given as follows:
– Comparing d_{low} and d_{high} to ζ, the validation criteria is based on comparing the absolute differences of the respective endpoints of the two coverage intervals:
ζ ≥ max(d_{low}, d_{high}) .
Then, if both d_{low}and d_{high} are no larger than ζ, the comparison is favorable, and the GUM uncertainty framework has been validated in this instance.
3 Results and discussion
The goal of this experimental study is to put into practice the theoretical model made above. To do so, we worked on a prismatic mechanical part, we measured the datum plane for fifteen repetitions in 10 points, to be able to evaluate experimentally the uncertainty matrix of the normal vector and compare it to the one obtained using the “nlinfit” function in Matlab, we then measured the substitute surface in 10 points. This method allows to evaluate the effect of the association criterion according to the method of least squares relative to the datum plane.
Laboratory: PCMT.
CMM used: Mitutoyo EuroC 544.
Software: Geopak.
Measuring volume: X = 500 mm, Y = 500 mm, Z = 400 mm.
Probe type: TP2.
Stylus diameter: 2 mm.
MPE: ET = ±(4 µm + L/200) with L in mm.
Temp: 20 ± 2°C.
The results are as follows:
The datum plane normal vector with its associated uncertainty matrix:
The tolerance plane two extreme points and their associated uncertainty matrix:
3.1 GUM application
The final uncertainty is obtained by propagating the elementary components through a linear approximation. Few simplifications can be done to calculate the Jacobian matrix:
The final form of the Jacobian matrix is given as follows:
The final form of the variance covariance matrix is given as follows:
which leads us to the parallelism error and its associated uncertainty:
3.2 Monte Carlo simulation
To begin the simulation, we generated the parallelism equation on the PFSoftAlgebraic Calculator software.
We then loaded the equation on the Monte Carlo editor software, based on the hypothesis that each of the parameters follow a normal distribution law with standard deviation σ=U/k from their respective matrices with a coverage factor k = 2.
After generating a sample of size 10^{6} we have obtained the following distribution function:
We obtain the following values:
Which leads us to the expanded parallelism associated uncertainty U = 2σ = 0.00459566 mm.
The numerical tolerance is ζ = 0.5E − 3mm, and (d_{low}, d_{high}) calculation is the following as they represent the limits for a 95.45% confidence interval () of the generated distribution:
After comparison, the numerical tolerance is significantly higher than d_{low} and d_{high}, meaning that the validation criterion max (d_{low}, d_{high}) ≤ ζ is verified, thus the validation of our uncertainty estimation.
3.3 Interlaboratory comparison
In order to determine whether our model is valid, an interlaboratory comparison must be conducted. ILCs aim to determine the efficacy and suitability of a certain analytical technique, in our case, can be used to validate our model's results.
Laboratory: Measurement Control Center (MCC).
CMM used: Zeis Duramax.
Software: Calypso.
Measuring volume: X = 500 mm, Y = 500 mm, Z = 500 mm.
Probe type: Vastxxttl3.
Stylus diameter: 2 mm.
Temp: 20 ± 2 °C.
p_{L} = 5.3 μm.
U_{L} = 3.3 μm.
While many methods of evaluation exist, the commonly used comparison method is the normalized error (E_{n}) formula, calculated using this equation [16,17]:
where p_{L} is the parallelism error measured by the participant laboratory and U_{L} its associated measurement uncertainty. if  E_{n} ≤ 1 the result is satisfactory.
After calculation, the normalized error is significantly inferior to one, we can then conclude that our model based on the GUM is valid, and that the CMM used in the PCMT laboratory is accurate.
4 Conclusion
The present paper shows a different approach for estimating the uncertainty of parallelism error of an industrial part using CMM. We began by setting the mathematical model of the parallelism of the surface according to the ISO 1101 specifications. We proceeded next to the inspecting the tolerance surfaces by probing a set of data points and estimating their measurement uncertainties according to the ISO 10360 instead of using the MPE given by the constructor. We then proceeded to the creation of the fitted least square datum plane and its associated parameters uncertainties using the algorithm we coded on Matlab, which allowed us to estimate and compare the associated uncertainties using the GUM method and Monte Carlo method, both showing remarkably close results (gap less than 10^{−4} mm). We finally proceeded to an interlaboratory comparison; The normalized error showed satisfactory results, which validates our developed model.
References
 R.V. Arencibia, C.C. Souza, H.L. Costa, A. PiratelliFilho, Simplified model to estimate uncertainty in CMM, J Braz. Soc. Mech. Sci. Eng. 37, 411–521 (2015) [Google Scholar]
 M.I. Khan, S. Ma, S. Saifullah, New method of plane to plane parallelism error measurement based on flatness plane, Adv. Mater. Res. 941–944, 2209–2218 (2014) [Google Scholar]
 J. Zha, H. Zhang, Y. Chen, A strategy to evaluate and minimize parallelism errors of a rotor system in a precision rotary table, Int. J. Adv. Manufactur. Technol. 106, 3641–3648 (2020) [Google Scholar]
 W. Jakubiec, W. Plowucha, First coordinate measurement uncertainty evaluation software fully consistent with the GPS philosophy, 12th CIRP Conference on Computer Aided Tolerancing, Proc. CIRP 10, 317–322 (2013) [Google Scholar]
 W. Plowucha, Uncertainty of coordinate measurement of geometrical deviations, 15th CIRP Conference on Computer Aided Tolerancing, Proc. CIRP 75, 361–366 (2018) [Google Scholar]
 M. Balasubramanian, R. Vignesh, R. Rajamani, Evaluation of uncertainty in angle measurement performed on a coordinate measuring machine, ICCAP 2021: Proceedings of the First International Conference on Combinatorial and Optimization (2021) [Google Scholar]
 A. Egidi, A. Balsamo, D. Corona, Monte Carlo method for evaluating the uncertainty of roundness measurements on a CMM − a Python implementation, INRIM Istituto Nazionale di Ricerca Metrologica, Report number: R.T. 14/2021 (2021) [Google Scholar]
 A. Jalid, S. Hariri, A. El Gharad, J.P. Senelaer, Comparison of the GUM and Monte Carlo methods on the flatness uncertainty estimation in coordinate measuring machine, Int. J. Metrol. Qual. Eng. 7, 302 (2016) [Google Scholar]
 A. Jalid, S. Hariri, N.E. Laghzale, Influence of sample size on flatness estimation and uncertainty in threedimensional measurement, Int. J. Metrol. Qual. Eng. 6, 102 (2015) [Google Scholar]
 JCGM 100:2008, Evaluation of measurement data − Guide to the expression of uncertainty in measurement (2008) [Google Scholar]
 JCGM 101, Evaluation of measurement data − Supplement 1 to the Guide to the expression of uncertainty in measurement − Propagation of distributions using a Monte Carlo method, BIPM Joint Committee for Guides in Metrology, Sevres (2008) [Google Scholar]
 International Organization for Standardization ISO 1101: 201702, Geometrical product specifications (GPS) − Geometrical tolerancing − Tolerances of form, orientation, location and runout, Norm (2017) [Google Scholar]
 K. Bahassou, A. Salih, A. Jalid, M. Oubrek, Modeling of the correction matrix for the calibration of measuring machines, IJMET 8, 862–870 (2017) [Google Scholar]
 K. Bahassou, A. Salih, A. Jalid, M. Oubrek, Measurement uncertainty on the correction matrix of the coordinate measuring machine, Int. J. Adv. Res. Eng. Technol. 10, 669–676 (2019) [Google Scholar]
 International Organization for Standardization ISO 103602: CMMs Used For Measuring Linear Dimensions, Norm (2009) [Google Scholar]
 A. Softić, N. ZaimovićUzunović, H. Basić, Proficiency testing and interlaboratory comparisons in laboratory for dimensional measurement, J. Trends Dev. Mach. Assoc. Technol. 16, 115–118 (2012) [Google Scholar]
 International Organization for Standardization ISO 13528:2015, Statistical Methods For Use In Proficiency Testing By Interlaboratory Comparison, Norm (2015) [Google Scholar]
Cite this article as: Nabil Habibi, Abdelilah Jalid, Abdelouahab Salih, Hassane Hanane, Estimation of parallelism measurement uncertainty according to the Geometrical Product Specifications standard using coordinate measuring machine, Int. J. Metrol. Qual. Eng. 14, 4 (2023)
All Figures
Fig. 1 Parallelism error. 

In the text 
Fig. 2 CMM used. 

In the text 
Fig. 3 Distribution law of MCM. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.