砂帶磨損對(duì)螺桿轉(zhuǎn)子磨削去除深度影響分析

摘要 為保證螺桿轉(zhuǎn)子表面磨削加工材料去除的均勻性，考慮砂帶磨料磨損對(duì)磨削去除深度的影響，建立螺桿轉(zhuǎn)子材料去除深度數(shù)學(xué)模型。根據(jù)彈塑性變形理論，結(jié)合磨削過(guò)程中磨粒的磨損規(guī)律，建立單顆磨粒材料去除率理論模型。將磨粒的磨損分為快速磨損階段和穩(wěn)定磨損階段2個(gè)階段?；趩挝唤佑|面積上磨粒的數(shù)目和磨粒的出刃規(guī)律，運(yùn)用數(shù)值積分和接觸壓力計(jì)算出磨粒的最大切削深度，利用積分建立基于砂帶磨損的宏觀材料去除率數(shù)學(xué)模型。該數(shù)學(xué)模型考慮了砂帶的磨損，其預(yù)測(cè)值和試驗(yàn)值相比，最大誤差為9.6%，最小誤差為4.1%。該試驗(yàn)數(shù)據(jù)充分驗(yàn)證了模型的有效性，可以為保障螺桿轉(zhuǎn)子廓形精度提供理論基礎(chǔ)。

關(guān)鍵詞 砂帶磨損；數(shù)學(xué)模型；去除深度；螺桿轉(zhuǎn)子；模型驗(yàn)證

中圖分類(lèi)號(hào) TG70；TG580 文獻(xiàn)標(biāo)志碼 A

文章編號(hào) 1006-852X（2024）02-0228-09

DOI 碼 10.13394/j.cnki.jgszz.2023.0094

收稿日期 2023-04-23 修回日期 2023-06-30

砂帶磨削幾乎可以加工工程上的所有材料，有著“萬(wàn)能磨削”之稱(chēng)[1]，但是在磨削過(guò)程中工件表面材料磨削去除深度存在差異，很大程度上會(huì)影響工件表面的形面精度。如何精準(zhǔn)地控制加工參數(shù)使工件在磨削過(guò)程中獲得均勻的材料去除率變得十分重要。預(yù)測(cè)材料去除率的模型較多，但總的來(lái)說(shuō)可分為基于加工理論、試驗(yàn)研究、設(shè)計(jì)實(shí)驗(yàn)和人工智能4類(lèi)[2]。

ZHE 等[3]根據(jù) Preston 方程及 Hertz 接觸理論預(yù)測(cè)磨削過(guò)程中工件材料的去除深度，以磨削壓力為重點(diǎn)研究對(duì)象，研究了砂帶磨削過(guò)程中接觸壓力對(duì)材料去除的影響。遲玉倫等[4]針對(duì)傳統(tǒng)材料去除理論模型未考慮工件表面不規(guī)則變形的問(wèn)題開(kāi)展了研究，利用分段函數(shù)建立了材料去除率理論模型。劉斐等[5]針對(duì)接觸輪在張緊力的作用下發(fā)生變形的問(wèn)題，采用復(fù)變函數(shù)解析法建立材料去除率模型。齊俊德等[6]研究了微觀單顆磨粒對(duì)材料去除率的影響規(guī)律，結(jié)合 Hertz接觸理論以及概率統(tǒng)計(jì)理論從單顆磨粒的材料去除率模型推廣到整體砂帶對(duì)工件的材料去除率模型，但此模型并未考慮砂帶磨損因素。 WANG 等[7]通過(guò) Hertz 接觸理論分析作用在工件表面的磨削力，采用 Hammann模型建立材料去除深度模型，結(jié)果表明非線性材料去除率模型的精度高于線性材料去除率模型的。

