求直廓环面蜗杆实际车削步骤

东海fyh126

lingcheng 发表于 2014-5-29 18:34
: [" T& }" A2 L& l' F我大四毕业论文，实在图书馆没东西。捉急啊

, ]  A: M+ w( E8 Q9 Q. s, P简单点说就是角相位和进刀量的函数关系问题，，，，，，，，

lingcheng

东海fyh126 发表于 2014-5-29 18:38
( x5 z9 c0 i* j" u简单点说就是角相位和进刀量的函数关系问题，，，，，，，，

1 H; g5 C. O8 B% S  ^0 ^8 H你也是大神啊，我有刀具轨迹的MATLAB代码，但是看不懂，尴尬

lingcheng

lingcheng 发表于 2014-5-29 14:53
% s1 _, Y* P/ l+ T大神求带啊，关于实例车削直廓环面蜗杆，包括齿根和齿顶

大神带我飞吧

lingcheng

function worm3
clear
clc
8 R) y5 r$ K4 c6 E, R1 J' J, la=240.65;%中心距
3 v/ E* b: {- r" \0 T8 S# l8 z: uz1=1;%蜗杆头数
z2=40;%齿轮齿数
! z3 q' G$ b$ _" Ni12=z2/z1;%传动比
zc=4;%蜗杆包围蜗轮齿数
d1=80;%蜗杆分度圆直径
d2=2*a-d1;%蜗轮分度圆直径
1 T( x: U3 }8 x* Mmt=d2/z2;%蜗轮端面模数
9 l  }9 _, d3 N0 y7 z1 S7 yda2=417.36;%蜗轮齿顶圆直径
df2=381.24;%蜗轮齿根圆直径
1 r7 F" y/ g8 r$ \4 i8 G# Q& V* t+ n: S0 Q: hha=(da2-d2)/2;%齿顶高
" d8 _# P* v, rhf=(d2-df2)/2;%齿根高
6 M8 ^: ^; @+ M- Y/ C  k9 tda1=d1+2*ha;%蜗杆顶圆直径
df1=d1-2*hf;%蜗杆根圆直径
/ z; K! y+ h/ V5 B3 R3 C( b, BRf1=a-df1/2;%蜗杆齿根圆弧半径
) D5 \* l! S) y; g/ @' GRa1=a-da1/2;%蜗杆齿顶圆弧半径
7 X  R/ D  @; O3 i4 R1 {gama=atan(d2/d1/i12);%蜗杆喉部分度圆导程角
1 k, i9 I9 A1 {$ P% p* N7 i1 rtau=2*pi/z2;%齿距角
$ }  E# g0 P0 [  J. S! Udb=153;%主基圆直径
alpha=asin(db/d2);%蜗轮分度圆压力角
phi_w=0.5*tau*(zc-0.45);%蜗杆包围蜗轮的工作半角
phi_0=alpha-phi_w;%工作起始角
Lw=d2*sin(phi_w);%蜗杆工作部分长度
. P: ]) P: D& [8 R! Vbeta=11/180*pi;%母平面倾斜角
, [" M, \- N* {; s' W- Mphi0min=phi_0;
phi0max=phi_0+2*phi_w;
u0=((d2/2)^2-(db/2)^2)^0.5;
phi0=linspace(phi0min,phi0max,300);
) _& N/ j* z) f( z6 DR=linspace(Ra1,Rf1,10);
phi1=phi0.*i12;
3 l; n) ]5 g  Q7 Fcath=2*alpha+162/z2/180*pi;
7 O5 Z& G$ U% |5 \  h
4 S+ T8 n1 [" B1 h2 Q, n  b& e* W%预设矩阵
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));%蜗杆乙面
* z! D  e: p, mx14=zeros(length(R),length(phi0));y14=zeros(length(R),length(phi0));z14=zeros(length(R),length(phi0));%甲面的螺旋线形式
; W) Z9 ^/ n! e) `; \+ p2 U" R' \8 Vtheta2=zeros(length(R),length(phi0));R2=zeros(length(R),length(phi0));Z2=zeros(length(R),length(phi0));%甲面的加工参数

# R* g  e' v2 K* nfigure(1)
for i=1:length(phi0)
& n' ^0 U' X1 v9 F; ^) c7 Y6 g  E    for j=1:length(R)
  i- R4 D, m* s+ S" @9 }2 _' e% ^        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)];
* G, B7 W$ ]7 k  I' U  C  m6 K' p        options=optimset('Display','off','TolFun',1e-7);
        x0=[u0 0];
% n; j' s  D: ~9 T        x=fsolve(f1,x0,options);
' N& T; z: Q6 t        u1(i,j)=x(1);
        v1(i,j)=x(2);
        x01(i,j)=u1(i,j);
9 u& G2 }- Z: a/ V) U% ]        y01(i,j)=v1(i,j).*sin(beta)-db/2;
8 b) ~2 n* L# Y" p  i        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));
( U3 |( i0 H4 B        z11(i,j)=-sin(phi0(i)).*x01(i,j)-cos(phi0(i)).*y01(i,j);
    end
    plot3(x11(i,,y11(i,,z11(i,,'r')
2 n7 \7 @# _# J0 m    hold on
' z0 e* e5 ~+ N+ {+ K5 K. `end
axis equal

