Galerkin Boundary Element Method for Direct Boundary Integral Eqution for Dirichlet Problem of 2-D Laplace Eqation

The direct boundary integral equation of two-dimensional Laplace equation for Dirichlet problem is(con-sidered).It is deduced by Green's formula and the fundamental solution.The most-used numerical method for solving(direct) boundary integral equation is collocation method,and seldom have been used the Galerkin scheme in this case.The direct boundary integral eqution is changed into the variational eqution.Using linear element,it is solved by Galerkin boundary method.In the variational eqution double integrations shall be carried out.The paper presents the analytical formula to calculate the inner integration and the Gaussian quadrature is used for the outer integration. The numerical experimentation proved thefaesibility and the efficiency.     

Garlerkin Boundary Element Method with Hyper Singular integral kernel

A Galerkin Boundary Elements was applied to solve the first kind of integral equation with hyper-singularity, which can be deduced from the direct boundary integral formula for the Neumann problem of Laplace equation. The concept of integration by parts in the sense of distributions was used. When boundary rotation is introduced, the two order derivatives of singular kernel are shifted to the boundary rotation of unknown function in the Galerkin variational formulation. While linear boundary elements are used for 2-dimensional problems, the boundary rotation on each element can be discretized into a constant vector, so that the integration can be performed in a simple way and the difficulty of numerical calculation for hyper-singularity is overcome. The results of numerical examples demonstrate that the scheme presented is practical and effective.     

A BOUNDARY ELEMENT METHOD FOR THE NEUMANN PROBLEM OF HELMHOLTZ EQUATION BY USE OF A DOUBLE LAYER POTENTIAL IN R~3

This paper presents a new numerical solution for Neumann problem of Helmholtz equation in R~3. The expression of the solution for this problem is obtained by use of a double layer potential and it leads to a Fredholm boundary integral equation of the first kind. Then, the existence and unicity of the integral equation which is equivalent to the boundary value problem are obtained in a suitable Sobolev space. Finally, a variational form which is equivalent to the integral equation is applied to the construction of a finite element method and the error estimate is given.     

Study on the Nonlinear Vibration System with Slow-changing Parameters

In this paper a posteriori error estimates for Galerkin approximation of general operator equations is firstly presented in the framework of Sobolev spaces. Then a practical posteriori error estimates formula for the adaptive boundary element method solving the acoustic scattering problem with a finite plane screen is obtained by triangulations. The mathematical model of this problem is the three dimensional Neumann boundary value problem of Helmholtz equation with finite plane boundary.     

Exact Integration of Constant Element of Elastomer in Boundary Element Method

The boundary integral in Boundary Element Method affects the precision and the speed of the method. If the boundary integral with constant element, the nonsingular integrals are popularly calculated by the Gauss numerical integral, and the singular integrals are popularly calculated by the analytical integral. This paper presents an alternative way with Gauss formula to transform the double integral in elastic problem on 3-d into the linear integrals on the boundary of each subdomains, so that all the singular integrals and nonsingular integrals are calculated by analytical method. The example indicates that this method makes the precision and the speed of BEM improve.     

A BOUNDARY ELEMENT METHOD FOR SOLVING DIRICHLET PROBLEM OF THE HELMHOLTZ EQUATION IN R~2

This paper presents a boundary element method for solving Dirichlet bou-ndary value problem of the Helmholtz equation in R~2.First,the existence andthe uniqueness of an extended solution for the problem are obtained.Then,thesolution is expressed in terms of simple layer potentials,and this expression,which is suitable to the interior as well as the exterior problem,leads to aboundary integral equation of the first kind.Finally,a finite element approachis applied to solve a variational form which is equivalent to the boundaryintegral equation.     

Galerkin Boundary Element Method for Two-dimension Laplace Equation of Neumman Condition

The authors apply the Galenkin variational equation to solve the integral equation with hyper singularity, which can be deduced from the double layer solution for Neumann problem of Laplace equation. The scheme of partial integration in the sense of distributions is introduced to reduce the hyper singularity integral into a weak one with the boundary rotation of unknown function. The numerical implementation with linear boundary elements is presented. The numerical examples illustrate the feasibility and efficiency of the method.     

