App Notes

Here you will find out latest app notes.

Using PSpice to Simulate the Discharge Behavior of Common Batteries

Category: 
Application Notes
Language: 
English

As the use of battery-operated electronic devices become more widespread, so does the need for simulation models used to analyze the operating characteristics of batteries. The most common batteries in use today are: non-rechargeable Alkaline cells, rechargeable Nickel-Cadmium (NICD) cells, Nickel-Metal-Hydride (NIMH) cells, and sealed Lead-Acid cells. This article presents PSpice behavioral models for simulating the four battery types mentioned above.

Battery Variables

All of the battery types modeled here share some common characteristics and deviations from ideal during discharge.

·    Only parameters that vary by 15% or more are considered in these models. Parameters that change less than 15% are not considered because the capacity of cells may vary from +/- 20%  to +/- 50% when shelf time, number of recharge cycles, and manufacturing variances are taken into account.

·    A reduced state of charge is specified for the start because the capacity of a completely charged cell decays with time. For Alkaline cells, this decay takes years to affect the usable capacity. For NICD and Lead-Acid batteries, the decay is 10% to 30% per month. 

Note: The term cell is used to indicate a single energy source. The term battery is used to indicate a power source composed of a single cell or several cells.

·    The usable capacity of a cell varies depending on the discharge rate. At very low discharge rates (<) 100 hours, all batteries are very efficient. At very fast discharge rates (< 10 hours), the batteries are not as efficient and the usable capacity is lost.

·     For pulsed loads with cycle times greater than 10 seconds, a cell gives more total capacity than under a constant load. The rest portion of the pulsed load allows the battery chemistry to recover some of the lost capacity. But, as the pulsed load cycle time becomes less than 1 second, the cell does not have enough time to recover and usable capacity is not increased. In these cases, the RMS value of the pulsed discharge current should be used in the simulation.

·    The affect of cell temperature on cell resistance is not modeled because, for the models, the change in resistance versus temperature falls below the 15% change threshold. The cell temperature changes may be accounted for by adjusting the parameters passed to the various cell subcircuits.

Note: Cell temperature affects both the cell resistance and usable capacity. Low cell temperatures reduce the usable capacity; only a slight decrease is noted at high temperatures. 

    Cell resistance is a function of the cell’s state of charge. Although there is a negligible effect on Lead-Acid and NICD types, Alkaline cells show a 2:1 to 4:1 increase in cell resistance from full charge to full discharge. Still, cell resistance is fairly flat and constant until 80% discharged, when the resistance increases sharply. The sharp fall in cell voltage during discharge can be looked upon as a large increase in cell resistance.

    Cell discharge voltage versus temperature is modeled only for NICD subcircuit because  the open circuit cell voltage variation with temperature even for a 0°C to 60°C range is much less than the difference in actual cell discharge voltage for the other types.

Behavioral Modeling

Figure 1 shows the results of discharging seven identically rated NICD cells to see how well their capacity track. These are old cells, , used every week for 1 to 2 years, and exhibit a 2:1 spread in measured capacity. Alkaline and Lead-Acid batteries have similar variations even between new cells.

The results indicate that there is little practical value in overly accurate models. Therefore, only those battery characteristics that show a 10 to 15% or greater change during discharge are modeled.

dcrgbehbat_fig1.png

                               Figure 1: AA Cell Discharge rate

 

As shown in Figure 2 and Figure 3, the batteries are modeled using the following functional blocks:

  • · Capacitor representing the A-H capacity of the cell.
  • · Discharge rate normalizer to determine the lost capacity at high discharge rates.
  • · A circuit to discharge the A-H capacity of the cell.
  • · Cell voltage versus state-of-charge lookup table.
  • · Cell resistance.
  • · The thermal effects of the cell under high discharge rates for NICD batteries.

To model a cell, perform the following steps:

  • · Generate a parameterized PSpice lookup table
  • · Model discharge current sense and cell resistance
  • · Modeling the State-of-Charge (SOC)

Link SOC with cell voltage to determine output

Generating a Parameterized PSpice Lookup Table

To model a cell, start by generating a parameterized PSpice lookup table Voltage Controlled Voltage Source (VCVS) using the following steps:

  1. Measure several actual discharge curves on a computerized constant-current load analyzer [1] at a low rate (20 to 200 hours) to get an actual voltage versus capacity curve.
  2. Make a single curve either by averaging several curves or by picking a typical curve from the data.
  3. Convert the data into a parameterized PSpice lookup table VCVS. This models the cell’s output voltage versus the state-of-charge at low discharge rates. A simplified VCVS definition is

E_Cell +OUT -OUT TABLE {V(x)} = (0,1.5) (0.5,1.3) (1.0,0.0)

 

Where:

  • E_Cell signifies the PSpice call to a VCVS named E_Cell
  • +OUT and -OUT are the output nodes of the VCVS   
  • TABLE is the PSpice behavioral modeling TABLE directive
  • {V(x)} is the controlling voltage for the table
  • (0,1.5) (0.5,1.3) (1.0,0.0) are the table pairs that are output to +OUT and -OUT based on the value of V(x). If V(x) is 0, signifying 0% discharge, E_Cell will have the first table pair value of 1.5 Volts. If the cell is 50% discharged, the second table pair will be used and so on. For in-between discharge values, PSpice uses linear interpolation between the table pairs.

