Monitoring of the Looseness in Cargo Bolts under Random Excitation Based on Vibration Transmissibility

Bolted joints are widely used in industrial applications and joint failure can cause a disastrous accident if loosening happens. Bolt loosening detection can be made by regular manual inspection or by using sensors based on different physical principles, such as acoustoelastic effect, piezoelectric active sensing, and electromechanical impedance. Compared with the above methods, vibration based bolt looseness monitoring using accelerometers is appealing for its economy and convenience for measurement. In this paper, cargo bolts looseness monitoring under random excitation is studied based on vibration transmissibility, which overcomes the drawback of commonly used vibration methods in finding local bolt looseness. Vibration transmissibility analysis only uses two vibration transducers to monitor bolt group looseness, where the vibration signal below the cargo bolts is used as the “input” and the other one above the cargo bolts is used as the “output.”)ere are 12 bolts in the cargo bolts studied in this paper, providing an essential clamping force to fix cargo during transportation. Six kinds of bolt group looseness with an increasing degree are simulated in the experiment. )e experimental analysis shows that variation of the spectral moment can be used to monitor the global variation of the torque wrench exerted on the cargo bolts.)e early stage of the bolt group looseness is that some one or two bolts begin to loose; however, the spectrum moment factor is insensitive to the local bolt looseness in the bolt group. To address this issue, the eigensystem realization algorithm (ERA) based on random input and output is utilized to find the subtle eigenvalue variation of the system matrix, which is neglected by the frequency transmissibility function. )e experimental results show the effectiveness of the proposed method for detecting local bolt looseness.