Finite Element Method of 3D Stable Heat Conduction of Steel Fluid in Continuous Casting Moulds

The heat transfer of steel fluid in continuous casting mould is a stable process and can be depicted with three-dimensional stable heat conduction equation depending on tension speed. The corresponding finite element equation, including the first, second and third boundary conditions, is deduced out with Galerkin residual method. The coded FEM program is used to analyze the temperature distribution of Q235 steel in continuous billet casting mould. The method proposed is a foundation of thermo-mechanical coupled analysis for the formation of solidified shells and stress in the shells in continuous casting.     

三维轴对称功能梯度材料瞬态热传导问题的自然单元法

为了更有效地求解三维轴对称功能梯度材料瞬态热传导问题,对无网格自然单元法应用于此类问题进行了研究,并发展了相应的计算方法。基于几何形状和边界条件的轴对称性,三维的轴对称问题可降为二维平面问题。为了简化本质边界条件的施加,轴对称面上的温度场采用自然邻近插值进行离散。功能梯度材料特性的变化由高斯点的材料参数进行模拟。时间域上,采用传统的两点差分法进行离散求解,进而得到瞬态温度场的响应。数值算例结果表明,提出的方法是行之有效的,理论及方法不仅拓展了自然单元法的应用范围,而且对三维轴对称瞬态热传导分析具有普遍意义。     

Application of Symmetric Galerkin Boundary Element Method in Limit Analysis of Plane Problems

Based on Melan's theorem, the symmetric Galerkin boundary element method (SGBEM) is used to discretize two dimensional structures and a computational formulation of structural limit analysis is established. The self equilibrium stress field is constructed by linear combination of several basic vectors, which are the stress differences between different iteration steps at the same increment using the traditional elastoplastic incremental method. Then the complex method is used to solve the nonlinear programming directly, so that the lower bound load multiplier of two dimensional structures is obtained. The validation of the present method has been confirmed by some numerical examples.     

A NEW BOUNDARY INTEGRAL EQUATION METHOD FOR SOLVING EXTERIOR BOUNDARY VALUE PROBLEMS OF THREE-DIMENSIONAL HELMHOLTZ EQUATION

This paper pressnts a new boundary integral equation method for solving exteri-or boundary value problems of three-dimensional Heimheltz equation by using the multiple reciproc-ity method.Firstly,integral representations of the solution in an exterior domain as well as on itsboundary,which have the peculiarity that integral kernels are infin ite seriesea developed from thenormal fundamental solution of Laplace equation and independent of the wavenumber,are given andproved under the Dirichlet condition.Then,based on the representation of the solution on the bound-ary,boundary integral equations for solving the Dirichlet and the Neumann boundary value prob-lems are obtained,and remarks for some problems concerned with solving these integral equationsnumerically are made.Finally, the advantages of the proposed method,as compared with the conven-tional boundary element methods,are summarized.     

Analyzing the StiffnesS of Foundation Plate By BEM

This.paper derives the basic solution of elastic half-space foundation plate byHankel transform,Then its boundary integral equation of distortion is estabished,The paper discuss-es bou ndary integral equation on plate of fixxing and/or simple support boundary.Then a systemof algebraic equations for the boundary values are established.A numerical example is given , its re-sult shows that this Method is correct.     

Exact Integration of Biharmonic Equation In Boundary Element Method

The boundary integral in Boundary Element Method effects the precision and the speed of the method. The boundary integrals are composed of the normal integrals and singular integrals. The normal integrals are popularly calculated by exact integral, and the singular integrals by the Gauss numerical integral. The singular integrals are low in precision when the source points approach the element. This paper presents an alternative way to transform the double integral in Biharmonic Equation on 3-d into the linear integrals on the boundary of each subdomain, so that all the singular integrals and nonsingular integrals are calculated by analytical method. It makes the precision and the speed of BEM improve.     

Boundary Element Methods with Overlapping Domain Decomposition for Parabolic Problem

