考虑离心体力的直线单元的解析法

对于体积力问题，可将体积力的域内积分转换为边界上的积分运算[1]，这将减少域内划分单元，使边界元法处理体积力时更简秉方便，本文在文[1]和[2]的基础上，提出用单元的解析积分取代数值积分，计算离心体力问题，推导了相应的直线单元即常单元和线性单元的解析式，算例表明，本文方法数值精度高，尤其在造近边界的内点时，效果迹较好。     

三维Helmholtz积分方程外问题几乎奇异积分的半解析算法

边界元法求解声场Helmholtz外问题时,由于简单闭合曲面外的无限域边界积分方程与原边值问题(外问题)不完全等价,从而会在某些激励波数(与相应内问题的特征波数重合)下不能获得唯一解。文章引入一种计算几乎奇异积分的半解析算法,结合CHIEF点法,在较宽的波数范围内计算了声场外问题近场和远场内的声压。计算结果表明,该算法不仅有效地克服了频域内解的非唯一问题,而且与单纯的CHIEF点法相比能够显著提高计算精度。     

边界元分析薄板自由振动问题的一个近似方法

本文提出边界元法分析域内具有支承及集中质量的薄板自由振动问题的近似方法，该方法在利用基本解的基础上，将域内积分化为边界积分来处理，节省了工作量，文中计算实例结果表明，该方法的精度满足实际工程的要求。．     

弹性动力学的双互易杂交边界点法

将双互易法同杂交边界点法相结合,提出了求解弹性动力问题的新型数值方法------双互易杂交边界点方法. 该算法在求解弹性动力问题时,将控制方程非齐次项的域内积分转化为边界积分. 该方法将问题的解分为通解和特解两部分,通解使用杂交边界点法求得,特解则使用局部径向基函数插值得到,从而实现了使用静力问题的基本解来求解动力问题. 计算时仅仅需要边界上离散点的信息,无论积分还是插值都不需要网格,域内节点仅用来插值非齐次项,因此该算法仍是一种边界类型的无网格方法. 数值算例表明,该方法后处理简单,计算精度高,适合于求解弹性动力问题.      

弹性薄板弯曲问题的双互易杂交边界点法

基于一种板的修正变分泛函,将杂交边界点法与双互易法结合,用于薄板弯曲问题的分析。该方法将问题的解分为齐次方程的通解和非齐次的特解两部分,特解采用径向基函数插值得到,而通解则使用杂交边界点法求解。在杂交边界点法用于求解通解的列式过程中,边界变量采用移动最小二乘近似,域内变量则采用基本解插值。与有限元法相比,该方法仅需要边界上离散点的信息,无论插值还是积分都不需要网格,域内点仅用来插值非齐次项,因而仍是一种纯边界类型的无网格方法。数值算例表明,本文方法能以很少的计算自由度获得与其它方法同样的计算精度,且具有前后处理简单、收敛速度快等优点,适合于求解工程中各种薄板的弯曲问题。     

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

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

任意支承条件下连续板系结构的有限板块解法

本文基于薄板小挠度弯曲理论,构造出板元内部解析、边界挠度和边界法向弯矩以带补充项的付氏级数逼近、同时考虑域内多点支承作用的板元位移函数,给出了一个适用于任意支承条件下连续板系结构的有限板块法求解格式。数值计算结果表明:本文的方法具有良好的计算精度和计算效率,适于工程应用。     

任意支承条件下连续板系结构的有限板块解法

本文基于薄板小挠度弯曲理论，构造出板元内部解析，边界挠度和边界法向弯矩以带补充项的付氏级数逼近，同时考虑域内多点支承作用的板元位移函数，给出了一处适用于任意支承条件下连续板系结构的有限板块法求解格式。数值计算结果表明：本文的方法具有良好的计算精度和计算效率，适于工程应用。     

用杂交边界点法分析简支薄板的热弯曲问题

