Assessment of the dynamic response of unreinforced masonry structures using a macroelement modeling approach

The seismic performance of unreinforced masonry structures is strongly associated with the interaction between in‐plane and out‐of‐plane mechanisms. The seismic response of these structures has been thoroughly investigated by means of experimental testing, analytical procedures, and computational approaches. Within the framework of the numerical simulations, models based on the finite element method provide a good prediction of the seismic performance of unreinforced masonry structures. However, they usually require a high computational cost and advanced user expertise to define appropriate mechanical properties and to interpret the numerical results. Because of these limitations, simplified models for practical applications have been developed during the last decades. Despite this, a great number of these models focus mostly on the evaluation of the in‐plane response, assuming box (or integral) behavior of the structure. In this paper, a simplified macroelement modeling approach is used to simulate the seismic response of 2 masonry prototypes taking into consideration the combined in‐plane and out‐of‐plane action. The numerical investigations were performed in the static and dynamic fields by using pushover analyses and nonlinear dynamic analyses respectively. The latter is a novel implementation of a model previously developed for static analysis. The results obtained from this study are in good agreement with those provided by a detailed nonlinear continuum FE approach, demonstrating the applicability of this macroelement model with a significant reduction of the computational cost.     

地铁隧道地震反应非线性分析

本文对动力有限元法中的一系列问题进行了讨论,对土体动力本构模型和土-结构接触面模型的设定提出了建议．在此基础上以天津地铁1、2号线工程为实例,分别对水平地震作用下圆形盾构隧道的动位移、动应力与动应变以及地表位移等动力反应进行了非线性分析,得出了一些有益的规律和结论,从而为地铁隧道的抗震设计与研究提供参考和依据。     

H形断面钢筋混凝土立体剪力墙的非线性动力响应分析

为了建立统一的动力响应分析模型,本文以NUPEC振动台试验的H形断面钢筋混凝土立体剪力墙为研究对象进行了三维非线性有限元动力响应分析。根据分析结果与试验结果的比较可知,在RC剪力墙到达最大承载力之前由简化模型和一般模型得到的动力响应特性与试验结果吻合较好,荷载-变形关系能很好模拟试验结果。但是,最大承载力之后,由于混凝土开裂、损伤、劣化的急剧发展,较难模拟混凝土开裂、裂缝的开闭及滑移等非线性特性,分析得到的加速度衰减较慢、位移响应较小。基于上述研究成果探讨并提出了进一步改善非线性有限元动力响应分析精度的建议。     

基于重力卫星数据监测地表质量变化的三维点质量模型法

利用卫星重力测量手段监测全球质量变化取得了巨大成功,本文基于牛顿万有引力定律在三维空间直角坐标系中导出利用重力卫星观测数据监测全球质量变化的三维点质量模型法,该方法可直接利用重力卫星的轨道和星间观测数据或时变重力场模型计算全球质量变化,由于利用卫星观测数据计算地表质量变化的向下延拓过程以及观测数据噪声的影响,需要采用合适的空间约束方程或正则化技术对解算结果进行约束或平滑处理.利用合成全球质量变化模型模拟一个月的GRACE双星轨道和星间距离变率数据计算全球质量变化,对三维点质量模型法进行分析验证,采用零阶Tikhonov正则化技术处理病态问题.结果表明,三维点质量模型法可有效用于重力卫星观测数据监测全球质量变化,为利用重力卫星观测数据监测全球质量变化提供一种可选的途径.

     

大型渡槽有限条法动力建模研究

根据渡槽结构特点，采用有限条划分渡槽槽身，推导出渡槽槽身结构的单元刚度矩阵、相容质量矩阵的显式表达式；渡槽支架采用空间梁单元模拟；联接槽身和支架的盆式橡胶支座采用弹性元件单元模拟。编程程序具体计算了某大型渡槽的模态，并用ANSYS对大型渡槽进行模态分析，计算结果表明：该模型计算的渡槽结构的固有频率和ANSYS计算结果接近。模型可用于大型渡槽的动力分析。     

中高层钢筋混凝土异形柱框架结构弹塑性时程分析

