Analysis for the Vibration Mechanism of the Spillway Guide Wall Considering the Associated-Forced Coupled Vibration

Featured Application: The ﬂow-induced vibration of a hydraulic structure is a threat to the safety of hydraulic engineering projects. Therefore, a new vibration mechanism is presented, which can be applied to control hydraulic structural vibration in structural design and maintenance. Abstract: During the ﬂood discharge in large-scale hydraulic engineering projects, intense ﬂow-induced vibrations may occur in hydraulic gates, gate piers, spillway guide walls, etc. Furthermore, the vibration mechanism is complicated. For the spillway guide wall, existing studies on the vibration mechanism usually focus on the vibrations caused by ﬂow excitations, without considering the inﬂuence of dam vibration. According to prototype tests, the vibrations of the spillway guide wall and the dam show synchronization. Thus, this paper presents a new vibration mechanism of associated-forced coupled vibration (AFCV) for the spillway guide wall to investigate the dynamic responses and reveal coupled vibrational properties and vibrational correlations. Di ﬀ erent from conventional ﬂow-induced vibration theory, this paper considers the spillway guide wall as a lightweight accessory structure connected to a large-scale primary structure. A corresponding simpliﬁed theoretical model for the AFCV system is established, with theoretical derivations given. Then, several vibrational signals measured in di ﬀ erent structures in prototype tests are handled by the cross-wavelet transform (XWS) to reveal the vibrational correlation between the spillway guide wall and the dam. Afterwards, mutual analyses of numeral simulation, theoretical derivation, and prototype data are employed to clarify the vibration mechanism of a spillway guide wall. The proposed mechanism can give more reasonable and accurate results regarding the dynamic response and amplitude coe ﬃ cient of the guide wall. Moreover, by changing the parameters in the theoretical model through practical measures, the proposed vibration mechanism can provide beneﬁts to vibration control and structural design.


Introduction
In recent years, comprehensive hydraulic engineering projects with high dams and large reservoirs have been rapidly developing, especially in China, such as Three Gorges, Jinping, and Xiluodu. Simultaneously, flow-induced vibrations (FIV) of hydraulic structures, which frequently occur in these projects, have posed a significant threat to engineering safety [1]. To investigate the mechanism of the FIV, researchers have proposed several types of source and mechanism of the FIV. Blevins [2] transform (WT) technique and its modifications have become widely utilized in signal processing. Extended from WT, the cross-wavelet transform (XWT) is a new method to evaluate the relationship between two signals. This method is more common in meteorological [23], hydrological [24] and medical [25] fields, and is demonstrated to be an efficient method of analyzing the vibrational correlation between two structures [26].
In conclusion, the studies on the FIV of spillway guide walls are inadequate, as the conventional methods usually focus on the vibrations caused by flow excitations, but the influence of dam vibration is always ignored, leading to the inaccuracy and irrationality of vibration analyses. In this paper, according to prototype tests and structural characteristics of a typical guide wall, a new perspective of associated-forced coupled vibration (AFCV) is proposed. Considering the lightweight characteristic of the guide wall, it can be deemed to be an accessory structure installed in the largescale primary structure, constituting the associated vibration (AV) system. Then, considering the connection between the guide wall and the foundation, the AV system was modified to be an AFCV system, with the establishment of a two-degrees-of-freedom theoretical model. The proposed theoretical model is subsequently verified based on prototype tests. By employing the signal processing method XWT and numerical simulation, the vibration mechanism can be clarified in detail utilizing the theoretical model, and simultaneously demonstrates the reasonability and correctness of the theoretical model. In practice, this investigation aims at contributing to the design, monitoring, and maintenance of spillway guide walls, and is beneficial to the vibrational control of similar lightweight hydraulic structures.

Prototype Tests for the Spillway Guide Wall Vibration
Three Gorges Hydropower Station (TGHS) is currently the largest hydropower station in the world, during the flood discharge, its structural safety is of utmost concern. The left spillway guide wall of TGHS is located between the left powerhouse and the spillway dam, as shown in Figure 1. During the flood discharge, the action area of fluctuation is large, and the action point is high, making the FIV a severe threat to the guide wall. To investigate the vibrational characteristics and mechanism of this guide wall, prototype tests of structural displacement were implemented. In the tests, several displacement sensors were installed on 1# spillway dam block and eight guide wall blocks (1#-8#). There are four measuring points on the dam block: 1# and 2# measuring points were located at the upstream side, 3# and 4# measuring points were located at the downstream side. Horizontal displacement sensors were arranged for each of the measuring points, and two vertical displacement sensors were arranged for 1# and 3# measuring points. Similarly, on each guide wall block, the upstream measuring point 1# had one horizontal displacement sensor, and the downstream measuring point 2# had one horizontal displacement sensor and one vertical displacement sensor. All of the measuring points were tested simultaneously under the same In the tests, several displacement sensors were installed on 1# spillway dam block and eight guide wall blocks (1#-8#). There are four measuring points on the dam block: 1# and 2# measuring points were located at the upstream side, 3# and 4# measuring points were located at the downstream side. Horizontal displacement sensors were arranged for each of the measuring points, and two vertical displacement sensors were arranged for 1# and 3# measuring points. Similarly, on each guide wall block, the upstream measuring point 1# had one horizontal displacement sensor, and the downstream measuring point 2# had one horizontal displacement sensor and one vertical displacement sensor. All of the measuring points were tested simultaneously under the same working condition (listed in the Appendix A). The layout of the measuring points is illustrated in Figure 2. For the low dominant frequency of the vibration on the dam and the guide wall, the sampling frequency is 50 Hz.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 4 of 25 working condition (listed in the Appendix). The layout of the measuring points is illustrated in Figure  2. For the low dominant frequency of the vibration on the dam and the guide wall, the sampling frequency is 50 Hz. The signals of typical measuring points for vibration displacements on the 1# dam block and each of the guide wall blocks is processed by the Fourier Transform (FT). The FT spectrums of horizontal and vertical displacements of the 1# dam block and typical guide wall blocks are respectively shown and compared in Figure 3 and  From the FT spectrums, it can be observed that there are several similar dominant frequencies of horizontal displacements between the 1# dam block and the guide wall blocks in a broad range of 0-15 Hz. However, for vertical displacements, similar dominant frequencies are concentrated in a narrow range of 0-5 Hz, which is the range of the dominant frequency of flow excitation [5,27]. For further investigation of this phenomenon, the dominant frequencies of horizontal displacements on the 1# dam block and guide wall blocks are listed in Table 1.   As shown in Table 1, apart from the dominant frequencies of flow fluctuations (less than 5 Hz), some of the dominant frequencies of the horizontal displacements of guide wall blocks tend to match that of 1# dam block in an approximate range: the range of 5.58-5.81 tends to match 6.04; 7.14-7.22 tends to match 7.11; 16.26 tends to match 16.27. The vibrational components with frequencies higher than 5 Hz are considered to represent the structural resonances of certain modes under flow excitations because the fundamental frequencies of structures are usually higher than 5 Hz, which will be analyzed in detail in Section 4.2. It can be inferred that the horizontal vibrational characteristics of the dam block and the guide wall blocks are relevant to some extent, which can postulate that the vibration propagated from the dam may participate in the guide wall vibration for certain defined frequencies. To verify the postulation and judge the vibrational relationship between 1# dam block and guide wall blocks, the theoretical analyses of the guide wall must be proposed.

