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

東海fyh126

lingcheng 發(fā)表于 2014-5-29 18:34
; R; `  ]% O/ y6 z9 q我大四畢業(yè)論文，實(shí)在圖書館沒東西。捉急啊

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

lingcheng

東海fyh126 發(fā)表于 2014-5-29 18:38
. H5 B9 o: _% K% B: {1 g4 E: t簡(jiǎn)單點(diǎn)說就是角相位和進(jìn)刀量的函數(shù)關(guān)系問題，，，，，，，，

你也是大神啊，我有刀具軌跡的MATLAB代碼，但是看不懂，尷尬
0 Y- y) R9 Y7 H1 B# i- H

lingcheng

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

3 ]/ v+ S: }/ |/ d$ m6 t: O5 _" p大神帶我飛吧
) L' ?4 B! p0 V& S" [3 ?

lingcheng

function worm3
clear
clc
a=240.65;%中心距
; H/ |9 e' Z5 E! o( n; f2 mz1=1;%蝸桿頭數(shù)
z2=40;%齒輪齒數(shù)
- _+ K; j. |9 Q6 ]; d7 A0 \! v1 ii12=z2/z1;%傳動(dòng)比
$ q/ k" ^9 w% D$ y+ czc=4;%蝸桿包圍蝸輪齒數(shù)
2 M% a$ y; L1 v( }2 v1 Md1=80;%蝸桿分度圓直徑
d2=2*a-d1;%蝸輪分度圓直徑
% V, H& l4 H# B7 X3 A. Zmt=d2/z2;%蝸輪端面模數(shù)
da2=417.36;%蝸輪齒頂圓直徑
df2=381.24;%蝸輪齒根圓直徑
ha=(da2-d2)/2;%齒頂高
hf=(d2-df2)/2;%齒根高
0 {9 m4 [! D5 \, A0 ~da1=d1+2*ha;%蝸桿頂圓直徑
- ~4 O5 b6 Q! y% d# p! ]df1=d1-2*hf;%蝸桿根圓直徑
# H( t" F/ h9 t% s( NRf1=a-df1/2;%蝸桿齒根圓弧半徑
1 T4 f0 n  \% KRa1=a-da1/2;%蝸桿齒頂圓弧半徑
gama=atan(d2/d1/i12);%蝸桿喉部分度圓導(dǎo)程角
tau=2*pi/z2;%齒距角
db=153;%主基圓直徑
! Q$ S$ a- t# F4 U8 B# ]. \1 zalpha=asin(db/d2);%蝸輪分度圓壓力角
7 I. x6 l7 M) L. {) r3 J0 |phi_w=0.5*tau*(zc-0.45);%蝸桿包圍蝸輪的工作半角
' _# S9 R8 F6 B* a$ |8 X' Aphi_0=alpha-phi_w;%工作起始角
Lw=d2*sin(phi_w);%蝸桿工作部分長(zhǎng)度
beta=11/180*pi;%母平面傾斜角
phi0min=phi_0;
phi0max=phi_0+2*phi_w;
u0=((d2/2)^2-(db/2)^2)^0.5;
' T8 q% r! @& y  q( A. ophi0=linspace(phi0min,phi0max,300);
R=linspace(Ra1,Rf1,10);
! k( L8 _, u9 y) j2 E6 c2 vphi1=phi0.*i12;
cath=2*alpha+162/z2/180*pi;

; S1 V3 V5 R9 E2 T8 f0 E%預(yù)設(shè)矩陣
" s( B2 H  D' B6 O6 Su1=zeros(length(phi0),length(R));v1=zeros(length(phi0),length(R));%甲面參數(shù)
( f* V  z  ?: Q* p! h; Yx01=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));%甲面的加工參數(shù)

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));%乙面的平面
# ]5 D% o0 E' J9 Z5 J* xx02=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));%蝸桿乙面
7 W6 o8 h0 q: A) x- ax14=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ù)
4 r% v5 q* y1 [  d9 `( h1 C4 ^
figure(1)
for i=1:length(phi0)
    for j=1:length(R)
8 d6 @1 s+ W) M% N, w  B$ B        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)];
# @  L, G' H: m5 a8 k        options=optimset('Display','off','TolFun',1e-7);
        x0=[u0 0];
        x=fsolve(f1,x0,options);
7 B! E5 y0 C( J7 x& {+ R1 p        u1(i,j)=x(1);
        v1(i,j)=x(2);
        x01(i,j)=u1(i,j);
& a! r! }, {) X) l; B        y01(i,j)=v1(i,j).*sin(beta)-db/2;
( Q5 ^9 t7 @" z4 t        z01(i,j)=v1(i,j).*cos(beta);
6 {& i: Z6 I; V- ~2 K        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));
7 u1 w5 u! ]: {. ^+ U1 a  k- M        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));
( l, b) b1 ]: T8 d9 P" g$ `        z11(i,j)=-sin(phi0(i)).*x01(i,j)-cos(phi0(i)).*y01(i,j);
    end
    plot3(x11(i,,y11(i,,z11(i,,'r')
, ?0 Z" B/ ~3 f" Y; z* Z! d    hold on
% m! l  U. |2 h0 |# o2 }0 _7 }9 Nend
axis equal

0 q& A/ f3 z1 A4 dfigure(1)
: T3 O% ?0 M6 i8 _  W+ cfor 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];
) h2 j: g& S6 x6 l! Q+ H2 o        x=fsolve(f2,x0,options);
        u2(i,j)=x(1);
: X& L. @$ }0 c2 _        v2(i,j)=x(2);
& ~& a1 @$ |  }  M        x03(i,j)=u2(i,j);
        y03(i,j)=v2(i,j).*sin(beta)+db/2;
. K/ S; K% q( V  Z- r8 K        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 a1 p% h+ J9 J. O$ f' 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));
3 q" J# \$ Y9 ^; t% 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));
+ t6 k. G. j9 S4 ~        z12(i,j)=-sin(phi0(i)).*x02(i,j)-cos(phi0(i)).*y02(i,j);
2 B+ J' l. K5 P6 D# f) |# H  j4 q) G; s    end
5 A; C" {$ T# E% C6 p! Q5 n8 ^7 o0 d& t    plot3(x12(i,:),y12(i,:),z12(i,:))
    hold on
4 I8 i; C, G, tend
% y2 R: P4 m4 Uaxis equal
. M; D6 l! s% |. t1 Y  s6 K
figure(2)
for j=1:length(R)
1 x8 Z8 c7 p8 n. f4 e    plot3(x11(:,j),y11(:,j),z11(:,j),'r')
. r9 s' v1 _# I0 s. p* [    hold on
end
6 F- ?# s+ _1 j. w. d( ~axis equal
figure(2)
for j=1:length(R)
    plot3(x12(:,j),y12(:,j),z12(:,j))
    hold on
end
axis equal
! N6 I. A, g- V8 `
x13=x11';%轉(zhuǎn)換為螺旋線的格式
2 ]- X5 ]; h# h9 c: ^y13=y11';
z13=z11';
x14=x12';
: x9 T/ p4 X5 M$ {$ Y) ny14=y12';
z14=z12';
; o" s$ e2 U- d/ k# I& t  M
for i=1:length(R)
( _' I, @/ k% j' o2 U* Z" E    for j=1:length(phi0)
        theta1(i,j)=cal_theta(x13(i,j),y13(i,j))/pi*180;
6 @1 D/ w; B: S- b' w2 |        R1(i,j)=(x13(i,j)^2+y13(i,j)^2)^0.5;%車床X方向?yàn)橹睆街?font class="jammer">+ h( ~3 T( }) L0 I3 j# ]
        Z1(i,j)=z13(i,j);
        theta2(i,j)=cal_theta(x14(i,j),y14(i,j))/pi*180;
( g+ k; T7 Z% C: E* @& I) k+ \( O        R2(i,j)=(x14(i,j)^2+y14(i,j)^2)^0.5;
        Z2(i,j)=z14(i,j);
    end
end
' B/ V- j7 Y9 u8 ^( `
8 {$ U1 V+ t. `5 R2 l, m%角度轉(zhuǎn)換為絕對(duì)式
for  i=1:length(R)
" p! K$ H# Q7 e; `' B    for j=1length(phi0)-1)
! y7 s3 w. d* r5 \6 @* Z        if theta1(i,j+1)<theta1(i,j)
            for k=(j+1):length(phi0)
% I* b! ~' K" m: m% m" ]                theta1(i,k)=theta1(i,k)+360;
, V3 a9 }3 X1 f4 t+ k8 ?2 U; V            end
        end
: W* }2 g$ ^! L& P! m) v! e/ g         if theta2(i,j+1)<theta2(i,j)
5 _$ u0 P7 l/ @  v- y            for k=(j+1):length(phi0)
. d; J2 m# x- |2 i+ o                theta2(i,k)=theta2(i,k)+360;
            end
        end
