1 基本概念
1.1 定义与作用
形函数是有限元分析中用来描述单元内部物理量分布的插值函数。它通过节点上的自由度，近似重建单元内任意位置的位移、应变、温度或其他场变量。
在工程计算中，形函数不仅承担“连接节点与单元内部”的作用，还用于建立位移场、推导单元刚度矩阵，并为后续的应力应变计算提供基础。
1.2 形函数的数学表达
一般情况下，单元内某个场变量可表示为各节点自由度与对应形函数的组合，即通过节点值加权得到连续的近似场。其常见形式为节点自由度乘以形函数后求和。
在局部坐标系中，形函数通常写成自然坐标的函数。对于不同单元类型，这些函数可以是线性的、多项式的、分片定义的，也可以采用更高阶插值形式。
1.3 与有限元法的关系
有限元法的核心思想之一，就是把连续体离散为若干单元，再用形函数在单元内部进行插值。形函数决定了单元如何表达未知量，因此是有限元离散的基础工具。
从求解流程看，形函数影响位移模式、单元刚度、载荷等效方式以及误差分布。可以说，没有形函数，就无法完成标准意义上的有限元建模与计算。
1.4 形函数的基本性质
1.4.1 节点插值性
节点插值性是指形函数在本节点取值为1，在其他节点取值为0。这样可保证节点自由度能够被准确传递到单元内部。
这一性质使得有限元解能够直接与节点未知量对应，便于组装整体方程，也便于后处理中的结果解释。
1.4.2 完备性
完备性是指形函数所构成的插值空间能够再现某些基本函数形式，通常至少应能表示常数项，较高阶单元还应能表示线性项及更复杂多项式项。
完备性不足会导致单元即使在网格细化后仍难以正确逼近真实解，从而影响计算精度和收敛表现。
1.4.3 连续性
连续性描述单元之间场变量及其导数的衔接情况。不同问题对连续性要求不同，例如位移型问题通常要求位移连续，而某些高阶理论还可能要求更高阶导数连续。
形函数的连续性等级直接决定单元之间能否平滑拼接，也影响模型是否适用于弯曲、板壳等对导数敏感的力学问题。
1.4.4 局部支撑性
局部支撑性是指某个形函数通常只在与其相关的少数单元内非零，而在远离节点的区域取零。
这种性质使有限元矩阵具有稀疏结构，显著提高了存储与求解效率，也便于采用大规模数值计算方法。
2 形函数的类型
2.1 一维形函数
一维形函数主要用于杆、梁等一维理想化结构。它们通常定义在单元轴线方向，并以节点坐标作为插值依据。
这类形函数形式较为直观，常见于线性杆单元、梁单元以及某些传热和扩散问题的离散表达。
2.1.1 线性形函数
线性形函数以节点坐标为基础，单元内部变量随位置作一次变化。其优点是表达简单、计算效率高。
由于阶次较低，它对曲线变化较强的场量逼近能力有限，但在粗略分析和初步设计中仍较常见。
2.1.2 二次形函数
二次形函数在单元中引入中间节点或附加控制点，可表示更平滑的曲线变化。与线性形函数相比，它通常能提供更好的精度。
在弯曲或梯度变化明显的问题中，二次插值往往更适合工程应用，也更容易获得较好的收敛速度。
2.1.3 高次形函数
高次形函数采用更高阶多项式进行插值，能够更准确地拟合复杂分布。随着阶次提高，单元内部变化的表达能力增强。
不过，高次单元对数值积分、网格质量和计算稳定性的要求也更高，因此在实际使用中需要平衡精度与成本。
2.2 二维形函数
二维形函数用于平面问题、板壳分析以及曲面单元描述。它们通常在二维自然坐标空间中建立，并根据节点布局定义插值关系。
其形式受到单元几何形状影响较大，常见类型包括三角形单元、四边形单元以及等参单元。
2.2.1 三角形单元形函数
三角形单元形函数常用于不规则区域的网格划分，适合复杂边界和局部加密。它的几何适应性较强，建模方便。
根据阶次不同，三角形单元可采用线性或更高阶插值。高阶三角形单元通常可改善应力梯度的表达能力。
2.2.2 四边形单元形函数
四边形单元形函数多用于规则或半规则区域，能够在较少单元数下取得较好的逼近效果。
这类单元常与自然坐标系配合使用，通过双线性或高阶插值描述面内变化，在工程结构中应用十分广泛。
2.2.3 等参单元形函数
等参单元是指几何形状与场变量都使用同一组形函数描述。这样既能表达单元边界形状，也能插值内部未知量。
这种统一处理方式便于建立复杂几何模型，并使坐标映射与位移插值保持一致，成为现代有限元的重要基础形式。
2.3 三维形函数
三维形函数用于实体结构分析，可描述体单元内部在三个方向上的场变量变化。其计算通常比一维和二维情形更复杂。
常见三维单元包括四面体、六面体，以及用于过渡连接的棱柱和金字塔单元。
2.3.1 四面体单元形函数
四面体单元便于处理复杂三维几何，尤其适合自动网格划分。其节点布局简单，适应性较强。
在实际应用中，四面体单元常用于不规则体域的离散，但若阶次较低，有时在精度与刚度方面不如六面体单元理想。
2.3.2 六面体单元形函数
六面体单元通常具有较好的数值性能和变形表达能力，尤其在规则网格中表现稳定。
其形函数多采用多维多项式插值，能够较有效地描述实体内部的应力与位移变化，因而在精细分析中常被优先考虑。
2.3.3 棱柱与金字塔单元形函数
棱柱与金字塔单元常用于不同网格类型之间的过渡区域。它们在几何适应性和单元连接上具有特殊作用。
这类单元的形函数构造相对复杂，但在混合网格中非常实用，能够提升整体网格的连续性与可生成性。
3 形函数的构造方法
3.1 多项式插值法
多项式插值法是构造形函数的基本思路之一，即用某一阶多项式逼近单元内的未知场。通过确定系数，可满足节点插值条件。
该方法理论清晰，适合从数学上推导单元公式，也便于扩展到不同阶次和不同维度。
3.2 拉格朗日插值法
拉格朗日插值法通过构造满足节点克罗内克性质的一组插值多项式，直接形成形函数表达式。
它是有限元中最常见的构造方式之一，优点是形式规范、实现简洁，特别适合位移型单元。
3.3 赫尔米特插值法
赫尔米特插值法不仅考虑节点处的函数值，还常同时考虑导数值，从而提高插值的光滑性。
这种方法适用于对转角、斜率或高阶连续性有要求的问题，在梁、板等弯曲分析中较有价值。
3.4 分片插值与样条方法
分片插值将整体曲线或曲面划分为若干区段，在每一段内分别构造插值函数，再通过连接条件保证整体光滑。
样条方法通常能够提供较好的连续性和更平顺的几何描述，在需要较高光滑度的场合有一定优势。
4 在不同结构单元中的应用
4.1 杆单元中的形函数
杆单元主要描述轴向拉压变形，形函数通常用于插值轴向位移。由于问题维度较低，其形式相对简单。
在杆系结构中，形函数决定了轴力与轴向位移之间的关系，也影响单元刚度矩阵的构造方式。
4.2 梁单元中的形函数
梁单元除位移外，往往还涉及转角或斜率自由度，因此形函数不仅要描述平移，还要协调弯曲变形。
梁单元中的插值要求通常高于杆单元，特别是在弯曲分析中，需要兼顾位移连续与转角连续。
4.2.1 位移型梁单元
位移型梁单元以节点位移为主要未知量，通过形函数重构梁轴线的挠度曲线。
这种单元形式便于理解和实现，是工程分析中常见的梁模型基础。
4.2.2 转角自由度与弯曲插值
当梁单元包含转角自由度时，形函数需要同时控制位移和斜率，使弯曲曲线具有合适的光滑性。
这类插值通常比纯位移插值更复杂，但能更准确地表达梁在弯矩作用下的变形特征。
4.3 板壳单元中的形函数
板壳单元用于描述薄板和壳体结构，其形函数通常要兼顾面内与面外变形，有时还要处理转角耦合问题。
由于板壳理论对连续性要求较高，相关形函数设计往往更注重光滑过渡和弯曲性能。
4.4 实体单元中的形函数
实体单元用于三维连续体分析，形函数在三个方向上同时插值位移分量。
这类单元可以处理复杂应力状态，但对网格质量和积分精度较为敏感，因此形函数选择直接关系到整体分析效果。
5 形函数的计算与实现
5.1 自然坐标系与局部坐标系
在有限元中，单元通常在自然坐标系或局部坐标系下定义形函数。这样可以统一不同几何单元的表达方式。
采用标准化坐标后，单元内部的数学形式更规整，便于公式推导和程序实现。
5.2 雅可比矩阵与坐标映射
雅可比矩阵用于描述自然坐标与实际坐标之间的映射关系。它在形函数导数转换中起着关键作用。
通过坐标映射，可以把规则参考单元上的插值结果转化到实际单元几何中，从而处理不规则形状与变形网格。
5.3 导数计算与应变矩阵
形函数的一阶导数常用于建立应变位移关系，并进一步形成应变矩阵。
在结构力学问题中，导数计算是否准确，直接影响应力结果、刚度矩阵以及后续的误差评估。
5.4 数值积分中的应用
形函数通常与数值积分配合使用，以计算单元刚度、质量矩阵和等效载荷。积分方案的选择与形函数阶次密切相关。
5.4.1 高斯积分点设置
高斯积分点的数量和位置通常根据单元类型及插值阶次确定。积分点过少可能导致精度不足，过多则增加计算量。
合理设置高斯点有助于在效率与准确性之间取得平衡。
5.4.2 积分精度与单元阶次
单元阶次越高，积分所需精度通常也越高。若积分阶次不足，可能引起刚度误差或数值假象。
因此，在实际实现中，形函数阶次、积分阶次和单元公式常需统一考虑。
6 影响形函数性能的因素
6.1 单元阶次选择
单元阶次决定了形函数的表达能力。低阶单元计算简洁，高阶单元通常精度更好。
具体选择时，应结合问题的场变量变化规律、计算资源以及所需结果精度进行权衡。
6.2 单元形状与畸变
单元几何形状会影响形函数的实际表现。若单元严重畸变，插值质量和数值稳定性可能下降。
因此，良好的网格质量是发挥形函数性能的重要前提，尤其在复杂结构区域更为明显。
6.3 网格划分密度
网格越密，通常越有利于捕捉局部变化，但同时也会增加自由度数量和计算成本。
形函数本身的阶次与网格密度往往需要联合考虑，才能获得合适的逼近效果。
6.4 计算稳定性与收敛性
形函数的设计会影响方程组的条件性质、数值稳定性和收敛速度。某些不合理的插值形式可能导致锁死、振荡或误差放大。
在工程计算中，稳定性与收敛性往往与单元选择、积分策略及网格质量共同相关。
7 相关问题与扩展
7.1 兼容性条件
兼容性条件要求相邻单元之间的场变量能够协调衔接，避免出现不连续或不合理的位移跳变。
这是有限元单元能够拼接成整体模型的基本要求，也是形函数设计的重要约束之一。
7.2 锁死现象与改进
锁死现象是指某些单元在特定受力条件下表现出异常刚硬，导致变形被过度抑制。
常见改进思路包括优化单元插值、调整积分方式或采用更合适的单元公式，以减轻不合理约束带来的影响。
7.3 超收敛与误差估计
在部分位置上，有限元结果可能比整体平均精度更高，这种现象常被称为超收敛。
误差估计则用于判断计算结果的可靠性，并为网格加密或阶次提升提供依据。形函数特性会影响超收敛区域的出现与误差分布。
7.4 p型、h型与hp型改进
p型改进是提高单元阶次，h型改进是细化网格，hp型改进则结合两者优点。
不同改进策略对应不同的形函数配置方式，常用于提升复杂问题中的计算精度与效率。
8 常见误解与术语辨析
8.1 形函数与基函数
形函数可以看作有限元中的一种基函数，但并非所有基函数都属于形函数。形函数强调节点插值性质，而基函数概念更为宽泛。
在有限元语境中，二者常被混用，但严格来说，形函数是与节点自由度直接对应的插值基。
8.2 形函数与插值函数
形函数本质上属于插值函数的一种，但其使用场景更明确，主要服务于有限元离散。
插值函数是更通用的数学概念，涵盖范围远大于有限元中的单元插值表达。
8.3 形函数与试函数
试函数是变分法或弱形式中用于近似未知解的函数，形函数则是构造试函数的具体工具之一。
在有限元中，试函数往往由形函数和节点自由度组合而成，因此两者关系密切但并不等同。
8.4 工程软件中的形函数设置
工程软件通常会根据单元类型自动选取默认形函数，用户在建模时更多是通过选择单元阶次、类型和网格方式来间接影响它。
某些软件还提供高阶单元、减缩积分或特殊插值选项，便于在精度、稳定性和计算量之间进行调整。