# Application of Improved Combined Deterministic-Stochastic Subspace Algorithm in Bridge Modal Parameter Identification

Authors: Peng Wen; Inamullah Khan; Jie He; Qiaofeng Chen
Journal: Shock and Vibration
 (2021)

---

## Abstract

Modal parameter identification is considered to be one of the most important tasks in structural health monitoring because it provides a reliable reference for structural vibration control, damage severity, and operational state. Moreover, at present, the combined deterministic-stochastic subspace algorithm is cogitated as one of the key algorithms in the modal parameter identification, which is why it is widely used in the modal parameter identification of bridge structures. In this paper, a novel method is proposed, which is a time-domain identification algorithm, based on sliding window-fuzzy C-means clustering algorithm-combined with deterministic-stochastic subspace identification (SC-CDSI), to achieve online intelligent tracking and identification of modal parameters for nonlinear time-varying structures. First of all, to realize the online tracking and identification process, it is necessary to divide the input and output signal of the nonlinear time-varying structure by windowing; for that, to determine the window function, window size and window step length according to the characteristics of the signal are analyzed. Secondly, in order to satisfy the intelligent identification of effective modals in stability diagram, the fuzzy C-means clustering algorithm is kept as a base, whereas frequency, damping ratio, and modal shapes serve as clustering elements, applied to fuzzy C-means clustering algorithm, and then the intelligent selection of effective modals is achieved. Finally, a shaking table test bridge is used as a modal parameter identification in lab, and its results are compared with the MIDAS finite element results. The compared results show that the proposed SC-CDSI identification algorithm can accurately achieve the intelligent identification of online tracking of the structural frequency, and the identification results are reliable to be used in real-life bridge structures.


---

## Body

## 1. Introduction

Stochastic subspace identification (SSI) [1] is one of the most widely used modal parameter identification algorithms, which can effectively identify the natural frequencies, damping ratios, and modal shapes of bridge structures under environmental excitation [2]. Along with the popularization and development of SSI algorithm, it has also got some limitations. According to a previous study [3], a corresponding improved algorithm for these limitations in SSI was proposed. However, the input signal of recognition system is environmental excitation (the signal is unknown); thus, a clear incentive-system-response model cannot be established in the decomposition process of the modal parameters and can lead to an expansion of the system's eigenvalue solution space. Thus, false modals can exist in the final recognition results, and thereby the accuracy of the recognition results will eventually be reduced.The combined deterministic-stochastic subspace identification (CDSI) method [4] is established on the basis of SSI algorithm theory. CDSI uses the structure's input signals (excitation) and output signals (response) as the input of the system for parameters identification. They can provide effective remedy for the limitations present in SSI algorithm, which in turn can improve the accuracy of the process. Additionally, the advantages of the CDSI algorithm are that it can fully account the inputs, outputs, and noise of the structure, which is why it is popular among the users. Belleri et al. [5] applied CDSI to concrete frame structures to perform an accurate identification of modal parameters in different damage states under earthquake excitation; Bakir [6] applied CDSI to identify frame structures under earthquake excitation. These recognition results show that this approach can effectively identify the natural frequency of the structure, and the identified modal shapes are reliable. Huang et al. [7] applied CDSI to multistory steel frame shaking table test to identify the different states of modal parameters. Thus, it can be concluded that CDSI can be used to the determine the structural damage state.In addition, Damore et al. [8] adopted the CDSI method to accurately identify the natural vibration frequency of the shaking table test model, which is a 6-layer steel frame.Thai et al. [4] conducted an experiment on the Walnut Creek Bridge and revealed that the combined deterministic-stochastic subspace method could accurately obtain the frequency information of the structure and the identification result was more ideal than that of the stochastic subspace identification.Moaveni and Asgarieh [9] applied the identification method of definite combined deterministic-stochastic subspace to identify the time-varying nonlinear structural system and compared its identification results with the results of wavelet transformation and finite element analysis.Lopes Dos Santos [10] introduced a recursive subspace identification algorithm, which is similar to the bilinear system identification algorithm based on Picard iterative process, and this iterative method effectively overcomes the shortcomings of LPV subspace algorithms.Reynders and Roeck [11] improved the deterministic-stochastic subspace algorithm and applied the new algorithm to identify the modal parameters of the Z24 bridge on the A1 highway between Bern and Zurich. As compared to the traditional algorithm, the new algorithm has the advantage of faster operation rate and gives more accurate identification results.SSI algorithm is widely used in modal parameters identification of bridge structures, but CDSI algorithm is rarely used in modal parameter identification of bridge structures. The main reason is that the core algorithm of CDSI is a linear algorithm [12], which is not suitable for nonlinear time-varying structures like bridge structures. Based on this approach, the sliding window technique [13] is nested within CDSI algorithm to perform the time-varying modal parameter tracking and identification. Additionally, to avoid manual involvement of distinguishing true and false modes, a fuzzy C-means (FCM) algorithm [14] is adopted for intelligent screening of modal parameter results to complete online tracking and intelligent identification of modal parameters.


## 2. Basic Principle of CDSI Algorithm

The flowchart of the CDSI algorithm can be seen in Figure1, where uk is the input signal of the system; B is input matrix; wk represents various types of noise derived from surrounding circumstances; xk is state vector of the system; A is state matrix; C is output matrix; vk is measured noise derived from sensors, D is midway transfer matrix, and yk is output data of the system.Figure 1
Flowchart of the combined deterministic-stochastic system.A combined deterministic-stochastic state equation is obtained and can be written as(1)xk+1=Ax+Buk+wk,yk=Cxk+|Duk+vk.Now, to determine the natural frequency, damping ratio, and modal shape, a stepwise equation will be presented in detail:Step 1: according to the input and output signals, the Hankel matrix is established, whereUP is past input data and Uf is future output data.(2)U02i−1=u0u1⋯uj−1u1u2⋯uj⋯⋯⋯⋯ui−1ui⋯ui+j−2uiui+1⋯ui+j−1ui+1ui+2⋯ui+j⋯⋯⋯⋯u2i−1u2i⋯u2i+j−2=UpUf.Step 2: to determine the state sequence, it can be written as(3)Xi=xixi+1⋯xi+j−2xi+j−1.Step 3: now to set up a linear equation containing a state matrix (A), an input matrix (B), an output matrix (C), and a midway transfer matrix (D), it can be written as (4)X˜i+1Yi+1=ABCDX˜iUi/i+ρwρv.In the above formula,ρw and ρw are residuals, where both are orthogonal to X˜i and X˜i+1.Step 4: now eigenvalue decomposition [15] of matrix A is performed and can be written as (5)A=eAcΔt,λi=eλicΔt,λic=InλiΔt,λic,λ¯ic=−ξiωi±jωi1−ξi2.Step 5: according to formula (6), the natural frequency value fj−fj+1/fj×100%<Δf����, damping ratio (MAC), and modal shape (Φ) of the structure are performed.(6)wi=ai2+bi2rads,fi=ai2+bi22πHz,Φ=CΨ.In the above formula,ψ is the eigenvector matrix of the system, and it is also a complex number.


## 3. Application of Sliding Window Technique

### 3.1. Addition of Windows Analysis

In the process of modal parameter identification of the response signal, the efficiency of parameter identification can be improved by controlling the length of the signal; therefore, in practical application, the length of the signal can be controlled by windowing truncation processing. Additionally, the length of signal can be controlled by windowing and truncating process. However, the window functions will directly correspond to the signal’s spectrum in the process on windowing. Different window forms cause different levels of spectrum leakage, and that is why the final frequency resolution effect will be different. Based on these conditions, it is particularly important to choose the type of window functions, the size of the windows, and the step length of sliding window.


### 3.2. Window Function Selection

Currently, commonly used window functions [16] are rectangular window, Hanning window, Hamming window, flat-topped window, Gaussian window, and so forth. Comparing the advantages and disadvantages of each window function, a rectangular window was selected as the windowing function of the response signal in this paper. Since the rectangular window suits better for our analysis due to its simplicity and accuracy in frequency identification, it is used in this paper. In addition, rectangular window is a zero-power window of time variable; it also has the advantage of main lobe concentration and can also accurately identify the frequency results, but its accuracy of amplitude identification is a little bit low, which does not have much impact on the analysis of this paper. Therefore, if it is to be only used for accurate identification of frequency main lobe and does not require the accuracy of amplitude recognition, which is the case in a practical application, such as identifying the natural frequency of the structure, a rectangular function can be preferably selected for windowing of signal. The rectangular window function can be written as(7)wt=1.


### 3.3. Window Size Determination

The window size and the quantity of initial data modus operand once entered into system each time will definitely affect the accuracy of the recognition results up to certain extent. If the window size is too small, the accurate stability diagram cannot be identified due to insufficient input signals and thus will easily result in mode omissions. However, if window size is too large, there will be too many false modes into the stability diagram because of too many input signals and will thereby reduce the efficiency of parameters identification process. Based on this scenario, some scholars [17] believe that the characteristics of the structure itself, the sampling frequency of the signal, and the noise level of the signal should be comprehensively considered. However, according to a previous study, it is not feasible to comprehensively determine the size of the window based on these factors, because there are too many factors involved, and each structure has its own unique characteristics. Therefore, the size of window can be defined by referring to the principle of order determination of stability diagram. For that, specific analysis steps are required, which are as follows:(1)
The true orderNof the structure is determined based on the order determination [18] of the stability diagram in the stochastic subspace method.(2)
Assume that the initial window data length isL1, corresponding window time is t1, and sampling frequency is cf; the relationship among these three is L1=t1×cf.(3)
It is assumed that the lag time between previous and latter windows isT; that is, the data length corresponding to the ith window is Li=t1+i−1×T×cf.(4)
The corresponding stability diagrams of each window are identified by CDSI algorithm; in case the system order is N, the numbers of corresponding stable pointsFi=fi1,fi2,fi3 to the first 3 modes of each stability diagram are to be calculated, because the order of the system is N. The maximum number of stable points corresponding to each mode is N/2; thus the percentages BFi=2fi1/N,2fi2/N,2fi3/N of the first three modes of stable points can be calculated.(5)
TheBFi numerical value changes with the signal length; due to this limitation, the percentage shows a gradual increase at first and then keeps a steady change, as shown in Figure 2. It is considered that once the percentage exceeds 90%, the corresponding signal length is considered to be a suitable window size.Figure 2
Percentage trend chart of stability points.


### 3.4. Window Step Size

Currently, the window size is determined, and still further work is required to determine the sliding step size of the window, that is, the data points moving forward by each window to achieve real-time update of the input signal. Assuming that the quantity of updated data of signal isV, the relationship between the data replacement of each window and the input and output of the system are shown in Figure 3.Figure 3
Schematic diagram of step length.According to the definition of the Hankel matrix [19], the length of past input signal is equal to the length of future output signal. Based on this situation, Vwill be defined as the data quantity corresponding to 1/2 window; that is, the overlap between windows is 1/2 window.


## 3.1. Addition of Windows Analysis

In the process of modal parameter identification of the response signal, the efficiency of parameter identification can be improved by controlling the length of the signal; therefore, in practical application, the length of the signal can be controlled by windowing truncation processing. Additionally, the length of signal can be controlled by windowing and truncating process. However, the window functions will directly correspond to the signal’s spectrum in the process on windowing. Different window forms cause different levels of spectrum leakage, and that is why the final frequency resolution effect will be different. Based on these conditions, it is particularly important to choose the type of window functions, the size of the windows, and the step length of sliding window.


## 3.2. Window Function Selection

