1、1单摆的非线性动力学分析张 亚 兵（兰州交通大学 车辆工程专业，甘肃 兰州，730070）摘 要: 研究单摆的运动，从是否有无阻尼和驱动力方面来分析它们对单摆运动的影响。对于小角度单摆的运动，从单摆的动力学方程入手，借助李雅普诺夫一次近似理论，推导出单摆的运动稳定性情况。再借助绘图工具matlab，对小角度和大角度单摆的运动进行仿真，通过改变参数，如阻尼大小、驱动力大小等绘出单摆运动的不同相图，对相图进行分析比较，从验证单摆运动的稳定性情况。关键词: 单摆；振动； 阻尼； 驱动力 Abstract: The vibration of simple pendulum is studied by
2、analyzing whether or not damp and drive force its influence of the simple pendulum. For small angle pendulum motion, pendulum dynamic equation from the start, with an approximate Lyapunov theory of stability of motion is derived pendulum situation. Drawing tools with help from matlab, small angle an
3、d wide-angle pendulum motion simulation, by changing the parameters, such as damping size, drive size draw simple pendulum of different phase diagram, analysis and comparison of the phase diagram, from the verification the stability of the situation pendulum movement.Key words: simple pendulum; vibr
4、ation; damp; drive force1 引言单摆是一种理想的物理模型 1，单摆作简谐振动（摆角小于 5）时其运动微分方程为线性方程，可以求出其解析解，而当单摆做大幅度摆角运动时，其运动微分方程为非线性方程，我们很难用解析的方法讨论其运动，这个时候可以用 MATLAB 软件对单摆的运动进行数值求解，并可以模拟不同情况下单摆的运动。随着摆角的减小，摆球的运动速率将越来越大，而加速度将单调下降,至时, 加速度取极小值。本文从动力学的角度详细考察了这一过程中摆球的0非线性运，得出了在运动过程中 的关系。.,t图 1 单摆模型 22 单摆的线性情况2.1 线性单摆的无阻尼振动如图所示，忽略细
5、绳重量，也不计小球受到的空气阻力，则上诉单摆可看成理想单摆，对其进行受力分由牛顿第二定律得： sinmga（1）因为 22dtlt )(ls（2）把（2）代入（1）式可得 0sin2mgdtl（3）将（3）两端同除以 可得 ml si2lt（4）令 ，其中 为自然频率.则（4）可变为lg00sin2dt（5）当 很小时, i故有, 02dt（6）解此方程得: titieCet 0021)(（7）若 为实数,则有 ,即*titititi eCe 0000 212*1 （8）所以, , *2C*12C（9）令 , .ieA21ie2则有 )cos()( 0)()( 00 tAt titi（10）)
6、cos()(0tAt（11）从能量守恒方面考虑:3可变形为 02dt 02dt（12）令 ，则有 t 20（13）两边同时乘以 ，得到 d 20d（14）在对两边求积分， d20（15）积分结果为 E2021（16）令 (动能), (势能).21T201V则有 ,机械能守恒.EV为椭圆方程：202图 2 无阻尼单摆的相平面轨迹图3 有阻尼和有驱动力单摆的运动分析有阻尼和有驱动力单摆的运动方程为(17).2sincosft在任意大振幅下，方程(17)的解变得十分复杂，下面利用计算机模拟，分别讨论单摆运动随初值的变化和其混沌运动。3.1 初值不同所产生的 曲线t为简单计，设 ，当t=0时，两振动初
7、始条件相差极小，0.1,2/3f4有 (18).11.2(0), ()0取 ，对(6)式在初始值(18)式下利用MATLAB绘图，其 变化曲线如ts t图4所示(其中实线为 ，虚线为 )。1t2t图3 有阻尼和驱动力的 图t由图4可以看出，当 时，两条曲线重合，两个解 不能分辨；025ts12()t、但当 时，两条曲线不再重合，两个解 、完全不一样，这种混25ts12()t、沌运动对初始条件的敏感性称为蝴蝶效应。3.2 振幅不同所产生的相图.为简单计，设方程(17)中，除驱动参数 取变值外，其余参数不变，即f。0.25,/3对(17)式在初始值式 下利用MATLAB作计算模拟绘图，相图.(0)
8、, ()0如图3所示。.当 时，振荡周期 等于外加周期力的周期T， ，1.06f 2/3T应应单周期解，其相图 如图4示。.5图4当 时， ，对应三倍周期解，其相图 如图 5 示。1.07f3T .图5当周期强迫力的振幅达到某一临界值 时， ，系统运动出1.684f2T现混沌，其相图 如图6示。.图64 无阻尼、无驱动力单摆的运动稳定性当单摆的阻尼因数为0时，即当单摆既无驱动力又无阻尼时，单摆的运动方程为 0sinmgl此线性系统的本征方程和本征值分别为 12i,本征值为纯虚根，线形方程的零解是为稳定的。将初值设定为：初始角度 =1.8，角速度为 0，此时的阻尼因数驱动力都为 0，利用 mat
9、lab，作出此时的单摆运动的相图为图 8 所示。6-0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4-0.4-0.3-0.2-0.100.10.20.30.4度度度度度图7 阻尼为0驱动力为0单摆的小摆角运动相图由图可知，无阻尼、无驱动力的单摆运动的相图是一个极限环。因此这种单摆的运动是稳定的。由此可以验证，当角度很小时，方程的推导是正确的，即线性时是稳定的。当增大初始角度，令=45时保持其他条件不变，利用matlab绘出大角度单摆无阻尼无驱动里的相图，如图9。-0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8-0.8-0.6-0.4-0.200
10、.20.40.60.8度度度度度图7 阻尼为0驱动力为0单摆的大摆角运动相图由图可知，当摆角变大时，无阻尼、无驱动的单摆运动的相图也是一个极限环，也就是说，此时的单摆运动也是稳定的。7综上所述，由方程推导的结果以及利用工具 matlab 绘图所得到的结论都是相同的，即不论摆角的大小，无阻尼、无驱动力的单摆的运动都是稳定的5 结论本文以单摆为研究对象。研究了实际条件下的单摆问题，即含阻尼的受迫运动。在受迫阻尼运动中，单摆的运动反映出如下特征，即:(l)描述运动特征的动力学方程是非线性的；(2)这些非线性的方程是“确定性的” ，不包含任何随时间变化的随机项；(3)在某些情况下，系统运动轨道的时间行
11、为存在对初始条件的敏感性，初始条件的微小差异可能导致结果的变化非常大；(4)整个系统长期行为的全局特征与初始条件无关。从本文的分析中，结合各种情况下的位移时间图像及相图，我们可以得出以下结论：（1）在摆角较小即 时，单摆其相轨迹是围绕原点的椭圆曲线，即我们5所熟知的简谐运动相轨迹，此时单摆在平衡位置附近作简谐运动；随着摆角的增加，单摆作非线性振动，周期与摆角 有关，周期随摆角 的增大而增大。因此，当摆角较大时，单摆运动存在若干稳定点和鞍点(不稳定点)，振动曲线和周期是非线性振动的结果；（2）在有阻尼和有驱动力的情况下，非线性单摆的振动对初始条件非常敏感，称为蝴蝶效应； 在有阻尼和有驱动力的情况
12、下，若取某一参数变值由小到大，其余参数不变，非线性振动的相图会出现由单周期解倍周期解四周期解混沌单周期解，如此反复。参考文献1王海期. 非线性振动M. 高等教育出版社，1992，262-2782刘延柱，陈立群. 非线性振动M. 高等教育出版社，2001，163-1693张德丰. MATLAB 数值计算方法M. 机械工业出版社，2010，225-2394李欣亚，张明路. 机械振动M. 清华大学出版社，2009，587-5895孙中奎，徐伟，杨晓丽. 求解强非线性动力系统响应的一种新方法J. 动力学与控制学报，第 3 卷第 2 期，2005，29-356陈文涛，龚善初. 单摆振动分析J. 湖南理工
13、学院学报，第 21 卷第 1 期，2008，66-707 Kunihik Kaneko.Oscillation and doubling of torus J. Progress of Theoretical Physics,1984,72(2).8 张伟,霍拳忠,李骊.非线性振动系统的异宿轨道分叉、次谐分叉和混沌J.应用数学和力学,1992 ,13:199.9 谢柏松.单摆运动的同宿轨道分叉、次谐分叉和混沌J.北京师范大学学报,2000,36(5):631.10 郎和.保守单摆系统中的混沌运动J.西北师范大学学报 ,2002 ,38 (4):108.11 李元杰.单摆的规则运动及混沌运动的研究J .大学物理,1998 ,17(9):6.8非线性系统控制理论单摆的动力学分析9学院：机电工程学院专业：车辆工程姓名：张 亚 兵学号：0211362
