Transactions of the Korean Society for Noise and Vibration Engineering
[ Article ]
Transactions of the Korean Society for Noise and Vibration Engineering - Vol. 30, No. 6, pp.546-552
ISSN: 1598-2785 (Print) 2287-5476 (Online)
Print publication date 20 Dec 2020
Received 04 Sep 2020 Revised 22 Oct 2020 Accepted 16 Nov 2020
DOI: https://doi.org/10.5050/KSNVE.2020.30.6.546

하이브리드 모델링 기법과 베이지안 최적화를 통한 엔진 마운트 응답개선

강성호* ; 양성수* ; 김성준** ; 주형준** ; 강연준
Engine Mount Response Improvement Using Hybrid Modeling and Bayesian Optimization
Seong Ho Kang* ; Sung Soo Yang* ; Seoung Joon Kim** ; Hyung Jun Ju** ; Yeon June Kang
*Member, School of Mechanical & Aerospace Engineering, Seoul National University, Student
**Hyundai Motor Group, Powertrain Structure CAE Team, Researcher

Correspondence to: Member, School of Mechanical & Aerospace Engineering, Seoul National University, Professor E-mail : yeonjune@snu.ac.kr # A part of this paper was presented at the 49th International Congress and Exposition on Noise Control Engineering

‡ Recommended by Editor Jae Young Kang




© The Korean Society for Noise and Vibration Engineering

Abstract

This paper presents a method for predicting and improving the response of the engine mount using hybrid modeling technique and Bayesian optimization. The hybrid modeling technique based on the FRF based substructuring (FBS) theory can improve the accuracy of results and reduce the required time and cost compared to full-FE analysis using both experimental and finite element results. Using FBS theory, a transfer function of complex system can be changed to transfer function of several simple substructures. At this time, the hybrid modeling technique replaces the transfer function of a specific substructure with the result of FE analysis. Using this method, it is possible to predict the effect of the design change in the entire system by changing the CAD model design for the substructure applied to FE analysis. Furthermore, an optimal design is proposed by applying the Bayesian optimization technique.

Keywords:

Frequency Response Function, Substructure, FRF Based Substructuring, Hybrid Modeling Method, Engine Mount, Alternator-bracket, Bayesian Optimization

키워드:

주파수응답함수, 부분구조물, 주파수응답함수를 이용한 부분구조 합성법, 혼합 부분구조 합성법, 엔진마운트, 알터네이터 브라켓, 베이지안 최적화

1. 서 론

차량은 서스펜션, 차량 바디, 파워트레인 등 여러 부분구조물들의 조합으로 이루어진 복잡한 시스템이다. 이러한 복잡한 시스템을 FE 해석을 통해 동특성을 분석할 경우 해석에 소요되는 시간이 오래 걸리며 각 부분구조물들 사이의 결합특성에 의해 그 정확도가 떨어지게 된다. 이러한 문제를 개선하는데 지금까지 frequency response function based substructuring (FBS) 이론이 유용하게 사용되었다. FBS이론은 전체시스템의 전달함수를 간단한 각 부분구조물의 전달함수들 조합으로 표현 가능케 한다. FBS 이론을 기반으로 하는 하이브리드 모델링 기법은 시험을 통해 얻은 전달함수와 FE 해석을 통해 얻은 전달함수를 복합해 전체시스템의 전달함수를 예측하는 기법으로, 결합된 전체구조물을 full-FE 해석으로만 진행하는 것에 비해 실제 시험 값을 반영하기 때문에 결과의 정확도가 높으며 소요되는 시간을 단축시킬 수 있다. 또한 FE 해석이 적용되는 부분구조물의 설계변경 시 부분구조물에만 대한 CAD 모델링과 FE 해석 결과로 전체구조물의 전달함수를 예측할 수 있기 때문에 여러 가지 가능성에 대해 많은 데이터를 얻을 수 있으며 이 결과를 바탕으로 머신러닝, 딥러닝 등이 적용 가능하다.

이러한 연구는 차량 개발 초기 단계에 부분적인 설계변경에 대해서 실제 물리적 모델제작 없이 최적의 설계변경 방향을 제시 할 수 있다. 이 연구는 엔진 알터네이터의 진동에 의한 엔진 마운트에서의 응답 개선을 목표로 하이브리드 모델링 기법을 적용하였으며 이를 통해 얻어낸 마운트의 응답 값들을 베이지안 최적화를 이용해 최적의 설계변경 방향을 제시한다.


2. 하이브리드 모델링 기법

2.1 FBS 이론

FBS 기법은 여러 개의 부분 구조물이 결합된 복잡한 시스템의 전달함수를 결합되는 지점에서의 displacement compatibility와 force equilibrium 성질을 이용해 부분구조물들의 전달함수로 표현하는 이론이다. 식 (1)은 두 개의 부분구조물 A, B가 결합되어 전체구조물을 이룰 때의 FBS 이론 식을 나타낸 것이며(1~3), Fig. 1은 각 부분구조물 A, B가 결합되어 하나의 전체구조물을 이루는 모습을 보여준다.

HSoai(a)HSoac(a)HScai(a)HScac(a)HSoac(b)HSoai(b)HScac(b)HScai(b)HScbi(a)HScbc(a)HSobi(a)HSobc(a)HScbc(b)HScbi(b)HSobc(b)HSobi(b)=HAoai(a)HAoac(a)HAcai(a)HAcac(a)[0][0][0][0][0][0][0][0]HBcbc(b)HBcbi(b)HBobc(b)HBobi(b)- HAoacaHAcaca-HBcbcb-HBobcbC[HAcai(a)HAcac(a)-HBcbcb-HBcbi(b)(1) 
C=HAcac(a)+HBcbc(b)+[Kc]-1-1(2) 
Fig. 1

A two-substructure coupled system

여기서 [Kc]-1항은 부분구조물 A와 B가 강체결합을 한다고 가정했을 때 0으로 수렴한다.

2.2 알터네이터의 ODS와 Waterfall 분석

엔진을 이루는 여러 가지 부품 중 알터네이터의 경우 상대적으로 중량이 무거우며 엔진바디에 브라켓으로 연결되어있다. 따라서, 알터네이터에서 발생하는 진동은 엔진 마운트를 통해 결국 차량 바디로 전달된다. 즉, 알터네이터의 진동으로 인한 엔진 마운트의 응답을 줄이면 차량 바디로 전달되는 진동을 줄일 수 있다. 이 논문에서는 알터네이터의 진동 특성과 관심주파수대역 설정을 위해 엔진의 operational deflection shape (ODS) 분석을 진행하였다(4). 엔진 RPM을 700부터 4000까지 증가시킨 결과를 기반으로 알터네이터의 진동 크기를 확인하였으며 이러한 분석을 통해 400 Hz 부근 대역에서 알터네이터의 Z축 방향 움직임이 큰 것을 확인하였다. Fig. 2는 알터네이터의 operational analysis를 토대로 얻어진 waterfall 결과이다. ODS와 waterfall 분석을 통해 하이브리드 모델링 기법 적용 시 개선시켜야 할 엔진 마운트의 관심 주파수대역을 350 Hz ~ 500 Hz로 설정하였으며, Z축 방향의 움직임을 줄이기 위해 알터네이터 브라켓의 Z축 방향 구조강성을 증가시키기로 결정하였다. 알터네이터 브라켓의 경우 cantilever beam의 형상을 띄기 때문에 굽힘에 대한 저항력을 증가시켜 그 움직임을 억제할 수 있다.

Fig. 2

Waterfall result of alternator in Z-direction

2.3 하이브리드 모델링 기법의 신뢰성 판단

하이브리드 모델링 기법을 적용하기 전에 부분구조물의 전달함수로 계산된 전체구조물의 전달함수가 실제 전달함수를 적절하게 예측하는지와 Fig. 1에서 부분구조물 B를 FE 결과로 대체하였을 때 얻어진 하이브리드 모델링의 결과가 full-FE 해석으로만 분석하였을 때 얻어진 결과에 비해 개선되었는지 확인이 필요하다. 부분구조물 A의 경우 엔진바디와 마운트가 포함된 구조물이며, 부분구조물 B의 경우 알터네이터와 알터네이터 브라켓이 포함된 구조물을 의미한다. Fig. 3은 부분구조물의 시험 결과를 이용해 계산된 전체시스템의 전달함수와 전체구조물을 실제 측정한 결과를 비교한 것이다. Fig. 3을 보면 관심 주파수대역인 350 Hz ~ 500 Hz에서 예측한 전달함수 레벨과 피크가 유사하며 이를 통해 하이브리드 모델링 기법에 이용되는 식 (1)의 유효성을 확인 할 수 있다. 식 (1)의 유효성을 확인하였으면 부분구조물의 시험 값을 FE 해석 결과로 대체하여 하이브리드 모델링 기법으로 예측한 결과의 신뢰성을 파악하여야 한다. Fig. 4는 실제 측정결과, full-FE 해석결과, 하이브리드 모델링 기법으로 얻어진 [HS]o(a),i(b)를 비교한 것이다. 이를 통해, 하이브리드 모델링 기법 적용 시 full-FE 해석으로 얻은 [HS]o(a),i(b)에 비해 실제 측정결과를 더 잘 묘사한다. 이는 하이브리드 모델링의 경우 전체구조물의 전달함수를 계산할 때 실제 측정값이 사용되기 때문이다. Table 1을 보면 full-FE 해석과 하이브리드 모델링 기법 적용 시 소요되는 시간을 비교하였다. Full-FE 해석은 [HS]o(a),i(b)를 얻기 위해 총 11시간이 소요되나 하이브리드 모델링 기법을 이용할 경우 이의 0.76 %에 해당하는 5분 정도가 소요된다. 전체구조물의 전달함수를 예측하는데 짧은 시간이 소요된다는 장점은 부분적으로 FE 해석이 적용된 부분구조물의 설계변경에 대한 다양한 사례 연구를 할 수 있다. 이는 실제 물리적 모델 없이 많은 경우에 대해 [HS]o(a),i(b)를 예측 할 수 있으며 3장에 소개될 베이지안 최적화를 적용하여 알터네이터 브라켓의 구조강성 변경을 최적화 시키기 위한 초석이 된다.

Fig. 3

Comparison of experimental and calculated [HS]o(a),i(b)in Z-direction

Fig. 4

Comparison of experimental, hybrid modeling and full-FE analysis [HS]o(a),i(b)in Z-direction

Comparison of the time required to calculate [HS]o(a),i(b)

2.4 하이브리드 모델링 기법 적용

앞서 하이브리드 모델링 기법의 유용성을 확인하였으므로 FE 해석이 적용된 부분구조물 B에서 알터네이터 브라켓의 구조강성을 변경하여 전체구조물의 진동 응답 레벨이 어떻게 변하는지 예측하였다. 구조강성 변경 대상은 알터네이터 브라켓이며 ODS분석을 통해 400 Hz 대역 부근에서 Z축 방향의 움직임이 큰 것을 확인하였으므로 Z축 방향의 움직임을 억제하기 위해 가상의 rod를 덧대었다. Fig. 5는 기존 브라켓과 구조강성 변경한 브라켓을 나타낸 것이다. Rod는 총 9개가 사용되었으며 모두 Z축 방향의 움직임 억제를 고려하여 위치를 설정하였고, 지름은 7.5 mm, 10 mm, 12.5 mm, 재질은 기존 브라켓과 동일한 aluminum alloy 그리고 steel인 경우에 대해 분석하였다. Table 2는 각 재질의 물성치를 나타낸 것이다. Table 3는 aluminum alloy rod를 사용하였을 때 각 rod의 지름에 대해 관심 주파수 대역에서 전체 응답레벨 변화와 기존 브라켓에서의 전체 응답레벨을 비교한 것이며, Fig. 6은 aluminum alloy rod의 두께 변화에 따라 전체구조물의 [HS]o(a),i(b)를 예측한 것이다.

Fig. 5

Alternator bracket

Material property of aluminum and steel

Comparison of overall response level (applying aluminum alloy rod, 350 Hz to 450 Hz)

Fig. 6

Comparison of [HS]o(a),i(b)in Z-direction according the diameter of aluminum alloy rod

Table 4는 동일한 두께 조건으로 steel rod를사용하였을 때 전체 응답레벨을 비교한 것이며, Fig. 7은 steel rod의 두께 변화에 따라 전체구조물의 [HS]o(a),i(b)를 예측한 것이다. 이를 통해 브라켓 구조강성 변경 시 엔진 마운트의 응답 레벨을 개선 시킬 수 있음을 확인하였다.

Comparison of overall response level (applying steel rod, 350 Hz to 450 Hz)

Fig. 7

Comparison of [HS]o(a),i(b)in Z-direction according the diameter of steel rod


3. 베이지안 최적화

3.1 베이지안 최적화 이론

베이지안 최적화는 지능적 샘플링 과정을 통해 주어진 모델의 최적 파라미터를 탐색하는 알고리즘이다. 특히 목적함수를 모르거나 수학적으로 표현하기 어려운 black-box optimization 문제를 효율적으로 해결할 수 있으며, 최근 기계학습, 딥러닝 등 인공지능 알고리듬의 hyperparameter를 효과적으로 튜닝할 수 있어 크게 각광받고 있다(5). 베이지안 최적화의 핵심 아이디어는 최소한의 탐색으로 최적값을 찾는 것이다. 이를 위해 현재까지 관측된 정보를 기반으로 미지의 목적함수를 확률적으로 추정하는 surrogate model을 이용한다. 이 연구에서는 식 (3)과 같이 베이지안 선형 회귀 모델을 이용하였다.

y=wTϕx+ϵ,(3) 

여기서 xRd는 input 벡터, ϕ:RdRl는 feature map, wRl은 weight 벡터, ϵ은 N(0,σ2)을 따르는 노이즈이다. 이 연구에서는 식 (4)와 같이 feature map의 inner product가 width가 η인 가우시안 커널에 해당하도록 정의하였다(6).

ϕxTϕx'=exp-x-x'2η2.(4) 

학습된 surrogate model로부터 빠르게 최적점을 탐색하기 위해서는 최적점이 존재할 가능성이 높은 위치를 탐색하는 함수가 필요하며, 이는 acquisition function으로 정의된다. 이 연구에서는 가장 널리 사용되는 expected improvement(EI) 함수를 이용하였고 이는 식 (5)와 같이 정의된다(7).

EIx= Emaxfx-fx+, 0.(5) 

여기서 f(x+)는 현재까지 조사한 관측값 중 최대 함수값을 의미한다. EI를 이용한 acquisition function은 현재까지 탐색된 데이터를 기반으로 최적 포인트가 존재할 가능성이 높은 위치를 추정하며, 해당 포인트의 관측을 통해 surrogate model이 갱신된다. 이러한 반복적 학습의 과정을 거쳐 효과적으로 global optimum값을 찾아낼 수 있다.

3.2 베이지안 최적화를 이용한 알터네이터 브라켓 구조강성 최적화

앞서 9개의 rod를 이용해 브라켓의 구조강성을 변경할 경우 중량이나 비용적인 측면에서 다소 불리하다. 따라서, 적은 rod의 개수로도 효율적인 엔진 마운트 응답 저감을 위해 베이지안 최적화 기법을 활용하였다. 9개의 rod 중 임의의 4개 rod를 선정한다고 가정하면 126개의 독립적인 경우를 가지는 조합이 나온다. 이 논문에선 각 조합을 구분하기 위해 임의로 No. 1 ~ No. 126까지의 번호를 부여하였다. 이 경우의 수를 모두 시행하기엔 소요시간이 지나치게 많이 필요하기 때문에 베이지안 최적화를 활용해 최소의 조합만을 검증하여 최적의 조합을 찾아내었다. 최초 임의의 5개 조합을 input으로 하여 surrogate model을 학습한 뒤 acquisition function을 통해 다음 5개 조합을 추천받는다. 추천받은 5개 조합의 응답레벨을 다시 input으로 하여 같은 학습과정을 반복한다. Fig. 8은 최초 설정한 임의의 5개 조합(random input 1)에 대해서 5번의 학습을 통해 가장 낮은 응답 레벨을 갖는 조합을 확인한 것이며 필요한 검증 데이터의 수는 126개 중 약 20 % 정도인 25개의 데이터에 불과하다. Fig. 9은 input으로 사용될 임의의 조합을 다르게 설정하여(random input 2) 학습을 진행한 결과이다. Random input 1과 비교하였을 때 전체 조합수의 약 16 % 수준인 20개 조합의 검증만으로도 최적의 조합을 찾을 수 있었다. 104.7 dB를 가지는 최적의 조합 No. 12의 검증을 위해 126개의 조합을 모두 하이브리드 모델링 기법을 통해 확인하였고, 그 결과 가장 낮은 응답레벨을 갖는 최적의 조합임을 확인 할 수 있었다. Fig. 10은 최적의 조합인 No. 12의 rod 위치를 묘사한 것이다.

Fig. 8

Learning steps to find the optimal combination (using first random input)

Fig. 9

Learning steps to find the optimal combination (using second random input)

Fig. 10

Optimal modified bracket (No.12)


4. 결 론

하이브리드 모델링 기법을 적용함으로써 FE 해석으로만 전달함수를 계산하는 것에 비해 정확도가 높아지며 소요되는 시간을 줄일 수 있어 여러 가지 구조변경에 따른 전체구조물의 전달함수 변화 예측이 가능하다. 이는 많은 데이터를 빠르게 얻을 수 있다는 장점이 있어 머신러닝, 딥러닝과 같은 인공지능 학습과 연결시킬 수 있는 장점이 있다. 이 논문에서는 ODS 분석을 통해 엔진 마운트 응답에 영향을 미치는 구조물과 관심 주파수대역을 설정하였다. 이후 FBS이론을 기반으로 하는 하이브리드 모델링을 엔진구조에 적용하였고 알터네이터 브라켓의 CAD 모델 변경을 통해 여러 가지 경우에 대해 엔진 마운트의 응답을 예측하고 개선하였다. 이후 머신러닝의 일부인 베이지안 최적화를 통해 엔진 마운트의 응답을 개선시키기 위한 알터네이터 브라켓의 최적화된 구조강성을 제시했으며 이러한 결과를 통해 향후 설계변경시 어떤 부분의 강성 및 구조를 중점적으로 고려해야 하는지 빠르게 파악할 수 있다. 이러한 접근은 비단 엔진뿐만 아니라 차량 전체 시스템에 적용시킬 수 있으며, 차량의 NVH문제 해결 및 차량개발의 선행단계에서 개발에 필요한 소요시간과 비용을 크게 절감시킬 수 있을 것이다(8,9).

기 호 설 명

[HS] : 전체구조물의 전달함수
[HA] : 부분구조물 A의 전달함수
[HB] : 부분구조물 B의 전달함수
[Kc] : 부분구조물 사이의 결합부 강성
c(a) : 부분구조물 A의 결합지점
c(b) : 부분구조물 B의 결합지점
i(b) : 부분구조물 B의 가진점
ϵ : 베이지안 선형 회귀 모델의 노이즈
o(a) : 부분구조물 A의 응답점
ϕ(x) : 베이지안 선형 회귀 모델의 feature map
w : 베이지안 선형 회귀 모델의 weight 벡터

Acknowledgments

이 연구는 현대자동차 및 서울대학교 정밀기계설계공동연구소의 지원을 통해 진행되었으며 INTER-NOISE 2020에 발표한 내용을 토대로 작성되었다.

References

  • Tsai, J. S. and Chou, Y.-F., 1988, The Identification of Dynamic Characteristics of a Single Bolt Joint, Journal of Sound and Vibration, Vol. 125, No. 3, pp. 487~502. [https://doi.org/10.1016/0022-460X(88)90256-8]
  • Ren, Y. and Beards, C. F., 1995, On Substructure Synthesis with FRF Data, Journal of Sound and Vibration, Vol. 185, No. 5, pp. 845~866. [https://doi.org/10.1006/jsvi.1995.0421]
  • Liu, L., 2003, A Frequency Response Function-based Inverse Substructuring Approach for Analyzing Vehicle System NVH Response, Ph.D. Dissertation, The Univerity of Alabama, Alabama, USA.
  • Jin, S. and Pai, P. F., 2000, Locating Structural Defects Using Operational Deflection Shapes, Journal of Intelligent Material Systems and Structures, Vol. 11, No. 8, pp. 613~630. [https://doi.org/10.1106/F1LC-VJB5-HRCY-49NW]
  • Snoek, J., Larochelle, H. and Adams, R. P., 2012, Practical Bayesian Optimization of Machine Learning Algorithms, Proceedings of Advances in Neural Information Processing Systems 25, pp. 2951~2959.
  • Rahimi, A. and Recht, B., 2008, Random Features for Large-scale Kernel Machines, Proceedings of Advances in neural information processing systems 25, pp. 1177~1184.
  • Jones, D. R., Schonlau, M. and Welch, W. J., 1998, Efficient Global Optimization of Expensive Black-box Functions, Journal of Global optimization, Vol. 13, No. 4, pp. 455~492. [https://doi.org/10.1023/A:1008306431147]
  • Cuppens, K., Sas, P. and Hermans, L., 2001, Evaluation of the FRF Based Substructuring and Modal Synthesis Technique Applied to Vehicle FE Data, Proceedings of the International Seminar on Moal Analysis 3, pp. 1143~1150.
  • Liu, L. and Lim, T. C., 2003, Application of FRF-based Inverse Substructuring Analysis to Vehicle NVH Problems, SAE Transactions, Vol. 112, pp. 2006~2012. [https://doi.org/10.4271/2003-01-1607]

Seong Ho Kang received a B.S. degree from Ulsan National Institute of Science and Technology in 2017. He is currently a M.S. candidate in Acoustic and Vibration Lab. at Seoul National University. His research interests include Hybrid CAE/CAT modeling method.

Sung Soo Yang received a B.S. degree from Sung Kyun Kwan University in 2018. He is currently a Ph.D. candidate in Acoustic and Vibration Lab. at Seoul National University. His research interests include Acoustic materials.

Yeon June Kang received a B.S. and M.S. degree from Seoul National University in 1998, 1990. He then receive his Ph.D. degree from Purdue University in 1994. He is currently a professor at the department of mechanical engineering at Seoul National University in Seoul, South Korea, since 1997. His research include Acoustical Materials, Automotive NVH and Sound quality.

Fig. 1

Fig. 1
A two-substructure coupled system

Fig. 2

Fig. 2
Waterfall result of alternator in Z-direction

Fig. 3

Fig. 3
Comparison of experimental and calculated [HS]o(a),i(b)in Z-direction

Fig. 4

Fig. 4
Comparison of experimental, hybrid modeling and full-FE analysis [HS]o(a),i(b)in Z-direction

Fig. 5

Fig. 5
Alternator bracket

Fig. 6

Fig. 6
Comparison of [HS]o(a),i(b)in Z-direction according the diameter of aluminum alloy rod

Fig. 7

Fig. 7
Comparison of [HS]o(a),i(b)in Z-direction according the diameter of steel rod

Fig. 8

Fig. 8
Learning steps to find the optimal combination (using first random input)

Fig. 9

Fig. 9
Learning steps to find the optimal combination (using second random input)

Fig. 10

Fig. 10
Optimal modified bracket (No.12)

Table 1

Comparison of the time required to calculate [HS]o(a),i(b)

Analysis method Required time for calculate
[HS]o(a),i(b)
Full-FE 39 600 s
Hybrid modeling 300 s

Table 2

Material property of aluminum and steel

Material Young’s modulus [GPa] Poisson’s ratio Structuredamping coefficient Massdensity [kg/mm3]
Aluminum alloy 71 0.33 0.03 2.7e-06
Steel 210 0.29 0.02 7.9e-06

Table 3

Comparison of overall response level (applying aluminum alloy rod, 350 Hz to 450 Hz)

Diameter of rod [mm] Overall response level (original bracket) [dB] Overall response level (modified bracket) [dB]
7.5 119.4 114.7
10 111.4
12.5 108.8

Table 4

Comparison of overall response level (applying steel rod, 350 Hz to 450 Hz)

Diameter of rod [mm] Overall response level (original bracket) [dB] Overall response level (modified bracket) [dB]
7.5 119.4 110.2
10 108.0
12.5 105.7