于艷濤，王甲興，徐永福，崔鳳奎

1洛陽(yáng)格瑞德重型機(jī)械有限公司 河南洛陽(yáng) 471006

2第一拖拉機(jī)股份有限公司大拖公司 河南洛陽(yáng) 471003

3中信重工機(jī)械股份有限公司 河南洛陽(yáng) 471039

4河南科技大學(xué)機(jī)電工程學(xué)院 河南洛陽(yáng) 471003

花 鍵軸在裝備制造領(lǐng)域應(yīng)用廣泛，常用的加工方法包括銑削法、滾切法、磨削法、冷搓、冷擠壓及冷滾打等[1]。其中，高速冷滾打花鍵成形技術(shù)是利用常溫下金屬材料具有的塑性變形，通過(guò)高速旋轉(zhuǎn)的滾打輪周期性擊打毛坯軸，使其產(chǎn)生塑性凸起的齒形工藝過(guò)程[2]。冷滾打與傳統(tǒng)切削相比具有明顯的優(yōu)點(diǎn)：效率高，節(jié)省原材料，制造成本低，齒形表面完整性好，齒面強(qiáng)度與耐磨性顯著提高，產(chǎn)品使用壽命[3]。此外，冷滾打還克服了冷擠壓、冷搓等工藝只適用于小模數(shù)、短花鍵制造的缺點(diǎn)。隨著高速冷滾打花鍵加工技術(shù)逐步取代傳統(tǒng)加工方法，探討冷滾打花鍵成形機(jī)理有著重要的科學(xué)意義和工程應(yīng)用價(jià)值。

1 冷滾打花鍵制造工藝

兩組冷滾打輪分別安裝在關(guān)于花鍵對(duì)稱的兩套動(dòng)力軸上，動(dòng)力軸的間距可根據(jù)花鍵毛坯尺寸進(jìn)行調(diào)整。花鍵加工時(shí)，冷滾打輪繞動(dòng)力軸中心線高速旋轉(zhuǎn)，周期性擊打軸坯，使其外部金屬產(chǎn)生塑性流動(dòng)，在滾打輪外形限制下，逐步形成均勻的花鍵齒形[4]。當(dāng)滾打輪軋入軸坯時(shí)，兩者之間產(chǎn)生摩擦力，使?jié)L打輪發(fā)生自轉(zhuǎn)，這樣有助于延長(zhǎng)其使用壽命。連續(xù)分度冷滾打花鍵加工示意如圖 1 所示。

圖1 冷滾打花鍵加工示意Fig.1 Sketch of cold rolling spline machining

2 花鍵冷滾打成形數(shù)值仿真

2.1 花鍵毛坯尺寸的確定

冷滾打花鍵軸坯直徑的估算方法有：①以分度圓作為軸坯直徑；② 按照成形前后軸截面面積相等進(jìn)行計(jì)算；③按照齒形面積計(jì)算公式轉(zhuǎn)換；④ 通過(guò)反復(fù)試驗(yàn)確定合適的軸坯尺寸。因花鍵在滾打過(guò)程中為塑性變形，假設(shè)塑性變形時(shí)體積不可壓縮，且忽略滾打過(guò)程中金屬微量的軸向流動(dòng)，可認(rèn)為工件滾打前后的橫截面積不變，即工件被滾打輪齒頂擠壓凹陷部分的面積等于工件被擠壓漲出部分的面積。由此，軸坯直徑

式中：df為漸開(kāi)線花鍵的齒根圓直徑，mm；Z為漸開(kāi)線花鍵的齒數(shù)；rb為漸開(kāi)線基圓半徑，mm；α a為花鍵齒頂圓壓力角，(°)；αf為花鍵齒根圓壓力角，(°)；ra為花鍵齒頂圓半徑，mm；θn為漸開(kāi)線與分度圓交點(diǎn)處的展角，(°)；θa為漸開(kāi)線與齒頂圓交點(diǎn)處的展角，(°)；rf為花鍵齒根圓半徑，mm；θf(wàn)為漸開(kāi)線與齒報(bào)圓交點(diǎn)處的展角，(°)。

2.2 工藝參數(shù)及運(yùn)動(dòng)設(shè)定

以東方紅拖拉機(jī)半軸花鍵加工為例，運(yùn)動(dòng)特性按照實(shí)際加工進(jìn)行設(shè)定。冷滾打輪繞其回轉(zhuǎn)中心以2 000 r/min 高速旋轉(zhuǎn)，放開(kāi)其沿自身軸線轉(zhuǎn)動(dòng)方向的自由度，以便在滾打摩擦力作用下可自由旋轉(zhuǎn)；工件繞其軸線以 125 r/min 連續(xù)旋轉(zhuǎn)，并且沿軸以 90 mm/min 速度進(jìn)給，其工藝參數(shù)如表 1 所列。

表1 半軸花鍵加工工藝參數(shù)Tab.1 Parameters of half-shaft spline machining

2.3 DEFORM 成形仿真流程及流動(dòng)應(yīng)力模型

DEFORM-3D 軟件的模塊結(jié)構(gòu)是由前處理器、模擬處理器和后處理器三大模塊組成。成形仿真系統(tǒng)的建立是將彈塑性有限元理論、剛塑性成形工藝學(xué)、計(jì)算機(jī)圖形處理技術(shù)等相關(guān)理論和技術(shù)進(jìn)行有機(jī)結(jié)合的過(guò)程[5]，其仿真流程如圖 2 所示。

圖2 DEFORM 成形仿真流程Fig.2 Process flow of DEFORM formation simulation

DEFORM 軟件中有能量定律的經(jīng)驗(yàn)公式和材料實(shí)測(cè)數(shù)據(jù) 2 種流動(dòng)應(yīng)力模型。其中，后者根據(jù)試驗(yàn)中實(shí)測(cè)數(shù)據(jù)得到，由于它遵循材料的真實(shí)變化規(guī)律，能夠達(dá)到較高的精度。對(duì)于實(shí)測(cè)數(shù)據(jù)模型，如果在模擬過(guò)程中，材料的狀況超出了應(yīng)變、應(yīng)變速率或溫度定義的范圍時(shí)，程序?qū)⒆詣?dòng)根據(jù)最近的兩點(diǎn)數(shù)據(jù)進(jìn)行插值計(jì)算，得到相關(guān)數(shù)據(jù)[6]。表格式實(shí)測(cè)數(shù)據(jù)模型常用下式表達(dá)：

2.4 仿真前處理

2.4.1 材料模型

滾打輪的材料為 W2Mo9Cr4VCo8，假設(shè)其為完全的剛性體，不發(fā)生變形。花鍵坯軸的材料為 40Cr，通過(guò)拉伸試驗(yàn)、彈塑性理論計(jì)算以及冷滾打過(guò)程中的材料強(qiáng)化性質(zhì)，可得到常溫下 40Cr 應(yīng)力-應(yīng)變曲線，如圖 3 所示，40Cr 的性能參數(shù)如表 2 所列。

圖3 40Cr 應(yīng)力-應(yīng)變?cè)囼?yàn)曲線Fig.3 Stress-strain test curve of 40Cr

表2 40Cr 性能參數(shù)Tab.2 Performance parameters of 40Cr

2.4.2 網(wǎng)格劃分及約束

在三維金屬成形仿真時(shí)，常用的單元類型有六面體單元和四面體單元。為了提高計(jì)算效率，在工件的成型位置進(jìn)行局部網(wǎng)格細(xì)化。使用 DEFORM-3D 提供的四面體單元自動(dòng)劃分，花鍵坯軸的網(wǎng)格劃分如圖4 所示。

圖4 網(wǎng)格劃分Fig.4 Mesh division

冷滾打是在常溫條件下進(jìn)行，且加工過(guò)程中有冷卻液，所以在數(shù)值模擬中暫不考慮物體之間的熱傳遞，僅考慮的是物體之間的摩擦條件和接觸條件。

3 仿真結(jié)果分析

冷滾打花鍵成形過(guò)程主要發(fā)生在滾打輪與花鍵毛坯接觸后，即將開(kāi)始階段設(shè)置在滾打輪與花鍵毛坯接觸位置，按照設(shè)計(jì)工藝參數(shù)進(jìn)行數(shù)值模擬，得到了不同時(shí)刻的花鍵外形輪廓，如圖 5 所示。

圖5 花鍵冷滾打成形過(guò)程Fig.5 Process of spline cold rolling formation

從冷滾打塑性成形數(shù)值模擬結(jié)果來(lái)看，當(dāng)滾打輪高速軋入軸坯后，迫使金屬沿滾打輪外形向兩側(cè)自由流動(dòng)。隨著周期擊打深入，滾打輪底部的金屬逐步被擠壓到齒頂，這樣就形成了完整的花鍵齒形輪廓。通過(guò)測(cè)量花鍵齒形輪廓，發(fā)現(xiàn)兩側(cè)齒面并不完全對(duì)稱，主要是因?yàn)檫B續(xù)分度冷滾打與軸坯的單向連續(xù)轉(zhuǎn)動(dòng)形成了干涉，造成轉(zhuǎn)動(dòng)方向軸坯會(huì)更早與滾打輪接觸，其受到的擠壓力更大，所以該面的齒形更早完成。在實(shí)際加工中，可以通過(guò)調(diào)整冷滾打輪安裝偏角減輕運(yùn)動(dòng)干涉，來(lái)獲得花鍵左右齒面對(duì)稱誤差很小的齒廓。

4 試驗(yàn)驗(yàn)證

為了驗(yàn)證數(shù)值仿真的準(zhǔn)確性，在ZRMe9 型花鍵冷滾打機(jī)床上，按照數(shù)值仿真的工藝參數(shù)加工制造了1 件花鍵軸，冷滾打花鍵軸如圖 6 所示。利用 Serein-CMM FUNCTION1000 型三坐標(biāo)測(cè)量?jī)x，選擇花鍵上任一個(gè)完整齒形，取其垂直于軸線截面內(nèi)的 15 個(gè)均布點(diǎn)進(jìn)行測(cè)量。在后處理中，同樣選取數(shù)值模擬齒形上 15 個(gè)均布點(diǎn)的坐標(biāo)。利用 MATLAB 數(shù)據(jù)分析工具進(jìn)行擬合，分別得到了冷滾打加工和數(shù)值模擬齒廓曲線，如圖 7 所示。由圖 7 不難看出，數(shù)值仿真與加工齒廓曲線基本吻合。

圖6 冷滾打花鍵軸Fig.6 Cold rolling spline shaft

圖7 數(shù)值仿真與加工齒廓曲線Fig.7 Numerical simulation and machined tooth profile curve

5 結(jié)語(yǔ)

通過(guò)采用 DEFORM 有限元模擬技術(shù)，對(duì)花鍵冷滾打過(guò)程進(jìn)行了數(shù)值仿真，得到了花鍵齒廓塑性成形過(guò)程。通過(guò)分析對(duì)比數(shù)值模擬花鍵與實(shí)際加工花鍵，驗(yàn)證了數(shù)值模擬的正確性，為后續(xù)冷滾打花鍵成型機(jī)理的研究奠定了基礎(chǔ)。