10节点四面体单元(TET10)
10节点四面体单元(TET10)的理论基础
TET10 — 3次元自动网格划分的主角
老师，TET4是"不要使用"，那TET10可以积极使用吗？
TET10是3次元结构分析的实务标准。几乎能对任意形状进行自动网格划分，精度也很高。可以说现代3次元FEM的大部分都是用TET10进行的。
TET10的结构
TET10和TET4有什么不同？
在TET4的4个顶点基础上，在每条边的中点加6个中间节点，共10个节点的四面体。每个节点3个自由度，共30个自由度。
形状函数是二次多项式：
二次的形状函数，所以单元内位移呈二次变化。应变线性变化。
是的。这是TET4（应变恒定）的决定性差异。弯曲产生的应力梯度可以用一个单元来表示。
精度的理论背景
为什么二次单元精度更高？
用泰勒展开就容易理解。当用多项式逼近真实位移场时：
- TET4（一次单元）近似到常数项和一次项 → 二次以上项是误差
- TET10（二次单元）近似到常数项、一次项和二次项 → 三次以上项是误差
收敛速度（网格细化时的误差减小速率）：
| 量 | TET4 | TET10 |
|---|---|---|
| 位移误差 | $O(h^2)$ | $O(h^3)$ |
| 应力误差 | $O(h)$ | $O(h^2)$ |
| 收敛速度 | 慢 | 快 |
TET10的应力误差以 $h^2$ 减小…网格减半时误差变为1/4。TET4减半只能减到1/2。
因此TET10即使网格较粗也能获得充分的精度。同样自由度数的情况下，TET10比TET4精度高得多。
数值积分
TET10的数值积分如何？
TET10的形状函数是二次多项式，B矩阵是一次。$B^T D B$ 是二次多项式，需要4点Gauss积分来准确积分。
| 积分格式 | 积分点数 | 精度 | 用途 |
|---|---|---|---|
| 4点（完全积分） | 4 | 准确积分二次多项式 | 标准 |
| 1点（低减积分） | 1 | 精度下降 | 特殊用途 |
TET10可以用低减积分吗？
几乎不会。与TET4不同，TET10的完全积分（4点）通常不会出现体积锁定。低减积分只会降低精度而没有优势。
不过Abaqus的C3D10M（Modified）采用特殊的积分方法。带hourglass控制的低减积分，提高了接触问题的稳定性。通常的C3D10和C3D10M的使用区分我后面会详细说明。
曲面的近似精度
TET10能准确表现曲面吗？
因为中间节点在边上，能将边变成曲线。通过将中间节点吸附到CAD曲面，能比TET4更准确地逼近圆柱和球面等曲面。
比如用TET4对圆孔建模时，圆会被近似为多边形。TET10的话，每条边都变成曲线，沿着圆。应力集中的评估精度直接关系到这一点。
如果不把中间节点吸附到CAD面会怎样？
中间节点会线性排列的"直边TET10"，曲面的近似精度会降低到TET4的水平。大多数前处理器会自动吸附到CAD面，但要确认设置。
总结
整理TET10的理论。
要点：
- 10个节点、二次形状函数、单元内应变线性变化 — 相比TET4质的提升
- 收敛快 — $O(h^2)$ 收敛应力。TET4的2倍阶数
- 自动网格容易 — 适配几乎任何3D形状
- 曲面表示准确 — 中间节点CAD吸附后能二次逼近曲面
- 3次元FEM的实务标准 — "不确定就用TET10"
TET10是兼有"自动网格容易"和"足够精度"的，实用上最平衡的单元呀。
正是如此。并非完美（HEX20效率更好），但包含网格生成在内的总体成本中TET10最优。
TET10二次形状函数推导
10节点四面体单元（TET10）具有4个顶点和6条边中点，共10个节点，形状函数为完全二次多项式（10项）。使用体积坐标（L1,L2,L3,L4，L1+L2+L3+L4=1）时，顶点形状函数为Li(2Li-1)，边中点为4LiLj。此定式由1960年代后期Zienkiewicz等在威尔士大学斯旺西校区系统化，因与任意形状部件自动网格划分（Delaunay分割）的契合度高而广泛普及。
10节点四面体单元(TET10)的数值计算方法
C3D10 vs. C3D10M
Abaqus有C3D10和C3D10M，有什么不同？
这是实务中最重要的使用区分之一。
| 特性 | C3D10 | C3D10M |
|---|---|---|
| 积分 | 4点完全积分 | 改进低减积分（hourglass控制） |
| 接触稳定性 | 略有不稳定 | 稳定 |
| 体积锁定 | $\nu > 0.45$ 时要注意 | 已对策 |
| 计算成本 | 低 | 略高 |
| 推荐场景 | 无接触，线性弹性 | 接触有，不可压缩材料，通用 |
C3D10M更通用，那总是用C3D10M就行？
Abaqus手册推荐C3D10M作为通用选择。C3D10的接触面容易出现棋盘纹样（应力振荡）。C3D10M解决了这个问题。
不过C3D10M有一个hourglass模式，在极少数情况下可能成为问题。C3D10有时会给出更"直爽"的结果。
各求解器的单元名
| 求解器 | 标准TET10 | 改进TET10 | 推荐 |
|---|---|---|---|
| Abaqus | C3D10 | C3D10M, C3D10MH | C3D10M |
| Nastran | CTETRA(10) | — | CTETRA(10) |
| Ansys | SOLID187 | — | SOLID187 |
| LS-DYNA | ELFORM=17 | — | ELFORM=17 |
Nastran和Ansys没有改进版？
Nastran的CTETRA(10)和Ansys的SOLID187内置了各自的稳定化方法，相当于C3D10M级别的对策已融入。名称没有分开，但内部上已改进了。
网格生成的最佳实践
生成TET10网格的要点？
1. 中间节点的CAD吸附 — 曲面上的中间节点投影到CAD面。忽略这点会降低曲面近似精度
2. 尺寸控制 — 应力集中部（孔、倒角、缺口）设置单元尺寸的局部控制
3. 单元品质确认 — 雅可比行列式 > 0.3、纵横比 < 5、最小角 > 10°
4. 网格渐变 — 细网格区域到粗网格区域的过渡要平缓。尺寸比相邻单元1.5倍以下为理想
网格渐变率多少合适？
相邻单元的尺寸比1.5倍以下推荐。超过2倍会在接口处精度下降。许多前处理器的自动网格有渐变率自动控制，但极端尺寸差时需要手动调整。
自适应网格细化
有自动细化网格的功能吗？
自适应网格细化（AMR）是基于误差估计对网格自动细分的方法。TET10因四面体分割容易，与AMR的相性好。
步骤：
1. 初始网格求解
2. 计算误差指标（应力的单元间不连续、ZZ误差估计等）
3. 对误差大的区域细化网格
4. 重新求解
5. 重复至收敛
Abaqus或Ansys能用AMR吗？
Ansys Workbench的Convergence功能最方便。指定目标应力变化率，会自动细化网格确认收敛。Abaqus的 *ADAPTIVE MESH 主要支持流体，结构的h-refinement多数手工。
总结
整理TET10实现细节。
要点：
- Abaqus选C3D10M — 接触稳定性、体积锁定对策
- 中间节点CAD吸附是必须 — 关系到曲面精度
- 网格渐变率1.5倍以下 — 尺寸激变会降低精度
- 自适应细化 — Ansys的收敛功能最方便
- 4点Gauss积分是标准 — TET10低减积分的优势少
TET10的4点Gauss积分
TET10标准使用4点或5点Gauss积分，能完全积分3次多项式。Keast（1986年）提出的4点积分格式使用体积坐标(0.1381966,0.1381966,0.1381966)等4点，能准确积分到4次多项式。11点格式精度进一步提高但计算成本约2.75倍。AbaqusC3D10M作为Modified TET，内置中点移动流控制以减轻体积锁定。
10节点四面体单元(TET10)的实务应用
TET10的实务工作流程
请教我用TET10的3次元分析的典型工作流程。
1. CAD导入＋简化 — 去除不必要的形状特征
2. 网格尺寸的规划 — 决定关注部位的必要尺寸
3. TET10自动网格生成 — 设置尺寸控制
4. 品质检查 — 雅可比行列式、纵横比
5. 求解分析
6. 结果评估 — 用非平均化应力分布检查品质
7. 网格收敛性验证 — 2水准应力变化5%以内
网格尺寸的确定方法
单元尺寸怎么决？
根据关注部位的形状特征：
| 形状特征 | 推荐单元尺寸（TET10） |
|---|---|
| 圆孔周围 | 孔半径 / 8 以下 |
| 倒角 | 倒角半径 / 3 以下 |
| 板厚方向 | 板厚 / 3 以下（最低2单元） |
| 螺栓孔 | 孔径 / 6 以下 |
| 焊接焊脚 | 脚长 / 3 以下 |
| 通用远场区域 | 最大尺寸 / 10 左右 |
倒角半径 / 3 的话，R5倒角要1.7 mm以下单元尺寸？
是的。倒角的应力集中需要充分的网格密度。倒角省略为尖边时会成为应力奇点，疲劳评估必要部位一定要保留倒角。
HEX单元的使用区分
TET10和HEX单元（六面体）怎么选择？
| 判断标准 | TET10 | HEX8/HEX20 |
|---|---|---|
| 复杂形状 | ○（自动网格） | ×（需手工网格） |
| 网格生成时间 | 短 | 长（数小时～数日） |
| 自由度效率 | 低（HEX的2～5倍） | 高 |
| 接触面稳定性 | 略低 | 高 |
| 不可压缩材料 | C3D10MH应对 | C3D8H应对 |
| 挤压形状 | △（非效率） | ○（sweep网格） |
"挤压形状"HEX有利呢。
2次元断面挤压形状（管、H形梁等）能用sweep网格高效生成HEX单元。TET10对这类形状建网格会不必要地增加自由度。
实务方法：
- 单纯形状·挤压形状 → HEX（sweep网格）
- 复杂形状·有机形状 → TET10（自动网格）
- 混合情况 → HEX和TET10混合网格（过渡用金字塔单元）
结果的品质指标
TET10分析结果的品质怎么判定？
非平均化应力的不连续最可靠。
对比相邻单元的应力值：
- 差5%以内 → 网格充分
- 差5～15% → 略粗。关注部位考虑细化
- 差15%以上 → 网格不足。需细化
大多数后处理器能显示非平均化应力？
Abaqus/CAE, Ansys Workbench, HyperView等主要后处理器都对应。Abaqus/CAE中"Average: 75%"等设置能改变平均化阈值。100%是完全平均化，0%是完全非平均化。
实务检查清单
请给TET10分析的检查清单。
"非TET4"是第一个确认项。和TET4页一样呢。
正因为那样重要。用TET10、检查网格品质、验证收敛性。守住这3步，3次元FEM的结果就信得过。
TET10自动网格生成工作流程
在复杂铸件（发动机缸盖等）FEM中，从CATIA/SOLIDWORKS形状自动生成TET10网格是业界标准。Altair HyperMesh的AutoMesh功能仅设置目标单元尺寸，即可用Delaunay-Advancing Front法自动生成TET10。2mm设定下发动机缸盖全体（体积约2L）网格生成在2024年16核工作站上约5分钟完成。
10节点四面体单元(TET10)的软件比较
TET10单元在各求解器中的特点
各求解器的TET10单元有差异吗？
基本精度基本相同，但附加功能有差异。
| 功能 | Abaqus C3D10M | Nastran CTETRA(10) | Ansys SOLID187 |
|---|---|---|---|
| 接触稳定性 | 最高（改进积分） | 良好 | 良好 |
| 体积锁定对策 | C3D10MH（混合） | 有限 | Mixed u-p对应 |
| 大变形 | NLGEOM对应 | SOL 400对应 | NLGEOM对应 |
| 自适应细化 | 有限 | 有限 | Convergence功能 |
| 陽解法 | C3D10M（Explicit对应） | — | — |
Abaqus的C3D10M在陽解法也能用，很稀有呢。
TET10通常在陽解法中难用（稳定时间增分会变小）。C3D10M的改进积分考虑了陽解法稳定性，碰撞分析也能用TET10。
网格生成工具
TET10网格的品质很大程度取决于网格生成工具，对吧？
| 工具 | TET10品质 | 特点 |
|---|---|---|
| Abaqus/CAE | 良好 | Advancing front算法 |
| Ansys Meshing | 优秀 | Patch conforming + Independent。品质指标表示充实 |
| HyperMesh | 良好 | TET-remesh功能。品质改善工具充实 |
| GMSH | 良好 | 开源。Netgen/Delaunay算法 |
| Simlab | 优秀 | Altair出品。自动网格品质高 |
Ansys Meshing的Patch conforming和Independent的区别？
通常推荐Patch conforming。形状复杂导致网格生成失败时试试Independent。
选型指南
TET10的选型指南？
TET10在任何求解器里信度都高呢。
是的。TET10是各求解器最用力开发的单元之一。自动网格用户大多用TET10，求解器商也在品质提升上投入。
TET10主要求解器实现一览
对应TET10的单元为Nastran的CTETRA(10)、Abaqus的C3D10/C3D10M、Ansys的SOLID187、Altair OptiStruct的CTETRA10。Abaqus的C3D10M用改进的中点位置定式来减轻体积锁定，对橡胶和生物软组织分析特别有效。Siemens NX Nastran 2023版TET10的并行求解器效率比前版本提升35%，超1000万单元的模型也能在实用时间内求解。
10节点四面体单元(TET10)的前沿研究
TET10的前沿研究
关于TET10有前沿的研究吗？
自动网格品质提升和新计算方法是主流。
p-法对TET10的高阶化
p-法是什么？
相对通常的h-法（细化网格提高精度），p-法是不改网格而提高单元多项式次数来提高精度。TET10（p=2）提升为p=3, p=4等。
p-法的优点：
- 无需重新生成网格
- 指数级收敛速度（光滑解）
- 应力集中部精度飙升
商用求解器能用p-法？
Simcenter Nastran的p-法单元（SOL 147/148）和MSC Nastran的p-单元对应。Ansys也有p-法选项。不过p-法与通常h-法工作流不同，普及有限。
四面体重新网格划分（TET-remesh）
大变形分析中单元变形时，分析途中对TET10网格自动重新划分的方法。ALE（Arbitrary Lagrangian-Eulerian）法的结构版。
在锻造或挤压成形这样的大变形问题里有效。Abaqus的ALE自适应网格或Defrom的自动重网格已实用化。
通过机器学习优化网格质量
有用AI改善网格生成的研究？
用GNN和CNN：
- 预测最优网格密度分布 — 从过往分析学习，首次就配置恰当密度
- 网格品质事前检查 — AI检测生成网格的品质问题
- 自适应细化指导 — AI推荐哪里应该细化
首次就配置对的网格的话，反复验证收敛性就不用了呢。
理想的话是如此。但现在还没到"AI网格设计超过人的设计能力"的阶段。研究很活跃，数年内有可能有商用的AI网格助手内置于CAE软件。
Poly-TET（多面体单元）
不拘泥四面体，用任意多面体单元的方法也在研究。VEM（Virtual Element Method）或Polygonal FEM这个系统。TET10的"对网格品质敏感"的弱点，想用单元形状自由度来解决。
总结
整理TET10的前沿研究。
TET10是"今日标准"，将来有可能通过AI和p-法"自动获得最优精度"的时代来临。
TET10与HEX20精度比较研究
TET10与HEX20（20节点六面体）的精度比较是FEM研究的经典课题。1990年代的参数研究（Cifuentes & Kalbag, 1992）表明同自由度数下HEX20在应力精度上约2倍优位，但实际复杂形状上自动TET10网格的通用性更高。2000年代以后，改进TET10（Abaqus公司的C3D10M等）体积锁定对策施行，与HEX20的精度差缩小到5%以内。
10节点四面体单元(TET10)的故障排查
TET10的故障
TET10也会出故障？比TET4好得多吧。
TET10优秀但没有零故障。要知道的问题有。
中间节点位置不正
中间节点的问题具体是什么？
中间节点偏离CAD面或相邻单元的中间节点不一致时出现。
症状：
- 曲面上应力呈锯齿状
- 单元雅可比行列式变负（单元翻转）
- 位移局部异常
对策：
- 确认前处理器的"中间节点CAD吸附"设置
- 品质检查中确认没有负雅可比行列式单元
- 有问题时局部重新生成网格
负雅可比行列式的单元能计算吗？
大多数求解器把负雅可比行列式的单元作为错误拒绝。Nastran的FATAL错或Abaqus的"Negative eigenvalue"会出现。一个负雅可比行列式单元就会让分析无法进行。
自由度数过多
TET10网格的自由度数太多，内存不足了。
TET10是TET4的2.5倍自由度（10节点 vs 4节点），HEX8的约3倍。大规模模型中自由度数成为问题。
对策：
- 对称条件活用 — 自由度减到1/2～1/8
- 子模型化 — 全体粗网格，局部精细
- 迭代法求解器 — 比直接法内存效率好
- HEX和TET10混合网格 — 单纯形状改HEX降低自由度
- 自适应网格 — 从粗网格开始，只对必要部位细化
混合网格的接接怎么做？
TET10和HEX20接接要用金字塔单元（PYRA13/14）夹在中间。金字塔单元有四角形面（HEX侧）和三角形面（TET侧）。主要前处理器能自动生成。
接触面不稳定性（C3D10的情况）
C3D10做接触分析应力会棋盘状。
C3D10的接触面稳定性问题。4点完全积分的C3D10中接触面压力振荡。
对策：
- 切换到C3D10M — 改进积分接触面稳定
- 接触面网格细化 — 缓解振荡
- 接触算法改变 — 罚函数法→Lagrange乘数法
C3D10M是接触问题的铁板呢。
Abaqus中包含接触的TET10分析没有理由不用C3D10M。接触不稳定是可预测的问题。
应力精度低于HEX
同样问题用HEX20和TET10解，TET10精度更低。
同自由度数下HEX20精度更高是事实。TET10要达到HEX20同等精度，需要自由度增加2～5倍。
但包含网格生成时间的总成本中TET10往往有利。HEX网格要花1天，TET10自动网格几分钟完成。
"HEX精度更高"和"TET10更实用"矛盾不？
不矛盾。精度效率（单位自由度精度）HEX优。工作流效率（网格生成＋分析＋评估总时间）TET10优。哪重视取决于项目。
总结
整理TET10故障处理。
C3D10M的接触稳定性和中间节点CAD吸附是2大关键点呢。
掌握这两点，TET10不会出大问题。TET10是实务中最信赖的3D单元。
TET10中点逸脱和负雅可比行列式
TET10的辺中点置于辺长25%以下或75%以上位置时会出现负雅可比行列式。HyperMesh的Quality Check中提取Jacobian<0.01的单元并修正是实务定石。Nastran中输出Fatal 2186。TET10网格自动生成时曲面吻合时最容易在曲率大的部位发生负雅可比行列式，HyperMesh的Local Remesh（局部重网格）重新生成最高效。
相关话题
价值
更详细
错误