Next Article in Journal
Prompt-Based Graph Convolution Adversarial Meta-Learning for Few-Shot Text Classification
Previous Article in Journal
Real-Time Interactive Parallel Visualization of Large-Scale Flow-Field Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Measurement Uncertainty Analysis of the Stitching Linear-Scan Method for the Measurable Dimension of Small Cylinders

1
School of Mechanical and Electrical Engineering, Lanzhou University of Technology, Lanzhou 730050, China
2
Shenzhen International Graduate School, Tsinghua University, Shenzhen 518131, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(16), 9091; https://doi.org/10.3390/app13169091
Submission received: 7 July 2023 / Revised: 4 August 2023 / Accepted: 8 August 2023 / Published: 9 August 2023

Abstract

:
A stitching linear scan method is proposed for roundness and diameter measurement of small cylindrical workpieces instead of the conventional rotary scan method due to the crucial alignments of eccentricity and inclination. To verify the reliability of the proposed method, by which the coordinates of the cross-sectional circle of a small cylinder is divided into several equal parts to be obtained and reconstructed, the diameter and roundness measurement uncertainties of the small cylinders with a diameter 1.5 mm are evaluated to be 0.047 μm and 0.095 μm, respectively, which can meet the uncertainty target of 0.1 μm. To investigate the measurable dimension by the proposed method, measurement uncertainty analysis of the small cylinders with various dimensions has been conducted according to the previous evaluation, since all the procedures are the same except for the uncertainty of X coordinates, which changes with the measured dimension’s change. The results show that the small cylinders with a diameter range from 0.01 mm to 50 mm can be measured by the proposed method when the position error θZ is reduced to the corresponding value. There is no measuring limitation set by the proposed model theoretically in the case of θX = θZ = 0.1°, while the machine has a measuring limitation.

1. Introduction

Fine mechanics, such as RV reducers, play a significant role in industrial robotics and machine tool spindles. There are many cylindrical parts with various diameters employed in RV reducers. The performance and life of an RV reducer are affected by the quality of these cylindrical parts very much [1,2,3,4,5]. The parameters of roundness and diameter are always used to evaluate the quality of the cylindrical parts. The more precise the parameter is, the better the quality is. To enhance the quality of the cylindrical parts, precision measurement is necessary [6,7,8,9,10,11]. Conventionally, roundness measurement can be carried out by the rotary scan method with a roundness measuring machine [12,13,14,15]. In recent years, some alternative roundness measurement methods have been proposed, such as the multiple probe method and online method [16,17,18,19]. However, due to the crucial alignments (inclination and eccentricity) of the rotary scan method, precision roundness measurement become difficult when the diameter of the measured workpiece becomes smaller. To address this issue, the stitching linear scan method is proposed [20,21,22]. In this method, the coordinates of the cross-sectional circle of the cylindrical workpiece, which is attached to a round magnetic jig and mounted on a V-groove, are divided into several equal arcs to be scanned linearly by a profilometer. A set of arc coordinate data can be obtained after one linear scan on the surface of the cylindrical workpiece. The measured cylindrical workpiece is rotated by an equal angle to a new measuring position by rotating the magnetic jig and then scanned again. Another set of arc data can be obtained. We then repeat the procedure several times till the entire cross-sectional circle of the cylindrical workpiece is scanned completely, then a series of arcs can be obtained. The arc profiles, which can be characterized according to the obtained arcs, are used to stitch into a roundness profile. Through the stitching angle error compensation, overlap parts integration, and filter processing, an accurate, integral, smooth, continuous roundness profile can be obtained [20,21,22].
Although small cylindrical parts can be measured by the stitching linear scan method, the measurable range of this method has not been analyzed. In this paper, measurement uncertainty is used to verify the reliability of the stitching linear scan method for the roundness measurement of small cylindrical parts with various diameters. Due to ambiguous random error and system error, the concept of measurement uncertainty was proposed in 1993. The measured result can be reasonably and quantitively evaluated according to the GUM (Guide to the Expression of Uncertainty in Measurement), which is actually an error propagation evaluation [23,24,25]. Measurement uncertainty is necessary for multiple disciplines, such as air gauge back-pressure uncertainty estimation for advanced test rigs and precision positioning measurement uncertainty [26,27]. In this paper, the analysis of measurement uncertainty is presented and the measurable range of the cylindrical workpiece can be confirmed.

2. Principles and Experiment

