Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Interferometric method of measuring external cylindrical surfaces

Open Access Open Access

Abstract

A new method of measuring external cylindrical surfaces is proposed, making use of a ring beam transform unit. The unit is composed of three cone mirrors, with which a parallel beam can be transformed into annular and convergent beams. The method has advantages over classic methods: It avoids physical contact with the cylinder surface and is fast and stitching-free. In experiments using a dynamic interferometer to demonstrate the feasibility of the method, the RMS difference of the axis contour from results using a Luphoscan is 0.0282 µm, while the roundness difference of the circumferential contour is 0.0962 µm.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. INTRODUCTION

Cylindrical surfaces are widely used in mechanical and optical systems such as aerostatic cylindrical bearings, error calibration of rotary stages, and interferometry systems. In most such systems, the profile of the cylindrical surface directly affects the performance of the system. Traditional instruments to measure cylindrical surfaces, e.g., a coordinate measurement machine, cylindricity measuring instrument, and contour graph with V-block, make contact with the surface, possibly resulting in surface damage in the measurement process [1,2]. Moreover, because these instruments scan the cylinder circle by circle along its axis, rotary state errors may be introduced into the results. With the development of computer-generated holograms (CGH) and subaperture stitching technology, interferometry methods have been proposed to measure cylindrical surfaces. Peng et al. sought to achieve high-precision measurements, including the introduction of Legendre–Fourier polynomials, to offset the form error by adopting a cylindrical projection to establish one-to-one mapping between the subaperture data and points on the cylinder surface [35]. However, stitching technology is inefficient because dozens of subapertures are needed, owing to the limited $F$-number of CGH. Accordingly, other methods have been proposed to achieve a higher measurement efficiency without sacrificing accuracy. Weckenmann et al. proposed a grazing-incidence interferometer equipped with CGH to measure cylindrical surfaces. This incorporated parallel data acquisition of the cylindrical surfaces and achieved excellent radial resolution [6]. White light interferometers have been used with conic mirrors to measure both continuous and stepped cylindrical surfaces, achieving an uncertainty as small as 1 µm [7,8]. Also, deflectometry in cylindrical coordinates has been proposed to measure the cross section of cylindrical parts and determine deviations from roundness [9].

In this paper, a modified Twyman–Green interferometry system is proposed to realize efficient measurement of the external cylindrical surface. The key element of the system is a ring beam transform unit (RBTU) composed mainly of three conic mirrors, which transform the parallel beam into annular and convergent beams, so that information on the cylindrical surface can be reflected back to form useful interferometric fringes. The system is realized with dynamic interferometry, and the results demonstrate the feasibility and validity of the proposed method.

2. METHOD

A. Measurement System

The setup for the proposed method is shown in Fig. 1. It consists of a laser source, beam expander, reference mirror, RBTU, detector, and other lenses. The beam expander is optional because the size of the beam determines the height of the cylinder accessible for a single measurement. The laser beam is split by the polarizing beam splitter (PBS) into measuring and reference beams. The measuring beam is then transformed into annular and convergent beams by means of three reflections from conic mirrors A, B, and C. Reflected by the sample cylinder, which is set to be coaxial with the measured system, the measuring beam returns to the PBS carrying surface information and forms interference fringes with the reference beam on the detector.

 figure: Fig. 1.

Fig. 1. Schematic of measuring system.

Download Full Size | PDF

Using the interference fringes and phase-stepping interferometry, the surface information can be obtained. Let ($x,y$) be the coordinates of a point on an interference fringe; the height ($x, y$) is the surface deviation at ($x,y$), and ${{ I}_1} - {{ I}_4}$ are the intensities at ($x, y$) for phase shifts of 0, ${0.5}\pi$, $\pi$, and ${1.5}\pi$, respectively. ${{I}_1} - {{ I}_4}$ can be obtained by the four interference fringes captured by the camera, and the surface deviation can be described by

$$\text{Height}(x,y) = \frac{\lambda}{{4\pi}}\arctan\left[{\frac{{{I_4} - {I_2}}}{{{I_1} - {I_3}}}} \right],$$
where $\lambda$ is the wavelength of the laser (632.8 nm in this paper) [10]. In this way, all the surface information in the measured range can be obtained.

B. Surface Reconstruction

