PharmaPy 干法制粒二次开发教程（9）：粉碎模块①——磨机类型、筛网角色与 Vogel–Peukert 破碎概率核

PharmaPy 干法制粒二次开发教程9粉碎模块①——磨机类型、筛网角色与 Vogel–Peukert 破碎概率核版本声明块工具/软件本项目二次开发的破碎动力学函数本系列二次开发示例非官方 API基于 Vogel–Peukert2003破碎概率核语言/环境Python 3.9 numpy/scipyscipy.special.expm1用于数值稳定本文目标先认清磨机是干什么的、筛网在中间扮演什么——再落地可复用的破碎概率函数并做一次单粒 Monte Carlo 演示感受一粒能不能碎是随机事件。一句话结论制药整粒用低能磨摇摆/锥式/锤式把辊压得到的 ribbon 打回颗粒筛网孔径定上限粒径而叶轮尖端速度主要决定细粉量Vogel–Peukert2003把破碎当作随机事件单颗粒破碎概率S 1−exp{−f_Mat·k·x·(W_m,kin−W_m,min)}其中x为初始粒径、W_m,min为物料破碎阈值能——本项目将其实现为纯函数并用 Monte Carlo 演示接口接续篇 8 的ribbon_psd。〇、本篇要解决的认知问题为什么整粒都用低速低能磨筛网到底是滤掉大小还是决定破碎程度摇摆、锥式、锤式三种磨各自的破碎机制是什么粒度控制变量分别是什么破碎为什么是个概率而非必然事件Vogel–Peukert 核里每一项的物理含义是什么粗看粉碎概率像公式可它到底怎么和颗粒粒径分布PSD对接破碎后子粒去哪了篇 8 给我交出来的 ribbon PSD我要怎么接进来作为粉碎模块的输入一、机制解析磨机只做一件事把大颗粒变成小颗粒筛网决定能多碎多少先破除一个常见的误解筛网不是挡筛网配筛孔大小剩下的自己掉下去这么简单。在锥式/摇摆磨里物料是被转子或摆动锤强行压过筛网孔破碎主要发生在过筛挤压 剪切这一下撞击里而不是筛网后端选粒径。所以筛网的真实角色是同时决定上限粒径几何与破碎程度力学孔小了颗粒要过孔就得先被撞碎到能穿过去于是平均粒径和细粉都跟着变。这也是为什么筛网孔径 转子速度两大变量可以解耦调节孔径定能过的最大尺寸转速定每次撞击的能量进而定破碎概率。制药整粒的工业经验高度一致为了少产生细粉应优先用低速的低能磨格转子转速通常1000 rpmHancock Vendola 2008 综述比较了摇摆、锥式、锤式三类结论见底账 C1。细粉的定义本系列统一为125 µm质量分数见篇 16。三种磨的破碎机理差异底账 C1 汇总磨机核心机理粒度控制变量特点摇摆式颗粒机oscillating摆动转子把物料挤出/压过筛网筛网、转子速、摆角作用温和适合脆性 ribbon锥式整粒机conicalQuadro Comil 类锥形筛网内旋转楔形叶轮剪切推进筛网尺寸/类型、叶轮尖端速度给出一致上限粒径、收窄 PSD锤式磨hammerFitzMill 类高速锤/刀冲击把物料击碎锤型、转速、筛网对脆性料效率高细粉偏多其中锥式整粒机Comil 类在制药 RC-milling 里用得最多因为它在一致性上限粒径 收窄 PSD上表现好。决定 Comil 出料 PSD 的最大因素 叶轮尖端速度 筛网孔径/几何Hancock Vendola 2008底账 C1。这给了我们一个工程指向做粉碎动力学模型时筛网尺寸进上限粒径叶轮尖端速度进破碎能量核。破碎为什么是概率一粒颗粒被压过筛孔或撞上锤头它有某个概率碎、也可能不碎。碎片数量、大小也随机。把这种随机性数学化的最小框架就是破碎概率核给定初始粒径x与单次碰撞的比冲击动能W_m,kin算这一粒碎不碎/碎多少的概率。这个概率再乘上单位时间的碰撞次数强度就是后面篇 10 破碎 PBM 里的**选择函数selection function**雏形。所以篇 9 这一步的价值是先把单次事件概率这一个最原子的问题讲清楚为篇 10 的群体化铺路。Vogel–Peukert2003物料破碎特性不是一个常数而是一组属性Vogel 与 Peukert 把多种材料的单颗粒冲击破碎数据拟合出一条主曲线mastercurve破碎概率统一写作底账 C3S 1 − exp{ − f_Mat · k · x · (W_m,kin − W_m,min) }S破碎概率0~1无量纲。x初始粒径 [m]。k冲击次数。W_m,kin质量比冲击动能 [J/kg] 单粒动能/质量。W_m,min材料破碎阈值能[J/kg]低于它几乎不碎。f_Mat物料破碎特性参数“易碎性”单位 [m·(J/kg)⁻¹·s⁻¹] 量纲上配合 k 与 x 使指数无量纲。指数里指数项写成expm1exp(x)−1当W_m,kin刚超过阈值时数值更稳。关键直觉x越大、f_Mat越大、超过阈值越多破碎概率就越接近 1——大颗粒在同样动能下更容易被撞碎这正对应破碎概率随粒径增大的物理事实。需要强调的边界真实性铁律f_Mat、W_m,min具体数值必须针对物料实测/标定不同材料差异很大本文不编造一组通用常数演示仅用示意参数底账 J 未核实项Vogel–Peukert 曾被改造成描述冲击破碎程度t_10ACS Engineering Au 2023本系列只用其概率形式。与篇 8 的交接点ribbon_psd作为粉碎输入篇 8 的RollerCompactor交出了ribbon_psd进料 ribbon 的x/distrib网格。粉碎模块的输入正是这一份 ribbon PSD——把每一档粒径的颗粒当作一批母粒各自按S(x)的概率破碎。本系列用工程选型见铁律取几何网格~30–50 bin本篇起统一在x_gr np.geomspace(10, 2000, num40)上工作。二、完整代码与逐行剖析下面把篇 8 的 ribbon PSD 接进来实现vogel_peukert_probability(...)含数值稳定的expm1再用一次单粒 Monte Carlo演示碎与不碎是随机事件、但统计上随粒径/动能单调上升。整段都是可运行的示例f_Mat/W_m,min为示意参数正文注明以实测/文献为准。importnumpyasnpfromscipy.specialimportexpm1# # 破碎概率核Vogel Peukert (2003) 示意参数以实测/文献为准# S 1 - exp{ -f_Mat * k * x * (W_m,kin - W_m,min) }# expm1(u)exp(u)-1 提升 W_m,kin 刚过阈值时的数值稳定性# defvogel_peukert_breakage_probability(x,W_m_kin,k1.0,f_Mat1.5e4,W_m_min8.0):单颗粒破碎概率。 x : 初始粒径 [m] W_m_kin : 质量比冲击动能 [J/kg] k : 冲击次数默认1 f_Mat : 物料破碎特性参数 [示意] W_m_min : 材料破碎阈值能 [J/kg] [示意] ifW_m_kinW_m_min:return0.0# 低于阈值概率恒为0exponent-f_Mat*k*x*(W_m_kin-W_m_min)return-expm1(exponent)# 1 - exp(u) 用 expm1(u)exp(u)-1 -expm1(-a)1-e^{-a}# ---- 接续篇8从 ribbon_psd 取进料逐档算破碎概率 ----defattach_ribbon_psd(x_gr,distrib_ribbon,W_m_kin15.0):把篇8的 ribbon PSDx/distrib 网格变成每档的破碎概率向量。x_cnp.sqrt(x_gr[:-1]*x_gr[1:])# bin 中心几何中点µmx_mx_c*1e-6# µm - mSnp.array([vogel_peukert_breakage_probability(xi,W_m_kin)forxiinx_m])# 归一保留母粒总质量份额参考 - 篇10 PBM 里做质量守恒M_totalnp.sum(distrib_ribbon)# 示意粒子数密度归一returnx_c,S,M_total# ---- 单粒 Monte Carlo 演示统计概率随粒径/动能的单调性 ----rngnp.random.default_rng(42)N_trials20000forWin[10.0,15.0,30.0]:brokennp.zeros(3)x_samplenp.array([100e-6,500e-6,1500e-6])# 0.1/0.5/1.5 mmforj,xiinenumerate(x_sample):svogel_peukert_breakage_probability(xi,W)broken[j]np.mean(rng.random(N_trials)s)# 蒙特卡洛统计print(fW_m,kin{W:5.1f}J/kg 破碎概率估计: f{broken[0]:.4f}/{broken[1]:.4f}/{broken[2]:.4f}(小/中/大粒))运行趋势应满足可作验收依据同一W下粒径越大概率越高broken[2] ≥ broken[1] ≥ broken[0]同一粒径下W越大概率越高。这些行为断言正是篇 10 里做 PSD 演化正确性的伏笔。逐行要点expm11−exp(u)里的1−exp(u)直接写会在大负数时舍入-expm1(-a) 1−e^{-a}更稳。工程上可以不用但这是数值稳定的好习惯呼应篇 10 的质量守恒。attach_ribbon_psd把篇 8 的ribbon_psd转成每档破碎概率。x_c用几何中点sqrt(x_gr[:-1]*x_gr[1:])因为对数网格下几何中点才是 bin 代表粒径呼应底部网格共轭铁律。Monte Carlo 那段用rng.random(N_trials) s统计频率逼近概率证明了S的概率语义——这与后面篇 10 的确定性破碎 PBM用S_i作速率系数是同一物理的两种视角先建立直觉。三、常见报错与排查概率算出来超过 1 或为负。现象vogel_peukert_breakage_probability返回 1.3 或 -0.5。根因W_m_kin − W_m_min是负的而你直接1−exp(正数)溢出成大于 1或f_Mat*x数字太大导致 exp 溢出。解法先判W_m_kin W_m_min直接返回 0本篇已内置并检查量纲——x必须以米别把 µm 直接代入。Monte Carlo 频率和理论概率对不上尤其大粒径。现象broken[2]明显偏离解析S。根因N_trials太小尾部大粒度概率接近 1随机波动被放大或每次循环都新建rng导致可复现性丢失。解法固定seed、N_trials提到2e4且对p≈1的档位不要用贝努利估计而直接读S本身。ribbon PSD 尺寸/量纲错导致概率全为 0 或全为 1。现象x_c用 µm 但函数内部按米算所有概率塌缩到某一边。根因attach_ribbon_psd里的x_m x_c*1e-6没生效或x_c已经错了。解法打印x_c的量级核对µm 应在 10~2000确认换算后才进核函数并在函数里加assert np.all(x_m0)。筛网上限粒径没用上碎出来的大粒超出筛孔。现象破碎后仍能看到 筛孔的颗粒。根因篇 9 只给了概率核还没接筛网几何约束——那是篇 11选择/破碎函数 筛网分级d_i和篇 12Tromp 分级的职责。解法明确本篇边界不要在本篇强行塞筛网把上限粒径作为输入留到篇 11 的网格内分类处理。四、动手练习概率单调性验收跑上面 Monte Carlo 段把x_sample换成 5 档50/200/500/1000/2000 µmW_m_kin[10,20,40]分别打印概率矩阵。判定按行看粒径升→概率升按列看动能升→概率升任一档位违例即为失败同时解析S与 Monte Carlo 频率偏差0.02。阈值能边界设W_m_kin恰等于W_m_min与W_m_min*1.01打印概率。判定恰好等于阈值时为0略超阈值时概率从极小值平滑起步用expm1稳定证明没有跳变。量纲自检把x手误传入 µm 版本观察概率塌缩。判定加一行assert报错让你抓住 µm/m 未换算这类单位显式换算铁律的违例。五、小结与下一篇预告本篇把粉碎战场铺开三种低能磨机理与筛网角色并把破碎的随机本质用 Vogel–Peukert 概率核 单粒 Monte Carlo 讲透。你应记住三个可迁移的心智筛网同时管上限粒径与破碎程度、Comil 类锥磨用叶轮尖端速度定细粉、破碎概率S(x,W_m,kin)是单次事件视角——它即将被剂伯爵化成篇 10 破碎 PBM 里决定 PSD 演化的选择函数。下一篇篇 10《粉碎②离散破碎 PBM 数值实现》——把单个颗粒的碎不碎升级为整批 PSD 如何随停留时间演化落籽 Hounslow 几何网格与 Kumar–Ramkrishna1996固定枢轴法用solve_ivp积分dN_i/dtΣ_{j≤i}S_j·B_ij·N_j − S_i·N_i并以质量守恒容差作为验收硬标准。本篇认知问题回显FAQQ1为什么要用低速低能磨筛网是滤大小还是定破碎A1制药整粒为少产细粉细粉125µm本系列口径而倾向低速低能磨常用转子1000 rpm。筛网不是单纯挡大小——在锥式/摇摆磨里物料是被转子强行压过筛孔所以筛网同时决定上限粒径几何与破碎程度力学同理叶轮尖端速度主要决定细粉量两者可解耦调节。Q2摇摆/锥式/锤式磨有什么机理差异A2摇摆磨靠摆动转子把物料挤出/压过筛网温和、适合脆性 ribbon锥式磨Comil 类靠锥形筛网内旋转叶轮剪切推进粒径控制最稳、能收窄 PSD工业用得最多锤式磨FitzMill 类靠高速锤头冲击对脆性料效率高但细粉偏多。决定 Comil 出料的最大因素是叶轮尖端速度 筛网孔径/几何。Q3破碎为什么是概率Vogel–Peukert 核每项是什么A3单颗粒撞上锤头碎不碎、碎成几块本质随机所以用概率描述。Vogel–Peukert2003写为S1−exp{−f_Mat·k·x·(W_m,kin−W_m,min)}x初始粒径、k冲击次数、W_m,kin单粒质量比冲击动能、W_m,min材料破碎阈值能、f_Mat物料破碎特性参数。指纹在指数项——粒径越大、动能超过阈值越多、材料越易碎概率越接近 1。Q4单粒概率怎么和 PSD 对接A4把 ribbon PSD 的每一档颗粒当作一批母粒各自按S(x)的概率破碎碎出的子粒再落回更小的档位——这就是篇 10 破碎 PBM 的出生/死亡结构Σ_{j≤i}S_j·B_ij·N_j − S_i·N_i。本篇的S(x,W_m,kin)就是决定演化的选择函数雏形。Q5篇 8 的 ribbon PSD 怎么接进粉碎模块各自按S(x)的概率破碎碎出的子粒再落回更小的档位——这就是篇 10 破碎 PBM 的出生/死亡结构Σ_{j≤i}S_j·B_ij·N_j − S_i·N_i。本篇的S(x,W_m,kin)就是决定演化的选择函数雏形。Q5篇 8 的 ribbon PSD 怎么接进粉碎模块A5用attach_ribbon_psd(x_gr, distrib_ribbon, W_m_kin)读篇 8 的ribbon_psd按 bin 几何中心sqrt(x_gr[:-1]*x_gr[1:])把粒径换算成米后再逐档算破碎概率返回(x_c, S, M_total)。M_total用于篇 10 的质量守恒核对x_c用几何中点是对数网格的 bin 代表粒径网格共轭铁律。