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

High-performance optical coherence velocimeter: theory and applications

Open Access Open Access

Abstract

We proposed a high-performance optical coherence velocimeter (OCV) based on broadband optical interference which achieves spatial resolution from interference cancellation or enhancement of different components of the broadband light. There is a challengeable issue for OCV that the interference fringes become blurred when the velocity of detected object is relatively large, hindering the pace of OCV application in high-velocity field. To resolve this, the relationship between blurry coefficient and OCV system parameters (e.g., exposure time, central wavelength, bandwidth of source) was derived. It was found that blurry coefficient changed with oscillatory decay form and reached the minimum at each order blurry velocity. It showed that maximum measurable velocity of OCV systems could reach 10th order blurry velocity. The measurement of vibration of the loudspeaker driven by a function signal generator was employed to experimentally verify the velocity measurement performance of the system. The experiment demonstrated that the developed OCV can provide large velocity measurement ranges from static to 25.2 mm/s with nanometer-level precision and maximum measurable vibration frequency of up to 50 kHz. However, in theory, the theoretical maximum measurable velocity can be up to 1.06 m/s for current OCV configuration. The OCV has high precision, large dynamic range, and high-velocity measurement capability, making it attractive for applications in mechanical structure vibration monitoring and acoustic measurement.

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

1. Introduction

Mechanical vibration measurement has always been a focus of industrial production, structural health monitoring, MEMS (Micro-Electro-Mechanical system), rotating machinery, hearing and acoustics. The optical non-contact measuring device claimed particular advantages over traditional instrumentation, such as accelerometers or strain gauges etc., particular for measurements on hot, light, or rotating structures where traditional contacting instrumentation would change dynamics or be difficult to attach. Laser Doppler Vibrometry (LDV) is the most common optical vibration measurement technique [1,2]. It detects the Doppler frequency shift that occurs when the light is scattered by the moving surface which is directly proportional to the surface velocity. Usually, the frequency shift caused by the velocity of 0.1 m/s is almost 9 orders of magnitude lower than the frequency of laser [1]. However, this poses a high demand for the monochromatic nature of the laser source and the stability of the system. In our present work, we proposed Optical Coherence Velocimeter (OCV) which is a new optical non-contact velocity measuring technology based on the principle of broadband optical coherence which achieve spatial resolution from interference cancellation or enhancement of different components of the broadband light. Compared with the LDV, the OCV velocity measurement range and the working distance of the OCV are smaller, however the OCV has the advantages of low cost, small size, and high precision of velocity measurement, which makes the performance better for micro structure vibration measurement and acoustic measurement.

In fact, the OCV is developed from the Optical Coherence Tomography (OCT) technology. OCT is a powerful 3-D inner structure optical imaging technology proposed in the 90s of last century. It has been rapidly used in ophthalmology field [3] and other medical fields [4] because it can provide an inner observation of micron-structure of organ tissue. The development trend of OCT technology is higher imaging spatial resolution (~1μm), higher sensitivity (more than 90 dB) and faster imaging speed (real-time imaging), which mainly depends on the improvement of the system’s principle [5] and the application of better performance devices [6,7]. With the improvement of the system performance, the application fields of system are no longer limited in the field of medical imaging, such as industrial composites material performance testing [8], film coating thickness on-line monitoring [9], car paint monitoring [10], art heritage protection [11], and paper industry [12] and so on. In recent years, the OCT is also beginning to be applied in the field of displacement measurement. In our previous work, our research group enhanced the spatial resolution of OCT to nanoscale level by combining OCT with spectrum correction method [13] and PZT (Piezoelectric ceramic transducer)’s hysteresis curve was measured [14]. It shows great potential for displacement measurement of optical coherence technology. Even more, we have proposed an optical coherence vibration tomography (OCVT) technology [15,16] which can obtain both the inner structure and the vibration parameters of the sample simultaneously. But it is a pity that just some simple attempts and only low-speed (about 0.5 mm/s) and low-frequency (about 0~500 Hz) motion can be measured, limiting the further application of OCVT. It is noted here that the Doppler OCT [17,18] also a velocity measurement technology which can measure fluid flow in biological tissues. However, the technology is an approximate estimation of fluid flow and the maximum measurable velocity estimation limited by phase-shift measurements between axial scans, regardless of the method (color Doppler or velocity variance). In the present work, we propose the OCV technology which is based on the principle of the broadband optical coherence. The phenomenon that the interference fringes become blurred or disappear when measuring a high velocity moving object, was introduced and analyzed. Some solutions to improve the velocity measurement performance of the system, are given. The system developed in the work can provide a large velocity measurement range from static to 1.06 m/s with very high precision. The measured maximum vibration frequency can be 50 kHz. OCV technology has excellent velocity measurement performance and can enrich the optical non-contact velocity measurement method together with LDV technology.

2. Principle and set-up of optical coherence velocimeter (OCV)

The traditional OCT system is based on the broadband light interference principle and its basic structure is the Michelson interference structure. The back-scattered light from the target is mixed with a mutually coherent reference beam to produce interference. The result of the interference is that the intensity of the merged beam changes periodically in the spectral domain. This is so-called frequency domain optical coherence tomography (FD-OCT). The frequency of the interference fringes corresponds to the spatial position information of the target, so OCT provides a spatial resolution accuracy based on wavelength. The OCV system has the same structure as the typical OCT system, consisting of four parts: sample, interference structure, spectrometer, and signal analysis device (PC). There are two interference structures: free space structure and optical fiber structure. A simplest OCV system will include a 2*2 optical fiber as the interference structure, with a commercial spectrometer and signal processor (e.g. computer or DSP). The free-space system is slightly more complicated than the optical fiber one, but the bandwidth of the light is not limited by optical fiber, and the system has a better signal than the optical fiber one.

