Intrinsic Characteristics of Forward Simulation Modeling Electric Vehicle for Energy Analysis

: The forward method for modeling electric vehicles is one of the most suitable for estimating energy consumption in different imposed driving cycles. However, a detailed description of the methodology used for the development of electric vehicle models is necessary and is scarce in the current literature. To ﬁll this gap, this study focuses on highlighting the intrinsic characteristics through a theoretical study with a mathematical model, complemented by demonstrative simulations in Matlab/Simulink. The results show that the forward method can be estimated more accurately based on the energy consumption of the electric vehicle. Moreover, this paper aims to be explicitly descriptive for the development of more complex electric vehicle models to incorporate real driving cycles, being able to size the drivetrain of the vehicle itself or develop ecological routes.


Introduction
Numerical simulation based on mathematical modeling is a general method to evaluate the energy consumption of an electric vehicle (EV) on a certain route. The numerical methods are widely applicable, from statistical models to evaluate the CO 2 emissions, kinematic models for the simulation of microscopic traces, or complex dynamic models to determine the influence of components [1]. There is a variety of computational tools that allow validation of the operation of EVs, focusing on responses to driving cycles [2]. The simulation process can be "backward" or "forward", or even a combination between them. The former is relatively easier to use/implement, while the latter is more complex when adjusting to real driving cycles. Although there have been numerous studies that address the "forward" method for the development of EV models, these studies do not tackle a clear and concise description of the components involved in the development of models and estimate their energy consumption based on the imposed driving cycle. The study of the intrinsic characteristics of the forward method creates a starting point to develop computational models based on artificial neural networks (ANN), with the aim of determining the relationships between a series of factors that affect the energy consumption of EVs as a function of input parameters.
There is a wide variety of EV models available in the literature (e.g., see [3][4][5][6][7][8][9] and references therein), which, depending on the required detail of each component, can be stationary, quasi-stationary, or dynamic. Ref. [3] presents the modeling and simulation of a hybrid EV, while [4] focuses on the acoustic simulation of engines. In [5], the authors present a flat dynamic simulation of vehicles based on a real driving cycle. In [6,7], the modeling and estimation of the energy consumption of an EV is carried out, with rechargeable storage systems [6] and accurate EV energy model estimation [7]. Likewise, in [8,9], the authors focus on the modeling of high-performance powertrains, considering simulation and validation [8], as well as transmission parameters [9]. Most EV models

•
Highlighting the intrinsic characteristics of the forward method; • Providing a systematic description of the blocks involved, together with their equations and necessary considerations for the development of the model; • Modeling, simulation, and validation of an electric vehicle by the forward method; • Energy analysis of the electric vehicle before an Urban Dynamometer Driving Schedule (UDDS) driving cycle.
The remainder of the paper is organized as follows: Section 2 presents the proposed materials and methods, including a mathematical representation of components using the forward method. Section 3 presents a case study with the results of the study. Finally, the main conclusions are duly drawn in Section 4.

Materials and Methods
The methodology used in this paper is schematically depicted in Figure 1. Firstly, the necessary background is presented, including a brief introduction to the forward method, which allows us to clearly identify its intrinsic characteristics. Next, this method is used to perform extensive simulations on a benchmark EV model, including the mathematical notations and model description in Matlab-Simulink, in order to demonstrate the consistency of the mathematical modeling.

Background
The backward simulation calculates the power demand of various components according to the driving cycle, to be posteriorly distributed among power sources in order to meet the vehicle's demand. Finally, the energy consumption of the motor and batteries is obtained [16]. In contrast, in the forward simulation, the model determines the energy demand in the process of the driving cycle, to transmit the energy demand to the vehicle controller. This controller receives the signal and considers the current state of the vehicle components to optimize and determine power allocation, together with the output of each power source. The movement flow is transferred from the driving system, passing through the transmission to the tire module. Then, by means of equations that model the driving module, the real speed that is transmitted to the driver is calculated to form a closed cycle. Figure 2 represents a simple scheme of the simulation process by the "forward" method. The driving cycle provides the speed to pass through the driver model, which controls the longitudinal interfaces of the vehicle (accelerator and brake pedal), based on the difference between the object and the vehicle speed. The battery supplies the motor with energy and the motor provides a torque that, through the kinematic chain, transfers the movement to the wheels, where the traction maintains the propulsion force. The position returns to the driving cycle to find a target speed, closing the calculation loop. Electricity 2022, 3, FOR PEER REVIEW 3 notations and model description in Matlab-Simulink, in order to demonstrate the consistency of the mathematical modeling.

