齊永超 呂玉山 李興山 許劉宛 衣軍任

(沈陽(yáng)理工大學(xué)機(jī)械工程學(xué)院 遼寧沈陽(yáng) 110159)

本文作者針對(duì)凹坑結(jié)構(gòu)化表面，提出拓?fù)淠ハ鞯募庸し椒?，依?jù)工件的結(jié)構(gòu)化表面特征設(shè)計(jì)結(jié)構(gòu)化砂輪，建立磨削運(yùn)動(dòng)模型，并進(jìn)行磨削過(guò)程仿真研究。

1 結(jié)構(gòu)化表面拓?fù)涮卣鞣治?/h2>

1.1 結(jié)構(gòu)單元形狀拓?fù)涮卣鞣治?/h3>

研究表明，凹坑型結(jié)構(gòu)化表面上結(jié)構(gòu)單元的形狀尺寸和排布影響著減阻減磨性能[11-12]。土壤動(dòng)物蜣螂頭部和上唇基分布有圓形及橢圓形的凹坑，具有很好的減黏降阻作用[13]。文中即以橢球冠及球冠形凹坑為例，建立公式(1)所示的結(jié)構(gòu)單元數(shù)學(xué)模型。

(1)

式中：n=h/c，c為橢球冠所在橢圓z軸半徑；l、w、h分別為結(jié)構(gòu)單元長(zhǎng)度、寬度和深度。

當(dāng)l>w時(shí),結(jié)構(gòu)單元形狀為橢球冠形,當(dāng)l=w時(shí)，結(jié)構(gòu)單元形狀為球冠形。圖1為結(jié)構(gòu)單元模型圖。

圖1 凹坑結(jié)構(gòu)單元模型

由數(shù)學(xué)模型(1)可知，參數(shù)l、w、h決定結(jié)構(gòu)單元的尺寸形狀，由此提取其為拓?fù)涮卣鲄⒘?，建立結(jié)構(gòu)單元拓?fù)涮卣骶仃嚍?/p>

(2)

1.2 結(jié)構(gòu)單元排布拓?fù)涮卣鞣治?/h3>

依據(jù)生物表面結(jié)構(gòu)單元排布方式，建立陣列、錯(cuò)位、葉序3種排布規(guī)律的數(shù)學(xué)模型，圖2為單元排布示意圖。定義單元長(zhǎng)度所在的工件與砂輪相對(duì)運(yùn)動(dòng)方向?yàn)楣ぜ潭ㄗ鴺?biāo)系xw軸方向，寬度方向?yàn)閥w軸方向，垂直工件表面方向?yàn)閦w軸方向，被磨削平面為zw=0平面。xw方向排布Iw列，yw方向排布Jw行。則陣列排布數(shù)學(xué)表達(dá)式：

圖2 結(jié)構(gòu)單元排布模型

(3)

錯(cuò)位排布數(shù)學(xué)表達(dá)式為

(4)

式中：i=1,2,3，......,Iw；j=1,2,3,......,Jw;(xi，yj)為任意結(jié)構(gòu)單元中心點(diǎn)坐標(biāo);Tx、Ty分別為xw、yw方向排布周期；φx為xw方向相位差。

葉序排布數(shù)學(xué)表達(dá)式為

(5)

式中：R為Van Iterson模型[14-15]中母體圓柱半徑；c為葉序排布生長(zhǎng)序數(shù)；葉序角θ=137.508°,k=1,2,3，......,Kw，為結(jié)構(gòu)單元序數(shù)；m為控制第k個(gè)結(jié)構(gòu)單元在xw軸方向上的排布序數(shù)；φy為yw方向相位差。

由公式(3)—(5)可知，結(jié)構(gòu)化表面的排布特征由排布周期與相位差決定，陣列排布由排布周期Tx和Ty控制結(jié)構(gòu)單元排布間距，行列之間整齊排列，相位差為0；錯(cuò)位排布xw方向與yw方向排布周期分別為Tx、Ty,偶數(shù)行與奇數(shù)行之間相位差為φx；葉序排布xw方向有相位差φx，yw方向有相位差φy。則提取3種排布的周期與相位建立特征向量矩陣為

[Twp]=[TxTyφxφy]T

(6)

2 結(jié)構(gòu)化砂輪設(shè)計(jì)

2.1 凹坑結(jié)構(gòu)化表面創(chuàng)成機(jī)制

圖3所示為單元長(zhǎng)度方向磨粒簇中心截面磨削過(guò)程示意圖，磨料簇在砂輪周向間隔排布。在工件上弧AB段為磨削階段，砂輪上A′B′段為磨粒簇與工件接觸磨削出的凹坑單元；BC段為空閑階段，砂輪上B′C′段不磨削工件，在結(jié)構(gòu)單元間形成平臺(tái)間隔。

圖3 截面磨削過(guò)程示意

中心截面各磨粒運(yùn)動(dòng)軌跡如圖4所示，磨削過(guò)程中磨料簇內(nèi)各個(gè)磨粒相對(duì)工件做延長(zhǎng)外擺線運(yùn)動(dòng)，各磨粒運(yùn)動(dòng)軌跡在工件表面包絡(luò)形成結(jié)構(gòu)單元輪廓。磨料簇最低點(diǎn)A′從切入點(diǎn)開(kāi)始進(jìn)入磨削階段，由于磨料簇高度逐漸增大，高點(diǎn)磨粒從切入點(diǎn)左側(cè)開(kāi)始切削，形成超越切入點(diǎn)的長(zhǎng)度l1，且存在某一點(diǎn)磨粒，超越切入點(diǎn)距離最大為l1max。l2為磨料簇與工件接觸時(shí)長(zhǎng)內(nèi)工件進(jìn)給的距離，l3為超越切出點(diǎn)長(zhǎng)度，與l1相等。設(shè)定砂輪磨料簇高度變化為f(α),磨料簇長(zhǎng)度為L(zhǎng)，最大高度差為H，周向排布間隔φs，砂輪基體半徑為rs，平均磨粒直徑為dg，磨削深度為ap，則長(zhǎng)度l1為

圖4 凹坑中心截面形貌

(7)

式中:α為進(jìn)入磨削狀態(tài)的磨料簇的圓心角度;β為不同α值時(shí)，對(duì)應(yīng)位置磨粒與平面開(kāi)始接觸時(shí)的角度;ps為3種排布的速度比。從而得到結(jié)構(gòu)單元長(zhǎng)度為

(8)

結(jié)構(gòu)單元間的排布周期為

(9)

由運(yùn)動(dòng)學(xué)關(guān)系可知，當(dāng)Tx=l時(shí)，結(jié)構(gòu)單元相接，凹坑相連;Tx>l時(shí)，結(jié)構(gòu)單元相離，可形成凹坑型結(jié)構(gòu)化表面。

2.2 工件與砂輪拓?fù)溆成潢P(guān)系模型

設(shè)定工件固定坐標(biāo)系為工件三維拓?fù)淇臻g，砂輪固定坐標(biāo)系為砂輪三維拓?fù)淇臻g，工件表面結(jié)構(gòu)單元與砂輪磨料簇的拓?fù)鋵?duì)應(yīng)關(guān)系中，在寬度方向與深度方向?yàn)榈戎底儞Q，在長(zhǎng)度方向，由公式(8)可知，受磨料簇形狀、速度比、磨削深度等影響，發(fā)生拉伸和彎曲變形，且結(jié)構(gòu)單元集合點(diǎn)與磨料簇集合點(diǎn)間存在一一映射關(guān)系。設(shè)定砂輪磨料簇拓?fù)涮卣骶仃嚺c排布特征矩陣為分別為

[Tsd]=[LWH1]T

(10)

[Tsp]=[σxσyδxδy]T

(11)

式中:L、W、H分別為磨料簇的長(zhǎng)度、寬度與高度;σx、σy、δx、δy分別為周向與軸向排布周期和相位差。

由這種拓?fù)渥儞Q關(guān)系，建立工件結(jié)構(gòu)單元與砂輪磨粒簇之間的拓?fù)溆成潢P(guān)系矩陣[Cd]與排布特征拓?fù)溆成渚仃嘯Cp]：

(12)

(13)

則工件與砂輪之間的拓?fù)溆成潢P(guān)系為

(14)

2.3 結(jié)構(gòu)化砂輪數(shù)學(xué)模型

由公式(14)得到砂輪的單元特征與排布特征，在長(zhǎng)度方向?qū)⒔Y(jié)構(gòu)單元長(zhǎng)半軸等分為S份，由公式(1)得到磨料簇離散數(shù)學(xué)模型-f(s·L/(2S),s·W/(2S))，s=1,2,......,S，得到磨料簇?cái)?shù)學(xué)模型中各磨粒點(diǎn)坐標(biāo)(xs,ys,zs)。設(shè)定陣列與錯(cuò)位排布砂輪表面磨粒簇周向排列Is列，軸向排列Js行，i=1,2,......,Is,j=1,2,......,Js。對(duì)于陣列排布有周向周期σx、軸向周期σy、相位差為0，由此建立陣列排布磨料簇砂輪數(shù)學(xué)模型：

(15)

錯(cuò)位排布奇數(shù)行與偶數(shù)行間有周向的相位差δx=σx/2，j為奇數(shù)時(shí)，數(shù)學(xué)模型與公式(15)相同，j為偶數(shù)時(shí)砂輪數(shù)學(xué)模型為

(16)

葉序排布在砂輪表面排布Ks個(gè)磨粒簇，k=1,2,3,......,Ks，周向相位差δx=θ，軸向相位差δy為葉序生長(zhǎng)序數(shù)c，則數(shù)學(xué)模型為

(17)

由拓?fù)溆成潢P(guān)系公式(14)和砂輪的數(shù)學(xué)模型(15)—(17)設(shè)計(jì)出砂輪，使用MatLab建立圖5所示的陣列、錯(cuò)位、葉序排布砂輪模型圖，砂輪由砂輪基體與規(guī)則排布的磨料簇組成。

圖5 磨料簇砂輪模型

3 磨粒運(yùn)動(dòng)軌跡模型

磨料簇砂輪逆磨平面結(jié)構(gòu)化表面，在磨削運(yùn)動(dòng)仿真中工件靜止，砂輪相對(duì)工件運(yùn)動(dòng)，砂輪逆時(shí)針勻速轉(zhuǎn)動(dòng)，轉(zhuǎn)速為ωs=2πns/60，砂輪相對(duì)工件勻速進(jìn)給，進(jìn)給速度為vw，砂輪與工件坐標(biāo)系之間的距離為e=r+H-ap。數(shù)學(xué)模型以時(shí)間t為變量，控制磨削運(yùn)動(dòng)過(guò)程，圖6為磨削過(guò)程示意圖。

圖6 磨削過(guò)程示意

定義砂輪固定坐標(biāo)系Ss(xs,ys,zs)，砂輪旋轉(zhuǎn)坐標(biāo)系Sn(xn,yn,zn),工件進(jìn)給坐標(biāo)系Sf(xf,yf,zf)，工件絕對(duì)坐標(biāo)系Sw(xw,yw,zw)。通過(guò)坐標(biāo)變換實(shí)現(xiàn)從砂輪固定坐標(biāo)到工件固定坐標(biāo)系的轉(zhuǎn)換，建立砂輪固定坐標(biāo)系到砂輪旋轉(zhuǎn)坐標(biāo)的變換矩陣[Msn]，砂輪旋轉(zhuǎn)坐標(biāo)系到工件進(jìn)給變換矩陣[Mnf]，工件進(jìn)給坐標(biāo)系到工件固定坐標(biāo)系變換矩陣[Mfw]，分別為

(18)

(19)

(20)

由變換矩陣(18)—(20)得到磨料簇砂輪上各個(gè)磨粒在不同時(shí)刻運(yùn)動(dòng)軌跡的數(shù)學(xué)模型：

[Q]=[Mfw]×[Mnf]×[Msn]×[P]

(21)

4 仿真結(jié)果與分析

根據(jù)砂輪數(shù)學(xué)模型和運(yùn)動(dòng)軌跡模型使用MatLab進(jìn)行仿真。仿真中砂輪基體半徑為50 mm，軸向?qū)挾葹?3 mm，磨粒平均直徑為95 μm，磨料簇弧長(zhǎng)為30 mm,磨料簇寬度為2 mm，磨料簇最大高度差為20 μm，設(shè)定砂輪轉(zhuǎn)速為1 000 r/min，通過(guò)控制工件進(jìn)給速度來(lái)控制轉(zhuǎn)速比，研究不同的排布規(guī)律、速度比與磨削深度對(duì)結(jié)構(gòu)化表面形貌的影響。

4.1 排布對(duì)結(jié)構(gòu)化表面形貌的影響

對(duì)于陣列排布砂輪，磨料簇軸向排布4行，周期3 mm，周向排布4列；錯(cuò)位排布軸向7行，周期1.5 mm，周向排布4列，偶數(shù)與奇數(shù)行間相位差1.5 mm；葉序排布軸向相位差0.8 mm，周向相位差為葉序角θ，排布12個(gè)磨料簇。設(shè)置3種排布的速度比ps分別為12.9、9、9，磨削深度ap=20 μm。仿真形貌分別如圖7所示，通過(guò)改變結(jié)構(gòu)單元的排布拓?fù)涮卣鱗Tw]，可得到不同排布的結(jié)構(gòu)化表面形貌，實(shí)現(xiàn)結(jié)構(gòu)化表面形貌的多樣化。此時(shí)陣列排布為相接的臨界狀態(tài)，而錯(cuò)位排布與葉序排布可在更小的速度比下，形成相離狀態(tài)的凹坑結(jié)構(gòu)化表面，具有更高的表面占用率。

圖7 不同排布形貌仿真結(jié)果

4.2 速度比對(duì)結(jié)構(gòu)化表面形貌的影響

以陣列排布為例，研究速度比對(duì)結(jié)構(gòu)化表面形貌的影響。設(shè)置磨削深度為20 μm不變，圖7(a)和圖8(a)、(b)是轉(zhuǎn)速比ps分別在12.9、14和16時(shí)工件表面的形貌，結(jié)構(gòu)單元的長(zhǎng)度分別為3.18、3.22、3.38 mm，結(jié)構(gòu)單元的排布間距分別為3.18、3.5、4 mm。由仿真結(jié)果可知，當(dāng)l=Tx時(shí)，凹坑單元相接處于臨界狀態(tài)，隨著速度比進(jìn)一步增大，結(jié)構(gòu)單元深度h與寬度w不變，結(jié)構(gòu)單元長(zhǎng)度l與周期Tx均增大。且進(jìn)一步研究得，在ps較小時(shí)，長(zhǎng)度l增長(zhǎng)速率較小，ps較大時(shí)增長(zhǎng)速率較大且趨于線性增長(zhǎng)；排布周期Tx與轉(zhuǎn)速比ps成正比例關(guān)系，比例系數(shù)為σx/(2π)，且增大的速率大于單元長(zhǎng)度的增長(zhǎng)速率。

圖8 不同轉(zhuǎn)速比形貌仿真結(jié)果

4.3 磨削深度對(duì)結(jié)構(gòu)化表面形貌的影響

在速度比一定時(shí)，研究磨削深度對(duì)結(jié)構(gòu)化表面形貌的影響。設(shè)置速度比為12.9，圖9(a)、(b)和圖7(a)是磨削深度ap分別在10、15和20 μm時(shí)工件表面的形貌，單元的長(zhǎng)度l分別為2.25、2.77、 3.22 mm，單元的排布間距Tx均為3.22 mm。隨著磨削深度的增大，單元排布周期Tx值不變，結(jié)構(gòu)單元長(zhǎng)度l與寬度w均增大，且增長(zhǎng)速率受單元形狀影響均逐漸減小。

圖9 不同磨削深度形貌仿真結(jié)果

5 結(jié)論

(1)基于磨削運(yùn)動(dòng)學(xué)關(guān)系與點(diǎn)集拓?fù)淅碚?，提出拓?fù)淠ハ鹘Y(jié)構(gòu)化減阻表面的方法，該方法可根據(jù)工件與砂輪的拓?fù)溆成潢P(guān)系設(shè)計(jì)磨料簇砂輪，實(shí)現(xiàn)結(jié)構(gòu)化表面形貌多樣化，從而加工出特定單元形狀與排布規(guī)律的凹坑型結(jié)構(gòu)化表面，達(dá)到預(yù)期的減阻、減磨性能。

(2)通過(guò)建立工件與砂輪拓?fù)淇臻g的映射關(guān)系矩陣，設(shè)計(jì)出磨粒簇砂輪并建立磨削運(yùn)動(dòng)軌跡模型，進(jìn)行拓?fù)淠ハ鬟^(guò)程仿真，仿真結(jié)果表明：工件與砂輪的速度比ps與磨削深度ap是影響結(jié)構(gòu)單元形貌的主要因素。磨削深度一定時(shí)，工件與砂輪的速度比增大，結(jié)構(gòu)單元的長(zhǎng)度l逐漸增大，排布周期Tx隨速度比線性增大，當(dāng)結(jié)構(gòu)單元長(zhǎng)度與排布周期Tx相等時(shí)，凹坑單元相接；單元寬度w與周期Ty不變。在速度比一定時(shí)，磨削深度越大，結(jié)構(gòu)單元的長(zhǎng)度與單元的寬度均越大，且結(jié)構(gòu)單元的排布周期與磨削深度無(wú)關(guān)。

(3)仿真結(jié)果驗(yàn)證了拓?fù)淠ハ鹘Y(jié)構(gòu)化表面方法的可行性，后續(xù)將進(jìn)行磨料簇砂輪的制造與磨削實(shí)驗(yàn)，對(duì)上述結(jié)論進(jìn)行實(shí)驗(yàn)驗(yàn)證。