The schematic diagram of free-space and fiber-optic OCV system are shown in Fig. 1(a) and Fig. 1(b) respectively. The basic principle of these two systems are similarly the Michelson interference structure, therefore we take free-space OCV system as an example to explain the system in detail. The output of a super-luminescent diode (SLD, EBS300003-01, Exalos) with a central wavelength of 870 nm and bandwidth of 100 nm corresponding to a theoretical axial resolution of 3.3μm in air, was collimated by lens L1 (f = 25 mm). The source was split into a reference beam and a detection beam by a beam-splitter BS (50:50). In this experiment, the detected object was a silicon mirror that was driven by a loudspeaker which was controlled by the signal from a function signal generator. Both reflection beams of reference beam and detection beam were combined to a signal beam after passing through the beam-splitter BS again. In particular, all components of the OCV system (not including the PC and function signal generator) are fixed on an optical platform with vibration isolation to eliminate the interference of environmental vibrations. In the meanwhile, the OCV system is installed in a room at constant temperature and humidity condition to ensure the sub-nanometer displacement accuracy of the OCV system. A homemade spectrometer consisting of grating (2400 line/mm), lens (f = 50 mm) and linear CCD camera (P4-CM-04K10D-00-R, 100k fps, 4k*2 pixel, DALSA) was employed to detect the signal beam. The components of the spectrometer are encapsulated in an opaque box to block the interference of ambient light source. Interference occurs when both the reference beam and detection beam are spatially matched in orientation and their optical lengths are matched within the coherence length of the light source. Generally, the detection range of the system is determined by the spectral resolution of the spectrometer and the central wave-length of the spectrum, so the detection range of this system was 7.2 mm corresponding to the spectral resolution of 0.025 nm of the homemade spectrometer when the center wavelength is 870 nm. The traditional OCT acquires the spatial position information by doing the fast Fourier transform (FFT) to the spectral interferogram. However, due to signal leakage effects, the calculation of amplitude, phase and frequency by the FFT of the measured finite-length spectral interferogram is normally different from the real one. It means that the OCT system can only identify the integral change of the interferogram fringes, and the original accuracy of the spatial resolution will be several microns. A spectrum correction theory for sparse spectrum was used as an enhanced signal processing method for OCV system. It could improve the spatial resolution of the system from micron to nanometer level of 0.05 nm. The relevant principles and conclusion have been reported in the previous work [14,15]. OCV is the derivation of displacement with respect to time to obtain the velocity of the target. The electronic clock provides an extremely high time accuracy, so the OCV system possess a very high velocity measurement accuracy because of its sub-nanometer displacement measurement accuracy.

 figure: Fig. 1

Fig. 1 Schematic diagram of optical coherence velocimeter (OCV) with free space reference structure (a) and fiber-optic reference structure (b). The vibration structure (LS, loud speaker) was controlled by a function signal generator (FSG). SLD, Super luminescent diode; PC, personal computer; BS, Beam-splitter; L, Lens; M, Mirror; SM, silicon mirror; FO, 2*2 fiber optic; G, Grating; Detector, Linear array camera

Download Full Size | PDF

3. Analysis of interference fringes becoming blurred at high velocity

In the traditional applications, OCT technology is used for static or quasi-static object imaging. The spectral interferogram of the interferometer can be expressed as:

I(λ)=S(λ)[Ir+Is+n2IrInscos(22πλΔLn)+nm2ImsInscos(22πλ(ΔdnΔdm))]
where S(λ) is the spectral density function of light source;λ is the wavelength; Ir,Is,InsandIms are the intensity of the reference beam, the sample beam, the reflection light from n and m layer of the sample, respectively. ΔLn means the optical path difference between nth layer of the sample and the reference mirror;Δdn means the optical path length from nth layer within the object to the surface of object. The first two terms in this equation are the DC intensities returning from the reference and detection arms. The third term carries the wanted information such as surface position and the inner structure information of the sample. The last term, the auto-correlation term, is usually negligible because of the value of this term is too small comparing to another term.

In the OCV system, the detected object is no longer static but dynamic. For a vibrating object, the optical path difference ΔL is no longer constant. The optical path difference ΔL becomes a function of timet. Since the signal collected by the camera, the signal receiver, is the result of light intensity integral in the exposure time, the spectral interferogram of OCV system can be expressed as:

I(λ)=S(λ)t1t2(Ir+Is+2IrIscos(22πλ(ΔL(t)+ν(t)t)))dt

wheret1,t2 represent the beginning and the end time of the camera exposure. ν(t)means the velocity of the vibration object. In order to simplify the analysis process, we assume that the object only has the surface layer along the axial direction and does not consider the auto-correlation term. In the process of the object vibration, for a sample with multi-layered structure, the distance between layers is constant, so making the above assumption does not affect the subsequent analysis results.

The Eq. (2) can be rewritten as:

I(λ)=S(λ)*(t2t1)+S(λ)*2IrIs*λ4πν(t2+t12)*[sin(4πλ(ΔL(t1)+ν(t2+t12)*t2)sin(4πλ(ΔL(t1)+ν(t2+t12)*t1)]=S(λ)*(t2t1)[(Ir+Is)+2IrIs*λ(2πν(t2+t12)*(t2t1))*sin(2πν(t)*(t2t1)λ)*cos(4πλ(ΔL(t1)+ν(t2+t12)*(t2+t12)))]
We can find that the intensity of the signal is proportional to the camera’s exposure time. Comparing Eq. (1) with Eq. (3), the first two terms which represent the DC terms have the same expression. But there are two changes in signal term. First, the optical path length which represent the position of object is changed fromΔL to ΔL(t1)+ν(t2+t12)*((t2+t1)/2), which means that the instantaneous position measured by OCV system is the position of the intermediate moment of the camera’s exposure time. Second, the intensity of the signal item is no longer simply determined by the strength of beam returned by the reference arm and the sample arm, a coefficient item which seems like of the form sin(β)/β is added to determine the intensity of the signal item together (The parameter β=2πv(t)(t2t1)/λ). We call the coefficient item as blurry coefficient and denote it by the symbol C=blur|sin(β)/β|.

The blurry coefficient item sin(β)/βis a classical mathematical expression. The result of this function oscillates like a sine function, but the amplitude of the oscillation decreases with the increment of the variableβ . The variableβ is related to the wavelengthλ, exposure timet and the vibration velocityν. For a certain system, the wavelength of light source is constant, so the variableβ is proportional to the exposure time t and vibration velocityν. The relationship between blurry coefficient and the vibration velocity is shown in Fig. 2. The position and negative signs of velocity represent the direction of motion. When the velocity is close to 0, the value of the blurry coefficient is close to 1. At this situation, we can omit the blurry coefficient, when the interferogram signal expression is for static or quasi-static OCT system. When the vibration velocity of the object increases, the blurry coefficient begins to decay rapidly. At some velocity value, e.g. the parameter β=nπ, the blurry coefficient becomes zero, which means that the contrast of interference fringe is zero, and the system cannot obtain the desired information. Over that velocity, the blurry coefficient starts rising up to a peak and then falling to zero again. The peaks are getting smaller and smaller and finally become to zero. It should be noted here that the sign of sin(β)/βdoes not represent whether the interference fringes are enhanced or canceled; only the magnitude of sin(β)/β represents the attenuation coefficient of the interference. Therefore, the blurry coefficient defined as the absolute value of the sin(β)/β. The velocities when blurry coefficients become 0 are called blurry velocities: first blurry velocity, second blurry velocity and so on.

 figure: Fig. 2

Fig. 2 The relationship between blurry coefficient and vibration velocity when the center wavelength is 800nm.

Download Full Size | PDF

νblur=nλ2(t2t1),n=1,2,3...

When the velocity of moving object is close to the blurry velocities, the interference fringes obtained by system become blurred and the system loses measurement ability. For real vibration measurement application, high velocity measurement capability is an important parameter of the system. It means that the larger the blurry velocity, the better the performance of the OCV system. According to Eq. (4), we could know that increasing the wavelength or reducing the exposure time could increase the system’s blurry velocities. However, changing the wavelength of the light source has a limited increase to blurry velocities, so the most effective way to increase the blurry velocities is to reduce the camera’s exposure time. In general, the system with better performance needs more high-speed, more sensitive camera and stronger optical density light source. Figure 2, the dash line is the relationship between blurry coefficient and velocity when exposure time is 100μs, and the solid line corresponds to 40μs . It can be seen that the smaller exposure time corresponds to greater blurry velocity of the system which allows to measure the object vibrating at higher velocity.

In the above analysis, the system will lose its measurement ability when the velocity of moving object is close to the first blurry velocity. In the experiment, we found that even if the velocity of vibration object is first-order blurry velocity, the interference fringes obtained by the developed OCV system will still not be completely blurred. The reason for such phenomena is due to the polychromatism of the broadband light source we used. According to Eq. (4), the blurry velocity is related to the wavelength. Assuming that the velocity of the vibration object ν is the first-order blurry velocity of the wavelengthλ1, the velocity ν is not first-order blurry velocity for the other wavelengths (e.g.λ2)., Furthermore, the farther the wavelength λ2 is away from the wavelengthλ1, the farther the vibration velocity deviates from the blurry velocity of the wavelengthλ1. The relationship between blurry coefficient and vibration velocity for different wavelength is shown in Fig. 3. We could find that different wavelengths correspond to different blurry velocities; also, large wavelength interval corresponds to large blurry velocity difference.

 figure: Fig. 3

Fig. 3 The velocity-independent blurry coefficient of different wavelengths.

Download Full Size | PDF

If a broadband light source is employed in the OCV system, the interference fringe contrast will be obvious even if the interference fringes of a part wavelength range will be blurred due to the corresponding velocity being located at blurry velocity. Hence, there is always some part of the entire wavelength range that will have obvious interference fringes if a broadband light source is used. For a narrow-band source (for example λwidth<0.025λ0) OCV system, the maximum velocity-measuring capability is the first blurry velocity. However, the velocity-measuring capability of broadband system can pass through the blurry velocities. Therefore, increasing the bandwidth of the light source is an alternative effective way to increase the velocity-measuring capability of the system. In the broadband OCV system, the interference signal collected is the result of a combination of all the wavelength components. The wavelength-independent blurry coefficient can be obtained by integrating the blurry coefficient by wavelength. For broadband source system, the blurry coefficient Cblur can be calculated as:

Cblur=λ0λ0+S(λ)*|sin(β)/β|dλS(λ)=2ln2πλwidth*exp[(2ln2λλ0λwidth)2]
where S(λ) means the normalized power spectral density function of Gaussian light source, λ0is the central wavelength of source and λwidth is the FWHM (full width at half maximum amplitude) which means the effective bandwidth of source.

It is difficult to get the numerical solution of an integral formula like Eq. (5). But we still have a way to get an approximate solution of the equation. We can discuss it in two situations.

The first case, when the wavelength λ varies within the bandwidth and the value ofβ changes a little, the situation corresponds to low vibration velocity. The change of the value ofβwithin entire bandwidth range of light source Δβnπλwidth/λ0at each order blurry velocity, wheren is the order number of blurry velocity. When Δβ is small (for example Δβ<π/6), the below formula can be established.

sin(β)β=sin(nπλ0/λ)(nπλ0λ)=sin(nπ+nπ(λ0λ)/λ)(nπλ0/λ)λλ0λ0

The blurry coefficient Cblur of each wavelength of nth order blurry velocity can calculated as

Cblur1=λ0λ0+2ln2πλwidth*exp[(2ln2λλ0λwidth)2]*|λλ0λ0|dλ=λwidth2ln2πλ0λ0λ0+exp[(2ln2λλ0λwidth)2]d(2ln2λλ0λwidth)2=λwidth2ln2πλ0=0.3388*λwidthλ0

From the above result, we can infer that in order to increase the blurry coefficient value at each order blurry velocity the effective way is to increase the ratio of the bandwidth and the central wavelength of the light source. Another interesting conclusion is that the blurry coefficient for each order blurry velocity is equal according to Eq. (7). However, in reality, as the order increases, the blurry coefficient decreases. The reason is that Eq. (7) is based on Eq. (6) to calculate the blurry coefficient. As the order of blurry velocity increases, that is, the Δβ (the change of the value of β within entire bandwidth range of light source) becomes larger, Eq. (6) will be false. The left part of Eq. (6) will gradually become smaller than the right part. This will result in that the calculated result of blurry coefficient Cblur become smaller. After performing simulations using light sources of various bandwidths, we find that this reduction of blurry coefficient of each order blurry velocity is only related to the change of the value of β within entire bandwidth range of light source. When Δβ equal to π/2, the value of blurry coefficient Cblurdecreases to 87% of the value calculated using Eq. (7) whilst when Δβequal toπ, the value of Cblur decreases to 60% of the value calculated using Eq. (7).

The second case, when the wavelengthλ varies within the bandwidth and the value of β changes by more thanπ, this situation corresponds to a large vibration velocity. At this situation, the value of βis much bigger than Δβ. The value of sin(β) changes very large in whole wavelength range. But the value of β changes slowly compared to sin(β) and it can be replaced by the β value at central wavelengthλ0. The light source can be considered as a rectangular spectrum light source whose normalized power spectral density function is: S(λ)=1λwidth,λ0λwidth2λλ0+λwidth2 . The blurry coefficient Cblurcan be simply approximated as

Cblur2λ0λwidth2λ0+λwidth21λwidth*sin(2πν*Δt/λ)*λ02πνΔtdλ=1λwidth*λ02πνΔt*λ0λwidth2λ0+λwidth2sin(2πν*Δt/λ)dλλ02πνΔt2π=1β*2π

In this case, as the velocity of object increases, the value of blurry coefficient Cblur will become steady and the oscillation of the value will disappear quickly. As can be seen from Eq. (8), the value of blurry coefficient Cblur is proportional to the central wavelength of the light source and inversely proportion to the vibration velocity and integration time of the spectrometer.

The above two cases are connected when Δβ within entire bandwidth of light source equal toπ. At this time, β=π*λ0/λwidth. So,

Cblur2=1β*2π=2π2λwidthλ0=0.2026*λwidthλ0Cblur1=60%*0.03388*λwidthλ0=0.2033*λwidthλ0Cblur2Cblur1
Therefore, we provide a simple way to calculate the value of blurry coefficient Cblur at each order blurry velocity. Figure 4 is the simulation results of the velocity dependent blurry coefficient with different bandwidths light sources. The velocity dependent blurry coefficient oscillates with the cycle of blurry velocity at low velocity, and the oscillation gradually becomes smaller and closer to the steady-state curve as the velocity becomes larger. The narrower the bandwidth of the light source is, the larger the amplitude of the oscillation is, and the smaller the value of the blurry coefficient of each order blurry velocity is. The steady-state curve is related to the central wavelength of the light source but independent of the bandwidth. However, broadband source’s blurry coefficient becomes steady more quickly. For the light source with a bandwidth of 160nm, the blurry coefficient becomes steady after five oscillations, corresponding toΔβ within entire bandwidth range equal to π at fifth-order blurry velocity. For the light source with a bandwidth of 80nm, the number of oscillation is 10, and Δβ within entire bandwidth range equal to π at tenth-order blurry velocity. Therefore, the discussion of above two cases are correct and the simulation results are consistent with the theoretical analysis results.

 figure: Fig. 4

Fig. 4 The velocity dependent blurry coefficient of system with different bandwidths.

Download Full Size | PDF

Here, we discuss the low-cost OCV system design, including choosing the bandwidth of light source and the achievable measurement capability of the system. It is related to the SSNR (Static Signal to Noise Ratio) which means the SNR of interference signal when the vibrator remains stationary. In order to achieve low-cost OCV system, the criteria for OCV system design are that the SNR reduction caused by blurry coefficient must be smaller than the SSNR; the bandwidth used by the system is as small as possible; the achievable velocity measurement capability is as large as possible. Therefore, it is generally required that the blurry coefficient of first-order blurry velocity should be greater than twice of SSNR, that is,

Cblur=0.3388*λwidthλ02*1SSNR

In this way, even if the blurry coefficient is close to the steady state and decreases to 60% of the first-order blurry coefficient, the interference signal still exists. The system has a high velocity measurement capability. If the SSNR of the system is large enough, using a narrow-bandwidth light source can also achieve high velocity measurement capability. Therefore, smoothing and brightening the surface of the vibration object to improve the SSNR is also an effective way to improve the system’s velocity measurement capability.

In the real engineering applications, sometimes the SSNR is unknown in advance, we generally recommend that the bandwidth of the light source should be greater than 1/10 of the central wavelength. In this case, the value of β changesπ/10 at first-order blurry velocity in the entire bandwidth range, and the value of blurry coefficient is 0.03388 corresponding 30 dB reduction of interference intensity. Even at tenth-order blurry velocity, the value of blurry coefficient is 0.02 corresponding 34 dB reduction. It satisfies most of the engineering applications. Hence, we could define the tenth-order blurry velocity as the maximum velocity measurement capability of the OCV system.

νOCV=5*λ0t2t1

4. Experiment and discussion

We Use the home-made OCV system to measure the vibration of loudspeaker which was excited by the signal generated from a function signal generator (Agilent 33220). The spectrum correction theory was used to improve the spectral resolution of the OCV system. The signal processing method based on spectrum correction method can be referred to reference [14] and [19]. In our previous work in Reference [14], we have experimentally verified the measurement accuracy of the OCV system by measuring the dynamic curves of a PZT stack actuator (Thorlabs) compared with a commercial STIL MICROMESURE system. MICROMESURE is a modular measurement system dedicated to high resolution 3D microtopography and to shape and texture analysis. The axial resolution of the MICROMESURE system is 8 nm. The results measured by the two systems show that the performance OCV system is much higher than the one of MICROMESURE system. It has been demonstrated that the developed OCV system has sub-nanometer measurement accuracy of 0.05 nm [14].The Fig. 5 (a) is OCV measured results with and without spectrum correction method for the sinusoidal vibration with amplitude in micrometers. From the figure, we could find that the spectrum-corrected OCV had better resolution. Furthermore, Fig. 5 (b) shows the he OCV measured results with spectrum correction method for the vibration amplitudes in nanometers. Combined with spectrum correction method, the displacement resolution of OCV system has been greatly improved to nanometer level resolution.

 figure: Fig. 5

Fig. 5 (a) The comparison of OCV measured results with and without spectrum correction method for the vibration with amplitude in micrometers; (b) the OCV measured results with spectrum correction method for the vibration amplitude in nanometers. For clarity, the curves are offset vertically.

Download Full Size | PDF

We choose a sinusoidal vibration to verify the high-velocity measurement performance of the developed OCV system. Firstly, a sinusoidal signal with amplitude of 6 v and frequency of 25 Hz generated by function signal generator is used to excite the loudspeaker vibration. The diaphragm of the loudspeaker vibrated in sinusoidal form under the driven of signals. The advantage of choosing this form of motion is that the velocity of object changes smoothly and the velocity at each time can be simply calculated by the current time and the maximum velocity; also, no additional velocity measuring device is needed. Figure 6 (a) is the loudspeaker’s vibration measured by the developed OCV system. The dashed line is the displacement of the vibrating loudspeaker diaphragm. The amplitude of vibration is 318μm, and the vibration period is 0.04 s corresponding to the frequency of 25 Hz. The velocity of vibration could be obtained by taking the derivative of the displacement with respect to time. The solid line in Fig. 6 (a) is the velocity curve of the vibration. The velocity curve is in cosine form with a maximum velocity of 25.2 mm/s.

 figure: Fig. 6

Fig. 6 (a) The displacement and velocity of the sinusoidal vibration measured by the OCV system. (b) The change of blurry coefficient in a motion cycle under different exposure time conditions.

Download Full Size | PDF

Figure 6 (b) is the blurry coefficient of the sinusoidal motion measured by OCV system with different exposure time. In order to keep the consistency of the experimental conditions, we inserted an attenuator with different extinction ratio into light path and adjusted the exposure time of the camera to make the same intensity of light received. The extinction ratios of 50%, 30%, 7.5% and 3.5% corresponded to the exposure time of 11.5μs, 19.5μs, 80μsand 165μsrespectively. In the present work, the central wavelength of the light source is 850 nm. The bandwidth is 100 nm slightly larger than 1/10 of central wavelength. According to Eq. (4), the first order blurry velocity of OCV system at 11.5 μsand the exposure time is 36.9 mm/s, this is greater than the maximum velocity 25.2 mm/s of sinusoidal motion. Therefore, in the whole motion cycle, the blurry coefficient of the system is less than 10 dB, the interference fringes collected by system are still clear. When the exposure time is 19.5μs, the first-order blurry velocity of the system is 21.8 mm/s. Considering that the velocity varies according to cosine form with a maximum velocity of 25.2 mm/s, the first-order blurry velocity occurs at 0.0035 s, 0.0165 s, 0.0235 s, 0.0365 s for the four exposure times. The blurry coefficient was reduced to about −33dB at first order blurry velocity and increased rapidly when it exceeded the first-order blurry velocity. When the exposure time is 80μs, the first-order blurry velocity of system is 5.3 mm/s. The motion cycle crosses four order blurry velocity, the blurry coefficient correspondingly had four decreasing minimum values. When the exposure time is 165μs, the first-order blurry velocity is 2.58 mm/s. The maximum velocity of motion is close to the 10th-order blurry velocity.It can be seen from the figure that the coefficient does have nine minimum values. It is proved that the OCV system has the measuring ability at 10th-order blurry velocity. It is noted that as the exposure time increases, the number of sampling points of a motion cycle is reduced. The exposure time of 11.5μscorresponds to 7844 points. For the exposure time of 165μs, the sampling point is only 480 points. Therefore, when the exposure time is large, the value of blurry coefficient is smoothed because there are not enough sampling points to describe the rapid change of blurry coefficient.

In summary, the velocity measurement capability of the system is inversely proportional to the exposure time; the shorter the exposure time, the higher the velocity measurement performance; also, when the bandwidth is wide enough, the velocity measurement performance of the system can reach 10th order blurry velocity.

In order to give a better understanding of the high-velocity measurement performance of the homemade OCV system, high-frequency vibration measurement results are shown in Fig. 7. Figure 7 (a) is the OCV measurement result of the vibration excited by the periodic swept frequency signal which sweep from 10 Hz to 500 Hz in 0.5 s. From the OCV measurement results, it can be seen that the vibration structure is basically vibrated according to the frequency of the swept signal. However, it still had residual oscillation in the low-frequency part which is originated from the last swept period. Figure 7 (b) is the OCV measurement results of vibration excited at 5 kHz, 10 kHz and 20 kHz. The measured frequencies by the OCV system are 4.9998 kHz, 9.9999 kHz and 19.9997 kHz, and the corresponding amplitudes are 680 nm, 70 nm and 46 nm. It is noted here that there a low frequency component in the vibration signal measured. That is due to the resonance of the loud speaker. These results demonstrated the excellent performance of the OCV system in high-velocity measurement field.

 figure: Fig. 7

Fig. 7 (a) The OCV measurement result of vibration excited by swept frequency signal with a range of 10-500 Hz (b) The OCV measurement result of 5 kHz, 10 kHz and 20 kHz acoustic signal.

Download Full Size | PDF

When the OCV system is applied in vibration measurement, it is more concerned with the maximum measurable frequency that the system can achieve. The displacement function of the vibration can be simply assumed as s=A*sin(2π*f*t), where A means the amplitude of vibration and thef means the frequency of the vibration. The maximum velocity of the vibration is:νmax=A*2πf. At the same time, the sampling of the system should satisfy Nyquist sampling law. Therefore, the maximum measurable frequency of the system should satisfy the following formula:

{f<fsampling2A*2πf<νOCV=5*λ0Δt

In this work, the central wavelength of the OCV system is 850 nm, the minimum exposure time is 4μs, and the maximum sampling frequency is 100 kHz. It can be calculated that the maximum measurable velocity of the OCV system is 1.06 m/s in theory. Under this condition, the maximum measurable frequency of the system is related to the amplitude of the vibration. The relationship can be shown in Fig. 8. The solid curve is the relation between maximum measurable frequency and amplitude calculated by Eq. (12). From this curve, we can obtain the maximum measurable frequency of the system from the amplitude of vibration. In real engineering applications, we can’t accurately know the amplitude of vibration, but we can estimate the approximate order of vibration amplitude. A so-called conservative curve was proposed (shown as the dashed line in Fig. 8) to estimate the maximum measurable frequency of vibration if the approximate amplitude’s order is known. Take an example, if we know that the amplitude order of vibration is a few microns, we can assume that the amplitude of vibration is 10 microns to calculate the maximum measurable frequency. Based on these two curves, the OCV system offers three application prospects: The first application is acoustic measurement, the frequency range of sound is from 20 Hz to 20 kHz, the measurable range of the amplitude of sound vibration can cover from 1.0e-9 m to 1.0e-5 m. The second application is mechanical structure vibration measurement, the frequency of mechanical structure is usually from several hertz to several thousand hertz, the measurable range of the amplitude of mechanical structure vibration can cover from sub-micron to several millimeters. The third application is displacement measurement, which is used to measure the displacement of low-velocity objects. Furthermore, the measurement range of the OCV system is very wide, from nanometers to centimeters.

 figure: Fig. 8

Fig. 8 The relationship between maximum measurable frequency and amplitude of vibration

Download Full Size | PDF

5. Conclusions

The optical coherence velocimeter (OCV) is an optical non-contact vibration measurement technology which is based on the broadband optical interference. It can provide nanometer level displacement measurement accuracy and the corresponding velocity measurement accuracy. When the target is moving at high speed, the interference fringes will be blurred and the OCV will lose the velocity measurement ability. Theoretical analysis showed that the exposure time of the camera, the central wavelength and the bandwidth of light source are the main factors that affect the maximum velocity measurement ability of the OCV system. The shorter the exposure time and the larger the central wavelength, the greater the maximum measurable velocity of the OCV system. The blurry coefficient was proposed to describe the contrast reduction of the interference fringes changes with oscillatory decay form as the velocity of the target increasing. It was found that the blurry coefficient changed with oscillatory decay form and reached the minimum at each order blurry velocity. The bandwidth of the light source can increase the trough of each order blurry velocity. It is strongly recommended that the bandwidth of the light source exceeds the 1/10 of the central wavelength, so that the maximum measurable velocity of the system can reach 10th order blurry velocity. A loudspeaker driven by function signal generator was used to verify the above conclusion. Subsequently, we analyzed the relationship between vibration amplitude, vibration frequency and the maximum measurable velocity. The experiment demonstrated that the OCV system developed in the work can provide a large velocity measurement ranges from static to 25.2 mm/s with nanometer-level precision and the maximum measurable vibration frequency of up to 50 kHz. Due to the limit of our current experimental facility, we cannot provide larger velocity vibration for the validation of the maximum measurable velocity of our OCV system. However, in theory, the maximum measurable velocity of the developed OCV system can be up to 1.06 m/s for current configuration (e.g., the central wavelength of 850 nm, the minimum exposure time of 4μs, and the maximum sampling frequency of 100 kHz). The optical coherence velocimeter has high precision, large dynamic range, and high-velocity measurement capability, making it attractive for applications in the field of mechanical structure vibration monitoring and acoustic measurement.

Funding

National Natural Science Foundation of China (NSFC) (51675103), State Key Laboratory of Mechanical System and Vibration (MSV-2018-07), Shanghai Natural Science Fund (18ZR1414200).

Disclosures

The authors declare that there are no conflicts of interest related to this article.

References

1. S. J. Rothberg, M. S. Allen, P. Castellini, D. Di Maio, J. J. J. Dirckx, D. J. Ewins, B. J. Halkon, P. Muyshondt, N. Paone, T. Ryan, H. Steger, E. P. Tomasini, S. Vanlanduit, and J. F. Vignola, “An international review of laser Doppler vibrometry: Making light work of vibration measurement,” Opt. Lasers Eng. 99, 11–22 (2017). [CrossRef]  

2. E. Yoshikawa and T. Ushio, “Wind ranging and velocimetry with low peak power and long-duration modulated laser,” Opt. Express 25(8), 8845–8859 (2017). [CrossRef]   [PubMed]  

3. A. Rashno, D. D. Koozekanani, P. M. Drayna, B. Nazari, S. Sadri, H. Rabbani, and K. K. Parhi, “Fully automated segmentation of Fluid/Cyst regions in optical coherence tomography images with diabetic macular edema using neutrosophic sets and graph algorithms,” IEEE Trans. Biomed. Eng. 65(5), 989–1001 (2018). [PubMed]  

4. S. Chen, J. Yi, W. Liu, V. Backman, and H. F. Zhang, “Monte Carlo investigation of optical coherence tomography retinal oximetry,” IEEE Trans. Biomed. Eng. 62(9), 2308–2315 (2015). [CrossRef]   [PubMed]  

5. J. F. de Boer, B. Cense, B. H. Park, M. C. Pierce, G. J. Tearney, and B. E. Bouma, “Improved signal-to-noise ratio in spectral-domain compared with time-domain optical coherence tomography,” Opt. Lett. 28(21), 2067–2069 (2003). [CrossRef]   [PubMed]  

6. M. Choma, M. Sarunic, C. Yang, and J. Izatt, “Sensitivity advantage of swept source and Fourier domain optical coherence tomography,” Opt. Express 11(18), 2183–2189 (2003). [CrossRef]   [PubMed]  

7. L. An, P. Li, T. T. Shen, and R. Wang, “High speed spectral domain optical coherence tomography for retinal imaging at 500,000 A‑lines per second,” Biomed. Opt. Express 2(10), 2770–2783 (2011). [CrossRef]   [PubMed]  

8. J. T. Oh and S. W. Kim, “Polarization-sensitive optical coherence tomography for photoelasticity testing of glass/epoxy composites,” Opt. Express 11(14), 1669–1676 (2003). [CrossRef]   [PubMed]  

9. M. Juuti, H. Tuononen, T. Prykari, V. Kontturi, M. Kuosmanen, E. Alarousu, J. Ketolainen, R. Myllylä, and K.-E. Peiponen, “Optical and terahertz measurement techniques for flat-faced pharmaceutical tablets: a case study of gloss, surface roughness and bulk properties of starch acetate tablets,” Meas. Sci. Technol. 20(1), 150301 (2009). [CrossRef]   [PubMed]  

10. J. Zhang, B. M. Williams, S. Lawman, D. Atkinson, Z. Zhang, Y. Shen, and Y. Zheng, “Non-destructive analysis of flake properties in automotive paints with full-field optical coherence tomography and 3D segmentation,” Opt. Express 25(16), 18614–18628 (2017). [CrossRef]   [PubMed]  

11. G. Latour, J. P. Echard, B. Soulier, I. Emond, S. Vaiedelich, and M. Elias, “Structural and optical properties of wood and wood finishes studied using optical coherence tomography: application to an 18th century Italian violin,” Appl. Opt. 48(33), 6485–6491 (2009). [CrossRef]   [PubMed]  

12. T. Prykari, J. Czajkowski, E. Alarousu, and R. A. Myllyla, “Optical coherence tomography as an accurate inspection and quality evaluation technique in paper industry,” Opt. Rev. 17(3), 218–222 (2010). [CrossRef]  

13. K. Ding, M. Xie, and X. Zhang, “Phase difference correction method for phase and frequency in spectral analysis,” Mech. Syst. Signal Process. 14(5), 835–843 (2000). [CrossRef]  

14. S. Zhong and Q. Zhang, “Enhanced optical coherence vibration tomography for sub-nanoscale-displacement-resolution calibration of piezoelectric actuators,” Sens. Actuators A Phys. 233, 42–46 (2015). [CrossRef]  

15. J. Zhong, S. Zhong, Q. Zhang, and L. Yao, “Two-dimensional optical coherence tomography for real-time structural dynamical characterization,” Opt. Lasers Eng. 66, 74–79 (2015).

16. S. Zhong, H. Shen, and Y. Shen, “Real-time monitoring of structural vibration using spectral-domain optical coherence tomography,” Opt. Lasers Eng. 49(1), 127–131 (2011). [CrossRef]  

17. V. X. D. Yang, M. L. Gordon, B. Qi, J. Pekar, S. Lo, E. Seng-Yue, A. Mok, B. C. Wilson, and I. A. Vitkin, “High speed, wide velocity dynamic range Doppler optical coherence tomography (Part I): System design, signal processing, and performance,” Opt. Express 11(7), 794–809 (2003). [CrossRef]   [PubMed]  

18. Y. Wang, D. Huang, Y. Su, and X. S. Yao, “Two-dimensional phase unwrapping in Doppler Fourier domain optical coherence tomography,” Opt. Express 24(23), 26129–26145 (2016). [CrossRef]   [PubMed]  

19. K. Ding, D. Cao, and W. Li, “An approach to discrete spectrum correction based on energy centroid,” Key Eng. Mater. 321, 1270–1273 (2006).

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 (8)

Fig. 1
Fig. 1 Schematic diagram of optical coherence velocimeter (OCV) with free space reference structure (a) and fiber-optic reference structure (b). The vibration structure (LS, loud speaker) was controlled by a function signal generator (FSG). SLD, Super luminescent diode; PC, personal computer; BS, Beam-splitter; L, Lens; M, Mirror; SM, silicon mirror; FO, 2*2 fiber optic; G, Grating; Detector, Linear array camera
Fig. 2
Fig. 2 The relationship between blurry coefficient and vibration velocity when the center wavelength is 800nm.
Fig. 3
Fig. 3 The velocity-independent blurry coefficient of different wavelengths.
Fig. 4
Fig. 4 The velocity dependent blurry coefficient of system with different bandwidths.
Fig. 5
Fig. 5 (a) The comparison of OCV measured results with and without spectrum correction method for the vibration with amplitude in micrometers; (b) the OCV measured results with spectrum correction method for the vibration amplitude in nanometers. For clarity, the curves are offset vertically.
Fig. 6
Fig. 6 (a) The displacement and velocity of the sinusoidal vibration measured by the OCV system. (b) The change of blurry coefficient in a motion cycle under different exposure time conditions.
Fig. 7
Fig. 7 (a) The OCV measurement result of vibration excited by swept frequency signal with a range of 10-500 Hz (b) The OCV measurement result of 5 kHz, 10 kHz and 20 kHz acoustic signal.
Fig. 8
Fig. 8 The relationship between maximum measurable frequency and amplitude of vibration

Equations (12)

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

I (λ) = S (λ) [ I r + I s + n 2 I r I ns cos( 2 2π λ Δ L n ) + nm 2 I ms I ns cos( 2 2π λ (Δ d n Δ d m ) )]
I (λ) = S (λ) t 1 t 2 ( I r + I s +2 I r I s cos(2 2π λ (Δ L (t) + ν (t) t)))dt
I (λ) = S (λ) *( t 2 t 1 )+ S (λ) *2 I r I s * λ 4π ν ( t 2 + t 1 2 ) *[sin( 4π λ (Δ L ( t 1 ) + ν ( t 2 + t 1 2 ) * t 2 )sin( 4π λ (Δ L ( t 1 ) + ν ( t 2 + t 1 2 ) * t 1 )] = S (λ) *( t 2 t 1 )[( I r + I s )+2 I r I s * λ (2π ν ( t 2 + t 1 2 ) *( t 2 t 1 )) *sin( 2π ν (t) *( t 2 t 1 ) λ )*cos( 4π λ (Δ L ( t 1 ) + ν ( t 2 + t 1 2 ) *( t 2 + t 1 2 )))]
ν blur = nλ 2( t 2 t 1 ) , n=1,2,3...
C blur = λ 0 λ 0 + S (λ) *| sin(β)/β | dλ S (λ) = 2 ln2 π λ width *exp[ (2 ln2 λ λ 0 λ width ) 2 ]
sin(β) β = sin(nπ λ 0 /λ) ( nπ λ 0 λ ) = sin(nπ+ nπ( λ 0 λ)/λ ) ( nπ λ 0 /λ ) λ λ 0 λ 0
C blur 1 = λ 0 λ 0 + 2 ln2 π λ width *exp[ (2 ln2 λ λ 0 λ width ) 2 ] *| λ λ 0 λ 0 |dλ = λ width 2 ln2 π λ 0 λ 0 λ 0 + exp[ (2 ln2 λ λ 0 λ width ) 2 ] d (2 ln2 λ λ 0 λ width ) 2 = λ width 2 ln2 π λ 0 =0.3388* λ width λ 0
C blur 2 λ 0 λ width 2 λ 0 + λ width 2 1 λ width *sin(2πν*Δt/λ) * λ 0 2πνΔt dλ = 1 λ width * λ 0 2πνΔt * λ 0 λ width 2 λ 0 + λ width 2 sin(2πν*Δt/λ) dλ λ 0 2πνΔt 2 π = 1 β * 2 π
C blur 2 = 1 β * 2 π = 2 π 2 λ width λ 0 =0.2026* λ width λ 0 C blur 1 =60%*0.03388* λ width λ 0 =0.2033* λ width λ 0 C blur 2 C blur 1
C blur =0.3388* λ width λ 0 2* 1 SSNR
ν OCV =5* λ 0 t 2 t 1
{ f< f sampling 2 A*2πf< ν OCV =5* λ 0 Δt
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.