The Guide Wall Vibration under FSI
It is noted that the guide wall is a hydraulic structure that is partially submerged in water, as shown in Figure 5. For such structures, the vibration is nonlinear and under the effect of the fluidstructure interaction (FSI) [27,28], which means the structural loads and responses are mutually From the FT spectrums, it can be observed that there are several similar dominant frequencies of horizontal displacements between the 1# dam block and the guide wall blocks in a broad range of 0-15 Hz. However, for vertical displacements, similar dominant frequencies are concentrated in a narrow range of 0-5 Hz, which is the range of the dominant frequency of flow excitation [5,27]. For further investigation of this phenomenon, the dominant frequencies of horizontal displacements on the 1# dam block and guide wall blocks are listed in Table 1.   As shown in Table 1, apart from the dominant frequencies of flow fluctuations (less than 5 Hz), some of the dominant frequencies of the horizontal displacements of guide wall blocks tend to match that of 1# dam block in an approximate range: the range of 5.58-5.81 tends to match 6.04; 7.14-7.22 tends to match 7.11; 16.26 tends to match 16.27. The vibrational components with frequencies higher than 5 Hz are considered to represent the structural resonances of certain modes under flow excitations because the fundamental frequencies of structures are usually higher than 5 Hz, which will be analyzed in detail in Section 4.2. It can be inferred that the horizontal vibrational characteristics of the dam block and the guide wall blocks are relevant to some extent, which can postulate that the vibration propagated from the dam may participate in the guide wall vibration for certain defined frequencies. To verify the postulation and judge the vibrational relationship between 1# dam block and guide wall blocks, the theoretical analyses of the guide wall must be proposed.

The Guide Wall Vibration under FSI
It is noted that the guide wall is a hydraulic structure that is partially submerged in water, as shown in Figure 5. For such structures, the vibration is nonlinear and under the effect of the fluidstructure interaction (FSI) [27,28], which means the structural loads and responses are mutually From the FT spectrums, it can be observed that there are several similar dominant frequencies of horizontal displacements between the 1# dam block and the guide wall blocks in a broad range of 0-15 Hz. However, for vertical displacements, similar dominant frequencies are concentrated in a narrow range of 0-5 Hz, which is the range of the dominant frequency of flow excitation [5,27]. For further investigation of this phenomenon, the dominant frequencies of horizontal displacements on the 1# dam block and guide wall blocks are listed in Table 1. As shown in Table 1, apart from the dominant frequencies of flow fluctuations (less than 5 Hz), some of the dominant frequencies of the horizontal displacements of guide wall blocks tend to match that of 1# dam block in an approximate range: the range of 5.58-5.81 tends to match 6.04; 7.14-7.22 tends to match 7.11; 16.26 tends to match 16.27. The vibrational components with frequencies higher than 5 Hz are considered to represent the structural resonances of certain modes under flow excitations because the fundamental frequencies of structures are usually higher than 5 Hz, which will be analyzed in detail in Section 4.2. It can be inferred that the horizontal vibrational characteristics of the dam block and the guide wall blocks are relevant to some extent, which can postulate that the vibration propagated from the dam may participate in the guide wall vibration for certain defined frequencies.
To verify the postulation and judge the vibrational relationship between 1# dam block and guide wall blocks, the theoretical analyses of the guide wall must be proposed.

The Guide Wall Vibration under FSI
It is noted that the guide wall is a hydraulic structure that is partially submerged in water, as shown in Figure 5. For such structures, the vibration is nonlinear and under the effect of the fluid-structure interaction (FSI) [27,28], which means the structural loads and responses are mutually influential. To investigate the effect of FSI, the coupled equations of motion (EOM) of FSI are discussed below.  Suppose that the fluid around the guide wall is homogenous, incompressible, non-viscous and irrotational. Therefore, the hydrodynamic pressure and velocity potential of the fluid satisfy the following relationships: where p denotes the hydrodynamic pressure field; φ denotes the velocity potential; ∇ is the del operator. Employing the Galerkin method to discretize Equation (1), the discrete fluid motion equation can be deduced as follows [29,30]: where ρ denotes the fluid density; r denotes the displacement vector; 0 q denotes the input excitation vector; and where N and s N denote the shape functions of the fluid and fluid-structure interface, respectively; r S , f S , and l S denote the areas of infinite boundary, free surface, and the fluid-structure interface, respectively; Ω denotes the volume of the fluid; w C denotes the compression wave velocity; g denotes the gravity acceleration; Λ denotes the coordinate transformation matrix; the superscript T Suppose that the fluid around the guide wall is homogenous, incompressible, non-viscous and irrotational. Therefore, the hydrodynamic pressure and velocity potential of the fluid satisfy the following relationships: where p denotes the hydrodynamic pressure field; φ denotes the velocity potential; ∇ is the del operator. Employing the Galerkin method to discretize Equation (1), the discrete fluid motion equation can be deduced as follows [29,30]: where ρ denotes the fluid density; r denotes the displacement vector; q 0 denotes the input excitation vector; and where N and N s denote the shape functions of the fluid and fluid-structure interface, respectively; S r , S f , and S l denote the areas of infinite boundary, free surface, and the fluid-structure interface, respectively; Ω denotes the volume of the fluid; C w denotes the compression wave velocity; g denotes the gravity acceleration; Λ denotes the coordinate transformation matrix; the superscript T denotes the transpose of the corresponding matrix. It is noted that the matrix H represents the "stiffness" of the fluid, and the matrix B represents the coupled characteristics on the fluid-structure interface. The EOM of a structure in the fluid is given as the following: where x is the displacement vector; M s , C, and K denote the mass, stiffness and damping matrices of the structure, respectively; f p denotes the hydrodynamic nodal forces on the fluid-structure interface; f 0 denotes the external excitation apart from f p . Gathering the contributions of each fluid element, f p can be obtained as the following: and Equation (5) can be rewritten as follows: It is noted that the dynamic behavior of a fluid-structure system considering FSI can be obtained as the combination of Equations (3) and (7). The combined EOM is given as follows: As an arbitrary load can be Fourier transformed into the sum of harmonic excitations with different frequencies, the input excitation q 0 can be considered as a harmonic force, Matching the form of the input excitation, the solution of Equation (7) can be expressed as follows: As mentioned above, the fluid is considered to be incompressible ( C w → ∞ ), which means A → 0 ; the effect of the free surface is usually neglected, so E → 0 . Moreover, the external excitation f 0 is equal to zero in this case as the excitation is flow excitation only. Substituting into Equation (8), the equation can be obtained as follows: The displacement x 0 can be obtained as follows: Equation (12) can be rewritten as follows: .. where Equation (13) is a common form of EOM. It is noted that P is a coupled input excitation. In practical prototype tests, the value of P can be measured as the fluctuation pressure acting on the fluid-structure interface, which is an actual value that can take all the coupled characteristics into consideration. Therefore, the effect of FSI can be concluded as M a , which is defined as the added mass in existing studies. The value of added mass is difficult to calculate. In existing studies, several models have been established to calculate the approximate value of the added mass. For large hydraulic structures such as dams, the most classic method is the Westergaard method [31,32]. The added mass can be calculated as follows: where ρ denotes the water density; h denotes the depth of water above the base; z denotes the distance from the surface. Concluded from practical hydraulic projects [5], the value in Equation (15) is larger than reality. In the engineering field, the added mass is usually considered as follows: With an added mass, the effect of FSI can be taken into consideration, and a more simplified theoretical model of the guide wall can be established.