Introduction
Bolted joints are widely used in mechanical and civil structures, such as machine tools, industrial robots, transport machines, bridges, and steel towers. A significant advantage of a bolted joint over other joint types, such as welded and riveted joints, is that they are capable of being disassembled. is feature, however, can cause joint failure leading to a disastrous accident if vibrational loosening happens, which is induced by flight/road loads or other environment uncertainties [1]. Hence, loosening detection of bolted joints is critical for structure safety, and it is highly desirable to achieve real-time online monitoring for bolts loosening detection under environmental and service conditions. e quality of a bolted connection is determined by bolt preload, which is usually controlled by the torque wrench method in engineering applications. However, it is difficult to achieve high precision of bolt preload because of the variation of friction between the bolt threads and the nut. Besides, the joints themselves may be located in inaccessible areas, which make regular inspection practically impossible.
As a nondestructive testing technology, acoustoelastic effect based methods can be applied to determine the bolt preload stress by measuring the time difference of reflecting waves [2]. However, the time difference is only in the range of tens of nanoseconds, and the data acquisition process requires a high sampling rate for the data acquisition system. Besides, most references based on the acoustoelastic effect use oscilloscopes to capture the signals, which are not suitable for the in situ bolt health monitoring [3]. e piezoelectric active sensing methods have been studied for bolt loosening detection. In the active sensing approach, one lead-zirconate-titanate (PZT) wafer is used as an actuator to generate an ultrasonic wave, and the other one is used as a sensor to record the wave that propagates through the bolted connection. e variation of the received wave energy [4][5][6], signal delay [7], spectral moment [8], and wavelet energy [9] have been experimentally investigated as an index for the bolting state. In the active sensing method, although the frequency for data acquisition can be lower than those used in the acoustoelastic effect based methods, it is still in the ultrasonic frequency band, which is quite a burden for data acquisition equipment in engineering applications.
Besides the abovementioned active sensing methods, the electromechanical impedance (EMI) techniques have also been studied for bolt loosening detection [10][11][12][13][14]. EMI measures the electrical impedance of a PZT actuator-sensor patch bonded to a bolted structure. Due to electromechanical impedance coupling, the change in the structure properties such as stiffness and damping due to bolt loosening can be detected through the variation of the measured electrical impedance. e impedance method is attractive for local dynamics of the bolt joint with a relatively lower sensing frequency than the piezoelectric active sensing methods [3]. However, for effective monitoring of bolts group joints, a lot of PZT actuator-sensor patches are needed for bolt loosening detection. Besides, the impedance signal can also be altered by other ambient variations such as temperature changes, making it susceptible to false alarms [10].
Vibration based bolt looseness monitoring using accelerometers is appealing for its economy and convenience for measurement. At present, much research about tool vibration based monitoring and fault diagnosis has been studied, which can be basically divided into two categories: (1) signal processing diagnostic method (2) and intelligent diagnostic method based on machine learning. Jiang and Yin generalized the recent advances in the multivariate statistical analysis based approaches for fault diagnosis and prognosis and developed a MATLAB toolbox for system fault detection (DB-KIT) [15]. Zhao et al. introduced an enhanced empirical wavelet transform for vibration based monitoring [16]; then they studied innovative machine learning methods for vibration based fault diagnosis [17] and prognosis [18]. e vibration based method is also effective for monitoring the connections with many bolted joints. e basic idea behind this technology is that any local damage somehow changes the dynamics of the structure and that the modal parameters identified from the structural response can reflect this change [19]. However, the measured vibration responses always reflect the global behavior of the whole structure. If there is another damage in the system that is under monitoring, or if the temperature and boundary conditions change, the bolt looseness monitoring based on vibration response will be inaccurate. To reduce the environmental effects, Cha et al. carried out a series of research on vision-based methods for bolt looseness detection, such as the Hough transform with support vector machines [20], the Viola-Jones algorithm with support vector machines [21], or convolutional neural networks [22,23], and these studies can effectively reduce the impact of environmental changes on detection accuracy.
Despite the limitations of the abovementioned vibration based techniques for bolt looseness monitoring, if there is a controlled input near the bolted structure, the vibration response can effectively reflect the bolt looseness. Razi et al. utilized an electric impact hammer to excite the bolted structure and analyzed the measured vibration using the empirical mode decomposition to establish an energy-based damage index [24]. With the impact hammer test, He and Zhu experimentally verified their study, which used the changes in natural frequencies to detect loosening of bolted connections in a full-size steel pipeline with bolted flanges [25]. Oregui et al. assessed the health condition of insulated rail joints in railway tracks by examining their dynamic response to impact excitation, where damage index is established by different frequency bands of vibration signals [26]. Wei et al. also utilized the exciting hammer to provide a stable excitation pulse signal and proposed a wavelet based damage index to monitor rail fastening status via a remote measurement system [27]. Kong et al. measured the responses of bolted joints, which were excited by impact, and employed the spectral features and the decision tree method to identify the bolt looseness level [28]. Using random excitation as a controlled input can be found in [29], where the authors studied statistical properties in time domain and signal power in different frequency bands for the measured vibration data, by which the bolt loosening can be detected.
In many practices for online monitoring of structural health, the controlled input is difficult to be achieved. us, transmissibility-based structural health monitoring has attracted much attention by utilizing the output-to-output relationship [30,31]. However, bolt looseness detection based on vibration transmissibility is rare in the literature, which is due to the fact that the responses on both sides of bolted connections may be uncorrelated at all. Fortunately, the cargo bolts in this study are another case: the motion of the whole structure is mainly in the vertical direction, and the responses below and above the bolted connection are obviously correlated.
In this paper, the vibration signal below the cargo bolts is used as the "input," and the vibration signal above the cargo bolts is used as the "output"; the concept of vibration transmissibility is used for bolt looseness monitoring. e frequency transmissibility function is calculated, and the derived spectrum moment factor is found to be related to the variation of torque level of the whole bolt group. However, the spectrum moment factor is insensitive to the local bolt looseness in the bolt group. To address this issue, the algorithm for the identification of a balanced state-space representation is utilized to find the subtle eigenvalue variation of the system matrix, which is neglected by the frequency transmissibility function. e first step of the proposed algorithm is the use of Markovsky's algorithm [32] to estimate a finite part of the impulse response from input-output data under random environmental excitation during cargo transportation. Once the impulse response is computed, the second step is to obtain a balanced state-space model via the eigensystem realization algorithm (ERA) [33], by which the eigenvalues of the identified system matrix can be calculated as the bolt looseness indicator. e organization of this paper is as follows: the second section describes the studied cargo bolts and bolt looseness simulation; the third section shows the effectiveness of the frequency transmissibility function on bolt group looseness while also showing its insensitivity to local bolt looseness; and the fourth and fifth sections give the solution based on identification in the state-space domain.