目前有很多學(xué)者利用仿真軟件將磨粒簡(jiǎn)化成球體或錐體等，對(duì)工件進(jìn)行切削仿真，觀察磨削參數(shù)對(duì)材料去除率的影響[8]。王榮全[9]對(duì)鋼軌打磨的磨削過(guò)程進(jìn)行了建模和實(shí)驗(yàn)研究，分析了砂帶磨粒對(duì)鋼軌的磨削作用，并以砂帶磨損帶來(lái)的影響作為下一步的研究目標(biāo)。杜瑤等[10]將砂帶的磨削過(guò)程分為3個(gè)階段，分別探究了不同磨削階段下砂帶磨削中密度纖維板的表面粗糙度以及材料去除率，此研究針對(duì)的是磨削平面，未對(duì)磨削曲面進(jìn)行研究。KARKALOS 等[11]建立了2個(gè) 金剛石磨粒磨削單晶銅模型，通過(guò)研究不同工藝參數(shù)下磨屑的形貌，分析得出最佳磨削參數(shù)，此研究也未考慮砂帶磨損帶來(lái)的影響。楊雅琦等[12]在磨削木材過(guò)程中結(jié)合砂帶磨損機(jī)理對(duì)磨削后的木材表面粗糙度進(jìn)行分析，同時(shí)也探究了工藝參數(shù)對(duì)砂帶磨損的影響規(guī)律。

為實(shí)現(xiàn)磨削時(shí)工件表面材料去除均勻性，相關(guān)學(xué)者采用各種方式進(jìn)行研究，例如，通過(guò)分析工件與磨具之間的接觸應(yīng)力或者通過(guò)實(shí)驗(yàn)仿真的方式進(jìn)行理論研究，將實(shí)驗(yàn)與智能算法相結(jié)合以建立材料去除率預(yù)測(cè)模型，但是這些研究通常局限在平面磨削與微觀的點(diǎn)接觸磨削上，由于螺旋曲面的測(cè)量以及壓力分析較為復(fù)雜，因此對(duì)復(fù)雜曲面的研究很少。當(dāng)磨削螺桿轉(zhuǎn)子類(lèi)零件時(shí)，磨削裝置與螺桿轉(zhuǎn)子之間的壓力隨著磨削路線的變化而發(fā)生變化，從而難以掌控磨削的深度。因此研究復(fù)雜磨削路線中磨料磨損對(duì)材料去除深度的影響顯得尤為重要，本文根據(jù)單磨粒的切削過(guò)程，分析不同切削深度下單磨粒的材料去除機(jī)理（如圖1所 示），并由此建立基于砂帶磨損的材料去除率數(shù)學(xué)模型。

1 砂帶磨粒磨損規(guī)律

1.1 磨粒模型簡(jiǎn)化及測(cè)量

工藝參數(shù)、砂帶材料以及工件材料共同決定了砂帶壽命，砂帶選型主要是挑選適合磨削多頭螺桿轉(zhuǎn)子的砂帶，即尋求一種磨削能力較強(qiáng)、服役壽命長(zhǎng)、磨削后能夠滿足螺桿轉(zhuǎn)子表面粗糙度要求的砂帶。目前棕剛玉、陶瓷、碳化硅、鋯剛玉4種類(lèi)型的砂帶在磨削作業(yè)時(shí)被廣泛使用，其具體的磨削屬性見(jiàn)表1。磨粒的幾何形狀會(huì)影響砂帶的磨削性能。

根據(jù)表1可以得出鋯剛玉砂帶的磨削效率較高且 耐磨性能最好。經(jīng)過(guò)實(shí)驗(yàn)發(fā)現(xiàn)，粒度為178μm 的鋯剛 玉砂帶可以滿足螺桿轉(zhuǎn)子的加工需求，通過(guò)超景深顯 微鏡觀察發(fā)現(xiàn)該類(lèi)型砂帶磨粒成圓頂圓錐狀，由圓頂 半徑 R、錐頂角α和攻角β表示，磨粒簡(jiǎn)圖如圖2所示，磨粒磨損高度示意圖如圖3所示。

式中：h'為磨粒圓頂?shù)饺清F體的高度，即虛擬值； R 為磨粒頂端半徑；Δh 為砂帶的磨損高度；H 為磨粒的 出刃高度。

采用超景深顯微鏡對(duì)砂帶磨損前后的高度進(jìn)行測(cè)量，超景深顯微鏡的測(cè)量界面如圖4所示。

砂帶基體為長(zhǎng)方體，其周長(zhǎng) l 為1400 mm，寬 b 為10 mm，假設(shè)底面圓形的直徑為 dm，磨粒的平均粒徑可由式（3）計(jì)算得出：

