低剪切条件下中药浸膏粉离散元仿真参数标定方法的优化  

Optimization of Discrete Element Simulation Parameter Calibration Method for Traditional Chinese Medicine Extract Powder Under Low Shear Conditions

作　　者：唐雪芳 李焕正 梁子辰 刘亦菲 刘颖[1] 徐芳芳[2,3] 徐冰[1] TANG Xuefang;LI Huanzheng;LIANG Zichen;LIU Yifei;LIU Ying;XU Fangfang;XU Bing(Beijing Key Laboratory for Production Process Control and Quality Evaluation of Chinese Medicine,Beijing University of Chinese Medicine,Beijing 102400,China;State Key Laboratory on Technologies for Chinese Medicine Pharmaceutical Process Control and Intelligent Manufacture,Nanjing University of Chinese Medicine,Nanjing 210023,China;Jiangsu Kanion Pharmaceutical Co.Ltd.,Lianyungang 222001,China)

机构地区：[1]北京中医药大学,中药生产过程控制与质量评价北京市重点实验室,北京102400 [2]南京中医药大学,中药制药过程控制与智能制造技术全国重点实验室,南京210023 [3]江苏康缘药业股份有限公司,江苏连云港222001

出　　处：《中国实验方剂学杂志》2025年第9期211-218,共8页Chinese Journal of Experimental Traditional Medical Formulae

基　　金：北京中医药大学基本科研业务(揭榜挂帅)项目(2023-JYB-JBZD-060);中药制药过程控制与智能制造技术全国重点实验室开放基金项目(SKL2024Z0205);国家工信部2023年产业基础再造和制造业高质量发展专项(TC2308068)。

摘　　要：目的:提高低剪切条件下离散元法模拟中药粉体体系加工过程的准确性。方法:以通塞脉片浸膏粉和芪葛颗粒浸膏粉为研究对象,采用休止角(AOR)测试法和剪切池测试法分别测定了2种物料的AOR和有效内摩擦角。基于Hertz-Mindlin with JKR V2接触模型和颗粒缩放理论,以颗粒-颗粒恢复系数(A)、颗粒-颗粒静摩擦系数(B)、颗粒-颗粒滚动摩擦系数(C)、颗粒-钢恢复系数(D)、颗粒-钢静摩擦系数(E)、颗粒-钢滚动摩擦系数(F)和Johnson-Kendall-Roberts(JKR)表面能(G)为试验因素,先以AOR为参考值对通塞脉片浸膏粉的仿真接触参数进行单一参考值标定,再以AOR和有效内摩擦角为参考值对通塞脉片浸膏粉和芪葛颗粒浸膏粉的仿真接触参数进行联合标定。然后采用Plackett-Burman试验筛选出对模拟参考值影响显著的关键接触参数,利用最陡爬坡试验确定关键接触参数最优范围,最后通过响应面试验设计建立关键接触参数与模拟参考值之间的回归模型,基于回归模型和满意度函数法对关键接触参数进行标定。结果:单一参考值下通塞脉片浸膏粉离散元接触参数A、B、C、D、E、F和G的最佳组合是0.100、0.718、0.616、0.100、0.400、0.250和0.075 J·m^(-2),经验证,其仿真AOR和仿真有效内摩擦角的相对误差分别为0.10%和-8.64%;联合参考值下通塞脉片浸膏粉离散元接触参数A、B、C、D、E、F和G的最佳组合为0.100、0.682、0.598、0.100、0.521、0.294和0.075 J·m^(-2),经验证,其仿真AOR和仿真有效内摩擦角的相对误差分别为0.10%和-0.18%;联合参考值下芪葛颗粒浸膏粉离散元接触参数A、B、C、D、E、F和G的最佳组合为0.150、0.370、0.330、0.150、0.500、0.500和0.100 J·m^(-2),经验证,其仿真AOR和仿真有效内摩擦角的相对误差分别为2.70%和-1.30%。与单一参考值法比较,联合标定方法不仅增加了表征颗粒-设备相互作用的关键接触参数而且更加准确�Objective:To improve the accuracy of discrete element method in simulating the processing of traditional Chinese medicine(TCM)powder system under low shear conditions.Methods:In this study,extract powders of Tongsaimai tablets and Qige granules were used as the research objects,the angle of repose(AOR)and effective angle of internal friction of the two materials were determined by AOR test method and shear cell test method.Based on the Hertz-Mindlin with JKR V2 contact model and particle scaling theory,taking the particle-particle restitution coefficien(t A),particle-particle static friction coefficient(B),particle-particle rolling friction coefficien(t C),particle-steel restitution coefficien(t D),particle-steel static friction coefficient(E),particle-steel rolling friction coefficient(F)and Johnson-Kendall-Roberts(JKR)surface energy(G)as test factors,the simulated contact parameters of Tongsaimai tablets extract powder were first calibrated with a single reference value using AOR as the reference value,and then the simulated contact parameters of Tongsaimai tablets extract powder as well as Qige granules extract powder were co-calibrated with AOR and effective angle of internal friction as the joint reference value,respectively.Then,Plackett-Burman design was used to screen the critical contact parameters that have a significant effect on the simulated reference value,and the steepest ascent design was used to determine the optimal range of the critical contact parameters,finally,the regression model between the critical contact parameters and the simulated reference values was established through the design of the response surface test,and the critical contact parameters were calibrated based on the regression model and the desirability function approach.Results:The optimal combination of discrete elemental contact parameters A-G for Tongsaimai tablets extract powder under a single reference value was 0.100,0.718,0.616,0.100,0.400,0.250 and 0.075 J·m^(-2),which was validated to have relative errors of 0.10%an

关 键 词：离散元法 休止角 有效内摩擦角 参数标定 关键接触参数 低剪切 响应面优化 

分 类 号：R277[医药卫生—中医学] R943[化学工程—制药化工] TQ460.6