A new instantaneous frequency extraction method for nonstationary response signals in civil engineering structures

Transcription

1 Original Article A new instantaneous frequency extraction method for nonstationary response signals in civil engineering structures Journal of Low Frequency Noise, Vibration and Active Control 2018, Vol. 37(4) ! The Author(s) 2018 DOI: / journals.sagepub.com/home/lfn Jing-Liang Liu 1,2, Jin-Yang Zheng 1, Xiao-Jun Wei 2, Fei-Yu Liao 1 and Yong-Peng Luo 1 Abstract Civil engineering structures in essence belong to time-varying and nonlinear structures and the resultant dynamic responses are usually nonstationary signals. A new method is proposed to extract instantaneous frequency of such nonstationary signals. The new method combines the extended analytical mode decomposition, normalization scheme, and improved arccosine function algorithm. In this method, a multicomponent response signal is first decomposed into several mono-components signals by the extended analytical mode decomposition theorem, and then the extracted mono-components are demodulated into purely frequency modulated signals by a repeated normalization scheme. After that, the improved arccosine function is presented to compute the instantaneous phase and instantaneous frequency. A numerical example and an experimental study of a cable with time-varying tension force are provided to illustrate the working of the method. The results demonstrate that the proposed method can extract instantaneous frequencies from multicomponent signals more accurately than Hilbert Huang transform. Keywords Instantaneous frequency, civil engineering structures, extended analytical mode decomposition, normalization, improved arccosine function Introduction Practical civil engineering structures essentially belong to time-varying and nonlinear structures and the resultant dynamic responses are often nonstationary due to the variation of dynamic properties under service loads or extreme loads. However, traditional signal processing techniques, including time- and frequency-domain analyses, are based on the assumption that the response signals are stationary and linear. 1 For example, the Fourier transform can only give meaningful interpretation to linear and stationary signals, but its application to nonstationary signals of time-varying or nonlinear structures is often problematical. 2 In contrast, time frequency analysis (TFA) method is potential to process nonstationary signals and characterize their time-varying features. Various TFA methods, including short-time Fourier transform, 3 quadratic distributions, 4 empirical mode decomposition (EMD)-based Hilbert Huang transform (HHT), 5,6 ensemble empirical modal decomposition (EEMD), 7 and continuous wavelet transform, 8,9 have been proposed in recent years to deal with nonstationary signals. However, for some civil engineering structures, the response signals may be composed of closely spaced or 1 College of Transportation and Civil Engineering, Fujian Agriculture and Forestry University, Fuzhou, China 2 School of Engineering, University of Warwick, Coventry, UK Corresponding author: Xiao-Jun Wei, School of Engineering, University of Warwick, Coventry CV4 7AL, UK. X.Wei.3@warwick.ac.uk Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License ( creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (

2 Liu et al. 835 low-frequency components. In such cases, some classic TFA methods such as Hilbert transform (HT) cannot work well. To address this issue, it is essential to introduce effective methods for signal decomposition. 10 EMD is designed to extract intrinsic mode functions from response signals adaptively by cubic spline interpolation. However, EMD-based methods are empirical methods and still face several challenges in some engineering applications. For example, EMD is difficult to decompose closely spaced frequency components. 11 By injecting a generated white noise into the original signal, EEMD may be able to separate closely spaced components but the corresponding calculation amount increases dramatically. In addition, the signal decomposition results using EEMD are sensitive to parameter selection, even pseudo-components may be produced if improper parameters are selected. 12 The wavelet transform as an advanced TFA technique is suitable for nonstationary signal processing and it has received wide acceptance in the dynamic signal analysis of time-varying structures in recent years. 13,14 However, Kijewski and Kareem 15 argued that civil engineering structures usually possess low-frequency motions and thus require finer frequency resolution. Moreover, parameter selection, which remained not well resolved, is critical for wavelet transform and it may significantly affect instantaneous frequency (IF) extraction results. Recently, Daubechies et al. 16 proposed synchrosqueezing wavelet transform (SWT) to reshape time frequency curves and reconstruct individual mono-components. Modified algorithms and successful applications about SWT can be found in Thakur et al., 17 Li and Liang, 18 and Cao et al. 19 It should be noted here that there is an assumption to hold by default for the SWT, that is the components of response signals should be asymptotic and well separated in frequency domain. However, some civil engineering structures such as large span cable-stayed bridges with dense stayed cables may break the assumption of well-separated frequencies and asymptotic signals. To address this issue, a new signal decomposition theorem called analytical mode decomposition (AMD) 20 was presented and on this basis, the author of this paper and his cooperators extends AMD theorem to solve the problem of modal mixing. 21 Numerical simulation results indicate that each frequency-modulated (FM) individual component between any two bisecting frequencies can be analytically extracted by extended AMD. However, the mono-components decomposed by the extended AMD are generally amplitude-modulated (AM) and frequency-modulated (FM) functions, and the amplitude change rate of the mono-component signal is usually not much slower than its phase change rate. Therefore, the significance of the demodulation process is highlighted. The local mean decomposition (LMD) method developed by Smith 22 is a self-adaptive algorithm to decompose AM FM signals into a small set of product functions, each of which is the product of an envelope signal and a FM signal. The LMD avoids the limitation of negative frequencies occurring in HT, but end effect and mode mixing problems remain two main limitations in LMD, which can be alleviated by proper parameter selection with regard to boundary condition, envelope estimation, and sifting stopping criterion. 23 Huang et al. 10 proposed a normalization scheme to demodulate signals through empirical amplitude method. The normalization scheme is exactly the process of demodulation based on iterative applications of cubic spline fitting and thus the AM and FM parts of the decomposed components are separated uniquely. With this repeated demodulation process, we can conduct HT or other TFA methods on the FM part alone, avoiding the difficulty stated in the Bedrosian theorem and the Nuttall theorems. 2 The normalization scheme not only successfully extract FM parts of original signals but also converge very fast in the repeated demodulation process. 10 Nevertheless, the normalization process could cause some deformation of the original signal, but the amount of the deformation is negligible, for there are rigid controlling points for the periodicity provided by the zero-crossing points in addition to the extrema. The next step is to extract IFs from FM signals. The most common methods for calculating instantaneous phase or IF include HT, Teager energy operator (TEO), zero-crossing, and arccosine function method. 2,10 Among them, HT is the most popular and direct method to extract IF. However, the IF extracted by HT may be negative at times, which is extremely hard to interpret. In addition, HT is very sensitive to random noises and hence leads to an IF curve with blurred and distorted lines. The TEO is a method to compute IF totally based on differentiations. A distinct advantage of the TEO is its excellent localization property and no integral transform needed as in HT. The shortcoming of this method is that TEO is based on a linear model for a single harmonic component only. Therefore, the approximation produced by the TEO will deteriorate and even break down when either the amplitude is a function of time or the waveforms have any interweave modulations or harmonic distortions. 10 The zero-crossing method has long been used to compute the mean frequency for narrow band signals. But unfortunately, the calculation results are relatively crude for the frequency should be constant over the period between two zero-crossings. The arccosine function is potential to calculate the phase from FM signals, and then the IF can be expressed as the derivative of the instantaneous phase with respect to time. Although the arccosine function is able to calculate the phase successfully, its range is strictly restricted to ½0; pš. Therefore, the

3 836 Journal of Low Frequency Noise, Vibration and Active Control 37(4) arccosine function cannot be directly applied to the purely FM signal whose phase range actually is the set of real numbers. 24 To address this issue, a new method called improved arccosine function is proposed in this paper to solve the instantaneous phase and IF more accurately based on the division of intervals along time axis. The main limitation of improved arccosine function method is the distortion of the instantaneous phase, which can be greatly reduced by a moving average algorithm. 23 With the moving average method, the instantaneous phase curves beforehand extracted by the improved arccosine function are smoothed without great distortions. The contribution of this paper is to propose a new IF extraction method, which is different from the traditional HHT. The proposed new method, which combines the extended AMD, normalization scheme, and improved arccosine function, is composed of three steps accordingly. First, the extended AMD is employed to decompose a response signal, which may include closely spaced or low-frequency components, into several mono-components. Second, the normalization scheme is used to demodulate the decomposed mono-component signal and thus the AM and FM parts are completely separated. Finally, the improved arccosine function is presented to compute the instantaneous phase. To reduce the distortion of the instantaneous phase caused by segmentation error, a moving average algorithm is introduced to smooth the target instantaneous phase curves. Then, the IF can be solved by taking a derivative of instantaneous phase with respect to time. The main difference between the proposed method and HHT is that the decomposition is not empirical but analytical. In addition, a major highlight of the new method is the presentation of the improved arccosine function. Theory and methodology The method proposed in this paper combines the techniques of extended AMD, normalization scheme, and improved arccosine function. It allows to decompose and demodulate multicomponent signals composed of closely spaced frequency or low-frequency components into mono-components with separated AM and FM parts, whereby the improved arccosine function can be used to achieve good time frequency resolutions. The flowchart of this new method is illustrated in Figure 1. The extended AMD theorem is first introduced to decompose a multicomponent response signal into several mono-components by selecting proper bisecting frequencies. A normalization scheme is then employed to demodulate the decomposed mono-component signals into purely FM. After that, the improved arccosine function is presented to characterize the IFs of the purely FM signals with a better time frequency resolution. Original signal If the signal is a multicomponent signal No Yes Step 1: Decompose the signal by the extended AMD theorem Step 2: Demodulate the extracted components by a normalization scheme Step 3: IF extraction by an improved arccosine function Stop Figure 1. The flowchart of the new method. AMD: analytical mode decomposition; IF: instantaneous frequency.

4 Liu et al. 837 Signal decomposition by the extended AMD theorem For nonstationary AM FM signals, the frequency of each individual component varies with time. Instead of selecting constant bisecting frequencies, time-varying bisecting frequencies estimated from a preliminary SWT are selected, and each FM component between any two bisecting frequencies can then be analytically extracted by the extended AMD theorem. Let x(t) denote a real response signal of r individual components with frequencies: x 1 ðtþ; x 2 ðtþ;...; x r ðtþ, all positive, in L 2 ( 1, þ1) of the real time variable t. It can be decomposed into r signals x ðdþ whose frequency ranges satisfy jx 1 ðþj t < x b1 ðtþ, x br 1 ð Þ ðtþ, and x br 1 ð Þ ðtþ < jx r ðtþj. That is xt ðþ¼ Xr x b1 ðtþ < jx 2 ðþj t < x b2 ðtþ;..., i¼1 x br 2 ð i ðtþ (i ¼ 1, 2,..., r) Þ ðtþ < jx r 1 ðþj t < x ðdþ i ðtþ (1) Here, x i ðtþ is a frequency variable for the decomposed signal x ðdþ i ðtþ, and x bi ðtþ 2ðx i ðtþ; x iþ1 ðtþþði ¼ 1; 2;...; r 1Þ is time-varying bisecting frequencies. Then, each individual signal can be determined by and x ðdþ 1 ¼ s 1 ðþ; t...; x ðdþ i ðþ¼ t s i ðþ t s i 1 ðþ; t...; xr d ðþ ðtþ ¼xt ðþ s r 1 ðtþ (2) s i ðþ¼ t sin½h bi ðtþšhfxðþcos½h t bi ðþšg t cos½h bi ðþšhfx t ðþsin t ½h bi ðþ t Šg ði ¼ 1; 2;...; r 1Þ (3) where Hfg returns the HT the function inside the square bracket and h bi ðþ¼ t R t 1 x biðsþds is the phase function of the ith bisecting frequency. x i ðtþ and x bi ðtþ can be estimated from a previous SWT T x ðx l ; bþ. 16,17 More details about the proofs and numerical studies of the extended AMD theorem can be found in Wang et al. 21 Demodulation by a normalization scheme For a mono-component signal x 1 ðþdecomposed t by the extended AMD theorem, the analytic signal of x 1 ðþcan t be defined as a complex-valued function which has the original signal as its real part and the HT as its imaginary part zt ðþ¼x 1 ðþþ t i H½x 1 ðþ t Š ¼ A 1 ðþe t i/ 1 t (4) p where H½Šrepresents the HT of a signal in the square bracket, i ¼ ffiffiffiffiffiffi qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1, A1 ðþ¼ t ½x 1 ðþ t Š 2 2 þ H½x 1 ðþ t Š and / 1 ðþ¼ t arctan H½x 1 ðþ t Š=x 1 ðþ t are instantaneous amplitude and instantaneous phase, respectively. Then, x1 ðþ t can be expressed in equation (5) through the definition of the analytic signal ðþ x 1 ðþ¼ t A 1 ðþcos/ t 1 ðþ t (5) Supposing the Bedrosian theorem is satisfied, the HT for the right side of equation (5) can be rewritten as H A 1 ðþcos/ t 1 ðþ t ¼ A1 ðþh t cos/ 1 ðþ t (6) Actually equation (6) can hold if the amplitude varies so slowly that the frequency spectra of the envelope and the carrier waves are disjoint according to the Bedrosian theorem. 2 In other words, the signal to be analyzed should be asymptotic or purely FM, but this assumption is usually broken for response signals in practical civil engineering structures. To address this issue, the signal normalization scheme is used to demodulate the extracted component signals into AM and FM parts uniquely. First, we calculate absolute values of the decomposed monocomponent signal and then identify all extrema of these absolute values, followed by defining the envelope

5 838 Journal of Low Frequency Noise, Vibration and Active Control 37(4) Figure 2. The multicomponent signal xt ðþ¼ ( cos2pt0 t 5 cos4pt 5< t 10 with its zero-crossing points and extrema. with a spline through all these extrema and designating it as e 1 ðþ. t As such, the normalization d 1 ðþcan t be expressed as d 1 ðþ¼ t x 1 ðþ=e t 1 ðþ t (7) However, e 1 ðþcould t have values smaller than the decomposed signal x 1 ðþat t some temporal locations and cause the normalized function to have amplitude greater than unity when the fitted envelope e 1 ðþbased t on the extrema has cuts through x 1 ðþ. t These conditions may be rare but are not impossible. Therefore, the normalization should be further conducted recursively as equation (8) d 2 ðþ¼ t d 1 ðþ=e t 2 ðþ; t...; d m ðþ¼ t d m 1 ðþ=e t m ðþ t (8) This iteration process does not stop until the normalized function values are all unity. Set m as the number of total iterations, and the final demodulated FM part D 1 ðþand t AM part A 1 ðþof t x 1 ðþare t given by ( D 1 ðþ¼ t d m ðþ t A 1 ðþ¼ t x 1 ðþ=d t 1 ðþ t So far the demodulation process is successfully realized and the FM signal separated from mono-component signal x 1 ðþcan t be used in the next step. IF extraction by improved arccosine function One way to express the nonstationarity is to find instantaneous phase and IF using HT. With HT, the instantaneous phase is solved by the definition of analytical signal expressed in equation (4). Then, the IF can be calculated by taking the derivative of instantaneous phase with respect to time. However, a purely oscillatory function with a zero reference level is a necessary condition for the abovementioned HT method. Otherwise, the distortion of phase function will happen and thus cause the deviation of the extracted IF from its true value. 5 Moreover, the sensitivity to noises limits its practical application. An alternative way to compute the phase from the FM signal is the arccosine function. The arccosine function can calculate the phase successfully, but its range is strictly restricted to ½0; pš. Therefore, the arccosine function cannot be applied to the purely FM signal whose phase range actually is a set of real numbers. As such, we present a new method called improved arccosine function based on the division of intervals along time axis to solve the instantaneous phase more accurately than the standard arccosine function. Typically the time domain graph of a purely FM signal x 1 ðþis t composed of several waves as shown in Figure 2, with each having a wave crest and trough. We can find two zero-crossing points in each wave by using the operator sgnð Þ. Set n as the total time points of the purely FM signal x 1 ðþ, t and define a new time series st ð i Þ in equation (10) st ð i Þ ¼ x 1 ðt i 1 Þ:x 1 ðt i Þ ði ¼ 2;...; n (9) (10)

6 Liu et al. 839 Then, the operator sgnð Þ is expressed as equation (11) 8 >< sgnðst ð i ÞÞ ¼ >: 1; 1; st ð i Þ > 0 0; st ð i Þ ¼ 0 st ð i Þ < 0 (11) If sgnðst ð i ÞÞ ¼ 1, it implies that there is at least one time instance between two points t i and t i 1 at which the value of x 1 ðþequals t zero. Then the time point t i can be approximated as zero-crossing point with i as the index, provided that the sampling frequency for the signal x 1 ðþis t large enough. In other words, the point t i can be directly defined as the zero-crossing point when sgnðst ð i ÞÞ ¼ 0. In this way, we can find all p zero-crossing points and define the indexes of them as z 1 ; z 2 ;...; z p. Then the time domain graph of a purely FM signal x 1 ðþcan t be divided into p þ 1 intervals named as [0; t z1 ], [t z1 ; t z2 ],...,[t zp ; t n ], and each individual maximum/minimum extremum in the intervals can be extracted, with the indexes of the extrema are named as c 1, c 2,..., c p, c pþ1. By using these extracted extrema, the time axis is separated into p þ 1 intervals and defined as [0; t c1 ], [t c1 ; t c2 ],..., [t cp 1 ; t cp ], [t cp ; t cpþ1 ]. ( Taking the signal x 1 ðþ¼ t cos2pt 0 t 5 for example, we can find two extrema points and one zerocrossing point at any time interval [t ci ; t ciþ1 ] after assigning the operator sgnð Þ; which is shown in Figure 2. cos4pt 5 < t 10 Assuming x 1 ðþis t monotonic in each individual interval [t ci ; t ciþ1 ], we can draw a conclusion that the phases of x 1 ðþin t the first interval fall in the range of ½0; pš. Likewise, the phases of x 1 ðþin t the second interval are located in ½p; 2pŠ. By analogy, the phases of x 1 ðþin t the ith interval will definitely be ½ði 1Þp; ipš. Therefore, the instantaneous phase / 1 ðþof t x 1 ðþcan t be accurately computed as follows 8 >< / 1 ðþ¼ t >: arccos½x 1 ðþ t Š t 2½0; t z1 Š ð 1Þ i 1 i arccos½x 1 ðþ t Šþ 2p floor 2 t 2½t zi ; t ziþ1 Š where floorðþ represents the element inside the parenthesis to the nearest integer toward minus infinity. However, the distortions of the phases usually occur at the locations near extrema points. The reason for this phenomenon is that the extrema points defined in the improved arccosine function algorithm can only be infinitely close to the true extrema points due to the discretization of signals and thus the occurrence of distortions near the center of the extrema is not a surprise. In addition, noises have an impact on the calculation of instantaneous phase, but the effect is small because most high-frequency noises are already removed by the extended AMD filter. In order to solve the problem of phase distortions, it is necessary to find a useful mathematical tool to smooth the whole instantaneous phase curve. Moving average method is one of useful methods for smoothing the computed instantaneous phases. A simple moving average method is actually a common average of n consecutive data points in time series with each point equally weighted. In this paper, only five-point moving average method is used and the formula is defined as follows (12) f a ðt i Þ ¼ ft ð i 2Þþ ft ð i 1 Þþ ft ð i Þþ ft ð iþ1 Þþ ft ð iþ2 Þ 5 (13) where ft ð i Þ stands for the value of the frequency at t i and f a ðt i Þ represents the averaged or smoothed frequency at t i. It is noted that the smoothed frequencies at the first two and the last two points cannot be directly calculated by equation (13). One possibility is to update the old value at abovementioned locations with the counterparts at neighboring locations. By using the moving average method, the computation stability and identification accuracy of instantaneous phases along the whole time axis can be greatly enhanced. Thus, the IF can be further computed by ft ðþ¼ xðþ t 2p ¼ 1 d/ 1 ðþ t 2p dt (14)

7 840 Journal of Low Frequency Noise, Vibration and Active Control 37(4) in which xðþand t / 1 ðþrepresent t the circular frequency and the instantaneous phase of the decomposed component signal, respectively. Since the variables ft ðþand / 1 ðþare t often discretized in practical engineering, the discrete version of equation (14) is expressed as Numerical examples ft ð i Þ ¼ 1 / 1 ðt i Þ / 1 ðt i 1 Þ (15) 2p t i t i 1 In this section, a numerical simulation case is employed to verify the effectiveness and accuracy of the proposed method. In addition, an index of accuracy (IA) is presented to quantify the accuracy of the IF identification. A multicomponent signal with two closely spaced AM FM components is considered yt ðþ¼y 1 ðþþ t y 2 ðtþ (16) where y 1 ðþ¼½3 t þ cosðptþšcos½8pt þ 2sinðptÞŠ and y 2 ðþ¼½4 t þ cosðptþšcos½13pt þ sinð2ptþš. The theoretical frequencies of the two components are closely spaced and they are expressed as f 1 ¼ 4 þ cosðptþhz and f 2 ¼ 6:5 þ cosðptþhz, respectively. The sampling frequency is 100 Hz and the sampling time is 10 s. A simulated 20% Gaussian white noise is added to consider the influence of noises. The noise intensity is defined by the following signal-to-noise ratio (SNR) SNR ¼ 10log 10 A 2 signal A 2 noise ðin dbþ (17) where A signal and A noise refer to the root mean square value of the signal and the noise, respectively, and the noise level means the ratio of A 2 noise to A2 signal. The waveform of the noisy multicomponent signal to be analyzed is plotted in Figure 3. The SWT with Morlet transform as parent function is first performed on the noisy signal to obtain the wavelet scalogram of yðtþ as presented in Figure 4. By selecting 5 Hz as the bisecting frequency, two individual components can be extracted by the extended AMD, which are shown in Figure 5. Then the normalization scheme is used to separate the Figure 3. The simulated AM FM multicomponent signal. Figure 4. Wavelet scalogram of the AM FM multicomponent signal.

8 Liu et al. 841 Figure 5. Extracted mono-components: (a) y 1 and (b) y 2. AMD: analytical mode decomposition; EMD: empirical mode decomposition. Figure 6. IFs extracted by the new method and HHT: (a) y 1 and (b) y 2. HHT: Hilbert Huang transform. decomposed mono-components into AM and FM parts until the maximum is close to 1. The demodulation results are also illustrated in Figure 5. One can see from Figure 5 that the decomposed mono-components extracted by the extended AMD agree with theoretical results well while the counterparts extracted by EMD do not. Moreover, the periodic change of waveform of the purely FM signals obtained by the normalization scheme is also in good accordance with that of theoretical counterparts. Finally, the instantaneous phases and IFs are computed by the improved arccosine function and then smoothed by the moving average method. Figure 6 presents the results of IF extraction using the new method and HHT, which indicates that the new method can provide a more accurate IF curve than HHT. To quantify the accuracy of the IF identification, an IA is defined as the root-mean-squared value over the total time duration of the IF and expressed as equation (18) vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi Z T ½f d ðtþ f e ðtþš 2 0 IA ¼ u Z t T ½f e ðtþš 2 0 (18)

9 842 Journal of Low Frequency Noise, Vibration and Active Control 37(4) Table 1. IA of IF extraction of the noisy signal with AM FM components (%). Methods New method HHT IA IA AM: amplitude modulated; FM: frequency modulated; HHT: Hilbert Huang transform; IA: index of accuracy; IF: instantaneous frequency. Figure 7. The cable test setup. where f d ðtþ is the extracted IF, and f e ðtþ is the exact/theoretical IF. IA 1 and IA 2 represent the IF identification accuracies of components y 1 ðtþ and y 2 ðtþ, respectively, The smaller the value of IA, the better accuracy of IF extraction. According to equation (18), the values of IA using HHT and the new method are computed and then presented in Table 1. We can see from Table 1 that the value of IA using the new method is much smaller than that of the HHT, which verifies the superiority of the new method, especially when response signals are contaminated by random noises. Experimental case study IF identification of a cable with linearly varying tension forces To further validate the accuracy of the proposed method, a cable composed of pieces of 7U5 steel wires is considered. The cable has elastic modulus E ¼ MPa, area of cross section A ¼ m 2, and density of unit length q ¼ 1.1 kg/m. The stiffness of the cable is changed due to the applied time-varying tension force so that the natural frequency of the cable is time dependent. The cable is fixed at one end, and the other end is connected to an MTS loading system. The total length of cable between two ends is 4.55 m. A accelerometer is installed vertically at the midpoint of the cable. The test setup is shown in Figure 7. During the cable test, an initial constant pretension force was first applied by MTS load actuator to the cable. When the test setup and data collection were ready, the cable tension force was changed continuously using the

10 Liu et al. 843 MTS load system. At the same time, the impact hammer was used to generate free vibration and the vertical acceleration responses were recorded at a sampling frequency of 600 Hz. The tension force with linear or sinusoidal change was considered during the test. The theoretical IFs of the cable were obtained by solving eigenvalues and eigenvectors of vibration equations, assuming that the parameters of the cable keep invariant over a relatively short time interval, which is designated as the time frozen method. 9 For simplicity, only the fundamental frequency was considered at this case. The theoretical fundamental frequencies of the cable at different fixed tension forces are listed in Table 2. All these fundamental frequencies were used as theoretical IFs at different constant tension forces to compare with the extracted IFs under time-varying tension forces. Table 2. The theoretical fundamental frequency of the cable under different constant tension forces. Tension (kn) Frequency (Hz) Tension (kn) Frequency (Hz) Tension (kn) Frequency (Hz) Figure 8. Measured cable tension forces with linear variation. Figure 9. Measured cable acceleration responses with linearly varying tension forces. Figure 10. Wavelet scalogram of cable acceleration responses with linearly varying tension forces.

11 844 Journal of Low Frequency Noise, Vibration and Active Control 37(4) For the first case, the initial tension force of the cable was set to 20 kn, and then the tension force increased linearly at the rate of 1.67 kn/s using the MTS load system. The duration of data acquisition was set to be 6 s. The measured tension force and acceleration responses were shown in Figures 8 and 9, respectively. By performing SWT on acceleration responses with Morlet wavelet as parent function, the wavelet scalogram was obtained and displayed in Figure 10, which indicates that there were two IF trajectories in the time frequency plane. Extended AMD was used to extract the fundamental frequency component and then the normalization scheme was introduced to demodulate it. The results of component decomposition and demodulation are presented in Figure 11, which shows that the periodic change of waveform of the FM signal obtained by the normalization scheme is in good accordance with that of theoretical counterparts. Finally, the improved arccosine Figure 11. Extracted mono-components of cable acceleration responses with linearly varying tension forces. AMD: analytical mode decomposition; EMD: empirical mode decomposition. Table 3. The theoretical fundamental frequency of the cable with linearly varying tension forces. Time (s) Tension (kn) Theoretical frequency (Hz) Time (s) Tension (kn) Theoretical frequency (Hz) Figure 12. Extracted IF of the cable with linearly varying cable tension forces. (a) Compared with HHT and theoretical IFs and (b) compared with theoretical IFs alone. HHT: Hilbert Huang transform.

12 Liu et al. 845 function and moving average algorithm were implemented to identify the instantaneous phases and IFs of the cable and the theoretical fundamental frequencies solved by the time frozen method are shown in Table 3. Figure 12(a) presented the theoretical fundamental frequencies and the extracted IFs using the proposed new method and HHT. It can be seen from Figure 12 that the new method can extract IFs of the cable subjected to linearly increased tension forces more accurately than HHT. To compare the performance of the new method with that of HHT, the values of IA are calculated according to equation (18) and hence presented in Table 4. It can be concluded from Table 4 that the proposed new method is much better than HHT as for the IF extraction of a cable with linearly varying tension forces. IF identification of a cable with sinusoidally varying tension forces In the second case, the initial tension force of the cable was set to 22 kn, and then the tension force varied sinusoidally, which is shown in Figure 13. The acceleration response signal with a duration of 6 s was measured at the middle of the cable and shown in Figure 14. The SWT with Morlet wavelet as parent function was first conducted on the acceleration response signal and the corresponding wavelet scalogram is plotted in Figure 15. As shown in Figure 15(a), the measured acceleration response signal contains several components. Figure 15(b) illustrates that the IF of the fundamental frequency component oscillates sinusoidally. Again, the extended AMD was used to decompose the original multicomponent response signals into several mono-components and the normalization scheme was then used to demodulate the fundamental frequency component. The results of extended AMD and normalization scheme are displayed in Figure 16. It can be seen from Figure 16 that there is no obvious change of the periodicity of waveform before and after the normalization scheme. Finally, the improved arccosine function was employed to extract the IFs of the Table 4. IA of IF extraction with linearly varying cable tension forces (%). Methods New method HHT IA HHT: Hilbert Huang transform; IA: index of accuracy; IF: instantaneous frequency. Figure 13. Measured cable tension forces with sinusoidal variation. Figure 14. Measured cable acceleration responses with sinusoidally varying tension forces.

13 846 Journal of Low Frequency Noise, Vibration and Active Control 37(4) Figure 15. Wavelet scalogram of cable acceleration responses with sinusoidally varying tension forces. Figure 16. Extracted mono-components of cable acceleration responses with sinusoidally varying tension forces. AMD: analytical mode decomposition; EMD: empirical mode decomposition. Table 5. The theoretical fundamental modal frequency of the cable with sinusoidally varying tension forces. Time (s) Tension (kn) Theoretical frequency (Hz) Time (s) Tension (kn) Theoretical frequency (Hz) Figure 17. Extracted IF of the cable with sinusoidally varying tension forces. (a) Compared with HHT and theoretical IFs and (b) compared with theoretical IFs alone. HHT: Hilbert Huang transform.

14 Liu et al. 847 Table 6. IA of IF extraction of the cable with sinusoidal varying tension forces (%). Methods New method HHT IA HHT: Hilbert Huang transform; IA: index of accuracy; IF: instantaneous frequency. cable, which were then smoothed by the moving average algorithm. The theoretical fundamental modal frequencies were solved using the time frozen method and shown in Table 5. The comparison of theoretical fundamental frequencies and identification results presented in Figure 17 indicates that the identified IFs using the new method are in reasonable accordance with theoretical results. However, the classical HHT method cannot effectively extract IFs because of its inherent defects and the noise contamination. Moreover, the IA presented in Table 6 shows again that the new method provides better time frequency resolutions than HHT on IF extraction of a cable with sinusoidal tension forces. Conclusions This paper presents a new IF extraction method for nonstationary response signals in civil engineering structures. The proposed method is a combination of the extended AMD, normalization scheme, and improved arccosine function. The extended AMD is able to decompose a multicomponent signal, which may consist of closely spaced modal or low-frequency components, into mono-component signals. The normalization scheme demodulates the extracted mono-component signals into purely FM signals. Then the improved arccosine function with the moving average algorithm as a smoother is presented to enhance the time frequency resolutions of IF curves. A numerical example of a signal with closely spaced AM FM components and an experimental case study on a cable with time-varying tension forces are employed to verify the effectiveness and accuracy of the proposed method. The results demonstrate that the proposed new method can generate more accurate IF identification results than HHT. As the new method only needs the output response data, therefore, it has the potential to trace the instantaneous properties of time-varying structures under ambient excitations. Acknowledgements We thank Dr Chao Wang of National University of Singapore and Professor Weixin Ren of Hefei University of Technology for providing us the data of a cable test with time-varying tension forces. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (Grants No ), Natural Science Foundation of Fujian Province (Grants No. 2016J05111), the Outstanding Youth Fund of Fujian Agriculture and Forestry University (Grants No. XJQ201728), China Postdoctoral Science Foundation under Grants No. 2018M632561, and the China Scholarship Council (No ). References 1. Lei Y, Lin J, He Z, et al. A review on empirical mode decomposition in fault diagnosis of rotating machinery. Mech Syst Signal Process 2013; 35: Huang NE and Wu Z. A review on Hilbert-Huang transform: method and its applications to geophysical studies. Rev Geophys 2008; 46: RG Griffin DW and Lim JS. Signal estimation from modified short-time Fourier transform. IEEE Trans Acoust Speech Signal Process 1984; 32: Cohen L. Time-frequency distributions: a review. Proc IEEE 1989; 77: Huang NE, Shen Z, Long SR, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis. Proc R Soc Lond Ser A Math Phys Eng Sci 1998; 454:

15 848 Journal of Low Frequency Noise, Vibration and Active Control 37(4) 6. Feng Z, Liang M and Chu F. Recent advances in time-frequency analysis methods for machinery fault diagnosis: a review with application examples. Mech Syst Signal Process 2013; 38: Wu Z and Huang NE. Ensemble empirical mode decomposition: a noise-assisted data analysis method. Adv Adapt Data Anal 2009; 1: Kijewski-Correa T and Kareem A. Nonlinear signal analysis: time-frequency perspectives. J Eng Mech 2007; 133: Wang C, Ren WX and Wang ZC. Instantaneous frequency identification of time-varying structures by continuous wavelet transform. Eng Struct 2013; 52: Huang NE, Wu Z, Shen Z, et al. On instantaneous frequency. Adv Adapt Data Anal 2009; 1: Zheng J, Pan H, Yang S, et al. Adaptive parameterless empirical wavelet transform based time-frequency analysis method and its application to rotor rubbing fault diagnosis. Signal Process 2017; 130: Chen X and Cui B. Efficient modeling of fiber optic gyroscope drift using improved EEMD and extreme learning machine. Signal Process 2016; 128: Hou ZK, Hera A and Shinde A. Wavelet-based structural health monitoring of earthquake excited structures. Comput Aided Civil Eng 2006; 21: Qin Y, Mao Y and Tang B. Multicomponent decomposition by wavelet modulus maxima and synchronous detection. Mech Syst Signal Process 2017; 91: Kijewski T and Kareem A. Wavelet transforms for system identification in civil engineering. Comput Aided Civil Eng 2003; 18: Daubechies I, Lu JF and Wu HT. Synchrosqueezed wavelet transforms: an empirical mode decomposition-like tool. Appl Comput Harmonic Anal 2011; 30: Thakur G, Brevdo E, Fuckar NS, et al. The synchrosqueezing algorithm for time-varying spectral analysis: robustness properties and new paleoclimate applications. Signal Process 2013; 93: Li C and Liang M. A generalized synchrosqueezing transform for enhancing signal time-frequency representation. Signal Process 2012; 92: Cao H, Xi S and Chen X. Zoom synchrosqueezing transform and iterative demodulation: methods with application. Mech Syst Signal Process 2016; 72 73: Chen GD and Wang ZC. A signal decomposition theorem with Hilbert transform and its application to narrowband time series with closely-spaced frequency components. Mech Syst Signal Process 2012; 28: Wang ZC, Ren WX and Liu JL. A synchrosqueezed wavelet transform enhanced by extended analytical mode decomposition method for dynamic signal reconstruction. J Sound Vib 2013; 332: Smith JS. The local mean decomposition and its application to EEG perception data. J R Soc Interface 2005; 2: Liu Z, Jin Y, Zuo MJ, et al. Time-frequency representation based on robust local mean decomposition for multicomponent AM-FM signal analysis. Mech Syst Signal Process 2017; 95: Guo SM, Tsai JSH, Chen CY, et al. An improved empirical mode decomposition based on time scale allocation method and the 2D mode mixing phenomenon judgement. Adv Data Sci Adapt Anal 2017; 9: