Parameter Calibration of Discrete Element Model for Simulation of Crushed Corn Stalk Screw Conveying

doi:10.13304/j.nykjdb.2022.0725

Abstract

Abstract:

In order to further reveal the conveying mechanism of crushed corn stalks feed in the screw conveyor and improve the accuracy of the parameters used in discrete element simulation research during the screw conveying process of crushed corn stalks， the crushed corn stalks were divided into straw core， straw leaf and straw skin. The discrete element simulation parameters were calibrated. First， the average and range of the intrinsic parameters and contact parameters of core， leaf and skin were obtained through physical tests， and the contact parameters with a relatively large range were screened by Plackett-Burman test using their significance. The test results showed that the factors that had significant effect on the simulated stacking angle were the core-leaf rolling friction coefficient， the skin-leaf rolling friction coefficient， and the leaf-leaf static friction coefficient. Then， the evaluation index is set as the three significant parameters of the relative error value between the test and the actual stacking angle to carry out the steepest climbing test， and the value range of the three significant parameters was optimized and narrowed， and finally the accumulation angle and the second-order regression model between the angle and the saliency parameter were established by Box-Behnken test. The target value of the accumulation angle was set to 30.4° measured by the physical test， and the optimal value of the saliency parameters were obtained by optimization， which the leaf-core rolling friction coefficient， skin-leaf rolling friction coefficient and leaf-leaf static friction coefficient were 0.325， 0.377 and 0.411. And the simulated stacking angle obtained when each parameter was optimal was 30.77°. In order to further test the accuracy of the parameter values， the t test showed P was more than 0.05， indicating that there was no significant difference in the value of the stacking angle between the simulation and the physical test， which could verify the reliability of the optimal parameter values. The results showed that it was feasible to use the above physical tests and optimization tests to measure and calibrate discrete element simulation parameters， and the calibrated parameters could provide reference for other simulation tests of crushing corn stalks.

Key words: crushed corn stalk, discrete element, stacking angle, parameter calibration

摘要：

为进一步揭示揉碎玉米秸秆饲料在螺旋输送器中的输送机理，提高揉碎玉米秸秆螺旋输送过程中离散元仿真研究所用参数的准确度，将揉碎玉米秸秆分为秸秆穰、秸秆叶、秸秆皮，采用物理试验和仿真优化设计相结合的方法对离散元仿真参数进行标定。首先，通过物理试验得到穰、叶、皮各项本征参数和接触参数的平均值和范围，将范围相对较大的接触参数通过Plackett-Burman试验利用其显著性进行筛选。结果表明，对仿真堆积角影响显著的因素有穰-叶滚动摩擦系数、皮-叶滚动摩擦系数、叶-叶静摩擦系数。对这3个显著性参数进行最陡爬坡试验，将其取值范围进行优化缩小，通过Box-Behnken试验建立堆积角与显著性参数的二阶回归模型，堆积角的目标值设置为物理试验测得的30.4°，寻优得到最佳的显著性参数的取值：叶-穰滚动摩擦系数为0.325，皮-叶滚动摩擦系数为0.377，叶-叶静摩擦系数为0.411，在此组合下仿真堆积角为30.77°。为进一步检验参数取值的准确性，通过t检验得到P>0.05，表明仿真与物理试验堆积角的值无明显差异，验证了最佳参数取值的可靠性。由此表明，应用上述各物理试验和优化试验来测定及标定离散元仿真参数是可行的，同时标定的参数为揉碎玉米秸秆的其他仿真试验提供了参考。

关键词: 揉碎玉米秸秆, 离散元, 堆积角, 参数标定

CLC Number:

TH132.1

Hongbo WANG, Zhipeng FAN, Wulantuya, Chunguang WANG, Zhe MA. Parameter Calibration of Discrete Element Model for Simulation of Crushed Corn Stalk Screw Conveying[J]. Journal of Agricultural Science and Technology, 2023, 25(3): 96-106.

王洪波, 樊志鹏, 乌兰图雅, 王春光, 马哲. 揉碎玉米秸秆螺旋输送仿真离散元模型参数标定[J]. 中国农业科技导报, 2023, 25(3): 96-106.

Figures/Tables 18

References 25

1	周腰华,周洋,赖晓璐,等.玉米秸秆综合利用技术模式研究[J].玉米科学, 2019,27(6):186-190.
	ZHOU Y H, ZHOU Y, LAI X L, et al.. Research on the technology model of comprehensive utilization of corn stalk [J]. Corn Sci., 2019, 27(6):186-190.