dm =68· Mv ?1.4"""""""""""""""""""""" （3）

式中：Mv 為磨粒的平均粒徑，mm。

建立砂帶表面形貌需要對(duì)其表面坐標(biāo)進(jìn)行采集，其中采樣間隔會(huì)影響砂帶表面形貌，測(cè)量間距越小越 接近真實(shí)的砂帶表面形貌，但這會(huì)導(dǎo)致測(cè)量數(shù)據(jù)點(diǎn)的 增加，進(jìn)而增加工作量，BLUNT 等[13]提出了適合砂帶 測(cè)量的采樣間隔，如式（4）所示。根據(jù)本文選擇的砂帶，采樣間隔應(yīng)為40μm。砂帶的植砂密度通常保持在砂帶表面積的40%～60%，為了重建砂帶表面圖像，采用超景深顯微鏡獲取砂帶表面形貌坐標(biāo)點(diǎn)，利用編程軟件對(duì)坐標(biāo)點(diǎn)進(jìn)行三維重建，砂帶磨損前后的表面形貌如圖5、圖6所示。

1.2 磨粒磨損規(guī)律

根據(jù)圖5、圖6可以得出，砂帶磨損前后磨粒出刃 高度不一致，測(cè)量發(fā)現(xiàn)未磨損的砂帶磨粒出刃分布規(guī) 律近似于高斯分布，磨損后的砂帶磨粒的磨損高度趨 于一致。將砂帶磨削25 min，每隔1 min 測(cè)量砂帶磨粒 的出刃高度和螺桿轉(zhuǎn)子表面粗糙度，測(cè)量結(jié)果如圖7、圖8所示。

根據(jù)圖7、圖8可以得出，砂帶磨粒的出刃高度在 磨削前5 min 時(shí)磨損較快，且砂帶磨粒出刃高度不一致，螺桿轉(zhuǎn)子表面粗糙度為0.7～1.0μm，此時(shí)稱(chēng)為快速磨 損階段。在磨削5～20 min 內(nèi)砂帶的磨損較為穩(wěn)定，砂帶磨粒出刃高度較為一致，螺桿轉(zhuǎn)子表面粗糙度為0.5 左右，此時(shí)稱(chēng)為穩(wěn)定磨損階段。在磨削20 min 后采用 粗糙度測(cè)量?jī)x測(cè)出的工件磨削表面的粗糙度不滿足驗(yàn) 收標(biāo)準(zhǔn)，故在確定磨削螺桿轉(zhuǎn)子表面材料去除率時(shí)，僅需計(jì)算磨削20 min 內(nèi)的材料去除率。

2 壓力仿真及磨粒切削仿真

2.1 接觸壓力仿真

在計(jì)算磨削材料去除率時(shí)，需確定砂帶與螺桿轉(zhuǎn) 子之間的接觸壓力。由于螺桿轉(zhuǎn)子型面復(fù)雜，解析方 式不易獲得接觸區(qū)域壓力分布情況，故采用有限元仿 真分析的方式獲取螺桿轉(zhuǎn)子與接觸輪之間的接觸壓力分布，具體如圖9所示，工件與接觸輪的材料屬性見(jiàn)表2。

在進(jìn)行網(wǎng)格劃分時(shí)，將接觸處的網(wǎng)格尺寸設(shè)置為0.1 mm ×0.1 mm，在接觸輪磨削時(shí)將正壓力設(shè)置為80 N。由以上條件，可以得到接觸區(qū)域各節(jié)點(diǎn)處的接觸應(yīng)力值。根據(jù)圖9接觸輪與螺桿轉(zhuǎn)子之間的壓力仿真結(jié)果可知，該參數(shù)下接觸輪的最大壓力約為0.866 MPa。

2.2 磨粒切削仿真

砂帶磨削過(guò)程中，磨粒逐漸發(fā)生磨損，其頂端半徑會(huì)逐漸增大，經(jīng)過(guò)測(cè)量發(fā)現(xiàn)初始時(shí)磨粒的圓頂半徑約 為10μm，在磨削5 min 后磨粒的圓頂半徑約為30μm，當(dāng)磨削25 min 時(shí)，磨粒的頂端半徑約為55μm，因此將 磨粒簡(jiǎn)化為圓頂圓錐進(jìn)行磨削仿真，根據(jù)其磨削力的 變化，計(jì)算磨損系數(shù)，磨削仿真如圖10至圖12所示。磨削力變化規(guī)律如圖13所示，參數(shù)設(shè)置如下：磨粒彈 性模量為300 GPa，泊松比為0.24，切削速度為10 m/s，工件采用六面體進(jìn)行網(wǎng)格劃分。

從圖10至圖12可以得出，磨粒的圓頂半徑越大工件與磨粒之間的應(yīng)力越小，因?yàn)槟チA頂半徑越大磨削時(shí)磨粒切削深度降低，磨粒與工件的接觸面積減小且應(yīng)力降低。由圖13可知，磨粒的磨削力隨著頂端半徑的增大而增大，因此砂帶磨鈍后切削能力下降。根據(jù)圖13將砂帶的磨削分為2部分，并根據(jù)磨削力的變化計(jì)算砂帶磨粒的磨損系數(shù)ξ， 即當(dāng)前磨損形態(tài)磨粒的磨削力與初始時(shí)未磨損磨粒磨削力的比值：

式中：Fc 為當(dāng)前磨損形態(tài)磨粒的磨削力；Fq 為未磨損磨粒的磨削力。

3 材料去除數(shù)學(xué)模型建立

3.1 磨粒切削時(shí)的彈塑性變形

當(dāng)最高磨粒的出刃高度小于螺桿轉(zhuǎn)子的最大彈性變形量時(shí)，砂帶表面的磨粒與螺桿轉(zhuǎn)子之間只發(fā)生彈性變形。當(dāng)單顆磨粒受到的壓力為 F0時(shí)，螺桿轉(zhuǎn)子發(fā)生彈性變形量，如式（6）～（9）所示[14]，單顆磨粒與工件作用示意圖如圖14所示。

式中：Pv 為單顆磨粒受到的平均壓強(qiáng)；Rm 為磨粒底端半徑；F0為施加壓力；E 為綜合彈性模量。

?=π（1） 16（9）RmE（F02）2）3（1）

式中：δ為彈性變形切入深度。

式中：δmax 為工件在發(fā)生彈性變形時(shí)，磨粒的最大切入深度。

式中：E1、E2分別為工件和磨粒的彈性模量；v1、v2分別為工件和磨粒的泊松比。經(jīng)過(guò)計(jì)算發(fā)現(xiàn)工件發(fā)生彈性變形時(shí)磨粒的切入深度為7.63×10-5μm，遠(yuǎn)小于磨粒塑性變形時(shí)磨粒的切入深度，所以忽略工件的彈性變形。

3.2 最高磨粒刺入深度求解

當(dāng)砂帶磨粒所受的切向力足夠大時(shí)，磨粒會(huì)對(duì)螺桿轉(zhuǎn)子表面產(chǎn)生耕犁或切削作用，此時(shí)磨粒受到的壓強(qiáng)通常取 Pv=HB/2。由于磨粒的形狀為圓頂圓錐，其切入螺桿轉(zhuǎn)子表面的部分分為球體和錐體2部分，球體的切入表面積 S1如式（11）所示，錐體的切入表面積 S2如式（13）所示。

式中：S1為圓頂部分與螺桿轉(zhuǎn)子的接觸面積。

式中：Fs1為磨粒圓頂對(duì)螺桿轉(zhuǎn)子產(chǎn)生的塑性變形力。

式中：S2為磨粒與螺桿轉(zhuǎn)子的接觸面積。

式中：Fs2為磨粒對(duì)螺桿轉(zhuǎn)子產(chǎn)生的塑性變形力。

采用超景深顯微鏡對(duì)砂帶磨粒單位面積的數(shù)目和出刃高度進(jìn)行測(cè)量，砂帶磨粒與螺桿轉(zhuǎn)子相互作用的情況如圖15所示，以磨粒平均出刃高度作為原點(diǎn)建立直角坐標(biāo)系，其中ω’為磨粒最高出刃高度。

可通過(guò)柱狀圖描繪出砂帶表面磨粒出刃高度的分布密度曲線，從而對(duì)其整體分布規(guī)律進(jìn)行評(píng)估，通過(guò)測(cè)量，砂帶磨粒出刃高度分布如圖16所示。

通過(guò)圖16可以計(jì)算出磨粒的平均高度為177.5μm，標(biāo)準(zhǔn)差為37.8。

式中：Nz 為單位面積上的磨粒總數(shù)占比率；Nm 為測(cè)量區(qū)域內(nèi)磨?？倲?shù)；M1、N1分別為測(cè)量的行數(shù)和列數(shù)；ΔL 為測(cè)量的步長(zhǎng)。

f（h）= e?（h?u）2，|h|lt;3σ""""" （16）

式中：σ為出刃高度均方差。

通過(guò)分析得出，砂帶單位面積上磨粒對(duì)螺桿轉(zhuǎn)子產(chǎn)生的壓力 Pz 主要由磨粒的彈性變形力與塑性變形力2部分組成，設(shè)最高磨粒的切入深度為ω′，可得出：

h0=3σ?ω′"""""""""""""""""" （17）

當(dāng)最高磨粒的切入深度ω′小于彈性變形的壓入深度，即3σ-ω′≤h0≤3σ-ω′+δmax 時(shí)，此時(shí)單位面積上磨粒對(duì)螺桿轉(zhuǎn)子產(chǎn)生的彈性變形力為：

Ft = lh 0+?max Nz F0f（h）dh"""""""" （18）

當(dāng)工件發(fā)生塑性變形時(shí)，假設(shè)最高磨粒的切入深度為δmaxlt;ω′lt;R，砂帶單位面積上磨粒對(duì)螺桿轉(zhuǎn)子產(chǎn)生的塑性變形力為：

F1= lh NzFS1f（h）dh"""""""" （19）

假設(shè)最高磨粒的切入深度為 Rlt;ω'lt;3σ， 砂帶單位面積上磨粒對(duì)螺桿轉(zhuǎn)子產(chǎn)生的塑性變形力為：

