Similar articles

Abstract: This paper presents a torque distribution method for an individual-wheel drive vehicle, in which each wheel is controlled individually by its own electric motor. The terminal sliding mode technique is employed for the motion control so as to track the desired vehicle motion obtained by interpreting the driver’s commands. Thus, finite-time convergence of the system’s dynamic errors can be achieved on the terminal sliding manifolds, as compared to the well-used linear sliding surface. By considering nonlinear constraints of the tire adhesive limits, a simple yet effective distribution strategy is introduced to allocate the motion control efforts to each of the four wheels. Through the use of a high-fidelity CarSim full-vehicle model, vehicle stability and handling performance of the proposed controller is evaluated in both open- and closed-loop simulations.

基于终端滑模控制的四轮独立驱动汽车转矩分配方法

研究目的：四轮独立驱动汽车四个车轮的电机转矩可以正、反向输出，有必要对其底盘集成控制系统进行针对性设计以保证各车轮间的协调运作。集成控制系统一般可以分为上层运动控制器和下层力分配器。对于运动控制器的设计，同类研究一般采用滑模控制方法来处理汽车运动的非线性特征，其中终端滑模控制因具有高的稳态精度和有限时间收敛的特点而成为研究热点。对于力分配器的设计，通常方法不能在保证运算效率的同时考虑到执行器的约束，从而很难应用于实际。本文采用终端滑模控制方法来设计运动控制器，通过分析驾驶员操作行为从而更好地追踪理想的车辆运动目标；并提出一种简单有效的转矩分配控制策略，通过考虑轮胎附着极限从而将运动总力分配至四个车轮上。
创新要点：本文创新性地将终端滑模控制应用到底盘集成控制系统以实现车辆纵向、侧向及横摆运动的联合控制；本文提出了一种新颖的转矩分配控制策略，将复杂的有约束控制分配问题分解至若干个简单的无约束分配子问题。
研究方法：本文采用分层式协调控制方案（图1），应用非奇异和全局快速终端滑模控制方法设计运动控制器，提出一种基于伪逆矩阵的有约束转矩分配策略，应用驾驶员最优预瞄加速度模型来描述和分析人-车闭环系统的运动响应，通过MATLAB/Simulink和CarSim的联合仿真对所设计的四轮独立驱动汽车转矩分配方法进行对比验证。所设定的三个仿真工况包括：开环方向盘角阶跃输入（图6、7及表2）、闭环双移线工况（图8-10）和闭环对开路面制动（图11-15），分别用以测试车辆横摆、侧向及纵向方向上的动力学响应。
重要结论：本文提出的基于终端滑模控制的四轮独立驱动汽车转矩分配方法将车辆的稳态转向特性由不足转向转变为中性转向，驾驶员从而能更容易地操纵车辆而不需要对汽车的非线性响应做出额外的转向补偿。与此同时，该分配方法在不影响车辆侧向稳定性的前提下可以准确地响应驾驶员的加速/制动意图。仿真结果表明终端滑模控制器相比传统滑模控制方法在车辆纵向、侧向、横摆方向上的运动控制效果均有一定程度地提高，而有约束的力分配器更可以明显地提高车辆的操纵性和稳定性。综合评价表明，所提出的基于终端滑模控制的四轮独立驱动汽车转矩分配方法获得了最佳的控制性能，符合设计要求。
汽车动力学；底盘集成控制；驾驶员模型；终端滑模控制；控制分配

1.  Introduction

 With the rapid development of electric motors and electronic communications, there is a trend in the automotive industry towards the use of electrified propulsion systems together with electronic control systems to achieve fuel efficiency and to improve handling performance. Almost all automakers throughout the world have been studying and promoting electric vehicles (EVs) for energy and environmental considerations. Individual-wheel drive (IWD) is a promising technique for EVs in which all wheels are controlled independently by their own in-wheel/hub motors (Hori, 2004).

 Because of the inherent characteristics of electric motors in which the torque can be negative, it is necessary for the global chassis control (GCC) system to coordinate the individual traction and braking actions of the four wheels. There are two main tasks in the GCC system (Johansen and Fossen, 2013). The first one is to compute a vector of virtual inputs to be applied to the vehicle so as to meet the motion control objectives of a given application. The other is to allocate these efforts to the individual actuators such that the total forces and moments generated by all actuators amount to the motion control efforts. In the design of the motion controller, the nonlinearity of the vehicle dynamics is usually handled by employing the sliding mode control (SMC) method (Mokhiamar and Abe, 2004; Li et al., 2008; Wang and Longoria, 2009). A terminal sliding mode control (TSMC) is a type of sliding control, in which the finite-time convergence and quick responsiveness can be achieved on the terminal sliding manifold (Liu and Wang, 2012). It has been successfully applied to applications of vehicle longitudinal (Kim et al., 2007) and lateral (Ren et al., 2011) control, but no study has been previously made regarding its effects on the combined motion control. In the matter of control allocation, different optimization methods have been proposed (Mokhiamar and Abe, 2004; Ono et al., 2006; Li et al., 2008; Wang and Longoria, 2009; Chen and Wang, 2011), although none of them are computationally efficient. An alternative plan is to employ the pseudo-inverse method to map the relationship between the motion control efforts and the actuator efforts (Weiskircher and Müller, 2012), but this does not consider the actuator saturation and consequently has limited practical utility.

 In this paper, a novel torque distribution controller is developed for an IWD EV (Fig. 1). First, the bicycle model and the optimal preview acceleration (OPA) driver model is used to provide the steering wheel angles δ _SW in the open- and closed-loop maneuvers to compute the desired vehicle states V _xd, , and ψ _d. Second, a terminal sliding mode controller (TSMCr) is designed to generate the motion control efforts, i.e., F _xd, F _yd, and M _zd, to track the desired motion commands, in which the steady-state tracking errors can be eliminated in finite time. Third, while the longitudinal tire forces of the IWD vehicles can be controlled individually, for the trade-off between computational precision and real-time performance, this paper employs a simplified tire model to effectively calculate lateral tire forces with acceptable accuracy. Thus, nonlinear constraints of the tire adhesive limits can be determined accordingly. Last, by considering these constraints, a distributor effectively allocates the motion control efforts to each of the four wheels based on a novel strategy, which divides the task into several unconstrained distributions. It should be noted that for concerns of cost, maintenance and safety, no four wheel steering (4WS) or active front steering (AFS) system is involved and the only controlled variables are the four wheels’ traction/braking torques.

2.  System modeling

2.1.  Vehicle dynamics model

 For the model-based controller design, the vehicle is modeled as a 3-DOF (three degrees of freedom) double track model. Its coordinate system is shown in Fig. 2. (1) , (2) , where V_x and V_y are the longitudinal and lateral velocities of the vehicle at cog (center of gravity), respectively; ψ is the yaw angle, the first- and second-order time derivatives of which are the yaw rate and yaw angular acceleration, respectively; C _D is the drag coefficient factor; A _f is the vehicle front area; ρ is the air density; m is the vehicle mass; I_z is the moment of inertia around the vertical axis; the motion control efforts F _xd, F _yd, and M _zd are the desired longitudinal force, lateral force, and yaw moment, respectively; F _xij is the tire longitudinal force, which is controlled timely (due to the fast response of the electric motor) and precisely (via an accurate measurement of the motor torque); F _yij is the tire lateral force, which can be calculated by a tire model; the first subscript of (·)_ij denotes the left (l) or right (r) side of the vehicle, and the second denotes the front (f) or rear (r) axle. The corresponding matrices are (3) , (4) , where l _f and l _r are the distances from cog to front and rear axles, respectively; t _f and t _r are the front and rear wheel treads, respectively; δ _lf and δ _rf are respectively the steering angles of the left-front and right-front wheels, where the Ackerman steering principle is applied, expressed as follows: (5) , where the minus (plus) sign is for the left (right) wheel; L is the wheelbase of the vehicle, L=l _f+l _r; δ _f is the front wheel steering angle.




 A simple yet powerful tool to describe the vehicle lateral dynamics is the well-known bicycle model, as shown in Fig. 2. The matrix form of this model is expressed as Eq. (6), in which the left and right side wheels are combined to one wheel, and the lateral tire forces F _yf and F _yr are proportional to the tire sideslip angles α _f and α _r, respectively. (6) , where C _αf and C _αr are the tire cornering stiffnesses of the front and rear axles, respectively.

 The relationship between the vehicle sideslip angle β and the lateral velocity can be expressed as (7) .

2.2.  Wheel model

 With the use of electric motors to drive all the wheels, an IWD EV can control the longitudinal tire forces individually based on the following equation: (8) , where T _ij is the traction/braking torque exerted by the electric motor; T _bij is the friction torque of the braking system and is assumed to be known; J _ij is the wheel spin inertia; is the wheel spin acceleration; R _ij is the wheel effective rolling radius; f _ij is the rolling resistance coefficient; F _zij is the tire vertical load, which can be expressed as follows: (9) where the minus (plus) sign is for the left (right) tire; g is the gravitational constant (9.8 m/s²); m _s is the sprung mass of the vehicle; k _f and k _r are the lateral weight-shift distribution ratios on the front and rear axles, respectively; h _s is the sprung mass height; a_x and a_y are the longitudinal and lateral accelerations at cog, respectively, and are given by (10)

2.3.  Tire model

 A simple tire model proposed by Sakai et al. (2002) is used to depict the relationship between the lateral force and the sideslip angle in the tire stable and monotonic region. By using the arctangent function to fit the experimental curve, the tire model can be formulated as follows: (11) where μ is the tire-road friction coefficient; C_α is the tire cornering stiffness; α is the tire sideslip angle, which is obtained by Eq. (12). Here, the value of the fitting constant i is changed to ‘2.9’ to improve the fitting accuracy at relatively large sideslip angles. (12) where σ _ij is the angle between the velocity of wheel center and the longitudinal axle of the vehicle body. They are given by (13) where the minus (plus) sign is for the left (right) wheel.

3.  Driver-vehicle system

3.1.  Reference model

 This study assumes that the driver’s task is only to keep the vehicle tracking at given speed profiles. So the desired longitudinal velocity is independent from the steering task: (14) where V _xd is the desired longitudinal velocity of the vehicle at cog; V _xr is the longitudinal velocity error; V _x0 is the initial longitudinal velocity of the vehicle at cog; a _xd is the driver acceleration or deceleration command.

 For stability consideration, the undesired vehicle sideslip motion that generates V_y should be suppressed to the minimum level, which can be effectively controlled by systems such as 4WS and AFS. However, these complex systems increase the overall cost and require additional maintenance. Moreover, since the vehicle’s lateral dynamics are greatly influenced by the active intervention of the AFS system, non-cooperative behaviors may occur between the human driver and the AFS system when they hold different targets (Na and Cole, 2013). This can cause an inexperienced driver to get stressed and increase his/her risk of a serious road accident. Therefore, we study the IWD EV with a conventional front-wheel-steering configuration and use a torque distribution approach to control the vehicle’s motion.

 The steady-state yaw rate response of the bicycle model in Eq. (6), i.e., steady-state yaw rate gain, , is obtained by (15) where K _s is defined as the understeer gradient, which usually has a negative value. It can be determined that the understeer gradient is actually a time-varying coefficient with the tire cornering stiffness varying at different vertical loads, road conditions, and slip ratios. Also, a non-zero K _s makes a nonlinear function of V_x . Furthermore, the driver tends to maneuver the vehicle according to its linear response, so that hereafter the value of K _s is set to be ‘0’, which implies that the steady-state yaw rate gain is V_x /L.

3.2.  Driver model

 The experiments should contain the open-loop maneuvers in order to verify the tracking performance of the proposed controller. Nevertheless, the controller should also be tested in the closed-loop maneuvers, so that the driver’s reactions to the vehicle dynamic responses can be studied and analyzed. This study employs the OPA driver model to mimic the human driver behavior (Guo and Guan, 1993). By using the schematic representation of Fig. 3 and assuming the vehicle proceeds with the current longitudinal and lateral velocities, the previewed lateral error between the desired path and the actual vehicle position ε is defined by (16) , where T _p is the preview time; r(t+T _p) is the desired path after T _p; y is the lateral offset of vehicle at cog.




 Fig. 4 shows the block diagram of this single-point preview model, which exerts the steering input in order to provide the optimal lateral acceleration and to minimize the previewed lateral error. Eq. (17) indicates that the OPA driver model predicts the vehicle’s future position based on the dynamic response of the aforementioned bicycle model. (17) where G _ay is the steady-state steering gain, which reflects the yaw dynamics of the internal model anticipated by the driver and is proportional to the aforementioned steady-state yaw rate gain; is the optimal front wheel steering angle. The driver’s actual input to the vehicle is the steering wheel angle δ _SW, given by (18) , where i _SW is the steering system ratio; T _C is the driver correction time; T _N is the driver action delay; T _D is the driver brain response delay; and s is the Laplace operator.




 For the lane-keeping task, the driver exerts steering effort not only to eliminate the deviation error, but also to minimize the orientation error with respect to the road. Moreover, it was discovered that accumulation of the yaw rate tracking error could cause the undesired vehicle yaw angle which greatly affects the vehicle heading and path-following accuracy (Wang and Longoria, 2009). Thus, a yaw angle control effort is considered to mitigate this effect. To sum up, the yaw error ψ _r and the yaw rate error are given by (19) where ψ _rd is the road path heading angle after T _P, as shown in Fig. 3.

4.  Terminal sliding mode control (TSMC)

 The global fast TSMC proposed by Park and Tsuji (1999) is employed to design the sliding surfaces s ₁ and s ₂, while the nonsingular TSMC technique proposed by Feng et al. (2002) is applied to design the sliding surface s ₃. The sliding surfaces are defined by (20) , (21) , (22) , where α _ieq and β _ieq are the strictly positive constants of TSMC (i=1, 2, 3); p_i and q_i are the positive odd numbers and p_i >q_i . In particular, for i=3, we also have: 1<p ₃/q ₃<2.

 Because the longitudinal and lateral velocities follow the dynamics of the first-order system in Eq. (1), the control laws u ₁ and u ₂ are directly derived from the sliding surfaces by letting s ₁ and s ₂ be equal to ‘0’, expressed as Eq. (23). Meanwhile, the yaw motion follows the dynamics of a second-order system, thus the control law is obtained from the differential equation, given by Eq. (24). (23) (24) where u _3n is the approaching law for the sliding surface s ₃. Here, the nonlinear reaching law u _3n proposed by Ren et al. (2011) rather than a linear one is adopted to achieve the finite-time convergence of and when , and thus allows us to have a better control performance. (25) , where α _3n and β _3n are strictly positive constants; p _3n and q _3n are positive odd numbers and p _3n>q _3n.

 To sum up, the control laws are given by (26) , (27) , (28) .

 The global fast sliding surface for the controllers u ₁ and u ₂ is expressed as Eq. (29), of which the convergence time is t _s when the initial state x(0)≠0 reaches x=0 which is given by Eq. (30). (29) (30) .

 By analyzing the solutions of the differential equation, expressed as Eq. (31), the finite-time convergence for the nonsingular TSMCr u ₃ occurs from the initial state s ₃(0)≠0 to s ₃=0. (31) .

 Hence, we have (32) .

 Then it is straightforward to determine that the yaw angle error ψ _r will reach ‘0’ along the sliding mode in finite time (Liu and Wang, 2012).

 For stability analysis, the Lyapunov functions are selected as (33)

 Then we have: (34) , (35) , (36) .

 Then the sliding conditions are satisfied and the sliding surfaces are attractive. Note that if the desired forces and moment cannot be generated due to the friction circle limitations of the four wheels, the stability of the proposed TSMCrs cannot be guaranteed. The friction circle of each tire would directly affect the range of the control efforts. The inequality constraint of the motion control efforts is written as (37) , where W _r is a constant weight balancing forces/moment, which is chosen as the reciprocal of the typical length between cog and road wheels (Ono et al., 2006). When the control efforts are calculated to be greater than this constraint, the motion control efforts would be limited by a ratio of τ _u. Note that this loose bound is to ensure that the desired motions are theoretically realizable before the efforts are distributed to the actuators, a precise limitation of the bound will be developed and handled in the torque distribution strategy, which will be illustrated in the next subsection.

5.  Constrained torque distribution

 The motion control efforts produced by TSMCr need to be distributed to the electric motors driving the four wheels. The functional relationship between the motion control efforts and the controlled variables can be compactly written as Eq. (38), which serves as a redundant linear constraint for the torque distribution problem. (38) where B is a 3×4 matrix representing the effectiveness of the actuators. There are 3 motion control objectives but only 4 controlled variables, which has limited potential for an optimal distribution. Furthermore, optimization approaches usually require large computational efforts and are not suitable for real-time applications. Thus, assuming that the pseudo-inverse B ⁺ exists, this study distributes the control efforts by calculating the torques with u _act= B ⁺ F _ud. The tire adhesive limit should also be considered as an additional constraint, which can be described by the friction ellipse concept (39) , where the values of the coefficients C _xij and C _yij are set to ‘1’ due to the use of the simple tire model in Eq. (11), which defines the friction limitation as a circle rather than an ellipse. However, it is still difficult and time-consuming to distribute torques directly under these nonlinear constraints of four tires.

 This study presents an efficient way to divide this problem into no more than 5 unconstrained distributions (Fig. 5). Note that for small lateral accelerations, all the 3 control efforts desired by TSMCr can be achieved; when F _yd is large, this force should be ignored, where the sideslip motion is suppressed mainly by the driver’s steering effort rather than by the controller itself; when 3 out of 4 wheels reach their adhesive limits, F _xd is also ignored and only M _zd is left to be achieved, which is thought to be more critical to the vehicle’s handling stability.




 Here, for notation simplification, u _{act_i} and F _{ud_j} can be expressed as Eqs. (40) and (41), respectively; B _ij is the corresponding matrices of u _{act_i} and F _{ud_j} . Particularly, B ₁ is the corresponding matrix of u _act and F _{ud_1}. The superscript (∙)⁺ denotes the pseudo-inverse operator. (40) (41)

 Once the tire longitudinal forces are obtained, the traction/braking torque of each wheel is calculated by Eq. (7).

6.  Results and discussion

 The proposed torque distribution method is designed in the MATLAB/Simulink environment. As depicted in Fig. 1, an A-Class Hatchback full-car model built in CarSim is employed here for the controller validation. Some of the vehicle’s parameters are listed in Table 1. Three simulation cases including both open- and closed-loop maneuvers are conducted. The responses of the vehicles with TSMCr and the traditional sliding mode controller (SMCr) proposed by Wang et al. (2009) are compared with each other. The unconstrained distribution approaches are also added for comparison, which does not take into consideration the tire adhesive limits.

6.1.  Open-loop step-steer maneuver

 In the open-loop simulation, the driver model is not involved, and the yaw error ψ _r is thus neglected. The vehicle is subjected to a step-like steering input at 1.0 s, with a front steering angle δ _f=2° (δ _SW=32.8°). The speed remains constant at 80 km/h and the road friction coefficient is 0.85. The vehicle dynamic responses are shown in Figs. 6 and 7 and Table 2.










 From Fig. 6, it can be seen that since the target yaw rate represents the neutral steering response, the conventional IWD vehicle without control tends to understeer and is not able to reach the yaw rate desired by the driver. For the vehicles using the unconstrained distribution methods, as soon as the steering commands are applied, their states oscillate in the time domain but can still be stabilized to the target values. This can be explained by the fact that in the high speed cornering case with a relatively large steering angle, the tire forces distributed by the unconstrained distribution method are beyond the tire adhesive limits, such that one or more tires are saturated and they cannot provide enough forces within the tire adhesive limits which can amount to the total yaw moment required by the motion controller. This problem can be effectively solved by the proposed constrained distribution strategy. Both vehicles with the constrained methods can accurately respond to the desired yaw rate, Table 2 shows that the constrained TSMCr provides the best control performance with the minimum error.

 As shown in Fig. 7, the suppression of the vehicle sideslip angles by these distribution methods is restricted due to the limited impacts on the lateral motion by the longitudinal tire forces, thus each controller acts more like a joint system of cruise control (CC) and direct yaw control (DYC), which controls the longitudinal and yaw motions simultaneously. However, the vehicles with the unconstrained methods lose their stabilities, while the constrained ones can still control the vehicle sideslip motion to be stable, whereas their sideslip angles are slightly larger than that of the conventional vehicle. The root mean square (RMS) values of the sideslip angles can be found in Table 2. Note that the counter steering effect of the conventional vehicle is eliminated by the proposed method, which does not produce the reversed sideslip angle and has a faster lateral response than without control.

6.2.  Closed-loop maneuver

 To further evaluate the control effects, closed-loop simulations are conducted and performed as the handling tests, in which the OPA model provides the driver’s steering inputs for the lane-keeping tasks. Note that the road path heading angle ψ _rd used to calculate the yaw error ψ _r in Eq. (19) is obtained by some road information acquisition algorithm, which is not presented in this study.

6.2.1.  Double-lane-change

 The first closed-loop simulation is a double-lane-change (DLC) maneuver, with a constant longitudinal velocity of 100 km/h and a road friction coefficient of 0.5. The results are shown in Figs. 8 to 10.










 It is shown in Fig. 8 that in the attempt to control the conventional vehicle to track the given path accurately, the OPA driver model cannot provide the appropriate steering effort due to the mismatch between the characteristics of the neutral steer anticipated by the driver and the handling response of the vehicle without control, and this therefore causes large lateral deviations from the ideal path. However, vehicles with both the unconstrained approaches do not lose lateral stability as they do in the step-steer maneuver and show fairly good path-following performance. Despite that, the constrained approaches provide good path-tracking performance. This is primarily because the constrained distribution methods can generate the motion control efforts determined respectively by TSMCr and SMCr as much as possible (especially when the vehicles are cornering with large lateral acceleration) within the tire adhesive limits. The RMS value of the lateral offset of the constrained TSMCr is 0.1693 m, compared to that of the constrained SMCr which is 0.2023 m.

 Fig. 9 shows the driver-vehicle system responses depicted in the phase plane, where the controller with the tighter trajectory but higher value of yaw rate response is supposed to maximize the stability margin and the vehicle responsiveness. It can be attributed to the aforementioned mismatch effect that the trajectory of the response of the vehicle without control appears much looser. Its maximum yaw rate and sideslip angle as well as those of the vehicles with unconstrained methods are located in the fourth quadrant of the plot, which indicates these vehicles tend to be less stable and perform poorly before entering into the second straight section of the DLC path. At the same time, the vehicle without control needs extra steering effort to be stabilized and is difficult for the driver to manipulate, as shown in Fig. 10. Both the constrained distribution methods attain good control performance, whereas the constrained TSMCr gives better handling performance (the tightest trajectory in Fig. 9) and increases maneuverability (the least steering effort in Fig. 10).

6.2.2.  Hard braking through a μ-split surface

 The second closed-loop simulation is the μ-split braking, which is mainly for the longitudinal dynamics validation. As shown in Fig. 11, the road is divided along a straight line with two different surfaces with μ=0.8 on the left side and μ=0.2 on the right side. The driver attempts to decelerate the vehicle from 120 km/h to 0 as quickly as possible and also steers the vehicle in order to stay in this straight lane.




 During the simulation testing, it was discovered that neither the two unconstrained methods nor the even torque distribution could stabilize the vehicle in this type of maneuver. Hence, the common and the laterally-stable anti-lock brake systems (ABS) are employed here to compare their control effects with each other (Mitschke and Wallentowitz, 1972). Figs. 11 and 12 show the results. It can be seen that although the maximum lateral offset of the laterally-stable ABS is only 0.0027 m, compared to that of the common one, which is 15.4727 m, it still takes 14.2 s for the former to stop the vehicle at 229.7168 m, more than twice as much braking time and distance as those of the common ABS. This is because, to avoid the additional yaw-moment when braking the vehicle, the laterally-stable ABS makes the longitudinal tire forces on both sides be equal to the lesser of the right side, and thus decreases its braking performance. For the common ABS, the vehicle drifts towards the μ-high side of the road, which causes shorter braking distance and time, meanwhile its maximum lateral velocity reaches up to 80 km/h. It can be concluded that both methods fail to balance the trade-off between the lateral stability and the braking efficiency.




 The constrained methods with TSMCr and SMCr are validated in Figs. 13 to 15. Fig. 13 shows that the proposed distribution strategy can fully utilize the road adhesion capability and both controllers provide vehicles with good braking efficiency (the vehicle with TSMCr stops in 8.25 s at 142.0866 m and the one with SMCr stops in 8.35 s at 142.7459 m) by not losing lateral stability (the maximum lateral offsets of TSMCr and SMCr are 0.5049 m and 0.6074 m, respectively) while the latter yields a slightly worse performance. The reason can be explained that with the same initial conditions the state errors of longitudinal and lateral velocities of the TSMCr converge to 0 more rapidly than those of the SMCr do, as shown in Fig. 14. Therefore, the vehicle with TSMCr has a shorter braking time with less lateral deviations which hence results in a shorter braking distance. Fig. 15 also suggests that due to its faster convergence of the yaw rate error, TSMCr offers a better tracking performance than SMCr.

7.  Conclusions

 In this paper, a novel torque distribution approach for an IWD vehicle is introduced in the design of the TSMC based controller to adjust and coordinate the torques from four electric motors. The proposed controller alters the handling characteristics of the vehicle from understeering to neutral steering so that the driver can manipulate the vehicle more easily without having his/her steering effort compensated for nonlinearities in vehicle yaw response. It also operates as a speed controller, which can accurately execute the driver’s acceleration and deceleration commands without degrading steering performance. Although no active steering system is used here, lateral stability is guaranteed and the sideslip motion is, to some extent, suppressed by the controller. Additionally, the constrained distribution methods deal with nonlinear constraints of the tire adhesive limit efficiently and effectively.

 Simulation results indicate that due to its fast and finite time convergence, and high steady-state precision, TSMCr provides more improved control effects than SMCr does. Moreover, the constrained distribution approach also significantly improves vehicle stability and handling performance. On an overall evaluation, the constrained TSMC based controller meets the design requirements and shows the best performance among all schemes during our study.

 Note that an IWD vehicle is a highly redundant system, which is more inclined to take the risk of actuator failure. Integrating fault tolerant control into this chassis control system will be performed in future work.

---

* Project supported by the Open Fund Project of the State Key Laboratory of Automotive Simulation and Control (No. 20120111), and the China Scholarship Council (No. 201306170075)