基于Pro/ENGINEER的渐开线变位直齿圆柱齿轮实体模型的二次开发

sunshine1026

Pro/ENGINEER是由美国PTC公司开发的三维CAD软件，在工业造型设计和模具设计制造行业中得到广泛应用。运用Pro/ENGINEER软件，设计人员可以在计算机上实现虚拟现实设计和制造。 在Pro/ENGINEER软件中，Program模块是自动化零件与组合件设计的一项重要工具，使用者可以经由非常简易且高级的程序语言来控制特征的出现与否、尺寸的大小、零件的出现与否、零件的个数等。当零件或组合件的Program设计完成后，往后读取此零件或组合件时，其各种变化情况即可利用问答式的方式得到不同的形状，以实现产品设计的要求。
$ W# g1 W# u: r- C5 y齿轮虽然在机械行业中有广泛的使用，但Pro/ENGINEER软件不能直接生成。通过复杂的造型设计虽可以生成齿轮，但变化齿数、模数和变位系数后，又需进行复杂的设计计算和造型形状判断，工作量很大，齿轮的三维造型设计成为造型设计中望而却步的工作。本文基于渐开线齿轮的生成原理，结合Program程序，研制出变位直齿圆柱齿轮实体造型的自动化设计程序，用户只要输入齿轮的设计参数（齿数z1、z2，模数m，变位系数x1、x2，齿槽圆角半径、齿宽、轴孔直径、键槽高度和宽度）就可以快速地生成齿轮实体，使用简单方便，该设计软件可以广泛地应用于网页设计和机械产品造型设计中。
2 齿轮齿槽形状的精确确定
齿形的渐开线形状可通过高级Sweep特征设计精确地设计出来，但齿轮齿槽的形状与齿轮齿数、模数、变位系数和压力角有关。具体地说，下列两种情况决定齿槽的形状：
6 R. R, f6 p2 |; I! C

& C$ ^8 S4 @6 F* n; | i7 z7 _2 w

图1 第一种情况的齿槽形状

图2 第二种情况的齿槽形状

即：齿根圆的半径小于渐开线的基圆半径时和齿根圆的半径大于渐开线的基圆半径时。为了叙述方便，下文中，分别称前者为a情况，后者为b情况。因齿轮的啮合只在渐开线部分，因此对于a情况，设计时可在基圆的以下部分加与渐开线相切的直线以连接齿根圆，如图1所示。而对于b情况则无需加相切直线，仅截取渐开线的一部分即可，如图2所示。
+ y Z- D3 I1 \1 [$ R3 程序设计时需解决的几个关键问题
) I3 a4 R9 {2 H3 m8 \6 u3.1 啮合角的计算
) D) h2 D& o: x" g- i5 J5 h) V变位齿轮的齿顶圆直径和啮合角有关，因此设计时必须先由输入的齿轮设计参数经过复杂的数学计算得到啮合角值。
我们知道，外齿合时的啮合角计算公式为：

! R9 d0 D0 ]* t: v4 h) d) Q8 |9 J可见无法由显式直接表达出来，因此也就无法直接计算出来。虽然用C语言可以数值计算出，但Program自带的程序不支持C语言，幸运的是利用Program自带的解方程程序可以解出，具体的设计程序是：
8 C9 O3 j1 a$ Y+ q$ s: T8 M
上式中必须赋初值，否则数值解会不确定，考虑到啮合角一般大于10度，所以初值设为10即可。由于Program程序中三角函数内的角度是度，因此公式右边最后的应处理成弧度，以便和tan同量纲。
' v4 T) ?2 f' t+ @* @& P3.2 a种情况和b种情况齿轮的调用
对于Pro/ENGINEER软件，在零件设计模式下是无法实现两种不同零件的调用的，但在装配模式下可以实现这种调用功能，具体的方法是：在.asm下新建文件，坐标系统为mmns，把a、b两种情况的齿轮造型在基准缺省方式下直接装配在一起，这时的装配件是a和b两种情况齿轮的重叠，Program编程时，通过判断齿轮的齿根圆半径是否大于渐开线的基圆半径来决定采用a种情况还是b种情况的齿轮。若齿根圆半径小于基圆半径，则采用a种情况的齿轮 ；否则，则采用b种情况的齿轮。
4 齿轮的自动化设计程序流程图
4.1 a、b两种情况时的齿轮自动化设计流程图
图3是a种情况时齿轮的自动化设计流程图。b种情况的齿轮自动化设计流程和上图类似，只是齿槽设计时改为b种情况即可。
/ J+ U; B/ O% Q9 c, a4.2 综合a种和b种两种情况时的齿轮自动化设计流程图
综合a和b情况后的自动化设计流程图如图4所示。
$ \* K* ]6 i# x, |( Z

