固液兩相磨粒流研拋螺旋齒輪的數(shù)值模擬研究

衛(wèi)麗麗，李俊燁，李丹妮，周立賓

（長春理工大學(xué) 機(jī)電工程學(xué)院，長春 130022）

固液兩相磨粒流研拋螺旋齒輪的數(shù)值模擬研究

衛(wèi)麗麗，李俊燁，李丹妮，周立賓

（長春理工大學(xué) 機(jī)電工程學(xué)院，長春 130022）

為了探討固液兩相磨粒流研拋螺旋齒輪質(zhì)量的影響，以螺旋齒輪為對象，進(jìn)行數(shù)值模擬分析。通過分析可知：從軸向來看，在同一入口速度條件下，隨流體不斷流入，受曲面形狀的影響，流道截面積不斷變化，近壁面流體跟隨曲面形狀時(shí)刻做上下起伏、旋轉(zhuǎn)運(yùn)動(dòng)，運(yùn)動(dòng)方向和流動(dòng)狀態(tài)不斷變化，流體能量逐漸損失，拋光效果也隨之減弱；從徑向來看，流道越窄或流道曲面彎曲形狀變化越均勻，流道截面積越均勻，磨粒流研拋工件表面質(zhì)量均勻性越好。為促進(jìn)固液磨粒流超精密加工技術(shù)的不斷發(fā)展提供了技術(shù)支持。

兩相流；磨粒流研拋；螺旋齒輪；數(shù)值模擬

螺旋齒輪因其具有振動(dòng)小、噪聲低、節(jié)能高等優(yōu)點(diǎn)廣泛應(yīng)用于各種精密儀表等機(jī)構(gòu)傳動(dòng)中，其齒面的形狀精度與表面質(zhì)量對傳動(dòng)過程中的影響尤為明顯，需高精度的精密齒輪滿足上述要求，普通加工方法難以獲得高精度的螺旋齒輪，因此需用磨粒流拋光技術(shù)進(jìn)行后續(xù)拋光，從而達(dá)到其精度要求［1-5］。

磨粒流加工技術(shù)是以顆粒為磨削刀具，以流體為載體，通過流體在零部件表面的流動(dòng)來帶動(dòng)顆粒在工件表面產(chǎn)生滑移摩擦和碰撞，達(dá)到去除材料的目的。磨粒流拋光裝置示意圖如圖1所示，上下液壓缸在壓力推動(dòng)下使得活塞上下移動(dòng)，從而推動(dòng)磨料在工件表面產(chǎn)生相對運(yùn)動(dòng)，工件表面產(chǎn)生較大的剪切力，從而實(shí)現(xiàn)對加工表面材料的微量去除，達(dá)到拋光的目的［6-13］。

圖1 磨粒流拋光裝置示意圖

1 固液兩相磨粒流拋光螺旋齒輪流場數(shù)值模擬

1.1 三維模型建立與網(wǎng)格劃分

由于螺旋齒輪常用于一些精密傳動(dòng)中，所以對齒輪的加工要求較高，其齒面表面質(zhì)量直接影響著齒輪的傳動(dòng)精度，因此需要對齒輪齒面進(jìn)行精密加工。根據(jù)螺旋齒輪相關(guān)技術(shù)要求，在三維軟件solid?works中建立三維模型，螺旋齒輪模型具體參數(shù)為模數(shù)為1，齒數(shù)為16，壓力角20°，螺旋角45°，螺旋齒輪模型如圖2所示。

圖2 螺旋齒輪三維模型

模型選擇根據(jù)固液兩相磨粒流加工特性，并結(jié)合螺旋齒輪的工作性能，需要對螺旋齒輪的齒面進(jìn)行加工，在進(jìn)行數(shù)值計(jì)算時(shí)需對數(shù)學(xué)模型進(jìn)行選擇。在磨粒流在拋光葉輪過程中，流體要流經(jīng)葉輪的葉片表面，由于葉片的幾何形狀為彎曲的，流體的運(yùn)動(dòng)也是不穩(wěn)定的，主要選取的模型主要有混合相（Mixure）方程、能量（Energy）方程、湍流方程及離散相方程［14］。湍流方程選擇標(biāo)準(zhǔn)k-ε模型，近壁面處理方式選擇標(biāo)準(zhǔn)壁面方程，假定流體滿足連續(xù)條件，其中模型經(jīng)驗(yàn)系數(shù)C1ε=1.44，C2ε=1.92，Cmu= 0.09，k對應(yīng)的普朗特?cái)?shù)σk=1.0，ε對應(yīng)的普朗特?cái)?shù)σε=1.3。

