毕业设计论文-运用Matlab分析机械振动

毕业 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 论文 政研论文下载论文大学下载论文大学下载关于长拳的论文浙大论文封面下载 -运用Matlab 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 机械振动 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 目: 运用Matlab分析机械振动 学院, 物理与电子科学学院 班级, 2008级物理二班 姓名, 李钊 指导教师, 陈翠萍 职称, 实验师 完成日期, 2012 年 5 月 18 日 运用Matlab分析机械振动 作者:李钊 指导老师:陈翠萍 (山西大同大学物理与电子科学学院，山西大同037009) 摘要: 日常生活中所说的振动是一种周期性的运动。所谓周期性运动是指在时间上具有重复性或往复性的一种运动，是遍及自然界及社会科学界的一种运动方式。在物理学中，广义地说凡是描述物质运动的物理量，在某一数值附近做周期性的变化，都叫做振动.本文主要是例举了关于振动的典型实例，用Matlab语言编制计算机程序进行仿真以达到研究简谐振动以及振动的合成，振动的能量等目的，并在文章的最后简要地介绍了共振的危害与应用。 关键字: Matlab语言 演示 振动 周期性 频率 合成 能量 共振 混沌现象 目 录 一(振动的概念与分类 ...................................................... 1 1.1狭义的振动 ......................................................... 1 1.1.1简谐振动的概念............................................................................................................. 1 1.1.2简谐振动的动力学特征 ............................................................................................... 1 1.2广义的振动及振动的分类 ............................................. 2 二(不同类型的振动合成及运用Matlab模拟演示 ............................... 3 2.1振动方向相同，频率相同的简谐振动的合成 ............................. 3 2.2 Matlab模拟振动方向相同，频率略有差异振动合成的“拍”现象 .......... 4 2.3二个振动方向互相垂直的简谐振动的合成 ............................... 7 三(运用Matlab演示典型实例分析简谐振动的能量转换 ......................... 7 3.1简谐振动的系统机械能 ............................................... 7 3.2 小球单摆的能量分析及Matlab演示: ................................. 7 3.3Matlab演示弹簧振子: ............................................... 8 3.4简谐振动的能量曲线 ................................................ 12 四(阻尼振动、受迫振动及位移共振 ......................................... 15 4.1阻尼振动 .......................................................... 15 4.2受迫振动 .......................................................... 18 4.3位移共振 .......................................................... 21 五、“不守规矩”的摆?混沌行为: .......................................... 22 5.1什么是“混沌”现象 ................................................ 22 5.2依赖初值的两种情况 ................................................ 22 六、振动的危害 ........................................................... 22 6.1生产中接触到的振动源 .............................................. 22 6.2振动引起的共振在历史事件中引起的危害 .............................. 23 6.3振动对人体器官的影响 .............................................. 23 七、共振创造了世界 ....................................................... 24 致谢 ..................................................................... 26 一(振动的概念与分类 1.1狭义的振动 1.1.1简谐振动的概念 狭义的振动指的是机械振动，即力学系统中的振动。 振动(或机械振动)指的是物体在平衡位置附近往复运动。(琴弦、钟鼓、机械钟 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 的摆轮、发动机座、高耸的烟囱和固体晶格点阵的分子和原子都在振动) 振动是以波的形式传播的，机械振动的传播即机械波。 简谐振动是指质点在线性回复力作用下围绕平衡位置的运动。 1.1.2简谐振动的动力学特征 (1)简谐振动过程中应该掌握的一些基本概念: 振幅A:简谐运动物体离开平衡位置最大位移的绝对值A。 振动的周期T:物体做一次完全振动所经历的时间。 频率f:单位时间内物体所作的完全振动的次数。 圆频率ω:一秒钟对应的圆心角。一次全振动对应的圆心角就是2π(即360度)。 相位:当振幅和频率一定时决定振动物体在任意时刻相对平衡位置的位移和速度的物理量。 相位概念的重要性还在于比较两个简谐振动之间在“步调”上的差异。设有两个同频率的谐振动，它们的振动表达式为: ,,,cos(，)xAt,111 , ,cos(，)xA,t,222, ,,,(,,,)t，(,,,)它们的相位差为: 2121 (2)以弹簧振子为例描述简谐振动的特点: 弹簧自由伸展时，滑块的位置为原点O(即平衡位置)， 表示位移: 1 平衡位置O点线 X X V X X 图1-1-2 弹簧阵子 2dx由牛顿第二定律知: m,f,,kxx2dt 由上式易知简谐振动的动力学方程 : 2 dx2，,x,0 02 dt x,Acos(,t，,)其解的形式为: 1.2广义的振动及振动的分类 从广义上说振动是指描述系统状态的参量(如位移、电压)在其基准值上下交替变化的过程。 电磁振动习惯上称为振荡。 力学系统能维持振动，必须具有弹性和惯性: ?由于弹性，系统偏离其平衡位置时，会产生回复力，促使系统返回来位置; ?由于惯性，系统在返回平衡位置的过程中积累了动能，从而使系统越过平衡位置 2 向另一侧运动。 正是由于弹性和惯性的相互影响，才造成系统的振动。 (1)按系统运动自由度分，有单自由度系统振动(如钟摆的振动)和多自由度系统振动。 有限多自由度系统与离散系统相对应，其振动由常微分方程描述; 无限多自由度系统与连续系统(如杆、梁、板、壳等)相对应，其振动由偏微分方程描述。 (2)方程中不显含时间的系统称自治系统;显含时间的称非自治系统。 (3)按系统受力情况分，有自由振动、衰减振动和受迫振动。 (4)按弹性力和阻尼力性质分，有线性振动和非线性振动。 (5)振动还可分为确定性振动和随机振动，后者无确定性规律。(如车辆行进中的颠簸) 二(不同类型的振动合成及运用Matlab模拟演示 2.1振动方向相同，频率相同的简谐振动的合成 ,,,cos(，)xAt,111 ,合简谐振动的分振动方程为: ,cos(，)xA,t,222, 振动矢量合成如图示: , A, A2, A12 10xxx21X 图2-1振动矢量合成 x,x，x? 用旋转矢量合成 12 22,,,AAA2AAcos(),，，,121221,,,, ,A,A，A,,? 合振幅矢量 AsinAsin,，121122,tan,,AcosAcos,,，1122, 3 合振动保持原振动方向不变。 x,Acos(,t，,)? 合振动方程 由此易知:一个质点同时参与两个振动方向相同、频率相同的简谐振动，合振动仍为简谐振动。 Matlab模拟编程如下: %两个同方向同频率的简谐振动的合成 clear %清除变量 a1=input('请输入第一个振动的振幅:’); %第一个振动的振幅 %a1=0.03; %参考值 第一个振动的初相 phi1=input('请输入第一个振动的初相的度数:');% %phi1=0; %参考值 phi1=phi1*pi/180; %化为弧度 a2=input('请输入第二个振动的振幅:'); %第二个振动的振幅 %a2=0.04; %参考值 phi2=input('请输入第二个振动的初相的度数:'); %phi1=0;90; phi2=phi2*pi/180; wt=linspace(0,4*pi); x1=a1*cos(wt+phi1); x2=a2*cos(wt+phi2); x=x1+x2; figure plot(wt,x1,'-.',wt,x2,'--',wt,x,'LineWidth',2)%画振动曲线 set(gca,'XTick',(0:8)*pi/2) %设置横坐标刻度 grid on %加网格 fs=16; %字体大小 title('同一直线上简谐振动的合成','FontSize',fs)%显示标题 2.2 Matlab模拟振动方向相同，频率略有差异振动合成的“拍”现象 二个振动方向相同、频率略有差异的简谐振动，其合振动不为简谐振动，产生“拍”现象( 4 拍频为(为两分振动频率) ,,,,,,,,2112 “拍”现象合振动图像如下所示: x1 t x2 t x3 t ,-,21 图2-2“拍”现象振动合成 Matlab演示“拍”现象: %拍的形成 clear %清除变量 d=10; %分母 %d=15; %分母t=0:0.01:60; %时间向量 w1=pi/2; %第一个角频率 dw=pi/d; %角频率之差 w2=w1+dw; %第二个角频率 x1=cos(w1*t); %第一个位移 x2=cos(w2*t); %第二个位移 x=x1+x2; %合位移 figure %创建图形窗口 subplot(3,1,1) %选择子图 plot(t,x1,t,x2,'--','LineWidth',2) %画位移曲线 grid on %加网格 5 leg1='\itx\rm_1/\itA\rm=cos\it\omega\rm_1\itt';%第一个图例字符串 leg2='\itx\rm_2/\itA\rm=cos\it\omega\rm_2\itt';%第二个图例字符串 legend(leg1,leg2) %图例 tit=['(\it\omega\rm_1=\pi/2,\Delta\it\omega\rm=\pi/',num2str(d),')'];%标题一部分 fs=16; %字体大小 title(['拍的形成' tit],'FontSize',fs) %加标题 %xx1=cos(dw*t/2); %调幅线 xx1=cos((w2-w1)*t/2); %调幅线(同上) %xx2=cos((w1+dw/2)*t); %无调幅的振动线 xx2=cos((w2+w1)*t/2); %无调幅的振动线(同上) subplot(3,1,2) %选择子图 plot(t,xx2,t,xx1,'r--','LineWidth',2) %画曲线 grid on %加网格 leg1='cos(\it\omega\rm_2'; %图例字符串的第一部分 leg2='\it\omega\rm_1)\itt\rm/2'; %图例字符串的第二部分 legend([leg1,'+',leg2],[leg1,'-',leg2])%图例 subplot(3,1,3) %选择子图 plot(t,x1+x2,t,2*xx1,'r--',t,-2*xx1,'m--','LineWidth',2)%画曲线 grid on %加网格 xlabel('\itt\rm/s','FontSize',fs) %标记横坐标 ylabel(['\itx/A\rm=\itx\rm_1/\itA\it+\itx\rm_2/\itA'],'FontSize',fs)%标记纵坐标 6 2.3二个振动方向互相垂直的简谐振动的合成 (1)若二振动频率相同，合振动轨迹一般为一椭圆. (2)若二振动频率成整数比，合振动轨迹为规则的稳定的闭合曲线，称利萨如图(但 若不成整数比，轨迹为不闭合的复杂曲线. 三(运用Matlab演示典型实例分析简谐振动的能量转换 3.1简谐振动的系统机械能 弹簧振子或扭摆振动系统中线性回复力为弹性力(或力矩)，它们是保守力(或力 矩)，所以简谐振动系统的总机械能守恒。 简谐振动的总机械能是简谐振动的动能与势能之和 现以单摆、弹簧振子为例讨论振动系统的动能和势能的变化。 3.2 小球单摆的能量分析及Matlab演示: l 2,单摆的周期: g ,h,l(1,cos,)最高点与最低点的高度差: 最高点时动能为0，最低点时势能为0 E,mg,h,mgl(1,cos,)所以振动的能量为: 抽象的问题具体化，运用Matlab演示单摆如下: %制作动画 %挂摆横梁 plot([-0.2;0.2],[0;0],'color','y','linestyle','-',... 'linewidth',10); %画初始位置的单摆 g=0.98; %重力加速度,可以调节摆的摆速 l=1; theta0=pi/4; x0=l*sin(theta0); y0=(-1)*l*cos(theta0); axis([-0.75,0.75,-1.25,0]); 7 axis('off'); %不显示坐标轴 %创建摆锤 head=line(x0,y0,'color','r','linestyle','.',... 'erasemode','xor','markersize',40); %创建摆杆 body=line([0;x0],[0;y0],'color','b','linestyle','-',... 'erasemode','xor'); %摆的运动 t=0; dt=0.01; while 1 t=t+dt; theta=theta0*cos(sqrt(g/l)*t); x=l*sin(theta); y=(-1)*l*cos(theta); set(head,'xdata',x,'ydata',y); set(body,'xdata',[0;x],'ydata',[0;y]); drawnow; end 3.3Matlab演示弹簧振子: 11222E,kx,kAcos(,t，,)弹簧振子的弹性势能: P22 112222E,mv,m,Asin(,t，,) 弹簧振子的动能: k22 2k,m,,k？, m 122？E,kAsin(,t，,) k2 8 12EEEkA,，,弹簧振子的总机械能: kp2 122EmA,, 2 1122EkxE,kA,因为 ,均较易进行计算，所以计算动能时常用P22 E,E,E kp 综上所述易知: ?任何简谐振动系统的机械能均可用下式计算 122E,kAcos(ωt，,)p2 122E,kAsin(ωt，,)k2 12EEEkA,，,kp2 ?简谐振动过程中，系统机械能守恒，但动能和弹性势能相互转化。 12EEEkA,，,kpE,Epk2 ?简谐振动的振幅与机械能的关系。 E2E20AA,,kk 理想弹簧振子的简谐振动Matlab编程如下: %理想弹簧阵子简谐运动 %Clear rectangle('position',[12,8.5,2,0.3],'FaceColor',[0.5,0.3,0.4]); axis([0,15,-1,10]); 9 %画顶板 hold on plot([13,13],[7,8.5],'r','linewidth',2); %画直线 y=2:.2:7; M=length(y); x=12+mod(1:M,2)*2; x(1)=13; x(end-3:end)=13; D=plot(x,y); %弹簧 C=0:.1:2*pi;r=0.35; t1=r*sin(C); F1=fill(13+r*cos(C),2+t1,'r'); % 球 10 set(gca,'ytick',[0:2:9]); set(gca,'yticklabels',num2str([-1:3]')); plot([0,15],[3.3,3.3],'black'); H1=plot([0,13],[3.3,3.3],'y'); % 句柄[黄线] Q=plot(0,3.8,'color','r'); % 运动曲线; td=[];yd=[]; T=0; text(2,9,'理想中的弹簧振子简谐振动','fontsize',16); set(gcf,'doublebuffer','on'); while T<12; pause(0.2); Dy=(3/2-1/2*sin(pi*T))*1/2; Y=-(y-2)*Dy+7; 11 Yf=Y(end)+t1; td=[td,T];yd=[yd,Y(end)]; set(D,'ydata',Y); set(F1,'ydata',Yf,'facecolor',rand(1,3)); set(H1,'xdata',[T,13],'ydata',[Y(end),Y(end)]); set(Q,'xdata',td,'ydata',yd) ; T=T+0.1; End 3.4简谐振动的能量曲线 ?能量曲线: 12EEEkA,，, 总机械能: kp2 12 2kx, 弹性势能能: Ep2 E,E,E 动能: kp 能量曲线如下图所示: EEEpk E Ek x EP 图3-4简谐振动的能量曲线 ?弹性势能与动能的平均值: 简谐振动中势能与动能的平均值相等且等于总机械能的一半。 TT1111222,,,，,,EEdtkAcos(t)dtkAPP,,00TT24 TT1111222,,,，,,EEdtkAsin(t)dtkAkk,,00TT24 ?以弹簧振子为例运用Matlab模拟能量曲线: %弹簧振子的动能,势能和机械能曲线 clear %清除变量 n=4; %周期的个数 t=linspace(0,2*pi)*n; %时间向量 x=cos(t); %振子位置 v=-sin(t); %速度 ek=v.^2; %动能 ep=x.^2; %势能 figure %建立图形窗口 subplot(2,1,1) %子图 13 plot(t,x,t,v,'--','LineWidth',2) %画位移和速度曲线 grid on %加网格 axis tight %紧贴坐标轴 fs=16; %字体大小 title('简谐振动的位移和速度','FontSize',fs)%显示标题 xlabel('\it\omegat','FontSize',fs) %标记横轴 \itx/A','速度\itv/\omegaA')%图例 legend('位移 set(gca,'XTick',(0:2*n)*pi) %设置横坐标刻度线 text(0,0,'\it\omega\rm=(\itk/m\rm)^{1/2}','FontSize',fs)%显示角频率 subplot(2,1,2) %子图 plot(t,ek,'--',t,ep,'-.',t,ek+ep,'LineWidth',2)%画能量曲线 grid on %加网格 axis tight %紧贴坐标轴 title('简谐振动的能量','FontSize',fs) %显示标题 xlabel('\it\omegat','FontSize',fs) %标记横轴 legend('动能\itT/E\rm_0','势能\itV/E\rm_0','机械能\itE/E\rm_0')%图例 text(0,0.5,'\itE\rm_0=\itkA\rm^2/2','FontSize',fs)%显示能量单位 set(gca,'XTick',(0:2*n)*pi) %设置横坐标刻度线 14 四(阻尼振动、受迫振动及位移共振 4.1阻尼振动 以上讨论均假设质点或刚体的振动不受任何阻力，由于能量守恒，它们将永远振动下去。然后事实上，振动系统都受阻力作用，如无外界能量补偿，振动幅将不断减小而归于静止。振动系统因受阻力作振幅减小的运动，叫做阻尼振动。 dx 设 阻力Fv,,,,,,xxdt 2dxdx由牛顿第二定律得 kx m,,,,2dtdt k,2 令 ,,,,0m2m 2dxdx2由此可得其动力学方程 ，2,，,x,002dtdt 根据阻尼因数之不同，可将此方程解出三种可能的运动状态: x , (1)欠阻尼状态 t ,，0,，,0 ,,t,x,Aecos(,t，,) 得质点的解 图4-1-1欠阻尼状态 x (2)临界阻尼 t ,,,,,00 ,,tx,(c，ct)e得质点的解 12 图4-1-2临界阻尼状态 x (3)过阻尼状态 ,，0,，,0 t 2222,(,,,,,)t,(,，,,,)t00x,ce，ce 得质点的解 12 图4-1-3过阻尼状态 15 运用Matlab语言模拟阻尼运动: %阻尼运动的类型 clear %清除变量 t=0:0.25:20; %固有角频率与时间的乘积w0t向量(约化时间向量) %b=0:0.5:1.5; %阻尼因子与固有角频率的倍数向量(约化阻尼因子向量) b=0:0.25:1.25; %阻尼因子与固有角频率的倍数向量(约化阻尼因子向量) n=length(b); %曲线条数 b(b==1)=1+eps; %将1值改为1加小量 [B,T]=meshgrid(b,t); %约化阻尼因子和约化时间矩阵 W=sqrt(1-B.^2); %准角频率矩阵 X=exp(-B.*T).*(cos(W.*T)+B./W.*sin(W.*T));%位移函数矩阵 V=-exp(-B.*T).*sin(W.*T)./W; %速度函数矩阵 %A=sqrt(B.^2-1); %参数矩阵 %X=exp(-B.*T).*((A+B).*exp(A.*T)+(A-B).*exp(-A.*T))/2./A;%位移函数矩阵(效果相同) %X=exp(-B.*T).*(cosh(A.*T)+B./A.*sinh(A.*T));%位移函数矩阵(效果相同) %V=-exp(-B.*T).*sinh(A.*T)./A; %速度函数矩阵(效果相同) f1=figure; %创建图形窗口 %plot(t,X,'LineWidth',2) %画位移曲线族 plot(t,X(:,1),'o-',t,X(:,2),'d-',t,X(:,3),'s-',t,X(:,4),'p-',... t,X(:,5),'h-',t,X(:,6),'<-') %画位移曲线族 fs=16; %字体大小 xlabel('\it\omega\rm_0\itt','FontSize',fs)%标记横坐标 ylabel('\itx/A','FontSize',fs) %标记纵坐标 title('质点在不同阻尼下的运动曲线','FontSize',fs)%标题 legend([repmat('\it\beta/\omega\rm_0:',n,1),num2str(b')])%加图例 txt='\it\beta/\omega\rm_0 小于1为欠阻尼,等于1为临界阻尼,大于1为过阻尼';%文本 text(0,-0.7,txt,'FontSize',fs) %显示文本 grid on %加网格 16 f2=figure; %创建图形窗口 %plot(t,V,'LineWidth',2) %画速度曲线族 plot(t,V(:,1),'o-',t,V(:,2),'d-',t,V(:,3),'s-',t,V(:,4),'p-',... t,V(:,5),'h-',t,V(:,6),'<-') %画位移曲线族 xlabel('\it\omega\rm_0\itt','FontSize',fs)%标记横坐标 ylabel('\itv/\omega\rm_0\itA','FontSize',fs)%标记纵坐标 title('质点在不同阻尼下的速度曲线','FontSize',fs)%标题 grid on %加网格 legend([repmat('\it\beta/\omega\rm_0:',n,1),num2str(b')])%加图例 pause %暂停,可取图 X1=[]; %位移矩阵清空 V1=[]; %速度矩阵清空 X2=[]; %位移矩阵清空 V2=[]; %速度矩阵清空 for i=1:n %按曲线循环 [tm,XV]=ode45('P5_7fun',t,[1;0],[],b(i));%计算位移和速度 X1=[X1,XV(:,1)]; %连接位移矩阵 V1=[V1,XV(:,2)]; %连接速度矩阵 s=['D2x+',num2str(2*b(i)),'*Dx+x'];%微分方程字符串 sx=dsolve(s,'x(0)=1','Dx(0)=0'); %微分方程的符号解 sv=diff(sx); %求速度的符号解 x=subs(sx,'t',t); %位移 v=subs(sv,'t',t); %速度 X2=[X2,x']; %连接位移矩阵 V2=[V2,v']; %连接速度矩阵 end %结束循环 figure(f1) %重开图形窗口 hold on %保持图像 plot(t,X1,'.',t,X2,'*') %画位移曲线 figure(f2) %重开图形窗口 17 hold on %保持图像 plot(t,V1,'.',t,V2,'*') %画速度曲线 %阻尼运动的二阶微分方程的函数 function f=fun(t,x,flag,b) f=[ x(2); %速度 -2*b*x(2)-x(1)]; %加速度 4.2受迫振动 dx设质点受到三种力:弹性力-kx，阻尼力，周期性外力，亦称驱动力-,dt F(t),Fcos,t0 根据牛顿第二定律得，受迫振动的动力学方程: 2dxdx m,,kx,,，Fcos,t02dtdt Fk,20令 ,，，f ,2,,,00mmm 2dxdx2得 ，2,，,x,fcos,t002dtdt ,,t,x,Aecos(,t，,)，Acos(,t，,)解方程得: 0 开始时，受迫振动的振幅较小，经过一定时间后，阻尼振动即可忽略不计，质点进 行由上式第二项所决定的与驱动力同频率的振动，称受迫振动的稳定振动状态，可表示 18 2dxdx2x,Acos(,t，,)，2,，,x,fcos,t如下:将此式带入方程得 0002dtdt f,,,20,A,， ,tan02222222,,,(,,,)，4,,00 x o t 图4-2 受迫振动自暂态发展为稳定振动.本图所示初始条件为t=0,,,初始条v,0x,000x 件影响暂态过程，不影响稳态振动 Matlab编程演示受迫振动如下: %物体在平衡点从静止开始的受迫振动曲线(用解析式) clear %清除变量 b=input('请输入约化阻尼因子(0~1):'); %键盘输入约化阻尼因子 %b=0.1; %参考值 if b<=0|b>=1 return,end %不符合条件则不向下执行程序 w=sqrt(1-b^2); %约化阻尼圆频率 s=['请输入约化驱动力圆频率(约化阻尼圆频率为',num2str(w),'):'];%提示字符串 W=input(s); %键盘输入约化驱动力圆频率 %W=2;6;1;0.6; %参考值 if W==1 W=1-eps;end %如果为1则改小一点 tm=30; %最大时间 t=0:0.001:tm; %时间向量 19 a1=sqrt(w^2*(W^2-1)^2+b^2*(W^2+1)^2)/w/((W^2-1)^2+4*b^2*W^2);%阻尼振幅 phi=atan2(b*(W^2+1),w*(W^2-1)); %阻尼振动初相 a2=1/sqrt((W^2-1)^2+4*b^2*W^2); %等幅振动振幅 PHI=atan2(-2*b*W,1-W^2); %等幅振动初相 x1=a1*exp(-b*t).*cos(w*t+phi); %阻尼振动函数 x2=a2*cos(W*t+PHI); %等幅振动函数 x=x1+x2; %合成振动 xm=max(abs(x)); %最大值 %----------------------------------------------------------- %创建图形窗口 figure subplot(3,1,1) %选子图 plot(t,x1,'LineWidth',2) %画曲线 grid on %加网格 axis([0,tm,-xm,xm]) %设置曲线范围 fs=12; %字体大小 title('减幅振动的位移时间曲线','FontSize',fs)%标题 ylabel('\itx\rm_1/\itA\rm_0','FontSize',fs)%标记纵坐标 txt=['\it\beta/\omega\rm_0=',num2str(b)];%阻尼因子字符串 txt=[txt,',\it\omega/\omega\rm_0=',num2str(w)];%连接阻尼圆频率 text(0,xm,txt,'FontSize',fs) %标记阻尼因子和阻尼圆频率 subplot(3,1,2) %选子图 plot(t,x2,'LineWidth',2) %画曲线 grid on %加网格 axis([0,tm,-xm,xm]) %设置曲线范围 title('等幅振动的位移时间曲线','FontSize',fs)%标题 ylabel('\itx\rm_2/\itA\rm_0','FontSize',fs)%标记纵坐标 txt=['\it\Omega/\omega\rm_0=',num2str(W)];%驱动力圆频率字符串 text(0,xm,txt,'FontSize',fs) %标记驱动力圆频率 subplot(3,1,3) %选子图 plot(t,x,'LineWidth',2) %画曲线 20 grid on %加网格 axis([0,tm,-xm,xm]) %设置曲线范围 txt='\itA\rm_0=\itF\rm_0/\itm\omega\rm_0^2';%振幅文本 text(0,xm,txt,'FontSize',fs) %标记振幅文本 title('受迫振动的位移时间曲线','FontSize',fs)%标题 xlabel('\it\omega\rm_0\itt','FontSize',fs)%标记横坐标 ylabel('\itx/A\rm_0','FontSize',fs) %标记纵坐标 4.3位移共振 任何物体产生振动后，由于其身的构成、大小、形状等物理特性，原先以多种频率开始的振动，渐渐会固定在某一频率上振动，这个频率叫做该物体的"固有频率"，因为它与该物体的物理特性有关。所以当人们从外界再给这个物体加上一个振动(称为策动)时，如果策动力的频率与该物体的固有频率正好相同，物体振动的振幅达到最大，这种现象叫做"共振"。 任何物体自身存在振动，当一个物体接受到另一个物体的振荡频率时，又巧好与这个物体的振荡频率相同时会产生共振。共振危害极大可以使大桥、房屋以及其它的建筑物瞬间倒塌，甚至还危及到人类的心脏使血管破裂而亡。 下面简要介绍一下位移共振: 当驱动力频率取某值时，振幅获得极大值。振动系统受迫振动时，其振幅达到极大 f0值的现象叫做位移共振。将式用微分法关于极大值的判据，可求A,022222(,,,)，4,,0 22出共振驱动力的圆频率为:亦称位移共振条件。 ,,,,2,r0 21 Ao ,,0 ,r 图4-3由于阻尼存在，位移共振时受迫振动的频率不等于驱动力频率，仅当阻尼无限小时，共振频率无限接近固有频率，产生极激烈的共振 五、“不守规矩”的摆?混沌行为: 5.1什么是“混沌”现象 23dxdx，,,x，x,,cos,t2动力学方程形如:的摆动现象称为混沌摆。与受迫dtdt3x振动有所区别的是出现了非线性项——，因此这扭摆称为受周期驱动的非线性振子。非线性系统的混沌行为是自然界普遍和重要的物理现象 振子的运动是“循规蹈矩”的，但与之相对比的是，这里的扭摆当参量取某些值时，却表现出“不守规矩”的行为。长期行为不可预测，呈现随机性。扭摆的这种行为称为混沌行为或混沌运动，也称作混沌。 5.2依赖初值的两种情况 依赖初值的两种情况:牛顿力学研究的运动是依赖初值的。当今。混沌研究表明，运动对初值的依赖本应该分为两类: 一类是一般的依赖初值;这类运动，给定初值，之前和之后的运动是完全确定的，这类运动是可“重现”，可“预报”的。 另一类是敏感地依赖初值;“运动依赖敏感初值”意味着初值有预测不出的微小微差，随着时间演化，两次运动间的偏离却可以观测得到，并且还变化不定，也就是说对于敏感依赖初值的运动。即便在实验中给予“相同的”初值(注意: “相同的”表示两次初值相差如此之小，以致测量不出两次初值的不同)两次运动也并不重复，即敏感依赖初值的运动是不可重复出现的，不可预报的。 六、振动的危害 6.1生产中接触到的振动源 (1)凿岩机、铆钉机、风铲等风动工具; ( 2 ) 砂轮机、抛光机、研磨机、电钻、电锯、林业用油锯等电动工具; 22 (3)摩托车、内燃机车、船舶等运输工具; (4)收割机、拖拉机、脱粒机等农业机械。 6.2振动引起的共振在历史事件中引起的危害 1831年，一队骑兵通过英国曼彻斯特附近的一座吊桥。突然，随着一声巨响，大桥断裂崩塌了，人与马纷纷坠入河中，导致死伤惨重。无独有偶，过了半个多世纪，1906年，俄国首都彼得格勒也有一支全副武装的沙皇军队，步伐整齐的通过爱记华特大桥时，大桥莫名其妙的倒塌了，当时情景狼狈不堪。 两件事情发生的方式大同小异，人们对此觉得特别震惊，并随即对此进行了调查。通过当时一大批顶尖物理学家的研究发现，在没有敌人破坏，又不是桥的质量问题时，肇事者正是这些受害者自己。由于他们步伐整齐产生的周期频率碰巧接近桥的固有频率，激起了大桥的共振，结果造成了桥断人亡的大事故～ 6.3振动对人体器官的影响 ( 1) 条件反射潜伏期改变;交感神经功能亢进;血压不稳、心律不齐等; 如触觉、温热觉、痛觉等皮肤感觉功能下降，特别是振动感觉，最早出现迟钝。 (2)40,300Hz的振动能引起周围毛细血管张力和形态的变化，症状表现为末梢血管痉挛、脑血流图异常等;心脏方面会出现心律不齐、心动过缓等病症。 ( 3) 振动会引起肌电图异常、握力下降、肌肉萎缩、肌纤维颤动和疼痛等病症。 (4)小于40Hz的大振幅振动容易引起骨关节的变动，通过X光，可清晰的观察到骨贸形成、骨质疏松、骨关节变形和骨关节坏死等症状。 ( 5)振动引起的听力变化，则以125,250Hz频段的听力下降最为明显，然而在初期仍以高频段听力损失为主，随后才会出现低频段听力下降。 振动和噪声有联合作用。 (6)长期处于振动工作环境的人可会患发局部振动病。局部振动病是以末梢循环障碍为主的疾病，也会影响肢体神经及运动功能。发病部位多表现于上肢末端，典型性为发作性手指变白(亦称白指)。 我国1957年就将局部振动病定为职业病。 ( 7 )众所周知，振动频率、加速度与振幅才是影响振动作用的因素。然而人体只会对频率在1,1000Hz范围的振动产生振动感觉。振动频率在发病过程中起重要作用。 频率在30,300Hz范围内的振动主要是引起末梢血管痉挛，白指。 频率相同时，危害随着加速度的增大而增大。而振幅大，频率低的振动主要作用于 23 前庭器官，且会使内脏发生移位。 频率一定时，振幅越大对机体影响越大。这是寒冷振动病发病的重要外部病因之一，寒冷会导致血流量减少，血液循环发生改变，造成局部供血不足，促进振动病发生。接触振动时间越长，振动病发病率越高。增长工间休息时间对预防振动病具有积极意义。 人对振动的敏感程度与身体部位有关。人体站立时对垂直振动敏感;躺卧时对水平振动敏感。有的作业要采取强制体位，甚至胸腹部或下肢部位贴振动体，这样造成的振动危害就更大了。加工部件硬度大时，人体所受危害大也会相应增大，冲击力大的振动容易导致骨关节发生病变。 七、共振创造了世界 有利必有弊，世界上的任何事物都逃不过自然规律，振动不只是能给人们带来危害，而且只要人们运用得当，振动还能给人们带来莫大的好处。共振筛和垂直输送器就是利用振动的共振原理最典型的例子。 振动中的共振是物理学上的一个运用频率非常高的专业术语。 共振在声学中称为“共鸣”，它是指物体因共振而发声的现象; 共振在电学中称为“谐振”。 产生共振的重要的条件之一是要有弹性，而且一件物体受外来的频率作用时，它的频率要与后者的频率相同或相接近。从总体来说，宇宙上大多数物质是具有弹性的，大到行星小到原子，几乎都能以一个或多个固有频率来振动。因此也可以这么说，共振现象是宇宙间最普遍和最频繁的现象之一。大家都知道宇宙是在一次剧烈的大爆炸后产生的，而促使这次大爆炸的根本原因之一就是共振。因而从某种程度来说，是共振产生了宇宙和世界万物，没有共振就没有世界。 共振在日常生活中的应用也是不可或缺的一部分，例如: (1)石油化工存在很多共振的地方，例如精馏塔的塔盘共振; (2)音响设备中扬声器纸盆的振动，各种弦乐器中音腔在共鸣箱中的振动等利用“力学共振” ; (3)电磁波的接收和发射利用了“电磁共振”; (4)医疗技术中非常普及的“核磁共振”等等。 参考文献: 24 [1] 漆安慎 杜婵英 编《力学》第二版 北京，高等教育出版社，2005年6月(2008重印) [2] 程守洙 江之永 主编 胡盘新 汤毓骏 宋开欣 修订《普通物理》第三版 北京，高等教育出 版社，1998年6月 [3] 东南大学等七所工科院校 编，《物理学》，【M】， 北京，高等教育出版社，2008年12月 [4] 陈怀琛 编《MATLAB》第三版 西安电子科技大学出版社，2007年7月(2010.6重印) [5]《大科技》杂志 : 《共振，能把地球一烈为二》 [6] 李 斌 《共振的危害与利用》 [7] 蔡圣君 《共振的运用与防止》 Using Matlab Analyzing Mechanical Vibration Abstract: The daily life of the said vibration is a cyclical movement. The so-called periodic motion is to point to in time is to summarize repetitive or a sport. In physics, general description of all material movement physical quantities, in a numerical around cyclical change, all are called vibration.This paper mainly introduced the vibration of the typical example, using Matlab language program by computer simulation to study simple harmonic vibration and vibration synthesis and vibrational energy purposes, and at the end of the article briefly introduces the harm and application of resonance. Key words: Matlab language Demo vibration of periodic resonance frequency synthesis energy chaos 25 致谢 本论文从立题到论文撰写整个过程都是在指导教师陈翠萍老师的悉心指导下完成的。特别是指导教师在学习上，思想上都给予了我极大的关怀和帮助，在传授我知识的同时，更注重培养我解决问题的思路和方法及创新能力，为我今后学习和工作打下了坚实的基础并开阔了我的视野。指导教师敏捷的思维和孜孜不倦的探索精神是我永远学习 我将最诚挚的的榜样，我所取得的每一点进步无不凝聚着指导教师陈翠萍老师的心血,谢意奉献给我的指导教师陈翠萍老师。 26