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.     

The comparative performance of seismic response spectrum combination rules in building analysis

The peak dynamic responses of two mathematical models of a fifteen-storey steel moment resisting frame building subjected to three earthquake excitations are computed by the response spectrum and time history methods. The models examined are: a ‘regular’ building in which the centres of stiffness and mass are coincident resulting in uncoupled modes with well-separated periods in each component direction of response; and an ‘irregular’ building with the mass offset from the stiffness centre of the building causing coupled modes with the translational modes having closely spaced periods. Four response spectrum modal combination rules are discussed and are used to predict the peak responses: (1) the square root of the sum of the squares (SRSS) method; (2) the double sum combination (DSC) method; (3) the complete quadratic combination (CQC) method; and (4) the absolute sum (ABS) method. The response spectrum results are compared to the corresponding peak time history values to evaluate the accuracy of the different combination rules. The DSC and the CQC methods provide good peak response estimates for both the regular and irregular building models. The SRSS method provides good peak response estimates for the regular building, but yields significant errors in the irregular building response estimates. The poor accuracy in the irregular building results is attributable to the effects of coupled modes with closely spaced periods. It is concluded that the DSC and CQC methods produce response estimates of equivalent accuracy. Both methods are recommended for general use. In addition to the DSC and CQC rules, the SRSS method is recommended for systems where coupled modes with closely spaced periods do not dominate the response.     

Conservatism in summation rules for closely spaced modes

It is shown that the method recommended by the Nuclear Regulatory Commission to be used to combine spectral response in the case of closely spaced modes is unnecessarily conservative for certain systems. Closely spaced modes arise in structures from symmetry and where there is a light appendage with a frequency close to one of the natural frequencies of the structure. In the former case, the closely spaced modes do not involve significant interaction between components of the system and the Nuclear Regulatory Commission Guide is reasonable. The latter case, that is when there are closely spaced modes where interaction of components occurs as in the examples of light appendages and torsionally unbalanced buildings, must be treated by consideration of the interacting components. The approach proposed here is that the modes that are not closely spaced be treated by modal analysis and the closely spaced modes, in the case of two closely spaced modes, be treated as a coupled two-degree-of-freedom system. If this is done, the beat phenomenon, the most important characteristic of the interaction, is evident, as is the associated result that the peak response of the coupled system is developed much later than the peak responses obtained in the individual modes. It is shown that the square root of the sum of the squares procedure underestimates, as expected, the response for undamped and very lightly damped systems, but for damped systems the square root of the sum of the squares method can be extremely conservative. It follows that the other methods specified by the Nuclear Regulatory Commission for closely spaced modes must be even more conservative.     

Dynamic characteristics of non-classically damped structures

The equations of motion of a structure in undamped modal coordinates may have non-zero off-diagonal terms in the damping matrix. Although these terms are commonly neglected, studies have shown that they may have a significant influence on the response to dynamic loads. In this paper, two independent criteria are developed to determine when these damping terms will affect the structure's modal properties and response. It is found that even small off-diagonal damping values can be significant if the structure has closely spaced natural frequencies. To quantify and understand the influence of these damping terms, closed-form analytical expressions are derived for the modal properties and harmonic and stochastic response of structures with closely spaced natural frequencies. One conclusion is that off-diagonal damping terms will decrease a modal damping ratio for each pair of closely spaced modes. This is significant, since a response analysis performed by neglecting these off-diagonal terms will underestimate the true response.     

有关Kanai-Tajimi模型的统计特征分析

选择金井清过滤白噪声模型,采用留数定理推导出K-T谱激励的自相关函数表达式,并在此基础上利用相关函数法分析出K-T谱激励下的平稳反应,并采用功率谱方法加以验证。这些结果可为结构随机地震反应时域分析和抗震可靠性评估提供基础。     

Probabilistic approach for modal identification using non‐stationary noisy response measurements only

