A quadrature formula for integrals with nearby singularities

The purpose of this paper is to propose a new quadrature formula for integrals with nearby singularities. In the boundary element method, the integrands of nearby singular boundary integrals vary drastically with the distance between the field and the source point. Especially, field variables and their derivatives at a field point near a boundary cannot be computed accurately. In the present paper a quadrature formula for ??‐isolated singularities near the integration interval, based on Lagrange interpolatory polynomials, is obtained. The error estimation of the proposed formula is also given. Quadrature formulas for regular and singular integrals with conjugate poles are derived. Numerical examples are given and the proposed quadrature rules present the expected polynomial accuracy. Copyright © 2006 John Wiley & Sons, Ltd.     

Evaluations of hypersingular integrals using Gaussian quadrature

A Gaussian quadrature formula for hypersingular integrals with second‐order singularities is developed based on previous Gaussian quadrature formulae for Cauchy principal value integrals. The formula uses classical orthonormal polynomials, and the formula is then specialized to the case of Legendre and Chebyshev polynomials. Numerical experiments are carried out using the current formula and a previous formula developed by Kutt. It is found that the two methods generally give similar results, and in some cases the current method works better. It has also been shown that the current method allows the choice of an appropriate weight which can increase the convergence rate and the accuracy of the results. Copyright © 1999 John Wiley & Sons, Ltd.     

Direct Gauss quadrature formulae for logarithmic singularities on isoparametric elements

Weakly singular logarithmic integrals occur in two-dimensional (2D) BEM problems when the integration is over an element which includes the current source node. For boundary elements with curved geometry such as the quadratic and cubic element then numerical integration must be used. Recently developed direct Gauss quadrature schemes which implicitly consider the integral to include a sum of singular and non-singular terms appear to be superior to conventional schemes which represent these separate terms explicitly. This paper discusses these direct quadratures and introduces new Gauss schemes which integrate exactly the logarithmic singularities on any one of the 3 nodes of a quadratic element using a single formula. This new quadrature may also be used for linear and constant elements. Similar quadrature formulae for the 4 singularities on a cubic element are also discussed. This new approach is both accurate and simple, reducing the size of the computer program and allowing the use of the same quadrature for several isoparametric types.     

A modified Gauss quadrature formula with special integration points for evaluation of Quasi-singular integrals

It is well known in the boundary element method that integration rules fail when the integrand presents a nearby singularity. This drawback arises when the field point is near the source point, i.e. in the case of a domain with very narrow boundaries or when the field point where we try to calculate stresses or any other field variables, is near the boundaries.In the present paper a quadrature formulas for ℓ isolated singularities near the integration interval, based on ordinary or special Langrange interpolatory polynomials, is obtained. This interpolatory formulas present similarities with known formulas for the numerical evaluation of singular integrals. Quadrature formulas for regular and singular integrals with conjugate poles are also derived. Numerical examples are given and the proposed quadrature rules present the expected polynomial accuracy.     

The evaluation of improper integrals encountered in the use of conformal transformation

Methods are described and reviewed for the accurate numerical evaluation of improper integrals encountered in conformal transformation solutions involving boundaries of relatively complicated shape. Four methods are reviewed for the solution of integrals containing end-point singularities. Two new methods are discussed for the solution of integrals containing both end point singularities and simple poles within the range of integration. One method uses a combination of a simple recursive formula and the coefficients of a Chebyshev series and a second method involves subtracting out the singularity and the use of Gauss–Jacobi quadrature. Both methods can give results of high accuracy and an upper limit of the error is readily found. A numerical example is taken which is typical of the application to practical problems and this brings out a comparison of the two methods.     

Complex hypersingular integrals and integral equations in plane elasticity

Summary Complex hypersingular (finite-part) integrals and integral equations are considered in the functional class of N. Muskhelishvili. The appropriate definition is given. Three regularization (equivalence) formulae follow from this definition. They reduce hypersingular integrals to singular ones and allow to derive hypersingular analogues for Sokhotsky-Plemelj's formulae and for conditions that are necessary and sufficient for the function to be piecewise holomorphic. Two approaches to get and investigate complex hypersingular equations follow from these results: one of them is based on the equivalence formulae; as to the other, it is based on above-mentioned conditions. As an example, authors' equation for plane elasticity is studied. The existence of a unique solution is stated and some advantages over singular equations are outlined. To solve hypersingular equations the quadrature rules are presented. The accuracy of different quadrature formulae is compared, the examples being used. They confirm the need to take into account asymptotics and to carry out a thorough analytical investigation to get safe numerical results.     

On non‐linear transformations for accurate numerical evaluation of weakly singular boundary integrals

This paper presents a study of the performance of the non‐linear co‐ordinate transformations in the numerical integration of weakly singular boundary integrals. A comparison of the smoothing property, numerical convergence and accuracy of the available non‐linear polynomial transformations is presented for two‐dimensional problems. Effectiveness of generalized transformations valid for any type and location of singularity has been investigated. It is found that weakly singular integrals are more efficiently handled with transformations valid for end‐point singularities by partitioning the element at the singular point. Further, transformations which are excellent for CPV integrals are not as accurate for weakly singular integrals. Connection between the maximum permissible order of polynomial transformations and precision of computations has also been investigated; cubic transformation is seen to be the optimum choice for single precision, and quartic or quintic one, for double precision computations. A new approach which combines the method of singularity subtraction with non‐linear transformation has been proposed. This composite approach is found to be more accurate, efficient and robust than the singularity subtraction method and the non‐linear transformation methods. Copyright © 2001 John Wiley & Sons, Ltd.     

Application of sigmoidal transformations to weakly singular and near‐singular boundary element integrals

Accurate numerical determination of line integrals is fundamental to reliable implementation of the boundary element method. For a source point distant from a particular element, standard Gaussian quadrature is adequate, as well as being the technique of choice. However, when the integrals are weakly singular or nearly singular (source point near the element) this technique is no longer adequate. Here a co‐ordinate transformation technique, based on sigmoidal transformations, is introduced to evaluate weakly singular and near‐singular integrals. A sigmoidal transformation has the effect of clustering the integration points towards the endpoints of the interval of integration. The degree of clustering is governed by the order of the transformation. Comparison of this new method with existing co‐ordinate transformation techniques shows that more accurate evaluation of these integrals can be obtained. Based on observations of several integrals considered, guidelines are suggested for the order of the sigmoidal transformations. Copyright © 1999 John Wiley & Sons, Ltd.     

An efficient numerical method for the solution of sliding contact problems

In this paper, an efficient numerical method to solve sliding contact problems is proposed. Explicit formulae for the Gauss–Jacobi numerical integration scheme appropriate for the singular integral equations of the second kind with Cauchy kernels are derived. The resulting quadrature formulae for the integrals are valid at nodal points determined from the zeroes of a Jacobi polynomial. Gaussian quadratures obtained in this manner involve fixed nodal points and are exact for polynomials of degree 2n ? 1, where n is the number of nodes. From this Gauss–Jacobi quadrature, the existing Gauss–Chebyshev quadrature formulas can be easily derived. Another apparent advantage of this method is its ability to capture correctly the singular or regular behaviour of the tractions at the edge of the region of contact. Also, this analysis shows that once if the total normal load and the friction coefficient are given, the external moment M and contact eccentricity e (for incomplete contact) in fully sliding contact are uniquely determined. Finally, numerical solutions are computed for two typical contact cases, including sliding Hertzian contact and a sliding contact between a flat punch with rounded corners pressed against the flat surface of a semi‐infinite elastic solid. These results provide a demonstration of the validity of the proposed method. Copyright © 2005 John Wiley & Sons, Ltd.     

A mapping method for numerical evaluation of two-dimensional integrals with 1/r singularity

Singular integrals occur commonly in applications of the boundary element method (BEM). A simple mapping method is presented here for the numerical evaluation of two-dimensional integrals in which the integrands, at worst, are O(1/r) (r being the distance from a source to a field point). This mapping transforms such integrals over general curved triangles into regular 2-D integrals. Over flat and curved quadratic triangles, regular line integrals are obtained, and these can be easily evaluated by standard Gaussian quadrature. Numerical tests on some typical singular integrals, encountered in BEM applications, demonstrate the accuracy and efficacy of the method.     

Semi‐sigmoidal transformations for evaluating weakly singular boundary element integrals

Accurate numerical integration of line integrals is of fundamental importance to reliable implementation of the boundary element method. Usually, the regular integrals arising from a boundary element method implementation are evaluated using standard Gaussian quadrature. However, the singular integrals which arise are often evaluated in another way, sometimes using a different integration method with different nodes and weights. Here, a co‐ordinate transformation technique is introduced for evaluating weakly singular integrals which, after some initial manipulation of the integral, uses the same integration nodes and weights as those of the regular integrals. The transformation technique is based on newly defined semi‐sigmoidal transformations, which cluster integration nodes only near the singular point. The semi‐sigmoidal transformations are defined in terms of existing sigmoidal transformations and have the benefit of evaluating integrals more accurately than full sigmoidal transformations as the clustering is restricted to one end point of the interval. Comparison of this new method with existing coordinate transformation techniques shows that more accurate evaluation of weakly singular integrals can be obtained. Based on observation of several integrals considered, guidelines are suggested for the type of semi‐sigmoidal transformation to use and the degree to which nodes should be clustered at the singular points. Copyright © 2000 John Wiley & Sons, Ltd.     

A direct approach for boundary integral equations with high‐order singularities

Boundary integral equations with extremely singular (i.e., more than hypersingular) kernels would be useful in several fields of applied mechanics, particularly when second‐ and third‐order derivatives of the primary variable are required. However, their definition and numerical treatment pose several problems. In this paper, it is shown how to obtain these boundary integral equations with still unnamed singularities and, moreover, how to efficiently and reliably compute all the singular integrals. This is done by extending in full generality the so‐called direct approach. Only for definiteness, the method is presented for the analysis of the deflection of thin elastic plates. Numerical results concerning integrals with singularities up to order r⁻⁴ are presented to validate the proposed algorithm. Copyright © 2000 John Wiley & Sons, Ltd.     

Singular integral equations of the second kind with generalized Cauchy‐type kernels and variable coefficients

A numerical solution method is presented for singular integral equations of the second kind with a generalized Cauchy kernel and variable coefficients. The solution is constructed in the form of a product of regular and weight functions. The weight function possesses complex singularities at the ends of the interval. The parameters defining the power of these singularities are obtained by solving for the characteristic equations. A Gauss–Chebychev quadrature formula is utilized in the numerical solution of the integral equations. Benchmark examples are considered in order to illustrate the validity of the solution method. Copyright © 1999 John Wiley & Sons, Ltd.     

Direct computation of Cauchy principal value integrals in advanced boundary elements

This paper presents a new method for the direct computation of strongly singular integrals existing in the Cauchy principal value sense. It can be usefully applied in the solution by the direct boundary element method of many different problems. Initially, some considerations are provided in order to summarize the state-of-the-art on this issue. Then the features of the proposed method are reported. The procedure allows the direct calculation of Cauchy principal value integrals with first-order singularity and it is applicable even in advanced boundary element methods employing high-order elements. It requires only the use of standard Gaussian quadrature formulae plus the computation of a logarithmic term. Some examples show the effectiveness and efficiency of the procedure.     

Analysis and evaluation of singular integrals by the invariant imbedding approach

The invariant imbedding formulae for evaluating singular integrals on flat domains are used to analyse the singular integrals which arise in the Boundary Element Method. Although all orders of singularity are considered, the primary objective of the paper is the evaluation of Cauchy principal value integrals (CPV). Conditions for boundedness of the CPV aré derived, and evaluation techniques are presented for situations in which the singular point lies on, or is in close proximity to, an element boundary. Specific cases analysed include straight and curved element boundaries, and corners. Other numerical aspects of the formulae are also discussed. By way of example, the analysis techniques are applied to the volume integrals which arise in continuum, inelastic problems.     

The general type of finite-part singular integrals and integral equations with logarithmic singularities used in fracture mechanics

Summary A new method is proposed, by using some special quadrature rules, for the numerical evaluation of the general type of finite-part singular integrals and integral equations with logarithmic singularities. In this way the system of such equations can be numerically solved by reduction to a system of linear equations. For this reduction, the singular integral equation is applied to a number of appropriately selected collocation points on the integration interval, and then a numerical integration rule is used for the approximation of the integrals in this equation. An application is given, to the determination of the intensity of the logarithmic singularity in a simple crack inside an infinite, isotropic solid.With 1 Figure     

运动介质中奇异边界元积分式的精确求解

采用边界元方法求解与运动介质相关声学问题时,难点之一是如何精确计算场点与源点重合所导致的奇异积分式。论文提出一种将具有奇性的单元面积分式拆分为奇性和非奇性积分部分分别进行计算的新方法。对奇性积分部分,经过严格的数学推导给出解析解;而对非奇性积分部分则通过高斯积分法处理。新方法可有效地提高边界元计算精度和效率,对运动介质中的有关声学问题的边界元数值计算具有重要意义。     

A sinh transformation for evaluating nearly singular boundary element integrals

An implementation of the boundary element method requires the accurate evaluation of many integrals. When the source point is far from the boundary element under consideration, a straightforward application of Gaussian quadrature suffices to evaluate such integrals. When the source point is on the element, the integrand becomes singular and accurate evaluation can be obtained using the same Gaussian points transformed under a polynomial transformation which has zero Jacobian at the singular point. A class of integrals which lies between these two extremes is that of ‘nearly singular’ integrals. Here, the source point is close to, but not on, the element and the integrand remains finite at all points. However, instead of remaining flat, the integrand develops a sharp peak as the source point moves closer to the element, thus rendering accurate evaluation of the integral difficult. This paper presents a transformation, based on the sinh function, which automatically takes into account the position of the projection of the source point onto the element, which we call the ‘nearly singular point’, and the distance from the source point to the element. The transformation again clusters the points towards the nearly singular point, but does not have a zero Jacobian. Implementation of the transformation is straightforward and could easily be included in existing boundary element method software. It is shown that, for the two‐dimensional boundary element method, several orders of magnitude improvement in relative error can be obtained using this transformation compared to a conventional implementation of Gaussian quadrature. Asymptotic estimates for the truncation errors are also quoted. Copyright © 2004 John Wiley & Sons, Ltd.     

Error estimation of quadrature rules for evaluating singular integrals in boundary element problems

The efficient numerical evaluation of integrals arising in the boundary element method is of considerable practical importance. The superiority of the use of sigmoidal and semi‐sigmoidal transformations together with Gauss–Legendre quadrature in this context has already been well‐demonstrated numerically by one of the authors. In this paper, the authors obtain asymptotic estimates of the truncation errors for these algorithms. These estimates allow an informed choice of both the transformation and the quadrature error in the evaluation of boundary element integrals with algebraic or algebraic/logarithmic singularities at an end‐point of the interval of integration. Copyright © 2000 John Wiley & Sons, Ltd.     

Direct evaluation of singular integrals in elastoplastic analysis by the boundary element method

This paper presents a new method for the direct and accurate evaluation of strongly singular integrals in the sense of Cauchy principal values and weakly singular integrals over quadratic boundary elements in three-dimensional stress analysis and quadratic internal cells in two-dimensional elastoplastic analysis by the boundary element method. A quadratic triangle polar co-ordinate transformation technique is applied to reduce the order of singularity of the singular integrals. Next, a form of Stokes' theorem is introduced in order to remove the singularity in the Cauchy principal value integrals; therefore, the evaluation of these integrals can be carried out by standard Gaussian quadrature. Numerical examples of 2-D elastoplastic problems and a 3-D elastic problem show the effectiveness and efficiency of the method.