凸轮机构matlab程序.docx
《凸轮机构matlab程序.docx》由会员分享,可在线阅读,更多相关《凸轮机构matlab程序.docx(8页珍藏版)》请在冰豆网上搜索。
凸轮机构matlab程序
凸轮轮廓程序:
>>e=20;
s0=77.46;
a1=0:
pi/36:
pi/3;
s1=50*[3*a1/pi-sin(6*a1)/(2*pi)];
x1=(s0+s1).*sin(a1)+e*cos(a1);
y1=(s0+s1).*cos(a1)-e*sin(a1);
k1=150/pi*[1-cos(6*a1)];
>>i1=[(k1-e).*sin(a1)+(s0+s1).*cos(a1)].*[(k1-e).*(k1-e)+(s0+s1).*(s0+s1)].^(-1/2);
>>j1=[-(k1-e).*cos(a1)+(s0+s1).*sin(a1)].*[(k1-e).*(k1-e)+(s0+s1).*(s0+s1)].^(-1/2);
>>x10=x1-10*j1;
>>y10=y1-10*i1;
a2=pi/3:
pi/36:
pi;
s2=50;
x2=(s0+s2).*sin(a2)+e*cos(a2);
y2=(s0+s2).*cos(a2)-e*sin(a2);
k2=0;
>>i2=[(k2-e).*sin(a2)+(s0+s2).*cos(a2)].*[(k2-e).*(k2-e)+(s0+s2).*(s0+s2)].^(-1/2);
j2=[-(k2-e).*cos(a2)+(s0+s2).*sin(a2)].*[(k2-e).*(k2-e)+(s0+s2).*(s0+s2)].^(-1/2);
x20=x2-10*j2;
y20=y2-10*i2;
>>a3=pi:
pi/36:
4*pi/3;
>>s3=50*[1-3*(a3-pi)/pi+sin(6*(a3-pi))/(2*pi)];
>>x3=(s0+s3).*sin(a3)+e*cos(a3);
>>y3=(s0+s3).*cos(a3)-e*sin(a3);
>>k3=50*[-3/pi+3/pi*cos(6*(a3-pi))];
>>i3=[(k3-e).*sin(a3)+(s0+s3).*cos(a3)].*[(k3-e).*(k3-e)+(s0+s3).*(s0+s3)].^(-1/2);
>>j3=[-(k3-e).*cos(a3)+(s0+s3).*sin(a3)].*[(k3-e).*(k3-e)+(s0+s3).*(s0+s3)].^(-1/2);
>>x30=x3-10*j3;
>>y30=y3-10*i3;
>>a4=4*pi/3:
pi/36:
2*pi;
>>s4=0;
>>x4=(s0+s4).*sin(a4)+e*cos(a4);
>>y4=(s0+s4).*cos(a4)-e*sin(a4);
>>k4=0;
>>i4=[(k4-e).*sin(a4)+(s0+s4).*cos(a4)].*[(k4-e).*(k4-e)+(s0+s4).*(s0+s4)].^(-1/2);
>>j4=[-(k4-e).*cos(a4)+(s0+s4).*sin(a4)].*[(k4-e).*(k4-e)+(s0+s4).*(s0+s4)].^(-1/2);
>>x40=x4-10*j4;
>>y40=y4-10*i4;
>>plot(x10,y10,'-g*',x20,y20,'-r*',x30,y30,'-b*',x40,y40,'-k*')
>>plot(x10,y10,'-g*',x20,y20,'-r*',x30,y30,'-b*',x40,y40,'-k*',x1,y1,'-g*',x2,y2,'-r*',x3,y3,'-b*',x4,y4,'-k*')
>>
凸轮压力角:
>>x=0:
0.01:
pi/3;
>>t=5.76*sin(6*x)+6*x.*sin(6*x)+(2-2*pi/15)*cos(6*x)-2+2*pi/15;
>>x
x=
Columns1through6
00.01000.02000.03000.04000.0500
Columns7through12
0.06000.07000.08000.09000.10000.1100
Columns13through18
0.12000.13000.14000.15000.16000.1700
Columns19through24
0.18000.19000.20000.21000.22000.2300
Columns25through30
0.24000.25000.26000.27000.28000.2900
Columns31through36
0.30000.31000.32000.33000.34000.3500
Columns37through42
0.36000.37000.38000.39000.40000.4100
Columns43through48
0.42000.43000.44000.45000.46000.4700
Columns49through54
0.48000.49000.50000.51000.52000.5300
Columns55through60
0.54000.55000.56000.57000.58000.5900
Columns61through66
0.60000.61000.62000.63000.64000.6500
Columns67through72
0.66000.67000.68000.69000.70000.7100
Columns73through78
0.72000.73000.74000.75000.76000.7700
Columns79through84
0.78000.79000.80000.81000.82000.8300
Columns85through90
0.84000.85000.86000.87000.88000.8900
Columns91through96
0.90000.91000.92000.93000.94000.9500
Columns97through102
0.96000.97000.98000.99001.00001.0100
Columns103through105
1.02001.03001.0400
>>t
t=
Columns1through6
0.00000.34610.69251.03791.38091.7202
Columns7through12
2.05462.38252.70283.01413.31503.6042
Columns13through18
3.88044.14244.38894.61874.83075.0237
Columns19through24
5.19675.34875.47885.58615.66985.7292
Columns25through30
5.76365.77255.75555.71225.64225.5455
Columns31through36
5.42195.27155.09444.89084.66104.4055
Columns37through42
4.12483.81963.49063.13872.76472.3699
Columns43through48
1.95521.52191.07140.60510.1243-0.3692
Columns49through54
-0.8740-1.3882-1.9102-2.4381-2.9701-3.5044
Columns55through60
-4.0388-4.5716-5.1008-5.6243-6.1401-6.6464
Columns61through66
-7.1410-7.6221-8.0877-8.5359-8.9648-9.3727
Columns67through72
-9.7578-10.1184-10.4528-10.7596-11.0372-11.2843
Columns73through78
-11.4996-11.6820-11.8304-11.9439-12.0215-12.0627
Columns79through84
-12.0668-12.0335-11.9623-11.8531-11.7059-11.5207
Columns85through90
-11.2979-11.0378-10.7410-10.4081-10.0399-9.6374
Columns91through96
-9.2016-8.7339-8.2356-7.7081-7.1530-6.5722
Columns97through102
-5.9674-5.3405-4.6937-4.0291-3.3489-2.6554
Columns103through105
-1.9511-1.2383-0.5195
>>x=0.46:
0.001:
0.47;
>>t=5.76*sin(6*x)+6*x.*sin(6*x)+(2-2*pi/15)*cos(6*x)-2+2*pi/15;
>>x
x=
Columns1through6
0.46000.46100.46200.46300.46400.4650
Columns7through11
0.46600.46700.46800.46900.4700
>>t
t=
Columns1through6
0.12430.07550.0266-0.0225-0.0717-0.1210
Columns7through11
-0.1704-0.2199-0.2696-0.3193-0.3692
>>x=0.462:
0.0001:
0.463;
>>t=5.76*sin(6*x)+6*x.*sin(6*x)+(2-2*pi/15)*cos(6*x)-2+2*pi/15;
>>x
x=
Columns1through6
0.46200.46210.46220.46230.46240.4625
Columns7through11
0.46260.46270.46280.46290.4630
>>t
t=
Columns1through6
0.02660.02170.01680.01190.00700.0021
Columns7through11
-0.0028-0.0077-0.0126-0.0176-0.0225
>>z=[00.4625];
>>y=[150/pi*(1-cos(6*z))-20].*[45.8+50*(3*z/pi-sin(6*z)/(2*pi))];
>>y=[150/pi*(1-cos(6*z))-20]./[45.8+50*(3*z/pi-sin(6*z)/(2*pi))];
>>y
y=
-0.43671.1121
>>a=atan(y)
a=
-0.41170.8384
>>a*180/pi
ans=
-23.590048.0382
>>
以上是计算推程压力角的临界值。
回程与推程相同,故临界压力角也相同。
近休止区临界压力角为23.6度;远休止区为11.8度。
所以整个凸轮运动过程中的临界压力角为48度。
许用压力角为50度(该凸轮为封闭力系)。
故满足要求。
(注:
可编辑下载,若有不当之处,请指正,谢谢!
)