求直廓環面蝸桿實際車削步驟

東海fyh126 · 發表于 2014-5-29 18:38:36

lingcheng 發表于 2014-5-29 18:34
0 ~9 M& [/ G* V5 Q' ?, f我大四畢業論文，實在圖書館沒東西。捉急啊

簡單點說就是角相位和進刀量的函數關系問題，，，，，，，，

lingcheng · 發表于 2014-5-29 18:40:53

東海fyh126 發表于 2014-5-29 18:38
4 I* H) k5 z' n3 j. j6 p簡單點說就是角相位和進刀量的函數關系問題，，，，，，，，

你也是大神啊，我有刀具軌跡的MATLAB代碼，但是看不懂，尷尬

lingcheng · 發表于 2014-5-29 18:57:36

lingcheng 發表于 2014-5-29 14:53 8 ]* u; ]$ e* ]" K
大神求帶啊，關于實例車削直廓環面蝸桿，包括齒根和齒頂

大神帶我飛吧

lingcheng · 發表于 2014-5-29 18:58:45

function worm3
clear
clc
a=240.65;%中心距
z1=1;%蝸桿頭數
z2=40;%齒輪齒數
i12=z2/z1;%傳動比
zc=4;%蝸桿包圍蝸輪齒數
d1=80;%蝸桿分度圓直徑
d2=2*a-d1;%蝸輪分度圓直徑
mt=d2/z2;%蝸輪端面模數
da2=417.36;%蝸輪齒頂圓直徑
df2=381.24;%蝸輪齒根圓直徑
ha=(da2-d2)/2;%齒頂高
hf=(d2-df2)/2;%齒根高
da1=d1+2*ha;%蝸桿頂圓直徑
df1=d1-2*hf;%蝸桿根圓直徑
Rf1=a-df1/2;%蝸桿齒根圓弧半徑
Ra1=a-da1/2;%蝸桿齒頂圓弧半徑
gama=atan(d2/d1/i12);%蝸桿喉部分度圓導程角
tau=2*pi/z2;%齒距角
db=153;%主基圓直徑
alpha=asin(db/d2);%蝸輪分度圓壓力角
phi_w=0.5*tau*(zc-0.45);%蝸桿包圍蝸輪的工作半角
phi_0=alpha-phi_w;%工作起始角
Lw=d2*sin(phi_w);%蝸桿工作部分長度
beta=11/180*pi;%母平面傾斜角
phi0min=phi_0;
phi0max=phi_0+2*phi_w;
u0=((d2/2)^2-(db/2)^2)^0.5;
phi0=linspace(phi0min,phi0max,300);
R=linspace(Ra1,Rf1,10);
phi1=phi0.*i12;
cath=2*alpha+162/z2/180*pi;

%預設矩陣
u1=zeros(length(phi0),length(R));v1=zeros(length(phi0),length(R));%甲面參數
x01=zeros(length(phi0),length(R));y01=zeros(length(phi0),length(R));z01=zeros(length(phi0),length(R));%甲面的平面
x11=zeros(length(phi0),length(R));y11=zeros(length(phi0),length(R));z11=zeros(length(phi0),length(R));%蝸桿甲面
x13=zeros(length(R),length(phi0));y13=zeros(length(R),length(phi0));z13=zeros(length(R),length(phi0));%甲面的螺旋線形式
theta1=zeros(length(R),length(phi0));R1=zeros(length(R),length(phi0));Z1=zeros(length(R),length(phi0));%甲面的加工參數

u2=zeros(length(phi0),length(R));v2=zeros(length(phi0),length(R));%乙面參數
x03=zeros(length(phi0),length(R));y03=zeros(length(phi0),length(R));z03=zeros(length(phi0),length(R));%乙面的平面
x02=zeros(length(phi0),length(R));y02=zeros(length(phi0),length(R));z02=zeros(length(phi0),length(R));%乙面的平面
x12=zeros(length(phi0),length(R));y12=zeros(length(phi0),length(R));z12=zeros(length(phi0),length(R));%蝸桿乙面
x14=zeros(length(R),length(phi0));y14=zeros(length(R),length(phi0));z14=zeros(length(R),length(phi0));%甲面的螺旋線形式
theta2=zeros(length(R),length(phi0));R2=zeros(length(R),length(phi0));Z2=zeros(length(R),length(phi0));%甲面的加工參數

