求直廓環(huán)面蝸桿實(shí)際車削步驟

東海fyh126

lingcheng 發(fā)表于 2014-5-29 18:34
# I- Y) m1 o3 f+ ?2 y我大四畢業(yè)論文，實(shí)在圖書館沒東西。捉急啊

簡單點(diǎn)說就是角相位和進(jìn)刀量的函數(shù)關(guān)系問題，，，，，，，，

lingcheng

東海fyh126 發(fā)表于 2014-5-29 18:38
簡單點(diǎn)說就是角相位和進(jìn)刀量的函數(shù)關(guān)系問題，，，，，，，，

7 Y/ W( {8 P2 m4 o: w, P' {你也是大神啊，我有刀具軌跡的MATLAB代碼，但是看不懂，尷尬
9 a4 p; |" u: {+ q+ J8 D0 y

lingcheng

lingcheng 發(fā)表于 2014-5-29 14:53
大神求帶啊，關(guān)于實(shí)例車削直廓環(huán)面蝸桿，包括齒根和齒頂

6 e6 }3 D+ Y$ i- |0 |% F) H大神帶我飛吧
: h, M( g2 e& Z; l

lingcheng

function worm3
0 @+ S* v( E, w% S+ eclear
/ U1 e. n% l  s% kclc
1 Q( i4 c  A$ y6 F9 ja=240.65;%中心距
. R' p6 q- d/ }  F  {% o/ m  Xz1=1;%蝸桿頭數(shù)
z2=40;%齒輪齒數(shù)
i12=z2/z1;%傳動比
zc=4;%蝸桿包圍蝸輪齒數(shù)
d1=80;%蝸桿分度圓直徑
d2=2*a-d1;%蝸輪分度圓直徑
mt=d2/z2;%蝸輪端面模數(shù)
da2=417.36;%蝸輪齒頂圓直徑
- ?$ W6 L; u  b! D( j5 z% ldf2=381.24;%蝸輪齒根圓直徑
. _% S" L' ?; C. H+ mha=(da2-d2)/2;%齒頂高
hf=(d2-df2)/2;%齒根高
  H, T) y5 s7 F% f8 D& I; yda1=d1+2*ha;%蝸桿頂圓直徑
df1=d1-2*hf;%蝸桿根圓直徑
: @3 F+ H6 `& o1 T- iRf1=a-df1/2;%蝸桿齒根圓弧半徑
Ra1=a-da1/2;%蝸桿齒頂圓弧半徑
; a: m7 L8 B- e, ngama=atan(d2/d1/i12);%蝸桿喉部分度圓導(dǎo)程角
tau=2*pi/z2;%齒距角
3 v& K0 C9 a" H: g9 x; f- Zdb=153;%主基圓直徑
" i7 ]& ]) ?" a9 y5 qalpha=asin(db/d2);%蝸輪分度圓壓力角
phi_w=0.5*tau*(zc-0.45);%蝸桿包圍蝸輪的工作半角
phi_0=alpha-phi_w;%工作起始角
1 A7 ~( \4 ?8 h% ]0 z! t. q' \1 A' }Lw=d2*sin(phi_w);%蝸桿工作部分長度
( o' _: [' B. vbeta=11/180*pi;%母平面傾斜角
9 w5 Q! t8 x" a: j7 y: X* Y& Gphi0min=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;
* n; F2 ]8 k5 e7 ]9 [: Q2 ]- E9 P" pcath=2*alpha+162/z2/180*pi;

. ?7 K! H: I+ a; B%預(yù)設(shè)矩陣
u1=zeros(length(phi0),length(R));v1=zeros(length(phi0),length(R));%甲面參數(shù)
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));%蝸桿甲面
5 V3 p6 T( M6 H1 Yx13=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));%甲面的加工參數(shù)
  k! A% W, @$ ^0 j3 x! `# Q
/ U! V' B7 J9 ^u2=zeros(length(phi0),length(R));v2=zeros(length(phi0),length(R));%乙面參數(shù)
x03=zeros(length(phi0),length(R));y03=zeros(length(phi0),length(R));z03=zeros(length(phi0),length(R));%乙面的平面
& t) X; ~' b4 G% G  Jx02=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));%蝸桿乙面
0 N9 V/ d+ g( r* |* r: Y9 Jx14=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));%甲面的加工參數(shù)

$ ]$ U3 q* w, m+ @, d/ vfigure(1)
for i=1:length(phi0)
/ l) x; d/ r% r- g4 @% Y) o    for j=1:length(R)
( K7 c* B' y* [) y% Z, i. 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];
( I% u1 j- f) L! O2 w        x=fsolve(f1,x0,options);
2 M: h2 A+ t2 z' M" O4 B        u1(i,j)=x(1);
        v1(i,j)=x(2);
        x01(i,j)=u1(i,j);
- `7 ?. D* a- C+ u6 x( _        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));
  n' n2 M5 ]. h8 s( Q! A0 ?( [( h        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);
' l# t& q" |, ]" I; U5 h    end
    plot3(x11(i,,y11(i,,z11(i,,'r')
    hold on
: s/ o: h6 ^& J, X- T( Kend
axis equal

figure(1)
( W' R; s4 k0 ~) f' m6 [4 d$ a2 \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);
. f: v/ x% q5 M+ h        x0=[u0 0];
: d" D' M, V$ @! K$ y3 {        x=fsolve(f2,x0,options);
" j; d/ c9 Y0 O3 n        u2(i,j)=x(1);
) _8 o* C' C7 _& V. M1 |* n        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);
& ?1 Q9 w! i( P# W0 U6 F8 \5 f        x02(i,j)=x03(i,j).*cos(cath)+y03(i,j).*sin(cath);
' k5 P/ |2 p/ V' J; r- L        y02(i,j)=-x03(i,j).*sin(cath)+y03(i,j).*cos(cath);
; P5 p" k2 i' e# |" o$ x0 Y        z02(i,j)=z03(i,j);
$ p$ ^3 i" R' X/ v        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));
8 v7 S) `+ [5 V7 ?3 g- ^        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
0 ?4 q& ~3 C+ [8 }5 Q( ^; S    plot3(x12(i,:),y12(i,:),z12(i,:))
/ r: P. Y6 r+ h    hold on
5 _2 d/ H% r+ @% E6 qend
5 ~8 v/ Q- o4 R+ V9 @1 C8 L: zaxis equal

figure(2)
for j=1:length(R)
/ z0 f5 G$ \/ q! g) X& ]7 }    plot3(x11(:,j),y11(:,j),z11(:,j),'r')
& u& g) t- @/ O8 Z3 i    hold on
end
, u" q- l" V: t& u0 ?3 caxis equal
figure(2)
6 c$ a( a; H/ @1 a( `( }# f+ ufor j=1:length(R)
    plot3(x12(:,j),y12(:,j),z12(:,j))
$ N* u, C2 V- _    hold on
3 o3 A- v- z2 |& J! Q" Iend
axis equal

4 U- t' ~& Z& D. f, M3 n5 D1 |5 ^& Ax13=x11';%轉(zhuǎn)換為螺旋線的格式
y13=y11';
z13=z11';
8 Q! ^4 E3 ]) |9 Z" nx14=x12';
y14=y12';
' a, S7 p) I  O7 [0 n; _z14=z12';

' h9 Z) A& l& tfor 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方向?yàn)橹睆街?br /> . M  N+ x6 Z7 @  {9 [) Y        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);
$ B& z8 J0 `0 s2 |1 u    end
end
: t2 T6 N( i% E( e) H
4 j, D0 g: f6 U( Q' B* u9 m%角度轉(zhuǎn)換為絕對式
for  i=1:length(R)
0 [9 ]. D, Z/ H+ {    for j=1length(phi0)-1)
        if theta1(i,j+1)<theta1(i,j)
8 G( Q" q5 d. |$ Z( h$ H, }3 M            for k=(j+1):length(phi0)
                theta1(i,k)=theta1(i,k)+360;
            end
% L1 R2 @' l7 U# J0 f        end
, D: C' X& f, K, }         if theta2(i,j+1)<theta2(i,j)
: b7 X. [& D+ L8 R' j            for k=(j+1):length(phi0)
                theta2(i,k)=theta2(i,k)+360;
6 g) n8 _5 ?: j) Z  b" m6 X            end
        end
    end
end

! y' ^+ M; X& G4 V2 zassignin('base','X1',R1)
' E# {9 k! v7 u# Y* gassignin('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');
4 M: Z3 R* s9 ~4 s% fprintf(fid,'O0001\n');
% for  i=1:length(R)
%     for j=1:length(phi0)
- j0 a- {3 I3 `%         fprintf(fid,'X=%0.4f\t',R1(i,j))
%         fprintf(fid,'Z=%0.4f\t',Z1(i,j))
. `3 Z) T+ c5 ?% K# A( u7 N& }%         fprintf(fid,'C=%0.4f\n',theta1(i,j))
%     end
) ?  Y" ?9 Q& U. f% end


' ~) o5 _  k% E, [, F* \

function theta=cal_theta(x1,y1)
if (x1>0)&&(y1<0)
    theta=acos(x1./(x1^2+y1^2)^0.5);%電機(jī)轉(zhuǎn)是角度還是弧度？
; U+ ~& I2 ]+ r4 \- gelseif (x1<0)&&(y1<0);
theta=acos(x1./(x1^2+y1^2)^0.5);
. r+ E; Z  o/ J* p, E0 F- x1 nelseif (x1<0)&&(y1>0);
! X& s6 A1 A6 ~- {0 @) }3 y    theta=2*pi-acos(x1./(x1^2+y1^2)^0.5);
elseif (x1>0)&&(y1>0);
  A+ ?! j! X, G. ^; i! P    theta=2*pi-acos(x1./(x1^2+y1^2)^0.5);
9 n$ l# |" D8 T+ c& X  R4 u( Celseif (x1>0)&&(y1>0);
    theta=2*pi-acos(x1./(x1^2+y1^2)^0.5);
  V, v) D. G  f. f4 s2 {elseif (x1>0)&&(y1==0);
    theta=0;
elseif (x1<0)&&(y1==0);
    theta=pi;
( R+ ~' \0 S. |" Celseif (x1==0)&&(y1<0);
& m2 o3 K& |8 @. E( u    theta=pi/2;
3 w) P5 Y! x9 o' z3 w: n3 o9 d1 w+ relseif (x1==0)&&(y1>0);
    theta=3*pi/2;
% m; n6 j; @, {" F, Z" Dend

東海fyh126

我也看不懂，，，，，

jxcrxssh

看了頭暈，有這么復(fù)雜，算了路過！

car1982

不懂
- T8 D8 Z! y  C1 {

/ J# e' {' g( U* c5 z1 X9 e) j+ n