MATLAB与系统仿真 - 图文

五、结合《机械工程设计分析M文件集(第2版)》中的例子，同时结合自己的本专业所学到的知识，设计一个MATLAB程序用来解决与自己学到的机械学科知识有关的实际问题。（本题10分）

【解】

机械运动中，槽轮机构再各种自动机械中的运用很广泛，槽轮机构是间歇性运动，槽轮再运转中有较大的动载荷，槽轮的运动过程分析也是一项很重要的工作，运动参数的设计是必不可少的。下面在MATLAB系统中就槽轮机构的运动进行分析：

新建一个M-file文件，输入下面的程序，保存为dfgsg.mdl % ?ú2????ú11???ˉ·???

dr=pi/180.0; % ???èó????èμ?×a???μêy % ?ú??2×a??·??§:-f20

lmd=sin(pi/z); % ?????ú±ú2ó??ú?ü1μ?3¤?è±è bc=10; % ?-?·2?3¤ cz=-f20/dr; % ?-?·3??μ zz=f20/dr; % ?-?·???μ

i=1; % ?ù?Y2?3¤±??ˉμ????ˉ2?êy???ócsDDêy??êy?÷ for f2=cz:bc:zz % ????2???????ò??￠àà???ù?è?￠àà???ó?ù?è wy=atan(lmd*sin(f2*dr)/(1+lmd*cos(f2*dr)));

sd=lmd*(cos(f2*dr)+lmd)/(1+2*lmd*cos(f2*dr)+lmd^2);

jsd=lmd*sin(f2*dr)*(1-lmd^2)/(1+2*lmd*cos(f2*dr)+lmd^2)^2; switch z % ???óc(i,:)±íê?μúiDDμ??÷áD?a?? case 4,c4(i,:)=[f2 wy/dr sd jsd]; case 6,c6(i,:)=[f2 wy/dr sd jsd]; case 8,c8(i,:)=[f2 wy/dr sd jsd]; case 10,c10(i,:)=[f2 wy/dr sd jsd]; end i=i+1; end end

% ê?3??ú2????ú11???ˉ2?êy ['??2?êy z=4']

[' ?ú??×a??',' 2???????ò?',' ???ù?è',' ???ó?ù?è'] % ???óc(:,j)±íê?μújáDμ??÷DD?a?? [c4(:,1),c4(:,2),c4(:,3),c4(:,4)] ['??2?êy z=6']

[' ?ú??×a??',' 2???????ò?',' ???ù?è',' ???ó?ù?è'] % ???óc(:,j)±íê?μújáDμ??÷DD?a?? [c6(:,1),c6(:,2),c6(:,3),c6(:,4)]

20

['??2?êy z=8']

[' ?ú??×a??',' 2???????ò?',' ???ù?è',' ???ó?ù?è'] % ???óc(:,j)±íê?μújáDμ??÷DD?a?? [c8(:,1),c8(:,2),c8(:,3),c8(:,4)] ['??2?êy z=10']

[' ?ú??×a??',' 2???????ò?',' ???ù?è',' ???ó?ù?è'] % ???óc(:,j)±íê?μújáDμ??÷DD?a?? [c10(:,1),c10(:,2),c10(:,3),c10(:,4)] %

% ????2????ú11???ˉ2?êy?ú??

figure(1); % éú3é2??????ˉ??í?′°?ú subplot(2,2,1); % ????μú1??×ó′°?ú plot(c4(:,1),c4(:,3),c4(:,1),c4(:,4)) % ????z= 4μ???í? title('?ú2???2?êy z=4') % ±ê×￠×ó′°?ú??3? axis([-3*pi/4/dr 3*pi/4/dr -1 1]) % ?¨ò?×?±ê?á·??§ grid % ?¤????

text(-85,-0.2,'\\epsilon/\\omega^{2}') % ±ê×￠àà???ó?ù?è??í? text(10,0.55,'\\omega/\\omega') % ±ê×￠àà???ù?è??í? ylabel('2??????ˉ??í?') % ?¨ò?×Y×?±ê?á??3? %

subplot(2,2,2); % ????μú2??×ó′°?ú plot(c6(:,1),c6(:,3),c6(:,1),c6(:,4)) % ????z= 6μ???í? title('?ú2???2?êy z=6')

axis([-3*pi/4/dr 3*pi/4/dr -0.6 0.6]) grid

