20节点六面体单元(HEX20)
20节点六面体单元(HEX20)的理论基础
HEX20 — 最高精度的3维单元
HEX20是FEM中3维单元精度最高的吗？
作为实用单元来说精度最高。20节点的二次六面体单元，每个节点3个自由度，共60个自由度。HEX8（24个自由度）的2.5倍DOF，但精度高出数个数量级。
形状函数
HEX20有8个顶点节点和12个边中点节点。形状函数是三次曲面（Serendipity）型：
顶点节点：
边中点节点（例如: $\xi_i = 0$）：
Serendipity型是什么？与Lagrange型不同吗？
Lagrange型（27节点HEX）在面上和内部也有节点，而Serendipity型（20节点HEX）仅在边上有中间节点，面中心和内部没有节点。Serendipity型DOF较少，精度与Lagrange型相近，实务中应用最广。
HEX20的精度
HEX20的精度水平如何？
与TET10对比：
| 相同DOF数下的精度 | HEX20 | TET10 |
|---|---|---|
| 位移 | 极高 | 高 |
| 应力 | 极高 | 高 |
| 效率（单位DOF精度） | 最高 | HEX20的50~70% |
收敛速度：
- 位移: $O(h^4)$（HEX20）vs. $O(h^3)$（TET10）
- 应力: $O(h^3)$（HEX20）vs. $O(h^2)$（TET10）
HEX20收敛速度快1个数量级。即使网格粗一半也能达到同等精度…。
所以在能构建HEX网格的形状上，HEX20最高效。问题在于网格生成的工作量。\"精度效率 vs. 网格生成成本\"的权衡，这就是HEX20与TET10选择的本质。
积分方案
HEX20的积分：
| 积分 | Gauss点 | 特性 |
|---|---|---|
| 完全积分（3×3×3） | 27 | 最高精度。锁定风险小 |
| 减缩积分（2×2×2） | 8 | 避免剪切锁定。沙漏模式1个 |
HEX20的减缩积分是2×2×2 = 8点？与HEX8的完全积分相同。
对的。HEX20的减缩积分（C3D20R）用8个积分点，精度远高于HEX8的完全积分。C3D20R是3D FEM中效率最高的单元之一。
C3D20R有沙漏问题吗？
HEX20的减缩积分仅有1个沙漏模式。与HEX8的12个相比极少，实务中基本无问题。只是单个单元的补丁测试需谨慎。
何时使用HEX20
什么时候应该用HEX20？
反过来HEX20不利的场合是？
总结
整理一下HEX20的理论。
要点：
- 20节点、Serendipity型、二次六面体 — 实用3D单元最高精度
- C3D20R（减缩积分）最推荐 — 8个积分点，最高效
- 收敛快 — 比TET10快1个数量级
- 网格生成是课题 — 自动HEX困难。适合Sweep/映射网格
- 不适合大变形和显式法 — HEX8R或TET10更稳定
HEX20是"只要能建就最强的单元"。
正是。\"最高精度\"与\"能建出来\"之间隔着网格生成这堵墙。越过这堵墙的工程师，HEX20是最强武器。
二次六面体单元的理论优越性
20节点六面体单元（Serendipity单元）由Ergatoudis等人在1966年提出，是等参数化表述的代表作。完整的二次多项式包含在其形状函数中，使单元内应力精度接近一阶单元的8倍。1970年代阿波罗计划后，美国航天局的航天器结构分析采用了该单元，一度被称为\"黄金标准单元\"。
20节点六面体单元(HEX20)的数值计算方法
C3D20R — 最高效的3D单元
C3D20R（Abaqus中HEX20的减缩积分）被称为\"最高效\"的理由是什么？
有3个原因：
1. 精度/DOF比最高 — 相同DOF数下精度是TET10的2~3倍
2. 无剪切锁定 — 减缩积分已规避
3. 体积锁定也很轻 — $\nu = 0.49$ 左右无问题
为什么推荐C3D20R而不是C3D20（完全积分）？
C3D20用27个积分点，计算成本高，而且当$\nu$较大时会出现体积锁定迹象。C3D20R只需8个积分点，计算快，锁定少。实务中C3D20R绝对占主导。
求解器的单元命名
| 变体 | Abaqus | Nastran | Ansys |
|---|---|---|---|
| 完全积分 | C3D20 | CHEXA(20) | SOLID186(full) |
| 减缩积分 | C3D20R | — | SOLID186(red.) |
| 混合型 | C3D20H, C3D20RH | — | u-P对应 |
Nastran的CHEXA(20)没有减缩积分选项？
Nastran的CHEXA(20)默认使用2×2×2（减缩）积分。切换到完全积分有选项，但默认就是减缩。即Nastran的CHEXA(20) ≈ Abaqus的C3D20R。
中间节点的处理
HEX20的中间节点在边的中点，在曲面上要投影到CAD面，对吧？
与TET10相同。将曲面上的中间节点投影到CAD面，能实现对曲面的二次近似。但HEX20有特殊注意点：
中间节点位置偏离太远会导致单元退化。曲率大的部分若把中间节点投影到CAD面，使其偏离边中点太远，雅可比行列式会变负。
如何应对？
对策：
- 曲率大的部分网格要细致，使中间节点偏移相对较小
- 必须进行雅可比行列式检查
- 用预处理软件的\"mid-node projection limit\"设置限制投影量
HEX20的应用限制
HEX20有不擅长的问题吗？
大变形分析中HEX20处于劣势。中间节点缺乏物理约束，大变形时容易飘向\"古怪位置\"。
金属成形等大变形问题应：
HEX20主要用于线性~中等非线性（接触、小变形塑性）。
总结
整理HEX20的数值方法。
要点：
- C3D20R（减缩积分）最推荐 — 8个积分点最高效
- Nastran的CHEXA(20)默认减缩积分 — 相当于C3D20R
- 留意中间节点投影量 — 曲率大部分易产生负雅可比
- 大变形用HEX8 — HEX20是线性~中非线性用
- 完全积分C3D20无需用 — 通常C3D20R就够
27点vs8点Gauss积分的对比
HEX20单元完全积分需3×3×3=27个Gauss点。计算成本高。1974年Barlow点理论指出14点（减缩积分）也能保精度，ANSYS SOLID186用3×3×3，而Abaqus C3D20R采用2×2×2（8点）的减缩积分，计算时间缩短约40%，精度基本不变。
20节点六面体单元(HEX20)的实务应用
HEX20的实务应用
HEX20在实务中的应用场景是什么？
HEX20大显身手的场景：
压力容器的应力分类
ASME Div. 2的应力分类需板厚方向应力分布精度高。HEX20仅需板厚方向2个单元就能准确捕捉应力梯度（C3D20R）。TET10则需4个以上。
基准验证
FEM分析精度检验（NAFEMS基准等）常用HEX20作为参考单元。先用HEX20求\"最准确的结果\"，再拿TET10或壳单元结果对比。
接触分析的精密评估
螺栓连接部或轴承接触压力分布精密评估时，HEX20接触面比TET10更稳定光滑。粗网格也能得到相对准确的压力分布。
网格划分策略
如何高效生成HEX20网格？
| 方法 | 适用形状 | 效率 |
|---|---|---|
| Sweep网格 | 拉伸体（管、轴） | 极高 |
| 旋转网格 | 旋转体（法兰、圆盘） | 极高 |
| 映射网格 | 矩形区域 | 高 |
| Multi-zone | 稍复杂形状分区 | 中等 |
| 手动（HyperMesh等） | 任意形状 | 低（工作量大） |
旋转体用HEX20效率最高。
压力容器或汽轮机转子这样的旋转体，可从2D断面旋转生成HEX20。这种情况下HEX20精度效率最大化。
从TET10到HEX20的转换判断
用TET10做过的分析想用HEX20验证，怎么办？
实务上用\"部分子模型\"方法最适合：
1. 用TET10做全体模型分析
2. 提取关注部位的边界位移
3. 将位移作为边界条件加到HEX20子模型
4. 比较HEX20与TET10结果
不用全部改成HEX20？
全改成HEX20的网格生成成本太大。\"必需精度的部分用HEX20\"这样的分部法，成本-效益最优。
实务检查清单
请给出HEX20的检查清单。
HEX20是\"能建就最强\"。检查清单也很简明。
HEX20自身故障少。问题大多出在网格生成阶段。网格一旦成功，HEX20是最可靠的3D单元。
涡轮增压器外壳分析
汽车涡轮增压器压气机外壳形状复杂，用HEX20约15万单元的模型，与四面体单元60万单元精度相当，但计算时间缩短为1/4。IHI技术报告（2019）验证了这一点。因此设计初期优化循环普遍选用HEX20。
20节点六面体单元(HEX20)的软件比较
HEX20在各求解器中的特性
各求解器的HEX20有区别吗？
基本精度相同。差异在附加功能和大规模性能。
| 功能 | Abaqus C3D20R | Nastran CHEXA(20) | Ansys SOLID186 |
|---|---|---|---|
| 减缩积分 | 明确选择 | 默认 | KEYOPT(2)选择 |
| 混合型 | C3D20RH | — | u-P对应 |
| 温度依赖 | 完全支持 | 支持 | 完全支持 |
| 大规模并行 | MPI支持 | MPI + GPU | MPI + GPU |
| 复合材料 | — | — | SOLID186 Layered |
Ansys的SOLID186 Layered是什么？
Ansys的SOLID186带Layered选项，支持层状结构（复合材）。通常的HEX20再加积层方向和各层材料定义，便能处理复合材的实体建模。
网格生成工具的比较
| 工具 | HEX20网格质量 | 特点 |
|---|---|---|
| HyperMesh | 最高（手动控制） | HEX网格的专业工具 |
| Ansys Meshing Multi-zone | 良好（半自动） | GUI易操作 |
| Abaqus/CAE Sweep | 良好 | Sweep网格稳定 |
| Cubit/Trelis | 良好 | 研究用。高级分割功能 |
| Gmsh | 基本 | TET→HEX转换有限 |
HEX网格生成中HyperMesh最强？
手动HEX网格的品质和灵活性，HyperMesh压倒优势。汽车或航空的碰撞模型几百万个HEX单元手工作业，用的就是HyperMesh的技能。
选型指南
总结一下？
HEX20是\"能建就最强\"这个结论没变。
是的。HEX20的单元技术本身在各求解器中都很成熟。瓶颈永远是网格生成。
HEX20各求解器的对应状况
Abaqus C3D20/C3D20R、ANSYS SOLID186、NX Nastran CHEXA（20节点）、LS-DYNA ELFORM=3都是20节点六面体，但节点序编号不同。Abaqus与ANSYS节点顺序有差异，混合模型转换时用HyperMesh的OptiStruct转块出错过。2015年前后这是常见FAQ。
20节点六面体单元(HEX20)的先进研究
HEX20的先进研究
HEX20有最新研究吗？
HEX20本身是成熟单元，但网格生成自动化和高阶化有新研究。
HEX27(Lagrange型) vs. HEX20(Serendipity型)
HEX27与HEX20比怎样？
HEX27是在面中心和体内部加节点的完全二次六面体。Serendipity型的HEX20相比，形状函数\"更完美\"，数学上更理想。
但实用上差异小。HEX27的DOF是27×3=81，比HEX20的60多35%。精度提升微小，不值得DOF增加。所以HEX20（Serendipity）是实务标准。
谱单元法
谱单元法是高阶多项式（p=4~8）的FEM变种。在HEX单元上放置高阶GLL（Gauss-Lobatto-Legendre）点。地震波传播模拟（SPECFEM3D等）的标准做法。
相当于把HEX20的p=2提升到p=8这样？
是的。谱单元自然得到对角质量矩阵，显式时间积分高效。地震学大规模波动模拟基本唯一选择。
IGA(等几何分析)的HEX
等几何分析（IGA）用NURBS/B-spline基替代Lagrange/Serendipity基。IGA的HEX单元：
- 与CAD形状完全一致（网格形状近似误差为0）
- $C^1$ 及以上连续（通常FEM是 $C^0$）
- 应力场光滑
应力的光滑性大有优势。
普通HEX20应力在单元间不连续（$C^0$），但IGA有 $C^1$ 以上连续，单元间应力光滑。应力集中部评估精度提高。
总结
整理HEX20的先进研究。
HEX20是成熟单元，高阶化和IGA开拓新可能。
非线性变形与HEX20的局限
HEX20在大变形分析时体积锁定倾向强。1984年Simo等提出F-bar法，独立处理体变形，克服锁定。Abaqus的C3D20H采用混合变形-压力（Hybrid），在橡胶、金属塑性大变形中精度比C3D20R高10~20%。
20节点六面体单元(HEX20)的故障排除
HEX20的故障
HEX20也会出故障吗？
HEX20是最稳定单元之一，但网格质量引发的问题仍有。
雅可比行列式为负（单元退化）
雅可比负是什么时候？
中间节点位置不当时发生。HEX20中间节点应在边中点附近，但曲面投影若使其偏离太远，雅可比会变负。
经验规则：中间节点应在边长的25%~75%范围内。超出这个范围雅可比易负。
如何应对？
完全积分的体积锁定
HEX20的完全积分（C3D20）会体积锁定吗？
当$\nu > 0.49$时可能发生。C3D20的27个积分点过多，体积拘束过剩。
应对：
- 改用C3D20R（减缩积分） — 8个积分点，拘束缓解
- 用C3D20RH（混合+减缩） — 非压缩材最稳定选择
结果是C3D20R最安全。
是的。C3D20（完全积分）除了\"害怕减缩积分沙漏\"才用，而HEX20的减缩积分沙漏实质无害。C3D20R作默认选择没问题。
DOF过多
HEX20模型DOF太多，内存不足。
HEX20每单元60个自由度，HEX8（24个自由度）的2.5倍。
应对：
\"全体HEX8I + 关注部位HEX20子模型\"是高效方案。
正是。HEX20精度需求通常有限。全体改HEX20网格成本巨大，用子模型优化最实际。
总结
整理HEX20的故障处理。
HEX20故障少于HEX8。单元品质高。
正确。HEX20的故障大多出在网格生成阶段。网格成功，HEX20是最可靠3D单元。
中点节点的品质检查
HEX20的中点节点偏离边中点超过25%则雅可比行列式变负，导致分析错误。商用网格划分工具有\"Jacobian Ratio\"指标，推荐范围1.0~0.6，特别是从导入的STEP数据自动网格时必检。HyperMesh从2000年代起已将该指标纳入标准品质检查。
相关主题
细节
错误