高压水射流喷丸的模拟方法

1.本发明属于高压水射流的技术领域，具体涉及一种高压水射流喷丸的模拟方法。


背景技术：

2.在机械工程材料行业中，零部件的使用寿命、可靠性以及产品的竞争性等一直是研究的重点。其中最关键的参数之一便是材料的疲劳强度，疲劳破坏也一直是零部件失效的主要形式。通常来说，改善零部件疲劳强度需要从改善材料表面性能角度进行研究，喷丸强化技术便是其中最有效的方式之一。喷丸工艺可以使材料表面产生塑性变形，提高材料表面的残余压应力，进而提高材料的疲劳性能以及抗应力腐蚀的能力。
3.高压水射流技术就是利用高压水发生装置，将水加压至几百个大气压以上，再通过有细小孔径的喷射装置转换为高速的微细"水射流"。这种"水射流"的速度极高，具有巨大的冲击能量，可以实现切割、清洗、喷丸和除锈等多种任务，其环保高效的特点，被广泛应用于航空航天、汽车、机械、煤矿等行业。高压水射流喷丸就是将高压水射流技术应用于材料表面的喷丸强化方面，以提高材料表面性能的工艺。
4.目前高压水射流喷丸技术在国外有着成熟的应用，国内由于技术等原因还处于发展阶段。近年来随着计算流体动力学(cfd)的发展，关于此类问题已不再局限于通过试验方法来研究，数值模拟技术作为一种成本更低的方案也已成为研究的主要手段之一。
5.目前科研人员常采用abaqus软件，将模型简化为一段高速高压水柱对靶向材料的冲击，利用cel方法进行数值模拟。但是此方案仅能模拟一段小水柱的冲击，不能建立从高压水射流入口到靶体的完整模型。因此，研发一种能完整模拟高压水射流从入口冲击到靶向材料的全过程模型，具有重要的现实和未来意义。


技术实现要素：

6.本发明的目的在于针对现有技术的不足之处，提供一种高压水射流喷丸的模拟方法，能完整模拟高压水射流从入口冲击到靶向材料的全过程模型，既可以实时监测流场中各个位置的物理信息，也可以得到更真实的结论。
7.为解决上述技术问题，本发明采用如下技术方案：
8.一种高压水射流喷丸的模拟方法，包括如下步骤：
9.步骤一：构建高压水射流喷丸的流体域和固体域；
10.步骤二：对流体域进行网格划分，将固体域导入静态结构分析模块的几何结构模块中；
11.步骤三：设置流体域中的气、液两相材料及属性，建立流体输运方程；
12.步骤四：设置流体域的边界条件、选择计算方法和离散方法；
13.步骤五：调整瞬态计算的参数，获得高压水射流喷丸作用下的流场信息；
14.步骤六：读取流场信息，输出靶体所受冲击最大时刻的压力信息；
15.步骤七：将压力信息传输给静态结构分析模块，并将压力加载于靶体受冲击的表
面；
16.步骤八：设置靶体材料，并对其进行网格划分，设置相应的边界条件，最后求解靶体受冲击的情况；
17.步骤九：根据步骤八求解的靶体受冲击情况，得到靶体受到高压水射流喷丸作用下，表面性能的变化。
18.进一步地，流体域包括喷嘴、射流核心区与冲击壁面，固体域包括靶体。
19.进一步地，步骤二中利用icem软件对流体域进行六面体结构网格划分，并对喷嘴和射流核心区域进行网格加密，再转换成非结构网格，并利用ansys-workbench平台导入fluent软件中。
20.进一步地，步骤三中的气相为空气，液相设置为可压缩液体，并采用tait状态方程对其进行建模；其中，tait方程是在等温条件下建立了液体密度与压力之间的非线性关系，其简化形式可以写为：
[0021][0022]
式中：
[0023]
k＝k0+nδp；
[0024]
δp＝p-p0；
[0025]
式中，p0为参考压力；ρ0为参考压力下的液体密度；k0为参考压力下液体的弹性模量；n为密度指数；p为液体绝对压力；ρ为实际压力下的液体密度；k为实际压力下液体的弹性模量；δp表示液体表压力。
[0026]
进一步地，步骤三中，选择多相流模型、湍流模型和壁面函数，建立流体输运方程。
[0027]
进一步地，多相流模型选择为vof模型，湍流模型选择为k-ε模型，壁面函数选择为标准壁面函数swf。
[0028]
进一步地，步骤四中，边界条件设置为压力入口-压力出口，出口与外界相连；求解器设置为压力基求解器，压力-速度耦合算法采用simple分离算法，空间和时间离散格式为一阶迎风格式。
[0029]
进一步地，步骤八具体包括：
[0030]
利用瞬态结构分析中的网格划分功能，对靶体进行网格划分，设置全局网格尺寸，并对受冲击表面进行面加密；
[0031]
将靶体底面固定，并加载外部压力数据至受冲击表面，然后对靶体的应力应变情况进行求解。
[0032]
与现有技术相比，本发明的有益效果为：本发明实现了对高压水射流喷丸工艺流体域的仿真，具体地，本发明利用ansys-workbench平台以及多相流流固耦合技术，将高压水射流喷丸中流体仿真与固体仿真进行了结合，建立了对高压水射流喷丸数值模拟的全过程模型；在流体仿真中，利用cfd数值模拟方法对高压水射流喷丸流体域进行分析，可以实时监测流场中各个位置的物理信息；在固体仿真中，可以实现与流体仿真的耦合，无需借助实验手段，便可得到更真实的结论。
附图说明
[0033]
图1为本发明实施例高压水射流喷丸的模拟方法的流程图；
[0034]
图2为本发明实施例中高压水射流喷丸流体域和固体域的参数化模型示意图；
[0035]
图3为本发明实施例中高压水射流喷丸流体域和固体域的三维模型示意图；
[0036]
图4为本发明实施例中高压水射流喷丸流体域的六面体非结构网格示意图；
[0037]
图5为本发明实施例中流体域截面液相体积分数云图；
[0038]
图6为本发明实施例中高压水射流喷丸固体域的网格示意图；
[0039]
图7为本发明实施例中靶体所受冲击载荷作用下的应力云图。
具体实施方式
[0040]
下面将结合本发明实施例对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0041]
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
[0042]
下面结合具体实施例对本发明作进一步说明，但不作为本发明的限定。
[0043]
如图1所示，本发明实施例提供一种高压水射流喷丸的模拟方法，包括以下步骤：
[0044]
步骤一：利用三维建模软件构建高压水射流喷丸的流体域和固体域；
[0045]
该步骤中，利用workbench平台内的spaceclaim三维绘图软件绘制流体域和固体域的模型，二者之间受高压水射流冲击的表面为流固耦合交界面，压力数据将通过这个交界面从流体传递至固体，其参数化模型如图2所示；其中流体域主要包括喷嘴、射流核心区与冲击壁面三个部分，喷嘴直径d＝0.3mm，靶距s＝1mm。固体域为靶体，在本实施例中，其尺寸形状为直径30mm、高5mm的圆柱体；流体域和固体域的参数化模型如图2所示，两者的三维模型示意图如图3所示。
[0046]
步骤二：利用workbench平台，将流体域导入网格划分软件并对流体域进行网格划分，再将其导入至cfd软件中，此外，将固体域导入静态结构分析模块的几何结构模块中；
[0047]
在该步骤中，利用workbench平台，将三维流体域导入icem网格划分软件中，抑制固体域，仅对流体域进行网格划分；具体地，cfd软件为fluent，利用icem软件对流体域进行六面体结构网格划分，并对喷嘴和射流核心区域进行网格加密，再转换成非结构网格，如图4所示，并利用ansys-workbench平台导入fluent软件中；
[0048]
通过workbench平台将三维模型数据导入静态结构分析中，并抑制流体域。
[0049]
步骤三：设置流体域中的气、液两相材料及属性，并选择多相流模型、湍流模型和壁面函数，从而建立流体输运方程；
[0050]
流体域内涉及气、液两相，气相为空气，液相为水；在本实施例中，设置空气密度为1.225kg/m3，粘度为1.7894
×
10-5
kg/(m
·
s)；水设置为可压缩液体，并采用tait状态方程对其进行建模，tait方程的简化形式可以写为：
[0051]
[0052]
式中：k＝k0+nδp；δp＝p-p0；
[0053]
其中：p0为参考压力，取值101325pa；ρ0为参考压力下的液体密度，取值998.2kg/m3；k0为参考压力下液体的弹性模量，取值2.2
×
109pa；n为密度指数，取值7.15；p为液体绝对压力；ρ为实际压力下的液体密度；k为实际压力下液体的弹性模量；δp表示液体表压力。水的粘度设置为0.001003kg/(m
·
s)；
[0054]
多相流模型选择为vof模型，湍流模型选择为standard k-ε模型，壁面函数选择为标准壁面函数(swf)，建立流体输运方程；
[0055]
其中，液相体积分数的连续性方程和多相流的动量方程为：
[0056][0057][0058]
其中，t是时间；ρ
l
是液相水的密度；α
l
是待求解的液相水体积分数；u
l
是液相水在流场中的速度；ρ是气液混合相密度；u是气液混合相速度；μ是湍流项动力粘度系数；g是重力加速度；f是体积力；。
[0059]
湍流模型的输运方程组为：
[0060][0061][0062]
式中，μ
t
是涡流粘度，其计算式为：
[0063][0064]
其中：t是时间；ρ是气液混合相密度；k是湍流动能；ε是湍流耗散率；u是流体速度；μ是湍流项动力粘度系数；sk、s
ε
是自定义源项；gk是由于层流速度梯度而产生的湍流动能；gb是由于浮力而产生的湍流动能；ym是在可压缩湍流中，过渡的扩散产生的波动；c
1ε
、c
2ε
、c
3ε
、c
μ
是模型常量，其值分别为1.44，1.92，1，0.09；σk、σ
ε
分别为k和ε的prandtl数，其值分别为1.0，1.3。
[0065]
步骤四：设置流体域的边界条件、选择计算方法和离散方法；
[0066]
依据相关文献边界条件设置为压力入口-压力出口，参考压力设置为101325pa，入口表压设置为60mpa，出口与外界相连，表压设置为0mpa。求解器设置为压力基求解器，压力-速度耦合算法采用simple分离算法，空间和时间离散格式为一阶迎风格式。
[0067]
步骤五：调整瞬态计算的参数，获得高压水射流喷丸作用下的流场信息；
[0068]
在本实施例中，瞬态计算的收敛残差设置为0.001，时间步长为1
×
10-7
s，迭代时间步数设为100，每时间步最大迭代步数设置为50，求解可得到高压水射流经过喉管冲击到材料表面的云图，如图5所示。
[0069]
步骤六：利用cfd-post软件读取流场信息，输出靶体所受冲击最大时刻的压力文件；
[0070]
在本实施例中，利用workbench平台内的cfd-post后处理软件读取在fluent软件中计算得到的高压水射流喷丸流场参数，并将流体域中的壁面所受冲击载荷文件导出，文件格式为*.axdt文件。
[0071]
步骤七：利用workbench中的外部数据软件读入压力文件，并传输给静态结构分析模块，再将压力加载于靶体受冲击的表面；
[0072]
在本实施例中，通过workbench平台内的外部数据软件可以读取步骤七中壁面所受冲击载荷的*.axdt格式文件，其中几何模型数据已在步骤二中被导入，再设置固定支撑，即可对模型进行求解。
[0073]
步骤八：设置靶体材料，并对其进行网格划分，设置相应的边界条件，最后求解；
[0074]
本实例的固体材料选用al6061-t6，材料属性：密度为2900kg/m3，杨氏模量为69gpa，泊松比为0.33；
[0075]
利用workbench瞬态结构分析中的mesh网格划分功能，对靶体进行网格划分，设置全局网格尺寸为0.5mm，并对受冲击表面进行面加密，最终得到的网格节点数为2538659，如图6所示；
[0076]
设置的边界条件为：将靶体底面固定、并加载外部压力数据至受冲击表面，根据网格划分结果和设置的边界条件，然后对靶体的应力应变等情况进行求解。
[0077]
步骤九：基于多相流流固耦合仿真结果，得到靶体在高压水射流喷丸作用下，表面性能的变化。在本实施例中，靶体所受冲击载荷作用下的等效应力云图如图7所示。
[0078]
以上仅为本发明较佳的实施例，并非因此限制本发明的实施方式及保护范围，对于本领域技术人员而言，应当能够意识到凡运用本发明说明书内容所作出的等同替换和显而易见的变化所得到的方案，均应当包含在本发明的保护范围内。