平面偏振光场光弹性实时相移法

光弹性法可广泛用于二维、三维应力或应变场的研究。根据偏振相移法的原理,提出一种新的积分相移法,并设计了用于全场等倾线和等差线实时分析的光弹性仪。在正交偏振光场下,按一定速度同步旋转起偏镜和分析镜,CCD摄像机可连续获取四幅相移光弹图像,由CCD像敏单元积分获得其光强值,通过所提出的相移法运算,可以迅速地得出全场等倾角和等差线的相位值。该方法和装置不仅可用于光学原件及玻璃制品的高速检测,而且也可用于随时间缓慢变化的光弹性分析中。     

Measurement of surface topography by RGB Shadow-Moiré with direct phase demodulation

In this paper we present the application of a direct demodulation method for the measurement of surface topography by means of Shadow-Moiré. In our set-up, we use three LEDs (with green, red and blue peak wavelengths) to illuminate the grating. Due to the different position of these light sources, a polychromatic Shadow-Moiré fringe pattern is produced, which can be described as the superposition of three monochromatic (red, green and blue) fringe patterns. Taking the image of this polychromatic fringe pattern with a RGB CCD camera, we get a monochromatic fringe pattern stored at each RGB channel of the CCD. The direct demodulation algorithm employed uses these fringe patterns to calculate the wrapped phase map. After unwrapping the phase map using a standard multi-grid technique, we implemented an automatic procedure to detect the area of interest of the phase map by removing low modulation zones and to calculate the absolute value of the phase. In this way it is possible to determine the topography of a surface with a single RGB snapshot maintaining a simple experimental set-up, which is an important feature, especially for the study of dynamic phenomena such as deformations. We present the experimental results obtained after measuring different objects with both smooth and rough surface textures.     

Digitally whole-field analysis of isoclinic parameter in photoelasticity by four-step color phase-shifting technique

This paper presents a whole-field method for automatic determination of the isoclinic parameter φ in photoelasticity. The method determines φ in the true phase interval (-π/2<φ+π/2) based on the four-step color phase-shifting approach and a new unwrapping method based on simple logic operations, a binary image and a gray-scale mask. A plane polariscope with a white light source is used for capturing raw photoelastic fringe images. The theoretical validation of the method is done with the problem of a circular disk under compression. Comparisons between theoretical and experimental results are also demonstrated. The experimental confirmation for the proposed method is tested with the photoelastic model having an isotropic point(s) and/or a singular point(s). The results show the method permits the reliable isoclinic parameter to be determined in the true phase interval with almost no defects of the isochromatic parameter.     

Whole field evaluation of stress components in digital photoelasticity—Issues, implementation and application

In photoelasticity, the method of obtaining the individual values of principal stresses/normal stresses separately is referred to as stress separation. Shear difference is one of the widely used techniques for stress separation in photoelasticity and one needs the value of fringe order and the isoclinic angle free of noise at every pixel over the domain. For accurate parameter determination, a ten-step phase shifting approach which uses a plane polariscope for isoclinic determination and a circular polariscope for isochromatic determination is proposed. A new quality guided approach for isoclinic unwrapping is developed. Isochromatic phasemap free of ambiguous zones is obtained by a new methodology and is unwrapped by a quality-guided approach. Whole field evaluation of stress components and its representation is then presented. The models used in this study are intentionally subjected to moderate loads showing a high level of isochromatic–isoclinic interaction. In view of this, the isoclinic data has several kinks which is found to cause streak formation in the whole field representation of separated stress components. An outlier smoothing algorithm is proposed for getting a smooth variation of the digital photoelastic parameters over the domain. Use of such smoothed data for stress separation has removed the streaks and has also greatly improved the accuracy of the separated stress components.     

Limitation of carrier fringe methods in digital photoelasticity

The carrier fringes method has been proposed in digital photoelasticity in combination with techniques such as Fourier transform and phase shifting method, without considering the influence of the isoclinics on the isochromatic patterns analysis. Unlike other optical methods as moiré and holographic interferometry, in photoelasticity the light intensity emerging from a circular polariscope is related to both the isochromatic retardation and the isoclinic parameter. As it is shown by the theoretical analysis, owing to the misalignment between the principal stresses in the model and in the carrier, the computed retardation is affected by an error which is the same for all photoelastic methods based on the use of carrier fringes. Consequently, the photoelastic analysis carried out by methods that use carrier fringes cannot be applied as a full-field technique. In detail, numerical simulations show that the retardation error is comparable (less than 0.05 fringe orders) with that of other photoelastic methods provided that the misalignment between the principal stresses in the model and in the carrier is less than 30°. On the contrary, in the model zones where the misalignment is higher than 30°, the retardation measurement can be affected by non negligible errors (up to 0.25 fringe orders).     

Pixel-by-pixel absolute phase retrieval using three phase-shifted fringe patterns without markers

This paper presents a method that can recover absolute phase pixel by pixel without embedding markers on three phase-shifted fringe patterns, acquiring additional images, or introducing additional hardware component(s). The proposed three-dimensional (3D) absolute shape measurement technique includes the following major steps: (1) segment the measured object into different regions using rough priori knowledge of surface geometry; (2) artificially create phase maps at different z planes using geometric constraints of structured light system; (3) unwrap the phase pixel by pixel for each region by properly referring to the artificially created phase map; and (4) merge unwrapped phases from all regions into a complete absolute phase map for 3D reconstruction. We demonstrate that conventional three-step phase-shifted fringe patterns can be used to create absolute phase map pixel by pixel even for large depth range objects. We have successfully implemented our proposed computational framework to achieve absolute 3D shape measurement at 40 Hz.     

Real-time phase-shift ing techniques for determining the photoelastic parameters: a theoretical comparison

Among data acquisition techniques in digital photoelasticity, the integrated phase shifting technique (IPST) can real-time analyze the photoelastic parameters at a video rate of the high speed CCD camera. In this paper, fourteen algorithms are described by different configurations of the rotating an analyzer at a constant rate and an output quarter-wave plate at another constant rate. The theoretical comparisons of the algorithms are given by the simulated phase distributions of the isochromatic and isoclinic parameters of the disk under two cases that the load keeps unchangeable or linearly increasing in exposure time of the camera. Then a guideline is given to alleviate the influence of the load changing with time on the IPST.     

Phase shifting photoelasticity in white light

The availability of image acquisition systems has led to the development of digital photoelasticity both in monochromatic and white light. In particular white light has been used mainly with the following methods: Spectral Content Analysis, RGB photoelasticity and phase shifting photoelasticity. Phase shifting photoelasticity in the colour domain has been used effectively for the determination of the isoclinic parameter in order to reduce the influence of the isochromatic fringes. The method has been also proposed for the determination of the isochromatic fringe order. This paper concerns the general characteristics of phase shifting photoelasticity in the colour domain. Special attention is drawn towards the influence of spectral content of the light source, the spectral response of the colour camera filters, the dispersion of the birefringence and the error of quarter wave plates.     

Numerical analysis of phase-stepping interferometric photoelasticity for plane stress separation

A phase-stepping method for interferometric photoelasticity is proposed to determinate whole-field plane stress components. Isopachic phase can be obtained by rotating polarizer, second quarter-wave plate and analyzer at definite optical arrangements. On the other hand, isochromatic and isoclinic phases can be determined in a circular polariscope arrangement. Furthermore, a load-stepping method is adopted to overcome the wrong mathematic sign of the isochromatics in the ambiguity regions and the influence of initial interferometric fringes on the isopachics. Light intensities and phase-stepping formula for the proposed method are derived using Jones calculus. Simulation of a circular disk under diametral compression demonstrates the feasibility of the proposed method.     

The influence of the quarter wave plates in automated photoelasticity

During the last decades, several methods have been proposed to automate photoelastic analyses. Some procedures are based on the circularly polarised light by using quarter wave plates. However, quarter wave plates are typically matched for a specific wavelength, and an error is introduced at different wavelengths. The error of quarter wave plates affects the measurement of isochromatic and isoclinic data. In this paper, the influence of the errors of quarter wave plates in some of the most common automated photoelastic methods is reviewed. The errors in the photoelastic data are given and the procedures to reduce, or eliminate, them are also suggested.     

