Shear hysteretic elasto-plastic earthquake response of soil systems

An analytical/numerical procedure is developed to calculate the shear hysteretic elasto-plastic earthquake response of soil systems. Dynamic level ground response as well as earth-dam response is evaluated on the basis of the analytical one-dimensional (ID) shear-beam and shear-wedge models respectively. Soil is modelled as an elasto-plastic hysteretic material following the principles of flow or incremental plasticity. Stiffness and strength degradation of this material may be specified as a function of accumulated plastic strain and associated pore-pressure build-up. The entire numerical procedure may be executed on a personal computer and consequently qualifies as a versatile tool in conducting preliminary design calculations or parametric-type investigations. Using this procedure, the responses of the Mexico City SCT site during the September 19, 1985 earthquake, and the California Wildlife site during the November 24, 1987 earthquake are computed and compared to those actually observed.     

Seismic analysis of structures with a fractional derivative model of viscoelastic dampers

Viscoelastic dampers are now among some of the preferred energy dissipation devices used for passive seismic response control. To evaluate the performance of structures installed with viscoelastic dampers, different analytical models have been used to characterize their dynamic force deformation characteristics. The fractional derivative models have received favorable attention as they can capture the frequency dependence of the material stiffness and damping properties observed in the tests very well. However, accurate analytical procedures are needed to calculate the response of structures with such damper models. This paper presents a modal analysis approach, similar to that used for the analysis of linear systems, for solving the equations of motion with fractional derivative terms for arbitrary forcing functions such as those caused by earthquake induced ground motions. The uncoupled modal equations still have fractional derivatives, but can be solved by numerical or analytical procedures. Both numerical and analytical procedures are formulated. These procedures are then used to calculate the dynamic response of a multi-degree of freedom shear beam structure excited by ground motions. Numerical results demonstrating the response reducing effect of viscoelastic dampers are also presented.     

A new spatial coherence model and analytical coefficients for multi-support response spectrum combination

In this paper, a new spatial coherence model of seismic ground motions is proposed by a fi tting procedure. The analytical expressions of modal combination (correlation) coeffi cients of structural response are developed for multi-support seismic excitations. The coeffi cients from both the numerical integration and analytical solutions are compared to verify the accuracy of the solutions. It is shown that the analytical expressions of numerical modal combination coeffi cients are of high accuracy. The results of random responses of an example bridge show that the analytical modal combination coeff icients developed in this paper are accurate enough to meet the requirements needed in practice. In addition, the computational effi ciency of the analytical solutions of the modal combination coeff icients is demonstrated by the response computation of the example bridge. It is found that the time required for the structural response analysis by using the analytical modal combination coeffi cients is less than 1/20 of that using numerical integral methods.     

Numerical artifacts associated with Rayleigh and modal damping models of inelastic structures with massless coordinates

The paper presents a detailed reexamination of the effects of three damping models on the inelastic seismic response of structures with massless degrees of freedom. The models considered correspond to (a) Rayleigh damping based on current properties (tangent stiffness), (b) Rayleigh damping based on initial properties, and (c) modal damping. The results suggest that some nonzero damping forces/moments at massless DOFs obtained in multistory frames for the case of Rayleigh damping with tangent stiffness may be numerical artifacts rather than a deficiency of the damping model. The results also indicate that significant artificial numerical oscillations in the velocities of the secondary components of MDOF structures are introduced when modal damping or mass-proportional damping is used.     

Computational aspects of a seismic performance test method using on-line computer control

The basic computational procedure of a newly developed seismic performance test méthod is presented and examined in this paper. This method uses a direct step-by-step integration technique to compute the displacement response of a test specimen subjected to a numerically specified seismic excitation record, utilizing the non-linear restoring forces actually developed by the specimen during the experiment. Due to the limitations of the experimental procedure and conditions, special analytical assumptions and numerical methods must be adopted. The adequacy of such analytical techniques is evaluated and possible computational errors are identified. An improved numerical integration algorithm is also proposed for this specific application. It has been found that the method can be a reliable and powerful experimental tool provided appropriate analytical and numerical criteria, as discussed in this paper, are observed.     

Approximate soil-structure interaction analysis by a perturbation approach: The case of stiff soils

An approximate solution of the classical eigenvalue problem governing the vibrations of a structure on an elastic soil is derived through the application of a perturbation analysis. For stiff soils, the full solution is obtained as the sum of the solution for a rigid-soil and small perturbing terms related to the inverse of the soil shear modulus. The procedure leads to approximate analytical expressions for the system frequencies, modal damping ratios and participation factors for all system modes that generalize those presented by other authors for the fundamental mode. The resulting approximate expressions for the system modal properties are validated by comparison with the corresponding quantities obtained by numerical solution of the eigenvalue problem for a nine-story building. The accuracy of the proposed approach and of the classical normal mode approach is assessed through comparison with the exact frequency response of the test structure.     

面向设计应用的地震动空间相干函数模型

本文对现有的常用地震动空间相干模型进行了总结,提出了一个新的面向工程抗震设计应用的形式统一的地震动空间相干函数模型,在此基础上推导出了多点地震反应谱和功率谱计算所需要的振型组合系数的解析表达式,避免了耗费时间的数值积分运算。本文模型与计算方法使多点地震激励下结构响应的计算时间减低至积分方法的1/20以下,使多点地震反应谱方法和多点地震功率谱方法在计算时间方面实用化。     

Static equivalent method for the kinematic interaction analysis of single piles

This paper presents a static equivalent approach to estimate the maximum kinematic interaction effects on piles subjected to lateral seismic excitation. Closed-form expressions are reported for the evaluation of the maximum free-field soil movements and for the computation of maximum pile shear force and bending moments. Firstly, modal analysis, combined with a suitable damped response spectrum, is used to evaluate the maximum free-field response. Secondly, the pile is schematised as a Winkler's beam subjected to equivalent static forces defined according to soil vibration modal shapes and amplitude. The method may be applied by using response spectra suggested by National Standards or those obtained with accelerograms. The procedure proposed may be conveniently implemented in simple spreadsheets or in commercial finite element programs and easily used by practicing engineers. Method accuracy is demonstrated by comparing the results with those obtained with a more rigorous model. Good results may be achieved by considering only the first soil vibration mode making the procedure straightforward for practical design purposes.     

Approximate modal decomposition of inelastic dynamic responses of wall buildings

Two approximate methods for decomposing complicated inelastic dynamic responses of wall buildings into simple modal responses are presented. Both methods are based on the equivalent linear concept, where a non‐linear structure is represented by a set of equivalent linear models. One linear model is used for representing only one vibration mode of the non‐linear structure, and its equivalent linear parameters are identified from the inelastic response time histories by using a numerical optimizer. Several theoretical relations essential for the modal decomposition are derived under the framework of complex modal analysis. Various numerical examinations have been carried out to check the validity of the proposed modal decomposition methods, and the results are quite satisfactory in all cases. Fluctuating bending moment and shear at any location along the wall height contributed by each individual vibration mode can be obtained. Modal contributions to shear and flexural strength demands, as well as the corresponding modal properties, under various seismic loading conditions can also be identified and examined in detail. Furthermore, the effects of higher vibration modes on seismic demands of wall buildings are investigated by using the modal decomposition methods. Several new insights into the complicated inelastic dynamics of multi‐story wall buildings are presented. Copyright © 2004 John Wiley & Sons, Ltd.     

MODAL ANALYSIS OF SOFT-SOIL SITES INCLUDING RADIATION DAMPING

For the one-dimensional analysis of soft-soil layers on an elastic half-space, a general form of analytical solution is developed for converting radiation damping due to energy leaking back to the half-space into equivalent modal damping, allowing the modal analysis technique to be extended to a site where radiation damping has to be accounted for. Closed-form solutions for equivalent modal damping ratios and effective modal participation factors are developed for a single layer with a shear wave velocity distribution varying from constant to linearly increasing with depth. Compact and recursive forms of solutions for equivalent modal damping ratios are developed for a system with an arbitrary number of homogeneous layers on an elastic half-space. Comparisons with numerical solutions show that the modal solutions are accurate. The nominal frequency of a site, i.e. the inverse of four times the total shear wave travel time through the layers, is an important parameter for estimating the high mode frequencies. A parameter study shows that for the same impedance ratio of the bottom layer to the elastic half-space, a system of soil layers with an increasing soil rigidity with depth has, in general, larger peak modal amplifications at the ground surface than does a single homogeneous layer on an elastic half-space, while a system with a decreasing soil rigidity with depth has smaller modal peak amplifications. © 1997 by John Wiley & Sons, Ltd.     

成层地基一维土层对地震的随机反应分析

首先基于改进的一维剪切梁模型，对成层土层推导了确定自振频率、振型函数、参与系数及稳态动力响应的封闭型解析表达式，首次证明了成层土层振型函数的正交性，然后在此基础上，利用随机振动理论，研究了成层土层对地震的随机动力反应问题，关于基岩输入地震加速度的功率谱密度函数，考虑了两种形式：白噪声谱和过滤白噪声谱。数值计算结果表明：对这两种谱，土层的最大期望反应是不相同的；平稳输入与输出过高估计了土层的随机反应。     

Seismic response analysis of multidrum classical columns

This paper presents a numerical investigation on the seismic response of multidrum classical columns. The motivation for this study originates from the need to understand: (a) the level of ground shaking that classical multidrum columns can survive, and (b) the possible advantages or disadvantages of retrofitting multidrum columns with metallic shear links that replace the wooden poles that were installed in ancient times. The numerical study presented in this paper is conducted with the commercially available software Working Model 2D?, which can capture with fidelity the sliding, rocking, and slide‐rocking response of rigid‐body assemblies. This paper validates the software Working Model by comparing selected computed responses with scarce analytical solutions and the results from in‐house numerical codes initially developed at the University of California, Berkeley, to study the seismic response of electrical transformers and heavy laboratory equipment. The study reveals that relative sliding between drums happens even when the g‐value of the ground acceleration is less than the coefficient of friction, µ, of the sliding interfaces and concludes that: (a) typical multidrum classical columns can survive the ground shaking from strong ground motions recorded near the causative faults of earthquakes with magnitudes M_w=6.0–7.4; (b) in most cases multidrum classical columns free to dislocate at the drum interfaces exhibit more controlled seismic response than the monolithic columns with same size and slenderness; (c) the shear strength of the wooden poles has a marginal effect on the sliding response of the drums; and (d) stiff metallic shear links in‐between column drums may have an undesirable role on the seismic stability of classical columns and should be avoided. Copyright © 2005 John Wiley & Sons, Ltd.     

Alternative method to locate centre of rigidity in asymmetric buildings

Response of asymmetric buildup under earthquake excitation often involves lateral vibration coupled with torsional vibration. Floor slab is, in general, assumed as rigid along the in‐plane direction. Building code provisions to account for the torsional effect in static force procedure are based on centre of rigidity or shear centre of the building. A convenient procedure is developed here to locate the centre of rigidity or shear centre, which can be implemented, using any standard building analysis software. The procedure is applicable for orthogonal as well as non‐orthogonal building systems and accounts for all possible definitions of static eccentricity to compute the design response. An irregular building is analysed to illustrate the proposed methodology. Significant variation in member force resultants is observed due to different definitions of static eccentricity. Finally, a mathematical proof is presented to substantiate the applicability of the proposed procedure to a non‐orthogonal building. Copyright © 2006 John Wiley & Sons, Ltd.     

Simplified seismic analysis of asymmetric building systems

The paper reviews the uncoupled modal response history analysis (UMRHA) and modal pushover analysis (MPA) procedure in the analysis of asymmetric structures. From the pushover curves in ADRS format, showing the relationships of base shear versus roof translation and base torque versus roof rotation, a bifurcating characteristic of the pushover curves of an asymmetric structure is observed. A two‐degree‐of‐freedom (2DOF) modal stick is constructed using lump mass eccentrically placed at the end of beam which is connected with the column by a rotational spring. By converting the equation of motion of a whole structure into 2DOF modal equations, all of the elastic properties in the 2DOF modal sticks can be determined accurately. A mathematical proof is carried out to demonstrate that the 2DOF modal stick is consistent with the single‐degree‐of‐freedom (SDOF) modal stick at elastic state. The bifurcating characteristic of modal pushover curves and the interaction of modal translation and rotation can be considered rationally by this 2DOF modal stick. In order to verify the effectiveness of this proposed 2DOF modal stick, a two‐storey asymmetric building structure was analysed by the UMRHA procedure incorporating this novel 2DOF modal sticks (2DMPA) and conventional SDOF modal sticks (SDMPA), respectively. The analytical results are compared with those obtained by nonlinear response history analysis (RHA). It is illustrated that the accuracy of the rotational response histories obtained by 2DMPA is much better than those obtained by SDMPA. Consequently, the estimations of translational response histories on flexible side (FS) and stiff side (SS) of the building structure are also improved. Copyright © 2006 John Wiley & Sons, Ltd.     

