病死动物化制搅拌过程流固耦合数值仿真与优化

第43卷第3期Vol.43No.32022年3月
Mar.2022
中国农机化学报
Journal of Chinese Agricultural Mechanization
DOI:10.13733/j.jcam.issn.2095⁃5553.2022.03.012
病死动物化制搅拌过程流固耦合数值仿真与优化*俞国红1，薛向磊1，郑航1，牛瑞征2，蒋建东2
(1.浙江省农业科学研究院农业装备研究所，杭州市，310021；
2.浙江工业大学特种装备与先进加工技术教育部重点试验室，杭州市，310014)
摘要：针对病死动物干化化制过程中物料组分复杂多变，流动不均衡等问题，提出一种新型双螺带式搅拌器。

采用ANSYS 软件进行了化制罐内物料与搅拌器流固耦合建模与仿真分析，研究不同搅拌速度下化制容器内的物料流动速度与分布规律，并结合干化容器内物料流线分布进行了搅拌器结构优化。

试验和数值优化结果表明，双螺带式搅拌器在10r/min 的搅拌转速下，化制容器内超过80%的物料流动速度均达到0.419m/s 以上。

相比单螺带搅拌器提升物料流动性能20%，一种新型双螺带式搅拌器提升了病死动物干化化制处理的效率和质量。

关键词：化制过程；搅拌器；流固耦合；流动特性；数值仿真中图分类号：S811.5
文献标识码：A
文章编号：2095⁃5553(2022)03⁃0092⁃07
俞国红,薛向磊,郑航,牛瑞征,蒋建东.病死动物化制搅拌过程流固耦合数值仿真与优化[J].中国农机化学报,2022,43(3):92-98
Yu Guohong,Xue Xianglei,Zheng Hang,Niu Ruizheng,Jiang Jiandong.Numerical simulation and optimization of fluid ⁃solid coupling in the mixing process of sick and dead animals [J].Journal of Chinese Agricultural Mechanization,2022,43(3):92-98
0引言
内化制法是将病死动物尸体及其废弃物破碎后投
入到高温高压容器中进行灭菌的工艺方法[1]，化制后的物料进行烘干和油脂分离，得到油脂和骨肉粉等副产品，因处理周期短，灭菌效果好，副产物质量高等优点化制法已成为国内外病死动物无害化处理的主流方法。

化制法经过国内外学者不断探索，已逐渐发展成熟，并得出相关结论：Gwyther 等[2]认为在处理因高致病性疾病死亡的畜禽时温度要求达到133℃，压力要求达到0.3MPa ，处理时间不少20min 的条件；Auvermann 等[3]认为影响化制过程的主要因素包括时间和温度、空气压力以及动物尸体被分割的大小，最好控制在30~50mm 之间；王德胜等[4]对现有干化化制无害化工艺进行利弊分析，通过层次分析法对不同工艺进行评价，得出综合效益最好的干化化制工艺，并对猪肉化制过程物料的特性进行研究，从而优化化制参数。

目前市场上多采用干化化制机完成化制处理过程，化制机中间设计有用于加热搅拌的搅拌单元，搅拌单元包括搅拌轴体以及轴体上安装的具有分解和输送功能的搅拌叶片，搅拌装置承担了病死动物化制处理工艺近70%的工作量，且化制罐内的物料流动直接影响着物料的加热速率和灭菌效果，因此搅拌器进行物料的搅拌直接影响到化制处理的效率和质量，而国内外学者针对干化化制设备研究鲜有报道：王长利等[5]开发了一套采用干化化制法对病死动物进行无害化处理的装置及其处理方法，该装备置处于试验阶段未见推广。

陶兰兰等[6]针对化工设备搅拌器类型与结构开展了研究，结果表明采用螺带式搅拌器更适合于黏度较高的流体化制过程。

陆旭等[7]针对搅拌轴的可靠性优化进行了研究，对空心轴的选用进行了论证，得出了空心轴适用的轴径范围。

张鑫[8]针对离心泵叶片与泵腔的流场与结构开展了研究，苏腾龙等[9]开展柔性浆流场特性的流固耦合数值模拟和试验，二者指出了搅拌器工作过程物料流动
收稿日期：2021年7月15日
修回日期：2021年8月25日
*基金项目：
省长基金项目(2020R30A60C01)——现代化规模猪场建设标准与智能化集成研究与示范第一作者：俞国红，男，1969年生，浙江杭州人，高级工程师；研究方向为农业机械设计与应用。

E⁃mail:Yuguohong69@163．com 通讯作者：蒋建东，男，1974年生，四川广安人，教授；研究方向为机械动力学、机电系统设计与控制。

E⁃mail ：****************.cn
流固耦合特性直接影响搅拌效果。

诸多学者针对搅拌过程中物料流固耦合特性的分析进行了探索[10-20],验证了基于有限元分析原理运用ANSYS 软件中的Fluent 模块研究搅拌器内部的流体特性的可行性。

通过试验测定干化化制过程不同时间阶段的物料特性，进行干化化制设备的有限元数值建模，通过流固耦合分析开展化制过程中物料流动特性数值仿真，开展螺带式搅拌器优化设计与分析，旨在提升病死动物干化化制处理的效率和质量。

1病死动物干化化制设备机理分析
根据农业农村部提出的病死动物无害化处理技术
规范规定：干化处理时物料中心温度不低于140℃，容器内绝对压力不低于0.5MPa ，处理时间需要达到4h 以上；设计新型病死动物高温高压化制处理系统组成如图1所示。

化制设备搅拌器工作环境有以下限制：化制罐为卧式反应容器，搅拌轴轴向跨度较大；容器直径较大，需要的叶片搅拌范围较大；物料黏度大，搅拌速度慢；螺带搅拌器搅拌时可以使流体产生径向流、轴向流和
环向流，搅拌效果好，且适用于黏度较高的流体。

常用的单螺带搅拌器，如图2(a)所示，当螺带的外径D j 为1400mm 时，螺带宽度b 为115mm ，螺距S 为1400mm ，螺带厚度δ为12mm 。

通过运用Fluent 流体仿真软件对单螺带搅拌器搅拌作用下的物料进行流场仿真，分析单螺带搅拌器的搅拌效果，如图2(b)、图2(c)所示。

5
4
3
2
1
76图1
新型无害化设备处理系统
Fig.1
New innocuous equipment processing system
1.输送泵
2.破碎组件
3.化制烘干组件
4.加热组件
5.螺旋输送机
6.粉碎机
7.油脂分离机
(a)单螺带搅拌器模型图
(b)纵截面速度分布云图(c)横截面速度分布云图
0.5001.000 m 0.2500.750
Veloc ity/(m·s -1)
2.0951.8861.6761.4671.2571.0480.8380.6290.4190.2100.000
x
0.5001.000 m
0.2500.750
Y Y
Z
Z X
X
图2
单螺带搅拌效果图
Fig.2
Stirring effect diagram of single ribbon
(a)双螺带搅拌器模型图
(b)纵截面速度分布云图(c)横截面速度分布云图
0.5001.000 m 0.2500.750
00.5001.000 m 0.2500.750
Y Y Z
Z X
X Velocity /(m·s -1)
2.0951.8851.6761.4661.2571.0470.8380.6280.4190.2090.000
Velocity /(m·s -1)
2.0951.8851.6761.4661.2571.0470.8380.628
0.419
0.209
0.000yz
yz
图3
双螺带搅拌效果图
Fig.3
Stirring effect diagram of double ribbon
搅拌轴周围大范围区域的物料流动速度较小，搅拌效果并不理想，且搅拌器左右两端的物料流动性相差较大，物料向搅拌器的右端聚集，会增大右侧搅拌器和罐体的负载，因此需要对单螺带搅拌器进行改进设计。

