设有液压质量控制系统(HMS)的底层柔性结构动力特性及其地震反应的振型分解法

本文对设有液压质量控制系统(HMS)的底层柔性结构体系的动力特性进行了试验研究和分析,通过对试验结果的分析,证明了设有液压质量控制系统(HMS)的底层柔性结构控制体系近似存在经典振型,其阻尼阵满足正交性条件,所以,可用振型分解法计算受控结构的地震反应。     

Damped vibration mode superposition method for dynamic response analysis

In the dynamic response analysis of extremely complex structural systems in which the damping characteristics of each element are independent, the damping matrix is not always diagonalized by the use of undamped free vibration mode shapes. In the present paper, a mode-superposition method by the use of damped free vibration mode shapes is developed for such structural systems. It is also shown that the Fast Fourier Transform (FFT) procedures, that are available for the dynamic response analysis of linear structural systems, are used effectively in this mode-superposition method with good accuracy.     

Efficient mode superposition algorithm for seismic analysis of non-linear structures

A step-by-step integration method is proposed to compute within the framework of the conventional mode superposition technique the response of bilinear hysteretic structures subjected to earthquake ground motions. The method is computationally efficient because only a few modes are needed to obtain an accurate estimate of such a response, and because it does not require the use of excessively small time steps to avoid problems of accuracy or stability. It is developed on the basis that the non-linear terms in the equations of motion for non-linear systems may be considered as additional external forces, and the fact that by doing so such equations of motion can be interpreted as the equations of motion of an equivalent linear system, excited by a modified ground motion. These linear equations are then subjected to a conventional modal decomposition and transformed, as with linear systems, into a set of independent differential equations, each representing the system's response in one of its modes of vibration. To increase the efficiency of the method and account properly for the participation of higher modes, these independent equations are solved using the Nigam-Jennings technique in conjunction with the so-called mode acceleration method. The accuracy and efficiency of the method is verified by means of a comparative study with solutions obtained with a conventional direct integration method. In this comparative study, including only a few modes, the proposed method accurately predicts the seismic response of three two-dimensional frame structures, but requiring only, on an average, about 47 per cent less computer time than when the direct integration method is used.     

弹性能量半径演化谱及其在多自由度体系地震反应计算中的应用

目前用于结构抗震设计的反应谱仅能反映峰值反应,无法体现反应值随时间的变化。文中提出一种弹性能量半径演化谱,可反映线性单自由度体系弹性能量(即动能与弹性势能之和)随地震持时的变化,且其峰值近似等于结构峰值位移。文中给出了利用地震动演化功率谱得到该演化谱的方法并进一步发展了一种计算线性多自由度体系地震位移反应的新方法。通过两座框架结构的地震反应计算,将新方法与传统振型组合法及时程分析法的计算结果进行对比,发现对于振型稀疏结构,新方法计算结果与SRSS法接近;而对于振型密集结构,新方法计算结果较CQC法更精确,且避免了CQC法相关系数的复杂计算。     

Dynamic response of an elastic-viscoplastic system in modal co-ordinates

A dynamic analysis of elastic–viscoplastic systems, incorporating the modal co-ordinate transformation technique, is presented. The formulation results in uncoupled incremental equations of motion with respect to the modal co-ordinates. The elastic–viscoplastic model adopted allows the analysis not to involve yielding regions and loading/unloading processes. An implicit Runge–Kutta scheme together with the Newton–Raphson method are used to solve the non-linear constitutive equations. Stability and accuracy of the numerical solution are improved by utilizing a local time step sub-incrementing procedure. Applications of the analyses to multi-storey shear buildings show that good results can be obtained for the maximum displacement response by including only a few lower modes in the computation, but the prediction of the ductility factor response tends to underestimate the peak values when too few modes are used. In addition, stable and valid results can be obtained even with a sizable time step increment.     

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.     

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

A method, based on the Hilbert–Huang spectral analysis, has been proposed by the authors to identify linear structures in which normal modes exist (i.e., real eigenvalues and eigenvectors). Frequently, all the eigenvalues and eigenvectors of linear structures are complex. In this paper, the method is extended further to identify general linear structures with complex modes using the free vibration response data polluted by noise. Measured response signals are first decomposed into modal responses using the method of Empirical Mode Decomposition with intermittency criteria. Each modal response contains the contribution of a complex conjugate pair of modes with a unique frequency and a damping ratio. Then, each modal response is decomposed in the frequency–time domain to yield instantaneous phase angle and amplitude using the Hilbert transform. Based on a single measurement of the impulse response time history at one appropriate location, the complex eigenvalues of the linear structure can be identified using a simple analysis procedure. When the response time histories are measured at all locations, the proposed methodology is capable of identifying the complex mode shapes as well as the mass, damping and stiffness matrices of the structure. The effectiveness and accuracy of the method presented are illustrated through numerical simulations. It is demonstrated that dynamic characteristics of linear structures with complex modes can be identified effectively using the proposed method. Copyright © 2003 John Wiley & Sons, Ltd.     

非经典振型分解法求混凝土房屋及其上钢塔的地震响应

钢结构与混凝土结构阻尼比不同,混凝土房屋与其顶上钢塔组成了非比例阻尼结构系统。本文用非经典振型分解法求解该类结构系统的线弹性地震响应,发现只用前几阶振型响应迭加的结果即可逼近直接积分法的精确度。     

Evaluation of seismic site factors in the Mississippi Embayment. I. Estimation of dynamic properties

Deep deposits of the Mississippi Embayment, overlying the New Madrid Seismic Zone, present unique challenges for the estimation of local site effects on propagated ground motion. Limited information is available on the dynamic properties of these deposits. This paper develops generalized depth dependent modulus degradation and damping curves specific to the embayment deposits. Depth dependent small strain damping is estimated using weak motion recordings during the Enola earthquake, 2001. Modulus degradation and damping curves are further constrained using limited laboratory test data of embayment soils at low confining pressure. At high confining pressures embayment specific data is unavailable and test data from outside the embayment is used in developing the dynamic properties. The representative modulus degradation and damping curves are used in non-linear and equivalent linear one-dimensional site response analyses. A companion paper describes a large-scale probabilistic seismic hazard analysis study in the Mississippi Embayment that integrates non-linear site effects.     

Dynamic property and seismic performance of square-based membranous latticed shells

The basic natural vibration property of the square-based single layer membranous latticed shell is analyzed using the subspace iteration method. Then, the seismic performances of the structure are studied in detail by the mode-superposition response spectrum method, and the reasonable combination number of vibration modes in seismic analysis by the mode-superposition response spectrum method is investigated. In addition, the formulas of seismic time-history response of a long-span spatial structure are deduced considering geometrical non-linearity, and the corresponding computer program is developed. The results by the response spectrum method and those by the time-history method are, then, compared. The researches show that, the seismic performance of the square-based membranous latticed shell is comparatively well, and the reasonable combination number of vibration modes is 15 in the vertical seismic analysis of this kind of structure by the mode-superposition response spectrum method, whereas, it is 60 in the horizontal seismic analysis. Moreover, the results of the vertical seismic analysis of the structure by the mode-superposition response spectrum method are on the safe side, whereas, in the horizontal seismic analysis of the structure, the results must be checked by the time-history method at the end. All this work will provide theoretical guidance to the actual engineering application of this structure.     

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.     

Finite element analysis of water flow in variably saturated soil

A two-dimensional Galerkin finite element model for water flow in variably saturated soil is presented. A fourth-order Runge-Kutta time integration method is employed which allows use of time steps at least 2 times greater than for a traditional finite difference approximation of time derivatives. For short total simulation times computer execution costs for the Runge-Kutta method are greater than for the finite difference approximation due to the start up cost of the Runge-Kutta method, but for longer simulation times the Runge-Kutta method requires considerably less computational effort even when automatic time-step adjustment is used with the finite difference procedure. A comparison of the method of influence coefficients and 2 × 2 Gaussian integration to compute element matrices indicates that the influence coefficient method reduces total execution time to 60% of that required for numerical quadrature. Computed pressure heads using the influence coefficient method and numerical integration are found to be in close agreement with each other even under conditions of highly non-linear soil properties in a heterogeneous domain. Fluxes computed by the two methods are also generally in close agreement except under extremely non-linear conditions when some deviations were observed at short simulation times.     

Seismic response of offshore structures in random seas

A study of the dynamic response of offshore structures to simultaneous loadings by random earthquake ground motions and random sea waves is presented. Emphasis is placed on the evaluation of dynamic soil-structure interaction effects and also on the evaluation of non-linear hydrodynamic damping effects due to sea waves for the seismic response. The structure is discretized using the finite element method. Sea waves are represented by Bretschneider's power spectrum and the Morison equation defines the wave forcing function. The Tajimi-Kanai power spectrum is used for the horizontal ground acceleration due to earthquakes. The governing equations of motion are obtained by the substructure method. Response analysis is carried out using the frequency-domain random vibration approach. It is found that the first few vibrational modes contribute significantly to the dynamic response. The response due to earthquake loadings is larger when the soil-structure interaction effects are considered. The hydrodynamic damping forces are higher in random seas than in still water and sea waves reduce the seismic response of offshore structures. Studies on the first passage probabilities of response indicate that small sea waves enhance the reliability of offshore structures against earthquake forces.     

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.     

Elimination of spurious higher-mode response in pseudodynamic tests

In a pseudodynamic test, errors in restoring-force feedback are introduced into numerical computations. Some of these errors can excite the higher-frequency response of the specimen. In this paper, the use of viscous and numerical dampings to eliminate spurious higher-frequency effects is studied. Since the tangent stiffness of a non-linear specimen cannot be measured accurately, initial-stiffness-dependent viscous damping is considered. In addition, an explicit integration algorithm with desired numerical damping properties is proposed and examined. The analytical and numerical studies presented indicate that viscous-damping properties can be substantially changed by non-linear deformations. For this reason, the use of numerical damping appears to be more advantageous.     

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

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

The ρ-family of algorithms for time-step integration with improved numerical dissipation

The dynamic analysis of complex non-linear structural systems by the finite element approach requires the use of time-step algorithms for solving the equations of motion in the time domain. Both an implicit and an explicit version of such a time-step algorithm, called the ρ-method, the parameter ρ being used for controlling numerical damping in the higher modes, are presented in this paper. For the implicit family of algorithms unconditional stability, consistency, convergence, accuracy and overshoot properties are first discussed and proved. On the basis of the algorithmic damping ratio (dissipation) and period elongation (dispersion) the ρ-method is then compared with the well-known implicit algorithms of Hilber, Newmark, Wilson, Park and Houbolt. An explicit version of the algorithm is also derived and briefly discussed. This shows numerical properties similar to the central difference method. Both versions of the algorithm have been implemented in a general purpose computer program which has been often used for both numerical tests and practical applications.     

Appropriate viscous damping for nonlinear time‐history analysis of base‐isolated reinforced concrete buildings

There is no consensus at the present time regarding an appropriate approach to model viscous damping in nonlinear time‐history analysis of base‐isolated buildings because of uncertainties associated with quantification of energy dissipation. Therefore, in this study, the effects of modeling viscous damping on the response of base‐isolated reinforced concrete buildings subjected to earthquake ground motions are investigated. The test results of a reduced‐scale three‐story building previously tested on a shaking table are compared with three‐dimensional finite element simulation results. The study is primarily focused on nonlinear direct‐integration time‐history analysis, where many different approaches of modeling viscous damping, developed within the framework of Rayleigh damping are considered. Nonlinear direct‐integration time‐history analysis results reveal that the damping ratio as well as the approach used to model damping has significant effects on the response, and quite importantly, a damping ratio of 1% is more appropriate in simulating the response than a damping ratio of 5%. It is shown that stiffness‐proportional damping, where the coefficient multiplying the stiffness matrix is calculated from the frequency of the base‐isolated building with the post‐elastic stiffness of the isolation system, provides reasonable estimates of the peak response indicators, in addition to being able to capture the frequency content of the response very well. Furthermore, nonlinear modal time‐history analyses using constant as well as frequency‐dependent modal damping are also performed for comparison purposes. It was found that for nonlinear modal time‐history analysis, frequency‐dependent damping, where zero damping is assigned to the frequencies below the fundamental frequency of the superstructure for a fixed‐base condition and 5% damping is assigned to all other frequencies, is more appropriate, than 5% constant damping. Copyright © 2013 John Wiley & Sons, Ltd.     

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.