地震作用下水-轴对称柱体相互作用的子结构分析方法

针对水-轴对称柱体动力相互作用问题,提出了一种地震作用下水-结构相互作用的时域子结构分析方法。基于三维不可压缩水体的波动方程和边界条件,利用分离变量法将其转换为环向解析、竖向和径向数值的二维模型;基于比例边界有限元推导了截断边界处无限域水体的动力刚度方程,并将水体内域有限元方程和人工边界处的动水压力进行耦合,从而得到结构表面的动水压力方程;将轴对称柱体结构的有限元方程与动水压力方程耦合,从而得到水-轴对称柱体结构系统的时域有限元方程;数值算例验证该文提出的水-轴对称动力相互作用的子结构方法,结果表明:该文方法具有很高的精度和计算效率。通过对水中轴对称结构地震响应和自振频率的分析表明:地震动水压力对结构自振频率和动力响应的影响随水深的增加而增大。     

A new integral identity for potential polygonal domain problems described by parametric linear functions

The paper presents a modification of the classical integral identity for two-dimensional potential boundary value problems with a linear segments boundary. The modification consists of describing integral contours in the integral identity by means of parametric linear functions. As a result of the modification, it was possible to obtain a new integral identity, which included the geometry of the boundary in its subinterval functions. The identity can be used for finding solutions in polygonal domains under the condition that the solutions previously obtained on the boundary are arrived at by the new parametric integral equation system. The proposed method makes it possible to obtain solutions of domain problems with no need for the discretization of the boundary geometry. The effectiveness of the method is illustrated by three testing examples.     

Direct method of solution for general boundary value problem of the Laplace equation

A general boundary value problem for two-dimensional Laplace equation in the domain enclosed by a piecewise smooth curve is considered. The Dirichlet and the Neumann data are prescribed on respective parts of the boundary, while there is the second part of the boundary on which no boundary data are given. There is the third part of the boundary on which the Robin condition is prescribed. This problem of finding unknown values along the whole boundary is ill posed. In this sense we call our problem an inverse boundary value problem. In order for a solution to be identified the inverse problem is reformulated in terms of a variational problem, which is then recast into primary and adjoint boundary value problems of the Laplace equation in its conventional form. A direct method for numerical solution of the inverse boundary value problem using the boundary element method is presented. This method proposes a non-iterative and unified treatment of conventional boundary value problem, the Cauchy problem, and under- or over-determined problems.     

Expansion formulae for wave structure interaction problems with applications in hydroelasticity

Alternate derivations of the expansion formulae for wave structure interaction problems are obtained in case of water of infinite depth and utilized to analyze the hydroelastic behavior of large floating structures. Considering the boundary value problem associated with Laplace equation having higher order boundary condition on the horizontal boundary and a Dirichlet type boundary condition on the vertical boundary in a quarter plane, Fourier sine transform is applied in the horizontal direction to convert the problem to a Sturm-Liouville type boundary value problem associated with non-homogeneous ordinary differential equation (ODE) in the transformed variable. Finally, inverting the transformed functions and applying the regularity criterion of the transformed function, the required expansion formula is derived. The expansion formula thus derived is extended to deal with similar boundary value problems having Neumann type boundary condition. The expansion formulae are applied to (i) analyze oblique scattering of flexural gravity waves by an articulated floating elastic plate and (ii) study the effect of compression on the oblique scattering of flexural gravity waves by a line discontinuity in a large floating ice sheet in water of infinite depth, which find applications in marine technology and arctic engineering, respectively. The present derivations of the expansion formulae are very simple and straightforward and can be easily used to study a large class of problems in the area of fluids and structures in mathematical physics and engineering.     

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 solution of the Dirichlet initial boundary value problem for the heat equation in exterior 3-dimensional domains using integral equations

A numerical method for the Dirichlet initial boundary value problem for the heat equation in the exterior and unbounded region of a smooth closed simply connected 3-dimensional domain is proposed and investigated. This method is based on a combination of a Laguerre transformation with respect to the time variable and an integral equation approach in the spatial variables. Using the Laguerre transformation in time reduces the parabolic problem to a sequence of stationary elliptic problems which are solved by a boundary layer approach giving a sequence of boundary integral equations of the first kind to solve. Under the assumption that the boundary surface of the solution domain has a one-to-one mapping onto the unit sphere, these integral equations are transformed and rewritten over this sphere. The numerical discretisation and solution are obtained by a discrete projection method involving spherical harmonic functions. Numerical results are included.     

Boundary integral equation solution of moving boundary phase change problems

Boundary integral equation methods are presented for the solution of some two-dimensional phase change problems. Convection may enter through boundary conditions, but cannot be considered within phase boundaries. A general formulation based on space-time Green's functions is developed using the complete heat equation, followed by a simpler formulation using the Laplace equation. The latter is pursued and applied in detail. An elementary, noniterative system is constructed, featuring linear interpolation over elements on a polygonal boundary. Nodal values of the temperature gradient normal to a phase change boundary are produced directly in the numerical solution. The system performs well against basic analytical solutions, using these values in the interphase jump condition, with the simplest formulation of the surface normal at boundary vertices. Because the discretized surface changes automatically to fit the scale of the problem, the method appears to offer many of the advantages of moving mesh finite element methods. However, it only requires the manipulation of a surface mesh and solution for surface variables. In some applications, coarse meshes and very large time steps may be used, relative to those which would be required by fixed grid domain methods. Computations are also compared to original lab data, describing two-dimensional soil freezing with a time-dependent boundary condition. Agreement between simulated and measured histories is good.     

A dual-reciprocity boundary element method for a class of elliptic boundary value problems for non-homogeneous anisotropic media

A dual-reciprocity boundary element method is proposed for the numerical solution of a two-dimensional boundary value problem (BVP) governed by an elliptic partial differential equation with variable coefficients. The BVP under consideration has applications in a wide range of engineering problems of practical interest, such as in the calculation of antiplane stresses or temperature in non-homogeneous anisotropic media. The proposed numerical method is applied to solve specific test problems.     

The boundary element method applied to the analysis of two-dimensional nonlinear sloshing problems

This paper deals with an application of the boundary element method to the analysis of nonlinear sloshing problems, namely nonlinear oscillations of a liquid in a container subjected to forced oscillations. First, the problem is formulated mathematically as a nonlinear initial-boundary value problem by the use of a governing differential equation and boundary conditions, assuming the fluid to be inviscid and incompressible and the flow to be irrotational. Next, the governing equation (Laplace equation) and boundary conditions, except the dynamic boundary condition on the free surface, are transformed into an integral equation by employing the Galerkin method. Two dynamic boundary condition is reduced to a weighted residual equation by employing the Galerkin method. Two equations thus obtained are discretized by the use of the finite element method spacewise and the finite difference method timewise. Collocation method is employed for the discretization of the integral equation. Due to the nonlinearity of the problem, the incremental method is used for the numerical analysis. Numerical results obtained by the present boundary element method are compared with those obtained by the conventional finite element method and also with existing analytical solutions of the nonlinear theory. Good agreements are obtained, and this indicates the availability of the boundary element method as a numerical technique for nonlinear free surface fluid problems.     

Regularized MFS-Based Boundary Identification in Two-Dimensional Helmholtz-Type Equations

We study the stable numerical identification of an unknown portion of the boundary on which a given boundary condition is provided and additional Cauchy data are given on the remaining known portion of the boundary of a two-dimensional domain for problems governed by either the Helmholtz or the modified Helmholtz equation. This inverse geometric problem is solved using the method of fundamental solutions (MFS) in conjunction with the Tikhonov regularization method. The optimal value for the regularization parameter is chosen according to Hansen's L-curve criterion. The stability, convergence, accuracy and efficiency of the proposed method are investigated by considering several examples.     

粘弹性层状半空间上刚体的垂直振动

本文较深入地研究了层状粘弹性半空间上刚性圆形基础的垂直振动问题。首先作者采用独创的方法求解弹性动力学基本方程，并借助于Ｔｈｏｍｓｏｎ Ｈａｓｋｅｌ的传递矩阵方法得出层状地基表面位移的Ｈａｎｋｅｌ变换式。然后按弹性力学混合边值问题建立层状粘弹性地基上刚体垂直振动的对偶积分方程，并化为第二类Ｆｒｅｄｈｏｌｍ积分方程。文末给出了数值算例。     

Finite Deflection of Slender Cantilever with Predefined Load Application Locus using an Incremental Formulation

In this paper, a class of problems involving space constrained loading on thin beams with large deflections is considered. The loading is such that, the locus of the force application point moves along an arbitrarily predefined path, fixed in space. Both linear elastic as well as elastic-perfectly plastic materials are considered. A simplification is realized using the moment-curvature relationship directly. The governing equation obtained is highly non-linear owing to inclusion of both material and geometric non-linearity. A general algorithm is described to solve the governing equation using an incremental formulation coupled with Runge Kutta 4^th order initial value explicit solver. Additionally, the presented method is capable of handling unloading and reverse loading conditions. An example problem where the load application point locus is an inclined straight line is solved to demonstrate the performance of the method. It is found that, the force response due to the inclined locus is stiffer than the vertical locus. This response is akin to dry friction condition on a vertical locus case.     

Implicit formulation with the boundary element method for nonlinear radiation of water waves

An accurate and efficient numerical method is presented for the two-dimensional nonlinear radiation problem of water waves. The wave motion that occurs on water due to an oscillating body is described under the assumption of ideal fluid flow. The governing Laplace equation is effectively solved by utilizing the GMRES (Generalized Minimal RESidual) algorithm for the boundary element method (BEM) with quadratic approximation. The intersection or corner singularity in the mixed Dirichlet–Neumann problem is resolved by introducing discontinuous elements. The fully implicit trapezoidal rule is used to update solutions at new time-steps, by considering stability and accuracy. Traveling waves generated by the oscillating body are absorbed downstream by the damping zone technique. To avoid the numerical instability caused by the local gathering of grid points, the re-gridding technique is employed, so that all the grids on the free surface may be re-distributed with an equal distance between them. The nonlinear radiation force is evaluated by means of the acceleration potential. For a mixed Dirichlet–Neumann problem in a computational domain with a wavy top boundary, the present BEM yields numerical solutions for the quadratic rate of convergence with respect to the number of boundary elements. It is also demonstrated that the present time-marching and radiation condition work successfully for nonlinear radiation problems of water waves. The results obtained from this study concur reasonably well with other numerical computations.     

Boundary Integral Equation for Eddy-Current Problems with Dirichlet Boundary Condition

A boundary integral equation (BIE) is developed for the eddy-current (EC) problems with Dirichlet boundary condition by considering the difference between the field without cracks and the one with cracks. Once the field and its normal derivative are given for the structure in the absence of cracks, normal derivative of the scattered field on the surface can be calculated by solving this integral equation numerically. For infinite-domain problems, this equation is more efficient than the conventional BIE due to a smaller computational region needed. Four kinds of two-dimensional EC problems have been solved using this integral equation. The surface impedance for different cases is presented in this paper. Numerical results are compared with analytical solutions and published numerical results. There are good agreements between them. Also, this concept can be extended to three-dimensional problems with other boundary conditions.     

Relaxation procedures for an iterative MFS algorithm for two-dimensional steady-state isotropic heat conduction Cauchy problems

We investigate two algorithms involving the relaxation of either the given Dirichlet data (boundary temperatures) or the prescribed Neumann data (normal heat fluxes) on the over-specified boundary in the case of the alternating iterative algorithm of Kozlov et al. [26] applied to two-dimensional steady-state heat conduction Cauchy problems, i.e. Cauchy problems for the Laplace equation. The two mixed, well-posed and direct problems corresponding to each iteration of the numerical procedure are solved using a meshless method, namely the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method. For each direct problem considered, the optimal value of the regularization parameter is chosen according to the generalized cross-validation (GCV) criterion. The iterative MFS algorithms with relaxation are tested for Cauchy problems associated with the Laplace operator in various two-dimensional geometries to confirm the numerical convergence, stability, accuracy and computational efficiency of the method.     

Stable Boundary and Internal Data Reconstruction in Two-Dimensional Anisotropic Heat Conduction Cauchy Problems Using Relaxation Procedures for an Iterative MFS Algorithm

We investigate two algorithms involving the relaxation of either the given boundary temperatures (Dirichlet data) or the prescribed normal heat fluxes (Neumann data) on the over-specified boundary in the case of the iterative algorithm of Kozlov91 applied to Cauchy problems for two-dimensional steady-state anisotropic heat conduction (the Laplace-Beltrami equation). The two mixed, well-posed and direct problems corresponding to every iteration of the numerical procedure are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method. For each direct problem considered, the optimal value of the regularization parameter is chosen according to the generalized cross-validation (GCV) criterion. The iterative MFS algorithms with relaxation are tested for over-, equally and under-determined Cauchy problems associated with the steady-state anisotropic heat conduction in various two-dimensional geometries to confirm the numerical convergence, stability, accuracy and computational efficiency of the method.     

Indirect Trefftz method for boundary value problem of Poisson equation

Trefftz method is the boundary-type solution procedure using regular T-complete functions satisfying the governing equation. Until now, it has been mainly applied to numerical analyses of the problems governed with the homogeneous differential equations such as the two- and three-dimensional Laplace problems and the two-dimensional elastic problem without body forces. On the other hand, this paper describes the application of the indirect Trefftz method to the solution of the boundary value problems of the two-dimensional Poisson equation. Since the Poisson equation has an inhomogeneous term, it is generally difficult to determine the T-complete function satisfying the governing equation. In this paper, the inhomogeneous term containing an unknown function is approximated by a polynomial in the Cartesian coordinates to determine the particular solutions related to the inhomogeneous term. Then, the boundary value problem of the Poisson equation is transformed to that of the Laplace equation by using the particular solution. Once the boundary value problem of the Poisson equation is solved according to the ordinary Trefftz formulation, the solution of the boundary value problem of the Poisson equation is estimated from the solution of the Laplace equation and the particular solution. The unknown parameters included in the particular solution are determined by the iterative process. The present scheme is applied to some examples in order to examine the numerical properties.     

Large eddy simulation of turbulent flows in external flow field using three-step FEM–BEM model

An innovative computational model is presented for the large eddy simulation (LES) modeling of multi-dimensional unsteady turbulent flow problems in external flow field. Based on the LES principles, the model uses a pressure projection method to solve the Navier–Stokes equations in transient condition. The turbulent motion is simulated by Smagorinsky sub-grid scale (SGS) eddy viscosity model. The momentum equation of the flow motion is solved using a three-step finite element method (FEM). The external flow field is simulated using a boundary element method (BEM) by solving a pressure Poisson equation that assumes the pressure as zero at the infinity. Through extracting the boundary effects on a specified finite computational domain, the model is able to solve the infinite boundary value problems. The present model is used to simulate the flows past a two-dimensional square rib and a three-dimensional cube at high Reynolds number. The simulation results are found to be reasonable and comparable with other models available in the literature even for coarse meshes.     

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.     

Boundary element analysis of nonlinear water wave problems

A new boundary element technique based on Green's formula is applied to the analysis of nonlinear water wave problems. The problems are formulated mathematically as two-dimensional nonlinear initial-boundary value problems in terms of a velocity potential, assuming the fluid to be inviscid and incompressible, and the flow to be irrotational. In the present paper, two kinds of wave-making problems are analysed: (1) a tsunami generated by sea bed elevation; (2) generation, propagation and run-up on a vertical wall of a solitary wave. Numerical results obtained by the present method are compared with available experimental data and analytical solutions. Excellent agreements are obtained.