电涡流阻尼下落的终端速度测量：以壁厚–温度二维扫描检验薄壁标度并反演铜的电阻率温度系数

标签：优先级:最高 · 桌面实验/电磁族 · 预算档位:≤700元 · Arduino+视频测量 · 技能取向:动手+数据拟合

1 · 研究问题

钕磁铁在非铁金属管中下落的终端速度（terminal velocity）v_t，是否在薄壁极限下严格满足 v_t ∝ ρ_e/(w)（ρ_e 为管材电阻率、w 为壁厚），该标度在壁厚增大到何处失效？进一步：能否把该装置当作"无电极电阻率计"，仅凭 v_t(T) 反演出铜的电阻率温度系数（temperature coefficient of resistivity）α，并与文献值 0.0039 K⁻¹ 在误差棒内一致？

2 · 研究背景与空白

磁铁穿过导体管时，管壁感生电涡流（eddy current），其磁场对磁铁产生阻尼力，磁铁在几厘米内即达终端速度，典型量级 0.03–0.2 m/s，肉眼与手机 60 fps 视频即可分辨。薄壁偶极子模型给出 v_t = C  m g  ρ_e  a⁴/(μ₀² m_d² w)（a 为管半径、m_d 为磁偶极矩），是少数中学装置能逼近解析理论的电磁问题。

已有工作：MacLatchy 等（Am. J. Phys., 1993）建立经典测量方案；Donoso 等（Eur. J. Phys., 2009）给出薄壁解析模型与实验对照；2021 年有论文报告了铜管温度对阻尼力的影响（终端速度随管温升高而增大），并指出此前文献只扫过厚度与材质、未扫温度。COMSOL 官方也有此算例。可见单变量版本已被反复做过。

空白在：没有公开工作在同一装置上做壁厚 × 温度二维扫描并用二维数据联合反演 α、同时标定薄壁标度的失效边界（厚壁下 v_t 对 w 的依赖偏离 1/w）。该空白属"参数维度 + 评价维度"型，方法完全公开、结论正负都成立，适合一年期学生课题。

3 · 可检验假设

• H1：对壁厚 w/a ≤ 0.15 的管，v_t 与 1/w 的双对数斜率为 -1.00 ± 0.10（95% CI）；当 w/a ≳ 0.3 时斜率显著偏离 −1（CI 不含 −1）。
• H2：在 0–80 ℃ 范围内由 v_t(T) 反演的铜电阻率温度系数 α 与文献值 0.0039 K⁻¹ 之差不超过 15%，且误差棒足以分辨 0 与 0.0039（即能否证"无温度效应"）。

4 · 量化验收标准

1. 方法学校验（硬门槛）：用同一视频测量链对无磁性等尺寸铜柱做自由落体，拟合 g 与 9.79 m/s² 偏差 ≤ 2%；再用单摆重测 g 交叉验证，偏差 ≤ 1%。不过关则后续全部数据无效。
2. 终端速度判据预先写死：取管中段 60% 区间，位置–时间线性拟合残差的 RSD ≤ 2% 方判定"已达终端"；每个 (w,T) 参数点重复下落 ≥ 10 次，报均值 ± 标准误。
3. 标度指数拟合用正交距离回归（orthogonal distance regression, ODR）与 Theil–Sen 稳健回归交叉验证，自助抽样（bootstrap，≥2000 次）给 95% CI；不得只报双对数最小二乘。
4. 主要误差源清单必须量化：释放初态（用非磁性导管定心）、管温不均（管外保温棉 + 上下两只 DS18B20，温差 > 2 ℃ 的轮次作废）、磁矩随温度漂移（磁铁不入水浴，只加热管子，并用手机磁力计核查磁铁温度前后 B 变化 < 1%）。
5. "差异不显著也是有效结论"条款：若某两档壁厚的 v_t 差异不显著，必须给出功效说明——误差棒需小到能分辨按模型预言的差值的 1/2。
6. 全部脚本（Tracker 导出 + Python 拟合）与原始视频开源，第三方可一键重跑。

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

用途	来源 / 工具	估价
铜管 4 根（同内径 ≈16 mm，壁厚 1/1.5/2/3 mm，长 0.5 m）+ 铝管 2 根	五金市场/淘宝定制切割	≈350 元
NdFeB 磁铁（D12–15 mm 圆片，叠 2–3 片）	淘宝	≈60 元
Arduino Nano + 2×DS18B20 温度探头 + 光电门（对射式，备选计时方案）	淘宝	≈80 元
保温水浴（家用保温桶 + 冰/热水）	家用	0 元
视频测量	手机 60/240 fps + Tracker（免费开源）/OpenCV	0 元
拟合与统计	Python（scipy.odr、numpy bootstrap），纯 CPU 秒级	0 元
对照基准	Donoso 2009 薄壁公式、铜 α 文献值（仅用于校验对比，不计入本项目贡献）	—

合计约 500–700 元。安全：磁铁夹伤与低温/热水烫伤为仅有风险，操作时戴棉手套、热水 ≤ 80 ℃。

6 · 方法路径

1. 搭装置 + 完成第 4.1 条方法学校验（自由落体 g、单摆 g）。
2. 标定终端速度测量链：同一条件连拍 15 次，确认 RSD ≤ 3%，否则先改进释放机构。
3. 室温下扫壁厚（4 档铜 + 2 档铝），ODR 拟合双对数斜率与失效边界。
4. 定一根薄壁铜管扫温度（0/20/40/60/80 ℃ 五点），反演 α。
5. 用铝管重复温度扫描做独立交叉校验（铝 α ≈ 0.0043 K⁻¹）。
6. 汇总二维数据，与薄壁/厚壁模型全局对照，写清残差结构。
7. 录制答辩视频：装置全景、带刻度尺的下落慢动作、温控读数同框，保留全部原始片段。

7 · 新颖性边界

本课题不声称发现新电磁物理，也不提出新模型；v_t 的薄壁公式、温度效应的存在均已发表（MacLatchy 1993；Donoso 2009；2021 年温度论文）。已有工作各自只扫单一变量。本项目的主结论是：(i) 同一装置上壁厚标度的定量失效边界；(ii) 由 v_t(T) 反演 α 的精度评估——把演示实验升级为"无电极电阻率测量法"并给出其不确定度预算。若反演值与文献一致，是该方法可用性的独立支持；若系统偏离，偏离量本身即厚壁/端部效应的量度。两种结果都成立。

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

• 第 3 周末：方法学校验未达标（g 偏差 > 2%）→ 停止扫描，先排查帧率标定与透视畸变；连续 2 周不过则改用光电门计时替代视频（备选方案已购）。
• 第 8 周末：单条件重复性 RSD > 5% → 降级路径 A：放弃温度轴，专做室温壁厚 + 材质（铜/铝/黄铜）标度与失效边界，主结论框架（标度检验 + 不确定度预算）不变。
• 第 16 周末：若温度扫描中磁矩漂移无法控制在 1% 内 → 降级路径 B：改为"冰水 vs 室温"两点式差分测量，只检验 H2 的方向与量级，不报 α 的精确值。
• 预算裁剪顺序：先砍铝管（−100 元），再砍 3 mm 厚壁管；两者砍掉后 H1 的失效边界结论弱化为"未观察到失效"，其余结论保留。

---