Currently, commonly used window functions [16] are rectangular window, Hanning window, Hamming window, flat-topped window, Gaussian window, and so forth. Comparing the advantages and disadvantages of each window function, a rectangular window was selected as the windowing function of the response signal in this paper. Since the rectangular window suits better for our analysis due to its simplicity and accuracy in frequency identification, it is used in this paper. In addition, rectangular window is a zero-power window of time variable; it also has the advantage of main lobe concentration and can also accurately identify the frequency results, but its accuracy of amplitude identification is a little bit low, which does not have much impact on the analysis of this paper. Therefore, if it is to be only used for accurate identification of frequency main lobe and does not require the accuracy of amplitude recognition, which is the case in a practical application, such as identifying the natural frequency of the structure, a rectangular function can be preferably selected for windowing of signal. The rectangular window function can be written as(7)wt=1.


## 3.3. Window Size Determination

The window size and the quantity of initial data modus operand once entered into system each time will definitely affect the accuracy of the recognition results up to certain extent. If the window size is too small, the accurate stability diagram cannot be identified due to insufficient input signals and thus will easily result in mode omissions. However, if window size is too large, there will be too many false modes into the stability diagram because of too many input signals and will thereby reduce the efficiency of parameters identification process. Based on this scenario, some scholars [17] believe that the characteristics of the structure itself, the sampling frequency of the signal, and the noise level of the signal should be comprehensively considered. However, according to a previous study, it is not feasible to comprehensively determine the size of the window based on these factors, because there are too many factors involved, and each structure has its own unique characteristics. Therefore, the size of window can be defined by referring to the principle of order determination of stability diagram. For that, specific analysis steps are required, which are as follows:(1)
The true orderNof the structure is determined based on the order determination [18] of the stability diagram in the stochastic subspace method.(2)
Assume that the initial window data length isL1, corresponding window time is t1, and sampling frequency is cf; the relationship among these three is L1=t1×cf.(3)
It is assumed that the lag time between previous and latter windows isT; that is, the data length corresponding to the ith window is Li=t1+i−1×T×cf.(4)
The corresponding stability diagrams of each window are identified by CDSI algorithm; in case the system order is N, the numbers of corresponding stable pointsFi=fi1,fi2,fi3 to the first 3 modes of each stability diagram are to be calculated, because the order of the system is N. The maximum number of stable points corresponding to each mode is N/2; thus the percentages BFi=2fi1/N,2fi2/N,2fi3/N of the first three modes of stable points can be calculated.(5)
TheBFi numerical value changes with the signal length; due to this limitation, the percentage shows a gradual increase at first and then keeps a steady change, as shown in Figure 2. It is considered that once the percentage exceeds 90%, the corresponding signal length is considered to be a suitable window size.Figure 2
Percentage trend chart of stability points.


## 3.4. Window Step Size

Currently, the window size is determined, and still further work is required to determine the sliding step size of the window, that is, the data points moving forward by each window to achieve real-time update of the input signal. Assuming that the quantity of updated data of signal isV, the relationship between the data replacement of each window and the input and output of the system are shown in Figure 3.Figure 3
Schematic diagram of step length.According to the definition of the Hankel matrix [19], the length of past input signal is equal to the length of future output signal. Based on this situation, Vwill be defined as the data quantity corresponding to 1/2 window; that is, the overlap between windows is 1/2 window.


## 4. Application of FCM Clustering Algorithm

SSI algorithm and CDSI algorithm involve manual participation to select effective modals [20] in stability diagram. To avoid subjectivity caused by manual participation in effective modal parameter identification, an intelligent screening method for effectiveness modal selection based on FCM clustering algorithm is proposed. The basic principle and steps of FCM clustering algorithm and then the intelligent identification of effective modes will be determined by combining the FCM algorithm with CDSI algorithm.


### 4.1. The Basic Principle of FCM Clustering Algorithm

Compared with K-means clustering and pedigree clustering algorithm, FCM clustering [21] has many advantages.The clustering principle of FCM is based on the degree of membership of the data points to identify its classification; therefore, it is a type of a fuzzy clustering algorithm [22] and can perform clustering of different dimensions of data. To achieve this, the data set Xis divided into class c using FCM algorithm, assuming that X consists of N groups sample; that is, X=x1,x2,…,xN, where each sample consists of n-dimensional vectors; that is,xj=xj1,xj2,…,xjn. The membership matrix U=uij of the data setsX is defined as a matrix of C×N, where C is the number of classifications, and uij is the jth sample in X; the membership belongs to grade ith, and it must meet the following conditions:(8)∑i=1cuij=1,∀j=1,2,…,N,0≤uij≤1,1≤i≤c,1≤j≤N,whereas the formulation of uij is performed as follows:(9)uij=1∑k=1cdij/dki2/m−1,i=1,2,…,c;j=1,2,…,n.Fuzzy indexm∈1,∞ and dij∈xj−vi are represented in the Euclidean distance between the data object xj and clustering center vi of ith cluster, where v1,v2,…,vi represents the clustering center set and its formulation is performed as follows: (10)vi=∑j=1nuijm×xj∑j=1nuijm.The objective function value of FCM algorithm isJu,v; its purpose of iterative process is to find the membership matrix and clustering center to minimize the function Ju,v.(11)Ju,v=∑j=1n∑i=1cuijmd2xj,vi.The specific steps of FCM clustering algorithm are as follows:(1)
The amount of classificationc and the maximum number of iterations m (fuzzy index) are determined,(2)
Randomly select a data set from the data setsX as the center point of the initial clustering,(3)
The membership matrixU is performed according to formula (9),(4)
The clustering centerv is performed according to formula (10),(5)
According to formula (11), membership matrix and clustering center are calculated, so that the function Ju,v has the minimum value and thus the clustering process got completed.


### 4.2. Combination of FCM Algorithm and CDSI Algorithm

To effectually realize the intelligent identification of modal in the stability diagram, a new parameter identification algorithm (C-CDSI) is proposed by combining FCM algorithm and CDSI algorithm. The specific steps for achieving such algorithm are as follows.


#### 4.2.1. Identification of Parameters of Each Window

The parameter resultsXi=Fi,wi,Φi corresponding to each window are identified by CDSI algorithm, where F is frequency value, w is damping ratio, and Φ is modal shape. Each parameter result is a matrix with m×n, where the value m is equal to 1/2 of the real order.(12)Fi=fim,n,wi=ξim,n,Φi=ψim,n.


