Next Article in Journal
Analysis of Discrete-Time Queues with Branching Arrivals
Next Article in Special Issue
Ordinary Kriging Interpolation Method Combined with FEM for Arch Dam Deformation Field Estimation
Previous Article in Journal
Polynomial Chaos Expansion-Based Enhanced Gaussian Process Regression for Wind Velocity Field Estimation from Aircraft-Derived Data
Previous Article in Special Issue
The Optical Path Method for the Problem of Oblique Incidence of a Plane Electromagnetic Wave on a Plane-Parallel Scatterer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structural Modal Parameter Identification Method Based on the Delayed Transfer Rate Function under Periodic Excitations

1
State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210024, China
2
College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing 210024, China
3
National Engineering Research Center of Water Resources Efficient Utilization and Engineering Safety, Hohai University, Nanjing 210024, China
4
College of Civil and Transportation Engineering, Hohai University, Nanjing 210024, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(4), 1019; https://doi.org/10.3390/math11041019
Submission received: 16 January 2023 / Revised: 11 February 2023 / Accepted: 14 February 2023 / Published: 16 February 2023
(This article belongs to the Special Issue Mathematical Modeling and Numerical Analysis for Applied Sciences)

Abstract

:
The dynamic response transfer rate function (TRF) is increasingly used in the field of structural modal parameter identification because it does not depend on the white noise assumption of the excitation. In this paper, the interference of periodic excitation on structural modal parameter identification using TRF is analyzed theoretically for a class of civil engineering structures with obvious periodic components in excitation, and then an identification method of structural real modal parameters is proposed. First, a delayed TRF is constructed, and the pseudo-frequency response function is further obtained to identify the periodic spurious poles of the whole system. Then, the effective identification of the real modal parameters of the structure is achieved by comparing the system poles identified via conventional TRF. Finally, the feasibility and robustness of the proposed method were verified using a calculation example with four-degrees-of-freedom system. In addition, the modal parameters of a structure under periodic excitation were effectively identified by taking a pumping station as an example, and the results show that the method accurately identified the structural modal parameters when the excitation contained periodic components, which has wider prospects for technical applications.

1. Introduction

In large civil engineering projects, the modal parameters of a structure are important pieces of information reflecting the dynamic characteristics of the structure, and accurate recording of the modal parameters is a prerequisite for effective structural health evaluation and damage diagnosis [1,2,3]. The operating modal analysis method can effectively obtain the modal parameters of a structure without affecting the normal operation of the structure, based only on the dynamic response of the structure under the environmental excitation [4,5]. Its low cost of operation and harmlessness to structures have led to its development as an effective tool for structural modal analysis of large civil engineering projects [6,7,8]. It is worth noting that the modal analysis method based on the transfer rate function (TRF) is an emerging method of structural operational modal analysis developed in recent decades which does not depend on the white noise excitation assumption of the load, greatly expands the scope of engineering applications of the operational modal analysis, and has received wide attention [9,10,11]. The commonly used TRFs take the form of response transmissibility [12,13] and power spectral density transmissibility of the structure [14]. According to the nature of the TRF intersection at the system poles, the identification of the modal parameters can be achieved by constructing the quasi-frequency response function matrix of the structure [15].
As mentioned above, the TRF has special properties. Devriendt et al. [16] successfully identified the modal parameters of a pedestrian bridge, extended the scalar transmissibility to the transmissibility matrix, and proposed a multivariate transmissibility [17] and pseudo-transmissibility [18] method for operational modal analysis. Weitjens et al. [19] proposed a transmissibility operational modal analysis method for multiple-input–multiple-output systems using multi-reference-point least-squares complex frequency domain estimation to solve for the system poles and derive the system modal parameters based on the steady-state diagram. Zhang et al. [20] constructed rational functions by applying transmissibility for two different loading cases and identified the modal parameters using orthogonal polynomial fitting. Araujo et al. [21] used a combination of transmissibility and blind source separation techniques to analyze structural modes for a low signal-to-noise ratio, mode-dense case. Li et al. [22] applied the vibration response transmissibility to analyze structural operating modes. It can be seen that the TRF has been widely used in the modal identification of various civil engineering structures.
In industrial plants, hydraulic (thermal) powerhouses, and other types of engineering structure, the ambient excitation contains an obvious cyclical component that is subject to the influence of rotating mechanical equipment during operation. Conventional modal analysis methods for TRFs do not account for the presence of significant periodic components in the ambient excitation [23], making it difficult to effectively identify false modes introduced by periodic excitation. Simultaneously, the periodic load components of such structures have a high energy proportion, which tends to obscure the real modal information of the structure [24]. In addition, there are some methods to identify false modes, such as sparse decomposition techniques and statistical methods [25], methods based on fuzzy clustering [26], and modal order-determination indexes [27]. Although some attempts have been made to address this issue, these methods still struggle to accurately identify the false modes introduced by periodic excitation [28].
For this purpose, based on analysis of the mechanism of the influence of periodic excitation on modal identification, this paper proposes identification of the false modes introduced by periodic excitation by establishing a delayed TRF. Simultaneously, the conventional TRF modal recognition method is combined to identify the real poles of the system to effectively obtain the real modal parameters of the structure. Finally, it was verified that the method is effective in identifying the real modal parameters of a structure when the ambient excitation contains a periodic component, and that the method has good applicability at various amplitudes and frequencies of periodic excitation, using a technical calculation example of a four-degrees-of-freedom (4-DOF) system.
The remainder of this paper is organized as follows. The basic theories of modal parameter identification, the periodic excitation effect, and the delayed TRF method are presented in Section 2. The example analysis of this study, including the numerical calculation of a 4-DOF system and modal parameter identification of a pumping station, is described in Section 3. Finally, the conclusions are given in Section 4.