As is known to all, the surface form of a measured workpiece can be modeled and reconstructed by coordinates. In this paper, the coordinates of the cross-sectional circle of a small cylindrical workpiece are divided into several equal arcs to be scanned by a profilometer. The diameter and center coordinates of these arcs can be fitted by the least square method. The arc profiles can be characterized according to the arc data. The roundness profile of the small cylindrical part can be obtained by stitching these arc profiles one by one in sequence.
As shown in Figure 1a, a small cylindrical workpiece with a diameter of 1.5 mm and length of 7.8 mm is attached to the circular magnetic jig, which is marked by eight equal lines on the surface. This kind of combination is mounted on the V-groove. As shown in Figure 1b, the stylus of the profilometer is brought to scan on the surface of the small cylindrical workpiece and returned to the initial position, after which the first arc coordinate data can be obtained. The small cylindrical workpiece is rotated by 45° by rotating the circular magnetic jig by one equal part. The stylus scans on the surface of the workpiece again and returns, then the second set of arc coordinate data can be obtained. We repeat the procedure seven times, then eight sets of arc coordinate data can be obtained. Since the maximum measuring inclination of the stylus is ±45°, only 85° arc coordinate data are extracted for the following data processing.
The eight arc profiles with respect to the arc coordinate can be characterized by a series of data processes. An inaccurate roundness profile can be formed by stitching these eight arc profiles together. Regarding the stitching procedure, as shown in Figure 2a,b, arc profiles 1 and 2 are at the initial position. As shown in Figure 2c, the first arc profile is kept static and the second arc profile is rotated by 45°, since there is a 45° angular displacement between adjacent arc profiles according to the measuring procedure. As shown in Figure 2d, the second arc profile is kept static, and the third arc profile is rotated by 45°, while the rest of the arc profiles are stitched in this manner. However, there are some stitching angle errors between adjacent arc profiles, since the small cylindrical workpiece is rotated manually, namely, the angular displacement between adjacent arc profiles is not always 45°. Therefore, stitching angle error compensation is necessary. An accurate, integral, smooth, continuous roundness profile can be obtained after stitching error compensation, overlap part combination, and low-pass filtering.

3. Measurement Uncertainty Analysis of the Linear Scan Method

3.1. Mathematical Modeling

The uncertainty of roundness and diameter measurements of a small cylinder with a diameter of 1.5 mm and length of 7.8 mm by the stitching linear scan method is carried out at first. Mathematical modelling was conducted to calculate the measurement uncertainty of the diameter and roundness. The arc coordinate (xi, zi) in the rectangular coordinate system can be obtained by a profilometer with its stylus. The mathematical model used for radius calculation and the uncertainty factors deriving from the measurements are shown in Figure 3a. The radius of each arc can be fitted by the least square method. Radius Ri of an arbitrary measuring point can be expressed by Equation (1), where (xi, zi) is the arc coordinate.
R i = x i 2 + z i 2
As derived from Equation (1), the combined standard uncertainty u(Ri) of the measured arc radius can be expressed as follows:
u 2 R i = R i x i u 2 x i + R i z i u 2 z i = x i 2 z i 2 + x i 2 u 2 x i + z i 2 z i 2 + x i 2 u 2 ( z i )
The uncertainties in the figure are as follows.
(1)
Uncertainty of output zi in the Z-axis direction.
u(ecalibration_Z): Uncertainty of stylus calibration in the Z-axis direction.
u(eresolution): Uncertainty due to stylus resolution.
u(erepeat): Uncertainty of repeatability.
(2)
Uncertainty of output xi in the X-axis direction
u(ecalibration_X): Uncertainty of stylus calibration in the X-axis direction.
u(ealignment_Z): Uncertainty due to position error around the Z-axis of the workpiece.
u(ealignment_X): Uncertainty due to position error around the X-axis of the workpiece.
In summary, the combined standard uncertainty zi of the Z-axis output and the combined standard uncertainty of the X-axis coordinate xi are expressed by Equations (3) and (4), respectively.
u 2 z i = u 2 e c a l i b r a t i o n _ z + u 2 e r e s o l u t i o n + u 2 e r e p e a t
u 2 x i = u 2 e c a l i b a t i o n _ X + u 2 e a l i g n m e n t _ Z + u 2 e a l i g n m e n t _ X
For the calculation of roundness, it is assumed that there is a deviation from the approximate center of the circle, rather than the radian of a perfect circle, and the approximate center of the circle is set as the origin. The output of zi can be obtained by Equation (5), where Ri is the radius of the arbitrary point of the measured arc and R0 is the peak radius.
z i = R 0 R i 2 x i 2
The radius Ri after the stitching process is calculated according to Equation (6).
u z q = u r θ m a x + u r θ m i n
Δri is the radial deviation at any measured point and Δr0 is the radial deviation at the arc apex.
R i = R ¯ + r i ( R ¯ r i )
R 0 = R ¯ + r 0 ( R ¯ r 0 )
We square both sides of Equation (5), substitute this into Equations (7) and (8), and rearrange.
r i = 2 R ¯ z i z i 2 x i 2 2 R ¯
From Equation (9), the combined standard uncertainty uri) of the radial deviation of any measurement point can be expressed as:
u 2 r i = r i R ¯ 2 u 2 R ¯ + r i x i 2 u 2 x i + r i z i 2 u 2 z i = x i 2 + z i 2 2 R ¯ 2 2 u 2 R ¯ + x i R ¯ 2 u 2 ( x i ) + 1 z i R 2 u 2 ( z i )
The uncertainty factors shown in Figure 3b are similar to the radius measurement model stated above. Therefore, the combined standard uncertainty zi for the Z-axis output and the combined standard uncertainty for the X-axis coordinate xi are expressed by Equations (6) and (7), respectively. The uncertainty of the radius after the stitching process is expressed by Equation (11), where u(Ri) is the combined standard uncertainty of the radius obtained from Equation (6) and N is the dividing number.
u R ¯ = u ( R i ) N = u ( R i ) 8

