UDK 519.61/.64:621.77 Original scientific article/Izvirni znanstveni članek ISSN 1580-2949 MTAEC9, 47(1)53(2013) CALCULATION OF THE LUBRICANT LAYER FOR A COARSE SURFACE OF A BAND AND ROLLS IZRAČUN SLOJA MAZIVA NA GROBI POVRŠINI TRAKU IN VALJEV Dušan Curcija1, Franc Vodopivec2, Ilija Mamuzic1 JCroatian Metallurgical Society, Berislaviceva 6, 10000 Zagreb, Croatia 2Institute of Metals and Technology, Lepi pot 11, 1000 Ljubljana, Slovenia mamuzic@simet.hr Prejem rokopisa — received: 2012-07-03; sprejem za objavo - accepted for publication: 2012-09-26 The effect of the average roughness of a lubricated band caused by dressing processes is analysed by applying the Reynolds differential equation for lubrication with the incorporated average roughness and evolution in the Fourier series to the third member. The analysis has shown that the average roughness has two effects on the lubricant-layer thickness in the entering section of the deformation zone. For a small surface roughness, the nominal lubricant-layer thickness decreases slowly (if the process is treated as occurring on a smooth surface) and the thickness grows again with an increase in the roughness. The basis for the analysis was the numerical Monte-Carlo method and the developed approximate analytical solution was in acceptable agreement with the numerical method. Keywords: surface roughness, lubricant-layer thickness, Reynolds equation, Monte-Carlo method, Fourier series Analiziran je vpliv povprečne hrapavosti mazanega traku pri procesih dresiranja. Podlaga analize je Reynoldsova diferencialna enačba za mazanje z vključeno povprečno hrapavostjo in obravnavo s Fourierovo vrsto do tretjega člena. Analiza je pokazala, da ima povprečna hrapavost dva učinka na debelino plasti maziva v vhodnem preseku področja deformacije. Pri majhni začetni hrapavosti se nominalna debelina plasti maziva počasi zmanjšuje (če se proces obravnava, kot da poteka na gladki površini) in znova raste, če se povečuje hrapavost. Podlaga za analizo je bila numerična metoda Monte Carlo, razvita pa je bila tudi približna analitična rešitev, ki se zadovoljivo ujema z numerično. Ključne besede: hrapavost površine, debelina plasti maziva, Reynoldsova enačba, metoda Monte Carlo, Fourierova vrsta 1 INTRODUCTION This technology is strongly associated with the quality of technological lubricants as it: • diminishes the contact friction, • removes the heat, cools the tool and diminishes the wear, • diminishes the deformation resistance and the deformation work, • diminishes the sticking to the tool and keeps the surface of the product clean. The basic groups examined in this work1-3 are: • liquid emulsions, • fats and compounds, • consistent lubricants, • transparent/glass lubricants, • powder lubricants and • metallic lubricants. Technological lubricants must meet a series of requirements, beginning with a high lubricity - the ability to form a flat, firm layer separating the contact surfaces - then there are thermal consistency and stability that prevent the damaging effect of the product corrosion, the properties not posing any health and environmental risks, etc. The liquid emulsions, whose compounds are mixtures of vegetable and mineral oils, are especially used in the cold rolling of 0.3-0.4 mm thick sheets and strips. In the cold rolling of sheets and strips, the dressing process is also used with an application of liquid lubricants to reduce undulation. 2 MATHEMATICAL MODELLING Mathematical modelling is a requirement of today's metallurgy4 5 and it is also used in the field of plastic deformation of metals. For an analysis of smooth surfaces6 7 the following equation is used: dp 6 ß(v o +v R) dx 12 ßQ £ 2(x) £ 2(x) £(x ) Q(x) = {udy=dß &£ 3(x)+ (1) (x) (2) The geometry of the lubricant contact8 and the length of the lubricant wedge are described with the relations (3), (4) and (5): £(x) = £ 0 + R cos a-, /1—I sin a — r (3) 2 Figure 1: Model of the tribomechanical system; 1- lubricant layer £(x) - nominal thickness for smooth surfaces, 2- band - in dressing processes the adhering angle a is low, 3- average band roughness d(x) -casual sheet roughness, 4- roll defined by surface smoothness. In Table 1 the roughness is Rz = 8 pm. Slika 1: Model tribomehanskega sistema; 1- plast maziva e(x) - nominalna debelina za gladke površine, 2- trak - pri procesu dresiranja pri majhnem kotu stika a, 3- povprečma hrapavost traku d(x) - slučajna hrapavost traku, 4- valj z definirano gladkostjo površine. V Tabeli 1 je njegova hrapavost Rz = 8 pm. a = R Î1-I cos a - -R+-R I - sin a 2 3 4 x ax x -(x) = -0 -ax3 0 2R 2R2 8 R3 (4) (5) I dp \ dx„ )= 6 o + Vr ) 1 1 I -2(xn)/ 1 \ \-3x)/ L -( xn )) = -( x)+ô( x) (6) (7) Reflexion of sheet roughness is added, as £o, to the lubricant wedge (4). The calculation is possible only with numerical mathematical methods and, in the program MATHEMATICA, the numerical method Monte Carlo was used. In the theoretical calculations regarding the model of the average roughness, the following function developed to the third term of Fourier series was applied: 4 1 1 o(x) = — (sin x+-sin3x+-sin5x) R (8) n 3 5 z 3 RESULTS AND DISCUSSION In Table 1 the standard values of geometrical, rheo-logical and kinematic characteristics of the processes of theoretical investigations are given according to the Russian-Ukrainian1011 authors. Table 1: Standard lubricant characteristics for theoretical calculations Tabela 1: Standardne značilnosti maziva za teoretične izračune For the average sheet roughness9, the mathematical relation in accordance with Figure 1 is: r h\ i Parameter y- piezo coefficient of viscosity po- rolling pressure vR- circumferential roll speed vo- sheet speed R- roll radius ^o- lubricant dynamic viscosity ^ = ^o exp (y * po) Barussa equation a- gripping angle £a- lubricant thickness on sheet A- technological parameter Value 2.18E-7 20E6 10 6 0.35 (0.25) 0.024-0.048 0-0.02 0.001-0.00001 1965512 (3934525) Unit Pa- Pa m/s m/s Pa s rad m- A = (1-exp(-y * p0)/6jU0y(v0+VR)) Rz - Rz = 1-10 ^m The parameters in Table 1 are of two groups: 1- lubricant rheological characteristics (po, y) 2- geometrical characteristics of the technological process (R, a, Rz) 3- kinematics (vo vr) m m Table 2: Lubricant-layer exit results (pm) Tabela 2: Izhodni rezultati za plast maziva (pm) 0 1 2 3 4 5 6 7 0 10.102 10.095 10.087 10.075 10.059 10.039 9.92 9.207 1 10.003 9.996 9.988 9.976 9.959 9.934 9.819 9.099 2 9.941 9.935 9.926 9.914 9.898 9.872 9.757 9.031 3 9.92 9.914 9.905 9.893 9.876 9.851 9.735 9.008 4 9.94 9.933 9.925 9.913 9.896 9.87 9.755 9.03 5 9.998 9.991 9.983 9.971 9.954 9.929 9.814 9.093 6 10.09 10.083 10.075 10.063 10.046 10.021 9.907 9.193 7 10.21 10.204 10.195 10.184 10.167 10.143 10.03 9.322 8 10.354 10.348 10.339 10.328 10.312 10.287 10.175 9.475 9 10.517 10.51 10.502 10.491 10.474 10.45 10.339 9.645 10 10.693 10.686 10.678 10.667 10.651 10.627 10.517 9.828 11 12 4- compounds (A, roughness space angle) The solutions of differential equation (6) are partially given in Table 2. The examined roughness is classified1213 in 10 vertical classes and the band profile roughness in 32 horizontal classes. In principle, with a decreasing band-lubricant thickness (£a in Figure 1) the lubricant thickness in the entering section of the metal deformation zone is also decreased (£0). As shown in14, the lubricant wedge has the ideal geometry and can give economic savings of the lubricant in the metalworking technology. The numerical integration of equation (6) was checked with the approximate1517 analytical solutions possible in the case of practical interest, which is found in equations (9), (10a)-(10e) and (11). Equation (9) is the simplest analytical solution that does not consider the thickness of the band lubricant layer, £a >> £0. With a clear complexity, equation (11) corrects this deficiency: 315 AR3 a7 -168R2 a4 -1824(52 = 0 1-exp(-yp0) (9) £ 0 = 0.5Ra2 ..A = 6 /j.y(v o + Vr ) W = A- 3a'R -- 4- Rjl - I cos a+ | - R(sin a—a) (10a) -W, = 3<5' 1 6 TOR - 1 6- Rj 1 - I cos a - -£a + I - R(sin a-a) W = £ Q+3£0Ô2 3 £ 2+6d2 R3 5aR r n ^t ]5 5 R 1 - | cos a + - R(sin a— a)| M 1 R R J R5 R^j1 -cos a - £a + - R(sin a—a) W1 + W2 + W3 * (W4 + W5) =0 (10b) (10c) (10d) (10e) (11) In Table 3 approximate numerical and analytical solutions are compared. The approximate numerical solutions can be compared with numerical integration only for the entering roughness profile, thus, at the entering section of the deformation zone with x = 0. It is clear from Table 3 that the simple analytical form of equation (9) with numerous approximations describes well the lubricant layer for the case of a lubricant excess on the sheet and the rolls. Table 3: Comparison of approximate analytical and numerical Monte-Carlo solutions for one point of the graph crossing from Figure 2 Tabela 3: Primerjava približnih analitičnih in numeričnih rešitev Monte Carlo za eno točko prereza grafa na sliki 2 Case conditions Approximate analytical solutions, eq. (11) and (9) Monte-Carlo method, eq. (6) x = 0 (initial roughness profile) Rz= 1 Mm R ~ 6 5 a = 0.00918759 rad A = 1965512 m-1 R = 0.35 m £a = 0.001 m £0 = 14.721 [m (11) £0 = 14.771 [m (9) £a = 0.0001 m £0 = 13.834 [m (11) £a = 0.001 m £0 = 14.772 [m £a = 0.0001 m £0 = 13.761 [m x = 0 (initial roughness profile) Rz= 10 |im Rz « 6 5 az= 0.0092867 rad A = 1965512 m-1 R = 0.35 m £a = 0.001 m £0 = 15.024 [m (11) £0 = 15.092 [m (9) £a = 0.0001 m £0 = 13.511 [m (11) £a = 0.001 m £0 = 15.077 [m £a = 0.0001 m £0 = 13.429 [m x = 0 (initial roughness profile) Rz= 10 |m Rz « 6 5 az= 0.00840867 rad A = 3934525 m-1 R = 0.25 m £a = 0.001 m £0 = 8.776 [m (11) £0 = 8.838 [m (9) £a = 0.0001 m £0 = 8.464 [m (11) £a = 0.001 m £0 = 8.755 [m £a = 0.0001 m £0 = 8.429 [m Table 4: Effect of the two-sided roughness of the sheet and rolls, congruous for £o Tabela 4: Vpliv dvostranske hrapavosti traku in valjev, kongruenten za £o x = 0 (initial roughness profile) = 10 |m, average roughness, horizontal (transversal) = 8 |m, longitudinal roll roughness = 6 d (GOST 2789-73) a = 0.00840867 rad A = 3934525 m-1 R = 0.25 m R R R Monte-Carlo method £a = 0.0001 m £0 = 9.299 [m £0 = 8.429 [m One-sided Two-sided roughness of the roughness of the sheet sheet and roll ea = 0.0001 m £a = 0.0001 m £0 = 8.429 |im e0 = 8.919 |im Rz ^ 0 £a = 0.0001 m £0 = 7.877 |im The longitudinal band profile on abscissa in shown in 66 classes and on ordinate in 11 classes for roughness (0-10 |im). It is useful to calculate the lubricant thickness £0 in the range of 8.5-12.5 |im in the area of I-I. Q, K and W designations connect the specific areas of the network diagram with the contour plot (an aircraft picture of the network diagram). B and C are the left and right sides of the band roughness defined as a sine evolution function in the Fourier series: B in the range of (w-2w) rad and C in (0-w) rad. Line P in Figure 2 represents the nominal lubricant-layer thickness on side C, thus, by having the thickness for Rz ~ 6 |im, an equivalent to the lubricant-layer thickness on a smooth surface is obtained. Side B does not have this property. 4 R R y Figure 2: Effect of the average sheet roughness and smooth rolls on £0 Slika 2: Vpliv povprečne hrapavosti traku in gladkih valjev na £0 In Figure 3 both sides of the roll longitudinal roughness C from Figure 2 are shown. The average roughness conserves the same properties as in Figure 2. The longitudinal roughness profile in the range of classes 33 to 66 gives a more stable hydrodynamic lubrication, while for classes 1 to 33 the hydrodynamic lubrication is already seriously impaired by the low roughness of the band and rolls. The lubricant layer decreases rapidly and spreads to fractal areas. A stable lubrication can be achieved on small band segments and around class 4 of the longitudinal sheet profile and around classes 10 and 30. The complex shapes of the lubrication space are probably determined by the band and roll roughness in the entering section of the deformation zone that determines a different lubrication layer than in the case of smooth-sheet and roll surfaces. 4 CONCLUSIONS Based on the results of theoretical analyses of the effect of the band roughness on the lubrication dressing processes, the following conclusions are proposed: • The average band roughness has a critical value when it starts to affect positively the lubricant layer with its increase in comparison with a smooth surface. Up to line P in Figure 2, the lubricant layer has a tendency to increase and to decrease the formation of sunk baskets in area Q. The theoretical explanation for this is that the surface roughness determines the shape of the lubricant layer for every value of Rz. This is the range of a stable lubrication. • If congruous roll roughness is added to the average band roughness, forming a longitudinal roll roughness with the positive side in the range of (0-^), the thickness of the lubricant layer in the entering section of the band deformation zone will increase its longitudinal profile from class 33 to 66 (Figure 3 and Table 4) and will approach the boundary lubrication. • The developed approximate analytical solutions agree with the numerical integration of equation (6) and ensure a reliable approach to the analysis. • If the technological process was performed with a nominal lubricant-layer thickness marked with line P in Figure 2 the best roll rhythm would be obtained without significant fluctuations of the lubricant thickness, especially in the case of the boundary-lubrication proximity. This includes the control of the roll roughness. 5 SYMBOLS AND FIGURES Figure 3: Effect of the average sheet roughness and longitudinal roll roughness on £o (Table 4) Slika 3: Vpliv povprečne hrapavosti traku in longitudinalne hrapavosti valjev na £0 (tabela 4) Symbol Unit Comment £o m, (¡¡m) Lubricant thickness in the entering section of the deformation zone (Figure 1) e(x) m Lubricant thickness in the range of [-a : 0], Figure 1, equations (3) and (5) £a m Lubricant thickness ahead of the entering section of the deformation zone a m Length of the lubricant wedge (Figure 1), equation (4) a rad Band dressing angle VR m/s Circumferential roll speed VT m/s Mandrel speed R m Roll radius Rz m Roughness of the band surface, equation (8) d2 Dispersion roughness of the sheet and rolls according to equation (9) dx Casual lubricant thickness depending on the band roughness (and rolls) < > Operative mathematical expectation x, y Descartes coordinates Q(x) - Volume use of lubricant (on the band perimeter) V0 Pa s Lubricant dynamic viscosity by the rolling pressure v Pa s Lubricant dynamic viscosity by the air pressure u m/s Lubricant rate on the abscissa V m2/N Piezo coefficient of lubricant viscosity P Pa Rolling pressure Q m2/s Use of lubricant on the mandrel perimeter - a one-dimensional model dp/dx Pa/m Pressure gradient in the lubricant layer, equation (1) sin a rad Marking the trigonometric function for the griping alpha angle H m Enter band thickness h m Exit band thickness A m-1 Technological parameter: A = [1- exp(-yp)] / [6^0y(vR + v0) ] exp, n 2.718 Base of natural logarithm (3.141) 14 1-1 Reference 1 |m 10-6 m Micrometre |m Band- and roll-roughness classes L |m Longitudinal holding-band profile Q, K, W Markers for Figure 2 6 REFERENCES 11. Mamuzic, V. M. Drujan, Teorija, Materijali, Tehnologija celičnih cijevi, Hrvatsko Metalurško Društvo, Zagreb 1996, 428-435 2 D. Curčija, Mater. Tehnol., 37 (2003) 5, 237-254 3 D. Curčija, I. Mamuzic, F. Vodopivec, Metalurgija, 45 (2006) 3, 250 (Summary) 4 A. I. Gubin, B. B. Veselovskiy, D. Curčija, A. A. Kochubey, Mathematical simulation and choice of optimum thermal models of continuous events, Metalurgija, 47 (2008) 3, 255 (Summary) 5 Iu. V. Brazaluk, O. O. Kochubey, D. Curčija, M. V. Polyakov, D. V. Yevdokymov, On a mathematical model of particle in liquid metal, Metalurgija, 47 (2008) 3, 256 (Summary) 6D. Curčija, I. Mamuzic, Mater. Tehnol., 43 (2009) 1, 23-30 7 D. Curčija, I. Mamuzic, Goriva i maziva, 48 (2009) 1, 3-28 8 D. Curčija, I. Mamuzic, Mater. Tehnol., 42 (2008) 2, 59-63 9 D. Curčija, I. Mamuzic, Metalurgija, 44 (2005) 4, 295-300 10 O. P. Maksimenko, N. P. Podberezniij, Izvestija Černaja metallurgija, 73 (2003) 10, 12-16 11 O. P. Maksimenko, A. A. Semenča, Sučasni problemi metalurgii, 8 (2005), 99-103 12 O. P. Maksimenko, O. E. Lejko, Sučasni problemi metalurgii, 8 (2005), 93-99 13 P. I. Klimenko, Sučasni problemi metalurgii, 8 (2005), 44-49 14S. M. Ionov, V. I. Kantorovič, S. A. Sepovalov, A. N. Krjukov, Sučasni problemi metalurgii, 8 (2005), 224-228 15 D. M. Me, S. P. Liu, J. F. Zheng, Met. Form. Technol., 20 (2002) 5, 29-32 16 Y. T. Keun, B. H. Lee, R. H. Wagner, J. Mater. Process. Technol., 130 (2002), 60-63 17 D. Curčija, I. Mamuzic, Goriva i maziva, 46 (2007) 1, 23-44