单摆的混沌运动.docxVIP

广西民族学院学报 (自然科学版) 第 9 卷第 2 期 JOURNAL OF GUANGXI UNIVERSITY FOR NATIONAL ITIES Vol. 9 No. 2 2003 年 5 月 ( Natural Science Edition) May. 2003 文章编号 :1007 - 0311( 2003) 02 - 0021 - 05 单摆的混沌运动Ξ 杨青勇 ( 南宁地区教育学院 物理系 ,广西 南宁 530001) 摘 要 : 采用相图方法和庞加莱截面法描述单摆的复杂运动 ,研究单摆运动中的分岔 、混沌等非线性特 征 1 关键词 : 单摆 ; 相图 ;庞加莱映射 ;混沌 中图分类号 : O313 文献标识码 : A 单摆是物理学中最简单的模型之一 , 传统的力学 教材一般只讨论单摆在摆幅很小的条件下作简谐振 动 、阻尼振动和受迫振动的特征. 事实上 , 如果不限制  0 0 D 0令 ω2 = g/ l ,τ = ω t ,Ω = ω 0 0 D 0 2 F = F , 则 (1) 式变为 : 2  γ 2 mω0  , f = 其摆幅 , 单摆在周期性策动力的作用下 , 其运动将有 m lω0 m g 2θ θ 意想不到的复杂性. 本文将从单摆的动力学方程出 发 ,采用相图 、庞加莱截面等描述方法研究单摆的复 d dτ2 = - 2βd τd τ - sinθ + f cosΩτ (2) 杂运动. 1 单摆模型的动力学方程 我们把传统的单摆模型一般化 :单摆的摆线换成 引入新变量 ω,φ 将 (2) 式化成自治方程形式 : · θ= ω · θ= - 2βω - sinθ + f cosφ φ′= Ω  (3) 质量可忽略不计的刚性杆 , 摆角θ的取值范围不受限 · 制 , 设摆长为 l , 摆球质量为 m , 沿切向受阻力 - γθl ( γ 为阻尼系数) 、重力 的分力 - m gsinθ以及周期策 动力 Fcos ωD t 作用 , 由牛顿第二定律 ∑Fτ = m aτ 得此单摆所满足的动力学方程为 这是一个反映单摆运动所遵循的动力学规律的 不显含时间的微分方程组. (3) 式中有 3 个可调参 量 :β, Ω 和 f , 每个参量的改变都会引起解的变化. 可以通过控制β,Ω, f 参量的变化去考察其解的类型 和结构的变化 , 从而得出反映系统运动特征的信息. ¨ · · m θl = - rθl - m gsinθ+ FcosωDt (1) 为使 (1) 式各物理量无量纲化 , 作如下标度变 换 : 2 单摆运动的相图及庞加莱截面 描述方法[1 ][2 ] Ξ 收稿日期 :2003 - 03 - 16 . 作者简介 :杨青勇(1962 - ) ,男 ,广西马山人 ,南宁地区教育学院物理系讲师. 广西民族学院学报 (自然科学版) 2003 年 5 月 第 9 卷 由于 (3) 式含有非线性项 , 一般而言 , 不能用解 析法求解. 对于这类微分方程 , 法国数学家庞加莱在 十九世纪末创建了一种微分方程的定性理论 , 发明了 相图和拓扑学方法 , 在不求出解的情况下 , 通过直接 考察微分方程的系数及其本身的结构去研究它的解 的性质. 相图方法是非线性动力学最基本的研究方 法. 描写动力学系统状态的空间称为相空间 , 相空间 的每一点都代表系统的一个状态 , 反映系统状态演化 过程的图像就是相图. 下面我们用相图方法以及与之 相联系的庞加莱截面法描述单摆的运动. 2. 1 无阻尼无策动力的单摆的相图 无阻尼 、无策动力相当于 β = 0 , f = 0 的情况 , 此时单摆是一保守系统 , 由 (3) 式得 ·  图 2 无阻尼无策动力的单摆相图 图 2 即为反映了单摆在不同初始值下对应运动 · 状态变化情况的相图. 当 H = 0 ,θ = 0 时 ,θ= 0 , 对 应于相图中的不动点 0 , 为单摆的稳定平衡点 ; 当 H ω′= θ dτ = - sinθ (4) 较小时 , 轨线为圆形闭合回线 , 表明单摆作简谐振动 ; 当 H 增大时 , 轨线将偏离圆形 , 系统运动表现出非线 · · · ·把dθ = dθ dθ = θdθ代入 (4) 式 , 整理得 · dτ dθ dτ dθ · · θdθ= - sinθdθ (5) (1) 若 θ很小 , sinθ≈θ, 在此情况下我们容易得 到 (4) 、( 5) 的解 , 它归结为简谐振动情形. 为便于从 相图的角度来分析其运动特征 ,