The Studied Cargo Bolts and Looseness
Simulation in Experiment e cargo bolts are used in a clamping support structure, which is designed to fix cargo during transportation, as shown in Figure 1. e enlarged view and cross section view of the structure are shown in Figures 2 and 3. e clamping support structure has 8 components as shown in Figure 3. e vertical braces (components 2 and 3) are fixed on the base (component 1). e vertical braces are also connected with two half-cylindrical clamps (components 4 and 5) by the cargo bolts (component 6 with 12 bolts, each side 6). When the cargo bolts are tightened, the simulated cargo (components 7 and 8, two half steel tubes) will be fixed by clamping force. e side view of the cargo bolts is also shown in Figure 2. Two vibration accelerometers below and above the cargo bolts are fixed on the structure. e studied cargo bolts in Figure 2 have twelve bolts, evenly spaced along the plane connecting boundary. Under stochastic base vibration excitation, the early stage of the bolt group looseness is that some one or two bolts begin to loose, which increases dynamic loads on adjacent bolts, leading to the increased likelihood of bolt looseness. When more bolts loose, the previous loose bolts may become tight again due to the elasticity of the connecting boundary. However, the pretightening force of all the bolts in the bolt group will decrease. us, six damage cases of the bolt group looseness are simulated in the experiment: ① 2 bolts at one corner are loosened while the other 10 bolts are tightened; ② 4 bolts at the four corners are loosened while the other 8 bolts are tightened; ③ 4 bolts at corners are tightened, while the left 8 bolts are loosened; ④ the torque wrench levels of all the 12 bolts are decreased with a slight degree; ⑤ the torque wrench levels of all the 12 bolts are decreased with a medium degree; ⑥ the torque wrench levels of all the 12 bolts are decreased with a serious degree. e illustration of the first two damage cases is shown in Figure 4.

The Proposed Algorithm for Cargo Bolts Looseness Monitoring
As discussed in Section 2, the motion of the cargo bolts is mainly in the vertical direction; thus, the concept of vibration transmissibility is suitable for bolt looseness monitoring, where the vibration signal below the cargo bolts can be used as the "input" and the vibration signal above the cargo bolts are used as the "output." In the following proposed algorithm, the frequency transmissibility function is calculated, and the derived spectrum moment factor is utilized to monitor the torque level variation of the whole bolt group. A brief theory description is shown in Section 3.1 and corresponding experimental results are shown in Section 4.
To further reflect the local bolt looseness in the bolt group, the algorithm for the identification of a balanced state-space representation is utilized to find the subtle eigenvalue variation of the system matrix. In the proposed state-space method, the first step is to estimate a finite part of  Shock and Vibration 3 the impulse response from input-output data under random environmental excitation during cargo transportation; the second step is to obtain a balanced state-space model via the eigensystem realization algorithm (ERA), by which the eigenvalues of the identified system matrix can be calculated as the local bolt looseness indicator. e theoretical background is introduced in Section 3.2 and verification is shown in section 4.2. Figure 2, two vibration accelerometers below and above the cargo bolts are fixed on the structure. In this paper, the vibration signal below the cargo bolts is used as the "input" and the vibration signal above the cargo bolts is used as the "output"; vibration transmissibility analysis is used for bolt looseness monitoring.

