王樹(shù)強(qiáng)，陳釗，陳昊雷，周游

（沈陽(yáng)化工大學(xué)機(jī)械與動(dòng)力工程學(xué)院，遼寧 沈陽(yáng) 110142）

機(jī)械鍍鋅是在常溫常壓條件下利用化學(xué)吸附沉積和物理碰撞使金屬鋅粉在工件表面成膜的綠色表面處理工藝[1-3]。典型的機(jī)械鍍鋅工藝是把經(jīng)過(guò)前處理的工件放入轉(zhuǎn)動(dòng)的鍍筒中，根據(jù)預(yù)定的鍍層厚度加入適量鋅粉、活性劑、沖擊介質(zhì)和水，形成碰撞和搓碾的流態(tài)環(huán)境，在藥劑和機(jī)械碰撞的共同作用下工件表面逐漸形成鍍層。整個(gè)鍍層的形成過(guò)程可以簡(jiǎn)化為“聚集→吸附→沉積→鑲嵌→致密化”[4]。離散單元法(DEM)是一種用以解決非連續(xù)介質(zhì)問(wèn)題的數(shù)值模擬方法[5]。它以求解空間離散元為單元，能夠求解所有的元在任意時(shí)刻的位移、速率、加速度等物理量，進(jìn)一步預(yù)測(cè)離散群體行為，被廣泛應(yīng)用在巖土、建筑、農(nóng)業(yè)、醫(yī)藥等領(lǐng)域顆粒體系的過(guò)程模擬、裝備設(shè)計(jì)等方面[6-8]。

為了使模擬數(shù)據(jù)更接近實(shí)際，本文采用離散單元法中的Hertz-Mindlin with JKR Cohesion接觸模型從微觀角度模擬分析了機(jī)械鍍鋅層形成過(guò)程中的顆粒運(yùn)動(dòng)規(guī)律，研究了在不同沖擊速率下顆粒群的運(yùn)動(dòng)規(guī)律，并通過(guò)實(shí)驗(yàn)驗(yàn)證得出了最佳沖擊速率。

1 模型的建立

1.1 仿真模型的選擇

在機(jī)械鍍鋅的鍍層致密化過(guò)程中，鋅顆粒間存在相互粘結(jié)的現(xiàn)象。常用的Hertz-Mindlin接觸模型不考慮顆粒之間的粘結(jié)力，僅考慮彈性變形；Hertz-Mindlin with Bonding接觸模型常用于模擬眾多小顆粒之間相互粘結(jié)成為大顆粒，并在外部作用力下小顆粒間的粘結(jié)發(fā)生斷裂、破碎的過(guò)程，該模型僅適用于有限粘結(jié)力大小的粘結(jié)[9]。而Hertz-Mindlin with JKR Cohesion接觸模型是一種建立在Hertz理論上的粘結(jié)性顆粒接觸模型，該模型考慮了顆粒間的粘結(jié)力、顆粒運(yùn)動(dòng)等因素，非常適用于顆粒間因存在液橋力、范德華力等而發(fā)生粘結(jié)和團(tuán)聚的情形[10]。在機(jī)械鍍鋅過(guò)程中，鋅粉顆粒受液橋力、鋅粉間作用力、滾筒內(nèi)動(dòng)力場(chǎng)和活性劑的影響，在鍍液中形成細(xì)小藻團(tuán)狀聚集體，并在玻璃珠的沖擊作用下被撞擊和粘結(jié)到工件表面，形成一層致密的鍍層[11]。因此，Hertz-Mindlin with JKR Cohesion接觸模型適用于模擬機(jī)械鍍鋅的鍍層形成過(guò)程。其JKR Cohesion簡(jiǎn)化模型如圖1所示，模型中的法向彈性力能夠反映顆粒間的粘彈性特征，由表面能和法向重疊量決定，具體見(jiàn)式(1)-(4)。

圖1 JKR粘結(jié)模型示意圖Figure 1 Schematic diagram of JKR cohesion model

式中FJKR為JKR法向彈性力(單位：N)，γ為表面能(單位：J/m2)，E*為等效彈性模量(單位：Pa)，δ為2個(gè)接觸鋅粉顆粒之間的切向重疊量(單位：m)，α為2個(gè)接觸鋅粉顆粒之間的法向重疊量(單位：m)，R*為等效接觸半徑(單位：m)，E1、E2為接觸鋅粉顆粒的彈性模量(單位：Pa)，μ1、μ2為接觸鋅粉顆粒的泊松比，R1、R2為接觸鋅粉顆粒的接觸半徑(單位：m)。

1.2 仿真參數(shù)選擇

設(shè)置模型的仿真參數(shù)如下：鋅粉，密度7 100 kg/m3，泊松比0.25，剪切模量3.94 × 1010Pa；Q235鋼基板，密度7 850 kg/m3，泊松比0.3，剪切模量8.0 × 1010Pa；玻璃作為沖擊介質(zhì)，密度2 500 kg/m3，泊松比為0.25，剪切模量為2.2 × 1010Pa。

為了研究沖擊速率對(duì)鍍層形成及鋅粉運(yùn)動(dòng)的影響，設(shè)置模型中的鋅粉為靜態(tài)生成和自然填充，顆粒直徑在 6 ～ 8 μm 范圍內(nèi)隨機(jī)分布；工件模型設(shè)為 200 μm × 200 μm × 200 μm 的立方體，空間網(wǎng)格尺寸為 9 μm。如圖2所示，在模擬過(guò)程中將初始模型分為5層，每層40 μm高。

圖2 模型初始狀態(tài)示意圖Figure 2 Sketch showing the initial state of model

2 結(jié)果與討論

2.1 運(yùn)用EDEM離散元軟件模擬機(jī)械鍍鋅的可行性分析

從圖3可以看出，無(wú)論是運(yùn)用EDEM離散元軟件模擬所得，還是機(jī)械鍍所得真實(shí)鋅層，01層的鋅顆粒均被嚴(yán)重壓縮。模擬結(jié)果的01層厚度約為12 μm，真實(shí)鋅層的01層厚度約為16 μm。模擬和真實(shí)機(jī)械鍍鋅層的相對(duì)致密度分別為76.6%和80.1%?？梢?jiàn)運(yùn)用EDEM離散元軟件分析機(jī)械鍍鋅的鍍層成型過(guò)程可行。

圖3 機(jī)械鍍鋅層截面的模擬示意圖(a)和真實(shí)形貌(b)Figure 3 Simulated sketch (a) and actual image (b) of cross-section of mechanically plated zinc coating

2.2 單一沖擊速率下鋅粉的運(yùn)動(dòng)特性分析

從圖4和圖5可知，機(jī)械鍍鋅的鍍層形成過(guò)程可分為3個(gè)階段。

圖4 在1.0 m/s沖擊速率下鋅粉群在縱向(a)和橫向(b)上的移動(dòng)速率變化Figure 4 Vertical (a) and horizontal (b) movement rates of zinc particles at an impact velocity of 1.0 m/s

圖5 在1.0 m/s沖擊速率下不同時(shí)間段的鋅顆粒間作用力Figure5 Interparticle forces between zinc particles at different processing stages when the impact velocity is 1.0 m/s

(1) 沖擊前期(0 ～ 0.07 ms)：屬于鋅粉群下落階段。在縱向上，01層鋅粉群最先受到?jīng)_擊力并獲得最大縱向移動(dòng)速率，沖擊能量減弱，05層鋅粉群最后獲得最小的縱向速率；在橫向上，鋅粉群未受到?jīng)_擊力的作用，橫向移動(dòng)速率無(wú)變化。

(2) 沖擊中期(0.07 ～ 0.16 ms)：屬于鋅粉群間隙壓縮階段。在縱向上，05層部分鋅粉與基體接觸，使得鋅粉群的縱向移動(dòng)速率減小，產(chǎn)生的反作用力不斷減弱，故05層的鋅粉縱向移動(dòng)速率變化最大，01層的鋅粉縱向移動(dòng)速率變化最小，最終各層的縱向移動(dòng)速率達(dá)到相對(duì)穩(wěn)定的狀態(tài)；在橫向上，鋅粉群仍未受到作用力，橫向移動(dòng)速率保持不變。

(3) 沖擊后期(0.16 ～ 0.21 ms)：屬于鋅粉群重排變形階段。在沖擊力和鋅粉間作用力(見(jiàn)圖5c和圖5d)的共同作用下，所有鋅粉都進(jìn)入波動(dòng)狀態(tài)，發(fā)生重排和壓縮變形。

根據(jù)上述分析可知，在沖擊前期和中期主要在縱向上對(duì)鋅粉間隙進(jìn)行壓縮，沖擊后期鋅粉在縱向和橫向上都有波動(dòng)，鋅粉顆粒間的重排和壓縮變形不斷進(jìn)行。

2.3 不同沖擊速率下鋅粉的運(yùn)動(dòng)特性分析

為研究沖擊速率對(duì)沖擊效果的影響，在保證沖擊行程相同的條件下，改變沖擊速率和沖擊時(shí)間(見(jiàn)表1)，對(duì)鍍層的形成進(jìn)行模擬分析，結(jié)果見(jiàn)圖6。

表1 沖擊速率與對(duì)應(yīng)的沖擊時(shí)間Table 1 Correspondence between impact velocity and impact time

圖6 不同沖擊速率下對(duì)應(yīng)的縱向(a)和橫向(b)速率的變化Figure 6 Variation of vertical (a) and horizontal (b) velocity under different impact velocities

圖6中的3組數(shù)據(jù)與圖4具有相似的變化規(guī)律，但各自的移動(dòng)速率波動(dòng)差異較大。在圖6a中，沖擊速率為1.0 m/s和1.4 m/s時(shí)，顆粒群的縱向移動(dòng)速率波動(dòng)最大；在圖6b中，沖擊速率為1.0 m/s時(shí)，顆粒群的橫向移動(dòng)速率波動(dòng)最大。波動(dòng)程度越大，鋅顆粒間的重排效果越好。故沖擊速率為1.0 m/s時(shí)的重排效果最佳，所得鋅鍍層性能應(yīng)該最好。

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

為驗(yàn)證模擬分析的可靠性，以?xún)?nèi)徑16 mm、外徑28 mm的圓形Q235鋼墊片作為基材，玻璃珠作為沖擊介質(zhì)，1 200目鋅粉和相關(guān)表面活性劑為實(shí)驗(yàn)材料，在自主設(shè)計(jì)的50L型機(jī)械鍍鋅機(jī)上進(jìn)行機(jī)械鍍鋅試驗(yàn)，并按照標(biāo)準(zhǔn)ISO 2409:2020Paints and Varnishes — Cross-cut Test對(duì)鍍鋅件進(jìn)行劃格試驗(yàn)。

從圖7可知，當(dāng)沖擊速率為0.6 m/s時(shí)，沖擊力不足，鋅粉過(guò)度聚集，重排效果不佳，鋅鍍層表面存在凸點(diǎn)，附著力為1級(jí)；當(dāng)沖擊速率為1.0 m/s時(shí)，鋅鍍層表面光滑平整，附著力為0級(jí)；當(dāng)沖擊速率為1.4 m/s時(shí)，鋅鍍層整體光滑平整，但由于沖擊速率過(guò)大，局部出現(xiàn)劃痕和點(diǎn)狀脫落，附著力為1級(jí)?？梢?jiàn)較佳的沖擊速率為1.0 m/s，與模擬分析結(jié)果吻合。

圖7 不同沖擊速率下機(jī)械鍍鋅試樣劃格試驗(yàn)的結(jié)果Figure 7 Cross-cut test results of zinc coatings mechanically plated at different impact velocities

3 結(jié)論

(1) 在機(jī)械鍍鋅過(guò)程中，位于上層的鋅粉發(fā)生形變和重排，位于中下層的鋅粉則以重排為主。

(2) 在機(jī)械鍍鋅的前期，縱向上鋅粉逐層下落，橫向上無(wú)變化；中期鋅粉的縱向空隙得以壓縮，橫向上無(wú)變化；后期鋅粉在縱、橫方向上都發(fā)生重排和壓縮變形。

(3) 沖擊速率為1.0 m/s時(shí)，所得鋅鍍層光滑、平整，附著力為0級(jí)，模擬分析結(jié)果與試驗(yàn)結(jié)果基本吻合。