单摆测量重力加速度试验的误差分析

1、单摆测量重力加速度实验的误差分析吉恒云南省通海县第二中学，云南，玉溪 652701 ）单摆实验是普通物理的基本实验之一 , 同时也是必做实验之一。其原理简 单、易懂，原则上只要在同一地点进行实验，都应得到相同结果，但在实际操作 过程中一些不可避免的因素会影响实验结果的精确度。 为提高实验的精确度， 减 小各种不可避免因素给实验结果带来的影响， 本文从以下几方面着手对此实验进 行分析和研究。首先，对摆角进行分析， 因为随摆角大小的变化， 摆遵循的运动规律是不一 样的。在实验原理中， 一般是把它理想化地当作简谐运动来处理， 让其满足简谐 运动的运动方程， 然后来求解其周期公式， 事实上这是有条件限

2、制的。 因此本文 采用了增维精细积分的方法来讨论单摆在什么样的摆角情况下才能够做线性动 力学分析，也就是单摆满足简谐运动运动规律的摆角范围。其次，单摆摆长的测量也是引起实验误差的原因之一。 本文就单摆摆长的不 同测量方法带来的 B类标准不确定度（由实验仪器的精确度引进）进行计算、分 析、比较，以选取最佳测量方法。 1单摆测量重力加速度的实验原理如图 1所示，单摆就是用一根不可伸长的轻线悬挂一个小球 , 使其可绕摆的 支点 O做摆动, 当小球作摆角很小的摆动时就是一个单摆。设小球的质量为 m , 其质心到支点 o的距离为 l （摆长） 。建立自然坐标系，根据受力分析，作用在小球上的切向力的大小为

3、 mgsin 当 很小时 , 有 sin , 此时切向力的大小近似为 mg 法向，绳的张力和重力的分力相平衡。根据牛顿第二 运动定律 , 质点动力学方程为：ma t mgd2因at l ddt2 ，代入上式得1)d2 dt 2上式即为单摆的运动微分方程。对上式移项得到d2dt2若令则有其解为2)d2dt202 03)A cos 0t式中 A 与是待定常数。由于单摆运动的周期性，应有A cos 0 tA cos 0 t0T又因余弦函数的周期为 2，故 0T 2将式( 2)代入得T24)g42l实验时 , 若只测量一个时间周期，则测量误差相对较大（5）因此一般采用测量连续摆动 n 个周期的时间 t

4、, 此时（ 5）式变为：2 n2lt26）以上为基本实验原理，理想情况（即忽略复摆，空气阻力，空气浮力等因素 对实验的影响） 下只要在同一地点进行实验， 都应得到相同结果， 但在实际操作 过程中一些不可避免的因素会影响实验结果的精度。2. 实验改进摆角范围讨论在上述实验原理中， 谈到当摆角很小时， 可将 sin 近似为 处理，但未确切 说明到底要多小；同时对 角增大了又是怎样的情况也没有考虑进来。因此，若 要建立单摆模型，让单摆摆动起来，摆角范围的大小就是首要解决的问题。记p则单摆非线性动力学方程变为7）gsinl8）根据增维思想 , 通过一维变量 x 1和 x 0的引入, 可将式（ 8）写成

5、矩阵形式的齐次状态方程p00g sinlpx000x0 1 09）T010TV p x T， Vp x T ， H00gsinl000于是式（ 8）简写为V HV根据矩阵分析理论 ,其解为VHte V0（10）其中, V0为初始状态向量, V 为任意时刻的状态向量。为得解的递推表达式 , 将时间轴等分为 N段, 即t0,t1,t2, ,tj,tj 1, tN 时间步长为t t j 1 t j则由t j时刻的状态向量 Vj, 可求tj 1时刻的状态向量 Vj 1。即HtVj 1 e V j（11）Ht其中指数矩阵 e 采用高精度、高效率的精细积分法。以上是采用高精度、高效率的增维精细积分推导出的

6、求解单摆非线性动力摆 角和摆速（取摆角初始速度都为零） 的迭推公式。 通过单摆线性动力响应与非线 性动力响应的对比分析 , 得出了能简化为线性动力问题的最大摆角值及最大摆角 和摆长对非线性动力问题摆角相位、摆速相位的变化规律 , 即当初始摆角小于 0 10o 时, 采用非线性动力方程求解的摆角时程曲线、摆角速度时程曲线分别 与同一问题线性动力方程求得的摆角时程曲线、 摆角速度时程曲线几乎完全重合 达到同一周期两者间的相位差很小 , 可以忽略不计。但是当初始摆角超过 0 10o 后, 随着初始摆角的增大 , 采用非线性动力方程求得的摆角、摆角速度 达到某一平衡位置所需的时间分别比采用线性动力方程

7、所得的摆角、 摆角速度达 到同一周期平衡位置所需的时间越来越长。因此 , 当初始摆角大于 0 10o 后, 不能采用简化的线性动力方程来求解单摆动力响应 , 否则所得结果与实际存在较大的差距 , 而应采用非线性动力方程求解。 所以，若要采用简谐运动方程来求 解单摆周期公式，实验中单摆的摆角应该小于 10o ，误差也才能因此而减小。 换言之，只有在单摆摆动的角度小于 10o 的情况下，单摆的微分方程才可表 示为式（ 1）的形式，才可求解得T 2 l（ 12）g（ ）相应地，重力加速度的表达式为 g 4 2 l2 。此即为实验原理中介绍的周期公式T2和重力加速度公式的来由。摆长的测量方法选择在测量

8、摆长的过程中， 方法选择的不同， 误差的来源以及传递就不同， 引起 的误差大也就不一样。 在选用的长度测量工具都为米尺和游标卡尺的情况下 （即 仪器的精确度等级相同， 其极限误差值相同） ，测量方法对由实验仪器引入的 B类标准不确定度 U （B）的影响是不同的。实验中用毫米刻度的米尺测量摆线长度， 用游标卡尺测量摆球直径， 其中米 尺的极限误差为 = 0.1mm ，游标卡尺的极限误差为 = 0.02mm摆长的测量方法 一般如下三种。如图 2所示， l1 为支点到摆球上端的距离， l2 为支点到摆球下端的距离， D为摆球的直径方法1：l l1 l2 /2其B类标准不确定度为：UB(20.13) (20.13)0.0408mm方法2：l l1 d2图2 单摆摆长B类标准不确定度为：UB0.0954 mm方法3：l l2 d 2UB0.