Dynamic soil pressures on rigid cylindrical vaults

A critical evaluation is made of the dynamic pressures and the associated forces induced by ground shaking on an upright, circular, rigid vault that is embedded in a uniform viscoelastic stratum of constant thickness and infinite extent in the horizontal plane. Both the vault and the stratum are presumed to be supported on a non-deformable base undergoing a space-invariant, uniform horizontal motion. The effects of both harmonic and earthquake-induced excitations are examined. Simple approximate expressions for the responses of the system are formulated, and comprehensive numerical data are presented which elucidate the underlying response mechanisms and the effects and relative importance of the various parameters involved. The parameters investigated include the height to radius ratio for the vault, the conditions at the vault-medium interface, and the material properties of the stratum. In addition to valuable insights into the response of the particular system being examined, the results presented provide a conceptual framework for the analysis and interpretation of solutions for more involved systems as well.     

Inversion of 2D spectral induced polarization imaging data

Laboratory measurements of various materials suggest that more information can be obtained by measuring the in‐phase and out‐of‐phase potentials at a number of frequencies. One common model used to describe the variation of the electrical properties with frequency is the Cole‐Cole model. Apart from the DC resistivity (ρ) and chargeability (m) parameters used in conventional induced‐polarization (IP) surveys, the Cole‐Cole model has two additional parameters, i.e. the time (τ) and relaxation (c) constants. Much research has been conducted on the use of the additional Cole‐Cole parameters to distinguish between different IP sources. Here, we propose a modified inversion method to recover the Cole‐Cole parameters from a 2D spectral IP (SIP) survey. In this method, an approximate inversion method is initially used to construct a non‐homogeneous starting model for the resistivity and chargeability values. The 2D model consists of a number of rectangular cells with constant resistivity (ρ), chargeability (m), time (τ) and relaxation (c) constant values in each cell. A regularized least‐squares optimization method is then used to recover the time and relaxation constant parameters as well as to refine the chargeability values in the 2D model. We present results from tests carried out with the proposed method for a synthetic data set as well as from a laboratory tank experiment.     

A new method to discriminate between a valid IP response and EM coupling effects

The problem of discrimination between a valid induced polarization (IP) response and electromagnetic (EM) coupling effects is considered and an effective solution is provided. First, a finite dimensional approximation to the Cole‐Cole model is investigated. Using the least‐squares approach, the parameters of the approximate model are obtained. Next, based on the analysis of overvoltage, a finite dimensional structure of the IP model is produced. Using this overvoltage‐based structure, a specific finite dimensional approximation of the Cole‐Cole model is proposed. Summarizing the analysis of the finite dimensional IP model, it is concluded that the proposed IP model, which fits the field data much better than the traditional Cole‐Cole model, is essentially an RC‐circuit. From a circuit‐analysis point of view, it is well known that an electromagnetic effect can be described by an RL‐circuit. The simulation results on experimental data support this conception. According to this observation, a new method to discriminate between a valid IP response and EM coupling effects is proposed as follows: (i) use a special finite dimensional model for IP–EM systems; (ii) obtain the parameters for the model using a least‐squares approach; (iii) separate RC‐type terms and RL‐type terms – the first models the IP behaviour, the latter represents the EM part. Simulation on experimental data shows that the method is very simple and effective.     

Modelling tidal signals enhanced by a submarine spring in a coastal confined aquifer extending under the sea

Submarine springs discharge offshore groundwater from confined aquifers extending under the sea. The effects of these springs on the propagation of tidal oscillations in coastal confined aquifers are not known. This paper presents an approximate analytical solution of tidal head fluctuations in a confined aquifer with one submarine spring. The aquifer is assumed to extend in all directions infinitely. The spring is represented by a permeable round column on the seabed, which penetrates completely the impermeable layer overlying the confined aquifer. The error of the approximate solution is negligible if the distance from the spring to the coastline is much greater than the radius of the permeable column representing the spring. Through a hypothetical example, we demonstrate that it is possible to identify the spring’s location using tidal signals observed from inland wells. Tidal groundwater head fluctuations from three inland observation wells at least are needed to determine the 5 model parameters, including the location (2 parameters), the radius of the permeable column representing the spring, the diffusivity of the aquifer, and the tidal loading efficiency of the system.     

On tuned mass dampers for reducing the seismic response of structures

This paper presents an energy‐based theoretical model for a two degree‐of‐freedom mechanical system. After a general formulation in Appendix A, the model is specialized to study tuned mass dampers as a means to substantially increase modal damping in order to induce a consequential decrease of the seismic response of the structures thus provided. Although approximate since it neglects coupling due to damping, it is shown that the model yields a first‐order approximation to the exact frequencies, providing values of optimum damping that closely match exact results proposed by others. In view of this, it is proposed that the model be applied through an iterative numerical procedure that identifies the pertinent optimum parameters. It is also shown that for certain particular benchmark cases the model provides closed‐form equations for the parameters defining the dynamic states related to these special conditions. Despite its approximate nature the model presented in this paper is rational, and due to its explicit consideration of energy balance and overall simplicity, it provides a convenient platform for the study of tuned mass dampers, as well as for other methods of structural passive control. Copyright © 2005 John Wiley & Sons, Ltd.     

Vibration of viscoelastic foundations

Simple approximate solutions are presented for the steady-state response of a massless, rigid disk, which is supported at the surface of a linear viscoelastic halfspace and is excited by a harmonically varying horizontal or vertical force, or by a harmonically varying moment about an axis in the plane of the disk. The halfspace is idealized either as a standard Voigt solid or as a constant hysteretic solid. The solution for the massless disk is then applied to the analysis of foundations with mass, and for each mode of excitation, comprehensive numerical data are presented to illustrate the effects of viscoelastic action on the response of both the massless disk and of foundations with mass It is shown that the principal effects of material damping are to increase the damping capacity of the foundation and to reduce its ‘natural’ frequency. In the mathematically equivalent representation of the elastically supported massless disk as a spring-dashpot combination without mass, it is shown that the effective stiffness of the spring may become negative over significant ranges of the exciting frequency, and that it is preferable to model the disk-halfspace system as an oscillator with mass, taking its spring constant equal to the static stiffness of the original system.     

Introducing new design spectra derived from Italian recorded ground motions 1972 to 2017

The shape of design spectra, traditionally based on regions characterized by constant displacement, constant velocity, and constant acceleration, has been discussed from a conceptual point of view by Calvi (2018). In the same study, a formulation for the definition of the design spectra relying on four parameters was proposed. Predictive models are proposed herein to calculate these four parameters, conditional on magnitude and distance. These models were developed using a large number of recorded ground motions in Italy, and they allow defining combined spectral acceleration versus spectral displacement plots. Such design spectra are shown to reasonably interpolate the experimental data, resulting in acceleration and displacement demand that approximate the response spectra resulting from +1σ results obtained from recorded ground motions. While it is recognized that numerous additional parameter should be considered (eg, focal depth and fault distance, site amplification), it is also concluded that this approach to define the seismic demand is promising toward a rationalization of seismic design. A thorough application of the approach developed and preliminary tested in this work may result in a re‐visitation of seismic design approaches and, ultimately, in a more efficient use of the available resources.     

Analytical solutions to nonlinear conservative oscillator with fifth-order nonlinearity

This paper describes analytical and numerical methods to analyze the steady state periodic response of an oscillator with symmetric elastic and inertia nonlinearity. A new implementation of the homotopy perturbation method (HPM) and an ancient Chinese method called the max-min approach are presented to obtain an approximate solution. The major concern is to assess the accuracy of these approximate methods in predicting the system response within a certain range of system parameters by examining their ability to establish an actual (numerical) solution. Therefore, the analytical results are compared with the numerical results to illustrate the effectiveness and convenience of the proposed methods.     

Evaluation of interaction effects on the system period and the system damping due to foundation embedment and layer depth

On the basis of some simplifying assumptions, a parametric analysis is made of the interaction effects on the effective period and damping of structures with embedded foundation in a soil layer. A simplified three-dimensional interaction model is used, in which the depth of a cylindrical foundation, the degree of contact between the ground and the footing walls and the depth of a homogeneous stratum over rigid rock are considered variable. The soil is replaced with impedance functions that are taken from a data base obtained with an appropriate numerical technique, so that suitable springs and -pots dependent on the excitation frequency are used. The system period and system damping are determined from the steady-state response of an equivalent single oscillator with flexible base subjected to a harmonic motion with constant amplitude, by equating its resonant response with that of a replacement oscillator with rigid base excited with the same motion. The influence of the foundation embedment and soil layer is investigated for several depths of both the footing and the stratum.It is confirmed that the system period decreases and the system damping increases with the foundation embedment only for sidewalls extending along the entire foundation depth. For embedded footings without sidewall or with sidewall in null contact with the surrounding soil, the effective system parameters behave opposite to those corresponding to the interface condition of total contact. Also, the system damping increases significantly with the layer depth, while the system period is practically insensitive to variations of this characteristic parameter. Finally, introducing additional permissible simplifications, an improved approximate solution for the effective period and damping of coupled systems is presented, which differs from previous analogous approximations in that damping factors of second order are not neglected and the foundation depth is explicitly considered.     

Higher-order approximation techniques for estimating stochastic parameter of a sediment transport model

