Multi-region adaptive finite element-boundary element method for elasto-plastic analysis

In this paper, a multi-region adaptive finite element–boundary element (FEM-BEM) method for elasto-plastic analysis is presented. The method is suitable for several plasticity models, i.e. von-Mises, Tresca, Mohr–Coulomb, Drucker–Prager, Hill, and Hoffman yield criteria. The domain of the original problem is divided into the finite element and the boundary element sub-domains: the FEM is utilized in regions where plastic material behaviour is expected to develop, whereas the complementary linear elastic region is approximated using the symmetric Galerkin BEM. The adaptive method estimates the finite element and boundary element sub-domains, automatically generates the corresponding meshes and adapts the sub-domains, according to the state of computation.     

Inverse analysis FOR two-dimensional structures using the boundary element method

Inverse analysis is currently an important subject of study in several fields of science and engineering. The identification of physical and geometric parameters using experimental measurements is required in many applications. In this work a boundary element formulation to identify boundary and interface values as well as material properties is proposed. In particular the proposed formulation is dedicated to identifying material parameters when a cohesive crack model is assumed for 2D problems. A computer code is developed and implemented using the BEM multi-region technique and regularisation methods to perform the inverse analysis. Several examples are shown to demonstrate the efficiency of the proposed model.     

Lead field computation for the electrocardiographic inverse problem--finite elements versus boundary elements

In order to be able to solve the inverse problem of electrocardiography, the lead field matrix (transfer matrix) has to be calculated. The two methods applied for computing this matrix, which are compared in this study, are the boundary element method (BEM) and the finite element method (FEM). The performance of both methods using a spherical model was investigated. For a comparable discretization level, the BEM yields smaller relative errors compared to analytical solutions. The BEM needs less computation time, but a larger amount of memory. Inversely calculated myocardial activation times using either the FEM or BEM computed lead field matrices give similar activation time patterns. The FEM, however, is also capable of considering anisotropic conductivities. This property might have an impact for future development, when also individual myocardial fiber architecture can be considered in the inverse formulation.     

Fast multipole method based solution of electrostatic and magnetostatic field problems

Abstact Applications of boundary element methods (BEM) to the solution of static field problems in electrical engineering are considered in this paper. The choice of a suitable BEM formulation for electrostatics, steady current flow fields or magnetostatics is discussed from user's point of view. The dense BEM matrix is compressed with an enhanced fast multipole method (FMM) which combines well-known BEM techniques with the FMM approach. An adaptive grouping scheme for problem oriented meshes is presented along with a discussion on the influence of the mesh to the efficiency of the FMM. The computational costs of the FMM algorithm are analyzed for typical problems in practice. Finally, some electrostatic and magnetostatic numerical examples demonstrate the simple usability and the efficiency of the FMM. Communicated by: U. Langer     

CUDA架构下的三维弹性静力学边界元并行计算

针对传统边界元法计算量大、计算效率低的问题,以三维弹性静力学的边界元法为对象,将基于CUDA的GPU并行计算应用到其边界元计算中,提出了基于CUDA架构的GPU并行算法.该算法首先对不同类型的边界元系数积分进行并行性分析,描述了相关的GPU并行算法,然后阐述了边界元方程组的求解方法及其并行策略.实验结果表明,文中算法较传统算法具有显著的加速效果.     

层次式直接边界元计算VLSI三维互连电容

文中将Ａｐｐｅｌ处理多体问题的层次式算法思想实现于直接边界元法,用以计算ＶＬＳＩ三维互连寄生电容。直接边界积分方程同时含有边界上的电势与法向电场强度,能比间接边界元法更方便地处理多介质及有限介质结构,直接边界元法的层次式计算涉及对三种边界（强加边界、自然边界与介质交界面）及两种积分核（１／ｒ与１／ｒ＾３）的处理,显著区别于基于间接边界元法、仅处理强加边界与一种分核的层次式算法。文中以边界元的层次划     

Stability improvement to BEM/FEM coupling scheme for 2D scalar wave problems

The stability problem appeared in boundary element method/finite element method (BEM/FEM) coupling is discussed in this paper. As the response at time t_n+1 relates to the excitations and responses at all previous times, i.e. response history, BEM is easier to be unstable compared with FEM. The Newmark method for FEM is unconditionally stable, oscillations appeared at any time would decrease step by step as time goes by. But the oscillation history caused by FEM may be big enough to cause stability problems to the BEM scheme which although may be stable when used independently. A new procedure is used in this paper to reduce the oscillation history caused by FEM so that it will not cause stability problem to BEM scheme and further to the coupling BEM/FEM scheme. Numerical examples show that the proposed procedure can improve significantly to the stability of the coupling BEM/FEM scheme and cause little numerical damping.     

