The Korean Society Of Automotive Engineers
[ Article ]
Transactions of the Korean Society of Automotive Engineers - Vol. 28, No. 1, pp.43-52
ISSN: 1225-6382 (Print) 2234-0149 (Online)
Print publication date 01 Jan 2020
Received 19 Aug 2019 Revised 02 Nov 2019 Accepted 04 Nov 2019
DOI: https://doi.org/10.7467/KSAE.2020.28.1.043

무향 칼만 필터를 이용한 타이어 수직하중 변화 추정

김석원1) ; 정용우1) ; 김진성1) ; 이승희1) ; 정정주*, 2)
1)한양대학교 전기공학과
2)한양대학교 전기생체공학부
Estimation of Vertical Load Variation on a Tire Using Unscented Kalman Filter
Suk Won Kim1) ; Yong Woo Jeong1) ; Jin Sung Kim1) ; Seung-Hi Lee1) ; Chung Choo Chung*, 2)
1)Department of Electrical Engineering, Hanyang University, Seoul 04763, Korea
2)Department of Electrical and Biomedical Engineering, Hanyang University, Seoul 04763, Korea

Correspondence to: *E-mail: cchung@hanyang.ac.kr

Copyright Ⓒ 2020 KSAE / 170-06
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License(http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium provided the original work is properly cited.

Abstract

In this study, we present an estimator of vertical tire forces based on Unscented Kalman Filter(UKF) to improve the performance of the chassis control. When the vehicle is maneuvered to change lanes, the vertical load transfer caused by the centrifugal force can affect cornering stability and rollover. In the previous study, a Kalman Filter with a fixed gain, obtained from infinite horizon, has been conventionally used. However, this estimator does not fully compensate for the model’s uncertainty on slippery road conditions. To reflect the model’s uncertainty effectively, this paper is proposing a new estimator of vertical tire forces with UKF. Moreover, to reflect roll behavior on the low road friction coefficient accurŀately, we used roll rate information. The effectiveness and the performance of the proposed method were validated by using CarSim for vehicle dynamics.

Keywords:

Roll dynamis, Chassis control, Unscented kalman filter, Vertical tire forces, Lateral load transfer, Road friction coefficient

키워드:

롤 역학, 샤시 제어, 무향 칼만 필터, 타이어 수직력, 횡방향 하중이동, 도로 마찰 계수

1. 서 론

최근 지능형 자동차에 관한 연구가 활발히 진행되면서, 정확한 자세 제어를 위한 능동안전시스템(Active Safety System)의 중요성이 커지고 있다. 능동안전과 관련된 대표적인 섀시 제어로는 Electronic Stability Control(ESC), Torque Vectoring System(TVS), Continuous Damping Control(CDC), Active Front Steering(AFS)과 같은 것들이 있다.1) 이들은 차량의 선회안정성을 향상시키기 위해 개발되었기 때문에, 차량 모션과 관련된 정보가 필요하다.

자동차는 선회운동시 원심력에 의해 수직하중의 이동이 발생하게 되는데, 이는 차량의 선회 안정성 및 조정성에 중요한 영향을 미친다. 타이어 수직력은 차량의 내부센서(In-Vehicle-Sensor)에서 직접 측정이 되지 않기 때문에, 이를 계산 하기 위한 여러 방법들이 제안되었다.2-4) 먼저 전 차량 모델을(Full Car Model) 이용하는 방법이 논의되었지만, 해당 모델은 서스펜션 역학(Suspension Dyna-mics) 및 타이어의 비선형적인 특성 때문에 효과적인 제어전략의 수립이 어려웠다.5) 이러한 이유들 때문에, 간략화된 모델을 사용하여 타이어 수직력을 추정하는 방법들이 연구되었다. 윤장열 등6)과 조완기 등7)은 차량의 Inertial Measurement Unit(IMU) 센서와 3-DOF Vehicle Model을 바탕으로 한 Steady-State Kalman Filter(KF)를 설계하였다. 하지만 이러한 추정방법은 차량의 실시간 정보를 업데이트 하지 못하기 때문에, 다양한 주행 환경 변화를 정확히 반영하지 못한다. 윤성환 등8)은 Plane Roll Model의 파라미터를 변경하였으나, 실제 다양한 도로 환경에서의 차량 파라미터는 알기 어렵기 때문에 해당 방식은 그 정확성이 떨어진다.

본 연구에서는 타이어 수직력 추정기의 성능을 개선하기 위해 보다 향상된 추정 알고리즘을 제안한다. 먼저 차량의 IMU센서에서 측정되는 횡가속도 및 요레이트 정보로만으로는 미끄러운 도로 환경에서 발생되는 선회 특성을 정확히 반영하지 못하기 때문에, 롤 레이트를 측정값으로 추가하였다. 롤 레이트는 차량의 수직방향 고유 주파수 특성으로 구성되어 있으므로, 다양한 도로 노면 환경에서의 차량 응답 특성을 얻을 수 있다.9) KF의 경우는 Ricatti equation에서 계산된 고정된 게인값을 사용하기 때문에, 주행중인 차량 정보를 실시간으로 반영하지 못한다. 본 연구에서는 선형 모델을 사용하였기 때문에 Extended Kalman Filter(EKF)를 사용하지 않았다.6) 그러므로 Unscented Kalman Filter(UKF)에 기반한 추정기 알고리즘을 사용하여, 게인을 지속적으로 업데이트 한다. 이를 통해 미끄러운 주행 환경에서의 차량 정보를 실시간으로 반영할 수 있다.

제안된 방법의 검증을 위해 차량 동역학 Simulation 프로그램인 CarSim을 사용하였다. Dry Asphalt에서는 이전에 방법과 제시된 방법간의 성능차이가 유의미하지 않았다. 하지만 Ice Road 상에서는 Tire slip angle이 2도 이상 발생되면서 Lateral Force가 Saturation 영역을 넘어가게 된다. 이로 인해 비선형성이 증가하여 두 가지 방법의 성능차이가 크게 발생했다. 따라서 제안된 방법은 낮은 도로 마찰 계수 상황에서 유효한 방법임이 검증되었다.

본 논문은 총 5개의 장으로 이루어져 있으며, 그 구성은 다음과 같다. 2장, 3장에서는 차량의 타이어 수직력을 모델링하며 추정기 알고리즘을 설명한다. 4장에서는 제안된 추정기를 설계하여, 차량 동역학 해석프로그램인 CarSim으로 그 성능을 검증하였다.


2. 차량의 수직역학 모델링

본 장에서는 차량의 수직역학을 모델링하여 타이어 수직력을 추정한다. 먼저 2.1절에서는 타이어 수직력 추정기 구조도를 제시하며, 2.2절에서는 차량의 수직역학을 모델링하는 과정에 대해 설명한다.

2.1 타이어 수직력 추정기 구조도

본 연구에서 제시한 타이어 수직력 추정기 구조도는 Fig. 1과 같다. 해당 구조도에서는 입력으로 종방향 속도와 조향각이 주어졌을 때, 횡 가속도 와 요 레이트, 롤 레이트가 계산 및 측정되는 과정을 보여준다.6) 이후 해당 데이터들을 입력으로 하는 추정기(KF/UKF)를 통해 횡 가속도와 롤 각을 추정한다. 이렇게 추정된 값들을 이용하여 횡방향 하중이동을 구하고, 차량의 종 가속도로부터 종방향 하중이동을 구한다. 이를 바탕으로 각 바퀴의 타이어 수직력을 추정한다.

Fig. 1

Structure of vertical tire forces estimation

2.2 수직역학 모델링

이번절에서는 타이어 수직력 추정에 사용될 Bicycle Model, Simple Roll Model, Lateral/Longitudinal Load Transfer Model에 대해 설명한다.

2.2.1 Bicycle Model

차량의 횡방향 움직임에 대한 동역학을 묘사하기 위해, 차량이 고속으로 주행하는 상황을 가정한다. 이로 인해 선회반경은 차량의 축거 보다 훨씬 크고, 조향각의 범위가 작은 범위에 존재하면, 전후륜의 내측과 외측을 아래 Fig. 2와 같이 하나의 휠로 나타낼 수 있다.2,10)

Fig. 2

Bicycle model

또한 도로의 경사각이 없을 때 y 축에 대하여 뉴턴의 제2 법칙을 적용하면, 다음 식 (1)과 같이 표현된다.

may=Fyf+Fyr(1) 

선회하는 차량의 횡 가속도는 y 축에 의한 모션과 구심가속도에 의한 성분을 더하여, 식 (2)와 같이 나타낼 수 있다.

ay=y¨+VxΨ˙(2) 

이를 이용하여 차량의 횡방향 병진운동을 정리하면, 식 (3)과 같다.

my¨+VxΨ˙=Fyf+Fyr(3) 

또한 z 축에 대한 모멘트 평형식을 구하면 식 (4)와 같이 표현할 수 있다.

IzΨ˙=lfFyf-lrFyr(4) 

자동차가 고속으로 주행 중인 상황에서는, 차량의 진행 방향과 타이어의 진행방향이 일치하지 않게 된다. 이로 인해 타이어에는 슬립이 발생된다. 이때 슬립각이 작은 범위에 있다면 Fig. 3과 같이 횡력과 선형적인 관계를 갖고 있으며, 이는 식 (5)와 같이 표현될 수 있다.11,12)

Fig. 3

Relation between slip angle and lateral force

Fy=Cα α(5) 

작은 각 범위에서는 전후륜의 슬립각을 아래 식 (6)식 (7)로 근사화할 수 있다.

αf=y˙+lfΨ˙Vx(6) 
αr=y˙-lrΨ˙Vx(7) 

이상에서 식 (3)~(7)을 정리하면 아래 식 (8)과 같은 상태 공간 방정식의 형태로 나타낼 수 있다.

상태변수 x=y˙Ψ˙ 에 대하여 입력 u=δf 라 하면,

x˙=Ax+Bu

여기서

A=a11a12a21a22B=b11b21a11=-2Cαf+CαrmVx,a12=2lrCαr-lfCαfmVx-Vx,a21=-2Cαflf-CαrlrIzVx,a22=-2l2rCαr+l2fCαfIzVxb11=2Cαfm,b21=2lfCαfIz이다.(8) 
2.2.2 Simple Roll Model

앞의 Fig. 2에서 제시된 Bicycle Model은 롤이 발생하지 않는, 평면 움직임의 영역만을 다루었다. 그러나 실제 차량은 선회시 횡 가속도에 의해 아래 Fig. 4와 같이 롤 이 발생된다.6,13,14)

Fig. 4

Simple roll model

이때 차량은 요철이 없는 평지를 주행한다고 가정하면, 서스펜션의 역학을 2개의 스프링과 댐퍼로 단순화 할 수 있다. 이때 횡 가속도에 의해 발생된 롤 모멘트를 고려하여, 롤 축에 대한 토크평형식을 구하면 아래 식 (9)와 같이 표현할 수 있다. 이때 롤 각이 작다고 가정하면 (sinϕϕ, cosϕ≃1) 식 (10)과 같이 정리될 수 있다.15)

Ixϕ¨=-msayhrccosϕ+msghrcsinϕ-Krollϕ-Crollϕ˙(9) 
ϕ¨=-msayhrcIx+msghrc-KrollϕIx-Crollϕ˙Ix(10) 
2.2.3 Lateral Load Transfer Model

아래 Fig. 5에서는 좌/우 수직력과 횡 방향힘, 롤 영향을 고려하여 뉴턴의 제2법칙을 적용한다. 이를 통해 좌/우 수직하중의 변화를 아래 식 (11)로 표현할 수 있다.11,16,17)

Fig. 5

Latereal load transfer model

Fziw2-Fzow2-Fyhr=Ixϕ¨+Crollϕ˙+Krollϕ(11) 
ΔFz,LA=Fzi-Fzo=2Fyhr+2Krollϕw(12) 

이때 차량이 정상상태에서 선회한다고 가정하면, 식 (12)로 정리할 수 있다(ϕ¨=ϕ˙=ϕ). 즉 차량에서 발생되는 수직하중의 변화는 선회력에 의한 횡방향 하중 이동과 롤 역학에 의한 하중이동의 합으로 정의할 수 있다.

2.2.4 Longitudinal Load Transfer Model

차량이 일정한 속도로 주행하면서 선회를 하면, Fig. 6과 같이 종방향 가속도에 변화가 발생한다.7,18,19)

Fig. 6

Longitudinal load transfer model19)

이는 전후륜 수직하중의 변화를 일으키기 때문에, 타이어 수직력 추정시 고려될 필요가 있다. 차량의 무게중심과 롤센터 사이에서 발생되는 종방향 힘과 전후륜 수직력 사이의 토크 평형식을 고려하면, 식 (13)과 같이 표현할 수 있다.

FzrL-FzfL+Fxhrc=Ixϕ¨+Crollϕ˙+Krollϕ(13) 
ΔFz,LONG=Fzr-Fzf=-FxhrcL=-maxhrcL(14) 

이때 차량이 정상상태에서 주행하고 있다고 가정하면(ϕ¨=ϕ˙=ϕ), 이를 식 (14)와 같이 정리할 수 있다.4)

2.2.5 각 바퀴의 수직력 계산

자동차의 타이어 수직력을 계산하기 위해서는 정적인 상태에서의 하중을 기준으로 종/횡방향 수직하중 변화를 고려해야 한다. 여기서 정적인 상태의 하중은 차량을 저속으로 주행시킨 상황에서 계측한 값을 기준으로 삼았다. 앞의 2.2.3절과 2.2.4절에서 구한 횡방향 하중이동, 종방향 하중이동의 변화를 고려한 각 바퀴의 수직력은 Fig. 7과 같고, 다음 식 (15)와 같이 정리할 수 있다.7,20)

Fig. 7

Calculation of vertical tire force of each wheel20)

Fz,fL=12Ff+ΔFzf,LA-ΔFz,LONGFz,fR=12Ff-ΔFzf,LA-ΔFz,LONGFz,rL=12Fr+ΔFzr,LA+ΔFz,LONGFz,rR=12Fr-ΔFzr,LA+ΔFz,LONG(15) 

3. 타이어 수직력 추정기 설계

3장에서는 앞의 2장에서 사용된 모델들을 바탕으로 3-DOF 차량동역학을 모델링하고, 4장에서 설계될 2가지의 추정기 알고리즘을 설명한다.

3.1 3-DOF Vehicle Model

앞의 2.2절에서 사용된 Bicycle Model과 Simple Roll Model을 결합하면, 식 (16)과 같이 3-DOF Vehicle Model로 표현할 수 있다. 이때 상태변수 x=y˙Ψ˙ϕϕ˙T에 대하여 입력 u=δf 라 하면, 출력 y=ayΨ˙ϕ˙T 으로 쓸 수 있고 이를 상태 공간방정식으로 정리하면 식 (16)과 같다.6)

윤장열 등6)과 조완기 등7)은 측정값으로 횡 가속도, 요 레이트만을 사용한 반면, 본 연구에서는 롤 레이트를 측정값에 포함시켰다. 해당 값은 차량의 수직방향 응답을 해석하는데 있어서 매우 중요한 요소로 쓰인다.

x˙=a11a1200a21a22000001a41a42a43a44x+b11b210b41uy=c11c120001000001x+d1100u(16) 

여기서

a11=-2Cαf+CαrmVx,a12=2lrCαr-lfCαfmVx-Vxa21=-2Cαflf-CαrlrIzVx,a22=-2l2rCαr+l2fCαfIzVxa41=mshrcIx-2Cαf+CαrmVx,a42=mshrcIx2lrCαr-lfCαfmVxa43=-Kroll+mghrcIx,a44=-CrollIxb11=2Cαfm,b21=2lfCαfIz,b41=2CαfmshrcIxmc11=-2Cαf+CαrmVx,c12=2lrCαr-lfCαfmVx,d11=2Cαfm

3.2 Steady-State Kalman Filter

이경수는 타이어 수직력 추정시 Steady-State Kalman Filter (KF)를 사용하였다.6-8) 해당 추정기의 경우 시불변 시스템에서는 Error covariance P가 정상상태에 빠르게 도달할 수 있기 때문에, 정상상태 Covariance를 미리 계산하여 게인 K 값을 쉽게 결정할 수 있었다. 즉 식 (18)과 같이 Ricatti equation 을 통해 게인을 빠르게 계산할 수 있는 장점이 있었다. 이때 주요 알고리즘을 정리하면 다음과 같다.21-23) 식 (17)과 같이 어떤 Plant가 주어졌을 때,

xk+1=Φxk+Γuky~k=Hxk(17) 
P=ΦPΦT-ΦPHTHTHT+R-1HPΦT+ΥQΥT(18) 
L=PHTHTHT+R-1(19) 
x^k+1=Φx^k+Γuk+ΦLy~k-Hx^k(20) 

식 (18)은 해당하는 Error covariance P 를 Discrete-time-algebraic Ricatti equation 통해 구하는 과정을 나타내고 있다. 이를 통해 얻은 해 P식 (19)의 게인 L 을 계산하게 되고, 결론적으로 위의 식 (20)과 같이 상태 변수를 추정 할 수 있다. 이러한 방법은 적은 계산만으로 게인을 선정할 수 있다는 장점이 있는 반면에, 시스템이 Steady state 하지 않은 상황에서 Plant의 변화를 실시간으로 반영하지 못하는 단점이 있다.23) 따라서 추정성능을 개선시키기 위해서는 게인을 실시간으로 업데이트 하는 알고리즘이 필요하다.

3.3 Unscented Kalman Filter

