Stepping Simulation of an Example Code (Simulating Battery Pack using Drive Cycle) stuck at 0% with no indication of Simulation Rate.

[Reetesh96]

Hi. I was trying to run one of the Example Codes present in the liionpack Website (Link: https://liionpack.readthedocs.io/en/latest/examples/05%20Drive%20cycles/).

On trying to run the Code, the simulation doesn't start at all and returns a warning: "warn("Matrix is exactly singular", MatrixRankWarning)". How do I rectify this?

The Code that I ran:

`import liionpack as lp
import liionpack as lp
import pybamm
import numpy as np
import os
import pandas as pd
import matplotlib.pyplot as plt

pybamm.set_logging_level("INFO")

drive_cycle = pd.read_csv(
"C:/Users/32231/Desktop/BattBMSModSim/US06Data.csv", comment="#", header=None
).to_numpy()

with plt.rc_context(lp.lp_context()):
plt.figure(figsize=(8, 6))
plt.plot(drive_cycle[:, 0], drive_cycle[:, 1])

experiment = pybamm.Experiment(
operating_conditions=["Run US06 (A)"],
drive_cycles={"US06": drive_cycle},
period="1 second",
)

parameter_values = pybamm.ParameterValues("Chen2020")

netlist = lp.setup_circuit(Np=4, Ns=1, Rb=1e-4, Rc=1e-3, Ri=5e-2, V=4.0, I=1.0)

output_d = lp.solve(
netlist=netlist,
parameter_values=parameter_values,
experiment=experiment,
initial_soc=0.75,
)

lp.plot_pack(output_d)
lp.show_plots()`

The content in the terminal that I am getting on running the code:

PS C:\Users\32231\Desktop\BattBMSModSim> & "C:/Program Files/Python39/python.exe" c:/Users/32231/Desktop/BattBMSModSim/DriveCycleSample Traceback (most recent call last): File "c:\Users\32231\Desktop\BattBMSModSim\DriveCycleSample", line 11, in <module> drive_cycle = pd.read_csv( File "C:\Users\32231\AppData\Roaming\Python\Python39\site-packages\pandas\util\_decorators.py", line 211, in wrapper return func(*args, **kwargs) File "C:\Users\32231\AppData\Roaming\Python\Python39\site-packages\pandas\util\_decorators.py", line 331, in wrapper PS C:\Users\32231\Desktop\BattBMSModSim> & "C:/Program Files/Python39/python.exe" c:/Users/32231/Desktop/BattBMSModSim/DriveCycleSample 2023-02-09 09:42:17.848 - [INFO] parameter_values.process_model(425): Start setting parameters for ElectrodeSOH model 2023-02-09 09:42:17.856 - [INFO] parameter_values.process_model(527): Finish setting parameters for ElectrodeSOH model 2023-02-09 09:42:17.856 - [INFO] discretisation.process_model(148): Start discretising ElectrodeSOH model 2023-02-09 09:42:17.862 - [INFO] discretisation.process_model(260): Finish discretising ElectrodeSOH model 2023-02-09 09:42:17.863 - [INFO] base_solver.solve(704): Start solving ElectrodeSOH model with Algebraic solver (lm) 2023-02-09 09:42:17.863 - [INFO] base_solver.set_up(111): Start solver set-up 2023-02-09 09:42:17.895 - [INFO] base_solver.set_up(236): Finish solver set-up 2023-02-09 09:42:17.901 - [INFO] base_solver.solve(938): Finish solving ElectrodeSOH model (success) 2023-02-09 09:42:17.902 - [INFO] base_solver.solve(939): Set-up time: 32.529 ms, Solve time: 5.376 ms (of which integration time: 4.250 ms), Total time: 37.905 ms 2023-02-09 09:42:17.933 - [INFO] base_model._build_model(550): Start building Single Particle Model 2023-02-09 09:42:17.964 - [INFO] base_battery_model.build_model(995): Finish building Single Particle Model 2023-02-09 09:42:17.965 - [INFO] parameter_values.process_model(425): Start setting parameters for Single Particle Model 2023-02-09 09:42:17.998 - [INFO] parameter_values.process_model(527): Finish setting parameters for Single Particle Model 2023-02-09 09:42:17.999 - [INFO] discretisation.process_model(148): Start discretising Single Particle Model 2023-02-09 09:42:18.001 - [INFO] discretisation.remove_independent_variables_from_rhs(1123): removing variable Discharge capacity [A.h] from rhs 2023-02-09 09:42:18.096 - [INFO] discretisation.process_model(260): Finish discretising Single Particle Model 2023-02-09 09:42:18.097 - [INFO] base_solver.set_up(111): Start solver set-up 2023-02-09 09:42:18.109 - [INFO] base_solver.set_up(236): Finish solver set-up Stepping simulation: 0%| | 0/601 [00:00<?, ?it/s]C:\Users\32231\AppData\Roaming\Python\Python39\site-packages\liionpack\netlist_utils.py:524: RuntimeWarning: divide by zero encountered in divide g[R_map] = 1 / value[R_map] # conductance = 1 / R C:\Users\32231\AppData\Roaming\Python\Python39\site-packages\scipy\sparse\linalg\_dsolve\linsolve.py:276: MatrixRankWarning: Matrix is exactly singular warn("Matrix is exactly singular", MatrixRankWarning) C:\Users\32231\AppData\Roaming\Python\Python39\site-packages\liionpack\netlist_utils.py:524: RuntimeWarning: divide by zero encountered in divide g[R_map] = 1 / value[R_map] # conductance = 1 / R C:\Users\32231\AppData\Roaming\Python\Python39\site-packages\scipy\sparse\linalg\_dsolve\linsolve.py:276: MatrixRankWarning: Matrix is exactly singular warn("Matrix is exactly singular", MatrixRankWarning)

[TomTranter]

This runs ok for me with latest versions of pybamm and liionpack. Check these are up to date and your drive cycle is not corrupted some how. I noticed you changed the file name. I used the one in the PyBaMM folder input directory