Abstract

The damping is important for forward and inverse structural dynamic analysis, and damping identification has become a hot issue in structural health monitoring recently. The dynamic responses of the structure can be measured in practice, and the structural parameter usually can be identified by inverse response sensitivity analysis. To reduce the measurement noise effect and enhance the effectiveness of the response sensitivity method, an enhanced sensitivity analysis method was proposed to identify the structural damping based on the Principal Component Analysis (PCA) method. The measured acceleration responses were analyzed by PCA method, and the updated analytical responses and the response sensitivities were projected into the subspace determined by the first-order principal component. The projection equations were adopted to identify the parameters of damping model. The proposed damping identification method was numerically validated with a planar truss structure at first, and then the experimental study was conducted with a steel planar frame structure. It shows that the proposed method is effective in identifying the parameters of damping model with better accuracy compared with the conventional acceleration response sensitivity method, and it is also robust to the sensor placement and measurement noise.

1. Introduction

The damping ratio is a dimensionless measure and a measure of describing how rapidly the oscillations of a structural system decay from one bounce to the next, which is a significant factor when analyzing the structural dynamic behavior dominated by energy dissipation [1]. Unlike the mass or stiffness that can be measured or determined by static test, the damping cannot be determined by measurement or static test [2], but the damping characteristic is very important in structural health monitoring. Accurate damping matrix construction is a determining factor in analyzing structural dynamic responses and predicting energy dissipation behavior, which makes the damping estimation become a key issue for the structural design, dynamic response analysis, structural health monitoring, etc.

In order to construct accurate structural damping matrix, damping identification has been studied by some researchers for different systems, including rotor systems [3], mistuned blisks [4] and monopile foundation [5]. For the structural system, the damping identification method can be classified into frequency domain method [68], time domain method [911], and time-frequency domain method [12, 13]. The half-power bandwidth method [6, 7] in frequency domain is usually adopted in the dynamic test, and Wang [8] studied the errors of calculating damping effect between the classical and the third-order half-power method. Li and Law [9] proposed a time domain damping ratio identification method based on the acceleration response sensitivity, and the proposed method is validated by numerical study and experimental study [11]. The time domain response sensitivity method was combined with iterative regularization method to identify damping ratios [10]. The wavelet analysis method [12, 13] and Hilbert–Huang Transform method [14] have also been adopted to identify the damping ratio.

Sensitivity analysis can be used to estimate the system output variation due to a perturbation in the system parameters by means of partial derivatives [15, 16], and forward sensitivity analysis has been used in many applications [1719]. Inverse sensitivity-based method with model updating is usually based on a first-order Taylor series that minimizes an error function to assess the system parameter perturbation, and time domain sensitivity method has been adopted in structural parameter identification widely [2022]. The time domain response sensitivity method has advantages including no requirement of computing the higher order system model parameters, obtaining the responses easily, and providing more identification equations, which makes it a good tool for damping identification. Time domain response sensitivity method has also been studied and gained significant attention in damping identification [911].

The time domain response sensitivity method also has the disadvantage of being sensitive to the measurement noise [23, 24], and it is also important to enhance the response sensitivity for the structural parameter identification [2426]. The Principal Component Analysis (PCA) technique, also known as Karhunen-Loeve transform or proper orthogonal decomposition [27], decomposes data series through orthogonal linear transformation to get the principal components, the first few of which contain more information of the parameter variation and less random noise information. The relationship between the system parameters and its output has been broadly studied with a combination of sensitivity analysis and Principal Component Analysis (PCA) [28, 29], and it has been proved that PCA method can improve the response sensitivity for structural damage identification [25, 26] with subspace projection method.

With the rapid development of measurement technique, the time domain responses of the structure can be obtained, and the acceleration responses of the structure can be measured easily. This paper will propose a damping identification method based on the acceleration response measurement. The inverse acceleration response sensitivity method for damping ratio identification is revisited, and the model updating method is also briefly reviewed in this paper at first. The measured acceleration responses are decomposed by PCA method, and time domain response sensitivity equations of damping ratio identification are projected into the finer subspace. With the iterative procedure of model updating and subspace projection, the enhanced sensitivity method is proposed to identify the parameters of damping model. The proposed method is described in detail, and it is validated with simulation studies on a plane truss structure, in which different sensor placements and different measurement noise levels are studied. A seven-storey steel frame was designed and manufactured in the laboratory, and hammer test was performed. The proposed acceleration response sensitivity enhancement method is used to identify the damping ratio of the steel frame structure.