Note: The actual lookup tables are composed of 30 or more pairs of data to provide finer granularity of the resulting discharge voltage curve.

 

Modeling Discharge Current Sense and Cell Resistance

To model the discharge current sense and the cell resistance, add a zero-valued voltage source in series with the output voltage. The cell resistance is modeled as a simple resistor for NICD or Lead-Acid cells and as a more complex variable resistance that depends on the cell’s state of charge for Alkaline cells.

Modeling the State-of-Charge

To model the state-of-charge (SOC), use a simple and appropriately sized capacitor as the charge storage element that simulates the available charge of the cell. Size this capacitor so that it has a value of 1 Volt at 100% cell capacity and 0.5 Volts at 50% cell capacity. This capacitor is given the following value at the start of the simulation by PSpice’s Parameterization function:

    C_CellCapacity 50 0 {3600*CAPACITY*FudgeFactor}

The capacitor, C_CellCapacity, is connected between nodes 50 and 0 and is given a value of the Amp-hour capacity of the cell times a conversion from hours to seconds (3,600 seconds = 1 hour) times a fudge factor (FudgeFactor). If a cell has a 10 Amphour capacity, C_CellCapacity equals10 * 3,600 or 36,000 Farads; this is a big capacitor, but a workable value that is easy to understand.

FudgeFactor adjusts for the difference in the manufacturer’s listed Amp-hour capacity (i.e., some cutoff voltage with some capacity remaining at the cutoff) and the simulated capacity of 0 Volts output at 0% remaining capacity. To correct for this, and still allow the model user to use the manufacturer’s listed capacity, a FudgeFactor value of 1.01 to 1.1 is included.

The actual usable capacity of a cell depends on the rate at which it is being discharged. Most manufacturers list the capacity at the most favorable rate, usually, at greater than 20 hours discharge. At any faster rate, the cell is less efficient and results in a nonlinear function of the discharge rate. This must be characterized as a lookup table at many discharge rates. This inefficiency is modeled as a VCVS in series with the output voltage of the battery state-of-charge node (the voltage on C_CellCapacity). This VCVS subtracts a given amount of capacity from the cell during discharge. The amount subtracted depends on the rate at which the cell is being discharged.

To determine the rate at which the cell is being discharged, it is convenient to normalize the discharge rate in Amps to a more conventional cell rate called the C rate. The C rate is defined as the capacity of the cell in Amp-hours when it is discharged completely in one hour. This normalization makes it easy to determine the cell inefficiency at different rates, and between different cell sizes, because it converts discharge in Amps to discharge in “C” units of the battery capacity at one hour. This conversion is done in the model by the VCVS, E_Rate, as follows.

E_Rate RATE 0 VALUE = {I(V_Sense) / CAPACITY}

E_Rate is the sensed discharge current in Amps divided by the Amp-hour capacity of the cell. The node, RATE, is the instantaneous rate at which the cell is being discharged, as shown in  Figure 2.

This instantaneous rate information can almost be fed directly to E_Lost_Rate to determine the actual available capacity. But, when the discharge is a low duty cycle, high value pulsed load, the cell supplies a large initial current which decays in seconds to a lower value. For pulsed loads, the cell recovers between pulses and delivers a higher proportion of its capacity than a cell under constant discharge. The delayed rate is modeled by an RC lowpass filter, as shown Figure 2 as  R1 and C1. The exact value of the RC time constant depends on the type and size of cell being simulated. E_Lost_Rate is built like the E_Cell table as follows.

     E_Lost_Rate 50 SOC TABLE { V(x) } = (0.0,0.0) (1.5,0.5)

The table entries indicate the capacity unavailable from the cell at high discharge rates. The table entry shows that at a discharge rate of 0, the cell loses 0% of its capacity (first entry). If the discharge rate is 1.5 times the rated capacity of the cell (1.5 C), the cell loses 50% of its capacity (second entry in the table).

In Figure 2, the State-Of-Charge (SOC) node is the subtraction of the voltage on the capacitor C_CellCapacity and E_Lost_Rate. The SOC node represents the capacity in the cell for a given discharge rate during the simulation. G_Discharge discharges C_CellCapacity at the cell rate. The voltage on node 50 relates to the capacity remaining in the cell if the discharge rate is low enough to actually run the cell dry. At low discharge rates, these two nodes are the same; at high discharge rates, node SOC is at a lower potential than node 50. If, at the end of a high discharge rate the cell reverts to a low discharge, nearly the entire rated capacity can be recovered from the cell. At the high discharge rate, approximately 60% of the cell’s rated capacity can be used.

 

Linking SOC with Cell Voltage to Determine Output

You can now determine the output by linking  the state of charge with the cell voltage to determine an output. The state of charge is 1 Volt for 100%, while the cell voltage table is just the opposite. To make the cell voltage correct, the state-of-charge voltage must be inverted as shown in Figure 2. The schematic in Figure 2 represents the subcircuit model common to all batteries without temperature effect.

 

