Анализ и построение функции Боде преобразователей постоянного напряжения с помощью простейших эквивалентных схем малого сигнала

№ 2’2022
PDF версия
В статье рассматривается анализ и построение функции Боде импульсных преобразователей постоянного напряжения, эквивалентные схемы которых приводятся к простейшим схемам низких порядков. Используются условия, при которых аналитические выражения для этих функций получаются простыми, что удобно применять при проектировании преобразователей постоянного напряжения.

Объединенная эквивалентная трансформаторная схема малого сигнала понижающего-повышающего преобразователя постоянного напряжения представлена в работе [1] и имеет вид, показанный на рис. 1.

Объединенная эквивалентная трансформаторная схема малого сигнала понижающего-повышающего преобразователя

Рис. 1. Объединенная эквивалентная трансформаторная схема малого сигнала понижающего-повышающего преобразователя

Передаточная функция Gvg(s) между выходом и входом преобразователя определяется установкой вариаций рабочего цикла Формула равного нулю и последующим нахождением передаточной функции от Формула до Формула:

Формула

Передаточная функция (1) описывает, как изменения или возмущения приложенного входного напряжения Формула приводят к нарушениям выходного напряжения Формула. Это важно при проектировании регулятора выходного напряжения в преобразователях постоянного напряжения (в конвертерах). Например, в автономном источнике питания входное напряжение vg(t) содержит нежелательные четные гармоники напряжения сети переменного тока. Передаточная функция Gvg(s) используется для определения влияния этих гармоник на выходное напряжение v(t) преобразователя.

Передаточная функция Gvg(s), которая описывает связь между управлением рабочим циклом Формула и выходным напряжением, находится путем установки изменений входного напряжения Формула на ноль и последующего получения модели эквивалентной схемы для определения Формула в функции Формула от [1, 2, 3, 4]:

Формула

Эта передаточная функция описывает, как изменения управляющего входа Формула влияют на выходное напряжение v(t). В системе регулятора входного напряжения Gvg(s) является ключевым компонентом и оказывает существенное влияние на характеристики регулятора.

Выходной импеданс Zвых(s) определяется при условии, что переменные Формула и Формула равны нулю. Zвых(s) описывает, как изменения в токе нагрузки влияют на выходное напряжение. Количественная оценка этой величины также важна для регуляторов напряжения.

Входное сопротивление Zвх(s) преобразователя играет существенную роль, когда на входе питания преобразователя добавлен фильтр электромагнитных помех (ФЭП). Отметим, что величины Zвх(s) и Zвых(s) фильтра в целом влияют на передаточную функцию преобразователя Gvg(s), что необходимо учитывать при проектировании полупроводниковых преобразователей.

Правила построения амплитудных и фазовых асимптот рассматриваются ниже, включая два типа особенностей, которые часто появляются в передаточных функциях преобразователя: резонансы и нули в правой части полуплоскости. Диаграммы Боде для передаточной функции при малосигнальном анализе повышающего преобразователя приведены в работах [3, 4, 5, 6, 7, 8].

Трудность, с которой сталкиваются при анализе схем, — это сложность модели схемы: практическая схема может содержать сотни элементов, и, следовательно, их анализ может привести к сложным выводам, неразрешимым уравнениям и множеству алгебраических ошибок. Некоторые инструменты, предназначенные для проектирования сложной преобразовательной системы, описаны в работах [1, 2, 3, 4, 5]. Представление передаточных функций в нормализованной форме напрямую раскрывает важные особенности реакции схемы. Аналитические выражения для этих функций и асимптот приводят к простым уравнениям, которые можно использовать при проектировании преобразователей.

Хорошо разделенные корни полиномов передаточной функции могут быть аппроксимированы простым способом. Ниже описывается графический метод построения графиков Боде передаточной функции по существу, путем проверки. Этот метод может:

  • уменьшить объем алгебры и связанной с ней алгебраических ошибок;
  • привести к более глубокому пониманию поведения схемы;
  • привести к пониманию, необходимому для получения подходящих приближений, которые сделают уравнения легко описываемыми.

