大摆角单摆运动周期的三个公式.docVIP

 大摆角单摆运动周期的三个公式 谭志中 (南通师范学院物理系 ,江苏 南通 226007) 要 : 建构出一类非线性动力学方程的近似处理方法 ,对单摆的动力学方程进行了修正 ,得到 单摆运动周期的三个新结论 。 : 非线性 ;单摆 ;周期 : O314 文献标识码 :A 文章编号 :1003 - 7551 (2003) 03 - 0028 - 03 引言 单摆运动问题是一个古老的问题 。众所周知 , 对于摆长为 L , 最大摆角为 θ0 的单摆作无阻尼运动时 , 其受力分析如图 1 所示 。根据转动定理可得到单摆的 动力学方程 : 1 d2θ g d t2 = - L sinθ ( 1 - 1) 式 (1 - 1) 不是一 个 简 单 的 微 分 方 程 , 而 是 一 个 非 线 性 微 分 方 程 , 它 能 产 生“混 沌[1 2 ]”现象 。非线性微分方程的解通常比较复杂 。当单摆摆角比较小 (θ0 50 ) 时 , 各种教材都给出了统一的近似处理方法 , 即对式 (1) 中取近似值 sinθ≈θ得 图 1 d2θ 2 d t2 = - ω0θ ( 1 - 2) g 其中 ω0 = 。显然式 ( 1 - 2) 是一个简单的二阶线性微分方程 , 由式 ( 1 - 2) 可以解得单摆的运动学方程 L ( 简谐运动方程) θ=θ0cos (ω0 t + φ0) (1 - 3) 当然 , 由于式 ( 1 - 2) 是一个近似微分方程 , 因而式 ( 1 - 3) 也是一个近似的简谐运动方程 , 当 θ0 50 时 , 其相应的运动周期公式为 : π T0 = 2 = 2π L (1 - 4) ω0 g 对于摆长为 L , 最大摆角为 θ0 的单摆 , 文 [ 3 ]和文 [ 4 ]分别从机械能守恒定律的角度和相图关系的角 度求出了精确的单摆运动周期公式 : π/ 2 2 T0 dφ π Θ 0 T = (1 - 5) 1 - sin2 (θ 0/ 2) ·sin2φ 式 ( 1 - 5) 适用于任意摆角情形下的单摆运动 , 但这一公式是用椭圆积分表示的 , 过于繁杂 , 应用时需 要查完椭圆积分表 , 实用性不强 。 本文对非线性微分方程 (1 - 1) 进行了一定的研究 , 通过视变量为常量的思维策略 , 得到比较理想的大 摆角情形下的三个简单 、实用 、精度较高的单摆运动周期公式 。 3 收稿日期 :2003 - 07 - 10 211 三个周期公式 对于非线性动力学方程 (1 - 1) 的解 , 本文另通过视变量为常量的方法 , 给出了单摆作无阻尼运动时大 摆角单摆运动周期的一个公式及其两个推论 。我们得到的结论为 T0 T = ( 2 - 1) cos (μθ0) 其中 μ = 0. 497 ,θ0 是指摆线偏离平衡位置的最大偏角 。 根据单摆周期公式 (2 - 1) 可以得到如下推论 T0 (2 - 2) T = 1 - μ1θ2 0 T = T0 (1 + μ1θ2 + μθ4) (2 - 3) 0 2 0 其中 μ1 = 0. 062 ,μ2 = 0. 0044 。 从第三节的相对误差对比表中将会看到式 (2 - 1) 、(2 - 2) 、(2 - 3) 是比较精确的周期公式 。 212 从纯数学的角度推导 sinθ的近似值 θ sinθ= 2sin (θ/ 2) cos (θ/ 2) ≈2 × cos (θ/ 2) =θcos (θ/ 2) ≈θcos (μθ0) ( 2 - 4) 2 此处将变量 cos (θ/ 2) 视为常量 cos (μθ0) , 其中 θ0 是指摆线偏离平衡位置的最大偏角 ,μ 是修正常数 。 将 (2 - 4) 代入 (1 - 1) 得 d2θ 2 d t2 = - ωθ (2 - 5) gcos (μθ0) 其中 ω = = ω μθ cos 0 , 式 ( 2 - 5) 即是对动力学方程 ( 1 - 2) 的修正 , 因为式 ( 2 - 5) 同样是 0 L 一个简单的二阶线性微分方程 , 由式 ( 2 - 5) 很容易解得单摆的周期公式 ( 2 - 1) , 即 T = 2π = 2π T0 L = ω gcos (μθ0) cos (μθ0) 其中 μ是一种修正常数 , 修正的方法是将式 (2 - 1) 与标准式 ( 1 - 5) 在取不同角度时进行比较 。本文 取修正值为 μ = 0. 497 。 213 对推论的推导 证明 : ( 根据单摆周期公式 2 - 1 进行推导) 根据泰勒级数展开式可得 cosλμθ= 1 - λ (μθ) 2 + λ(3λ- 2) (μθ) 4 -