검색 전체 메뉴
PDF
맨 위로
OA 학술지
미세 그루브가 있는 무한폭 Slider 베어링의 윤활해석:제2보 - 그루브 깊이의 영향 Lubrication Analysis of Infinite Width Slider Bearing witha Micro-Groove: Part 2 - Effect of Groove Depth
  • CC BY
  • CC BY
ABSTRACT

It is currently well known that surface textures act as lubricant reservoirs, entrap wear debris, and hydrodynamic bearings, which can lead to certain increases in load-carrying capacities. Until recently, the vast majority of research has focused on parallel sliding machine components such as thrust bearings, mechanical face seals, piston rings, etc. However, most sliding bearings have a convergent film shape in the sliding direction and their hydrodynamic pressure is mainly generated by the wedge action. Following the first part of the present study that investigates the effect of groove position on the lubrication performances of inclined slider bearings, this paper focuses on the effects of groove depths and film thicknesses. Using a commercial computational fluid dynamics (CFD) code, FLUENT, the continuity and Navier–Stokes equations are numerically analyzed. The results show that the film thickness and groove depth have a significant influence on the pressure distribution. The maximum pressure occurs at the groove depth where the vortex is found and, as the depth increases, the pressure decreases. There is also a groove depth to maximize the supporting load with the film thickness. The friction force acting on the slider decreases with deeper grooves. Therefore, properly designed groove depths, depending on the operating conditions, can improve the load-carrying capacity of inclined slider bearings as compared to the bearings without a groove.


