Efficient analysis of dynamic response of linear systems

After reviewing briefly a recently proposed procedure for evaluating the dynamic transient response of a classically damped linear system from its corresponding steady-state response, a modified procedure is presented which also appears to be highly efficient for non-classically damped systems of the type encountered in studies of soil-structure interaction. The concepts involved are developed by reference to viscously damped single-degree-of-freedom systems, and numerical solutions are included to illustrate the accuracy and efficiency of the proposed procedure and its superiority over the classical Discrete Fourier Transform approach.     

Mode-superposition methods in dynamic analysis of classically and non-classically damped linear systems

Mode-superposition analysis is an efficient tool for the evaluation of the response of linear systems subjected to dynamic agencies. Two well-known mode-superposition methods are available in the literature, the mode-displacement method and the mode-acceleration method. Within this frame a method is proposed called a dynamic correction method which evaluates the structural response as the sum of a pseudostatic response, which is the particular solution of the differential equations, and a dynamic correction evaluated using a reduced number of natural modes. The greater accuracy of the proposed method with respect to the other methods is evidenced through extensive numerical tests, for classically and non-classically damped systems.     

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.     

Mode‐acceleration approach to seismic response of multiply‐supported secondary systems

A mode‐acceleration approach has been proposed for estimating the seismic response of a linear, classically‐damped, multiply‐supported secondary system within the framework of a power spectral density function (PSDF)‐based stochastic approach, while the primary system is linear and classically‐damped. Response transfer functions have been formulated in terms of chosen numbers of fixed‐base modes of the primary and secondary systems. The proposed approach does not involve the determination of combined system properties, and is applicable to the secondary systems with high mass ratios also. Through a few example primary–secondary systems and an example band‐limited white noise excitation, it has been shown that this approach leads to reasonably accurate results when only a few primary and secondary modes are to be considered. The proposed formulation has been used to obtain input data for a decoupled response spectrum analysis of secondary systems. This data accurately accounts for the effects of interaction between the primary and secondary systems. It is shown to lead to substantial reductions in the errors associated with the envelope spectrum method in the case of moderately heavy to heavy secondary systems and when the spatial coupling does not play a major role. Copyright © 2002 John Wiley & Sons, Ltd.     

Time delayed control of structural systems

Time delays are ubiquitous in control systems. They usually enter because of the sensors and actuators used in them. Traditionally, time delays have been thought to have a deleterious effect on both the stability and the performance of controlled systems, and much research has been done in attempting to eliminate them, compensate for them, or nullify their presence. In this paper we take a different view. We investigate whether purposefully injected time delays can be used to improve both the system's stability and performance. Our analytical, numerical, and experimental investigation shows that this can indeed be done. Analytical results of the effects of time delays on collocated and non‐collocated control of classically damped and non‐classically damped systems are given. Experimental and numerical results confirm the theoretical expectations. Issues of system uncertainties and robustness of time delayed control are addressed. The results are of practical value in improving the performance and stability of controllers because these characteristics (performance and stability) improve dramatically with the intentional injection of small time delays in the control system. The introduction of such time delays constitutes a ‘minimal change’ to a controller already installed in a structural system for active control. Hence, from a practical standpoint, time delays can be implemented in a nearly costless and highly reliable manner to improve control performance and stability, an aspect that cannot be ignored when dealing with the economics and safety of large structural systems subjected to strong earthquake ground shaking. Copyright © 2003 John Wiley & Sons, Ltd.     

Design of friction damped structures using lateral force procedure

A trilinear model is used to simulate the seismic resisting mechanism of a single‐degree‐of‐freedom friction‐damped system to reflect the situation in which both dampers and frame members lose their elastic resistance. The seismic response of the friction‐damped system is normalized with respect to the response of its corresponding linear system by an approach that incorporates a credible equivalent linearization method, a damping reduction rule and the algebraic specification of the design spectrum. The resulting closed‐form solutions obtained for the normalized response are then used to define a force modification factor for friction‐damped systems. This force modification factor, together with the condensation procedure for multi‐degree‐of‐freedom structures, enables the establishment of a quasi‐static design procedure for friction‐damped structures, which is intended for the benefit and use of structural practitioners. A curve‐fitting technique is employed to develop an explicit expression for the force modification factor used with the proposed design procedure; it is shown that this simplification results in satisfactory accuracy. Finally, a design example is given to illustrate the validation of the proposed design procedure. Copyright © 2000 John Wiley & Sons, Ltd.     

Improved dynamic correction method in seismic analysis of both classically and non‐classically damped structures

A step‐by‐step approximate procedure taking into consideration high‐frequency modes, usually neglected in the modal analysis of both classically and non‐classically damped structures, is presented. This procedure can be considered as an extension of traditional modal correction methods, like the mode‐acceleration method and the dynamic correction method, which are very effective for structural systems subjected to forcing functions described by analytical laws. The proposed procedure, herein called improved dynamic correction method, requires two steps. In the first step, the number of differential equations of motion are reduced and consequently solved by using the first few undamped mode‐shapes. In the second step, the errors due to modal truncation are reduced by correcting the dynamic response and solving a new set of differential equations, formally similar to the original differential equations of motion. The difference between the two groups of differential equations lies in the forcing vector, which is evaluated in such a way as to correct the effects of modal truncation on applied loads. Copyright © 2001 John Wiley & Sons, Ltd.     

On a modal damping identification model for non-classically damped linear building structures subjected to earthquakes

A simple modal damping identification model developed by the present authors for classically damped linear building frames is extended here to the non-classically damped case. The modal damping values are obtained with the aid of the frequency domain modulus of the roof-to-basement transfer function and the resonant frequencies of the structure (peaks of the transfer function) as well as the modal participation factors and mode shapes of the undamped structure. The assumption is made that the modulus of the transfer function of the non-classically damped structure matches the one of the classically damped structure in a discrete manner, i.e., at the resonant frequencies of that function modulus. This proposed approximate identification method is applied to a number of plane building frames with and without pronounced non-classical damping under different with respect to their frequency content earthquakes and its limitations and range of applicability are assessed with respect to the accuracy of both the identified damping ratios and that of the seismic structural response obtained by classical mode superposition and use of those identified modal damping ratios.     

Partial eigensolution of damped structural systems by Arnoldi's method

An efficient numerical algorithm is developed to solve the quadratic eigenvalue problems arising in the dynamic analysis of damped structural systems. The algorithm can even be applied to structural systems with non-symmetric matrices. The algorithm is based on the use of Arnoldi's method to generate a Krylov subspace of trial vectors, which is then used to reduce a large eigenvalue problem to a much smaller one. The reduced eigenvalue problem is solved and the solutions are used to construct approximate solutions to the original large system. In the process, the algorithm takes full advantage of the sparseness and symmetry of the system matrices and requires no complex arithmetic, therefore, making it very economical for use in solving large problems. The numerical results from test examples are presented to demonstrate that a large fraction of the approximate solutions calculated are very accurate, indicating that the algorithm is highly effective for extracting a number of vibration modes for a large dynamic system, whether it is lightly or heavily damped.     

Response spectrum analysis for non-classically damped linear system with multiple-support excitations

A new response spectrum method, which is named complex multiple-support response spectrum (CMSRS) method in this article, is developed for seismic analysis of non-classically damped linear system subjected to spatially varying multiple-supported ground motion. The CMSRS method is based on fundamental principles of random vibration theory and properly accounts for the effect of correlation between the support motions as well as between the modal displacement and velocity responses of structure, and provides an reasonable and acceptable estimate of the peak response in term of peak seismic ground motions and response spectra at the support points and the coherency function. Meanwhile, three new cross-correlation coefficients or cross covariance especially for the non-classically damped linear structures with multiple-supports excitations are derived under the same assumptions of the MSRS method of classically damped system. The CMSRS method is examined and compared to the results of time history analyses in two numerical examples of non-classically damped structures in consideration of the coherences of spatially variable ground motion. The results show that for non-classically damped structure, the cross terms representing the cross covariance between the pseudo-static and dynamic component are also quite small just as same as classically damped system. In addition, it is found that the usual way of neglecting all the off-diagonal elements in transformed damping matrix in modal coordinates in order to make the concerned non-classically damped structure to become remaining proportional damping property will bring some errors in the case of subjected to spatially excited inhomogeneous ground motion.     

Asymmetric one‐storey elastic systems with non‐linear viscous and viscoelastic dampers: Simplified analysis and supplemental damping system design

Investigated is the accuracy in estimating the response of asymmetric one‐storey systems with non‐linear viscoelastic (VE) dampers by analysing the corresponding linear viscous system wherein all non‐linear VE dampers are replaced by their energy‐equivalent linear viscous dampers. The response of the corresponding linear viscous system is determined by response history analysis (RHA) and by response spectrum analysis (RSA) extended for non‐classically damped systems. The flexible and stiff edge deformations and plan rotation of the corresponding linear viscous system determined by the extended RSA procedure is shown to be sufficiently accurate for design applications with errors generally between 10 and 20%. Although similar accuracy is also shown for the ‘pseudo‐velocity’ of non‐linear VE dampers, the peak force of the non‐linear VE damper cannot be estimated directly from the peak damper force of the corresponding linear viscous system. A simple correction for damper force is proposed and shown to be accurate (with errors not exceeding 15%). For practical applications, an iterative linear analysis procedure is developed for determining the amplitude‐ and frequency‐dependent supplemental damping properties of the corresponding linear viscous system and for estimating the response of asymmetric one‐storey systems with non‐linear VE dampers from the earthquake design (or response) spectrum. Finally, a procedure is developed for designing non‐linear supplemental damping systems that satisfy given design criteria for a given design spectrum. Copyright © 2003 John Wiley & Sons, Ltd.     

Steady-state and transient responses of non-classically damped linear systems

A recently proposed procedure for interrelating the steady-state and transient responses of multi-degree-of-freedom, classically damped linear systems is extended to non-classically damped systems. The extension is formulated for baseexcited systems, and it is illustrated by simple examples.     

Earthquake analysis of concrete gravity dams including dam-water-foundation rock interaction

A general procedure for analysis of the response of concrete gravity dams, including the dynamic effects of impounded water and flexible foundation rock, to the transverse (horizontal) and vertical components of earthquake ground motion is presented. The problem is reduced to one in two dimensions, considering the transverse vibration of a monolith of the dam. The system is analysed under the assumption of linear behaviour for the concrete, foundation rock and water. The complete system is considered as composed of three substructures—the dam, represented as a finite element system, the fluid domain, as a continuum of infinite length in the upstream direction, and the foundation rock region as a viscoelastic half-plane. The structural displacements of the dam are expressed as a linear combination of Ritz vectors, chosen as normal modes of an associated undamped dam-rock system. The effectiveness of this analytical formulation lies in its being able to produce excellent results by considering only a few Ritz vectors. The generalized displacements due to earthquake motion are computed by synthesizing their complex frequency responses using Fast Fourier Transform procedures. The stress responses are calculated from the displacements. An example analysis is presented to illustrate results obtained from this analytical procedure. Computation times for several analyses are presented to illustrate the effectiveness of the procedure.     

Identification of modal parameters of non‐classically damped linear structures under multi‐component earthquake loading

A method for parametric system identification of classically damped linear system in frequency domain is adopted and extended for non‐classically damped linear systems subjected up to six components of earthquake ground motions. This method is able to work in multi‐input/multi‐output (MIMO) case. The response of a two‐degree‐of‐freedom model with non‐classical damping, excited by one‐component earthquake ground motion, is simulated and used to verify the proposed system identification method in the single‐input/multi‐output case. Also, the records of a 10 storey real building during the Northridge earthquake is used to verify the proposed system identification method in the MIMO case. In this case, at first, a single‐input/multi‐output assumption is considered for the system and modal parameters are identified, then other components of earthquake ground motions are added, respectively, and the modal parameters are identified again. This procedure is repeated until all four components of earthquake ground motions which are measured at the base level of the building are included in the identification process. The results of identification of real building show that consideration of non‐classical damping and inclusion of the multi‐components effect of earthquake ground motions can improve the least‐squares match between the finite Fourier transforms of recorded and calculated acceleration responses. Copyright © 2006 John Wiley & Sons, Ltd.     

New approaches for non‐classically damped system eigenanalysis

This paper presents three new approaches for solving eigenvalue problems of non‐classically damped linear dynamics systems with fewer calculations than the conventional state vector approach. In the latter, the second‐order differential equation of motion is converted into a first‐order system by doubling the size of the matrices. The new approaches simplify the approach and reduce the number of calculations. The mathematical formulations for the proposed approaches are presented and the numerical results compared with the existing method by solving a sample problem with different damping properties. Of the three proposed approaches, the expansion approach was found to be the simplest and fastest to compute. Copyright © 2005 John Wiley & Sons, Ltd.     

非比例阻尼线性体系地震反应计算的振型分解反应谱法

以非比例阻尼线性体系地震反应计算实数形式的一般解答为基础,推导得到了非比例阻尼线性体系水平地震作用计算的多种形式,建立了非比例阻尼线性体系地震反应计算振型分解反应谱法的基本过程与步骤。最后,以一个五层剪切型结构为例,通过与各种常用直接积分方法计算结果的比较,证实了本文非比例阻尼线性体系地震反应计算实数形式的一般解答的高精度与可靠性。通过对多种形式地震作用所得地震效应的比较,证实了非比例阻尼线性体系地震反应振型分解反应谱方法的可靠性及可行性。     

Dynamics of flexible systems with partial lift-off

The dynamic behaviour of a simplified model of a multi-storey building, supported by an elastic foundation and allowed to uplift, is examined. The building is modelled by an n-degree-of-freedom oscillator, while the foundation is represented by a viscously damped two-spring model which permits uplift. This model has been shown in previous studies to be an accurate approximation to the more realistic but more complex Winkler foundation. Approximate values for the characteristic frequencies of the interacting system are presented and a simple, first-mode solution is developed. The response of the system is non-linear and the apparent fundamental period increases with the amount of lift-off. In contrast to the first mode of the superstructure which participates strongly in the interaction, the second and higher modes of the building are not affected significantly by either the interaction with the soil or the uplift. The study shows that lift-off results in larger rocking motion of the structure, but it is not clear from the analysis and the example whether the interfloor displacements are consistently increased or decreased, since this appears to depend on the properties of the system and the excitation.     

Quantification of ground-motion parameters and response spectra in the near-fault region

This study focuses on the characteristics of near-fault ground motions in the forward-direction and structural response associated with them. These ground motions are narrow-banded in nature and are characterized by a predominant period at which structures excited by them are severely affected. In this work, predominant period is defined as the undamped natural period of a single-degree-of-freedom (SDOF) oscillator at which its 5% damped linear elastic pseudo-spectral velocity (PSV) contains a clear and dominant peak. It is found that a linear relationship exists between predominant period and seismic moment. An empirical equation describing this relationship is presented by using a large set of accelerograms. Attenuation equations are developed to estimate peak ground velocity (PGV) as a function of earthquake magnitude and source-to-site distance. In addition, a predictive equation for spectral shapes of PSV (i.e., PSV normalized by PGV) is presented as a continuous function of the undamped natural period of SDOF oscillators. The model is independent of PGV, and can be used in conjunction with any available PGV attenuation relation applicable to near-fault ground motion exhibiting forward-directivity effects. Furthermore, viscous damping of the SDOF is included in the model as a continuous parameter, eliminating the use of so-called damping correction factors. Finally, simple equations relating force reduction factors and displacement ductility of elasto-plastic SDOF systems are presented.     

Active control of MDOF structures with supplemental electrorheological fluid dampers

A method for design of an active control system for multistorey structures using Electrorheological (ER) dampers is presented. Incorporated at various levels of a structural frame, ER dampers are used to improve the response of the structure during earthquakes. Optimal control theory was used to design the ER devices. The aim of the design is to find the most suitable combination of the minimum required forces produced by the ER dampers to obtain the optimal structural response. The mechanical response of ER fluid dampers is regulated by an electric field, depending on the displacements and velocities of the frame. Numerical analysis of an ER damped seven-storey structure is represented as an example. Significant improvement of the structural response was obtained using optimal active controlled ER dampers compared to passive controlled and uncontrolled structures. Copyright © 1999 John Wiley & Sons, Ltd.     

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.