边界元法在环境声学中的应用

边界元法是边界积分方程的数值解法 ,是随着计算机技术的发展而出现的。建立声学边界积分方程分两种方法 :直接法与间接法。本文介绍了边界元法在环境声学中的应用 ,如声屏障和不同情况下道路周围的声场分布、复杂气象条件对声传播的影响的问题等。由于边界元法是半解析半数值解法。在解边界积分方程时会遇到解的存在与唯一性问题。     

薄壁杆件翘曲剪应力的边界元精确积分解法

用非连续边界元对薄壁杆件的约束扭转进行了分析,推导出了求解边界点二次翘曲函数值的边界积分方程,给出了边界积分方程数值求解时积分计算的精确表达式。数值算例表明:利用边界积分方程方法分析薄壁杆件的约束扭转问题时效率和精度高,同时采用精确积分可以有效的处理"边界层效应"问题。     

基于失效边界特性的体系可靠度计算方法

根据各失效模式验算点分布特性的不同,对体系失效边界的特性进行了补充分析。得到了多失效模式下边界不光滑程度较高的原因。另外分析了Gauss-Hermite积分方法求解体系失效概率的适用性问题。然后基于这些特性,提出了一种结合降维分解和梯形积分技术来计算体系失效概率的改进方法。实例分析表明该方法具有较高的精度和效率。     

边界元法分析功能梯度涂层材料

在常规边界元法中引入几乎奇异积分的解析算法,使边界元法可以分析涂层结构的强度.计算了在赫兹压力作用下,各向同性涂层和功能梯度涂层两种涂层结构中的Tresca应力分布,绘制了应力等值线图.计算发现使用各向同性涂层时,Tresca应力的最大值出现在涂层和基体交界面上,且在交界面上应力存在明显的不连续性.通过边界元法分析,发现采用功能梯度涂层,可以降低最大的Tresca应力值,削减交界面上的应力不连续性.     

多角形域上第一类边界积分方程的高精度配置法

本文提出了一种多角区域上的第一类边界积分方程的高精度算法离散之前采用特殊周期变换,消去边界积分方程未知函数在积分端点的奇异性,然后使用常元配置法求解.该方法在内点获得超收敛O(h3).此外,通过Richardson整体外推,可进一步提高内点解的精度.实际计算结果表明,该方法优于同一问题的Galerkin方法甚至机械求积法.     

基于分布源边界点法的近场声全息技术

分布源边界点法是一种新型的声辐射数值计算方法,通过分布在振动体内部一系列的特解源间接地构造出声源与声场之间的传递矩阵,从而有效地避开了采用边界元法求解声辐射问题时所存在的复杂的变量插值、奇异积分的处理、特征波数处解的非唯一性等问题.本文首先建立了分布源边界点法实现近场声全息变换的模型,并提出采用Tikhonov正则化方法稳定重建过程,抑制了重建误差的影响,随后通过实验研究进一步验证了文中方法的可行性和正确性.     

二维正交各向异性结构弹塑性问题的边界元分析

给出了二维正交各向异性结构弹塑性问题的边界元分析方法, 包括相应边界积分方程、内点应力公式、边界元求解格式以及弹塑性应力计算方法。在弹塑性分析中, 引入了Hill-Tsai 屈服准则, 采用初应力法和切向预测径向返回法确定实际应力状态。通过具体算例分析了二维正交各向异性结构的弹塑性应力和塑性区分布情况, 部分数值结果与已有结果进行了比较, 两者基本吻合。结果表明, 本文中给出的边界元法可以有效地用于求解二维正交各向异性结构的弹塑性问题。      

平面电磁弹性固体的虚边界元——最小二乘配点法

从电磁弹性固体平面问题的基本方程出发,依据弹性力学虚边界元法的基本思想,利用电磁弹性固体平面问题的基本解,提出了电磁弹性固体平面问题的虚边界元——最小二乘配点法。电磁弹性固体的虚边界元法在继承传统边界元法优点的同时,有效地避免了传统边界元法的边界积分奇异性的问题。由于仅在虚实边界选取配点,此方法不需要网格剖分,并且不用进行积分计算。最后给出了一些具体算例,并和已有的解析解进行了对比,结果表明提出的虚边界元方法有很高的精度。     

求解夹杂问题的有限部积分──边界元法

利用Kelvin解及有限部积分的概念和方法,导出求解含夹杂二维有限弹性体的超奇异积分方程,继而使用有限部积分与边界元结合的方法,为其建立了数值求解方法,即有限部积分与边界元法.最后计算了若干典型数值例子夹杂端部的应力强度因子.      

