对边界法的改进（Ⅱ）

就传统边界元法在计算靠近边界的内点位（或场分量）值方面所存在的问题，提出了一种改进方法，即边界单元外移法，它将距所求内点最近的边界单元外移一适当距离，以减小基本解异性对求解稳定性和精度的不良影响。推导了外移单元上等效源密度的计算公式，给出了一个算例。     

SolvingPolyzoneProblemsofIndividualMediaUsingBoundaryElementMethod

文中在常规间接边界单元法的基础上，建立了多区域耦合边界单元法，以便求解矿山岩体中普遍存在的分区介质问题。用常规边界单元法解多域问题，要给定内外边界条件，但是，内部边界条件确定是困难的。本文试图用多域耦合边界单元法解决这一困难，即只需给定计算区域最外围边界条件，而各子域之间的边界条件借助于协调关系确定，并且各子域可用不同的材料参数     

极坐标系中具有旋转自由度的平面扇形单元

通过引进结点刚体转角，提出两具有平面内旋转自由度的平面４结点扇形单元，数值计算表明这两种单元具有很高的计算精度，为具有圆形边界的结构有限元分析提供了有效实用的单元。     

用混合边界元法解岩体不连续问题

介绍了一种混合边界元法——耦合直接边界元法与位移间断法（Ｔｈｅ Ｄｉｓｐｌａｃｅｍｅｎｔ Ｄｉｓｃｏｎｔｉｎｕｉｔｙ Ｍｅｔｈｏｄ），该方法分别采用边界单元和位移间断 单元模拟结构的边界及各种不连续构造，适于求解工程中含有软弱结构面岩 体的应力分析问题．由于在位移间断单元上只定义了位移差为基本未知量， 总体方程阶数较一般边界元法显著降低，从而节省了机时和存储量．     

对三柱体绕流流场和压力场的研究

采用边界元方法对三柱体绕流问题进行了研究,得到了速度场、内点的压力场及边界的压力分布;为了提高计算的精度,在边界积分离散化时采用线单元法,并用混合单元法对边界角点进行了处理,计算表明结果是合理的。     

边界积分方程——边界元法在三维断裂力学中的应用

本文简要地评述了边界积分方程——边界元法在三维断裂力学中的应用,也回顾了边界积分方程及边界单元的离散形式,并采用常应力插值和线性位移插值函数的三角形单元,编制了计算程序,对一个具有径向内表面椭圆裂纹的受内压圆筒进行了 K_1值的计算。计算结果与交替法所得结果进行了比较。计算表明:在三维断裂力学领域内,边界积分方程——边界元法是十分有效的,可以用比较粗糙的边界元网格得到精确的结果。计算是在 TQ-16电子计算机上进行的。     

二阶椭圆型偏微分方程有限分析解法的精度研究

该文对单元内二阶椭圆型偏微方程有限分析解法的精度作了深入研究。首先应用有限分析方法对求得单元内二阶椭圆型方程的有限分析解，其后通过泰勒级数展开方法，对单元内二阶椭圆型偏微分方程的有限分析解作了展开，引入算子将复杂方程进一步简化，证明了在任意边界近似函数条件下单元内二阶椭圆型偏微分方程的有限分析解的精度为步长ｈ的Ｏ（ｈ４）。结果为有限分析法解工程椭圆型偏微分方程问题，在解法精度上提供了保障。     

有限元法罚单元中的内外边界元方法

论证了边界元法及如何应用有限元法罚单元中的内、外边界法计算已知位移(变形)反求载荷(反力)类问题。它为工程界处理复杂的边界条件问题提供了一简单可行的方案。     

三维弹塑性摩擦接触多极边界元法

三维弹塑性摩擦接触问题是多重非线性问题，对接触区表面和塑变区的离散，须划分大量单元进行大规模运算才能获得接触位移、面力及应力场的准确信息。传统边界元法由于离散自由度所需内存大，CPU计算时间冗长，完整解题运算规模受到限制。本文在三维弹性多极留数边界元法的基础上，开发研制三维弹性数学规划型摩擦接触多极边界元法及源程序，建立三维弹塑性摩擦接触多极边界元法并研制其源程序，更新了课题组开发的原传统的三维弹塑性摩擦接触边界元法。数值试验表明，本法使计算机内存量减少近百倍，从而使细划分单元的大规模运算成为可能，并提高了计算精度。     

基于图像八叉树的三维比例边界有限元多面体网格生成算法

基于图像八叉树方法,提出了平衡八叉树和多面体网格修剪相结合的三维比例边界有限元多面体网格算法,该算法根据结构尺寸建立恰好完全包含整个结构的立方体网格,再建立结构图像像素信息,根据像素信息,按照2∶1的平衡分割原则递归地进行等分,完成平衡八叉树网格生成。结构内单元网格完全保留,结构外单元网格删除,对于结构边界单元网格,提出采用面-立方体相交判断方法进行边界单元与结构边界相交面的筛选,搜寻结构单元与结构边界表面相交点,通过有序连接相交点形成边界单元切割面,再结合边界单元其他几个面,构成裁剪后的多面体单元。数值算例结果表明,基于本文算法生成的比例边界有限元网格计算结果具有较好的精度和边界适应性。     

轴对称异质体引起波散射问题的边界元解

通过将全空间内轴对称异质体边界载荷沿周向作Ｆｏｕｒｉｅｒ级数展开，然后利用边界元法在母线方向求解边界积分方程，使弹性波散射问题的维数由三维降到一维。该方法充分利用旋转体轴对称的几何特点，采用环壳单元使计算量较普通边界元法更小，并且收敛速度较其他边界元法更快。     

轴对称异质体引起波散射问题的边界元解

通过将全空间内轴对称异质体边界载荷沿周向作Ｆｏｕｒｉｅｒ级数展开，然后利用边界元法在母线方向求解边界积分方程，使弹性波散射问题的维数由三维降到一维．该方法充分利用旋转体轴对称的几何特点，采用环壳单元使计算量较普通边界元法更小，并且收敛速度较其他边界元法更快．     

基于无厚度Goodman单元的粘弹性人工边界模拟

基于粘弹性人工边界,推导了三维一致粘弹性人工边界单元的刚度矩阵和阻尼矩阵.利用边界单元厚度对计算结果影响不大,引入了一种新的一致粘弹性单元的模拟方法,即无厚度Goodman单元.基于无厚度Goodman单元对一致粘弹性边界单元进行模拟,并通过半无限空间和均匀半空间的算例进行验证.结果表明,使用无厚度Goodman单元实现的一致粘弹性边界单元简单、实用性强,并且可以得到足够的精度.相比于有厚度的边界单元,无厚度Goodman单元符合边界单元无厚度的实际情况,理论完善,且对于可能存在的边界不规则情况适用性强.     

湍流数值模拟用于洁净隧道污染控制的研究(Ⅰ)——层流罩下工作台位置与污染边界曲线的关系

在直接求解N-S方程进行湍流数值计算的基础上,采用数值计算与实验测量相结合的方法,对单侧布置层流罩的侧回风洁净隧道内,层流罩下工作台摆放位置与污染边界曲线的关系进行了研究．研究结果表明：采用等涡线图可以非常直观地定性分析射流边界层的发展情况,进而对污染边界的发展进行预测；工作台的适当外移,有助于扩大层流罩下的洁净空间．     

面向对象的边界元程序设计

按照面向对象的程序设计方法，遵循边界元分析的本质，建立了有关描述边界元模型的类，用链表方式实现结点、配置点、边界单元和内部单元的数据存放、用多态性实现单元的自由链接，方便的实现了单元增减、复连通区域和同一程序解不同问题等功能。采用VC 编制了边界元配置法的数值计算程序，并给出了三维Laplace方程在球体区域上的算例。并将此方法的计算结果与精确解进行了比较，计算结果吻合良好。     

动力边界单元法及其在土一结构动力耦合效应分析中的应用

本文研究了动力边界单元法及其在土一结构动力耦合效应分析中的应用。文中利用Helmholtz势函数定理详细地推导了动力边界单元法的基本解,并且通过Laplace积分变换,应用直接列法建立了弹性动力的边界积分方程以及相应的数值求解公式。编制了动力边界单元法的计算机程序DYBEMI。最后对于不同激振频率和不同埋置深度的基础振动动力响应等算例进行了计算与分析。结果表明动力边界单元法在土一结构动力耦合效应问题中的分析是一种有效的方法。     

不规则边界面上离散节点三角形单元的自动连结

三维边界单元法解空间弹性力学问题,只需将边界面划分为单元。本文对几何形状不规则或边界节点分布无规律的边界面,给出了将边界节点自动连结成三角形单元网格的一种方法。     

大冲角翼栅绕流流场及压力场的边界元方法

采用粘性流场边界元方法，对大冲角的翼栅绕流问题进行了研究。得到了速度场、内点的压力分布及翼栅边界面上的压力分布。为了提高计算精度，在对积分方程离散化时，加密了刻分网格，并且用混合单元方法对角点进行处理。计算方法简洁结果合理，适应性较强，显示了本文方法的优越性。     

三维边界元单元自动剖分的一种新算法

提出了一种空间表面三角边界元单元自动剖分的新算法，给出节点及单元形成的算法及有关规则。最后，给出了典型计算实例。     

不同间距的五圆柱体绕流的流体动力特性研究

利用边界元方法对不同间距的五圆柱体绕流问题进行研究——即边界元方法在水利方面的应用,得到了不同间距的五圆柱体绕流的速度场、内点的压力场、边界的压力分布以及各圆柱体所受到的黏性摩擦阻力。为了提高计算的精度,在边界积分离散化时采用线单元法,并用混合单元法对边界角点进行了处理。计算结果合理,处理问题方便、简捷,方法适应性强。