静态有限元分析

我第一次做 静态有限元分析是算一个压力容器的筒体应力，用一阶四面体自动网格，最大应力 140 MPa，评审一句”你这网格够吗”让我重做，换成二阶单元并加密后应力变 95 MPa——原来是坏单元在应力集中处造了个假峰。那次让我把 静态有限元分析 的第一步从”跑求解器”改成”先看网格报告”。

一、为什么 静态有限元分析先卡网格

线性四面体（C3D4）在工程应力场里基本是劝退级，易体积自锁、应力平滑差。我做结构件默认二阶四面体（C3D10）或六面体扫掠，薄壁件优先中面+壳单元。曾经一个齿轮用一阶四面体，齿根应力比二阶高 60%，完全失真。做 静态有限元分析时，单元类型比加密整体网格更该先定，坏单元不会因加密消失。

二、核心原理：离散与积分阶次

静态有限元把结构用形函数插值，单元内应变靠 B 矩阵从节点位移反推。一阶四面体单积分点会出现体积自锁，刚度虚高；二阶单元配足够积分点才稳。关键认知：应力是位移的导数，对网格更敏感，所以应力集中处必须局部加密而非全局加密。我用应力而非位移作为收敛指标，因为位移收敛了应力可能还没。

三、关键技术要点：应力集中与子模型

静态有限元分析最致命的是应力集中（倒角、孔边、焊缝）被粗糙网格抹掉或夸大。我处理孔边应力默认用子模型（submodeling）：全局粗算定位，局部切出细分网格映射边界，得到可信峰值。另一个坑是单元长宽比过高，薄壁转角处单元畸变导致应力尖峰。做之前我必查单元质量三项：长宽比、雅可比、翘曲度。

四、实操流程：从 CAD 到可信应力

我的标准动作：CAD 清理 → 简化非关键特征 → 单元选型（二阶/六面体优先）→ 局部加密应力集中区 → 约束加载 → 求解 → 网格无关性验证（尺寸减半应力变化<3%）→ 必要时子模型。最容易被砍的是无关性验证，但 静态有限元分析 的结论稳健性全靠它，我默认保留。

补一个实算：压力容器筒体，全局网格最大应力 140 MPa（含坏单元假峰），网格无关性后降到 98 MPa，子模型处理接管处应力集中得到 165 MPa 真实峰值。如果只交第一版，你会要么被假峰吓到过度设计，要么漏掉真实峰值低估风险。这之后我做应力分析，无关性+子模型是标配。

五、常见踩坑：约束与单位

静态有限元分析常见错是过约束制造假安全，或欠约束导致刚体位移求解报错。单位制混乱更是老生常谈，mm 模型按 m 给密度结果差千倍。还有人把接触当绑定处理所有连接，漏掉摩擦和预紧，连接刚度虚高。这些错都不在求解器，在建模前没想清楚。

六、复盘：有限元是流程不是黑箱

回过头看，静态有限元分析最值钱的是用可控的流程把”应力多少、可不可信”量化到能写进报告，但前提是网格和收敛做扎实。我接项目默认先问应力集中区在哪、要不要子模型，再定网格策略。把单元、无关性、子模型三件事做对，应力才敢拿去评审。被证明有用的，是愿意为一个孔边多切一层子模型的扎实。下次有人给我一张”最大应力”的图，我会先问——你的网格无关性做没做。