盛佳俊¹， 郑正鼎¹， 周子杨¹， 吴鹏民²

1.武汉科技大学 冶金装备及其控制教育部重点实验室， 机械传动与制造工程湖北省重点实验室， 湖北 武汉 430081；

2.五冶集团上海有限公司焦化工程建造标准研究所， 上海 201999

引用格式：

盛佳俊， 郑正鼎， 周子杨， 等. 基于离散元法的耐火泥颗粒仿真参数标定［J］. 中国粉体技术， 2026， 32（4）： 1-13.

Sheng Jiajun， Zheng Zhengding， Zhou Ziyang， et al. Calibration of simulation parameters for refractory mortar particles based on discrete element method［J］. China Powder Science and Technology， 2026， 32（4）： 1-13.

DOI：10.13732/j.issn.1008-5548.2026.04.011

收稿日期： 2025-09-22， 修回日期： 2026-04-02，上线日期： 2026-06-02。

基金项目：国家自然科学基金项目，编号：52505516。

第一作者：盛佳俊（1997—），男，硕士研究生，研究方向为机械设计。E-mail：shengjj97@163.com。

通信作者：郑正鼎（1995—），男，讲师，博士，硕士生导师，研究方向为先进制造装备设计与工艺开发。E-mail：zdzheng@wust.edu.cn。

摘要：【目的】结合耐火泥颗粒的微观参数及形态特征，以堆积角为响应值，建立离散元参数标定方法，为耐火泥泵送系统的CFD-DEM优化设计提供高保真输入参数。【方法】首先选择粒径分别为20、300、550 μm的耐火泥颗粒作为研究对象，建立类圆形、类块形和类锥形3种形态的耐火泥颗粒离散元模型，然后通过实际实验获取堆积角的数值，再通过Plackett-Burman试验设计筛选出待标定参数中的显著影响因素，接着利用最陡爬坡试验缩小显著因素的合理取值范围，最后通过Box-Behnken试验进行响应面分析得到二阶回归方程模型，确定反映耐火泥颗粒堆积特性的显著影响因素最优参数组合。【结果】耐火泥颗粒的堆积角平均值为41.87°；在7个待标定参数中，对堆积角的极其显著影响因素为颗粒表面能、颗粒-颗粒间滚动摩擦系数，显著影响因素为颗粒-颗粒间静摩擦系数，不显著影响因素为颗粒-颗粒间恢复系数、颗粒-管道间恢复系数、颗粒-管道间静摩擦系数以及颗粒-管道间滚动摩擦系数；通过Box-Behnken响应面分析法确定了最优参数组合，即颗粒表面能为0.004 J/m²、颗粒-颗粒间静摩擦系数为1.096、颗粒-颗粒间滚动摩擦系数为0.17，不显著影响因素取中值；仿真试验堆积角平均值为42.23°，每组仿真结果与实验结果相对误差均小于5%，总样本相对误差为0.68%。【结论】根据耐火泥颗粒离散元参数标定方法确定的最优参数组合进行仿真试验，获得的堆积角与实验值拟合度高，能够精准预测耐火泥颗粒的实际堆积行为。

关键词：耐火泥颗粒；堆积角；离散元；参数标定

Abstract

Objective By integrating the micro-parameters and morphological characteristics of refractory mortar particles and taking the angle of repose as the response value， a discrete element parameter calibration method is established to provide high-fidelity input parameters for the CFD-DEM optimization design of refractory mortar pumping systems.

Methods Refractory mortar particles with particle sizes of 20， 300， and 550 μm were selected as research objects， and discrete element models of refractory mortar particles with three morphologies， namely quasi‑circular， quasi‑block， and quasi‑conical， were established. Numerical values of the angle of repose were obtained through physical experiments. Significant influencing factors among the parameters to be calibrated were screened out using Plackett-Burman experimental design. The reasonable value ranges of the significant factors were narrowed down by the steepest ascent experiment. Finally， response surface analysis was carried out via Box-Behnken experiments to obtain a second-order regression equation model， and the optimal parameter combination of significant influencing factors reflecting the accumulation characteristics of refractory mortar particles was determined.

Results and Discussion The average angle of repose of refractory mortar particles measured through multiple groups of physical stacking experiments was 41.87°. The results of the Plackett-Burman experiment showed that among the seven parameters to be calibrated， particle surface energy and the rolling friction coefficient between particles were highly significant influencing factors on the angle of repose， the static friction coefficient between particles was a significant influencing factor， and the coefficient of restitution between particles， the coefficient of restitution between particles and the pipe wall， the static friction coefficient between particles and the pipe wall， and the rolling friction coefficient between particles and the pipe wall were insignificant influencing factors. The steepest ascent experiment further narrowed the reasonable value ranges of particle surface energy， the rolling friction coefficient between particles， and the static friction coefficient between particles. The optimal parameter combination was determined by the Box‑Behnken response surface methodology， namely the particle surface energy was 0.004 J/m²， the static friction coefficient between particles was 1.096， and the rolling friction coefficient between particles was 0.17， with the insignificant influencing factors set at their median values. Based on the optimal parameter combination， discrete element simulation verification was performed using the EDEM software， and the average angle of repose obtained from the simulation experiments was 42.23°. The relative error between each simulation result and the experimental result was less than 5%， and the relative error of the total samples was 0.68%.

Conclusion Simulation experiments are carried out using the optimal parameter combination determined by the discrete element parameter calibration method for refractory mortar particles. The resulting angle of repose shows a high fitting degree with the experimental value， enabling accurate prediction of the actual accumulation behavior of refractory mortar particles. This provides high-precision simulation input parameters for the CFD-DEM optimization design of pumping systems.

Keywords： refractory mortar particle； angle of repose； discrete element； parameter calibration

参考文献（References）

［1］吴鹏民. 焦炉绿色智能砌筑发展方向及相关技术探讨［J］. 建筑科技， 2024， 8（1）： 43-46， 53.

Wu Pengmin. Discussion on the development direction of green intelligent masonry for coke ovens and related technologies［J］. Building Technology， 2024， 8（1）： 43-46， 53.

［2］Roussel N， Geiker M R， Dufour F， et al. Computational modeling of concrete flow： general overview［J］. Cement and Concrete Research， 2007， 37（9）： 1298-1307.

［3］Roussel N， Gram A.Simulation of fresh concrete flow： State-of-the Art report of the RILEM technical committee 222-SCF［M］. Dordrecht： Springer Netherlands， 2014..

［4］Sun Zhenjiao， Chen Lianjun， Ma Guanguo， et al. Flow characteristics of moist-mixed materials for shotcrete： from experiment to CFD-DEM simulation［J］. Powder Technology， 2023， 428： 118821.

［5］Coetzee C J， Scheffler O C. Review： the calibration of DEM parameters for the bulk modelling of cohesive materials［J］. Processes， 2023， 11（1）： 5.

［6］Ma Haozhou， Wang Xuewen， Li Bo， et al. Calibration of discrete element microparameters of coal based on the response surface method［J］. Particulate Science and Technology， 2022， 40（5）： 543-557.

［7］Xia Rui， Li Bo， Wang Xuewen， et al. Measurement and calibration of the discrete element parameters of wet bulk coal［J］. Measurement， 2019， 142： 84-95.

［8］Zhan Yijian， Gong Jian， Huang Yulin， et al. Numerical study on concrete pumping behavior via local flow simulation with discrete element method［J］. Materials， 2019， 12（9）： 1415.

［9］Zhang Hong， Zheng Mingxin， Rong Yao， et al. Calibration and sensitivity analysis of macro and meso parameters of discrete element model for coal measure soil［J］. Journal of Physics： Conference Series， 2023， 2519（1）： 012036.

［10］刘佳， 连峰， 刘治， 等. 影响砂砾土堆积角形成的离散元关键参数［J］. 中国粉体技术， 2025， 31（1）： 98-109.

Liu Jia， Lian Feng， Liu Zhi， et al. Key discrete element parameters influencing angle of repose formation of gravel soil［J］. China Powder Science and Technology， 2025， 31（1）： 98-109.

［11］Feng Y T， Owen D R J. Discrete element modelling of large scale particle systems I： exact scaling laws［J］. Computational Particle Mechanics， 2014， 1（2）： 159-168.

［12］He Ping， Fan Yiwei， Pan Banglong， et al. Calibration and verification of dynamic particle flow parameters by the back-propagation neural network based on the genetic algorithm： recycled polyurethane powder［J］. Materials， 2019， 12（20）： 3350.

［13］吴震， 王利强， 徐立敏， 等. 基于静态和动态休止角的钛白粉离散元仿真参数标定［J］. 中国粉体技术， 2023， 29（4）： 108-119.

Wu Zhen， Wang Liqiang， Xu Limin， et al. Discrete element parameters calibration of titanium dioxide based on static and dynamic repose angles［J］. China Powder Science and Technology， 2023， 29（4）： 108-119.

［14］蔡文源， 王利强， 徐立敏. 基于静态和动态休止角的超细碳酸钙离散元参数标定［J］. 中国粉体技术， 2024， 30（4）： 81-93.

Cai Wenyuan， Wang Liqiang， Xu Limin. Discrete elemental parameter calibration of ultrafine calcium carbonate based on static and dynamic angle of repose［J］. China Powder Science and Technology， 2024， 30（4）： 81-93.

［15］韩伟， 王绍宗， 张倩， 等. 基于JKR接触模型的微米级颗粒离散元参数标定［J］. 中国粉体技术， 2021， 27（6）： 60-69.

Han Wei， Wang Shaozong， Zhang Qian， et al. Discrete element parameter calibration of micron sized powder particles based on JKR contact model［J］. China Powder Science and Technology， 2021， 27（6）： 60-69.

［16］GB/T 22459.5—2008 耐火泥浆 第5部分：粒度分布（筛分析）试验方法［S］.

GB/T 22459.5—2008 Refractory mortars - part 5： determination of grain size distribution （sieve analysis）： ［S］.

［17］Li Zhengquan， Chen Huimin， Wu Yukun， et al. CFD-DEM analysis of hydraulic conveying of non-spherical particles through a vertical-bend-horizontal pipeline［J］. Powder Technology， 2024， 434： 119361.

［18］Ma Huaqing， Liu Chang， Wang Wenrui， et al. DEM-FEM investigation of the particle transport process in a flexible tube［J］. Powder Technology， 2025， 455： 120776.

［19］Johnson K L， Kendall K， Roberts A D. Surface energy and the contact of elastic solids［J］. Proceedings of the Royal Society of London A Mathematical and Physical Sciences， 1971， 324（1558）： 301-313.

［20］Chen Xizhong， Elliott J A. On the scaling law of JKR contact model for coarse-grained cohesive particles［J］. Chemical Engineering Science， 2020， 227： 115906.

［21］Asaf Z， Rubinstein D， Shmulevich I. Determination of discrete element model parameters required for soil tillage［J］. Soil and Tillage Research， 2007， 92（1/2）： 227-242.

［22］Shi Chong， Yang Wenkun， Yang Junxiong， et al. Calibration of micro-scaled mechanical parameters of granite based on a bonded-particle model with 2D particle flow code［J］. Granular Matter， 2019， 21（2）： 38.

［23］罗帅， 袁巧霞， Gouda Shaban， 等. 基于JKR粘结模型的蚯蚓粪基质离散元法参数标定［J］. 农业机械学报， 2018， 49（4）： 343-350.

Luo Shuai， Yuan Qiaoxia， Gouda Shaban， et al. Parameters calibration of vermicomposting nursery substrate with discrete element method based on JKR contact model［J］. Transactions of the Chinese Society for Agricultural Machinery， 2018， 49（4）： 343-350.

［24］向伟， 吴明亮， 吕江南， 等. 基于堆积试验的黏壤土仿真物理参数标定［J］. 农业工程学报， 2019， 35（12）： 116-123.

Xiang Wei， Wu Mingliang， Lyu Jiangnan， et al. Calibration of simulation physical parameters of clay loam based on soil accumulation test［J］. Transactions of the Chinese Society of Agricultural Engineering， 2019， 35（12）： 116-123.

［25］Ucgul M， Fielke J M， Saunders C. Three-dimensional discrete element modelling of tillage： determination of a suitable contact model and parameters for a cohesionless soil［J］. Biosystems Engineering， 2014， 121： 105-117.

［26］Alizadeh M， Asachi M， Ghadiri M， et al. A methodology for calibration of DEM input parameters in simulation of segregation of powder mixtures： a special focus on adhesion［J］. Powder Technology， 2018， 339： 789-800.