Inelastic displacement demand of bridge columns considering shear–flexure interaction

Reinforced concrete bridge columns exhibit complex hysteretic behavior owing to combined action of shear, bending moment, and axial force under multi‐directional seismic shakings. The inelastic displacement of columns can be increased by shear–flexure interaction (SFI). This paper develops a simple yet reliable demand model for estimating the inelastic displacement and ductility based on the nonlinear time history analyses of 24 full‐size columns subject to a suite of near‐fault ground motions. A coupled hysteretic model is used to simulate the shear‐flexure interactive (SFI) behavior of columns and the accumulated material damage during loading reversals, including pinching, strength deterioration, and stiffness softening. Guided by rigorous dimensional analysis, the inelastic displacement responses of bridge columns are presented in dimensionless form showing remarkable order. A dimensionless nonlinearity index is derived taking into account of the column strength, ground motion amplitude, and softening or hardening post‐yield behavior. Strong correlation is revealed between the normalized inelastic displacement and the dimensionless structure‐to‐pulse frequency, the dimensionless nonlinearity index as well as the aspect ratio. Two regressive equations for displacement and ductility demands are proposed and validated against the simulation results. The SFI effects are discussed and included explicitly through the aspect ratio in the proposed model. This study offers a new way to realistically predict the inelastic displacement of columns directly from structural and ground motion characteristics. Copyright © 2010 John Wiley & Sons, Ltd.     

Optimum lateral load pattern for seismic design of elastic shear‐buildings incorporating soil–structure interaction effects

Recently, several new optimum loading patterns have been proposed by researchers for fixed‐base systems while their adequacy for soil–structure systems has not been evaluated yet. Through intensive dynamic analyses of multistory shear‐building models with soil–structure interaction subjected to a group of 21 artificial earthquakes adjusted to soft soil design spectrum, the adequacy of these optimum patterns is investigated. It is concluded that using these patterns the structures generally achieve near optimum performance in some range of periods. However, their efficiency reduces as soil flexibility increases especially when soil–structure interaction effects are significant. In the present paper, using the uniform distribution of damage over the height of structures, as the criterion, an optimization algorithm for seismic design of elastic soil–structure systems is developed. The effects of fundamental period, number of stories, earthquake excitation, soil flexibility, building aspect ratio, damping ratio and damping model on optimum distribution pattern are investigated. On the basis of 30,240 optimum load patterns derived from numerical simulations and nonlinear statistical regression analyses, a new lateral load pattern for elastic soil–structure systems is proposed. It is a function of the fundamental period of the structure, soil flexibility and structural slenderness ratio. It is shown that the seismic performance of such a structure is superior to those designed by code‐compliant or recently proposed patterns by researchers for fixed‐base structures. Using the proposed load pattern in this study, the designed structures experience up to 40% less structural weight as compared with the code‐compliant or optimum patterns developed based on fixed‐base structures. Copyright © 2012 John Wiley & Sons, Ltd.     

R/C frame–infill interaction model and its application to Indonesian buildings

This paper proposes a new analytical model for masonry‐infilled R/C frames to evaluate the seismic performance considering R/C frame–infill interactions. The proposed analytical model replaces masonry infill with a diagonal compression strut, which represents distributed compression transferred between frame and infill interfaces. The equivalent strut width is presented as a function of the frame–infill contact length, which can be evaluated by static equilibriums related to compression balance and lateral displacement compatibility at the frame–infill interfaces. The proposed analytical model was verified through comparisons with experimental results obtained for several brick masonry‐infilled R/C frames representing a typical R/C building with nonstructural masonry infill in Indonesia. As a result, good agreements were observed between the experimental and analytical values of the lateral strength and ductility of the infilled frames. The seismic performances of two earthquake‐damaged R/C buildings with different damage conditions were evaluated considering infill effects by applying the proposed analytical model. Consequently, the nonstructural brick masonry infill significantly affected the seismic resistances of the buildings, which seemed to lead to differing levels of damage for each building. These results indicate that the proposed analytical model can be an effective tool for more precisely screening earthquake‐vulnerable existing R/C buildings in Indonesia. Copyright © 2016 The Authors. Earthquake Engineering & Structural Dynamics Published by John Wiley & Sons Ltd.     

