毕业论文-基于边界元的煤岩体稳定性分析

中国矿业大学本科生毕业论文姓名金红全学号10094715学院力学与建筑工程学院专业工程力学论文题目基于边界元的煤岩体稳定性分析指导教师高亚楠职称讲师二一三年六月徐州摘要在数学模拟方法中，以往只有“有限元方法”可以灵活的处理煤岩体复杂的非线性性状及地质的不连续性。近十余年来，边界元法已有很大发展，提供了求解煤岩体力学问题的另一种数值方法，无论是处理复杂的非线性因素，边界元法都具有它的特色，是一种简单而效率高的数值方法。煤层巷道的稳定性受地应力影响极大，分析不同地应力分布情况下煤岩巷道的稳定性对煤炭开采具有较大的工程意义和指导作用。本文首先介绍了边界元的发展及基本原理，随后采用边界元算法，计算了不同地应力分布情况下煤岩巷道的位移场、应力场分布及其演化规律。计算表明在地应力与水平方向呈45时，巷道变形较为明显，容易失稳，因此，在实际开挖时，应尽量避免45角方向开挖深部煤岩层巷道。关键字边界元法；煤岩体稳定性；二维弹性问题；间接边界元法；直接边界元法；ABSTRACTTHEFINITEELEMENTMETHODWASTHEONLYWAYTOWELLSOLVETHECOMPLEXNONLINEARBEHAVIOROFCOALANDROCKMASSANDGEOLOGICALDISCONTINUITYAMONGALLTHEMATHEMATICALSIMULATIONMETHODSDURINGRECENTTENYEARS,ASASIMPLEANDHIGHEFFICIENTNUMERICALMETHOD,BOUNDARYELEMENTMETHODHASBEENDEVELOPEDGREATLY,WHICHPROVIDESANALTERNATIVEMETHODTODEALWITHCOALANDROCKMASSMECHANICSANDSHOWSITSSPECIALTYINDEALINGWITHTHECOMPLEXNONLINEARELEMENTTHESTABILITYOFCOALROADWAYINFLUENCEDBYCRUSTALSTRESSISGREATANALYSISTHESTABILITYOFTHECOALROCKUNDERDIFFERENTCRUSTALSTRESSDISTRIBUTIONOFCOALMININGHASGREATENGINEERINGSIGNIFICANCEANDGUIDINGROLETHISPAPERFIRSTINTRODUCESTHEDEVELOPMENTOFTHEBOUNDARYELEMENTANDTHEBASICPRINCIPLETHENITUSESTHEBOUNDARYELEMENTMETHODCOALROCKISCALCULATEDUNDERTHEDIFFERENTINSITUSTRESSDISTRIBUTIONOFDISPLACEMENTFIELD,STRESSFIELDDISTRIBUTIONANDITSEVOLUTIONREGULARITYCALCULATIONSHOWSTHATINCRUSTALSTRESSANDTHE45TOTHEHORIZONTALDIRECTIONTHEROADWAYDEFORMATIONISRELATIVELYOBVIOUSTHEROADWAYISEASYTOLOSESTABILITYTHEREFORE,THEACTUALEXCAVATIONSHOULDAVOIDA45ANGLEDIRECTIONROADWAYEXCAVATIONOFDEEPCOALSTRATAKEYWORDSBOUNDARYELEMENTMETHODSTABILITYOFCOALANDROCKMASSTWODIMENSIONALELASTICPROBLEMSINDIRECTBOUNDARYELEMENTMETHODDIRECTBOUNDARYELEMENTMETHOD目录1边界元发展概况111发展历史1111萌芽与奠基期195019782112方法完善与初步应用期197819903113理论完善与广泛应用期1990至今512我国边界元法研究概况513岩体工程中边界元的发展62边界元法基础821概述8211边界积分方程的建立8212对于一般问题的推广10213BETTI定理、KELVIN解及SOMIGLIANA等式1122本章小结133几种常见的边界元方法1431几种常见的问题14311半平面问题14312FLAMANT问题1632几种常见的方法18321位势问题的边界元法简介18322位势问题的间接法边界积分方程19323虚应力法建立的边界积分方程21324位移间断法建立的边界积分方程2233本章小结234案例分析2441案例设计及软件介绍2442计算模型2543参数设定及云图分析2544计算结果及分析2845本章小结355结论与展望3651结论3652展望36参考文献37翻译部分39致谢461边界元发展概况11发展历史力学的问题应用到实际工程时，我们可以将其归结为求解定解问题的数学问题。不过，只有比较简单的问题可以通过数学推导得出问题的解析解，一般遇到的复杂问题都需要用实际有效的数值模拟方法来求解问题。边界积分方程边界元方法建立在有限元方法发展之后，成为了解决实际工程问题的一种有效而且比较广泛的数值模拟理论。边界元方法降低了问题的维数使这种方法更优于其他数值模拟的理论，这种方法通过设定边界上未知的一些量把他们作为基本量，根据问题的需要可以通过设定的量求解出求解域内的未知量。求解弹性力学问题时，边界元方法的解可以精确的满足求解域内的偏微分方程，因此边界元法优于有限元法的方面是边界元法的解具有较高的精度。与此同时，在某一些相关情况下，例如线弹性体的应力集中问题，有应力奇异性的弹性裂纹问题，考虑脆性材料中裂纹扩展的结构软化分析，局部进人塑性的弹塑性局部应力问题以及弹性接触的问题等，边界元方法比有限元方法更加有效、更加精确。边界元方法的这些特点，使边界元法得到了广泛的应用，并且受到了力学研究领域的重视。相对于边界元方法而言，比有限元方法有很多优势首先，在解决问题时，边界元方法通过降低问题维数使求解问题更加简单，比如把比较复杂的空间问题简化为平面问题求解，这样有效地降低了问题的难度。其二，有限元法采用的是将求解区域进行离散化，边界元方法则是只把边界区域离散，这样划分单元远小于区域划分，使得求解问题时需要求解的方程和求解所需要的数据都减少了很多，既节约了时间又降低了难度，然后，边界元方法不需要建立型函数，边界元方法直接建立在控制微分方程，以及对应的边界条件上，即使泛函无法建立，边界元方法也可以求解，而有限元就不可以了，所以这种方法可以求解无限域内的一些问题。最后，边界元法有较高的精度，因为他离散解析相结合，通过基本解来解决问题，精度提高。有优点也会有缺点，边界元方法在求解问题时需要知道基本解，也需要建立GREEN函数并且一些问题的基本解释不知道的，如非线性问题，因此，边界元方法不能求解。不过，边界元方法的使得在许多领域难以求解的问题得到解决。被广泛应用到了许多领域。而且在力学当面取得了广泛的应用，也得到了丰富的成果，边界元的发展起步比较晚，但是边界元方法的发展取得了很大的进展，因为边界元方法解决问题的优势，很多人都为边界元的发展做出了贡献。常规边界元的弱点更为明显，举例如下边界元方法求解问题时，所需要求解的代数方程是满阵，而且是非对称的，对于1000个边界点的二维问题，大概需要计算2000个自由度，这样对应的系数矩阵就有20000000个元素，其中元素是双精度的，因此，每个系数矩阵占有800M的储存空间。可想而知，边界元方法计算的难度，这样边界元方法只能够对某个关键构件进行分析计算，难以分析复杂结构。总体来开，边界元有优点也有缺点，而且相对于有限元来说，在实际应用中也发挥自己的优势。边界元有限元耦合法能够综合两者的优势，经过发展具有良好的应用前景。边界元法的发展可分如下几个时期111萌芽与奠基期1950197850年代初期，MUSKHELISHVILI1953将积分方程的方法用于结构力学的分析，KELLOGG1953用积分方程的方法求解LAPLACE的问题，这便是边界元方法的前身。现代边界积分方程方法与FREDHOLM的工作有着直接关系，他讨论了建立在离散技术上的求解方法。关于间接边界元法的概念是JASWON，HESS和SYMM等研究形成的。关于直接边界元方法，曾出现在KUPRADZE的著作中，但更多的早期工作是RIZZO和CRUSE用边界积分方程方法求解弹性动力学问题和经典的弹性力问题。在这一时期，RICHARDSHAW对波的传播问题的边界积分方程方法进行了广泛的研究。1960年，他完成了博士学位论文，并在其后发表了两篇重要论文，提出了有任意形状障碍的声波脉冲的瞬态散射问题的边界积分方程法。除此之外，他还对弹性动力学的间接边界积分公式、特征值问题、流固藕合问题、三维散射问题、渐近膨胀解和扩散问题等进行了研究。1963年，JASWON和PONTER讨论了扭转问题的积分方程方法，第一次利用了边界值和法向导数的积分关系。同年，JASWON对LAPLACE方程由势理论建立了边界积分方程的数值方法，为间接边界元法的提出作出了重要贡献。其后，JASWON等人建立了平面弹性静力学的边界积分方程，提出了数值求解的有效途径，并首次用边界积分方程方法求解了板弯曲问题。1966年，SYMM建立了保角映射下的边界积分方程。1969年他发展了边界积分方程在势问题包括热传导分析方面的应用。1967年，RIZZO运用了BETTISOMIGHANA公式建立弹性静力学问题的边界积分公式，指出了边界位移和面力的函数关系，这是文献中最早的一篇关于直接边界元方法的论文。虽然这些公式的数学理论源于KAPRADZE的著作，但是RIZZO以一种简明的形式提出了与当今边界元法有着密切联系的公式。1967年，CRUSE完成了直接边界元方法若干重要问题的推导，随后，CRUSE与RIZZO和SHIPPY配合，对这些边界积分公式进行了数值求解，相继提出了直接边界元法的若干重要论文。边界元法实施的困难之一是积分奇异性的处理。SYMM在70年代对二维势问题的边界积分方程中的积分奇异性问题进行了研究，并发展了计算软件。1973年，BREBBIA、WATSON等将边界积分方程应用于应力分析问题。1975年，LACHAT完成了他的博士论文，第一次使用高次单元求解三维弹性静力学问题，彻底解决了边界积分方程中的奇异积分问题，大大提高了计算精度，为边界元法的发展作出了非常重要的贡献。1974年，CRUSE首先使用了曲面元建立了三维弹性应力分析的边界积分方程的新模式，为几何区域的更准确描述，提高边界元法的精度做了重要工作。CRUSE还讨论了由边界面力获得表面应力、体积力向边界力转换技术、断裂力学问题以及对特殊形状的裂纹采用特殊的应力函数等。这些成果对现代边界元法的发展起了重要作用。1976年，CROUCH建议用位移不连续法DISPLACEMENTDISCONTINUITYMETHOD，IEDDM求解平面弹性问题，这是一种间接边界元法，它以单元均匀位移不连续位移分量为未知数，可以很便利地求解岩石力学问题，因而，CROUCH被公认为是间接边界元法的开创者。1977年，CRUSE就固体力学的边界积分方程法，包括直接法和间接法的数学基础发表论文，是该方面最早、全面的系统性的理论著作。同年，SYMM将直接边界元法应用于有界面的多介质问题，是非均质问题最早的具有开拓性的贡献。1978年BRADY与BRAY提出了一种四级QUADRUPOLES虚载荷用于模拟矿山薄层采场的变形，这种方法后来被确认为应力不连续法FIETITIOUSSTRESSMETHOD，IEFSM。DDM和FSM均可用于模拟裂缝或夹层，两种方法实质上相同，但DDM法更适合于裂隙、断层的模拟。边界元法BOUNDARYELEMENTMETHOD，IEBEM这一名称是CRUSE于1973年首先提出，但之后的有关文章包括CRUSE自己也没有再使用这一提法，而用的是边界积分方程法BOUNDARYINTEGRALEQUATIONMETHOD，IEBIEM。1977年，BREBBIA和BANERJEE重新使用了边界元法这个名称，边界元法从此有了明确的定义。1978年，由BREBBIA编著的第一本边界元法专著出版，对边界元法的发展有着极为重要的意义，其重要性在于它指出了边界元法与其他数值方法特别是有限元法的关系，提出了如何用加权余量法来建立边界积分方程，初步形成了边界元法的理论体系，确立了边界元法作为一种数值方法的地位，标志着边界元法从此进入了系统性的研究时期。112方法完善与初步应用期197819901978年，第一届边界元法国际会议在英国南安普敦SOUTHAMPTON大学举行。此后，边界元法国际会议几乎每年一次在世界各地举行，迄今己举行了28次。大量论文和专著先后面世，发展之快、水平之高是前所未有的。1984年，边界元法国际性刊物ENGINEERINGANALYSISJOUMAL创刊，它主要致力于边界元法研究新进展的宣传，为边界元法的发展起了重要的推动作用。从这些会议文集和各种刊物，如ENGINEERINGANALYSISJOURNAL、COMPUTERANDSTRUETURES、INTJOUMALOFNUMERIEALMETHODSINENGINEERING、COMPUTATIONALMECHANIES和COMPUTERMETHODSINAPP1IEDMECHANICSANDENGINEERING等登载的论文以及BREBBIA和BANERJEE等人的专著来看，这一时期边界元法的发展可归结为以下三个方面在数学领域虽然边界元法的发展是由计算机的迅速发展和广泛应用带来的，但也与近代数学理论的发展关系密切。边界元法在数学领域的研究，不仅解决了由于积分奇异性造成的困难，同时又对收敛性、误差分析以及各种不同的边界元法形式的统一进行了数学分析，为边界元法的可行性和可靠性提供了基本的理论。WENDIAND是研究边界元法数学理论的主要学者，其主要工作有数值积分的渐近误差分析，有限元和边界元耦合方法的误差分析，边界元法解的稳定性，弹性力学和流体流动问题边界元法的数学理论，断裂力学边界元法的误差分析等。HSIAO等人在边界元法特别是有限元和边界元耦合方法的数学理论方面作了许多工作。近年来，一些学者将有限元理论中的区域分解方法引入边界元法，讨论了其数学理论。总的说来，边界元法数学理论的研究还落后于方法和应用的研究，与有限元法数学理论的研究尚有一定的差距，有待进一步研究和发展。在方法与应用领域70年代以前，边界元法的研究仅仅限于解决以下几个方面的问题流体力学、弹性静力学、板弯曲问题、波的传播、势问题、断裂力学等，而且对于一些问题的研究也只是初步的尝试。现在，边界元法的发展已涉及很多工程和科学领域，几乎可以解决所有有限元法能够解决的问题。对于线性问题，边界元法的应用己经规范化；对非线性问题，其方法亦趋于成熟。边界元法在线性问题方面的研究和应用包括弹性力学、瞬态弹性动力学、断裂动力学、断裂力学、稳态弹性动力学、动态板弯曲问题、板弯曲问题、壳体分析、温度场和弹性热应力、壳的动态响应分析、势问题包括热传导、散射、扩散、势流、静电分析等、瞬态势问题、波的传播、稳态势问题、流体力学、流体动力学、声学、反问题等。边界元法在非线性问题方面的研究与应用已经涉及非弹性力学包括塑性、弹粘塑性、弹塑性、蠕变等、非弹性动力学、非弹性断裂动力学、非弹性断裂力学、非弹性壳体分析、非线性势问题、弹性有限变形、非线性断裂力学、含时间的非线性势问题、非线性瞬态热分析、材料非线性热分析、非线性板壳分析、岩土力学、非线性瞬态波的传插、非弹性有限变形等。相对来说，边界元法与其他数值方法主要是有限元法的耦合方法发展比较缓慢。目前，边界元与有限元耦合方法的研究和应用主要涉及以下几个方面非线性问题、断裂力学；弹塑性力学、岩土力学、势问题、流固祸合问题、弹性力学、土动力学、热分析、电力工程等。在工程和工业技术领域，边界元法的应用已涉及到采矿、土建、结构优化、机械、电力、地震、航空、地质、汽车、水工、桥梁等诸多方面。在应用软件领域边界元法作为一种数值方法，其应用要通过计算程序来实现。这种计算程序作为应用软件，是随着边界元法的发展而发展的。BREBBIA在边界元法的第一本专著中就附有简单的计算程序，这个程序为其后许多计算程序的研制建立了模式，为边界元法应用软件的发展和边界元法的应用都起到了良好的推动作用。1978年以后，随着边界元法国际会议在世界各地逐年举行，陆续有边界元法应用软件的新成果问世。1982年，在第四届边界元法国际会议上，英国南安普敦大学的DANSON介绍了他们研制的边界元分析程序包BEASY，这是国际上第一个边界元法大型软件。1985年以来边界元技术国际会议在世界各地举行，它着重于边界元计算技术的研究和应用、包括工程应用、计算技术和工业应用等，为边界元应用软件的发展起到促进作用。但是，随着计算机技术的迅速发展，计算机软件已成为商品，稍稍复杂一点的计算程序都不会无偿地在文献中发表。现在，以边界元法为内容的部分书籍和文献中附有简单的程序，是以教学为目的，着重说明边界元法的基本理论和方法，供读者学习边界元法时参考，只具有初等实用价值，所解的问题也只是简单的线性问题。1988年，MAEKERLE和BREBBIA在文献中从软件的来源、类型、应用范围、前后处理、元素库、材料性质、特殊功能和硬件准备等多方面对135个边界元法应用软件进行了归纳。基本上反映了这一时期边界元法应用软件的发展水平和趋势。可以看出，这一时期边界元法应用软件的发展取得了一定的成绩，但与应用于各个领域的边界元法本身的发展及计算机软件技术的发展尚有一定距离。现在，边界元法应用软件己由原来的解决单一问题的计算程序向具有前后处理功能、可以解决多种问题的边界元法程序包发展，已经形成的较大程序包有BEASY英国、CASTOR法国、BETSY德国、SURFES日本、EZBEA美国等。其中，BEASY在当今边界元应用软件市场可谓独占鳌头，在英国、美国、法国和日本等国的大学、研究所和公司得到了一定的应用。但是，可以说，边界元法程序包无论在质量上，还是数量上均与有限元法程序包的发展有一定差距。要想使边界元法象有限元法那样得到广泛应用，还必须发展各种各样的高质量的通用程序包。113理论完善与广泛应用期1990至今90年代以后，边界元法在以下几个方面得到进一步发展（1）数学理论的完善边界元法像有限元法那样在收敛性、误差分析和各种不同的边界元形式的统一等方面形成较规范的数学理论。（2）方法与应用的开拓随着边界元研究的深入，解决各种非线性问题的边界元法己有不同程度地发展和完善，边界元法的应用范围进一步拓宽。（3）应用软件的开发随着边界元方法和理论的完善，已有更多功能齐全的边界元法通用程序包问世，促使边界元法得到更广泛的应用。预计不久的将来，边界元法的应用将与有限元法互为补充。但尚难以并驾齐驱，有待进一步发展。12我国边界元法研究概况我国的边界元法研究起步于20世纪70年代末期，清华大学的杜庆华在推动我国的边界元法研究方面起了重要作用。在他的创导下，从1985年开始每3年举行一次全国工程中边界元法叙述会议。前6次依次在重庆、南宁、武汉、南京、北京、重庆举行，第7次于2003年在秦皇岛和全球华人边界元和无网格法学术讨论会共同举行，第8次于2006年在合肥与亚太国际工程数值方法学术会议共同举行，第9次于2009年10月在南京与亚太国际工程数值方法学术会议共同举行。由杜庆华和日本境界要素法研究会的田中正隆（MTANAKA）和小林昭一（SKOBAYASHI）发起，从1987年开始每3年举行两次中日边界元法学术会议。举办的依次为日本的轻井泽（1987），北京清华大学（1988），日本东京的八王子（1990）北京清华大学（1991），日本北海道的札幌1992,上海同济大学（1993），日本九州的福冈（1995），北京清华大学（1998）。除此之外，1986年杜庆华在清华大学组织召开了一次边界元法的国际会议，2002年，姚振汉和MALIABADI联合主持了第3届国际边界元技术学术会议（BETEQ2002）。20世纪八九十年代国内在边界元法及其应用研究方面成果较突出的有清华大学工程力学系的杜庆华、姚振汉、岑志章的研究组（固体力学线性、非浅性问题）；西安交通大学工程力学系的嵇醒、楼志文等（固体力学线性、非浅性问题），嵇醒后来转到同济大学；上海交通大学工程力学系的张永文（断裂力学），鲁船敬（流体力学）等，华中理工大学工程力学系的黄玉盈（流固耦合）、钟伟芳（弹性波传播）等，武汉大学数学系的冯振兴（可动边界问题）；西北工业大学的叶天麒（固体力学）；中国科技大学黄茂光的研究组（板壳问题）；合肥工业大学土木系王有成的研究组（板壳问题）等。中国科学院数学与系统科学研究院的余德浩则师承冯康，在自然边界元法的数学理论方面成果突出。我国先后出版了多本有关边界元法的专著，其中面向工科专业的代表性的有杜庆华等1989年在高等教育出版社出版的边界积分方程边界元法的力学基础与工程应用，嵇醒等1997年在同济大学出版社出版的边界元法进展及通用程序等；面向数学专业的有祝家麟1991年在科学出版社出版的椭圆边值问题的边界元分析，余德浩1993年在科学出版社出版的自然边界元方法的数值理论等；与工程密切结合的有申光宪2005年在科学出版社出版的多极边界元方法和轧制工程等。此外还有秦荣关于样条边界元法的专著、孙焕纯关于无奇异边界元法的专著、田宗若的有关复合材料分析的边界元法的专著、杨德全的包含流体力学边界元法内容的教科书，以及一些针对一定应用背景的书籍等。大陆学者在边界元法的工程应用方面做出突出贡献的可以提到清华大学张楚汉运用动力边界元法与断裂力学原理提出重力坝地震断裂与拱坝裂缝扩展模型，在边坡与地下工程研究方面提出了时域边界元与离散元耦合模型；太原重机学院黄庆学、燕山大学申光宪等将边界元法用于延长型轧机轴承寿命研究，取得突出成果。与国际同行相比，我国边界元研究较大的差距在于边界元软件的开发及其在工程应用中的推广，国内至今还没有得到较广泛的实际应用的边界元法应用软件。13岩体工程中边界元的发展岩体工程领域遇到的问题比较适合用边界元的方法，对于一般的岩体工程问题，设定好边界条件后就可以用边界元方法进行解决，但是，一般情况下，普通的边界元方法只可以解决理想的弹性问题，并且要求岩体的材料是各向同性的，然而，在实际情况下，岩石是不能满足理想状态的，因此，边界元方法不能广泛的应用于岩石工程领域，经过几十年的发展，边界元方法只是在理论领域有了较大的发展，而在工程实际应用方面相对落后，即使有研究实际工程问题的文献，也只是简单的介绍，所以边界元的发展在工程领域较慢。不过也有一些学者有较大的贡献。比如，姚振汉与杜庆华总结边界元应用在最近发展的研究中介绍了二维问题含缺陷弹性体移动接触情况下的边界元方法，而且简单介绍了边界元快速算法，使边界元在岩石工程领域展示了广阔的前景，一些学者研究洞室开挖等实际问题并且编写相应的软件，也有学者研究巷道稳定性，并进行实际参数的测定，通过边界元计算软件进行计算得出正确的结论，例如，曾留伟将边界元方法与电脑计算软件相结合，求解了三维状态下地下洞室的开挖问题，找到了洞室变化的规律，经过总结并加以整理得出了洞室围岩的稳定性。刘兴业利用边界元方法对围岩的应力以及变形进行研究，给出了测定围岩稳定性的几种方法。虽然，边界元方法在工程领域发展迟缓，经过近期多位学者的努力，也取得了较大的成就，系统的得出了一些计算方法和基本理论，为以后的发展铺平了道路。2边界元法基础21概述本章通过介绍边界元法的基本思想及其基本公式，结合地下矿山实例进行模拟计算与分析，获得了满足采矿工程问题需要的、有实际意义的结果。借助计算机软件ROCSCIENCE软件，通过对深部煤层巷道的岩体边界稳定性进行分析，通过选择施压角度、各点线受力分析、水平及垂直方向的应力位移，采用数值统计得出受力位移曲线，进而对煤岩体稳定性进行分析。物理问题中，我们经常遇到的是位势相关问题，也就是我们经常求解的调和方程，在求解调和方程时，我们通常建立所求问题的边界积分方程，然后再将建立的边界积分方程进行推广，本章节，我们首先推导边界积分方程，然后将这种研究方法推广，在一般的偏分方程方面进行研究，最后，以位势问题作为我们推导边界方程的实例，通过位势问题介绍边界元方法。本章设计思路首先建立边界积分方程，然后对一般问题进行推广，最后介绍BETTI定理、KELVIN解及SOMIGLIANA等式，从而对边界元方程进行理论研究。211边界积分方程的建立首先利用空间域或者上任意一点的调和方程的解，V积分方程为（即边界变量的表达式）,SSPQPQDSPVN,SSIIINQ（21）这样我们就建立了边界积分方程，只需要将奇异点P（基本解的奇异点）从域或者内趋于边界点P，如图11所示。V于是得到,SSPQCQDSPN（22）或者用指标符号写成,SSIIICPPQPQDSQP（23）其中,00LIMLIMSSIIPQCPDSPQNDSN（24）是如图11所示为把不满足GREEN等式应用条件的奇异点从积分域中除去而S作的以为半径的球面，是小球面各点指向P点的外法线方向。NQ不难得证，对于光滑边界，而对角点P则，为角1/2CP/4C点处在域内一侧度量的立体角（即内角）的弧度数。如果对最一般的情况把边界S分为给定函数值的部分和给定法向导数的部分S，则（22）式还可以改写成GSGSSSSSSPQPQCPPQDPQDSNN（25）其中和为边界给定量。对于边界每点有一个边界未知量上的或者G上GS的。而这一边界积分方程的奇异点也可以选取在边界上的任意点，对于这样QP的问题，在数学上是是定的。建立了这样的边界积分方程，这种方法通过设定边界上未知的一些量把他们作为基本量，根据问题的需要可以通过设定的量求解出求解域内的未知量，这样我们逐步进行求解通过设定边界上未知的一些量把他们作为基本量，根据问题的需要可以通过设定的量求解出求解域内的未知量。通过推导我们可以看出，调和方程的边界积分方程可以从GREEN等式出发建立起来。从更有普遍意义（即适用于许多不同问题边界积分方程的建立）的角度来看，建立边界积分方程的过程可以看成算子的基本方法不需建立权函数作为算子基本解，按照加权余量法的基本方法建立边界积分方程；建立方程后应用GAUSS方法把边界积分和域内积分进行相关链接，进而得到，求解域内域点P的解，然后我们再用边界变量表示；再得到方程最后将基本解点P（所谓奇异点）趋于边界点，得P到边界积分方程。这个过程可以简称为用加权余量法建立边界积分方程。对于二维位势问题，可以用同样的方法由二维LAPLACE方程出发来导出相应的边界积分方程。得出域内解的积分公式为SSPQPQDN（26）相应的边界积分方程为SSPQCPPQDN（27）或者写成GSSSSPQPQCPPQDPQDNN（28）其中用表示二维域的边界，基本解为1LN2,SAPQRPQ（29）其中是任意选定的一个距离标尺，具有长度的量纲，数值可以是1或者根据A实际问题的几何尺度选定某个数。这里的之所以要有长度量纲是因为对数是个超A越函数，该函数的自变量必须是无量纲的。将（28）式用指标符号表示，即,GSSSSCPPQPQNDQPQPNQD（210）其中，对于二维问题的指标用希腊字母表示，其取值范围为1,2。212对于一般问题的推广对于二维域上的一般二阶线性偏微分方程0,LX（211）其中微分算子为LABC（212）即,0,LABCXA（213）引入权函数，按照加权余量格式，令,LDD（214）利用GAUSS公式可得,LDANBANDABCD（215）其中可定义,ABCM（216）微分算子称为算子的共轭算子，一般情况下为了进行相应于前面所述的推ML导应使权函数满足方程P（217）即取共轭算子的基本解为权函数。对于特例ML（218）称算子为自共轭算子，对于二阶线性微分算子而言，自共轭的条件即是LA常系数，且。显然调和算子就是一例，因此可取算子的基本解作为权函数0BL来推导边界积分方程。对于非自共轭情况，应取共轭算子的基本解为权函数。213BETTI定理、KELVIN解及SOMIGLIANA等式为了以比较直观的方式来建立弹性静力学的边界积分方程，可以将BETTI功互等定理（简称BETTI定理）作为推导的出发点，我们简单介绍一下互等定理，为后面的推导提供理论依据。弹性体一般有2种平衡状态，造成的，和他的，造成的。1IFIT1IU2IFIT2IU接下来我们计算，在消耗的功，我们可以采用平衡方程或者GAUSS公式IFITIU就可以求出相应的值21212121,IIIJIIJVSSVFUDTDNDUD,IJIJVJU（219）2121,IJIJVVDD利用应力应变的关系可以写出211212,IJIJKLKLIJEU（220）于是又（219）式不难得出21212112,IIIJIJVSVVFUDTUDDDIISFUTU（221）此式即为BETTI定理，它可以叙述如下假如同一弹性体承受两组体积力和表面力的作用，那么第一组力，在由第二组力所引起的位移上所做的功等于1IFIT2IU第二组力，在第一组力所引起的位移上所做的功。在实际应用中，这两2IFIT1IU组荷载与相应的变形状态中通常一组是待求的真实状态，而另一组是为求解方便而引起的辅助状态。我们为了弹性静力学的边界积分方程，需要在点处作用一个单位集中力，而P这个力引起的辅助状态变形后，满足一定的方程为，这个方程为,0SSIKJIJKIJGUPQUPQ（222）这个辅助问题称为KELVIN问题，是有经典的解析解的。设集中力沿方向作用于坐标原点（图22）。它FOZO在边界上满足的条件是原点处应力奇异性与集中力是保持一致的；所有的应力分量在无穷远处为0。我们对集中力的解释为取一个小球洞，在其表面上有一个载荷系，这个载荷系的极限就是集中力。由于该问题的轴对称性质，可以在圆柱坐标系下采用LOVE应变函数求解。得到的解为,Z23321,0,1616ZFZFZUUGRGRR（223）推广到一般情况并采用指标符号，在我们所取的任何点（称源点）沿方向PIX作用时，对于单位的力，在空间域内的任意一个点（称场点）处引起的方向位QJ移分量可表示为,1346SIJIJIJUPQRGR（224）其中222133RXPXQPXP（225）这是点与点之间的距离。在无穷远处，在点，应力和位移都是有奇异性QPP的，并且我们所求解的应力状态为0，位移也趋于0，并且应力奇异性与集中力相似，在空间域内的任意一个点（称场点）处引起的方向位移分量。JX我们如果假设在互等定理的等式（221）中取得KELVIN解（225），并且他的变形状和我们前面设定的是差不多的，则2SJIJUQP（226）紧接着我们有限域内推导待解问题，可以看到点在域内，点为域内的VPV任何一点。对于这种状态我们可以得到而且在点处仅有沿方向作用的单位集IX中力，域内不可能有分布的力，而在边界则应该作用与（224）S式的位移场相对应的面力，若将边界上任意场点记作，则可得2JTQ2SJIJTQTP,213128IJIJIJJIRRNRN（227）代入（221）式可以得到（22811SSIIJJIJJVSUPPQFDVUPQTDS1SIJJSTPQUDS）对于待解问题的真实解略去上标（1），即可写成SSIIJJIJJVSFT（229）SIJJSTPQUDS这就是弹性理论的SOMIGLIANA等式，对无体积力情况则简化为SSIIJJIJJSSPTQTPQUDS（230）由（229）式我们得到了以下几种结论，弹性理论的解的性质大都可以这样确定边界各点的位移与面力全部已经确定，则域内任意点的位移也都随之确定。IUIT这是求解弹性力学问题的理论方法，即把边界上未知的量全部求解出来。对于边界元的方法我们在下一章进行讨论22本章小结通过推导我们可以看出建立边界积分方程，然后对一般问题进行推广，最后介绍BETTI定理、KELVIN解及SOMIGLIANA等式，从而对边界元方程进行理论研究，这样求解一般问题是，求解思路与方法就比较清晰。本章在列出弹性力学的微分提法即其偏微分方程边值问题的基础上，介绍利用赋予力学意义的数学公式（即为公式化的力学规律）来推导边界积分方程由BETTI功互等定理出发，利用KELVIN基本解，导出SOMIGLIANA等式，最终得到弹性力学的边界积分方程。3几种常见的边界元方法31几种常见的问题下面简单介绍几种边界元方法适用的问题311半平面问题在平面应变问题的各向同性弹性半平面内任意一点作用单位集中力所引起的位移和面力可表示成KELVIN解与辅助解之和，即,SKSCSKSCUTT（31）叠加辅助解后满足半平面的直线边界无面力作用的边界条件。其中下标代表在点作用方向的单位集中力引起的点方向的位移分量。PXQX的具体公式为SCU222111412121241212214223483LN32348LNSCDSCDSCDSCDRCXKRCXRRUUKR224RCXRR（32）其中（参看图31）1122/1/121,0ARCTN,8GDRXQPXQPRXXK相应的面力分量为CSCSTN（33）为与上述位移辅助解对应的应力辅助解，其具体公式列出如下CS2221C11246RC13XCSSXRCXRRK2221C112246SSSXCXCXRRRR2221C112246RCX3CSSRXCXRK221C212466SSCXRCXRRR2211CC1212463XCSSSCXRCXRRKR2221C122461SSRCXCXRR（34）其中。1/4K当点为半平面边界线上的点时，上述由KELVIN解和辅助解叠加所得基本P2OX解即FLAMANT解，其公式可以写成1121LN,SFDUR22SK11,SFDR221LNSU,SFRTR（35）其中。这种解可用于上有外力作用的部分。1/2DKG2OX相应的边界积分方程为000,SSSCPUPQFDUPQTDTPQUDP00SFSFSFQ,SSUQTTU（36）其中为半平面内部的边界，例如空洞边界，0为半平面边界上有面力作用的部分。如图32以上对于半平面问题的基本解公式都是针对平面应变问题给出的，对于平面应力问题只要将公式中的泊松比改为即可。对于半空间问题，/1如果在板空间的平面边界有部分受外力作用，也可以做类似处理。312FLAMANT问题我们首先建立一个各向同性的半平面，在这个半平面上作用一个集中力F，这样的问题就是我们所研究的FLAMANT问题，这一节我们就来用这种方法下的奇异解构造边界方程的边界元方法。这种问题一般应用到固体力学领域，对于FLAMANT问题的说明我们可以通过右图可以看出FY代表沿Z轴作用的线荷载，而它的单位为N/M。这样的求解问题，我们可以看成是平面求解问题，很显然，边界元法降低了求解问题的维数。它的求解域为。,0XY可以得到在时半平面的应力为0Y（37）位移为（其中是常量）LYX2FU12ARCTNGYXYX122YLY（38）得到在半平面边界上，向量的分量为和，的这样IJITNXYXYTYT我们可以得到单位外法线向量为，而它的两个方向的分量分别为为和J0XN。通过方程37我们可以得到，在原点不可以求解外，对于半平面上1YN0XY的任意一点，应力大小为0。也就是说，在原点处，半平面受到的力是单一的。也就是前面所说的集中力。这样我们就可以建立下面的函数该函数定义如下ARCTNTAYYARKXX（39）为了得到在方程（38）中位移分量，我们就必须得到反正切函数的值。ARCTN/YX其中表示的主值，有TANYARCXRCTAN/2T/2AYX（310）K由变量X和Y决定，即1,XYK当当当（311）因此的值位于和之间。对半无限区域，该函数的值在ARCTN/YX0Y和0之间，即有，ARCTN/0XY（312）特别的，当Y0时，有XYYX如果如果（313）把这些结果代入（38）的第一个方程可得12,04XYUFXG（314）可以看出半平面表面位移在X方向的分量为一个定值，不随着方向去改变。也就是说，（如图31中所示）半平面上的点都会向着偏离O点的方向移动。0YF相面对38中的另外的一个方程进行讨论，令，可以得到0Y1YYUFLNXG（315）在（315）方程中得到位移在力作用点为正无穷大，在远离这个作用YU0Y点时越来越小。从而能够看出边界上，施加载荷时的一些性质。虽然我们是在平,XY面设定的一些平面应变条件，平面应变条件适用于法向平行于轴的任意平面，但Z是线型力在方向是无限长的，我们是没有办法改变的。Z对FLAMANT问题，方程315中的是用来量化值的常量。由方程315知，LYU当时。,0XLY0YU（316）32几种常见的方法下面简单介绍几种边界元方法适用的方法321位势问题的边界元法简介边界元法是基于边界积分方程，采用与有限元法类似的分元、离散思想建立起来的。由于边界积分方程是定义在边界上的，因此边界元法只要在边界上划分单元。以二维位势问题为例，首先假设把整个边界划分为个单元，即在每个单EN1ENJJ元上对边界变量采用一定的插值函数（也称形函数）插值，例如可采用只保证相邻单元间未知量本身连续的LAGRANGE插值11MLLN11MLLNN（317）其中是单元上的局部坐标，一般情况下；M为单元节点数，L为单元,内节点序号，为相应的插值函数。需要提到与一般有限元不同的是，这里还可MLN以采用相邻单元间未知量本身都不连续的常值单元。原则上说，只有未知量是必须利用型函数插值的，边界给定量可以不进行插值离散化，而直接利用精确给定函数。此外，一般说来对不同边界未知量还可以采用不同的插值规律。但为简单起见，下面不妨对未知量、给定量两者都进行离散插值，且对函数值和法向导数值采用同样的插值函数。将（317）式代入边界积分方程（27）得1EJNMMSLLJLCPQPQDN1EJSLLJLNN（318）将离散插值以后的弱解代入此方程，将得到误差1EJNMMSLLJLPCQPQDN1EJSLLJLN（319）采用加权余量法的配点格式，即权函数采用各节点处的DIRACDELTA函数，得,0,12,IIPDPN（320）其中为总节点数。由此即得N1EJMMSLLJLCPQNPQDN10,12,EJNSLLJLINN（321）注意到单元号，单元内节点号和节点整体编号的关系，并区分节点未知量JLI和节点给定量，实际上（321）式刚好是对于个节点未知量的个线性代数方程。N方程的未知量系数给定和给定量系数中包含有积分JMSLNPQDJLN（322）它们叫做核函数与形函数乘积的积分。核函数即与基本解有关的量，如，它们是离散化前积分方程中的积分核。这些积分一般用数SPQSPQN值积分法积分，而且当核函数的奇异点在积分域上时这些积分是奇异积分。奇异IP积分在计算式要作特殊处理，具体处理方法在后面关于边界元数值方法的后面章节进行介绍。一旦通过数值积分计算出方程组的系数矩阵和右端列阵，可得一满阵的线性代数方程组。由此解出边界节点未知量，也就是解得了边界未知量。由于这里的，同时作为互相独立的边界变量，因此得到的的，具有NN大致相同的精度。当要求域内点的未知量时可以根据域内未知量用边界变量表示的积分表达式通过数值积分求出。322位势问题的间接法边界积分方程在此之前介绍的都是直接法的边界积分方程，直接以原始物理问题或微分方程边值问题中的边界未知量作为基本未知量，即边界积分方程的求解未知量。直接法是一种直观且有效的方法，但不是唯一的建立边界积分方程的方法。下面将简要介绍几种别的方法。位势问题的求解是要求出满足边界条件的调和函数。根据数学的位势理论，调和函数可用在边界连续分布的源所产生的单层势或双层势来表示，而根据满足条件的要求就可确定这种源在边界的分布，从而把要求的调和函数完全确定。为了得到求解NEUMANN问题的积分方程，可以假设所求的调和函数表示为源密度为的单层势QPDP（323）其中的核函数是对数势，也就是二维问题的基本解PSPQ1LN2,RPQ（324）对于三维问题，（323）式中的应改为，核函数是NEWTON势SP14,PRQ（325）也就是三位位势问题的基本解。SPQ（326）将（323）式在边界场点对外法向求导、再代入边界条件，即得Q12PQDPGQN（327）这就是对源密度的积分方程，由于是定义在边界的，所以这也是边界积分方程。对三维问题只要在此方程中将改为。S为得到DIRICHLET问题的边界积分方程，经典的位势方法是假设可以表示为密度为的源的双层势PQDPN（328）将（328）式中的内点趋于边界点，再代入边界条件，可得Q12PQQDPN（329）这就是对于源密度的边界积分方程。对于三维问题只要把改为即可。S方程（327）和（329）都是第二类FREDHOLM积分方程，由此将源密度或解出之后，就可由（323）和（328）式求域内任意点的。此外对于DIRICHLET问题也可采用单层势（323），得到边界积分方程QPQDP（