An Analytical Model for Defect Depth Estimation Using Pulsed Thermography

The use of pulsed thermography as a non-destructive evaluation tool for damage monitoring of composite materials has dramatically increased in the past decade. Typically, optical flashes are used as external heating sources, which may cause poor defect definition especially for thicker materials or multiple delaminations. SMArt thermography is a new alternative to standard pulsed thermography as it overcomes the limitations on the use of external thermal sources. Such a novel technology enables a built-in, fast and in-depth assessment of both surface and internal material defects by embedding shape memory alloy wires in traditional carbon fibre reinforced composite laminates. However, a theoretical model of thermal wave propagation for SMArt thermography, especially in the presence of internal structural defects, is needed to better interpret the observations/data measured during the experiments. The objective of this paper was to develop an analytical model for SMArt thermography to predict the depth of flaws/damage within composite materials based on experimental data. This model can also be used to predict the temperature contrast on the surface of the laminate, accounting for defect depth, size and opening, thermal properties of material and defect filler, thickness of the component, and intensity of the excitation energy. The results showed that the analytical model gives good predictions compared to experimental data. This paper is one of the first pioneering work showing the use thermography as a quantitative non-destructive tool where defect size and depth could be assessed with good accuracy.


Introduction
Flash thermography is the most commonly used thermographic non-destructive evaluation (TNDE) technique (e. g. [1,2]). Flash thermograhy is an active TNDE approach [3], which consists in using a number of high intensity flash lamps to heat the surface of the part under inspection. The result is an instantaneous rise in temperature of the surface followed by a cooling as the heat applied to the surface is conducted into the part. The cooling of the surface is monitored by an infrared (IR) camera and the video images are stored in a personal computer (PC) for viewing at the end of the test.
Detection of delamination defects between plies is of the greatest importance for composite laminates, as they can significantly reduce the strength and performance of structures under compressive in-plane loading, eventually giving rise after bucking to global plate instability [4][5][6].
Adopting flash thermograhy defect sizing from the thermal images has been performed employing both analytical, experimental and numerical techniques by [7][8][9][10][11]. Especially the quantitative prediction of defect depth using flash thermography is an important on-going research topic. The methods reported in literature for this purpose can be distinguished in empirical and experimental. The main techniques were reviewed and contrasted by [12,13], who used samples with flat-bottom holes as simulated defects in their experiments. Experimental studies [14] have showed that delaminations with a ratio of diameter to depth of the delamination less than 0.5 cannot be detected using thermography.
An analytical model for flash thermography has been developed in [15] which has been found to provide good results when compared against both experimental and finite difference simulations [16]. This model is based on the impulse heating transient response of an adiabatically isolated layer [17,18] where the thermal response function of the layer is given by the thermal wave interference expression obtained in [19].
Flash thermography presents a series of complications caused by the use of external lamps.
• The positioning of the lamp(s) from the part under inspection can have a significant impact on test results. • The heating from a single lamp will be non-uniform, peaking at the center. A more uniform field of heating can be produced by two or more lamps imposing potentially a substantial increase in the cost of the test equipment. • Lamp filaments and tubes can continue to emit infrared radiation (lamp after glow) for a significant amount of time (seconds) after the lamp has been switched off which may be detected by the IR camera as it is reflected from the test piece or other surfaces into the camera lens. The result can be the generationof an unwanted IR fog that can obscure required IR images of the surface of the part under inspection. The effect is extremely important for parts whose defect response times are very short (<1s) as on this time scale after glow will be very significant. • The optical characteristics of the surface of the material or its coating heavily impact the results of the inspection. Many materials for example are translucent, i. e. light is absorbed in them over a range of depths (and consequently the heat source obtained by flash excitation is not confined to the surface; it is distributed over a range of depths), but their optical absorption coefficients are not readily available.
In this paper a novel approach to pulsed thermography, which does not require external lamps, but is instead based on SMArt composites, is considered. A SMArt composite is a new kind of smart multifunctional material obtained by embedding shape memory alloy (SMA) wires within traditional carbon fibre-reinforced composite (CFRP) composites [20].
Smart materials are materials designed to adapt to the environment in a controlled fashion following the modification of one or more of their properties by external stimuli [21].
Multifunctional materials lead to optimal system performance, by combining different functions into a single material, otherwise not possible through independent subsystem optimization [22].
Specifically, the embodiment of SMA within the lay-up of CFRP laminates has been widely studied as a valid manufacturing procedure to enhance the impact resistance of traditional composite structures due to their unique physical properties such as shape memory effect and superelasticity [23]. The electrical variation of embedded SMA wires to the strain distribution within hybrid glass-fiber reinforced plastic (GFRP) laminates has been correlated by [24], while a hybrid CFRP/SMA composite with damages suppression function that is enabled by activating the martensite-austenite transformation was manufactured by [25].
However, to this date, only few works have investigated the possibility of exploiting the presence of an internal hybrid grid to enable the laminate to have multiple additional features, hence resulting in a real multifunctional system characterized by specific properties that go further than the traditional load-bearing functionality. Based on these premises, a SMA based multifunctional composite (named SMArt) was developed exploiting their intrinsic thermoelectric properties for applications such as sensing, thermography and de-icing.
In particular SMArt thermography is a form of material enabled thermography [26] where the embedded SMA wires are used as heat sources to generate in situ power resistive heating (ohmic heating or joule heating) that can be used to perform pulsed thermography. Although good results in term of damage detection have been already obtained with this new technique, there is still the need for a rigorous theoretical frame in which this phenomenon must be included. Such a model should not be limited to the inclusion of SMA, but should be extendable to any kind of conductive wire embedded into the traditional lay-up of a composite laminate.
In this paper an analytical model for pulsed thermography using SMArt composites will be developed specifically for in-plane, delamination-like defects. This model can be used to predict the temperature contrast on the surface of the laminate, accounting for 1. defect depth, size and opening; 2. thermal properties of material and defect filler; 3. thickness of the component; 4. depth of the SMA wire and intensity of the excitation energy.
Moreover, using this model an analytical expression for evaluating the defect depth will be obtained.
Finally, the results of the analytical model and defect depth evaluation are contrasted against experimental results.
"Methods" describes the theory behind the analytical model and its application to evaluating analytically the defect depth, while "Results" compares the model against some experimental data. "Conclusions" provides some conclusions from this work and summarizes the major findings.

