1. 引言
近年来我国淡水鱼加工市场迅速增长,淡水鱼的前处理工序繁多,其中去鳞加工是鱼原料加工的关键工序,同时也是整个工序中的技术难点 [1]。随着淡水鱼加工量的大幅增长,具体高质量和大批量加工能力的去鳞方法和设备成为行业的亟需势。目前加工企业广泛采用的鱼去鳞方式主要有滚筒式去鳞 [2]、辊刷式去鳞 [3] 和水射流式去鳞 [4]。相对于其他两种方式,水射流式去鳞具有去鳞率高、加工速度快的优点。但仍存在一定不足,水射流式去鳞的参数较多,而鱼类原料大小种类不一,水射流去鳞的优化参数组合困难。当水射流去鳞工艺参数设置不合理时,将出现去鳞效果不佳、鱼体表皮损伤、鱼尾/鱼鳍损伤等问题,影响鱼去鳞加工的合格率和外观品质等。水射流去鳞效果与鱼体受到射流打击力特性密切相关,因此有必要分析水射流参数与鱼体受到射流打击力的关系。
Fluent是分析水射流仿真的常用软件,经常用于工程、农业、食品等领域 [5] [6] [7] [8] 物料表面清洗去渍过程的仿真分析。李仁杰等 [9] 对比鱼去鳞过程的试验和仿真结果,验证了采用Fluent模拟鱼去鳞过程的仿真精度,为本文研究提供了方法参考。
本文基于Fluent建立水射流去鳞的流固耦合分析模型,分析水射流参数变化对鱼体表面应力分布的影响趋势,建立鱼去鳞效果的力学特性评价指标,探讨去鳞工艺参数优化的准则,为鱼类原料水射流式去鳞工艺优化提供参考。
2. 模型与方法
以鱼体受水射流冲击后的应力分布作为主要参考来评定鱼去鳞过程的受力情况。因此本文采用流固耦合计算方法,该方法一般用在流体与固体间存在强烈的相互作用时,对流体流场与固体应力应变的考察,本文采用单向流固耦合进行仿真分析。单向流固耦合分析过程主要包括:流场仿真和结构仿真,每个过程都包括建立几何模型,网络划分,参数设置,求解和后处理。
2.1. 模型建立
本文在研究鱼体定向上选取了鲤鱼作为研究对象,鲤鱼长轴方向上呈长纺锤形,横截面呈楔形,鱼背厚度大于腹背厚度,鱼体重心靠近鱼头部分,本文鲤鱼模型参照李仁杰 [10] 等方法经过对众多鲤鱼样本的体型参数进行测量,经过对大量鲫鱼样本的体型参数进行测量,得到鲫鱼形态参数如表1所示;根据表1利用Solid Works建立鲫鱼近似三维模型如图1所示,通过企业调研和多种喷嘴的比较发现扇形喷嘴喷射比较均匀,适用于去鳞工作,根据鲤鱼的体型特征,25˚扇形喷嘴与其它相比,可以较好的对鲤鱼鱼体进行覆盖。所以选择扇型喷嘴作为喷嘴模型,其孔径为2.4 mm,喷射角度为25˚。以上述鲤鱼和喷嘴模型为基础建立Fluent流场仿真模型如图2所示。
![](Images/Table_Tmp.jpg)
Table 1. Morphological parameter measurement of crucian
表1. 鲫鱼形态参数测量
![](//html.hanspub.org/file/1-2570296x11_hanspub.png)
Figure 2. Fluent flow field simulation model
图2. Fluent流场仿真模型
2.2. 网络划分
采用Ansys Fluent自带的网格划分软件Icem对射流流场进行六边形网格进行划分,除各个入口、出口交接处做加密化处理,其它流场部分采用自动网络划分 [11],流场模型的网格划分如图3所示。
2.3. 参数设置
选用标准k-ε湍流模型对整个气液两相流场进行模拟仿真,k值为0.3,ε值为38.08选择稳流,设置空气为基本相,水为第二相;选择非稳态离散求解器,入口边界定义为压力入口,设定水体积分数为1;外流场边界条件也采用压力边界,出口压力为一个标准大气压,设定回流体积分数为0;喷嘴壁面定义为无滑移壁面,采用PISO算法对压力和速度进行耦合求解计算 [6]。
2.4. 水射流去鳞原理
水射流鱼去鳞的原理如图4所示。上方喷嘴喷射经由水泵加压后的压力水以一定角度打击鱼体表面,随着鱼被水平输送,水射流的剪切作用使鱼鳞从鱼体表面上剥离 [12]。当输送速度一定时,鱼去鳞效果与水射流参数(射流压力、靶距和入射角)有关。
![](//html.hanspub.org/file/1-2570296x13_hanspub.png)
Figure 4. Principle diagram of water jet descaling
图4. 水射流鱼去鳞原理图
3. 仿真与结果分析
3.1. 仿真方案
以射流压力、射流靶距和射流入射角3个参数分别进行单因素仿真试验。试验方案如表2所示。
![](Images/Table_Tmp.jpg)
Table 2. Single factor test protocol
表2. 单因素试验方案
3.2. 仿真结果
当射流入射角和靶距一定时(α = 20˚、L = 100 mm),射流压力8 MPa、9 MPa、10 MPa、11 MPa应力云场如图5所示,从图5可以看出,随着射流压力的增加,鱼体所受的最大等效应力上升,这是因为压力增加,鱼体受到的打击力变大,从而最大等效应力变大;当射流入射角和射流压力一定时(α = 20˚、P = 9 MPa),射流靶距100 mm、130 mm、160 mm、190 mm应力云场如图6所示,从图6可以看出,随着射流靶距的增加,鱼体所受的最大等效应力下降,这是因为由于空化作用,随射流靶距减少,射流打击力会增大,最大等效应力大;当射流压力和靶距一定时(P = 9 MPa、L = 100 mm),射流入射角20˚、25˚、30˚、35˚应力云场如图7所示,从图7可以看出,随着射流入射角的增加,最大等效应力先变大后变小,存在一个最大值,因为当射流角度变小时,虽然射流打击力的垂直分量变小,但射流打击力的水平分量变大。
![](//html.hanspub.org/file/1-2570296x14_hanspub.png)
Figure 5. Stress cloud diagram under different jet pressure
图5. 不同射流压力下应力云图
![](//html.hanspub.org/file/1-2570296x15_hanspub.png)
Figure 6. Stress cloud diagram under different target distances
图6. 不同靶距下应力云图
![](//html.hanspub.org/file/1-2570296x16_hanspub.png)
Figure 7. Stress cloud diagram at different jet incident angles
图7. 不同射流入射角下应力云图
3.3. 力学性能指标建立
水射流去鳞过程可视为射流系统相对静止,鱼体沿着其轴向水平线性往复移动,在此过程中水射流部分的打击在鱼体表面进行去鳞作业。水射流去鳞的力学性能可用鱼体受力均匀性、最大打击力、有效打击比3个指标来衡量 [13]。
受力均匀性采用鱼体表面的应力最大的径向切线上的Von-Mises等效应力分布曲线表示鱼去鳞过程受力的均匀性,其中该曲线的变化趋势越平缓,鱼体表面在去鳞过程中受力的均匀性越好。
最大打击力采用鱼体表面受到的最大Von-Mises等效应力表示最大打击力。最大Von-Mises等效应力变化应折中,使得去鳞率程度高且对鱼体的表面损伤小。
有效打击比在射流压力一定的情况下,以鱼体受到的打击力与射流参数方案能提供的总打击力比值来衡量该方案的射流利用水平,有效打击比越大说明射流参数方案的能耗有效利用率越高。按下式计算有效打击比。
(1)
式中,
——有效打击比,%;
——鱼体受到的打击力,N;
——射流参数方案提供的总打击力,N。
3.4. 受力均匀性的影响因素分析
固定射流入射角和射流压力(α = 20˚、P = 9 MPa)、射流入射角和射流靶距(α = 20˚、L = 100 mm)、射流压力和射流靶距(P = 9 MPa、L = 100 mm),分别建立鱼体表面径向等效应力分布曲线随射流靶距、射流压力及射流入射角变化曲线如图8、图9、图10所示。
由图8、图9、图10可知,射流压力、射流入射角和射流靶距的等效应力分布曲线均呈中间向两端平稳降低的规律,射流压力、射流入射角和射流靶距的变化基本不影响鱼体表面径向等效应力分布趋势。选取中间2/3段(−30 mm, 30 mm)曲线的应力数据,射流压力一定时(P = 9 MPa),射流靶距100 mm、130 mm、160 mm、190 mm时的等效应力最大变化率分别为12.9%,9.4%,7.6%,4.5%,在给定的水射流参数取值范围内,随着射流靶距减少,鱼体受力均匀性下降;射流靶距一定时(射流靶距L = 100 mm),射流压力8 MPa、9 MPa、10 MPa、11 MPa时的等效应力最大变化率分别为15.4%,12.7%,15.2%,15.3%,在给定的水射流参数取值范围内,随着压力变化,鱼体受力均匀性基本不变;射流压力、射流靶距一定时(P = 9 MPa、L = 100 mm),射流入射角20˚、25˚、30˚、35˚时的等效应力最大变化率分别为12.7%,25.0%,25.3%,26.5%,在给定的水射流参数取值范围内,随着射流入射角的增加,鱼体受力均匀性下降。其中,射流方案P = 9 MPa、L = 100 mm、β = 20˚的受力均匀性最好。因此设定射流参数时,在满足提供足够的水射流喷射打击力的同时,应尽量采用较小的射流靶距和射流压力,以提高鱼体受力的均匀性。
![](//html.hanspub.org/file/1-2570296x21_hanspub.png)
Figure 8. Radial equivalent stress distribution curve of fish body surface under different target distances
图8. 不同射流靶距下鱼体表面径向等效应力分布曲线
![](//html.hanspub.org/file/1-2570296x22_hanspub.png)
Figure 9. Radial equivalent stress distribution curve on fish body surface under different inlet pressures
图9. 不同射流压力下时鱼体表面径向等效应力分布曲线
![](//html.hanspub.org/file/1-2570296x23_hanspub.png)
Figure 10. Radial equivalent stress distribution curve on the fish body surface under different jet incident angles
图10. 不同射流入射角下时鱼体表面径向等效应力分布曲线
3.5. 最大打击力的影响因素分析
固定射流入射角(α = 20˚)的条件下,分别建立最大Von-Mises等效应力随射流压力和射流靶距的变化曲线,如图11、图12所示。
从图11、图12可知,在射流入射角一定时,最大Von-Mises等效应力随射流靶距增大下降,在相同射流靶距下,射流压力越大,最大Von-Mises等效应力越大;最大Von-Mises等效应力随射流压力增加上升,在相同射流压力下,射流靶距越小,最大Von-Mises等效应力越大。因此采用较小射流靶距和较大射流压力的参数组合可获得采用较大射流靶距和较小射流压力组合同样的射流打击力效果,增大射流压力将提高射流打击力,为了保证足够的打击力,应选择尽量小的射流靶距。
![](//html.hanspub.org/file/1-2570296x24_hanspub.png)
Figure 11. Variation trend of maximum Von-Mises equivalent stress with target distance
图11. 最大Von-Mises等效应力随射流靶距的变化趋势
![](//html.hanspub.org/file/1-2570296x25_hanspub.png)
Figure 12. The variation trend of the maximum Von-Mises equivalent stress with the inlet pressure
图12. 最大Von-Mises等效应力随射流压力的变化趋势
3.6. 有效打击比的影响因素分析
固定射流压力(P = 9 MPa)的条件下,分别建立有效打击比随射流靶距和射流入射角的变化曲线,如图13、图14所示。
由图13、图14可知,在射流压力一定时,有效打击比随着射流靶距和射流入射角的增加,有效打击比都是下降趋势;射流距离固定时,射流入射角越小,有效打击比越大;射流入射角固定时,射流靶距越小,有效打击比越大;其中仿真方案P = 9 MPa、L = 100 mm、β = 20˚的有效打击比最大,即为97.3%,说明该方案的能耗有效利用率最大。
![](//html.hanspub.org/file/1-2570296x26_hanspub.png)
Figure 13. The effective blow ratio vs. jet distance
图13. 有效打击比随射流距离的变化曲线
![](//html.hanspub.org/file/1-2570296x27_hanspub.png)
Figure 14. The effective blow ratio vs. jet incident angle
图14. 有效打击比随射流入射角的变化曲线
3.7. 射流参数优化
水射流最佳工艺要求是鱼体要有比较好的受力均匀性、比较高的有效打击比、适中的最大打击力;在选取水射流去鳞最优参数组合时,应将受力均匀性、有效打击比和最大打击力3项指标进行综合考虑,综合上述3个指标,通过对各参数组合下的力学特性分析,在试验给定的射流参数范围内,当射流压力9 MPa,射流入射角20˚,射流靶距100 mm,鱼体受力均匀性、有效打击比、最大打击力最优,因此在试验给定射流参数范围内,选择该组参数作为水射流去鳞的最佳射流组合较为合适。
4. 结论
1) 根据水射流去鳞原理,基于Fluent建立水射流去鳞的流固耦合分析模型,实现了连续的鱼腹背定向清洗。
2) 经单因素研究,采用较小的射流靶距和射流压力,可以提高鱼体受力的均匀性;采用较小射流靶距和较大射流压力的参数组合,可获得采用较大射流靶距和较小射流压力组合同样的射流打击力效果;有效打击比随着射流靶距和射流入射角的增加均呈下降趋势。
3) 当射流压力9 MPa,射流入射角20˚,射流靶距100 mm时,鱼体受力均匀性、有效打击比和最大打击力最优,为鱼类原料水射流式去鳞设备提供了参考。
4) 本研究的局限性在于:单一对鲤鱼模型为仿真实验对象,后续可增加其它鱼类对比研究,流场模拟采用的是扇形喷嘴单点喷射方式,后续可对空间多点式水射流去鳞技术进行研究。
基金项目
国家重点研发计划项目(2018YFD0400800);辽宁省自然科学基金项目(20180551017和20180550792)。
NOTES
*通讯作者。