System identification of linear structures based on Hilbert–Huang spectral analysis. Part 1: normal modes

Based on the Hilbert–Huang spectral analysis, a method is proposed to identify multi‐degree‐of‐freedom (MDOF) linear systems using measured free vibration time histories. For MDOF systems, the normal modes have been assumed to exist. In this method, the measured response data, which are polluted by noises, are first decomposed into modal responses using the empirical mode decomposition (EMD) approach with intermittency criteria. Then, the Hilbert transform is applied to each modal response to obtain the instantaneous amplitude and phase angle time histories. A linear least‐square fit procedure is proposed to identify the natural frequency and damping ratio from the instantaneous amplitude and phase angle for each modal response. Based on a single measurement of the free vibration time history at one appropriate location, natural frequencies and damping ratios can be identified. When the responses at all degrees of freedom are measured, the mode shapes and the physical mass, damping and stiffness matrices of the structure can be determined. The applications of the proposed method are illustrated using three linear systems with different dynamic characteristics. Numerical simulation results demonstrate that the proposed system identification method yields quite accurate results, and it offers a new and effective tool for the system identification of linear structures in which normal modes exist. Copyright © 2003 John Wiley & Sons, Ltd.     

Response of non-classically damped structures in the modal subspace

The evaluation of the dynamic response of non-classically damped linear structures requires the solution of an eigenproblem with complex eigenvalues and modal shapes. Since in practice only a small number of complex modes are needed, the complex eigenvalue problem is solved in the modal subspace in which the generalized damping matrix is not uncoupled by classical real modes. It follows that the evaluation of the structural response requires in both cases the determination of complex modes by numerical techniques, which are not as robust as techniques currently used for the solution of the real eigenvalue problem, and the use of complex algebra. In the present paper an unconditionally stable step-by-step procedure is presented for the response of non-classically damped structures in the modal subspace without using complex quantities. The method is based on the evaluation of the fundamental operator in approximated form of the numerical procedure. In addition, the method can be easily modified to incorporate the modal superposition pseudo-static correction terms.     

Modal decomposition method for stationary response of non-classically damped systems

The stationary response of multi-degree-of-freedom non-classically damped linear systems subjected to stationary input excitation is studied. A modal decomposition procedure based on the complex eigenvectors and eigenvalues of the system is used to derive general expressions for the spectral moments of response. These expressions are in terms of cross-modal spectral moments and explicitly account for the correlation between modal responses; thus, they are applicable to structures characterized with significant non-classical damping as well as structures with closely spaced frequencies. Closed form solutions are presented for the important case of response to white-noise input. Various quantities of response of general engineering interest can be obtained in terms of these spectral moments. These include mean zero-crossing rate and mean, variance and distribution of peak response over a specified duration. Examples point out several instances where non-classical damping effects become significant and illustrate the marked improvement of the results of this study over conventional analysis based on classical damping approximations.     

A subspace modal superposition method for non-classically damped systems

For structures with non-proportional damping, complex eigenvectors or mode shapes must be used in order to decoe the equations of motion. The resulting equations can then be solved in a systematic way. The necessity of solvie complex eigenvalue problem of a large system remains an obstacle for the practical application of the method. This stres utilizes the fact that in practice only a small number of the complex modes are needed. Therefore, these complex modes be approximated by a linear combination of a small number of the undamped modes, which can be obtained by established methods with less cost. An additional eigenvalue problem is then solved in a subspace with a much sm dimension to provide the best combination coefficient for each complex mode. The method of solution for the decoue equations is then carried over, using the approximate complex modes expressed in undamped mode shapes, to resue simple formulas for the time- and frequency-domain solution. Thus, an efficient modal superposition method is develoe for non-proportionally damped systems. The accuracy of this approximate method is studied through an example. Comparing the frequency response result using the approximate method with that using the exact complex modes, found that the error is negligible.     

Hilbert-Huang变换在密频结构阻尼识别中的应用

Hilbert—Huang变换是一种新的数据处理方法，由经验模分解(Empirical Mode Decomposition)技术及Hilbert变换两部分组成。本文研究此方法对于密频结构阻尼识别的应用。首先对于两自由度系统模型，说明该方法用于阻尼识别的步骤。进而研究存在频率密集现象的高层建筑的阻尼识别问题。上述结果与理论值及由半功率带宽法的识别值进行了比较，对比显示Hilbert．Huang方法较传统方法具有良好的识别密频结构阻尼的性能，适用于大型结构的系统识别。     

Simplifications of CQC method and CCQC method

