Next Article in Journal
End-to-End Encrypted Message Distribution System for the Internet of Things Based on Conditional Proxy Re-Encryption
Next Article in Special Issue
Research on Electric Oil–Pneumatic Active Suspension Based on Fractional-Order PID Position Control
Previous Article in Journal
The Development of a Measuring System for Intraoral SpO2
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Adaptive Unscented Kalman Filter for the Estimation of the Vehicle Velocity Components, Slip Angles, and Slip Ratios in Extreme Driving Manoeuvres

1
Centre for Automotive Engineering, University of Surrey, Guildford GU2 7XH, UK
2
McLaren Automotive, Woking GU21 4YH, UK
3
Tenneco Automotive, 3800 Sint-Truiden, Belgium
4
Independent Researcher, 21100 Varese, Italy
5
Departement of Mechanical and Aerospace Engineering, Politecnico di Torino, 10129 Turin, Italy
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(2), 436; https://doi.org/10.3390/s24020436
Submission received: 20 October 2023 / Revised: 2 December 2023 / Accepted: 26 December 2023 / Published: 10 January 2024
(This article belongs to the Special Issue Vehicle Sensing and Dynamic Control)

Abstract

:
This paper presents a novel unscented Kalman filter (UKF) implementation with adaptive covariance matrices (ACMs), to accurately estimate the longitudinal and lateral components of vehicle velocity, and thus the sideslip angle, tire slip angles, and tire slip ratios, also in extreme driving conditions, including tyre–road friction variations. The adaptation strategies are implemented on both the process noise and measurement noise covariances. The resulting UKF ACM is compared against a well-tuned baseline UKF with fixed covariances. Experimental test results in high tyre–road friction conditions show the good performance of both filters, with only a very marginal benefit of the ACM version. However, the simulated extreme tests in variable and low-friction conditions highlight the superior performance and robustness provided by the adaptation mechanism.

1. Introduction

The latest generation of vehicle dynamics controllers, including stability, traction, active steering, and active suspension controllers, as well as automated driving systems, demands increasing amounts of data from the vehicle and its surroundings [1,2,3,4,5]. However, specific vehicle variables, such as the sideslip angle, tyre slip ratios and slip angles, vehicle speed, and tyre–road friction level, cannot be robustly and cost-effectively measured for mass production implementations, and therefore must be estimated. While the data measurements and estimations used since the 1990s for vehicle stability control (VSC) [6] are sufficient for rather simple rule-based algorithms, operating only in emergency scenarios, the progress in high-performance vehicle controllers, implying more frequent and smooth interventions, has increased the demand for more accurate information on the vehicle variables [7,8]. For example, reliable vehicle speed and slip ratio estimations are crucial for effective wheel slip control through actuators capable of continuously and seamlessly modulating the wheel torque levels [9].
In this context, machine learning methods have gained significant attention [10,11,12,13,14,15,16]. These algorithms offer enhanced adaptability to unpredictable situations as they are not based on predefined models, and tend to exhibit increased robustness in scenarios where the accurate modelling of the system dynamics is challenging. Although the resulting agents typically demand less computational power than model-based approaches, their practical implementation requires a very substantial amount of training data [16,17,18]. As mentioned in [18], ‘constructing a complete dataset that can cover all driving conditions is extremely difficult if not impossible’. As a consequence, although machine learning approaches are likely to be more diffusely implemented for state estimation in the near future, Kalman filters are still widely adopted tools to achieve the real-time estimation of vehicle dynamics variables, starting from the measurements from multiple on-board sensors [19].
Several algorithmic solutions have been proposed to account for the vehicle system nonlinearities in the internal model of the filter. This can be typically achieved through extended Kalman filters (EKFs), which imply a linearisation at each time step, or unscented Kalman filters (UKFs), which directly use nonlinear models. Comparative studies of EKF and UKF performance for the same vehicle model, e.g., see [16,20,21,22,23,24], highlight the superiority of UKFs, towards which there has been a progressive shift [17].
In terms of vehicle models embedded in the filters, the literature has assessed a wide range of options. For instance, within a complex estimation architecture for automated vehicles, reference [25] implemented a Kalman filter for sideslip angle estimation, using a 2-degree-of-freedom (2-DoF) linear bicycle model, considering the yaw and sideslip dynamics. This setup performed well during normal driving; however, it lost accuracy at high lateral accelerations or in critical transients. A possible solution to such limitation is the adoption of a double-track vehicle model (DTM), which has become rather widespread [26]. For example, the DTMs within the filters in [21,27,28] include the lateral tyre force nonlinearities, which are modelled through the Pacejka magic formula or the Dugoff formulation. The EKF implementation in [29] and the UKFs in [30,31,32] augmented the DTM by including the longitudinal vehicle dynamics and wheel dynamics, which resulted in 7-DoF vehicle models, enabling tyre slip ratio estimation.
In general, there are limited studies assessing Kalman filter performance on low-friction surfaces. References [29,32] are rare exceptions, but their setups are constrained by fixed settings of the filter parametrisation. Additionally, the tests in [29,32] omitted high wheel slip cases, such as wheel spinning and locking, analysing only the sideslip angle estimation.
In most vehicle dynamics filters, vehicle and tyre parameters are assumed to remain constant over time. However, some implementations deal with parametric uncertainties by embedding parallel estimators for the main vehicle parameters, or by concurrently estimating vehicle states and parameters. For example, the Kalman filter in Wenzel et al. [33] was supported by a second Kalman filter, which enabled the update of the vehicle mass, yaw mass moment of inertia, and longitudinal position of the centre of gravity. The algorithms in [30,31,32,34,35] varied the peak tyre–road friction coefficient and/or the cornering stiffness during their operation.
Similarly to the vehicle and tyre parameters, the EKF/UKF covariance matrices usually remain constant over time, and can be tuned through Bayesian optimisations [36] to find desirable trade-offs for a catalogue of manoeuvres. However, the effectiveness of the resulting filter is limited to specific conditions, and significant issues can occur during real-world operation. An effective approach to improve predictability and handle situations where the uncertainty characteristics of the system change over time is to employ adaptive covariance matrices [37,38,39,40,41,42,43,44]. The adaptation approaches can either be based on theoretical algorithms, e.g., using fading factors to give higher weight to the last filter iterations, or on covariance definition rules, e.g., imposed through fuzzy logic algorithms or other heuristics benefitting from the understanding of vehicle dynamics.
References [37,38,39,40] present adaptation strategies for the measurement noise covariances. Specifically, an adaptation scheme limited to the wheel speed measurement noise covariances was used in the EKF in [37], to target high wheel slip conditions, based on empirical rules. In such conditions, although the wheel speed and longitudinal acceleration measurements processed by the typical vehicle stability control unit are accurate, reliable, and without outliers, it is usually better to vary the weight associated with the angular wheel speed and longitudinal acceleration measurements, to prioritise the latter for vehicle speed estimation. Among the theory-based algorithms, references [38,39] used the moving window estimation method, while reference [40] relied on a maximum a posteriori algorithm. However, the assessment was limited to high-friction conditions, and no adaptation was applied to the process noise covariances.
By employing a comparable method to [40], reference [41] adapted the process and measurement noise covariances. Although this algorithm can offer stability in scenarios with noisy or incomplete data, its complexity could become critical when applied to filters characterised by vehicle models with higher number of degrees of freedom (DoFs). Reference [42] discusses a normalized innovation squared (NIS) algorithm to adaptively change both process and measurement noise covariances, by comparing the difference between predicted and observed measurements to the predicted covariance of measurement noise. However, alongside its computational complexity, the NIS approach requires precise tuning, and lacks robustness when encountering model inaccuracies. Reference [43] proposes an exponential decay adaptive unscented filter, in which the exponential decay factor helps assign different weights or decay rates to older and newer measurements, allowing the filter to adaptively weigh recent information more heavily than older data. Nevertheless, the algorithm lacks flexibility in accommodating rapid variations in state variables, resulting in compromised estimation accuracy, as shown for a sinusoidal steering manoeuvre. The authors conclude that “…however, the average weighting method and the general exponential weighting method are not flexible enough. As a result, when the state variable changes rapidly, the noise cannot be accurately estimated, and the state variable cannot be evaluated more accurately…. The proposed method has some improvement effect compared with the original method, but the effect is not significant enough in some conditions, which is a defect of the proposed method”. In summary, the adaptation solutions in [41,42,43] primarily relied on the measurement noise level, and were not particularly effective for swift variations of states or tyre–road friction parameters, or for incorporating physics-based heuristics. Additionally, the proposed filters used simplified models, neglecting the wheel dynamics, which are of the essence for tyre slip dynamics estimation in extreme conditions. Finally, the evaluation was limited to simulations with constant or slowly varying torque demands, without significant wheel slip events.
Interestingly, van Aalst et al. [44] combined parameter and covariance adaptations. In fact, the cornering stiffness values are adapted with a random walk model, while the filter covariances are adjusted according to the operating conditions of the vehicle, by using vehicle-dynamics-derived heuristics. The resulting scheme can thus effectively account for uncertainties and nonlinearities across a wide range of excitation inputs. However, the algorithm, based on an EKF approach, is not designed for longitudinal speed and slip ratio estimation, and is assessed only in high-friction conditions.
A recent trend is the use of combinations of state estimators, e.g., through so-called fusion architectures that embed both kinematic and dynamic vehicle models. In fact, a kinematic model is not affected by uncertain vehicle parameters or tyre behaviour, and can be particularly effective in transient conditions and for large excitation magnitudes. On the contrary, during normal vehicle operation, dynamic model-based filters provide better performance. For example, the solution in [45] combined a consensus Kalman filter, using a dynamic vehicle model, and a kinematic model-based Kalman filter. The architecture was applied to an automated vehicle, with the two filters receiving heading and velocity errors determined through a global navigation satellite system (GNSS) and a reduced inertial navigation system (R-INS). A weighting scheme based on the lateral excitation level fused the two outputs, bringing improved estimation accuracy during critical driving conditions. In the architecture in [46], an algorithm adaptively determined the cornering stiffness of a single-track model for sideslip angle estimation during quasi-steady-state cornering, through a proportional integral (PI) controller. Additionally, a kinematic approach in the form of an EKF was employed for transients and large excitation magnitudes. Similarly to [45], a weighting strategy combined the filter outputs. Whilst the resulting fusion architectures were assessed through a wide range of experimental data, neither paper embedded tests on low- μ surfaces, or incorporated the wheel dynamics, i.e., the tyre slip ratios could not be accurately estimated during wheel slip control transients. Furthermore, these approaches are more complex to tune than a single estimator, and thus, instead of fusing two estimators to maximise performance, a covariance adaptation strategy on a single estimator may be more practical.
In summary, the available literature lacks:
  • Adaptive dynamic model-based Kalman filters that concurrently vary the process and measurement noise covariances through vehicle-dynamics-derived heuristics. These should be defined as a function of relevant error variables, based on the difference between the measured and estimated outputs, with the specific scope of enhancing performance in highly dynamic conditions, including significant longitudinal tyre slip variations induced by wheel torque or tyre–road friction transients. These scenarios, unlike the current adaptive implementations, require internal models accounting for the wheel dynamics.
  • The assessment of the sideslip angle, velocity and tyre slip estimation performance benefit of adaptive Kalman filters in both high- and low-friction conditions, including μ -jump tests, and manoeuvres with very high levels of longitudinal and lateral slip.
