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

東海fyh126

lingcheng 發表于 2014-5-29 18:34
我大四畢業論文，實在圖書館沒東西。捉急啊

6 A2 d( C! w6 F0 v& p7 o簡單點說就是角相位和進刀量的函數關系問題，，，，，，，，

lingcheng

東海fyh126 發表于 2014-5-29 18:38
5 S/ E$ C! {1 }& Q簡單點說就是角相位和進刀量的函數關系問題，，，，，，，，

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

lingcheng

lingcheng 發表于 2014-5-29 14:53
大神求帶啊，關于實例車削直廓環面蝸桿，包括齒根和齒頂

0 B8 k0 c2 x) H+ X+ ]" E6 y大神帶我飛吧

lingcheng

function worm3
clear
# N, q! c" `9 _% h9 I6 Dclc
a=240.65;%中心距
1 K' ]3 R7 T- }8 Q) F# p, Rz1=1;%蝸桿頭數
. p" f. h+ W: Q) k) r3 oz2=40;%齒輪齒數
i12=z2/z1;%傳動比
zc=4;%蝸桿包圍蝸輪齒數
d1=80;%蝸桿分度圓直徑
d2=2*a-d1;%蝸輪分度圓直徑
' @6 n1 Y1 s; pmt=d2/z2;%蝸輪端面模數
da2=417.36;%蝸輪齒頂圓直徑
df2=381.24;%蝸輪齒根圓直徑
5 D0 q; Z0 }6 V! [ha=(da2-d2)/2;%齒頂高
hf=(d2-df2)/2;%齒根高
3 R) x2 V$ x( ^5 Ada1=d1+2*ha;%蝸桿頂圓直徑
1 T# {3 D; m! i+ Pdf1=d1-2*hf;%蝸桿根圓直徑
8 J! b- }8 s* x- A, P+ nRf1=a-df1/2;%蝸桿齒根圓弧半徑
1 c: y9 G. e& u$ VRa1=a-da1/2;%蝸桿齒頂圓弧半徑
+ @9 b- t2 ?2 Cgama=atan(d2/d1/i12);%蝸桿喉部分度圓導程角
tau=2*pi/z2;%齒距角
1 T/ A4 o2 n" o$ C& f  r2 n5 q$ j  rdb=153;%主基圓直徑
alpha=asin(db/d2);%蝸輪分度圓壓力角
phi_w=0.5*tau*(zc-0.45);%蝸桿包圍蝸輪的工作半角
phi_0=alpha-phi_w;%工作起始角
Lw=d2*sin(phi_w);%蝸桿工作部分長度
4 R. g/ s% b0 {; qbeta=11/180*pi;%母平面傾斜角
2 v  @  z" v* M8 d$ M) L7 Fphi0min=phi_0;
$ P7 K+ U) n) b; Ophi0max=phi_0+2*phi_w;
) |* X8 d+ T7 ^) E8 g% Ju0=((d2/2)^2-(db/2)^2)^0.5;
) W4 b5 w$ x. v0 _phi0=linspace(phi0min,phi0max,300);
R=linspace(Ra1,Rf1,10);
! @( p1 v$ }! yphi1=phi0.*i12;
cath=2*alpha+162/z2/180*pi;

%預設矩陣
4 f1 b8 n: k' J! Ku1=zeros(length(phi0),length(R));v1=zeros(length(phi0),length(R));%甲面參數
  @1 `& v7 T4 o' q4 ux01=zeros(length(phi0),length(R));y01=zeros(length(phi0),length(R));z01=zeros(length(phi0),length(R));%甲面的平面
+ A0 E5 S! z/ e% n5 c* [4 Rx11=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));%甲面的加工參數
' O; D2 L& R! K
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));%乙面的平面
! |5 b  z5 x( r5 ax02=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));%蝸桿乙面
6 S" ~( \' |  zx14=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));%甲面的加工參數
0 w: ]% F) I0 _* A' E/ G! r
9 @/ x/ j8 z3 {  ^' _, W7 Efigure(1)
6 {, E  d6 F: b6 m/ L1 Kfor i=1:length(phi0)
/ j: A1 }/ d* R, ^5 M( f0 {% M% G    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)];
0 m4 _  t: r, Q6 {3 B        options=optimset('Display','off','TolFun',1e-7);
        x0=[u0 0];
& N% m) R1 L! J2 j5 w        x=fsolve(f1,x0,options);
        u1(i,j)=x(1);
        v1(i,j)=x(2);
        x01(i,j)=u1(i,j);
/ ?, \& K, Y* Q        y01(i,j)=v1(i,j).*sin(beta)-db/2;
& `, B- F4 p1 o/ t* h8 U% t; {        z01(i,j)=v1(i,j).*cos(beta);
  M$ a+ q4 K* R  v        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));
