Для решения этих задач была создана универсальная подпрограмма WOD, с обращения к которой начинались программы идентификации разработанные для моделей с одним гс — контуром поляризации, с двумя гс — контурами в форме (3.17) с чисто активной нагрузкой в форме (3.1) с учетом индуктивности цепи разряда, а также программы расчета обусловленности матриц наблюдений и погрешности идентификации. Программа написана для работы в диалоговом режиме и рассчитана на обработку информации по всем аккумуляторам, включенным в последовательную батарею. Исходные данные, вводимые с клавиатуры дисплея; количество аккумуляторов в батарее, массивы разрядного тока и напряжений на каждом аккумуляторе и всей батарее в целом, снятые в произвольных точках кривых iH (f) и Амн (/) (к — номер аккумулятора), число дискрет между точками вводимых массивов (вспомогательные массивы). Подпрограмма осуществляет компенсацию индуктивного выброса по описанной методике, рассчитывает нагрев нагрузки и среднее ее Г значение на время теста. На выходе подпрограммы формируется массив напряжений и тока с равномерным шагом дискретизации, характерным для работы с АП. Формирование производится в двух вариантах при линейной и экспоненциальной интерполяции между точками исходных массивов. Первоначальная обработка осциллограмм производилась вручную с помощью увеличительного проекционного аппарата, что уменьшало погрешности считывания. Результаты идентификации 23-х аккумуляторов 6ЭМ145 при различной степени заряженности 8 представлены в табл. 3.1. В той же таблице в графах 14, 15 даны значения косвенных параметров то и т£ которые предлагается использовать для оценки степени заряженности аккумуляторных батарей. Параметр т0 представляет собой следующее то ГФС«с Т.qc Га + Гп +Гф м3 [Ф] 112 |
компенсации индуктивного выброса, который заключается в следующем (рис.2.20). На зарегистрированных кривых Ли„(/) и iH(t) фиксируется точка после окончания переходного процесса. На рис. 2.23 она соответствует пятой дискрете. Вторая точка берется при удвоении времени (десятая дискрета). По измеренному значению ино вычисляется сопротивление нагрузки в 5-й и 10-й дискретах. '-,(5) = /»(5) г>(Ю) = Ц-°-Ац"(10) /„(10) Предположим, что сопротивление нагрузки за время от t=0 до t=10 увеличивается за счет нагрева примерно по линейному закону. Тогда, используя линейную экстраполяцию, находим г„(0) = 2г„(5)-г„(10); /„(0) = г„(0) Аналогично может быть представлен процесс обработки осциллограмм для ввода их в ПК. Необходимо было, во-первых, обеспечить максимальную точность при дискретном представлении сигналов и потому уменьшить ша дискретизации D (а значит и увеличивать размерность вводимых массивов); во-вторых, иметь возможность изменять шаг дискретизации внутри самого массива для более детальной проработки начальных участков, кривых, что обеспечивало необходимую точность идентификации; и в третьих, для обработки машинных программ реального диагностического устройства преобразовывать исходный массив к виду, характерному для работы с аналого-цифровыми преобразователями. Для решения этих задач была создана универсальная подпрограмма WOD, с обращения к которой начинались все программы идентификации, разработанные для моделей с одним гс контуром поляризации, с двумя гс контурами в форме (2.29) с чисто активной нагрузкой в форме (2.13) с учетом индуктивности цепи разряда, а также программы расчета обусловленности матриц наблюдений и погрешности идентификации. 107 Рис. 2.20. Осциллограммы кривых &uH(f) и Программа написана для работы в диалоговом режиме и рассчитана на обработку информации по всем аккумуляторам, включенным в последовательную батарею. Исходные данные, вводимые с клавиатуры дисплея; количество аккумуляторов в батарее, массивы разрядного тока и напряжений на каждом аккумуляторе и всей батарее в целом, снятые в 108 произвольных точках кривых /„(/) и AwH(Z) (к-номер аккумулятора), число дискрет между точками вводимых массивов (вспомогательные массивы). Подпрограмма осуществляет компенсацию индуктивного выброса по описанной методике, рассчитывает нагрев нагрузки и среднее ее значение на время теста. На выходе подпрограммы формируется массив напряжений и тока с равномерным шагом дискретизации, характерным для работы с АП. Формирование производится в двух вариантах при линейной и экспоненциальной интерполяции между точками исходных массивов. Первоначальная обработка осциллограмм производилась вручную с помощью увеличительного проекционного аппарата, что уменьшало погрешности считывания. В экспериментальных исследованиях были использованы ПК и диалоговый вычислительный комплекс ДВК-2. Результаты идентификации 23-х аккумуляторов 6ЭМ145 при различной степени заряженности б представлены в табл. 2.7. В той же таблице в графах 14,15 даны значения косвенных параметров т0 и те которые предлагается использовать для оценки степени заряженности аккумуляторных батарей. Параметр т0 представляет собой следующее т _ гФсдс _ ^дс Г„+Г„+Гф мз где Т = гф • cqc постоянная времени одиночного контура актива активационной поляризации. Как показал анализ, наблюдается достаточно хорошая корреляция между начальной емкостью аккумуляторной батареей Отах (или „ S бтах. «номинальной») степенью заряженности е =---------------) и параметром т0. Quito Коэффициент корреляции близок к единице (0,989). На рис.2.21 показаны экспериментальные точки этой зависимости и прямая регрессия вида: т0 = -0,79 + 5,36s0 прямая регрессия, (2.69) 109 |