Modeling of Battery for EV using EMTP/ATPDraw

  • Kim, Jun-Hyeok (Dept. of Electrical and Electronic Engineering, Sungkyunkwan University) ;
  • Lee, Soon-Jeong (Dept. of Electrical and Electronic Engineering, Sungkyunkwan University) ;
  • Kim, Eung-Sang (Korea Electrotechnology Research Institute) ;
  • Kim, Seul-Ki (Korea Electrotechnology Research Institute) ;
  • Kim, Chul-Hwan (Dept. of Electrical and Electronic Engineering, Sungkyunkwan University) ;
  • Prikler, Laszlo (Dept. Electric Power Engineering, Budapest University of Technology and Economics)
  • Received : 2013.03.19
  • Accepted : 2013.09.26
  • Published : 2014.01.01


As environmentally friendly energy takes center stage, interests for Electric Vehicles/Plug in Hybrid Electric Vehicles (EVs/PHEVs) are getting increase. With this trend, there is no doubt EVs will take large portion to penetrations of total cars. Therefore, accurate EV modeling is required. Battery is one of the main components with the power system view of aspect. Hence, in this paper, reviews and discussions of some types of batteries for EV are contained by considering energy density and weight of the batteries. In addition, simulations of Li-Ion battery are accomplished with various variables such as temperature, capacity fading and charge/discharge current. It is confirmed that temperature is the main factor of capacity fading. Validation of the modeled battery is also conducted by comparing it with commercialized battery.

1. Introduction

These days most of countries and companies focus on eco-friendly growth. With this circumstance, Electric Vehicles (EVs) take center stage. EVs can reduce the amount of consumption of fuel fossil, and thus, it can reduce the amount of CO2 emission. However, this is not the only reason why the interests about EVs are growing. EVs can operate as distributed generation as well. That is to say, it can operate as power system device for compensating insufficient power or for SMART GRID. For these reasons, therefore, it is clear that EVs will take large portion of total complement of car. In this state of affairs, analyses of effect on distribution system with EVs are essential. And accurate battery modeling should be considered as first priority for the analyses.

Since future studies of this paper are analyzing the effects on grid when EVs are connected, EMTP/ATPDraw which is one of the widely used programs for analyzing transient phenomena is used. However, there is no such battery model in that program. Hence modeling of battery with EMTP/ATPDraw is conducted in this paper.

There are some types of battery modeling, specifically: experimental, circuit-based and mathematical. Experimental modeling and circuit-based modeling cannot estimate State-of-Charge (SoC). However, estimations of SoC are very important part of battery modeling. Hence, mathematical modeling is used in this paper. This paper is divided into four main parts. In the section 2, brief introductions of batteries are contained. And battery for modeling in this paper is selected based on section 2. In the section 3, theoretical backgrounds of battery modeling are presented. The section 4 contains simulation result of modeled battery with various aspects such as internal resistance, charging/discharging current and capacity fading. And it is verified in section 5. This paper ends with conclusion in the last section.


2. Battery Selection

There are some types of batteries for EV such as Lead-Acid, Lithium-Ion, Nickel-Cadmium and Nickel-Metal-Hydride. To model battery, selection of battery is one of the most crucial factors that have to be considered. Brief characteristics of each type of batteries are indicated in Table 1 [1-5].

Table 1.Brief characteristics of batteries for EV

However the first priority considering EV is how long EV can drive with the battery. Therefore, weight and energy density of the battery should be considered first. Brief data of weights and costs of batteries are summarized in Table 2 [6].

Table 2.Weights and costs of batteries those are applicable to electric vehicles

The more energy the battery has, the more driving distance EV gets. And in the same way, the less weight the battery has, the more driving distance EV gets. Hence, the battery which has high energy density and low weight should be contemplated to use for EV. In this respect, Li-Ion batteries are getting the limelight because of its high energy density and relatively low weight as shown in Table 2. Therefore, Li-Ion battery is used in this paper.


3. Battery Modeling

3.1 EMTP/ATPDraw

In this paper EMTP/ATPDraw is used for modeling the battery. The Electromagnetic Transient Program (EMTP) is the tool used to simulate transient electromagnetic (EM) phenomena, and it is one of the most widely used programs throughout electric utilities [7]. And ATPDraw is a graphical, mouse-driven pre-processor to the ATP version of the Electromagnetic Transients Program (EMTP) [8].

MODELS in ATP is a general-purpose description language supported by an extensive set of simulation tools for the representation and study of time-variant systems [8]. The MODELS provides the monitoring and controllability of power system as well as some other algebraic and relational operations for programming. With some compromised functions such as repetition, conditional path selection, and user defined functions, it is also called a new TACS [9].

As stated, there is no doubt that EV will be brought out. However, without analyses for possible effect on power system, it cannot be commercialized. In this paper, EMTP/ ATPDraw which has advantages of analyzing transient phenomenon is used. As mentioned, however, there is no battery model with the program. Thus, modeling of battery with the program is indispensable not only analyzing battery itself but also analyzing the effects on grid when EVs are connected.

3.2 Battery models

There are some types of battery modeling such as experimental model, circuit-based model and mathematical model. However, experimental model does not properly describe cell dynamics. Hence, circuit-based model and mathematical model are discussed in this paper [10].

3.2.1 Circuit-based model

Since this method can represent electrical characteristics of battery, circuit-based model can possibly be considered as suitable modeling method. The most commonly used circuit-based battery model is Thevenin battery model as in Fig. 1.

Fig. 1.Thevenin battery model

Where,Voc Open-circuit voltage (V)R0 Internal resistance (Ω)R1 Overvoltage resistance (Ω)C1 Capacitance of the battery (F)Vbatt Battery voltage (V)

The main disadvantage of this model is all components in this equivalent circuit are assumed to be constant. However, it varies depending on the conditions of battery [11]. Furthermore, circuit-based models including thevenin model do not take account of SoC.

3.2.2 Mathematical model

In 1965, Shepherd [12] suggested battery model mathematically, as in (1). It describes electrochemical behavior of the battery directly [13].

Where,E0 No-load voltage (V)K Polarization voltage (V)Q Battery capacity (Ah)A Exponential zone amplitude (V)B Exponential zone time constant inverse (Ah)-1R Internal resistance (Ω)it Actual battery charge (Ah)i Battery current (A)Vbatt Output voltage of the battery (V)

However, this model cannot reflect performance of the battery because it assumes internal resistance as constant. As mentioned, it is impossible to describe exact performance of the battery with fixed internal components.

Generic battery model [14] progresses with Shepherd’s model, as in (2-3).

In Generic battery model, K(Q/(Q−it)i term in Shepherd’s model is replaced with K(Q/(Q−it) to prevent algebraic loop and simulation instability. Nonetheless, this model still assumed internal resistances as constant.

In reality, however, internal resistances are varied depending on the value of SoC. On the other hand, in [15] all equations are functions of SoC. That makes possible internal resistances have various values according to the value of SoC. Therefore, equations below [15] are used, in this paper.

Where,VOC Open-Circuit Voltage (V)RTransient_S, CTransient_S Short-time constantsRTransient_L, CTransient_L Long-time constantsRSeries Immediate voltage drop constant

In above Eqs. (4-9), RSeries is responsible for immediate voltage drop of the step response. RTransient_S and CTransient_S are responsible for short-time constant of the step response. RTransient_L and CTransient_L are responsible for long-time constant of the step response [15]. The main advantage of using these equations is the fact that internal resistances can be varied depending on the value of SoC.

3.3 Measurement methods for state-of-charge

There are three main methods to measure the SoC: voltage measurement method, impedance measurement method and coulomb counting method [16].

Voltage Measurement Method (VMM) is the simplest method to calculate SoC. It uses a relationship between output voltage of the battery and SoC. However, the relationship can be affected by temperature and charge/ discharge rate of the battery. Hence, although it is easy to calculate SoC, the result is rather roughly accurate.

Impedance Measurement Method (IMM) uses relationship between internal impedance and SoC. As stated above, internal impedance is varied depending on conditions of the battery. And SoC is one of those conditions. Thus, if the value of internal impedance is identified, it is easy to calculate the value of SoC. However, it is not widely used because of the difficulties to measure internal resistances while the battery is charged or discharged.

Current is used in Coulomb Counting Method (CCM). TCurrent is used in Coulomb Counting Method (CCM). The unit of the energy in the electric charge is coulombs and that is equal to the integration of overtime-current. Therefore, from this point, the remaining capacity in the battery can be calculated by measuring the current which is flowing into or leaving from the battery. To calculate accurate SoC in CCM, self-discharge current and columbic efficiency should be taken into account. But still it is useful and is making acceptably accurate result. So, in this paper, CCM is used with ignoring self-discharge current and columbic efficiency. Calculation of SoC using CCM is following (10)~(11).

Where,Q Battery capacity (Ah)i Battery current (A)SoC0 Initial value of SoC

3.4 Capacity fading

“Capacity fading” is one of the main characteristics of battery. Capacity fading refers to the phenomenon that usable capacity of the battery decreases due to run-time, temperature and charging/discharging cycle. Generally, it is considered that the battery is available until it has 80% of its initial capacity. Hence, with this, consideration for the effect of capacity fading is important factor in battery modeling [17].

There are three main reasons causing capacity fading which are temperature, time and cycle of the battery. In (12), capacity loss by temperature and time is expressed [18].

Where,Q Capacity loss (%)A Pre-exponential factorEa Activation energy (J)R Gas constantT Temperature (K)t Timez Exponent of time

Among above parameters, A and Ea are adjustable parameters. In this paper, they were set at 1.1443×106, 4.257×104 and 0.5 respectively [18, 19].

Capacity fading caused by cycling is indicated in (13) [20, 21].

Where,Qloss_cycle Capacity loss by cycling (%)IPR Number of cycle converted by Ah-process

In (13), capacity fading due to cycle is expressed as similar as (12). The number of cycles is converted to Ah-processed (IPR), the integral of absolute value of current over time [21]. Ah-variable represents the amount of charge transported during cycling. Advantage of using Ahprocess is allowances for quantifying and correlating the capacity fading behaviors with different C-rates [20].


4. Modeling of the Battery using EMTP/ATPDRAW

4.1 State-of-charge

As stated, CCM is used in this paper. The equations for calculation are following (10)~(11).

In this paper, battery capacity Q is set at 2.3Ah. And it is assumed that the modeled battery has negative current when it is charged, on the other hand, has positive current when it is discharged. According to the equations, SoC should have linear characteristic if the battery current set as constant. Therefore, as the battery current increase, the time which SoC is reached to 1 will be shorten.

4.2 Internal resistances

As SoC is changed, internal resistances are changed as well. Repeatedly, this characteristic is one of the most important factors. Simulation result of the internal resistances is presented in Table 3 and Fig. 2.

Table 3.Battery internal resistances

Fig. 2.Simulation result of internal resistances

As shown in Table 3 and Fig. 2, using (5), (6) and (8), internal resistances can be calculated and it rapidly fluctuates with SoC while SoC increases from 0 to 0.1. And it is gradually changed until it reaches each of their minimum values. This implies most of power losses are occurred when charging or discharging of the battery begins or ends, respectively.

4.3 Charging characteristics

For charging mode, it is assumed that No-load is connected. This means, in this simulation, only Voc is concerned. Simulation is conducted using (4) and simulation result is shown in Fig. 3.

Fig. 3.Simulation result of Voc when the battery is charged

In this simulation, battery capacity is set at 2.3Ah and charging current is set at 1A as constant value.

Generally, Open-Circuit Voltage Voc sharply changes when charging is started. In other words, when SoC is roughly between 0 to 0.2, Voc rapidly changes. And it gradually changes as SoC is reached to 1. This tendency is indicated in Fig. 3. It is noticed Voc starts at 2.6V and reaches at 4.1V in Fig. 3. That is to say, cut-off voltage and nominal voltage of modeled battery is 2.6V and 4.1 V, respectively.

4.4 Discharging characteristics

4.4.1 Effect of discharging current

Discharging mode is also following (4-9). The main difference between charging mode and discharging mode is consideration of internal resistances. In discharging mode, it should be taken into account. Theoretically, as discharge current increase, power losses due to internal resistances should be larger. The reason is power loss follows “P = I2R”. Simulation results when the discharging current varies are presented in Fig. 4.

Fig. 4.Simulation result of output voltage of the modeled battery when discharge current caries 0.1C to 1C

In Fig. 4, C means the discharge current which can make battery fully discharged in 1 hour. As expected, power consumed by internal resistances is grown as discharge current increase.

4.4.2 Capacity fading

Effect of capacity fading can be analyzed by (12) which is based on Arrhenius equation. In this simulation, irreversible capacity loss according to time is analyzed with three different temperatures.

Theoretically, as temperature higher, chemical reactions in battery are become faster. For this reason, the performance of battery can provide better performance, however, at the same time, it causes extra power loss. In other words, high temperature causes reduction of the battery life. As a rule of thumb, reaction rate of the battery is become doubles as the temperature of battery increase 10 ℃. For instance, if specific battery can be used in an hour at 30℃, it can be used in two hours at 20℃.

Time is also one of the main factors causing capacity fading. As mentioned previously, battery makes output voltage using chemical reaction. However, although users do not use the battery, it does not mean its chemical reactions stop either. Thus, as time goes by, the battery life is kept reducing.

For the reasons above, it is expected that life reduction or capacity loss of the battery becomes larger as temperature and time increase.

As shown in Table 4 and Fig. 5, power loss due to time and temperature gets worse as the time and temperature increase. In addition, it is noticed that higher temperature causes severely more power loss in much shorter time than those of lower temperatures. At 15000 days, for example, 1.41 percent of power losses are occurred at 5℃, however 4.88 and 14.36 percent of power losses are occurred at 25 ℃ and 45 ℃, respectively. This simulation result verifies theoretical expectation, as stated above.

Table 4.Power loss due to time and temperature

Fig. 5.Simulation results of capacity loss versus time

As charging/discharging repeated, irreversible capacity loss is occurred as well and it is called “cycling fading”. Simulation results of capacity loss due to cycling are presented in Table 5 and Fig. 6.

Table 5.Power loss due to cycle number and temperature

Fig. 6.Simulation results of capacity loss due to charging/discharging cycle (conducted with 1C)

In this simulation, 1C is used and it is shown that capacity losses by cycling are significantly worse when it is conducted in high temperature. Percentage capacity loss is 1.537% at 5℃. On the other hand, it causes 15.552% of capacity loss at 45℃.

Capacity fading by temperature, time and cycle, by the way, cannot be separated. Therefore, it should be considered at once. Simulation results considering those three variables are shown in Fig. 7.

Fig. 7.Simulation results of total capacity loss (25℃ and 45℃)

This simulation is conducted with 25℃ and 45℃, and assumed charging/discharging cycle is implemented one time per day. As shown in Fig. 7, it causes about 14 and 41 percent of power loss at 25℃ and 45℃, respectively. It suggests that temperature is the most crucial factor affect power loss of the battery.


5. Validation of Modeled Battery

Comparisons between modeled battery and commercialized battery “ANR26650” are presented in Fig. 8 ~ Fig. 10. In these comparisons, capacity of modeled battery Q is set at 2.3Ah and discharge current is provided sequentially.

Fig. 8.Comparison of battery output voltage in modeled battery and commercialized battery at 25℃

Fig. 9.Comparison of capacity loss by cycling in modeled battery and commercialized battery at 25℃

Fig. 10.Comparison of battery output voltage in modeled battery and commercialized battery at 45℃

Owing to the difference in open-circuit voltage, modelled battery and ANR26650 have not exactly same result, however, those two batteries have considerably similar pattern, as shown in Fig. 8. And analogousness of capacity loss by cycling is also shown in Fig. 9 and Fig. 10.


6. Conclusion

Since the battery for EV should have high energy density and low weight, Li-Ion battery is chosen for this paper. And mathematical battery model is used because of some explained disadvantages of circuit-based model and experimental model. It is introduced that the three main factors for battery modeling are SoC, internal resistances and capacity fading. With mathematical modeling, SoC and varying internal resistances are taken into account. And it is shown that capacity fading can be caused by temperature, time and cycling. The higher temperature batteries have the more power losses it gets. Likewise, the longer time and larger cycle numbers cause more power losses on battery. Therefore, it can lead reduction of battery life as presented in simulated results. Especially, it is shown that temperature is the most crucial factor among stated variables.

Validation of modeled battery is presented by comparing commercialized battery (ANR26650) and modeled battery. Simulation results verify its validity.

Without sufficient researches of possible effect caused by connected EVs in power system, it is impossible to employ EV technologies. Especially, transient phenomena on power system are more likely to be affected by EVs. Therefore, for analyzing battery characteristics and the effects on grid when EVs are connected are mainly discussed. To be specific, EMTP/ATPDraw which is one of the most appropriate programs for analyzing transient phenomena and simulations on power systems is used for that.

However, despite its advantages in analyzing power system such as transient phenomena, EMTP/ATPDraw does not have any battery models. Hence, modeling battery using EMTP/ATPDraw is worth conducting. So with verified battery in this paper, it is expected that the effects on grid when EVs are connected will be researched for further studies


Supported by : Korea Institute of Energy Technology Evaluation and Planning (KETEP)


  1. Isidor Buchmann, "Batteries in a Portable World", 2001.
  2. D. A. Corrigan, "Introduction to NiMH Battery Technology," June 21, 2002.
  3. Dr. Forrest A.Trumbore ed, "Modern Battery Technology" Center for Professional Advancement, 1995.
  4. M. A. Fetcenko, S. Venkatesan, S. R. Ovshinsky, "Selection of Metal Hydride Alloys For Electrochemical Applications," Electrochemical SoCiety Proceedings Volume 92-5, pp. 141-167Dr. Forrest A.Trumbore ed, "Modern Battery Technology" Center for Professional Advancement, 1995.
  5. R.F. Nelson, "High-power batteries for the new 36/42 V automotive systems", Journal of Power Sources 107 (2002) 226-239
  6. Henry Oman, "On-board energy and power management on electric vehicles: effect of battery type", The AIAA/IEEE/SAE Digital Avionics Systems Conference, November 1998.
  7. Chul-Hwan Kim, Myung-Hee Lee, Raj K. Aggarwal, Allan T. Johns, "Educational Use of EMTP MODELS for the Study of a Distance Relaying Algorithm for Protecting Transmission Lines", IEEE Transactions on Power Systems, Vol. 15, No. 1, February 2000.
  8. ATPDraw for Windows 3.1x/95/NT User's Manual.
  9. Keon-Woo Park, Hun-Chul Seo, Chul-Hwan Kim, Chang-Soo Jung, Yeon-Pyo Yoo, Yong-Hoon Lim, "Analysis of the Neutral Current for Two-Step-Type Poles in Distribution Lines", IEEE Transactions on power delivery, Vol. 24, No. 3, July 2009.
  10. Olivier Tremblay, Louis-A. Dessaint, "Experimental Validation of a Battery Dynamic Model for EV Applications", World Electric Vehicle Journal, Vol. 3, May 2009.
  11. H. L. Chan, D. Sutanto, "A New Battery Model for use with Battery Energy Storage Systems and Electric Vehicles Power Systems", Power Engineering SoCiety Winter Meeting, 2000.
  12. C. M. Shepherd, "Design of Primary and Secondary Cells" Journal of The Electrochemical SoCiety, Vol. 112, pp. 657-664, 1965.
  13. Hanlei Zhang and Mo-Yuen Chow, "Comprehensive Dynamic Battery Modeling for PHEV Applications", Power and Energy SoCiety General Meeting, 2010.
  14. Olivier Tremblay, Louis-A. Dessaint, Abdel-Illah Dekkiche, "A Generic Battery Model for the Dynamic Simulation of Hybrid Electric Vehicles", Vehicle Power and Propulsion Conference, 2007.
  15. Min Chen, Gabriel A. Rincon-Mora, "Accurate Electrical Battery Model Capable of Predicting Runtime and I-V Performance", IEEE Transactions on Energy Conversion, June 2006.
  16. Vipin Prajapati, Herbert Hess, Edward James William II, Vishu Gupta, Mattew Huff, Milos Manic, Freeman Rufus, Ash Thakker, Justin Govar, "A Literature Review of State of Charge Estimation Techniques Applicable to Lithium Poly-Carbon Monoflouride (LI/CFx) Battery", IEEE International Conference on Power Electronics (IICPE), 2010.
  17. O. Erdinc, B. Vural, M. Uzunoglu, "A dynamic lithium-ion battery model considering the effects of temperature and capacity fading", IEEE International Conference on Clean Electrical Power, June 2009.
  18. R. Spotnitz, "Simulation of capacity fade in lithiumion batteries", Journal of Power Sources 113(2003) 72-80.
  19. A123 Systems, "ANR26650M1A Cells Data Sheet".
  20. John Wang, Ping Liu, Jocelyn Hicks-Garner, Elena Sherman, Souren Soukiazian, Mark Verbrugge, Harshad Tataria, James Musser, Peter Finamore, "Cycle-life model for graphite-LiFePO4 cells", Journal of Power Sources, 196 (2011) 3942-3948.
  21. P. Ramadass, Bala Haran, Ralph White, Branko N. Popov, "Mathematical modeling of the capacity fade of Li-ion cells", Journal of Power Sources 123(2003) 230-240.

Cited by

  1. Coordinated Control Algorithm for Distributed Battery Energy Storage Systems for Mitigating Voltage and Frequency Deviations vol.7, pp.3, 2016,
  2. Modeling and Analysis of a Low-Voltage DC Distribution System vol.4, pp.4, 2015,
  3. Estimating Lithium-Ion Battery State of Charge and Parameters Using a Continuous-Discrete Extended Kalman Filter vol.10, pp.8, 2017,
  4. Novel performance of ultrathin AlPO4 coated H2Ti12O25 Exceeding Li4Ti5O12 in cylindrical hybrid supercapacitor vol.273, 2015,
  5. Investigation on Li-ion battery charging polarization characteristics and influence factors vol.22, pp.9, 2016,
  6. New Reclosing Technique in Distribution System with Battery Energy Storage System vol.30, pp.1, 2016,
  7. New Configuration and Novel Reclosing Procedure of Distribution System for Utilization of BESS as UPS in Smart Grid vol.9, pp.4, 2017,
  8. Simulation of electrochemical behavior in Lithium ion battery during discharge process vol.13, pp.1, 2018,
  9. Fabrication and electrochemical properties of cylindrical hybrid supercapacitor using H2Ti12O25 as anode material vol.143, 2015,
  10. Range Extension of Light-Duty Electric Vehicle Improving Efficiency and Power Density of IPMSM Considering Driving Cycle vol.65, pp.12, 2016,
  11. New adaptive reclosing technique using second-order difference of THD in distribution system with BESS used as uninterruptible power supply vol.90, 2017,
  12. Novel adaptive reclosing scheme using wavelet transform in distribution system with battery energy storage system vol.97, 2018,
  13. Optimal Distributed Nonlinear Battery Control vol.5, pp.3, 2017,
  14. Analysis of peak shaving effect of demand power using Vehicle to Grid system in distribution system vol.7, pp.1, 2017,
  15. Reclosing scheme using synchronism checking for utilization of BESS in distribution system 2017,
  16. A Study on the Voltage Sag During the EVs Charging Considering Domestic Data vol.29, pp.1, 2015,
  17. Development of Reclosing Method in a Distribution System with Distributed Generation and Battery Energy Storage System vol.11, pp.6, 2018,