F2= lh0（3） R NzFS2f（h）dh

磨粒單位面積所受的外力為：

結(jié)合式（18）～式（20），通過(guò)仿真求出單位面積上的壓強(qiáng)分布，進(jìn)而可以得到最高磨粒的實(shí)際切削深度。

3.3 材料去除模型建立

在磨削過(guò)程中，假設(shè)在單位時(shí)間 t 內(nèi)磨粒與螺桿轉(zhuǎn)子接觸掃掠過(guò)的體積為微元 Q。砂帶磨粒進(jìn)入磨削區(qū)域后對(duì)微元持續(xù)地進(jìn)行材料去除，直至微元處的體積被磨粒切削掉為止。微元 Q 的總體積為接觸區(qū)域內(nèi) 所有磨粒切削材料體積的總和，將微元 Q 進(jìn)行放大，dx、 dy、dz 分別為微元 Q 的長(zhǎng)寬高，微元簡(jiǎn)圖如圖17所示。

由圖17可得，微元的表面積為：

S 3= dldx""""""""""""""""""""""""" （22）

式中：dl 為磨粒單位時(shí)間內(nèi)砂帶旋轉(zhuǎn)磨粒經(jīng)過(guò)微元的實(shí)際長(zhǎng)度。

dl = VS dt"""""""""""""" （23）

V′= Vc ± wr"""nbsp;"""""""""""""""""" （24）

式中：VS 為砂帶的選擇速度；V'為砂帶與螺桿轉(zhuǎn)子的相對(duì)速度；wr 為工件旋轉(zhuǎn)速度，同向?yàn)樨?fù)，逆向?yàn)檎?/p>

單位時(shí)間內(nèi)微元處參與磨削的磨粒數(shù)為：

N2= S 1Nz lh0（3）σ f（h）dh""""""""" （25）

當(dāng)磨粒的切削深度小于球冠半徑時(shí)，磨粒頂部球冠切削橫截面積為：

ΔA1=πR2 R2 sinαcosβ"""""" （26）

單顆磨粒的球冠切削體積為：

V1=ΔA1 dy"""""""""""""" （27）

當(dāng)磨粒的切削深度大于球冠半徑時(shí)，磨粒切入的橫截面積為：

ΔA2=πR2 R2 sin αcos β+

（h? h0+ h′）2tan α? Rsin α（h′+ R ? Rsin α）（28）

磨粒的去除體積為：

V2=ΔA2 dy""""nbsp;"""""""" （29）

微元處磨?？傮w材料去除體積為：

V =dxdydz = h 0+R e?h2/2σ2ΔA1 dh+

lh0（3） R e?h2/2σ2ΔA2 dh）dy

在微元處磨粒的磨削深度為：

H1= d（d）y（z）= d（N）x（2）dy2（V0）= vf（ξ）2π（Nz）σ lh 0+R e?h2/2σ2（ΔA1）dh+lh0（3） R e?h2/2σ2ΔA2 dh）dh（31）

將式（31）簡(jiǎn)化為：

H1= lh 0+R e?h2/2σ2（ΔA1）dh+lh0（3） R e?h2/2σ2ΔA2 dh）dh

對(duì)于砂帶磨削穩(wěn)定期，其磨粒刺入深度區(qū)域一致，因此將磨粒形貌化作一致，計(jì)算出單顆磨粒的材料去 除量，結(jié)合接觸面積上的總磨粒數(shù)即可計(jì)算總材料去 除量，假設(shè)磨粒從某一位置切削到另一位置走過(guò)的距 離為 L1，則磨粒切削掉的材料總體積應(yīng)為：

V0=¨Δ（Δ）A2（A1）·（·）L1（L1），，（） ap（ap） R（R）"""""""" （33）

式中：ap 為穩(wěn)定磨削期間磨粒的切削深度。

當(dāng)磨粒與螺桿轉(zhuǎn)子接觸的時(shí)間無(wú)限小時(shí)，單顆磨粒切削的材料量為：

dap = v′dt"""""""""""""""""""""" （34）

若已知砂帶磨損的起始位置，則可計(jì)算磨削穩(wěn)定期的磨削深度，具體如下：

