Step-by-Step Small-Signal Modeling and Control of a Light Hybrid Electric Vehicle Propulsion System

: This paper develops step-by-step a complete electric model of a light hybrid electric vehicle propulsion system. This model includes the vehicle mass, the radius and mass of the wheels, the aerodynamic proﬁle of the vehicle, the electric motor and the motor drive, among other elements. Each element of the model is represented by a set of equations, which lead to getting an equivalent electric circuit. Based on this model, the outer and inner loop compensators of the motor drive control circuit are designed to provide stability and a fast dynamic response to the system. To achieve this, the steady-state equations and the small-signal model of the equivalent electric circuit are also obtained. Furthermore, as these elements are the main load of the power distribution system of the fully electric and light hybrid electric vehicle, the input impedance model of the set composed of the input ﬁlter, the motor drive, the motor, and the vehicle is presented. This input impedance is especially useful to get the system stability of the entire power distribution system.


Introduction
The presence of electric vehicles (EV) in the automotive industry is growing rapidly [1]. According to the International Energy Agency, 125 million electric vehicles are expected to be manufactured by 2030 [2]. This growing market requires the optimization of the vehicle design and, therefore, an optimized propulsion system. If the system is too small, the driving range decreases; however, if the system is too big, the cost of the vehicle increases. Thus, the correct sizing is needed to properly design the vehicle.
Power sources such as fuel cells, batteries, and supercapacitors have to manage a power that depends on the power consumption of the vehicle. This power consumption is related to two key factors: the shape and weight of the vehicle and the driving profile. Figure 1 shows, as an example, the block diagram of a hybrid vehicle power distribution system and the main load, composed of the electric motor, the vehicle, the motor drive, and the filter. The vehicle is powered by an electric motor, which needs a motor drive to control its speed and a filter to diminish the propagation of current ripple effects towards the system. The primary energy source, e.g., a fuel cell, needs a DC-DC power converter that adapts and controls its power. The secondary energy source, e.g., a battery, can be connected either directly to the DC bus or via a DC-DC power converter.
The vehicle block diagram in Figure 1 shows the different elements in charge of transforming the power sources energy into force and linear speed in the vehicle. This transformation depends on many factors, such as the mass of the vehicle, the radius and the mass of the wheels, and the aerodynamic profile of the vehicle, among others. Besides, the terrain conditions have a considerable influence, The electric motor model is another key element. State of the art shows that some studies use DC motors [16,17] while other studies use AC motors, e.g., induction motors in [18][19][20], permanent magnet synchronous motors in [18,[21][22][23] and, finally, switched reluctance motors in [17,23,24]. The procedure used in each type of motor is similar, although the set of equations and control techniques could be different between them.
In addition, the electric motor needs a motor drive to control the desired vehicle speed. This motor drive must provide to the entire system good stability and dynamic response, which is achieved with the outer speed loop compensator and the inner current loop compensator, as shown in Figure 2.  The electric motor model is another key element. State of the art shows that some studies use DC motors [16,17] while other studies use AC motors, e.g., induction motors in [18][19][20], permanent magnet synchronous motors in [18,[21][22][23] and, finally, switched reluctance motors in [17,23,24]. The procedure used in each type of motor is similar, although the set of equations and control techniques could be different between them.
In addition, the electric motor needs a motor drive to control the desired vehicle speed. This motor drive must provide to the entire system good stability and dynamic response, which is achieved with the outer speed loop compensator and the inner current loop compensator, as shown in Figure 2. aerodynamic profile of the vehicle, among others. Besides, the terrain conditions have a considerable influence, e.g., the slope of the terrain and the friction between wheels and road. Another factor could be the environmental conditions such as wind or wet terrain. Different approaches can be used to model the vehicle. Reviewing the literature, some authors make use of several specific software. ADVISOR® (2002, Nat. Renewable Energy Lab, Golden, CO, USA) is used in [3] to model the vehicle with the aim of developing a power management strategy for an urban driving profile. MATLAB® (2017a, MathWorks, Natick, MA, USA) toolbox PSAT (2.1.11, University College Dublin, Dublin, Ireland) is used in [4] to model Zero-Emission Buses, and SIMULINK® (2019, MathWorks, Natick, MA, USA) is used in [5] to simulate a parallel hybrid powertrain and analyze the fuel consumption.
Another approach is to analyze the forces affecting the vehicle [6][7][8][9][10][11][12][13][14][15]. This approach is used in [6] to develop a sizing method for Hybrid Electric Vehicles (HEV), in [7] to analyze the sensitivity of the performance of a Battery Electric Vehicle (BEV), and in [8,9] to develop the energy consumption model of a plug-in HEV. In addition, this modeling approach is used in [10] to develop a technical assessment of utilizing an electrical variable transmission system in HEV, in [11] this approach is used to model and simulate an electric propulsion system of an EV, and in [12] is used to validate a HEV dynamic model for energy management including tire pressure distribution. Finally, in [13], this modeling approach is used to analyze the drivability of plug-in HEV, including the gear transmissions effects and, in [14,15] is used to design a control strategy for HEV with continuously variable transmission. These models are known as backward-facing vehicle models, which assume that the required demand is always met. The electric motor model is another key element. State of the art shows that some studies use DC motors [16,17] while other studies use AC motors, e.g., induction motors in [18][19][20], permanent magnet synchronous motors in [18,[21][22][23] and, finally, switched reluctance motors in [17,23,24]. The procedure used in each type of motor is similar, although the set of equations and control techniques could be different between them.
In addition, the electric motor needs a motor drive to control the desired vehicle speed. This motor drive must provide to the entire system good stability and dynamic response, which is achieved with the outer speed loop compensator and the inner current loop compensator, as shown in Figure 2.   On the other hand, the dynamic of each power converter included in the power distribution system is affected by the other power converters. Thus, the input impedance of the filter-motor drive-electric motor-vehicle combination must be considered to ensure system stability. As shown in Figure 1, the primary energy source, e.g., a fuel cell, is connected to the DC bus via a DC-DC power converter. The secondary energy source, the battery, is also connected to the DC bus via a bidirectional DC-DC power converter. In this type of architecture, stability cannot be guaranteed by default, unlike those architectures with the battery connected directly to the bus. For this reason, the impedance of the loads connected to the DC bus, in this case, the vehicle, motor, motor-drive, and filter, including the speed and current compensator, must be calculated, as shown in Figure 2.
Therefore, one of the main goals of the paper is to develop step-by-step a complete electric model of a light hybrid electric vehicle propulsion system, considering the motor, motor drive, and the light vehicle (including the vehicle mass, the radius and the mass of the wheels, the aerodynamic profile of the vehicle, among others). The model parameter values affect the relationship between the power demanded to the sources and the vehicle speed.
The second goal is to obtain the analytical expression of the small-signal model of the entire system composed of the electric motor, the motor drive (with two control loops), and the vehicle. This analytical expression shows the relationship among the wheel speed, the wheel torque, the electric motor torque, the electric motor current, and the current demanded to the DC bus, in closed-loop.
This model allows developing the electrical simulation and the optimal design of the vehicle power distribution system, as well as to get a controlled speed of the vehicle and a stable system.
Although the modeling procedure can be extended to other types of vehicles, this paper focuses on light vehicles, those whose weight is less than 1000 kg and which are mainly used in cites or in industrial speed-limited areas.
The third main goal is to obtain the input impedance model of the set composed of the input filter, the motor drive, the motor, and the vehicle, which is needed to assure the stability of the entire system. This analytical model results are compared to those obtained from the switched circuit and the averaged circuit simulated with PSIM ® (2015, Powersim, Rockville, MD, USA).
As this paper is focused on the motor and the vehicle modeling along with their integration in the entire system, the modeling of other parts of the overall system, e.g., the thermal management, battery management system, fuel cell, etc., are not considered.
The paper is organized as follows. First, Section 2 describes the vehicle and the electric motor model. With this information, the control loop to manage the speed vehicle is developed in Section 3. Next, the electrical input impedance analysis is shown in Section 4. Finally, Section 5 summarizes the main obtained conclusions.

Time Domain Modeling of the Vehicle and the Electric Motor
This section presents both the vehicle model and the motor model that allow the power profile to be calculated from a specific driving profile. The vehicle model is presented first, then the motor model is analyzed, and finally, the electric equivalent of both the electric motor and the vehicle is presented.
The vehicle model is based on the analysis of the applied forces on the vehicle and its impact on the electric motor operating conditions. These forces produce a torque on the shaft of the wheel that becomes a torque in the motor shaft.
In order to move the vehicle, the motor has to overcome a resistance torque, T R (t) (N·m), which is composed of two main parts, the torque needed to move the vehicle (acceleration torque) T Accel (t), and the torque to overcome the friction and the aerodynamic drag forces (rolling torque) T Roll (t) (1). The acceleration torque is composed of the torque needed to move the wheels T wm (t) and the torque needed to move the vehicle T vm (t).
T Accel (t) = T wm (t) + T vm (t) (2) The torque to move the wheels is defined in (3) where J wm is the moment of inertia of the wheel (kg·m 2 ), v is the speed of the vehicle (m/s), and, m w and r are the mass (kg) and radius (m) of the wheels.
The torque that moves the vehicle is defined in (4) where m is the mass of the vehicle (kg), and ω is the angular speed of the wheel.
A mass concentrated in the wheel surface is supposed in (3) and (4), i.e., the torque needed in the shaft to move the vehicle mass, m, or the wheel mass, m w , is equivalent to the torque needed to move a point particle of mass m or m w at a distance r of the shaft.
The rolling torque T Roll (t) is due to the friction of the wheels T fr (t) and the aerodynamic resistance T aer (t), both defined in (7) and (10), respectively.
T f r (t) = m·g· f r· cos β·r (7) In (7), g is the gravitational constant (m/s 2 ), β is the slope of the terrain ( • ), and fr is the friction coefficient defined in (8), where C sta and C dyn are the static and dynamic friction coefficient.
Then, replacing (9) and (10) in (6), the rolling torque results: T Roll (t) = [C sta ·m·g· cos β·r] + C dyn ·m·g· cos β·r 2 ·ω(t) + 1 2 ·ρ·C x ·S·r 3 ·ω 2 (t) (11) Using the angular speed instead of the linear speed and replacing (5) and (11) in (1), (12) is obtained: dt + [C sta ·m·g· cos β·r] + C dyn ·m·g· cos β·r 2 ·ω(t) Energies 2019, 12, 4082 5 of 20 Equation (12) has four terms. The first term is the acceleration torque (needed to move the mass of the vehicle and wheels), the second and third ones are the static and dynamic friction torque, respectively, and the last one is the aerodynamic torque. The forces that cause these torques are represented in Figure 3. Next, the electric motor model is developed. This motor transforms the electric power into mechanic power in the motor shaft. The DC motor is chosen to explain the procedure step-by-step, although other types of electric motors could be chosen using the same procedure. Figure 4 shows the model of the electric motor, being ω(t) the angular speed of the wheels (rad/s), TM(t) the motor torque (N·m) on the wheels, iarm(t) the current that goes through the rotor (A) and varm(t) the voltage between the motor connections (V). Rarm, Larm, and E(t) are, respectively, the electrical resistance (Ω), inductance (H), and counter-electromotive force (V) of the electric motor. The equations of a DC motor are (13) and (14) where K1 and K2 are constant parameters of the motor and, ωmot(t) and Tmot(t) are the angular speed and torque of the DC motor. The electric flux ϕ is assumed constant, which lets replacing K1 and K2, for Kω and KT, which are the speed constant (V/(rad/s)) and the torque constant (N·m/A).
The ratio between the motor shaft and the wheel shaft, rt, allows obtaining the motor torque, TM(t), and the angular speed of the wheel, ω(t). Replacing (15) and (16) into (13) and (14) derives in (17) and (18), which includes the efficiency of the vehicle transmission system and motor, ηT.
The Equation (19) describes the DC electric motor model. On the other hand, in steady-state, TR(t) has to be the same as the motor torque TM(t) that moves the vehicle, as can be seen in (20). Therefore, the Equations (17), (18), (19) and (20) are the four basic equation used to represent a DC electric motor. Next, the electric motor model is developed. This motor transforms the electric power into mechanic power in the motor shaft. The DC motor is chosen to explain the procedure step-by-step, although other types of electric motors could be chosen using the same procedure. Figure 4 shows the model of the electric motor, being ω(t) the angular speed of the wheels (rad/s), T M (t) the motor torque (N·m) on the wheels, i arm (t) the current that goes through the rotor (A) and v arm (t) the voltage between the motor connections (V). R arm , L arm , and E(t) are, respectively, the electrical resistance (Ω), inductance (H), and counter-electromotive force (V) of the electric motor. Next, the electric motor model is developed. This motor transforms the electric power into mechanic power in the motor shaft. The DC motor is chosen to explain the procedure step-by-step, although other types of electric motors could be chosen using the same procedure. Figure 4 shows the model of the electric motor, being ω(t) the angular speed of the wheels (rad/s), TM(t) the motor torque (N·m) on the wheels, iarm(t) the current that goes through the rotor (A) and varm(t) the voltage between the motor connections (V). Rarm, Larm, and E(t) are, respectively, the electrical resistance (Ω), inductance (H), and counter-electromotive force (V) of the electric motor.  (13) and (14) where K1 and K2 are constant parameters of the motor and, ωmot(t) and Tmot(t) are the angular speed and torque of the DC motor. The electric flux ϕ is assumed constant, which lets replacing K1 and K2, for Kω and KT, which are the speed constant (V/(rad/s)) and the torque constant (N·m/A).
The ratio between the motor shaft and the wheel shaft, rt, allows obtaining the motor torque, TM(t), and the angular speed of the wheel, ω(t). Replacing (15) and (16) into (13) and (14) derives in (17) and (18), which includes the efficiency of the vehicle transmission system and motor, ηT.
The Equation (19) describes the DC electric motor model. On the other hand, in steady-state, TR(t) has to be the same as the motor torque TM(t) that moves the vehicle, as can be seen in (20). Therefore, the Equations (17), (18), (19) and (20) are the four basic equation used to represent a DC electric motor.  (13) and (14) where K 1 and K 2 are constant parameters of the motor and, ω mot (t) and T mot (t) are the angular speed and torque of the DC motor. The electric flux φ is assumed constant, which lets replacing K 1 and K 2 , for K ω and K T , which are the speed constant (V/(rad/s)) and the torque constant (N·m/A).
The ratio between the motor shaft and the wheel shaft, rt, allows obtaining the motor torque, T M (t), and the angular speed of the wheel, ω(t). Replacing (15) and (16) into (13) and (14) derives in (17) and (18), which includes the efficiency of the vehicle transmission system and motor, η T .
The Equation (19) describes the DC electric motor model. On the other hand, in steady-state, T R (t) has to be the same as the motor torque T M (t) that moves the vehicle, as can be seen in (20). Therefore, the Equations (17), (18), (19) and (20) are the four basic equation used to represent a DC electric motor.
Summarizing, the electrical and mechanical equations of the electric motor are (17), (18), (19) and (12). In the model developed, two considerations are assumed, flat terrain and the linear friction coefficient.
The addends of (12) are shortened in (21). Then, the equivalent circuit of the electric motor is shown in Figure 5. where: T R1 = C sta ·m·g· cos β·r (23) Equation (21) can be interpreted as an electric circuit, where T M (t) can be interpreted as a current, as shown in Figure 5. Summarizing, the electrical and mechanical equations of the electric motor are (17), (18), (19) and (12). In the model developed, two considerations are assumed, flat terrain and the linear friction coefficient.
The addends of (12) are shortened in (21). Then, the equivalent circuit of the electric motor is shown in Figure 5. where: Equation (21) can be interpreted as an electric circuit, where TM(t) can be interpreted as a current, as shown in Figure 5.

Control Loop Design
A control circuit is designed with a speed outer loop and a current inner loop to control the motor operation. The vehicle and motor small-signal model must be obtained to achieve this design.

Vehicle Small-Signal Model
The small-signal model of the vehicle and the motor is obtained by linearizing the terms included in (21) at a working point W, which is the angular speed of the wheels. Concerning the notation, henceforth, the terms â is a perturbation of a variable a. To obtain the linearization of the acceleration torque at W, (26) is replaced in (21): which gives (27): where (27) is composed of the motor torque at the working point W, , and the perturbation of the motor torque, ( ), defined in (28) and (29) respectively:

Control Loop Design
A control circuit is designed with a speed outer loop and a current inner loop to control the motor operation. The vehicle and motor small-signal model must be obtained to achieve this design.

Vehicle Small-Signal Model
The small-signal model of the vehicle and the motor is obtained by linearizing the terms included in (21) at a working point W, which is the angular speed of the wheels. Concerning the notation, henceforth, the terms â is a perturbation of a variable a. To obtain the linearization of the acceleration torque at W, (26) is replaced in (21): which gives (27): where (27) is composed of the motor torque at the working point W, T M W , and the perturbation of the motor torque,T M (t), defined in (28) and (29) respectively: In (28), dW/dt = 0, because W is not a function of time. In (29),ω 2 (t) is neglected because the terms of the second degree or higher are not considered. Therefore, the motor torque can be expressed, as shown in (30): Finally, applying the Laplace transformation to (29) gives (31): where J T and Rol are defined in (22) and (32), considering (24) and (25):

Vehicle and Motor Control Loop
The control loop together with the blocks diagram explained in Figure 2 and Figure 6 is shown in Figure 7. Ri(s) is the transfer function of the current compensator, Mod(s) includes both the modulator and the motor drive power stage, and Rv(s) is the transfer function of the speed compensator. Finally, the terms Ksc(s) and Ksω(s) are the gain related to the internal current loop sensor and external speed loop sensor, respectively. The reference signal is the driving profile, ωref(s).

Vehicle and Motor Control Loop
The control loop together with the blocks diagram explained in Figures 2 and 6 is shown in Figure 7. R i (s) is the transfer function of the current compensator, Mod(s) includes both the modulator and the motor drive power stage, and R v (s) is the transfer function of the speed compensator. Finally, the terms K sc (s) and K sω (s) are the gain related to the internal current loop sensor and external speed loop sensor, respectively. The reference signal is the driving profile, ω ref (s). The next step is to simplify the blocks diagram shown in Figure 7. First, the inner loops must be grouped. To achieve that, the blocks of the right side are shifted as shown in Figure 8.
Finally, the two blocks displaced before, now go back to the initial position, as shown in Figure  10.  The next step is to simplify the blocks diagram shown in Figure 7. First, the inner loops must be grouped. To achieve that, the blocks of the right side are shifted as shown in Figure 8. The next step is to simplify the blocks diagram shown in Figure 7. First, the inner loops must be grouped. To achieve that, the blocks of the right side are shifted as shown in Figure 8.
Finally, the two blocks displaced before, now go back to the initial position, as shown in Figure  10.   The next step is to simplify the blocks diagram shown in Figure 7. First, the inner loops must be grouped. To achieve that, the blocks of the right side are shifted as shown in Figure 8.
Finally, the two blocks displaced before, now go back to the initial position, as shown in Figure  10. Finally, the two blocks displaced before, now go back to the initial position, as shown in Figure 10.
Finally, the two blocks displaced before, now go back to the initial position, as shown in Figure  10.  T Next, the inner loop is simplified in Figure 11, where BCi(s) is the closed-loop transfer function of the current inner loop (41).
Next, the inner loop is simplified in Figure 11, where BCi(s) is the closed-loop transfer function of the current inner loop (41).
From the block diagram shown in Figure 11, the transfer function of the speed open loop Tω(s) can be calculated to design the speed compensator Rv(s).
Based on the Equations (40) and (42), the inner loop compensator and the outer loop compensator can be obtained.
As an example, the data to calculate both compensators, which are obtained using the design software SmartCtrl® (version, manufacturer, city, state abbreviation if USA and Canada, Country) [25], are shown in Table 1. In this case, the gain of the current loop Ksc(s) and speed loop Ksω(s) are assumed constant. The vehicle and motor data are summarized in Table 2 and Table 3, shown in Section 4. The driving profile shown in Figure 12 is used to evaluate the designed control. It corresponds to a European driving profile, ECE-15, applied to light vehicles. From the block diagram shown in Figure 11, the transfer function of the speed open loop T ω (s) can be calculated to design the speed compensator R v (s).
Based on the Equations (40) and (42), the inner loop compensator and the outer loop compensator can be obtained.
As an example, the data to calculate both compensators, which are obtained using the design software SmartCtrl ® (version, manufacturer, city, state abbreviation if USA and Canada, Country) [25], are shown in Table 1. In this case, the gain of the current loop K sc (s) and speed loop K sω (s) are assumed constant. The vehicle and motor data are summarized in Tables 2 and 3, shown in Section 4.  [26].

Variable Value
Wheel radius, r 0.27 m Wheel total mass, m w 60 kg Vehicle mass, m 800 kg Aerodynamic coefficient, C dyn 0.31 Frontal surface of the vehicle, S 1.7 m 2 Driver and filter efficiency, η DF 99% Vehicle transmission system and motor efficiency, η T 75% The driving profile shown in Figure 12 is used to evaluate the designed control. It corresponds to a European driving profile, ECE-15, applied to light vehicles.
From the block diagram shown in Figure 11, the transfer function of the speed open loop Tω(s) can be calculated to design the speed compensator Rv(s).
Based on the Equations (40) and (42), the inner loop compensator and the outer loop compensator can be obtained.
As an example, the data to calculate both compensators, which are obtained using the design software SmartCtrl® (version, manufacturer, city, state abbreviation if USA and Canada, Country) [25], are shown in Table 1. In this case, the gain of the current loop Ksc(s) and speed loop Ksω(s) are assumed constant. The vehicle and motor data are summarized in Table 2 and Table 3, shown in Section 4. The driving profile shown in Figure 12 is used to evaluate the designed control. It corresponds to a European driving profile, ECE-15, applied to light vehicles.   Figure 13 shows the results of applying the ECE-15 driving profile, shown in Figure 12, to the motor-vehicle system. Figure 13a shows the wheel angular speed that must be achieved, ω ref , and the simulated, ω. Figure 13b shows both the torque in the wheel, T M , and the torque in the motor, T mot . The difference between them is due to the ratio between the motor shaft and the wheel shaft, rt. Figure 13c shows the current through the rotor, i arm , and the demanded current to the DC bus. As the motor input voltage is lower than the DC bus voltage, the current through the rotor is higher than the current demanded to the DC bus.
As can be seen in Figure 13, the speed increments or decrements imply positive or negative motor torque, which involve positive or negative currents in the rotor. Thus, using the carried out modeling a driving profile can be transformed into a current profile that has to be supplied by the power distribution system from the DC bus. motor-vehicle system. Figure 13 a shows the wheel angular speed that must be achieved, ωref, and the simulated, ω. Figure 13 b shows both the torque in the wheel, TM, and the torque in the motor, Tmot. The difference between them is due to the ratio between the motor shaft and the wheel shaft, rt. Figure 13 c shows the current through the rotor, iarm, and the demanded current to the DC bus. As the motor input voltage is lower than the DC bus voltage, the current through the rotor is higher than the current demanded to the DC bus.  Table 1 to Table 3.
As can be seen in Figure 13, the speed increments or decrements imply positive or negative motor torque, which involve positive or negative currents in the rotor. Thus, using the carried out modeling a driving profile can be transformed into a current profile that has to be supplied by the power distribution system from the DC bus.

Electrical Impedance Analysis
In this section, the small-signal input impedance of the system composed of the vehicle, the electric motor, the motor drive, and the filter is obtained. This input impedance can affect system stability. Figure 14 shows the analyzed system, divided into three main blocks. The first one is the motor drive that enables the vehicle to run and brake, which includes an inner loop to control iarm and outer loop to control ω. The second one is the input circuit of the electric motor (shown as an Electric motor). Finally, the last one is the vehicle equivalent circuit. In addition, Figure 14 shows the input filter of the motor drive-electric motor-vehicle system.
The outer loop compares the angular speed of the motor, ω(t), with the reference signal, which depends on the driving profile. The inner loop uses the output of the outer loop as a reference signal and generates the modulator input signal.
In this analysis, it is assumed that the motor speed dynamics is much slower than the motor current dynamics. This hypothesis allows us to assume that the counter-electromotive force E(t) is constant. Therefore, the reference signal of the current loop, which is the output of the speed loop, Figure 13. Example of the main Motor-vehicle waveforms obtained with an ECE-15 driving profile and the data included in Table 1 to Table 3.

Electrical Impedance Analysis
In this section, the small-signal input impedance of the system composed of the vehicle, the electric motor, the motor drive, and the filter is obtained. This input impedance can affect system stability. Figure 14 shows the analyzed system, divided into three main blocks. The first one is the motor drive that enables the vehicle to run and brake, which includes an inner loop to control i arm and outer loop to control ω. The second one is the input circuit of the electric motor (shown as an Electric motor). Finally, the last one is the vehicle equivalent circuit. In addition, Figure 14 shows the input filter of the motor drive-electric motor-vehicle system.  The averaged model of the motor drive shown in Figure 15 is calculated to obtain the electric motor impedance, where vcontrol(t) is the output signal of the inner loop compensator, and vtri(t) is the modulator signal. The outer loop compares the angular speed of the motor, ω(t), with the reference signal, which depends on the driving profile. The inner loop uses the output of the outer loop as a reference signal and generates the modulator input signal.
In this analysis, it is assumed that the motor speed dynamics is much slower than the motor current dynamics. This hypothesis allows us to assume that the counter-electromotive force E(t) is constant. Therefore, the reference signal of the current loop, which is the output of the speed loop, can also be assumed constant. Taking this into account, there is no dependence between the electric motor circuit and the vehicle equivalent circuit from the point of view of the current loop operation.
The averaged model of the motor drive shown in Figure 15 is calculated to obtain the electric motor impedance, where v control (t) is the output signal of the inner loop compensator, and v tri (t) is the modulator signal. The averaged model of the motor drive shown in Figure 15 is calculated to obtain the electric motor impedance, where vcontrol(t) is the output signal of the inner loop compensator, and vtri(t) is the modulator signal. The relation between the bus power, PBus(t), and the motor power, Pmot(t), is given by (43), where ηDF is the efficiency of both motor drive and filter (not included Figure 5): where Pmot(t) and PBus(t) are: The averaged output voltage vAB and the averaged input current ig are: where d is the duty cycle, defined as: The relation between the bus power, P Bus (t), and the motor power, P mot (t), is given by (43), where η DF is the efficiency of both motor drive and filter (not included Figure 5): where P mot (t) and P Bus (t) are: The averaged output voltage v AB and the averaged input current i g are: where d is the duty cycle, defined as: Afterward, once the motor drive averaged model is known, the motor drive-electric motor-vehicle system averaged model can be obtained. The model, shown in Figure 16, uses the previously claimed assumptions.

= (48)
Afterward, once the motor drive averaged model is known, the motor drive-electric motorvehicle system averaged model can be obtained. The model, shown in Figure 16, uses the previously claimed assumptions. The obtained model is nonlinear. Therefore, the Equations (46) and (47) must be linearized at the working points Vcontrol, VBus, and Iarm. Figure 17 depicts the small-signal model of the motor drive-electric motor-vehicle system without the control loop. In Figure 17, Z'motor(s) is the motor drive-electric motor-vehicle system input impedance. The block diagram of the motor drive-electric motor-vehicle system small-signal model, together with the control loop, is shown in Figure 18. The equations used to obtain the block diagram are (36), (49), (50), where the counter-electromotive force E(t) is surmised constant. The obtained model is nonlinear. Therefore, the Equations (46) and (47) must be linearized at the working points V control , V Bus , and I arm .î Figure 17 depicts the small-signal model of the motor drive-electric motor-vehicle system without the control loop.

= (48)
Afterward, once the motor drive averaged model is known, the motor drive-electric motorvehicle system averaged model can be obtained. The model, shown in Figure 16, uses the previously claimed assumptions. The obtained model is nonlinear. Therefore, the Equations (46) and (47) must be linearized at the working points Vcontrol, VBus, and Iarm. Figure 17 depicts the small-signal model of the motor drive-electric motor-vehicle system without the control loop. In Figure 17, Z'motor(s) is the motor drive-electric motor-vehicle system input impedance. The block diagram of the motor drive-electric motor-vehicle system small-signal model, together with the control loop, is shown in Figure 18. The equations used to obtain the block diagram are (36), (49), (50), where the counter-electromotive force E(t) is surmised constant. In Figure 17, Z' motor (s) is the motor drive-electric motor-vehicle system input impedance. The block diagram of the motor drive-electric motor-vehicle system small-signal model, together with the control loop, is shown in Figure 18 As stated above, the reference of the current is assumed constant. This way, a simplified model is shown in Figure 19. As stated above, the reference of the current is assumed constant. This way, a simplified model is shown in Figure 19. As stated above, the reference of the current is assumed constant. This way, a simplified model is shown in Figure 19. As the goal is to obtain the transfer function that relates the current that goes through the rotor iarm, the duty cycle d, and bus voltage vBus, the blocks diagram shown in Figure 20 can be considered. The two diagrams placed on the left side in Figure 20 are named Giv(s) and Gdv(s). Their transfer functions appear in (51) and (52). They are obtained from the blocks diagram in Figure 19. where: Replacing (53) in (52) gives: Figure 19. Simplified block diagram of the motor drive-electric motor-vehicle system with the inner control loop and without current reference signal variation.
As the goal is to obtain the transfer function that relates the current that goes through the rotor i arm , the duty cycle d, and bus voltage v Bus , the blocks diagram shown in Figure 20 can be considered. Figure 18. Block diagram of the motor drive-electric motor-vehicle system, including the inner control loop.
As stated above, the reference of the current is assumed constant. This way, a simplified model is shown in Figure 19. As the goal is to obtain the transfer function that relates the current that goes through the rotor iarm, the duty cycle d, and bus voltage vBus, the blocks diagram shown in Figure 20 can be considered. The two diagrams placed on the left side in Figure 20 are named Giv(s) and Gdv(s). Their transfer functions appear in (51) and (52). They are obtained from the blocks diagram in Figure 19. where: Replacing (53) in (52) gives: The two diagrams placed on the left side in Figure 20 are named Giv(s) and Gdv(s). Their transfer functions appear in (51) and (52). They are obtained from the blocks diagram in Figure 19. where: Replacing (53) in (52) gives: Once Giv(s) and Gdv(s) are known, the input impedance of the motor drive-electric motor-vehicle system represented in Figure 17, Z' motor (s), can be obtained using the block diagram shown Figure 20.
Last, to obtain the input impedance, the working points of D, and I arm must be calculated. To achieve it, the steady-state circuit shown in Figure 21 is used together with the electric motor mechanical equations.
Last, to obtain the input impedance, the working points of D, and Iarm must be calculated. To achieve it, the steady-state circuit shown in Figure 21 is used together with the electric motor mechanical equations. On the other hand, using the motor torque Equation (59), Iarm is derived in (60).
Finally, D can be obtained, as shown in (61).
The vehicle and electric motor data must be known to develop the impedance analysis. The data used are shown in Table 2 and Table 3, which correspond to a light vehicle and an electric motor model GEM 5BC449JB6007.  Equations (56) and (57) are obtained from the static circuit of Figure 21.
E is defined in (58), where W is the angular speed at the working point: On the other hand, using the motor torque Equation (59), I arm is derived in (60).
Finally, D can be obtained, as shown in (61).
The vehicle and electric motor data must be known to develop the impedance analysis. The data used are shown in Tables 2 and 3, which correspond to a light vehicle and an electric motor model GEM 5BC449JB6007.
If the electric motor allows energy recovery, the motor has two modes. In motor mode, the vehicle is accelerating, and thus, the motor demands energy from the system. In generator mode, the vehicle is decelerating, and therefore, the motor supplies energy to the system. In both modes, the model and behavior of the motor drive, the electric motor, and the control loops are the same. Thus the impedance calculation procedure remains equal.
The working point W is chosen using a standard driving profile. This profile depends on the vehicle and the driving pattern used. As commented before, the driving profile chosen is the Urban Driving Cycle (UDC or ECE-15), which represents the typical driving conditions of busy European cities, with 50 km/h maximum speed. The working points in both motor and generator mode are shown in Figure 22, and the corresponding data is shown in Table 4. the impedance calculation procedure remains equal. The working point W is chosen using a standard driving profile. This profile depends on the vehicle and the driving pattern used. As commented before, the driving profile chosen is the Urban Driving Cycle (UDC or ECE-15), which represents the typical driving conditions of busy European cities, with 50 km/h maximum speed. The working points in both motor and generator mode are shown in Figure 22, and the corresponding data is shown in Table 4.  Next, the frequency response of the small-signal input impedance obtained by three different methods are shown in Figure 23 a and b for motor mode, and in Figure 23 c and d for generator mode. The first method is the frequency response obtained with the averaged circuit, shown in Figure 16, by using Equations (46) and (47). The second method is the frequency response obtained from the switched circuit, shown in Figure 14. The third method is the frequency response obtained from the small-signal analytical model, shown in Figure 17, by using Equations (49) and (50). Next, the frequency response of the small-signal input impedance obtained by three different methods are shown in Figure 23a,b for motor mode, and in Figure 23c,d for generator mode. The first method is the frequency response obtained with the averaged circuit, shown in Figure 16, by using Equations (46) and (47). The second method is the frequency response obtained from the switched circuit, shown in Figure 14. The third method is the frequency response obtained from the small-signal analytical model, shown in Figure 17, by using Equations (49) and (50). In both motor mode and generator mode, the input impedance of the converter presents two parts. Below 1 kHz, the impedance can be represented as a resistance, and above this frequency the input impedance has inductive behavior. Note that the Z´motor(s) does not include the typical input capacitor, which in this case belongs to the filter block. The only difference between motor and generator mode is the phase value, so the motor mode has negative impedance, and the generator mode has positive impedance.
The mean absolute percentage error (MAPE) is measured using the natural impedance values (Ω) instead of dB. On the one hand, by comparing the switched circuit frequency response and the averaged circuit frequency response the MAPE reaches a value of 2.17%. On the other hand, by In both motor mode and generator mode, the input impedance of the converter presents two parts. Below 1 kHz, the impedance can be represented as a resistance, and above this frequency the input impedance has inductive behavior. Note that the Z' motor (s) does not include the typical input capacitor, which in this case belongs to the filter block. The only difference between motor and generator mode is the phase value, so the motor mode has negative impedance, and the generator mode has positive impedance.
The mean absolute percentage error (MAPE) is measured using the natural impedance values (Ω) instead of dB. On the one hand, by comparing the switched circuit frequency response and the averaged circuit frequency response the MAPE reaches a value of 2.17%. On the other hand, by comparing the small-signal analytical model frequency response and the averaged circuit frequency response, the higher MAPE value is 3.98%.
However, as the transmission of high-frequency ripples to the system must be limited, a filter is placed between the motor drive and the DC bus, as Figure 14 illustrates. Therefore, the impedance including the filter, Z motor (s), shown in Figure 24, must be obtained.
electric motor-vehicle system in both motor mode, (a) gain and (b) phase; and generator mode, (c) gain and (d) phase.
In both motor mode and generator mode, the input impedance of the converter presents two parts. Below 1 kHz, the impedance can be represented as a resistance, and above this frequency the input impedance has inductive behavior. Note that the Z´motor(s) does not include the typical input capacitor, which in this case belongs to the filter block. The only difference between motor and generator mode is the phase value, so the motor mode has negative impedance, and the generator mode has positive impedance.
The mean absolute percentage error (MAPE) is measured using the natural impedance values (Ω) instead of dB. On the one hand, by comparing the switched circuit frequency response and the averaged circuit frequency response the MAPE reaches a value of 2.17%. On the other hand, by comparing the small-signal analytical model frequency response and the averaged circuit frequency response, the higher MAPE value is 3.98%.
However, as the transmission of high-frequency ripples to the system must be limited, a filter is placed between the motor drive and the DC bus, as Figure 14 illustrates. Therefore, the impedance including the filter, Zmotor(s), shown in Figure 24, must be obtained.
The goal of the filter is to reduce the propagation of the high-frequency switching ripple of the motor drive-electric motor-vehicle system towards the rest of the powertrain. The data related to the filter are shown in Table 5. Table 5. Input filter data of the motor drive-electric motor-vehicle system.

Variable
Value The frequency responses of the small-signal input impedance, Z motor (s), are shown in Figure 25a,b for motor mode, and in Figure 25c,d for generator mode. Figure 25 shows the averaged circuit, shown in Figure 16, using Equations (46) and (47); the switched circuit, shown in Figure 14; and the small-signal analytical model, shown in Figure 17, using Equations (49) and (50).
In both modes, a resonance and an anti-resonance can be observed. The resonance is imposed by L fm and C 2m , 1.59 kHz, and the anti-resonance depends on the L fm and C 1m , 35.58 kHz.
At low frequency, the filter does not affect, but above 159 Hz the input impedance is determined by the input filter. circuit, shown in Figure 16, using Equations (46) and (47); the switched circuit, shown in Figure 14; and the small-signal analytical model, shown in Figure 17, using Equations (49) and (50).
In both modes, a resonance and an anti-resonance can be observed. The resonance is imposed by Lfm and C2m, 1.59 kHz, and the anti-resonance depends on the Lfm and C1m, 35.58 kHz.
At low frequency, the filter does not affect, but above 159 Hz the input impedance is determined by the input filter.  The higher MAPE obtained by comparing the switched circuit frequency response and the averaged circuit frequency response is 2.43%. On the other hand, the MAPE reaches a maximum value of 4.6% when comparing the small-signal analytical model frequency response and the averaged circuit frequency response.
Finally, as can be noticed in Figures 23 and 25, the obtained analytical model results fit with the simulation results obtained by PSIM ® [27], and therefore, the analytical models are validated.

Conclusions
This paper develops, step-by-step and in detail, the main analyses that must be addressed to obtain a complete electric model of a light hybrid electric vehicle propulsion system, as well as the input impedance expressions necessary to evaluate the system stability.
The complete electric model of a light hybrid electric vehicle propulsion system includes the most important parameters that affect the relationship between the vehicle speed and the energy demanded from the power sources. These parameters are the vehicle mass, the radius and the mass of the wheels, the aerodynamic profile of the vehicle, the electric motor, the motor drive, and the motor and transmission system efficiency, among others.
In addition, the paper shows the procedure to obtain the analytical expression of the small-signal model of the complete analyzed system compound by the vehicle, the electric motor, the motor drive (including two control loops), and the input filter. The inner and outer open-loop gains, which determine the current and speed compensator included in the motor drive, has been obtained from the small-signal model. These loops provide good stability and dynamic response, which helps the motor to achieve the desired speed. With this information, the simulation results show the relationship among speed vehicle, the wheel torque, the electric motor torque, the current demanded by the DC bus, and the electric motor current, in closed-loop.