米糠和碎米的接触参数测量与离散元仿真标定

doi:10.13304/j.nykjdb.2022.0835

中国农业科技导报 ›› 2024, Vol. 26 ›› Issue (2): 127-136.DOI: 10.13304/j.nykjdb.2022.0835

• 智慧农业农机装备 • 上一篇

米糠和碎米的接触参数测量与离散元仿真标定

陈林¹(), 余南辉¹, 王立宗¹, 范吉军²(), 雷港¹, 刘晓鹏³, 周龙¹, 周劲²

^1.武汉轻工大学机械工程学院，武汉 430023
^2.武汉轻工大学电气与电子工程学院，武汉 430023
^3.武汉轻工大学动物科学与营养工程学院，武汉 430023

收稿日期:2022-10-02 接受日期:2023-04-22 出版日期:2024-02-15 发布日期:2024-02-04
通讯作者: 范吉军
作者简介:陈林 E-mail：1575810724@qq.com；
基金资助:
湖北省重点研发计划项目(2022BBA0047);2021年湖北省揭榜制粮食科技项目(鄂粮函〔2021〕18号);2022年湖北省揭榜制粮食科技计划项目(鄂粮函〔2022〕21号)

Measurement of Contact Parameters and Discrete Element Simulation Calibration of Rice Bran and Broken Rice

Lin CHEN¹(), Nanhui YU¹, Lizong WANG¹, Jijun FAN²(), Gang LEI¹, Xiaopeng LIU³, Long ZHOU¹, Jin ZHOU²

^1.School of Mechanical Engineering，Wuhan Polytechnic University，Wuhan 430023，China
^2.School of Electrical and Electronic Engineering，Wuhan Polytechnic University，Wuhan 430023，China
^3.School of Animal Science and Nutritional Engineering，Wuhan Polytechnic University，Wuhan 430023，China

Received:2022-10-02 Accepted:2023-04-22 Online:2024-02-15 Published:2024-02-04
Contact: Jijun FAN

摘要/Abstract

摘要：

为确定米糠和碎米分离装备离散元仿真所需的接触参数，利用台架试验与仿真试验对米糠、碎米颗粒接触参数进行测量和标定。分别采用自由落体碰撞法、碰撞扩散法、斜面滑动法和斜面滚动法标定了碎米与亚克力间的碰撞恢复系数为0.45、静摩擦系数为0.38、滚动摩擦系数为0.17，米糠与亚克力间的碰撞恢复系数为0.15，静摩擦系数为0.76。基于休止角测量试验，结合MATLAB软件图像处理技术获取了颗粒堆休止角；通过最陡爬坡试验和二次正交旋转组合试验，建立了影响因素与休止角之间的二次回归模型，并以休止角的相对误差最小为约束条件，对回归方程寻优求解，得到了最佳参数组合，其中，碎米间静摩擦系数0.60、滚动摩擦系数0.40、碰撞恢复系数0.61，米糠间静摩擦系数0.81、滚动摩擦系数0.18、碰撞恢复系数0.15。研究结果为米糠和碎米分离装备离散元仿真参数设置提供了参考。

关键词: 米糠, 碎米, 离散元, 接触参数, 标定

Abstract:

In order to determine the contact parameters required for the discrete element simulation of rice bran and broken rice separation equipment， the contact parameters of rice bran and broken rice particles were measured and calibrated through bench test and simulation test. By free-fall impact method， impact diffusion method， inclined plane sliding method and inclined plane rolling method respectively， it was calibrated that the collision recovery coefficient between broken rice and acrylic was 0.45， static friction coefficient was 0.38， rolling friction coefficient was 0.17， the collision recovery coefficient between rice bran and acrylic was 0.15， static friction coefficient was 0.76. Based on the measurement test of repose angle and MATLAB image processing technology， the repose angle of particle reactor was obtained. Through the steepest climb test and quadratic orthogonal rotation combination test， a quadratic regression model between the influencing factors and repose angle was established， and the minimum relative error of the repose angle was taken as the constraint condition to optimize the regression equation. The optimal parameter combination was obtained： static friction coefficient 0.60， rolling friction coefficient 0.40， collision recovery coefficient 0.61 for broken rice； static friction coefficient was 0.81， the rolling friction coefficient was 0.18， and the collision recovery coefficient was 0.15 for rice bran. The results provided reference for the discrete element simulation parameter setting of rice bran and broken rice separation equipment.

Key words: rice bran, broken rice, discrete element method, contact parameters, calibration of parameters

中图分类号:

S226.5

陈林, 余南辉, 王立宗, 范吉军, 雷港, 刘晓鹏, 周龙, 周劲. 米糠和碎米的接触参数测量与离散元仿真标定[J]. 中国农业科技导报, 2024, 26(2): 127-136.

Lin CHEN, Nanhui YU, Lizong WANG, Jijun FAN, Gang LEI, Xiaopeng LIU, Long ZHOU, Jin ZHOU. Measurement of Contact Parameters and Discrete Element Simulation Calibration of Rice Bran and Broken Rice[J]. Journal of Agricultural Science and Technology, 2024, 26(2): 127-136.

图/表 19

表1 材料本征参数

Table 1 Material parameter

材料

Material

密度Density/（kg·m^-3）

泊松比

Poisson’s ratio

图1 米糠和碎米等效离散元模型

Fig. 1 Equivalent discrete element model of rice bran and broken rice

图2 碎米碰撞恢复系数测量原理及试验

Fig. 2 Measurement principle and test of collision recovery coefficient of broken rice

图3 米糠与亚克力间碰撞恢复系数测量试验

Fig. 3 Measurement test of collision recovery coefficient between rice bran and acryli

图4 静摩擦系数测量原理与仿真试验

Fig. 4 Measurement principle and simulation test of static friction coefficient

图5 滚动摩擦系数测量试验

Fig. 5 Test of rolling friction coefficient

图6 休止角测量试验

Fig. 6 Test of angle of repose

图7 MATLAB软件图像处理

Fig. 7 Image processing by MATLAB software

图8 碎米碰撞恢复系数与反弹高度拟合曲线

Fig. 8 Fitting curve of collision recovery coefficient and rebound height of broken rice

图9 米糠碰撞恢复系数与扩散直径拟合曲线

Fig. 9 Fitting curve of collision recovery coefficient and diffusion diameter of rice bran

图10 碎米静摩擦系数与斜面倾角拟合曲线

Fig. 10 Fitting curve of the static friction coefficient and the inclination angle of the inclined plate of broken rice

图11 碎米滚动摩擦系数与水平滚动距离拟合曲线

Fig. 11 Fitting curve of rolling friction coefficient and horizontal rolling distance of broken rice

表2 最陡爬坡试验设计与结果

Table 2 Steepest climb test design and results

序号Number	试验因素Test factor			试验结果Test result
序号Number	x₅	x₆	x₇	θ_y ´/（°）	$σ$ /%
1	0.35	0.10	0.25	20.5	41.3
2	0.45	0.20	0.35	25.2	27.8
3	0.55	0.30	0.45	27.8	20.2
4	0.65	0.40	0.55	33.1	5.2
5	0.75	0.50	0.65	41.3	18.3
6	0.85	0.60	0.75	43.3	24.1

表2 最陡爬坡试验设计与结果

Table 2 Steepest climb test design and results

序号Number	试验因素Test factor			试验结果Test result
序号Number	x₅	x₆	x₇	θ_y ´/（°）	$σ$ /%
1	0.35	0.10	0.25	20.5	41.3
2	0.45	0.20	0.35	25.2	27.8
3	0.55	0.30	0.45	27.8	20.2
4	0.65	0.40	0.55	33.1	5.2
5	0.75	0.50	0.65	41.3	18.3
6	0.85	0.60	0.75	43.3	24.1

表3 仿真试验因素编码

Table 3 Simulation test factor coding

编码Code	试验因素Test factor
编码Code	x₅	x₆	x₇
-1.682	0.48	0.23	0.38
-1	0.55	0.30	0.45
0	0.65	0.40	0.55
1	0.75	0.50	0.65
1.682	0.82	0.57	0.72

表4 二次正交旋转组合试验设计与结果

Table 4 Experimental design and results of quadratic orthogonal rotation combination

序号Number	试验因素Test factor			Y/%
序号Number	A	B	C	Y/%
1	-1	-1	-1	20.2
2	1	-1	-1	11.2
3	-1	1	-1	8.5
4	1	1	-1	20.6
5	-1	-1	1	15.2
6	1	-1	1	13.7
7	-1	1	1	5.5
8	1	1	1	18.3
9	-1.682	0	0	9.9
10	1.682	0	0	15.3
11	0	-1.682	0	26.0
12	0	1.682	0	24.3
13	0	0	-1.682	3.6
14	0	0	1.682	0.6
15	0	0	0	2.6
16	0	0	0	3.0
17	0	0	0	2.5
18	0	0	0	3.6
19	0	0	0	2.5

表5 二次回归模型方差分析

Table 5 ANOVA of quadratic regression model

方差来源 Soruce of variation	平方和 Mean square	自由度 Freedom	F值 F value	P值 P value
R² =0.995 3	R_Adj²=0.990 7	Predicted R² =0.964 3	CV=7.28%	Adep precision=38.901 0
模型Model	1226.53	9	213.39	<0.000 1^**
A	45.33	1	70.98	<0.000 1^**
B	5.75	1	9.00	0.015 0^*
C	9.59	1	15.02	0.003 8^*
AB	169.28	1	265.07	<0.000 1^**
AC	11.52	1	18.04	0.002 2^**
BC	0.245 0	1	0.383 6	0.551 0
A²	170.84	1	267.51	<0.000 1^**
B²	868.12	1	1 359.34	<0.000 1^**
C²	0.4161	1	0.651 6	0.440 4
残差Residual	5.75	9
失拟项Lack of fit	4.86	5	4.35	0.089 5
纯误差Pure error	0.892 0	4
总和Sum	1 232.28	18

表6 二次优化后回归模型方差分析

Table 6 ANOVA of regression model after second optimization

方差来源 Soruce of variation	平方和 Mean square	自由度 Freedom	F值 F value	P值 P value
R² =0.998 8	R_Adj²=0.997 6	Predicted R² =0.985 3	CV=8.38%	Adep precision =42.844 5
模型Model	1 200.53	7	205.72	<0.000 1**
A	40.37	1	48.43	<0.000 1**
B	7.71	1	9.24	0.011 2*
C	12.08	1	14.49	0.002 9**
AB	156.65	1	187.90	<0.000 1**
AC	8.40	1	10.08	0.008 8**
A²	173.99	1	208.71	<0.000 1**
B²	885.88	1	1 062.64	<0.000 1**
残差Residual	9.17	11
失拟项Lack of fit	8.28	7	5.30	0.063 1
纯误差Pure error	0.892 0	4
总和Sum	1 209.70	18

图12 碎米堆积角

Fig. 12 Broken rice stacking angle

图13 米糠堆积角

Fig. 13 Rice bran stacking angle

参考文献 29

1	赵志浩,邓媛元,魏振承,等.大米适度加工和副产物综合利用现状及展望[J].广东农业科学,2020,47(11):144-152.
	ZHAO Z H, DENG Y Y, WEI Z C, et al.. Current status and prospect of appropriate processing of rice and comprehensive utilization of by-products [J]. Guangdong Agric. Sci., 2020,47(11):144-152.
2	邹昱成,鞠兴荣,何荣,等.食品配料米糠的研究进展[J].中国粮油学报,2018,33(7):123-129.
	ZOU Y C, JU X R, HE R, et al.. Research progress for rice bran as a food ingredient [J]. J. Chin. Cereals Oils Assoc., 2018, 33(7):123-129.
3	李永祥,李阳,王明旭,等.基于EDEM的大产量原粮振动清理筛筛分效果分析[J].中国油脂,2018,43(5):157-160.
	LI Y X, LI Y, WANG M X, et al.. Screening effect of large production raw grain vibration cleaning sieve based on EDEM [J]. China Oils Fats, 2018, 43(5):157-160.
4	党贺,赵宝生,马学东,等.CFD-DEM耦合的旋风式谷物清选装置的模拟研究[J].中国农机化学报,2020,41(4):86-91.
	DANG H, ZHAO B S, MA X D, et al.. Simulating research on cyclone grain cleaning device based on CFD-DEM coupling [J]. J. Chin. Agric. Mechan., 2020, 41(4):86-91.
5	王立宗,余南辉,范吉军,等.上吸风碾米机环形风道内颗粒运动特性[J].包装与食品机械,2022,40(2):25-30.
	WANG L Z, YU N H, FAN J J, et al.. Motion characteristics of particles in annular duct of top-suction rice mill [J]. Packag. Food Mach., 2022, 40(2):25-30.
6	HU M J, et al.. Measurement and calibration of the discrete element parameters of coated delinted cotton seeds [J]. Agriculture, 2022, 12(2) :286-286.
7	LI H, ZENG R, NIU Z, et al.. A calibration method for contact parameters of maize kernels based on the discrete element method [J/OL]. Agriculture, 2022, 12(5): 664 [2023-10-26]. .
8	刘羽平,张拓,刘妤.稻谷颗粒模型离散元接触参数标定与试验[J].中国农业科技导报,2019,21(11):70-76.
	LIU Y P, ZHANG T, LIU Y. Calibration and experiment of contact parameters of rice grain based on discrete element method [J]. J. Agric. Sci. Technol., 2019,21(11):70-76.
9	闫建伟,魏松,胡冬军,等.白萝卜种子颗粒模型离散元接触参数标定与试验[J].中国农业科技导报,2022,24(5):119-128.
	YAN J W, WEI S, HU D J, et al.. Parameter calibration of radish seeds with different filling particle radius by DEM [J]. J. Agric. Sci. Technol., 2022,24(5):119-128.
10	李飞翔,王鹏,王云飞,等.基于堆积试验的玉米包衣种子离散元参数标定[J].中国农业科技导报,2022,24(7):97-107.
	LI F X, WANG P, WANG Y F, et al.. Calibration of discrete element parameters of corn coated seeds based on stacking test [J]. J. Agric. Sci. Technol., 2022,24(7):97-107.
11	田剑锋,石林榕,杨小平,等.党参种子的离散元仿真参数标定与试验验证[J].干旱地区农业研究,2022,40(2):240-249.
	TIAN J F, SHI L R, YANG X P, et al.. Calibration and experimental validation of discrete element simulation parameters for Codonopsis pilosula seed [J]. Agric. Res. Arid Areas, 2022, 40(2):240-249.
12	张春,杜文亮,陈震,等.荞麦米筛分物料接触参数测量与离散元仿真标定[J].农机化研究,2019,41(1):46-51.
	ZHANG C, DU W L, CHEN Z, et al.. The measurement of contact parameters of buckwheat rice screening material and discrete element simulation calibration [J]. J. Agric. Mechan. Res., 2019, 41(1):46-51.
13	刘凡一,张舰,李博,等.基于堆积试验的小麦离散元参数分析及标定[J].农业工程学报,2016,32(12):247-253.
	LIU F Y, ZHANG J, LI B, et al.. Calibration of parameters of wheat required in discrete element method simulation based on repose angle of particle heap [J]. Trans. Chin. Soc. Agric. Eng., 2016, 32(12):247-253.
14	朱康,惠延波,周颖,等.基于X射线断层扫描的小麦籽粒模型构建及离散元参数标定[J].中国粮油学报,2022,37(9):218-223.
	ZHU K, HUI Y B, ZHOU Y, et al.. Gonstruction of wheat grain model and calibration of discrete element parameters based on X-ray tomography [J]. J. Chin. Cereals Oils Assoc., 2022,37(9):218-223.
15	Balevičius R, Sielamowicz I, Mróz Z, et al.. Investigation of wall stress and outflow rate in a flat-bottomed bin: a comparison of the DEM model results with the experimental measurements [J]. Powder Technol., 2011, 214(3): 322-336.
16	于庆旭,刘燕,陈小兵,等.基于离散元的三七种子仿真参数标定与试验[J].农业机械学报,2020,51(2):123-132.
	YU Q X, LIU Y, CHEN X B, et al.. Calibration and experiment of simulation parameters for Panax notoginseng seeds based on DEM [J]. Trans. Chin. Soc. Agric. Mach., 2020, 51(2):123-132.
17	吴孟宸,丛锦玲,闫琴,等.花生种子颗粒离散元仿真参数标定与试验[J].农业工程学报,2020,36(23):30-38.
	WU M C, CONG J L, YAN Q, et al.. Calibration andexperiments for discrete element simulation parameters ofpeanut seed particles [J]. Trans. Chin. Soc. Agric. Eng., 2020, 36(23):30-38.
18	王立宗.上吸风环形风道碾米机的设计与实验研究[D].武汉:武汉轻工大学,2021.
	WANG L Z. Design and experimental study of up-suction rice mill with annular air duct [D]. Wuhan: Wuhan Polytechnic University, 2021.
19	王立军,刘天华,冯鑫,等.农业和食品领域中颗粒碰撞恢复系数的研究进展[J].农业工程学报,2021,37(20):313-322.
	WANG L J, LIU T H, FENG X, et al.. Research progress of the restitution coefficients of collision of particles in agricultural and food fields [J]. Trans. Chin. Soc. Agric. Eng., 2021,37(20):313-322.
20	王成军,李耀明,马履中,等.小麦籽粒碰撞模型中恢复系数的测定[J].农业工程学报,2012,28(11):274-278.
	WANG C J, LI Y M, MA L Z, et al.. Experimental study on measurement of restitution coefficient of wheat seeds in collision models [J]. Trans. Chin. Soc. Agric. Eng., 2012, 28(11): 274-278.
21	冯斌,孙伟,石林榕,等.收获期马铃薯块茎碰撞恢复系数测定与影响因素分析[J].农业工程学报,2017,33(13):50-57.
	FENG B, SUN W, SHI L R, et al.. Determination of restitution coefficient of potato tubers collision in harvest and analysis of its influence factors [J]. Trans. Chin. Soc. Agric. Eng., 2017, 33(13): 50-57.
22	张胜伟,张瑞雨,陈天佑,等.绿豆种子离散元仿真参数标定与排种试验[J].农业机械学报,2022,53(03):71-79.
	ZHANG S W, ZHANG R Y, CHEN T Y, et al.. Discretc element simulation paramcter caliraion and seelingverification test of mung bean seeds [J]. Trans. Chin. Soc. Agric. Eng., 2022,53(3):71 -79.
23	JÓZEF H, MICHAŁ B, RAFAŁ M, et al.. Determination of the restitution coefficient of seeds and coefficients of visco-elastic Hertz contact models for DEM simulations [J]. Biosys. Eng., 2017, 161:106-119.
24	FANG W Q, WANG X Z, et al.. Review of material parameter calibration method [J/OL]. Agriculture, 2022, 12(5): 706 [2023-10-26]. .
25	韩燕龙,贾富国,唐玉荣,等.颗粒滚动摩擦系数对堆积特性的影响[J].物理学报,2014,63(17):173-179.
	HAN Y L, JIA F G, TANG Y R, et al.. Influence of granular coefficient of rolling friction on accumulation characteristics [J]. Acta Physiol. Sin., 2014,63(17):173-179.
26	贾富国,韩燕龙,刘扬,等.稻谷颗粒物料堆积角模拟预测方法[J].农业工程学报,2014,(11) :254-260.
	JIA F G, HAN Y L, LIU Y, et al.. Simulation predictionmethod of repose angle for rice particle materials [J]. Trans. Chin. Soc. Agric. Eng., 2014(11):254-260.
27	GUO Z, CHEN X, LIU H, et al.. Theoretical and experimental investigation on angle of repose of biomass-coalblends [J]. Fuel, 2014, 116(15):131-139.
28	吴爱祥,孙业志,刘湘平.散体动力学理论及其应用[M].北京:冶金工业出版社, 2002:1-364.
29	WANG J, XU C, XU Y, et al.. InfluencingFactors analysis and simulation calibration of restitution coefficient of rice grain [J/OL]. Appl. Sci., 2021, 11(13):5884 [2023-10-26]. .

[1]	王洪波, 樊志鹏, 乌兰图雅, 王春光, 马哲. 揉碎玉米秸秆螺旋输送仿真离散元模型参数标定[J]. 中国农业科技导报, 2023, 25(3): 96-106.
[2]	于淼, 周海宾, 丁京涛, 程红胜, 沈玉君, 范盛远, 张曦, 王健, 徐鹏翔, 程琼仪. 基于EDEM的餐厨垃圾组成颗粒间接触参数标定[J]. 中国农业科技导报, 2023, 25(12): 111-120.
[3]	马紫涛, 赵智豪, 全伟, 石方刚, 高晨, 吴明亮. 基于EDEM的水稻残茬秸秆离散元仿真参数标定[J]. 中国农业科技导报, 2023, 25(11): 103-113.
[4]	单发科, 康朔, 朱建锡, 王永维, 王俊. 基于EDEM的粉垄和旋耕作业混肥效果研究[J]. 中国农业科技导报, 2023, 25(11): 90-102.
[5]	周婷, 孙松林, 朱海英, 彭才望. 含水率对黑水虻生物转化猪粪有机肥黏结流动的影响[J]. 中国农业科技导报, 2023, 25(10): 126-136.
[6]	李飞翔, 王鹏, 王云飞, 葛越锋, 唐凯怿, 李得志. 基于堆积试验的玉米包衣种子离散元参数标定[J]. 中国农业科技导报, 2022, 24(7): 97-107.
[7]	宋世圣, 孙松林, 方芹, 彭才望, 周婷, 朱海英. 黑水虻生物转化餐厨垃圾有机肥离散元模型参数标定[J]. 中国农业科技导报, 2022, 24(6): 123-132.
[8]	闫建伟, 魏松, 胡冬军, 刘启合, 张富贵. 白萝卜种子颗粒模型离散元接触参数标定与试验[J]. 中国农业科技导报, 2022, 24(5): 119-128.
[9]	杨贵川, 张富贵, 郑乐, 王震, 孔曼曼, 章鑫鹏. 半夏块茎物理特性研究及离散元仿真参数标定[J]. 中国农业科技导报, 2022, 24(10): 99-108.
[10]	全伟, 吴明亮, 官春云, 罗海峰. 油菜钵苗移栽机成穴器外形优化试验研究[J]. 中国农业科技导报, 2021, 23(10): 97-106.
[11]	刘妤,刘羽平,张拓. 基于离散元与多体动力学的微耕机旋耕刀轴负荷分析[J]. 中国农业科技导报, 2020, 22(11): 79-86.
[12]	向伟1,2,吴明亮1*,吕江南2,马兰2,全伟1,刘佳杰2,肖乐3. 基于EDEM的油菜移栽成穴装置作业性能仿真与试验研究[J]. 中国农业科技导报, 2019, 21(7): 70-81.
[13]	刘羽平,张拓,刘妤*. 稻谷颗粒模型离散元接触参数标定与试验[J]. 中国农业科技导报, 2019, 21(11): 70-76.
[14]	. 利用信息技术改造传统农业——智能化农业信息技术应用示范工程[J]. , 2001, 3(2): 49-53.

米糠和碎米的接触参数测量与离散元仿真标定

Measurement of Contact Parameters and Discrete Element Simulation Calibration of Rice Bran and Broken Rice

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表 19

参考文献 29

相关文章 14

编辑推荐

Metrics