Co-simulation Methodologies for Hybrid and Electric Vehicle Dynamics

In recent decades, full electric and hybrid electric vehicles have emerged as an alternative to conventional cars due to a range of factors, including environmental and economic aspects. These vehicles are the result of considerable efforts to seek ways of reducing the use of fossil fuel for vehicle propulsion. Sophisticated technologies such as hybrid and electric powertrains require careful study and optimization. Mathematical models play a key role at this point. Currently, many advanced mathematical analysis tools, as well as computer applications have been built for vehicle simulation purposes. Given the great interest of hybrid and electric powertrains, along with the increasing importance of reliable computer-based models, the author decided to integrate both aspects in the research purpose of this work. Furthermore, this is one of the first final degree projects held at the ETSII (Higher Technical School of Industrial Engineers) that covers the study of hybrid and electric propulsion systems. The present project is based on MBS3D 2.0, a specialized software for the dynamic simulation of multibody systems developed at the UPM Institute of Automobile Research (INSIA). Automobiles are a clear example of complex multibody systems, which are present in nearly every field of engineering. The work presented here benefits from the availability of MBS3D software. This program has proven to be a very efficient tool, with a highly developed underlying mathematical formulation. On this basis, the focus of this project is the extension of MBS3D features in order to be able to perform dynamic simulations of hybrid and electric vehicle models. This requires the joint simulation of the mechanical model of the vehicle, together with the model of the hybrid or electric powertrain. These sub-models belong to completely different physical domains. In fact the powertrain consists of energy storage systems, electrical machines and power electronics, connected to purely mechanical components (wheels, suspension, transmission, clutch…). The challenge today is to create a global vehicle model that is valid for computer simulation. Therefore, the main goal of this project is to apply co-simulation methodologies to a comprehensive model of an electric vehicle, where sub-models from different areas of engineering are coupled. The created electric vehicle (EV) model consists of a separately excited DC electric motor, a Li-ion battery pack, a DC/DC chopper converter and a multibody vehicle model. Co-simulation techniques allow car designers to simulate complex vehicle architectures and behaviors, which are usually difficult to implement in a real environment due to safety and/or economic reasons. In addition, multi-domain computational models help to detect the effects of different driving patterns and parameters and improve the models in a fast and effective way. Automotive designers can greatly benefit from a multidisciplinary approach of new hybrid and electric vehicles. In this case, the global electric vehicle model includes an electrical subsystem and a mechanical subsystem. The electrical subsystem consists of three basic components: electric motor, battery pack and power converter. A modular representation is used for building the dynamic model of the vehicle drivetrain. This means that every component of the drivetrain (submodule) is modeled separately and has its own general dynamic model, with clearly defined inputs and outputs. Then, all the particular submodules are assembled according to the drivetrain configuration and, in this way, the power flow across the components is completely determined. Dynamic models of electrical components are often based on equivalent circuits, where Kirchhoff’s voltage and current laws are applied to draw the algebraic and differential equations. Here, Randles circuit is used for dynamic modeling of the battery and the electric motor is modeled through the analysis of the equivalent circuit of a separately excited DC motor, where the power converter is included. The mechanical subsystem is defined by MBS3D equations. These equations consider the position, velocity and acceleration of all the bodies comprising the vehicle multibody system. MBS3D 2.0 is entirely written in MATLAB and the structure of the program has been thoroughly studied and understood by the author. MBS3D software is adapted according to the requirements of the applied co-simulation method. Some of the core functions are modified, such as integrator and graphics, and several auxiliary functions are added in order to compute the mathematical model of the electrical components. By coupling and co-simulating both subsystems, it is possible to evaluate the dynamic interaction among all the components of the drivetrain. ‘Tight-coupling’ method is used to cosimulate the sub-models. This approach integrates all subsystems simultaneously and the results of the integration are exchanged by function-call. This means that the integration is done jointly for the mechanical and the electrical subsystem, under a single integrator and then, the speed of integration is determined by the slower subsystem. Simulations are then used to show the performance of the developed EV model. However, this project focuses more on the validation of the computational and mathematical tool for electric and hybrid vehicle simulation. For this purpose, a detailed study and comparison of different integrators within the MATLAB environment is done. Consequently, the main efforts are directed towards the implementation of co-simulation techniques in MBS3D software. In this regard, it is not intended to create an extremely precise EV model in terms of real vehicle performance, although an acceptable level of accuracy is achieved. The gap between the EV model and the real system is filled, in a way, by introducing the gas and brake pedals input, which reflects the actual driver behavior. This input is included directly in the differential equations of the model, and determines the amount of current provided to the electric motor. For a separately excited DC motor, the rotor current is proportional to the traction torque delivered to the car wheels. Therefore, as it occurs in the case of real vehicle models, the propulsion torque in the mathematical model is controlled through acceleration and brake pedal commands. The designed transmission system also includes a reduction gear that adapts the torque coming for the motor drive and transfers it. The main contribution of this project is, therefore, the implementation of a new calculation path for the wheel torques, based on performance characteristics and outputs of the electric powertrain model. Originally, the wheel traction and braking torques were input to MBS3D through a vector directly computed by the user in a MATLAB script. Now, they are calculated as a function of the motor current which, in turn, depends on the current provided by the battery pack across the DC/DC chopper converter. The motor and battery currents and voltages are the solutions of the electrical ODE (Ordinary Differential Equation) system coupled to the multibody system. Simultaneously, the outputs of MBS3D model are the position, velocity and acceleration of the vehicle at all times. The motor shaft speed is computed from the output vehicle speed considering the wheel radius, the gear reduction ratio and the transmission efficiency. This motor shaft speed, somehow available from MBS3D model, is then introduced in the differential equations corresponding to the electrical subsystem. In this way, MBS3D and the electrical powertrain model are interconnected and both subsystems exchange values resulting as expected with tight-coupling approach.When programming mathematical models of complex systems, code optimization is a key step in the process. A way to improve the overall performance of the integration, making use of C/C++ as an alternative programming language, is described and implemented. Although this entails a higher computational burden, it leads to important advantages regarding cosimulation speed and stability. In order to do this, it is necessary to integrate MATLAB with another integrated development environment (IDE), where C/C++ code can be generated and executed. In this project, C/C++ files are programmed in Microsoft Visual Studio and the interface between both IDEs is created by building C/C++ MEX file functions. These programs contain functions or subroutines that can be dynamically linked and executed from MATLAB. This process achieves reductions in simulation time up to two orders of magnitude. The tests performed with different integrators, also reveal the stiff character of the differential equations corresponding to the electrical subsystem, and allow the improvement of the cosimulation process. When varying the parameters of the integration and/or the initial conditions of the problem, the solutions of the system of equations show better dynamic response and stability, depending on the integrator used. Several integrators, with variable and non-variable step-size, and for stiff and non-stiff problems are applied to the coupled ODE system. Then, the results are analyzed, compared and discussed. From all the above, the project can be divided into four main parts: 1. Creation of the equation-based electric vehicle model; 2. Programming, simulation and adjustment of the electric vehicle model; 3. Application of co-simulation methodologies to MBS3D and the electric powertrain subsystem; and 4. Code optimization and study of different integrators. Additionally, in order to deeply understand the context of the project, the first chapters include an introduction to basic vehicle dynamics, current classification of hybrid and electric vehicles and an explanation of the involved technologies such as brake energy regeneration, electric and non-electric propulsion systems for EVs and HEVs (hybrid electric vehicles) and their control strategies. Later, the problem of dynamic modeling of hybrid and electric vehicles is discussed. The integrated development environment and the simulation tool are also briefly described. The core chapters include an explanation of the major co-simulation methodologies and how they have been programmed and applied to the electric powertrain model together with the multibody system dynamic model. Finally, the last chapters summarize the main results and conclusions of the project and propose further research topics. In conclusion, co-simulation methodologies are applicable within the integrated development environments MATLAB and Visual Studio, and the simulation tool MBS3D 2.0, where equation-based models of multidisciplinary subsystems, consisting of mechanical and electrical components, are coupled and integrated in a very efficient way.

Electric current of an electrified jet issuing from a long metallic tube.

This paper presents an analysis of the transport of electric current in a jet of an electrically conducting liquid discharging from a metallic tube into a gas or a vacuum, and subject to an electric field due to a high voltage applied between the tube and a far electrode. The flow, the surface charge and the electric field are computed in the current transfer region of the jet, where conduction current in the liquid becomes surface current due to the convection of electric charge accumulated at its surface. The electric current computed as a function of the flow rate of the liquid injected through the tube increases first as the square root of this flow rate, levels to a nearly constant value when the flow rate is increased and finally sets to a linear increase when the flow rate is further increased. The current increases linearly with the applied voltage at small and moderate values of this variable, and faster than linearly at high voltages. The characteristic length and structure of the current transfer region are determined. Order-of-magnitude estimates for jets which are only weakly stretched by the electric stresses are worked out that qualitatively account for some of the numerical results.

Development of a grid connected micro wind generator. A practical activity for the course on electric generation with wind energy

This paper describes a practical activity, part of the renewable energy course where the students have to build their own complete wind generation system, including blades, PM-generator, power electronics and control. After connecting the system to the electric grid the system has been tested during real wind scenarios. The paper will describe the electric part of the work surface-mounted permanent magnet machine design criteria as well as the power electronics part for the power control and the grid connection. A Kalman filter is used for the voltage phase estimation and current commands obtained in order to control active and reactive power. The connection to the grid has been done and active and reactive power has been measured in the system.

Organisational learning, strategic rigidity and technology adoption: Implications for electric utilities and renewable energy firms

This paper examines the implications of strategic rigidness for technology adoption behaviours among electric utilities. Such behaviours lead to heterogeneity in firm performance and consequently affect the electric utility industry. The paper's central aim is to identify and describe the implications of strategic rigidness for a utility firm's decision making in adopting newer renewable energy technologies. The findings indicate that not all utility firms are keen to adopt these new technologies, as these firms have traditionally been operating efficiently with a more conventional and mature technological arrangement that has become embedded in the organisational routine. Case studies of Iberdrola S.A. and Enel S.p.A. as major electric utilities are detailed to document mergers and acquisitions and technology adoption decisions. The results indicate that technology adoption behaviours vary widely across utility firms with different organisational learning processes and core capabilities.

Ion correlations due to high-frequency electric field and their effect on the nonlinear plasma conductivity

The influence of a strong, high‐frequency electric field on the ion‐ion correlations in a fully ionized plasma is investigated in the limit of infinite ion mass, starting with the Bogoliubov‐Born‐Green‐Kirkwood‐Yvon hierarchy of equations; a significant departure from the thermal correlations is found. It is shown that the above effect may substantially modify earlier results on the nonlinear high‐frequency plasma conductivity.

Electrostatic plasma instabilities excited by a high‐frequency electric field

The electrostatic plasma waves excited by a uniform, alternating electric field of arbitrary intensity are studied on the basis of the Vlasov equation; their dispersion relation, which involves the determinant of either of two infinite matrices, is derived. For ω0 ≫ ωpi (ω0 being the applied frequency and ωpi the ion plasma frequency) the waves may be classified in two groups, each satisfying a simple condition; this allows writing the dispersion relation in closed form. Both groups coalesce (resonance) if (a) ω0  ≈  ωpe/r (r any integer) and (b) the wavenumber k is small. A nonoscillatory instability is found; its distinction from the DuBois‐Goldman instability and its physical origin are discussed. Conditions for its excitation (in particular, upper limits to ω0,k, and k⋅vE,vE being the field‐induced electron velocity), and simple equations for the growth rate are given off‐resonance and at ω0  ≈  ωpi. The dependence of both threshold and maximum growth rate on various parameters is discussed, and the results are compared with those of Silin and Nishikawa. The threshold at ω0  ≈  ωpi/r,r  ≠  1, is studied.

Finite Element Analysis Model of a Contactless Transformer for Battery Chargers in Electric Vehicles

A contactless transformer model is proposed in this paper using Finite Element Analysis (FEA). This model can be used to simulate Inductive Coupling Power Transfer (ICPT) systems with good accuracy of the transformer and reduce the fabrication time of these systems. The model not only takes into account the geometry of the windings but also the frequency effects in them. As the transformer does not have a magnetic core, it is complicated to model because the flux is expanded in the area around the windings. In order to obtain a very accurate model, it is necessary to use a 2D/3D field solver.

Magnetic and electric current morphoiogy in the plasma ejected by a laser-irradiated foil

A quasisteady model for the plasma ablated from a thick foil by a laser pulse, at low $lln $ and R /A i within a low, narrow range, is given (4, is absorbed intensity, /zL wavelength, R focalspot radius). An approximate analytical solution is given for the two-dimensional plasma dynamics. At large magnetic Reynolds number Rm, the morphology of the magnetic field shows features in agreement with recent results for high intensities. Current lines are open: electric current flows toward the spot near its axis, then turns and flows away. The efficiency of converting light energy into electric energy peaks at Rm- 1, both the validity of the model. and accuracy of the solution are discussed, The neighborhood of the spot boundary is analyzed in detail by extending classical Prandtl-Meyer results.

Electrostatic plasma instabilities excited by a high-frequency electric field. AT(30-1)-1238-MATT-669

An analysis of the electrostatic plasma instabilities excited by the application of a strong, uniform, alternating electric field is made on the basis of the Vlasov equation. A very general dispersion relation is obtained and discussed. Under the assumption W 2 O » C 2 pi. (where wO is the applied frequency and wpi the ion plasma frequency) a detailed analysis is given for wavelengths of the order of or large compared with the Debye length. It is found that there are two types of instabilities: resonant (or parametric) and nonresonant. The second is caused by the relative streaming of ions and electrons, generated by the field; it seems to exist only if wO is less than the electron plasma frequency wpe. The instability only appears if the field exceeds a certain threshold, which is found.

Optical behaviour of cholesteric blue phases in electric fields

The optical behaviour of cholesteric mixtures of negative dielectric anisotrony under electric fields is reported. A mixture of S 311~ (31.35 %) + N 5 was employed. AC voltages (f = 1000 Hz) betweeen 0 and 150 volts were applied. Cells 23 micron thick, with internal SnO2 electrodes, were used.

Strategic rigidity and foresight for technology adoption among electric utilities

The variation in the adoption of a technology as a major source of competitive advantage has been attributed to the wide-ranging strategic foresight and the integrative capability of a firm. These possible areas of competitive advantage can exist in the periphery of the firm's strategic vision and can get easily blurred as a result of rigidness and can permeate in the decision-making process of the firm. This article explores how electric utility firms with a renewable energy portfolio can become strategically rigid in terms of adoption of newer technologies. The reluctance or delay in the adoption of new technology can be characterized as strategic rigidness, brought upon as a result of a firm's core competence or core capability in the other, more conventional technology arrangement. This paper explores the implications of such rigidness on the performance of a firm and consequently on the energy eco-system. The paper substantiates the results by emphasizing the case of Iberdrola S.A., an incumbent firm as a wind energy developer and its adoption decision behavior. We illustrate that the very routines that create competitive advantage for firms in the electric utility industry are vulnerable as they might also develop as sources of competitive disadvantage, when firms confront environmental change and uncertainty.

Evaluation of the Magnetic Field Generated by the Inverter of an Electric Vehicle

In hybrid and electric vehicles, passengers sit very close to an electric system of significant power, which means that they may be subjected to high electromagnetic fields. The hazards of long-term exposure to these fields must be taken into account when designing electric vehicles and their components. Among all the electric devices present in the power train, the electronic converter is the most difficult to analyze, given that it works with different frequencies. In this paper, a methodology to evaluate the magnetic field created by a power electronics converter is proposed. After a brief overview of the recommendations of electromagnetic fields exposure, the magnetic field produced by an inverter is analyzed using finite element techniques. The results obtained are compared to laboratory measurements, taken from a real inverter, in order to validate the model. Finally, results are used to draw some conclusions regarding vehicle design criteria and magnetic shielding efficiency.

Dual Active Bridge based Battery Charger for Plug-in Hybrid Electric Vehicle with Charging Current Containing Low Frequency Ripple

High power density is strongly preferable for the on-board battery charger of Plug-in Hybrid Electric Vehicle (PHEV). Wide band gap devices, such as Gallium Nitride HEMTs are being explored to push to higher switching frequency and reduce passive component size. In this case, the bulk DC link capacitor of AC-DC Power Factor Correction (PFC) stage, which is usually necessary to store ripple power of two times the line frequency in a DC current charging system, becomes a major barrier on power density. If low frequency ripple is allowed in the battery, the DC link capacitance can be significantly reduced. This paper focuses on the operation of a battery charging system, which is comprised of one Full Bridge (FB) AC-DC stage and one Dual Active Bridge (DAB) DC-DC stage, with charging current containing low frequency ripple at two times line frequency, designated as sinusoidal charging. DAB operation under sinusoidal charging is investigated. Two types of control schemes are proposed and implemented in an experimental prototype. It is proved that closed loop current control is the better. Full system test including both FB AC-DC stage and DAB DC-DC stage verified the concept of sinusoidal charging, which may lead to potentially very high power density battery charger for PHEV.

Nationwide impact and vehicle to grid application of electric vehicles mobility using an activity based model

This paper describes the impact of electric mobility on the transmission grid in Flanders region (Belgium), using a micro-simulation activity based models. These models are used to provide temporal and spatial estimation of energy and power demanded by electric vehicles (EVs) in different mobility zones. The increment in the load demand due to electric mobility is added to the background load demand in these mobility areas and the effects over the transmission substations are analyzed. From this information, the total storage capacity per zone is evaluated and some strategies for EV aggregator are proposed, allowing the aggregator to fulfill bids on the electricity markets.

Flexible low-cost system to test batteries and ultracapacitors for electric and hybrid vehicles in real working conditions

Batteries and ultracapacitors for hybrid and electric vehicles must satisfy very demanding working conditions that are not usual in other applications. In this sense, specific tests must be performed in order to draw accurate conclusions about their behaviour. To do so, new advanced test benches are needed. These platforms must allow the study of a wide variety of energy storage systems under conditions similar to the real ones. In this paper, a flexible, low-cost and highly customizable system is presented. This system allows batteries and ultracapacitors to be tested in many and varied ways, effectively emulating the working conditions that they face in an electric vehicle. The platform was specifically designed to study energy storage systems for electric and hybrid vehicles, meaning that it is suitable to test different systems in many different working conditions, including real driving cycles. This flexibility is achieved keeping the cost of the platform low, which makes the proposed test bench a feasible alternative for the industry. As an example of the functionality of the platform, a test consisting of a 17-minute ARTEMIS urban cycle with a NiMH battery pack is presented.