title: Workbench 三维锚杆拉拔模拟操作流程
date: 2024-01-15
tags:

Workbench 三维锚杆拉拔模拟操作流程

[!tip] 前提知识
本文基于 [[锚杆仿真要点]] 中所述的二维分析逻辑扩展而来，建议先阅读该文档了解基础概念。

进行三维（3D）抗浮锚杆拉拔模拟时，整体的分析逻辑（材料本构、两步加载、位移控制）与二维模型一致，但计算资源消耗会呈指数级增加，且网格划分与接触收敛的难度会大幅提升。

三维模型通常用于分析二维无法涵盖的复杂工况，例如：群锚效应、倾斜拉拔（抵抗风荷载等产生的斜向扬压力）、非均匀土层分布、或复杂的上部结构连接。

[!warning] 计算成本警告
3D 模型的计算时间可能是 2D 模型的 10-100 倍，请确保硬件配置足够（建议 16 核以上 CPU，32GB+ 内存）。

以下是针对 ANSYS Workbench模拟的核心操作指南与避坑要点：

三维拉拔目录

在 3D 非线性分析中，几何不能直接“画完就导”，必须经过精心切分。

利用对称性（极度重要）：

[!important] 核心优化策略
如果是单根垂直锚杆且土体均匀，绝对不要建立完整的 360° 圆柱模型。
- 强烈建议建立 1/4 或 1/2 对称模型。这不仅能节省 75% 的网格量和计算时间，还能消除刚体位移，帮助模型更容易收敛。
几何切块 (Geometry Slicing)：
- 为了后续能划分出高质量的六面体网格，必须在 SpaceClaim 中对土体进行切分。通常采用“同心圆柱+放射状”的切法，将大块土体切分为多个可以被“扫掠 (Sweep)”的规则块。
拓扑共享设置：
- 土体内部的各个切块之间，设置为 Share Topology = Share（共享节点，保证应力连续）。
- 锚杆/注浆体 与土体之间，必须设置为 Share Topology = None（不共享节点），因为我们需要在这里手动建立摩擦接触面，允许它们发生相对滑动。

在涉及岩土本构（如 Mohr-Coulomb）和摩擦接触的 3D 仿真中，网格质量直接决定了计算能否算下去。

摒弃四面体，拥抱六面体 (Hexahedral Mesh)：

[!danger] 网格类型警告
自动生成的四面体网格（Tetrahedrons）在处理不可压缩材料的大变形时极易发生”体积自锁 (Volumetric Locking)”，导致土体异常坚硬，计算结果严重失真且极难收敛。
- 必须使用 MultiZone（多区划分） 或 Sweep（扫掠） 方法，强制生成纯六面体单元（底层对应 APDL 中的 Solid185 或 Solid186 单元）。
接触面网格匹配：
- 插入 Sizing（尺寸控制），确保注浆体外表面和土体内孔表面的网格尺寸尽可能一致。网格尺寸差异过大会导致接触面法向力计算振荡。

三维曲面的接触比二维线接触复杂得多，容易发生初始穿透。

建立接触： 选择注浆体外圆柱面为 Contact，土体内孔面为 Target。类型设定为 Frictional（摩擦）。
算法 (Formulation)： 改为 Augmented Lagrange（增广拉格朗日）。这种算法对法向刚度的敏感度较低，更适合大变形滑移。
法向刚度 (Normal Stiffness)： 初始默认因子为 1.0，如果在拉拔瞬间直接报错不收敛，通常是因为接触刚度太大。将其改为手动输入（Manual），设定为 0.1 甚至 0.01。

[!tip] 收敛技巧
如果遇到不收敛问题，优先调小法向刚度因子，这是最容易解决的问题。
初始干涉处理 (Interface Treatment)： 设置为 Adjust to Touch（调整以接触）。这能自动消除由于网格离散化带来的微小几何间隙或穿透，防止计算第一步就因为找不到接触而发散。

对称边界： 如果你使用了 1/4 模型，必须在切开的两个平面上施加 Frictionless Support（无摩擦约束）。这在力学上完美等效于对称边界（限制了法向位移，允许面内滑动）。
侧面与底面：
- 土体底面：Fixed Support（全固定）。
- 土体外侧圆柱面：施加 Cylindrical Support（圆柱坐标约束），仅固定径向 (Radial)，释放轴向 (Axial) 和环向 (Tangential)。
两步加载法：
- 同样需要分为两步：Step 1 开启重力建立初始地应力，Step 2 施加向上的强制位移 (Displacement) 进行拉拔。

三维非线性岩土接触问题的计算矩阵非常庞大。

开启大变形： Large Deflection = ON。
增加子步数： 在 Step 2 中，初始子步数 (Initial Substeps) 建议设为 100 左右，最大子步数 (Max Substeps) 设为 1000。
并行计算 (HPC)： 在 Workbench 主界面 -> Tools -> Options -> Solution Process 中，务必将核心数 (Default Number of Cores) 调到你电脑的最大物理核心数（例如 8核或 16核）。如果仅用单核计算，3D 模型可能需要几天时间才能跑完。

[!warning] 性能优化
切勿忽视 HPC 设置！单核计算 3D 非线性问题可能需要几天时间。