2 A1 q7 j3 h! N  xfigure(1)
, y* E2 D7 I( ~! Ofor i=1:length(phi0)
    for j=1:length(R)
4 D3 M' i3 v3 |4 Z% u& 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)];
" L9 L! d0 N% l; p. P        options=optimset('Display','off','TolFun',1e-7);
2 j' h- l! X$ H, u- i; }) o- S7 z        x0=[u0 0];
        x=fsolve(f2,x0,options);
        u2(i,j)=x(1);
        v2(i,j)=x(2);
% b6 A3 v# y4 n2 T& q        x03(i,j)=u2(i,j);
$ y# U3 P  B' J% N6 b* b        y03(i,j)=v2(i,j).*sin(beta)+db/2;
        z03(i,j)=v2(i,j).*cos(beta);
  }' C4 H/ I' [0 ^        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);
) k. s0 W. H7 i+ c; q+ c: W        z02(i,j)=z03(i,j);
7 t  Q* U1 X2 g/ r' w7 _; q# }        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));
  a) |) Y: ]1 v# \+ w5 {$ v2 k        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);
6 q  u  [  Y; W: o% E: ~    end
    plot3(x12(i,:),y12(i,:),z12(i,:))
    hold on
: c) A8 @/ i9 Fend
axis equal

" b9 q& p% N8 i0 Cfigure(2)
# P) Y7 N2 _0 afor j=1:length(R)
    plot3(x11(:,j),y11(:,j),z11(:,j),'r')
' u# |9 C5 Q* a$ P" L  A6 A6 y- A8 A    hold on
end
axis equal
figure(2)
for j=1:length(R)
    plot3(x12(:,j),y12(:,j),z12(:,j))
    hold on
end
axis equal
- @4 F) J) }) f' g& P
. p1 y+ J& d  d" [/ Ax13=x11';%转换为螺旋线的格式
y13=y11';
+ f; ^+ V, T8 g2 wz13=z11';
x14=x12';
' ^4 L+ t$ U! _4 ey14=y12';
1 g) U. S, W. F1 I7 D) ~6 q3 H" Z8 f, xz14=z12';
, T) H0 l" Y& [
9 g: I1 I  v' {; C" }( A. Kfor i=1:length(R)
    for j=1:length(phi0)
9 J. n! l  |6 o" \* k) ^1 x. c        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);
+ N0 W+ [. ~' Z* ]% ^        theta2(i,j)=cal_theta(x14(i,j),y14(i,j))/pi*180;
; U3 w2 \- X+ f9 _: w) C  W        R2(i,j)=(x14(i,j)^2+y14(i,j)^2)^0.5;
        Z2(i,j)=z14(i,j);
9 j; {# _% R$ \% M- z    end
end

%角度转换为绝对式
/ C4 s, m; m2 [% n, dfor  i=1:length(R)
: O  I5 a" B9 O, [, v- W6 L: f    for j=1length(phi0)-1)
# |) K) C6 E- ?5 b+ m0 z. T4 c8 W6 U        if theta1(i,j+1)<theta1(i,j)
            for k=(j+1):length(phi0)
                theta1(i,k)=theta1(i,k)+360;
$ l, Z3 P5 ^6 {6 J            end
        end
# w* X0 i" F6 m/ _' U3 D1 t         if theta2(i,j+1)<theta2(i,j)
            for k=(j+1):length(phi0)
& D: |/ N% E/ ~: Z  H) i* y                theta2(i,k)=theta2(i,k)+360;
            end
        end
    end
end

assignin('base','X1',R1)
% [/ o& Q- J2 ]( K/ Eassignin('base','Z1',Z1)
assignin('base','C1',theta1)
+ }. J3 k3 \0 Z3 c; J' rassignin('base','X2',R2)
  C6 _. j) m; ]assignin('base','Z2',Z2)
assignin('base','C2',theta2)
& }% A7 b2 Y2 y3 J9 B
%G代码的输出
% fid=fopen('turncode.txt','wt');
% fprintf(fid,'O0001\n');
% for  i=1:length(R)
%     for j=1:length(phi0)
& Y4 E1 k5 e) I( H; `%         fprintf(fid,'X=%0.4f\t',R1(i,j))
%         fprintf(fid,'Z=%0.4f\t',Z1(i,j))
/ c% s. c, B; ]$ q0 H' c%         fprintf(fid,'C=%0.4f\n',theta1(i,j))
%     end
% end
2 {5 R* [' E0 L0 B: J) e! T
4 j, o  K2 \! O
  u6 L& \* x& C/ s4 v( ]9 J

function theta=cal_theta(x1,y1)
if (x1>0)&&(y1<0)
( |6 r# O0 r0 O    theta=acos(x1./(x1^2+y1^2)^0.5);%电机转是角度还是弧度？
4 m1 @% Z9 C, g# U, Jelseif (x1<0)&&(y1<0);
1 d  [; Z- \" C0 U9 ]% _theta=acos(x1./(x1^2+y1^2)^0.5);
( j3 G; J0 l: Z- u' `elseif (x1<0)&&(y1>0);
    theta=2*pi-acos(x1./(x1^2+y1^2)^0.5);
elseif (x1>0)&&(y1>0);
. x' ?$ l1 q9 D" J    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);
- o$ w& Y0 _% x9 W! ?! w& x" Relseif (x1>0)&&(y1==0);
    theta=0;
+ @5 p2 O$ }' ~8 u( ^, `elseif (x1<0)&&(y1==0);
    theta=pi;
elseif (x1==0)&&(y1<0);
    theta=pi/2;
4 C) ]4 S) o- ^3 \/ Telseif (x1==0)&&(y1>0);
: _4 P3 @& d) l$ I3 |! N    theta=3*pi/2;
end

东海fyh126

我也看不懂，，，，，

jxcrxssh

看了头晕，有这么复杂，算了路过！
/ o$ t0 K& G+ q( l' A

car1982

不懂


$ U; _3 y9 g7 J" s7 Y% _