Three-dimensional nonlinear analysis for seismic soil–pile-structure interaction

A three-dimensional method of analysis is presented for the seismic response of structures constructed on pile foundations. An analysis is formulated in the time domain and the effects of material nonlinearity of soil on the seismic response are investigated. A subsystem model consisting of a structure subsystem and a pile-foundation subsystem is used. Seismic response of the system is found using a successive-coupling incremental solution scheme. Both subsystems are assumed to be coupled at each time step. Material nonlinearity is accounted for by incorporating an advanced plasticity-based soil model, HiSS, in the finite element formulation. Both single piles and pile groups are considered and the effects of kinematic and inertial interaction on seismic response are investigated while considering harmonic and transient excitations. It is seen that nonlinearity significantly affects seismic response of pile foundations as well as that of structures. Effects of nonlinearity on response are dependent on the frequency of excitation with nonlinearity causing an increase in response at low frequencies of excitation.     

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.     

Influence of frequency‐dependent soil–structure interaction on the fragility of R/C bridges

Bridge performance under earthquake loading can be significantly influenced by the interaction between the structure and the supporting soil. Even though the frequency dependence of the interaction mentioned in this study has long been documented, the simplifying assumption that the dynamic stiffness is dominated by the mean or predominant excitation frequency is still commonly made, primarily as a result of the associated numerical difficulties when the analysis has to be performed in the time domain. This study makes use of the advanced lumped parameter models recently developed 1 in order to quantify the impact of the assumption on the predicted fragility of bridges mentioned in this study. This is achieved by comparing the predicted vulnerability for the case of a reference, well studied, actual bridge using both conventional, frequency‐independent, Kelvin–Voigt models and the aforementioned lumped parameter formulation. Analysis results demonstrate that the more refined consideration of frequency dependence of soil–structure interaction at the piers and the abutments of a bridge not only leads to different probabilities of failure for given intensity measures but also leads to different hierarchy and distribution of damage within the structure for the same set of earthquake ground motions even if the overall probability of exceeding a given damage state is the same. The paper concludes with the comparative assessment of the effect for different soil conditions, foundation configurations, and ground motion characteristics mentioned in this study along with the relevant analysis and design recommendations. Copyright © 2016 John Wiley & Sons, Ltd.     

Pipe–soil interaction for segmented buried pipelines subjected to dip faults

This paper describes an investigation of pipe–soil interaction equations suggested by currently used pipeline seismic design codes and the applicability of these equations to segmented pipelines. The results of computer‐aided analyses were compared to results obtained in full‐scale experiments on a segmented ductile iron pipeline 93 mm in diameter and 15 m in length. The pipeline was installed 600 mm below the ground surface in a sandy soil compacted to two different subgrade reaction values. The type of fault considered was a reverse fault with an intersection angle of 60° with the pipeline, and the fault movement was a total of 350 mm in three same steps in the fault trace direction. The findings of this study demonstrate the necessity of considering the nature of soil behavior in pipe–soil interaction equations and the effects of connection joints on the integrated response of pipelines to fault‐induced ground deformations. A new combination of equations constituting a direction‐wise selection from among the equations proposed by currently used guidelines is introduced as a new series to describe pipe–soil interaction for segmented pipelines and is verified using the results of full‐scale experiments. Copyright © 2014 John Wiley & Sons, Ltd.     

A coupling model of FE–BE–IE–IBE for non-linear layered soil–structure interactions

