HUNGARIAN JOURNAL OF INDUSTRY AND CHEMISTRY Vol. 49(1) pp. 37–46 (2021) hjic.mk.uni-pannon.hu DOI: 10.33927/hjic-2021-06 ELECTRIC VEHICLE MODELLING AND SIMULATION OF A LIGHT COMMERCIAL VEHICLE USING PMSM PROPULSION AMINU BABANGIDA*1 AND PÉTER TAMÁS SZEMES1 1Department of Mechatronics Engineering, University of Debrecen, Debrecen, HUNGARY Even though the Internal Combustion Engine (ICE) used in conventional vehicles is one of the major causes of global warming and air pollution, the emission of toxic gases is also harmful to living organisms. Electric propulsion has been developed in modern electric vehicles to replace the ICE. The aim of this research is to use both the Simulink and Simscape toolboxes in MATLAB to model the dynamics of a light commercial vehicle powered by electric propulsion. This research focuses on a Volkswagen Crafter with a diesel propulsion engine manufactured in 2020. A rear-wheel driven electric powertrain based on a Permanent Magnet Synchronous Motor was designed to replace its front-wheel driven diesel engine in an urban environment at low average speeds. In this research, a Nissan Leaf battery with a nominal voltage of 360 V and a capacity of 24 kWh was modelled to serve as the energy source of the electric drivetrain. The New European Driving Cycle was used in this research to evaluate the electric propulsion. Another test input such as a speed ramp was also used to test the vehicle under different road conditions. A Proportional Integral controller was applied to control the speed of both the vehicle and synchronous motor. Different driving cycles were used to test the vehicle. The vehicle demonstrated a good tracking capability in each type of test. In addition, this research determined that the fuel economy of electric vehicles is approximately 19% better than that of conventional vehicles. Keywords: MATLAB, New European Driving Cycle, Permanent Magnet Synchronous Motor, Propor- tional Integral, Volkswagen, Internal Combustion Engine 1. Introduction Environmental effects such as air pollution and global warming which are harmful to our health are primarily caused by internal combustion engines (ICE) in conven- tional vehicles. “In recent decades, the research and de- velopment activities related to transportation have em- phasized the development of high-efficiency, clean and safe transportation.” [1] Nowadays, electric vehicles are being developed to reduce these toxic effects and achieve safer transportation networks. In Ref. [2], Un-Noor et al. stated that the development of electric vehicles is im- mensely beneficial to our environment since this will lead to a reduction in the adverse effects of greenhouse gas emissions. After reviewing the literature on Electric Ve- hicles (EV), they drew up different design and develop- ment processes in terms of vehicle modelling, EV con- figurations, battery management and electrical machine drives. In EVs, ICE propulsion has been replaced by elec- tric propulsion, consisting of electric motor drives, en- ergy sources and other auxiliaries. Therefore, a reason- able amount of effort has been made in the field of in- dustrial automation to make the transition from vehicles powered by traditional ICEs to those driven by EVs [3]. However, Shariff et al. [4] stated that ‘Greenhouse gas *Correspondence: aminubabangida24@gmail.com emission and the increased cost of petroleum products are the major factors that need a shift from internal combus- tion engines to Electric Vehicles.” Electric vehicles are solutions to this problem [5]. In this paper, a Volkswagen Crafter with a 2.0 diesel TDI CR engine manufactured in 2020 is examined by focusing on the replacement of its ICE with rear-wheel-driven electric motor propulsion. According to our literature review, numerous papers have focused on electric vehicles, e.g., Wahono et al. [6] who compared three forms of range extender engines for electric cars based on simulations to overcome some dis- advantages (such as the weight) of EVs over conventional vehicles driven by ICEs. A further study by Marmaras et al. [7] simulated the driver behaviour of EVs in road transport networks and electrical grids. The EV studied was modelled to investigate its integration in both the electrical grid and road transport networks. Although a multi-agent platform was used to model driver behaviour, a fleet of 1000 EVs were used as a case study where known and unknown profiles were chosen to explore the results. The cars were considered to be typical smart mechatronics systems. New trends in smart systems are summarised in Ref. [8]. Since it is difficult to formulate the exact dynamic equations of a car, Aracil et al. [9] pro- posed the Hardware-In-the-Loop simulation as an opti- https://doi.org/10.33927/hjic-2021-06 mailto:aminubabangida24@gmail.com 38 BABANGIDA AND SZEMES Figure 1: Battery electric vehicle mal solution. They discovered that electric vehicles have direct and indirect implications on road transport networks and elec- trical grids. Another study [10] presented a “ride com- fort performance evaluation on EV conversion via simu- lations.” This study aimed to investigate the ride comfort of a vehicle before being converted into an electric ve- hicle. The study considered a full car model with 7 de- grees of freedom (DOF). The two results were validated by evaluating the performance of the vehicle before and after replacing its conventional ICE with an electric mo- tor. 2. Electric powertrain The electric powertrain of a Battery Electric Vehicle (BEV) consists of two components, namely the electrical and mechanical parts. A schematic diagram of the gen- eral layout of the electric powertrain of a BEV is depicted in Fig. 1 where M denotes the electric motor, while the thin and thick lines represent the electrical and mechani- cal parts, respectively. 2.1 Electrical parts The electrical parts consist of the battery, DC-DC con- verter, inverter and the controller of the electrical ma- chine, which all depend on the electrical machine applied as the EV. 2.2 Mechanical parts The mechanical parts of the electric powertrain con- sist of the transmission system, axles, wheels and chas- sis. In this paper, a three-phase Permanent Magnet Syn- chronous Motor (PMSM) was applied in a simplified energy-equivalent PMSM model that makes use of the losses resulting from the detailed PMSM model. 3. Simulation of the dynamic system of the electric powertrain A simplified closed-loop representation of the electric powertrain applied in our electric Crafter (e-Crafter) is presented in Fig. 2, which consists of a Nissan Leaf bat- tery as well as the PMSM drives, transmission and chas- sis subsystems adapted from Refs. [11] and [12], where Vref denotes the speed reference and Vehspd stands for the vehicle speed. Figure 2: Simplified EV powertrain Figure 3: Transmission system [1] 3.1 Vehicle transmission system A single-speed transmission system consists of various elements such as the gearbox, torque converter and the final drive. As is described in Ref. [1], the torque con- verter couples the gearbox to the vehicle and the gearbox contains the appropriate gear ratios, where U denotes the velocity required and V represents the actual velocity as a result of the final drive. However, Fig. 3 shows a simplified single speed trans- mission used in modeling our electric vehicle. The pres- ence of the torque converter in the transmission system of Fig. 3 clearly indicates that it is automatic transmission. 3.2 Vehicle modelling Vehicle dynamics is the study of the motion of a ve- hicle and is comprised of three categories, that is, lon- gitudinal, lateral and vertical dynamics. In this paper, the longitudinal dynamics of the car are modelled in a MATLAB/Simscape/Simulink environment. “In practical terms, a vehicle not only travels on a level road but also up and down the slope of a roadway as well as around corners.” [13] A simplified model and more detailed de- scription of dynamical behaviour can be found in Ref. [8]. To model the vehicle dynamics, it is necessary to de- scribe the forces acting on the vehicle using a Free Body Diagram as shown in Fig. 4. The tractive force acting on the chassis can be de- scribed by [3, 14] Ft = Fad + Frr + Fhc + FA, (1) where Fad denotes the wind resistance, which depends on the density of air ρ, surface area of the front of the vehicle Af , drag coefficient Cd and its speed V , calculated from Fad = 0.5ρCdAfV 2. (2) The rolling resistance Frr depends on the weight of the vehicle w (w = mg), rolling resistance coefficient Crr and the angle of inclination α: Frr = wCrr cos α. (3) Hungarian Journal of Industry and Chemistry ELECTRIC VEHICLE MODELLING AND SIMULATION 39 Figure 4: Vehicle dynamics [3] Table 1: Vehicle specifications Parameters Specifications Vehicle Mass 3500 kg Centre of Gravity 0.254 m Front Axle 1 m Rear Axle 1.346 m Rolling Resistance 0.013 Drag Coefficient 0.3 Air Density 1.225 kg/m3 Gravity 9.81 m/s2 The forces resulting from the grade resistance and resis- tance to acceleration are given as Fhc = w sin α, (4) FA = 1.04 ma, (5) where m denotes the mass of the vehicle in kg, 1.04 is its inertia and a stands for its acceleration. The tractive power and energy needed to propel the vehicle are given by the following equations: P = FV, (6) E = Pt, (7) where F denotes the tractive force in newtons and t repre- sents the time in seconds. However, the vehicle dynamics system was simulated using the parameters as specified in Table 1. 3.3 Tyre Dynamics using the Magic Formula “The Tire-Road Interaction (Magic Formula) block mod- els the longitudinal forces at the tire-road contact patch using the Magic Formula of Pacejka.” [25] In this pa- per, both the tyres attached to the front and rear axles of the vehicle were modelled using the Magic Formula. The tyre coefficients used were B, C, D and E. The values of these coefficients, adapted from Ref. [25], are shown in Table 2. Table 2: Tyre specifications [25] Surfaces Constant Coefficients B C D E Dry Tarmac 10 1.9 1 0.97 Wet Tarmac 12 2.3 0.82 1 Snow 5 2 0.3 1 Ice 4 2 0.1 1 Figure 5: Control loop by applying a PI Controller 3.4 Speed controller A PI controller was developed to control the speed of both the motor and vehicle. This PI controller was im- plemented in the energy-equivalent model of our PMSM. PID controllers are used in many industrial applications because of their simple structure and robustness [15]. Since noise is a measured parameter, the derivative part is not usually used [15]. The general representation of the PI controller is presented in Fig. 5. 3.5 Integral performance criteria “Criteria based on disturbance rejection and system ro- bustness are proposed to assess the performance of PID controllers.” [16] “A two-block structured singular value measures the robustness, and the disturbance rejection is measured by the minimum singular value of the in- tegral gain matrix.” [16] In this paper, five criteria used in a closed-loop control system are employed to assess the performance of our PI controller. They are customar- ily calculated for different control setpoints such as step input and ramp input. In this research, the performance of our controller was assessed using various test inputs. “It is well-known that a well-designed control system should meet the disturbance attenuation, setpoint track- ing, robust stability, and robust performance.” [16] “The first two requirements are traditionally referred to as ’per- formance’ and the third, ’robustness’ of a control sys- tem.” [16] The following are the criteria stated in Ref. [16]: IAE = ∫ ∞ 0 |e(t)|dt, (8) ITAE = ∫ ∞ 0 t |e(t)|dt, (9) ISE = ∫ ∞ 0 e2(t) dt, (10) 49(1) pp. 37–46 (2021) 40 BABANGIDA AND SZEMES Figure 6: Structure of the PMSM [17] ITSE = ∫ ∞ 0 te2(t) dt, (11) where IAE stands for the Integral Absolute Error, ITAE the Integral Time Absolute Error, ISE the Integral Square Error, and ITSE the Integral time Squared Absolute Er- ror [16]. These performance indexes were used to tune our PID controller. 4. Permanent magnet synchronous motor A PMSM, which is widely used to overcome the disad- vantages of a Brushless DC Motor (BLDC) [17], is pro- posed in this research. Virani et al. employed the Field- Oriented Control (FOC) approach to control the speed and torque of the PMSM of an electric car [17]. More- over, Espina et al. in their review of Speed Anti-Windup PI strategies for Field-Oriented Control of Permanent Magnet Synchronous Motors emphasized that PMSMs are gaining popularity when compared to other AC Mo- tors due to their higher efficiency, lower inertia as well as reduction in weight and volume [18]. This study suggests PMSM has advantages in EV ap- plications over other types of electric motors. The gen- eral structure of the PMSM motor is presented as shown in Fig. 6. The rotor having a permanent magnet mounted on it creates a rotating magnetic field, which in turn pro- duces a sinusoidal electromagnetic field. 4.1 PMSM Mathematical Model The modelling of a PMSM was carried out based on the following assumptions [17, 26]: 1. There is distribution of the sinusoidal Magnetomo- tive Force (MMF) in the air gap. 2. Restriction in the saliency according to the rotor po- sition. 3. Ignoring the hysteresis and saturation. 4. Assuming a balanced 3-phase supply voltage. 5. Assuming that the back EMF (electromotive force) is sinusoidal. The 3-phase supply voltage is given by [17] VA = PΨA + IARS, (12) VB = PΨB + IBRS, (13) VC = PΨC + ICRS, (14) where IA, IB and IC denote the phase currents, VA, VB, and VC represent the phase voltages, ΨA, ΨB and ΨC stand for the flux linkages, and RS and P are the phase resistance and "derivative operator," respectively [17]. However, using the reference frame dq, the model of the PMSM can be represented in the rotating reference frame dq as Vq = RSIq + ωrλd + Pλp, (15) Vd = RSId −ωrλp + Pλq, (16) λq = LqIq and λd = LdId + λr, (17) Vq = RSIq + ωr (LdId + λr) + PLqIq, (18) Vd = RSId −ωrLqIq + P(LdId + λr), (19) Therefore, the torque developed by the PMSM is given by Te = 3 2 P 2 (ΨrIq + (Ld −Lq)IdIq) , (20) where P denotes the number of poles of the ma- chine [17]: the electric torque derived in Eq. 20 is di- vided into two components, namely the “mutual reac- tance torque” [17] and “reluctance torque”, the latter re- sults from the difference in reluctance between the q- and d-axes [17]. However, for the PMSM, when Lq = Ld = LS, the torque generated by the PMSM is [17] Te = 3 2 P 2 (ΨrIq) . (21) The three-phase voltages of the detailed PMSM model are presented in Fig. 7. Since it can be seen that the simu- lation runs with a stop time of 0.2 s, which is very slow, it runs more slowly than in real time. To resolve this issue, a methodology was adapted using this detailed model of the PMSM to obtain the electrical losses of an energy- equivalent model. The corresponding 3-phase currents obtained by sim- ulating the detailed three-phase PMSM are shown in Fig. 8. In this model, the torque induced by the structure of the chassis is not included in the modelling process, namely the cogging torque [19]. Flux harmonics are also present since the magnet in the PMSM is composed of “neodymium, iron and boron”. [19] Therefore, its mag- netic flux density is usually affected by variations in tem- perature. [19] Similarly, there parametric uncertainties are present due to mechanical and electrical parameters. In terms Hungarian Journal of Industry and Chemistry ELECTRIC VEHICLE MODELLING AND SIMULATION 41 Figure 7: Three-phase voltages Figure 8: Three-phase currents of mechanical ones, the inertia of the PMSM is uncer- tain due to its changing behaviour under different operat- ing conditions. Regarding electrical ones, the stator resis- tance, which is a function of the temperature, influences the control of the current loop performance [19]. The electrical losses resulting from the simulation of the PMSM are presented in Table 3. These losses were obtained by using the corresponding torque and speed vectors, moreover, this method was adopted from Ref. [11]. 4.2 Simplified PMSM model In this paper, a three-phase model of the PMSM designed was converted into an equivalent energy model based Table 3: PMSM electrical losses Speed Torque experimental electrical losses (kW) [rpm] [Nm] 99 9 0.0811 0.321 0.943 2.146 4.379 454 45 0.0621 0.299 0.925 2.191 4.567 803 80 0.0454 0.281 0.9 2.217 4.796 1149 114 0.0298 0.262 0.858 2.217 4.950 1499 149 0.016 0.251 0.873 2.23 4.998 Figure 9: PMSM equivalent energy model Table 4: Motor specifications Parameters Specifications Maximum Power 80 kW Maximum Torque 280 Nm Time Constant 0.02 s Series Resistance 0 Rotor Inertia 3.9×10−4 kg m2 Rotor Damping 10−5 Nm/(rad/s) on the electrical losses obtained from the detailed three- phase model of the PMSM. The motor parameters used during the simulation are presented in Table 4. From Fig. 9, the transfer function from TL to ω is given by G(s) = 1 sJ + B , (22) H(s) = ω ωr = F(s)G(s) 1 + F(s)G(s) = α/s 1 + α/s , (23) where TL denotes the load torque, ω represents the achieved speed, ωr refers to the reference speed, J stands for the total of the moments of inertia, B is the coefficient of viscosity and α denotes the bandwidth of the speed control [20] F(s) = α s ( 1 sJ + B ) . (24) Transforming Eq. 3 into the PI form [14] yields F(s) = αJ + αB s = Kp + Ki s . (25) Therefore, in the case of the simplified equivalent PMSM Model, only the outer loop was analysed, which uses a PI controller to control the speed of the motor. Generally, a three-phase PMSM, using the FOC (Field-Oriented Con- trol) strategy, consists of two control loops with two PIs in the inner loop to control the current vectors. The PI controller can be mathematically represented as u(t) = kpe(t) + ki ∫ e(t) dt, (26) where kp and ki denote the proportional and integral co- efficients, respectively and e(t) represents the error be- tween the reference and the feedback signal [21]. 49(1) pp. 37–46 (2021) 42 BABANGIDA AND SZEMES Table 5: Battery parameters [13] Chemistry Symbol Cell Voltage Specific Energy Cycle Life Specific Power Self-Discharge (V) (Wh/kg) (W/kg) (per month) Lead-Acid PbA 2 35 ≈ 500 250 − 500 5 Nickel-metal hydride NiMH 1.2 30 − 100 > 1000 200 − 600 > 10 Lithium-ion Li-ion 3.8 80 − 160 > 1000 250 − 600 < 2 Lithium-titanate LTO 2.5 50 − 100 > 20, 000 N/A N/A Alkaline ZnMnO2 1.5 110 N/A N/A < 0.3 Table 6: Battery specifications Parameters Specifications Battery Nominal Voltage 360 V Battery Capacity 24 kWh Battery Charge 66.2 Ah Energy Density 140 Wh/kg Power Density 2.5 kW/kg Battery Power 90 kW 5. Nissan Leaf Battery In this research, a Nissan Leaf Battery manufactured in 2011 was used to design the electric vehicle. The bat- tery has a nominal voltage of 360 V and a capacity of 24 kWh [13]. It consists of 48 modules and 4 cells (2 in par- allel, 2 in series), amounting to 192 cells [14]. The total voltage of the battery pack is approximately 403.2 V. The possible arrangement of the cells in a battery pack has been studied in general by Emadi [22]. Table 5 shows the batteries available along with their chemistries. In this paper, a lithium-ion battery was used. The specific energy defines the energy stored in the bat- tery per unit of weight. The cell voltage is 4.2 V when fully charged and 2.5 V when discharged. While the battery specifications for a Nissan Leaf manufactured in 2011 are detailed in Table 6, this pa- per used a model of a built-in battery simulated using Simulink based on these specifications to simulate our traction battery. 5.1 Battery modelling Generally, our battery was analysed using an equiva- lent circuit model from the literature to study its be- haviour mathematically. This equivalent circuit was mod- elled based on a Nissan Leaf to carry out the analysis be- fore being compared with the built-in battery simulated using MATLAB. The cell voltage of the battery as stud- ied in Ref. [23] is given by Vb = V 0 r − RT nF ln QR, (27) where R denotes the ideal gas constant, T represents the temperature and QR stands for the reaction quotient, which is a function of the concentrations of the reactants. Figure 10: Battery equivalent circuit model However, since each battery is associated with an ohmic drop [23], the equation can be modified as Vb = V 0 r − RT nF ln QR −RbIb. (28) The above equation can then be modified to include the battery capacity as Vb (Ib,y) = V 0 r −A ln (By)−Ky−Fe G(y−y3)−RbIb, (29) where y denotes a variable that can be related to the ca- pacity, DoD (depth of discharge), SoC (state of charge) or the cell energy, y3 represents the value at which the ex- ponential decay begins, and A, B, K, F and G stand for values determined by curve fitting [23]. The above equa- tion can be expressed in terms of the DoD as Vb(Ib, DoD) = V 0 r −A ln(B DoD)− −K DoD −FeG(DoD−DoD3) −RbIb (30) In terms of the no-load, the ohmic drop can be expressed as Vb(nl)(DoD) = V 0 r −A ln(B DoD)−Fe G (DoD−DoD3). (31) The aforementioned equations can be modelled to repre- sent the equivalent circuit shown in Fig. 10. Therefore, the above mathematical model of the battery can be represented by the simplified battery model presented in Fig. 10 based on the already built-in Simulink model. 5.2 Battery parameters The design considerations of a battery ensure it functions safely and reliably. Therefore, the battery management Hungarian Journal of Industry and Chemistry ELECTRIC VEHICLE MODELLING AND SIMULATION 43 Figure 11: Simulink model of the EV system takes into account three parameters, namely the total voltage of the battery pack, the total temperature of an individual cell in the pack and the total current, before calculating the State of Charge (SoC), State of Health (SoH), Safe Operating Envelope (SOE) and faults. The SoC, which is expressed as a percentage, determines the amount of voltage also as a percentage. In this paper, the SoC was found to be 99%, that is, almost fully charged. A friction brake was used to stop our vehicle. The bat- tery was naturally recharged due to regenerative braking. The SoH represents the battery’s capacity relative to its capacity when initially installed. Finally, the SOE shows the amount of current that can be charged or discharged at any given time. 6. Matlab model The complete model of our electric vehicle simulated in a Simulink environment is depicted in Fig. 11. The model consists of the Nissan Leaf battery, the PMSM drives, the transmission systems and the vehicle subsystems. 7. Results The "New European Driving Cycle (NEDC)" used in this research to test our vehicle and other driving cycles to comply with energy consumption and emissions reduc- tion targets is presented in Fig. 12. In this paper, only three cycles of 200 seconds in duration were used in the "NEDC." Other driving cycles such as the "Urban Dy- namometer Driving Schedule (UDDS)," were used to test our electric vehicle. The NEDC, adopted from the literature, is mainly used to determine the consumption of electric vehicles, gas emissions, etc. The reference speed and the speed achieved by the vehicle over 200 seconds are both rep- resented in Fig. 13. The power consumption of the bat- tery in kilowatts is shown in Fig. 14. The Nissan Leaf can provide over 90 kW of power, 50 kW of which was con- sumed by our electric vehicle. The energy consumption of the battery in kWh is presented in Fig. 15. Given that the SoC achieved by our battery was 99%, i.e., it was almost fully charged, the operating tempera- ture used by our Battery Management System was ne- Figure 12: NEDC driving cycle Figure 13: Speed achieved by the vehicle Figure 14: Power consumption of the battery Figure 15: Energy consumption of the battery 49(1) pp. 37–46 (2021) 44 BABANGIDA AND SZEMES Figure 16: Battery current Figure 17: Battery SoC Figure 18: Achieved motor speed Figure 19: Motor torque Figure 20: Power consumption of the motor Figure 21: Energy consumption of the motor and battery glected. In a future study, a complete Battery Manage- ment System will be implemented to determine and en- sure safe operating conditions for our battery system. The battery current is represented in Fig. 16. The bat- tery draws a maximum current of 75 A when the ve- hicle is accelerating and a minimum current of −55 A when decelerating. The SoC of the battery as a percent- age, which is approximately maintained at 99% due to natural regenerative braking, is presented in Fig. 17: The reference speed of the motor was 4250 rpm, which accurately tracks the actual rate, and is shown in Fig. 18. The maximum torque of 280 Nm, as seen in Table 4 of the motor specifications, is shown in Fig. 19, while the power consumption of the motor, which is approxi- mately 49 kW, is presented in Fig. 20. The losses can be analysed in terms of the energy transmitted from the stor- age system, namely from the battery to the wheels of the vehicle. The energy consumption of both the battery and motor are plotted in Fig. 21. The difference between the plots in Fig. 21 is the en- ergy lost when transmitting from the energy supply sys- tem, that is, from the battery to our PMSM. It can be seen that: if Eb = 0.3305 kWh is the energy consumed by our Nissan Leaf battery and Em = 0.2599 kWh is the en- ergy consumed by our PMSM, then the energy efficiency, η = Em/Eb = 78.63 %. Therefore, our electric vehicle is approximately 19% Hungarian Journal of Industry and Chemistry ELECTRIC VEHICLE MODELLING AND SIMULATION 45 Table 7: PI controller parameters Gains Performance Indices Kp Ki IAE ISE ITAE ITSE 60 60 4.1141 × 10−6 1.0935 × 10−12 2.8676 × 10−5 5.1996 × 10−12 60 40 1.353 × 10−4 1.0372 × 10−9 0.0011 6.3622 × 10−9 5 20 7.9202 × 10−10 3.1365 × 10−20 7.9202 × 10−9 3.1363 × 10−19 more efficient than conventional vehicles. However, Du et al. determined that an EV can be up to approximately 15% more efficient than conventional vehicles driven by ICEs in terms of fuel consumption (fuel economy) [24]. The speed of the vehicle when tested by a ramp input is depicted in Fig. 22. In this case, the vehicle accelerates from rest to a final speed of 5.55 m/s (20 km/hr). The PMSM implemented as a result of a ramp test signal is presented in Fig. 23 below. The consumption of vehicles was less during this test to prove the perfor- mance of our motor and electric vehicle under these ideal conditions. 7.1 Settings of the PI controller based on the performance matrix The settings of the PI controller based on the performance matrix obtained by conducting several experiments until optimized gains had been achieved is presented in Table Figure 22: Speed of the vehicle due to a ramp input Figure 23: Speed of the motor due to a ramp input 7. It can be seen from this table that minimum values of these performance indices were required to obtain these optimal gains. Moreover, the system is stable and the con- troller somewhat reliable. However, a future study may propose intelligent tuning techniques to achieve a realis- tic level of performance by accurately taking into consid- eration the model uncertainties and parameter variations due to operating conditions, e.g., temperature and humid- ity amongst other considerable factors. 8. Conclusions The modelling and simulation of an electric powertrain based on the PMSM Motor of a light commercial vehicle, a VW Crafter manufactured in 2020, has been presented. A PI-based classical control algorithm of the outer con- trol loop of the PMSM was used to control the speed of the vehicle and motor. PMSM-based electrical-rear-wheel-driven traction of a VW Crafter manufactured in 2020 was modelled. This motor was chosen due to its high degree of efficiency over other electrical traction machines. This motor propelled the vehicle during its motoring action and while recharg- ing the battery when acting as a generator. A Nissan Leaf battery with a rating of 360 V manu- factured in 2011 and an energy supply to the system of 24 kWh was used. The battery was used to supply energy to propel the vehicle. This research determined that the fuel consumption decreased by a significant percentage by replacing conventional vehicles driven by ICEs with electric vehicles. Acknowledgments This research work is supported by TKP2020-NKA-04. Project no. TKP2020-NKA-04 has been implemented with the support provided from the National Research, Development and Innovation Fund of Hungary, financed under the 2020-4.1.1-TKP2020 funding scheme. REFERENCES [1] Ehsani, M.; Gao, Y.; Longo, S.; Ebrahimi, K.M.: Modern Electric, Hybrid Electric, and Fuel Cell Vehicles, Third ed. (CRC Press) 2018 DOI: 10.1201/9780429504884 [2] Un-Noor, F.; Padmanaban, S.; Mihet-Popa, L.; Mol- lah M. N.; Hossain, E.: A comprehensive study of key electric vehicle (EV) components, tech- nologies, challenges, impacts, and future direction 49(1) pp. 37–46 (2021) https://doi.org/10.1201/9780429504884 https://doi.org/10.1201/9780429504884 46 BABANGIDA AND SZEMES of development, Energies, 2017, 10(8), 1–82 DOI: 10.3390/en10081217 [3] Mohd, T. A. T.; Hassan, M. K.; Aziz, WMK. A.: Mathematical modeling and simulation of an elec- tric vehicle, J. Mech. Eng. & Sci., 2015, 8, 1312– 1321, DOI: 10.15282/jmes.8.2015.6.0128 [4] Shariff, S. M.; Iqbal, D.; Saad Alam, M.; Ahmad, F.: A state of the art review of electric vehicle to grid (V2G) technology, IOP Conference Series: Materi- als Science and Engineering, 2019, 561(1), 012103 DOI: 10.1088/1757-899X/561/1/012103 [5] Cioroianu, C. C.; Marinescu, D. G.; Iorga, A.; Sibiceanu, A. R.: Simulation of an electric vehi- cle model on the new WLTC test cycle using AVL CRUISE software IOP Conference Series: Materi- als Science and Engineering, 2017, 252(1), 012060 DOI: 10.1088/1757-899X/252/1/012060 [6] Wahono, B.; Nur, A.; Santoso, W. B.; Praptijanto, A.: A comparison study of range-extended engines for electric vehicle based on vehicle simulator, J. Mech. Eng. Sci., 2016, 10(1), 1803–1816 DOI: 10.15282/jmes.10.1.2016.5.0173 [7] Marmaras, C.; Xydas, E.; Cipcigan, L.: Simula- tion of electric vehicle driver behavior in road transport and electric power networks, Transp. Res. Part C Emerg. Technol., 2017, 80, 239–256 DOI: 10.1016/j.trc.2017.05.004 [8] Fodor, D.; Enisz, K.: Vehicle dynamics based ABS ECU verification on real-time hardware-in-the-loop simulator International Power Electronics and Mo- tion Control Conference and Exposition, 2014, 1247–1251 DOI: 10.1109/EPEPEMC.2014.6980683. [9] Aracil, C.; Sziebig, G.; Korondi, P.; Oh, S.; Tan, Z.; Ruderman, M.; He, W.; Ding, L.; Luo, H.; Yin, S.: Toward smart systems: Their sensing and control in industrial electronics and applications, IEEE Indus- trial Electronics Magazine, 2021, 15(1), 104–114 DOI: 10.1109/MIE.2020.3042171. [10] Abu Bakar, S. A.; Muhamad Said, M. F.; Aziz, A. A.: Ride comfort performance evaluations on electric vehicle conversion via simulations, ARPN J. Eng. Appl. Sci., 2015, 10(17), 7794–7798 http://www.arpnjournals.com/jeas [11] Ito, I.: Battery Electric Vehicle Model in Simscape, 2021 https://github.com [12] MathWorks Student Competitions Team (2020). Student Competitions - Physical Modeling Train- ingm, MATLAB Central File Exchange, 2020 https://www.mathworks.com [13] Yang, S.; Lu, Y.; Li, S.: An overview on vehicle dy- namics, Int. J. Dyn. Control, 2013, 1, 385–395 DOI: 10.1007/s40435-013-0032-y [14] Saleem, A.; Liu, N.; Junjie, H.; Iqbal, A.; Hayyat, M. A.; Mateen, M.: Modelling of an elec- tric vehicle for tractive force calculation along with factors affecting the total tractive power and energy demand, 2020 3rd International Con- ference on Computing, Mathematics and Engi- neering Technologies (iCoMET), 2020, 1–5 DOI: 10.1109/iCoMET48670.2020.9073845 [15] Dogruer, T.; Tan, N.: Design of PI controller, using optimization method in fractional order control sys- tems, IFAC-PapersOnLine, 2018, 51(4), 841–846 DOI: 10.1016/j.ifacol.2018.06.124 [16] Tan, W.; Marquez, H. J.; Chen, T.: Perfor- mance assessment of PID controllers, Con- trol Intell. Syst., 2004, 32(3), 158–166 DOI: 10.2316/journal.201.2004.3.201-1309. [17] Virani, V.; Arya, S.; Baria, J.: Modelling and con- trol of PMSM drive-by field oriented control for HEV, Advances in Power Generation from Renew- able Energy Sources (APGRES), 2019, 1–11 DOI: 10.2139/ssrn.3442515. [18] Espina, J.; Arias, A.; Balcells, J.; Ortega, C.: Speed anti-windup PI strategies review for field-oriented control of permanent magnet synchronous ma- chines, Compatibility and Power Electronics, 2009, 279–285 DOI: 10.1109/CPE.2009.5156047. [19] Yang, J.; Chen, W. H.; Li, S.; Guo, L.; Yan, Y.: Disturbance/uncertainty estimation and attenuation techniques in PMSM drives - A survey, IEEE Trans. Ind. Electron., 2017, 64(4), 3273–3285 DOI: 10.1109/TIE.2016.2583412. [20] Han, W.: Simulation model development of elec- tric motor and controller, Master’s thesis, Chalmers University of Technology, Gothenburg, Sweden, 2017 https://publications.lib.chalmers.se [21] Lina, W.; Kun, X.; De Lillo, L.; Empringham, L.; Wheeler, P.: PI controller relay auto-tuning using delay and phase margin in PMSM drives, Chi- nese J. Aeronaut., 2014, 27(6), 1527–1537, DOI: 10.1016/j.cja.2014.10.019 [22] Emadi, A.: Advanced electric drive vehicles, (CRC Press, Tailor & Francis Group, London–New York), 2015 ISBN: 978-1-4665-9770-9 [23] Al azze, Q.: Field-oriented control of permanent magnet synchronous motors based on DSP con- troller, Master’s thesis, Southern Illinois University, Edwardsville, USA, 2014 [24] Du, G.; Cao, W.; Hu, S.; Lin, Z.; Yuan, T.: Design and assessment of an electric vehicle powertrain model based on real-world driving and charging cycles,IEEE Trans. Vehicular Tech., 2019, 68(2), 1178–1187 DOI: 10.1109/TVT.2018.2884812 [25] Pacejka, H.: Tire and vehicle dynamics, (Butterworth-Heinemann, Elsevier Ltd., Oxford, UK), 2005 ISBN: 978-0-0809-7017-2 [26] Szalay, I.; Fodor, D.; Kohlrusz, G.: Modeling of slotless surface-mounted PM synchronous motor for sensorless applications, 2014 IEEE Interna- tional Electric Vehicle Conference (IEVC), 2014, 1– 5 DOI: 10.1109/IEVC.2014.7056198 Hungarian Journal of Industry and Chemistry https://doi.org/10.3390/en10081217 https://doi.org/10.3390/en10081217 https://doi.org/10.15282/jmes.8.2015.6.0128 https://doi.org/10.1088/1757-899X/561/1/012103 https://doi.org/10.1088/1757-899X/252/1/012060 https://doi.org/10.15282/jmes.10.1.2016.5.0173 https://doi.org/10.15282/jmes.10.1.2016.5.0173 https://doi.org/10.1016/j.trc.2017.05.004 https://doi.org/10.1016/j.trc.2017.05.004 https://doi.org/10.1109/EPEPEMC.2014.6980683 https://doi.org/10.1109/MIE.2020.3042171 http://www.arpnjournals.com/jeas/research_papers/rp_2015/jeas_0915_2643.pdf https://github.com/mathworks/Simscape-Battery-Electric-Vehicle-Model/releases/tag/1.0.0 https://www.mathworks.com/matlabcentral/fileexchange/52631-student-competitions-physical-modeling-training https://doi.org/10.1007/s40435-013-0032-y https://doi.org/10.1007/s40435-013-0032-y https://doi.org/10.1109/iCoMET48670.2020.9073845 https://doi.org/10.1109/iCoMET48670.2020.9073845 https://doi.org/10.1016/j.ifacol.2018.06.124 https://doi.org/10.2316/journal.201.2004.3.201-1309 https://doi.org/10.2316/journal.201.2004.3.201-1309 https://doi.org/10.2139/ssrn.3442515 https://doi.org/10.2139/ssrn.3442515 https://doi.org/10.1109/CPE.2009.5156047 https://doi.org/10.1109/TIE.2016.2583412 https://doi.org/10.1109/TIE.2016.2583412 https://publications.lib.chalmers.se/records/fulltext/253368/253368.pdf https://doi.org/10.1016/j.cja.2014.10.019 https://doi.org/10.1016/j.cja.2014.10.019 https://doi.org/10.1109/TVT.2018.2884812 https://doi.org/10.1109/IEVC.2014.7056198 Introduction Electric powertrain Electrical parts Mechanical parts Simulation of the dynamic system of the electric powertrain Vehicle transmission system Vehicle modelling Tyre Dynamics using the Magic Formula Speed controller Integral performance criteria Permanent magnet synchronous motor PMSM Mathematical Model Simplified PMSM model Nissan Leaf Battery Battery modelling Battery parameters Matlab model Results Settings of the PI controller based on the performance matrix Conclusions