On the uncertainty of measurements by measuring the form of a surface/Mootemaaramatuse uurimine pinna kontuuri mootmisel.1. INTRODUCTION
Uncertainty of measurement is, by its definition, a parameter, associated with the result of the measurement and characterizing the dispersion of the values that can be attributed to the measurand [1,2]. It reflects the lack of exact knowledge about the value of the measurand. Thus owing to the uncertainty, arising from random effects and from imperfect correction of the result for systematic effects, the result of the measurement after correction for recognized systematic effects is still only a rough estimate of the true value of the measurand. For this reason, each measuring result should be associated with information about the uncertainty, identifying the possible dispersion of the true value of the measurand. In metrology laboratories, mostly standardized procedures are used in evaluating measuring uncertainty. However, these procedures require extended statistical and mathematical knowledge, which usually is not available in industry.
In [3-6], surface roughness was measured with a roughness measuring instrument. The uncertainty of measurement results could be estimated only by the uncertainty contribution of the measuring instrument. This forms about 10-15% of the total indication. Besides that, no other contributors of uncertainty were considered to estimate the measurement results. In [7-9], step height was measured with a surface roughness measuring instrument. To evaluate the measurement results, in addition to the uncertainty, contributed by the measuring instrument, the uncertainty caused by the measurer was considered. As a result, the measurement results became more reliable. In above-mentioned papers, however, other contributions to uncertainty--the stylus radius, measurement force, surface angle, etc.--were ignored.
The current research has the aim to consider all the possible uncertainty contributors, essential in estimating a surface contour of a complicated form.
2. EVALUATION OF THE UNCERTAINTY OF MEASUREMENTS
Uncertainty in measurement comprises, in general, many components. Some of these components may be evaluated from the statistical distribution of the results of a series of measurements and can be characterized by experimental standard deviations (type A evaluation of uncertainty). The other components, which can also be characterized by standard deviations, are evaluated from assumed probability distributions, based on experience or other information (type B evaluation of uncertainty).
According to the reference document , the first step in determining the uncertainty of a measurement is to calculate the model function f that shows the relationship between the input quantities ([X.sub.1], [X.sub.2],..., [X.sub.N]) and the quantity to be measured Y:
Y = [florin] ([X.sub.1], [X.sub.2],..., [X.sub.i],..., [X.sub.N]) (1)
Model function [florin] represents the procedure of measurement and the method of evaluation. It describes how values of the output quantity Y are obtained from the values of the input quantities [X.sub.i] In most cases it is an analytical expression, but it may also be a group of such expressions which include corrections and correction factors for systematic effects, thereby leading to a more complicated relationship that cannot be written down as an explicit function. In this case, [florin] may be determined numerically.
An estimate of the measurand , Y the output estimate denoted by , y is obtained from Eq. (1) using input estimates [x.sub.i] for the values of the input quantities [X.sub.i]:
y = [florin]([x.sub.1], [x.sub.2],..., [x.sub.i],..., [x.sub.N]). (2)
First, standard uncertainties u ([x.sub.i]) of all input estimates [x.sub.i] should be evaluated. For an input estimate of [X.sub.i] obtained from the statistical analysis of a series of observations (type A evaluation of uncertainty), standard deviation of the mean value of [X.sub.i] is calculated as
S([bar.X.sub.i]) = S([X.sub.ij])/[square root of ([n.sub.i])](3)
Variance [s.sub.2], ([x.sub.i,j]) for non-correlated input values is calculated as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)
If some of the input quantities are correlated, the correlation should be considered in Eq. (4).
Standard uncertainty ( ) i u x is equal to the standard deviation of the mean
u(xi) = s([[bar.x].sub.i]) (5)
Input values are the best estimates that were corrected in terms of all effects, significant for the model. If that was not the case, the necessary corrections were introduced as separate input quantities.
Due to insufficient knowledge, estimations of the input quantities are not exact, leading to uncertainty, characterized by standard deviation of the output quantity Y. The calculation of the output quantity is performed applying the law of propagation of variances to Eq. (1):
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)
Therefore it is necessary to know standard deviations, called standard uncertainties, of each of the input quantities u([[x].sub.i]). Depending on how the standard uncertainty is estimated, the set of input quantities may be divided into two categories [1,2]:
Evaluation Type A
Standard uncertainty of input quantities can be evaluated in the course of statistical analysis of a series of observations.
Evaluation Type B
The standard uncertainty of input quantities is evaluated by means of tools different from the statistical analysis of a series of observations. In this case, the information can come from the following sources: calibration certificates, handbooks, producers' specifications and hypotheses on the density function of the input quantity.
Calculating standard uncertainty of the output quantity on the basis of the law of propagation of variances, the expanded uncertainty of measurement U can be obtained when multiplying the standard uncertainty by a coverage factor : k
U = ku(y). (7)
The value of k depends on the probability distribution of the output quantity y and on the level of confidence. The assigned expanded uncertainty corresponds to a coverage probability of approximately 95%. In most of the calibrations the output distribution can approach a normal distribution with k = 2.
3. CONTOUR MEASUREMENT METHODS
Complicated surface contours were measured using the surface texture measuring system "Perthometer Concept" produced by company MAHR . A simplified schema of the system is shown in Fig. 1 and a photo in Fig. 2.
The Institute of Mechatronics has owned this system for about one year. "Perthometer Concept" is a modular computer-controlled station for measuring and analysing roughness, contour and topography of surfaces. Its software runs under Windows. Operation is therefore quickly learned, easy to understand and compatible with other Windows' applications. PCV 200 contour drive unit with an exchangeable tracing arm was used in our research. The high-precision PCV 200 contour drive unit is a long-distance instrument for the assessment of radii, distances, angles and straightness deviations. The smooth traverse and the computer-assisted error correction guarantee reproducible measurements with utmost vertical and horizontal resolution in a measuring field of 200 x 50 mm. PCV 200 contour drive unit allows automatic lowering and lifting of the tracing arm with programmable speed and quick positioning. The measuring force can be adjusted from 2 to 120 mN. Rigid design and unique material provide a highly dynamic construction. Drive unit has programmable measuring routines including lowering, lifting and positioning of the tracing arm and selectable measuring speeds.
[FIGURE 1 OMITTED]
[FIGURE 2 OMITTED]
4. MEASUREMENT MODEL
Proceeding from Eq. (2), the measurement model can be expressed as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (8)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (10)
where [delta][x.sub.MI] is the correction from the measuring instrument, [delta][x.sub.r] is the correction from the stylus radius, [delta][x.sub.F] is the measurement force correction, [delta][x.sub.cv] is the surface curvature correction, [delta][x.sub.cc] is the surface concavity correction and [delta][x.sub.ang] is the correction from the surface angle.
Now, we can express the measurement model by the following equation:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)
5. COMBINED UNCERTAINTY OF A MEASUREMENT RESULT
The standard uncertainty to be ascribed to the estimate y of an output quantity , Y which is evaluated from the estimates of a number of input quantities, is named [1,2] combined standard uncertainty. By introducing this concept, it is possible to distinguish the uncertainty of the output quantity from the uncertainty of other quantities that occur in the measurement model. However, the uncertainty of an input quantity is, in its turn, often obtained from a relevant measurement model, which means that during the evaluation process it was itself determined with an uncertainty. Similarly, we can use the output from the measurement model as an input for a measurement task. The concept of combined standard uncertainty is therefore only of limited use. The symbol u(y) is used for the standard uncertainty to be ascribed to the estimate , y regardless of the way in which the uncertainty has been evaluated. The combined standard uncertainty is the positive square root of the combined variance, which is the weighted sum of the experimental variances and covariances of all input quantities considered in the measurement model. The experimental variances and covariances are obtained from the experimental standard deviations u([x.sub.i]), associated with the estimates [x.sub.i] of the input quantities [X.sub.i]. In our case the combined standard uncertainty is determined as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)
Standard uncertainties of input quantities from different sources were determined. The following results were obtained experimentally and their standard uncertainties were calculated applying type B method of uncertainty evaluation.
Indication (contour) x
Indication in this case is a contour we can see on the screen of the computer (Fig. 3). Standard uncertainty of the indication can be determined according to the printer resolution. The current printer resolution is [Delta]x = [+ or -]1 [micro]m:
u(x) = [Delta]x/[square root of 3] = 0.6 [micro]m.
[FIGURE 3 OMITTED]
Measuring instrument correction [delta][x.sub.MI]
This correction was not found in the calibration certificate, but it was mentioned, that the indication can change within the limits of [[Delta].sub.MI] = [+ or -]0.5 [micro]m:
u([delta][x.sub.MI]) = [[Delta].sub.MI]/[square root of 3] [congruent to] 0.3 [micro]m.
Stylus radius correction [delta][x.sub.r]
Research results indicated that stylus radius correction does not remarkably affect the contour measurements. Thus we have assumed [delta][x.sub.r] [congruent to] 0 and u([delta][x.sub.r]) [congruent to] 0.
Measuring force correction [delta][x.sub.F]
Measuring force correction and its standard uncertainty can be calculated as follows. From the Hertz formula the elastic deformation can be calculated. The worst situation, sphere-sphere, was assumed. The correction value is to be considered equal to zero and its standard uncertainty can be calculated according to the following equation:
u([delta][x.sub.F]) = [[Delta].sub.F]/[square root of 3] [congruent to] 0.3 [micro]m,
where [[Delta].sub.F] = [+ or -]0.6 [micro]m.
Surface complexity corrections [delta][x.sub.cv], [delta][x.sub.cc]
The corrections due to the surface curvature and concavity are assumed to be equal to zero: [delta][x.sub.cv] [congruent to] [delta][x.sub.cc] [congruent to] 0.
The standard uncertainty of these corrections can be calculated from the following equations:
u([delta][x.sub.cv]) = [[Delta].sub.cv]/[square root of 3] [congruent to] 0.9 [micro]m,
u([delta][x.sub.cc]) = [[Delta].sub.cc]/[square root of 3] [congruent to] 0.9 [micro]m,
where [[Delta].sub.cv] . and [[Delta].sub.cc] . have been found experimentally: [[Delta].sub.cc] = [[Delta].sub.cc] = [+ or -] 1.5 [micro]m.
Correction of the surface angle [delta][x.sub.ang]
Correction of the surface angle [delta][x.sub.ang] = 0 and its standard uncertainty can be calculated from the equation:
u([delta][x.sub.ang]) = [[Delta].sub.ang]/[square root of 3] [approximately equal to] 1.7 [micro]m,
where ang . was experimentally determined during the research applying the angle standards: [[Delta].sub.ang] = [+ or -] 2.9 [micro]m.
The above-mentioned quantities and their estimations are presented in Table 1. From Eq. (12) we have
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]
Hence, the expanded uncertainty can be given, according to Eq. (7), as
U = ku(y) = 2 * 2.3 = 4.6 [micro]m [congruent to] 5 [micro]m.
A measurement model and method for calculating the expanded uncertainty of the surface measurement has been elaborated. It is possible to give an estimation to the surface elements obtained in the printout after measuring the contour of a complicated surface. It has been shown how to estimate the variation range and to analyse the limits within which the numerical values of the surface contour can change. Finally, the quality of the measured values can be evaluated and the measurement uncertainty of the latter can be estimated.
This work was supported by Estonian Ministry of Education targeted financing grant No. 14250503 and Estonian Science Foundation grant No. 6172.
Received 2 November 2006, in revised form 29 January 2007
[1.] Guide to the Expression of Uncertainty in Measurement (first edition 1993, corrected and reprinted 1995). International Organization for Standardization, Geneva, 1995.
[2.] Laaneots, R. and Mathiesen, O. An Introduction to Metrology. TUT Press, Tallinn, 2006.
[3.] Torims, T. Research on Machining of Diesel Engine Crankshaft Journals Surfaces. Riga Technical University, Riga, 2005.
[4.] Rudzitis, J., Torims, T., Konrads, G. and Shperbergs, J. Three-dimensional roughness analyses for air compressor cylinders. In Proc. 5th International Conference of DAAAM Baltic Industrial Engineering. Tallinn, 2006, 163-165.
[5.] Rudzitis, J., Shiron, E., Skurba, M. and Torims, T. Classification of rough surfaces. In Scientific Proceedings of Riga Technical University. Series 6. Machine Science and Transport. Production Engineering, vol. 2. RTU, Riga, 2001, 53-57.
[6.] Rudzitis, J., Torims, T., Gernish, E. and Konrads, G. Rough surfaces contact examination. In Scientific Proceedings of Riga Technical University. Series 6. Machine Science and Transport. Production Engineering, vol. 14. RTU, Riga, 2004, 19-24.
[7.] Hoffmann, K.-P. and Laaneots, R. Pindepaksusetalonide kalibreerimine. EVS Teataja, Tallinn, 1999, No. 7/8, 13-17.
[8.] Hasche, K., Hermann, K., Mirande, W., Pohlenz, F. and Seemann, R. Development of scanning force microscopy with subnanometric capability. In Proc. Fourth Seminar on Quantitative Microscopy QM 2000. Semmering, 2000 (PTB-Bericht PTB-F-39), 219-225.
[9.] Hoffmann, K.-P., Ahbe, T., Herrmann, K., Hasche, K., Pohlenz, F. and Sun, J. Development and calibration of standards for the coating thickness in the range of micrometer and nanometer. Surface Coating Technol., 2003, 169-170, 732-734.
[10.] Perthometer Concept. Operating Instructions. MAHR GmbH, Gottingen. Germany, 2001.
Indrek Abiline, Rein Laaneots, Maidu Nanits and Jurgen Riim
Institute of Mechatronics, Tallinn University of Technology, Ehitajate tee 5, 19086 Tallinn, Estonia; email@example.com
Table 1. Estimates of the input quantities and their uncertainties Quantity Estimate Standard Dispersion [X.sub.i] [x.sub.i] uncertainty [u.sup.2] u([x.sub.i]), ([x.sub.i]) [micro]m x contour 0.6 0.36 [[delta]x.sub.MI] 0 0.3 0.09 [[delta]x.sub.r] 0 0 0 [[delta]x.sub.F] 0 0.3 0.09 [[delta]x.sub.cv] 0 0.9 0.81 [[delta]x.sub.cc] 0 0.9 0.81 [[delta]x.sub.ang] 0 1.7 2.89 [Sigma] 5.05