The surface deviation results appear on a circular plane since the purpose is to measure surface information on a cylindrical surface; hence, surface reconstruction is needed to establish the relationship between the surface deviation and the cylindrical surface. Moreover, the eccentricity and tilt of the cylindrical surface relative to the RBTU directly affect the surface deviation, which needs to be corrected in the reconstruction process in order to obtain a more accurate result. As shown in Fig. 2, the reconstruction can be divided into three steps: (1) transforming the measured results to the 3D cylindrical surface; (2) fitting the surface with a spatial cylinder to remove the tilt and eccentricity between the sample and the RBTU; and (3) computing the distance between the measured points and the axis of the fitted cylindrical surface to check on the cylindricity.

 figure: Fig. 2.

Fig. 2. Surface reconstruction process. (a) Transforming the surface deviation result to a 3D cylindrical surface. (b) Fitting the surface with a spatial cylinder. (c) Recomputing the cylindricity.

Download Full Size | PDF

According to the geometric relationships, the radial direction of the results corresponds to the axial direction of the cylinder. Let O be the center of the surface deviation results. Then, (${x_0}, {y_0}$) are the pixel coordinates of a point P, $\theta$ is the angle between OP and OX, as shown in Fig. 2(a), $d$ is the pixel size, and ${{\rm P}_1}$ is the projection of point P on the cylinder. The location of ${{\rm P}_1}$ can be described by

$$\left[ \begin{array}{c} x \\ y \\ z \\ \end{array} \right]=\left[ \begin{array}{c} (\text{Height}({{x}_{0}},{{y}_{0}})+\text{R)}\cdot \text{cos}\theta \\[4pt] (\text{Height}({{x}_{0}},{{y}_{0}})+\text{R)}\cdot \text{sin}\theta \\[4pt] d\cdot \sqrt{x_{0}^{2}+y_{0}^{2}} \\ \end{array} \right],$$
where $R$ is a radius artificially added to reconstruct the cylinder shape. It has no effect on the evaluation of the cylindricity because it is equivalent to the radial offset of the whole surface. The theoretical radius of the measured cylinder surface is usually designated $R$.
 figure: Fig. 3.

Fig. 3. Sample experimental device. (a) Experimental system. (b) Layout of RBTU. (c) Size and measured region of cylindrical surface. (d) Image of cylindrical surface sample.

Download Full Size | PDF

Let $q$ designate the axial line of a cylinder. ${{\rm P}_0}$ (${x_1}, {y_1}, {z_1}$) is a point on $q$, and ($l, m, n$) specify the direction vector of $q$, with $r$ being the radius of the cylinder. According to the definition of a cylinder, the equation of a cylinder can be described by Eq. (3).

Combining the point cloud found by Eq. (2) and the equation of a cylinder, the parameters of the cylinder can be obtained by minimizing Eq. (4) according to the nonlinear least-square method [11]. Then, the surface information corrected for eccentricity and tilt can be obtained by computing the distance between each point and the axis of the cylinder:

$$\begin{split}&{(x - {x_1})^2} + {(y - {y_1})^2} + {(z - {z_1})^2} \\&\quad- \frac{{{{(l(x - {x_1}) + m(y - {y_1}) + n(z - {z_1}))}^2}}}{{{l^2} + {m^2} + {n^2}}} = {r^2},\end{split}$$
in which ${x_1}$, ${y_1}$, ${z_1}$, $l$, $m$, $n$, and $r$ are parameters to be estimated:
$$\sum\limits_i {{{\left(\sqrt {{{({x_i} - {x_1})}^2} + {{(y - {y_i})}^2} + {{({z_i} - {z_1})}^2} - \frac{{{{(l({x_i} - {x_1}) + m({y_i} - {y_1}) + n({z_i} - {z_1}))}^2}}}{{{l^2} + {m^2} + {n^2}}}} - r\right)}^2}} .$$

3. EXPERIMENT AND DISCUSSION

A. Experimental Setup

To verify the practicality of the proposed method, experiments based on dynamic interferometry were conducted. As shown in Fig. 3(a), the system included a dynamic interferometer (PhaseCam 6000 made by 4D Technology Corporation), a ring beam transform unit (RBTU), and several multiaxis adjustment stages. The camera of PhaseCam 6000 has an effective image size of ${996} \times {984}$ pixels, and the pixel size is 9.03 µm. The effective radius of the laser beam was 8.5 mm, which means the maximum theoretical measured range was 4.25 mm. To minimize the influence of the geometric position between three conic mirrors, the RBTU is specially designed, as shown in Fig. 3(b). Conic mirror A is glued to an optical flat with antireflection film on both sides, and the position can be adjusted by four screws. The relative position between conic mirrors B and C is ensured by the conic surface concentric to them. The cylindrical surface sample is an aluminum alloy part manufactured by single-point diamond turning with a diameter of 30 mm, in which an arc with 2 µm width and 5 mm length is manufactured as the measure region, as shown in Figs. 3(c) and 3(d).

