管道应力分析，看这一篇就够了——ASME B31.3 管道应力分析工具书

在工程仿真中，我们经常遇到计算域边界随时间变化的问题：内燃机活塞运动、阀门开闭、飞行器投放物分离、流体与结构相互作用……这些问题的共同特点是——网格必须随之运动或变形。Ansys Fluent 的动网格（Dynamic Mesh）模型正是为此而生。本文基于 Fluent 2025R2 官方用户手册，系统梳理动网格的完整设置流程、三种核心更新方法（平滑、动态层铺、重构）及其详细参数、内燃缸专用功能、六自由度（6 DOF）求解器、隐式更新、接触检测、事件定义以及稳态动网格应用。一、动网格基本设置流程（共13步）无论是瞬态还是稳态问题，设置动网格的标准步骤如下。注意：此流程仅包含动网格模型本身必需的步骤，其他模型（湍流、传热、多相流等）按常规设置即可。启用瞬态或稳态流动Setup → General 瞬态问题：选择Transient；稳态问题：参见后文“稳态动网格应用”章节。 准备相邻运动区域的网格拓扑 若单元区域相对运动，边界必须拥有独立的面区域且不共享节点。可通过断开连接或分别生成网格再合并。 设置单元区域条件和边界条件Setup → Cell Zone Conditions ，Setup → Boundary Conditions动态区域运动属性会自动设定壁面速度。若壁面实际不运动则需声明动态区域并排除运动影响。 启用动网格模型并设置参数Setup → Dynamic Mesh → Dynamic Mesh 勾选 Dynamic Mesh，选择 Options（如 In‑Cylinder、Six DOF、Implicit Update）和 Mesh Methods（平滑、层铺、重构）。 创建动态区域Setup → Dynamic Mesh → Create/Edit... 为运动/变形区域指定运动类型（刚体、变形、用户自定义等）。 （可选）显示预设运动Setup → Dynamic Mesh → Display Zone Motion... 定义事件（仅瞬态）Setup → Dynamic Mesh → Events... 控制特定时刻的动作（创建/删除滑移界面、激活区域等）。 保存 Case 和 Data 文件File → Write → Case & Data ⚠️ 重要：每次保存 Data 时必须同时保存 Case 文件，因为网格位置信息存储在 Case 文件中。 预览网格运动Setup → Dynamic Mesh → Preview Mesh Motion... 指定压力-速度耦合方案 瞬态计算推荐使用PISO 算法。 设置自动保存Solution → Calculation Activities → Autosave （可选）设置图形动画Calculation Activities Task Page 使用“动网格库朗数”指导求解 场变量Moving Mesh Courant Number指示网格运动在一个时间步内扫过的单元数量。 二、动网格更新方法总览Fluent 提供三类网格更新方法，可组合使用：方法 特点 适用网格类型 平滑（S moothing） 节点移动，拓扑不变 任何类型 动态层铺（Layering） 增加/删除网格层 棱柱形（六面体/楔形/四边形） 重构（Remeshing） 局部重新划分网格 三角形/四面体为主 ⚠️ 注意：平滑、层铺和重构与网格自适应不兼容。三、平滑方法（S moothing Methods）平滑方法通过移动内部节点来“吸收”边界运动，节点数量和连接关系保持不变。启用方式：在 Dynamic Mesh Task Page 中勾选S moothing，然后点击Settings...选择具体方法。3.1 基于扩散的平滑（Diffusion‑Based S moothing）控制方程为扩散方程：∇·(γ ∇ u̅) = 0其中 u̅ 为网格位移速度，γ 为扩散系数。边界条件来自用户指定的边界运动，在变形边界上网格运动与边界相切（法向速度分量为零）。扩散系数的两种形式基于边界距离（Boundary Distance）：γ = 1/dα，α≥0。α=0 均匀扩散；α越大，近边界区域保留越好，远处吸收更多运动。💡 对于旋转运动，推荐以 α=1.5 作为起始值。 基于单元体积（Cell Volume）：γ = 1/Vα，V为归一化单元体积。α越大，大单元吸收更多运动，小单元质量保持更好。 数值求解方法对于无非多面体单元或悬挂节点的网格：采用有限元离散，直接获得各节点的位移速度。对于包含多面体或悬挂节点的网格：采用有限体积法，单元中心的解通过反距离加权插值到节点。可通过 TUI 命令强制使用有限体积法：/define/dynamic-mesh/controls/s moothing-parameters/diffusion-fvm? yes线性求解器选项（AMG Stabilization）：CG（共轭梯度法）：默认，利用矩阵对称性，速度快；若出现负体积可增加迭代次数至200~500。 GMRES（广义最小残差法）：更稳健，尤其适用于高长宽比网格，但内存和计算时间更高。 BCGSTAB（双共轭梯度稳定法）：类似 CG，发散时回退到 GMRES。 适用性：支持所有单元类型，但对多面体不如 RBF 稳健；计算成本高于弹簧平滑，但网格质量更好；平移运动处理优于旋转运动（旋转运动推荐 RBF）；与边界层平滑方法、面区域重构方法不兼容；不支持变形共形周期边界或固体移动壁面。3.2 基于弹簧的平滑（Spring‑Based S moothing）将网格边理想化为相互连接的弹簧网络。边界节点位移产生与位移成正比的力。F̅ᵢ = Σ kⱼ(Δxⱼ-Δxᵢ)kⱼ = k_fac / √|xᵢ-xⱼ|k_fac 为弹簧常数因子，范围 0~1。平衡时净力为零，得到迭代方程：Δx̅ᵢ^{m+1} = (Σ kⱼ Δx̅ⱼ^{m}) / (Σ kⱼ)收敛判据：Δx_rms^{m} / Δx_rms^{1} < Convergence Tolerance。弹簧常数因子的影响：默认值1有一定阻尼，内节点位移可能超出运动边界（产生畸变单元）；值0无阻尼，内节点位移更好地跟随边界，保持原始分布。适用性：适用于运动以平动为主、主要垂直于边界的场景；非四面体/三角形单元区域需谨慎，多面体单元尤其容易产生高畸变 → 推荐 RBF；不支持变形共形周期边界或固体移动壁面。3.3 基于线弹性固体的平滑（Linearly Elastic Solid Based S moothing）将网格视为线弹性固体，控制方程：∇·σ=0，σ=2με+λ tr(ε)I，ε=½(∇d̅+(∇d̅)ᵀ)。泊松比 ν = λ/(2(λ+μ))，范围 -1.0~0.5，仅支持常数材料属性。边界条件：若指定了变形边界的几何（平面、圆柱等），则法向变形为零，切向应力为零；若几何为unspecified（仅此方法和 RBF 支持），则边界也可在法向变形，所有方向应力为零。适用性：计算成本比扩散平滑更高，但对某些网格和运动能更好地保持质量；支持 2D 三角形/四边形，3D 四面体/六面体/楔形/金字塔；不支持多面体单元、悬挂节点；与边界层平滑、面区域重构不兼容；不支持变形共形周期边界或固体移动壁面。3.4 径向基函数平滑（Radial Basis Function S moothing, RBF）RBF 求解一个位移场，该位移场在指定运动的边界上满足位移值，然后插值到所有其他节点。主要优势：特别适合旋转运动，可吸收比扩散/弹簧平滑更多的旋转量。支持悬挂节点和多面体单元。扩散平滑对多面体不稳健，RBF 则很稳健。唯一支持以下场景的平滑方法：平滑区域内包含变形共形周期边界；平滑区域是带有移动壁面的固体（如烧蚀、侵蚀/堆积问题）。 变形边界几何选项：可指定几何（平面、圆柱等）让节点沿几何切向滑动；也可选择unspecified，节点可任意方向移动。对于非共形界面，若两侧界面均为 unspecified 的变形动态区域，Fluent 会保证两个界面保持连接（仅 RBF 支持，且要求界面两侧为同一单元区域）。高级设置：Relative Convergence Tolerance（控制位移场求解精度）；Verbosity=1 在控制台打印绝对容差；Local S moothing：仅对运动边界附近的单元进行平滑，远离边界的网格不受影响，可大幅减少计算量。若启用局部平滑，通常还需要同时启用重构，因为运动被吸收在更小的区域中。边界层单元会在局部平滑之前独立平滑，以保持形状。注意：RBF 在分布式并行计算中的扩展性不如其他方法；与边界层平滑方法不兼容。3.5 从参考位置平滑（S moothing from a Reference Position）对于扩散、线弹性固体或 RBF 平滑，可以指定从固定的参考位置进行平滑。这对于周期性或准周期性运动（如涡轮机械叶片颤振）非常有用，可使每个循环的网格质量保持一致。启用方式：在 Mesh S moothing Parameters 对话框中勾选S moothing From Reference Position。注意：若同时启用 Layering 或 Remeshing，或激活了自动自适应判据，则该选项不可用；不要以 legacy 格式保存 case 文件，否则重启后参考位置会不同。3.6 拉普拉斯平滑（Laplacian S moothing）最简单的方法：将每个节点移动到其相邻节点的几何中心。x̅ᵢ^{m} = (Σ xⱼ^{m})/n_i， x̅ᵢ^{m+1} = x̅ᵢ^{m}(1-β) + x̅ᵢ^{m}ββ 为拉普拉斯节点松弛因子。Fluent 仅在改善网格质量（畸变率降低）时才更新节点位置。可用于变形边界（3D 三角形单元，2D 线性单元），也用于 2.5D 重构方法中。3.7 边界层平滑方法（Boundary Layer S moothing Method）用于使边界层网格随相邻运动/变形面区域一起变形。启用方式：在 Dynamic Mesh Zones 对话框中为面区域勾选Deform Adjacent Boundary Layer with Zone。工作原理：边界层中每个单元的节点坐标与底层面对应节点使用相同的位移向量。适用性：支持 3D 楔形/六面体边界层，2D 四边形边界层；仅与弹簧平滑兼容，不可与扩散、线弹性固体或 RBF 平滑同时使用；主要用于平动，大幅旋转会导致边界层单元畸变。四、动态层铺（Dynamic Layering）适用于棱柱形网格（3D 多面体/六面体/楔形，2D 四边形）。根据运动边界相邻层的高度来添加或删除网格层。层分裂条件（膨胀时）：hmin> (1+αs)·hideal层合并条件（压缩时）：hmin< αc·hideal其中 αs为分裂因子，αc为合并因子，hideal为理想单元高度（可常数或随时间和曲轴角变化）。分裂选项：高度法（Height Based）：产生高度为 hideal的新层 + 高度为 h - hideal的剩余层。 比率法（Ratio Based）：局部新/旧层高比恒为 αs。 适用条件：运动边界相邻单元必须是棱柱形（即使单元区域内混合其他形状也可以）；网格层必须由单侧面区域完全包围（除非使用滑移界面）；不能与自适应同时使用（几乎在所有情况下）。五、重构（Remeshing）当边界位移相对于局部单元尺寸过大时，仅靠平滑会导致网格质量恶化甚至负体积。Fluent 会聚合不满足畸变率或尺寸判据的单元，并局部重新划分网格。5.1 统一重构（Unified Remeshing）—— 推荐默认启用 Unified Remeshing 后，算法结合多种重构方法的优点，并默认尝试保持初始网格尺寸分布，简化设置，提高稳健性（尤其适合并行计算）。高级设置（点击 Advanced...）：Maximum Cell Skewness：3D 默认 0.9，2D 默认 0.7 Maximum Face Skewness：用于边界上的面重构 Retain Size Distribution：默认启用，若禁用则使用初始网格中的最小/最大长度尺度作为标记依据 尺寸控制（Sizing Controls）和棱柱控制（Pris m Controls）可通过 TUI 命令创建，用于精细控制网格尺寸分布或边界层参数（第一层高度、增长率、层数等）。在 Watertight Geometry 网格工作流中创建的尺寸/棱柱控制会自动导入统一重构。5.2 基于方法的重构（Methods‑Based Remeshing）若选择 Methods‑Based Remeshing，可单独启用以下子方法：局部单元重构（Local Cell Remeshing）：标记单元的条件（畸变率>阈值，或长度尺度超出最小/最大范围，或高度不满足）。若局部单元重构无法充分降低最大单元畸变率，则自动触发单元区域重构。尺寸重构间隔（Size Remeshing Interval）控制基于尺寸和高度的标记频率。 局部面重构（Local Face Remeshing）：仅适用于3D。对变形边界上的三角形面进行标记和重构，也可重构边界层中的楔形单元。需在 Mesh Method Settings 中启用 Local Face，并在 Dynamic Mesh Zones 中为每个边界启用 Local 选项。 单元区域重构（Cell Zone Remeshing）：对整个单元区域完全重新划分网格，可选择同时重构所有相邻变形动态面区域。支持2D三角形、3D四面体以及3D边界层中的楔形单元。边界层检测自动进行（可手动通过 TUI 指定第一层高度、增长率和层数）。限制：不能重构有悬挂节点的区域。 面区域重构（Face Region Remeshing）：主要用于内燃缸配置，对运动面区域（如活塞顶）与相邻变形面区域交界处的面进行重构。当膨胀层高度达到最大长度尺度时，一层变为两层；压缩到最小长度尺度时，两层合并为一层。推荐：若运动范围很大（超过缸径一半），最好将网格体积分解，使边界层完全位于刚体运动区域内。 2.5D 表面重构（2.5D Surface Remeshing）：专门用于可拉伸的棱柱网格（如齿轮泵）。本质上是一个2D三角形网格沿法向拉伸形成3D网格。三角形表面网格在一侧重构和平滑，变化沿拉伸方向传递到另一侧。运动面区域采用刚体运动，拉伸面设置为变形区域（启用重构+平滑），对面设置为仅平滑。只支持 Laplacian 平滑，不与其他平滑方法兼容。拉伸必须沿直线且垂直于变形区域，横截面沿拉伸方向恒定。并行时推荐使用垂直于拉伸面的分区方法（如 Cartesian-Y 或 Cartesian-Z）。 六、尺寸函数（Sizing Function）用于局部重构时，基于边界单元尺寸生成背景网格，从而控制网格尺寸的空间分布。背景网格：Fluent 在每个流体区域周围画一个大约两倍于区域尺寸的包围盒，根据最短特征长度和用户指定的 Resolution 细分，创建背景网格。在每个背景网格顶点 I 计算尺寸函数值：SF_I = (Σ D_J^{-1}·Δs_J) / (Σ D_J^{-1})然后通过拉普拉斯平滑，并插值到任意单元形心 P：size_p = size_b × γγ = 1 + α·d_b^{1+2β} (α>0) 或 γ = 1 + α·d_b^{1/(1-β)} (α<0)参数：α（Variation）：-1 ~ ∞，控制远离边界时单元尺寸的变化幅度。α=0.5 表示最大尺寸为边界尺寸的1.5倍。 β（Rate）：-0.99 ~ 0.99，控制尺寸变化速率。正值变化更平缓，负值变化更剧烈。 d_b：归一化距离（最大值为1）。 Resolution：背景网格分辨率（2D默认3，3D默认1）。 启用尺寸函数后，单元被标记为需要重构的条件是：size ∉ [4/5·γ·size_b, 5/4·γ·size_b]。七、体积网格更新过程总结Fluent 根据启用的模型和单元形状自动决定每个区域使用哪种更新方法：四面体/三角形区域 → 平滑 + 局部重构棱柱形区域（六面体/楔形） → 动态层铺拉伸棱柱（2.5D） → 2.5D表面重构注意：重构可能导致并行计算中各分区网格略有不同，从而引起非常微小的解差异。若模型包含大量壳传导区域，每次网格更新会删除并重建所有壳，非常耗时。八、时间推进的瞬态考虑二阶时间格式可与重构、层铺和平滑一起使用。层铺事件和网格平滑期间保持二阶精度；重构事件发生时，时间推进精度自动降为一阶。若重构非常频繁（每个时间步都发生），实际就是一阶格式。可通过 TUI 命令define/dynamic-mesh/transient-settings/verbosity查看当前使用的时间格式。九、特征检测（Feature Detection）在 3D 中，重构和/或平滑可能导致角部被倒角。若某些角部（特定角度范围）对流动很重要，可启用特征检测来保留它们。在 Dynamic Mesh Zones 对话框的 Geometry Definition 选项卡中，勾选 Feature Detection 并设置 Feature Angle：180°（默认）不保留任何角部，0°保留所有角部。适用性：仅与面区域重构配合使用，特征边上的节点不允许平滑。十、内燃缸专用设置（In‑Cylinder）启用 In‑Cylinder 选项后，需要设置：曲轴转速（Crank Shaft Speed）、起始曲轴角（Starting Crank Angle）、曲轴周期（Crank Period，通常720°四冲程或360°二冲程）、曲轴角步长（Crank Angle Step Size，默认0.5°）。内置活塞运动函数 piston‑full：活塞位置作为曲轴角的函数ps= √[(rc+L)² - x_offset²] - rc·cosθc- √[L² - (rc·sinθc+x_offset)²]其中 rc曲柄半径，L 连杆长度，x_offset 活塞销偏移（有正负），θc= θs+ t·Ωshaft。限制参数：Piston Stroke Cutoff（活塞行程截止值，控制层铺何时开始）；Minimum Valve Lift（最小气门升程，防止气门与阀座接触时单元退化为扁平）。气门升程剖面：使用 Profile 文件定义，必须包含 angle 和 lift 字段。示例：((ex-valve 5 point) (angle 0 180 270 360 720) (lift 0.05 0.05 1.8 0.05 0.05))滑移接口事件：气门关闭时（升程<最小值）删除滑移接口，并将接口区域转换为壁面；气门打开时重新创建滑移接口。这些通过 Dynamic Mesh Events 实现。涡流/滚流输出：可启用 Write In‑Cylinder Output，设置输出频率、涡流中心计算方法（重心/固定/UDF）、涡流轴、滚流轴、交叉滚流轴等。十一、六自由度（Six DOF）求解器计算作用在刚体上的气动力和力矩（通过压力与剪切应力积分），加上重力、弹簧力、推力等外部载荷，求解刚体运动。适用于外挂物投放、弹射分离等。启用：Dynamic Mesh Task Page → Options → Six DOF设置刚体属性：在 Six DOF Properties 对话框中为每个运动物体创建属性集，必须定义：Mass（质量）、Inertia Tensor（惯性张量，相对于重心和刚体方向）。可选用：One DOF Translation（单自由度平动，可定义弹簧刚度、预紧力、运动范围）、One DOF Rotation（单自由度转动，需定义转动惯量）。运动历史输出：勾选 Write Motion History，为每个运动物体生成历史文件。配合隐式更新：若六自由度运动依赖于流场，建议启用隐式网格更新。若仍不稳定，可在边界区域的 Solver Options 选项卡中启用 Solution Stabilization。十二、隐式网格更新（Implicit Update）对于网格运动依赖于流场的应用（如六自由度、流固耦合），隐式更新允许在一个时间步内多次更新网格，基于收敛中的流场解逐步调整网格位置，实现更强的耦合，提高鲁棒性。启用：Dynamic Mesh Task Page → Options → Implicit Update设置（Options 对话框 → Implicit Update 选项卡）：Update Interval：一个时间步内网格更新的迭代频率 Motion Relaxation ω（0~1）：节点位置松弛，xk= ω·xcomputed,k+ (1-ω)·xk-1 限制：不能与密度基显式求解器、稳态问题、内燃缸同时使用；使用 UDF 计算运动时，UDF 必须在每次调用时使用当前流场，不能存储之前的信息。十三、接触检测与流控（Contact Detection）用于检测运动表面是否与其他表面发生接触（距离小于阈值）。接触信息可在 UDF 中用于约束运动，或控制接触区域的流动。启用：Dynamic Mesh Task Page → Options → Contact Detection设置：从 Face Zones 列表中选择参与检测的面（仅壁面和 Rigid Body/User‑Defined 动态区域）；Proximity Threshold：触发接触的距离阈值；UDF：接触时调用的函数（可选）。流控方法：Contact Zones（接触区法）：将接触区域内的单元分离到一个新的单元区域，可设置多孔介质阻力等。完全兼容所有求解器，但不能实现完全流动阻断；若接触区域快速变化，计算开销较大。 Contact Marks（接触标记法）：在接触区域边界施加零质量通量数值边界条件，完全阻断流动。计算效率高，无泄漏，可减少窄缝中的网格需求。限制：仅压力基求解器；不支持 Green‑Gauss Node Based 梯度、QUICK 格式、NITA 格式；不推荐 outflow 边界。可启用 Solution Stabilization 改善稳定性。 显示接触标记：通过 TUIdefine/dynamic-mesh/controls/contact-parameters/render-contactcells?启用，然后在 Contours 中选择 Cell Info... → Contact Cell Mark 显示。十四、动态网格事件（Dynamic Mesh Events）用于在瞬态模拟的特定时刻（时间或曲轴角）执行预定动作。内燃缸模拟中尤其重要（气门开闭、滑移接口创建/删除）。设置：Dynamic Mesh Events 对话框。增加 Number of Events，勾选事件，输入名称，指定触发时间（At Time）或曲轴角（Act Crank Angle），点击 Define... 选择事件类型并设置参数。点击 Preview... 可播放事件以验证。可用事件类型：Change Zone Type、Copy Zone BC、Activate/Deactivate Cell Zone、Create/Delete Sliding Interface、Change Motion Attribute、Change Time Step Size、Change Under‑Relaxation Factors、Insert/Remove Boundary Zone Layer、Insert/Remove Interior Zone Layer、Insert/Remove Cell Layer、Execute Command、Replace Mesh、Inert EGR Reset、Diesel Unsteady Flamelet Reset 等。内燃缸事件循环：事件在满足 θ_event = θ_c ± n·θ_period 时重复执行。十五、稳态动网格应用动网格也可用于稳态问题，例如：在不同阀门位置达到稳态后进行比较，或流固耦合界面达到稳态后更新网格。与瞬态的主要区别：使用CG_MOTION UDF指定运动（dtime 默认为1，若需位移1 mm，设速度为1e-3 m/s）。不能使用瞬态 profile。可通过 Journal 文件或 Execute Commands 手动触发网格更新。可启用额外的重构（基于畸变率，在边界移动后执行）。不可用：In‑Cylinder、Implicit Update、Dynamic Mesh Events。 更新方式：GUI：Solution → Run Calculation → Update Dynamic Mesh...；TUI：solve → mesh-motion。建议：六自由度稳态问题可能需要调整伪时间步长，确保在更新前稳态解已较好收敛。十六、重要注意事项汇总每个时间步必须同时保存 Case 和 Data 文件。重构可能导致并行计算中各分区微小差异。重构事件发生时二阶时间精度自动降为一阶；若重构非常频繁，建议直接用一阶格式。计算前务必使用 Preview Mesh Motion 检查是否出现负体积或严重畸变。周期性运动启用 Smoothing From Reference Position 可保持循环间网格质量一致。边界层平滑方法主要用于平动，大幅旋转会使边界层单元畸变。需要完全阻断流动时优先选择 Contact Marks 方法。多面体网格推荐 RBF 平滑，不要用弹簧或扩散（除非已验证可行）。层铺与自适应几乎在所有情况下都不能同时使用。设置自动保存频率，避免因意外丢失长时间计算结果。来源：仿真与工程