Critical damping of structures with elastically supported visco‐elastic dampers

This paper presents a new formulation for critical damping of structures with elastically supported visco‐elastic dampers.Owing to the great dependence of damper performance on the support stiffness, this model is inevitable for reliable modelling of structures with visco‐elastic dampers. It is shown that the governing equation of free vibration of this model is reduced to a third‐order differential equation and the conventional method for defining the critical damping for second‐order differential equations cannot be applied to the present model. It is demonstrated that the region of overdamped vibration is finite in contrast to that (semi‐infinite) for second‐order differential equations and multiple critical damping coefficients exist. However, it turns out that the smaller one is practically meaningful. Copyright © 2001 John Wiley & Sons, Ltd.     

The effects of reservoir geometry on the seismic response of gravity dams

Conventional seismic analysis of gravity dams assumes that the behaviour of the dam–water–soil system can be represented using a 2‐D model since dam vertical contraction joints between blocks allow them to vibrate independently from each other. The 2‐D model assumes the reservoir to be infinite and of constant width, which is not true for certain types of reservoirs. In this paper, a boundary element method (BEM) model in the frequency domain is used to investigate the influence of the reservoir geometry on the hydrodynamic dam response. Important conceptual conclusions about the dam–reservoir system behaviour are obtained using this model. The results show that the reservoir shape influences the seismic response of the dam, making it necessary to account for 3‐D effects in order to obtain accurate results. In particular, the 3‐D pressure and displacement responses can be substantially larger than those computed with the 2‐D model. Copyright © 2007 John Wiley & Sons, Ltd.     

Seismic response of intake towers including dam–tower interaction

The seismic response of the intake–outlet towers has been widely analyzed in recent years. The usual models consider the hydrodynamic effects produced by the surrounding water and the interior water, characterizing the dynamic response of the tower–water–foundation–soil system. As a result of these works, simplified added mass models have been developed. However, in all previous models, the surrounding water is assumed to be of uniform depth and to have infinite extension. Consequently, the considered added mass is associated with only the pressures created by the displacements of the tower itself. For a real system, the intake tower is usually located in proximity to the dam and the dam pressures may influence the equivalent added mass. The objective of this paper is to investigate how the response of the tower is affected by the presence of the dam. A coupled three‐dimensional boundary element‐finite element model in the frequency domain is employed to analyze the tower–dam–reservoir interaction problem. In all cases, the system response is assumed to be linear, and the effect of the internal fluid and the soil–structure interaction effects are not considered. The results suggest that unexpected resonance amplifications can occur due to changes in the added mass for the tower as a result of the tower–dam–reservoir interaction. Copyright © 2008 John Wiley & Sons, Ltd.     

Representation of radiation damping by fractional time derivatives

When modelling unbounded domains, formulation of a matrix‐valued force–displacement relationship which can take radiation damping into account is of major importance. In this paper, a method to describe the dynamic stiffness by a system of fractional differential equations in the time‐domain is presented. Here, a doubly asymptotic rational approximation of the low‐frequency force–displacement relationship is used, whereas a direct interpretation of the asymptotic part as a fractional derivative is possible. The numerical solution of the corresponding system of fractional differential equations is demonstrated using the infinite beam on elastic foundation as an example. Copyright © 2003 John Wiley & Sons, Ltd.     

An efficient approach for frequency-domain and time-domain hydrodynamic analysis of dam–reservoir systems

An efficient procedure is developed for the hydrodynamic analysis of dam–reservoir systems. The governing equations of hydrodynamic pressure in the frequency as well as time domain are derived in the framework of the scaled boundary finite element method. The water compressibility and absorption of reservoir sediments can be conveniently taken into consideration. By extending the reservoir to infinity with uniform cross-section, only the dam–reservoir interface needs to be discretized to model the fluid domain, and the hydrodynamic pressure in the stream direction is solved analytically. Several numerical examples including a gravity dam with an inclined upstream face and an arch dam with a reservoir of arbitrary cross-section are provided to demonstrate the computational efficiency and accuracy of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydrodynamic effects in a reservoir with semicircular cross-section and absorptive bottom