The adjustment in the measurement process contains two main steps. (1) The sample is removed from the RBTU and the relative position of the RBTU and the dynamic interferometer is adjusted to make the optic axis parallel to the cylinder axis as well as to keep the laser at the center of the RBTU. In this step, the RBTU reflects the rays to the dynamic interferometer, so that the interference fringes indicate the relative position of the RBTU and the dynamic interferometer. (2) The sample is put back into the RBTU and adjusted, with the RBTU and dynamic interferometer fixed, until reasonable interference fringes appear.

B. Results and Discussion

Figure 4 shows the interference fringes of an experiment. There is a fuzzy part in the center due to (1) the top of conic mirror A in Fig. 1 can never be manufactured to its theoretical shape; (2) the lateral resolution along the angle theta is radius-dependent, and the effective surface information decreases with the decreases of radius. MATLAB is adopted to achieve the surface deviation computation, unwrapping, and surface reconstruction.

 figure: Fig. 4.

Fig. 4. Interference fringes of an experiment.

Download Full Size | PDF

Figure 5(a) shows the original surface deviation unwrapped by phase-stepping interferometry based on the interference fringes in Fig. 4. The axial resolution is the pixel size, while the lateral resolution is the quotient of the radius of the cylindrical sample and the radius in terms of the number of pixels on the camera sensor according to Fig. 2. In order to obtain sufficient sampling density and a complete cylindrical surface, an annular region between two concentric circles is established as the region of interest, in which the radius of the inner circle is 150 pixels and the radius of external circle is 450 pixels; hence, the axial measure range is $({450} - {150}) \times {9.03} = {2.709}\;{\rm mm}$. The reconstructed surface corrected for eccentricity and tilt of region of interest is shown in Fig. 5(b), in which the colored bar indicates the radius deviation of the measure point relative the estimated $r$ in Eq. (3), and the maximum surface deviation is 1.3 µm within the region of interest.

 figure: Fig. 5.

Fig. 5. Measured results. (a) Original results. (b) Reconstructed results.

Download Full Size | PDF

 figure: Fig. 6.

Fig. 6. Measured results with Luphoscan 420 and comparisons. (a) Photo of the Luphoscan 420. (b) Point cloud and comparison line.

Download Full Size | PDF

As shown in Fig. 6(a), a Luphoscan 420 manufactured by Taylor Hobson is adopted to measure the cylindrical surface sample in order to verify the measured results. Figure 6(b) shows the results measured in this way, in which the colored bar is the radius deviation obtained by Luphoscan and reveals a maximum surface deviation of 2.2 µm over a 5 mm length region.

 figure: Fig. 7.

Fig. 7. Results comparison. (a) Axial line comparison results. (b) Circumferential line comparison results.

Download Full Size | PDF

The proposed interferometry system and Luphoscan provide relative measurements not tied to the actual radius of the cylindrical surface. The contour lines adopted to compare the two methods must therefore be adjusted to the same reference value, which does not affect the comparison of surface deviations. Figure 7(a) shows comparative results for the axial contour line, for which the contours are integrally shifted and rotated. The RMS difference of the two contours is 0.0288 µm, indicating that the contours from the two methods match well. Figure 7(b) shows comparative results for the circumferential contour line, in which the two contours are adjusted to the reference radius of 1.65 µm. The roundness determined by Luphoscan is 0.1833 µm; it is 0.2795 µm for the proposed interferometric method.

4. CONCLUSIONS

A new interferometry method to measure radial deviation of an external cylindrical surface is presented. It makes use of a ring beam transform unit composed of three conic mirrors. An experiment using a dynamic interferometer demonstrates the feasibility of the method. In the experiment, the axial contour differs from the Luphoscan results by 0.0282 µm, and the roundness difference of the circumferential contour is 0.0962 µm. The proposed method can obtain the 360° surface map in a single measurement without physical contact, which improves measurement efficiency and avoids surface damage. However, the surface deviation of three conic mirrors is directly reflected in the results, though it is much smaller than the deviation of sample, which may be studied further to remove system errors and improve measurement accuracy.

Funding

Science Challenging Program, China Academy of Engineering Physics (JCKY2016212A506-0106).

Acknowledgment

The authors thank Shanghai Engineering Research Center of AI & Robotics and Engineering Research Center of AI & Robotics, Ministry of Education, for their support.

Disclosures

The authors declare no conflicts of interest.