1.2 計(jì)算模型邊界條件的設(shè)置

（1）入口邊界條件

固相設(shè)置：固體顆粒相為SiC顆粒，同樣采用速度入口條件，初始速度大小與液相相同，由于在磨粒流加工過程中，是以液壓油作為載體，SiC顆粒作為實(shí)質(zhì)的磨削刀具，所以SiC顆粒的濃度大小會影響磨削效果，若顆粒相濃度過低，則顆粒對壁面的碰撞幾率較小，加工效果不明顯，且加工效率較低，但是顆粒相濃度也不能過高，雖然在理論上顆粒濃度越大，與工件壁面碰撞幾率就越大，磨削效率也越高，但是由于流道內(nèi)部結(jié)構(gòu)復(fù)雜，當(dāng)濃度過大時(shí)，會使得流體運(yùn)動(dòng)受阻運(yùn)行不通暢，導(dǎo)致磨粒大量沉積，從而影響流體的湍流狀態(tài)，因此應(yīng)設(shè)置合理的顆粒相的體積分?jǐn)?shù)，最大限度提高加工效率，因此設(shè)置顆粒相體積分?jǐn)?shù)為0.2。

（2）出口邊界條件

由于在對流體仿真計(jì)算之前并不清楚出口處的壓力與速度，出口與外界直接相連，流動(dòng)狀態(tài)假定為完全發(fā)展的湍流，所以設(shè)定出口端邊界條件為自由出口。

（3）壁面邊界條件

由于在磨粒流加工時(shí)是對工件壁面進(jìn)行加工，在加工過程中由于工件是固定不動(dòng)的，對壁面的切削力主要來自于在外界壓力驅(qū)動(dòng)下流體所負(fù)載的磨粒流經(jīng)工件待加工表面與壁面產(chǎn)生的相對滑移碰撞，對于固定的工件待加工表面來說，磨粒與壁面的相對滑移即為粒子本身的運(yùn)動(dòng)，因此壁面邊界條件應(yīng)選擇無滑移的壁面邊界條件。

1.3 螺旋齒輪仿真結(jié)果分析

通過設(shè)置不同速度對螺旋齒輪近壁面的流場狀態(tài)進(jìn)行數(shù)值模擬，經(jīng)計(jì)算大約在125步時(shí)得到收斂，可得在對螺旋齒輪流場的仿真計(jì)算時(shí)模型選擇與參數(shù)設(shè)置正確，為更好分析磨粒流體對螺旋齒輪齒面的影響，通過后處理得到不同速度下螺旋齒輪近壁面的流場狀態(tài)，從而得出磨粒流對螺旋齒輪齒面的拋光效果。由于螺旋齒輪的傳動(dòng)主要靠兩輪齒之間的嚙合來傳遞力矩，所以需要進(jìn)行精密加工的主要部位是輪齒的齒面，因此對螺旋齒輪進(jìn)行分析時(shí)，主要對輪齒近壁面的流場進(jìn)行分析，通過對螺旋齒輪在不同入口速度條件下進(jìn)行流體仿真計(jì)算，得到的靜壓云圖如圖3所示。

圖3 不同入口速度下的靜壓云圖

從圖3不同入口速度下的靜壓云圖可以看出，在同一入口速度下，輪齒上端所受靜壓較大，這是因?yàn)樵谘剌S線方向，流體前行時(shí)會受到齒輪彎曲輪齒影響，由于輪齒形狀呈現(xiàn)一定的彎曲傾斜性，當(dāng)磨料流體沿軸向運(yùn)動(dòng)時(shí)，會受到輪齒螺旋齒面的阻擋，外界壓力會首先作用于輪齒上端，因此輪齒上端齒面的所承受的靜壓也較大；當(dāng)磨料流體流入兩螺旋齒之間后，隨兩輪齒間不斷深入，由于彎曲齒面的阻擋，會造成一定的壓力損失，導(dǎo)致輪齒齒面近壁面的靜壓自上而下逐漸較?。浑S入口速度不斷增加，螺旋輪齒近壁面的靜壓也在不斷增加，意味著流動(dòng)的流體對壁面產(chǎn)生的有效壓力越大。未受擾動(dòng)的流體靜壓與流體完全受阻時(shí)的總壓之差即為動(dòng)壓，與之相對應(yīng)的不同入口速度條件下的動(dòng)壓如圖4所示。

圖4 不同入口速度下的動(dòng)壓云圖

通過在不同速度條件下對螺旋齒輪不同部位進(jìn)行選點(diǎn)取值，得到了不同速度條件下不同部位的動(dòng)壓力分布，如表1所示。

如圖4所示的動(dòng)壓云圖和表1動(dòng)壓分布表可以看出在螺旋輪齒上端湍動(dòng)能大于下端，且在齒頂部位湍動(dòng)能更加明顯，在輪齒的兩側(cè)面湍動(dòng)能明顯大于齒根部，這是因?yàn)樵趶较蛏淆X頂端與外界約束裝置內(nèi)壁面距離最小，所形成的流道面積最??；同時(shí)由于輪齒為螺旋形狀，流體流經(jīng)該區(qū)域時(shí)受到流道形狀的改變，導(dǎo)致阻力變大，動(dòng)能轉(zhuǎn)化為壓力能，因此齒頂處的壓力能會大于齒根處，而且當(dāng)初始碰到齒輪輪齒時(shí)，由于流體的運(yùn)動(dòng)方向會由平行于軸向方向轉(zhuǎn)為與軸線成一定角度的螺旋流向，運(yùn)動(dòng)方向的突然改變且受到流道壁面阻擋，會導(dǎo)致此時(shí)流體所受的阻力猛然增加，壓力能瞬時(shí)增大，因此會在上端輪齒齒頂處形成的動(dòng)壓最大，對輪齒上端齒頂處拋光效果最好。

表1 不同入口速度條件下不同部位的動(dòng)壓分布表

圖5 不同入口速度下的湍動(dòng)能云圖

由圖5不同入口速度下的湍動(dòng)能云圖可以看出在齒輪上端磨料流體剛進(jìn)入兩齒之間的流道時(shí)，湍動(dòng)能最大，隨流體在兩螺旋輪齒之間不斷前行，湍動(dòng)能明顯減弱；同樣是因?yàn)榱鞯佬螤畹脑蚴艿阶枇?dǎo)致能量降低，但在各輪齒頂端湍動(dòng)能基本一致，這是由于與外界約束裝置內(nèi)壁面形成的流道是一樣的，所以在同一橫截面上不同輪齒同一部位湍動(dòng)能也基本一致，隨入口速度的不斷增加，湍動(dòng)能在不斷增加，所以可以通過增加入口速度來提升輪齒齒面的表面質(zhì)量。

圖6 不同入口速度下的湍流強(qiáng)度云圖

通過在不同入口速度條件下的湍流云圖進(jìn)行匯總，得出了螺旋齒輪不同部位的湍流強(qiáng)度，具體分布如表2所示。

表2 不同入口速度條件下不同部位的湍流強(qiáng)度分布表

通過圖6所示的湍流強(qiáng)度云圖和表2所示的動(dòng)壓分布表可以看出，隨入口速度不斷增強(qiáng)，湍流強(qiáng)度不斷增強(qiáng)，湍流強(qiáng)度增加的梯度卻有逐漸削弱的趨勢，即當(dāng)湍流強(qiáng)度增大到一定值后，入口速度的再次增加，湍流強(qiáng)度有趨于一個(gè)恒定值而不再增加的趨勢；輪齒頂端的湍流強(qiáng)度大于輪齒兩側(cè)及齒根部位，說明流道截面積的縮小有利于增強(qiáng)湍流強(qiáng)度，流體的無序性運(yùn)動(dòng)使得其附載的磨粒對壁面進(jìn)行無規(guī)則的碰撞或在壁面上進(jìn)行摩擦滑移，從而使得表面紋理更加光整，輪齒表面質(zhì)量更高。

