Next Article in Journal
The Prediction of Intrinsically Disordered Proteins Based on Feature Selection
Previous Article in Journal
An Improved Genetic Algorithm for Emergency Decision Making under Resource Constraints Based on Prospect Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Speed Planning and Friction Coefficient Estimation Algorithm for Intelligent Electric Vehicles

School of Mechanical Engineering, Southeast University, Nanjing 211189, China
*
Author to whom correspondence should be addressed.
Algorithms 2019, 12(2), 44; https://doi.org/10.3390/a12020044
Submission received: 16 December 2018 / Revised: 31 January 2019 / Accepted: 12 February 2019 / Published: 20 February 2019

Abstract

:
To improve the safety of intelligent electric vehicles and avoid side slipping on curved roads with changing friction coefficients, an integrated speed planning and friction coefficient estimation algorithm is proposed. With this algorithm, the speeds of intelligent electric vehicles can be planned online using estimated road friction coefficients to avoid lane departures. When a decrease in the friction coefficient is detected on a curved road with a large curvature, the algorithm will plan a low and safe speed to avoid side slipping. When a normal friction coefficient is detected, the algorithm will plan a higher speed for normal driving. Simulations using MATLAB and CarSim have been performed to demonstrate the effectiveness of the designed algorithm. The simulation results suggest that the proposed algorithm is applicable to speed planning on curved roads with changing friction coefficients.

1. Introduction

With the rapid development of vehicle technology [1,2,3], considerable attention has been devoted to intelligent vehicles [4,5,6]. Various technologies have been investigated to satisfy the demanding requirements of commercialization. In particular, lane keeping systems (LKSs) and cruise control are two fundamental technologies for intelligent vehicles. A number of studies have been performed in these critical areas [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22]. A nested proportional-integral-derivative (PID) algorithm with two closed loops was presented for lane keeping in [7]. A super-twisting sliding mode scheme was suggested in [8] to eliminate vehicle lateral tracking error and handle the chatting phenomenon. Three kinds of sliding mode control algorithms were suggested and compared with regard to the actuator saturation scenario in [9]. In [10], a control algorithm consisting of inner and outer loops in the control system was proposed for an LKS based on sliding mode control. In [11], the authors addressed the problem of integrated direct yaw moment control and front steering for an LKS with a transfer function. In [23], a proportional-type space-learning control scheme was constructed for vehicle path tracking, and a stability analysis of the proposed scheme was also performed. In [12], to achieve guaranteed lane keeping and speed tracking performance, a motion controller was developed based on fuzzy model predictive control. In [13], to consider nonlinear vehicle dynamics, a cruise control algorithm was proposed using nonlinear model predictive control. In [14], a self-learning cruise control scheme, in which the system parameters can be adjusted online, was developed for speed tracking. In [15], using a fractional linear feedforward prefilter, a cruise control method was designed considering an acceleration signal. In [16], a cruise controller was constructed considering energy consumption. In [17], the development of cruise control at a motor company was described. In [18], a crone feedback control strategy was extended to a cruise control application. In [19], a super-twisting sliding mode algorithm with adaptive gains was used in the design of a cruise controller. In [20], a cruise controller was proposed based on a fuzzy PID control method. In [21], a piecewise method was used in the development of a cruise control system. In [22], an online learning fuzzy scheme was considered when developing a cruise control system. In [24], an indirect field-oriented control algorithm was proposed for electric vehicles considering the motor uncertainty. Simulations and experiments have shown that the algorithms proposed in these studies are effective for lane keeping or cruise control on normal roads; thus, these investigations have been useful in the development of intelligent vehicles.
Most of these studies have focused on normal roads with constant friction coefficients. However, due to changes in weather or the traffic environment, the friction coefficient may be reduced in real situations. For instance, water or ice on the road surface may significantly decrease the friction coefficient, causing the handling performance of a vehicle on a curved road to degrade. At high speeds, the lane keeping performance may not be maintained due to the limited adhesion on a curved road with a low friction coefficient; however, this concern has seldom been considered in previous studies on lane keeping and cruise control technologies. Thus, it is essential to develop an algorithm that can be used for online speed planning based on estimated friction coefficients. On a curved road with a high curvature and a low friction coefficient, such an algorithm can plan a lower speed to ensure safety.
The aim of this paper is to design an integrated speed planning and friction coefficient estimation algorithm to address changing friction coefficients on curved roads to avoid potential accidents. Two key points should be considered in regard to such road conditions. First, the friction coefficients should be estimated online to provide information about the road surface. Second, the vehicle speed should be planned based on these estimated friction coefficients. On a normal road with a high friction coefficient, a constant desired speed is preferred. When a decrease in the friction coefficient is detected, a safe speed should be planned to avoid side slipping and accidents on curved roads with large curvatures. In this paper, a friction coefficient estimation algorithm is proposed based on torque injection [25,26]. By means of torque injection, the slip ratios of the tires can be increased, which is essential for accurately estimating the friction coefficient. The friction coefficient is estimated based on the slip ratios and side-slip angles over a given sampling period. With the aim of calculating a safe vehicle speed, a speed planning algorithm is proposed that considers a changing friction coefficient on a curved road. The speed profile is planned to prevent potential accidents caused by high speeds and the saturation of friction on curved roads with changing friction coefficients, which is a topic that has seldom been discussed in previous studies.
By combining the proposed speed-planning and friction coefficient estimation algorithms, an integrated speed-planning and friction coefficient estimation (ISPFCE) algorithm is proposed. This algorithm consists of two components, as shown in Figure 1. Note that the sensor system, motion controller, and electric vehicle are also included in this figure to illustrate the overall relationships in the system. The first component of the algorithm is the friction coefficient estimation, in which the friction coefficient is estimated through optimization. The second component is the speed planning, which is designed to calculate a safe speed considering the desired speed and safety constraints. A motion controller is used for lane keeping and cruise control; however, this controller is not the focus of this study. Thus, the designed algorithm can be utilized in combination with various motion controllers based on previously-designed lane-keeping and cruise control algorithms to improve their performance.
The remainder of this paper is organized as follows. In Section 2, the friction estimation algorithm based on torque injection is proposed. In Section 3, the speed-planning method considering the road curvature and friction coefficient is presented. Simulations conducted to verify the proposed algorithm are reported in Section 4. Finally, conclusions are presented in Section 5.

2. Friction Coefficient Estimation