This paper addresses the problem of identification of the modal parameters for a structural system using measured non‐stationary response time histories only. A Bayesian time‐domain approach is presented which is based on an approximation of the probability distribution of the response to a non‐stationary stochastic excitation. It allows one to obtain not only the most probable values of the updated modal parameters and stochastic excitation parameters but also their associated uncertainties using only one set of response data. It is found that the updated probability distribution can be well approximated by a Gaussian distribution centred at the most probable values of the parameters. Examples using simulated data are presented to illustrate the proposed method. Copyright © 2002 John Wiley & Sons, Ltd.     

Rosenblueth's modal combination rule for systems with non-classical damping

A simple rule is derived to combine, within the framework of a complex mode superposition, the maximum modal responses of systems such as soil-structure and structure-equipment systems, for which closely spaced natural frequencies are likely, and for which, because of the large difference in the damping values of their various components, the assumption of an orthogonal damping matrix may lead to significant errors. The rule constitutes the generalization of Rosenblueth's rule for systems with closely spaced natural frequencies and classical modes, and is expressed in terms of their complex mode shapes and natural frequencies. Its derivation is based on the theory of a complex modal analysis for systems with non-classical modes of vibration and on Rosenblueth's original derivation. As in this original derivation, earthquake ground motions are modelled as a stationary white noise process, but the formulae obtained under this assumption are modified later on to account for the transient nature of actual earthquakes. A numerical example is presented to illustrate the application of the rule, and a comparative study with numerical integration solutions is performed to assess its accuracy. In this comparative study, it predicts the numerical integration solutions with an average error of 0.3 per cent.     

The complete SRSS modal combination rule

The complete Square‐Root‐of‐Sum‐of‐Squares (c‐SRSS) modal combination rule is presented. It expresses the structural response in terms of uncoupled SDOF modal responses, yet accounting fully for modal response variances and cross‐covariances. Thus, it is an improvement over the classical SRSS rule which neglects contributions from modal cross‐covariances. In the c‐SRSS rule the spectral moments of the structural response are expressed rigorously in terms of the spectral moments of uncoupled modal responses and of some coefficients that can be computed straightforwardly as a function of modal frequencies and damping, without involving the computation of cross‐correlation coefficients between modal responses. An example shows an application of the c‐SRSS rule for structural systems with well separated and closely spaced modal frequencies, subjected to wide‐band and narrow‐band excitations. Comparisons with response calculations using the SRSS and the Complete Quadratic Combination rules are given and discussed in detail. Based on the c‐SRSS rule a response spectrum formulation is introduced to estimate the maximum structural response. An example considering a narrow‐band excitation from the great Mexico earthquake of September 19, 1985, is given and the accuracy of the response spectrum formulation is examined. Copyright © 2010 John Wiley & Sons, Ltd.     

Response of secondary systems in structures subjected to transient excitation

An analytical method, based on matrix perturbation theory, is developed whereby a simple estimate can be obtained of the maximum dynamic response of lightly damped, light equipment (modelled as a n(²)-degree-of-freedom system) attached to a structure (modelled as a n(¹)-degree-of-freedom system) subjected to ground motion or impact. A natural frequency of the equipment is considered close or equal to a natural frequency of the structure. It is assumed that the information available to the designer is a time history of the ground motion or impact, or an associated design spectrum; the fixed base modal properties of the structure; and the fixed base modal properties of the equipment. The method employed avoids the direct conventional analysis of a n(²) + n(¹)-degree-of-freedom system either by modal or by matrix time-marching methods; as well as errors in estimates of peak response due to the possible unreliability of numerical schemes because of the lightness of the equipment, or due to uncertainty as to the appropriate procedure for summing the contributions of the two closely spaced modes which occur in the system. The proposed procedure is demonstrated for an example equipment-structure system. Computed results based on the method are in close agreement with results obtained through a Newmark time-integration scheme.     

Blind modal identification of output‐only structures in time‐domain based on complexity pursuit