A coupling model of Finite Elements (FEs), Boundary Elements (BEs), Infinite Elements (IEs) and Infinite Boundary Elements (IBEs) is presented for analysis of soil–structure interaction (SSI). The radiation effects of the infinite layered soil are taken into account by FE–IE coupling, while the underlying bed rock half-space is discretized into BE–IBE coupling whereby the non-horizontal bed rock surface can be accounted for. Displacement compatabilities are satisfied for all types of aforementioned elements. The equivalent linear approach is employed for approximation of nonlinearity of the near field soil. This model has some advantages over the current SSI program in considering the bed rock half-space and non-vertical wave incidence from the far field. Examples of verification demonstrate the applicability and accuracy of the method when compared with the FLUSH program. Finally, the effects of the relative modulus ratio E_r/E_s of rock and soil and the incident angles of non-vertical waves on the responses of the structure and the soil are examined. Copyright © 1999 John Wiley & Sons, Ltd.     

A discontinuous Galerkin method for seismic soil–structure interaction analysis in the time domain

A technique for modeling transient wave propagation in unbounded media is extended and applied to seismic soil–structure interaction analysis in the time domain. The technique, based on the discontinuous Galerkin method, requires lower computational cost and less storage than the boundary element method, and the time‐stepping scheme resulting from Newmark's method in conjunction with the technique is unconditionally stable, allowing for efficient and robust time‐domain computations. To extend the technique to cases characterized by seismic excitation, the free‐field motion is used to compute effective forces, which are introduced on the boundary of the computational domain containing the structure and the soil in the vicinity of the structure. A numerical example on a dam–foundation system subjected to seismic excitation demonstrates the performance of the method. Copyright © 2003 John Wiley & Sons, Ltd.     

Prediction of free field vibrations due to pile driving using a dynamic soil–structure interaction formulation

This paper presents a numerical model for the prediction of free field vibrations due to vibratory and impact pile driving. As the focus is on the response in the far field where deformations are relatively small, a linear elastic constitutive behaviour is assumed for the soil. The free field vibrations are calculated by means of a coupled FE–BE model based on a subdomain formulation. First, the case of vibratory pile driving is considered, where the contributions of different types of waves are investigated for several penetration depths. In the near field, the soil response is dominated by a vertically polarized shear wave, whereas in the far field, body waves are importantly attenuated and Rayleigh waves dominate the ground vibration. Second, the case of impact pile driving is considered. A linear wave equation model is used to estimate the impact force during the driving process. Apart from the response of a homogeneous halfspace, it is also investigated how the soil stratification influences the ground vibration for the case of a soft layer on a stiffer halfspace. When the penetration depth is smaller than the layer thickness, the layered medium has no significant influence on ground vibrations. However, when the penetration depth is larger than the layer thickness, the influence of the layered medium becomes more significant. The computed ground vibrations are finally compared with field measurements reported in the literature.     

Influence of dynamic soil–structure interaction on the nonlinear response and seismic reliability of multistorey systems

A set of reinforced concrete structures with gravitational loads and mechanical properties (strength and stiffness) representative of systems designed for earthquake resistance in accordance with current criteria and methods is selected to study the influence of dynamic soil–structure interaction on seismic response, ductility demands and reliability levels. The buildings are considered located at soft soil sites in the Valley of Mexico and subjected to ground motion time histories simulated in accordance with characteristic parameters of the maximum probable earthquake likely to occur during the system's expected life. For the near‐resonance condition the effects of soil–structure interaction on the ductility demands depend mainly on radiation damping. According to the geometry of the structures studied this damping is strongly correlated with the aspect ratio, obtained by dividing the building height by its width. In this way, for structures with aspect ratio greater than 1.4 the storey and global ductility demands increase with respect to those obtained with the same structures but on rigid base, while for structures with aspect ratio less than 1.4 the ductility demands decrease with respect to those for the structures on rigid base. For the cases when the fundamental period of the structure has values very different from the dominant ground period, soil–structure interaction leads in all cases to a reduction of the ductility demands, independently of the aspect ratio. The reliability index β is obtained as a function of the base shear ratio and of the seismic intensity acting on the nonlinear systems subjected to the simulated motions. The resulting reliability functions are very similar for systems on rigid or on flexible foundation, provided that in the latter case the base rotation and the lateral displacement are removed from the total response of the system. Copyright © 2006 John Wiley & Sons, Ltd.     

