Advanced
Effect of Contact Stiffness on Brake Squeal Analysis Using Analytical FE Squeal Model
Effect of Contact Stiffness on Brake Squeal Analysis Using Analytical FE Squeal Model
Transactions of the Korean Society for Noise and Vibration Engineering. 2014. Oct, 24(10): 749-755
Copyright © 2014, The Korean Society for Noise and Vibration Engineering
  • Received : July 04, 2014
  • Accepted : September 26, 2014
  • Published : October 20, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
재 영 강
Corresponding Author ; Member, Division of Automotive & Mechanical Engineering, Kongju National University E-mail :jkang@kongju.ac.krTel : +82-41-521-9263

Abstract
The analytical-finite element(FE) squeal model is applied to investigate the squeal propensity associated with contact stiffness of the disc brake system. The system contact stiffness is incorporated into the perturbed equations of motion in the analytical manner where the brake components are modeled by FE method. The results show that the contact stiffness of the friction material and the contact stiffness between the pads and caliper are the influential factors on the squeal propensity. Particularly, the modal instability of the 3200 Hz squeal mode drastically changes with respect to the contact stiffness between the pads and caliper.
Keywords
1. 서 론
브레이크 스퀼(squeal) 분야는 자동차 및 진동 분야에서 꾸준히 연구되어 왔다. 특히 자동차 브레이크 스퀼 연구시 시스템 형상의 복잡성으로 유한요소법(finite element method)을 이용하여 선형불안정성을 계산하는 복소수 고유치 해석(complex eigenvalue analysis)을 주로 이용하여 왔다 (1~6) .
하지만 상용유한요소 소프트웨어를 이용하여 연구할 때는 항상 제약이 따르게 된다. 브레이크 스퀼 해석은 일반적으로 상용 유한요소 소프트웨어를 이용하여 수행되고 있으나, 모델링 기법이나 알고리즘에 대해 원개발자 이외의 일반 연구자들이 접근하기 힘들고 해석시간(computation time)이 상대적으로 길어서 시스템 파라미터 연구를 수행하는데 어려움이 있다. 따라서 실제 시스템의 스퀼 시험 결과와의 오차를 근본적으로 유추하기가 어렵다. 이러한 유한요소기법 자체의 한계 및 문제점들을 극복하기 위해서는 현재 상용소프트웨어의 스퀼 알고리즘의 문제점을 파악하여야 한다.
첫째 기존 상용소프트웨어가 제공하는 스퀼 해석모델은 접촉 모델이 단순화되어 있다. Kang (7~9) 은 이를 개선한 정교 마찰 접촉 및 디스크 회전모델을 수학-유한요소 융합모델로 제시하여 이를 극복할 수 있음을 보였다.
둘째 기존 상용소프트웨어를 이용하여 브레이크 시스템을 구축할 때 각 단품 간 접촉되는 지점에서의 접촉강성이나 접촉 경계면 조건을 어떻게 설정할지 기준이 없다. 특히 접촉강성은 접촉투과법(penalty method)에서 제공하는 기준값(default)을 이용하기 마련이다. Nam 등은 캘리퍼 접촉강성에 대한 스퀼 연구를 기존 상용소프트웨어를 이용하여 수행하였다. 하지만 해석 시간의 제약으로 인해서 몇 개의 접촉 강성 조건에 대해서만 조사하였고 캘리퍼 모델은 고려하지 않았다 (10) .
이 연구에서는 수학-유한요소 융합모델 (7~9) 을 이용 하여 단품 간 접촉부 강성에 따른 스퀼 성향을 연구하였다. 기존 상용소프트웨어는 해석 시간 때문에 시스템 파라미터 연구에 제약이 있지만, 스퀼 융합 모델은 접촉모델링의 정교성 이외에도 해석 시간이 대단히 짧다는 장점이 있어서 수만 번의 해석을 단 몇 시간 내에 수행할 수 있다. 따라서 이러한 스퀼 융합모델을 이용하여 디스크 브레이크의 접촉부 강성의 스퀼 영향도를 연구하고자 하였다.
2. 운동방정식
디스크 브레이크 모델을 구축하기 위해서 디스크, 패드 및 캘리퍼의 유한요소법(FEM)으로 모델링을 구성한다. Fig. 1 에서 보듯이 각각의 단품을 이용하여 브레이크 시스템을 구축할 때 서로 접촉 결합이 이루어진다. 이러한 접촉 결합 상황을 구현할 때 접촉투과법(penalty method)를 이용하여 두 접촉 파트간 접촉강성을 삽입하여 모델링 하여야 한다. 만일 패드와 디스크간 접촉강성(contact stiffness, kc )없이 강체 접촉한다면 미소 진동에 의한 접촉을 표현하기 힘들고, 접촉강성에 의한 모드연성이 발생하지 아니하여 모드연성(mode-coupling) 스퀼 메커니즘을 표현할 수 없게 된다 (11~12) . 특히 캘리퍼와 패드의 사이드(side)부에서 브레이크 작동 시 접촉이 발생하는데 이를 접촉투과법(penalty method)으로 접촉모델링 할 경우 접촉부에 대한 사이드 접촉강성( kb )을 정의하기가 모호하지만 스퀼 특성에 큰 영향을 줄 수 있다. 이러한 접촉강성부의 스퀼에 대한 영향도를 살펴보고자 Fig. 1 Fig. 2 와 같이 패드와 캘리퍼 간 접촉강성 및 패드와 디스크 간 접촉강성을 시스템 파라미터로 선정하여 운동방정식을 구축한다.
PPT Slide
Lager Image
Configuration of FE brake model with contact stiffness between pads and caliper
Fig. 2 와 같이 브레이크 시스템 단위에서 디스크는 Z축에 대해 일정한 속도( Ω )로 회전하고, 브레이크 패드 뒷면의 하중( No )에 의해 디스크와 패드간 마찰이 발생한다. 우선 디스크 및 패드의 접촉면에서의 변위벡터를 다음과 같이 표현한다.
PPT Slide
Lager Image
Coordinate systems and contact stiffness of friction material
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서
PPT Slide
Lager Image
,
PPT Slide
Lager Image
,
PPT Slide
Lager Image
는 디스크-패드 접촉면에서 각각 디스크의 반지름방향, 접선방향 및 수직방향의 디스크의 n 번째 모드형상 함수이고
PPT Slide
Lager Image
,
PPT Slide
Lager Image
,
PPT Slide
Lager Image
는 패드에 대한 모드형상 함수이다. 또한 e r , e θ , e z 는 각각 디스크의 반지름방향, 접선방향 및 수직방향의 단위방향벡터이다. 참고문헌 (7~9) 에 따르면 정지좌표계( x)의 패드와 회전좌표계(
PPT Slide
Lager Image
)의 디스크간 접촉응력을 좌표계를 통일하여 표현할 수 있다. 이때 디스크-패드 접촉면에서 접촉강성 kc 에 의한 접촉강성에너지가 다음과 같이 표현 된다.
PPT Slide
Lager Image
여기서 Ac 는 디스크와 패드간 접촉면을 의미한다.
또한 패드와 캘리퍼 간 접촉점에서의 접촉강성연성을 추가적으로 표현하고자 한다. 캘리퍼와 패드의 접촉점에서의 변위벡터를 다음과 같이 표현할 수 있다.
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서,
PPT Slide
Lager Image
,
PPT Slide
Lager Image
,
PPT Slide
Lager Image
는 접촉점 bi ( b 1 , b 2 , b 3 , b 4 )에서 각각 Z , X ,Y 방향의 캘리퍼의 n 번째 모드형상 함수를 의미한다. 또한
PPT Slide
Lager Image
는 접촉점 bi 에서 각각 Z , X ,Y 방향의 패드의 n 번째 모드형상 함수이다. e x , e y 는 각각 X ,Y 방향의 단위방향벡터이다. 캘리퍼와 패드간 접촉강성 kb 에 의해 발생하는 강성에너지는 다음과 같다.
PPT Slide
Lager Image
여기서 패드-디스크 접촉면에서 발생하는 마찰응력은 쿨롱의 법칙으로 표현할 수 있다.
PPT Slide
Lager Image
이때 V rel 는 디스크 및 패드 접촉점에서의 상대 미끄럼속도이다. 이러한 마찰응력벡터의 접선방향성분은 마찰연성(friction-coupling)을 발생시키며 선형화에 의해서 마찰연성행렬로 표현되고 방사방향성분은 시스템 감쇠역할을 하는 마찰댐핑(friction damping) 행렬로 전개가 된다 (11) . 이때 마찰연성행렬은 비대칭으로 시스템 불안정성을 직접 발생시킬 수 있으므로 스퀼해석에 매우 중요하다.
선형화된 형태로 운동방정식을 전개하면 다음과 같이 나타낼 수 있다.
PPT Slide
Lager Image
여기서 F d ( = F dT ) 는 마찰 댐핑(friction damping) 행렬이고, [ ω 2 ] , K c ( = K cT ) , K b ( = K bT ) , F c ( ≠ F cT ) 는 각각 단품 고유진동수, 패드-디스크 접촉강성, 패드-캘리퍼 접촉강성, 마찰연성 행렬을 의미한다.
이러한 시스템 브레이크 유한요소모델을 단품간 접촉강성으로 구성하여 시스템을 구성하면 시스템 접촉강성에 따른 스퀼 영향도를 폭넓게 연구할 수 있다. 특히 시스템 접촉강성을 특정값으로 지정하여 해석할 때 숨어있는 브레이크 스퀼 특성들을 파악하기가 대단히 어렵다. 하지만 이 연구에서 제안하는 방식으로 시스템 접촉강성에 따른 스퀼 민감도를 분석하면 전체적인 스퀼 특성 및 잠재적인 스퀼위험도 등을 파악하는데 큰 도움이 된다.
위 운동방정식 (8)은 선형미분방정식들의 집합으로 디스크 시스템의 동적 안정성(dynamic stability)을 다음과 같이 판별할 수 있다. 이러한 시스템의해를 다음과 같이 가정할 수 있다.
PPT Slide
Lager Image
위 식 (9)를 운동방정식 (8)에 대입하여 정리하면 다음과 같이 구할 수 있다.
PPT Slide
Lager Image
여기서 U 는 시스템 고유벡터로 0이 되면 무의미하므로 시스템 행렬의 역행렬이 존재하지 않을 조건은 다음과 같다.
PPT Slide
Lager Image
위 특성방정식(characteristic equation)의 해를 구하면 총 시스템 자유도 N tot ( = N b + 2 N p + N c )의 개수만큼의 해 λ 를 구할 수 있고 이를 시스템 고 유치(system eigenvalue)라고 한다. 이때 디스크 마찰진동시스템의 고유치는 다음과 같은 복소수 형태로 나타나게 되어 시스템방정식 (8)의 일반해를 묘사하게 된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서 n 번째 모드의 고유치의 실수부 Re( λn ) 가 양수이면 해당 모드에 의해 선형적으로 발산하므로 시스템의 선형안정성(linear stability)의 지표가 되며 브레이크 시스템의 마찰에 의한 소음발현 시점 (onset of squeal)을 나타낸다 (3~6) .
3. 해석결과
선형운동방정식 (8)의 고유치(eigenvalue)를 계산하기 위해서 계산 코드(solver)를 개발하여 계산하였다. 이때 선형운동방정식 (8)은 접촉부 강성 kc kb 에 직접 영향을 받고 있으며, 이를 두 접촉강성에 대해 안정성맵(stability map)를 구할 수 있다. 이를 통해 FE 스퀼해석에서 접촉강성부의 모델링에 대한 스퀼 민감도 및 중요성에 대해 분석하고자 한다. 패드 뒷면에 작용하는 정하중은 N o = 2000 N 이고 속도는 Ω = 6 rad / s 기준으로 해석한다. 또한 패드 마찰재의 기준 접촉강성 kc,nom 은 1.6×10 12 N ⋅m −3 이고 캘리퍼와 패드간 기준 접촉강성 kb,nom 은 3.0×10 7 N ⋅m −1 이며 마찰계수는 0.35 로 일정하다.
Fig. 3 은 마찰계수가 0.35로 상수일 때 패드와 캘리퍼 간 접촉강성 kb 에 따라 특정 구간에서 양의 실수부 ( Re( λ ) > 0 )가 발생하고 있음을 보여준다. 특히 2300 Hz 대역과 3200 Hz 대역에서 스퀼이 발생할 수 있음을 보여주며, Fig. 4 는 스퀼주파수에 상응하는 실수부의 크기를 축적하여 보여주고 있다. 이 는 이러한 모드들이 각각 패드 좌우방향(Y)의 변위에 민감하여 좌우방향의 접촉강성 kb 에 대한 주변 모드와의 마찰연성이 특정 구간에서 커지면서 모드 불안정성이 발생하고 있음을 의미한다.
PPT Slide
Lager Image
Frequency loci with respect to contact stiffness kb [%](= kb / kb,nom × 100[%]) ; dot represents Re (λ ) > 0
PPT Slide
Lager Image
Positive real parts of eigenvalues corresponding to system frequency of Fig. 3
Fig. 5 는 패드와 캘리퍼 간 접촉강성 kb 이 상기 두 스퀼주파수의 실수부에 미치는 영향도를 보여준다. 여기서 접촉강성 kb 이 변하면서 실수부가 급격히 변하는 것을 알 수 있으며, 이는 접촉강성이 스퀼해석 결과에 크게 영향을 주고 있음을 의미한다. 특히 3200 Hz 모드의 접촉강성 kb 에 대한 스퀼 성향 민감도가 매우 큼을 알 수 있다. 이는 3200 Hz 모드의 스퀼 특성을 분석할 때 패드와 캘리퍼 간 접촉부 모델링에 따라 큰 폭의 해석결과 편차를 발생시키게 된다.
PPT Slide
Lager Image
Real part loci with respect to contact stiffness kb / kb,nom × 100[%] up to 4 kHz
Fig. 6(a) kb 를 정해두고 패드와 디스크간 접촉강성 kc 에 따른 스퀼 발생 주파수를 보여준다. 여기서도 마찬가지로 2300 Hz 및 3200 Hz 구간에서 불안 정성이 발생하고 있음을 알 수 있다. Fig. 6(b) 는 스퀼주파수에 상응하는 양의 실수부를 보여주고 있다. 여기서 접촉강성 kc 이 상대적으로 적을 때 3200 Hz 구간의 불안정성이 사라지고 있다.
PPT Slide
Lager Image
Frequency loci with respect to contact stiffness kc[%] ( = kc / kc,nom × 100[%] ) at kb [%], = 80% (a) frequency loci, (b) positive real parts; dot represents Re (λ ) > 0
이렇듯 접촉강성의 변화에 따라 주요 스퀼모드의 2300 Hz 모드와 3200 Hz 모드의 스퀼성향이 변하고 있음을 알 수 있다. 이러한 접촉부 강성들의 변화에 대한 스퀼 영향도를 살펴보기 위해서 kb 를 100등분하고 kc 를 100등분하여 각각의 ( kb , kc )조합에 따른 스퀼해석을 총 10,000회 실시한다. 일반적인 유한요소 상용소프트웨어로 이렇게 많은 해석을 수행하는 것은 사실상 불가능하다. 하지만 이 연구에서 제안하는 모델로 계산하면 인텔 코어 i7 4GB RAM 기준으로 5시간 정도 소요되어 매우 효율적인 스퀼해석기법이다. Fig. 7 에서 보듯이 스퀼모드의 양의 실수부가 접촉부 강성에 매우 민감하게 변하고 있음을 알 수 있다. 특히 패드와 캘리퍼 간 접촉강성 kb 의 변화에 따라 특정 구간에서 양의 실수부가 급격히 증가함을 알 수 있다. 이는 시스템상으로 kb 가 스퀼에 민감한 값으로 설정되어 있다면 스퀼 특성이 악화될 수 있다. 또한 패드 및 디스크 간 접촉강성 kc 이 70 % 이상 증가함에 따라 3200 Hz 모드의 불안정성이 발생함을 알 수 있다. 해석 측면에서 보면 스퀼 해석 모델을 구축할 때 kb kc 를 실제 시스템과 일치시키지 못하면 약간의 편차가 매우 큰 결과 오차를 만들 수 있음을 암시하고 있다.
PPT Slide
Lager Image
Real parts with respect to two contact stiffness kb / kb,nom × 100[%] and kc / kc,nom × 100[%] , (a) 3D view, (b) top view
4. 결 론
이 연구에서는 디스크 브레이크 스퀼 해석을 유한요소법을 이용하여 수행할 때 접촉부 강성이 스퀼 해석 결과에 큰 영향을 줄 수 있음을 보여주었다. 이 연구에 사용된 디스크 브레이크 시스템에서는 2300 Hz 및 3200 Hz 근방의 특정 모드에서 불안정성이 발생하였다. 또한 패드와 디스크간 접촉강성과 패드와 캘리퍼 간 접촉강성에 의해서 스퀼 모드의 양의 실수부가 크게 변할 수 있음을 보여주었다. 특히 패드와 캘리퍼 간 접촉강성이 3200 Hz 스퀼 모드의 불안정성에 큰 영향을 주었다.
브레이크 스퀼 해석에서는 오직 스퀼 주파수만 중요하다. 따라서 스퀼 모드가 접촉강성에 큰 영향을 받고 있다는 점은 스퀼 해석 시 접촉강성 모델링이 매우 중요하다는 점을 보여주고 있다.
향후 연구에서는 이러한 해석결과를 바탕으로 브레이크 다이나모메터를 이용하여 접촉강성에 따른 스퀼 특성을 실험적으로 연구할 예정이다.
Acknowledgements
이 논문은 2014년도 정부(미래창조과학부)의 재원으로 한국연구재단의 지원을 받아 수행된 기초연구 사업임(No. NRF-2014R1A1A1002092).
BIO
Jaeyoung Kang is an associate professor in the Department of Mechanical and Automotive Engineering, Kongju National University. He received his Ph.D. degree in Mechanical Engineering from Purdue University in 2008. His research interests include friction noise, structural vibration and nonlinear dynamics.
References
Kinkaid N. M. , O’Reilly O. M. , Papadopoulos P. 2003 Automotive Disc Brake Squeal Journal of Sound and Vibration 267 105 - 166    DOI : 10.1016/S0022-460X(02)01573-0
Ouyang H. , Nack W. , Yuan Y. , Chen F. 2005 Numerical Analysis of Automotive Disc Brake Squeal: a Review International Journal of Vehicle Noise and Vibration 1 207 - 231    DOI : 10.1504/IJVNV.2005.007524
Nack W. , Joshi A. M. 1995 Friction Induced Vibration: Brake Moan SAE Warrendale, PA
Nack W. 2000 Brake Squeal Analysis by Finite Elements International Journal of Vehicle Design 23 263 - 275
Bajer A. , Belskyl V. , Zeng L. 2003 Combining a Nonlinear Static Analysis and Complex Eigenvalue Extraction in Brake Squeal Simulation SAE Warrendale, PA
Bajer A. , Belskyl V. , Kung S. 2004 The Influence of Friction-induced Damping and Nonlinear Effects on Brake Squeal Analysis SAE Warrendale, PA
Kang J. , Krousgrill C. M. , Sadeghi F. 2009 Wave Pattern Motion and Stick-slip Limit Cycle Oscillation of a Disc Brake Journal of Sound and Vibration 325 552 - 564    DOI : 10.1016/j.jsv.2009.03.030
Kang J. 2009 Squeal Analysis of Gyroscopic Disc Brake System based on Finite Element Method International Journal of Mechanical Science 51 284 - 294    DOI : 10.1016/j.ijmecsci.2009.02.003
Kang J. 2009 Linear Stability Analysis of a Rotating Disc Brake for Squeal Noise Transactions of the Korean Society for Noise and Vibration Engineering 19 (10) 1092 - 1098    DOI : 10.5050/KSNVN.2009.19.10.1092
Nam J. , Kang J. 2013 Brake Squeal Analysis with Respect to Caliper Contact Stiffness Transactions of the Korean Society for Noise and Vibration Engineering 23 (8) 717 - 724    DOI : 10.5050/KSNVE.2013.23.8.717
Kang J. , Krousgrill C. M. , Sadeghi F. 2009 Comprehensive Stability Analysis of Disc Brake:Gyroscopic, Negative Slope and Mode-Coupling Instability Journal of Sound and Vibration 324 387 - 407    DOI : 10.1016/j.jsv.2009.01.050
Kang J. 2012 Effect of Contact Area Variation on Squeal Propensity of Disc Modes Proc. IMechE Part C: J. Mechanical Engineering Science 226 2467 - 2473    DOI : 10.1177/0954406211435872