앞절에서는 정상상태의 Ricatti equation을 이용한 게인을 사용하여 추정기를 설계하였다. 하지만 이러한 방법은 차량 모션이 비선형적으로 변하는 상황에서는 정확성이 떨어지는 단점이 있다. 본 연구에서는 이러한 추정방법의 단점을 보완하기 위해, Sample에 기반하여 Covariance 를 업데이트 하는 Unscented Kalman Filter(UKF) 알고리즘을 제시한다.21-24) UKF는 주어진 상태변수의 적은 샘플을 가지고 Weighted된 Sigma Points를 선정한다. 이렇게 선정된 값들을 바탕으로, 가우시안 분포를 Unscented Transform(UT)하여 Mean과 Covariance를 구한다.23) 임의의 변수 x 에 대해 차원이 M인 비선형 함수가 주어질 때, x는 아래식 (21)~식 (23)와 같이 Weighted 된 Sample Points 로 표현할 수 있다.

χ=xi02M,W=Wi02Mi=0,2M(21) 
x0=x¯+s0xi=x¯+si,xi+M=x¯-si+M for i=0,..M(22) 
W0=κM+κW0<i2M=κ2M+κ(23) 

Covariance P식 (24)와 같이 대각 행렬로 표현할 수 있다고 가정하면, 시그마 포인트 χ식 (25)에서와 같이 Scaling parameter κ에 의해 결정될 수 있다. 또한 λ는 시그마포인트의 Spread를 결정하며 작은 양수 값으로 알려져 있다.

P=diagσ12,.,σM2(24) 
S=si02MSST=M+κPsj=σjM+κ1+λ2 for j=1,.M(25) 

다음 식들은 모델에 기반한 Prediction 과정을 나타낸다. 식 (26)에서는 예측한 평균을 계산하며, 식 (27)에서 예측한 공분산을 얻는다.

x¯k+1=i=02MWix¯k+1i(26) 
P¯k+1=i=02MWix¯k+1i-x¯k+1x¯k+1i-x¯k+1T+Qk(27) 

각각의 예측된 포인트들이 모델을 통해 식 (28)과 같이 Instantiate 되어, 식 (29)에서는 예측한

Y¯k+1=y¯k+1i02M=hχ¯k+1(28) 
y¯k+1=i=02MWiy¯k+1i=i=02MWihx¯k+1i(29) 

출력값을 계산하게 된다. 이때의 Innovation covariance와 Cross-correlation 행렬은 각각 식 (30)식 (31)과 같다.

P¯xy=i=02MWix¯k+1i-x¯k+1y¯k+1i-y¯k+1T(30) 
P¯yy=i=02MWiy¯k+1i-y¯k+1y¯k+1i-y¯k+1T+Rk(31) 

모델로부터의 계산과정 이후 측정값과의 Correction 과정을 정리하면 아래와 같다. 이때의 상태 추정은 식 (32)와 같고, 게인과 A priori covariance는 각각 식 (33), (34)로 정리할 수 있다.23,24)

x^k=x¯k+Lkyk-hkx¯k(32) 
Lk=P¯xyP¯yy-1(33) 
Pk=P¯k-LkP¯yyLkT(34) 

4. 시뮬레이션 및 결과

4.1 시뮬레이션 조건

본 연구에서 제시한 타이어 수직력 추정기 알고리즘을 검증하기 위해 MATLAB/CarSim을 이용하였다. 그리고 다음과 같은 주행 시나리오를 상정하였다. 시나리오는 Single Lane Change 상황(시속 70 km/h)과 Double Lane Change(시속 50 km/h) 상황으로 분류하였다. 그리고 이를 다양한 노면 마찰계수 상황에서 시뮬레이션을 진행하였다. 이때 True 값은 CarSim에서 계측된 값을 의미한다.

E-Class SUV Parameter

4.1.1 시뮬레이션 조건 및 도로 환경

(1) Steering wheel angle

Fig. 8

Steering wheel angle (a) Single lane change, (b) Double lane change

(2) Vehicle trajectory

Fig. 9

Vehicle trajectory

4.2 시뮬레이션 결과

4.2.1 Single lane change (Mu=1.0)
4.2.2 Single Lane Change (Mu=0.2)

제시된 Single Lane Change에서 Fig. 10과 같이 도로 노면 마찰계수가 1.0인 상황에서는 Steady-State-KF 와 UKF의 성능은 유사하였다. 그러나 Fig. 11에서와 같이 도로 노면마찰계수가 0.2인 상황에서는 UKF의 추정값이 Ground Truth와 더 일치하는 경향을 보였다.

Fig. 10

Vertical tire forces estimation (a) Front right, (b) Front left, (c) Rear right, (d) Rear left [kN]

Fig. 11

Vertical tire forces estimation (a) Front right, (b) Front left, (c) Rear right, (d) Rear left [kN]

4.2.3 Double Lane Change(Mu=1.0)
4.2.4 Double Lane Change (Mu=0.2)

제시된 Dobule Lane Change에서 Fig. 12에서와 같이 도로 노면 마찰계수가 1.0인 상황에서는 Steady-State-KF와 UKF의 성능은 유사하였다. 그러나 Fig. 13에서와 같이 도로 노면마찰계수가 0.2인 상황에서는 UKF의 추정값이 Ground Truth와 더 일치하는 경향을 보였다. 즉 Ice Road 상에서는 Tire slip angle이 2도 이상 발생하면서, Lateral Force가 Saturation 영역을 넘어가게 된다. 이로 인해 비선형성이 증가하여 두 가지 방법의 성능차이가 크게 발생했다.

Fig. 12

Vertical tire forces estimation (a) Front right, (b) Front left, (c) Rear right, (d) Rear left [kN]

Fig. 13

Vertical tire forces estimation (a) Front right, (b) Front left, (c) Rear right, (d) Rear left [kN]

4.3 차량의 Roll rate 비교

Fig. 14는 CarSim에서 계측된 롤 레이트를 비교하였다. 마찰계수가 1인 상황에서는 롤 레이트가 사인파 형태로 움직인 반면, 마찰계수가 0.2인 경우에는 불규칙적으로 움직였다. 후자의 경우와 같이 비정상적인 선회상황에서는 타이어가 그립을 잃게 된다. 따라서 UKF에 기반한 추정기를 통해 롤 레이트의 변화된 특성을 반영하여, 타이어 수직력을 추정해야 한다.

Fig. 14

Measured roll rate (a) Single Lane Change (b) Double Lane Change

4.4 RMSE (Root Mean Square Error) 비교

아래 Table 2 ~ Table 5를 통해 제안된 Lane Change/ Double Lane Change 상황에서 UKF가 추정한 타이어 수직력이 KF가 추정한값보다 Ground Truth와 의 RMSE가 적은 것으로 나타났다.

RMSE of vertical tire forces : Lane change (Mu=1.0)

RMSE of vertical tire forces : Lane change (Mu=0.2)

RMSE of vertical tire forces : Double lane change (Mu=1.0)

RMSE of Vertical tire forces : Double lane change(Mu=0.2)


5. 결 론

본 연구에서는 UKF를 이용한 타이어 수직력 추정기 알고리즘을 제시하였다. 제안한 타이어 수직력 추정기를 검증하기위해 Single Lane Change, Double Lane Change 주행상황을 상정하였고, 각각 아스팔트 도로와 빙판길 도로 조건으로 분류하여 시뮬레이션을 진행하였다. 시뮬레이션 결과 도로 마찰계수가 낮은 상황에서는 비정상적인 선회상황이 발생하여 Covariance를 지속적으로 업데이트해야할 필요성이 있었다. 또한 제시된 추정알고리즘은 Mu값이 0.2~1.0 사이에서 유효한 결과를 보인 것으로 관찰되었다. 즉 미끄러운 도로에서 타이어 수직력 추정시 Covariance를 업데이트 하는 UKF 알고리즘을 적용하여 추정성능을 높일 필요가 있다. 이를 통해 정확한 타이어 수직력 정보가 필요한 여러 샤시 제어시스템에 적용한다면 그 성능이 향상될 것으로 기대된다.

Nomenclature

Cαf : cornering stiffness of front tires, N/rad
Cαr : cornering stiffness of rear tires, N/rad
y˙ : lateral velocity at c.g of vehicle, m/s
Ψ˙ : yaw rate, rad/s
lf : distance from c.g to front tires, m
lr : distance from c.g to rear tires, m
Vx : longitudinal velocity at c.g of vehicle, m/s
w : track width, m
Ff : static vertical force of front tire, N
Fr : static vertical force of rear tire, N
Fy, f : lateral force of the rear tire, N
Fy, f : lateral force of the front tire, N
ax : longitudinal acceleration, m/s
ay : lateral acceleration, m/s
IZ : yaw moment of inertia, Kg-m2
δf : front steering wheel angle, rad
m : total mass of vehicle, kg
mS : sprung mass, kg
ϕ : roll angle, rad
ϕ˙ : roll rate, rad/s
Kroll : roll stiffness, Nm/rad
Croll : roll damping coefficient, Nms/rad
hrc : height of the roll center from the c.g, m
Ix : inertia moment about its x-axis, Kg-m2

Acknowledgments

A part of this paper was presented at the KSAE 2018 Fall Conference and Exhibition

본 연구는 2019 년도 산업통상자원부 및 산업기술평가관리원(KEIT) 연구비 지원(20000293, 환경센서와 차량내부센서 융합기반 노면상태 검출 기술 개발) 과 산업통상자원부 및 중소벤처기업부의 연구비 지원(S2367759, 통합 이종카메라 시스템을 이용한 운전자지원시스템 개발) 에 의하여 연구되었음.

References

  • W. Cho, J. Choi, C. Kim, S. Choi and K. Yi, “Unified Chassis Control for the Improvement of Agility, Maneuverability, and Lateral Stability,” IEEE Transactions on Vehicular Technology, Vol.61, No.3, pp.1008-1020, 2012. [https://doi.org/10.1109/TVT.2012.2183152]
  • H. B. Pacejka, Tyre and Vehicle Dynamics, Butterworth-Heinemann, Oxford, 2002.
  • S. Song, M. C. K. Chun, J. Huissoon and S. L. Waslander, “Pneumatic Trail Based Slip Angle Observer with Dugoff Tire Model,” IEEE Intelligent Vehicles Symposium Proceedings, pp.1127-1132, 2014. [https://doi.org/10.1109/IVS.2014.6856538]
  • A. Rezaeian, R. Zarringhalam, S. Fallah, W. Melek, A. Khajepour, S. -Ken Chen, N. Moshchuck and B. Litkouhi, “Novel Tire Force Estimation Strategy for Real-time Implementation on Vehicle Applications,” IEEE Transactions on Vehicular Technology, Vol.64, No.6, pp.2231-2241, 2015. [https://doi.org/10.1109/TVT.2014.2345695]
  • J. Daknlallah, S. Glaser, S. Mammar and Y. Sebsadji, “Tire-road Forces Estimation Using Extended Kalman Filter and Sideslip Angle Evaluation,” American Control Conference, pp.4597-4602, 2008. [https://doi.org/10.1109/ACC.2008.4587220]
  • J. Yoon, D. Kim and K. Yi, “Design of a Rollover Index based Vehicle Stability Control Scheme,” Vehicle System Dynamics, Vol.45, No.5, pp.459-475, 2007. [https://doi.org/10.1080/00423110701245165]
  • W. Cho, J. Yoon, S. Yim, B. Koo and K. Yi, “Estimation of Tire Forces for Application to Vehicle Stability Control,” IEEE Transactions on Vehicular Technology, Vol.59, No.2, pp.638-649, 2010. [https://doi.org/10.1109/TVT.2009.2034268]
  • S. H. Yun, J. G. Yang, C. S. Kim and K. S. Huh, “A Study of Vertical Load Estimation for 4-wheel Torque Vectoring Control,” KSAE Fall Confenrence Proceedings, pp.211-213, 2016.
  • S. W. Kim, Y. W. Jeong, J. S. Kim, S. -H. Lee and C. C. Chung, “Vehicular Vertical Tire Foces Estimation Using Unscented Kalman Filter,” Asian Control Conference, pp.325-330, 2019.
  • R. Rajamani, Vehicle Dynamics and Control, Springer, New York, 2012. [https://doi.org/10.1007/978-1-4614-1433-9]
  • T. D. Gillespie, Fundamentals of Vehicle Dynamics, Society of Automotive Engineers, Warrendale, 1992. [https://doi.org/10.4271/R-114]
  • R. N. Jazar, Vehicle Dynamics: Theory and Applications, Springer, New York, 2008. [https://doi.org/10.1007/978-0-387-74244-1]
  • A. Hac, T. Brown and J. Martens, “Detection of Vehicle Rollover,” SAE 2004-01-1757, 2004. [https://doi.org/10.4271/2004-01-1757]
  • M. Doumiati, A. Charara, A. Victorino and D. Lechner, Vehicle Dynamics Estimation Using Kalman Filtering: Experimental Validation, Wiley, Hoboken, 2013. [https://doi.org/10.1002/9781118578988]
  • S. Bae, J. R. Yun, J. M. Lee and T. O. Tak, “An Effect of the Complexity in Vehicle Dynamics Models on the Analysis of Vehicle Dynamics Behaviors: Model Comparison and Validation,” Transactions of KSAE, Vol.8, No.6, pp.267-278, 2000.
  • D. C. Barton and J. D. Fieldhouse, Automotive Chassis Engineering, Springer, New York, 2018.
  • J. Reimpell, H. Stoll and J. W. Betzler, The Automotive Chassis: Engineering Principles, SAE International, Woburn, 2001.
  • M. Meywerk, Vehicle Dynamics, Wiley, West Sussex, 2015.
  • GMC Terrain SLE/SLT, https://www.gmc.com/suvs/terrain-smallsuv, , 2019.
  • Auto-Tech-Review, Technology, https://autotechreview.com/technology/tech-update/sensors-ecus-actuators-to-boost-suspension-performance, , 2019.
  • J. L. Crassidis and J. L. Junkins, Optimal Estimation of Dynamic Systems, Chapman and Hall, Florida, 2004. [https://doi.org/10.1201/9780203509128]
  • V. Kordic, Kalman Filter, InTech, Vukovar, 2010. [https://doi.org/10.5772/233]
  • S. W. Kim, Vertical Tire Forces Estimation and Active Steering Control Considering Roll Dynamics of Vehicle, M. S. Thesis, Hanyang University, Seoul, 2019.
  • S. W. Kim, Y. W. Jeong, J. S. Kim, S. -H. Lee and C. C. Chung, “Design of Roll Angle Estimator for Vehicles Using Unscented Kalman Filter,” KSAE Fall Conference Proceedings, pp.495-500, 2018.

Fig. 1

Fig. 1
Structure of vertical tire forces estimation

Fig. 2

Fig. 2
Bicycle model

Fig. 3

Fig. 3
Relation between slip angle and lateral force

Fig. 4

Fig. 4
Simple roll model

Fig. 5

Fig. 5
Latereal load transfer model

Fig. 6

Fig. 6
Longitudinal load transfer model19)

Fig. 7

Fig. 7
Calculation of vertical tire force of each wheel20)

Fig. 8

Fig. 8
Steering wheel angle (a) Single lane change, (b) Double lane change

Fig. 9

Fig. 9
Vehicle trajectory

Fig. 10

Fig. 10
Vertical tire forces estimation (a) Front right, (b) Front left, (c) Rear right, (d) Rear left [kN]

Fig. 11

Fig. 11
Vertical tire forces estimation (a) Front right, (b) Front left, (c) Rear right, (d) Rear left [kN]

Fig. 12

Fig. 12
Vertical tire forces estimation (a) Front right, (b) Front left, (c) Rear right, (d) Rear left [kN]

Fig. 13

Fig. 13
Vertical tire forces estimation (a) Front right, (b) Front left, (c) Rear right, (d) Rear left [kN]

Fig. 14

Fig. 14
Measured roll rate (a) Single Lane Change (b) Double Lane Change

Table 1

E-Class SUV Parameter

Parameter Description Value Unit
Caf Cornering stiffness of front tires 129487 N/rad
Car Cornering stiffness of rear tires 129487 N/rad
lf Distance from c.g to front tires 1.18 m
lr Distance from c.g to rear tires 1.77 m
Iz Yaw moment of inertia of vehicle 2488 kg-m2
m Total mass of vehicle 1862 kg
Kroll Roll stiffness 135206 Nm/rad
Croll Roll damping coefficient 20651 Nms/rad
ms Sprung mass 1592 kg
hκ Height of the roll center from the c.g 0.385 m
Ix Moment of inertia about its roll axis 614 kg-m2
w Track width 1.575 m
hr Height of the roll center 0.334 m
h Height of the c.g 0.719 m

Table 2

RMSE of vertical tire forces : Lane change (Mu=1.0)

Force FR FL RR RL
Filter KF UKF KF UKF KF UKF KF UKF
RMSE 4.8 2.4 4.8 2.3 4.8 1.9 4.9 1.9

Table 3

RMSE of vertical tire forces : Lane change (Mu=0.2)

Force FR FL RR RL
Filter KF UKF KF UKF KF UKF KF UKF
RMSE 3.8 2.5 3.8 2.4 4.5 2.6 4.4 2.5

Table 4

RMSE of vertical tire forces : Double lane change (Mu=1.0)

Force FR FL RR RL
Filter KF UKF KF UKF KF UKF KF UKF
RMSE 5.4 1.7 5.4 1.5 5.4 2.2 5.4 2.2

Table 5

RMSE of Vertical tire forces : Double lane change(Mu=0.2)

Force FR FL RR RL
Filter KF UKF KF UKF KF UKF KF UKF
Unit:[kN]
RMSE 3.7 2.2 3.7 2.1 4.8 2.9 4.8 2.8