2	赵福迪,魏思明,黄文明.玉米秸秆机械化利用情况与发展途径分析[J].农机使用与维修, 2019(7):48-48.
3	乌兰图雅,青林,王春光.揉碎玉米秸秆螺旋-气力耦合输送装置设计[J].农业工程学报, 2019,35(6):29-38.
	WULANTUYA, QING L, WANG C G. Design of screw-pneumatic coupling conveying device for crushing corn stalks [J]. Trans. Chin. Soc. Agric. Eng., 2019, 35(6):29-38.
4	DHEERAJ M, ABHISHEK S, HARSH P, et al.. A review of granular flow in screw feeders and conveyors [J]. Powder Technol., 2020, 366:369-381.
5	扈伟昊,杨发展,赵国栋,等.基于离散元法的立式旋耕刀工作参数分析与优化[J].中国农机化学报, 2022,43(10):25-32, 41.
	HU W H, YANG F Z, ZHAO G D, et al.. Analysis and optimization of working parameters of vertical rotary blade based on discrete element method [J]. China J. Agric. Mach. Chem., 2022, 43(10):25-32, 41
6	刘元义,于圣洁,胥备,等.基于离散元的设施农业就地翻土犁的研究与试验[J/OL].吉林大学学报(工学版),2022 [2022-10-09]. .
	LIU Y Y, YU S J, XU B, et al. Research and test of on-site plow for facility agriculture based on discrete element [J/OL]. J. Jilin Univ. (Eng.), 2022 [2022-10-09]. .
7	李修银,廖敏,杨杰,等.羌活籽粒和珍珠岩的离散元参数标定及排种验证[J].中国农机化学报, 2022,43(8):40-46.
	LI X Y, LIAO M, YANG J, et al.. Discrete element parameter calibration and seed metering verification of Qianghuo grain and perlite [J]. China J. Agric. Mach. Chem., 2022, 43 (8):40-46.
8	陈晨,马学东,刘飞宇,等.偏心摆振下圆簸箕中谷物分离的离散元模拟研究[J].农机化研究,2023,45(5):13-19.
	CHEN C, MA X D, LIU F Y, et al.. Discrete element simulation study on grain separation in circular dustpan under eccentric sway [J]. Res. Agric. Mech., 2023, 45(5):13-19.
9	季雨.基于CFD-DEM螺旋输送机内颗粒流动特性数值模拟研究[D].大庆:东北石油大学, 2020.
	JI Y. Numerical simulation of particle flow characteristics in screw conveyor based on CFD-DEM [D]. Daqing: Northeast Petroleum University, 2020.
10	宁廷州,李偎,缑亚楠.基于EDEM的倾斜螺旋输送机输送能力研究[J].牡丹江师范学院学报(自然科学版),2019(4):38-41.
	NING T Z, LI C, JIN Y N. Research on conveying capacity of inclined screw conveyor based on EDEM [J]. J. Mudanjiang Norm. Univ. (Nat. Sci.), 2019(4):38-41.
11	张胜伟,张瑞雨,陈天佑,等.绿豆种子离散元仿真参数标定与排种试验[J].农业机械学报, 2022, 53(3):71-79.
	ZHANG S W, ZHANG R Y, CHEN T Y, et al.. The discrete element simulation parameter calibration and seeding test of mung bean seeds [J]. J. Agric. Mach., 2022, 53(3):71-79.
12	丁文波,朱继平,陈伟,等.基于的青稞接触参数仿真标定[J].中国农机化学报,2021,42(9):114-121.
	DING W B, ZHU J P, CHEN W, et al.. Simulation calibration of contact parameters of highland barley based on EDEM [J]. Chin. J. Agric. Mach., 2021, 42(9):114-121.
13	赵吉坤,宋武斌,李晶晶.基于的水稻秸秆建模及力学性能分析[J].土壤通报,2020,51(5):1086-1093.
	ZHAO J K, SONG W B, LI J J. Modeling and mechanical properties analysis of rice straw based on EDEM [J]. Soil Bull., 2020, 51(5):1086-1093.
14	康宏彬,刘铭,王雷,等.基于的马铃薯收获机分离输送装置仿真分析[J].农机化研究,2022,44(5):1-8.
	KANG H B, LIU M, WANG L, et al.. Simulation analysis of separation and conveying device of potato harvester based on EDEM [J]. Res. Agric. Mech., 2022, 44(5):1-8.
15	刘瑞,李衍军,刘忠军,等.包衣玉米种子离散元参数分析与标定[J].农业机械学报,2021,52(S1):1-8.
	LIU R, LI Y J, LIU Z J, et al.. Discrete element parameter analysis and calibration of coated corn seeds [J]. J. Agric. Mach., 2021, 52(S1):1-8.
16	马彦华,宋春东,宣传忠,等.苜蓿秸秆压缩仿真离散元模型参数标定[J].农业工程学报,2020,36(11):22-30.
	MA Y H, SONG C D, XUAN C Z, et al.. Parameter calibration of discrete element model for alfalfa straw compression simulation [J]. Trans. Chin. Soc. Agric. Eng., 2020, 36(11):22-30.
17	张涛,刘飞,赵满全,等.玉米秸秆接触物理参数测定与离散元仿真标定[J].中国农业大学学报,2018,23(4):120-127.
	ZHANG T, LIU F, ZHAO M Q, et al.. Measurement of contact physical parameters of corn stalks and discrete element simulation calibration [J]. J. China Agric. Univ., 2018, 23(4):120-120.
18	赖庆辉,袁海阔,胡子武,等.三七种苗物料特性研究及离散元法参数标定[J].扬州大学学报(农业与生命科学版), 2018,39(2):74-79.
	LAI Q H, YUAN H K, HU Z W, et al.. Study on material properties of Panax notoginseng seedlings and parameter calibration by discrete element method [J]. J. Yangzhou Univ. (Agric. Life Sci.), 2018, 39(2):74-79.
19	中华人民共和国国家卫生和计划生育委员会. 食品安全国家标准食品中水分的测定: . [S].北京:中国标准出版社, 2016.
20	于克强.转轮式全混合日粮混合机混合机理分析及试验研究[D].吉林:东北农业大学, 2015.
	YU K Q. Mixing mechanism analysis and experimental research of rotary total mixed ration mixer [D]. Jilin: Northeast Agricultural University, 2015.
21	ZHOU L, YU J Q, LIANG L S, et al.. DEM parameter calibration of maize seeds and the effect of rolling friction [J]. Processes, 2021, 9(6):914-914.
22	胡广锐,卜令昕,张恩宇,等.苹果两向异性力学特性和跌落试验研究[J].农机化研究, 2021, 43(4):154-160.
	HU G R, BU L X, ZHANG E Y, et al.. Study on anisotropic mechanical properties and drop test of apples [J]. Res. Agric. Mech., 2021, 43(4):154-160.
23	WANG X W, MA H Z, LI B, et al.. Review on the research of contact parameters calibration of particle system [J]. J. Mech. Sci. Technol., 2022, 36:1-1.
24	贾富国,韩燕龙,刘扬,等.稻谷颗粒物料堆积角模拟预测方法[J].农业工程学报, 2014, 30(11):254-260.
	JIA F G, HAN Y L, LIU Y, et al.. Simulation and prediction method of accumulation angle of rice granular materials [J]. Trans. Chin. Soc. Agric. Eng., 2014, 30(11):254-260.
25	张锋伟,宋学锋,张雪坤等.玉米秸秆揉丝破碎过程力学特性仿真与试验[J].农业工程学报,2019,35(9):58-65.
	ZHANG F W, SONG X F, ZHANG X K, et al. Simulation and test of mechanical properties in the process of corn straw kneading and crushing [J]. Trans. Chin. Soc. Agric. Eng., 2019, 35(9):58-65

对象 Object	静摩擦系数 Static friction coefficient			滚动摩擦系数 Rolling friction coefficient			碰撞恢复系数 Collision restitution coefficient
对象 Object	穰 Straw core	叶 Leaf	皮 Skin	穰 Straw core	叶 Leaf	皮 Skin	穰 Straw core	叶 Leaf	皮 Skin
穰 Straw core	2.25~3.08			0.23~0.36			0.36~0.42
叶 Leaf	0.75~1.04	0.31~0.40		0.27~0.36	0.29~0.38		0.29~0.38	0.29~0.40
皮 Skin	0.78~1.04	0.31~0.40	0.62~0.75	0.31~0.38	0.29~0.42	0.27~0.32	0.36~0.43	0.36~0.41	0.33~0.43
45钢 45 steel	1.00~1.07	0.36~0.40	0.31~0.34	0.31~0.42	0.18~0.25	0.25~0.29	0.37~0.46	0.36~0.43	0.40~0.45