A three-dimensional dam-reservoir system under seismic load is analysed. The dam is assumed to be rigid. The reservoir is an infinite channel with semi-circular cross-section. The exact analytical solution, based on the assumption of potential fluid motion is presented, as well as numerical results for selected parameters.The most significant parameters are: the direction and frequency content of the seismic input; the radiation damping at the reservoir bottom; and the compressibility of the fluid. The response of the system depends strongly on the direction of the input ground motion. This is shown by the transfer functions as well as by the pressure time histories due to two earthquakes with different frequency content. The energy absorption at the reservoir bottom is important. A simple plane-wave model shows, that even for a rock foundation, the amount of transmitted energy can reach up to 80%. For comparison the case without bottom absorption is also shown. Compressbility has to be included to capture the resonance effects. The exact analytical solution is also used to verify numerical results obtained by a new method that combines a finite element model with a rigorous radiation boundary for the infinite channel in the time domain.     

Direct finite element method for nonlinear earthquake analysis of 3‐dimensional semi‐unbounded dam–water–foundation rock systems

A direct finite element method for nonlinear earthquake analysis of 2‐dimensional dam–water–foundation rock systems has recently been presented. The analysis procedure uses standard viscous‐damper absorbing boundaries to model the semi‐unbounded foundation‐rock and fluid domains and specifies the seismic input as effective earthquake forces at these boundaries. Presented in this paper is a generalization of the direct finite element method with viscous‐damper boundaries to 3‐dimensional dam–water–foundation rock systems. Step‐by‐step procedures for determining the effective earthquake forces starting from a ground motion specified at a control point on the foundation‐rock surface is developed, and several numerical examples are computed and compared with independent benchmark solutions to demonstrate the effectiveness of the analysis procedure for modeling 3‐dimensional systems.     

2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method

To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion.     

Direct finite element method for nonlinear analysis of semi‐unbounded dam–water–foundation rock systems

A direct finite element method is presented for nonlinear earthquake analysis of interacting dam–water–foundation rock systems. The analysis procedure applies viscous damper absorbing boundaries to truncate the semi‐unbounded fluid and foundation‐rock domains and specifies at these boundaries effective earthquake forces determined from the design ground motion defined at a control point on the free surface. The analysis procedure is validated numerically by computing the frequency response functions and transient response of an idealized dam–water–foundation rock system and comparing with results from the substructure method. Because the analysis procedure is applicable to nonlinear systems, it allows for modeling of concrete cracking, as well as sliding and separation at construction joints, lift joints, and at concrete–rock interfaces. Implementation of the procedure is facilitated by commercial finite element software with nonlinear material models that permit modeling of viscous damper boundaries and specification of effective earthquake forces at these boundaries. Copyright © 2017 John Wiley & Sons, Ltd.     

Nonlinear earthquake analysis of high arch dam–water–foundation rock systems

A nonlinear finite element model for earthquake response analysis of arch dam–water–foundation rock systems is proposed in this paper. The model includes dynamic dam–water and dam–foundation rock interactions, the opening of contraction joints, the radiation damping of semi‐unbounded foundation rock, the compressibility of impounded water, and the upstream energy propagating along the semi‐unbounded reservoir. Meanwhile, a new equivalent force scheme is suggested to achieve free‐field input in the model. The effects of the earthquake input mechanism, joint opening, water compressibility, and radiation damping on the earthquake response of the Ertan arch dam (240 m high) in China are investigated using the proposed model. The results show that these factors significantly affect the earthquake response of the Ertan arch dam. Such factors should therefore be considered in the earthquake response analysis and earthquake safety evaluation of high arch dams. Copyright © 2011 John Wiley & Sons, Ltd.     

二维土坝的两相动力分析

提出了一种分析饱和土坝动力反应的方法,考虑了土坝的两相介质特性,在固液耦联动力方程的基础上,选取固相位移,液相位移、孔隙水压作为场变量,采用伽辽金加权残数法进行有限元空间离散化,然后在时域上采用Ｗｉｌｓｏｎ－θ法进行逐步积分。该方法不仅能计算出固相位移和液相位移,而且能直接得到孔隙水压的反应过程。文中以一饱和土坝模型进行算例分析,并与将其作为单相介质时的结果进行了比较。该法可用于分析饱和介质的地震     

A FE–BE method for dynamic analysis of dam–foundation–reservoir systems in the time domain

