[发明专利]基于近壁区理论的软性磨粒流光整加工控制方法

摘要

基于近壁区理论的软性磨粒流光整加工控制方法，包括以下步骤：建立软性磨粒流运动的湍流数学模型：采用可实现k-ε双方程湍流模型对流体进行分析；采用欧拉多相流模型方程作为软性磨粒流的控制方程；软性磨粒流的数值模拟求解：网格划分；初始条件和边界条件；离散相模型初始化；方程离散化；给定求解控制参数；迭代运算求解。本发明具有能够提供光整加工时的理论支持，明确加工机理，提高抛光效率的优点。

主权项

1.基于近壁区理论的软性磨粒流光整加工控制方法，包括以下步骤：1)、建立软性磨粒流运动的湍流数学模型：(1.1)、采用可实现k-ε双方程湍流模型(Realizable k-ε Model)对流体进行分析，k和ε的输运方程如下：∂(ρk)∂t+∂(ρkui)∂xi=∂∂xj[(μ+μtσk)∂k∂xj]+Gk-ρϵ---(1)]]>∂(ρϵ)∂t+∂(ρϵui)∂xi=∂∂xj[(μ+μtσϵ)∂ϵ∂xj]+ρC1Eϵ-ρC2ϵ2k+vϵ---(2)]]>其中：ρ为流体密度；x_i，x_j为各坐标分量；σ_k，σ_ε为湍动能k和耗散率ε的湍流普朗特数(Prandtl Number)σ_k＝1.0，σ_s＝1.2；C₂＝1.9；C1=max(0.43,ηη+5);---(3)]]>η=(2EijgEij)12kϵ,]]>Eij=12(∂ui∂xj+∂uj∂xi);]]>C_k表示由于平均速度梯度引起的湍动能产生；μ为分子粘性系数，μ_t为湍流粘性系数，μt=ρCμk2ϵ;---(4)]]>系数Cμ=1A0+ASU*kϵ;---(5)]]>(1.2)、采用欧拉多相流模型方程作为软性磨粒流的控制方程，欧拉模型的第i相质量守恒方程为：第i相的动量守恒方程为：式(4)，(5)中，α_i为体积分数，ρ_i是第i相的流体密度，v_i是第i相的速度，表示从第i相到第j相的质量传递；是i，j两相之间的相互作用力，p是所有相共同受到的压力，是i，j两相之间的速度；是外部体积力，是i相作用于j相的升力，表示为：是i相粒子对j相加速粒子施加的虚拟质量力，表示为：FVvm,i=0.5ρiαj(dividt-djvjdt)---(9)]]>为压力应变张量，表示为：τi==αiμi(▿viV+▿viVT)+αi(λi-23μi)▿·viVI=---(10)]]>2)、软性磨粒流的数值模拟求解：(2.1)网格划分：根据软性磨粒流流道中区域的相异性，在不同的计算区域采用结构化和非结构化网格交替划分的方法来对数值模拟中的流道进行网格划分；(2.2)初始条件和边界条件：设置软性磨粒流在流道中运动的初始条件和边界条件；(2.3)离散相模型(Discrete Phase Model DPM)初始化：采用DPM对软性磨粒流在流道内的轨迹进行仿真模拟，对DPM进行初始化；(2.4)方程离散化：利用数值模拟方法，将作为计算域的连续流道壁面离散成为相互独立的网格，将其视为未知变量，并建立以未知变量为参数的代数方程，通过对求解代数方程组求解来得到节点值，而计算区域内的其他位置上的值则根据响应节点值来确定；(2.5)给定求解控制参数：针对代数离散方程组，对流体的物理性质，软性磨粒流的湍流模型参数，迭代精度，计算步长和频率等进行设定；(2.6)迭代运算求解：对数值运算中所需要的解进行多次迭代计算，在运算结果的残差收敛后，迭代过程结束。