圖7 不同入口速度下的壁面剪切力云圖

通過在不同速度條件下的壁面剪切力云圖進(jìn)行匯總，螺旋齒輪不同部位所受的壁面剪切力如表3所示。

表3 不同入口速度條件下不同部位的壁面剪切力分布表

通過圖7所示的壁面剪切力云圖和表3所示的壁面剪切力分布表可以看出輪齒上端壁面剪切力較大，當(dāng)磨料流體初始接觸螺旋齒輪時(shí)，能量較為充足，首先作用于上端輪齒，所以輪齒齒面的上部分加工效果更為明顯，但在兩輪齒之間的齒槽處，壁面剪切力卻很小，這是由于流道過大造成的，當(dāng)流體快速通過該流道時(shí)，會受到流道形狀影響，導(dǎo)致流體流經(jīng)齒槽處時(shí)壓力較小，流體對壁面的貼近力不明顯，從而使粒子與齒槽處壁面碰撞的機(jī)會大大減少，齒槽處加工效果弱于齒頂處，且隨入口速度不斷增大壁面剪切力也在逐步增大，壁面剪切力與速度成正相關(guān)，所以可以通過增加入口速度的方式來提高加工磨粒流體附帶的磨料粒子對壁面的碰撞機(jī)會，從而提高磨粒流拋光質(zhì)量。

2 結(jié)論

（1）固液兩相磨粒流研拋的螺旋齒輪是曲面形狀，需要待加工壁面與外部裝置相互配合組成封閉的型腔流道，故在進(jìn)行固液兩相磨粒流數(shù)值分析時(shí)需對待加工曲面外部進(jìn)行流體包覆。

（2）從磨粒流研拋螺旋齒輪的數(shù)值云圖可知，隨著入口速度的不斷增加，螺旋齒輪的靜壓、動(dòng)壓、湍動(dòng)能、湍流強(qiáng)度、壁面剪切力都隨之增加，因而可以通過增加入口速度來提高磨粒流研拋異形曲面的表面質(zhì)量。在同一入口速度條件下，從軸向方向來看，當(dāng)磨料流體初始進(jìn)入流道時(shí)，相關(guān)各參數(shù)都相對較大，說明加工效果也最好，隨流體不斷前行，受曲面形狀和流道截面積不斷變化的影響，近壁面流體隨著曲面結(jié)構(gòu)形狀時(shí)刻做上下起伏和旋轉(zhuǎn)運(yùn)動(dòng)，運(yùn)動(dòng)方向和流動(dòng)狀態(tài)不斷變化，流體能量逐漸損失，拋光效果也隨之減弱，在出口處磨粒流研拋?zhàn)畈?；從徑向來看，流道越窄或流道曲面彎曲形狀變化越均勻，各?xiàng)參數(shù)的變化幅度也較小，流道截面積越均勻，磨粒流研拋異形曲面獲得的表面質(zhì)量均勻性越高。

［1］ 李俊燁，胡敬磊，張心明，等.固液兩相流拋光內(nèi)齒輪齒面夾具的設(shè)計(jì)和分析［J］.制造業(yè)自動(dòng)化，2016（11）：66-69+74.

［2］ 唐進(jìn)元，劉艷平.直齒面齒輪加載嚙合有限元仿真分析［J］.機(jī)械工程學(xué)報(bào)，2012，48（5）：124-131.

［3］ 姚鵬華.基于ANSYS的汽車變速器斜齒輪的有限元分析［J］.機(jī)械研究與應(yīng)用，2016，29（6）：18-20.

［4］ 尚春民，范景峰，張心明.兩種齒輪泵脈動(dòng)特性對比分析［J］.長春理工大學(xué)學(xué)報(bào)：自然科學(xué)版，2004，27（3）：41-43.