The response-spectrum mode superposition method is widely used for seismic response analyses of linear systems. In using this method, the complete quadratic combination （CQC） is adopted for classically damped linear systems and the complex complete quadratic combination （CCQC） formula is adopted for non-classically damped linear systems. However, in both cases, the calculation of seismic response analyses is very time consuming. In this paper, the variation of the modal correlation coefficients of displacement, velocity and displacement-velocity with frequency and damping ratios of two modes of interest are studied, Moreover, the calculation errors generated by using CQC and square-root-of-the-sum-of-thesquares （SRSS） methods （or CCQC and CSRSS methods） for different damping combinations are compared. In these analyses, some boundary lines for classically and non-classically damped systems are plotted to distinguish the allowed minimum frequency ratio at given geometric mean of the damping ratios of both modes if their relativity is neglected. Furthermore, the simplified method, which is a special mode quadratic combination method considering only relativity of adjacent modes in CQC method and named simplified CQC or partial quadratic combination （PQC） method for classically damped linear system, is proposed to improve computational efficiency, and the criterion for determination of how many correlated modes should be adopted is proposed. Similarly, the simplified CCQC or complex partial quadratic combination （CPQC） method for the non-classically damped linear system and the corresponding criterion are also deduced. Finally, a numerical example is given to illustrate the applicability, computational accuracy and efficiency of the PQC and CPQC methods.     

Dynamic characteristics and seismic response of adjacent buildings linked by discrete dampers

Coupling adjacent buildings using discrete viscoelastic dampers for control of response to low and moderate seismic events is investigated in this paper. The complex modal superposition method is first used to determine dynamic characteristics, mainly modal damping ratio and modal frequency, of damper-linked linear adjacent buildings for practical use. Random seismic response of linear adjacent buildings linked by dampers is then determined by a combination of the complex modal superposition method and the pseudo-excitation method. This combined method can effectively and accurately determine random seismic response of non-classically damped systems in the frequency domain. Parametric studies are finally performed to identify optimal parameters of viscoelastic dampers for achieving the maximum modal damping ratio or the maximum response reduction of adjacent buildings. It is demonstrated that using discrete viscoelastic dampers of proper parameters to link adjacent buildings can reduce random seismic responses significantly. Copyright © 1999 John Wiley & Sons Ltd.     

A method for tuning tuned mass dampers for seismic applications

This paper studies tuned mass dampers (TMDs) resulting in high modal damping for mechanical systems incorporating such devices for the purpose of seismic response reduction. Focusing on the determination of damping and tuning, the proposed methodology identifies a point of multiplicity of complex eigenvalues and eigenvectors, resulting in different parameters for TMDs according to their location with respect to such multiplicity condition. It is shown that significant equal modal damping and average modal damping can be induced by properly tuning highly damped TMDs, obtaining parameters intrinsic to the mechanical systems, and excitation independent. Further, it is shown that the methodology yields, as particular cases, two proposals by others using TMDs for the same purpose of seismic response abatement. Copyright © 2012 John Wiley & Sons, Ltd.     

Complex complete quadratic combination method for damped system with repeated eigenvalues

A new response-spectrum mode superposition method, entirely in real value form, is developed to analyze the maximum structural response under earthquake ground motion for generally damped linear systems with repeated eigenvalues and defective eigenvectors. This algorithm has clear physical concepts and is similar to the complex complete quadratic combination (CCQC) method previously established. Since it can consider the effect of repeated eigenvalues, it is called the CCQC-R method, in which the correlation coefficients of high-order modal responses are enclosed in addition to the correlation coefficients in the normal CCQC method. As a result, the formulas for calculating the correlation coefficients of high-order modal responses are deduced in this study, including displacement, velocity and velocity-displacement correlation coefficients. Furthermore, the relationship between high-order displacement and velocity covariance is derived to make the CCQC-R algorithm only relevant to the high-order displacement response spectrum. Finally, a practical step-by-step integration procedure for calculating high-order displacement response spectrum is obtained by changing the earthquake ground motion input, which is evaluated by comparing it to the theory solution under the sine-wave input. The method derived here is suitable for generally linear systems with classical or non-classical damping.     

Eigenproperties of non-classically damped primary structure and equipment systems by a perturbation approach

Closed-form expressions are obtained to calculate the approximate complex eigenvalues and eigenvectors of a system composed of a non-classically damped primary structure and a single degree of freedom oscillator. The expressions are obtained through a systematic second order perturbation analysis of a transformed eigenvalue problem of the combined system. The possibility of tuning between the structure and equipment is considered. The dynamic properties of the combined system are derived in terms of the complex eigenvalues and eigenvectors of the supporting structure and the frequency, mass and damping ratio of the equipment. Examples demonstrating the accuracy of the expressions for the eigenvalues and eigenvectors are presented. These eigenproperties are used for generation of floor response spectra for non-classically damped structures to incorporate the dynamic interaction effects between the structure and equipment.     

