企业官网建设流程全解析

直接压片单元模型开发教程（14）：力传递比与弹出力——上下冲力、模壁摩擦与压制功

版本声明块

• 工具/软件：本系列自研TabletPress单元类的力学表征模块（本系列二次开发示例类，非官方 API）；表征量清单依 USP-NF⟨1245⟩，力传递比定义依 Hussain(1991)，实测锚点依 Roopwani(2019)。
• 语言/环境：Python 3.9 + numpy（<2）+ scipy（<1.14）+ matplotlib
• 本文目标：把"力从上面打进去、有多少传到下面、脱模又要花多少力"这条链条算清楚，量化润滑剂与模壁摩擦的作用。

一句话结论：力传递比定义为R = 最大下冲力 / 最大上冲力（Hussain 1991），Roopwani(2019) 实测：未处理加巴喷丁77.7%±1.2、共处理85.8%±2.5、外加润滑93.5%±0.8，同时三者弹出应力分别为 13.9 / 7.1 / 1.1 MPa（随共处理/润滑降低）；力在上冲与下冲之差消耗于模壁摩擦；USP-NF⟨1245⟩ 把轴向/径向应力、残余模壁应力、压制功（work of compaction）、弹出曲线与 SRS 列为标准表征量。硬脂酸镁常用0.5–1.5% w/w（★ 行业站，以官方文档为准）。

〇、本篇要解决的认知问题

1. 力传递比R = 最大下冲力/最大上冲力为什么永远小于 100%？差值去哪了？
2. 为什么"未处理加巴喷丁 77.7%"、“共处理 85.8%”、"外加润滑 93.5%"三种状态 R 逐个升高？
3. 弹出力/弹出应力为什么与 R 反着走（R 高时弹出应力低到 1.1 MPa）？
4. 硬脂酸镁加多少？加多了除了脱模变好，还会不会带来坏处？
5. 怎么把上下冲力、摩擦功、弹出力与压制功写成一个能量守恒、能自检的可运行模块？

一、机制解析

1.1 力在上冲与下冲之间"漏"掉了

一片药从上面受压、下面顶着，直觉上会以为"上下受力相等"。但实测不是——下冲感受到的峰值力总小于上冲，因为中间隔着一层模壁摩擦：片剂被压紧后侧向膨胀、紧紧顶住冲模壁，粉体/片剂沿模壁下滑时摩擦力把一部分轴向力"吃掉"了。底账明确把这个比值定义为：

R = 最大下冲力 / 最大上冲力 （Hussain 1991）

于是有：

F_upper ──►[ 片剂 ]──► F_lower │ └─ 模壁摩擦把 (F_upper − F_lower) 消耗掉 R = F_lower / F_upper < 1 （R 越大，力传递越充分、摩擦损失越小）

一个自洽的能量视角：上冲做的功 = 下冲接收的功 + 模壁摩擦耗散的功（忽略其他损失）。这就是为什么"摩擦功"可以从上下冲力差积分出来——本篇代码正是这么做的。

1.2 实测锚点：三种表面状态，三档 R

Roopwani(2019) 的加巴喷丁实验给了三条可直接引用的锚点，是理解 R 的最佳教材：

三条数据讲了一个完整故事：R 与弹出应力同向于"摩擦大小"，但方向相反——摩擦越小，R 越高（越接近 1），弹出应力越低（越容易脱模）。外加润滑能把 R 提到 93.5%、把弹出应力压到 1.1 MPa，说明润滑对"力传递效率 + 脱模"是双重收益。

注意"共处理 vs 外加润滑"的机制不同：共处理是在颗粒层面改变表面/塑性（内润滑或表面改性），外加润滑是在冲模壁上喷涂润滑剂。两者都降低摩擦，但一个是材料属性、一个是过程操作——放大时前者随处方走、后者随设备走，不可互换。

1.3 弹出力与模壁摩擦：脱模那一下要花多少力

弹出力（ejection force）的物理来源很清晰（底账）：脱模时片剂与冲模壁摩擦。影响因素列为四类：模壁摩擦、润滑剂浓度、片剂密度、材料变形类型。把弹出应力量化，就是上面那张表的第三列：13.9 / 7.1 / 1.1 MPa。

从应力到力，需要接触面积。片剂与模壁的接触面积可近似为圆柱侧面积：

A_contact = π · D · t （D=片径, t=片厚） F_eject = σ_eject · A_contact

这是一个几何上自然的近似（侧面积），配合实测的σ_eject就能算出脱模力。但必须声明：弹出力的精确建模依赖残余模壁应力与摩擦系数，具体形式以官方文档/原文为准——本篇把它作为工程近似使用，用于比较不同润滑条件下的相对大小（这恰恰是它的合理用途）。

润滑剂方面，底账记录硬脂酸镁（magnesium stearate）常用 0.5–1.5% w/w（★ 行业站）。这属于常见范围，正式使用须以官方文档为准。还要记住一个工程常识：润滑剂（尤其疏水的硬脂酸镁）加多会降低片剂强度、延缓崩解——所以 R 升高的收益是有代价的，"把润滑加到最大"不是最优解。

1.4 USP-NF⟨1245⟩ 的力学表征清单

USP-NF⟨1245⟩ 明确列出压缩表征涉及的力学量，本篇覆盖其中相当一部分：

• 轴向应力 / 径向应力：轴向即上下冲方向的应力，径向即侧向、与模壁摩擦相关；
• 残余模壁应力：脱模瞬间仍留在模壁上的应力，直接关联弹出力；
• 压制功（work of compaction）：上冲力沿位移积分得到的能量；
• 弹出曲线（ejection curve）：脱模过程中力随位移/时间的曲线；
• 应变速率敏感性（SRS）：与 dwell time、转速耦合（第 7、18 篇）。

把这些量对齐到模型里，就得到本篇的"力-能量流水线"：

F_upper, y(t) F_upper − F_lower 上冲力-位移曲线 ──► 压制功 W=∫F dy │ │ ▼ │ 模壁摩擦功 W_fric = ∫(F_upper−F_lower)dy ▼ R = F_lower/F_upper ──► 力传递效率诊断 │ ▼ 弹出力 F_eject = σ_eject · (πDt) ──► 脱模难易 / 润滑效应

1.5 与前后篇的连接

第 13 篇讲"两段压力怎么把片子压低"，本篇讲"这个力在上下冲之间怎么分配、脱模要花多少"——两者共同构成压片机的"力平衡"。第 15 篇的片重会改变受压时的片剂密度，从而改变模壁摩擦与 R；第 19 篇 QbD 则把润滑水平（或 R）与破碎力一并作为响应变量。所以 R 不只是"诊断指标"，还是可标定的模型参数。

1.6 为什么 R 是有用的诊断量：它把"摩擦"变成可读数字

力传递比最实用的价值，是把看不见的"模壁摩擦"压缩成一个可比较、可追踪的数字。它的用法可以归纳成三条：

1. 处方/表面状态的指纹：同一台机器、同一个压力，R 从 77.7% 跳到 93.5%，一定是摩擦变了（共处理或外加润滑）。R 的稳定性因此可以当"批次一致性/润滑均匀性"的在线征兆。Roopwani(2019) 三档值（77.7%±1.2、85.8%±2.5、93.5%±0.8）正是这种"指纹"的标定样本。
2. 设备与模具状态的告警：模壁磨损、冲头/冲模配合变差、清洁残留都会改变摩擦，进而改变 R。持续监控 R 的漂移，比事后发现"脱模力升高、片剂顶裂"要早。
3. 放大与转移的比较基线：从压片机 A 换到 B，若同一处方在相同 dwell time 下 R 明显不同，说明两者的模壁摩擦/几何状态不可直接类比，工艺窗口不能原样平移。

要把 R 用成诊断量，必须先把测量条件说清楚：冲头直径、压力水平、转速（关联 dwell time）、片重与片厚。否则两个 R 之间没有可比性——这也是本系列反复强调的"元数据与数值同行"。

1.7 三类功与两条曲线：USP-NF⟨1245⟩ 的力学表征怎么读

