首页 / 专利库 / 电脑图像 / 三维网格 / 一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法

一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法

阅读:1039发布:2020-06-15

专利汇可以提供一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法专利检索,专利查询,专利分析的服务。并且本 发明 涉及了一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法,本发明计算求解过程包括:(1)内腔流道模型建立;(2)内腔流道模型网格划分;(3)启动求解器;(4)求解器选择;(5)导入网格进行检查网格;(6)固液两相磨粒流参数设置;(7)流场参数的设置;(8)流场初始化、计算及结果检查;(9)后处理;(10)不同亚格子模型的模拟结果分析:采用大涡模拟的数值研究方法分析四种常见的亚格子模型对磨粒流数值计算结果的影响,对不同亚格子模型的静态压 力 、动压、速度、壁面剪切力和涡旋分布等进行了数值分析;(11)最佳亚格子模型:分析不同亚格子模型 对流 场的预测能力以及扑捉细节的能力,得到最佳亚格子模型。,下面是一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法专利的具体信息内容。

1.一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法,其特征在于;其具体步骤如下:
(1)内腔流道模型建立:选择直管摆线型流道为研究对象,工件流道的最大直径是5mm,流道沟槽处的半圆直径为1mm,距离流道中心最近处的圆弧半径为1.5mm,利用SolidWorks软件对直管摆线型流道进行三维建模,采用不同的亚格子模型对直管摆线型流道进行大涡数值模拟;
(2)内腔流道模型网格划分:结构化网格比非结构化网格收敛速度快,计算过程不容易出现错误,因此数值计算之前需要进行网格无关性检验,对模型以六面体划分网格,流道的网格划分采用最小网格单元数量,直管摆线型流道的网格数量为32.1万;
(3)启动求解器;
(4)求解器选择:3D双精度非稳态压基求解器;
(5)导入网格进行检查网格;
(6)固液两相磨粒流参数设置:磨粒流为固液两相流,采用多相流模型为混合(Mixture)模型,磨粒流的固相为颗粒,密度为3170kg/m3,液相为油,密度为1260kg/m3;
(7)流场参数的设置:
(a)模型选择及条件计算控制参数设定:内腔流道为复杂曲面,为了提高求解精度,采用的湍流模型是为LES模型,离散方法采用有限体积法进行离散,动量方程采用边界中心差分格式,体积分数采用一阶迎格式,瞬态方程采用二阶隐式格式,结合压力-速度耦合的SIMPEC算法对流场进行求解;
(b)入口边界条件:选择压力入口为边界条件;
(c)出口边界条件:由于在固液两相磨粒流加工机床的特性,压力出口为标准大气压
(8)流场初始化、计算及结果检查;
(9)后处理:应用后处理软件处理结果;
(10)大涡模拟中应用不同亚格子模型的模拟结果分析;采用大涡模拟的数值研究方法分析四种常见的亚格子模型对磨粒流数值计算结果的影响,对不同亚格子模型的静态压力、动压、速度、壁面剪切力和涡旋分布等进行了数值分析;
(11)最佳亚格子模型:分析不同亚格子模型对流场的预测能力以及扑捉细节的能力,得到针对此种流道的最佳亚格子模型。
2.根据权利要求1所述的一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法,其特征在于;分别采用四种亚格子模型分别为:Smagorinsky、WALE、WMLES和KET模型对直流道进行数值计算,设置压力进口为2Mpa;
(1)四种亚格子模型的静态压力呈现相同的变化趋势,在进口处压力最大,在出口处压力最小,压力分布呈现一个很明显的阶梯状,说明磨粒流的研抛效果由强到弱,但四种亚格子模型的静态压力数值大小存在差异,在相同的初始条件下采用不同的亚格子模型得到不同的计算结果,说明不同的亚格子模型有不同的数值计算精度;KET模型动态大值要大于Smagorinsky模型、WALE模型和WMLES模型,说明KET模型的捕捉能力要优于其余三个亚格子模型,再结合磨粒流的研抛规律,压力越大表面质量越好,因此可认为KET模型是符合磨粒流研抛的亚格子模型;
(2)四种亚格子模型的动态压力有相同的变化趋势,在流道中间位置,动态压力是一个逐渐变大的过程,在进口处动态压力较小,在出口处动态压力达到最大;而在壁面处动态压力则是一个缓慢减小的过程,说明磨粒流对壁面的研抛效果逐渐减弱,四种亚格子模型的数值基本一致,但仍存在细小偏差;KET亚格子模型的计算结果最大根据磨粒流抛光质量的影响规律,入口压力在某个范围内和抛光质量成正比,压力越大,磨粒流研抛质量越好,因此可以认为KET亚格子模型是较合适的;
(3)四种亚格子模型的速度分布也十分相近,速度越大,磨粒流与壁面之间的作用效果越强,磨粒流对壁面的剪切力越大,磨粒流对工件内表面的抛光质量越好,四种亚格子模型的速度分布和数值大小十分相近,KET亚格子模型的数值相对其他三个亚格子模型更大,WALE亚格子模型和KET亚格子模型比其他两个亚格子模型先达到最大值,Smagorinsky模型和WMLES模型则稍微有点延后,说明WALE和KET亚格子模型能够捕捉到微小的瞬态流场变化,有较高的数值模拟精度;
(4)磨粒流研抛效果取决于磨粒流对壁面的作用效果,也就是壁面剪切力,可通过对壁面剪切力的分析来判断不同亚格子模型的数值模拟精度,四种亚格子模型在进口处的壁面剪切力要大于出口处,随着远离进口,壁面剪切力越来越小,表明磨粒流对进口处的抛光效果好,对出口处的抛光效果较差;KET亚格子模型的数值最大,该分布趋势和压力、速度和动态压力的分布一致;Smagorinsky亚格子模型和WMLES亚格子模型的壁面数值变化跨度大,没有一个逐步递减的过程,表明这两个亚格子模型的捕捉能力较差,模拟精度较低;WALE亚格子模型和KET亚格子模型剪切力数值变化平缓,说明这两个亚格子模型的数值模拟精度较高,能够实时监测流体的瞬态运动变化,捕捉流场的流动细节;
(5)分析流体的运动轨迹,从涡旋的大小和疏密程度可以推测该位置的动能大小和能量传递,从而可以推测磨粒流的对该位置的研抛效果,涡旋的分布基本集中在位置2和4中间位置,尤其是大涡,而位置1和3的分布较小,这是因为流道中部为主流区域,大涡能够得到发展;大涡携带能量,大涡分布越多,携带能量越多,动能越大,对工件壁面的研抛效果越好;大涡都是在主流区域形成,即位置2和4之间,在位置1和3等四个位置出现小涡,这说明主流区域的能量较大,流体动能发展的比较充分,流体对位置2和4的作用效果相比位置1和
3要明显,磨粒流对位置2和4的研抛效果要优于位置1和3;根据流线分布可知WALE亚格子模型和KET亚格子模型的涡旋分布比其他两个亚格子模型分布更密集且涡旋也更大,表明WALE亚格子模型和KET亚格子模型捕捉流场的变化能力更强,对流体的动态变化和流动细节更加精确的模拟;在这四个亚格子模型,WALE亚格子模型和KET亚格子模型的模拟结果出现较多的小涡,其余两个亚格子模型则分布较少,表明WALE亚格子模型和KET亚格子模型对流体的微弱流动把握的更加细致,模拟结果更精准。