2. Methodology

2.1. Damping Identification Based on Acceleration Response Sensitivity

The equation of motion of a damped linear structure can be written aswhere , , and are the mass, damping, and stiffness matrices of the structural system, respectively. is the vector of excitations on the structure and is the mapping matrix for the excitations. , , and are vectors of the acceleration, velocity, and displacement responses, respectively.

Performing differentiation to both sides of (1) with respect to damping ratio , which is the critical parameter in structural damping model, we havewhere is the th damping ratio of the structural system, and , , and are the acceleration, velocity, and displacement sensitivity vectors, respectively, which can be determined by Newmark- method solving (2).

The acceleration sensitivity vector corresponding to th damping ratio can be rewritten as . All the sensitivity vectors are assembled asThe identification equation for all the damping ratios of a structure can be represented asThe higher order term can be omitted in (4). With an iterative method the damping ratio perturbations can be determined from (4), and Tikhonov regularization is used for optimizing the following objective function in the th iteration aswhere is the regularization parameter in the th iteration obtained with the L-curve method [30].

After is solved, the damping matrix is updated with where is the assumed initial damping ratio vector.

Then after recalculating the structural responses and the sensitivity matrix, the vector for the next identification iteration is obtained until the given convergence criterion is met:Tol values are selected to meet the challenge in convergence of the identified results with measurement noise effect, and then the final value of the damping ratio can be obtained.

2.2. Acceleration Response Sensitivity Enhancement by PCA

The initial structural damping ratio vector is assumed as , and the initial analytical acceleration responses corresponding to in the initial stage can be determined from (1) with the Newmark-β method. The measured acceleration responses can be obtained with the accelerometers in the structure. Based on the initial analytical acceleration response vectors and measured acceleration response vectors , the identification equation for the first identification iteration can be written asThe PCA is applied to the measured acceleration responses , and the covariance matrix of is , when written in the form of spectral decomposition as where is the vector of eigenvalues of the covariance matrix; is the corresponding eigenvector matrix which defines an orthogonal subspace; p is the number of measured acceleration responses here.

The principal components of measured acceleration responses are then obtained as The initial analytical responses can be projected into the subspace constructed by the eigenvector matrix asThe same projection can be conducted to the sensitivity vectors, so the projection of (8) with the same subspace can be written asThe orthogonal vectors insure each equation in (12) is independent in the subspace. The dimension of the projection sensitivity matrix for each principle component is , where is the length of the measured acceleration and is the number of the damping ratios. If there are principal components selected for determining the damping ratios, the projection sensitivity matrix dimension becomes .

The principal components, which contain less measured noise information and most of the structural damping perturbation information, will be selected, and the selection of principal components will be discussed in next section. Then the damping ratio perturbation can be obtained by Tikhonov regularization method as (5) with the projection sensitivity equation in (12).

Similar to the conventional response sensitivity method [9], the model updating procedure is conducted when each damping ratio perturbation is determined in the th iterative step. After each iteration step, the damping ratio is updated with (6), and the analytical acceleration and sensitivity matrix are updated. The procedure of the proposed method for damping ratios identification is presented as the following 5 steps.

Step 1. The measured acceleration responses are analyzed with the PCA method, and the principal components, which contain less measured noise information and more damping ratio perturbation information, are selected.

Step 2. Based on the finite element model and initial damping ratio, the corresponding analytical acceleration responses and sensitivity matrix are calculated with the Newmark- method, and then the identification equation in (12) is determined based on (1)-(12).

Step 3. The identification equation is solved with Tikhonov regularization method to obtain the damping ratio perturbations.

