Direct estimation of the seismic demand and capacity of oscillators with multi‐linear static pushovers through IDA

SPO2IDA is introduced, a software tool that is capable of recreating the seismic behaviour of oscillators with complex quadrilinear backbones. It provides a direct connection between the static pushover (SPO) curve and the results of incremental dynamic analysis (IDA), a computer‐intensive procedure that offers thorough demand and capacity prediction capability by using a series of nonlinear dynamic analyses under a suitably scaled suite of ground motion records. To achieve this, the seismic behaviour of numerous single‐degree‐of‐freedom (SDOF) systems is investigated through IDA. The oscillators have a wide range of periods and feature pinching hysteresis with backbones ranging from simple bilinear to complex quadrilinear with an elastic, a hardening and a negative‐stiffness segment plus a final residual plateau that terminates with a drop to zero strength. An efficient method is introduced to treat the backbone shape by summarizing the analysis results into the 16, 50 and 84% fractile IDA curves, reducing them to a few shape parameters and finding simpler backbones that reproduce the IDA curves of complex ones. Thus, vast economies are realized while important intuition is gained on the role of the backbone shape to the seismic performance. The final product is SPO2IDA, an accurate, spreadsheet‐level tool for performance‐based earthquake engineering that can rapidly estimate demands and limit‐state capacities, strength reduction R‐factors and inelastic displacement ratios for any SDOF system with such a quadrilinear SPO curve. Copyright © 2006 John Wiley & Sons, Ltd.     

Imaging Pathways in Fractured Rock Using Three‐Dimensional Electrical Resistivity Tomography

Major challenges exist in delineating bedrock fracture zones because these cause abrupt changes in geological and hydrogeological properties over small distances. Borehole observations cannot sufficiently capture heterogeneity in these systems. Geophysical techniques offer the potential to image properties and processes in between boreholes. We used three‐dimensional cross borehole electrical resistivity tomography (ERT) in a 9 m (diameter) × 15 m well field to capture high‐resolution flow and transport processes in a fractured mudstone contaminated by chlorinated solvents, primarily trichloroethylene. Conductive (sodium bromide) and resistive (deionized water) injections were monitored in seven boreholes. Electrode arrays with isolation packers and fluid sampling ports were designed to enable acquisition of ERT measurements during pulsed tracer injections. Fracture zone locations and hydraulic pathways inferred from hydraulic head drawdown data were compared with electrical conductivity distributions from ERT measurements. Static ERT imaging has limited resolution to decipher individual fractures; however, these images showed alternating conductive and resistive zones, consistent with alternating laminated and massive mudstone units at the site. Tracer evolution and migration was clearly revealed in time‐lapse ERT images and supported by in situ borehole vertical apparent conductivity profiles collected during the pulsed tracer test. While water samples provided important local information at the extraction borehole, ERT delineated tracer migration over spatial scales capturing the primary hydrogeological heterogeneity controlling flow and transport. The fate of these tracer injections at this scale could not have been quantified using borehole logging and/or borehole sampling methods alone.     

Historical seismicity, palaeoseismicity and seismic risk in Western Macedonia, Northern Greece

Western Macedonia, Northern Greece, was a seismically quiescent region for one or more centuries, and was regarded as a nearly aseismic, rigid block inside a broad zone of distributed continental deformation and faulting, and a region of minimum seismic risk. Consequently, the May 13, 1995 destructive earthquake (M = 6.6) which hit this assumed aseismic zone was a surprise for scientists, government and population.However, historical and archaeoseismic evidence, as well as coastal change data indicate that the assumed aseismic region of Western Macedonia has been affected in the last 2,000 years by at least seven, and possibly nine destructive earthquakes. One of these earthquakes occurred in circa 1700, and probably had the same epicentre with, but higher magnitude than the 1995 shock.The earthquake in circa 1700 is deduced from historical data and is modelled on the base of a swarm of church repairs which is explained as post-seismic recovery of the broader Kozani area: except for certain well known cases of towns or areas in which religious privileges were granted, large scale repairs or reconstruction of churches during the Ottoman period were possible only after Sultan's permissions, usually following earthquakes and other calamities.It can hence be concluded that some, at least, of the apparently aseismic regions inside broad zones of distributed seismicity are hit by stronger shocks, but with longer (200 years or more) recurrence intervals than their adjacent zones. Consequently, the seismic risk of the apparently aseismic regions is certainly not low, especially since relatively long periods of seismic quiescence lead to constructions vulnerable to earthquakes.     

Buckling and vibration analysis of flexible beams resting on an elastic half-space

The present paper deals with the buckling and vibration of prismatic beams resting on an elastic half-space. The computational procedure developed herein utilizes the advantages of both an analytical approach and a finite element scheme. This is accomplished by deriving exact frequency and axial force dependent stiffness matrices using the analytical solutions of the governing differential equation as ‘shape functions’. The major advantages of the proposed approach over previous ones are pointed out and discussed in detail. Numerical results demonstrating the performance of the proposed method are presented in the final part of the paper.     

Modelling of R/C members accounting for shear failure localisation: Finite element model and verification

Reinforced concrete (R/C) frame buildings designed according to older seismic codes represent a large part of the existing building stock worldwide. Their structural elements are often vulnerable to shear or flexure‐shear failure, which can eventually lead to loss of axial load resistance of vertical elements and initiate vertical progressive collapse of a building. In this study, a computationally efficient member‐type finite element model for the hysteretic response of shear critical R/C frame elements up to the onset of axial failure is presented; it accounts for shear‐flexure interaction and considers, for the first time, the localisation of shear strains, after the onset of shear failure, in a critical length defined by the diagonal failure plane. Its predictive capabilities are verified against experimental results of column and frame specimens and are shown to be accurate not only in terms of total response, but also with regard to individual deformation components. The accuracy, versatility, and simplicity of this finite element model make it a valuable tool in seismic analysis of complex R/C buildings with shear deficient structural elements.     