2. Basic Theories

2.1. Modal Parameter Identification Based on Acceleration TRF

For a viscous damping system with N degrees of freedom, the forced vibration equation is written as follows:
M y ¨ t + C y ˙ t + K y t = f t
where M, C, and K are the matrices of mass, damping, and stiffness, respectively; f t is the excitation vector; and y t is the displacement response vector.
Under the condition that the initial system states are 0 and the system excitation is smooth, the Laplace transformation of Equation (1) yields
Y s = H s F s = k = 1 Q H k s F k s
where Y s and F s are the Laplace transformations of y t and f t , respectively; s is the Laplace variable; H s is the transfer function matrix of the system, H s = s 2 M + s C + K 1 ; Q is the number of load degrees of freedom; and H k s is the row vector in H s corresponding to the load degrees of freedom k.
The acceleration response TRF t i j A s is defined as the ratio of the acceleration response of degrees of freedom i and j in the Laplace domain. It can be proven [18] that under the condition that initial displacement and initial velocity are zero, the acceleration transmissibility of the structure is equal to the displacement transmissibility, that is
t i j A s = A i s A j s = t i j D s = Y i s Y j s = k = 1 Q H i k s F k s k = 1 Q H j k s F k s
where Y i s and A i s denote the Laplace transformations of the displacement response y i t and acceleration response a i t of DOF i, respectively. H i k s denotes the i-th row and k-th column element of H s .
For the sake of simplicity, the transmissibility of the acceleration response is denoted by t i j s in the following. From the form of Equation (3), it can be stated that the TRF t i j s under multipoint excitation conditions is related to both the system properties H s and the excitation F s . Therefore, t i j s differs for different loading conditions.
According to the theory of modal analysis and the residue theorem, the following relationship can be obtained:
lim s λ m s λ m H s = ψ m l m T
where λ m is the pole of the m-th order mode; ψ m and l m T are the vibration vector and the modal participation vector corresponding to λ m , respectively.
From Equations (3) and (4), the following equation can be obtained when s λ m :
lim s λ m t i j s = lim s λ m k = 1 Q s λ m H i k s F k s lim s λ m k = 1 Q s λ m H j k s F k s = k = 1 Q ψ m , i l m , k F k λ m k = 1 Q ψ m , j l m , k F k λ m = ψ m , i ψ m , j k = 1 Q l m , k F k λ m k = 1 Q l m , k F k λ m = ψ m , i ψ m , j
From Equation (5), it can be seen that for any load excitation condition, the TRF t i j s of the system acceleration response at the natural frequency pole λ m always remains the same. Based on this property, the identification of structural modal parameters can be carried out.
The acceleration response sequence is divided into N f periods and the TRF t i j f k s under each period is calculated separately. By choosing a specific output DOF j and a time period f s , the matrix T H s can be constructed as follows:
T H s = 1 i = 1 , i j N 0 f k = 1 , f k f s N f t i j f k s t i j f s s t 1 j 1 s t 1 j 2 s t 1 j N f s t 2 j 1 s t 2 j 2 s t 2 j N f s t N 0 j 1 s t N 0 j 2 s t N 0 j N f s
where the superscript of t i j f k s indicates a different time period.
From Equation (5), the denominator of T H s tends to 0 when s λ m . That is, the system pole λ m is also the pole of T H s . It can be proven [29] that T H s is equivalent to the transfer function matrix H s of the system in the neighborhood of the system poles, and T H s is said to be the matrix of the proposed frequency response function. Therefore, conventional frequency domain class methods, such as frequency and spatial domain decomposition (FSDD) [30], can be applied to identify the modal parameters, including the system intrinsic frequency, modal vibration pattern, and damping ratio.
First, the singular value decomposition of T H ( j ω ) at each discrete frequency point ω = ω i is performed as follows:
T H ( j ω i ) = U i S i U i H
where U i is the matrix of singular value vectors; S i is a diagonal matrix of singular values.
The enhanced proposed frequency response function can be constructed as follows:
T H e n ( j ω i ) = Re [ u i H T H ( j ω i ) u i ]
where u i denotes the first singular value vector in U i . The T H e n ( j ω i ) obtained at this point is the scalar. At the m-th order mode, the following exist:
T H e n ( j ω ) ω ω m = 2 Re q m j ω λ m = 2 q m σ m σ m 2 + ω m 2 + ω 2 2 ω ω m
For all frequency points ( ω k ( k = 1 , 2 , , L ) ) near the natural frequency ω m , σ m and ω m are obtained by least-squares fitting, and then the modal parameters are calculated as follows:
f m = σ m 2 + ω m 2 ξ m = σ m σ m 2 + ω m 2
After the value of the natural frequency f m is obtained, the corresponding mode vibration pattern u m can be obtained from the matrix of singular value vectors U i ( ω i = 2 π f m ) determined by Equation (7).