Background
The backward simulation calculates the power demand of various components according to the driving cycle, to be posteriorly distributed among power sources in order to meet the vehicle's demand. Finally, the energy consumption of the motor and batteries is obtained [16]. In contrast, in the forward simulation, the model determines the energy demand in the process of the driving cycle, to transmit the energy demand to the vehicle controller. This controller receives the signal and considers the current state of the vehicle components to optimize and determine power allocation, together with the output of each power source. The movement flow is transferred from the driving system, passing through the transmission to the tire module. Then, by means of equations that model the driving module, the real speed that is transmitted to the driver is calculated to form a closed cycle. Figure 2 represents a simple scheme of the simulation process by the "forward" method. The driving cycle provides the speed to pass through the driver model, which controls the longitudinal interfaces of the vehicle (accelerator and brake pedal), based on the difference between the object and the vehicle speed. The battery supplies the motor with energy and the motor provides a torque that, through the kinematic chain, transfers the movement to the wheels, where the traction maintains the propulsion force. The position returns to the driving cycle to find a target speed, closing the calculation loop.   Figure 2. Schematic representation of the forward simulation process.

Mathematical Modeling of the Forward Method
This section provides the mathematical foundations of the forward method, as well as the explanation of the blocks depicted in Figure 2.

Driving Cycle Model
To estimate the SOC of on-board batteries and EV driving range, simulation models are based on internationally legislated driving cycles, and some represetative examples are as follows: Federal Test Procedure 72/75 (FTP 72/75), the New European Driving Cycle (NEDC) (see Figure 3a), the Japanese Cycle 08 (JC08), the Worldwide Harmonized Light Duty Cycle (WLTP), the 15-mode test cycle of the Economic Commission for Europe (ECE 15), and the UDDS (see Figure 3b), among many others. All of them differ in the sections of acceleration, braking, speed, cycle duration time, etc. [22][23][24][25][26].

Mathematical Modeling of the Forward Method
This section provides the mathematical foundations of the forward method, as well as the explanation of the blocks depicted in Figure 2.

Driving Cycle Model
To estimate the SOC of on-board batteries and EV driving range, simulation models are based on internationally legislated driving cycles, and some represetative examples are as follows: Federal Test Procedure 72/75 (FTP 72/75), the New European Driving Cycle (NEDC) (see Figure 3a), the Japanese Cycle 08 (JC08), the Worldwide Harmonized Light Duty Cycle (WLTP), the 15-mode test cycle of the Economic Commission for Europe (ECE 15), and the UDDS (see Figure 3b), among many others. All of them differ in the sections of acceleration, braking, speed, cycle duration time, etc. [22,23,24,25,26]. In 2019, Europe adopted the WLTP as standard, thus rendering obsolete the NEDC. This cycle tests vehicles on longer routes, higher accelerations, and higher speeds; however, it is still believed that it does not fit reality. EVs are initially tested on a dynamometric bench in different driving cycles; one of the most important is the UDDS, which simulates driving in urban areas intermittently, starts the vehicle and accelerates it to the maximum speed allowed and stops it again. It is used to measure the energy consumption in the city. In the development of EV models for simulation, the UDDS cycle is conventionally used, by which the EV is subjected to multiple test cycles until the battery is completely discharged and stops by itself, thus providing preliminary autonomy. This study applies this driving cycle and not the obsolete NEDC cycle or the WLTP homologation cycle.

Driver Model
The controller module is suitable for the forward simulation model, by allowing the simulation model to form a closed-loop system, for which the PID controllers (Proportional, Integral, Derivative) are used, as observed in Figure 4. The forward method requires a driver model that acts as an accelerator/brake pedal actuator to follow the target speed (imposed driving cycle). For this, various controller models can be used; one of the most common is the PID controller, whose scheme is shown in Figure 4 [27]. The PID controller is based on repetitive tests and experience to determine the output values. The control system has a minimum static error between the actual simulation speed and the target speed; this error is necessary to maintain the torque of the drive wheels, because the changes in the simulation speed are always a function of the changes in the target speed [28]. If the PID parameters are properly calibrated, the delays will be slight and the In 2019, Europe adopted the WLTP as standard, thus rendering obsolete the NEDC. This cycle tests vehicles on longer routes, higher accelerations, and higher speeds; however, it is still believed that it does not fit reality. EVs are initially tested on a dynamometric bench in different driving cycles; one of the most important is the UDDS, which simulates driving in urban areas intermittently, starts the vehicle and accelerates it to the maximum speed allowed and stops it again. It is used to measure the energy consumption in the city. In the development of EV models for simulation, the UDDS cycle is conventionally used, by which the EV is subjected to multiple test cycles until the battery is completely discharged and stops by itself, thus providing preliminary autonomy. This study applies this driving cycle and not the obsolete NEDC cycle or the WLTP homologation cycle.

Driver Model
The controller module is suitable for the forward simulation model, by allowing the simulation model to form a closed-loop system, for which the PID controllers (Proportional, Integral, Derivative) are used, as observed in Figure 4. The forward method requires a driver model that acts as an accelerator/brake pedal actuator to follow the target speed (imposed driving cycle). For this, various controller models can be used; one of the most common is the PID controller, whose scheme is shown in Figure 4 [27]. The PID controller is based on repetitive tests and experience to determine the output values. The control system has a minimum static error between the actual simulation speed and the target speed; this error is necessary to maintain the torque of the drive wheels, because the changes in the simulation speed are always a function of the changes in the target speed [28]. If the PID parameters are properly calibrated, the delays will be slight and the motor drive torque is set efficiently; otherwise, the delays will be longer and could even cause a system discharge.  The transfer function of the PID controller is given by Equation (1).
where k p is the proportional coefficient, k i is the integral coefficient, and k d is the derivative coefficient.
In [29], the authors compared the conventional PID control with the linear square scheme (LQR), determining that LQR may outperform the conventional PID control. In particular, the scheme shown in Figure 5  With these settings, it is observed in Figure 6 that the LQR controller improves the overshoot and settling time compared to the PID controller that takes values of P = 0.00012, I = 0.00047, and D = 0.000068 [30]. However, the PID controller has a shorter response time, being better behaved with fast speed transitions, maintaining a minimum error between the target and real speeds. On the basis of this conclusion, it is enough to apply a PID controller that performs the function of the pedal (acceleration/brake) to follow the target speed profile.    The transfer function of the PID controller is given by Equation (1).
where k p is the proportional coefficient, k i is the integral coefficient, and k d is the derivative coefficient.
In [29], the authors compared the conventional PID control with the linear square scheme (LQR), determining that LQR may outperform the conventional PID control. In particular, the scheme shown in Figure 5 Figure 6 that the LQR controller improves the overshoot and settling time compared to the PID controller that takes values of P = 0.00012, I = 0.00047, and D = 0.000068 [30]. However, the PID controller has a shorter response time, being better behaved with fast speed transitions, maintaining a minimum error between the target and real speeds. On the basis of this conclusion, it is enough to apply a PID controller that performs the function of the pedal (acceleration/brake) to follow the target speed profile. The transfer function of the PID controller is given by Equation (1).
where k p is the proportional coefficient, k i is the integral coefficient, and k d is the derivative coefficient.
In [29], the authors compared the conventional PID control with the linear square scheme (LQR), determining that LQR may outperform the conventional PID control. In particular, the scheme shown in Figure 5  With these settings, it is observed in Figure 6 that the LQR controller improves the overshoot and settling time compared to the PID controller that takes values of P = 0.00012, I = 0.00047, and D = 0.000068 [30]. However, the PID controller has a shorter response time, being better behaved with fast speed transitions, maintaining a minimum error between the target and real speeds. On the basis of this conclusion, it is enough to apply a PID controller that performs the function of the pedal (acceleration/brake) to follow the target speed profile.     The transfer function of the PID controller is given by Equation (1).
where k p is the proportional coefficient, k i is the integral coefficient, and k d is the derivative coefficient.
In [29], the authors compared the conventional PID control with the linear square scheme (LQR), determining that LQR may outperform the conventional PID control. In particular, the scheme shown in Figure 5  With these settings, it is observed in Figure 6 that the LQR controller improves the overshoot and settling time compared to the PID controller that takes values of P = 0.00012, I = 0.00047, and D = 0.000068 [30]. However, the PID controller has a shorter response time, being better behaved with fast speed transitions, maintaining a minimum error between the target and real speeds. On the basis of this conclusion, it is enough to apply a PID controller that performs the function of the pedal (acceleration/brake) to follow the target speed profile.

Brake Model
The braking force is proportional to the torque generated at the wheels. In order to stop the tire, a force must be applied that is contrary to the rotation movement, which is calculated in (2). In fact, Equation (2) calculates the necessary braking for just one wheel, it being necessary to multiply this result by the total number of wheels in order to obtain the necessary total braking force in the vehicle [31].
where N f is the braking torque (N), r d is the radius of the tire (m). The ratio that determines the recoverable energy by means of regenerative braking is general and applies to any driving cycle [32]. The regenerative braking energy efficiency η rb is defined by Equation (3) [33,34].
where E recoverable [kWh ] is the energy recovered during braking and E available [kWh] is the maximum energy available to be recovered during braking. The available energy is calculated using Equation (4).
where P Wheels (−) (t) is the negative part of the power in the wheels in [kW].

Electric Motor Model
The output power of the alternating current electric motor depends on its efficiency, which is modeled as a quasi-static efficiency map that depends on the output speed and torque, with the efficiency of the inverter considered constant. With these premises, the output power (P em ) can be calculated by Equation (5) [35].
where P inv is the power of the inverter (W), W em is the angular speed of the electric motor (rad/s), T em is the torque of the electric motor (Nm), and η em is the electric motor efficiency (%). During simulation, the electric motor controller receives a power demand and adjusts the torque output to supply it. Since the current speed of the motor cannot change instantaneously, the theoretical torque is calculated by the controller based on the current speed of the motor and the power demand, and then finds the current efficiency of the motor from the efficiency map. However, it is also possible to resort to specific parameters of the electric motor by replacing the efficiency map with values of maximum torque and maximum power-see Equation (5)-and considering the motor loss constants (k c , k i , and k w ) for the loss of power, power output is limited by the capacity of the electric motor, with the following equations: P in = T·W + P loss (8)

Transmission Model
The drivetrain system transmits the movement obtained by the electric motor to the wheels as shown in Figure 7, which is modeled by (11) and (12).
where η g is the transmission efficiency (%), F t is the pulling force (N), G is the transmission ratio [36].

Transmission Model
The drivetrain system transmits the movement obtained by the electric motor to the wheels as shown in Figure 7, which is modeled by (11) and (12).
where η g is the transmission efficiency (%), F t is the pulling force (N), G is the transmission ratio [36]. Equation (11) is used when the electric motor is delivering mechanical power (T el-em and Equation (12) is used when the electric motor is receiving mechanical power (T em-el [37]. Likewise, the relationship between the linear speed "v" of the electric vehicle and the angular speed Wem of the motor in (rad/s) can be obtained by Equation (13).
Finally, the transmission force is obtained with Equation (14).
where T Loss is the friction losses associated with the torque.

Battery Model
The battery modeling can be addressed in various ways, the most used being the electric and physical electrochemical models [38,39]. The equivalent electrical model is simple to build and adjust, giving good results in estimating voltage, current, SOC, and heat rate [40]. On the other hand, it is possible to express the physical model of the battery based on electrochemical phenomena. However, this method can be complex and computationally expensive, and is also commonly used to evaluate its behavior outside the predetermined temperature range of the batteries. Another key aspect of this modeling is the ability to use models of physics-based (as opposed to empirical) aging that give insight into the growth rate of the cathode film, the lithium coating, and the insulation of the material due to cracking [41,42]. In summary, the objective is to estimate the energy consumption of the EV throughout a driving cycle but not to estimate the aging of each battery cell or the degradation of the performance of the system over time. Thus, it could be more feasible to use the equivalent electrical model for battery [43]. Through the equivalent electrical model, this type of charging and discharging behavior is represented in an Equation (11) is used when the electric motor is delivering mechanical power (T el−em ) and Equation (12) is used when the electric motor is receiving mechanical power (T em−el ) [37]. Likewise, the relationship between the linear speed "v" of the electric vehicle and the angular speed W em of the motor in (rad/s) can be obtained by Equation (13).
Finally, the transmission force is obtained with Equation (14).
where T Loss is the friction losses associated with the torque.

Battery Model
The battery modeling can be addressed in various ways, the most used being the electric and physical electrochemical models [38,39]. The equivalent electrical model is simple to build and adjust, giving good results in estimating voltage, current, SOC, and heat rate [40]. On the other hand, it is possible to express the physical model of the battery based on electrochemical phenomena. However, this method can be complex and computationally expensive, and is also commonly used to evaluate its behavior outside the predetermined temperature range of the batteries. Another key aspect of this modeling is the ability to use models of physics-based (as opposed to empirical) aging that give insight into the growth rate of the cathode film, the lithium coating, and the insulation of the material due to cracking [41,42]. In summary, the objective is to estimate the energy consumption of the EV throughout a driving cycle but not to estimate the aging of each battery cell or the degradation of the performance of the system over time. Thus, it could be more feasible to use the equivalent electrical model for battery [43]. Through the equivalent electrical model, this type of charging and discharging behavior is represented in an electrical circuit using the resistive battery model, as shown in Figure 8. The battery model in the simulation is not subjected to a deep discharge (below 20%) nor until all energy is exhausted (0%); it is only subjected to the analysis of the state of charge at the end of the UDDS cycle. Electricity 2022, 3, FOR PEER REVIEW 8 electrical circuit using the resistive battery model, as shown in Figure 8. The battery model in the simulation is not subjected to a deep discharge (below 20%) nor until all energy is exhausted (0%); it is only subjected to the analysis of the state of charge at the end of the UDDS cycle. The maximum battery power denoted by pmax points out the maximum power that the storage set can exchange with the system, and is given by Equation (15).
where V oc is the open circuit voltage (V), V ve is the minimum motor controller voltage or the minimum battery voltage, and R int is the charge and discharge resistance. The charge or discharge current is calculated with Equation (16).
where P is the charging/discharging power of the batteries. The SOC of the batteries can be estimated using Equation (17).
where Ah max is the maximum capacity of the battery bank and Ah used is the used capacity of the battery bank.

EV Global Model
Longitudinal Dynamics For EV modeling, the inertia and resistances of the vehicle must be represented. Longitudinal dynamics analyzes the behavior of the vehicle when it is subjected to speed fluctuations in the longitudinal plane, neglecting the lateral acceleration [44]. In this case, the vehicle moves in a straight line, accelerating, braking, or changing the slope. Longitudinal dynamics is affected by different resistances, such as rolling, aerodynamics, and slope resistance [45]. To analyze these factors, a vehicle moving on an inclined plane and in a straight line is considered.
The aerodynamic drag is a force acting on objects moving in a fluid opposite to their motions, which is calculated with Equation (18) [27].
where ρ is the density of air, A is the frontal area of the EV, C d is the drag coefficient that depends on the shape of the vehicle for a particular vehicle C d and varies between 0.3 and 0.4, and v is the speed of the EV. The maximum battery power denoted by p max points out the maximum power that the storage set can exchange with the system, and is given by Equation (15).
where V oc is the open circuit voltage (V), V ve is the minimum motor controller voltage or the minimum battery voltage, and R int is the charge and discharge resistance. The charge or discharge current is calculated with Equation (16).
where P is the charging/discharging power of the batteries. The SOC of the batteries can be estimated using Equation (17).
where Ah max is the maximum capacity of the battery bank and Ah used is the used capacity of the battery bank.

EV Global Model Longitudinal Dynamics
For EV modeling, the inertia and resistances of the vehicle must be represented. Longitudinal dynamics analyzes the behavior of the vehicle when it is subjected to speed fluctuations in the longitudinal plane, neglecting the lateral acceleration [44]. In this case, the vehicle moves in a straight line, accelerating, braking, or changing the slope. Longitudinal dynamics is affected by different resistances, such as rolling, aerodynamics, and slope resistance [45]. To analyze these factors, a vehicle moving on an inclined plane and in a straight line is considered.
The aerodynamic drag is a force acting on objects moving in a fluid opposite to their motions, which is calculated with Equation (18) [27].
where ρ is the density of air, A is the frontal area of the EV, C d is the drag coefficient that depends on the shape of the vehicle for a particular vehicle C d and varies between 0.3 and 0.4, and v is the speed of the EV. Two other parameters are of importance when modeling the longitudinal dynamics: the rolling and slope resistance. The first coefficient refers to the resistance force of the rotational movement of the tire in contact with the surface, which is denoted by F rr and can be calculated with Equation (19), while the slope resistance is denoted by F rp and is given by (20). Finally, it is worth mentioning that the pulling force F t is given by (21) [27,46].
where m is the mass of the EV, g is the gravity acceleration, C r is the coefficient of friction, θ is the angle of inclination between the EV and the surface, and a is the acceleration of the EV.

Multibody
In this section, the multibody modeling that represents each individual component of the EV based on their lateral and vertical physical properties is analyzed. Three degrees of freedom are considered for this theoretical study; however, it is worth mentioning that other, more sophisticated models can represent up to 94 degrees [47]. To express the complete block of the multibody, in this paper, the mathematical models of the vehicle frame, suspension, and wheels are considered.
The frame vehicle is modeled as a rigid body on which the translational displacement in the three directions of the fixed coordinate system (x, y, z) and the rotary movement around the three axes (ϕ, χ, ψ) are represented, in order to obtain the output vector X = [x y z ϕ χ ψ] T . These movements occur due to the effect of rolling resistance F fr R ij , tire forces F fr x ij , F fr y ij , and tire forces on the frame suspension F fr Zs ij [48]. For the calculation of the angular displacement, Euler's Law is applied. In addition, the angular accelerations of the vehicle can be obtained through the centroidal inertia matrix and the moments exerted by the forces on the body [49]. To obtain the angular accelerations of EV (w p ), Equation (22) is applied.
where · E is the Euler derivate function, ϕ represents the rotation on the x-axis, χ represents the rotation on the y-axis, and ψ represents the rotation on the z-axis. For the suspension model, the position vectors r cg s ij must initially be established for each vertex R s ij with respect to the coordinate system of the center of gravity, with the aim of determining the overall displacement, as shown in Figure 9. The reaction force of the suspension on the frame can be calculated by using Equation (23).
Finally, the mathematical modeling of the wheels complements the multibody model. In this sense, the behavior of the wheels significantly affects the performance of the vehicle due to the horizontal forces and moments generated [50]. In addition, the wheels support other types of efforts by having to hold up the weight of the vehicle, as well as the forces of aerodynamic load and mass transfer during braking or turning of the vehicle [49]. Experimental methods are used to simplify and model the tire. One of the most used and well-known is the Pacejka model. This model relates the longitudinal force (F x ), the lateral force F y , and the alignment moment (M z ) with the drift angle (α) and the drift ratio (κ) of the wheels by means of a set of formulas to interpolate all data. First, the centers of rotation of the wheels Q ij are defined in the horizontal plane with respect to the origin of the frame coordinate system by means of the position vectors q f r ij . Figure 10 details the dimensions of the vehicle [50]. Finally, the mathematical modeling of the wheels complements the multibody model. In this sense, the behavior of the wheels significantly affects the performance of the vehicle due to the horizontal forces and moments generated [50]. In addition, the wheels support other types of efforts by having to hold up the weight of the vehicle, as well as the forces of aerodynamic load and mass transfer during braking or turning of the vehicle [49]. Experimental methods are used to simplify and model the tire. One of the most used and well-known is the Pacejka model. This model relates the longitudinal force , the lateral force , and the alignment moment with the drift angle and the drift ratio of the wheels by means of a set of formulas to interpolate all data. First, the centers of rotation of the wheels ( ) are defined in the horizontal plane with respect to the origin of the frame coordinate system by means of the position vectors . Figure  10 details the dimensions of the vehicle [50].   Finally, the mathematical modeling of the wheels complements the multibody model. In this sense, the behavior of the wheels significantly affects the performance of the vehicle due to the horizontal forces and moments generated [50]. In addition, the wheels support other types of efforts by having to hold up the weight of the vehicle, as well as the forces of aerodynamic load and mass transfer during braking or turning of the vehicle [49]. Experimental methods are used to simplify and model the tire. One of the most used and well-known is the Pacejka model. This model relates the longitudinal force , the lateral force , and the alignment moment with the drift angle and the drift ratio of the wheels by means of a set of formulas to interpolate all data. First, the centers of rotation of the wheels ( ) are defined in the horizontal plane with respect to the origin of the frame coordinate system by means of the position vectors . Figure  10 details the dimensions of the vehicle [50].  Lastly, it is necessary to rotate the horizontal forces of the front wheels by the angle of twist δ w around the z-axis, as shown in Equation (24). To sum up, the different sub-blocks involved in EV modeling are depicted in Figure 11, for simplicity.
where F fr su ij represents the position vector of each component applied from the suspension on the chassis and the wheel, respectively, F fr w ij represents the position vector of each component applied from the wheels on the chassis, and Ω ij is the angular velocity of the wheels. where F su ij fr represents the position vector of each component applied from the suspension on the chassis and the wheel, respectively, F w ij fr represents the position vector of each component applied from the wheels on the chassis, and Ω is the angular velocity of the wheels. Figure 11. EV model with the corresponding subsystems.

Simulation
A modeling of the EV is developed with the forward method approach (Figure 12), which comprises all the blocks described throughout the study, using longitudinal dynamics for vehicle modeling. The cycle used is the UDDS (see Figure 3b) and the input values are collected in Table 1.  Figure 11. EV model with the corresponding subsystems.

Simulation
A modeling of the EV is developed with the forward method approach (Figure 12), which comprises all the blocks described throughout the study, using longitudinal dynamics for vehicle modeling. The cycle used is the UDDS (see Figure 3b) and the input values are collected in Table 1.

Results and Discussion
This section presents the results and discussion of the simulations performed on Matlab-Simulink, in which all the parameters and forces interacting on the EV and the energy of the vehicle at the end of the imposed driving cycle can be estimated, as observed in Figure 13. The EV profile is similar to the UDDS route, which is passed to the driver model. The PID controller converts the error into signals for the accelerator and brake pedal. Therefore, its representation is observed in detail throughout the driving cycle, as seen in Figure 14. Figure 15 shows the SOC of the battery bank throughout the driving cycle. The initial state of charge was considered equal to 80.7% and, after completing the route (1369 s), declined to 74.3%, being therefore reduced by 6.4%. Figure 16 shows the battery voltage variation thorughout the driving cycle. The average voltage required is 340 V, with the maximum capacity of the battery being 350 V, a capacity that is not exceeded during the driving cycle.

Results and Discussion
This section presents the results and discussion of the simulations performed on Matlab-Simulink, in which all the parameters and forces interacting on the EV and the energy of the vehicle at the end of the imposed driving cycle can be estimated, as observed in Figure 13. The EV profile is similar to the UDDS route, which is passed to the driver model. The PID controller converts the error into signals for the accelerator and brake pedal. Therefore, its representation is observed in detail throughout the driving cycle, as seen in Figure 14. Figure 15 shows the SOC of the battery bank throughout the driving cycle. The initial state of charge was considered equal to 80.7% and, after completing the route (1369 s), declined to 74.3%, being therefore reduced by 6.4%. Figure 16 shows the battery voltage variation thorughout the driving cycle. The average voltage required is 340 V, with the maximum capacity of the battery being 350 V, a capacity that is not exceeded during the driving cycle.

Results and Discussion
This section presents the results and discussion of the simulations performed on Matlab-Simulink, in which all the parameters and forces interacting on the EV and the energy of the vehicle at the end of the imposed driving cycle can be estimated, as observed in Figure 13. The EV profile is similar to the UDDS route, which is passed to the driver model. The PID controller converts the error into signals for the accelerator and brake pedal. Therefore, its representation is observed in detail throughout the driving cycle, as seen in Figure 14. Figure 15 shows the SOC of the battery bank throughout the driving cycle. The initial state of charge was considered equal to 80.7% and, after completing the route (1369 s), declined to 74.3%, being therefore reduced by 6.4%. Figure 16 shows the battery voltage variation thorughout the driving cycle. The average voltage required is 340 V, with the maximum capacity of the battery being 350 V, a capacity that is not exceeded during the driving cycle.

Results and Discussion
This section presents the results and discussion of the simulations performed on Matlab-Simulink, in which all the parameters and forces interacting on the EV and the energy of the vehicle at the end of the imposed driving cycle can be estimated, as observed in Figure 13. The EV profile is similar to the UDDS route, which is passed to the driver model. The PID controller converts the error into signals for the accelerator and brake pedal. Therefore, its representation is observed in detail throughout the driving cycle, as seen in Figure 14. Figure 15 shows the SOC of the battery bank throughout the driving cycle. The initial state of charge was considered equal to 80.7% and, after completing the route (1369 s), declined to 74.3%, being therefore reduced by 6.4%. Figure 16 shows the battery voltage variation thorughout the driving cycle. The average voltage required is 340 V, with the maximum capacity of the battery being 350 V, a capacity that is not exceeded during the driving cycle.     Figure 17 shows the battery current variation thorughout the route. The average current required is 50 A; the current is proportional to the speed, and if the speed increases, more current is demanded. In Figure 18, a comparison of the variation in the speed and the variation in the battery power is shown. It can be seen that the velocity changes are proportional to the battery power. Table 2 shows the energy results of the EV using the Electricity 2022, 3 214 method proposed for the UDDS cycle; for a distance of approximately 12 km, the battery has been discharged by almost 25%.   Figure 17 shows the battery current variation thorughout the route. The average current required is 50 A; the current is proportional to the speed, and if the speed increases, more current is demanded. In Figure 18, a comparison of the variation in the speed and the variation in the battery power is shown. It can be seen that the velocity changes are proportional to the battery power. Table 2 shows the energy results of the EV using the method proposed for the UDDS cycle; for a distance of approximately 12 km, the battery has been discharged by almost 25%.     Figure 17 shows the battery current variation thorughout the route. The average current required is 50 A; the current is proportional to the speed, and if the speed increases, more current is demanded. In Figure 18, a comparison of the variation in the speed and the variation in the battery power is shown. It can be seen that the velocity changes are proportional to the battery power. Table 2 shows the energy results of the EV using the method proposed for the UDDS cycle; for a distance of approximately 12 km, the battery has been discharged by almost 25%.     Undoubtedly, this is a method used by many studies for the modeling and simulation of an EV for charging and discharging estimation throughout the driving cycle. This study focuses on the model development process in general, and not on developing each component in its complex form, as in Ref. [35], where electric motors are designed and analyzed for applications in plug-in EVs, or in [51,52], where the regenerative brake is studied and modeled. Similarly, in Ref. [53], the authors model other complementary and important components, such as the converter DC/DC and inverter. After understanding the process of forward modeling, it is possible to include it in the integral simulation chain, replacing the basic components, including energy management systems, by analyzing the effect of temperature on a system's energy storage [54]; thus, it will be possible to obtain more accurate and realistic estimates.
The multibody modeling of the vehicle in the actual literature is neglected, not providing much difference compared to the longitudinal modeling, when estimating the energy consumption. The multibody modeling of the vehicle depends on the results that are required; on the other hand, this modeling is ideal for analyzing the dynamic performance of each component (suspension, chassis, and tire) and even for analyzing the lateral behavior of the steering system. The multibody modeling involves sub-models that depend on the degrees of freedom, its modeling being much more complex. In Ref. [47], the authors analyze a model with 94 degrees of freedom. Longitudinal dynamics is favorable for estimating energy consumption, because it can include a slope profile together with the speed profile to be input parameters to the driver's model, and it can also incorporate the consumption of auxiliary elements of the EV, such as electric motors for the steering system, air conditioning, pumps, etc.
When comparing the characteristics of the backward and forward methods, a peculiarity of the latter stands out, since the acceleration of the vehicle is expressed as a function of the accelerator position, unlike the backward method, where the force required to accelerate the vehicle is determined only from the speed cycle. For this reason, the forward method is a more complex and accurate method, which serves to develop real-time eco-driving systems and eco-routes that can improve energy efficiency to extend the EV's travel range, as studied in [43]. Table 3 compares autonomy with the simulation of a more accurate model from the Matlab/Simulink library, adjusting the input values of Table 1, and also with [14], a reference used for the validation of the model. The model developed in this study has higher energy consumption per hour than the Matlab models and [14]. These EV models (Matlab, Ref. [14]) are more complex and apply control algorithms for the drivetrain and for battery energy management; it incorporates components such as BMS, sensors, regenerative braking, etc., leading to greater accuracy in estimating energy consumption. The battery charge status for the EV model from the matlab/simulink library is shown in Figure 19 and for the EV model from [14] it is shown in Figure 20. The three models are developed by the same method (forward), but there is a large difference in energy consumption, especially in the model developed in this study. It should be noted that this study focuses on highlighting the development process of the EV model by the forward method, but not on developing a model that improves consumption autonomy. For the estimation of autonomy, values from the United States Environmental Protection Agency (US EPA) are used. The model developed in this study has higher energy consumption per hour than the Matlab models and [14]. These EV models (Matlab, Ref. [14]) are more complex and apply control algorithms for the drivetrain and for battery energy management; it incorporates components such as BMS, sensors, regenerative braking, etc., leading to greater accuracy in estimating energy consumption. The battery charge status for the EV model from the matlab/simulink library is shown in Figure 19 and for the EV model from [14] it is shown in Figure 20. The three models are developed by the same method (forward), but there is a large difference in energy consumption, especially in the model developed in this study. It should be noted that this study focuses on highlighting the development process of the EV model by the forward method, but not on developing a model that improves consumption autonomy. For the estimation of autonomy, values from the United States Environmental Protection Agency (US EPA) are used.

Conclusions
A detailed study of the modeling and simulation of the EV is presented considering the forward method for its development, highlighting the intrinsic characteristics. One of the aspects that clearly differs compared to the backward method is the model of the driver, who acts as an acceleration and braking pedal to follow the imposed target speed. Driver modeling can be achieved by different controllers: the classic PID, Pole Location Controller, Linear Quadratic Regulator (LQR), and Observer Based Controller (OBC). LQR is one of the controllers that best adjusts the throttle position; therefore, it is compared with the PID controller, looking at their response time to a simple step. The LQR controller improves the overshoot, while properly adjusting the PID controller reduces the response and stabilization time; for this study, a PID controller is used, where the results do not show a greater error between the target speed and the real speed of the vehi-

Conclusions
A detailed study of the modeling and simulation of the EV is presented considering the forward method for its development, highlighting the intrinsic characteristics. One of the aspects that clearly differs compared to the backward method is the model of the driver, who acts as an acceleration and braking pedal to follow the imposed target speed. Driver modeling can be achieved by different controllers: the classic PID, Pole Location Controller, Linear Quadratic Regulator (LQR), and Observer Based Controller (OBC). LQR is one of the controllers that best adjusts the throttle position; therefore, it is compared with the PID controller, looking at their response time to a simple step. The LQR controller improves the overshoot, while properly adjusting the PID controller reduces the response and stabilization time; for this study, a PID controller is used, where the results do not show a greater error between the target speed and the real speed of the vehicle. The use of one controller or another depends on the results required and the type of vehicle.
The results of the forward method can estimate the energy consumption of the EV, as well as the energy losses in the components. In this particular case, the input data (Table 2) were obtained as a result after fulfilling the driving cycle UDDS 1396 s, with a distance traveled of 11.99 km. The SOC dropped from 80.7% to 74.3%, with energy consumption of 2209 kWh and energy losses in some components, such as the electric motor, of 0.429 kWh, the transmission of 0.221 kWh, and the battery of 0.0324 kWh.
The simulation was also compared with other more precise models ( Table 3); one of them was an existing model in the Matlab library (Powertrain Blockset version 1.10), and we also compared it with the reference [14], which uses the AUTONOMIE simulator. For all of them, the same values are adjusted. From the outset, where there is a notable difference is in the autonomy; the model developed in this study consumes more energy, having an autonomy of 4.45 kWh/100 km, and the one that consumes the least is the one that uses the AUTONOMIE simulator, with a consumption of 0.97 kWh/100 km.
Future research can focus on these concepts to develop computational models based on artificial neural networks (ANN) to determine the relationships between a series of factors that affect the energy consumption of EVs depending on the input factors, where it is determined a weight for each factor based on its relative importance using training algorithms. ANN can also be used to predict driving behavior by classifying driving patterns using Global Positioning System data, being an extremely important approach method as it is data-driven and self-adaptive.