Step 4. The damping ratio vector is updated, and with the new damping ratios, the damping matrix of the structural model is updated and the analytical responses and corresponding response sensitivity matrix are recalculated. The identification equation in (12) is determined based on (1)-(12) with the updated damping ratios.

Step 5. Step 3 to Step 4 are repeated, and the iterative process will stop until the convergence criterion defined in (6) is satisfied.

The flowchart of the proposed damping identification method is shown in Figure 1.

3. Numerical Simulation Study

3.1. Numerical Simulation Case

A two-dimensional truss structure as shown in Figure 2 serves for the numerical study. This structure is modeled using 61 truss finite elements, and each node consists of 2 Dofs including the translational displacements in horizontal and vertical directions. This structure is pin-supported at Node 1 and roller-supported at Node 25, and the node number and element number are also shown in Figure 2. Each vertical and horizontal member is of 1.0-meter length, and the cross-sectional area of all members is 0.0016 m2. The Rayleigh damping model is adopted, which has two damping ratios needed to be identified.

The mass density of material and the elastic modulus of material are 7.8×103 kg/m3 and 2.06 GPa, respectively. The system mass and stiffness matrices can be established based on the geometry and material properties of the model. The first eight natural frequencies of the structure are 1.5792 Hz, 5.2627 Hz, 7.0661 Hz, 12.6623 Hz, 18.4986 Hz, 21.2952 Hz, 29.0158 Hz, and 32.3490 Hz, respectively.

It assumes that the real damping ratio coefficients associated with the first and second modes are 1% and 2% in the Rayleigh damping model, respectively. The truss structure is subjected to vertical and horizontal external dynamic forces as shown in (13), and the dynamic forces are modeled as follows to simulate excitations over a relatively wide range of frequencies.The accelerometers are set to collect acceleration responses of the structure, and the duration of measurement is 1.0 second with 1000 Hz sampling rate. The “measured responses” are simulated by adding random components to the “measured responses without noise effect” aswhere is the “measured response without noise effect” which can be determined by Newmark- method with the real damping ratio in numerical simulation; is the noise level; is the standard deviation of the “measured response without noise effect”; and is a vector of random values with zero mean and unit standard deviation.

It is impractical to measure the responses of all the Dofs in a real structure, and the measurement location will affect the identification result. This work will not discuss the sensor placement effect, and six arbitrarily selected sets of the sensors, which are shown in Table 1, are adopted to validate that the proposed damping identification method is robust to sensor placement.

3.2. Damping Identification without Measurement Noise

Both of the initial damping ratio coefficients are assumed to be 3%. The structural responses are simulated with the Newmark- method from the structure with the real damping ratios and initial damping ratios, in which the responses with the real damping ratios are considered as “measured responses without noise effect” and the responses with the initial damping ratios are considered as “initial analytical responses”. The measurement noise is not considered in this case, so the “measured responses without noise effect” are adopted to identify the damping ratio in this section to validate the correctness of the proposed method.

The “measured responses without noise effect” from different sensor placement are decomposed by PCA method, and the corresponding “analytical responses” are projected into the subspace constructed by the principal components as (11). The damping ratios are identified with single principal component based on the “measured responses without noise effect”, and the identified results from the proposed method with different sensor placements are listed in Table 2. It shows that the identified results match the real value perfectly indicating the correctness and good accuracy of the proposed method in damping identification, and it can be concluded that without considering the measurement noise the satisfactory identification results can be achieved based on every single component.

The PCA decompose time domain data series into principal components via orthogonal linear transformation, and the components corresponding to larger eigenvalues usually contain more information of the original data series. In order to study the perturbation information of damping ratio in the principal components of the measured responses, a ratio is defined in (15) to present the relative perturbation information of damping ratio in different principal components.where is the norm of the difference between the jth principal component of the “measured responses without noise effect” and the corresponding projected component of “initial analytical responses”, and is the norm of the difference between the “measured responses without noise effect” and “initial analytical responses”.

All the ratios with different sensor placements are shown in Table 3. The ratios of the lower order principal component are usually larger than that of higher order principal components. It indicates that the lower order principal component contains more perturbation information of damping ratios, which makes the lower order principal component better for damping ratio identification than higher order components.

