齿轮论坛 www.gearbbs.net

 找回密码
 注册

手机号码,快捷登录

手机号码,快捷登录

搜索
查看: 3028|回复: 4
收起左侧

Matlab 编程摆线轮不知哪里错了,望高手止点

[复制链接]
发表于 2016-9-8 15:51 | 显示全部楼层 |阅读模式

本帖子中包含更多资源。

您需要 登录 才可以下载或查看,没有账号?注册

x
程序如下, 而且没有生成一个完整的齿圈。

%% parameters
rp=34;      %  针齿中心圆半径

rrp=11.1;                %  砂轮圆弧半径 (针齿半径)

a=3.8;                    %  偏心距

zp=7;                     %  针轮齿数

delta_rp=-0.043;          % 移距修形
delta_rrp=-0.033;         % 等距修形

delta=0;                  % 转角修形

iH= 1;                    % 摆线轮和针轮的相对传动比

phi=[0:1:180];            % degree

%%
intermediate parameters
K1=a*zp/(rp+delta_rp);
for j=1:length(phi)
Sr=1+K1^2-2*K1*cos(phi(j)*pi/180);   

  
xc(j)=(rp+delta_rp-(rrp+delta_rrp)*Sr^(-1/2)) ...  
     *cos(((1-iH)*phi(j)-delta)*pi/180)-a/(rp+delta_rp) ...
     *(rp+delta_rp-zp*(rrp+delta_rrp)*Sr^(-1/2))*cos((iH*phi(j)+delta)*pi/180);
yc(j)=(rp+delta_rp-(rrp+delta_rrp)*Sr^(-1/2)) ...     
     *sin(((1-iH)*phi(j)-delta)*pi/180)+a/(rp+delta_rp) ...
     *(rp+delta_rp-zp*(rrp+delta_rrp)*Sr^(-1/2))*sin((iH*phi(j)+delta)*pi/180);

end
plot(xc,yc,'Linewidth',2)
title('shape reform')
xlabel('x_c')
ylabel('y_c')

data=[xc' yc'];

save e:\kala_xy.txt data -ascii

发表于 2016-10-7 06:34 | 显示全部楼层
参数有问题吧
发表于 2016-12-17 20:36 | 显示全部楼层
plot函数,乘除法使用.*和./代替*和/。试试吧
发表于 2017-3-27 21:04 | 显示全部楼层
公式错了 ...
发表于 2017-12-28 13:37 | 显示全部楼层
语句有问题
您需要登录后才可以回帖 登录 | 注册

本版积分规则

QQ|小黑屋|手机版|Archiver|齿轮论坛

GMT+8, 2024-4-27 12:25 , Processed in 0.161786 second(s), 9 queries , MemCache On.

Powered by Discuz! X3.4 Licensed

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表