大学物理实验转动惯量测量:三线摆与扭摆原理、Python数据处理及误差分析

📅 发布时间:2026/9/18 18:55:54
大学物理实验转动惯量测量:三线摆与扭摆原理、Python数据处理及误差分析
简介这份大学物理实验报告文档面向高校理工科学生及物理实验初学者围绕“测量刚体的转动惯量”这一经典力学实验展开帮助读者完成从实验原理推导到数据处理与结论撰写的完整流程。文档系统梳理了实验目的、转动定律MIβ的推导、m–1/t²与r–1/t两种曲线改直作图法以及刚体转动仪、滑轮、秒表、砝码等仪器的使用要点并给出质量分布对转动惯量影响的定性观察思路。资源包共1个doc文件大小约1.73MB内容涵盖实验步骤、数据记录表格与处理示例可直接作为实验预习、报告撰写或复习参考。目前已有637人学习下载适合需要规范实验报告格式、掌握作图法求转动惯量的读者借鉴使用。1. 从一份 .doc 实验报告说起转动惯量测量到底在测什么大学物理实验里「测量刚体的转动惯量」几乎是每个工科生都绕不开的一课。多数人拿到手的是一份 Word 文档模板标题写着「大学物理实验报告测量刚体的转动惯量.doc」要求填数据、算结果、写误差分析。但真正让人卡住的往往不是报告格式而是背后那套物理量到底怎么来的转动惯量 J 为什么不能像质量一样直接称出来三线摆和扭摆测的到底是什么公式里的周期 T 为什么要测几十次。转动惯量描述的是刚体绕轴转动时的惯性大小单位是 kg·m²它等于各质元质量与到转轴距离平方的乘积之和。同一个物体绕不同轴的转动惯量完全不同这就是平行轴定理存在的意义。实验的核心任务就是通过测量周期、几何尺寸这些「好测的量」反推出这个「不好测的量」。这份报告适合两类人一类是正在做实验、需要把数据处理和误差分析写清楚的学生另一类是想用代码把重复计算自动化、顺便验证理论公式的工程实践者。下面从原理、装置、公式推导一路讲到用 Python 批量处理数据、排查系统误差把这份 .doc 里该有的东西真正算明白。2. 三线摆与扭摆测转动惯量的原理和公式推导2.1 三线摆的简谐振动模型三线摆由上下两个圆盘和三条等长悬线组成下盘绕中心竖直轴扭转一个小角度后释放做近似简谐的扭转振动。设下盘质量 m₀悬线长 L悬点到盘心距离 r下盘半径 R当扭转角 θ 很小时系统势能变化与 θ² 成正比恢复力矩 M -kθ于是转动周期T 2π√(J / k)通过几何关系把 k 用 m₀g、r、R、L 表示得到下盘空盘转动惯量J₀ (m₀ g R r / (4π² L)) · T₀²放上待测刚体后总转动惯量 J_total 用同样的公式只是质量换成 m₀m、周期换成 T₁待测物转动惯量 J J_total − J₀。这里的关键近似是 θ 足够小一般要求扭转角小于 5°否则 sinθ ≈ θ 的线性化失效周期会随振幅漂移。2.2 扭摆的角谐振动与弹簧扭转常数扭摆是把物体挂在一条金属丝的扭转弹性上绕丝轴扭转后释放。恢复力矩来自金属丝的扭转M -DθD 为扭转常数。周期同样满足 T 2π√(J / D)。要测 J必须先标定 D用一个已知转动惯量的标准体常见是两段对称圆柱或圆环测出周期反解 D再换待测物。装置恢复力矩来源待标定量适用对象三线摆重力势能变化无几何量直接算圆盘、圆环、不规则薄板扭摆金属丝扭转弹性扭转常数 D细长杆、小尺寸刚体落体法砝码重力牵引无带轮轴的规则体2.3 平行轴定理在实验中的用法如果待测刚体绕的轴不通过质心而是平行于过质心的轴、距离为 d则J J_c m d²实验里常把两个相同圆柱对称放在圆盘上改变它们到中心的距离 d测总转动惯量用 J 对 d² 作图斜率就是 2m两个圆柱截距是 2J_c。这是验证平行轴定理最直观的做法也是报告里误差分析的重点d 的测量误差会被平方放大。注意三线摆公式里的 R 和 r 都是「悬点到轴心」的距离不是直径量具用游标卡尺时读数要除以 2这一步错会导致结果差 4 倍。3. 用 Python 处理实验数据从原始周期到转动惯量3.1 数据录入与单位统一实验原始数据通常是「测 30 个周期用时 t重复 5 次」这种形式直接算单周期误差大。先把数据整理成结构化格式统一到国际单位制。import numpy as np # 三线摆几何参数单位统一为 m、kg m0 0.355 # 下盘质量 kg R 0.0625 # 悬点到盘心距离 m r 0.0325 # 悬点到轴心距离 m L 0.500 # 悬线长 m g 9.794 # 本地重力加速度 m/s^2 # 空盘5 次测 30 个周期总时间s t0_30 np.array([42.15, 42.18, 42.12, 42.20, 42.16]) T0 t0_30.mean() / 30 # 单周期 u_T0 t0_30.std(ddof1) / 30 # 周期标准不确定度 # 加待测圆环后 t1_30 np.array([48.90, 48.85, 48.95, 48.88, 48.92]) T1 t1_30.mean() / 30 u_T1 t1_30.std(ddof1) / 30 print(fT0 {T0:.5f} s, u {u_T0:.5f}) print(fT1 {T1:.5f} s, u {u_T1:.5f})逻辑说明用 30 个周期的总时间除以 30 得到单周期比直接测一个周期精度高一个量级。std(ddof1)是样本标准差除以 √n 后得到平均值的标准不确定度这里 n5所以是除以 30 再除以 √5 的等价写法代码里直接对总时间求标准差再除 30对应单周期的不确定度。参数说明ddof1表示贝塞尔校正g取当地值不同城市差异可达 0.03 m/s²报告里应注明取值来源。3.2 转动惯量计算与不确定度传递def J_disk(m, R, r, L, T, g9.794): 三线摆测转动惯量返回 J 及其相对不确定度 J (m * g * R * r) / (4 * np.pi**2 * L) * T**2 return J J0 J_disk(m0, R, r, L, T0) m_ring 0.210 J_total J_disk(m0 m_ring, R, r, L, T1) J_ring J_total - J0 # 相对不确定度主要来自 T平方项和几何量 rel_T 2 * u_T1 / T1 rel_geo np.sqrt((0.0001/R)**2 (0.0001/r)**2 (0.0005/L)**2) u_J J_ring * np.sqrt(rel_T**2 rel_geo**2) print(fJ0 {J0:.6e} kg·m^2) print(fJ_ring {J_ring:.6e} ± {u_J:.2e} kg·m^2)逻辑说明转动惯量正比于 T²所以周期的不确定度贡献要乘 2。几何量用游标卡尺测R、r 的分辨力约 0.02 mm取 0.1 mm 作为 B 类不确定度L 用米尺取 0.5 mm。合成时按相对不确定度平方和开根。参数说明m_ring是待测圆环质量用电子天平测如果圆环内外径已知还可以用理论值 J ½m(R₁²R₂²) 做对照偏差在 5% 以内算正常。3.3 用最小二乘验证平行轴定理把两个圆柱放在不同距离 d 处测总转动惯量做 J–d² 线性拟合。from numpy.polynomial import polynomial as P d np.array([0.00, 0.02, 0.04, 0.06, 0.08]) # 圆柱到中心距离 m J_meas np.array([1.82e-3, 2.05e-3, 2.74e-3, 3.86e-3, 5.41e-3]) coef np.polyfit(d**2, J_meas, 1) # 一次拟合斜率2m截距2J_c slope, intercept coef m_cyl slope / 2 J_c intercept / 2 print(f拟合斜率 {slope:.4f} kg, 反推单圆柱质量 {m_cyl:.4f} kg) print(f截距 {intercept:.4e}, 单圆柱质心转动惯量 {J_c:.4e} kg·m^2)逻辑说明平行轴定理 J J_c md²两个圆柱时总 J 2J_c 2md²所以 J 对 d² 作图为直线斜率是 2m截距是 2J_c。把拟合出的质量和天平称量值对比就能判断实验系统是否可靠。参数说明np.polyfit返回高次到低次系数这里一次拟合返回 [斜率, 截距]。d 的测量要用卡尺测圆柱中心到转轴的距离不是边缘距离。4. 实验操作步骤与常见系统误差排查4.1 三线摆的规范操作流程调节底座水平用水平仪放在下盘上调三个底脚螺丝直到气泡居中。这一步不做扭转时会引入摆动分量周期测不准。测量几何量悬线长 L 用米尺测上下盘悬点间距R 和 r 用游标卡尺测悬点到盘心的距离各测三次取平均。测空盘周期轻轻扭转下盘约 3°5°等振动稳定后开始计时用光电门或秒表测 30 个周期重复 5 次。放待测物把圆环或圆柱对称放在下盘上确认质心与盘心重合重复测周期。记录环境温度、气压对金属丝扭转常数有影响扭摆实验尤其要记。4.2 周期测量误差的三个来源误差来源表现处理办法扭转角过大周期随振幅变化控制 θ 5°从大角度往小角度测看是否漂移下盘晃动周期读数跳动大调水平、等稳定后再计时计时起停偏差人为反应时间约 0.2 s测 30 个周期摊薄或用光电门提示如果 5 次测的 30 周期总时间极差超过 0.3 s先别急着算回去检查水平 and 悬线是否等长数据本身不可靠时任何公式都救不回来。4.3 扭摆标定扭转常数的实操扭摆要先标定 D。取一个已知转动惯量的标准圆环J_std 由几何算得测周期 T_std则 D 4π²J_std / T_std²。再换待测物测 TJ D T² / (4π²)。金属丝的 D 会随使用次数变化同一批实验里标定一次即可但跨天实验要重新标。# 扭摆标定与测量 J_std 1.25e-4 # 标准圆环转动惯量 kg·m^2 T_std 1.842 # 标准圆环周期 s D 4 * np.pi**2 * J_std / T_std**2 T_unknown 2.315 # 待测物周期 s J_unknown D * T_unknown**2 / (4 * np.pi**2) print(f扭转常数 D {D:.4e} N·m/rad) print(f待测物 J {J_unknown:.4e} kg·m^2)逻辑说明D 的单位是 N·m/rad标定后对待测物直接套用周期公式。注意 T 的平方关系周期测偏 1%J 偏 2%。参数说明J_std由圆环质量 and 内外径算出J ½m(R₁²R₂²)标准体的转动惯量精度直接决定整个实验的准确度优先选加工精度高的标准件。5. 把报告写成可复现文档数据表、误差分析与结论写法5.1 数据表该放哪些列一份能复现的实验报告数据表至少包含测量次数、30 周期总时间、单周期、平均值、标准差。几何量单独一张表列出每次测量值和平均值。不要只写最终结果中间量缺失别人无法核对。次数t₃₀/sT/s偏差/s142.151.40500.0007242.181.40600.0017342.121.4040−0.0003442.201.40670.0024542.161.40530.00105.2 误差分析要落到具体项误差分析不是写「存在误差」四个字。要分 A 类统计和 B 类仪器逐项给出数值。比如周期的不确定度来自计时反应时间B 类约 0.2 s 摊到 30 个周期是 0.0067 s和重复性A 类标准差除以 √5。几何量的不确定度来自卡尺分辨力。最后合成相对不确定度和理论值偏差对比判断是系统误差还是随机误差主导。5.3 用理论值做交叉验证圆环绕中心轴的转动惯量理论值 J ½m(R₁²R₂²)。把实验值和理论值比较若偏差超过 10%优先查悬点到轴心距离是否量错、待测物是否偏心、扭转角是否过大。偏差在 5% 以内报告结论就可以写「实验结果与理论值吻合」。# 圆环理论值对照 m_ring 0.210 R1, R2 0.050, 0.070 # 内径、外径 m J_theory 0.5 * m_ring * (R1**2 R2**2) J_exp 3.42e-4 # 前面实验算出的值 print(f理论值 {J_theory:.4e}, 实验值 {J_exp:.4e}, 偏差 {abs(J_exp-J_theory)/J_theory*100:.2f}%)逻辑说明理论值只依赖质量和几何尺寸是独立于周期测量的另一条路径。两条路径结果接近说明周期测量和几何测量没有大的系统偏差。参数说明R₁、R₂ 是圆环内外半径用卡尺测内径和外径后除以 2如果圆环是矩形截面公式要换成对应形式。5.4 报告文档的排版技巧回到那份 .doc 本身公式用 Word 自带公式编辑器或 LaTeX 插件别用图片贴公式否则改数据时对不上。数据表用三线表单位写在表头。不确定度保留一到两位有效数字结果值和不确定度末位对齐。最后附上原始数据照片或电子表格别人拿到这份 .doc 能直接复算才算合格。本文还有配套的精品资源点击获取