# Real-Time and Quantitative Measurement of Crack-Tip Stress Intensity Factors Using Digital Holographic Interferometry

Authors: Haiting Xia; Rongxin Guo; Feng Yan; Heming Cheng
Journal: Advances in Materials Science and Engineering
 (2018)

---

## Abstract

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.


---

## Body

## 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 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 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 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.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.Establishingx 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 Figure2, 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.Figure 2
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.


## 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 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.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.Establishingx 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 Figure2, 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.Figure 2
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 Figure2. 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 Validation

### 4.1. Experimental Results of Phase Difference

By using the image-plane reconstruction algorithm given in Section2.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.Figure 3
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.
(a)
(b)
(c)
(d)In order to unwrap the phase data in Figure3, 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.Figure 4
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.
(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”), 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.Figure 5
The unwrapped phase of the object wave asσ00 = 2 MPa obtained from (a) Gold, (b) Quality, (c) Flynn, and (d) Lp.
(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. Figure6 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.Figure 6
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.
(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 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.Table 1
The comparison between the measured values and theoretical values ofKI.
Tensile load (MPa)
MeasuredKI (MPa·m0.5)
TheoreticalKI (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.0In 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.


## 4.1. Experimental Results of Phase Difference

By using the image-plane reconstruction algorithm given in Section2.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.Figure 3
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.
(a)
(b)
(c)
(d)In order to unwrap the phase data in Figure3, 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.Figure 4
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.
(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”), 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.Figure 5
The unwrapped phase of the object wave asσ00 = 2 MPa obtained from (a) Gold, (b) Quality, (c) Flynn, and (d) Lp.
(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. Figure6 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.Figure 6
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.
(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 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.Table 1
The comparison between the measured values and theoretical values ofKI.
Tensile load (MPa)
MeasuredKI (MPa·m0.5)
TheoreticalKI (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.0In 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 ofKI 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.



---

*Source: 1954573-2018-05-16.xml*