USP-NF⟨1245⟩ 把一系列力学量列为压缩表征内容，本篇把其中最核心的几条串成"能量账本"：

读这张表的要点是：压制功是"打进去"的能量，摩擦功是被模壁消耗的部分，弹出功/弹出力是"拿出来"的代价。三者量纲都是"力×位移"或"力"，但反映的是压缩循环三个不同阶段。把它们分开记，才能诊断问题出在哪里——例如"压制功正常但弹出力异常高"往往指向脱模侧（模壁摩擦、润滑不足、片剂密度过高），而不是主压本身。

一个容易被忽略的关系是：R 越低，说明同一上冲力下传下去的力越少，摩擦功占比越高（因为摩擦功正比于F_upper − F_lower）。这把 R 与能量账本直接挂钩：R 不只是"力的比值"，还是"能量在哪被吃掉"的指示器。

1.8 润滑的收益与代价：为什么 0.5–1.5% 是"折中区间"

底账记录硬脂酸镁常用 0.5–1.5% w/w（★ 行业站，以官方文档为准）。理解这个区间为什么存在，要同时看收益与代价两面：

• 收益：润滑降低模壁摩擦 → R 升高、弹出力/弹出应力下降 → 顶裂、脱模困难的风险降低。Roopwani(2019) 的外加润滑档（R=93.5%、σ_eject=1.1 MPa）展示了收益上限。
• 代价：润滑剂（尤其疏水的硬脂酸镁）是"惰性"成分，会降低片剂强度（削弱颗粒间键合），高用量还会延缓崩解（第 11 篇的紧凑性/成形性也因此受扰）。

所以工程上不是"润滑加到最大"，而是在"脱模够顺 + 片强够高 + 崩解达标"三者之间取折中，0.5–1.5% 就是这样一种常见折中区间。注意它属于★级行业常见范围，实际用量必须以处方研究/官方文档为准，不能把它当"合规默认值"。

1.9 压制功、弹出功与应变速率敏感性：能量视角的收束

USP-NF⟨1245⟩ 把SRS（应变速率敏感性）与压制功、弹出曲线并列为表征量。用能量视角可以把它们串成一条主线：

• 压制功（∫F dy）随转速/加载速率变化：加载越快，材料可变形的时间越短，塑性变形减少，达到同一密度往往需要更大的力——这正是 SRS 的表现（第 7 篇的 dwell time、第 18 篇的转速放大深挖）。
• 弹出功/弹出力反映脱模阶段的摩擦代价，与润滑和片剂密度强相关。
• 两者合起来给出"这一片消耗多少能量、其中多少耗在摩擦/脱模上"的诊断。

工程含义是：比较两台机器、两个转速下的压制功与弹出力，必须在相同 dwell time 或相同加载速率下进行，否则差异里混入了 SRS，无法判断是设备差异还是速率差异。把这层口径写进模型文档，是让R、W、F_eject这些量"可比"的前提。

还要强调一点能量账本的用法：W_comp（压制功）应等于"下冲接收的功 + 摩擦耗散的功 + 弹性储能（卸压时部分回弹释放）"的近似平衡。本模型只显式算了前两项，把弹性回弹归入误差项——所以当W_fric/W_comp异常大（例如超过物理合理范围）时，往往说明 R 数据有误或存在未建模的损失，应回查测量，而不是相信模型输出。

二、完整代码与逐行剖析

下面的脚本把上冲力-位移曲线积分成压制功、用上下冲力算 R、用实测弹出应力算弹出力，并对照三档实测锚点做诊断。

importnumpyasnpfromscipy.integrateimporttrapezoid# scipy>=1.6：trapezoid；旧版为 np.trapzimportmatplotlib.pyplotasplt plt.rcParams["font.sans-serif"]=["SimHei","Microsoft YaHei","DejaVu Sans"]plt.rcParams["axes.unicode_minus"]=False# ============================================================# 1. 力传递比 R = 最大下冲力 / 最大上冲力（Hussain 1991）# ============================================================defforce_transmission(F_upper_max,F_lower_max):"""返回 R（无量纲）。输入单位任意但须一致（推荐 N）。"""ifF_upper_max<=0:raiseValueError("上冲力必须为正")returnF_lower_max/F_upper_max# ---- 用 Roopwani(2019) 三档实测锚点自检 ----ANCHORS={"未处理加巴喷丁":dict(R=0.777,sd=0.012,sigma_eject=13.9),"共处理":dict(R=0.858,sd=0.025,sigma_eject=7.1),"外加润滑":dict(R=0.935,sd=0.008,sigma_eject=1.1),}F_upper_demo=10000.0# N（演示）forname,ainANCHORS.items():F_lower=a["R"]*F_upper_demoprint(f"{name}: R={a['R']*100:.1f}%±{a['sd']*100:.1f}% "f"→ 下冲力≈{F_lower:.0f}N（上冲{F_upper_demo:.0f}N）")# ============================================================# 2. 上冲力-位移曲线 → 压制功 W = ∫F dy（trapezoid 数值积分）# 合成一条“加载-卸载”示意曲线（演示；位移-时间曲线以第 4 篇为准）# ============================================================y_up=np.linspace(0.0,4.0,200)# 位移 mm（压缩段）F_up=10000.0*np.sin(np.pi*y_up/(2*4.0))# 示意：平滑升到峰值再回落W_comp=trapezoid(F_up,y_up)# N·mm = mJprint(f"\n压制功（上冲）W ≈{W_comp:.1f}N·mm ={W_comp/1000:.2f}J")# ============================================================# 3. 模壁摩擦功：由上下冲力差积分# F_lower(y) = R · F_up(y)；摩擦功 = ∫(F_up − F_lower)dy# ============================================================deffriction_work(y,F_up,R):"""返回 (摩擦功, 下冲力数组)。R 为力传递比。"""F_low=R*F_upreturntrapezoid(F_up-F_low,y),F_low W_fric,F_low=friction_work(y_up,F_up,R=0.777)print(f"模壁摩擦功（R=0.777）≈{W_fric:.1f}N·mm（占压制功{W_fric/W_comp*100:.1f}%）")# ============================================================# 4. 弹出力：F_eject = σ_eject · (π·D·t)（工程近似，以官方/原文为准）# ============================================================defejection_force(sigma_eject_MPa,D_mm,t_mm):"""弹出力 N。σ_eject 为实测弹出应力；A_contact = π·D·t 为侧面积近似。"""A_contact=np.pi*D_mm*t_mm# mm²returnsigma_eject_MPa*A_contact# ND_mm,t_mm=10.0,4.0forname,ainANCHORS.items():Fe=ejection_force(a["sigma_eject"],D_mm,t_mm)print(f"{name}: 弹出应力={a['sigma_eject']}MPa → 弹出力≈{Fe:.1f}N")# ============================================================# 5. 润滑对 R 的工程近似：R 随润滑水平单调升并饱和（示意见，待标定）# ============================================================defR_from_lubricant(c_pct,R0=0.777,c_ref=1.0):"""c_pct: 硬脂酸镁 % w/w（常用 0.5–1.5，以官方文档为准）。工程近似。"""returnR0+(1.0-R0)*(1.0-np.exp(-np.maximum(c_pct,0.0)/c_ref))c_grid=np.linspace(0.0,1.5,50)R_grid=R_from_lubricant(c_grid)print(f"\n润滑 0.0%→R={R_from_lubricant(0.0)*100:.1f}% ; "f"1.5%→R={R_from_lubricant(1.5)*100:.1f}%（工程近似，以官方文档为准）")# ============================================================# 6. 绘图：R–润滑曲线 + 三档锚点# ============================================================plt.figure(figsize=(7.4,5.0))plt.plot(c_grid,R_grid*100,lw=2,color="#2980b9",label="R(c) 工程近似")forname,ainANCHORS.items():plt.axhline(a["R"]*100,ls=":",lw=1.1,color="#7f8c8d")plt.text(1.52,a["R"]*100,f"{name}{a['R']*100:.1f}%",va="center",fontsize=9,color="#34495e")plt.xlabel("硬脂酸镁 % w/w（常用 0.5–1.5，以官方文档为准）")plt.ylabel("力传递比 R [%]")plt.title("润滑水平与力传递比：R 单调升并饱和（工程近似，锚点来自 Roopwani 2019）")plt.legend(frameon=False,loc="lower right");plt.tight_layout()plt.savefig("force_transmission.png",dpi=150)plt.show()

逐段剖析：

• force_transmission：定义直白，但加了"上冲力必须为正"的守卫——现实数据里若传感器符号接反（下冲力为负），R 会算成负值而不报错。这类物理合理性守卫在力学科目里很值钱。
• ANCHORS字典：把三档实测值作为变量而非硬编码，既方便打印诊断，也方便后面把仿真 R 与实测 R 直接对比。注意每次引用都带±sd——统计口径要随数值一起走。
• trapezoid：∫F dy用梯形法。W = trapezoid(F_up, y_up)单位为 N·mm = mJ。注意scipy.integrate.trapezoid与旧版np.trapz的差异——scipy<1.14 里两者都在，但新代码建议用 scipy 版本。用示意正弦曲线只是为了演示积分流程，真实曲线取实测力-位移数据。
• friction_work：基于"下冲力 = R×上冲力"的线性近似，把摩擦功写成上下冲力差对位移的积分。它直观展示了"R 越低、摩擦功占比越高"——这就是 R=77.7% 时摩擦吃掉近 1/4 能量的定量表达。
• ejection_force：用侧面积πDt把实测弹出应力换算成力。注释写明"工程近似，以官方/原文为准"，避免被当成普适公式。三档对比能直接看出"润滑把脱模力打了下来"。
• R_from_lubricant：把"润滑↑→摩擦↓→R↑"做成饱和指数，起点锚到 0.777（未处理）。它保证 R<1 且单调——物理上合理。c 的单位 % w/w，注释标注 0.5–1.5% 的常见范围与"以官方文档为准"。
• 绘图：把三档实测 R 画成水平参考线，与模型曲线放在一起，便于看出模型在什么润滑水平能覆盖到"外加润滑 93.5%"。

反直觉默认值与坑（表格）

三、常见报错与排查

现象一：force_transmission返回负值或 >1。
根因：上下冲力符号约定不一致（一个传感器正压为负）或单位不同。解法：加"上冲力>0"守卫；对数据先做符号统一（向下压为正），并画一次上下冲力时间序列目视核对；R>1 在物理上不可能，应视为数据错误。

现象二：NameError: name 'trapezoid' is not defined或np.trapz 弃用警告。
根因：scipy 版本差异——scipy.integrate.trapezoid在较新版本才有，旧代码用np.trapz；而 numpy<2 环境下np.trapz可用但属于旧接口。解法：优先from scipy.integrate import trapezoid，若报错改np.trapz；无论如何参数顺序必须是(被积函数数组, 自变量数组)，写反会把 W 算成完全错误的值。

现象三：压制功数值比预期大/小几个数量级。
根因：单位混用——y 用 mm 而力用 N 得 N·mm（=mJ），若误以为 J 就会差 1000 倍；或把 kN 当 N。解法：显式记录"F 用 N、y 用 mm、W 单位 N·mm"，并打印W/1000换算成 J；底账反复强调 kN（力）与 MPa（应力）不可混用，此处同理。

现象四：弹出力与文献对不上。
根因：弹出力是过程量，随残余模壁应力、摩擦系数、片型变化；本篇用πDt侧面积是工程近似。解法：把弹出力用于同类条件下的相对比较（如比较不同润滑），不做跨设备绝对断言；精确形式写"以官方文档/原文为准"。

现象五：把"外加润滑 93.5%"当处方目标 R。
根因：混淆"共处理"与"外加润滑"两种降摩擦机制，且忽略高润滑对片强/崩解的副作用。解法：明确 R 是"表面状态 × 处方 × 设备"的函数；把 93.5% 作为该实验条件下外加润滑的实测值，工艺设计时结合片剂强度与崩解综合取润滑水平。

四、动手练习

1. 摩擦功占比自证：固定F_up与y，把 R 从 0.777 扫到 0.935，算W_fric/W_comp。判定：R 越高、摩擦功占比越低（趋近 0），且 R=0.777 时占比明显更大——你能用"力传递越充分、摩擦损耗越小"解释。
2. 弹出力三档对比：用ejection_force对 13.9/7.1/1.1 MPa、D=10 mm、t=4 mm 算力。判定：三个力随 σ_eject 成比例，外加润滑那档最小；你能说明这与 R 最高（93.5%）是同一物理（摩擦低）的两面。
3. 润滑敏感性：改c_ref（饱和速度）重画 R–c 曲线，看模型在 c=1.5% 时能否接近 93.5%。判定：曲线单调有界（R<1）、随 c_ref 减小更快饱和；你能解释c_ref是待标定参数而非材料常数。

五、小结与下一篇预告

本篇把压片机的"力平衡"算通了：力传递比R = 最大下冲力/最大上冲力（Hussain 1991）永远小于 1，差额消耗于模壁摩擦；Roopwani(2019) 的三档锚点（77.7%±1.2 / 85.8%±2.5 / 93.5%±0.8）显示 R 与弹出应力（13.9/7.1/1.1 MPa）反着走，降摩擦（共处理或外加润滑）双赢；硬脂酸镁常用 0.5–1.5% w/w（以官方文档为准），但要权衡片强与崩解。你还把 USP-NF⟨1245⟩ 列出的压制功、弹出曲线等表征量落成了∫F dy、σ_eject·πDt的可运行实现。

力与密度讲完，接下来换一条线：片剂到底装了多少料——这直接决定含量与主压缩力的起点。

下一篇：第 15 篇《片重差异与模具填料——die fill 模型与 USP⟨905⟩ 均匀度判定》——片重模型W_t = f(ρ_b, D, DP, TS, c)、转台速度比桨叶速度更主导填料、存在最佳中等转速（90 rpm 均匀性下降），以及 USP⟨905⟩ 的 AV 判定器。

本篇认知问题回显（FAQ）

Q1：力传递比为什么永远小于 100%？差值去哪了？
A1：R = 最大下冲力/最大上冲力（Hussain 1991）。片剂被压紧后侧向膨胀顶住模壁，沿模壁下滑的摩擦把一部分轴向力消耗掉，故下冲峰值力总小于上冲，R<1。从上冲做功减去下冲接收的功，就是模壁摩擦耗散的能量，可用∫(F_upper−F_lower)dy积分得到。

Q2：为什么未处理加巴喷丁、共处理、外加润滑三档 R 逐个升高？
A2：Roopwani(2019) 实测 R 分别为 77.7%±1.2、85.8%±2.5、93.5%±0.8。共处理在颗粒层面改善表面/塑性、外加润滑在模壁涂覆润滑剂，二者都降低摩擦，使力传递更充分、R 更接近 1。注意两者机制不同（材料属性 vs 过程操作），放大时不可互换。

Q3：弹出力为什么与 R 反着走？
A3：同一份摩擦决定了两者：摩擦越小，R 越高（传递越充分），脱模阻力越小、弹出应力越低。Roopwani(2019) 三档弹出应力 13.9/7.1/1.1 MPa 正好与 R=77.7%/85.8%/93.5% 反序对应，是"降摩擦双赢"的直接证据。

Q4：硬脂酸镁加多少？加多了有什么坏处？
A4：底账记录常用0.5–1.5% w/w（★ 行业站，以官方文档为准）。它通过降低模壁摩擦提升 R、降低弹出力，但润滑剂（尤其疏水硬脂酸镁）加多会降低片剂强度、延缓崩解，因此不是越多越好，需与片强/崩解综合权衡。

Q5：怎么把上下冲力、摩擦功、弹出力、压制功写成能量自洽的模块？
A5：用上冲力-位移曲线积分得压制功W=∫F dy（trapezoid，参数顺序 F、y）；由F_lower=R·F_upper把上下冲力差积分得摩擦功；由F_eject=σ_eject·(πDt)（侧面积近似、以官方/原文为准）算弹出力；并用R(c)饱和指数描述润滑效应。整体保证 R<1、摩擦功占比随 R 升高而降低，逻辑自洽。