返回   MATLAB中国论坛|MATLAB爱好者之家—不仅仅是MATLAB! > 基础科学-不仅仅是MATLAB! > 数学


数学 A discussion board for Mathematics.



回复
 
LinkBack 主题工具 显示模式
旧 2013-05-30, 05:24 PM   #1
初级会员
 
注册日期: 2013-05-30
帖子: 1
感谢他人: 0
有 0 帖获得 0 感谢
声望力: 0
danxiaojuan 正向着好的方向发展
默认 求大神看看我的程序哪里不对

小女菜鸟一枚,只懂皮毛。这是根据立铣刀铣削力模型编出来的程序,要出铣削力的仿真图形。一直调不出来,望各位大神赐教!!!
clear
ap=3;ae=8;
c=0.1;
n=10000;
D=16;N=2;beta=pi*30/180;
fi_st=0;fi_ex=pi/3;
Ktc=675.5;Krc=277.9;Kac=195.8;
Kte=6.6;Kre=0.43;Kae=7.5;
deta_fi=pi/100;
deta_ap=0.01;
fi_p=2*pi/N;
K=2*pi/deta_fi;
L=ap/deta_ap;
for i=1:K
Fx(i)=0;Fy(i)=0;Fz(i)=0;Ft(i)=0;fi(i)=0;
fi(i)=fi_st+i*deta_fi;
for k=1:N;
fi_1=fi(i)+(k-1)*fi_p;
fi_2=fi_1;
for j=1:L;
a(j)=j*deta_ap;
fi_2=fi_1-2*tan(beta)*a(j)/D;
if fi_st<=fi_2<=fi_st;
h(j)=c*sin(fi_2);
deta_Ft=deta_ap*(Ktc*h(j)+Kte);
deta_Fr=deta_ap*(Kre*h(j)+Kre);
deta_Fa=deta_ap*(Kac*h(j)+Kae);
deta_Fx=-deta_Ft*cos(fi_2)-deta_Fr*sin(fi_2);
deta_Fy=deta_Ft*sin(fi_2)-deta_Fr*cos(fi_2);
deta_Fz=deta_Fa;
Fx(i)=Fx(i)+deta_Fx;
Fy(i)=Fy(i)+deta_Fy;
Fz(i)=Fz(i)+deta_Fz;
else ap>ap_st;
ap=ap+deta_ap;
end
F(i)=(Fx(i)^2+Fy(i)^2)^(1/2);
Tc(i)=Ft(i)*D/2;
end
end
end
plot(i,Fx);
plot(i.Fy);
plot(i.Fz);
plot(i,Tc)
danxiaojuan 当前离线   回复时引用此帖
回复

书签

主题工具
显示模式

发帖规则
不可以发表新主题
可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码
Trackbacks are 禁用
Pingbacks are 禁用
Refbacks are 启用



所有时间均为北京时间。现在的时间是 11:43 AM


Powered by vBulletin
版权所有 ©2000 - 2017,Jelsoft Enterprises Ltd.
陕ICP备07001583号

SEO by vBSEO ©2009, Crawlability, Inc.