Experimental and Numerical Study on Characterization and Evaluation of Surface Cracks with Wireless Ultrasonic Sensor

Adopting both wireless ultrasonic sensing and numerical simulation techniques, this research investigates the interaction between Rayleigh wave and artificial surface cracks of varying depths. When analyzing experimental ultrasonic data collected by a wireless sensing node, the signals are enhanced through a two-step procedure including signal reconstruction and envelope extraction.)e waveforms are interpreted in detail by analyzing wave components through time-of-flight technique. A finite element (FE) model is devised to properly simulate the experimental testing. )e simulated waveforms are consistent with experimental results and corroborate the analysis and explanations of experimental waveforms. Based on both experimental and numerical waveform analysis, a relationship between ultrasonic characteristic parameter and crack size is established for the quantitative estimation purpose. )e proposed model shows a good agreement with data from both test and literatures.


Introduction
e service life of major civil engineering structures can range from decades to hundreds of years. During this period, adverse effects, such as environmental erosion, material aging, and long-term loading effect, can lead to accumulation of damage and structural deterioration [1][2][3]. Under extreme circumstances, unexpected failure of aging and deteriorating infrastructures can threat public safety. In order to prevent such catastrophes, real-time information on structural condition should be collected. To this end, numerous structural sensing techniques have been developed and studied.
Among various techniques, nondestructive testing (NDT) has been widely used in practice. NDT technique can reveal both internal and surface defects without introducing further damage to the structure and usually with little interruption to the normal operation. Correlated with changes of different physical phenomena (e.g., pressure, heat, and scattering of elastic waves), NDT can provide valuable data for evaluating structural integrity, safety, and reliability.
Among different NDT methods, Rayleigh wave-based ultrasonic testing method is particularly suited for detecting the presence, orientation, distribution, and geometric shape of surface cracks. A primary advantage of the Rayleigh wave is that during propagation it suffers less spatial attenuation than body waves. Most energy of Rayleigh wave propagates within a wavelength from the free surface of a tested specimen [4]. Moreover, Rayleigh wave is sensitive to the surface defects with a depth less than a few times of the wavelength [5]. e Rayleigh wave velocity is constant in the same medium, and the wavelength is only related to frequency. erefore, the detection depth can be easily controlled by varying the frequency of the wave. Because of all the features mentioned above, Rayleigh wave is ideal for nondestructive evaluation with low-power wireless sensing devices.
However, because the interaction mechanism between Rayleigh wave and crack is not yet fully understood, the detection and sizing of surface cracks with Rayleigh wave is still a research subject of great interest. A number of experimental research studies on this topic have been reported.
Viktorov [4] observed a scattering phenomenon of Rayleigh wave at a notch in an aluminum bar and proposed a crack size estimation method involving reflection and transmission coefficients. After calculation, he also pointed out that there was a loss of complementarity between the two coefficients due to the conversion and energy loss at the interface near the notch. On the basis of Viktorov's work, Domarkas et al. [6] studied the scattering from surface cracks at a fixed incidence angle with Rayleigh waves of different frequencies.
ey explained the variation of reflection coefficient as resonances at the crack and used it for crack dimension estimation with a satisfactory precision. Rokhlin and Kim [7] conducted a detailed experimental study of scattering of Rayleigh waves at pit-induced cracks in an AI2024-3 steel plate. ey studied the reflected waveforms under the influences of crack depths and tension stepup loads, compared the experimental data with theoretical results, and explored the potential application for crack sizing in a fatigue test. Hernandez-Valle et al. [8] adopted laser-induced Rayleigh wave for studying branched defects in aluminum samples. Branch was added to a slot to simulate stress corrosion defects, and a connection was successfully established between characteristics of collected signals and the geometry of cracks by amplitude and frequency spectrum analysis.
In parallel to experimental study on the interaction between Rayleigh wave and surface defects, a great number of articles have also been published on aspects of theoretical modeling and numerical simulations. Comparing experimental studies, theoretical and numerical research can provide detailed motion prediction of all particles in the body, providing guidance towards improving future experiments and data analysis. For example, Kino [9] and Auld [10] derived mathematical formula by reciprocity theorem, which could be used to calculate the transmission coefficient from one transducer to another through an arbitrarily shaped flaw. Achenbach and Gautesen [11] extended geometrical theory of diffraction from the optics research to investigate the diffraction of elastic waves by a crack. ey study detailed the propagation process from the source to the receiver through a crack and successfully formulated the characteristic interference pattern of diffracted wave using both frequency and angle. Valle et al. [12] established a finite element (FE) model to simulate the propagation of guided circular wave in a cracked hollow cylinder. ey showed the feasibility to accurately size and locate a radial crack by analyzing simulated results using Auld's formula and timefrequency representations. Rosli et al. [13] performed another study on detecting surface-breaking defects of a wide range of depths and angles. ey built a 2D FE model to simulate Rayleigh wave interactions and revealed the trends that transmission coefficients vary with different geometric parameters of defects, which agreed well with experimental results.
Traditional cabled equipment for Rayleigh wave testing is usually bulky, expensive, and inconvenient to deploy in the field. e equipment is better suited for short-term inspection than long-term deployment in the field. As reported by Straser and Kiremidjian [14], for a general bridge testing system, cable installation can take up to 75% of total time and 25% of total cost. On the contrary, wireless smart sensors with attractive characteristics, such as compactness, ability of wireless data transmission, low cost, power efficiency, and on-board processing capabilities [15], can be a promising alternative to traditional cabled equipment. e compact size and wireless communication make deployment much more convenient in the field. e low-cost attribute makes it feasible and economical to deploy a dense array of sensors on a civil structure, significantly increasing the capability for accurate damage localization. On-board processing feature enables computation of damage-sensitive parameters in real time, immediately upon data collection. Designed with these characteristics, Martlet is a wireless sensing platform codeveloped by the University of Michigan, Georgia Institute of Technology and Michigan Technological University [16]. Martlet adopts a dual-core Texas Instruments processor with a hardware floating-point unit, which supports high-speed data acquisition at up to 3 MHz sample rate. It also incorporates an extensible platform to support heterogeneous sensing for both global-and localbased damage detection. A customized daughter board was designed for ultrasonic measurement and integrated with the Martlet mother board to compose a wireless ultrasonic sensing node. Architectural design and a preliminary performance validation of the wireless ultrasonic sensing node was illustrated thoroughly in a preceding article [17].
is paper aims to investigate the interaction of Rayleigh wave with artificial surface cracks of varying depths, through both wireless sensing experiments and FE analyses. e rest of the paper is organized as follows. Section 2 illustrates the experimental setup for the ultrasonic measurement using a wireless sensing node. In Section 3, proper selection of key parameters (e.g., material properties, excitation wave generation, meshing size, time step, and boundary condition) is discussed through finite element modeling and correlating simulation with experimental results. Section 4 provides analysis and comparison of Rayleigh wave scattering from the perspective of waveform and transmission coefficients. e study will help improve quantitative size estimation of surface defects and ultrasonic waveform analysis in future fatigue crack detection. Section 5 presents the summary and conclusions.

Experimental Setup
e wireless ultrasonic sensing node consists of a Martlet mother board and an ultrasonic daughter board that can be plugged-in on top ( Figure 1). A detailed description of the hardware architecture has been previously reported [17] and thus will not be repeated herein. In short, the ultrasonic daughter board contains modules for both ultrasonic signal excitation and signal receiving. e wireless ultrasonic sensing node can perform signal conditioning to the received ultrasonic signal prior to digitization. e excitation module generates a series of square wave bursts at 500 kHz with peak-to-peak amplitude up to 40 V (an amplitude of 18 V is used in this study). In the receiving module, a band-pass filter and an amplifier are used to increase the signal-to-noise ratio of the received ultrasonic signal. Considering the frequency of interest is around 500 kHz, the high-pass and low-pass cutoff frequencies are set at 159 Hz and 2 MHz, respectively, to suppress background noise while avoiding waveform distortion. In addition, the gain for output signal amplification is selectable among 10 dB, 20 dB, and 30 dB using a rotary switch on the ultrasonic daughter board (30 dB by default in the experiment). e accuracy of the generation and recording of ultrasonic waveforms are verified by comparing to those generated/received by cabled commercial equipment. Figure 2 presents the thick steel specimen plate used in this study, which is 241.3 mm × 152.4 mm × 25.4 mm (thickness). Four artificial cracks with different depths (0.51, 1.27, 2.29, and 3.05 mm) were created using an electrical discharge machine to simulate surface cracks with different levels of damage (depths). e artificial cracks have the same length of 9.3 mm and opening width of 0.58 mm. e distances between adjacent artificial cracks are about 25 mm. e experimental setup is shown schematically in Figure 3. e ultrasonic measurements were performed using a pair of ultrasonic transducers (Model: WC50-0.5, Ultran Group) with a central frequency of 500 kHz. Each transducer was attached to a Plexiglas wedge for Rayleigh wave generation and detection separately, and the distance between the source and receiver is 25.4 mm in a typical pitch-catch arrangement. Light lubrication oil was used as acoustic couplant at the interfaces between transducer/wedge and wedge/specimen surfaces. When the experiment begins, the server wakes up the wireless ultrasonic sensing node that commands the ultrasonic daughter board. e ultrasonic daughter board then generates five cycles of 18 V and 500 kHz square wave at every burst and repeats the same bursts every 260 μs. e bursts are fed into the transmitting ultrasonic transducer to launch Rayleigh waves in the specimen. Consequently, the propagating waves are detected by the receiving ultrasonic transducer, and the ultrasonic signals are collected by the signal receiving module in the ultrasonic daughter board. After filtering and amplifying, analog output signals are digitized at a sampling rate 2 MHz. e digitized data of the received ultrasonic signal can be wirelessly transferred to the server. Using the amplitude of the received Rayleigh wave, which propagated through the artificial crack, it is expected to distinguish the depth of each artificial crack. Before reviewing the experimental results, the next section first describes finite element modeling that attempts to simulate this ultrasonic testing.

Finite Element Modelling
Rayleigh wave has perpendicular and parallel components, which indicates it is a two-dimensional propagation wave [4]. In the material with characteristics of homogeneous, isotropic, and linearly elastic (e.g., steel in this paper), the interaction of plane Rayleigh wave in a steel specimen with an artificial crack can be analyzed using a two-dimensional model [18], as presented in Figure 4. e established FE model attempts to simulate the experimental testing by characterizing the interaction between Rayleigh wave and artificial cracks under the plane-strain condition. To fulfill this end, PLANE42 element is adopted in commercial  software ANSYS, which is a plane element and defined by four nodes each with translation along both x and z directions. As shown in Figure 4, artificial cracks have a rectangular shape with different depths and their surfaces have stress-free boundary conditions. e height of the FE model is the same as the specimen thickness (25.2 mm).
Considering both be consistent with the experiment and convenience for modeling, the length of the model is taken to be 100.8 mm. e coordinate origin of the model is set at the bottom center, and the coordinate orientations are shown in Figure 4 as well.

Material Properties and Excitation Signal Generation.
Since the material properties have a significant influence on the simulation results [19], they should be chosen as close as possible to the actual properties of the specimen. According to the handbook published by American Society of Metals [20], the density, elastic modulus, and Poisson's ratio of simulated steel material are 7.85 × 10 3 kg/m 3 , 200 GPa, and 0.24, respectively. Moreover, considering reference data provided in literature [21] and a measurement carried out with a specimen of similar material [22], the propagation velocities of compression, shear, and Rayleigh waves in the model are selected as 5900 m/s, 3200 m/s, and 2989 m/s, respectively. e excitation signal is generated at the position labeled as "source" in Figure 4, and the surface distance from the source to the artificial crack center is 12.7 mm. Two approaches are frequently used to shape waveforms of excitation signals. One is to approximately apply a time-varying force waveform at the "source" FE node based on a certain function (e.g. half-sine function) [23]. e other is to specify the displacement waveform of the "source" node and generate a displacement record closely matching the motion equation of Rayleigh wave in homogeneous elastic materials. e actual motion that occurs in an experimental testing is subject to the electromechanical coupling and sensitivity coefficients of the particular ultrasonic transducer used in the experiment. To identify an appropriate waveform that models the motion of "source" node, Treeby et al. [24] applied a semisinusoidal window to a series of sinusoidal signals to generate a toneburst. ey compared the simulated wave with the theoretical solution in both time and frequency domains and demonstrated the high accuracy. In this paper, the second approach is adopted to generate the excitation signal, and the waveform of a five-cycle 500 kHz tone-burst with normalized amplitude is presented in Figure 5.

Finite Element Meshing.
Meshing size in finite element modeling can significantly influence simulation accuracy, which is affected by the wavelength and frequency of interest. According to Mirahmadi and Honarvar [25], the maximum meshing size should be less than 1/10 of the shortest propagating wavelength. For a 500 kHz Rayleigh wave, the smallest wavelength is approximately 5.98 mm; thus, the maximum element size should be less than 0.598 mm. In addition, Valliappan and Murti [26] provided another empirical equation to estimate the element size g: where ζ is an empiric constant related to material properties, which takes a value of 0.25 for homogeneous materials; V s is the shear wave velocity; and f is the maximum Rayleigh wave frequency of interest. By equation (1) and material properties discussed in Section 3.1, the maximum element size is calculated and should be smaller than 1.6 mm. Considering the calculations above, the element sizes along x-and z-axis are set at 0.3 mm and 0.4 mm in this study. As a result, the undamaged FE model has 21,424 elements and 21,888 nodes, and the number is reduced by 2-18 for a model with an artificial crack of varying depth.

