Abstract

This research proposes a correlation coefficient for detecting and evaluating defects in beams, which brings about a positive outcome in terms of accuracy and efficiency. This parameter surpasses other parameters, such as natural frequency and damping coefficient, thanks to its sensitivity to structural changes. Our results show that although the damping coefficient had more variation than the natural frequency value in the same experiment, its changes were insufficient and unstable at different levels of defects. In addition, the proposed correlation coefficient parameter has a linear characteristic and always changes significantly according to increasing levels of defects. The results outweigh damping coefficient and natural frequency values. Furthermore, this value is always sensitive to measurement channels, which could be an important factor in locating defects in beams. The testing index is statistically evaluated by a normal distribution of the amplitude value of vibration measurement signals. Changes and shifts in this distribution are the basis for evaluating beam defects. Thus, the suggested parameter is a reliable alternative for assessing the defects of a structure.

1. Introduction

Mathematically, research on detecting defects in bridge beams and other construction components has taken place in two major areas. The first study uses mathematical models to simulate the structures of mechanical systems. Researchers do not focus on characteristics of defects but on their general specifications in terms of kinematic and dynamic responses. The defects caused by changes in geometric dimensions, geometric shapes, or materials might include cracks, such as horizontal cracks [14] or vertical cracks [5, 6], geometric shapes, [7, 8] such as circles, ovals, or obliques, a reduction of mechanical characteristics, such as elastic, torsional, or bending modulus [911] or a combination of the above factors. These studies use certain parameters to evaluate the overall stiffness level, EJx. The second research area, which seeks to obtain a solution to the inverse problem, is aimed at determining structural characteristics by using algorithms to retrieve relevant information. Some of these algorithms also evaluate the performance of a structure under actual working conditions. In order to increase the accuracy of parameters collected from different measurement signals, researchers have developed different mathematical methods, for example, statistical methods [12, 13] or combination with optimization algorithms [1416], algorithms of artificial neural networks [1719], wavelet transformation [20, 21], or neuron-fuzz [15, 22]. Figure 1 shows the steps involved in detecting and forecasting defects in a structure.

The procedure goes from the lowest to highest level, in detail:(i)Detecting: to inspect and confirm the presence of defects(ii)Locating: to localize defects(iii)Evaluating: to assess the impact of these defects on the structure and determine their degrees in order to organize periodic tests and inspection(iv)Forecasting: to promptly forecast defects, minimizing loss and damage to personnel and assets

As detecting is very important, it has attracted the interest of researchers all over the world. Researchers have mostly opted for relatable signals, such as vibration measurement signals, deformation measurement signals, and transposition signals (deflection).

Some outstanding examples using vibration measurement signals are found in [23, 24]. Parameters related to these signals, such as natural frequencies, mode shapes, vibration amplitudes, and damping coefficients, were commonly used in [25, 26]; among these, natural frequency was one of the most popular. Many of the above studies have shown that natural frequency does not express changes in the overall stiffness of mechanical systems very well, as it is not sensitive enough in practice. For example, in [27], when the beam had three cracks that accounted for 50% of its thickness, the natural frequency was reduced by less than 5%. Salawu [28] measured the status of a concrete bridge before and after the trimmer beam was replaced and realized that the frequency of the first six segments changed by 1.7% on average. The results from Salawu show that natural frequency can only detect relatively large defects but cannot locate them due to the same changes in frequency [29]. Other researchers, such as Loland [30], Fox [31], Meneghetti and Maggiore [32], and Choy et al. [33], also had the same viewpoint [3438]. Therefore, we must discover new characteristics that are sensitive enough to exploit.

This article proposes a new characteristic built on actual vibration signals and, specifically, on the correlation coefficient. Changes in this parameter are the basis for evaluating the existence and development of defects or locating defects inside structures. Using test beam models, this study not only determines the natural frequencies and damping coefficients of mechanical systems suffering from different extents of defects but also compares their changes to the new proposed parameter using mathematical statistics. In addition, it also defines the normal distribution of vibration measurement signals through experimentation. Most importantly, this new parameter can be applied to complex structural models.

2. Theoretical Basis

2.1. Theoretical Vibration of Beam Model

Any vibration in a mechanical system is considered to be a linear combination of basic vibrations. In this investigation, a bar type is used with a basic load-bearing state of compression and tension at the center, bending, torsion, or a combination of different kinds of vibration. However, every vibration can be composed of each type of basic vibration. They are usually represented by the harmonic vibration formula:in which ωNi is the ith natural frequency. Vectors following the order of Wi are called ith private vectors and are assigned so that

These values depend on the geometric dimensions of bridge span, material characteristics, and sectional structure. Technically, natural frequency has a special feature: when the frequency of a constraining force is equal to the natural frequency value, the amplitude of the corresponding response will sharply increase due to resonance. Determining the natural frequency of basic vibrations of elastic beams applies the formulas below.(i)Natural frequency of the bending vibration of a simply supported beam:in which l is the length of the structure, EJz is its bending stiffness, A is its sectional area, and ρ is its net weight.(ii)Natural frequency of torsional vibration:in which GJz is the convolution stiffness and J0 is the area of the cross-section.(iii)Natural frequency of torsional-bending vibration:in which , , and .

When , the torsional-bending frequency (5) is determined byin which and .

Equations (3), (4), and (6) show that if the mechanical system vibrates, the impact of impeding factors will make the amplitude decrease over time. The vibration energy in the system is gradually lost until it vanishes (becomes equal to 0). The decreasing energy directly affects the overall stiffness of the beams. However, the natural frequency value rarely changes with respect to the stiffness reduction of structures, which means the energy loss is unable to change the natural frequency value. To evaluate the energy loss mentioned above, researchers focus on examining the damping coefficient. When a mechanical system vibrates under bending, torsion, or a torsional-bending state, the vibration equation is

When a mechanical system suffers from damping vibrations, the natural frequency is reduced. To evaluate the new vibration state (commonly called vibration damping), studies have examined the damping frequency of mechanical systems. ωd is the damped natural frequency. Therefore,

To determine the damping coefficient ζ, experimental methods record the decreasing vibration graphs of the system. q1 and q2 represent two adjacent peaks of the decreasing vibration graph (Figure 2). At the recorded time, the two peaks are t1 and t2. We thus have a ratio:

As q1 and q2 are two adjacent peaks of the graph, t2=t1+T d:

Noting equation (10), when applying the damped vibration equation to equation (8), we get the following equation:

Equation (9) then becomes equation (12). This is the basis for determining energy loss through vibration amplitude:

The natural logarithm of the above equation is called the logarithmic decrement and is symbolized as δ. This coefficient expresses the loss of energy during vibration as

Provided that , the final expression becomes .

If , then the above expression can be written as .

In order to increase the accuracy of determining the ζ value, instead of two adjacent peaks, we propose to take two peaks at m intervals. Thus, the ratio between q1 and qm is

We combine equation (14) with equation (13) to determine the damping coefficient under all circumstances:

However, the results of calculating the damping coefficient ζ usually contain errors. During actual experimentation, we have identified two common cases:(a)The calculated results have higher values than the reality. This case usually occurs when the test beam is overloaded at a low velocity. In this circumstance, the deflection of the beam suddenly increases to maximum, while the vibration decreases quickly.(b)The calculated results have smaller values than the reality. This case usually occurs when the beam suffers from forced vibration, with its frequency approaching its natural frequency. The load and velocity on the beam are sufficient to make the bridge vibrate drastically. At that point, the forced vibration lasts longer than the decreasing vibration. As such, the test beam vibrates harshly, and forced vibration in beams creates large vibration amplitudes. This phenomenon causes a deterioration in the bridge capability and destroys the structure.

Therefore, calculating the actual damping value is relatively difficult with remaining deviations. This is why we propose a new parameter to evaluate the energy loss of structural vibration. This parameter must simplify both calculation and measurement.

2.2. Proposal of Correlation Coefficient Model

Correlation coefficients are commonly used in cases of studying the relation between two parameter sets with the same number of samples X and Y. There are different types of calculations based on usage. The most popular correlation coefficient is the Pearson correlation coefficient, which is calculated by dividing the covariances of two variables by their multiplied standard deviations. This calculation is shown asin which E is the operator of expected value and Cov is the covariance. Another formula is widely applied according to equation (16):

As , and so is Y.

On the contrary,

Therefore, equation (16) is rewritten as

The correlation coefficient shall only be defined as above if the standard deviation is bounded and not equal to zero. Mathematically, the correlation coefficient is equal to one in the case of covariant linear correlation and −1 in the case of contravariant linear correlation. Other values ranging from [−1, 1] show the level of linear dependence between variables. The closer the correlation coefficients are to −1 and 1, the stronger the correlation between variables. However, the values in mechanical fields only range between [0, 1]. It can be concluded that the correlation coefficient is equal to 1 when no energy loss occurs and equal to 0 when no energy exists after vibration transmission.

3. Experimental Model

The experimental model was created at the Laboratory of Applied Mechanics (LAM), Ho Chi Minh City University of Technology (Figure 3), utilizing three groups of equipment: a model group, excitation-creating group, and signal-receiving group.

The model group includes a steel plate with a 0.9 m length, 0.01 m width, and 0.005 m height (see Table 1) to simulate the load-bearing beam. The steel beam rests on two trimmer beams, as shown in Figure 4(a). Additionally, a deceleration part is installed at each head of the test beam to make the model more realistic. When a vehicle approaches, passes through the beam, and finishes its movement, it does not make the steel beam become tight. To simulate a complicated force acting on the test beam, the study installed one more thin homogeneous plate that creates eccentric vibrations. Thus, this experiment can change not only the running velocity of vehicles but also the rotating velocity of vehicle engines (Figure 4(b)).

The excitation-creating group consists of one engine driven by a belt-drive system. This system lets the vehicle run over the beam in one direction at a steady velocity, as shown in Figure 5(a). When the vehicle moves, the velocity is controlled by two inverter machines that can change the rotating velocity of the engine and the running velocity of the vehicles (Figure 5(b)).

The signal-receiving group measures various kinds of parameters, including the velocity, acceleration, deformation, and transposition-meter measurement signals. The transposition meter is directly attached to the middle point of the beam (Figure 6(a)). For vibration signals, parameters from velocity and acceleration sensors were collected. There are four acceleration sensors and four velocity sensors attached to the beam at four locations. The measurement points are distributed evenly on the beam, as shown in Figure 6(b). At each point, we arranged one velocity sensor and one acceleration sensor to simultaneously measure both signals. Finally, a strain gauge was installed at the center of the beam, as shown in Figure 6(c).

The experiment evaluated six different velocities: V50 = 18.84 cm/s, V60 = 25.12 cm/s, V70 = 28.26 cm/s, V80 = 31.4 cm/s, V90 = 35.3 cm/s, and V100 = 39.2 cm/s. Each velocity responded to eight different levels of defects created in the beam (according to Table 2). We treated cracks as defects. At each stage of the experiment, the cracks were deepened, until the beam was nearly destroyed. The moving loads from left abutments to right abutments on the beams and sample parameters are shown in Table 3.

4. Results and Discussion

4.1. Stability of Natural Frequency

A power spectrum for the test beam was created, ranging from nondefective to different levels of defects. The results, shown in Figure 7, include Figure 7(a): a free vibration of the intact beam, Figure 7(b): a free vibration of the beam with a 0.4 mm deep and 0.3 mm wide crack, Figure 7(c): a free vibration of the beam with a 0.8 mm deep and 0.3 mm wide crack, Figure 7(d): a free vibration of the beam with a 1.2 mm deep and 0.3 mm wide crack, Figure 7(e): a free vibration of the beam with a 1.6 mm deep and 0.3 mm wide crack, Figure 7(f): a free vibration of the beam with a 2.4 mm deep and 0.3 mm wide crack, Figure 7(g): a free vibration of the beam with a 2.6 mm deep and 0.3 mm wide crack, and Figure 7(h): a free vibration of the beam with two defects of 0.3 mm width wherein the first crack is 2.6 mm deep and the second crack is 0.8 mm deep.

Changes in the natural frequency value were too small (not exceeding 4%) during the experiment. They were not adequate for determining the existence of defects in the beam. Furthermore, under free vibration, not only one natural frequency but also others nearby on the spectrum appeared. These frequencies had almost identical values (with only a 0.05 Hz differential), and the amplitudes of different spectra corresponded to different levels of cracks. The natural frequency tended to continuously reduce together with the levels of cracks. However, this reduction was too small in comparison with the changes in the defects. From the intact level (beam with no crack) to the maximum level of defect (H26), the value of the frequency was reduced by no more than 0.5 Hz. Therefore, changes in natural frequency cannot evaluate, locate, or determine the level of defects in the beam.

Figures 8(a)8(h) show the frequency values of the power spectra under forced action when the vehicle model passes through the beam. These values correspond to the levels of defects (as shown in Table 2) at different velocities (V50 = 18.84 cm/s, V60 = 25.12 cm/s, V70 = 28.26 cm/s, V80 = 31.4 cm/s, V90 = 35.3 cm/s, and V100 = 39.2 cm/s). When the vehicle moved on the beam, together with a basic natural frequency (determined when the beam vibrates freely), there were more forced frequencies, such as f1 = 6.2 Hz, f2 = 7.8 Hz, f3 = 9.4 Hz, and f4 = 10.8 Hz, and higher frequencies beyond the natural frequency value. Corresponding to different levels of defects, the natural frequency remains nearly unchanged, which is similar to the free vibration of beams. This result accords with prior research [27, 28, 39, 40]. Most researchers conclude that changes in natural frequency are not sensitive enough to determine a partial breakdown of the test or actual beams. The problem is the accuracy of measurement equipment: it is unable to spot these small differences. Furthermore, frequency is an overall quantity: it is not sensitive to partial changes. Thus, natural frequency does not provide enough information on structural changes.

4.2. Changes of Damping Coefficient

This research also studied changes in damping components at different levels of defects. In general, damping coefficients collected from vibrations of sensors always vary in a similar way. The results show that(i)From defect level H0 to level H12, the damping coefficient tended to increase and reached a maximum value of 0.3. This change was relatively larger and more distinct than natural frequency values, as shown in Figures 914. However, under different circumstances of defects, such as H16, the value was reduced to 0.15, lower than that of the intact beam (with no defect). Further investigation showed that this coefficient increased slightly until level H26, where the beam was entirely destroyed. Moreover, there was no difference in damping coefficients at each measurement point (sensor location) according to Table 3.(ii)The velocity of the experiment did not significantly affect the value of the damping. Figures 914 show the damping coefficient values at six different velocities corresponding to each level of defects.

In general, the change in the damping coefficient is insufficient and unstable. Therefore, we can hardly recognize the existence of defects, especially in beams suffering from major defects. In addition, the changes in the damping coefficients at four measurement points were relatively the same, which made it difficult to locate cracks in the beam. Although it performs better than natural frequency, the damping coefficient is still an impractical, low-sensitivity parameter. This conclusion is accepted in other studies [25, 26], in which the damping coefficient value is affected by the mechanical characteristics of the material. The structural changes in the experimental models only have impacts on small parts, not the whole beam structure.

4.3. Changes in Correlation Coefficient
4.3.1. Changes in Correlation Coefficient

Due to insufficient damping coefficients and natural frequencies, we propose to use correlation coefficients. Both Figure 15 and Table 2 show the relation between two measurements (Figure 4(a)) channels through correlation coefficients. The new proposed parameter has always shown a significant linear change corresponding to the increase in defect degree. According to Figure 15, when correlation coefficients decreased, the level of cracks was found to have increased. They decreased within the range of [1, 0]. If the correlation coefficient is equal to 1, there is perfect correlation. This usually happens to correlation coefficients of measurements at the same point (autocorrelation). On the contrary, when the correlation coefficient between two points is equal to 0, there is no binding relation. In terms of energy distribution in beams, the reduction of correlation coefficients implies energy loss at each level of defects. The correlation coefficient is influenced by velocities during the experiment. This means that the velocity in the experimental model scatters the vibrating energy. It changes the fluctuating values of consecutive elements. As defined in the draft, the more the value of the velocity increased, the more the value of the correlation coefficient decreased. In the manuscript, we want to define the relationship between these two parameters (correlation coefficient value and experimental velocity). Hence, the existence of defects is the main reason for this loss. In each case, we have the following comments:(i)Figure 15 shows changes in correlation coefficients between channel 1 and the remaining channels. The value of channel 1 had significant changes compared with other channels (100%), from a nondefective level to defect level H26. This change is the rule and remains steady throughout the extent of defects. It also surpasses the results of the damping coefficient and natural frequency values. This research will delve into the effective usage of this parameter at different levels of defects.(ii)Figures 1618 show changes in correlation values in channel 2 compared with the remaining channels as well as in channel 3 with the remaining channels. The results indicate the symmetric characteristics of the two measurement channels throughout different levels. In other words, the distribution of measurement channels in the model can affect this value during evaluation. Therefore, this value is sensitive to the measurement channel and crucial for locating defects in the beam. Similar conclusions can be made upon investigating the correlation values of channels R1-4 (between channel 1 and channel 4) and channels R4-1 (between channel 4 and channel 1).(iii)When the velocity changed, the correlation coefficients sped up rapidly (Figures 1518). This suggests that the relation between correlation coefficients and velocity corresponds linearly to each level of defects. The impact of velocity on this parameter increased the feasibility of detecting defects, as it can show the overall condition of a structure under acting force.

4.3.2. The Relation between Correlation Coefficients and Defects

Corresponding to different levels of defects, correlation coefficients are distributed in pairs as follows:(i)Regarding the value of correlation coefficients between two consecutive channels (R1-2 and R3-4) without cracks, Figure 19 shows almost no connection with the levels of cracks in the beam. When the cracks were at the highest level (H26) and V100 velocity, the changes in channel R1-2 were approximately 40% in comparison with the intact beam. The change was nearly 38% in channel R3-4. The distribution of the correlation coefficient rule in these similar channels is equivalent and relatively linear at each level of cracks. This is the standard input for evaluating the presence of cracks in beams.(ii)Regarding the correlation coefficient values of two consecutive channels (R2-3) with cracks (Figure 20), the correlation value changed quickly according to the increase in the cracks. This change was significant (over 90%) between the intact beams after they reached the level (H26) at V100 velocity. In the first (R1-2 and R3-4) and second (R2-3) circumstances with the same level of defects, the changes were sufficient for determining the location of the cracks in the beam.(iii)Regarding the correlation coefficients of two channels that were far away from each other (R1-3 and R2-4) (see Figure 21) with defects in between, their values among these channels were relatively different. Under forced vibration, the vibration energy was transmitted from channel 1 to channel 2, passing through the defect and finally to channel 3. As analyzed above, the energy was partially lost during transmission from one point to another according to the direction of vibration. It experienced maximum loss when passing by the defect location. However, the correlation coefficients of (R1-3) had lower energy loss because the distance between channel 1 and the defect location was large. When it passed the defect location, the remaining energy was not large enough to result in such a significant loss. Contrary to (R1-3), the correlation value of (R2-4) clearly increased at each level of defects. This meant that the distance between the measurement point of channel 2 and the defect location was shorter, causing more significant losses. We can immediately notice the impact of defects on the correlation coefficients.

Figure 22 shows the changes in the correlation coefficient values at different levels of defects in the beam. This is evidence that affirms the sensitivity of this parameter to defects. These changes, compared with the changes in natural frequency values and damping coefficients, are relatively large. Therefore, they support the application of this coefficient for detecting and evaluating defects in structures.

4.3.3. The Impact of Velocity on Correlation Coefficients and the Development of Defects

Velocity has a significant impact on changes in correlation coefficient values at different levels of defects in beams. In order to simplify the evaluation process, this study only investigated a case of two consecutive channels with the presence of (R2-3) or the absence of (R1-2 and R3-4) defects. The results show that(i)For intact beams (beams without defects), as shown in Figure 23(a), or beams with some defects (H4, H8, H12), as in Figure 23(b), the impact of velocity was insignificant. The defect was not enough to reduce the correlation coefficients. The impact of velocity shows that the energy consumed when passing through the defect location is small.(ii)On the contrary, at the high scale of defects, the impact of velocity was more significant. Figures 23(e)23(h) show that the higher the velocity, the more energy in the beam. However, the higher the levels of defects, the more energy is lost. Thus, there is a large difference in energy loss between intact beams and defective beams, which is an important point to evaluate in the study of defects in beams.

5. Conclusion

In this research, we proposed a new parameter—the correlation coefficient—to detect and evaluate defects in beams. It is much more sensitive to defects than the natural frequency or the damping coefficient. This research has raised four prominent points:(i)The natural frequency values in the experimental model were not sensitive enough to detect defects in the beam. The changes in this value were too small compared to the development of defects. Therefore, natural frequency measurements hardly meet the requirements of detecting defects in structures.(ii)Although the damping coefficient had more variation than the natural frequency value in the same experiment, its changes were also insufficient and unstable at different levels of defects. Therefore, it is difficult to apply the damping coefficient method to detect the existence of defects, especially in beams suffering from major defects. In addition, changes in the damping coefficients at four measurement channels were relatively the same, so this characteristic makes it hard to locate cracks in the beam. Thus, the damping coefficient is still a low-sensitivity parameter that is difficult to apply.(iii)The new proposed parameter has a linear characteristic and always significantly changes according to increasing levels of defects. The results show that the measurability of the changes using this method outweighs those from damping coefficient and natural frequency values. Furthermore, this value is always sensitive to measurement channels, which could be an important factor in locating defects in beams.(iv)However, this research has limited value for investigating the relation between correlation coefficients and changes in beam stiffness due to cracks. Future research needs to be conducted to look into the quantitative evaluation of correlation coefficients to clarify that relation in the monitoring and identification of defects in structures.

Data Availability

All data generated or analysed during this study are included in this published article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors wish to thank various people for their contribution towards this work. Mr. Lam Q. Tran and Mr Vien N. Le contributed for technical support on this project. Prof. Nhi K. Ngo, staff of Applied Mechanics (LAM), Ho Chi Minh City University of Technology, contributed in collecting the plant data and all the technicians helped in handling the instruments. The authors would like to specially thank Prof. Nhi K. Ngo, research project supervisor, for her professional guidance and valuable support and constructive recommendations on this project.