1 历史与发展
有限元法的形成并非源于单一学科，而是数学、力学与计算技术长期交汇的结果。它最初服务于复杂结构分析，随后逐步扩展到热学、流体、电磁学等多个领域，最终发展为通用型数值方法。随着工程问题规模增大、模型复杂度提高，有限元法也从手工推导和小规模计算，演变为高度依赖计算机的系统化分析框架。
1.1 早期数学与工程背景
在有限元法出现之前，工程计算主要依赖解析法、经验公式和近似理论。对于几何复杂、边界条件繁多或材料行为非线性的问题，传统方法往往难以给出可用解。与此同时，变分法、插值理论和数值分析等数学分支不断成熟，为后来离散化思想的建立提供了基础。
1.2 有限元法的形成
有限元法的核心观念是把连续问题拆分为若干局部问题，再通过整体组装恢复系统行为。这一思路在20世纪中叶逐渐明确，并在航空航天、土木和机械设计中得到实践验证。其成熟过程体现了工程需求对数学方法的反向推动。
1.2.1 结构力学中的萌芽
早期的结构分析中，研究者已尝试将复杂构件划分为较小部分进行近似计算。桁架、框架和薄板问题尤其适合这种处理方式，因为它们本身具有明显的离散连接特征。此类尝试虽未形成统一理论，却为后来的单元法提供了直观原型。
1.2.2 变分原理与弱形式的影响
变分原理使得力学问题能够通过能量最小化或驻值条件来表达，而弱形式则降低了对解的光滑性要求。两者共同促成了有限元法的数学基础，使复杂偏微分方程能够转化为更适合计算机处理的积分方程组。这一转化是有限元法广泛适用的关键原因之一。
1.3 计算机时代的推动
电子计算机的出现极大改变了有限元法的发展速度。矩阵组装、方程求解和网格管理等过程高度依赖计算能力，只有在计算机可高效执行之后，有限元法才真正具备大规模应用的条件。随后，算法优化、图形界面和自动网格生成工具进一步降低了使用门槛。
1.4 现代有限元法的发展趋势
现代有限元法正向更高精度、更强耦合和更大规模方向演进。自适应网格、并行计算、非线性材料模型、多尺度分析等技术不断完善，使其能够处理更复杂的工程场景。同时，等几何分析、数据驱动方法和多物理场集成也在扩展其应用边界。
2 基本原理
有限元法的基本思想是将一个连续区域划分为有限数量的子区域，并在每个子区域内用简单函数描述未知量的变化。通过在单元层面建立局部方程，再将这些方程按节点关系组装，可得到整个系统的近似解。该方法兼具数学严谨性与工程可操作性，因此成为求解复杂边值问题的常用工具。
2.1 连续体离散化思想
连续体离散化是有限元法最核心的步骤。原本在空间上连续分布的位移、温度或电势等物理量，被转换为有限个节点上的离散未知量。随着单元数量增加、尺寸减小，离散模型通常会逐步逼近原始连续模型。
2.2 单元划分与节点
在有限元模型中，单元是基本计算单元，节点则是单元之间交换信息的位置。节点上通常存储位移、温度等自由度，而单元内部的场变量由节点值通过插值得到。单元划分质量直接影响计算精度、稳定性和效率。
2.2.1 网格的概念
网格是单元和节点在空间中的整体排列方式。它不仅决定几何形状如何被近似，也影响误差分布、收敛速度和求解难度。合适的网格通常需要兼顾单元尺寸、形状质量以及局部细化需求。
2.2.2 单元类型
单元类型按照维度和几何形式可分为一维、二维和三维等多种类别。不同单元适合不同问题：细长构件常用梁单元，薄壁结构常用壳单元，实体材料分析则多用三维单元。选型不当可能导致计算效率下降或结果失真。
2.3 形函数与插值
形函数用于描述单元内部物理量如何由节点值插值得到。它们通常具有局部支持性和一定的连续性要求，是连接离散节点与连续场变量的桥梁。形函数的阶次越高，通常可获得更好的近似能力，但计算代价也会增加。
2.4 弱形式与加权残值法
弱形式通过对原方程进行积分处理，降低了对解的导数阶数要求。加权残值法则是在近似解代入方程后，要求残差在某种加权意义下为零，从而获得离散方程。二者在有限元法中紧密相关，构成了从微分方程到代数方程的重要过渡。
2.4.1 变分法基础
变分法关注函数的极值问题，常用于建立物理系统的能量表达。对许多力学系统而言，平衡状态对应某种泛函的驻值，因此可通过变分过程导出控制方程。有限元法在处理结构问题时，常直接借助这一理论基础。
2.4.2 积分形式的建立
将控制方程乘以权函数并在区域内积分，可得到等价的积分表达。此过程把局部微分条件转化为全局平均意义上的约束，允许近似解在单元内部存在一定误差。积分形式也更适合数值积分与矩阵化处理。
2.5 边界条件处理
边界条件决定了求解问题是否唯一以及解的物理意义。有限元法通常需要同时处理几何边界、约束边界和载荷边界，并将它们准确映射到离散模型中。边界条件设置不合理，往往会导致整体结果偏离实际。
2.5.1 位移边界条件
位移边界条件也称本质边界条件，通常直接规定某些节点的位移、温度或电势值。由于这类条件直接约束未知量本身，通常在组装方程时通过消元或罚函数等方式施加。
2.5.2 力边界条件
力边界条件属于自然边界条件，表现为表面力、热流或电流密度等外部作用。它们通常以边界积分的形式进入系统方程，反映外界对计算域的输入或交换。
3 数学基础
有限元法建立在偏微分方程、线性代数、数值积分和误差理论等多个数学分支之上。它的有效性不仅取决于模型本身，还依赖于离散化、求解和误差控制的数学性质。理解这些基础，有助于判断算法适用范围及结果可信度。
3.1 偏微分方程与边值问题
许多工程现象都可归结为偏微分方程所描述的边值问题，例如弹性平衡、热扩散和势场分布。此类问题通常具有复杂边界和多尺度特征，解析求解困难，而有限元法则通过离散近似为其提供了实用解法。
3.2 线性代数基础
有限元离散后通常得到规模庞大的线性方程组或非线性方程组。矩阵的结构、求解效率和数值稳定性，决定了计算可行性。线性代数因此成为有限元法不可或缺的基础支撑。
3.2.1 稀疏矩阵
由于每个单元只与相邻节点发生耦合，有限元方程矩阵通常具有稀疏特征。这种结构既减少了存储需求，也为专门的求解算法提供了优化空间。稀疏矩阵技术是大规模有限元计算的重要组成部分。
3.2.2 特征值问题
在振动、屈曲和稳定性分析中，有限元模型常转化为特征值问题。通过求解特征值和特征向量，可以获得固有频率、失稳模态等关键信息。这类问题对矩阵精度和算法鲁棒性要求较高。
3.3 数值积分
许多单元矩阵和载荷向量都需要通过积分计算获得。由于解析积分常不可行，实际中常采用高斯积分等数值积分方法。积分精度与单元阶次、几何复杂度和材料非线性程度密切相关。
3.4 收敛性与误差分析
收敛性描述了当网格逐渐细化时，数值解逼近真实解的趋势。误差分析则试图区分离散误差、截断误差和建模误差等来源。二者共同用于评估结果的可靠程度。
3.4.1 网格收敛
网格收敛研究的是网格尺寸减小时，解是否稳定接近某一极限值。一般而言，更细的网格可提高精度，但在高梯度区域或奇异点附近，仍需特别处理才能获得可靠收敛。
3.4.2 近似误差来源
误差主要来自单元插值能力有限、网格质量不足、边界条件近似以及数值积分误差等方面。若材料模型或几何简化过强，还会引入额外的模型误差。对这些来源进行辨识，是改进仿真的基础。
4 有限元建模流程
有限元分析通常遵循从问题抽象到结果解释的完整流程。每一步都可能影响最终结论，因此建模并不只是“划网格并求解”，而是一个连续的工程判断过程。合理的流程设计能够显著提升结果质量与计算效率。
4.1 问题定义与理想化
建模首先要明确分析目的，例如强度评估、温度分布还是固有频率预测。随后需要对真实对象进行理想化处理，去除对结果影响较小的细节，同时保留关键特征。理想化程度过高或过低，都会影响模型的有效性。
4.2 几何建模
几何建模用于建立待分析对象的空间形状。实际工程中，模型常需经过简化、修补和特征删除，以适配网格生成和求解需求。几何质量直接关系到后续单元划分的顺畅程度。
4.3 网格划分
网格划分是将几何模型转化为离散计算域的关键步骤。单元大小、形状和分布应结合分析目标进行设计，重要区域通常需要加密，而低敏感区域可适当粗化，以平衡精度和效率。
4.3.1 结构化网格
结构化网格具有规则的拓扑关系，节点编号和单元连接清晰，便于存储和高效计算。它常用于几何规则、边界平滑的问题，但在复杂外形上生成难度较大。
4.3.2 非结构化网格
非结构化网格对复杂几何更具适应性，常用于工程实际中的不规则模型。其单元连接更灵活，但数据管理和局部质量控制通常更为复杂。
4.4 材料参数设置
材料参数决定模型对载荷的响应方式，包括弹性模量、泊松比、导热系数、密度等。对于各向异性、温度相关或非线性材料，还需定义更复杂的本构关系。参数选取是否准确，往往直接决定仿真结果的可信度。
4.5 求解器选择
求解器决定方程组如何被求解。对于小规模线性问题，直接法较为稳妥；对于超大规模系统或非线性问题，则常需迭代法或专用算法。选择时需要考虑模型规模、非线性强度和计算资源。
4.6 后处理与结果解释
后处理包括应力云图、变形图、温度场和误差指标等结果展示。更重要的是对结果进行物理解释，判断峰值是否合理、局部异常是否源自网格或边界设定，并结合工程经验形成结论。
5 单元类型
单元类型决定了有限元模型中基本近似单元的几何表达方式和自由度配置。不同单元在精度、计算量和适用场景方面各有侧重，因此需要根据问题特点进行选择。
5.1 一维单元
一维单元通常用于描述细长结构，如杆、梁和索等。其计算模型较为简洁，适合轴向、弯曲或扭转等主导行为明显的对象。尽管形式简单，但在许多工程问题中非常实用。
5.2 二维单元
二维单元主要用于平面应力、平面应变和薄板类问题。它们能够在较低计算成本下反映平面内的力学或热学行为，因此常见于结构局部分析和截面研究。
5.2.1 三角形单元
三角形单元具有良好的几何适应性，适合复杂边界和不规则区域。其网格生成相对容易，但低阶三角单元的精度往往不如较高阶四边形单元。
5.2.2 四边形单元
四边形单元在规则区域中通常具有较好的数值性能和更高的插值效率。对于弯曲较少或可规则分区的几何，四边形网格往往能提供更平滑的结果。
5.3 三维单元
三维单元用于实体结构分析，能够直接描述真实的空间应力、应变和传热过程。它们是最通用也最耗费计算资源的一类单元，适用于复杂载荷和空间耦合问题。
5.3.1 四面体单元
四面体单元适合复杂三维几何，自动网格生成能力强，使用广泛。对于初步分析和复杂外形模型，它往往是最容易实施的选择，但在精度与单元数量之间需要权衡。
5.3.2 六面体单元
六面体单元通常在数值精度和计算效率上表现较好，尤其适合规则实体或可映射分区的模型。它对网格质量要求较高，但在高质量网格下常能获得优良结果。
5.4 特殊单元
特殊单元针对特定结构特点建立，如薄壁、细长或层状构件。它们通过简化几何维度或强化某类自由度表达，提高了某些问题的计算效率。
5.4.1 壳单元
壳单元适用于厚度远小于其他尺寸的结构，如金属板壳和薄壁容器。它能在保留面内与弯曲行为的同时减少实体单元数量。
5.4.2 板单元
板单元主要用于薄板弯曲问题，重点描述法向挠度与弯曲应力。与壳单元相比，它对面内效应的表达通常更简化。
5.4.3 梁单元
梁单元适合线状结构的受力分析，可描述轴力、剪力、弯矩和扭矩等效应。它在结构设计和整体刚度分析中常被广泛采用。
6 常见应用领域
有限元法的应用范围十分广泛，几乎覆盖所有涉及连续介质和复杂边界的工程与科学问题。不同领域虽物理量不同，但计算框架高度相似，均依赖离散化、组装和求解过程。
6.1 结构分析
结构分析是有限元法最经典的应用方向之一，主要研究构件或结构在外载作用下的应力、应变和变形。其结果可用于强度校核、刚度评估和安全裕度判断。
6.1.1 静力分析
静力分析关注系统在静态载荷下的平衡状态，如重力、压力和固定载荷等。通过有限元法可获得位移场、应力分布及关键危险部位。
6.1.2 动力分析
动力分析考虑时间变化载荷或惯性效应，常用于振动、冲击和疲劳相关研究。此类问题通常更复杂，需要处理质量、阻尼和时间积分等因素。
6.2 热分析
热分析用于研究温度场、热流和热应力等问题。它可应用于电子器件散热、热处理过程以及高温构件设计等场景。若材料性质随温度变化，计算过程通常会更复杂。
6.3 流体分析
在流体分析中，有限元法可用于求解速度场、压力场和流动稳定性等问题。虽然流体计算还常与其他离散方法并行发展，但有限元法在复杂边界和耦合问题中具有明显优势。
6.4 电磁场分析
电磁场分析主要用于求解电势、磁场强度和涡流分布等。有限元法能够处理不规则几何、复杂材料和边界条件，因此在电机、传感器和电磁装置设计中十分常见。
6.5 多物理场耦合
多物理场耦合指多个物理过程相互影响，如热-结构耦合、电-热耦合或流固耦合。有限元法可通过统一离散框架描述这些相互作用，是复杂系统仿真的重要基础。
7 算法与求解技术
有限元离散后形成的方程组往往规模巨大且结构复杂，因此需要高效算法支持。求解技术不仅影响速度，也决定模型能否在合理时间内完成计算。
7.1 直接法
直接法通过矩阵分解等方式一步得到方程解，通常稳定可靠，适合中小规模问题。其缺点是内存占用较高，在超大规模模型中可能受限。
7.2 迭代法
迭代法通过逐步逼近求解结果，适合稀疏且规模很大的线性系统。它通常比直接法更节省存储，但对初值、预条件和收敛条件较为敏感。
7.2.1 共轭梯度法
共轭梯度法常用于对称正定矩阵系统，具有较好的效率和存储性能。它在结构分析等问题中应用广泛，尤其适合稀疏大系统。
7.2.2 Krylov子空间方法
Krylov子空间方法是一类常见迭代技术，适用于多种非对称或大型稀疏系统。其性能通常与预条件器密切相关，在复杂工程计算中非常重要。
7.3 非线性求解
非线性问题广泛存在于材料塑性、大变形、接触和几何非线性分析中。此类问题通常不能一次求解，需要借助迭代策略逐步逼近平衡状态。
7.3.1 牛顿迭代法
牛顿迭代法通过线性化非线性方程并逐步修正未知量，具有较快的局部收敛速度。它是有限元非线性求解中最常用的方法之一。
7.3.2 路径跟踪方法
路径跟踪方法用于处理载荷变化过程中可能出现的分叉、失稳或极限点问题。它能够帮助模型跨越复杂响应区间，获得更完整的平衡路径。
7.4 瞬态问题求解
瞬态问题涉及时间演化过程，如热扩散、动态响应和时变场问题。求解时需要同时考虑空间离散与时间离散，时间步长的选择会明显影响精度与稳定性。
7.5 并行计算与加速技术
随着模型规模持续增大，并行计算已成为有限元求解的重要手段。多核处理、分布式计算和图形处理器加速能够显著缩短分析时间，也推动了更大规模、更复杂模型的应用。
8 误差、稳定性与验证
有限元结果的可信度取决于误差控制、数值稳定性以及是否经过充分验证。实际工程中，结果不仅要“算得出”，还要“算得对”，这需要系统的检验机制。
8.1 模型误差
模型误差来自对真实系统的简化与假设，例如忽略某些细节、采用理想材料模型或简化边界条件。它并不属于数值计算本身，但往往对最终结论影响最大。
8.2 数值误差
数值误差主要包括离散误差、舍入误差和求解误差。它们通常与网格大小、算法精度和计算机浮点运算有关。通过改善单元质量和求解设置，可在一定程度上减小这类误差。
8.3 稳定性分析
稳定性分析用于判断数值过程是否会在扰动下保持可靠。若算法对微小变化过于敏感，结果可能出现发散或异常振荡。稳定性良好的方法更适合工程应用。
8.4 验证与确认
验证与确认是确保模型可信的重要步骤。验证关注“数值实现是否正确”，确认则关注“模型是否真实反映物理现象”。二者结合，才能提高有限元分析的可用性。
8.4.1 基准算例
基准算例是具有已知解析解或高可信参考结果的测试问题。通过与这些标准问题比较，可以检查程序实现、单元性能和求解精度。
8.4.2 实验对比
实验对比是将仿真结果与实际测量数据进行比较。若二者在可接受范围内一致，说明模型参数、边界条件和数值设置较为合理。
9 软件与工具
有限元法的应用离不开专门软件支持。现代软件通常集成建模、网格生成、求解和后处理功能，并提供脚本接口以适应不同复杂度的工程任务。
9.1 商业软件
商业有限元软件通常具有成熟的界面、稳定的求解器和完善的技术支持，适合工程项目和工业流程。它们在通用功能、可靠性和工作效率方面优势明显。
9.2 开源软件
开源软件为研究和教学提供了灵活平台，用户可查看源码、修改算法并进行二次开发。它们常用于学术验证、方法研究和定制化仿真。
9.3 前处理与后处理工具
前处理工具负责几何导入、网格划分和参数设置，后处理工具则用于查看结果和生成报告。二者构成有限元工作流的重要两端，直接影响使用体验和分析效率。
9.4 脚本化与二次开发
脚本化使有限元流程更加自动化，适合重复任务、参数扫描和定制算法。二次开发则能将有限元平台与特定工程需求结合，提高灵活性。
9.4.1 自动建模
自动建模可减少重复操作，适合系列化产品或标准化结构的批量分析。它通常依赖脚本、模板和参数化几何。
9.4.2 批量计算
批量计算用于一次运行多个工况或参数组合，常见于优化设计、灵敏度分析和统计研究。它对计算管理和资源调度能力要求较高。
10 相关方法与扩展
有限元法并非唯一的数值离散方法，其他方法在某些场景中具有互补优势。了解这些方法有助于在不同问题类型中选择更合适的求解框架。
10.1 有限差分法
有限差分法通过网格点上的差分近似导数，形式简单、实现直接，适合规则区域和较标准的偏微分方程问题。与有限元法相比，它在复杂几何适应性上通常较弱。
10.2 有限体积法
有限体积法强调守恒量在控制体上的平衡，常用于流体力学和传热问题。它在守恒性方面表现突出，因此与有限元法在不同应用领域中各有优势。
10.3 边界元法
边界元法只对边界进行离散，适合无限域或半无限域问题。由于维数降低，它在特定场景中计算效率较高，但对适用方程类型有一定限制。
10.4 谱元法
谱元法结合了有限元的分区思想和谱方法的高阶逼近能力，通常具有较高精度。它适用于对解精度要求极高的科学计算问题。
10.5 等几何分析
等几何分析将几何建模与数值分析统一在同一套函数表示下，常与CAD模型直接衔接。它有助于减少几何近似误差，并提高设计到分析的衔接效率。
11 典型工程案例
有限元法在工程实践中通常用于评估结构安全、优化设计和预测失效风险。典型案例覆盖从大型土木结构到微观材料尺度，体现了其跨尺度适应能力。
11.1 桥梁与建筑结构
在桥梁和建筑结构分析中，有限元法可用于评估承载能力、变形控制和局部应力集中。对于复杂连接、不同材料组合和非规则几何，它尤其适合进行整体性能预测。
11.2 机械零部件强度分析
机械零部件常面临交变载荷、接触和疲劳等问题，有限元法可帮助识别危险区域并优化结构尺寸。齿轮、轴类和壳体等部件都常见此类分析。
11.3 复合材料仿真
复合材料具有各向异性和层合结构特点，有限元法可用于研究层间应力、分层风险和整体刚度。其分析过程往往比均质材料更复杂，对材料模型要求更高。
11.4 高温与热应力问题
在高温环境中，温度梯度会引起显著热应力，进而影响结构寿命。有限元法可同时处理温度场与应力场，帮助评估热变形、热疲劳及相关失效风险。
11.5 微结构与材料尺度分析
在微结构层面，有限元法可用于研究颗粒分布、孔隙结构和局部应力集中对宏观性能的影响。随着多尺度建模方法的发展，它在材料设计与性能预测中的作用日益突出。