根据试验，造成搅拌轴周围搅拌效果较差的原因
是螺带与搅拌轴距离较大，为改善搅拌轴周围物料的搅拌效果，本文在搅拌轴和大直径螺带之间加装旋向相反的小直径螺带(图3）：外径700mm ，宽度55mm ，螺距700mm ，螺带厚度8mm 。

改进的搅拌器使物料在小螺旋带的转动下回流，有效改善物料的流动性，减小搅拌器和化制罐的负载。

对改进后的双螺带搅拌器搅拌作用下的物料进行流场仿真，对比图2、图3可知，单螺带搅拌器搅拌作用
下约有60%左右的物料流动速度达到0.3m/s以上，双螺带搅拌器搅拌作用下有80%以上的物料流动速度均达到了0.419m/s以上，双螺带搅拌器搅拌作用下，物料向搅拌器一端聚集的现象得到明显缓解，因此双螺带搅拌器的搅拌效果明显优于单螺带搅拌器，有效改善了搅拌器的搅拌效果。

2化制容器内搅拌过程物料流动模型与特性分析
病死动物化制处理的物料主要是经过破碎的动物尸体及其废弃物，以病死猪化制处理为例：物料中肉块所占体积比约为70%，骨头约占20%，其余为毛发、粪便、血水等杂物。

化制物料组分复杂，黏度较大，流动性差，对搅拌器的强度和刚度要求较高，为优化改进型化制搅拌器工作参数，本文对化制容器内搅拌过程物料流动模型与特性进行分析。

由于双螺旋搅拌器搅拌叶片由旋向相反大小螺带组成，在其搅拌作用下，物料的的流动趋势可分解为轴向流、径向流和环向流，如图4所示。

(a)主视图(b)侧视图
图4搅拌过程中物料的流动模型
Fig.4Flow model of material during mixing
1.层流区域
2.湍流区域
3.轴向流
4.径向流
5.环向流
物料升温过程中，不同液面高度的物料温度有一定差别，径向流可以促进不同液面高度的物料相互混合流动，加快整体物料的升温速度。

轴向流是指物料从化制罐的一端向另一端的流动，由于搅拌器大小螺带的表面积不同，大螺旋带输送物料的能力比小螺旋带强，导致物料在大螺旋带的推动作用下不断向化制罐的一端积聚，部分物料未能及时被小螺旋带输送，继续在挤压作用下，借助自身的流动性向另一端流动，导致物料积聚部位的搅拌器和罐体承受较大压力。

轴向流动靠近化制罐内壁和搅拌轴，流动性较差，可视为层流区域，对物料的升温速度影响较小，因此应尽量减小物料的轴向流动趋势。

分析搅拌器的结构特点可得：径向流对搅拌器转速的敏感性比轴向流要强，当搅拌器转速提升后，径向流流速的提升较轴向流明显。

由于径向流动区域位于罐体内壁部位和搅拌轴周围的两个轴向流动区域之间，根据伯努利原理搅拌器转速增加，引起两个轴向流动区域速度较小，压强较大，而径向流动区域速度较大，压强较小，轴向流动区域的部分物料会有向径向流动区域运动的趋势，如图5所示。

因此可推断提升搅拌器转速可以削弱轴向流，增强径向流，即增强涡旋
现象。

图5提高转速后物料的流动趋势变化
Fig.5Change in the flow tendency of the material
after increasing the speed
综上所述，轴向流对物料升温速度影响不大，且会对化制罐和搅拌器产生额外应力，而径向流可以有效改善物料流动性，加快物料升温速度。

以本文改进双螺带搅拌器为例，其转速从10m/s变为20m/s，大螺带外沿物料流动速度由
0.362
m/s增大为0.905m/s，物料流动效果提升了250%。

3基于ANSYS的化制过程物料流动流固耦合数值仿真试验
基于以上结论，本文建立化制过程物料流动模型，采用ANSYS进行化制过程物料流动流固耦合数值仿真试验。

3.1物料参数测定
将猪肉块（约5cm×5cm×5cm）放在高温高压容器中，使用集热式恒温加热磁力搅拌器，以导热油为传热介质，对高温高压容器进行加热，保证容器内的最终温度达到150℃，分别测量加热时间为0min、30min、60min时物料的近似黏度，模拟化制处理前中后期物料的状态，对三个时间点的物料流动状态进行分析。

取化制过程中的三个时间点，对物料黏度进行测定，不同加热时间的物料状态如图6所示。

（a）t=0min（b）t=30min（c）t=60min
图6不同加热时间的物料状态
Fig.6Material status for different heating times
使用黏度计测得加热时间t 分别为0min 、30min 、60min 时，物料近似黏度分别为83.2Pa·s 、18.7Pa•s 、7.4Pa•s 。

可见化制处理前期的物料近似黏度较高，随着化制处理的进行，物料中的固态油脂在高温条件下变为液态油脂，物料的近似黏度迅速下降。

显然搅拌器在化制处理的前期所受的载荷最大，因此只需校核最大黏度条件下即化制过程开始时的搅拌器强度。

3.2
ANSYS 流固耦合数值建模
新型高压化制设备的化制罐和搅拌器的结构3D 模型如图7所示，本文针对部分结构特征进行简化，同时由于搅拌器搅拌物料时，物料对搅拌器叶片的阻力较大，导致搅拌器发生一定极小程度的变形，对物料流动状态的影响可以忽略，本文采用ANSYS Workbench 仿真软件建立单向流固耦合分析有限元网格模型，单元数量为622247，节点数量为113448，平均网格偏斜率为0.23，
网格质量良好。

(a)
模型图
(b)网格图
图7
化制罐及搅拌过程模型Fig.7
Finite element model of chemical tank and
mixing process
3.3物料流动性分析
使用CFD-post 模块进行结果后处理，创建yz 面
和xy 面显示物料在流域纵截面和横截面上的速度分布(图8)。

除去封头顶部、罐体内表面以及搅拌轴附近的少量物料，其他物料流动速度均达到0.362m/s 以上，且左右两端的物料流动状态基本一致，因此化制处理前期，在双螺带搅拌器搅拌作用下，罐内物料的流动性良好。

试验表明：随着化制处理的进行，物料黏度逐渐降低，物料流动性逐渐提高。

在螺带式搅拌器搅拌作用下，物料流动的速度矢量图、物料速度矢量在xy 和yz 坐标面上的分量如图9所示，由图可知物料高速流动区域主要分布在大螺带附近，螺带附近物料的速度方向与螺旋带螺旋线的切线方向相同。

速度方向与
螺旋带螺旋线的切线方向相同。

图9（a ）为物料速度矢量在yz 坐标面上的分量，物料在yz 坐标面上的速度分量均小于0.3m/s ，且在大螺旋带和小螺旋带之间存在多处涡旋现象，有利于物料的混合搅拌。

图9（b ）为物料速度矢量在xy 坐标面上的分量，物料在yz 坐标面上的速度分量分布不均匀性较强，平均流动速度较大。

因此，化制罐内的物料在搅拌器搅拌作用下同时存在轴向流、径向流和环向流，轴向流和径向流的流动速度相对较小，环向流流速较大，整体搅拌效果良好。

0.500 1.000 m 0.2500.750
Velocity /(m·s -1)
1.8091.6281.4471.2661.0850.9050.7240.5430.3620.1810.000
yz
Velocity /(m·s -1)1.8091.6281.4471.2661.0850.9050.7240.5430.3620.1810.000
xy
0.500 1.000 m
(a)流域纵截面速度云图
(b)流域横截面速度云图图8物料流动速度分布云图
Fig.8
Material flow velocity distribution
0.500 1.000 m 0.2500.750
Velocity /(m·s -1)1.8091.6281.4471.2661.0850.9050.7240.5430.3620.1810.000
yz
Velocity /(m·s -1)
1.8091.6281.4471.2661.0850.9050.7240.5430.3620.1810.000
xy
0.500 1.000 m 0.250
0.750
(a)yz 面上速度分量
(b)xy 面上的速度分量图9物料流动速度矢量图Fig.9
Material flow speed vector
0.500 1.000(m)0.2500.750
Velocity / m·s -1
1.8091.6281.4471.2661.0850.9050.7240.5430.3620.1810.000
xy
Velocity / m·s -1 1.8091.6281.4471.2661.0850.9050.7240.5430.362
0.181
0.000xy
00.500 1.000(m)0.2500.750
(a)n =10r/min
(c)n =30r/min
(b)n =20r/min
(d)n =40r/min
图10不同转速下物料流动速度矢量图Fig.10
Material flow speed vector at different speeds
3.4不同转速下的物料流动性分析
为了验证前文中所提出的提高搅拌器转速可以削
弱轴向流，增强径向流的推论，选取4组不同的转速值10r/min 、20r/min 、30r/min 、40r/min 分别进行流场分析，得到不同搅拌轴转速下的物料流动速度矢量在xy
和yz 坐标面上的分量，如图10所示。

不同转速下物料流动速度矢量图对比，当搅拌器转速为10r/min 时物料的流动趋势主要为轴向流动和环向流动，径向流动趋势较小，涡旋现象不明显。

随着物料的径向流动趋势不断增大，涡旋现象逐渐增强。

从数值仿真看：在搅拌器转速10r/min 、20r/min 、30r/min 、40r/min 增加过程中，轴向流动平均速度分别达到0.174m/s ，0.362m/s ，0.645m/s ，1.158m/s ，而径向截面最大速度分别达到0.221m/s ，0.398m/s ，0.612m/s ，0.937m/s 。

以搅拌器转轴为中心轴，不同半径圆柱面区域的物料径向流动平均速度变化情况如图11所示，由图可
以看出：不同半径圆柱面区域的物料径向流动平均速度变化情况，随着搅拌器转速的提高，物料的径向截面流动速度不断增大，且涡旋区域的影响范围不断扩大，有利于不同层面的物料进行混合。

搅拌器转中心径向575mm 处物料流动最大速度分别为：0.362m/s ，0.905m/s ，1.451m/s ，2.123m/s ，显著改善了物料的流动性，提升搅拌器搅拌效果。

E / m ·s -1
/mm
图11不同位置的物料径向流动速度变化曲线Fig.11
Radial flow velocity curve of materials at
different positions
4试验
根据仿真优化分析结果，进行改进搅拌器温升试
验（图12），通过监控界面观察并记录化制罐内温度值。

(a)实物图(b)显示界面图
#
K NJO
(c)温度变化图图12验证试验Fig.12
Verification test
试验结果表明，搅拌轴的搅拌速度不同，物料的升温和降温速度也不相同，当搅拌轴转速为5r/min 时，物料从22℃升至140℃用了约50min 时间，而搅拌轴转速为20r/min 时同样的过程只用了36min ，可见随着搅拌轴转速的提高，物料的升温和降温速度也有所加快。

同时与单螺带搅拌器化制过程温升对比，双螺带搅拌器温升时间缩短，证明其搅拌效果和化制罐内物料流动性能得到了提升。

5结论
1）提出一种改进病死动物化制设备的新型双螺
带式搅拌器结构，较好的提升化制过程中的物料搅拌效果，相比单螺带搅拌器提升物料流动性能20%。

2）通过化制罐搅拌过程中物料流动性流固耦合数值仿真分析试验，研究10r/min 、20r/min 、30r/min 、
40r/min 不同搅拌速度下化制容器内的物料流动速度与分布规律，仿真及试验表明：随着搅拌器转速的提高，物料的径向流动速度不断增大，且涡旋区域的影响范围不断扩大，有利于不同层面的物料进行混合，提高物料的流动性，改善搅拌器搅拌效果
参考文献
[1]戴冕.病死动物无害化处理工艺及关键设备研究[D].武汉：
武汉理工大学,2018.
Dai Mian.Study on the harmless treatment process and key equipment of dead animals [D].Wuhan:Wuhan University of Technology,2018.
[2]Gwyther C L,Williams A P,Golyshin P N,et al.The environmental and biosecurity characteristics of livestock carcass disposal method:A review [J].Waste Management,2011,31(4):767-778.
[3]Kalbasiashtari A,Schutz M M,Auvermann B W.Carcass rendering systems for farm mortalities:A review [J].Journal of Environmental Engineering &Science,2008,7(3):199-211.[4]王德胜.病死动物干化化制工艺研究及优化[D].杭州：浙江工业大学,2019.
Wang Desheng.Research and optimization of dry rendering process of diseased and dead animals [D].Hangzhou:Zhejiang University of Technology,2019.
[5]陶兰兰,季承玺.关于化工设备搅拌器的选择及设计研究[J].
化工管理,2018,490(19):154.
[6]陆旭.搅拌设备空心搅拌轴的设计开发[J].石油和化工设备,2016(11):35-36.
[7]王长利,王金茹,刘春茹,等.病死动物干化化制法无害化处理装置及其处理方法[P].中国:104368588.2014-10-13.[8]张鑫.基于流固耦合的低比转速离心泵流动诱导振动研究[D].杭州：浙江理工大学,2016.
Zhang Xin.Investigations on fluid⁃structure interaction of flow⁃induced vibration of low specific speed centrifugal pumps [D].Hangzhou:Zhejiang Sci⁃Tech University,2016.
[9]苏腾龙.柔性桨流场特性的流固耦合数值模拟和试验研究[D].济南:山东大学,2018.
Su Tenglong.Fluid⁃slid interaction simulation and experiment study on the flow field characteristics of flexible ⁃blade impel⁃ler [D].Jinan:Shandong University ，2018.
[10]Khan H M.Traditional methods of carcass disposal:A
review[J].Journal of Dairy Veterinary&Animal Research, 2017,5(1):15-16.
[11]Xi Jialin,Dong Hongmin,Zhu Zhiping,et al.Research
progress on composting treatment of dead livestock[J].
Journal of Agricultural Science&Technology,2010. [12]Hou.Numerical methods for fluid⁃structure interaction:A
review[J].Communications in Computational Physics, 2012,12(2):337-377.
[13]胡效东,王超,王灏,等.基于流-固耦合理论的搅拌反应
器机械特性[J].济南大学学报(自然科学版),2017,31(2): 150-158.
Hu Xiaodong,Wang Chao,Wang Hao,et al.Mechanical properties of stirred reactors based on Fluid⁃Structure interaction theory[J].Journal of University of Jinan(Science and Technology),2017,31(2):150-158.
[14]蔡琳,江寒冰.搅拌机叶片的流固耦合计算机仿真分析[J].
制造业自动化,2012,34(2):69-71.
Cai Lin,Jiang Hanbing.Analysis of FSI for mixer blades based on computer simulation[J].Manufacturing Automa⁃tion,2012,34(2):69-71.
[15]巩亚东,高志州,温雪龙.一种发酵机搅拌装置的流固耦
合有限元仿真研究[J].东北大学学报(自然科学版),2017, 38(11):1584-1589.[16]胡效东,王海月,刘馥瑜,等.流固耦合的搅拌器振动特性[J].
机械设计与研究,2017,33(5):63-67.
Hu Xiaodong,Wang Haiyue，Liu Fuyu,et al.Study on vibration characteristics for agitator with fluid⁃structure coupling[J].Machine Design and Research,2017,33(5): 63-67.
[17]包凯,刘清友.基于workbench流固耦合节流阀数值分析与
优化设计[J].机械设计与研究,2014,30(4):116-119.
Bao Kai,Liu Qingyou.The numerical analysis and optimization design for throttle valves based on fluid⁃soild coupling workbench[J].Machine Design and Research, 2014,30(4):116-119.
[18]刘厚林,徐欢,吴贤芳,等.基于流固耦合的导叶式离心泵
强度分析[J].振动与冲击,2013,32(12):27-30.
[19]刘欣.柔性Rushton搅拌桨的振动特性研究[D].济南：山东
大学,2020.
Liu Xin.Study on vibration characteristics of flexible⁃blade rushton impeller[D].Jinan:Shandong University,2020. [20]王海月.搅拌反应器振动模态分析与试验研究[D].青岛：
山东科技大学,2018.
Wang Haiyue.Vibration modal analysis and experimental study of stirrde reactor[D].Qingdao:Shandong University of Science and Technology,2018.
Numerical simulation and optimization of fluid-solid coupling in the mixing process of
sick and dead animals
Yu Guohong1,Xue Xianglei1,Zheng Hang1,Niu Ruizheng2,Jiang Jiandong2
(1.Institute of Agricultural Equipment,Zhejiang Academy of Agricultural Sciences,Hangzhou,310021,China;
2.Key Laboratory of Special Purpose Equipment and Advanced Manufacturing Technology,Ministry of
Education,Zhejiang University of Technology,Hangzhou,310014,China)
Abstract:To solve the problems of complex and changeable material components and uneven flow in the drying process of sick and dead animals,a new type of double ribbon agitator was proposed.The fluid⁃solid coupling of the material in the chemical tank and the agitator was carried out using ANSYS Software.Modeling and simulation analysis was used to investigatethe material flow rate and distribution mechanism in the chemical vessel under different stirring speeds and combine the material flow line distribution in the drying vesselso as to optimize the structure of the agitator.Experimental and numerical optimization results showed that the maximum deformation of the double ribbon agitator was0.21mm at a stirring speed of10r/min,and flow velocity of more than80%of the materials in the chemical system reached0.419m/pared with the single ribbon agitator,the material flow performance wasincreased by20%,and the efficiency and quality of the drying process of sick and dead animals were improved.
Keywords:chemical process;beater;fluid⁃structure coupling;flow characteristics;numerical simulation。