［5］ Kenda J，Duhovnik J，Tav?ar J，et al.Abrasive flow ma?chining applied to plastic gear matrix polishing［J］.Inter?national Journal of Advanced Manufacturing Technolo?gy，2014，71（71）：141-151.

［6］ 李敏，袁巨龍，吳喆，等.復(fù)雜曲面零件超精密加工方法的研究進(jìn)展［J］.機(jī)械工程學(xué)報(bào)，2015（5）：178-191.

［7］ 李琛，計(jì)時(shí)鳴，譚大鵬，等.軟性磨粒流加工特性及近壁區(qū)域微切削機(jī)理［J］.機(jī)械工程學(xué)報(bào)，2014，50（9）：161-168.

［8］ 丁金福，劉潤之，張克華，等.磨粒流精密光整加工的微切削機(jī)理［J］.光學(xué)精密工程，2014，22（12）：3324-3331.

［9］ 計(jì)時(shí)鳴，李琛，譚大鵬，等.基于Preston方程的軟性磨粒流加工特性［J］.機(jī)械工程學(xué)報(bào)，2011，47（17）：156-163.

［10］ 張雷，吳紹菊，李俊燁，等.基于多物理耦合場的固液兩相磨粒流伺服閥閥芯噴嘴拋光研究［J］.長春理工大學(xué)學(xué)報(bào)：自然科學(xué)版，2016，39（4）：87-89.

［11］ 喬澤民，李俊燁，張雷，等.磨粒流研磨液顆粒特性的耗散粒子動(dòng)力學(xué)數(shù)值模擬［J］.長春理工大學(xué)學(xué)報(bào)：自然科學(xué)版，2016，39（5）：61-64.

［12］ Walia R S，Shan H S，Kumar P.Morphology and in?tegrity of surfaces finished by centrifugal force assisted abrasive flow machining［J］.International Journal of Advanced Manufacturing Technology，2008，39（11）：1171-1179.

［13］ Mali H S，Manna A.Optimum selection of abrasive flow machining conditions during fine finishing of Al/ 15wt%SiC-MMC using Taguchi method［J］.Interna?tional Journal of Advanced Manufacturing Technology，2010，50（9-12）：1013-1024.

［14］ 宋林森，史國權(quán)，李占國.ANSYS在激光打孔溫度場仿真中的應(yīng)用［J］.長春理工大學(xué)學(xué)報(bào)：自然科學(xué)版，2006，29（4）：19-21.

Study on Numerical Simulation of Helical Gear with Solid Liquid Two Phase Abrasive Flow

WEI Lili，LI Junye，LI Danni，ZHOU Libin
（School of Mechatronical Engineering，Changchun University of Science and Technology，Changchun 130022）

In order to investigate the effect of solid-liquid two-phase abrasive flow on the quality of helical gear，the helical gear is taken as the object，the numerical simulation analysis is carried out.Through analysis we can see：from the horizontal point of view，at the same inlet velocity conditions，with the fluid flowing in，influenced by the surface shape of the channel area is con?stantly changing，the near wall fluid surface shape follow ups and downs and rotary motion，changing the direction of movement and the flow state of fluid energy loss，the polishing effect also weakened；from the radial point of view，the narrower the flow path or the more uniform the curved shape of the runner surface，the more uniform the cross-sectional area of the flow channel，the higher the uniformity of the surface quality of the workpiece.Which can provide theoretical basis for the continuous improve?ment of the ultra-precision machining technology.

two-phase；abrasive flow polishing；helical gear；numerical simulation

TH117.1

A

1672-9870（2017）03-00054-05

2017-03-17

國家自然科學(xué)基金資助項(xiàng)目（51206011）；吉林省科技發(fā)展計(jì)劃資助項(xiàng)目（20160101270JC，20170204064GX）；吉林省教育廳項(xiàng)目（吉教科合字［2016］第386號）

衛(wèi)麗麗（1992-），女，碩士研究生，E-mail：916033207@qq.com

李丹妮（1982-）女，碩士，副研究員，E-mail：ldn@cust.edu.cn