Sideband Peak Count in a Vibro-Acoustic Modulation Method for Crack Detection

: This paper presents a new method of signal processing for vibro-acoustic modulation (VAM) methods in order to detect damage accumulation in steel samples. Damage in the tested samples was produced by cycle loading, which, with a small amplitude, was used as a pump wave to modulate an ultrasonic probe wave. Multiple sideband peaks were observed, which were used to characterize the modulation effect. We propose the effectiveness sideband peak number (SPN) method as an indicator of any damage accumulation when the load cycle is applied. Moreover, after comparing the SPN with the previously used modulation index (MI), we concluded that, for some of the samples, the SPN provided better damage indication than the MI. The presented results can be explained by a simple model of bilinear crack nonlinearity. This model demonstrates that the amplitude dependences of the sideband components on the pump and the probe wave amplitudes are very different from the quadratic crack model that is usually used for MI test explanation.

Among the numerous NA approaches, one of the most investigated is the vibroacoustic modulation (VAM) technique that is also referred to as nonlinear wave modulation spectroscopy and dynamic acoustic-elastic testing (DAET) [9]. This method is based on the modulation of high-frequency probe acoustic waves by low-frequency vibration pump waves. One explanation of this modulation is that a low-frequency wave changes the sound speed in the tested material, leading to a change in the amplitude and phase modulation of the probe wave. This phenomenon of sound speed change in a solid material under applied pressure is called acoustoelasticity. The first paper on this subject was published in 1958 [10]. Toupin and Bernstein later presented the theory of this effect on solid material [11].
The two kinds of approaches to VAM differ by way of low-frequency vibration generation. The first approach is that vibrations are excited by harmonic sources. The first paper on this matter was [12], while one of the newest papers is [13]. The other sources of harmonic vibration excitation are piezoceramic transducers [14,15], speakers [16,17], lasers [18], and loading machines [19]. The second and simpler VAM approach of lowfrequency excitation by applying mechanical impact was first suggested by Korotkov and Sutin [20] (see also [21][22][23]). This impact excitation does not require additional equipment and can avoid additional nonlinearity produced by contact between low-frequency sources and the tested material. The majority of VAM techniques consider the level of sideband components with sum and difference frequencies as an indicator of nonlinearity. In many publications [24][25][26], the ratio of sideband peaks to the probe wave amplitude is used as a Acoustics 2022, 4 75 modulation indicator. This approach is based on the assumption of the quadratic type of nonlinearity due to a crack presence.
In many VAM experiments, the generation of multiple sideband components was observed. These sideband components have frequencies f ± nF, where f is the probe frequency, F is the pump frequency, and n is the integer number. Hu et al. [14] demonstrated the generation of multiple sideband peaks in VAM experiments in order to detect fatigue cracks in aluminum. The authors did not count the sideband peaks but used Hilbert-Huang transforms to demodulate the received signal. The generation of multiple sideband peaks in the VAM method was presented in [27]. This paper investigated the VAM application for delamination detection in glass-fiber-reinforced polymer laminates, and sideband peak amplitudes were used as features in a machine learning classification algorithm.
In the work presented here, we investigated multiple sideband peak generation using cycled loads for the damage accumulation and for modulation of the ultrasonic wave, according to how it was made in [19]. As a result, there was an increase in the number of sideband peaks due to damage accumulation, allowing the peak numbers to be utilized as damage indicators. The result was compared with the more commonly used modulation index (MI).
In previously conducted VAM experiments [19,[24][25][26][28][29][30][31], the MI was calculated on the basis of first sideband peaks with frequencies f ± F and numerous sideband peaks with n > 1 were not taken into account. The goal of our work is to investigate multiple sideband peak generation and to consider whether these peaks can provide more information about damage accumulation for the same A108 samples that were investigated in [19,29,30].
It is worth mentioning that a slightly different sideband peak count method was used earlier [8,28,32,33]. In this method, the excitation of the samples was conducted with multifrequency signals, and frequency components generated due to crack nonlinearity were counted. These results were explained by the quadratic nonlinearity that is not applicable to our experiments. Instead, the multiple sideband peak generation that we observed can be explained by the simple bilinear model, which predicts the unusual dependence of sideband peaks on the amplitudes of the probe and pump waves.

Configuration of the Experimental Setup
The experimental scheme is the same as that described in [19], which is illustrated in Figure 1. The tested samples of our VAM experiment were 10 inch A108 steel bars with a cross-section 1 × 1 /8 inch cross-section. The same samples were used in earlier tests [19,29,30] that allowed a comparison of our SPN method with the previous MI method. The sample mechanical properties are illustrated in Table 1. At the center of each specimen, a circular notch with a 1 /4 inch diameter was drilled to control the location of the expected fatigue crack. The specimen was installed in a material test system (MTS), and two piezoelectric discs were attached to the samples as ultrasonic transmitters and signal receivers. Two processes, namely, damaging and measuring, were performed alternatively on the tested samples until each sample came to a complete break. During the degradative process, the cyclic tensile loads, with a target setpoint of 10.5 kN and amplitude of 10 kN, were applied to the sample by the MTS. The lifetime of the six tested samples varied from 45,000 to 120,000 cycles of such fatigue loads.
The MTS was utilized during the VAM measuring process to generate the pump wave with a frequency of 10 Hz and an amplitude of 0.6 kN. The probe wave was generated by an NI USB 6361 data acquisition (DAQ) board amplified 50 times by a piezo driver power amplifier and was sent to the sample by one of the piezoelectric discs illustrated in Figure 1. Thus, in the presence of damage in the sample, the probe and pump wave interact with each other at the notch area.
Another piezoceramic disc, which received a modulated signal, went through a bandpass filter with a passband of 110-220 kHz and was sent back to the DAQ board. Finally, the digital signal was sent to the computer, where spectral signal analysis was conducted.
Because of the dependence of the modulation on the probe frequency wave, we conducted tests for five probe waves frequencies. Since single frequency use does not provide reliable results, the frequency response is acquired by a sweeping frequency test of the sample at its intact stage with the pump wave applied [31,34]. Several criteria need to be considered when choosing the probe wave frequency. Firstly, the magnitude of the carrier should be strong enough on an excellent signal-to-noise ratio. Secondly, the magnitude of the sidebands in the intact sample should be lower than −50 dB for the Two processes, namely, damaging and measuring, were performed alternatively on the tested samples until each sample came to a complete break. During the degradative process, the cyclic tensile loads, with a target setpoint of 10.5 kN and amplitude of 10 kN, were applied to the sample by the MTS. The lifetime of the six tested samples varied from 45,000 to 120,000 cycles of such fatigue loads.
The MTS was utilized during the VAM measuring process to generate the pump wave with a frequency of 10 Hz and an amplitude of 0.6 kN. The probe wave was generated by an NI USB 6361 data acquisition (DAQ) board amplified 50 times by a piezo driver power amplifier and was sent to the sample by one of the piezoelectric discs illustrated in Figure 1. Thus, in the presence of damage in the sample, the probe and pump wave interact with each other at the notch area.
Another piezoceramic disc, which received a modulated signal, went through a bandpass filter with a passband of 110-220 kHz and was sent back to the DAQ board. Finally, the digital signal was sent to the computer, where spectral signal analysis was conducted.
Because of the dependence of the modulation on the probe frequency wave, we conducted tests for five probe waves frequencies. Since single frequency use does not provide reliable results, the frequency response is acquired by a sweeping frequency test of the sample at its intact stage with the pump wave applied [31,34]. Several criteria need to be considered when choosing the probe wave frequency. Firstly, the magnitude of the carrier should be strong enough on an excellent signal-to-noise ratio. Secondly, the magnitude of the sidebands in the intact sample should be lower than −50 dB for the minimum influence of structural nonlinearity and maximum sensitivity of the VAM. Thirdly, in order to avoid the node or the antinode of the modes of the sample, the frequency should be selected from peaks or troughs of the frequency response curve.

Results
The MI was computed according to the VAM technique described in [35] in order to evaluate the remaining fatigue life using VAM signals. The MI was usually calculated from the level of probe frequency (A 0 ) and both side components A + and A − having frequencies f + F and f − F, where f is the probe frequency and F is the pump frequency. The VAM signals were recorded at every cycle point for five probe frequencies in the range of 150 to 210 kHz and a pump frequency of 10 Hz. The recorded signals were studied in the frequency domain by applying the frequency Fourier transform (FFT).
Examples of the signal spectra for Sample #3 for the frequencies around the probe frequency f = 175 kHz are presented in Figure 2. These spectra are shown for different numbers of cycles, including the level of sideband peaks that were observed with cycle number increase.
frequency should be selected from peaks or troughs of the frequency response curve.

Results
The MI was computed according to the VAM technique described in [35] in order to evaluate the remaining fatigue life using VAM signals. The MI was usually calculated from the level of probe frequency ( 0 ) and both side components + and − having frequencies + and − , where is the probe frequency and is the pump frequency. The VAM signals were recorded at every cycle point for five probe frequencies in the range of 150 to 210 kHz and a pump frequency of 10 Hz. The recorded signals were studied in the frequency domain by applying the frequency Fourier transform (FFT).
Examples of the signal spectra for Sample #3 for the frequencies around the probe frequency f = 175 kHz are presented in Figure 2. These spectra are shown for different numbers of cycles, including the level of sideband peaks that were observed with cycle number increase. The number of cycles to failure was 52,106 (see Table 2).
Representative results of the MI were averaged for five tested probe frequencies. The tested samples had a different number of cycles to failure and we used this number for normalization of load cycles in the figures below. This normalized load cycle number can be considered as the sample lifetime.
The MI dependences for three samples are shown in Figure 3.  Table 2). The nonlinear effects in VAM are usually described by the following equation: Representative results of the MI were averaged for five tested probe frequencies. The tested samples had a different number of cycles to failure and we used this number for normalization of load cycles in the figures below. This normalized load cycle number can be considered as the sample lifetime.
The MI dependences for three samples are shown in Figure 3.
There are three stages apparent in the diagram of the MI measurement that are depicted by a change of inclination in the trend. At the first stage, the MI starts with a slight difference between −60 dB and −45 dB with increasing cyclic load. At this stage, there is no visible crack. The second stage is the trending line before growing, where invisible damage occurs in the samples. We refer to this stage as "initial damage accumulation". Figure 2a,b display the spectra of Sample #3. These spectra of MI signals remain unchanged until the cycle load ratio achieves an approximate 0.9 MI increase that is evidenced in the spectrum in Figure 2c due to the presence of a crack. This MI increase is the third stage, where the damage becomes visible until the end of the measurement. There are three stages apparent in the diagram of the MI measurement that are depicted by a change of inclination in the trend. At the first stage, the MI starts with a slight difference between −60 dB and −45 dB with increasing cyclic load. At this stage, there is no visible crack. The second stage is the trending line before growing, where invisible damage occurs in the samples. We refer to this stage as "initial damage accumulation". Figure 2a,b display the spectra of Sample #3. These spectra of MI signals remain unchanged until the cycle load ratio achieves an approximate 0.9 MI increase that is evidenced in the spectrum in Figure 2c due to the presence of a crack. This MI increase is the third stage, where the damage becomes visible until the end of the measurement.   Table 2).
The sideband peaks were counted for frequencies f + nF, where n is an integer number. Red dots mark the counted peaks for these frequencies. Sideband peak numbers (SPNs) were calculated on the right side of the sideband because the spectrum is symmetrical regarding the probe frequency f. As a result, the SPN expanded due to damage accumulation, as seen in Figure 5.   There are three stages apparent in the diagram of the MI measurement that are depicted by a change of inclination in the trend. At the first stage, the MI starts with a slight difference between −60 dB and −45 dB with increasing cyclic load. At this stage, there is no visible crack. The second stage is the trending line before growing, where invisible damage occurs in the samples. We refer to this stage as "initial damage accumulation". Figure 2a,b display the spectra of Sample #3. These spectra of MI signals remain unchanged until the cycle load ratio achieves an approximate 0.9 MI increase that is evidenced in the spectrum in Figure 2c due to the presence of a crack. This MI increase is the third stage, where the damage becomes visible until the end of the measurement.   Table 2).
The sideband peaks were counted for frequencies f + nF, where n is an integer number. Red dots mark the counted peaks for these frequencies. Sideband peak numbers (SPNs) were calculated on the right side of the sideband because the spectrum is symmetrical regarding the probe frequency f. As a result, the SPN expanded due to damage accumulation, as seen in Figure 5.  Table 2).
The sideband peaks were counted for frequencies f + nF, where n is an integer number. Red dots mark the counted peaks for these frequencies. Sideband peak numbers (SPNs) were calculated on the right side of the sideband because the spectrum is symmetrical regarding the probe frequency f. As a result, the SPN expanded due to damage accumulation, as seen in Figure 5.
Notably, from a comparison of Figures 3 and 5, both techniques, namely, MI and SPN, reflect damage accumulation in a similar way. Hence, both parameters can be used for damage and crack assessment. Acoustics 2022, 4 FOR PEER REVIEW 6 Figure 5. SPN versus the number of normalized cycles.
Notably, from a comparison of Figures 3 and 5, both techniques, namely, MI and SPN, reflect damage accumulation in a similar way. Hence, both parameters can be used for damage and crack assessment.
However, unlike Samples #1-3, the tests with Samples #4-6 showed an MI and SPN behavior difference. The examples of the modulation spectra in Figure 6 show that the level of sideband peaks was not increased with damage accumulation, but the number of sideband peaks was increased, as seen in Figure 7. However, unlike Samples #1-3, the tests with Samples #4-6 showed an MI and SPN behavior difference. The examples of the modulation spectra in Figure 6 show that the level of sideband peaks was not increased with damage accumulation, but the number of sideband peaks was increased, as seen in Figure 7. Notably, from a comparison of Figures 3 and 5, both techniques, namely, MI and SPN, reflect damage accumulation in a similar way. Hence, both parameters can be used for damage and crack assessment.
However, unlike Samples #1-3, the tests with Samples #4-6 showed an MI and SPN behavior difference. The examples of the modulation spectra in Figure 6 show that the level of sideband peaks was not increased with damage accumulation, but the number of sideband peaks was increased, as seen in Figure 7.  Notably, from a comparison of Figures 3 and 5, both techniques, namely, MI and SPN, reflect damage accumulation in a similar way. Hence, both parameters can be used for damage and crack assessment.
However, unlike Samples #1-3, the tests with Samples #4-6 showed an MI and SPN behavior difference. The examples of the modulation spectra in Figure 6 show that the level of sideband peaks was not increased with damage accumulation, but the number of sideband peaks was increased, as seen in Figure 7.  In the tests of these samples, the SPN demonstrated much better performance for damage detection than the MI. The dependence of MI for this case in Figure 8 does not display a significant improvement compared to the MI method with crack accumulation.
Acoustics 2022, 4 FOR PEER REVIEW 7 In the tests of these samples, the SPN demonstrated much better performance for damage detection than the MI. The dependence of MI for this case in Figure 8 does not display a significant improvement compared to the MI method with crack accumulation.   Another advantage of an SPN calculation compared to MI is that the SPN is less dependent on the choice of probe frequency. Figure 10 shows the MI and SPN measured for Sample #3 for all five tested frequencies. It is clearly seen that the SPN has significantly fewer fluctuations for different frequencies compared to the MI.  Figure 9 shows an evaluation of Samples #4-6 using the SPN method. The result presents a noticeable trend toward an increase during the fatigue test. After 80% of their lifetime, a significant increase was observed in the sideband values. The results of a robustness test confirm that the SPN method can provide a consistent measure of minor damage to the specimens.
Acoustics 2022, 4 FOR PEER REVIEW 7 In the tests of these samples, the SPN demonstrated much better performance for damage detection than the MI. The dependence of MI for this case in Figure 8 does not display a significant improvement compared to the MI method with crack accumulation.   Another advantage of an SPN calculation compared to MI is that the SPN is less dependent on the choice of probe frequency. Figure 10 shows the MI and SPN measured for Sample #3 for all five tested frequencies. It is clearly seen that the SPN has significantly fewer fluctuations for different frequencies compared to the MI. Another advantage of an SPN calculation compared to MI is that the SPN is less dependent on the choice of probe frequency. Figure 10 shows the MI and SPN measured for Sample #3 for all five tested frequencies. It is clearly seen that the SPN has significantly fewer fluctuations for different frequencies compared to the MI. Acoustics 2022, 4 FOR PEER REVIEW 8 Figure 10. Variations of parameter MI (upper) and SPN (lower) for Sample #4 for different probe frequencies.

Bilinear Crack Model of Multiple Sideband Peak Interpretation
As mentioned in the introduction, various VAM research studies have observed the presence of multiple sideband peaks. These multiple sideband peaks can be explained by complex (nonquadratic) crack nonlinearity and/or higher-order interactions [36].
A much simpler and more reliable explanation can be made on the basis of a bilinear crack model. Many publications consider the nonlinear acoustic effect in the structures with closing (breathing) cracks, and some examples are presented in [3,[37][38][39]. The bilinear model in [40] displayed multiple sideband peaks in a cracked beam due to the simulated effect of stochastic vibrations [41].
In this work, we considered the generation of sideband peaks according to the simplest crack model. Figure 11 shows the force-displacement spring relationship of a bilinear model in a section medium. Thus, the crack is closed by positive applied force (compression) and opened by negative applied force (extension), which increases the stiffness coefficient.

Bilinear Crack Model of Multiple Sideband Peak Interpretation
As mentioned in the introduction, various VAM research studies have observed the presence of multiple sideband peaks. These multiple sideband peaks can be explained by complex (nonquadratic) crack nonlinearity and/or higher-order interactions [36].
A much simpler and more reliable explanation can be made on the basis of a bilinear crack model. Many publications consider the nonlinear acoustic effect in the structures with closing (breathing) cracks, and some examples are presented in [3,[37][38][39]. The bilinear model in [40] displayed multiple sideband peaks in a cracked beam due to the simulated effect of stochastic vibrations [41].
In this work, we considered the generation of sideband peaks according to the simplest crack model. Figure 11 shows the force-displacement spring relationship of a bilinear model in a section medium. Thus, the crack is closed by positive applied force (compression) and opened by negative applied force (extension), which increases the stiffness coefficient. Acoustics 2022, 4 FOR PEER REVIEW 9 Figure 11. Schematic illustration of opening-closing action of cracks in beams.
The relationship of the normalized displacement on the normalized applied force is shown in Figure 12. Some essential features of the multifrequency sideband peak generation can be explained by the simple bilinear crack model, which has various stiffnesses for positive and negative force.
Consider excitation frequencies, which are much lower than resonance frequencies, and small nonlinearity k(q) = K0 + Dk(q), where the nonlinear part of the stiffness coefficient Equation (2) can then be written in the following form: For the bilinear case, we consider the variation of the stiffness as follows: (4) Figure 12 shows the case for = 0.5. This strong nonlinearity is only used for the visibility of the stiffness variation. For demonstration of multiple sideband peak generation, = 0.1 was used. Illustration of the sideband components for both models was conducted for the case where the input stress is the sum of the probe wave and the pump wave with the same amplitudes P1 and P2. Figure 11. Schematic illustration of opening-closing action of cracks in beams.
The relationship of the normalized displacement on the normalized applied force is shown in Figure 12.
Acoustics 2022, 4 FOR PEER REVIEW 9 Figure 11. Schematic illustration of opening-closing action of cracks in beams.
The relationship of the normalized displacement on the normalized applied force is shown in Figure 12. Some essential features of the multifrequency sideband peak generation can be explained by the simple bilinear crack model, which has various stiffnesses for positive and negative force.
Consider excitation frequencies, which are much lower than resonance frequencies, and small nonlinearity k(q) = K0 + Dk(q), where the nonlinear part of the stiffness coefficient Equation (2) can then be written in the following form: For the bilinear case, we consider the variation of the stiffness as follows: (4) Figure 12 shows the case for = 0.5. This strong nonlinearity is only used for the visibility of the stiffness variation. For demonstration of multiple sideband peak generation, = 0.1 was used. Illustration of the sideband components for both models was conducted for the case where the input stress is the sum of the probe wave and the pump wave with the same amplitudes P1 and P2. Some essential features of the multifrequency sideband peak generation can be explained by the simple bilinear crack model, which has various stiffnesses for positive and negative force. m ..
Consider excitation frequencies, which are much lower than resonance frequencies, and small nonlinearity k(q) = K 0 + Dk(q), where the nonlinear part of the stiffness coefficient Dk(q) << K 0 .
Equation (2) can then be written in the following form: For the bilinear case, we consider the variation of the stiffness as follows: (4) Figure 12 shows the case for α = 0.5. This strong nonlinearity is only used for the visibility of the stiffness variation. For demonstration of multiple sideband peak generation, α = 0.1 was used.
Illustration of the sideband components for both models was conducted for the case where the input stress is the sum of the probe wave and the pump wave with the same amplitudes P 1 and P 2 . σ = P 1 cos ωt + P 2 cos Ωt, where the ω = 2π f and Ω = 2πF. Figure 13 represents the sideband component for the bilinear and quadratic model, calculated using the following parameters and values: frequency of the pump wave, 10 Hz; frequency of the probe wave, 200 kHz; nonlinear parameter, a = 0.1; P 1 = P 2 = 1.
Acoustics 2022, 4 FOR PEER REVIEW 10 where the = 2 and = 2 . Figure 13 represents the sideband component for the bilinear and quadratic model, calculated using the following parameters and values: frequency of the pump wave, 10 Hz; frequency of the probe wave, 200 kHz; nonlinear parameter, a = 0.1; P1 = P2 = 1.
(a) (b) Figure 13. Simulation spectrum (a) of the bilinear and (b) quadratic crack models.
The dependences of the first sideband components on the probe and pump wave amplitudes are presented in Figure 14. The solid curved line represents the bilinear model in both figures, while the straight dashed line represents the quadratic model. Interestingly, both dependences are practically the same and reveal a difference with regard to the quadratic model. Figure 15 reveals examples of the spectra calculated from the bilinear crack models. A significant difference in variation of sideband components was observed for highorder components. Figure 15 shows examples of the calculated spectra. The dependences of the first sideband components on the probe and pump wave amplitudes are presented in Figure 14. The solid curved line represents the bilinear model in both figures, while the straight dashed line represents the quadratic model. Interestingly, both dependences are practically the same and reveal a difference with regard to the quadratic model. Figure 15 reveals examples of the spectra calculated from the bilinear crack models.
Acoustics 2022, 4 FOR PEER REVIEW 10 where the = 2 and = 2 . Figure 13 represents the sideband component for the bilinear and quadratic model, calculated using the following parameters and values: frequency of the pump wave, 10 Hz; frequency of the probe wave, 200 kHz; nonlinear parameter, a = 0.1; P1 = P2 = 1.
(a) (b) Figure 13. Simulation spectrum (a) of the bilinear and (b) quadratic crack models.
The dependences of the first sideband components on the probe and pump wave amplitudes are presented in Figure 14. The solid curved line represents the bilinear model in both figures, while the straight dashed line represents the quadratic model. Interestingly, both dependences are practically the same and reveal a difference with regard to the quadratic model. Figure 15 reveals examples of the spectra calculated from the bilinear crack models. A significant difference in variation of sideband components was observed for highorder components. Figure 15 shows examples of the calculated spectra. Analysis of the simple bilinear crack model shows that this model can explain the multiple sideband peak generation. In addition, this model demonstrates that the strong dependence of the sideband peaks on the probe and pump wave amplitudes could be a topic of future experimental investigations.

Conclusions
In this study, a new method based on sideband peak count was suggested for estimating the nonlinear effects of the vibro-acoustic modulation technique of crack detection. This method allows for the consideration of the multiple peaks generated in VAM due to nonquadratic crack nonlinearities. An experimental comparison between the earlier developed modulation index and the new sideband peak number parameters demonstrates that, in some of the tests, the latter reflects damage more effectively than the classical MI. A simple bilinear crack model was applied for an interpretation of multiple sideband peak generation This model indicates the effect of the unusual dependence of the sideband peaks on the probe and pump wave amplitudes, which is a potential topic for future experimental investigations. Funding: The work of Abdullah Alnutayfat was financially supported by Prince Sattam Bin Abdulaziz University, Saudi Arabia. sudden death in June 2021 shocked and saddened the Stevens community. Sophia, a dedicated professor in the Department of Civil, Environmental, and Ocean Engineering at the Stevens Institute of Technology, was passionate about research and teaching. She was honored with numerous awards and recognitions for her teaching and research. Sophia initiated this work and led this research until her death. The memory of Sophia will always be in our hearts.

Conflicts of Interest:
The authors declare no conflicts of interest. A significant difference in variation of sideband components was observed for highorder components. Figure 15 shows examples of the calculated spectra.
Analysis of the simple bilinear crack model shows that this model can explain the multiple sideband peak generation. In addition, this model demonstrates that the strong dependence of the sideband peaks on the probe and pump wave amplitudes could be a topic of future experimental investigations.

Conclusions
In this study, a new method based on sideband peak count was suggested for estimating the nonlinear effects of the vibro-acoustic modulation technique of crack detection. This method allows for the consideration of the multiple peaks generated in VAM due to nonquadratic crack nonlinearities. An experimental comparison between the earlier developed modulation index and the new sideband peak number parameters demonstrates that, in some of the tests, the latter reflects damage more effectively than the classical MI. A simple bilinear crack model was applied for an interpretation of multiple sideband peak generation This model indicates the effect of the unusual dependence of the sideband peaks on the probe and pump wave amplitudes, which is a potential topic for future experimental investigations.