共查询到14条相似文献,搜索用时 218 毫秒
1.
分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效. 相似文献
2.
二维位势边界元法高阶单元几乎奇异积分半解析算法 总被引:1,自引:1,他引:0
准确计算几乎奇异积分是边界元法难题之一。目前,对于一般的高阶单元的几乎奇异积分尚缺乏通用高效的计算方法。本文在单元局部坐标系中表征了二维高阶单元的几何特征,提出了源点相对高阶单元的接近度概念。针对二维位势边界元法的3节点二次等参单元,构造出与单元积分核具有相同几乎奇异性的近似奇异核函数。从二维位势几乎奇异积分单元积分核中扣除近似奇异核函数,把几乎奇异积分项转换为规则积分和奇异积分两部分之和,规则积分部分用常规Gauss数值积分计算,奇异积分部分由导出的解析公式计算,从而建立了二维位势问题高阶单元几乎强奇异和超奇异积分的半解析算法。算例结果表明了本文半解析算法的有效性和计算精度。 相似文献
3.
4.
三维变系数热传导问题边界元分析中几乎奇异积分计算 总被引:2,自引:2,他引:0
在边界积分的数值计算过程中,当源点离积分单元很近时,边界积分就会具有几乎奇异性,此时不能直接用高斯数值积分公式计算几乎奇异积分。本文以三维非均质热传导问题为例,介绍了一种计算几乎奇异边界积分的新方法。首先,采用Newton-Raphson迭代算法确定积分单元上离源点最近的点;然后,将积分单元上任意一点的坐标在最近点处展开成泰勒级数,并计算源点到积分单元任意点的距离;最后,将距离函数代入几乎奇异边界积分中,并运用指数变换方法导出积分单元上几乎奇异积分的计算公式。文中给出了两个非均质热传导问题的算例来验证所述方法的正确性、有效性和稳定性。 相似文献
5.
正交各向异性位势问题边界元法中几乎奇异积分的解析算法 总被引:3,自引:0,他引:3
几乎奇异积分的计算困难阻碍了边界元法的工程应用。本文针对二维正交各向异性位势问题边界元法中近边界点的几乎奇异积分,采用分部积分法,导出一种直接的解析计算公式。该解析公式可以精确计算线性单元上的几乎奇异积分。对二次单元,可将其细分为几个线性元,采用该解析公式近似计算其边界积分。当内点离当前积分单元较远时,仍保持常规高斯数值积分模式;而当内点离其较近时,因常规高斯积分结果失效,则采用该解析积分取代高斯数值积分。数值算例证明了该算法的有效性和精确性。二次元计算结果比线性元计算结果更精确。 相似文献
6.
7.
边界元法分析狭长体结构 总被引:3,自引:0,他引:3
针对边界元法分析狭长结构时遇到的几乎奇异积分难以计算的困难,将几乎奇异积分划分为两种类型,分别通过分部积分交换把引起积分几乎奇异的参量移至积分号之外,从而建立了一个新的正则化算法,解决了边界积分方程中几乎奇异积分的计算难题。文中用边界元法计算了弹性力学平面问题的狭长结构,算例证明了本法的有效性。 相似文献
8.
9.
薄体位势问题边界元法中的解析积分算法 总被引:1,自引:0,他引:1
薄体结构的数值分析是边界元法的难点问题之一。该文导出了一种完全解析积分算法,用这种算法计算了薄体平面位势问题边界元法中出现的几乎弱奇异、强奇异和超奇异积分。当边界离散为一系列线性单元,边界积分方程离散计算的积分可归纳为三种形式。对薄体问题,源点与积分单元距离通常相距很近,这些积分产生显著几乎奇异性,直接采用常规高斯积分不能有效计算。为此该文导出了这些几乎奇异积分的全解析计算公式。按源点与单元的距离是否为零,公式分两种情况。新算法采用全解析积分公式处理几乎奇异积分,首先精确计算出薄体问题边界未知位势和法向位势梯度,然后再进一步计算了域内点的物理参量。算例表明该文算法可处理狭长比为1.E-08的薄体问题,显示了边界元法分析薄体问题具有独特的优势。 相似文献
10.
针对边界元法中高阶单元中几乎奇异积分计算难题,解剖了二维边界元法高阶单元的几何特征,定义源点相对高阶单元的接近度。将高阶单元上奇异积分核函数用近似奇异函数逼近,从而分离出积分核中主导的奇异函数部分,其奇异积分核分解为规则核函 数和奇异核函数两项积分之和。规则核函数用常规高斯数值积分,再对奇异核函数积分导出解析公式,从而建立了一种新的半解析法,用于高阶边界单元上几乎强奇异和超奇异积分计算。给出3个算例,采用边界元法高阶单元的半解析法计算了弹性力学薄体结构和近边界点位移/应力,并与线性边界元正则化算法结果作了比较,结果表明提出的二次元的半解析算法更加有效。特别是分析薄体结构,采用正则化算法的线性边界元分析比有限元有显著优势,而用提出的二次边界元半解析算法分析比其线性元的有效接近度又减小了4个量级。 相似文献
11.
12.
分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效. 相似文献
13.
精确高效地计算近奇异积分,对边界元法的成功实施至关重要,也是边界元法在实际工程计算中面临的主要障碍之一。论文提出了一种基于球面细分技术的近奇异积分计算方法,可以精确计算任意基本解类型、任意单元形状和任意源点位置的近奇异积分。该方法首先通过计算源点到单元的最近最远距离,来确定球面细分的初始半径和终止半径;然后通过一系列半径呈指数级增长的球面来分割积分单元,得到一系列三角形和四边形子单元;最后把细分后得到的子单元变成弧形状,即三角形和四边形子单元分别变成扇形和环形子单元。由于球面细分是直接在三维笛卡尔坐标系下进行的,所以它适用于任何类型的单元。此外,由于基本解主要是源点到场点距离的函数,因此在同等精度下,近奇异积分在子单元的环向上所需要的高斯积分点数将大大减少。在径向方向上,由于球半径系列呈指数级变化,各个子块可以做到等精度高斯积分。数值算例表明,与传统近奇异积分计算方法相比,论文提出的方法更加稳定,精度更高。 相似文献
14.
边界积分方程中近奇异积分计算的一种变量替换法 总被引:2,自引:0,他引:2
准确估计近奇异边界积分是边界元分析中一项很重要的课题,其重要性仅次于对奇异积分的处理.
近年来已发展了许多方法,都取得了一定程度的成功,但这个问题至今仍未得到彻
底的解决. 基于一种新的变量变换的思想和观点,提交了一种通用的积分变换法,
它非常有效地改善了被积函数的震荡特性,从而消除了积分的近奇异性,在不增加计算量的情况
下, 极大地改进了近奇异积分计算的精度. 数值算例表明,其算法稳定,效率高,
并可达到很高的计算精度,即使区域内点非常地靠近边界,仍可取得很理想的结果. 相似文献