3.2. Measurement Uncertainty Evaluation

The standard uncertainty of each uncertainty factor can be analyzed and obtained according to GUM.
  • Uncertainty coefficient of the Z-axis coordinate zi.
(1)
Uncertainty of stylus calibration in Z-axis direction u(ecalibration_Z).
The calibration of the employed Form Talysurf PGI420 is conducted with the dedicated master sphere shown in Figure 4a. The Pt value, which is the minimum and maximum width after linear expansion of the arc measured at the known radius values shown in Figure 4b, includes the shape error of the master ball, the shape error of the stylus tip diameter, the electrical noise of the system, and mechanical and environmental disturbances. The Pt value = 0.0337 μm exists in the radial direction, but assuming that it exists in the Z-axis direction and the X-axis direction, respectively, the uncertainty of the probe calibration in the Z-axis direction can be expressed by the following equation.
u e c a l i b r a t i o n _ Z = 0.0337   μ m = 33.7 n m
(2)
Uncertainty due to stylus resolution u(eresolution).
The stylus of the Form Talysurf PGI420 has a vertical resolution of 3.2 nm. Since the uncertainty due to the resolution of the probe is considered to be a square distribution of ±1.6 nm, it is obtained by the following equation using B-type evaluation.
u e r e s o l u t i o n = 3.2 / 2 3 = 0.92   n m
(3)
Uncertainty of repeated measurements u(erepeat).
The uncertainty due to the measurement error can be evaluated by the results of the repeated measurements of the small cylinder with a diameter 1.5 mm. Since it has a repeatability of 59.93 nm from the results of 10 repeated measurements, the uncertainty of repeated measurements can be obtained by the following formula using A type evaluation.
u e r e p e a t = 59.93 10 = 18.95   n m
2.
Uncertainty coefficient of the X-axis coordinate xi.
(1)
Uncertainty of stylus calibration in the X-axis direction u(ecalibration_X).
When picking up the output zi in the Z-axis direction as the uncertainty coefficient, it is assumed that Pt values exist in the Z-axis and X-axis directions, respectively, so the uncertainty of the stylus calibration in the X-axis direction is the same as that in the X-axis direction, expressed by the following equation.
u e c a l i b r a t i o n _ X = 33.7   n m
(2)
Uncertainty due to position error around the Z-axis of the workpiece u(ealignment_Z).
The uncertainty due to the workpiece position error can be obtained according to B type evaluation by using the value of the calculated measurement error, which can be obtained from the mathematical modeling. The influence of the position error around the Z-axis of the measured cylinder is shown in Figure 5a, in which D is the diameter of the workpiece, R is the radius, (xedge, zedge) is the end point coordinate of the obtained arc, and θZ is the angle of the position error around the Z-axis. The position error around the Z-axis causes the measurement result of the geometric circle of the geometric cylinder to appear as an ellipse, with its main axis in the X-axis direction. Since the arc of φ = 85°, that is, the arc of height h from the vertex shown in Equation (17), is extracted regardless of the size of the workpiece diameter, an error occurs in the X-axis coordinate. If the error angle is ±θZ° and this error has a rectangular distribution, the uncertainty due to the position error around the Z-axis can be represented by the following equation.
h = R × ( 1 cos φ 2 )
u e a l i g n m e n t _ Z = x e d g e cos θ Z 3 x e d g e = R sin φ 2 ( 1 cos θ Z 3 1 )
In the actual measurement, the alignment of the rotating stage and the scanning of the stylus were repeated to converge the obtained arc length to the minimum, and the measurement was performed with the error angle θZ being as small as possible. Assuming that the error angle at this time is ±1°, the uncertainty due to the attitude error around the Z-axis can be obtained as follows, by substituting the radius after stitching into Equation (18).
u e a l i g n m e n t _ Z = 0.74940 × sin 85 ° 2 1 cos 1 ° 3 1 = 0.00002568   m m = 25.68   n m
(3)
Uncertainty due to position error around the X-axis of the workpiece u(ealignment_Z).
The influence of the position error around the X-axis of the measured cylinder is shown in Figure 5b, in which D is the workpiece diameter, R is the radius, the coordinates (xedge, zedge) are the end point of the obtained arc, and θX is the angle of position error, which causes measurements of geometric circles to geometric cylinders to appear as ellipses, with the major axis along the Z-axis. Since the arc with height h from the vertex is extracted in the same way as when there is a posture error around the Z-axis, an error occurs in the X-axis coordinate. Assuming that the error angle is θX, which forms a rectangular distribution, the uncertainty due to the position error around the X-axis is expressed by the following equation.
u e a l i g n m e n t _ X = x e d g e x e d g e cos θ X 3 = R sin φ 2 ( 1 cos θ X 3 )
In the actual measurement, the manual stage was moved in the Y-axis direction while the stylus was in contact with the flat surface of the mounting table, and alignment was performed using the tilting stage so that the change in displacement output was minimized. As a result, the error angle θX was measured as small as possible. Assuming that the error angle at this time is ±0.1°, the uncertainty due to the attitude error around the X-axis is obtained as follows, by substituting the radius after stitching processing into Equation (20).
u e a l i g n m e n t _ X = 0.74940 × sin 85 ° 2 1 cos 1 ° 3 = 0.00000257   m m = 0.26   n m
Table 1 and Table 2 show a summary of the calculation results for each uncertainty. From these results, the combined standard uncertainty of the Z-axis direction output zi and the combined standard uncertainty of the X-axis coordinate xi are expressed by Equation (21).
u z i = u 2 e c a l i b r a t i o n _ Z + u 2 e r e s o l u t i o n + u 2 e r e p e a t = 33.7 2 + 0.92 2 + 18.95 2 = 38.45   n m
u x i = u 2 e c a l i b r a t i o n _ X + u 2 e a l i g n m e n t _ Z + u 2 e a l i g n m e n t _ X = 33.7 2 + 25.68 2 + 0.26 2 = 42.37   n m
Using the calculated uncertainties, the combined standard uncertainty u(Ri) of the radius of the measured arc can be obtained from Equation (23). u(Ri) can be obtained by Equation (24).
u r θ = u m θ = 26.92   n m
u R i = x i 2 z i 2 + x i 2 u 2 x i + z i 2 z i 2 + x i 2 u 2 z i = 0.50575 2 0.55280 2 + 0.50575 2 × 42.37 2 + 0.55280 2 0.55280 2 + 0.50575 2 × 38.45 2 = 40.28   n m
The uncertainty of the radius after the stitching process can be obtained by substituting the values obtained in Equation (13) into Equation (11) to obtain the following equation.
u R ¯ = u ( R i ) 8 = 40.28 8 = 14.24   n m
The uncertainty of the diameter measurement is twice the uncertainty of the radius.
u D = 2 × u R i = 2 × 14.24 = 28.48   n m
The expanded uncertainty can be calculated when k = 2.
U D = k u D = 2 × 28.48 = 56.96   n m
As can be seen above, the uncertainty of the diameter measurement of the 1.5 mm small cylindrical workpiece using the linear scanning method was estimated to be ±47.24 nm, achieving a target measurement uncertainty within ±0.1 μm. Table 3 summarizes the results of the standard uncertainty.
Using the calculated uncertainty, the combined standard uncertainty uri) of the radial deviation of an arbitrary measured point on the arc can be obtained from Equation (10). Substituting the combination of (xi, mi), the coordinates of the endpoint of the arc (xedge, zedge) = (0.50575, 0.55280) and the radius after the stitching process, for which u obtained from the measurement uri) is maximum, the following equation is obtained.
u r i = x i 2 + z i 2 2 R ¯ 2 2 u 2 R ¯ + x i R ¯ 2 u 2 x i + 1 z i R ¯ 2 u 2 z i = 0.50575 2 + 0.55280 2 2 × 0.74940 2 2 × 14.24 2 + 0.50575 0.74940 2 × 42.37 2 + 1 0.55280 0.74940 2 × 38.45 2 = 34.84   n m
Using the same assumptions as in Equations (6) and (29), the combined standard uncertainty u(Δzq) of roundness can be expressed as follows.
u z q = 2 u ( r ( θ ) )
u z q = 2 u ( r i )
Substituting the uncertainty calculated in Equation (26), the combined standard uncertainty u(Δzq) is obtained from Equation (27), as shown in Equation (31) below.
u z q = 2 u r i = 2 × 34.84 = 69.68   n m
The expanded uncertainty can be obtained when k = 2 shown in Equation (32).
U z q = k u z q = 2 × 69.68 = 139.36   n m
As can be seen above, the uncertainty in the roundness measurement of a small cylindrical workpiece with a diameter of 1.5 mm was estimated to be ±95.44 nm using the linear scanning method, achieving a target measurement uncertainty of ±0.1 μm or less. Table 3 summarizes the results of the standard uncertainty.