Higher-order approximation techniques for estimating stochastic parameter of the non-homogeneous Poisson (NHP) model are presented. The NHP model is characterized by a two-parameter cumulative probability distribution function (CDF) of sediment displacement. Those two parameters are the temporal and spatial intensity functions, physically representing the inverse of the average rest period and step length of sediment particles, respectively. Difficulty of estimating the parameters has, however, restricted the applications of the NHP model. The approximation techniques are proposed to address such problem. The basic idea of the method is to approximate a model involving stochastic parameters by Taylor series expansion. The expansion preserves certain higher-order terms of interest. Using the experimental (laboratory or field) data, one can determine the model parameters through a system of equations that are simplified by the approximation technique. The parameters so determined are used to predict the cumulative distribution of sediment displacement. The second-order approximation leads to a significant reduction of the CDF error (of the order of 47%) compared to the first-order approximation. Error analysis is performed to evaluate the accuracy of the first- and second-order approximations with respect to the experimental data. The higher-order approximations provide better estimations of the sediment transport and deposition that are critical factors for such environment as spawning gravel-bed.     

Higher-order approximation techniques for estimating stochastic parameter of a sediment transport model

Higher-order approximation techniques for estimating stochastic parameter of the non-homogeneous Poisson (NHP) model are presented. The NHP model is characterized by a two-parameter cumulative probability distribution function (CDF) of sediment displacement. Those two parameters are the temporal and spatial intensity functions, physically representing the inverse of the average rest period and step length of sediment particles, respectively. Difficulty of estimating the parameters has, however, restricted the applications of the NHP model. The approximation techniques are proposed to address such problem. The basic idea of the method is to approximate a model involving stochastic parameters by Taylor series expansion. The expansion preserves certain higher-order terms of interest. Using the experimental (laboratory or field) data, one can determine the model parameters through a system of equations that are simplified by the approximation technique. The parameters so determined are used to predict the cumulative distribution of sediment displacement. The second-order approximation leads to a significant reduction of the CDF error (of the order of 47%) compared to the first-order approximation. Error analysis is performed to evaluate the accuracy of the first- and second-order approximations with respect to the experimental data. The higher-order approximations provide better estimations of the sediment transport and deposition that are critical factors for such environment as spawning gravel-bed.     

Effects of wave passage on the relevant dynamic properties of structures with flexible foundation

An evaluation of the wave passage effects on the relevant dynamic properties of structures with flexible foundation is presented. A simple soil–structure system similar to that used in practice to take into account the inertial interaction effects by the soil flexibility is studied. The kinematic interaction effects due to non‐vertically incident P, SV and Rayleigh waves are accounted for in this model. The effective period and damping of the system are obtained by establishing an equivalence between the interacting system excited by the foundation input motion and a replacement oscillator excited by the free‐field ground motion. In this way, the maximum structural response could be estimated from standard free‐field response spectra using the period and damping of the building modified by both the soil flexibility and the travelling wave effects. Also, an approximate solution for the travelling wave problem is examined over wide ranges of the main parameters involved. Numerical results are computed for a number of soil–structure systems to identify under which conditions the effects of wave passage are important. It comes out that these effects are generally negligible for the system period, but they may significantly change the system damping since the energy dissipation within the soil depends on both the wave radiation and the diffraction and scattering of the incident waves by the foundation. Copyright © 2001 John Wiley & Sons, Ltd.     

Theoretical and experimental study of time domain-induced polarization in water-saturated sands

A theoretical model of spectral-induced polarization (IP) of sand is presented. In the proposed model, contacts of sand grains and intergrain solution-filled space are considered as electrical current passages of varying thickness, which differ in values of ion transport number. Ion-selective narrow passages are considered as active zones, large passages as passive. The proposed model describes spectral IP characteristics for the medium where the length of passive zones is much greater than the length of active ones. The model is called short narrow pores (SNP) model. The SNP model predicts a growth of IP time constant with increase of length of ion-selective zone. Both the time domain and frequency domain parameters are described. The parameters of Cole–Cole model corresponding to the SNP model were also found.The behaviour of model parameters is compared with experimental data obtained on natural and sieved sands using time domain technique. The natural sand spectra correspond neither to the simple SNP model nor simple Cole–Cole model with single time constant because the lengths of ion-selective zones vary, reflecting the grain-size distribution.The spectra of sieved sand compared with the theoretical SNP spectra reveal close correspondence between experimental data and theoretical parameters. For four sieved sands, both the theoretical and experimental data show that the time constant of the IP is proportional to the square of the average grain size.     

Variance component estimation in errors-in-variables models and a rigorous total least-squares approach

