平面四桿機構的設計與運動分析Matlab代碼

1、平面四桿機構的設計與運動分析M代碼平面四桿機構的設計M代碼-A=cos(50*pi/180),cos(50-35)*pi/180),1;cos(75*pi/180),cos(75-80)*pi/180),1;cos(105*pi/180),cos(105-125)*pi/180),1;B=cos(35*pi/180);cos(80*pi/180);cos(125*pi/180);P=ABm= P(1)n=-m/ P(2)l=sqrt(mA2+nA2+1-2*n*P(3)運行設計結果顯示：Command Window A= cos r：50*pi/180) cos (50-35) #pi/180

2、) j 1; cos r75*pi/180) cos (75-B=cos f35*pi/180):cos (80*pi/180);cosf125*pi/180):P=ABm二 P (1)n=-m/ P(2)l=sqrt1m 2+n 2+l-2*n* P(3)P =1.5815-1.26401.0235m 二1.5815n =1.25121 =1.5829平面四桿機構的運動分析M代碼-%參數(shù)賦值clc,clear10=1.2512;11=1.0;12=1.5829;13=1.5815;M=-1;Omiga1=10;Theta1=0:0.01:360;Theta1=Theta1*pi/180;%求

3、解各個構件位移、速度、加速度A=2*11*12*sin(Theta1);B=2*12*(11*cos(Theta1)-10);C=11A2+12A2+10A2-13A2-2*11*10*cos(Theta1);E=2*11*13*sin(Theta1);F=2*13*(11*cos(Theta1)-10);G=12A2-11A2-13A2-10A2+2*11*10*cos(Theta1);Theta3=2*atan(E+M*sqrt(E.A2+ F2- GA2)./(F-G);Theta2=2*atan(A+M*sqrt(A.A2+B.A2-C.A2)./(B-C);Omiga2=Omiga1

4、*11*sin(Theta1-Theta3)./(12*sin(Theta3-Theta2);Omiga3=Omiga1*11*sin(Theta1-Theta2)./(13*sin(Theta3-Theta2);A1fa3=(Omiga1A2*11*cos(Theta1-Theta2)+Omiga2.A2*12-Omiga3.A2*13.*cos(Theta3-Theta2)./ (13*sin(Theta3-Theta2);A1fa2=(-Omiga1A2*11*cos(Theta1-Theta3)+Omiga3.A2*13-Omiga2.A2*12.*cos(Theta2-Theta3)

5、./ (12*sin(Theta2-Theta3);%繪圖Theta1=Theta1*180/pi;Subp1ot(3,1,1)p1ot(Theta1,Theta3*180/pi),grid onx1abe1(曲柄轉角(人。);y1abe1(CD桿角位移(人。);Subp1ot(3,1,2)p1ot(Theta1,Omiga3),grid onx1abe1(曲柄轉角(人。);y1abe1(CD 桿角速度(rad/s);Subp1ot(3,1,3)p1ot(Theta1,A1fa3),grid onx1abe1(曲柄轉角(人。);y1abe1(CD 桿角加速度(rad/sA2);by Xu jianpingFS-Ed_i Inseirt Tclos- Dmsktop w-ndclw Hfp-PCD桿角位移。）CDGfex_l_一口口.口i-2CD OCD CD CDg_ks15Q 愣口 2S3S350 4S soCD桿角速度（ad/s）cnunCD CD CDC