Theoretical Model
The hydraulic structures, including dams, spillway guide walls, hydraulic gates, are complex multi-degrees-of-freedom (MDOF) systems with complicated connections and damping conditions. For ease of calculation, the complex MDOF system should be simplified appropriately. This paper mainly concerns the vibrational response of the guide wall and the dam so that the system can be transformed into a two-degrees-of-freedom system. It is worth mentioning that the guide wall has a lighter weight compared to the weight of the large-scale concrete dam. To investigate the simplified form of the vibration system in this circumstance, the existing research of the hydraulic gate, another type of lightweight hydraulic structure, can provide a reference.
To investigate the hydraulic gate vibration, Liang [12] established a two-degrees-of-freedom theoretical model, shown in Figure 6. The model is inspired by the dynamic vibration absorber (DVA), which is commonly employed in the transportation field to control the vibration of vehicles [33]. The two-degrees-of-freedom system consists of two parts: the lightweight accessory structure and the heavyweight primary structure. While the primary structure is under the external excitation P 1 , the vibration component induced by P 1 is transferred from the primary structure to the accessory structure through the simplified stiffness K 2 and damping C 2 , leading to an intensity increase in the accessory structure. This kind of vibration can be defined as associated vibration (AV), and the two-degrees-of-freedom system can be consequently named as the AV system. In practice, the hydraulic gate can be treated as the accessory structure, and the dam can be treated as the primary structure. By analyzing the dynamic responses of the accessory structure in this model, the mechanism of intense vibration that occurred in the surface outlet radial gate of the Jinping arch dam in China has been successfully elucidated.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 8 of 25 The value of added mass is difficult to calculate. In existing studies, several models have been established to calculate the approximate value of the added mass. For large hydraulic structures such as dams, the most classic method is the Westergaard method [31,32]. The added mass can be calculated as follows: where ρ denotes the water density; h denotes the depth of water above the base; z denotes the distance from the surface.
Concluded from practical hydraulic projects [5], the value in Equation (15) is larger than reality. In the engineering field, the added mass is usually considered as follows: With an added mass, the effect of FSI can be taken into consideration, and a more simplified theoretical model of the guide wall can be established.

Theoretical Model
The hydraulic structures, including dams, spillway guide walls, hydraulic gates, are complex multi-degrees-of-freedom (MDOF) systems with complicated connections and damping conditions. For ease of calculation, the complex MDOF system should be simplified appropriately. This paper mainly concerns the vibrational response of the guide wall and the dam so that the system can be transformed into a two-degrees-of-freedom system. It is worth mentioning that the guide wall has a lighter weight compared to the weight of the large-scale concrete dam. To investigate the simplified form of the vibration system in this circumstance, the existing research of the hydraulic gate, another type of lightweight hydraulic structure, can provide a reference.
To investigate the hydraulic gate vibration, Liang [12] established a two-degrees-of-freedom theoretical model, shown in Figure 6. The model is inspired by the dynamic vibration absorber (DVA), which is commonly employed in the transportation field to control the vibration of vehicles [33]. The two-degrees-of-freedom system consists of two parts: the lightweight accessory structure and the heavyweight primary structure. While the primary structure is under the external excitation 1 P , the vibration component induced by 1 P is transferred from the primary structure to the accessory structure through the simplified stiffness 2 K and damping 2 C , leading to an intensity increase in the accessory structure. This kind of vibration can be defined as associated vibration (AV), and the twodegrees-of-freedom system can be consequently named as the AV system. In practice, the hydraulic gate can be treated as the accessory structure, and the dam can be treated as the primary structure. By analyzing the dynamic responses of the accessory structure in this model, the mechanism of intense vibration that occurred in the surface outlet radial gate of the Jinping arch dam in China has been successfully elucidated.  It is noted that the spillway guide wall has a lighter weight than the dam, and can also be considered as an accessory structure connected to the dam. However, the theoretical model for the vibration mechanism of the spillway guide wall is quite different from that for the hydraulic gate vibration mechanism. The reason is that the hydraulic gate is rigidly connected to the dam, while the guide wall is connected to not only the dam but also the foundation. Moreover, the connection between the guide wall and the dam is usually in a flexible form, such as sealing devices, diving joint, and the connection between the guide wall and the foundation is similar to the dam-foundation connection. With regard to these characteristics, the former model can be modified as shown in Figure 7, with the constraint between the guide wall and the foundation taken into account. Note that only horizontal constraints are considered in the theoretical model. The reasons are that in prototype results, the vertical displacements of the dam and the guide wall show less synchronization for structural resonance frequencies, which is analyzed in Section 2. As shown in Figure 8, the vertical displacement of the guide wall is measured as being smaller, which, in practice, is not a major threat to safety [13].
Appl. Sci. 2019, 9, x FOR PEER REVIEW 9 of 25 It is noted that the spillway guide wall has a lighter weight than the dam, and can also be considered as an accessory structure connected to the dam. However, the theoretical model for the vibration mechanism of the spillway guide wall is quite different from that for the hydraulic gate vibration mechanism. The reason is that the hydraulic gate is rigidly connected to the dam, while the guide wall is connected to not only the dam but also the foundation. Moreover, the connection between the guide wall and the dam is usually in a flexible form, such as sealing devices, diving joint, and the connection between the guide wall and the foundation is similar to the dam-foundation connection. With regard to these characteristics, the former model can be modified as shown in Figure  7, with the constraint between the guide wall and the foundation taken into account. Note that only horizontal constraints are considered in the theoretical model. The reasons are that in prototype results, the vertical displacements of the dam and the guide wall show less synchronization for structural resonance frequencies, which is analyzed in Section 2. As shown in Figure 8, the vertical displacement of the guide wall is measured as being smaller, which, in practice, is not a major threat to safety [13]. The EOMs for this modified two-degrees-of-freedom system are as follows [34,35]: where 1 x ,  1 x , and  1 x denote the displacement, velocity, and acceleration of the primary structure, respectively; 2 x ,  2 x , and  2 x denote the displacement, velocity, and acceleration of the accessory structure, respectively; 1 P and 2 P denote the time-variant external exciting force applied to the primary structure and the accessory structure, respectively; 1 M and 2 M denote the coupled mass of the accessory structure and the accessory structure, respectively, which are the sums of the structural and added masses of the corresponding structure, i.e.,  It is noted that the spillway guide wall has a lighter weight than the dam, and can also be considered as an accessory structure connected to the dam. However, the theoretical model for the vibration mechanism of the spillway guide wall is quite different from that for the hydraulic gate vibration mechanism. The reason is that the hydraulic gate is rigidly connected to the dam, while the guide wall is connected to not only the dam but also the foundation. Moreover, the connection between the guide wall and the dam is usually in a flexible form, such as sealing devices, diving joint, and the connection between the guide wall and the foundation is similar to the dam-foundation connection. With regard to these characteristics, the former model can be modified as shown in Figure  7, with the constraint between the guide wall and the foundation taken into account. Note that only horizontal constraints are considered in the theoretical model. The reasons are that in prototype results, the vertical displacements of the dam and the guide wall show less synchronization for structural resonance frequencies, which is analyzed in Section 2. As shown in Figure 8, the vertical displacement of the guide wall is measured as being smaller, which, in practice, is not a major threat to safety [13]. The EOMs for this modified two-degrees-of-freedom system are as follows [34,35]: where 1 x ,  1 x , and  1 x denote the displacement, velocity, and acceleration of the primary structure, respectively; 2 x ,  2 x , and  2 x denote the displacement, velocity, and acceleration of the accessory structure, respectively; 1 P and 2 P denote the time-variant external exciting force applied to the primary structure and the accessory structure, respectively; 1 M and 2 M denote the coupled mass of the accessory structure and the accessory structure, respectively, which are the sums of the structural and added masses of the corresponding structure, i.e.,  The EOMs for this modified two-degrees-of-freedom system are as follows [34,35]: ..
x 1 denote the displacement, velocity, and acceleration of the primary structure, respectively; x 2 , .
x 2 denote the displacement, velocity, and acceleration of the accessory structure, respectively; P 1 and P 2 denote the time-variant external exciting force applied to the primary structure and the accessory structure, respectively; M 1 and M 2 denote the coupled mass of the accessory structure and the accessory structure, respectively, which are the sums of the structural and added masses of the corresponding structure, i.e., M j = M sj + M aj ( j = 1, 2); K 1 and C 1 denote the stiffness and damping of the primary structure, respectively; K 2 and C 2 represent the stiffness and damping between the accessory structure and the primary structure, respectively; K 3 and C 3 represent the stiffness and damping between the accessory structure and the foundation, respectively. As K 2 and C 2 approach zero, the connection to the primary structure tends to disappear; therefore, the accessory structure becomes an independent structure subjected to the forced vibration (FV) only. Similarly, while K 3 and C 3 approach zero, the accessory structure becomes an absolute accessory unit of the dam, and the associated vibration (AV) will be generated. Therefore, the dynamic system shown in Figure 7 can be named as the associated-forced coupled vibration (AFCV) system, and in this AFCV system, K 2 and C 2 represent the factors of AV; K 3 and C 3 represent the factors of FV.
In this study, the exciting forces are flow excitations measured on the fluid-structure interface, which can be Fourier transformed as the sum of a series of harmonic forces in the engineering vibrational field. Hence, P 1 and P 2 can be expressed as harmonic forces [4,5,12,33,36] with forms as below: where F 1 and F 2 denote the amplitudes of corresponding external exciting forces; i is the imaginary unit; ω 0 is the frequency of excitation, which is considered as an equalized value in Equations (18) and (19), since both excitations are induced from supercritical flow.
To investigate the mechanism of guide wall vibration, the amplitude coefficient of the accessory structure is calculated to analyze its dynamic response in the AFCV system, and the aforementioned two kinds of vibrations are taken into account and a comparison is made to explore their contribution to the coupled vibration.

Parameter Sensitivity Analyses
To investigate the influence of every parameter on the amplitude coefficient given in Equation (25), sensitivity analyses for each parameter are performed, and reasonable values are set for all the parameters except the target one beforehand. Referring to the minimum value specified in the code for seismic design of hydraulic structures [37], the damping ratio of the structure made of concrete is considered to be 0.05, i.e., ξ 1 = ξ 2 = ξ 3 = 0.05. The excitation amplitude ratio γ is considered to be 1.5, which is concluded from previous flow excitation studies [13,38,39]. According to the practical structures of TGHS [40] and the added mass of each part of structures, the coupled mass ratio µ is calculated to be 31. The parameter λ 13 is considered to be 0.18, which is deduced from the coupled mass ratio. It is noted that the values of these parameters are only referential values for qualitative studies of each parameter, and the values can be changed to apply this analysis to other practical projects. Based on these predefined values, the dynamic response (amplitude coefficient) with the variation of λ 12 is shown in Figures 9-11.
As shown in Figure 9, the influence of the parameter λ 12 on the dynamic response is investigated for the case of λ 12 ≤ λ 13 = 0.2. Two peaks appear when the external excitation frequency is roughly equal to the natural frequency of the primary structure and the accessory structure in the AV system (i.e., ω 0 = ω 1 and ω 0 = ω 2 ). With the increase of λ 12 , the peak value of ω 0 = ω 1 decreases, while the peak value of ω 0 = ω 2 shows no apparent change. The peak of ω 0 = ω 1 demonstrates that the vibration of the accessory structure synchronizes with the primary structure, i.e., AV is generated; and the peak value of ω 0 = ω 1 represents the intensity of AV. The significant decrease in the amplitude coefficient when ω 0 ω 1 can be explained as the reduction of vibrational synchronization between the vibrations of the primary structure and the accessory structure. It is noticeable that the excitation frequency of the first peak is slightly greater than ω 1 , which is considered to be induced by the constraint between the accessory structure and the foundation.
As shown in Figure 10, while 0.2 < λ 12 ≤ 5, there are still two peaks on the curve, but the peak corresponding to ω 0 = ω 1 tends to disappear with the increase of λ 12 . Different from Figure 9, the second peak no longer occurs when ω 0 = ω 2 , but is generated when ω 0 = ω 3 , and its value increases with an increasing λ 12 . This phenomenon can be illustrated as the recession of AV: as λ 12 increases, which means the frequency ω 2 gradually approaches and exceeds ω 1 , the vibration synchronization is weakened, causing a further decrease in the vibration induced by the primary structural vibration. Simultaneously, with a larger ω 2 and a smaller ω 3 , the FV of the accessory structure subjected to external excitation becomes more significant so that the peak value of ω 0 = ω 3 increases with the increase in λ 12 . ω ω demonstrates that the vibration of the accessory structure synchronizes with the primary structure, i.e., AV is generated; and the peak value of = 0 1 ω ω represents the intensity of AV. The significant decrease in the amplitude coefficient when ≠ 0 1 ω ω can be explained as the reduction of vibrational synchronization between the vibrations of the primary structure and the accessory structure. It is noticeable that the excitation frequency of the first peak is slightly greater than 1 ω , which is considered to be induced by the constraint between the accessory structure and the foundation.  As shown in Figure 10, while . , and its value increases with an increasing 12 λ . This phenomenon can be illustrated as the recession of AV: as 12 λ increases, which means the frequency 2 ω gradually approaches and exceeds 1 ω , the vibration synchronization is weakened, causing a further decrease in the vibration induced by the primary structural vibration. Simultaneously, with a larger 2 ω and a smaller 3 ω , the FV of the accessory structure subjected to external excitation becomes more significant so that the peak value of 0 increases with the increase in 12 λ . Another situation where 12 5 λ > is shown in Figure 11. In this situation, the curve becomes single-peaked. The peak of 0 1 ω ω = mentioned above has entirely disappeared, which indicates that the vibration of the accessory structure has totally separated from that of the primary structure and ω ω demonstrates that the vibration of the accessory structure synchronizes with the primary structure, i.e., AV is generated; and the peak value of = 0 1 ω ω represents the intensity of AV. The significant decrease in the amplitude coefficient when ≠ 0 1 ω ω can be explained as the reduction of vibrational synchronization between the vibrations of the primary structure and the accessory structure. It is noticeable that the excitation frequency of the first peak is slightly greater than 1 ω , which is considered to be induced by the constraint between the accessory structure and the foundation.  As shown in Figure 10, while . , and its value increases with an increasing 12 λ . This phenomenon can be illustrated as the recession of AV: as 12 λ increases, which means the frequency 2 ω gradually approaches and exceeds 1 ω , the vibration synchronization is weakened, causing a further decrease in the vibration induced by the primary structural vibration. Simultaneously, with a larger 2 ω and a smaller 3 ω , the FV of the accessory structure subjected to external excitation becomes more significant so that the peak value of 0 increases with the increase in 12 λ . Another situation where 12 5 λ > is shown in Figure 11. In this situation, the curve becomes single-peaked. The peak of 0 1 ω ω = mentioned above has entirely disappeared, which indicates that the vibration of the accessory structure has totally separated from that of the primary structure and By the discrimination of the above three situations, the dynamic response with a various parameter has been adequately analyzed. Generally speaking, with the increase of 12 λ , AV decreases until it disappears. As 01 λ increases, the first peak roughly matches 0 1 ω ω = and the peak value significantly decreases with a rising 12 λ , and disappears when 12 5 λ > . Moreover, the second peak initially corresponds to = ω ω , but when the value of λ approaches that of λ , the abscissa of this Another situation where λ 12 > 5 is shown in Figure 11. In this situation, the curve becomes single-peaked. The peak of ω 0 = ω 1 mentioned above has entirely disappeared, which indicates that the vibration of the accessory structure has totally separated from that of the primary structure and becomes an independent FV. Moreover, the peak value of ω 0 = ω 3 marginally rises with a larger λ 12 .
By the discrimination of the above three situations, the dynamic response with a various parameter has been adequately analyzed. Generally speaking, with the increase of λ 12 , AV decreases until it disappears. As λ 01 increases, the first peak roughly matches ω 0 = ω 1 and the peak value significantly decreases with a rising λ 12 , and disappears when λ 12 > 5. Moreover, the second peak initially corresponds to ω 0 = ω 2 , but when the value of λ 12 approaches that of λ 13 , the abscissa of this peak remains almost unchanged with the variation of λ 12 and gradually approaches ω 0 = ω 3 .
As illustrated in Figure 12, the influence of the parameter λ 13 on the dynamic response is investigated. Figure 12a shows the dynamic response of the accessory structure in the case of λ 12 = 0.5. It is noteworthy that the peak value is remarkably increased with an increasing λ 13 in the case of ω 0 = ω 1 , denoting the dominance of AV; the other peak originally corresponds to ω 0 = ω 3 while λ 13 < λ 12 = 0.5, but becomes asymptotic to ω 0 = ω 2 as λ 13 gets larger. Moreover, Figure 12b,c are detailed to obtain more specific analyses for the influence of λ 13 on the dynamic response with different λ 12 . It is evident that similar conclusions can be drawn from Figure 12a,b. For the one-peak situation (Figure 12c), the peak of ω 0 = ω 3 abides by the same law as that in Figure 12a, before eventually becoming asymptotic to ω 0 = ω 2 when λ 13 > λ 12 . The phenomenon above can be interpreted as the growth of AV with the increase of λ 13 . When λ 13 is increased, which means that the natural frequency of the accessory structure in the independent FV system (as an independent structure) has been decreased. It can be inferred that the influence of AV will become salient. As illustrated in Figure 13, the influence of the damping ratio 1 ξ on the dynamic response is investigated. It can be discovered that the increase of 1 ξ suppresses the value of the first peak, but has the opposite effect on the value of the second peak. According to Figure 14, the value of the second peak decreases with the increase of the damping ratio 2 ξ , but the value of the first peak has not been Combined with the analyses of λ 12 , it can be concluded that λ 12 is a parameter determining the number of peaks, and λ 13 is a parameter determining peak values. Both λ 12 and λ 13 determine the vibration type. AV is dominant when λ 12 < λ 13 ; as λ 12 /λ 13 gets larger, the vibration gradually becomes the independent FV. The dividing point is λ 12 = λ 13 from the curves. Based on the above analyses, the parameter λ 12 is considered to be 0.5 in subsequent analyses to take the effects of both AV and independent FV into consideration, which can be seen from Figure 10.
As illustrated in Figure 13, the influence of the damping ratio ξ 1 on the dynamic response is investigated. It can be discovered that the increase of ξ 1 suppresses the value of the first peak, but has the opposite effect on the value of the second peak. According to Figure 14, the value of the second peak decreases with the increase of the damping ratio ξ 2 , but the value of the first peak has not been influenced. Figure 15 illustrates the influence of the damping ratio ξ 3 on the dynamic response, which is similar to that of ξ 2 , but the decrease influence of ξ 3 on the second peak is greater. In conclusion, the three damping ratios change the peak values, but they do not change the dominant vibration type. As illustrated in Figure 13, the influence of the damping ratio 1 ξ on the dynamic response is investigated. It can be discovered that the increase of 1 ξ suppresses the value of the first peak, but has the opposite effect on the value of the second peak. According to Figure 14, the value of the second peak decreases with the increase of the damping ratio 2 ξ , but the value of the first peak has not been influenced. Figure 15 illustrates the influence of the damping ratio 3 ξ on the dynamic response, which is similar to that of 2 ξ , but the decrease influence of 3 ξ on the second peak is greater. In conclusion, the three damping ratios change the peak values, but they do not change the dominant vibration type. Figure 13. The dynamic response of the accessory structure upon varying the parameter 1 ξ . Figure 13. The dynamic response of the accessory structure upon varying the parameter ξ 1 .
Appl. Sci. 2019, 9, x FOR PEER REVIEW 15 of 25 Figure 14. The dynamic response of the accessory structure upon varying the parameter 2 ξ . Figure 14. The dynamic response of the accessory structure upon varying the parameter ξ 2 .
As illustrated in Figure 16, as the parameter γ increases, the first peak value is significantly increased, while the second peak value is slightly reduced. The changes reflect the transformation of the vibration type. According to Figure 17, the increase of the parameter µ can result in a decrease in the value of first peak. Moreover, it causes a marginal increase in the second peak value but has no effect on the vibration type.  As illustrated in Figure 16, as the parameter γ increases, the first peak value is significantly increased, while the second peak value is slightly reduced. The changes reflect the transformation of the vibration type. According to Figure 17, the increase of the parameter μ can result in a decrease in the value of first peak. Moreover, it causes a marginal increase in the second peak value but has no effect on the vibration type.   As illustrated in Figure 16, as the parameter γ increases, the first peak value is significantly increased, while the second peak value is slightly reduced. The changes reflect the transformation of the vibration type. According to Figure 17, the increase of the parameter μ can result in a decrease in the value of first peak. Moreover, it causes a marginal increase in the second peak value but has no effect on the vibration type.

Vibration Correlation Analyses by XWT Method
As recommended in the introduction, cross-wavelet transform (XWT) is an efficient method to analyze the correlation between two signals. Suppose the wavelet coefficients

Vibration Correlation Analyses by XWT Method
As recommended in the introduction, cross-wavelet transform (XWT) is an efficient method to analyze the correlation between two signals. Suppose the wavelet coefficients W y (a, τ) of a finite time signal y(t) with where a represents the scaling parameter; t represents time; τ represents the location of t; ψ(t) denotes the mother wavelet; and ψ (t) extends the mother wavelet to a wavelet family. In practice, the mother wavelet ψ(t) is always considered as a Morlet wavelet: where ω is the nondimensional frequency, which is considered to be six to satisfy the admissibility condition [41]. For two signals, y(t) and z(t), the XWT is defined as follows: the asterisk * denotes the complex conjugate. In engineering practice, the XWT is always transferred to the form of the cross-wavelet scalogram (XWS) to analyze the non-stationary signal: XWS y,z (a, τ) = XW y,z (a, τ)XW y,z * (a, τ) = XW y,z (a, τ) 2 By defining XW y,z (a, τ) and XWS y,z (a, τ), the correlation between two signals in both time and frequency domains is revealed, and the value of XWS y,z (a, τ) represents the correlation level.
Applying the XWT technique for signals obtained from 1# dam block and eight guide wall blocks, the correlation level in the form of XWS is illustrated in Figure 18. The figures show that the vibration of the dam block and guide wall blocks demonstrate clear synchronized correlation for the frequency band of about 2.5-7 Hz. It can be interpreted that under these synchronized vibration frequencies, AV tends to occur. For the 1# and 2# guide wall block, the correlations of both time and frequency domains are weaker, and the dominant synchronized frequency is located around 2.8 Hz, which may represent the natural frequency of the large-sized dam. It can be deduced that the vibration of these guide wall blocks may derive from the dam block, i.e., AV is dominant in these guide wall blocks. Meanwhile, in downstream guide wall blocks (3#-8# blocks), the correlative frequency bands are broader, and the dominant synchronized frequency is gradually increased to about 6.2 Hz, which represents the natural frequency of the small-sized guide wall. The natural frequency of the dam becomes less evident, from which it can be inferred that the influence of AV is attenuated in downstream guide wall blocks. The guide wall vibration transforms into independent FV, and the broad correlative frequency bands are derived from the broad frequency band of flow fluctuation that acts on both the dam and the guide wall.
Meanwhile, in downstream guide wall blocks (3#-8# blocks), the correlative frequency bands are broader, and the dominant synchronized frequency is gradually increased to about 6.2 Hz, which represents the natural frequency of the small-sized guide wall. The natural frequency of the dam becomes less evident, from which it can be inferred that the influence of AV is attenuated in downstream guide wall blocks. The guide wall vibration transforms into independent FV, and the broad correlative frequency bands are derived from the broad frequency band of flow fluctuation that acts on both the dam and the guide wall. Combined with the theoretical model of the AFCV system in Section 3, the above phenomena can be elucidated. Generally speaking, the damping ratio of a built-up structure is fixed, which can infer that parameters 1 ξ and 3 ξ are constant; each guide wall block can be considered as being the same structures, therefore, the parameters 13 λ and μ are fixed in this case. For blocks which are located more downstream, the connecting effect to the dam is weakened, and the natural frequency of the accessory Combined with the theoretical model of the AFCV system in Section 3, the above phenomena can be elucidated. Generally speaking, the damping ratio of a built-up structure is fixed, which can infer that parameters ξ 1 and ξ 3 are constant; each guide wall block can be considered as being the same structures, therefore, the parameters λ 13 and µ are fixed in this case. For blocks which are located more downstream, the connecting effect to the dam is weakened, and the natural frequency of the accessory structure in the AV system has the same tendency, leading to an increase in λ 12 ; the longer propagation path facilitates the attenuation of vibration, which means an increase in ξ 2 ; the downstream blocks are located at the dissipation zone, and contact with more complex flow conditions and more severe pulsation, causing γ to decrease. Summarized from Figures 9-11, 14, and 16, it can be discovered that the above changes of λ 12 and γ sharply decline the amplitude coefficient for ω 0 = ω 1 , which is equivalent to the attenuation of the AV intensity; the change in ξ 2 increases the amplitude coefficient for ω 0 = ω 2 , which means that the independent FV is enhanced. As a result, the XWS patterns are more concentrated on lower frequencies, as shown in Figure 18a-c; the value of XWS in downstream blocks is remarkably higher than that in upstream blocks, as shown in Figure 18d-h. Hence, the actual vibrational circumstances validate the correctness of the theoretical model in Section 3.

Numerical Simulation
Based on the above analysis, it can be concluded that AV can be detected in the guide wall structure, and it indeed amplifies the vibrations of certain blocks on the guide wall during the discharge process. To explore the detailed mechanism of guide wall vibration, numerical simulations for the vibration mode of the dam and the guide wall are proposed using the commercial software ANSYS.
Compared with the operational modal parameter identification by eigensystem realization algorithm (ERA) on the same project [42], it is shown that the vibration appears as the vibration of an integral structure, rather than the independent vibration of each block, due to the transmission of shear forces through joints and thermal expansion (the tests were implemented in the summer). Therefore, the cementation between adjacent blocks is considered, with the integral numerical models built, and the frequencies of the dam and the guide wall are calculated for both dry and wet modes.
The numerical models and the corresponding material parameters (density ρ, elastic modulus E, and Poisson ratio ν) are shown in Figure 19. For the guide wall, the connection to the foundation is simplified as the all-degrees of freedom (all-DOF) constraint applied to the bottom, which corresponds to the constraints in the independent FV system; the connection to the dam is simplified as the all-DOF without the longitudinal displacement constraint for the diving joint between the guide wall and the dam, corresponding to the constraints in the AV system. For the AFCV system, both kinds of constraints are applied. For the dam, the all-DOF constraint is applied to the bottom, and the diving joint between the guide wall and dam is considered as the all-DOF without the longitudinal displacement constraint. For the wet mode, the range of water is considered to be twice as much as the height of the guide wall or the dam. guide wall and the dam, corresponding to the constraints in the AV system. For the AFCV system, both kinds of constraints are applied. For the dam, the all-DOF constraint is applied to the bottom, and the diving joint between the guide wall and dam is considered as the all-DOF without the longitudinal displacement constraint. For the wet mode, the range of water is considered to be twice as much as the height of the guide wall or the dam. Table 2 shows several natural frequencies of the dam for both dry and wet modes. There are several modes with frequencies that match actual dominant frequencies. It is noted that the dam is a structure with closely spaced natural frequencies, whether in dry or wet mode. That is, under the flow fluctuation with a broad frequency band, there is always a mode with a natural frequency which is close to the excitation frequency.    Table 2 shows several natural frequencies of the dam for both dry and wet modes. There are several modes with frequencies that match actual dominant frequencies. It is noted that the dam is a structure with closely spaced natural frequencies, whether in dry or wet mode. That is, under the flow fluctuation with a broad frequency band, there is always a mode with a natural frequency which is close to the excitation frequency. As mentioned above, by changing constraints, the vibration modes for the AV system, independent FV system and AFCV system can be respectively analyzed. Table 3 illustrates some typical natural frequencies of the guide wall. There are also several modes with frequencies that match the actual dominant frequencies. According to prototype tests, the root mean square (RMS) of the horizontal displacement reaches a maximum in the 6# guide wall block. Through numerical simulation, it can be found that for all three vibration systems, there are always several natural frequencies in the range of 2-7 Hz that are in accordance with actual vibrational characteristics. The corresponding vibration patterns, their orders, and natural frequencies are illustrated in Figure 20. Therefore, these vibration modes can be inferred to participate in the actual vibration.
According to prototype tests, the root mean square (RMS) of the horizontal displacement reaches a maximum in the 6# guide wall block. Through numerical simulation, it can be found that for all three vibration systems, there are always several natural frequencies in the range of 2-7 Hz that are in accordance with actual vibrational characteristics. The corresponding vibration patterns, their orders, and natural frequencies are illustrated in Figure 20. Therefore, these vibration modes can be inferred to participate in the actual vibration.

Vibration Mechanism Analyses for the TGHS Left Guide Wall
Based on the spectrum analyses, XWT correlation analyses, and numerical simulation, the detailed mechanism of the TGHS left guide wall vibration can be clarified with the proposed AFCV theory. The guide wall vibration consists of the AV induced by the dam vibration and the independent FV caused by flow fluctuation. The contribution of AV depends on the location of the guide wall block. Using the theoretical model detailed in Section 3, the mechanism can be verified as below.
According to the modal analyses, the dam is a structure with closely spaced modes, therefore, there is always a natural frequency of the dam that matches the excitation frequency. The concerned frequencies for all of the corresponding vibration systems have participated in the vibration, which means they also match the excitation frequencies. Referring to the corresponding natural frequencies in Figure 20, the parameter 12 λ can be considered as 0.8, and 13 λ can be considered as 1. Consequently, the dynamic response of the guide wall is shown in Figure 21.
For this case, parameters 2 ξ and γ are considered to be corresponding parameter values of the integral guide wall, as well as the parameter 12 λ , as mentioned above. For upstream guide wall blocks, which are further from the dissipation zone and closer to the dam, the flow fluctuations are weak, and the frequency bands are narrow. As mentioned above, the external excitation matches the natural frequency of the dam; therefore, 01 λ is located around 1. As shown in Figure 21, there is a peak value at 01 1 λ = , indicating that the guide wall vibration is significantly affected by the dam, i.e., AV is significant, corresponding to Figure 18a-c. While, for downstream guide wall blocks, which are closer to the dissipation zone, the flow here demonstrates strong fluctuation and a broad frequency bandwidth. The ratio between the frequencies of external excitation and the dam (i.e., 01 λ ) range from 0.2 to 5 (analogous to Figure 10). It is noted that the two peaks are located in this range, which

Vibration Mechanism Analyses for the TGHS Left Guide Wall
Based on the spectrum analyses, XWT correlation analyses, and numerical simulation, the detailed mechanism of the TGHS left guide wall vibration can be clarified with the proposed AFCV theory. The guide wall vibration consists of the AV induced by the dam vibration and the independent FV caused by flow fluctuation. The contribution of AV depends on the location of the guide wall block. Using the theoretical model detailed in Section 3, the mechanism can be verified as below.
According to the modal analyses, the dam is a structure with closely spaced modes, therefore, there is always a natural frequency of the dam that matches the excitation frequency. The concerned frequencies for all of the corresponding vibration systems have participated in the vibration, which means they also match the excitation frequencies. Referring to the corresponding natural frequencies in Figure 20, the parameter λ 12 can be considered as 0.8, and λ 13 can be considered as 1. Consequently, the dynamic response of the guide wall is shown in Figure 21. , the contribution from independent FV is more significant than that from AV.
Because there are also several natural frequencies of the dam which are in this range, it illustrates a broader correlation between the dam block and the guide wall blocks, as seen in Figure 18d-h. Therefore, by the verification of prototype tests, it can be demonstrated that the vibration mechanism of the TGHS left guide wall can be distinctly clarified by the theoretical model, in other words, the proposed AFCV theory is valid and correct. For this case, parameters ξ 2 and γ are considered to be corresponding parameter values of the integral guide wall, as well as the parameter λ 12 , as mentioned above. For upstream guide wall blocks, which are further from the dissipation zone and closer to the dam, the flow fluctuations are weak, and the frequency bands are narrow. As mentioned above, the external excitation matches the natural frequency of the dam; therefore, λ 01 is located around 1. As shown in Figure 21, there is a peak value at λ 01 = 1, indicating that the guide wall vibration is significantly affected by the dam, i.e., AV is significant, corresponding to Figure 18a-c. While, for downstream guide wall blocks, which are closer to the dissipation zone, the flow here demonstrates strong fluctuation and a broad frequency bandwidth. The ratio between the frequencies of external excitation and the dam (i.e., λ 01 ) range from 0.2 to 5 (analogous to Figure 10). It is noted that the two peaks are located in this range, which indicates that both AV and independent FV contribute to the guide wall vibration. Due to the greater peak value of ω 0 = ω 2 , the contribution from independent FV is more significant than that from AV. Because there are also several natural frequencies of the dam which are in this range, it illustrates a broader correlation between the dam block and the guide wall blocks, as seen in Figure 18d-h.

Conclusions
Therefore, by the verification of prototype tests, it can be demonstrated that the vibration mechanism of the TGHS left guide wall can be distinctly clarified by the theoretical model, in other words, the proposed AFCV theory is valid and correct.

Conclusions
In this paper, a new vibration mechanism is presented to analyze the flow-induced vibration (FIV) of the spillway guide wall in a hydropower station, with the new concepts of associated vibration (AV) and associated-forced coupled vibration (AFCV) defined. The prototype tests on the left guide wall of the Three Gorges Hydropower Station (TGHS) indicate that there are several similar dominant frequencies of the horizontal displacements between 1# dam block and guide wall blocks. To investigate the mechanism of the vibration synchronization, a two-degrees-of-freedom theoretical model of the AFCV system is proposed, considering the lightweight spillway guide wall as an accessory structure installed in both the large-scale primary structure (the dam) and the foundation. Afterwards, the practicability of this theoretical model is verified by mutual analyses. Using the cross-wavelet transform (XWT), the correlation between the vibrational signals measured on the dam block and guide wall blocks, respectively, is illustrated, which confirms the existence of AFCV. With the vibrational characteristics calculated by numerical simulation, the vibration mechanism is clarified by the theoretical model.
In brief, the guide wall vibration is not only a simple independent forced vibration (FV) but also a vibration influenced by the dam, of which the effect can be named as AFCV. This kind of vibration derives from the effect of the shear force through joints and thermal expansion. The contribution level of the AV on the guide wall depends on the relative position of the guide wall, the dam, and the position of the dissipation zone. When the guide wall is positioned closer from the dam and further to the dissipation zone, the guide wall vibration appears to possess synchronized characteristics with the dam, proving the existence of AV; the synchronized vibrational frequency band is narrow but concentrated, which indicates AV is a dominant vibration component in this case. On the contrary, when the guide wall is positioned further from the dam and closer to the dissipation zone, the guide wall vibration appears to be with a broad frequency band and possesses more independent FV characteristics, although the vibration intensity tends to increase because of the stronger fluctuation in the flow.
It is noted that all of the above changes in the vibrational characteristics can be elucidated by varying the parameters in the theoretical model, such as λ 12 , ξ 2 , γ, which indicates the validity of the theoretical model. Furthermore, by changing parameters using practical engineering measures, the model can provide a certain beneficial reference regarding vibration reduction.
Moreover, it should be pointed out that due to the complexity of the engineering structures and FSI effect, the mechanism of vibration propagation and the engineering measures to reduce propagation still require further investigation.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A the theoretical model. Furthermore, by changing parameters using practical engineering measures, the model can provide a certain beneficial reference regarding vibration reduction. Moreover, it should be pointed out that due to the complexity of the engineering structures and FSI effect, the mechanism of vibration propagation and the engineering measures to reduce propagation still require further investigation. the theoretical model. Furthermore, by changing parameters using practical engineering measures, the model can provide a certain beneficial reference regarding vibration reduction. Moreover, it should be pointed out that due to the complexity of the engineering structures and FSI effect, the mechanism of vibration propagation and the engineering measures to reduce propagation still require further investigation.