Для лучшего усвоения приведенных вопросов приведем краткий обзор графиков Боде для простых схем.

Величина передаточной функции G в безразмерной форме в децибелах имеет следующий вид:

Формула

Когда описываемый параметр имеет размерность, эту величину необходимо нормализовать с помощью деления на базовую величину. Например, для представления сопротивления в децибелах следует выполнить деление на величину базового сопротивления:

Формула

где Rбаз — величина базового сопротивления. Значение величины базового сопротивления выбирают произвольно, но нужно уточнить, какое базовое значение используется. Например, если ||Z|| = 5 Ом и выбрано Rбаз = 10 Ом, то ||Z||дБ = 20lg(5/10) = –6 дБ относительно 10 Ом.

В практических расчетах для удобства базовое сопротивление выбирается Rбаз = 1 Ом и считается, что импеданс при этом выражен в децибел-омах (дБОм). Таким образом, при импедансе 5 Ом его эквивалент равен ||Z||дБ = 20lg(5/1) = 14 дБ.

При оценке гармоник входного тока преобразователя часто выбирают базовый ток 1 мкА. Например, если ток равен 1000 мкА, то его эквивалент в дБ равен 60 дБ.

Рассмотрим простую однополюсную схему, простой RC-фильтр нижних частот, показанный на рис. 2.

Простой RC-фильтр нижних частот

Рис. 2. Простой RC-фильтр нижних частот

Передаточная функция определяется отношением выходного и входного напряжений:

Формула

Передаточная функция фильтра совпадает со следующей стандартной нормализованной формой для одного полюса:

Формула

где параметр w0 = 2pf0, находится приравниванием коэффициентов при s в формулах (5) и (6):

w0 = 1/RC.                          (7)

Поскольку R и C вещественные положительные величины, то w0 тоже вещественная положительная величина. Знаменатель уравнения (6) содержит корень при s = –w0 и, следовательно, G(s) содержит действительный полюс в левой половине комплексной плоскости. Чтобы найти величину и фазу передаточной функции, положим s = jw (j2 = –1). При s = jw из (4) получим:

Формула

Формула

и, соответственно, в децибелах получим вид:

Формула

Самый простой способ изобразить график Боде для G(jw) — исследовать асимптотическое поведение этой функции при высоких и низких частотах.

Для малой частоты w<< w0, или f<<f0 и

Формула

с учетом (10) и (11) получим:

Формула

При высоких частотах

Формула

Формула

Из (14) получим:

Формула

Отклонение точной кривой от асимптот можно найти, просто вычислив уравнение (9) на угловой частоте w = w0:

Формула

В децибелах величина равна:

Формула

На рис. 3 приведены графики функции Боде для амплитуды и фазы в соответствии со схемой на рис. 2, где красным цветом показан график асимптоты, который в частоте f0 имеет максимальное отклонение 3 дБ.

Резюме графика амплитуды и фазы Боде

Рис. 3. Резюме графика амплитуды и фазы Боде

Фаза передаточной функции G(jw) определяется из (6) и равна:

Формула

Фаза ∠G(jw) стремится к нулю на низкой частоте и к –90° на высокой частоте. При частоте f = f0G(jw) = –45°. Поскольку высокочастотная и низкочастотная фазовые асимптоты не пересекаются, нужна третья асимптота для аппроксимации фазы на угловой частоте f0. Один из способов сделать это проиллюстрирован на рис. 3, где наклон асимптоты выбран таким, чтобы он совпадал с наклоном фактической кривой при частоте f = f0. Можно показать, что при таком выборе частоты пересечения асимптот fa и fb имеют вид [2, 3]:

Формула

Более простой вариант, с приемлемой точностью соответствующий реальной кривой, это:

fa = f0/10;

fb = 10 f0.                           (20)