Automatic evaluation of photoelastic fringe constant by the nonlinear least-squares method

A nonlinear least-squares method was applied to determine the photoelastic fringe constant from isochromatic phase field obtained by digital phase-shifting photoelasticity. To perform automatic evaluation, not only the photoelastic fringe constant but the center coordinates of a disk under diametral compression are also determined simultaneously using the nonlinear least-squares method. A set of over-determined nonlinear equations is obtained to estimate the unknown parameters by an iterative procedure based on Newton–Raphson method. Simulation and experimental results show that the parameters are evaluated accurately from the isochromatic phase field obtained from the phase-shifting photoelasticity. The photoelastic fringe constant and the disk center coordinates can be evaluated easily, simply and automatically by the technique.     

Automated photoelasticity: weighted least-squares determination of field stresses

The three-wavelength approach to phase-stepping photoelasticity as developed by the author is extended to determine automatically full-field stress tensor values. The only need for the user to calibrate the results is to give the material fringe value and the value of a stress at a single point. Four phase-stepped images illuminated by three wavelengths of light, that differ by prescribed increments, are collected using a semi-circular polariscope and an RGB CCD camera. A ramped phase map for the isochromatic parameter (α) is produced in the wrapped range −π/2<απ/2 that can be calibrated automatically. The value of the isoclinic angle (θ) is determined in the wrapped range −π/4<θπ/4. A discrete cosine transform algorithm has been developed to separate the stresses into Cartesian components. A convenience of the method is that accurate results can be obtained using a least-squares error minimisation process. The results obtained from experimental testing of a disc-in-compression specimen using transmission photoelasticity presented in comparison with theoretical solutions demonstrate the accuracy of the new approach.     

Towards effective phase unwrapping in digital photoelasticity

Phase-shifting techniques in photoelasticity have a peculiar problem of isoclinic–isochromatic interaction, which lead to ambiguous zones in phase map. A comprehensive discussion on the reason for the formation of ambiguous zones is presented in this paper. A new method is developed to remove this ambiguity in the phase map. The new methodology is validated with theoretically and experimentally generated phase maps. A flexible algorithm for phase unwrapping is developed using tiling. The applicability of this is demonstrated using various benchmark problems consisting of simply and multiply connected bodies and slices cut from stress frozen models.     

基于正弦和三角波条纹投影的三维测量方法

在现有的针对复杂物体表面形貌的三维测量方法中, 为了完成绝对相位的测量, 通常需要处理至少6幅条纹图像, 限制了测量速度。提出了采用2幅正弦条纹和2幅三角波条纹图来获得物体三维形貌的方法。利用两步相移正弦条纹和两步相移三角条纹得到截断相位, 再利用两步相移三角波条纹得到条纹级次, 减少了投影条纹幅数, 提高测量速度。在得到条纹级次时, 计算三角波条纹强度调制和强度对比度, 与计算相位相比, 可以减少数据处理的时间, 进一步提高测量速度, 同时能减小物体表面反射率的影响, 提高了测量精度。测量最大高度为39 mm的阶梯状标准块, 得到的最大绝对误差和最大的RMS误差分别为0.045 mm和0.041 mm。验证了该方法的有效性和实用性, 在高速实时的复杂形貌三维测量中有广泛的应用前景。     

Automatic interferogram analysis techniques applied to quasi-heterodyne holography and ESPI

This paper describes significant developments in methods for the automatic, quantitative analysis of interferograms. All areas of analysis have been considered: fringe field generation, pre-processing, and phase unwrapping.

A new quasi-heterodyne holographic technique is described in which the image is reconstructed using a single beam. The errors in the reconstructed fringe field are mainly linear in form, and an error compensation scheme is proposed. The final error in the phase measurement using automatic analysis is λ/40.

The process of image smoothing by an averaging filter is considered to reduce the effects of random noise. It is shown that by measuring the signal-to-noise ratio of the fringe field an optimum degree of smoothing may be applied. This is demonstrated on holographic and electronic speckle pattern interferometry (ESPI) data.

Two methods for cosinusoidal fringe image combination are compared, using three or four fields. It is shown that an automatic analysis can be achieved using four phase stepped images.

