Abstract
This paper addresses the issue of controlling the drivetrain of electric vehicles. Taking into account both internal and external system disturbances, including the vehicle’s mass, rotational friction of the shafts, wind speed, vehicle aerodynamics, road type, and slope constraints, the controller’s task is to ensure robustness in vehicle behavior. The significant dynamics of these disturbances and uncertainties in vehicle parameters have a substantial impact on vehicle performance. To overcome these challenges, a nonlinear model of the entire controlled system is developed. Subsequently, a robust nonlinear controller is designed using the damping function version of the backstepping design technique to compensate for all uncertain terms. Within this framework, two primary control loops are established. Firstly, a speed control loop is implemented to achieve precise tracking of the driver’s speed reference. Secondly, the machine current is optimized to generate maximum torque. A formal analysis based on Lyapunov stability is conducted to describe the control system’s performance. Despite parameter uncertainties, it is demonstrated that all control objectives are asymptotically achieved. Ultimately, all control objectives are validated through simulation results using Matlab/Simulink, showcasing the efficiency and robustness of the proposed control technique.
Similar content being viewed by others
Introduction
As the world moves toward more sustainable transportation, electric vehicles (EVs) are taking center stage in automotive innovation. This shift is driven by growing environmental concerns, stricter emission regulations, and rapid advancements in battery technology1,2. By 2040, EVs could exceed 50% of new car sales, making their role in reducing carbon emissions and improving efficiency crucial3. This demand calls for robust powertrain systems that adapt to various driving needs4. As EVs shape the future of mobility, optimizing their performance through advanced control strategies is essential.
In the ___domain of automotive engineering, cruise control systems play a pivotal role in enhancing driving comfort, fuel efficiency, and overall road safety5,6. These systems are meticulously designed to maintain a desired vehicle speed, relieving the driver from constant throttle adjustments and allowing focus on critical driving tasks7,8. The need for precise and effective control in such systems is paramount. While traditional proportional-integral-derivative (PID) controllers are widely adopted for their simplicity9,10,11, they often struggle with the nonlinear dynamics and parameter uncertainties (e.g., wind speed, vehicle mass) inherent in electric vehicle (EV) powertrains, requiring frequent retuning. Advanced model predictive control (MPC) offers predictive capabilities12,13, but its reliance on precise models and high computational demands limits real-time applicability. The literature also explores alternative structures, including fractional-order PID (FOPID)14, model predictive12,13, deep reinforcement learning15, fuzzy logic16, and real PID plus second-order derivative17 controllers.”
Despite the satisfactory performance often delivered by these controllers, their design process can be time-consuming and may not consistently yield optimal results. To address these challenges, researchers have embraced metaheuristic algorithms, showing promise in solving complex optimization problems18. In the realm of vehicle cruise control systems, the application of metaheuristic algorithms has the potential to significantly enhance the time-___domain performance and stability of the control system. Consequently, this paper focuses on developing and applying an enhanced Runge Kutta (RUN) optimizer, named IRUN, to achieve superior control performance. The proposed IRUN optimizer integrates several advanced strategies, including quadratic interpolation, Laplacian segment mutation, Levy flight, and information-sharing-based local search mechanisms, to augment its effectiveness. Through the incorporation of these strategies, the resulting IRUN algorithm showcases improved optimization capabilities, making it well-suited for fine-tuning the controller.
The study focused on comprehensive modeling of the electric vehicle’s traction system, incorporating the intricate components such as the Lithium-ion battery, the Motor-Converter association, and the mechanical transmission. This modeling approach extended to encompass the vehicle’s body dynamics and longitudinal behavior. The core innovation lay in the development and application of a robust nonlinear controller, tailored to optimize the vehicle’s performance. The methodology involved meticulous design, implementation, and validation of this controller using available empirical data. The research sought to enhance understanding and efficiency in electric vehicle powertrain systems, striving to contribute to advancements in sustainable transportation technology.
In order to comprehensively evaluate the effectiveness of the proposed method, the reference speed is provided that has to press/release either the accelerator in order to reduce the error between the actual speed and the speed from a drive cycle The driver model provides the torque demand to match the drive cycle speed profile. Thereafter, from the driver set-point, the energy required to overcome the opposing forces acting on the vehicle namely Aerodynamic Drag Force,Rolling Resistance Force, represents Gradient Resistance Force and the Inertia Resistance Force, is computed. The backward approach considers a reference speed profile, as input, to determine the forces acting at the wheels and then processes backward through the powertrain.
Summarising the aforementioned literature and considering the discussion mentioned above, the main contributions of this paper can be stated as follows:
-
Precision in Vehicle Speed Tracking: This paper addresses the critical challenge of ensuring that the vehicle speed accurately follows the reference speed set by the driver. This is particularly noteworthy given the various disturbances encountered, including wheel friction, wind velocity, vehicle aerodynamics, vehicle mass, and variations in road type or slope. The proposed approach focuses on mitigating these disturbances to enhance the accuracy of speed tracking, contributing to a safer and more comfortable driving experience.
-
Disturbance Mitigation Strategies: Innovative strategies to mitigate disturbances, contributing to a robust, adaptable vehicle speed control system.
-
Optimization of Stator Current for Maximum Torque: optimizes stator current for maximum torque, enhancing overall vehicle performance and fuel efficiency.
-
Adaptability to Variable Road Conditions: Acknowledging the influence of road type and slope on vehicle dynamics, the proposed methods aim to optimize speed control in diverse road conditions. This adaptability is a key contribution, as it enhances the overall safety and stability of the vehicle across varying terrains.
These contributions have real-world benefits for EV powertrain control. Precise speed tracking improves energy efficiency and safety, whether in stop-and-go city traffic or on steep inclines, making EVs more reliable for daily drivers. Effective disturbance mitigation ensures stable performance despite external factors like wind or changing loads, reducing strain on powertrain components and extending battery life, which is key to lowering costs. Optimizing stator current for maximum torque enhances acceleration and hill-climbing, improving the driving experience and making EVs more competitive with traditional combustion vehicles. Altogether, these advancements make EV control systems more robust and adaptable, paving the way for next-generation autonomous and high-performance electric vehicles while accelerating the shift to sustainable transportation.
The rest of the paper is organized as follows: Section “Vehicle modelling” introduces the mathematic models of the whole system. In Section “Robust nonlinear controller design”, the control strategy using nonlinear robust backstepping controller is described in details. Section “System simulation” shows the simulation setup, results and analysis. Finally, concluding remarks are discussed in Section “Robustness of the proposed control”.
Vehicle modelling
The considered conversion system consists, as shown in Fig. 1, of a association of several modules. On the one side, there is a storage battery is used to supply (reversible energy source) the power conversion chain through a DC/AC inverter. The role of this latter is to generate an adequate power supply to the electromechanical converter (synchronous machine), associated in turn with a transmission system representing the differential and the various mechanical connections between the traction machine and the vehicle’s wheels. On the other side, at the end of the conversion chain, there is the mechanical load where the rolling environment of the vehicle is considered (wind, slope, load, type of road, vehicle shape...).
Thanks to their better mass/power ratio and their ability to develop a much higher power level and present a more satisfactory efficiency, permanent magnet synchronous motors (PMSM) are more suitable for electric vehicles. A control unit allowing the control of all the elements of the traction chain is necessary in order to satisfy the driver requirements.
Model of vehicle body and longitudnal dynamics
An electric vehicle is subject to forces that the traction system must overcome in order to move the vehicle forward or backward, These forces are distributed as shown in Fig. 2.
Several more or less complex mechanical models are developed in the literature19. In this paper, the model considering the vehicle as a massive point and which undergoes a time varying forces sum has been adopted20. The following opposing forces are considered in the proposed model as a common technique in the literature as shown in Fig. 2:
-
Gradient Resistance Force: due to the road inclination with regard to the horizontal plane.
-
Rolling Resistance Force: mainly due to the friction between the tyres and the road.
-
Aerodynamic Drag Force: due to the friction between the vehicle body and the air.
-
Inertia Resistance Force: related to the forces required for the linear acceleration of the vehicle and the increase of the rotational speed of the rotating components
Where:
\(F_{aero\ drag}\) denotes the Aerodynamic Drag Force, \(F_{rolling\ resistance}\) designates the rolling resistance force, \(F_{grade}\) represents the Gradient Resistance Force, \(F_{acc}\) is the Inertia Resistance Force.
Where:
\(T_{r\omega }\) is the resistant torque on wheels,\(R_{\omega }\) designates the wheel radius,V is the vehicle speed in (m/s),\(V_{w}\) is the wind speed, which has a positive sign when it is tailwind and a negative sign when it is headwind,M is the vehicle mass in (kg),g is the gravitational constant \((m/s^2)\), \(\theta\) is the road slope angle in radians, \(\rho\) is the air density in \((kg/m^3)\),\(S_f\) is the frontal surface of the vehicle, \(C_x\) is the drag coefficient,\(C_i\) is the inertia factor, which is assumed to be around 1.15 (a typical value),\(C_{rolling\ resistance}\) is the rolling resistance coefficient.
Model of mechanical trnasmission
In order to adjust the wheels speed \((\Omega _\omega )\) and/or torque \((T_{rm})\), a gearbox linked to the drive shaft with the wheel axis. This gearbox can be modeled by a product of reduction ratio (r) and efficiency representing the friction losses \((\eta _t)\)20.
where \(P_\omega\) output mechanical power on the wheel side, \(P_m\) motor mechanical power, \(T_{rm}\) motor resistive torque; \(\Omega _\omega\) wheel angular velocity; \(\Omega\) is the angular rotor speed.
where r reduction ratio; \(\eta _t\) gearbox efficiency;
Model of Lithium-ion battery
The modeling of batteries is an important and complex issue, i.e. a mathematical description of its chemical behavior. Indeed, it is required to model its charging and discharging behaviors according to the battery type and parameters variation i.e. state of charge (SOC), current, temperature, etc. Several battery models have been developed in previous works. A compromise between complexity, accuracy and model parameterization is considered when chosing a model : i) Thevenin-based, impedance-based, and runtime-based models. ii) Model using neural network21; iii) electrochemical models22,23,24; iv) Empirical models25,26.
In the present work, the modified Shepherd model is adopted25. It is a dynamic model that can be parameterized to generate the charge and discharge curves. This model is based on equations that describe the battery open circuit voltage as a function of the state of charge and the exchanged current, as shown in Fig. 3.
The open circuit voltage (OCV) in charge and discharge modes of a Lithium-ion battery, according to the shpherd model, is described in(5) (if \(i<0\)) and (6) (if \(i>0\)) respectively (see (7)(8)).
-
Charge mode \((i<0)\)
$$\begin{aligned} OCV=&E_0-K_i^*\dfrac{Q_n}{0.1Q_n-Q_b}+KQ_b\dfrac{Q_n}{0.1Q_n-Q_b}+A\exp (-BQ_b) \end{aligned}$$(5) -
Discharge mode \((i>0)\)
$$\begin{aligned} OCV=&E_0-K_i^*\dfrac{Q_n}{Q_n-Q_b}+KQ_b\dfrac{Q_n}{0.1Q_n-Q_b}+A\exp (-BQ_b) \end{aligned}$$(6)
where \(E_0\) is the constant voltage; \(Q_n\) is the nominal capacity; \(Q_b\) is the discharged capacity; K is the polarization constant; A is the amplitude of the exponential zone; B is the time constant inverse of the exponential zone; and \(i^*\) is filtred current with low dynamic. Several methods, for estimating the battery’s SoC, exist in the literature. These methods are mainly based on the electrical quantity’s measurement (voltage and current) and battery parameters. According to Shepherd model, assume that the initial SoC is known and the battery current is available. The simplest and robust method named Coulomb counting method (CCM) can be used. The Eqs. (7) and (8) allow to compute SoC(t) based on the battery current measurement i(t).
The battery voltage \(V_b\) depends on the current direction and can be described as follows :
Generally, the battery parameters quoted in (5), (6) and (9), namely \(R_b\) which is the internal resistance of the battery K, A, B and \(E_0\) can be determined empirically by OCV\((Q_b)\) curve especially from the discharge curve or from the datasheet provided by the battery manufacturer25,27. Figure 4 shows the characteristic of voltage variation as a function of charge capacity during a discharge case28. It is a non-linear characteristic that allows the battery-parameters to be determined as follows:
where \(V_f\) is the fully charged voltage; \(V_e\) the voltage at the end of exponential zone; \(V_n\) is the rated voltage; \(Q_n\) is the maximum capacity of the battery; \(Q_e\) is the capacity at the end of exponential zone; \(i_n\) is the rated current of the battery; \(\eta\) represents the battery efficiency; \(\alpha\) and \(\beta\) constants values are usually determined to improve the fit to a battery data26. *
Model of the motor-converter association
The structure of the sub-systém considered is described by the Fig. 5. It consist of converting electric energy provided by the battery into mechanical energy, using the permanent magnet synchronous motors (PMSM)29. After the well-known Park transformation, the mathematical model of the converter-machine association in the dq frame (more suitable for developing control laws), can be described by the Eqs. (15a)–(15c)30.
where \({R_s}\) and \({L_s}\) are the stator resistor and inductance; \({i_d}\) and \({i_q}\) denotes the stator currents in dq coordinates; \({J_m}\), \({F_m}\) and p are respectively the rotor inertia, viscous coefficient and number of poles pairs; \(T_{rm}\) denote the resistant torque on the motor which include all forces applied on the vehicle ( see Eq. (4) and (2)); \({k_m=p.\phi _r}\) is the motor coefficient, with \(\phi _r\) is the rotor flux norm; \(u_d\) and \(u_q\) are the DC/AC input control in the dq coordinates; \(V_b\) is the battery voltage described in (9).
In view from the motor axis, the resistant torque \(T_{rm}\) take the following expression, obtained by replacing the Eq. (2) in (2.2).
where \(K_{D\omega }\) is the wind disturbance coeficient, \(\delta _{DM}\) is the mass disturbance function, which take variation following the state of the road and the slope percentage, \(K_i\) is the inertia coeficient.
Remarks
-
\(\delta _{DM}\) is very dependent on the mass of the vehicle that its variation is generally very limited. In this work, we consider the parameters of the function \(\delta _{DM}\) as constant and defined by their average values.
-
As the dynamics of the aerodynamic force is slowly with respect to the speed, one can approximate \(F_{aer}\) to a linear function around the rated speed. Using Taylor young’s theorem, the aerodynamic torque takes the following expression:
where \(A_i,(i=1,..4\)) are a constants coefficients depending to the vehicle parameters, as defined bellow:
Using (16) and (18) it follows from (15) that the speed dynamics undergoes the Eq. (20), defined in the follo state space representation of the whole system including the vehicle environment, mechanical transmission and motor-converter association.
where J is the whole system inertia \(J=J_m+K_i M\), \(F_v\) is the vehicle friction including the motor friction, transmission friction, and the linear terme of wind friction by the vehicle, which is expressed by \(F_v=F_m+A_3V_\omega +A_2\) and \(\psi _w\) is the non linear influence of wind velocity on vehicle taking the expression \(\psi _w=A_1V_\omega ^2+A_4\). Now, let us introduce the state variables \(x_1=\Omega\), \(x_2=i_q\) and \(x_3=i_d\) and the inputs control \(u_1=u_q\) and \(u_2=u_d\). Then the state space representation of the traction system can be presented as:
Robust nonlinear controller design
Control objectives
This section focuses on developing a nonlinear controller which is capable of achieving the following operational control objectives:
-
Speed regulation : the vehicle speed must track, as accurately as possible the reference speed specified by the driver, despite the several disturbance, in particular the wheel friction, wind velocity, vehicle aerodynamic, vehicle mass and road type or slope inconvenience.
-
Regulate the d-axis current component to zero, in order to optimize the stator current in the machine and thus develop a maximum torque31.
In the following subsections, a nonlinear state feedback controller will be performed using the Robust Backstepping technique to achieve the above control objectives. The choice of backstepping technique in this study is motivated by its systematic, recursive design process. Backstepping allows for step-by-step stabilization of subsystems, such as speed and current loops, while explicitly addressing nonlinearities and uncertainties through Lyapunov-based stability guarantees32,33.
Vehicle model with parameter uncertainties
The disturbance undergoes by the vehicule are not supposed to be known \((F_v, V_\omega )\) or not accuratly mesurable (M). It’s just supposed to be delimited within a known intervals as shown below:
with \(M_0\) the nominal mass of the vehicle, \(F_{v0}\) and \(V_{\omega 0}\) are respectively the mean values of the vehicle friction and wind velocity, concerning \(\Delta _M\), \(\Delta _F\) and \(\Delta _\omega\) denotes the possibly varying uncertainties such that.
where \(\Delta ^{MIN}\) and \(\Delta ^{MAX}\) are known bounds, using (22) in (21a), the dynamics of the speed is described as follows:
With \(\phi _1\) and \(\Delta\) are the uncertain terms defined as:
Speed regulation loop
Now, the goal is to designe the control law of the vehicule speed reference tracking, which is foist by the driver. Basing on robust backstepping technique, where the machine parameters and the vehicle disturbances must be tolerated, \(u_1\) in (21b) stand as input control of the loop, so let consider the following speed tracking error.
By substituting (24) in (26), the error dynamics can be described as follows:
To get stabilizing the subsystem (24), let consider \(\alpha =\dfrac{3}{2} \dfrac{k_m}{J}x_2\) as a virtual control input for \(z_1\) dynamics, which must track to the stabilizing function \(\alpha _{ref}\) defined as mentioned below.
Note that \(|\phi _1|\) denotes the Euclidean norm of \(\phi _1\) and \(\varsigma _1>0\) is a negative scalar design parameter. The quantity \(k_1|\phi _1|^2z_1\) is a nonlinear damping term that is introduced to dominate the unknown term \(\phi _1\Delta\).
If \(\alpha =\alpha _{ref}\), one would have:
Then, if \(V_1=0.5z_1^2\) considered as the the Lyapunov function candidate, one would get the following time-derivative along the \(z_1\) -trajectory.
Note that in case of no variation in disturbances \((\Delta =0)\), the dynamic of the Lyapunov function (30) is defined as negative function \((\dot{V}_1<0)\); consequently the \(z_1\) error vanish in a finite time. However, to ensure regulation in all cases, let introduce a new error \(z_2=\alpha -\alpha _{ref}\). The virtual control input (28) should track \(\alpha _{ref}\), in order to ensure global loop regulation and eliminate the associated error \(z_1\). Nevertheless, using the fact that \(\alpha =z_2+\alpha _{ref}\) in (27), it follows that the \(z_1\)-dynamics undergos the following equation:
The next aim is to ensure the asymptotic stability of the subsystem (21a)-(21b); i.e. the two tracking errors \((z_1,z_2)\) must converge to zero in a finite time. Thus, the time derivative of\(z_2\)-error gives
with
To recap, the errors dynamics \(\dot{z_1}\) and \(\dot{z_2}\) finally follow up the equations rewritten below:
To get the stabilizing control law for the latter errors dynamics an augmented candidate Lyapunov function is defined as:
and using (36)-(37), the time derivative of \(V_2\) yields:
Substituting the Eq. (37) in (39) turns out to be:
According to (40), to ensure the convergence of the system in a defined time, the following control law \(u_1\) is proposed.
where \(\varsigma _2>0\) is a new positive design parameter for the second loop, and \(k_2 |\phi _2|^2 z_2\) is an introduced additional nonlinear damping term to overcome the second uncertain term \(\phi _2 \Delta\). Considering (41), the time derivative of the Lyapunov function \(\dot{V_2}\) becomes:
As can be easily deduced from (42), \(\dot{V_1}\) is defined as negative function \((\dot{V_1}<0)\) although the uncertain term varies or disappears \((\Delta =0)\). The stability analysis of the overall system (21a)-(21c) will be summarized in Theorem 1. Now, one determines the remaining control input \(u_2\).
d-Axis current regulation
Now the target is the control of d-axis current \(x_3\) whose dynamics is given by the Eq. (21c). For this let’s consider the tracking error \(z_3=x_3-i_{dref}\). Then, by time derivative of \(z_3\) and using (21c), one obtains the error dynamics.
To achieve objective (ii) and make \(x_3\) converges asymptotically to \(i_{dref}=0\), let’s introduce a new quadratic Lyapunov function \(V_3=0.5z_3^2\). The asymptotic stability is provided when \(\dot{V_3}\) is negative definite, given in (44).
Where \(\varsigma _3>0\) is a new design parameter. Now, the input control of the d-axis loop can be easily deduced from (44) as follows:
Finally combining (45) and (43) gives:
To finish the design of the controller, in the following theorem, the control closed loops are analysed.
Theorem 1
Consider the traction chain system described by the nonlinear model (21a)-(21b) and the control laws (28), (41) and (45), there exist a Lyapunov function candidate \(V=V_2+V_3\) and a positive constants \(\varsigma _i\) ,\((i=1,2,3)\)and \((k_1\),\(k_2)\) such that:
-
The resulting closed-loop system undergoes, in the \((z_1, z_2, z_3)\)-coordinates, the following equation:
$$\begin{aligned} \dot{z_1}&=-(\varsigma _1+k_1|\phi _1|^2)z_1-\phi _1\Delta +z_2 \end{aligned}$$(47a)$$\begin{aligned} \dot{z_2}&=-(\varsigma _2+k_2|\phi _2|^2)z_2-\phi _2\Delta \end{aligned}$$(47b)$$\begin{aligned} \dot{z_3}&=-\varsigma _3z_3 \end{aligned}$$(47c)with \((z_1=x_1-\Omega _{ref}), (z_2=\alpha -\alpha _{ref}),\) and \((z_3= x_3-i_{dref})\)
-
The error vector \(z(t)=(z_1,z_2,z_3)\) converges exponentially to a compact neighborhood of the origin [0 0 0]. Consequently, whatever the initial conditions, the observation error z(t) can be made arbitrarily small letting \((k_1,k_2)\) be sufficiently larges.
Proof
-
(i).
Equation (47) is specifically obtained by substituting (36) into (37) and from Eqs. (27) and (46).
-
(ii).
To prove the second part (ii), let us consider the Lyapunov candidate function :
$$\begin{aligned} V=V_2+V_3=0.5(z^2_1+z^2_2+z^2_3) \end{aligned}$$(48)It easily follows that :
$$\begin{aligned} \dot{V}=&-\Big (\varsigma _1z^2_1+\varsigma _2z^2_2+\varsigma _3z^2_3+(k_1|\phi _1|^2)z^2_1+\phi _1\Delta z_1)+(k_2|\phi _2|^2)z^2_2+\phi _2\Delta z_2)-z_1z_2\Big ) \end{aligned}$$(49)\(\dot{V}\) can be bounded as follows:
$$\begin{aligned}&\dot{V}<-\left(\varsigma _1z^2_1+\varsigma_2z^2_2+\varsigma _3z^2_3-|z_1||z_2|+k_1 \left(|\phi_1||z_1|-\dfrac{||\Delta ||_\infty }{2k_1}\right)\right.\\&\qquad\left. +k_2 \left(|\phi_2||z_2|-\dfrac{||\Delta ||_\infty }{2k_2}\right)-\dfrac{||\Delta||_\infty }{4k_1}-\dfrac{||\Delta ||_\infty }{4k_2}\right)\end{aligned}$$(50)The inequality (50) can be rewritten under the following form :
$$\begin{aligned} \dot{V} \le -\varsigma V +\dfrac{||\Delta ||_\infty }{2\varsigma k} \end{aligned}$$(51)where \(||\Delta ||_\infty\) denotes the \(L_\infty\) norm and \(\varsigma = \min (\varsigma _1, \varsigma _2, \varsigma _3)\) and \(k = \min (k_1, k_2)\). Recalling that \(V=0.5(z_1^2+z_2^2+z_3^2)\), which implies that the vector \((z_1 (t)\),\(z_2 (t)\) and \(z_3 (t))\) converges exponentially to the following compact:
$$\begin{aligned} z^2_1+z^2_2+z^2_3 \le \dfrac{||\Delta ||_\infty }{\varsigma k} \end{aligned}$$(52)It is clear that the dimensions of such a compact are inversely proportional to \(\varsigma k\).
\(\square\)
System simulation
Protocol of simulation
In this section, investigates the performance of the proposed controller the module of the EV diagrams described by Fig. 6 is established in the MATLAB/Simulink environment.
To evaluate the performance of the proposed controllers, we take the wind speed profile shown in Fig. 7, and the vehicle mass shown in Fig. 8 and select a profile over a wide range of variation of the longitudinal dynamics of the vehicle, The simulation parameters are given in Table 1.
Remark 1
The control parameters of the robust backstepping controller were tuned using a trial-and-error method in MATLAB/Simulink to achieve optimal performance. This iterative process aimed to minimize key metrics such as convergence time and steady-state error. While trial-and-error proved effective for this study, alternative optimization methods could be explored for more systematic parameter tuning, such as trajectory and velocity planning using quartic bezier curves34 or QPSOMPC-based chassis coordination control35, which provide advanced optimization frameworks for vehicle control applications.
Simulation results
Use of the robust backstepping control technique, has demonstrated remarkable resilience to variations in the vehicle’s longitudinal dynamics. Figure 9a vividly illustrates the speed dynamics, demonstrating exceptional responsiveness with a very fast response time. Figures 9b and c offer zoomed views to show the response of the proposed control to rapid changes in the speed reference signal. This is a key achievement, as it signifies that the speed of the electric vehicle adeptly adheres to the reference trajectory.
In Fig. 10, we can observe the torque response, showcasing the adaptability of our control system to variations in parameters impacting longitudinal dynamics. The ability of torque to fluctuate appropriately in response to changes in vehicle mass, friction, and wind speed further exemplifies the system’s robustness. Moreover, Fig. 11 provides a snapshot of another essential control objective. Here, the \(i_{sq}\) current, which represents the machine’s torque, exhibits a commendable performance. Notably, the \(i_{sd}\) current remains consistently at zero, underscoring the effectiveness of our control strategy in both static and dynamic current responses. In sum, our control system, based on the backstepping robust control technique, excels in maintaining high-speed dynamics with a short response time, as demonstrated in Fig. 9, the speed tracking performance is indeed impressive. Figures 10 and 12 showcase the system’s adaptability, and Fig. 11 attests to the successful management of static and dynamic current responses. These results affirm the robustness and practicality of our approach, which holds promise for real-world applications.
Figure 13 demonstrates the behavior of the electrical angle \(\theta\) over time, indicating precise and stable control as highlighted at the bottom of this figure. The three-phase stator currents (\(i_{sabc}\)) are shown in Fig. 14. The zoomed-in region displays near-perfect sinusoidal waveforms for the stator currents, confirming harmonic integrity and minimal distortion during steady-state operation. The system demonstrates precise and reliable tracking of both the electrical position \(\theta\) and the stator currents (\(i_{sabc}\)) , maintaining stability and accuracy even during dynamic operating conditions.
Robustness of the proposed control
To verify the robustness of the proposed control against variations in the system characteristics, the dynamic behavior of the system was studied by introducing a \(30\%\) uncertainty on the nominal values of the machine parameters, including the stator resistance and inductance (\(R_s, L_s\)), the viscous friction coefficient (\(F_m\)), and the rotor inertia (\(J_m\)). These variations were applied to simulate realistic conditions where the machine parameters may differ from their nominal values due to wear, temperature variations, or other external factors. The performance of the proposed control was evaluated in this context, and the results are presented in Figs. 15 and 16. These figures illustrate the evolution of the machine states, including the rotor angular velocity (\(\Omega\)) and the stator currents in the three-phase reference frame (\(I_{abc}\)).
During the time interval \(\left[ 2, 2.2 \right] s\) , the actual values of the machine parameters were used in the simulated model, while the nominal values were used in the control. Despite this discrepancy, after a brief transient period, the control demonstrated the ability to maintain stable and accurate performance, ensuring that the machine states converge to their reference values.
The proposed control demonstrates stability and effectiveness even under significant variations in machine parameters, confirming its robustness against system uncertainties. The results indicate that it accurately tracks speed and current references despite disturbances caused by parameter variations, as highlighted in the zoomed view in Fig. 15. Moreover, these variations do not significantly impact overall system performance, validating the robustness of the control design.
Comparative analysis
To rigorously validate the proposed robust backstepping controller, a comparative analysis was conducted against a standard backstepping controller and a conventional PID controller, with results summarized in Table 2. All controllers were simulated in MATLAB/Simulink under identical conditions, including the reference speed profile (Figure 7) and vehicle mass variations (Figure 8). The proposed robust backstepping controller achieves a root mean square error (RMSE) of 0.12rad/s, significantly outperforming the standard backstepping controller (0.25rad/s) and PID controller (0.35rad/s), reflecting superior speed tracking precision. Its convergence time to within \(2\%\) of the reference speed is 0.15s, compared to 0.22s for standard backstepping and 0.28 s for PID, indicating faster response dynamics. The peak overshoot is reduced to \(1.5\%\) versus \(2.8\%\) and \(4.2\%\), respectively, while the steady-state error drops to 0.03rad/s from 0.07rad/s and 0.10rad/s. Moreover, under varying wind speed and mass conditions, the proposed controller maintains robust performance with RMSE values of 0.14rad/s and 0.13rad/s, compared to 0.28rad/s and 0.26rad/s for standard backstepping, and 0.40rad/s and 0.38rad/s for PID. These results, highlight the proposed approach’s enhanced accuracy, responsiveness, and resilience to parameter uncertainties.
Conclusion
In conclusion, the primary objective of this study was to develop a precise computer-based model for Electric Vehicle (EV) energy consumption during various driving cycles. A comprehensive forward vehicle simulation model was evaluated through a simulation using MATLAB/Simulink, encompassing the intricate powertrain system and longitudinal vehicle dynamics. This model integrated the Thevenin equivalent circuit battery model, inverter, and permanent magnet synchronous motor (PMSM), contributing to a holistic roepresentation of the EV’s energy dynamics. Additionally, the resistance forces opposing the vehicle’s motion were accurately characterized within the longitudinal vehicle dynamics, further enhancing the model’s accuracy. To regulate the vehicle’s speed, a robust nonlinear controller, specifically the backstepping robust technique, was adeptly designed and incorporated into an elaborated driver model. A formal stability analysis employing Lyapunov tools provided assurance that the system operates in a stable manner, crucial for real-world applications. The results obtained from this comprehensive model undeniably demonstrate its effectiveness in achieving the defined control objectives with remarkable performance and robustness. Looking ahead, future work could focus on integrating real-time adaptive algorithms to enhance the controller’s responsiveness to unpredictable environmental conditions, such as sudden weather changes or road surface variations. Additionally, experimental validation on a physical EV prototype would help bridge the gap between simulation and practical implementation. Investigating the scalability of the proposed approach across different EV architectures (e.g., hybrid or multi-motor systems) could further broaden its applicability.
Data availibility
The data that supports the findings of this study are available from the corresponding author on request.
Abbreviations
- CCM:
-
Coulomb counting method
- DC/AC:
-
Direct current/Alternating current
- EV:
-
Electric vehicle
- OCV:
-
Open circuit voltage
- PID:
-
Proportional-integral-derivative
- PMSM:
-
Permanent magnet synchronous motors
- RUN:
-
Runge kutta
- SOC:
-
State of charge
- RMSE:
-
Root mean square error
- MPC:
-
Model predictive control
References
Gnanavendan, S. et al. Challenges, solutions and future trends in ev-technology: A review. IEEE Access 12, 17242–17260. https://doi.org/10.1109/ACCESS.2024.3353378 (2024).
Ashfaq, S. et al. Comparing the role of long duration energy storage technologies for zero-carbon electricity systems. IEEE Access 12, 73169–73186. https://doi.org/10.1109/ACCESS.2024.3397918 (2024).
Dalkic-Melek, G., Saltik, E. C. & Tuydes-Yaman, H. Electric vehicle (ev) market penetration in countries with rising motorization rates. Int. J. Civ. Eng. https://doi.org/10.1007/s40999-024-01039-z (2024).
Samiullah, M., Amrr, S. M., Al-Awami, A. T., AbouOmar, M. S. & El Ferik, S. Robust high gain control for uncertain hybrid powertrain systems in fuel cell electric vehicles. Transp. Res. Procedia 84, 400–407. https://doi.org/10.1016/j.trpro.2025.03.089 (2025).
Vasebi, S., Hayeri, Y. M. & Saghiri, A. M. A literature review of energy optimal adaptive cruise control algorithms. IEEE Access https://doi.org/10.1109/ACCESS.2023.3241140 (2023).
Liu, M., Chen, W., Huang, J. & Ning, Y. Kalman filter-based model predictive control for an adaptive cruise control system considering measurement noise. SAE Int. J. Connect. Autom. Veh. 3, 53–66. https://doi.org/10.4271/12-03-01-0005 (2020).
Theodosis, D., Karafyllis, I. & Papageorgiou, M. Cruise controllers for lane-free ring-roads based on control lyapunov functions. J. Franklin Inst. 360, 6131–6161. https://doi.org/10.1016/j.jfranklin.2023.02.003 (2023).
Ekinci, S. et al. Revolutionizing vehicle cruise control: An elite opposition-based pattern search mechanism augmented info algorithm for enhanced controller design. Int. J. Comput. Intell. Syst. 16, 129. https://doi.org/10.1007/s44196-023-00304-8 (2023).
Chaturvedi, S. & Kumar, N. Design and implementation of an optimized pid controller for the adaptive cruise control system. IETE J. Res. https://doi.org/10.1080/03772063.2021.2012282 (2021).
Izci, D. & Ekinci, S. A novel hybrid aso-nm algorithm and its application to automobile cruise control system. In Proceedings of 2nd International Conference on Artificial Intelligence: Advances and Applications: ICAIAA 2021, 333–343, https://doi.org/10.1007/978-981-16-6332-1_29 (Springer, 2022).
Izci, D., Ekinci, S., Kayri, M. & Eker, E. A novel improved arithmetic optimization algorithm for optimal design of pid controlled and bode’s ideal transfer function based automobile cruise control system. Evol. Syst. 13, 453–468 (2022).
Zhang, R. et al. Constrained hybrid optimal model predictive control for intelligent electric vehicle adaptive cruise using energy storage management strategy. J. Energy Storage 65, 107383. https://doi.org/10.1016/j.est.2023.107383 (2023).
Landolfi, E., Salvi, A., Troiano, A. & Natale, C. Hardware-in-the-loop validation of an adaptive model predictive control on a connected and automated vehicle. Int. J. Adapt. Control Signal Process. https://doi.org/10.1002/acs.3583 (2023).
Pradhan, R. & Pati, B. B. Optimal fopid controller for an automobile cruise control system. In 2018 International Conference on Recent Innovations in Electrical, Electronics & Communication Engineering (ICRIEECE), 1436–1440 (IEEE, 2018).
Selvaraj, D. C., Hegde, S., Amati, N., Deflorio, F. & Chiasserini, C. F. A deep reinforcement learning approach for efficient, safe and comfortable driving. Appl. Sci. 13, 5272 (2023).
Zhu, M. & Tan, G. Research on cooperative adaptive cruise control (cacc) based on fuzzy pid algorithm (Tech. Rep, SAE Technical Paper, 2023).
Ekinci, S., Izci, D., Abualigah, L. & Zitar, R. A. A modified oppositional chaotic local search strategy based aquila optimizer to design an effective controller for vehicle cruise control system. J. Bionic Eng. 20(4), 1828–1851 (2023).
Izci, D. & Ekinci, S. An improved run optimizer based real pid plus second-order derivative controller design as a novel method to enhance transient response and robustness of an automatic voltage regulator. E-Prime-Adv. Electr. Eng., Electron. Energy 2, 100071 (2022).
Sanguesa, J. A., Torres-Sanz, V., Garrido, P., Martinez, F. J. & Marquez-Barja, J. M. A review on electric vehicles: Technologies and challenges. Smart Cities 4, 372–404 (2021).
Vafamand, N., Arefi, M. M., Khooban, M. H., Dragičević, T. & Blaabjerg, F. Nonlinear model predictive speed control of electric vehicles represented by linear parameter varying models with bias terms. IEEE J. Emerg. Sel. Top. Power Electron. 7, 2081–2089 (2018).
Dang, X. et al. Open-circuit voltage-based state of charge estimation of lithium-ion battery using dual neural network fusion battery model. Electrochim. Acta 188, 356–366 (2016).
Drummond, R. & Duncan, S. R. Structural identifiability of a pseudo-2d li-ion battery electrochemical model. IFAC-PapersOnLine 53, 12452–12458 (2020).
Li, J., Adewuyi, K., Lotfi, N., Landers, R. G. & Park, J. A single particle model with chemical/mechanical degradation physics for lithium ion battery state of health (soh) estimation. Appl. Energy 212, 1178–1190 (2018).
Zhang, C. et al. A generalized soc-ocv model for lithium-ion batteries and the soc estimation for lnmco battery. Energies 9, 900 (2016).
Raszmann, E., Baker, K., Shi, Y. & Christensen, D. Modeling stationary lithium-ion batteries for optimization and predictive control. In 2017 IEEE Power and Energy Conference at Illinois (PECI), 1–7 (IEEE, 2017).
Hemi, H., M’Sirdi, N. K. & Naamane, A. A new proposed shepherd model of a li-ion open circuit battery based on data fitting. In IMAACA 2019 (2019).
Li, S. & Ke, B. Study of battery modeling using mathematical and circuit oriented approaches. In 2011 IEEE Power and Energy Society General Meeting, 1–8 (IEEE, 2011).
El Mezdi, K., El Magri, A., Watil, A., El Myasse, I. & Bahatti, L. Integrated control and energy flow management for hybrid grid-connected photovoltaic/wind systems with battery storage using fuzzy logic controllers. IFAC-PapersOnLine 58, 442–447. https://doi.org/10.1016/j.ifacol.2024.07.522 (2024).
Mansouri, A. et al. Designing control systems and multimode power management of wecs including energy storage. IFAC-PapersOnLine 58, 7–12. https://doi.org/10.1016/j.ifacol.2024.07.451 (2024).
Mansouri, A. et al. Theoretical development and experimental evaluation of a nonlinear observer for sensorless wecs. IFAC-PapersOnLine 58, 1–6. https://doi.org/10.1016/j.ifacol.2024.07.450 (2024).
Han, Z., Xu, N., Chen, H., Huang, Y. & Zhao, B. Energy-efficient control of electric vehicles based on linear quadratic regulator and phase plane analysis. Appl. Energy 213, 639–657 (2018).
El Myasse, I. et al. Improvement of real-time state estimation performance in hvdc systems using an adaptive nonlinear observer. IFAC J. Syst. Control. 27, 100244. https://doi.org/10.1016/j.ifacsc.2024.100244 (2024).
El Myasse, I. et al. Design and implementation of advanced output feedback control for vsc-hvdc transmission system. IFAC-PapersOnLine 58, 521–526. https://doi.org/10.1016/j.ifacol.2024.07.535 (2024).
Chen, T., Cai, Y., Chen, L. & Xu, X. Trajectory and velocity planning method of emergency rescue vehicle based on segmented three-dimensional quartic bezier curve. IEEE Trans. Intell. Transp. Syst. 24, 3461–3475. https://doi.org/10.1109/TITS.2022.3224785 (2023).
Chen, T., Xu, X., Cai, Y., Chen, L. & Li, K. Qpsompc-based chassis coordination control of 6widagv for vehicle stability and trajectory tracking. J. Franklin Inst. 362, 107458. https://doi.org/10.1016/j.jfranklin.2024.107458 (2025).
Miri, I., Fotouhi, A. & Ewin, N. Electric vehicle energy consumption modelling and estimation-a case study. Int. J. Energy Res. 45, 501–520. https://doi.org/10.1002/er.5700 (2021).
Funding
Open access funding provided by Manipal Academy of Higher Education, Manipal
Author information
Authors and Affiliations
Contributions
All authors contributed equally to the main contributor to this paper. All authors read and approved the final paper.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
El Myasse, I., Lmouradi, M., El Magri, A. et al. Enhancing electric vehicle powertrain energy efficiency using robust nonlinear control approaches. Sci Rep 15, 19973 (2025). https://doi.org/10.1038/s41598-025-04950-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-025-04950-0