图3 a种情况的齿轮自动化设计流程

图4 综合a、b两种情况后的齿轮自动化设计流程

5 源程序
# \4 D# K- _) p4 F5.1 a种情况的齿轮自动化设计Program原程序
VERSION
REVNUM 23833
LISTING FOR PART GEAR_A
0 J- w Q `: w$ k5 Y$ S- P/ e: JINPUT
1 L* g' O j/ @2 W J: |& ^Z1 NUMBER /*齿轮1齿数
"Enter Z1: "
8 k7 z Y, V' v# r" h. R3 fZ2 NUMBER /*齿轮2齿数
"Enter Z2: "
M NUMBER /*模数
+ E' ]6 r. u) d. _5 @"Enter M: "
5 Y, }6 K7 j+ n/ w5 e1 A( bX1 NUMBER /*齿轮1变位系数
" r5 u* H. c/ U$ N"Enter x1:"
X2 NUMBER /*齿轮2变位系数
. G# ^& }' p4 D. h$ z0 o: B% k"Enter x2:"
ALPHA NUMBER /*压力角
"Enter the pressure angle: "
# p4 l) `$ W: \6 u+ o7 U8 H% o# WGEAR_WIDTH NUMBER /*齿轮宽度
"Enter the gear width: "
RAD_FILLET NUMBER /*齿槽圆角半径
"Enter rad_fillet:"
" n5 S8 g8 \6 W" O# l# T* HSHAFT_DIA NUMBER /*轴孔直径
( \7 j: h: Q) a5 V% `/ O' O"Enter shaft_dia:"
KEY_WIDTH NUMBER /*键槽宽度
"Enter key_width:"
9 |& V, B$ o# [KEY_HEIGHT NUMBER /*键槽高度
"Enter key_height:"
: x5 A+ F l+ p$ ZEND INPUT
RELATIONS
INVA=TAN(ALPHA)-ALPHA*PI/180
/*标准压力角时的渐开线函数值
6 r: U! h: O. c Y ~INVAP=2*(X1+X2)*TAN(ALPHA)/(Z1+Z2)
, v5 H6 K$ [/ j1 v% c　　　　+INVA
　　　/*计算啮合角时的渐开线函数值
7 [3 X* e6 {: O. x. m0 L! Y0 JAP=10
SOLVE /*解方程程序
INVAP=TAN(AP)-AP*PI/180
& ?7 C1 I5 j& D& YFOR AP
& m) i9 \1 L) V+ ~* d: `( b5 ]4 vHA=1 /*齿顶高系数
/ h# j1 q4 c/ g' z& NIF M<1 /*径向间隙系数的确定
C=0.35
ELSE
C=0.25
ENDIF
% F$ |$ Y/ l$ Z0 n7 r$ ~. ^R1 = 0.5*M*Z1 /*齿轮1分度圆半径
R2 = 0.5*M*Z2 /*齿轮1分度圆半径
D = R1+R2 /*标准中心距
RD = D*COS(ALPHA)/COS(AP)
/*实际中心距
) ^1 `0 k C5 r6 ^- nRB1 = R1*COS(ALPHA) /*基圆半径
RF1 = R1-M*(HA+C-X1) /*齿根圆半径
RA1 = RD-R2+M*(HA-X2) /*齿顶圆半径
5 V% Y7 e- J9 D, Y( H, _THICK1 = M*(0.5*PI+2*X1*TAN(ALPHA))
/*分度圆上齿厚（弧长）
" ]: _ h$ s0 Q' K+ DTOOTH_THICK1=THICK1*COS(ALPHA)
　　+2*RB1*INVA /*基圆上齿厚（弧长）
9 `# H# W5 ~& a) yANG_TOOTH_THICK1=TOOTH_THICK1/RB1*180/PI /*基圆上齿厚角（度）
7 f* }, z; L' \& K- pD0=RA1 /*内部参数赋值
~8 b( t$ s& ]- |D1=GEAR_WIDTH
, ?, d6 W( V b9 u! @% l* WD2=360/Z1
D3=ANG_TOOTH_THICK1
D4=180/Z1
D5=RB1
D13=RF1
D14=RAD_FILLET
D15=360/Z1
6 Z! l- U9 w3 ^( A$ y$ y( i2 HP16=Z1
) |& z( f( {) L* d) LD36=SHAFT_DIA
D37=KEY_HEIGHT
D38=KEY_WIDTH
. }5 _9 q n5 g, W% K! B7 CEND RELATIONS
…………………………
5.2 b种情况的齿轮自动化设计Program原程序
+ L: N+ s9 m$ c$ n% n. LVERSION
REVNUM 60662
LISTING FOR PART GEAR_B
INPUT
9 N! ~% Y4 m: g, G; \9 IZ1 NUMBER /*齿轮1齿数
* R/ i0 E9 Y9 o1 z) \# J"Enter Z1: "
+ Z2 i5 G! J; c& V3 {5 NZ2 NUMBER /*齿轮2齿数
"Enter Z2: "
! g1 {8 h: x; ]! qM NUMBER /*模数
"Enter M: "
X1 NUMBER /*齿轮1变位系数
) K& q8 @8 ?$ P4 u ?+ ^"Enter x1:"
X2 NUMBER /*齿轮2变位系数
: s' B7 M- U! P) V5 ^" @"Enter x2:"
) v0 Q* Z" P) M% E h' AALPHA NUMBER /*压力角
"Enter the pressure angle: "
$ j5 v- S# _! C: d: vGEAR_WIDTH NUMBER /*齿轮宽度
"Enter the gear width: "
RAD_FILLET NUMBER /*齿槽圆角半径
5 z+ d7 f9 P9 g2 H" b1 }, X"Enter rad_fillet:"
4 I/ u8 x! o3 s4 h- i: lSHAFT_DIA NUMBER /*轴孔直径
3 v7 G* X1 J! J3 j9 O5 d"Enter shaft_dia:"
KEY_WIDTH NUMBER /*键槽宽度
. a( w+ G) l2 n( C- `"Enter key_width:"
5 ^- z0 `% V2 C" X5 q. GKEY_HEIGHT NUMBER /*键槽高度
"Enter key_height:"
END INPUT
, F* P# P9 t9 GRELATIONS
* t6 X8 ~0 x- k" JINVA=TAN(ALPHA)-ALPHA*PI/180
9 w$ n. O$ v' T0 @) _/*标准压力角时的渐开线函数值
M. Y9 v* [: c: \' W0 G6 j8 AINVAP=2*(X1+X2)*TAN(ALPHA)/(Z1+Z2)
　　　　+INVA
; `& A9 l* D8 S( y/ _8 l　　　/*计算啮合角时的渐开线函数值
. ]) r" z% L/ i8 V5 s; ?AP=10
SOLVE /*解方程程序
( B# H; G. P7 b3 z. bINVAP=TAN(AP)-AP*PI/180
' b) r: G V! W' A5 r3 VFOR AP
HA=1 /*齿顶高系数
IF M<1 /*径向间隙系数的确定
C=0.35
ELSE
C=0.25
ENDIF
' O4 l% O' W: Q+ W1 A3 GR1 = 0.5*M*Z1 /*齿轮1分度圆半径
/ y5 |5 w! o% E, P4 YR2 = 0.5*M*Z2 /*齿轮1分度圆半径
D = R1+R2 /*标准中心距
RD = D*COS(ALPHA)/COS(AP)
. |! o0 s+ {- r2 F5 S* z/*实际中心距
RB1 = R1*COS(ALPHA) /*基圆半径
2 v, d% |& v- t5 [. qRF1 = R1-M*(HA+C-X1) /*齿根圆半径
RA1 = RD-R2+M*(HA-X2) /*齿顶圆半径
: ], i; Y- s. H) t7 ?, i3 C O& STHICK1 = M*(0.5*PI+2*X1*TAN(ALPHA))
/*分度圆上齿厚（弧长）
0 w! B/ c, q- J+ e# B/ K5 bTOOTH_THICK1=THICK1*COS(ALPHA)
1 M2 i' P7 K- ^) [; H+ t　　+2*RB1*INVA /*基圆上齿厚（弧长）
/ V8 a' E D# b8 }* s7 J& Y) t& f gANG_TOOTH_THICK1=TOOTH_THICK1/RB1*180/PI /*基圆上齿厚角（度）
: }9 ]* R! K1 AD0=RA1 /*内部参数赋值
) p9 H9 L2 C n" w& I1 {3 [D1=GEAR_WIDTH
0 |( \! J2 v$ |% p( @ pD2=360/Z1
, e; P$ F+ v) E2 R& {( K Z) m3 BD3=ANG_TOOTH_THICK1
% g. Y0 N8 V1 c* [5 z! uD4=180/Z1
D5=RB1
D13=RF1
& U9 a. O! O) eD14=RAD_FILLET
D15=360/Z1
P16=Z1
D116=SHAFT_DIA
, `/ J6 O4 l8 j# k0 tD117=KEY_HEIGHT
D118=KEY_WIDTH
END RELATIONS
# e. {. K" i7 o/ y& v1 d……………………