FEM模型审查清单
理论与物理
概述 — 为何需要检查清单
模型审查检查清单具体检查什么内容？是在执行分析之前查看的东西吧？
是的。这是FEM模型质量保证的一部分，是在执行分析前系统性地检查网格质量、边界条件、材料属性、载荷条件的清单。NASA的"FEM最佳实践"（NASA-STD-5002A）规定了包含35项的检查清单，在航空航天结构分析中必须执行此清单。
有35项这么多吗！但我很少听说我们自己的团队会执行如此严格的检查...
遗憾的是，在许多实际工作中，往往变成"划分网格后扔给求解器，收敛了就OK"。但实际上，收敛 ≠ 结果正确。例如，在汽车碰撞模拟中，用Tied Contact近似螺栓连接但未经验证就继续推进，导致与实物试验结果完全不符的情况并不少见。
在自己的团队中，制定最低限度的检查清单并强制要求审查者签字确认是最佳实践。由另一位工程师而非分析者本人进行审查，可以防止因先入为主导致的错误。
原来如此，双重检查机制很重要啊。具体需要查看哪些项目呢？
大致可分为6个类别：
- 网格质量检查 — 纵横比、偏斜度、雅可比等指标
- 边界条件验证（BC Verification） — 约束的合理性、是否存在过约束
- 材料属性验证 — 杨氏模量、泊松比的值与单位制的一致性
- 载荷条件审查 — 载荷工况的完备性、载荷值的依据
- 自由-自由模态检查 — 通过无约束模态分析验证连接性
- 反力平衡检查 — 确认 $\sum F = 0$、$\sum M = 0$
检查清单的数学基础
检查清单给人一种"操作手册"的印象，它也涉及数学公式吗？
当然。许多检查项目都有定量的判定标准。例如反力平衡检查，用以下公式评估：
通常要求这个 $\varepsilon_{bc}$ 低于0.1%。此外，还需确认力矩平衡：
必须同时确认力和力矩两者平衡才行啊。即使力平衡了，如果力矩有偏差也没意义。
没错。例如，在悬臂梁末端施加集中载荷 $P$ 时，固定端的反力应为 $R_y = P$（力平衡），反力矩应为 $M = P \cdot L$（力矩平衡）。通过FEM输出结果来确认这一点。如果数值有偏差，则说明存在边界条件设置错误或载荷定义遗漏。
关于网格质量，整体刚度方程为：
这个刚度矩阵 $[K]$ 是由各单元的刚度矩阵 $K_e$ 组装而成。单元刚度通过数值积分计算：
这里 $|J(\xi_g)|$ 是雅可比行列式。单元发生畸变时，雅可比会退化，积分精度下降。所以在网格质量检查中要确认雅可比比。如果出现负雅可比，则该单元发生翻转，物理上无意义。
啊，了解了公式背后的原理，就明白为什么会有"纵横比小于5"之类的标准了。这不仅仅是经验法则，而是直接关系到数值积分的精度。
NASA的质量保证体系
NASA的检查清单具体是什么内容？您刚才说总共有35项对吧？
详情可参考NASA-STD-5002A《结构分析与试验要求》及其补充文档NASA-HDBK-5014。按类别总结主要项目如下：
| 类别 | 检查项目（示例） | 项目数 |
|---|---|---|
| 几何保真度 | 与CAD的几何差异、简化合理性的记录 | 4 |
| 网格质量 | 纵横比、偏斜度、雅可比、连接性 | 7 |
| 材料属性 | 杨氏模量/泊松比的来源、温度依赖性、单位制 | 5 |
| 边界条件 | 约束自由度数量、对称条件、接触定义 | 6 |
| 载荷条件 | 载荷工况完备性、安全系数的应用 | 5 |
| 分析执行 | 收敛判定标准、求解器警告的确认 | 4 |
| 结果验证 | 反力平衡、能量平衡、自由-自由检查 | 4 |
体系真完善啊。航空航天以外的领域也能用吗？
当然可以。汽车、核电、医疗器械等行业中，也有很多参考NASA框架制定自己检查清单的例子。例如，汽车行业Euro NCAP的碰撞分析指南中有类似的要求。关键在于确定"自己团队的分析中风险最高的是哪个部分"，并以此为重点制定检查清单。
反力平衡式的物理意义
- $\sum R_{reaction}$（反力总和）：所有约束点计算出的反力总和。无约束的自由度上反力为零。可直接从求解器输出文件（.f06, .sta, .rst等）获取。
- $\sum F_{applied}$（施加载荷总和）：用户定义的所有载荷的总和。分布载荷时按积分值计算。包含重力时需加上总质量 $\times g$。
- $\varepsilon_{bc}$（平衡误差）：此值较大时，表明可能存在边界条件定义遗漏、载荷作用面指定错误或数值误差。
雅可比与单元质量的关系
- 雅可比行列式 $|J|$ 表示等参映射的局部"畸变程度"。正方形→正方形的映射中 $|J| = $ 常数（理想值）
- $|J| \leq 0$ 的单元是"翻转单元"，物理上无意义。许多求解器会在分析前将其作为错误检测出来
- $|J|$ 的最小/最大比（雅可比比）低于0.3时，精度下降显著。六面体单元推荐0.5以上
数值解法与实现
网格质量检查
网格质量检查具体要看哪些指标、依据什么标准呢？不同软件输出的指标五花八门，很困扰。
主要有5个质量指标。这与求解器无关，是通用的：
| 指标 | 定义 | 理想值 | 合格标准 | 不合格会怎样 |
|---|---|---|---|---|
| 纵横比 | 最长边/最短边 | 1.0 | < 5.0 | 应力精度下降，特别是弯曲问题 |
| 雅可比比 | $|J|_{min} / |J|_{max}$ | 1.0 | > 0.3 | 积分精度下降，若为负值则是翻转单元 |
| 翘曲角 | 面的非平面度 | 0° | < 15° | 壳单元刚度评估错误 |
| 偏斜度 | 内角与理想角的偏差 | 0 | < 0.7 | 收敛性变差 |
| 锥度比 | 对边长度比的偏差 | 0 | < 0.5 | 应力分布畸变 |
所有单元都必须满足标准吗？复杂形状的话，好像总有一部分会超标...
问得好。现实目标是99%以上的单元满足标准。如果存在不合格单元，需确认该单元不在关注区域（应力集中区、破坏判定区域）。如果报告应力的区域存在不合格单元，则该结果不可信。
实际工作中，输出"质量直方图"来确认分布是常规做法。Hypermesh中使用Quality Index面板，Ansys中使用Mesh Metrics可以一键查看。
边界条件验证（BC Verification）
边界条件检查，感觉只要看看"这里固定了"不就行了吗？还需要确认那么多东西吗？
边界条件是FEM分析中最容易出错的地方。必须确认以下几点：
- 是否存在过约束：静定问题的3D实体，约束基本是6个自由度（3平动+3转动）。约束超过此数会产生人为应力
- 约束的位置和方向：螺栓孔内表面该固定的地方是否固定了外表面？圆柱坐标系设置是否正确？
- 对称条件的合理性：对称面上的面外位移是否被约束？反对称载荷下是否采用了反对称边界条件？
- 接触定义的合理性：Tied/Bonded接触的主从面选择、摩擦系数设置
过约束的话，固定得越多不是越安全吗？这样不行吗？
这是常见的误解。过约束会阻碍结构本应发生的变形，从而计算出现实中不存在的应力。例如，将热膨胀的管道两端全部自由度固定，会在实际法兰可能滑动的部位产生巨大的虚拟应力。结果被判定为"危险"，但实际上完全没问题——这是假阳性，会误导设计判断。
自由-自由模态检查
自由-自由模态检查是什么？"自由-自由"就是无约束的意思吧？为什么特意要解除约束呢？
正是如此。这是一种不施加任何边界条件而进行模态分析（特征值分析）的测试。目的有三个：
- 模型连接性验证：如果所有部件正确连接，则刚体模态只有6个（3平动+3转动）。如果出现7个以上，则说明有部件浮空（连接遗漏）
- 刚体模态频率确认：理论上6个刚体模态的固有频率为零。FEM中由于数值误差不会完全为零，但 $f_{rigid} < 0.01\,\text{Hz}$ 是目标。如果超过1Hz就有问题
- 质量合理性确认：观察自由-自由状态下第7个及之后的弹性模态振型，确认变形模式在物理上是否合理
诶，刚体模态出现7个就是连接遗漏？实际中有因此得救的案例吗？
在实际工作中经常遇到。特别是装配体模型有100个以上部件时，只要有一处Tied Contact定义遗漏，那个部件就会浮空。静力分析中可能会因"奇异矩阵"错误而停止，但根据载荷工况有时也可能偶然收敛。自由-自由检查是预先暴露这种"静默缺陷"的非常有效的方法。
方法很简单。Nastran中使用SOL 103且不指定SPC，Abaqus中使用*FREQUENCY步骤且不设边界条件。计算成本也低于常规分析的十分之一。
实践指南
材料属性验证
材料属性检查，不就是输入杨氏模量和泊松比吗？我以为直接用MIL手册的值就行了...
比起数值本身，单位制的一致性才是最可怕的。实际工作中最常见的致命错误是"弄错杨氏模量的单位"。例如：
| 单位制 | 杨氏模量（钢） | 密度（钢） | 常见错误 |
|---|---|---|---|
| SI (m, kg, s, Pa) | 2.1 × 1011 Pa | 7850 kg/m³ | 直接输入GPa值 → 偏差1000倍 |
| mm, t, s, MPa | 2.1 × 105 MPa | 7.85 × 10-9 t/mm³ | 密度仍输入kg/m³的值 |
| mm, kg, ms (碰撞) | 210 GPa → 210 kN/mm² | 7.85 × 10-6 kg/mm³ | 忘记时间单位ms |
哇，mm-t-s系的密度数值好小啊。如果弄错了会怎样？
密度弄错会导致惯性力、重力载荷、固有频率全部出错。静力分析中若无重力可能不会察觉，但在动力分析或模态分析中，结果会变得极其离谱。检查清单中需确认：
- 记录输入材料属性的来源（手册名称、试验报告编号）
- 计算模型总质量，与实物重量比较（误差2%以内为目标）
- 若有温度依赖性，确认使用温度范围内的数据是否正确输入
- 复合材料情况下，确认铺层顺序和纤维取向角是否与设计图纸一致
载荷条件审查
载荷条件审查...
なった
詳しく
報告