基于EDEM的水稻残茬秸秆离散元仿真参数标定

doi:10.13304/j.nykjdb.2022.1047

摘要/Abstract

摘要：

为了提高离散元法模拟水稻秸秆在切碎还田过程中的准确性，以水稻收获后田间的残茬秸秆为研究对象，通过斜面试验测定秸秆与钢、秸秆与秸秆之间的接触参数，基于离散元仿真模型，通过堆积角试验和剪切试验，以秸秆模型颗粒之间的粘结半径、切向临界应力、法向临界应力，秸秆之间的恢复系数、静摩擦系数、动摩擦系数为试验因素，以秸秆的抗剪力与径向堆积角为评价指标，利用Box-Behnken 试验进行多目标优化仿真标定。结果表明，秸秆颗粒之间的粘结半径为1.06 mm，秸秆法向临界应力为4.77×10¹⁰ Pa、切向临界应力为4.67×10⁶ Pa；秸秆与秸秆之间的碰撞恢复系数为0.21、静摩擦系数为0.19、动摩擦系数为0.09。最优参数组合条件下的仿真试验与物理试验对比验证表明，仿真试验与物理试验中秸秆的抗剪力、径向堆积角的相对误差分别为1.7%、2.8%。研究结果可为秸秆与秸秆之间、秸秆与农机具之间的离散元仿真分析提供参考。

关键词: 水稻, 秸秆, 离散元法, 试验

Abstract:

In order to improve the accuracy of the parameters during the process of rice straw shredding and returning to the field in discrete element simulation， the stubble straw in the field after rice harvest were taken as the research object， and the contact parameters between straw and steel， straw and straw were measured through the slope test.Then based on simulation model and through the stacking angle and shear test， the bonding radius between straw model particles， tangential critical stress， normal critical stress， coefficient of restitution， static friction coefficient and dynamic friction coefficient were taken as the test factors， taking the shear resistance and radial stacking angle of straw as evaluation indicators， a multi-objective optimization simulation calibration experiment were conducted using the Box-Behnken experimental scheme.The results showed that the bonding radius between straw particles was 1.06 mm， and the normal critical stress of straw was 4.77×10¹⁰ Pa， tangential critical stress was 4.67×10⁶ Pa， the collision restitution coefficient between straw and straw was 0.21， the static Friction coefficient was 0.19， and the rolling friction coefficient was 0.09. The results of the simulation test results under the optimal parameter combination conditions and the physical test showed that the relative errors of the shear resistance and radial stacking angle of straw in the simulation test and physical test were 1.7% and 2.8%， respectively. This study could provide reference for the discrete element simulation analysis between straw and straw， as well as between straw and agricultural machinery.

Key words: rice, straw, discrete element method, test

中图分类号:

S126

马紫涛, 赵智豪, 全伟, 石方刚, 高晨, 吴明亮. 基于EDEM的水稻残茬秸秆离散元仿真参数标定[J]. 中国农业科技导报, 2023, 25(11): 103-113.

Zitao MA, Zhihao ZHAO, Wei QUAN, Fanggang SHI, Chen GAO, Mingliang WU. Calibration of Discrete Element Parameter of Rice Stubble Straw Based on EDEM[J]. Journal of Agricultural Science and Technology, 2023, 25(11): 103-113.

图/表 14

参考文献 30

1	刘淑军, 李冬初, 黄晶, 等. 1988—2018年中国水稻秸秆资源时空分布特征及还田替代化肥潜力[J]. 农业工程学报, 2021,37(11):151-161.
	LIU S J, LI D C, HUANG J,et al.. Temporal and spatial distribution characteristics of rice stalk resources and its potential of synthetic fertilizers substitution returning to farmland in China from 1988 to 2018 [J]. Trans. Chin. Soc. Agric. Eng., 2021,37(11):151-161.
2	侯杰, 谢方平, 王修善, 等. 水稻茎秆接触物理参数测定与离散元仿真标定[J]. 江西农业大学学报, 2022,44(3):747-758.
	HOU J, XIE F P, WANG X S, et al.. Measurement of contact physical parameters of flexible rice straw and discrete element simulation calibration [J]. Acta Agric. Univ. Jiangxiensis, 2022,44(3):747-758.
