西工大机械原理大作业报告Word下载.docx
《西工大机械原理大作业报告Word下载.docx》由会员分享，可在线阅读，更多相关《西工大机械原理大作业报告Word下载.docx（34页珍藏版）》请在冰豆网上搜索。
（1）推杆运动规律
匀速时：
远休：
等加速：
等减速：
近休：
（2）轮廓线方程
（a）理论廓线
（b）内包络线
（c）外包络线
3.计算程序
%------------------------------绘制凸轮轮廓曲线------------------------------
title（'
绘制凸轮轮廓曲线'
）
holdon
%-----------凸轮部件几何数据输入--------
e=10;
%偏心距
rb=35;
%基圆半径
rr=15;
%滚子半径
%--------------运动规律输入-------------
a1=150;
%等速截止角
a2=180;
%远休截止角
a3=300;
%匀加匀减截止角
h=15;
%行程
%--------------计算运动规律-------------
s0=sqrt（rb^2-e^2）;
aymax=0;
%最大压力角
ymaxa=0;
%最大压力角转过角度
roumin=10000;
%最小曲率半径
roumina=0;
%最小曲率半径对应角度
ay30=0;
%30度对应压力角
dzh=pi/180;
%角度转弧度
a3=（a3-a2）/2;
a2=a2-a1;
vmax=h/a3;
fora=0.0001:
0.0001:
360
f=round（10000*a）;
ifa<
=a1%等速
s（f）=h*a/a1;
dsa（f）=h/a1;
dssaa（f）=0;
elseifa>
a1&
a<
=（a1+a2）%远休
s（f）=h;
dsa（f）=0;
（a1+a2）&
=（a1+a2+a3）%回程等加速
dsa（f）=-（a-a1-a2）/a3*vmax;
s（f）=h+dsa（f）*（a-a1-a2）/2;
dssaa（f）=-vmax/a3;
（a1+a2+a3）&
a<
=（a1+a2+2*a3）%回程等减速
dsa（f）=-（a1+a2+2*a3-a）/a3*vmax;
s（f）=-（a1+a2+2*a3-a）*dsa（f）/2;
dssaa（f）=vmax/a3;
（a1+a2+2*a3）&
=360%近休
s（f）=0;
dsa（f）=0;
end
%-----------计算轮廓曲线---------
xb（f）=-（s（f）+s0）*sin（a*dzh）-e*cos（a*dzh）;
%理论廓线
yb（f）=（s（f）+s0）*cos（a*dzh）-e*sin（a*dzh）;
dxb（f）=-（s0+s（f））*cos（a*dzh）-（dsa（f）-e）*sin（a*dzh）;
dyb（f）=-（s0+s（f））*sin（a*dzh）+（dsa（f）-e）*cos（a*dzh）;
dxxbb（f）=（s0+s（f）-dssaa（f））*sin（a*dzh）-（2*dsa（f）-e）*cos（a*dzh）;
dyybb（f）=-（s0+s（f）-dssaa（f））*cos（a*dzh）-（2*dsa（f）-e）*sin（a*dzh）;
x（f）=xb（f）-rr*dyb（f）/sqrt（dxb（f）*dxb（f）+dyb（f）*dyb（f））;
%内包络线
y（f）=yb（f）+rr*dxb（f）/sqrt（dxb（f）*dxb（f）+dyb（f）*dyb（f））;
xw（f）=xb（f）+rr*dyb（f）/sqrt（dxb（f）*dxb（f）+dyb（f）*dyb（f））;
%外包络线
yw（f）=yb（f）-rr*dxb（f）/sqrt（dxb（f）*dxb（f）+dyb（f）*dyb（f））;
%-----------计算曲率半径-----------
rou（f）=（（dxb（f）*dxb（f）+dyb（f）*dyb（f））^（3/2））/（dxb（f）*dyybb（f）-dyb（f）*dxxbb（f））;
ifrou（f）<
roumin
roumin=rou（f）;
roumina=a;
%-----------计算压力角------------
ay（f）=atan（（dsa（f）-e）/（sqrt（rb*rb-e*e）+s（f）））;
ifay（f）<
aymax
aymax=-ay（f）/dzh;
ymaxa=a;
%----------查看各点坐标---------------
a30=-ay（30000）/dzh;
fori=1:
1:
360%用zb查看实际曲线各个点坐标
zb（i,1）=x（i*10000）;
zb（i,2）=y（i*10000）;
360%用zbT查看理论曲线各个点坐标
zbT（i,1）=xb（i*10000）;
zbT（i,2）=yb（i*10000）;
%------------------画图------------------
plot（x,y,'
-'
'
Color'
[100],'
LineWidth'
2）;
%内包络线，红色
plot（xb,yb,'
[010],'
1）%理论廓线，绿色
plot（xw,yw,'
[001],'
2）%外包络线，蓝色
alpha=0:
pi/50:
2*pi;
xj=rb*cos（alpha）;
yj=rb*sin（alpha）;
plot（xj,yj,'
.'
[0.80.30.9],'
0.5）%基圆
xg=-e+rr*cos（alpha）;
yg=s0+rr*sin（alpha）;
plot（xg,yg,'
[000],'
1.5）%滚子
gx=[-e,-e];
gy=[s0s0+50];
plot（gx,gy,'
1.5）%推杆
xp=e*cos（alpha）;
yp=e*sin（alpha）;
plot（xp,yp,'
:
'
[0.90.60.3]）%偏心圆
axis（[-120120-8585]）
axisequal
text（31,106,'
30°
角对应压力角='
[000]）
text（92,106,num2str（a30）,'
text（50,99,'
最大压力角='
text（92,99,num2str（aymax）,'
text（20,92,'
最大压力角对应转角='
text（92,92,num2str（ymaxa）,'
text（35,85,'
最大小曲率半径='
text（92,85,num2str（roumin）,'
text（13,78,'
最小曲率半径对应转角='
text（92,78,num2str（roumina）,'
text（60,68,'
内包络线'
[100]）
text（60,61,'
外包络线'
[001]）
text（60,54,'
理论轮廓线'
[010]）
text（60,47,'
滚子轮廓线'
text（60,40,'
偏心距圆轮廓线'
[0.90.60.3]）
text（60,33,'
基圆轮廓线'
[0.80.30.9]）
text（-2,0,'
O'
4.计算结果及分析
（x,y）坐标如下图，本程序虽然在画图时取了3,600,000个点进行计算，但是为了表示结果，在程序中控制每10000个点输出一个，即1°
一个点，共360个点。
（程序的67-72行）
其他计算结果显示在figure1中，见“5.凸轮理论廓线及实际廓线图”。
以下计算结果均为h=15的结果
理论轮廓线坐标（分了四栏）：
-10.585593699678533.4613719078855
-11.171452874641333.3714705779453
-11.757397471818833.2712821587502
-12.343246349534033.1607763070274
-12.928817333490033.0399258779396
-13.513927273394332.9087069518753
-14.098392100204532.7670988602346
-14.682026883976032.6150842101964
-15.264645892294132.4526489084537
-15.846062649272332.2797821839047
-16.426089995096832.0964766092865
-17.004540146100831.9027281217397
-17.581224755346331.6985360422921
-18.155********6531.4839030942496
-18.728541511356831.2588354204836
-19.2987946998670.023*********
-19.866524554522030.7774376610052
-20.431540837202230.5211370987833
-20.993653119592130.2544608844998
-21.552670846766729.9774324787985
-22.108403401124729.6900788418560
-22.660660166647029.3924304426620
-23.209250593458029.0845212671200
-23.7539842626713
升级会员 