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.

Impedance of foundations using the Boundary Integral Equation Method

Dynamic soil-structure interaction has been for a long time one of the most fascinating areas for the engineering profession. The building of large alternating machines and their effects on surrounding structures as well as on their own functional behavior, provided the initial impetus; a large amount of experimental research was done,and the results of the Russian and German groups were especially worthwhile. Analytical results by Reissner and Sehkter were reexamined by Quinlan, Sung, et. al., and finally Veletsos presented the first set of reliable results. Since then, the modeling of the homogeneous, elastic halfspace as a equivalent set of springs and dashpots has become an everyday tool in soil engineering practice, especially after the appearance of the fast Fourier transportation algorithm, which makes possible the treatment of the frequency-dependent characteristics of the equivalent elements in a unified fashion with the general method of analysis of the structure. Extensions to the viscoelastic case, as well as to embedded foundations and complicated geometries, have been presented by various authors. In general, they used the finite element method with the well known problems of geometric truncations and the subsequent use of absorbing boundaries. The properties of boundary integral equation methods are, in our opinion, specially well suited to this problem, and several of the previous results have confirmed our opinion. In what follows we present the general features related to steady-state elastodynamics and a series of results showing the splendid results that the BIEM provided. Especially interesting are the outputs obtained through the use of the so-called singular elements, whose description is incorporated at the end of the paper. The reduction in time spent by the computer and the small number of elements needed to simulate realistically the global properties of the halfspace make this procedure one of the most interesting applications of the BIEM.

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.

Dynamic earth pressures against a retaining wall caused by Rayleigh waves

An approximate procedure for studying harmonic soil-structure interaction problems is presented. The presence of Rayleigh waves is considered and the resulting governing equations of the dynamic soil-structure system are solved in the time domain. With this method the transient and steady states of a vibratory motion and also the nonlinear behaviour of the soil can be studied. As an example, the dynamic earth pressure against a rigid retaining wall is investigated. The loads are assumed to be harmonic Rayleigh waves with both static and dynamic surface surcharges. The dependence of the results on the excitation frequency is shown.

Simplified longitudinal seismic analysis of buried tunnels

In this paper an analytical static approach to analyse buried tunnels under seismic surface waves (Rayleigh and Love waves), propagating parallel to the tunnels axis, is provided. In the proposed method, the tunnel is considered as a beam on elastic foundation by using a Winkler model to represent the subgrade reaction and the soil-structure interaction. The seismic load is imposed by giving at the base of the soil springs a determined configuration corresponding to the free-field motion. From the solution of the differential governing equations of the problem, results are obtained in form of relative displacements between points of tunnel, and therefore the seismic bending moments and shearing forces, acting on the tunnel cross section, can be computed.

Revisiting the effect of foundation embedment on elastic settlement: A new approach

The effect of foundation embedment on settlement calculation is a widely researched topic in which there is no scientific consensus regarding the magnitude of settlement reduction. In this paper, a non-linear three dimensional Finite Element analysis has been performed with the aim of evaluating the aforementioned effect. For this purpose, 1800 models were run considering different variables, such as the depth and dimensions of the foundation and the Young’s modulus and Poisson’s ratio of the soil. The settlements from models with foundations at surface level and at depth were then compared and the relationship between them established. The statistical analysis of this data allowed two new expressions, with a mean maximum error of 1.80%, for the embedment influence factor of a foundation to be proposed and these to be compared with commonly used corrections. The proposed equations were validated by comparing the settlements calculated with the proposed influence factors and the true settlements measured in several real foundations. From the comprehensive study of all modelled cases, an improved approach, when compared to those proposed by other authors, for the calculation of the true elastic settlements of an embedded foundation is proposed.

A simple method to predict elastic settlements in foundations resting on two soils of differing deformability

This paper shows the analysis results obtained from more than 200 finite element method (FEM) models used to calculate the settlement of a foundation resting on two soils of differing deformability. The analysis considers such different parameters as the foundation geometry, the percentage of each soil in contact with the foundation base and the ratio of the soils’ elastic moduli. From the described analysis, it is concluded that the maximum settlement of the foundation, calculated by assuming that the foundation is completely resting on the most deformable soil, can be correlated with the settlement calculated by FEM models through a correction coefficient named “settlement reduction factor” (α). As a consequence, a novel expression is proposed for calculating the real settlement of a foundation resting on two soils of different deformability with maximum errors lower than 1.57%, as demonstrated by the statistical analysis carried out. A guide for the application of the proposed simple method is also explained in the paper. Finally, the proposed methodology has been validated using settlement data from an instrumented foundation, indicating that this is a simple, reliable and quick method which allows the computation of the maximum elastic settlement of a raft foundation, evaluates its suitability and optimises its selection process.

Cálculo del asiento elástico tridimensional de cimentaciones de rigidez variable en terrenos con una capa rígida inclinada

Las fórmulas basadas en la teoría de la elasticidad son ampliamente utilizadas para el cálculo de asientos de cimentaciones, ya que la totalidad de la normativa geotécnica recomienda su empleo. No obstante, estos métodos no cubren todas las situaciones geotécnicamente posibles ya que frecuentemente las condiciones geológicas son complejas. En este trabajo se analiza la influencia de la presencia de una capa rígida inclinada en los asientos elásticos de una cimentación superficial. Para ello se han resuelto 273 modelos tridimensionales no lineales de elementos finitos, variando los parámetros clave del problema: la inclinación y la profundidad de la capa rígida y la rigidez de la cimentación. Finalmente, se ha realizado un análisis estadístico de los resultados de los modelos y se ha propuesto una fórmula que puede ser utilizada en el cálculo de asientos por métodos elásticos, para tener en consideración la presencia de una capa rígida inclinada en profundidad.

Updating of an instrumented building model considering amplitude dependence of dynamic resonant properties extracted from seismic response records

The authors would like to express their gratitude to their supporters. Drs Jim Cousins, S.R. Uma and Ken Gledhill facilitated this research by providing access to GeoNet seismic data and structural building information. Piotr Omenzetter’s work within the Lloyd’s Register Foundation Centre for Safety and Reliability Engineering at the University of Aberdeen is supported by Lloyd’s Register Foundation. The Foundation helps to protect life and property by supporting engineering-related education, public engagement and the application of research.

Updating of an instrumented building model considering amplitude dependence of dynamic resonant properties extracted from seismic response records

The authors would like to express their gratitude to their supporters. Drs Jim Cousins, S.R. Uma and Ken Gledhill facilitated this research by providing access to GeoNet seismic data and structural building information. Piotr Omenzetter’s work within the Lloyd’s Register Foundation Centre for Safety and Reliability Engineering at the University of Aberdeen is supported by Lloyd’s Register Foundation. The Foundation helps to protect life and property by supporting engineering-related education, public engagement and the application of research.

Seismic Response Analysis of a Full-Scale Base-Isolated Structure via Measurements and Modeling

The full-scale base-isolated structure studied in this dissertation is the only base-isolated building in South Island of New Zealand. It sustained hundreds of earthquake ground motions from September 2010 and well into 2012. Several large earthquake responses were recorded in December 2011 by NEES@UCLA and by GeoNet recording station nearby Christchurch Women's Hospital. The primary focus of this dissertation is to advance the state-of-the art of the methods to evaluate performance of seismic-isolated structures and the effects of soil-structure interaction by developing new data processing methodologies to overcome current limitations and by implementing advanced numerical modeling in OpenSees for direct analysis of soil-structure interaction.

This dissertation presents a novel method for recovering force-displacement relations within the isolators of building structures with unknown nonlinearities from sparse seismic-response measurements of floor accelerations. The method requires only direct matrix calculations (factorizations and multiplications); no iterative trial-and-error methods are required. The method requires a mass matrix, or at least an estimate of the floor masses. A stiffness matrix may be used, but is not necessary. Essentially, the method operates on a matrix of incomplete measurements of floor accelerations. In the special case of complete floor measurements of systems with linear dynamics, real modes, and equal floor masses, the principal components of this matrix are the modal responses. In the more general case of partial measurements and nonlinear dynamics, the method extracts a number of linearly-dependent components from Hankel matrices of measured horizontal response accelerations, assembles these components row-wise and extracts principal components from the singular value decomposition of this large matrix of linearly-dependent components. These principal components are then interpolated between floors in a way that minimizes the curvature energy of the interpolation. This interpolation step can make use of a reduced-order stiffness matrix, a backward difference matrix or a central difference matrix. The measured and interpolated floor acceleration components at all floors are then assembled and multiplied by a mass matrix. The recovered in-service force-displacement relations are then incorporated into the OpenSees soil structure interaction model.

Numerical simulations of soil-structure interaction involving non-uniform soil behavior are conducted following the development of the complete soil-structure interaction model of Christchurch Women's Hospital in OpenSees. In these 2D OpenSees models, the superstructure is modeled as two-dimensional frames in short span and long span respectively. The lead rubber bearings are modeled as elastomeric bearing (Bouc Wen) elements. The soil underlying the concrete raft foundation is modeled with linear elastic plane strain quadrilateral element. The non-uniformity of the soil profile is incorporated by extraction and interpolation of shear wave velocity profile from the Canterbury Geotechnical Database. The validity of the complete two-dimensional soil-structure interaction OpenSees model for the hospital is checked by comparing the results of peak floor responses and force-displacement relations within the isolation system achieved from OpenSees simulations to the recorded measurements. General explanations and implications, supported by displacement drifts, floor acceleration and displacement responses, force-displacement relations are described to address the effects of soil-structure interaction.