This paper targets the identified gap, with the following contributions:
  • A UKF implementation with vehicle-dynamics-based adaptive formulations for the following: (a) the wheel speed measurement covariances, based on variables that are robustly representative of the longitudinal tyre slip condition; (b) the tyre–road friction coefficient process noise covariance, based on error variables depending on the longitudinal and lateral accelerations; and (c) the process noise covariances of the yaw rate and the longitudinal and lateral velocity components, based on the estimation errors with respect to the available relevant measurements.
  • The experimental validation of the UKF with adaptive covariance matrices, referred to as UKF ACM, along extreme high-friction manoeuvres, including significant longitudinal and lateral accelerations.
  • The validation of UKF ACM through a high-fidelity and experimentally validated model, in conditions with very low tyre–road friction, including μ -jumps.
  • The comparison of UKF ACM with a baseline UKF with fixed covariance values that are well calibrated.
The remainder is organised as follows: Section 2 introduces the case study vehicle and simulation models; Section 3 presents the proposed UKF implementations; Section 4 deals with the case study manoeuvres and relevant key performance indicators; Section 5 reports and critically analyses the experimental and simulation results; finally, Section 6 draws the main conclusions.

2. Case Study Vehicle and Associated Models

2.1. Reference Vehicle

The case study was an Audi e-tron sport utility vehicle prototype, used as one of the demonstrators of the European Horizon 2020 EVC1000 project. The vehicle has two on-board electric motors, one per axle, connected to the wheels via a single-speed transmission and mechanical differential. The main vehicle parameters are in Table 1. Figure 1 shows a photograph of the vehicle, with a Kistler Correvit S-350 sensor mounted on the rear left door. The sensor provides accurate measurements of the longitudinal and lateral components of vehicle velocity, and thus sideslip angle, which are recorded for validating the values estimated by the filters. The vehicle also includes the typical factory-installed measuring equipment, e.g., an inertial measurement unit (IMU), as well as steering wheel angle, individual wheel speed, accelerator pedal position, and tandem master cylinder pressure sensors. The considered filters only used the conventional production vehicle measurements as inputs.

2.2. High-Fidelity Vehicle Simulation Model

A high-fidelity vehicle simulation model was implemented with the CarMaker tool by IPG, for the simulation-based assessment of the UKF performance. The model includes the six DoFs of the sprung mass, the four DoFs associated with the unsprung mass motions, and the four wheel rotations, and considers the associated dynamic couplings. The suspension elasto-kinematic characteristics and nonlinearities associated with the shock absorbers and bump stops were also accounted for. Tyre behaviour was emulated with the Pacejka magic formula [47] including relaxation effects. The model parameters were directly provided by the involved car manufacturer, which was part of the EVC1000 initiative.

2.3. Internal Model of the Filters

The internal model embedded in the UKFs was specifically designed to achieve its intended objectives, i.e., it includes only the essential DoFs for estimating the sideslip angle, vehicle speed, slip angles, and slip ratios, in order to achieve computational efficiency and real-time implementability. Therefore, the selected formulation is a nonlinear 7-DoF model, considering the longitudinal, lateral, and yaw dynamics of the vehicle, in addition to the rotational dynamics of each of the four wheels. The schematic in Figure 2 highlights the key parameters and variables, together with their sign conventions, according to the ISO standard [48]. The internal model was developed by Mazzilli et al. [49]; hence, its detailed derivation will not be included here.
The primary governing equations of motion are the following:
  • Longitudinal force balance
    a x = v ˙ x v y r = 1 m t o t i = 1 2 F x , i c o s δ i F y , i s i n δ i + i = 3 4 F x , i F d r a g
  • Lateral force balance
    a y = v ˙ y + v x r = 1 m t o t i = 1 2 F x , i s i n δ i + F y , i c o s δ i + i = 3 4 F y , i
  • Yaw moment balance
    r ˙ = 1 I z [ a i = 1 2 ( F x , i s i n δ i + F y , i c o s δ i ) + ( t f 2 y C O G ) ( F y , 1 s i n δ 1 F x , 1 c o s δ 1 )            ( t f 2 + y C O G ) ( F y , 2 s i n δ 2 F x , 2 c o s δ 2 ) ( t r 2 y C O G ) F x , 3        + ( t r 2 + y C O G ) F x , 4 b i = 3 4 F y , i + i = 1 4 M z , i ]
  • Wheel moment balance
    I w e q , i ω ˙ i = M E M , i M B , i F x . i R l , i M y , i
    where the increasing values of the integer i , with i = 1,…, 4, used as a subscript, refer to the front left, front right, rear left, and rear right vehicle corners; a x and a y are the longitudinal and lateral accelerations; v x and v y are the longitudinal and lateral components of the vehicle velocity; F x , i and F y , i are the longitudinal and lateral tyre forces; δ i is the steering angle of the i -th wheel; r is the yaw rate; F d r a g is the aerodynamic drag force; y C O G indicates the lateral position of the centre of gravity with respect to the plane of symmetry of the vehicle; M z , i is the tyre self-aligning moment; I w e q , i is the equivalent mass moment of inertia of the wheel; ω ˙ i is the angular wheel acceleration; M E M , i is the electric motor torque level referred to the i -th corner, computed from the measured motor current, the gearbox and final drive ratios, as well as the respective efficiencies; M B , i indicates the braking torque at the individual corner, which is estimated from the measured tandem master cylinder pressure; R l , i is the laden tyre radius; and M y , i is the rolling resistance moment.
The model computes the longitudinal and lateral load transfers through steady-state equations based on a x and a y . The vehicle sideslip angle at the centre of gravity, β C O G , which is one of the key variables to be estimated, is given by:
β C O G = tan 1 v y v x
Similarly, the classical kinematic equations are used for the tyre slip ratios, σ x , i , and slip angles, α i , which are modified according to the analyses in [50,51], to ensure stability and realism also during low-speed operation. The longitudinal and lateral tyre forces are obtained through a simplified version of the magic formula (version 5.2, see [47]).
In summary, the formulation of the internal model of the filter is totally independent from the CarMaker model in Section 2.2. The very substantial differences in key aspects of the models (e.g., number of DoFs, load transfer formulations, magic formula version) also imply different model parametrisations. To make the UKF performance assessment more realistic and conservative, the model calibrations, based on the available vehicle data, did not involve any further optimisation to match the collected experimental results. The diversity in model arrangements and tuning setups makes the following simulation analyses realistic and indicative of the operational robustness of the UKFs.

2.4. Experimental Validation of the Models