ap = l 2 Ldl""""""""""" （35）

穩(wěn)定磨削期間總的材料磨削量為：

Nz · SξPvΔA1 ll2 Ldl， ap lt; R

Q總=〈 N · Sξ（H）P（B）ΔA2" l2（l1）""""""""""""" （36）

綜上所述，砂帶磨損對(duì)磨削螺桿轉(zhuǎn)子材料去除率的求解流程如圖18所示。

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

4.1 試驗(yàn)設(shè)備及流程

本文采用自主研發(fā)的磨削裝置對(duì)螺桿轉(zhuǎn)子凹面進(jìn)行試驗(yàn)，首先設(shè)置試驗(yàn)參數(shù)，然后將磨削裝置的接觸輪調(diào)整到螺桿轉(zhuǎn)子凹面，使其沿著螺旋線進(jìn)行運(yùn)動(dòng)，最后采用線性激光測(cè)量?jī)x對(duì)螺桿轉(zhuǎn)子磨削前后的表面輪廓進(jìn)行測(cè)量，一共測(cè)量8組試驗(yàn)。

4.2 螺桿轉(zhuǎn)子磨削材料去除率模型試驗(yàn)驗(yàn)證

使用 LJ-X Navigator 軟件采集 LJ-X 系列測(cè)得的圖 像及坐標(biāo)點(diǎn)，并將采集的結(jié)果通過(guò) LJ-X Navigator 作為 數(shù)字日記進(jìn)行輸出，打磨之前將線性激光位移傳感器 安裝在夾具上，使線性激光垂直打在螺桿轉(zhuǎn)子截面輪 廓上，調(diào)整夾具至合適的位置，使線性激光測(cè)量?jī)x測(cè)量 出的坐標(biāo)位置在有效范圍內(nèi)。測(cè)出的部分?jǐn)?shù)據(jù)點(diǎn)見(jiàn)表3， X 坐標(biāo)代表螺桿輪廓數(shù)據(jù)點(diǎn)的橫向位置，Y 坐標(biāo)代表螺 桿輪廓數(shù)據(jù)點(diǎn)的縱向位置，磨削前后需要將坐標(biāo)點(diǎn)相 匹配，匹配原則如式（37）所示，實(shí)際螺桿轉(zhuǎn)子廓形與 理論螺桿轉(zhuǎn)子部分廓形如圖19所示。

mind = è （（xi（′）? xi ）2+（yi（′）? yi ）2）

式中：mind 為磨削前后坐標(biāo)位置的最小距離。

挑選采集的數(shù)據(jù)點(diǎn)并將其進(jìn)行擬合，在同一坐標(biāo)系中將磨削前螺桿轉(zhuǎn)子曲面數(shù)據(jù)點(diǎn)坐標(biāo)減去磨削后螺桿轉(zhuǎn)子曲面坐標(biāo)即可得到實(shí)際磨削深度，砂帶實(shí)際磨削深度與計(jì)算值對(duì)比結(jié)果見(jiàn)表4。

通過(guò)對(duì)比螺桿轉(zhuǎn)子磨削前后表面輪廓測(cè)量得到的磨削深度與數(shù)學(xué)模型計(jì)算的磨削深度發(fā)現(xiàn)，試驗(yàn)8組數(shù)據(jù)中最大誤差為9.6%，最小誤差為4.1%。分析誤差主要來(lái)源于兩方面，一方面是磨粒磨削初始時(shí)高度的分布狀態(tài)并非完全符合高斯分布，在磨損過(guò)程中對(duì)工件進(jìn)行切削的磨粒頂部不是規(guī)則的球形。另一方面則是由于螺桿轉(zhuǎn)子與接觸輪實(shí)際接觸區(qū)域的壓力分布與有限元分析存在誤差。但模型的整體準(zhǔn)確率保持在10%以內(nèi)，由此證明了建立的螺桿轉(zhuǎn)子材料去除率數(shù)學(xué)模型的可行性。

5 結(jié)論

（1）對(duì)砂帶表面磨粒出刃高度進(jìn)行測(cè)量，探究砂帶磨粒的磨損規(guī)律，利用計(jì)算機(jī)編程建立砂帶表面三維數(shù)字化模型，解決了砂帶表面形貌難以重塑的問(wèn)題。

（2）利用有限元仿真軟件對(duì)螺桿轉(zhuǎn)子與接觸輪接觸區(qū)域的壓力進(jìn)行仿真，獲得接觸區(qū)域壓力分布情況，并利用單顆磨粒磨削仿真的方式明確磨削歷程中磨削力的變化規(guī)律。

（3）將砂帶磨粒簡(jiǎn)化為圓頂圓錐，利用磨粒的彈塑性變形、數(shù)值分析等理論建立了基于砂帶磨損的材料去除率數(shù)學(xué)模型。通過(guò)試驗(yàn)驗(yàn)證，該模型的誤差在10%以內(nèi)，可以為螺桿轉(zhuǎn)子砂帶磨削表面的均勻去除提供理論基礎(chǔ)。

參考文獻(xiàn)：

［1］ 黃云.砂帶磨削技術(shù)的研究現(xiàn)狀和發(fā)展方向簡(jiǎn)介[J].金剛石與磨料磨具工程，2020，40（3）：1-4.

HUANG Yun. Brief introduction of research status and development direction of belt grinding technology [J]. Diamond amp; Abrasives Engineering，2020，40（3）：1-4.

［2］ 潘杰， 陳凡， 楊煒， 等.基于 SPSO–BP 神經(jīng)網(wǎng)絡(luò)的自適應(yīng)拋光工藝參數(shù)匹配[J].表面技術(shù)，2022，51（8）：387-399.

PAN Jie， CHEN Fan， YANG Wei， et al. Adaptive polishing process parameter matching based on SPSO-BP neural network [J]. Surface technology，2022，51（8）：387-399.

［3］ HE Z， LI J Y， LIU Y M， et al. Investigating the effects of contact pressure on rail material abrasive belt grinding performance [J]. International Journal of Advanced Manufacturing Technology，2017，93（10）：1-8.

［4］ 遲玉倫， 俞鑫， 劉斌， 等.基于軸承套圈毛坯表面分析的磨削材料去除率模型與應(yīng)用實(shí)驗(yàn)研究[J].表面技術(shù)，2023，52（4）：338-353，373.

CHI Yulun， YU Xin， LIU Bin， et al. Research on grinding material removal rate model and application experiment based on bearing ring blank surface analysis [J]. Surface technology，2023，52（4）：338-353，373.

［5］ 劉斐， 王偉， 王雷， 等.接觸輪變形對(duì)機(jī)器人砂帶磨削深度的影響[J].機(jī)械工程學(xué)報(bào)，2017，53（5）：86-92.

LIU Fei， WANG Wei， WANG Lei， et al. Effect of Contact Wheel’s Deformation on Cutting Depth for Robotic Belt Grinding [J]. Journal of mechanical engineering，2017，53（5）：86-92.

［6］齊俊德， 陳冰.考慮單磨粒作用的砂帶磨削機(jī)理模型[J].金剛石與磨料磨具工程，2020，40（3）：13-20.

QI Junde， CHEN Bing. Mechanism model of belt grinding considering single abrasive action [J]. Diamond amp; Abrasives Engineering，2020，40（3）：13-20.

［7］WANG Y J， HUANG Y， CHEN Y， et al. Model of an abrasive belt grinding surface removal contour and its application [J]. The International Journal of Advanced Manufacturing Technology，2016，82（9/10/11/12）：2113-2122.

［8］劉偉， 劉仁通， 鄧朝暉， 等.單顆磨粒磨削仿真研究進(jìn)展[J].宇航材料工藝，2018，48（4）：1-8.

LIU Wei， LIU Rentong， DENG Chaohui， et al. Research progress no single abrasive grain grinding simulation [J]. Aerospace materials amp; technology，2018，48（4）：1-8.

［9］王榮全.面向鋼軌打磨的砂帶磨削過(guò)程建模與實(shí)驗(yàn)研究[D].北京：北京交通大學(xué)， 2016.

WANG Rongquan. The modeling and experimental research of belt- grinding process in rail grinding[D]. Beijing： Beijing Jiaotong of University， 2016.

［10］杜瑤， 田彪， 張健， 等.中密度纖維板磨削時(shí)的砂帶磨損研究及壽命評(píng)判[J].北京林業(yè)大學(xué)學(xué)報(bào)，2021，43（10）：126-134.

DU Yao， TIAN Biao， ZHANG Jian， et al. Research on abrasive belt wear and working life evaluation [J]. Journal of Beijing forestry university，2021，43（10）：126-134.

［11］ KARKALO N E， MARKOPOULOS A P， KUNDRáK J. Molecular dynamics model of nano-metric peripheral grinding [J]. Procedia Cirp，2017（58）：281-286.

［12］楊雅琦， 汪生瑩， 張亦鳴， 等.棕剛玉磨料砂帶磨削木材的磨損機(jī)制及其壽命[J].林業(yè)工程學(xué)報(bào)， 2020， 5（6）：121-127.

YANG Yaqi， WANG Shengying， ZHANG Yiming， et al. Wear mechanism and life of brown corundum abrasive belt in wood sanding process[J]， Journal of forestry engineering， 2020， 5（6）：121-127.

［13］ BLUNT L， EBDON S. The application of three-dimensional surface measurement techniques to characterizing grinding wheel topography [J].International Journal of Machine Tools and Manufacture，1996，36（11）：1027-1126.

［14］ JOURANI A， HAGèGE B， BOUVIER S， et al. Influence of abrasive grain geometry on friction coefficient and wear rate in belt finishing [J]. Tribology International，2013（59）：30-37.

作者簡(jiǎn)介

陳斯睿，男，1998年生，碩士研究生。主要研究方向：復(fù)雜曲面精密加工。

Email： a13897990910@163.com

通信作者：孫興偉，女，1970年生，教授，博士研究生導(dǎo)師。主要研究方向：數(shù)控裝備及理論。

Email：ww7247@126.com

Analysis of the influence of abrasive belt wear on the removal depth ofgrinding screw rotor

CHEN Sirui1，2，SUN Xingwei1，2，YANG Heran1，2，LIU Yin1，2

（1. School of Mechanical Engineering， Shenyang University of Technology， Shenyang 110870， China）

（2. Key Laboratory of Numerical Control Manufacturing Technology for Complex Surfaces ofLiaoning Province， Shenyang 110870， China）

Abstract To ensure the uniformity of material removal during surface grinding of screw rotors， a mathematical model for the material removal depth of screw rotors is established， considering the impact of abrasive wear on grinding re- moval depth. Based on the theory of elastic-plastic deformation and the wear pattern of abrasive particles during grind- ing， a theoretical model for the removal rate of single abrasive material is established. The wear of abrasive particles is divided into two stages： the rapid wear stage and the stable wear stage. Based on the number of abrasive particles per unit contact area and the rule of abrasive cutting edge， the maximum cutting depth of abrasive particles is calculated by numerical integration and contact pressure， and the mathematical model of macroscopic material removal rate based on abrasive belt wear is established. The model consideres the wear of the sand belt. The maximum error of the predicted value is 9.6%， and the minimum error is 4.1% compared with the experimental value. The experimental data fully veri- fy the effectiveness of the model， and can provide a theoretical basis for ensuring the profile accuracy of the screw rotor.

Key words belt wear；mathematical model；removal depth；screw rotor；model verification