Как видно на рис. 3, где эта асимптота сравнивается с реальной кривой, полюс вызывает изменение фазы в диапазоне частот примерно на две декады. Наклон асимптоты в этом диапазоне частот составляет –45° на декаду. В частотах излома fa и fb фактическая фаза отклоняется от асимптот примерно на 5,7°.

Обычно рекомендуется выражать однополюсные передаточные функции в нормализованной форме.

Уравнение (6) легко интерпретируется в w0 = 2pf0 благодаря его нормализованной форме. На низких частотах, где член (s/w0) мал по величине, передаточная функция приблизительно равна 1. На высоких частотах (s/w0>>1) передаточная функция приблизительно равна (s/w0)–1 или (f/f0)–1 (w0 = 2pf0). Таким образом, передаточная функция записывается непосредственно в терминах ее основных характеристик, то есть ее асимптот и угловых частот.

 

Схема с одиночным нулевым откликом

Схема с одним нулевым откликом содержит корень в числителе передаточной функции и может быть записана в следующей нормализованной форме:

Формула

Формула

Фаза задается:

Формула

При низких частотах фаза стремится к нулю, а при высоких частотах к +90°. Вне интервала частот f0/10<f<10f0 фазовая асимптота имеет наклон +45°/декада.

Нули правой полуплоскости

При малосигнальном анализе нули в передаточных функциях часто встречаются в правой полуплоскости:

Формула

Корень уравнения (24) положительный и, следовательно, лежит в правой половине комплексной s-плоскости. Передаточная функция имеет абсолютную величину аналогично (21):

Формула

Следовательно, только по абсолютному значению невозможно различить, находится ноль в правой полуплоскости или в левой полуплоскости.

Фаза задается ∠G(jw) = –arctg(w/w0), и это совпадает с фазой схемы с одиночным нулевым откликом по формуле (18). То есть при нуле в правой полуплоскости отклик по амплитуде будет таким же, как и при нуле в левой полуплоскости, но фазовая характеристика аналогична полюсу. Амплитудно-фазовые характеристики приведены на рис. 4.

Амплитудно-фазовые характеристики при нуле в правой полуплоскости

Рис. 4. Амплитудно-фазовые характеристики при нуле в правой полуплоскости

Две другие формы (перевернутый полюс и перевернутый ноль) возникают из-за инверсии оси частот.

Перевернутый полюс имеет передаточную функцию:

Формула

Преобразуя (25), получим:

Формула

Отметим, что перевернутый полюс имеет усиление на высоких частотах, равное 1, а асимптота на низких частотах имеет наклон +20 дБ/декада.

Перевернутый ноль имеет передаточную функцию:

Формула

Преобразуя (27), получим:

Формула

Примером использования передаточной функции этого типа является пропорционально-интегральный контроллер при построении связи с конструкцией контура обратной связи [3, 4].

Отметим, что перевернутый ноль имеет усиление на высоких частотах, равное 1, а асимптота на низких частотах имеет наклон –20 дБ/декада.

 

Квадратичный полюсный отклик (резонанс)

Рассмотрим теперь передаточную функцию G(s) двухполюсного фильтра частот (рис. 5). Отметим, что большинство понижающих-повышающих преобразователей содержат фильтр такого типа [2, 3, 8].

Двухполюсный фильтр частот

Рис. 5. Двухполюсный фильтр частот

Передаточная функция цепи (рис. 5) имеет следующий вид:

Формула

Передаточная функция (29) в знаменателе содержит многочлен второго порядка. Обозначая a1 = L/R, a2 = LC в (29), получим:

Формула

Чтобы построить график Боде этой передаточной функции, попробуем разложить знаменатель на два его корня:

Формула

Использование квадратичной формулы приводит к следующему выражению для корней:

Формула

Формула

Передаточная функция по уравнению (29) может быть записана в следующей стандартной нормализованной форме:

Формула

где w0 = 1/√LC — собственная частота последовательного LC-контура, Q = R/r — добротность контура, r = L/C — волновое сопротивление контура.

Абсолютное значение G(s) определяется выражением:

Формула

На низких частотах w/w0<<1 и, следовательно, ||G(s)|| 1 при w<< w0.

При высоких частотах w/w0>>1 член (w/w0)4 доминирует в выражении внутри радикала уравнения (35). Следовательно, высокочастотная асимптота имеет вид:

Формула

Таким образом, высокочастотная асимптота имеет наклон –40 дБ/декада, асимптоты пересекаются в точке f = f0 и не зависят от параметра Q. Параметр Q влияет на отклонение фактической кривой от асимптот в окрестности частот f0 (рис. 6).

Графики Боде передаточной функции для схемы второго порядка

Рис. 6. Графики Боде передаточной функции для схемы второго порядка

Основные особенности графика Боде передаточной функции второго порядка приведены на рис. 6. Как видно из этих графиков, при увеличении добротности контура Q отклонение фактических кривых от асимптот в окрестности частот f0 увеличивается.

Фаза функции G(s) определяется:

Формула

Фаза стремится к 0° на низкой частоте и к –180° на высокой. При частоте f = f0 фаза равна –90°. Как показано на рис. 6, увеличение значения добротности Q вызывает более резкое изменение фазы между асимптотами 0° и 180°. Поэтому необходима асимптота на средней частоте, чтобы аппроксимировать фазовый переход в окрестности угловой частоты f0. Как и в случае реального одиночного полюса, можно выбирать наклон этой асимптоты идентичным наклону реальной кривой при частоте f = f0. Можно показать, что такой выбор приводит к следующим частотам изломов асимптот [3]:

Формула

Лучшим выбором, который согласуется с приближением (20), используемым для реального однополюсного соединения, является:

Формула

При таком выборе асимптота средних частот имеет наклон Q –180°/декада.

 

Приближение при низких значениях добротности Q

Как упоминалось выше, когда корни полиномиального уравнения второго порядка знаменателя (30) действительны, можно разложить знаменатель на множители и построить диаграмму Боде, применив асимптоты для действительных полюсов. В дальнейшем используется следующая нормализованная форма:

Формула

Это особенно желательно, когда угловые частоты w1 и w2 хорошо разделены по значению [3, 4, 7]. Сложность данной процедуры заключается в выражении угловых частот w1 и w2 в зависимости от элементов схемы R, L, C, что часто приводит к сложным и непонятным выражениям, особенно когда схема содержит много элементов. Даже в случае простой схемы (рис. 5), передаточная функция которой задается уравнением (30), эта обычная квадратная формула приводит к следующей сравнительно сложной формуле для угловых частот:

Формула

Данное уравнение практически не дает представления о том, как угловые частоты зависят от значений элементов. Например, можно показать, что, когда угловые частоты w1 и w2 хорошо разделены по значениям, они могут быть выражены с высокой точностью гораздо более простыми соотношениями [3, 4]:

Формула

В этом случае w1, по существу, не зависит от значения емкости C, а w2 — от значения индуктивности L, однако уравнение (38) очевидно показывает, что обе угловые частоты зависят от всех значений элементов. Простые выражения (39) намного предпочтительнее, чем уравнения (38), и могут быть легко получены с использованием приближения низкой добротности.

Предположим, что передаточная функция была выражена в стандартной нормализованной форме (30). Для Q ≤ 0,5 воспользуемся формулой корней квадратного уравнения, чтобы записать действительные корни много­члена знаменателя уравнения (34):

Формула

Формула

Угловую частоту w2 можно выразить следующим образом:

Формула

где:

Формула

Отметим, когда Q<< 0,5, тогда 4Q2<< 1 и можно с высокой точностью определить, что когда Q быстрее уменьшается ниже Q< 1, тогда F(Q) стремится к 1 (рис. 7).

Зависимость функции F(Q) при малых значениях добротности

Рис. 7. Зависимость функции F(Q) при малых значениях добротности

Используя (43), (45) и (46), получим следующее выражение для w1:

Формула

Так как F(Q) стремится к 1, когда Q<< 0,5, то, следовательно, из (47) получим:

w1 = Qw0, для Q<< 0,5.           (48)

Асимптотические величины для случая малых добротностей представлены на рис. 8; при Q< 0,5 два полюса в w0 разделяются на реальные полюса: один действительный полюс находится на угловой частоте w1< w0, а другой — на угловой частоте w2 > w0.

Асимптотические величины для случая малых добротностей

Рис. 8. Асимптотические величины для случая малых добротностей

Угловые частоты легко аппроксимируются с помощью (45) и (48). Для схемы на рис. 5 параметры Q и w0 задаются формулой (34). Для случая Q< 0,5 можно получить следующие простые аналитические выражения для угловых частот w1 и w2, используя (45) и (48):

Формула

Таким образом, приближение при низких значениях добротности позволяет получить простые аналитические выражения для определения угловых частот.

Рассмотренный способ получения передаточных функций и угловых частот можно обобщить и применять для определения параметров более сложных схем.

Приближение низкого значения добротности можно использовать, чтобы найти приближенное аналитическое выражение для корней n-порядка:

P(s) = 1 + a1s + a2s2 + … + ansn.       (50)

Желательно разложить полином P(s) на множители в виде:

P(s) = (1 + t1s)(1 + t2s)…(1 + tns).      (51)

Постоянные коэффициенты времени t1, t2tn можно связать с исходными коэффициентами a1, a2an и произвести умножение в уравнении (51), в результате получим:

a1 = t1 + t2 +… + tn,

a2 = t1(t2 +…+ tn) + t2(t3 +…+ tn) +…,

a3 = t1t2(t3 +…tn) + t2t3(t4 + … + tn) + ….,

an = t1t2t3 ++ tn.                   (52)

Общее решение этой системы уравнений сводится к точному разложению полинома произвольной степени на множители, что является безнадежной задачей. Тем не менее уравнение (52) дает способ аппроксимации корней. Предположим, что все константы t1, t2tn действительны и хорошо разделены по значению. Без ограничения общности далее можно предположить, что константы расположены в порядке убывания величины:

|t1|>>| t2|>>| t3|>>…>>| tn|.          (53)

С учетом (53) из (52) получим:

a1t1,

a2t1t2,

a3t1t2t3,

an t1t2t3 tn.                            (54)

Из (54) определяются константы времени t1, t2, t3tn в следующем виде:

Формула

Соответственно, если:

Формула

то полином P(s) можно разложить на множители в следующем виде:

Формула

При нарушении первого неравенства выражения (56) первые два корня следует оставить в квадратичной форме.

В качестве примера рассмотрим схему фильтра электромагнитных помех (ФЭМП). На рис. 9 показана упрощенная схема входного каскада ФЭМП.

Упрощенная схема входного каскада ФЭМП

Рис. 9. Упрощенная схема входного каскада ФЭМП

Как известно, такие фильтры обычно размещаются на входе мощных преобразователей (конвертеров) для ослабления коммутационных гармоник, присутствующих во входном токе конвертера. Используя известные методы расчета электрических схем, можно легко получить для фильтра схемы на рис. 9 передаточную функцию в следующем виде:

Формула

Как видно из (58), знаменатель передаточной функции содержит полином третьей степени с постоянными коэффициентами:

Формула

С учетом (56), (57) знаменатель передаточной функции (58) можно представить в следующем виде:

Формула

Согласно выражению (54), приближение справедливо при условии, что:

Формула

Эти неравенства не могут быть удовлетворены без удовлетворения условия L1 >>L2. При L1 >>L2 условия (61) можно упростить:

Формула

Используя (62) из (60), получим:

Формула

Таким образом, в случае если передаточная функция содержит три хорошо разделенных полюса, уравнения (60) и (63) представляют собой приблизительные аналитические факторизации знаменателя уравнения (58). Необходимость выполнения условий (61) и (62) для получения простых аналитических выражений в зависимости от величин элементов L1, L2, R и C позволяет с помощью (63) выбирать значения параметров элементов, чтобы получить заданные частоты полюсов.

Когда второе неравенство уравнения (61) нарушается, то второй и третий корни следует оставить в квадратичной форме:

Формула

Это приближение оправдано при условии, что:

Формула

Выражение (65) можно упростить при условии:

Формула

Заметим, что больше не требуется выполнять условие RC >> L2/R. Уравнение (66) позволяет дополнительно упростить (64) и привести к виду:

Формула

Таким образом, для этого случая передаточная функция содержит низкочастотный полюс, хорошо отделенный от высокочастотной квадратичной пары полюсов.

На рис. 10 приведены графические зависимости функции Боде для схемы, показанной на рис. 9, при условии (66) согласно выражениям (67) и (60) соответственно. В графиках (рис. 10) красным цветом показаны приближенные кривые согласно (58) при знаменателе (67). Как видно из этих графиков, при соблюдении условий (66) отклонения между реальной и приближенной функциями незначительны.

Графические зависимости функции Боде при условии (66)

Рис. 10. Графические зависимости функции Боде при условии (66)

В случае невыполнения первого неравенства в условии (61) первый и второй корни следует оставить в квадратичной форме:

Формула

Это приближение оправдано при условии, что:

Формула

Или в более упрощенном виде:

Формула

Для этого случая передаточная функция содержит низкочастотную квадратичную пару полюсов, которые хорошо отделены от высоко­частотного реального полюса.

Если приведенные выше приближения не наблюдаются, то все три корня приблизительно равносильны по величине и в этом случае необходимо найти другие способы для упрощенного получения и расчета передаточной функции преобразователей [6–8].

 

Выводы

  1. Представление передаточных функций в нормализованной форме напрямую раскрывает важные особенности реакции схемы. Аналитические выражения для этих функций и асимптот приводят к простым уравнениям, которые можно использовать при проектировании преобразователей.
  2. Рассмотренные выше способы для получения передаточных функций и угловых частот можно обобщить и применять для более сложных схем, для этого необходимо принимать соответствующие приемлемые допущения.
Литература
  1. Петросян Н. Н., Бегоян К. В., Кароян Г. С. Анализ повышающих-понижающих конвертеров с помощью эквивалентной трансформаторной схемы малого сигнала // Силовая электроника. 2019. №. 2.
  2. Hsu S., Brown A., Rensink L., Middlebrook R. D. Modeling and Analysis of Switching Dc-to-Dc Converters in Constant-Frequency Current Programmed Mode. IEEE Power Electronics Specialists Conference, 1979 Record.
  3. Ericson R. W., Maksimovic D. Fundamentals of Power Electronics. Springer Science & Business Media, 2007.
  4. Петросян Н. Н., Мелконян В., Бегоян К. В. Анализ работы цепей управления и построение простой модели малого сигнала повышающе-понижающего конвертера с помощью алгебраического подхода // Силовая электроника. 2019. №6.
  5. Middlebrook R. D. Topics in Multiple-Loop Regulators and Current-Mode Programming. IEEE Power Electronics Specialists Conference, 1985 Record.
  6. Middlebrook R. D. Modeling Curreent Programmed Buck and Boost Regulatorts // IEEE Power Electronics Specialists Conference. 1989. Vol. 4. No. 1.
  7. Vergnese G., Bruzos C., Mahabir K. Averaged and Sampled-Data Models for Current Mode Control: A Reexamination. Power Electronics Specialists Conference, 1989 Record.
  8. Brown A., Middlebrook R. D. Sampled-Data Modeling of Switching Regulators. IEEE Power Electronics Specialists Conference, 1981 Record.

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *