?? 03-連桿機構(gòu)的運動精度綜合-1.m__.htm
字號:
<html><head><title>MATLAB.rar 03-連桿機構(gòu)的運動精度綜合-1.m</title>
<LINK href="/inc/read_style.css" type=text/css rel=stylesheet></head>
<body>
<p><a href=http://www.pudn.com>www.pudn.com</a> > <a href="http://www.pudn.com/downloads111/sourcecode/graph/detail461284.html">MATLAB.rar</a> > 03-連桿機構(gòu)的運動精度綜合-1.m</p><!-- saved from http://www.pudn.com -->
<script src="/inc/gg_read1.js"></script><BR>
<pre name="code" class="m">
disp ' ****** 曲柄滑塊機構(gòu)的運動學(xué)分析 *******'
SC=input(' 輸入滑塊行程的均值(mm) SC = ');
P=input(' 輸入曲柄軸心至滑銷最遠(yuǎn)距離(mm) P = ');
E=input(' 輸入機構(gòu)偏心距的均值(mm) E = ');
RL=input(' 輸入曲柄與連桿長度比的均值 RL = ');
DR=input(' 輸入曲柄長度偏差(mm) DR = ');
DL=input(' 輸入連桿長度偏差(mm) DL = ');
DE=input(' 輸入機構(gòu)偏心距偏差(mm) DE = ');
N=input(' 輸入曲柄轉(zhuǎn)速(r/min) N = ');
L=sqrt((P-SC)^2-E^2)/(1-RL);
fprintf(1,' 連桿長度的均值(mm) L = %3.6f \n',L)
R=RL*L;
fprintf(1,' 曲柄長度的均值(mm) R = %3.6f \n',R)
CR=DR/3;CL=DL/3;CE=DE/3;
EL=E/L;
fprintf(1,' 偏心距與連桿長度比的均值(mm) EL = %3.6f \n',EL)
fprintf(1,' 曲柄長度的標(biāo)準(zhǔn)離差(mm) CR = %3.6f \n',CR)
fprintf(1,' 連桿長度的標(biāo)準(zhǔn)離差(mm) CR = %3.6f \n',CL)
fprintf(1,' 偏心距的標(biāo)準(zhǔn)離差(mm) CE = %3.6f \n',CE)
W=pi*N/30;
fprintf(1,' 曲柄的角速度(mm) W = %3.6f \n',W)
CRL=sqrt((R*CL)^2+(L*CR)^2)/L^2;
fprintf(1,' 曲柄與連桿長度比的標(biāo)準(zhǔn)離差 CRL = %3.6f \n',CRL)
CEL=sqrt((E*CL)^2+(L*CE)^2)/L^2;
fprintf(1,' 偏心距與連桿長度比的標(biāo)準(zhǔn)離差 CEL = %3.6f \n',CEL)
theta=0:10:360;
hd=theta.*pi/180;
% 計算滑塊位移、速度、加速度的均值
S=R.*(1-cos(hd)-EL.*sin(hd)+0.5.*RL.*sin(hd).^2);
V=R.*W.*(sin(hd)-EL.*cos(hd)+0.5.*RL.*sin(2.*hd));
A=R.*W^2.*(cos(hd)+EL.*sin(hd)+RL.*cos(2.*hd));
figure(1);
subplot(1,3,1);
plot(theta,S,'r')
title('\bf \mus 線圖')
subplot(1,3,2);
plot(theta,V,'k')
title('\bf \muv 線圖')
xlabel('\bf 曲柄轉(zhuǎn)角\theta(度)')
subplot(1,3,3);
plot(theta,A,'b')
title('\bf \mua 線圖')
% 計算滑塊位移、速度、加速度的標(biāo)準(zhǔn)離差
CS=sqrt((1-cos(hd)+(0.5.*RL.*sin(hd)-EL).*sin(hd)).^2.*CR^2+(0.5.*(CRL.*sin(hd)).^2-CEL^2).^2.*(R.*sin(hd)).^2);
CV=W.*sqrt((sin(hd)-EL.*cos(hd)+0.5.*RL.*sin(2.*hd)).^2.*CR^2+(0.5.*(CRL.*sin(hd)).^2-(CEL.*cos(hd)).^2)*R^2);
CA=W^2.*sqrt((cos(hd)+EL.*sin(hd)+RL.*cos(2.*hd)).^2.*CR^2+(R.*CRL.*sin(hd)).^2+(RL.*CEL.*cos(2.*hd)).^2);
figure(2);
subplot(1,3,1);
plot(theta,CS,'r')
title('\bf \sigmas 線圖')
subplot(1,3,2);
plot(theta,CV,'k')
title('\bf \sigmav 線圖')
xlabel('\bf 曲柄轉(zhuǎn)角\theta(度)')
subplot(1,3,3);
plot(theta,CA,'b')
title('\bf \sigmaa 線圖')
% 計算滑塊位移、速度、加速度的偏差
DS=3.*CS;DV=3.*CV;DA=3.*CA;
% 計算滑塊位移、速度、加速度的最大值和最小值
SM=S+DS;SN=S-DS;
VM=V+DV;VN=V-DV;
AM=A+DA;AN=A-DA;
% 計算滑塊位移、速度、加速度的差值
SD=2.*DS;VD=2.*DV;AD=2.*DA;
figure(3);
subplot(1,3,1);
plot(theta,SD,'r')
title('\bf \Deltas 線圖')
subplot(1,3,2);
plot(theta,VD,'k')
title('\bf \Deltav 線圖')
xlabel('\bf 曲柄轉(zhuǎn)角\theta(度)')
subplot(1,3,3);
plot(theta,AD,'b')
title('\bf \Deltaa 線圖')
</pre>
<script src="/inc/gg_read2.js"></script><BR>
<script src="http://s117.cnzz.com/stat.php?id=1236358&web_id=1236358&show=pic" language="JavaScript" charset="gb2312"></script>
</body></html>
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -