简介这份大学物理实验报告文档面向高校理工科学生及物理实验初学者围绕“测量刚体的转动惯量”这一经典力学实验提供从实验目的、原理推导到数据处理与结论的完整记录。内容涵盖刚体转动定律MIβ的验证思路以及作m–1/t²图、r–1/t图两种曲线改直法求转动惯量的具体方法并涉及塔轮、伸杆、配重物组成的刚体在砝码拖动下的运动分析还包含质量分布对转动惯量影响的定性观察。资源包共1个doc文件大小约1.73MB结构完整、公式与步骤清晰可直接用于实验预习、报告撰写参考或数据处理对照。目前已有637人学习下载适合需要快速掌握实验原理、理清作图法处理数据思路、完成规范实验报告的学生参考使用。1. 从一份 .doc 实验报告说起转动惯量测量到底在测什么很多人拿到《大学物理实验报告测量刚体的转动惯量.doc》这类文档第一反应是抄数据交作业。但如果你真去实验室转一圈会发现这套装置测的是转动惯量 I验证的却是 M Iβ 这条转动定律。塔轮、伸杆、配重物组成一个绕竖直轴转动的刚体砝码下落拖动细线细线绕在塔轮上整个系统的运动被简化成砝码匀加速下落 刚体匀角加速转动两个耦合过程。文档里给出的核心近似式 mgr 2hI/(rt²) 就是全部实验的落脚点r、h、t、m 四个量里前三个用卡尺、钢尺、秒表直接量m 是砝码质量I 就是待求量。适合刚接触转动定律的本科生也适合想用最小二乘或 Python 复现数据处理的人。真正难的不是公式而是把曲线改直这件事做对——把 m 对 1/t² 作图、把 r 对 1/t 作图两条直线的斜率分别对应两个不同的 I 表达式斜率取错、单位漏换结果能差一个数量级。2. 转动定律的数学推导与两种曲线改直方案2.1 从牛顿第二定律到转动定律的联立砝码受重力 mg 和绳张力 T下落加速度 a运动方程是 mg − T ma。刚体受绳的力矩 Tr 和轴摩擦力矩 Mf转动方程是 Tr − Mf Iβ。绳与塔轮无相对滑动所以 a rβ。砝码从静止下落高度 h有 h at²/2。四个方程联立消去 T、a、β得到m(g − a)r − Mf 2hI/(rt²)当 m 远小于刚体质量、a 远小于 g 时g − a ≈ gMf 相对 Tr 可忽略于是近似为 mgr 2hI/(rt²)。这一步近似是整个实验的误差来源之一文档里砝码质量不宜太大的注意事项就是为它服务的。2.2 方案 A固定 r 和 h作 m–1/t² 图固定塔轮半径 r 和下落高度 h把近似式改写为 m K₁·(1/t²)其中 K₁ 2hI/(gr²)。实验中换不同砝码质量 m测对应下落时间 t作 m 对 1/t² 的图。如果图是直线转动定律就被验证了。斜率 K₁ 拿到后反解 I K₁gr²/(2h)。2.3 方案 B固定 m 和 h作 r–1/t 图固定砝码质量 m 和下落高度 h把近似式改写为 r K₂·(1/t)其中 K₂ √(2hI/(mg))。换不同塔轮半径 r测同一砝码的下落时间 t作 r 对 1/t 的图。斜率 K₂ 拿到后反解 I K₂²mg/(2h)。两种方案测的是同一个刚体理论上 I 应该一致。文档给出的结果一个是 1.78×10⁻³ kg·m² 量级另一个是 1.087×10⁻² kg·m² 量级差了近六倍——这不是实验做错了而是两种方案里刚体的配重位置、塔轮半径选取不同加上近似条件满足程度不一样。下面这张表把两种方案的差异列清楚对比项方案 Am–1/t²方案 Br–1/t固定量r、hm、h自变量砝码质量 m塔轮半径 r因变量1/t²1/t斜率K₁ 2hI/(gr²)K₂ √(2hI/(mg))反解 II K₁gr²/(2h)I K₂²mg/(2h)主要误差源砝码质量小、a≈g 近似换半径时滑轮高度未同步调提示方案 B 换塔轮半径时必须重新调滑轮高度保证拉线与塔轮平面共面否则绳的张力方向偏离切向力矩 Tr 直接算错。3. 用 Python 复现数据处理斜率拟合与 I 的计算3.1 把实验数据整理成可计算的数组假设你从实验记录里拿到方案 A 的一组数据r 3.00 cm 0.0300 mh 0.800 mg 取 9.80 m/s²砝码质量从 5 g 到 35 g每档测三次取平均。下面用 numpy 做最小二乘拟合直接算出斜率 K₁ 和 I。import numpy as np # 方案 A 数据砝码质量 m(kg) 与三次下落时间(s) m_list np.array([0.005, 0.010, 0.015, 0.020, 0.025, 0.030, 0.035]) t_trials np.array([ [8.42, 8.45, 8.40], [6.01, 5.98, 6.03], [4.92, 4.90, 4.95], [4.28, 4.30, 4.26], [3.84, 3.86, 3.82], [3.52, 3.50, 3.54], [3.27, 3.25, 3.29], ]) t_mean t_trials.mean(axis1) # 每档取平均 x 1.0 / t_mean**2 # 自变量 1/t^2 y m_list # 因变量 m # 最小二乘拟合 y K1 * x强制过原点 K1 np.sum(x * y) / np.sum(x * x) r 0.0300 # 塔轮半径 m h 0.800 # 下落高度 m g 9.80 # 重力加速度 m/s^2 I_A K1 * g * r**2 / (2 * h) print(f斜率 K1 {K1:.6e}) print(f方案 A 转动惯量 I {I_A:.6e} kg·m^2)这段代码的关键点有三个一是t_trials.mean(axis1)对每档三次计时取平均降低随机误差二是x 1.0 / t_mean**2完成曲线改直把非线性关系转成线性拟合三是强制过原点的斜率公式sum(x*y)/sum(x*x)因为理论上 m 0 时 1/t² 0截距应为零。如果拟合出明显非零截距说明系统误差轴摩擦、绳未拉直不可忽略。3.2 方案 B 的拟合与两种结果对比方案 B 换半径测时间自变量是 1/t因变量是 r斜率 K₂ 反解 I 时带平方对时间测量误差更敏感。# 方案 B 数据塔轮半径 r(m) 与同一砝码下落时间(s) r_list np.array([0.015, 0.020, 0.025, 0.030, 0.035, 0.040]) t_B np.array([6.85, 5.42, 4.51, 3.86, 3.38, 3.01]) x_B 1.0 / t_B y_B r_list K2 np.sum(x_B * y_B) / np.sum(x_B * x_B) m_fixed 0.020 # 固定砝码质量 kg I_B K2**2 * m_fixed * g / (2 * h) print(f斜率 K2 {K2:.6e}) print(f方案 B 转动惯量 I {I_B:.6e} kg·m^2)参数说明m_fixed是方案 B 里固定不变的砝码质量必须和实验时实际用的砝码一致h沿用方案 A 的高度前提是两次实验高度确实没变。如果两个 I 值差得离谱先检查单位——r 用米、m 用千克、h 用米任何一处用厘米或克结果都会错 10² 到 10⁶ 倍。3.3 用残差判断数据是否可信拟合完不能只看斜率还要看残差。残差是实测值减拟合值的差如果残差随自变量单调变化说明模型有系统偏差。y_pred K1 * x residual y - y_pred for mi, ri in zip(m_list, residual): print(fm{mi*1000:.0f}g 残差{ri:.6e} kg)正常情况下残差应在零附近随机分布量级远小于砝码质量。如果 5 g 那档残差特别大多半是下落时间太短、秒表启停不同步导致的这也是文档里计时与启动一致这条注意事项的由来。4. 装置调节与误差控制从轴尖接触到秒表启停4.1 转轴铅直与轴尖点接触转动惯量测量的最大系统误差来自轴摩擦力矩 Mf。文档要求轴尖与轴槽尽量为点接触目的是把滑动摩擦变成近似点接触减小 Mf。实际操作时先把刚体放上轴用手轻拨让它自由摆动如果摆动几次就停说明摩擦大能转十几圈才停才算调好。转轴不铅直会让重力产生附加力矩砝码下落时刚体还会进动时间测量直接失真。4.2 拉线缠绕与滑轮共面拉线必须平行缠绕在塔轮上不能重叠。重叠会让有效半径 r 随缠绕层数变化而公式里 r 是常数这一条被破坏整个推导就不成立。滑轮高度要调到拉线与塔轮轴垂直、与滑轮面共面否则绳的张力方向不沿切向力矩变成 Tr·cosθθ 一偏误差就进来了。换塔轮半径时这一步必须重做不能偷懒。4.3 秒表启停与砝码释放的同步计时误差是随机误差的主要来源。常见做法是一手扶住砝码到起始高度另一手拿秒表释放砝码的同时按下秒表砝码触地瞬间再按一次。更稳的办法是用光电门替代秒表但文档里的装置只有秒表那就靠多测几次取平均。每档测三次是最低要求条件允许测五次把标准差算出来超过 0.1 s 的数据直接剔除重测。4.4 砝码质量上限与 a≪g 条件近似式 mgr 2hI/(rt²) 成立的前提是 a ≪ g。砝码越重a 越大近似误差越大。文档建议砝码从 5 g 到 35 g就是为了把 a 控制在 g 的百分之几以内。如果你用 100 g 砝码下落时间可能不到 1 s秒表根本反应不过来而且 a 接近 gg − a 不能再近似成 g公式直接失效。注意实验思考题里问的随机误差和系统误差随机误差主要来自计时和拉线平行度系统误差主要来自轴摩擦和空气阻力。前者靠多次测量平均压制后者只能靠装置调节尽量减小无法完全消除。5. 进阶技巧用光电门和线性回归把误差压到 1% 以内如果实验室条件允许把秒表换成光电门是提升精度最直接的手段。光电门记录砝码通过两个固定位置的时间间隔配合已知距离算出平均速度再外推得到加速度时间分辨率能到毫秒级。具体做法是在砝码下落路径上装两个光电门间距 Δs 固定记录通过时间 Δt则 a 2Δs/Δt² 的修正形式考虑初速度不为零时用 v₂² − v₁² 2aΔs。拿到 a 之后不再依赖 a ≪ g 近似直接用完整式 m(g − a)r − Mf 2hI/(rt²) 反解 I把近似误差也消掉。数据处理上除了强制过原点的拟合还可以做带截距的线性回归用截距反推 Mf。以方案 A 为例拟合 m K₁·(1/t²) b截距 b 对应的是摩擦力矩的等效质量项b 越小说明装置调得越好。用 scipy 的linregress一行就能拿到斜率和截距再算 R²R² 低于 0.995 就说明数据点偏离直线太多得回去查是哪一档测坏了。最后一个容易被忽略的技巧把两种方案测出的 I 值放在一起比较时先确认两次实验的配重物位置是否完全相同。配重物在伸杆上的位置直接改变刚体的质量分布I 对配重位置极其敏感——这正是文档里观察刚体质量分布对转动惯量的影响那一步要定性说明的内容。配重往外移I 增大下落时间变长配重往内收I 减小下落变快。把配重位置、塔轮半径、砝码质量三个变量在实验记录里写清楚后面复现和排错才有依据。本文还有配套的精品资源点击获取