Helmholtz声学边界积分方程中奇异积分的计算

提出了一种非等参单元的四边形坐标变换,它将积分的曲面单元映射为另一四边形单元,通过两次坐标变换引入的雅可比行列式可以消除Helmholtz声学边界积分方程中的弱奇异型O(1／r))积分．而且利用δr／δn以及坐标变换可以同时消除坐标变换无法消除的Cauchy型(O(1／r^2))奇异积分,并给出了消除奇异性的详细证明．该方法给Helmholtz声学边界积分方程中的弱奇异积分与Cauchy奇异积分的计算以及编程提供了极大便利。     

结构振动辐射声场的预估——边界积分方程中奇异积分的间接处理

本文讨论了利用边界积分方程和边界元技术计算结构的稳态外辐射声场的方法,同时,对边界元方法所固有的奇异数值积分提出了一种简单方便的间接处理方法。计算实例证明所编计算程序和奇异积分处理方法是成功的。利用该程序在已知结构表面振速分布的条件下,可以求出该结构在自由声场中的声功率、表面辐射效率以及声场中任意点的声压值和相位。对一个实际钢质空心封闭圆筒作了计算与实测的比较,结果显示了该方法可应用于实际结构或机器的前景。     

用边界元法计算声辐射时高次奇异积分的处理方法

边界元法应用于计算辐射声场时，由于奇积分的存在，会影响以计算结果的精度，本文描述了处理带有１／ｒ奇异积分和１／ｒ＾２二次奇异积分处理方法，包括数学证明和数值积分方法，计算结果表明这种方法能提高精度。     

用样条边界元计算振动体的三维稳态声辐射

本文采用三次Ｂ样条边界单元计算振动体的三维稳态声辐射。实际计算表明：采用样条边界元可获得较好的数值计算结果。此外，本文在计算声场内点声压的Ｈｅｌｍｈｏｌｔｚ边界积分公式的基础上，导出了计算声场内点质点振速和声强的边界积分公式。文中还给出了应用本文方法计算的算例结果。     

粘弹性分析的复变边界积分方程法

本文应用弹性－粘弹性对应原理提出了基于复位势基本解的二维粘弹性分析的复变边界积分方程方法，给出了所有基本关系式，编制了相应的计算程序并给出计算实例，与已有工作相比，本方法具有公式统一，程序简洁通用，边界单元数少和效率高等特点。     

边界元法分析中的边界切向应力差分算法

在用边界元法作弹性应力分析中,不能直接计算出弹性体边界切向应力。本文在边界元法分析的基础上,用差分法计算边界切向应力。推导出常边界单元情况下边界切向应力的差分公式。计算表明文中所述方法是可行的,并且简单实用。所研究的方法和公式也适用于高次边界单元的边界切向应力的计算。     

组装结构中的混合边界积分方程模型

本文基于混合应力概念，对结构组装中的荷载传递问题，运用混合边界积分方程方法，建立了比较简的力学模型方程系统，给出了组合方程的离散形式，标定了各个量值的力学含义及数学求解途径，指出了研究结果的意义。     

边界积分方程的一种高阶估算技术

基于边界积分方程技术讨论了边界函数未知时的数值估算方法，推导出了均匀各向同性体普遍适用的ＢＩＥＭ公式，给出了积分方程边界值估算法技术的离散格式及数值计算形式，最后分析了矩形及三角形板数值做算例子。     

模态边界元法

本文介绍了用有限区域弹性体代替传统的无限区域弹性体建立边界方程基本解的思想。文中以有限弹性体模态的线性组合表示稳态基本解,建立了积分与激励频率无关的稳态边界元的直接表达式和间接表达式。其次文中提出了求解弹性系统固有频率和振型的设想,导出了边界元法特征方程。最后,文中以一维固有特性问题为例对原理、方法的正确性和精确度进行了检验,得到了不同边界条件下的固有频率和振型。文中所给数据表明与理论值吻合良好。     

Improved numerical method for the Traction Boundary Integral Equation by application of Stokes' theorem

This paper concerns the direct numerical evaluation of singular integrals arising in Boundary Integral Equations for displacement (BIE) and displacement gradients (BIDE), and the formulation of a Traction Boundary Integral Equation (TBIE) for solving general elastostatic crack problems. Subject to certain continuity conditions concerning displacements and tractions at the source point, singular integrals in the BIE and the BIDE corresponding to coefficients of displacement and displacement gradients at the source point are shown to be of a form that allows application of Stokes' theorem. All the singular integrals in 3-D BIE and BIDE are reduced to non-singular line integrals, and those in 2-D BIE and BIDE are evaluated in closed form. Remaining terms involve regular integrals, and no references to Cauchy or Hadamard principal values are required. Continuous isoparametric interpolations used on continuous elements local to the source point are modified to include unique displacement gradients at the source point which are compatible with all local tractions. The resulting numerical BIDE is valid for source points located arbitrarily on the boundary, including corners, and a procedure is given for constructing a TBIE from the BIDE. Some example solutions obtained using the present numerical method for the TBIE in 2-D and 3-D are presented. © British Crown Copyright 1997/DERA.     

Boundary integral equation for tangential derivative of flux in Laplace and Helmholtz equations

In this paper, the boundary integral equations (BIEs) for the tangential derivative of flux in Laplace and Helmholtz equations are presented. These integral representations can be used in order to solve several problems in the boundary element method (BEM): cubic solutions including degrees of freedom in flux's tangential derivative value (Hermitian interpolation), nodal sensitivity, analytic gradients in optimization problems, or tangential derivative evaluation in problems that require the computation of such variable (elasticity problems in BEM). The analysis has been developed for 2D formulation. Kernels for tangential derivative of flux lead to high‐order singularities (O(1/r³)). The limit to the boundary analysis has been carried out. Based on this analysis, regularization formulae have been obtained in order to use such BIE in numerical codes. A set of numerical benchmarks have been carried out in order to validate theoretical and practical aspects, by considering known analytic solutions for the test problems. The results show that the tangential BIEs have been properly developed and implemented. Copyright © 2005 John Wiley & Sons, Ltd.     

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

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

Comparison of regularization methods for solving the Cauchy problem associated with the Helmholtz equation

In this paper, several boundary element regularization methods, such as iterative, conjugate gradient, Tikhonov regularization and singular value decomposition methods, for solving the Cauchy problem associated to the Helmholtz equation are developed and compared. Regularizing stopping criteria are developed and the convergence, as well as the stability, of the numerical methods proposed are analysed. The Cauchy problem for the Helmholtz equation can be regularized by various methods, such as the general regularization methods presented in this paper, but more accurate results are obtained by classical methods, such as the singular value decomposition and the Tikhonov regularization methods. Copyright © 2004 John Wiley & Sons, Ltd.     

Boundary knot method for some inverse problems associated with the Helmholtz equation

The boundary knot method is an inherently meshless, integration‐free, boundary‐type, radial basis function collocation technique for the solution of partial differential equations. In this paper, the method is applied to the solution of some inverse problems for the Helmholtz equation, including the highly ill‐posed Cauchy problem. Since the resulting matrix equation is badly ill‐conditioned, a regularized solution is obtained by employing truncated singular value decomposition, while the regularization parameter for the regularization method is provided by the L‐curve method. Numerical results are presented for both smooth and piecewise smooth geometry. The stability of the method with respect to the noise in the data is investigated by using simulated noisy data. The results show that the method is highly accurate, computationally efficient and stable, and can be a competitive alternative to existing methods for the numerical solution of the problems. Copyright © 2005 John Wiley & Sons, Ltd.     

上半平面上的第一类椭圆型方程组的自然积分方程

主要用Fourier变换法讨论第一类椭圆型方程组在上半平面的自然积分方程和Poisson积分公式。     

An Artificial Boundary Condition for a Class of Quasi-Newtonian Stokes Flows

An artificial boundary condition method, derived in terms of infinite Fourier series, is applied to solve a class of quasi-Newtonian Stokes flows. Based on the natural boundary reduction involving an artificial condition on the artificial boundary, the coupled variational problem and its numerical solution are obtained. The unique solvability of the continuous and discrete formulations are discussed, and the error analysis for the problem is also considered. Finally, an a posteriori error estimate for the corresponding problem is provided.     

体积源边界点法及其在声辐射计算中的应用

利用作者提出的体积源边界点法，对结构振动声辐射的计算问题进行了研究。给出了该法在全频率范围内存在唯一解的数学证明；通过诸多不同边界曲面和不同边界量分布的声辐射算例，从计算精度、计算稳定性、对振动体任意表面几何形状的适应性以及克服解的非唯一性等方面，对该方法的有效性进行了检验。     

Dispersion analysis and error estimation of Galerkin finite element methods for the Helmholtz equation

When applying numerical methods for the computation of stationary waves from the Helmholtz equation, one obtains ‘numerical waves’ that are dispersive also in non-dispersive media. The numerical wave displays a phase velocity that depends on the parameter k of the Helmholtz equation. In dispersion analysis, the phase difference between the exact and the numerical solutions is investigated. In this paper, the authors' recent result on the phase difference for one-dimensional problems is numerically evaluated and discussed in the context of other work directed to this topic. It is then shown that previous error estimates in H¹-norm are of nondispersive character but hold for medium or high wavenumber on extremely refined mesh only. On the other hand, recently proven error estimates for constant resolution contain a pollution term. With certain assumptions on the exact solution, this term is of the order of the phase difference. Thus a link is established between the results of dispersion analysis and the results of numerical analysis. Throughout the paper, the presentation and discussion of theoretical results is accompanied by numerical evaluation of several model problems. Special attention is given to the performance of the Galerkin method with a higher order of polynomial approximation p(h-p-version).