有限元分析是什么？从网格、载荷到结果验证

有限元分析（FEA，Finite Element Analysis）是把连续物体或物理场离散成有限个单元，结合材料、载荷和边界条件，求得位移、温度等量的近似分布，再解释其工程意义的分析方法。彩色云图只是结果的一种显示方式，分析的核心是建立什么模型、求解什么量，以及怎样确认这个结果值得使用。

在制造设计中，它可以帮助理解支架受力后的变形、零件内部的应力分布、温度变化或振动特性。有限元分析通常属于工业软件中的研发设计与工程仿真工作；软件能够完成大量计算，但不会替分析者判断实际零件如何安装、承受什么载荷。

有限元分析、有限元法与CAE是什么关系

有限元法（FEM，Finite Element Method）指底层数值方法；有限元分析FEA则指运用这种方法开展建模、求解与结果解释的工作。CAE是Computer-Aided Engineering，即计算机辅助工程，覆盖的工程分析活动更广，不应把所有CAE仿真都称为有限元分析。

COMSOL的有限元方法说明从物理方程出发解释了它的作用：复杂几何和边界条件下，往往难以得到解析表达式，于是用离散模型求近似解。结构问题可求位移并进一步得到应变、应力；传热问题可求温度场和热流。它并不是只会做“受力云图”的工具。

一个有意义的分析应先提出问题。例如，“这个支架会不会太软”需要关注特定位置的位移；“是否存在局部屈服风险”需要选择适用的应力与材料判断条件；“是否接近共振”则涉及固有频率和实际激励。问题不同，需要的分析类型、输入与验收指标也不同，不能从同一张静力云图同时推导所有结论。

单元、节点和方程怎样连接起来

“有限元”中的“元”是单元。它可以是表示细长构件的线单元、描述薄壁的壳单元，也可以是填充体积的实体单元。节点承载模型中的未知量；在结构分析里，这些未知量通常包括位移，某些单元还包括转动。节点数不等于自由度数，具体取决于单元形式与求解的物理量。

在常见的节点型有限元中，单元内部的场用形函数从节点值近似表示。例如，两个节点之间的轴向位移可以先用线性变化描述。相邻单元共享节点时，也共享相应未知量，因此单元之间必须协调，而不是各自随意移动、最后拼成一张图片。

在线性静力结构分析中，组装后的关系常写成K u＝f：K是整体刚度矩阵，u是节点位移等未知量组成的向量，f是等效节点载荷向量。每个单元贡献自己的关系，再按节点连接汇总到整体方程。施加已知位移条件后，求解其余位移；应变来自位移的空间变化，应力再由材料关系得到。Jönköping University的有限元教学章节给出了这一组装与求解过程。

这也解释了为什么位移看起来平滑，不代表应力已经足够准确：应力计算涉及位移的变化率，对局部网格与处理方式可能更敏感。热传导等其他物理问题也有各自的方程与未知量，不能把结构刚度方程无条件套到全部分析中。

计算前必须定义哪些条件

CAD描述形状，有限元模型还必须描述材料、连接和工况。可以按下面的顺序检查输入是否足以回答问题。

输入	需要说明的内容	遗漏或简化可能改变什么
分析目标与几何	关注位置、工况、保留的孔和圆角、梁壳实体理想化	删掉关键圆角，可能改变局部应力；用过于简化的模型可能遗漏目标现象
材料	弹性、塑性或其他适用模型，以及对应温度和状态的参数	刚度、变形和失效判断；弹性模量不能替代屈服强度
载荷	大小、方向、施加区域、时间变化；力、压力或体载荷的含义	合力、力矩和载荷进入结构的路径
约束与连接	哪些方向限制移动、支承刚度、接触或连接方式	反力如何产生、各零件怎样传力、是否发生相对滑动

把螺栓安装面全部设为完全固定，相当于作出很强的支承假设；实际连接如果能够转动、滑移或变形，结果就可能不同。对常规静力模型，约束不足会留下刚体运动，使问题无法按预期求解；约束过多则可能把真实结构算得过硬。目的应是表达实际连接，不能为了消除报错任意补几个固定面。

力和压力也不能混用。总力施加在一个面上，与给这个面的每个节点重复施加同样大小的力，不是同一个载荷。分析者应核对最终合力、方向和作用区域，确认外载经过接触或连接传到预期支座，而不是被某个意外约束承担。

分析类型同样是条件的一部分。材料进入塑性、大转动使几何关系改变或接触状态变化时，简单线性假设可能不够；惯性影响明显时，需要考虑动态响应。材料参数与几何单位必须一致，例如采用N和mm时，弹性模量可用N/mm²，不能把以Pa表示的数字原样填入而不换算。

两段拉杆实例：从节点位移算到应力

以图中的一维均匀拉杆为例：总长L＝1000 mm，截面积A＝100 mm²，弹性模量E＝200000 N/mm²；左端轴向位移为0，右端承受向右F＝1000 N的力。假设材料线弹性、变形很小，只研究轴向拉伸，忽略自重、温度变化和夹持处的三维局部效应。以下是可手算的教学模型，不是企业试验或某款软件实测结果。