3	邱进, 吴明亮, 方友祥, 等. 水稻秸秆利用研究现状与发展趋势[J]. 当代农机, 2015(4):72-75.
	QIU J, WU M L, FANG Y X,et al.. Research status and development trend of rice straw scales [J].Contemporary Farm Mach., 2015(4):72-75.
4	逯克鑫, 屈佳伟, 于晓芳, 等. 沙壤土条件下秸秆还田减施氮肥对不同玉米品种氮效率的影响[J]. 内蒙古农业大学学报(自然科学版), 2022, 43(5): 37-45.
	LU K X, QU J W, YU X F,et al.. Effect of straw returning and nitrogen fertilizer reduction on nitrogen efficiency of different maize varieties in sandy ioam soils [J]. J. Inner Mongolia Agric.Univ. (Nat. Sci. ), 2022, 43(5): 37-45.
5	张博睿, 蔡影, 刘文雯, 等. 水稻秸秆还田对小麦季农田径流氮磷及COD的影响[J]. 安徽科技学院学报, 2022,36(1):18-24.
	ZHANG B R, CAI Y, LIU W W,et al.. Rice straw return measures on nitrogen, phosphorus and COD of wheat season in farm runoff [J]. J. Anhui Sci. Technol. Univ., 2022,36(1):18-24.
6	崔思远, 朱新开, 张莀茜, 等. 水稻秸秆还田年限对稻麦轮作田土壤碳氮固存的影响[J]. 农业工程学报, 2019,35(7):115-121.
	CUI S Y, ZHU X K, ZHANG C Q,et al.. Effects of years of rice straw retention on soil carbon and nitrogen sequestration in rice-wheat system [J].Trans. Chin. Soc. Agric. Eng., 2019,35(7):115-121.
7	曹阳, 陈新兵, 沙之敏, 等. 水稻秸秆还田量对土壤重金属及小麦产量的影响[J]. 上海交通大学学报(农业科学版), 2019,37(4):6-11.
	CAO Y, CHEN X B, SHA Z M,et al.. Effect of rice straw returning on accumulation of heavy metals in soil and yield of wheat [J]. J. Shanghai Jiaotong Univ. (Agric. Sci.), 2019,37(4):6-11.
8	曹若梅. 水稻秸秆还田对土壤含水量的影响[J]. 安徽农业科学, 2019,47(18):85-87.
	CAO R M. Effects of returning rice straw to field on soil moisture content [J]. J. Anhui Agric. Sci., 2019,47(18):85-87.
9	邱进, 吴明亮, 官春云, 等. 动定刀同轴水稻秸秆切碎还田装置结构设计与试验[J]. 农业工程学报, 2015,31(10):11-19.
	QIU J, WU M L, GUAN C Y, et al.. Design and experiment of chopping device with dynamic fixed knife coaxial for rice straw [J].Trans. Chin. Soc. Agric. Eng., 2015,31(10):11-19.
10	薛惠岚,薛少平,杨青,等. 秸秆粉碎覆盖与施肥播种联合作业的实现及机具设计[J]. 农业工程学报, 2003, 19(3):104-107.
	XUE H L, XUE S P, YANG Q, et al.. Implementation of combined work of straw crushed for mulching and seeding with fertilizer and design of the machine [J].Trans. Chin. Soc. Agric. Eng., 2003, 19(3):104-107.
11	贾洪雷, 赵佳乐, 姜鑫铭, 等. 行间免耕播种机防堵装置设计与试验[J]. 农业工程学报, 2013,29(18):16-25.
	JIA H L, ZHAO J L, JIANG X M, et al.. Design and experiment of anti-blocking mechanism for inter-row no-tillage seeder [J]. Trans. Chin. Soc. Agric. Eng., 2013, 29(18): 16-25.
12	陈晓光, 武涛, 张增学, 等. 基于EDEM甘蔗切割器入土切割仿真试验分析[J]. 农机化研究, 2021,43(11):7-13.
	CHEN X G, WU T, ZHANG Z X, et al.. Simulation test analysis of cutting into the soil based on EDEM sugarcane cutter [J]. J. Agric. Mech. Res., 2021,43(11):7-13.
13	周里群, 许欣, 关汗青, 等. 基于离散元法的沥青混凝土振动切削过程的数值模拟[J]. 机械工程学报, 2017,53(22):166-175.
	ZHOU L Q, XU X, GUAN H Q,et al.. Numerical simulation of vibration cutting for asphalt concrete paving based on discrete element method [J]. J. Mech. Eng.,2017,53(22):166-175.
14	赵淑红, 高连龙, 袁溢文, 等. 基于离散元法的深松作业玉米秸秆运动规律[J]. 农业工程学报, 2021,37(20):53-62.
	ZHAO S H, GAO L L, YUAN Y W, et al.. Maize straw motion law in subsoiling operation using discrete element method [J].Trans. Chin. Soc. Agric. Eng.,2021,37(20):53-62.
15	耿端阳, 孙延成, 牟孝栋, 等. 基于差速辊的青贮玉米籽粒破碎仿真试验及优化[J]. 吉林大学学报(工学版), 2022,52(3):693-702.
	GENG D Y, SUN Y C, MOU X D, et al.. Simulation test and optimization of grain breakage of silage maize based on differential roller [J]. J. Jilin Univ. (Eng. Technol.) 2022,52(3):693-702.
16	张锋伟, 宋学锋, 张雪坤, 等. 玉米秸秆揉丝破碎过程力学特性仿真与试验[J]. 农业工程学报, 2019,35(9):58-65.
	ZHANG F W, SONG X F, ZHANG X K, et al.. Simulation and experiment on mechanical characteristics of kneading and crushing process of corn straw [J]. Trans. Chin. Soc. Agric. Eng., 2019,35(9):58-65.
17	TEKESTE M Z, BALVANZ L R, HATFIELD J L, et al.. Discrete element modeling of cultivator sweep-to-soil interaction: worn and hardened edges effects on soil-tool forces and soil flow [J/OL]. J. Terramechanics, 2019, 82:1-11 [2022-10-30]. .
18	ABDIEL R L B, ULRICH H, THAI S D, et al.. Hypoplastic particle finite element model for cutting tool-soil interaction simulations: numerical analysis and experimental validation [J]. Underg. Space, 2018,3(1):61-71.
19	马彦华, 宋春东, 宣传忠, 等. 苜蓿秸秆压缩仿真离散元模型参数标定[J]. 农业工程学报, 2020,36(11):22-30.
	MA Y H, SONG C D, XUAN C Z, et al.. Parameters calibration of discrete element model for alfalfa straw compression simulation [J]. Trans. Chin. Soc. Agric. Eng., 2020,36(11):22-30.
20	张涛, 刘飞, 赵满全, 等. 玉米秸秆接触物理参数测定与离散元仿真标定[J]. 中国农业大学学报, 2018,23(4):120-127.
	ZHANG T, LIU F, ZHAO M Q, et al..Determination of corn stalk contact parameters and calibration of discrete element method simulation [J]. J. China Agric. Univ., 2018,23(4):120-127.
21	王万章, 刘婉茹, 袁玲合, 等. 基于EDEM的收获期小麦植株离散元参数标定[J]. 河南农业大学学报, 2021,55(1):64-72.
	WANG W Z, LIU W R, YUAN L H, et al..Calibration of discrete element parameters of wheat plants at harvest period based on EDEM [J]. J. Henan Agric. Univ., 2021,55(1):64-72.
22	ANNOUSSAMY M, RICHARD G, RECOUS S, et al.. Change in mechanical properties of wheat straw due to decomposition abd moisture [J]. Appl. Eng. Agric., 2000,16(6):657-664.
23	刘凡一. 清选装置中小麦颗粒和短茎秆离散元建模研究[D]. 杨凌:西北农林科技大学, 2018.
	LIU F Y. Discrete element modelling of the wheat particles and short straw in cleaning devices [D]. Yangling: Northwest A&F Universit, 2018.
24	JIA H L, DENG J Y, DENG Y L, et al.. Contact parameter analysis and calibration in discrete element simulation of rice straw [J]. Int. J. Agric. Biol. Eng, 2021,14(3):72-81.
25	MAK J, CHEN Y, SADEK M A. Determining parameters of a discrete element model for soil-tool interaction [J]. Soil Till. Res., 2012,118:117-122.
26	丁启朔, 任骏, BELAL Eisa Adam, 等. 湿粘水稻土深松过程离散元分析[J]. 农业机械学报, 2017,48(3):38-48.
	DING Q S, REN J, BELAL E A, et al.. DEM analysis of subsoiling process in wet clayey paddy soil [J] Trans. Chin. Soc. Agric. Eng, 2017,48(3):38-48.
27	沈成, 李显旺, 张彬, 等. 苎麻茎秆台架切割试验与分析[J]. 农业工程学报, 2016,32(1):68-76.
	SHEN C, LI X W, ZHANG B, et al..Bench experiment and analysis on ramie stalk cutting [J] Trans. Chin. Soc. Agric. Eng., 2016,32(1):68-76.
28	王洪波,樊志鹏,乌兰图雅,等.揉碎玉米秸秆螺旋输送仿真离散元模型参数标定[J].中国农业科技导报,2023,25(3):96-106.
	WANG H B, FAN Z P, Wulantuya, et al..Parameter calibration of discrete element model for simulation of crushed corn stalk screw conveying [J]. J. Agric. Sci. Technol. China,2023,25(3): 96-106.
29	陈声显. 玉米秸秆力学模型及压缩成型设备研究[D].长春:吉林大学, 2011.
	CHEN S X. Research on mechanical model and compression molding equipment of corn straw [D]. Changchun: Jilin University, 2011.
30	赵吉坤,宋武斌,李晶晶.基于EDEM的水稻秸秆建模及力学性能分析[J].土壤通报,2020,51(5):1086-1093.
	ZHAO J K, SONG W B, LI J J.Modeling and mechanical analysis of rice straw based on discrete element mechanical model [J]. Chin. J. Soil Sci., 2020,51(5):1086-1093.

材料 Material	参数 Parameter	数值 Value
水稻秸秆 Rice straw	泊松比 Poisson ratio	0.4
	剪切模量Modulus of shear /MPa	1
	密度 Density/（kg·m^-3）	126.4
钢 Steel	泊松比 Poisson ratio	0.3
	剪切模量 Modulus of shear/MPa	79 000
	密度 Density/（kg·m^-3）	7 800

材料 Material	参数 Parameter	数值 Value
水稻秸秆 Rice straw	泊松比 Poisson ratio	0.4
	剪切模量Modulus of shear /MPa	1
	密度 Density/（kg·m^-3）	126.4
钢 Steel	泊松比 Poisson ratio	0.3
	剪切模量 Modulus of shear/MPa	79 000
	密度 Density/（kg·m^-3）	7 800

序号 No.	法向临界应力 Normal critical stress/（10¹⁰ Pa）	切向临界应力 Tangential critical stress/（10⁶ Pa）	粘结半径 Stalk bonding radius/mm	抗剪力 Shear resistance/N	相对误差 Relative error/%
1	7	7	0.50	14.3	39.9
2	6	6	0.75	20.4	14.2
3	5	5	1.00	26.2	10.1
4	4	4	1.25	21.3	10.5
5	3	3	1.50	13.7	42.4

序号 No.	法向临界应力 Normal critical stress/（10¹⁰ Pa）	切向临界应力 Tangential critical stress/（10⁶ Pa）	粘结半径 Stalk bonding radius/mm	抗剪力 Shear resistance/N	相对误差 Relative error/%
1	7	7	0.50	14.3	39.9
2	6	6	0.75	20.4	14.2
3	5	5	1.00	26.2	10.1
4	4	4	1.25	21.3	10.5
5	3	3	1.50	13.7	42.4

试验因素Test factor	低水平Low level （-1）	中水平Middle level （0）	高水平High level （1）
x₁：粘结半径 Stalk bonding radius/mm	0.75	1.00	1.25
x₂：法向临界应力Normal critical stress/（10¹⁰ Pa）	4	5	6
x₃：切向临界应力Tangential critical stress/（10⁶ Pa）	4	5	6
A：碰撞恢复系数 Collision restitution coefficient	0.21	0.31	0.41
B：静摩擦系数 Static friction coefficient	0.13	0.23	0.33
C：动摩擦系数 Rolling friction coefficient	0.05	0.15	0.25