The CarMaker and internal models have been experimentally validated through series of quasi-steady-state and transient tests in high-friction conditions, carried out during the EVC1000 project. Further experimental data was gathered on the very low-friction ice tracks at the Ebbenjarka testing facility in Sweden, which was characterised by variability and unpredictability of the surface conditions.
Figure 3 reports a sample of the validation results, referring to (a) a quasi-steady-state 40 m radius skid pad manoeuvre on dry tarmac; (b) a sinusoidal steering manoeuvre from an initial speed of ~80 kmh−1 and constant accelerator pedal position; and (c) a sinusoidal steering manoeuvre at ~40 kmh−1 on an icy surface.
The three sets of data—the experimental results and the simulation results from the two models—highlight that the CarMaker model provided a consistently good match for experiments (a–c), and therefore could be considered a reliable tool for UKF performance assessment across a wide range of conditions, including very low tyre–road friction. On the contrary, although performing well in the high- μ tests, the UKF models severely underestimated the sideslip angle magnitude in the low-friction test (c), which highlights the necessity of the UKF ACM.

3. Adaptive Unscented Kalman Filter Architecture

3.1. Filter Architecture and Strategy

Figure 4 is an overview schematic of the filter architecture, which was implemented and assessed in real-time on a dSPACE MicroAutoBox II unit (900 MHz, 16 Mb flash memory), with a 10 ms sampling time.
The system state, input, and output (corresponding to the available measurements) vectors, respectively, x , u , and y , are:
x =   v x    v y     r     ω i     μ ¯ m a x u =   δ i    M E M , i     M B , i    y =   a x    a y     r     ω i  
As per [49], the average peak tyre–road friction coefficient ( μ ¯ m a x ) of the four tyres is one of the states, and is estimated via a random walk model in the time update step. μ ¯ m a x is used as an additional degree of freedom to provide more accurate output estimations through the UKF, such as those of the vehicle speed and sideslip angle, as opposed to attaining an accurate friction coefficient estimation.
The presented estimators are based on the UKF algorithm [24], which operates according to two main steps:
  • The time update (or prediction) step, in which a predicted state vector x k , also known as the a priori state vector, where k is the time step, is computed by using the previous state estimate and the internal nonlinear vehicle model. x k is calculated by estimating the mean of so-called sigma points, which approximate the mean and covariance of the system, representing the state estimate and its associated uncertainty. In this step, the process noise covariance matrix, Q , is used to influence the propagation of the generated sigma points. High values of the elements of Q indicate significant uncertainty and low confidence in the internal model dynamics, which consequently affects the a priori prediction [52].
  • The measurement update (or correction) step, in which x k is corrected to produce the a posteriori state vector estimate, x k + , by multiplying the error between the predicted measurements and the real measured data, also known as the innovation, e y , k , by the Kalman gain, K k , the calculation of which is beyond the scope of this paper:
    x k + = x k + K k e y , k
The measurement noise covariance matrix, R , which represents the uncertainty in the vehicle sensors, is used to influence the spread of a second set of sigma points, generated similarly to those in step 1., and is also directly involved in the calculation of K k . Within the UKF algorithms, larger values of the elements of R imply reduced significance associated with the measurements.
A third covariance matrix, namely the state covariance matrix P , which quantifies the uncertainty in the state vector x k , is used in both steps of the UKF algorithm. While P is a full matrix that is updated iteratively at each time step according to the UKF equations, Q and R are diagonal matrices to reduce the computational load, under the assumption that the process and measurement noise sources are independent from each other. Hence, each element of Q and R correlates to the uncertainty of a state or measurement. In the specific implementation, Q and R are given by:
Q = d i a g [ Q v x Q v y Q r Q ω 1 Q ω 2 Q ω 3 Q ω 4 Q μ ¯ m a x ]
R = d i a g [ R a x R a y R r R ω 1 R ω 2 R ω 3 R ω 4 ]
where the subscripts in the notations of the elements of Q indicate the related state, while the subscripts of the components of R refer to the relevant measurement. Q and R are typically kept constant. In fact, well-calibrated fixed-covariance UKFs are sufficient in the majority of real-world cases, but there are exceptional conditions where this is not true, and thus an ACM setup is beneficial to robustness and convergence. In this study, two UKF variants are presented: (i) the baseline UKF, simply referred to as UKF, with fixed values of the matrices; and (ii) UKF ACM, including the adaptation of the elements of Q and R indicated in the boxes in (8) and (9).
The process and measurement noise covariances of the baseline UKF were obtained via an optimisation routine, minimising a cost function based on the error between the UKF estimated outputs and the measurements, see [49], which was computed from test results obtained along the high-friction handling circuit of the Lommel proving ground (Belgium).

3.2. Wheel Speed Measurement Noise Covariance Adaptation

In UKFs with constant covariances, the accuracy of the v x estimation typically worsens as the magnitude of the longitudinal tyre slip ( σ x ) increases. In fact, when the tyres tend to lock or spin, the wheel speed measurements, despite remaining accurate and robust, are not effective for straightforward vehicle speed estimation. In these conditions, it is more effective to prioritise (i) the speed estimate from the nonlinear internal vehicle model of the time update step, which embeds the tyre slip dynamics of each corner, if the tyre–road friction factor is correctly estimated; and/or (ii) the a x measurement, which, through integration in the time domain, could ideally directly bring the speed profile.
To achieve the desired effect, UKF ACM includes an adaptation of the R matrix values for the four wheel speed measurements, i.e., R ω i , at each time step, according to:
R ω i = R ω , m i n + R ω , m a x R ω , m i n 1 ϵ ω i 2
where R ω , m i n and R ω , m a x are the minimum and maximum values for R ω i , and ϵ ω i is an adaptive weight.
The weight ϵ ω i varies according to the following bivariate normal distribution:
ϵ ω i = e 1 2   a x , ω i a x ¯ σ a 2 + v x , ω i v ^ x , ω ¯ σ v 2
where a x is the measured longitudinal acceleration; v ^ x , ω is a dead reckoning estimated longitudinal speed based on the average of the four wheel speeds; σ a and σ v are tuneable standard deviation parameters for the longitudinal acceleration and speed; and a x , ω i and v x , ω i are the longitudinal acceleration and velocity of the respective corner, obtained from the measured i -th wheel speed and transformed to the vehicle centre of gravity, under the simplifying assumption of zero longitudinal tyre slip:
v x , ω 1 = r t f 2 + ω 1 R 1 cos δ 1 v x , ω 2 = r t f 2 + ω 2 R 2 cos δ 2 v x , ω 3 = r t r 2 + ω 3 R 3 v x , ω 4 = r t r 2 + ω 4 R 4 a x , ω i = v ˙ x , ω i
where R i is the rolling radius of the respective tyre. The bar notation ‘ ¯ ’ in the numerator of (11) denotes a moving average computed on four samples with 10 ms discretisation. The term a x , ω i a x ¯ is an indicator of the absolute rate of variation of the tyre slip in the considered corner, since the reference is represented by the longitudinal vehicle acceleration measurement, while the term v x , ω i v ^ x , ω ¯ highlights the wheels with significantly different longitudinal slip from the average value among the corners. Hence, from the adaptation mechanism perspective, the absence of criticality corresponds to conditions of rather uniform longitudinal tyre slip distribution among the corners, and the absence of significant longitudinal slip rates. During the filter implementation phase, it was extensively verified through simulations and experiments that the combination of the two terms is effective and robust in detecting the critical corners from the viewpoint of the slip ratio dynamics.
Figure 5 shows the adopted ϵ ω i distribution as a function of a x , ω i a x ¯ and v x , ω i v ^ x , ω ¯ . ϵ ω i varies between 0, indicating the presence of critical slip dynamics, and 1, indicating the absence of detected longitudinal slip criticalities. Based on (10), ϵ ω i = 1 corresponds to R ω , m i n , i.e., the condition of maximum weight given to the wheel speed measurement, while ϵ ω i = 0 implies R ω , m a x , i.e., the minimum significance given to the wheel speed signal. The use of a standard deviation approach in (11) enables the capture of the range of potential values and the likelihood of their occurrence within a normal distribution framework. This is particularly useful in the presented adaptation algorithm, in which the two variables are interrelated and exhibit a certain level of correlation in their behaviour [53].
During the filter calibration phase, the parametrisation of the standard deviation coefficients σ a and σ v was carried out through brute-force testing in manoeuvres implying different levels of excitation of the longitudinal vehicle dynamics. R ω , m i n was set to be equal to the R ω value that was optimised for the baseline UKF, providing desirable performance for a very wide range of operating vehicle conditions. Instead, R ω , m a x was determined via a dedicated optimization process. By using the value from van Aalst’s research [37] as the initial condition, a multi-objective genetic algorithm was executed with constant R ω i values, equal among the corners, i.e., R ω 1 = = R ω 4 = R ω . The routine minimises an objective function, J b , which includes six terms, accounting for the relevant normalized root mean square error ( N R M S E ) values between the measurements and the filter estimation outputs:
a r g   m i n R ω   J b = 1 N s b p = 1 N s e w , b , p e n o r m , b 2     with   b = a x , v x , ω 1 , ω 2 , ω 3 , ω 4
where the index b —used both as a main variable and as a subscript—refers to the longitudinal acceleration a x , longitudinal speed v x , and wheel speeds ω i ; N s is the number of samples within the considered manoeuvre; e n o r m , b is the normalisation factor for the corresponding error variable; and e w , b , p is the weighted estimation error of the sample p , which is given by:
e w , b , p = W 1 , b ( b ^ p b p ) i f b ^ p b p < e t h , 1 , b W 2 , b ( b ^ p b p ) i f   e t h , 1 , b b ^ p b p < e t h , 2 , b W 3 , b ( b ^ p b p ) i f b ^ p b p e t h , 2 , b
where e t h , 1 , b and e t h , 2 , b , with e t h , 1 , b >   e t h , 2 , b , are predefined b ^ p b p thresholds; and the scaling factors W 1 , b , W 2 , b , and W 3 , b , with W 3 , b W 2 , b W 1 , b , are set to penalise excessive error values. The optimisation round only covered manoeuvres with extremely critical steady-state and transient longitudinal tyre slip conditions, and the resulting optimal R ω value was set as R ω , m a x . Further trial-and-error calibration of R ω , m a x was implemented at the end of the overall filter adaptation mechanism tuning, by considering the adaptations of the other UKF covariances.

3.3. Process Noise Covariance Adaptation

The process noise covariances Q r , Q v x , and Q v y are characterised by a similar adaptation structure to the wheel speed measurement noise covariances. During the development phase, it was verified that these adaptations tend to improve estimation robustness, especially in extremely critical conditions. For example, this applies to specific counter-steering and drifting manoeuvres, in which the yaw rate and front steering angles often have opposite signs. Hence, the internal model tends to be less accurate, and the UKF must increase its reliance on measured values.
The covariances vary between minimum and maximum set values as functions of ϵ j :
Q j = Q j , m i n + Q j , m a x Q j , m i n 1 ϵ j 2
where the subscript j , with j = r , v x , v y , indicates the relevant state; the minimum and maximum covariance values, Q j , m i n and Q j , m a x , are computed with a similar procedure to the one used for R m i n and R m a x ; and ϵ j is calculated based on a Gaussian distribution law:
ϵ j = e 1 2 l ^ l ¯ σ l 2
where the notation l , with l = r , a x , a y , indicates the relevant measured variable; l ^ is the value of the same variable estimated by the UKF; and σ l is a tuneable standard deviation parameter. Whilst r is used for the computation of ϵ r , the accelerations a x and a y are adopted for ϵ v x and ϵ v y , instead of v x and v y , as it is preferable to build the adaptation upon the estimation error with respect to measured values such as a x and a y .
Extensive testing revealed a significant influence of the measurement covariance for μ ¯ m a x , i.e., Q μ ¯ m a x , on the adaptability and rate of change of μ ¯ m a x , especially in scenarios characterised by high levels of lateral and/or longitudinal tyre slip. Consequently, for this covariance, a different adaptation strategy was developed, using a two-dimensional look-up table, based on two variables related to the lateral and longitudinal dynamics: (a) the moving average of the error between the lateral acceleration value estimated by the UKF and the one measured by the IMU, i.e., e a y = a ^ y a y ¯ ; and (b) the time derivative of the moving average of the longitudinal acceleration error, e ˙ a x = d a ^ x a x ¯ / d t , which is a good indicator of friction-related criticalities of the longitudinal vehicle dynamics.
The Q μ ¯ m a x look-up table was designed to enhance the responsiveness of the state μ ¯ m a x to the actual tyre–road friction level, and to induce rapid reactions of the random walk model. In the generation process, initial one-dimensional maps were generated through brute-force testing, and then were further refined through a gradient-descent optimisation algorithm for each of the purely longitudinal and lateral test scenarios in Section 4.1. The objective was to minimise a cost function based on the relevant key performance indicators (KPIs) in Section 4.2, for the prevalent dynamics excited in the specific manoeuvre. The attained optimal Q μ ¯ m a x values were subsequently combined via interpolation to construct the two-dimensional look-up tables. These were further optimised along manoeuvres with combined high slip conditions. The tuning process was repeated for alternative look-up table input variables. e a y and e ˙ a x emerged as the optimal choices. Figure 6 shows the resulting map implemented in UKF ACM.

4. Test Scenarios and Key Performance Indicators

4.1. Manoeuvres for Performance Assessment

The resulting UKFs were tested: (i) experimentally, in high-friction conditions, at the Lommel proving ground; and (ii) through the high-fidelity CarMaker model, in manoeuvres with very low tyre–road friction levels, which could not be performed within the available experimental facility. Both (i) and (ii) include very challenging conditions from the vehicle dynamics viewpoint. The test scenarios were the following:
  • Test scenario 1: experimental 60 m radius skid pad test, in which the vehicle was slowly accelerated while the driver applied steering angle corrections for tracking the reference trajectory, until the car reached the cornering limit, and could no longer stay within the reference lane.
  • Test scenario 2: experimental handling circuit lap with the vehicle pushed to its limit by a professional test driver. This test was designed to stress the vehicle at high lateral and longitudinal accelerations, and targeted the UKF performance assessment in peak acceleration conditions.
  • Test scenario 3: simulated acceleration and braking test on a very low-friction surface ( μ 0.3). The manoeuvre involved the vehicle accelerating from a standstill to a speed of ~ 70 kmh−1, followed by heavy braking whereby the anti-lock braking system (ABS) module was activated throughout. A conventional ABS algorithm was chosen, which uses a control law based on the combination of longitudinal tyre slip and wheel deceleration [54]. The ABS was fed with the true values of the relevant variables from the high-fidelity vehicle simulation model, and the two filters received the same inputs from the model, i.e., the simulation results were not affected by the presence of the filters. On the contrary, the traction controller was purposely kept inactive during the acceleration phase, which thus implied significant wheel spinning. The overall test targeted the assessment of the vehicle speed estimation performance in extremely critical conditions.
  • Test scenario 4: simulated acceleration test with μ -jump, i.e., with a sudden transition from μ 0.3 to μ 1, which was followed by an immediate acceleration at the vehicle’s maximum capability once the rear wheels had crossed over to the higher friction surface. With the electric motors installed on this vehicle, this equated to a rise in a x from 0 to 6 ms−2 in under 0.3 s.
  • Test scenario 5: simulated slow sinusoidal steering (with a 0.25 Hz frequency and 45 deg amplitude) manoeuvre at μ 0.3 from an initial speed of 70 kmh−1, with the torque demand set to the constant value that would maintain the entry speed if the vehicle was operated in straight line. The steering wheel angle amplitude, δ s w a , m a x , was set to 45 deg, corresponding to a steady-state a y of ~6 ms−2. The test focused on the sideslip angle estimation performance in very low-friction conditions.
Although test scenario 3 involved very low-friction conditions, the filters were initialised with μ ¯ m a x , 0 = 1, corresponding to dry tarmac operation; similarly, to make the manoeuvre more challenging following the tyre-road friction transition, in test scenario 4 the initialisation was set to μ ¯ m a x , 0 = 0.3. In all other scenarios, the initialisation was set to μ ¯ m a x , 0 = 1.
Additionally, to ascertain the robustness of the estimators with respect to realistic variations of typical parameters, a sensitivity analysis was conducted on UKF and UKF ACM with ± 10% independent changes of the vehicle mass m t o t , yaw mass moment of inertia I z , and tyre cornering stiffness C , for which the changes were made to the equivalent scaling factor of the magic formula model, λ k y . Due to the inability to change the experimental vehicle demonstrator, the analysis was performed by varying the parameters within the internal model of the filters. For consistency, the same method was applied to the simulation-based analyses.

4.2. Key Performance Indicators

The performance of UKF and UKF ACM was analysed through the root mean square (hence, the subscript ‘rms’ in the following notations) values of the error between the UKF estimation outputs and the measured data. Hence, the resulting KPIs are (i) e β C O G , r m s , the rms value of the estimation error of the sideslip angle at the centre of gravity; (ii) e v x , r m s , the rms value of the longitudinal speed estimation error; (iii) e σ x , F , r m s , the rms value of the average estimation error on the slip ratios of the front tyres; and (iv) e σ x , R , r m s , the same as (iii), but applied to the slip ratios across the rear axle.
Moreover, for each indicator (i–iv), the relative percentage reduction, U K F   A C M | U K F , of the estimation error brought by UKF ACM with respect to UKF is computed for nominal internal model parameters. The robustness analyses also include consideration of the following indicators:
U K F , r p | U K F = e x , r m s , U K F e x , r m s , U K F , r p e x , r m s , U K F 100 U K F   A C M , r p | U K F = e x , r m s , U K F e x , r m s , U K F   A C M , r p e x , r m s , U K F 100
where e x , r m s , U K F indicates the KPI value for UKF, related to the estimation of the variable x (with x =   β C O G , v x , σ x , F , σ x , R ), for nominal parametrisation of the internal model of the filter; e x , r m s , U K F , r p is the KPI value for UKF, referred to the estimation of the variable x , for a variation of the internal model parameter indicated by the subscript ‘ r p ’ (i.e., robustness parameter, defined such that r p =   m t o t , I z , C ) in the robustness analysis; and e x , r m s , U K F A C M , r p is the equivalent of e x , r m s , U K F , r p for UKF ACM. Hence, U K F , r p | U K F and U K F A C M , r p | U K F are the percentage variations of the estimation performance of UKF and UKF ACM with respect to the benchmarking UKF with nominal parametrisation.

5. Results

5.1. Assessment for Nominal Vehicle Parameters