. W# z1 f& G8 p8 y( E    end
end
. A) p$ r' ~) S6 x
) v6 V0 y* S2 {; cassignin('base','X1',R1)
assignin('base','Z1',Z1)
" o6 S( G8 _! `% o$ qassignin('base','C1',theta1)
# z& i! B. N/ `6 d5 e  lassignin('base','X2',R2)
  Z; V8 e7 [5 z, K9 L4 Vassignin('base','Z2',Z2)
assignin('base','C2',theta2)
+ a) I. g3 h- m
& u* j; Z" l" A' [5 a% F# j%G代碼的輸出
% fid=fopen('turncode.txt','wt');
% fprintf(fid,'O0001\n');
# l* M; R: I8 P7 @& q% for  i=1:length(R)
%     for j=1:length(phi0)
; _+ n) w' Y1 \2 y3 V5 s%         fprintf(fid,'X=%0.4f\t',R1(i,j))
  B1 Z# k2 s5 D. K8 I2 k% X% z%         fprintf(fid,'Z=%0.4f\t',Z1(i,j))
, A& U9 g+ I1 ]3 R1 j" G. L%         fprintf(fid,'C=%0.4f\n',theta1(i,j))
. Y0 w7 S" v' u1 N4 P6 \4 |; h%     end
6 {7 W% B. q7 o( U( U' R( |- Y$ \% end
( X* J2 ?: f8 x! Q) z. ]8 J0 Y

% q- V) ]/ ~& K

function theta=cal_theta(x1,y1)
if (x1>0)&&(y1<0)
    theta=acos(x1./(x1^2+y1^2)^0.5);%電機(jī)轉(zhuǎn)是角度還是弧度？
elseif (x1<0)&&(y1<0);
6 p( R! d1 ?5 y; T* ytheta=acos(x1./(x1^2+y1^2)^0.5);
& ~% u3 s0 M; {1 T6 Eelseif (x1<0)&&(y1>0);
$ s4 @% M  X$ l  Y6 h0 e. ?; |    theta=2*pi-acos(x1./(x1^2+y1^2)^0.5);
6 j' r2 e. j, E( k6 Telseif (x1>0)&&(y1>0);
4 Q1 {% I$ t0 Z8 m, l% S' a    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);
) Q' j% g& g# ?0 x% H# C: _1 i5 }elseif (x1>0)&&(y1==0);
, ~3 Q( n) c8 \0 f: r) r    theta=0;
elseif (x1<0)&&(y1==0);
1 B* e8 A/ O) j0 w: ?0 y    theta=pi;
7 m+ n0 h. D- I& @elseif (x1==0)&&(y1<0);
* w3 _" ]& R9 t9 k$ Q! a' T# G    theta=pi/2;
, B! q: u5 R7 ^. \( Kelseif (x1==0)&&(y1>0);
    theta=3*pi/2;
end

東海fyh126

我也看不懂，，，，，

jxcrxssh

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

car1982

不懂

' x6 s9 b! K( ?5 U