#### 4.2.2. Frequency Clustering

Accordingly, the frequency, damping ratio, and modal shape ofXi=Fi,wi,Φi and Xj=Fj,wj,Φj are clustered by FCM algorithm. In order to achieve this, the following steps will be followed:(1)
The clustering data set is determined, and it is assumed that the modal frequencies off1m,n and f2m,n are clustered.(2)
The number of clustering classifications is determined; because there are two types of modal frequency of the data set for each clustering, the number of classifications is set asc=2.(3)
The maximum iterations are determined; because the rows in parameter results matrix are true orderN of 1/2 of the system, the maximum iteration is set as m=N/2.(4)
The first-order frequency resultsf1m,1 in f1 and the first-order frequency results f2m,1 in f2 are clustered. If these are in the same category, the average value will be the new first-order frequency result; that is, f1−2m,1=f1m,1+f2m,1/2. All of these modes in the frequency matrix will be clustered and so on. Then, a new frequency matrix f1−2 will be established by the clustered new modes with that of the nonclustered modes.


#### 4.2.3. Damping Ratio Clustering

The damping ratio clustering identification of corresponding two-order modes inξ1m,n and ξ2m,n is performed. To verify the authenticity of the clustered two-order modes in frequency matrix f1−2, the following steps will be followed:(1)
Assuming thatf1m,1 and f2m,1 in f1−2 are similar items, the clustering is completed based on Section 4.2.2.(2)
If these are similar terms, the average value will be a new first-order damping ratio result;ξ1−2m,1=ξ1m,1+ξ2m,1/2.(3)
If these are not similar terms,f1−2m,1 will be split into the initial items f1m,1 and f2m,1 , and then a new frequency matrix q1−2 will be reconstructed.(4)
A new damping ratio matrixξ1−2 is formed from all the clustered damping ratios and the nonclustering damping ratios.


#### 4.2.4. Modal Mode Clustering

The modal shape clustering of the corresponding two-order modes inψ1m,n and ψ2m,n is performed; to verify the authenticity of the clustering results of frequency cluster in q1−2 and damping ratio cluster in ψ2m,n, the following steps will be followed:(1)
Assuming thatf1m,1 and f2m,1 in q1−2 are similar items, the clustering identification of ψ1m,1 and ψ2m,1 is completed based on Section 4.2.2.(2)
If these are similar terms, the average value will be a new first-order modal shape result;ψ1−2m,1=ψ1m,1+ψ2m,1/2.(3)
If these are not similar terms,f1−2m,1 will be split into the initial items f1m,1 and f2m,1, and then a new frequency matrix Q1−2 will be reconstructed. Additionally, ξ1−2m,1 will be split into the initial items ξ1m,1 and ξ2m,1, and then a new damping ratio matrix ζ1−2 will be reconstructed.(4)
A new modal shape matrixψ1−2 is formed from all the clustered modal shapes and the nonclustering modal shapes.


#### 4.2.5. Clustering of Each Window Parameter Successively

A new set of parameter resultsX1−2=Q1−2,ξ1−2,ψ1−2 are established by Q1−2, ζ1−2, and ψ1−2, and the same modal cluster between X1−2 and X3 based on Section 4.2.2 to Section 4.2.4 is completed to obtain a new clustering result X2−3. Based on this approach, the parameter result clustering among all windows is hereby termed as completed. Assuming that the total number of windows is N, the final clustering parameter result will be XN−1−N.The effective modal parameters identification (the frequency, damping ratio, and modal shape are stable modals) inXN−1−N can be performed by counting each cluster modal of XN−1−N, which is made up of the corresponding modals of windows. If modal’s clustering window is greater than 0.8N, it will be an effectiveness modal and will be drawn in stability diagram.Based on the above descriptions, a new modal parameters online intelligent identification algorithm (SC-CDSI) that combines sliding window technology, FCM clustering algorithm, and CDSI algorithm is proposed. The specific steps of this method are shown in the flowchart diagram in Figure4.Figure 4
Flowchart of SC-CDSI.


## 4.1. The Basic Principle of FCM Clustering Algorithm

Compared with K-means clustering and pedigree clustering algorithm, FCM clustering [21] has many advantages.The clustering principle of FCM is based on the degree of membership of the data points to identify its classification; therefore, it is a type of a fuzzy clustering algorithm [22] and can perform clustering of different dimensions of data. To achieve this, the data set Xis divided into class c using FCM algorithm, assuming that X consists of N groups sample; that is, X=x1,x2,…,xN, where each sample consists of n-dimensional vectors; that is,xj=xj1,xj2,…,xjn. The membership matrix U=uij of the data setsX is defined as a matrix of C×N, where C is the number of classifications, and uij is the jth sample in X; the membership belongs to grade ith, and it must meet the following conditions:(8)∑i=1cuij=1,∀j=1,2,…,N,0≤uij≤1,1≤i≤c,1≤j≤N,whereas the formulation of uij is performed as follows:(9)uij=1∑k=1cdij/dki2/m−1,i=1,2,…,c;j=1,2,…,n.Fuzzy indexm∈1,∞ and dij∈xj−vi are represented in the Euclidean distance between the data object xj and clustering center vi of ith cluster, where v1,v2,…,vi represents the clustering center set and its formulation is performed as follows: (10)vi=∑j=1nuijm×xj∑j=1nuijm.The objective function value of FCM algorithm isJu,v; its purpose of iterative process is to find the membership matrix and clustering center to minimize the function Ju,v.(11)Ju,v=∑j=1n∑i=1cuijmd2xj,vi.The specific steps of FCM clustering algorithm are as follows:(1)
The amount of classificationc and the maximum number of iterations m (fuzzy index) are determined,(2)
Randomly select a data set from the data setsX as the center point of the initial clustering,(3)
The membership matrixU is performed according to formula (9),(4)
The clustering centerv is performed according to formula (10),(5)
According to formula (11), membership matrix and clustering center are calculated, so that the function Ju,v has the minimum value and thus the clustering process got completed.