figure(1)
for i=1:length(phi0)
for j=1:length(R)
      f1=@(x)[(a-((x(1)*cos(phi0(i))-(x(2)*sin(beta)-db/2)*sin(phi0(i))-a)^2+(x(2)*cos(beta))^2)^0.5)^2+(x(1)*sin(phi0(i))+(x(2)*sin(beta)-db/2)*cos(phi0(i)))^2-R(j)^2;(x(1)*(cos(beta)/i12+sin(beta)*cos(phi0(i)))+(db/2*sin(phi0(i))-a)*sin(beta))/sin(phi0(i))-x(2)];
      options=optimset('Display','off','TolFun',1e-7);
      x0=[u0 0];
      x=fsolve(f1,x0,options);
      u1(i,j)=x(1);
      v1(i,j)=x(2);
      x01(i,j)=u1(i,j);
      y01(i,j)=v1(i,j).*sin(beta)-db/2;
      z01(i,j)=v1(i,j).*cos(beta);
      x11(i,j)=-cos(phi1(i)).*cos(phi0(i)).*x01(i,j)+cos(phi1(i)).*sin(phi0(i)).*y01(i,j)-sin(phi1(i)).*z01(i,j)+a.*cos(phi1(i));
      y11(i,j)=sin(phi1(i)).*cos(phi0(i)).*x01(i,j)-sin(phi1(i)).*sin(phi0(i)).*y01(i,j)-cos(phi1(i)).*z01(i,j)-a.*sin(phi1(i));
      z11(i,j)=-sin(phi0(i)).*x01(i,j)-cos(phi0(i)).*y01(i,j);
