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

Elliptically polarized light photoelasticity based on LCD

Open Access Open Access

Abstract

We propose a three-wavelength elliptically polarized light photoelasticity method for high efficiency and low-cost stress measurement. By illuminating the sample with two different forms of elliptically polarized light for each wavelength sources, twelve images are acquired. From these images, phase delay and the principal internal stress difference are precisely computed using developed algorithms. Our proposed method based on an LCD panel has the unrivalled advantage that elliptically polarized light can be automatically adjusted, which reduces the mechanical rotation of the system, in contrast to the traditional six-step phase-shifting photoelasticity method, which requires manual rotation for circularly polarized light. In addition, the system has the potential to theoretically expand the area of illumination infinitely, thereby expanding the measurement area. The viability of the suggested methods is confirmed with numerical simulation and stress measurement.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Stress-induced birefringence is an important characteristic that is widely observed in optical materials such as silicon, plastic, and glass [13]. To investigate the form of stress, various methods have been developed based on effects of stress on materials and components [37]. A majority of these techniques are primarily based on stress-strain and photoelasticity measurements [8] and also incorporates a variety of other approaches such as x-ray diffraction [9,10], Stoney curvature [1113], and micro-Raman spectroscopy [1416]. The development of photoelasticity approaches, such as the striped glass slide [17], RGB photoelasticity [18], and phase-shifting [1924] procedures, has sparked renewed interest due to the accessibility of cutting-edge computational methods and advanced image processing techniques. Some methods for measuring large samples appear as the sample size rises, such as the deformation loading method [25], Layer-by-layer milling method [26], and large diameter linear polarization measuring instrument [27]. However, these large-sample measurement approaches result in spoiling the sample or low accuracy. Phase-shifting, as a non-destructive stress measurement method, has a straightforward measuring process and small relative error. A phase-shifting approach to measuring stress in large samples would be useful in a range of commercial applications.

Traditional 2D stress models consider two major stresses to be orthogonal on the measurement plane, with the difference between them representing the level of stress at a given location [28,29] as

$$\delta = \frac{{2\pi Cd({{\sigma_1} - {\sigma_2}} )}}{\lambda }$$

C denotes the stress optical coefficient of the sample, which is determined by the model material and the wavelength λ of the illuminating light. The thickness of the sample perpendicular to the direction of light propagation is represented by d, and $({{\sigma_1} - {\sigma_2}} )$ is the principal stress difference at the given point.

Phase shifting is a technique that involves introducing a specific phase shift by rotating a wave plate, which in turn modulates the phase information in the domain. Phase-shifting techniques, such as six-step phase-shifting [30], eight-step phase-shifting [31], and others, are currently mature technologies. In the majority of the existing phase-shifting photoelasticity investigations, a circularly polarized light is produced by two quarter-wave plates, which is then used to illuminate the sample and take six intensity images. The principle of classical six-step phase-shifting method is schematically shown in Fig. 1. The light system is symmetrical and circularly polarized light generated by a polarizer ${P_1}$ and a quarter-wave plate ${Q_1}$ illuminates the sample, and the linearly polarized light carrying stress information is captured by a CCD after passing through the quarter-wave plate ${Q_2}$ and polarizer ${P_2}$. The intensity at the detection plane is derived from the Jones matrix as

$$I = {I_b} + \frac{1}{2}{I_0}[{1 \pm sin 2({\beta - \gamma } )cos\delta \mp cos2({\beta - \gamma } )sin2({\theta - \gamma } )sin\delta } ]$$
where θ denotes the angle between the principal stress direction and the x-axis, γ represents the angle between fast axis of the second wave plate ${Q_2}$ and the x-axis, β corresponds to the angle between the vibration direction of the analyzer ${P_2}$ and the x-axis. The addition and subtraction symbols in the Eq. (2) correspond to left-handed circularly polarized light and right-handed circularly polarized light, respectively (the upper part represents left-handed light, and the lower half represents right-handed light). Ib denotes the light intensity of the background, and I0 represents the intensity of the illumination. Using the six combinations of $\beta $ and $\gamma $ as shown in Table 1, phase shift $\delta $ can be precisely calculated according to Eqs. (3) and (4).
$$\theta = 0.5 \times ta{n^{ - 1}}\left( {\frac{{{I_3} - {I_5}}}{{{I_6} - {I_4}}}} \right)$$
$$\delta = \frac{\lambda }{{2\pi }} \times ta{n^{ - 1}}\left[ {\frac{{({{I_3} - {I_5}} )sin2\theta + ({{I_6} - {I_4}} )cos2\theta }}{{{I_1} - {I_2}}}} \right]$$

 figure: Fig. 1.

Fig. 1. Traditional photoelasticity measurement method based on circularly polarized light. ${P_1}$, ${P_2}$: Polarizers; ${Q_1}$, ${Q_2}$: Quarter-wave plates; $s$: Sample.

Download Full Size | PDF

Tables Icon

Table 1. The principle expression of classical six-step phase-shifting

In practical applications, classical phase-shifting methods are required to be illuminated with circularly polarized light, which requires a quarter-wave plate in front of the sample. As a result, circularly polarized light is obtained by rotating the plate. And in order to measure large samples, the aperture of the plate will increase, resulting in a sharp increase in the cost of the system.

This study proposes a simple and adaptable optical system for a stress assessment using an elliptically polarized light photoelasticity approach. By using liquid crystal display (LCD) instead of quarter-wave plates, the system can automatically adjust the form of elliptically polarized light, greatly reducing costs and providing a theoretical basis for large sample measurement. The illumination comprises of three closer wavelengths and an LCD panel directs two forms of elliptically polarized light towards the sample in a vertical position. A CCD then records images of linearly polarized light that has been modulated by a quarter wave plate and a polarizer. The principal stress difference $({{\sigma_1} - {\sigma_2}} )$ is then precisely determined from a set of 12 phase shifted images using developed algorithms. The experimental results with a polycarbonate sample show that the photoelasticity method of elliptically polarized light is in good agreement with the stress calculation results of the traditional six-step phase-shifting photoelasticity method. Moreover, this approach demonstrates twice the efficiency. The photoelasticity method using elliptically polarized light delivers greater precision than typical procedures, with an error value of only 0.36 × 105, according to numerical simulation results.