KEYWORD
슬라이더 베어링 , 표면조직가공 , 나비에-스톡스 방정식 , 유체윤활 , 수치해석
  • 1. 서 론

    기계의 상대운동면에서 발생하는 마찰과 마멸(Wear)을 저감시킬 수 있는 가장 보편적인 방법은 적절한 점도의 윤활유를 사용하여 운전에 적합한 방법으로 윤활하는 것이다[1]. 하지만 유체의 전단으로 인한 점성마찰은 피할 수 없기 때문에 효율향상에는 한계가 있다. 현재 기계의 효율과 내구수명을 보다 향상시키기 위하여 전세계에서 집중적으로 연구중인 대표적인 3가지의 핵심분야는 새로운 저마찰 윤활제, 저마찰/저마멸 코팅(Coating) 기술 및 Surface texturing 기술의 개발과 이들의 실제적용이다[2].

    이 중에서 20여년 전에 Etsion[3]에 의해서 체계적인 연구가 시작된 Surface texturing은 윤활면에 미세한 크기의 그루브(Groove), 딤플(Dimple) 등을 가공하여 마찰 저감과 함께 하중지지능력까지 향상시킬 수 있는 최신의 표면처리기술이다[4-6]. Brizmer 등[7]은 윤활면이 평행한 경우에도 부분적인 Texturing은 상당한 정도의 하중지지가 가능함을 레이놀즈 방정식에 대한 수치해석으로 규명하였다. 현재에는 이 기술의 적용으로 경계윤활 상태로 운전된다고 알려진 평면 스러스트 베어링(Plain thrust bearing), 메카니컬 시일(Mechanical seal), 피스톤링(Piston ring) 등의 수명을 크게 향상시켰다. 특히, 유막두께에 비하여 Texturing 깊이가 상대적으로 얕을수록 캐비테이션(Cavitation) 발생이 강하며, 이것으로 인한 입구부에서의 윤활유 흡입(Inlet suction)이 압력발생에 주로 기여하는 것으로 인정되고 있다[8-9]. 이와 같이 윤활성능에 영향을 미치는 중요한 인자중의 하나인 Texturing 깊이에 대하여 많은 연구가 수행되고 있다[3, 5-6, 10-14]. 특히, Jeong & Park[12], Park & Kim[14]은 구형딤플(Dimple)의 깊이와 온도경계조건이 평행 스러스트 베어링의 열유체(Thermohydrodynamic: THD) 윤활성능에 미치는 영향을 조사하였다. 이와 같이 지금까지의 Texturing 깊이에 대한 연구는 윤활면이 평행한 문제에 거의 국한되었기에 유막이 축소되는(Convergent film) 대부분의 미끄럼 베어링에 대한 연구는 크게 미흡한 실정이다.

    한편, Texturing한 경우에는 하지 않은 경우의 유막두께에 비하여 상대적으로 깊고 유막이 급격하게 변하므로 유체의 관성에너지를 무시할 수가 없다. 따라서 레이놀즈 방정식 보다는 Navier-Stokes 방정식과 연속방정식을 사용하여 수치해석하는 것이 보다 적절하며[15], 이러한 경우에는 전산유체역학(Computational fluid dynamics: CFD) 해석방법을 사용하는 것이 현실적이다. Sahlin 등[11]은 평행한 윤활면에서 그루브의 폭과 깊이가 윤활성능에 미치는 영향을 조사하기 위하여 Navier-Stokes 방정식을 해석하였다. 베어링의 하중지지능력이 최대인 그루브 깊이가 존재하며, 와류(Vortex)는 이 부근의 깊이에서 발생하였다. Cupillard 등[16]은 입구부에 위치한 사각형 그루브 3개의 깊이가 경사진 슬라이더 베어링의 윤활성능에 미치는 영향을 조사하였다. 그루브가 없는 경우에 대한 상대적인 하중지지능력은 유막의 수축비(Convergence ratio)가 아주 작을 경우에는 그루브가 깊어질수록 향상되는 반면에 큰 경우에는 모든 깊이에서 감소하였다. 그리고 깊은 그루브에서는 와류가 발생하였다. Papadopoulos 등[17]은 CFD 해석방법과 유전 알고리즘(Genetic algorithm)을 사용하여 다수의 그루브가 있는 경사진 3D 슬라이드 베어링 모델에서 최적인 그루브 깊이 등을 구하였다. 이상에서 고찰에서와 같이 Texturing 깊이가 경사진 슬라이더 베어링의 윤활성능에 미치는 영향에 대한 연구결과는 크게 부족한 실정이다.

    본 논문에서는 Surface texturing한 슬라이더 베어링의 윤활성능을 향상시키기 위한 연구의 일환으로 그루브의 위치에 대한 제1보[18]에 이어서 깊이에 의한 영향을 상용 CFD 프로그램을 사용하여 조사하고자 한다.

    2. 수치해석

       2-1. 해석 모델

    Fig. 1에는 본 논문에서 사용한 무한폭 슬라이더 베어링에 대한 수치해석모델을 개략적으로 나타내었다[18]. 한 개의 그루브가 있는 패드(Pad)(길이 L)에 대하여 슬라이더는 상대속도 U로 운동한다. 여기서 a는 입구에서 사각형 단면의 그루브(폭 c, 깊이 hp)가 중앙에 위치한 미 세 베어링 셀(Cell)(길이 b)까지의 거리이다. N은 결과의 설명을 위하여 도입한 그루브의 위치(Position)를 나타내는 숫자이다. 해석에서 사용한 위치는 N = 4로 고정하였으며, 이 경우에 a = L-3b이다.

       2-2. 지배방정식 및 경계조건

    유막내의 윤활유 유동이 정상상태, 비압축성, 층류유동인 경우에 대한 연속방정식과 Navier-Stokes 방정식은 각각 다음과 같다.

    image
    image

    베어링의 입·출구부에는 식 (3)의 경계압력이 작용하며, 모든 베어링면에서는 점착조건(No-slip condition)을 만족한다.

    image

       2-3. 수치해석

    본 논문에서 수치해석에 사용한 격자계는 참고문헌[18]과 거의 동일하며 수치해석은 상용 CFD 프로그램인 FLUENT[19]를 사용하였다. Table 1에는 해석에서 사용한 베어링의 상세사양, 운전조건 및 윤활유의 물성자료를 나타내었다.

    [Table 1.] Specification of inclined slider bearing and operating condition

    label

    Specification of inclined slider bearing and operating condition

    본 논문에서는 그루브의 폭과 위치, 유막의 기울기 (h0-hn)를 고정하고 그루브 깊이 hp 및 입·출구부 유막두께 h0, hn 의 변화에 따른 영향을 조사하였다.

    3. 결과 및 고찰

    Figs. 2~3에는 h0 = 2.0 μm, hn = 1.0 μm인 경우에 hp 의 변화에 따른 압력과 속도의 분포를 각각 나타내었다. Fig. 4Fig. 3의 그루브내에서의 속도분포를 확대한 그림이다. 이상의 그림들은 이해를 쉽게 하기 위하여 유막두께 방향으로 크게 확대하였다. 이들과 동일한 순서로 Figs. 5~7h0 = 3.0 μm, hn = 2.0 μm인 경우에 대한 그림으로 Figs. 2~4 보다 1.0 μm 두껍지만 유막의 기울기는 동일하다. 따라서 유막의 쐐기효과(Wedge effect)가 상대적으로 많이 약화되므로 발생압력은 크게 감소하였다. 그루브가 없는 경우(hp = 0 μm)보다 그루브가 특정깊이인 경우에 상대적으로 고압이 발생하는 점이 특징적이다. 베어링 입구부에서는 역류가 전혀 발생하지 않지만 그루브 내에서는 유막두께와 그루브 깊이에 따라서 와류(Vortex)가 발생할 수 있음을 알 수 있다(Figs. 4 & 7). 즉, h0 = 2.0 μm, hn = 1.0 μm인 경우에는 hp = 0.5 μm에서 와류가 나타나지만 h0 = 3.0 μm, hn = 2.0 μm인 경우에는 hp = 1.0 μm에서도 이것이 관찰되지 않았다. 이상의 결과에서 그루브에 의한 압력상승이 그루브내에서의 와류발생과 밀접한 관련이 있을 것으로 추정된다[16]. 참고로 와류 발생이 없는 Fig. 7(a)에서의 (−)속도는 사각형인 그루브 바닥의 양쪽 모서리부에서 발생하는 아주 작은 모서리 유동(Corner flow)으로 인한 결과이다.

    Figs. 8~9에는 hp에 따른 그루브 중앙(X = 0.35)에서의 속도분포와 길이방향 압력분포를 h0 = 2.0 μm, hn = 1.0 μm 와 h0 = 3.0 μm, hn = 2.0 μm에 대하여 각각 나타내었다. Fig. 8에서 속도분포가 모두 오목한(Concave) 것은 이 위치보다 하류가 고압임을 의미한다(Fig. 2 & 5). 그루브 바닥에서 유동방향의 속도구배가 0인 hp 가 각각 0.3 μm (Fig. 8(a)), 1.0 μm(Fig. 8(b))이며, 이 보다 더 깊은 경우에 나타나는 (−)속도는 와류가 발생한다는 증거이다. 그루브가 없는 경우 보다 있는 경우에 더 높은 압력이 발생할 수 있음을 Fig. 9에서 보여주고 있다. 특히, 최고압력이 발생하는 hpFig. 8에서 고찰한 그루브 바닥에서의 속도구배가 0인 결과와 일치하였다. 이는 그루브내에서 슬라이더의 운동으로 유동하는 유체의 관성에너지 일부가 압력으로 변환되기 때문으로 추정된다. 이 보다 그루브가 깊어질수록 와류(Fig. 4 & 7)의 성장에 더 많은 관성에너지가 소모되므로 압력은 점차 낮아진다.

    Fig. 10에는 지지하중과 슬라이더에 작용하는 마찰력에 미치는 hp의 영향을 나타낸 것으로 각 유막두께에서 그루브가 없는 경우의 결과로 무차원화하였다. 그루브가 깊어질수록 증가하던 지지하중은 최대값에 도달한 다음에는 거의 직선적으로 감소하였다. 이러한 결과는 Fig. 9에서 고찰한 것과 같이 와류의 생성 및 이의 크기 증가에 따른 관성에너지의 손실에 기인하는 것으로 추정된다. 특히, 지지하중이 최대로 되는 그루브 깊이가 유막 두께의 크기에 따라서 달라지는 점이 특징적일 뿐만 아니라 이의 상대적인 효과에도 큰 차이가 있다. 한편, 마찰력은 유막두께의 크기에 거의 관계없이 그루브가 깊어질수록 감소하는 것은 Fig. 8에 나타낸 것과 같이 속도구배의 감소로 인한 것으로 이해된다.

    이상의 결과에서 경사진 슬라이더 베어링인 경우에도 그루브의 깊이를 적절하게 설계하면 상당한 정도의 하중지지능력 향상이 가능함을 확인하였다. 따라서 슬라이더 베어링의 윤활성능을 향상시키기 위해서는 다양한 Surface texturing 관련인자와 운전조건 등에 대한 추가 연구가 요구된다.

    4. 결 론

    본 논문에서는 그루브의 깊이가 Surface texturing한 슬라이더 베어링의 윤활성능에 미치는 영향을 조사하기 위하여 미세 베어링 모델을 상용 CFD 프로그램으로 수치해석하였다. 해석결과를 요약하면 다음과 같다.

    1. 그루브내에서 와류가 발생하려는 깊이에서 최고압력이 발생하였으며, 이 보다 더 깊어지면 압력은 감소하였다.

    2. 유막두께의 크기에 따라서 지지하중이 최대인 그루브 깊이가 존재하였다.

    3. 운전조건에 따라서 그루브 깊이를 적절하게 설계하면 그루브가 없는 경우보다 하중지지능력을 향상시킬 수 있다.

    4. 그루브는 슬라이더에 작용하는 마찰력을 감소시키며 그루브가 깊을수록 이의 효과는 크게 된다.

    5. 그루브 깊이가 경사진 슬라이더 베어링의 윤활성능에 미치는 영향이 상당하다는 결과를 처음으로 제시하였다.

    본 논문의 결과는 각종 미끄럼 베어링의 윤활성능향상을 위한 설계자료로 사용될 수 있으며, 다양한 조건에 대한 추가연구가 요구된다.

    Nomenclature

    a: Distance to the bearing cell, see Fig. 1 b: Bearing cell size (μm), see Fig. 1 c: Groove width (μm), see Fig. 1 h0: Inlet film thickness (μm), see Fig. 1 hn: Outlet film thickness (μm), see Fig. 1 hp: Groove depth (μm), see Fig. 1 L: Bearing length N: Number indicating groove position, see Fig. 1 p: Pressure (Pa) P0, Pn: Boundary pressure (Pa) U: Sliding speed of slider (m/s) : Flow velocity vector x, y: Cartesian coordinates X: Dimensionless x-coordinate, x/L η: Absolute oil viscosity (kg/m ∙ s) ρ: Oil density (kg/m3)

참고문헌
  • 1. Cameron A. 1965 The principle of lubrication google
  • 2. Holmberg K., Andersson P., Erdemir A. 2012 Global energy consumption due to friction in passenger cars [Tribol. Int.] Vol.47 P.221-234 google cross ref
  • 3. Etsion I. 2005 State of the art in laser surface texturing [ASME J. Tribol.] Vol.127 P.248-253 google cross ref
  • 4. Sudeep U., Tandon N., Pandey R. K. 2015 Performance of lubricated rolling/sliding concentrated contacts with surface textures: A review [J. Tribol.] Vol.137 P.031501 google cross ref
  • 5. Ibatan T., Uddin M. S., Chowdhury M. A. K. 2015 Recent development on surface texturing in enhancing tribological performance of bearing sliders [Sur. Coating Technol.] Vol.272 P.102-120 google cross ref
  • 6. Gropper D., Wang L., Harvey T. J. 2016 Hydrodynamic lubrication of textured surfaces: A review of modeling techniques and key findings [Tribol. Int.] Vol.94 P.509-529 google cross ref
  • 7. Brizmer V., Kligerman Y., Etsion I. 2003 A laser surface textured parallel thrust bearing [Tribol. Trans.] Vol.46 P.397-403 google cross ref
  • 8. Olver A. V., Fowell M. T., Spikes H. A., Pegg I. G. 2006 ‘Inlet suction’, a load support mechanism in non-convergent, pocketed, hydrodynamic bearings [Proc. IMechE., J. Eng. Tribol.] Vol.220 P.105-108 google cross ref
  • 9. Fowell M., Olver A. V., Gosman A. D., Spikes H. A., Pegg I. 2007 Entrainment and inlet suction: Two mechanisms of hydrodynamic lubrication in textured bearings [ASME J. Tribol.] Vol.129 P.221-230 google cross ref
  • 10. Brajdic-Mitidieri P., Gosman A. D., Ioannnides E., Spikes H. A. 2005 CFD analysis of a low friction pocketed pad bearing [ASME J. Tribol.] Vol.127 P.803-812 google cross ref
  • 11. Sahlin F., Glavatskih S. B., Almqvist T., Larsson R. 2005 Two-dimensional CFD-analysis of micro-patterned surfaces in hydrodynamic lubrication [ASME J. Tribol.] Vol.127 P.96-102 google cross ref
  • 12. Jeong Y. H., Park T. J. 2014 THD analysis of surface textured parallel thrust bearing: Effect of dimple radius and depth [J. Korean Soc. Tribol. Lubr. Eng.] Vol.30 P.303-310 google cross ref
  • 13. Caramia G., Carbone G., De Palma P. 2015 Hydrodynamic lubrication of micro-textured surfaces: Two dimensional CFD-analysis [Tribol. Int.] Vol.88 P.162-169 google cross ref
  • 14. Park T. J., Kim M. G. 2017 Thermohydrodynamic lubrication of surface-textured parallel slider bearing: effect of dimple depth [J. Korean Soc. Tribol. Lubr. Eng.] Vol.33 P.288-295 google
  • 15. Rom M., Muller S. 2018 An effective Navier-Stokes model for the simulation of textured surface lubrication [Tribol. Int.] Vol.124 P.247-258 google cross ref
  • 16. Cupillard S., Cervantes M. J., Glavatskih S. 2008 Pressure buildup mechanism in a textured inlet of a hydrodynamic contact [ASME J. Tribol.] Vol.130 P.1-10 google
  • 17. Papadopoulos C. I., Efstathiou E. E., Nikolakopoulos P. G., Kaiktsis L. 2011 Geometry optimization of textured three-dimensional micro-thrust bearings [J. Tribol.] Vol.133 P.041702 google cross ref
  • 18. Park T. J., Jang I. G. 2019 Lubrication analysis of infinite width slider bearing with a micro-groove: Part 1 - Effect of groove position [Tribol. Lubr.] Vol.35 P.376-381 google
  • 19. 2011 ANSYS FLUENT User Guide, Release 14.0 google
이미지 / 테이블
  • [ Fig. 1. ]  Schematic of 2D inclined slider bearing model.
    Schematic of 2D inclined slider bearing model.
  • [ ] 
  • [ ] 
  • [ ] 
  • [ Table 1. ]  Specification of inclined slider bearing and operating condition
    Specification of inclined slider bearing and operating condition
  • [ Fig. 2. ]  Contour plot of pressure distribution: h0 = 2.0 μm, hn = 1.0 μm. hp is (a) 0 μm, (b) 0.5 μm, (c) 1.0 μm.
    Contour plot of pressure distribution: h0 = 2.0 μm, hn = 1.0 μm. hp is (a) 0 μm, (b) 0.5 μm, (c) 1.0 μm.
  • [ Fig. 3. ]  Contour plot of velocity distribution: h0 = 2.0 μm, hn = 1.0 μm. hp is (a) 0 μm, (b) 0.5 μm, (c) 1.0 μm.
    Contour plot of velocity distribution: h0 = 2.0 μm, hn = 1.0 μm. hp is (a) 0 μm, (b) 0.5 μm, (c) 1.0 μm.
  • [ Fig. 4. ]  Contour plot of velocity distribution in groove: h0 = 2.0 μm, hn = 1.0 μm. hp is (a) 0.5 μm, (b) 1.0 μm.
    Contour plot of velocity distribution in groove: h0 = 2.0 μm, hn = 1.0 μm. hp is (a) 0.5 μm, (b) 1.0 μm.
  • [ Fig. 5. ]  Contour plot of pressure distribution: h0 = 3.0 μm, hn = 2.0 μm. hp is (a) 0 μm, (b) 1.0 μm, (c) 1.5 μm.
    Contour plot of pressure distribution: h0 = 3.0 μm, hn = 2.0 μm. hp is (a) 0 μm, (b) 1.0 μm, (c) 1.5 μm.
  • [ Fig. 6. ]  Contour plot of velocity distribution: h0 = 3.0 μm, hn = 2.0 μm. hp is (a) 0 μm, (b) 1.0 μm, (c) 1.5 μm.
    Contour plot of velocity distribution: h0 = 3.0 μm, hn = 2.0 μm. hp is (a) 0 μm, (b) 1.0 μm, (c) 1.5 μm.
  • [ Fig. 7. ]  Contour plot of velocity distribution in groove: h0 = 3.0 μm, hn = 2.0 μm. hp is (a) 1.0 μm, (b) 1.5 μm.
    Contour plot of velocity distribution in groove: h0 = 3.0 μm, hn = 2.0 μm. hp is (a) 1.0 μm, (b) 1.5 μm.
  • [ Fig. 8. ]  Velocity distribution at X = 0.35. (a) h0 = 2.0 μm, hn = 1.0 μm, (b) h0 = 3.0 μm, hn = 2.0 μm.
    Velocity distribution at X = 0.35. (a) h0 = 2.0 μm, hn = 1.0 μm, (b) h0 = 3.0 μm, hn = 2.0 μm.
  • [ Fig. 9. ]  Pressure distribution. (a) h0 = 2.0 μm, hn = 1.0 μm, (b) h0 = 3.0 μm, hn = 2.0 μm.
    Pressure distribution. (a) h0 = 2.0 μm, hn = 1.0 μm, (b) h0 = 3.0 μm, hn = 2.0 μm.
  • [ Fig. 10. ]  Effect of groove depth on the dimensionless (a) load-carrying capacity, (b) friction force.
    Effect of groove depth on the dimensionless (a) load-carrying capacity, (b) friction force.
(우)06579 서울시 서초구 반포대로 201(반포동)
Tel. 02-537-6389 | Fax. 02-590-0571 | 문의 : oak2014@korea.kr
Copyright(c) National Library of Korea. All rights reserved.