A hybrid time frequency domain formulation for non-linear soil–structure interaction

Methods that combine frequency and time domain techniques offer an attractive alternative for solving Soil–Structure-interaction problems where the structure exhibits non-linear behaviour. In the hybrid-frequency-time-domain procedure a reference linear system is solved in the frequency domain and the difference between the actual restoring forces and those in the linear model are treated as pseudo-forces. In the solution scheme explored in this paper, designated as the hybrid-time-frequency-domain (HTFD) procedure, the equations of motion are solved in the time domain with due consideration for non-linearities and with the unbounded medium represented by frequency-independent springs and dampers. The frequency dependency of the impedance coefficients is introduced by means of pseudo-forces evaluated in the frequency domain at the end of each iteration. A criterion of stability for the HTFD approach is derived analytically and its validity is sustained numerically. As is often the case, the criterion takes the form of a limit of unity on the spectral radius of an appropriately defined matrix. Inspection of the terms in this matrix shows that convergence can be guaranteed by suitable selection of the reference impedance. The CPU times required to obtain converged solutions with the HTFD are found, in a number of numerical simulations, to be up to one order of magnitude less than those required by the alternative hybrid-frequency-time-domain approach. © 1998 John Wiley & Sons, Ltd.     

Comparison of genetic algorithms and shuffled complex evolution approach for calibrating distributed rainfall–runoff model

Three methods, Shuffled Complex Evolution (SCE), Simple Genetic Algorithm (SGA) and Micro‐Genetic Algorithm (µGA), are applied in parameter calibration of a grid‐based distributed rainfall–runoff model (GBDM) and compared by their performances. Ten and four historical storm events in the Yan‐Shui Creek catchment, Taiwan, provide the database for model calibration and verification, respectively. The study reveals that the SCE, SGA and µGA have close calibration results, and none of them are superior with respect to all the performance measures, i.e. the errors of time to peak, peak discharge and the total runoff volume, etc. The performances of the GBDM for the verification events are slightly worse than those in the calibration events, but still quite satisfactory. Among the three methods, the SCE seems to be more robust than the other two approaches because of the smallest influence of different initial random number seeds on calibrated model parameters, and has the best performance of verification with a relatively small number of calibration events. Copyright © 2009 John Wiley & Sons, Ltd.     

A geomorphology‐based integrated stream–aquifer interaction model for semi‐gauged catchments

Many of the existing stream–aquifer interaction models available in the literature are very complex with limited applicability in semi‐gauged and ungauged catchments. In this study, to estimate the influent and effluent subsurface water fluxes under limited geo‐hydrometeorological data availability conditions, a simple stream–aquifer interaction model, namely, the variable parameter McCarthy–Muskingum (VPMM) hillslope‐storage Boussinesq (hsB) model, has been developed. This novel model couples the VPMM streamflow transport with the hsB groundwater flow transport modules in online mode. In this integrated model, the surface water–groundwater flux exchange process is modelled by the Darcian approach with the variable hydraulic heads between the river stage and groundwater table accounting for the rainfall forcing. Considering the exchange fluxes in the hyporheic zone and lateral overland flow contribution, this approach is field tested in a typical 48‐km stretch of the Brahmani River in eastern India to simulate the streamflow and its depth with the minimum Nash–Sutcliffe efficiency of 94% and 88%; the maximum root mean square error of 134 m³/s and 0.35 m; and the minimum index of agreement of 98% and 97%, respectively. This modelling approach could be very well utilized in data‐scarce world‐river basins to estimate the stream–aquifer exchange flux due to rainfall forcings.     

Fluid–structure interaction analysis of 3-D rectangular tanks by a variationally coupled BEM–FEM and comparison with test results