Boundary element method is a numerical method for solving partial differential equations. There are several formulations of boundary element method (BEM) applied to solve a parabolic differential equation.The approach,which employs time- dependent fundamental solution,allows longer time steps in time integration than other approaches,and this can cut down on time for computer implementation with high precision.Domain decomposition method,which decompose the domain that a given problem is to be solved into subdomains,has the advantages of reducing the large problem into smaller ones and reducing the complex problem into simpler ones,and allows parallel computing.An overlapping domain decomposition method is applied combining a boundary element formulation with time-dependent fundamental solution to solve a diffusion equation. Firstly, by domain decomposition, the problem divided into two problems on subdomains, and then the initial-Boundary problems are solved by boundry element method on each subdomain.Some numerical examples are presented to illustrate feasibility and efficiency of the method. The numerical experiments show that the convergence rate of the method is dependent with the overlapping degree of the subdomains.     

THE THREE-DIMENSIONAL CALCULATING MODEL OF INVERSE PROBLEM IN ECG AND A SPECIAL BOUNDRYVALUE PROBTEM

In this paper, the mathematical model and the three-dimensional finite clement formulation of inverse problem in electrocardiography (ECG) under a general inhomogcneous and anisotropic configuration of the torso conductor is focused on. The following points are especially considered: 1. The mathamatical model is described by a differential equation instead of the Fredholm integral equation. 2. A kind of special boundary value problem is defined as ill-posed boundary value problem. 3. A group of finite element formulation, in which the method that treats anisotropic medium, is built up. 4. The calculating method of ill-posed boundary value problem is discussed.     

Complex Variable Function Solution of Stress and Displacement of Surrounding Rock Buried Deep Horseshoe-Shaped Tunnel Excavation

An analytical solution was presented for the stresses and displacements around horseshoe-shaped tunnel using Cauchy integral method and Optimization theory in an elastic half-plane. Since the depth of a tunnel was larger than the size of the tunnel, gravity loads were simplified as uniform normal loads along far-field boundary. Finite element model was used to verify the accuracy of analytical solutions. Comparison between the results of numerical analysis by finite element method and those from the closed solutions indicates that the closed solution is reliable and applicable for the stress and displacement field around horseshoe-shaped tunnel tunnels at great depth.     

Study on Strengthening and Toughening Mechanism of Nb V Bearing High Strength Steels

It is an important subject to solve numerically the problems of anisotropic unsteady seepage flow.This paper presents a numerical method for computing the anisotropic_orthotropic unsteady seepage flow with BEM.The paper introduces the so_called boundary_only technique to solve time dependent problem.The Volume domain integral is voided by using the Green formula again.Several numerical examples are presented.     

Remarks for Various Boundary IntegralEquations Reduced from the ExteriorNeumann Problem of Helmholtz Equation

This paper presents a discussion on various boundary integral equations reduced from the exterior Neumann problem of Helmholtz equation.The author analyses how the famous difficulty that some equations have no unique solution when the wave number k is an eigenvalue of an interior problem is arised in the course of reducing these equations from Helmholtz representations,and proposes a method of overcoming the difficulty,that is,introducing a direct boundary integral equation which has unique solution for all wave numbers k and is equivalent to the original boundary value problem.Besides,advantages and shortcomings for these integral equations are estimated respectively.     

Numerical Analysis for Torsion of Prismatic Bars with Sectorial Cross-section

In mechanical design and application, in order to calculate torsion angle and shear stresses of prismatic bars with seetorial cross-section undergoing extemal couples, the numerical method is presented. The torsion equation is nonhomogeneous partial differential equation. First, using the method of separation of variables, torsion stress function is acquired in polar coordinate. Then, the method of boundary collocation is improved to calculate the undetermined parameters. Finally,approximate numerical solutions of stress function and shear stresses in cross-section are obtained. It is given the several calculation results of shear stresses of prismatic bars with different vertex angles. These results show that the method has some precision and application feasibility in engineering design. The method of separation of variables is combinied with the method of boundary collocation simplified calculation process.     

Effects of Thorax Muscle Fiber on Body Surface Potentials

A numerical method calculating the forward problem of electrocardiogram (ECG) has been presented. By using this method, the torso part including anisotropic muscle layers has been discretized by means of the finite element method (FEM). The rest of torso have been divided into the surface elements by the boundary element method (BEM). The effects of the skelltel muscle layers over thorax wall, have been discussed in terms of a three-dimensional torso models.