初紅艷 陳立博 安 然 蔡力鋼

北京工業(yè)大學(xué)先進(jìn)制造與智能技術(shù)研究所，北京，100124

0 引言

隨著社會(huì)的發(fā)展，人們對(duì)印刷品的質(zhì)量要求也越來(lái)越高。為了提高印刷質(zhì)量，國(guó)內(nèi)外學(xué)者進(jìn)行了諸多研究。BHOOMKAR等[1]提出用“世界一流制造”(WCM)的方法來(lái)提高膠印機(jī)的印刷質(zhì)量，采用油墨溫度控制裝置保持油墨溫度恒定以保證印刷質(zhì)量；ROBERT等[2]研究了現(xiàn)代柔性印刷機(jī)中溫度對(duì)油墨黏性及油墨轉(zhuǎn)移率的影響。目前，國(guó)內(nèi)外學(xué)者對(duì)普通膠印機(jī)墨輥的熱效應(yīng)研究較少，而隨著技術(shù)的發(fā)展，膠印機(jī)的印刷速度越來(lái)越高，在印刷過(guò)程中，隨著印刷速度的提高，油墨、墨輥的溫度也隨之升高。溫度能顯著地改變油墨的黏度，進(jìn)而影響油墨在墨輥間的流動(dòng)與傳遞，因此對(duì)印刷過(guò)程中油墨的溫升情況進(jìn)行研究很有必要。

彈性流體動(dòng)力潤(rùn)滑簡(jiǎn)稱(chēng)彈流潤(rùn)滑，它主要研究點(diǎn)線(xiàn)接觸中彈性表面間的流體動(dòng)力潤(rùn)滑問(wèn)題，如接觸區(qū)流體壓力與流體通道厚度等。膠印機(jī)輸墨系統(tǒng)中，在正常工況下，一軟一硬兩墨輥相互擠壓對(duì)滾，油墨從墨輥間的墨路通道通過(guò)，軟輥在擠壓中發(fā)生彈性變形，且接觸區(qū)相對(duì)較小，該模型可近似地看成一個(gè)彈流潤(rùn)滑系統(tǒng)，對(duì)接觸區(qū)的溫度分析需要用到熱彈流潤(rùn)滑理論。

學(xué)者針對(duì)大剛度接觸表面間的熱彈流潤(rùn)滑進(jìn)行了研究。王大偉等[3]在綜合考慮載荷、速度及潤(rùn)滑劑熱效應(yīng)等影響因素的基礎(chǔ)上，基于Ree-Eyring型非牛頓流體模型建立了滾動(dòng)軸承熱彈流潤(rùn)滑的數(shù)學(xué)模型，求得了壓力、溫度及油膜厚度的分布，探討了速度、載荷、滾子半徑等參數(shù)對(duì)滾動(dòng)軸承潤(rùn)滑特性的影響；BRACCESI等[4]應(yīng)用熱彈流潤(rùn)滑相關(guān)理論分析了滾珠絲杠中的潤(rùn)滑油溫度變化；武斌[5]建立了低速重載齒輪傳動(dòng)微觀熱彈流潤(rùn)滑的數(shù)學(xué)模型，并針對(duì)該數(shù)學(xué)模型中形式相當(dāng)復(fù)雜的二階非線(xiàn)性偏微分方程組的特征，運(yùn)用多重網(wǎng)格法求解了接觸區(qū)油膜的溫升情況。

本文結(jié)合印刷機(jī)實(shí)際情況，對(duì)一軟一硬兩個(gè)接觸表面之間的熱彈流潤(rùn)滑問(wèn)題進(jìn)行詳細(xì)研究。

1 熱彈流潤(rùn)滑理論

熱彈流潤(rùn)滑理論的研究對(duì)象是點(diǎn)線(xiàn)接觸摩擦副中的熱現(xiàn)象和熱效應(yīng)，包括接觸區(qū)內(nèi)油膜溫度分布及溫度場(chǎng)對(duì)壓力分布、油膜形狀與厚度等的影響。接觸區(qū)內(nèi)溫度升高會(huì)使?jié)櫥瑒┑奈锢砘瘜W(xué)特性發(fā)生較大的變化，因此對(duì)接觸區(qū)內(nèi)溫度場(chǎng)的研究意義重大。等溫彈流潤(rùn)滑理論中忽略了熱效應(yīng)的影響，對(duì)于速度較小或滑動(dòng)速度較小時(shí)的情況較為適用。當(dāng)高速滾動(dòng)或有較大滑動(dòng)時(shí)，就必須考慮熱效應(yīng)所帶來(lái)的影響。在求解熱彈流潤(rùn)滑問(wèn)題時(shí)，需在等溫彈流潤(rùn)滑的基礎(chǔ)上聯(lián)立能量方程、熱界面方程，即需聯(lián)立求解雷諾方程、彈性變形方程、能量方程、熱界面方程、運(yùn)動(dòng)方程。

1.1 雷諾方程

雷諾方程可描述狹小間隙中黏性流體的流動(dòng)，是流體力學(xué)中Navier-Stokes方程的一種特殊形式。在線(xiàn)彈流潤(rùn)滑情況下，當(dāng)潤(rùn)滑油在兩固體表面之間進(jìn)行層流流動(dòng)時(shí)，所形成的油膜壓力分布可由下式來(lái)描述：

(1)

式中，u為卷吸速度；η為黏度；p為油膜壓力；ρ為潤(rùn)滑油密度；t為時(shí)間；h為膜厚；x方向?yàn)闈?rùn)滑油流動(dòng)方向。

式(1)為線(xiàn)接觸條件下的雷諾方程，適用于求解各種非穩(wěn)態(tài)工況下流體潤(rùn)滑問(wèn)題。

1.2 彈性變形方程

線(xiàn)接觸問(wèn)題通?？梢暈闊o(wú)限長(zhǎng)柱體的接觸。在膠印機(jī)的實(shí)際工況下，兩墨輥接觸區(qū)相對(duì)于墨輥長(zhǎng)度特別小，因此適用于線(xiàn)接觸問(wèn)題。在無(wú)限平面上的任意分布力p(x)在M點(diǎn)產(chǎn)生的彈性變形δ(x)見(jiàn)圖1。

圖1 分布力作用下的變形Fig.1 Deformation under the distributed force

在彈流潤(rùn)滑中，兩個(gè)被潤(rùn)滑油膜隔開(kāi)的表面都是彈性表面。由于受到油膜壓力的作用，兩表面受到數(shù)值相等而方向相反的法向分布力的作用，從而產(chǎn)生了彈性變形，將彈性變形與兩柱體間隙厚度相疊加，就可得到線(xiàn)接觸彈流問(wèn)題的油膜幾何方程(即彈性變形方程)：

(2)

式中，hc為中心膜厚，此時(shí)hc即為原點(diǎn)處的實(shí)際膜厚；E′為當(dāng)量彈性模量；R為當(dāng)量曲率半徑。

1.3 能量方程

在熱彈流潤(rùn)滑理論中，當(dāng)發(fā)生高速滾動(dòng)或滾動(dòng)中伴有較大滑動(dòng)速度時(shí)，就必須考慮潤(rùn)滑中的熱效應(yīng)和溫度場(chǎng)的影響。在潤(rùn)滑過(guò)程中，油膜受到黏性剪切和壓縮作用而發(fā)熱，造成溫度升高，同時(shí)，所產(chǎn)生的熱量通過(guò)對(duì)流和熱傳導(dǎo)而散失，這樣，就需要能量方程來(lái)計(jì)算潤(rùn)滑油膜的溫度場(chǎng)[6]。

如果不考慮體積力和熱輻射的影響，根據(jù)能量守恒定律，黏性流體流動(dòng)時(shí)的能量方程有如下普遍形式：

(3)

(4)

式中 ，T為溫度;cp為質(zhì)量定壓熱容；Φ為耗散功；k為熱導(dǎo)率；u、v、w分別為直角坐標(biāo)系下流體沿x、y、z方向的速度；y為油膜流動(dòng)法向；z為垂直于xy平面的方向。

(5)

1.4 熱界面方程

在求解潤(rùn)滑膜的能量方程時(shí)，需要根據(jù)固體的熱傳導(dǎo)情況來(lái)確定潤(rùn)滑膜與固體相接觸的界面上的邊界條件。

一般形式的固體熱傳導(dǎo)方程為

(6)

式中，c為質(zhì)量熱容。

由于接觸體處于運(yùn)動(dòng)狀態(tài)，故接觸表面的溫度計(jì)算可以歸結(jié)為具有移動(dòng)熱源的半無(wú)限大體的熱傳導(dǎo)問(wèn)題。忽略x和y方向的熱傳導(dǎo)，最終可得兩個(gè)界面上邊界條件：

(7)

式中，K1、K2分別為兩個(gè)界面的傳熱系數(shù)；T0為初始溫度。

1.5 運(yùn)動(dòng)方程

在求解能量方程時(shí)，需要先確定潤(rùn)滑油膜的速度場(chǎng)，這時(shí)就需要根據(jù)給定的壓力分布和油膜厚度來(lái)求解運(yùn)動(dòng)方程以確定速度場(chǎng)。

對(duì)于線(xiàn)接觸問(wèn)題，運(yùn)動(dòng)方程為

(8)

2 油墨溫度場(chǎng)的熱彈流潤(rùn)滑模型建立

2.1 方程的量綱一化

2.2 方程所需參數(shù)的確定

為了給下文油墨溫度分析提供數(shù)據(jù)支持及分析載荷、速度、材料等參數(shù)變化對(duì)油墨溫度的影響趨勢(shì)，選取了包含實(shí)際工況參數(shù)在內(nèi)的較大范圍的參數(shù)。

在調(diào)節(jié)印刷壓力時(shí)，往往以?xún)奢仈D壓時(shí)的接觸寬度作為判斷依據(jù)，北人集團(tuán)公司生產(chǎn)的某對(duì)開(kāi)膠印機(jī)的串墨輥的推薦接觸半寬為4～5 mm，根據(jù)有限元軟件仿真模擬反推可知載荷力的大小約為100 N，本文取100～500 N的載荷范圍。

在實(shí)際印刷過(guò)程中，正常的印刷速度為0.36～1.5萬(wàn)張每小時(shí)，對(duì)應(yīng)的最小墨輥角速度約為10 rad/s，本文取10～30 rad/s的速度范圍。

依據(jù)墨輥實(shí)際情況確定各項(xiàng)參數(shù)，兩輥半徑均為32.5 mm，上下兩輥相互擠壓對(duì)滾，見(jiàn)圖2。

圖2 擠壓對(duì)滾模型Fig.2 Squeezing rolling model

實(shí)際工況中，膠印機(jī)的硬輥和軟輥均采用鋼制輥芯，但鋼芯直徑一般為墨輥直徑的一半，本文重點(diǎn)分析墨輥擠壓接觸區(qū)的油墨溫度，接觸區(qū)距離鋼制輥芯較遠(yuǎn)，所以忽略鋼質(zhì)輥芯的影響。設(shè)定硬輥(硬塑膠輥)彈性模量E1為4 GPa，泊松比μ1為0.33；軟輥(橡膠輥)彈性模量E2為7.84 MPa，泊松比μ2為0.4。印刷油墨的初始環(huán)境黏度為42 Pa·s；印刷油墨的環(huán)境密度為1.5 g/cm3；通過(guò)實(shí)驗(yàn)測(cè)得，油墨黏壓系數(shù)為2.57×10-8m2/N。橡膠的熱導(dǎo)率為0.2 W/(m·K)，質(zhì)量定壓熱容為1 600 J/(kg·K)，密度為1 050 kg/m3。塑料的熱導(dǎo)率為0.12 W/(m·K)，質(zhì)量定壓熱容為1 540 J/(kg·K)，密度為1 050 kg/m3。

2.3 方程的求解

與等溫彈流潤(rùn)滑相比，熱彈流潤(rùn)滑的基本問(wèn)題是要根據(jù)能量方程求解潤(rùn)滑膜的溫度分布，同時(shí)，在求解雷諾方程時(shí)，將黏度作為溫度和壓力的函數(shù)，而由能量方程可知，溫度場(chǎng)的計(jì)算又取決于壓力的分布和界面邊界條件的確定，這樣，在熱彈流潤(rùn)滑的計(jì)算中，除去等溫彈流潤(rùn)滑所需的方程外，還需要聯(lián)立求解能量方程和熱界面方程。

由于熱彈流潤(rùn)滑數(shù)值計(jì)算的工作量巨大，故侯克平等[8]提出了熱彈流潤(rùn)滑的簡(jiǎn)化求解方法，其基本思想是忽略熱效應(yīng)對(duì)壓力分布和膜厚的影響，在等溫條件下求得不同工況下彈流潤(rùn)滑壓力分布和油膜厚度，然后在此基礎(chǔ)上，直接通過(guò)求解能量方程、運(yùn)動(dòng)方程和熱界面方程及潤(rùn)滑油的黏溫方程來(lái)確定油膜的溫度分布。

在計(jì)算時(shí)，由于能量方程中的溫度、黏壓溫方程中的黏度及彈性變形量都隨著壓力變化，因此，首先應(yīng)給出一個(gè)初始的壓力分布，本文采用ANSYS軟件計(jì)算所得的接觸應(yīng)力作為初始?jí)毫Ψ植?；其次，?yīng)給出一初始溫度分布，本文采用數(shù)值等于室溫的均勻溫度場(chǎng)；計(jì)算黏度值后，代入能量方程求得溫度，利用新的溫度值修正黏度，反復(fù)迭代直至兩次計(jì)算的溫度差小于0.01 ℃，即可看作溫度收斂。求解過(guò)程見(jiàn)圖3。

圖3 溫度求解流程Fig.3 Temperature solving process

3 油墨層溫度分析

3.1 接觸區(qū)油墨層壓力及厚度分布

代入載荷100 N進(jìn)行運(yùn)算，接觸區(qū)油墨壓力分布及油墨墨層厚度分布見(jiàn)圖4。圖4中，橫坐標(biāo)為接觸區(qū)節(jié)點(diǎn)編號(hào)，接觸區(qū)中心處在橫坐標(biāo)中70號(hào)節(jié)點(diǎn)所在位置，可以看出，油墨墨層壓力最大值出現(xiàn)在接觸區(qū)中心稍靠近出口處，最大值約為1.44 MPa。油墨墨層厚度在接觸區(qū)前逐漸減小，在接觸區(qū)內(nèi)相對(duì)穩(wěn)定，在近出口處出現(xiàn)了頸縮的現(xiàn)象，頸縮處油墨厚度最小，接觸區(qū)中心墨層厚度約為0.073 mm。

(a)壓力分布 (b)厚度分布圖4 壓力及厚度分布Fig.4 The distribution of pressure and thickness

3.2 接觸區(qū)油墨層溫度分布

沿油墨層厚度方向取5層，所得油墨各層溫度分布見(jiàn)圖5。

圖5 墨層厚度方向各層溫度Fig.5 Layers temperature

當(dāng)只討論油墨中間層溫度時(shí)，可以看出隨著油墨逐漸進(jìn)入擠壓區(qū)，油墨的溫度不斷升高，在接觸區(qū)墨層壓力較大的區(qū)域，油墨的溫度變化也較大，最后隨著遠(yuǎn)離最大壓力區(qū)，油墨溫度逐漸降低。這是由于高壓力使得油墨黏度增加，從而能量方程中的耗散功一項(xiàng)增大，使得高壓力區(qū)油墨溫度的值也較高。

比較不同位置的溫升，可以發(fā)現(xiàn)，中間層油墨溫度最高，100 N時(shí)，中間層最大溫升值為2.231 ℃，靠近兩墨輥處的溫度較低，塑料表面層最大溫升值為1.851 ℃，橡膠表面層最大溫升值為1.679 ℃。這是由于換熱的影響，油墨所產(chǎn)生的熱量傳遞給了兩側(cè)墨輥，從而導(dǎo)致了靠近墨輥處的油墨溫度降低。比較橡膠表面層和塑料表面層的墨層溫度，可以看出塑料表面層的溫度更高，這是由于塑料的熱導(dǎo)率較小，因此傳播熱的能力較低，這樣，塑料表面就容易聚集較多的熱量，從而產(chǎn)生較高的溫升。

3.3 載荷對(duì)接觸區(qū)溫升的影響

保持其他參數(shù)不變，只改變載荷力大小，計(jì)算接觸區(qū)最大溫升，設(shè)定載荷分別為100 N、150 N、200 N、250 N、300 N，油墨中間層最大溫升見(jiàn)圖6。

圖6 不同載荷下的溫升Fig.6 Temperature rise under different loads

載荷增加即接觸區(qū)的壓力增大，隨著壓力的增大,油墨黏度增大，同時(shí)能量方程中的耗散功一項(xiàng)增大，因此，隨著載荷的增大，油墨中間層最大溫升提高，溫升由載荷為100 N時(shí)的2.231 ℃，升高到300 N時(shí)的3.41 ℃ 。

3.4 速度對(duì)接觸區(qū)溫升的影響

保持其他參數(shù)不變，只改變兩輥對(duì)滾速度大小，計(jì)算接觸區(qū)最大溫升，設(shè)定對(duì)滾速度分別為10 rad/s、15 rad/s 、20 rad/s、25 rad/s、30 rad/s，接觸區(qū)油墨中間層最大溫升見(jiàn)圖7。

圖7 不同速度下的溫升Fig.7 Temperature rise under different speeds

由圖7可以看出，隨著兩輥對(duì)滾速度的增加，兩輥間接觸區(qū)油墨中間層的溫度變化不大，即對(duì)滾速度對(duì)溫升的影響較小。

3.5 材料對(duì)接觸區(qū)溫升的影響

由于軟輥的材質(zhì)多為橡膠，且不同橡膠的彈性模量等參數(shù)相差較大，為了分析不同材質(zhì)下的油墨層的最大溫升，本節(jié)分別選用彈性模量為8 MPa、16 MPa、24 MPa、32 MPa的橡膠進(jìn)行溫度場(chǎng)分析。

保持其他參數(shù)不變，只改變橡膠材質(zhì)，接觸區(qū)油墨中間層最大溫升見(jiàn)圖8?？梢钥闯?，隨著軟輥彈性模量的增大，油墨中間層的最大溫升提高，由8 MPa時(shí)的2.23 ℃升高到32 MPa時(shí)的2.94 ℃。這是由于軟輥彈性模量的增大導(dǎo)致了當(dāng)量彈性模量的增大，當(dāng)量彈性模量的增大導(dǎo)致了接觸區(qū)壓力分布的變化，壓力的增大使得油墨層溫度提高。

圖8 不同軟輥材料下的溫升Fig.8 Temperature rise under different materials

4 結(jié)論

(1)比較沿厚度方向不同位置處的墨層溫度變化，結(jié)果顯示中間層油墨溫度最高，越靠近兩墨輥處溫度越低。這是由于靠近兩輥表面的油墨受換熱的影響，流失了一部分熱量，且熱導(dǎo)率低的材料表面溫度較高。

(2)隨著載荷的增加，油墨中間層最大溫升提高。這是由于載荷增加即接觸區(qū)的壓力增大，隨著壓力的增大,油墨黏度增大，同時(shí)能量方程中的耗散功一項(xiàng)增大，進(jìn)而使接觸區(qū)最大溫升升高。

(3)隨著兩輥對(duì)滾速度的增加，兩輥間接觸區(qū)油墨中間層的溫度變化不大，即對(duì)滾速度對(duì)溫升的影響較小。

(4)隨著軟輥彈性模量的增大，兩輥間接觸區(qū)油墨中間層的溫升逐漸提高。