2. Theory

The principle of the proposed elliptically polarized light photoelasticity is illustrated in Fig. 2. The system is illuminated by an LED light source that emits three wavelengths (${\lambda _1},\; {\lambda _2},\; \; and\; {\lambda _3}$), while a Liquid Crystal Display (LCD) panel modulates the linearly polarized light from polarizer ${P_1}$ to an elliptically polarized light. Q is a wide-band quarter-wave plate, which modulates elliptically polarized light back to linearly polarized light, with the fast-axis of Q parallel to the major axis of elliptically polarized light. Polarizer ${P_2}$ is used to construct bright and dark fields. Since the wavelengths of illuminations are closer, respective refractive indices are almost equal. The outer polarization film must be removed from the LCD in the system in order to guarantee that elliptically polarized light is modulated. Liquid Crystals (LCs) possess both the fluidity of a liquid and the anisotropy of crystalline materials. The ordered arrangement of liquid crystal molecules results in an intermediate state having qualities of both crystal and liquid. Particularly cholesteric liquid crystals have exceptional optical rotation, making them perfect for producing elliptically polarized light [3234]. These liquid crystal molecules collectively rotate when a voltage is applied, with the angle of rotation being affected by the voltage [35,36]. By altering the degree of rotation of the liquid crystal molecules, one can produce a variety of elliptically polarized light, each with a unique angle that places its major axis differently than linearly polarized light. The elliptically polarized light emerging from the LCD can be expressed as:

$${e_x} = asin \left( {wt + \frac{\pi }{2}} \right)$$
$${e_y} = bsin({wt} )$$
${e_x}$ and ${e_y}$ are components along major and minor axis of elliptically polarized light with a and b as respective amplitudes. w is the vibration frequency of the light wave and t is the time. The x-axis of the coordinate system is set to overlap with the major axis of elliptically polarized light. The sample is transparent and the stress induced birefringence modulates elliptically polarized light to carry the stress information which is expressed by
$$e_x^{\prime} = a\sin \left( {wt + \frac{\pi }{2}} \right)co{s^2}\theta + \; bsin({wt} )sin\theta cos\theta + a\sin \left( {wt + \frac{\pi }{2} - \delta } \right)si{n^2}\theta - bsin({wt - \delta } )cos\theta sin\theta $$
$$e_y^{\prime} = a\sin \left( {wt + \frac{\pi }{2}} \right)cos\theta sin\theta + \; bsin({wt} )si{n^2}\theta - a\sin \left( {wt + \frac{\pi }{2} - \delta } \right)sin\theta cos\theta + bsin({wt - \delta } )co{s^2}\theta $$

 figure: Fig. 2.

Fig. 2. Schematic of elliptically polarized light photoelasticity measurement. $LED({{\lambda_1}/{\lambda_2}/{\lambda_3}} )$: LEDs with ${\lambda _1} = 591nm$, ${\lambda _2} = 613nm$, ${\lambda _3} = 631nm$; ${L_1}$, ${L_2}$: Lens; LCD: 7-inch LCD panel with 1 piece of polarization film on the side facing ${L_1}$; $s$: Sample; Q: Wide-band quarter-wave plate; P: Polarizer; CCD: Charge coupled device.

Download Full Size | PDF

$\theta $ represents the angle between principal stress direction and the x-axis. The orientation of analyzer ${P_2}$ is responsible for a bright or dark field illumination on the CCD plane. In particular, when the transmittance direction of polarizer ${P_2}$ is parallel to the linearly polarized light, a bright field is produced, and when it is perpendicular to the linearly polarized light, a dark field is produced. The intensities of bright and dark field images captured by CCD are represented by

$${I_ \bot } = {I_b} + {\sin ^2}\left( {\frac{\delta }{2}} \right)[{Acos4\theta + B} ]$$
$${I_\parallel } = {I_b} + {a^2} + {b^2} - {\sin ^2}\left( {\frac{\delta }{2}} \right)[{Acos4\theta + B} ]$$
where $A = \frac{{2{a^2}{b^2}}}{{{a^2} + {b^2}}} - \frac{{({{\textrm{a}^2} + \; {b^2}} )}}{2}$, $B = \frac{{2{a^2}{b^2}}}{{{a^2} + {b^2}}} + \frac{{({{\textrm{a}^2} + \; {b^2}} )}}{2}$. The background light ${I_b}$ present in both bright and dark field images is removed by subtracting Eq. (9) from Eq. (10) for the same wavelength of illumination as
$$I = {I_\parallel } - {I_ \bot } = {a^2} + {b^2} - 2{\sin ^2}\left( {\frac{\delta }{2}} \right)[{Acos4\theta + B} ]$$

In general stress models [37], the principal stress difference $({{\sigma_1} - {\sigma_2}} )$ is often regarded as the most crucial parameter, as it represents the magnitude of stress. Since C, d, and $\lambda $ are all known in Eq. (1), the phase delay $\delta $ now becomes an important parameter to determine. By using light sources of three different wavelengths for lighting and measurement, the determined value of $\delta $ becomes more accurate and Table 2 lists the final intensity calculation based on Eq. (11) for each wavelength.

Tables Icon

Table 2. Calculated intensity of each wavelength with elliptically polarized light obtained by LC rotation angles 0° and 45°

LC rotation angle, which represents the deflection angle of the major axis of the elliptically polarized light passing through the liquid crystal, is determined by calibrating the system prior to the experiment.