Output only modal identification and structural damage detection using time frequency ＆ wavelet techniques

The primary objective of this paper is to develop output only modal identifi cation and structural damage detection.Identif ication of multi-degree of freedom(MDOF) linear time invariant(LTI) and linear time variant(LTV—due to damage) systems based on Time-frequency(TF) techniques—such as short-time Fourier transform(STFT),empirical mode decomposition(EMD),and wavelets—is proposed.STFT,EMD,and wavelet methods developed to date are reviewed in detail.In addition a Hilbert transform(HT) approach to determine ...     

Free‐vibration analysis of a three‐dimensional soil–structure system

A procedure which involves a non‐linear eigenvalue problem and is based on the substructure method is proposed for the free‐vibration analysis of a soil–structure system. In this procedure, the structure is modelled by the standard finite element method, while the unbounded soil is modelled by the scaled boundary finite element method. The fundamental frequency, and the corresponding radiation damping ratio as well as the modal shape are obtained by using inverse iteration. The free vibration of a dam–foundation system, a hemispherical cavity and a hemispherical deposit are analysed in detail. The numerical results are compared with available results and are also verified by the Fourier transform of the impulsive response calculated in the time domain by the three‐dimensional soil–structure–wave interaction analysis procedure proposed in our previous paper. The fundamental frequency obtained by the present procedure is very close to that obtained by Touhei and Ohmachi, but the damping ratio and the imaginary part of modal shape are significantly different due to the different definition of damping ratio. This study shows that although the classical mode‐superposition method is not applicable to a soil–structure system due to the frequency dependence of the radiation damping, it is still of interest in earthquake engineering to evaluate the fundamental frequency and the corresponding radiation damping ratio of the soil–structure system. Copyright © 2001 John Wiley & Sons, Ltd.     

Optimal design of viscoelastic dampers using eigenvalue assignment

In this study a procedure for determining the optimum size and location of viscoelastic dampers is proposed using the eigenvalue assignment technique. Natural frequencies and modal damping ratios, required to realize a given target response, are determined first by the convex model. Then the desired dynamic structural properties are realized by optimally distributing the damping and stiffness coefficients of viscoelastic dampers using non‐linear programming based on the gradient of eigenvalues. This optimization method provides information on the optimal location as well as the magnitude of the damper parameters. The proposed procedure is applied to the retrofit of a 10‐story shear frame, and to a three‐dimensional structure with an asymmetric plan. The analysis results confirm that the responses of model structures retrofitted by the proposed method correspond well with the given target response. Copyright © 2003 John Wiley & Sons, Ltd.     

两种识别结构模态参数时域法的对比研究

模态参数是有效评估结构安全状况的关键参数,在结构抗震加固和健康诊断领域得到广泛应用。与频域法相比较,时域法直接利用实测的振动信号识别模态参数,不需要进行频域变换,减少数据处理带来的误差,并且可以实现大型结构的在线识别,真实地反应结构的现状。以同济大学12层钢筋混凝土标准框架振动台模型试验完整数据为对象,在详细介绍ITD法和复指数法2种时域法理论的基础上,通过编程选取结构不同测点的振动加速度时程数据,识别了小震和强震工况下12层钢筋混凝土框架模型振动台试验模型的模态频率和阻尼比,并结合移动谱识别结构模态参数的时变特性。结果表明:ITD法和复指数法可有效地识别结构的模态参数,自振频率的识别精度较高,而阻尼比的离散度较大;小震工况频率变化值不大,而强震工况频率值较初始时刻有明显的下降,这与试验现象是吻合的,进一步说明移动谱与这2种时域法相结合可以反应结构在塑性阶段的参数时变特性。     

基于加窗Hilbert变换的复偏振分析方法及其应用

考虑到传统Hilbert变换所隐含的无限长序列假设的实际局限性, 本文将多窗分析法引入到短时序列Hilbert变换中, 通过构建复协方差矩阵和对该矩阵的特征值与特征向量的求解, 获得三分量地震记录的时变偏振参数;基于不同波型的偏振特性与实测偏振参数, 采用线性_余弦权重函数或高斯权重函数自适应空间滤波, 识别与分离具有特定偏振特性的不同地震响应. 针对天然地震三分量记录的处理结果表明, 该方法在识别、分离具有特定偏振特性的不同地震响应方面具有一定的潜力.     

Seismic response analysis of highway overcrossings including soil–structure interaction