2.2. The Effect of Periodic Excitation on Modal Identification

The environmental excitation acting on large civil engineering structures is commonly assumed to be random white noise excitation since it cannot be measured directly, and this is referred to as the Class A condition. Under some special conditions, the structure’s ambient excitation also contains a periodic component, referred to as the Class B condition. The input–output relationship of the system in the frequency domain for these two cases is shown in Figure 1.
The input–output relationship of the structural system under the Class A condition is expressed as follows:
Y ( s ) = H ( s ) E
In classical control theory [31], any colored noise signal F k ( s ) can be considered as the output after feeding ideal white noise into a specific filter G k ( s ) , that is
F k ( s ) = G k ( s ) E k
For a structural system in the Class B condition, the input–output relationship is expressed as follows:
Y ( s ) = H ( s ) F ( s ) = H ˜ ( s ) E
where E is a complex constant vector consisting of E N i , and E = E 1 , E 2 , , E N i T .
Thus, the structure in the Class B condition can be regarded as an equivalent system H ˜ ( s ) , consisting of the filter G ( s ) and the original system H ( s ) :
H ˜ ( s ) = H ( s ) G ( s )
where G ( s ) is the diagonal matrix composed of G N i ( s ) , and G ( s ) = d i a g G 1 ( s ) , G 2 ( s ) , , G N i ( s ) .
Therefore, the inherent modal information of the structure can be identified using the dynamic response in the Class A condition when the load information is unknown, while the identification result of the Class B condition additionally contains information about the periodic excitation. The effect of periodic excitation on the modal analysis of TRF operation is demonstrated by the following theoretical derivation.
For an engineering structure under the Class B condition, the ambient excitation in the Laplace domain can be expressed as follows:
F s = P s + R s
where P s is periodic excitation and R s is random excitation.
As the system excitation is smooth, the spectral characteristics of the periodic and random excitations in a finite time period l can be expressed as follows:
P l s = a r , s = γ r 0 , s = o t h e r w i s e
R l s = r a n d μ s , σ s 2
where γ r is the rth frequency in the periodic excitation; a r is the amplitude vector of the periodic excitation; μ s and σ s 2 are the mean and variance of the random excitation amplitude at s, respectively.
When there is a significant periodic component in the dynamic response, the amplitude of the periodic component of the structural excitation is much larger than that of the random excitation component. The following equation is then available for the load degree of freedom k:
R k l γ r a r k = 0
According to Equation (3), the TRF of acceleration in any finite length time slot l can be obtained as follows:
t i j l s = k = 1 Q H i k s F k l s k = 1 Q H j k s F k l s
From Equations (15)–(19), when s γ r , the following exist:
lim s γ r t i j l s = lim s γ r k = 1 Q H i k s F k l s lim s γ r k = 1 Q H j k s F k l s = k = 1 Q H i k γ r P k l γ r + R k l γ r k = 1 Q H j k γ r P k l γ r + R k l γ r   = k = 1 Q H i k γ r a r k + o a r k k = 1 Q H j k γ r a r k + o a r k k = 1 Q a r k H i k γ r k = 1 Q a r k H j k γ r
The final expression of Equation (20) is a constant. Therefore, when the external load consists of random excitation and stronger periodic excitation, it is known from Equations (5) and (20) that the TRF t i j l s at any time period l tends towards a constant value at the system pole λ m and at the periodic excitation frequency γ r .
In addition, from Equations (15)–(19), the following equation is obtained when s λ m , γ r :
t i j l s = k = 1 Q H i k s F k l s k = 1 Q H j k s F k l s = k = 1 Q H i k s P k l s + R k l s k = 1 Q H j k s P k l s + R k l s = k = 1 Q H i k s 0 + R k l s k = 1 Q H j k s 0 + R k l s = k = 1 Q H i k s R k l s k = 1 Q H j k s R k l s
Equation (21) shows that at other frequency locations, t i j l s is also random due to the randomness of R l s at different time periods.
The above derivation shows that the TRF has similar properties at the periodic excitation frequency γ r and at the system pole λ m . γ r is the real periodic pole of the equivalent system, so the false mode cannot be rejected by conventional methods and will interfere with the identification of the real modal parameters of the structure.

2.3. Delayed TRF Approach for Periodic False Mode Identification

