Анализ и синтез широтно-импульсных систем

УДК 68i.5

А. Н. Ловчиков, Е. Е. Носкова АНАЛИЗ И СИНТЕЗ ШИРОТНО-ИМПУЛЬСНЫХ СИСТЕМ

Рассматриваются разработанные авторами методики анализа и синтеза автоматических систем с широтно-импульсной модуляцией.

Ключевые слова: широтно-импульсные системы, широтно-импульсная модуляция, автоматические системы.

Системы автоматического управления с широтноимпульсной модуляцией (ШИМ) относят к классу нелинейных импульсных. Подобные системы известны достаточно давно и в настоящее время широко применяются, однако разработанных методик анализа и синтеза, особенно для инженерных приложений, практически нет.

Предлагаемые методики анализа и синтеза систем с ШИМ позволяют решить большинство из тех задач, которые при этом возникают: исследовать устойчивость системы и в малом, и в большом; изучить динамические характеристики системы с учетом инерционности ШИМ-ключа; дать оценку увеличения или уменьшения запасов устойчивости при работе системы в нелинейных режимах, объяснить причину возникновения в системе многочастотных колебаний; синтезировать корректирующие устройства, обеспечивающие необходимые запасы устойчивости, построить управление, обеспечивающее системе с ШИМ оптимальность по быстродействию.

Разработанные методики анализа и синтеза основываются на представлении системы с ШИМ как существенно нелинейной системы автоматического управления [1], включающей в себя релейно-гистерезисный элемент с зоной нечувствительности, отражающий реальные характеристики любого ключевого элемента, непрерывную нелинейную часть, содержащую звено запаздывания, отражающее инерционные свойства ключа, на которую воздействует постоянная составляющая и периодическое (пилообразное, треугольное или синусоидальное) внешнее воздействие (рис. 1). Передаточная функция ^ос(р) вводится при коррекции динамических характеристик системы.

Такой подход дает возможность применять методы исследования существенно нелинейных систем и на их основе разработать методики анализа и синтеза систем с ШИМ, позволяющий эффективно решать задачи обеспечения устойчивости для линейных и нелинейных режимов, более корректно сформулировать условия перехода от нелинейной к квазинепрерывной модели, решать задачи построения многомодульных широтно-импульсных систем (ШИС).

Рис. 1. Структурная схема широтно-импульсной системы

Система автоматического регулирования (см. рис. 1), относится к классу нелинейных систем, так как содержит

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

В подавляющем большинстве практических случаев передаточная функция линейной части W^p) в ШИС обладает свойством фильтра [1 ]:

\wjjnM)\<<\wjjw)\, n = 2, 3,... .

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

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

Уравнение динамики системы имеет вид [1]

X + e^ wjp) FX = W» Чп -/вын^)

гдеf^t) = в sin( w2 • t).

Решение уравнения находится в виде

X=XB+X,

B м’ v 7

где X = X + X„; X - медленно меняющаяся составляю-

м U A м

щая; X0 - постоянная составляющая; XA = At • sin(rot • t) -переменная составляющая, характеризующая автоколебания; XB = A2 • sin(w2 • t + ф) - переменная составляющая, характеризующая вынужденные колебания в системе; Oj - частота автоколебаний, ю2 = 2л / T - вынужденная частота; T - период переключения.

При использовании принятых обозначений уравнение (1) разбивается на два уравнения соответственно для медленно меняющихся и колебательных составляющих: Q(p) • Xw + R(p)e^ •Woc (p) • F0 =

= (p)•иоп • Q(p) ,

[Q( p) + R( p^W^ (p) • (q + q'p / W2) +

+Q

где принято WRJp>) = R(p)/Q(p).

Нелинейность и коэффициенты гармонической линеаризации определятся из выражений

F[X] = F0(A2,Xм) + [q(A2,XM) + q'(A2) • p/W2] • Xb,

ч(А2, Xм) =

и

А2 -л

1 - + . 1 - (Ь + Хм )2

а2

Ч' (А2) =-

А,'

2 •и1 • Ь

А2 - л '

Р ° (А2 , Хм ) = ^ ^

2 2-л

. Ь + Хм . Ь - Xм ^

агс8іп--------- - агсєіп

А2

а2

2 " *2 0

Амплитуда А2 определяется из равенства

Ш Ю2)2

Ю(/Ю2) + КО2) • в-"”2 • Кс(]”) • (Ч + Іч')2

Решая (3), можно определить функцию смещения

ПАХ) = ФХМ)

для заданной нелинейности -Р[Х], которая в дальнейшем используется при определении автоколебательного режима и решения уравнения для постоянных составляющих. Для этого линеаризуем функцию смещения Ф(Хм) в пределах рабочего режима системы

Ф(Xм ) = Ф0 + Ка • Ха

где Ка =

СФ

С учетом (4), (5) и (6) система уравнений для определения возможности возникновения автоколебательного режима имеет вид

б(0)Х0+Я(0) Гос(0)Ф» = Жо(0)и^(0),

Q(p)xA+к(р)е-тос(р)кхА = о.

После подстановки в (7) р = }^х определяютсяХ0, А1 и w1, т. е. условия возникновения или отсутствия автоколебаний в системе, на основании чего делается заключение об устойчивости системы. Изменяя Жлч(р) или вводя корректирующие устройства в цепь обратной связи, т. е. также изменяя Ж (р) и используя вышеизложенную методику, можно синтезировать рассматриваемую ШИС исходя из необходимых запасов устойчивости.

Наличие вынужденных колебаний в системе приводит к эффекту вибрационной линеаризации нелинейности [3]. В этом случае релейно-гистерезисный элемент приобретает квазинепрерывную нелинейную характеристику, зависящую от параметров непрерывной линейной части и нелинейности. Однако такая нелинейная характеристика с большим основанием может быть линеаризована в рабочей точке, чем релейно-гистерезисная. Если в выражении для коэффициента Р0(4) принять Х .=Х, то согласно (6) коэффициент передачи линеаризованной нелинейности определится выражением

( \

1 1

к = и.

А 2л

^4 - (Ь + Xо)2 ^ - (Ь -Xо):

где А2 - амплитуда вынужденных колебаний; Х0 - постоянная составляющая на входе нелинейности. Для предварительных расчетов значение коэффициента КА приближенно определяется исходя из требований к синтезируемой системе. Если принять Ж (р) = 1, то при заданных

статической точности Х0 и пульсациях Д и на стадии предварительного проектирования А2 определяется прибли-

Таким образом, ШИС можно представить в виде квазилинейной системы с передаточной функцией разомкнутого контура Ж (р), определяемой выражением

Жр(р) = Ка Жлч(р) Жо(р) в-*.

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

Рассмотренная структурная схема является классической при решении задачи построения системы, оптимальной по быстродействию [4]. В соответствии с теорией оптимального управления, произведение егХр Ж (р) является объектом управления (ОУ), ^Х] выдает управляющее воздействие, ограниченное сверху и и снизу 0. Сама функция ЩХ] имеет релейную характеристику, что, в соответствии с принципом максимума Л. С. Понтрягина, необходимо при построении оптимальной по быстродействию системы и наличии ограничений на величину управляющего воздействия.

На основании вышеизложенного задача синтеза ШИС по критерию оптимального быстродействия сводится к определению условий переключения релейного регулирующего элемента ^Х], т. е. нахождению уравнения линии переключения [3; 4]. Применительно к схеме (см. рис. 1) задача синтеза предполагает нахождение такой передаточной функции Ж (р), которая обеспечит оптимальное быстродействие системе.

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

Для простейших случаев уравнение оптимальной фазовой траектории можно получить аналитически. В большинстве практических случаев этого сделать не удается. Методика построения оптимального по быстродействию ШИС в этом случае выглядит следующим образом. По известной передаточной функции объекта управления рассчитываются фазовые траектории и определяется оптимальная. Оптимальная траектория аппроксимируется. Полученное уравнение аппроксимации принимается за уравнение линии переключения теперь уже квазиоптимальной по быстродействию ШИС.

Следует отметить, что сложность решения задачи построения оптимальной по быстродействию системы в значительной степени зависит от вида передаточной функции объекта управления. Изложенная выше методика более справедлива для систем второго порядка без учета запаздывания. Для этого случая важную роль будет играть вид корней характеристического уравнения объекта управления. В случае простых корней, согласно теореме об «-интервалах, оптимальная система отрабатывает возмущение за два интервала постоянства. Если корни комплексные, то количество интервалов постоянства зависит от уровня возмущающего воздействия. Последнее утверждение справедливо и для систем более высокого поряд-

ка, но для таких систем сложно, а часто и невозможно использовать метод фазовых траекторий.

Предложенная выше методика построения оптимальной по быстродействию системы справедлива для больших отклонений и не учитывает наличия внешнего периодического воздействия. Однако и в этом случае, как показано в [3], можно использовать метод фазовых траекторий для анализа поведения ШИС в режиме вынужденных колебаний.

Рассмотрим применение вышеизложенных методик на примере анализа и синтеза импульсного стабилизатора напряжения (ИСН) понижающего типа (рис. 2).

Нелинейная модель ИСН может быть представлена в виде схемы (см. рис. 1). Линейной частью ИСН является однозвенный фильтр с передаточной функцией вход фильтра - выход фильтра с учетом параметров нагрузки:

К

Ж (р) =------г:-------гг,

л.ч^' а2 • р + а1 • р +1

где Кс = Ян/(Ин + г) - коэффициент передачи; Ин - сопротивление нагрузки, г1 - активное сопротивление обмотки дросселя фильтра; Ь - индуктивность дросселя фильтра; С - емкость конденсатора фильтра; а2 = Кс • Ь• С; а1 = (Ь + И • г1 • С)/(Ян + г1). Нелинейность Е[Х] и звено чистого запаздывания имитируют

