Система автоматического регулирование температуры теплоносителя зерносушилки

Курсовой проект - Компьютеры, программирование

Другие курсовые по предмету Компьютеры, программирование

Скачать Бесплатно!
Для того чтобы скачать эту работу.
1. Пожалуйста введите слова с картинки:

2. И нажмите на эту кнопку.
закрыть



, входящих в уравнения отражен в описании САР. Параметры Тс, Тя, Тм, Тя - постоянные времени, с; kc, kf, kA, ky, k0 - коэффициенты передачи. Выходной сигнал усилителя ограничен уровнем 220 В. Заданная температура теплоносителя θС = 120 С.

 

3. Построение структурной схемы нескорректированной САУ и определение передаточных функций её звеньев

 

1)Найдем передаточную функцию объекта управления (камера смешивания)

 

 

 

 

2) Найдем передаточную функцию электродвигателя

 

 

 

3) Найдем передаточную функцию усилителя

 

 

 

4) Найдем передаточную функцию датчика температуры

 

 

 

5) Сравнивающий орган

 

 

Итоговая схема САР

 

 

 

4. Оценка точности и анализ устойчивости исходной системы

 

Приведем схему САР к единичной обратной связи

 

 

 

Разложим звенона более простые звенья, для этого определим колебательное оно или апериодическое

Характеристическое уравнение знаменателя

 

 

 

Если

 

,то

 

звено апериодическое, иначе колебательное

 

Тм2 >4ТмТя

0,0182 > 4∙0,018∙0,002

0,000324 > 0,000144

 

звено апериодическое.

 

Т1 = 0,0157

Т2 = 0,0023

 

Оценка точности

а) ошибка по задающему воздействию

 

 

б) ошибка по возмущению

 

 

Оценка устойчивости

Для оценки устойчивости используем пакет программ MATLAB и его модуль Control System Toolbox

Задаем передаточную функцию разомкнутой системы

 

>> W=tf([4.5],[0.000061 0.000467 3.75 3.31 1 0])

 

Transfer function:

 

>> step(W);grid on

 

 

Строим ЛАЧХ и ЛФЧХ.

Определяем запасы по фазе и амплитуде

>> margin(W);grid on

 

 

Запас по амплитуде достаточен Gm=14,1 dB и входит в рекомендуемый запас 6-20dB.

Запас по фазе достаточен Pm=40,5 deg и входит в рекомендуемый запас 30-60deg.

Определяем устойчивость по критерию Найквиста

>> nyquist(W);grid on

 

 

 

 

Охватывает точку (-j;0), значит по критерию Найквиста разомкнутая система не устойчива.

Построим замкнутую систему н. основе разомкнутой

>> Z=feedback(W,1)

Transfer function:

 

>> step(Z);grid on

 

 

 

Строим ЛАЧХ и ЛФЧХ.

Определяем запасы по фазе и амплитуде

>> margin(Z);grid on

 

 

Определяем устойчивость по критерию Найквиста

>> nyquist(Z);grid on

 

 

 

Не охватывает точку (-j;0), значит по критерию Найквиста замкнутая система устойчива. Годограф замкнутой системы проходит дальше от точки (-j;0), чем годограф разомкнутой системы, значит замкнутая система более устойчива.

 

5. Построение желаемой ЛАЧХ

 

Построим ЛАЧХ исходной нескорректированной системы

ν =1, k=4.5, ТД=1.3, TС=1.3, T1=0.0157, T2=0.0023

20lgK=20lg4.5=13.06 дБ

L(ω)=20lgA(ω)=20lg|W(jω)|=20lg4.5 -20lg

-20lg-20lg-20lg

lg w0= -0.11

lg w1= -0.11

lg w2= 1.8

lg w3= 2.6

 

Первая точка ЛАЧХ 20lgK=13.06 дБ

Начальный наклон при ν =1 равен

Звенья апериодические и каждое из них изменяет наклон на

Для построения желаемой ЛАЧХ необходимо найти желаемый передаточный коэффициент:

Построение желаемой ЛАЧХ

Находим частоту среза по графику 5.24

При σ=15% Рмах=1

 

ωcp ≥ ωc1

ωcp = 100

 

Строим на миллиметровке ЛАЧХ исходной и желаемой системы, находим ЛАЧХ корректирующего устройства как результата вычитания ЛАЧХ исходной системы из ЛАЧХ желаемой.

По графику желаемой ЛАЧХ находим

 

lgKж=1.9 дБ

 

 

6. Оценка показателей качества желаемой системы

 

Записываем передаточную функцию желаемой ЛАЧХ и вычисляем запасы по фазе и амплитуде, устойчивость по критерию Найквиста в пакете MATLAB.

 

lgω4 = 2,6

ω4= 398c-1

 

Transfer function:

 

>> step(W);grid on

 

 

Строим ЛАЧХ и ЛФЧХ.

Определяем запасы по фазе и амплитуде

>> margin(W);grid on

 

 

Построенная MATLAB ЛАЧХ совпадает с построенной на миллиметровке, запасы по фазе и амплитуде улучшились по сравнению с исходной MATLAB ЛАЧХ

Определяем устойчивость по критерию Найквиста

>> nyquist(W);grid on

Определяем устойчивость по критерию Найквиста

>> nyquist(W);grid on

 

 

Диаграмма желаемой системы проходит дальше от точки (-j;0), чем исходной, поэтому она более устойчива.

 

7. Синтез последовательного корректирующего устройства

 

В результате вычислений в предыдущем пункте мы получили исходную и желаемую ЛАЧХ. ЛАЧХ корректирующего устройства получается в результате вычитания

 

 

По графику строим передаточную функцию корректирующего устройства.

 

 

По графику ЛАЧХ корректирующего устройства на миллиметровке находим

 

20lgKку=28 дБ

lgKку=1.4 дБ

Kку=25.1

>> W=tf([0.65,43.44,65.5,25.1],[0.000000011,0.000014,0.0067,1])

 

 

Transfer function:

 

0.65 s^3 + 43.44 s^2 + 65.5 s + 25.1

1.1e-008 s^3 + 1.4e-005 s^2 + 0.0067 s + 1

 

Отобразим её

step(W);grid on

Строим ЛАЧХ и ЛФЧХ.

Определяем запасы по фазе и амплитуде

 

 

Строим ЛАЧХ и ЛФЧХ.

Определяем запасы по фазе и амплитуде

 

 

 

Построенная MATLAB ЛАЧХ КУ совпадает с построенной на миллиметровке

 

8. Реал

s