多极边界元法积分异性创新探析
导读：面对现实生活中工程问题的复杂性，直接求解控制方程解析解是非常困难的，这就促进了各种数值方法的飞速发展。由本站硕士论文中心整理。
第一章绪论
在现实代工程实际问题中，数值仿真技术正发挥越来越重要的作用，随着力学的深入，计算机功能的日新月异，计算机辅助设计已经成为现代工程规划设计不可缺少的手段.随着各种数值方法和软件的出现，并应用于各种工程问题中，成功的解决了工程中许多复杂问题.然而现代工程中有一些不连续、不均匀，多场祸合等关键问题依然存在，对解决工程问题造成了严重障碍，要想成功更好的解决工程问题，必须克服这一重重的困难，只有解决这些不确定的因素，才能更好的为工程服务。
由于现实生活中工程问题的复杂性，在处理工程问题时往往需要将某些工程问题转化为满足某种控制方程和一定边界条件的数学问题来处理，然而，遗憾的是直接求解控制方程解析解是非常困难的，这就促进了各种数值方法的飞速发展。有限元法曾一度在数值方法的发展史上起着极为重要的作用。但是，就像其他任何事物一样，亦存在其不足的一面，比如有限单元法处理问题时自由度数目庞大，要求计算机容量大，且计算精度低等不足，这就导致了边界元法的发展，在实际的工程问题中，两种方法相互补充，发挥各自的优势，可以经济，有效的解决实际问题。
1. 1边界元法发展历史、现状和动向
对于边界元法的基本理论研究己有一百多年的历史了，但是将边界元法作为一种数值方法提出则在二十世纪五六十年代左右，二十世纪50年代初期，Mskhel于1953年首次将边界积分方程法用于分析结构力学，同年，Kellog用边界积分方程法对拉普拉斯问题进行了求解，这便是边界元法的前身。
关于边界元法的全面研究始于20世纪五六十年代，1978年英国南安普敦大学以C. A.Brebbia "’为首的一个研究小组对60年代起在边界元方面做的工作进行了汇总，出版了((The Boundary Element Method for Engineers》一书，正式标志着边界元法作为一种数值计算方法的诞生，“边界元法(Boundary Element Method)”一词也正式被使用。
边界元法归纳为两类，即直接边界元法和间接边界元法，关于间接边界元法的概念是Jaswon, I-less, Symm和Massonnet等提出来的m。关于直接边界元法的概念，出现在Kupradze的著作中，关于直接边界元法的应用，早期工作是Rizz。和Cruse `"，用边界积分方程方法求解经典的弹性力学问题和弹性动力学问题。在这一时期，Richard Shaw对声波传播问题的边界积分方程方法进行了广泛的研究。1960年，他完成了博士学位论文，并在其后发表了两篇重要论文提出了有任意形状障碍的声波脉冲的瞬态散射问题的边界积分方程法。另外，他还对弹性动力学间接边界积分公式、三维散射问题、流固祸合问题、特征值问题、扩散问题和渐近膨胀解等问题进行了研究。1963年，Jaswon和Ponte:.r讨论了扭转问题的边界积分方程法，第一次利用了边界值和法向导数的积分关系。同年，,Jaswon对laplace方程由势理论建立了边界积分方程的数值方法，为间接边界元法的提出作出了重要贡献。其后，Jaswon等人建立了平面弹性静力学的边界积分方程，提出了数值求解的有效途径，并首次用边界积分方程法求解了板弯曲问题。1966年，Symm建立了保角映射下的边界积分方程 1969年他发现了边界积分方程在势问题包括热传导分析问题方面的应用。1967年，Rizz。运用Betti-Somi.gliang公式建立了弹性静力学问题的边界积分公式，指出了边界位移和面力的函数关系。虽然这些公式的数学理论源于Kaprudze的著作，但是Rizzo以一种简明的形式提出了与当今边界元法有着密切联系的边界积分公式。其后，Cruse, Rizzo和Shippy对这些边界积分公式进行了数值求解。对于边界元法求解实施的最大困难是边界积分奇异性的处理。Symm在70年代对二维势问题的边界积分方程中的积分奇异性问题进行了研究，并开发了相应的计算机软件。1973年，Brebbia, waston等对边界元的奇异积分进行了讨论和归纳。1975年，在他的博士学位论文中首次使用高次单元求解三维弹性静力学问题，彻底解决了边界积分方程中的奇异积分问题，大大提高了数值计算的精度，从而为边界元法的发展作出了非常重要的贡献。1974年，Cruse使用曲面元建立了三维弹性应力分析的边界积分方程的新模式，为区域的更准确描述，提高边界元法的计算精度做了重要工作。Gruse还就如何由边界面力获得表面应力、体积力如何向边界力转换、断裂力学问题以及对特殊形状的裂纹采用何种特殊的应力函数等。这些成果对现代边界元法的发展起了重要的推动作用。1977年，Symm等将直接边界元法应用于所有界面的多介质问题，是这方面最早的贡献。1978年，由英国南安普顿大学教授Brebbia编著的世界第一本边界元法专著出版，标志着边界元法基本理论初步形成，同时也确立了边界元法作为一种数值方法的地位，那以后边界元法进入了系统性的研究和发展时期。同年，在英国南安普顿大学召开了国际第一届边界元法会议。此后，边界元法国际会议每年一次轮流在世界各地举行，迄今为止已举行了好几十次。每次会议都有大量论文和专著问世，这期间边界元法的发展速度是很快的，发展水平也是非常高的。为了更好的一记录边界元法的研究成果并对边界元法的研究进展进行宣传，1984年，边界元法国际性刊物创刊，它主要就是致力于边界元法研究成果和新进展的宣传，这为边界元法的开拓和发展起了重要的推动作用。
自1978年以后的30多年来，随着计算机功能的日新月异，一代代高性能高容量计算机的陆续问世，边界元法逐渐被许多方面的工程所采用，比如应用于非线性科学工程领域的数值解析，成功地解决了如涡轮机叶片冷却孔群优化布设、内燃机曲轴强度分析、车L钢机机架强度及刚度分析、海洋石油钻井台柱防腐蚀优化设计等复杂工程问题。如今，边界元法的发展巳涉及工程和科字的很多领域，几乎可以解决所有的有限元法能够解决的问题。纵观30年来的边界元法国际国内研究成果，可以看到边界元法主要在以下几个方面得到了进一步发展和成熟:
1.数学理论的完善。
如今，边界元法像有限元法那样在收敛性、误差分析和各种不同的边界元形式的统一等方面形成了规范的数学理论。边界元法的数学理论突破了有限元法数学理论的框架而发展。
2.方法与应用的开拓。
随着边界元法研究的深入，解决各种非线性问题的边界元法将不断发展和完善，边界元法的应用范围将进一步拓宽。
3.应用软件的开发
随着边界元法理论的完善，将有更多功能齐全的边界元法通用程序包问世，促使边界元法得到更广泛的应用。边界元法与有限元法互为补充，共同发展。这一时期边界元法的发展将力求赶上有限元法的发展水平，并将与有限元法一起在工程和工业方面得到广泛应用。
我国关于边界元法的研究开始于1978年。30年来，随着国际边界元法的发展，我国的边界元法研究也取得了很大成绩。目前，我国己举办过几次边界元法国际会议、比如第三届亚太国际工程中计算方法学术会议和第九届全国工程中边界元、无网格等数值方法学术会议均在我国召开，这些会议的召开进一步促进了我国边界元法研究工作的展开和发展。我国边界元法的研究主要在以下几个方面展开并取得了巨大的进展:
一数学方面:从边界元法取得的成果看，我国研究者在边界元法数学理论方面做了大量的研究工作。主要有:中科院余得浩教授建立了自然边界元法的数学理论;重大学祝家麟教授创立了椭圆型边界积分方程的边界元法数学理论;另外其对边界元法中的区域分解方法;收敛性和误差分析等都做了讨论，比并取得了良好的成果。
二.应用方面:我国学者边界元法的应用方面了很多研究工作，我国边界元法主要在力学方面得到了广泛应用，主要有:弹性力学;弹性动力学;断裂力学;动态断裂力学;流体力学;非弹性力学;弹塑性有限变形;复合材料力学;复合材料界面断裂力学、岩土力学、有限元与边界元稿合方法和反问题等。此外，边界元法在其他方面的应用也有不少涉及，比如在板壳问题，势问题，波的传播等方面。
三.应用软件方面:伴随着边界元法的研究，对一求解以上各种问题的边界元法，我国学者都发展了相应的计算机软件，有的己经应用于解决工程实际问题，比如用FoRTRA}语一言编写的处理二维，三维问题的各种源程序。并取得了良好的效果。这些程序基本是分别针对某一类具体的工程问题而编写的，其没有完善的前后处理，在应用上显得不如有限元法通用源程序方便。还有待于进一步完善。
国内对边界元法的研究相对国外较晚。清华大学杜庆华教授于1978年开始研究边界元法，并作了广泛的宣传，促进了边界元法的发展。后来冯康、胡海昌、何广乾等一批国内著名的专家开始研究边界元，这使国内边界元法研究工作得到了迅猛的发展。
国内对于边界元法的研究工作起源于固体力学，那以后才一转入非线性领域的。边界元法在固体力学上的研究主要有，清华大学岑章志博士在其博士论文中采用了初应力形式的边界积分方程，分析了三维弹塑性问题，这在奇异积分的计算以及加速收敛方法等方面取得很好的效果。西安交通大学陈政清博士等提出弹塑性大变形边界元法论，并完成了圆柱形拉伸试件颈缩的定量数值模拟。燕山大学肖宏博卜oz’建立了三维弹塑性有限形变边界元法，并给出轧制过程边界元法源程序和板带轧制过程变形一面力一应力场，还对奇异问题做了很好的处理。
边界元法边界条件的非线性使其矛盾非常突出，如果再叠加弹塑性和摩擦非线性，求解的难度将更加复杂。另外边界元法计算效率不适应大规模工程运算课题的需要，这一直困扰着边界元法的进步。为此，提高边界元计算效率和扩大运算规模，对边界元法的发展至关重要。在边界元法求解过程中普遍注意到，整个运算量主要集中在积分的计算和方程组的求解两个方面。因此，寻求方程组系数积分运算和方程组求解的快速有效的方法是扩大边界元法解题规模的关键。在边界元法发展的初期，由于计算机条件的限制，求解问题的规模很有限。为了处理大规模问题，通常采用多子域法。即把计算区域分成多个子域，分别对各个子域建模后，通过结交面的位移连续和面力平衡条件形成块稀疏矩阵，这适合求解细长域问题。1996年R. A. Bialecki讨论了多子域的内外存交换求解器，并给出了在CRAY-YMP机上进行的4万自由度稳态热传导问题算例。
近十多年来随着计算机的迅速发展，边界元法并行计算开始应用于工程实际问题，国内外研究者结合不同的并行环境，针对不同问题对边界元法并行算法的展开了研究。如1992年E. M. Daoudi "n’和1995年M. Kreienmeyer "}'’在微机Transputer环境下进行了边界元并行算法的研究。西北工业大学A_卜天麒在SGI PowerChallenge工作站上用Krylov子空间法并行求解了不连续边界元法矩阵方程。清华大学的姚振汉课题组在基于Linux网络操作系统的PVM系统上，用8台微机并行求解了石油钻探螺纹偏心钻挺的应力场，完成了5万多自由度规模运算，所需时间为43小时41分。燕山大学申光宪课题组在基于Windows操作系统的PVM系统上对系统方程组的求解进行了并行计算，在7台节点机上成功地模拟了二辊冷轧机的轧制过程。
边界元法是吸收有限元法离散技术发展起来的继有限元之后的又一种数值计算方法，目前在国外边界元法在很多领域都得到了广泛应用，并取得了良好的效果，然而在我国尽管边界元理论方面取得了不少成就，但在应用软件方面还很缺乏，使得边界元法方面的应用还落后于一些发达国家，目前我国有少数处理二维，三维问题的边界元法实用程序。这些软件尽管也能解决一些工程问题，但是缺乏完善的前后处理功能，在应用上元显得不如有限元法软件那样方便。还有待于完善。
1. 2边界元法的一些优缺点
边界元法作为边界积分方程的离散数值解法有它特有的优点，同时也在数值上引出了一些新的问题，由于采用边界上的积分形式，降低了问题的空间维数，一个三维弹性静力学IA题利用边界积分形式则变为二维问题。平面问题，轴对称问题，柱体扭转问题以及平板弯曲问题利用边界积分形式都只需要利用一维的边界离散，即变为一维问题。三维边界元法在边界离散处理时只需要用二维曲面单元和一维线性单元，降低问题的维数无疑将使未知数数目显著降低，以弹性分析为例，对一些示范性的具体事例，为了达到一定的精度，在有限元法中常用含有数千个未知数的方程组，而用边界元法求解同样问题只需用含有一，二百个未知数的方程组就能实现。另外当边界元法应用于无限域的情况时是方便的，因为这时，为了建立域内量与边界量之间的关系，本来就己经选用了适用于无限域的基本解。边界元法可以直接用边界量作为求解变量来建立边界积分方程，也可以从虚设的源密度函数或虚设边界的虚载荷作为求解变量来形成边界积分方程这是非常方便的，然而边界元法形成的刚度矩阵一般是满阵的，求解满阵方程组计算量比较大，而有限元法一般导致对称的刚度矩阵，而且在适当的节点编号下可得到非零元素集中于对角线附近的带状稀疏矩阵，一般说来边界元在计算上具有较高的精度的特点，但是必须考虑利用最后形成的系数矩阵，其矩阵元素的求值有大量的积分运算，特别是采用高度协调的边界元时所涉及到的积分运算量是比较大的，同时还应该指出，在动力问题和进入塑性状态的计算中，除了边界积分外还有域内的积分，这样，有时一还需要将域内积分化为边界积分或者同时使用边界元法和有限元法。由边界元法和有限元法各自的优缺点可以看出，在实际工程中，有限元法和边界元法应该互为补充，取长补短，充分发挥各自的优势。更好的运用于工程实践，为其服务。
1. 3边界元法的研究进展及未来发展方向
全面了解边界元法的进展情况并看到这种数值方法的潜力及宽广的应用背景首先见于1963年对势问题的应用。经过对弹性力学和塑性力学方面的一系列问题尝试后，在七十年代中期可以看到边界元法工程应用的诸多方面。自1978年以后，每年举行的国际边界元学术会议，大体上记录了边界元法基本理论方面的研究成果和边界元技术在工程上的具体应用进展情况。国际边界元会议除了反映固体力学边界元法研究外，提到的其他方面研究有:位势流及波的传播，如工ngham于1985年的研究工作“R,，非均匀椭圆型问题，抛物线形问题，如J. Wu于1985年的工作“])J，各项异性介质，如Kobayashi于1986年的研究工作，粘弹性，粘塑性，如Kobayashi于1982年的研究工作和TELLS于1982年的研究工作。岩土地基结构和土壤相互作用，开挖，深孔的渗压固化作用也一直受到注意，如Alarcon, Kuroki于1982年的研究工。此外，电磁场的边界元分析，藕合技术及形状优化’特征值分析，平板弯曲等也受到工程界的重视。
近年来在我国举行过几次边界元法学术会议，标志着我国在边界元法方面的研究已经上升到了新的阶段。目前对边界元法的研究己经从简单的线弹性力学问题扩展到弹塑性，非线性，带时间变量等问题。30年来，边界元法在各个专门问题中的研究进展是很大的。这里的进展有Baner jee于 1980年提出的二维非线性位势问题研究，Cruse于1980年提出的二维及三维断裂分析, Shaw于1980年提出的波的边界元法‘, Rizz。和Shippy于1980年提出的热弹性问题， Tottenham,于1980年提出的板壳边界元法等以Wexle:对电机工程问题的研究，Zienkiewicz小组对有限元一边界元藕合问题的研究Banneer及Shaw的化及融化分析, Bannejee及Butterfield的多孔介质分析, Liu及Ligget、的水波分析J. C. Wu `的一般粘性流分析，Mukher jee 对金属非弹性变形与时间有关问题的研究，Kobayashi 的特征值分析和行波对于隧道空洞瞬态应力的分析，Eatockteylor 的水动力载荷的边界塑性研究，Mukher jee 对大应变大变形的塑性和粘塑性的研究及非线性断裂的研究，田中正隆对薄板大挠度的研究,， Wu,对N-S流动与时间有关的分析Liu和Liggett ’对非线性水波以及有关晃荡问题的边界元法研究，另外对于非线性位势问题，波的传播问题，特别对接触问题和流体与结构的相互作用做了专门的论述。由于边界元方法与有限元法相比有突出的优点引起了愈来愈多的科技工作者的注意，可以预计这种方法将会在各个方面得到迅速的发展。但任何一种方法都不可能十全十美，对边界元法来讲也不例外。它在很多方面还不完善，还有待于进一步研究和发展，综合边界元技术的潜在优点和目前存在的不足，可以预测边界元法在未来的研究工作中主要从以下几个方面进行:
( 1)关于边界元法的基本理论，包括对基本解、方法的收敛性及误差分析等各方面进行理论研究。
( 2)在数值计算方法和程序设计方面，如何在最经济的条件下获得高精度的解，它包括对各利‘奇异积分的有效处理和采用何种最佳方案来形成有效的计算机软件。.
(3)如何处理满的、非对称系数矩阵，探讨这类线性方程组的有效而迅速的解法。
(4)扩大边界元法的应用，使其广泛应用于诸如与时间相关的问题，非线性问题等。
(5)将边界元法应用于结构形状优化等能充分发挥边界元优势的领域。
(6)与有限元法取长补短，充分发挥其优势，使其形成一种更为有效的数值分析工具。
1. 4弹性接触问题边界元法的发展
接触这个概念对我们来说并不陌生，在现实生活中接触的例子几乎处处可见，比如在机械行业领域内机械设备，零构件之间广泛存在着接触，铰接结构中的销轴与耳片间的接触，齿轮传动中的轮齿接触、盘轴配合中的联接面间的接触，螺纹联接中联接件间的接触等。构件在接触过程中，在接触区内会产生应力集中现象，若两构件间存在相对运动，还会造成接触面间的摩擦磨损，这将促使疲劳裂纹早期形成，大大降低零构件的使用寿命。因而在机械设计中对接触问题的研究是十分必要的。对于接触问题，在大多数情况下，由于接触件的几何形状和所受载荷都很复杂，很难得到理论解，因此对接触问题大多采用数值解法。处理接触问题的数值计算方法主要有有限元法和边界元法两种，在20世纪60年代左右有限元法开始应用于处理接触问题，边界元法运用于接触问题分析是在20世纪80年代，1980Andersson等人首先用边界元法分析弹性接触问题，其根据全量理论进行推导，建立了二维无摩擦接触问题的边界积分方程，后来又从增量理论出发，导出了元增量求解的迭代过程，将边界元法推广到了二维弹性带摩擦接触问题中。Paris’采用增量迭代技术，并用不连续单元处理了接触区面力不连续问题。之后汪骏书、郑家栋分别对二维弹性无摩擦和带摩擦的接触问题进行了研究和拓展。对于三维弹性接触问题的研究，始于20世纪80年代中期，Hasegawa分析了刚性轴横向插入空心圆柱体时轴与孔之间的三维弹性无摩擦接触问题o Liu Shubi。等人用全量法推导了无摩擦三维弹性接触问题的边界积分方程，并用凝聚技术得至al仅包含接触边界上的未知量的迭代方程，大大提高了迭代效率。温卫东研究了带摩擦的弹性接触问题，利用增量理论建立了相应的边界积分方程和迭代方程，解决了带摩擦情况下切向局部坐标的选取问题。与有限元法相比，边界元法求解接触问题有如下优点:
(1)能直接求得边界上的位移和应力，数值精度较高。面上的位移和应力分布规律，边界元法只须离散边界，只对区域的边界及接触区域划分单元，不仅使数据准备工作量减少，而且在获得相同精度的情况下，使方程未知量数目大大减少，从而提高求解效率。
(2)可任意选取内部点的数目及位置，使之更适合于求解所关心的局部危险点的应力。而在有限元计算中，当单元数及高斯点数确定后，计算点的位置及数目一般就不能任意变动。
(3)用边界元法求解时，弹性体边界上的作用力和位移都是边界积分方程中的直接未知量，它们在求解时具有相同的精度。因此，接触问题最关心的接触压力与位移也具有相同的精度。
(4)边界元法为弹性体边界上法向作用力和切向作用力的祸合提供了可能性，这是我们分析无摩擦、有摩擦接触问题的关键。
如今，边界元法求解弹性接触问题的理论及算法正逐步走向成熟，并且应用于更多的工程实际问题。
1. 5边界元法源程序简介
任何一种数值方法，其应用都要通过计算机程序来实现，而作为一种应用软件，要使工程技术人员能方便地使用，就必须具有可以解决多种工程问题的大型程序包，几十年来，大量专用和通用的有限元法程序包研制成功，使得有限元法成为目前应用最为广泛的数值方法。要使边界元法像有限元法那样得到广泛应用，必须开发与有限元法程序包能够相媲美的大型边界元法程序包。目前国外已出现了一系列的边界元法程序包，但由于其结构和功能等方面的原因，还需要在推广应用的过程中不断改进和完善。
纵观现今国内外边界元法应用程序软件包，相对比较完善的主要大型程序包有BEASY, CASTOR, SURFES, EZBEA和HIBEM等几个，这些程序包中，只有BEASY含有独立于计算机程序模块的大型前后处理程序系统，其他程序包则是在计算程序模块的基础上附加以简单的前后处理功能模块而形成的，其功能，应用范围等都受到了一定条件的限制。
边界元法程序包是伴随着边界元理论的发展而发展的，1978年，第一届边界元法国际会议在英国的南安普顿大学举行，在那以后边界元法得到了突飞猛进的研究和发展，一系列计算机程序陆续面世，英国南安普顿大学D. J. Danson在1982年举行的边界元会议上宣布了他研制成功的边界元法程序包，这在国际上是首次出现的边界元法程序包，对边界元法来说，这具有跨时代的意义，标志着边界元法的研究进入了一个崭新的阶段，在随后的几次边界元法国际会议以及Brebbia, Baner jee "''等人发表的专著中陆续有关于边界元法程序包系统研究和应用的论文发表，同时也有不少其他边界元法程序包问世。从1985年起，在每年举行的边界元法国际会议上，都会有大量关于边界元法理论研究和工业应用方面的成果出现，这对边界元法计算机源程序的发展起到了巨大的推动作用，从边界元法程序包的发展可以看出其发展现状主要体现在以下几个方方面:
1.落后于边界元法的理论发展水平，目前边界元法己经运用于一些工程领域内，并成功的解决了一些工程实际问题，然而，与有限元法相比，边界元法应用的范围还很狭窄，还需要扩大应用范围，如今，边界元法对于线性问题的分析己经成熟，对非线性问题的研究也取得了巨大的进步，也逐步走向成熟。目前的大型边界元法程序包基本局限于解决静态线性问题，不能对动态，几何非线性和双重非线性问题进行处理。
2.落后于软件技术的发展水平，现今世界上知名的有限元法程序包有几十个，它们是经过长期的实践考验而得出的，功能是比较齐全和可靠的，基本上代表了当今大型工程应用软件的技术发展水平。先进的有限元法程序包具有便利的前处理和完备的后处理功能，方便的联机询问功能，重开始功能，错误诊断能力，用户子程序接口技术及子结构技术等。相比之下，边界元法在应用程序软件方面，无论质量还是数量都显得不如有限元，边界元法的程序包还需在实践中加以考验，还有待于进一步完善和处理。
1. 6本章小结
本章主要介绍了边界元法的简介以及在国际国内的发展过程，介绍了边界元多年来的研究进展和未来的发展方向，最后还对弹性接触问题的边界元法以及边界元法源程序的发展历程做了介绍，通过这些提高了对边界元法的认识，对进一步了解边界元法，研究边界元法具有十分重要的意义。
参考文献
[1].Mskhel. Corrected Boundary Integral Equation in Planar Thermoelasto-Plasticity. Int. J.Solids Struct. 1953,13(4) :356-359a
[2] Kellog. A New Version of the Fast Multipole Method for the Laplace Equation in Three Dimensions. Acta Numerical. 1953, 6:229-270
[3] C. A. Brebbia. The Boundary Element Method for Engineers. Pentech Press. London, 1978.
[4] Jaswon.Hess and symm 0困)Algorithm for Stokes AndLaplace Interaction of Particles.Phys. Fluid. 1996,8:1990-2010
[5] Kupradze,B,n(1963).memoo nomeuauan Teopuu ynpyzocru.
[6] Rizzo. Cruse A. Disc Based Block Elimination Technique Used for the Solution of Nonsymmetrical Fully Populated Matrix Systems Encountered in the Boundary Element Method. Proceeding of International Symposium on Recent Development in Boundary Element Method. Southampton University. 1978.391-404
[7] Jaswon.Ponter, M.H.Aliabadi. Out-of-Core Solver for Large, Multi-Zone Boundary Element Matrices. International Journal for Numerical Methods in Engineering.1995,38:1507-1533
[8] Symm. Crony. A Block Equation Solver for Large Unsymmetric Matrices Arising in the Boundary Element Method. International Journal for Numerical Methods in Engineering.1982, 18:997-1017
[9] Brebbia, waston. Linear Equations Solver for Large Block Matrices Arising in Boundary Element Methods. In: C.A. Brebbia, W.L. Wendland, G. Kuhn, Eds. Boundary Elements IX,Vol 1，Computational Mechanics Publications. Southampton: Springer, 1987:1256-1274
[10]岑章志.三维弹塑性分析的有限元一边界元藕合方法.〔清华大学工学博士学位论文〕.1984:36-70
[11]陈正清.弹塑性体有限形变问题的边界元法及拉伸颈缩分析.[西安交通大学工学博士学位论文〕.1987:50-86
[12]肖宏.用三维弹塑性边界元法模拟轧制过程.〔燕山大学工学博士论文」‘1993:78-89
[13]姚振汉，杜庆华.No. 4/5, pp89}93, 2001，
[14」申光宪.边界元法.边界元法应用的若干近期研究及国际新进展清华大学学报.4J匕京:机械工业出版社，1998:138-152
[15]尹欣.三维弹性问题边界元法并行计算及其工程应用.仁清华大学工学博士学位论文〕.2000:1一16
[16]杜庆华.边界积分方程方法一边界元法.北京:高等教育出版社，1989:145-169International Conference on Boundary Element Techniques(论文集)，Bei jing, 2002.清华大学出版社:77-82
[17]杜庆华.关于连续介质力学工程数值分析方法.全国有限元会议.北京，1978:352-359a
[18]汪骏书，用边界元法接弹性接触问题，第一届工程中的边界元法会议论文集，重庆，pp95}100, 1985. 10
[19〕郑家栋，解弹性接触问题边界积分方程一边界元法，固体力学学报，No. 2, pp253}257, 198
[20]姚寿广.边界元数值方法及其工程应用.北京:国防工业出版社，1995
中文摘要 5-6
ABSTRACT 6-7
第一章 绪论 10-19
1.1 边界元法发展历史、现状和动向 10-14
1.2 边界元法的一些优缺点 14
1.3 边界元法的研究进展及未来发展方向 14-16
1.4 弹性接触问题边界元法的发展 16-17
1.5 边界元法源程序简介 17-18
1.6 本章小结 18-19
第二章 多极边界元法 19-28
2.1 多极边界元法简介 19
2.2 多极边界元法的理论基础 19-24
2.3 曲面上的多极展开法 24-25
2.4 多极边界元法的原理分析 25-26
2.5 广义极小残余算法 26-27
2.6 本章小结 27-28
第三章 弹性力学问题多极边界元法 28-42
3.1 弹性多极边界元法问题的建立 28-29
3.2 弹性多极边界元法问题的基本定理与等式 29-30
3.3 弹性力学多极边界元问题的基本解 30-31
3.4 多极边界元法边界积分方程的建立 31-33
3.5 多极边界元法边界积分方程的离散化处理 33-35
3.6 多极边界元法角点问题的处理 35-36
3.7 多极边界元法表面应力的求解 36-39
3.8 多极边界元法内点位移和应力的求解 39-41
3.9 多极边界元法的数值精度分析 41
3.10 本章小结 41-42
第四章 弹性接触问题多极边界元法 42-49
4.1 接触问题的基本关系 42-44
4.2 接触问题边界积分方程的建立 44-45
4.3 接触区方程组建立及求解 45-48
4.4 本章小结 48-49
第五章 多极边界元法积分奇异性的处理 49-61
5.1 边界元法奇异积分的研究现状与动向 49-50
5.2 多极边界元法积分奇异性处理方法 50-60
5.2.1 子单元分割法 50-55
5.2.2 内蕴变换法 55-60
5.3 本章小结 60-61
第六章 多极边界元法源程序设计介绍 61-70
6.1 弹性多极边界元法程序设计介绍 61-63
6.1.1 弹性多极边界元法主要符号说明 61-62
6.1.2 弹性多极边界元法程序结构介绍 62-63
6.1.3 弹性多极边界元法源程序设计流程图 63
6.2 弹性接触多极边界元法程序设计介绍 63-66
6.2.1 主要符号说明 63-64
6.2.2 弹性接触多极边界元法程序结构介绍 64-65
6.2.3 弹性接触多极边界元法程序流程图 65-66
6.3 算例 66-68
6.4 本章小结 68-70
结论 70-72
参考文献 72-76
攻读硕士期间发表论文 76-78
致谢 78
您可能有工程硕士学位论文方面的购买需求，请到工程论文硕士论文频道选取：