基于Matlab的二自由度系统模拟实验.docVIP

PAGE PAGE 2 基于Matlab的两自由度振动系统模拟实验报告 实验目的 深入了解两自由度振动系统的模态正交性。 　　2. 掌握Matlab编程基本语言和两自由度系统的响应模态求解方法。  二、 实验原理 如图1所示的系统，设时，两个圆盘恰处于平衡位置, 设。 图1 振动系统模型图 设为的转角，则描述了系统的运动情况。故该系统的自由度为两个。 当不考虑图示系统的阻尼和外界激励时，根据牛顿运动定律,其运动的微分方程为： (1) 其中， 由，可得频率方程如下： (2) 求解式(2)得, (3) 系统振型为，，其振型图如下： 图2 系统振型图 令, 则 , (4) 式(4)为振动系统的模态正交特性。 当考虑图1系统的外界激励时，即时,其运动的微分方程为： (5) 为了求出方程(5)的稳态解,可令,将代入式(5),同时方程两边同乘以,并利用振动系统的模态正交特性,有 (6) 利用式(6)很容易求得, 再由求得。 模拟实验设备和软件 计算机和Matlab6.5版本软件系统 四、实验步骤 、熟悉Matlab语言的基本语句 和基本功能函数 二)、 按照基本原理进行编程和调试 三) 、设置振动系统基本参数进行模拟和运行 五、实验结果和分析 (a) 若,,,,求图1所示的系统固有频率和响应历程图 计算系统固有频率和响应历程图的程序如下: I2=1. kt2=25. F1=5. F2=10. omika=4.5 omika1=sqrt(2*(2-sqrt(2)))*sqrt(kt2/I2)/2. omika2=sqrt(2*(2+sqrt(2)))*sqrt(kt2/I2)/2. U=[1.,1.;sqrt(2),-sqrt(2)]; fs=100. dt=1/100; t=0:dt:1023*dt; q1=(((F1+sqrt(2)*F2)/(4*I2))/(omika1^2-omika^2))*sin(omika*t); q2=(((F1-sqrt(2)*F2)/(4*I2))/(omika2^2-omika^2))*sin(omika*t); for i=1:1024 x1(i)=U(1,1)*q1(i)+U(1,2)*q2(i); x2(i)=U(2,1)*q1(i)+U(2,2)*q2(i); end subplot(2,1,1) plot(t,x1); axis([0,11,-1.,1.]) xlabel(time(s)) ylabel(x1(t)) title(the response of two freedom degree system) subplot(2,1,2) plot(t,x2); xlabel(time(s)) ylabel(x2(t)) axis([0,11,-1.,1.]) end 程序运行可得到系统固有频率分别为2.7060(rad/s), 6.5328(rad/s),系统响应历程见图3。从图3可以看到与相位相同,且的幅值较大。 图3 两自由度系统的响应 (b) 若,,,,求图1所示的系统固有频率和响应历程图 仅需在上面程序基础上改变, 即设置kt2=14.2 程序运行可得到系统固有频率分别为2.0394 (rad/s), 4.9235 (rad/s),系统响应历程见图4。从图4可以看到与相位差180度, 且的幅值较大。 图4 两自由度系统的响应