正交各向异性位势问题边界元法中几乎奇异积分的解析算法

几乎奇异积分的计算困难阻碍了边界元法的工程应用。本文针对二维正交各向异性位势问题边界元法中近边界点的几乎奇异积分,采用分部积分法,导出一种直接的解析计算公式。该解析公式可以精确计算线性单元上的几乎奇异积分。对二次单元,可将其细分为几个线性元,采用该解析公式近似计算其边界积分。当内点离当前积分单元较远时,仍保持常规高斯数值积分模式;而当内点离其较近时,因常规高斯积分结果失效,则采用该解析积分取代高斯数值积分。数值算例证明了该算法的有效性和精确性。二次元计算结果比线性元计算结果更精确。     

薄体位势问题边界元法中的解析积分算法

薄体结构的数值分析是边界元法的难点问题之一。该文导出了一种完全解析积分算法，用这种算法计算了薄体平面位势问题边界元法中出现的几乎弱奇异、强奇异和超奇异积分。当边界离散为一系列线性单元，边界积分方程离散计算的积分可归纳为三种形式。对薄体问题，源点与积分单元距离通常相距很近，这些积分产生显著几乎奇异性，直接采用常规高斯积分不能有效计算。为此该文导出了这些几乎奇异积分的全解析计算公式。按源点与单元的距离是否为零，公式分两种情况。新算法采用全解析积分公式处理几乎奇异积分，首先精确计算出薄体问题边界未知位势和法向位势梯度，然后再进一步计算了域内点的物理参量。算例表明该文算法可处理狭长比为1．E-08的薄体问题，显示了边界元法分析薄体问题具有独特的优势。     

三维变系数热传导问题边界元分析中几乎奇异积分计算

在边界积分的数值计算过程中,当源点离积分单元很近时,边界积分就会具有几乎奇异性,此时不能直接用高斯数值积分公式计算几乎奇异积分。本文以三维非均质热传导问题为例,介绍了一种计算几乎奇异边界积分的新方法。首先,采用Newton-Raphson迭代算法确定积分单元上离源点最近的点;然后,将积分单元上任意一点的坐标在最近点处展开成泰勒级数,并计算源点到积分单元任意点的距离;最后,将距离函数代入几乎奇异边界积分中,并运用指数变换方法导出积分单元上几乎奇异积分的计算公式。文中给出了两个非均质热传导问题的算例来验证所述方法的正确性、有效性和稳定性。     

二维热弹性力学边界元法中几乎奇异积分的正则化

针对二维热弹性力学边界元法中近边界点的几乎强奇异和超奇异积分,采用一种通用算法,将其实施正则化．该方法适用于线性单元,与近边界点邻近的单元上的积分采用正则化积分公式计算,远处单元的积分仍保持常规高斯积分．算例证明了该法的有效性和精确性．     

三维边界元法高阶元几乎奇异积分半解析法

分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.      

A NEW SEMI-ANALYTIC ALGORITHM OF NEARLY SINGULAR INTEGRALS IN HIGH ORDER BOUNDARY ELEMENT ANALYSIS OF 2D ELASTICITY

针对边界元法中高阶单元中几乎奇异积分计算难题,解剖了二维边界元法高阶单元的几何特征,定义源点相对高阶单元的接近度。将高阶单元上奇异积分核函数用近似奇异函数逼近,从而分离出积分核中主导的奇异函数部分,其奇异积分核分解为规则核函 数和奇异核函数两项积分之和。规则核函数用常规高斯数值积分,再对奇异核函数积分导出解析公式,从而建立了一种新的半解析法,用于高阶边界单元上几乎强奇异和超奇异积分计算。给出3个算例,采用边界元法高阶单元的半解析法计算了弹性力学薄体结构和近边界点位移/应力,并与线性边界元正则化算法结果作了比较,结果表明提出的二次元的半解析算法更加有效。特别是分析薄体结构,采用正则化算法的线性边界元分析比有限元有显著优势,而用提出的二次边界元半解析算法分析比其线性元的有效接近度又减小了4个量级。     

二维位势边界元法高阶单元几乎奇异积分半解析算法

准确计算几乎奇异积分是边界元法难题之一。目前,对于一般的高阶单元的几乎奇异积分尚缺乏通用高效的计算方法。本文在单元局部坐标系中表征了二维高阶单元的几何特征,提出了源点相对高阶单元的接近度概念。针对二维位势边界元法的3节点二次等参单元,构造出与单元积分核具有相同几乎奇异性的近似奇异核函数。从二维位势几乎奇异积分单元积分核中扣除近似奇异核函数,把几乎奇异积分项转换为规则积分和奇异积分两部分之和,规则积分部分用常规Gauss数值积分计算,奇异积分部分由导出的解析公式计算,从而建立了二维位势问题高阶单元几乎强奇异和超奇异积分的半解析算法。算例结果表明了本文半解析算法的有效性和计算精度。     

三维Laplace方程边界元中线性单元的精确积分法

边界元中的边界积分计算影响计算精度和计算速度。非奇异积分一般采用数值积分，当配置点接近积分单元时，计算精度降低。未知函数线性插值得到的解是连续解，但计算难度增大。本文采用积分区域变换，将三维Laplace问题的二维积分化为一维积分，这样奇异积分和非奇异积分能采用精确积分的方法计算，使求解精度，计算速度都得到提高。     

A NEW SEMI-ANALYTIC ALGORITHM OF NEARLY SINGULAR INTEGRALS IN HIGH ORDER BOUNDARY ELEMENT ANALYSIS OF 3D POTENTIAL

分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.     

二维正交异性位势问题高阶边界元几乎奇异积分半解析算法

边界元法中高阶单元上的几乎奇异积分一直难以计算。针对正交各向异性位势问题,提出一个半解析算法准确计算了其高阶单元上的几乎奇异积分。首先将正交各向异性材料中源点到单元的距离函数在局部坐标系下渐近展开,采用级数展开式构造出与奇异积分核函数具有相同奇异性的可积近似核函数;然后利用扣除法的思想,原奇异积分核减去近似积分核后再加回,几乎奇异积分便转换为规则部分和奇异部分之和,规则积分采用Gauss数值积分计算,奇异积分由文中推导出解析公式计算。通过两个正交各向异性的热传导算例表明,本文建立的高阶单元半解析算法能准确高效地计算近边界内点位势和位势梯度。     

二维正交各向异性位势问题的高阶单元快速多极边界元法

研制了一种适用于二维正交各向异性位势问题的高阶单元(线性单元和二次单元)快速多极边界元法. 在快速多极边界元法中, 源点对于远场区域的积分采用快速多极展开式计算, 而对于近场区域的积分则直接进行计算. 高阶单元的使用使得近场积分, 尤其是奇异积分和几乎奇异积分的计算更加复杂. 通过引入复数表达对其进行简化, 若边界采用线性单元插值, 近场积分可直接解析计算; 若采用二次单元插值, 则给出一个半解析算法计算近场积分. 高阶单元奇异积分和几乎奇异积分计算难题的解决, 使得高阶单元快速多极边界元法不仅能够计算一般结构, 也能被应用于超薄体结构, 拓宽了高阶单元快速多极边界元法的适用范围. 数值算例表明, 若计算精度一定, 高阶单元快速多极边界元法较常值单元快速多极边界元法使用的单元数量显著减少, 且高阶单元快速多极边界元法计算时间与自由度数量成线性关系, 其计算效率仍处于$O(N)$量级, 因此高阶单元快速多极边界元法可更加高效求解大规模问题.      

REGULARIZATION OF NEARLY SINGULAR INTEGRALS IN THE BOUNDARY ELEMENT METHOD OF POTENTIAL PROBLEMS

IntroductionAsanimportantnumericalmethod ,BoundaryElementMethod (BEM)hasbeenappliedinmanyareas[1].However,theBEMhasthedifficultiesofcalculatingsingularintegralsatnodesonboundaryoratinteriorpointsveryclosetotheboundary .TheaccuracyoftheBEMdependsontheprecisionofthecalculatedvaluesofthesingularintegrals,toagreatdegree.Manyresearchersdevotethemselvestothetreatmentofthesingularintegrals[2～3],whicharereviewedindetailbyRef.[4] .Ageneralregularizationalgorithmofevaluatingthephysicalquantitiesa…     

边界积分方程中近奇异积分计算的一种变量替换法

准确估计近奇异边界积分是边界元分析中一项很重要的课题,其重要性仅次于对奇异积分的处理. 近年来已发展了许多方法,都取得了一定程度的成功,但这个问题至今仍未得到彻 底的解决. 基于一种新的变量变换的思想和观点,提交了一种通用的积分变换法, 它非常有效地改善了被积函数的震荡特性,从而消除了积分的近奇异性,在不增加计算量的情况 下, 极大地改进了近奇异积分计算的精度. 数值算例表明,其算法稳定,效率高, 并可达到很高的计算精度,即使区域内点非常地靠近边界,仍可取得很理想的结果.