Detection of Bolt Looseness Based on Average Autocorrelation Function

Bolted connections are widely used in multiple engineering fields including aerospace and mechanical engineering due to their numerous advantages like the ability to bear relatively heavy loads, low costs, easy installation, and implementation. Bolt looseness may lead to costly disasters in industries and some cases of injuries. A new method for bolt looseness detection based on average autocorrelation function is proposed. It does not usually directly extract the looseness damage feature from the original response signal to detect the bolt looseness, but it uses the average autocorrelation function value at time lag T� 0 of the vibration pixel displacement signal to establish the looseness damage index. In terms of structural arrangement of this paper, firstly, the theoretical background of the proposed method is given. *en, an experimental system for bolt looseness detection based on computer vision is designed, and a verification experiment is carried out with the bolted connection plate as the experimental object. *e results show that the proposed method can effectively obtain the location of the looseness damage of the bolted connection plate, which provides a new technical reference for the online monitoring of the looseness damage of bolted connection plate.


Introduction
Bolted connections are widely used in multiple engineering fields including aerospace and mechanical engineering, because they have numerous advantages such as the ability to bear relatively heavy loads, low cost, and easy installation and implementation [1]. However, bolt looseness, induced by fatigue, impact, and thermal loads, may lead to costly disasters in industries and some cases of injuries [2,3]. erefore, for maintaining the healthy and stable operation of equipment, how to detect the bolt looseness simply and effectively as early as possible has attracted more and more attention by many researchers [4].
Nowadays, the methods of bolt looseness detection proposed by many researchers mainly include impedancemethod methods [5][6][7], vibration-based methods [8][9][10], and ultrasonic-based methods [11][12][13]. In the impedancebased methods, the piezoelectric transducers (PZT) are usually attached to the surface of a target structure being inspected [14]. Representatively, Lee [15] attempted a novel technique to detect bolt looseness using a transfer impedance technique for monitoring the structural health of a bolt joint. In addition, Sun and Zhang [16] proposed a monitoring method based on an impedance method that was researched for bolt looseness of steel structure. However, the impedance-based measurement method has high cost of the installation sensors. Moreover, in the vibration-based methods, representatively, Razi et al. [17] introduced a vibration-based health monitoring strategy for detecting the loosening of bolts in a pipeline's bolted flange joint. ey conducted numerical and experimental studies to verify the feasibility of the proposed method. Huda et al. [18] investigated a vibration testing and health monitoring system based on an impulse response excited by laser ablation. Furthermore, they also estimated the bolt loosening by detecting fluctuations of the high frequency response with their health monitoring system. In addition, Esmaeel et al. [19] explored a novel vibration-based damage detection methodology using the empirical mode decomposition (EMD) to calculate the energy damage index (EDI) in common industrial bolted joints. Additionally, in the ultrasonic-based methods, representatively, Zhang et al. [20] presented a bolt looseness recognition method based on the subharmonic resonance analysis. Furthermore, they also carried out analytical prediction and numerical simulations to verify the validity of the loosening detection method for bolted joint structures. Amerini and Meo [21] developed a tightening/loosening state index based on the first-order acoustic moment for the linear acoustic/ultrasound method. In addition, Fierro and Meo [22] proposed an in situ structural health monitoring approach based on the evaluation of linear and nonlinear modulated acoustic moments for the assessment of the loosened state of bolts. However, the ultrasonic-based method lacks accuracy due to the use of linear phenomena such as reflection and scattering [23]. e vibration-based method is simple in principle and easy to operate, and its measurement principles are mainly divided into two aspects: the first aspect is that the stiffness of the bolted connection is reduced by bolt looseness, which will result in the changes of the dynamic characteristic parameters (natural frequencies, damping ratio, and mode shapes) of the structure, and the second aspect is that the various signal processing methods are used to process the signals obtained from the sensors attached to the structure to extract the damage feature information of the structure. He and She [24] developed a new identification method for bolt looseness in wind turbine towers for preventing the collapse of wind turbine towers. ey found that the natural frequencies, damping ratio, and vibration characteristics are not very sensitive to the looseness of flange bolts, so they used the first-order phase difference curve to find the bolt looseness of the flanges. Vibration tests and analyses were carried out on six wind turbine towers with loose bolts, which verified their feasibility of the proposed method. Zhou et al. [25] formulated an energy-based damage index based on the high-frequency intrinsic mode functions. ey conducted a vibration test on the frame structure, and experimental results showed that the energy-based damage index can accurately detect bolt looseness damage. Chen et al. [26] proposed a looseness diagnosis method for connecting bolt of fan foundation based on sensitive mixeddomain features of excitation response and manifold learning. ey verified the feasibility of the proposed method with experimental results. Pnevmatikos et al. [27] selected the difference of the wavelet coefficients of the acceleration response of the healthy and loosened connection structure as an indicator of the damage and presented an application of wavelet analysis for damage detection of steel frames with bolted connections. e feasibility of the wavelet approach to damage detection was verified through experiments. He and Zhu [28] proposed using structural natural frequency changes to detect the loosening of bolted connections. ey took the steel pipe structure with bolted flanges as the experimental object and successfully detected the position of the loose bolts. Milanese et al. [29] explored two signal processing techniques to use in assessing the connection strength. e first method relied on basic statistical properties of the measured strains and their time derivatives, and the second method was based on power of the signal in different frequency bands. ey also verified the feasibility of the proposed methods by experimental studies.
In general, the vibration-based methods mentioned above in [28,29] can be used to identify bolt looseness. It is mainly based on the dynamic characteristic parameters or the signal processing methods to establish the damage index (DI). However, when the degree of bolt looseness damage is small, the change of the structural dynamic characteristic parameters is not very obvious. Further, the signal processing methods proposed in [9,25,26] can be usually implemented to detect the bolt looseness only when the input is confirmed. Actually, when the device is in operation, it means that the detailed information of the input is unknown, so how to recognize looseness based on ambient excitation is an important key affair. Since the excitation of the structure is unknown under operating conditions, taking into account the relevant theory of Natural Excitation Technique (NExT) in [30,31], the vibration response characteristics of the structure are analyzed under ambient excitation. Moreover, it is found that the autocorrelation function of the structural vibration displacement response is related to the structural modal parameters. Based on this reason, this paper proposes using the autocorrelation function of vibration response to establish damage index (DI) for detecting the bolt looseness. e remaining sections of this paper are as follows. e theoretical background of the response analysis under natural excitation, the definition of looseness damage index, and the proposed method with procedures are introduced in Section 2. e design principle of the experimental system and the experimental method are described in Section 3. e results are discussed in Section 4 to demonstrate the feasibility of the proposed method. At last, some conclusions are drawn in Section 5.

Response Analysis under Natural Excitation.
For a structural dynamics system with n degrees of freedom (DOFs), the dynamic equation describing the characteristics of the system can be expressed as [30,31] M € where M, C and K are the mass matrix, damping matrix, and stiffness matrix of the structure with the dimension of n × n, respectively. x (t) is the displacement vector of the structure response, and f (t) is the excitation vector of the structure. Because the response of the linear time-invariant system can be expressed as the product of the mode shapes and the modal coordinates, the structural response can be written as where Φ is the modal matrix composed of the modal vectors, q (t) is the vector of modal coordinates, q r (t) is the r th modal coordinate value, φ r is the r th mode shape, and n is the number of DOFs.
Since real normal modes are assumed, M, C, and K are diagonalized. Next, substituting equation (2) into equation (1) will result in the following [30].
where ω nr is the r th modal frequency, ζ r is the r th modal damping ratio, m r is the r th modal mass, and r � 1,. . ., n. Assume that, under the zero initial condition, we use Duhamel's integral to solve equation (3), and it will be obtained that Substituting equation (4) into equation (2) results in the following.
From equation (5), it can be obtained that if the excitation is f k (t) at point k on the system, the vibration response x ik (t) at point i can be expressed as where x ik (t) is the vibration displacement response at point i, f k (t) is the excitation at point k, φ ir is the r th mode shape value at point i, φ kr is the r th mode shape value at point k, and n is the number of DOFs. g r (t) is the impulse response function related to the r th modal parameters of the system. In addition, g r (t) is where m r , ζ r , and ω nr are the r th modal mass, modal damping, and modal frequency, respectively. ω dr is the r th damped modal frequency, . Based on the definition of autocorrelation function, the autocorrelation function R iik (T) of the vibration displacement response at point i with the time lag T is where E[·] is the expectation operator. Substituting equation (6) into equation (8) will obtain the following.
where φ ir is the r th mode shape value at point i, φ kr is the r th mode shape value at point k, φ is is the s th mode shape value at point i, and φ ks is the s th mode shape value at point k. n is the number of DOFs. g r (t) is the impulse response function related to the r th modal parameters. g s (t) is the impulse response function related to the s th modal parameters.
Assuming that the excitation f k (t) at point k is Gaussian white noise, the autocorrelation function of excitation can be expressed as where α k is the coefficient depending on exciting point k, and δ(t) is the pulse function. Substituting equation (10) into equation (9), and calculating the first integration by using the definition of δ(t), note that Because −∞ ≤ τ ≤ t, if we define λ � t − τ, the limits of λ will be 0 to ∞. en, the above equation (11) can be rewritten as Due to the definition of g r (t) from equation (7), equation (12) can be deformed into [30].
where A r iik and B r iik are functions of the modal parameters.
As described in Ref. [30], define the variable c rs as follows: tan c rs � I rs J rs . (16) Substitute equations (14) and (15) into equation (13), and simplify it to produce the autocorrelation function.
Note that the inner summation on s and k is just the summation of constants multiplied by the sine function. In addition, its frequency ω dr is fixed, and phase angel c rs is variable. erefore, the inner summation on s and k can be rewritten as a new sine function with a new phase angel θ r and a new constant multiplier G ir . where Based on equation (18), the above equation (17) can be further rewritten as where G ir is the coefficient depending on r th modal parameters and the response point i. θ r is the phase angle depending on r th modal parameters. Consequently, it can be found from equation (19) that the autocorrelation function of the structural vibration displacement response is related to the structural modal parameters.

