CC BY
268
УДК 541.136:519.876.5
ОЦЕНКА ОСТАТОЧНОЙ ЁМКОСТИ ЛИТИЙ-ИОННОГО АККУМУЛЯТОРА В РЕЖИМЕ РЕАЛЬНОГО ВРЕМЕНИ
М. М. Хандорин, В. Г. Букреев1 и
ОАО «Научно-производственный центр «Полюс» 634050, Россия, Томск, пр. Кирова, 56в 1 Национальный исследовательский Томский политехнический университет 634050, Россия, Томск, пр. Ленина, 30
и E-mail: polus@online.tomsk.net, bukreev@tpu.ru Поступила в редакцию 30.07.14 г.
Предложены алгоритмы оценки остаточной и полной ёмкости литий-ионного аккумулятора на основе модели Тевенина. Рассмотрен вариант определения параметров модели с применением сплайн-интерполяции дискретных экспериментальных данных. Выполнен анализ возможности реализации предложенных алгоритмов в приложениях реального времени на микроконтроллере.
Ключевые слова: литий-ионный аккумулятор, модель Тевенина, интегрирующий счётчик ампер-часов, сплайн-интерполяция, микроконтроллер.
LITHIUM-ION BATTERY CAPACITY ESTIMATION AT THE REAL-TIME
M. M. Khandorin, V. G. Bukreev1 и
Joint Stock Company «Research and Production Center «Polus», 56 v, Kirov ave., Tomsk, 634050, Russia 1 National Research Tomsk Polytechnic University, 30, Lenin ave., Tomsk, 634050, Russia
и E-mail: polus@online.tomsk.net, bukreev@tpu.ru Received 30.07.14
At this paper algorithms of accessible and full capacity estimation are presented. The model discrete parameter interpolation, based on spline function is shown. Presented algorithms are analyzed for use in real-time applications.
Key words: lithium-ion battery, integrating ampere-hour counter, Thevenin model, spline-interpolation, microcontroller, Mathlab, Simulink.
ВВЕДЕНИЕ
Важной характеристикой литий-ионного аккумулятора (ЛИА) является его остаточная ёмкость, для определения которой часто используется программная или аппаратная реализация счётчика ампер-часов [1]. Как правило, алгоритм работы такого счётчика включает установку начального значения остаточной ёмкости при условии достижения полной заряженности и фиксированного напряжения аккумулятора. При этом предполагается, что первоначально аккумулятор имеет остаточную ёмкость, указанную в технической документации, и один из алгоритмов определения полной ёмкости аккумулятора использует момент его полного разряда [2].
Несмотря на то что программная реализация счётчика ампер-часов не требует значительных вы-
числительных ресурсов, при его эксплуатации в составе измерительной системы оценки состояния аккумулятора выявлен ряд существенных недостатков. Во-первых, аккумулятор не всегда заряжается до максимального значения напряжения, что обусловливает невозможность определения момента начальной установки счётчика ампер-часов и, как следствие, накопление ошибки. Во-вторых, счётчик ампер-часов не учитывает часть ёмкости, которую теряет ЛИА в составе батареи из-за саморазряда при хранении и вследствие протекания балансировочных токов аккумуляторов. Те же недостатки проявляются и при использовании счётчика ампер-часов для оценки полной ёмкости ЛИА. Кроме того, остаточная ёмкость аккумулятора экспоненциально зависит от тока нагрузки согласно закону Пеккерта [1]. Невозможность точного учёта этого эффекта приводит к значительным ошибкам определения остаточ-
© ХАНДОРИН М. М., БУКРЕЕВ В. Г., 2014
ной ёмкости, особенно при больших токах разряда аккумулятора.
Определение ёмкости аккумулятора с помощью счётчика ампер-часов используется для лабораторных испытаний ЛИА, а для работы в составе систем электропитания космических аппаратов чаще применяются другие методы, основанные на математических моделях аккумулятора (например, фильтр Калмана, структуры на основе нечёткой логики и нейронных сетей) и обладающие гораздо большей точностью [1, 2].
Однако сдерживающим фактором применения таких методов для систем измерения остаточной ёмкости аккумулятора в реальном времени является необходимость значительных вычислительных ресурсов [1, 3].
Исходя из этого предлагается вычислительный алгоритм, позволяющий на основе степени заряжен-ности ЛИА оценить в реальном времени остаточную и полную ёмкость аккумулятора. Данный алгоритм также исключает недостатки, присущие классическому определению этих параметров с использованием счётчика ампер-часов.
ПОСТРОЕНИЕ МОДЕЛИ ЛИТИЙ-ИОННОГО АККУМУЛЯТОРА НА ОСНОВЕ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ
Поведение ЛИА на малом отрезке времени достаточно адекватно описывается математической моделью Тевенина [3, 4]:
ип(0 = - UC + Сп 'лиаСО;
илиА(0 = Eo - Un(t) - RoimA(t),
(1)
где ип(г) - поляризационное напряжение ЛИА; Яп, Сп, Ео, Яо - поляризационное сопротивление, поляризационная ёмкость, ЭДС, активное сопротивление соответственно; г'шдХО - ток разряда-заряда аккумулятора; иЛИА(г) - напряжение на выходе аккумулятора.
Малый отрезок времени используется для того, чтобы пренебречь влиянием степени заряжен-ности на параметры модели (1). При описании работы аккумулятора на больших интервалах времени необходимо учитывать функциональные зависимости параметров Ео(БОС), Яо(БОС), Яп(БОС), Сп(Б ОС) от степени заряженности Б ОС аккумулятора [3, 4]. Эти функции могут быть экспериментально определены в дискретных точках разрядной характеристики с промежуточной интерполяцией значений [2, 4]. Иногда [1, 4] для этой цели применяется кусочно-линейная аппроксимация, однако для уменьшения объёма экспериментальных данных це-
(2)
лесообразно интерполировать дискретные значения E0i (SOCi), Ro' (SOCi), Rni (SOCi), (SOCi) кубическим сплайном. На каждом интервале между двумя интерполяционными узлами [S OC—, S OCi] промежуточные значения параметров вычисляются по формулам [5]:
Cn (S OC) = aCi + bQ (S OC) + Cf (S OC)2 +
+ Ц- (SOC)3; Rn (S OC) = aRi + bRi (SOC) + CR (S OC)2 +
+ djt (SOC)3; E0i (S OC) = aEi + bEi (S OC) + C-f- (S OC)2+ + % (SOC)3;
CR0. т
R0i (S OC) = aR0. + bR0. (S OC) + ^ (S OC)2 +
dR0■ о
+ -£ (SOC)3,
где aCi, bCi, CCi, dCi, aR0., bRoi, CR0i, dR0i ,aEi, bEi, CEi, dEi - коэффициенты сплайна, определяемые предварительно в диапазоне изменения характеристик или в режиме реального времени.
Таким образом, для описания функцийE0(S OC), R0(SOC), Ra(S OC), Ca(S OC) получаем семейство сплайнов, определяемых степенью заряженности SOC. Вычисление значений параметров по дифференцируемой сплайн-функции позволяет более точно решить дифференциальное уравнение (1).
Так как начальный и конечный участки разрядной характеристики c экспоненциальной зависимостью параметров модели от остаточной ёмкости занимают около 20% ёмкости аккумулятора [6], то для обеспечения адекватности модели (1) параметры^ (SOCi), Rot (SOCi), Rm (SOCi), CUi (S OCi) определяются не менее чем в 20 опорных точках. Такое количество экспериментальных данных позволяет применить равномерную сетку дискретных значений S OCi при расчёте коэффициентов сплайна, обеспечивающих необходимую точность определения параметров аккумулятора [6]. Следует отметить, что значения параметров Ra(SOC), Cn(SOC) в модели Тевенина для разряда и заряда аккумулятора будут отличаться при одинаковых значениях S OC, поэтому их описание будет требовать, как минимум, двух семейств сплайнов.
Для имитационного моделирования переходной характеристики аккумулятора формируется импульс тока ¿лиа(0 разряда (нагрузки), а затем равный по значению импульс тока заряда. Далее измеряются значения напряжения Цлиа(0, которые используются для определения параметров модели Те-венина.
Расчёт напряжения Um_A(t) аккумулятора производится путём решения системы уравнений (1) ме-
тодом Рунге-Кутта. Обозначив результат этих вычислений функцией ¥гк(Яп, Сп, Ео,Яо, /диа!, к), запишем среднеквадратическое отклонение:
п 2
г = ^ (илИА; - Рук(Яп, Сп,Ео,Яо, /ЛИА;, , (3)
I
где илиа; - значение напряжения аккумулятора в дискретных точках.
Определение значений параметров модели в конечном итоге сводится к оптимизации многомерной функции суммы квадратов отклонений. Для оптимизации многомерной функции, включающей изменяемые параметры модели аккумулятора, целесообразно использовать алгоритм координатного спуска (рис. 1) [7].
Рис. 1. Алгоритм определения параметров модели Тевенина
Алгоритм реализует подбор параметров модели за несколько итераций, каждая из которых состоит из четырёх этапов. На каждом из них последовательно оптимизируется один из параметров - Ео, Яо,
Сп, Яп, при этом остальные параметры фиксируются.
Алгоритм подстройки параметров модели Теве-нина под экспериментальный аккумулятор предполагает остановку при уменьшении темпа сокращения среднеквадратической ошибки. Как только изменение среднеквадратической ошибки между последующей итерацией тп+\ и текущей гп оказывается меньше порогового значения, равного 1 ■ 10-1о, то вычислительный процесс завершается. В результате моделирования в среде МаАаЬ Simulink установлено, что такая процедура позволяет достичь суммарной среднеквадратической ошибки г <0.1 (при использовании встроенной модели ЛИА в среде МаЙаЬ для имитации реальной батареи).
Функция (3) может иметь несколько экстремумов, и для сходимости алгоритма к области глобального минимума необходимо определить направление поиска, которое обеспечит его нахождение с приемлемой точностью. Для ускорения оптимизации предлагается использовать значения координат предыдущего интерполяционного узла в качестве начальной точки поиска.
ОЦЕНКА СТЕПЕНИ ЗАРЯЖЕННОСТИ ЛИА
Модель Тевенина достаточно адекватно описывает процессы в ЛИА только для некоторого известного значения степени заряженности, и очевидно, что по мере протекания разрядно-зарядного тока степень заряженности аккумулятора будет изменяться. Данные изменения будут приводить к появлению разности е(0 между напряжением Цлиа(0 на выходе аккумулятора и оценкой Удиа (0 этого напряжения по модели Тевенина:
e(t) = ^лиа(t) - ^лид(г).
(4)
Так как напряжение на выходе аккумулятора Цпиа(0 будет связано со степенью заряженности SOC через параметры модели Тевенина, соотношение (4) можно использовать для их определения.
Введём поправочный коэффициент для коррекции степени заряженности ЛИА по мере накопления ошибки по напряжению. Таким образом, степень за-ряженности в любой момент времени будет описываться соотношением
SOCi = STOOCo + ^ ô j (e),
j=o
(5)
где S ОСо - начальная степень заряженности, Ь^ (е) -поправочный коэффициент в предыдущие моменты времени.
Если степень заряженности Б ОС в модели ЛИА не превышает реальную Б ОС, т. е. БОС < БОС, то ошибка принимает отрицательное значение: е(0 < 0, при БОС > БОС имеем положительные значение: е(г) > 0. Таким образом, используя поправочный коэффициент 6(8) при вычислении значений^ ОС, можно получить примерное равенство: БОС~БОС. Примем пропорциональную зависимость между поправочным коэффициентом и ошибкой по напряжению иЛИА(0 аккумулятора:
6(8) = кг(г), (6)
где к - нормировочный множитель.
ОЦЕНКА ОСТАТОЧНОЙ ЁМКОСТИ ЛИТИЙ-ИОННОГО АККУМУЛЯТОРА
Для оценки остаточной ёмкости ЛИА 0ост необходимо знать его полную ёмкость и сте-
пень заряженности. При этом остаточная ёмкость может быть вычислена по уравнению
бост — QmaxS ОС.
(7)
Полная ёмкость ЛИА, как правило, известна в начальный момент его использования и по мере эксплуатации будет снижаться. В реальной системе электропитания происходят неполные циклы разряда-заряда аккумулятора, поэтому построение алгоритма определения полной ёмкости при частичном заряде и разряде представляет собой важную задачу. Так, алгоритм Impedance Track [8] позволяет определить полную ёмкость ЛИА по линейному участку разрядной характеристики. Исходными данными для работы этого алгоритма является график зависимости степени заряженности аккумулятора SOC от напряжения E0 без нагрузки. Используя эту характеристику, можно вычислить S OC следующим образом.
Предположим, что в некоторый начальный момент времени известны значение ёмкости Qi и соответствующая ей степень заряженности S OC\. Далее при разряде аккумулятора счётчик ампер-часов подсчитывает изменение ёмкости аккумулятора (Qi — — Q2) на интервале времени (t2 — ti). Уменьшение остаточной ёмкости ЛИА приведёт к снижению степени заряженности до значения SOC2. Предполагая линейность характеристик ЛИА на интервале (t2 — —1\), полную электрохимическую ёмкость аккумулятора можно определить по выражению
Qmax —
(Qi - Q2)
(8)
(БОСг - БОС2).
Наибольшая точность вычисления Qmж будет обеспечиваться на линейном участке разрядной харак-
теристики при изменении относительных значений S ОС в пределах от 0.6 до 0.4.
Таким образом, объединяя модель Тевенина, алгоритмы определения степени заряженности и полной ёмкости ЛИА, можно вычислить остаточную ёмкость аккумулятора (рис. 2).
Рис. 2. Блок-схема алгоритма оценки остаточной ёмкости ЛИА
МОДЕЛИРОВАНИЕ ПОСТРОЕННЫХ АЛГОРИТМОВ В СРЕДЕ МА1ЪАВ SIMULINK
Для проверки работоспособности алгоритма оценки остаточной ёмкости было проведено ими-
тационное моделирование в среде МаИаЪ Simulink, где в качестве эквивалента реального аккумулятора использовалась модель батареи ЛИГП-10 ТУ 3482-063-20503890-2006. Параметры модели приведены в табл. 1.
Таблица 1
Параметры модели батареи ЛИГП-10
Параметр Значение
Номинальное напряжение, В 3.6
Полная ёмкость, А-ч 10
Максимальная ёмкость, А-ч 10.35
Напряжение при полном заряде, В 3.9
Номинальный ток разряда, А 2
Внутренне сопротивление, Ом 0.02
Ёмкость при номинальном напряжении, А-ч 4.5
Граница экспоненциальной зоны по напряжению, В 3.8
Граница экспоненциальной зоны по ёмкости, А-ч 0.5
Вычислительный алгоритм решения дифференциальных уравнений для модели Тевенина выполнен на стандартных элементах среды Matlab, прочие алгоритмы моделировались логическими блоками. Параметры модели определялись на основе экспериментальных данных при помощи специальной программы, реализующей оптимизационный алгоритм координатного спуска. Значения этих параметров при различной степени заряженности ЛИА приведены в табл. 2.
Первый численный эксперимент был поставлен для проверки адекватности оценки значений степени заряженности ЛИА при различных начальных условиях SOCo (рис. 3).
В результате исследований установлено, что алгоритм обеспечивает погрешность определения SOC на приемлемом уровне, при этом точность практически не зависит от начальных условий.
Для проверки алгоритма оценки полной ёмкости аккумулятора на основе изменения степени за-ряженности и показаний счётчика ампер-часов была построена его логическая модель в среде Matlab Simulink. При этом используется имитационная модель аккумулятора и модель Тевенина из предыдущего численного эксперимента.
100
50
-Истинное значение S OC
---Оценка S OC при S OC0 = 50%
.......Оценка S OC при S OC0 = 0%
100
200
300
400
500
Рис. 3. Определение степени заряженности ЛИА
На рис. 4 представлен график работы алгоритма при определении значения полной ёмкости ЛИА, реальное значение которой равно 10.35 А ч.
F
<
и 10.4
10.2
10.0
9.8
—|-1-1-1-1-1-1-1-1-1—
Истиное значениезначение Сщ^ = Ю.35 А-ч
|_
i__
Момент определения Сщах при другом-токе разряда
Начальное значение Сщах — Ю
_i_
Момент определения Сщ^
_L
_L
_L
_L
_L
5000 10000 15000 20000 25000
t, с
Рис. 4. Определение значения полной ёмкости ЛИА
Таблица 2
Значения параметров модели Тевенина
0
0
t, с
SOC E0, В R0, Ом Зарядное сопротивление Rth, Ом Разрядное сопротивление Rth, Ом Зарядная ёмкость Cth, Ф Разрядная ёмкость Cth, Ф
0.2 2.96227 0.01216 0.03414 0.12636 344.81133 94.77973
0.3 3.35284 0.01480 0.03390 0.08289 327.88814 139.40513
0.4 3.54723 0.01560 0.03680 0.06201 291.63949 180.17834
0.5 3.66364 0.01577 0.04175 0.04985 251.40727 215.55111
0.6 3.74112 0.01572 0.04913 0.04189 211.68956 245.26821
0.7 3.79638 0.01545 0.06050 0.03635 172.10408 267.24373
0.8 3.83778 0.01489 0.07964 0.03247 131.81754 277.53342
0.9 3.86981 0.01348 0.11822 0.03011 89.88268 262.67406
1.0 3.93204 0.01000 0.23429 0.03085 47.65809 250.45772
В первый момент времени полная ёмкость будет равна произвольному начальному значению. Это продолжается до тех пор, пока аккумулятор не разрядится до 60% (до начала линейного участка), после начинается работа алгоритма определения Стах. К моменту достижения степени заряженности ЛИА 40% работа алгоритма завершается. При следующем определении Стах алгоритм может сформировать значение, отличающееся от предыдущего на 0.1-0.5 А ч. Это связано с тем, что разрядная характеристика аккумулятора имеет небольшую нелинейность, следовательно, изменение тока разряда будет отражаться на результате.
С использованием имитационной модели в среде МайаЬ и модели Тевенина был проведён численный эксперимент для оценки точности алгоритма определения остаточной ёмкости, структура которого была построена из логических блоков среды МайаЬ. Входным сигналом для имитационной модели ЛИА является ток разряда или заряда, изменяющийся случайным образом один раз за 10 с. При этом оценка остаточной ёмкости аккумулятора обеспечивается с точностью не хуже 2.5% (рис. 5).
г, с
Рис. 5. Сравнение оценочного значения остаточной ёмкости ЛИА (сплошная линия) с действительным (штриховая линия)
В первоначальный момент времени работы алгоритма (примерно 500 с) погрешность весьма велика, так как выполняется приближение расчётного значения степени заряженности к действительному значению. Далее, до момента времени 6000 с, погрешность остаётся также значительной из-за неточного определения полной ёмкости аккумулятора. После завершения работы алгоритма определения полной ёмкости аккумулятора погрешность вычисления остаточной ёмкости сокращается до минимального значения и составляет примерно 2.5%.
ОЦЕНКА ВОЗМОЖНОСТИ РЕАЛИЗАЦИИ АЛГОРИТМА НА МИКРОКОНТРОЛЛЕРЕ В РЕАЛЬНОМ ВРЕМЕНИ
Так как большинство микроконтроллеров обладает ограниченными вычислительными ресурсами, представляется важным вопрос об объёме памяти и процессорном времени, необходимых для выполнения алгоритма.
Объём памяти для хранения коэффициентов сплайнов при 16-битной разрядности каждого коэффициента определяется исходя из следующего: четыре коэффициента в каждом сплайне, четыре параметра модели Тевенина, десять узлов сплайна, два режима работы (разряд и заряд), что составляет 320 байт. Хранение промежуточных данных при решении дифференциальных уравнений методом Рун-ге-Кутта потребует 16 байт ОЗУ. С учётом необходимого резерва вычислительных алгоритмов ответственных микропроцессорных систем общий объём оперативной памяти не будет превышать 320 байт. Затраты процессорного времени для одной итерации алгоритма оценки остаточной ёмкости ЛИА и метода Рунге-Кутта на языке Си с использованием 16-разрядной целочисленной арифметики будут составлять 10 мс.
ЗАКЛЮЧЕНИЕ
Работоспособность и эффективность предлагаемой методики оценки остаточной ёмкости ЛИА доказана путём численных экспериментов в среде МаНаЬ Simulink на имитационной модели аккумулятора ЛИГП-10. Алгоритм позволяет оценивать степень заряженности ЛИА с погрешностью не более 2%, а остаточную ёмкость - не более 2.5%.
По сравнению с алгоритмом, основанным на фильтре Калмана [3], предложенный алгоритм требует значительно меньших вычислительных затрат, точность определения с его помощью степени за-ряженности ЛИА не зависит от задания начальных условий.
Кроме того, разработанный алгоритм позволяет оценивать полную электрохимическую ёмкость ЛИА в процессе эксплуатации, что открывает путь не только для предупреждения аварийных ситуаций, но и для управления ресурсом ЛИА.
Высокая точность и низкие требования к аппаратно-программным ресурсам микроконтроллера делают полученную методику оценки остаточной ёмкости ЛИА перспективной для внедрения в аппаратуре различного назначения.
СПИСОК ЛИТЕРАТУРЫ
1. Chang W. Y. The state of charge estimating methods for battery // Hindawi Publishing Corporation ISRN Applied Mathematics. 2013. № 1. P. 1-7.
2. Ng K. S, Moo C. S. Chen Y. P., Hsieh Y. C. Enhanced Coulomb counting method for estimating state-of-charge and state-of-health of lithium-ion batteries. // Applied Energy. 2009. Vol. 86, № 9. P. 1506-1511.
3. Hongwen H., Rui X., Jinxin F. Evaluation of Lithium-Ion Battery Equivalent Circuit Models for State of Charge Estimation by an Experimental Approach. // Energies. 2011. № 4. P. 582-598.
4. Rahmoun A., Biechl H. Modelling of Li-ion batteries using equivalent circuit diagrams. // Przegl^d elektrotechniczny. 2012. № 7b. P. 152 - 156.
5. Завьялов Ю. С., Квасов Б. И., Мирошниченко В. Л. Методы сплайн-функций. М. : Наука, 1980.
6. Таганова А. А., Бубнов Ю. И., Орлов С. Б. Герметичные химические источники тока. Элементы и аккумуляторы, оборудование для испытаний и эксплуатации. СПб. : Химиздат, 2005.
7. Bunday B. D. Basis optimisation methods. London : Edward Arnold Publ., 1984.
8. Wen S. Impedance Track Gas Gauge for Novices // Texas Instruments Application Report [офиц. сайт]. URL: http://www.ti.com/lit/an/slua375/slua375.pdf (дата обращения 10.06.2014).
REFERENCES
1. Chang W. Y. The state of charge estimating methods for battery. Hindawi Publishing Corporation ISRN. Applied Mathematics. 2013, vol. 1, pp. 1-7.
2. Ng K. S., Moo C. S., Chen Y. P., Hsieh Y. C. Enhanced Coulomb counting method for estimating state-of-charge and state-of-health of lithium-ion batteries. Applied Energy. 2009, vol. 86, no. 9, pp. 1506-1511.
3. Hongwen H., Rui X., Jinxin F. Evaluation of lithium-ion battery equivalent circuit models for state of charge estimation by an experimental approach. Energies, 2011, no. 4, pp. 582-598.
4. Rahmoun A., Biechl H. Modelling of Li-ion batteries using equivalent circuit diagrams. Przeglqd elektrotechniczny [Electrical Review]. 2012, no. 7b, pp. 152 - 156.
5. Zavyalov Yu. S., Kvasov B. I., Miroshnichenko V. L. Metodysplayn-funktsiy. [Methods of spline functions]. Moscow, Nauka Publ., 1980 (in Russian).
6. Taganova A. A., Bubnov Yu. I., Orlov S. B. Germetichnyye khimicheskiye istochniki toka. Elementy i akkumulyatory, oborudovaniye dlya ispytaniy i ekspluatatsii. [Sealed power sources. Elements and battery testing equipment and operation]. S. Peterburg, Khimizdat Publ., 2005 (in Russian).
7. Bunday B.D. Basis optimisation methods. London. Edward Arnold Publ., 1984.
8. Wen S. Impedance Track Gas Gauge for Novices, Texas Instruments Application Report. Avialable at: http://www.ti.com/lit/ an/slua375/slua375.pdf (accessed 10 June 2014).
СВЕДЕНИЯ ОБ АВТОРАХ
Хандорин Михаил Михайлович - инженер-конструктор, ОАО «НПЦ «Полюс». Служебный телефон 8 (382-2) 55-46-94, e-mail: polus@online.tomsk.net
Букреев Виктор Григорьевич - д-р техн. наук, профессор, Национальный исследовательский Томский политехнический университет, кафедра «ЭПЭО». Служебный телефон +79138542196, e-mail: bukreev@tpu.ru
