Skip to main content

ORIGINAL RESEARCH article

Front. Control Eng., 27 October 2021
Sec. Control and Automation Systems
This article is part of the Research Topic Linear Parameter Varying Systems Modelling, Identification and Control View all 6 articles

Linear Parameter Varying Path Tracking Control for Over-Actuated Electric Vehicles

Alex GimondiAlex GimondiMatteo Corno
Matteo Corno*Sergio Matteo SavaresiSergio Matteo Savaresi
  • Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Milano, Italy

This paper discusses a multi-layer linear parameter varying path tracking controller for autonomous full electric vehicles with 4 wheel drive. Many approaches for path tracking are present in the literature, but the majority of them utilize only the steering wheel. Nowadays, the massive introduction of electric vehicles opens new opportunities for vehicle dynamics control; multiple electric motors, via differential torque, offer an alternative way to influence vehicle motion. The multi-layer structure, managing separately steering wheel and electric motors, permits the exploitation of pre-existing control systems and maintains a straightforward tuning process. Furthermore, the robustness of the path tracking controller is enhanced by introducing tire nonlinearities. Linear parameter varying models offer a direct way to incorporate tire characteristics in the design process. Exploiting H technique the scheduled controller can be seamlessly synthezised. The proposed strategy has been compared in simulation to a benchmark integrated, i.e., a unique controller that manages all the actuators, approach. During double lane changes beyond the vehicle limits, the multi-layer controller guarantees stability and excellent tracking with an average error of 12 and 24 cm on high and low grip respectively.

1 Introduction

Two major trends are driving the automotive industry nowadays. On one hand, self-driving cars promise to reduce road congestion (Stogios et al., 2019), increase safety (Fagnant and Kockelman, 2015) and reduce pollution (Huang et al., 2018). The literature on the topic is very rich. Autonomous driving systems are usually developed around three interacting but separate tasks: sensing and perception, planning, and control (Amer et al., 2017). This paper focuses on the last one: our goal is to track a known, obstacle-free, reference trajectory. Furthermore, trajectory tracking comprises two objectives: to follow a desired path and a reference speed. We circumscribe our work to lateral control: a separate longitudinal controller will track the target velocity. On the other hand, electric vehicles are being massively introduced in the market due to their ecological potentialities (Zinnari et al., 2020; Frieske et al., 2013). Besides environmental aspects, electric vehicles have many advantages; in the context of this work, we are mainly interested in the versatility they offer in terms of drivetrain design. Among the possible drivetrain configurations, the one with 4 electric motors (4WD), one per wheel, offers interesting opportunities for vehicle dynamics control. For example, the possibility to generate a differential torque between left and right wheels permits to influence the vehicle lateral dynamics: this is called torque vectoring (TV). Mainly, TV systems regulate the vehicle yaw rate to improve safety and drivability (Lucchini et al., 2020).

In our work, we consider an autonomous electric vehicle with 4WD. Within this framework, two possibilities exist to control the lateral dynamics: integrated approaches, i.e. techniques that coordinately actuate the steering angle and the motors; or multi-layer approaches, i.e., two controllers that manage respectively the steering angle and the motors.

Integrated approaches are widely discussed in the literature (Goodarzi et al., 2008; Yakub and Mori, 2015; Falcone et al., 2008; Peng et al., 2020). Regarding multi-layer strategies, the single control problems have been investigated in depth. Concerning path tracking, the literature presents many interesting approaches ranging from: Model Predictive Control (MPC) (Guo et al., 2020; Katriniok et al., 2013; Beal and Gerdes, 2013; Falcone et al., 2007; Liniger et al., 2017), sliding mode (Tagne et al., 2013), nonlinear control (Peters et al., 2011), to potential fields (Rossetter et al., 2004). MPC is considered the state of the art. Nevertheless, its use is sensitive to model accuracy and requires a high computational demand. Therefore, MPC implementation entails a non-negligible amount of effort. For TV literature analysis the reader can refer to De Novellis et al., 2014); Kaiser et al., 2014).

However, the interaction of the two control systems is hardly discussed in the literature. Only Chatzikomis et al. (2018) investigate the topic: integrated approaches are compared to a blind (each subsystem is tuned independently) multi-layer structure. The authors highlight the fact that integrated approaches are not robust to friction variation. From a pure control performance standpoint, an integrated approach should be better. However, automotive systems are often built on pre-existing structures and control layers, often designed by different teams, if not different suppliers. Redesigning the entire vehicle dynamics control from scratch would be expensive and discard years of experience and certification efforts. For these reasons, it is interesting and could be beneficial to study the effect of multi-layer approaches.

In this context, the paper designs a nested control structure. The inner control loop uses the electric motors to regulate the yaw rate dynamics. The outer control loop takes care of controlling the steering angle to track a desired trajectory. The core of the scheme is a quasi-linear parameter varying (q-LPV) H controller. It incorporates in an easy way the tire nonlinearities and guarantees a straightforward tuning process. The multi-layer q-LPV/H, compared to fixed structure and integrated approaches, showed excellent performance. The proposed solutions are validated in simulation with IPG CarMaker on ISO double lane change (DLC) at different velocities and road conditions. This work extends our previous papers (Roselli et al., 2017; Corno et al., 2021) from several standpoints.

1. Differently from the cited literature, we consider an over-actuated vehicle. The control system thus also solves a control allocation problem.

2. The role of the preview distance is better formalized than in Corno. In fact, in Section 3, we introduce the course rate preview control that, to the best of our knowledge, has not been discussed before.

3. We propose a formal way of considering a hierarchical structure by designing a multi-layer controller that explicitly considers the torque vectoring system,

4. We improve the robustness of the path tracking controller by using the q-LPV control structure to consider the tire nonlinearities.

Linear parameter varying (LPV) systems (Mohammadpour and Scherer, 2012) have been extensively used in automotive applications (Li et al., 2021). They have been employed to consider the velocity dependency of the vehicle dynamics (Corno et al., 2021), to schedule a path following controller based on vehicle states (Alcalá et al., 2017), to manage the priority between different actuators (Fergani et al., 2016), or, as in our case, to consider tire nonlinearities (Németh et al., 2016). In particular, since the scheduling parameters are functions of the model state, the system is classified as q-LPV.

To sum up, the contributions of this work are:

• the design of the multi-layer course rate preview structure tuned exploiting the LPV formulation. Especially, the novelty consists in considering the TV system in the designing of the higher level controller;

• a comparison between integrated and multi-layer approaches, highlighting the effectiveness of the proposed approach.

The paper is organized as follows. In Section 2, the models used for control design are presented. Section 3 describes the control structures: a multi-layer H is a necessary step to attain the q-LPV controller. Furthermore, some insight on the parameter choice are given. Section 4 validates the proposed strategies and compares their peculiarities. Finally, Section 5 summarizes the concluding remarks.

2 Vehicle Dynamics Modeling

We describe the vehicle lateral dynamics with the classical single-track model (Abe, 2015):

may=Fyfcosδw+FyrIzṙ=Fyflf+Fyrlr+Mz(1)

where ay is the lateral acceleration, Fy f,r are respectively the front and rear lateral tire forces, r is the yaw rate, m the vehicle mass, Iz is the vehicle moment of inertia, lf,r are the distances from the centre of gravity (CoG) respectively to the front and rear axle, Mz is an external yaw moment, and δw is the steering angle at the wheel.

The tire forces are influenced by many parameters and variables, e.g. vertical load Fz, camber angle γ, lateral slip α, longitudinal slip λ, and friction μ. In general:

Fy=fFz,α,λ,γ,μ.(2)

There exist different ways to describe the nonlinear function f(⋅), e.g., Pacejka (2005); Kiencke and Nielsen (2000); Fiala (1954). Among them, we have preferred the Pacejka’s magic formula for its accuracy. Figure 1 shows an example of f(⋅) for different road surfaces. In this phase, the longitudinal slip is set to 0, and, for the sake of simplicity, also the camber angle is neglected. For a fixed Fz, (Eq. 2) can be linearized around a given lateral slip angle ᾱ:

Fyαα=ᾱ=fᾱ+Cᾱαᾱ(3)

where C represents the generalized1 cornering stiffness defined as:

Cᾱ=fαα=ᾱ(4)

for the sake of clarity, subscripts f,r for α, Fy, C have been dropped out in the previous relations. To obtain the linear single-track model, we can write the lateral acceleration as:

ay=vβ̇+r(5)

where v is the vehicle velocity and β is the sideslip angle, i.e., the angle between the single-track longitudinal axis and the velocity vector. Furthermore, with the hypothesis of small angles, we can express the wheel side slips as (Kiencke and Nielsen, 2000):

αf=β+lfrvδwαr=βlrrv.(6)

FIGURE 1
www.frontiersin.org

FIGURE 1. Tire lateral forces for different friction coefficients at Fz = 6000 N.

Substituting (Eq. 3) (linearized around ᾱ=0), (Eq. 5), and (Eq. 6) in (Eq. 1), the linear single-track model arises:

β̇ṙ=Cf+Crmv1+CrlrCflfmv2CrlrCflfIzCflf2+Crlr2Izvβr+Cfmv0CflfIz1IzδwMz(7)

where Cf = Cf(0) and Cr = Cr(0).

In order to control the vehicle position, the previous model is augmented with its relative motion with respect to a known trajectory (Kapania and Gerdes, 2015). Figure 2 summarizes the main variables. In particular, the dynamics of the heading difference Δψ, i.e., the angle between the vehicle and the trajectory orientation, is

Δψ̇=rρv(8)

where ρ is the trajectory curvature. The lateral error dynamics become:

ėcg=vβ+Δψ.(9)

FIGURE 2
www.frontiersin.org

FIGURE 2. Augmented single-track variables.

Eventually, we obtain the augmented single-track state space model:

β̇ṙΔψ̇ėcg=Cf+Crmv1+CrlrCflfmv200CrlrCflfIzCflf2+Crlr2Izv000100v0v0βrΔψecg+Cfmv00CflfIz1Iz000v000δwMzρ(10)

where the vector δwMzρ represents the model inputs. δw and Mz are control variables; whereas ρ is a measurable disturbance (usually in autonomous driving the path is known in advance).

The vehicle model is completed by the actuators dynamics. We consider a fully electric vehicle with 4WD and a steer-by-wire system. Each wheel is connected to an electric motor able to generate both positive and negative torques. Figure 3 plots the static characteristics of the electric motors. A first-order low-pass filter with a bandwidth of 20 Hz and a time delay of 10 ms complete the motor dynamics model.

FIGURE 3
www.frontiersin.org

FIGURE 3. Front and rear motors characteristics.

Besides, the steering actuator we consider is composed of an electric motor linked to the steering column via a belt transmission. The system is modeled with a second-order transfer function

Gsteers=ωsteers2+2ξsteerωsteer+ωsteer2eTacts(11)

where ωsteer = 2π4.1, ξsteer = 0.1 and Tact = 0.08 have been experimentally identified (Roselli et al., 2017; Corno et al., 2021).

3 Course Rate Preview Control

In order to control the lateral error, we propose a nested loop approach. The inner loop tracks a desired course rate; it actuates both the steering angle and the motors torque and exploits a preview term usually available in autonomous driving applications. The outer loop corrects the desired course rate based on the actual path tracking error.

We have chosen a multi-layer structure to manage the over-actuated configuration. The course rate controller sets the steering angle; in addition, a standard Torque Vectoring (TV) system generates the motor torques. The scheme in Figure 4 shows the general control structure.

FIGURE 4
www.frontiersin.org

FIGURE 4. Multi-layer course rate preview control scheme. The resulting steering wheel command is fed both to the vehicle and to a TV system that generates a suitable yaw moment.

Firstly, for the sake of clarity, we briefly describe the TV block. Afterwards, the course rate preview is designed. Initially, we have adopted the H technique to tune the inner loop. Then, exploiting the same structure, we have applied the H on a LPV version of the single-track in which the tire nonlinearities have been introduced.

3.1 Torque Vectoring

A TV system is generally composed of three blocks (De Novellis et al., 2014): a reference generator, a high-level controller; and a torque allocator.

The reference generator has the objective of providing a suitable yaw rate reference (rref). It is defined as in Rajamani (2011):

rref=vL+v2Kusdesδw,if vL+v2Kusdesδw<0.85μgv0.85μgvsignδw,if vL+v2Kusdesδw0.85μgv(12)

where L is the wheelbase, g is the gravitation acceleration, Kus des is the desired understeering coefficient. By changing Kus des we can impose the desired steering behavior; in our case study, we require a neutral steering vehicle, thus Kus des is set to 0.

The high-level controller tracks the desired yaw rate reference by generating a yaw moment. We design this controller using model-based techniques on the nominal transfer function from the yaw moment to the yaw rate obtained from (Eq. 7) at 90 km/h. A PI controller achieves a bandwidth of 1.5 Hz with a phase margin of 80°.

The torque allocator solves the control allocation problem. It distributes the driver torque request to obtain the desired yaw moment. To formulate the optimization problem, we can write a torque balance at vehicle CoG:

Mztot=cosδwFxfrFxfldf2+sinδwFxfr+Fxfllf+cosδwFyfr+Fyfllf+sinδwFyflFyfrdf2+FxrrFxrldr2Fyrl+Fyrrlr(13)

where df,r are respectively the front and rear track widths, and Fx f,r are the longitudinal tire forces. With the hypothesis of small δw (Eq. 13) becomes:

Mztot=FxfrFxfldf2+Fyfr+Fyfllf+FxrrFxrldr2Fyrl+Fyrrlr.(14)

To influence the lateral dynamics we can exploit the longitudinal forces that are directly controllable with the electric motors. Thus, the yaw moment that can be generated is:

Mz=FxfrFxfldf2+FxrrFxrldr2(15)

Supposing that a linear relation exists between torque and force: Tw=FxRw, where Tw is the torque generated by each electric motor at wheel side, and Rw is the wheel radius; (Eq. 15) can be written as

u=ATv(16)

where A=df2Rwdf2Rwdr2Rwdr2Rw and Tv=TflTfrTrlTrr. Note that possible inaccuracies introduced by the hypothesis done are compensated by the closed loop control system.

To compute the optimum torque allocation a minimization problem is set:

minimizeTiJTisubject toATv=MzTihbiTilbi(17)

where hbi and lbi are respectively the highest and lowest motor torque for each corner. The cost function J(Ti)=(Td4Tfl)2+(Td4Tfr)2+(Td4Trl)2+(Td4Trr)2 aims at matching the total driver request Td. We formalize the problem as a constrained optimization. This provides a tool to prioritize the lateral stability of the vehicle, while still considering the driver’s torque request. If both the driver’s torque request and the yaw moment request cannot be met, the yaw moment takes priority. Furthermore, note that in this implementation we uniformly distributed the torque among the four wheels. The formulation is of course open to other torque distributions. The control allocation is solved online as a weighted least-square problem as in Harkegard (2002).

3.2 H Course Rate Preview

The proposed path tracking controller is based on the course rate, defined as ϕ=r+β̇; and it is composed of three elements: the course rate loop, the reference pre-filter and the lateral error loop.

The course rate loop is an H controller tuned on the generalized plant of Figure 5, where:

w represents the exogenous input (ϕref);

u is the control variable (δw cmd);

z represents the weighted exogenous outputs ([w1 w2 w3]);

y represents the measured output (eϕ).

FIGURE 5
www.frontiersin.org

FIGURE 5. H control scheme.

P(s) is the generalized plant model including the vehicle dynamics, the actuators, and the shaping functions.

In the multi-layer approach, the course rate loop controls a vehicle that is equipped with a TV system. There are two ways to consider this in the design. If the underlying control system is protected, one could design the controller on an equivalent single-track model that describes the expected closed-loop dynamics. Alternatively, if one has the access to the controller, it (along with its states) can be embedded in the generalized plant. We followed this second approach, augmenting (Eq. 7) with the TV controller:

β̇ṙζ̇=Cf+Crmv1+CrlrCflfmv20CrlrCflfIzCflf2+Crlr2+kpvIzvkiIz010βrζ+Cfmv0CflfIzkpvIzL+Kusdesv20vIzL+Kusdesv2δwδwcmd(18)

where ζ represents the TV integrator state, and kp,i the TV controller parameters. Note that the TV is fed with the steering wheel angle command. The objective of the H controller is to find a Kϕ(s) such that the H norm of the transfer function from w to z is minimized. The three exogenous outputs are defined as w1 = W1(s)eϕ, w2 = W2(s)δw cmd and w3 = W3(s)ϕ, where W1,2,3(s) are frequency weights designed to shape the closed-loop response. In particular, W1(s) and W3(s) have been chosen to guarantee a bandwidth of 1 Hz and a good reference tracking at low frequency, W2(s) has been chosen to limit the control effort at high frequency (Skogestad and Postlethwaite, 2007; Corno et al., 2021). Figure 6 shows the inverse of the weights and the obtained sensitivity (S=(I+GδwcmdϕKϕ)1), control sensitivity (Q = KϕS), and complementary sensitivity (T=SGδwcmdϕKϕ) functions. Gδwcmdϕ is the transfer function from the steering wheel command to the course rate obtainable from the block scheme of Figure 4.

FIGURE 6
www.frontiersin.org

FIGURE 6. Bode magnitude plots of the sensitivity (S), the control sensitivity (Q) and the complementary sensitivity (T) functions. The dashed lines are the inverse of the weights used to shape the plant response. Multi-layer H case.

The reference pre-filter, Fρ(s), generates the reference for the course rate loop. The reference depends on the curvature of the target path. Ideally Gϕrefϕ(s), i.e., the transfer function from ϕref to ϕ, should be 1. Thus, the simplest choice would be: Fρ(s)=Gϕrefϕ(s)1. In practice, we are interested in having a unitary transfer function only up to a certain bandwidth; we assume that most manoeuvrers excite the course rate up to 1 Hz (Roselli et al., 2017). Thanks to this hypothesis, we approximate Gϕrefϕ(s) as a time delay (Figure 7) which makes Fρ(s) a time advance that provides the future reference course rate to the controller; no casualty issue arises since the trajectory is known in autonomous driving applications.

FIGURE 7
www.frontiersin.org

FIGURE 7. Bode plot of Gϕrefϕ(s) (continuos) an its approximation (dashed).

The lateral error loop: finally, the external loop can be closed to guarantee a tracking also in the case of small errors in the course rate control. Kecg(s) is tuned based on the transfer function from δe to ecg in which the course rate loop is considered. The controller, tuned to guarantee null steady state error, yields a bandwidth of 0.47 Hz with a phase margin of 67°. For the sake of completeness, the gain of Kecg(s) is scheduled according to the vehicle velocity:

Kecg0=Kecg,90090/3.6v(19)

where Kecg,90 is the controller transfer function tuned at 90 km/h.

In the following, the pre-filter and the lateral error controller are kept fixed for all the proposed configurations. Since the minimum bandwidth requirement is the same between the controllers, the higher levels do not need to be changed.

3.3 Q-LPV/H Course Rate Preview

The course rate controller considers a linear single-track model. When a vehicle performs emergency obstacle avoidance manoeuvrers, a linear model may not accurately describe the vehicle dynamics. If reliable estimations of the vehicle side slip angle and friction coefficient (as the one provided by Grip et al. (2008) are available, one can improve the performance and robustness of the controller by considering the tire nonlinearities in the design. The LPV framework provides a straightforward and systematic way to do so. Model (18) can be expressed as a q-LPV system by choosing the parameter vector

p=CfαfmvCfαflfIzCrαrmvCrαrlrIz.(20)

By linearizing the single-track around an arbitrary ᾱfr we obtain:

ẋ=A0+A1p1+A2p2+A3p3+A4p4x+B0+B1p1+B2p2+B3p3+B4p4ũ+E(21)

where:

A0=0100kpIzkiIz010A1=1lfv0000000A2=0001lfv0000A3=1lrv0000000A4=0001lrv0000B0=000kpIzvL+Kusdesv20vL+Kusdesv2B1=100000B2=001000B3=000000B4=000000E=ffᾱf+frᾱr+Cfᾱf+Crᾱrmvlfffᾱflrfrᾱr+lfCfᾱflrCrᾱrIz0x=βrζũ=δwδwcmd.

Matrix E can be regarded as a constant input to the system that will be compensated by guaranteeing a high low frequency gain in the controller.

Some remarks concerning the choice of the parameters are due. Firstly, note that the actual physical scheduling parameters are the front and rear generalized cornering stiffnesses. However, in defining the LPV system, we have employed a larger parameter set. The combined parameter vector p has the advantage of better scaling the numerical range of the parameter variability thus improving the numerical properties of the resulting Linear Matrix Inequalities (LMI). By fixing as higher bounds (C̄f,r) the cornering stiffnesses on dry asphalt and as lower bounds (C̲f,r) their half, with the proposed solution p1,3 range from 1 to 3 and p2,4 from 50 to 100. If we had chosen p = [Cf(αf) Cr(αr)] the parameters would have spaced from to 103 to 2 ⋅ 103. By reducing the parameter variation magnitude, we obtain a more numerically robust problem that results in a lower performance index when numerically solved.

Secondly, we highlight the fact that there exists a trade-off between parameter set amplitude and performance. When it comes to decide the parameters range, one may be tempted to consider the absolute maximum and minimum values of the generalized cornering stiffnesses as bounds. It this way, all the possible system trajectories would have been considered. However, synthesizing a controller that stabilizes the system in all the theoretically reachable conditions would hardly satisfy the limitation imposed on the actuator usage generating an inadequate solution. Conversely, if we shrink too much the parameters space the advantages of scheduling vanishes.

The last remark regards the vehicle speed. We have decided to not schedule the system according to the velocity. Indeed, considering that its dynamics depend on the vehicle mass and the generalized cornering stiffness variation depends on fα, it is sensible to assume that the velocity varies slower than Cf,r and its impact can be accounted for with a slower and simpler controller scheduling. In the next section this issue will be further discussed.

Based on the above LPV model, we design a self-scheduled H controller using the same structure and weights rationale of the previous fixed structure controller. Figure 8 exemplifies the closed-loop transfer functions for a particular p. The pre-filter and the lateral error loop tuning are unchanged.

FIGURE 8
www.frontiersin.org

FIGURE 8. Bode magnitude plots of the sensitivity (S), the control sensitivity (Q) and the complementary sensitivity (T) functions. The q-LPV system has been evaluated for p̄=(34Cfmv,34CflfIz,34Crmv,34CrlrIz). The dashed lines are the inverse of the weights used to shape the plant response. Multi-layer q-LPV/H case.

4 Validation and Discussion

The analysis has been carried out in simulation exploiting IPG CarMaker2, a multi-body vehicle simulation software. It considers critical elements such as tire nonlinear characteristics (Magic Formula 6.1), suspension dynamics and kinematics and load transfer phenomena. The controlled vehicle is a SUV with a 4WD drivetrain, Table 1 lists its main parameters.

TABLE 1
www.frontiersin.org

TABLE 1. Vehicle parameters.

In order to assess the course rate preview effectiveness a DLC manoeuvre is executed at different velocities and with different road friction coefficients. For a more objective analysis three performance indexes are considered:

• root mean square of the lateral error: RMSecg=1T0Tecg(γ)dγ;

• steering actuator usage: SAU=1T0T|δwcmd(γ)|dγ;

• yaw moment usage: YMU=1T0T|Mz(γ)|dγ.

First of all, we analyze the behavior of the course rate preview. It is compared against a lateral controller based on Kapania and Gerdes (2015). The benchmark features a feedback contribution based on the look-ahead error and a feedforward term; TV is kept unchanged and active. The controller tuning minimizes the RMSecg. Figure 9 shows a DLC on high grip at 80 km/h; which is within the vehicle handling limits. The figure plots the reference and executed trajectory, the course rate and the error at the centre of gravity. From the figure, one appreciates how the control scheme perfectly tracks the reference trajectory. The RMSecg of the multi-layer H is about 3 cm and the peak error is 10 cm. The benchmark solution shows a RMSecg of 7 cm and a peak error of 20 cm. The course rate preview working mechanism is well depicted in the course rate plot. The required course rate (dotted black) is fed in advance as reference (dashed black), thanks to the pre-filter the course rate loop lag is attenuated. The curvature preview term is more effective than the use of a look-ahead distance. Indeed, ϕmultilayer is in advance with respect to ϕbenchmark.

FIGURE 9
www.frontiersin.org

FIGURE 9. DLC at 80 km/h on μ = 1. Comparison between the multi-layer H and the proposed benchmark solution.

For a comprehensive analysis, the H and the q-LPV/H are compared to an integrated approach. Thus, an additional controller is introduced as benchmark: a MISO q-LPV/H. Indeed, the multi-layer approach has the advantage of keeping the actuators management separated. This is convenient from the architecture point of view, but may be suboptimal from the performance standpoint. To design the integrated approach, many choices of the controllable outputs can be done. We decided to consider the system as multi-input ([δw cmd, Mz]) single-output ([ϕ]). The lateral error dynamics can be considered slower than the course rate one. Utilizing the same approach of the course rate preview: the inner loop, a MISO q-LPV/H, computes the steering angle and the yaw moment, the pre-filter and lateral error loop are left unchanged. The q-LPV/H approach previously described has been integrated adding a weight on the yaw moment and the MISO q-LPV/H controllers is synthesized.

Table 2 resumes the simulation results. We consider two manoeuvres both beyond the handling limits of the vehicle. The perfect tracking of the DLC at 100 km/h on high grip would require a maximal lateral acceleration of 12.4 m/s2, the DLC on μ = 0.4 yields a peak lateral acceleration of 7.4 m/s2. Given the tires limits, the maximal available lateral acceleration is 1 g and 0.4 g on high and low grip respectively. From the table, the following remarks are due:

• at 100 km/h, the multi-layer approaches maximize the usage of the yaw moment, fully exploiting the over-actuated architecture. The MISO q-LPV/H exploits both the actuators, but its lateral error is beyond 20 cm. The multi-layer q-LPV/H achieves the best performance considering all the indexes.

• at 80 km/h, despite the high steering wheel usage, the multi-layer approach keeps the vehicle stable. Comparing the q-LPV approaches, the MISO q-LPV/H reaches slightly better performance. Anyway, by looking at the whole tests the multi-layer q-LPV/H demonstrates a better average behavior.

TABLE 2
www.frontiersin.org

TABLE 2. Performance indexes of the DLC for different speeds and friction coefficients.

We analyze also the time domain behaviour of the controllers. Figure 10 shows the DLC at 100 km/h on high grip. In particular, we compare the multi-layer q-LPV/H and MISO q-LPV/H approaches. It is evident that the first one has better performance. From the course rate plot, it can be observed that the MISO controller has a lower bandwidth with respect to the multi-layer approach and struggles to follow the reference. Anyway, in the design of the MISO controller we have requested the same minima requirements. However, LMI solution returns a performance index of 10. Tuning the multi-input single-output controller is not straightforward. The complexity of the MISO setup effectively prevents it to attain its full potential.

FIGURE 10
www.frontiersin.org

FIGURE 10. DLC at 100 km/h on μ = 1. Multi-layer q-LPV/H vs MISO q-LPV/H.

To discuss the validity of the hypothesis done to obtain the controller we further analyze the manoeuvre. First of all, the single-track model was derived with the hypothesis of small angles. From the top plot of Figure 11 it can be observed that the steering wheel angle at the wheels reaches non negligible values. Thus, the control is robust with respect to high steering commands. We highlight that the most critical non-linearity in the single-track mode is the tire linearization. Indeed, the bottom plot shows that we cover both the linear and nonlinear regions. Considering this behavior, as we did with LPV formulation, is paramount. Lastly, the comment about prioritizing the cornering stiffness scheduling with respect to the velocity is motivated. From the figure it is clear that the vehicle speed is varying slower than αf,r, and as consequence slower than the generalized cornering stiffness. Thus, the velocity dependency is considered with more classical gain scheduling techniques, e.g. linearization scheduling (Rugh and Shamma, 2000).

FIGURE 11
www.frontiersin.org

FIGURE 11. DLC at 100 km/h on μ = 1. Multi-layer q-LPV/H vs MISO q-LPV/H. Detail of the steering wheel angle at the wheel, the vehicle velocity and the front/rear wheel slip angles.

Figure 12 compares the multi-layer and MISO approaches on low grip. The interesting fact is that both approaches obtain a very similar behavior. Thus, the multi-layer strategy reaches the performance of an integrated method.

FIGURE 12
www.frontiersin.org

FIGURE 12. DLC at 80 km/h on μ = 0.4. Multi-layer q-LPV/H vs MISO q-LPV/H.

Finally, Figure 13 shows the RMSecg of the DLCs for different lower bounds choice. As previously pointed out, there exists a trade-off between parameters range and performance. On one hand, the wider the parameters space the more difficult is the synthesis of the controller resulting in unsatisfactory performance. On the other hand, a narrow parameter space nullifies the effectiveness of the scheduling strategy, i.e., guarantee stability and leverage actuators peculiarities (this is evident from Table 2).

FIGURE 13
www.frontiersin.org

FIGURE 13. Sensitivity of the RMSecg with respect to the lower bounds of the parameters space (C̲f,r). The lower bounds are expressed as a percentage of the nominal cornering stiffness.

5 Conclusion

In this paper, we have proposed a multi-layer q-LPV/H control scheme. The multi-layer approach, separating the actuators management, exploits possible pre-existing control systems and simplifies the design process. We highlight the fact that the multi-layer course rate preview is aware of the underlying TV system since it is considered in the tuning process. Furthermore, it has been compared to MISO technique demonstrating excellent performance.

Besides, the introduction of tire nonlinearities, exploiting LPV formulation, proved to increase the conditions in which a DLC is feasible. The choice of the parameters and their range have been discussed.

The next steps consist in introducing the velocity as an additional scheduling variable and validating experimentally the proposed strategies.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

AG contributed to the design of control law, simulation and writing of the original draft. MC helped with the methodology, guidance, review and editing. SS helped with review and editing. All authors contributed to the article and approved the submitted version.

Conflict of Interest

The 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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Footnotes

1usually, the term cornering stiffness refers to the linearized curve for α = 0

2https://ipg-automotive.com/

References

Abe, M. (2015). Vehicle Handling Dynamics: Theory and Application. Butterworth-Heinemann.

Google Scholar

Alcalá, E., Puig, V., Quevedo, J., and Escobet, T. (2017). Gain Scheduling Lpv Control Scheme for the Autonomous Guidance Problem Using a Dynamic Modelling Approach. London: IET. https://arxiv.org/abs/1712.00390.

Google Scholar

Amer, N. H., Zamzuri, H., Hudha, K., and Kadir, Z. A. (2017). Modelling and Control Strategies in Path Tracking Control for Autonomous Ground Vehicles: A Review of State of the Art and Challenges. J. Intell. Robot Syst. 86, 225–254. doi:10.1007/s10846-016-0442-0

CrossRef Full Text | Google Scholar

Beal, C. E., and Gerdes, J. C. (2013). Model Predictive Control for Vehicle Stabilization at the Limits of Handling. IEEE Trans. Contr. Syst. Technol. 21, 1258–1269. doi:10.1109/TCST.2012.2200826

CrossRef Full Text | Google Scholar

Chatzikomis, C., Sorniotti, A., Gruber, P., Zanchetta, M., Willans, D., and Balcombe, B. (2018). Comparison of Path Tracking and Torque-Vectoring Controllers for Autonomous Electric Vehicles. IEEE Trans. Intell. Veh. 3, 559–570. doi:10.1109/TIV.2018.2874529

CrossRef Full Text | Google Scholar

Corno, M., Panzani, G., Roselli, F., Giorelli, M., Azzolini, D., and Savaresi, S. M. (2021). An Lpv Approach to Autonomous Vehicle Path Tracking in the Presence of Steering Actuation Nonlinearities. IEEE Trans. Contr. Syst. Technol. 29, 1766–1774. doi:10.1109/TCST.2020.3006123

CrossRef Full Text | Google Scholar

De Novellis, L., Sorniotti, A., Gruber, P., and Pennycott, A. (2014). Comparison of Feedback Control Techniques for Torque-Vectoring Control of Fully Electric Vehicles. IEEE Trans. Veh. Technol. 63, 3612–3623. doi:10.1109/tvt.2014.2305475

CrossRef Full Text | Google Scholar

Fagnant, D. J., and Kockelman, K. (2015). Preparing a Nation for Autonomous Vehicles: Opportunities, Barriers and Policy Recommendations. Transportation Res. A: Pol. Pract. 77, 167–181. doi:10.1016/j.tra.2015.04.003

CrossRef Full Text | Google Scholar

Falcone, P., Borrelli, F., Asgari, J., Tseng, H. E., and Hrovat, D. (2007). Predictive Active Steering Control for Autonomous Vehicle Systems. IEEE Trans. Contr. Syst. Technol. 15, 566–580. doi:10.1109/TCST.2007.894653

CrossRef Full Text | Google Scholar

Falcone, P., Eric Tseng, H., Borrelli, F., Asgari, J., and Hrovat, D. (2008). Mpc-based Yaw and Lateral Stabilisation via Active Front Steering and Braking. Vehicle Syst. Dyn. 46, 611–628. doi:10.1080/00423110802018297

CrossRef Full Text | Google Scholar

Fergani, S., Sename, O., and Dugard, L. (2016). An LPV/H∞ Integrated Vehicle Dynamic Controller Integrated Vehicle Dynamic Controller. IEEE Trans. Veh. Technol. 65, 1880–1889. doi:10.1109/TVT.2015.2425299

CrossRef Full Text | Google Scholar

Fiala, E. (1954). Seitenkraften am rollenden Luftreifen. VdI 96, 973–979.

Google Scholar

Frieske, B., Kloetzke, M., and Mauser, F. (2013). “Trends in Vehicle Concept and Key Technology Development for Hybrid and Battery Electric Vehicles,” in 2013 World Electric Vehicle Symposium and Exhibition (EVS27), Barcelona, Spain, November 17–20, 2013, 1–12. doi:10.1109/evs.2013.6914783

CrossRef Full Text | Google Scholar

Goodarzi, A., Sabooteh, A., and Esmailzadeh, E. (2008). Automatic Path Control Based on Integrated Steering and External Yaw-Moment Control. Proc. Inst. Mech. Eng. K: J. multi-body Dyn. 222, 189–200. doi:10.1243/14644193jmbd120

CrossRef Full Text | Google Scholar

Grip, H. F., Imsland, L., Johansen, T. A., Fossen, T. I., Kalkkuhl, J. C., and Suissa, A. (2008). Nonlinear Vehicle Side-Slip Estimation with Friction Adaptation. Automatica 44, 611–622. doi:10.1016/j.automatica.2007.06.017

CrossRef Full Text | Google Scholar

Guo, N., Zhang, X., Zou, Y., Lenzo, B., and Zhang, T. (2020). A Computationally Efficient Path-Following Control Strategy of Autonomous Electric Vehicles with Yaw Motion Stabilization. IEEE Trans. Transp. Electrific. 6, 728–739. doi:10.1109/TTE.2020.2993862

CrossRef Full Text | Google Scholar

Harkegard, O. (20022002). Efficient Active Set Algorithms for Solving Constrained Least Squares Problems in Aircraft Control Allocation. Proc. 41st IEEE Conf. Decis. Control. 2, 1295–1300. doi:10.1109/CDC.2002.1184694

CrossRef Full Text | Google Scholar

Huang, K., Yang, X., Lu, Y., Mi, C. C., and Kondlapudi, P. (2018). Ecological Driving System for Connected/automated Vehicles Using a Two-Stage Control Hierarchy. IEEE Trans. Intell. Transport. Syst. 19, 2373–2384. doi:10.1109/tits.2018.2813978

CrossRef Full Text | Google Scholar

Kaiser, G., Liu, Q., Hoffmann, C., Korte, M., and Werner, H. (2014). Lpv Torque Vectoring for an Electric Vehicle with Experimental Validation. IFAC Proc. Volumes 47, 12010–12015. 19th IFAC World Congress. doi:10.3182/20140824-6-ZA-1003.00163

CrossRef Full Text | Google Scholar

Kapania, N. R., and Gerdes, J. C. (2015). Design of a Feedback-Feedforward Steering Controller for Accurate Path Tracking and Stability at the Limits of Handling. Vehicle Syst. Dyn. 53, 1687–1704. doi:10.1080/00423114.2015.1055279

CrossRef Full Text | Google Scholar

Katriniok, A., Maschuw, J. P., Christen, F., Eckstein, L., and Abel, D. (2013). “Optimal Vehicle Dynamics Control for Combined Longitudinal and Lateral Autonomous Vehicle Guidance,” in 2013 European Control Conference (ECC), Orlando, FL, December 12–15, 2011, 974–979. doi:10.23919/ecc.2013.6669331

CrossRef Full Text | Google Scholar

Kiencke, U., and Nielsen, L. (2000). Automotive Control Systems: For Engine, Driveline, and Vehicle. Heidelberg: Springer.

Google Scholar

Li, P., Nguyen, A.-T., Du, H., Wang, Y., and Zhang, H. (2021). Polytopic Lpv Approaches for Intelligent Automotive Systems: State of the Art and Future Challenges. Mech. Syst. Signal Process. 161, 107931. doi:10.1016/j.ymssp.2021.107931

CrossRef Full Text | Google Scholar

Liniger, A., Domahidi, A., and Morari, M. (2017). Optimization-based Autonomous Racing of 1: 43 Scale Rc Cars. ArXiv abs/1711.07300.Hoboken, NJ: Wiley

Google Scholar

Lucchini, A., Formentin, S., Corno, M., Piga, D., and Savaresi, S. M. (2020). Torque Vectoring for High-Performance Electric Vehicles: An Efficient Mpc Calibration. IEEE Control. Syst. Lett. 4, 725–730. doi:10.1109/LCSYS.2020.2981895

CrossRef Full Text | Google Scholar

Mohammadpour, J., and Scherer, C. W. (2012). Control of Linear Parameter Varying Systems with Applications. Springer Science & Business Media.

Google Scholar

Németh, B., Gáspár, P., and Bokor, J. (2016). “Lpv-based Integrated Vehicle Control Design Considering the Nonlinear Characteristics of the Tire,” in 2016 American Control Conference (ACC), Boston: IEEE, 6893–6898. doi:10.1109/ACC.2016.7526758

CrossRef Full Text | Google Scholar

Pacejka, H. (2005). Tire and Vehicle Dynamics. Elsevier.

Google Scholar

Peng, H., Wang, W., An, Q., Xiang, C., and Li, L. (2020). Path Tracking and Direct Yaw Moment Coordinated Control Based on Robust Mpc with the Finite Time Horizon for Autonomous Independent-Drive Vehicles. IEEE Trans. Veh. Technol. 69, 6053–6066. doi:10.1109/TVT.2020.2981619

CrossRef Full Text | Google Scholar

Peters, S., Frazzoli, E., and Iagnemma, K. (20112011). Differential Flatness of a Front-Steered Vehicle with Tire Force Control. IEEE/RSJ Int. Conf. Intell. Robots Syst., 298–304. doi:10.1109/iros.2011.6094800

CrossRef Full Text | Google Scholar

Rajamani, R. (2011). Vehicle Dynamics and Control. Springer US. Mechanical Engineering Series.

Google Scholar

Roselli, F., Corno, M., Savaresi, S. M., Giorelli, M., Azzolini, D., Irilli, A., and Panzani, G. (2017). “H∞ Control with Look-Ahead for Lane Keeping in Autonomous Vehicles,” in 2017 IEEE Conference on Control Technology and Applications (CCTA), Kohala Coast, Hawaii, August 27–30, 2017, 2220–2225. doi:10.1109/CCTA.2017.8062781

CrossRef Full Text | Google Scholar

Rossetter, E. J., Switkes, J., and Gerdes, J. C. (2004). Experimental Validation of the Potential Field Lanekeeping System. Int. J. Automotive Technol. 5, 95–108. doi:10.1109/ACC.2003.1240417

CrossRef Full Text | Google Scholar

Rugh, W. J., and Shamma, J. S. (2000). Research on Gain Scheduling. Automatica 36, 1401–1425. doi:10.1016/S0005-1098(00)00058-3

CrossRef Full Text | Google Scholar

Skogestad, S., and Postlethwaite, I. (2007). Multivariable Feedback Control: Analysis and Design, Vol. 2. Hoboken, NJ: Wiley.

Google Scholar

Stogios, C., Kasraian, D., Roorda, M. J., and Hatzopoulou, M. (2019). Simulating Impacts of Automated Driving Behavior and Traffic Conditions on Vehicle Emissions. Transportation Res. D: Transport Environ. 76, 176–192. doi:10.1016/j.trd.2019.09.020

CrossRef Full Text | Google Scholar

Tagne, G., Talj, R., and Charara, A. (2013). “Higher-order Sliding Mode Control for Lateral Dynamics of Autonomous Vehicles, with Experimental Validation,” in 2013 IEEE Intelligent Vehicles Symposium (IV), Gold Coast, Australia, June 23–26, 2013, 678–683. doi:10.1109/IVS.2013.6629545

CrossRef Full Text | Google Scholar

Yakub, F., and Mori, Y. (2015). Comparative Study of Autonomous Path-Following Vehicle Control via Model Predictive Control and Linear Quadratic Control. Proc. Inst. Mech. Eng. D: J. automobile Eng. 229, 1695–1714. doi:10.1177/0954407014566031

CrossRef Full Text | Google Scholar

Zinnari, F., Strada, S., Tanelli, M., Formentin, S., Savaia, G., and Savaresi, S. M. (2020). “Mining the Electrification Potential of Fuel-Based Vehicles Mobility Patterns: a Data-Based Approach,” in 2020 IEEE International Conference on Human-Machine Systems (ICHMS), September 7–9, 2020, 1–S. doi:10.1109/ICHMS49158.2020.9209557

CrossRef Full Text | Google Scholar

Keywords: path tracking, linear parameter varying (LPV) systems, torque vectoring (TV), autonomous vehicles, electric vehicles, control

Citation: Gimondi A, Corno M and Savaresi SM (2021) Linear Parameter Varying Path Tracking Control for Over-Actuated Electric Vehicles. Front. Control. Eng. 2:750190. doi: 10.3389/fcteg.2021.750190

Received: 30 July 2021; Accepted: 11 October 2021;
Published: 27 October 2021.

Edited by:

Fernando A.C.C. Fontes, University of Porto, Portugal

Reviewed by:

Fotis Nicholas Koumboulis, National and Kapodistrian University of Athens, Greece
Basilio Lenzo, University of Padua, Italy

Copyright © 2021 Gimondi, Corno and Savaresi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Matteo Corno, bWF0dGVvLmNvcm5vQHBvbGltaS5pdA==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.