RESEARCH ARTICLE

Extended model predictive control scheme for smooth path following of autonomous vehicles

  • Qianjie LIU 1 ,
  • Shuang SONG 1 ,
  • Huosheng HU 2 ,
  • Tengchao HUANG 1 ,
  • Chenyang LI 1 ,
  • Qingyuan ZHU , 1
Expand
  • 1. Department of Mechanical and Electrical Engineering, Xiamen University, Xiamen 361102, China
  • 2. School of Computer Science and Electronic Engineering, University of Essex, Colchester CO4 3SQ, UK

Received date: 12 Jul 2021

Accepted date: 22 Sep 2021

Published date: 15 Mar 2022

Copyright

2022 Higher Education Press 2022.

Abstract

This paper presents an extended model predictive control (MPC) scheme for implementing optimal path following of autonomous vehicles, which has multiple constraints and an integrated model of vehicle and road dynamics. Road curvature and inclination factors are used in the construction of the vehicle dynamic model to describe its lateral and roll dynamics accurately. Sideslip, rollover, and vehicle envelopes are used as multiple constraints in the MPC controller formulation. Then, an extended MPC method solved by differential evolution optimization algorithm is proposed to realize optimal smooth path following based on driving path features. Finally, simulation and real experiments are carried out to evaluate the feasibility and the effectiveness of the extended MPC scheme. Results indicate that the proposed method can obtain the smooth transition to follow the optimal drivable path and satisfy the lateral dynamic stability and environmental constraints, which can improve the path following quality for better ride comfort and road availability of autonomous vehicles.

Cite this article

Qianjie LIU , Shuang SONG , Huosheng HU , Tengchao HUANG , Chenyang LI , Qingyuan ZHU . Extended model predictive control scheme for smooth path following of autonomous vehicles[J]. Frontiers of Mechanical Engineering, 2022 , 17(1) : 4 . DOI: 10.1007/s11465-021-0660-4

1 Introduction

Autonomous vehicles have wide application potentials in intelligent transportation and military fields [13]. In general, they require several well-designed capabilities such as motion control, path planning and following, positioning, mapping, and decision making [4]. Path following is controlling the steering action such that the vehicle can move along the reference lane path with good driving safety and riding comfort according to the motion planning and vehicle driving state feedback. To achieve path following accuracy and driving stability simultaneously, path following control should effectively improve the optimal path following performance of autonomous vehicles under different driving conditions [5]. Interaction or integration between the vehicle and the road dynamics can be very substantial to enhancing the overall performance of path following. Many classical control approaches have been deployed in the path following of autonomous vehicles, which can be classified into four categories.
(1) Geometric and kinematic control. It is based on the geometric and kinematic vehicle model. Pure pursuit and Stanley methods are the most popular geometric controllers among the existing methods, which are developed for DARPA Grand Challenge vehicles [68]. These types of controllers are suitable for smooth path following at low speeds because they neglect the vehicle nonlinear dynamics that is very crucial at high speeds [9].
(2) Model-free control. No mathematical models are used, and vehicle dynamics is treated as a black box. Simple PID, fuzzy, and neural network control methods are designed for automatic steering control based on tracking errors [1013]. Moreover, several adaptive composite control methods, such as fuzzy-PID, neural network PID, adaptive fuzzy, and neural controllers [1416], are usually applied to perform parameter adjustment.
(3) State feedback control. The controller is designed by using accurate system models and implement the rigorous operations in the explicit control laws [2,17,18]. For instance, linear quadratic regulator, sliding mode control, and H infinity control methods are several common state feedback control techniques, where Lyapunov theory is generally used to evaluate the nominal stability of such systems [19,20]. The said methods cannot effectively consider the vehicle and road constraints, resulting in actuator saturation or even instable dynamics.
(4) Model predictive control (MPC). It can predict the multistep output behaviors of nonlinear system dynamics in the future and correct the following error of the system by solving the constrained optimal control problem in the moving horizon [21,22]. It can also explicitly tackle the multiconstraint problem of MIMO systems under sophisticated dynamic environments. Therefore, it has been widely used in the path following of autonomous vehicles [23,24].
The key tasks of MPC include accurate dynamic modeling, multiconstraint design, and optimization solution. Taking the high-fidelity dynamic model as a predictive model can enhance the predictive ability for the future dynamic behaviors of the vehicle and ensure a good maneuverability. The kinematic model without tires and lateral dynamics is generally used in low-speed driving scenarios [25]. Most vehicle lateral dynamics is modeled by a bicycle model due to its simplicity and high practicability, which is beneficial to improving the discretization and prediction of MPC [26]. Full vehicle dynamic model is also applied to simulate the lateral and longitudinal forces generated on each wheel during maneuvering [27]. It can completely describe the nonlinear characteristics of vehicle slip and yaw but also greatly increases its computational cost.
To date, different error performance evaluation criteria have been considered in path following controllers to track lateral error, heading error, preview look-ahead distance, and sideslip angle error and their derivatives [17,28,29], which are focused on vehicle dynamics for improving the path following performance. The multiconstraint design is also required to tackle the boundary problem from the vehicle motion and driving environment under several unstable maneuver conditions. A wide range of lateral dynamic, actuator, and road environmental constraints are encountered in path following, which are imposed on the control and output variables [3033]. Finally, the designed MPC control objective forms an online optimization problem. Furthermore, three kinds of metaheuristic optimization algorithms are presented to solve the optimization of the MPC controllers used for mobile robot path tracking control [34]. CVXGEN, particle swarm optimization, and genetic algorithm are also used in iterative operations [3537].
However, the current path following methods of autonomous vehicles mainly focus on the lateral dynamic control, but the roll dynamics on the road bank is often overlooked. As roll motion is related to vehicle model fidelity and affects the safety and ride comfort of vehicles, roll effects should be considered in the vehicle dynamic modeling for path following problem. The road centerline derived by minimizing the lateral error and heading error is strictly tracked in most path following cases, and many control techniques have been applied to address the path planning and path following problems through two independent controllers separately [32,33,38].
Nevertheless, the reference path generated by the road centerline or path planner may not be smooth enough for or even conflicts with vehicle dynamics, which may lead to overshoot, oscillation, or even instability for autonomous vehicles following it. To solve these problems, the improved multiconstraint design and optimization solution with the integration of vehicle and road dynamics should be considered in path following operations, which are the main motivations for the research work reported in this paper. The main contributions of this paper include the following:
(1) The vehicle lateral and roll constraint design for stability and environmental envelopes is further concretely derived in MPC controller formulation.
(2) An extended MPC scheme considering the driving path features is created to tackle the interaction of vehicle and road dynamics in path following operations.
(3) The proposed scheme can obtain the smooth transition on the optimal drivable path and improve the stability and maneuverability of autonomous vehicles.
The rest of the paper is organized as follows. In Section 2, vehicle dynamic modeling for path following problems, in which road curvature and inclination are used to describe the vehicle lateral and roll dynamics accurately, is described. In Section 3, a multiconstraint MPC formulation is explicitly designed by considering sideslip and rollover envelopes, road boundaries, and actuator saturation, in which the extended MPC scheme considering driving path features is presented and differential evolution (DE) optimization algorithm is used for its solution. In Section 4, Carsim–Simulink joint simulations and hardware in the loop (HIL) test are conducted to verify the effectiveness of the proposed extended MPC path tracking method. Finally, a brief conclusion and future prospects are provided in Section 5.

2 Vehicle dynamic modeling for path following problem

Building a high-precision vehicle model is the key task to determining the overall performance of the path following controller, which requires accurately capturing the lateral and roll dynamics [39]. However, the flat road is considered in most modeling works for the path following problem, and the vehicle roll effect from the suspension system is generally ignored in most modeling works for the path following problem, which is not conducive to the future dynamic prediction and constraint exertion under complex road and high-speed conditions. The road bank must be taken to enhance the model fidelity of the path tracking controller. Figure 1 shows the schematic diagram of vehicle dynamics. The vehicle can be reduced to a two-tire configuration at the front and the rear owing to symmetric dynamics of the left and right tracks. Considering the rollover instability of the vehicle, the dynamic features of sprung mass and unsprung mass under road inclination conditions are also discussed.
Fig.1 Vehicle dynamic modeling with road bank: (a) lateral dynamics, (b) roll dynamics.

Full size|PPT slide

In the assumption of small angles, the differential dynamic equations can be described by
{m(v ˙y+vxr)=m s hφ ¨+2Fy f+2 Fy r mg φr, Iz r ˙=2Fyflf2 Fy rlr,Ix φ ¨= msgh(φ +φ r)+msh( v ˙y+vxrh φ ¨)Kφφ Cφ φ ˙, e ˙y=vxeψ+vy,e ˙ψ=r κ vx,
where m and ms are the total vehicle mass and the sprung vehicle mass, respectively, lf and lr are the distances from the center of gravity of the vehicle to the front and rear axles, respectively, Ix and Iz are the roll moment of inertia and the yaw moment of inertia, respectively, Kφ and Cφ represent the combined roll stiffness coefficient and the roll damping coefficient, respectively, h is the distance from the spring mass to the roll center, r is the vehicle yaw rate, κ is the reference road curvature, vx and vy are the longitudinal and lateral velocities in the x and y directions, respectively, Fyf and Fyr are the lateral forces for the front and rear tires, respectively, φ and φr are the roll angle of the vehicle and the road bank angle, respectively, ey and eψ are the lateral and heading errors of the vehicle, respectively, which are obtained by calculating the difference between road path and output vehicle trajectory, and g is the gravitational acceleration.
The nonlinear dynamics of the vehicle is mainly reflected in the lateral tire force; therefore, a reasonable linearization of the tire model is approximated based on a tire brush model, where the lateral forces of the front and rear tires are modeled as affine [40]:
Fyj={Cdj( αj+α hj) C0 j, αj <αhj, Cjα j, αhj α jαhj, Cdj( αjαhj)+C0j,αj> α hj,
where j = {f, r} represents the front and rear tires, αj is the tire slip angle, αhj is the saturation slip angle, and Cj, Cdj, and C0j are the equivalent tire cornering stiffness identified from experimental data or more complex tire models.
The tire slip angles can also be linearized in the assumption of small angles.
{ αf= t an 1( vy+ lfrvx )δf=vy+l f rvxδf,αr= t an 1( vylrrvx)=vy lrrvx,
where δf is the front steer angle.
The tire forces symmetrically construct three regions of operations, namely, negative saturation, linear saturation, and positive saturation. Figure 2 shows the phase plane trajectories of the nonlinear vehicle dynamics generated by Eqs. (1)–(3) under different initial conditions. Two distinct areas with unstable and stable trajectories are plotted by red and blue marks separately. The stable operating region of the tire is determined by analyzing the lateral dynamic behaviors in the phase plane. Beside the stable equilibrium at the origin, two unstable equilibria clearly appear at approximately (αf, αr) = ±(0.12, 0.15) rad.
Fig.2 Phase plane trajectories of nonlinear vehicle dynamics: (a) lateral dynamic, (b) tire slip angle.

Full size|PPT slide

Let x=[vyrφ φ ˙ey eψ] T be the state vector, and δf is the control input u. In addition, γ = [φr κ]T is considered the external disturbance input. A constant vehicle longitudinal velocity is also assumed. As a result, the vehicle model with lateral and roll dynamics is described below:
M x ˙= Amx+Bumu+ B v mγ,
where M is system matrix, and
M=[ m0 00msh 00I z0 0 0 00 0 10 0 0msh00Ix+m s h20 00 000100000 01],
Am is system state matrix, Bum is system input matrix, and Bvm is system disturbance matrix.
Equation (4) can be further rewritten as
x ˙= Acx+Bucu+ B v cγ,
where Ac is state matrix, Ac = M−1Am, Buc is input matrix, Buc = M−1Bum, and Bvc is disturbance matrix, Bvc = M−1Bvm.
The following output map for yaw angle and lateral position states are defined.
y=C c x=[ 000010 000001]x,
where y is output vector, and Cc is output matrix.

3 Design of MPC-based path following control

In this section, combined with the vehicle model captured the lateral and roll dynamics, the MPC controller formulation with multiple constraints is first performed to deal with the prediction behavior, optimization, and constraint envelopes under the road bank and curvature disturbances. Then, an extended MPC-based control scheme considering the integration of vehicle and road dynamics is presented for the smooth path following of autonomous vehicles. Figure 3 shows a block diagram of the extended MPC-based control scheme.
Fig.3 Extended MPC-based control scheme.

Full size|PPT slide

3.1 MPC controller formulation

Considering the constantly changing road surface and traffic circumstances, a remarkable process in the MPC controller design for path following is to predict the future behaviors of the vehicle at each sample time Ts correspondingly. The path following problem is formulated by using a discrete-time prediction model and solved in receding horizon. Equations (5) and (6) need to be transformed into a discrete state–space model at sample time. The discrete matrices for state variables, control input, external inputs, and output variables are described by
{ A d =eA c Ts, Bud=k Ts (k+ 1)Tse Ac[(k+ 1) Tsτ]dτ B u c, Bvd= kTs (k+ 1) TseA c [(k+ 1)Tsτ]dτ B v c .
where Ad is discrete state matrix, Bud is discrete input matrix, Bvd is discrete disturbance matrix, and k is discrete time.
Therefore, the discrete state space model is obtained as
{ x( k+1)=Ad x(k) +B u du(k)+B v dγ(k) , y(k)=Cd x(k) ,
where Cd is discrete output matrix, and Cd = Cc.
An integrator is introduced to reduce the control error and enhance the precise execution. The control input of the path following controller is obtained by the incremental form u.
{ ξ (k+1)= Aξ(k) +B u Δu(k) +B v γ (k),y(k) =Cξ(k),
where augmented state vector ξ (k)= [x(k) u( k 1)] and input increment Δu(k) =u(k)u(k1), A is augmented state matrix, A=[ AdBud0I], Bu is augmented input matrix, Bu=[ B u dI], Bv is augmented disturbance matrix, Bv=[ B v d0], and C is augmented output matrix, C=[ Cd0].
This work assumes that the prediction horizon is Np, the control horizon is Nc (Nc < Np), and state variable ξ(k) captures the current plant dynamic behaviors by measurement or estimation. According to the current vehicle state at time k, the future vehicle dynamics performed at times k+1 to k+Np could be predicted by
{ ξ (k+1)= Aξ(k )+ BuΔ u(k) +B v γ (k),ξ(k+ 2)=A2ξ(k )+ ABuΔu(k) +B u Δu(k+1) +ABv γ (k)+ Bv γ (k),ξ(k+ Nc)=AN cξ(k)+AN c 1 BuΔu(k) + +BuΔu(k+N c 1)+A Nc1 B v γ (k) + + Bvγ(k), ξ(k+Np)=AN pξ(k)+AN p 1 BuΔu(k) + +AN p NcBuΔu(k+N c 1) + ANp1 Bv γ (k)+ + Bvγ(k),
where ∆u(k+Nc) = ∆u(k+Nc+1) = ∙∙∙ = ∆u(k+Np−1) = 0.
Control input remains constant when control horizon exceeds the Nc time steps. By using successive substitution, u(k+Nc−1) = u(k+Nc) = ∙∙∙ = u(k+Np−1) is assumed. In this case, the output states at times k+1 to k+Np could be computed by combining Eqs. (8)–(10).
y( k+i)=Cξ(k+i) , i=1 ,2,...,Np .
Furthermore, the predicted sequences over a horizon Np and control sequences over a horizon Nc are defined as
Xp (k)= [ξ(k+1) ξ (k+2) ξ(k+ Np)], Y p (k)= [y( k+1) y(k+2)y(k+N p)], Δ Uc (k)= [Δ u(k) Δu(k+1)Δu(k+N c 1)],
where Xp(k) is Np-step state prediction at time k, Yp(k) is Np-step output prediction at time k, and ∆Uc(k) is Nc-step input variations at time k.
By combining Eqs. (11) and (12), the predicted state and output vector is computed as
{ X p (k)= A~x ξ(k )+ B~uxΔUc (k)+ B~ vxγ (k), Y p (k)= A~ ξ (k)+ B~ u ΔUc (k)+ B~ v γ (k),
where A~x, B ~ux, and B~ vx are the predicted state matrix, input matrix, and disturbance matrix, respectively, and A ~, B ~u, and B ~v are the predicted state matrix, input matrix, and disturbance matrix for output system, respectively.
Equation (13) shows the prediction function, in which the future output trajectories of the model are predicted by using current state ξ(k), input variations ∆Uc(k), and external disturbance input γ(k).
The control sequence can be obtained by solving the optimization objective. The goal of autonomous vehicle lateral control is to track the desired road centerline accurately in the existing research [30,38], and conventional cost function Jc applied to the above prediction model is generally selected as
Jc= i=1Np y(k+i) y r ef( k+i)Q2+ i=0 Nc1 Δu(k+i)R2+λε2,
where yref(k+i) is the desired reference output, ε is a slack variable, Q and R are the corresponding output and input weight matrices, λ is slack weight, and the large weights will provide a more important priority to perform the corresponding optimization objective. The first term describes the predictive error between the prediction and reference outputs. The second term is the appropriate control inputs for the steering system. The third item reflects the optimal relaxation ability of the system, which avoids the infeasibility solution.
By substituting predicted output Eq. (13) into cost function Eq. (14), a quadratic expression in terms of the input variations ∆Uc(k) can be obtained by
Jc=VkTHkVk+GkVk+Pk,
where Hk, Vk, Gk, and Pk are quadratic transformation matrices that converts Eq. (14) into Eq. (15).
The output deviation E(k) of the free response system in the predictive domain is equivalently arranged into
E(k) =A~ξ(k) +B ~v γ (k) Yref( k) ,
where Yref is reference output matrix.
Considering the limited work conditions encountered in the control input, control increment, and predictive output at each time step, the path following problem is enforced to satisfy the following control constraints, which are concretely described in multiconstraint design.
{Δ umin( k+i)Δu(k+i) Δ umax( k+i),i=0,1, ...,Nc1, u m in( k+i) u(k+i) u m ax( k+i),i =0,1,...,Nc1, y m in( k+i) y(k+i) y m ax( k+i),i=1,2, ...,Np,
where ∆umin and ∆umax are minimum and maximum of input increment, respectively, umin and umax are minimum and maximum of control input, respectively, and ymin and ymax are minimum and maximum of output vector, respectively.
To solve the control increment of the objective function in the moving control horizons, the following optimization problem is performed by
m inΔUc,εJc (ξ(k), u(k) ,ΔUc (k), ε ), s .t .C1Δ δ f minΔ Uc (k) C1Δ δ f max, C1δf minC1u( k1 )+ C2ΔUc (k) C1δ f max, C3yh minA ~ h ξ (k)+ B~ u hΔUc (k)+ B~ v hγ(k) C3yh max, C4ys minεA ~ s ξ (k)+ B~ u sΔUc (k)+ B~ v sγ(k) C4y s max+ε ,
where C1=[111]N c ×1, C2=[1 0 0110111]N c × Nc, C3= [ I Ny h INy hINyh]Np×1, C4=[I Nys IN y s I Nys] Np×1, δfmin and δfmax are the minimum andmaximum of steer angle in the steering actuator, ∆δfmin and ∆δfmax are the minimum and maximum of steer angle increment, yhmin and yhmax are the lower and upper bounds on the output variables to be hard constrained, ysmin and ysmax are the lower and upper bounds on the output variables to be soft constrained, A~h, B~uh, and B~ vh are coefficient matrices for hard constraints, and A~s, B~us, and B~ vs are coefficient matrices for soft constraints, which can be computed by the defined multiple constraints.
The generated optimization problem consisting of input variations ∆Uc(k) and the linear constraints can be solved by
m inΔUc,ε[ Δ U c (k)ε]T Hk[ Δ Uc (k)ε]+ Gk[ Δ Uc (k)ε], s.t .[ C1Δ δ f min 0][ Δ U c (k)ε] [C1Δ δ f maxM],[ C20 C 2 0 B~ u h0 B ~ uh0 B~us1 B~ u s1][ ΔUc (k)ε][ C1δ f max C1u( k1 )C1δ f max+C1u( k1 )C3yh max A~ h ξ (k) B~ v hγ(k) C3yh min+A ~ h ξ (k)+ B~ v hγ(k) C4y s max A~ s ξ (k) B~ v sγ(k) C4ys min+A ~ s ξ (k)+ B~ v sγ(k)] .

3.2 Multiconstraint design for stability and environmental envelopes

Stability and environmental envelopes in the above optimization problem are further given by the multiconstraint design. Given that an excessively large tire sideslip angle will emerge the vehicle slip instability shown in Fig. 2, the tire slip angle imposes a limit angle αt within the stable region to avoid vehicle slip and ensure the effectiveness of the linearized tire model. Vehicle slip is generally characterized by the envelopes of the lateral velocity and the yaw rate, which can be derived from the saturation of the rear tires [35]. Therefore, the rear tire sideslip angle must satisfy the following constraint by combining Eq. (3).
vy lrrvx α t .
An additional bound on yaw rate provides the maximum steady-state condition of Eq. (1), and this requires the yaw rate to yield
r +g vx φ r min( Cf αt(1+lf/ lf lrl r)m vx, Cr αt(1+lr/ lr lfl f)m vx)=R .
To couple the constraint Eqs. (20) and (21) with the optimization problem, it is rewritten as
E1 ξ(k )+ F1γ(k) M1,
where E 1= [ 1vxl r vx0 0000 10000], F1=[ 00 g vx0], and M1=[ αtR].
Equation (22) can ensure that the vehicle meets the slip constraint and allow the sideslip angle to enter the nonlinear region for a short time to reduce vehicle instability. Moreover, producing a large lateral acceleration in high-speed emergency obstacle avoidance is easy, resulting in vehicle rollover. Thus, the rollover stability of the vehicle should also be considered in lateral motion control. Lateral load transfer rate (LTR) is a widely used rollover index that can dynamically represent tire grounding force and lateral motion states, and rollover warning and active safety control based the LTR are generally used in anti rollover control [41]. However, conducting emergency control when the vehicle is about to roll over is often very late. Considering the critical unstable dynamics of high-speed vehicles, a rollover constraint is applied to prevent rollover risk, and the rollover threshold is limited within a suitable range to prevent the tires from rolling off the ground. Estimating LTR is very difficult because normal force sensors are expensive. It is difficult to be accurately captured in the complex road surface. An alternative equivalent expression for LTR, which depends on the roll states and the vehicle parameters, can be obtained by
LTRd=2Tr[ msgh(φ +φ r)+msh(v ˙y+vxr hφ ¨) Ixφ ¨mg],
where LTRd is the equivalent LTR, andTr is the vehicle track.
To prevent rollover hazard, the rollover characterization index is constrained by
LTRdmaxLTR d LTR dm ax,
where LTRdmax is the maximum of LTRd.
Using the state space model, the rollover constraint matrix is expressed as
E2 ξ(k )+ Ed2ξ ˙ (k)+ F2 γ(k ) M2,
where E 2= [02 msh vxT r mg2 mshT r m 000], E d 2= [ 2mshTrmg 00 ms h2IxTrmg 00], F2=[ 00 2 mshT r m 0], and M2 = LTRdmax.
Considering the obstacles and the road boundary, a collision-free optimal trajectory should be generated in a given feasible driving area. According to the envelope dimension of the vehicle and the drivable road region in the path following problem, lateral road environment envelopes on vehicle trajectory can expressed as a combination of the lateral deviation thresholds at current time k, which is represented as time-varying constraints on ey.
{ ey lk=ey m axk 12wd s , eyrk=ey m ink+ 12w+ds,
where eyl represents the left road boundary, eyr represents the right road boundary, w is the vehicle width, and ds is the desired minimum distance for the vehicle to an obstacle or road edge.
Further considering the vehicle length for feasible longitudinal road region, the vehicle envelope constraint is represented by the lateral position of the front and rear axles, which should satisfy the following conditions [30].
{ yF= y0+l f (ψ +β), yR= y0 lr(ψ+β),
where yF and yR are the lateral position of the front and rear axles, respectively, y0 is the lateral position of the vehicle origin, ψ is the yaw angle, and β is the vehicle sideslip angle.
By combining Eqs. (26) and (27), the feasible road region obtained from the lateral and longitudinal vehicle envelopes is given by
En vminey kEnvmax,
where Envmin and Envmax are minimum and maximum feasible road regions, and Envmin=max( ey rklf (ψ +β), ey rk+l r (ψ +β)) and Envmax=min( ey lk lf (ψ +β), eylk+ lr (ψ +β)).

3.3 Extended model predictive steering control scheme

The conventional MPC controller with quadratic programming cannot solve the complicated path following optimization problem, which considerably limits the flexible application of MPC. To solve this problem, an extended MPC scheme is proposed here for the smooth steering control of autonomous vehicles. Considering stability and comfort in the control process, three indictors are deployed to generate a smooth path following performance:
(1) Shortest tracking path: The path length is as small as possible, which is conducive to save driving time and fuel consumption.
(2) Minimum curvature: It requires optimal path smoothness to avoid the occurrence of dangerous conditions caused by transient large curvature path.
(3) Path heading following: It makes the tracked path heading consistent with the road centerline to reduce path deviation. Thus, the optimized cost function Je is designed by
Je=i=1 Np1 Γκρ( k+i)2+ i=1Np1 ΓsS(k+i)2+ i=1 Np1 Γψ eψ( k+i)2+ i=0 Nc1 ΓuΔu(k+i)2,
where ρ(k+i) is the driving path curvature of the Np-step at current time k, S(k+i) is the driving path length of the Np-step at current time k, and Гκ, Гs, Гψ, and Гu are the corresponding weight factors.
Given the global vehicle pose (x0, y0, ts) as transient variables in the time step, the corresponding path curvature and length in the prediction step can be solved analytically by the prediction model Eq. (13). The global vehicle pose is added as
{y ˙0= vxsinψ +v ycosψ,x ˙0=vxcosψ vysinψ.
To enhance the road environment envelope in the optimized cost for drivable road region, the constraint Eq. (28) is converted into a penalty function J( Env) of feasible road region in this paper.
J(Env) =i=1 Np1{ min( ey k+iEnv m in,0)2 + max( ey k+iEnvmax,0)2}.
Combined with Eqs. (8)–(13) and (20)–(31), a new prediction model (Ap, Bp, Cp, Dp) is formed, and the extended path following optimization problem can be described as
m inΔUc (k)Je+σ J( Env),
subject to corresponding constraint Eqs. (20)–(28) and
ξp (k +1)=A p ξp (k)+ BupΔu(k) +B v pγ(k) ,
ρ(k+i)=x ˙0(k+i) y ¨0(k+i)y ˙0 (k+i)x ¨0 (k+i)[(x ˙0(k+i))2+ (y ˙0 (k+i))2 ] 32,
S( k+i)= Δx(k+i) 2+ Δy (k +i) 2,
eψ( k+i)=ψ( k+i) ψ ref( k+i),
where σ is penalty factor, ∆x and ∆y are the distance between two adjacent points in x and y directions, and ψref is the reference heading angle.
The integration and multiple constraints of vehicle and road dynamics makes the path following problem extend the derivation of intermediate variables and further non convex. To implement the extended MPC controller of path following for iteratively searching the large spaces of candidate solutions, a DE algorithm with strong robustness is used to solve the path following optimization problem, which is suitable for accelerating execution time. DE algorithm can iteratively process the constrained optimization problem by mutation, crossover, and selection operations [42]. Figure 4 describes the optimization procedure of DE algorithm in path following, and the detailed process is given as follows.
Fig.4 Optimization procedure of DE algorithm in path following.

Full size|PPT slide

Step 1. The control increment for the front wheel steering angle is described as an encoding particle, and the initial population of control increments is randomly generated within the restrained range. In this operation, several parameters are configured, including population size, mutation weight, and crossover probability. The population is randomly initialized by
{Wi Δ δ f ij=Δδf min+rand× (Δδf maxΔδf min)}, (i=1 ,2,...,P;j=1 ,2,...,D),
where Wi is the initial population, Δδfij is the initial individual, rand is a random number with a uniform probability distribution in [0, 1], P is the population size, and D is the number of the objective control inputs, which are generated for each control horizon in the optimization problem, that is, D = Nc.
Step 2. The optimization problem defined in Eq. (32) with multiple constraints is regarded as the major fitness function. As mentioned, the path curvature, length, and heading error in the prediction horizon can be calculated by using the obtained initialized population in the control step. The objective fitness function value corresponding to every reshaped population is recorded.
Step 3. A differential variable is generated by subtracting two different random objective control inputs. Mutation and crossover operations are conducted to generate new experimental input individuals in the search space by updating the differential variables. Furthermore, self-adaptive robust operations are considered in this process. For each objective control input vector, the corresponding differential mutation is expressed by
ViG+1=W r1G +η r ( W r2G Wr3G ),
where ViG+1 is the (G+1)th-generation mutation control inputs, the individual serial numbers r1, r2, and r3 are different and randomly generated, and ηr is a robust mutation factor.
To enhance the diversity of the population, a crossover operation between initial individual Δδfij and its mutation vij is also performed as follows:
ui jG +1= { vijG+1,ifrandCR o rj=rd,Δ δ fijG,otherwise,
where uij is the new experimental input individual, CR is the crossover probability within [0, 1], and rd is a randomly generated integral number in [1, 2, ..., D].
Step 4. The fitness of each control input in the updated population is evaluated. The experimental control input sequence Ui is selected as the offspring when its fitness value is better than that of the objective input individual. Otherwise, the objective input individual remaining becomes the offspring. The evolution is repeated for a fixed iteration or ended when the search converges with a given accuracy. The best input individual is used to work out the optimization problem. The selection method is
WiG+1={ UiG+1,iff( UiG+1)< f( WiG), WiG,otherwise,
where f is the fitness function defined in Eq. (32).
Step 5. The solved MPC controller is implemented with the first control input of the generated optimal objective sequence under the constraints. Then, the optimization problem is solved forwardly by updating these control input sequences.

4 Experimental results and analysis

To validate the feasibility and the effectiveness of the developed extended MPC controller for path following of autonomous vehicles, Carsim–Simulink cosimulation and HIL test are performed under different driving conditions. The parameters characterizing the vehicle dynamic model are listed as below:
m = 1530 kg,
ms = 1370 kg,
Ix = 671.3 kg·m2,
Iz = 2315.3 kg·m2,
lf = 1.11 m,
lr = 1.67 m,
h = 0.52 m,
Tr = 1.55 m,
Kφ = 183791 N·m/rad,
Cφ = 4904 N·m·s/rad,
Cf = 66800 N/rad,
Cr = 62700 N/rad.
The main parameters used in the MPC controller by trial adjustment are shown as below:
Ts = 0.02 s,
Np = 20,
Nc = 5,
δmin = –0.52 rad,
δmax = 0.52 rad,
δmin = –0.12 rad/s,
δmax = 0.12 rad/s,
eymin = –3 m,
eymax = 3 m,
αt = 0.1 rad,
LTRdmax = 0.7,
Г κ = 1500,
Гs = 20,
Гψ= 350,
Гu = 250.

4.1 Simulation results

In the simulation test, the nominal path and the road surface are presented in Fig. 5, which shows the road curvature, heading, and road elevation in the path station. To verify the accuracy of rollover index LTRd, LTR is used as nominal reference object. Figure 6 presents a comparison of LTR, LTRd, and LTRd without considering φr at the speed of 30 m/s, which clearly shows an evident error in LTRd without considering φr. As the rollover index LTRd can effectively capture LTR, it can be effectively used to characterize the rollover dynamics of the vehicle and further exert its constraint.
Fig.5 Nominal path and road bank: (a) nominal path, (b) road bank.

Full size|PPT slide

Fig.6 Comparison of LTR, LTRd, and LTRd without considering φr.

Full size|PPT slide

Figure 7 presents the simulation results of path following for considering multiple constraints and road surface in the conventional MPC controller, where the road centerline is considered the expected trajectory. More specifically, Fig. 7(a) shows that the controller with multiple constraints can effectively track the nominal path, and the path following controller with multiple constraints and road surface can remarkably further reduce tracking error. The controller without multiple constraints fails to follow the drivable path after approximately 28 s. Figure 7(b) shows that yaw angle following performs well on the controller with multiple constraints. Moreover, Fig. 7(c) shows that the lateral stability can effectively avoid the occurrence of unstable sideslip behavior in three turns. The apparent oscillating steering angle appears in Fig. 7(d) due to unstable road dynamics and mismatched trace path. To eliminate the overshoots of the conventional MPC controller, smooth path following is further considered in the proposed extended MPC controller.
Fig.7 Simulation results of path following for considering multiple constraints and road bank: (a) lateral error, (b) yaw error, (c) lateral velocity versus yaw rate, and (d) steer angle.

Full size|PPT slide

Table 1 summarizes the root mean square (RMS) and maximum (Max) values of path following for conventional MPC controller. The controller considering constraints and road surface achieves smaller RMS and Max values of path following than other controllers, and its maximum yaw error and lateral velocity are reduced by 82.9% and 43.8%, respectively, compared with the unconstrained one. Owing to the failure path following behavior shown in the unconstrained controller, its lateral error and steer angle are also far beyond the feasible tracking area.
Tab.1 RMS and Max values of path following for conventional MPC controller
Controller Lateral error/m Yaw error/rad Yaw rate/(rad·s−1) Lateral velocity/(m·s−1) Steer angle/rad
RMS Max RMS Max RMS Max RMS Max RMS Max
With φr and constraint 0.0170 0.1019 0.0078 0.0340 0.1487 0.3347 0.3688 1.1068 0.0245 0.0915
Without φr 0.0249 0.1415 0.0080 0.0356 0.1498 0.3262 0.3707 1.1351 0.0247 0.1001
Without constraint 0.0283 0.1990 0.1672 0.3529 0.4438 1.9691
Figure 8 presents a comparison of the simulation results of path following for conventional and extended MPC controllers at 20 m/s, where road boundary is the enveloped region by the left and right road boundaries eyl and eyr in Eq. (26). Although the conventional MPC controller can make it track the reference road centerline, the steering angle and the path curvature show more prominent oscillations than that of the extended MPC controller when the vehicle is turning because the road dynamic information is not considered in the conventional MPC controller, which will scarify the comfort and the stability of the vehicle. The autonomous vehicle ignores road trafficability by tracking the lane center line and is not available for the driver’s driving habits. To improve the quality of path following, the driving path information of the vehicle should be considered in the path following optimization problem.
Fig.8 Simulation results of path following for conventional and extended MPC controllers at 20 m/s: (a) path following, (b) steer angle, (c) path curvature, (d) yaw rate, and (e) rollover threshold.

Full size|PPT slide

More specifically, Fig. 8(a) shows that the vehicle’s driving path is constantly adjusted in the drivable road region instead of strictly tracking the road center line. The MPC controller drives the vehicle to perform the optimal transition steering behavior when the vehicle is making a turn and keep heading straight following the driving deviation from the lane center line at the end of the turn, which can preferably complete the transition steering of the next turning.
Figures 8(b) and 8(c) show that the steering and curvature profiles for the extended MPC controllers are much smoother than the ones for the conventional MPC controller, especially at road corners. Its curvature is remarkably improved compared with the original one, which achieves the flexible smooth path following. In Figs. 8(d) and 8(e), the path following of the extended MPC has a higher stable yaw rate and a smaller rollover threshold than those of the conventional MPC, which can satisfy the better yaw stability and antirollover performance. Consequently, the extended MPC scheme can make full use of vehicle maneuverability and provide a more satisfactory path following performance than the conventional one.
Figure 9 shows the comparisons of the path following results for conventional and extended MPC controllers at 30 m/s. The steering performance of both schemes deteriorate at different degrees. Such situations may occur due to poor lateral stability and more sensitive steering at high-speed corners. The extended MPC scheme can achieve the acceptable path following performance and provide more superior path curvature and stable lateral dynamics. Figure 9(a) shows that the large transition steering behavior is required at high speed, which results in an increase in the steering angle performed by the autonomous vehicle, as displayed in Fig. 9(b). These results further confirm that the extended MPC controller can ensure the good effectiveness and the robustness of smooth path following for autonomous vehicles.
Fig.9 Simulation results of path following for conventional and extended MPC controllers at 30 m/s: (a) path following, (b) steer angle, (c) path curvature, (d) yaw rate, and (e) rollover threshold.

Full size|PPT slide

The RMS and Max values of path following for conventional and extended MPC controllers under different speeds are illustrated in Table 2. The extended MPC controller can effectively reduce the RMS and Max values of path following. The RMS values of the extended MPC controller are slightly reduced compared with the traditional one, whereas its Max values show a substantial reduction. This result is not unexpected because the extended controller plays a major role in suppressing the undesirable oscillation for path following. According to the maximum steering angle and path curvature, the smoothness of the extended controller is improved by more than 40% than that of the traditional one, and its lateral stability represented by yaw rate and LTRd is improved by approximately 19% on average.
Tab.2 RMS and Max values of path following for conventional and extended MPC controllers
Controller Speed/(m·s−1) Steer angle/rad Path curvature/(10−3 m−1) Yaw rate/(rad·s−1) LTRd
RMS Max RMS Max RMS Max RMS Max
Extended MPC 20 0.0173 0.0319 4.73 8.3 0.0949 0.1658 0.0822 0.1941
30 0.0235 0.0464 4.74 8.5 0.1447 0.2795 0.2201 0.4778
Conventional MPC 20 0.0180 0.0557 4.86 16.8 0.0965 0.2112 0.0977 0.2739
30 0.0247 0.0915 4.92 14.8 0.1510 0.3347 0.2565 0.5281

4.2 HIL experimental results

For further verifying the performance of the proposed extended MPC scheme, an HIL test experiment based on Simulink and dSPACE system is carried out, as shown in Fig. 10. It mainly consists of the dSPACE system, steering system with servo motor, loading cylinder, sensors, monitor, and power supply. In this system, the hardware bench system realizes the real-time steering control of the autonomous vehicle and feedback of the sensor measurement data, and the loading cylinder is used to simulate steering resistance. The dSPACE system compiles and embeds the vehicle–road dynamic model and control algorithm. Then, different strategies are generated to calculate steering angle command and send it to the actuator driver in real time.
Fig.10 HIL test system based on Simulink and dSPACE system.

Full size|PPT slide

A curve test driving at 30 m/s is performed on the HIL system, and the experiment results are shown in Fig. 11. A comparison is made in terms of path following, steering angle, path curvature, and lateral yaw dynamics with conventional and extended MPC controllers. The vehicle model and control parameters configured in the HIL system are not completely consistent with those set in the simulation environment. The steering angle is fed back by the sensor measurement, and other results are output by the embedded vehicle dynamics model, where the disturbance response mainly appearing in the steering angle is caused by communication delay and measurement noise. Figure 11(a) shows that the path trajectory with the conventional MPC control is strictly limited in the road centerline, which fails to use other drivable road areas comprehensively, whereas the extended MPC control strategy enables the autonomous vehicle to follow the optimal route within the road lane. Figures 11(b) and 11(c) show that the extended MPC control displays a smoother path following performance than the conventional one. Figure 11(d) shows that the lateral dynamics for the extended MPC control is maintained in a smaller envelope. By contrast, the oscillation of the conventional MPC controller at the corner conduces to the divergence of lateral dynamics. Apparently, the extended MPC control is more suitable for smooth path following of autonomous vehicles.
Fig.11 Experiment comparison results performed in HIL system: (a) path following, (b) steer angle, (c) path curvature, and (d) lateral velocity versus yaw rate.

Full size|PPT slide

5 Conclusions

An extended MPC controller, which considers vehicle and road dynamics, is developed in this paper to find autonomously the optimal driving path in the drivable road regions for smooth path following of autonomous vehicles. A vehicle dynamic model is created for accurately capturing the lateral motion and roll dynamics of the vehicle. Multiple constraints are concretely derived in the MPC controller formulation for stability and environmental envelopes. Then, the extended MPC method solved by DE optimization algorithm is proposed and takes the vehicle and road dynamics as a combined problem for smooth path following. Simulation and HIL experiments are conducted to verify the feasibility and the performance of the proposed extended MPC scheme. The results show that the proposed approach can obtain the optimal path following on the drivable road and satisfy the lateral dynamic stability constraints of the vehicle. It can perform smooth transition steering on the road turns and maintain the path heading safely instead of crossing the road area in the road straight section. Smooth path following performance can be further enhanced for broader road drivable area and in coordination with human drivers. Moreover, the proposed method can effectively address the undesirable oscillations generated in constant speed turning. In the future work, to handle the complex turn conditions under variable speed operation, the longitudinal control with time-varying velocity and nonlinear tire model will be considered for the path following control with sharp turn behavior. Real vehicle experiments will also be conducted in a dynamic environment.

Nomenclature

Cdj, Cj, C0j Cornering stiffness coefficient
Cφ Combined roll damping coefficient
CR Crossover probability within [0, 1]
ds Desired minimum safe distance
D Objective control input size
ey Lateral error
eyl Left road boundary
eyr Right road boundary
eψ Heading error
Envmax Maximum feasible road region
Envmin Minimum feasible road region
f Fitness function
Fyf Lateral force of front tire
Fyr Lateral force of rear tire
g Gravitational acceleration
h Distance from the spring mass to the roll center
Ix Roll moment of inertia
Iz Yaw moment of inertia
k Discrete time
Kφ Combined roll stiffness coefficient
lf Distance from center to front axle
lr Distance from center to rear axle
LTR Lateral load transfer rate
LTRd Equivalent LTR
LTRmax Maximum value of LTRd
m Total vehicle mass
ms Sprung vehicle mass
Nc Control horizon
Np Prediction horizon
P Population size
r Vehicle yaw rate
rand Random number with a uniform probability distribution in [0, 1]
rd Randomly generated integral number in [1, 2, ..., D]
S Driving path length
Tr Vehicle track
Ts Sample time
uij New experimental input individual
vij Mutation individual of Δδfij
vx Longitudinal velocity
vy Lateral velocity
w Vehicle width
x0 Longitudinal position of vehicle origin
x Distance between two adjacent points in x direction
y0 Lateral position of vehicle origin
yF Lateral position of front axle
yR Lateral position of rear axle
y Distance between two adjacent points in y direction
αhj Saturation slip angle
αj Tire slip angle
αf Front tire slip angle
αr Rear tire slip angle
αt Limit tire slip angle
β Vehicle sideslip angle
δf Front steer angle
δfmax Maximum value of δf
δfmin Minimum value of δf
ε Slack variable
φ Roll angle of vehicle
φr Road bank angle
ηr Robust mutation factor
κ Reference road curvature
λ Slack weight
σ Penalty factor
ρ Driving path curvature
ψ Yaw angle
ψref Reference heading angle
δfmax Maximum of steer angle increment
δfmin Minimum of steer angle increment
Δδfij Initial individual
Г κ Weight factor of ρ
Гs Weight factor of S
Гu Weight factor of u
Гψ Weight factor of eψ
A Augmented state matrix
Ac State matrix
Ad Discrete state matrix
Am System state matrix
A~ Predicted state matrix for output system
A ~h, B~ uh, B~ vh Coefficient matrices for hard constraints
A ~s, B~ us, B~ vs Coefficient matrices for soft constraints
A ~x Predicted state matrix
Bu Augmented input matrix
Buc Input matrix
Bud Discrete input matrix
Bum System input matrix
Bv Augmented disturbance matrix
Bvc Disturbance matrix
Bvd Discrete disturbance matrix
Bvm System disturbance matrix
B ~u Predicted input matrix for output system
B ~ux Predicted input matrix
B ~v Predicted disturbance matrix for output system
B ~vx Predicted disturbance matrix
C Augmented output matrix
Cc Output matrix
Cd Discrete output matrix
E Output deviation
E1, E2, Ed2 Constraint state matrix
F1, F2 Constraint input matrix
Hk, Vk, Gk, Pk Quadratic transformation matrices
Jc Conventional cost function
Je Optimized cost function
J(Env) Penalty function of feasible road region
M System matrix
M1, M2 Constraint output matrix
Q Output weight matrix
R Input weight matrix
u Control input
umax Maximum value of u
umin Minimum value of u
Ui Experimental control input sequence
u Input increment
umax Maximum value of u
umin Minimum value of u
Uc Input variation
Vi Mutation control input
Wi Initial population
x State vector
Xp State prediction
y Output vector
yhmax Hard constraint of upper output bound
yhmin Hard constraint of lower output bound
ymax Maximum value of y
ymin Minimum value of y
yref Desired reference output
ysmax Soft constraint of upper output bound
ysmin Soft constraint of lower output bound
Yp Output prediction
Yref Reference output matrix
γ Disturbance input
ξ Augmented state vector

Acknowledgements

This work was funded by the National Natural Science Foundation of China (Grant No. 52075461), the Key Project in Science and Technology Plan of Xiamen, China (Grant No. 3502Z20201015), and the Science and Technology Plan of Fujian Province of China (Grant No. 2021H6019).
1
LiK Q, GaoF, LiS E, ZhengY, GaoH. Robust cooperation of connected vehicle systems with eigenvalue-bounded interaction topologies in the presence of uncertain dynamics. Frontiers of Mechanical Engineering, 2018, 13( 3): 354– 367

DOI

2
SpielbergN A, BrownM, KapaniaN R, KegelmanJ C, GerdesJ C. Neural network vehicle models for high-performance automated driving. Science Robotics, 2019, 4( 28): aaw1975–

DOI

3
BadueC, GuidoliniR, CarneiroR V, AzevedoP, CardosoV B, ForechiA, JesusL, BerrielR, PaixãoT M, MutzF, dePaula Veronese L, Oliveira-SantosT, DeSouza A F. Self-driving cars: a survey. Expert Systems with Applications, 2021, 165 : 113816–

DOI

4
GonzálezD, PérezJ, MilanésV, NashashibiF. A review of motion planning techniques for automated vehicles. IEEE Transactions on Intelligent Transportation Systems, 2016, 17( 4): 1135– 1145

DOI

5
GuoJ H, LuoY G, LiK Q, DaiY. Coordinated path-following and direct yaw-moment control of autonomous electric vehicles with sideslip angle estimation. Mechanical Systems and Signal Processing, 2018, 105 : 183– 199

DOI

6
ThrunS, MontemerloM, DahlkampH, StavensD, AronA, DiebelJ, FongP, GaleJ, HalpennyM, HoffmannG, LauK, OakleyC, PalatucciM, PrattV, StangP, StrohbandS, DupontC, JendrossekL E, KoelenC, MarkeyC, RummelC, van NiekerkJ, JensenE, AlessandriniP, BradskiG, DaviesB, EttingerS, KaehlerA, NefianA, MahoneyP. Stanley: the robot that won the DARPA Grand Challenge. Journal of Field Robotics, 2006, 23( 9): 661– 692

DOI

7
SniderJ M. Automatic steering methods for autonomous automobile path tracking. Dissertation for the Doctoral Degree. Pittsburgh: Carnegie Mellon University Pittsburgh, 2009

8
Morales J, Martínez J L, Martínez M A, Mandow A. Pure-pursuit reactive path tracking for nonholonomic mobile robots with a 2D laser scanner. EURASIP Journal on Advances in Signal Processing, 2009, (1): 935237

9
ChenQ Y, SunZ P, LiuD X, LiX. Ribbon model based path tracking method for autonomous ground vehicles. Journal of Central South University, 2014, 21( 5): 1816– 1826

DOI

10
ZhangC Z, HuJ F, QiuJ B, YangW, SunH, ChenQ. A novel fuzzy observer-based steering control approach for path tracking in autonomous vehicles. IEEE Transactions on Fuzzy Systems, 2019, 27( 2): 278– 290

DOI

11
MarinoR, ScalziS, NettoM. Nested PID steering control for lane keeping in autonomous vehicles. Control Engineering Practice, 2011, 19( 12): 1459– 1467

DOI

12
MohammadzadehA, TaghavifarH. A robust fuzzy control approach for path-following control of autonomous vehicles. Soft Computing, 2020, 24( 5): 3223– 3235

DOI

13
HuangW, WongP K, WongK I, VongC M, ZhaoJ. Adaptive neural control of vehicle yaw stability with active front steering using an improved random projection neural network. Vehicle System Dynamics, 2021, 59( 3): 396– 414

DOI

14
ZhaoP, ChenJ J, SongY, TaoX, XuT, MeiT. Design of a control system for an autonomous vehicle based on adaptive-PID. International Journal of Advanced Robotic Systems, 2012, 9( 2): 44–

DOI

15
AbdelhakimG, AbdelouahabH. A new approach for controlling a trajectory tracking using intelligent methods. Journal of Electrical Engineering & Technology, 2019, 14( 3): 1347– 1356

DOI

16
GuoJ H, HuP, LiL H, WangR. Design of automatic steering controller for trajectory tracking of unmanned vehicles using genetic algorithms. IEEE Transactions on Vehicular Technology, 2012, 61( 7): 2913– 2924

DOI

17
KritayakiranaK, GerdesJ C. Using the centre of percussion to design a steering controller for an autonomous race car. Vehicle System Dynamics, 2012, 50(sup1 S1): 33– 51

18
KapaniaN R, GerdesJ C. Design of a feedback-feedforward steering controller for accurate path tracking and stability at the limits of handling. Vehicle System Dynamics, 2015, 53( 12): 1687– 1704

DOI

19
ChenW W, ZhaoL F, WangH R, HuangY. Parallel distributed compensation/H control of lane-keeping system based on the Takagi-Sugeno fuzzy model. Chinese Journal of Mechanical Engineering, 2020, 33( 1): 61–

DOI

20
CalzolariD, SchurmannB, AlthoffM. Comparison of trajectory tracking controllers for autonomous vehicles. In: Proceedings of IEEE 20th International Conference on Intelligent Transportation Systems (ITSC). Yokohama: IEEE, 2017

21
FalconeP, BorrelliF, AsgariJ, TsengH E, HrovatD. Predictive active steering control for autonomous vehicle systems. IEEE Transactions on Control Systems Technology, 2007, 15( 3): 566– 580

DOI

22
FalconeP, BorrelliF, TsengH E, AsgariJ, HrovatD. Linear time-varying model predictive control and its application to active steering systems: stability analysis and experimental validation. International Journal of Robust and Nonlinear Control, 2008, 18( 8): 862– 875

DOI

23
PengH N, WangW D, AnQ, XiangC, LiL. Path tracking and direct yaw moment coordinated control based on robust MPC with the finite time horizon for autonomous independent-drive vehicles. IEEE Transactions on Vehicular Technology, 2020, 69( 6): 6053– 6066

DOI

24
SuhJ, ChaeH, YiK. Stochastic model-predictive control for lane change decision of automated driving vehicles. IEEE Transactions on Vehicular Technology, 2018, 67( 6): 4771– 4782

DOI

25
LiZ J, DengJ, LuR Q, XuY, BaiJ, SuC Y. Trajectory-tracking control of mobile robot systems incorporating neural-dynamic optimized model predictive approach. IEEE Transactions on Systems, Man, and Cybernetics: Systems, 2016, 46( 6): 740– 749

DOI

26
AmerN H, ZamzuriH, HudhaK, KadirZ A. Modelling and control strategies in path tracking control for autonomous ground vehicles: a review of state of the art and challenges. Journal of Intelligent & Robotic Systems, 2017, 86( 2): 225– 254

DOI

27
SongY T, ShuH Y, ChenX B. Chassis integrated control for 4WIS distributed drive EVs with model predictive control based on the UKF observer. Science China. Technological Sciences, 2020, 63( 3): 397– 409

DOI

28
TagneG, TaljR, ChararaA. Design and comparison of robust nonlinear controllers for the lateral dynamics of intelligent vehicles. IEEE Transactions on Intelligent Transportation Systems, 2016, 17( 3): 796– 809

DOI

29
ChowdhriN, FerrantiL, IribarrenF S, ShyrokauB. Integrated nonlinear model predictive control for automated driving. Control Engineering Practice, 2021, 106 : 104654–

DOI

30
GuoH Y, LiuJ, CaoD P, ChenH, YuR, LvC. Dual-envelop-oriented moving horizon path tracking control for fully automated vehicles. Mechatronics, 2018, 50 : 422– 433

DOI

31
BealC E, GerdesJ C. Model predictive control for vehicle stabilization at the limits of handling. IEEE Transactions on Control Systems Technology, 2013, 21( 4): 1258– 1269

DOI

32
JiJ, KhajepourA, MelekW W, HuangY. Path planning and tracking for vehicle collision avoidance based on model predictive control with multiconstraints. IEEE Transactions on Vehicular Technology, 2017, 66( 2): 952– 964

DOI

33
LiX H, SunZ P, CaoD P, LiuD, HeH. Development of a new integrated local trajectory planning and tracking control framework for autonomous ground vehicles. Mechanical Systems and Signal Processing, 2017, 87 : 118– 137

DOI

34
MerabtiH, BelarbiK, BouchemalB. Nonlinear predictive control of a mobile robot: a solution using metaheuristcs. Journal of the Chinese Institute of Engineers, 2016, 39( 3): 282– 290

DOI

35
BrownM, FunkeJ, ErlienS, GerdesJ C. Safe driving envelopes for path tracking in autonomous vehicles. Control Engineering Practice, 2017, 61 : 307– 316

DOI

36
GuoH Y, LiuF, XuF, ChenH, CaoD, JiY. Nonlinear model predictive lateral stability control of active chassis for intelligent vehicles and its FPGA implementation. IEEE Transactions on Systems, Man, and Cybernetics: Systems, 2019, 49( 1): 2– 13

DOI

37
DuX X, HtetK K K, TanK K. Development of a genetic-algorithm-based nonlinear model predictive control scheme on velocity and steering of autonomous vehicles. IEEE Transactions on Industrial Electronics, 2016, 63( 11): 6970– 6977

DOI

38
GuoH Y, ShenC, ZhangH, ChenH, JiaR. Simultaneous trajectory planning and tracking using an MPC method for cyber-physical systems: a case study of obstacle avoidance for an intelligent vehicle. IEEE Transactions on Industrial Informatics, 2018, 14( 9): 4273– 4283

DOI

39
CaoH T, SongX L, ZhaoS, BaoS, HuangZ. An optimal model-based trajectory following architecture synthesising the lateral adaptive preview strategy and longitudinal velocity planning for highly automated vehicle. Vehicle System Dynamics, 2017, 55( 8): 1143– 1188

DOI

40
Di CairanoS, TsengH E, BernardiniD, BemporadA. Vehicle yaw stability control by coordinated active front steering and differential braking in the tire sideslip angles domain. IEEE Transactions on Control Systems Technology, 2013, 21( 4): 1236– 1248

DOI

41
XiaoL J, WangM, ZhangB J, ZhongZ. Vehicle roll stability control with active roll-resistant electro-hydraulic suspension. Frontiers of Mechanical Engineering, 2020, 15( 1): 43– 54

DOI

42
LiuQ J, ChenW, HuH S, ZhuQ, XieZ. An optimal NARX neural network identification model for a magnetorheological damper with force-distortion behavior. Frontiers in Materials, 2020, 7 : 10–

DOI

Outlines

/