基于Pro/ENGINEER的漸開線變位直齒圓柱齒輪實體模型的二次開發(fā)

sunshine1026

Pro/ENGINEER是由美國PTC公司開發(fā)的三維CAD軟件，在工業(yè)造型設(shè)計和模具設(shè)計制造行業(yè)中得到廣泛應(yīng)用,。運用Pro/ENGINEER軟件,，設(shè)計人員可以在計算機上實現(xiàn)虛擬現(xiàn)實設(shè)計和制造。    在Pro/ENGINEER軟件中,，Program模塊是自動化零件與組合件設(shè)計的一項重要工具,，使用者可以經(jīng)由非常簡易且高級的程序語言來控制特征的出現(xiàn)與否、尺寸的大小,、零件的出現(xiàn)與否,、零件的個數(shù)等。當(dāng)零件或組合件的Program設(shè)計完成后,，往后讀取此零件或組合件時,，其各種變化情況即可利用問答式的方式得到不同的形狀，以實現(xiàn)產(chǎn)品設(shè)計的要求,。
    齒輪雖然在機械行業(yè)中有廣泛的使用,，但Pro/ENGINEER軟件不能直接生成。通過復(fù)雜的造型設(shè)計雖可以生成齒輪,，但變化齒數(shù),、模數(shù)和變位系數(shù)后，又需進行復(fù)雜的設(shè)計計算和造型形狀判斷,，工作量很大,，齒輪的三維造型設(shè)計成為造型設(shè)計中望而卻步的工作。本文基于漸開線齒輪的生成原理,，結(jié)合Program程序,，研制出變位直齒圓柱齒輪實體造型的自動化設(shè)計程序，用戶只要輸入齒輪的設(shè)計參數(shù)（齒數(shù)z1,、z2,，模數(shù)m，變位系數(shù)x1,、x2,，齒槽圓角半徑、齒寬,、軸孔直徑,、鍵槽高度和寬度）就可以快速地生成齒輪實體，使用簡單方便,，該設(shè)計軟件可以廣泛地應(yīng)用于網(wǎng)頁設(shè)計和機械產(chǎn)品造型設(shè)計中,。
2 齒輪齒槽形狀的精確確定
: w( W; m) O6 |$ Q    齒形的漸開線形狀可通過高級Sweep特征設(shè)計精確地設(shè)計出來，但齒輪齒槽的形狀與齒輪齒數(shù),、模數(shù),、變位系數(shù)和壓力角有關(guān)。具體地說,，下列兩種情況決定齒槽的形狀：
  
6 S) ]: N, V6 K' v& F4 ]) M

圖1 第一種情況的齒槽形狀

圖2 第二種情況的齒槽形狀

    即：齒根圓的半徑小于漸開線的基圓半徑時和齒根圓的半徑大于漸開線的基圓半徑時,。為了敘述方便，下文中,，分別稱前者為a情況,，后者為b情況。因齒輪的嚙合只在漸開線部分,，因此對于a情況,，設(shè)計時可在基圓的以下部分加與漸開線相切的直線以連接齒根圓，如圖1所示,。而對于b情況則無需加相切直線,，僅截取漸開線的一部分即可，如圖2所示,。        
3 程序設(shè)計時需解決的幾個關(guān)鍵問題
    3.1 嚙合角的計算
7 C! ?% t/ O  |8 r  {    變位齒輪的齒頂圓直徑和嚙合角有關(guān),，因此設(shè)計時必須先由輸入的齒輪設(shè)計參數(shù)經(jīng)過復(fù)雜的數(shù)學(xué)計算得到嚙合角值。
    我們知道,，外齒合時的嚙合角計算公式為：

    可見無法由顯式直接表達出來,，因此也就無法直接計算出來。雖然用C語言可以數(shù)值計算出,，但Program自帶的程序不支持C語言,，幸運的是利用Program自帶的解方程程序可以解出，具體的設(shè)計程序是：

& w# ?# H& ~- B6 o9 P    上式中必須賦初值,，否則數(shù)值解會不確定,，考慮到嚙合角一般大于10度，所以初值設(shè)為10即可,。由于Program程序中三角函數(shù)內(nèi)的角度是度,，因此公式右邊最后的應(yīng)處理成弧度,，以便和tan同量綱。
    3.2  a種情況和b種情況齒輪的調(diào)用
  V$ d, J6 w2 |* O8 A0 f; l7 w    對于Pro/ENGINEER軟件,，在零件設(shè)計模式下是無法實現(xiàn)兩種不同零件的調(diào)用的,，但在裝配模式下可以實現(xiàn)這種調(diào)用功能，具體的方法是：在.asm下新建文件,，坐標(biāo)系統(tǒng)為mmns,，把a、b兩種情況的齒輪造型在基準(zhǔn)缺省方式下直接裝配在一起,，這時的裝配件是a和b兩種情況齒輪的重疊,，Program編程時，通過判斷齒輪的齒根圓半徑是否大于漸開線的基圓半徑來決定采用a種情況還是b種情況的齒輪,。若齒根圓半徑小于基圓半徑,，則采用a種情況的齒輪 ；否則,，則采用b種情況的齒輪,。
  K) L2 m) g) c; @8 {) Y4 齒輪的自動化設(shè)計程序流程圖
) Z+ W" `9 ~6 A    4.1 a、b兩種情況時的齒輪自動化設(shè)計流程圖
4 c1 }+ C( Z7 X: m. B* n    圖3是a種情況時齒輪的自動化設(shè)計流程圖,。b種情況的齒輪自動化設(shè)計流程和上圖類似,，只是齒槽設(shè)計時改為b種情況即可。
    4.2 綜合a種和b種兩種情況時的齒輪自動化設(shè)計流程圖
. K  N# F8 j& d( `  s" z4 X    綜合a和b情況后的自動化設(shè)計流程圖如圖4所示,。
" \8 F) b  D( ^3 j5 {  A% [

圖3 a種情況的齒輪自動化設(shè)計流程

圖4 綜合a,、b兩種情況后的齒輪自動化設(shè)計流程

5 源程序
, J2 e6 \' f, c9 y" F" p+ ~    5.1 a種情況的齒輪自動化設(shè)計Program原程序
VERSION
REVNUM 23833
LISTING FOR PART GEAR_A
" S/ m6 d* o; d/ V: }INPUT
! ^' g2 t4 x+ W0 L; o" m, iZ1 NUMBER    /*齒輪1齒數(shù)
"Enter Z1: "
0 }' J) v5 e+ b# k5 l+ f! T9 mZ2 NUMBER    /*齒輪2齒數(shù)
. z! A8 n# F( e0 }$ Q3 }) d  t"Enter Z2: "
M NUMBER    /*模數(shù)
"Enter M: "
+ U  c5 D5 E4 V. VX1 NUMBER   /*齒輪1變位系數(shù)
"Enter x1:"
X2 NUMBER   /*齒輪2變位系數(shù)
"Enter x2:"
ALPHA NUMBER  /*壓力角
9 ]# q# U) `# T"Enter the pressure angle: "
  G1 Q5 Y) W" z9 |: ~5 m! F: `7 K. `GEAR_WIDTH NUMBER  /*齒輪寬度
"Enter the gear width: "
+ j+ m3 ?3 V" `  c& K, qRAD_FILLET NUMBER   /*齒槽圓角半徑
7 s" z6 ?: g* `1 t( H; o- R"Enter rad_fillet:"
SHAFT_DIA NUMBER    /*軸孔直徑
"Enter shaft_dia:"
) J1 S: m* Z, r7 N3 z$ T& E) x, QKEY_WIDTH NUMBER   /*鍵槽寬度
"Enter key_width:"
KEY_HEIGHT NUMBER  /*鍵槽高度
"Enter key_height:"
END INPUT
RELATIONS
INVA=TAN(ALPHA)-ALPHA*PI/180  
      /*標(biāo)準(zhǔn)壓力角時的漸開線函數(shù)值
" S% n/ Z/ Y* K% g# t& c# xINVAP=2*(X1+X2)*TAN(ALPHA)/(Z1+Z2)
　　　　+INVA   
/ Q0 U5 C" N; G6 v& J& O4 u　　　/*計算嚙合角時的漸開線函數(shù)值
AP=10
SOLVE         /*解方程程序
# u4 L4 d* ?9 j; ~+ S  INVAP=TAN(AP)-AP*PI/180
FOR AP        
HA=1           /*齒頂高系數(shù)
- o& g( I' U3 {3 `+ e) VIF M<1          /*徑向間隙系數(shù)的確定
   C=0.35
ELSE
   C=0.25
) f" h% |/ W* EENDIF         
R1 = 0.5*M*Z1   /*齒輪1分度圓半徑
$ K' J; u  \3 C" _R2 = 0.5*M*Z2   /*齒輪1分度圓半徑
D = R1+R2      /*標(biāo)準(zhǔn)中心距
9 ^& D) D# b1 f, J8 F! GRD = D*COS(ALPHA)/COS(AP)
            /*實際中心距
RB1 = R1*COS(ALPHA)   /*基圓半徑
8 X0 ?" l9 z7 f3 D6 ~1 y& CRF1 = R1-M*(HA+C-X1)   /*齒根圓半徑
RA1 = RD-R2+M*(HA-X2)  /*齒頂圓半徑
THICK1 = M*(0.5*PI+2*X1*TAN(ALPHA))
              /*分度圓上齒厚（弧長）
$ s" Q3 P4 j1 l: E# e% |2 ]TOOTH_THICK1=THICK1*COS(ALPHA)
- X8 M- H. {% w3 ?' B3 H9 g' F  O　　+2*RB1*INVA  /*基圓上齒厚（弧長）
ANG_TOOTH_THICK1=TOOTH_THICK1/RB1*180/PI        /*基圓上齒厚角（度）
D0=RA1           /*內(nèi)部參數(shù)賦值
# {, t6 }$ ~) M, TD1=GEAR_WIDTH
D2=360/Z1
D3=ANG_TOOTH_THICK1
D4=180/Z1
; a: P6 N' g: g5 eD5=RB1
9 @5 ~1 g+ N: H1 d4 C) G. R2 lD13=RF1
& a* W% `; o5 I* hD14=RAD_FILLET
# X# n" l1 z* S+ C. PD15=360/Z1
P16=Z1
D36=SHAFT_DIA
D37=KEY_HEIGHT
D38=KEY_WIDTH
END RELATIONS
…………………………
5.2  b種情況的齒輪自動化設(shè)計Program原程序
VERSION
REVNUM 60662
LISTING FOR PART GEAR_B
INPUT
Z1 NUMBER    /*齒輪1齒數(shù)
" h# l, N" G; c9 r2 T0 V5 ]/ d3 s "Enter Z1: "
Z2 NUMBER    /*齒輪2齒數(shù)
"Enter Z2: "
M NUMBER    /*模數(shù)
6 C2 I' D( @+ z/ ` "Enter M: "
X1 NUMBER   /*齒輪1變位系數(shù)
1 E+ w! g& d, H( {% w  N3 T/ m( C "Enter x1:"
' t( d2 r% x! R6 s: n9 y( ?' g X2 NUMBER   /*齒輪2變位系數(shù)
"Enter x2:"
5 c# Y2 j" B7 j8 O ALPHA NUMBER  /*壓力角
"Enter the pressure angle: "
& C. T6 l2 I6 |9 J+ M9 l1 Y: p GEAR_WIDTH NUMBER  /*齒輪寬度
8 n8 _8 F7 z. U; V4 J "Enter the gear width: "
RAD_FILLET NUMBER   /*齒槽圓角半徑
+ x' x1 H9 A; }" H3 @# n1 ~ "Enter rad_fillet:"
5 ?  p. b5 Y! B- c% t. v SHAFT_DIA NUMBER    /*軸孔直徑
2 J4 y. T3 K& w5 e: I "Enter shaft_dia:"
+ P- ^9 l0 L& V4 u; v' l$ ]% o KEY_WIDTH NUMBER   /*鍵槽寬度
1 u. b+ \3 d" U "Enter key_width:"
" D/ z  g4 Q- t- t5 a4 o+ Y- @ KEY_HEIGHT NUMBER  /*鍵槽高度
5 q: K, L/ b' U "Enter key_height:"
9 I( O( O% S# ?  l2 U" I& ~END INPUT
RELATIONS
* d0 {3 A$ m& }5 }% @4 e9 T9 iINVA=TAN(ALPHA)-ALPHA*PI/180  
      /*標(biāo)準(zhǔn)壓力角時的漸開線函數(shù)值
2 m) ^1 Z! G: v* P$ E) A3 rINVAP=2*(X1+X2)*TAN(ALPHA)/(Z1+Z2)
　　　　+INVA   
　　　/*計算嚙合角時的漸開線函數(shù)值
AP=10
" B# c( }( ?( e( G( ?SOLVE         /*解方程程序
  INVAP=TAN(AP)-AP*PI/180
FOR AP         
. [. o" f  t5 I5 h0 U, HHA=1           /*齒頂高系數(shù)
' e/ \5 }( V9 c# J5 `7 MIF M<1          /*徑向間隙系數(shù)的確定
   C=0.35
# Z! J- [: M$ z, a1 tELSE
. v  f: w7 Y, u$ A9 ~   C=0.25
ENDIF         
R1 = 0.5*M*Z1   /*齒輪1分度圓半徑
% J" |# b8 T9 Y" UR2 = 0.5*M*Z2   /*齒輪1分度圓半徑
  Z# g/ ~4 e* x" h& D; t D = R1+R2      /*標(biāo)準(zhǔn)中心距
5 \+ Y/ S. V. A* Y8 |RD = D*COS(ALPHA)/COS(AP)
/ h4 G& t7 J$ \8 ?; c               /*實際中心距
: R- n& q9 ?, [2 `6 N2 xRB1 = R1*COS(ALPHA)   /*基圓半徑
6 y7 @2 [, P% HRF1 = R1-M*(HA+C-X1)   /*齒根圓半徑
1 n  E" q4 b' P2 u* Z- {RA1 = RD-R2+M*(HA-X2)  /*齒頂圓半徑
THICK1 = M*(0.5*PI+2*X1*TAN(ALPHA))
              /*分度圓上齒厚（弧長）
TOOTH_THICK1=THICK1*COS(ALPHA)
　　+2*RB1*INVA  /*基圓上齒厚（弧長）
+ s, x9 O( J1 ~$ A% Z+ _% |ANG_TOOTH_THICK1=TOOTH_THICK1/RB1*180/PI        /*基圓上齒厚角（度）
D0=RA1           /*內(nèi)部參數(shù)賦值
% s  r' ~/ r+ OD1=GEAR_WIDTH
- w% s1 ~. M% L- p$ Y- ]D2=360/Z1
1 I  m" k  S7 D) A5 ]2 q8 W. `D3=ANG_TOOTH_THICK1
D4=180/Z1
) l; y1 Q: O& B/ f0 rD5=RB1
9 H, g" V+ t6 V5 I7 m8 kD13=RF1
D14=RAD_FILLET
D15=360/Z1
' ?- O( j6 @5 C/ `P16=Z1
' ]) ^8 L9 `: A# H5 ZD116=SHAFT_DIA
5 d. C/ x* p+ [% N% hD117=KEY_HEIGHT
D118=KEY_WIDTH
1 V  ^& {1 }# [END RELATIONS
5 o( O4 a& d# s……………………