grandov.ru страница 1
скачать файл
Цель работы – определение областей возможных значений двух параметров, обеспечивающих устойчивую работу САР угловой скорости паровой турбины, методом D-разбиения; экспериментальная проверка теоретических результатов с помощью компьютерного моделирования.
Теоретическая часть

Описание работы САР угловой скорости паровой турбины

Функциональная схема (блок-схема) системы автоматического регулирования угловой скорости паровой турбины (Рис. 1.1) показана на Рис. 1.2


Рис. 1.1 Система автоматического регулирования угловой скорости паровой турбины


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


На Рис. 1.1 цифрами обозначены:

  1. Турбина

  2. Грузы

  3. Рычаг

  4. Задвижка

  5. Муфта

При отклонении угловой скорости турбины (Рис. 1.1) от заданного значения меняется центробежная сила грузов, в связи с чем изменяется положение муфты, которая рычагом приводит в действие исполнительный механизм – задвижку.

Данная САР является системой прямого действия.


Краткий вывод дифференциальных уравнений звеньев рассматриваемой САР

1) Регулируемый объект (паровая турбина)

Активный момент на валу турбины

, 

где – коэффициент, зависящий от конструкции паровой турбины, – величина открытия задвижки, – скорость движения пара.

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

, 

где – установившийся активный момент.

Уравнение  приведем к виду

. 

Момент сопротивления на валу турбины

, 

где – установившийся момент сопротивления на валу, - момент от изменения нагрузки.

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

, 

где – приведенный к валу турбины момент инерции всех вращающихся масс.

Подставляя в уравнение  зависимости  и  и учитывая то, что в установившемся режиме , получим уравнение турбины в приращениях

. 

Номинальными значениями для турбины будут , , . С учетом этого  примет вид



. 

Вводя следующие обозначения , , , , получим



Уравнение  характеризует турбину как систему первого порядка с передаточной функцией



, 

где - постоянная времени, - оператор Лапласа.

2) Чувствительный элемент (центробежный маятник)

Центробежный маятник показан на Рис. 1.3


Рис. 1.3 Центробежный маятник


На Рис. 1.3 цифрами обозначены

  1. Вал

  2. Демпфер

  3. Муфта

  4. Шары

  5. Пружина

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

Для составления дифференциального уравнения центробежного маятника используется уравнение Лагранжа II рода

, 

где – обобщенная координата (в данном случае в качестве выбирается выходная координата – перемещение муфты ), – кинетическая энергия, – обобщенная сила, которая в общем случае системы с одной степенью свободы имеет вид

, 

где – потенциальная энергия системы, - функция рассеяния энергии, - внешняя сила, приложенная к динамической системе. Считается, что .

Кинетическая энергия системы равна



где – кинетическая энергия шаров во вращательном движении вокруг оси , – кинетическая энергия шаров во вращательном движении относительно точек , – кинетическая энергия масс в поступательном движении вдоль оси .

Указанные составляющие кинетической энергии



;

;

,

где – масса частей, участвующих в поступательном движении вдоль оси .

Потенциальная энергия системы



где – потенциальная энергия масс, движущихся параллельно оси , – потенциальная энергия шаров, – потенциальная энергия пружины.

Указанные составляющие потенциальной энергии

;

;

,

где – коэффициент жесткости пружины, сила предварительного поджатия пружины.

Функция рассеяния энергии имеет вид

,

где – коэффициент вязкого трения.

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

, 

где – постоянная времени муфты маятника, – постоянная времени демпфера, – коэффициент статизма.

Уравнение  характеризует центробежный маятник как звено второго порядка – колебательное звено, передаточная функция которого имеет вид

. 

В МВТУ-3.6 коэффициенты колебательного звена имеет иной вид

, 

где – коэффициент усиления, – постоянная времени, – коэффициент относительного демпфирования.

Связь между коэффициентами в  с коэффициентами в  следующая



3) Исполнительный механизм (задвижка)

Уравнение движения задвижки как безынерционного звена



, 

где – передаточный коэффициент, – безразмерное перемещение точки А (Рис. 1.1), соответствующее стационарному режиму, – малое отклонение той же точки относительно .

4) Другие элементы в цепи обратной связи

Угловые скорости вала центробежного регулятора и вала турбины связаны между собой соотношением



Или, разделив на ,



, 

где – передаточное число связи вала турбины с валом регулятора.

В рассматриваемой системе автоматического регулирования использована жесткая отрицательная обратная связь. Рассматривая рычаг 3 (Рис. 1.3) как звено нулевого порядка, представим его уравнение в виде

. 
Теоретический расчет областей устойчивости методом D-разбиения

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