dcrgbehbat_fig2.png

Figure 2: Functional schematic  developed for all of teh modeled cell types, only minor changes are required to complete each detailed model type

 

Model Differences for Different Battery Types

Alkaline cells

Alkaline cell resistance is not fixed throughout the discharge.The resistance model is developed by determining the relationship of the output’s current and voltage, then linking this to the battery’s state of charge. The cell has a small resistance increase from 100% to 20% cell capacity, then increases to twice its initial value at 0% capacity.

When the discharge current approaches 100mA, the discharge capacity versus discharge rate produces a kink in the discharge curve. Below 100mA, the cell loses capacity gradually; above 100mA, the rate of lost capacity increases significantly. Although the 100mA discharge rate is the same for all of the cell sizes, the C rating is not. Because it is not possible to relate this kink to a specific C discharge rate, a separate E_Lost_Rate table must be developed for each cell size.

A separate subcircuit model is used to model the 9V Alkaline cell. The cell resistance change versus discharge state is more pronounced in this type of battery, and is modified accordingly in the model.

 

Nickel-cadmium cells

These cells are often used at very high discharge rates up to 20C. Discharging a fully charged cell in 5 or 6 minutes (10C rate) releases significant amounts of heat. To account for this, a thermal model is included in the cell subcircuit.

The thermal temperature rise of a cell per watt dissipated in free air is approximately Θcell = 13.41 * V-0.61 where V is the cell volume in cubic inches, and Θcell is the thermal rise of the cell in °C per watt dissipated.The thermal time constant for material of the density used in making NICD batteries is approximately 20 minutes per pound or expressed in more convenient terms, 2.65 seconds per gram.

These empirical relationships are used with the calculated cell power dissipation (Cell Discharge Current2 * Cell Resistance) to get a temperature rise and time constant model for the cell temperature. The cell temperature rise above ambient temperature is available at node CELL_TEMP in the NICD model. The temperature information is also used to add or subtract from the cell discharge voltage to account for the cell temperature E_Temp in the NICD model.

Another small modification is made to the NICD model to facilitate the direct entry of manufacturers’ rated capacity data. Most NICD batteries are not rated at their maximum capacity for low discharge rates. The norm is to rate them at the C to C/5 rate, leaving 30% more than the rated capacity if the cell is used at low discharge rates. To account for this difference a Voltage- Controlled Current Source (VCCS), G_LowRate is used to add a small amount of current to C_CellCapacity during discharge at rates less than 1C.

 

Nickel-Metal-Hydride cells

These cells are modeled like the NICD cells, but without the fast discharge thermal effects.

Lead-acid cells

Since these cells are almost universally used in batteries composed of 3, 6, or more cells (6 or 12 Volt batteries), the model is changed slightly. The single cell voltage is multiplied by the number of cells to get the total battery voltage.The schematic in Figure 3 refers to the schematic shown in Figure 2 implemented for NICD battery with an additionally added cell temperature block.

 

dcrgbehbat_fig3.png

                                 Figure 3: NICD battery discharge models

 

Using the Discharge Models

To use the battery models in a simulation, add the required number of cells to the circuit file, and pass the appropriate parameters to the subcircuit model. The manufacturers’ data sheets may be consulted for information on the cell parameters shown in Table 1.

 

Cell Model Parameter Meaning
ALKALINE CAPACITY Capacity of the cell measured at greater than 100 hour discharge time
  RESISTANCE Initial cell resistance in ohms
     
9V ALKALINE   No parameters are needed for a generic 9 Volt Alkaline battery
     
NICKELCADMIUM CAPACITY Amp-hour capacity of cell measured at the C/5 rate
  RESISTANCE Mid-discharge cell resistance in ohms
  CELLTEMP1 Initial cell temperature at start of simulation
  VOLUME1 Volume of cell in cubic inches
  WT1 The weight of the cell in grams
     
NICKEL METAL HYDRIDE CAPACITY Amp-hour capacity of cell measured at the C/5 rate
  RESISTANCE Mid-discharge cell resistance in ohms
     
LEAD-ACID CAPACITY Capacity of the battery measured at the 20 hour rate
  RESISTANCE Mid-discharge cell resistance in ohms
  CELLS The number of cells that make up the battery(6V = 3 cells, 12V = 6, etc.)


Note: These parameters are only needed at a discharge rate greater than 5C to account for cell temperature rise.

Since little standardization exists for rating methods between manufacturers, the popular typical parameters are summarized in Table 2 for the most common consumer batteries

 

