Vibration Measurement Using a Fringe Pattern in Reflective Monochromatic Interferometry
 Author: Kim Minsu, Yoon DoYoung, Pahk Heuijae
 Publish: Journal of the Optical Society of Korea Volume 19, Issue5, p494~502, 25 Oct 2015

ABSTRACT
This paper introduces methods to measure vibration using a fringe pattern. These methods use variations of a fringe pattern in reflective monochromatic interferometry, without additional components. With the proposed methods we measured the vibrations of four waveform with amplitude 100 nm. When the vibrational amplitude is greater than a quarter wavelength of the light employed, however, the measured results are distorted due to ambiguity. Thus we propose advanced methods to solve this problem, and also measure the vibrations of two waveformswith an amplitude of 1 μm. To verify the performance of the proposed methods, we compare the results to those from an accelerometer. Multifrequency vibrations of 1, 5, 10, and 20 Hz are measured by both techniques, and the results compared in the frequency domain.

KEYWORD
Interferometry , Vibration measurement , Fringe analysis

I. INTRODUCTION
Interferometry, which offers high resolution and nondestructive and highspeed measurements, is widely used in precision metrology, such as for liquid crystal displays (LCDs) [1], semiconductors [2, 3], and microelectromechanical systems (MEMS). Interferometry acquires interferograms as it moves along the
z axis, and surface shapes are reconstructed from the acquired interferograms. Therefore the sensitivity of interferometry to outofplane vibration is based on the principle of its acquisition and reconstruction process. Because of its sensitivity to vibrations, they should be evaluated and managed to enhance interferometric performance.Accelerometers are widely used in industry to measure vibrations. However, they have some disadvantages, one being that the accelerometer needs to be attached to a vibrating part. Furthermore, the mass of the sensor can distort vibration measurements.
The most extensively studied methods for measuring vibration are speckle pattern interferometry (SPI) [48] and digital holographic interferometry (DHI) [911], which have been used to measure the vibration of the interferometry system itself. These methods use a laser as a light source and measure the deformation of the sample. The light source is a limitation for these systems.
Techniques for vibration measurement with Michelson interferometry have been studied in various ways. A method for realtime measurement of outofplane vibration using an additional beam splitter (BS) and two gratings has been proposed [12]. This system requires additional optical components, and it was difficult to align the gratings and reference mirror. Some problems also arose when measuring linear and sharp vibrations.
Interferometry with a selfpumped phase conjugate mirror (SPCCM) has been used to measure vibration [13]. Vibration on a rough surface could be measured by this system, but the range of measurable amplitude was too narrow, because of the properties of the SPCCM. Dual Michelson interferometry using optical fibers has been studied [14, 15]. These systems have the advantage of being comparatively insensitive to external noise and long measuring distances. However, they require an additional laser and a complicated feedback loop to compensate for external noise [14]. Several other ways have been studied, including the use of Moiré interferometry with gratings [16, 17] and frequencyscanned interferometry [18]. A vibrationmeasurement method using Michelson interferometry without any additional component has been studied [19]. This method obtains the intensity as a function of time and converts it to the frequency domain by fast Fourier transform. This method has the advantage of a simple setup and low costbut it only measures vibrational frequency, not vibrational amplitude or phase error.
In this paper, we propose methods for measuring displacement and phase error caused by vibration. These methods use variations of the fringe pattern in reflective monochromatic interferometry, without additional components. These methods can measure environmental vibrations, along with those of the interferometry system. The most important advantage of the proposed methods is that they can be easily applied to any kind of interferometry.
II. BASICS OF INTERFEROMETRY
A basic Mirau interferometry approach is shown in Fig. 1. Light from an LED or a laser is incident upon BS1 and is reflected to the objective lens. Then, the reflected beam is split by BS2 into sample and reference arms. The split beams are reflected from the sample and reference mirror respectively, and then the reflected beams pass through BS2 and BS1. Finally, the beams are detected by a chargecoupled device (CCD) camera. These two beams travel different paths and create an interferogram.
The interferogram frommonochromatic interferometry can be described as:
where
I_{DC} is the average background intensity as the DC term,γ is the visibility,k is the wave number,d _{1} is the distance from BS2 to a reference mirror, andd _{2} is the distance from BS2 to the sample.The difference in the distances that the two beams travel, 2(
d _{1}−d _{2}), is referred to as the optical path difference (OPD). The intensity of the interferogram varies between bright and dark according to the OPD. The intensities are acquired as the OPD is changed via mechanical moving parts. Then, the surface profile is recovered from the acquired intensities by phase reconstruction. The most widely used reconstruction method, the SchwiderHariharan algorithm [20, 21], requires five intensities as OPD changes byπ /2. It reconstructs the phase as shown below:where
However, outofplane vibrationcauses a change in OPD during mechanical adjustment. In the presence of vibration of a single frequency, Eq. (1) can be reexpressed as:
where
α is the amplitude andω _{c} is the angular frequency of the vibration.Vibrations cause uncontrolled changes in OPD, and this can induce phase errors in the reconstruction process. When the surface profile is restored, phase errors introduce ripples into the measured results for a flat surface [22]. Thus, it is important to measure the variations of the interference fringe caused by vibrations to optimize an interferometry system and compensate for phase errors.
III. PROPOSED MEASUREMENT METHODS FOR VIBRATION
In interferometry, several fringe patterns appear when the sample has a tilted flat surface. Variations in fringe patterns are detected as they move laterally due to outofplane vibration.
The OPD is changed when the sample moves from position A to position B due to outofplane vibration, as shown in Fig. 2. If a dark fringe appears at point 1 when the sample is placed at position A, it shifts to point 2 as the sample moves to position B. Therefore, the fringe pattern shifts a distance
ε_{x} when the displacement due to vibration is equal toε_{z} . Thus outofplane vibration can be measured using these variations of the fringe pattern. In other words, the influence of vibration on interferometry is measured directly using the variation in the fringe patterns.Outofplane vibration is the main source of measurement error among various types, including inplane vibration and tilt vibration. Inplane vibration does not change the OPD; therefore, the proposed method cannot measure inplane vibration. On the other hand, tilt vibration alters thetilt angle
β and also affects the proposed methods. We can determineβ as a function of time from the fringe spacingl_{x} and the wavelength of lightλ . Outofplane vibration, however, is a major source of measurement error, while there is little tilt vibration, because a typical interferometry setup is constructed on a massive, damped air table. Thus, we consider only the case of outofplane vibration in this paper.3.1. Displacement Measurement Method
Images are captured by CCD at a fixed time interval without any mechanical moving process, while the sample is vibrated vertically. We call the first captured image the
reference image, and the otherstarget images. Vibrational displacement is calculated using the reference and target images.Figures 3(c)3(d) show the profiles that were extracted from reference (Fig. 3(a)) and target images (Fig. 3(b)) in the lateral direction. The displacement of the sample can be calculated by analyzing variations in the extracted profiles.
I_{DC} andγ are affected by the nonuniformity of the light source and detector. To simplify the calculations, in this paper we assume thatI_{DC} andγ are equal in all of the images in monochromatic interferometry. We should also assume that the displacement due to vibration is smaller than a quarter wavelength, because of the ambiguity problem in monochromatic interferometry [20]. Figure 4 compares two extracted profiles;ε_{x} is the distance that the fringe pattern (profile) shifted, andl_{x} is the fringe spacing.As shown in Fig. 2, the ratio of the half wavelength
λ /2 tol_{x} represents the tangent of the tilt angleβ of the sample. Likewise, the ratio of displacement due to vibrationε_{z} toε_{x} is expressed by the tangent ofβ . Thereforeε_{z} can be determined from the following relationship:In this paper
l_{x} andε_{x} are calculated by sinusoidal fitting of the extracted profile.l_{x} is the peakto peak distance in a single extracted profile, whileε_{x} is the distance between peak positions when comparing different extracted profiles.Figure 5 shows a procedure for calculating the displacement caused by vibration. Images are captured, profiles are extracted from reference and target images, and finally the displacement of vibration is calculated from all extracted profiles using Eq. (5).
3.2. Phase Error Measurement Method
The phase error due to vibration can be calculated using the fringe spacing and the distance that the fringe pattern shifts. However, if the flat part of the sample is not wide enough, or the tilt angle is very small, the fringe spacing and the distance that the fringe pattern shifts cannot be obtained. Therefore, the following method is required to measure the phase error due to vibration.
The phase error can be measured by the subtraction of one profile from another profile. Figure 6 shows that
I_{Diff} is the result of this subtraction, wheree is the phase difference between profiles in the horizontal direction, which should not be confused with the phase errorθ_{error} from OPD change in thez direction.As in the section above, we should assume that the phase difference between two profiles is less than
π , to avoid ambiguity.Let the position at which
I_{Diff} is zero bex _{1}. Thenx _{1} satisfies the following relationships:The two cosine values cannot be the same because
e is not zero; therefore the two cosine values have the same magnitude but different signs. Hence, the following relationship is satisfied:Similarly, where
I_{Diff} has its largest value, positionx _{2} is expressed as in Eq. (8) whene is between zero andπ . Ife is between −π and zero,x _{2} isx _{1}−π /2 rather thanx _{1}+π /2. To simplify the calculation ofe , one needs to assume that the magnitude ofe is between zero andπ .Using Eq. (8), the magnitude of the phase difference 
e  can be determined.The calculated value of 
e  is between 0 andπ , as assumed above. To increase the range ofe to 2π , the sign ofe is required; by including the sign ofe , the range of the phase difference is expanded from [0,π ] to [−π ,π ]. The sign is determined from the positions of the peaks of the profiles, as shown below:In addition, as shown in Fig. 2,
e is proportional to the ratio of the distance that the interference fringe pattern has shifted,ε_{x} , to the fringe spacing,l_{x} .The phase error due to vibration is proportional to the ratio of displacement
ε_{z} to the half wavelengthλ /2. As mentioned in Sec. 3.1, the ratio ofλ /2 tol_{x} and the ratio ofε_{z} toε_{x} are equal to the tangent of the tilt angleβ . Thus the phase errorθ_{error} due to vibration has the same value as the phase differencee of the profiles.IV. EXPERIMENTAL RESULTS FOR SMALLAMPLITUDE VIBRATION
The Mirau interferometry employed in this study is shown in Fig. 1. We use an optical bandpath filter on an LED light source that provides quasimonochromatic light. The beam after the bandpass filter has a center wavelength of 550 nm and bandwidth of 20 nm. We use a multiaxis piezo nanopositioning stage P518.ZCL (PI) as the outofplane vibration generator. A 33500B series waveform generator (Agilent) is used to control the waveform and amplitude of the vibration. A Vieworks CCD camera with 640×480 pixels is used to capture the images. A 10 of the vibration. A Vieworks CCD camera with 640×480 pixels is used to capture the images. A 10x objective lens is used so that the field of view (FOV) of the experimental system is 464×348 μm.
There should be at least one whole fringe pattern in each image, to measure the vibrational amplitude and phase error. The lower limit of the tilt angle
β is determined from the FOV and the wavelength of the light. When there are two full fringe patterns in the horizontal direction of the FOV,β is about 0.14°. Therefore, in the proposed system the vibrational amplitude and phase error can be measured whenβ is greater than 0.14°.To measure the vibration, 100 images are captured without a mechanical moving process; only variations of the fringe pattern due to outofplane vibration are obtained. A voltage of 15 mV is applied to the vibration generator to create vibrations of amplitude 100 nm. We measure four waveforms using the proposed methods. Images are captured at time intervals of 1/35 s during vibration, and the displacement caused by the vibration is calculated using Eq. (5). Figure 7 shows that the measured displacement has the same shape and normalized amplitude as the input voltage. These results confirm that the proposed technique for displacement is a convenient and effective approach for measureming any waveform from small vibrations.
In the same manner, we measure phase error from the captured images using Eqs. (9) and (10). Figure 8 shows a comparison of phase error to input voltage. The results confirm that the proposed technique for phase error is also useful for any waveform.
V. ADVANCED VIBRATION MEASUREMENT METHODS
5.1. Ambiguity Problem of the Proposed Methods
There is an ambiguity of 2
π in monochromatic interferometry. Likewise, the proposed methods are also ambiguous. As shown in Fig. 9, the measured displacement is distorted when the vibrational amplitude is greater than a quarter wavelength. Ambiguities in the methods proposed above can be removed by applying advanced methods.5.2. Advanced Method for Largeamplitude Vibration: Relative Method
In the methods proposed earlier, reference and target images are used to calculate displacement and phase error, whereas in a relative method, two consecutive target images are used. As in the earlier methods, profiles are extracted from images, and the relative displacement and phase error are calculated using profiles from consecutive captured images. Displacement and phase error are determined by the accumulation of relative values, as shown in Fig. 10. By using a relative method, largeamplitude vibration can be measured without ambiguity, if one assumes that the relative displacement between two consecutive images is smaller than a quarter wavelength. This assumption is valid for cases in which the frame rate of the camera is higher than the vibrational frequency, as the frame rate corresponds to the reciprocal of the time interval required to acquire the images. In other words, largeamplitude vibrations can be measured unambiguously if the frame rate is higher than the vibrational frequency.
[FIG. 10.] Procedure for using relative displacement to calculate the displacement due to vibration.
If we let the relative displacement of the i^{th} and (i+1)^{th} image be (Δ
ε_{z} )_{i} , Eq. (5) can be rewritten as follows:where (Δ
ε_{x} )_{i} is the relative distance the extracted profiles shifted from the i^{th} to the (i+1)^{th} image. Wavelengthl and fringe spacingl_{x} do not change in this process. Thus, calculation of displacement from these relative value is given by:In the same way, a phase error larger than 2
π can be determined from the relative phase error Δθ_{error} . The i^{th} relative phase error, (Δθ_{error} )_{i}, from the i^{th} and (i+1)^{th} images is expressed as:Thus the phase error can be determined by:
A voltage of 150 mV is used to create a vibration of amplitude 1 μm. 100 images are acquired at time intervals of 1/55 s. Figure 11 shows the two waveforms that are measured by the relative method, and a comparison of the measured results to the input voltage.
As shown in Fig. 11, the relative method solves the ambiguity problem in largeamplitude vibration. The measured results have shape and normalized amplitude identical to those of the input voltage. These results confirm that the advanced methods are convenient and effective for measuring any vibration.
VI. FREQUENCY ANALYSIS AND COMPARISON TO RESULTS FROM AN ACCELEROMETER
Frequencies of vibration are as important as vibrational amplitude, and can be determined by frequency analysis. Measured displacement can be converted to vibrational frequency by frequency analysis such as fast Fourier transform, and the main frequencies of vibration can be determined by finding peaks in the frequency domain. When a camera is used to measure vibration, the range of measurable frequencies depends on the frame rate of the camera,
i.e . the number of images it can acquire in one second; The greater the frame rate, the wider the range of measurable frequencies.To verify the performance of the proposed method, measurement results are compared to those from an accelerometer, which is used predominantly in industry. Multifrequency vibration is measured by both techniques, and the results compared in the frequency domain. The measured vibration is the sum of four sine waveforms with frequencies of 1, 5, 10, and 20 Hz respectively. The amplitudes of all of the waveforms are the same, and the vibration is expressed as:
300 images are captured at intervals of 1/110 s, and displacement is calculated by the relative method. Figure 12 shows the displacement measurements.
The measured displacement is analyzed in the frequency domain by fast Fourier Transformation (FFT). Figure 13 shows the four peaks, which are identified at 1, 5, 10, and 20 Hz.
The same vibration is measured with an accelerometer (DS 2000 series, Onosokki) and the average result for 10 s evaluated using Onosokki’s software. Figure 14 shows three peaks identified by the accelerometer at 5, 10, and 20 Hz. These peaks are obvious, but the peak at 1 Hz looks like noise. The failure to measure the peak at 1 Hz is caused by the limitations of the accelerometer at low frequencies.
Using the proposed method, the same results can be obtained for frequencies greater than 5 Hz, and better results can be achieved for lower frrequencies, than with an accelerometer. Thus the proposed method is useful for measuring vibrational frequencies.
VII. CONCLUSION
In this paper we have proposed methods for measuring displacement and phase error due to vibration. We developed and tested the methods to verify their performance. This work is summarized as follows.
(1) Our system consisted of the following components: an LED light source, an optical bandpass filter to provide quasimonochromatic light, an objective lens, and a CCD camera to capture images. This is a basic Mirau interferometry setup that is used extensively to measure surface shapes. Moving parts or additional optical components are not needed. The proposed methods only use variations in fringe patterns caused by changes in the OPD. Thus, the proposed methods can be applied to measurement of environmental vibration as well as the vibration of the system itself. A vibration generator was used to create outofplane vibration. (2) Images were captured by a CCD camera at the same time interval, without any mechnical moving process. We referred to the first image that was obtained as the reference image, and the others as target images. The profiles extracted from the reference and target images in the lateral plane were used to calculate displacement and phase error due to vibration. We assumed that the displacement caused by the vibration was smaller than a quarter wavelength. (3) The vibrations we measured had amplitudes of 100 nm and four different waveforms: ramp, sine, square, and triangle. The calculated displacement and phase error were compared to theinput voltage, and the results confirmed that the proposed methods are useful for measuring the waveforms of small vibrations. (4) In the case of largeamplitude vibrations, the measured results were distorted. To solve this problem, we proposed the use of an advanced method called the relative approach. The advanced method used the same system as the earlier methods. Two consecutive images were captured and used to calculate relative displacement and phase error. The relative values were calculated from two consecutive profiles that were extracted in the same way as in the earlier methods. Then displacement and phase error were determined by the accumulation of relative values. We assumed that the relative displacement of two consecutive images was less than a quarter wavelength. (5) The vibrations we measured had amplitudes of 1 μm and two waveforms: sine and triangle. Relative values were calculated and accumulated to determine displacement and phase error. The measured results were compared to theinput voltage. This comparison indicated that the proposed methods are useful for measuring any type of vibration. (6) To verify the performance of the proposed methods, the measurement results were compared in the frequency domain with those from an accelerometer. Multifrequency vibrations of 1, 5, 10, and 20 Hzwere measured using the proposed methods and an accelerometer. The displacement measured by the proposed method was analyzed in the frequency domain by fast Fourier Transformation; then we measured the same vibration with an accelerometer. Both the proposed method and the accelerometer identified peaks at 5, 10, and 20 Hz, but the peak at 1 Hz was only identified by the proposed method. These results confirmed that the proposed method is useful for measuring vibrational frequencies. Furthermore, the proposed method performed better than a conventional accelerometer for lowfrequency vibration.
The proposed methods can be used to measure displacement and phase error due to outofplane vibration. Additionally, vibrational frequencies can be determined by an analysis of the displacement measured in the frequency domain. Therefore, vibration can be directly measured and easily managed using interferometry in industrial applications. The most important advantage of the proposed method is that it can be applied easily to almost any kind of interferometry, without any additional components.

[FIG. 1.] Mirau interferometry: BS, beam splitter; d1, distance from BS2 to the reference mirror; d2, distance from BS2 to the sample.

[]

[]

[]

[]

[FIG. 2.] Shift of a fringe pattern due to outofplane vibration in the case of a tilted sample: lx, period of the fringe pattern in the horizontal direction; λ, wavelength; εz, displacement due to vibration; εx, distance the fringe pattern shifts.

[FIG. 3.] Reference and target images and extracted profiles in the horizontal direction: (a) reference image, (b) target image, (c) profile extracted from the reference image, (d) profile extracted from the target image.

[FIG. 4.] Comparison of the two profiles: Ir, profile extracted from the reference image; It, profile extracted from the target image; lx, fringe spacing; εx, distance the fringe pattern shifted.

[]

[FIG. 5.] Procedure for calculating the displacement due to vibration.

[FIG. 6.] Comparison and subtraction of profiles from reference and target images.

[]

[]

[]

[]

[]

[]

[]

[FIG. 7.] Measured displacements of vibration compared to input voltage: (a) ramped wave input voltage and displacement caused by vibration; (b) sine wave; (c) square wave; (d) triangle wave.

[FIG. 8.] Measured phase errors, comparing vibration to input voltage: (a) ramp wave input voltage and phase error due to vibration; (b) sine wave; (c) square wave; (d) triangle wave.

[FIG. 9.] Ambiguity of measurementfor large vibration.

[FIG. 10.] Procedure for using relative displacement to calculate the displacement due to vibration.

[]

[]

[]

[]

[FIG. 11.] Measured displacement and phase error compared to input voltage: (a) and (b), displacement; (c) and (d), phase error.

[]

[FIG. 12.] Measured displacement of multifrequency vibration calculated by the proposed method.

[FIG. 13.] Peaks in the frequency domain measured by the proposed method.

[FIG. 14.] Peaks in the frequency domain measured by an accelerometer.