Spectral Moment Factor. As shown in
At time instant k, let the sampled vibration signal below the cargo bolts be u(k), and let the sampled signal above the cargo bolts be y(k). e frequency transmissibility function (FTF) S(ω) can be calculated using the cross power spectrum G yu (ω) between u(k) and y(k) and using auto power spectrum G uu (ω) of u(k) [31]: To quantify the cargo bolts looseness level, the spectral moment factor is defined as shown in equation (2), where only one input and one output are considered for simplicity. Spectral moments can describe the spectral distribution of a statistical signal, such as structural responses [34], EEG signals [35], and wave energy [36]. Spectral moments have been adopted for assessing the loosing state of a bolted structure by the piezoelectric active sensing method [8], where the requirements for additional energy consumption and high sampling frequency are needed. In this paper, spectral moments of vibration transmissibility function are calculated and used as an indicator of the torque level variation of the whole bolt group.
Spectral moment with 0 and 1 order: Spectral moment factor:

Estimation of the Impulse Response from Input-Output Data under Random Excitation.
A state-space model for a dynamic system can be written as follows: e matrix A ∈ R n×n is the state transition matrix that completely characterizes the dynamics of the system by its eigenvalues.
e vectors u(k) ∈ R m×1 , y(k) ∈ R p×1 , and x(k) ∈ R n×1 are the m inputs, p outputs, and n states, respectively, of the system at time instant k. To consider process noise and measurement noise of the model, noise vectors w(k) and v(k) are added into the model.
Given the input vector u(k) and the output vector y(k), the state-space identification algorithms can be used to identify the system matrices A, B, and C, where the orthogonal or oblique projection of input and output data matrices are calculated through the QR decomposition, followed by Singular Value Decomposition (SVD) to determine the order, the observability matrix, and the state sequence. en, the extraction of the state-space model is achieved through the solution of a least-squares problem.
ere are various kinds of state-space identification algorithms in the literature, such as N4SID, MOESP, and CVA, and Van Overschee incorporated them with different userdefined weighting matrices [37]. If only the estimation of the matrices A and C is concerned, from which modal parameters can then be extracted, the algorithms are generally equivalent [38].
Markovsky proposed a computationally efficient algorithm for the computation of a finite part of the impulse responses directly from the input-output data [32]. Once the impulse response is computed, a balanced state model can be derived without returning to the original observed data.
Suppose we have N samples of input and output data: the problem is how to compute the first L samples of impulse response matrices h(k) ∈ R p×m , k ∈ 0, . . . , L − 1 { }. e common sense is that N > L and L > n max . n max is an upper bound of the system order n.
Similar to the definition of block Hankel matrices in [37], the following notation is used with f � f(0), . . . , f(N − 1) { }: H Equation (5) defines the block Hankel matrix of f having Δ block rows. en, through block Hankel matrices of input u and output y, matrices of U p , U f , Y p , and Y f can be defined as follows: and row dim denotes the number of block rows of a matrix or vector. l max is an upper bound of the system lag l.
e impulse response sequence with L samples can be defined as a matrix with pL rows and m columns: where T denote the transpose of a matrix in equation (7). Clearly, H pL×m is the result of the system input as follows: where I denotes the identity matrix with m rows and columns and 0 denotes the zero matrix with L(m − 1) rows and m columns in equation (8).
Under zero initial conditions, there exists a matrix G, such that Since G can be calculated through the first three submatrices in equation (9), the impulse sequences can be calculated as H pL×m � Y f G. e calculation can be implemented via QR decomposition, and the details can be referred to in [32].

ERA for System Matrix Estimation.
eoretically, the impulse sequences, or the Markov parameters, can be denoted using matrices A, B, and C defined in equation (4): ERA begins by forming the pα × mβ Hankel matrix H(k) composed of the impulse responses h(k): (11) It can be derived that the H(k) is connected with observability matrix P α , the system matrix A , and controllability matrix Q β : e ERA process starts with the factorization of the Hankel matrix H(0) by the use of SVD: where the columns of matrices U and V are orthonormal and Σ is a rectangular matrix: in which 0 denotes zero matrix with an appropriate dimension, and the matrix Σ r is where σ i is monotonically nonincreasing: If we define U r and V r as the first r columns of the orthogonal matrices U and V, then we have Because of measurement noise, nonlinearity, and computer round-off error, the Hankel matrix H(0) will usually be of full rank, which does not, in general, equal the true order of the system under test. By using SVD, noise and other interferences are removed through the truncation of nonzero small singular values.
A balanced realization can be obtained by choosing the truncated observability matrix and controllability matrix as follows: To compute the system matrix A, we can use the Hankel matrix H(1):

