一维振子运动的数值求解研究.docxVIP

?

?

一维振子运动的数值求解研究

?

?

张定梅

摘要：建立了一维振子运动的牛顿运动学方程，利用MATLAB软件数值研究了振子的动力学特性.通过调节参数实现了简谐振动、阻尼振动、受迫振动及复杂非线性振动.给出了各种振动的运动时间序列，讨论了参数对振子的振动周期及振幅的影响，并使用功率谱和吸引子图分析了振子产生的混沌的特性.与传统的解析求解方法相比，数值求解方法更加快捷、直观，避免了复杂的数学积分运算，并且通过运动轨迹图可快速掌握振子的运动规律.

关键词：一维振子;MATLAB;Runge-Kutta算法;混沌

：O322?：A?：1673-260X（2019）09-0019-03

1引言

机械振动是自然界中普遍存在的运动形式，也是大学物理课程中一项重要的教学内容.对于简单的简谐振动可以利用数学积分的方法给出解析解[1]，但对于存在复杂外部扰动的非线性振动系统，解析求解将变得非常困难，因此数值求解将是一种好的替代方式.关于振动问题数值求解的研究已见有关报道[2-4]，例如，文献[2]利用拉格朗日方法建立了复雜的六弹性振子的二维运动方程，并用maple软件进行了数值求解，给出振子运动时间序列图.文献[3]研究了各种组合弹簧振子中的非线性振动问题，给出了数值解和解析解并进行了对比.文献[4]研究了Duffing振子和VanderPol振子耦合之后的非线性的动力学行为，利用Simulink仿真给出了不同频率及耦合系数下耦合振子的相图和庞加莱截面图.在大学物理教学中，一维振子运动是一种最基本的教学内容，为了能够系统的、直观的研究一维振子在各种外部扰动下的动力学特性，本文建立了一维振子的牛顿运动学方程，采用四阶Runge-Kutta算法并利用MATLAB软件进行数值求解，通过合理调节参数，实现了简谐振动、阻尼振动、受迫振动及复杂非线性振动，给出了振子的运动时间序列图，并且得到了有趣的混沌现象.

2四阶Runge-Kutta算法

四阶Runge-Kutta算法是求解常微分方程及方程组常用的数值解法，它基于Taylor展开，要求被求解的函数具有较好的光滑性，并且其截断误差为o（h5），具有较高的求解精度.四阶Runge-Kutta公式的一般表达形式为：

其中：

3一维振子运动学方程

在实际生活中，一维振子除了受到回复力作用以外，往往还受到阻力和外部扰动，为了全面的描述一维振子的运动情况，本文建立的一维振子运动学方程如下[5-7]：

其中，m为振子的质量，-k1x为线性回复力，为与速度成正比的阻力，-k3x3为与位移有关的三阶非线性阻力项，F1cos1t+F2sin2t为周期性的外力.为简化运算，m取1kg，并且将二阶常微分方程化为一阶微分方程组求解，具体表达式如下：

4仿真结果及结果

本文采用四阶Runge-Kutta算法对方程（4）进行数值求解，求解步长为0.01s，采用MATLAB软件进行仿真.

4.1简谐振动

简谐振动是振子最简单的运动形式，此时的振子只受到回复力的作用，方程（4）的参数取值情况如图题所示.图1给出了振子在不同的回复系数下的振动曲线，如图1（a）所示，当A=1s-2时，此时振子的运动轨迹为一规则正弦曲线，振幅为0.110m，振动频率为0.16Hz.当A=3s-2时（图1（b）），由于回复力的增大，振幅变为0.063m，振动频率变为0.28Hz.当A=6s-2时（图1（c）），振幅减小为0.045m，振动频率增大为0.39Hz.从图中可以直观地看出随着回复系数的增大，振子的振动振幅慢慢缩小，同时振动频率慢慢变大.

4.2阻尼振动

阻尼振动是指振子受到外界阻力下的运动，本文中考虑阻力与振子的运动速度成正比的情况.在阻尼振荡下，方程（4）的参数取值情况如图题所示.图2（a）给出了阻尼系数B=0.01s-1时振子振动的时间序列，可以看出，由于此时阻尼系数较小，振子的振幅衰减较慢.当继续增加阻尼系数到0.03s-1，如图2（b）所示，振子振动到200s时振幅几乎为零.继续增大阻尼系数到0.06s-1（图2（c）），振子到150s就几乎停止振动了.因此，从图中可以直观地看出，随着阻尼系数的增大，振子的振幅衰减的越来越快.

4.3受迫振动

受迫振动指的是振子在外部周期力作用下的运动，振子在不受外力作用时，由于自身恢复力的作用会产生一个固有振荡频率，此频率与回复系数有关，本文中回复系数A取2s-2，相应的本征角频率为1.414Hz.图3考虑了三种不同的外部扰动频率下振子的运动情况.图3（a）为D=0.5ms-2，ω1=1Hz情况下的受迫振动振动曲线，此时外部扰动频率与振子的固有频率不相同，但差距较小，振子在72s后趋于规则振荡，振荡幅度为0.52m.当D=0.5ms-2，ω1=1.414Hz时（图3（b