该文利用杂交边界点法对简支薄板的热弹性弯曲进行了分析计算.采用薄板的热弹性理论,通过薄板的修正变分原理建立了各向同性薄板的边界局部积分方程,域内变量使用基本解插值,而边界上的变量则用移动最小二乘法近似.计算时仅需边界上离散点的信息,无论变量近似还是数值积分都不需要网格,因此该方法是一种纯边界类型无网格方法.数值算例表明,杂交边界点法在分析薄板的热弯曲问题时具有效率高、精度高和收敛性好等优点.     

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

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

一种新型的边界元法——边界轮廓法

利用传统边界元积分方程的被积函数的散度等于零的特性，提出一种新型的边界元法——边界轮廓法，使求解问题的维数降低两维。对线弹性平面问题，选择二次位移形函数，求得相应的位移和应力势函数，使二维问题的求解转化为边界点的数值计算，给出了边界点的位移和面力及域内点的应力和位移的计算公式。实例计算表明，该方法具有较高的精度。     

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

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

薄板结构的边界元优化设计

完成了一种用边界元法对弹性薄板结构进行优化设计的计算机分析及彩色绘图系统.本文论述了平板弯曲问题的边界元法及其优化设计.系统采用本文作者提出的三方程协调方案,并精确计算分布荷载的域内积分.本文完成了多个平板弯曲问题的边界元分析算例及板结构边界元优化设计的工程算例.结果表明,本系统方法先进、结果精确可靠,具有十分明显的工程实用价值.     

边界配位法数值计算的收敛性和可行性的一些问题

本文在一种裂纹体上得到了Williams 级数的收敛性结果,依此进行边界配置计算和讨论,提出了边界配位法在收敛域内外的计算可行性问题和"主核计算域"的概念.     

弹性地基板的无奇异边界元解法

本文用无奇异边界单元法,分析弹性地基上薄板的弯曲问题,考虑Winkler和双参数地基模型,选择第三类复变量的Bessel函数作为该问题的基本解,由此导出了一组权函数,在数值解法中,对面荷载积分项均统一化为边界积分,避免了在域内划分网格,文中还给出了一种考虑桩支承的边界元分析方法及域内点弯矩和接触应力的计算公式。     

域内具有支承及集中质量的薄板自由振动的边界元分析

本文基于薄板小挠度弯曲问题的基本解,建立了任意边界条件、域内具有支承及附有征意个集中质量的薄板自山振动的边界积分方程,文中计算了若干算例,其精度是实际工程中所允许的。     

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

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

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

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

八节点等参元滑动接触算法的改进

分析认为采用八节点等参元计算滑动接触问题时,往往存在较大的接触应力计算误差的主要原因在于这种单元具有两种不同的节点。为此,本文提出了双模拟接触边界法,即分别采用两种不同的节点,形成两个模拟接触边界,并限制“从接触边界”上的节点只能和“主接触边界”上其同类节点形成的模拟接触边界形成接触关系,从而有效的减小了采用八节点等参元计算滑动接触问题时的误差。     

SBFEM与非局部宏微观损伤模型相结合的准脆性材料开裂模拟

混凝土是一种被广泛应用于土木和水利工程中的准脆性材料, 在各种内外部因素的作用下, 开裂是混凝土结构最为普遍的破坏形式, 准确模拟结构的开裂过程, 对于结构的安全评估至关重要. 将比例边界有限元与非局部宏微观损伤模型相结合提出一种准脆性材料开裂模拟新方法. 以比例边界有限元子域的比例中心作为物质点, 通过两比例中心(物质点对)之间的物质键的正伸长率来定义微细观损伤, 将某点影响域内物质键的微细观损伤加权平均得到该点的宏观拓扑损伤. 再引入能量退化函数, 将宏观拓扑损伤嵌入到比例边界有限元的基本框架中. 充分利用比例边界有限元网格允许存在悬挂节点的优势, 采用四叉树网格离散技术进行快速、高质量的多级网格划分与过渡. 通过一个I型开裂与一个混合型开裂的两个典型算例, 验证了该方法可捕获结构裂纹扩展路径与荷载变形曲线. 与现有的方法相比, 本文的损伤模型可得到更准确的局部开裂损伤带, 结果更为合理, 且具有更高的计算精度和计算效率. 当损伤过程区网格尺寸小于影响域半径的1/5时, 计算结果不存在网格敏感性问题.