Open access peer-reviewed chapter

Powerful Multilevel Simulation Tool for HiL Analysis of Urban Electric Vehicle’s Propulsion Systems

Written By

Raul Octavian Nemeș, Mircea Ruba, Sorina Maria Ciornei and Raluca Maria Raia

Submitted: May 4th, 2021 Reviewed: May 24th, 2021 Published: June 24th, 2021

DOI: 10.5772/intechopen.98532

Chapter metrics overview

92 Chapter Downloads

View Full Metrics

Abstract

The general focus of the proposed chapter is to describe a complex yet transparent solution for advanced simulation analysis of urban electric vehicles propulsion unit. As general rule, precise and realistic results are obtained only when performing real-time simulations, engaging dedicated software for such applications. Hence, simulation of an electric vehicle as a complete solution can become rather difficult. The authors targeted advanced analysis of the propulsion unit, including the motor, the battery, the power converter, and its control. These are designed using multilevel models in Matlab/Simulink, referring to different complexity levels of each assembly. Another feature of the models is their organization, based on Energetic Macroscopic Representation (EMR), this easing the process of inter-connecting models correctly. Nevertheless, the mechanical, aerodynamical and road profile details are included using Amesim Software. All the simulations are performed on a real-time target, using a National Instruments PXIe embedded controller. The latter runs NI VeriStand software, allowing real-time communication between Amesim and Simulink offering in the same time possibility to read/write analog/digital IOs for external communication. This feature in fact is used when passing from modeling to Hardware in the Loop (HIL) analysis, replacing the simulated assembly with the actual one.

Keywords

  • real-time
  • Hardware in the Loop
  • EMR
  • multilevel
  • electric vehicle

1. Introduction

Electrification of urban transportation services, from high power solutions, such as busses and trains, to low power ones, such as small cargo delivery vehicles became highly important. The pressure of developing such solutions is continuously increasing due to the highly crowded urban scenarios in which merchandise and goods must be delivered as fast as possible. Investigating, researching and creating such solutions challenges many universities and OEM engineers. Nowadays, taking advantage of the achievements in the field of computer aided design, any serious study requires advanced simulations, fusing measured data, estimated and computed ones into powerful computation models, able to replicate behavior of the actual hardware. It is important to mention that at a certain level of complexity, studies can become rather difficult by means of large required computation time, powerful computers to perform and many model’s implementations time.

In order to achieve high precision and accurate phenomenon simulation, one must use several software interconnected, each performing analysis in their field of expertise. For example, when simulation of an urban electric vehicle (UEV), all the mechanical, aerodynamical and road conditioning factors can be simulated in Amesim Software, while all the power electronics, electrical drive assemblies and their advanced control can be simulated in Matlab/Simulink. Such a solution ensures high fidelity and correct analysis approach of each assembly of the UEV. However, a co-simulation of this magnitude becomes rather complicated, time and resource demanding, returning a non-real-time behavioral screening. Hence, one solution is to engage a third software, National Instruments VeriStand. This was created to integrate exported model from many simulation software and create one complex simulation solution. Another hands-on advantage is that using a PXIe embedded controller, VeriStand can upload the model and perform the simulation in real-time, allowing the user to observe the system behavior just like in the real life. The simulation time can be from seconds elongated to tens of hours. The data recording can be streamed and decimated as desired. Such a solution will be presented in this chapter.

Another important feature of this tool is given by the real-time target’s nature. This allows external communication via analog and digital I/O s simplifying and facilitating the transition from model-based analysis to hardware in the loop (HIL) based analysis.

In order to offer the possibility of selecting the accuracy and the magnitude of the model, the authors considered developing multilevel complexity model, all organized in Energetic Macroscopic Representation (EMR). The latter ensures the user with the correct understanding of the system’s power flow, creates a direct path when designing the control loops and facilitates the model exchange when passing from one level of complexity to another. It has to be mentioned that all the models presented in the following sections are designed in EMR philosophy. However, the interest of this chapter is not focused towards this special organization method, hence the reader is encouraged to lecture references [1, 2, 3].

Advertisement

2. Multilevel modeling of electric vehicle’s propulsion unit

The multilevel approach for developing simulation models for an electric vehicle’s propulsion system must be broken down into its main assemblies: the power supply (batteries), the electronic converter, the electrical machine and the general system controller (electronic control unit-ECU). To have a comprehensive analysis, each of these assemblies will be addressed separately in the following sub-sections. For each model, the corresponding energetic macroscopic representation (EMR) block will be presented, separately from its content. The latter will be in fact the modeling approach of each assembly. EMR is an energy flow based graphical organization philosophy easing the building process of complex electromechanical systems. The connection between blocks is always based on action-reaction principle according to the physical causality. In Figure 1, the basic elements used in EMR organization are depicted. These will be used in the following sections as some assemblies, for example, an electrical machine will include more than one such block in their graphical description.

Figure 1.

Main EMR pictograms.

Before diving in the presentation of the multilevel models, it is lucrative to explain the philosophy of building a control scheme using EMR. This is necessary because, in the present chapter, the control units will be also under discussion and their multilevel attribute will come from the level of compromises chosen when building its EMR representation. The control loop can be simply generated by inverting all the model blocks. Hence, the blue pictograms (the control) are a reversed model of the orange ones (the model). The main difference is that in the control loop, only the action path is reversed, while the reaction is neglected, as the latter one is necessary only for the model description and not for the control. In the above indicated references, more details about building EMR control loops are presented, hence the reader in encouraged to visit those.

2.1 Multilevel modeling of batteries

When it comes to model a power supply for an electric vehicle, one can choose between ideal or non-ideal approaches. The latter category can be divided into less or more complex mathematical description. In general, besides the ideal approach the non-ideal one include first and second order battery modeling [4, 5]. A source or a sink in EMR is represented by a green oval, outputting voltage, and inputting current or vice versa.

In Figure 2 the pictogram used for the battery as electrical source in a simulation is depicted. The content of this block however, can be adjusted to the designer’s needs. In Figure 2(a)(c), the ideal battery model is depicted (where the voltage is imposed invariant of any changes in the system), the first and second order models (where the battery voltage is dependent on the sourced/sinked current and the state of charge (SOC)). The latter models are called first and second order ones due to the number of parallel RC branches that compound their circuit.

Figure 2.

First (a) and second (b) order equivalent battery circuits and their EMR pictogram (c).

Vbat=Vocu1RsIbat;du1dt=1RpCpu1+1CpIbatE1

The mathematical translation of the circuits is detailed in Eq. (1) and Eq. (2). It is important to mention that the circuits parameter’s values for such models are identified using laboratory experiments and actual battery cells. Such approaches were carried out by the authors on a LG (LGABD11865) battery with a rated capacity of 3000 mAh, 3.75 V rated, 4.2 V maximum over charge voltage, 2.7 V minimum discharge voltage, 0.5–1 A charging current and 0.2–0.5 A discharging current.

Vbat=Vocu1u2RsIbat
du1dt=1R1C1u1+1C1Ibat
du2dt=1R2C2u1+1C2IbatE2

The complete process described step-by-step to identify the parameters of the battery cell needed to run the correct model is presented in detail in [4]. Even more, not only the experimental process itself is presented, but also very important post processing steps mandatory to reach high accuracy results.

The open circuit voltage (VOC) is identical for both models, being measured for the same conditions. The recorded function of the SOC is depicted in Figure 3.

Figure 3.

The open circuit cell voltage vs. SOC.

To help the reader emphasize an image about the precision of these models, a comparative analysis was carried out, superimposing their results for the same testing conditions. After the parameter’s identification process [5, 6, 7, 8], both models presented in Eq. (1) and Eq. (2) were built in Matlab/Simulink, as depicted in Figure 4-left.

Figure 4.

The Matlab/Simulink cell models (left) and the comparative analysis (right).

The considered test was to apply equal length current pulses with equal length relaxation time between them, discharging the cell from 100–10%. Initially, this test was applied to the actual cell, the results being also used to identify the necessary data for the simulation models. Imposing the same current variation to the two models, and comparing their voltage calculations with the measured one, one can conclude that the second order model reaches higher accuracy than the first order one. In the zoomed plot from Figure 4-right, it can be observed that the second order and measured voltages are in best agreement, while the first order’s one marks a slight deviation.

To quantify this difference to a certain palpable value, the root mean square error (RMSE) is computed using Eq. (3) and the instantaneous error between the two signals.

RMSE=1T0TVmeasuredtVmodelt2dtE3

The comparative results of the RMSE and the instantaneous error are depicted in Figure 5 (left and right). In both cases one can observe that the second order model yells lower disturbance compared to the first order one. More, in Figure 5-left shows that during the ongoing discharging process, the RMSE for the first order model increases gradually faster than the one of the second order one.

Figure 5.

The RMSE (left) and the instantaneous error (right) of the battery models.

To conclude this section, one can understand that increasing the rank of the battery model, the gap between real and simulated results is mitigated. The cost of this is the increased complexity of models, time-consuming parameters identification and the necessity of more data post processing operations. On the contrary, the first order model is quite simple to design, its parameters are simple to identify with classical calculations and there is no need for values post processing. If the compromise of accepting a small error in the battery model with the gain of simplicity, the first order model is the hands-on solution.

2.2 Multilevel modeling of power electronics’ losses

The general approach when dealing simulation approaches of electronic converters regard conventional ideal ones, with or without switching, so called smooth or switched models. However, none of these consider the losses that occur in a non-ideal converter. The largest losses are due to the switching components (power transistors and/or diodes), while the rest of the losses inside a converter can be neglected (those on the driver and measuring circuits). One reason why the losses are generally neglected overall is because their actual measurement is rather difficult and requires dedicated equipment. In doing so, one simple method is to consider a general efficiency of the converter and add it into its calculations. So, already some straightforward solutions for modeling electronic converters were mentioned, solutions that are available in nearly any reference the reader would search. To add value to the research and to the present chapter, the authors introduce 2 different approaches to compute the losses that occur on one MOSFET transistor. Of course, these can be then extrapolated towards the total number of power switches.

In Figure 6-bottom the EMR pictogram used for any electronic converter is presented. It inputs the DC link voltage (Vdc), the load currents (i_load) and the duty cycle (m_s). After the internal calculations, the converter outputs the resulting voltage (u) to be supplied to the load and its output currents (i_conv). The inputs and outputs do not change with the chosen method of internal calculations. Either is a smooth or switched model, these remain the same and valid even if we consider DC or AC converters.

Figure 6.

Smooth (a) and switched (b) converter architectures and their EMR pictogram (c).

In addition to the above depicted models, this block can include also the calculus for the converter losses. Even if in the standard pictogram there are no input/output ports considered for this, the information, being additional, can be manipulated to/from the block via wireless labels.

Taking into consideration of these losses is mandatory when one is interested in reaching optimal converter architectures [9, 10]. The following calculus can be adapted to any power transistor, the authors using the IRFP4668PBF MOSFET, having its specifications listed in Table 1.

ParameterSymbolValueUnit
Forward voltageV120V
Drain to source resistanceRds0.008Ω
Switching Frequencyfsw10000Hz
Output CapacitanceCoss810e-12F
Turn off timeToff64e-9s
Turn on timeTon41e-9s
Forward voltage dropVf1.3V
Reverse Leakage currentIr0.00025A
Reverse Recovery TimeTrr135e-9s
Max Recovery CurrentIrrm8.7A
Gate ResistanceRg11Ω
Plateau VoltageUplateau6.1V
Gate CapacitanceCGD1200e-12F
Gate CapacitanceCGD21200e-12F
Reverse recovery chargeQrr700e-9C
Rise Timetr105e-9s
Fall timetf74e-9s
Reverse Recovery CurrentIDoff8.7A
Average diode currentIFav20A
Diode on-state zero-current voltageud00.5V

Table 1.

The main specifications for IRFP4668PbF MOSFET.

The losses computed in any active switching electronic component can be broken down into: conduction losses, switching losses and blocking (leakage) losses, the later normally being neglected, as expressed in Eq. (4).

PlossT=PcondctionT+PswitchingTE4

As mentioned before, two approaches will be presented in detail. The first one computes the power losses during conduction and switching the MOSFET while the second one computes losses during switch-on and switch-off transients.

In the first approach, the conduction and switching losses are computed based on Eqs. (5) and (6), with regards to the semiconductor’s parameters and its operation. The latter is described via the duty cycle (D), the switching frequency (fsw) and the transistor’s output capacitance (Coss)

PconductionT=Id2RdsDE5
PswitchingT=VId2Ton+Tofffsw+CossV2fswE6

On the other hand, the losses occurring in the diode attached to the transistor are broken down into conduction, reverse and switching ones. The total losses due to the diode’s operation are:

PlossD=PcondctionD+PreverseD+PswitchingDE7

Each of the eq.’s (7) components are formulated in expressions (8)(10), with parameters valued and explained in Table 1:

Pcondction=IfVfDE8
Preverse=IrV1DE9
Pswitching=TrrIrrmVrfsw2E10

This model is quite simple and represents the hands-on solution for quantifying the converter losses without engaging too many mathematical tools.

The more complex model, the conduction losses are computed function of the on-state resistance (Rds) and the current passing through the component (id):

Pcond=Rds2idE11

The next component of the total losses, the switching one, is broken down into two states, according to the on and off transients. There are different opinions and mathematical interpretations on how to compute the voltage fall-time (tfu) to gain accuracy. Details are presented in Ref. [10]. Engaging a straight forward approach, the fall-time is considered to be composed of two parts in accordance with the variation of the gate capacitance:

tfu=tfu1+tfu22
tfu1=UDDRDSIdRGCGD1UdrvUplateauE12
tfu2=UDDRDSIdRGCGD2UdrvUplateau

In Eq. (12), the terms (Udrv) and (UDD) refer to the driver voltage that excites the transistor’s gate and the converter’s DC link voltage, respectively. The rest of the terms are to be found in Table 1.

The complementary model to compute the switch-off losses differ just by subtracting from Eq. (12) the driver voltage.

tru=tru1+tru22
tru1=UDDRDSIdRGCGD1UplateauE13
tru2=UDDRDSIdRGCGD2Uplateau

The required energy to turn on and off the transistor is calculated in Eq. (14) and (15).

EonT=UDDIDontru+tfu2+QrrUDDE14
EoffT=UDDIDofftru+tf2E15

With the above computed ingredients, one can calculate the power losses on the transistor, Eq. (16), and diode, Eq. (17).

PswT=EonT+EoffTfswE16
PswD=EonD+EoffDfswEonDfswE17

The losses in the MOSFET and the free-wheeling diode are computed as the sum of the conduction and switching losses.

PTotalT=PconductionT+PswitchingT
PTotalT=RDSonI2Drms+EonM+EoffMfswE18
PTotalD=PConductionD+PswitchingD
PTotalD=uD0IFav+RDI2Frms+EonDfswE19

It has to be mentioned that during operation, in one transistor or diode, the currents can be positive or negative. In an analytical approach, this would lead to decreased losses, hence, to correctly calculate those, the absolute current’s value is used.

Figure 7.

The instantaneous converter losses for the simple and complex models.

Comparing the results of the two approaches (Figure 7) one can conclude that the complex model reaches higher calculated values. This is more realistic as it contains more loss sources compared to the simple model. However, the simple model proves to offer a quite good quantification of the losses despite its ease in approach, reaching satisfactory results.

2.3 Multilevel modeling of permanent magnet synchronous machines

The general approach when modeling a permanent magnet synchronous machine (PMSM) is to use its mathematical representation in dq0 reference frame, fixed to the rotor.

vsdRsisdesd=ddtψds
vsqRsisqesq=ddtψqsE20
Tpmsm=32pΨdsiqsΨqsids
esd=ωΨqs,esq=ωΨds
Ψds=Ψmd+LdIds+LdqIsq
Ψqs=Ψmqd+LqIqs+LqdIdsE21
JdΩdt=TpmsmTloadE22

The variables in Eqs. (20) and (21) are: the phase resistance (Rs), the direct and quadrature inductances (Ld and Lq) and the cross-coupling inductances (Ldq and Lqd). The number of pole pairs is denoted with (p), while the rotational speed is marked by (ω) and the friction coefficient in marked by (B). The dq currents (Ids and Iqs) are computed by integration of the voltage expressions.

For a PMSM, the EMR model is constructed by coupling 3 pictograms, one accumulation element (for the windings) described by Eq. (20), one electromechanical conversion (for the torque generation), described by Eq. (21) and another accumulation element (for the mechanical equilibrium), described by Eq. (22), (visit Table 1 for clarifications). Such a model is depicted in Figure 8.

Figure 8.

The EMR for a PMSM.

Returning to Eq. (21), the dq magnetic flux equations give the opportunity to choose the level of modeling complexity. In this section, 3 different such levels will be detailed: based on constant dq inductances, based on variable dq inductances and a flux-linkage based model. The latter is the most complex and precise one.

To have a comprehensive analysis, like for the battery and the power inverter, for the PMSM the model’s parameters are listed in Table 2.

Electric and geometric parameters
DimensionValueUnit
1Pole pairs3
2Rotor speed2000[rot/min]
3Rated torque12[Nm]
4Rated current per phase21[A]
5Mechanical power2500[W]
6DC Voltage link120[V]

Table 2.

The PMSM under analysis.

The first level of complexity consists in using constant values for the d and q axes inductances. The identification of these is a very simple process. It can be carried out in finite element analysis (FEA) or on the real machine (or in both environments). It requires only fixing the rotor in 0 and 90 electrical degrees positions and applying in both one controlled voltage pulse. The time dependent slope of the current will return after some simple calculations the global machine inductance (L). Having this value, the d and q components are easy to be computed, using Eq. (23). Details about the identification process can be selected from references [11, 12].

Ld=23Lθel=0,Lq=23Lθel=90E23

The second level model regards detailed information about the machine under testing. The latter needs to be simulated in a FEA environment in order to compute the dq inductance’s variations function of the machine current. This is generally ranged from ±10 × In, where In is the machine’s rated current. The method itself that allows preliminary identification of the 3-phase magnetic fluxes is the so called “frozen permeability” approach. Applying it demands to set to 0 the permanent magnet’s remanent flux density, freezing its permeability. More details about this method are to be found in Ref.s [12, 13]. The resulted phase magnetic fluxes will be used to compute the dq inductances versus the machine current, applying Eq. (24).

Ψd=23Ψa12Ψb12Ψc,Ψq=13ΨbΨc.
Ld=Ψdid,Lqd=Ψqid,Lq=Ψqiq,Ldq=ΨdiqE24

The outcome of these calculations is stored in look up tables (LUT) and introduced in the machine model. It is clear that such a procedure is impossible to be implemented in a real testbench or it would require complex and expensive hardware. One can observe that using the calculus from Eq. (23), the direct as well as the cross-coupling inductances are computed. Their variation function of the ranged current is depicted in Figure 9.

Figure 9.

The direct and cross-coupling inductances calculated with frozen permeabilities.

The last and altogether the most complex model is based on using instead of inductances, the variation of the d and q currents function of the d and q magnetic fluxes. These are 2D LUTs as depicted in Figure 10.

Figure 10.

The 2D LUTs of idd, Ψq) and iqd, Ψq).

It is important to mention that building the 2D LUTs from Figure 10 is a quite complex and challenging process. Here to the designer needs to have the FEA model of the machine under study. Several simulations need to be accomplished, fetching initially the Ψd (id, iq) and Ψq (id, iq). From this information mathematical inversion methods need to be applied to get idd, Ψq) and iqd, Ψq). The complete process is explained in detail in [14].

In Figure 11 a simple comparative analysis of the results obtained when simulating all the 3 models for the same conditions are depicted. Focusing the attention on the variation of the q axis magnetic flux, it can be noticed that flux linkage and frozen permeability models do saturate, while the one based on constant inductances reaches large values. The latter is ideal; hence the machine will never saturate. In this condition, one can perform simulations within rated operation. However, when attempting to go in overload region the results will not reflect the real machine operation.

Figure 11.

The PMSM torque and Ψq during benchmarking the simulation levels.

The flux linkage and frozen permeabilities models are behaving quite in the same manner. However, the complexity of the flux linkage model is high, both in data post processing and design. The ultimate conclusion is that using the frozen permeabilities model will lead to quite satisfactory results reflecting closely the actual machine behavior.

2.4 Designing the control loop

Energetic Macroscopic Representation is one organization concept that besides helping the designer to build simulation programs that reflect the actual power flow of a real testbench, it also proves to be a very lucrative solution when it comes to building control loops. To do so, one has just to inverse the simulation model and keep only the action path, discarding the reaction one.

In Figure 12 the control loop can be observed being highlighted in light blue. Its inversed character compared to the model (in orange) is noticeable and each block contains the mathematical description of the inversed action.

Figure 12.

The complete PMSM simulation model in EMR organization.

The reversed accumulation block contains a regulator that outputs the reference torque to be developed function of the error between the actual and the imposed speed (see Eq. (25)).

ξnPMSM=nPMSM_refnPMSM
ξnPMSMPIkp_nki_nTPMSME25

The inversed block without controller computes the reference q axis current function of the refence torque above expressed. In the same time, the d axis current is imposed to 0, this for torque maximization purposes.

isq_ref=TPMSM32pΨPM+LdLqisd
isd_ref=0E26

This reference d and q currents are compared with the actual ones. The resulting errors are sent to regulators that output the reference d and q voltages. This process is carried out in the third control block (see Figure 12 from right to left).

ξiqd=isqd_refisqd
ξiqdPIkp_iki_iusqd
usqd_est=usqd+esqd_estE27

The resulted d and q voltage refences are converter to 3-phase quantities and sent to the last control block, that computes the duty-cycle (m) for each phase. The duty-cycle is sent to the PWM generator (in orange) closing the model’s loop. More details about generating the process model are indicated in Ref. [15].

Advertisement

3. Real-time model analysis

To develop an analysis tool able to replicate as close as possible the behavior of a real electric urban vehicle, it is clear that using complex simulation developments for the electrical assemblies will not be enough. However, building the mechanical models using advanced mathematical tools is highly time consuming and would require a lot of expertise and education in this direction. Hence, a wise solution is to engage a software created for such simulations, like Simcenter Amesim from Siemens.

In order to highlight the previously mentioned application with a practical example, an urban vehicle’s mechanical model is developed. The vehicle is an electrical tricycle used for cargo delivery purposes. Its main specifications are presented in Table 3.

ParameterSymbolValue
Maximum massm250 kg
Maximum speedv10 m/s
Wheel radiusRw12″
Maximum road slopeα15%
Rolling resistancef0.017
Front areaA1.2 m2
Air drag coefficientCs0.36
Gear ratioGr5

Table 3.

Tricycle specifications.

In Figure 13, the Amesim model is depicted, observing that none electrical assemblies are included in the model. Those are modeled in Matlab/Simulink as it will be presented in the following sections. The link between the two models will be performed using National Instruments VeriStand software, able to create an integrated project that runs a co-simulation on a real-time processor for the complete vehicle.

Figure 13.

The mechanical assemblies’ Amesim model.

The Matlab/Simulink model for the electrical assemblies is built putting together all the components detailed in Section 2 based on the hierarchy depicted in Figure 12. It has to be mentioned that the pictograms cannot be exactly replicated in Simulink as the inputs are on the left and the outputs always on the right. However, the organization is the same as it can be observed in Figure 14.

Figure 14.

The urban vehicle’s EMR in Matlab/Simulink.

Having all the models above presented, the goal is to create a virtual simulation platform that allows simple transition hardware in the loop (HiL) testing. The latter is the ultimate goal of testing any electromechanical system. It allows users to test certain physical assemblies while the rest of the system is virtual, running as simulation on a real-time target. This must ensure enough computation power and speed to cope with the demands of the actual hardware. National Instruments have in their portfolio hardware (NI PXIe embedded controller) and software (NI VeriStand) that are able to integrate into one real-time simulation all the above detailed models. The PXI computer, via its onboard field programable gate array (FPGA) enables the connection of the virtual model with the testbench using analog and digital IOs.

In Figure 15, the main components of the analysis platform are depicted. One can observe that VeriStand software that integrates the Amesim and Simulink models runs on the PXI and using the FPGA channels communicates with the actual hardware. The user has the possibility to run the virtual model and in parallel the actual hardware in a bidirectional communication network, performing as an entire entity. The real-time processor also features the possibility to run only the simulation of the entire vehicle, without any hardware connected to it. The difference between such a simulation and one running on a PC is that the latter would take tens of hours for a road cycle of 30 minutes while on the real-time computer it will take the elapsed 30 minutes.

Figure 15.

The platform’s hardware/software architecture.

Both scenarios will be presented in detail in the following section, comparing the results and discussing them.

Advertisement

4. System validation via urban scenario

Taking advantage of the flexibility of the above presented analysis platform one can perform a study using only simulation, hence a virtual vehicle or can combine virtual elements (mechanical assemblies) with real ones (the propulsion motor). For a robust simulation model, it is important to ensure that using it as reference will return results that mimic closely the actual hardware. In doing so, the platform depicted in Figure 15 was initially tested only for simulation, the complete system being entirely virtual and running on the PXIe under VeriStand software. The latter ensured the continuous communication of the Simulink model (running all electrical assemblies) with Amesim (running all mechanical systems).

It has to be mentioned that the chosen simulation used the most complex model assemblies from those presented in Section 2. In order to avoid redundancy, results for the rest of the levels are not presented in the chapter, however during the presentations in Section 2, the main differences between the models were already presented.

The second step was to keep in the virtual level the mechanical assemblies (in Amesim) and the battery model (in Simulink). The rest of the virtual electrical assemblies (the inverter, the motor and the controller) were replaced with their laboratory homologs. These remained however connected to the simulator via the analog/digital IOs presented earlier.

The outcome of this comparative analysis was more than satisfactory. Firstly, the reference torque was plotted versus the simulated and the measured ones. It can be observed that the agreement is very good proving the accuracy of the simulator. Knowing that in a PMSM the q current component is responsible for the torque production, for both measured and simulated cases this was recorded and depicted in Figure 16-right. The same conclusions as for the comparative analysis of the torque characteristics can be considered. In order to avoid redundancy, the d currents were not depicted as those values are forced to 0 at all times.

Figure 16.

The comparative analysis of the torque (right) and q current component (left).

The slight differences between the values plotted in Figure 16 are more due to measurements error and noise. No filtering was used what’s so ever in order to avoid any unnecessary postprocessing (Figure 17).

Figure 17.

The comparative analysis of the speed (right) and its zoomed section (left).

The slight differences between the values plotted in Figure 16 are more due to measurements error and noise. No filtering was used what’s so ever in order to avoid any possible influence over the quality results.

Smoother results can be reached if certain data sampling time or other such procedure is considered. However, the interest in this entire work was to prove that one can reach very accurate results when building responsible simulation programs. Also using a real-time processor can benefit the study with fast analysis and reliable results as well.

Advertisement

5. Conclusions

When it comes to develop simulation tools wise enough to reach the expectations of the fast-growing industry in the field of engineering, one needs to consider solutions that become hands-on. These must be flexible to changes, simple to implement and accurate when it comes to results. Using real-time processors to run these simulations offers reduced analysis time, accurate results and close to real study. The same processors when engaged, ease the transition from simulation to HiL testing, simply by replacing virtual components with real ones. Choosing between the complexity levels of the virtual models, allows the user to select the accuracy and the necessary time invested in the development. The system under analysis being compound of several assemblies, more complex models can be considered for those that are of interest, while the rest can be ideal ones.

In the present chapter the main assemblies of an electrical urban vehicle’s traction system are presented, offering the choice of complexity, mathematical description and EMR organization. The latter is introduced as graphical method for representing the elements of any simulation program by respecting the actual action/reaction physical and natural laws.

It has to be mentioned that the nature of the present chapter is more towards a review than of an academic lecture, hence the reader is encouraged to consult the indicated references that are guidance and complementary information to the one described in the previous pages.

Advertisement

Acknowledgments

This work was supported by a grant of the Ministry of Research, Innovation and Digitization, CNCS/CCCDI – UEFISCDI, project number PN-III-P2-2.1-PED-2019-4384, within PNCDI III.

This work was supported by a grant of the Ministry of Research, Innovation and Digitization, CNCS/CCCDI – UEFISCDI, project number PN-III-P1-1.1-TE-2019-0411, within PNCDI III.

References

  1. 1. A. Bouscayrol, R. Schoenfeld, G. Dauphin-Tanguy, G.-H. Geitner, X. Guillaud, A. Pennamen, and J. Hautier, “Different Energetic Descriptions for Electromechanical Systems,” in Proc. of the European Conference on Power Electronics and Applications (EPE), Sep. 2005, Dresden, Germany
  2. 2. K. Chen, A. Bouscayrol, and W. Lhomme, “Energetic macroscopic representation and inversion-based control: Application to an electric vehicle with an electrical differential,” Journal of Asian Electric Vehicles, vol. 6, no. 1, pp. 1097–1102, Jun. 2008
  3. 3. A. Bouscayrol, X. Guillaud, P. Delarue, and B. Lemaire-Semail, “Energetic macroscopic representation and inversion-based control illustrated on a wind-energyconversion system using hardware-in-the-loop simulation,” IEEE Transactions on Industrial Electronics, vol. 56, no. 12, pp. 4826–4835, Dec. 2009
  4. 4. Nemes, R.O., Maria Ciornei, S., Ruba, M., Martis, C., "Parameters Identification Using Experimental Measurements for Equivalent Circuit Lithium-Ion Cell Models", 11th International Symposium on Advanced Topics in Electrical Engineering, ATEE 2019; Bucharest; Romania; 28 March 2019 through 30 March 2019; DOI: 10.1109/ATEE.2019.8724878
  5. 5. Ruba, M., Ciornei, S., Nemeş, R., Martis, C.,"Detailed Design of Second Order Model of Lithium-Ion Battery Simulator Based on Experimental Measurements", 11th International Symposium on Advanced Topics in Electrical Engineering, ATEE 2019; Bucharest; Romania; 28 March 2019 through 30 March 2019; DOI: 10.1109/ATEE.2019.8724878
  6. 6. A. Nikolian, Firouz, Y. Gopalakrishnan, R. Timmermans, J. M. Timmermans, N. Omar, P. van den Bossche, J. van Mierlo, “Lithium ion batteriesDevelopment of advanced electrical equivalent circuit models for nickel manganese cobalt lithium-ion,”,Energies 2016
  7. 7. X. Zhang, J. Lu, S. Yuan, J. Yang, X. Zhou, “A novel method for identification of lithium-ion battery equivalent circuit model parameters considering electrochemical properties,” J. Power Sources 2017
  8. 8. M. Mathew, Q. H. Kong, J. McGrory, M. Fowler, Simulation of lithium ion battery replacement in a battery pack for application in electric vehicles, J. Power Sources 2017
  9. 9. D. Graovac and M. Pürschel, “MOSFET power losses calculation using the data- sheet Parameters,” 2006.https://www.infineon.com
  10. 10. S. Shuvo, E. Hossain, T. Islam, A. Akib, S. Padmanaban and M. Z. R. Khan, "Design and hardware implementation considerations of modified multilevel cascaded H-bridge inverter for photovoltaic system," in IEEE Access, vol. 7, pp. 16504-16524, 2019, doi: 10.1109/ACCESS.2019.2894757
  11. 11. D. E. Pinto, A. C. Pop, J. Kempkes, J. Gyselinck, “dq0-modeling of interior permanent-magnet synchronous Machines for High-Fidelity Model Order Reduction, Optim & Intl Aegean Conference on Electrical Machines and Power Electronics, May, 2017
  12. 12. D. Zarko, D. Ban, R. Klaric, “Finite element approach to calculation of parameters of an interior permanent magnet motor”, Automatika 46, 2005
  13. 13. J. Wu, J. Wang, C. Gan, Q. Sun and W. Kong, "Efficiency optimization of PMSM drives using field-circuit coupled FEM for EV/HEV applications," in IEEE Access, vol. 6, pp. 15192-15201, 2018, doi: 10.1109/ACCESS.2018.2813987
  14. 14. S. M. Ciornei, M. Ruba, R. O. Nemeş and C. Marţiş, "Multi-level models for a light electric vehicle propulsion system using EMR organisation," 2020 International Symposium on Power Electronics, Electrical Drives, Automation and Motion (SPEEDAM), 2020, pp. 357-362, doi: 10.1109/SPEEDAM48782.2020.9161855
  15. 15. M. Ruba, R. O. Nemes, S. M. Ciornei, C. Martis, A. Bouscayrol and H. Hedesiu, "Digital twin real-time FPGA implementation for light electric vehicle propulsion system using EMR organization," 2019 IEEE Vehicle Power and Propulsion Conference (VPPC), 2019, pp. 1-6, doi: 10.1109/VPPC46532.2019.8952428

Written By

Raul Octavian Nemeș, Mircea Ruba, Sorina Maria Ciornei and Raluca Maria Raia

Submitted: May 4th, 2021 Reviewed: May 24th, 2021 Published: June 24th, 2021