壳单元、四面体与六面体网格应该如何选择？

同一个零件用壳单元、四面体或六面体建模，结果可能明显不同，而每个模型都能正常求解。本文从关注结果、几何形态和建模成本出发，梳理三类单元的适用范围、常见陷阱与选择流程。

同一个零件，分别用壳单元、四面体和六面体建模，算出的位移、应力和固有频率可能明显不同，而每个模型都能正常求解。单元类型选错通常不会报错，而是带来系统性偏差：要么刚度偏大、应力偏低，要么关注部位根本表达不了真实的应力状态。

本文从几何形态、关注结果和建模成本三个角度，梳理三类单元的适用范围与常见陷阱。讨论以 Simcenter Nastran、Simcenter Optistruct 等隐式结构求解器为背景，方法同样适用于其他有限元软件。

本文要点

  • 先定要看的结果，再看几何：薄壁件优先用壳单元，中面质量决定模型质量；
  • 实心件能扫掠就用六面体，不能就用二阶四面体，一阶四面体只用于粗估或填充；
  • 壳与实体交界不能直接共节点，应使用 RBE3、MPC 等连接传递弯矩；
  • 下结论前至少加密一次做收敛性检查，并区分真实应力集中与应力奇异点。

一、选单元之前，先回答三个问题

  1. 要看什么结果？整体刚度、模态和变形，与圆角根部、焊缝、接触面附近的局部应力，对单元的要求完全不同。
  2. 几何是什么形态？厚度远小于面内尺寸的薄壁件、实心块体、细长杆件，对应的理想化方式不同。
  3. 能投入多少建模时间？四面体几乎全自动；六面体往往要先切分几何；壳单元需要抽中面、赋厚度。

单元类型是这三个问题折中的结果，脱离具体问题谈“最好的单元”没有意义。

二、壳单元：薄壁结构的首选

壳单元把三维薄壁零件理想化为“中面 + 厚度属性”。西门子 Simcenter 博客指出，薄壁零件通常用壳单元比实体单元更合适，前提是先在内外表面之间抽取中面。钣金件、冲压件、车身与机身蒙皮、注塑薄壁件都属于这一类。工程上常用的经验是：壁厚小于面内特征尺寸的约 1/10 时，优先考虑壳单元。处在两者之间的中厚件，可以用实体单元并保证厚度方向有足够层数，也可以对比壳与实体两种模型的结果。

壳单元的优势在于：

  • 自由度少：厚度方向不划分网格，同等精度下模型规模远小于实体网格。
  • 弯曲表达好：Simcenter Nastran 的 CQUAD4/CTRIA3 基于 Mindlin–Reissner 壳理论，可计入横向剪切变形（需在 PSHELL 中定义 MID3），薄板弯曲不依赖厚度方向的网格层数。
  • 便于优化：厚度是属性参数，在 Simcenter Optistruct 中可以直接作为尺寸优化的设计变量。
使用提示

  • 中面质量决定模型质量。厚度变化、加强筋和凸台处的中面断开或偏移，会直接改变刚度与质量分布。
  • 四边形优先，三角形远离关注区。西门子 Nastran 单元库说明，CTRIA3 可能表现出过大的刚度，尤其是薄膜应变；建议把它放在非关注区域，并尽量接近等边。
  • 局部应力会失真。壳模型表达不了厚度方向的应力，也没有真实的圆角：圆角根部的应力可能偏低，壳与壳的交线处又可能因中面重叠出现奇异的高应力。
  • 法向转动没有刚度。CQUAD4/CTRIA3 不为垂直于壳面的转动自由度提供弹性刚度，共面模型需用 AUTOSPC 或 PARAM,K6ROT 处理，避免刚度矩阵奇异。

三、四面体单元：自动化程度最高，但应默认用二阶

四面体几乎可以对任意复杂实体自动划分，是铸件、锻件、机加工件等实心几何最常用的单元。它的关键问题在阶次。

TET4

一阶四面体

常应变单元，单元内应变恒定。在弯曲和应力梯度明显的区域表现得过刚，位移和应力偏低，网格加密后收敛也慢；对橡胶、大塑性变形等近不可压缩问题，还会出现体积锁定。只用于粗略估算或非关注区的填充。

TET10

二阶四面体

增加了边中节点，应变可以在单元内线性变化，弯曲和应力集中的表达明显改善。在隐式结构分析中，应作为四面体网格的默认选择。

使用二阶四面体时还要注意：

  • 节点数增长快：同样的单元尺寸下，TET10 的节点数远多于一阶单元。大装配应控制全局尺寸，只在关注区加密。
  • 薄壁部位要保证层数：用实体单元表达薄壁时，厚度方向单元太少会让弯曲偏硬。二阶单元在厚度方向只有一层时，整体变形有时仍可接受，但应力梯度大的部位要加密，并做收敛性检查。
  • 几何清理仍然必要：微小倒角和碎面会生成极小或畸形的单元，影响结果精度并拖慢求解。

四、六面体单元：效率和稳定性最好，建模成本最高

在规则几何上，六面体可以用更少的节点达到同等精度。单元排列整齐，便于在厚度方向分层，也便于沿应力梯度方向控制尺寸。它常用于：

  • 可以扫掠或映射的几何，如挤出件、回转件、厚板与法兰；
  • 接触面、螺栓连接、焊缝等需要规则网格的局部区域；
  • 大变形、塑性等非线性分析，以及需要稳定应力结果的疲劳热点区。

代价是建模时间。复杂零件通常要先切分成可扫掠的子体，再逐块划分。做不到全六面体时，可以用楔形单元（Nastran 中为 CPENTA）处理三角形底面的扫掠区，用金字塔单元（Nastran 中为 CPYRAM）在六面体与四面体之间过渡。

使用提示一阶六面体（8 节点）同样有阶次问题：全积分形式在弯曲中容易出现剪切锁定，减缩积分形式可能出现沙漏模式。求解器一般通过改进的单元公式或积分方案缓解这些问题，但厚度方向仍要有足够的单元层数；必要时改用二阶六面体（20 节点）。

五、单元对照

单元 适合的几何 优势 主要风险 建模成本
壳单元（CQUAD4/CTRIA3 等） 薄壁、钣金、蒙皮 自由度少、弯曲好、厚度可优化 中面质量；缺少厚度方向应力 中（抽中面、赋厚度）
二阶四面体（TET10） 复杂实心件 几乎全自动 节点数多；薄壁层数不足 低
一阶四面体（TET4） 仅用于填充或粗估 划分最快 过刚、应力偏低 低
六面体（HEX8/HEX20） 可扫掠的规则几何；关注区局部 节点少、结果稳定 一阶的锁定或沙漏；切分费时 高

六、混合建模：把单元用在该用的地方

实际工程很少只用一种单元，常见做法有两种：

SHELL + SOLID

壳—实体混合

大面积薄壁用壳单元，厚重或关注的局部用实体单元。实体单元节点只有三个平动自由度，壳单元节点还有转动自由度，直接共节点会在交界处形成“铰接”，弯矩传不过去。交界处应使用 RBE3、MPC 等连接，或把壳单元延伸进实体内部。

HEX + TET

六面体为主、四面体填充

关注区和规则区域用六面体，复杂的过渡区用四面体，两者之间用金字塔单元或绑定接触连接。

经验提示混合模型最容易出问题的地方是交界。应单独检查交界处的位移和应力是否连续，交界本身不应成为应力最高点。

七、一套可执行的选择流程

前三步决定用什么单元，后三步决定结果能不能用。下面的选择器先给出起步单元，关注区仍按第 3 步单独处理。

交互工具

单元类型选择器

回答两个问题，得到起步单元；最终仍需网格质量检查和收敛性检查。

1　零件是薄壁件吗？
查看完整 6 步流程（文字版）
  1. 明确评价指标整体刚度、模态、局部应力、疲劳寿命还是接触压力。
  2. 判断几何形态薄壁优先用壳单元；实心件能扫掠就用六面体，不能就用二阶四面体。
  3. 处理关注区对关注区单独加密，必要时改用六面体或实体子模型。
  4. 检查网格质量雅可比、翘曲、长宽比、内角等指标要符合求解器和企业规范，不合格的单元不能留在关注区。
  5. 做收敛性检查至少加密一次，关键结果（最大应力、一阶频率等）变化很小时再下结论。
  6. 与试验或理论解对标条件允许时，用应变或模态测试数据验证模型。
经验提示第 5 步最容易被省略。加密后应力仍在持续上升，说明网格或单元类型还不足以表达该处的应力梯度；如果该处是尖角或点载荷引起的奇异点，应力会随加密无限增大，应修正几何或载荷施加方式，而不是继续加密。

八、三个常见误区

  • 把一阶四面体当默认：划分最快，但刚度偏大、应力偏低，结论容易偏于危险。
  • 只看网格是否“好看”：质量指标合格只说明数值上可用，不代表单元类型适合这个问题。
  • 壳与实体直接共节点：交界处弯矩传递失真，局部应力和整体刚度都会受影响。

单元类型选择没有万能答案，但有可重复的判断顺序：先定结果指标，再看几何形态，然后为关注区分配合适的单元，最后用收敛性检查和试验数据确认。对多数工程问题，薄壁件用壳单元；实心件能扫掠就用六面体，不能就用二阶四面体；关注区再局部加密或换用六面体。这是合理的起点。

如果您正在使用 Simcenter Hypermesh 搭建有限元模型，希望评估现有的网格与单元策略，可以了解数领科技的 工程仿真与物理测试 服务。

相关文章

继续阅读