Detection of the crack in an object is a critical problem for the health monitoring of a transparent object. The real-time and quantitative measurement of the crack-tip stress intensity factor (SIF) remains an open issue. In this paper, an approach for real-time and quantitative measurement for the SIFs of a Mode I crack is presented based on digital holographic interferometry (DHI). A transmission digital holographic system is established to measure the phase difference of an object wave during loading. The expression to achieve the SIF from the phase difference is formulated. To enhance the accuracy of measurement, calibrated phase unwrapping based on least-squares and iteration and median filtering is applied to retrieve the actual phase from the noisy wrapped one. The SIFs of the Mode I crack in a transparent polymethyl methacrylate (PMMA) specimen are measured by this approach. The results are compared with the theoretical ones to demonstrate the feasibility of the proposed approach.
National Natural Science Foundation of China11462009113620071. Introduction
Transparent materials are widely used in many engineering applications such as conductive materials [1–3], electrode materials [4–7], thermoelectric materials [8], encapsulation materials [9], and bulletproof materials [10, 11]. The fracture and fatigue problems are critical for many applications of the transparent materials. The crack is the key reason of materials failure. Usually the stress intensity factor (SIF) is the main quantity to evaluate the influence of the crack in the materials in elasticity. In the last decades, many measurement approaches were used to detect the SIF near the crack in transparent materials including photoelasticity [12, 13], coherent gradient sensing (CGS) [14–16], and digital coherent sensing (DGS) [17–19]. However, these approaches need birefringent materials, identification of the fringes, phase shifting, or fabrication of speckles which affect the application and accuracy of the measurement of the SIF.
As a noninvasive and full-field measurement method, digital holographic interferometry (DHI) has been widely used in many measurements including deformation [20–22], contour [23, 24], the field of fluid [25, 26], vibration [27, 28], and microscopy [29, 30]. By applying DHI in the measurement of transparent materials, one can get the stress-optic constant, stress field, and strain field [31, 32]. Especially, DHI can detect the damage and cracks in transparent materials which means it can be used in the health monitoring of a transparent structure and experimental study on the fracture and damage mechanics [33]. Usually, DHI obtains physical quantities through the phase difference of an object wave. As it is obtained through numerical operations based on the arctangent function, the phase is wrapped into the interval −π,π and needed to be unwrapped to retrieve the actual nonwrapped phase [34]. However, because stress near the crack tip is concentrated, the phase fringes are very close which results high speckle noise in the phase maps. The existence of high noise makes phase unwrapping very difficult and even to fail. Spatial filtering can be applied so as to reduce noise before phase unwrapping [35]. However, filtering will smear the phase jumps which will make unwrapping to fail. So, detection of the crack by DHI remains a qualitative analysis and does not realize the quantitative measurement which limits the application of DHI in the measurement of transparent materials.
In this paper, we aim to realize unwrapping the phase near the crack tip by using the calibrated phase unwrapping method based on least-squares and iterations [36]. Moreover, we present the scheme for the measurement of SIF near the crack tip in a transparent plate by DHI. In addition, the measured results are compared with the theoretical ones to validate the feasibility of the proposed approach. This paper is organized as follows: Section 2 first presents the principles of the proposed approach; Section 3 introduces the experimental scheme of the measurement for KI in transparent materials by the proposed approach; in Section 4, the measured results are given and validated by the comparison to the theoretical ones; and some conclusions are given in Section 5 finally.
2. Principles of the Proposed Approach2.1. Digital Holographic Interferometry
The principle diagram of transmission digital holographic interferometry [31] is shown in Figure 1. The input laser is divided into the object wave and reference wave by a beam splitter. The object wave and reference wave are both expanded and collimated by a spatial filter and a lens, respectively. The transparent object is placed in the collimated object wave and is imaged to the CCD target plane by an equivalent 4f system. Here, the negative lens in the system is used to adjust the field of view of interest. The collimated reference wave is refracted by a beam splitter and interfered with the object wave at the CCD plane. The principle axis of the reference wave forms a certain specified angle with that of the object wave which is normal to the CCD plane so that this system belongs to the image-plane off-axis holography. The CCD camera captures holograms at different loading states. The phase difference related to the stress in the specimen is reconstructed according to the following principles shown in Figure 1.
Principle diagram of the digital holographic system. BS1 and BS2 are the beam splitters; M1 and M2 are the mirrors; L1, L2, L3, and L4 are the positive lenses; CL is the negative lens; SF1 and SF2 are the spatial filters.
Establishing x and y coordinates in the CCD plane and z coordinate normal to the plane, the object wave Ox,y and the reference wave Rx,y can be expressed as(1)Ox,y=ox,yexpjkφx,y,Rx,y=rx,yexpjkθxx+θyy,where j=−1, k=2π/λ, λ is the wavelength of the laser, θx and θy are the angles between z axis and the projection of the reference wave vector on the xz plane and yz plane, respectively.
The hologram captured by the CCD is described as(2)Ix,y=o2x,y+r2x,y+ox,yexpjkφx,yrx,yexp−jkθxx+θyy+ox,yexp−jkφx,yrx,yexpjkθxx+θyy.
Performing Fourier transform on both sides of the above equation, we can get the spectrum of the hologram:(3)ℑIx,y=G0fx,fy+Gfx+θxλ,fy+θyλ+G∗fx−θxλ,fy−θyλ.
On the right side of this equation, the first item is the spectrum of the zero-order diffracted wave, the second is the spectrum of the object wave which is modulated by the reference wave and translated to the point (−θx/λ, −θy/λ) relative to the center of the spectrum plane, and the third is the spectrum of the conjugate object wave. So, by selecting the appropriate angles of θx and θy, we can separate the spectrum of the object wave Gfx,fy from the other spectrums. A rectangular window filter is used to extract Gfx,fy. Then, the inverse Fourier transform of Gfx,fy is carried out to reconstruct the complex modulated object wave fx,y as follows:(4)fx,y=ℑ−1Gfx,fy.
The phase difference of object waves at two loading states, due to the specimen’s deformation, is expressed as(5)ψx,y=anglef1x,yf2x,y,where ψx,y is the wrapped phase difference, f1x,y and f2x,y are the complex object waves at different loading states, and angle. is the function to calculate the phase angle of a complex value.
2.2. Phase Recovering Approach
In this paper, we use calibrated phase unwrapping algorithm based on the least-squares and iterations (CPULSI) [36] and median filtering to recover the true phase from the noisy wrapped phase. Here, φij represents the true phase and ψij (∈−π,π, i.e., modulo 2π) represents the wrapped phase at the grid point (i,j) of a phase map. Their relationship can be expressed by the wrapping operator W. which is defined as(6)Wφij=ψij=φij+2πkiji=0,1,…,M−1;j=0,1,…,N−1,where −π≤ψij≤π, kij is an integer, M and N are, respectively, the number of grid points with respect to the i and j indexes. The 1st order spatial wrapped phase derivatives are defined as(7)Δijx=Wψi+1j−ψiji=0,1,…,M−2;j=0,1,…,N−1,Δijx=0,otherwise,Δijy=Wψij+1−ψiji=0,1,…,M−1;j=0,1,…,N−2,Δijy=0,otherwise,where Δijx and Δijy are, respectively, the difference with respect to the i and j indexes.
In practice, the presence of noise will generate errors between noisy and noise-free phase derivatives and make unwrapping difficult, even to fail. Here, we proposed a calibration approach in [36] to calibrate the phase derivatives exhibiting large errors:(8)Δijx=sgnΔijxGxifΔijx≥Tx,Δijx=Δijx,otherwise,Δijy=sgnΔijyGyifΔijy≥Ty,Δijy=Δijy,otherwise,where sgn… is the signum function, Tx and Ty are the thresholds, and Gx and Gy are the calibrated phase derivatives. These parameters are defined as (E… means the statistical average)(9)Tx=EΔijx2−EΔijx2,Ty=EΔijy2−EΔijy2,Gx=1MN∑i=0i=M−1∑j=0j=N−1Δijx,Gy=1MN∑i=0i=M−1∑j=0j=N−1Δijy.
This calibration approach means that the phase derivatives whose values are larger than the standard deviation of phase derivatives are replaced by the average value of phase derivatives.
In the least-squares sense, the optimal solution φij can be obtained from the discrete Poisson equation with the Neumann boundary conditions:(10)φi+1j−2φij+φi−1j+φij+1−2φij+φij−1=ρij,where(11)ρij=Δijx−Δi−1jx+Δijy−Δij−1y.
The discrete Poisson equation can be solved by many methods such as fast Fourier transform (FFT), discrete cosine transform (DCT), or the multigrid method. In this algorithm, the DCT method is selected to solve the least-squares phase unwrapping problem.
Theoretically, the solution obtained from (10) is the exact one. However, there exist errors between the unwrapped phase and the true phase due to noise and to the smoothing performance of the least-squares method. In the proposed approach, iterations of unwrapped phase errors are utilized to seek more accurate results [36].
After unwrapping, median filtering is used to reduce the noise in the phase map. Median filtering is one of the most commonly used spatial filtering methods. It uses the median value in the kernel to substitute the value of center of the kernel [37]. So, it is easy to implement and efficient for impulsive noise and speckle noise. In addition, the phase map must be padded out of the frontiers with the length of the kernels symmetrically.
2.3. Calculation of the Stress Intensity Factor
For the uniform transparent plate specimen normal to the object wave principle axis, the relationship between the stress and phase difference of the object wave in plane stress state can be expressed as [31](12)σx+σy=λ2πhcΔφ,where σx and σy are the normal stresses in the plane, h is the thickness of the plate, c is the stress-optic constant of materials, and Δφ is the unwrapped denoised phase.
Taking into account the specimen with a Mode I through crack as shown in Figure 2, the stress near the crack tip can be expressed as [38](13)σx=KI2πrcosθ21−sinθ2sin3θ2,(14)σy=KI2πrcosθ21+sinθ2sin3θ2,where KI is the stress intensity factor (SIF) of the Mode I crack and r and θ are the polar coordinates relative to the crack tip.
Specimen shape and loading way.
By adding (13) and (14), we can get the equation to calculate KI from the stresses as(15)KI=σx+σy2πr2cosθ/2.
So, from (12) and (15), one can calculate the stress sum and SIF near the crack tip.
3. Experiment
In order to validate its feasibility, the proposed approach is used to measure the SIF near the crack tip in a transparent specimen. The specimen is fabricated with polymethyl methacrylate (PMMA). The shape and load of the specimen are shown in Figure 2. The size of the specimen is 400 mm × 100 mm × 4 mm. A through sharp crack is fabricated in the center of the specimen by the fatigue testing machine, and its length is 22 mm. The specimen is subjected to the vertical uniform tensile load σ00 increasing from 0 to 3 MPa which has been estimated to ensure linear elasticity. The elastic modulus of PMMA is E = 3.24 GPa, and Poisson’s ratio is ν=0.35. The wavelength of the laser is λ = 532 nm. A 1024 × 768 pixel CCD camera with a pixel size of 6.4 μm is used to record digital holograms.
4. Experimental Results and Validation4.1. Experimental Results of Phase Difference
By using the image-plane reconstruction algorithm given in Section 2.1, we can obtain the phase difference of any applied load relative to the initial load (σ00=0). The wrapped phase differences at different load are shown in Figure 3. From these wrapped phase maps, it can be seen that there exists high speckle noise, the fringes near the crack tip are very close, and the phase data inside the crack region is invalid. So, these phase maps are very difficult to unwrap.
The wrapped phase of the object wave through the PMMA specimen with a crack. (a) σ00 = 1 MPa, (b) σ00 = 1.5 MPa, (c) σ00 = 2 MPa, and (d) σ00 = 2.5 MPa.
In order to unwrap the phase data in Figure 3, we use a mask to shield the crack region; then, CPULSI is used to unwrap the masked phase maps. The unwrapped phase maps are shown in Figure 4, and we can see that all the phase maps are unwrapped successfully.
The unwrapped phase of the object wave through the PMMA specimen with a crack. (a) σ00 = 1 MPa, (b) σ00 = 1.5 MPa, (c) σ00 = 2 MPa, and (d) σ00 = 2.5 MPa.
In order to illustrate the necessity of CPULSI used in this paper, four other different classical phase unwrapping algorithms were selected: Goldstein branch cut algorithm (namely, “Gold”), quality-guided path following algorithm (“Quality”), Flynn’s minimum-weight-discontinuity algorithm (“Flynn”), and minimum Lp norm algorithm (“Lp”) [34]. Figure 5 shows the unwrapped phase maps of the object wave as σ00 = 2 MPa obtained from these four algorithms. Comparing Figure 4(c), we can see that the unwrapped result from CPULSI is better than those from the other algorithms obviously. In fact, the other algorithms cannot successfully unwrap the phase data in this study because there all exist discontinuities in the maps of Figure 5. So, it is necessary to use CPULSI in this study.
The unwrapped phase of the object wave as σ00 = 2 MPa obtained from (a) Gold, (b) Quality, (c) Flynn, and (d) Lp.
4.2. Calculation of the Stress Field and SIF near the Crack Tip
In order to reduce the noise in the phase maps, median filtering with the kernel size of 21 × 21 pixels is applied on the unwrapped phase maps. The denoised unwrapped phase is used to calculate stress sum near the crack tip. Figure 6 shows the stress sum as σ00 = 1 MPa, 1.5 MPa, 2 MPa, and 2.5 MPa. From these results, it can be seen that the maximum stress sums exist in the zone near the crack tip and those near the crack plane close to zero which is in accordance with the fracture theory.
The sum of stress (σx+σy) in the PMMA specimen with a crack. (a) σ00 = 1 MPa, (b) σ00 = 1.5 MPa, (c) σ00 = 2 MPa, and (d) σ00 = 2.5 MPa.
According to (15), we can calculate SIFs from stress sums at different loads. Because (13) and (14) only express the stresses infinitely close to the crack tip in the range of linear elastic, and the K-dominate zone must be determined before calculation. Here, we select the dashed sectorial domain in Figure 6(a) which is 1.5 mm ≤ r ≤ 2.5 mm and −45° ≤ θ ≤ 45° as the K-dominate zone is referring to some literatures [14, 15, 18, 19]. The average value of the SIFs in the K-dominate zone is set as the measured value. The measured SIFs at different loads are shown in Table 1.
The comparison between the measured values and theoretical values of KI.
Tensile load (MPa)
Measured KI (MPa·m0.5)
Theoretical KI (MPa·m0.5)
Relative error (%)
1.0
0.193
0.191
1.0
1.5
0.291
0.287
1.4
2.0
0.379
0.383
−1.0
2.5
0.503
0.479
5.0
In order to validate the accuracy of the measurement, the theoretical values of the SIFs are calculated by using the method in [39]. For the crack in Figure 2, KI can be calculated as(16)KI=Fσ00πa,where σ00 is the vertical uniform tensile load, a is the half length of the crack, and F is the calibration coefficient for the size of the specimen which can be expressed as(17)F=1−0.025ab2+0.06ab4secπa2b,where b is the half width of the specimen. The calculated theoretical values of SIFs at different loads are also shown in Table 1. We can obtain the relative errors by comparing the measured values and theoretical ones. From Table 1, it can be seen that the error of the measured SIF does not exceed 5% which means the results of the proposed approach are accurate.
5. Conclusions
In this paper, we address the real-time and quantitative measurement for the SIFs of the Mode I crack in the transparent plate by DHI. For the phase processing near the crack tip which is affected by the close fringes and high speckle noise, we use phase unwrapping algorithm CPULSI and median filtering to recover the phase maps. The calculation method of KI from the phase difference is formulated based on linear elastic fracture mechanics. The proposed approach is used to measure KI of the crack in the PMMA specimen. The measured results are validated in comparison with the theoretical ones. The proposed approach can be used for the real-time and quantitative test of transparent materials.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Natural Science Foundation of China (11462009 and 11362007).
ZhaoF.-G.KongY.-T.XuZ.-W.YaoX.ZuoB.LiW.-S.High-performance flexible transparent conductive films achieved by cooperation between 1D copper nanowires and 2D graphene materials201755509551610.1039/c7tc00860k2-s2.0-85021680587YangB.YaoC.YuY.LiZ.WangX.Nature degradable, flexible, and transparent conductive substrates from green and earth-abundant materials20177110.1038/s41598-017-04969-y2-s2.0-85022345527TangQ.ShenH.YaoH.JiangY.ZhengC.GaoK.Preparation of silver nanowire/AZO composite film as a transparent conductive material20174311106111310.1016/j.ceramint.2016.10.0482-s2.0-84995463485WuM.SunD.TanC.TianX.HuangY.Al-doped ZnO monolayer as a promising transparent electrode material: a first-principles study201710435910.3390/ma100403592-s2.0-85017405824WroblewskiG.SwatowskaB.PowroznikW.JakubowskaM.StapinskiT.Structural and optical properties of spray coated carbon hybrid materials applied to transparent and flexible electrodes201720177342467210.1155/2017/34246722-s2.0-85028684031WimalanandaM. D. S. L.KimJ.-K.LeeJ.-M.Toward the ultra-transparent electrode by using patterned silver nanowire and graphene layered material201712591910.1016/j.carbon.2017.09.0402-s2.0-85029515224TrungT. Q.LeeN.-E.Materials and devices for transparent stretchable electronics2017592202222210.1039/c6tc05346g2-s2.0-85014510365YangC.SouchayD.KneiβM.Transparent flexible thermoelectric material based on non-toxic earth-abundant p-type copper iodide thin film201781607610.1038/ncomms160762-s2.0-85022178980KangD. J.ParkG. U.ParkH. Y.ImH.-G.A robust transparent encapsulation material: Silica nanoparticle-embedded epoxy hybrid nanocomposite201714410711310.1016/j.compscitech.2017.03.0322-s2.0-85015915514BenitezT.GómezS. Y.de OliveiraA. P. N.TravitzkyN.HotzaD.Transparent ceramic and glass-ceramic materials for armor applications20174316130311304610.1016/j.ceramint.2017.07.2052-s2.0-85026410815ShimG.-I.KimS.-H.AhnD.-L.Experimental and numerical evaluation of transparent bulletproof material for enhanced impact-energy absorption using strengthened glass/polymer composite20169715016110.1016/j.compositesb.2016.04.0782-s2.0-84968756919KhaleghianS.EmamiA.TehraniM.SoltaniN.Analysis of effective parameters for stress intensity factors in the contact problem between an asymmetric wedge and a half-plane using an experimental method of photoelasticity20134344745310.1016/j.matdes.2012.07.0382-s2.0-84864814708EllingsenM. D.KhannaS. K.Experimental investigation of static interfacial fracture in orthotropic polymer composite bimaterials using photoelasticity2010132202100710.1115/1.40002232-s2.0-77955319888MaK.XieH.Mixed-mode fracture investigation of PMMA with initial single/double crack(s) interference using phase-shifted coherent gradient sensing method20175929630710.1016/j.polymertesting.2017.02.0152-s2.0-85013940020MaK.XieH.A novel phase shifting technique of coherent gradient sensing method for measuring crack-tip K-dominance20175781239124810.1007/s11340-017-0293-z2-s2.0-85020286455YaoX. F.YehH. Y.XuW.Fracture investigation at V-notch tip using coherent gradient sensing (CGS)20064351189120010.1016/j.ijsolstr.2005.03.0432-s2.0-29244436694SundaramB. M.TippurH. V.Dynamic mixed-mode fracture behaviors of PMMA and polycarbonate201717618621210.1016/j.engfracmech.2017.02.0292-s2.0-85015687913SundaramB. M.TippurH. V.Dynamics of crack penetration vs. branching at a weak interface: an experimental study20169631233210.1016/j.jmps.2016.07.0202-s2.0-84981312100PeriasamyC.TippurH. V.Measurement of crack-tip and punch-tip transient deformations and stress intensity factors using digital gradient sensing technique20139818519910.1016/j.engfracmech.2012.09.0112-s2.0-84876989813MontrésorS.PicartP.SakharukO.MuravskyL.Error analysis for noise reduction in 3D deformation measurement with digital color holography201734B9B1510.1364/josab.34.0000b92-s2.0-85018735735KhodadadD.SinghA. K.PedriniG.SjödahlM.Full-field 3D deformation measurement: comparison between speckle phase and displacement evaluation201655277735774310.1364/ao.55.0077352-s2.0-84988811675OliveiraG. N.RodriguesD. M. C.NunesL. C. S.dos SantosP. A. M.Digital Fourier transform holography applied to investigate mechanical deformation in polymers201250121798180310.1016/j.optlaseng.2012.06.0162-s2.0-84865490241LédlV.PsotaP.KavánF.MatoušekO.MokrýP.Surface topography measurement by frequency sweeping digital holography2017567808781410.1364/ao.56.0078082-s2.0-85030027833DekiffM.BerssenbrüggeP.KemperB.DenzC.DirksenD.Simultaneous acquisition of 3D shape and deformation by combination of interferometric and correlation-based laser speckle metrology201564825484010.1364/boe.6.0048252-s2.0-84961801015HesselingC.HomeyerT.PeinkeJ.GülkerG.Particle depth position detection by 2D correlation in digital in-line holography201641214947495010.1364/ol.41.0049472-s2.0-84995442527DesseJ.-M.PicartP.OlchewskyF.Quantitative phase imaging in flows with high resolution holographic diffraction grating20152318237262373710.1364/oe.23.0237262-s2.0-84957608091StetsonK. A.Effect of changing speckles in digital holography on measurements of static and vibratory displacements201655164485448910.1364/ao.55.0044852-s2.0-84973926814PoittevinJ.PicartP.GautierF.PezeratC.Quality assessment of combined quantization-shot-noise-induced decorrelation noise in high-speed digital holographic metrology20152324309173093210.1364/oe.23.0309172-s2.0-84962076658EldridgeW. J.SheinfeldA.RinehartM. T.WaxA.Imaging deformation of adherent cells due to shear stress using quantitative phase imaging201641235235510.1364/ol.41.0003522-s2.0-84962107006MalekM.KhelfaH.PicartP.MounierD.PoilâneC.Microtomography imaging of an isolated plant fiber: a digital holographic approach2016553A111A12110.1364/ao.55.00a1112-s2.0-84962272783XiaH.GuoR.YanF.Simultaneous measurement of stress-optic constant and stress field of transparent plate by digital holographic interferometry201612724119741198110.1016/j.ijleo.2016.09.1042-s2.0-84990038555XiaH.GuoR.FanZ.ChengH.YangB.Non-invasive mechanical measurement for transparent objects by digital holographic interferometry based on iterative least-squares phase unwrapping201252443944510.1007/s11340-011-9516-x2-s2.0-84858446698GuoR.XiaH.FanZ.YangB.HeT.Research on the interaction of collective cracks with digital holography2009303839in ChineseGhigliaD. C.PrittM. D.1998Hoboken, NJ, USAWileyMontresorS.PicartP.Quantitative appraisal for noise reduction in digital holographic phase imaging20162413143221434310.1364/oe.24.0143222-s2.0-84977544310XiaH.MontresorS.GuoR.Phase calibration unwrapping algorithm for phase data corrupted by strong decorrelation speckle noise20162425287132873010.1364/oe.24.0287132-s2.0-85006147433VillarS. A.TorcidaS.AcostaG. G.Median filtering: a new insight201758113014610.1007/s10851-016-0694-02-s2.0-85007438897LiZ. N.2012Beijing, ChinaBehang University Pressin ChineseChinese Aeronautical Establishment1993Beijing, ChinaScience Pressin Chinese