Definition of Looseness Damage Index.
e structural looseness damage often changes the physical properties of the system. e looseness of bolted connections essentially causes the stiffness of the structure to decrease and ultimately leads to changes in the dynamic characteristics before and after the bolts are loose. erefore, the bolt looseness can be detected based on the change of the dynamic characteristics caused by the looseness damage.
Considering that the autocorrelation function of each measuring point is related to the structural dynamic characteristic parameters in Section 2.1, when the time lag T is 0, equation (19) can be expressed as where G ir is the coefficient depending on r th modal parameters and the response point i, and φ ir is the r th mode shape value at point i. From equation (20), it is can be seen that the autocorrelation function value R ii (0) at time lag T � 0 is composed of parameters such as mode shape and modal frequency. erefore, the autocorrelation function value at time lag T � 0 of the vibration response can be used for damage detection.
To reduce the influence of noise on the data results, it is usually necessary to obtain multiple sets of data under the same environment or segment the data obtained at one time.
en, the results of multiple sets of data are averaged for reducing the influence of noise. For this reason, this paper proposes using the average autocorrelation function R ii (0) at time lag T � 0 to establish the looseness damage indicator, and its calculation equation is where K is the number of signal sets, R x k x k (0) is the autocorrelation function value at time lag T � 0 of the k th collected data, and R ii (0) is the average autocorrelation function value at time lag T � 0. erefore, the relative change degree of the average autocorrelation function value at time lag T � 0 is calculated, and the bolt looseness damage index can be established as where δ is taken as the bolt looseness damage index, R ok (0) is the average autocorrelation function value in tightened 4 Shock and Vibration state, and R damage (0) is the average autocorrelation function value in loose state. It can be seen from equation (22) that the larger the damage index δ, the greater the relative difference of the average autocorrelation function value before and after the bolts are loose, and the greater the possibility of bolt looseness. Apparently, the smaller the bolt looseness damage index δ, the smaller the relative difference of the average autocorrelation function value before and after the bolts are loose, and the less the possibility of bolt looseness. erefore, the bolt corresponding to the maximum looseness damage index can be found out as the loose bolt.