Boundary Conditions.
In the FE model, the nodes at the specimen bottom are fixed in both vertical and horizontal directions. Due to sufficient spacing between the artificial cracks and the boundaries of the test specimen, it should experience little interference from boundaries or a neighboring artificial crack. erefore, the wave can be modeled as propagating indefinitely along both positive and negative directions of x-axis. To achieve this effect in the FE model, an artificial boundary is set to offset the impact caused by the left and right boundaries. In this paper, artificial boundaries with viscous springs are adopted along the x-axis as the absorbing boundary condition ("Absorbing B. C." in Figure 4) for its simplicity, stability, and effectiveness [27]. By adding linear spring and damper elements along both horizontal and vertical directions at the nodes on an absorbing boundary, most of the propagating energy is absorbed to minimize reflection from these boundaries. e stiffness coefficient for the springs and the damping coefficient for the dampers are calculated as follows [28]: where K x and K z are stiffness coefficients along x and z directions, separately; α x and α z denote empirical artificial boundary parameters, which are set at 1.0 and 0.5 according to the recommendation [29]; G represents shear modulus of the material; R is the distance between the source and the boundary; C x and C z are damping coefficients along x-and zaxis, respectively; ρ denotes the medium density; and V p is the compression wave velocity.

Simulation Time
Step. e selection of time step should take into account both numerical accuracy and computing time. An empirical formula was given by Valliappan and Murti [26]: where τ represents the time step and g is the element size calculated from equation (1). As a result, the time step should be smaller than 0.27 μs. In order to be comparable with experimental results, the time step is set at 0.1 μs in our simulation. Based on the experimental data [30], the collected signal fades away in about 90 μs, and the total number of simulated time steps is determined as 900.

Data Processing Approaches for the Wireless Ultrasonic
Data. e sampling frequency used in the experiments was 2 MHz. Because the central frequency of the ultrasonic signal was at around 500 kHz, 4 data points were sampled for every cycle of the ultrasonic signal. Although the sampling frequency is relatively low, it meets the requirement of Nyquist-Shannon sampling theorem [31,32]. Signal reconstruction using upsampling and the cardinal sine function (a.k.a. sinc function) can be carried out to restore most of the details in the band-limited ultrasonic waveform [33]: where x re (t) is the reconstructed and upsampled signal; t represents time; x [i] is the i th sampled data point in the wireless ultrasonic data; n denotes the total number of data points; f s is the wireless sampling frequency. During the measurement, a slight dislocation of the transducer positions causes change in the received signal with a different waveform and different local amplitudes. e solid and dashed curves in Figure 6 illustrate two collected signals for which the receiving transducer position is slightly disturbed.
is phenomenon can affect peak amplitude determination, which, in turn, causes an error in artificial crack size estimation. In order to reduce the effect of transducer positions, an envelope detection process using discrete Hilbert transform is performed. e theoretical derivation of the algorithm was presented in a previous paper [17]. Generally, the envelope of an upsampled signal can be extracted using where x env re (t) is the corresponding envelope curve for a collected signal x re (t) and x re (t) denotes the Hilbert transform of x re (t), which can be obtained through a series of calculations involving convolution and Fourier transform. e solid curve in Figure 6 shows the envelopes of the two collected signals overlap.

Characterization of Wireless Ultrasonic Data.
e wireless data with a sampling frequency of 2 MHz were upsampled to 10 MHz with the algorithm proposed in the preceding section. e validation of the algorithm was conducted by a comparison between the upsampled wireless sensing signals and corresponding signals collected from commercial cabled equipment at same sampling frequency. Detailed analysis was reported in the earlier paper [17]. Figures 7(a)-7(e) show the ultrasonic waveforms and corresponding envelope curves, for each of the undamaged case and four artificial crack cases. e waveform fluctuates very little after 70 μs, and for ease of comparison between experimental and simulation waveforms in the following section, the gain of the collected signals (30 dB) is removed.
Using the propagation velocities of compression, shear, and Rayleigh waves discussed in Section 3.1, the propagation distance difference between waveforms of different arrival times can be calculated. Based on geometric analysis, the causes of waveforms at different arrival times can then be deduced. erefore, relative arrival times of different wave components are crucial in this study, while absolute arrival time is irrelevant. For comparison, in Figures 7(a)-7(e), all waveforms are time-aligned so that the first arriving peak with a significant positive amplitude is at around 10 μs.
In these figures, a very similar waveform trend can be observed. e first group of waves arrives approximately between 8 and 20 μs, which means the propagation distance is the shortest. e first waves packet contains the highest peak, and its amplitude decreases significantly for larger artificial crack depths. Since the amplitude has a positive correlation with energy, the first arriving component has the highest percentage of the total energy. Considering this component has the features of shortest propagation distance and the highest share of energy, based on Pertsch et al. [34], the first arrival can be attributed to transmitted Rayleigh wave, which travels along the horizontal surface and the edge of the artificial crack. It is marked as path #1 in Figure 7(f ), and its propagation distance is the straight-line distance between two sensors, which is 25.4 mm. e second group of waves arrives between approximately 20 and 30 μs. Based on the arriving time difference from path #1, the additional propagation distance should be about 29.89 mm (for Rayleigh wave) or 32 mm (for shear wave). erefore, the total propagation distance of the second group of waves is between 55.29 mm and 57.4 mm. Based on geometry relationship, the distance for a wave traveling from the source to the receiver through a single bottom reflection (path #2 in Figure 7(f )) is 56.80 mm, which confirms the second group to be bottom reflection waves. In addition, its amplitude is low and the arriving time slightly lengthens as the artificial crack depth increases. ese features match the description of creeping waves reported in literature [35]; thus the waves propagating through paths #3 and #4 in Figure 7(f ) may also contribute to the formation of second group of waves. e third group of waves arrives at approximately 30 μs. Following the procedures of time-difference analysis used previously, the traveling distance for these waves is around 100.13 mm. Considering the distance for one propagation is around 56 mm, the third group of waves is likely the double reflection shear waves (path #5 in Figure 7(f )). In addition, Pertsch [18] pointed out that a wave with a delay of 25 μs exists in the excitation signal, due to the reflection inside the wedge with identical dimensions. erefore, this "delayed" excitation signal travels through path #1 and should reach the receiver at around 35 μs. erefore, it can also contribute to the formation of third group of waves. e full route of this delayed wave should be a combination of paths #6 and #1 in Figure 7(f ).
Last but not least, the fourth group of waves arrives at approximately 45 μs, and the arriving time varies with the depth of the artificial crack. Based on the different arrival times, the propagation distances are roughly between 130 mm and 160 mm and match the geometric distance of a triple reflection. Since this group of waves helps little in the artificial crack size estimation, the origin of the fourtharriving waveform will not be discussed in detail in this paper.

FEM Simulation of the Interaction between Surface
Acoustic Wave and a Crack. A proper FE model should not only accurately reflect experimental conditions, but also provide simulation results that are close to experimental data. To validate the FE simulation, both undamaged and damaged scenarios are adopted. Figure 8   considering correlation coefficients are close to 1 and the amplitude of all waveforms is on the order of 0.01 V (much higher than the RMSE values), it can be concluded that the simulations well modeled the experiments. e selected waveform of the excitation signal is also deemed appropriate to simulate the ultrasonic source in the experiment. Hence, the simulation methods established in this study can be used for future analysis of wave interaction, for example, identifying the origin of waveforms with different arrival times.
If the waveforms in the time domain are examined carefully, they can be divided into multiple segments, from which several common features can be observed by comparing the waveforms. Firstly, it can be found that simulation waveforms agree well with the experimental results between 10 and 20 μs. is waveform segment is easily attributed to the Rayleigh wave propagation, based on previous time-of-flight analysis. In the meantime, it can also be observed that there is an approximately 1 μs delay between simulated waveforms and the experimental ones during the 20-40 μs (see dashed boxes in Figure 8). is phenomenon indicates that the velocity of the simulated wave may have a slight difference from the experiment. In the literature [36,37], it is generally acknowledged that wave velocity is only associated with material properties in a homogeneous linear elastic medium that is strain-/stress-free. Despite the careful selection of material properties, we can reasonably suspect that the minor time difference between simulated and experimental waveform is caused by a small difference between simulated and experimental material properties. Considering the similarity in analysis, scenario of cracks depths at 0.40 mm is selected as the representative for the study of the origin of wave segments with different arrival times.
At the 14 th μs in Figure 8(b), the peak of the transmitted Rayleigh wave passes through the receiver for the 0.40 mm artificial crack. Figure 9(a) shows the simulated wave propagation amplitude at this instant. In the figure, four significant waves can be observed, that is, incident wave, transmitted wave, bottom reflected wave, and artificial crack reflected wave. Among them, the propagation directions of incident and transmitted waves are opposite and are perpendicular to a virtual plane defined by x � − 12.6 mm (source position). e wave amplitudes are nearly symmetric about this plane, except that the peak amplitude of the transmitted wave is a little lower than that of the incident wave, which is due to the energy loss of the transmitted wave at the artificial crack. In addition, the transmitted wave shows a significant attenuation inside the specimen, which agrees with the feature of a Rayleigh wave. erefore, the simulation corroborates the conclusion of experimental analysis, in which the first wave group mainly consists of transmitted Rayleigh waves. In addition, the figure shows a bulk wave propagates toward the bottom, gets reflected by the bottom surface, and finally arrives back at the top surface.
is bottom reflected wave arrives at the receiver approximately between 20 μs and 30 μs, and its propagation feature matches the description of shear wave in Liu et al. [38]. Consequently, the bottom reflecting shear wave in Figure 9(a) supports the previous time-of-flight analysis of the experimental data collected by the receiver between 20 μs and 30 μs. Figure 9(b) shows another example instant from the same simulation model as in Figure 9(a). e instant corresponds to the 30 th μs in Figure 8(b), when the incident wave has propagated across the left boundary. From the figure it can be seen that the incident wave fades out and barely leaves any reflected waves. erefore, it can be concluded that the absorbing B. C. set in this study is functional and can help to minimize the reflection at the boundary. Furthermore, a bottom-surfacebottom reflected wave is shown in Figure 9(b) to reach the receiver around 35 μs after the data collection. e existence of this bottom double reflected wave in the simulation explains the experimental analysis of the wave between 30 and 45 μs in Figure 8(b).

Transmission Coefficients versus Crack
Size. Based on the waveform analysis above, it can be concluded that the firstarrival wave (i.e., transmitted Rayleigh wave) contains most of the energy in a collected signal and can be used to analyze the variation in artificial crack depth. Hence, based on the envelope amplitude of the transmitted Rayleigh wave, a transmission coefficient (T r ) is adopted for artificial crack size assessment and defined as where A env is the envelope amplitude of the transmitted Rayleigh wave from an artificial crack and A env undamaged is the envelope amplitude of propagated Rayleigh wave collected at the same distance from an undamaged area. e envelope amplitudes are obtained through discrete Hilbert transform as described in Section 4.1. Figure 10(a) illustrates the compare the experimental transmission coefficients versus normalized artificial crack depths (crack depth/wavelength of the 500 kHz Rayleigh wave) with these available in literature [4,[39][40][41][42]. Overall trend from these six studies is shown to match each other. Furthermore, the plot calculated from the wireless sensing data (in this study) is within the bounds provided by the references and close to the ones presented in Hévin et al. [41]. Using the wireless sensing data, a linear regression is performed between normalized crack depth and transmission coefficient, and the result is presented in Figure 10(b). e correlation coefficient and RMSE between experimental data and linear regression curve are 0.963 and 0.069, respectively. ese results show that using the transmission coefficient, the linear regression relationship can be used as a preliminary estimate for crack depth under similar testing conditions.

Conclusions
Aiming for the nondestructive quantitative assessment of crack sizes, this paper proposes a study on the interaction between Rayleigh wave and artificial surface cracks of different sizes using both experimental and numerical methods. A self-developed wireless ultrasonic sensing device is adopted for the ultrasonic measurement of surface cracks in a steel plate. In addition, for more effective identification purpose, a set of signal processing procedures is applied to the collected ultrasonic signals. rough wave propagation analysis, wave components with different arrival times are identified and explained. Emulating the experimental setup, an FE model is established. And a careful discussion is conducted on the key parameters, such as simulated material properties, method for Rayleigh wave excitation, finite element mesh size, boundary condition setup, and time step. e simulated waveforms are consistent with experimental results and corroborate the analysis of experimental waveforms. On the basis of analysis, transmission coefficient of the transmitted Rayleigh wave, which is sensitive to crack depth, is adopted to represent crack size. Furthermore, the linear regression results (between crack depth and transmission coefficient) have the potential for future quick estimation of crack depths under similar testing conditions.

Data Availability
All data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.