紡織鋼領(lǐng)類金剛石涂層耐沖擊性能仿真研究

蘇慧明， 郭飛飛， 楊軼飛， 王珂心

(西安工程大學(xué) 機電工程學(xué)院， 陜西 西安 710048)

環(huán)錠細紗機運轉(zhuǎn)時，其鋼絲圈繞鋼領(lǐng)高速旋轉(zhuǎn)，產(chǎn)生摩擦磨損的同時，紡線的張力和摩擦力也會使鋼絲圈偏移和跳躍[1]，由此產(chǎn)生的瞬時沖擊效應(yīng)不可忽視，因此提升鋼領(lǐng)的抗沖擊性能可以有效降低由于沖擊導(dǎo)致的鋼領(lǐng)失效問題，從而提升鋼領(lǐng)的使用壽命。運用表面處理技術(shù)是提高鋼領(lǐng)使用壽命的主要方法之一，包括光面鋼領(lǐng)、鍍鉻鋼領(lǐng)、化學(xué)鍍鎳磷鋼領(lǐng)、超聲波減磨鋼領(lǐng)、陶瓷鋼領(lǐng)以及DLC涂層鋼領(lǐng)[2-5]。類金剛石涂層是一種亞穩(wěn)態(tài)的非晶碳膜，碳原子以sp3和sp2鍵結(jié)合。因為其具有硬度高、耐磨性好、彈性模量高以及摩擦因數(shù)低等與金剛石類似的性能，很多場合代替金剛石使用以降低生產(chǎn)成本；類金剛石涂層的諸多優(yōu)點使其受到廣泛關(guān)注[6-7]。目前，國內(nèi)外學(xué)者對于DLC涂層耐沖擊性能的研究主要分為實驗探究和數(shù)值仿真模擬2種方式，且大多集中在實驗方面，施紹淼[8]研究了不同結(jié)構(gòu)類金剛石碳薄膜的動態(tài)沖擊摩擦和靜態(tài)疲勞載荷磨損性能和失效機理。王璋等[9]85通過動力學(xué)響應(yīng)考察Cr-DLC涂層的沖擊磨損行為，用沖擊力峰值和能量吸收率來表征抗沖擊性的優(yōu)劣。占智遠[10]設(shè)計并制備了類金剛石薄膜碳納米管陣列復(fù)合結(jié)構(gòu)，以能量耗散來表征其對基體抗沖擊性能的提升。對于DLC涂層應(yīng)用于紡織鋼領(lǐng)的數(shù)值模擬研究相對較少。陳崇亮[11]建立了連續(xù)沖擊情況下DLC涂層與基體剝離的物理模型，認為類金剛石涂層在沖擊載荷下極易產(chǎn)生疲勞，形成微小裂紋，進而發(fā)生涂層剝落，但其只是運用理論推理得出，并沒有進行相應(yīng)的仿真模擬。KANG等[12]將有限元分析模型與納米壓痕技術(shù)相結(jié)合的方法來確定多層FG Ti/TiC/DLC涂層材料的屈服強度,但其模型建立在納米壓痕技術(shù)之上，并不能代表沖擊時的動態(tài)響應(yīng)。課題組以DLC應(yīng)用于紡織鋼領(lǐng)涂層為背景，研究DLC涂層基體和純基體在不同沖擊荷載作用下的力學(xué)響應(yīng)和破壞形態(tài)。針對DLC涂層的特性，運用有限元軟件ABAQUS建立等效沖擊模型，通過對比2種鋼領(lǐng)基體材料受到不同沖擊載荷下的結(jié)構(gòu)瞬時動力響應(yīng)，考察可能存在的涂層破壞形態(tài)，為DLC涂層在鋼領(lǐng)領(lǐng)域的應(yīng)用提供一定的參考；為本質(zhì)上改善鋼領(lǐng)的服役壽命提供技術(shù)支撐，最終實現(xiàn)鋼領(lǐng)的長效運行。

1 鋼領(lǐng)/鋼絲圈沖擊有限元模型

1.1 鋼領(lǐng)/鋼絲圈結(jié)構(gòu)幾何模型