This paper presents a systematic procedure for the seismic response analysis of highway overcrossings. The study employs an elementary stick model and a more sophisticated finite element formulation to compute response quantities. All dynamic stiffnesses of approach embankments and pile groups are approximated with frequency‐independent springs and dashpots that have been established elsewhere. A real eigenvalue analysis confirms the one‐to‐one correspondence between modal characteristics obtained with the three‐dimensional finite element solutions and the result of the simpler stick‐model idealization. A complex eigenvalue analysis yields modal damping values in the first six modes of interest and shows that modal damping ratios assume values much higher than those used by Caltrans. The validity of the proposed method is examined by comparing the computed time response quantities with records from the Meloland Road and Painter Street overcrossings located in southern and northern California, respectively. The proposed procedure allows for inexpensive parametric analysis that examines the importance of considering soil–structure interaction at the end abutments and centre bent. Results and recommendations presented by past investigations are revisited and integrated in comprehensive tables that improve our understanding of the dynamic characteristics and behaviour of freeway overcrossings. The study concludes with a step‐by‐step methodology that allows for a simple, yet dependable dynamic analysis of freeway overcrossings, that involves a stick model and frequency‐independent springs and dashpots. Copyright © 2002 John Wiley & Sons, Ltd.     

Modal equations of linear structures with viscoelastic dampers

Several types of energy dissipation devices using viscoelastic materials have been proposed to reduce vibration in structures subjected to wind and earthquake excitations. At constant temperature and small strain levels, the mechanical behaviour of Viscoelastic (VE) materials can be described using linear operators. In general, the stiffness and damping matrices of structures using VE devices are frequency dependent; this implies that the classical second-order differential equations for the modal co-ordinates are not a complete model for this type of structures. In this paper, the concept of modal coupling in the frequency domain is addressed, expressions for diagonalizable frequency-dependent stiffness and damping matrices are given, and an iterative technique for the computation of the response of viscoelastic structures is studied. Necessary and sufficient conditions for convergence of the technique are given and numerical examples are developed to illustrate the application of the method.     

基于复阻尼理论的混合结构Rayleigh阻尼模型

混合结构的阻尼矩阵不满足经典阻尼条件,导致传统的模态叠加法无法适用。复阻尼理论无法适用于时域计算,其自由振动响应中存在发散现象。针对混合结构的阻尼矩阵非比例性和复阻尼理论的时域发散性,基于频域等效原则构建了求解Rayleigh阻尼系数的数学优化模型,进而得到与复阻尼理论等效的Rayleigh阻尼运动方程。算例分析表明:依据位移时程响应和结构等效阻尼比可证明Rayleigh阻尼运动方程的正确性。基于本文研究成果,等效复阻尼理论的混合结构Rayleigh阻尼运动方程可直接采用模态叠加法,结合其确定的结构等效阻尼比,为混合结构的振型分解反应谱法提供理论依据。     

基于小波变换的层间隔震结构地震响应分析

利用小波多分辨率分析将地震动加速度分解为多频段小波分量,并运用复模态方法推导其计算层间隔震体系在地震作用下的动力响应公式,讨论各频段地震信号及结构响应的能量分配。同时利用小波时频工具分析地震动能量在时频域内的分布对层间隔震结构响应的影响,进而为考察地震动非平稳性对层间隔震结构非线性分析的影响提供方法。利用小波分析的以上优势,对一典型层间隔震结构分别进行弹性和弹塑性分析,结果表明弹性体系在地震作用下的响应可由该地震波各小波分量的响应叠加而得,地震动能量在时间上的集中会对层间隔震结构响应产生不利影响。     

Towards a seismic design method for plane steel frames using equivalent modal damping ratios

A method for seismic design of plane steel moment resisting frames based on the use of equivalent modal damping ratios is developed. The method determines the design base shear of the structure through spectrum analysis using rationally obtained equivalent modal damping ratios instead of the crude strength reduction (behavior) factor. An equivalent linear structure, which retains the mass and initial stiffness of the original non-linear structure and takes into account geometrical non-linearity and inelasticity in the form of equivalent, time-invariant, modal damping ratios is established. The equivalent damping ratios for the first few significant modes are numerically computed by first iteratively forming a frequency response transfer function modulus until it satisfies certain smoothness criteria and then by solving a set of non-linear algebraic equations. Thus, design equations providing equivalent damping ratios as functions of period and allowable deformation and damage are constructed using extensive numerical data coming from plane steel moment resisting frames excited by various seismic motions. These equations can be used in conjunction with a design spectrum, appropriately constructed for high damping values, and modal synthesis tools to calculate the seismic design forces of the structure. The proposed method is illustrated by numerical examples. It is concluded that unlike the usual approach of seismic codes employing a single common value of the strength reduction factor value for all modes, the proposed approach working with deformation and damage dependent equivalent modal damping ratios leads to more accurate and deformation and damage controlled results.