The acceleration dynamic response of two DOFs in the numerical model of a multiple DOF system at time periods T1 and T2 is shown in Figure 2, and a set of delayed TRFs t i j 1 , 2 s and t i j 2 , 1 s can be constructed, defined as follows:
t i j 1 , 2 s = Y i 1 s Y j 2 s ,   t i j 2 , 1 s = Y i 2 s Y j 1 s
where Y 1 s and Y 2 s are the Laplace transformations of the acceleration response in time period 1 and time period 2, respectively; i and j are different response DOFs.
Due to the limited duration of time periods T1 and T2, the spectra of the periodic and random components of the system excitation during the two time periods satisfy Equations (16) and (17).
From Equations (4) and (15)–(19), when s λ m , the following exist:
lim s λ m t i j 1 , 2 s = lim s λ m k = 1 Q s λ m H i k s F k 1 s lim s λ m k = 1 Q s λ m H j k s F k 2 s = k = 1 Q ψ m , i l m , k P k 1 λ m + R k 1 λ m k = 1 Q ψ m , j l m , k P k 2 λ m + R k 2 λ m = k = 1 Q ψ m , i l m , k 0 + R k 1 λ m k = 1 Q ψ m , j l m , k 0 + R k 2 λ m = ψ m , i ψ m , j k = 1 Q l m , k R k 1 λ m k = 1 Q l m , k R k 2 λ m
Due to the randomness of R l s , generally R k 1 λ m R k 2 λ m , the delayed transmission rates at different combinations of time periods are not the same at the system’s intrinsic pole λ m , i.e.,
t i j 1 , 2 λ m t i j 2 , 3 λ m t i j k 1 , k λ m
Similarly, from Equations (4) and (15)–(19), the following equation is obtained when s γ r :
lim s γ r t i j 1 , 2 s = lim s γ r k = 1 Q H i k s F k 1 s lim s γ r k = 1 Q H j k s F k 2 s = k = 1 Q H i k γ r P k 1 γ r + R k 1 γ r k = 1 Q H j k γ r P k 2 γ r + R k 2 γ r = k = 1 Q H i k γ r a r k + o a r k k = 1 Q H j k γ r a r k + o a r k k = 1 Q a r k H i k γ r k = 1 Q a r k H j k γ r
Analogously, it can be obtained that
lim s γ r t i j 2 , 1 s k = 1 Q a r k H i k γ r k = 1 Q a r k H j k γ r
Comparing Equations (25) and (26), it is obtained that
t i j 1 , 2 γ r t i j 2 , 1 γ r
The value of Equation (26) is independent of the selected time period, so the delayed transmissibility at each combination of time periods is approximately equal at the periodic excitation frequency γ r , i.e.,
t i j 1 , 2 γ r t i j 2 , 3 γ r t i j k 1 , k γ r
It is clear that the delayed TRF t i j 1 , 2 s intersects at the periodic excitation frequency γ r and varies at the natural frequency λ m of the system. The corresponding proposed frequency response function is constructed according to Equation (6), and its peak indicates only periodic false modes. Thus, by comparing the modal parameter identification results of the conventional TRF and the delayed TRF, the effective identification of the real modal parameters of the structure under periodic excitation can be achieved.

3. Example Analysis

3.1. Numerical Calculation of a 4-DOF System

3.1.1. Dynamic Response of the System

The modal parameters were identified for a 4-DOF spring-damped-mass system as shown in Figure 3. Excitations F1 and F2 for 100 s were applied to degrees of freedom 1 and 2, and both loads contained both random and periodic excitations, as shown in Table 1. The parameters of the system were as follows: k1~k4 = 1000 N/m, m1~m4 = 1 kg, and c1~c4 = 0.5 N·s/m. The Newmark method was used to solve the acceleration dynamic response of the system with a sampling frequency of 200 Hz, where the acceleration responses of DOFs 3 and 4 in 40–50 s are shown in Figure 4. The overall calculation showed that the natural frequency from 1st to 4th order of the system was in the range of 0–10 Hz, so the sampling frequency satisfied Shannon’s sampling theorem ( f s 2 f max ).

3.1.2. Modal Parameter Identification Based on TRF

The acceleration dynamic response of 100 s was divided into 20 groups of equal length, the TRF t34 between DOF 3 and DOF 4 was estimated using the Hv method, and the obtained 20 groups of t34 are shown in Figure 5, which is represented by a curve with random colors. The red dashed line in Figure 5 represents the theoretical value of natural frequency in each order of the system, and the blue dashed line represents the periodic frequency in the excitation. The TRFs around the above frequencies are shown in Figure 6. It can be seen that the overall TRFs of the multiple dynamic responses had significant discretion due to the random components in the excitation, but always intersected at the natural and periodic excitation frequencies of the system.
The proposed frequency response function matrix was constructed according to Equation (6); Figure 7, which contains six peaks, shows element T H 1 , 1 ( s ) : the 4th-order intrinsic frequency peak and the 2nd-order periodic excitation peak.
The delayed TRF t 34 i j (i = 1, 2, …, 19; j = i + 1) of the acceleration response of DOF 3 and DOF 4 were constructed according to Equation (18), as shown in Figure 8 and Figure 9, and it can be seen that the delayed TRFs intersected only at the periodic excitation frequency.
At the same time, the quasi-frequency response function was constructed according to Equation (6), and the pole indication curve was obtained as shown in Figure 10. At this point, only the 2nd-order frequency peak of the periodic excitation is included in Figure 10.
The peak frequencies of the two curves in Figure 7 and Figure 10 were identified; comparison of the two could accurately distinguish the natural frequency of the system, and then the true modal vibration pattern and damping ratio of the structure could be identified. The theoretical and identified values of the system modal frequencies and damping ratios are shown in Table 2, and the theoretical and identified values of the modal vibration patterns are displayed in Figure 11. These identification results show that the TRF method has high accuracy for the identification of system modal parameters.