隨著數(shù)值仿真模擬技術(shù)的不斷完善，仿真模擬實驗已經(jīng)在很多領(lǐng)域得到應(yīng)用。仿真模擬實驗與實際實驗相比，具有不受空間和地點限制，試錯成本低，節(jié)約時間成本，以及提取相關(guān)數(shù)據(jù)極為方便的特點。同時，仿真模擬可以看到試樣內(nèi)部的力學(xué)響應(yīng)以及不易發(fā)現(xiàn)的損傷變形，從而不斷完善設(shè)計方案，在實際實驗條件比較極端或者不容易實現(xiàn)時，仿真模擬的優(yōu)勢會更加明顯。目前，大多數(shù)紡織機械的鋼領(lǐng)與鋼絲圈的相對線速度在18～40 m/s之間,有的甚至在40 m/s以上，鋼領(lǐng)與鋼絲圈之間的沖擊距離很短，但在如此高的線速度下，鋼領(lǐng)受到的沖擊依然不可小覷。

對環(huán)錠細紗機鋼領(lǐng)/鋼絲圈結(jié)構(gòu)進行簡化并構(gòu)建有限元沖擊模型。運用SoildWorks三維建模軟件建立簡化沖擊模型,模型主要包括圓球沖擊體、基體和涂層3部分。圓球沖擊體模擬鋼絲圈與鋼領(lǐng)接觸的尖端部分，將其視為圓形小球，半徑3 mm；涂層材料為類金剛石(DLC)，涂層厚度300 μm；基體為厚度5 mm，直徑?30 mm的圓盤。圓球沖擊方向與基體上表面呈90°，且過基體表面圓心軸線，2物體之間設(shè)置一微小間隔，以模擬實際沖擊距離，更好地觀察沖擊過程。將其導(dǎo)入ABAQUS有限元分析軟件，給予基體四周固定約束(U3=U2=U1=UR1=UR2=UR3)，即分別限制基體在X,Y,Z方向的移動和轉(zhuǎn)動；沖擊體只允許z向平移(U1=U2=UR3=0)，即限制沖擊體在X,Y方向的移動和Z方向的轉(zhuǎn)動。假定涂層與基體理想結(jié)合在一起(即完全結(jié)合)，施加綁定約束(Tie)，同時為防止沖擊引起的體接觸使2體之間發(fā)生網(wǎng)格竄動，在沖擊體表面與涂層上表面之間設(shè)置切向無摩擦接觸(tangential behavior frictionlness)，涂層上表面和基體間采用罰函數(shù)接觸法(penalty contact method)。根據(jù)鋼領(lǐng)實際運行速度，對沖擊體(圓球)施加初始動能即初速度，運用動能與速度關(guān)系公式ΔEk=1/2[m(v2-v′2)]，可實現(xiàn)動能與速度間的轉(zhuǎn)換；ΔEk為動能的變化量，m為物體的質(zhì)量，v和v′分別為沖擊運動前后的速度。圓球沖擊體初速度v選取10，20，30和40 m/s(意在確定圓球的初始動能，忽慮重力等影響因素)。在滿足計算結(jié)果精度的同時兼顧運算效率，對涂層和基體采用不同的網(wǎng)格密度，將基體、涂層和圓球進行網(wǎng)格分區(qū)?；w和涂層同時分為中心細化區(qū)域和外圍區(qū)域，中心細化區(qū)域取邊長為4 mm的正方形，劃分為網(wǎng)格較密的單元，外圍區(qū)域受到的球體沖擊作用較小，可視為無限域，故劃分為相對稀疏的單元；在基體上表面以及其與涂層交界處劃分較密網(wǎng)格，而基體下表面劃分為較稀疏網(wǎng)格，上密下疏，逐漸過渡。單元類型采用C3D8T，鋼領(lǐng)/鋼絲圈及仿真沖擊模型網(wǎng)格劃分如圖1所示。

圖1 鋼領(lǐng)/鋼絲圈及仿真沖擊模型網(wǎng)格劃分Figure 1 Grid of textile ring/ring traveler and simulation impact model

1.2 材料選取及力學(xué)性能

假定DLC涂層為理想彈塑性材料(elastic-perfectly plastic materia)[13]，圓球沖擊體施加剛體模型(rigid body)，其不作為考察對象且對結(jié)果的影響可以忽略，基體采用20#鋼材料常用的Johnson-Cook本構(gòu)模型[14]，表1給出具體材料力學(xué)屬性。

1.3 沖擊模型驗證

通過逐步細化基體網(wǎng)格密度探究仿真結(jié)果的穩(wěn)定性。圓球采用剛體模型，其網(wǎng)格密度隨基體網(wǎng)格密度逐步細化，對仿真結(jié)果的影響可以忽略；將基體中心細化部分邊長劃分單元個數(shù)作為變量，初始沖擊速度選取30 m/s ，取基體所受沖擊載荷力的最大值和沖擊下陷深度的最大值來探究網(wǎng)格密度對于仿真實驗結(jié)果的影響程度。擬合數(shù)據(jù)得到網(wǎng)格無關(guān)性仿真曲線圖，如圖2所示。

表1 沖擊模型材料力學(xué)屬性

圖2 網(wǎng)格無關(guān)性仿真曲線Figure 2 Grid-independent simulation graph

從圖2曲線可以看出2組數(shù)據(jù)均經(jīng)歷一定波動后趨于穩(wěn)定狀態(tài)；當中心細化區(qū)域邊長單元個數(shù)達到27后波動逐漸減小，在單元個數(shù)達到40及其以上時數(shù)據(jù)波動已控制在±0.4%以內(nèi)。由此得出：當中心細化區(qū)域邊長單元個數(shù)達到40以上時，盲目提高網(wǎng)格數(shù)量，只會單方面提高計算機的運算負擔，但對最終結(jié)果的影響可以忽略。本有限元模型選取中心細化區(qū)域邊長單元個數(shù)為62，節(jié)點總數(shù)845 711，單元總數(shù)814 680。

使用此有限元模型計算初始沖擊能(initial impact energy)分別為0.08，0.17，0.33和0.50 mJ時基體所受到的載荷力，將其與王璋等所做的實驗結(jié)果[9]80進行比較，擬合曲線如圖3所示。2條曲線皆隨著初始沖擊能的增加而不斷升高，且仿真模擬結(jié)果與實際沖擊實驗結(jié)果非常接近，偏差在5%以內(nèi)；仿真模擬值偏高是由于在實際實驗中空氣阻力或人為干擾等因素導(dǎo)致的能量流失 ，但整體上仿真模擬結(jié)果與實際實驗結(jié)果趨勢一致，由此驗證了此沖擊模型的正確性和仿真結(jié)果的有效性。

圖3 仿真計算與實驗結(jié)果對比Figure 3 Comparison between simulation calculation and experimental results

2 不同沖擊載荷下模型動態(tài)響應(yīng)分析

2.1 基體殘余應(yīng)力及沖擊坑深度分析

通過對比各沖擊速度下純基體和DLC涂層基體的表面和內(nèi)部殘余應(yīng)力分布及沖擊坑深度，來考察DLC涂層的抗沖擊性能以及對于鋼領(lǐng)基體的保護作用。2種基體的殘余應(yīng)力云圖如圖4所示，基體表面沖擊坑呈環(huán)狀分布，隨著沖擊速度的提升，沖擊坑直徑不斷加大，且在沖擊坑底部和邊緣部分存在應(yīng)力集中現(xiàn)象。純基體的沖擊坑相對較大，且殘余應(yīng)力值在表面分布相對均勻；DLC涂層基體的殘余應(yīng)力由沖擊坑底部向外圍逐步擴散，且沖擊坑底部為殘余拉應(yīng)力，沖擊坑四周內(nèi)壁為殘余壓應(yīng)力。

圖4 基體殘余應(yīng)力分布云圖Figure 4 Cloud map of matrix residual stress distribution

取沖擊速度為30 m/s時2種基體殘余應(yīng)力場截面分布云圖，如圖5所示。2種基體內(nèi)部的殘余應(yīng)力場分布有較大差異。純基體受圓球沖擊后殘余應(yīng)力主要分為2部分：基體表面受沖擊形成的沖擊坑隆起部分和沖擊坑下方區(qū)域。沖擊坑隆起部分的殘余應(yīng)力既有拉應(yīng)力也有壓應(yīng)力，最大值分別為379.2 MPa和-361.7 MPa；沖擊坑下方內(nèi)部區(qū)域分布著一部分殘余壓應(yīng)力，最大值為-164.1 MPa，距離沖擊坑底部1.69 mm。沖擊坑底部區(qū)域的應(yīng)力梯度變化較為緩和，但在沖擊坑隆起部分區(qū)域，應(yīng)力值變化突然增大，基體材料受沖擊發(fā)生塑性變形使其在沖擊的下方區(qū)域出現(xiàn)一定的殘余壓應(yīng)力。由于基體表面的形貌變化，一部分基體材料受沖擊變形堆積于沖擊坑邊緣處，同時圓球向下(z方向)沖擊時，也伴隨著其向沖擊坑四周(x-y方向)的擠壓作用，基體材料受到拉力和擠壓的共同作用，當沖擊結(jié)束離開后，堆積較多材料處就形成了較高的殘余應(yīng)力。

