Assessment of expected damage on buildings subjected to Lorca earthquake through an energy-based seismic index method and nonlinear dynamic response analyses

The city of Lorca (Spain) was hit on May 11th, 2011, by two consecutive earth-quakes of magnitudes 4.6 and 5.2 Mw, causing casualties and important damage in buildings. Many of the damaged structures were reinforced concrete frames with wide beams. This study quantifies the expected level of damage on this structural type in the case of the Lorca earth-quake by means of a seismic index Iv that compares the energy input by the earthquake with the energy absorption/dissipation capacity of the structure. The prototype frames investigated represent structures designed in two time periods (1994–2002 and 2003–2008), in which the applicable codes were different. The influence of the masonry infill walls and the proneness of the frames to concentrate damage in a given story were further investigated through nonlinear dynamic response analyses. It is found that (1) the seismic index method predicts levels of damage that range from moderate/severe to complete collapse; this prediction is consistent with the observed damage; (2) the presence of masonry infill walls makes the structure very prone to damage concentration and reduces the overall seismic capacity of the building; and (3) a proper hierarchy of strength between beams and columns that guarantees the formation of a strong column-weak beam mechanism (as prescribed by seismic codes), as well as the adoption of counter-measures to avoid the negative interaction between non-structural infill walls and the main frame, would have reduced the level of damage from Iv=1 (collapse) to about Iv=0.5 (moderate/severe damage)

A comparison among different spectrum compatible earthquake simulation methods

The need for the simulation of spectrum compatible earthquake time histories has existed since earthquake engineering for complicated structures began. More than the safety of the main structure, the analysis of the equipment (piping, racks, etc.) can only be assessed on the basis of the time history of the floor in which they are contained. This paper presents several methods for calculating simulated spectrum compatible earthquakes as well as a comparison between them. As a result of this comparison, the use of the phase content in real earthquakes as proposed by Ohsaki appears as an effective alternative to the classical methods. With this method, it is possible to establish an approach without the arbitrary modulation commonly used in other methods. Different procedures are described as is the influence of the different parameters which appear in the analysis. Several numerical examples are also presented, and the effectiveness of Ohsaki's method is confirmed.

The direct boundary element method: 2D site effects assessment on laterally varying layered media (methodology)

The Direct Boundary Element Method (DBEM) is presented to solve the elastodynamic field equations in 2D, and a complete comprehensive implementation is given. The DBEM is a useful approach to obtain reliable numerical estimates of site effects on seismic ground motion due to irregular geological configurations, both of layering and topography. The method is based on the discretization of the classical Somigliana's elastodynamic representation equation which stems from the reciprocity theorem. This equation is given in terms of the Green's function which is the full-space harmonic steady-state fundamental solution. The formulation permits the treatment of viscoelastic media, therefore site models with intrinsic attenuation can be examined. By means of this approach, the calculation of 2D scattering of seismic waves, due to the incidence of P and SV waves on irregular topographical profiles is performed. Sites such as, canyons, mountains and valleys in irregular multilayered media are computed to test the technique. The obtained transfer functions show excellent agreement with already published results.

Simplified computational method for non-linear seismic analysis of bridges

After the experience gained during the past years it seems clear that nonlinear analysis of bridges are very important to compute ductility demands and to localize potential hinges. This is specially true for irregular bridges in which it is not clear weather or not it is possible to use a linear computation followed by a correction using a behaviour factor. To simplify the numerical effort several approximate methods have been proposed. Among them, the so-called Dynamic Plastic Hinge Method in which an evolutionary shape function is used to reduce the structure to a single degree of freedom system seems to mantein a good balance between accuracy and simplicity. This paper presents results obtained in a parametric study conducted under the auspicies of PREC-8 european research program.

A scaled-time test method

The Pseudo-Dynamic Test Method (PDTM) is being developped currently as an alternative to the shaking table testing of large size models. However, the stepped slow execution of the former type of test has been found to be the source of important errors arising from the stress relaxation. A new continuous test method, wich allows the selection of a suitable time-scale factor in the response in order to control these errors, es proposed here. Such scaled-time response is theoretically obtained by simply augmenting the mass of the structure for wich some practical solutions are proposed.

Estudio probabilista del daño por choque según la metodología PEER. Aplicación a un paso superior sísmicamente aislado

En países de sismicidad moderada es habitual el uso de puentes cuyo tablero está unido mediante apoyos elastoméricos y/o deslizantes, tanto a las pilas como a los estribos. La principal vulnerabilidad sísmica de estos puentes está asociada a los choques. Dicha vulnerabilidad puede valorarse estimando la probabilidad de que se produzca un determinado daño en un periodo de tiempo. En este artículo se plantean las líneas actuales para estudiar este problema, y se aplican al estudio de la vulnerabilidad de un paso superior situado en la provincia de Granada, que puede considerarse como sísmicamente aislado. El marco del estudio es la metodología del Pacific Earthquake Engineering Research Center, cuya aplicación se ha adaptado al caso estudiado. En el artículo se describen cada una de las etapas de que consta, aunque se presta especial atención a la identificación y caracterización de los daños susceptibles de producirse. El modelo de la estructura, que debe permitir una adecuada caracterización de los choques, también será presentado en detalle. Dado que el puente es representativo de una gran cantidad de estructuras del inventario construido en España, los resultados son doblemente importantes. Por una parte, sirven de guía para la aplicación de esta metodología a estructuras similares. Y por otra, proporcionan valores de probabilidades que pueden servir de referencia.

Dynamic impedances of bridge abutments

The soil-structure interaction at bridge abutments may introduce important changes in the dynamic properties of short to medium span bridges. The paper presents the results obtained, through the use of the Boundary Element Method (B.E.M.) technique in several typical situations, including semiinfinite and layered media. Both stiffness and damping properties are included.

Shaking table tests of structures with hysteretic dampers: experimental results vs. prediction using non-linear static methods

The use of seismic hysteretic dampers for passive control is increasing exponentially in recent years for both new and existing buildings. In order to utilize hysteretic dampers within a structural system, it is of paramount importance to have simplified design procedures based upon knowledge gained from theoretical studies and validated with experimental results. Non-linear Static Procedures (NSPs) are presented as an alternative to the force-based methods more common nowadays. The application of NSPs to conventional structures has been well established; yet there is a lack of experimental information on how NSPs apply to systems with hysteretic dampers. In this research, several shaking table tests were conducted on two single bay and single story 1:2 scale structures with and without hysteretic dampers. The maximum response of the structure with dampers in terms of lateral displacement and base shear obtained from the tests was compared with the prediction provided by three well-known NSPs: (1) the improved version of the Capacity Spectrum Method (CSM) from FEMA 440; (2) the improved version of the Displacement Coefficient Method (DCM) from FEMA 440; and (3) the N2 Method implemented in Eurocode 8. In general, the improved version of the DCM and N2 methods are found to provide acceptable accuracy in prediction, but the CSM tends to underestimate the response.

Maximum Likelihood Estimation of Modal Parameters in Structures Using the Expectation Maximization Algorithm

This paper presents a time-domain stochastic system identification method based on maximum likelihood estimation (MLE) with the expectation maximization (EM) algorithm. The effectiveness of this structural identification method is evaluated through numerical simulation in the context of the ASCE benchmark problem on structural health monitoring. The benchmark structure is a four-story, two-bay by two-bay steel-frame scale model structure built in the Earthquake Engineering Research Laboratory at the University of British Columbia, Canada. This paper focuses on Phase I of the analytical benchmark studies. A MATLAB-based finite element analysis code obtained from the IASC-ASCE SHM Task Group web site is used to calculate the dynamic response of the prototype structure. A number of 100 simulations have been made using this MATLAB-based finite element analysis code in order to evaluate the proposed identification method. There are several techniques to realize system identification. In this work, stochastic subspace identification (SSI)method has been used for comparison. SSI identification method is a well known method and computes accurate estimates of the modal parameters. The principles of the SSI identification method has been introduced in the paper and next the proposed MLE with EM algorithm has been explained in detail. The advantages of the proposed structural identification method can be summarized as follows: (i) the method is based on maximum likelihood, that implies minimum variance estimates; (ii) EM is a computational simpler estimation procedure than other optimization algorithms; (iii) estimate more parameters than SSI, and these estimates are accurate. On the contrary, the main disadvantages of the method are: (i) EM algorithm is an iterative procedure and it consumes time until convergence is reached; and (ii) this method needs starting values for the parameters. Modal parameters (eigenfrequencies, damping ratios and mode shapes) of the benchmark structure have been estimated using both the SSI method and the proposed MLE + EM method. The numerical results show that the proposed method identifies eigenfrequencies, damping ratios and mode shapes reasonably well even in the presence of 10% measurement noises. These modal parameters are more accurate than the SSI estimated modal parameters.

The scaled-time test as an alternative to the pseudo-dynamic test

The stepped and excessively slow execution of pseudo-dynamic tests has been found to be the source of some errors arising from strain-rate effect and stress relaxation. In order to control those errors, a new continuous test method which allows the selection of a more suitable time scale factor in the response is proposed in this work. By dimensional analysis, such scaled-time response is obtained theoretically by augmenting the inertial and damping properties of the structure, for which we propose the use of hydraulic pistons which are servo-controlled to produce active mass and damping, nevertheless using an equipment which is similar to that required in a pseudo-dynamic test. The results of the successful implementation of this technique for a simple specimen are shown here.

Damage assessment on building structures subjected to the recent near-fault earthquake in Lorca (Spain)

The city of Lorca (Spain) was hit on May 11th 2011 by two consecutive earthquakes with 4.6 and 5.2 Mw respectively, causing casualties and important damage in buildings. Lorca is located in the south-east region of Spain and settled on the trace of the Murcia-Totana-Lorca fault. Although the magnitudes of these ground motions were not severe, the damage observed was considerable over a great amount of buildings. More than 300 of them have been demolished and many others are being retrofitted. This paper reports a field study on the damage caused by these earthquakes. The observed damage is related with the structural typology. Further, prototypes of the damaged buildings are idealized with nonlinear numerical models and their seismic behavior and proneness to damage concentration is further investigated through dynamic response analyses.

Elastodynamics

In this chapter some applications of boundary element techniques to dynamic problems are presented. First, the basic theory is briefly reviewed in order to provide the necessary background to interpret the numerical results (for a fuller account of elastodynamic theory we recommend a study of the specialized literature). The second part of the chapter is devoted to the numerical implementation of the BEM. The presentation is based on the steady-state solution because this is the area in which most experience exists. This is by no means a limitation of the BEM method, and the use of integral transformations to obtain transient solutions is a well established procedure. Finally, in the third part three examples are presented. The first example is the steady-state solution of a plate under cyclic forces with and without a crack. The second example relies on the determination of soil compliances necessary to study soil-structure interaction and the third example treats the problem of the influence of different incidence angles of incoming waves in foundations. The last two examples are relevant to earthquake engineering problems for which the BEM is very well suited.

Performance of masonry buildings during the 2011 Lorca earthquake

On Wednesday 11th May 2011 at 6:47 pm (local time) a magnitude 5.1 Mw earthquake occurred 6 km northeast of Lorca with a depth of around 2 km. As a consequence of the shallow depth and the small epicentral distance, important damage was produced in several masonry constructions and even led to the collapse of some of them. Pieces of the facades of several buildings fell down onto the sidewalk, being one of the reasons for the killing of a total of 9 people. The objective of this paper is to describe and analyze the failure patterns observed in unreinforced masonry buildings ranging from 3 to 8 floors in height. First, a brief description of the local building practices of masonry buildings is given. Then, the most important failure types of masonry buildings are described and discussed. After that, a more detailed analysis of one particular building is presented.

Performance of buildings with masonry infill walls during the 2011 Lorca earthquake

On Wednesday 11th May 2011 at 6:47 pm (local time) a magnitude 5.1 Mw earthquake occurred 6 km northeast of Lorca with a depth of around 5 km. As a consequence of the shallow depth and the small epicentral distance, important damage was produced in several masonry constructions and even led to the collapse of one of them. Pieces of the facades of several buildings fell down onto the sidewalk, being one of the reasons for the killing of a total of 9 people. The objective of this paper is to describe and analyze the failure patterns observed in reinforced concrete frame buildings with masonry infill walls ranging from 3 to 8 floors in height. Structural as well as non-structural masonry walls suffered important damage that led to redistributions of forces causing in some cases the failure of columns. The importance of the interaction between the structural frames and the infill panels is analyzed by means of non-linear Finite Element Models. The resulting load levels are compared with the member capacities and the changes of the mechanical properties during the seismic event are described and discussed. In the light of the results obtained the observed failure patterns are explained. Some comments are stated concerning the adequacy of the numerical models that are usually used during the design phase for the seismic analysis.

Foreword special issue Lorca’s earthquake

(This is an excerpt from the content) On May 11 2011 at 1705 hours, a small 4.5 Mw. magnitude earthquake struck the town of Lorca in south-eastern Spain. Other than alarmed citizens, only minor damage to buildings occurred due to this quake. Unfortunately at 1847 hours, a second shock registering a magnitude of 5.1 Mw. and very shallow (just around 2 km under the city) produced the largest seismic catastrophe registered in Spain in the last 120 years. This second shock is commonly referred to as “Lorca’s earthquake” and the following papers describe the context, circumstances and consequences of the event. Spain is a country of moderate seismic hazard in a global context. Before the Lorca earthquake, the most destructive earthquake in modern times was the so-called “Andalusian earthquake” (25th December 1884) that resulted in 750 fatalities and more than 1,500 injuries, reaching X in Mercalli’s intensity scale. Despite the lack of catastrophic events in the last 120 years, Spain has always had a scientific interest in seismic ...