3. FEA应变输出基础
在结构应变测试与FEA对标工作中，最容易出现的一类误解是：把有限元后处理里读到的“应变”当成和应变片一样直接测得的物理量。实际上，有限元求解器真正解出的基本未知量是节点位移，应变和应力都是在此基础上通过求导、插值和本构关系“算出来”的派生量。理解这条计算链，才能判断后处理中读到的应变值是否可信、是否与应变片测点可比。
一、从位移到应变：求导带来的精度损失
有限元位移法的基本方程是 [K]{u} = {F}，求解得到节点位移向量 {u}。单元内部的位移场由形函数插值得到：
u(x) = Σ Ni(x) · ui （节点位移插值）
ε(x) = du/dx = Σ (dNi/dx) · ui （对形函数求导）
σ(x) = D · ε(x) （本构关系，D为弹性矩阵）关键点在于：应变是位移的一阶导数。位移解本身在节点上是精确的（在收敛意义下），但求导会放大误差。因此应变、应力的精度天然低于位移精度，这也是为什么对标时应优先核对位移或整体刚度，再核对局部应变。
另一个直接后果是：应变在单元内部通常不是常数，而是随位置变化的。以四节点四边形单元（Q4）为例，其位移场是双线性的，应变场是双线性的导数，在单元内呈线性变化；而三节点三角形单元（CST）位移线性、应变为常数，整个单元只有一个应变值。这就引出了单元类型对strain结果的核心影响。
二、单元类型如何影响应变结果
不同单元类型的应变表达能力差异很大，直接影响后处理读数的物理含义：
- CST（常应变三角形）：应变为常数，网格稍粗就会严重“过刚”，应变普遍偏小。仅适合粗略估算，不建议用于应变对标。
- Q4 / 六面体线性单元：应变在单元内线性变化，但存在剪切自锁和沙漏风险，弯曲主导区域应变偏小，需要足够密的网格。
- Q8 / 二次单元：位移二次、应变线性，能较好捕捉弯曲和应力集中，是应变对标的首选。
- 壳单元：需区分中面薄膜应变与上下表面弯曲应变，应变片贴在表面，对标时必须取表面层应变而非中面值。
- 实体单元：表面节点的应变需通过单元内插值外推到表面，外推方式（单元内插值 vs 节点平均）会带来差异。
此外，积分点与节点的区别至关重要。求解器先在高斯积分点计算应变（精度最高），再外推并平均到节点。因此同一位置，积分点值、单元节点值、全局节点平均值可能不同。对标时应明确后处理读取的是哪一种。
三、后处理中的应变提取方式
常见的应变提取方式有三类，用途不同：
- 云图（Contour）：默认显示节点平均应变，视觉平滑，但会削低峰值。用于整体趋势判断，不用于精确对标。
- 探针 / 路径提取（Probe / Path）：在指定坐标或沿路径读取积分点或节点值，可保留峰值，适合与应变片测点一一对应。
- 单元/积分点列表：直接输出每个积分点的应变分量，精度最高，适合做定量对标与误差统计。
对标时还需注意应变分量方向。应变片只测栅长方向的平均应变，因此FEA中应提取与片栅方向一致的正应变分量（如 εx），而不是等效应变或主应变。若测点为三向应变花，则需提取 εx、εy、γxy 三个分量，再按应变花公式换算主应变，与FEA同位置同分量比较。
四、一个可复现的示例参数
假设一块 200 mm × 20 mm × 5 mm 的钢制悬臂梁，左端固定，右端施加 50 N 横向力，材料 E = 206 GPa、ν = 0.3。在距固定端 30 mm 的上表面贴一片沿梁轴向的应变片。用二次六面体单元（C3D20）建模，网格在贴片区域加密至 1 mm。后处理时：
- 读取该点上表面的 εx（轴向正应变），而非中面值；
- 使用探针而非云图，避免节点平均削峰；
- 记录积分点值与节点值两组数据，观察差异。
若积分点 εx ≈ 320 με，节点平均 εx ≈ 305 με，而应变片实测约 315 με，则说明节点平均确实削低了峰值，积分点值更接近实测。这类差异正是对标中需要量化和解释的内容，而不是简单归因于“仿真不准”。