3.3. Principal Components Selection for Damping Identification with Measurement Noise

In practice the measurement noise cannot be avoided, so the measurement noise effect should be considered in damping identification. All the principal components of the “measured responses without noise effect” can be used to identify the damping ratio with perfect accuracy, but the principal component will be affected by the measurement noise. In this section, the measurement noise information distribution in the principal component will be discussed, and the principal component selection criterion will be given.

Both of the initial damping ratio coefficients are also assumed to be 3% in this section. The structural responses are simulated with the Newmark-β method from the structure with the real damping ratios and initial damping ratios, in which random noise is added to the responses with the real damping ratios being considered as “measured responses”. In this section 10% measurement noise is considered in damping ratio identification and the measurement noise information distribution in different principal components is discussed.

The “measured responses” from different sensor placement are decomposed by PCA method, and the corresponding “measured responses without noise effect” are projected into the subspace constructed by the principal components as (11). In order to study the measurement noise information in the principal components of the measured responses, a ratio is defined in (16) to present the measurement noise information distribution in different principal components.where is the norm of the difference between the jth principal component of the “measured responses” and the corresponding projected component of “measured responses without noise effect”.

Ten different measurement noise series are added to the “measured responses without noise effect” to simulate the “measured responses” as in (14), and the ratio defined in (16) is calculated. The values of with different measurement noise series are shown in Figure 3; it seems that the noise information distribution in every principal component is stable for all the sensor placements, in which the first principal component always contains least noise information. The mean values for different sensor placements are shown in Table 4, and it can be seen that the lower order principal components are noted to be less affected by measurement noise.

In the last section, it has been proved that the lower order principal components of the “measured responses without noise effect” also usually contain more perturbation information of damping ratios. The “measured responses without noise effect” in (15) are replaced by “measured responses”, and the values of with different measurement noise series are shown in Figure 4. It can be seen that the perturbation information of damping ratio in every principal component is stable for all the sensor placements with measurement noise, and the damping perturbation information distribution is similar to that without measurement noise. The mean values for different sensor placements are shown in Table 5, and they are close to the ratio of the “measured responses without noise effect”.

Based on the above discussions the lower principal components contain more perturbation information of damping ratios and less measurement noise information, so the identification equation based on the lower order principal component will give more accurate identification result considering the measurement noise. The first-order principal component is adopted to identify the damping ratios, and ten different measurement noise series are considered for each sensor placement. The ten identification results are shown in Figure 5 with the proposed method and conventional sensitivity method, and it can be seen that the proposed method can give more accurate and stable identification result than the conventional method.

The mean values and coefficients of variation of the identification results with the proposed method are shown in Table 6, and those of the identification results with the conventional response sensitivity are shown in Table 7. Comparing the mean values in Table 6 with those in Table 7, the mean values of the identification results from proposed method are closer to the true value, so it can be concluded that the proposed method can give more accurate identification results.

Comparing the coefficients of variation in Table 6 with those in Table 7, the coefficients of variation, which can represent the standardized measure of dispersion, for identification results from the proposed method are much smaller than those in Table 7, so it shows that the proposed method can give more stable identification results.

3.4. Damping Identification with Different Measurement Noise Level

In this section 20% and 30% measurement noise levels are considered to investigate the performance of the proposed method under higher level measurement noise, and the identification results are compared with the conventional response sensitivity method. The evolution of damping ratios in these two noise level scenarios with different sensor placements is shown in Figure 6, and the finial identified results are shown in Table 8. The identification results in Figure 6 and Table 8 with 20% and 30% noise levels indicate that the sensor placement effect is not significant for the proposed method, but for the conventional method the difference of the identification results between different sensor placements is much larger.

The damping ratio identification results show that the proposed method can identify the damping ratio with acceptable accuracy even under 30% measurement noise, and the difference of identification results between different noise levels is very small, which means the proposed method is less affected by measurement noise. With the increasing of the measurement noise level, the results from the conventional method become worse, and with 30% measurement noise the identification error of the conventional method is too large to be accepted.

4. Experimental Study

4.1. Experimental Model and Finite Element Modeling

