Measurement uncertainty in process of line scales calibrating.
Runje, Biserka ; Medic, Srdjan
Abstract: The paper presents main characteristics of the device for
calibration of line scales and measurement uncertainty evaluation by GUM
and MCS method. As a part of research on the impact of measurement
uncertainty the following was investigated: the position of laser light
sources and optical components, minimizing Abbe's error (Bosse at
al. 2007), the determination of the middle line of line scales,
alignment of line scale and laser beam, straightness movement of table,
pitch, roll and yaw angles, environmental conditions affect the laser
wavelength and the geometry of device and the impact of losing focus
while moving of table. Measurement uncertainty evaluation has been
validated in comparison measuremens EURAMET Key Comparison, EURAMET.L-K7
"Calibration of line scales"
Key words: measurement uncertainty, line scale, length
1. INTRODUCTION
The Laboratory for Precise Measurement of Length, which is at the
same time the National Laboratory for Length (in text
'Laboratory') takes part in CIPM MRA comparisons of length
standards, which include line scales as very important standards of
length. Calibration of the line scales at the level of measurement
uncertainties of the order of value U = 0,1 [micro]m, k = 2. P = 95%
represents today still a world problem, although these levels of
measurement uncertainties are necessary in the context of ensuring the
traceability. So, the Laboratory started to design their own
optoelectronic system for the calibration of line scales.
2. MEASUREMENT DEVICE FOR CALIBRATING OF LINE SCALES
The measuring range of the device is 800 mm and it is primarily
intended for the calibration of line scales. The sighting process is
done by means of a microscope with a digital CCD camera Olympus DP 70
with 12, 5 Megapixels.
The microscope is fitted with lens of different magnification (10X,
20X, 50X). The lenses are selected in compliance with the object of
measurement.
The measuring system used is the laser interferometer (Reinshaw ML
10). The basis of the Renishaw Laser Interferometer system is He-Ne
Laser operating at a wavelength of 0,663 [micro]m. Measurement device
for calibrating of line scales is presented in Figure 1. In order to
achieve order in the above-mentioned measurement uncertainties, it is
necessary to use software in the process of detecting the line centre of
the measuring scale in reference to requirement limits (Beers and
Penzes,1999) The software solution functions in such a way that all the
pixels of a certain image are transmitted into a black & white
combination and then the position of the line centre is calculated by
arithmetic algorithms (Druzovec at. al. 2009).
The software solution provides the exact position of the line
centre in pixels. In order to convert the values in pixels into the
length values, it is necessary to calibrate the pixels size, i.e. to
find out the length value of every pixel.
[FIGURE 1 OMITTED]
3. CALCULATION OF THE MEASUREMENT UNCERTAINTY BY APPLYING GUM AND
MCS METHOD
The mathematical model of measurement has been given by expression
(1):
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)
where:
[N.sub.i]--Number of wavelengths
[lambda]--Laser wavelength
[n.sub.air.sup.-]--Refractive index of air
[delta][l.sub.ni]--Interferometer nonlinearity
[delta][l.sub.DP]--Deadpath influence
[delta][l.sub.li]--Interferometer cosine error
[delta][l.sub.Az]--Abbe offset in z and pitch
[delta][l.sub.Ay]--Abbe offset in y and yaw
L--Nominal length of line scale
[[alpha].sub.s]--Thermal exp. Coeficient
[DELTA][t.sub.s]--Deviation scale temperature from 20[degrees]C
[delta][l.sub.sh]--Scale alignement horizontaly
[delta][l.sub.sv.sup.-]--Scale alignement verticaly
[delta][l.sub.ai]--Scale support influence
[delta][E.sub.alg]--Line quality influence
[delta][e.sub.fok]--Focus loosing influence
[delta][l.sub.opt]--Uncertainty of measurement optics due to temp.
dev.
[delta][l.sub.sE]--Reproducibility of line detection
The yields of components of the standard uncertainty for the line
scale of 100 mm are presented in Table 1.
Calculation of the measurement uncertainty (validation) has also
been performed, by means of MCS method (JCGM 101:2008.) Probability
density function of the output value has been obtained by M = 100000
simulations. The probability density function g([x.sub.i]) has been
simulated by the MCS method based on the expression (1). Figures 2 and 3
show the probability density functions of the output value [L.sub.MS]
where the distance between spots of reference and reflected beams are s
= 2 mm and s = 5 mm respectively.
While the GUM method assumes normal distribution of the output
value, the MCS method yielded experimental distribution of the output
value that may more or less match the assumed normal distribution. The
form of the experimental curve will depend primarily on the probability
density function of the most significant input value (Medic et al.,
2003).
[FIGURE 2 OMITTED]
[FIGURE 3 OMITTED]
In this case, due to the dominant influence of interferometer
cosine error (Quenelle, 1983.) on the measurement uncertainty, the
normal distribution assumes, through length increase, the
characteristics of a trapezoid distribution (Fig. 3).
4. CONCLUSION
By designing the measurement system for calibration of precise line
scales, the Laboratory has opened the possibility of carrying out the
international comparisons in the field of line scales. Thus, the
Laboratory participated in the EUROMET project 882 "Calibration of
line scales", L-K7. Intercomparison results of measuring the length
of the 100 mm line scale are presented in Figure 4. Figure 4 shows that
obtained results of Laboratory have no significant deviation compared to
average of results of METAS, PTB and MIKES and that they have the same
trend.
[FIGURE 4 OMITTED]
The participation in this international comparison measurement was
representing a real validation of the device and evaluated measurement
uncertainty. The obtained results of this comparison will be good
indication about direction of future research in a way to reduce
measurement uncertainty in calibration of line scales.
6. REFERENCES
Bosse, H.; Flugge, J.; Koning, R. A. (2007), Method for the in situ determination of Abbe errors and their correction, Measurement science
and tehnology, 18, 476-481
Druzovec, M.; Acko B.; Godina, A.; Welzer T. (2009), Robust
algorithm for determining line centre within a video positional
measuring system, Optics and Laser sin Engineering 47
Beers, J. S; Penzes W. B. (1999), The NIST Length Scale
Interferometer, Journal of Research of the National Institute of
Standards and Technology
Medic, S.; Mudronja, V.; Runje, B. (2003), Examples of Applying
Monte Carlo Simulations In The Field of Measurement Uncertainties of The
Standard of Length // Proceedings of the XVII IMEKO World Congress.
Cavtat, Croatia
*** JCGM 101:2008, Evaluation of measurement data--Supplement 1 to
the "Guide to the expression of uncertainty in
measurement"--Propagation of distributions using a Monte Carlo
method
Tab. 1. Yields of components of standard uncertainty,
and sources of uncertainty
Amount
of
Source an d Component Stand.
of uncertaint
Uncertainty, [x.sub.i] Distr. y u([x.sub.i]
Abbe offset in z and R 16,8 not
catch, [delta]IAz
Abbe offset in y and R 4,3 nm
yaw, [delta][I.sub.Ay]
Laser Wavelength, R 0,03
[delta][lambda]
Air temperature, tair R 0,12[degrees]C
Air pressure, pair R 13 Pa
Relative humidity, RHair R 0,06
Edlen equation N 2 x [10.sup.-8]
uncertainty, [delta]nair
Deadtath, [delta]IDP R 1,8 nm
Interferometer U 3 nm
nonlincari, [delta]INL
Interferometer cosine R 0,48L
error, [delta][I.sub.Ii]
Deviation scale N 0,12[degrees]C
temperature from 20
[degrees]C, [DELTA]ts
Thermal exp. Coef, R 0,289 x [10.sup.-7]
[[alpha].sub.s],
[K.sup.-1]
Scale alignement hor., R 0,001L
[delta]ISh
Scale alignement vert, R 0,0023L
[delta]ISV
Scale support, [delta]lai R 0,0058L
Line quality, [delta]Ealg N 6,4 nm
Focus loosing, [delta]fok N 18 nm
Measurement optics, R 58 nm
[delta]lost
Interferometer R 0,003
resolution, N
Reproducibility of line N 11,6 mn
detection, [delta]ISE
Combined variance
Linearised expanded
measurement uncertainty
U,P = 95%, k = 2
Yield to
Source an d Component [c.sub.i] = [partial measure.
of derivative]dL/[partial uncertainty,
Uncertainty, [x.sub.i] derivative][x.sub.i] nm, L in mm
Abbe offset in z and 1 16,8
catch, [delta]IAz
Abbe offset in y and 1 4,3
yaw, [delta][I.sub.Ay]
Laser Wavelength, L 0,03L
[delta][lambda]
Air temperature, tair 9,5 x [10.sup.-7] 0,112 x L
L/[degrees]C
Air pressure, pair 2,7 x [10.sup.-7] L/Pa 0,035 x L
Relative humidity, RHair 8,5 x [10.sup.-7] L 0,050 x L
Edlen equation L 0,020 x L
uncertainty, [delta]nair
Deadtath, [delta]IDP 1 1,8
Interferometer 1 3
nonlincari, [delta]INL
Interferometer cosine 1 0,48 x L
error, [delta][I.sub.Ii]
Deviation scale 5 x [10.sup.-7] L/K 0,06 x L
temperature from 20
[degrees]C, [DELTA]ts
Thermal exp. Coef, L x 0,5 K 0,0145 x L
[[alpha].sub.s],
[K.sup.-1]
Scale alignement hor., 1 0,001 x L
[delta]ISh
Scale alignement vert, 1 0,0023 x L
[delta]ISV
Scale support, [delta]lai 1 0,0058 x L
Line quality, [delta]Ealg l 6,4
Focus loosing, [delta]fok 1 18
Measurement optics, 1 58
[delta]lost
Interferometer [lambda]/2 1
resolution, N
Reproducibility of line 1 11,6
detection, [delta]ISE
Combined variance [u.sup.2] = ([65.sup.2] +
[0,5.sup.2] x [L.sup.2]) nm, L in mm
Linearised expanded U = (130 + 0,66 x L) nm, L in mm
measurement uncertainty
U,P = 95%, k = 2