Next Article in Journal
Two-Degrees-of-Freedom PID Control with Kalman Filter for Engraving Machine System
Next Article in Special Issue
Study of Error Flow for Hydraulic System Simulation Models for Construction Machinery Based on the State-Space Approach
Previous Article in Journal
Research on a Variable-Stiffness Joint and Its Application in Actuators
Previous Article in Special Issue
Improved Vehicle Vibration Control through Optimization of Suspension Parameters Using the Response Surface Method and a Non-Linear Programming with a Quadratic Lagrangian Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Study on Coordinated Control of Path Tracking and Vehicle Stability for Autonomous Vehicles on Low-Friction Roads

1
Department of Electronic Engineering, Korea National University of Transportation, Chungju-si 27469, Republic of Korea
2
Research Center for Electrical and Information Technology, Seoul National University of Science and Technology, Nowon-gu, Seoul 01811, Republic of Korea
*
Author to whom correspondence should be addressed.
Actuators 2023, 12(11), 398; https://doi.org/10.3390/act12110398
Submission received: 5 October 2023 / Revised: 24 October 2023 / Accepted: 24 October 2023 / Published: 25 October 2023

Abstract

:
This paper presents a comparative study on coordinated control of path tracking and vehicle stability for autonomous vehicles on low-friction roads. Generally, a path-tracking controller designed on high-friction roads cannot provide good performance under low-friction conditions. To cope with the problem, a coordinated control between path tracking and vehicle stability has been proposed to date. In this paper, three types of coordinated controllers are classified according to the controller structure. As an actuator, front-wheel steering, four-wheel steering, and four-wheel independent braking and driving are adopted. A common feature of these types of controllers is that front steering and yaw moment control are adopted as control inputs. To convert the yaw moment control into tire forces generated by combinations of multiple actuators, a control allocation method is applied. For each type, a controller is designed and simulated using vehicle simulation software. From the simulation results, a performance comparison among those controller types is carried out. Through comparison, it is shown that there are small differences among those types of controllers in terms of path tracking.

1. Introduction

Autonomous driving has been regarded as a next-generation solution for future transportation in research groups and the automotive industry since 2010 because it can prevent several types of traffic accidents and make road safety and traffic flow better [1,2,3]. From the literature survey, it is widely known that a generic modular system architecture for autonomous driving comprises detection, localization, prediction, planning, and control [2]. Among those topics, path-tracking control (PTC) has been intensively studied [4,5]. Consequently, a great deal of papers have been published in the area of PTC [4,5,6,7,8,9,10].
The aim of PTC is to make a vehicle follow a target path. For this purpose, most of the path-tracking controllers proposed to date have tried to minimize the lateral offset error (ey) and heading error (eφ), which are derived from the current vehicle position and target path [8,9]. Generally, front-wheel steering (FWS) or four-wheel steering (4WS) has been adopted as a control input or an actuator of a path-tracking controller. In this paper, 4WS means the combination of FWS and rear-wheel steering (RWS). This controller has shown good performance for path tracking on high-friction roads because this is designed and tuned under high-friction conditions. However, a PTC designed on high-friction roads cannot provide good performance under low-friction conditions [11,12,13]. This is caused by the fact that the lateral tire force of a vehicle becomes small on low-friction roads. To cope with the problem, PTC has been integrated with vehicle stability or lateral stability control. Vehicle stability control (VSC) has been intensively studied since the middle of the 1990s [14,15,16,17]. In this paper, VSC stands for yaw rate tracking and lateral stability controls, where the former tries to minimize the yaw rate error (γe) and the latter tries to keep the side-slip angle (β) as small as possible. Generally, the yaw moment control (∆Mz) has been adopted as a control input of VSC. ∆Mz is a virtual input, which should be generated by multiple actuators installed in real vehicles. For this reason, ∆Mz is to be converted into the tire forces generated by multiple actuators. This procedure is called control allocation or yaw moment distribution. Actuators used to generate ∆Mz are FWS, 4WS, RWS, and four-wheel independent steering/driving/braking (4WIS/4WID/4WIB) [18,19]. Most control allocation methods are quadratic programming (QP) with equality or/and inequality constraints [11,12,13,18,19].
Numerous papers have been published on the coordinated control of path tracking and vehicle stability on low-friction road conditions [11,12,13,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40]. These papers can be classified into three types. Figure 1 shows the control structures of these types.
The first type of coordinated control method is to use two-level controllers: upper- and lower-level controllers [13,20,21,22,23,24,25,26,27,28,29]. As a control input, steering angles and ∆Mz are calculated by controllers such as model predictive control (MPC), PID control, sliding mode control (SMC), and H control in the upper-level controllers (UPLC). Among the control inputs, steering angles and ∆Mz are in charge of PTC and VSC in low-friction roads, respectively. For controller design, most of these types of controllers adopted a state-space model with ey, eφ, and its rates. ∆Mz is converted into tire forces generated by a single actuator or a combination of multiple ones such as FWS, RWS, 4WS, 4WIS, 4WIB, and 4WID in the lower-level controller (LWLC). For this purpose, a control allocation method was adopted. The steering angles of FWS, RWS, and 4WS, and braking and traction torques of 4WID and 4WIB are calculated from the converted tire forces. The effects of combinations of multiple actuators on path-tracking performance were deeply investigated in the previous research [13]. According to this work, there are few differences among the combinations of FWS, 4WS, 4WID, and 4WIB in terms of path-tracking performance. In other words, it is quite enough to use only FWS for path tracking. Let us denote this type of method as PTC#1.
The second type of coordinated control method is to use three-level controllers: upper-, middle-, and lower-level controllers [12,30,31,32,33,34,35,36,37]. In the UPLC, the front steering angle is calculated by a driver model or a controller. As a driver model, pure pursuit and Stanley methods have been widely used to date [4,11,12]. As a controller for steering angle generation, linear quadratic regulator (LQR), MPC, and SMC have been widely adopted with the state-space model whose state variables are ey and eφ. This is identical to PTC#1 with FWS. In the middle-level controller (MDLC), a reference yaw rate (γref) is calculated from the front steering angle (δf), and a controller is designed to calculate ∆Mz, which is needed to make a vehicle follow the reference yaw rate. Most of MDLCs adopted SMC in order to calculate ∆Mz. In the previous paper, the reference yaw rate is directly generated as a control input of the UPLC, H control [34]. In the LWLC, ∆Mz is converted into tire forces generated by multiple actuators. For this purpose, a control allocation method is applied. Let us denote this type of method as PTC#2. In view of PTC#1, it seems unnecessary to use the three-level controller of PTC#2 because PTC#1 has the upper- and middle-level controllers of PTC#2 in the UPLC. In other words, to generate ∆Mz, PTC#2 requires two controllers although PTC#1 does a single controller. If the path-tracking performance of PTC#2 is not superior to that of PTC#1 on low-friction roads, then PTC#2 is useless.
The third type of coordinated control method is identical to the second one except that the reference yaw rate is obtained from a target path in the UPLC [11,38,39,40]. In this type of method, a driver model or LQR is needed to generate δf. In the MDLC of this method, SMC has been applied to calculate ∆Mz, which is needed to make a vehicle follow the reference yaw rate. In the LWLC, a control allocation method was applied to convert ∆Mz into tire forces generated by multiple actuators such as RWS, 4WID, and 4WIB. Let us denote this type of method as PTC#3. In view of PTC#1, PTC#3 has the identical problem to PTC#2 because PTC#3 has a three-level structure.
The first feature of PTC#1, PTC#2, and PTC#3 is the control allocation in LWLC, which is applied to distribute ∆Mz into tire forces generated by combinations of multiple actuators. In fact, control allocation and LWLC are synonymous. For control allocation, an optimization method with QP has been applied to date. These are quite general in vehicle stability control or lateral stability control. For this reason, PTC#1, PTC#2, and PTC3 try to integrate path tracking and vehicle stability control for driving on low-friction roads.
The second feature of PTC#2 and PTC#3 is the calculation procedure of γref. PTC#1 does not use γref. In PTC#2 and PTC#3, γref for path tracking is calculated from δf or a target path in the MDLC. In PTC#1, PTC#2, and PTC#3, δf is generated using LQR. Because δf is already determined using LQR in the UPLC, it is excluded from the control allocation procedure. In other words, δf is not determined using the control allocation procedure. For this reason, the actuators used for control allocation are RWS, 4WID, and 4WIB.
The third feature of these methods is to use a preview or lookahead function along the heading direction [12,21,22,24,34,38,40]. Generally, the larger the preview distance the smaller the steering angle. Especially, the preview distance should be larger for lateral stability on low-friction roads. On the contrary, the preview distance should be smaller for path tracking on high-friction roads [12,41]. For this reason, the preview distance should be carefully tuned for path tracking and lateral stability if it is adopted [13].
The aim of this paper is to compare three types of coordinated controllers, PTC#1, PTC#2, and PTC#3, for PTC and VSC on low-friction road conditions. As a state variable, ey, eφ, β, and γ are selected, and these variables are included in the LQ objective function for PTC and VSC. In the state-space model, a preview function is introduced. In PTC#1, four input configurations, [δf], [δf δr], [δfMz] and [δf δrMz], are selected from three control inputs, δf, δr, and ∆Mz. If ∆Mz is selected as a control input, a control allocation method is adopted to distribute it to the tire forces generated using multiple actuators. As an actuator for control allocation, 4WID and 4WIB are adopted [13,39]. To compare those controllers, a simulation is conducted using the vehicle simulation package, CarSim. The double lane change maneuver for collision avoidance is adopted as a test scenario [11,12,20,26,27,37]. From the simulation results, it is shown that there are few differences among the three types of coordinated controllers on low-friction roads.
The contributions of this paper can be condensed as follows:
  • This paper compares three types of coordinated controllers for PTC and VSC using LQR with regard to path-tracking performance. Different from the previous works, β and γ are included as state variables instead of the rates of lateral offset and heading errors. A preview function is also introduced into the state-space model. With the state vector, it is easier to include the preview function in the state-space model.
  • Through comparison, the effects of input configurations and controller structures on path-tracking performance under low-friction conditions were investigated. With the investigation, which type of controller is more effective than others for path tracking is determined.
  • From the comparison results, the best control structure and input configuration for path tracking and vehicle stability on low-friction roads is recommended. If there are little differences among those types of controllers, the simplest one is the best.