REFERENCES

1. W. Q. Zhao, Z. Xue, J. B. Tan, and Z. B. Wang, “SSEST: A new approach to higher accuracy cylindricity measuring instrument,” Int. J. Mach. Tools Manuf. 46, 1869–1878 (2006). [CrossRef]  

2. Y. Chen, Y. Machida, Y. Shimizu, H. Matsukuma, and W. Gao, “A stitching linear-scan method for roundness measurement of small cylinders,” CIRP Ann. 67, 535–538 (2018). [CrossRef]  

3. J. Peng, Y. Yu, and H. Xu, “Compensation of high-order misalignment aberrations in cylindrical interferometry,” Appl. Opt. 53, 4947–4956 (2014). [CrossRef]  

4. J. Peng, Q. Wang, X. Peng, and Y. Yu, “Stitching interferometry of high numerical aperture cylindrical optics without using a fringe-nulling routine,” J. Opt. Soc. Am. A 32, 1964–1972 (2015). [CrossRef]  

5. J. Peng, Y. Yu, D. Chen, H. Guo, J. Zhong, and M. Chen, “Stitching interferometry of full cylinder by use of the first-order approximation of cylindrical coordinate transformation,” Opt. Express 25, 3092–3103 (2017). [CrossRef]  

6. A. Weckenmann, J. Bruning, S. Patterson, P. Knight, and J. Bryan, “Grazing Incidence Interferometry for High Precision Measurements of Cylindrical Form Deviations,” CIRP Ann. 50, 381–384 (2001). [CrossRef]  

7. A. G. Albertazzi and A. Dal Pont, “A white light interferometer for measurement of external cylindrical surfaces,” in 5th International Workshop on Automatic Processing of Fringe Patterns Fringe (Springer, 2005), pp. 632–639.

8. A. G. Albertazzi, M. R. Viotti, R. M. Miggiorin, and A. Dal Pont, “Applications of a white light interferometer for wear measurement of cylinders,” Proc. SPIE 7064, 70640B (2008). [CrossRef]  

9. M. J. F. Carvalho, C. L. N. Veiga, and A. Albertazzi, “Deflectometry in cylindrical coordinates using a conical mirror: principles and proof of concept,” J. Braz. Soc. Mech. Sci. Eng. 41, 380 (2019). [CrossRef]  

10. D. Malacara, “Phase shifting interferometry,” in Optical Shop Testing, 3rd ed. (Wiley, 2007), pp. 547–666.

11. G. Lukács, R. Martin, and D. Marshall, Faithful Least-Squares Fitting of Spheres, Cylinders, Cones and Tori for Reliable Segmentation, H. Burkhardt and B. Neumann, eds. (Springer Berlin Heidelberg, 1998), pp. 671–686.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. Schematic of measuring system.
Fig. 2.
Fig. 2. Surface reconstruction process. (a) Transforming the surface deviation result to a 3D cylindrical surface. (b) Fitting the surface with a spatial cylinder. (c) Recomputing the cylindricity.
Fig. 3.
Fig. 3. Sample experimental device. (a) Experimental system. (b) Layout of RBTU. (c) Size and measured region of cylindrical surface. (d) Image of cylindrical surface sample.
Fig. 4.
Fig. 4. Interference fringes of an experiment.
Fig. 5.
Fig. 5. Measured results. (a) Original results. (b) Reconstructed results.
Fig. 6.
Fig. 6. Measured results with Luphoscan 420 and comparisons. (a) Photo of the Luphoscan 420. (b) Point cloud and comparison line.
Fig. 7.
Fig. 7. Results comparison. (a) Axial line comparison results. (b) Circumferential line comparison results.

Equations (4)

Equations on this page are rendered with MathJax. Learn more.

Height ( x , y ) = λ 4 π arctan [ I 4 I 2 I 1 I 3 ] ,
[ x y z ] = [ ( Height ( x 0 , y 0 ) + R) cos θ ( Height ( x 0 , y 0 ) + R) sin θ d x 0 2 + y 0 2 ] ,
( x x 1 ) 2 + ( y y 1 ) 2 + ( z z 1 ) 2 ( l ( x x 1 ) + m ( y y 1 ) + n ( z z 1 ) ) 2 l 2 + m 2 + n 2 = r 2 ,
i ( ( x i x 1 ) 2 + ( y y i ) 2 + ( z i z 1 ) 2 ( l ( x i x 1 ) + m ( y i y 1 ) + n ( z i z 1 ) ) 2 l 2 + m 2 + n 2 r ) 2 .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.