橡胶气球的非单调压强–体积曲线与双气球失稳：以循环充放气量化 Mullins 预处理对失稳阈值的影响

标签：优先级:中高 · 桌面实验/软物质-热力学族 · 预算档位:≤400元 · Arduino压强传感器 · 技能取向:细致标定+非线性物理

1 · 研究问题

乳胶气球的压强–体积曲线 P(V) 的非单调形状（先峰后谷再上翘）经过多少次充放气循环（Mullins 预处理，Mullins effect）后收敛到稳定曲线——逐循环峰值压强的衰减是否呈指数收敛，特征循环数是多少？进一步：预处理程度不同的两只气球连通（双气球实验，two-balloon experiment）时，终态（等大/不等大）边界是否随预处理状态按 P(V) 曲线预言移动？

2 · 研究背景与空白

气球膜是超弹性（hyperelastic）橡胶球壳，Laplace 关系叠加应变硬化给出非单调 P(V)：这是双气球"大者更大"反直觉失稳的根源。Merritt & Weinhaus（Am. J. Phys. 46, 1978）给出双气球实验经典分析；Müller & Strehlow 系统测过气球压强曲线；已发表实验明确指出乳胶有显著迟滞——"气球须充放至少 10 次才能得到稳定曲线"，且隔天部分回复（Mullins 恢复）。即：曲线形状、双气球失稳、迟滞的存在均已发表。

空白在：文献把"预处理 ≥ 10 次"当经验操作规程，未见对收敛动力学（逐循环峰压衰减律、特征循环数、隔天回复量）的定量测量，也未见把预处理状态当自变量去移动双气球失稳边界的实验。这是"评价维度"型空白：把别人扔进 SOP 的步骤变成被测对象。器材只需压强传感器与注射器，正负结论均有效。

3 · 可检验假设

• H1：逐循环峰值压强 P_peak(n) 相对首循环的衰减满足指数收敛 P_peak(n)-P_∞∝ e^-n/n，特征循环数 n 的 95% CI 上界 ≤ 10（与经验规程一致）；静置 24 h 后回复量可测且 > 3 倍测量误差。
• H2：一只"新"气球与一只充分预处理气球连通时，气体净流向由两者当前 P(V) 曲线预测的压差方向决定，预测正确率 ≥ 90%（≥ 20 组）。

4 · 量化验收标准

1. 方法学校验（硬门槛）：MPX5010DP 用自制水柱压强计在 0.3–6 kPa 五点标定，偏差 ≤ 2%；体积测量（大注射器分度 + 排水法交叉）互差 ≤ 3%。同批气球在充分预处理后复现文献 P(V) 非单调形状（峰–谷结构清晰，峰压量级 3–6 kPa）。不过关则全部数据无效。
2. 充放气用 150 mL 注射器定容步进（每步 50 mL，静置 10 s 读稳态压），加载/卸载全程记录，迟滞环面积逐循环报告；每种气球型号 ≥ 5 只样本，报样本间 RSD。
3. 拟合：指数收敛用非线性最小二乘 + bootstrap 95% CI，并用 Theil–Sen 对 ln(P_peak-P_∞)–n 做稳健交叉验证；P_∞ 作为拟合参数并报告不确定度。
4. 误差源量化：温度（橡胶模量温敏，实验室 ±2 ℃ 内有效并记录）、气体温升（缓慢步进准等温，步进后 10 s 弛豫）、漏气（每循环终点闭阀 60 s 压降 < 1% 方有效）、气球批间差（同袋同色）。
5. "隔天回复不显著"亦为有效结论，须给出可分辨 3% 峰压变化的误差棒证明。
6. 全部 P–V 原始数据与脚本开源。

5 · 数据与工具（含器材清单与估价）

用途	来源 / 工具	估价
MPX5010DP ×2 + Arduino Nano + 三通/阀门/软管	淘宝	≈90 元
150 mL 注射器 ×2、透明水柱管（标定用）	淘宝	≈50 元
同批次乳胶气球 100 只（两种尺寸）	淘宝	≈60 元
排水法量筒/量杯、温度计	淘宝	≈60 元
手机延时 + 白板记录；Python 拟合（纯 CPU）	免费	0 元
对照基准	Merritt & Weinhaus 1978 双气球理论、已发表 P(V) 曲线与迟滞描述（仅校验对比，不计入贡献）	—

合计约 300–400 元。安全：气球爆破声与碎片，充气不超过标称直径 90%，戴护目镜。

6 · 方法路径

1. 搭气路，完成压强/体积标定与 P(V) 形状复现（第 4.1 条）。
2. 制定循环 SOP（步进速率、静置时间、终点体积固定），预实验确认漏气判据可执行。
3. 对 ≥ 5 只气球逐循环测 P(V) 至 20 循环，拟合收敛律与 n^*（H1 前半）。
4. 静置 24 h/72 h 复测，量化 Mullins 回复（H1 后半）。
5. 配对不同预处理程度的双气球连通实验 ≥ 20 组，逐组先用各自 P(V) 曲线做书面预测再开阀（预测判据预先写死），统计正确率（H2）。
6. 独立交叉校验：用视频测气球直径换算体积（球形近似，量化非球形误差），与注射器计量比对。
7. 视频：P(V) 采集全程延时、双气球开阀对决实拍（预测卡片同框）、装置讲解；成片 ≤ 3 min。

7 · 新颖性边界

本课题不声称发现非单调 P(V)、双气球失稳或 Mullins 效应——均已发表（Merritt & Weinhaus 1978；Müller & Strehlow；已发表实验明确记载"≥10 次循环 + 隔天部分回复"）。丘奖 2022 铜奖 "Exploring the Thermodynamic Properties of Rubber Band and Related Refrigerator Design"（上海平和）研究橡皮筋熵弹性制冷，与本题同属橡胶弹性但对象、观测量不同，须点名。本项目主结论是：把预处理从操作规程变为被测对象——收敛动力学的首个定量描述（n、回复量），以及预处理状态对双气球终态的可预测操纵。若 n 远大于 10，则直接修正流行实验规程；若双气球终态不可由 P(V) 预测，则说明速率/迟滞效应不可忽略——两种结果都有效。

8 · 决策门槛（go / no-go）

• 第 4 周末：压强标定或漏气判据不过 → 更换阀门与硅胶管；2 周不过则改用凡士林密封 + 缩短单循环时长的快测方案（准等温性下降须在误差预算中加项）。
• 第 9 周末：样本间 RSD > 15%（批间差淹没循环效应）→ 降级路径 A：全部结论改为单球内配对（每球自身第 n 循环 vs 第 1 循环），放弃跨球平均，主结论框架不变。
• 第 18 周末：双气球开阀后流动过快无法记录中间态 → 在连通管加毛细节流，把 H2 降级为仅判终态方向；判据不变。
• 预算裁剪顺序：先砍第二种气球尺寸，再砍第二只传感器（双气球压差改为交替单点测量，时间分辨率降低）。

---