正交各向异性厚板振动分析的边界积分方程法

本文采用正交各向异性厚板静力问题的基本解作为边界积分方程的核函数,利用加权残数法建立了正交各向异性厚板振动分析的边界积分方程。文中详细地讨论了边界积分方程的数值处理过程并给出了若干数值算例以论证本文方法的正确性。      

The first-kind and the second-kind boundary integral equation systems for solution of frictionless contact problems

An original approach to the numerical solution of displacement boundary integral equation (BIE) and traction hypersingular boundary integral equation (HBIE) by the boundary element method (BEM) for contact problems is given. The main point is to show, how the contact conditions are used to formulate the first-kind and the second-kind BIE systems in the case of frictionless two-body elastic contact. The solution of the first-kind BIE is performed by symmetric Galerkin BEM; the second-kind BIE is solved by an appropriate collocation BEM. The contact problem in itself is solved by the method of subsequent approximations of contact region. Both forms of BIE system are compared in several numerical examples. This comparison is made for different kinds of contact problem. The major emphasis is put on the evaluation of contact pressure. The obtained results are compared with referenced numerical and with the analytical ones.     

Rayleigh wave correction for the BEM analysis of two‐dimensional elastodynamic problems in a half‐space

A simple, elegant approach is proposed to correct the error introduced by the truncation of the infinite boundary in the BEM modelling of two‐dimensional wave propagation problems in elastic half‐spaces. The proposed method exploits the knowledge of the far‐field asymptotic behaviour of the solution to adequately correct the BEM displacement system matrix for the truncated problem to account for the contribution of the omitted part of the boundary. The reciprocal theorem of elastodynamics is used for a convenient computation of this contribution involving the same boundary integrals that form the original BEM system. The method is formulated for a two‐dimensional homogeneous, isotropic, linearly elastic half‐space and its implementation in a frequency domain boundary element scheme is discussed in some detail. The formulation is then validated for a free Rayleigh pulse travelling on a half‐space and successfully tested for a benchmark problem with a known approximation to the analytical solution. Copyright © 2004 John Wiley & Sons, Ltd.     

ON THE REMOVAL OF RIGID BODY MOTIONS IN THE SOLUTION OF ELASTOSTATIC PROBLEMS BY DIRECT BEM

A theoretical and numerical study of the removal of rigid body motions in the solution of the boundary form of Somigliana identity and of the corresponding discretized linear system of the direct BEM is presented. This study is based on the Fredholm theory of linear operators and mechanical aspects of the problem. Various methods suitable for implementation in BEM codes are analyzed and relations between apparently different methods are shown. The relation between global equilibrium conditions and solvability of the discretized linear system of the direct BEM is discussed.     

Coupling BEM,FEM and analytic solutions in steady-state potential problems

Problems solved by using different steady-state solution techniques in adjacent subregions are discussed. The computational domain typically consists of two subregions, with a linear boundary value problem in one of them. BEM or analytical methods are used to solve the problem in this subregion. Static condensation of the off-interfacial degrees of freedom in this subdomain produces a linear set of equations linking nodal potentials and fluxes on the interface. This set of equations is generated by solving a sequence of boundary value problems in the linear subregion. Access to the source version of the software used to solve these boundary value problems is not required. Thus, the condensation can be accomplished using any commercial BEM code. The resulting set of equation is then treated as a boundary condition attached to the second subregion. In the latter, any numerical technique can be used and both linear and nonlinear problems may be considered. The paper addresses coupling of BEM and FEM, BEM and BEM and analytical solutions with BEM and FEM. Numerical examples are included.     

BEM for crack-hole problems in thermopiezoelectric materials

The boundary element formulation for analysing interaction between a hole and multiple cracks in piezoelectric materials is presented. Using Green's function for hole problems and variational principle, a boundary element model (BEM) for a 2-D thermopiezoelectric solid with cracks and holes has been developed and used to calculate stress intensity factors of the crack-hole problem. In BEM, the boundary condition on the hole circumference is satisfied a priori by Green's function, and is not involved in the boundary element equations. The method is applicable to multiple-crack problems in both finite and infinite solids. Numerical results for stress and electric displacement intensity factors at a particular crack tip in a crack-hole system of piezoelectric materials are presented to illustrate the application of the proposed formulation.     

正交各向异性平面问题边界元——边境元素法研究

在文献[1]中,本文作者研究了正交各向异性平面问题边界元素法的有关基本理论和计算公式,在上述工作的基础上,本文进一步研究各向异性平面问题边界邻域的应力分析。当采用边界元素法分析应力时,由于边界积分的奇异性,边界邻域应力的计算结果往往存在一定误差。为解决此问题,本文提出一个基于修正余能原理的所谓边境元素,包括四节点边境元素、八节点边境元素和三节点边境元素等。在边界元素法求解的基础上,进一步利用本文所述边境元素法,得到了非常满意的计算结果。      

Analytical integration of the moments in the diagonal form fast multipole boundary element method for 3-D acoustic wave problems

A diagonal form fast multipole boundary element method (BEM) is presented in this paper for solving 3-D acoustic wave problems based on the Burton-Miller boundary integral equation (BIE) formulation. Analytical expressions of the moments in the diagonal fast multipole BEM are derived for constant elements, which are shown to be more accurate, stable and efficient than those using direct numerical integration. Numerical examples show that using the analytical moments can reduce the CPU time by a lot as compared with that using the direct numerical integration. The percentage of CPU time reduction largely depends on the proportion of the time used for moments calculation to the overall solution time. Several examples are studied to investigate the effectiveness and efficiency of the developed diagonal fast multipole BEM as compared with earlier p³ fast multipole method BEM, including a scattering problem of a dolphin modeled with 404,422 boundary elements and a radiation problem of a train wheel track modeled with 257,972 elements. These realistic, large-scale BEM models clearly demonstrate the effectiveness, efficiency and potential of the developed diagonal form fast multipole BEM for solving large-scale acoustic wave problems.     

The importance of diagonal dominance in the iterative solution of equations generated from the boundary element method

The unsymmetric matrix equations generated from the boundary element method (BEM) can be solved iteratively, with convergence to the correct solution guaranteed, if the boundary element system of equations can be first transformed into an equivalent, diagonally dominant system. A transformation is presented which selectively annihilates terms in the coefficient matrix of the system Ax = b until an equivalent, row diagonally dominant system, if available, is obtained. The new, row diagonally dominant system is well suited for use with Jacobi and Gauss-Seidel point iterative equation solvers. The diagonal dominizing transformation presented in this work is not suitable for large systems of equations but is useful as a research tool for studying the importance of diagonal dominance in the iterative solution of equations generated from the BEM. A simple Laplacian problem is used to examine the structure of the BEM equations and to introduce the diagonal dominizing transformation. The importance of diagonal dominance is shown by comparing iterative convergence of positive-definite, symmetric positive-definite and diagonally dominant systems of BEM equations obtained from a plane strain elasticity problem.     

Comparison of boundary collocation methods for singular and non-singular axisymmetric heat transfer problems

This paper presents a computational study of some boundary collocation solution methods for the Laplace equation in cylindrical coordinates with axisymmetry. The methods compared are (i) the direct boundary element method (BEM), (ii) the method of fundamental solutions (MFS) with fixed sources and (iii) the Trefftz method. Relative accuracy of these methods are compared for two test problems. The first problem is a simple problem of heat transfer through a cylindrical rod which is a standard benchmark problem in this field. The second problem deals with heat transfer in silicon melt for Czochralski (CZ) process which involves a singularity in the boundary conditions at the corner of the crystal–melt interface. All the three methods indicated above are highly successful for the simple (first) problem with MFS and Trefftz being simpler to implement than the BEM. However, the Trefftz method was not effective for the second problem due to the boundary singularity and the MFS showed oscillations near the singularity point. Hence the use of higher order non-conforming elements with accurate Gauss–Kronrod integration schemes in the direct BEM method was investigated. It was found that the boundary singularity does not deteriorate the accuracy of the results if this improved numerical integration procedure is used in the direct BEM. Hence higher order elements with Gauss–Kronrod integration schemes can be used for the solution of many free interface problems encountered in crystal growth.     

Analytical integral algorithm in the BEM for orthotropic potential problems of thin bodies

There exist nearly singular integrals for boundary layer effect problem and thin body effect problem in the boundary element method (BEM). A new completely analytical integral algorithm is proposed and applied to evaluate the nearly singular integrals in the BEM for two-dimensional orthotropic potential problems of thin bodies. The completely analytical integral formulas are derived with integration by parts for the linear boundary interpolation. The present algorithm applies these analytical formulas to deal with the nearly singular integrals. The unknown potentials and fluxes at boundary nodes are firstly calculated accurately and then the physical quantities at the interior points are computed. Two benchmark numerical examples on heat conduction demonstrate that the present algorithm can handle thin structures with the thickness-to-length ratio down to 1.E−08. This indicates that the BEM is especially accurate and efficient for numerical analysis of thin body problems.