end
plot3(x11(i,

,y11(i,

,z11(i,

,'r')
hold on
end
axis equal

figure(1)
for i=1:length(phi0)
for j=1:length(R)
      f2=@(x)[(a-(((x(1)*cos(cath)+(x(2)*sin(beta)+db/2)*sin(cath))*cos(phi0(i))-(-x(1)*sin(cath)+(x(2)*sin(beta)+db/2)*cos(cath))*sin(phi0(i))-a)^2+(x(2)*cos(beta))^2)^0.5)^2+((x(1)*cos(cath)+(x(2)*sin(beta)+db/2)*sin(cath))*sin(phi0(i))+(-x(1)*sin(cath)+(x(2)*sin(beta)+db/2)*cos(cath))*cos(phi0(i)))^2-R(j)^2;(x(1)*(cos(beta)/i12+sin(beta)*cos(phi0(i)-cath))+(-db/2*sin(phi0(i)-cath)-a)*sin(beta))/sin(phi0(i)-cath)-x(2)];
      options=optimset('Display','off','TolFun',1e-7);
      x0=[u0 0];
      x=fsolve(f2,x0,options);
      u2(i,j)=x(1);
      v2(i,j)=x(2);
      x03(i,j)=u2(i,j);
      y03(i,j)=v2(i,j).*sin(beta)+db/2;
      z03(i,j)=v2(i,j).*cos(beta);
      x02(i,j)=x03(i,j).*cos(cath)+y03(i,j).*sin(cath);
      y02(i,j)=-x03(i,j).*sin(cath)+y03(i,j).*cos(cath);
      z02(i,j)=z03(i,j);
      x12(i,j)=-cos(phi1(i)).*cos(phi0(i)).*x02(i,j)+cos(phi1(i)).*sin(phi0(i)).*y02(i,j)-sin(phi1(i)).*z02(i,j)+a.*cos(phi1(i));
      y12(i,j)=sin(phi1(i)).*cos(phi0(i)).*x02(i,j)-sin(phi1(i)).*sin(phi0(i)).*y02(i,j)-cos(phi1(i)).*z02(i,j)-a.*sin(phi1(i));
      z12(i,j)=-sin(phi0(i)).*x02(i,j)-cos(phi0(i)).*y02(i,j);
end
plot3(x12(i,:),y12(i,:),z12(i,:))
hold on
end
axis equal

figure(2)
for j=1:length(R)
plot3(x11(:,j),y11(:,j),z11(:,j),'r')
hold on
end
axis equal
figure(2)
for j=1:length(R)
plot3(x12(:,j),y12(:,j),z12(:,j))
hold on
end
axis equal

x13=x11';%轉換為螺旋線的格式
y13=y11';
z13=z11';
x14=x12';
y14=y12';
z14=z12';

for i=1:length(R)
for j=1:length(phi0)
      theta1(i,j)=cal_theta(x13(i,j),y13(i,j))/pi*180;
      R1(i,j)=(x13(i,j)^2+y13(i,j)^2)^0.5;%車床X方向為直徑值
      Z1(i,j)=z13(i,j);
      theta2(i,j)=cal_theta(x14(i,j),y14(i,j))/pi*180;
      R2(i,j)=(x14(i,j)^2+y14(i,j)^2)^0.5;
      Z2(i,j)=z14(i,j);
end
end

%角度轉換為絕對式
for  i=1:length(R)
for j=1

length(phi0)-1)
      if theta1(i,j+1)<theta1(i,j)
         for k=(j+1):length(phi0)
            theta1(i,k)=theta1(i,k)+360;
         end
      end
      if theta2(i,j+1)<theta2(i,j)
         for k=(j+1):length(phi0)
            theta2(i,k)=theta2(i,k)+360;
         end
      end
end
end

assignin('base','X1',R1)
assignin('base','Z1',Z1)
assignin('base','C1',theta1)
assignin('base','X2',R2)
assignin('base','Z2',Z2)
assignin('base','C2',theta2)

%G代碼的輸出
% fid=fopen('turncode.txt','wt');
% fprintf(fid,'O0001\n');
% for  i=1:length(R)
%    for j=1:length(phi0)
%       fprintf(fid,'X=%0.4f\t',R1(i,j))
%       fprintf(fid,'Z=%0.4f\t',Z1(i,j))
%       fprintf(fid,'C=%0.4f\n',theta1(i,j))
%    end
% end

function theta=cal_theta(x1,y1)
if (x1>0)&&(y1<0)
theta=acos(x1./(x1^2+y1^2)^0.5);%電機轉是角度還是弧度？
elseif (x1<0)&&(y1<0);
theta=acos(x1./(x1^2+y1^2)^0.5);
elseif (x1<0)&&(y1>0);
theta=2*pi-acos(x1./(x1^2+y1^2)^0.5);
elseif (x1>0)&&(y1>0);
theta=2*pi-acos(x1./(x1^2+y1^2)^0.5);
elseif (x1>0)&&(y1>0);
theta=2*pi-acos(x1./(x1^2+y1^2)^0.5);
elseif (x1>0)&&(y1==0);
theta=0;
elseif (x1<0)&&(y1==0);
theta=pi;
elseif (x1==0)&&(y1<0);
theta=pi/2;
elseif (x1==0)&&(y1>0);
theta=3*pi/2;
end

東海fyh126 · 發表于 2014-5-29 22:26:57

我也看不懂，，，，，

jxcrxssh · 發表于 2014-6-27 22:00:58

看了頭暈，有這么復雜，算了路過！

car1982 · 發表于 2014-6-28 09:32:24

不懂

		自動登錄	找回密碼
密碼			注冊會員

国产日韩欧美久久久精品图片|国产综合有码无码中文字幕|国产一区二区综合视频|国产亚洲精品电影网站在线观看|国产精品一区在线

求直廓環面蝸桿實際車削步驟