%% Batería ANR26650M1 - A123 Systems
clear
clc

Ns = 200; % Número de celdas en serie
Np = 10; % Número de celdas en paralelo

%  Parámetros de la batería:
Q = 2.3; % Capacidad máxima a la temperatura nominal (Ah)
Eo = 3.3678; % Tensión constante de circuito abierto a la temperatura nominal (V)

K = 0.0029; % Constante de polarización a la temperatura nominal (Ohmios)
Rohm = 0.01; % Resistencia óhmica (Ohmios)
Cpol = 0; % Capacidad de polarización (V/Ah)


% Constantes:
Alpha = 1093.6; % Exponente alfa de la constante de polarización (K)
Beta = 8390.6; % Exponente beta de la resistencia óhmica (K)

A = 0.3446; % Tensión exponencial (V)
B = 13.0435; % Capacidad exponencial (1/Ah)

Tf = 10; % Constante del filtro de corriente
Kc = 0.1; % Constante del término 2 para el cálculo de E durante la carga


% Temperaturas:
Ta = 25; % Temperatura ambiente (ºC)
To = 298.15; % Temperatura nominal a la que se han medido todos los parámetros que dependen de esta (K)

% Constantes del modelo térmico:
dE_dT = 0.000010001; % Variación de la tensión de circuito abierto con la temperatura (V/K)
dQ_dT = 0.0037; % Variación de la capacidad con la temperatura (Ah/K)
Rout = 0.6; % Resistencia térmica convectiva (K/W)
Tth = 1000; % Constante de tiempo térmica (s)


% Constantes del modelo de capacidad (cycling):
a = 2.091620116180383e-08; % Término cuadrático del polinomio de temperatura (1/K^2)
b = -1.217949508593764e-05; % Término lineal del polinomio de temperatura (1/K)
c = 0.001802766883936; % Término independiente del polinomio de temperatura 
d = -1.708177699275979e-06; % Constante del término exponencial del acoplamiento entre la temperatura y el Crate (1/K)
e = 0.055578784861400; % Constante del término exponencial del Crate

% Constantes del modelo de capacidad (calendar):
f = 5.980801169706742e06; % Constante del calendar ageing (Ah/días^0.5)
g = 0.689770643186819; % Constante del término exponencial del SOC
h = -6.464711916845157e03; % Constante del término exponencial de la temperatura (1/K)


% Valores críticos de la batería
Umax = 3.7*Ns; % Tensión máxima admisible en bornes de la batería (V)
Umin = 2*Ns; % Tensión mínima admisible en bornes de la batería (V)
Tcrit = 65; % Temperatura máxima admisible en cada celda de la batería (ºC)


% Consignas del ensayo:
Id = Q*3*Np; % Corriente en valor absoluto durante la descarga (A)
Ic = Q*3*Np; % Corriente en valor absoluto durante la carga (A)
SOC0 = 100; % Estado de carga inicial (%)
w = 1.5; % Número de ciclos diarios (ciclos/día)