3.1.3. Robustness Check of the Delayed TRF

The above arithmetic example confirms that the joint delayed TRF approach could effectively identify the intrinsic modal parameters of the system. Thus, considering that in some extreme cases, the periodic excitation energy is lower or closer to some natural frequencies, the following calculation example is set out in this section to test the detection effect of the periodic excitation frequency in such cases. The four sets of periodic components in the load excitation included 1.9 Hz (fn1 + 0.15 Hz) and 7.91 Hz (fn3 + 0.2 Hz), close to the natural frequency, and 3.1 Hz (fn1 + 1.35 Hz) and 6 Hz (fn2 + 0.97 Hz), far from the natural frequency. On this basis, the signal-to-noise ratios of the periodic and random components of the excitation were set to 20 dB, 0 dB, and −20 dB, representing cases of strong, moderate, and weak periodic excitation, respectively.
The conventional TRF and delayed TRF were used to establish the proposed frequency response function, respectively (Figure 12). From Figure 12, even when the intensity of the periodic component of the excitation was relatively weak or its frequency was closer to the natural frequency of the structure, the frequency of the periodic excitation could be accurately identified using the proposed frequency response function of the delayed TRF, so the method has strong robustness for identifying periodic false modes.

3.2. Modal Parameter Identification Example of a Pump Station Plant

The considered medium-sized pumping station adopts the vertical unit type (Figure 13) and the rated speed of the pump unit is 500 r/min. The static and dynamic loads resulting from pump operation are mainly carried by the two lower main beams. To identify the dynamic characteristics of the main beam in the vertical direction, sensors were placed on its lower part as shown in Figure 14. The acceleration response of the structure was measured under stable operation of the unit, with a sampling frequency of 1024 Hz. The acceleration response of each measurement point in 10 s is presented in Figure 15.
The quasi-frequency response was constructed using the conventional transmissibility and delayed transmissibility, as shown in Figure 16.
Comparing the poles in Figure 16a,b yielded the natural and periodic frequencies of the structured system, as shown in Table 3. The unit rotational frequency was 8.33 Hz, which was converted from the rated speed of the pump, 500 r/min. The frequency of each periodic excitation in Table 3 is a multiple of the rotational frequency.
The FSDD method was further applied to identify the modal vibration pattern of the main beam in each order, as shown in Figure 17. The MAC indexes of the modal vibration identification results are shown in Table 4, which shows that the modal vibration modes were approximately linearly independent of each other, verifying the accuracy of the modal vibration modes. The corresponding damping ratios of each order were 1.6%, 2.2%, and 2.8%. The above results show that the combined transmissibility and delayed transmissibility methods effectively identified the periodic false modes of the main beam of the pump pier and accurately identified the modal parameters of the structure.

4. Conclusions

To address the interference problem of periodic components in ambient excitation when identifying structural modal parameters, this paper proposes a method for identifying structural modal parameters based on the delayed TRF of the dynamic response, and presents a verification of the effectiveness and accuracy of the proposed method through numerical models and practical engineering calculations. The main conclusions are as follows:
Through theoretical derivation, the property was demonstrated that the delayed TRF tends to intersect at the periodic frequencies of the excitation and varies at the natural frequencies of a structure, and thus, the TRF can be used to indicate the periodic frequencies in the excitation. The joint application of the delayed TRF and the regular TRF could be used to identify the periodic false modes of the structure and thus extract the real mode parameters of the structure. The numerical calculation of the 4-DOF system confirmed that the method described in this paper can effectively eliminate the interference of periodic false modes and identify the real modal parameters of a system, and the applicability of the method was tested in the case of weak periodic excitation or periodic excitation with a frequency close to the natural frequency of the structure. Furthermore, based on the measured dynamic response data of a pump station structure, the modal parameters of the main beam of the pump pier were effectively identified.
In this study, the excitation generated by rotating mechanical equipment was simplified to ideal periodic excitation, and a modal recognition method based on the transfer rate function is proposed, which provided better recognition of the working condition of the unit’s smooth operation. However, the frequency or amplitude of the periodic excitation is constantly changing in the case of nonsmooth operating conditions, such as frequent changes in the load or input power of the unit. The effect of the periodic excitation on the identification of structural modal parameters and the effective identification method need to be further studied.

Author Contributions

Conceptualization, Y.X., D.Z. and C.S.; methodology, Y.X.; software, Y.X. and C.S.; validation, D.Z., S.Z. and H.G.; formal analysis, Y.X.; investigation, S.Z. and H.G.; resources, Y.X., D.Z. and C.S.; data curation, Y.X. and C.S.; writing—original draft preparation, Y.X.; writing—review and editing, D.Z., C.S. and H.G.; visualization, C.S. and S.Z.; supervision, D.Z. and H.G.; funding acquisition, D.Z., C.S. and H.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (52209159, 52179128, U2243223), the Fundamental Research Funds for the Central Universities (B220201042), National Key R&D Program (2018YFC1508603), Water Conservancy Science and Technology Project of Jiangsu (2022024), China Postdoctoral Science Foundation (2021M701044), Jiangsu Young Science and Technological Talents Support Project TJ-2022-076.

