第三章matlab在力学中的应用

3-3 机械振动物体在平衡位置附近的往返叫做振动或机械振动。

振动的传播称为波，机械振动的传播称为机械波。

振动和波动是涉及物理及众多领域的一种非常普遍而重要的运动形式，研究振动和波动的意义已远远超过了力学的范围。

本节利用MATLAB 来处理机械振动的一些问题。

简谐振动质点在线性回复力作用下围绕平衡位置的运动称为简谐远动，它是最基本的振动。

下面，我们通过两个例子来讨论简谐运动的动力学和运动学特征。

（1） 弹簧振子系统的简谐运动·题目（ex3311）设弹簧阵子系统由质量为m 的滑块和劲度系数为k 的弹簧所组成已知t=0时，m 在A 处，即x 0=A ，并由静止开始释放。

试研究滑块的运动规律。

·解题分析以x 表示质点相对原点的位移，线性回复力f=-kx 。

由牛顿第二定律以及题设条件，可写出弹簧振子的振动微分防尘及初始条件为22t 00(0)(0)0d x kx dt m x A dx v dt =+====滑块速度分别为22dx v dy d x a dt== 令2,kmω=用符号法求解上述微分方程，求出运动方程、速度和加速度，并绘制出,()x t v x a x ---相轨迹和曲线。

（2） 单摆·题目（ex3313）设单摆的摆长为l ，摆锤质量为，将摆锤拉开一角度θ，然后放开使其自由摆动。

在不计空气阻力的情况下，分小摆角和大摆角两种情况，讨论单摆的角位移θ随时间t 的变化规律。

·解题分析由牛顿第二定律，有222sin sin ,d g dt lθθωθω=-=-= 其中，g 为重力加速度。

① 小角摆动假定角位移很小，sin θ≈θ,上式为220d gdt lθθ+= ② 大角摆动222sin sin d g dt lθθωθ=-=-上式是非线性方程。

为了方便起见，将θ用y 来表示，上式又可以写为下列一阶微分方程组1221;sin()dy dy g y y dt dt l==- 用MATLAB 编程解此方程组。

取l=1m,g=s 2。

初始条件取为073ππθ=试取和，比较二者的运动规律。

简谐振动的合成（1）同方向简谐振动的合成 ·题目（ex3321）设一物体同时参加了再同一直线上的两个的简谐振动，其简谐振动分别表示为11112222cos()cos()x A t x A t ωαωα=+=+讨论不同频率、不同初位相时简谐振动的合成。

·解题分析 由题知，合振动为12111222cos()cos()cos()x x x A t A t A t ωαωαωα=+=+++=+其中，A =讨论：① 频率相同，初相位相同例如，取A 1=20,A 2=10，ω1=ω2=5, α1=α2=0,可以看出，合成的振动是振幅为A=A 1+A 2的同频率简谐振动。

② 频率相同，初相位不同例如，取A 1=20,A 2=10，ω1=ω2=5, α1=0,α2=2π/3，或普遍写为Δα=α2-α1=2k π，可以看出，合成的振动是振幅为A=A 1+A 2的同频率简谐振动。

若取A 1=20,A 2=10，ω1=ω2=5, α1=0,α2=-3π ,或普遍表示为Δα=α2-α1=±（2k+1）π ，则合振动振幅为A=A 1-A 2，是两个分振幅之差，合振动仍然为简谐振动。

③ 拍现象当频率差21ωω-很小时，则合振动振幅出现周期变化，此现象称为拍。

合振幅每变化一次周期叫一拍，单位时间内拍出现的次数叫拍频。

例如，取12121210.05,0.04,/3,4, 1.1 4.4,A A a a πωπωωπ=======运行上面程序便可显示出拍现象。

为了图形更好看一些，时间取t=0::10,运行结果如图所示。

为了得到拍效果的声音效果，可以修改以上程序一些参数。

例如，将程序写为clearA1=5;A2=5;w1=*pi;w2=*w1;a1=pi/3;a2=a1;t=0:40000;x1=A1*cos(w1*t+a1);x2=A2*cos(w2*t+a2);x=x1+x2;pause,sound(x1);pause(5),sound(x2);在命令窗口中键入上述程序后按回车键，便可听到三段声音，其中最后一段就是合成后拍的拍音。

（3） 两相互垂直简谐振动的合成 ·题目（ex3322）设一物体同时参加了垂直方向上的两个简谐振动，其简谐振动方程分别为()()111222cos ,cos x A t y A t ωαωα=+=+求其合成振动的轨迹。

·解题分析用计算机画其运动轨迹，不需解轨道方程，只要给定时间数组，算出相应点的x 和y 值，然后作y-x 图即可。

为了能使读者看到图形的绘制过程，程序中使用了for 循环语句以获得动画效果，这对理解图形是很有帮助的。

·程序(ex3322) 讨论：① 频率相同情况两振动方向垂直、频率相同的简谐振动的合振动轨迹为直线、圆或椭圆，轨迹的形状和运动方向由分振动的振幅和相位差决定。

当120,,2a a ππ-=时为直线；取/4,/2,3/4,πππ等时为圆或椭圆。

例如，取12124,0,2/2,a a ωωππ====和则分别得到图3-3-8所示的直线和椭圆。

② 频率不同情况 李萨如图形一般来说，在相互垂直分振动频率不同的情况下，合振动的轨迹不能形成稳定的图案，当如果分振动频率比21/2ωω=成简单整数比，则合振动的轨迹为稳定而闭合的曲线，曲线的花样和分振动的频率比、初相位有关，这些轨迹图形称为李萨如图。

图3-3-9给出了在取121212215,0,/4,2,4(/2)A A a a πωωωω=======即时，由上述程序运行的结果，读者可利用该程序方便地画出李萨如图形。

当两分振动频率比21/ωω是无理数时，合成的运动将永不重复已走过的路径，它的轨迹将逐渐密布在由振幅限定的矩形面积内。

图3-3-10所示的这种非周期运动称为准周期运动。

阻尼振动前面讨论的是不受阻力的自由简谐振动，而振动系统实际上都受阻力的作用，若无外界能量补充，振幅将逐渐减小为零。

这种振动系统因受阻力作振幅减小的运动称为阻尼振动。

·题目(ex3331)一弹簧阵子系统除了受到弹性力外，还受到与速度成正比的黏性阻力的作用，试比较无阻尼、阻尼系数分别为和ω三种情况的运动曲线。

·解题分析设质点在一条直线上，并选择质点的平衡位置为原点，弹簧振子系统中质点受阻力 dxf dtγ=-阻 （γ为阻力系数）和弹性力f kx =-弹 （k 为劲度系数）的作用，根据牛顿第二定律，有22220d x dx x dt dtβω++= 式中，2,,2k m mγωβωβ==和分别为圆频率和阻尼系数。

·程序（ex3331）运行结果如图3-3-11所示。

讨论：对于一定的振动系统，根据阻尼因数大小的不同，由方程可解出三种可能的运动状态，即欠阻尼状态()βω<、过阻尼状态()βω>和临界状态()βω=。

在上述程序中，将频率和阻尼因数写成输入参量，便可作出三种状态的运动图像。

·程序（ex33311）例如，取14.0,s ωβ-=分别取为,和，运行上述程序即可得到图3-3-12。

受迫振动 共振·题目（ex3341）振动系统在连续周期性外力作用下进行的振动叫受迫振动。

考虑质点受三种力作用，即弹性力、阻尼力和驱动力，试讨论受迫振动的运动特征和位移共振。

·解题分析弹性力和阻尼力在前面已经介绍过了，现设驱动力为周期性外力()0cos '.F t F t ω=由牛顿定律，受迫振动方程为2222cos 'd x dx x h t dt dtβωω++= 其中20,,2F k h m m mγωβ=== 初始条件为00,0dxt x dt===时， ·程序（ex3341）运行结构如图3-3-14所示。

讨论：① 迫振动的运动特征由上面的图3-3-14的下图可以看出，开始时，受迫振动的振幅较大，经过一定时间后，阻尼振动即可忽略不计，质点进行与驱动力同频率的振动，称为阻尼受迫振动的稳定状态。

② 位移共振上面的第一幅图给出了受迫振动的振幅随频率比1ωω的变化曲线，可以看出，当驱动力的频率1ω接近系统固有频率时，受迫振动的振幅急剧增大，发生共振现象。

当取不同的阻尼因数时，可以得到一条曲线。

阻尼过大时曲线较为平缓，而减小阻尼时，曲线就变得陡直。

将前面的程序改写为 ·程序（ex33411）分别取B 为,,和并运行该程序得图3-3-14。

读者也可以改变k 值。

看看受迫振动曲线有什么变化，并从物理上给以解释。

非线性振子 混沌现象简谐振动的特征是振子所受的恢复力与位移成正比，方向总指向平衡位置。

但是，有些振动尽管也是呈现周期性变化，而恢复力与位移的关系却是非线性的，平衡位置的均值也将向某一方向偏移。

这种在偏离线性关系的恢复力作用下的振动系统的振动称为非线性振动。

例如，位移较大的弹簧振动、大摆角单摆、复摆运动以及扭摆等都会出现非线性振动。

与线性振动的可重复性和可预知性不同，非线性振动在初始条件给定的情况下不再是可重复的和可预知的，出现了不确定性，这种现象称为混沌。

·题目（ex3351）设非线性振子受1122cos cos F t F t ωω和两个强迫力作用，振子的运动方程为231211222cos cos d x k x k x F t F t dtωω-+=+ 试用MATLAB 求上述非线性方程在两个初始条件相差很小时的解，并画出x-t 图形。

·解题分析 令12,,dxy x y dt==则上述方程可化为一阶微分方程组 123211211122cos cos dy y dtdy k y k y F t F tdtωω⎧=⎪⎪⎨⎪=-++⎪⎩ 为简单起见，可取1212120.5,1,1,0.8F F k k ωω======。

取两个相差微小的初始条件：第一个是x 0=0,v 0=;第二个是x 0=0,v 0=。

分别用这两个初始条件解上述方程，然后将结果绘制在同一幅图中，用以显示初始条件对运动的影响。

·程序（ex3351）运行结果如图3-3-15所示。

可以看出，随时间的演化，原来看似一致的两个运动逐渐分成两个独立的运动，显示出对初始条件的敏感。

混沌行为与粒子的随机运动不同，它告诉我们牛顿定律具有内在的随机性，而不再使我们以前所为的那种确定性理论。

另外，读者还可以选取不同的F 1、F 2以及圆频率值来观察振子的响应。

波动如果在空间某处发生的扰动以一定的速度向四处传播，则称这种传播着的扰动为波。

机械扰动在介质中的传播形成机械波，电磁扰动在真空或介质中的传播形成电磁波。