Effects of local site conditions on the seismic response of municipal solid waste landfills

Over the last decade, the seismic response of landfills made of municipal solid waste has drawn attention mainly due to the environmental and public-health issues that could be raised in the event of a failure. Nevertheless, there are several associated technical issues that have not been adequately investigated. One of these is the impact of local site conditions on the earthquake-induced accelerations and, thereby, on the seismic design of a landfill. This study presents the results of a parametric numerical simulation that has been performed in order to examine the effects of local site conditions on the dynamic response of a typical landfill. Emphasis is given on the special characteristics of ground motion, while the material nonlinearity of both soil and waste is taken into account by an equivalent-linear procedure. Results indicate that local site conditions may play a significant role in the seismic response of a landfill. However, this role cannot be judged a priori as beneficial or detrimental, as it depends not only on soil conditions and seismic excitation, but also on the material and geometric characteristics of the landfill.     

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.     

最小剪力系数及其调整方法对超高层建筑地震响应的影响

超高层结构地震剪力响应由振型分解反应谱法得到的结果经常不能满足规定的最小剪力系数要求。为此,文章简述剪力系数的概念和调整方法,以具有不同剪力系数的两个模型对比分析结构弹性、弹塑性地震响应差异,探讨剪力系数对超高层结构地震响应的影响。以通过强度和刚度调整使最小剪力系数满足规范要求的两个模型,分析不同调整方法引起的结构响应的合理性。结果表明:满足最小剪力系数的结构的弹性基底剪力大、层间位移角较小,结构的弹塑性位移响应也较小,受力状态优于不满足最小剪力系数的结构,安全性得到了提高。结构弹性倾覆力矩需求和弹塑性基底剪力按刚度调整大于按强度调整;结构弹塑性最大顶点位移和层间位移角响应相差不大,但出现刚度大\,层间位移角也大的与抗震理论相悖的情况;在满足抗震要求的情况下,构件的受力状态则是按强度调整更优,构件截面更加经济合理。     

An innovative methodology for seismic retrofitting of existing RC buildings by metal shear panels

In the present paper, the seismic upgrading of existing reinforced concrete (RC) structures by means of steel and pure aluminium shear panels is examined. After a preliminary experimental evaluation of the performance of the bare RC structure, a design approach based on the capacity spectrum method has been developed according to the procedure provided in the ATC 40 American guidelines. First, the geometrical configuration of the applied shear panels has been defined according to simplified analytical relationships, while appropriate steel members have been designed to allow the insertion of shear panels in the existing RC structure. Then, complex finite element models have been implemented in order to check the reliability of the proposed design procedure. Also, a numerical evaluation of the global response of the upgraded structure has been processed aiming at evaluating the interaction between the RC structure and the metal devices. Finally, the effectiveness of the applied shear panels has been proven by means of full‐scale experimental tests, which confirmed the significant improvement of the RC structure performance, in terms of strength, stiffness and deformation capacity. Copyright © 2008 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.     

挡土墙地震反应非线性波动模拟

本文运用解耦近场非线性波动数值模拟方法研究挡土墙地震反应，为反映墙土体系在地震作用下的位移机制，引入了Desai薄层单元模拟墙土间接触面，并采用双线型本构关系作为接触面单元和土体的非线性模型，在此基础上给出了解决P—SV问题的非线性显式有限元时域递推公式，为进一步发展非线性波动数值模拟技术提供了有益经验。为验证本文方法及适用性，将数值模拟结果与Zeng，X．和Madabhushi，X．P．G．等的离心机试验和弹塑性数值模拟结果进行对比。结果表明：墙土体系加速度、挡土墙顶底相对滑移、沉降和墙体倾角等同离心机试验模拟结果基本吻合，与弹塑性数值模拟结果相似。