4. 网格划分加速技术

在大型模型求解中，网格划分往往占据前处理阶段的大部分时间，而网格数量与质量又直接决定后续求解的规模与收敛性。加速网格划分并不是单纯把网格“画粗”，而是通过合理的全局与局部控制、合适的划分方法、平滑的尺寸过渡，以及可控的质量约束，在精度、稳定性和计算量之间取得平衡。本节围绕这四条主线展开。

全局控制与局部控制的分工

全局网格控制决定整个模型的基准尺寸，是网格规模的“总闸门”。在 Ansys Mechanical 的 Mesh 分支下，Defaults → Element Size 设定全局单元尺寸，Relevance Center（Coarse/Medium/Fine）与 Relevance（-100～100）共同调节整体疏密倾向。对于大型装配体，建议先用较粗的全局尺寸（例如结构特征尺寸的 1/10～1/20）快速生成一版网格，确认可划分性后再局部加密，而不是一开始就追求细网格。

局部控制则针对关键区域单独加密，常用手段包括：Body Sizing（整体尺寸）、Face Sizing（面尺寸）、Edge Sizing（边尺寸与分段数）、Inflation（边界层）、Contact Sizing（接触面尺寸）以及 Sphere of Influence（球形影响区）。局部尺寸必须小于全局尺寸才会生效，且应遵循“局部加密、远离区域保持粗”的原则。例如一个 500 mm 长的支架，全局尺寸取 10 mm，螺栓孔周边用 Face Sizing 取 2 mm，孔边再叠加 3 层 Inflation，这样既保证孔边应力梯度可捕捉，又不会让整体单元数失控。

扫掠与多区网格：用规则拓扑换速度

对几何规则的长条形、旋转体或可分解为规则块的实体，扫掠（Sweep）与多区（MultiZone）方法能生成六面体主导网格，单元数远少于四面体，且求解精度和收敛性通常更好。扫掠要求实体具有可识别的源面与目标面，且侧面为可映射的扫掠路径；多区方法则自动将实体分解为可扫掠块与自由块，对复杂但仍有规则特征的几何更宽容。

使用扫掠时，关键控制项是 Sweep → Src/Trg Selection（自动或手动指定源/目标面）、Sweep Num Divs（扫掠方向分段数）和 Free Face Mesh Type（源面用 Quad/Tri 或 All Tri）。分段数直接决定厚度方向单元层数，对弯曲问题一般建议不少于 2～3 层。多区网格则通过 Method → MultiZone 指定，可设置 Mapped Mesh 与 Free Mesh Type，并允许对个别块单独控制。若几何存在小圆角、短边或微小面，扫掠常失败，此时可先做几何清理（合并短边、去除装饰性圆角）再重试，这比强行用四面体加密更省时间。

尺寸过渡策略：让疏密变化“平滑”

局部加密后，如果细网格与粗网格直接相邻，会在过渡区产生大量畸形单元，既影响精度又拖慢求解。控制过渡的核心参数是 Growth Rate（增长率），它定义相邻单元尺寸的比值。默认值通常为 1.2，工程上建议控制在 1.1～1.3 之间：过小（如 1.05）会使单元数迅速膨胀，过大（如 1.5 以上）则过渡区单元形状变差。

过渡策略有两种典型做法。其一是“由细到粗”自然过渡：在关键区设定小尺寸，依靠 Growth Rate 逐层放大到全局尺寸，适合应力集中区。其二是“分区控制”：用 Face Sizing 或 Body Sizing 把模型划分为几个尺寸等级，例如孔边 2 mm、过渡区 4 mm、远场 10 mm，每级之间用 Growth Rate 衔接。对于大型装配体，还可以利用 Match Control 保证对称面或周期面上的网格一致，避免接触对两侧网格严重失配导致收敛困难。

网格质量与求解稳定性的平衡

网格质量指标中，最需要关注的是长宽比（Aspect Ratio）、雅可比比率（Jacobian Ratio）、偏斜度（Skewness）和正交质量（Orthogonal Quality）。对结构静力分析，一般要求雅可比比率大于 0.6、偏斜度小于 0.85、长宽比控制在 20 以内；接触面和应力集中区应更严格。质量检查可在 Mesh → Statistics → Mesh Metric 中查看，并利用 Mesh → Insert → Named Selection 定位最差单元。

质量与规模是一对矛盾：加密能提升精度，但单元数增加会显著拉长求解时间，且过度加密反而可能因单元过小、刚度矩阵病态而降低收敛性。实践中可遵循“三步走”：第一步用粗网格试算，确认边界条件、接触和载荷传递无误；第二步在应力梯度大的区域局部加密，观察关键结果（如最大等效应力、最大变形）是否趋于稳定；第三步做一次网格无关性验证，若结果变化小于 5%，即可认为当前网格密度足够。对于大型模型，还可开启 Use Adaptive Sizing 与并行网格划分（Tools → Solve Process Settings 中设置网格划分核数），把划分时间从小时级压缩到分钟级。

一个可复用的参数配置示例

以某大型焊接支架为例，整体尺寸约 800 mm，含 12 个螺栓孔和 4 条焊缝。可参考如下配置：全局 Element Size = 12 mm，Relevance Center = Medium；螺栓孔 Face Sizing = 2.5 mm，孔边 Inflation = 3 层、Growth Rate = 1.2；焊缝区 Body Sizing = 4 mm；规则筋板采用 Sweep，Sweep Num Divs = 4；全局 Growth Rate = 1.2。按此配置，单元数约 80 万～120 万（示例估算，实际取决于几何细节），相比全模型 2 mm 均匀网格可减少 60% 以上单元，同时孔边与焊缝仍保持足够分辨率。划分完成后重点检查 Jacobian Ratio 最小值与 Skewness 最大值，若最差单元出现在非关键区，可接受；若出现在孔边或焊缝，则应缩小该处尺寸或降低 Growth Rate 重新划分。

需要强调的是，上述数值属于教学示例，实际项目应结合几何特征尺寸、材料、载荷类型和求解器类型调整。网格加速的本质不是牺牲精度，而是把计算资源集中到真正影响结果的位置，让每一份单元都“花在刀刃上”。