Methods
Embedding an SMA wire at a depth h in a composite laminated plate of thickness L and conducting a impulse of resistive heating through this wire, it is possible to identify the presence of a defect using infrared imaging of the variations of the temperature contrast on the surface of the composite laminate. Here we are assuming an inplane, delamination-like circular defect of diameter D at a depth d from the surface of the material, as described in Fig. 1.
No assumptions are made for the laminated fiberreinforced composite in terms of number of layers and direction of the reinforcement fibers in each layer, as in the following (see "Derivation of an analytical model") only homogenized equivalent thermal properties are being considered. Fig. 1 Identification of an in-plane, delamination-like, circular defect in a composite laminate through pulsed thermography using an SMA wire as heat source. Note that an SMA wire (in blue) is positioned at a depth h from the surface of the composite layup, while a delaminationlike circular defect (in red), of diameter D, is positioned at a depth d from the surface. The Cartesian coordinate axes (x, y, z) are oriented such that x, y are in-plane coordinates, with the y-axis running along the SMA wire, and z is the transverse coordinate. T (h, t) is the temperature rise on the surface of the material at time t An analytical model for pulsed thermography using an SMA wire as heat source will be derived in the "Derivation of an analytical model". Then "Evaluation of Defect Depth" will explore the possibility of using this model to obtain an approximate expression for evaluating analytically the defect depth.

Derivation of an analytical model
Consider the laminated fiber-reinforced composite plate shown in Fig. 1, in which the origin of the coordinate system is placed at the center of the SMA wire, which will be assumed to all effects and purposes to act as a line source, and the components of the Cartesian coordinates are (x, y, z) where x, y are in-plane coordinates, and z is the transverse coordinate. Furthermore, it is assumed that the SMA wire runs along the y-axis.
The three-dimensional inhomogeneous linear problem of finding the temperature field T (r, t), which results from the transient heat conduction within the composite laminate due to an internal heat source g (r, t), is given in equation where r = xî + yĵ + zn, ⊂ R 2 is the domain of influence, S is the boundary of the composite laminate and ∇ 2 = ∂ 2 ∂x 2 î + ∂ 2 ∂y 2 ĵ + ∂ 2 ∂z 2 n is the Laplacian operator. In equation (2.1) k represents the thermal conductivity of the composite material (W m −1 K −1 ), whilst α = k ρc is the thermal diffusivity of the composite material (m −2 ) with ρ and c the density (kg m −3 ) and the specific heat capacity (W kg −1 K −1 ) of the composite material, respectively.
The thermal properties in equation (2.1) are the homogenized equivalent thermal properties of the laminated fiberreinforced composite. (In the case of unidirectionallyreinforced fiber composites, the physical properties are anisotropic, but it is common to assume these materials to be transversely isotropic. The anisotropy in thermal conductivity will be accounted for, in the following, using a specific parameter, the thermal diffusivity anisotropy, see equation (2.5) on page 13, whose effects have been studied by one of the authors in [15].) After resistive heating, if a uniform impulse of thermal energy is released instantaneously by the SMA wire at time t = 0 and x = z = 0, but along the entire y-axis in a circular cylindrical system of coordinates in the amount (per unit length) Q 0 Jm −1 , then it is a well known result (e. g. [27]) that the temperature is independent of y and the corresponding fundamental solution is given as where r = √ x 2 + z 2 and T (r, t) is the temperature rise at a distance r from the SMA wire at a time t after the impulse of heat energy has been released. We are interested in the cross-sectional propagation of the transient thermal wave (i. e. along the z-axis), so the corresponding one-dimensional (1D) field solution of problem equation (2.1) is independent of x and y.
When there is no defect in the component being inspected the heat released by the SMA wire will propagate straight through the laminated fiber-reinforced composite plate without impediments causing over time a temperature rise on the surface at a distance h from the wire. The temperature rise T nd (h, t) on the surface of the non-defective material will consist of two contributions.
1. The first contribution T F nd (h, t) is constituted by the forward wave, which reaches directly the surface, and all its subsequent reflections at the back face of the laminate of thickness L.
2. The second one T R nd (h, t) instead is composed by the reverse wave, which reaches the surface only after being reflected at the back face of the laminate, and all its subsequent reflections by the same back face.
Note that in theory there are an infinite number of reverberations, but in practice a summation over the first six terms is adequate. The 1D analytical model of the temperature rise (or background temperature) on the surface of the non-defective material T nd (h, t) is then given by equation (2.3).
Note that the in equation (2.3) the magnitude of the heating produced at the surface varies inversely with the thermal conductivity k of the composite material being impulse heated and that this heat decays with time t.
As a delamination-like defect lying beneath the surface in a plane parallel to the surface is being considered (see Fig. 1), then both the magnitude and the time dependence of the temperature rise on the laminate surface, T (h, t), will be altered by the defect, as the conduction of the heat released by the SMA wire to the surface will be reduced or blocked by the defect causing a slower temperature rise than indicated by equation (2.3). This results in the area of the surface over a defect appearing cooler in the IR video images collected during the test, enabling the detection of the defect.
As a first approximation, a region containing a delamination-like defect can be treated as a layer of thickness d, the depth of the defect, below the surface.
When a defect is present in the component under inspection, a part of the heat, generated by the the SMA wire and being conducted towards the surface will be blocked by the defect, where is the effective thermal reflectivity of the defect, assumed 100 % for a wide open defect. This applies to both direct and reverse waves and all their reflections. Then only a fraction (1 − ) of the heat will therefore be able to reach the surface. Another fraction of this heat reverberates at the surface after being reflected by the defect having a round-trip path of length 2d, with following terms having round-trip path lengths of 4d, 6d, etc. The 1D analytical model of the temperature rise T d (h, t) on the surface of the defective material is then given by equation (2.4).
Real defects are not infinite in size, but are finite in their lateral dimensions, so heat flowing in their vicinity cannot be assumed to be a 1D phenomenon. As a circular defect of diameter D has been assumed, then the thermal lateral diffusion around the defect will also contribute to the defective temperature rise over time on the surface of the plate [15]. The physical assumption here is that the lateral diffusion of heat from the edge of the defect to the center will cause an increase with time in the temperature on the surface of the laminate over the center of a defect. For a circular defect of diameter D, the diffusion distance is D/2.
Accounting for lateral diffusion, then the temperature rise on the surface of the defective material (equation (2.4)) can be rewritten as where A is the thermal diffusivity anisotropy of the composite material. This type of anisotropy assumes a prominent importance on the temperature contrast caused by defects for materials that are thermally anisotropic, such as composites, for which in-plane thermal conductivity typically exceeds through-the-thickness conductivity because of the layered structure of composites. For this reason, the inplane thermal diffusivity of composites is larger by a factor of ≈ 3-5 than the through-the-thickness thermal diffusivity. For thermally anisotropic materials the diffusivity should include the anisotropy factor A in the exponential term, while for the thermally isotropic ones A = 1 in equation (2.5). The exponential term in equation (2.5), which multiplies the 1D terms, accounts for the physics of the diffusion of heat from the edge of the circular defect to the center, a distance D/2 away. In practice, (equation (2.5)) can be used to provide the temperature rise at any distance from the edge of the defect. This fact will be used in "Results".
The temperature contrast T c (h, t) on the surface of the material is given by equation (2.6) The temperature contrast at the surface over a defective region is negative, it decreases with time until it reaches a minimum, after which it increases tending to zero. A thermal contrast image of the defect will be obtained, only if 1. the magnitude of the minimum temperature contrast is significantly above the noise level of the IR camera (typically 0.02K), and 2. the minimum contrast time is long enough to be recorded by the camera (typical frame rate 60Hz).
Examples of thermal images, along with a further discussion, will be presented in "Results". Effects of defect opening must also be accounted for, as real defects have a finite opening that may range from less than a μm to several mm. The thermal barriers presented by such defects are strongly dependent on defect opening and the thermal properties of the host material. Defects are usually treated as thermal contact resistances, R = l d /k d , where l d is the defect opening and k d is the thermal conductivity of the material filling the defect (usually it is assumed that defects are air filled).
In previous work [28] on thermal wave interference an expression was obtained for the thermal reflectivity of a defect represented as a thermal contact resistance.
In expression (equation (2.7)) k is the thermal conductivity of the host material and σ is the thermal wave number αt is equivalent to μ for transient thermal phenomena, where α is the thermal diffusivity of the host material. Making use of this in equation (2.7) leads to expression (equation (2.8)).
Expression (equation (2.8)) shows an effective thermal reflectivity of a defect that is a function of its thermal contact resistance R and the thermal properties k, ρ, c, of the host material. The results of using this expression have been found to be in good agreement with numerical modeling studies of the effects of defect opening on temperature contrast [10]. Expression (equation (2.8)) is used in (equations (2.4) and (2.5)) to compute the temperature contrast of defects of specified openings.
The temperature contrast can change dramatically depending on the defect opening, so this is a very important parameter whose effects will be analyzed in "Results".

Evaluation of Defect Depth
The possibility of obtaining from the model derived in "Derivation of an analytical model" an approximate expression for evaluating analytically the defect depth will be explored in the present subsection.
Considering in equation (2.3) only the first term in the summations in m, i. e. m = 0, and in equation (2.4) only the first term in the summations in m and two terms in the summations in n, i. e. m = 0 and n = 0, 1, then the temperature contrast T c (h, t) can be rewritten as Neglecting the second and the sixth terms on the righthand side of equation (2.9), an approximate expression for the temperature contrast is given by 16αAt (2.10) As discussed in "Derivation of an analytical model", the temperature contrast on the surface of the laminate over a defective area reaches a minimum T c = T c min at a certain time t = t min after the pulsed heat excitation has been released by the SMA wire.
Assuming known material and defect thermal properties, geometry of the sample and defect diameter and opening, and substituting with the expression in equation (2.8) and the values T c = T c min and t = t min derived from the experimental results in equation (2.10), then the following quadratic equation in the defect depth d, the only remaining unknown, is obtained with solutions 16πT cmin k (t min ) 2 q 1 q 2 q 3 + 4πR 2 T cmin kt min ζ 2 q 1 q 2 q 3 + 16πRT cmin k (t min ) (3/2) (2.14) Only the positive solution d 1 in equation (2.12) makes physical sense, and can be used to evaluate analytically the defect depth, assumed known all other parameters.

Results
In this section the results of the analytical model and the approximate expression for evaluating analytically the defect depth, which were presented in "Methods", are compared against experimental data. The analytical model was implemented in MATLAB [29], which is a numerical computing scripting language notably suitable for matrix manipulations.
The experimental results are for a four-ply laminated composite plate with orientations [0 • 90 • 90 • 0 • ] obtained from a composite T 700/M21 unidirectional prepreg with fiber volume fraction ≈ 57-59 %. Each lamina has a thickness of ≈ 0.1375mm. The dimensions of the plate were 10× 6 × 0.55mm 3 (width × height × thickness). A NiTi SMA wire with a diameter of 350μm was positioned between the 3rd and the 4th plies, while the presence of a defective area was modeled using a 11cm 2 Teflon patch of ≈ 0.05mm thickness. The inclusion of Teflon inserts is a well-known technique for introducing artificial in-plane delaminations in a laminate composite structure [20,26,30].
A pulse of electric potential difference V was applied at the ends of the SMA wire in order to induce resistive heating in the component under test for a short heating period. The effects of the heating were then captured by the thermal camera in a thermographic system. The thermal camera was an InSb electrically cooled infrared camera with a noise-equivalent temperature difference of ≈ 18-25mK and a resolution of 320 × 240 pixels (width × height). The camera was used at a frame rate of 50Hz.
The heating period was varied in the experiments between 10-50ms in steps of 10ms in order to characterize the impulsive nature of the analytical model. A V of 2V was applied in all tests. An electrical resistance of 2.5 was measured in the SMA wire. All tests were conducted at ambient temperature (≈ 25 • ). Figure 2 shows a comparison between the infrared images of the experimental data for the 30ms heating period at different time instants after the impulsive heating excitation (which in this case occurred at t = 3.1s).
The following signal processing techniques were adopted in order to visualize the defect in the thermal images. Fig. 2 between the pre-excitation and post-excitation thermal images. Table 1 Thermal properties of CFRP used in the analytical model

Background subtraction was applied to the images in
0.50 1200 1700 5.88 4 2. Thermal images were smoothed locally using simple moving averaging. 3. Moreover, the digital levels of the visualized thermal images were also rescaled based on a user selectable region of interest (ROI), which is represented by a green rectangle in the Fig. 2(a)-(d). Figure 2(a) refers to the thermal image at a time immediately after the pulsed heat excitation, while ( Fig. 2(b)) is for a time t=0.18s after the pulsed excitation where the temperature contrast reaches its minimum value. The presence of the patch is clearly visible in Fig. 2(b) as a break in the temperature rise in the thermal image. Figure 2(c) and (d) show how at times much later there is a temperature rise also on the surface area over the patch caused by heat diffusion both through and around the defective area. Figure 3(a) shows the results for the experimental data for the 30ms heating period. All graphs in Fig. 3(a) are plotted starting from when the SMA wire starts heating, after V has been applied, and were smoothed locally using simple moving averaging. The non-defective temperature rise is being measured on the surface of the plate in a small nondefective area directly on top of the SMA wire, so that it is uniform as much as possible. After the initial rise the nondefective temperature decreases because of the effects of thermal convection towards the air surrounding the sample and thermal conduction in-plane to the SMA wire. The temperature is given in Kelvin. The non-defective temperature increases by ≈ 1.50 Kelvin. The defective temperature rise is measured on the surface of the plate in a small defective area directly on top of the SMA wire, so that it is uniform as much as possible. The graph of the defective temperature rise takes into account non only the effects of thermal diffusion through the defect, but also the effects of lateral thermal diffusion around the defect. The defective temperature rise is of ≈ 0.60 Kelvin. The temperature contrast in Fig. 3(a) is negative reaching a minimum at 0.20s. The temperature contrast reduction is of ≈ 1.60 Kelvin. Figure 3(b) shows the results from the analytical model (equation (2.6)) corresponding to the experimental data in Fig. 3(a).
For comparison with the experimental data, the analytical model is computed using a plate thickness L of 0.55mm and a wire depth h of 0.4125mm (i. .e. 3 × 0.1375mm, which is the thickness of a single lamina). A defect depth d of 0.1375mm (i. .e. the thickness of a single lamina) and a defect diameter D of 10mm were assumed. Furthermore, to allow for imperfections in the bonding between patch and host material during the manufacturing process, it was assumed that the defect opening was 100mm and that the defect filler was air. The thermal properties utilized for the laminated fiber-reinforced composite plate are given in Table 1.
The values in Table 1 were derived by fitting the results of the analytical model to those of the experimental data. In particular, the procedure consisted in two stages.