Proposed Method and
Procedures. Based on the above analysis, a new method for bolt looseness detection based on average autocorrelation function is proposed. e flow chart of the proposed method, as shown in Figure 1, can be summarized as follows.
(1) For the undamaged structure, obtain the vibration signal x (m) of the location near the bolt Bm under ambient excitation (we denoted these bolts as B1, B2, . . ., B (L −1), and BL, m � 1, 2, . . . , L), and continuously collect K sets of data under the same environment. Take the k th set of data x k as an example, x k is [x k1 , x k2 , . . . , x (kN−1) , x kN ], and use equation (23) to calculate the autocorrelation function value at time lag T � 0.
(2) Average the autocorrelation function values of the K sets of data, and use equation (24) to calculate the average correlation function value at time lag T � 0.
(3) For the damaged structure, obtain the vibration signal y (m) of the location near the bolt Bm under ambient excitation, and continuously collect K sets of data under the same environment. Take the k th set of data y k as an example, y k is [y k1 , y k2 , . . . , y (kN−1) , y kN ], and use equation (25) to calculate the autocorrelation function value at time lag T � 0.
(4) Average the autocorrelation function values of the K sets of data, and use equation (26) to calculate the average correlation function value at time lag T � 0.
(5) Calculate the relative change degree δ of the average autocorrelation function value at time lag T � 0 before and after the bolts are loose.
where R (m) ok (0) is the average autocorrelation function value at time lag T � 0 before the bolt is loose,

Detecting System of Bolt Looseness Based on Computer
Vision. In order to further verify the feasibility of the proposed method in this paper, the sketch of the detecting system of bolt looseness is shown in Figure 2. e detecting system mainly consists of three parts: bolted connection plate, the shaker with its control system, and vision-measurement system. Further, the shaker with its control system consists of a shaker, a power amplifier, and a computer equipped with the LMS Test.Lab software. e visionmeasurement system consists of a camera, an optical lens, and a computer equipped with the image acquisition software.
In addition, Figure 3 shows the principle of visionmeasurement system. a is the image distance. b is the object distance, and O is the optical center of the lens. A is the spatial position of the object when the object is stationary, and B is the projection position of the object on the imaging plane when the object is stationary.
As shown in Figure 3, we assume that the direction of movement of the object is parallel to the imaging plane. When the time kΔt has passed, the object moves from A to A′, and its displacement in space is denoted as x. Further, the projection position of the object will move from B to B ′ , and its pixel displacement can be denoted as y. According to the similar triangle relationship in Figure 3 will find the following: If the time kΔt is, respectively, 1Δt, 2Δt, . . . , (N − 1)Δt, NΔt, the vibration signal of the object relative to the stationary position can be denoted as x (t).

