Shock and Vibration

Volume 2015, Article ID 846943, 9 pages

http://dx.doi.org/10.1155/2015/846943

## Estimation of Evolutionary Spectra of Monitored Seismic Ground Motions by Transformation of Correlation Functions

School of Civil Engineering and Mechanics, Huazhong University of Science and Technology, Wuhan 430074, China

Received 30 July 2014; Accepted 1 October 2014

Academic Editor: Bo Chen

Copyright © 2015 Zhifeng Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Spatially varying seismic ground motions recorded by distributed structural health monitoring systems (SHMS) can be used to improve the performances of civil engineering structures, necessitating estimation of the evolutionary power spectral density as an indispensable procedure for utilizing records of SHMS. This paper proposes a method for the estimation of evolutionary power spectral density of a nonstationary process by transforming the correlation functions of its sample time histories. First, the background of the theory of evolutionary power spectral density is reviewed in detail. Relationship between the EPSD and the correlation function of a reference stationary process is then established. Formulas are derived for estimating this correlation function directly from nonstationary sample time histories so that the EPSD can be obtained. The implementation procedure of the proposed method is also detailed. Finally, a numerical example is presented, which validates the proposed method, demonstrates its application for SHMS, and displays its capabilities by comparison with the traditional method.

#### 1. Introduction

With the emergence of large-scale infrastructures, such as long-span bridges, skyscrapers, and space structures, structural health monitoring system (SHMS), which aims at enhancing the safety and reliability and reducing the operation and maintenance costs, has been an active research regime [1–5]. One necessary application of the SHMS is to record and analyze structural responses under extreme natural effects such as earthquake and wind storms. Estimating the statistical characteristics of these natural effects is therefore an indispensable procedure for exploiting SHMS. Because these natural effects, especially the seismic ground motion records, usually display nonstationary features, methods for accurately estimating some useful statistical characteristics of nonstationary seismic ground motions are of considerable importance for the application of seismic records given by SHMS equipped on long-span civil engineering structures.

For nonstationary stochastic processes such as seismic ground motions or downburst winds, their time-frequency properties or time-dependent spectra attracted lots of researches in the past fifty years. Some existing approaches for depicting these properties are listed as follows: the evolutionary power spectral density proposed by Priestley [6], which is interpreted as energy distribution over frequency at different time for those can be regarded as quasi-stationary process; the short-time Fourier transform (STFT), obtained by using a sliding window for filtering [7]; the Wigner-Ville distribution (WVD), which is defined as the Fourier transform of an instantaneous correlation function [8]; the wavelet transform (WT) scalogram, which is represented by a group of varied sized wavelets [9]; the Hilbert spectrum acquired by Hilbert Huang transform (HHT), which is represented by Hilbert transform of intrinsic mode function (IMF) components [10]; the Chirplet transform (CT), which can be treated as the extension and combination of STFT and WT [11]. Among them, due to its physical comparability with the power spectral density of stationary processes and its mathematical completeness, the approach of evolutionary power spectral density (EPSD) is most widely used in earthquake engineering for characterizing seismic ground motions.

Although the concept of EPSD has been developed for a relatively long time period, the methods for estimating EPSD from sample time histories of seismic ground motion still remain only partially resolved. A variety of methods have been proposed for estimation of EPSD. Priestley proposed an original EPSD estimating method through windowed filtering [6]. Conte and Peng introduced multiple-windowed STFT method in the EPSD estimation [12]. Spanos and Failla proposed an innovative approach to estimate the EPSD of univariate nonstationary stochastic process via wavelet transform (WT) [13], which was later extended to multivariate cases by Huang and Chen [14]. Chen estimated the EPSD of downburst wind through time-varying autoregressive mode [15]. Kayhan proposed a data-adaptive method for EPSD estimation [16, 17]. Nevertheless, due to the basic difficulty of estimation of EPSD posed by the limited amount of sample time histories and by the Heisenberg uncertainty principle, none of these methods could perfectly be applicable for most situations. In addition, these methods only take advantage of the time series directly but not of the transformation of correlation function, because of the blurred relationship from nonstationary time-varying correlation function to EPSD.

This paper aims to propose an alternative method for EPSD estimation of a set of multivariate sample time histories, which yields the EPSD by transforming their estimated correlation functions. First, the background information of the theory of EPSD and its traditional estimation method is briefly introduced. On top of that, the relationship between the EPSD and the correlation function of a reference stationary process is established. A method is then derived for estimating this correlation function directly from nonstationary sample time histories, on the basis of which the EPSD can be obtained by transformation. Smoothing procedure necessary for eliminating random fluctuations is also presented. Then the implementation procedure of the proposed method of estimation of EPSD is detailed. Finally, a numerical example is presented in order to demonstrate the validity and the application of the proposed method.

#### 2. Background of Evolutionary Power Spectral Density (EPSD)

Nonstationary seismic ground motions recorded by a SHMS installed on long-span structures usually are modelled as a one-dimensional* N*-variates (1D-*N*v) nonstationary stochastic vector process consisting of* N* components , where the superscript () denotes the sample of a nonstationary stochastic vector process. For seismic ground motions, the mean values of are set equal to zero without loss of generality.

The component process of has the following representation [6]:
where is the time-frequency modulating function and
where both and denote Gaussian stochastic processes with orthogonal increments, which satisfy
where (^{∗}) denotes complex conjugation and is the stationary power spectral density of the stationary process to be modulated by 1; and therefore
where is the auto-EPSD of . Accordingly, the evolutionary power spectral density (EPSD) matrix of is given by
The elements of the evolutionary power spectral density matrix are defined as
where is the modulating function and is the cross stationary power spectral density to be modulated.

The time-varying correlation matrix of is given by The relationship between and is shown by the following transformation:

The theory of evolutionary processes proposed by Priestley [6] can be applied only for semistationary processes; that is, the nonstationary process can be regarded as “stationary” in a certain period of time . Generally, this time period is defined as follows [6]: Apparently, the value of cannot be determined, which is a priori.

##### 2.1. Priestley’s Method for Estimation of EPSD

For estimating the EPSD from sample time histories, traditionally the method proposed by Priestley [6] is invoked, which is briefly introduced here. This approach is showed by the following transforms: where is the estimation of EPSD, is a window function, and is a weight function, which are usually selected as

#### 3. EPSD Estimation via Estimated Correlation Function

To utilize the correlation functions to estimate EPSD, firstly, a distinct relationship between EPSD and the correlation function of a reference stationary process of the original nonstationary process is derived, by which the problem is transformed to estimate the correlation function of that specified reference stationary process. Consequently, a set of formulas is derived for a specially designed method to estimate that correlation function from sample time histories and thus enables the estimation of EPSD. The estimated EPSD from a limited number of samples may possess large random fluctuations, and a smoothing procedure is presented to address this. For facilitating the application of the proposed method, its implementation procedure is also presented.

##### 3.1. Transforming Correlation Function to EPSD

With reference to the spectral representation of an evolutionary process in 1, we can define a family of stochastic processes by fixing the time instant as follows: that is, for a specified time instant , is a stochastic process with respect to the pseudo time but not the real time , and can be regarded as stationary since does not vary with the pseudo time . Consequently, defining as the cross power density spectra between the stationary process and the stationary process , with defined as the correlation function between and , the relationship between and can be depicted as a Fourier transformation pair:

In this context, when the correlation function is calculated by estimating over a period* T*, for example, by
then can be calculated by 14. On the other hand, using 13 has the following representation:
By comparing 7 with 17, it can derived that once the pseudo time is replaced with real time , the power spectral density of can be equal to the EPSD of :
Jointly invoking 14–18, by estimation of the correlation function of the reference stationary process , the EPSD of can be readily obtained by Fourier transformation. However, since only the sample time histories of but not of are available, formulas for estimating from need to be developed.

##### 3.2. Formulas for Estimating Correlation Function

Recalling that the evolutionary processes are supposed to be semistationary over a time period of such that , the following approximate relationship can be presumed for EPSD: Utilizing this relationship and taking into consideration the relationship between EPSD and time-varying correlation function of 9, we can derive that Equation 20 indicates that, under the assumption of local stationarity, the correlation function of the reference process can be approximated by the nonstationary time-varying correlation functions in the time period . Therefore, the correlation function can be calculated approximately through the following formula: Consequently, can now be estimated by the Fourier transform of :

##### 3.3. Smoothing

For estimating EPSD of nonstationary processes from the ground motions recorded by SHMS, an inherent difficulty is that the amount sample of time histories is very limited: usually one set of sample time histories is available. This entails considerable part of the stochastic error in the estimated EPSD , which indicates that is highly unreliable and thus not applicable in practice. In order to reduce the fluctuations of , two weight functions and must be included to smooth in the time domain and in the frequency domain so that where represents a smoothed result of , which significantly decreases stochastic errors although it may increase the bias of estimated EPSD. In consideration of the constraints due to limited samples, can be used as an applicable estimation of EPSD.

##### 3.4. Implementation Procedure

In practice, suppose sets of sample time histories of the seismic ground motion recorded by SHMS are presented. For using the proposed method to estimate the EPSD, the following three steps are included.Estimate the correlation function of the reference stationary processes from the given sample time histories:
where , is the time step, and* P* is the total amount of samples.Performing the Fourier transform of , we obtain the unsmooth estimation of EPSD :
where and is the frequency step.Smoothing the unsmoothed estimation of evolutionary spectra in both the time domain and the frequency domain, we obtain the final result:

#### 4. Numerical Example

In this part, a numerical example is presented to demonstrate the capabilities of the proposed approach in estimating the EPSD of seismic ground motions. Due to the lack of multipoint measurements, the SHMS-recorded ground motions are represented by an ensemble of multivariate nonstationary time histories generated conforming to EPSD of a well-known measured seismic ground motion time history and the coherence function given by statistical results of measurements. Then the method proposed in this paper was adopted to estimate the ESPD of sample time histories. The estimated results are compared with the target EPSD of the generated samples as well as the corresponding results yielded by traditional estimation method.

##### 4.1. Sample Time Histories of Seismic Ground Motions

In this example, the ground motions were modeled as multivariate nonstationary stochastic vector processes according to the EPSD of the seismic ground motion record of El-Centro. Figure 1 shows the time history of El-Centro wave. The evolutionary spectrum of the El-Centro wave is denoted by as showed in Figure 2. The ground motions at three points are treated as trivariate nonstationary stochastic vector processes. The configuration of the three points is displayed in Figure 3. For simplicity, the effect of wave propagation is not considered in this example.