基于水平集的软性磨粒流流场测试方法技术

一种基于水平集的软性磨粒流流场测试方法，包括以下步骤：１）通过与流体Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程结合，建立软性磨粒流的水平集函数的运动方程；２）对空间作离散化和投影；３）对时间作离散化；４）水平集函数与流体体积模型函数的结合进行迭代；５）距离函数的重构；６）对软性磨粒流流场的测试：利用水平集方法的有效捕捉界面发生拓扑变化的原理，对流场特性进行较好地分析，找到最佳的加工状态。本发明专利技术提供一种提高计算精度、准确得到最佳加工状态的基于水平集的软性磨粒流流场测试方法。

【技术实现步骤摘要】

本专利技术涉及磨粒流加工领域，尤其是一种软性磨粒流的结构化表面无工具加工中的流场测试方法。
技术介绍
目前，模具领域对结构化表面的粗糙度要求越来越高，例如精密注塑模具，不但对注塑件表面对应的模具面的粗糙度有极高要求，甚至对于注塑流道的粗糙度也有很高要求。对于尺寸小或几何形态特殊的结构化表面难以使用工具进行接触式光整加工，如抛光或研磨，无论自动化或手工加工都是如此，这一问题目前尚无有效方法解决，且缺乏深入的研究。鉴于模具应用量大面广，其制造水平直接影响制造业的整体水平，因此对模具结构化表面光整加工方法的研究，具有重要的现实意义。 所谓的软性磨粒流是一种液-固两相磨粒流，具有弱黏性或无黏性，因此具有更好的流动特性并可实现湍流流动，基于软性磨粒流实现模具结构化表面无工具精密加工的原理是(1)通过在被加工的结构化表面附近配置约束模块，构成磨粒流约束流道，使被加工表面成为流道壁面的一部分；(2)以约束流道内流动的软性磨粒流替代加工工具实现对被加工表面的光整加工；(3)软性磨粒流的有效加工是在湍流状态下进行，它不是通过射流的形式强力冲击被加工表面，而是利用磨粒的微力微量切削的频繁作用实现表面的逐步光整，湍流流场中的磨粒运动的随机性实现了表面纹理无序化，直至实现结构化表面无工具镜面级加工。
技术实现思路
为了克服已有软性磨粒流流场的计算精度低、无法准确得到最佳加工状态的不足，本专利技术提供一种提高计算精度、准确得到最佳加工状态的。 本专利技术解决其技术问题所采用的技术方案是 一种，所述测试方法包括以下步骤 1)通过与流体Navier-Stokes方程结合，建立软性磨粒流的水平集函数的运动方程； 2)对空间作离散化和投影用一组有限个离散的点来代替原来的连续空间，同时经投影作用使得有散度的转化为无散度的量； 3)对时间作离散化； 4)水平集函数与流体体积模型(Volume of fluid Method，VOF)函数的结合进行迭代； 5)距离函数的重构重构就是在已知函数分段平均值的情况下，构造满足一定精度下的原函数分布，由于数值方法的内在效应，经过若干演算时间后，原距离函数已不再满足要求，通过求解初始问题的稳定解来实现； 6)对软性磨粒流流场的测试利用水平集方法的有效捕捉界面发生拓扑变化的原理，对流场特性进行较好地分析，找到最佳的加工状态。 进一步，所述软性磨粒流是松散磨粒与流体混合，构成液-固两相或气-固两相磨粒流。 再进一步，所述步骤1)中，利用连续介质表面Φ(x，t)来描述两相流的界面边界，计算域内其他计算点处Φ值则表示该点到相界面的距离，其定义为Φ＝Const，即t时刻时，边界零水平集Φ Γ＝{x|Φ(x，t)＝0}(1) 在磨粒运动过程中Φ遵循的守恒方程 式(2)中u为流场中的速度矢量； 为Φ的梯度 在此设Φ＜0为磨粒流区(l1)，Φ＞0为流体区(l2)，则有 两相流的边界区域可表示为 式(4)中 x∈Γ；σ为边界张力系数；Γ为边界区域；p为流体压强；μ为流体粘度；D/Dt为随体导数； 对于边界的速度单位矢量方向 及曲率κ可表示为 及 为了保持Φ的距离函数性质，利用SUSSMAN M等提出的方法，对方程(2)进行时间离散的迭代重新初始化，其解为 初始化条件d(x，0)＝Φ(x)； 式(6)中 当sign(0)＝0时，d(x，τ)与Φ(x，t)表示的零水平集等价；τ为重新初始化过程中迭代的虚拟时间，d为距离函数；根据方程(6)来简单地表示出系统的稳态d(x，τsteady)，即 更进一步，所述步骤4)中，根据水平集连续介质表面Φ的公式，定义体积模型F 式中 Ω为网格划分单元。 所述步骤6)中，以SIMPLEC算法求解，模拟分析流体特性。 本专利技术的技术构思为水平集(Level set Method，LSM)具有基于边界的处理结构拓扑变化能力和较好的鲁棒性，Burger M、Wang M Y等人使用响应泛函的灵敏度信息构造了水平集进化所需的速度场，并采用LSM进行了拓扑优化设计。近年来，LSM在图像分割、目标跟踪等图像处理技术方面应用广泛，CASELLES V、PARAGIOS N等人用该方法求解由活动轮廓模型得到的偏微分方程，从而实现对图像的分割。随着计算流体力学(Computational Fluid Dynamics，CFD)的快速发展，LSM逐渐扩展到多相流研究领域，其主要求解方法为投影法与压力耦合方程的半隐相容(Semi-Implicit Method for Pressure Linked Equations Consistent，SIMPLEC)算法。 所述的软性磨粒流是松散磨粒与流体混合，可构成液-固两相或气-固两相磨粒流，磨粒流的流体性质决定其可变化无形且无孔不入，因此，基于磨粒流形成了一些表面加工方法。软性磨粒流，具有弱黏性或无黏性，因此具有更好的流动特性并可实现湍流流动，在处于湍流状态下，可利用磨粒的微量切削的频繁作用，以约束流道内流动的软性磨粒流替代加工工具，实现对表面的无工具的光整加工。 所述的水平集方法把低维空间上的函数转化为高位空间分析，将两流体分界面用水平集函数(定义为距界面的代数距离)的零水平集隐含地给出。该水平集函数可以很容易地用一个简单的对流方程随着时间推进来更新，从而成功地描述界面的运动状态，而其结构边界演化原理能有效地进行的目标追踪和矢量分析。 所述的VOF模型通过研究网格单元中的流体和网格体积比函数F来确定自由面，追踪流体的变化，而非追踪自由液面上质点的运动。用体积率函数表示流体自由面的位置和流体所占的体积，其方法占内存小，是一种简单而有效的分析流场的方法。VOF方法可以处理自由面重入等强非线性现象，但在处理F的变化时稍显繁琐，有一定人为因素。 所述的水平集的特点决定其能够灵活地描述结构拓扑变化和边界捕获，但在求解过程中存在严重的质量损失，致使计算精度下降。水平集方法与VOF模型的结合(CLSVOF)，是由于VOF模型在边界区域是不连续的，且在处理复杂的几何问题时也较困难，但当与水平集方法结合后，能较好地解决这个问题.同时，也能很好地减小水平集方法的质量损失。同时在处理两相流流场的未知边界时，存在无法使流体的体积模型转化为可计算区域模型的问题，但利用水平集的方法便可以定位边界和体积模型，从而与VOF模型结合，实现对流体模型的模拟与计算。 本专利技术的有益效果主要表现在模具结构化表面软性磨粒流精密加工方法，利用磨粒流与加工表面接触时的壁面效应，形成磨粒对表面的微切削，从而实现表面光整加工；由于磨粒流可与待加工表面形成良好仿形接触，可进行无工具化精密加工，因此其在曲面和异型面加工中具有明显优势。当前关于固-液两相流的研究方法很多，例如VOF模型，但因为值不连续，所以在计算界面的位置和平均曲率(确定表面张力)时精度较低。而水平集方法凭借其在结构拓扑变换方面的优势，能更好地确定相界面的位置，同时与VOF模型的结合能有效地解决质量损失等问题，从而更准确地描述磨粒及磨粒群的运动状态和流体特性，有助于磨粒流本文档来自技高网...

【技术保护点】
一种基于水平集的软性磨粒流流场测试方法，其特征在于：所述测试方法包括以下步骤：１）通过与流体Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程结合，建立软性磨粒流的水平集函数的运动方程；２）对空间作离散化和投影：用一组有限个离散的点来代替原来的连续空间，同时经投影作用使得有散度的转化为无散度的量；３）对时间作离散化；４）水平集函数与流体体积模型函数的结合进行迭代；５）距离函数的重构：重构就是在已知函数分段平均值的情况下，构造满足一定精度下的原函数分布，由于数值方法的内在效应，经过若干演算时间后，原距离函数已不再满足要求，通过求解初始问题的稳定解来实现；６）对软性磨粒流流场的测试：利用水平集方法的有效捕捉界面发生拓扑变化的原理，对流场特性进行较好地分析，找到最佳的加工状态。

【技术特征摘要】

【专利技术属性】
技术研发人员：计时鸣，翁晓星，谭大鹏，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：86[中国|杭州]