将它分为两个500 mm长的杆单元，得到左、中、右三个节点，位移分别记作u₁、u₂、u₃。根据杆单元的力—位移关系，每段轴向刚度k＝EA/l＝200000×100÷500＝40000 N/mm；单元内力取决于两端位移差。

左端已知u₁＝0。中间节点没有外加力，两个单元在这里的作用要平衡，因此有40000×（2u₂－u₃）＝0；右端满足40000×（u₃－u₂）＝1000。联立可得u₂＝0.025 mm、u₃＝0.050 mm。中间节点只有一个位移值，两段不能各给它一个不一致的位置，这就是共享节点组装的具体含义。

每段伸长均为0.025 mm，因此轴向应变ε＝0.025÷500＝0.00005；轴向应力σ＝Eε＝10 N/mm²，也就是10 MPa。支座对杆的反力向左，大小为1000 N，与向右的外力平衡。整杆公式δ＝FL/EA也得到0.050 mm，σ＝F/A得到10 MPa，分别提供位移与应力的交叉检查。

这个特殊模型的位移本来就是沿长度线性变化，线性杆单元能够表示它，所以不需要靠大量加密才能得到上述解。它适合解释组装与核对，却不能证明复杂零件也能用两个单元算准。如果要研究真实夹具附近的应力集中、横向变形或连接接触，就需要能表达那些现象的模型。

网格收敛为什么不等于越细越好

网格收敛研究是在材料、载荷、几何和边界条件保持一致时，逐次改变网格，观察关注结果是否趋于稳定。例如，始终比较同一点的位移或同一区域、同一输出方式的应力，并记录单元类型、网格尺度及目标值。COMSOL的网格细化说明指出，不同评价量可能以不同速度收敛；不能只看总位移变化很小，就宣布所有局部应力也已收敛。

可用相邻两次结果的变化量衡量敏感性，但允许变化多少应由问题、输入精度和决策需要确定，不存在适合所有分析的统一百分比。除缩小单元外，也可提高单元阶次或在梯度较大的区域局部细化。网格形状、连接质量和单元适用性同样重要；Ansys关于网格检查的说明明确区分单元形状与连接问题，单元外观合格并不保证模型整体连接正确。

还要区分求解器收敛与网格收敛。前者通常描述当前离散方程的迭代残差等达到规定条件，后者关注换用更细网格后结果是否仍改变。一个粗网格上的非线性方程可以迭代收敛，仍不足以说明空间离散已经准确；模型即使在两方面都收敛，也可能用了错误的材料或约束。

局部应力奇异是另一个常见问题。理想尖角、某些点载荷或突变约束可能让数学模型中的峰值随网格加密继续升高，不趋向一个有限值。Ansys的应力奇异性说明给出了这些情况。此时无限加密不会自动解决问题，也不能一律删掉红色区域：若关注的正是该处强度，应重新检查真实圆角、接触面积、连接和材料假设，必要时采用适当的局部模型或评价方法。

有限圆角处的真实应力集中与理想奇异点不同。前者在适当模型和网格下可以趋于有限值，后者的峰值行为来自所作理想化。判断依据应包括几何、载荷方式及细化后的变化，而不是仅凭“颜色很红”。

怎样判断有限元结果能不能用于工程决策

阅读结果时，先确认显示的是总位移、某方向位移、主应力、等效应力还是其他量，以及单位、坐标系和工况。软件常放大变形以便观察，云图中的弯曲形状不一定按真实比例显示；不同图片的自动色标范围也可能不同，红色只代表当前显示范围的一端，不直接等于失效。

随后检查数量级、反力平衡、变形方向和载荷路径。拉杆例中，如果右端拉力向右而整体却缩短，或支座反力与1000 N外力明显不平衡，就应检查符号、单位和施载方式，而不是先增加网格。复杂模型还要关注接触状态、求解警告，以及关键结果对连接刚度、材料和载荷不确定性的敏感程度。

有限元教学资料对verification与validation的区分有助于理解两类证据：数值验证关注方程是否被正确、充分地求解，可通过解析解、基准题、平衡与网格研究检查；模型确认关注数学模型是否足以代表实际对象，需要与相关工况的测量或试验比较，并考虑量测不确定性。拉杆手算相符属于前一类，不等于已经证明真实夹具和实物都符合模型假设。

一份可用于决策的分析记录，应能让别人重建所研究的工况、几何与连接简化、材料来源、载荷和约束、单元与网格、关键结果及验证依据。若结论是“该工况下变形满足要求”，就要给出对应位置的位移与允许值；若要判断屈服、疲劳或稳定性，则还需要各自适用的材料数据、载荷历史和评价方法。静力应力没有超过某个数字，并不能自动证明疲劳寿命或其他工况也满足要求。

参考资料

• COMSOL：有限元方法与网格细化，离散近似、单元及收敛评价。
• Jönköping University：Finite element analysis与The rod equation，组装、边界条件、杆单元和结果核验。
• Ansys：Assessing stress singularity，应力峰值与网格细化的关系。
• Ansys：Meshing Your Solid Model，网格形状、连接及分析结果检查。