By coupling FEM and BEM, a numerical method was developed for dynamic response analyses of dam–foundation–reservoir systems in the time domain. During formulation, the weighted residual procedure was applied to the coupling of several equations of motion for solid and fluid in the FE and BE regions, and an algorithm similar to the Newmark beta procedure was finally obtained. The algorithm is advantageous in that it takes into account all the effects of dam–foundation, dam–reservoir and reservoir–foundation interactions, as well as of the absorption of both elastodynamic and hydrodynamic waves at the boundaries of the foundation and the reservoir. To demonstrate the validity of the present method, the impulsive response of a dam–foundation–reservoir system was calculated using the algorithm, and showed a good agreement with the existing results obtained by other researchers.     

Response of unbounded soil in scaled boundary finite‐element method

The scaled boundary finite‐element method is a powerful semi‐analytical computational procedure to calculate the dynamic stiffness of the unbounded soil at the structure–soil interface. This permits the analysis of dynamic soil–structure interaction using the substructure method. The response in the neighbouring soil can also be determined analytically. The method is extended to calculate numerically the response throughout the unbounded soil including the far field. The three‐dimensional vector‐wave equation of elasto‐dynamics is addressed. The radiation condition at infinity is satisfied exactly. By solving an eigenvalue problem, the high‐frequency limit of the dynamic stiffness is constructed to be positive definite. However, a direct determination using impedances is also possible. Solving two first‐order ordinary differential equations numerically permits the radiation condition and the boundary condition of the structure–soil interface to be satisfied sequentially, leading to the displacements in the unbounded soil. A generalization to viscoelastic material using the correspondence principle is straightforward. Alternatively, the displacements can also be calculated analytically in the far field. Good agreement of displacements along the free surface and below a prism foundation embedded in a half‐space with the results of the boundary‐element method is observed. Copyright © 2001 John Wiley & Sons, Ltd.     

Two-dimensional dynamic analysis of concrete gravity and embankment dams including hydrodynamic effects

An analysis procedure in the frequency domain is developed for determining the earthquake response of two-dimensional concrete gravity and embankment dams including hydrodynamic effects; responses of the elastic dams and compressible water are assumed linear. The dam and fluid domain are treated as substructures and modelled with finite elements. The only geometric restriction is that an infinite fluid domain must maintain a constant depth beyond some point in the upstream direction. For such an infinite uniform region, a finite element discretization over the depth is combined with a continuum representation in the upstream direction. The fluid domain model approximately accounts for interaction between the fluid and underlying foundation medium through a damping boundary condition applied along the reservoir bottom, while the dam foundation is assumed rigid. Several examples are presented to demonstrate the accuracy of the fluid domain model and to illustrate dam responses obtained from the analysis procedure.     

Nonstationary seismic responses of structure with nonlinear stiffness subject to modulated Kanai–Tajimi excitation

A simple structure under earthquake excitation is modeled as a single‐degree‐of‐freedom system with nonlinear stiffness subject to modulated Kanai–Tajimi excitation. The nonstationary responses including the nonstationary probability densities of the system responses and the statistical moments are obtained in semi‐analytical form. By applying the stochastic averaging method based on the generalized harmonic functions, the averaged Fokker–Planck–Kolmogorov(FPK) equation governing the nonstationary probability density of the amplitude is derived. Then, the solution of the FPK equation is approximately expressed by a series expansion in terms of a set of properly selected basis functions with time‐dependent coefficients. According to the Galerkin method, the time‐dependent coefficients are solved from a set of linear first‐order differential equations. Thus, the nonstationary probability densities of the amplitude and the state responses as well as the statistic moments of the amplitude are obtained. Finally, two types of the modulating functions, i.e. constant function and exponential function, are considered to give some semi‐analytical formulae. The proposed procedures are checked against the Monte Carlo simulation. The effects of the structure natural frequency and the intensity of the excitation as well as the ground stiffness on the system responses are discussed. It should be pointed out that the proposed method is good for broadband excitation and light damping. Copyright © 2011 John Wiley & Sons, Ltd.     

A non-reflecting boundary condition for the finite element modeling of infinite reservoir with layered sediment