A new algorithm to automatically unwrap the phase of complex fringe patterns is described. The fringe field is segmented into small rectangular areas, called tiles. This allows local data to be obtained on fringe consistency and density. A confidence tree can then be formed to produce an optimal solution for the whole field. Results are presented and discussed for both holographic and ESPI data.     

相移阴影莫尔条纹正交化解调技术

提出一种基于克莱姆正则化分析法的三帧自标定相移阴影莫尔三维轮廓技术.该技术首先采用移动光栅的方法获得相移条纹图,然后通过不同帧相移条纹图相减去除条纹图背景,进而结合克莱姆正交化法和最小二乘法,发展了一种相位解调方法,提取了测量相位.以五步Harlharan算法为参考,用不同算法对同一物体表面进行测量.结果表明,相对于典型的三步相移法和主量分析方法,提出的方法测量得到的相位误差最小(0.5rad),且简化了测量过程.     

Generalized frequency selection in multifrequency interferometry

We present a generalized frequency selection method for N-frequency interferometry to form an optimum geometric series at synthetic wavelengths. The absolute range that is measurable is bounded by the number of beat frequency operations, phase noise, and the number of wavelengths used to form the geometric series of synthetic wavelengths. Theoretical predictions are compared with experimental results from a full-field fringe projector. A comparison of this technique with the method of excess fractions shows orders-of-magnitude faster processing with similar measurement reliability.     

Conversion from phase map to coordinate: Comparison among spatial carrier, Fourier transform, and phase shifting methods

For full-field phase measurement methods, many algorithms have been developed to extract a phase map from fringe image(s). Both phase wrapping and unwrapping algorithms have been extensively investigated by many researchers, but few papers can be found on how to calculate the coordinates of surface points from a phase map. This paper focuses on algorithms that show how a phase map can be used to calculate coordinates. Details are given for single image methods such as Fourier transform, spatial carrier methods, and multiple image methods like traditional phase-shifting methods. Algorithms that can be used to convert a phase map to coordinates and some issues related to these conversion algorithms are discussed. An artifact is measured using these phase measurement methods. The results show that using the correct algorithm to convert a phase map to coordinates is a key to obtaining accurate measurement results.     

Linear least squares approach for evaluating crack tip fracture parameters using isochromatic and isoclinic data from digital photoelasticity

In this work, digital photoelasticity technique is used to estimate the crack tip fracture parameters for different crack configurations. Conventionally, only isochromatic data surrounding the crack tip is used for SIF estimation, but with the advent of digital photoelasticity, pixel-wise availability of both isoclinic and isochromatic data could be exploited for SIF estimation in a novel way. A linear least square approach is proposed to estimate the mixed-mode crack tip fracture parameters by solving the multi-parameter stress field equation. The stress intensity factor (SIF) is extracted from those estimated fracture parameters. The isochromatic and isoclinic data around the crack tip is estimated using the ten‐step phase shifting technique. To get the unwrapped data, the adaptive quality guided phase unwrapping algorithm (AQGPU) has been used. The mixed mode fracture parameters, especially SIF are estimated for specimen configurations like single edge notch (SEN), center crack and straight crack ahead of inclusion using the proposed algorithm. The experimental SIF values estimated using the proposed method are compared with analytical/finite element analysis (FEA) results, and are found to be in good agreement.     

Fractional Fringe Order Method Using Fourier Analysis for Absolute Measurement of Block Gauge Thickness

A fringe analysis method, which combines the Fourier transform with the fractional fringe order method for measuring the absolute thickness of a block gauge, is presented. An approximate integer part of the fringe order is estimated by mechanical measurements, and the fractional part is determined by Fourier analysis of the interferometric fringe patterns. The fringe patterns are obtained with a Michelson interferometer by illumination of several selected wavelengths, respectively. The use of the fractional fringe order method can determine accurately more than 2ϖ phase jumps. The measured thickness of a block gauge is 5999115.0 ± 0.7 nm, which is comparable with the standard thickness of 6 mm. The influence of wavelength and approximate integer part of fringe order on the measurement accuracy is discussed.