Bridging the Scale Gap from Leaf to Canopy in Biosphere-Atmosphere Gas and Particle Exchanges

Terrestrial ecosystems, occupying more than 25% of the Earth's surface, can serve as

`biological valves' in regulating the anthropogenic emissions of atmospheric aerosol

particles and greenhouse gases (GHGs) as responses to their surrounding environments.

While the signicance of quantifying the exchange rates of GHGs and atmospheric

aerosol particles between the terrestrial biosphere and the atmosphere is

hardly questioned in many scientic elds, the progress in improving model predictability,

data interpretation or the combination of the two remains impeded by

the lack of precise framework elucidating their dynamic transport processes over a

wide range of spatiotemporal scales. The diculty in developing prognostic modeling

tools to quantify the source or sink strength of these atmospheric substances

can be further magnied by the fact that the climate system is also sensitive to the

feedback from terrestrial ecosystems forming the so-called `feedback cycle'. Hence,

the emergent need is to reduce uncertainties when assessing this complex and dynamic

feedback cycle that is necessary to support the decisions of mitigation and

adaptation policies associated with human activities (e.g., anthropogenic emission

controls and land use managements) under current and future climate regimes.

With the goal to improve the predictions for the biosphere-atmosphere exchange

of biologically active gases and atmospheric aerosol particles, the main focus of this

dissertation is on revising and up-scaling the biotic and abiotic transport processes

from leaf to canopy scales. The validity of previous modeling studies in determining

iv

the exchange rate of gases and particles is evaluated with detailed descriptions of their

limitations. Mechanistic-based modeling approaches along with empirical studies

across dierent scales are employed to rene the mathematical descriptions of surface

conductance responsible for gas and particle exchanges as commonly adopted by all

operational models. Specically, how variation in horizontal leaf area density within

the vegetated medium, leaf size and leaf microroughness impact the aerodynamic attributes

and thereby the ultrane particle collection eciency at the leaf/branch scale

is explored using wind tunnel experiments with interpretations by a porous media

model and a scaling analysis. A multi-layered and size-resolved second-order closure

model combined with particle

uxes and concentration measurements within and

above a forest is used to explore the particle transport processes within the canopy

sub-layer and the partitioning of particle deposition onto canopy medium and forest

oor. For gases, a modeling framework accounting for the leaf-level boundary layer

eects on the stomatal pathway for gas exchange is proposed and combined with sap

ux measurements in a wind tunnel to assess how leaf-level transpiration varies with

increasing wind speed. How exogenous environmental conditions and endogenous

soil-root-stem-leaf hydraulic and eco-physiological properties impact the above- and

below-ground water dynamics in the soil-plant system and shape plant responses

to droughts is assessed by a porous media model that accommodates the transient

water

ow within the plant vascular system and is coupled with the aforementioned

leaf-level gas exchange model and soil-root interaction model. It should be noted

that tackling all aspects of potential issues causing uncertainties in forecasting the

feedback cycle between terrestrial ecosystem and the climate is unrealistic in a single

dissertation but further research questions and opportunities based on the foundation

derived from this dissertation are also brie

y discussed.

An experimental investigation into pile group-layered soil interaction

An interaction analysis of an axially loaded single pile and pile group with and without a pile cap in a layered soil medium has been investigated using the two-dimensional photoelastic method. A study of the pile or pile group behaviour has been made, varying the pile cap thickness as well as the embedded length of the pile in the hard stratum. The shear stress distribution along the pile-soil interface, non-dimensionalized settlement values of the single pile and the interaction factor for the pile group have been presented. Wherever possible, the results of the present analysis have been compared with available numerical solutions.

Computer simulation of underground blast response of pile in saturated soil

This paper treats the blast response of a pile foundation in saturated sand using explicit nonlinear finite element analysis, considering complex material behavior of soil and soil–pile interaction. Blast wave propagation in the soil is studied and the horizontal deformation of pile and effective stresses in the pile are presented. Results indicate that the upper part of the pile to be vulnerable and the pile response decays with distance from the explosive. The findings of this research provide valuable information on the effects of underground explosions on pile foundation and will guide future development, validation and application of computer models.

Effect of pile spacing on pile demand for small pile groups in laterally spreading soil

Liquefaction-induced lateral spreading has been responsible for widespread damage to pile foundations in many large earthquakes. The specification of inertial and kinematic pile and pile cap demands is a particularly challenging aspect of the analysis of pile foundations in laterally spreading soils. This paper presents and examines the results from a pair of dynamic centrifuge tests focusing on pile and pile cap demands for small pile groups with different pile spacings. Inertial and kinematic pile cap forces and lateral pile group interaction are examined with regard to the overturning mechanism that dominated the pile group response. © 2014 Taylor & Francis Group.