For test scenarios 1–5, the UKF and UKF ACM outputs are compared against the measured data or simulation results in the following Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11, while Table 2 reports the KPIs for all tests. For simplicity of notation, the symbol ‘ ^ ’ is omitted from the estimated variables.
Test scenario 1. The maximum measured lateral acceleration achieved in this test exceeded 10 ms−2, which implies a real tyre–road friction coefficient slightly above 1. Figure 7 reports the time histories of the main estimated variables and the respective measured profiles. The main useful filter outputs are β C O G and v x . Nevertheless, to assess the overall matching of the estimated vehicle dynamics with the measured behaviour, the figure—similarly to the following ones—also reports additional variables, in this case a x , a y , and r , which are commonly measured in production vehicles, and μ ¯ m a x , the auxiliary state that facilitates estimation convergence. Both UKFs estimate the states well in these quasi-steady-state conditions, with some minor differences. For example, UKF ACM outputs a μ ¯ m a x profile that is higher—by almost 0.1—than the one from UKF. The increased μ ¯ m a x provides modest improvements, of up to 4%, in the estimation of v x and β C O G . This is noticeable at ~34 s, where both UKF ACM estimations are clearly closer to the measured data, e.g., see the inset in the subplot of v x . In general, although not being cause for concern in this specific test, the discrepancies between the measurements and the UKF estimates were caused by the model assumptions (e.g., absence of the pitch and heave motions, and the vertical dynamics of the unsprung masses) and mismatches, which tend to be particularly evident in correspondence with the model nonlinearities associated with tyre, suspension, and aerodynamic behaviour.
Test scenario 2. The handling circuit excited the vehicle cornering response in quasi-steady-state and transient conditions, including significant interactions between the longitudinal and lateral dynamics. In fact, a x often exceeded 6 ms−2, and a y reached the cornering limit with values of over 10 ms−2. In this test, the estimation of β C O G was paramount for active safety purposes, as it would be used by vehicle stability controllers and other chassis controllers. The same time history plots as in scenario 1 are shown in Figure 8. Aside from a jump in μ ¯ m a x at t 80 s, which is attributed to a longitudinal acceleration spike associated with an instance of hard braking, the outputs from the UKF variants are substantially indistinguishable. Both UKFs accurately estimate the sideslip response, with e β , r m s values of ~0.5 deg. This level of accuracy was particularly impressive, considering the signal noise magnitude typical of the measurements from the production sensors of the specific vehicle demonstrator, feeding their inputs to the UKFs. Similarly to the first test scenario, UKF ACM brought a minor improvement in all KPIs, with a ~2% reduction of e β , r m s .
Test scenario 3. This was the first low-friction surface test, in which the μ ¯ m a x initialisation at 1 represented a worst-case scenario, given that the actual friction factor was ~0.3. The vehicle accelerated in the first 7.5 s, achieving an a x peak of ~3 ms−2, which implied significant wheel spinning given the absence of traction control, before braking to a standstill in approximately the same time, with ABS interventions throughout to prevent wheel locking. Figure 9 reports the time histories of the most relevant variables. On the a x profile, at the beginning of the test, there is a very sharp rise in the a x estimated by the filters, as the driver starts to accelerate and the algorithms have been initialised for high-friction conditions, i.e., they expect higher a x values than those actually measured, which are limited by the surface condition. The error between the accelerations and velocities is high, so ϵ ω i drops to its bottom boundary, and subsequently R ω i (see R ω 1 , n in the figure, i.e., the normalised R ω 1 profile for comparison purposes) increases, to indicate that UKF ACM must give low weight to the measured wheel speeds. This initial transient corresponds to large e ˙ a x , which, in turn, at t 1 s, causes a small increase in Q μ ¯ m a x , n , i.e., the normalised value of Q μ ¯ m a x , responsible for the steep drop in μ ¯ m a x . After the drop, as the estimated μ ¯ m a x matches the real friction level, the wheel speed measurements become relevant again. Unlike UKF, UKF ACM is then able to track the measured a x and v x for the remainder of the acceleration phase with great accuracy. For example, the maximum magnitude of the slip ratio estimation error amounts to ~0.6 for UKF ACM, against ~3.2 for UKF. Both UKF variants estimate v x well during the braking period, although UKF ACM is smoother, whereas UKF staggers during the ABS cycling phase. Correspondingly, the amplitude of the a x oscillations during the braking period is roughly halved by UKF ACM. Overall, across the test, UKF ACM reduces e v x , r m s by 94%, from 32.3 ms−1 to 1.9 ms−1, and both e σ x , F , r m s and e σ x , R , r m s by 87%.
The important conclusion is that UKF ACM provided decisively superior v x estimation performance in conditions of very high longitudinal slip. In this and the following tests, the magnitude of the covariance variations for Q μ ¯ m a x and R ω i was significantly higher than for Q r , Q v x and Q v y , so the latter are not shown in the figures, although they are contributing to the performance improvement. In the experiments for test scenarios 1 and 2, it was also verified that the adoption—within UKF—of a Q μ ¯ m a x value facilitating swift variations of μ ¯ m a x , i.e., corresponding to the top level of Q μ ¯ m a x visible in Figure 9 for test scenario 3, would provoke undesirable oscillations of the estimated variables, and unstable filter behaviour during high-friction operation, which confirmed the necessity of the proposed adaptation mechanism.
Test scenario 4. The μ -jump test helped ascertain how well the UKF could detect a swift change in surface, reflecting real-world scenarios, such as road contamination, black ice, or simply transitions from wet to dry tarmac, where accurate estimation, although not being safety-critical, is important not to unnecessarily limit vehicle performance. The results are presented in Figure 10. The road surface changes from low to high μ , and then the vehicle begins to accelerate, just before t = 2 s, with maximum torque demand, which results in an almost step-like increase in a x to ~6 ms−2. UKF fails to match the amplitude and resulting shape of the a x curve, as its model expects the vehicle to be in low-friction conditions, whereas UKF ACM is much more prompt. In fact, the rise in e ˙ a x leads to a spike in Q μ ¯ m a x which allows for μ ¯ m a x to jump from ~0.3 to ~0.75, whereas in UKF, μ ¯ m a x only reaches a value of ~0.38. Hence, UKF ACM better matches the v x measurement curve, with e v x , r m s reducing from ~5 ms−1 to ~1.4 ms−1, and the peak error value dropping from ~10 ms−1 to ~2.5 ms−1. With regards to e σ x , F , r m s and e σ x , R , r m s , UKF ACM reduces these KPIs by ~57%, with the maximum magnitude of the slip ratio estimation error dropping by ~45%, i.e., from 0.175 to 0.095.
Test scenario 5. In the low-friction sinusoidal steering test (Figure 11), μ ¯ m a x was subject to an initial reduction to reasonable levels for both UKF and UKF ACM, occurring at the very beginning of the simulation, when the vehicle accelerated in a straight line to reach the required initial condition, before the relevant section of the manoeuvre reported in the figure. Within the core part of the test, the vehicle exhibits high levels of lateral slip for relatively low a y . During each steering transient, UKF overestimates the lateral acceleration peak, and then immediately afterwards underestimates a y . The e a y oscillations during the transients provoke the Q μ ¯ m a x , n oscillations in the respective plot. These cause the UKF ACM μ ¯ m a x profile to dip as the vehicle slips laterally, which significantly improves the estimation of the peaks and overall shape of a y and β C O G . For example, the maximum magnitude of the β C O G estimation error is reduced from over 1 deg to less than 0.3 deg, a 70% improvement compared to UKF. In conclusion, this test demonstrated the enhanced performance of UKF ACM in its β C O G estimation on low- μ surfaces. Minor improvements, amounting to less than 5%, were also seen in e v x , r m s , e σ x , F , r m s , and e σ x , R , r m s .

5.2. Robustness Analysis

Figure 12 summarises the results of the robustness analysis along test scenarios 2 and 3 (see also Table A1 in the Appendix A for the full dataset). In the plots, the histograms refer to U K F   A C M | U K F , while the error bars represent the range of variation of U K F , r p | U K F and U K F   A C M , r p | U K F , across the sensitivity checks.
The histograms in subplot (a) show that, on the handling circuit in high-friction conditions and with nominal internal model parametrisation, the improvements brought by UKF ACM are marginal (~2% on average) but consistent across the KPIs. The top and bottom caps of the error bars of UKF ACM are all at a higher level than the respective ones for UKF, which corresponds to better performance for the considered parameter range. This is confirmed in Table A1, in which all UKF ACM KPIs are better than those for UKF with the same internal model mismatch. For example, UKF occasionally performs worse by nearly 20% in terms of e β , r m s with respect to its nominal case, particularly when altering the tyre cornering stiffness. The maximum e β , r m s performance decay of the UKF ACM with internal model mismatch with respect to the nominal UKF is only 7%. With the exception of that of e v x , r m s , UKF ACM has shorter error bars than UKF, which represents a smaller spread of values, and as such is symptomatic of robustness. Notably, the error bar for e v x , r m s is longer for UKF ACM, but the increased spread is only in the positive range, where the filter with covariance adaptation markedly improves the estimation.
For test scenario 3 (see subplot (b)), the spread of e σ x , F , r m s and e σ x , R , r m s is roughly the same for both UKF variants, while the one for e v x , r m s is lower for UKF ACM. Nevertheless, the difference in estimation accuracy is so wide that the error bars for the two filter variants do not have any overlap.

6. Conclusions

This paper presented a novel vehicle state estimator, i.e., the so-called UKF ACM, of the vehicle longitudinal and lateral components of the vehicle velocity, tyre slip ratios, as well as sideslip and slip angles. The novelty lies in the adaptive formulations of the process and measurement noise covariances. For fairness of assessment, the fixed covariances of a well-calibrated baseline UKF were also used as the baseline values for UKF ACM. The two UKF variants were compared in a variety of experimental and simulated manoeuvres, including very low tyre–road friction conditions as well as longitudinal and lateral dynamics.
Based on the results in Section 5, the main conclusions are the following:
  • In high-friction conditions near the limit of handling, the performance of the baseline UKF with fixed covariances was very similar to that of UKF ACM, with the latter providing small but still noteworthy benefits.
  • In extreme longitudinal slip cases on low μ with highly incorrect friction level initialisation within the filters, UKF ACM performed significantly better than UKF. In fact, unlike UKF, the variations in the process and measurement noise covariances Q μ ¯ m a x and R ω i , related to the friction random walk model and wheel speed measurements, enabled UKF ACM to promptly detect the actual tyre–road friction level, and achieve highly accurate speed and slip ratio estimation. Similarly, UKF ACM was very effective in identifying instantaneous and extreme changes in μ , with the related positive impact in terms of the resulting estimation.
  • In extreme lateral slip conditions on very low μ surfaces, the increased sensitivity of μ ¯ m a x of UKF ACM allowed it to outperform the baseline UKF, with improvements of over 50% in e β , r m s . This highlighted a clear safety improvement, as accurate sideslip angle estimation is necessary for typical vehicle chassis controllers.
  • UKF ACM has shown notable robustness with respect to UKF. In fact, when varying—within the internal models of the filters—parameters that would realistically change during real-world vehicle operation, (i) the UKF ACM KPI decay was maintained within the tolerable range of 10% of the original values for nominal conditions, while UKF experienced a maximum performance reduction that approached 20%; and (ii) the UKF ACM results were the same or better, e.g., by more than 80% in the high longitudinal tyre slip conditions of test scenario 3, than the corresponding UKF ones.
Future work will include (i) experimental testing on surfaces with low and variable tyre–road friction conditions; (ii) the production-oriented implementation of UKF ACM on automotive control hardware; (iii) the integration of machine learning algorithms in the UKF adaptation scheme; and (iv) the comparison of the proposed formulation with the most promising machine learning techniques for vehicle dynamics state estimation.

Author Contributions

Conceptualization, A.A., S.D.P., S.v.A., D.I., P.G. and A.S.; Methodology, S.D.P., S.v.A., K.P., D.I., P.G. and A.S.; Software, A.A., S.D.P., P.S., S.v.A. and D.I.; Validation, A.A., S.D.P., P.S., S.v.A., K.P. and D.I.; Formal analysis, A.A., S.D.P., P.S. and S.v.A.; Investigation, A.A., S.D.P., S.v.A., K.P., D.I. and P.G.; Resources, E.B.; Data curation, A.A., P.S. and K.P.; Writing—original draft, A.A. and S.D.P.; Writing—review & editing, P.S., E.B., P.G. and A.S.; Visualization, P.S.; Supervision, E.B., P.G. and A.S.; Project administration, E.B. and A.S.; Funding acquisition, E.B. and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Horizon 2020 Programme of the European Commission under grant agreement 824250 (EVC1000 project).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

Author Stefano De Pinto was employed by the company McLaren Automotive. Authors Sebastiaan van Aalst, Kylian Praet and Emilie Boulay were employed by the company Tenneco Automotive. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Appendix A

Table A1. KPIs for the robustness analysis for test scenarios 2 and 3.
Table A1. KPIs for the robustness analysis for test scenarios 2 and 3.
e β , r m s e v x , r m s e σ x , F , r m s e σ x , R , r m s
[deg][ms−1][−][−]
Test scenario 2—Experimental handling circuit
UKF0.5102.5150.01570.0173
UKF m t o t +10%0.5712.6440.01570.0175
UKF m t o t –10%0.5112.4910.01580.0171
UKF I z +10%0.5212.5150.01570.0173
UKF I z −10%0.5042.5140.01570.0173
UKF C +10%0.5372.4890.01550.0171
UKF C –10%0.6042.5350.01620.0177
UKF ACM0.5012.5040.01560.0171
UKF ACM m t o t +10%0.5412.6390.01570.0174
UKF ACM m t o t –10%0.5252.4210.01570.0170
UKF ACM I z +10%0.5242.5040.01570.0173
UKF ACM I z −10%0.5162.5040.01570.0173
UKF ACM C +10%0.5252.5030.01540.0170
UKF ACM C –10%0.5542.5290.01620.0177
Test scenario 3—Simulated acceleration and braking with μ 0.3
UKF-32.2541.39881.0325
UKF m t o t +10%-31.6601.35751.0047
UKF m t o t –10%-34.8761.44541.0576
UKF I z +10%-32.2531.39881.0325
UKF I z −10%-32.2761.40061.0337
UKF C +10%-32.2111.40151.0376
UKF C –10%-32.4401.41611.0441
UKF ACM-1.9250.18280.1361
UKF ACM m t o t +10%-3.1530.14890.1339
UKF ACM m t o t –10%-3.1950.19600.1527
UKF ACM I z +10%-2.4490.13570.1085
UKF ACM I z −10%-1.8950.17450.1296
UKF ACM C +10%-1.8860.17440.1305
UKF ACM C –10%-1.9940.18740.1395

References

  1. Mazzilli, V.; De Pinto, S.; Pascali, L.; Contrino, M.; Bottiglione, F.; Mantriota, G.; Gruber, P.; Sorniotti, A. Integrated chassis control: Classification, analysis and future trends. Annu. Rev. Control 2021, 51, 172–205. [Google Scholar] [CrossRef]
  2. Theunissen, J.; Tota, A.; Gruber, P.; Dhaens, M.; Sorniotti, A. Preview-based techniques for vehicle suspension control: A state-of-the-art review. Annu. Rev. Control 2021, 51, 206–235. [Google Scholar] [CrossRef]
  3. Ersal, T.; Kolmanovsky, I.; Masoud, N.; Ozay, N.; Scruggs, J.; Vasudevan, R. Connected and automated road vehicles: State of the art and future challenges. Veh. Syst. Dyn. 2020, 58, 672–704. [Google Scholar] [CrossRef]
  4. Olovsson, T.; Svensson, T.; Wu, J. Future connected vehicles: Communications demands, privacy and cyber-security. Commun. Transp. Res. 2022, 2, 100056. [Google Scholar] [CrossRef]
  5. Xu, Q.; Li, K.; Wang, J.; Yuan, Q.; Yang, Y.; Chu, W. The status, challenges, and trends: An interpretation of technology roadmap of intelligent and connected vehicles in China. J. Intell. Connect. Veh. 2022, 5, 1–7. [Google Scholar] [CrossRef]
  6. van Zanten, A.T.; Erhardt, R.; Pfaff, G. VDC, The Vehicle Dynamics Control System of Bosch. SAE J. Passeng. Cars Part 1 1995, 104, 1419–1436. [Google Scholar]
  7. Singh, K.B.; Arat, M.A.; Taheri, S. Literature review and fundamental approaches for vehicle and tire state estimation. Veh. Syst. Dyn. 2019, 57, 1643–1665. [Google Scholar] [CrossRef]
  8. Jin, X.; Yin, G.; Chen, N. Advanced Estimation Techniques for Vehicle System Dynamic State: A Survey. Sensors 2019, 19, 4289. [Google Scholar] [CrossRef]
  9. Savaresi, S.; Tanelli, M. Active Braking Control Systems Design for Vehicles, 1st ed.; Springer: London, UK, 2010. [Google Scholar]
  10. Lin, H.; Liu, Y.; Li, S.; Qu, X. How generative adversarial networks promote the development of intelligent transportation systems: A survey. IEEE/CAA J. Autom. Sin. 2023, 10, 1781–1796. [Google Scholar] [CrossRef]
  11. Wei, W.; Shaoyi, B.; Lanchun, Z.; Kai, Z.; Yongzhi, W.; Weixing, H. Vehicle sideslip angle estimation based on general regression neural network. Math. Probl. Eng. 2016, 2016, 3107910. [Google Scholar] [CrossRef]
  12. Melzi, S.; Sabbioni, E. On the vehicle sideslip angle estimation through neural networks: Numerical and experimental results. Mech. Syst. Signal Process. 2011, 25, 2005–2019. [Google Scholar] [CrossRef]
  13. Ghosh, J.; Tonoli, A.; Amati, N. Deep learning based virtual sensor for vehicle sideslip angle estimation: Experimental Results. SAE Tech. Pap. 2018, 1089. [Google Scholar] [CrossRef]
  14. Bonfitto, A.; Feraco, S.; Tonoli, A.; Amati, N. Combined regression and classification artificial neural networks for sideslip angle estimation and road condition identification. Veh. Syst. Dyn. 2020, 58, 1766–1787. [Google Scholar] [CrossRef]
  15. González, L.P.; Sánchez, S.S.; Garcia-Guzman, J.; Boada, M.J.L.; Boada, B.L. Simultaneous Estimation of Vehicle Roll and Sideslip Angles through a Deep Learning Approach. Sensors 2020, 20, 3679. [Google Scholar] [CrossRef]
  16. Bertipaglia, A.; Mol, D.; Alirezaei, M.; Happee, R.; Shyrokau, B. Model-based vs Data-driven Estimation of Vehicle Sideslip Angle and Benefits of Tyre Force Measurements. arXiv 2022, arXiv:2206.15119. [Google Scholar]
  17. Chindamo, D.; Lenzo, B.; Gadola, M. On the Vehicle Sideslip Angle Estimation: A Literature Review of Methods, Models, and Innovations. Appl. Sci. 2018, 8, 355. [Google Scholar] [CrossRef]
  18. Liu, J.; Wang, Z.; Zhang, L.; Walker, P. Sideslip Angle Estimation of Ground Vehicles: A Comparative Study. IET Control Theory Appl. 2020, 14, 3490–3505. [Google Scholar] [CrossRef]
  19. Chen, B.C.; Hsieh, F.C. Sideslip angle estimation using extended Kalman filter. Veh. Syst. Dyn. 2008, 46, 353–364. [Google Scholar] [CrossRef]
  20. Bersani, M.; Vignati, M.; Mentasti, S.; Arrigoni, S.; Cheli, F. Vehicle state estimation based on Kalman filters. In Proceedings of the AEIT International Conference of Electrical and Electronic Technologies for Automotive, Turin, Italy, 2–4 July 2019; pp. 1–6. [Google Scholar]
  21. Doumiati, M.; Victorino, A.C.; Charara, A.; Lechner, D. Onboard Real-Time Estimation of Vehicle Lateral Tire–Road Forces and Sideslip Angle. IEEE/ASME Trans. Mechatron. 2011, 16, 601–614. [Google Scholar] [CrossRef]
  22. Wang, P.; Fan, X.; Chen, X.; Yi, J.; He, S. UKF Estimation Method of Centroid Slip Angle for Vehicle Stability Control. Int. J. Control Autom. Syst. 2023, 21, 2259–2266. [Google Scholar] [CrossRef]
  23. Laenza, A.; Mantriota, G.; Reina, G. On the vehicle dynamics prediction via model-based observation. Veh. Syst. Dyn. 2023, 1–22, in press. [Google Scholar] [CrossRef]
  24. Wan, E.A.; Van Der Merwe, R. The unscented Kalman filter for nonlinear estimation. In Proceedings of the IEEE 2000 Adaptive Systems for Signal Processing, Communications, and Control Symposium, Lake Louise, AB, Canada, 1–4 October 2000; pp. 153–158. [Google Scholar]
  25. Xiong, L.; Xia, X.; Lu, Y.; Liu, W.; Gao, L.; Song, S.; Han, Y.; Yu, Z. IMU-Based Automated Vehicle Slip Angle and Attitude Estimation Aided by Vehicle Dynamics. Sensors 2019, 19, 1930. [Google Scholar] [CrossRef] [PubMed]
  26. Madhusudhanan, A.K.; Corno, M.; Holweg, E. Vehicle sideslip estimator using load sensing bearings. Control Eng. Pract. 2016, 54, 46–57. [Google Scholar] [CrossRef]
  27. Lenzo, B.; Ottomano, G.; Strano, S.; Terzo, M.; Tordela, C. A Physical-Based Observer for Vehicle State Estimation and Road Condition Monitoring. In Proceedings of the IOP Conference Series: Materials Science and Engineering, Rome, Italy, 22–24 July 2020. [Google Scholar]
  28. Villano, E.; Lenzo, B.; Sakhnevych, A. Cross-combined UKF for vehicle sideslip angle estimation with a modified Dugoff tire model: Design and experimental results. Meccanica 2021, 56, 2653–2668. [Google Scholar] [CrossRef]
  29. Li, L.; Jia, G.; Ran, X.; Song, J.; Wu, K. A variable structure extended Kalman filter for vehicle sideslip angle estimation on a low friction road. Veh. Syst. Dyn. 2014, 52, 280–308. [Google Scholar] [CrossRef]
  30. Antonov, S.; Fehn, A.; Kugi, A. Unscented Kalman filter for vehicle state estimation. Veh. Syst. Dyn. 2011, 49, 1497–1520. [Google Scholar] [CrossRef]
  31. Heidfeld, H.; Schünemann, M.; Kasper, R. UKF-based State and tire slip estimation for a 4WD electric vehicle. Veh. Syst. Dyn. 2019, 58, 1479–1496. [Google Scholar] [CrossRef]
  32. Wielitzka, M.; Busch, A.; Dagen, M.; Ortmaier, T.; Serra, G. Unscented Kalman filter for state and parameter estimation in vehicle dynamics. In Kalman Filters-Theory for Advanced Applications; IntechOpen: Rijeka, Croatia, 2018; pp. 56–75. [Google Scholar]
  33. Wenzel, T.A.; Burnham, K.J.; Blundell, M.V.; Williams, R.A. Dual extended Kalman filter for vehicle state and parameter estimation. Veh. Syst. Dyn. 2006, 44, 153–171. [Google Scholar] [CrossRef]
  34. Liao, Y.W.; Borrelli, F. An adaptive approach to real-time estimation of vehicle sideslip, road bank angles, and sensor bias. IEEE Trans Veh Technol. 2019, 68, 7443–7454. [Google Scholar] [CrossRef]
  35. Bechtoff, J.; Koenig, L.; Isermann, R. Cornering stiffness and sideslip angle estimation for integrated vehicle dynamics control. IFAC-PapersOnLine 2016, 49, 297–304. [Google Scholar] [CrossRef]
  36. Bertipaglia, A.; Shyrokau, B.; Alirezaei, M.; Happee, R. A two-stage bayesian optimisation for automatic tuning of an unscented kalman filter for vehicle sideslip angle estimation. IEEE Intell. Veh. Symp. 2022, 2022, 670–677. [Google Scholar]
  37. van Aalst, S. Virtual Sensing for Vehicle Dynamics. Ph.D. Thesis, KU Leuven, Leuven, Belgium, 2020. [Google Scholar]
  38. Chen, T.; Cai, Y.; Chen, L.; Xu, X.; Jiang, H.; Sun, X. Design of Vehicle Running States-Fused Estimation Strategy Using Kalman Filters and Tire Force Compensation Method. IEEE Access 2019, 7, 87273–87287. [Google Scholar] [CrossRef]
  39. Zhujie, S.; Fen, L.; Youqun, Z.; Shougang, M. Vehicle state and parameter estimation based on improved adaptive dual extended Kalman filter with variable sliding window. In Proceedings of the 2022 6th CAA International Conference on Vehicular Control and Intelligence (CVCI), Nanjing, China, 28–30 October 2022; pp. 1–6. [Google Scholar]
  40. Chen, Y.; Yan, H.; Li, Y. Vehicle State Estimation Based on Sage–Husa Adaptive Unscented Kalman Filtering. World Electr. Veh. J. 2023, 14, 167. [Google Scholar] [CrossRef]
  41. Wang, Y.; Li, Y.; Zhao, Z. State Parameter Estimation of Intelligent Vehicles Based on an Adaptive Unscented Kalman Filter. Electronics 2023, 12, 1500. [Google Scholar] [CrossRef]
  42. Pang, H.; Wang, P.; Wang, M.; Hu, C. On accurate estimation of vehicle lateral states based on an improved adaptive unscented Kalman filter. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2022, 2022, 09544070221132328. [Google Scholar] [CrossRef]
  43. Zhang, Y.; Li, M.; Zhang, Y.; Hu, Z.; Sun, Q.; Lu, B. An Enhanced Adaptive Unscented Kalman Filter for Vehicle State Estimation. IEEE Trans. Instrum. Meas. 2022, 71, 1–12. [Google Scholar] [CrossRef]
  44. van Aalst, S.; Naets, F.; Boulkroune, B.; De Nijs, W.; Desmet, W. An adaptive vehicle sideslip estimator for reliable estimation in low and high excitation driving. IFAC-PapersOnLine 2018, 51, 243–248. [Google Scholar] [CrossRef]
  45. Xia, X.; Hashemi, E.; Xiong, L.; Khajepour, A. Autonomous Vehicle Kinematics and Dynamics Synthesis for Sideslip Angle Estimation Based on Consensus Kalman Filter. IEEE Trans. Control Syst. Technol. 2023, 31, 179–192. [Google Scholar] [CrossRef]
  46. Li, X.; Xu, N.; Li, Q.; Guo, K.; Zhou, J. A fusion methodology for sideslip angle estimation on the basis of kinematics-based and model-based approaches. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2020, 234, 1930–1943. [Google Scholar] [CrossRef]
  47. Pacejka, H.B. Tire and Vehicle Dynamics, 2nd ed.; SAE International: Warrendale, PA, USA, 2006. [Google Scholar]
  48. ISO 7401; Road Vehicles—Lateral Transient Response Test Methods—Open-Loop Test Methods. ISO: Geneva, Switzerland, 2011.
  49. Mazzilli, V.; Ivone, D.; De Pinto, S.; Pascali, L.; Contrino, M.; Tarquinio, G.; Gruber, P.; Sorniotti, A. On the benefit of smart tyre technology on vehicle state estimation. Veh. Syst. Dyn. 2022, 60, 3694–3719. [Google Scholar] [CrossRef]
  50. Jung, S.; Kim, T.Y.; Yoo, W.S. Advanced slip ratio for ensuring numerical stability of low-speed driving simulation. Part I: Longitudinal slip ratio. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2019, 233, 2000–2006. [Google Scholar] [CrossRef]
  51. Kim, T.Y.; Jung, S.; Yoo, W.S. Advanced slip ratio for ensuring numerical stability of low-speed driving simulation: Part II—Lateral slip ratio. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2019, 233, 2903–2911. [Google Scholar] [CrossRef]
  52. Welch, G.; Bishop, G. An Introduction to the Kalman Filter; University of North Carolina at Chapel Hill: Chapel Hill, NC, USA, 2006; pp. 2–12. [Google Scholar]
  53. Lai, C.D.; Balakrishnan, N. Continuous Bivariate Distributions; Springer New York: New York, NY, USA, 2009. [Google Scholar]
  54. Burckhardt, M.; Reimpell, J. Fahrwerktechnik: Radschlupe-Regelsysteme; Vogel: Würzburg, Germany, 1993. [Google Scholar]
Figure 1. The case study Audi e-tron at the Ebbenjarka testing facility (Sweden), including the Kistler Correvit S-350 sensor, providing accurate measurements of the velocity components.
Figure 1. The case study Audi e-tron at the Ebbenjarka testing facility (Sweden), including the Kistler Correvit S-350 sensor, providing accurate measurements of the velocity components.
Sensors 24 00436 g001
Figure 2. Top view schematic of the internal model of the estimator with indication of the main variables and parameters.
Figure 2. Top view schematic of the internal model of the estimator with indication of the main variables and parameters.
Sensors 24 00436 g002
Figure 3. Examples of experimental validation of the CarMaker model and the internal UKF model on (a) a 40 m radius skid pad on dry tarmac; (b) a sinusoidal steering manoeuvre on dry tarmac, from ~80 kmh−1; and (c) a sinusoidal steering test at ~40 kmh−1 on an icy surface (t is time).
Figure 3. Examples of experimental validation of the CarMaker model and the internal UKF model on (a) a 40 m radius skid pad on dry tarmac; (b) a sinusoidal steering manoeuvre on dry tarmac, from ~80 kmh−1; and (c) a sinusoidal steering test at ~40 kmh−1 on an icy surface (t is time).
Sensors 24 00436 g003
Figure 4. Overview schematic of the presented UKF ACM state estimator.
Figure 4. Overview schematic of the presented UKF ACM state estimator.
Sensors 24 00436 g004
Figure 5. Plot of ϵ ω i as a function of a x , ω i a x ¯ and v x , ω i v ^ x , ω ¯ .
Figure 5. Plot of ϵ ω i as a function of a x , ω i a x ¯ and v x , ω i v ^ x , ω ¯ .
Sensors 24 00436 g005
Figure 6. Plot of Q μ ¯ m a x as a function of e a y and e ˙ a x .
Figure 6. Plot of Q μ ¯ m a x as a function of e a y and e ˙ a x .
Sensors 24 00436 g006
Figure 7. UKF comparison plots for an experimental skid pad manoeuvre (test scenario 1).
Figure 7. UKF comparison plots for an experimental skid pad manoeuvre (test scenario 1).
Sensors 24 00436 g007
Figure 8. UKF comparison plots for an experimental high-friction handling circuit lap (test scenario 2).
Figure 8. UKF comparison plots for an experimental high-friction handling circuit lap (test scenario 2).
Sensors 24 00436 g008
Figure 9. UKF comparison plots for an acceleration and braking test with μ 0.3 (test scenario 3). For the immediacy of notation, the front left and rear left vehicle corners are indicated with ‘FL’ and ‘RL’.
Figure 9. UKF comparison plots for an acceleration and braking test with μ 0.3 (test scenario 3). For the immediacy of notation, the front left and rear left vehicle corners are indicated with ‘FL’ and ‘RL’.
Sensors 24 00436 g009
Figure 10. UKF comparison plots for an acceleration test with μ -jump from low to high tyre–road friction conditions (test scenario 4).
Figure 10. UKF comparison plots for an acceleration test with μ -jump from low to high tyre–road friction conditions (test scenario 4).
Sensors 24 00436 g010
Figure 11. UKF comparison plots for the sinusoidal steering test with μ 0.3 (test scenario 5).
Figure 11. UKF comparison plots for the sinusoidal steering test with μ 0.3 (test scenario 5).
Sensors 24 00436 g011
Figure 12. Histograms of U K F   A C M | U K F , and error bars for U K F , r p | U K F and U K F   A C M , r p | U K F , for the considered robustness analyses, along (a) the experimental handling circuit (test scenario 2); and (b) the simulated acceleration and braking test in low-friction conditions (test scenario 3).
Figure 12. Histograms of U K F   A C M | U K F , and error bars for U K F , r p | U K F and U K F   A C M , r p | U K F , for the considered robustness analyses, along (a) the experimental handling circuit (test scenario 2); and (b) the simulated acceleration and braking test in low-friction conditions (test scenario 3).
Sensors 24 00436 g012
Table 1. Main parameters of the case study vehicle.
Table 1. Main parameters of the case study vehicle.
DescriptionSymbolValueUnit
Vehicle mass in real testing conditions m t o t 2843[kg]
Yaw mass moment of inertia I z 4124[kgm2]
Front semi-wheelbase a 1.47[m]
Rear semi-wheelbase b 1.46[m]
Front track width t f 1.60[m]
Rear track width t r 1.60[m]
Centre of gravity height h C O G 0.63[m]
Wheel radius R 0.38[m]
Table 2. KPIs for UKF and UKF ACM along the considered test scenarios.
Table 2. KPIs for UKF and UKF ACM along the considered test scenarios.
e β , r m s U K F   A C M | U K F e v x , r m s U K F   A C M | U K F e σ x , F , r m s U K F   A C M | U K F e σ x , R , r m s U K F   A C M | U K F
[deg][%][ms−1][%][−][%][−][%]
Test Scenario 1—Experimental skid pad
UKF0.429-1.531-0.0264-0.0287-
UKF ACM0.415−3.15%1.472−3.83%0.0257−2.54%0.0282−1.84%
Test Scenario 2—Experimental handling circuit
UKF0.510-2.515-0.0157-0.0173-
UKF ACM0.501−1.76%2.504−0.43%0.0156−0.83%0.0171−0.87%
Test Scenario 3—Simulated acceleration and braking with μ   0.3
UKF--32.254-1.3988-1.0325-
UKF ACM--1.925−94%0.1828−87%0.1361−87%
Test scenario 4—Simulated acceleration with μ -jump
UKF--4.993-0.0956-0.0928-
UKF ACM--1.373−72%0.0422−56%0.0401−57%
Test scenario 5—Simulated sinusoidal steering test with μ   0.3
UKF0.508-0.985-0.0327-0.0328-
UKF ACM0.233−54%0.948−3.71%0.0322−1.41%0.0323−1.40%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alshawi, A.; De Pinto, S.; Stano, P.; van Aalst, S.; Praet, K.; Boulay, E.; Ivone, D.; Gruber, P.; Sorniotti, A. An Adaptive Unscented Kalman Filter for the Estimation of the Vehicle Velocity Components, Slip Angles, and Slip Ratios in Extreme Driving Manoeuvres. Sensors 2024, 24, 436. https://doi.org/10.3390/s24020436

AMA Style

Alshawi A, De Pinto S, Stano P, van Aalst S, Praet K, Boulay E, Ivone D, Gruber P, Sorniotti A. An Adaptive Unscented Kalman Filter for the Estimation of the Vehicle Velocity Components, Slip Angles, and Slip Ratios in Extreme Driving Manoeuvres. Sensors. 2024; 24(2):436. https://doi.org/10.3390/s24020436

Chicago/Turabian Style

Alshawi, Aymen, Stefano De Pinto, Pietro Stano, Sebastiaan van Aalst, Kylian Praet, Emilie Boulay, Davide Ivone, Patrick Gruber, and Aldo Sorniotti. 2024. "An Adaptive Unscented Kalman Filter for the Estimation of the Vehicle Velocity Components, Slip Angles, and Slip Ratios in Extreme Driving Manoeuvres" Sensors 24, no. 2: 436. https://doi.org/10.3390/s24020436

APA Style

Alshawi, A., De Pinto, S., Stano, P., van Aalst, S., Praet, K., Boulay, E., Ivone, D., Gruber, P., & Sorniotti, A. (2024). An Adaptive Unscented Kalman Filter for the Estimation of the Vehicle Velocity Components, Slip Angles, and Slip Ratios in Extreme Driving Manoeuvres. Sensors, 24(2), 436. https://doi.org/10.3390/s24020436

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop