三维非规则非均匀边界元网格的简便的高精度算法

对三维直接边界元中一阶奇异积分、一阶近奇异积分以及非奇异积分进行统一处理，给出了一种提高积分计算精度的简便有效的方法，对非规则非均匀边界元网格可获得比一般方法高得多的计算精度，非常适合边界形状比较复杂的三维实际问题的边界元分析．     

pFFT快速边界元方法模拟三维声散射

研究了用pFFT快速边界元方法模拟声散射问题的关键技术。采用Burton—Miller方程消除了声学边界元方法中外问题解的不唯一现象。为此,文中研究了采用常量元时该方程中超奇异积分的计算方法。最后,通过对平面声波的刚性圆球声散射的数值模拟,验证了建立的声学pFFT快速边界元方法。     

非连续边界元积分的精确表达式及相关问题

以二维位势问题边界元分析为例,给出了利用线性非连续边界元离散边界积分方程时系数矩阵积分计算的精确表达式,通过和利用Gauss积分方法计算系数矩阵所得数值结果的比较表明：配位点选择不同对数值计算结果精度影响的主要原因是积分计算的精度,尤其当配位因子选择较大时,存在的准奇异积分（Nearly Singular Integrals)很难利用常规Gauss积分方法准确求得。     

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

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

三维常数势边界元中的精确积分

对三维问题边界元方法中应用最广泛的常数边界元的积分提出一种精确积分方法。借助于一个假想的闭合曲面,将特定的势场应用于边界积分方程,发现对于三维问题,常数势项的积分可以化作球面三角型的面积计算,而导数项的积分则可在平面域用极坐标进行。本文方法结果精确,公式简单,同一计算公式可以用来计算非奇异、几乎奇异和奇异积分,统一了积分算法。     

二维边界元奇异积分和多域缩聚法分析

基于基本解的一种新的表达式，对二维边界元分析中奇异积分的精确求解进行了讨论，从几何方面对基本解的奇异性进行了分析，给出了超参非连续元离散位势和弹性力学问题边界积分方程时奇异积分计算的精确式，从而为判断各种近似方法的优劣和间接方法的精度提供了依据，也为精确地分析了大规模问题提供了一条有效的途径。     

三维有限体平片裂纹的超奇异积分方程与边界元法

利用Ｓｏｍｉｇｌｉａｎａ公式及有限部积分的概念，导出了含任意平片裂纹三维有限体问题的超奇异积分方程组，并联合使用有限部积分与边界元法，建立了数值求解方法．在裂纹前沿附近单元，采用与理论分析一致的平方根位移模型，以提高数值结果的精度．最后计算了若干典型例子的应力强度因子．     

位势问题非连续边界元分析

本文推导了二维位势问题超参非连续边界元分析的公式,证明了利用曲边单元对位势问题进行分析时,Ｃａｕｃｈｙ型奇异积分可以直接确定,没有必要采用通常的常位势方法．通过数值算例对非连续边界元分析中最优配点问题进行了说明．     

内时弹塑性力学边界积分理论和边界元计算（二）

本文在文[1]的基础上应用边界积分方程求得了球壳和简体弹塑性问题的全量解析解。首先求出其增量形式的解,然后对内时标度积分求得其最终解。与经典解比较可知本文结果是较为精确和理想的。     

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

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

非连续边界元——有限元耦合方法分析

对边界元－有限元耦合方法进行了分析，采用非连续元离散边界积分方程，解决了耦合分析中自由度约束问题，给出了非连续边界元同有限元耦合的具体实施步骤，通过对二维弹性力学和流＝固耦合问题分析，表明了该文方法的有效性。     

边界元法中计算几乎奇异积分的一种无奇异算法

边界元法中存在几乎奇异积分的计算困难。引起边界单元上几乎奇异积分的因素是源点到其邻近单元的最小距离δ。本文拓展文[1]的思想,进一步采用分部积分将δ移出奇异积分式中积分核之外,转换后积分核是δ的正则函数。所以几乎强奇异和超奇异积分被化为无奇异的规则积分与解析积分的和,可由通常的Gauss数值积分解。文中应用此正则化技术求解了弹性力学平面问题的近边界点位移和应力。     

水闸闸室结构三维样条边界元分析

用三维样边界元法分析水闸闸室结构。底板，闸墩和载水墙等为其子结构，交通桥，工作桥和胸墙等处理为内部支撑。地基和边载可以是任意的，只要能给定地表位移面力关系。在各种工况下，不论是设置平板门还是弧形门，是平底板是反拱底板，即使在稀疏剖分下也能给出高精度的位移场，应力场和地基反力场。     

基于球面细分法的高效高精度近奇异积分计算

精确高效地计算近奇异积分,对边界元法的成功实施至关重要,也是边界元法在实际工程计算中面临的主要障碍之一.论文提出了一种基于球面细分技术的近奇异积分计算方法,可以精确计算任意基本解类型、任意单元形状和任意源点位置的近奇异积分.该方法首先通过计算源点到单元的最近最远距离,来确定球面细分的初始半径和终止半径;然后通过一系列半径呈指数级增长的球面来分割积分单元,得到一系列三角形和四边形子单元;最后把细分后得到的子单元变成弧形状,即三角形和四边形子单元分别变成扇形和环形子单元.由于球面细分是直接在三维笛卡尔坐标系下进行的,所以它适用于任何类型的单元.此外,由于基本解主要是源点到场点距离的函数,因此在同等精度下,近奇异积分在子单元的环向上所需要的高斯积分点数将大大减少.在径向方向上,由于球半径系列呈指数级变化,各个子块可以做到等精度高斯积分.数值算例表明,与传统近奇异积分计算方法相比,论文提出的方法更加稳定,精度更高.     

薄壳问题的三维虚边界元解法

直接从三维弹性力学微分方程出发,依据三维的Kelvin解,应用最小二乘法建立了三维虚边界元法解薄壳问题的一般方法。本方法的显著优点是：不论求解何种壳体问题,方法的思想是不变的,均以三维的Kelvin解来建立方程,而勿需对不同几何形状的壳体采用不同的基本解。文中给出了数值算例,以作为本方法的应用。本文方法与边界元直接法相比,优点在于无需处理奇异积分,且系数阵是对称的;再者,本文方法思想简单,程序实现容易。     

三维结构在循环加载条件下弹塑性分区边界元法及其应用

本文发展了一种三维结构在单调和循环加载条件下的弹塑性分析的分区边界元法并完成了相应软件的开发。采用了改进的双面理论模拟循环加载条件下的材料性质,引进了有效的奇异积分数值处理方法和推广的返回向量法(作为双面本构理论的近似积分方法,最后,运用该软件对常用的工程构件——吊钩进行了应力应变分析,并与相应的实验结果作了较为详细的比较分析。     

边界积分方程中超奇异积分的解法

本文对边界积分方程中所存在的超奇异积分的数值解法作了综述，并介绍了它的一些应用。     

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

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