Advanced
Simulation of Electric Vehicles Combining Structural and Functional Approaches
Simulation of Electric Vehicles Combining Structural and Functional Approaches
Journal of Electrical Engineering and Technology. 2014. May, 9(3): 848-858
Copyright © 2014, The Korean Institute of Electrical Engineers
  • Received : May 24, 2013
  • Accepted : October 22, 2013
  • Published : May 01, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
L.I. Silva
Grupo de Electrónica Aplicada, Dpto de Electricidad y Electrónica, Facultad de Ingeniería, Universidad Nacional de Rio Cuarto, Argentina
G.A. Magallán
Grupo de Electrónica Aplicada, Dpto de Electricidad y Electrónica, Facultad de Ingeniería, Universidad Nacional de Rio Cuarto, Argentina
P.M. De la Barrera
Grupo de Electrónica Aplicada, Dpto de Electricidad y Electrónica, Facultad de Ingeniería, Universidad Nacional de Rio Cuarto, Argentina
C.H. De Angelo
Grupo de Electrónica Aplicada, Dpto de Electricidad y Electrónica, Facultad de Ingeniería, Universidad Nacional de Rio Cuarto, Argentina
G.O. García
Grupo de Electrónica Aplicada, Dpto de Electricidad y Electrónica, Facultad de Ingeniería, Universidad Nacional de Rio Cuarto, Argentina

Abstract
In this paper the construction of a model that represents the behavior of an Electric Vehicle is described. Both the mechanical and the electric traction systems are represented using Multi-Bond Graph structural approach suited to model large scale physical systems. Then the model of the controllers, represented with a functional approach, is included giving rise to an integrated model which exploits the advantages of both approaches. Simulation and experimental results are aimed to illustrate the electromechanical interaction and to validate the proposal.
Keywords
1. Introduction
Transportation solutions have been thoroughly investigated to improve fuel economy and reduce vehicle emissions. Electric Vehicles (EV) have been proposed as possible solutions to meet these goals, see Fig. 1 .
PPT Slide
Lager Image
(a) Vehicle prototype; (b) Vehicle geometry in CAD design.
Modeling and Simulation are key tasks for testing novel concepts. They help in the design stage allowing the evaluation of different vehicle configurations, inexpensively and without implementing a prototype. A well detailed revision of vehicle models and their features is given in [1] . An overview of EVs with a focus on hybrid configurations and energy management is discussed in [2] .
Vehicle models are classified as steady-state, dynamic and quasi-steady-state. In steady-state models transients are neglected and require less computation power than dynamic models which are physic-based and reproduce transient evolutions. An intermediate solution is a quasi-steady-state model. Another classification is based on the calculation direction: Forward models follow the engine-to-wheel direction whereas backward models follow the wheel-to-engine relations. Generally, dynamic models are used in forward calculation and hardware-in-the-loop simulation [3] .
Once the model of the vehicle is defined according to the application, the representation of the model can be accomplished with a functional or structural approach [4] . In a functional approach the connections between subsystems are virtual links. It gives an overall description of the system and is used in applications oriented to its control. On the other hand, structural approaches are based on the topology of the real system and the links respond to a physical connection. Hence a direct correspondence between the system under analysis and its representation is generated. Structural approaches are useful for system design and its analysis in faulty conditions.
In applications such as those addressed in this paper, it is convenient to resort to energy-based approaches. They explicitly define the energy flow and/or power interaction between components. Some energy-based approaches are:
Bond graphs (BG) allows interconnecting sub-systems based on their power flow [5] and facilitates the intercomnection of multi domain systems [6] . BG contains dissipative, accumulative, conversion elements and junctions. Elements are interconnected by bonds that carry the exchange variables (general flow and effort) and allow an easy visualization of the physical system topology and power flow.
Multi Bond graphs (MBG) [7] is a vectorial extension of standard BG. It is also a non-causal approach suited to model large scale multi domain physical systems.
Power Flow Diagrams (PFD) [8] is based on splitting the bond to differentiate flow and effort variables to better interpret the system. New pictograms are defined to highlight the different elements in the system. PFDs have been used to study electric drives and traction systems [9] .
Energetic Macroscopic Representation (EMR) [10] makes apparent the power interchange between subsystems. It is composed of connected elements that highlight energy properties of the system. It allows determining systematically an inversion based control scheme. In [11] this approach is used to model and control HEVs.
After defining the vehicle model and its representation, different analyses can be conducted via simulation. Within the large variety of simulation tools devoted to EVs and HEVs, Advisor [12] and QSS-TB [13] are aimed to develop quasi steady-state models and have been used successfully to evaluate fuel consumption and its optimization [14] .
In [15] it is used Matlab/Simulink to construct dynamic models which is a flexible simulation tool where the user can define and program each component of the vehicle.
More specific tools incorporate libraries oriented to model EVs and HEVs such as Dymola [16] and PSAT [17] .
Finally there are specific softwares in which the modeling is carried out by interconnecting predefined blocks making the procedure faster but less flexible. In [18] it is used PSIM whereas in [19] it is used Carsim.
The goal of this paper is to develop a complete 3-D representation of a four wheels EV (see Fig. 1 ) using a structural energy-based approach whereas its control is represented using a functional description. In this way, the advantages of each approach are better exploited. This idea is inspired by [20] where a first coupling between a structural and a functional approach is proposed.
In this work a power system is simulated using PSIM (with a structural component library) and the control is designed using the EMR (functional approach). Also in [21] it is studied the advantage of using BG for representing the system and block diagrams for representing the controller.
Here the BG model is only used to obtain the dynamic equations of the powertrain but then the entire system (vehicle and control) is simulated in Matlab / Simulink using block diagrams (i.e. the structural approach is lost).
To meet the goal, in our proposal, the 3-D mechanics and the electric traction system of the vehicle are modeled with MBG and its control is described in block diagrams. In order to simulate the complete system behavior, the vehicle model and its control are integrated in the same simulation environment.
Models of the 3-D mechanical behavior of vehicles using BG were already presented in [22] . In [23] the translational and rotational dynamics are grouped and modeled together by mean of the MBG reducing significantly these models but the electromechanical interaction was neither modeled nor analyzed. To include this interaction, BG-based models of EVs and HEVs were presented in [24] . They include a BG model of the electric motor but the mechanical behavior is reduced to longitudinal dynamics only. Thus the proposed model in this work complements those models presented in recent literature.
Nomenclature
PPT Slide
Lager Image
2. Mechanical Characteristics and Modeling
In this Section, the model of the vehicle is created using Dymola which is an integrated environment for developing models in the open Modelica language [25] . Dymola allows generating structural representations which consist of connected components. These components have “connectors” that describe the interaction possibilities, e.g. an electrical pin, a mechanical flange, or an input/output signal. By connecting components a physical system model can be constructed. Dymola includes libraries to generate BG and MBG models. In 3-D mechanical applications, elements are represented compactly using MBG. Three single bonds are grouped into a Multi-Bond representing the rotational dynamics and another Multi-Bond is used for describing the translational dynamics.
Chassis and suspensions are modeled using the “Modelica library for MBGs and its application in 3-D mechanics” [26] . These components are internally represented with elements of the MBG library and are provided with external frames where the rotational and translational Multi-bonds are connected. The mechanical connection between two 3-D components is made via their frames that act as interfaces for exchanging power. As example, the upper part of Fig. 2 depicts the graphical representation of the connection of two 3-D components called “fixed translations” while the lower part shows that their internal models uses the MBG approach. The connection implies that both components have the same angular / linear velocities and positions at the connecting point.
PPT Slide
Lager Image
Connection of two different 3-D mechanical elements.
- 2.1. Chassis
The model of the chassis contains a sprung mass that behaves like a lumped mass for the translational dynamics whereas for the rotational dynamics it contains the inertia tensor thus acting like a rigid body. The laws for the 3-D body motion are given by the Euler’s equations
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
where F = [ Fx Fy Fz]T and T = [ Tx Ty Tz]T are the forces (in the global frame) and torques (in the local frame) applied over the rigid body, respectively. Also, m and a represent the mass and acceleration of the chassis. Moreover, L , 𝜔 and ICH represent the angular momentum, angular velocity and moment of inertia tensor of the chassis, respectively. Fig. 3 shows the 3-D body representation together with its internal representation in MBG.
PPT Slide
Lager Image
(a) Free body in 3-D mechanics, (b) Internal composition in MBG.
To complete the chassis model “Fixed translations” are needed. They are ideal undeformable rods considered massless. Forces and velocities are given in the global frame whereas torques and angular velocities in the local frame. The set of equations is given by
PPT Slide
Lager Image
PPT Slide
Lager Image
where R and r are the orthonormal orientation matrix and the vector that links both frames, respectively. The representation in 3-D mechanics is depicted in Fig. 2 and its internal MBG representation is presented in Fig. 4 . Eq. (4) is graphically expressed in the “0” of the upper part whereas Eq. (5) is in the “1” of the lower part.
PPT Slide
Lager Image
Internal representation of fixed translations in MBG.
Fig. 5 shows the complete chassis model. Here a set of fixed translations is used to compose the undeformable structure of the chassis that projects the efforts and flows from the sprung mass -located in the center of gravity toward the points where the suspensions are mechanically connected in the frames noted with ① , ②, ③ and ④.
PPT Slide
Lager Image
Complete chassis modeled with 3-D mechanics components.
- 2.2. Suspension
The suspension system is composed by four suspensions that couple each wheel to the chassis. Fig. 6(a) shows that each suspension contains an ideal spring and a damper (from the 3-D mechanics library) which are connected according to the real mechanical configuration.
PPT Slide
Lager Image
(a) 3-D mechanical model, (b) Internal composition in MBG.
Spring force ( Fk ), damper force ( FB ) and total suspension force ( Fsusp ) are given by
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
where velocities on frames 1 and 2 are v 1 and v 2 and Δx is the displacement respect to the rest position ( x 0 ).
The suspension stiffness and damper coefficient matrices are k S = − kS I 3x3 and B S = BS I 3x3 , respectively. Fig. 6(b) depicts the internal model in MBG.
- 2.3. Wheel and tires
Fig. 7(a) contains the Dymola internal model of the wheel that includes the calculation of the normal force at the contact patch ( FN ) based on a BG model of the vertical dynamics.
PPT Slide
Lager Image
(a) BG-MBG internal model of the wheel, (b) 3- Dmechanical representation.
The longitudinal / lateral forces ( Fx / Fy ) depend on the tire characteristics and certainly influence the dynamic behavior of the vehicle. These forces are determined using a 2-D LuGre friction model [27] that is widely used to estimate and control the traction contact force [28] . It allows representing different road conditions and the transient behavior of the traction forces with a concise model that considers the coupling between x and y directions. Rotational dynamics is represented in a BG model that receives the torque ( Te ) generated by the Induction Machine (IM) and calculates its angular speed ( 𝜔r ). The torque balance also includes: rolling resistance, inertia momentum of the wheel ( IW ) and the opposition torque due to the longitudinal force on the patch contact at distance rW . Then, translational dynamics along three axes are grouped into a multi-bond, as well as rotational dynamics.
These multi-bonds are used to incorporate the 3-D mechanical frame. Fig. 7(b) shows the Dymola 3-D mechanics model. It contains a mechanical frame to connect with the suspension, a single bond to interact with the IM and receives the steering signal ( δ ).
3. Electric Traction System
Two identical IM drive the rear wheels which are controlled independently. This feature allows the implementation of yaw stability control [29 , 30] using a simple strategy to estimate the sideslip angle [31] .
The set of equations concerning voltage and flux linkage in qd coordinates together with the equivalent circuit are given in [32] . The qd coordinates representation means that all variables are referred to an arbitrary reference frame rotating with angular velocity 𝜔 qd0 = dρ/dt . Angle ρ is the d-axis position. Electromagnetic torque ( Te ) as function of electrical variables, obtained with the power invariant transformation, is given by
PPT Slide
Lager Image
to represent the power interaction between electrical and mechanical domains modulated gyrators ( MGY ) were used, see Fig. 8 . From Fig. 8(a) the following equation can be derived
PPT Slide
Lager Image
BG diagram used to couple electrical and mechanical domains.
PPT Slide
Lager Image
Taking into account that the MGY ’s do not add or dissipate power and using the power conservation at the “1” the following equation can be obtained
PPT Slide
Lager Image
Replacing (10) into (11) an alternative expression for Te is
PPT Slide
Lager Image
It can be proved that Eq. (12) is equivalent to (9), thus the MGY presented in Fig. 8(b) can be used to couple the electrical and mechanical domains.
To keep consistency with the mechanical model, the IM equivalent circuit was represented in a MBG of cardinality two. The resulting model is presented in the compact diagram shown in Fig. 9(a) together with its Dymola representation in Fig. 9(b) . Variables and parameters are
PPT Slide
Lager Image
(a) MBG representation; (b) Wrapped model in Dymola.
PPT Slide
Lager Image
For the case of study, electric variables are expressed using a qd stationary frame ( 𝜔 qd0 = 0). Variables in this frame are noted as αβ variables.
4. Integrated Model
This Section is devoted to generate the integrated model which includes in the same simulation environment the complete model of the vehicle (with a structural representation in Section IV-A) and the controllers (with a functional representation in Section IV-B).
- 4.1 Complete model of the vehicle
In order to obtain the structural representation of the entire vehicle, the mechanical and electrical sub-systems were interconnected. Each suspension and the chassis are connected though a revolute joint, which permits relative rotation along one axis. It also allows the chassis pitch movement while each suspension remains in vertical position. Similarly the lower end of the suspension system connects to the wheel’s frame through a revolute joint. This connection allows the spinning of the wheel. Fig. 10 shows the complete connection of the left rear wheel and suspension. Connections of the other three wheels are analogous (front wheels do not have traction motor). It was also incorporated a non-lineal resistance “ rnl ” that models the dragging force along the longitudinal direction.
PPT Slide
Lager Image
Complete connection of Wheel and Suspension in Dymola.
- 4.2 Control Scheme
The most adequate approach used to describe a control scheme is a functional representation. This Section presents the scheme implemented in the prototype which gives a clear overall description. Additionally, each block is internally programmed giving hierarchical levels of descriptions.
In order to emulate a mechanical differential, an equal torque strategy is proposed. This strategy is implemented by the Electronic Differential System (EDS), see Fig. 11 , which determines the reference torque to be applied by each IM. Based on this reference torque a Field Oriented Control (FOC) was implemented on each motor, see Fig. 12 . A well detailed description of the EDS/FOC strategy is given in [33] . It is proposed a simple strategy. More complex controllers regarding an EDS are presented in [34] .
PPT Slide
Lager Image
Block Diagram representation of the EDS.
PPT Slide
Lager Image
Block Diagram representation of the FOC.
- 4.3 Integration of structural / functional representations
The final part of the proposal consists on coupling the vehicle model with its control in the same simulation environment in order to study the complete vehicle performance.
The physical system model, built in Dymola, is compiled by Modelica and exported to the Matlab/ Simulink environment. To this end, the model (see Fig. 10 ) is provided with the inputs / outputs necessary to interact in Simulink. Inputs are the steering angle ( δ ) and the reference αβ voltages applied to the IMs. Currents and rotor angular speed are the outputs used to implement both FOCs. The control algorithm is obtained in Simulink to close the control loop.
The integrated model which combines both approaches is displayed in Fig. 13 . It shows explicitly that the inputs are the driver commands and the outputs can be any physical variable of interest.
PPT Slide
Lager Image
The resulting model allows evaluating different controllers in a simulation environment, but using a complete vehicle model. Thus, any modification in the controller configuration or in its adjustment can be easily tested under simulation. Once the control algorithm has been tuned, it can be transferred to the digital controller using toolboxes that allow directly translating the control algorithm to a particular digital processor [35] . These toolboxes generate optimized C code for DSP / FPGA from plain Simulink diagrams.
Additionally, BG allows us to modify quickly the physical model of the vehicle before producing any change in the actual prototype. Thus, the impact of such changes (e.g. tires, suspensions) over the controller or the overall vehicle performance can be verified and evaluated immediately, at practically no cost. Once the modifications have been tested and the results are satisfactory, they can be implemented in the actual vehicle prototype.
- 5. Simulations Results and Experimental Validations
The first part (Section V-A) provides simulation and experimental results that validate the proposed model. Then in Section V-B only simulation results are given in order to highlight potential applications of the proposed methodology. These results are obtained using the scheme presented in Fig. 13 . Parameters of the mechanical and electrical submodels are shown in Table 1 .
Vehicle model parameters.
PPT Slide
Lager Image
Vehicle model parameters.
- 5.1. Model validation
5.1.1 Acceleration-Brake
This experiment is aimed to evaluate the electromechanical power interaction. The vehicle was conducted from rest to the maximum speed and then back to rest. This sequence was performed twice. Fig. 14 shows simulation and experimental results of IM1. Maximum 𝜔 1 ≈ 1050 rpm corresponds to Vx ≈ 30 km/h. In Fig. 15 the consumed electric power of IM1 is depicted. Its negative value during braking, evidences the regenerative braking capability. The correlation between experimental and simulation results ensures a correct prediction of the energy consumed for any driving cycle.
PPT Slide
Lager Image
Simulation/experimental results for Acceleration- Brake.
PPT Slide
Lager Image
Electric power during the Acceleration-Brake.
Along the complete maneuver the total energy demanded to drive the vehicle is ≈ 97.8 kJ whereas the recovered energy is ≈ 29.9 kJ. From these results it is possible to quantify the improvement of the energetic efficiency which is approximately 30.56%.
5.1.2 Turning Maneuver
Figs. 16 and 17 show simulation and experimental results of the steering angle and IMs variables obtained in this experiment.
PPT Slide
Lager Image
Steering angle evolution during the turning maneuver.
PPT Slide
Lager Image
Simulation and experimental results for the turning maneuver experiment: (a) Right and (b) left wheel motor.
The vehicle began stopped and then accelerated and remained at constant speed. Fig. 16 shows that the vehicle first turned right, then it turned left and came back to a straight direction. Fig. 17 shows that while the vehicle was turning right, the inner curve wheel speed (IM1) was smaller than that of the outer one (IM2) and then it changed as the vehicle direction changed. During the whole experiment, wheel torques (proportional to Iqs ) were the same for both motors as expected due to the implemented EDS algorithm.
The largest difference in Fig. 17 can be seen for Iqs (mean squared error of 3.96%). Simulation and experimental results during the turning maneuver presented EDS algorithm. The largest difference in Fig. 17 can be a mean squared error smaller than 4% for all variables. The inclusion of lateral dynamics allows reproducing turning maneuvers with high fidelity. This property is necessary to evaluate the vehicle maneuverability and stability under dangerous conditions. Furthermore, it helps to design vehicle yaw controllers and test their performance.
- 5.2. Potential applications
Results presented in this section emphasize the capability of the model to simulate dangerous situations and/or variables difficult or impossible to measure. Fig. 18 (a) shows an experiment where both left wheels were blocked during braking. The vehicle was driven at a constant speed of 70 km/h and after 0.5 seconds the brake was applied (mechanical brake on frontal wheels). At t = 0.6 s, road condition for left wheels was changed to extremely slippery, producing their blockage. Right wheels kept braking along producing a net torque. The vehicle continued to slow down as shown in Fig. 18(a) but the tangential velocities of the blocked and non-blocked wheels were absolutely different.
PPT Slide
Lager Image
(a) Velocities during blockage, (b) IMs energy flow.
Because of the wheel blockage, the regenerative braking capability of the IM was lost in the left rear wheel (continuous line in Fig. 18(b) ) whereas the IM located on right rear wheel (dashed line in Fig. 18(b) ) kept recovering energy.
An important yaw angle appeared ( Fig. 19(a) ) due to asymmetry on the longitudinal forces. This yaw rate
PPT Slide
Lager Image
produced load transference toward the left side and a negative roll angle occurred. Lateral slip is seen in Fig. 9(b) .
PPT Slide
Lager Image
(a) Roll and yaw evolution during blockage, (b) Lateral slip of rear and frontal wheels during blockage.
Variables of the FOC presented in Section Ⅳ-B are depicted in Fig. 20 . They are the measured Ia and the Park transformed Iqs in the IM connected to the blocked wheel.
PPT Slide
Lager Image
Currents behavior during blockage experiment.
Simulation results illustrate the wide range of applications possible to tackle. It is possible to assess the response during risky situations, power flow between electric and mechanical domains and variables regarding the control.
Moreover, the behavior of the tire-road interaction can be studied. These variables are needed to estimate and control the traction contact forces. In [28] the longitudinal dynamics of the proposed model was successfully used to this end.
6. Conclusion
A novel simulation methodology based on a combination of structural (for the physical system) and functional (for the controller) approaches oriented to developing electric vehicle (EV) concepts was proposed. A remarkable contribution of the proposal is that applying this combination it is possible to exploit the advantages of each approach.
The structural approach is based on Multi-Bond Graph (MBG) that allows representing each EV sub-system in a compact/modular way. This approach allows representing the power interactions between sub-systems independently of their physical nature. Moreover, the functional approach is accomplished using Block Diagrams (BD) that gives a clear overall description of the input/output relationships in the controller.
The proposed procedure consists on integrating both representations (represented in MBG and BD) into a single simulation environment. The resulting model allows evaluating and verifying rapidly any modification on the vehicle or in the controller, reducing time and development costs. The versatility of the proposal permits to implement other EV configurations and helps to understand the dynamics of electrically driven vehicles.
Experimental results were obtained using a real vehicle prototype. The high correlation between simulation and experimental results demonstrates the practical feasibility of the proposal.
BIO
L. I Silva He received his MSc in Space Sciences and Technologies from the Lulea Tekniska Universitet, Sweeden in 2007 and his PhD in Engineering Sciences from Universidad Nacional de Río Cuarto, Argentina in 2012. Since 2012 he is a postdoctoral student at Universidad Nacional de Río Cuarto. His research topics are oriented to modeling and simulation applied to electric and hybrid electric vehicles.
G. A Magallán He received the Electronic Engineer degree from the Universidad Tecnológica Nacional, Paraná, Argentina, in 2002 and the M.Sc. and Dr. Eng. degree from the Universidad Nacional de Río Cuarto, Río Cuarto, Argentina, in 2009 and 2010, respectively. He is also with the Consejo Nacional de Investigaciones Científicas y Técnicas, Buenos Aires, Argentina. His research interests include traction control of electric vehicles, dsp-based implementations, hybrid electric energy storage systems.
P. M. De la Barrera He received the B.Sc. and M.Sc. degrees in electrical engineering from Universidad Nacional de Rio Cuarto, Rio Cuarto, in 2003 and 2006, respectively, and the Ph.D. degree from Universidad Nacional del Sur, Bahía Blanca, Argentina, in 2009. Since 1998, he has been with the Grupo de Electrónica Aplicada, Universidad Nacional de Río Cuarto. He is also currently with Consejo Nacional de Investigaciones Científicas y Técnicas, Buenos Aires, Argentina. His research interests include modelling, control and fault diagnosis of electric machines, ac motor drives, electric vehicles, and renewable energy generation.
C. H. De Angelo He received the Electrical Engineer degree from the Universidad Nacional de Río Cuarto, Río Cuarto, Argentina, in 1999 and the Dr. Eng. degree from the Universidad Nacional de La Plata, La Plata, Argentina, in 2004. In 1994, he joined the Grupo de Electrónica Aplicada, Universidad Nacional de Río Cuarto. He is also with the Consejo Nacional de Investigaciones Científicas y Técnicas, Buenos Aires, Argentina. His research interests are fault diagnosis of electric machines, sensorless motor control, electric vehicles, and renewable-energy generation.
G. O. García He received the Electrical and Electronics Engineering degree from the Universidad Nacional de Córdoba, Córdoba, Argentina, in 1981 and the M.Sc. and Dr. degrees in electrical engineering from COPPE, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil, in 1990 and 1994, respectively. Since 1994, he has been with the Universidad Nacional de Río Cuarto, Río Cuarto, Argentina, where he is currently the Director of the Grupo de Electrónica Aplicada. He is also with the Consejo Nacional de Investigaciones Científicas y Técnicas, Buenos Aires, Argentina. His research interests include power electronics, electric vehicles, and renewable energy conversion.
References
Gao D. , Mi C. , Emadi A. 2007 “Modeling and simulation of electric and hybrid vehicles,” Proc. of the IEEE 95 (4) 729 - 745    DOI : 10.1109/JPROC.2006.890127
Kamil Çagatay Bayindir A. T. , Gözüküçük Mehmet Ali 2011 “A comprehensive overview of HEV: Powertrain configurations, powertrain control techniques and electronic control units,” Energy Conversion and Management 52 (2) 1305 - 1313    DOI : 10.1016/j.enconman.2010.09.028
Wang L. , Zhang Y. , Yin C. , Zhang H. , Wang C. 2012 “Hardware in the loop simulation for the design and verification of the control system of a series-parallel hybrid electric city-bus,” Simulation Modelling Practice and Theory 25 148 - 162    DOI : 10.1016/j.simpat.2012.02.010
Chan C. , Bouscayrol A. , Chen K. 2010 “Electric, hybrid, and fuel-cell vehicles: Architectures and modeling,” Vehicular Technology, IEEE Transactions on 59 (2) 589 - 598    DOI : 10.1109/TVT.2009.2033605
Gawthrop P. , Bevan G. 2007 “Bond-graph modeling,” Control Systems Magazine, IEEE 27 (2) 24 - 45
Karnopp D. C. , Margolis D. L. , Rosenberg R. C. 2006 System Dynamics: Modeling And Simulation of Mechatronic Systems Willey Intersciences New York, USA
Jang J. , Han C. 1998 “Proposition of a modeling method for constrained mechanical systems based on the vector bond graph,” Journal of the Franklin Institute 335 (3) 451 - 469    DOI : 10.1016/S0016-0032(96)00129-9
Schonfeld R. , Geitner G.-H. 2004 “Power flow and information flow in motion control system,” in Proc. EPE-PEMC
Geitner G.-H. 2006 “Power flow diagrams using a bond graph library under simulink,” in IEEE Industrial Electronics, IECON 2006-32nd Annual Conference on 5282 - 5288
Bouscayrol A. , Davat B. , de Fornel B. , Francois B. , Hautier J. P. , Meibody-Tabar F. , Pietrzak-David M. 2000 “Multimachine multiconverter system: application for electromechanical drives,” Eur. Phys. J., Appl. Phys. 10 (2) 131 - 147    DOI : 10.1051/epjap:2000124
Boulon L. , Hissel D. , Bouscayrol A. , Pape O. , Pera M. 2010 “Simulation model of a military HEV with a highly redundant architecture,” Vehicular Tech, IEEE Trans. on 59 (6) 2654 - 2663    DOI : 10.1109/TVT.2010.2045522
Wipke K. , Cuddy M. , Burch S. 1999 “ADVISOR 2.1: a user-friendly advanced powertrain simulation using a combined backward/forward approach,” Vehicular Technology, IEEE Transactions on 48 (6) 1751 - 1761    DOI : 10.1109/25.806767
Guzzella L. , Amstutz A. 1999 “CAE tools for quasistatic modeling and optimization of hybrid powertrains,” Vehicular Technology, IEEE Transactions on 48 (6) 1762 - 1769    DOI : 10.1109/25.806768
Somayajula D. , Meintz A. , Ferdowsi M. 2009 “Designing efficient hybrid electric vehicles,” Vehicular Technology Magazine, IEEE 4 (2) 65 - 72
Sekour M. , Hartani K. , Draou A. , Allali A. 2013 “Sensorless fuzzy direct torque control for high performance electric vehicle with four in-wheel motors,” Journal of Electrical Engineering and Technology 8 (3) 530 - 543    DOI : 10.5370/JEET.2013.8.3.530
Zupancic B. , Sodja A. 2012 “Computer-aided physical multi-domain modelling: Some experiences from education and industrial applications,” Simulation Modelling Practice and Theory
Borhan H. , Vahidi A. , Phillips A. , Kuang M. , Kolmanovsky I. , Di Cairano S. 2012 “MPC-based energy management of a powersplit hybrid electric vehicle,” Control Systems Technology, IEEE Trans-actions on 20 (3) 593 - 603    DOI : 10.1109/TCST.2011.2134852
Onoda S. , Emadi A. 2004 “PSIM-based modeling of automotive power systems: conventional, EVs, and HEVs,” Vehicular Technology, IEEE Transactions on 53 (2) 390 - 400    DOI : 10.1109/TVT.2004.823500
Wang R. , Wang J. 2011 “Fault-tolerant control with active fault diagnosis for 4W independently driven electric ground vehicles,” VT, IEEE Trans on 60 (9) 4276 - 4287
Delarue P. , Bouscayrol A. , Barrade P. 2010 “Energetic macroscopic representation and PSIM R_ simulation: Application to a DC / DC converter input filter stability,” in Vehicle Power and Propulsion Conference (VPPC), 2010 IEEE 1 - 6
Mihael Cipek J. P. , Pavkovic Danijel 2013 “A controloriented simulation model of a power-split hybrid electric vehicle,” Applied Energy 101 121 - 133    DOI : 10.1016/j.apenergy.2012.07.006
Bera T. , Bhattacharya K. , Samantaray A. 2011 “Evaluation of antilock braking system with an integrated model of full vehicle system dynamics,” Simulation Modelling Practice and Theory 19 (10) 2131 - 2150    DOI : 10.1016/j.simpat.2011.07.002
Silva L. , Magallan G. , De Angelo C. , Garcia G. 2008 “Vehicle dynamics using multi-bond graphs: Four wheel electric vehicle modeling,” in Industrial Electronics, 2008. IECON 2008. 34th Annual Conference of IEEE 2846 - 2851
Djeziri M. , Merzouki R. , Bouamama B. 2009 “Robust monitoring of an EV with structured and unstructured uncertainties,” Vehicular Tech, IEEE Trans. on 58 (9) 4710 - 4719    DOI : 10.1109/TVT.2009.2026281
2012 “Modelica Language Specification, version 3.3,” Available: https://www.modelica.org/documents
Zimmer D. , Cellier F. 2006 “The modelica multi-bond graph library,” in Proc. 5th Int. Modelica Conference 2 559 - 568
Velenis E. , Tsiotras P. , Canudas-de Wit C. , Sorine M. 2005 “Dynamic tyre friction models for combined longitudinal and lateral vehicle motion,” Vehicle System Dynamics 43 (1) 3 - 29
Magallan G. , De Angelo C. , Garcia G. 2011 “Maximization of the traction forces in a 2WD electric vehicle,” Vehicular Technology, IEEE Transactions on 60 (2) 369 - 380    DOI : 10.1109/TVT.2010.2091659
Kang J. , Yoo J. , Yi K. 2011 “Driving control algorithm for maneuverability, lateral stability, and rollover prevention of 4WD EV with independently driven front and rear wheels,” Vehicular Technology, IEEE Trans on 60 (7) 2987 - 3001    DOI : 10.1109/TVT.2011.2155105
Cho W. , Choi J. , Kim C. , Choi S. , Yi K. 2012 “Unified chassis control for the improvement of agility, maneuverability, and lateral stability,” Vehicular Technology, IEEE Transactions on 61 (3) 1008 - 1020    DOI : 10.1109/TVT.2012.2183152
Nam K. , Oh S. , Fujimoto H. , Hori Y. 2013 “Estimation of sideslip and roll angles of electric vehicles using lateral tire force sensors through rls and kalman filter approaches,” Industrial Electronics, IEEE Transactions on 60 (3) 988 - 1000    DOI : 10.1109/TIE.2012.2188874
Silva L. , Magallan G. , de la Barrera P. , De Angelo C. , Garcia G. 2010 “Modeling of EV dynamics with multi-bond graphs,” IEEE in Vehicle Power and Propulsion Conf. (VPPC’10) 1 - 7
Magallan G. A. , De Angelo C. H. , Garcia G. O. 2009 “A Neighborhood Electric Vehicle Development with Individual Traction on Rear Wheels,” Int. Journal of Electric and Hybrid Vehicles (IJEHV) 2 (2) 115 - 136    DOI : 10.1504/IJEHV.2009.029037
Sampaio R. , Hernandes A. , do Valle Magalha V. , Fernandes es , Becker M. , Siqueira A. 2012 “A new control architecture for robust controllers in rear electric traction passenger HEVs,” Vehicular Tech, IEEE Trans on 61 (8) 3441 - 3453    DOI : 10.1109/TVT.2012.2208486
Duma R. , Dobra P. , Abrudean M. , Dobra M. 2007 “Rapid prototyping of control systems using embedded target for TI c2000 DSP,” in Control & Automation. MED ’07 Conf. on 1 - 5