机械原理大作业.docVIP

机械原理作业 C程序代码 #include stdio.h #include math.h #define pi 3.14159 void main() { float a=0.1,b=0.3,e=0.03,w1=10; float m,t,L,angle,d,t1,t2,w2,a2,a3,v3; printf(please input jiaodu(example 30 °shuru 30）:\n); scanf(%f,d); d=d/180*pi; m=e-a*sin(d); t=m/b; L=asin(t); angle=L*180/pi; if(angle0) angle=angle+360; t1=a*w1*cos(d); t2=b*cos(L)*(-1); w2=t1/t2; a2=((a*w1*w1*sin(d)+b*w2*w2*sin(L))/(b*cos(L)))*(-1); v3=((a*w1*sin(d-L))/cos(L))*(-1); a3=((a*w1*w1*cos(L-d)+b*w2*w2)/cos(L))*(-1); printf(\n*****************************\n); printf(\nJiaoweiyi2=%.3f° \n\nJiaoshudu2=%.3f rad/s\n,angle,w2); printf(\nJiaojaishudu2=%.3f rad/s2\n\nshudu3=%.3f m/s\n\njiashudu3=%.3f m/s2\n,a2,v3,a3); printf(\n*****************************\n); } 运行结果图： 当输入的角度为30°时的结果如下： 当输入的角度为60°时的结果如下： 当输入的角度为120°时的结果如下： MATLAB图像编程源程序： clear; a=0.1;b=0.3;e=0.03;w1=10; d=0:360; l=d/180*pi; L= asin((e-a.*sin(l))/b); g=L*180/pi; w=a*w1*cos(l)./(b.*cos( asin((e-a.*sin(l))/b))*(-1)); t1=a.*w1.^2.*sin(d/180*pi); t2=b.*w.^2.*sin(asin((e-a.*sin(l))./b)); t=t1+t2; j=(b.*cos(asin((e-a.*sin(l))./b)))*(-1); a2=t./j; v3=((a*w1.*sin(d./180*pi-asin((e-a.*sin(l))/b)))./cos(asin((e-a.*sin(l))/b)))*(-1); a3=((a.*w1.^2.*cos(asin((e-a.*sin(l))./b)-l)+b.*w.^2)./cos(asin((e-a.*sin(l))/b)))*(-1); figure(1) subplot(5,1,1); plot(d,g); grid on axis on title(角位移θ2) xlabel(角度);ylabel(角位移θ2); title(角度与角位移θ2图像); subplot(5,1,2) plot(d,w);title(构件2的角速度w2图像);xlabel(角度);ylabel(角速度w2);grid on; subplot(5,1,3) plot(d,a2);title(构件2的角加速度a2图像);xlabel(角度);ylabel(角加速度a2);grid on; subplot(5,1,4) plot(d,v3);title(构件3的速度v3图像);xlabel(角度);ylabel(速度v3);grid on; subplot(5,1,5) plot(d,a3);title(构件3的加速度v3图像);xlabel(角度);ylabel(加速度v3);grid on; 所的的图像MATLAB软件显示：