Output‐only modal identification is needed when only structural responses are available. As a powerful unsupervised learning algorithm, blind source separation (BSS) technique is able to recover the hidden sources and the unknown mixing process using only the observed mixtures. This paper proposes a new time‐domain output‐only modal identification method based on a novel BSS learning algorithm, complexity pursuit (CP). The proposed concept—independent ‘physical systems’ living on the modal coordinates—connects the targeted constituent sources (and their mixing process) targeted by the CP learning rule and the modal responses (and the mode matrix), which can then be directly extracted by the CP algorithm from the measured free or ambient system responses. Numerical simulation results show that the CP method realizes accurate and robust modal identification even in the closely spaced mode and the highly damped mode cases subject to non‐stationary ambient excitation and provides excellent approximation to the non‐diagonalizable highly damped (complex) modes. Experimental and real‐world seismic‐excited structure examples are also presented to demonstrate its capability of blindly extracting modal information from system responses. The proposed CP is shown to yield clear physical interpretation in modal identification; it is computational efficient, user‐friendly, and automatic, requiring little expertise interactions for implementations. Copyright © 2013 John Wiley & Sons, Ltd.     

Stochastic analysis of structures and piping systems subjected to stationary multiple support excitations

A stochastic method has been developed for seismic analysis of structures and piping systems subjected to multiple support excitations. In either the time or the frequency domain, mean and extreme values of structural and piping system response can be found, including the effects of cross-correlations of modal response and cross-correlations of multiple support excitations. Stationary white noise and stationary filtered white noise ground excitations are used. A computer program has been developed to carry out the stochastic seismic analysis. Results for a realistic nuclear power plant structure and piping system with and without modal cross-correlations and support excitation cross-correlations are compared. From these results, it is concluded that neglecting cross-correlations can lead to large errors. The stochastic method reported is shown to be more accurate than the response spectrum method and more economical than the time-history method; therefore, it is recommended for seismic analysis of nuclear power plants.     

关于结构振型参与系数和振型贡献的分析

采用振型分解反应谱法求解多自由度弹性体系的地震反应时,为了在满足所需计算精度的前提下减少工作量,需要对振型数量进行合理的选择,而振型数的确定主要取决于结构各阶振型对总体反应的贡献。通过数学推导,对振型贡献及振型数量的选择问题进行了研究。首先,讨论了振型参与系数的性质,在此基础上给出了能够反映结构振型贡献参数的数学表达式,对这些参数的力学含义进行了解释,并给出了相关证明;其次,对有效质量法、振型位移控制法等基于不同振型贡献标准的确定振型数的方法进行了分析比较,指出了其合理性和不足。本文研究对进一步理解结构振型贡献和振型数的选择问题具有一定的理论意义。     

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

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

Random response of a single-degree-of-freedom vibro-impact system with clearance

The closed form solutions of the stationary random response of a single-degree-of-freedom vibro-impact system with clearance are formulated in this paper. The Hertz contact law from elasticity is used to model the contact phenomena between the mass and constraint during vibration. The excitation is assumed to be a stationary white Gaussian process with zero mean. Through solving the time-independent Fokker-Planck equation, the stationary responses are obtained analytically. The effects of contact stiffness and clearance on the response are discussed probabilistically. It is found that, when the clearance is about twice the square root of the mean square response of the corresponding linear system, the contact phenomena are almost negligible.     

空间网格结构多维多点随机地震响应分析的高效算法

将林家浩教授提出的“虚拟激励法”进一步推广应用于空间网格结构多维多点非平稳随机地震响应分析,推导了多维虚拟激励随机振动分析方法的理论公式,给出了峰值响应估计方法,并讨论了多维地震动的随机模型及参数选取,通过编制的专用计算机程序分析了网壳结构的随机地震响应。本方法自动包含了参振振型间及各输人地震分量间的相关项,计算精确、快速,非常适合分析频率密集型空间网格结构的随机地震响应,是一种高效的随机振动分析算法。     

Stochastic considerations in seismic analysis of structures

A method is presented for stochastic modelling of a design earthquake by a power spectral density function for seismic analysis of structures. The method can be adopted with information currently available in the form of design response spectra for earthquake motion. Accurate seismic responses of structures can be easily obtained using such stochastic models. The methods for accurate response analysis of structures with closely spaced modes and for generation of floor response spectra of a building using a prescribed ground response spectrum directly are also presented. The hypothesis that a design earthquake can be modelled by a power spectral density function is used only implicitly in developing these methods.     

