Numerical Simulation Analysis of an Oversteer In-Wheel Small Electric Vehicle Integrated with Four-Wheel Drive and Independent Steering

Similar to conventional vehicle, most in-wheel small EVs that exist today are designed with understeer (US) characteristic.They are safer on the road but possess poor cornering performance.With recent in-wheelmotor and steer-bywire technology, high cornering performance vehicle does not limit to sport or racing cars. We believe that oversteer (OS) design approach for in-wheel small EV can increase the steering performance of the vehicle. However, one disadvantage is that OS vehicle has a stability limit velocity. In this paper, we proposed a Four-Wheel Drive and Independent Steering (4WDIS) for in-wheel small EV with OS characteristic.The aim of implementing 4WDIS is to develop a high steer controllability and stability of the EV at any velocity. This paper analyses the performance of OS in-wheel small EV with 4WDIS by using numerical simulation. Two cornering conditions were simulated which are (1) steady-state cornering at below critical velocity and (2) steady-state cornering over critical velocity. The objective of the simulation is to understand the behavior of OS in-wheel small EV and the advantages of implementing the 4WDIS. The results show that an in-wheel small EV can achieve high cornering performance at low speed while maintaining stability at high speed.


Introduction
In recent years, concerns over environmental issues by petroleum-based transportations along with the issue of fossil fuel depletion around the world have led to renewed interest over electric vehicle (EV).Small EVs have become prominent as a means of transportation, especially in urban areas where there are limited spaces [1,2].The driving method of such an electric vehicle can be divided into two types.The most basic type is the central motor system, where a conventional combustion engine is replaced with an electric motor.The other type is the in-wheel motor system, where a driving motor is located near the hub of each wheel.We believe that the inwheel motor system has tremendous potential in the future of EVs for its compact size and high-energy efficiency.In past research, results show that the in-wheel motor system offers more freedom of movement for each wheel and can be governed separately by control systems such as the antilock brake system (ABS) and the traction control system (TCS) [3][4][5].
A few car manufacturers have marketed compact, cheap, and reliable in-wheel small electric vehicle that is aimed for daily usage.Similar to conventional vehicle, most in-wheel small EVs that exist today are designed with understeer (US) characteristic [6].Even though they are safer on the road, they have poor cornering performance.With in-wheel motor and steer-by-wire technology nowadays, high cornering performance vehicle does not limit only sport or racing cars.We believe that oversteer (OS) design approach for in-wheel small EV can increase the steering response and lower the steering input with shorter cornering radius.However, one disadvantage of OS characteristic is that it has a stability limit velocity [7].In this paper, we proposed a Four-Wheel Drive and Independent Steering (4WDIS) for in-wheel small EV that is designed with OS characteristic.The aim of implementing 4WDIS is to develop a high steer controllability and stability of an OS in-wheel small EV at any velocity.
This paper analyses the performance of 4WDIS on an OS in-wheel small EV by using numerical simulation.Two cornering conditions were simulated which are (1) steadystate cornering at below critical velocity and (2) steady-state cornering over critical velocity.The objective of the simulation is to understand the behavior of an OS in-wheel small EV and the advantages of implementing the 4WDIS.The 4WDIS is governed by an Intelligence Steering Control System (ISCS) that evaluates the dynamics of the EV, the surrounding condition (i.e., narrow space), the road surface (i.e., dry, wet and icy road), and vehicle stability (i.e., crosswind and soft collision) in response to the drivers driving input.The ISCS will control the four-wheel driving and braking with each in-wheel motor and steer angle by the steering actuators to ensure safe maneuverability.The results show that, by implementing 4WDIS, an OS in-wheel small EV can achieve high cornering performance and maintain stability at any speed.

Main Symbols
Nomenclature shows the constant parameters and variables with the meaning, value, and unit used in this paper, respectively.

Limitations of OS and US Characteristics
A basic understanding of OS and US characteristic is necessary before any control method can be implemented.Figure 1 shows the relationship of the steady-state yaw rotational speed to velocity, and Figure 2 shows the relationship of steady-state side-slip angle to velocity between OS and US characteristic during steady-state cornering at constant steer angle.Based on the figure, the OS characteristic is separated into stable region and unstable region depending on the velocity of the vehicle.When the vehicle velocity  is below the stability limit velocity   , OS vehicle can achieve steady-state cornering with a higher yaw rotational speed in comparison to US characteristic.A higher yaw rotational speed should also mean that the vehicle can achieve a shorter turning radius.The side-slip angle becomes negative as the absolute magnitude increases with velocity.This means that when the velocity increases, the vehicle will point to the inner side of the circular path of the cornering.However, in the case of the velocity  ≥   , the yaw rotational speed and the side-slip angle will produce an infinite value with regard to a constant steer angle.It does not mean that the vehicle cannot be driven above this velocity but it is important to note that the vehicle motion is unstable.For this particular reason, most vehicles are manufactured with US characteristic that has no stability limit velocity as a compensation to low steering response.Nonetheless, US vehicle also has a tendency to cause accidents at high speed due to its restrictive steering.

Four-Wheel Drive and Independent Steering (4WDIS) Approach
4.1.Steering Characteristics.Two-wheel steering (2WS) or front wheel steering system (FWS) is generally known in all moving vehicles.This system works by rotating the front wheel in an angle that produces yaw motion, which then push the rear wheels with no steer angle sideways.A 2WS vehicle with oversteer characteristic will have a shorter turning radius but it is difficult to stabilize with front wheel steering alone.Four-wheel steering (4WS) has been introduced and marketed by Honda since 1987 in their Prelude model and then was joined by other car manufacturers to promote stability for conventional vehicle.The trend was short-lived due to the mechanical constrain of the system [8].Nonetheless, with inwheel motor system and steer-by-wire technology, a small electric vehicle can have the freedom of design and layout to implement Four-Wheel Drive and Independent Steering (4WDIS).In this research, the main purpose of implementing 4WDIS is to have high cornering performance for OS inwheel small EV while maintaining stability at any speed [9].Considering that each wheel's driving torque can be governed separately and the direction of the wheels can be manipulated with certain control, we had concluded that the steering modes for 4WDIS as in Figure 3 and the descriptions are as follows: (a) Opposite Steering (Phase 1).The rear wheels will turn on the opposite direction as the front wheels.This will increase the vehicle yaw rotational speed, which tighten the turning radius.
(b) Parallel Steering (Phase 2).The rear wheels will rotate in the same direction parallel to the front.This steering does not generate yaw rotational speed but instead reduces the yaw rotational speed.
(c) "Zero-Radius" Steering (Phase 3).The front wheels have to be "toe-in" position and the rear wheels in "toe-out" position.The driving direction of the left and right is in opposite.This allows the vehicle to rotate at its yaw rotational center.

Intelligence Steering Control System (ISCS
).An Intelligence Steering Control System (ISCS) is used to govern the 4WDIS.The ISCS evaluates the dynamics of the EV, the surrounding condition (i.e., narrow space), the road surface (i.e., dry, wet, and icy road), and vehicle stability (i.e., crosswind and soft collision) with regard to the drivers steering input.After the evaluation, the ISCS will choose the suitable steering characteristics to ensure stability and maneuverability.
In this paper, the evaluated parameter focuses on the dynamic of the in-wheel small EV and dry road surface.Thus, when studying the cornering performance of OS in-wheel small EV, only the alternation between phases 1 and 2 as shown in Figure 3 can be observed in the results.

Analysis Vehicle Model
Toyota COMS AK10E-PC shown in Figure 4, which is an inwheel small electric vehicle manufactured by Toyota Auto Body Co., was used as an analysis vehicle model in the simulation.The vehicle uses a 2WD in-wheel motor and has a front wheel steering.However, in the simulation, we added two in-wheel motors to the front wheels and two-steering actuator at the rear wheels.Thus, the mass of the vehicle model is increased from the actual vehicle.Furthermore, the weight distribution of the vehicle model was modified so that the weight shifted to the rear to represents a highly OS inwheel small EV while the dimensions of the vehicle remain the same.

Simulation Block Diagram
Figure 5 shows the simulation block diagram for the vehicle model and the ISCS.The vehicle model, which is an OS in-wheel small EV equipped with 4WDIS, is calculated by a nonlinear system.Much of past researches adopted easy way of modelling a control system based on a vehicle with linear characteristics.By this assumption, other parameters that could influence the dynamic motion of the vehicle are neglected and the modelled control system is deemed to be low precision [10][11][12][13].In this paper, the purpose of adopting nonlinear system for calculating the vehicle model is that the system represents more precisely an actual vehicle.However, a feedback control cannot be executed directly to nonlinear system because of yaw rotational speed coupling with the velocity components as expressed in dynamic equation of motion in later section.As an option, the ISCS, which is a state observer unit that consist of a feedback controlled linear system, is used to control the vehicle model.Each model in the simulation block will be explained in detail in the following sections.

Nonlinear Dynamic Equations.
The vehicle model has six degrees of freedom: three translational along the -, -, and -axes and three rotational about the -, -, and axes.In this study, however, a planar motion is assumed so that a translation along the -axis and rotations about the and -axes are disregarded.Figure 6 depicts the force vector used to construct the dynamic equations of motion for the vehicle and for the yaw rotational speed.Unlike linear dynamic equation, the vehicle side-slip angle at the center of gravity and wheels steer angle is not considered as small.Thus, we took the consideration that the longitudinal and lateral acceleration are influenced by the yaw rotational speed, which can be seen in (1).
The assumptions mentioned above necessitate calculating the dynamic equations for longitudinal and lateral velocities, including the coupling of the yaw rotation speed with the velocities: (1) 7.2.Nonlinear Tire Characteristic.In the case of linear model, the forces acting on the vehicle is determined by the motions of the vehicle that are mainly side-slip angle and yaw rotational speed at the center of gravity.In our study, the acting forces are influenced by the wheels of the vehicle.
Based on Fiala's theory, each wheel of the vehicle are modelled as a brush tire [14].This theory assumes that with side-slip angle, the tire tread base deforms elastically in the lateral and longitudinal direction of the wheel.Due to this, the equations for friction and lateral forces are dependent on the value of wheels side-slip angle, slip ration, and friction coefficient.The deformation of the tire tread rubber is shown by the dimensionless variable   : where When   > 0, the contact surface comprises adhesive and slip region.Thus, longitudinal force and lateral force can be written as However, when   ≤ 0, the contact surface only holds the slip region.Therefore, the longitudinal force and lateral force can be written as Here, the integral  represents front and rear, while the integral  represents right and left.The side-slip angle   which is mainly   ,   ,   ,   is given as follows: The slip ratio   can be defined as the traction between the road and the tire surface and is written as Then, by using the slip ratio value, the coefficient of friction   can be approximated by the following equation: where  = 1.0 (dry asphalt) ,  = 0.2 (icy road) .

State Observer Unit
8.1.Linear Dynamic Equations.The ISCS consist of a state observer unit.This unit uses linear dynamic equations for the translational motion and yaw rotational motion such as It is to be noted that there is no difference in the characteristics of the left and right tires due to linearity; hence, a 2-wheel model is derived from the equations.The sideslip angle  * is assumed to be small, which then lead to the assumption that the direction perpendicular to the traveling direction of the vehicle almost coincides with the lateral International Journal of Vehicular Technology direction .Therefore, no coupling exists between the yaw rotational speed and the vehicle's translational velocity.
The linear lateral force is defined by Here, the cornering stiffness of the front and rear wheels is   and   .These values determine the sensitivity of the ISCS [15].In contrast, with the tire characteristics of the nonlinear model, the actual vehicle motion determines the linear model lateral forces, which are related to the sideslip angle, yaw rotational speed, and front and rear steering angles.The above equations can be substituted into the linear dynamic equations of motion (10) and the yaw rotational speed (11).

Linearized Differential Equation.
The linearized differential equation represents the linear vehicle dynamic equation of motion for the state observer.It is arranged as a set of firstorder ordinary differential equations in the vector state form: where The state vector  represents the side-slip angle and yaw rotational speed of the vehicle.The steering angles of the front and rear wheels are the inputs of this system.The value for the front wheel angle   is a driver-selected input whereas the rear wheel angle   is initially 0 degrees.During high-speed cornering, the rear wheel steering is used as the control input by introducing a feedback gain −  : 8.3.Optimal Control.In order to find the optimal control for the linearized differential equation, the following evaluation function  is employed: where in which  11 ,  22 , and  are appropriately chosen constant weighting matrices for the state vector , the first-order differential of state vector ẋ , and the control input   , respectively.Table 1 shows the characteristics of the weighting matrices value that control the state and control value.The control input has its own high gain  and thus the weighting matric  = 1.The weighting matrices of  11 and  22 were used to stabilized the state vector, which is yaw rotational speed and side-slip angle to prevent it from reaching infinite value by the oversteering characteristic of the vehicle.The control input weighting matric  = 1.
The optimal solution for  can be designed if  is a positive definite matrix and the state observer control input   is given by where With  =   ≥ 0 being the unique positive-semidefinite solution of the algebraic Riccati equation, International Journal of Vehicular Technology 7 The unique positive-semidefinite  can also be expressed as Thus, (20) can be expanded as Finally, the output of the vehicle model and the state observer was compared to form a close loop system.In general, there is no other method to directly measure the output of side-slip angle.However, yaw rotational speed of a vehicle can be calculated with gyro sensors.In this system, the estimated yaw rotational speed output of state observer unit and the measured yaw rotational speed output of the vehicle model are compared and multiplied by high gain .Then, this gain is fed back into the linear model.The linearized differential equation of ( 13) becomes The renewed state observer control input can also be redefined as follows: With this equation, the FFC for the 4WDIS in-wheel small EV can be changed to a FBC that gives greater precision.

Determine the Vehicle Model Steering Characteristic.
In the previous section, we mentioned that the vehicle model used the Toyota COMS specification and modified it to represent an OS in-wheel small electric vehicle in the simulation.
A simulation of a steady-state cornering at various velocities with regard to a constant front steer angle was done to determine the stability limit velocity of the vehicle model before 4WDIS could be implemented.The constant front steer angle is fix at 10 ∘ at  = 10 s.When the yaw rotational speed and the side-slip angle give an infinite value, the velocity before that moment is acknowledged as the stability limit velocity.

Implementation of 4WDIS with ISCS.
A 4WDIS for OS in-wheel small EV was constructed as shown is previous simulation block diagram.In order to investigate the performance of 4WDIS in depth, we have divided the simulation as follows: (1) Steady-state cornering below stability limit velocity.
Similar to previous simulation condition, a steady-state cornering of a 4WDIS in-wheel small EV at below stability limit velocity and over stability limit velocity, with regard to a constant front steer angle, was carried out.The constant front steer angle is fix at 10 ∘ at  = 10 s.The results of the steadystate yaw rotational speed and steady-state side-slip angle for 4WDIS are compared to the previous simulation.

Stability Velocity Limit of the OS Vehicle Model during 2WS Cornering.
A simulation of steady-state cornering was carried out to determine the stability limit velocity of the OS vehicle model.Figures 7 and 8 show the results of the steadystate yaw rotational speed and steady-state side-slip angle with regard to constant front wheel steer angle, respectively.Based on both of these graphs, in regard to constant front steer angle of 10 ∘ , the stability limit velocity   is 15.5 km/h.At this velocity, the vehicle is stable and can reach the maximum steady-state yaw rotational speed of 1.01 rad/s with a maximum absolute magnitude of the steady-state side-slip angle of 0.301 rad.When the velocity is above 15.5 km/h, the yaw rotational speed and the side-slip angle will give an infinite value.However, the vehicle can still be driven but it is very unstable.This shows that we have successfully modelled a highly OS in-wheel small EV.
Nonetheless, we believe that even in the stable region which is below the stability limit velocity, the OS small EV can be improved.In Figure 7, when the constant velocity International Journal of Vehicular Technology increased, the vehicle model took longer time to generate a constant yaw rotational speed.This is due to the increase in side-slip angle to the constant velocity.The relation between steady-state yaw rotational speed and velocity is plotted in Figure 9, while the relation between steady-state side-slip angle and velocity is plotted in Figure 10.

Steady-State Cornering below Stability Limit Velocity with 4WDIS.
Theoretically, OS vehicles can already produce high yaw rotational speed and stability when cornering below the stability limit velocity.Logically, any control system is unnecessary in the stable region.We can assume that the cornering performance in the stable region of the OS in-wheel small EV will be suppressed by implementing the 4WDIS.In this section, the control effect of 4WDIS implementation onto an OS in-wheel small EV during steady-state cornering below the stability limit velocity in regard to constant front steer angle 10 ∘ is examined.Figures 11 and  12 show the results of the steady-state yaw rotational speed and steady-state side-slip angle for the vehicle model with 4WDIS, respectively.The average time for the yaw rotational speed to reach a constant value was also reduced to 25 s in comparison to the 2WS vehicle model.Even with high yaw rotational speed, the steady-state side-slip angle also remains constant.The steady-state side-slip angle was in negative values which indicate that the vehicle was facing the inner circle of the cornering direction as the cornering initiated.For better comparison, the relation between steady-state yaw rotational speed to velocity and steady-state side-slip angle to velocity for vehicle model with 4WDIS and 2WS was plotted in Figures 13 and 14 increased from 10 km/h to 15 km/h, the steady-state yaw rotational speed gradually decreased.A similar pattern occurred for the side-slip angle if we observe in absolute magnitude.
In order to comprehend these matters, we have to examine the output of the ISCS. Figure 15 shows the rear wheel steering angle during the steady-state cornering of the 4WDIS vehicle model.As soon as the front wheel steer input initiated at  = 10 s, the rear wheel steer initiated at the opposite direction as the front wheel.The rear wheel steering angle maintained an opposite angle as the front wheel which increased the yaw rotational speed.As the velocity approach near the stability limit, the opposite angle of the rear wheel decreased.

Steady-State Cornering over Stability
Limit Velocity with 4WDIS.The stability limit velocity of the modelled OS in-wheel small EV has been determined in the previous simulation.The previous simulation also proves that the implementation of 4WDIS is able to improve the maneuverability of the OS small electric vehicle even at stable region.The more anticipated matter is the control of 4WDIS for the OS small electric vehicle during cornering above the stability velocity limit.The previous steady-state cornering simulation condition is repeated again for this vehicle model.
In this section, the control effect of the 4WDIS on an OS in-wheel small EV during steady-state cornering over the stability limit velocity in regard to constant front steer angle 10 ∘ is investigated.Figure 16 shows the result of the steady-state yaw rotational speed for the vehicle model with 4WDIS.From this figure, the yaw rotational speed decreased as the velocity of the vehicle model increased.However, if we observe closely, the time to produce a constant steadystate yaw rotational speed is instantaneous and maintains the average time of 25 s, especially at velocity range of 30 to 50 km/h.The results of steady-state side-slip angle for this simulation in Figure 17 show a reverse pattern from the yaw rotational speed.The side-slip angle increased as the velocity of the vehicle model increased.Based on these results, we can assumed that the 4WDIS control input suppressed the output more as the velocity of the vehicle model increased.
International Journal of Vehicular Technology In order to understand the control method of the 4WDIS by the ISCS, we have to examine the output of the system.Figure 18 shows the front and rear wheel steering angle during the steady-state cornering simulation of the 4WDIS vehicle model.In the first few second after the front wheel steering angle initiated, the ISCS produced a rear wheel steering angle in minus value that gave the instantaneous yaw rotational speed mentioned before.The minus angle of the rear wheel was smaller as the steady-state velocity increased because the excess yaw rotational speed produced by OS vehicle model also increased.As time progressed, the ISCS gradually increased the rear wheel steering angle to counter the excess yaw rotational speed and sustain a constant value.The angle of the rear wheel is higher for velocity range over 15 km/h which proves that our early assumption was correct.Based on these results, an OS in-wheel small EV can achieve a stable steady-state cornering even when the velocity is over the stability limit velocity with the implementation of 4WDIS.Simply put, the 4WDIS eliminate the stability velocity that was present in the OS small electric vehicle with 2WS.The 4WDIS control input suppressed the yaw rotational speed and side-slip angle during high speed by rotating the rear wheels at the same direction as the front wheels.Additionally, it can also improve the maneuverability of the vehicle at low speed by rotating the rear wheels at the opposite direction of the front wheels.

Conclusion
There are various methods to design an OS vehicle.In this paper, we chose the transfer of weight distribution method to model an OS in-wheel small EV.There are also a few ways to recognize an OS vehicle and one of them is that the vehicle will have a stability limit velocity.When a steady-state cornering simulation was carried out below the stability limit velocity   , the OS in-wheel small EV can achieve steadystate cornering with a higher yaw rotational speed with regard to constant low steering input.However, when the vehicle velocity  ≥   , the yaw rotational speed and the side-slip angle produced an infinite value in regard to the same constant steer angle.
The main objective of implementing 4WDIS onto an OS in-wheel small EV is to ensure stability during cornering at any velocity.The 4WDIS utilized the opposite at low speed to gain yaw motion and parallel steering to reduce the yaw motion of the vehicle.
In this paper, the Intelligence Steering Control System (ISCS) evaluated the dynamic of the in-wheel small EV and dry road surface.In the future, other parameters such as lateral disturbance or icy road surface which has effect on the driving or the vehicle will be taken into consideration in the simulation analysis.Wheel rotational speed, rad/s.

Figure 1 :
Figure 1: Steady-state yaw rotational speed to velocity between OS and US characteristic.

Figure 4 :
Figure 4: Analysis vehicle model of Toyota COMS.

Figure 7 :
Figure 7: Steady-state yaw rotational speed to time change with 2WS.

Figure 15 :
Figure 15: The front and rear wheel steering angle during the steadystate cornering at below stability limit velocity for the vehicle model with 4WDIS.

Figure 17 :
Figure 16: Steady-state yaw rotational speed to time change for the vehicle model with 4WDIS.

Figure 18 :
Figure18: The front and rear wheel steering angle during the steadystate cornering at over stability limit velocity for the vehicle model with 4WDIS.

Figure 19 :
Figure 19: (a) Steady-state yaw rotational speed to velocity for the vehicle model with 4WDIS.(b) Steady-state side-slip angle to velocity for the vehicle model with 4WDIS.
Constant Parameters Symbols, Meaning, Value, and Unit : Width of wheel interact surface, 0.1, m   : Front tread, 0.840, m   : Rear tread, 0.815, m : Yaw inertia moment at gravity point of vehicle, 1470, kgm 2   : Longitudinal tread rubber stiffness, 3.33 × 10 6 , N/m 3   : Lateral tread rubber stiffness, 3.33 × 10 6 , N/m 3   : Front wheel cornering stiffness, 2000.00,N/rad:Rear wheel cornering stiffness, 2612.62,N/rad:Length of wheel interact surface, 0.15, m   : Length from front axle wheel to gravity point of vehicle, 0.725, m   : Length from rear wheel axle to gravity point of vehicle, 0.555, m : Length from front axle to rear axle of the wheel, 1.28, m : Mass of the vehicle, 421.61, kg : Radius of the wheel, 0.23, m.Velocity at center gravity, m/s   : Stability limit velocity, m/s   ,   ,   ,   : Friction force for each tire, N   ,   ,   ,   : Lateral force for each tire, N   : Wheels load subject to load transfer, N : Side-slip angle of vehicle, rad  * : Side-slip angle of state observer, rad   ,   ,   ,   : Side-slip angle of each wheel, rad :Friction coefficient   ,   :Front and rear wheels steer angle, rad : :L o n g i t u d i n a l v e l o c i t y , m / s V:Lateral velocity, m/s :