轴承定位套热锻成形问题的无网格SPH法分析.docxVIP

轴承定位套热锻成形问题的无网格 SPH 法分析 轴承定位套热锻成形问题的无网格 SPH 法分析 杨垠耘，刘泓滨，刘欣玮 （昆明理工大学 机电工程学院，云南 昆明 650500） 摘要：基于配点型无网格理论和紧支试函数加权残量法理论基础，构造了满足一致性条件的核近似函数，在域内及域边界上节点满足平衡方程及力和位移边界条 件的基础上，寻找稳定方案， 并建立了三维弹塑性 SPH 模型。将所建模型运用到某重装企业的轴承定位套锻件生产研究中，对热锻过程中的材料特性参数、锻件尺寸参数等进行了 全面的 SPH 法和 FEM 法模拟，通过与实验数据的对比分析，验证了 SPH 模型的正确性和 SPH 法在分析金属塑性成形过 程的优越性，对实际生产具有一定的参考价值。 关键词： 锻造成形；轴承定位套； SPH 法； FEM 法；模拟 20 世纪 70 年代，随着计算机技术的发展和软件工程的兴起及有限单元 法（ Finite Element Method ， FEM ）非线性理论的不断成熟，模锻过程的计算机模拟仿真技术得到了空前发展，有限元法 也成为了分析金属塑性成形最有效的方法。但在处理大变形问题时，由于有限元法依赖于网格而使其计算受阻，需多次 进行网格划分，大大降低了计算效率和分析准确性 [1]。 为 解决此问题，很多国内外专家学者提出了用不依赖于网格的 无网格法进行金属塑性成形分析计算 [2] 。经过近三十年的发 展，无网格法已日 成熟，并大量运用在金属塑性成形 中。其中，由于光滑粒子 力学法（ Smooth Particle Hydrodynamic ，SPH）是一种 Lagrange 方法 [3] ， 就能 更好地避免欧拉网格与材料的界面 ，所以在流体力学分 析及高速撞 等方面得到了广泛的 用，近些年来在金属塑 性成形 域也已逐 得到 用。 Bonet 等将其 用于理想状 下的平面 粗和 称 造 研究 [4] ；国内李 生等将其用于平面 下微 材料的金属塑性 形 研 究[5] 。本文以某企 生 中的某型号 承定位套 例，建立 FEM 及 SPH 分析模型，分 运用两种模型 承定位套 成形 程 行模 仿真， 通 其仿真 果与 数 据的 比分析， SPH 模型的正确性和 越性。 1 SPH 基本原理 光滑 点流体 力学 （SPH）方法是采用核近似的配点型无网格方法，它已 被耦合到 示 力有限元分析程 序 PRONTO 和 LS-DYNA 中，成功求解高速碰撞、 爆炸等 ，充分 了有限元和 SPH 法各自的 [6] 。 在 SPH 中，求解域用 N 个 点 xI （I=1 ，2，?， N ）来离散。令运 方程在各 点 足， 并将 力 空 坐 的散度 σ ij ，得 式 中： σJ—— xI 的 力 量。 SPH 法中多采用差分法在 域中 行 分，利用式（ 1）可以由 刻的 力 量 σIn 求出 刻 tn 的加速度 VIn ，然后利用蛙跳（ leap-frog ）差分格式得到 刻 tn+1/2 的速度 vIn+1/2 和 刻 tn+1 的坐 xIn+1 ： 式中 步 tn+1/2=tn+1-tn 。 里加速度和坐 是在 刻 tn+1 算的，而速度 是在 刻 tn+1/2 算的。 刻 tn+1/2 的速度 vIn+1/2 可求得此 刻的速度梯度， 而求得 刻 的 率： ρ ˙+ρ ▽·v=0 利用中心差分格式， 得： 再 由本构方程可得到 刻 tn+1 的 力： 1.1 核函数 由 δ的 性 ，有：  如果用任意 ω（ x-ˉ）来代替 δ（ x-x ˉ）， 式（ 7）不再 格成立，函数的近似 ：  式中函数 ω（ x- ） 即 核函数。 由于在 SPH 法中，核函数具有 性、 称 性、 一性、半正定性、光滑性等特点，所以常采用 B 条 函数作 SPH 中常 的核函数 [7] ： 式中： r=‖ x-ˉ‖ /h， d 是空 数， 一化常数 G 等于（三 ） 。 1.2 定 方案 在 SPH 中，采用 特 的 定方案，以避免方程的病 。大多数 定方案都是在域内 点平衡方程、 界 Γ t 上 点 界条件和 界 Γu 上 点位移 界条件的矩 形 式上增加控制方程残差 。 未知函数在域内 行泰勒展 开，保留高 的有限点法更 理想： 1.3 形函数 在 SPH 中，将流体剖分成数个小体 元， 各体 元的 心 x1，x2，?， xN 作 点。 把各体 元的 量 mJ 予 点 xJ，便可以将 xJ 粒子或 点（ Particle）。各 点所 的体 VJ=mJ/ ρJ，其中 ρ J=ρ（xJ） xJ 的密度。因此： 可 ，SPH 的形函数其 就是