Intensities ${I_1}$, ${I_3}$, and ${I_5}$ are obtained with LC rotation angle of 0° (major axis of elliptically polarized light is parallel to the transmittance axis of the polarizer) and ${I_2}$, ${I_4}$, and ${I_6}$ are obtained with LC rotation angle of 45° (major axis of elliptically polarized light makes an angle of 45° with the polarizer transmittance axis).

The intensity equations in Table 2 are combined to eliminate the term related to the principal stress direction θ and obtain the phase retardation term as $\cos \delta $. For example, $cos{\delta _1}$ is obtained by combining equations of ${I_1}$ and ${I_2}$ and applying double angle formula as

$$cos{\delta _1} = \frac{{{A_2}({{I_1} - a_1^2 - b_1^2} )+ {A_1}({{I_2} - a_2^2 - b_2^2} )}}{{{A_2}{B_1} + {B_2}{A_1}}} + 1$$

Similar steps are followed to obtain $cos{\delta _2}$ from ${I_3}$ and ${I_4}$ and $cos{\delta _3}$ from ${I_5}$ and ${I_6}$ which are further used to obtain accurate stress difference $({{\sigma_1} - {\sigma_2}} )$ [38]. Let ${\delta _1} = 2\pi Cd({{\sigma_1} - {\sigma_2}} )/{\lambda _1}$, ${\delta _2} = 2\pi Cd({{\sigma_1} - {\sigma_2}} )/{\lambda _2},$ λ1<λ2, and ${\delta _1}({x,y} )- {\delta _2}({x,y} )< 2\pi $. Wavelength dependence of $\delta $ contributes to a small difference in the calculated values of ${\delta _1}$, ${\delta _2}$ and ${\delta _3}$. If $cos({{\delta_1}} )\ge cos({{\delta_2}} )$, $co{s^{ - 1}}[{\cos ({{\delta_1}} )} ]$ is a decreasing function in the range of [0, π], that is, ${\delta _1} ={-} co{s^{ - 1}}[{\cos ({{\delta_1}} )} ]$; on the contrary, $co{s^{ - 1}}[{\cos ({{\delta_1}} )} ]$ is an increasing function in the range of [0, π], that is, ${\delta _1} = {\cos ^{ - 1}}[{cos({{\delta_1}} )} ]$. Similarly, stress difference can also be obtained from ${\delta _2}$ and ${\delta _3}$. The specific formula is shown below