A 7-storey planar steel frame is adopted to verify the proposed damping identification method in the laboratory, and more details of the test can be found in the thesis of the second author [26]. The configuration of the test model is shown in Figure 7. The measured mass densities of the column and beam materials are 7850kg/m3 and 7764kg/m3, respectively, and the measured cross sections of the column and beam elements are 50.06 mm × 5.10 mm and 49.88 mm × 8.06 mm, respectively. Additional mass blocks have been placed on the 1/4 and 3/4 length along the beam members to simulate the inertia of floor slab in practice, and the weights of mass blocks are listed in Table 9.

B&K 3023 accelerometers shown in Figure 8(a) and B&K Nexus amplifiers shown in Figure 8(b) were adopted to collect the acceleration responses. Impact test was performed with SINOCERA LC-04A dynamic hammer shown in Figure 8(c) to apply the impact excitation, and the National Instrument data acquisition board shown in Figure 8(d) was adopted to record data for the test.

To reduce the discrepancy between the analytical model and the experimental model, the initial finite element model is updated before the damping ratios identification. The baseline finite element model of the test structure is obtained from model updating with the modal sensitivity approach. Young’s modulus for each element and the stiffness values for two restraints at the supports are updated, because the dimensions and mass can be measured in situ. The detailed procedure for the baseline model updating can be found in [26], in which the accuracy of the baseline finite element model is verified by the minimized errors in the Modal Assurance Criteria (MAC) values before and after updating.

4.2. Damping Identification of the Experimental Model

Based on the Rayleigh damping model assumption, the proposed method is adopted to identify the first two damping ratios of the test model. The measured acceleration responses of the test model under hammer impact force are adopted, and the hammer impact location and measurement locations are shown in Table 10. Both of the initial damping ratios for the first and second modes are set as 0.5%, and the analytical responses calculated by the assumed damping ratios and baseline finite element model are compared with the measured responses shown in Figures 9(a), 9(c), and 9(e).

Because the true values of the damping ratios of the test model are unknown, the relative errors between the measured responses and analytical responses with the identified damping ratios are selected to represent the identification error. The damping ratios are identified as and with the proposed method, and the analytical responses calculated by the identified damping ratios and baseline finite element model are compared with the measured responses in Figures 9(b), 9(d), and 9(f).

The relative errors between the initial analytical responses and measured responses are shown in Table 11, and the corresponding relative errors between the analytical responses with the identified damping ratios and measured responses are also shown in Table 11. It shows that with the identified damping ratios the analytical responses match the measured responses very well and using the identified damping ratios from the proposed approach to calculate the analytical response can give a smaller relative error in the dynamic response prediction compared with the experimental measured response.

5. Conclusions and Discussions

A damping identification method is proposed based on the combination of response sensitivity analysis and PCA method, in which the response sensitivity equation is projected into the subspace constructed by the first principal component of the measured responses. The first principal component of the measured responses contains more damping ratio variation information and less noise information, so the measurement noise effect is reduced significantly and the sensitivity of the measured response is enhanced simultaneously by the projection. This proposed method is validated in numerical studies with a planar truss structure with up to 30% measurement noise, and the enhanced sensitivity method performs better compared to the conventional sensitivity-based method for damping ratio identification. With a steel place frame, the proposed method is also validated experimentally, and the identification results also provide experimental evidences to support conclusions drawn above.

It should be noted that the Rayleigh damping model is adopted in this paper, and the damping ratios in Rayleigh damping model are identified based on the proposed method. A generalized projection method was proposed, so we can use similar procedure to identify the damping parameters in other types of damping model. It should also be noted that the identification results of the proposed method have some variations with different sensor placements under high level measurement noise, so the sensor placement effect and improved damping identification method should be studied later.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; and in the decision to publish the results.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The work described in this paper was supported by “Financial Grant from the China Postdoctoral Science Foundation” (Grant no. 2015M580267; Grant no. 2016T90296) and “National Natural Science Foundation of China” (Grant no. 51708159; Grant no. 51008094). The Ph.D. thesis of the second author [26] provided theoretical foundation for this paper.