#弹簧摆的混沌行为.pdf

1 1 2 2 2 ? ? ?? ? ? ? ?? ? ? ? 　第 37 卷 2001 年第 3 期 　　　　　西 北 师 范 大 学 学 报 (自然科学版) 　Vol 37 　2001 　No 3 　　　　　Journal of Northwest Normal University (Natural Science) 　 91 2 2 2 2 · 教学研究 · 弹簧摆的混沌行为 石玉仁 　薛具奎 (西北师范大学 　物理与电子工程学院 　甘肃 兰州 　730070) 摘 　要 : 对受迫阻尼弹簧摆系统的运动进行了分析 , 结合数值计算和理论分析判断了非线性方程的解 的各种性质 , 即系统运动的各种状态. 结果表明 , 在系统参数的某个范围内 , 弹簧摆的运动状态呈现 十分丰富的分岔 、混沌等现象. 关键词 : 弹簧摆 ; 分岔 ; 混沌 中图分类号 : O 322 　　　文献标识码 : A 　　　文章编号 : 1001 988 Ⅹ(2001) 03 0091 07 近 10 年迅速发展起来的混沌理论 , 为解决复杂的非线性问题带来了希望 , 也为经典的 牛顿力学注入了新的活力 , 将牛顿力学提高到了一个新的认识水平. 因此 , 改革传统的物理 学内容 , 将混沌理论引入大学物理学已势在必行. 对此 , 文献[ 1 ]作了很有意义的工作. 新 近出版的一些力学教材[2 ,3 ] 也作了很有益的尝试. 笔者对 1 个典型的非线性系统 ———弹簧振 子 , 作了较全面的理论和数值分析 , 展示了弹簧振子系统的复杂特征. 1 　模型及其动力学方程 在轻质弹簧下面挂一小球 , 小球质量为 m , 可看作质点 , 弹 簧的质量忽略不计 , 便组成一个弹簧摆系统. 设弹簧原长为 l0 , 倔强系数为 k , 在所研究的范围内遵从虎克定律 , 并且假定弹簧摆 在运动过程中不会弯曲. 在极坐标系 ( r ,θ) 下 , 小球运动过程中的受力情况如图 1 所示. 其中 k ( r - l0 ) 为弹性恢复力 ; γr 、γθr 分别为小球在 r 方向和θ方 向受到的空气阻力 , γ为阻尼系数 ; mg 为重力 ; F = AcosωD t 为 图 1 　弹簧摆受力分析 弹簧摆系统受到的一水平方向周期驱动力 , ωD 为驱动力圆频率 , A 为驱动力振幅. 则得小 球的运动方程为 m ( θr¨ + 2 θr ) = - γθr - mgsinθ + AcosωD tcosθ, 2 ( 1 ) m ( ¨r - θr ) = - k ( r - l0 ) - γr + mgcosθ + AcosωD tsinθ. 2 2 令 ω0 = g/ l0 , 两边同除以 m