$${cos^{ - 1}}[{cos ({{\delta_1}} )} ]= \left\{ {\begin{array}{{cc}} {\; \begin{array}{{cc}} {{cos^{ - 1}}[{cos ({{\delta_1}} )} ],}&{cos ({{\delta_1}} )< cos({{\delta_2}} )} \end{array}}\\ {\begin{array}{{cc}} { -{cos^{ - 1}}[{cos ({{\delta_1}} )} ],}&{cos ({{\delta_1}} )\ge cos({{\delta_2}} )} \end{array}} \end{array}} \right.$$
$${cos^{ - 1}}[{cos ({{\delta_2}} )} ]= \left\{ {\begin{array}{{cc}} {\; \begin{array}{{cc}} {{cos^{ - 1}}[{cos ({{\delta_2}} )} ],}&{cos ({{\delta_2}} )< cos({{\delta_3}} )} \end{array}}\\ {\begin{array}{{cc}} { -{cos^{ - 1}}[{cos ({{\delta_2}} )} ],}&{cos ({{\delta_2}} )\ge cos({{\delta_3}} )} \end{array}} \end{array}} \right.$$

Unwrapping the phase retardation and combining the results of above two calculations, the accurate stress difference $({{\sigma_1} - {\sigma_2}} )$ can be obtained by the following formula

$$({{\sigma_1} - {\sigma_2}} )= \left\{ {\begin{array}{{cc}} {\; \begin{array}{{cc}} {\frac{{{\delta_1}{\lambda_1}}}{{2\pi Cd}},}&{\frac{{{\delta_1}{\lambda_1}}}{{2\pi Cd}} > \frac{{{\delta_2}{\lambda_2}}}{{2\pi Cd}}} \end{array}}\\ {\; \begin{array}{{cc}} {\frac{{{\delta_2}{\lambda_2}}}{{2\pi Cd}},}&{\frac{{{\delta_1}{\lambda_1}}}{{2\pi Cd}} \le \frac{{{\delta_2}{\lambda_2}}}{{2\pi Cd}}} \end{array}} \end{array}} \right.$$

Flow-chart of elliptically polarized light photoelasticity measurement is schematically shown in Fig. 3. Two types of elliptically polarized light with three wavelengths are used to acquire a total of 12 images in bright and dark fields. Intensities in Table. 2 are calculated from Eq. (11), with the background light removed which is further used to obtain $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$. $\frac{{{\delta _1}{\lambda _1}}}{{2\pi Cd}}$ and $\frac{{{\delta _2}{\lambda _2}}}{{2\pi Cd}}$ are obtained after phase unwrapping. The principal stress difference $({{\sigma_1} - {\sigma_2}} )$ is determined by the larger value of $\frac{{{\delta _1}{\lambda _1}}}{{2\pi Cd}}$ and $\frac{{{\delta _2}{\lambda _2}}}{{2\pi Cd}}$.

 figure: Fig. 3.

Fig. 3. Flow-chart of elliptically polarized light photoelasticity measurement.

Download Full Size | PDF

3. Simulation

The proposed method is verified in a numerical simulation using a virtual polycarbonate sample in the setup in Fig. 2. The sample has a diameter of 25 mm and it is loaded at two opposite points on the same horizontal plane and the two principal stress components (${\sigma _1}$ and ${\sigma _2}$) and phase difference ($\delta $) are expressed using the following formula [39,40]

$$\left\{ {\begin{array}{{c}} {{\sigma_1} = \frac{{{\sigma_x} + {\sigma_y}}}{2} + \sqrt {{{\left( {\frac{{{\sigma_x} - {\sigma_y}}}{2}} \right)}^2} + \tau_{xy}^2} }\\ {{\sigma_2} = \frac{{{\sigma_x} + {\sigma_y}}}{2} - \sqrt {{{\left( {\frac{{{\sigma_x} - {\sigma_y}}}{2}} \right)}^2} + \tau_{xy}^2} }\\ {\delta = \frac{{2\pi Cd({{\sigma_1} - {\sigma_2}} )}}{\lambda }} \end{array}} \right.$$
where C is the stress-optical constant of the given material, d is the thickness of sample. ${\sigma _x}$ and ${\sigma _y}$ are two principal inner stress components along x and y-axis respectively, which are expressed by Eq. (16) [41,42]
$$\left\{ {\begin{array}{{c}} {{\sigma_x} ={-} \frac{{2P}}{{\pi d}}\left\{ {\frac{{({R - y} ){x^2}}}{{r_1^4}} + \frac{{({R + y} ){x^2}}}{{r_2^4}} - \frac{1}{{2R}}} \right\}}\\ {{\sigma_y} ={-} \frac{{2P}}{{\pi d}}\left\{ {\frac{{{{({R - y} )}^3}}}{{r_1^4}} + \frac{{{{({R + y} )}^3}}}{{r_2^4}} - \frac{1}{{2R}}} \right\}}\\ {{\tau_{xy}} ={-} \frac{{2P}}{{\pi d}}\left\{ {\frac{{{{({R + y} )}^2}x}}{{r_2^4}} - \frac{{{{({R - y} )}^2}x}}{{r_1^4}}} \right\}} \end{array}} \right.$$
$R$ is the radius of the round plate sample and P is the load pressure. Also, $r_1^2 = {x^2} + {({R - y} )^2}$, $r_2^2 = {x^2} + {({R + y} )^2}$.

Numerical simulation assumed the following values: ${\lambda _1} = 591nm$, ${\lambda _2} = 613nm$, ${\lambda _3} = 631nm$, $P ={-} 150N$, $d = 1.5mm$, $C ={-} 400 \times {10^{ - 13}}P{a^{ - 1}}$. Phase difference for three wavelengths $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$ are calculated from 12 captured images and are presented in Figs. 4(a1)-4(a3), which corresponds to Step 3 in Fig. 3. The two wrapped images calculated by Eq. (12) and Eq. (13) are shown in Figs. 4(b1) and 4(b2), which corresponds to Step 4 in Fig. 3. The reconstructed $({{\sigma_1} - {\sigma_2}} )$ measured by elliptically polarized light photoelasticity is shown in Fig. 4(d), which corresponds to Step 6 in Fig. 3.

Theoretical stress difference obtained with the parameters used in numerical simulation is shown in Fig. 4(e). Figures 4(c1)-4(c6) show images captured by classical six-step phase-shifting method according to Table 1, and the reconstructed $({{\sigma_1} - {\sigma_2}} )$ is shown in Fig. 4(f). Figure 4(g) and 4(h) records the difference between theoretical $({{\sigma_1} - {\sigma_2}} )$ and the results obtained in numerical simulation and classical method. The errors are of the order of 10−5 (0.36 × 10−5 for elliptically polarized light photoelasticity and 0.8 × 10−5 for classical six-step phase-shifting) which concludes that the elliptically polarized light system has better measurement accuracy than the classical six-step phase-shifting method. The grid distribution of elliptically polarized light output from the LCD panel as illumination will have an impact, even if the method has shown its ability to increase efficiency in a number of stress measurements. This is addressed in the experiment section.

 figure: Fig. 4.

Fig. 4. Numerical simulation of elliptically polarized light photoelasticity and classical six-step phase-shifting photoelasticity. (a1)-(a3) $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$; (b1)-(b2) wrapped ${\delta _1}$ and wrapped ${\delta _2}$; (c1)-(c6) intensity images of classical six-step phase-shifting photoelasticity; (d) reconstructed $({{\sigma_1} - {\sigma_2}} )$ by elliptically polarized light photoelasticity; (e) theoretical $({{\sigma_1} - {\sigma_2}} )$; (f) reconstructed $({{\sigma_1} - {\sigma_2}} )$ by classical six-step phase-shifting photoelasticity; (g) Difference of (d) and(e); (h) Difference of (f)-(e).

Download Full Size | PDF

4. Experiments

The proposed elliptically polarized light stress measurement system is shown in Fig. 5(a) along with classical six-step phase-shifting photoelasticity system in Fig. 5(b). The light source has three LEDs with wavelengths 591 nm, 613 nm, and 631 nm. The system also includes a CCD camera (Allied Vision, Germany), LCD panel (7 “3-in-1/HD, Xianglan, China) with polarization film removed and a computer to modulate elliptically polarized light. The polarizer and wide-band quarter-wave plate (Thorlabs, United States) are located on the optical axis of the system. A polycarbonate plate with a diameter of 25 mm and a thickness of 3 mm is used as a sample. Symmetrical point forces are applied at either end of the plate's horizontal diameter. The custom made three wavelength LED light sources are connected by a three-in-one optical fiber with a diameter of 1 mm. A specific controller is used to select the desired wavelength that is coupled to the fiber.

 figure: Fig. 5.

Fig. 5. Optical system for photoelasticity measurement. (a) elliptically polarized light photoelasticity and (b) classical six-step phase-shifting photoelasticity. ${\lambda _1} = 591nm$, ${\lambda _2} = 613nm$, ${\lambda _3} = 631nm$; The arrangement includes LCD (liquid crystal display), ${P_1}$, ${P_2}$, $P$(polarizers), ${Q_1}$, ${Q_2}$, $Q$(wideband quarter-wave plate), $s$(sample), ${L_1}$, ${L_2}$, $L$(lens), $FH$(fiber holder) and CCD camera.

Download Full Size | PDF

Stress-induced birefringence of the polycarbonate sample is measured by modulating different types of elliptically polarized light as outlined in Table 2, which corresponds to Step 2 in Fig. 3. To evaluate the accuracy of the elliptically polarized light photoelasticity system, the same sample is also measured using the classical six-step phase-shifting photoelasticity method. Figures 6(a1)-6(a3) show Step 3 where $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$ are calculated by elliptically polarized light photoelasticity. Wrapped ${\delta _1}$ and ${\delta _2}$ in Step 4 are shown in Fig. 6(b1) and 6(b2). Figures 6(c1)-6(c6) show intensity images obtained by classical six-step phase-shifting photoelasticity. Figure 6(d) and Fig. 6(e) are the reconstructed $({{\sigma_1} - {\sigma_2}} )$ via proposed elliptically polarized light photoelasticity and classical six-step phase-shifting photoelasticity, respectively. Both of them are evenly distributed and in good agreement with the numerical simulation results, as shown in Fig. 4 (d) and (f). Figure 6(f) shows the difference among reconstructed $({{\sigma_1} - {\sigma_2}} )$ obtained by either methods, where the maximum absolute difference is 2 Pa. Larger variations are concentrated on the edge of the sample, which may be caused by the aberration of the system or the grid formation of liquid crystal molecules.

 figure: Fig. 6.

Fig. 6. Experimental verification of elliptically polarized light photoelasticity. (a1) -(a3) $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$; (b1) -(b2) wrapped ${\delta _1}$ and wrapped ${\delta _2}$; (c1) –(c6) intensity images of classical six-step phase-shifting photoelasticity; (d) reconstructed $({{\sigma_1} - {\sigma_2}} )$ via elliptically polarized light photoelasticity; (e) reconstructed $({{\sigma_1} - {\sigma_2}} )$ via classical six-step phase-shifting photoelasticity; (f) difference of (d) and (e)

Download Full Size | PDF

The effectiveness of the suggested method for larger stress measurement is investigated, and the findings are displayed in Fig. 7. The difference obtained in reconstructed $({{\sigma_1} - {\sigma_2}} )$ by either method indicates that the proposed method also works well in larger stress measurement, where the maximum absolute difference is 3 Pa.

 figure: Fig. 7.

Fig. 7. Experimental verification of elliptically polarized light photoelasticity in higher pressure. (a1) -(a3) $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$; (b1) -(b2) wrapped ${\delta _1}$ and ${\delta _2}$; (c1) –(c6) intensity images of classical six-step phase-shifting photoelasticity; (d) reconstructed $({{\sigma_1} - {\sigma_2}} )$ via elliptically polarized light photoelasticity; (e) reconstructed $({{\sigma_1} - {\sigma_2}} )$ via classical six-step phase-shifting photoelasticity; (f) difference of (d) and (e)

Download Full Size | PDF

In common LCD panels, liquid crystals exist in small units, resulting in a grid-like distribution. The distance between the LCD panel and the sample, and the aberrations in the system form different morphologies of the grid, with a square grid being the most common form. The existence of such grids affects the transmittance of LCD panel and hinders the measurement of background light that needs to be measured before the sample is formally measured. Considering that the grid is a high frequency spatial information, an ideal low-pass filter is used to solve the above problem. Figures 8(a1)-8(a3) show $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$ without low-pass filtering, and Figs. 8(b1)-8(b3) show $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$ after low-pass filtering. Without a frequency domain filtering, LC grid reduces fringe sharpness and contrast, which results in an error in stress calculation, as shown in Fig. 8(c) and Fig. 8(d). The diameter of the elliptically polarized light is determined by the size of the LCD panel and hence the system can theoretically achieve large aperture illumination of the sample.

 figure: Fig. 8.

Fig. 8. The importance of frequency domain filtering. (a1) -(a3) $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$ without frequency domain filtering; (b1) –(b3) $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$ after frequency domain filtering; (c) Reconstructed $({{\sigma_1} - {\sigma_2}} )$ without frequency domain filtering; (d) Reconstructed $({{\sigma_1} - {\sigma_2}} )$ after frequency domain filtering.

Download Full Size | PDF

To demonstrate the suggested method's potential for industrial applications, a quantitative study is performed using a quarter-wave plate and a half-wave plate with a central wavelength of 633 nm. The internal phase delay of both plates was evenly distributed. Figure 9 (a) represents the measurement result of the half-wave plate, while Fig. 9 (b) illustrates the measurement result of the quarter-wave plate. Only 613 nm single wavelength was used to measure the above phase delay, and the mean values are 3.4221 rad and 1.6470 rad, respectively, while the theoretical results were 3.3618 rad and 1.6195 rad, respectively. The phase errors of the two plates were 1.8% and 1.7%, respectively. The photoelasticity approach using elliptically polarized light has the flexibility to replace the polarizer and CCD combo with a polarizing camera. Because the wave plate requires only two revolutions, the photoelasticity measurement efficiency is boosted by a factor of two when compared to the typical six-step phase-shifting photoelasticity.

 figure: Fig. 9.

Fig. 9. Internal phase delay measured for (a) half wave plate and (b) quarter wave plate using the proposed method.

Download Full Size | PDF

5. Conclusion

In this paper, an elliptically polarized light photoelasticity measurement system is proposed, which employs that elliptically polarized light for illumination generated by LCD panel. The principal stress difference $({{\sigma_1} - {\sigma_2}} )$ is calculated from twelve intensity images obtained by illumination with two types of elliptically polarized light with three wavelengths. This system has an efficient modulation mode and greatly reduce the system cost, which can effectively reduce the mechanical rotation, and make it suitable for many applications. Results obtained in numerical simulation and experiments demonstrate the accuracy in the measurements. With advanced electronic control modulation technology of liquid crystals, large scale stress measurements can find more practical applications in engineering and material science.

Funding

National Natural Science Foundation of China (61827816).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. A. Ajovalasit, G. Petrucci, and M. Scafidi, “A critical assessment of automatic photoelastic methods for the analysis of edge residual stresses in glass,” J. Strain Anal. Eng. Des. 49(5), 361–375 (2014). [CrossRef]  

2. M. Scafidi, G. Pitarresi, A. Toscano, G. Petrucci, S. Alessi, and A. Ajovalasit, “Review of photoelastic image analysis applied to structural birefringent materials: glass and polymers,” Opt. Eng. 54(8), 081206 (2015). [CrossRef]  

3. K. Ramesh and V. Ramakrishnan, “Digital photoelasticity of glass: A comprehensive review,” Opt. Lasers Eng. 87, 59–74 (2016). [CrossRef]  

4. F.W. Hecker and B. Morche, “Computer-Aided Measurement of Relative Retardations in Plane Photoelasticity,” Springer Netherlands (1986).

5. A. Ajovalasit, S. Barone, and G. Petrucci, “Towards Rgb Photoelasticity - Full-Field Automated Photoelasticity in White-Light,” Exp. Mech. 35(3), 193–200 (1995). [CrossRef]  

6. P. J. Withers and H. Bhadeshia, “Overview - Residual stress part 1 - Measurement techniques,” J. Mater. Sci. Technol. 17(4), 355–365 (2001). [CrossRef]  

7. Z. Lei, D. Yun, Y. Kang, and L. Shao, “A Review of Digital Photoelasticity,” J. Exp. Mech. 19(4), 393–402 (2004).

8. S. Xiao, B. Li, H. Cui, and J. Wang, “Sensitive measurement of stress birefringence of fused silica substrates with cavity ring-down technique,” Opt. Lett. 43(4), 843 (2018). [CrossRef]  

9. S. Labat, P. Gergaud, O. Thomas, B. Gilles, and A. Marty, “Interdependence of elastic strain and segregation in metallic multilayers: An x-ray diffraction study of (111) Au/Ni multilayers,” J. Appl. Phys. 87(3), 1172–1181 (2000). [CrossRef]  

10. L. N. Taylor, A. K. Brown, A. J. Pung, E. G. Johnson, and J. J. Talghader, “Continuous-wave laser damage of uniform and nanolaminate hafnia and titania optical coatings,” Opt. Lett. 38(21), 4292–4295 (2013). [CrossRef]  

11. G. C. A. M. Janssen, M. M. Abdalla, F. van Keulen, B. R. Pujada, and B. van Venrooy, “Celebrating the 100th anniversary of the Stoney equation for film stress: Developments from polycrystalline steel strips to single crystal silicon wafers,” Thin Solid Films 517(6), 1858–1867 (2009). [CrossRef]  

12. J. M. Pureza, M. M. Lacerda, A. L. De Oliveira, J. F. Fragalli, and R. A. S. Zanon, “Enhancing accuracy to Stoney equation,” Appl. Surf. Sci. 255(12), 6426–6428 (2009). [CrossRef]  

13. A. Fluri, D. Pergolesi, V. Roddatis, A. Wokaun, and T. Lippert, “In situ stress observation in oxide films and how tensile stress influences oxygen ion conduction,” Nat. Commun. 7(1), 10692 (2016). [CrossRef]  

14. B. P. Groth, S. M. Langan, R. A. Haber, and A. B. Mann, “Relating residual stresses to machining and finishing in silicon carbide,” Ceram. Int. 42(1), 799–807 (2016). [CrossRef]  

15. P. Safaie, A. Eshaghi, and S. R. Bakhshi, “Optical properties of oxygen doped diamond-like carbon thin films,” J. Alloys Compd. 672, 426–432 (2016). [CrossRef]  

16. Z. W. Zhao, F. Xue, Y. Z. Mao, N. Chen, and T. Y. Wang, “Effects of annealing on the residual stresses distribution and the structural properties of Si core fiber,” Opt. Fiber Technol. 41, 193–199 (2018). [CrossRef]  

17. W. Rupeng and F. P. Chiang, “The Image Photo-Carrier Theory And Its Application To The Determination Of Principal Stress Direction,” Proc. SPIE 814, 257–261 (1987). [CrossRef]  

18. B. N. Simon and K. Ramesh, “Colour adaptation in three fringe photoelasticity using a single image,” Exp. Technol. 35(5), 59–65 (2011). [CrossRef]  

19. E. A. Patterson and Z. F. Wang, “Towards full field automated photoelastic analysis of complex components,” Strain 27(2), 49–53 (1991). [CrossRef]  

20. G. Petrucci, “Full-field automatic evaluation of an isoclinic parameter in white light,” Exp. Mech. 37(4), 420–426 (1997). [CrossRef]  

21. E. A. Patterson and Z. F. Wang, “Simultaneous observation of phase-stepped images for automated photoelasticity,” J. Strain Anal. Eng. Des. 33(1), 1–15 (1998). [CrossRef]  

22. X. Ji, S. Wei, X. Zhang, and J. Chen, “Optical experiment testing for fracture toughness of directional plexiglass along different directions,” Proc. SPIE 7375, 73751M (2009). [CrossRef]  

23. X. Liu and S. Dai, “Unwrapping isochromatic by combining white light photoelasticity with phase shifting thechnology,” J. Mech. Str. 38(4), 711–715 (2016). [CrossRef]  

24. J. C. Brinez, A. R. Martinez, and J. W. Branch, “Computational hybrid phase shifting technique applied to digital photoelasticity,” Optik 157, 287–297 (2018). [CrossRef]  

25. M. Li, L. Xu, and D. Wang, Load method for thermal deformation load stress of overflow brick, Caihong Display Devices Co Ltd (2021).

26. K. Tolstikhin, “An approach to differentiation of non-smooth functions obtained during residual stress measurements by layer-removal method,” J. Eng. Math. 103(1), 87–95 (2017). [CrossRef]  

27. S. Han, P. Xiao, D. Wang, F. Lu, J. Wang, J. Qiu, X. Shi, and H. Wu, Large-view-field stress analyzer, China Building Materials Test & Certific (2016).

28. S. M. Crawford and H. Kolsky, “Stress Birefringence in Polyethylene,” Proc. Phys. Soc. B 64(2), 119–125 (1951). [CrossRef]  

29. T. C. Liu, J. P. Hu, L. L. Zhu, R. Y. Zhou, C. Zhang, L. Y. Gu, X. D. Sun, J. F. Yu, A. J. Zeng, and H. J. Huang, “Mechanical stress birefringence of optical plates,” Appl. Opt. 59(24), 7371–7375 (2020). [CrossRef]  

30. B. Tao, Y. Yuan, X. L. Zhou, and H. Chen, “An integrated solution for compression-molded glass lenses,” Optical Design and Testing Vii 10021 (2016).

31. Z. LiNa, J. I. XinHua, and C. JinLong, “Research on eight-step phase shifting technique in digital photoelasticity and its applications,” J. Mech. Str. 31(3), 378–381 (2009).

32. I. Shiyanovskaya, A. Khan, S. Green, G. Magyar, O. Pishnyak, D. Marhefka, and J. W. Doane, “Single-substrate encapsulated cholesteric LCDs: Coatable, drapable, and foldable,” J. Soc. Inf. Disp. 14(2), 181–186 (2006). [CrossRef]  

33. M. C. Yeh, S. H. Yang, and W. Lee, “Color tuning in thermo-sensitive chiral photonic liquid crystals based on the pseudo-dielectric heating effect,” J. Mol. Liq. 296, 112082 (2019). [CrossRef]  

34. J. D. Lin, Y. S. Zhang, J. Y. Lee, T. S. Mo, H. C. Yeh, and C. R. Lee, “Electrically Tunable Liquid-Crystal-Polymer Composite Laser with Symmetric Sandwich Structure,” Macromolecules 53(3), 913–921 (2020). [CrossRef]  

35. S. W. Oh, J. M. Baek, S. H. Kim, and T. H. Yoon, “Optical and electrical switching of cholesteric liquid crystals containing azo dye,” RSC Adv. 7(32), 19497–19501 (2017). [CrossRef]  

36. J. M. Baek, S. W. Oh, S. H. Kim, and T. H. Yoon, “Fabrication of an initially-focal-conic cholesteric liquid crystal cell without polymer stabilization,” Displays 52, 55–58 (2018). [CrossRef]  

37. E.E.G. Pericles and S. Theocaris, Matrix Theory of Photoelasticity. (Springer, 1979).

38. S. C. Wang, X. L. Gao, P. F. Zhu, X. L. He, Z. L. Jiang, Y. Kong, S. Y. Wang, and C. Liu, “Three-wavelength digital photoelasticity,” Appl. Opt. 62(4), 972–978 (2023). [CrossRef]  

39. C. A. Magalhaes and P. A. A. M. Jr, “New numerical methods for the photoelastic technique with high accuracy,” J. Appl. Phys. 112(8), 3654–3658 (2012). [CrossRef]  

40. B. Tao, Y. Liu, L. C. Huang, G. J. Chen, and B. J. Chen, “3D reconstruction based on photoelastic fringes,” Concurr. Comput. 34(1), (2022).

41. J. H. Ye, F. Q. Wu, and J. Z. Sun, “Estimation of the tensile elastic modulus using Brazilian disc by applying diametrically opposed concentrated loads,” Int. J. Rock Mech. Min. Sci. 46(3), 568–576 (2009). [CrossRef]  

42. A. M. Comanici and P. D. Barsanescu, “Modification of Mohr's criterion in order to consider the effect of the intermediate principal stress,” Int. J. Plast. 108, 40–54 (2018). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

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

Fig. 1.
Fig. 1. Traditional photoelasticity measurement method based on circularly polarized light. ${P_1}$, ${P_2}$: Polarizers; ${Q_1}$, ${Q_2}$: Quarter-wave plates; $s$: Sample.
Fig. 2.
Fig. 2. Schematic of elliptically polarized light photoelasticity measurement. $LED({{\lambda_1}/{\lambda_2}/{\lambda_3}} )$: LEDs with ${\lambda _1} = 591nm$, ${\lambda _2} = 613nm$, ${\lambda _3} = 631nm$; ${L_1}$, ${L_2}$: Lens; LCD: 7-inch LCD panel with 1 piece of polarization film on the side facing ${L_1}$; $s$: Sample; Q: Wide-band quarter-wave plate; P: Polarizer; CCD: Charge coupled device.
Fig. 3.
Fig. 3. Flow-chart of elliptically polarized light photoelasticity measurement.
Fig. 4.
Fig. 4. Numerical simulation of elliptically polarized light photoelasticity and classical six-step phase-shifting photoelasticity. (a1)-(a3) $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$; (b1)-(b2) wrapped ${\delta _1}$ and wrapped ${\delta _2}$; (c1)-(c6) intensity images of classical six-step phase-shifting photoelasticity; (d) reconstructed $({{\sigma_1} - {\sigma_2}} )$ by elliptically polarized light photoelasticity; (e) theoretical $({{\sigma_1} - {\sigma_2}} )$; (f) reconstructed $({{\sigma_1} - {\sigma_2}} )$ by classical six-step phase-shifting photoelasticity; (g) Difference of (d) and(e); (h) Difference of (f)-(e).
Fig. 5.
Fig. 5. Optical system for photoelasticity measurement. (a) elliptically polarized light photoelasticity and (b) classical six-step phase-shifting photoelasticity. ${\lambda _1} = 591nm$, ${\lambda _2} = 613nm$, ${\lambda _3} = 631nm$; The arrangement includes LCD (liquid crystal display), ${P_1}$, ${P_2}$, $P$(polarizers), ${Q_1}$, ${Q_2}$, $Q$(wideband quarter-wave plate), $s$(sample), ${L_1}$, ${L_2}$, $L$(lens), $FH$(fiber holder) and CCD camera.
Fig. 6.
Fig. 6. Experimental verification of elliptically polarized light photoelasticity. (a1) -(a3) $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$; (b1) -(b2) wrapped ${\delta _1}$ and wrapped ${\delta _2}$; (c1) –(c6) intensity images of classical six-step phase-shifting photoelasticity; (d) reconstructed $({{\sigma_1} - {\sigma_2}} )$ via elliptically polarized light photoelasticity; (e) reconstructed $({{\sigma_1} - {\sigma_2}} )$ via classical six-step phase-shifting photoelasticity; (f) difference of (d) and (e)
Fig. 7.
Fig. 7. Experimental verification of elliptically polarized light photoelasticity in higher pressure. (a1) -(a3) $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$; (b1) -(b2) wrapped ${\delta _1}$ and ${\delta _2}$; (c1) –(c6) intensity images of classical six-step phase-shifting photoelasticity; (d) reconstructed $({{\sigma_1} - {\sigma_2}} )$ via elliptically polarized light photoelasticity; (e) reconstructed $({{\sigma_1} - {\sigma_2}} )$ via classical six-step phase-shifting photoelasticity; (f) difference of (d) and (e)
Fig. 8.
Fig. 8. The importance of frequency domain filtering. (a1) -(a3) $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$ without frequency domain filtering; (b1) –(b3) $cos{\delta _1}$, $cos{\delta _2}$, and $cos{\delta _3}$ after frequency domain filtering; (c) Reconstructed $({{\sigma_1} - {\sigma_2}} )$ without frequency domain filtering; (d) Reconstructed $({{\sigma_1} - {\sigma_2}} )$ after frequency domain filtering.
Fig. 9.
Fig. 9. Internal phase delay measured for (a) half wave plate and (b) quarter wave plate using the proposed method.

Tables (2)

Tables Icon

Table 1. The principle expression of classical six-step phase-shifting

Tables Icon

Table 2. Calculated intensity of each wavelength with elliptically polarized light obtained by LC rotation angles 0° and 45°

Equations (17)

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

δ = 2 π C d ( σ 1 σ 2 ) λ
I = I b + 1 2 I 0 [ 1 ± s i n 2 ( β γ ) c o s δ c o s 2 ( β γ ) s i n 2 ( θ γ ) s i n δ ]
θ = 0.5 × t a n 1 ( I 3 I 5 I 6 I 4 )
δ = λ 2 π × t a n 1 [ ( I 3 I 5 ) s i n 2 θ + ( I 6 I 4 ) c o s 2 θ I 1 I 2 ]
e x = a s i n ( w t + π 2 )
e y = b s i n ( w t )
e x = a sin ( w t + π 2 ) c o s 2 θ + b s i n ( w t ) s i n θ c o s θ + a sin ( w t + π 2 δ ) s i n 2 θ b s i n ( w t δ ) c o s θ s i n θ
e y = a sin ( w t + π 2 ) c o s θ s i n θ + b s i n ( w t ) s i n 2 θ a sin ( w t + π 2 δ ) s i n θ c o s θ + b s i n ( w t δ ) c o s 2 θ
I = I b + sin 2 ( δ 2 ) [ A c o s 4 θ + B ]
I = I b + a 2 + b 2 sin 2 ( δ 2 ) [ A c o s 4 θ + B ]
I = I I = a 2 + b 2 2 sin 2 ( δ 2 ) [ A c o s 4 θ + B ]
c o s δ 1 = A 2 ( I 1 a 1 2 b 1 2 ) + A 1 ( I 2 a 2 2 b 2 2 ) A 2 B 1 + B 2 A 1 + 1
c o s 1 [ c o s ( δ 1 ) ] = { c o s 1 [ c o s ( δ 1 ) ] , c o s ( δ 1 ) < c o s ( δ 2 ) c o s 1 [ c o s ( δ 1 ) ] , c o s ( δ 1 ) c o s ( δ 2 )
c o s 1 [ c o s ( δ 2 ) ] = { c o s 1 [ c o s ( δ 2 ) ] , c o s ( δ 2 ) < c o s ( δ 3 ) c o s 1 [ c o s ( δ 2 ) ] , c o s ( δ 2 ) c o s ( δ 3 )
( σ 1 σ 2 ) = { δ 1 λ 1 2 π C d , δ 1 λ 1 2 π C d > δ 2 λ 2 2 π C d δ 2 λ 2 2 π C d , δ 1 λ 1 2 π C d δ 2 λ 2 2 π C d
{ σ 1 = σ x + σ y 2 + ( σ x σ y 2 ) 2 + τ x y 2 σ 2 = σ x + σ y 2 ( σ x σ y 2 ) 2 + τ x y 2 δ = 2 π C d ( σ 1 σ 2 ) λ
{ σ x = 2 P π d { ( R y ) x 2 r 1 4 + ( R + y ) x 2 r 2 4 1 2 R } σ y = 2 P π d { ( R y ) 3 r 1 4 + ( R + y ) 3 r 2 4 1 2 R } τ x y = 2 P π d { ( R + y ) 2 x r 2 4 ( R y ) 2 x r 1 4 }
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.