第四章计算实例
4 1对外槽轮机构不同的槽数进行分析
外槽轮机构的槽数2-4 6 8 10时计算从动件槽轮机构的运动 参数日，u，：，并且绘制槽轮的运动线图（粪角速度\／n、和类角加速 度；／”）。
4 1 1外槽轮机构槽数不同时编写的M文件 外槽轮机构运动分析的M文件和运行结果 >> %外槽轮机构运动分析
Dr=pi/180.0;%角度与弧度的转换系数
%销轮2转角范围:-f20 lmd=sin(pi/z);%计算曲柄2与机架1的长度此 bc=5;%循环步长 cz=-f20/dr;%循环初值 zz=f20/dr;%循环终值 i=1;%根据步长变化的运动参数矩阵cs行数计数器 for f2=cz:bc:zz %计算槽轮角位移,类角速度,类角加速度 wy=atan(lmd*sin(f2*d r）/（1-lmd'*cos(f2*d r)））． sd=lmd*(cos(f2*d r)-lmd)/(1 -2*lmd*cos(f2*d r)+lmd^2); jsd=-lmd*sin(f2*dr)*(1-lmd^2)/(1-2*lmd*cos(f2*dr)+lmd^2)^2; switch z %矩阵c(i,:)表示第i行的各列元素 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]; End i=i+1; End End %输出外槽轮机构运动参数 [＇轮槽数 z=4＇] [＇销轮转角＇，＇槽轮角位移＇，＇角速度＇，＇角加速度＇] % 矩阵c(:,j)表示第j列的各行元素 [c8(:,1),c8(:,2),c8(:,4)] [＇轮槽数 z=10＇] [＇销轮转角＇,＇槽轮角位移＇,＇角速度＇,＇角加速度＇] % 矩阵c(:,j)表示第j列的各行元素 [c10(:,1),c10(:,2),c10(:,3),c10(:,4)] % %绘制槽轮机构运动参数曲线 Figure(1); %生成槽轮运动线图窗口 Subplot(2,2,1); %选择第一个子窗口 Plot(c4(:,1),c4(:,3),c4(:,1),c4(:,4)) %绘制z=4的窗口 Title(外槽轮槽数 z=4’) %标注子窗口的名称 Axis([-pi/4/dr pi/4/dr -6 6]) %定义坐标轴范围 Grid %栅格线 Text(-2,4,2,’\\epsilon/\\omega^{2}’) %标注类角加速度线图 Text(20,1.6,’\\omega/\\omega’) %标注类角速度线图 Ylabel(‘槽轮运动线图’) %定义坐标轴名称 % Subplot(2,2,2); %选择第2个子窗口 Plot(c6(:,1),c6(:,3),c6(:,1),c6(:,4)) %绘制z=6的线图 Title(‘外槽轮槽数 z=6’) Axis([-pi/3/dr pi/3/dr -1.5 1.5]) Grid Text(10,-0.7,’\\epsilon/\\omega^{2}’) Text(30,0.7,’\\omega/\\omega’) Ylabel(‘槽轮运动线图’) % Subplot(2,2,3); %选择第3个子窗口 Plot(c8(:,1),c8(:,3),c8(:,1),c8(:,4)) %绘制z=8的线图 Title(‘外槽轮槽数 z=8’) Axis([-3*pi/8/dr 3*pi/8/dr -0.8 0.8]) Grid Text(12,-0.3,’\\epsilon/\\omega^{2}’) Text(40,0.4,’\\omega/\\omega’) Ylabel(‘槽轮运动线图’) % Subplot(2,2,4); %选择第4个子窗口 Plot(c10(:,1),c10(:,3),c10(:,1)(:,4)%绘制z=10的线图 Title(‘外槽轮槽数 z=10’) Axis([-2*pi/5/dr 3*pi/5dr -0.5 0.5]) Grid Text(15,-0.2,’\\epsilon/\\omega^{2}’) Text（40,0.3’\\omega/\\omega’） Ylabel(‘槽轮运动线图’) % Figure(2); %生成类线图窗口 Subplot(1,2,1); %选择第一个子窗口 Plot(c4(:,1),c4(:,3),c6(:,3),c8(:,1),c8(:,3),c10(:,1),c10(:,3) Title(‘\\omega/\\omega’) Axis([-f20/dr 20/dr -0.1 2.5]) Grid Text(-10,0.35,’z=10’) Text(-8,0.7,’z=8’) Text(-8,1.1’,z=6’) Text(10,2.1,’z=4’) Ylabel(‘槽轮类角速度线图’) % Subplot(1,2,2); %选择第二个子窗口 Plot(c4(:,1),c4(:,4),c6(:,4),c8(:,1),c8(:,4),c10(:,1),c10(:,4) Title(‘\\epsilon/\\omega^{2}’) Axis([-f20/dr f20/dr -5.5 5.5]) Grid Text(-50,0.2,’z=10’) Text(-30,0.9,’z=8’) Text(-25,1.6,’z=6’) Text(0,3.5,’z=4) Ylabel(槽轮类角加速度线图’) %