
- •Теория автоматического управления
- •Рецензент кандидат технических наук, доцент, заведующий кафедрой автоматики и телемеханики фгбоу впо «ВятГу» в. И. Семёновых
- •Введение
- •Исходные данные
- •Статическая характеристика тп
- •2. Задание для выполнения курсовой работы (содержание
- •3. Методические указания по выполнению курсовой работы
- •4. Пример выполнения курсового проекта Составление математического описания системы
- •4.1.1. Схема и исходные данные
- •4.1.2. Уравнения во временной области и их операторные преобразования. Нахождение передаточных функций для всех элементов системы
- •4.1.3. Полная и линеаризованная структурные схемы
- •4.1.4. Определение численных значений коэффициентов связи и постоянных времени неизменяемой части системы
- •4.2. Анализ установившегося режима системы
- •4.2.1. Составление структурной схемы для установившегося режима
- •4.2.2. Определение необходимого коэффициента передачи
- •4.2.3. Определение значения переменных для режима с заданной частотой nЗад
- •4.2.4. Определение статических отклонений n частоты вращения вала привода от заданного значения nЗад
- •4.3. Исследование динамики системы
- •4.3.1. Построение аппроксимированной лачх системы в разомкнутом состоянии и проверка устойчивости
- •4.3.2. Построение желаемой лачх системы в разомкнутом состоянии
- •4.3.3. Нахождение лачх последовательного корректирующего устройства, определение передаточной функции корректирующего устройства
- •4.3.4. Расчёт кривой переходного процесса на эвм
- •Оценка качества переходного процесса
- •Теория автоматического управления
2. Задание для выполнения курсовой работы (содержание
пояснительной записки)
Составление математического описания системы
1.1. Вычертить схему (рис. 1) и привести исходные данные выполняемого варианта работы.
1.2. Записать уравнения во временной области, привести их операторные преобразования и найти передаточные функции для всех элементов системы.
1.3. Составить полную и линеаризованную структурные схемы.
1.4. Определить численные значения коэффициентов связи и постоянных времени неизменяемой части системы.
Анализ установившегося режима
2.1. Составить структурную схему для установившегося режима.
2.2. Определить необходимый коэффициент передачи замкнутого контура системы kР обеспечивающий в замкнутой системе статическую ошибку регулирования nЗ не более 5 % от регулируемой частоты вращения, если задан диапазон регулирования D=10:1 при максимальной частоте вращения, равной номинальной nН.
Определить требуемый коэффициент усиления kI операционного усилителя У1.
2.3. Определить значение переменных (Ud,U3,U2,U1,UТГ,UЗАД) для режима работы с заданной частотой вращения nЗАД.
2.4. Определить статические отклонения n частоты вращения вала привода от заданного значения nЗАД в установившемся режиме в системе без обратной связи (nР) и в замкнутой системе (nЗ) при ступенчатом воздействии основных возмущений МС=МН и UC=0,1· Ud ЗАД . Дать сравнительную оценку точности поддержания частоты вращения в статике для разомкнутой и замкнутой системы.
Исследование динамики системы
3.1. Построить аппроксимированную логарифмическую амплитудно-частотную характеристику (ЛАЧХ) системы в разомкнутом состоянии. Проверить устойчивость замкнутой системы.
3.2. Построить желаемую ЛАЧХ системы в разомкнутом состоянии, исходя из следующих требований: перерегулирование скорости по управляющему воздействию не более 20…30 %; время переходного процесса – минимальное; и ослабление всех возмущений – максимальное во всем диапазоне регулирования.
3.3. Найти ЛАЧХ последовательного корректирующего устройства, определить передаточную функцию корректирующего устройства, выбрать схему и определить параметры Z1 и Z2 корректирующего устройства.
3.4. Рассчитать кривую переходного процесса, отклонение nЗ частоты вращения вала привода от заданного значения nЗАД при приложении заданного возмущения. Расчет производить моделированием САУ на ЭВМ. По полученной кривой оценить качество переходного процесса.
3. Методические указания по выполнению курсовой работы
К первой части
Математическое описание системы приводится на основе составления системы дифференциальных уравнений для элементов системы при общепринятых допущениях.
Электродвигатель постоянного тока независимого возбуждения с рабочей машиной при допущениях постоянного магнитного потока двигателя с компенсированной реакцией якоря двигателя, с абсолютно жесткими механическими соединениями и постоянном моменте инерции J привода описывается системой уравнений. Система состоит из дифференциальных уравнений баланса напряжений якорной цепи и движения привода, а так же соотношений между ЭДС Е и частотой вращения n, электромагнитным моментом М и током I якорной цепи
(1.1)
где CЕ – коэффициент, учитывающий магнитный поток и конструктивное исполнение электродвигателя.
Тиристорный преобразователь ТП с СИФУ описывается неоднородным дифференциальным уравнением первого порядка
,
(1.2)
где ТП – постоянная времени (ТП=0,05 с для всех вариантов);
kП – коэффициент определяемый статической характеристикой ТП.
Операционные усилители У1, У2 описываются передаточной функцией
,
(1.3)
где ZОС(р) и ZВХ(р) – операторные сопротивления цепей обратной связи и входного данного операционного усилителя.
При этом операторное сопротивление активной цепи равно R, емкостной – 1/СР , индуктивной – LР. При сложных соединениях этих сопротивлений применяют правила, аналогичные правилам расчета соединений действительных сопротивлений. Если на вход операционного усилителя (У2 в задании) прикладывается несколько различных воздействий U1, U2 и т. д. с входными сопротивлениями Z1, Z2 и т. д., то операционный усилитель описывается соотношениями
(1.4)
где U1, U2 … берутся со знаком, соответствующим полярности напряжения на входе операционного усилителя, т. е. при ZОС=Z1=Z2=… операционный усилитель может служить сумматором в алгебраическом смысле слова.
Тахогенератор представляется в виде линейного безинерционного усилительного звена с передаточным коэффициентом kТГ.
Фильтр на выходе тахогенератора можно рассматривать как отдельное звено с входным напряжением UТГ и выходным U1. Такое допущение основывается на том, что внутреннее сопротивление тахогенератора можно считать пренебрежимо малым, а нагрузочное сопротивление R2 на порядок больше внутреннего сопротивления фильтра. Полярность UТГ подбирается такой, чтобы в установившемся режиме сигнал обратной связи на входе У2 был обратным по знаку сигналу UЗАД. Постоянная времени фильтра ТФ определяется произведением СфRФ.
Рекомендуется после математического описания элемента системы во временной области приводить его изображение по Лапласу при нулевых начальных условиях, передаточную функцию и структурную схему – полную с графическим изображением статической характеристики, если таковая имеется, и линеаризованную – с указанием конкретной передаточной функции элемента.
При определении
передаточной функции электродвигателя
за выходную переменную следует принимать
частоту вращения вала n,
за входное задающее воздействие –
напряжение Ud
, в качестве возмущающих воздействий –
момент статического сопротивления МС
и отклонения напряжения сети
UС.
Следует учесть, что динамические свойства
электродвигателя характеризуются двумя
постоянными времени: электромагнитной
ТЯ=LЯЦ/RЯЦ
и электромеханической ТМ=J·RЯЦ/СЕ2.
Коэффициент передачи двигателя
относительно Ud
определяется соотношением kd=1/CE.
Сравнивая передаточные функции
электродвигателя по задающему воздействию
и колебательного звена, можно определить
коэффициент затухания
.
При 0<<1
двигатель представляется колебательным
звеном, при >1
– апериодическим звеном второго порядка.
Собственная частота d
колебаний электродвигателя равна
.
В целом рассматриваемая система электропривода представляет собой одноконтурную замкнутую САУ с последовательным корректирующим устройством.
Прямой канал полной структурной схемы САУ включает в себя последовательное соединение звеньев: сумматора У2 с входными воздействиями: задающим UЗАД, обратной связью U1 и выходным воздействием U2; коррекции У1 со входом U2 и выходом U3 ; тиристорного преобразователя, состоящего из инерционного звена 1-го порядка с коэффициентом усиления, равным 1, и безинерционного линейного звена с таблично заданной нелинейной статической характеристикой преобразователя (входным воздействием преобразователя служит U3 , а выходным – Ud); электродвигателя с указанными выше воздействиями Ud, n, UС, Мс.
Канал обратной связи состоит из усилительного звена тахогенератора ТГ с коэффициентом передачи kтг, входным воздействием n, выходным Uтг и из фильтра в виде инерционного звена 1-го порядка.
Линеаризованная структурная схема САУ получается, если полные переменные X(U3АД , U2 , U3 , Ud , UС, I, М, Мс, n, Uтг, U1 ) представить в виде суммы ХА+Х, где ХА – значение переменных в рабочей точке статических характеристик звеньев. После сокращения статических составляющих левой и правой части уравнений звеньев линеаризованные структурные схемы звеньев отразят зависимость между приращениями входных и выходных переменных ( Uзад, U2, …, U1 ). Форма записи передаточных функции линейных звеньев при этом не изменяется, а статические характеристики нелинейных безинерционных звеньев будут представлены в виде коэффициентов динамической линеаризации в рабочей точке.
При расчёте численных значений коэффициентов передачи и постоянных времени звеньев САУ следует приводить формулу в буквенных обозначениях параметров, затем то же выражение с подставленными в него значениями параметров САУ и после этого конечный результат с указанием размерности в системе СИ. Конструктивный коэффициент СЕ электродвигателя может быть рассчитан по уравнению баланса напряжений якоря двигателя в установившемся номинальном режиме, т. е.
,
(1.5)
где переменные Uн, Iн, nн определены данными электродвигателя (табл. 2), а сопротивление RЯД собственно цепи якоря двигателя составляет 0,5 RЯЦ сопротивления всей якорной цепи (табл. 2).
Из графической
зависимости Ud
(Uз)
(табл. 1), построенной в одинаковом
масштабе переменных обеих осей,
определяется коэффициент kп
динамической линеаризации статической
характеристики тиристорного
преобразователя.
,
или графически kп
равен тангенсу угла наклона касательной,
проведённой к статической характеристике
рабочей точки А. Рабочая точка А
определяется значением выпрямленного
напряжения Ud|A
в режиме идеального холостого хода
электродвигателя. Необходимо определить
значение kп
в трёх рабочих точках:
kп мин – при максимальном значении Ud в заданном диапазоне регулирования D, т. е. Ud мин=Uн/D, kп мин используется при определении коэффициента k1 корректирующего звена (п. 2.2);
k п макс – в точке с максимальным наклоном статической характеристики в пределах диапазона регулирования, kп макс используется для определения динамического сомножителя передаточной функции корректирующего звена (пп. 3.1…3.3);
kп зад – при заданном значении выпрямленного напряжения, т. е. при
Ud зад=СЕ·nзад, kп зад используется для построения ЛАЧХ САУ в заданном рабочем режиме.
Ко второй части задания
Структурная схема для установившегося режима составляется на основе уравнений элементов САУ в статике или на основе линеаризованной структурной схемы САУ формальным путём приравнивания оператора p к 0. Следует обратить внимание, что последним способом установившийся режим работы двигателя может быть описан лишь на основе полной передаточной функции двигателя по отношению к Ud и Мс с выходом по частоте вращения, т. е. с учётом внутренней обратной связи двигателя.
По полученной структурной схеме в установившемся режиме можно определить статическое отклонение частоты вращения n вала при приложении Мс=Мн в разомкнутой системе без обратной связи (без ТГ и фильтра) – nрМс и в замкнутой САУ nзМс. Нетрудно убедиться, что соблюдается равенство
,
(2.1)
kр=k1 · kп·kd·kтг , (2.2)
где kр – коэффициент передачи замкнутого контура САУ в разомкнутом состоянии.
Отклонение nрМс рассчитывается непосредственно из структурной схемы в установившемся режиме, при этом, согласно (1.1), Мс=Мн=СЕ·Iн.
Учитывая, что относительное падение частоты вращения в статике при приложении Мс имеет максимальную величину при минимальной частоте вращения в пределах заданного диапазона регулирования D и, исходя из требований п. 2 задания, имеем соотношение
nзМс0,05 nн / D. (2.3)
Тогда необходимый коэффициент передачи kр мин может быть найден из уравнения (2.1), а требуемый из условий статики коэффициент kI операционного усилителя У1 из уравнения (2.2) (учитывая, что kр мин в качестве сомножителя имеет коэффициент передачи тиристорного преобразователя, равный kп мин, найденный в первой части курсовой работы).
Рекомендуется после определения kI найти, согласно (п. 2.2), коэффициенты kр макс и kр зад, при этом в соотношении (2.2) меняется лишь значение коэффициента передачи тиристорного преобразователя (k п макс и
kп зад). Значения kр зад используется при выполнении пп. 2.4; kр макс – при выполнении пп.3.1…3.3 задания.
К третьей части задания
Как известно из теории автоматического управления, устойчивость замкнутой САУ может быть определена с помощью критерия Найквиста, который в случае рассматриваемой устойчивой САУ в разомкнутом состоянии и анализе её частотных характеристик в логарифмическом масштабе сводится к условию, чтобы
|р(с)|<180 , (3.1)
где частота среза с системы определяется значением частоты точки пересечения ЛАЧХ САУ в разомкнутом состоянии и оси абцисс (Lр(с)=0).
Показатели качества переходного процесса в замкнутой САУ зависят от запаса устойчивости по фазе (=180–|р(с)|) и по амплитуде L(L равно абсолютному значению Lр при частоте, где |р|=180). В частности это иллюстрируется на рис. 2, где представлены зависимости перерегулирования в замкнутой САУ по задающему воздействию и L.
Рис. 2. Кривые запаса устойчивости по модулю L и по фазе
от перерегулирования
В работе рассматривается минимально – фазовая система, т. е. система, состоящая в разомкнутом состоянии из звеньев, имеющих однозначную зависимость между АЧХ и ФЧХ. Для такой системы анализ устойчивости можно провести по аппроксимированной ЛАЧХ системы в разомкнутом состоянии без построения ФЧХ. Значение i при любой выбранной частоте i для минимально-фазовой системы может быть приближённо определено по усреднённому наклону ЛАЧХ в этой частоте и равно
iсрi·90, (3.2)
.
(3.3)
Усреднённый наклон срi определяется путём нахождения координаты LВi, отстоящей на одну декаду в сторону высоких частот, т. е. LВi=L10i, и координаты LНi, отстоящей на одну декаду в сторону низких частот, т. е. LНi=L0,1i.
Некорректированная САУ в разомкнутом состоянии состоит из последовательно соединённых звеньев:
– усилительного – усилителя У1 с коэффициентом передачи kI, найденным из условия статики во второй части задания;
– инерционного – тиристорного преобразователя с коэффициентом передачи kп и частотой сопряжения п=1/Тп;
– колебательного
– электродвигателя с коэффициентом
передачи kd
и собственной частотой колебаний
;
– усилительного – тахогенератора с коэффициентом передачи kтг;
– инерционного – фильтра с единичным коэффициентом передачи и частотой сопряжения ф=1/Тф.
При проверки устойчивости некорректированной САУ и построении желаемой ЛАЧХ системы в разомкнутом состоянии (пп. 3.1 и 3.2) коэффициенты передачи всех звеньев (в том числе и желаемой kI) могут быть объединены в соответствии с п. 2.2. в один коэффициент kр. Следует учитывать, что наихудшие условия, с точки зрения устойчивости, при заданных параметрах динамических звеньев будут находиться при максимальном значении kр в требуемом диапазоне регулирования, т. е. при k р макс=kр (см. указания ко второй части).
Процесс построения аппроксимированной ЛАЧХ некорректированной САУ в разомкнутом состоянии можно ускорить, если воспользоваться следующей методикой:
определить значение ординаты Lp(0), равное lgkр макс, и абсцисс частот сопряжения и колебаний lgп, lgd, lgф;
на низких частотах графика Lpнк() отложить ординату, равную Lp(0), и провести через эту ординату прямую с нулевым наклоном (параллельно оси абсцисс) до ближайшей меньшей собственной частоты одного из звеньев;
по мере увеличения частоты изменять наклон ЛАЧХ на –1 лог/дек в абсциссах, соответствующих частотам, сопряжения инерционных звеньев, и на –2 лог/дек в абсциссах, соответствующих собственной частоте колебательного звена.
На рис. 3 приведена наиболее типичная для заданных вариантов ЛАЧХ Lpнк() некорректированной САУ в разомкнутом состоянии, построенная по данной методике. Параметры вариантов системы таковы, что некорректированная САУ либо неустойчива, либо имеет малый запас устойчивости, не удовлетворяющий требованиям обеспечения заданного качества регулирования. Для примера (рис. 3) значение ФЧХ в частоте снк среза некорректированной САУ в соответствии с (3.2) и (3.3) составит Lв(снк)=–4лог; Lн(снк)=1,57 лог и pнк(снк) = [Lв(снк)–Lн(снк)] · 45= =(–4–1,57) · 45=–250,6 т. е. некорректированная САУ неустойчива, т. к. |pнк(снк)|>180.
Рис. 3. К построению ЛАЧХ: некорректированной САУ LРНК,
желаемой LРЖЕЛ и коррекции LК'
При построении желаемой ЛАЧХ системы в разомкнутом состоянии в соответствии с требованием обеспечения перерегулирования по задающему воздействию не более 30 % необходим запас по фазе 45 и по амплитуде L0,75 лог (рис. 2). Это соответствует фазе pжел(с)135 или, согласно (3.2) и (3.3), соотношению
|Lв(с) – Lн(с)|3 лог. (3.4)
Следует учитывать, что требование минимального времени регулирования выполняется при максимально возможной величине частоты с скорректированной САУ, а требование максимального ослабления возмущений выполняется при минимальном уменьшении координат Lp в области низких частот (до с).
Частоты сопряжения на участках желаемых ЛАЧХ LpЖЕЛ() с различными наклонами следует выбирать с таким расчётом, чтобы передаточная функция корректирующего устройства была реализована наиболее просто, что достигается, если частоты сопряжения участков желаемой ЛАЧХ выбраны равными частотам сопряжения звеньев исходной некорректированной системы.
В примере рис. 3, если стремиться в области существенных частот (области частот с) к желаемой ЛАЧХ, близкой к оптимуму с однократным интегрированием, удобно принять частоты сопряжения участков с наклоном –1 лог/дек и –2 лог/дек, равной d, участков с наклоном –2 лог/дек и –3 лог/дек, равной ф, и частоту сопряжения 4 участков с наклонами –3 лог/дек и –4 лог/дек достаточно высокой, чтобы она не влияла на запас устойчивости и L, т. е. 4>100. Для такой конфигурации LРЖЕЛ() (см. рис. 3) L определяется координатами LРЖЕЛ при абсциссе, лежащей на середине отрезка между lgd и lgф, т. е. должно соблюдаться условие
(или lgС lgФ–0,75). (3.5)
Аналогично определяется LH(С)=1 лог (поскольку наклон LРЖЕЛ() до частоты С равен –1 лог/дек, а lg kР МАКС>1):
LВ(С)= – (lg d – lg С) – 2(lg Ф – lg d) – 3(lg C – lg Ф+1).
Тогда согласно (3.4),
|Lв(с)– Lн(с)|=|lg d –2 lg С +lg Ф – 3 – 1| 3 ,
отсюда
.
(3.6)
Частота С желаемой ЛАЧХ принимается равной наименьшему из значений, рассчитанных по (3.5) и (3.6). Например, для рис. 3: по (3.5)
lg С 1,7–0,75=0,95; по (3.6) lg С (1,54+1,7–1)/2=1,145.
Выбираем lg С=0,95 (С=100,95=8,91 с–1).
По приведенной выше методике может быть определена частота С и построена желаемая ЛАЧХ LРЖЕЛ() и иной конфигурации.
Динамическая составляющая ЛАЧХ LК!() корректирующего устройства (коэффициент kI учтен в kР МАКС) получается путем вычитания зависимости LРНК() из LРЖЕЛ(), т. е.
LК1()= LРЖЕЛ()-LРНК().
Полная ЛАЧХ корректирующего устройства
LК()=lg k1 + LК1().
В качестве схемотехнической реализации корректирующего устройства может быть рекомендован операционный усилитель, включенный по схеме рис. 4 а. Используя методику, изложенную в первой части задания, нетрудно получить передаточной функции в виде
.
(3.7)
Рис. 4. К расчету корректирующего устройства: а – принципиальная схема; б…д – ЛАЧХ при различных параметрах
При различных вариациях параметров можно реализовать различные виды коррекции:
опережающую (рис. 4 б) – при С3=0;
запаздывающую (рис. 4 в) – при С2=0;
комбинированные: (рис. 4 г) – при T3>T4>T1>T2;
(рис. 4 д) – при T1>T2>T3>T4.
Постоянные времени числителя передаточной функции (3.7) соответствуют частотам сопряжения аппроксимированной ЛАЧХ LК() с положительным переходом, т. е. наклона с –1 на 0 или с 0 на +1, а знаменатель соответствует частотам сопряжения участков LК() с отрицательным переходом, т. е. с 0 на –1 или с +1 на 0, если следовать по LК() в сторону возрастания частот.
Рассматриваемому примеру соответствует коррекция вида, указанного на рис. 4 д. Частота сопряжения 1=1/Т1; 2=1/Т2; 3=1/Т3 и 4=1/Т4 легко получаются, если взять антилогарифм (десятичный) абсцисс, приведенных на рис. 3. По (3.7) можно составить систему пяти уравнений (четыре постоянные и k1=R6/R4 c шестью неизвестными параметрами). Поэтому один из параметров необходимо определить. Рекомендуется задать значение С2=1 мкФ, а остальные параметры получить из решения системы уравнений
(3.8)
.
Рассчитать кривую переходного процесса n(t), отклонение nЗ частоты вращения вала привода от заданного значения nЗАД при приложении заданного возмущения наиболее удобно с использованием программ MatLab и System View.
Показатели
качества переходного процесса tР=tПП
– время
регулирования или время переходного
процесса;
– перерегулирование; М
–
колебательность.
1. Временем регулирования tР называется наименьший промежуток времени, по истечению которого кривая переходного процесса регулируемой величины будет (при своём дальнейшем изменении) отклонятся от установившегося значения не более чем на Δ. Примем Δ=5 % от Δn∞.
2.
Перерегулированием называется отношение
разности между максимальным и
установившимся значением регулируемой
величины к её установившемуся значению,
умноженному на 100 %:
3. Колебательность М переходного процесса характеризуется числом колебаний за время переходного процесса.