首页|期刊导航|电工技术学报|曲面边界元中基于投影点单元细分的近奇异积分计算

曲面边界元中基于投影点单元细分的近奇异积分计算OA

Calculation of Nearly Singular Integrals in Curved Boundary Element Method Based on Element Subdivision at Projection Points

中文摘要英文摘要

边界元法中近奇异积分的处理一直是影响其计算精度的关键因素,尤其在薄型结构分析和复杂几何问题中表现尤为突出.近奇异积分由于其数学特性和数值计算的复杂性,处理难度往往高于奇异积分.在电场计算的间接曲面边界元法中,该文提出一种通用的近奇异积分数值计算方法.该方法通过将场点在积分单元上的投影点作为细分依据,在参数域内对单元进行自适应细分.该方法适用于任意场点位置下的平面单元和曲面单元.在多种单元类型和整体算例模型中计算验证,结果表明该方法相较于传统的高斯积分法显著提高了计算精度.通过对近奇异积分的精确处理,该方法有效提升了边界元法在实际工程问题分析中的计算精度,有利于边界元法在更广泛领域的应用.

The boundary element method(BEM)is a numerical computation method developed from the ideas of element division and discretization of the finite element method(FEM).The BEM is challenging to handle alongside other numerical methods,so it is important to leverage its strengths in high precision and computational efficiency.The singular/nearly singular integration problem,unique to BEM,has significantly affected its accuracy.In particular,the nearly singular integrals are often more difficult to handle than singular integrals due to their mathematical properties and the complexity of numerical computation. This paper proposes ageneralized numerical computation method for nearly singular integrals to improve the accuracy of BEM,based on the curved boundary element method(CBEM)and coordinate transformation.The method subdivides the element in a scaled form,using the field-point projections onto the integration element as the basis for subdivision.For slender elements that often appear in surface meshing,the scaling is optimized by considering the longer side's dimensions to balance sub-element dimensions in all directions,thereby improving computational accuracy.Since the subdivision process is carried out in the parameter domain,the method can be applied to all planar and curved elements.In addition,adaptive methods are used to determine the number of subdivisions,ensuring high accuracy at any field point location. In the single-element calculation,the planar triangular element and the spherical quadrilateral element are considered.The results show that,compared with the ordinary 4-point Gaussian integration and Gaussian point subdivision methods,the proposed method achieves higher accuracy and better adaptability to different field point locations,especially on curved elements.Through iterative convergence,the new method can effectively ensure a sufficient number of subdivisions.A model of an indoor ring conductor is constructed.In the comparative analysis,the proposed method effectively improves overall field accuracy and the maximum value.Subsequently,this paper takes the electric-field analysis problem of the actual converter valve tower as an example.The calculations are accelerated using the fast multipole algorithm,and the results are obtained on anordinary PC.The results show that the method achieves higher accuracy when using a similar number of Gaussian points,especially for the maximum field strength,which is of greater concern in engineering.In addition,the advantages of BEM in terms of computational time and cost are demonstrated in this calculation example. The element subdivision method based on projection points effectively improves the accuracy of the nearly singular integrals,thereby reducing the overall computational error of BEM.This improvement highlights the advantage of BEM,demonstrating great potential across a wide range of applications.Therefore,this method is expected to provide a more accurate and efficient solution for the computation of complex engineering problems.

段一伟;王泽忠

华北电力大学电气与电子工程学院 北京 102206华北电力大学电气与电子工程学院 北京 102206

信息技术与安全科学

边界元法近奇异积分单元细分曲面单元电场计算

Boundary element methodnearly singular integralselement subdivisioncurved elementelectric field calculation

《电工技术学报》 2026 (10)

3221-3229,9

10.19595/j.cnki.1000-6753.tces.250894

评论