A method for variance component estimation (VCE) in errors-in-variables (EIV) models is proposed, which leads to a novel rigorous total least-squares (TLS) approach. To achieve a realistic estimation of parameters, knowledge about the stochastic model, in addition to the functional model, is required. For an EIV model, the existing TLS techniques either do not consider the stochastic model at all or assume approximate models such as those with only one variance component. In contrast to such TLS techniques, the proposed method considers an unknown structure for the stochastic model in the adjustment of an EIV model. It simultaneously predicts the stochastic model and estimates the unknown parameters of the functional model. Moreover the method shows how an EIV model can support the Gauss-Helmert model in some cases. To make the VCE theory into EIV model more applicable, two simplified algorithms are also proposed. The proposed methods can be applied to linear regression and datum transformation. We apply these methods to these examples. In particular a 3-D non-linear close to identical similarity transformation is performed. Two simulation studies besides an experimental example give insight into the efficiency of the algorithms.     

地基土-桩基础-核电站辅助厂房结构相互作用体系的地震响应分析

本文研究了地基土-桩-核岛辅助厂房结构相互作用体系在地震激励时的“平-扭”耦联动力响应,在分析中分别采用M.Novak的近似理论,M.Novak常数法,M法,常数法计算了桩基础的阻抗函数,并考虑了在同一基础上的安全壳和另一辅助厂房对所研究的辅助厂房动力响应的影响。     

A web‐based methodology for the prediction of approximate IDA curves

A web‐based methodology for the prediction of approximate IDA curves, which consists of two independent processes, is proposed. The result of the first process is a response database of the SDOF model, whereas the second process involves the prediction of approximate IDA curves from the response database by using n‐dimensional linear interpolation. Such an approach enables user‐friendly prediction of the seismic response parameters with high accuracy. In order to demonstrate the capabilities of the proposed methodology, a web application for the prediction of the approximate 16th, 50th and 84th fractile responses of an RC structure was developed. For the presented case study, the response database was computed for a set of 30 ground motion records and the discrete values of six structural parameters. Very good agreement between the computed and the approximated IDA curves of the four‐storey RC building was observed. Copyright © 2012 John Wiley & Sons, Ltd.     

Knowledge,transparency, and refutability in groundwater models,an example from the Death Valley regional groundwater flow system

This work demonstrates how available knowledge can be used to build more transparent and refutable computer models of groundwater systems. The Death Valley regional groundwater flow system, which surrounds a proposed site for a high level nuclear waste repository of the United States of America, and the Nevada National Security Site (NNSS), where nuclear weapons were tested, is used to explore model adequacy, identify parameters important to (and informed by) observations, and identify existing old and potential new observations important to predictions. Model development is pursued using a set of fundamental questions addressed with carefully designed metrics. Critical methods include using a hydrogeologic model, managing model nonlinearity by designing models that are robust while maintaining realism, using error-based weighting to combine disparate types of data, and identifying important and unimportant parameters and observations and optimizing parameter values with computationally frugal schemes. The frugal schemes employed in this study require relatively few (10–1000 s), parallelizable model runs. This is beneficial because models able to approximate the complex site geology defensibly tend to have high computational cost. The issue of model defensibility is particularly important given the contentious political issues involved.     

Experimental verification of torsional response control of asymmetric buildings using MR dampers

This paper proposes a semiactive control system to reduce the coupled lateral and torsional motions in asymmetric buildings subjected to horizontal seismic excitations. Magnetorheological (MR) dampers are applied as semiactive control devices and the control input determination is based on a clipped‐optimal control algorithm which uses absolute acceleration feedback. The performance of this method is studied experimentally using a 2‐story building model with an asymmetric stiffness distribution. An automated system identification methodology is implemented to develop a control‐oriented model which has the natural frequencies observed in the experimental system. The parameters for the MR damper model are identified using experimental data to develop an integrated model of the structure and MR dampers. To demonstrate the performance of this control system on the experimental structure, a shake table is used to reproduce an El Centro 1940 N–S earthquake as well as a random white noise excitation. The responses for the proposed control system are compared to those of passive control cases in which a constant voltage is applied to the MR damper. Copyright © 2003 John Wiley & Sons, Ltd.     

Model tests on horizontal pile-to-pile interaction incorporating local non-linearity and resonance effects

Dynamic experimental studies on horizontal interaction factors for laterally loaded model soil–pile systems subjected to low-to-high amplitude pile head loading are reported, focusing on the effect of dynamics and local non-linearity. Results obtained from dynamic experiments with different amplitudes of harmonic lateral loading on instrumented model soil–pile systems conducted on a shaking table indicate that the resonance of soil–pile system shows a profound impact on the horizontal interaction factors. Furthermore, behavior of the soil–pile system becomes highly non-linear with increasing amplitude of loading, as the extent of local non-linearity around the pile increases. Consequently, group effects cannot be predicted by using established approximate equations based on the assumption of linear viscoelastic behavior. To this end, new semi-empirical equations are proposed for obtaining the horizontal interaction factors, as an extension of available approximate equations to incorporate the local non-linear behavior of soil–pile system for both small and large pile spacings.