работу импульсного ключа УТ, управляемого ШИМ. В общем случае ШИМ состоит из компаратора, на один из входов которого подается пилообразное напряжение (и ), а на второй - управляющее (сигнал ошибки) напряжение (иоп - и2). При анализе ИСН считалось, что на вход компаратора вместо и действует синусоидальное напряжение с такой же амплитудой. Величина амплитуды принималась исходя из допустимой статической ошибки менее 1 %. На основании вышеизложенного метода при разработке нелинейной модели ИСН напряжения ипил и и моделируются источником вынужденных колебаний/в (0 и постоянной составляющей соответственно.

Расчеты проводились при следующих исходных данных: частота преобразования 25 кГц; и = 40 В; иоп = 27 В; Ь = 10 мкГн; С = 5 мФ; г1 = 0,010 м; ( = 1 мкс; И = 10 м; В = 0,2 В; Ж (р) = К - передаточная функция звена обратной связи; К = 1; Ь - зона нечувствительности компаратора, Ь = 0,001.

Согласно [2], процессы в ИСН описываются уравнением динамики (1). Решение этого уравнения ищется в виде (2) из системы уравнений (3) для медленно меняющихся и колебательных составляющих. Решая (3), определяем функцию смещения (5) и линеаризуем ее в преде-

лах рабочего режима ИСН (6). Из рассчитанных функций смещения (рис. 3) при различных значениях параметров фгьтра (И-руКТУЕНОСТЖ - рис; 31 а и емкостях - рис. 3, б) следует, что значение КА в рабочей точке равно 14.. .15.

Уравнение (1) теперь примет вид (7) и неизвестные А , ю1 и Х0 определяются из (7) подстановкой р =_/'ю1:

1 - а2 •ю12 + КА • Кс • К • соб(ю1 •т) = 0, а1 •ю1 - КА • Кс • К • •т) = 0,

X0 + Кс • К • Ф0 = К иоп.

При решении (10)/ должно быть меньше 25 кГц. Тогда величина 2р/ •/ = ю14 мала и б1п(ю1 • /) » ю1 •т, а соб(ю1 •т) » 1. Отсюда искомое решение для автоколебаний системы, т.е. условие их существования определится из выражения

г1 • с + ЬУя - К • т • КА = 0.

При подстановке значений параметров в (11) следует, что в системе отсутствуют автоколебания. Однако если совместное запаздывание ключа и компаратора будет равно 4 мкс, то в системе возникнут автоколебания с частотой 2,7 кГц.

Для анализа динамических характеристик ИСН как квазилинейной системы передаточную функцию разомкнутого контура представим согласно (9), где КА = 15 (рис. 4). Построим логарифмическую амплитудную (рис. 4, а) и фазочастотную (рис. 4, б) характеристики ИСН (сплошная линия). Там же пунктирной линией изображены аналогичные характеристики скорректированной системы, когда

^о.с (Р) =

1+Т • Р

1 + Т2 • Р !

при Т1 = 1,25 мс, Т2 = 0,125 мс.

Параметры корректирующего устройства выбирались исходя из заданного запаса по фазе а = 60°, который обеспечивает корректирующее устройство на частоте среза системы.

В соответствии с изложенной выше методикой синтеза оптимальной по быстродействию широтно-импульсной системы, используя передаточную функцию (10), строим фазовые характеристики и определяем оптимальную траекторию (расчеты в пакете МаШсаф (рис. 5). Штриховой линией показаны фазовые характеристики исходной нескорректированной системы, сплошной линией - оптимальная траектория, штрих-пунктирной - ква-зиоптимальная (рис. 5, а). Поскольку аналитически определить выражение, описывающее оптимальную линию переключения, практически невозможно, была проведе-

на ее аппроксимация, которая позволила получить выражение для линии переключения X +1,01 •Ю-9 • X2 • sign(X) = 0. Квазиоптимальная линия переключения (см. рис. 5, а), имеет вид

X +1,3 -10-4 • X = 0. Переходные процессы в ИСН при включении напряжения питания для трех случаев: нескорректированной системы, с квадратичным (сплошная линия) и линейным (штрихпунктирная линия) законами управления электронным ключом ИСН - приведены на рис. 5, б.

а

Рис. 3. Функции смещения: а -

Справедливости ради следует отметить, что результаты расчетов (см. рис. 5) получены без учета нелинейной характеристики диода. В реальной схеме без корректирующего устройства наблюдается режим прерывистых токов и переходной процесс не является колебательным. Однако и в этом случае (рис. 6) (расчет с помощью пакета МютоСАР) введение коррекции с рассчитанными параметрами обеспечивает системе оптимальность по быстродействию, что подтверждает правильность предложенной методики синтеза.

б

б - при различных емкостях фильтра

Рис. 4. Характеристики ИСН: а - логарифмические; б - фазочастотные

5

а

к

я

I

С

я

X

о

о

X

э

а

Л

т

Л ! \ Л

) | ' . 1 \

Г 1 \ 1 1 \ ! '1 7 V /

і V

0 0.001 0.002 Время

а б

Рис. 5. Характеристики нескорректированной и оптимальной системы: а - фазовые характеристики; б - переходные процессы

40

30

20

10

Предложенные методики дают возможность с единых позиций, используя хорошо апробированные методы, с высокой точностью решать задачи анализа и синтеза ШИС: оценивать устойчивость системы и синтезировать ее, исходя из заданных запасов устойчивости, строить оптимальную по быстродействию систему. Данные методики хотя и позволяют уже сейчас решать большое количество задач, связанных с ШИС, однако требуют развития для решения более сложных проблем, например, проведения анализа и синтеза при наличии дополнительных нелинейностей в ШИС, построения оптимальной по быстродействию системы высокого порядка.

Библиографические ссылки

1. Бессекерский В. А., Попов Е. П. Т еория систем автоматического управления. М. : Профессия, 2003.

2. Волович Г. М. Динамика вентильных источников вторичного электропитания постоянного тока. М. : Энер-гоатомиздат, 1991.

3. Попов Е. П. Теория нелинейных систем автоматического регулирования и управления. М. : Наука, 1988.

4. Павлов А. А. Синтез релейных систем, оптимальных по быстродействию (метод фазового пространства). М. : Наука, 1966.

Рис. 6. Переходные процессы в ИСН

A. N. Lovchikov, E. E. Noskova ANALYSIS AND SYNTHESIS OF WIDE-PULSE SYSTEMS

In the article we examine and work out methods of analysis and synthesis of automatic systems with wide-pulse modulation.

Keywords: wide-pulse systems, wide-pulse modulation, automatic systems.

© Ловчиков А. Н., Носкова Е. Е., 2010