- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
直动从动件盘形凹轮机构的计算机辅助设计
直动从动件盘形凹轮机构的计算机辅助设计
—— Matlab 实现
雷昊
航天学院
探测制导与控制技术
1080420103
涡轮机构的计算机辅助程序
题目:
用解析法求出涡轮的轮廓线的方程,精确计算出轮廓各点的坐标,用计算机软件把轮廓线绘制出来。
二、设计思路
由于几乎本人用过的所有c的编译器对win7都不兼容,只能选用matlab来实现程序。
设计思路:
程序开始,由操作者输入偏心距E,基圆半径R,最大位移H,以及凹轮的转角范围J1、J2、J3、J4。程序读入输入的变量。
首先设置for循环,递增变量J的值判断J的值位于的区间,从而确定S的函数。(S1=(H/2)*(1-cos(pi*J/J1)); S2=H; S3=H-(H/2*(1-cos(pi*J/J3))); S4=0;)。
根据x=E*cos(J)+(sqrt(R*R-E*E)+S)*sin(J);
y=(sqrt(R*R-E*E)+S)*cos(J)-E*sin(J)来计算在不同角度J下,X,Y 的坐标值,运用打点函数plot(),来实现绘制函数图像。
由于题目要求要设置旋转的图像,所以,为实现此功能,还必须再在程序中设置一个for循环,选取I为变量,使I从0度以5度为步长增加至360度,每一部都绘制出一个图形,只不过,对应的角度不一样,在新图像生成时候,使用clf或者hold off来清除上一个图像,在用延时函数pause()来使图像暂留, 这样就有了图形旋转的效果。
三、程序:
实现s随着J的变化函数
以及s的导数随J的变换函数的程序
R=input(请输入基圆半径R的值:);
E=input(请输入偏心距E的值:);
H=input(请输入最大位移H的值:);
J1=input(请输入推程J1的值:);
J2=input(请输入远休止J2的值:);
J3=input(请输入回程J3的值:);
J4=input(请输入近休止J4的值:);
hold on
for J=0:pi/60:2*pi
S1=(H/2)*(1-cos(pi*J/J1));
SS1=((pi*H)*sin(pi*J/J1))/(2*J1);
S2=H;
SS2=0;
S3=H-(H/2*(1-cos(pi*J/J3)));
SS3=-(pi*H*sin(pi*J/J3))/(2*J3);
S4=0;
SS4=0;
if J(2*pi-J4)
S=S4;
SS=SS4;
elseif J(2*pi-J4-J3)J(2*pi-J4)
S=S3;
SS=SS3;
elseif J(2*pi-J4-J3-J2)J(2*pi-J4-J3)
S=S2;
SS=SS2;
else
S=S1;
SS=SS1;
end
plot(J,S);
plot(J,SS,b:*);
end
axis square
实现凸轮旋转的程序:
R=input(请输入基圆半径R的值:);
E=input(请输入偏心距E的值:);
H=input(请输入最大位移H的值:);
J1=input(请输入推程J1的值:);
J2=input(请输入远休止J2的值:);
J3=input(请输入回程J3的值:);
J4=input(请输入近休止J4的值:);
J=0:pi/10:2*pi;
x1=R*cos(J);
y1=R*sin(J);
s0=sqrt(R^2-E^2);
for i=1:360
if i=270
s(i)=0;
elseif i=180
s(i)=H-H*(1-cos(pi*i/J1))/2;
elseif i=90
s(i)=H;v(i)=0;a(i)=0;
else
s(i)=H*(1-cos(pi*i/J3))/2;
end
x(i)=E*cos(pi*i/180)+(s0+s(i))*sin(pi*i/180);
y(i)=(s0+s
文档评论(0)