The simulation of 3D elastic scattering produced by thin rigid inclusions using the traction boundary element method

A mixed formulation that uses both the traction boundary element method (TBEM) and the boundary element method (BEM) is proposed to compute the three-dimensional (3D) propagation of elastic waves scattered by two-dimensional (2D) thin rigid inclusions. Although the conventional direct BEM has limitations when dealing with thin-body problems, this model overcomes that difficulty. It is formulated in the frequency domain and, taking into account the 2-1/2D configuration of the problem, can be expressed in terms of waves with varying wavenumbers in the zdirection, k_z. The elastic medium is homogeneous and unbounded and it should be noted that no restrictions are imposed on the geometry and orientation of the internal crack.     

分布式多区域多能微网群协同AGC算法

综合能源多区域协同是电网发展趋势, 而核心问题是采用何种方法对多区域进行协同. 本文基于Q ($\sigma $)融入了资格迹及双重Q学习, 提出一种面向多区域多能微网群的多智能体协同控制算法, 即DQ ($\sigma ,\lambda $), 避免传统强化学习动作探索值高估的同时, 来获取分布式多区域的协同. 通过对改进的IEEE两区域负荷频率控制模型及三区域多能微网群自动发电控制(Automatic generation control, AGC)模型仿真, 结果表明, 与传统方法相比, 所提算法具有快速收敛性和更优动态性能, 能获得分布式多区域多能微网群的协同.     

并行计算水下大尺度弹性壳体的低频声散射

有限元与边界元耦合模型是研究水下弹性壳体目标低频声散射常用的数值方法。应用该模型计算大尺度弹性目标的声散射时需要大量的计算时间与存储空间,采用并行数值的方式可以解决这一问题。首先并行计算生成有限元矩阵和边界元矩阵,然后应用并行化的广义极小残差(GMRES)迭代算法求解大型非对称线性方程组。详细叙述了并行GMRES(m)迭代算法的执行过程,并以球壳的声散射计算为例分析了迭代步数对算法收敛情况的影响。最后计算了Benchmark目标模型的低频散射声场,分析了其收发分置散射目标强度以及表面声场的分布。     

求解二维随机多区域声波散射问题的快速多极边界元方法

快速多极算法(FMM)是求解大尺度边界元问题的一种很有效的快速算法.应用快速多极算法求解二维随机多区域声散射问题的边界积分方程.首先给出了求解该问题的边界积分方程,进而给出快速多极算法求解的算法实现过程以及积分算子的相应多极展开、局部展开和相应系数的转化关系式.最后通过对数值例子的计算表明快速多极算法在求解随机多区域声散射问题时的可行性及高效性,其求解存储量和计算量都是O(N).     

Boundary element solution for the Cauchy problem in linear elasticity using singular value decomposition

In this paper the singular value decomposition (SVD), truncated at an optimal number, is analysed for obtaining approximate solutions to ill-conditioned linear algebraic systems of equations which arise from the boundary element method (BEM) discretisation of an ill-posed boundary value problem in linear elasticity. The regularisation parameter, namely the optimal truncation number, is chosen according to the discrepancy principle. The numerical results obtained confirm that the SVD+BEM produces a convergent and stable numerical solution with respect to decreasing the mesh size discretisation and the amount of noise added into the input data.     

NiHu: An open source C++ BEM library

This paper introduces NiHu, a C++ template library for boundary element methods (BEM). The library is capable of computing the coefficients of discretised boundary integral operators in a generic way with arbitrarily defined kernels and function spaces. NiHu’s template core defines the workflow of a general BEM algorithm independent of the specific application. The core provides expressive syntax, based on the operator notation of the BEM, reflecting the mathematics behind boundary elements in the C++ source code. The customisable Component library contains elements specific to particular applications such as different numerical integration techniques and regularisation methods. The library can be used for creating a standalone C++ application using external open source libraries, or compiling a Matlab toolbox through the MEX interface. By massively exploiting C++ template metaprogramming, NiHu generates optimised codes for specific applications, including heterogeneous problems. The paper introduces the main concepts of the novel development, demonstrates its versatility and flexibility and compares the implementation’s performance to that of other open source projects.     

A new fast multipole boundary element method for two dimensional acoustic problems

A new fast multipole boundary element method (BEM) is presented in this paper for solving large-scale two dimensional (2D) acoustic problems based on the improved Burton–Miller formulation. This algorithm has several important improvements. The fast multipole BEM employs the improved Burton–Miller formulation, and successfully overcomes the non-uniqueness difficulty associated with the conventional BEM for exterior acoustic problems. The improved Burton–Miller formulation contains only weakly singular integrals, and avoids the numerical difficulties associated to the evaluation of the hypersingular integral, it leads to the numerical implementations more efficient and straightforward. Furthermore, the fast multipole method (FMM) and the approximate inverse preconditioned generalized minimum residual method (GMRES) iterative solver are adopted to greatly improve the overall computational efficiency. The numerical examples with Neumann boundary conditions are presented that clearly demonstrate the accuracy and efficiency of the developed fast multipole BEM for solving large-scale 2D acoustic problems in a wide range of frequencies.     

Analytical integral algorithm applied to boundary layer effect and thin body effect in BEM for anisotropic potential problems

A new completely analytical integral algorithm is proposed and applied to the evaluation of nearly singular integrals in boundary element method (BEM) for two-dimensional anisotropic potential problems. The boundary layer effect and thin body effect are dealt with. The completely analytical integral formulas are suitable for the linear and non-isoparametric quadratic elements. The present algorithm applies the analytical formulas to treat nearly singular integrals. The potentials and fluxes at the interior points very close to boundary are evaluated. The unknown potentials and fluxes at boundary nodes for thin body problems with the thickness-to-length ratios from 1E−1 to 1E−8 are accurately calculated by the present algorithm. Numerical examples on heat conduction demonstrate that the present algorithm can effectively handle nearly singular integrals occurring in boundary layer effect and thin body effect in BEM. Furthermore, the present linear BEM is especially accurate and efficient for the numerical analysis of thin body problems.     

A combined boundary element and finite element solution of slab and slab-on-girder bridges

The boundary element method (BEM) has been shown by many researchers to be an efficient numerical tool with which to analyse various engineering structures. In particular, the method has been extensively applied to a large number of plate bending problems. The distinct advantage of the method is in the reduction in dimensionality of the problem and, as a consequence of this, both computer time and data preparation are significantly reduced. In this paper, the suitability of a new method combining the advantages of both the BEM and the finite element method (FEM) is studied. The method is first applied to investigate the conventional plate bending problems. After the validity of the method is established it is then extended to analyse the more complicated problems of slabs and slab-on-girder bridges. Through a series of slab-on-girder bridge calculations, it is demonstrated that the method is not only accurate and fast converging but its ease of application and data preparation is not attainable by the ordinary FEM.     

Symmetric BEM formulations for 2-D steady-state and transit potential problems

This paper presents a new concept for symmetric boundary element method (SBEM) applicable to 2-D steady-state and transit potential problems. Two kinds of SBEM formulations are derived. Symmetry is obtained simply through matrix manipulation, and no hypersingularity appears. Therefore, SBEM is much easier than the traditional symmetric Galerkin BEM. Compared with the traditional asymmetric BEM, the present SBEM can reduce the computational cost for time domain problems only. However, when applied to BEM/FEM coupling procedure, SBEM can reduce the computational cost for both steady-state and time domain problems. Three numerical examples are included to illustrate the effectiveness and accuracy of the present formulations.     

Analysis of three-dimensional anisotropic heat conduction problems on thin domains using an advanced boundary element method

In this paper, an advanced boundary element method (BEM) is developed for solving three-dimensional (3D) anisotropic heat conduction problems in thin-walled structures. The troublesome nearly singular integrals, which are crucial in the applications of the BEM to thin structures, are calculated efficiently by using a nonlinear coordinate transformation method. For the test problems studied, promising BEM results with only a small number of boundary elements have been obtained when the thickness of the structure is in the orders of micro-scales (10^?6), which is sufficient for modeling most thin-walled structures as used in, for example, smart materials and thin layered coating systems. The advantages, disadvantages as well as potential applications of the proposed method, as compared with the finite element method (FEM), are also discussed.     

A semi-analytical algorithm for the evaluation of the nearly singular integrals in three-dimensional boundary element methods

The nearly singular integrals occur in the boundary integral equations when the source point is close to an integration element (as compared to its size) but not on the element. In this paper, the concept of a relative distance from a source point to the boundary element is introduced to describe possible influence of the singularity of the integrals. Then a semi-analytical algorithm is proposed for evaluating the nearly strongly singular and hypersingular integrals in the three-dimensional BEM. By using integration by parts, the nearly singular surface integrals on the elements are transformed to a series of line integrals along the contour of the element. The singular behavior, which appears as factor, is separated from remaining regular integrals. Consequently standard numerical quadrature can provide very accurate evaluation of the resulting line integrals. The semi-analytical algorithm is applied to analyzing the three-dimensional elasticity problems, such as very thin-walled structures. Meanwhile, the displacements and stresses at the interior points very close to its bounding surface are also determined efficiently. The results of the numerical investigation demonstrate the accuracy and effectiveness of the algorithm.