喬培君,金 秋
(1.武漢理工大學(xué) 汽車工程學(xué)院,湖北 武漢 430070;2.武漢理工大學(xué) 材料科學(xué)與工程學(xué)院,湖北 武漢430070;3.現(xiàn)代汽車零部件技術(shù)湖北省重點實驗室,湖北 武漢 430070)
零度弧齒錐齒輪是中點螺旋角為零度的弧齒錐齒輪,用于相交軸之間傳遞運動,具有傳動效率高、傳遞平穩(wěn)、工作可靠、噪音小等優(yōu)點,在輪船、汽車、礦山、農(nóng)機等方面發(fā)揮著極其重要的作用[1]。由于其齒面形狀是不規(guī)則的圓弧形、錐狀,因此其幾何建模比較困難,不利于其優(yōu)化設(shè)計與加工制造的實現(xiàn)。目前,零度錐齒輪實體建模方法有以下兩種:一種是基于球面漸開線形成原理,利用三維造型軟件進行實體建模[2],該方法將輪齒進行參數(shù)化設(shè)計,沒有涉及嚙合原理、機床加工體系,不利于齒形優(yōu)化設(shè)計與齒輪機床加工制造的進行;另一種方法是以空間嚙合原理與機床加工體系為基礎(chǔ)[3-4],利用微分幾何進行矩陣變換,得到齒面方程,由數(shù)據(jù)點擬合成齒面模型,該方法得到的齒面更加光滑精確。第二種方法通過將弧齒面在機床平面上進行旋轉(zhuǎn)投影,選取由頂錐面、根錐面、外錐面、內(nèi)錐面的投影線作為邊界構(gòu)成的區(qū)域來求解齒面點坐標,這種方法在確定區(qū)域邊界時較為復(fù)雜繁瑣,且存在誤差。筆者將第二種方法進行改進、延伸,通過設(shè)計全齒高等參數(shù)確定加工刀盤的刀刃高度和刀具參數(shù)t的范圍;依據(jù)產(chǎn)形原理,在“刀盤-搖臺坐標系轉(zhuǎn)換”示意圖中確定刀盤轉(zhuǎn)角φm的范圍,以此來確定邊界條件,因此該方法更為精確,曲面參數(shù)范圍的選取更為簡潔,建模過程的控制更為嚴格。根據(jù)銑齒加工原理,由刀盤軸截面方程及截面參數(shù)推導(dǎo)出零度錐齒輪的大小輪齒面方程,在確定刀盤轉(zhuǎn)角范圍的前提下,運用MATLAB計算齒坯凹凸面的點坐標,在Pro/E里建立了零度錐齒輪的精確幾何模型[5]。
零度弧齒錐齒輪是在銑齒機上加工的,如圖1所示,在機床的搖臺機構(gòu)上可以模擬出一個假想齒輪,某一時刻刀盤的一個旋轉(zhuǎn)錐面相當于假想齒輪的一個齒面。工件與假想齒輪以某傳動關(guān)系各自旋轉(zhuǎn)同時進行嚙合,刀盤此時在工件輪坯上切出一個齒槽。如此運動,刀盤就能加工出來輪齒齒面。這種加工方法稱為展成法,該假想齒輪稱為產(chǎn)形輪[6-8]。在大輪與小輪嚙合部位之間模擬出一個平面產(chǎn)形輪如圖2所示。當小輪與大輪以一定的傳動比做嚙合運動時,產(chǎn)形輪可以替代任意一個工件進行旋轉(zhuǎn)運動,而不影響兩個齒輪的整體嚙合效果。其中,大輪、小輪、產(chǎn)形輪共用一個錐頂O,重合部分是接觸點P,ω1、ω2、ωg分別為小輪、大輪和產(chǎn)形輪的轉(zhuǎn)速。
圖1 機床設(shè)計原理
圖2 產(chǎn)形輪展成原理
2.1.1 大輪刀盤切削錐面方程
圖3為加工零度弧齒錐齒輪大輪所用的雙面銑刀盤的截面示意圖,圖中,r2v為加工左旋大輪所用的刀盤半徑,α為刀面齒形角,值為20°。刀盤截面分為內(nèi)刀面、外刀面與齒槽底面,其中,內(nèi)刀面加工大輪的凸面,外刀面加工大輪的凹面。S1為大輪大端齒槽厚,ht為全齒高,ha2為大輪齒頂高,hf2為大輪齒根高,通過這些變量求出刀刃處刀頂距w2的值,t為刀盤齒形變量,可得刀盤截面方程為:
圖3 大輪所用刀盤軸截面
將刀盤截面繞Zt軸旋轉(zhuǎn)一周,得到刀盤切削錐面方程和刀盤切削曲面單位法矢量分別為:
(2)
(3)
2.1.2 大輪機床的坐標轉(zhuǎn)換關(guān)系
通過圖4機床坐標系的轉(zhuǎn)換位置圖可知,大輪的齒面方程與嚙合方程分別為:
r2b(φ2m,t,θ)=T2b2a(φ2)×T2a2m×
T2m20(φ2m)×T202×r2(t,θ)
(4)
n2b(φ2m,t,θ)=L2b2a(φ2)×L2a20×
L202m(φ2m)×L2m2×n2(t,θ)
(5)
圖4 加工大輪的機床轉(zhuǎn)換關(guān)系
式(4)中,T202為刀盤坐標系到搖臺坐標系的轉(zhuǎn)換矩陣;T2m20為搖臺坐標系到機床坐標系的轉(zhuǎn)換矩陣,T2a2m為機床坐標系到輔助坐標系的轉(zhuǎn)換矩陣,T2b2a為輔助坐標系到輪坯坐標系的轉(zhuǎn)換矩陣。式(5)中,Lij為各轉(zhuǎn)換矩陣提取前三行與前三列所構(gòu)成的子矩陣。其中,φ2=-i2×φ2m,i2為機床滾比,φ2m為產(chǎn)形輪的轉(zhuǎn)角。
大輪的嚙合方程為:
n2b(φ2m,t,θ)×V2b2m(φ2m,t,θ)=0
(6)
(7)
式(6)中,n2b為產(chǎn)形輪與大輪在嚙合點的公共單位法矢量;V2b2m為產(chǎn)形輪與大輪在嚙合點的相對運動速度[9];其求解方法見式(7),τ為時間變量。將式(4)與式(6)聯(lián)立,在MATLAB軟件中求解,可消除參數(shù)φ2m,得到由曲面參數(shù)t和θ表示的大輪齒面方程r2b(t,θ)。
小輪的齒面方程推導(dǎo)過程與大輪的齒面方程推導(dǎo)過程的原理是一樣的,但是小輪的加工情況比大輪復(fù)雜,右旋小輪使用單刀刀面法加工,即內(nèi)刀面與外刀面分別加工出小輪的凸面與凹面。加工小輪的刀盤軸截面如圖5所示,小輪的刀盤截面方程與切削錐面方程分別為:
(8)
(9)
圖5 小輪所用刀盤軸截面
那么小輪的齒面方程與嚙合方程為:
r1b(φ1m,t,θ)=T1b1a(φ1)×T1a1m×
T1m10(φ1m)×T101×r1(t,θ)
(10)
n1b(φ1m,t,θ)×V1b1m(φ1m,t,θ)=0
(11)
將式(10)與式(11)聯(lián)立,在MATLAB軟件中求解,可得到由參數(shù)和表示的小輪齒面方程。
大輪和小輪的主要參數(shù)與主要機床調(diào)整參數(shù)如表1和表2所示。
表1 齒輪副的基本參數(shù)
表2 機床調(diào)整參數(shù)
圖6 刀盤-搖臺坐標轉(zhuǎn)換關(guān)系
在MATLAB[10]里面編寫循環(huán)程序,設(shè)定t和θ的步長進行計算,得到齒面點在輪坯坐標系下的坐標值。
(1)以大輪為例將計算得到的點坐標值按照X,Y,Z排列成3列N行的形式,保存為Pro/E能識別的“ibl”文件。用記事本打開該文件進行以下的編輯操作:
open
arclength
begin section!1
begin curve!15.4608498e+01-2.8226705e-01 1.7203320e+02
5.8056706e+01-1.9763803e+00 1.7124564e+02
6.1504915e+01-3.6704936e+00 1.7045809e+02
……
begin curve!25.3586427e+01 -8.2940911e-01 1.6881336e+02
5.7023882e+01 -2.5292791e+00 1.6799194e+02
……
將該文件導(dǎo)入到Pro/E中,生成齒面網(wǎng)格線,通過邊界混合命令生成的輪齒齒面,最后形成應(yīng)當去除材料的部分曲面,編輯生成齒槽實體如圖7所示。
圖7 齒槽實體
(2)在Pro/E中建立齒坯實體模型。對齒坯模型進行去除材料操作,陣列形成大輪三維實體模型如圖8所示,用同樣的方法得到小輪的三維實體模型如圖9所示。
圖8 大輪實體模型
(3)將大小齒輪裝配后,在Pro/E里進行運動仿真,如圖10所示。
圖9 小輪實體模型
圖10 齒輪裝配體模型
(1)根據(jù)零度弧齒錐齒輪的銑齒加工原理,建立了齒輪加工的坐標轉(zhuǎn)換關(guān)系,由刀盤軸截面方程及截面參數(shù)推導(dǎo)出零度錐齒輪的齒面方程。
(2)基于展成法加工原理,確定銑刀盤轉(zhuǎn)角范圍,運用MATLAB編程計算齒面點的坐標值。
(3)將點坐標值以一定的格式導(dǎo)入到Pro/E中,生成齒廓曲線,通過邊界混合功能生成完整齒面,最后得到完整的實體模型,為零度弧齒錐齒輪的傳動分析提供了理論依據(jù)。在此基礎(chǔ)上,運用Pro/E建立了零度錐齒輪的大小齒輪模型,并進行裝配與運動仿真,為零度錐齒輪的傳動分析與加工制造提供理論依據(jù)。
[1] 曾韜.螺旋錐齒輪的設(shè)計與加工[M].哈爾濱:哈爾濱大學(xué)出版社,1989.
[2] 王犇.弧齒錐齒輪參數(shù)化設(shè)計及嚙合傳動規(guī)律研究[D].武漢:武漢理工大學(xué),2010.
[3] Litvin F L. Gear Geometry and Applied Theory[M]. New Jersey: Prentice Hall, 1994.
[4] 吳序堂.齒輪嚙合原理[M].西安:西安交通大學(xué)出版社,2009.
[5] 吳鑫輝.弧齒錐齒輪實體建模與齒面接觸分析[D].大連:大連交通大學(xué),2012.
[6] Zheng F Y,Hua L,Han X H. Linkage Model and Manufacturing Process of Shaping Non-circular Gears[J]. Mechanism and Machine Theory, 2016,96:192-212.
[7] Zheng F Y,Hua L,Han X H. Synthesis of Indexing Mechanisms with Non-circular Gears[J]. Mechanism and Machine Theory, 2016,105:108-128.
[8] Zheng F Y,Hua L,Han X H. The Mathematical Model and Mechanical Properties of Variable Center Distance Gears Based on Screw Theory[J]. Mechanism and Machine Theory, 2016, 101:116-139.
[9] Zheng F, Hua L, Chen D F, et al. Generation of Non-Circular Spiral Bevel Gears by Face-Milling Method[J]. Journal of Manufacturing Science & Engineering, 2016,138(8)198-210.
[10] 劉衛(wèi)國.MATLAB 程序設(shè)計與應(yīng)用[M].北京:高等教育出版社,2006.