Shock and Vibration 5
Correspondingly, the pixel displacement signal of the object relative to the stationary position can be expressed as y (t). en, combining equation (28) will result in the following: If the relative position of the camera and the object is fixed, both a and b will be constant. According to equation (29), we can find that the actual vibration displacement x (t) is linearly proportional to the pixel displacement y (t), and their scale factor is b/a. Based on the above analysis, when the bolt looseness damage index is calculated by equation (29), it can be seen that the scale factor b/a will be removed. Apparently, the damage index value calculated by using the pixel displacement is equal to the damage index value calculated by the actual displacement. erefore, we can directly use the pixel displacement signal to calculate the bolt looseness damage index. Figure 4 shows the experimental setup for bolt looseness detection, and the shaker and bolted connection plate are shown in detail in Figure 5. e bolted connection plate consists of an aluminum plate (size: 777 mm × 50 mm × 2 mm) and another aluminum plate (size: 510 mm × 50 mm × 2 mm) integrated by ten M6 bolts. In addition, the distance between the bolts is 50 mm, and we denoted these M6 bolts as B1, B2, . . ., B9, and B10. e experimental setup also consists of a shaker, a computer equipped with the LMS Test.Lab software, a power amplifier, and a vision-measurement system, as depicted in Figure 4. e vision-measurement system consists of a camera (IMI Tech/IMB-3213UP), an optical lens (TOKINA/TC1214-  Undamaged structure

Experimental Setup.
Calculate autocorrelation function value at T = 0 Damaged structure Calculate the average autocorrelation function value at T = 0 Obtain K sets of data 3MPG), and a computer equipped with the image acquisition software. e camera uses a CMOS sensor as the image receiver to collect 8-bit grayscale images at a speed of 300 fps. ese images are transferred to the computer through the USB3.0 interface to obtain the vibration video signal. Table 1 shows the specific parameters of visionmeasurement system. e shaker is controlled by a computer, which is equipped with the LMS Test. Lab software. e driving signal can be manually set as an ambient excitation to drive the vibration of the bolted connection plate. e bottom right of Figure 4 is a torque wrench, which can be used to adjust the pretension force of the bolt. ere are a total of ten M6 bolts on the bolted connection. Marked points are pasted near each bolt on the plate to establish eye-catching features, and we denoted these marked points as P1, P2, . . . , P9, and P10. As shown in Figure 6, we can see the position of the marked points from a partial image.

Experimental Method.
To avoid the bolt pretension force from exceeding the aluminum material yield strength, when we use the torque wrench to set the pretension force to 3 N·m, the bolt is considered to be perfect integration. When the pretension force is set to 0 N·m, the bolt is considered to be complete looseness. When the pretension force is set to 1.5 N·m, the bolt is considered to be half looseness.
e Gaussian white noise excitation with the frequency from 15 to 36 Hz and duration 2 minutes is taken as an input in the experiment. In addition, the frame rate of camera is set to 300 fps. When the bolted connection plate is vibrating, a sequence of images can be recorded by the camera, and there is no need for camera calibration in this process. After obtaining the sequence of images, python and open source computer vision (OpenCV) libraries are used to process these images for obtaining the vibration pixel displacement signal of each marked point. en, the relative difference of the average autocorrelation function value will be calculated. e following section gives detailed inspection results of different locations of bolt looseness.

Experiments for Bolt Looseness Location
e pixel displacement signal of each marked point in tightened state can be obtained, and 3 sets of data can be obtained continuously under the same environment. When the pretension force of bolt B2 is set to 0N·m (complete looseness) by a torque wrench, these pretension forces of other bolts are set to 3N·m. e pixel displacement signal of each point in the loose state can be obtained by the vision-measurement system, and 3 sets of data can be obtained continuously. Taking the first set of data as an example, vibration pixel displacements at P1 are shown in Figure 7, where Figure 7(a) is the sample pixel displacement for the undamaged structure, and Figure 7(b) is pixel displacement for the structure of bolt B2 looseness.
By the above vibration pixel displacements at P1, the autocorrelation functions can be calculated. Figure 8 shows the autocorrelation function at P1 before and after damage.
Because the autocorrelation function of the structural vibration displacement response is related to the structural modal parameters, the autocorrelation function value at time lag T � 0 of each marked point is used to calculate the detection index. Finally, the histogram of bolt looseness damage index is shown in Figure 9.
From Figure 9, one observes that the looseness damage index of bolt B2 is the largest, which indicates that the relative change of the average autocorrelation function value at time lag T � 0 of bolt B2 is relatively large, while the relative changes of the average autocorrelation function value at time lag T � 0 of other bolts are relatively small. Consequently, it can be considered that bolt B2 has looseness. In fact, the actual looseness is achieved by adjusting the pretension force of bolt B2 to 0 N·m. erefore, the detection result of the proposed method is consistent with the actual bolt looseness location.