对象 Object	静摩擦系数 Static friction coefficient			滚动摩擦系数 Rolling friction coefficient			碰撞恢复系数 Collision restitution coefficient
对象 Object	穰 Straw core	叶 Leaf	皮 Skin	穰 Straw core	叶 Leaf	皮 Skin	穰 Straw core	叶 Leaf	皮 Skin
穰 Straw core	2.25~3.08			0.23~0.36			0.36~0.42
叶 Leaf	0.75~1.04	0.31~0.40		0.27~0.36	0.29~0.38		0.29~0.38	0.29~0.40
皮 Skin	0.78~1.04	0.31~0.40	0.62~0.75	0.31~0.38	0.29~0.42	0.27~0.32	0.36~0.43	0.36~0.41	0.33~0.43
45钢 45 steel	1.00~1.07	0.36~0.40	0.31~0.34	0.31~0.42	0.18~0.25	0.25~0.29	0.37~0.46	0.36~0.43	0.40~0.45

代号 Code	仿真参数 Simulation parameter	水平 Level		代号 Code	仿真参数 Simulation parameter	水平 Level
代号 Code	仿真参数 Simulation parameter	-1	+1	代号 Code	仿真参数 Simulation parameter	-1	+1
A	穰-穰静摩擦系数 Straw corestraw core static friction coefficient	2.14	3.27	G	皮-皮静摩擦系数 Skinskin static friction coefficient	0.58	0.75
B	穰-穰滚动摩擦系数 Straw corestraw core rolling friction coefficient	0.23	0.47	H	皮-叶静摩擦系数 Skinlesf static friction coefficient	0.31	0.45
C	穰-皮静摩擦系数 Straw coreskin static friction coefficient	0.73	1.07	I	皮-叶滚动摩擦系数 Skinlesf rolling friction coefficient	0.29	0.45
D	穰-皮碰撞恢复系数 Straw coreskin collision restitution coefficient	0.33	0.47	J	叶-叶静摩擦系数 Lesflesf static friction coefficient	0.31	0.45
E	穰-叶静摩擦系数 Straw corelesf static friction coefficient	0.75	1.11	K	叶-叶碰撞恢复系数 Lesflesf collision restitution coefficient	0.27	0.40
F	穰-叶滚动摩擦系数 Straw corelesf rolling friction coefficient	0.25	0.42

代号 Code	仿真参数 Simulation parameter	水平 Level		代号 Code	仿真参数 Simulation parameter	水平 Level
代号 Code	仿真参数 Simulation parameter	-1	+1	代号 Code	仿真参数 Simulation parameter	-1	+1
A	穰-穰静摩擦系数 Straw corestraw core static friction coefficient	2.14	3.27	G	皮-皮静摩擦系数 Skinskin static friction coefficient	0.58	0.75
B	穰-穰滚动摩擦系数 Straw corestraw core rolling friction coefficient	0.23	0.47	H	皮-叶静摩擦系数 Skinlesf static friction coefficient	0.31	0.45
C	穰-皮静摩擦系数 Straw coreskin static friction coefficient	0.73	1.07	I	皮-叶滚动摩擦系数 Skinlesf rolling friction coefficient	0.29	0.45
D	穰-皮碰撞恢复系数 Straw coreskin collision restitution coefficient	0.33	0.47	J	叶-叶静摩擦系数 Lesflesf static friction coefficient	0.31	0.45
E	穰-叶静摩擦系数 Straw corelesf static friction coefficient	0.75	1.11	K	叶-叶碰撞恢复系数 Lesflesf collision restitution coefficient	0.27	0.40
F	穰-叶滚动摩擦系数 Straw corelesf rolling friction coefficient	0.25	0.42

参数 Parameter	调整后的偏差平方和 Adjusted sum of squared deviations	F值 F value	P值 P value
A	0.92	8.79	0.207 1
B	14.39	137.64	0.054 1
C	1.17	11.15	0.185 2
D	3.54	33.89	0.108 3
E	5.55	53.08	0.086 8
F	28.21	269.90^*	0.038 7
G	0.95	9.11	0.203 7
H	4.86	47.56	0.162 7
I	19.25	184.18^*	0.046 8
J	50.43	482.43^*	0.029 0
K	7.94	75.94	0.072 7