First
The through-the-thickness thermal diffusivity α was derived by fitting the temperature contrast  Fig. 4(b) with the defect opening 100μm the minimum temperature contrast and time are T c min = −1.66K and t min = 0.21s, which are respectively ≈ 3 % and ≈ 10 % greater than for the same curve in Fig. 4(a)  Fig. 5 Temperature contrast versus elapsed time computed using the approximate expression (equation (2.10)) for the same parameters utilized in Fig. 3(b). In this case the minimum temperature contrast T c min = −1.61K at time t min = 0.20s obtained experimentally to the one computed using equation (2.6) for a value of = 1, which corresponds to a null contribution of the defective temperature rise. Second The thermal anisotropy A was derived by fitting the defective temperature rise obtained experimentally to the one computed from equation (2.5). The fitting was performed at a distance D/2=2.5mm from the edge of the defect, where the lateral diffusion was seen to be especially strong in the thermal images.
The non-defective and defective temperatures in Fig. 3(a) are constructed using, respectively, m = 0, . . . , 10 in equation (2.3) and m = 0, . . . , 10 and n = 0, . . . , 10 in equation (2.4). The temperature contrast computed by the analytical model agrees quite well the the experimental results in Fig. 3(a), with the model predicting a trough in the temperature contrast of −1.60K at 0.19s.
In Table 2 the results predicted by the analytical model (M) (equation (2.6)) and its approximate expression (AM) (equation (2.10)) for the minimum contrast time t min are compared against the experimental ones for 10ms, 20ms, 30ms, 40ms and 50ms heating periods. Figure 4(a) shows the temperature contrast versus elapsed time curves computed using the analytical model (equation (2.6)) for four defect openings: wide open (2.5cm), 100μm, 10μm and 1μm, holding fixed all other parameters. These curves are compared with a noise margin of 0.05K that is taken as being the minimum temperature contrast necessary to produce a useful thermal image of a defect. The results change significantly depending on the defect opening being considered, as was also found in [10], with the minimum contrast temperature and time in these curves increasing with defect depth. Figure 4(b) shows the temperature contrast versus elapsed time curves for same parameters in Fig. 4(a) but doubling the defect depth to 0.2750mm. In this case for the defect opening 100mm the minimum temperature contrast is T c min = −1.66K at time t min = 0.21s which are respectively ≈ 3 % and ≈ 10 % greater than for the same curve in Fig. 4(a), which gives an indication of the parameter sensitivity of the minimum contrast temperature and position for changes in the defect depth. The sensitivity of the temperature contrast and position to changes in defect depth using the approximate expression (equation (2.10)) will be analyzed in Fig. 6.
In Fig. 4 defects with all of the four openings exceed the noise margin, indicating SMArt thermograhy to be suitable as a TNDE technique for the inspection, but the sensitivity of the contrast to defect opening is evident, with the temperature contrast curve for the 1mm opening size just above the noise margin. Figure 5 shows the temperature contrast versus elapsed time computed using the approximate expression (equation (2.10)) for the same parameters utilized in Fig. 3(b). Fig. 6 Sensitivity of minimum contrast temperature and time versus defect depth using the approximate expression (equation (2.10)) for the two defect depths in Fig. 4. Note how minimum contrast temperature and time are not very sensitive to variations of the defect depth using this approximate expression, as a ≈ 50 % variation in defect depth causes only a couple of percentage points variation in these parameters In this case the minimum temperature contrast is T c min = −1.61K at time t min = 0.20s. Plugging these values into the expression of d 1 given by equation (2.12) gives the correct answer for the defect depth 1.375 × 10 −4 m. Figure 6 shows the sensitivity of minimum contrast temperature and time to variations in the defect depth using the approximate expression (equation (2.10)) for the two defect depths in Fig. 4. Minimum contrast temperature and time were not found to be very sensitive to variations of the defect depth using this approximate expression, as a ≈ 50 % variation in defect depth causes only a couple of percentage points variation in these parameters. For this reason, small errors in the minimum contrast temperature and time might produce large errors in the defect depth when computed analytically using the expression (equation (2.12)). Figure 7 shows the results predicted for the defect depth d by d 1 in equation (2.12) from the experimental data for the 10ms, 20ms, 30ms and 40ms heating periods. These graphs are plotted against the width of the plate, as this is the dimension along which the SMA wire is running through the plate (i. e. the abscissas in Fig. 2). Contrast ( Fig. 7(a)-(d)) to Fig. 2(b) to have a visual comparison of the defect size and location. The differences between the maximum predicted value in Fig. 7(a)-(d), compared to the assumed value ≈ 0.1375mm, range between 1-7 %.
The defect depth was evaluated using the following procedure.
1. Temperature contrast curves were generated from the experimental data, using the same non-defective temperature rise computed at a small area located on top of the SMA wire, and were smoothed locally using simple moving averaging. 2. The defect depth was then evaluated at each location along the wire employing the expression for d 1 in equation (2.12) where the minimum temperature T c min and time t min are computed from the temperature contrast curve of the experimental data.
There will be small positive or negative values of contrast temperature for all non-defective regions compared to the reference one, because of local differences in the heat conduction phenomena and noise in the thermographic system. These can be safely discarded in three ways: (a) by setting to zero the defect depth when the contrast temperature is found to be positive; (b) by setting to zero the defect depth when the contrast temperature is found to be negative, but is far away from the location of the defect (which is assumed to have been previously estimated from the thermal images), and In fact, by construction, expression (equation (2.12)) will produce a real positive, physically meaningful result for the defect depth only for values of contrast temperature and time that fit closely the curve in Fig. 5. 3. Finally, previously discarded results for the defect depth were reconstructed using inpainting techniques, which aim at filling-in holes in digital data by propagating surrounding data (e.g. [31]), the only condition being that the defect depth on the boundaries of the plate should be zero.
In Figs. 5, 6 and 7 the defect diameter D is assumed to have been estimated precisely from the thermal images. In real applications a component would carry a grid of SMA wires, at known distances and depths, that can be used to estimate the extent of the damage, hence the value of the diameter to be used in the expression for d 1 in equation (2.12). A methodology for estimating the size of the internal damage based on the a priori knowledge of the inter-wire distance and length is discussed in [20].
Despite good agreement with the experimental data (errors ranging between 1-7 %), some limitation must be highlighted. The analytical expression for the defect depth assumes known material and defect thermal properties, geometry of the sample and defect diameter and opening. While the defect diameter can be estimated from the thermal images, it is not possible to estimate the defect opening. All that is feasible is to evaluate the defect depths for a range of defect openings for a specific application. Furthermore, minimum contrast temperature and time using the approximate expression were not found to be very sensitive to large variations of the defect depth, but highly dependent on defect opening.

Conclusions
This paper developed an analytical model to derive an approximate expression for evaluating the defect depth for pulsed thermography performed using embedded SMA wires in laminated composites. The model was compared against experimental data and showed a good agreement in predicting the general form of the contrast curves for images of defects and for predicting the time and magnitude of the peak in contrast. In addition, the results for the defect depth computed showed vey good agreement with the experimental data (errors ranging between 1-7 %). This paper showed that SMart thermography and the developed analytical model could be used as a quantitative nondestructive tool where defect size and depth could be estimated with good accuracy.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.