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.     

Modeling viscous damping in nonlinear response history analysis of buildings for earthquake excitation

The Rayleigh damping model, which is pervasive in nonlinear response history analysis (RHA) of buildings, is shown to develop ‘spurious’ damping forces and lead to inaccurate response results. We prove that a viscous damping matrix constructed by superposition of modal damping matrices—irrespective of the number of modes included or values assigned to modal damping ratios—completely eliminates the ‘spurious’ damping forces. This is the damping model recommended for nonlinear RHA. Replacing the stiffness‐proportional part of Rayleigh damping by the tangent stiffness matrix is shown to improve response results. However, this model is not recommended because it lacks a physical basis and has conceptual implications that are troubling: hysteresis in damping force–velocity relationship and negative damping at large displacements. Furthermore, the model conflicts with the constant‐damping model that has been the basis for fundamental concepts and accumulated experience about the inelastic response of structures. With a distributed plasticity model, the structural response is not sensitive to the damping model; even the Rayleigh damping model leads to acceptable results. This perspective on damping provides yet another reason to employ the superior distributed plasticity models in nonlinear RHA. OpenSees software has been extended to include a damping matrix defined as the superposition of modal damping matrices. Although this model leads to a full populated damping matrix, the additional computational demands are demonstrated to be minimal. Copyright © 2015 John Wiley & Sons, Ltd.     

Problems encountered from the use (or misuse) of Rayleigh damping

Rayleigh damping is commonly used to provide a source of energy dissipation in analyses of structures responding to dynamic loads such as earthquake ground motions. In a finite element model, the Rayleigh damping matrix consists of a mass‐proportional part and a stiffness‐proportional part; the latter typically uses the initial linear stiffness matrix of the structure. Under certain conditions, for example, a non‐linear analysis with softening non‐linearity, the damping forces generated by such a matrix can become unrealistically large compared to the restoring forces, resulting in an analysis being unconservative. Potential problems are demonstrated in this paper through a series of examples. A remedy to these problems is proposed in which bounds are imposed on the damping forces. Copyright © 2005 John Wiley & Sons, Ltd.     

Discussion of ‘Modelling viscous damping in nonlinear response history analysis of buildings for earthquake excitation’ by Anil K. Chopra and Frank McKenna

This discussion deals with recommendations in the paper on appropriate damping formulations for use in nonlinear response history analysis of buildings. Concern over potentially excessive damping forces and moments should extend beyond the damping moments produced by the stiffness proportional part of Rayleigh damping that corresponds to rotational springs used to explicitly model plastic hinges. The key to an appropriate damping formulation for nonlinear analysis is a realistic mechanism that allows all damping forces and moments to be meaningfully assessed. Then features can be added to keep these forces and moments within reasonable bounds. Copyright © 2016 John Wiley & Sons, Ltd.     

多级串联非比例阻尼隔震结构地震响应分析

建立了多级串联非比例阻尼隔震结构动力分析模型;引用分区瑞利阻尼模型将非比例阻尼矩阵分解为瑞利阻尼矩阵和体现非比例阻尼的余项阻尼矩阵,推导出结构的阻尼矩阵;并编制了MATLAB动力时程分析程序,对一实际隔震工程进行地震响应分析.结果表明:随着下部结构刚度的增加,结构的层剪力比和隔震层位移响应峰值均趋向于基础隔震结构的对应值;当下部结构为一层,且层间刚度大于上部结构底层层间刚度4～6倍时,可以近似按基础隔震结构进行动力分析.     

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

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

Evaluation of Rayleigh damping and its influence on engineering demand parameter estimates

The effects of Rayleigh damping model on the engineering demand parameters of two steel moment‐resisting frame buildings were evaluated. Two‐dimensional models of the buildings were created and response history analysis were conducted for three different hazard levels. The response history analysis results indicate that mass‐proportional damping leads to high damping forces compared with restoring forces and may lead to overestimation of floor acceleration demands for both buildings. Stiffness‐proportional damping, on the other hand, is observed to suppress the higher‐mode effects in the nine‐story building resulting in lower story drift demands in the upper floors compared with other damping models. Rayleigh damping models, which combine mass‐proportional and stiffness‐proportional components, that are anchored at reduced modal frequencies lead to reasonable damping forces and floor acceleration demands for both buildings and does not suppress higher‐mode effects in the nine‐story building. Copyright © 2012 John Wiley & Sons, Ltd.     

Implementation and application of the unconditionally stable explicit parametrically dissipative KR‐α method for real‐time hybrid simulation

In real‐time hybrid simulations (RTHS) that utilize explicit integration algorithms, the inherent damping in the analytical substructure is generally defined using mass and initial stiffness proportional damping. This type of damping model is known to produce inaccurate results when the structure undergoes significant inelastic deformations. To alleviate the problem, a form of a nonproportional damping model often used in numerical simulations involving implicit integration algorithms can be considered. This type of damping model, however, when used with explicit integration algorithms can require a small time step to achieve the desired accuracy in an RTHS involving a structure with a large number of degrees of freedom. Restrictions on the minimum time step exist in an RTHS that are associated with the computational demand. Integrating the equations of motion for an RTHS with too large of a time step can result in spurious high‐frequency oscillations in the member forces for elements of the structural model that undergo inelastic deformations. The problem is circumvented by introducing the parametrically controllable numerical energy dissipation available in the recently developed unconditionally stable explicit KR‐α method. This paper reviews the formulation of the KR‐α method and presents an efficient implementation for RTHS. Using the method, RTHS of a three‐story 0.6‐scale prototype steel building with nonlinear elastomeric dampers are conducted with a ground motion scaled to the design basis and maximum considered earthquake hazard levels. The results show that controllable numerical energy dissipation can significantly eliminate spurious participation of higher modes and produce exceptional RTHS results. Copyright © 2014 John Wiley & Sons, Ltd.     

结构阻尼与材料阻尼的关系

针对固体材料在弹性阶段的内阻尼,采用材料粘弹性本构关系,通过有限单元法建立了考虑剪切变形影响的杆单元阻尼矩阵,得到了材料拉压粘滞系数和剪切粘滞系数之间的关系。我们可以通过实验测试材料的损耗因子或粘滞阻尼系数(材料阻尼),像结构质量矩阵、刚度矩阵一样通过比较确切的计算得到弹性阶段结构内阻尼的阻尼矩阵(结构阻尼),而且此阻尼矩阵既基于材料阻尼实验测试,又便于数学处理,且物理意义明确。其次,本文分析了材料阻尼对结构阻尼比影响,得到了材料损耗因子与结构模态阻尼比间的关系,并通过10层钢筋混凝土剪切型框架结构进一步给出了具体的数值结果。     

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.     

Effects of the improper modeling of viscous damping on the first-mode and higher-mode dominated responses of base-isolated buildings

In many finite element platforms, a classical global damping matrix based on the elastic stiffness of the system (including isolators) is usually developed as part of the solution to the equations of motion of base-isolated buildings. The conducted analytical and numerical investigations illustrate that this approach can lead to the introduction of unintended damping to the first and higher vibration modes and the spurious suppression of the respective structural responses. A similar shortcoming might be observed even when a nonclassical damping model (ie, an assembly of the superstructure and isolation system damping sub-matrices) is used. For example, the use of Rayleigh damping approach to develop the superstructure damping sub-matrix can lead to the undesired addition of damping to the isolated mode arising from the mass-proportional component of the superstructure damping. On the other hand, the improper use of nonclassical stiffness-proportional damping (eg, determining the proportional damping coefficient, β_k , based on the first mode) can result in assigning significant damping to the higher-modes and the unintended mitigation of the higher-mode responses. Results show that a nonclassical stiffness-proportional model in which β_k is determined based on the second modal period of a base-isolated building can reasonably specify the intended damping to the higher modes without imparting undesirable damping to the first mode. The nonclassical stiffness-proportional damping can be introduced to the numerical model through explicit viscous damper elements attached between adjacent floors. In structural analysis software such as SAP2000®, the desired nonclassical damping can be also modeled through specifying damping solely to the superstructure material.     

Constant-ductility inelastic displacement ratios for displacement-based seismic design of self-centering structures

This paper focuses on constant-ductility inelastic displacement ratios of self-centering single-degree-of-freedom (SDF) systems with two different levels of energy dissipation capacity, in the presence of 5% viscous damping ratio. A statistical analysis is developed considering an earthquake database composed of 228 ground motions recorded in California with magnitudes greater than six and organized for NEHRP soil class, ground motion duration, and peak ground acceleration. The response of self-centering SDF systems with large variability of initial periods, ductility levels, and postyield stiffness ratios is investigated and compared with the responses of SDF systems with bilinear plastic, Clough, and Takeda hysteresis. The inelastic demand variation with soil class, initial period, postyield stiffness ratio, unloading stiffness degradation, ductility level, and hysteretic behavior is highlighted. Simple and conservative analytical estimates of constant-ductility inelastic displacement ratios for mean and 90th percentile values in terms of initial period, ductility level, and postyield stiffness ratio are proposed to allow the extension of the Displacement-Based Design via Inelastic Displacement Ratio (C_μDBD) to self-centering structural systems.     

公路简支梁桥的地震能量响应分析

利用有限元软件SAP2000建立了某公路简支梁桥的有限元模型,以7条典型强震记录为输入,研究了公路简支梁桥的地震能量响应及其分配规律。结果表明:①地基柔性效应对公路简支梁桥的地震能量响应及其分配规律的影响较小;②当场地土质变软时,地震总输入能、结构阻尼耗能和结构阻尼耗能比均呈递增趋势,而结构滞回耗能和结构滞回耗能比则不断减小,即地基土体作为桥梁动力系统的一部分,增大了系统阻尼,并分担了部分非弹性变形;③随着PGA增大,输入结构的地震能量也增加,导致塑性铰的非弹性变形增加,即结构滞回耗能和结构阻尼耗能增大。     

Constant-ductility strength demand spectra for seismic design of structures

In displacement-based seismic design, constant-ductility strength demand spectra （CDSDS） are very useful for preliminary design of new structures where the global displacement ductility capacity is known. The CDSDS can provide the required inelastic lateral strength of new structures from the required elastic lateral strength. Based on a statistical study of nonlinear time-history for an SDOF system, the mean CDSDS corresponding to four site conditions are presented and approximate expressions of the inelastic spectra are proposed, which are functions of the structural period and ductility level. The effects of site conditions, structural period, level of ductility, damping and post-yield stiffness of structures on CDSDS are also investigated. It is concluded that site conditions, ductility level and structural period have important effects on the CDSDS and damping, post-yield stiffness effects are rather complex and of minor importance. The damping, post-yield stiffness effects depend on both the level of ductility and the natural period of structures.     

Practical causal hysteretic damping

A number of experiments indicate that the internal damping corresponding to the energy dissipation of many materials is essentially frequency independent. Accordingly, an analysis model that can express such characteristics (called a hysteretic damping model) in the time domain is needed. Although a great number of investigations into this subject have been carried out, there are a few practical methods. In this paper, a simple hysteretic damping model which satisfies the causality condition is presented using an extension of the complex stiffness transfer method that the author has proposed. Compared with the energy proportional damping model and the Biot model, the applicability and the efficiency of this model to time history response analyses were confirmed well by example problems. Copyright © 2006 John Wiley & Sons, Ltd.     

高速铁路客站模态特征分析与阻尼模型的注记

针对"房桥合一"高速铁路客站质量、刚度分布严重不均匀、不同阻尼比构件繁多等特点,对其模态特征和用于时程法阻尼模型的确定进行了研究。通过位能加权激励原理和位能公式的阐述、天津西站II区的模态分析与振型分解时程法的应用,并以此为标准进行了5种不同瑞利阻尼比例系数时程法的响应比较。结果表明:位能加权法激励为各模态的振型响应,求得阻尼为振型阻尼,结构模态质量累计数突变发生在第9阶,选择第1,9阶振型确定的瑞利阻尼比例系数较合理。振型分解时程法的振型阻尼可基于振型响应的位能加权法确定,直接时程法的瑞利阻尼宜选择第1阶与模态质量累计数发生突变的振型来确定,可供结构设计与分析参考。     

小应变硬化模型(HSS)在Rayleigh波作用下场地响应分析中的应用

为研究冲击荷载或地震作用下产生的,以Rayleigh波为主的面波对浅层地表土体动力响应特征以及数值模拟中土层阻尼的设置方法,以厦门地区浅层的素填土及粉质黏土为研究对象,采用有限元动力分析,土体本构采用小应变硬化模型(HSS),利用模型本身的滞回环特性,输入变化的小应变参数,考察HSS模型的小应变参数对场地动力响应的影响,并与土体采用摩尔-库伦模型结合Rayleigh阻尼("MC+Rayleigh阻尼")的计算结果进行对比。研究表明:当采用带有滞回环的HSS模型时,波速随初始剪切模量Gref0的增大而增大,但振幅减小,残余变形量也有所减小;小应变参数γ0.7对波的影响较小;HSS模型能够给出残余变形量,而"MC+Rayleigh阻尼"由于本构模型为理想弹塑性模型,在卸载重加载条件下表现为纯弹性行为,无法反映出卸载重加载过程中塑性应变的积累及其累积阻尼效应;但HSS模型还不能够全面反映循环加载作用下塑性体积应变的累积,因此在考虑滞回阻尼的基础上,仍然建议借助Rayleigh阻尼来更加全面地模拟土体的实际阻尼特性。     

Implementation of configuration dependent stiffness proportional damping for the dynamics of rigid multi-block systems

The distinct element method (DEM) has been used successfully for the dynamic analysis of rigid block systems. One of many difficulties associated with DEM is modeling of damping. In this paper, new procedures are proposed for the damping modeling and its numerical implementation in distinct element analysis of rigid multi-block systems. The stiffness proportional damping is constructed for the prescribed damping ratio, based on the non-zero fundamental frequency effective during the time interval while the boundary conditions remain essentially constant. At this time interval, the fundamental frequency can be estimated without complete eigenvalue analysis. The damping coefficients will vary while the damping ratio remains the same throughout the entire analysis. A new numerical procedure is developed to prevent unnecessary energy loss that can occur during the separation phases. These procedures were implemented in the development of the distinct element method for the dynamic analyses of piled multi-block systems. The analysis results for the single-block and two-block systems were in a good agreement with the analytic predictions. Applications to the seismic analyses of piled fourblock systems revealed that the new procedures can make a significant difference and may lead to much-improved results.     

基于楼层剪力的等效阻尼比计算方法研究

针对现有附加有效阻尼比计算方法存在的问题,本文从能量的角度揭示了阻尼比对结构影响的机理。从结构设计的角度,提出一种在时程分析下基于楼层剪力的消能减震结构等效阻尼比计算方法。对布置黏滞阻尼器和软钢阻尼器的消能减震模型,采用本文提出的等效阻尼比计算方法,建立等效结构进行结构响应对比。结果表明,由该计算方法得到的等效阻尼比能够准确地评估阻尼器在结构中的耗能效果,建立的等效结构能够准确反映消能减震结构实际情况。基于楼层剪力的等效阻尼比计算方法通过等效结构楼层剪力大于或等于消能减震结构楼层剪力判断迭代完成,该方法计算过程不涉及阻尼器参数及结构形式,适用于所有阻尼器类型与结构类型。计算得到的等效结构进行设计能够确保结构设计的安全。     

瑞利阻尼对砂土场地非线性地震响应的影响分析

在地震荷载作用下,自由场地会产生土体侧向变形和地表响应放大现象。由于土体的高度非线性,计算自由场地地震响应时,不同的阻尼比及剪切模量取值是造成其计算结果与试验结果相差较大的原因之一。目前动力计算常采用瑞利阻尼方法,其系数取值会在一定程度上影响计算结果。选用两模态简化瑞利阻尼系数计算方法,分析土体阻尼比及控制频率的取值对计算结果的影响,对比离心机模型试验,利用开源有限元平台OpenSees,采用适合于土体动力分析的多屈服面本构模型(PDMY),建立剪切梁模型模拟三维自由场地,并分析瑞利阻尼参数对自由场地地震响应和侧向变形计算结果的影响。结果表明,针对相对密度为60%的Nevada干砂,阻尼比为4%、控制频率比为5时,场地响应计算结果与试验结果较为符合。综合分析显示场地非线性响应时域计算时,应特别注意选用的瑞利阻尼参数值。