Основываясь на функциональной схеме, структурная схема регулирования будет выглядеть следующим образом


Рис. 1.4 Структурная схема регулирования

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



.

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



. 

Отсюда следует характеристическое уравнение, соответствующее исследуемой САР,



или


. 

Подставив в выражение  и разделив его на два уравнения, получим



Откуда найдем



Задавая значения в интервале по формулам  можно построить кривую D-разбиения в плоскости параметров и . Каждая точка кривой D-разбиения определяет такие значения параметров, при которых характеристическое уравнение имеет чисто мнимые корни. При переходе через кривую D-разбиения по крайней мере один вещественный корень или пара комплексно-сопряженных корней характеристического уравнения пересекает мнимую ось. Штриховка D-кривой производится слева по направлению движения при движении от к , если главный определитель системы  , и справа, если .

Криволинейная граница D-разбиения должна быть дополнена “особыми прямыми”. Такие уравнения получаются подстановкой в уравнение  частот , и других , при которых решение системы  имеет вид . В рассматриваемой системе отсутствуют.

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



.

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

. 

Уравнение особой прямой при можно получить, приравняв нулю коэффициент при старшем члене характеристического уравнения 



;

или


. 

Штриховка “особых прямых” наносится таким образом, чтобы одновременно заштрихованные или одновременно незаштрихованные стороны D-кривой и прямой располагались навстречу друг другу.



Практическая часть

Исходные данные

Неизменяемые параметры

с, с, с, , .

варьируемые параметры

, .
На Рис. 1.5 показана структурная схема системы автоматического регулирования угловой скорости паровой турбины, реализованная в МВТУ-3.6.

Рис. 1.5 Структурная схема


Кривая D-разбиения, построенная по формулам  при , а также две особые прямые, полученные согласно формулам  и , показаны на Рис. 1.6

Рис. 1.6 D-разбиение


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

Коэффициенты колебательного звена (Рис. 1.5) определяются по формулам , где с для точек , и



; c; .

Получим графики переходных процессов для точек , и , показанных на Рис. 1.6.




  1. Точка : c2, .

Рис. 1.7 График переходного процесса в точке


Переходный процесс асимптотически устойчивый.


  1. Точка : c2, .

Рис. 1.8 График переходного процесса в точке


Переходный процесс неасимптотически устойчивый.


  1. Точка : c2, .

Рис. 1.9 График переходного процесса в точке


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






0.90

0.6452

0.95

0.681

1.00

0.7169

1.05

0.7527

1.10

0.7886


Рис. 1.10 Рабочая характеристика паровой турбины

Как видно из графика (Рис. 1.10), рабочая характеристика линейная, что и ожидалось согласно дифференциальному уравнению .

Заключение. В данной лабораторной работе методом D-разбиения были определены области возможных значений двух параметров – коэффициента усиления и муфты маятника , обеспечивающих устойчивую работу САР угловой скорости паровой турбины. Построением в МВТУ-3.6 графиков переходных процессов для трех точек из плоскости D-разбиения была подтверждена корректность определения данных областей. Была также построена рабочая характеристика паровой турбины, которая оказалась линейной.

Данная САР является системой со статическим регулированием, так как равновесие системы имеет место при различных значениях регулируемой величины и каждому значению регулируемой величины соответствует единственное определенное положение регулирующего элемента.





скачать файл



Смотрите также:
Теоретическая часть Описание работы сар угловой скорости паровой турбины
94.64kb.
Основы катехизаторской работы Вероника Карасевич 3 курс
317.92kb.
Кинетическое описание dutpase фермента метаболизма нуклеотидов бактерии Escherichia coli
171.58kb.
Теплотехнические устройства
231.3kb.
Теоретическая часть Основные термины и положения в windows для работы с файлами и папками
68.38kb.
8. теоретическая механика
184.54kb.
Глава Теоретическая часть
951.8kb.
Программа учебной дисциплины «теоретическая механика» Специальность: 130102 «Технологии геологической разведки»
157.25kb.
A и b положительные константы. Найдите уравнение кривой, описывающей траекторию частицы, и изобразите ее график. Определите зависимости от времени абсолютных величин скорости и ускорения частицы
107.08kb.
Изучению курса предшествуют следующие дисциплины: курс основного иностранного языка, история языка, риторика, теоретическая фонетика, теоретическая грамматика, лексикология
113.57kb.
Ченцов Е. Л. (1), Клочкова В. Г. (1), Мирошниченко А. С. (2), Мансет Н
22.54kb.
Каждое описание содержит основные теоретические положения, относящиеся к данному разделу дисциплины, описание лабораторной установки, порядок проведения работы и перечень контрольных вопросов
528.8kb.