单摆初始摆角对单摆运动的影响.doc

单摆初始摆角对单摆运动的影响翁丰壕(南昌航空大学飞行器学院330063)摘要:在大摆角情况下,单摆运动方程sinθ≈θ中假设不成立,所以,近似后得到的解在大摆角的情况下也不准确。单摆的运动方程是二阶常微分方程,本文利用MATLAB中ode函数直接对其求解,得到各初始摆角情况下的运动图像。分析其运动规律,表明摆角在某个临界值附近时,单摆的运动情况会有很大差别,即混沌现象。关键词:摆角、周期、混沌现象单摆模型是一个形状、大小都可以看成质点的小球系在不计伸长和质量的摆线上的理想模型。它是处理摆动问题中必不可少的模型,作为一种经典物理模型,已在大多大学物理教材中有涉及。在处理单摆模型时只讨论了小角度摆动时的近似运动状况,当单摆做大摆角运动时,运动微分方程为非线性方程,由于计算上的困难没有涉及。,其求解的精确度在全局下未必能够得到保证。学过常微分方程我们知道,对于高阶的常微分方程,可以降价为一阶微分方程组,进而可以使用数值分析中的相关理论得到较为精确的数值解。因此可以使用以上的方法对大摆角运动进行精确的求解[3]。本文将在以下几个方面对任意角度单摆问题进行研究:分析相空间轨迹随能量的变化;研究单摆运动周期与最大摆角的关系。无阻尼单摆无驱动力、无阻尼的单摆的运动可由的方程(1)描述,其中θ为铅直方向与摆线的夹角,l为摆长,g为重力加速度。θ+glsinθ=0(1)数学理论表明,方程(1)没有严格的解析解,一般的情况下,sinθ≈θ的近似不成立。利用数值分析中的龙格库塔法可以求出方程(1)的数值解,参数取值为g=10m/s2,l=1m,m=1kg。单摆的周期性运动的图像如图1所示。可以看出,与小角度摆动不同的是,初始摆角影响单摆的周期,且初始摆角越大,周期越长,当初始摆角为π时,单摆的周期趋于无穷大,即系统出现了相变现象。图1不同初始摆角情况下单摆运动图像文献[1]中给出了无驱动力、无阻尼单摆的周期公式为:T=0π/2dφ1-sin2θ2sin2φ(2)令R=TT0,称做周期比,T0=2πlg,为小角度下周期。则R表示为:R=2π0π/2dφ1-sin2θ2sin2φ(3)画出R随初始摆角θ的关系图如下:图2周期比R与初始摆角θ关系图可以看出随着单摆摆角的增大,它的周期比也会增加,即周期也会增大,而且,增加的速度也在增大,与图1的结果吻合。单摆初始摆角越大,其机械能越大。从能量守恒角度来看,单摆遵守机械能守恒定律:E=12m(lθ)2+mgl(1-cosθ)(4)θ=±2Eml2-2gl1-cosθ(5)画出速度与摆角的关系图(称之为相图[1],[5]),由图3可得,,相轨迹是一条闭合曲线;当能量大于临界能量值时,相轨迹是非闭合的,分为上下两支,按顺时针方向转(角速度小于零)和逆时针方向转动(角速度大于零).当能量E=2mgl时,单摆处于往复摆动与单向转动之间的“临界状态”.该情况下单摆的运动对能量参数十分敏感,当能量略小于E=2mgl时发生往复摆动,当能量略大于E=2mgl时为