Data Availability Statement

Not applicable.

Acknowledgments

The support of Hohai University, China is also gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, J.J.; Yun, C.B. Damage diagnosis of steel girder bridges using ambient vibration data. Eng. Struct. 2006, 28, 912–925. [Google Scholar] [CrossRef]
  2. Han, J.; Zheng, P.; Wang, H. Structural modal parameter identification and damage diagnosis based on Hilbert-Huang transform. Earthq. Eng. Eng. Vib. 2014, 13, 101–111. [Google Scholar] [CrossRef]
  3. Hakim, S.J.S.; Razak, H.A.; Ravanfar, S.A. Fault diagnosis on beam-like structures from modal parameters using artificial neural networks. Measurement 2015, 76, 45–61. [Google Scholar] [CrossRef]
  4. Foti, D.; Giannoccaro, N.I.; Vacca, V.; Lerna, M. Structural operativity evaluation of strategic buildings through finite element (FE) models validated by operational modal analysis (OMA). Sensors 2020, 20, 3252. [Google Scholar] [CrossRef]
  5. Tsuchimoto, K.; Narazaki, Y.; Spencer Jr, B.F. Development and validation of a post-earthquake safety assessment system for high-rise buildings using acceleration measurements. Mathematics 2021, 9, 1758. [Google Scholar] [CrossRef]
  6. Magalhães, F.; Cunha, Á.; Caetano, E. Vibration based structural health monitoring of an arch bridge: From automated OMA to damage detection. Mech. Syst. Signal Process. 2012, 28, 212–228. [Google Scholar] [CrossRef]
  7. Yan, Z.; Liu, H. SMoCo: A powerful and efficient method based on self-supervised learning for fault diagnosis of aero-engine bearing under limited data. Mathematics 2022, 10, 2796. [Google Scholar] [CrossRef]
  8. Pan, C.; Ye, X.; Mei, L. Improved automatic operational modal analysis method and application to large-scale bridges. J. Bridg. Eng. 2021, 26, 04021051. [Google Scholar] [CrossRef]
  9. Yan, W.J.; Zhao, M.J.; Sun, Q.; Ren, W.X. Transmissibility-based system identification for structural health monitoring: Fundamentals, approaches, and applications. Mech. Syst. Signal Process. 2019, 117, 453–482. [Google Scholar] [CrossRef]
  10. Sun, Q.; Yan, W.J.; Ren, W.X.; Liu, L.L. Application of transmissibility measurements to operational modal analysis of railway, highway, and pedestrian cable-stayed bridges. Measurement 2019, 148, 106880. [Google Scholar] [CrossRef]
  11. Liu, T.; Xu, H.; Ragulskis, M.; Cao, M.; Ostachowicz, W. A data-driven damage identification framework based on transmissibility function datasets and one-dimensional convolutional neural networks: Verification on a structural health monitoring benchmark structure. Sensors 2020, 20, 1059. [Google Scholar] [CrossRef] [Green Version]
  12. Devriendt, C.; Guillaume, P. The use of transmissibility measurements in output-only modal analysis. Mech. Syst. Signal Process. 2007, 21, 2689–2696. [Google Scholar] [CrossRef]
  13. Lage, Y.E.; Maia, N.M.M.; Neves, M.M.; Ribeiro, A.M.R. Force identification using the concept of displacement transmissibility. J. Sound Vib. 2013, 332, 1674–1686. [Google Scholar] [CrossRef]
  14. Araújo, I.G.; Laier, J.E. Operational modal analysis using SVD of power spectral density transmissibility matrices. Mech. Syst. Signal Process. 2014, 46, 129–145. [Google Scholar] [CrossRef]
  15. Araújo, I.D.G. Transmissibility-based operational modal analysis: Unified concept and its application. Mech. Syst. Signal Process. 2022, 178, 109302. [Google Scholar] [CrossRef]
  16. Devriendt, C.; Guillaume, P. Identification of modal parameters from transmissibility measurements. J. Sound Vib. 2008, 314, 343–356. [Google Scholar] [CrossRef]
  17. Devriendt, C.; De Sitter, G.; Guillaume, P. An operational modal analysis approach based on parametrically identified multivariable transmissibilities. Mech. Syst. Signal Process. 2019, 24, 1250–1259. [Google Scholar] [CrossRef]
  18. Devriendt, C.; Weijtjens, W.; De Sitter, G.; Guillaume, P. Combining multiple single-reference transmissibility functions in a unique matrix formulation for operational modal analysis. Mech. Syst. Signal Process. 2013, 40, 278–287. [Google Scholar] [CrossRef]
  19. Weijtjens, W.; De Sitter, G.; Devriendt, C.; Guillaume, P. Operational modal parameter estimation of MIMO systems using transmissibility functions. Automatica 2014, 50, 559–564. [Google Scholar] [CrossRef]
  20. Zhang, Y.N.; Wang, T.; Xia, Z.P. Transmissibility based operational modal analysis. J. Vib. Meas. Diagn. 2015, 35, 945–949+995, (In Chinese with English abstract). [Google Scholar]
  21. Araújo, I.G.; Sánchez, J.A.G.; Andersen, P. Modal parameter identification based on combining transmissibility functions and blind source separation techniques. Mech. Syst. Signal Process. 2018, 105, 276–293. [Google Scholar] [CrossRef]
  22. Li, X.Z.; Dong, X.J.; Yue, X.B.; Huang, W.; Peng, Z.K. Dynamic characteristics of vibration response transmissibility and its application in operational modal analysis. J. Vib. Shock. 2019, 38, 62–70, (In Chinese with English abstract). [Google Scholar]
  23. Mohanty, P.; Rixen, D.J. Operational modal analysis in the presence of harmonic excitation. J. Sound Vib. 2004, 270, 93–109. [Google Scholar] [CrossRef]
  24. Mohanty, P.; Rixen, D.J. Modified ERA method for operational modal analysis in the presence of harmonic excitations. Mech. Syst. Sign. Process. 2006, 20, 114–130. [Google Scholar] [CrossRef]
  25. He, X.H.; Hua, X.G.; Chen, Z.Q.; Huang, F.L. EMD-based random decrement technique for modal parameter identification of an existing railway bridge. Eng. Struct. 2011, 33, 1348–1356. [Google Scholar] [CrossRef]
  26. Carden, E.P.; Brownjohn, J.M. Fuzzy clustering of stability diagrams for vibration-based structural health monitoring. Comput. Aided Civ. Infrastruct. Eng. 2008, 23, 360–372. [Google Scholar] [CrossRef]
  27. Lian, J.; Li, H.; Zhang, J. ERA modal identification method for hydraulic structures based on order determination and noise reduction of singular entropy. Sci. China Ser. E Technol. Sci. 2009, 52, 400–412. [Google Scholar] [CrossRef]
  28. Dong, X.; Lian, J.; Yang, M.; Wang, H. Operational modal identification of offshore wind turbine structure based on modified stochastic subspace identification method considering harmonic interference. J. Renew. Sustain. Energy 2014, 6, 033128. [Google Scholar] [CrossRef]
  29. Heylen, W.; Lammens, S.; Sas, P. Modal Analysis Theory and Testing; Katholieke Universiteit Leuven: Leuven, Belgium, 1997; Volume 200. [Google Scholar]
  30. Zhang, L.M.; Wang, T.; Tamura, Y. A frequency-spatial domain decomposition (FSDD) method for operational modal analysis. Mech. Syst. Signal Process. 2010, 24, 1227–1239. [Google Scholar] [CrossRef]
  31. Pintelon, R.; Schoukens, J. System Identification: A Frequency Domain Approach; John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar]
