三维问题边界元法中几乎奇异积分的正则化算法

当源点靠近边界单元时,边界积分方程通常存在几乎奇异积分的计算难题.基于三角形单元,将源点到单元的距离与单元特征长度比值定义为接近度,用于度量边界单元中积分奇异性的程度.将单元上的面积分在局部的极坐标系ρθ下表示,利用一些初等函数的积分公式,获得对变量ρ作单层积分的解析表达式.几乎强奇异和超奇异面积分被转化为沿单元围道上一系列线积分,而Gauss数值积分能够有效计算这些线积分.应用该算法分析三维弹性薄壁结构获得了成功.     

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

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

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

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

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

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

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

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

三维边界元法中高阶单元上的几乎奇异积分

三维边界元分析中,高阶几何单元上的几乎奇异积分计算是一个重要而且困难的问题,该文对此进行了研究。使用8节点四边形和6节点三角形曲面单元来描述几何边界;构造了新的距离函数;拓展原有的指数函数非线性变换到三维边界元法中,利用拓展的变换来消除被积函数的几乎奇异性。数值算例表明,该算法稳定,效率高,即使计算点到实际边界的距离很小,依然可获得令人满意的数值解。     

关于边界元法中奇异积分的处理

关于边界元法中奇异积分的处理臧跃龙，嵇醒（西安交通大学，７１００４９）（上海同济大学，上海２０００９２）关键词边界元法，积分奇异性，奇异性的消除１引言与有限元法不同，边界元法的数值积分通常带有对数或一、二阶奇异性．如二维问题，基本解带有对数奇异性，其...     

一种处理弹性动力边界元法中奇异积分的方法

利用边界元法求解瞬态弹性动力学问题时,时域基本解函数的分段连续性和奇异性为该问题的求解带来很大的困难。为了解决时域基本解中的奇异性问题,本文依据柯西主值的定义,对经过时间解析积分之后的时域基本解进行奇异值分解,将其分成奇异和正则积分两部分;其中正则部分可通过采用常规高斯积分方法来计算,而奇异部分具有简单的形式,可以利用解析积分计算。经过上述操作之后,就可以达到直接消除时域基本解中奇异积分的目的。和传统方法相比,本文方法并不依赖静力学基本解来消除奇异性,是一种直接求解方法。最后给定两个数值算例来验证本文提出方法的正确性和可行性,结果表明使用本文算法可以解决弹性动力学边界积分方程中的奇异性问题。     

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

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

边界元法计算近边界点参量的一个通用算法

针对边界元法存在近力界点参量计算的困难,给出了一个通用性方法,将近边界点到边界单元的距离参数通过分部积分变换到积分式之外,从而计算出二维问题近边界点参量的几乎强奇异和超奇异积分,由此,对任何近边界点参量,提出了整套计算方案,算例证明了本法的有效性。     

导数场边界积分方程分析近边界应力分布

研究二维弹性力学问题边界积分方程,通过分部积分变换消除了常规导数边界积分方程中的超奇异积分,获得仅含强奇异积分的应力自然边界积分方程.对于近边界应力的计算,进一步运用正则化算法解析计算其中的几乎强奇异积分.较常规边界元法相比,应力自然边界积分方程可以求解离边界更加接近的内点应力值.算例证明了文中方法的可应用性和有效性.     

边界元法分析狭长体结构

针对边界元法分析狭长结构时遇到的几乎奇异积分难以计算的困难,将几乎奇异积分划分为两种类型,分别通过分部积分交换把引起积分几乎奇异的参量移至积分号之外,从而建立了一个新的正则化算法,解决了边界积分方程中几乎奇异积分的计算难题。文中用边界元法计算了弹性力学平面问题的狭长结构,算例证明了本法的有效性。     

二维弹性问题边界元法中边界层效应问题的变换法

基于间接规则化边界积分方程,有效估计奇异边界积分,准确求得边界量,为场变量的计算奠定了基础。在计算场变量时,针对二维弹性力学边界元法中出现的几乎奇异积分,本文采用一类非线性变量替换法,有效地改善了被积函数的震荡特性,从而消除了核积分的几乎奇异性;在不增加计算量的情况下,极大地改进了几乎奇异积分计算的精度,成功地求解了弹性体近边界点上的力学参量,避免了边界层效应。此外,本文引入一种精确几何单元逼近,对于圆弧边界,这样的插值逼近几乎是精确的,提高了计算精度。数值算例表明,本文算法稳定,效率高,并可达到很高的计算精度,即使场点非常靠近边界,如场点到积分单元的距离小到纳米级,仍可避免边界层效应现象。