浙江农业学报 ›› 2026, Vol. 38 ›› Issue (7): 1471-1480.DOI: 10.3969/j.issn.1004-1524.20250269
收稿日期:2025-04-01
出版日期:2026-07-25
发布日期:2026-08-20
作者简介:史荣旭,研究方向为农业机械装备。E-mail:shirongxu928@163.com
通讯作者:
*宋学锋,E-mail: songxf@gsau.edu.cn
基金资助:
SHI Rongxu(
), SONG Xuefeng*(
), ZHANG Fengwei, DAI Fei
Received:2025-04-01
Published:2026-07-25
Online:2026-08-20
摘要:
为了提高玉米秸秆加工环节中离散元仿真研究所用参数的准确度,本研究以甘肃地区种植的大京九26玉米秸秆为研究对象,采用物理试验和仿真分析相结合的方法,进行玉米秸秆双层柔性模型的参数标定研究。首先,应用Plackett-Burman试验对玉米秸秆双层柔性模型的初始参数进行筛选,方差分析结果表明,玉米秸秆内穰弹性比、表皮弹性比和表皮塑性比对玉米秸秆轴向压缩破裂临界载荷有显著(p<0.05)影响。其次,以破裂临界载荷为评价指标,采用最陡爬坡和Box-Behnken试验建立了破裂临界载荷与上述3项参数的二次多项式回归模型,以物理试验得到的实际破裂临界载荷2 235.76 N为目标值,对各参数进行寻优,得到最优参数组合为内穰弹性比0.859、表皮弹性比0.848和表皮塑性比0.869。最后,在该参数组合下构建玉米秸秆双层柔性模型,进行玉米秸秆轴向压缩、径向压缩和剪切试验对比,结果表明,轴向压缩仿真试验与物理试验的破裂临界载荷相对误差为0.55%,径向压缩仿真试验与物理试验的破裂临界载荷相对误差为6.87%,剪切仿真试验与物理试验的临界剪切力相对误差为8.50%,且秸秆受力变化趋势基本一致,说明最优参数组合具有可行性与准确性。本研究标定的玉米秸秆双层柔性模型可为秸秆加工装备的设计优化提供参考。
中图分类号:
史荣旭, 宋学锋, 张锋伟, 戴飞. 基于离散元法的玉米秸秆双层柔性模型参数标定与试验[J]. 浙江农业学报, 2026, 38(7): 1471-1480.
SHI Rongxu, SONG Xuefeng, ZHANG Fengwei, DAI Fei. Parameter calibration and experiment of maize straw double-layer flexible model based on discrete element method[J]. Acta Agriculturae Zhejiangensis, 2026, 38(7): 1471-1480.
图3 纤维状柔性颗粒(左)与壳状柔性颗粒(右)的内力图 Fn为关节的法向力;Ft为关节的切向力;MT为关节的扭矩;MB1和MB2为关节的弯曲力矩。
Fig.3 Schematic of the internal force on fibrous flexible particles (left) and shell-shaped flexible particles (right) Fn is the normal force for joints; Ft is the tangential force for joints; MT is the torque for joint; MB1 and MB2 are the bending moments for the joint.
| 序号 No. | 各因素的设定水平Setting level for different factors | 破裂临界载荷/N Critical fracture load/N | ||||||
|---|---|---|---|---|---|---|---|---|
| A | B | C | D | E | F | G | ||
| 1 | -1 | -1 | -1 | 1 | -1 | 1 | 1 | 2 030.14 |
| 2 | 1 | -1 | 1 | 1 | -1 | 1 | 1 | 2 164.87 |
| 3 | -1 | -1 | 1 | -1 | 1 | 1 | -1 | 2 088.23 |
| 4 | 1 | -1 | 1 | 1 | 1 | -1 | -1 | 2 418.19 |
| 5 | -1 | 1 | 1 | 1 | -1 | -1 | -1 | 2 051.58 |
| 6 | -1 | 1 | 1 | -1 | 1 | 1 | 1 | 1 943.17 |
| 7 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | 1 902.05 |
| 8 | 1 | -1 | -1 | -1 | 1 | -1 | 1 | 2 173.64 |
| 9 | 1 | 1 | -1 | -1 | -1 | 1 | -1 | 1 852.11 |
| 10 | -1 | 1 | -1 | 1 | 1 | -1 | 1 | 2 340.54 |
| 11 | 1 | 1 | 1 | -1 | -1 | -1 | 1 | 1 855.68 |
| 12 | 1 | 1 | -1 | 1 | 1 | 1 | -1 | 2 335.36 |
表1 Plackett-Burman试验设计与结果
Table 1 Design and result of Plackett-Burman test
| 序号 No. | 各因素的设定水平Setting level for different factors | 破裂临界载荷/N Critical fracture load/N | ||||||
|---|---|---|---|---|---|---|---|---|
| A | B | C | D | E | F | G | ||
| 1 | -1 | -1 | -1 | 1 | -1 | 1 | 1 | 2 030.14 |
| 2 | 1 | -1 | 1 | 1 | -1 | 1 | 1 | 2 164.87 |
| 3 | -1 | -1 | 1 | -1 | 1 | 1 | -1 | 2 088.23 |
| 4 | 1 | -1 | 1 | 1 | 1 | -1 | -1 | 2 418.19 |
| 5 | -1 | 1 | 1 | 1 | -1 | -1 | -1 | 2 051.58 |
| 6 | -1 | 1 | 1 | -1 | 1 | 1 | 1 | 1 943.17 |
| 7 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | 1 902.05 |
| 8 | 1 | -1 | -1 | -1 | 1 | -1 | 1 | 2 173.64 |
| 9 | 1 | 1 | -1 | -1 | -1 | 1 | -1 | 1 852.11 |
| 10 | -1 | 1 | -1 | 1 | 1 | -1 | 1 | 2 340.54 |
| 11 | 1 | 1 | 1 | -1 | -1 | -1 | 1 | 1 855.68 |
| 12 | 1 | 1 | -1 | 1 | 1 | 1 | -1 | 2 335.36 |
| 序号 No. | 各因素的设定值 Setting value of each factor | 破裂临界载荷/N Critical fracture load/N | 相对误差/% Relative error/% | ||
|---|---|---|---|---|---|
| A | D | E | |||
| 1 | 0.70 | 0.70 | 0.70 | 1 627.36 | 27.21 |
| 2 | 0.75 | 0.75 | 0.75 | 1 810.82 | 19.01 |
| 3 | 0.80 | 0.80 | 0.80 | 1 934.51 | 13.47 |
| 4 | 0.85 | 0.85 | 0.85 | 2 209.24 | 1.37 |
| 5 | 0.90 | 0.90 | 0.90 | 2 351.53 | 5.18 |
| 6 | 0.95 | 0.95 | 0.95 | 2 446.78 | 9.44 |
表2 最陡爬坡试验设计与结果
Table 2 Design and results of the steepest ascent experiment
| 序号 No. | 各因素的设定值 Setting value of each factor | 破裂临界载荷/N Critical fracture load/N | 相对误差/% Relative error/% | ||
|---|---|---|---|---|---|
| A | D | E | |||
| 1 | 0.70 | 0.70 | 0.70 | 1 627.36 | 27.21 |
| 2 | 0.75 | 0.75 | 0.75 | 1 810.82 | 19.01 |
| 3 | 0.80 | 0.80 | 0.80 | 1 934.51 | 13.47 |
| 4 | 0.85 | 0.85 | 0.85 | 2 209.24 | 1.37 |
| 5 | 0.90 | 0.90 | 0.90 | 2 351.53 | 5.18 |
| 6 | 0.95 | 0.95 | 0.95 | 2 446.78 | 9.44 |
| 序号 No. | 各因素的设定水平 Setting level of each factor | 破裂临界载荷/N Critical fracture load/N | |||
|---|---|---|---|---|---|
| A | D | E | |||
| 1 | 1 | 0 | -1 | 2 140.77 | |
| 2 | 0 | 0 | 0 | 2 205.37 | |
| 3 | 1 | 1 | 0 | 2 300.97 | |
| 4 | 0 | 1 | 1 | 2 320.36 | |
| 5 | 0 | 0 | 0 | 2 209.85 | |
| 6 | 0 | -1 | -1 | 1 951.44 | |
| 7 | -1 | 0 | -1 | 2 109.01 | |
| 8 | -1 | 1 | 0 | 2 211.26 | |
| 9 | 0 | 0 | 0 | 2 216.11 | |
| 10 | 0 | 0 | 0 | 2 224.95 | |
| 11 | -1 | 0 | 1 | 2 179.65 | |
| 12 | 0 | -1 | 1 | 2 041.99 | |
| 13 | 1 | -1 | 0 | 2 029.04 | |
| 14 | 0 | 0 | 0 | 2 242.57 | |
| 15 | 0 | 1 | -1 | 2 232.54 | |
| 16 | 1 | 0 | 1 | 2 235.85 | |
| 17 | -1 | -1 | 0 | 1 988.69 | |
表3 Box-Behnken试验设计与结果
Table 3 Design and results of Box-Behnken experiment
| 序号 No. | 各因素的设定水平 Setting level of each factor | 破裂临界载荷/N Critical fracture load/N | |||
|---|---|---|---|---|---|
| A | D | E | |||
| 1 | 1 | 0 | -1 | 2 140.77 | |
| 2 | 0 | 0 | 0 | 2 205.37 | |
| 3 | 1 | 1 | 0 | 2 300.97 | |
| 4 | 0 | 1 | 1 | 2 320.36 | |
| 5 | 0 | 0 | 0 | 2 209.85 | |
| 6 | 0 | -1 | -1 | 1 951.44 | |
| 7 | -1 | 0 | -1 | 2 109.01 | |
| 8 | -1 | 1 | 0 | 2 211.26 | |
| 9 | 0 | 0 | 0 | 2 216.11 | |
| 10 | 0 | 0 | 0 | 2 224.95 | |
| 11 | -1 | 0 | 1 | 2 179.65 | |
| 12 | 0 | -1 | 1 | 2 041.99 | |
| 13 | 1 | -1 | 0 | 2 029.04 | |
| 14 | 0 | 0 | 0 | 2 242.57 | |
| 15 | 0 | 1 | -1 | 2 232.54 | |
| 16 | 1 | 0 | 1 | 2 235.85 | |
| 17 | -1 | -1 | 0 | 1 988.69 | |
图6 物理试验与仿真试验的结果 FSac和FEac分别为轴向压缩试验中仿真与物理试验的破裂临界载荷;FSrc和FErc分别为径向压缩试验中仿真与物理试验的破裂临界载荷;FSs和FEs分别为剪切试验中仿真与物理试验的最大剪切力。
Fig.6 Results of physical and simulation test FSac and FEac represent the critical fracture loads obtained from simulation and physical test in axial compression test, respectively; FSrc and FErc represent the critical fracture loads obtained from simulation and physical test in radial compression test, respectively; FSs and FEs represent the maximum shear forces obtained from simulation and physical test in shear test, respectively.
| [1] | 窦婉毓, 颜双双, 孙仁华, 等. 基于21世纪以来我国作物产量、人口与政策分析未来秸秆还田的必然性[J]. 中国农业大学学报, 2025, 30(12): 291-305. |
| Dou W Y, Yan S S, Sun R H, et al. The inevitability of returning straw to the field in the future based on crop yield, population, and policy in China since the 21st century[J]. Journal of China Agricultural University, 2025, 30(12): 291-305. | |
| [2] | 赵伟, 黄红英, 蒋建雄, 等. 全秸秆纤维育苗穴盘的制备及其性能分析[J]. 农业工程学报, 2025, 41(2): 233-241. |
| Zhao W, Huang H Y, Jiang J X, et al. Preparation and performance analysis of whole straw fiber seedling trays[J]. Transactions of the Chinese Society of Agricultural Engineering, 2025, 41(2): 233-241. | |
| [3] | Li H T, Miao Z, Zhang B. Review and prospect of comprehensive straw utilization and government policy in China[J]. Chinese Journal of Population, Resources and Environment, 2022, 20(4): 402-406. |
| [4] | 牟孝栋. 基于碟盘辊的全株青贮玉米破碎机理研究与试验[D]. 淄博: 山东理工大学, 2022. |
| Mu X D. Study and experiment on crushing mechanism of whole silage maize based on dish roller[D]. Zibo: Shandong University of Technology, 2022. | |
| [5] | 刘城宇, 贺正楚, 卢小龙. 可持续发展目标下农作物秸秆收集运输的碳减排优化分析[J]. 农业工程学报, 2022, 38(10): 239-248. |
| Liu C Y, He Z C, Lu X L. Optimization analysis of carbon emission reduction from crop straw collection and transportation under the sustainable development goals[J]. Transactions of the Chinese Society of Agricultural Engineering, 2022, 38(10): 239-248. | |
| [6] | 李胜男, 纪雄辉, 邓凯, 等. 区域秸秆资源分布及全量化利用潜力分析[J]. 农业工程学报, 2020, 36(12): 221-228. |
| Li S N, Ji X H, Deng K, et al. Analysis of regional distribution patterns and full utilization potential of crop straw resources[J]. Transactions of the Chinese Society of Agricultural Engineering, 2020, 36(12): 221-228. | |
| [7] | 霍丽丽, 赵立欣, 孟海波, 等. 中国农作物秸秆综合利用潜力研究[J]. 农业工程学报, 2019, 35(13): 218-224. |
| Huo L L, Zhao L X, Meng H B, et al. Study on straw multi-use potential in China[J]. Transactions of the Chinese Society of Agricultural Engineering, 2019, 35(13): 218-224. | |
| [8] | 王韦韦, 蔡丹艳, 谢进杰, 等. 玉米秸秆粉料致密成型离散元模型参数标定[J]. 农业机械学报, 2021, 52(3): 127-134. |
| Wang W W, Cai D Y, Xie J J, et al. Parameters calibration of discrete element model for corn stalk powder compression simulation[J]. Transactions of the Chinese Society for Agricultural Machinery, 2021, 52(3): 127-134. | |
| [9] | 马彦华, 宋春东, 宣传忠, 等. 苜蓿秸秆压缩仿真离散元模型参数标定[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]. Transactions of the Chinese Society of Agricultural Engineering, 2020, 36(11): 22-30. | |
| [10] | 樊琦. 玉米秸秆破碎力学性能试验研究[D]. 呼和浩特: 内蒙古农业大学, 2019. |
| Fan Q. Study on mechanical properties of corn stalk crushing[D]. Hohhot: Inner Mongolia Agricultural University, 2019. | |
| [11] | 赵玉清, 赵勇, 杨志, 等. 青贮玉米秸秆径向和轴向剪切特性试验研究[J]. 饲料研究, 2021, 44(17): 93-96. |
| Zhao Y Q, Zhao Y, Yang Z, et al. Experimental study on the radial and axial shear characteristics of silage corn stalk[J]. Feed Research, 2021, 44(17): 93-96. | |
| [12] | 陈争光, 王德福, 李利桥, 等. 玉米秸秆皮拉伸和剪切特性试验[J]. 农业工程学报, 2012, 28(21): 59-65. |
| Chen Z G, Wang D F, Li L Q, et al. Experiment on tensile and shearing characteristics of rind of corn stalk[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012, 28(21): 59-65. | |
| [13] | 郭思佳, 刘新柱, 刘春山, 等. 玉米秸秆压缩特性试验研究[J]. 农机化研究, 2025, 47(2): 153-158. |
| Guo S J, Liu X Z, Liu C S, et al. Experiment study on compression characteristics of corn straw[J]. Journal of Agricultural Mechanization Research, 2025, 47(2): 153-158. | |
| [14] | Xia Y D, Klinger J, Bhattacharjee T, et al. The elastoplastic flexural behaviour of corn stalks[J]. Biosystems Engineering, 2022, 216: 218-228. |
| [15] | 张李娴. 玉米秸秆力学特性的离散元建模方法研究[D]. 杨凌: 西北农林科技大学, 2017. |
| Zhang L X. Research on the discrete element modeling method of corn stalk’s mechanical characteristics[D]. Yangling: Northwest A & F University, 2017. | |
| [16] | Zhang S L, Zhao H B, Wang X Z, et al. Discrete element modeling and shear properties of the maize stubble-soil complex[J]. Computers and Electronics in Agriculture, 2023, 204: 107519. |
| [17] | Liu L, Wang X L, Zhang X C, et al. Determination and verification of parameters for the discrete element modelling of single disc covering of flexible straw with soil[J]. Biosystems Engineering, 2023, 233: 151-167. |
| [18] | Leblicq T, Smeets B, Ramon H, et al. A discrete element approach for modelling the compression of crop stems[J]. Computers and Electronics in Agriculture, 2016, 123: 80-88. |
| [19] | 张锋伟, 宋学锋, 张雪坤, 等. 玉米秸秆揉丝破碎过程力学特性仿真与试验[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]. Transactions of the Chinese Society of Agricultural Engineering, 2019, 35(9): 58-65. | |
| [20] | 刘禹辰, 张锋伟, 宋学锋, 等. 基于离散元法玉米秸秆双层粘结模型力学特性研究[J]. 东北农业大学学报, 2022, 53(1): 45-54. |
| Liu Y C, Zhang F W, Song X F, et al. Study on mechanical properties for corn straw of double-layer bonding model based on discrete element method[J]. Journal of Northeast Agricultural University, 2022, 53(1): 45-54. | |
| [21] | 史瑞杰, 戴飞, 赵武云, 等. 胡麻茎秆离散元柔性模型建立与接触参数试验验证[J]. 农业机械学报, 2022, 53(10): 146-155. |
| Shi R J, Dai F, Zhao W Y, et al. Establishment of discrete element flexible model and verification of contact parameters of flax stem[J]. Transactions of the Chinese Society for Agricultural Machinery, 2022, 53(10): 146-155. | |
| [22] | Shi L R, Yang X P, Zhao W Y, et al. Investigation of interaction effect between static and rolling friction of corn kernels on repose formation by DEM[J]. International Journal of Agricultural and Biological Engineering, 2021, 14(5): 238-246. |
| [23] | Dai F, Song X F, Zhao W Y, et al. Motion simulation and test on threshed grains in tapered threshing and transmission device for plot wheat breeding based on CFD-DEM[J]. International Journal of Agricultural and Biological Engineering, 2019, 12(1): 66-73. |
| [24] | Dai F, Song X F, Guo W J, et al. Simulation and test on separating cleaning process of flax threshing material based on gas-solid coupling theory[J]. International Journal of Agricultural and Biological Engineering, 2020, 13(1): 73-81. |
| [25] | Mao H P, Wang Q R, Li Q L. Modelling and simulation of the straw-grain separation process based on a discrete element model with flexible hollow cylindrical bonds[J]. Computers and Electronics in Agriculture, 2020, 170: 105229. |
| [26] | 夏俊芳, 张朋, 苑宏文, 等. 水稻秸秆离散元柔性模型参数标定与试验验证[J]. 农业机械学报, 2024, 55(9): 174-184. |
| Xia J F, Zhang P, Yuan H W, et al. Calibration and verification of flexible rice straw model by discrete element method[J]. Transactions of the Chinese Society for Agricultural Machinery, 2024, 55(9): 174-184. | |
| [27] | Freireich B, Litster J, Wassgren C. Using the discrete element method to predict collision-scale behavior: a sensitivity analysis[J]. Chemical Engineering Science, 2009, 64(15): 3407-3416. |
| [28] | El-Emam M A, Shi W D, Zhou L. CFD-DEM simulation and optimization of gas-cyclone performance with realistic macroscopic particulate matter[J]. Advanced Powder Technology, 2019, 30(11): 2686-2702. |
| [1] | 任宁, 俞国红, 郑航, 陈志东. 基于离散元法的茶园土壤参数标定[J]. 浙江农业学报, 2025, 37(6): 1353-1359. |
| [2] | 杨锦皓, 江洁, 刘行, 李福强. 基于离散元法的水稻种子参数标定与旱作排种仿真[J]. 浙江农业学报, 2025, 37(11): 2376-2386. |
| [3] | 熊厚全, 江洁. 黄色赤红壤与开沟部件互作的离散元参数标定及试验[J]. 浙江农业学报, 2024, 36(6): 1400-1412. |
| [4] | 孔德, 叶自然, 谭向峰, 代梦迪, 赵先亮, 孔德栋. 生菜种子物理接触参数测定与离散元仿真标定[J]. 浙江农业学报, 2024, 36(4): 932-942. |
| [5] | 彭才望, 周婷, 孙松林, 谢烨林, 魏源. 基于堆积试验的黑水虻离散元仿真参数标定与分析[J]. 浙江农业学报, 2022, 34(4): 814-823. |
| [6] | 王锋, 张锋伟, 戴飞, 张陆海, 赵伟, 杨小平. 双层平筛式半夏收获机设计与试验[J]. 浙江农业学报, 2021, 33(10): 1946-1955. |
| [7] | 李为宁, 柏宣丙, 李兵. 滚筒式茶叶提香机结构参数优化[J]. 浙江农业学报, 2020, 32(2): 348-358. |
| [8] | 邱秀文, 周桂香, 王慧娟, 杨丽丽. 纤维素分解菌的分离筛选及其对秸秆的分解特性[J]. 浙江农业学报, 2017, 29(4): 637-643. |
| [9] | 邵泱峰1,梅洪飞1,潘忠潮1,刘欢2,王超琦2. 玉米秸秆还田对土壤有机碳、微生物功能多样性及甘蓝产量的影响[J]. 浙江农业学报, 2016, 28(5): 838-. |
| 阅读次数 | ||||||
|
全文 |
|
|||||
|
摘要 |
|
|||||