Shock and Vibration 5
en, the state transition matrix A can be calculated as

Physical Mode Determination through Stabilization
Diagram and Clustering. e estimated system matrix A can be decomposed by eigendecomposition as A � ΨΛΨ − 1 , where Ψ is the eigenvector matrix and Λ � diag(λ di ) is the diagonal eigenvalue matrix. e matrix Λ contains the (n/2) discrete-time eigenvalues λ di and λ * di in complex conjugated pairs, which are directly related to the frequency and damping properties of the system. e discrete-time eigenvalues are first converted to continuous time eigenvalues λ ci as λ ci � (ln(λ di )/dt) , where dt is the time step of the digital data acquisition system. e undamped natural frequencies f ni and damping ratios ζ i can then be easily calculated from the conjugate pair of complex-valued eigenvalues: λ ci , λ * ci � − 2πζ i f ni ± j(2πf ni ) ����� 1 − ζ 2 i . As the true system order is often unknown, a common practice is to calculate the modal parameters for increasing model orders n. Besides, the weakly excited modes, which appear only at high model orders, can be detected. But if the model order n is higher than the true system order, the irrelevant noise will be modeled. ese irrelevant noise modes can be visually removed in the stabilization diagram, in which the eigenfrequencies are plotted against the increasing model orders. e task can be very tedious in cases when there are closely coupled modes, heavily damped modes, and so forth.
is creates a need for automating this selection process, thus avoiding user judgment by formulating a more reliable and intelligent process. A clustering algorithm is adopted to select the stable modes from the large set of parameters identified at various model orders. In this approach, one single set of modal estimates (a frequency and the associated damping ratio) is chosen among all the available parameters identified, and all sets that are closer to this point than a certain "radius" are clustered based on a distance measure between frequencies and damping estimates. e distance between mode i and j is a Euclid distance that incorporates the eigenfrequency and damping ratio difference: When certain clusters are selected, the variance of the clusters should be checked. When the variances of eigenfrequencies or damping ratios in each cluster are larger than a user-defined threshold, a conclusion of unsuccessful clustering will be made, and a new modal identification procedure should be restarted.

Flowchart of the Proposed Algorithm for Cargo Bolts
Looseness Monitoring. Based on the theoretical description in Sections 3.1 and 3.2, the flowchart of the cargo bolts looseness monitoring is shown in Figure 5. e proposed algorithm has four major steps: (1) Cargo bolts are tightened for transportation. Because the cargo weight and bolt connection stiffness influence the calculated frequency transmissibility function, calibration is needed to calculate a baseline and threshold for a specific cargo with the tightened bolt connections. e algorithm for calibration is the same in the following monitoring procedure. e recalibration process for each transportation also helps to reduce the environmental impact, because the time for each transportation is relatively short.
(2) During cargo transportation, vibration signals below and above the bolted connection are measured as the "input" and "output" data for transmissibility analysis. e measurement is repeated for a constant duration. In this application, the duration is set as four seconds. e most time-consuming part is the estimation of the impulse response from inputoutput data, which is around 0.7 seconds in the experimental data processing. Other parts of the monitoring algorithm are less than 0.3 seconds. If the consuming times for consecutive data acquisition, display, and storage are included, the four-second duration is reasonable for online monitoring. (3) During each four-second cycle, the spectral moment factor is calculated to monitor the torque level variation of the whole bolt group. A significant decrease of the spectral moment factor reflects a severe bolt loosening. When the factor is below the preset threshold, an emergency stop should be made for cargo bolts retightening. (4) During each four-second cycle, the subtle eigenvalue variation of the system is calculated as the local bolt looseness indicator. When an abnormal frequency appears, while the spectral moment factor shows no noticeable change, a warning will arise in the monitoring system. In this case, the transportation goes on as usual until the next scheduled stop, where bolts connection state will be checked and retightened.

Spectral Moment Factor as Looseness Indicator of the Cargo Bolts.
With the whole structure under vertical random excitation in the experiment, the calculated FTFs for cargo bolts under two states are shown in Figures 6 and 7. It can be seen in the figures that the spectral components are basically the same, while the peak positions decreased when bolt looseness happens. is is because of the decreased connection stiffness due to bolt looseness.
Variation of the defined spectral moment factor with the looseness level of cargo bolts is shown in Figure 8. e looseness level 0 means the cargo bolts in normal condition. Looseness levels 1-6 correspond with the six damage cases defined in Section 2. Because the spectral components are mainly located in 0-1000 Hz, as shown in Figures 6 and 7, H is set as 1000 Hz when calculating the spectral moments in equation (2).
It can be seen from Figure 8 that the decreased value of spectral moment factor corresponds with the increase of looseness level for the cargo bolts. e vibration transmissibility analysis can be an effective tool for monitoring the cargo bolts looseness. However, the spectral moment factor is insensitive to local bolt looseness of damage case ①. To further increase the sensitivity of detecting local bolt looseness, vibration transmissibility analysis in the statespace domain is introduced and verified in the following sections. e experiment has been carried out three times, and the results have good repeatability.

Local Bolt Looseness Detection via Subtle Eigenvalue
Variation of the System Matrix

Simulation Demonstration.
A simulation example is given at first to give an explanation of the subsequent experimental results. Figure 9 shows a two-DOF system with base random excitation. e acceleration measured on the base plate is used as the system input, while the acceleration measured on the top mass element is used as the system output. e two-DOF system used in this simulation is tuned to have two close spaced natural frequencies with a relatively high damping ratio. e two natural frequencies are set as 1 Hz and 1.025 Hz, and the damping ratio is set as 0.016. Uncorrelated white noises w(k) and v(k) shown in equation (4)      include uncertainty effects. Equation (4) describes a generalized state-space model with multi-input and multioutput. In this study, the input and output are all simplified to one channel. Figure 10 shows the estimated impulse response of the two-DOF system from the two measured accelerations under random excitation. e sampling frequency is 25 Hz, and signal-noise-ratio (SNR) is set as 20 dB. e calculated FTF using equation (1) is shown in Figure 11. Only one peak frequency with 1 Hz can be found, while the other natural frequency of 1.025 Hz is lost. e lack of frequency resolution is due to the effect of damping, which cannot be solved by simply increasing the acquired data length. e result using state-space domain identification is shown in Figure 12. e identified eigenfrequencies are confirmed to be a true physical mode because they are stable when the system order is increasing. e identified eigenfrequencies are 1 and 1.024 Hz, which are very close to their actual value.
When selecting the physical mode from the candidate eigenfrequencies, which are calculated from the system matrix A with a relatively high order, cluster analysis in the frequency-damping map is used. e distribution of identified frequency-damping is shown in Figure 13, where the noise effect is compared. When there is no noise, the identified mode is quite repeatable. When SNR is 20 dB, the identified results are very close to their actual value. As the SNR is decreasing, the identified results deviate from the actual value, but the error is still small. For example, the identified eigenfrequencies are 1.006 Hz and 1.035 Hz, respectively, and the error is less than 1%.
A comparative study is made for different impulse responses estimation. In this paper, the impulse responses h(k) used in ERA are estimated using the state-space model as shown in Section 3.2.1, while the commonly used method for h(k) estimation is the inverse Fourier transform (IFFT) from the frequency transmissibility function S(ω).
e results using h(k) are shown in Figures 14 and 15, which can be compared with Figures 12 and 13, where the length of impulse responses and parameters used in ERA is the same. Because there is variance in power spectra estimation for random signals using equation (1) [39], the estimated h(k) by IFFT of S(ω) has a larger error than h(k) using the state-space model. ere are three closely spaced frequencies for different model orders, as shown in Figure 14.
Cluster analysis is carried out in the identified frequencydamping map of Figure 15. ey are identified as one cluster, and 1.008 Hz with damping ratio 0.045 is identified using equation (21). e identified result cannot actually separate the two closely spaced modes of 1 Hz and 1.025 Hz.       Figure 16 and Figure 17. It can be seen from Figure 16 that the estimated base frequency using FTF is 202.5 Hz for damage case ①, while the base frequency of healthy cargo bolts is 203.5 Hz shown in Figure 6. e difference is quite a nuance and will be neglected in actual monitoring. When the bole looseness level is increased from case ① to case ②, the base frequency in the calculated FTF is decreased to 197 Hz. e results of Figures 16 and 17 can explain the result of the spectral moment factor in Figure 8, where the spectral moment factor for case ① does not changes while the spectral moment factor for case ② has a visible decrease.
Analysis using the algorithm in Section 3.2 is carried out in this section. Firstly, the estimated impulse response for the cargo bolts in damage case ① is shown in Figure 18. To avoid noise effect, impulse response with very small amplitudes is cut off. e sampling frequency is 3840 Hz. e second step is to identify the system matrix A and calculate its eigenvalues using the estimated impulse response. From Figure 19 and Figure 20, the identified eigenfrequencies for the cargo bolts in damage case ① can be confirmed to be 196 Hz, 203 Hz, and 230 Hz. e identified base frequency is 196 Hz, which cannot be detected using FTF in Figure 16.      Shock and Vibration e results of Figure 20 show that the identified base mode with 196 Hz is quite concentrated in the frequencydamping map, which shows relatively high energy compared with the scattered result of the third mode with 230 Hz. e emerge of 196 Hz is induced by asymmetric stiffness, which is caused by the local bolt looseness in the corner as shown in Figure 4. However, the base mode of 196 Hz is too close to the second mode of 203 Hz, which covers the base mode in the frequency domain using FTF in Figure 16. e results of Figure 19 and Figure 20 show that the analysis of vibration transmissibility in the state-space domain can discover the subtle eigenvalue variation of the system matrix.
When the bolt looseness level is increased from case ① to case ②, the identification results are shown in Figures 21  and 22. e base mode decreases from 203 Hz to 197 Hz due to the stiffness loss. As the damage case ② is four-bolt looseness at four corners, shown in Figure 4, the asymmetric effect does not exist. us, there are only two strong and distinct modes, and the frequency-damping is quite concentrated in Figure 22.
A comparative study is also made for case ① in this section for different impulse responses estimation, as conducted in the simulation Section 4.2.1. e results using h(k) are shown in Figures 23 and 24, which can be compared with Figures 19 and 20, where the calculation parameters are the same. From the comparison of Figures 19 and 23, it can be seen that the estimated h(k) by IFFT of S(ω) has a larger error than h(k) using the state-space model. ere are also three frequency lines in Figure 23, but the lines have more variance than Figure 19. Cluster analysis is carried out in the identified frequency-damping map of Figure 24, where it can be seen that the points of the first mode have a large dispersion area. e first two modes are identified as one cluster, and only two modes are identified with frequencies 199.8 Hz and 231.1 Hz.

Conclusion
is paper gives a case for vibration based bolt group looseness monitoring. e method uses the concept of vibration transmissibility. e motion of the whole structure studied in this paper is mainly in a vertical direction; thus, the use of the responses below and above the bolted connection as the "input" and "output" data for transmissibility analysis is feasible for bolt looseness detection. e frequency transmissibility function is calculated, and the derived spectrum moment factor is found to be related to the variation of torque level of the whole bolt group, where the decreased value of spectral moment factor corresponds with the increase of looseness level for the cargo bolts. e early stage of the bolt group looseness is that some one or two bolts begin to loose; however, the spectrum moment factor is insensitive to the local bolt looseness in the bolt group. To address this issue, the algorithm for the identification of a balanced state-space representation is further utilized to find the subtle eigenvalue variation of the system matrix, which is neglected by the frequency transmissibility function. e experimental results can discover the subtle eigenvalue variation of the system matrix and show the effectiveness of the proposed method for detecting local bolt looseness.

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

Conflicts of Interest
e authors declare that they have no conflicts of interest.