Model predictive control allocation for stability improvement of four-wheel drive electric vehicles in critical driving condition

Transcription

1 IET Control Theory & Applications Research Article Model predictive control allocation for stability improvement of four-wheel drive electric vehicles in critical driving condition Haiyan Zhao 1,2, Bingzhao Gao 1, Bingtao Ren 2, Hong Chen 1,2, Weiwen Deng 1 ISSN Received on 4th November 2014 Revised on 25th July 2015 Accepted on 24th August 2015 doi: /iet-cta State Key Laboratory of Automotive Simulation and Control, Jilin University, Changchun , Jilin, People s Republic of China 2 Department of Control Science and Engineering, Jilin University (Campus NanLing), Changchun , Jilin, People s Republic of China chenh@jlu.edu.cn Abstract: To improve the vehicle stability of an electric vehicle (EV) with four in-wheel motors, the authors investigate the use of a non-linear control allocation scheme based on model predictive control (MPC) for EVs. Such a strategy is useful in yaw stabilisation of the vehicle. The proposed allocation strategy allows a modularisation of the control task, such that an upper level control system specifies a desired yaw moment to work on the EVs, while the control allocation is used to determine control inputs for four driving motors by commanding appropriate wheel slips. To avoid unintended side effects, skidding or discomforting the driver in critical driving condition, the MPC method, which permits us to consider constraints of actuating motors and slip ratio, is proposed to deal with this challenging problem. An analytical approach for the proposed controller is given and applied to evaluate the handing and stability of EVs. The experimental results show that the designed MPC allocation algorithm for motor torque has better performance in real time, and the control performance can be guaranteed in the real-time environment. Nomenclature β, γ sideslip angle/yaw rate δ f, δ front wheel steering angle from controller/driver M z direct yaw moment input V, M vehicle velocity/mass F yf, F yr lateral force of front/rear tires h, R e height of centre of gravity/wheel radius I ω moment of inertia of the wheels a x, a y longitudinal/lateral acceleration I z yaw moment of inertia L f, L r distance from the centre of gravity to front/rear axle L distance from the front axle to rear axle d tread of vehicle C f, C r cornering stiffness at front/rear tire T i motor torque on the wheel K a, K b curve fitting coefficients α f, α r tire sideslip angles of the front/rear F xi wheel longitudinal force F zi static loads at front and rear wheels 1 Introduction In recent years, electric vehicles (EVs) driven by electric motors have recently been widely recognised and adopted in many countries as an effective mean to reduce pollution [1, 2]. Compared with traditional internal combustion engine drive vehicles (ICVs), motor torque generation is fast and accurate [3, 4]. This advantage enables us to realise high performance control of EVs. However, in four-wheel drive (4WD) case, due to four driving motors, control system has four permitted independent control inputs [5], which are all available to control yaw rate and sideslip angle. If we design control strategy for 4WD EVs, one issue faced is how to deal with such actuator redundancy, which means how to allocate torque for four driving motors [6]. The subject of resolving the redundancy of over-actuated systems has been an active field of research, see [7] for a comprehensive overview. Most existing algorithms neglect the actuator dynamics or deal with the actuator dynamics separately, thereby, a static relationship between actuator dynamic output and plant inputs is assumed [8]. Especially, some constrained control allocation methods are based on some scaling of the unconstrained optimal control allocation, such that the resulting control allocation is projected onto the boundary of the set of attainable generalised forces. The static control allocation (SCA) method starts with the unconstrained control allocation computed using some pseudo-inverse. Following the seminal work of [9], several algorithms for linear control allocation and adaptive control allocation have been developed [10 13]. In addition, actuator dynamics can be considered by model predictive control (MPC) allocation as shown by the authors [14, 15]. The main objective of this paper is to present a control allocation scheme for EVs, which permits us to consider the constraints of actuating motors, and simultaneously consider the characteristic of wheel slip. Building on the results of [16], this paper provides results for an MPC-based dynamic control allocation algorithm for asymptotically stable tracking of non-constant input commands. Here, we also propose an integrated yaw stability control algorithm for 4WD EV. The whole control strategy consists of two control loops. The upper controller is concise and comparable to the practice structure of the modern automotive control [16], which benefits to bridge the gap between advanced model-based control and engineering implementation. The main contribution is a real-time MPC allocation strategy to provide the additional yaw torque. The MPC problem herein is posed as a sequential quadratic programming problem with the constraints [17, 18]. In order to control tire slip on low-friction coefficient road, the soft constraint of tire longitudinal slip ratio is translated into a penalty function to enforce the slip changing within the stable slip range. To evaluate the improvement in performance achieved when EV drives in critical conditions, we also compare the MPC allocation scheme with the SCA scheme. The simulation results also extend those of [16] in that controller evaluation method is used, thereby increasing the importance of accounting for actuator dynamics, and highlighting the requirement for sophisticated control allocation. At last, to improve the real-time performance of controller based on MPC for torque allocation, a scheme for the controller implementation based on FPGA (field-programmable gate array) is adopted The Institution of Engineering and Technology 2015

2 3 System model for controller design To design control law, a non-linear model based on vehicle dynamics is developed in this section, which mainly consists of vehicle, tire and in-wheel motor. The symbols in modelling procedure and the related physical meanings are listed in Nomenclature section. To design the controller, a widely used simplified single track vehicle model is considered to capture the essential vehicle lateral steering dynamics as follows Fig. 1 Vehicle dynamics control structure 2 Integrated yaw stability control algorithm The concept of the proposed control algorithm is shown in Fig. 1. The upper control loop is used for calculating virtual control effectors M z, which is required for realising steady and dynamic characteristics of control system. The control allocation task is to distribute this generalised force to the individual tyres, in an optimal manner by mapping M z to physical control inputs for motor. The controlled output variables are γ and β, and a reference model is used for calculating desired input commands γ r and β r in accordance with driver s commands. A reference model is designed based on the linear two-degreesof-freedom vehicle model, which is formed as follows β = 2C f (δ β {L f γ/v }) + 2C r ( β +{L r γ/v }) γ, MV (1a) γ = 2L f C f (δ β {L f γ/v }) 2L r C r ( β +{L r γ/v }), I z (1b) where the responses of γ and β to the front wheel steering input δ are second-order systems. From the perspective of frequency response analysis, we derive the transfer function from steering angle δ to γ as mentioned in [19] (see (2)) where ω n is the natural frequency, ξ is the damping coefficient and G γ (s) is the steady gains transfer function. To achieve an ideal closed-loop performance, the reference model is defined as follows ωd 2 γ r (s) = G kγ (s) s 2 + 2ω d ξ d s + ωd 2 δ(s), (3) where ω d = k ωn ω n, ξ d = k xi ξ and G kγ (s) = k g G γ (s). To improve the yaw rate phase delay and make response quick, we should increase the bandwidth of system, which means frequency factor k ωn > 1. To reduce the volatility of second-order system response, the damping coefficient should be larger than the previous, as damping factor k xi > 1. Due to the conflict of rapidity and stationarity performance exists, k ωn and k xi need to coordinate. In addition, γ r should be considered by its adhesion saturation value γ r μg/v, and the sideslip angle β is desired to be zero to improve handling stability, i.e. β r = 0. It should be noted that the modularisation offered by the control allocation would to a large extent simplify the design of such a controller. β = F yf + F yr γ, MV (4a) γ = L f F yf L r F yr + M z. I z (4b) Assuming that the longitudinal speed V is constant, only the lateral and yaw motion of the vehicle is taken into account. 3.1 Tire model A widely used semi-empirical tyre model to calculate steady-state tyre force and moment characteristics for use in vehicle dynamics studies is based on the so called Magic formula [20]. Tire lateral forces can be calculated as follows F yi = D y sin(c y arctan(b y α i E y (B y α i arctan B y α i ))), i = f, r (5) where B y, C y, D y and E y all depend on the normal vertical forces F z,sothef y is a function of α and F z. To simplify the model and keep tire non-linear characteristic, we expand the magic formula of F y with the Taylor expansion method, and obtain a cubic polynomial as follows F yf = 2C f (1 K a αf 2 )α f, F yr = 2C r (1 K b αr 2 )α r, (6a) (6b) where K a, K b are coefficients, which can be obtained by the Taylor expansion of the Magic formula model of tires. Tire side slip angles of the front and rear wheels are calculated as α f = β + (L f /V )γ δ and α r = β (L r /V )γ. The practical longitudinal slip ratio κ is defined as follows κ = v sx v x = v x ωr e v x, (7) where v sx, v x are the slip velocity and the wheel velocity, respectively. To calculate the required motor torque, a torque balance equation is formed for each wheel. The motion dynamics of a single wheel are calculated as follows 1 4 M V = F x, (8a) J ω = T e T b R e F x, (8b) where F x is a tire longitudinal force and T e is the motor drive torque in each wheel. According to (7) and (8), we assume that γ(s) δ(s) = (2C f L f /I z )s +[{4C f C r (L f + L r )}/MVI z ] s 2 + ([{2(C f Lf 2 + C rlr 2)}/VI z]+[{2(c f + C r )}/MV ])s +[{4C f C r (L f + L r ) 2 2MV 2 (C f L f C r L r )}/MV 2 I z ], ωn 2 = G γ (s) s 2 + 2ω n ξ s + ωn 2. (2) The Institution of Engineering and Technology

3 wheel longitudinal velocity v x is equal to vehicle longitudinal velocity V, so the wheel slip ratio dynamics is derived to κ = ωr ev ωr e V V 2 = (T e R e F x T b )R e JV ( = R2 e JV κ + 1 (1/4)MV ωr ef x (1/4)MV 2 ) F x + R e JV (T e T b ). Here the mechanical braking mode is not considered, so T b = 0. To calculate force easily, the formula can be simplified as, F x = C k κ, where C k is decided by tire vertical load F z as follows C k = F z (p k1 + p k2 df z ) e p k3 dfz λ k, df z = F z F z0 F z0, F z0 = 4000 N, (9) (10) where the data of factors p k1, p k2 and p k3 are derived from [20]. As mentioned above, the slip derivative dynamics can easily be formulated as follows κ = ( R2 e JV κ + 1 (1/4)MV ) C k κ + R e JV T e. (11) It is important to note that, although all of the variables mentioned above are accessible in a simulation environment, they must be measured or estimated in real life. 3.2 In-wheel motor model We assume that there is a motor torque controller for each wheel and the ideal motor closed-loop dynamics is simplified as follows T e = 1 τs + 1 T c, (12) where τ is the closed-loop response time which is a control characteristic of motor torque controller, and T c is the motor driving torque command. Due to motor torque saturation, the constraints of the control inputs T c can be written as follows T emax T c T emax, (13) where T emax is the maximum torque output of the motor. Then, the total drive torque T t from driver pedal command determines the driving torque of four motors, which is an equality constraint as follows T t = T cfl + T cfr + T crl + T crr. (14) 4 Yaw moment control allocation If we design control strategy for 4WD EVs, one issue faced is how to deal with such actuator redundancy, which leads to a control allocation problem. Due to the strong non-linearities and dynamic constraints, the use of MPC techniques will generally be desired for lateral vehicle control. 4.1 Design of optimisation function Combine (9), (12) and (13), we get the state-space model as follows ( ẋ i = R2 e JV x ) i + 1 C ki x i + R e {(1/4)MV } JV u i, y = d 2 ( C kflx fl + C kfr x fr C krl x rl + C krr x rr ), i = fl, fr, rl, rr (15a) where x =[κ fl, κ fr, κ rl, κ rr ] T, u =[T cfl, T cfr, T crl, T crr ] T, y = M z. (15b) The main task is to determine control input commands for four motors to (i) keep the tire from slipping or locking; (ii) obtain the maximum tire load ratio; (iii) energy consumption as small as possible; and (iv) improve yaw rate stability. The sample time k is chosen as k = int(t /T s ), where T is running time and T s is a fixed sampling step. The system model (15) can be described as a discrete non-linear form as follows x(k + 1) = f k (x(k), u(k)) T s + x(k), y(k) = C x(k), (16a) (16b) where output matrix is C = d [ ] Ckfl C kfr C krl C krr. 2 According to the main principle of MPC, p is defined as predictive horizon to extend the forecast of the system future output and m is defined as control horizon and set m = p for simplification. Usually horizon p m, and here p is simply set to be equal to m. The values of p and m determined by trial and error not only ensure control performance, but also reduce optimisation burden. The future control input sequence U(k) is optimisation vector in the following cost function and the predictive control output sequence is Y(k), which can be defined as follows u(k k) u(k + 1 k) U(k) =. u(k + m 1 k) r(k) r(k) R(k) =.. r(k) p 1 y(k + 1 k) y(k + 2 k), Y(k) =,. m 1 y(k + p k) p 1 where each vector of U(k) is an array of control inputs u, and each vector of Y(k) is an array of system outputs y, R(k) is defined as output reference sequence containing m arrays of virtual control input M z which keeps the same in each predictive horizon. The main control requirement is to ensure the EV good handling and stability on cornering, i.e. to make actual yaw rate γ and sideslip angle β track the required reference output γ r and β r quickly. Optimum torque allocation control is used to convert virtual control effector to physical control inputs for available sets of motors. So we define the first cost function as follows J 1 = Y(k) R(k) 2 Q p = [(M zr (k + i k) M z (k)) 2 Q], (17) i=1 where M zr is the reference control input, Q is the positive weight factors for adjusting tracking performance. As the tyre longitudinal slip κ directly reflects the tyre slip performance, we should restrict κ with a soft constraint of system. However, this output constraint of κ in this paper is a non-linear constraint for the optimisation problem, and it will affect the speed of solving. Hence, we consider a penalty function for anti-slip instead of slip soft constraints to facilitate fast solution, which enforces κ to remain within the stable slip range and is calculated as follows κ i (k) κ max, κ i (k) κ max, κ max E i (k) = 0, κ max κ i (k) κ max, (18) κ i (k) + κ max, κ i (k) κ max, κ max where the subscript i denotes fl, fr, rl, rr for each wheel The Institution of Engineering and Technology 2015

4 In addition, we also consider energy consuming problem of the motors. As the larger motor torque means more energy provided from the battery, we should make the sum of torque command squares Tc 2 as small as possible to save energy while ensuring stable vehicle motion. Therefore, we define an additional cost function in a form as follows J 2 = U(k) 2 R m 1 = (T cfl (k + j j) 2 + T cfr (k + j j) 2 j=0 + T crl (k + j j) 2 + T crr (k + j j) 2 ) R, (19) where R = diag(r 0, R 0, R 0, R 0 ) is positive weight factor to build a suitable form of the normalised function for T ci, and affect the gaps among four motor torque commands. To reduce the rate of change of the control action and ensure driving of the vehicle, a cost function is considered to control the rate of change of the control action as follows where P is a positive penalty matrix. It should be noted that terminal penalty matrix P can be determined by offline, which can reduce the computing time of controller. Equation (21) guarantees closed-loop stability by extending the finite prediction horizon to infinite horizon. 4.2 Total control objective To this end, the total control objective is written as follows J mpc (x(k), U(k)) = J 1 + J 2 + J 3 + J 4, = Y(k) R(k) 2 Q + U(k) 2 R subject to (16), (18) and + U(k) 2 S + x(k + p k) 2 P, (22) T emax T ci (k + j k) T emax, i = fl, fr, rl, rr, j = 0, 1, 2,..., m 1. (23) J 3 = U(k) 2 S m 1 = ( T cfl (k + j j) 2 + T cfr (k + j j) 2 j=0 + T crl (k + j j) 2 + T crr (k + j j) 2 ) S, (20) where the control input changing sequence U(k) is calculated by u(k) = u(k) u(k 1) and S = diag(s 0, S 0, S 0, S 0 ) is positive weight factor to coordinate the rate of change of control inputs T ci. To extend the finite prediction horizon to infinite horizon, a terminal cost is needed to penalise the states at the end of the finite horizon. This quadratic terminal cost in (21) controls the non-linear EV system by local linear state feedback to guarantee closed-loop stability [18] J 4 = x(k + p k) 2 P = x(k + p k)t P x(k + p k), (21) T t = T cfl (k + j k) + T cfr (k + j k) + T crl (k + j k) + T crr (k + j k), j = 0, 1, 2,..., m 1. (24) The objective function in (22) includes four conflicting contributions. The relative importance of the objective function contributions is controlled by setting the time dependent weight matrices Q, R, S, P, these are chosen to be positive definite. 5 Simulation results For verifying the derived model and testing the designed controller, a more realistic electric vehicle system simulation model is established in AMESim environment as shown in Fig. 2. The AMESim model include body model, electric wheel model, in-wheel model, steering model, suspension model, battery model and so on. Fig. 2 Electric vehicle model based on AMESim The Institution of Engineering and Technology

5 5.1 Simulation and parameter tuning Table 1 Parameters of virtual vehicle model In this study, the simulation time T = 15 s, the fixed sample step T s = 0.01 s. According to the data of vehicle, we set the factors of expected response γ r, β r of ideal vehicle, weight factors, prediction horizon and control horizon as follows k ωn = 1.9, k xi = 1.3, k g = 0.83, ω d = 8.47, ξ d = Q = 8550, R 0 = 7500, S 0 = 1050, p = m = 3. (25) Based on the weight factors and data of vehicle as mentioned above, the terminal penalty matrix P is calculated by solving the Lyapunov equation (21) as follows P = Setting these parameters requires multiple sets of tests and once determined, do not need to be updated online to adapt to different scenarios. The weight factor Q plays regulatory roles in the following the expectations to maintain vehicle stable. R 0 affects the gaps among four motor torque commands. With the increasing of R 0, the differences between four torques become smaller. S 0 is positive weight factors to coordinate the rate of change of control inputs δ f and T c. After several commissioning tests, enlarging S 0 can make change of control inputs smoothly. In the following sections, the responses are compared with that in the case without any control and the case with SCA. Parameters Symbol Unit Value Symbol Unit Value I z kgm m kg d m h m L f m L r m C r N/rad C f N/rad 23,540 L mh 0.05 R om Re m 0.29 ρ f Wb 0.1 κ max 0.08 T emax Nm 500 of the actual vehicle model used in the simulation are shown in Table Simulation results and analysis The standard double lane-change (DLC) manoeuvre, ISO :1999(E), is established in AMESim. The following scenario is simulated: the first one (0 70 m) is the beginning phase with speed of 80 km/h in straight road topology, the second phase ( m) presents the double lane road, in the third phase (after 200 m), the EV moving up the sloped road of 17% (high friction surface) and 8% (low friction surface). From the above simulation scenarios, all the dynamic characteristics of the vehicle can be adequately described, and the yaw stability of EV and energy performance of the motor can be fully verified High friction surface μ = 0.8: The trajectory responses and detailed dynamic response are shown in Fig. 3. EV Fig. 3 Dynamic response of EV: μ = 0.8 a Sideslip angle and yaw rate b Displacement X Y Z of EV c Torque distribution: MPC allocation d Torque distribution: static control 2692 The Institution of Engineering and Technology 2015

6 with SCA and proposed controls can pass the test track with similar responses. Compared to the SCA control, the proposed control shows slightly smaller yaw rate tracking error and sideslip angle deviation. Without control, the vehicle can also driving normally just shows the most oscillatory responses. As shown in Fig. 3b, the vehicle begin climbing at 10 s, the vehicle velocity will be reduced and the motor torque begin to increase for overcoming the influence of slope. Thus, the vibration of the curve at 10 s is large. From Figs. 3c and d, we can see that the motor torque keep close to saturation because of the effort of the control input constraints. Simultaneously, the optimal performance can be obtained and the proposed control algorithm shows smaller energy consumption and good climbing ability Low friction surface μ = 0.3: To further evaluate the anti-slip performance of the proposed method, the EV model is driven on low-friction coefficient roads through DLC manoeuvre, which is closed-loop test manoeuvres used to evaluate the lateral dynamic of a vehicle. The simulation is carried out at an initial velocity of about 60 km/h. The road surface is assumed to be wet road with a friction coefficient of 0.3. The trajectory responses and detailed dynamic response are shown in Fig. 4. From Fig. 4a, we can seen that the vehicle without control shows more oscillatory responses and loses the yaw motion stability. The proposed control can track the desired yaw rate while maintaining the sideslip angle within a certain range. Although the SCA control can stabilise the yaw motion, its yaw rate tracking error and sideslip angle deviation are larger than those of the proposed control. It can be seen from Fig. 4b, vehicle trajectory slightly deviates from the desired trajectory. The control input as active steering angle δ f is shown in Figs. 4c and d. It can be seen, MPC allocation strategy can better satisfy the driver s operation intention and has the better handling performance than SCA strategy. In addition, the result Fig. 5a illustrates that tire longitudinal slip ratio spread quickly and the κ on no-control condition is beyond 0.3, which means that the ground cannot provide enough static friction force for tire. Compared with this, the tire longitudinal slip ratio κ is controlled in the range of 0.02 by MPC allocation and control performance is better than SCA strategy. Distributed torque command T c as functions of time is shown in Fig. 5b. Compared to the SCA strategy, the proposed control shows smaller energy consumption and the proposed controller has corrected the steering angle δ from driver and the motor torque T c is under motor constraint T emax. 5.3 Controller evaluation This section deals with three major performances, the controller s performance, the driver s mental workload and the motor actuator s workload as determining factors for vehicle handling and stability qualities. The task performance J1 e is calculated by the following equation as the integral of the square of an error between the desired path and the actual state of the vehicle t J1 e = β 0(e 2 + e2 γ ) dt (26) where e β = β β r, e γ = γ γ r. Apparently, the smaller the value of J1 e is, the higher is the task performance. The whole control performance is improved by coordination of the upper and lower controllers, i.e. the value of J1 e is associated with control objectives J 1, J 2, J 3, J 4. Fig. 4 Dynamic response of EV: μ = 0.3 a Sideslip angle and yaw rate b Displacement X Y Z of EV c Steering angle d Local expansion of steering angle The Institution of Engineering and Technology

7 Fig. 5 Response of slip ratio and torques: μ = 0.3 a Tire longitudinal slip ratio b Motor torques Accordingly, the magnitude of the driver s physical workload can be determined by the integral of the square of the steering wheel angle as follows [21] J e 2 = t 0(δ δ f ) 2 dt. (27) This means that the smaller the value of J2 e is, the lighter is the driver s physical workload. Moreover, the value of J2 e reflect the control performance of the upper controller. It is known that in general the frequent adjustment for motors will cause the motor to overheat and result in thermal protection. Simultaneously, motor actuator wear is inevitable. In addition, the level of smooth operation for motor actuator can be determined by the integral of the square of the rate of torque change as follows t J3 e = fl 0( T 2 + T fr 2 + T rl 2 + T rr 2 ) dt, (28) where T ij, ij = fl, fr, rl, rr. is the rate of torque change. This means that the smaller the value of J3 e is, the lighter the frequent operation or the motor actuator wear is. From the perspective of energy saving, the value of J3 e is associated with the selection of weight coefficient of J 3. It is worth emphasising that these evaluation criteria can be used for designing the model predictive control (MPC) allocation controller. According to different control requirements, the weight coefficient of the control objectives J 1, J 2, J 3, J 4 need to be correspondingly adjusted. We randomly select 50 sets of control parameters within a certain range for experiments under double lane manoeuvre. Table 2 shows the values of the J1 e, J 2 e, J 3 e. From Table 2, we can know that the two kinds of controllers can help drivers successfully complete the test in double lane condition, but the MPC allocation (MPCA) shows the best performances. On the one hand, improve the control performance by 47%; on the other hand, improve driver s physical workload by 41%. In particular, Table 2 Performance analysis of the considered controller Parameter MPCA SCA Improvement ratio, % J1 e, rad s J2 e, rad2 s J3 e, Nm the MPC allocation strategy can relieve the motor actuator wear and improve the smooth operation by 86%. Thus, MPC allocation plays an important role for the vehicle s handling and stability performance. 5.4 Experimental results To solve the real-time problem of MPC in miniaturised device and embedded system, the MPC controller is built on an FPGA chip. The FPGA implementation allows the particle swarm optimisation (PSO) algorithm to exploit its parallel search capabilities. Therefore, the idea of using PSO to solve non-linear optimisation problem in MPC will be a good alternation. The PSO is employed to achieve realtime operation due to its naturally parallel capabilities. The proposed FPGA-based MPC-PSO controller consists of random number generator, fixed-point arithmetic, PSO solver and universal asynchronous receiver/transmitter (UART) communication interface. The PSO algorithm has three loops which account for most of the computation at each solution; (i) update the particles from generation to generation; (ii) update every particle in one generation; and (iii) update one-dimension of a particle. The FPGA chip is selected as Altera Stratix III EP3SL150F1152 and the clock frequency is set to 50 MHz. The FPGA implementation architecture of the designed MPC-PSO controller for torque allocation control is shown in Fig. 6a. For validating the FPGA-based MPC controller, the testing platform is based on an FPGA board and an xpc-target real-time system is built up and shown in Fig. 6b. The FPGA board is 2694 The Institution of Engineering and Technology 2015

8 Fig. 6 Hardware implementation a Implementation architecture b Test platform the hardware implementation of the MPC-PSO controller. The EV model is established in host PC and then is downloaded to the target PC after compiling. PC1 runs the FPGA design tools (Quartus II and Catapult C) and is used to design the MPC-PSO controller. The result data exchanged between the FPGA-based controller and the EV plant on target PC are transferred through UART RS232 serial link. The synthesis results for the FPGA implementation of the PSO-v1 architecture are presented in Table 3. It is worthy noting that the resource consumption and performance before optimisation and after optimisation are different and the speed of solving the problem is improved greatly at the cost of Table 3 Synthesis results of the PSO-v1 architecture Architecture Cycles Time, ms Slack Area (LUTs) PSO-v1 344, ,939 more resource consumption of the FPGA. In this paper, we carry out before optimisation. Experimental parameters are the same as the condition when μ = 0.3 without sloped road. The real-time simulation results are Fig. 7 Experimental results: μ = 0.3 a Sideslip angle b Yaw rate c Slip ratio response d Torque response The Institution of Engineering and Technology

9 shown in Fig. 7. It can be seen that the designed shifting controller meets the demand of real time and control. However, compared with the simulation results shown in Section 5.2, there is a certain tracking deviation. First, the optimisation solver is NAG-E04WD when offline simulation, the optimisation tools PSO is used to design the controller when FPGA hardware implementation. Thus, different solvers lead to a certain bias. Second, and noise always exists when the FPGA communicates with the xpc-target. In addition, the quantisation errors result in low measuring precision of the system output and the controller output. Last, the time needed by once solution of controller calculations is close to the sampling time, thus the speed of the controller still needs to be improved. 6 Conclusions and future work This paper proposed the hierarchical control system of the front steering angle compensation and the traction force distribution for 4WD EVs based on optimal control theory. Simulation study clarified that the vehicle yaw motion and the vehicle attitude, sideslip motion, can be much improved by the proposed hierarchical control system. The upper control loop uses a tripe-step non-linear control method to specify control dynamic which is required for keeping yaw stability. The lower control loop is used to determine control inputs for driving motors by MPC allocation method. Finally, the designed controller is evaluated through simulation model established in the AMESim environment and compared with static allocation method under the same condition. The results show that the proposed control strategy can improve the stability of 4WD EV, especially EV drives in critical conditions. 7 Acknowledgments This work was supported by the 973 Program (no. 2012CB821202), the National Nature Science Foundation of China (no , no ) and Jilin Provincial Science Foundation of China (no. 3D513U315420). 8 References 1 Mutoh, N: Driving and braking torque distribution methods for front- and rear-wheel-independent drive-type electric vehicles on roads with low friction coefficient, IEEE Trans. Ind. Electron., 2012, 59, (10), pp Njeh, M., Cauet, S., Coirault, P. et al.: H control strategy of motor torque ripple in hybrid electric vehicles: an experimental study, IET Control Theory Appl., 2011, 5, (1), pp Poussot-Vassal, C., Sename, O., Dugard, L. et al.: Vehicle dynamic stability improvements through gain-scheduled steering and braking control, Veh. Syst. Dyn., 2011, 49, (10), pp Wang, R., Zhang, H., Wang, J: Linear parameter-varying-based fault-tolerant controller design for a class of over-actuated non-linear systems with applications to electric vehicles, IET Control Theory Appl., 2014, 8, (7), pp Novellis, L.D., Sorniotti, A., Gruber, P.: Wheel torque distribution criteria for electric vehicles with torque-vectoring differentials, IEEE Trans. Veh. Technol., 2014, 63, (4), pp Tjonnas, J., Johansen, T.A.: Stabilization of automotive vehicles using active steering and adaptive brake control allocation, IEEE Trans. Control Syst. Technol., 2010, 18, pp Johansen, T.A., Fossen, T.I: Control allocation a survey, Automatica, 2013, 49, (5), pp Luo, Y., Serrani, A., Yurkovich, S. et al.: Model predictive dynamic control allocation with actuator dynamics, Proc. American Contribution Conf., Karlsruhe, Germany, 2004, vol. 2, pp Durham, W.C: Computationally efficient control allocation, AIAA J. Guid. Control Dyn., 2001, 24, (3), pp Liao, F., Lum, K.Y., Wang, J.L. et al.: Adaptive control allocation for nonlinear systems with internal dynamics, IET Control Theory Appl., 2010, 4, (6), pp Schofield, B., Hagglund, T: Optimal control allocation in vehicle dynamics control for rollover mitigation, Proc. American Contribution Conf., 2008, pp dela Pena, D.M., Alamo, T., Ramrez, D. et al.: Min-max model predictive control as a quadratic program, IET Control Theory Appl., 2007, 1, (1), pp Chen, Y., Wang, J: Fast and global optimal energy-efficient control allocation with applications to over-actuated electric ground vehicles, IEEE Trans. Control Syst. Technol., 2012, 51, (6), pp Vermillion, C., Sun, J., Butts, K: Model predictive control allocation for overactuated systems stability and performance, Proc. 46th IEEE Conf. on Decision and Control, New Orleans, LA, USA, 2007, pp Bächle, T.B., Graichen, K., Buchholz, M. et al.: Slip-constrained model predictive control allocation for an all-wheel driven electric vehicle, Proc. of the 19th IFAC World Congress, Cape Town, South Africa, 2014, pp Zhao, H., Gao, B., Ren, B.: Integrated control of in-wheel motor electric vehicles using a triple-step nonlinear method, J. Franklin Inst., 2015, 352, pp Chen, H., Allg"ower, F.: A quasi-infinite horizon nonlinear model predictive control scheme with guaranteed stability, Automatica, 1998, 34, (10), pp Chen, H.: Model predictive control, System and Control Series, (Science Press, Beijing, China, 2013, 1st edn.). 19 Doumiati, M., Sename, O., Dugard, L. et al.: Integrated vehicle dynamics control via coordination of active front steering and rear braking, Eur. J. Control, 2013, 19, (1), pp Pacejka, H.B.: Tyre and vehicle dynamics (Elsevier, London, UK, 2005, 2nd edn.). 21 Acarman, T: Nonlinear optimal integrated vehicle control using individual braking torque and steering angle with on-line control allocation by using statedependent Riccati equation technique, Dyn. Syst. Meas. Control, 2000, 122, (3), pp The Institution of Engineering and Technology 2015