Arch dam–fluid interactions: By fem–bem and substructure concept

Arch dams can be conveniently analysed by the finite element method. For dam–fluid interaction problems, the fluid domain may be more conveniently handled by the boundary element method as a substructure first before connecting to the dam substructure. The added-mass matrix calculated from the fluid domain is symmetrized and lumped first so that the banded and symmetrical characteristics of the finite element method are retained. In the boundary element formulation, a mirror image method and quadratic elements are used for computational efficiency and accuracy. The strong singular terms are handled by using a solution which satisfies the governing equation and the free surface boundary condition. Infinite boundary conditions at the upstream of the reservoir can be reasonably approximated from the fundamental solution with accurate results, if the interior pressure distribution in the fluid domain is neglected. Numerical solutions on hydrodynamic pressure distribution and the natural frequencies of the dam–reservoir system with various water levels are obtained and compared with available analytical and experiment results.     

Quasi‐dual reciprocity boundary‐element method for incompressible flow: Application to the diffusive–advective equation

This work presents a new boundary‐element method formulation called quasi‐dual reciprocity formulation for heat transfer problems, considering diffusive and advective terms. The present approach has some characteristics similar to those of the so‐called dual‐reciprocity formulation; however, the mathematical developments of the quasi‐dual reciprocity approach reduces approximation errors due to global domain interpolation. Some one‐ and two‐dimensional examples are presented, the results being compared against those obtained from analytical and dual‐reciprocity formulations. The method convergence is evaluated through analyses where the mesh is successively refined for various Peclet numbers, in order to assess the effect of the advective term. Copyright © 2003 John Wiley & Sons, Ltd.     

Dynamic analysis of dam–reservoir-foundation interaction in time domain

In this paper, a time domain dynamic analysis of the dam–reservoir-foundation interaction problem is developed by coupling the dual reciprocity boundary element method (DRBEM) for the infinite reservoir and foundation domain and the finite element method for the finite dam domain. An efficient coupling procedure is formulated by using the substructuring method. Sharans boundary condition at the far end of the infinite fluid domain is implemented. To verify the proposed scheme, numerical examples are carried out and compared with available exact solutions and finite–finite element coupling results for the problem of the dam–reservoir interaction. Finally, a complete dam–reservoir-foundation interaction problem is solved and its solution is compared with previously published results.The author is thankful to the anonymous reviewer of this paper for his suggestions and comments, which improved considerably the present paper.     

A new boundary condition for energy radiation in covered reservoirs using BEM

In this paper, a new boundary condition accounting for energy radiation at the far end of covered reservoirs is proposed. Using boundary element modelling (BEM), the boundary condition is investigated through analysis of the hydrodynamic pressure within a two-dimensional ice-covered reservoir impounded by a gravity dam. The proposed boundary condition accounts for reservoir bottom absorption effects and the presence of an ice cover at infinity. Seismic excitation is introduced by subjecting the dam and the reservoir to a horizontal harmonic ground motion. The effectiveness of the proposed method and the accuracy of the boundary condition are examined through a parametric study. The boundary condition is shown to be accurate even when placed near the dam upstream face, and the results obtained are in excellent agreement when compared to those from a mathematical model developed by the authors in a previous work. Some fundamental aspects of hydrodynamic pressure within ice-covered reservoirs are also discussed.     

Coupling FEM and symmetric BEM for dynamic interaction of dam–reservoir systems

A numerical model coupling boundary and finite elements suitable for dynamic dam–reservoir interaction is presented herein. This model involves standard finite element idealization of the dam structure displacements and a new symmetric boundary element formulation of the unbounded reservoir domain leading to an equivalent symmetric stiffness matrix for the discretized pressure field. These two basic parts of the computation are directly coupled by imposing an equilibrium condition at the fluid–structure interface, then the resulting algebraic system is reduced by localizing the coupled terms in the global mass matrix such as usually achieved in the added-mass formulation. Finally, the performance and the accuracy of this model are examined by comparing its results to those obtained from three other numerical models.     

高拱坝气幕隔震效果的三维数值模拟