说明书全文

一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法

技术领域

[0001] 本发明涉及一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法,属于固液两相磨粒流精密加工技术领域。

背景技术

[0002] 大涡模拟的关键是亚格子模型的选取,不同的亚格子模型适应不同的流场需求,亚格子模型在燃烧、噪声、绕流等领域都进行了相关研究,主要是考察涡旋分布的分离和附着;而在磨粒流研抛领域,一般是考察速度、压等参数,因为最终的结果是看磨粒流研抛效果,这些参数与磨粒流的研抛情况有直接关系,当然涡旋分布也会对磨粒流的研抛质量起着一定作用;大涡模拟方法为瞬态求解,可以捕捉整个流场的瞬态变化,获得磨粒流在研抛过程中的最大值,不同的亚格子模型会得到不同的数值计算结果;通过对流体的压力、动态压力、速度和壁面剪切力等参数以及流线分布图进行分析,观察不同亚格子模型之间的差异和优势,确定合适的亚格子模型用于磨粒流研抛数值模拟;利用最佳亚格子模型进行数值模拟,对磨粒流的研抛效果进行合理预测,这一发明为磨粒流数值模拟提供一定的理论和技术支持。

发明内容

[0003] 本发明的目的在于提供一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法,该方法的取得能为磨粒流数值模拟提供重要的技术指导。
[0004] 为了实现上述目的,一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法的技术方案如下:
[0005] (1)内腔流道模型建立:选择直管摆线型流道为研究对象,工件流道的最大直径是5mm,流道沟槽处的半圆直径为1mm,距离流道中心最近处的圆弧半径为1.5mm,利用SolidWorks软件对直管摆线型流道进行三维建模,采用不同的亚格子模型对直管摆线型流道进行大涡数值模拟;
[0006] (2)内腔流道模型网格划分:结构化网格比非结构化网格收敛速度快,计算过程不容易出现错误,因此数值计算之前需要进行网格无关性检验,对模型以六面体划分网格,流道的网格划分采用最小网格单元数量,直管摆线型流道的网格数量为32.1万;
[0007] (3)启动求解器;
[0008] (4)求解器选择:3D双精度非稳态压力基求解器;
[0009] (5)导入网格进行检查网格;
[0010] (6)固液两相磨粒流参数设置:磨粒流为固液两相流,采用多相流模型为混合(Mixture)模型,磨粒流的固相为颗粒,密度为3170kg/m3,液相为油,密度为1260kg/m3;
[0011] (7)流场参数的设置:
[0012] (a)模型选择及条件计算控制参数设定:内腔流道为复杂曲面,为了提高求解精度,采用的湍流模型是为LES模型,离散方法采用有限体积法进行离散,动量方程采用边界中心差分格式,体积分数采用一阶迎格式,瞬态方程采用二阶隐式格式,结合压力-速度耦合的SIMPEC算法对流场进行求解;
[0013] (b)入口边界条件:选择压力入口为边界条件;
[0014] (c)出口边界条件:由于在固液两相磨粒流加工机床的特性,压力出口为标准大气压
[0015] (8)流场初始化、计算及结果检查;
[0016] (9)后处理:应用后处理软件处理结果;
[0017] (10)大涡模拟中应用不同亚格子模型的模拟结果分析;采用大涡模拟的数值研究方法分析四种常见的亚格子模型对磨粒流数值计算结果的影响,对不同亚格子模型的静态压力、动压、速度、壁面剪切力和涡旋分布等进行了数值分析;
[0018] (11)最佳亚格子模型:分析不同亚格子模型对流场的预测能力以及扑捉细节的能力,得到针对此种流道的最佳亚格子模型。
[0019] 进一步,采用四种亚格子模型分别为:Smagorinsky、WALE、WMLES和KET模型对直流道进行数值计算,设置压力进口为2Mpa;
[0020] (1)四种亚格子模型的静态压力呈现相同的变化趋势,在进口处压力最大,在出口处压力最小,压力分布呈现一个很明显的阶梯状,说明磨粒流的研抛效果由强到弱,但四种亚格子模型的静态压力数值大小存在差异,在相同的初始条件下采用不同的亚格子模型得到不同的计算结果,说明不同的亚格子模型有不同的数值计算精度;KET模型动态大值要大于Smagorinsky模型、WALE模型和WMLES模型,说明KET模型的捕捉能力要优于其余三个亚格子模型,再结合磨粒流的研抛规律,压力越大表面质量越好,因此可认为KET模型是符合磨粒流研抛的亚格子模型;
[0021] (2)四种亚格子模型的动态压力有相同的变化趋势,在流道中间位置,动态压力是一个逐渐变大的过程,在进口处动态压力较小,在出口处动态压力达到最大;而在壁面处动态压力则是一个缓慢减小的过程,说明磨粒流对壁面的研抛效果逐渐减弱,四种亚格子模型的数值基本一致,但仍存在细小偏差;KET亚格子模型的计算结果最大根据磨粒流抛光质量的影响规律,入口压力在某个范围内和抛光质量成正比,压力越大,磨粒流研抛质量越好,因此可以认为KET亚格子模型是较合适的;
[0022] (3)四种亚格子模型的速度分布也十分相近,速度越大,磨粒流与壁面之间的作用效果越强,磨粒流对壁面的剪切力越大,磨粒流对工件内表面的抛光质量越好,四种亚格子模型的速度分布和数值大小十分相近,KET亚格子模型的数值相对其他三个亚格子模型更大,WALE亚格子模型和KET亚格子模型比其他两个亚格子模型先达到最大值,Smagorinsky模型和WMLES模型则稍微有点延后,说明WALE和KET亚格子模型能够捕捉到微小的瞬态流场变化,有较高的数值模拟精度;
[0023] (4)磨粒流研抛效果取决于磨粒流对壁面的作用效果,也就是壁面剪切力,可通过对壁面剪切力的分析来判断不同亚格子模型的数值模拟精度,四种亚格子模型在进口处的壁面剪切力要大于出口处,随着远离进口,壁面剪切力越来越小,表明磨粒流对进口处的抛光效果好,对出口处的抛光效果较差;KET亚格子模型的数值最大,该分布趋势和压力、速度和动态压力的分布一致;Smagorinsky亚格子模型和WMLES亚格子模型的壁面数值变化跨度大,没有一个逐步递减的过程,表明这两个亚格子模型的捕捉能力较差,模拟精度较低;WALE亚格子模型和KET亚格子模型剪切力数值变化平缓,说明这两个亚格子模型的数值模拟精度较高,能够实时监测流体的瞬态运动变化,捕捉流场的流动细节;
[0024] (5)分析流体的运动轨迹,从涡旋的大小和疏密程度可以推测该位置的动能大小和能量传递,从而可以推测磨粒流的对该位置的研抛效果,涡旋的分布基本集中在位置2和4中间位置,尤其是大涡,而位置1和3的分布较小,这是因为流道中部为主流区域,大涡能够得到发展;大涡携带能量,大涡分布越多,携带能量越多,动能越大,对工件壁面的研抛效果越好;大涡都是在主流区域形成,即位置2和4之间,在位置1和3等四个位置出现小涡,这说明主流区域的能量较大,流体动能发展的比较充分,流体对位置2和4的作用效果相比位置1和3要明显,磨粒流对位置2和4的研抛效果要优于位置1和3;根据流线分布可知WALE亚格子模型和KET亚格子模型的涡旋分布比其他两个亚格子模型分布更密集且涡旋也更大,表明WALE亚格子模型和KET亚格子模型捕捉流场的变化能力更强,对流体的动态变化和流动细节更加精确的模拟;在这四个亚格子模型,WALE亚格子模型和KET亚格子模型的模拟结果出现较多的小涡,其余两个亚格子模型则分布较少,表明WALE亚格子模型和KET亚格子模型对流体的微弱流动把握的更加细致,模拟结果更精准。
[0025] 本发明的有益效果在于:本发明提供了一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法,磨粒流抛光质量受到多种因素的影响,例如挤出压力、介质粘度、磨粒尺寸等,为了更精确的描述磨粒流加工过程,对数值模拟过程中的亚格子模型进行优选,不同的模型对流场模拟的精确程度不同,利用最佳亚格子模型进行数值模拟,对磨粒流的研抛效果进行合理预测,能够实时监测流体的瞬态运动变化,模拟更加精准。附图说明
[0026] 图1一种固液两相磨粒流加工内腔流道亚格子模型数值选择方法流程图
[0027] 图2直管摆线型流道模型;
[0028] 图3流道截面位置划分;

具体实施方式

[0029] 以下结合附图来对本发明进行详细的描绘。然而应当理解,附图的提供仅为了更好的理解本发明,它们不应该理解成对本发明的限制。
[0030] (1)内腔流道模型建立:选择直管摆线型流道为研究对象,工件流道的最大直径是5mm,流道沟槽处的半圆直径为1mm,距离流道中心最近处的圆弧半径为1.5mm,利用SolidWorks软件对直管摆线型流道进行三维建模,采用不同的亚格子模型对直管摆线型流道进行大涡数值模拟;
[0031] (2)内腔流道模型网格划分:结构化网格比非结构化网格收敛速度快,计算过程不容易出现错误,因此数值计算之前需要进行网格无关性检验,对模型以六面体划分网格,流道的网格划分采用最小网格单元数量,直流道的网格数量为32.1万;
[0032] (3)启动求解器;
[0033] (4)求解器选择:3D双精度非稳态压力基求解器;
[0034] (5)导入网格进行检查网格;
[0035] (6)固液两相磨粒流参数设置:磨粒流为固液两相流,采用多相流模型为混合(Mixture)模型,磨粒流的固相为碳化硅颗粒,密度为3170kg/m3,液相为油,密度为1260kg/m3;
[0036] (7)流场参数的设置:
[0037] (a)模型选择及条件计算控制参数设定:内腔流道为复杂曲面,为了提高求解精度,采用的湍流模型是为LES模型,离散方法采用有限体积法进行离散,动量方程采用边界中心差分格式,体积分数采用一阶迎风格式,瞬态方程采用二阶隐式格式,结合压力-速度耦合的SIMPEC算法对流场进行求解;
[0038] (b)入口边界条件:选择压力入口为边界条件;
[0039] (c)出口边界条件:由于在固液两相磨粒流加工机床的特性,压力出口为标准大气压;
[0040] (8)流场初始化、计算及结果检查;
[0041] (9)后处理:应用后处理软件处理结果;
[0042] (10)大涡模拟中应用不同亚格子模型的模拟结果分析;采用大涡模拟的数值研究方法分析四种常见的亚格子模型对磨粒流数值计算结果的影响,对不同亚格子模型的静态压力、动压、速度、壁面剪切力和涡旋分布等进行了数值分析;
[0043] (11)最佳亚格子模型:分析不同亚格子模型对流场的预测能力以及扑捉细节的能力,得到针对此种流道的最佳亚格子模型;
[0044] 采用四种亚格子模型分别为:Smagorinsky、WALE、WMLES和KET模型对直流道进行数值计算,设置压力进口为2Mpa;
[0045] (1)四种亚格子模型的静态压力呈现相同的变化趋势,在进口处压力最大,在出口处压力最小,压力分布呈现一个很明显的阶梯状,说明磨粒流的研抛效果由强到弱,但四种亚格子模型的静态压力数值大小存在差异,在相同的初始条件下采用不同的亚格子模型得到不同的计算结果,说明不同的亚格子模型有不同的数值计算精度;KET模型动态大值要大于Smagorinsky模型、WALE模型和WMLES模型,说明KET模型的捕捉能力要优于其余三个亚格子模型,再结合磨粒流的研抛规律,压力越大表面质量越好,因此可认为KET模型是符合磨粒流研抛的亚格子模型;
[0046] (2)四种亚格子模型的动态压力有相同的变化趋势,在流道中间位置,动态压力是一个逐渐变大的过程,在进口处动态压力较小,在出口处动态压力达到最大;而在壁面处动态压力则是一个缓慢减小的过程,说明磨粒流对壁面的研抛效果逐渐减弱,四种亚格子模型的数值基本一致,但仍存在细小偏差;KET亚格子模型的计算结果最大根据磨粒流抛光质量的影响规律,入口压力在某个范围内和抛光质量成正比,压力越大,磨粒流研抛质量越好,因此可以认为KET亚格子模型是较合适的;
[0047] (3)四种亚格子模型的速度分布也十分相近,速度越大,磨粒流与壁面之间的作用效果越强,磨粒流对壁面的剪切力越大,磨粒流对工件内表面的抛光质量越好,四种亚格子模型的速度分布和数值大小十分相近,KET亚格子模型的数值相对其他三个亚格子模型更大,WALE亚格子模型和KET亚格子模型比其他两个亚格子模型先达到最大值,Smagorinsky模型和WMLES模型则稍微有点延后,说明WALE和KET亚格子模型能够捕捉到微小的瞬态流场变化,有较高的数值模拟精度;
[0048] (4)磨粒流研抛效果取决于磨粒流对壁面的作用效果,也就是壁面剪切力,可通过对壁面剪切力的分析来判断不同亚格子模型的数值模拟精度,四种亚格子模型在进口处的壁面剪切力要大于出口处,随着远离进口,壁面剪切力越来越小,表明磨粒流对进口处的抛光效果好,对出口处的抛光效果较差;KET亚格子模型的数值最大,该分布趋势和压力、速度和动态压力的分布一致;Smagorinsky亚格子模型和WMLES亚格子模型的壁面数值变化跨度大,没有一个逐步递减的过程,表明这两个亚格子模型的捕捉能力较差,模拟精度较低;WALE亚格子模型和KET亚格子模型剪切力数值变化平缓,说明这两个亚格子模型的数值模拟精度较高,能够实时监测流体的瞬态运动变化,捕捉流场的流动细节;
[0049] (5)分析流体的运动轨迹,从涡旋的大小和疏密程度可以推测该位置的动能大小和能量传递,从而可以推测磨粒流的对该位置的研抛效果,涡旋的分布基本集中在位置2和4中间位置,尤其是大涡,而位置1和3的分布较小,这是因为流道中部为主流区域,大涡能够得到发展;大涡携带能量,大涡分布越多,携带能量越多,动能越大,对工件壁面的研抛效果越好;大涡都是在主流区域形成,即位置2和4之间,在位置1和3等四个位置出现小涡,这说明主流区域的能量较大,流体动能发展的比较充分,流体对位置2和4的作用效果相比位置1和3要明显,磨粒流对位置2和4的研抛效果要优于位置1和3;根据流线分布可知WALE亚格子模型和KET亚格子模型的涡旋分布比其他两个亚格子模型分布更密集且涡旋也更大,表明WALE亚格子模型和KET亚格子模型捕捉流场的变化能力更强,对流体的动态变化和流动细节更加精确的模拟;在这四个亚格子模型,WALE亚格子模型和KET亚格子模型的模拟结果出现较多的小涡,其余两个亚格子模型则分布较少,表明WALE亚格子模型和KET亚格子模型对流体的微弱流动把握的更加细致,模拟结果更精准。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