Cell Type Capacity(Amp-hour) Resistance(Ohms) Volume(in3) Weight(gm) Cell Number
Alkaline Cells    
Typical subcircuit call for a single Alkaline N cell:
X1 +node -node SOC RATE ALKALINE PARAMS: CAPACITY=0.9,RESISTANCE=5
N 0.0 0.8 - - -
AAA 1.2 0.6 - - -
AA 2.5 0.3 - - -
C 7.5 0.2 - - -
D 16.4 0.07 - - -
Nickel-Cadmium Cells (Standard)1
Typical subcircuit call for a single NICD N cell:
X1 +node -node SOC RATE CELL_TEMP NICD PARAMS: CAPACITY=0.15, RESISTANCE=5, VOLUME=0.2, WT=9
N             0.15 0.027 0.2 9 -
AAA 0.18 0.021 0.24 10 -
AA 0.48 0.012 0.48 24 -
SUBC 1.2 0.005 1.1 50 -
C2 1.8 0.0045 1.6 80 -
D2 4.0 0.0035 3.4 160 -
Nickel-Cadmium Cells (Standard)1
Typical subcircuit call for a single NICD N cell:
X1 +node -node SOC RATE CELL_TEMP NICD PARAMS: CAPACITY=0.15, RESISTANCE=5, VOLUME=0.2, WT=9
AA 1.1 0.03 - - -
4/5A 1.5 0.02(estimated) - - -
Lead-Acid Cells
Typical subcircuit call for a 6 Volt, 1.3 Amp-hour Lead-Acid battery:
X1 +node -node SOC RATE LEADACID PARAMS: CAPACITY=1.3,RESISTANCE=0.06,CELLS=3
6V-1.3A-hr 1.3 0.06 - - 3
6V-4.0A-hr 4.0 0.025 - - 3
6V-6.5A-hr 6.5 0.02 - - 3
6V-10A-hr 10 0.015 - - 3
12V-1.3A-hr 1.3 0.12 - - 6
12V-4.0A-hr 4.0 0.05 - - 6
12V-6.5A-hr 6.5 0.04 - - 6
12V-10A-hr 10 0.03 - - 6

 

1. Use the Volume and Weight parameters for NICD cells only to simulate a discharge rate grater than approximately 5C and when the temperature profile of the cell is needed.
2. Although real NICD C and D cells can be purchased, most consumer C and D batteries are actually SUB C cells in a big empty can.

 

Example Circuit of NICD AA and N Cell Discharge Test

To demonstrate the use of a cell in an actual circuit a real AA NICD cell is discharged into a constant resistance of 2Ohms. 

Then the circuit shown in Figure 4 is used with PSpice to simulate the 2Ohm and 5Ohm discharge. Running this circuit with the capacity of the NICD model normalized to the actual capacity of the measured cell, shows results which compare very favorably with the real behavior.

 

dcrgbehbat_fig4.png

         Figure 4: Schematic for NICD AA and N cell Discharge test

 

dcrgbehbat_fig5.png

                           Figure 5: Simulation result of NICD AA and N cell Discharge

 

N cell is simulated by changing the properties of the cell. The .IC statement sets the initial conditions, and must be set for every subcircuit used. V(X1.50) sets the initial charge on node 50 of the X1 subcircuit. This is the voltage on the battery Amp-hour capacity model which simulates the initial state of charge. Setting this node to 1 Volt equals an initial state of charge of 100%. Likewise, 0.8 Volts would represent an 80% initial state of charge.

The next initial condition ( V(X1.60)=0) sets the voltage on the delayed lost rate calculator to zero. This allows the voltage on capacitor C1 (internal to the subcircuit) to start at 0 Volts as it would if the discharge current was zero before the simulation started. Another way to achieve this result is to switch on the discharge currents just after the simulation starts. This automatically sets the delayed lost rate voltage to zero at the start of the simulation.

Note:To perform simulation of other cell models use the .olb file.

 

Temperature Effects

During a fast discharge, the cell temperature of a NICD can change by 25°C or more. The effect of cell temperature on voltage is accounted for in the NICD model by changing the cell voltage based on the calculated temperature. The other models do not incorporate temperature effects directly into the simulation.

The major temperature influence on the cell is capacity. This may be accounted for by adjusting this parameter at the start of a simulation. The following equations give the new capacity for each cell type at any discharge temperature from 0 to 60°C based on the initial capacity at 25°C.

1.Alkaline Cells:
    NewCapacity = OldCapacity * (0.85 + 8.64E-3 * T - 1.05E-4 * T*T)

2.Nickel-Cadmium Cells:
If T > 25°C     NewCapacity=OldCapacity
If T < 25°C     NewCapacity=OldCapacity*(0.815+7.5E-3*T)

3.Nickel-Metal-Hydride Cells:
    NewCapacity=OldCapacity*(0.913+1.1E-2*T-3.0E-4*T*T)

4.Lead-Acid Cells:
    NewCapacity=OldCapacity*(0.84 + 7.96E-3*T-6.07E-5*T*T)

 

To use the equations, simply plug in the 25°C initial capacity for OldCapacity and the new discharge temperature for T into the proper equation. If needed, these equations can also be built into the subcircuit models.

 

dcrgbehbat_fig6.png

Figure 6: Schematic for Alkaline AA Cell Discharge test

dcrgbehbat_fig7.png

                                  Figure 7: Simulation result of Alkaline AA Cell

 

dcrgbehbat_fig8.png

                  Figure 8: Simulation result of Alkaline AA Cell at different temperatures

Figure 8 shows the simulation results for the discharging rate of Alkaline AA cell at 0, 25, and 100°C temperatures corresponding to the green,red and blue color traces respectively.

 

Limitations of the Models

The circuits presented in this chapter trade off accuracy with simplicity and simulation time. There are several cases where the simulated and actual results vary significantly.

  • Rest time of Alkaline batteries not considered: In practice, Alkaline cells are sometimes used in 2 to 8 hour shifts followed by rest time for the remainder of the day. This rest time allows the cell to recover part of its discharge capacity. This phenomenon is not specifically modeled, and when simulated with the alkaline subcircuit, the observed capacity may be up to 25% short over a cell’s actual performance. Using the Alkaline models in these patterns will give a conservative estimate of capacity.
  • Internal resistance at low discharge rates to low terminals not considered: When a battery is discharged to a low terminal voltage level and then disconnected from the load, the battery voltage will recover to some higher level in an hour or so. This phenomenon is accounted for in these models when they are discharged at high current levels then left to rest. When the models are used at low discharge rates to low terminal voltages however, they do not show this voltage recovery. The battery chemistry tries to make a voltage potential difference even if only a few molecules of unused material remain. In this state of discharge, the internal resistance of the cell can be an order of magnitude, or more, than its initial value. If any load is reconnected, the terminal voltage will quickly collapse again to zero.
  • Usage beyond standard cutoff voltages not accounted for: The models were designed to be used with the standard cutoff voltages as specified by the battery makers. For NICD batteries, this is 0.8 to 1.1 Volts per cell. For Alkaline cells, the cutoff is 0.8 to 1.2 Volts per cell. For Lead-Acid cells, the cutoff voltage is typically 1.5 to 1.7 Volts per cell. Usage beyond these limits should be studied carefully because they are not specifically accounted for in the modeling process.

  • Models not designed for parallel connection: These models are not designed to be connected in parallel.This is not acceptable in consumer design anyway, because there is no way to guard against the end user putting an Alkaline cell in parallel with a NICD or Carbon cell. If these different types of cells are connected together, the charge most likely will not equalize. This results in overcharging and leakage of the weak cell,  causing damage. 

 

Simulation Speed

The goal of simulation is to obtain results faster than can be achieved with the hardware, or to measure behaviors which cannot easily be accessed in the hardware. The following notes should help when making speed/accuracy trade-offs.

  • Be careful with the models that you attach to these batteries during a simulation. Simulate the power drain from your circuit, not the transistor-level circuit itself. These models have been tried and verified with many different discharge regimens, and are believed to be accurate enough to allow finding maximum or minimum battery life. The trends these models simulate are believed to be basically accurate, even though the absolute capacity simulated may be 20% or so off.
  • Do not simulate pulsed current loads with cycle times less than 5 seconds. Using short cycle time pulsed currents may make the simulation run slower than real time. To speed up the simulation with fast pulsed loads, use the RMS average of the pulsed current. This will provide you with a ballpark answer.
  • Use a minimum of semiconductor models hooked up to these models. Semiconductors contain many internal nonlinear equations that must be solved for each time point, thus slowing simulation time.
  • If you experience convergence problems, use the .IC directive to set the initial voltages on critical nodes in your circuit.

 

Annexure

The following are the five different Cell discharge models.These models are presented in the library file battery.lib and its associated symbols are presented in the battery.olb file.These .lib and .olb files can be used to perform discharge test of other batteries.

** PSpice Alkaline battery discharge model **
    .SUBCKT ALKALINE +OUTPUT -OUTPUT SOC RATE + PARAMS: CAPACITY=1, RESISTANCE=1
    C_CellCapacity 50 0 { 3600 * CAPACITY * 1.01 }
    R1 50 0 1G
    E_Resistance 20 10 VALUE = {I(V_Sense) * RESISTANCE * V(Cell_Res)}
    E_Cell_R Cell_Res 0 TABLE { V(50) } = (0,2) (0.2,1) (1,1)
    R3 Cell_Res 0 1G
    V_Sense -OUTPUT 20 0
    E_Invert Invert 0 TABLE { V(SOC) } = (0,1) (1,0)
    R4 Invert 0 1G
    E_Rate RATE 0 VALUE = { I(V_Sense)/CAPACITY }
    R2 RATE 60 10 
    C1 60 0 1
    G_Discharge SOC 0 VALUE = { I(V_Sense) } 
    E_Lost_Rate 50 SOC TABLE { V(60) } =
* * Use one of the following tables!!! * *
;----- Use this table for N cells -----
;+ (0.0,0.0) (0.019,0.056) (0.043,0.13) (0.072,0.28)
;+ (0.12,0.39) (0.21,0.58) (0.31,0.69)
;----- Use this table for AAA and AA cells -----
;+ (0.0,0.0) (0.018,0.08) (0.043,0.14) (0.08,0.2)
;+ (0.14,0.3) (0.26,0.48) (0.4,0.6)
;----- Use this table for C cells -----
;+ (0.0,0.0) (0.17,0.13) (0.035,0.31) (0.055,0.45)
;+ (0.093,0.53) (0.17,0.65) (0.27,0.73)
;----- Use this table for D cells -----
;+ (0.0,0.0) (0.0091,0.091) (0.017,0.15) (0.032,0.36)
;+ (0.058,0.42) (0.079,0.61) (0.13,0.73) (0.18,0.82)
    E_Cell +OUTPUT 10 TABLE { V(Invert) }=
+ (0.000E+00 1.528E+00) (2.320E-03 1.511E+00) (4.640E-03 1.500E+00)
+ (9.280E-03 1.481E+00) (1.392E-02 1.468E+00) (1.856E-02 1.457E+00)
+ (2.552E-02 1.442E+00) (3.248E-02 1.430E+00) (3.944E-02 1.419E+00)
+ (4.872E-02 1.406E+00) (5.800E-02 1.394E+00) (6.728E-02 1.380E+00)
+ (7.656E-02 1.370E+00) (1.206E-01 1.326E+00) (2.691E-01 1.230E+00)
+ (5.522E-01 1.126E+00) (8.213E-01 1.021E+00) (9.025E-01 9.901E-01)
+ (9.257E-01 9.792E-01) (9.443E-01 9.676E-01) (9.559E-01 9.564E-01)
+ (9.628E-01 9.445E-01) (9.698E-01 9.299E-01) (9.744E-01 9.181E-01)
+ (9.791E-01 9.043E-01) (9.814E-01 8.937E-01) (9.837E-01 8.800E-01)
+ (9.860E-01 8.654E-01) (9.884E-01 8.470E-01) (9.907E-01 8.040E-01)
+ (9.930E-01 6.417E-01) (9.953E-01 3.795E-01) (9.976E-01 3.354E-01)
+ (1.0 0.0)
    .ENDS

 

** PSpice 9 Volt Alkaline battery discharge model **
    .SUBCKT ALK_9V +OUTPUT -OUTPUT SOC RATE + PARAMS: CAPACITY=0.565
    E_Rate RATE 0 VALUE = { I(V_Sense)/CAPACITY }
    R2 RATE 60 10
    C1 60 0 1
    G_Discharge SOC 0 VALUE = { I(V_Sense) }
    E_Lost_Rate 50 SOC TABLE { V(60) } = + (0.0,0.0) (0.025,0.009) (0.046,0.080) (0.088,0.14) (0.18,0.21) (0.71,0.45)

    C_CellCapacity 50 0 { 3600 * CAPACITY * 1.06 }
    R1 50 0 1G
    E_Resistance 20 10 VALUE = { I(V_Sense) * 2.0 * V(Cell_Res) }
    E_Cell_R Cell_Res 0 TABLE { V(50) } = (0,4) (0.2,2) (1,1)
    R3 Cell_Res 0 1G
    V_Sense -OUTPUT 20 0
    E_Invert Invert 0 TABLE { V(SOC) } = (0,1) (1,0)
    R4 Invert 0 1G

 

E_Cell +OUTPUT 10 TABLE { V(Invert) } =
+ (0.00 9.18) (0.05 8.82) (0.10 8.62)
+ (0.15 8.41) (0.20 8.30) (0.25 8.21)
+ (0.30 8.09) (0.35 7.99) (0.40 7.95)
+ (0.45 7.89) (0.50 7.79) (0.55 7.66)
+ (0.60 7.55) (0.70 7.18) (0.75 6.96)
+ (0.80 6.58) (0.85 6.12) (0.90 5.42)
+ (0.95 4.51) (1.00 0.00)
    .ENDS

 

** PSpice Nickel-Cadmium battery discharge model **
    .SUBCKT NICD +OUTPUT -OUTPUT SOC RATE CELL_TEMP + PARAMS: CAPACITY=1, RESISTANCE=1, CELLTEMP=25, VOLUME=1, WT=1
    E_Rate RATE 0 VALUE = { I(V_Sense)/CAPACITY }
    R2 RATE 60 1 
    C1 60 0 3
    E_LowRate LowRate 0 TABLE {V(RATE)} = (0,0) (0.001,0.2) (0.1,0.2) (1,0)
    R3 LowRate 0 1G
    G_LowRate 0 50 VALUE = { V(LowRate)*I(V_Sense) }
    G_Discharge SOC 0 VALUE = { I(V_Sense) } 
    E_LostRate 50 SOC TABLE { V(60) } = (1.0,0) (10,0.25)
    C_CellCapacity 50 0 { 3600 * CAPACITY * 1.03 }
    R1 50 0 1G
    E_Temp_Rise 70 0 VALUE = {PWR(I(V_Sense),2)*RESISTANCE*(13.4 PWR(VOLUME,-0.6065))}
    V_Ambient 80 70 { CELLTEMP }
    R_Thermal CELL_TEMP 80 { 2.65 * WT }
    C2 CELL_TEMP 0 1
    R_Cell 20 30 { RESISTANCE }
    E_Temp 10 20 TABLE { V(CELL_TEMP) } = (0,-0.025) (25,0) (60,-0.100)
    V_Sense -OUTPUT 30 0
    E_Invert Invert 0 TABLE { V(SOC) } = (0,1) (1,0)
    R4 Invert 0 1G
    E_Cell +OUTPUT 10 TABLE { V(Invert) } =

 

+ (0.0000000000E+00 1.3148600000E+00) (1.7391197842E-03 1.3114600000E+00)
+ (8.6956352158E-03 1.3008400000E+00) (1.7391252284E-02 1.2910200000E+00)
+ (3.1304265000E-02 1.2794000000E+00) (4.8695517284E-02 1.2685600000E+00)
+ (6.4347649784E-02 1.2608600000E+00) (9.2173675215E-02 1.2504800000E+00)
+ (1.3739093478E-01 1.2401000000E+00) (2.9217314000E-01 1.2300000000E+00)
+ (4.9738998228E-01 1.2199000000E+00) (6.0347665207E-01 1.2099000000E+00)
+ (7.3738938358E-01 1.1909200000E+00) (7.8782396620E-01 1.1801600000E+00)
+ (8.2086740543E-01 1.1705000000E+00) (8.4695429293E-01 1.1602800000E+00)
+ (8.6608462870E-01 1.1501000000E+00) (8.7999764142E-01 1.1403000000E+00)
+ (8.9043241457E-01 1.1309800000E+00) (8.9912806793E-01 1.1212400000E+00)
+ (9.0782372129E-01 1.1090600000E+00) (9.1304106250E-01 1.1001000000E+00)
+ (9.1825840370E-01 1.0899000000E+00) (9.2347583564E-01 1.0784600000E+00)
+ (9.2695405706E-01 1.0705000000E+00) (9.3217141642E-01 1.0586000000E+00)
+ (9.3564971043E-01 1.0508000000E+00) (9.3912795000E-01 1.0430000000E+00)
+ (9.4434529120E-01 1.0303200000E+00) (9.4782360336E-01 1.0207400000E+00)
+ (9.5130184293E-01 1.0102600000E+00) (9.5304094457E-01 1.0046200000E+00)
+ (9.5478006435E-01 9.9866000000E-01) (9.5651918413E-01 9.9248000000E-01)
+ (9.5825830392E-01 9.8596000000E-01) (9.5999749629E-01 9.7910000000E-01)
+ (9.6173659793E-01 9.7178000000E-01) (9.6347571771E-01 9.6386000000E-01)
+ (9.6521483750E-01 9.5520000000E-01) (9.6695395728E-01 9.4556000000E-01)
+ (9.6869307707E-01 9.3460000000E-01) (9.7043217870E-01 9.2192000000E-01)
+ (9.7217137108E-01 9.0686000000E-01) (9.7391049086E-01 8.8908000000E-01)
+ (9.7564961064E-01 8.6722000000E-01) (9.7738873043E-01 8.3990000000E-01)
+ (9.7912783206E-01 8.0636000000E-01) (9.8086695185E-01 7.6520000000E-01)
+ (9.8260607163E-01 7.1436000000E-01) (9.8434519142E-01 6.6000000000E-01)
+ (9.8608438379E-01 6.0778000000E-01) (9.8782348543E-01 5.5698000000E-01)
+ (9.8956260521E-01 5.0776000000E-01) (9.9130172500E-01 4.5810000000E-01)
+ (9.9304084478E-01 4.0860000000E-01) (9.9477996457E-01 3.5850000000E-01)
+ (9.9651906620E-01 3.0526000000E-01) (9.9825825858E-01 2.4604600000E-01)
+ (9.9999737836E-01 1.8616600000E-01) (1.0000000000E+00 0.0000000000E+00)
    .ENDS

 

* PSpice Lead-Acid battery discharge model **
    .SUBCKT LEADACID +OUTPUT -OUTPUT SOC RATE + PARAMS: CAPACITY=1, RESISTANCE=1, CELLS=3
    E_Rate RATE 0 VALUE = { I(V_Sense)/CAPACITY }
    R2 RATE 60 60 
    C1 60 0 1
    G_Discharge SOC 0 VALUE = { I(V_Sense) }
    E_Lost_Rate 50 SOC TABLE { V(60) } = + (0.05,0.0) (0.089,0.11)(0.16,0.20)(0.62,0.39)(0.8,0.47) (1.6,0.44)
    C_CellCapacity 50 0 { 3600 * CAPACITY * 1.15 }
    R1 50 0 1G
    R_Cell 10 20 { RESISTANCE }
    E_Battery +OUTPUT 10 VALUE = { V(Cell_V) * CELLS }
    V_Sense -OUTPUT 20 0
    E_Invert Invert 0 TABLE { V(SOC) } = (0,1) (1,0)
    R4 Invert 0 1G
    R5 Cell_V 0 1G
    E_Cell Cell_V 0 TABLE { V(Invert) } =

+ (0.000E+00 2.171E+00) (5.222E-04 2.149E+00) (1.828E-03 2.128E+00)
+ (1.263E-01 2.101E+00) (4.908E-01 2.001E+00) (6.385E-01 1.949E+00)
+ (7.459E-01 1.900E+00) (7.834E-01 1.875E+00) (8.117E-01 1.850E+00)
+ (8.313E-01 1.826E+00) (8.436E-01 1.801E+00) (8.517E-01 1.773E+00)
+ (8.556E-01 1.750E+00) (8.591E-01 1.724E+00) (8.616E-01 1.702E+00)
+ (8.646E-01 1.676E+00) (8.677E-01 1.648E+00) (8.707E-01 1.623E+00)
+ (8.732E-01 1.600E+00) (8.850E-01 1.499E+00) (8.965E-01 1.401E+00)
+ (9.000E-01 1.333E+00) (1.000E+00 0.000E+00)
    .ENDS

 

** PSpice Nickel-Metal-Hydride battery discharge model **
    .SUBCKT NIMH +OUTPUT -OUTPUT SOC RATE + PARAMS: CAPACITY=1, RESISTANCE=1
    E_Rate RATE 0 VALUE = { I(V_Sense)/CAPACITY }
    R2 RATE 60 1
    C1 60 0 3

 

    E_LowRate LowRate 0 TABLE {V(RATE)} = (0,0) (0.001,0.15) (0.1,0.1) (0.2,0)
    R3 LowRate 0 1G
    G_LowRate 0 50 VALUE = { V(LowRate)*I(V_Sense) }
    G_Discharge SOC 0 VALUE = { I(V_Sense) }
    E_LostRate 50 SOC TABLE { V(60) } = (0.2,0.0) (1.0,0.15) (5,0.2)
    C_CellCapacity 50 0 { 3600 * CAPACITY * 1.01 }
    R1 50 0 1G
    R_Cell 20 30 { RESISTANCE }
    V_Sense -OUTPUT 30 0
    E_Invert Invert 0 TABLE { V(SOC) } = (0,1) (1,0)
    R4 Invert 0 1G
    E_Cell +OUTPUT 20 TABLE { V(Invert) } =
+(0.0000E+00 1.3346E+00) (7.0989E-03 1.3244E+00) (1.6327E-02 1.3144E+00)
+(2.9283E-02 1.3042E+00) (4.2593E-02 1.2942E+00) (6.8859E-02 1.2841E+00)
+(1.3008E-01 1.2733E+00) (4.3605E-01 1.2633E+00) (5.1165E-01 1.2532E+00)
+(5.8033E-01 1.2432E+00) (6.4635E-01 1.2331E+00) (7.0190E-01 1.2231E+00)
+(7.5834E-01 1.2130E+00) (8.0324E-01 1.2030E+00) (8.3075E-01 1.1929E+00)
+(8.5116E-01 1.1828E+00) (8.6820E-01 1.1727E+00) (8.8310E-01 1.1627E+00)
+(8.9641E-01 1.1527E+00) (9.0848E-01 1.1425E+00) (9.1860E-01 1.1324E+00)
+(9.2730E-01 1.1223E+00) (9.3475E-01 1.1122E+00) (9.4167E-01 1.1021E+00)
+(9.4841E-01 1.0919E+00) (9.5480E-01 1.0817E+00) (9.6013E-01 1.0716E+00)
+(9.6439E-01 1.0615E+00) (9.6776E-01 1.0515E+00) (9.7060E-01 1.0407E+00)
+(9.7291E-01 1.0299E+00) (9.7486E-01 1.0190E+00) (9.7663E-01 1.0080E+00)
+(9.7823E-01 9.9782E-01) (9.8001E-01 9.8706E-01) (9.8196E-01 9.7630E-01)
+(9.8391E-01 9.6612E-01) (9.8586E-01 9.5606E-01) (9.8799E-01 9.4542E-01)
+(9.9012E-01 9.3524E-01) (9.9225E-01 9.2518E-01) (9.9420E-01 9.1498E-01)
+(9.9580E-01 9.0400E-01) (9.9687E-01 8.9186E-01) (9.9740E-01 8.7990E-01)
+(9.9775E-01 8.6280E-01) (9.9793E-01 8.4818E-01) (9.9811E-01 8.2718E-01)
+(9.9828E-01 7.9518E-01) (9.9846E-01 7.4066E-01) (9.9864E-01 6.4712E-01)
+(9.9882E-01 5.1380E-01) (9.9899E-01 3.3476E-01) (1.0000E+00 0.0000E+00)
    .ENDS

 

References

1 Hageman, S.C. “PC-based power supply tester slashes setup time to minutes,” Personal Engineering & Instrumentation News, December 1990, Pages 65-66.

2. MicroSim, Circuit Analysis Reference Manual, Version 5.3,1993, MicroSim Corporation, Irvine, CA. 

3. DURACELL, Alkaline Dioxide Batteries, DURACELL INC., Bethel, CT.

4. SANYO, Various literature on Nickel-Cadmium batteries,Sanyo Electric Inc., Little Ferry, NJ.

5. Panasonic, Sealed Lead-Acid Batteries - Technical Handbook, Panasonic Industrial Co. Secaucus, NJ.

6. Panasonic, Batteries, Panasonic Industrial Co. Secaucus,NJ.

7. GATES, Rechargeable Batteries - Applications Handbook,1992, ISBN 0-7506-9228-6, Butterworth-Heinemann,Stoneham, MA

8. General Electric, Nickel-Cadmium Battery Application Engineering Handbook, 1975, General Electric Company.

9. Power-Sonic, Battery Selector Guide, Power-Sonic Corp.,Redwood City, CA.

 

© Copyright 2016 Cadence Design Systems, Inc. All rights reserved. Cadence, the Cadence logo, and Spectre are registered trademarks of Cadence Design Systems, Inc. All others are properties of their respective holders.

Download PSpice and try it for free! Download Free Trial
Cadence