. H5 f6 Z* D# g, q        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);
7 w  S( O/ H3 u8 Z7 Z- T    end
% |  }; u  F8 L/ I, r    plot3(x11(i,,y11(i,,z11(i,,'r')
    hold on
, y* c. |( u/ {; i$ pend
3 r& s. x7 q1 N; J! L5 l- Uaxis equal

# x# V/ _, {; m7 s9 H+ Xfigure(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);
! r9 m( c9 O4 `1 B2 y3 r/ ?        v2(i,j)=x(2);
        x03(i,j)=u2(i,j);
! s& Z8 l$ ?1 ~6 A- R        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);
2 }) O$ r7 E0 ]+ S$ o8 b0 A2 B- J        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
9 ?; }9 _* a5 A& Y' I5 o" m    plot3(x12(i,:),y12(i,:),z12(i,:))
    hold on
& l& {' N. J4 B4 E- B8 P, p+ send
7 a1 A; U/ k* Z$ kaxis equal

% X2 |& M* D, b9 f3 c. r' w) H' g" Xfigure(2)
for j=1:length(R)
    plot3(x11(:,j),y11(:,j),z11(:,j),'r')
    hold on
end
5 |, k8 b: ~; V5 H8 Gaxis equal
$ F5 m% R% S$ x- ?figure(2)
for j=1:length(R)
6 m3 P0 P. W5 T' I    plot3(x12(:,j),y12(:,j),z12(:,j))
    hold on
  }0 C: y+ ]6 `end
axis equal

1 Z6 s4 J1 H( N4 Y; ex13=x11';%轉換為螺旋線的格式
y13=y11';
5 B3 w9 j  P; G/ hz13=z11';
x14=x12';
6 ^$ {5 e$ C) B( Y# g4 v; j8 gy14=y12';
; \& [# W* R9 [! P5 P& F) k! H% az14=z12';

for i=1:length(R)
; ~+ h' f, e: d1 [, t2 m' u" j    for j=1:length(phi0)
        theta1(i,j)=cal_theta(x13(i,j),y13(i,j))/pi*180;
1 [3 D* A6 G+ V0 U; D! G. N1 M        R1(i,j)=(x13(i,j)^2+y13(i,j)^2)^0.5;%車床X方向為直徑值
* t1 a8 c. z7 V% V. O1 _! Q        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
+ R- Y7 W1 p; P4 Oend

%角度轉換為絕對式
5 o9 s; d8 h% v" d9 ifor  i=1:length(R)
! P; g. F  u8 k9 d0 E7 X: y- W" o    for j=1length(phi0)-1)
9 g0 o: w, _- X, a        if theta1(i,j+1)<theta1(i,j)
            for k=(j+1):length(phi0)
                theta1(i,k)=theta1(i,k)+360;
            end
        end
1 }) k+ M9 c5 ?5 M         if theta2(i,j+1)<theta2(i,j)
# I9 U* G/ C& [3 I3 P            for k=(j+1):length(phi0)
                theta2(i,k)=theta2(i,k)+360;
9 K4 [& u8 t# A: u            end
3 ?( F! W, U8 }0 V9 r  Q        end
; R0 x4 u- @, b( V  ~    end
end
# G. ], s+ Y. D, W
assignin('base','X1',R1)
. I: X! R) P; ?% f/ g" bassignin('base','Z1',Z1)
assignin('base','C1',theta1)
4 s/ J- C4 I7 I. s" j1 iassignin('base','X2',R2)
* H) r. o9 h$ {+ dassignin('base','Z2',Z2)
assignin('base','C2',theta2)

%G代碼的輸出
% fid=fopen('turncode.txt','wt');
3 P( m# M( |5 M  ?7 t% fprintf(fid,'O0001\n');
% for  i=1:length(R)
7 l7 X1 Y9 b1 e4 U8 @%     for j=1:length(phi0)
%         fprintf(fid,'X=%0.4f\t',R1(i,j))
$ i% X* L6 j, V3 f7 z3 [2 O%         fprintf(fid,'Z=%0.4f\t',Z1(i,j))
' E/ f: e) ?# _5 T$ [. Q# F%         fprintf(fid,'C=%0.4f\n',theta1(i,j))
: M! l4 K2 T7 L%     end
% end




5 h; ~7 N7 z, I4 B7 sfunction theta=cal_theta(x1,y1)
! B, I% N) Q1 m; o$ oif (x1>0)&&(y1<0)
    theta=acos(x1./(x1^2+y1^2)^0.5);%電機轉是角度還是弧度？
. h- ?) v. i" @. h! G) c4 Eelseif (x1<0)&&(y1<0);
  [% G/ l# F, ~7 X8 H" x8 z; Ktheta=acos(x1./(x1^2+y1^2)^0.5);
elseif (x1<0)&&(y1>0);
    theta=2*pi-acos(x1./(x1^2+y1^2)^0.5);
/ K- y2 D5 Y& g8 Yelseif (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);
9 q# W! A9 ^" o( C* Q- u8 A" k! ~# T% zelseif (x1>0)&&(y1==0);
" I( }/ `$ Y" t2 ?/ X! a+ s. V    theta=0;
& h: u2 x; ^; yelseif (x1<0)&&(y1==0);
# R9 J4 j- T' H7 d# {    theta=pi;
elseif (x1==0)&&(y1<0);
2 s" h$ _2 `/ W/ f3 A    theta=pi/2;
elseif (x1==0)&&(y1>0);
    theta=3*pi/2;
end

東海fyh126

我也看不懂，，，，，

jxcrxssh

看了頭暈，有這么復雜，算了路過！
: y/ A" U9 M; M6 u! U" D7 K

car1982

不懂