对锦屏一级拱坝进行了气幕隔震效果的三维数值模拟，其中气幕单元采用拉格朗日法中的位移格式，库水单元采用欧拉法中的压力场格式，坝体与气幕的交界面自动满足位移协调条件，气幕与库水交界面满足力的平衡条件。数值仿真结果表明：气幕显著降低了上游坝面的动水压力，气幕厚度越大，隔震效果越好。与无气幕时相比，气幕厚度为1m、3m、5m时，动水压力分别降低了77.1%、83.3%、87.7%，位移反应也相应降低。因此，气幕隔震技术在水工建设中有很好的推广前景。     

An improved boundary element formulation for wave propagation problems

In this paper the dual reciprocity formulation for scalar wave equations and elastodynamic problems developed by Nardini & Brebbia is extended to the problem of waves propagating in an infinite domain by applying the Sommerfeld's radiation condition on a suitable artificial boundary. The free surface condition of first order can also be taken into consideration. To validate the present scheme, some examples have been worked out and compared with analytical solutions.     

An efficient time‐domain perfectly matched layers formulation for elastodynamics on spherical domains

Many practical applications require the analysis of elastic wave propagation in a homogeneous isotropic media in an unbounded domain. One widely used approach for truncating the infinite domain is the so‐called method of perfectly matched layers (PMLs). Most existing PML formulations are developed for finite difference methods based on the first‐order velocity‐stress form of the elasticity equations, and they are not straight‐forward to implement using standard finite element methods (FEMs) on unstructured meshes. Some of the problems with these formulations include the application of boundary conditions in half‐space problems and in the treatment of edges and/or corners for time‐domain problems. Several PML formulations, which do work with FEMs have been proposed, although most of them still have some of these problems and/or they require a large number of auxiliary nodal history/memory variables. In this work, we develop a new PML formulation for time‐domain elastodynamics on a spherical domain, which reduces to a two‐dimensional formulation under the assumption of axisymmetry. Our formulation is well‐suited for implementation using FEMs, where it requires lower memory than existing formulations, and it allows for natural application of boundary conditions. We solve example problems on two‐dimensional and three‐dimensional domains using a high‐order discontinuous Galerkin (DG) discretization on unstructured meshes and explicit time‐stepping. We also study an approach for stabilization of the discrete equations, and we show several practical applications for quality factor predictions of micromechanical resonators along with verifying the accuracy and versatility of our formulation. Copyright © 2014 John Wiley & Sons, Ltd.     

Plane stress and plate bending coupling in BEM analysis of shallow shells

In this paper the shear deformable shallow shells are analysed by boundary element method. New boundary integral equations are derived utilizing the Betti's reciprocity principle and coupling boundary element formulation of shear deformable plate and two‐dimensional plane stress elasticity. Two techniques, direct integral method (DIM) and dual reciprocity method (DRM), are developed to transform domain integrals to boundary integrals. The force term is approximted by a set of radial basis functions. Several examples are presented to demonstrate the accuracy of the two methods. The accuracy of results obtained by using boundary element method are compared with exact solutions and the finite element method. Copyright © 2000 John Wiley & Sons, Ltd.     

A coupled BEM and arbitrary Lagrangian–Eulerian FEM model for the solution of two‐dimensional laminar flows in external flow fields

This paper describes a new computational model developed to solve two‐dimensional incompressible viscous flow problems in external flow fields. The model based on the Navier–Stokes equations in primitive variables is able to solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the pressure projection method. The external flow field is simulated using the boundary element method by solving a pressure Poisson equation that assumes the pressure as zero at the infinite boundary. The momentum equation of the flow motion is solved using the three‐step finite element method. The arbitrary Lagrangian–Eulerian method is incorporated into the model, to solve the moving boundary problems. The present model is applied to simulate various external flow problems like flow across circular cylinder, acceleration and deceleration of the circular cylinder moving in a still fluid and vibration of the circular cylinder induced by the vortex shedding. The simulation results are found to be very reasonable and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

Seismically induced, non-stationary hydrodynamic pressure in a dam-reservoir system

Stochastic seismic analysis of hydrodynamic pressure in a dam-reservoir system is presented in this paper. The analysis is conducted assuming infinite reservoir compressible fluid and modeling seismic acceleration as a normal zero-mean stochastic process obtained by Penzien filter. The non-homogeneous boundary conditions associated to the problem have been incorporated into the equation of pressure wave scattering in the form of a forcing function turning the non-homogeneous boundary value problem into an homogeneous one. Solution obtained via modal analysis in time-domain is coupled with the use of classical Itô stochastic differential calculus to characterize the stochastic hydrodynamic pressure field. Both cases of hydrodynamic pressure acting along the upstream face of the dam in presence of stationary and non-stationary seismic accelerations have been considered.     

