Research Article  Open Access
Haiting Xia, Rongxin Guo, Feng Yan, Heming Cheng, "RealTime and Quantitative Measurement of CrackTip Stress Intensity Factors Using Digital Holographic Interferometry", Advances in Materials Science and Engineering, vol. 2018, Article ID 1954573, 8 pages, 2018. https://doi.org/10.1155/2018/1954573
RealTime and Quantitative Measurement of CrackTip Stress Intensity Factors Using Digital Holographic Interferometry
Abstract
Detection of the crack in an object is a critical problem for the health monitoring of a transparent object. The realtime and quantitative measurement of the cracktip stress intensity factor (SIF) remains an open issue. In this paper, an approach for realtime 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 leastsquares 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.
1. 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 fullfield 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 stressoptic 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 leastsquares 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 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 Approach
2.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 imageplane offaxis 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.
Establishing and coordinates in the CCD plane and coordinate normal to the plane, the object wave and the reference wave can be expressed aswhere , , is the wavelength of the laser, and are the angles between axis and the projection of the reference wave vector on the plane and plane, respectively.
The hologram captured by the CCD is described as
Performing Fourier transform on both sides of the above equation, we can get the spectrum of the hologram:
On the right side of this equation, the first item is the spectrum of the zeroorder diffracted wave, the second is the spectrum of the object wave which is modulated by the reference wave and translated to the point (, ) 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 and , we can separate the spectrum of the object wave from the other spectrums. A rectangular window filter is used to extract . Then, the inverse Fourier transform of is carried out to reconstruct the complex modulated object wave as follows:
The phase difference of object waves at two loading states, due to the specimen’s deformation, is expressed aswhere is the wrapped phase difference, and are the complex object waves at different loading states, and 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 leastsquares and iterations (CPULSI) [36] and median filtering to recover the true phase from the noisy wrapped phase. Here, represents the true phase and (, i.e., modulo ) represents the wrapped phase at the grid point () of a phase map. Their relationship can be expressed by the wrapping operator which is defined aswhere , is an integer, and are, respectively, the number of grid points with respect to the and indexes. The 1st order spatial wrapped phase derivatives are defined aswhere and are, respectively, the difference with respect to the and indexes.
In practice, the presence of noise will generate errors between noisy and noisefree 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:where is the signum function, and are the thresholds, and and are the calibrated phase derivatives. These parameters are defined as ( means the statistical average)
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 leastsquares sense, the optimal solution can be obtained from the discrete Poisson equation with the Neumann boundary conditions:where
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 leastsquares 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 leastsquares 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]where and are the normal stresses in the plane, is the thickness of the plate, is the stressoptic 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]where is the stress intensity factor (SIF) of the Mode I crack and and are the polar coordinates relative to the crack tip.
By adding (13) and (14), we can get the equation to calculate from the stresses as
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 increasing from 0 to 3 MPa which has been estimated to ensure linear elasticity. The elastic modulus of PMMA is = 3.24 GPa, and Poisson’s ratio is . 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 Validation
4.1. Experimental Results of Phase Difference
By using the imageplane reconstruction algorithm given in Section 2.1, we can obtain the phase difference of any applied load relative to the initial load (). 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.
(a)
(b)
(c)
(d)
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.
(a)
(b)
(c)
(d)
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”), qualityguided path following algorithm (“Quality”), Flynn’s minimumweightdiscontinuity algorithm (“Flynn”), and minimum Lp norm algorithm (“Lp”) [34]. Figure 5 shows the unwrapped phase maps of the object wave as = 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.
(a)
(b)
(c)
(d)
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 = 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.
(a)
(b)
(c)
(d)
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 Kdominate zone must be determined before calculation. Here, we select the dashed sectorial domain in Figure 6(a) which is 1.5 mm ≤ ≤ 2.5 mm and −45° ≤ ≤ 45° as the Kdominate zone is referring to some literatures [14, 15, 18, 19]. The average value of the SIFs in the Kdominate zone is set as the measured value. The measured SIFs at different loads are shown in Table 1.

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, can be calculated aswhere is the vertical uniform tensile load, is the half length of the crack, and is the calibration coefficient for the size of the specimen which can be expressed aswhere 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 realtime 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 from the phase difference is formulated based on linear elastic fracture mechanics. The proposed approach is used to measure 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 realtime 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).
References
 F.G. Zhao, Y.T. Kong, Z.W. Xu, X. Yao, B. Zuo, and W.S. Li, “Highperformance flexible transparent conductive films achieved by cooperation between 1D copper nanowires and 2D graphene materials,” Journal of Materials Chemistry C, vol. 5, pp. 5509–5516, 2017. View at: Publisher Site  Google Scholar
 B. Yang, C. Yao, Y. Yu, Z. Li, and X. Wang, “Nature degradable, flexible, and transparent conductive substrates from green and earthabundant materials,” Scientific Reports, vol. 7, no. 1, 2017. View at: Publisher Site  Google Scholar
 Q. Tang, H. Shen, H. Yao, Y. Jiang, C. Zheng, and K. Gao, “Preparation of silver nanowire/AZO composite film as a transparent conductive material,” Ceramics International, vol. 43, no. 1, pp. 1106–1113, 2017. View at: Publisher Site  Google Scholar
 M. Wu, D. Sun, C. Tan, X. Tian, and Y. Huang, “Aldoped ZnO monolayer as a promising transparent electrode material: a firstprinciples study,” Materials, vol. 10, no. 4, p. 359, 2017. View at: Publisher Site  Google Scholar
 G. Wroblewski, B. Swatowska, W. Powroznik, M. Jakubowska, and T. Stapinski, “Structural and optical properties of spray coated carbon hybrid materials applied to transparent and flexible electrodes,” Journal of Nanomaterials, vol. 2017, Article ID 3424672, 7 pages, 2017. View at: Publisher Site  Google Scholar
 M. D. S. L. Wimalananda, J.K. Kim, and J.M. Lee, “Toward the ultratransparent electrode by using patterned silver nanowire and graphene layered material,” Carbon, vol. 125, pp. 9–19, 2017. View at: Publisher Site  Google Scholar
 T. Q. Trung and N.E. Lee, “Materials and devices for transparent stretchable electronics,” Journal of Materials Chemistry C, vol. 5, no. 9, pp. 2202–2222, 2017. View at: Publisher Site  Google Scholar
 C. Yang, D. Souchay, M. Kneiβ et al., “Transparent flexible thermoelectric material based on nontoxic earthabundant ptype copper iodide thin film,” Nature Communications, vol. 8, p. 16076, 2017. View at: Publisher Site  Google Scholar
 D. J. Kang, G. U. Park, H. Y. Park, and H.G. Im, “A robust transparent encapsulation material: Silica nanoparticleembedded epoxy hybrid nanocomposite,” Composites Science and Technology, vol. 144, pp. 107–113, 2017. View at: Publisher Site  Google Scholar
 T. Benitez, S. Y. Gómez, A. P. N. de Oliveira, N. Travitzky, and D. Hotza, “Transparent ceramic and glassceramic materials for armor applications,” Ceramics International, vol. 43, no. 16, pp. 13031–13046, 2017. View at: Publisher Site  Google Scholar
 G.I. Shim, S.H. Kim, D.L. Ahn et al., “Experimental and numerical evaluation of transparent bulletproof material for enhanced impactenergy absorption using strengthened glass/polymer composite,” Composites Part B, vol. 97, pp. 150–161, 2016. View at: Publisher Site  Google Scholar
 S. Khaleghian, A. Emami, M. Tehrani, and N. Soltani, “Analysis of effective parameters for stress intensity factors in the contact problem between an asymmetric wedge and a halfplane using an experimental method of photoelasticity,” Materials and Design, vol. 43, pp. 447–453, 2013. View at: Publisher Site  Google Scholar
 M. D. Ellingsen and S. K. Khanna, “Experimental investigation of static interfacial fracture in orthotropic polymer composite bimaterials using photoelasticity,” Journal of Engineering Materials and Technology, vol. 132, no. 2, p. 021007, 2010. View at: Publisher Site  Google Scholar
 K. Ma and H. Xie, “Mixedmode fracture investigation of PMMA with initial single/double crack(s) interference using phaseshifted coherent gradient sensing method,” Polymer Testing, vol. 59, pp. 296–307, 2017. View at: Publisher Site  Google Scholar
 K. Ma and H. Xie, “A novel phase shifting technique of coherent gradient sensing method for measuring cracktip Kdominance,” Experimental Mechanics, vol. 57, no. 8, pp. 1239–1248, 2017. View at: Publisher Site  Google Scholar
 X. F. Yao, H. Y. Yeh, and W. Xu, “Fracture investigation at Vnotch tip using coherent gradient sensing (CGS),” International Journal of Solids and Structures, vol. 43, no. 5, pp. 1189–1200, 2006. View at: Publisher Site  Google Scholar
 B. M. Sundaram and H. V. Tippur, “Dynamic mixedmode fracture behaviors of PMMA and polycarbonate,” Engineering Fracture Mechanics, vol. 176, pp. 186–212, 2017. View at: Publisher Site  Google Scholar
 B. M. Sundaram and H. V. Tippur, “Dynamics of crack penetration vs. branching at a weak interface: an experimental study,” Journal of the Mechanics and Physics of Solids, vol. 96, pp. 312–332, 2016. View at: Publisher Site  Google Scholar
 C. Periasamy and H. V. Tippur, “Measurement of cracktip and punchtip transient deformations and stress intensity factors using digital gradient sensing technique,” Engineering Fracture Mechanics, vol. 98, pp. 185–199, 2013. View at: Publisher Site  Google Scholar
 S. Montrésor, P. Picart, O. Sakharuk, and L. Muravsky, “Error analysis for noise reduction in 3D deformation measurement with digital color holography,” Journal of the Optical Society of America B, vol. 34, pp. B9–B15, 2017. View at: Publisher Site  Google Scholar
 D. Khodadad, A. K. Singh, G. Pedrini, and M. Sjödahl, “Fullfield 3D deformation measurement: comparison between speckle phase and displacement evaluation,” Applied Optics, vol. 55, no. 27, pp. 7735–7743, 2016. View at: Publisher Site  Google Scholar
 G. N. Oliveira, D. M. C. Rodrigues, L. C. S. Nunes, and P. A. M. dos Santos, “Digital Fourier transform holography applied to investigate mechanical deformation in polymers,” Optics and Lasers in Engineering, vol. 50, no. 12, pp. 1798–1803, 2012. View at: Publisher Site  Google Scholar
 V. Lédl, P. Psota, F. Kaván, O. Matoušek, and P. Mokrý, “Surface topography measurement by frequency sweeping digital holography,” Applied Optics, vol. 56, pp. 7808–7814, 2017. View at: Publisher Site  Google Scholar
 M. Dekiff, P. Berssenbrügge, B. Kemper, C. Denz, and D. Dirksen, “Simultaneous acquisition of 3D shape and deformation by combination of interferometric and correlationbased laser speckle metrology,” Biomedical Optics Express, vol. 6, pp. 4825–4840, 2015. View at: Publisher Site  Google Scholar
 C. Hesseling, T. Homeyer, J. Peinke, and G. Gülker, “Particle depth position detection by 2D correlation in digital inline holography,” Optics Letters, vol. 41, no. 21, pp. 4947–4950, 2016. View at: Publisher Site  Google Scholar
 J.M. Desse, P. Picart, and F. Olchewsky, “Quantitative phase imaging in flows with high resolution holographic diffraction grating,” Optics Express, vol. 23, no. 18, pp. 23726–23737, 2015. View at: Publisher Site  Google Scholar
 K. A. Stetson, “Effect of changing speckles in digital holography on measurements of static and vibratory displacements,” Applied Optics, vol. 55, no. 16, pp. 4485–4489, 2016. View at: Publisher Site  Google Scholar
 J. Poittevin, P. Picart, F. Gautier, and C. Pezerat, “Quality assessment of combined quantizationshotnoiseinduced decorrelation noise in highspeed digital holographic metrology,” Optics Express, vol. 23, no. 24, pp. 30917–30932, 2015. View at: Publisher Site  Google Scholar
 W. J. Eldridge, A. Sheinfeld, M. T. Rinehart, and A. Wax, “Imaging deformation of adherent cells due to shear stress using quantitative phase imaging,” Optics Letters, vol. 41, no. 2, pp. 352–355, 2016. View at: Publisher Site  Google Scholar
 M. Malek, H. Khelfa, P. Picart, D. Mounier, and C. Poilâne, “Microtomography imaging of an isolated plant fiber: a digital holographic approach,” Applied Optics, vol. 55, no. 3, pp. A111–A121, 2016. View at: Publisher Site  Google Scholar
 H. Xia, R. Guo, F. Yan et al., “Simultaneous measurement of stressoptic constant and stress field of transparent plate by digital holographic interferometry,” OptikInternational Journal for Light and Electron Optics, vol. 127, no. 24, pp. 11974–11981, 2016. View at: Publisher Site  Google Scholar
 H. Xia, R. Guo, Z. Fan, H. Cheng, and B. Yang, “Noninvasive mechanical measurement for transparent objects by digital holographic interferometry based on iterative leastsquares phase unwrapping,” Experimental Mechanics, vol. 52, no. 4, pp. 439–445, 2012. View at: Publisher Site  Google Scholar
 R. Guo, H. Xia, Z. Fan, B. Yang, and T. He, “Research on the interaction of collective cracks with digital holography,” Laser Journal, vol. 30, pp. 3839, 2009, in Chinese. View at: Google Scholar
 D. C. Ghiglia and M. D. Pritt, TwoDimensional Phase Unwrapping: Theory, Algorithms, and Software, Wiley, Hoboken, NJ, USA, 1998.
 S. Montresor and P. Picart, “Quantitative appraisal for noise reduction in digital holographic phase imaging,” Optics Express, vol. 24, no. 13, pp. 14322–14343, 2016. View at: Publisher Site  Google Scholar
 H. Xia, S. Montresor, R. Guo et al., “Phase calibration unwrapping algorithm for phase data corrupted by strong decorrelation speckle noise,” Optics Express, vol. 24, no. 25, pp. 28713–28730, 2016. View at: Publisher Site  Google Scholar
 S. A. Villar, S. Torcida, and G. G. Acosta, “Median filtering: a new insight,” Journal of Mathematical Imaging and Vision, vol. 58, no. 1, pp. 130–146, 2017. View at: Publisher Site  Google Scholar
 Z. N. Li, Applied Fracture Mechanics, Behang University Press, Beijing, China, 2012, in Chinese.
 Chinese Aeronautical Establishment, Handbook of StressIntensity FactorsRevised Edition, Science Press, Beijing, China, 1993, in Chinese.
Copyright
Copyright © 2018 Haiting Xia et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.