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

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.


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.

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", "rl", and "rr" denote the values corresponding to the front left, front right, rear left, and rear right tires, respectively.F xi and α i (i = f l, f r, rl, and rr) 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.

Tire Force and Effective Radius Estimation
The estimated vertical load force on each tire can be expressed as [27][28][29]: 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]: where i = f l, f r, rl, and rr; Rwi is the estimated effective radius of the corresponding tire; R g is the geometric radius of the tire; and d ri = Fzi 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: where i = f l, f r, rl, and rr; 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]: where Fy f and Fyr are the estimated lateral forces on the front and rear axles, respectively, and I z denotes the moment of inertia around the vertical axle.

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]: and: where v xi and θ i (i = f l, f r, rl, and rr) 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: and: where i = f l, f r, rl, and rr and ω i is the rolling speed of the corresponding wheel, which can be obtained from an angular speed sensor on that wheel.

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.

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: 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: 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 Figures 5 and 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: where i = f l, f r, rl, and rr 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 Figures 7 and 8, respectively.To facilitate the discussion of the estimation of the friction coefficient, the tire force use ratio is defined as follows: The relationships between the tire force use ratio and the slip ratio and side-slip angle with different friction coefficients are plotted in Figures 9 and 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: 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.

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: (14) and: To facilitate the calculation, the following optimization problem is defined to estimate the friction coefficient of each front tire: min : 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, µ eU is the threshold for accepting an estimated friction coefficient as reliable, µ eL 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.

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 Figures 12 and 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,sum denotes the resultant longitudinal drive force on the four tires; T inj 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 Figures 12 and 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 cyc 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: where t is the current time, T inj is the injected torque, and K T is a constant parameter used to calculate T inj .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 inj and T inj may be different from one cycle to another.The torque injected into the front wheels is T inj .The torque injected into the rear wheels is −T inj .

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.

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 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: 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 ) can be calculated as: 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.

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.

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: where R min 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: where f sS < 1 is a safety factor for avoiding skid conditions.To avoid rollover, the lateral acceleration should satisfy: The limit on the speed for avoiding rollover can thus be defined as: where f sR < 1 is a safety factor for avoiding rollover conditions.Thus, the maximum safe speed on a path segment can be expressed as:

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: where v smax,L and v smax,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: where f sL < 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: 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: 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: where v dec = t t f sL • µ • g.

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: and: By substituting Equation (19) into Equations ( 30) and (31), the speed-planning constraints can be rewritten as: and: Equations ( 32) and ( 33) can be further rewritten as: where . .

Objective Function
To reduce the error between the desired and planned speeds, the following objective function is defined: where . The optimal planned speed can be obtained by minimizing F(e).

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) (36) subject to: 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.

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: where v s and v e are the planned speeds at p s and p e , respectively.

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 Figures 20 and 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 Figures 22-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 Figures 22 and 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 Figures 25 and 26.The details of the friction estimation and speed planning are further plotted in Figures 27-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.

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.

Figure 1 .
Figure 1.Structure of a vehicle system including the proposed algorithm.

Figure 2 .
Figure 2. Structure of the friction coefficient estimation algorithm.

FFigure 4 .
Figure 4. Structure of the Magic tire model.

FFigure 5 .
Figure 5. Structure of the Magic x function.

FFigure 6 .
Figure 6.Structure of the Magic y function.

FFigure 7 .
Figure 7. Structure of the Magic x,i function.

FFigure 8 .
Figure 8. Structure of the Magic y,i function.

Figure 10 .
Figure 10.Tire force use ratio as a function of the side-slip angle.

9 Figure 11 .
Figure 11.Effectiveness factor as a function of the tire slip ratio for different friction coefficients.

Figure 22 .Figure 23 .
Figure 22.Lane-keeping performance of vehicles controlled by different controllers: (a) Controller A and (b) Controller B.

Figure 24 .
Figure 24.Path tracking errors: (a) Controller A and (b) Controller B.
Tire force use ratio as a function of the slip ratio.