text(-85,-0.3,'\\epsilon/\\omega^{2}') text(10,0.4,'\\omega/\\omega') ylabel('2??????ˉ??í?') %

subplot(2,2,3); % ????μú3??×ó′°?ú plot(c8(:,1),c8(:,3),c8(:,1),c8(:,4)) % ????z= 8μ???í? title('?ú2???2?êy z=8')

axis([-3*pi/4/dr 3*pi/4/dr -0.40 0.40]) grid

text(-85,-0.3,'\\epsilon/\\omega^{2}') text(10,0.32,'\\omega/\\omega') ylabel('2??????ˉ??í?') %

subplot(2,2,4); % ????μú4??×ó′°?ú plot(c10(:,1),c10(:,3),c10(:,1),c10(:,4)) % ????z=10μ???í? title('?ú2???2?êy z=10')

axis([-3*pi/4/dr 3*pi/4/dr -0.35 0.35]) grid

text(-80,-0.22,'\\epsilon/\\omega^{2}')

21

text(10,0.28,'\\omega/\\omega') ylabel('2??????ˉ??í?') %

figure(2); % éú3éàà??í?′°?ú subplot(1,2,1); % ????μú1??×ó′°?ú

plot(c4(:,1),c4(:,3),c6(:,1),c6(:,3),c8(:,1),c8(:,3),c10(:,1),c10(:,3))

title('\\omega/\\omega')

axis([-3*pi/4/dr 3*pi/4/dr -0.05 0.45]) grid

text(-12,0.43,'z=4') text(-12,0.35,'z=6') text(-12,0.29,'z=8') text(-15,0.21,'z=10') ylabel('2???àà???ù?è??í?') %

subplot(1,2,2); % ????μú2??×ó′°?ú

plot(c4(:,1),c4(:,4),c6(:,1),c6(:,4),c8(:,1),c8(:,4),c10(:,1),c10(:,4))

title('\\epsilon/\\omega^{2}')

axis([-3*pi/4/dr 3*pi/4/dr -0.6 0.6]) grid

text(40,0.22,'z=10') text(-60,-0.15,'z=8') text(-110,-0.45,'z=6') text(60,0.05,'z=4')

ylabel('2???àà???ó?ù?è??í?')

运行该程序可以得到以下参数依据： ans =

轮槽数 z=4

ans =

销轮转角 槽轮角位移 角速度 角加速度

ans =

-135.0000 -45.0000 0 -1.0000 -125.0000 -44.2584 0.1371 -0.6104 -115.0000 -42.4270 0.2229 -0.3936 -105.0000 -39.8961 0.2795 -0.2656

22

-95.0000 -36.8948 0.3184 -0.1858 -85.0000 -33.5652 0.3460 -0.1337 -75.0000 -30.0000 0.3660 -0.0981 -65.0000 -26.2621 0.3808 -0.0728 -55.0000 -22.3962 0.3918 -0.0542 -45.0000 -18.4349 0.4000 -0.0400 -35.0000 -14.4035 0.4060 -0.0287 -25.0000 -10.3217 0.4101 -0.0193 -15.0000 -6.2060 0.4128 -0.0111 -5.0000 -2.0708 0.4141 -0.0036 5.0000 2.0708 0.4141 0.0036 15.0000 6.2060 0.4128 0.0111 25.0000 10.3217 0.4101 0.0193 35.0000 14.4035 0.4060 0.0287 45.0000 18.4349 0.4000 0.0400 55.0000 22.3962 0.3918 0.0542 65.0000 26.2621 0.3808 0.0728 75.0000 30.0000 0.3660 0.0981 85.0000 33.5652 0.3460 0.1337 95.0000 36.8948 0.3184 0.1858 105.0000 39.8961 0.2795 0.2656 115.0000 42.4270 0.2229 0.3936 125.0000 44.2584 0.1371 0.6104 135.0000 45.0000 0 1.0000

ans =

轮槽数 z=6

ans =

销轮转角 槽轮角位移 角速度 角加速度

ans =

-120.0000 -30.0000 0.0000 -0.5774 -110.0000 -29.5432 0.0870 -0.4274 -100.0000 -28.3345 0.1516 -0.3188 -90.0000 -26.5651 0.2000 -0.2400 -80.0000 -24.3737 0.2366 -0.1822 -70.0000 -21.8622 0.2645 -0.1390

23