Bolt B5 Looseness.
To further verify the accuracy of the proposed method in this paper, correspondingly, the Shock and Vibration pretension force of bolt B5 is set to 0 N·m (complete looseness) by a torque wrench, and pretension force of other bolts are set to 3 N·m. e pixel displacement signal of each point can be obtained by the vision-measurement system. en, according to the flow chart of the proposed method described above, the histogram of the bolt looseness damage index can be shown in Figure 10.
As shown in Figure 10, the looseness damage index of bolt B5 is the largest, which indicates that the relative change of the average autocorrelation function value at time lag T � 0 of bolt B5 is relatively large. Consequently, it can be considered that bolt B5 has looseness. In fact, the actual looseness is achieved by adjusting the pretension force of bolt B5 to 0 N·m. erefore, the detection result of the proposed method is consistent with the actual bolt looseness location.

Bolts B2 and B5
Looseness. When the pretension force of bolts B2 and B5 is set to 0 N·m (complete looseness) by the torque wrench, these pretension forces of other bolts are set to 3 N·m. e pixel displacement signal of each point can be  Figure 6: e position of the marked points. obtained by the vision-measurement system. According to the method flow described above, the histogram of the bolt looseness damage index can be shown in Figure 11. From Figure 11, it can be seen that the looseness damage index of bolt B5 is the largest, followed by bolt B2, which indicates that the relative change of the average autocorrelation function value at time lag T � 0 of bolts B2 and B5 is relatively large. Consequently, it can be considered that bolts B2 and B5 have looseness. In fact, the actual looseness is achieved by adjusting the pretension force of bolts B2 and B5 to 0 N·m. erefore, the detection result of the proposed method is consistent with the actual bolt looseness location.

Experiments for the Influence of Bolt Looseness Degree.
To further explore the applicability of the proposed method, a verification experiment is carried out to study the influence of different looseness degrees. e pretension force of bolt B5 is set to 0N·m, 1.5N·m, and 1.8N·m to simulate, respectively, 100% looseness, 50% looseness, and 40% looseness. e procedure of experiment is carried out following the same procedure described in Section 3.2. Figure 12 shows the relationship of damage index for different looseness degrees.
As shown in Figure 12, the damage index value changes with the change of the looseness degree. It also can be seen that the damage index value of bolt B5 in the state of complete looseness is larger than the damage index value in the state of 50% looseness. It shows that the smaller the pretension force, the greater the looseness degree. And the greater the relative difference of its autocorrelation function at time lag T � 0, the greater the damage index. In addition, from Figure 12, it also can be seen that the damage index is basically proportional to the level of the looseness.    Shock and Vibration

Conclusions
In this research, a new method for bolt looseness detection based on average autocorrelation function was proposed. It did not usually directly use the original response signal to extract the looseness damage features but used the average autocorrelation function value at time lag T � 0 of the original signal to establish the looseness damage index. By comparing the relative changes of each bolt looseness damage index, the bolt corresponding to the maximum looseness damage index could be denoted as the loose bolt. An experimental system for bolt looseness detection based on computer vision was designed. Validation experiments were carried out to demonstrate the feasibility of the proposed method, in which different bolt looseness locations and degrees were evaluated. Last but not least, this proposed method could be extended to other complex equipment structures with bolted connections.
Data Availability e data and the MATLAB programs 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.