机械振动模型.PPT

* * 机械振动的模型 有一个弹簧,它的上端总是固定,下端挂个质量为m的物体,当物体处于静止的状态时,作用在物体上的重力与弹性力大小相等,方向向反.这个位置就是物体的平衡位置. 当物体处于平衡位置的时候,受到向下的重力为mg,弹簧的向上弹力k△l的作用,其中是k弹簧的弹性系数,是弹簧的受重力mg作用后向下拉伸的长度,有k△l=mg. 为了研究物体的运动规律,选取平衡位置为坐标原点, 取x轴垂直向下.从而当物体处于平衡位置的时候,有x=0,但物体受到的外力F(t)作用时,从平衡位置开始运动,x(t)代表物体在t时的位置,当物体开始运动时,受到下面4个力的作用. (1)物体的重力W=mg,方向是向下的，与坐标轴的方向一致. (2)弹簧的弹力R,当△l+x0时,弹力与x轴的方向相反R=-k(△l+x), 当△l+x0时,弹力的与x轴的方向相同,即R=-k(△l+x).因此弹簧的弹力总有R=-k(△l+x). (3)空气阻力D,物体在运动过程中总受到空气或其他介质的阻力作用,使震动逐渐减弱,阻力的大小与物体的运动速度成正比,方向与运动的方向相反,该阻力系数为c,在时刻物体运动速度为diff(x(t),t) ,因此D=-c*diff(x(t),t) (4)物体在运动过程中还受到随时间变化的外力F(t)作用,方向可能不确定.依赖的F(t)正负. .因此物体的运动满足二阶线形微分方程 m*diff(diff(x(t),t),t)+c*diff(x(t),t)+k*x(t)=F(t). (A).我们首先研究没有空气阻力和外力作用的弹簧振动.这个时 候方程变为 diff(x(t),t),t)+k*x(t)=0 令w0^2=k/m.方程的解为 x(t) = _C1*sin(wo*t)+_C2*cos(wo*t) C1,C2为常数.令sin(theta)=C1/ sqrt(C1^2+C2^2), cos(theta) =C2/ sqrt(C1^2+C2^2)^0.5. 取A= sqrt(C1^2+C2^2),theta=arctan(C1/C2), X(t)= sqrt(C1^2+C2^2)*(C1/ sqrt(C1^2+C2^2)*cos(w0*t) + C2/ sqrt(C1^2+C2^2)*sin(w0*t) ). = Asin(w0*t+theta) 令w0=4,theta=0,A=1.我们可以看到物体的运动是周期震动的.物体的振幅是A,theta是初相位,振幅和初相位都依赖于初始条件的选择. (B).现在我们考虑有空气阻力而没有外力作用的弹簧震动.这个时候微分方程变为 m*diff(diff(x(t),t),t)+c*diff(x(t),t)+k*x(t)=0; 上面方程的特征方程为 m*u^2+c*u+k=0; 特征根为 u1 = 1/2/m*(-c+sqrt(c^2-4*m*k)), u2 = 1/2/m*(-c-sqrt(c^2-4*m*k)); 方程的解为 x(t) = _C1*exp(1/2/m*(-c+sqrt(c^2-4*m*k))*t)+_C2*exp(-1/2*(c+sqrt(c^2-4*m*k))/m*t) 这里我们可以把方程的通解写成 x(t)= A*exp(-c*t/(2*m))*sin(u*t+theta); 我们令A=c=m=1,theta=0,u=10;那么方程的图形为 我们可以看出弹簧的震动已经不是周期的了,震动的最大偏离A*exp(-c*t/(2*m))随着时间的增加而不断的减小,如图所 给的,最后趋于平衡位置x=0. (C)假如物体在运动的过程中既有空气的阻力又有周期外力F(t) 作用,假设F(t)=Fo^cos(w*t),则弹簧震动的方程为 m*diff(diff(x(t),t))+c*diff(x(t),t)+k*x= Fo^cos(w*t). 这里求的方程的一个特解为 g(t)= F0*sin(w*t+theta)/[(k-m*w^2)^2+c^2*w^2]^0.5 这里tg(theta)=(k-m*w^2)/(c*w) 因此方程的通解为x(t)= A*exp(-c*t/(2*m))*sin(u*t+theta)+ F0*sin(w*t+theta)/ sqrt[(k-m*w^2)^2+c^2*w^2]; 令A=w=k=m=w=c=F0=1,u=10; 这里方程变为 x(t)= exp(-t/2)*sin(10*t)+ sin(t); 这里我