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 Cd and the rolling resistance coefficient μrr are varying because of wind and road conditions, respectively. In this paper, the resistance in the armature winding (Ra) of the DC motor, the aerodynamic drag coefficient Cd, the rolling resistance coefficient μrr, and the mass of the passengers (ΔM) 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.
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]:
(1)F=μrrmg︸1+12ρACdv2︸2+mgsinϕ︸3+mdvdt︸4,
where μrr is the rolling resistance coefficient, m is the mass of the EV, g is the gravity acceleration, ρ is the air density, A is the frontal area of the vehicle, Cd is the drag coefficient, v is the driving velocity of the vehicle, and ϕ is the hill climbing angle. In this paper, the aerodynamic drag coefficient Cd and the rolling resistance coefficient μrr are assumed to be varying with time. In [1] m is a constant which is a very stringent assumption. In the proposed work m includes the mass of vehicle M and the mass of passengers ΔM; that is, m=M+ΔM. Thus m 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 F will produce a counterproductive torque to the driving motor, which is represented by the following relationship:
(2)TL=FrG,
where TL is the torque produced by the driving motor, r is the tyre radius of the EV, and G 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]
(3)Jdωdt=Lafi2-Bω-TL,(La+Lfield)didt=V-(Ra+Rf)i-Lafiω,
where J is the inertia of the motor, including the gearing system and the tyres, ω is the motor angular speed, i is the armature current (also field current), La,Ra,Lfield, and Rf are the armature inductance, armature resistance, field winding inductance, and field winding resistance, respectively, B is the viscous coefficient, TL corresponds to the external torque, V is the control input voltage, and Laf 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 Ra 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]
(4){J+mr2G2}dωdt=Lafi2-Bω-rG{μrrmg+12ρACdv2+mgsinϕ}(La+Lfield)didt=V-(Ra+Rf)i-Lafiω.
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 v = 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.
Parameters of the EV system [11].
Motor
Vehicle
La+Lfield (mH)
6.008
M (kg)
800
Ra+Rf (Ω)
0.12
A (m2)
1.8
B (N.M.s.)
0.0002
ρ (kg/m3)
1.25
J (kg m2)
0.05
Cd
0.3
Laf (mH)
1.766
ϕ (°)
0
V (V)
0~48
μrr
0.015
i (A)
78 A (250 max)
r (m)
0.25
ω (r/min)
2800 (v = 25 km/h)
G
11
ΔM = 150 (kg) for t ≤ 195 and t≥ 780 (kg), ΔM = 220 for t > 195, and t≤ 585 and ΔM = 300 (kg) for t > 585 and t < 780.
Open-loop response of the EV system.
The relation between the driving velocity of the vehicle v and the motor angular speed ω is given as
(5)v=rGω,
where r is the tyre radius of the EV and G 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]:
(6)Ti+1(x-)=2x-Ti(x-)-Ti-1(x-),T0(x-)=1,
where Ti(x-) is a Chebyshev polynomial, i is the order of Chebyshev polynomials chosen, and here x- is scalar quantity. T1(x-) can be chosen as x-, 2x-, 2x--1 or 2x-+1. In this paper T1(x-) is chosen as x-. For example, an enhanced pattern using the Chebyshev polynomials for x-∈ℜ2 is obtained as
(7)ϕ-(x-)=[1T1(x-1)T2(x-1)⋯T1(x-2)T2(x-2)⋯]T,
where Ti(x-j) is a Chebyshev polynomial, i is the order of the selected Chebyshev polynomial, and j=1, 2. ϕ-(x-) 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].
Structure of CNN.
On the basis of approximation property of CNN [19], a general nonlinear function y-(x-) can be approximated by CNN as
(8)y-(x-)=W-Tϕ-+ε-,
where ε- is the CNN functional reconstruction error vector and ∥ε-∥≤ε-N which is bounded, W- is the optimal weight matrix, and ϕ- denotes the Chebyshev polynomial basis function. The output of the CNN is given by
(9)y-^=W-^Tϕ-,
where W-^ is the estimate of the optimal weight matrix W-.
3. Problem Statement
The complete dynamics in (4) can be described as
(10)X˙=f(X)+g(X)u,y=h(X),
where
(11)X=[x1x2]=[ωi],f(X)=[-Ra+RfLa+Lfieldx2-LafLa+Lfieldx1x21J+m(r2G2)]1J+m(r2/G2){(μrrmg+12ρACdr2G2x12+mgsinϕ)Lafx22-Bx1-rG×(μrrmg+12ρACdr2G2x12+mgsinϕ)}-Ra+RfLa+Lfieldx2-LafLa+Lfieldx1x21J+m(r2G2)],(12)g(X)=[01La+Lfield],h(X)=x1.
Assuming the hill climbing angle ϕ to be zero, the nonlinear functions f(X) and g(X) become
(13)f(X)=[K1m+K2{K3x22-K4x1-K5x12-K6m}-K7x2-K8x1x2],g(X)=[0K9],h(X)=x1,
where K1=G2/r2, K2=(G2/r2)J, K3=Laf, K4=B, K5=(1/2)ρACd(r3/G3), K6=(r/G)μrrg, K7=(Ra+Rf)/(La+Lfield), K8=Laf/(La+Lfield), and K9=1/(La+Lfield).
As a result, the EV system becomes
(14)x˙1=K1m+K2{K3x22-K4x1-K5x12-K6m},x˙2=-K7x2-K8x1x2+K9u,y=x1.
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 y to track a specified reference trajectory yd in the presence of time varying mass “m” and varying armature winding resistance (Ra), aerodynamic drag coefficient Cd, and the rolling resistance coefficient μrr; that is,
(15)limt→∞(y-yd)=0.
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:
(16)z1=x1,z2=x˙1,
the system (14) becomes
(17)z˙1=z2,z˙2=-2K1K3K7m+K2x22-2K1K3K8m+K2z1x22-K1K4m+K2z2-2K1K5m+K2z1z2+2K1K3K9m+K2x2u,y=z1,
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
(18)e1=z1-z1d=z1-yd(19)e2=z2-z2d.
The time derivative of e1 can be found as
(20)e˙1=z˙1-y˙d.
From (17) and (19), we have (20) as
(21)e˙1=e2+z2d-y˙d.
By selecting z2d=y˙d-k-1e1, (21) yields
(22)e˙1=e2-k-1e1,
where k-1 is the positive constant.
Differentiating (19),
(23)e˙2=z˙2-z˙2d.
Substituting z˙2 from (17), (23) becomes
(24)e˙2=-2K1K3K7m+K2x22-2K1K3K8m+K2z1x22-K1K4m+K2z2-2K1K5m+K2z1z2+2K1K3K9m+K2x2u-z˙2d.
To stabilize (24), the control effort u is chosen as
(25)u={z˙2d-e1+2K1K3K7m+K2x22+2K1K3K8m+K2z1x22+K1K4m+K2z2+2K1K5m+K2z1z2-k-2e2}×1{(2K1K3K9)/(m+K2)}x2,
where k-2is the positive constant.
Substituting (25) in (24) gives
(26)e˙2=-e1-k-2e2.
To prove the convergence of the EV system a Lyapunov function is chosen as
(27)L1=12e12+12e22.
The time derivative of (27) is
(28)L˙1=e1e˙1+e2e˙2.
Using (22) and (26) we get
(29)L˙1=-k-1e12-k-2e22.
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:
(30)z˙1=z2,z˙2=f(x,z)+g(x,z)u,y=z1.
Here the nonlinear functions f(x,z) and g(x,z) are unknown as m, Ra, Cd, and μrr are varying with time.
The error dynamics is defined as
(31)e1=z1-z1d=z1-yde2=z2-z2d.
The CNN based backstepping controller is designed in the following two steps.
Step 1. The time derivative of e1 is given by
(32)e˙1=z˙1-y˙d.
Using (30) and (31), (32) becomes
(33)e˙1=e2+z2d-y˙d.
Stabilize (33) by choosing
(34)z2d=y˙d-k1e1,
where k1is the positive constant.
Substituting (34) in (33), e˙1 becomes
(35)e˙1=e2-k1e1.
Step 2. The time derivative of e2 can be obtained as
(36)e˙2=z˙2-z˙2d.
By using (30), (36) results in
(37)e˙2=f(x,z)+g(x,z)u-z˙2d.
The unknown nonlinear functions f(x,z) and g(x,z) will be estimated by CNN 1 and CNN 2, respectively. The nonlinear functions f(x,z) and g(x,z) can be represented by a CNN as
(38)f=W1Tϕ1+ε1,g=W2Tϕ2+ε2,
where ε1 and ε2 are the bounded CNN approximation errors, W1 and W2 are the optimal weight matrices, and ϕ1 and ϕ2 are the basis functions.
The estimate f^ of f and g^ of g can be written as
(39)f^=W^1Tϕ1,g^=W^2Tϕ2,
where W^1 is the estimate of the W1 and W^2 is the estimate of the W2.
Adding and subtracting g^(x,z)u and rearranging, (37) we have
(40)e˙2=f(x,z)+{g(x,z)-g^(x,z)}u-z˙2d+g^(x,z)u.
To stabilize (40), the control effort u is chosen as
(41)u=1g^(x,z){-f^(x,z)+z˙2d-k2e2-e1},
where k2 is the positive constant. The block diagram of the overall system is presented in Figure 4.
Block diagram of CNN based backstepping controller.
Substituting (41) in (40), e˙2 becomes
(42)e˙2={f(x,z)-f^(x,z)}+{g(x,z)-g^(x,z)}u-k2e2-e1.
Define the estimation error as
(43)f~=f-f^,g~=g-g^.
Using (38) and (39) in (43) gives
(44)f~=W~1Tϕ1+ε1,g~=W~2Tϕ2+ε2,
where W~1=W1-W^1 and W~2=W2-W^2 are weight errors.
Now by applying (44) in (42) e˙2 becomes
(45)e˙2=W~1Tϕ1+ε1+(W~2Tϕ2+ε2)u-k2e2-e1.
Two standard assumptions, which are commonly used in the neural networks literature, are given below [26].
Assumption 1.
The optimal weights W1 and W2 are bounded by known positive values so that
(46)∥W1∥F≤W1M,∥W2∥F≤W2M.
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 ∥•∥F denotes the Frobenius norm; that is, given a matrix A, the Frobenius norm is defined by
(47)∥A∥F2=tr(ATA)=∑i,jaij2.
Assumption 2.
The desired trajectory yd 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 AnalysisTheorem 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,
(48)W^˙1=η1ϕ1e2-ρ-η1∥e2∥W^1,(49)W^˙2=η2ϕ2ue2-ρ-η2∥e2∥|u|W^2,
where η1 and η2 are the learning rate and ρ- is damping coefficient, then the weight errors W~1=W1-W^1, W~2=W2-W^2 and the errors e1 and e2 are uniformly ultimately bounded (UUB).
Proof.
Consider the Lyapunov function
(50)L2=12e12+12e22+12tr{W~1Tη1-1W~1}+12tr{W~2Tη2-1W~2}.
The time derivative of (50) is
(51)L˙2=e1e˙1+e2e˙2+tr{W~1Tη1-1W~˙1}+tr{W~2Tη2-1W~˙2}.
Now substitute e˙1 and e˙2 from (35) and (45), respectively, and perform a simple manipulation (i.e., using xTy=tr{xTy}=tr{yxT}, for placing weight matrices inside a trace operator). Then we have
(52)L˙2=-k1e12+e2(ε1+ε2u-k2e2)+tr{W~1T(ϕ1e2+η1-1W~˙1)}+tr{W~2T(ϕ2e2+η2-1W~˙2)}.
With the adaptation law given in (48) and (49), (52) becomes
(53)L˙2=-k1e12+e2(ε1+ε2u-k2e2)+tr{W~1T(ρ-∥e2∥W^1)}+tr{W~2T(ρ-∥e2∥|u|W^2)}.
Apply the following inequality [27]:
(54)tr[W~T(W-W~)]=〈W~,W〉F-∥W~∥F2≤∥W~∥FWM-∥W~∥F2.
And assume that the upper bounds are as follows:
(55)∥ε1∥≤ε1M,∥ε2∥≤ε2M,|u|≤uS.
Now we can express (53) as
(56)L˙2≤-k1∥e1∥2+∥e2∥(ε1M+ε2MuS-k2∥e2∥)+ρ-∥e2∥{(∥W~1∥FW1M-∥W~1∥F2)}+ρ-∥e2∥uS{(∥W~2∥FW2M-∥W~2∥F2)}.
As the first term is always negative, now defining nonnegative term and completing the square terms in (56) yield
(57)L˙2≤-∥e2∥[k2∥e2∥-ε1M-ε2MuS(∥W~2∥F2-W2M2)2+ρ-(∥W~1∥F2-W1M2)2-ρ-W1M24+ρ-uS(∥W~2∥F2-W2M2)2-ρ-uSW2M24(∥W~2∥F2-W2M2)2]
which is negative as long as either (58), (59), or (60) holds
(58)[ε1M+ε2MuS+ρ-(W1M2/4+uSW2M2/4)]k2<∥e2∥,(59)W1M2+(W1M24)+ε1Mρ-≤∥W~1∥F.
Or
(60)W2M2+uS(W2M24+ε2Mρ-)≤∥W~2∥F.
Thus, L˙2 is negative outside a compact set. According to a standard Lyapunov theorem extension [28], this demonstrates uniform ultimate boundedness of weight errors W~1, W~2 and errors e1, e2.
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
(61)e(t)=z1(t)-yd(t),e˙(t)=z˙1(t)-y˙d(t)=z2(t)-y˙d(t).
And the filtered tracking error is defined as
(62)r(t)=e˙(t)+Λe(t),
where Λ is the positive constant.
Differentiating (62) and rearranging it, we have
(63)r˙(t)=Λr(t)-y¨d(t)+z˙2(t)-Λ2e(t).
Substituting z˙2(t) from EV system (17), (63) becomes
(64)r˙(t)=Λr(t)-y¨d(t)+X1x+F1xu-Λ2e(t),
where X1x=-(2K1K3K7/(m+K2))x22−(2K1K3K8/(m+K2))z1x22−((K1K4)/(m+K2))z2−((2K1K5)/(m+K2))z1z2 and F1x=((2K1K3K9)/(m+K2))x2.
System (64) may be written as
(65)r˙(t)=Λr(t)+F1xu(t)+h(x),
where
(66)h(x)=X1x-y¨d(t)-Λ2e(t).
Now we define an auxiliary control input u-(t), which is to be optimized in the next subsection as
(67)u-(t)=h(x)+F1xu(t)
with u(t) as the control input. The closed-loop system becomes
(68)r˙(t)=Λr(t)+u-(t).
5.1. Optimal Controller Using Hamilton-Jacobi-Bellman (H-J-B) Optimization
The augmented system [27] is achieved using (62) and (68)
(69)[e˙r˙]=[-ΛI0Λ][er]+[0I]u-(t)
or with shorter notation
(70)z~˙(t)=Az~(t)+Bu-(t)
with z~(t) being defined as z~(t)=[e(t)r(t)]T, A=[-ΛI;0Λ], and B=[0I]T. A quadratic performance measure J(u-) is as follows:
(71)J(u-)=∫t0∞L(z~,u-)dt
with the Lagrangian
(72)L(z~,u-)=12z~T(t)Qz~(t)+12u-T(t)Ru-(t)=12[er][Q11Q12Q12TQ22][er]+12u-TRu-.
The objective is to find the auxiliary control input u-(t) that minimizes the quadratic performance measure J(u-) subject to the constraints imposed by (69) which will be denoted by u-*(t) in the presence of known nonlinearities.
A necessary and sufficient condition for u-*(t) is that there exist a function V=V(z~,t), which is accredited as the value function and satisfies the H-J-B equation [29]
(73)∂V(z~,t)∂t+minu-[H(z~,u-,∂V(z~,t)∂t,t)]=0,
where the Hamiltonian of optimization is described as
(74)H(z~,u-,∂V(z~,t)∂t,t)=L(z~,u-)+∂V(z~,t)∂tz~˙
and V(z~,t) satisfies the partial differential equation
(75)-∂V(z~,t)∂t=L(z~,u-*)+∂V(z~,t)∂tz~˙.
The minimum is achieved for u-(t)=u-*(t), and the Hamiltonian is then given by
(76)H*=minu-[L(z~,u-)+∂V(z~,t)∂tz~˙]=H(z~,u-*,∂V(z~,t)∂t,t)=-∂V(z~,t)∂t.
Lemma 4 (see [27]).
The function V composed of z~ and K satisfies the H-J-B equation
(77)V=12z~TPz~=12z~T[K00I]z~,
where Λ and K in (62) and (77) respectively can be found from the Riccati differential equation
(78)PA+ATPT-PBR-1BTP+P˙+Q=0.
The optimal control u-*(t) that minimizes (71) subject to (70) is
(79)u-*(t)=-R-1BTPz~=-R-1r(t).
Let Q, R be chosen such that
(80)Q=[Q11Q12Q12TQ22]>0,R-1=Q22
with Q12+Q12T<0 [27]. Then the Λ and K required in Lemma 4 can be calculated as given below
(81)K=-12(Q12+Q12T)>0,(82)ΛTK+KΛ=Q11
with (82) solved for Λ.
Using (67) and (79), the input to the EV system is given as
(83)u*(t)=(inv(Fx1)){u-*(t)-h(x)},
where h(x) is given by (66) and is assumed to be known. The following subsection details the stability analysis.
5.2. Stability Analysis
Suppose that K and Λ exist that satisfy Lemma 4, and in addition, there exist constants k1′ and k2′ such that 0<k1′<k2′<∞, and the spectrum of P is bounded in the sense that k1′I<P<k2′I on (t0,∞). Then using the feedback control (79) into (70) results in the controlled nonlinear system becomes
(84)z~˙(t)={A-BR-1BTP}z~(t).
The function V(z~,t) is chosen as a suitable Lyapunov function; it shows that dV/dt<0 for all ∥z~∥≠0. The time derivative of V(z~,t) is given by
(85)dV(z~,t)dt=∂V(z~,t)∂t+∂V(z~,t)∂tz~˙.
From the solution of the H-J-B equation (75) and using (85), it results that
(86)dV(z~,t)dt=-L(z~,u-*).
Substituting optimal control law (79) into (72), the time derivative of V(z~,t) becomes
(87)dV(z~,t)dt=-12{z~TQz~+(BTPz~)TR-1(BTPz~)}<0just∀t>0,z~≠0.
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 h(x) is treated as an unknown nonlinear function. The function h(x) is estimated using CNN neural network. The nonlinear function h(x) can be represented by a CNN as
(88)h(x)=WTφ(x)+ε(x)∥ε(x)∥≤εM,
where φ(x) is a basis function for the CNN. The block diagram in Figure 5 demonstrates the CNN neural controller based on H-J-B optimization.
Block diagram of CNN based optimal adaptive controller.
The estimate h^(x) of h(x) can be written as
(89)h^(x)=W^Tφ(x).
Using (67), (79), and (89), the input is given by
(90)u(t)=(inv(Fx1)){u-*(t)-W^Tφ(x)-ν(t)},
where ν(t) is a robustifying term, which is given by
(91)ν(t)=-kzr(t)∥r(t)∥
with kz≤bd and r(t) being defined as the filtered tracking error in (62). Using (88) and (90), (65) becomes
(92)r˙(t)=Λr(t)+u-*(t)+W~Tφ(x)+ε(x)-ν(t),
where W~=W-W^ is the weight-estimation error. Using (92) in (70) yields
(93)z~˙(t)=Az~(t)+B[u-*(t)+W~Tφ(x)+ε(x)-ν(t)]
with A, B, and z~ being given in (69) and (70).
Using the optimal control law (79) into (93), we have
(94)z~˙(t)=(A-BR-1BTP)z~(t)+B{W~Tφ(x)+ε(x)-ν(t)}.
Theorem 5.
Suppose the optimal control law u-*(t) given by (79) minimizes the quadratic performance measure J(u-) given in (71). If the weights of the CNN are updated according to adaptive learning law given by
(95)W^˙=φ(x)z~TPBΓ-k∥z~∥W^
with Γ>0 and k>0, then the errors e(t), r(t), and W~(t) are uniformly ultimately bounded.
Proof.
Consider the following Lyapunov function:
(96)L3=12z~T[K001]z~+12tr(W~TΓ-1W~),
where K is positive given by (81). The time derivative L˙3 of (96) becomes
(97)L˙3=z~TPz~˙+12z~TP˙z~+tr(W~TΓ-1W~˙).
Substituting z~˙ from (94), (97) becomes
(98)L˙3=z~TPAz~-z~TPBR-1BTPz~+12z~TP˙z~+z~TPB{W~Tφ(x)+ε(x)-ν(t)}+tr(W~TΓ-1W~˙).
Using z~TPAz~=(1/2)z~T{ATP+PA}z~ and from the Riccati equation (78), we obtain
(99)12ATP+12PA+12P˙=-12Q+12PBR-1BTP.
Then applying (99) in (98) and performing a simple manipulation for placing weight matrices inside a trace operator, we have
(100)L˙3=-12z~TQz~-12z~TPBR-1BTPz~+z~TPB{ε(x)-ν(t)}+tr{W~T(Γ-1W~˙+φz~TPB)}.
Now substitute the robustifying term (91), the adaptive learning law (95), and the following inequality:
(101)tr[W~T(W-W~)]=〈W~,W〉F-∥W~∥F2≤∥W~∥FWM-∥W~∥F2.
The time derivative L˙3 becomes
(102)L˙3≤-12∥z~∥2{λmin(Q)+λmin(R-1)}+∥z~∥εM+k∥z~∥(∥W~∥FWM-∥W~∥F2).
Completing the square terms yields
(103)L˙3≤-12∥z~∥[(∥W~∥F-12WM)2∥z~∥{λmin(Q)+λmin(R-1)}+k(∥W~∥F-12WM)2-εM-14kWM2]
which is guaranteed to be negative as long as either (104) or (105) holds
(104)(εM+(1/4)kWM2){λmin(Q)+λmin(R-1)}≤∥z~∥,(105)εMk+14WM2+12WM≤∥W~∥F.
Thus, L˙3 is negative outside a compact set. According to a standard Lyapunov theory extension [28], this demonstrates uniform ultimate boundedness of e(t), r(t), and W~(t).
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 k-1= 15 and k-2 = 15. The controller design parameters for CNN based backstepping are chosen as k1 = 0.8 and k2 = 0.35. For update of parameters in (48) and (49), η1, η2, and ρ- are chosen as 1, 0.0001, and 0.15, respectively. For approximating f(x,z) the input to the CNN 1 is z1, z2, and for approximating g(x,z) the input to the CNN 2 is z1, z2. The order of Chebyshev polynomial is chosen as 1 for both z1 and z2. The parameters of CNN are initialized to zero. Thus ϕ1 and ϕ2 have dimension (5 × 1), where ϕ1 = ϕ2 = [1z12z12-1z22z22-1]T. W1 and W2 have dimension of (5 × 1). For optimal controller the design parameters are chosen as Λ = 1.25, K = 4. For CNN based optimal adaptive controller parameter kz is chosen as 0.0001. For update of parameters in (95), Γ and k are chosen as 0.01. For approximating h(x) the input to the CNN is z1, z2, e, yd, y˙d, and y¨d. 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 W have dimension of (13 × 1). The initial conditions for [z1(0),z2(0)]T = [0.01,0.1]T.
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.
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 (v = 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.
Performance and tracking error of conventional backstepping controller.
Performance and tracking error of CNN based backstepping controller.
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.
Comparative results of driving cycle test.
Controller
Amp-hour consumption
Nonlinear optimal [11]
4.48 km/11.97 AH
Nonlinear robust [11]
4.4825 km/10.78 AH
CNN based backstepping
4.4834 km/10.75 AH
CNN based optimal adaptive
4.4795 km/10.70 AH
RMS tracking error of controllers designed.
Controller
RMS value of tracking errors
Conventional backstepping
0.0301
CNN based backstepping
0.0279
Optimal
0.0165
CNN based optimal adaptive
0.0040
Performance and tracking error of optimal controller.
Performance and tracking error of CNN based optimal adaptive controller.
Results of NEDC test of [11].
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.
Traction force.
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.
Control effort u(t).
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.
LarminieJ.LowryJ.2003West Sussex, UKJohn Wiley & SonsChengY.Van MierloJ.Van Den BosscheP.LataireP.Energy sources control and management in hybrid electric vehiclesProceedings of the 12th International Power Electronics and Motion Control Conference (EPE-PEMC '06)September 2006Portoroz, Slovenia5245302-s2.0-4254911545110.1109/EPEPEMC.2006.283214OpilaD. F.WangX.McGeeR.GillespieR. B.CookJ. A.GrizzleJ. W.An energy management controller to optimally trade off fuel economy and drivability for hybrid vehicles2012206149015052-s2.0-8155523561710.1109/TCST.2011.2168820MurpheyY. L.ParkJ.ChenZ.KuangM. L.MasrurM. A.PhillipsA. M.Intelligent hybrid vehicle power control-part I: machine learning of optimal vehicle power20126135193530MurpheyY. L.ParkJ.KiliarisL.KuangM. L.MasrurM. A.PhillipsA. M.WangQ.Intelligent hybrid vehicle power control-part II: online intelligent energy management2013626979ZhangH.TolbertL. M.OzpineciB.Impact of SiC devices on hybrid electric and plug-in hybrid electric vehicles20114729129212-s2.0-7995293660610.1109/TIA.2010.2102734DixonJ.NakashimaI.ArcosE. F.OrtúzarM.Electric vehicle using a combination of ultracapacitors and ZEBRA battery20105739439492-s2.0-7684908663010.1109/TIE.2009.2027920YeM.BaiZ.CaoB.Robust control for regenerative braking of battery electric vehicle2008212110511142-s2.0-5734918525610.1049/iet-cta:20070333PahlevaninezhadM.DasP.DrobnikJ.MoschopoulosG.JainP. K.BakhshaiA.A nonlinear optimal control approach based on the control-lyapunov function for an AC/DC converter used in electric vehicles20128596614WangR.WangJ.Passive actuator fault-tolerant control for a class of overactuated nonlinear systems and applications to electric vehicles201362972985HuangQ.HuangZ.ZhouH.Nonlinear optimal and robust speed control for a light-weighted all-electric vehicle2009344374442-s2.0-6304910606410.1049/iet-cta.2007.0367PooraniS.Udaya KumarK.RenganarayananS.Intelligent controller design for electric vehicle4Proceedings of the 57th IEEE Semiannual Vehicular Technology Conference (VTC '03)April 2003Jeju, Korea244724502-s2.0-0042469094KhatunP.BinghamC. M.SchofieldN.MellorP. H.Application of fuzzy control algorithms for electric vehicle antilock braking/traction control systems2003525135613642-s2.0-014190636710.1109/TVT.2003.815922SchwarzerV.GhorbaniR.Drive cycle generation for design optimization of electric vehicles2013628997KachroudiS.GrossardM.AbrougN.Predictive driving guidance of full electric vehicles using particle swarm optimization20126113093919NamatameA.UedaN.Pattern classification with Chebyshev neural networks199232331LeeT.-T.JengJ.-T.The Chebyshev-polynomials-based unified model neural networks for function approximation19982869259352-s2.0-003228767910.1109/3477.735405PatraJ. C.KotA. C.Nonlinear dynamic system identification using Chebyshev functional link artificial neural networks20023245055112-s2.0-003668527210.1109/TSMCB.2002.1018769PurwarS.KarI. N.JhaA. N.On-line system identification of complex systems using Chebyshev neural networks2007713643722-s2.0-3375096104210.1016/j.asoc.2005.08.001PatraJ. C.Chebyshev neural network-based model for dual-junction solar cells20112611321392-s2.0-7995194575910.1109/TEC.2010.2079935PurwarS.KarI. N.JhaA. N.Adaptive output feedback tracking control of robot manipulators using position measurements only2008344278927982-s2.0-3864909657310.1016/j.eswa.2007.05.030ZouA.-M.KumarK. D.HouZ.-G.Quaternion-based adaptive output feedback attitude control of spacecraft using chebyshev neural networks2010219145714712-s2.0-7795634476410.1109/TNN.2010.2050333ZouA.-M.KumarK. D.HouZ.-G.LiuX.Finite-time attitude tracking control for spacecraft using terminal sliding mode and chebyshev neural network20114149509632-s2.0-7996070047910.1109/TSMCB.2010.2101592WangD.HuangJ.Neural network-based adaptive dynamic surface control for a class of uncertain nonlinear systems in strict-feedback form20051611952022-s2.0-1384429925710.1109/TNN.2004.839354KhalilH. K.20023rdUpper Saddle River, NJ, USAPrentice HallKwanC.LewisF. L.Robust backstepping control of nonlinear systems using neural networks20003067537662-s2.0-003431343410.1109/3468.895898KimY. H.LewisF. L.Optimal design of CMAC neural-network controller for robot manipulators200030122312-s2.0-003374286010.1109/5326.827451NarendraK. S.AnnaswamyA. M.A new adaptive law for robust adaptation without persistent excitation19873221341452-s2.0-0023288277LewisF. L.SyrmosV. L.19952ndNew York, NY, USAJohn Wiley & Sons