用杆模型编制了钢筋混凝土异形柱框架结构的弹塑性时程分析程序FEANT。采用该程序对一中高层大开间钢筋混凝土异形柱框架结构振动台试验模型进行了计算分析，计算值与试验值吻合较好，能满足工程需要。     

Nonlinear surface ground analysis via statistical approach

A statistical approach is proposed for nonlinear surface ground analysis. In contrast to the conventional method which deals with only a single ground motion for equivalent linearization of soil properties, a design response spectrum defined at the upper level (bottom of the surface ground) of an engineering bedrock can be handled as the target design earthquake in the present paper. The effective shear strain in each soil layer is evaluated by means of a statistical procedure in which the mean peak shear strain is computed in terms of its standard deviation and the corresponding peak factor. The stiffness and damping ratio of each soil layer are obtained iteratively from the nonlinear relation of stiffness reduction factors and damping ratios with respect to the strain level. After the evaluation of the equivalent stiffness and damping ratio of every soil layer, the ground surface response spectrum is transformed from the design response spectrum defined at the upper level of the engineering bedrock via the one-dimensional wave propagation theory. The reliability and accuracy of the proposed analysis method is examined through the comparison with the results by the conventional method (represented by the program) for many simulated spectrum-compatible ground motions.     

Evaluation of substructuring method for seismic soil-structure interaction analysis of bridges

This paper evaluates the commonly used substructuring method for analysis of bridge systems where the bridge is divided into two sub-systems: the bridge superstructure and the substructure including the pile foundations, abutments, and soil. Modeling of the soil-structure interaction (SSI) in the system is simplified by replacing the pile foundations, abutments, and soil with sets of independent equivalent linear springs and dashpots at the base of the superstructure. The main objective of the paper is to examine how well the substructuring method simulates the seismic response of a bridge system. The baseline data required for the evaluation process is derived from analyzing a fully-coupled continuum bridge model, already validated for the instrumented two-span Meloland Road Overpass. The same bridge system is also simulated using the substructuring method. The results from both approaches are compared, and it is shown that the differences between them can be significant. The substructuring method consistently overestimates the pier base shear forces and bending moments and the pier top deflections. Moreover, the spectral response of the bridge structure is mispredicted. The analyses are repeated for a three-span bridge system subjected to several ground motions, leading to a similar observation as before. Hence, the current state of practice for simulating seismic SSI in bridges using the substructure model is shown to be too simplified to capture the major mechanisms involved in SSI.     

A method for nonlinear dynamic response analysis of semi-infinite foundation soil

This paper presents development of a special finite difference method for the nonlinear dynamic response analysis of semi-infinite foundation soil. Semi-infinite domain is mapped into the finite domain using special mapping. For the region of engineering interest, mapping is isometrical, and for far field, shrink mapping which transforms an infinite interval into a finite interval is adopted. Using linear and nonlinear constitutive models, the responses of semi-infinite foundation soil are computed using a proposed method with a small mesh model and an extensive mesh model. Surface loadings or incident earthquake waves are applied to the models in the computations. Good agreements were obtained among the theoretical and computed results of the two models and the effectiveness of the proposed method was demonstrated.     

Nonlinear analysis of masonry structures using fiber-section line elements

For seismic analysis of unreinforced masonry (URM) buildings characterized by a box-like behavior, a widely accepted model is based on the equivalent frame idealization of walls. The equivalent frame model uses 1D elements to represent the vertical piers and horizontal spandrels which are connected by rigid nodes. The mechanical characterization of the elements is one of the crucial aspects to predict reasonably the building seismic behavior. Through the comparison with pseudo-static and dynamic experimental tests performed on two-story full-scale buildings, this paper validates the frame modeling in the OpenSees framework, which includes a fiber-section force-based beam element for the axial-flexural behavior, coupled with a cyclic shear-deformation phenomenological law.     

Incremental dynamic analysis with consideration of modeling uncertainties

Incremental dynamic analysis (IDA) has been extended by introducing a set of structural models in addition to the set of ground motion records which is employed in IDA analysis in order to capture record‐to‐record variability. The set of structural models reflects epistemic (modeling) uncertainties, and is determined by utilizing the latin hypercube sampling (LHS) method. The effects of both aleatory and epistemic uncertainty on seismic response parameters are therefore considered in extended IDA analysis. The proposed method has been applied to an example of the four‐storey‐reinforced concrete frame, for which pseudo‐dynamic tests were performed at the ELSA Laboratory, Ispra. The influence of epistemic uncertainty on the seismic response parameters is presented in terms of summarized IDA curves and dispersion measures. The results of extended IDA analysis are compared with the results of IDA analysis, and the sensitivity of the seismic response parameters to the input random variable using the LHS method is discussed. It is shown that epistemic uncertainty does not have significant influence on the seismic response parameters in the range far from collapse, but could have a significant influence on collapse capacity. Copyright © 2008 John Wiley & Sons, Ltd.     

土石坝地震响应的非线性剪切条模型与对比分析

对一维剪切条计算模型进行改进,提出了土石坝非线性地震反应的简化计算方法。首先将坝体沿坝高离散为一系列的具有不同剪切刚度与阻尼比等参数特性的层状体系,建立了各层的振动控制方程及其边值条件,进而采用数学物理方程方法进行了求解,确定了体系的振动特性,并根据振型叠加原理和Duhamel积分确定了坝体地震反应的线弹性解。采用等价线性化方法考虑坝料的动力非线性性质,通过对线弹性地震响应的反复迭代计算,使得各层土的模量和阻尼比与其相应的剪应变水平相协调,确定出与非线性坝体系统相等效的线性解答,并将所得到的地震响应作为非线性地震响应的近似解。最后,以均质坝和心墙坝作为算例进行了具体的数值计算,将所得结果与有限元数值解进行对比分析,论证了所提方法的适用性和合理性。     

Approximate incremental dynamic analysis using the modal pushover analysis procedure

Incremental dynamic analysis (IDA)—a procedure developed for accurate estimation of seismic demand and capacity of structures—requires non‐linear response history analysis of the structure for an ensemble of ground motions, each scaled to many intensity levels, selected to cover the entire range of structural response—all the way from elastic behaviour to global dynamic instability. Recognizing that IDA of practical structures is computationally extremely demanding, an approximate procedure based on the modal pushover analysis procedure is developed. Presented are the IDA curves and limit state capacities for the SAC‐Los Angeles 3‐, 9‐, and 20‐storey buildings computed by the exact and approximate procedures for an ensemble of 20 ground motions. These results demonstrate that the MPA‐based approximate procedure reduces the computational effort by a factor of 30 (for the 9‐storey building), at the same time providing results to a useful degree of accuracy over the entire range of responses—all the way from elastic behaviour to global dynamic instability—provided a proper hysteretic model is selected for modal SDF systems. The accuracy of the approximate procedure does not deteriorate for 9‐ and 20‐storey buildings, although their dynamics is more complex, involving several ‘modes’ of vibration. For all three buildings, the accuracy of the MPA‐based approximate procedure is also satisfactory for estimating the structural capacities for the limit states of immediate occupancy, collapse prevention, and global dynamic instability. Copyright © 2006 John Wiley & Sons, Ltd.     

城市地下变电站非线性地震反应分析

为研究城市地下变电站结构及其内部电气设备的动力特性及地震反应规律,以西安市快速轨道交通二号线张家堡地下变电站为工程背景,依据黏弹性人工边界理论及有限单元法,建立考虑土体-结构-设备动力相互作用的三维有限元计算模型,首先通过半空间自由场模型验证了地震动输入方法的有效性,继而开展地下变电站系统的模态分析和动力时程分析,得到了该类结构的动力特性及地震响应规律。研究结果表明,结构上覆土体及电气设备在地震作用下产生的惯性力是结构构件应力增长的主要原因,8度罕遇地震作用下地下变电站内电气设备的加速度放大系数明显超出规范要求,需要采取相应保护措施。     

结构动力反应分析的三阶显式方法

本通过对传统动力反应分析方法的总结，阐明了建立隐式和显式方法的一般思路及数学本质，提出了使用系统位移反向向量三阶导数的隐工和实用显式积分方法－3阶显式方法，分析了该显式方法的精度和稳定性，并对建立更高阶隐式和显式方法以及方法的精度和稳定性作了初步讨论。最后，通过算例对本方法、献[1]方法和经典的常平均加速度法（隐式方法、视为精确解）的精度和稳定性进行了比较分析。结果表明，本方法具有明显的优点。     

Physical model for dynamic analysis of structures with FPS isolators

This article presents a physical model for frictional pendulum isolators (FPS) that is ready to be implemented in most commercial software. The model is capable of accounting for effects such as large deformations, sticking, and uplift and impact by sensing the normal loads in the isolators through a gap element. Sticking has been incorporated into the model by extending the Park–Wen hysteretic model to the case of large deformations. The proposed model has been tested against a theoretically ‘exact’ formulation leading to essentially identical results. To facilitate its use, the physical FPS model has been cast into a typical non‐linear structural element format, i.e. with deformation as input and restoring force as output. Examples of a building and a bridge have been chosen to show the potential of the element and to provide further insight into the earthquake response of structures with FPS isolators; in particular, in aspects such as the orientation in placement of the isolator, sticking, P? Δ, and other large deformation effects. Copyright © 2003 John Wiley & Sons, Ltd.     

单向张弦梁结构风振响应的时域分析

张弦梁屋盖结构自重轻,地震作用下的动力问题不是结构设计的控制因素,但是在风荷载作用下结构的动力响应较为显著,在设计中应给予足够重视。本文首先利用Matlab编程工具,模拟结构的风速时程曲线,然后利用有限元软件中的瞬态动力学分析模块对这种结构进行动力时程分析,讨论其在脉动风荷载作用下的动力响应,并将结果和频域法脉动风荷载的分析结果进行对比,验证了计算结构的可靠性。得到一些有意义的结论,供张弦梁结构的工程应用参考。     

Post‐quake decision analysis using dynamic programming

We introduce a general decision analysis procedure based on stochastic dynamic programming in the post‐quake aftershock environment. The damage sustained by the building due to the mainsheet, the time‐varying aftershock rates and the potential for further damage progression in the post‐quake environment are all factors taken into consideration in the proposed methodology. This procedure enables the optimal decision after the mainshock to be selected based on the minimization of expected financial losses, subject to a constraint on a minimal level of individual life‐safety, using a consistent probabilistic framework to explicitly quantify the uncertainties in the variables. Copyright © 2008 John Wiley & Sons, Ltd.     

闸门-地基非线性动力相互作用分析

采用大型有限元软件对地基、闸门进行整体建模,讨论了在动力相互作用的分析中,简单截断人工边界的合理取值范围以及土体非线性的塑性性质对分析结果的影响;比较了在深层搅拌桩加固的复合地基中,按桩、土分别独立建模与按常规的平均法进行分析的差异。分析结果表明,不同的分析模型,人工边界的合理取值范围是不一样的;土体的滤波作用及其塑性的非线性性质对分析结果有重大的影响;对深层搅拌桩加固的复合地基,应该按实际把桩、土分别独立建模进行分析。分析结果对工程实际有一定的指导意义。     

不同设防标准RC框架结构基于易损性分析的抗震性能评估

借助非线性动力时程分析,对严格按照规范Ⅶ、Ⅷ、Ⅸ度设计的5个三跨6层钢筋混凝土框架结构开展易损性分析,建立了基于峰值加速度的易损性曲线。从易损性的角度对不同设防标准RC框架结构的抗震性能做了定量评价,并探讨了设防标准对RC框架结构易损性的影响。分析表明,对应于设防小震、中震及大震水平的峰值加速度,结构“小震不坏”、“中震可修”和“大震不倒”的失效概率均在18%以内,可认为结构满足三水准的性态控制目标。随着结构设防标准的提高,其易损性随之降低,相同峰值加速度对应的各个破坏状态的超越概率均有所降低。此外,将框架结构的设防烈度提高1度,其“大震不倒”的失效概率会明显减小。而将框架结构的设防烈度降低1度,其“大震不倒”的失效概率会显著增加,最高可达4倍。