单摆非线性动力学.doc

单摆的非线性动力学分析张 亚 兵（兰州交通大学 车辆工程专业，甘肃 兰州，730070）摘 要: 研究单摆的运动，从是否有无阻尼和驱动力方面来分析它们对单摆运动的影响。对于小角度单摆的运动，从单摆的动力学方程入手，借助李雅普诺夫一次近似理论，推导出单摆的运动稳定性情况。再借助绘图工具matlab，对小角度和大角度单摆的运动进行仿真，通过改变参数，如阻尼大小、驱动力大小等绘出单摆运动的不同相图，对相图进行分析比较，从验证单摆运动的稳定性情况。关键词: 单摆；振动； 阻尼； 驱动力 Abstract: The vibration of simple pendulum is studied by 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 and 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; vibration; damp; drive force1 引言 单摆是一种理想的物理模型1，单摆作简谐振动（摆角小于5）时其运动微分方程为线性方程，可以求出其解析解，而当单摆做大幅度摆角运动时，其运动微分方程为非线性方程，我们很难用解析的方法讨论其运动，这个时候可以用MATLAB软件对单摆的运动进行数值求解，并可以模拟不同情况下单摆的运动。随着摆角的减小，摆球的运动速率将越来越大，而加速度将单调下降,至时, 加速度取极小值。本文从动力学的角度详细考察了这一过程中摆球的非线性运，得出了在运动过程中的关系。 图1 单摆模型 2单摆的线性情况2.1线性单摆的无阻尼振动如图所示，忽略细绳重量，也不计小球受到的空气阻力，则上诉单摆可看成理想单摆，对其进行受力分由牛顿第二定律得： （1）因为 （2）把（2）代入（1）式可得 （3）将（3）两端同除以可得 （4）令，其中为自然频率.则（4）可变为 （5）当很小时,故有, （6）解此方程得: （7）若为实数,则有,即 （8）所以, , （9）令,.则有 （10） （11）从能量守恒方面考虑: 可变形为 （12）令，则有 （13）两边同时乘以，得到 （14）在对两边求积分， （15）积分结果为 （16）令(动能),(势能).则有,机械能守恒.为椭圆方程： 9图2 无阻尼单摆的相平面轨迹图3 有阻尼和有驱动力单摆的运动分析有阻尼和有驱动力单摆的运动方程为 (17)在任意大振幅下，方程(17)的解变得十分复杂，下面利用计算机模拟，分别讨论单摆运动随初值的变化和其混沌运动。3.1 初值不同所产生的曲线为简单计，设，当t=0时，两振动初始条件相差极小，有 (18)取，对(6)式在初始值(18)式下利用MATLAB绘图，其变化曲线如图4所示(其中实线为，虚线为)。图3 有阻尼和驱动力的图由图4可以看出，当时，两条曲线重合，两个解不能分辨；但当时，两条曲线不再重合，两个解、完全不一样，这种混沌运动对初始条件的敏感性称为蝴蝶效应。3.2 振幅不同所产生的相图为简单计，设方程(17)中，除驱动参数取变值外，其余参数不变，即。对(17)式在初始值式下利用MATLAB作计算模拟绘图，相图如图3所示。当时，振荡周期等于外加周期力的周期T，应应单周期解，其相图如图4示。图4当时，对应三倍周期解，其相图如图5示。图5当周期强迫力的振幅达到某一临界值时，系统运动出现混沌，其相图如图6示。图64 无阻尼、无驱动力单摆的运动稳定性当单摆的阻尼因数为0时，即当单摆既无驱动力又无阻尼时，单摆的运动方程为此线性系统的本征方程和本征值分别为本征值为纯虚根，线形方程的零解是为稳定的。将初值设定为：初始角度=1.8，角速度为0，此时的阻尼因数驱动力都为0，利用matlab，作出此时的单摆运动的相图为图8所示。图7 阻尼为0驱动力为0单摆的小摆角运动相图由图可知，无阻尼、无驱动力的单摆运动的相图是一个极限环。因此这种单摆的运动是稳定的。由此可以验证，当角度很小时，方程的推导是正确的，即线性时是稳定的。当增大初始角度，令=45时保持其他条件不变，利用matlab绘出大角度单摆无阻尼无驱动里的相图，如图9。图7 阻尼为0驱动力为0单摆的大摆角运动相图由图可知，当摆角变大时，无阻尼、无驱动的单摆运动的相图也是一个极限环，也就是说，此时的单摆运动也是稳定的。综上所述，由方程推导的结果以及利用工具matlab绘图所得到的结论都是相同的，即不论摆角的大小，无阻尼、无驱动力的单摆的运动都是稳定的5 结论本文以单摆为研究对象。研究了实际条件下的单摆问题，即含阻尼的受迫运动。在受迫阻尼运动中，单摆的运动反映出如下特征，即:(l)描述运动特征的动力学方程是非线性的；(2)这些非线性的方程是“确定性的”，不包含任何随时间变化的随机项；(3)在某些情况下，系统运动轨道的时间行为存在对初始条件的敏感性，初始条件的微小差异可能导致结果的变化非常大；(4)整个系统长期行为的全局特征与初始条件无关。从本文的分析中，结合各种情况下的位移时间图像及相图，我们可以得出以下结论：（1）在摆角较小即时，单摆其相轨迹是围绕原点的椭圆曲线，即我们所熟知的简谐运动相轨迹，此时单摆在平衡位置附近作简谐运动；随着摆角的增加，单摆作非线性振动，周期与摆角有关，周期随摆角的增大而增大。因此，当摆角较大时，单摆运动存在若干稳定点和鞍点(不稳定点)，振动曲线和周期是非线性振动的结果；（2）在有阻尼和有驱动力的情况下，非线性单摆的振动对初始条件非常敏感，称为蝴蝶效应； 在有阻尼和有驱动力的情况下，若取某一参数变值由小到大，其余参数不变，非线性振动的相图会出现由单周期解倍周期解四周期解混沌单周期解，如此反复。参考文献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. 湖南理工学院学报，第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 谢柏松.单摆运动的同宿轨道分叉、次谐