A variationally coupled BEM–FEM is developed which can be used to analyse dynamic response, including free-surface sloshing motion, of 3-D rectangular liquid storage tanks subjected to horizontal ground excitation. The tank structure is modelled by the finite element method and the fluid region by the indirect boundary element method. By minimizing a single Lagrange function defined for the entire system, the governing equation with symmetric coefficient matrices is obtained. To verify the newly developed method, the analysis results are compared with the shaking-table test data of a 3-D rectangular tank model and with the solutions by the direct BEM–FEM. Analytical studies are conducted on the dynamic behaviour of 3-D rectangular tanks using the method developed. In particular, the characteristics of the sloshing response, the effect of the rigidity of adjacent walls on the dynamic response of the tanks and the orthogonal effects are investigated. © 1998 John Wiley & Sons, Ltd.     

Non-singular time domain BEM with applications to 3D inertial soil–structure interaction

A direct time domain boundary element method is presented based on the Stokes fundamental solutions, discretized in both time and space, and an efficient time step-by-step solution that minimizes the accumulation of errors. A non-singular numerical integration procedure, in the Cauchy sense, is proposed for the generation of the associated influence matrices. This methodology is shown to be efficient for the solution of a number of computationally intensive problems in the area of soil–structure interaction. In addition, an algorithm for the direct calculation of the response of massive foundations to externally applied forces and/or obliquely incident seismic waves is introduced. The accuracy and computational efficiency of the proposed methodologies is established through a number of comparison studies.     

A comparison of fibre‐optic distributed temperature sensing to traditional methods of evaluating groundwater inflow to streams

There are several methods for determining the spatial distribution and magnitude of groundwater inputs to streams. We compared the results of conventional methods [dye dilution gauging, acoustic Doppler velocimeter (ADV) differential gauging, and geochemical end‐member mixing] to distributed temperature sensing (DTS) using a fibre‐optic cable installed along 900 m of Ninemile Creek in Syracuse, New York, USA, during low‐flow conditions (discharge of 1·4 m³ s^?1). With the exception of differential gauging, all methods identified a focused, contaminated groundwater inflow and produced similar groundwater discharge estimates for that point, with a mean of 66·8 l s^?1 between all methods although the precision of these estimates varied. ADV discharge measurement accuracy was reduced by non‐ideal conditions and failed to identify, much less quantify, the modest groundwater input, which was only 5% of total stream flow. These results indicate ambient tracers, such as heat and geochemical mixing, can yield spatially and quantitatively refined estimates of relatively modest groundwater inflow even in large rivers. DTS heat tracing, in particular, provided the finest spatial characterization of groundwater inflow, and may be more universally applicable than geochemical methods, for which a distinct and consistent groundwater end member may be more difficult to identify. Copyright © 2011 John Wiley & Sons, Ltd.     

H-adaptivity applied to liquefiable soil in nonlinear analysis of soil–pile interaction

This paper is concerned with application of the h-adaptive finite element method to dynamic analysis of a pile in liquefiable soil considering large deformation. In finite element analysis of pile behavior in liquefiable soil during an earthquake, especially considering large deformation of liquefied ground, error due to discretization in the zone near the pile becomes very large. Our purpose was to refine the approximation of the finite element method. The updated Lagrangian formulation and a cyclic elasto-plastic model based on the kinematic hardening rule were adopted to deal with the nonlinearity of the soil. The mixed finite element and finite difference methods together with the u-p formulation and Biot's two-phase mixture theory were used. To improve the accuracy and increase the efficiency of finite element analysis, an h-adaptive scheme that included a posteriori error estimation and h-version mesh refinement was applied to the analysis. The calculated results of effective stress were smoothed locally by the extrapolation method and smoothed stress was used to calculate the L₂ norm of the effective stress error in the last step of the calculation of each time increment. The mesh was refined by a fission procedure based on the indication of the error estimate As a numerical example, a soil–pile interaction system loaded cyclically was analyzed by our method.