简介:这份大学物理实验报告文档面向高校理工科学生及物理实验初学者,围绕“测量刚体的转动惯量”这一经典力学实验,提供从实验目的、原理推导到数据处理与结论的完整记录。报告以塔轮、伸杆及配重物组成的刚体为研究对象,通过砝码下落时间间接求解转动惯量,并验证刚体转动定律,同时涵盖作m–1/t²图与r–1/t图两种曲线改直的数据处理方法,以及质量分布对转动惯量影响的定性观察。资源包共1个doc文件,大小约1.73MB,内容结构清晰,包含实验仪器、操作步骤、数据记录与处理结果,可直接作为实验预习、报告撰写或复习参考。目前已有637人学习下载,适合需要规范实验流程、理解转动定律与作图法处理数据的学习者参考借鉴。
1. 从一份 .doc 实验报告说起:转动惯量测量到底在测什么
大学物理实验里,「测量刚体的转动惯量」几乎是每个工科生都绕不开的一课。多数人拿到手的任务,是写一份.doc格式的实验报告:填数据表、算平均值、做误差分析。但真正把这份报告写扎实的人不多,因为大部分人只把它当成「套公式算数」,没搞清转动惯量到底在描述什么。
转动惯量是刚体绕轴转动时惯性大小的量度,记作 (J) 或 (I),单位是 (\text{kg}\cdot\text{m}^2)。它和质量的地位类似:平动里质量越大越难加速,转动里转动惯量越大越难改变角速度。区别在于,转动惯量不只取决于质量,还取决于质量相对转轴的分布——同样一块铁,做成细杆和做成圆盘,绕不同轴的转动惯量能差好几倍。
这份实验报告要解决的核心问题就是:给定一个刚体(常见是三线摆上的圆盘、扭摆上的金属杆,或气垫导轨上的滑块),怎么用可测量的物理量把它的转动惯量算出来,并评估结果可信度。适合正在做这个实验的学生,也适合需要复现实验数据处理流程、把纸质报告电子化的助教和工程师。
2. 转动惯量的理论模型与三线摆测量原理
2.1 转动惯量的定义式与平行轴定理
转动惯量的离散定义式是:
[ J = \sum_i m_i r_i^2 ]
连续体则写成积分:
[ J = \int r^2 , dm ]
其中 (r) 是质量元 (dm) 到转轴的垂直距离。这个式子说明,离轴越远的质量贡献越大,而且是按平方增长。所以设计飞轮时把质量往边缘堆,就是为了在总质量不变的前提下把转动惯量做大。
当刚体绕不同轴转动时,用平行轴定理换算:
[ J = J_C + m d^2 ]
(J_C) 是绕质心轴的转动惯量,(m) 是总质量,(d) 是两轴间距。实验里测出的往往是绕某个非质心轴的转动惯量,报告里要换算回质心轴,这一步经常被漏掉,导致结果系统性偏大。
2.2 三线摆测转动惯量的力学推导
三线摆由上下两个圆盘和三条对称悬线组成,下盘绕中心竖直轴做小角度扭转振动。设下盘质量 (m_0),悬线长 (L),上下盘悬点所在圆半径分别为 (r)、(R),扭转周期为 (T),则下盘绕中心轴的转动惯量为:
[ J_0 = \frac{m_0 g R r}{4 \pi^2 L} T^2 ]
这个公式的推导思路是:扭转时下盘重心升高,重力势能转化为动能,在小角度近似下把运动方程线性化,得到简谐振动形式,周期与转动惯量的平方根成正比。所以只要测准 (T),就能反推 (J_0)。
放上待测物体后,系统总转动惯量变为 (J_0 + J_x),周期变为 (T_1),于是:
[ J_x = \frac{(m_0 + m_x) g R r}{4 \pi^2 L} T_1^2 - J_0 ]
注意这里质量项要用总质量,因为待测物体也参与了重力势能变化。很多报告在这里直接减 (J_0) 而忘了改质量,误差能到百分之几。
2.3 扭摆法与气垫导轨法的适用边界
三线摆适合测圆盘、圆环这类轴对称刚体,操作简单但对悬线对称性要求高。扭摆法用一根金属丝提供恢复力矩,公式为:
[ J = \frac{k}{4 \pi^2} T^2 ]
(k) 是金属丝的扭转系数,需要先用已知转动惯量的标准件标定。扭摆法适合测形状不规则的小刚体,但标定环节引入的误差会传递到结果里。
气垫导轨法把滑块放在导轨上,通过弹簧或滑轮系统施加力矩,测角加速度反推转动惯量。它适合验证平行轴定理,因为滑块位置可以精确移动,(d) 能直接读出。缺点是导轨水平度、气源稳定性都会影响结果。
| 方法 | 适用对象 | 主要误差源 | 典型相对误差 |
|---|---|---|---|
| 三线摆 | 圆盘、圆环 | 悬线不等长、周期计时 | 2%~5% |
| 扭摆 | 不规则小刚体 | 金属丝标定、阻尼 | 3%~8% |
| 气垫导轨 | 滑块、验证平行轴定理 | 导轨水平度、气膜厚度 | 1%~3% |
选哪种方法,取决于待测物体的形状和实验室现有器材。报告里如果只写一种方法,建议在误差分析里说明为什么没选另外两种。
3. 用 Python 处理实验数据并生成报告表格
3.1 从原始计时数据算周期与转动惯量
实验现场通常用光电门或秒表记录多次全振动时间。假设测了 10 个周期共 5 组数据,每组 10T,先算平均周期,再代入公式。下面是一段可直接运行的 Python 代码:
import numpy as np # 原始数据:每组测10个周期的总时间,单位秒 t_10T = np.array([12.34, 12.31, 12.36, 12.33, 12.35]) # 计算单周期 T = t_10T / 10 T_mean = T.mean() T_std = T.std(ddof=1) # 样本标准差 # 三线摆参数 m0 = 0.355 # 下盘质量 kg g = 9.794 # 本地重力加速度 m/s^2 R = 0.0850 # 下盘悬点半径 m r = 0.0425 # 上盘悬点半径 m L = 0.5000 # 悬线长 m # 下盘转动惯量 J0 = (m0 * g * R * r) / (4 * np.pi**2 * L) * T_mean**2 # 不确定度传递:J 与 T^2 成正比,相对不确定度是 T 的两倍 u_J0 = J0 * 2 * (T_std / T_mean) print(f"平均周期 T = {T_mean:.4f} s") print(f"周期标准差 = {T_std:.4f} s") print(f"下盘转动惯量 J0 = {J0:.6f} kg·m^2") print(f"J0 的不确定度 = {u_J0:.6f} kg·m^2")这段代码的逻辑是:先把 10T 换算成单周期,用样本标准差估计计时离散程度,再代入三线摆公式。关键参数说明:ddof=1表示样本标准差,实验数据一般用样本而非总体;g取当地值,不要直接用 9.8,否则会引入约 0.06% 的系统偏差;R、r、L都要用米,量具读数换算时注意游标卡尺和米尺的分度值。
3.2 加装待测物体后的增量法计算
放上待测圆环后,总质量变为 (m_0 + m_x),周期变为 (T_1)。增量法代码:
# 加装待测圆环后的数据 t_10T_loaded = np.array([15.21, 15.18, 15.24, 15.19, 15.22]) T1 = t_10T_loaded / 10 T1_mean = T1.mean() T1_std = T1.std(ddof=1) mx = 0.120 # 待测圆环质量 kg # 系统总转动惯量 J_total = ((m0 + mx) * g * R * r) / (4 * np.pi**2 * L) * T1_mean**2 # 待测圆环转动惯量 Jx = J_total - J0 # 不确定度合成 u_Jx = np.sqrt((J_total * 2 * T1_std / T1_mean)**2 + u_J0**2) print(f"加装后平均周期 T1 = {T1_mean:.4f} s") print(f"系统总转动惯量 = {J_total:.6f} kg·m^2") print(f"待测圆环转动惯量 Jx = {Jx:.6f} kg·m^2") print(f"Jx 的不确定度 = {u_Jx:.6f} kg·m^2")这里用增量法而不是直接套圆环理论公式,是因为实验目的是「测量」而非「验证」。增量法把下盘本身的转动惯量作为已知量扣除,剩下的就是待测物体的贡献。注意质量项必须用总质量,这是前面提到的常见错误点。不确定度按平方和开根合成,因为 (J_x = J_{\text{total}} - J_0),两者独立。
3.3 用 python-docx 把结果写回 .doc 报告
实验报告要求交.doc或.docx,手动填表容易抄错。用python-docx可以把计算结果直接写进模板:
from docx import Document doc = Document("实验报告模板.docx") # 定位到数据表格,假设是第2个表格 table = doc.tables[1] # 填入周期数据 for i, val in enumerate(t_10T): table.cell(i + 1, 1).text = f"{val:.2f}" # 填入计算结果 table.cell(6, 1).text = f"{T_mean:.4f}" table.cell(7, 1).text = f"{J0:.6f}" table.cell(8, 1).text = f"{Jx:.6f}" table.cell(9, 1).text = f"{u_Jx:.6f}" doc.save("实验报告_已填数据.docx")逻辑说明:doc.tables[1]取第二个表格,索引从 0 开始,具体位置要看模板结构;cell(row, col)的行列也从 0 开始,第 0 行通常是表头。参数上,f"{val:.2f}"保留两位小数,与原始计时精度一致;转动惯量保留六位小数是因为数值本身在 (10^{-3}) 量级,保留太少会丢有效数字。
提示:运行前先备份模板,
python-docx会覆盖原文件。如果模板里有合并单元格,cell()的索引可能和视觉行列不一致,建议先用for row in table.rows: print([c.text for c in row.cells])打印结构确认。
4. 误差来源排查与实验报告常见扣分点
4.1 周期计时误差为什么被放大两倍
从公式 (J \propto T^2) 可以看出,周期测偏 1%,转动惯量就偏 2%。这是误差传递里最容易被忽视的一条。假设用秒表手动计时,反应时间约 0.2 s,测 10T 时相对误差是 (0.2/12 \approx 1.7%),传到 (J) 上就是 3.4%。改用光电门可以把计时误差降到毫秒级,相对误差降到 0.1% 以下。
所以报告里如果写「用秒表测量」,误差分析必须体现这个放大效应。具体做法是对 (J = kT^2) 两边取对数再微分:
[ \frac{\Delta J}{J} = 2 \frac{\Delta T}{T} ]
这条式子直接写进报告,比笼统说「计时误差较大」有说服力得多。
4.2 悬线不等长与下盘偏心的影响
三线摆公式假设三条悬线等长且上下盘严格水平。实际装配时,悬线长度差 1 mm,在 500 mm 的线长下相对偏差 0.2%,但会引起下盘倾斜,扭转时出现摆动分量,周期测量会偏大。排查方法是:静止时用水平仪检查下盘,手动扭转后观察是否只有纯转动、没有前后晃动。
下盘偏心则会让质心偏离转轴,等效于给系统附加了一个单摆运动。偏心量 (e) 带来的周期修正约为:
[ T' \approx T \left(1 + \frac{e^2}{2Rr}\right) ]
如果 (e = 2) mm,(Rr \approx 3.6 \times 10^{-3} , \text{m}^2),修正量约 0.06%,看起来不大,但和计时误差叠加后就不可忽略了。报告里可以写「通过调整悬线使下盘居中,偏心量控制在 1 mm 以内」。
4.3 报告里必须交代的不确定度合成步骤
很多报告只给最终不确定度,不写合成过程,这是扣分重灾区。完整的 A 类不确定度(统计)和 B 类不确定度(仪器)要分开算:
# A类不确定度:多次测量的标准差 u_A = T_std / np.sqrt(len(T)) # B类不确定度:仪器误差限除以根号3(均匀分布) delta_inst = 0.01 # 光电门分辨力 0.01 s u_B = delta_inst / np.sqrt(3) # 合成标准不确定度 u_T = np.sqrt(u_A**2 + u_B**2) # 传到转动惯量 u_J = J0 * 2 * u_T / T_mean print(f"u_A = {u_A:.5f} s, u_B = {u_B:.5f} s, u_T = {u_T:.5f} s") print(f"转动惯量合成不确定度 = {u_J:.6f} kg·m^2")参数说明:delta_inst取仪器分辨力或最大允许误差,光电门一般取分辨力;除以 (\sqrt{3}) 是假设误差在区间内均匀分布,这是大学物理实验的通用约定。最后结果写成 (J = (J_0 \pm u_J) , \text{kg}\cdot\text{m}^2),并注明置信概率约 68%。
注意:如果 (u_J / J_0) 超过 5%,先别急着写结论,回头检查是不是漏了质量项或周期单位换算错了。单位错误是实验报告里最高频的低级失误,10T 的时间如果误当成 T 代入,结果会差 100 倍。
5. 把测量结果和理论值对比:验证平行轴定理的进阶做法
5.1 用滑块移动验证平行轴定理
平行轴定理 (J = J_C + md^2) 是转动惯量实验里最适合做线性拟合的验证点。做法是:在气垫导轨上放一个标准滑块,先测质心轴的 (J_C),再把滑块平移距离 (d),测新的转动惯量。理论上 (J) 对 (d^2) 作图应该是一条直线,截距是 (J_C),斜率是 (m)。
import numpy as np from scipy import stats # 滑块平移距离 d 和对应测得的转动惯量 J d = np.array([0.00, 0.05, 0.10, 0.15, 0.20]) # m J = np.array([1.20e-3, 1.35e-3, 1.80e-3, 2.55e-3, 3.60e-3]) # kg·m^2 # 对 d^2 做线性拟合 slope, intercept, r_value, p_value, std_err = stats.linregress(d**2, J) print(f"拟合斜率 = {slope:.4f} kg(理论值应为滑块质量)") print(f"拟合截距 = {intercept:.6f} kg·m^2(理论值应为 J_C)") print(f"相关系数 R = {r_value:.5f}")逻辑说明:linregress返回斜率、截距和相关系数。斜率对应质量 (m),截距对应质心轴转动惯量 (J_C)。如果 (R > 0.999),说明平行轴定理在这个实验精度内成立。参数上,d要用米,且从质心位置开始量;J的单位统一用 (\text{kg}\cdot\text{m}^2),否则斜率的量纲对不上。
5.2 理论值与测量值的偏差怎么解释
以圆环为例,理论转动惯量 (J_{\text{th}} = \frac{1}{2}m(R_1^2 + R_2^2)),其中 (R_1)、(R_2) 是内外半径。把测量值 (J_x) 和理论值对比,算相对偏差:
[ E = \frac{|J_x - J_{\text{th}}|}{J_{\text{th}}} \times 100% ]
如果 (E) 在 5% 以内,通常认为实验合格。超过 5% 时,按以下顺序排查:先看质量称量是否包含夹具,再看半径测量是否用了内径还是外径,最后检查周期是否混入了摆动分量。报告里不要只写「误差在允许范围内」,要具体指出主要误差来源和占比。
| 对比项 | 测量值 | 理论值 | 相对偏差 |
|---|---|---|---|
| 圆环转动惯量 | (J_x) | (\frac{1}{2}m(R_1^2+R_2^2)) | 计算得出 |
| 平行轴定理斜率 | 拟合斜率 | 滑块质量 (m) | 计算得出 |
| 平行轴定理截距 | 拟合截距 | (J_C) | 计算得出 |
5.3 一个容易被忽略的技巧:先标定再测量
不管是三线摆还是扭摆,系统本身的转动惯量 (J_0) 或扭转系数 (k) 都是通过几何参数算出来的,而这些参数(悬线长、悬点半径)的测量误差会直接进入结果。更稳的做法是:用一个已知转动惯量的标准件先标定系统,把 (J_0) 或 (k) 当成拟合参数反推,而不是硬算。
具体操作是:测标准件(比如一个已知质量的均匀圆盘)的周期 (T_s),由 (J_s = \frac{1}{2}m_s R_s^2) 算出理论值,反推系统常数 (C = J_s / T_s^2)。之后测任何物体都用 (J = C T^2),这样悬线长和半径的测量误差被吸收进 (C) 里,不再单独传递。代价是需要一个精度足够的标准件,但实验室通常都有。
这个技巧在报告里可以作为「改进方案」写,比单纯重复教材步骤更能体现对误差传递的理解。最后把标定后的 (C) 和直接计算的 (C) 对比,如果两者相差在 1% 以内,说明几何参数测量可靠;差得多,就要回头检查量具和读数。
本文还有配套的精品资源,点击获取