平稳随机地震地面运动过程模型及其统计特征

地震地面运动过程具有强烈的随机性,应用随机理论对实际工程结构进行地震可靠性分析和抗震设计与加固时都需要建立合理的随机地震地面运动模型,本文选择3种典型的随机地震动模型,即理想白噪声模型、金井清模型和改进的金井消模型,分析了它们的物理概念、频域特征以及适用范围。引入状态向量,建立状态方程．通过复振型叠加法分析了地震地面运动过程的时域统计特性,推导出3种随机地震动模型的相关函数的解析表达式．这些结果可为结构随机地震反应时域分析和抗震可靠性评估提供基础。     

Modal combination rules for multicomponent earthquake excitation

A response spectrum method for dynamic analysis of linear structures subjected to multicomponent seismic input is developed. The method is based on elementary concepts of stationary random vibration and assumes the existence of a set of principal directions along which components of ground motion are uncorrelated. Modal combination rules in terms of ground response spectra are developed for the mean and standard deviation of peak responses and for rootmean-square responses. These rules account for correlations between modal responses of the structure, as well as correlations between the input components. When the position of principal directions is unknown, two alternative rules are proposed: one uses the direction which is most critical for the response quantity of interest, and the other considers the direction as a random variable. The proposed method is simple for practical implementation and gives more accurate results than other existing methods.     

A simplified multisupport response spectrum method

A simplified multisupport response spectrum method is presented.The structural response is a sum of two components of a structure with a first natural period less than 2 s.The first component is the pseudostatic response caused by the inconsistent motions of the structural supports,and the second is the structural dynamic response to ground motion accelerations.This method is formally consistent with the classical response spectrum method,and the effects of multisupport excitation are considered for any modal response spectrum or modal superposition.If the seismic inputs at each support are the same,the support displacements caused by the pseudostatic response become rigid body displacements.The response spectrum in the case of multisupport excitations then reduces to that for uniform excitations.In other words,this multisupport response spectrum method is a modification and extension of the existing response spectrum method under uniform excitation.Moreover,most of the coherency coefficients in this formulation are simplified by approximating the ground motion excitation as white noise.The results indicate that this simplification can reduce the calculation time while maintaining accuracy.Furthermore,the internal forces obtained by the multisupport response spectrum method are compared with those produced by the traditional response spectrum method in two case studies of existing long-span structures.Because the effects of inconsistent support displacements are not considered in the traditional response spectrum method,the values of internal forces near the supports are underestimated.These regions are important potential failure points and deserve special attention in the seismic design of reticulated structures.     

Defining equivalent stationary PSDF to account for nonstationarity of earthquake ground motion

This paper presents three approaches to defining the stationary power spectrum density function (PSDF) of strong ground acceleration, for prediction of structural response corresponding to the strong-motion stationary part of the input excitation. The first approach defines the PSDF in terms of the Fourier amplitude spectrum and a stationary duration of ground acceleration. The PSDF obtained by this approach predicts accurately the response of structures with low to intermediate natural periods. In the second approach, we introduce the concept of stationary duration of response, which is defined as a function of the natural period and damping ratio of the oscillator. Using this approach, it is possible to get accurate estimates of response amplitudes for the broad range of natural periods. However, it is not convenient in practical applications to deal with several stationary durations for a given input excitation. Further, to evaluate these durations it is necessary to specify both the Fourier and the response spectra of ground accelerations; whereas the common engineering practice is to specify the response spectrum only. Therefore, the third approach suggests the use of the response ‘spectrum compatible’ PSDF. The paper presents several improvements in the general methodology used for this purpose. The improvements mainly relate to using more accurate peak factors and to using the transient nature of response. The spectrum compatible PSDFs, as evaluated in the present study, provide realistic specification of strong ground motion for stochastic seismic response analyses of structures.