The friction coefficient estimation algorithm, which serves as the basis for speed planning, is proposed in this section. As shown in Figure 2, the proposed algorithm has four components. The first component is the tire force and effective radius estimation, in which the longitudinal, lateral, and vertical forces, as well as the effective radius of each tire are estimated. The second component is the calculation of the tire side-slip angles and slip ratios, which are calculated for each tire. The third component is the calculation of the friction coefficient, which produces the final estimated friction coefficient. The last component is torque injection [25,26], which involves calculating the drive torque to be injected.
In this study, an electric vehicle with four-wheel steering and independent four-wheel drive is considered. The structure of the vehicle used for estimating the friction coefficient is shown in Figure 3. The vertical tire force and side-slip angle of each tire can be calculated using this model. In this figure, l f and l r represent the distances from the front and rear axles, respectively, to the center of gravity (CG). l s is equal to half of the front and rear wheel treads. v x and v y represent the longitudinal and lateral velocities, respectively, at the CG. ω r is the yaw rate. The subscripts “ f l ”, “ f r ”, “ r l ”, and “ r r ” denote the values corresponding to the front left, front right, rear left, and rear right tires, respectively. F x i and α i ( i = f l , f r , r l , and r r ) are the longitudinal tire force and side-slip angle, respectively, of the corresponding tire. δ f and δ r are the steering angles of the front and rear wheels, respectively.

2.1. Tire Force and Effective Radius Estimation

The estimated vertical load force on each tire can be expressed as [27,28,29]:
F ^ z f l = m g l r 2 L m a x h c 2 L m a y h c 2 l s F ^ z f r = m g l r 2 L m a x h c 2 L + m a y h c 2 l s F ^ z r l = m g l f 2 L + m a x h c 2 L m a y h c 2 l s F ^ z r r = m g l f 2 L + m a x h c 2 L + m a y h c 2 l s ,
where a x and a y represent the longitudinal and lateral accelerations, respectively, at the CG, which can be obtained by means of integrated GPS/INS sensors [30]; L = l f + l r ; h c and m are the height of the CG and the total mass, respectively; and g denotes the acceleration due to gravity.
Under normal conditions, the effective radius of each tire is mainly determined by the vertical force on that tire and can be estimated as [31]:
R ^ w i = R g d r i 3 ,
where i = f l , f r , r l , and r r ; R ^ w i is the estimated effective radius of the corresponding tire; R g is the geometric radius of the tire; and d r i = F ^ z i k z is the deflection of the tire under the vertical force acting on it, where k z is the vertical stiffness of the tire.
Based on the dynamics of each wheel, the longitudinal force on each wheel can be estimated as:
F ^ x i = T i J w ω ˙ i R ^ w i + F ^ z i f r ,
where i = f l , f r , r l , and r r ; J w and f r are the rotational inertia and rolling resistance coefficient, respectively, of each wheel, which are constant factors; and T i and ω ˙ i are the drive torque and angular acceleration, respectively, of the corresponding wheel. The drive torque can be obtained from the motor control system, and the angular acceleration can be obtained from the wheel angular acceleration sensors.
Due to the coupling of the lateral forces on the tires on the same axle, it is difficult to calculate the lateral force on each tire separately. However, the lateral forces on the front and rear axles can be estimated based on the vehicle dynamics [32]:
F ^ y f = 1 l f + l r ( l f m a y I z ω ˙ r ) F ^ y r = 1 l f + l r ( l r m a y I z ω ˙ r ) ,
where F ^ y f and F ^ y r are the estimated lateral forces on the front and rear axles, respectively, and I z denotes the moment of inertia around the vertical axle.

2.2. Calculation of the Tire Side-Slip Angles and Slip Ratios

The speed and direction of each wheel can be calculated as follows [33,34,35,36,37]:
v x f l = v x ω r l s cos δ f + v y + ω r l f sin δ f v x f r = v x + ω r l s cos δ f + v y + ω r l f sin δ f v x r l = v x ω r l s cos δ r + v y ω r l r sin δ r v x r r = v x + ω r l s cos δ r + v y ω r l r sin δ r
and:
θ f l = arctan v y + ω r l f v x ω r l s θ f r = arctan v y + ω r l f v x + ω r l s θ r l = arctan v y ω r l r v x ω r l s θ r r = arctan v y ω r l r v x + ω r l s ,
where v x i and θ i ( i = f l , f r , r l , and r r ) are the speed and direction, respectively, of each wheel. The longitudinal speed, lateral speed, and yaw rate can be obtained by means of integrated GPS/INS sensors [30]. The front and rear steering angles can be obtained from angle sensors in the steering system.
The side-slip angle and slip ratio of each tire can be expressed as follows:
α i = δ i + θ i
and:
s i = v x i ω i R ^ w i max v x i , ω i R ^ w i ,
where i = f l , f r , r l , and r r and ω i is the rolling speed of the corresponding wheel, which can be obtained from an angular speed sensor on that wheel.

2.3. Calculation of the Friction Coefficient

To obtain a reliable estimate of the friction coefficient, the characteristics of a tire are analyzed based on the Pacejka tire model [38], and the effectiveness of the corresponding friction coefficient estimation is discussed. Then, an optimization-based estimation algorithm is proposed.

2.3.1. Analysis of the Tire Model

The Magic tire model [38] is a commonly-used tire model in vehicle dynamics investigations, and it has also been found to be suitable for experimental applications. Thus, in this study, the Magic tire model is used to estimate the friction coefficient:
F x = D x sin C s arctan B x s E x B x s arctan B x s F y = D y sin C y arctan B y α E B y α arctan B y α ,
where B x , C x , D x , E x , B y , C y , D y , and E y are the parameters of the tire model. The above functions constitute a common simplified description of the actual Magic tire model. The parameters of the tire model are mainly affected by the vertical force and the friction coefficient [38]. However, the friction coefficient does not explicitly appear in Equation (9), which makes this model difficult to use for friction coefficient estimation. The structure of the Magic tire model is visualized in Figure 4. This model consists of four parts [38], labeled from A–D. Parts A and C are used to calculate the longitudinal and lateral forces under conditions of pure longitudinal slip and pure side slip, respectively. Parts B and D are used to calculate the longitudinal and lateral forces, respectively, under combined slip conditions based on the calculation results from Parts A and C. The corresponding detailed calculations can be found in [38].
To facilitate the estimation of the friction coefficient, the Magic tire model can be further rewritten as follows:
F x = Magic x F z , s , α , μ F y = Magic y F z , s , α , μ ,
where F x and F y are the longitudinal and lateral tire forces, respectively. The functions defined above provide a simple way to represent the Magic tire model for the purpose of friction coefficient estimation. The structures of these functions are visualized in Figure 5 and Figure 6, respectively.
For friction coefficient estimation, the vertical force, side-slip angle, and slip ratio are known and are viewed as parameters instead of inputs in the defined functions. Once the vertical force, side-slip angle, and slip ratio of a tire have been obtained, the relationships between F x , F y and μ can be further rewritten as follows:
F x = Magic x , i μ F y = Magic y , i μ ,
where i = f l , f r , r l , and r r and the functions Magic x , i and Magic y , i represent the relationships between F x , F y , and μ when the vertical force, side-slip angle, and slip ratio are known for tire i. The structures of the functions defined above are illustrated in Figure 7 and Figure 8, respectively.
To facilitate the discussion of the estimation of the friction coefficient, the tire force use ratio is defined as follows:
μ u = F x 2 + F y 2 F z .
The relationships between the tire force use ratio and the slip ratio and side-slip angle with different friction coefficients are plotted in Figure 9 and Figure 10. Under conditions of a low slip ratio or low side-slip angle, the friction coefficient has only a small effect on the tire force use ratio, as shown in Figure 9. Thus, the friction coefficient also exerts only a small effect on the tire force under such conditions, which causes friction coefficient estimation based on the tire force to be unreliable. For instance, if the slip ratio and side-slip angle are both zero, then the longitudinal and lateral forces are also zero regardless of the friction coefficient. Indeed, the friction coefficient has no effect on the longitudinal and lateral forces under such conditions. Therefore, it is impossible to estimate the friction coefficient accurately based on the longitudinal and lateral tire forces under conditions of a low slip ratio and a low side-slip angle.
To measure the effectiveness of the estimation, an effectiveness factor is defined as follows:
μ e = Magic x , i μ + Δ μ Magic x , i μ 2 + Magic y , i μ + Δ μ Magic y , i μ 2 / F z Δ μ ,
where μ e is the effectiveness factor, which indicates the effectiveness of the estimation, and Δ μ is a small number. μ e can be viewed as reflecting the effect of the friction coefficient on the tire force use ratio. When μ e is small, a change in the friction coefficient has only a small effect on the tire force, whereas a large value of μ e indicates that this effect is large. Thus, friction coefficient estimation based on the tire force is inaccurate when μ e is small. Because the estimated friction coefficients serve as the basis for speed planning, unreliable friction coefficient estimates should not be used. Only reliable estimates should be applied for speed planning.

2.3.2. Process of Estimation

To obtain a reliable estimate of the friction coefficient, the following estimation process is considered. The longitudinal tire force can be estimated from Equation (3), and the lateral force on each axle can be estimated from Equation (4). Because the lateral force on each axle is the total lateral force on both tires on the same axle, under ideal conditions, the following equations should be satisfied:
Magic x , f l μ ^ f = F ^ x f l Magic x , f r μ ^ f = F ^ x f r
and:
Magic y , f l μ ^ f + Magic y , f r μ ^ f = F ^ y f .
To facilitate the calculation, the following optimization problem is defined to estimate the friction coefficient of each front tire:
min : F M ( μ ^ f ) = w x ( Magic x , f l μ ^ f F ^ f l ) 2 + w x ( Magic x , f r μ ^ f F ^ f r ) 2 + w y ( Magic y , f l μ ^ f + Magic y , f r μ ^ f F ^ y f ) 2
where μ ^ f is the estimated friction coefficient of each front tire. The meaning of Equation (16) is that if the estimated friction coefficient is equal to the actual value, then F M ( μ ^ f ) = 0 . Thus, the friction coefficient of the tires on the front axle can be estimated by means of the optimization problem defined in Equation (16). The friction coefficient of the tires on the rear axle can also be estimated using the same method.
At each sampling instant, two estimated friction coefficients and their corresponding effectiveness factors can be obtained according to Equations (13) and (16). Due to the uncertainty of the model and the error in the numerical computations, the values from several successive sampling instants are used for estimation. The pseudocode for obtaining a reliable estimate of the friction coefficient is presented in Algorithm 1. In this algorithm, n t is the total number of sampling instants used for estimation in each sampling window, m e is the total number of effective estimates, μ e U is the threshold for accepting an estimated friction coefficient as reliable, μ e L is the threshold for accepting a friction coefficient estimate from a single sampling instant for use in the final estimation, and μ ^ is the final output of the friction coefficient estimation process.
Algorithm 1 Friction coefficient estimation.
Input: μ e i , j , μ ^ i , j ( i = f and r; j = 1 , 2 , , n t )
Output: μ ^
1: procedure
2:   if max ( μ e f , 1 , μ e r , 1 ) > μ e U then
3:      μ ^ s = 0
4:      m e = 0
5:     for j = 1 n t do
6:       for i = f , r do
7:         if μ e i , j > μ e L then
8:            μ ^ s = μ ^ s + μ ^ i , j
9:            m e = m e + 1
10:         end if
11:       end for
12:     end for
13:      μ ^ = μ ^ s / m e
14:     return μ ^
15:   else
16:    A reliable estimate cannot be found in this sampling window.
17:   end if
18: end procedure

2.4. Torque Injection

The relationship between the effectiveness factor and the tire slip ratio is illustrated in Figure 11. As shown in this figure, as the slip ratio increases, the effectiveness factor also increases. Thus, the friction coefficient estimation is more reliable with a higher slip ratio. By contrast, when a vehicle is traveling on a flat road at a constant speed, the tire slip ratio of each tire is low, which may lead to inaccurate estimation, as discussed previously. Thus, to obtain reliable estimates, the torque injection method is applied for friction coefficient estimation. The front wheels are injected with drive torques, and the rear wheels are injected with equal braking torques. The slip ratios of the front and rear wheels become larger as the injected torques increase. Thus, the estimation of the friction coefficient becomes more reliable as the slip ratio of each tire increases.
The torque injection method can be used to force tires to work at high slip ratios while having only a small effect on the longitudinal dynamics of the vehicle [25]. However, note that torque injection also leads to tire wear and power loss. Thus, the need for reliable estimation should be balanced with the disadvantages of torque injection. The principles of and procedure for the torque injection method are illustrated in Figure 12 and Figure 13, respectively. In these figures, T c , f and T c , r denote the normal front and rear drive torques, respectively, without torque injection; F x , s u m denotes the resultant longitudinal drive force on the four tires; T i n j is the injected drive torque; and T f and T r are the actual drive torques on the front and rear wheels, respectively. As indicated in Figure 12 and Figure 13, in the torque injection method, the front wheels are injected with drive torques, and the rear wheels are injected with equal braking torques. Thus, the total injected drive torque has no effect on the vehicle’s motion, but higher slip ratios of the front and rear wheels can be obtained, which are essential for accurately estimating the friction coefficient.
In this paper, an adaptive torque injection method is proposed for friction coefficient estimation. The torque injection signal is illustrated in Figure 14. In this figure, t c y c denotes the constant length of the torque injection cycle. A shorter cycle improves the update frequency for road friction estimation, but increases the power loss, and vice versa. t 0 denotes the start time of each cycle. t e f f denotes the time when a reliable friction coefficient is obtained in each cycle.
The torque injected in each cycle is calculated as follows:
T i n j = K T · t i n j t i n j = t t 0 ,
where t is the current time, T i n j is the injected torque, and K T is a constant parameter used to calculate T i n j . At the beginning of each cycle, the injected torque is zero.
As the injected torque increases, the slip ratios of the front and rear wheels also increase. Thus, the effectiveness factor will also increase. Once the effectiveness factor reaches the specified threshold, a reliable estimate of the friction coefficient can be obtained; the time at which this occurs is denoted by t e f f . The injected torque is then set to zero for the time remaining in the current cycle, as indicated in Figure 14. t i n j and T i n j may be different from one cycle to another. The torque injected into the front wheels is T i n j . The torque injected into the rear wheels is T i n j .

3. Speed Planning

Generally, it is preferable to control a vehicle such that it moves at a constant desired speed for the purposes of comfort and economy. However, for a vehicle on a slippery road with a large curvature, traveling at a high speed may lead to a situation that could cause a serious traffic accident, such as lane departure. Thus, the vehicle’s speed should be planned online in accordance with the estimated friction coefficient and road shape.
In this section, an optimization-based algorithm for speed planning is proposed. First, the optimization variables, constraints, and objective function are designed for the speed planning task. Then, the optimization problem is constructed. Finally, a method of calculating the optimal safe speed is presented. Speed planning plays an important role in the vehicle control system, which serves as the link between the friction coefficient estimation and the motion controller.

3.1. Optimization Variables

To simplify the speed-planning problem, the path in front of the considered vehicle is separated into several segments of the same length by means of nodes. These segments and nodes are called path segments and path nodes for convenience. The speeds of the vehicle at each end of each of these segments are subjected to optimization. The travel speeds between the ends of each segment can be calculated through linear interpolation. For each segment, the path node at which the vehicle enters the segment is called the start node of the segment, and the path node at which the vehicle departs from the segment is called the end node of the segment. Note that a path can be separated into any number of segments for speed planning. However, the more segments are used for speed planning, the greater is the computational burden. Because the friction coefficient may change along the path, it is inappropriate to plan the entire path along the vehicle’s route simultaneously. To simplify the problem, the vehicle’s speed is planned successively for road sections of a certain length in front of the vehicle. As an example, the path in front of a vehicle is shown in Figure 15. The section of the path within a certain distance in front of the vehicle can be divided into n segments by n + 1 nodes, as shown in Figure 16. These nodes are labeled from p 0 p n . Note that p 0 is located at the CG of the vehicle. The path segment between p i 1 and p i is denoted by s i . Each of these path segments has a length of l p .
The planned speeds at p i ( i = 1 , 2 , , n ) are represented by v i ( i = 1 , 2 , , n ), respectively. Under the assumption that the speed of the vehicle varies linearly between each pair of adjacent nodes, the planned speeds of the vehicle along the path are further illustrated in Figure 17. As shown, once the planned speeds at these nodes have been determined, the planned speed profile along the path can be obtained. The error between the planned speed at each node and the desired speed is defined as:
e i = v i v d ,
where v d is the desired speed. Clearly, if e i = 0 for all i = 1 , 2 , , n , then the planned speed is equal to the desired speed along the entire path. Once e i ( i = 1 , 2 , , n ) has been determined, v i ( i = 1 , 2 , , n ) can be calculated as:
v i = v d + e i .
As shown in Figure 18, the planned speeds can also be determined from the speed errors at the nodes. Thus, the speed errors at the nodes, denoted by e = e 0 e 1 e 2 e n 1 e n T , are defined as the decision variables for speed planning. Note that because the speed of the vehicle cannot track a step reference, the speed at node p 0 is not designated as one of the decision variables.

3.2. Constraints

This subsection presents the speed-planning constraints. First, constraints on the speed on a path segment are designed by considering several dangerous driving conditions. Then, constraints on the speed at a node are discussed. Finally, constraints for speed planning along the entire considered path are proposed.

3.2.1. Constraints on the Speed on a Path Segment

Consider the speed plan for a path segment. Traveling at an excessively high speed on the segment may lead to a skid or rollover. To avoid a skid, the lateral acceleration should satisfy:
a y = v 2 R m i n μ g ,
where R m i n is the minimum radius of the path segment. Thus, a constraint on the planned speed as determined by the lateral acceleration can be defined as:
v S m a x = f s S μ g R m i n ,
where f s S < 1 is a safety factor for avoiding skid conditions.
To avoid rollover, the lateral acceleration should satisfy:
m g l s m a y h c = m v 2 R m i n h c .
The limit on the speed for avoiding rollover can thus be defined as:
v R m a x = f s R g l s R m i n h c ,
where f s R < 1 is a safety factor for avoiding rollover conditions.
Thus, the maximum safe speed on a path segment can be expressed as:
v s m a x = min v S m a x , v R m a x .

3.2.2. Constraints on the Speed at a Node

Generally, each node is adjacent to two path segments. To ensure safety, the maximum allowable planned speed at a node should be equal to the minimum value between the maximum safe speeds on the two adjacent segments, as expressed below:
v n m a x , i = min v s m a x , i , v s m a x , i + 1 ,
where v s m a x , L and v s m a x , R are the maximum safe speeds on the segments preceding and following the node, respectively. For the last node p n , only the constraints on the path segment s n need to be considered.
Due to the limits of the road friction, the longitudinal acceleration should be limited to avoid tire slip. Because the lateral velocity of the vehicle is small, the derivative of the speed is simply used to approximate the longitudinal acceleration:
v ˙ x a x f s L · μ · g ,
where f s L < 1 is a safety factor for limiting the longitudinal acceleration of the vehicle.
For a given path segment, the travel time of the vehicle on that segment is limited by the maximum safe speed on that segment; this limitation can be expressed as:
t t < l p v s m a x .
The planned speed profile along a segment is obtained through linear interpolation from the planned speeds at the end points. Thus, the derivative of the speed is constant in each segment. A constraint on the change in speed in a segment can be calculated as follows:
| v e v s | t t · f s L · μ · g ,
where v s and v e are the planned vehicle speeds at the start and end points, respectively.
Equation (28) can be further rewritten as:
v d e c v e v s v d e c ,
where v d e c = t t f s L · μ · g .

3.2.3. Constraints on the Planned Speed Profile along the Path

In summary, according to Equations (26)–(29), the constraints on the planned speeds at all nodes can be expressed as:
v m i n v 1 v n m a x , 1 v m i n v 2 v n m a x , 2 v m i n v n v n m a x , n
and:
v d e c 1 v 1 v 0 v d e c 1 v d e c 2 v 2 v 1 v d e c 2 v d e c n v n v n 1 v d e c n .
By substituting Equation (19) into Equations (30) and (31), the speed-planning constraints can be rewritten as:
v m i n v d e 1 v n m a x , 1 v d v m i n v d e 2 v n m a x , 2 v d v m i n v d e n v n m a x , n v d
and:
v d e c 1 + e 0 e 1 v d e c 1 + e 0 v d e c 2 e 2 e 1 v d e c 2 v d e c n e n e n 1 v d e c n .
Equations (32) and (33) can be further rewritten as:
b L Ae b U b l b e b u b ,
where e = e 1 e 2 e n , b L = v m i n v d v m i n v d v m i n v d , b U = v n m a x , 1 v d v n m a x , 2 v d v n m a x , 3 v d , b l b = v d e c 1 + e 0 v d e c 2 v d e c n , b u b = v d e c 1 + e 0 v d e c 2 v d e c n , and A = 1 0 0 0 0 1 1 0 0 0 0 1 1 0 0 0 0 0 1 0 0 0 0 1 1 .

3.3. Objective Function

To reduce the error between the desired and planned speeds, the following objective function is defined:
F e = 1 2 e T He ,
where H = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 . The optimal planned speed can be obtained by minimizing F ( e ) .

3.4. Construction of the Optimization Problem

Considering the discussed constraints, the optimization problem for planning the optimal safe speed profile is constructed as follows:
min : F ( e )
subject to:
b L Ae b U b l b e b u b .
The optimization problem defined above is a typical quadratic programming problem, which can be solved with qpOASES [39]. By solving this optimization problem, the speed errors at all nodes can be obtained. The planned speed at each node along the path can then be calculated in accordance with Equation (19). Using the conventional linear interpolation method, the planned speed profile on each segment along the path can also be obtained.

3.5. Calculation of the Planned Speed

Once the optimal speed planning solution has been found, the planned speed at any point can be calculated through simple linear interpolation. Consider a path section with a single segment, as shown in Figure 19. p s and p e are the start and end points, respectively, of this segment. p v is the position of the vehicle. l v is the length of the path between p s and p v . s s is the length of the path between p s and p e . The planned speed of the vehicle at p v can be calculated as:
v p = l v l p v s + l p l v l p v e ,
where v s and v e are the planned speeds at p s and p e , respectively.

4. Numerical Simulation

To verify the effectiveness of the proposed ISPFCE algorithm, a numerical simulation was conducted with CarSim and Simulink. The vehicle model in CarSim is a 27-degree-of-freedom full-vehicle model that is equipped with nonlinear tire models. The measurement noise characteristics of realistic sensors (such as those in the RT4000 navigation system from Oxford Technical Solutions Ltd., Bicester, U.K.) were reproduced in the simulation. The proposed algorithm was implemented using Simulink. The simulation parameters are listed in Table 1.
In the simulation, the initial speed of the vehicle at the origin point, which was also considered to be the desired speed, was set to 23 m/s. The shape and curvature of the road path are shown in Figure 20 and Figure 21, respectively. The minimum radius of the path was 187.5 m ( 1 / R = 0.0053(1/m)), as shown in Figure 21. The friction coefficient was initially 0.85; subsequently, it decreased to 0.2 and then returned to 0.85 in a stepwise manner, as shown in Figure 21.
As discussed previously, the proposed ISPFCE algorithm can be applied in combination with previously-investigated motion controllers. In the simulation, two motion controllers were designed for comparison, one based on fuzzy model predictive control [12] (Controller A) and one based on a linear quadratic regulator [40] (Controller B). The performances of these motion controllers with and without the proposed ISPFCE algorithm are presented for comparison.
The lane-keeping performance and speeds of vehicles controlled by Controllers A and B with and without ISPFCE are illustrated in Figure 22, Figure 23 and Figure 24. As shown in Figure 23, without the proposed ISPFCE algorithm, Controllers A and B attempt to control the corresponding vehicles such that they travel at the desired speed along the entire road. However, as discussed above, due to the limited road friction at a tight bend, it is impossible for the vehicles to pass through the bend at the desired speed. Figure 22 indicates that the vehicles controlled by Controllers A and B cannot remain on the path at the bend, meaning that the vehicles cannot stay within their lane on the road. The large lateral position errors shown in Figure 24 indicate that the vehicles leave the road during the simulation.
By contrast, with the proposed ISPFCE algorithm, the speeds of the vehicles controlled by Controllers A and B are reduced to low values before the tight bend with the corresponding low friction coefficient, as indicated in Figure 23. Thus, these vehicles can move along the reference path with small tracking errors, as shown in Figure 22 and Figure 24. This comparison demonstrates that the proposed ISPFCE algorithm can be applied with different motion controllers to handle roads with changes in their friction coefficients. Notably, traveling at a low speed at all times would lead to increased time consumption and decreased traffic flow. Thus, once a high friction coefficient is again detected, the vehicle speed should be increased. As shown in Figure 23, the vehicle speeds return to the desired value once a high friction coefficient is again detected, demonstrating that the proposed ISPFCE algorithm can handle various changes in the friction coefficient on a curved road.
The control inputs to the vehicles controlled by Controllers A and B with ISPFCE are shown in Figure 25 and Figure 26. The details of the friction estimation and speed planning are further plotted in Figure 27, Figure 28, Figure 29, Figure 30, Figure 31, Figure 32, Figure 33 and Figure 34. As shown in Figure 21, the road friction coefficient changes from 0.85 to 0.2 at approximately 550 m and changes back to 0.85 at approximately 2900 m. Due to the time elapsed during each torque injection cycle, there is some delay before the estimated friction coefficient reaches the real value, as illustrated in Figure 28. The results of online speed planning are presented in Figure 29. As shown in Figure 31, the initial speed plan is Plan A. Because the estimated friction coefficient is high, the ISPFCE algorithm assumes that the vehicle can move at the desired high speed along the path. Once a decrease in the friction coefficient is detected, the ISPFCE algorithm determines that the desired high speed is not safe on the curved road. Thus, speed Plan B is generated, as shown in Figure 30. The relationship between speed Plan B and the road curvature is plotted in Figure 32. As shown in Figure 32, the planned vehicle speed is reduced to a low value before the tight bend. Thus, the vehicle can follow the path safely. While the vehicle is still on a section of the road with low friction, it should follow the path at a low speed. Thus, speed Plan C is generated. As shown in Figure 33, the planned speed is very low on these low-friction bends. When a high friction coefficient is again detected, the ISPFCE algorithm determines that it is safe for the vehicle to again move along the reference path at the desired speed. Thus, speed Plan D is generated, as plotted in Figure 34, in which the planned speed is increased to the desired value to satisfy the set requirements.
These simulation results demonstrate that the proposed ISPFCE algorithm can handle changes in the friction coefficient on a curved road. Once a decrease in the friction coefficient that may lead to dangerous conditions is detected, a low, safe vehicle speed can be planned to pass through the tightly-curved road section. When the friction coefficient again increases to a high value, the ISPFCE algorithm can plan a return to the desired speed along the road. Moreover, these simulation results also demonstrate the flexibility of using the ISPFCE algorithm in combination with different motion controllers for general applications.

5. Conclusions

In this paper, an integrated speed-planning and friction coefficient estimation (ISPFCE) algorithm is proposed. The friction coefficient is estimated online by means of torque injection. If a decrease in the friction coefficient is detected before a tight bend, the speed can be adjusted to a safe value before an accident could potentially occur. A simulation based on CarSim and Simulink demonstrates the effectiveness of the proposed algorithm. The present paper reveals that a sudden decrease in the friction coefficient is dangerous for vehicles on curved roads, potentially leading to serious traffic accidents. Conventional lane-keeping systems and cruise control systems without friction estimation may fail to handle such a decrease in the friction coefficient properly, particularly on a curved road with a large curvature. One feasible solution is to estimate the friction coefficient online and reduce the speed before such a bend with a low friction coefficient, which is the key purpose of the algorithm proposed in this paper.
As discussed above, speed-planning and friction coefficient estimation are critical for ensuring the safety of vehicles. In the proposed ISPFCE algorithm, only speed planning for the host vehicle is considered. In fact, however, it is also useful to estimate the friction coefficients and planned safe speeds for other vehicles. Thus, given a suitable vehicular network, it is also essential to send relevant information to other vehicles that cannot estimate the friction coefficient for planning safe speeds. The question of how to use the information collected by the host vehicle to improve the safety of other nearby vehicles will be investigated in the future. Due to the practical limitations concerning available equipment and experimental sites, only a numerical simulation is reported in this study. Further experimental investigations will also be conducted to improve the performance of the ISPFCE algorithm.

Author Contributions

Conceptualization, C.B.; methodology, C.B.; validation, C.B., T.Z., G.Y. and L.X.; writing, original draft preparation, C.B. and T.Z.; writing, review and editing, T.Z. and L.X.; supervision, G.Y.

Funding

This research was funded by the National Natural Science Foundation of China, Grant Numbers 51575103 and U1664258; the National Key R&D Program in China, Grant Numbers 2016YFB0100906 and 2016YFD0700905; the Foundation of State Key Laboratory of Automotive Simulation and Control, Grant Number 20160112; and the Scientific Research Foundation of Graduate School of Southeast University and Southeast University Excellent Doctor Degree Thesis Training Fund, Grant Numbers YBJJ1703 and YBJJ1704.

Acknowledgments

The authors wish to express their gratitude to the anonymous reviewers and the Associate Editor for their rigorous comments during the review process.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ISPFCEintegrated speed planning and friction coefficient estimation
CGcenter of gravity

References

  1. D’Ambrosio, C.; Sbarra, G.; Tiberti, M.; Verrelli, C.M.; Consolini, L. A new spatial learning control for autonomous vehicles: Experimental results. In Proceedings of the Annual Conference on System of Systems Engineering, Paris, France, 19–22 June 2018; pp. 68–73. [Google Scholar]
  2. Antonelli, D.; Pasquale, L.; Salvatore, A.; Tiberti, M.; Verrelli, C.M. Electric vehicles under slip constraints: Experimental results. In Proceedings of the AEIT International Annual Conference, Bari, Italy, 3–5 October 2018; pp. 1–6. [Google Scholar]
  3. Marino, R.; Pasquale, L.; Scalzi, S.; Verrelli, C.M. Automatic Rotor Speed Reference Generator for Electric Vehicles Under Slip Constraints. IEEE Trans. Intell. Transp. Syst. 2015, 16, 3473–3478. [Google Scholar] [CrossRef]
  4. Sc., A.E.D.; Arem, B.V.; Amditis, A.; Lytrivis, P.; Portouli, E.; Papp, Z.; Schutter, B.D.; Ploeg, J.; Baskar, L.D.; Naus, G. Handbook of Intelligent Vehicles; Springer: Berlin, Germany, 2012; pp. 1217–1267. [Google Scholar]
  5. Wang, J.; Rajamani, R. Should adaptive cruise-control systems be designed to maintain a constant time gap between vehicles? IEEE Trans. Veh. Technol. 2004, 53, 1480–1490. [Google Scholar] [CrossRef]
  6. Qu, T.; Chen, H.; Cao, D.; Guo, H. Switching-based stochastic model predictive control approach for modeling driver steering skill. IEEE Trans. Intell. Transp. Syst. 2015, 16, 365–375. [Google Scholar] [CrossRef]
  7. Marino, R.; Scalzi, S.; Netto, M. Nested PID steering control for lane keeping in autonomous vehicles. Control Eng. Pract. 2011, 19, 1459–1467. [Google Scholar] [CrossRef]
  8. Tagne, G.; Talj, R.; Charara, A. Higher-Order Sliding Mode Control for Lateral Dynamics of Autonomous Vehicles, with Experimental Validation. In Proceedings of the Intelligent Vehicles Symposium, Gold Coast, QLD, Australia, 23–26 June 2013; pp. 678–683. [Google Scholar]
  9. Yang, J.; Zhao, J.; Lei, J. Research on three kinds of lane keeping method based on anti-saturation and adaptive method. Int. J. Hybrid Inf. Technol. 2015, 8, 317–326. [Google Scholar] [CrossRef]
  10. Lei, J.; Wu, H.; Yang, J.; Zhao, J. Sliding mode lane keeping control based on separation of translation and rotation movement. Opt. Int. J. Light Electron Opt. 2016, 127, 4369–4374. [Google Scholar] [CrossRef]
  11. Raksincharoensak, P.; Nagai, M.; Shino, M. Lane keeping control strategy with direct yaw moment control input by considering dynamics of electric vehicle. Veh. Syst. Dyn. 2006, 44, 192–201. [Google Scholar] [CrossRef]
  12. Bian, C.; Yin, G.; Zhang, N.; Xu, L. Takagi-sugeno fuzzy model predictive controller design for combining lane keeping and speed tracking of four wheels steering and four wheels drive electric vehicle. In Proceedings of the Chinese Control and Decision Conference, Chongqing, China, 28–30 May 2017; pp. 4067–4072. [Google Scholar]
  13. Luo, Y.; Chen, T.; Li, K. Nonlinear model predictive cruise control of hybrid electric vehicle. Chin. J. Mech. Eng. 2015, 51, 11–21. [Google Scholar] [CrossRef]
  14. Wang, J.; Xu, X.; Liu, D.; Sun, Z.; Chen, Q. Self-learning cruise control using kernel-based least squares policy iteration. IEEE Trans. Control Syst. Technol. 2014, 22, 1078–1087. [Google Scholar] [CrossRef]
  15. Taymans, A.; Melchior, P.; Malti, R.; Aioun, F.; Guillemard, F.; Servel, A. Cruise control of an electric vehicle through fractional linear feedforward prefiltering of an acceleration reference signal. IFAC-PapersOnLine 2017, 50, 12569–12574. [Google Scholar] [CrossRef]
  16. Schwickart, T.; Voos, H.; Hadji-Minaglou, J.R.; Darouach, M.; Rosich, A. Design and simulation of a real-time implementable energy-efficient model-predictive cruise controller for electric vehicles. J. Frankl. Inst.-Eng. Appl. Math. 2015, 352, 603–625. [Google Scholar] [CrossRef]
  17. Eichelberger, A.H.; McCartt, A.T. Toyota drivers’ experiences with dynamic radar cruise control, pre-collision system, and lane-keeping assist. J. Saf. Res. 2016, 56, 67–73. [Google Scholar] [CrossRef] [PubMed]
  18. Morand, A.; Moreau, X.; Melchior, P.; Moze, M.; Guillemard, F. CRONE cruise control system. IEEE Trans. Veh. Technol. 2016, 65, 15–28. [Google Scholar] [CrossRef]
  19. Reichhartinger, M.; Horn, M. Robust cascaded automatic cruise control of electric vehicles. Int. J. Robust Nonlinear Control 2016, 26, 662–680. [Google Scholar] [CrossRef]
  20. Qiu, C.; Liu, C.; Shen, F.; Chen, J. Design of automobile cruise control system based on Matlab and fuzzy PID. Trans. Chin. Soc. Agric. Eng. 2012, 28, 197–202. [Google Scholar]
  21. Lee, Y.; Ryoo, Y.J. Design of cruise control system for electric vehicle using piece-wised control. Adv. Intell. Syst. Comput. 2014, 269, 37–44. [Google Scholar] [CrossRef]
  22. Onieva, E.; Godoy, J.; Villagra, J.; Milanes, V.; Perez, J. On-line learning of a fuzzy controller for a precise vehicle cruise control system. Expert Syst. Appl. 2013, 40, 1046–1053. [Google Scholar] [CrossRef] [Green Version]
  23. Consolini, L.; Verrelli, C.M. Learning control in spatial coordinates for the path-following of autonomous vehicles. Automatica 2014, 50, 1867–1874. [Google Scholar] [CrossRef]
  24. Marino, R.; Scalzi, S.; Verrelli, C.M. Speed and traction control in electric vehicles with induction motors. In Proceedings of the Electrical Systems for Aircraft, Railway and Ship Propulsion, Bologna, Italy, 16–18 October 2012. [Google Scholar]
  25. Chen, Y.; Wang, J. Adaptive vehicle speed control with input injections for longitudinal motion independent road frictional condition estimation. IEEE Trans. Veh. Technol. 2011, 60, 839–848. [Google Scholar] [CrossRef]
  26. Deur, J.; Pavković, D.; Burgio, G.; Hrovat, D. A model-based traction control strategy non-reliant on wheel slip information. Veh. Syst. Dyn. 2011, 49, 1245–1265. [Google Scholar] [CrossRef]
  27. Xiong, L.; Yu, Z.; Wang, Y.; Yang, C.; Meng, Y. Vehicle dynamics control of four in-wheel motor drive electric vehicle using gain scheduling based on tyre cornering stiffness estimation. Veh. Syst. Dyn. 2012, 50, 831–846. [Google Scholar] [CrossRef]
  28. Zhenta, L. Optimal allocation of tire force based on vehicle stability. J. Mech. Electr. Eng. 2015, 32, 433–438. [Google Scholar] [CrossRef]
  29. Osborn, R. TaehyunShim. Independent control of all-wheel-drive torque distribution. Veh. Syst. Dyn. 2006, 44, 529–546. [Google Scholar] [CrossRef]
  30. Zhang, H.; Huang, X.; Wang, J.; Karimi, H.R. Robust energy-to-peak sideslip angle estimation with applications to ground vehicles. Mechatronics 2015, 30, 338–347. [Google Scholar] [CrossRef]
  31. Jazar, R.N. Vehicle Dynamics: Theory and Application; Springer: Berlin, Germany, 2008. [Google Scholar]
  32. Villagra, J.; D’Andréanovel, B.; Fliess, M.; Mounier, H. A diagnosis-based approach for tire-road forces and maximum friction estimation. Control Eng. Pract. 2011, 19, 174–184. [Google Scholar] [CrossRef]
  33. Wang, J.; Longoria, R.G. Coordinated vehicle dynamics control with control distribution. In Proceedings of the American Control Conference, Minneapolis, MN, USA, 14–16 June 2006; pp. 5348–5353. [Google Scholar]
  34. Bian, C.; Yin, G.; Hu, M.; Ning, Z.; Nan, C. Fault tolerant control design for all wheels independent drive and steering electric vehicle with singular effective drive wheel. In Proceedings of the Vehicle Power & Propulsion Conference, Hangzhou, China, 17–20 October 2016. [Google Scholar]
  35. Guvenc, B.; Acarman, T.; Guvenc, L. Coordination of steering and individual wheel braking actuated vehicle yaw stability control. In Proceedings of the IEEE Intelligent Vehicles Symposium, Columbus, OH, USA, 9–11 June 2003; pp. 288–293. [Google Scholar]
  36. Wang, R.; Wang, J. Fault-tolerant control for electric ground vehicles with independently-actuated in-wheel motors. J. Dyn. Syst. Meas. Control Trans. ASME 2012, 134, 021014. [Google Scholar] [CrossRef]
  37. Feiqiang, L.; Jun, W.; Zhaodu, L. On the vehicle stability control for electric vehicle based on control allocation. In Proceedings of the IEEE Vehicle Power and Propulsion Conference, Harbin, China, 3–5 September 2008; pp. 1–6. [Google Scholar]
  38. Pacejka, H.B. Tire and Vehicle Dynamics, 2nd ed.; Elsevier: New York, NY, USA, 2006. [Google Scholar]
  39. Ferreau, H.J.; Kirches, C.; Potschka, A.; Bock, H.G.; Diehl, M. qpOASES: A parametric active-set algorithm for quadratic programming. Math. Program. Comput. 2014, 6, 327–363. [Google Scholar] [CrossRef]
  40. Tavan, N.; Tavan, M.; Hosseini, R.; Tavan, N.; Tavan, M.; Hosseini, R. An optimal integrated longitudinal and lateral dynamic controller development for vehicle path tracking. Lat. Am. J. Solids Struct. 2014, 12, 1006–1023. [Google Scholar] [CrossRef]
Figure 1. Structure of a vehicle system including the proposed algorithm.
Figure 1. Structure of a vehicle system including the proposed algorithm.
Algorithms 12 00044 g001
Figure 2. Structure of the friction coefficient estimation algorithm.
Figure 2. Structure of the friction coefficient estimation algorithm.
Algorithms 12 00044 g002
Figure 3. Vehicle structure model used for estimation.
Figure 3. Vehicle structure model used for estimation.
Algorithms 12 00044 g003
Figure 4. Structure of the Magic tire model.
Figure 4. Structure of the Magic tire model.
Algorithms 12 00044 g004
Figure 5. Structure of the Magic x function.
Figure 5. Structure of the Magic x function.
Algorithms 12 00044 g005
Figure 6. Structure of the Magic y function.
Figure 6. Structure of the Magic y function.
Algorithms 12 00044 g006
Figure 7. Structure of the Magic x , i function.
Figure 7. Structure of the Magic x , i function.
Algorithms 12 00044 g007
Figure 8. Structure of the Magic y , i function.
Figure 8. Structure of the Magic y , i function.
Algorithms 12 00044 g008
Figure 9. Tire force use ratio as a function of the slip ratio.
Figure 9. Tire force use ratio as a function of the slip ratio.
Algorithms 12 00044 g009
Figure 10. Tire force use ratio as a function of the side-slip angle.
Figure 10. Tire force use ratio as a function of the side-slip angle.
Algorithms 12 00044 g010
Figure 11. Effectiveness factor as a function of the tire slip ratio for different friction coefficients.
Figure 11. Effectiveness factor as a function of the tire slip ratio for different friction coefficients.
Algorithms 12 00044 g011
Figure 12. Principles of torque injection.
Figure 12. Principles of torque injection.
Algorithms 12 00044 g012
Figure 13. Procedure for torque injection.
Figure 13. Procedure for torque injection.
Algorithms 12 00044 g013
Figure 14. Sawtooth torque injection signal.
Figure 14. Sawtooth torque injection signal.
Algorithms 12 00044 g014
Figure 15. An example path.
Figure 15. An example path.
Algorithms 12 00044 g015
Figure 16. Division of the example path.
Figure 16. Division of the example path.
Algorithms 12 00044 g016
Figure 17. Speeds at nodes.
Figure 17. Speeds at nodes.
Algorithms 12 00044 g017
Figure 18. Speed errors at nodes.
Figure 18. Speed errors at nodes.
Algorithms 12 00044 g018
Figure 19. A path segment.
Figure 19. A path segment.
Algorithms 12 00044 g019
Figure 20. Road path.
Figure 20. Road path.
Algorithms 12 00044 g020
Figure 21. Curvature and friction coefficient of the road.
Figure 21. Curvature and friction coefficient of the road.
Algorithms 12 00044 g021
Figure 22. Lane-keeping performance of vehicles controlled by different controllers: (a) Controller A and (b) Controller B.
Figure 22. Lane-keeping performance of vehicles controlled by different controllers: (a) Controller A and (b) Controller B.
Algorithms 12 00044 g022
Figure 23. Speeds of vehicles controlled by different controllers: (a) Controller A and (b) Controller B.
Figure 23. Speeds of vehicles controlled by different controllers: (a) Controller A and (b) Controller B.
Algorithms 12 00044 g023
Figure 24. Path tracking errors: (a) Controller A and (b) Controller B.
Figure 24. Path tracking errors: (a) Controller A and (b) Controller B.
Algorithms 12 00044 g024
Figure 25. Steering angles of vehicles controlled by different controllers: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Figure 25. Steering angles of vehicles controlled by different controllers: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Algorithms 12 00044 g025
Figure 26. Drive torques of vehicles controlled by different controllers: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Figure 26. Drive torques of vehicles controlled by different controllers: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Algorithms 12 00044 g026
Figure 27. Injected drive torques: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Figure 27. Injected drive torques: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Algorithms 12 00044 g027
Figure 28. Estimated friction coefficients: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Figure 28. Estimated friction coefficients: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Algorithms 12 00044 g028
Figure 29. Speed plans: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Figure 29. Speed plans: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Algorithms 12 00044 g029
Figure 30. Piecewise speed plans: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Figure 30. Piecewise speed plans: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Algorithms 12 00044 g030
Figure 31. Speed Plan A: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Figure 31. Speed Plan A: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Algorithms 12 00044 g031
Figure 32. Speed Plan B: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Figure 32. Speed Plan B: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Algorithms 12 00044 g032
Figure 33. Speed Plan C: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Figure 33. Speed Plan C: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Algorithms 12 00044 g033
Figure 34. Speed Plan D: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Figure 34. Speed Plan D: (a) Controller A + ISPFCE and (b) Controller B + ISPFCE.
Algorithms 12 00044 g034
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParameterValueParameterValue
m1412 kg R g 0.3 m
l f 1.016 m l r 1.564 m
l s 0.77 m h c 0.54 m

Share and Cite

MDPI and ACS Style

Bian, C.; Zhu, T.; Yin, G.; Xu, L. Integrated Speed Planning and Friction Coefficient Estimation Algorithm for Intelligent Electric Vehicles. Algorithms 2019, 12, 44. https://doi.org/10.3390/a12020044

AMA Style

Bian C, Zhu T, Yin G, Xu L. Integrated Speed Planning and Friction Coefficient Estimation Algorithm for Intelligent Electric Vehicles. Algorithms. 2019; 12(2):44. https://doi.org/10.3390/a12020044

Chicago/Turabian Style

Bian, Chentong, Tong Zhu, Guodong Yin, and Liwei Xu. 2019. "Integrated Speed Planning and Friction Coefficient Estimation Algorithm for Intelligent Electric Vehicles" Algorithms 12, no. 2: 44. https://doi.org/10.3390/a12020044

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