High‐order Higdon‐like boundary conditions for exterior transient wave problems

Recently developed non‐reflecting boundary conditions are applied for exterior time‐dependent wave problems in unbounded domains. The linear time‐dependent wave equation, with or without a dispersive term, is considered in an infinite domain. The infinite domain is truncated via an artificial boundary ??, and a high‐order non‐reflecting boundary condition (NRBC) is imposed on ??. Then the problem is solved numerically in the finite domain bounded by ??. The new boundary scheme is based on a reformulation of the sequence of NRBCs proposed by Higdon. We consider here two reformulations: one that involves high‐order derivatives with a special discretization scheme, and another that does not involve any high derivatives beyond second order. The latter formulation is made possible by introducing special auxiliary variables on ??. In both formulations the new NRBCs can easily be used up to any desired order. They can be incorporated in a finite element or a finite difference scheme; in the present paper the latter is used. In contrast to previous papers using similar formulations, here the method is applied to a fully exterior two‐dimensional problem, with a rectangular boundary. Numerical examples in infinite domains are used to demonstrate the performance and advantages of the new method. In the auxiliary‐variable formulation long‐time corner instability is observed, that requires special treatment of the corners (not addressed in this paper). No such difficulties arise in the high‐derivative formulation. Published in 2005 by John Wiley & Sons, Ltd.     

Hydrodynamic pressure on an accelerating dam and criteria for cavitation

Summary The effect of finite reservoir on the hydrodynamic pressure due to horizontal as well as vertical ground excitations has been studied. It is found that for horizontal accelerations the hydrodynamic pressure force decreases as the size of the reservoir decreases. The effect of vertical acceleration on the pressure force on a dam is simply to adjust the hydrostatic pressure by replacing the gravitational constant by an effective gravitational acceleration and this is true for any arbitrary shapes of the reservoir. A simple criterion has been presented in this paper which would enable dam engineers to determine whether a given earthquake could cause cavitation at the dam-water interface or not.     

Non‐reflecting artificial boundaries for transient scalar wave propagation in a two‐dimensional infinite homogeneous layer

This paper presents an exact non‐reflecting boundary condition for dealing with transient scalar wave propagation problems in a two‐dimensional infinite homogeneous layer. In order to model the complicated geometry and material properties in the near field, two vertical artificial boundaries are considered in the infinite layer so as to truncate the infinite domain into a finite domain. This treatment requires the appropriate boundary conditions, which are often referred to as the artificial boundary conditions, to be applied on the truncated boundaries. Since the infinite extension direction is different for these two truncated vertical boundaries, namely one extends toward x →∞ and another extends toward x→‐ ∞, the non‐reflecting boundary condition needs to be derived on these two boundaries. Applying the variable separation method to the wave equation results in a reduction in spatial variables by one. The reduced wave equation, which is a time‐dependent partial differential equation with only one spatial variable, can be further changed into a linear first‐order ordinary differential equation by using both the operator splitting method and the modal radiation function concept simultaneously. As a result, the non‐reflecting artificial boundary condition can be obtained by solving the ordinary differential equation whose stability is ensured. Some numerical examples have demonstrated that the non‐reflecting boundary condition is of high accuracy in dealing with scalar wave propagation problems in infinite and semi‐infinite media. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical and spectral investigations of Trefftz infinite elements

The numerical and spectral performance of novel infinite elements for exterior problems of time‐harmonic acoustics are examined. The formulation is based on a functional which provides a general framework for domain‐based computation of exterior problems. Two prominent features simplify the task of discretization: the infinite elements mesh the interface only and need not match the finite elements on the interface. Various infinite element approximations for two‐dimensional configurations with circular interfaces are reviewed. Numerical results demonstrate the good performance of these schemes. A simple study points to the proper interpretation of spectral results for the formulation. The spectral properties of these infinite elements are examined with a view to the representation of physics and efficient numerical solution. Copyright © 1999 John Wiley & Sons, Ltd.     

Modal analysis of plates using the dual reciprocity boundary element method

This paper presents a new method for determining the natural frequencies and mode shapes for the free vibration of thin elastic plates using the boundary element and dual reciprocity methods. The solution to the plate's equation of motion is assumed to be of separable form. The problem is further simplified by using the fundamental solution of an infinite plate in the reciprocity theorem. Except for the inertia term, all domain integrals are transformed into boundary integrals using the reciprocity theorem. However, the inertia domain integral is evaluated in terms of the boundary nodes by using the dual reciprocity method. In this method, a set of interior points is selected and the deflection at these points is assumed to be a series of approximating functions. The reciprocity theorem is applied to reduce the domain integrals to a boundary integral. To evaluate the boundary integrals, the displacements and rotations are assumed to vary linearly along the boundary. The boundary integrals are discretized and evaluated numerically. The resulting matrix equations are significantly smaller than the finite element formulation for an equivalent problem. Mode shapes for the free vibration of circular and rectangular plates are obtained and compared with analytical and finite element results.     

Weak forms for modelling of rotationally symmetric,multilayered structures,including anisotropic poro‐elastic media

A weak form of the anisotropic Biot's equation represented in a cylindrical coordinate system using a spatial Fourier expansion in the circumferential direction is presented. The original three dimensional Cartesian anisotropic weak formulation is rewritten in an arbitrary orthogonal curvilinear basis. Introducing a cylindrical coordinate system and expanding the circumferential wave propagation in terms of orthogonal harmonic functions, the original, geometrically rotationally symmetric three dimensional boundary value problem, is decomposed into independent two‐dimensional problems, one for each harmonic function. Using a minimum number of dependent variables, pore pressure and frame displacement, a computationally efficient procedure for vibro‐acoustic finite element modelling of rotationally symmetric three‐dimensional multilayered structures including anisotropic porous elastic materials is thus obtained. By numerical simulations, this method is compared with, and the correctness is verified against, a full three‐dimensional Cartesian coordinate system finite element model. Copyright © 2012 John Wiley & Sons, Ltd.     

Diagonalization procedure for scaled boundary finite element method in modeling semi‐infinite reservoir with uniform cross‐section

To improve the ability of the scaled boundary finite element method (SBFEM) in the dynamic analysis of dam–reservoir interaction problems in the time domain, a diagonalization procedure was proposed, in which the SBFEM was used to model the reservoir with uniform cross‐section. First, SBFEM formulations in the full matrix form in the frequency and time domains were outlined to describe the semi‐infinite reservoir. No sediments and the reservoir bottom absorption were considered. Second, a generalized eigenproblem consisting of coefficient matrices of the SBFEM was constructed and analyzed to obtain corresponding eigenvalues and eigenvectors. Finally, using these eigenvalues and eigenvectors to normalize the SBFEM formulations yielded diagonal SBFEM formulations. A diagonal dynamic stiffness matrix and a diagonal dynamic mass matrix were derived. An efficient method was presented to evaluate them. In this method, no Riccati equation and Lyapunov equations needed solving and no Schur decomposition was required, which resulted in great computational costs saving. The correctness and efficiency of the diagonalization procedure were verified by numerical examples in the frequency and time domains, but the diagonalization procedure is only applicable for the SBFEM formulation whose scaling center is located at infinity. Copyright © 2009 John Wiley & Sons, Ltd.     

Explicit finite element perfectly matched layer for transient three‐dimensional elastic waves

The use of a perfectly matched layer (PML) model is an efficient approach toward the bounded‐domain modelling of wave propagation on unbounded domains. This paper formulates a three‐dimensional PML for elastic waves by building upon previous work by the author and implements it in a displacement‐based finite element setting. The novel contribution of this paper over the previous work is in making this finite element implementation suitable for explicit time integration, thus making it practicable for use in large‐scale three‐dimensional dynamic analyses. An efficient method of calculating the strain terms in the PML is developed in order to take advantage of the lack of the overhead of solving equations at each time step. The PML formulation is studied and validated first for a semi‐infinite bar and then for the classical soil–structure interaction problems of a square flexible footing on a (i) half‐space, (ii) layer on a half‐space and (iii) layer on a rigid base. Numerical results for these problems demonstrate that the PML models produce highly accurate results with small bounded domains and at low computational cost and that these models are long‐time stable, with critical time step sizes similar to those of corresponding fully elastic models. Copyright © 2008 John Wiley & Sons, Ltd.