The design of seismic resistant concrete gravity dam necessitates accurate determination of hydrodynamic pressure developed in the adjacent reservoir. The hydrodynamic pressure developed on structure is dependent on the physical characteristics of the boundaries surrounding the reservoir including reservoir bottom. The sedimentary material in the reservoir bottom absorbs energy at the bottom, which will affect the hydrodynamic pressure at the upstream face of the dam. The fundamental parameter characterizing the effect of absorption of hydrodynamic pressure waves at the reservoir bottom due to sediment is the reflection coefficient. The wave reflection coefficient is determined from parameters based on sediment layer thickness, its material properties and excitation frequencies. An analytical or a closed-form solution cannot account for the arbitrary geometry of the dam or reservoir bed profile. This problem can be efficiently tackled with finite element technique. The need for an accurate truncation boundary is felt to reduce the computational domain of the unbounded reservoir system. An efficient truncation boundary condition (TBC) which accounts for the reservoir bottom effect is proposed for the finite element analysis of infinite reservoir. The results show the efficiency of the proposed truncation boundary condition.     

Influence of time-domain dam–reservoir interaction on cracking of concrete gravity DAMS

Based on a non-linear dam-reservoir interaction model, a study investigating the earthquake response of concrete gravity dams is presented. For the propagation of cracks in unreinforced mass concrete, a discrete crack approach formulation based on the finite element method is applied. A special crack element is used to follow a fictitious crack in order to account for a zone of microcracks developing at the crack tip. The reservoir is modelled using the boundary element method. At a fictitious boundary dividing the irregular finite part of the reservoir from the regular infinite part, the loss of energy due to pressure waves moving away towards infinity is taken into account rigorously. Analyses are performed on the tallest non-overflow monolith of the Pine Flat Dam located in Kern County, California. The interaction of a dam, which may exhibit cracks in mass concrete, with a reservoir domain of arbitrary geometry extending to infinity is studied. Some main parameters are investigated. The importance of tools capable of handling the non-linear dam-reservoir interaction is emphasized.     

Systematic lumped‐parameter models for foundations based on polynomial‐fraction approximation

Based on the approximation by polynomial‐fraction, a series of systematic lumped‐parameter models are developed in this paper for efficiently representing the dynamic behaviour of unbounded soil. Concise formulation is first employed to represent the dynamic flexibility function of foundation with a ratio of two polynomials. By defining an appropriate quadratic error function, the optimal coefficients of the polynomials can be directly solved from a system of linear equations. Through performing partial‐fraction expansion on this polynomial‐fraction and designing two basic discrete‐element models corresponding to the partial fractions, systematic lumped‐parameter models can be conveniently established by connecting these basic units in series. Since the systematic lumped‐parameter models are configured without introducing any mass, the foundation input motion can be directly applied to these models for their applications to the analysis of seismic excitation. The effectiveness of these new models is strictly validated by successfully simulating a semi‐infinite bar on an elastic foundation. Subsequently, these models are applied for representing the dynamic stiffness functions for different types of foundation. Comparison of the new models with the other existing lumped‐parameter models is also made to illustrate their advantages in requiring fewer parameters and featuring a more systematic expansion. Copyright © 2002 John Wiley & Sons, Ltd.     

Analytical and numerical analysis of pressure drawdown in a poroelastic reservoir with complete overburden effect considered

Starting from an analytical reservoir model that incorporates full interaction with an elastic overburden, a new hybrid mathematical approach is developed by combining two numerical discretization methods. A tabular reservoir (petroleum reservoir or an aquifer) in an infinite or semi-infinite domain is viewed as a macroscopic displacement discontinuity, allowing use of the efficient displacement discontinuity mathematical method to calculate stresses and displacements that arise because of pressure changes. A 3-D finite element method using a poroelastic formulation is used to discretize the reservoir itself. By coupling the displacement discontinuity and finite element methods, a 3-D large-scale poroelastic reservoir can be simulated within an infinite or semi-infinite domain. The numerical model has been verified through comparison to known solutions, and some time-dependent pressure drawdown problems are analyzed. Results indicate that including the complete overburden (reservoir surroundings) response has a significant effect on pressure drawdown in a poroelastic reservoir during pumping, and should be incorporated in appropriate applications such as well test equations and subsidence analyses.