Modelling of R/C members accounting for shear failure localisation: Hysteretic shear model

Reinforced concrete (R/C) frame buildings designed according to older seismic codes represent a large part of the existing building stock worldwide. Their structural elements are often vulnerable to shear or flexure‐shear failure, which can eventually lead to loss of axial load resistance of vertical elements and initiate vertical progressive collapse of a building. In this study, a hysteretic model capturing the local shear response of shear‐deficient R/C elements is described in detail, with emphasis on post‐peak behaviour; it differs from existing models in that it considers the localisation of shear strains after the onset of shear failure in a critical length defined by the diagonal failure planes. Additionally, an effort is made to improve the state of the art in post‐peak shear response modelling, by compiling the largest database of experimental results for shear and flexure‐shear critical R/C columns cycled well beyond the onset of shear failure and/or up to the onset of axial failure, and developing empirical relationships for the key parameters defining the local backbone post‐peak shear response of such elements. The implementation of the derived local hysteretic shear model in a computationally efficient beam‐column finite element model with distributed shear flexibility, which accounts for all deformation types, will be presented in a companion paper.     

Analytical modelling of near‐source pulse‐like seismic demand for multi‐linear backbone oscillators

Nonlinear static procedures, which relate the seismic demand of a structure to that of an equivalent single‐degree‐of‐freedom oscillator, are well‐established tools in the performance‐based earthquake engineering paradigm. Initially, such procedures made recourse to inelastic spectra derived for simple elastic–plastic bilinear oscillators, but the request for demand estimates that delve deeper into the inelastic range, motivated investigating the seismic demand of oscillators with more complex backbone curves. Meanwhile, near‐source (NS) pulse‐like ground motions have been receiving increased attention, because they can induce a distinctive type of inelastic demand. Pulse‐like NS ground motions are usually the result of rupture directivity, where seismic waves generated at different points along the rupture front arrive at a site at the same time, leading to a double‐sided velocity pulse, which delivers most of the seismic energy. Recent research has led to a methodology for incorporating this NS effect in the implementation of nonlinear static procedures. Both of the previously mentioned lines of research motivate the present study on the ductility demands imposed by pulse‐like NS ground motions on oscillators that feature pinching hysteretic behaviour with trilinear backbone curves. Incremental dynamic analysis is used considering 130 pulse‐like‐identified ground motions. Median, 16% and 84% fractile incremental dynamic analysis curves are calculated and fitted by an analytical model. Least‐squares estimates are obtained for the model parameters, which importantly include pulse period T_p. The resulting equations effectively constitute an R ? μ ? T ? T_p relation for pulse‐like NS motions. Potential applications of this result towards estimation of NS seismic demand are also briefly discussed. Copyright © 2016 John Wiley & Sons, Ltd.     

Sand fabric evolution effects on drain design for liquefaction mitigation

This paper revisits the seminal work of Seed and Booker (1977) [21] on the design of infinitely permeable drains for liquefaction mitigation. It is shown that their basic mathematical assumption for the rate of earthquake-induced excess pore pressure generation overlooks sand fabric evolution effects during cyclic loading and eventually leads to underestimation of the drain effectiveness. This is because such effects cause peak excess pore pressures to be attained at the early stages of partially drained shaking, followed by a gradual attenuation even if shaking continues undiminished, a response feature not predicted by the original formulation. In addition, special emphasis is given to the analytical relation describing the excess pore pressure build-up until liquefaction in undrained tests. This relation was considered unique in the original work, for reasons of simplicity, thus neglecting sand fabric evolution effects that may differentiate it for various sands, densities and loading conditions. Hence, a revised analytical formulation is proposed, which takes into account both above effects of sand fabric evolution. The paper provides a quantitative assessment of their influence on drain effectiveness and establishes a new set of charts for drain design. Experimental measurements from shaking table tests, as well as robust numerical simulations are shown, which underline the necessity for the revised solution and design charts.     

Turbidite stacking patterns in salt‐controlled minibasins: Insights from integrated analogue models and numerical fluid flow simulations

The sea floor of intraslope minibasins on passive continental margins plays a significant role in controlling turbidity current pathways and the resulting sediment distribution. To address this, laboratory analogue modelling of intraslope minibasin formation is combined with numerical flow simulations of multi‐event turbidity currents. This approach permits an improved understanding of evolving flow–bathymetry–deposit interactions and the resulting internal stacking patterns of the infills of such minibasins. The bathymetry includes a shelf to slope channel followed by an upper minibasin, which are separated by a confining ridge from two lower minibasins that compares well with analogous bathymetries reported from natural settings. From a wider range of numerical flow experiments, a series of 100 consecutive flows is reported in detail. The turbidity currents are released into the channel and upon reaching the upper minibasin follow a series of stages from short initial ponding, ‘filling and spilling’ and an extended transition to long retrogradational ponding. Upon reaching the upper minibasin floor, the currents undergo a hydraulic jump and therefore much sediment is deposited in the central part of the minibasin and the counterslope. This modifies the bathymetry such that in the fill and spill stage, flow stripping and grain‐size partitioning cause some finer sediment to be transported across the confining ridge into the lower minibasins. Throughout the basin infill process, the sequences retrograde upstream, accompanied by lateral switching into locally formed depressions in the upper minibasin. After the fill and spill stage, significant deposition occurs in the channel where retrograding cyclic steps with wavelengths of 1 to 2 km develop as a function of pulsating flow criticality. These results are at variance with conventional schemes that emphasize sequential downstream minibasin filling through ponding dominated by vertical aggradation. Comparison of these results with published field and experimental examples provides support for the main conclusions.