## 4.2. Combination of FCM Algorithm and CDSI Algorithm

To effectually realize the intelligent identification of modal in the stability diagram, a new parameter identification algorithm (C-CDSI) is proposed by combining FCM algorithm and CDSI algorithm. The specific steps for achieving such algorithm are as follows.


### 4.2.1. Identification of Parameters of Each Window

The parameter resultsXi=Fi,wi,Φi corresponding to each window are identified by CDSI algorithm, where F is frequency value, w is damping ratio, and Φ is modal shape. Each parameter result is a matrix with m×n, where the value m is equal to 1/2 of the real order.(12)Fi=fim,n,wi=ξim,n,Φi=ψim,n.


### 4.2.2. Frequency Clustering

Accordingly, the frequency, damping ratio, and modal shape ofXi=Fi,wi,Φi and Xj=Fj,wj,Φj are clustered by FCM algorithm. In order to achieve this, the following steps will be followed:(1)
The clustering data set is determined, and it is assumed that the modal frequencies off1m,n and f2m,n are clustered.(2)
The number of clustering classifications is determined; because there are two types of modal frequency of the data set for each clustering, the number of classifications is set asc=2.(3)
The maximum iterations are determined; because the rows in parameter results matrix are true orderN of 1/2 of the system, the maximum iteration is set as m=N/2.(4)
The first-order frequency resultsf1m,1 in f1 and the first-order frequency results f2m,1 in f2 are clustered. If these are in the same category, the average value will be the new first-order frequency result; that is, f1−2m,1=f1m,1+f2m,1/2. All of these modes in the frequency matrix will be clustered and so on. Then, a new frequency matrix f1−2 will be established by the clustered new modes with that of the nonclustered modes.


### 4.2.3. Damping Ratio Clustering

The damping ratio clustering identification of corresponding two-order modes inξ1m,n and ξ2m,n is performed. To verify the authenticity of the clustered two-order modes in frequency matrix f1−2, the following steps will be followed:(1)
Assuming thatf1m,1 and f2m,1 in f1−2 are similar items, the clustering is completed based on Section 4.2.2.(2)
If these are similar terms, the average value will be a new first-order damping ratio result;ξ1−2m,1=ξ1m,1+ξ2m,1/2.(3)
If these are not similar terms,f1−2m,1 will be split into the initial items f1m,1 and f2m,1 , and then a new frequency matrix q1−2 will be reconstructed.(4)
A new damping ratio matrixξ1−2 is formed from all the clustered damping ratios and the nonclustering damping ratios.


### 4.2.4. Modal Mode Clustering

The modal shape clustering of the corresponding two-order modes inψ1m,n and ψ2m,n is performed; to verify the authenticity of the clustering results of frequency cluster in q1−2 and damping ratio cluster in ψ2m,n, the following steps will be followed:(1)
Assuming thatf1m,1 and f2m,1 in q1−2 are similar items, the clustering identification of ψ1m,1 and ψ2m,1 is completed based on Section 4.2.2.(2)
If these are similar terms, the average value will be a new first-order modal shape result;ψ1−2m,1=ψ1m,1+ψ2m,1/2.(3)
If these are not similar terms,f1−2m,1 will be split into the initial items f1m,1 and f2m,1, and then a new frequency matrix Q1−2 will be reconstructed. Additionally, ξ1−2m,1 will be split into the initial items ξ1m,1 and ξ2m,1, and then a new damping ratio matrix ζ1−2 will be reconstructed.(4)
A new modal shape matrixψ1−2 is formed from all the clustered modal shapes and the nonclustering modal shapes.


### 4.2.5. Clustering of Each Window Parameter Successively

A new set of parameter resultsX1−2=Q1−2,ξ1−2,ψ1−2 are established by Q1−2, ζ1−2, and ψ1−2, and the same modal cluster between X1−2 and X3 based on Section 4.2.2 to Section 4.2.4 is completed to obtain a new clustering result X2−3. Based on this approach, the parameter result clustering among all windows is hereby termed as completed. Assuming that the total number of windows is N, the final clustering parameter result will be XN−1−N.The effective modal parameters identification (the frequency, damping ratio, and modal shape are stable modals) inXN−1−N can be performed by counting each cluster modal of XN−1−N, which is made up of the corresponding modals of windows. If modal’s clustering window is greater than 0.8N, it will be an effectiveness modal and will be drawn in stability diagram.Based on the above descriptions, a new modal parameters online intelligent identification algorithm (SC-CDSI) that combines sliding window technology, FCM clustering algorithm, and CDSI algorithm is proposed. The specific steps of this method are shown in the flowchart diagram in Figure4.Figure 4
Flowchart of SC-CDSI.


## 4.2.1. Identification of Parameters of Each Window

The parameter resultsXi=Fi,wi,Φi corresponding to each window are identified by CDSI algorithm, where F is frequency value, w is damping ratio, and Φ is modal shape. Each parameter result is a matrix with m×n, where the value m is equal to 1/2 of the real order.(12)Fi=fim,n,wi=ξim,n,Φi=ψim,n.


## 4.2.2. Frequency Clustering

Accordingly, the frequency, damping ratio, and modal shape ofXi=Fi,wi,Φi and Xj=Fj,wj,Φj are clustered by FCM algorithm. In order to achieve this, the following steps will be followed:(1)
The clustering data set is determined, and it is assumed that the modal frequencies off1m,n and f2m,n are clustered.(2)
The number of clustering classifications is determined; because there are two types of modal frequency of the data set for each clustering, the number of classifications is set asc=2.(3)
The maximum iterations are determined; because the rows in parameter results matrix are true orderN of 1/2 of the system, the maximum iteration is set as m=N/2.(4)
The first-order frequency resultsf1m,1 in f1 and the first-order frequency results f2m,1 in f2 are clustered. If these are in the same category, the average value will be the new first-order frequency result; that is, f1−2m,1=f1m,1+f2m,1/2. All of these modes in the frequency matrix will be clustered and so on. Then, a new frequency matrix f1−2 will be established by the clustered new modes with that of the nonclustered modes.


## 4.2.3. Damping Ratio Clustering

The damping ratio clustering identification of corresponding two-order modes inξ1m,n and ξ2m,n is performed. To verify the authenticity of the clustered two-order modes in frequency matrix f1−2, the following steps will be followed:(1)
Assuming thatf1m,1 and f2m,1 in f1−2 are similar items, the clustering is completed based on Section 4.2.2.(2)
If these are similar terms, the average value will be a new first-order damping ratio result;ξ1−2m,1=ξ1m,1+ξ2m,1/2.(3)
If these are not similar terms,f1−2m,1 will be split into the initial items f1m,1 and f2m,1 , and then a new frequency matrix q1−2 will be reconstructed.(4)
A new damping ratio matrixξ1−2 is formed from all the clustered damping ratios and the nonclustering damping ratios.


## 4.2.4. Modal Mode Clustering

The modal shape clustering of the corresponding two-order modes inψ1m,n and ψ2m,n is performed; to verify the authenticity of the clustering results of frequency cluster in q1−2 and damping ratio cluster in ψ2m,n, the following steps will be followed:(1)
Assuming thatf1m,1 and f2m,1 in q1−2 are similar items, the clustering identification of ψ1m,1 and ψ2m,1 is completed based on Section 4.2.2.(2)
If these are similar terms, the average value will be a new first-order modal shape result;ψ1−2m,1=ψ1m,1+ψ2m,1/2.(3)
If these are not similar terms,f1−2m,1 will be split into the initial items f1m,1 and f2m,1, and then a new frequency matrix Q1−2 will be reconstructed. Additionally, ξ1−2m,1 will be split into the initial items ξ1m,1 and ξ2m,1, and then a new damping ratio matrix ζ1−2 will be reconstructed.(4)
A new modal shape matrixψ1−2 is formed from all the clustered modal shapes and the nonclustering modal shapes.


## 4.2.5. Clustering of Each Window Parameter Successively

A new set of parameter resultsX1−2=Q1−2,ξ1−2,ψ1−2 are established by Q1−2, ζ1−2, and ψ1−2, and the same modal cluster between X1−2 and X3 based on Section 4.2.2 to Section 4.2.4 is completed to obtain a new clustering result X2−3. Based on this approach, the parameter result clustering among all windows is hereby termed as completed. Assuming that the total number of windows is N, the final clustering parameter result will be XN−1−N.The effective modal parameters identification (the frequency, damping ratio, and modal shape are stable modals) inXN−1−N can be performed by counting each cluster modal of XN−1−N, which is made up of the corresponding modals of windows. If modal’s clustering window is greater than 0.8N, it will be an effectiveness modal and will be drawn in stability diagram.Based on the above descriptions, a new modal parameters online intelligent identification algorithm (SC-CDSI) that combines sliding window technology, FCM clustering algorithm, and CDSI algorithm is proposed. The specific steps of this method are shown in the flowchart diagram in Figure4.Figure 4
Flowchart of SC-CDSI.


## 5. Shaking Table Test Bridge

In order to verify the reliability of the improved combined deterministic-stochastic subspace algorithm (SC-CDSI) as mentioned above, a shaking table test of large cable-stayed bridge was performed to analyze its modal parameters variations against the earthquake excitations.


### 5.1. General Introduction of Model Bridge

The model bridge tested in lab was a prototype of actual cable-stayed bridge (scaled down to a ratio of 1 : 20). The test bridge span in the laboratory was 6.5 + 19 + 6.5 m. The bridge type, its structure layout, floor plan of main girder, and lab model images are shown in Figure5; as far as the construction materials of lab bridge are concerned, cable tower and its piers were manufactured from normal reinforcement (round steel) with a diameter of 6 mm and M15 microconcrete, respectively, shear reinforcements were made from 10# steel wires, both main girder and steel box girder were manufactured from a steel plate with a thickness of 5 mm, and the cable is manufactured from a wire rope with a diameter of 10 mm.Figure 5
Bridge model. (a) Elevation and plan (unit: cm). (b) Reality image of bridge.
(a)(b)


### 5.2. Location of Acceleration Sensor (Accelerometer)

To perform shaking table test for this bridge, four acceleration sensors were installed at the bottom of foundation to monitor the input signals of structure, and eleven accelerometers were installed on the main girder to collect the response signals of structure, as shown in Figure6.Figure 6
Arrangement of acceleration sensors (unit: cm).The sensors used in the cable-stayed bridge are piezoelectric accelerometers, because its working principle is to measure the vibration accelerations by utilizing the characteristics that some special sensitive cores will generate a charge signal proportional to the vibration acceleration. The characteristics of these sensors include the following: wide measuring frequency range, large measuring range, small size, light weight, simple and solidity structure, less external interference, and no need for any external power supply to generate charge signals.


### 5.3. Earthquake Input

The inputs of earthquake excitation are shown in Table1. It can be noticed that the peak acceleration is 0.1 g, the direction of earthquake wave is transverse excitation, and the types of waveforms are Chichi wave and site wave. The acceleration time-history curve of Chichi earthquake wave in C1 working condition is shown in Figure 7. The sampling frequency of sensors is 256 Hz.Table 1
Working conditions of test table.
No.NameDirection of earthquake waveEarthquake waveformPeak acceleration (g)1C1Transverse excitationChichi0.12S1Transverse excitationSite wave0.1Figure 7
Time-history curve of acceleration (C1).


### 5.4. MIDAS Finite Element Analysis Results

A simulation analysis model of the bridge structure was established using MIDAS 2017. Bridge’s piers, towers, and main girders are all simulated using nonlinear beam-column elements, and cables are simulated using truss elements; bearings are all pot-type rubber bearings, and zero-length elements are adopted. The bridge model is shown in Figure8. The eigenvalue analysis of the bridge was performed, and its first five frequencies are obtained, as shown in Table 2. Figure 9 shows the modal shapes from the first mode to fifth mode.Figure 8
Bridge model (MIDAS).Table 2
Natural frequency (unit: Hz).
OrdersThe first orderThe second orderThe third orderThe fourth orderThe fifth orderFrequency1.142.513.774.638.83Figure 9
Modal shape of the first five orders (MIDAS). (a) Vertical bending mode of the first order (f1 = 1.14 Hz). (b) Vertical bending of the second order (f2 = 2.51 Hz). (c) Vertical bending of the third order (f3 = 3.77 Hz). (d) Vertical bending of the fourth order (f4 = 4.63 Hz). (e) Vertical bending of the fifth order (f5 = 8.83 Hz).
(a)(b)(c)(d)(e)


### 5.5. Parameter Identification of Each Window

According to the sliding window technology and CDSI algorithm, the modal parameters of the cable-stayed bridge under two types of earthquake excitations are identified. Assume that the initial window time is 1 s, but the appropriate window size is given to be 4 s; additionally, the corresponding sliding window step size is 2 s, which means that window data is normally updated after every 2 s. The signals of each window are identified, and the stability results corresponding to each window in two types of working conditions are obtained. Figure10 shows the results of the first two stability diagrams in C1 working condition.Figure 10
Stability diagram (working condition: C1). (a) Stability diagram (time: 0 s–4 s). (b) Stability diagram (time: 2 s–6 s).
(a)(b)The specific values of the frequency in each window in operating conditions of C1 and S1 are identified, and the trend chart of the first five changing frequency values with time is given as shown in Figure11. Analyzing the frequency changes of each order in the figure, it is shown that the frequency value of the first five orders of bridge structure is stable in both conditions; that is, the bridge structure’s natural frequency is stable under the action of two types of earthquake waves. It is also shown that the cable-stayed bridge shows no damage in first two types of conditions.Figure 11
Frequency trend chart. (a) Frequency trend chart (working condition: C1). (b) Frequency trend chart (working condition: S1).
(a)(b)


### 5.6. Intelligent Identification of Modal Effectiveness

According to the proposed FCM clustering algorithm inSection 3.2, a clustering analysis of each window’s parameter results which is under two types of conditions is performed to reach intelligent identification of modal effectiveness in the stability diagram. Figure 12 is the corresponding effective stability diagram of C1 and S1 working conditions.Figure 12
Effective stability diagram. (a) Effective stability diagram (working condition: C1). (b) Effective stability diagram (working condition: S1).
(a)(b)Comparing the MIDAS finite element results with the first five orders of frequency results which are under the conditions of C1 and S1, as shown in Table3, the percent difference between the two results in the table and the errors between the effectiveness modal frequency values obtained by the SC-CDSI algorithm and the finite element results difference are within 7%; thus, frequency identification results of the proposed SC-CDSI algorithm are hereby considered reliable.Table 3
Comparative analysis results of frequency.
ResultsThe first orderThe second orderThe third orderThe fourth orderThe fifth orderMIDAS results (Hz)1.042.563.774.548.43Working condition: C1 (Hz)1.112.443.874.688.74Percentage difference 1 (%)−6.34.9−2.6−3.0−3.5Working condition: S1 (Hz)1.092.423.894.728.84Percentage difference 2 (%)−4.65.8−3.1−3.8−4.6


### 5.7. Modal Shape Results of First Three Orders

Based on these approaches, the SC-CDSI algorithm not only has verified accurately the true modes in the stability diagram but also has identified the modal shapes of structure accurately. Additionally, the first three mode shapes of cable-stayed bridge are obtained, as shown in Figure13, the X-coordinate in the figure represents the 11 measurement points of sensor, and the Y-coordinate is the mode shape after 0-1 planning; the dotted line is the identification result of the proposed algorithm; the solid line is the mode shape diagram, which is obtained by MIDAS. By comparing these shapes to the first three mode shapes selected from Figure 9, it can be seen that the modal modes and finite element results obtained by the SC-CDSI algorithm have more than 90% similarity, which is why the modal shapes obtained by SC-CDSI algorithm are reliable.Figure 13
Modal shape of the first three orders. (a) Vibration diagram of the first order. (b) Vibration diagram of the second order. (c) Vibration diagram of the third order.
(a)(b)(c)


## 5.1. General Introduction of Model Bridge

The model bridge tested in lab was a prototype of actual cable-stayed bridge (scaled down to a ratio of 1 : 20). The test bridge span in the laboratory was 6.5 + 19 + 6.5 m. The bridge type, its structure layout, floor plan of main girder, and lab model images are shown in Figure5; as far as the construction materials of lab bridge are concerned, cable tower and its piers were manufactured from normal reinforcement (round steel) with a diameter of 6 mm and M15 microconcrete, respectively, shear reinforcements were made from 10# steel wires, both main girder and steel box girder were manufactured from a steel plate with a thickness of 5 mm, and the cable is manufactured from a wire rope with a diameter of 10 mm.Figure 5
Bridge model. (a) Elevation and plan (unit: cm). (b) Reality image of bridge.
(a)(b)


## 5.2. Location of Acceleration Sensor (Accelerometer)

To perform shaking table test for this bridge, four acceleration sensors were installed at the bottom of foundation to monitor the input signals of structure, and eleven accelerometers were installed on the main girder to collect the response signals of structure, as shown in Figure6.Figure 6
Arrangement of acceleration sensors (unit: cm).The sensors used in the cable-stayed bridge are piezoelectric accelerometers, because its working principle is to measure the vibration accelerations by utilizing the characteristics that some special sensitive cores will generate a charge signal proportional to the vibration acceleration. The characteristics of these sensors include the following: wide measuring frequency range, large measuring range, small size, light weight, simple and solidity structure, less external interference, and no need for any external power supply to generate charge signals.


## 5.3. Earthquake Input

The inputs of earthquake excitation are shown in Table1. It can be noticed that the peak acceleration is 0.1 g, the direction of earthquake wave is transverse excitation, and the types of waveforms are Chichi wave and site wave. The acceleration time-history curve of Chichi earthquake wave in C1 working condition is shown in Figure 7. The sampling frequency of sensors is 256 Hz.Table 1
Working conditions of test table.
No.NameDirection of earthquake waveEarthquake waveformPeak acceleration (g)1C1Transverse excitationChichi0.12S1Transverse excitationSite wave0.1Figure 7
Time-history curve of acceleration (C1).


## 5.4. MIDAS Finite Element Analysis Results

A simulation analysis model of the bridge structure was established using MIDAS 2017. Bridge’s piers, towers, and main girders are all simulated using nonlinear beam-column elements, and cables are simulated using truss elements; bearings are all pot-type rubber bearings, and zero-length elements are adopted. The bridge model is shown in Figure8. The eigenvalue analysis of the bridge was performed, and its first five frequencies are obtained, as shown in Table 2. Figure 9 shows the modal shapes from the first mode to fifth mode.Figure 8
Bridge model (MIDAS).Table 2
Natural frequency (unit: Hz).
OrdersThe first orderThe second orderThe third orderThe fourth orderThe fifth orderFrequency1.142.513.774.638.83Figure 9
Modal shape of the first five orders (MIDAS). (a) Vertical bending mode of the first order (f1 = 1.14 Hz). (b) Vertical bending of the second order (f2 = 2.51 Hz). (c) Vertical bending of the third order (f3 = 3.77 Hz). (d) Vertical bending of the fourth order (f4 = 4.63 Hz). (e) Vertical bending of the fifth order (f5 = 8.83 Hz).
(a)(b)(c)(d)(e)


## 5.5. Parameter Identification of Each Window

According to the sliding window technology and CDSI algorithm, the modal parameters of the cable-stayed bridge under two types of earthquake excitations are identified. Assume that the initial window time is 1 s, but the appropriate window size is given to be 4 s; additionally, the corresponding sliding window step size is 2 s, which means that window data is normally updated after every 2 s. The signals of each window are identified, and the stability results corresponding to each window in two types of working conditions are obtained. Figure10 shows the results of the first two stability diagrams in C1 working condition.Figure 10
Stability diagram (working condition: C1). (a) Stability diagram (time: 0 s–4 s). (b) Stability diagram (time: 2 s–6 s).
(a)(b)The specific values of the frequency in each window in operating conditions of C1 and S1 are identified, and the trend chart of the first five changing frequency values with time is given as shown in Figure11. Analyzing the frequency changes of each order in the figure, it is shown that the frequency value of the first five orders of bridge structure is stable in both conditions; that is, the bridge structure’s natural frequency is stable under the action of two types of earthquake waves. It is also shown that the cable-stayed bridge shows no damage in first two types of conditions.Figure 11
Frequency trend chart. (a) Frequency trend chart (working condition: C1). (b) Frequency trend chart (working condition: S1).
(a)(b)


## 5.6. Intelligent Identification of Modal Effectiveness

According to the proposed FCM clustering algorithm inSection 3.2, a clustering analysis of each window’s parameter results which is under two types of conditions is performed to reach intelligent identification of modal effectiveness in the stability diagram. Figure 12 is the corresponding effective stability diagram of C1 and S1 working conditions.Figure 12
Effective stability diagram. (a) Effective stability diagram (working condition: C1). (b) Effective stability diagram (working condition: S1).
(a)(b)Comparing the MIDAS finite element results with the first five orders of frequency results which are under the conditions of C1 and S1, as shown in Table3, the percent difference between the two results in the table and the errors between the effectiveness modal frequency values obtained by the SC-CDSI algorithm and the finite element results difference are within 7%; thus, frequency identification results of the proposed SC-CDSI algorithm are hereby considered reliable.Table 3
Comparative analysis results of frequency.
ResultsThe first orderThe second orderThe third orderThe fourth orderThe fifth orderMIDAS results (Hz)1.042.563.774.548.43Working condition: C1 (Hz)1.112.443.874.688.74Percentage difference 1 (%)−6.34.9−2.6−3.0−3.5Working condition: S1 (Hz)1.092.423.894.728.84Percentage difference 2 (%)−4.65.8−3.1−3.8−4.6


## 5.7. Modal Shape Results of First Three Orders

Based on these approaches, the SC-CDSI algorithm not only has verified accurately the true modes in the stability diagram but also has identified the modal shapes of structure accurately. Additionally, the first three mode shapes of cable-stayed bridge are obtained, as shown in Figure13, the X-coordinate in the figure represents the 11 measurement points of sensor, and the Y-coordinate is the mode shape after 0-1 planning; the dotted line is the identification result of the proposed algorithm; the solid line is the mode shape diagram, which is obtained by MIDAS. By comparing these shapes to the first three mode shapes selected from Figure 9, it can be seen that the modal modes and finite element results obtained by the SC-CDSI algorithm have more than 90% similarity, which is why the modal shapes obtained by SC-CDSI algorithm are reliable.Figure 13
Modal shape of the first three orders. (a) Vibration diagram of the first order. (b) Vibration diagram of the second order. (c) Vibration diagram of the third order.
(a)(b)(c)


## 6. Conclusion

Based on the limitations of deterministic-stochastic subspace algorithm, which is not suitable for nonlinear time-varying structures like bridge, a sliding window technique is proposed in combination with CDSI algorithm to implement the real-time tracking and identification of modal parameters for bridge structures. Additionally, considering that the identification of modal effectiveness in stability diagram needs manual participation, a fuzzy C-means (FCM) clustering algorithm is proposed to identify true and false modes in stability diagram to implement intelligent screening of modal parameter results. Finally, based on the shaking table test of a large cable-stayed bridge in lab, the proposed SC-CDSI algorithm is applied to identify the modal parameters of bridge, and its results are compared with the MIDAS finite element results. The compared results show that the proposed SC-CDSI algorithm not only accurately identifies the natural frequencies values of bridge structure but also accurately identifies the modal shape results and can accurately realize the online intelligent tracking and identification of modal parameters for real-life bridge.



---

*Source: 8855162-2021-03-11.xml*