平面連桿機(jī)構(gòu)設(shè)計(jì)方案-文庫(kù)吧
2025-04-10 01:08 本頁(yè)面
【正文】 0至360度，步長(zhǎng)為pi/6th34=zeros(length(th2),2)。 %建立一個(gè)N行2列的零矩陣，第一列存放options=optimset(39。display39。,39。off39。)。 %θ_3,第二列存放θ_3for m=1:length(th2) %建立for循環(huán)，求解θ_3，θ_4th34(m,:)=fsolve(39。fourbarposition39。,[1 1],… %調(diào)用fsove函數(shù)求解關(guān)于θ_3，θ_4options,th2(m),L2,L3,L4,L1)。 %的非線性超越方程，結(jié)果保存在th34中endy=L2*sin(th2)+L3*sin(th34(:,1)39。)。 %連桿3的D端點(diǎn)Y坐標(biāo)值x=L2*cos(th2)+L3*cos(th34(:,1)39。)。 %連桿3的D端點(diǎn)X坐標(biāo)值xx=[L2*cos(th2)]。 %連桿3的C端點(diǎn)X坐標(biāo)值yy=[L2*sin(th2)]。 %連桿3的C端點(diǎn)Y坐標(biāo)值figure(1)plot([x。xx],[y。yy],39。k39。,[0 L1],[0 0],… %繪制連桿3的幾個(gè)位置點(diǎn)39。k^39。,x,y,39。ko39。,xx,yy,39。ks39。)title(39。連桿3的幾個(gè)位置點(diǎn)39。)xlabel(39。水平方向39。)ylabel(39。垂直方向39。)axis equal%XY坐標(biāo)均衡th2=[0:2/72:2]*pi。 %重新細(xì)分曲柄輸入角度θ_2，步長(zhǎng)為5度th34=zeros(length(th2),2)。options=optimset(39。display39。,39。off39。)。for m=1:length(th2)th34(m,:)=fsolve(39。fourbarposition39。,[1 1],…options,th2(m),L2,L3,L4,L1)。endfigure(2)plot(th2*180/pi,th34(:,1),th2*180/pi,th34(:,2)) %繪制連桿3的角位移關(guān)于曲柄2的角位移圖plot(th2*180/pi,th34(:,1)*180/pi,…th2*180/pi,th34(:,2)*180/pi) %繪制搖桿4的角位移關(guān)于曲柄2的角位移圖axis([0 360 0 170])%確定XY邊界值grid%圖形加網(wǎng)格xlabel(39。主動(dòng)件轉(zhuǎn)角\theta_2(度)39。)ylabel(39。從動(dòng)件角位移(度)39。)title(39。角位移線圖39。)text(120,120,39。搖桿4角位移39。)text(150,40,39。連桿3角位移39。)w2= %設(shè)定曲柄角速度f(wàn)or i=1:length(th2)A=[L3*sin(th34(i,1)) L4*sin(th34(i,2))?！璍3*cos(th34(i,1)) L4*cos(th34(i,2))]。B=[w2*L2*sin(th2(i))。 w2*L2*cos(th2(i))]。w=inv(A)*B。w3(i)=w(1)。w4(i)=w(2)。endfigure(3)plot(th2*180/pi,w3,th2*180/pi,w4)。 %繪制角速度線圖axis([0 360 175 200]) text(50,160,39。搖桿4角速度(\omega_4)39。)text(220,130,39。連桿3角速度(\omega_3)39。)gridxlabel(39。主動(dòng)件轉(zhuǎn)角\theta_2(度)39。)ylabel(39。從動(dòng)件角速度(rad\cdot s^{1})39。)title(39。角速度線圖39。)for i=1:length(th2)C=[L3*sin(th34(i,1)) L4*sin(th34(i,2))?！璍3*cos(th34(i,1)) L4*cos(th34(i,2))]。D=[w2^2*L2*cos(th2(i))+w3(i)^2*L3*cos(th34(i,1))w4(i)^2*L4*cos(th34(i,2))。...w2^2*L2*sin(th2(i))+w3(i)^2*L3*sin(th34(i,1))w4(i)^2*L4*sin(th34(i,2))]。a=inv(C)*D。a3(i)=a(1)。a4(i)=a(2)。endfigure(4)plot(th2*180/pi,a3,th2*180/pi,a4)。 %繪制角加速度線圖axis([0 360 70000 65000])text(50,50000,39。搖桿4角加速度(\alpha_4)39。)text(220,12000,39。連桿3角加速度(\alpha_3)39。)gridxlabel(39。從動(dòng)件角加速度39。)ylabel(39。從動(dòng)件角加速度(rad\cdot s^{2})39。)title(39。角加速度線圖39。)disp 39。曲柄轉(zhuǎn)角連桿轉(zhuǎn)角搖桿轉(zhuǎn)角連桿角速度搖桿角速度連桿加速度搖桿加速度39。ydcs=[th239。*180/pi,th34(:,1)*180/pi,th34(:,2)*180/pi,w339。,w439。,a339。,a439。]。disp(ydcs) 程序運(yùn)行結(jié)果輸出 * * * * * * 平面四桿機(jī)構(gòu)的運(yùn)動(dòng)分析 * * * * * *曲柄轉(zhuǎn)角 連桿轉(zhuǎn)角搖桿轉(zhuǎn)角連桿角速度搖桿角速度連桿加速度搖桿加速度 +004 *16
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1