This paper consists of five sections. A design procedure of three types of coordinated controllers is described in Section 2. In Section 3, several measures for path-tracking performance are presented. The simulation is carried out and the simulation results are discussed in terms of the performance measures in Section 4. The conclusion is provided in Section 5.

2. Design of Controllers for Coordinated Control of Path Tracking and Vehicle Stability

2.1. Derivation of State-Space Model

Following the literature of PTC, a 2-DOF bicycle model is selected as vehicle one. Figure 2 shows the coordinates and variables of the bicycle model and the target path used for path tracking [12,13]. This model describes the lateral and yaw motions of a vehicle under the assumption that vx is constant. With γ and β as a state variable, the equations of motions for the model are obtained as (1) [42,43]. At the front and rear wheels, the slip angles, αf and αr, are defined as (2) under the assumption that tan−1(θ) ≈ θ for small θ. The lateral tire forces, Fyf and Fyr, at the front and rear wheels are calculated as (3) under the assumption that the lateral tire force is linear with respect to the slip angle. By combining (1), (2), and (3), the linear state-space equations for the bicycle model are derived as (4).
m v x β ˙ + γ = F y f + F y r I z γ ˙ = l f F y f l r F y r + Δ M z
α f = δ f tan 1 β + l f γ v x δ f β l f γ v x α r = δ r tan 1 β l r γ v x δ r β + l r γ v x
F y f = 2 C f α f , F y r = 2 C r α r
β ˙ = ς 1 m v x β + ς 2 m v x 2 1 γ + C f m v x δ f + C f m v x δ r γ ˙ = ς 2 I z β + ς 3 I z v x γ + l f C f I z δ f l r C r I z δ r + Δ M z I z ς 1 = C f C r , ς 2 = C f l f + C r l r , ς 3 = l f 2 C f l r 2 C r
In the literature of PTC, the lateral offset and heading errors, ey and eφ, have been defined at point P in Figure 2. In this paper, a lookahead or preview function is introduced to improve the path-tracking performance when calculating those errors [12,21,22,24,34,38,40]. The lookahead or preview distance Lp is calculated using (5). In (5), tp is the preview interval or preview time, which acts as a velocity gain. In Figure 2, point C is located at the center of gravity (C.G.) of a vehicle. With Lp, the lookahead point Q and the point R on the target path are obtained. At point R, the time-derivative of ey and eφ is calculated as (6) under the assumption that eφ is smaller than 10° and its consequence is sineφeφ.
L p = t p v x
e ˙ y = v x sin e φ v x β L p γ v x e φ v x β L p γ e ˙ φ = φ ˙ d φ ˙ = v x χ γ
With the state variables of ey, eφ, β, and γ, the state vector, x, the disturbance, w, and the control input, u, are defined as (7). From (4), (6), and (7), the state-space equation for PTC and VSC is obtained as (8) [21,22,24]. In (8), the system, disturbance, and input matrices, A, B1, and B2, are defined in (9).
x = e y e φ β γ T , w = χ , u = δ f δ r Δ M z T
x ˙ = A x + B 1 w + B 2 u
A = 0 v x v x L p 0 0 0 1 0 0 ς 1 m v x ς 2 m v x 2 1 0 0 ς 2 I z ς 3 I z v x , B 1 = 0 v x 0 0 , B 2 = 0 0 0 0 0 0 C f m v x C r m v x 0 l f C f I z l r C r I z 1 I z
The control input u defined in (7) is composed of δf, δr, and ∆Mz. From those elements of the control input u, four control inputs, u1, u2, u3, and u4, are set as (10), which correspond to the input configurations, IC#1, IC#2, IC#3, and IC#4, respectively [13]. According to four control inputs, u1, u2, u3, and u4, the input matrices, B21, B22, B23, and B24, are defined in (10) from (9) [13]. In (10), B2(k) represents the k-th column of the matrix B2 in (9). As shown in (10), δf is available for all ICs. δr is available for IC#2, IC#3, and IC#4. For IC#3 and IC#4 with ∆Mz, 4WIB and 4WID are available when generating ∆Mz. ∆Mz in IC#3 and IC#4 should be converted into δr and the braking and traction torques (TBi and TDi) at each wheel using a control allocation procedure. For IC#4, RWS, i.e., δr, is not used when converting ∆Mz into tire forces because it is already included in the control input, u4. In this paper, it is assumed that TBi and TDi are generated via 4WIB and 4WID, respectively. In this paper, IC#1 and IC#2 are used as a baseline of PTC#1. Moreover, IC#1 is used as a steering angle generator for PTC#2.
u 1 = δ f B 21 = B 2 1 for I C # 1 u 2 = δ f δ r T B 22 = B 2 1 B 2 2 for I C # 2 u 3 = δ f Δ M z T B 23 = B 2 1 B 2 3 for I C # 3 u 4 = δ f δ r Δ M z T B 24 = B 2 1 B 2 2 B 2 3 for I C # 4

2.2. Design of LQR for PTC#1 and PTC#2

LQ objective functions, J1, J2, J3, and J4, for input configurations, IC#1, IC#2, IC#3, and IC#4, of (10) are provided in (11), respectively. As shown in (11), the first two terms of J0 including ey and eφ represent the path-tracking objective, and the last two terms of J0 including β and γ represent the vehicle stability objective, which can be divided into lateral stability and yaw rate tracking. In this manner, PTC and VSC are integrated into a single control framework. J1, J2, J3, and J4 can be converted into the vector matrix form of (12). In (12), the matrices Q and Ri are provided as (13). In (11) and (13), the weight ρi is set using Bryson’s rule, (14), where ξi is the maximum allowable value of the corresponding term [44]. The actual tuning procedure for path tracking and vehicle stability is carried out by adjusting ξi. The control input ui of IC#i in LQR is calculated as (15), where Si is the solution of the Riccati equation for IC#i.
J 0 = 0 ρ 1 e y 2 + ρ 2 e φ 2 + ρ 3 β 2 + ρ 4 γ 2 d t J 1 = J 0 + 0 ρ 5 δ f 2 d t for I C # 1 J 2 = J 0 + 0 ρ 5 δ f 2 + ρ 6 δ r 2 d t for I C # 2 J 3 = J 0 + 0 ρ 5 δ f 2 + ρ 7 Δ M z 2 d t for I C # 3 J 4 = J 0 + 0 ρ 5 δ f 2 + ρ 6 δ r 2 + ρ 7 Δ M z 2 d t for I C # 4
J i = = 0 x u i T Q 0 0 R i x u i d t , i = 1 , , 4
Q = d i a g ρ 1 , ρ 2 , ρ 3 , ρ 4 , R 1 = ρ 5 for I C # 1 R 2 = d i a g ρ 5 , ρ 6 for I C # 2 R 3 = d i a g ρ 5 , ρ 7 for I C # 3 R 4 = d i a g ρ 5 , ρ 6 , ρ 7 for I C # 4
ρ i = 1 ξ i 2
u i = K i x = R i 1 B 2 i T S i x , i = 1 , , 4

2.3. Derivation of the Reference Yaw Rate for PTC#2 and PTC#3

There are two methods that have been used to derive the reference yaw rate, γref. The first method uses the FWS angle and the bicycle model provided in (4), and the second is to use a target path. In the first method, γref is calculated to γs as provided in (16) [12,18,19,31,32,33,35,36,37,42]. The front steering angle, δf, is obtained from the LQR with IC#1 of PTC#1. In (16), Kγ is the steady-state yaw rate gain. In the previous study, it was shown that Kγ can be tuned regardless of Cf, Cr, m, and vx for the purpose of enhancing path-tracking performance [11]. Generally, γs is bounded by (17) [43].
γ s = 2 C f C r l f + l r v x 2 C f C r l f + l r 2 + m v x 2 l r C r l f C f δ f = K γ δ f
γ s 0.85 μ g v x
In the second method, γref is obtained from a target path without a vehicle model [11,38,39,40]. The typical method for this purpose is to derive γp from Figure 1 as provided in (18) [38]. Another method is to use a polynomial interpolation, which is adopted in this paper [11,12,39]. Figure 1 shows the desired and target paths and three points, C, D, and Q. Let the coordinate values of the points, C, D, and Q, be (x0,y0), (x1,y1), and (x2,y2), respectively. In Figure 1, the desired path is marked as a dashed bold line, connecting points C and D. To derive γref for PTC, the curvature of the desired path should be calculated. Under the assumption that point D is close enough to the vehicle, the desired path is represented with a second-order polynomial, as provided in (19). Three coefficients, a1, a2, and a3, of the polynomial are determined using (20) from the positions of the points C and D and the slope at the point C. With the coefficients of the polynomial, the curvature χd at point C is calculated as in (21). Finally, γref for PTC is calculated into γd as in (22) [11,39]. In (22), the preview interval Kp plays an identical role to Kγ in (16). In this paper, Kp is set to 1 and γd is also limited by (17).
γ p = 2 t p tan 1 e y L p β
y ( x ) = a 1 x 2 + a 2 x + a 3
a 1 a 2 a 3 = x 0 2 x 0 1 x 1 2 x 1 1 2 x 0 1 0 1 y 0 y 1 y ( x 0 )
χ d = y x 0 1 + y x 0 2 3 / 2 = 2 a 1 1 + 2 a 1 x 0 + a 2 2 3 / 2
γ d = K p v x χ d

2.4. Design of Yaw Rate Tracking Controller for PTC#2 and PTC#3

To make γ follow the reference yaw rates, γs or γd, it is necessary to calculate ΔMz. The yaw rate error, γe, is defined as (23), where this is set to the error surface, s. To make the surface asymptotically converge to zero, the condition of (24) should be satisfied [11,12,18,19,39]. By combining (1), (23), and (24), ∆Mz is computed as (25). In this paper, Kc is set to 5.
s = γ e = γ γ r e f
s ˙ = K c s K c > 0
Δ M z = I z γ ˙ d l f F y f + l r F y r I z K c γ γ r e f

2.5. Control Allocation with WLS for PTC#1, PTC#2, and PTC#3

The control yaw moment, ΔMz, is obtained commonly from PTC#1, PTC#2, and PTC#3. Once ∆Mz is obtained, it should be converted into the longitudinal of wheels, which are generated using the actuators, i.e., 4WIB and 4WID, via a control allocation. For control allocation, a weighted least square (WLS)-based method is adopted in this paper [20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37].
Figure 3 shows tire forces at each wheel and ∆Mz, when ∆Mz is positive [13,18,19,45]. In Figure 3, ∆Fyf and ∆Fyr are the lateral tire forces of the front and rear wheels, which are generated via FWS and RWS, respectively. ∆Fx1, ∆Fx2, ∆Fx3, and ∆Fx4 are the longitudinal tire forces generated via 4WID and 4WIB. If ∆Fi is positive, it is generated as the traction torque TDi via 4WID. Otherwise, it is generated as the braking torque TBi via 4WIB. These six tire forces should be determined to generate ∆Mz. The WLS-based method is used for this purpose.
From the geometric relationship shown in Figure 3, the equilibrium condition among ∆Mz and tire forces is obtained as (26) [11,12,18,19,39]. In (26), the elements of the vector p are provided in (27).
p 1 p 2 p 3 p 4 p 5 p 6 p Δ F y f Δ F y r Δ F x 1 Δ F x 2 Δ F x 3 Δ F x 4 q = Δ M z
p 1 = 2 l f cos δ f , p 2 = 2 l r cos δ r , p 3 = l f sin δ f + t f cos δ f , p 4 = l f sin δ f t f cos δ f , p 5 = l r sin δ r + t r cos δ r , p 6 = l r sin δ r t r cos δ r
The objective function of WLS is defined with the quadratic form as in (28). In (28), μFzi represents the radius of the friction circle at the wheel i. Assuming that the road surface is flat, it can be easily estimated using acceleration signals measured by sensors [46]. In (28), κ is the vector of virtual weights, κi, which is used to select the combination of actuators [39,47,48]. The equality constraint (26) is converted into the quadratic form (29). In the previous study, (26) must be satisfied in order to generate ∆Mz. The objective function (28) and the constraint (29) are combined into a single objective one (30), with the Lagrange multiplier, ζ. In (30), ζ should be set to 1 or higher. In this paper, ζ is set to 10. Otherwise, the constraint (26) is not satisfied. The optimum of (30) is algebraically obtained as (31) by differentiating (30) with respect to q [13].
J E = κ 1 Δ F y f 2 + κ 3 Δ F x 1 2 μ F z 1 2 + κ 1 Δ F y f 2 + κ 4 Δ F x 2 2 μ F z 2 2 + κ 2 Δ F y f 2 + κ 5 Δ F x 3 2 μ F z 3 2 + κ 2 Δ F y f 2 + κ 6 Δ F x 4 2 μ F z 4 2 = p T W p W = diag 1 μ F z 1 2 + 1 μ F z 2 2 , 1 μ F z 3 2 + 1 μ F z 4 2 , 1 μ F z 1 2 , 1 μ F z 2 2 , 1 μ F z 3 2 , 1 μ F z 4 2 κ κ = diag κ 1 κ 2 κ 3 κ 4 κ 5 κ 6
J E C = p q Δ M z T p q Δ M z
J C A = J E + ζ J E C = q T W q + ζ p q Δ M z T p q Δ M z
q o p t = ζ W + η p T p 1 p T Δ M z
When applying WLS for control allocation with ∆Mz, arbitrary actuator combinations composed of FWS, RWS, 4WS, 4WIB, and 4WID can be set according to a particular input configuration. For example, if u4 of IC#4 in (10) is selected as a control input, only 4WID and 4WIB are available to be used in generating ∆Mz because δf and δr of FWS, RWS, and 4WS are already determined using LQR. To reflect this fact, the virtual weights, κi, can be set according to a particular combination of multiple actuators, FWS, RWS, 4WS, 4WIB, and 4WID [39,47,48].
The vectors of virtual weights corresponding to steering actuators, i.e., FWS, RWS, and 4WS, are provided in (32). Initially, all the virtual weights are set to 1. In (32), ε is 10−4 and the terms denoted by have no effect on the steering actuators. The first and second terms of κ in (32) correspond to ∆Fyf and ∆Fyr, respectively, which are converted into the front and rear steering angles. If the first term of κ is set to ε, then only ∆Fyf is generated via WLS. This is shown in the first row of (32).
FWS : κ = diag ε 1 RWS : κ = diag 1 ε 4 WS : κ = diag ε ε
The vectors of virtual weights representing 4WIB and 4WID are provided in (33). Figure 3 shows the fact that only ∆Fx2 and ∆Fx4 should be generated if ∆Mz is negative and 4WIB is available. This is represented by the second row of the first term, 4WIB in (33). If 4WIB and 4WID are available for generating ∆Mz, no constraints are imposed on the longitudinal forces, which is represented by the last row of (33), i.e., 4WID + 4WIB. The virtual weights provided in (32) and (33) can be set according to combinations of actuators available for the purpose of generating ∆Mz.
4 WID : κ = diag 1 ε 1 ε if Δ M z > 0 κ = diag ε 1 ε 1 if Δ M z < 0 4 WIB : κ = diag ε 1 ε 1 if Δ M z > 0 κ = diag 1 ε 1 ε if Δ M z < 0 4 WID + 4 WIB : κ = diag ε ε ε ε
In this paper, TBi and TDi at the wheel i are calculated from ∆Fx1, ∆Fx2, ∆Fx3, and ∆Fx4, determined via the WLS-based method, i.e., qopt in (31). Following the sign of ∆Fxi, TBi, and TDi at the wheel i are calculated as (34), respectively [13].
T B i if Δ F x i < 0 T D i if Δ F x i > 0 = h r w i Δ F x i υ i , ω i , i = 1 , , 4
δf and δr are calculated from ∆Fyf and ∆Fyr obtained in qopt and the definitions of (2) and (3). The definition of the lateral tire force (3) is rewritten as (35). In (35), σ is the parameter representing the magnitude of the cornering stiffness, Ci [34]. In this paper, σ is set to 6. For 4WS, δf and δr are calculated as (36) by combining (2) with (35) [13,49]. In these formulations, β is needed. However, it is hard to measure. For this reason, an extended Kalman Filter is adopted to estimate β [50] in this paper.
ϕ i = Δ F y i σ C i , i = f , r
δ f = α f + β + l f γ v x = ϕ f + β + l f γ v x δ r = α r + β l r γ v x = ϕ r + β l r γ v x

3. Performance Measures for Path Tracking and Vehicle Stability

In the literature of PTC, ey and eφ have been used to represent the path-tracking performance. Different from the literature on PTC, a double lane change maneuver for collision avoidance is adopted as a target path in this paper [11,12,13,20,26,27,37]. This path has been used to measure the reachability and agility of a path-tracking controller [12,13]. Figure 4 shows the target path and vehicle trajectory. From the points on the target path and the vehicle trajectory in Figure 4, five measures standing for path-tracking performance are defined as (37): the peak’s center offset, ∆X, the peak’s lateral offset, ∆Y, the percentage over-shoot, OS%, the response delay, ∆DX, and the settling delay, ∆SX. These measures are illustrated in Figure 4. In (37), the subscripts X and Y represent the x- and y-positions of the point , respectively. Basically, the smaller the absolute magnitudes of these measures, the better the path-tracking performance.
Δ X = D X A X = D X 73.20 Δ Y = D Y A Y = D Y 3.53 Δ D X = E X B X = E X 91.50 Δ S X = G X C X = G X 190.00 O S % = F Y 1.65 1.65 + 3.53 × 100
In Figure 4, ∆X and ∆Y represent the agility and the reachability, respectively. In this paper, if ∆Y is larger than −0.05 m, it is regarded that the path-tracking performance of a controller is satisfactory. OS% stands for lateral damping, which represents agility. In this paper, it is satisfactory if OS% is less than 16%, which is equal to the overshoot of 0.85 m. ∆DX stands for the response delay, which stands for the agility of the longitudinal motion. ∆SX can be regarded as the settling time, which stands for the convergence speed of the vehicle’s lateral motion to a target value of y-position. In this paper, it is satisfactory if ∆SX is less than 16 m. Generally, if OS% is near 1% under the condition that ∆Y is larger than −0.05 m, the best performance is achieved. Detailed explanations of these measures can be found in [12,13].
As mentioned earlier, VSC stands for yaw rate tracking and lateral stability. Generally, yaw rate tracking performance is measured with the yaw rate error γe and lateral stability performance is measured with β. It has been regarded that lateral stability performance is satisfactory if β is less than 3° [51]. Let us denote the maximum absolute value of β to be MASSA. The side-slip angular rate, combined with β, has been also used as a measure of lateral stability. Let us denote the maximum absolute value of the side-slip angular rate to be MASSAR [12].

4. Simulation and Discussion

4.1. Simulation Environment

A simulation was conducted to compare three controllers, PTC#1, PTC#2, and PTC#3, on low-friction roads in terms of path-tracking performance. These controllers were implemented using MATLAB/Simulink (2019a) and connected with the vehicle simulation software CarSim 8 [52]. The configuration of MATLAB/Simulink and CarSim is provided in Figure 1. In Figure 1, the autonomous electrical vehicles with FWS/4WS, 4WID, and 4WIB are replaced by CarSim, and the other parts are implemented with MATLAB/Simulink. A target path in Figure 4 was selected as a test scenario. In the simulation, the built-in F-segment sedan model provided in the CarSim software was selected [52]. From the model, the parameters of the bicycle model are provided in Table 1. In Table 1, the cornering stiffness, Cf and Cr, are obtained from simulation in CarSim by comparing β and γ of CarSim and the state-space model under a step steering maneuver. The steering actuators of FWS and 4WS were modeled as the first-order system with a time constant of 0.05. 4WIB and 4WID were modeled as the first-order system with a time constant of 0.1. The initial vx and μ were set to 60 km/h and 0.4, respectively. The speed controller provided in CarSim was applied upon driving to maintain a constant speed, which means that Kγ is maintained as constant as possible.
In real vehicles, the magnitude of ∆Mz is limited by the maximum tire forces and road surface friction. To reflect this fact, it is necessary to set a bound on ∆Mz. For this purpose, ∆Mz is limited to a certain value after calculating it using LQR. For example, IC#4 of PTC#1 can use the smallest ∆Mz among the input configurations because only the longitudinal tire forces are available for control allocation. A detailed description of how to select the limit of ∆Mz can be found in Reference [13].
The tuning was carried out by trial and error such that ∆Y is larger than −0.05 m and OS% is near 0.9, which gives the smallest ∆DX and ∆SX. The most important tuning parameters are the velocity gain tp, ζ1 and ζ2 in (14). If RWS is adopted, the weight on RWS, i.e., ζ6 in (14), should be set to a small value such that β is smaller than 2.5°.

4.2. Simulation with PTC#1

The first simulation was carried out with PTC#1. For IC#1, IC#2, IC#3, and IC#4, several actuator combinations composed of RWS, 4WS, 4WID, and 4WIB were used for control allocation. The tuning parameters of PTC#1 are tp and ζi as provided in (5) and (14), respectively. The simulation results of PTC#1 are summarized in Table 2. The velocity gains and weights of the LQ objective function are provided in Table 3. As shown in Table 2, there are few differences among the input configurations and actuator combinations in terms of path-tracking performance. A notable feature is that the use of RWS increases MASSA. This is caused by the fact that RWS steers the rear wheels along the identical direction to the front wheels. These results are identical to those of [13].
For IC#3, ∆Mz is added to IC#1, i.e., FWS. As shown in Table 2, the ∆X of IC#3 is reduced by adding ∆Mz to IC#1. Especially, RWS and 4WID have a large effect on the reduction in ∆X. Among them, 4WID is preferred to RWS because 4WID does not increase MASSA. On the contrary, IC#4 adds ∆Mz to IC#2, i.e., 4WS. As shown in Table 2, there are few effects of adding ∆Mz to IC#2 on path-tracking performance although ∆X is slightly decreased by adding ∆Mz. In Table 2, the actuator combinations with RWS in IC#3 are equivalent to those of IC#4. As shown in Table 2, the former is better than the latter in terms of path-tracking performance. This fact means that the control allocation with multiple actuators for 4WS is not effective. Instead, it is desirable to add ∆Mz to IC#1 and use control allocation with multiple actuators. Besides the increase in MASSA and the reduction in ∆X, there are few effects of adding ∆Mz to IC#1 or IC#2. Figure 5, Figure 6 and Figure 7 show the trajectories and side-slip angles of every input configuration and actuator combination provided in Table 2. As shown in these figures, there are little differences among trajectories. On the contrary, only the side-slip angles are different from one another. From these results, it can be concluded that it is desirable to not use the multiple actuators needed to generate ∆Mz for path tracking.

4.3. Simulation with PTC#2 and PTC#3

The second simulation was carried out with PTC#2 and PTC#3. In PTC#2, the front steering angle needed to calculate the reference yaw rate is generated using the IC#1 of PTC#1. This is identical to PTC#3 except the reference yaw rate is calculated from a target path. With the reference yaw rate, ∆Mz is calculated using SMC. For this reason, RWS, 4WID, and 4WIB were used to generate ∆Mz in PTC#2 and PTC#3. Because PTC#2 and PTC#2 use the IC#1 of PTC#1 to generate δf, the tuning parameters of PTC#2 and PTC#3 are identical to those of PTC#1.
Table 4 and Table 5 summarize the simulation results of PTC#2 and PTC#3, respectively. Table 6 shows the best and worst results for PTC#1, PTC#2, and PTC#3. In Table 5, IC#1 and IC#2 of PTC#1 are selected as a baseline, and the gray rows indicate the worst results. Table 7 and Table 8 show the velocity gains and weights of LQ objective functions for PTC#2 and PTC#3, which correspond to Table 4 and Table 5, respectively. Figure 8, Figure 9, Figure 10 and Figure 11 show the trajectories and the side-slip angles of each actuator combination of PTC#2 and PTC#3, respectively. As shown in Table 4, the use of RWS increases MASSA, which is identical to PTC#1. Instead, ∆X is significantly reduced by RWS. As shown in Table 5, this tendency can also be identified via PTC#3. Without RWS, the 4WIB of PTC#2 provided the best result among PTC#1, PTC#2, and PTC#3. This fact can be identified in Figure 8. With RWS, the 4WID of PTC#2 provided a good result although MASSA was increased via RWS. By comparing the results among Table 2, Table 4 and Table 6, it can be identified that PTC#2 is superior to PTC#1 in terms of path-tracking performance. By comparing PTC#1 and PTC#3 as shown in Table 1, Table 5 and Table 6, PTC#1 and PTC#3 have no superiority over each other in terms of path-tracking performance.
Although PTC#2 provided better results than PTC#1, there are few differences among PTC#1, PTC#2, and PTC#3, as shown in the trajectories of Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11. The side-slip angles are less than 3° for all controllers, which is acceptable in terms of lateral stability. PTC#1 and PTC#3 require a two-level controller. Especially, IC#1 and IC#2 of PTC#1 do not require any control allocation schemes. On the contrary, PTC#2 requires a three-level controller including a yaw rate tracking controller and control allocation. As shown in Section 2, the yaw rate tracking controller and control allocation require several variables that are difficult to measure. For this reason, IC#1 and IC#2 of PTC#1 are the best choices for designing a path-tracking controller because they are the simplest among the controllers. Moreover, IC#1 of PTC#1 is preferred to IC#2 because it gives a smaller side-slip angle than IC#2.

5. Conclusions

The aim of this study is to compare three types of coordinated controllers, i.e., PTC#1, PTC#2, and PTC#3 on low-friction roads. With the state-space model derived from a 2-DOF bicycle model and a target path, the three controllers were designed. To calculate the steering angles and yaw moment control, LQR was designed. For yaw rate tracking, SMC was designed. As an actuator for path tracking and yaw moment generation, FWS, RWS, 4WS, 4WID, and 4WIB were adopted. For yaw moment distribution or control allocation with the multiple actuators, the WLS-based method was applied. To check the path-tracking performance of each controller, a simulation was carried out using the vehicle simulation package CarSim. For simulation, double lane change for collision avoidance was adopted as a test scenario. From the simulation results, it was shown that PTC#2 is superior to PTC#1 and PTC#3, and that there are little differences among PTC#1 and PTC#3. For this reason, PTC#3 is not recommended as a path-tracking controller. However, there are small differences between PTC#1 and PTC#2 in terms of path-tracking performance except for the peak’s lateral offset and the settling delay. PTC#1 requires a two-level controller and PTC#2 uses a three-level control structure. Moreover, IC#1 and IC#2 of PTC#1 do not require yaw rate tracking and control allocation. As provided in Section 2, the yaw rate tracking controller and control allocation require several variables that are difficult to measure. Because the simplest has been shown to be the best, it is recommended that PTC#1 is the best controller for coordinated control of path tracking and vehicle stability control on low-friction roads. This conclusion was also drawn in the previous research [12]. As pointed out in the previous research, there are no significant differences among actuator combinations [13]. In other words, FWS or 4WS itself is quite enough for coordinated control of path tracking and vehicle stability control. Another conclusion from the simulation results is that it is important to tune the controller gain or select a controller type that is effective on both high- and low-friction roads.

Author Contributions

M.P. conceptualized the main idea and designed this study. S.Y. participated in formulating the idea, as well as validating the proposed method and results. S.Y. implemented the methodology and drafted the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Ministry of Trade, Industry and Energy (MOTIE, Korea). (Project Name: Development of Fail Operation technology in Lv.4 autonomous driving systems/Project Number: 20018055). Also, this work is supported by the Korea Agency for Infrastructure Technology Advancement (KAIA) grant funded by the Ministry of Land, Infrastructure and Transport (Grant RS-2021-KA160501).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

4WSFour-wheel steering
4WIBFour-wheel independent braking
4WIDFour-wheel independent drive
FWSFront-wheel steering
MASSAMaximum absolute side-slip angle
RWSRear-wheel steering
WLSWeighted least square
Cf, CrCornering stiffness of front and rear tires (N/rad)
ey, eφLateral offset error (m) and heading error (rad)
Fxi, FziLongitudinal and vertical tire forces of i-th wheel (N)
Fyf, FyrFront and rear lateral tire forces in the 2-DOF bicycle model (N)
h(ω)A function of rotor speed representing the capacity curve of an electric motor
IzYaw moment of inertial (kg·m2)
JiLQ objective function for the input configuration IC#i
KiGain matrix of LQR for input configuration IC#i
LpLookahead distance (m)
LWLCLower-level controller
lf, lrDistance from cog to front and rear axles (m)
mVehicle total mass (kg)
MDLCMiddle-level controller
OS%Percentage of overshoot in the lower lane of the target path
pVector used for the equality constraint in the WLS-based method
qVector of tire forces as a solution of the WLS-based method
rwiRadius of i-th wheel (m)
TBi, TDiBraking and traction torques applied at i-th wheel (N·m)
UPLCUpper-level controller
tf, trHalf of track widths of front and rear axles (m)
tpPreview time for lookahead distance
vx, vyLongitudinal and lateral velocities of cog of a vehicle (m/s)
WWeighting matrix of WLS-based method
X(∗), Y(∗)X- and y-positions of the point ∗ on the target path and vehicle trajectory
Yref(X)y-position of the target path with respect to X
αf, αrTire slip angles of front and rear wheels (rad)
βSide-slip angle of C.G. of a vehicle (rad) = tan−1(vy/vx) ≈ (vy/vx)
δf, δrFront and rear steering angles (rad)
ε10−4 used as a virtual weight in WLS-based method
FxiControl longitudinal and lateral forces generated by an actuator (N)
MzControl yaw moment as a control input in LQR (N·m)
X, ∆YDifferences between x- and y-positions at the peak points of the target path
DX, ∆SXResponse and settling delays of vehicle trajectory with respect to target path
γ, γrefMeasured and referenced yaw rates (rad/s)
γs, γp, γdReference yaw rates obtained from steering angle and target path (rad/s)
ζTuning parameter on relaxation term of equality constraint
χCurvature at a particular point on a target path.
κVirtual weight on the longitudinal and lateral tire forces
κVector of virtual weights
ξiThe maximum allowable value of i-th term in the LQ objective function
ζVector of the maximum allowable values
φHeading angle of a vehicle
φdDesired heading angle obtained with lookahead
Ψref(X)Heading angle of the target path with respect to X
μTire-road friction coefficient
ωiRotational speed of i-th wheel (rad/s)
ρiWeight on i-th term in LQ objective function
τTime constant of the first-order system used for reference yaw rate
ʋiRatio of reduction gear of i-th wheel

References

  1. Montanaro, U.; Dixit, S.; Fallaha, S.; Dianatib, M.; Stevensc, A.; Oxtobyd, D.; Mouzakitisd, A. Towards connected autonomous driving: Review of use-cases. Veh. Syst. Dyn. 2019, 57, 779–814. [Google Scholar] [CrossRef]
  2. Yurtsever, E.; Lambert, J.; Carballo, A.; Takeda, K. A survey of autonomous driving: Common practices and emerging technologies. IEEE Access 2020, 8, 58443–58469. [Google Scholar] [CrossRef]
  3. Omeiza, D.; Webb, H.; Jirotka, M.; Kunze, M. Explanations in autonomous driving: A survey. IEEE Trans. Intell. Transp. Syst. 2022, 23, 10142–10162. [Google Scholar] [CrossRef]
  4. Sorniotti, A.; Barber, P.; De Pinto, S. Path tracking for automated driving: A tutorial on control system formulations and ongoing research. In Automated Driving; Watzenig, D., Horn, M., Eds.; Springer: Cham, Switzerland, 2017. [Google Scholar]
  5. Amer, N.H.; Hudha, H.Z.K.; Kadir, Z.A. 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. 2017, 86, 225–254. [Google Scholar] [CrossRef]
  6. Paden, B.; Cap, M.; Yong, S.Z.; Yershov, D.; Frazzoli, E. A survey of motion planning and control techniques for self-driving urban vehicles. IEEE Trans. Intell. Veh. 2016, 1, 33–55. [Google Scholar] [CrossRef]
  7. Bai, G.; Meng, Y.; Liu, L.; Luo, W.; Gu, Q.; Liu, L. Review and comparison of path tracking based on model predictive control. Electronics 2019, 8, 10. [Google Scholar] [CrossRef]
  8. Yao, Q.; Tian, Y.; Wang, Q.; Wang, S. Control strategies on path tracking for autonomous vehicle: State of the art and future challenges. IEEE Access 2020, 8, 161211–161222. [Google Scholar] [CrossRef]
  9. Rokonuzzaman, M.; Mohajer, N.; Nahavandi, S.; Mohamed, S. Review and performance evaluation of path tracking controllers of autonomous vehicles. IET Intell. Transp. Syst. 2021, 15, 646–670. [Google Scholar] [CrossRef]
  10. Stano, P.; Montanaro, U.; Tavernini, D.; Tufo, M.; Fiengo, G.; Novella, L.; Sorniotti, A. Model predictive path tracking control for automated road vehicles: A review. Annu. Rev. Control 2022, 55, 194–236. [Google Scholar] [CrossRef]
  11. Jeong, Y.; Yim, S. Integrated path tracking and lateral stability control with four-wheel independent steering for autonomous electric vehicles on low friction roads. Machines 2022, 10, 650. [Google Scholar] [CrossRef]
  12. Lee, J.; Yim, S. Comparative study of path tracking controllers on low friction roads for autonomous vehicles. Machines 2023, 11, 403. [Google Scholar] [CrossRef]
  13. Park, M.; Yim, S. Comparative study on effects of input configurations of linear quadratic controller on path tracking performance under low friction condition. Actuators 2023, 12, 151. [Google Scholar] [CrossRef]
  14. Ahangarnejad, A.H.; Radmehr, A.; Ahmadian, M. A review of vehicle active safety control methods: From antilock brakes to semiautonomy. J. Vib. Control 2021, 27, 1683–1712. [Google Scholar] [CrossRef]
  15. Zhang, L.; Zhang, Z.; Wang, Z.; Deng, J.; Dorrell, D.G. Chassis coordinated control for full x-by-wire vehicles—A review. Chin. J. Mech. Eng. 2021, 34, 42. [Google Scholar] [CrossRef]
  16. Peng, H.; Chen, X. Active safety control of x-by-wire electric vehicles: A survey. SAE Int. J. Veh. Dyn. Stab. NVH 2022, 6, 115–133. [Google Scholar] [CrossRef]
  17. Indu, K.; Kumar, M.A. Electric vehicle control and driving safety systems: A review. IETE J. Res. 2023, 6, 482–498. [Google Scholar] [CrossRef]
  18. Yim, S. Comparison among active front, front independent, 4-wheel and 4-wheel independent steering systems for vehicle stability control. Electronics 2020, 9, 798. [Google Scholar] [CrossRef]
  19. Nah, J.; Yim, S. Vehicle stability control with four-wheel independent braking, drive and steering on in-wheel motor-driven electric vehicles. Electronics 2020, 9, 1934. [Google Scholar] [CrossRef]
  20. Hang, P.; Chen, X.; Luo, F. Path-Tracking Controller Design for a 4WIS and 4WID Electric Vehicle with Steer-by-Wire System (No. 2017-01-1954); SAE Technical Paper: Warrendale, PA, USA, 2017. [Google Scholar]
  21. Guo, J.; Luo, Y.; Li, K. An adaptive hierarchical trajectory following control approach of autonomous four-wheel independent drive electric vehicles. IEEE Trans. Intell. Transp. Syst. 2018, 19, 2482–2492. [Google Scholar] [CrossRef]
  22. Guo, J.; Luo, Y.; Li, K.; Dai, Y. Coordinated path-following and direct yaw-moment control of autonomous electric vehicles with sideslip angle estimation. Mech. Syst. Signal Process. 2018, 105, 183–199. [Google Scholar] [CrossRef]
  23. Hang, P.; Chen, X.; Luo, F. LPV/H controller design for path tracking of autonomous ground vehicles through four-wheel steering and direct yaw-moment control. Int. J. Automot. Technol. 2019, 20, 679–691. [Google Scholar] [CrossRef]
  24. Peng, H.; Wang, W.; An, Q.; Xiang, C.; Li, L. 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. 2020, 69, 6053–6066. [Google Scholar] [CrossRef]
  25. Wang, Y.; Shao, Q.; Zhou, J.; Zheng, H.; Chen, H. Longitudinal and lateral control of autonomous vehicles in multi-vehicle driving environments. IET Intell. Transp. Syst. 2020, 14, 924–935. [Google Scholar] [CrossRef]
  26. Wu, H.; Li, Z.; Si, Z. Trajectory tracking control for four-wheel independent drive intelligent vehicle based on model predictive control and sliding mode control. Adv. Mech. Eng. 2021, 13, 16878140211045142. [Google Scholar] [CrossRef]
  27. Xiang, C.; Peng, H.; Wang, W.; Li, L.; An, Q.; Cheng, S. Path tracking coordinated control strategy for autonomous four in-wheel-motor independent-drive vehicles with consideration of lateral stability. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2021, 235, 1023–1036. [Google Scholar] [CrossRef]
  28. Barari, A.; Afshari, S.S.; Liang, X. Coordinated control for path-following of an autonomous four in-wheel motor drive electric vehicle. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2022, 236, 6335–6346. [Google Scholar] [CrossRef]
  29. Wang, W.; Zhang, Y.; Yang, C.; Qie, T.; Ma, M. Adaptive model predictive control-based path following control for four-wheel independent drive automated vehicles. IEEE Trans. Intell. Transp. Syst. 2022, 23, 14399–14412. [Google Scholar] [CrossRef]
  30. Jin, L.; Gao, L.; Jiang, Y.; Chen, M.; Zheng, Y.; Li, K. Research on the control and coordination of four-wheel independent driving/steering electric vehicle. Adv. Mech. Eng. 2017, 9, 1687814017698877. [Google Scholar] [CrossRef]
  31. Ren, Y.; Zheng, L.; Khajepour, A. Integrated model predictive and torque vectoring control for path tracking of 4-wheeldriven autonomous vehicles. IET Intell. Transp. Syst. 2019, 13, 98–107. [Google Scholar] [CrossRef]
  32. Guo, L.; Ge, P.; Yue, M.; Li, J. Trajectory tracking algorithm in a hierarchical strategy for electric vehicle driven by four independent in-wheel motors. J. Chin. Inst. Eng. 2020, 43, 807–818. [Google Scholar] [CrossRef]
  33. Ahn, T.; Lee, Y.; Park, K. Design of integrated autonomous driving control system that incorporates chassis controllers for improving path tracking performance and vehicle stability. Electronics 2021, 10, 144. [Google Scholar] [CrossRef]
  34. Liang, Y.; Li, Y.; Zheng, L.; Yu, Y.; Ren, Y. Yaw rate tracking-based path-following control for four-wheel independent driving and four-wheel independent steering autonomous vehicles considering the coordination with dynamics stability. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2021, 235, 260–272. [Google Scholar] [CrossRef]
  35. Xie, J.; Xu, X.; Wang, F.; Tang, Z.; Chen, L. Coordinated control based path following of distributed drive autonomous electric vehicles with yaw-moment control. Control Eng. Pract. 2021, 106, 104659. [Google Scholar] [CrossRef]
  36. Tong, Y.; Jing, H.; Kuang, B.; Wang, G.; Liu, F.; Yang, Z. Trajectory tracking control for four-wheel independently driven electric vehicle based on model predictive control and sliding model control. In Proceedings of the 2021 5th CAA International Conference on Vehicular Control and Intelligence (CVCI), Tianjin, China, 29–31 October 2021. [Google Scholar]
  37. Du, Q.; Zhu, C.; Li, Q.; Tian, B.; Li, L. Optimal path tracking control for intelligent four-wheel steering vehicles based on MPC and state estimation. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2022, 236, 1964–1976. [Google Scholar] [CrossRef]
  38. Sun, X.; Wang, Y.; Hu, W.; Cai, Y.; Huang, C.; Chen, L. Path tracking control strategy for the intelligent vehicle considering tire nonlinear cornering characteristics in the PWA form. J. Frankl. Inst. 2022, 359, 2487–2513. [Google Scholar] [CrossRef]
  39. Jeong, Y.; Yim, S. Path tracking control with four-wheel independent steering, driving and braking systems for autonomous electric vehicles. IEEE Access 2022, 10, 74733–74746. [Google Scholar] [CrossRef]
  40. Hu, H.; Bei, S.; Zhao, Q.; Han, X.; Zhou, D.; Zhou, X.; Li, B. Research on trajectory tracking of sliding mode control based on adaptive preview time. Actuators 2022, 11, 34. [Google Scholar] [CrossRef]
  41. Tarhini, F.; Talj, R.; Doumiati, M. Adaptive look-ahead distance based on an intelligent fuzzy decision for an autonomous vehicle. In Proceedings of the 2023 IEEE Intelligent Vehicles Symposium (IV), Anchorage, AK, USA, 4–7 June 2023. [Google Scholar]
  42. Wong, H.Y. Theory of Ground Vehicles, 3rd ed.; John Wiley and Sons, Inc.: New York, NY, USA, 2001. [Google Scholar]
  43. Rajamani, R. Vehicle Dynamics and Control; Springer: New York, NY, USA, 2006. [Google Scholar]
  44. Bryson, A.E.; Ho, Y.C. Applied Optimal Control; Hemisphere: New York, NY, USA, 1975. [Google Scholar]
  45. Wang, J.; Wang, R.; Jing, H.; Chen, N. Coordinated active steering and four-wheel independently driving/braking control with control allocation. Asian J. Control 2016, 18, 98–111. [Google Scholar] [CrossRef]
  46. Rezaeian, A.; Zarringhalam, R.; Fallah, S.; Melek, W.; Khajepour, A.; Chen, S.-K.; Moshchuck, N.; Litkouhi, B. Novel tire force estimation strategy for real-time implementation on vehicle applications. IEEE Trans. Veh. Technol. 2015, 64, 2231–2241. [Google Scholar] [CrossRef]
  47. Yim, S.; Choi, J.; Yi, K. Coordinated control of hybrid 4WD vehicles for enhanced maneuverability and lateral stability. IEEE Trans. Veh. Technol. 2012, 61, 1946–1950. [Google Scholar] [CrossRef]
  48. Yim, S. Coordinated control with electronic stability control and active steering devices. J. Mech. Sci. Technol. 2015, 29, 5409–5416. [Google Scholar] [CrossRef]
  49. Yakub, F.; Mori, Y. Comparative study of autonomous path-following vehicle control via model predictive control and linear quadratic control. Proc. Inst. Mech. Eng. Part D J Automob. Eng. 2015, 229, 1695–1714. [Google Scholar] [CrossRef]
  50. Kim, H.H.; Ryu, J. Sideslip angle estimation considering short-duration longitudinal velocity variation. Int. J. Automot. Technol. 2011, 12, 545–553. [Google Scholar] [CrossRef]
  51. National Highway Traffic Safety Administration. FMVSS No. 126, Electronic Stability Control Systems: NHTSA Final Regulatory Impact Analysis; National Highway Traffic Safety Administration: Washington, DC, USA, 2007. [Google Scholar]
  52. Mechanical Simulation Corporation. VS Browser: Reference Manual, The Graphical User Interfaces of BikeSim, CarSim, and TruckSim; Mechanical Simulation Corporation: Ann Arbor, MI, USA, 2009. [Google Scholar]
Figure 1. Structures of PTC#1, PTC#2, and PTC#3.
Figure 1. Structures of PTC#1, PTC#2, and PTC#3.
Actuators 12 00398 g001
Figure 2. Coordinates and variables of the 2-DOF bicycle model and the target path.
Figure 2. Coordinates and variables of the 2-DOF bicycle model and the target path.
Actuators 12 00398 g002
Figure 3. Coordinate system on vehicle geometry with tire forces and yaw moment control.
Figure 3. Coordinate system on vehicle geometry with tire forces and yaw moment control.
Actuators 12 00398 g003
Figure 4. Measures for path-tracking performance.
Figure 4. Measures for path-tracking performance.
Actuators 12 00398 g004
Figure 5. Simulation results with IC#1 and IC#3 of PTC#1. (a) Trajectories. (b) Side-slip angles.
Figure 5. Simulation results with IC#1 and IC#3 of PTC#1. (a) Trajectories. (b) Side-slip angles.
Actuators 12 00398 g005aActuators 12 00398 g005b
Figure 6. Simulation results with IC#2 and IC#3 of PTC#1. (a) Trajectories. (b) Side-slip angles.
Figure 6. Simulation results with IC#2 and IC#3 of PTC#1. (a) Trajectories. (b) Side-slip angles.
Actuators 12 00398 g006
Figure 7. Simulation results with IC#2 and IC#4 of PTC#1. (a) Trajectories. (b) Side-slip angles.
Figure 7. Simulation results with IC#2 and IC#4 of PTC#1. (a) Trajectories. (b) Side-slip angles.
Actuators 12 00398 g007aActuators 12 00398 g007b
Figure 8. Simulation results with IC#1 and PTC#2. (a) Trajectories. (b) Side-slip angles.
Figure 8. Simulation results with IC#1 and PTC#2. (a) Trajectories. (b) Side-slip angles.
Actuators 12 00398 g008aActuators 12 00398 g008b
Figure 9. Simulation results with IC#2 and PTC#2. (a) Trajectories. (b) Side-slip angles.
Figure 9. Simulation results with IC#2 and PTC#2. (a) Trajectories. (b) Side-slip angles.
Actuators 12 00398 g009
Figure 10. Simulation results with IC#1 and PTC#3. (a) Trajectories. (b) Side-slip angles.
Figure 10. Simulation results with IC#1 and PTC#3. (a) Trajectories. (b) Side-slip angles.
Actuators 12 00398 g010aActuators 12 00398 g010b
Figure 11. Simulation results with IC#2 and PTC#3. (a) Trajectories. (b) Side-slip angles.
Figure 11. Simulation results with IC#2 and PTC#3. (a) Trajectories. (b) Side-slip angles.
Actuators 12 00398 g011
Table 1. Parameter of 2-DOF bicycle model obtained from F-segment sedan in CarSim.
Table 1. Parameter of 2-DOF bicycle model obtained from F-segment sedan in CarSim.
ParameterValueParameterValue
ms1823 kgIz6286 kg·m2
Cf42,000 N/radCr62,000 N/rad
lf1.27 mlr1.90 m
tf0.80 mtr0.80 m
Table 2. Summary of the simulation results for IC#1, IC#2, IC#3, and IC#4 of PTC#1.
Table 2. Summary of the simulation results for IC#1, IC#2, IC#3, and IC#4 of PTC#1.
Actuator
Combinations
X
(m)
Y
(m)
OS%DX
(m)
SX
(m)
MASSA
(deg)
MASSAR
(deg/s)
IC#1FWS1.570.0021.08.984.840.5813.13
IC#24WS1.610.0130.99.125.061.6714.73
IC#3FWS + 4WID1.19−0.0350.98.914.640.5913.34
+4WIB1.29−0.0241.08.804.630.590.59
+4WID + 4WIB1.30 −0.0151.08.714.590.5913.22
FWS + RWS1.25 −0.0030.98.744.581.3613.84
+4WID1.35 0.0001.08.874.701.3514.02
+4WIB1.590.0001.09.175.011.1513.11
+4WID + 4WIB1.60 0.0040.99.134.991.1313.30
IC#4+4WID1.45−0.0020.99.215.072.0915.30
+4WIB1.480.0041.09.105.042.0815.08
+4WID + 4WIB1.490.0100.99.045.021.9615.07
Table 3. Velocity gains and weights of LQ objective function for IC#1, IC#2, IC#3, and IC#4 of PTC#1.
Table 3. Velocity gains and weights of LQ objective function for IC#1, IC#2, IC#3, and IC#4 of PTC#1.
Actuator
Combinations
tpξ1ξ2ξ3ξ4ξ5ξ6ξ7
IC#1FWS0.600.56005.0000.3010.000.05--
IC#24WS0.600.55000.7000.3010.00.050.005-
IC#3FWS + 4WID1.190.05180.0050.100.100.03-50.0
+4WIB1.290.05320.0050.100.100.03-50.0
+4WID + 4WIB1.300.05720.0050.100.100.03-50.0
FWS + RWS1.250.06620.0050.100.100.03-50.0
+4WID1.350.06440.0050.100.100.03-50.0
+4WIB1.590.05850.0050.100.100.03-50.0
+4WID + 4WIB1.600.05980.0050.100.100.03-50.0
IC#4+4WID1.450.05660.0050.100.100.030.00150.0
+4WIB1.480.05690.0050.100.100.030.00150.0
+4WID + 4WIB1.490.06060.0050.100.100.030.00150.0
Table 4. Summary of the simulation results for PTC#2.
Table 4. Summary of the simulation results for PTC#2.
Actuator
Combinations
X
(m)
Y
(m)
OS%DX
(m)
SX
(m)
MASSA
(deg)
MASSAR
(deg/s)
+4WID1.50 −0.0120.98.704.530.5814.13
+4WIB0.52 −0.0360.96.662.400.6612.56
+4WID + 4WIB1.15 −0.0281.08.895.390.6413.39
+RWS0.390.0200.97.804.261.8412.41
+4WID0.42 0.0260.97.663.881.6712.39
+4WIB0.68 0.0230.98.024.241.5412.08
+4WID + 4WIB0.65 0.0150.78.294.881.4411.47
Table 5. Summary of the simulation results for PTC#3.
Table 5. Summary of the simulation results for PTC#3.
Actuator
Combinations
X
(m)
Y
(m)
OS%DX
(m)
SX
(m)
MASSA
(deg)
MASSAR
(deg/s)
+4WID1.75−0.0331.09.074.760.6813.59
+4WIB1.56−0.0291.07.704.020.7312.42
+4WID + 4WIB1.67−0.0241.09.345.700.6212.28
+RWS1.08−0.0341.08.415.591.7514.35
+4WID0.93−0.0340.98.124.801.5913.93
+4WIB1.510.0040.99.105.821.5113.85
+4WID + 4WIB1.01−0.0201.08.604.881.4613.80
Table 6. Summary of the best results for PTC#1, PTC#2, and PTC#3.
Table 6. Summary of the best results for PTC#1, PTC#2, and PTC#3.
PTCActuator
Combinations
X
(m)
Y
(m)
OS%DX
(m)
SX
(m)
MASSA
(deg)
MASSAR
(deg/s)
PTC#1IC#1: FWS1.570.0021.08.984.840.5813.13
IC#2: RWS1.610.0130.99.125.061.6714.73
IC#3: FWS + RWS1.25−0.0030.98.744.581.3613.84
IC#4: FWS + RWS + 4WID1.45−0.0020.99.215.072.0915.30
PTC#2FWS + 4WIB0.52−0.0360.96.662.400.6612.56
FWS + 4WID1.5−0.0120.98.704.530.5814.13
PTC#3FWS + RWS + 4WID0.93−0.0340.98.124.801.5913.93
FWS + RWS + 4WIB1.510.0040.99.105.821.5113.85
Table 7. Velocity gains and weights of LQ objective functions for PTC#2.
Table 7. Velocity gains and weights of LQ objective functions for PTC#2.
Actuator
Combinations
tpξ1ξ2ξ3ξ4ξ5
+4WID0.0650.05390.0050.050.500.03
+4WIB0.0650.07400.0050.050.500.03
+4WID + 4WIB0.0650.06260.0050.050.500.03
+RWS0.0600.07000.0050.100.500.03
+4WID0.0600.08180.0050.100.500.03
+4WIB0.0600.06560.0050.100.500.03
+4WID + 4WIB0.0600.07250.0050.100.500.03
Table 8. Velocity gains and weights of LQ objective functions for PTC#3.
Table 8. Velocity gains and weights of LQ objective functions for PTC#3.
Actuator
Combinations
tpξ1ξ2ξ3ξ4ξ5
+4WID0.070.07340.0100.050.500.03
+4WIB0.070.11720.0150.050.500.03
+4WID + 4WIB0.070.11230.0150.100.500.03
+RWS0.070.09690.0100.100.500.03
+4WID0.070.10120.0100.100.500.03
+4WIB0.070.09350.0100.100.500.03
+4WID + 4WIB0.070.10450.0100.100.500.03
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Park, M.; Yim, S. Comparative Study on Coordinated Control of Path Tracking and Vehicle Stability for Autonomous Vehicles on Low-Friction Roads. Actuators 2023, 12, 398. https://doi.org/10.3390/act12110398

AMA Style

Park M, Yim S. Comparative Study on Coordinated Control of Path Tracking and Vehicle Stability for Autonomous Vehicles on Low-Friction Roads. Actuators. 2023; 12(11):398. https://doi.org/10.3390/act12110398

Chicago/Turabian Style

Park, Manbok, and Seongjin Yim. 2023. "Comparative Study on Coordinated Control of Path Tracking and Vehicle Stability for Autonomous Vehicles on Low-Friction Roads" Actuators 12, no. 11: 398. https://doi.org/10.3390/act12110398

APA Style

Park, M., & Yim, S. (2023). Comparative Study on Coordinated Control of Path Tracking and Vehicle Stability for Autonomous Vehicles on Low-Friction Roads. Actuators, 12(11), 398. https://doi.org/10.3390/act12110398

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

Article Metrics

Back to TopTop