Figure 1. Input–output relation diagram of two types of structure: (a) random excitation; (b) contains periodic excitation.
Figure 1. Input–output relation diagram of two types of structure: (a) random excitation; (b) contains periodic excitation.
Mathematics 11 01019 g001
Figure 2. Acceleration responses of measurement points in different time periods.
Figure 2. Acceleration responses of measurement points in different time periods.
Mathematics 11 01019 g002
Figure 3. Four-degree-of-freedom spring-damping-mass system.
Figure 3. Four-degree-of-freedom spring-damping-mass system.
Mathematics 11 01019 g003
Figure 4. Acceleration dynamic response of the system. (a) DOF 3; (b) DOF 4.
Figure 4. Acceleration dynamic response of the system. (a) DOF 3; (b) DOF 4.
Mathematics 11 01019 g004
Figure 5. Transmissibility t34 of response in a finite time period (20 groups).
Figure 5. Transmissibility t34 of response in a finite time period (20 groups).
Mathematics 11 01019 g005
Figure 6. TRFs near each characteristic frequency. (a) Natural frequency 1; (b) natural frequency 2; (c) natural frequency 3; (d) natural frequency 4; (e) excitation frequency 1; (f) excitation frequency 2.
Figure 6. TRFs near each characteristic frequency. (a) Natural frequency 1; (b) natural frequency 2; (c) natural frequency 3; (d) natural frequency 4; (e) excitation frequency 1; (f) excitation frequency 2.
Mathematics 11 01019 g006
Figure 7. Natural frequency indication curve of T H 1 , 1 ( s ) .
Figure 7. Natural frequency indication curve of T H 1 , 1 ( s ) .
Mathematics 11 01019 g007
Figure 8. Delayed TRF t 34 i j of response in a finite time period (19 groups).
Figure 8. Delayed TRF t 34 i j of response in a finite time period (19 groups).
Mathematics 11 01019 g008
Figure 9. Local delayed TRF. (a) Excitation frequency 1; (b) excitation frequency 2.
Figure 9. Local delayed TRF. (a) Excitation frequency 1; (b) excitation frequency 2.
Mathematics 11 01019 g009
Figure 10. Natural frequency indication curve based on delayed transmissibility.
Figure 10. Natural frequency indication curve based on delayed transmissibility.
Mathematics 11 01019 g010
Figure 11. Identification results of the modal vibration pattern. (a) 1st order; (b) 2nd order; (c) 3rd order; (d) 4th order.
Figure 11. Identification results of the modal vibration pattern. (a) 1st order; (b) 2nd order; (c) 3rd order; (d) 4th order.
Mathematics 11 01019 g011
Figure 12. The virtual frequency response function with different period excitation amplitudes. (a) Conventional transmissibility (SNR = 20 dB); (b) delayed transmissibility (SNR = 20 dB); (c) conventional transmissibility (SNR = 0 dB); (d) delayed transmissibility (SNR = 0 dB); (e) conventional transmissibility (SNR = −20 dB); (f) delayed transmissibility (SNR = −20 dB).
Figure 12. The virtual frequency response function with different period excitation amplitudes. (a) Conventional transmissibility (SNR = 20 dB); (b) delayed transmissibility (SNR = 20 dB); (c) conventional transmissibility (SNR = 0 dB); (d) delayed transmissibility (SNR = 0 dB); (e) conventional transmissibility (SNR = −20 dB); (f) delayed transmissibility (SNR = −20 dB).
Mathematics 11 01019 g012
Figure 13. Schematic diagram of the pumping unit of a medium-sized pumping station.
Figure 13. Schematic diagram of the pumping unit of a medium-sized pumping station.
Mathematics 11 01019 g013
Figure 14. Arrangement of measuring points of the main beam of the pier base. (a) Diagram of measuring point arrangement; (b) upstream-side main beam acceleration measurement point.
Figure 14. Arrangement of measuring points of the main beam of the pier base. (a) Diagram of measuring point arrangement; (b) upstream-side main beam acceleration measurement point.
Mathematics 11 01019 g014
Figure 15. Acceleration response of each measurement point. (a) Measurement point 1; (b) measurement point 2; (c) measurement point 3; (d) measurement point 4; (e) measurement point 5.
Figure 15. Acceleration response of each measurement point. (a) Measurement point 1; (b) measurement point 2; (c) measurement point 3; (d) measurement point 4; (e) measurement point 5.
Mathematics 11 01019 g015
Figure 16. System pole recognition based on the virtual frequency response function. (a) Conventional TRF; (b) delayed TRF.
Figure 16. System pole recognition based on the virtual frequency response function. (a) Conventional TRF; (b) delayed TRF.
Mathematics 11 01019 g016
Figure 17. Identification results of vibration modes of the main beam. (a) 1st order; (b) 2nd order; (c) 3rd order.
Figure 17. Identification results of vibration modes of the main beam. (a) 1st order; (b) 2nd order; (c) 3rd order.
Mathematics 11 01019 g017
Table 1. Two dynamic load components.
Table 1. Two dynamic load components.
LoadF1F2
Component 1Periodic excitation 1 (3.1 Hz)Periodic excitation 2 (6 Hz)
Component 2White noise 1White noise 1
1 The white noise components of F1 and F2 are uncorrelated with each other, and the signal-to-noise ratio of the periodic component relative to the white noise is 10 dB.
Table 2. Identification results of the modal frequency and damping ratio of the 4-DOF system.
Table 2. Identification results of the modal frequency and damping ratio of the 4-DOF system.
Modal
Number
Frequency (Hz)Damping Ratio (%)
Theoretical ValueIdentified ValueError Rate (%)Theoretical ValueIdentified ValueError Rate (%)
11.7481.7500.100.270.283.70
25.0335.030−0.060.790.77−2.53
37.7117.693−0.231.211.242.48
49.4599.450−0.091.491.522.01
Table 3. Frequency identification results of the main beam.
Table 3. Frequency identification results of the main beam.
Frequency TypeFrequency Value (Hz)
Natural frequency34, 43.5, 61.5
Periodic frequency8.33, 16.67, 25, 41.7, 50, 58.3, 66.7, 83.3
Table 4. MAC index of modal shape recognition results.
Table 4. MAC index of modal shape recognition results.
Order of Modes1st Order2nd Order3rd Order
1st order1.0000.0050.001
2nd order0.0051.0000.002
3rd order0.0010.0021.000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, Y.; Zheng, D.; Shao, C.; Zheng, S.; Gu, H. Structural Modal Parameter Identification Method Based on the Delayed Transfer Rate Function under Periodic Excitations. Mathematics 2023, 11, 1019. https://doi.org/10.3390/math11041019

AMA Style

Xu Y, Zheng D, Shao C, Zheng S, Gu H. Structural Modal Parameter Identification Method Based on the Delayed Transfer Rate Function under Periodic Excitations. Mathematics. 2023; 11(4):1019. https://doi.org/10.3390/math11041019

Chicago/Turabian Style

Xu, Yanxin, Dongjian Zheng, Chenfei Shao, Sen Zheng, and Hao Gu. 2023. "Structural Modal Parameter Identification Method Based on the Delayed Transfer Rate Function under Periodic Excitations" Mathematics 11, no. 4: 1019. https://doi.org/10.3390/math11041019

APA Style

Xu, Y., Zheng, D., Shao, C., Zheng, S., & Gu, H. (2023). Structural Modal Parameter Identification Method Based on the Delayed Transfer Rate Function under Periodic Excitations. Mathematics, 11(4), 1019. https://doi.org/10.3390/math11041019

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop