Integrated Yaw Stability Control of Electric Vehicle Equipped with Front/Rear Steer-by-Wire Systems and Four In-Wheel Motors

: This paper presents the integrated motion control method for an electric vehicle (EV) equipped with a front/rear steer-by-wire (SbW) system and four in-wheel motor (IWM). The proposed integrated motion control method aims to maintain stable cornering. To maintain vehicle agility and stability, the lateral force and yaw rate commands of the vehicle are generated by referring to the neutral steering characteristics. The driver’s driving force command, the lateral force command based on the bicycle model, and the yaw moment generated by the high-level controller are distributed into the driving force of each wheel and the lateral force of the front and rear wheels by the yaw moment distribution. Finally, the distributed forces are directly controlled by a low-level controller. To directly control the forces, a driving force observer and a lateral force observer were introduced via driving force estimation in the IWMs and rack force estimation in the SbW system. The control performance is veriﬁed through computer simulations.


Literature Review
EV has attracted considerable attention owing to concerns regarding air pollution and resource depletion. EV's drivetrain has many advantages over an internal combustion engine. In addition, the SbW system has attracted interest in recent years. The SbW system has many advantages, as the steering wheel is disconnected from the rack-and-pinion system. These advantages can be summarized as follows [1]: • Improved configuration of a vehicle's steering functionality; • Improved design flexibility for a vehicle; • Reduced vibration, which enhances the driver's comfort; • Ability of the rack-and-pinion system to receive road surface information independently of the driver's operation.
Active motion control, which ensures the safety of passengers and vehicles, has been studied extensively in recent decades. Various vehicle motion control systems, such as four-wheel steering (4WS), active front steering (AFS), active rear steering (ARS), and direct yaw moment control (DYC), have been developed to maintain vehicle stability, even in emergencies. In addition, vehicle motion control systems for longitudinal motion, such as the anti-lock brake system (ABS) and traction control system (TCS), have also been studied.
In [2], two new methods were proposed on the basis of DOB and sliding mode control (SMC) for ABS control even in the case of uncertainties in the road conditions and the height of the center of gravity (CG). Thus, the wheel slip ratio could be maintained at the desired value. The vehicle braking performance was improved on the basis of iterative learning control theory by integrating the vehicle model, slip ratio model, and vehicle speed observer [3]. A TCS for the IWM system was proposed as a driving torque control method for EV on the basis of a driving torque observer. Driving torque control improves vehicle stability by preventing excessive wheel slip on the basis of the relationship between the driving torque and the slip ratio [4]. SMC-based wheel slip control was designed without a tire-road friction observer; tire-road friction is the most difficult to predict or measure among the various model uncertainties. The controller was verified using a vehicle equipped with a front IWM [5].
The additional yaw moment required for yaw rate control can be generated via active steering as well as on the basis of the wheel braking force and driving force. The purpose of DYC is to create an additional yaw moment to implement the yaw motion intended by the driver. A method for applying adaptive sliding mode control was presented to compensate for the nonlinear characteristics, disturbances, and model uncertainties of tires, and an experiment was performed to confirm the control performance [6]. In [7], to improve the maneuverability and stability of yaw motion during cornering, a robust DYC design based on the vehicle side slip angle was adopted by employing DOB. As EVs employ IWMs, they are particularly suitable for active torque distribution control because all the wheels can generate the torque independently [8]. In [9], the yaw moment control was designed by considering the saturation of the tire and applying variable cornering stiffness. In [10], smooth SMC was introduced for yaw rate control of target vehicles with electronic four-wheel drive (e-4WD) coexisting with IWMs in the front and a mechanical drivetrain in the rear, and torque distribution based on daisy-chaining allocation was proposed.
AFS, which can improve the handling characteristics of a vehicle while driving, has attracted considerable attention in the automotive industry. AFS is capable of generating additional lateral forces on the front wheels by controlling the steering angle. Most AFS control systems are based on the feedback of the yaw rate, side slip angle, and lateral acceleration. When the lateral acceleration is relatively small, the lateral force has a linear relationship with the wheel side slip angle and the 4WS and AFS can be controlled relatively ideally. When the lateral acceleration is large, the tire force is saturated and the 4WS and AFS cannot control the vehicle. The effects of various types of torque distribution on vehicle performance were analyzed under various steering and straight driving conditions, and it was confirmed that the strategy of torque distribution based on the driving conditions increases vehicle stability [11]. In [12], feedback control based on the yaw rate and lateral acceleration through an SbW system was proposed, and the AFS of the SbW system was shown to be more effective than DYC. Active steering control of an SbW system was proposed to adjust the handling characteristics of the vehicle according to the driver's preference, and its effect was verified experimentally [13]. In [14], because the lateral tire force measurements were obtained using multi-sensing hub units developed by NSK Ltd., direct lateral tire force control could be realized.
In [15], an integrated control strategy based on 4WS and DYC was proposed. An adaptive sliding mode control that is robust against model uncertainty and external disturbances was designed to calculate the desired yaw moment and front wheel steering angle [16]. The Takagi-Sugeno fuzzy control technique was introduced to solve problems related to model uncertainty and nonlinearity, and the lateral stability of the vehicle was improved by integrating AFS and DYC [17]. In [18], four critical conditions were introduced, and an integrated motion control plan consisting of three layers was presented. The top-level controller calculates the yaw moment of the vehicle, the middle-level force allocation controller distributes the braking force or driving force and steering angle of each wheel, and the bottom-level slip ratio controller generates the driving force or braking force of each wheel.
There are many difficulties in realizing the integrated control method, and first, high performance of the vehicle control unit is required. In [19], an asymptotic optimality approach with low memory requirements and low computational complexity was demonstrated on the basis of the Lyapunov function and cascade theory. Second, a robust control method against model uncertainty and external disturbances of the vehicle is required. In [20], a nonlinear disturbance observer based on second-order sliding mode control was proposed to overcome the chatter in SMC. In [21], backstepping control and second-order SMC were proposed to overcome the chatter in sliding mode control, nonlinear characteristics of the vehicle, model uncertainty, and external disturbances. As the vehicle motion is predominantly governed by the forces generated between the tires and the road, real-time knowledge of the tire forces is crucial when predicting vehicle motion; accordingly, several methods have been studied for real-time tire-road force estimation [22][23][24][25][26].

Contributions of the Article
This paper proposes advanced yaw rate control and lateral force estimation of EV equipped with 4WS of the SbW systems and 4WD of the IWM systems. Integrated motion control of EVs equipped with 4WS of the SbW system and 4WD of the IWM system can generate the independent driving force of four wheels and lateral force of front and rear wheels, and appropriate control input must be made for driver controllability and vehicle stability. The existing vehicle motion control method generates four independent driving forces and steering angles of the front and rear wheels. It is difficult to consider the vehicle's stability in this method. In addition to this, studies such as model predictive control methods are not only difficult to perform real-time control, but also difficult to implement vehicle integrated motion control algorithms such as vehicle side slip angle, roll angle, and pitch angle control. Therefore, the proposed integrated controller implements a 1 ms loop control algorithm by composing three layers of a high-level controller, a yaw moment distribution, and a low-level controller considering vehicle stability and handling characteristics. The contributions of the article are as follows. • In contrast to the commonly used lateral force estimation that relies on the tire model, in this study, the lateral tire force is estimated using the DOB of the SbW system and the kinematics of the suspension. • By minimizing the workload on each tire, the stability of the vehicle is improved. • The 3-DOF motion control of the vehicle is possible through the tire vector control, which directly controls the longitudinal and lateral forces of the tire.
The remainder of this paper is organized as follows. Section 2 introduces the nonlinear 3-DOF model of 4WS and the three layers for vehicle motion control. Section 3 describes a high-level controller for yaw stability control based on a sliding mode controller. In addition, yaw moment distribution for distributing the lateral force of the front/rear SbW system and the driving force of each wheel is introduced. Section 4 describes the low-level controller, driving force controller, and lateral force controller. Section 5 presents the evaluation results of the proposed controller obtained by a computer simulation. Finally, Section 6 concludes the paper.

Vehicle Modeling
In this section, a 3-DOF vehicle model is introduced to describe the motion of 4WS and IWM driven EV. In addition, the three layers for vehicle motion control are presented.

3-DOF Vehicle Model
The planar vehicle model shown in Figure 1 includes the longitudinal, lateral, and yaw motions, which are described as follows: where m is the total mass of the vehicle, I z is the yaw moment of inertia at the CG, γ is the yaw rate, and a x and a y are the longitudinal and lateral accelerations of the vehicle, respectively. Further, M z denotes the yaw moment generated by the driving force and the lateral force, d denotes the model uncertainty and disturbance. F x and F y are the longitudinal and lateral forces. The first subscripts of F x and F y represent the front and rear axles ( f , r), whereas the second subscripts represent the left and right (l, r). δ f and δ r are the front and rear steering angles at the wheel, i.e., δ f l = δ f r = δ f , δ rl = δ rr = δ r . l f and l r are the distances from the CG to the front and rear axles, and d f and d r are the track widths of the front and rear axles. The pitch, roll, and vertical motions are ignored in this model, and assuming that the difference between the left and right lateral forces is negligible, the yaw motion can be expressed as four driving forces and two lateral forces.

Three Layers Control System
In this study, the vehicle motion controller consists of three layers, and the overall algorithm of the vehicle motion controller is shown in Figure 2. The yaw rate command (γ d ) is from a reference model based on the driver's steering angle (δ) and vehicle speed (v x ), i.e., the bicycle model for steady-state cornering. In addition, the lateral acceleration command (a * y ) is designed as a reference model, and the lateral force command is generated by lateral dynamics, i.e., F * y = ma * y . The desired yaw rate and lateral acceleration for a given steering angle and vehicle speed are as follows: where τ r is the cutoff frequency of the desired model filter, l is the wheelbase, ρ rp is the rack-and-pinion gear ratio, and ρ wr is the rack-and-wheel gear ratio. k us is the vehicle stability factor that describes the steering characteristics. The vehicle characteristics are classified into understeer, oversteer, and neutral-steer according to the sign of l f C f − l r C r in k us .

High-Level Controller
The high-level controller, combined with a yaw stability controller and a yaw moment observer, is shown in Figure 3. The purpose of yaw stability control is to improve the handling performance of the vehicle and maintain its stability during cornering. The yaw stability controller is designed by employing the SMC method, which is known as a robust control method. To overcome chatter, a saturation function is employed instead of a sign function in the SMC method. Although the saturation function reduces the chatter, the tracking performance is degraded. A trade-off between chatter reduction and tracking performance can be achieved by employing a yaw moment observer (YMO) and adjusting the thickness of the saturation function.

Yaw Moment Observer
YMO is designed using the control input M z of the yaw motion and the measured yaw rate, and the nominal yaw model is given by: where I zn is the yaw moment of the nominal inertia at the CG and P n (s) is the transfer function of the nominal model. As shown in Figure 3, it is important to set an appropriate Q-filter cutoff frequency because the performance of YMO depends on it. The Q-filter is designed as a low-pass filter satisfying the robust stability condition. In general, the Q-filter in YMO is designed as a first-order low-pass filter as follows: where ω Q is the cutoff frequency.

Design of Yaw Stability Controller
The yaw stability controller is designed in the SMC method and consists of two steps. First, the sliding surface is designed so that the system trajectory obtains the desired characteristics along the defined surface. Then, a discontinuous controller is designed so that the system trajectory reaches the sliding surface within a finite time. Sliding mode control has a fast response, excellent transient performance, and robustness against disturbances and model uncertainty. The sliding surface (S(t)) and reaching condition (Ṡ(t)) are given by: where k y 1 , k y 2 , and k y 3 are the controller gain and e = γ d − γ. The saturation function is defined as sat . The control input of the proposed SMC is given by the following equation: Using M z , the closed-loop error dynamics can be expressed as follows: Theorem 1. Suppose that d of yaw motion is bounded, i.e., |d| < |d + |, and k y 1 , k y 2 , and k y 3 are selected according to: If the above-mentioned conditions are satisfied, the system is uniformly stable.
Proof. The following positive Lyapunov function V is used to prove that the closed-loop control system is uniformly stable: The time derivative of V is given by: Because the discontinuous saturation function is used, it is divided into the range when S ≥ φ and S < φ. When S ≥ φ,V is expressed as: When S < φ,V is expressed as: Defining k y 3 > D, this system is stable at φD φk y 2 +k y 2 < S. Thus, this system is uniformly stable. YMO improves the control performance of the system by reducing the magnitude of D.

Yaw Moment Distribution
In this study, the yaw moment distribution is considered as the workload of the tire. The forces acting on the tire, i.e., the driving force and lateral force, must be within the friction circle to stably drive the vehicle [15]. For this reason, the yaw moment distribution designs a cost function that minimizes the workload of the tire, which is given by: where η is the tire workload and µ is the coefficient of friction between the tire and the road surface. The vector x is denoted by a diagonal matrix, expressed as µ ij and F zij . The lateral force is estimated as the sum of the left and right wheels, and the lateral force of each wheel is assumed as follows: The vertical tire force F zij for weight transfer owing to the longitudinal and lateral acceleration of the vehicle is given by: where g is the gravitational acceleration and h is the height of the CG. The weight transfer is the complexity of the suspension dynamics and disturbance; therefore, only the longitudinal and lateral accelerations are considered. The constraints are as follows: where the coefficient matrix is denoted by A and b is the command of the vehicle motion control. For the yaw moment distribution where the cost function J has the minimum value, the solution x * is as follows. For further details, refer to Appendix A.
The distributed driving force and lateral force on each wheel of the yaw moment distribution are controlled by the driving force controller and lateral force controller of the low-level controller.

Low-Level Controller
This section describes a driving force controller and a lateral force controller that provide feedback of the driving force and lateral force, as shown in

Driving Force Observer
The driving force observer is designed to estimate the friction force between the tire and the road. Assuming that the current control of the IWM is sufficiently fast and accurate, the driving force can be estimated using the nominal model and angular velocity of the wheel. The wheel moment equation of the quarter-car model is given by: where J t is the total inertia of the wheel, ω w is the wheel angular velocity, g pg is the gear ratio, T m is the motor torque, F x is the driving force, and R is the wheel radius. The subscript i is the front and rear wheels and j is the left and right wheels. The state observer is given by: where ω m is the angular velocity of the IWM.

Driving Force Controller
The driving force controller based on the estimated force feedback is shown in Figure 4. The driving force controller is designed with a multiple-SMC method. The driving force controller consists of a DFO and a wheel speed control. The slip ratio is generated by controlling the driving force. The relationship between the wheel speed and the slip ratio 7 is given by: where y is the slip ratio when the vehicle decelerates, and the slip ratio and y are approximately equal when the vehicle accelerates in the vicinity of slip ratio = 0. The first-order dynamics of the driving force is expressed as follows: where τ is the relaxation time constant, and D sij is the driving stiffness of each wheel. The sliding surface (S 1 ) and reaching condition (Ṡ 1 ) of the outer loop is given by: where e = F * xij − F xij . The control input y ij is given by: where D s is the nominal driving stiffness. Using y ij , the closed-loop error dynamics can be written as follows: Theorem 2. Suppose that d 1 is bounded, i.e., |d 1 | < |d + 1 |, and the gains k 1 and k 2 are selected according to: If the above-mentioned conditions are satisfied, the system is uniformly stable.
Proof. The following positive Lyapunov function V 1 is used to prove that the closed-loop control system is uniformly stable: The time derivative of V 1 is given by: Defining k 2 > D 1 , this system is stable at φ 1 D 1 k 2 < S 1 . Thus, this system is uniformly stable. The dynamics of wheel speed control is expressed as follows: where T d is the driving torque. The sliding surface (S 2 ) and reaching condition (Ṡ 2 ) of the inner loop is given by: where e = ω * mij − ω mij . The control input T mij is given by: where J n is the nominal inertia. Using T mij the closed-loop error dynamics can be written as follows:Ṡ Theorem 3. Suppose that d 2 is bounded, i.e., |d 2 | < |d + 2 |, and the gains k 3 and k 4 are selected according to: If the above-mentioned conditions are satisfied, the system is uniformly stable.
Proof. The following positive Lyapunov function V 2 is used to prove that the closed-loop control system is uniformly stable: The time derivative of V 2 is given by: This system is stable at J ti D 2 J ni k 4 < S 2 . Thus, this system is also uniformly stable.

Lateral Force Observer
In this study, the rack and pinion are installed in the front and rear steering systems, and the motor is connected to the pinion gear. The lateral force, which is the sum of the left and right forces, can be estimated from the rack force acting on the pinion gear and the kinematics of the suspension. DOB is employed for estimating the rack force, and assuming that the variation of the ratio of the transport distance of the rack gear to the wheel angle is negligible, i.e.,ρ wr ≈ 0, the kingpin torque can be estimated using the rack force. The kingpin torque is generated by the driving force, vertical force, and lateral force, and the estimates of the vertical force and driving force are presented in Sections 3.3 and 4.1, respectively. The dynamics of the pinion gear is given by: where J r and B r are the steering system inertia and friction, respectively. Further, θ, T p , and T dis are the rotation angle, motor torque, and disturbance of the pinion gear. The subscripts i is the front and rear wheels. The kingpin torque observer is shown in Figure 7. The estimated kingpin torque is given by: whereT dis is the estimated disturbance. J n and B n are the nominal inertia and friction of the steering system, respectively. The Q-filter is designed as a second-order low-pass filter, which is given by: where ω q and ζ are the cutoff frequency and damping ratio, respectively. The kingpin torque generated by the driving force, lateral force, and vertical force is analyzed by focusing on the geometric parameters. The steering geometry is defined in Figure 8. The kingpin axis vector of the front left steering geometry is obtained using the kingpin angle (Ψ) and caster angle (Φ) of the two-dimensional plane in Figure 8, and the kingpin axis vector ( a k ) is given by: where L x and L y are the kingpin offset and caster trail at the wheel center, respectively. The direction of the kingpin axis vector ( u) is given by: The vector ( r) perpendicular to the kingpin axis vector and passing through the force application point is given by: where (r x , r y , r z ), is the point of application of the force at the tire and l, m, n is the direction vector of r. The intersection (p) of r and a k is as follows. For further details, refer to Appendix B and Figure 8.
The moment arm (R k ) between the force acting on the tire and the kingpin axis vector ( a k ) is given by: The kingpin torque can be calculated using the driving force, lateral force, and vertical force acting on the tire and the moment arm and is expressed as: The lateral force acting on the left and right tire is given by: The sum of the lateral forces of the left and right wheels is given by: where ∆F zi and ∆F xi are the differences between the left and right vertical and driving forces, respectively.

Lateral Force Controller
The lateral force controller, which directly controls the lateral force, is shown in Figure 6. It consists of a lateral force observer and pinion motor position control that controls the pinion gear angle of the rack and pinion system. The lateral force control is designed on the basis of a linear tire model and first-order dynamics. The lateral force of the linear tire model uses the vehicle velocity, yaw rate, side slip angle (β), and steering angle. In this study, the position of the pinion gear is controlled and the angle of the pinion gear is regarded as the only control input to control the lateral force. The linear lateral tire model with for first-order dynamics is given by: where τ y is the nominal relaxation time constant, C is the nominal cornering stiffness, and α is the side slip angle at the wheel. The dynamics of the lateral force are given by: The subscript i is the front and rear wheels. The lateral force generated by the side slip angle and yaw rate is regarded as the disturbance because the side slip angle, which affects the lateral force, is not measurable owing to the sensor cost and the yaw rate is strongly affected by the time-varying parameter. Thus, the dynamics of the lateral force can be rewritten as:Ḟ where d yi includes external disturbances and the effects of time-varying parameters, the side slip angle and yaw rate. The nominal model of the lateral force control is given by: The lateral force control system consists of a feedforward controller, feedback controller, and DOB. As stated previously, d yi plays a crucial role in the lateral force control. DOB is employed to compensate for d yi . The feedforward controller (C l f f ) in the lateral force control system is designed as the inverse of the nominal model.
The feedback controller (C l f b ) is designed with a pole-zero cancellation method based on the nominal model, and it is expressed as: where k P and k I are the proportional and integral gains, respectively. The pinion motor position control consists of a feedforward controller and feedback controller. The feedforward controller is designed on the basis of a lateral force command, and the feedback controller is designed using a pole-zero cancellation method based on the nominal model of the steering system. The relationship between the feedforward controller and the lateral force is given by: where r kp is the ratio of the lateral force and the force acting on the pinion gear, which is set to 78 in this study. Because the torque of the pinion motor is governed by the lateral force generated during driving, the feedforward controller is designed as the kinematics between the lateral force command and the pinion motor, assuming that the difference between the actual lateral force and the estimated lateral force is negligible. The nominal model of the steering system is given by: where J rn and B rn is the inertia and friction coefficient of the nominal model, respectively, and T pm is the torque of the pinion motor. The feedback controller for pinion motor position control is given by: where k pP and k pD are the conventional proportional and differential gains of pinion motor position control, respectively.

Simulation Environment
The computer simulations were performed using the vehicle simulation software CarSim, MATLAB/Simulink, and Amesim. CarSim uses mathematical models to implement vehicle motion and provides the necessary data for integrated vehicle motion control. MATLAB/Simulink is used to implement an integrated vehicle motion controller, while Amesim designs the E-powertrain, reducer, and SbW system. The SbW system used in this study was modeled on the basis of the actual model of the dual-motor steering system in the author's research group, and it is shown in Figure 9a. The performance of the rack force estimator was confirmed experimentally, as shown in Figure 9b, using a load cell sensor and a load motor.

Simulation Results
The parameters of the vehicle and controller are listed in Table 1.  The E-powertrain was a 4WD system consisting of an IWM system and a 4WS consisting of a front and rear SbW system. The simulation was performed in two ways: double lane change (DLC) and random steering and drving force commands (Random). The Random simulation use compared not only the proposed method and PID method for the higher-level controller but also the tire workload when four driving forces and two lateral forces were distributed using the yaw moment and the pseudo-inverse method. The percentage of the root mean square (RMS) errors of estimated forces are listed in Table 2. The percentage of RMS error is given by: where F act andF are the sensor values of CarSim and estimated force, respectively.

Results for Double Lane Change
The first simulation performed a DLC. The coefficient of friction between the road surface and the tire was 0.9 under dry asphalt road conditions. In this simulation, the steer signal was given at 7 s. The results are shown in Figures 10 and 11. The performance of the driving force observer and lateral force observer was satisfactory at a vehicle speed of 30 m/s. The trajectory of the vehicle is shown in Figure 10a, and the result of the lateral acceleration close to the maximum coefficient of friction between the road surface and the tire is shown in Figure 10b. Figure 10c shows that the high-level controller, yaw moment distribution, and low-level controller introduced in this study effectively track the yaw rate according to the driver's steering angle. The yaw moment generated by the high-level controller and the workload ratio calculated from the CarSim sensor are shown in Figure 10d,e, respectively. In Figure 10d, a larger yaw moment cannot be generated even though the vehicle yaw rate error increases because of the saturation function of the yaw stability controller. As shown in Figure 10e, the tire workload ratio confirms that the maximum friction force of each wheel is produced.  The driving force and lateral force distributed in the yaw moment distribution are shown in Figures 12a-d and 11c,d; the performance of the low-level controller can also be confirmed. In this simulation, the steering angle command and the actual steering angle at the wheel are shown in Figure 12e, and it can be confirmed that an additional steering angle is generated for the direct lateral force control. The rear wheel steering angle is shown in Figure 11a. The influence of the kingpin torque and vertical force on the lateral force estimation is shown in Figure 11b, and the influence of the driving force on the lateral force estimation in this simulation is negligible.

Results for Random
The second simulation performed random steering and random driving force. It verified the performance of yaw stability control and lateral force estimation according to the load transfer and driving force when the vehicle is accelerating and decelerating. The steering signal was given at 4 s. The simulation results are shown in Figures 13-15. The vehicle was driven between 18.5 m/s and 23 m/s, as shown in Figure 13a. The accel-eration of the vehicle is shown in Figure 13b; the vehicle motion control according to the acceleration and deceleration of the vehicle was successful. As mentioned above, the results of the proposed method and the PID method for the high-level controller are shown in Figure 13c, and the error of each method is shown in Figure 13d.    To verify the vehicle's stability, tire workloads were compared, and the proposed algorithm results are shown in Figure 13e. Figure 14a shows the tire workload means of the distributed algorithm with four driving forces and two lateral forces using the proposed yaw moment distribution and pseudo-inverse method. It can be seen that the tire workload means of the proposed algorithm had a significantly smaller value than the pseudo-inverse at 11, 16.3, and 17.5 s. The standard deviation of the tire workload of each wheel is shown in Figure 14b; this shows that the proposed algorithm has similar values for the load on each wheel.
The performance of the low-level controller, according to the driver's steering and driving force command, is shown in Figures 14c and 15a,d,e. It can be confirmed that the performance of the low-level controller was satisfactory. In Figure 15b, it can be seen that the steering angle is generated differently from the driver's steering command because of the direct lateral force control. The steering angle according to the direct lateral force control of the rear wheel is shown in Figure 15c.

Conclusions
This paper proposed vehicle motion control based on three layers, namely a high-level controller, yaw moment distribution, and a low-level controller, to improve the handling of EV equipped with a front/rear SbW system and 4WD of the IWM system. The control method was to improve the stability of the vehicle by minimizing the tire workload for the 3-DOF motion of the vehicle, and this was realized by directly controlling the driving force and the lateral force. The high-level controller generates a yaw moment to control the yaw rate of the vehicle, and the yaw moment distribution was introduced to distribute the driving force of each wheel and the lateral force of the front and rear wheels. The low-level controller directly controlled the distributed driving force and lateral force by driving force control and lateral force control, respectively. In this paper, we designed a driving force observer and a lateral force observer that do not use a tire model to estimate the driving force and lateral force. The control performances were verified via computer simulation using CarSim, MATLAB/Simulink, and Amesim. The simulation was conducted with double lane change as well as random steering and driving force commands. In future works, the states of the vehicle are estimated using the estimated driving force and lateral force.
Author Contributions: Y.S. and K.C. have equally contributed to this paper, performed the numerical simulations, and proposed control algorithms. K.N. contributed considerable effort in conducting the HILs test and organizing the manuscript as a corresponding author. All authors have read and agreed to the published version of the manuscript.
Funding: This paper was supported by the Korea Institute for Advancement of Technology (KIAT) grant funded by the Korea Government (MOTID)(P0008473, HRD Program for Industrial Innovation) and and in part by the Yeungnam University Research Grant 220A380108.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix B
The vector r perpendicular to the kingpin axis and passing through the point of force acting on the tire is equal to the direction of the moment arm of the kingpin torque. The intersection point (p) of vector ( r) and vector ( a k ) is the support point of the moment arm, shown in Figure 8. The x, y, and z of p are as follows: The direction vector ( v) of r is as follows: The vectors a k and the vectors r are perpendicular, and the vectors u and the vectors v must be orthogonal. u · v = 1 tan(Φ + 90) , 1 tan(Ψ + 90) , 1 · l, m, n = 0