#### Abstract

Two nonlinear controllers are proposed for a light-weighted all-electric vehicle: Chebyshev neural network based backstepping controller and Chebyshev neural network based optimal adaptive controller. The electric vehicle (EV) is driven by DC motor. Both the controllers use Chebyshev neural network (CNN) to estimate the unknown nonlinearities. The unknown nonlinearities arise as it is not possible to precisely model the dynamics of an EV. Mass of passengers, resistance in the armature winding of the DC motor, aerodynamic drag coefficient and rolling resistance coefficient are assumed to be varying with time. The learning algorithms are derived from Lyapunov stability analysis, so that system-tracking stability and error convergence can be assured in the closed-loop system. The control algorithms for the EV system are developed and a driving cycle test is performed to test the control performance. The effectiveness of the proposed controllers is shown through simulation results.

#### 1. Introduction

In recent years due to stricter emission standards and global fuel supply issues, researchers in the automobile industry are paying more attention to zero-polluting electric vehicles (EVs) for energy conservation and environmental protection. EV has emerged as a promising alternative to improve fuel economy while meeting the tightened emission standards [1]. EV is used in many applications, particularly for the purpose of patrolling and other short-range transportations. A lot of work has been reported in the literature for reducing the cost, extending the driving range, and sophisticated energy management strategies to improve the performance and usage of energy [2–6]. Direct current (DC) power is supplied from the battery, and therefore EVs driven by DC motor are a favorable selection. The control of DC motor is simple and it can provide comparatively larger startup torque. In addition to the primary function of propulsion, the DC motor can also be used effectively as the braking device because of its fast torque response characteristics and capability of regeneration [7, 8].

The dynamics of EV are inherently nonlinear and it is well known that nonlinear control techniques are superior to the conventional linear control methods [9, 10]. The performance of nonlinear control techniques, specifically, the differential-geometric approach, to the control of EV is presented in [11]. Due to the complex operating conditions of EV, intelligent or fuzzy control is suggested in [12, 13]. In [12], the fuzzy logic controller (FLC) is modeled to be capable of increasing the initial torque required for the electric vehicle drive with variable speed characteristics and high efficiency. In [13], fuzzy logic based controller to control the wheel slip for electric vehicle antilock braking systems (ABSs) is developed. A methodology to generate stochastic drive cycles for the design and control optimization of EVs is detailed in [14]. An algorithm for determining online energy based driving guidance for an EV using particle swarm optimization is developed in [15].

Neural networks have been used for approximation of nonlinear systems, for classification of signals and for associative memory. Chebyshev neural network (CNN) has been shown to be able to approximate any continuous functions over a compact set to arbitrary accuracy [16–18]. CNN is a functional link neural network (FLN) based on Chebyshev polynomials. The efficacy of CNN in the areas of online system identification [19, 20] and tracking controller for nonlinear systems has been established [21–23]. The proposed controller does not necessitate exact knowledge of the unknown nonlinearities. The CNN is used for estimating the unknown nonlinearities of the system. The adaptation laws for the CNN weights are such that they guarantee the stability of the system. The tracking error mainly depends on the CNN, feedback functions to be used for the weight adaptation law, and other design parameters.

It is not possible to precisely model the dynamics of an EV as some parameters may vary with time/conditions. For example, the resistance in the armature winding of the DC motor changes as the temperature varies and the aerodynamic drag coefficient and the rolling resistance coefficient are varying because of wind and road conditions, respectively. In this paper, the resistance in the armature winding () of the DC motor, the aerodynamic drag coefficient , the rolling resistance coefficient , and the mass of the passengers () are considered to be varying with time resulting in unknown nonlinearities. The aim of this paper is to design CNN based backstepping controller and CNN based optimal adaptive controller for EV in the presence of unknown nonlinearities and test the performance of the overall system on NEDC drive cycle test.

The paper is organized as follows. In Section 2, the description of the complete EV system and the structure of CNN will be presented. The problem statement is introduced in Section 3. The design of conventional backstepping controller and CNN based backstepping controller is described in Section 4. In Section 5, we give an optimal control design for EV systems using the H-J-B equation followed by a CNN based optimal adaptive controller. Section 6 validates the performance of the proposed controllers through simulations and drive cycle test followed by conclusion.

#### 2. EV System Description and CNN Structure

An EV system dynamics mainly comprises two parts: the vehicle dynamics and dynamics of the motor system, as shown in Figure 1. Motor system is connected to EV system through transmission unit, which includes the gearing system. In the actual EVs, the driver provides the command signal through the accelerator/brake pedal in the form of acceleration/deceleration to the controller of the propulsion system. The DC motor is used in the proposed EV system for propulsion and DC motor system is connected to EV system through transmission unit, which includes the gearing system. Accordingly, the speed of DC motor is controlled so as to control the actual EV system.

##### 2.1. Vehicle Dynamics

The major factors that affect the vehicle dynamics are road condition, aerodynamic drag, hill climbing, acceleration, and so forth. After these factors are taken into account, vehicle dynamics can be written as follows [1]: where is the rolling resistance coefficient, is the mass of the EV, is the gravity acceleration, is the air density, is the frontal area of the vehicle, is the drag coefficient, is the driving velocity of the vehicle, and is the hill climbing angle. In this paper, the aerodynamic drag coefficient and the rolling resistance coefficient are assumed to be varying with time. In [1] is a constant which is a very stringent assumption. In the proposed work includes the mass of vehicle and the mass of passengers ; that is, . Thus is varying with time and not a constant.

In the vehicle dynamics (1), the first term assimilates to the rolling resistance force; the second term assimilates to the aerodynamic drag force; the third term assimilates to the hill climbing force; and the fourth term assimilates to the acceleration force. This resultant traction force will produce a counterproductive torque to the driving motor, which is represented by the following relationship: where is the torque produced by the driving motor, is the tyre radius of the EV, and is the gearing ratio.

##### 2.2. Motor Dynamics

The EV is driven by a DC motor and the dynamics of which are given by [11] where is the inertia of the motor, including the gearing system and the tyres, is the motor angular speed, is the armature current (also field current), , and are the armature inductance, armature resistance, field winding inductance, and field winding resistance, respectively, is the viscous coefficient, corresponds to the external torque, is the control input voltage, and is the mutual inductance between the armature winding and the field winding, generally nonlinear because of saturation. In this paper the resistance in the armature winding of the DC motor is considered to be varying as the armature winding resistance of the DC motor changes as the temperature varies.

##### 2.3. Complete Dynamics and Open-Loop Response

With vehicle dynamics (1) and motor dynamics (3), the complete dynamics of the EV system becomes [11]

By using (4), the open-loop response of a light-weighted all-electric vehicle is studied. The parameters used on a light-weighted all-electric vehicle are specified in Table 1. The simulation result of open-loop response is shown in Figure 2. The plot shows the full power speed characteristics. As given in Table 1 the desired nominal speed is = 25 km/hr whereas in open-loop conditions the speed is beyond 40 km/hr which is not acceptable. Therefore, it is required to design proper controller.

The relation between the driving velocity of the vehicle and the motor angular speed is given as where is the tyre radius of the EV and is the gearing ratio.

##### 2.4. Structure of Neural Network

In this paper, a single layer CNN is considered for the NN structure. CNN consists of a functional expansion (FE) block and a single-layer perceptron network. The purpose of the FE block is to increase the dimension of the input pattern so as to improve representation of the input pattern in a higher dimensional space. Chebyshev expansions are frequently used for approximations to functions as they are much more efficient than other power series expansion of the same degree. Among orthogonal polynomials, the Chebyshev polynomials which are derived from the solution of the Chebyshev differential equation occupy an important place, since, in the case of a broad class of functions, expansions in Chebyshev polynomials converge more rapidly than expansions in other set of polynomials. Hence, we consider the Chebyshev polynomials as basis functions for the neural network.

The Chebyshev polynomials can be generated by the following recursive formula [17]: where is a Chebyshev polynomial, is the order of Chebyshev polynomials chosen, and here is scalar quantity. can be chosen as , , or . In this paper is chosen as . For example, an enhanced pattern using the Chebyshev polynomials for is obtained as where is a Chebyshev polynomial, is the order of the selected Chebyshev polynomial, and , . denotes the Chebyshev polynomial basis function.

Referring to Figure 3, the architecture of the CNN consists of two parts [17], namely, numerical transformation part and learning part. The numerical transformation is the FE of the input pattern consisting of a finite set of Chebyshev polynomials. Consequently, the Chebyshev polynomial basis can be considered as a new input vector. The learning part involves functional-link neural network based on Chebyshev polynomials. The CNN is a single-layered neural network, and, in general, its learning is fast [16, 17].

On the basis of approximation property of CNN [19], a general nonlinear function can be approximated by CNN as where is the CNN functional reconstruction error vector and which is bounded, is the optimal weight matrix, and denotes the Chebyshev polynomial basis function. The output of the CNN is given by where is the estimate of the optimal weight matrix .

#### 3. Problem Statement

The complete dynamics in (4) can be described as where Assuming the hill climbing angle to be zero, the nonlinear functions and become where , , , , , , , , and .

As a result, the EV system becomes

The aim of this paper is to compare and test the performance of CNN based backstepping controller with a CNN based optimal adaptive controller, which forces the plant output to track a specified reference trajectory in the presence of time varying mass “” and varying armature winding resistance (), aerodynamic drag coefficient , and the rolling resistance coefficient ; that is,

#### 4. Nonlinear Backstepping Controller Design

The nonlinear backstepping controller is designed in the following ways.

##### 4.1. Conventional Backstepping Controller

In this subsection, the steps involved in the development of backstepping controller for system (14) are discussed.

By selecting the following state transformation: the system (14) becomes which is in the strict feedback form [24, 25]. In this case, all the nonlinear functions are considered known. The error dynamics is defined as The time derivative of can be found as From (17) and (19), we have (20) as By selecting , (21) yields where is the positive constant.

Differentiating (19), Substituting from (17), (23) becomes To stabilize (24), the control effort is chosen as where is the positive constant.

Substituting (25) in (24) gives

To prove the convergence of the EV system a Lyapunov function is chosen as The time derivative of (27) is Using (22) and (26) we get Thus, it can be easily seen that the system is globally asymptotically stable.

##### 4.2. CNN Based Backstepping Controller

The EV system in (17) is rewritten in the following form for the implementation of the CNN based backstepping controller: Here the nonlinear functions and are unknown as , , , and are varying with time.

The error dynamics is defined as The CNN based backstepping controller is designed in the following two steps.

*Step 1.* The time derivative of is given by
Using (30) and (31), (32) becomes
Stabilize (33) by choosing
where is the positive constant.

Substituting (34) in (33), becomes

*Step 2.* The time derivative of can be obtained as
By using (30), (36) results in

The unknown nonlinear functions and will be estimated by CNN 1 and CNN 2, respectively. The nonlinear functions and can be represented by a CNN as where and are the bounded CNN approximation errors, and are the optimal weight matrices, and and are the basis functions.

The estimate of and of can be written as where is the estimate of the and is the estimate of the .

Adding and subtracting and rearranging, (37) we have To stabilize (40), the control effort is chosen as where is the positive constant. The block diagram of the overall system is presented in Figure 4.

Substituting (41) in (40), becomes

Define the estimation error as Using (38) and (39) in (43) gives where and are weight errors.

Now by applying (44) in (42) becomes

Two standard assumptions, which are commonly used in the neural networks literature, are given below [26].

*Assumption 1. *The optimal weights and are bounded by known positive values so that
We only need to know that ideal weights exist to prove the convergence analysis. The exact value of the ideal weights need not be known. The symbol denotes the Frobenius norm; that is, given a matrix , the Frobenius norm is defined by

*Assumption 2. *The desired trajectory and its derivatives up to second order are bounded.

Based on the above Assumptions 1 and 2, the stability analysis is given in Section 4.3.

##### 4.3. Stability Analysis

Theorem 3. *Consider the EV system (30) and control input (41) satisfying Assumptions 1 and 2. If the weights of the CNN 1 and CNN 2 are updated according to adaptation law given in (48) and (49), respectively,
**
where and are the learning rate and is damping coefficient, then the weight errors , and the errors and are uniformly ultimately bounded (UUB).*

*Proof. *Consider the Lyapunov function
The time derivative of (50) is
Now substitute and from (35) and (45), respectively, and perform a simple manipulation (i.e., using , for placing weight matrices inside a trace operator). Then we have
With the adaptation law given in (48) and (49), (52) becomes

Apply the following inequality [27]:
And assume that the upper bounds are as follows:
Now we can express (53) as
As the first term is always negative, now defining nonnegative term and completing the square terms in (56) yield
which is negative as long as either (58), (59), or (60) holds
Or

Thus, is negative outside a compact set. According to a standard Lyapunov theorem extension [28], this demonstrates uniform ultimate boundedness of weight errors , and errors , .

#### 5. Nonlinear Optimal Controller Design Using CNN

The following approach is used to design the nonlinear optimal controller for the EV system (17).

The tracking errors are defined as And the filtered tracking error is defined as where is the positive constant.

Differentiating (62) and rearranging it, we have Substituting from EV system (17), (63) becomes where −−− and .

System (64) may be written as where

Now we define an auxiliary control input , which is to be optimized in the next subsection as with as the control input. The closed-loop system becomes

##### 5.1. Optimal Controller Using Hamilton-Jacobi-Bellman (H-J-B) Optimization

The augmented system [27] is achieved using (62) and (68) or with shorter notation with being defined as , , and . A quadratic performance measure is as follows: with the Lagrangian

The objective is to find the auxiliary control input that minimizes the quadratic performance measure subject to the constraints imposed by (69) which will be denoted by in the presence of known nonlinearities.

A necessary and sufficient condition for is that there exist a function , which is accredited as the value function and satisfies the H-J-B equation [29] where the Hamiltonian of optimization is described as and satisfies the partial differential equation The minimum is achieved for , and the Hamiltonian is then given by

Lemma 4 (see [27]). *The function composed of and satisfies the H-J-B equation
**
where and in (62) and (77) respectively can be found from the Riccati differential equation
**
The optimal control that minimizes (71) subject to (70) is
*

Let , be chosen such that with [27]. Then the and required in Lemma 4 can be calculated as given below with (82) solved for .

Using (67) and (79), the input to the EV system is given as where is given by (66) and is assumed to be known. The following subsection details the stability analysis.

##### 5.2. Stability Analysis

Suppose that and exist that satisfy Lemma 4, and in addition, there exist constants and such that , and the spectrum of is bounded in the sense that on . Then using the feedback control (79) into (70) results in the controlled nonlinear system becomes

The function is chosen as a suitable Lyapunov function; it shows that for all . The time derivative of is given by From the solution of the H-J-B equation (75) and using (85), it results that Substituting optimal control law (79) into (72), the time derivative of becomes

The time derivative of the Lyapunov function is negative definite, implying global exponential stability.

##### 5.3. CNN Based Optimal Adaptive Controller

In Section 5.1, the nonlinear function (66) is assumed to be known. This assumption is relaxed and is treated as an unknown nonlinear function. The function is estimated using CNN neural network. The nonlinear function can be represented by a CNN as where is a basis function for the CNN. The block diagram in Figure 5 demonstrates the CNN neural controller based on H-J-B optimization.

The estimate of can be written as Using (67), (79), and (89), the input is given by where is a robustifying term, which is given by with and being defined as the filtered tracking error in (62). Using (88) and (90), (65) becomes where is the weight-estimation error. Using (92) in (70) yields with , , and being given in (69) and (70).

Using the optimal control law (79) into (93), we have

Theorem 5. *Suppose the optimal control law given by (79) minimizes the quadratic performance measure given in (71). If the weights of the CNN are updated according to adaptive learning law given by
**
with and , then the errors , , and are uniformly ultimately bounded.*

*Proof. *Consider the following Lyapunov function:
where is positive given by (81). The time derivative of (96) becomes
Substituting from (94), (97) becomes
Using and from the Riccati equation (78), we obtain
Then applying (99) in (98) and performing a simple manipulation for placing weight matrices inside a trace operator, we have
Now substitute the robustifying term (91), the adaptive learning law (95), and the following inequality:
The time derivative becomes
Completing the square terms yields
which is guaranteed to be negative as long as either (104) or (105) holds

Thus, is negative outside a compact set. According to a standard Lyapunov theory extension [28], this demonstrates uniform ultimate boundedness of , , and .

#### 6. Simulation Results

The drive cycle tests that are currently used for light-weighted EVs are new European driving cycle (NEDC), Federal Test Procedure (FTP-75), and JC08. The NEDC is used in Europe, and the low powered EV version of this cycle is used in India. The FTP 75 cycle is used in USA and the JC08 in Japan. In order to show the validity of the proposed controllers, the NEDC is used for testing the performance. The NEDC is a driving cycle consisting of four repeated ECE-15 driving cycles and an extra-urban driving cycle (EUDC) [1]. The maximum speed of NEDC is 120 km/h but it is scaled to 50 km/h when applied in this paper [11].

The simulation is implemented in MATLAB 7.8.0 (2009a) with m-file programming. The controller design parameters for conventional backstepping are chosen as = 15 and = 15. The controller design parameters for CNN based backstepping are chosen as = 0.8 and = 0.35. For update of parameters in (48) and (49), , , and are chosen as 1, 0.0001, and 0.15, respectively. For approximating the input to the CNN 1 is , , and for approximating the input to the CNN 2 is , . The order of Chebyshev polynomial is chosen as 1 for both and . The parameters of CNN are initialized to zero. Thus and have dimension (5 × 1), where = = . and have dimension of (5 × 1). For optimal controller the design parameters are chosen as = 1.25, = 4. For CNN based optimal adaptive controller parameter is chosen as 0.0001. For update of parameters in (95), and are chosen as 0.01. For approximating the input to the CNN is , , , , , and . The order of Chebyshev polynomial is chosen as 1 for all inputs to the CNN. The parameters of CNN are initialized to zero. Thus and have dimension of (13 × 1). The initial conditions for = .

Performance of designed controllers for mass variation as given in (106) is considered. Passengers mass is increased/decreased at different point of time in the driving cycle.The variation in armature winding resistance of the DC motor due to temperature changes, the variation in the aerodynamic drag coefficient, and the variation in the rolling resistance coefficient are considered as given in (107), (108), and (109), respectively. The variation in mass, armature winding resistance, aerodynamic drag coefficient, and rolling resistance coefficient are assumed to be known in conventional backstepping and optimal controller. The variation in mass, armature winding resistance, aerodynamic drag coefficient, and rolling resistance coefficient are considered unknown in CNN based backstepping controller and CNN based optimal adaptive controller.

Choosing

The drive cycle test performances and tracking errors for conventional backstepping controller and CNN based backstepping controller are shown in Figures 6 and 7, respectively. It is clear that the conventional backstepping controller has better tracking performance than the CNN based backstepping controller, in the range of speed below designed nominal speed ( = 25 km/hr). The CNN based backstepping controller performs much better in high speed range. The amp-hour consumption for CNN based backstepping controller is 4.4834 km/10.75 AH.

**(a)**

**(b)**

**(a)**

**(b)**

The drive cycle test performances and tracking errors for optimal controller and CNN based optimal adaptive controller are shown in Figures 8 and 9, respectively. The amp-hour consumption for CNN based optimal adaptive controller is 4.4795 km/10.70 AH. It can be seen that the optimal controller has better tracking performance in the speed range below designed nominal speed than conventional backstepping controller and CNN based backstepping controller. However, for the high speed range the performance of CNN based backstepping controller is comparable to the optimal controller. CNN based optimal adaptive controller gives the best tracking performance on the entire speed range of the drive cycle test as compared to CNN based backstepping controller designed in current work and the results presented in [11] which are reproduced in Figure 10. The comparative result of the driving cycle test is presented in Table 2 and to provide detailed quantitative analysis of the designed controllers, root mean squared (RMS) tracking error of controllers is presented in Table 3. The proposed CNN based optimal adaptive learning shows both robustness and adaptation to changing system dynamics and unknown nonlinearities.

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

The traction force is required to propel the EV in forward direction. It is produced by the DC motor torque and transferred through transmission unit, which includes the gearing system and finally drive the vehicle. While the vehicle is in motion, there are forces that try to stop its movement. These forces usually include rolling resistance and aerodynamic drag force. The simulation result of traction force of the EV is shown in Figure 11.

In present work, all the parameters are specified on a light-weighted all-electric vehicle and are given in Table 1. The mass of the vehicle is considered as 800 kg. This includes the body mass of the vehicle, accessories, weight of the DC motor, and batteries. For modeling light-weighted EVs, we have considered four batteries each of 12 volts and this restricts the control effort in the range of 0 ~ 48 V (see Table 1). The control effort is shown in Figure 12.

#### 7. Conclusion

Two nonlinear controllers for a light-weighted all-electric vehicle, CNN based backstepping controller and CNN based optimal adaptive controller, are presented in this paper. The unknown nonlinearities in EV system arise due to varying mass of passengers, varying resistance in the armature winding of the DC motor and variation in aerodynamic drag coefficient and the rolling resistance coefficient are estimated by CNN. The CNN weights are updated online according to the adaptive-learning algorithm, which is obtained from Lyapunov stability analysis, so that system-tracking stability and error convergence can be assured in the closed-loop system. The salient feature of the proposed design methodologies demonstrates that the control objective is obtained with unknown nonlinear dynamics of the EV system. The NEDC is used for testing the performance of the proposed controllers. It is shown that the tracking performance of the controllers designed in this paper is satisfactory in both the cases. The test results for CNN based backstepping controller and the CNN based optimal adaptive controller have better tracking performance than that reported by Huang et al. [11] and amp-hour consumption is also less than the nonlinear controllers described in [11]. From simulation results and the comparative and quantitative results presented in Tables 2 and 3, respectively, it is clear that the CNN based optimal adaptive controller gives better performance as compared to other controllers.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

The authors would like to thank the editor and the anonymous reviewers for their valuable comments and constructive suggestions which have helped in improving the quality of the paper.