1 基本概念
1.1 定义与核心思想
边界元法是一种求解偏微分方程的数值方法，其基本思路不是对整个求解区域进行离散，而是将问题转化为边界上的积分方程，再对边界进行划分和计算。由于只需处理区域边界，所需自由度通常明显少于体离散方法。
这种方法的关键在于借助问题的基本解，把场变量表示为边界量的组合，从而把原本的空间问题压缩为低一维的边界问题。对于二维区域，计算对象通常落在一维边界；对于三维区域，则落在二维曲面上。
1.2 与其他数值方法的关系
边界元法与有限元法、有限差分法同属偏微分方程的离散求解工具，但侧重点不同。前两者通常直接对区域内部进行网格化，而边界元法更强调边界信息的利用。它常被视为一种基于积分表示的数值分析框架。
在实际工程中，边界元法既可以独立使用，也可以与其他方法结合，以弥补各自的适用范围与计算特性差异。
1.2.1 与有限元法的比较
有限元法擅长处理材料非均匀、几何复杂以及强非线性问题，因为它对区域内部的局部建模能力较强。边界元法则在均匀介质、无限域或外部场问题中更具优势，尤其适合只需关注边界响应的场景。
两者在离散方式上也有明显差别。有限元法需要对体域剖分，未知量分布于内部单元节点；边界元法只需在边界上布置单元和未知量，因此模型规模往往更小。不过，边界元法通常会形成全局稠密矩阵，计算与存储代价可能随规模快速上升。
1.2.2 与有限差分法的比较
有限差分法主要通过网格点上的差商近似导数，结构直观，便于在规则区域实现。边界元法则不依赖对整个区域的规则网格化，而是利用积分方程直接关联边界数据与场量。
相较之下，有限差分法在规则网格和简单边界上更高效，但对复杂边界适应性较弱。边界元法对几何外形的表达通常更自然，尤其在外场传播、辐射与远场衰减等问题中，表现更为突出。
1.3 适用问题类型
边界元法适用于线性、均匀介质中的边值问题，尤其是外部域、半无限域和无限域问题。例如，静电势、稳态热传导、弹性位移场、声辐射和电磁散射等，都是其常见应用对象。
当问题的控制方程存在已知基本解，且边界条件能较方便地写入积分形式时，边界元法更容易发挥优势。若系统包含强烈非线性、材料参数显著变化或体力项复杂分布，则通常需要额外处理或与其他方法联合使用。
2 数学基础
2.1 偏微分方程与边界条件
边界元法所处理的对象通常是椭圆型、抛物型或双曲型偏微分方程。为了使问题具有唯一解，必须配合适当的边界条件，如狄利克雷条件、诺依曼条件或混合边界条件。
在物理解释上，边界条件决定了已知量与未知量在边界上的分布方式。例如，在位势问题中，边界上可给定势值或法向通量；在弹性问题中，则可指定位移或表面力。边界元法的数学建模往往正是围绕这些边界量展开。
2.2 格林函数与基本解
格林函数和基本解是边界元法的重要理论基础。基本解表示单位点源在无限空间中引起的响应，它是构造积分表示的核心工具。借助基本解，可以把区域内任一点的场量写成边界上的积分形式。
格林函数则可看作满足特定边界条件的响应函数。实际求解中，基本解常用于建立边界积分方程，并将原问题的微分形式转化为积分形式，从而为边界离散奠定基础。
2.3 积分方程理论
边界元法本质上依赖积分方程理论。通过对控制方程施加加权函数或利用格林第二恒等式，可把区域内的偏微分方程转换为边界积分方程。这样一来，未知量不再分布于整个体域，而主要集中在边界上。
这类积分方程通常包含核函数、边界未知量以及已知项。其求解可形成线性代数系统，之后通过数值线性代数方法得到边界变量，再进一步恢复区域内的场分布。
2.3.1 直接法与间接法
直接法是最常见的形式，它直接以位移、速度、势函数及其法向导数等物理量作为未知量建立方程。这种方法物理意义清楚，便于与工程边界条件对应。
间接法则引入虚拟密度或势层分布，通过单层势、双层势等表示构造解。它在某些理论推导和特殊边界处理上较灵活，但未知量的物理含义通常不如直接法直观。
2.3.2 单层势与双层势表示
单层势表示把边界看作分布源的集合，源强沿边界积分后产生场响应。该表示形式常用于位势和传热问题中，结构相对直接。
双层势表示则可理解为沿边界分布一组偶极子或法向导数型源项。它在处理某些边界积分公式时非常重要，尤其在构造第二类积分方程时常见。实际应用中，两者也可结合使用，以增强数值稳定性和灵活性。
2.4 边界积分公式
边界积分公式是边界元法的核心表达式之一，它将区域内点的场量与边界上的值联系起来。对于不同类型的控制方程，积分公式的具体形式会有所差异，但基本结构通常由基本解、边界变量及其导数构成。
通过该公式，求解域中的任一点可以由边界数据重建。这种“由边界推出内部”的思路，使得边界元法在外场计算和远场预测中具有天然优势。
3 离散化与求解过程
3.1 边界划分与单元选取
边界元法首先需要将边界划分为若干离散单元。二维问题中常用线单元，三维问题中则常用面单元。单元形状可根据几何复杂程度选取为直线段、三角形、四边形或更高阶曲面片。
单元划分的密度通常与几何曲率、边界条件变化以及解的梯度有关。在应力集中、尖角或边界不连续处，往往需要更细的离散，以提高局部精度。
3.2 插值函数与形函数
在每个边界单元内，未知边界量通常通过插值函数进行近似。形函数用于将单元节点上的变量延拓到单元内部，从而建立连续或分段连续的近似表示。
低阶插值实现简单，适合快速计算；高阶插值则能更好地描述边界变量变化，特别是在曲率较大或解变化较快的区域。不过，高阶模型也会带来更复杂的积分与更高的实现成本。
3.3 系统方程的构建
完成边界离散后，边界积分方程会被转换为线性代数方程组。未知量通常包括边界上的场变量及其共轭量，系统矩阵由各边界单元之间的影响系数组成。
由于边界元法中的每个边界单元通常都会影响到其他单元，形成的矩阵一般是全耦合的。与局部耦合的体离散方法相比，这类系统在规模扩大时更容易出现计算负担。
3.3.1 矩阵元素计算
矩阵元素表示一个边界单元对另一个单元或观测点的影响强度。其计算通常涉及对基本解及形函数的积分。对于规则单元，可以采用解析积分、数值积分或半解析方法。
在远场单元之间，积分核通常较平滑，可用常规高斯积分处理；而在观测点接近源单元时，积分核会变得尖锐甚至奇异，需要更精细的数值策略。
3.3.2 奇异积分处理
奇异积分是边界元法中的关键技术问题。由于基本解在源点附近往往具有对数型或幂型奇异性，常规数值积分可能失效。为此，通常需要采用坐标变换、奇异性分离、主值积分或解析积分等方法。
对于近奇异积分，还常借助局部细分、自适应积分或专门的奇异处理公式，以提高精度和稳定性。奇异积分处理的好坏，往往直接影响整体结果质量。
3.4 数值求解与后处理
完成系统方程组后，可使用直接法或迭代法求解边界未知量。得到边界解后，再通过积分公式恢复域内任意点的响应，如位移、应力、势函数、温度或声压等。
后处理阶段通常包括场量可视化、误差评估和工程量提取。对于结构分析，还可能进一步计算应力集中系数、反力分布和边界通量；对于声学与电磁问题，则常关注散射截面、辐射功率和场强分布。
4 经典边界元模型
4.1 位势问题边界元法
位势问题通常对应拉普拉斯方程或泊松方程的某些简化情形。边界元法在这类问题中形式最为清晰，常用于静电势、稳态扩散和理想流动等模型。
在均匀介质中，位势场往往能直接由边界单层势或双层势构造，因此模型结构相对简单，也便于用来说明边界元法的基本原理。
4.2 弹性力学边界元法
弹性力学中的边界元法主要用于求解位移场、应力场及表面力分布。其基本解来源于弹性理论中的格林函数，能够描述点力作用下的位移响应。
这类方法适合处理应力集中、孔洞、裂纹和无限域弹性体等问题。由于边界上力与位移之间存在明确的共轭关系，边界元法在弹性接触和断裂分析中也具有较高实用价值。
4.3 热传导边界元法
热传导问题中，边界元法常用于稳态导热分析，特别是在材料热导率均匀且边界条件明确时。温度和热流之间的关系可通过热传导基本解写成边界积分形式。
在热设计中，它常用于研究散热器、电子器件外壳以及大尺度热扩散区域。对于无限域或半无限域中的温度场，边界元法往往比体网格方法更简洁。
4.4 声学边界元法
声学边界元法主要用于求解亥姆霍兹方程及相关波动问题，尤其适合辐射、散射和噪声传播分析。声场在外部流体域中的传播特性，使其与边界积分方法高度契合。
在车辆噪声、机械振动声辐射和建筑声学中，边界元法常被用来预测远场声压和声强分布。对于开放边界问题，它能自然地处理辐射条件，而不必显式截断无限远空间。
4.5 电磁场边界元法
电磁场边界元法通常面向高频散射、导体表面电流和开域辐射问题。它通过电场和磁场的边界积分表示，将复杂场域问题转化为边界上的未知量求解。
在天线设计、电磁兼容分析和目标散射研究中，这类方法有较高应用价值。不过，电磁问题往往方程耦合更强，数值实现也更容易受到频率、尺度和网格质量的影响。
5 算法与实现技术
5.1 线性与非线性问题处理
边界元法在经典形式中更适合线性问题。若材料行为、边界条件或控制方程出现非线性，通常需要采用迭代线性化、增量求解或与其他数值方法耦合。
非线性处理常见于接触、塑性、热物性变化和大变形等情形。此时，边界元法往往不再单独完成全部计算，而是与域离散技术协同工作。
5.2 时间相关问题
对于瞬态问题，边界元法可以在时域或频域中建立求解框架。时间依赖性会显著增加问题复杂度，因为场量不仅与空间位置有关，还与时间演化过程相关。
5.2.1 时域边界元法
时域边界元法直接在时间变量上离散，能较直观地描述瞬态响应。它适合冲击载荷、瞬态波传播和时变边界激励等问题。
这种方法的优点是物理意义清楚，但计算中往往要处理历史项，随着时间步增长，累积代价也会增加。因此其实现通常比稳态问题更复杂。
5.2.2 频域边界元法
频域边界元法通过对时间变量进行变换，把瞬态问题转化为频率域中的稳态型问题。它特别适用于周期激励、稳态振动和声学辐射分析。
在频域中，每个频点通常独立求解，便于并行处理。但若需要完整的时域响应，则还要再进行反变换或频率合成。
5.3 快速算法与加速技术
由于边界元法会生成稠密矩阵，直接求解在大规模问题中代价较高，因此需要借助快速算法降低时间和存储消耗。相关技术已经成为现代边界元计算的重要组成部分。
5.3.1 快速多极子方法
快速多极子方法通过把远距离单元的相互作用聚合处理，减少直接计算的次数，从而显著降低复杂度。它特别适合处理大规模散射和辐射问题。
该方法将局部作用与远场作用分离，对远场部分进行层次化近似，是边界元法实现高效计算的重要手段之一。
5.3.2 迭代求解策略
对于大型线性系统，迭代求解常比直接分解更经济。配合预条件技术，可以在较少内存占用下求得可接受精度的解。
常见策略包括共轭梯度类方法、GMRES 等通用迭代框架。其效果取决于矩阵性质、预条件器质量以及收敛判据设置。
5.4 稳定性与收敛性分析
边界元法的稳定性与收敛性取决于离散阶次、单元划分、积分精度和方程形式。若积分误差过大或边界处理不当，可能导致解振荡、数值病态或收敛缓慢。
在实际应用中，通常需要通过网格加密、误差估计与对照试验来检验结果可靠性。稳定且可收敛的模型，往往需要在精度和计算成本之间取得平衡。
6 工程应用
6.1 结构分析
在结构工程中，边界元法常用于分析受力构件、壳体和开域弹性体的响应。它能够较有效地预测边界位移、表面应力和局部集中现象。
当结构外部环境影响较明显时，如土-结构相互作用或无限域弹性问题，边界元法可减少人为截断边界带来的误差。
6.2 断裂与接触问题
边界元法在裂纹扩展、应力强度因子计算和接触边界分析中具有较高实用性。由于裂纹问题本质上往往与边界奇异性有关，因此特别适合采用边界积分形式研究。
接触问题中，边界条件通常会随接触状态变化而更新，数值上需要反复迭代。边界元法虽能提供较清晰的边界力学图像，但实现时需要仔细处理接触约束和非光滑边界。
6.3 噪声与振动分析
在噪声控制和振动工程中，边界元法常用于分析声源辐射、结构振动向空气中的传递以及外部声场分布。它能较自然地描述开放空间中的波传播。
对于车辆、机械设备和建筑构件的振动噪声预测，这一方法常与模态分析或有限元结果结合，以获得更完整的系统响应。
6.4 热管理与传热设计
在电子散热、设备热设计和热防护研究中，边界元法常用于稳态或准稳态热分析。它可在较少离散量的情况下估计复杂外形表面的温度分布。
对于外部冷却环境或半无限导热介质，边界元法尤其适合用来评估热扩散和边界换热效应，减少对大规模体网格的依赖。
6.5 电磁兼容与散射问题
在电磁兼容、电磁散射和导体外场分析中，边界元法可以直接描述场在开放空间中的传播与反射。对于导体表面电流、散射波与远场特性，边界积分框架较为合适。
这类应用中常关注屏蔽效果、散射截面以及外部辐射模式。边界元法因此成为高频电磁工程中的重要工具之一。
7 优势与局限
7.1 计算维度降低的优势
边界元法最突出的优点之一是维度降低。由于只离散边界，模型规模通常比体离散方法更小，尤其在三维问题中，这种优势更为明显。
在边界变化相对集中的问题里，减少的自由度可以显著降低建模和求解负担，使计算更聚焦于关键区域。
7.2 无限域问题处理能力
对于无限域和外部域问题，边界元法无需像体网格方法那样人为截断远场边界，因此更容易满足辐射条件和远场衰减规律。
这使它在声学、波动传播和开放边界弹性分析中具有独特地位，常能得到较自然的外场描述。
7.3 对基本解的依赖
边界元法的适用性在很大程度上依赖于基本解的可得性。若控制方程缺少方便使用的解析基本解，或基本解形式过于复杂，方法实现就会受到限制。
因此，该方法在均匀线性系统中更成熟，而在复杂介质和复杂方程中常面临理论与数值上的双重障碍。
7.4 非线性与非均匀介质中的困难
当材料性质随空间位置变化明显，或控制方程具有强非线性时，边界元法通常难以单独处理全部物理效应。此时，积分表示可能不再保持简单闭合，需要引入域积分或迭代策略。
这些问题会削弱边界元法“只算边界”的简洁性，使其优势不如在线性均匀问题中那样突出。
7.5 大规模计算的存储与效率问题
边界元法形成的系统矩阵通常是稠密的，这意味着存储需求和运算量都可能迅速增长。对于大规模三维模型，直接法往往难以承受。
因此，虽然边界元法在理论上具有高效性，但在工程大模型中仍需要依赖压缩存储、快速算法和并行计算来维持可用性。
8 发展历程
8.1 理论起源
边界元法的思想源于积分方程理论、势理论和格林函数方法。早期研究主要集中在解析解和边值问题的积分表示上，为后来的数值化奠定了基础。
随着数学物理和工程分析的发展，人们逐渐认识到边界积分形式在处理外场和无限域问题时的便利性，这推动了其从理论工具走向计算方法。
8.2 计算机辅助发展
计算机技术的引入，使边界积分方程能够被离散并转化为可求解的代数系统。由此，边界元法开始从手工推导的分析方法，演变为可用于工程实际的数值工具。
在这一阶段，单元技术、矩阵组装和数值积分方法逐步成熟，边界元法的应用范围也随之扩展到结构、声学和热学等多个领域。
8.3 现代高效边界元技术
随着问题规模不断扩大，传统边界元法在存储和计算上的不足愈发明显。为此，快速多极子、层次矩阵、预条件迭代等高效技术相继出现。
这些进展使边界元法能够处理更大规模的三维问题，也推动了其在高频波动、复杂散射和大范围场计算中的应用。
8.4 与其他方法的融合发展
现代边界元法越来越常与有限元法、域积分法和多尺度方法结合使用。通过耦合框架，可以兼顾局部非线性、复杂材料和开放边界的计算需求。
这种融合趋势使边界元法不再是独立的孤立工具，而是成为多种数值技术体系中的关键组成部分。
9 相关概念
9.1 边界积分方程
边界积分方程是将偏微分方程转化为边界上积分关系的方程形式，是边界元法的理论核心。它描述边界未知量与边界条件之间的耦合。
9.2 域积分法
域积分法是直接在求解区域内部建立积分表示或数值积分的计算方法。与边界元法相比，它通常更适合处理非均匀介质和体源项较强的问题。
9.3 有限元-边界元耦合法
有限元-边界元耦合法是一种将有限元法与边界元法结合的计算框架，常用于既包含复杂内部结构又涉及无限域外场的问题。它兼顾了体域建模与边界辐射处理的优势。
9.4 无网格方法
无网格方法是一类不依赖传统单元网格连接关系的数值技术，强调节点和核函数的近似表示。它与边界元法在降低网格依赖方面有相近目标，但实现路线和理论基础不同。