3.3. Variation of Uncertainty Due to Change in Workpiece Diameter

For the purpose of examining the cylindrical work diameter suitable for roundness measurement using the profilometer, which is the stitching linear scan method, a simulation is performed to calculate the measurement uncertainty when the work diameter is changed. Of the uncertainty factors in (1) to (3) of 1 and (1) to (3) of 2, shown in the previous section, the ones that change significantly depending on the workpiece diameter are the uncertainty u(ealignment_Z) due to the workpiece orientation error around the Z-axis and the uncertainty u(ealignment_X) due to the workpiece orientation error around the X-axis. Figure 6a,b show the graphs of the uncertainty u(xi) of the X-axis coordinate when the workpiece diameter is varied from 0.01 mm to 50.00 mm using Equations (13) and (14). From Figure 6a, as the diameter increases, u(ealignment_Z) increases and becomes almost the same value as u(xi), but the amount of change in u(ealignment_X) is small. In other words, u(ealignment_Z) is the dominant factor for changes in u(xi). This is because the attitude error angles θZ = ±1° and θX = ±0.1° are set in consideration of the alignment method and stylus resolution, and θZ is estimated to be larger than θX. Figure 7a shows a graph showing changes in diameter expanded uncertainty U(D) and roundness expanded uncertainty U(Δzq) created by substituting u(xi), which varies depending on the workpiece diameter, and other standard uncertainties. It can be seen that both U(D) and U(Δzq) increase as the workpiece diameter increases due to the influence of u(xi). Figure 7b shows an enlarged graph of the work diameter range from 0.01 mm to 10 mm. From this graph, when the attitude error angle θZ = ±1° and θX = ±0.1°, it can be read that the conditions for the workpiece diameter that can theoretically achieve the target measurement uncertainty within ±0.1 μm are a U(D) of 5.58 mm or less and a U(Δzq) of 2.11 mm or less.
In the previous simulation, the position error angles were set to θZ = ±1° and θX = ±0.1°, but improving the alignment increases the upper limit of the workpiece diameter that satisfies the target measurement uncertainty. There is room for improvement in the position error angle θZ around the Z-axis. Therefore, Figure 8a,b show the changes in the expanded uncertainty of diameter U(D) and the expanded uncertainty of roundness U(Δzq) when the attitude angle error occurs in a rectangular distribution of θZ = ±0.1°, ±0.3°, ±0.5°, ±0.6°, ±0.7°, ±0.8°, ±0.9°, ±1°, ±1.5° and ±2°. Table 4 summarizes the upper limit of the workpiece diameter that can achieve the target uncertainty of diameter and roundness within ±0.1 μm at each attitude error angle θZ. From Table 4, we see that it is desirable to converge θZ as small as possible by more accurate alignment when measuring cylindrical workpieces with a large diameter. It can be confirmed that the proposed method can cover all diameters of less than ϕ 3 mm, which are difficult to measure using the rotational scanning measurement method, as long as θZ is kept smaller than ±0.8°. Furthermore, if θZ can be made smaller, the proposed method can be used for workpieces with a larger diameter. If alignment is possible up to θZ = θX = ±0.1°, there is no limitation on the workpiece diameter theoretically. However, as the workpiece diameter increases, the measurement range in the Z-axis and X-axis directions expands, so it is thought that the upper limit of the workpiece diameter will be reached due to the increase in uncertainty and the limitation of the measurement range of the measuring equipment used.

3.4. Discussion

For the precision roundness measurement of cylindrical parts, the conventional rotary scan method cannot meet the requirements. Therefore, the stitching linear scan method has been proposed. Since the repeatability of measurement is very important, large numbers of experiments are necessary [28]. For practical applications, some cylindrical parts, such as the motor cylinders [29], can be measured by the stitching linear scan method. Meanwhile, the reliability of the proposed method can be verified in practical application. Contact measurements with a stylus or probe may not be accurate due to the measuring force and deformations, therefor it is necessary to develop non-contact measurements with optical linear encoder technology [30,31]. Furthermore, the proposed method and the verification manner can be applied for the measurement or testing of cylindrical parts or round parts made of other materials, such as functionally graded plates, which play a significant role in building materials [32,33]. Measurement uncertainty is necessary for testing, inspection and precision measurement. For the verification of the reliability of the proposed method, the measurement of uncertainty is carried out in this paper. The upper limitation of the proposed method can be investigated by the analysis and applications in a variety of fields.

4. Conclusions

The roundness profile of the small cylindrical parts can be reconstructed by coordinate modeling in the stitching linear scan method. To verify the reliability of the method, the measurement uncertainty of a small cylinder with a diameter 1.5 mm and length 7.8 mm was evaluated to be 0.095 μm, which satisfies the target uncertainty of 0.1 μm. According to the analysis, it can be determined that the uncertainty due to the position error around the Z-axis is the main uncertainty source of the X coordinate, after which the measurement uncertainties of the workpieces with a diameter from 0.01 mm to 50 mm can be analyzed. As a result, all the small cylinders with the dimensions mentioned above can be measured with a measurement uncertainty less than 0.1 μm when the alignment is improved properly. Theoretically, there is no limitation to the stitching linear scan method, while the employed profilometer has a measuring limitation. Therefore, a machine with a large measuring range will be used for investigating the limitation of the stitching linear scan method for roundness measurements of cylindrical workpieces. In addition, the measuring efficiency by the stitching linear scan method should be enhanced as well.

Author Contributions

J.Z.: methodology; resources; formal analysis; writing—original draft. L.Z.: investigation; data curation; review and editing. D.W.: review and editing; project administration. B.S.: data curation; visualization. Q.L.: investigation; visualization. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Defense Basic Scientific Research Program of China (Grant numbers JCKY2019427D002).

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qiu, Z.; Xue, J. Review of Performance Testing of High Precision Reducers for Industrial Robots. Measurement 2021, 183, 109794. [Google Scholar] [CrossRef]
  2. Yu, Z.; Hao, Z. Calibration of the Angle Measurement Error Caused by the Industrial Reducer Performance Test Instrument Torsional Deformation. Sci. Rep. 2022, 12, 21742. [Google Scholar] [CrossRef] [PubMed]
  3. Xu, L.X.; Chen, B.K.; Li, C.Y. Dynamic Modelling and Contact Analysis of Bearing-Cycloid-Pinwheel Transmission Mechanisms Used in Joint Rotate Vector Reducers. Mech. Mach. Theory 2019, 137, 432–458. [Google Scholar] [CrossRef]
  4. Yue, H.; Wu, X.; Shi, Z.; Zhang, Y.; Ye, Y.; Zhang, L.; Fu, Y. A Comprehensive Cycloid Pin-Wheel Precision Reducer Test Platform Integrated with A New Dynamic Measurement Method of Lost Motion. Metrol. Meas. Syst. 2022, 29, 207–229. [Google Scholar] [CrossRef]
  5. Xie, Y.H.; Xu, L.X.; Deng, Y.Q. A Dynamic Approach for Evaluating the Moment Rigidity and Rotation Precision of a Bearing-Planetary Frame Rotor System Used in RV Reducer. Mech. Mach. Theory 2022, 173, 104851. [Google Scholar] [CrossRef]
  6. Gao, W. Precision Nanometrology; Springer Series in Advanced Manufacturing; Springer London: London, UK, 2010; ISBN 978-1-84996-253-7. [Google Scholar]
  7. Gao, W. (Ed.) Metrology; Precision Manufacturing; Springer: Singapore, 2019; ISBN 978-981-10-4937-8. [Google Scholar]
  8. Gao, W. Surface Metrology for Micro-and Nanofabrication; Elsevier: Oxford, UK, 2020; ISBN 9780128178508. [Google Scholar]
  9. Gao, W. Precision Nanometrology and Its Applications to Precision Nanosystems. Int. J. Precis. Eng. Manuf. 2005, 6, 14–20. [Google Scholar]
  10. Wei, G.; Yuki, S.; Hane, K.; Soyama, H.; Adachi, K. Measurement and Instrumentation; Elsevier: Amsterdam, The Netherlands, 2017. [Google Scholar]
  11. Gao, W.; Shimizu, Y. Optical Metrology for Precision Engineering; De Gruyter: Berlin, Germany, 2021; ISBN 9783110542363. [Google Scholar]
  12. ISO/TC 213; ISO 12181-1; 2003 Geometrical Product Specifications (GPS)-Roundness-Part 1: Vocabulary and Parameters of Roundness. American National Standards Institute: Washington, DC, USA, 2007.
  13. Taylor Hobson Ltd. A Guide to the Measurement of Roundness. Introduction to Roundness. Available online: http://www.tarkkuustuonti.fi/Kampanjat/Brochure_Roundness_Booklet.pdf (accessed on 6 July 2023).
  14. Taylor Hobson Ltd. Roundness Measurement Equipment|Form Measurement|Cylindricity Measuring Instrument|Roundness Tester. Available online: https://www.taylor-hobson.com/products/roundness-form (accessed on 6 July 2023).
  15. Sui, W.; Zhang, D. Four Methods for Roundness Evaluation. Phys. Procedia 2012, 24, 2159–2164. [Google Scholar] [CrossRef] [Green Version]
  16. Bai, J.; Wang, Y.; Wang, X.; Zhou, Q.; Ni, K.; Li, X. Three-Probe Error Separation with Chromatic Confocal Sensors for Roundness Measurement. Nanomanufacturing and Metrology. Nanomanufacturing Metrol. 2021, 4, 247–255. [Google Scholar] [CrossRef]
  17. Gao, W.; Kiyono, S. On-Machine Roundness Measurement of Cylindrical Workpieces by the Combined Three-Point Method. Measurement 1997, 21, 147–156. [Google Scholar] [CrossRef]
  18. Cai, Y.; Xie, B.; Ling, S.; Fan, K.-C. On-Line Measurement Method for Diameter and Roundness Error of Balls. Nanomanufacturing Metrol. 2020, 3, 218–227. [Google Scholar] [CrossRef]
  19. Gao, W.; Haitjema, H.; Fang, F.Z.; Leach, R.K.; Cheung, C.F.; Savio, E.; Linares, J.M. On-Machine and in-Process Surface Metrology for Precision Manufacturing. CIRP Ann. 2019, 68, 843–866. [Google Scholar] [CrossRef] [Green Version]
  20. Chen, Y.L.; Machida, Y.; Shimizu, Y.; Matsukuma, H.; Gao, W. A Stitching Linear-Scan Method for Roundness Measurement of Small Cylinders. CIRP Ann. 2018, 67, 535–538. [Google Scholar] [CrossRef]
  21. Li, Q.; Shimizu, Y.; Saito, T.; Matsukuma, H.; Gao, W. Measurement Uncertainty Analysis of a Stitching Linear-Scan Method for the Evaluation of Roundness of Small Cylinders. Appl. Sci. 2020, 10, 4750. [Google Scholar] [CrossRef]
  22. Li, Q.; Shimizu, Y.; Saito, T.; Matsukuma, H.; Cai, Y.; Gao, W. Improvement of a Stitching Operation in the Stitching Linear-Scan Method for Measurement of Cylinders in a Small Dimension. Appl. Sci. 2021, 11, 4705. [Google Scholar] [CrossRef]
  23. JCGM 100:2008; GUM 1995 with Minor Corrections-Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement; JCGM: Sèvres, France, 2008; Volume 50.
  24. JCGM 200:2012; International Vocabulary of Metrology–Basic and General Concepts and Associated Terms (VIM)-2008 Version with Minor Corrections. BIPM: Sèvres, France, 2012; Volume 3.
  25. Wen, X.-L.; Zhao, Y.-B.; Wang, D.-X.; Pan, J. Adaptive Monte Carlo and GUM Methods for the Evaluation of Measurement Uncertainty of Cylindricity Error. Precis. Eng. 2013, 37, 856–864. [Google Scholar] [CrossRef]
  26. Jermak, C.J.; Jakubowicz, M.; Derezynski, J.; Rucki, M. Uncertainty of the Air Gauge Test Rig. Int. J. Precis. Eng. Manuf. 2017, 18, 479–485. [Google Scholar] [CrossRef]
  27. Gao, W.; Kim, S.W.; Bosse, H.; Haitjema, H.; Chen, Y.L.; Lu, X.D.; Knapp, W.; Weckenmann, A.; Estler, W.T.; Kunzmann, H. Measurement Technologies for Precision Positioning. CIRP Ann. Manuf. Technol. 2015, 64, 773–796. [Google Scholar] [CrossRef]
  28. Skrzek, T.; Rucki, M.; Górski, K.; Matijošius, J.; Barta, D.; Caban, J.; Zarajczyk, J. Repeatability of High-Pressure Measurement in a Diesel Engine Test Bed. Sensors 2020, 20, 3478. [Google Scholar] [CrossRef]
  29. Jermak, C.J.; Dereżyński, J.; Rucki, M. Measurement System for Assesment of Motor Cylinder Tolerances and Roundness. Metrol. Meas. Syst. 2023, 25, 103–114. [Google Scholar] [CrossRef]
  30. Gurauskis, D.; Przystupa, K.; Kilikevičius, A.; Skowron, M.; Matijošius, J.; Caban, J.; Kilikevičienė, K. Development and Experimental Research of Different Mechanical Designs of an Optical Linear Encoder’s Reading Head. Sensors 2022, 22, 2977. [Google Scholar] [CrossRef]
  31. Kowalik, M.; Rucki, M.; Paszta, P.; Gołębski, R. Plastic Deformations of Measured Object Surface in Contact with Undeformable Surface of Measuring Tool. Meas. Sci. Rev. 2016, 16, 254–259. [Google Scholar] [CrossRef] [Green Version]
  32. Tounsi, A.; Tahir, S.; Al-Osta, M.; Do-Van, T.; Bourada, F.; Bousahla, A.A.; Tounsi, A. An Integral Quasi-3D Computational Model for the Hygro-Thermal Wave Propagation of Imperfect FGM Sandwich Plates. Comput. Concr. 2023, 32, 61–74. [Google Scholar]
  33. Tounsi, A.; Mostefa, A.H.; Attia, A.; Bousahla, A.A.; Bourada, F.; Tounsi, A.; Al-Osta, M.A. Free Vibration Investigation of Functionally Graded Plates with Temperature-Dependent Properties Resting on a Viscoelastic Foundation. Struct. Eng. Mech. 2023, 86, 1–16. [Google Scholar]
Figure 1. (a) Linear scan with a stylus; (b) small cylinder mounted on the V-groove and scanned linearly.
Figure 1. (a) Linear scan with a stylus; (b) small cylinder mounted on the V-groove and scanned linearly.
Applsci 13 09091 g001
Figure 2. The stitching procedure: (a) profile 1; (b) profile 2; (c) the stitching of profiles 1 and 2; (d) the stitching of the rest of the arc profiles.
Figure 2. The stitching procedure: (a) profile 1; (b) profile 2; (c) the stitching of profiles 1 and 2; (d) the stitching of the rest of the arc profiles.
Applsci 13 09091 g002
Figure 3. (a) Schematic of the radius measurement and the uncertainty components; (b) schematic of the roundness measurement and the uncertainty components.
Figure 3. (a) Schematic of the radius measurement and the uncertainty components; (b) schematic of the roundness measurement and the uncertainty components.
Applsci 13 09091 g003
Figure 4. (a) Standard ball; (b) calibration result.
Figure 4. (a) Standard ball; (b) calibration result.
Applsci 13 09091 g004
Figure 5. (a) Influence of attitude error of roll workpiece around Z-axis; (b) influence of attitude error of roll workpiece around X-axis.
Figure 5. (a) Influence of attitude error of roll workpiece around Z-axis; (b) influence of attitude error of roll workpiece around X-axis.
Applsci 13 09091 g005
Figure 6. Variation of u(xi) in accordance with the diameter of workpiece: (a) influence of uealignment_Z); (b) influence of uealignment_X).
Figure 6. Variation of u(xi) in accordance with the diameter of workpiece: (a) influence of uealignment_Z); (b) influence of uealignment_X).
Applsci 13 09091 g006
Figure 7. Variation of U(D) and U(Δzq) in accordance with diameter of workpiece: (a) diameter ф 0.01–50 mm; (b) diameter ф 0.01–10 mm.
Figure 7. Variation of U(D) and U(Δzq) in accordance with diameter of workpiece: (a) diameter ф 0.01–50 mm; (b) diameter ф 0.01–10 mm.
Applsci 13 09091 g007
Figure 8. (a) Variation in U(D) with the diameter of the workpiece; (b) variation in U(Δzq) in accordance with the diameter of workpiece.
Figure 8. (a) Variation in U(D) with the diameter of the workpiece; (b) variation in U(Δzq) in accordance with the diameter of workpiece.
Applsci 13 09091 g008
Table 1. Uncertainty budget of the output in the Z-axis direction.
Table 1. Uncertainty budget of the output in the Z-axis direction.
Source of
Uncertainty
SymbolTypeCoverage FactorStandard UncertaintySensitivity Coefficient|ci| × u(xi)
nm
Calibration
of probe
u(ecalibration_Z)A33.7133.7
Resolutionu(eresolution)B 3 0.9210.92
Repeatabilityu(erepeat)A18.95118.95
Combined standard
uncertainty
u(zi)38.45
Table 2. Uncertainty budget of the X-axis coordinate.
Table 2. Uncertainty budget of the X-axis coordinate.
Source of
Uncertainty
SymbolTypeCoverage FactorStandard
Uncertainty
Sensitivity Coefficient|ci| × u(xi)
nm
Calibration
of probe
u(ecalibration_X)A33.7133.7
Attitude error
around
Z-axis
u(ealignment_Z)B 3 25.68125.68
Attitude error
around
X-axis
u(ealignment_X)B 3 0.2610.26
Combined standard uncertaintyu(xi)42.37
Table 3. Uncertainty budget of diameter measurement by the proposed method.
Table 3. Uncertainty budget of diameter measurement by the proposed method.
Source of
Uncertainty
SymbolTypeCoverage FactorStandard UncertaintySensitivity Coefficient|ci| × u(xi)
nm
Output in
Z-axis direction
u(zi)38.45138.45
Coordinate of
X-axis
u(xi)42.37142.37
Combined
standard uncertainty
u(Ri)40.28
Radius after
stitching process
u ( R ¯ ) 14.24
Diameteru(D)28.48
Expanded
uncertainty (k = 2)
U(D)56.96
Table 4. Limits of diameter that can be achieved.
Table 4. Limits of diameter that can be achieved.
Attitude Error Angle
θZ deg.
0.10.30.50.60.70.80.911.52
Work Diameter Upper Limit mm (U(D))-61.6522.3115.5011.398.726.895.582.481.39
Work Diameter Upper Limit mm (U(Δzq))-23.348.445.864.313.302.602.110.930.52
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, J.; Zhang, L.; Wu, D.; Shen, B.; Li, Q. Measurement Uncertainty Analysis of the Stitching Linear-Scan Method for the Measurable Dimension of Small Cylinders. Appl. Sci. 2023, 13, 9091. https://doi.org/10.3390/app13169091

AMA Style

Zhao J, Zhang L, Wu D, Shen B, Li Q. Measurement Uncertainty Analysis of the Stitching Linear-Scan Method for the Measurable Dimension of Small Cylinders. Applied Sciences. 2023; 13(16):9091. https://doi.org/10.3390/app13169091

Chicago/Turabian Style

Zhao, Jiali, Liang Zhang, Dan Wu, Bobo Shen, and Qiaolin Li. 2023. "Measurement Uncertainty Analysis of the Stitching Linear-Scan Method for the Measurable Dimension of Small Cylinders" Applied Sciences 13, no. 16: 9091. https://doi.org/10.3390/app13169091

APA Style

Zhao, J., Zhang, L., Wu, D., Shen, B., & Li, Q. (2023). Measurement Uncertainty Analysis of the Stitching Linear-Scan Method for the Measurable Dimension of Small Cylinders. Applied Sciences, 13(16), 9091. https://doi.org/10.3390/app13169091

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop