Двухфазный повышающий конвертер как звено системы автоматического управления
Импульсные преобразователи повышающего типа необходимы в тех случаях, когда первичный источник электропитания (аккумулятор, солнечная батарея, генератор постоянного напряжения) не обеспечивает необходимое напряжение на нагрузке.
При малых уровнях потребляемой мощности достаточно использовать однотактный повышающий конвертер (рис. 1), в котором на интервале включенного состояния транзистора γ = tи/T энергия запасается в дросселе, а при выключении транзистора она через диод VD передается в конденсатор С и на нагрузку. В результате напряжение на нагрузке становится выше входного напряжения [1, 2].
Если коммутируемые токи превышают допустимый ток одного транзистора, то можно использовать параллельное включение приборов. Однако при этом возникают проблемы не только одновременного включения, но и подбора транзисторов по прямому падению напряжения. Для минимизации этих проблем используется многофазный режим работы преобразователей [1–3].
Двухфазный повышающий конвертер — это два повышающих импульсных преобразователя, питающихся от одного источника и работающих на общую нагрузку с временным сдвигом на полпериода (рис. 2).
Основной задачей преобразователей (рис. 1, 2) является обеспечение стабильного выходного напряжения при изменении в широких пределах входного напряжения и выходного тока. Стабилизация выходного напряжения реализуется введением обратных компенсационных связей. Для определения параметров цепей обратных связей при заданных статических и динамических характеристиках необходимо составить расчетную схему и определить передаточные функции вход/выход и управление/выход.
Передаточные функции однотактных конвертеров на основе усредненных моделей получены в работах [1, 4 и 5]. В работе [2] были определены динамические параметры двухфазного повышающего импульсного преобразователя, представленного структурной схемой, приведенной на рис. 3.
На схеме символом D обозначены малые отклонения величин Uвх.ср(t), iвх.ср(t), uн.ср(t), iн.ср(t), γ(t) от их стационарных значений. В соответствии со структурной схемой, приведенной в работе [2], была получена следующая передаточная функция по управляющему воздействию:
где G(p) = 2/(Lp + r); Z(p) = 1/(pC).
При определении передаточной функции Wγ(p) в каноническом виде в работе [2] не было учтено сопротивление нагрузки Rн и эквивалентное последовательное сопротивление выходного конденсатора rC. Поскольку повышающий конвертер не должен работать в режиме холостого хода, то необходимо учитывать Rн. С учетом Rн и допущением того, что
соотношения для Wγ(p) будут следующими:
где
Если же учесть еще и сопротивление потерь rC выходного конденсатора C, то
и можно получить более точную передаточную функцию по входному воздействию:
где
В отличие от формулы (1), передаточная функция (2) имеет ноль в левой комплексной полуплоскости, обусловленный сопротивлением rC потерь конденсатора C.
Проверим достоверность выражений (1) и (2) и принятых при их выводе допущений, используя альтернативный способ определения динамических параметров. Например, передаточная функция может быть составлена на основании переходной характеристики (кривой разгона) [6].
Рассмотрим применение этого метода для определения динамических параметров двухфазного повышающего преобразователя. Для получения кривой разгона будем использовать компьютерное моделирование в Simulink MATLAB (рис. 4) для следующих исходных данных:
L = 60 мкГн; С = 410 мкФ; r = 0,01 Ом; rC = 0,003 Ом; Rн = 1,9 Ом; U = 175 B; γ = 0,253; f = 20 кГц.
После выхода преобразователя на установившийся режим работы система управления (блок Subsystem на рис. 4) скачком изменяла коэффициент заполнения g, и на осциллографе ScopeURH фиксировался переходный процесс изменения напряжения нагрузки. Нормированная кривая разгона напряжения нагрузки представлена на рис. 5.
На первом периоде работы преобразователя после действия скачка напряжения управления на увеличение длительности включенного состояния tи транзистора напряжение нагрузки сначала уменьшается, так как длительность выключенного состояния транзистора, при котором энергия, накопленная в дросселе, передается в нагрузку, становится меньше, чем в предыдущем периоде работы транзисторов. Эта особенность работы повышающего преобразователя хорошо видна на кривой разгона, соответствующей переходной характеристике неминимально-фазового звена.
Неминимально-фазовый четырехполюсник можно заменить двумя последовательно соединенными четырехполюсниками, один из которых является минимально-фазовым, а другой — четырехполюсником чисто фазового сдвига [6]. Длительность интервала, когда напряжение на нагрузке меньше предыдущего значения напряжения, составляет всего один период (рис. 5), поэтому можно аппроксимировать кривую разгона переходной характеристикой колебательного звена, а этот интервал учесть временем запаздывания tз.
Передаточная функция колебательного звена имеет следующий вид:
Динамическими параметрами колебательного звена являются постоянная времени Т — величина, обратная собственной частоте колебательного звена, и относительный коэффициент затухания ξ. Эти параметры по кривой разгона можно определить несколькими методами (рис. 6).
В первом методе параметры Т и x определяются по выражениям [7]
Во втором методе [8]
где
В третьем методе параметры колебательного звена определяются по амплитудно-фазовой характеристике (АФХ) по следующим простым соотношениям:
где ωТ — частота, при которой АФХ пересекает мнимую ось.
Для построения АФХ кривую разгона аппроксимируем отрезками прямых (рис. 6 б) [9]. Затем определяем коэффициенты ci, i = 1,n, характеризующие углы наклона n прямолинейных отрезков:
c1 = h1/t1; c2 = (h2 – h1) / (t2 – t1); …; cn = (hn – hn–1) / (tn – tn–1).
Значения вещественной U и мнимой V частей частотной передаточной функции W(jw) при каждой частоте ω рассчитываем по формулам
где n — число аппроксимирующих отрезков;
τ1i = (ti – ti-1) / 2 и τ2i = (ti + ti-1) / 2
На рис. 7 показана АФХ двухфазного повышающего преобразователя.
На рис. 8 приведены нормированная кривая разгона напряжения нагрузки и аппроксимирующие кривые колебательного звена второго порядка. Все три метода дают близкие результаты.
В первом методе: t1 = 2,73×10–4 c, t3 = 7,397×10–4 c и T = 1,4329×10–4 c, ξ = 0,2638.
Второй метод дает A1 = 0,697 B, A2 = 0,321 B и по формулам (4) получаем T = 1,4737×10–4 c, ξ = 0,1225.
В третьем методе по АФХ (рис. 7) определяем ωт = 6,6596×103 рад/с и V(ω0) = –4,2169. По соотношениям (5) вычисляем T = 1,5016×10–4 c, ξ = 0,1186.
Из рис. 8 следует, что для двухфазного повышающего преобразователя наименьшее отклонение аппроксимирующей функции от экспериментальной кривой разгона дает третий метод, поэтому преобразователь может быть представлен по напряжению нагрузки колебательным звеном с передаточной функцией
и звеном чистого запаздывания e–pt3.
В таблице приведены результаты расчета динамических параметров двухфазного повышающего преобразователя, полученные различными методами.
Метод определения параметров |
T, c |
ξ |
Расчет по формуле (1) без учета rC |
0,14812×10–3 |
0,1074 |
Расчет по формуле (2) c учетом rC |
0,14824×10–3 |
0,1115 |
Моделирование по формуле (5) |
0,1502×10–3 |
0,1186 |
В таблице небольшие расхождения результатов, полученных по формулам (2) и (5), связаны с тем, что в формуле (2) не учитывались параметры ключевых приборов.
На рис. 9 б показаны результаты вычислений откликов (рис. 9 а) при подаче единичного скачка на входы передаточных функций, полученных по формулам (1) и (2), в сравнении с экспериментальной кривой разгона двухфазного повышающего преобразователя. Наилучшее приближение к экспериментальной кривой разгона дает передаточная функция с учетом сопротивления потерь конденсатора. Но это улучшение незначительно, и, следовательно, для принятых исходных данных значением rC можно пренебречь, поскольку учет его усложняет модель.
На рис. 10 приведены логарифмические амплитудно- и фазочастотные характеристики двухфазного повышающего преобразователя, построенные по кривой разгона напряжения нагрузки с использованием передаточной функции (2) и передаточной функции колебательного звена с параметрами, определенными по формулам (5).
Полученные частотные характеристики практически совпали. Это свидетельствует о том, что для приведенных исходных данных динамические параметры и передаточные функции двухфазного повышающего конвертера по напряжению нагрузки могут быть определены как по формуле (2), так по более простому выражению (6).
Динамические параметры могут быть определены и по току дросселей, а поскольку в двухфазном повышающем преобразователе два дросселя, то удобнее определять передаточную функцию по входному току id. После выхода преобразователя на установившийся режим работы система управления изменяла коэффициент заполнения g, и на осциллографе ScopeId (рис. 4) фиксировался переходный процесс изменения тока, потребляемого от источника питания. Нормированная кривая разгона входного тока id представлена на рис. 11, а АФХ, построенная по кривой разгона, приведена на рис. 12. Логарифмические АЧХ и ФЧХ конвертера показаны на рис. 13.
Из приведенных характеристик следует, что передаточная функция по току id наряду с колебательным должна содержать еще и форсирующее звено
Этот факт подтверждается и результатами аппроксимации кривой разгона колебательным звеном (рис. 14). Как и в случае с кривой разгона по напряжению нагрузки, параметры колебательного звена определялись тремя методами по формулам (3–5). Следует отметить, что применение первого метода в данном случае неприемлемо, так как коэффициент затухания оказывается меньше нуля.
Параметры T и ξ колебательного звена принимаем равными полученным по второму методу, поскольку длительность периода и отношение A1/A2 совпадают с аналогичными параметрами кривой разгона h(t):
T = 1,4752×10–4 c и ξ = 0,1223.
Для определения постоянной времени форсирующего звена в качестве первого приближения используем значение производной тока dh/dt = 2,1×104 A/c в начале кривой h(t) (рис. 15).
После подстановки в передаточную функцию (7) принятых значений T, ξ и начального значения t = dh/dt × T2 = 4,57×10–4 выполним проверку правильности значения постоянной времени форсирующего звена — подадим единичный скачок на функцию (7) и сравним отклик с экспериментальной кривой разгона. Из рис. 16 следует, что имеет место практически полное совпадение модели (синяя кривая) с экспериментальной кривой h(t) (красный пунктир) и нет необходимости выполнять итерационную процедуру уточнения постоянной времени форсирующего звена.
Для проверки правильности определения значения постоянной времени τ форсирующего звена можно использовать альтернативный алгоритм.
Передаточная функция системы из K последовательно соединенных звеньев равна произведению передаточных функций отдельных звеньев системы:
АЧХ последовательного соединения звеньев равна произведению АЧХ отдельных звеньев, а ФЧХ — сумме ФЧХ соединенных звеньев:
При последовательном соединении звеньев с известными АФХ результирующая АФХ может быть получена путем сложения углов, под которыми направлены радиусы-векторы АФХ отдельных звеньев при заданной частоте, а результирующий радиус-вектор будет равен произведению модулей звеньев на этой же частоте. Поскольку двухфазный повышающий преобразователь по входному току id может быть представлен в виде последовательного соединения двух звеньев — колебательного и форсирующего, то его передаточная функция (7) может быть записана в виде произведения
W(p) = Wк(p)Wф(p),
где Wк(p) и Wф(p) — передаточные функции колебательного и форсирующего звеньев соответственно. Тогда по известным АФХ системы и колебательного звена можно получить АФХ форсирующего звена. Для заданной частоты ω выполним вычитание из угла результирующего радиуса-вектора угла радиуса-вектора колебательного звена, а модуль форсирующего звена получим делением модуля результирующего радиуса-вектора на модуль радиуса-вектора колебательного звена. На рис. 17 показаны АФХ двухфазного повышающего преобразователя (рис. 12), построенная по экспериментальной кривой разгона входного тока id, АФХ колебательного звена с параметрами T = 1,4752×10–4 c и ξ = 0,1223 и АФХ форсирующего звена с постоянной времени t = 4,57×10–4.
На рис. 18 представлены частотные характеристики форсирующего звена, полученные в результате вычислений.
АФХ и частотные характеристики форсирующего звена отличаются от идеальных, что свидетельствует о том, что передаточная функция должна содержать звено более высокого порядка. Но поскольку различия наблюдаются только на частотах, близких к частоте управления, то неидеальностью можно пренебречь и определить постоянную времени форсирующего звена по ФЧХ (рис. 18 а). А именно, для угла 45° находим соответствующее значение частоты, по которой получим:
1/τ = ω = 2,21×103 с–1,
откуда τ = 4,52×10–4 с, что достаточно близко к значению постоянной времени, определенному по производной тока в начале кривой разгона.
Выводы
- Определены динамические характеристики и получены передаточные функции двухфазного повышающего конвертера по напряжению нагрузки и по входному току. Достоверность полученных передаточных функций подтверждена моделированием процессов в системе Simulink MATLAB.
- При определении передаточных функций конвертера необходимо учитывать сопротивление нагрузки. Учет сопротивления потерь в конденсаторе фильтра позволяет получить более точную передаточную функцию, но существенно не влияет на динамические параметры преобразователя.
- На основе полученных передаточных функций можно выполнить синтез параметров регуляторов обратной связи двухфазного повышающего конвертера.
- Четти П. Проектирование ключевых источников электропитания. М.: Энергоатомиздат, 1990.
- Белов Г. А., Павлова А. А., Серебрянников А. В. Динамические модели силовой части двухфазного повышающего импульсного преобразователя / Динамика нелинейных электротехнических и электронных систем: Материалы 8-й Всерос. науч.-техн. конф. – Чебоксары: изд-во Чуваш. ун-та, 2009.
- Ромаш Э. М., Драбович Ю. И., Юрченко Н. Н., Шевченко П. Н. Высокочастотные транзисторные преобразователи. М.: Радио и связь, 1988.
- Севернс Р., Блум Г. Импульсные преобразователи постоянного напряжения для систем вторичного электропитания. М.: Энергоатомиздат, 1988.
- Коршунов А. Анализ способов стабилизации выходного напряжения повышающего импульсного преобразователя постоянного тока // Компоненты и технологии. 2007. № 2.
- Бессонов Л. А. Линейные электрические цепи. М.: Высш. шк., 1968.
- Власов-Власюк Б. М. Экспериментальные методы в автоматике. М.: Машиностроение, 1969.
- Гольдфарб Л. С., Балтрушевич А. В., Нетушил А. В. и др. Теория автоматического управления. М.: Высш. шк., 1976.
- Пономарев В. Ф., Воеводина В. В. Расчет систем автоматического регулирования технологических объектов. Калининград: Калининградский технический институт рыбной промышленности и хозяйства, 1972.