边界元法在断裂力学中的研究综述
摘要：边界元法在域内采用基本解，只在边界上进行离散，代数方程组的未知数少，对应力变化剧烈的地方能得到较好计算结果。本文简要介绍了国内外利用边界元法研究断裂力学中裂纹问题的现状，并对研究中的一些关键问题进行了探讨。
关键词：边界元法；裂纹；断裂力学；特殊单元法
引言
在断裂力学中，由于裂纹尖端附近的应力场存在奇异性,以致直接应用常规数值方法分析断裂力学问题的效果往往较差,因此需要结合断裂力学的特点发展更有效的数值计算方法.
边界元法是在经典的积分方程的基础上，吸收了有限元法的离散技术而发展起来的计算方法[1]。边界元法在域内采用基本解，只在边界上进行离散，因此实际上是将问题降维处理，如果是各维尺度相近的大型问题，代数方程组的未知数将按指数规律减少，这无疑将大大减少准备工作、存贮量与机时[1]。另外，计算误差只来源于边界，区域内由解析公式计算，这就具有解析-数值计算的特点，有较高精度，对应力变化剧烈的地方能得到较好的结果，在边界上也能保持其精度，这些是有限元法所做不到的。这些特点，对边界元法应用在线弹性断裂力学问题上的应用是很有利的。
本文首先对边界元法在断裂力学中研究现状作一简介，在此基础上提出研究中存在的一些关键问题进行了初步探讨。
1．边界元法在断裂力学中研究现状
断裂力学研究的裂纹问题关键是确定应力强度因子（sif）。应力强度因子(sif)通常用来表征裂纹尖端附近区域应力场的强弱，通过它可以把构件几何形状、裂纹形状、尺寸及应力联系起来，并以它为基础来定义材料断裂的临界参数，从而把裂纹对构件断裂的影响进行定量计算。
用边界元解决裂纹问题,一般可以归纳为以下几个关键步骤:1)、建立边界积分方程；2)、选择单元模式；3)、处理裂纹尖端及其他边界奇异性；4)、实施数值或精确积分；5)、解最终线性代数方程组；6)、计算应力强度因子[2]。
样条边界元法
摘要：
1.样条边界元法的概述
2.样条边界元法的基本原理
3.样条边界元法的应用领域
4.样条边界元法的优缺点分析
正文：
样条边界元法是一种数值分析方法，主要用于求解偏微分方程的边界值问题。它以边界上的函数值为基本变量，通过构造合适的边界元函数，将边界值问题转化为求解边界元函数的线性或非线性方程组。这种方法具有广泛的应用领域，如工程技术、物理学、地球物理学等。下面将对样条边界元法的基本原理、应用领域以及优缺点进行详细介绍。
一、样条边界元法的基本原理
样条边界元法的基本思想是：首先对求解域进行边界离散，然后将边界离散点上的函数值作为基本变量，构造边界元函数。接着通过引入边界元系数，将边界值问题转化为求解边界元函数的线性或非线性方程组。最后，利用数值方法求解该方程组，得到边界值问题的解。
样条边界元法的基本原理主要包括以下几个方面：
1.边界离散：将求解域的边界进行离散，形成一组边界节点。这些边界节点通常是等分布的，可以根据实际问题进行调整。
2.构造边界元函数：根据边界节点上的函数值，构造边界元函数。样条边界元法通常采用 Bézier 曲线或 NURBS 曲线作为边界元函数的基函数。
3.引入边界元系数：在构造边界元函数的基础上，引入边界元系数，将边界值问题转化为求解边界元函数的线性或非线性方程组。
4.数值求解：利用数值方法（如有限元法、直接解法等）求解边界元方程组，得到边界值问题的解。
二、样条边界元法的应用领域
样条边界元法作为一种高效的数值分析方法，在多个领域具有广泛的应用，如：
1.工程技术：在机械结构分析、流体力学、热传导等方面有广泛应用。
2.物理学：在场论、量子力学、相对论等领域有重要应用。
3.地球物理学：在地震波传播、地壳构造、重力场模拟等方面具有重要意义。
三、样条边界元法的优缺点分析
样条边界元法具有以下优点：
1.灵活性：样条边界元法可以根据求解域的特点，灵活地选择边界元函数，提高了求解效率。
word格式-可编辑-感谢下载支持
对有限元法、有限差分法、边界元法和模拟电荷法的粗略总结:
有限元法（finite element method）：将连续的求解域离散为一组单元的组合体，用在每个单元内假设的近似函数来分片的表示求解域上待求的未知场函数，近似函数通常由未知场函数及其导数在单元各节点的数值插值函数来表达。从而使一个连续的无限自由度问题变成离散的有限自由度问题。缺点是有限元必须同时对所有域内节点和边界节点联立求解，待求未知数多，要求解的方程规模大，导致输入数据多，计算的准备工作量大。
有限差分法（finite difference method）：直接从微分方程出发，将求解区域划分为网格，近似地用差分、差商代替微分、微商，于是无限度的问题化成有限自由度的问题。这种方法在解决规则边界的问题时极为方便，但是正是由于这种限制而增加了它的局限性，即对于非规则边界的问题适用性较差。
边界元法（boundary element method）：边界元法是在有限元法之后发展起来的一种较精确有效的工程数值分析方法 。它以定义在边界上的边界积分方程为控制方程，通过对边界分元插值离散，化为代数方程组求解。它与基于偏微分方程的区域解法相比，由于降低了问题的维数，而显著降低了自由度数，边界的离散也比区域的离散方便得多，可用较简单的单元准确地模拟边界形状，最终得到阶数较低的线性代数方程组。又由于它利用微分算子的解析的基本解作为边界积分方程的核函数 ，而具有解析与数值相结合的特点，通常具有较高的精度。特别是对于边界变量变化梯度较大的问题 ，如应力集中问题 ，或边界变量出现奇异性的裂纹问题，边界元法被公认为比有限元法更加精确高效。由于边界元法所利用的微分算子基本解能自动满足无限远处的条件，因而边界元法特别便于处理无限域以及半无限域问题。边界元法的主要缺点是它的应用范围以存在相应微分算子的基本解为前提，对于非均匀介质等问题难以应用，故其适用范围远不如有限元法广泛，而且通常由它建立的求解代数方程组的系数阵是非对称满阵，对解题规模产生较大限制。对一般的非线性问题，由于在方程中会出现域内积分项，从而部分抵消了边界元法只要离散边界的优点。
- 1 - 边界元法在流—固相互作用分析中的应用
边界元法是一种相对新的数值解法，最近几十年来，越来越多的研究者开始应用它。在实际工程中，流固相互作用一直是一个挑战。它涉及到大量复杂的物理现象，包括流体流动、固体变形、换热和化学反应等等。边界元法能够有效地解决这些问题。
流固相互作用涉及到非线性的热传导、粘性流体力学和换热等。为了准确地解决这些问题，研究者必须掌握物理现象的本质，并正确估算参数，以及选择合适的有限元函数。边界元法是一种有助于准确描述复杂物理现象的计算方法。它基于弱解的离散技术，能够有效解决工程中的复杂问题。
要想有效地应用边界元法，需要正确建模问题，并考虑到外部环境的影响。为此，研究者必须将实际物理现象简化成数学模型，进而采用有限元法进行求解。由于实际系统的复杂性，有时周边环境对解决问题也会产生影响，因此研究者需要在模型中考虑这些因素。
边界元弗洛伊德模型是一种常用的并能够有效求解各种复杂流固相互作用问题的计算方法。它利用有限元方法离散描述物理现象，并使用复杂形式的迭代算法解决问题。在实际工程中，这一模型能够模拟流体流动、固体变形、换热等复杂的过程。因此，边界元弗洛伊德模型是流固相互作用分析中的一个被广泛应用的有效计算方法。
边界元法在流固相互作用分析中的应用取得了显著成果。经过多年的研究，它已被用于求解各种工程问题，如潜艇体结构的振动分析、内部流动模型、汽车发动机模型、汽车变速器模型、油藏模型、水力 - 2 - 学模型等。近年来，边界元法在流固相互作用分析中的应用也越来越广泛，可以模拟更复杂的物理现象，如固体植入体的热学性能模拟、多相流模拟、多场耦合分析等。
随着计算机技术和软件研发的不断改善，越来越多的工程问题被改用边界元法解决。计算机实验表明，边界元法是一种高效而可靠的数值解法，在大多数情况下可以准确地预测实际系统的运行情况。