涂覆DLC涂層后，殘余拉應(yīng)力主要集中在沖擊坑底部，呈橢圓狀向四周輻射，其最大值為761.1 MPa，最大壓應(yīng)力出現(xiàn)在基體上表面沖擊坑周圍區(qū)域，最大值為-402 MPa?；w上表面沖擊坑周圍區(qū)域的應(yīng)力梯度變化較大，但整體的應(yīng)力分布相對平緩。 基體沖擊坑周圍區(qū)域同樣受材料堆積影響而形成殘余壓應(yīng)力， 而在沖擊坑中心，基體同時受到涂層擠壓、材料堆積和涂層與其之間的結(jié)合力影響，導(dǎo)致其產(chǎn)生較高的殘余應(yīng)力。

在不同沖擊載荷下,涂覆DLC涂層的基體沖擊坑深度明顯降低，且從圖6基體沖擊坑深度截面云圖可以看出，沖擊坑邊緣的基體材料隆起、材料堆積和應(yīng)力集中現(xiàn)象也明顯改善。隨著沖擊動能的增加，2種基體的沖擊坑深度均不斷增加，但DLC涂層基體的沖擊坑深度總體小于純基體，且隨著初始動能的上升，差值趨于下降，這反映了DLC涂層對于沖擊坑深度的減緩效果，且隨著沖擊動能升高，減緩效果在降低。由此可知在鋼領(lǐng)表面涂覆DLC涂層可有效地改善表面形貌變化和基體沖擊坑外圍的應(yīng)力集中現(xiàn)象，提升鋼領(lǐng)的抗沖擊性能；當沖擊動能不斷升高，涂層對于基體的保護作用在降低，當沖擊坑底部的應(yīng)力值超過涂層與基體的粘結(jié)強度后，會引起涂層與基體分離形成裂紋，同時在沖擊坑隆起部分因拉伸應(yīng)力作用產(chǎn)生環(huán)形裂紋，出現(xiàn)涂層剝落現(xiàn)象，這對于涂層與基體間的粘結(jié)強度來說是很大的考驗。

圖6 基體沖擊坑深度截面云圖及時程曲線Figure 6 Matrix impact crater depth section cloud map and time curve

2.2 應(yīng)變與載荷力分析

通過對比各沖擊速度下2種基體的截面應(yīng)變云圖以及時程曲線，分析基體內(nèi)部應(yīng)變值與其所受載荷力的關(guān)系，進而探究DLC涂層的抗沖擊性能機理。由圖7時程曲線可以看出2種基體的最終應(yīng)變值均為負數(shù)值，說明了沖擊造成的主要是壓應(yīng)變，且2種基體皆產(chǎn)生不同程度的塑性變形。由圖7基體截面應(yīng)變云圖可以看出純基體和DLC涂層基體的應(yīng)變分布基本相同，壓應(yīng)變主要在沖擊坑底部，拉應(yīng)變集中于沖擊坑邊緣隆起部分；2種基體拉應(yīng)變最大值分別為0.142 90和0.008 06，壓應(yīng)變最大值分別為0.387 4和0.055 6。DLC涂層基體的應(yīng)變值遠小于純基體，且在沖擊過后恢復(fù)一定應(yīng)變值，而純基體恢復(fù)的量很小，這反映了DLC涂層不僅分擔了一部分的應(yīng)變，還使基體受到的應(yīng)變進一步降低。

圖7 基體截面應(yīng)變云圖及時程曲線Figure 7 Strain cloud map and time curve of matrix section

由圖8載荷力時程曲線圖分析可知2種試樣的載荷力時程曲線都經(jīng)歷了先升高再降低為0的過程。DLC涂層基體在4種沖擊速度下的載荷力最大值均大于純基體，這與其沖擊坑底部的應(yīng)力值較大相符合。DLC涂層延長了沖擊接觸時間，使應(yīng)變發(fā)生相對緩和，涂層的彈性模量大于基體，這使得涂覆涂層的基體不易變形，而載荷力受接觸面積和自身動能影響，當動能相同時，DLC涂層的變形較小，意味著圓球與涂層的接觸面積較小，其瞬時載荷力也就相應(yīng)變大；這也同時解釋了基體沖擊坑底部高應(yīng)力集中現(xiàn)象，因此選用鋼領(lǐng)涂層材料時，在其他材料參數(shù)一定的情況下，應(yīng)盡量提高涂層的彈性模量。

圖8 載荷力時程曲線Figure 8 Load force time history curve

2.3 圓球動能及基體總能對比

取圓球的動能變化和2種基體所受總熱力學(xué)能數(shù)據(jù)匯總得到圖9和圖10?；w總熱力學(xué)能包括可恢復(fù)彈性應(yīng)變能、非彈性過程的能量耗散(主要指塑性應(yīng)變)、黏彈性或者蠕變過程的能量耗散和偽應(yīng)變能。仿真沖擊過程時間間隔較短，其中不涉及固液之間的轉(zhuǎn)換和應(yīng)力的長時間加載條件， 故此處不考慮黏彈性或者蠕變過程的能量耗散和偽應(yīng)變能，只探究圓球的全過程動能變化和基體所受總熱力學(xué)能的全時程變化情況。該處理方式可以很好地反映整個沖擊過程中圓盤和基體之間的能量轉(zhuǎn)化及DLC涂層在沖擊過程中對基體材料起到的保護作用。從圓球的全過程動能曲線圖中可以看出有DLC涂層的基體和純基體在受到圓球各初始動能(初速度)沖擊時，都表現(xiàn)出動能先下降至0后又上升的趨勢，這符合實際實驗中圓球沖擊基體然后再反彈的現(xiàn)象，同時DLC涂層基體在受到圓球同等初速度沖擊后所回升的動能值要普遍大于純基體。圓球?qū)幽軅鬟f給涂層表面，再通過涂層傳導(dǎo)給基體，將動能轉(zhuǎn)化為材料的熱力學(xué)能儲存在材料內(nèi)部，而后其中一部分熱力學(xué)能重新轉(zhuǎn)化為圓球的動能；純基體在沖擊過程中發(fā)生較多塑性變形，吸收的能量較大，導(dǎo)致圓球沖擊后動能也相對較低。而隨著沖擊動能的增加，2種試樣的總熱力學(xué)能差值在逐漸縮小，故沖擊動能的增加會使DLC涂層對基體的保護作用減弱，王璋等[9]83所做的涂層動力學(xué)響應(yīng)實驗印證了此觀點。涂層和基體內(nèi)部儲存的熱力學(xué)能會越積越多，當達到裂紋萌生或擴展到臨界熱力學(xué)能時，就會引起涂層裂紋的產(chǎn)生或者不斷擴展，最終導(dǎo)致涂層失效。

圖9 圓球動能曲線Figure 9 Sphere kinetic energy curve

圖10 總熱力學(xué)能曲線Figure 10 Matrix total internal energy curve

3 結(jié)論

課題組利用ABAQUS有限元分析軟件建立一種針對環(huán)錠細紗機紡織鋼領(lǐng)系統(tǒng)簡化沖擊模型，對DLC涂層基體和純基體進行仿真沖擊模擬實驗，通過改變沖擊體初速度的方式來滿足不同沖擊載荷的要求，用對比基體殘余應(yīng)力應(yīng)變、沖擊坑深度、最大載荷力以及系統(tǒng)能量變化的方式，來分析DLC涂層的抗沖擊性能，同時考察DLC涂層在實際工況下經(jīng)受沖擊后其表面內(nèi)部可能存在的一些應(yīng)力集中區(qū)域和涂層破壞形式，得到以下結(jié)論：

1) DLC涂層可以有效降低沖擊坑的直徑和深度，并改善鋼領(lǐng)的表面形貌和殘余應(yīng)力集中現(xiàn)象，在沖擊坑底部和邊緣隆起部分有較多殘余應(yīng)力集中，會引起涂層生成裂紋，導(dǎo)致涂層剝落失效，應(yīng)盡量選用較高粘結(jié)強度的DLC涂層；

2) DLC涂層可以顯著降低基體應(yīng)變值，并延長沖擊接觸時間，較高的彈性模量使基體應(yīng)變降低，從而提升涂層的抗沖擊性能，在選用DLC涂層材料時，在其他材料參數(shù)一定的情況下，應(yīng)盡量提高涂層的彈性模量；

3) 從系統(tǒng)能量的角度考察了DLC涂層的抗沖擊性能，涂覆DLC涂層后，基體產(chǎn)生更少的塑性變形，能量耗散也相對變小，有更多的能量轉(zhuǎn)化為基體的熱力學(xué)能和圓球的動能，基體總熱力學(xué)能和圓球動能在沖擊結(jié)束后的回復(fù)量顯著升高；

4) 選用DLC涂層材料時，較高的涂層彈性模量以及較大的膜基粘結(jié)強度對于鋼領(lǐng)整體抗沖擊性能的提升至關(guān)重要。