CC BY
94
УДК 621.355: 519.713
М.М. Хандорин, В.Г. Букреев
Методика оценки остаточной емкости литий-ионной аккумуляторной батареи при изменении температуры аккумулятора
Дана методика оценки остаточной емкости при изменении температуры литий-ионной аккумуляторной батареи. Для описания процессов заряда-разряда использована модель Тевенина, учитывающая изменение температуры батареи. Представлены результаты экспериментальной проверки предложенной методики. Ключевые слова: литий-ионный аккумулятор, модель Тевенина, сплайн-интерполяция, линейная интерполяция, микроконтроллер. ао1: 10.21293/1818-0442-2017-20-2-120-122
В современной космической технике все шире применяются литий-ионные аккумуляторные батареи. Разработчиков космических аппаратов привлекают их высокая энергоемкость и длительный ресурс эксплуатации. Для прогнозирования достаточности запаса энергии на борту наиболее важны степень заряженности SOC(t) (State of Charge), остаточная емкость Q\(t) в ампер-часах и коэффициент технической годности SOH(t) (State of Health). Особенно это актуально для космических аппаратов с переменной электрической нагрузкой в течение сеанса, так как при изменении нагрузки в небольших пределах энергобаланс аппарата, как правило, может быть рассчитан заранее [1, 2].
Разработанная методика обладает низкой вычислительной сложностью и может применяться в том числе в системах электроснабжения малых космических аппаратов, аналогичных [3].
Краткое описание методики оценки остаточной емкости
Способ оценки остаточной емкости батареи был ранее описан авторами [4, 5], здесь приведем лишь краткие сведения.
Для оценки остаточной емкости аккумулятора используется модифицированная модель Тевенина [4, 6], записанная следующим образом:
Uак (t) =1 (¿ак (t) R S +Eo ~ ^к (t)) , (1)
Т
где RS и т - полное внутреннее сопротивление и постоянная времени аккумулятора; E0 - его ЭДС.
Определение параметров модели Тевенина на всей разрядной характеристике батареи проводится путем интерполяции данных параметров. В качестве узлов интерполяции используются фиксированные значения параметров в точках разрядной характеристики SOC(t). Для интерполяции E0 применяется кубический сплайн, в то же время интерполяция параметров RS, Т основана на линейной функции. Отказ от применения кубического сплайна для интерполяции параметров RS, т связан с большим разбросом значений между узлами, что ведет к снижению точности и появлению недопустимой ошибки.
Вычисление степени заряженности в j-й момент времени базируется на выражении
SO Cj = SO Cj_i + k ((ак j._i -Сак j_i ), (2) измеренное напряжение аккумулятора;
где иак
иак - оценка его напряжения на основе модели;
к - коэффициент пропорциональности, определяемый из соотношения
1
к=-
10Um
(3)
где U_
напряжение полностью заряженного
аккумулятора.
Способ определения параметров модели аккумулятора по реакции на импульсное воздействие достаточно подробно описан в [4].
Для модели (1) используется линейная зависимость параметров Е0, т и полной емкости аккумулятора Q0 от текущей температуры Т. Температурная же зависимость полного сопротивления представляется нелинейной функцией [7]
Rs (T) = Rs (Tref )exp
( ( Kr
1
Y\
T
ref
(4)
где Тгеу - нормальная температура аккумулятора; КК - коэффициент, определяемый на основе экспериментальных данных.
Вычисление Кк проводится по формуле
ln
RS (Tmin)
Kr =-
RS (Tref )
(5)
T
1 m
T
ref
где Ттш - минимально допустимая температура батареи.
Таким образом, для учета изменений характеристик батареи во всем температурном диапазоне достаточно определить параметры модели Тевенина при двух предельных значениях температуры.
М.М. Хандорин, В.Г. Букреев. Методика оценки остаточной емкости литий-ионной аккумуляторной батареи
121
Описание алгоритма вычисления параметров модели Тевенина с учетом температуры
Алгоритм оценки степени заряженности и остаточной емкости батареи в целом повторяет ранее опубликованный алгоритм [5]. При этом существенным изменением, сделанным в настоящей работе, является обновление параметров модели Тевенина в каждом цикле измерений с учетом текущей температуры батареи.
В памяти контроллера сохраняются параметры модели Тевенина при нормальной Tref и при минимальной Tmin температуре аккумулятора. Расчет параметров модели Тевенина для текущей температуры проводится в три этапа:
- вычисление коэффициентов линейных функций для параметров E0, т, Q0;
- вычисление коэффициента KR по формуле (5);
- вычисление текущих значений E0(T), т(Т), Rs(T), Qo(T).
Полученные значения параметров модели Теве-нина используются в алгоритме оценки остаточной емкости батареи.
При нагрузке номинальным током батарея не может быть полностью разряжена, и часть емкости Qa(I) не используется. Это происходит потому, что часть напряжения батареи падает на ее внутреннем сопротивлении RS, и попытка продолжить разряд приведет к падению напряжения на ней ниже минимально допустимого уровня Umin. Для вычисления Qa(I) применяется следующий алгоритм:
1) из уравнения (1) вычисляется напряжение батареи при текущем токе разряда для различных значений SOC;
2) по полученным значениям выполняется бинарный поиск напряжения Umin для определения минимально возможной степени заряженности SOCmin с данным значением тока I;
3) вычисляется Qa(I) по формуле
Qa = SOCmin Qo. (6)
Обновление параметров модели при деградации
При изменении параметров аккумулятора в процессе эксплуатации параметры модели могут быть обновлены повторением полного цикла заряда-разряда или его части. В последнем случае зависимость E0(SOC) не обновляется, а определение степени заряженности для узлов RS(SOC), t(SOC) происходит на основе измерения текущей ЭДС батареи.
Описание экспериментальной установки и методики проведения эксперимента
Структура и характеристики экспериментальной установки в полной мере представлены ранее [4]. При этом было модернизировано программное обеспечение микроконтроллера для обеспечения работы с моделями, учитывающими изменение температуры батареи.
Методика проведения эксперимента включает четыре этапа:
- вычисление параметров модели батареи при температуре +30 °С путем измерения ее реакции на импульсное воздействие током;
- определение параметров модели батареи при температуре 0 °С;
- разряд батареи номинальным током при температуре +10 °С, при этом точность оценки степени заряженности определяется по эталонному счетчику ампер-часов. Такой температурный режим выбран для проверки точности алгоритмов интерполяции параметров;
- сбор статистической информации в течение семи циклов заряда-разряда батереи для проверки точности методики.
Результаты эксперимента
Оценка параметров модели батареи при высокой температуре завершилась в соответствии с ожидаемыми значениями. Оценка ее параметров при низкой температуре показала, что параметры Я5, т невозможно интерполировать при помощи сплайн-функций, так как из-за больших перепадов и т между узлами интерполяции возникают броски в сплайн-функции.
Погрешность оценки остаточной емкости в начальный момент времени работы алгоритма незначительно превышает 10% (рис. 1, 2). По мере работы алгоритма погрешность снижается и остается на уровне 6%.
30 25 20 1 < t?,, А ч
-Ох
I J 10 5
■ьХ, "ЧА Ç
5000 10000 15000
Рис. 1. Зависимость остаточной емкости батареи от времени (сплошной линией показана оценка по методике, пунктирной - оценка по эталонному счетчику ампер-часов)
емкости батареи AQ^(t)
Для накопления статистики проведено несколько циклов заряда-разряда батареи, в ходе которых определялась погрешность оценки остаточной емко-
сти и степени заряженности. Было произведено семь циклов заряда-разряда: один - при температуре минус 10 °С; один - при +38 °С; два - при +30 °С; один - при 0 °С; два - при +10 °С.
Вся разрядная характеристика батареи подразделялась на три участка по степени заряженности, на каждом из которых была рассчитана погрешность оценки Q1(t), приведенная к шкале (см. таблицу).
муляторной батареи для применения на космическом аппарате / М.М. Хандорин, В.Г. Букреев // Авиакосмическое приборостроение. - 2017. - № 3. - С. 45-55.
5. Хандорин М.М. Оценка остаточной емкости литий-ионного аккумулятора в режиме реального времени / М.М. Хандорин., В.Г. Букреев // Электрохимическая энергетика. - 2014. - Т. 14, № 2. - С. 65-69.
6. Hongwen H. Evaluation of Lithium-Ion Battery Equivalent Circuit Models for State of Charge Estimation by an Experimental Approach / H. Hongwen, X. Rui, F. Jinxin // Energies. - 2011. - № 4. - P. 58 -598.
7. Saw L.H. Electro-thermal analysis of Lithium Iron Phosphate battery for electric vehicles / L.H. Saw, K. Soma-sundaram, Y. Ye et al. // Journal of Power Sources. -Vol. 249. - P. 231-238.
Погрешность оценки остаточной емкости Q1(t)
по участкам разрядной характеристики
Интер- Среднеквадратическое Среднеквадратическое
вал SOC отклонение, А-ч отклонение, %
1-0,8 1,6412 5,86
0,8-0,2 1,7074 6,09
0,2-0 1,5961 5,70
Заключение
В результате проведения экспериментов установлено, что погрешность оценки остаточной емкости батареи не превышает ±10% шкалы на всем протяжении разрядной характеристики. Таким образом, предложенная методика оценки остаточной емкости вполне работоспособна. Она позволяет учитывать текущий ток батареи и влияние температурного воздействия. Причины появления высокой погрешности оценки остаточной емкости в начальный момент времени являются предметом дальнейших исследований.
Литература
1. Романенко А.С. Энергобаланс в автономной системе электроснабжения с регулятором мощности на основе матричного преобразователя тока солнечной батареи / А.С. Романенко, А.В. Кобзев, В.Д. Семенов // Доклады ТУСУРа. - 2013. - № 4 (30). - C. 115-119.
2. Шиняков Ю.А. Энергетический анализ структурных схем систем электроснабжения автоматических космических аппаратов // Изв. Том. политехи. ун-та. - 2006. -Т. 309, № 8. - С. 152-155.
3. Романенко А. С. Модель автономной системы электроснабжения с секционированной солнечной батареей в пакете MatLab / А.С. Романенко, В. Д. Семенов // Доклады ТУСУРа. - 2011. - № 2 (24), ч. 1. - C. 269-274.
4. Хандорин М.М. Экспериментальная проверка алгоритма оценки остаточной емкости литий-ионной акку-
Хандорин Михаил Михайлович
Инженер-конструктор АО «НПЦ «Полюс» (Томск),
аспирант каф. ЭПЭО
Тел.: +7 (382-2) 55-46-94
Эл. почта: info@polus-tomsk.ru
Букреев Виктор Григорьевич
Д-р техн. наук, профессор каф. электропривода
и электрооборудования (ЭПЭО)
Энергетического института
Национального исследовательского
Томского политехнического университета (НИТПУ)
Тел.: +7-913-854-21-96
Эл. почта: bukreev@tpu.ru
Khandorin M.M., Bukreev V. G.
Method to estimate the available capacity of a Lithium-ion battery when the battery temperature changes
The article is devoted to the method used to evaluate the available capacity for spacecraft lithium-ion battery when the battery temperature changes. For the charge and discharge process description, when a battery temperature changes, the Tevenin model is used. The experimental test method result is presented.
Keywords: lithium-ion battery, Thevenin model, spline-interpolation, linear interpolation, microcontroller.
