Идентификация параметра зазора при регулировании положения

СОДЕРЖАНИЕ

ВВЕДЕНИЕ

1 РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ НЕИЗМЕНЯЕМОЙ ЧАСТИ ПРИВОДА

2 СИНТЕЗ АЛГОРИТМА РЕГУЛЯТОРА СКОРОСТИ

2.1 Общие сведения

2.2 Определение дискретной частотной характеристики неизменяемой части

2.3 Определение желаемой дискретной частотной характеристики

2.4 Определение дискретной передаточной функции регулятора скорости

3 ПРОВЕДЕНИЕ ИССЛЕДОВАНИЯ МОДЕЛИ НА ЭВМ БЕЗ ВОЗМУЩЕНИЯ И ПРИ НАЛИЧИИ ВОЗМУЩЕНИЙ

3.1 Определение реакции системы без возмущений

3.2 Определение реакции системы под нагрузкой

3.3 Работа привода под нагрузкой

4 РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ ПОДСИСТЕМЫ ИДЕНТИФИКАЦИИ

ВЫВОДЫ

ПЕРЕЧЕНЬ ССЫЛОК

Реферат

Курсовая работа содержит 27 страниц, 24 рисунка, 4 источника.

Объект исследования – цифровой электропривод.

Цель работы – разработать математическую модель подсистемы идентификации зазора в технологическом перемещении рабочего органа по критерию точности.

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

ПЕРЕДАТОЧНАЯ ФУНКЦИЯ, ЦИФРОВОЙ ЭЛЕКТРОПРИВОД, СИНТЕЗ, РЕГУЛЯТОР, ИДЕНТИФИКАЦИЯ, ПЕРЕХОДНОЙ ПРОЦЕСС, ЧАСТОТНАЯ ХАРАКТЕРИСТИКА, МОМЕНТ ИНЕРЦИИ

Введение

Курсовая работа по дисциплине «Идентификация и моделирование технологических объектов» отражает основные предпосылки для выполнения идентификации параметров элементов цифровых электроприводов, управляемых ЭВМ.

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

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

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

 

1.     Разработка математической модели

неизменяемой части привода

Неизменяемая часть состоит из фиксатора, широтно-импульсного преобразователя, двигателя и датчика обратной связи.

Структурная схема цифрового электропривода представлена на рисунке 1.1

Рисунок 1.1 – Структурная схема цифрового электропривода.

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

- ПФ ШИП , (1.1)

где  (запаздывания регулятора/период дискретности).

- ПФ двигателя , (1.2)

- ПФ датчика обратной связи , (1.3)

- ПФ фиксатора . (1.4)

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

(1.5)

где . (1.6)

Дискретная ПФ неизменяемой части:

, ,

где    и т.д.

. (1.7)

В найденные формулы подставим числа, заданные в условии:

Т0=10 мс,

τр=0.2Т0=2 мс,

 мс.

Коэффициент передачи цифрового электропривода , т.к. .

Отсюда , .

Двигатель ПБВ 132 имеет характеристики:

- номинальный момент 35 Нм,

- номинальная скорость 600 об/мин,

- момент инерции якоря 0.189,

- Тя 14.2 с.

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

, .

Значит, дискретная передаточная функция цифрового электропривода будет иметь вид:

2                   Синтез алгоритма регулятора скорости 2.1 Общие сведения

Синтез регулятора базируется на двух частотных характеристиках:

-     дискретной частотной характеристике неизменяемой части привода ( ДЧХ НЧ)

-     желаемой дискретной частотной характеристике

ДЧХ параметрически оптимизируемого регулятора определяется из соотношения:

 (2.1)

где  - желаемая ДЧХ разомкнутого контура в области низких и средних частот.

 - характеризует запаздывание, связанное с дискретными операциями.

2.2  Определение дискретной частотной характеристики неизменяемой части

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

 (2.2)

 (2.3)

Заменяя в дискретной передаточной функции z на w и w на  можно перейти к дискретной частотной характеристике

 (2.4)

Переход от z к w осуществляется с помощью этого выражения.

По выражению (1.7) можно записать:

 (2.5)

,(2.6)

, (2.7)

, (2.8)

, (2.9)

В нашем случае , тогда.

Окончательно получим ДЧХ неизменяемой части в виде:

 (2.10)

2.3 Определение желаемой дискретной частотной характеристики

Тип желаемой частотной характеристики зависит от требуемого порядка астатизма. Существует несколько видов типовых ЛЧХ цифровых электроприводов (чаще рассматривается пять видов типовых ЛЧХ). Все типовые ЛЧХ имеют в области средних частот асимптоты с наклоном -20дб/дек и отличаются друг от друга лишь наклоном в низкочастотной области. В области высоких частот формы ЛЧХ приняты одинаковыми. На рисунке 2.1 приведены пять типов желаемых ЛЧХ.

 

Рисунок 2.1 – Типовые ЛЧХ цифровых электроприводов.

Передаточные функции, соответствующие приведенным на рисунке 2.1 желаемым ЛЧХ, соответственно имеют вид:

, (2.11)

, (2.12)

, (2.13)

, (2.14)

. (2.15)

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

При проектировании или исследовании приводов станков с ЧПУ должны быть известны параметры точности при максимальных значениях скоростей и ускорений. В связи с этим входное воздействие необходимо подобрать так, чтобы оно было эквивалентно заданным параметрам. В качестве типового воздействия для электроприводов станков обычно используется гармонический сигнал.

, (2.16)

, (2.17)

. (2.18).

где  - обобщенная координата, для которой заданы:

 - максимальная скорость;

 - допустимое ускорение;

 - скоростная ошибка.

Для заданных по условию курсовой работы значений , ,  определим эти параметры:

, ,

 с-1,

  м

Таким образом, эквивалентное гармоническое воздействие будет иметь вид:

 (2.19)

При этом должно выполняться условие, что ошибка будет меньше или примерно равна отношению  к модулю частотной характеристики:

, (2.20)

. (2.21)

В низкочастотной области , следовательно:

 (2.22)

 дб

Определим,  - отношение допускаемого ускорения к ошибке. Для создания запаса кривую поднимают на 3 дб вверх  с-1. Без подъема кривой на 3 дб вверх  с-1.

Коэффициент передачи системы:

 дб.

, (2.23)

где М – коэффициент колебательности, в моем случае М=1.2. Тогда:

 с-1

Частота, при которой участок высокочастотной ЛЧХ переходит в наклон 0 дб/дек определится:

 с-1 (2.24).

Для обеспечения требуемой точности и коэффициента колебательности системы, а в конечном итоге, требуемой устойчивости, необходимо обязательное соблюдение условия, чтобы участок ЛЧХ, соответствующий наклону 0 дб/дек лежал не ниже значения, вычисленного по следующей формуле:

 дб.

Таким образом, результирующую ЛАЧХ мы получим после подъема исходной ЛАЧХ на 4.682 дб, что показано на рисунке 2.3

Определим коэффициент передачи “новой” системы.

 дб.

 - такое значение без подъема будет у ЛАЧХ при наклоне 0 дб/дек

Рисунок 2.2 – Желаемая ЛАЧХ цифрового электропривода.

  2.4 Определение дискретной передаточной функции регулятора скорости

Желаемая дискретная передаточная функция представлена выражением (предпосылки для ее записи представлены в предыдущем подпункте) (2.25)

(2.25)

 (2.26)

. (2.27)

Согласно формулам (2.10) и (2.25) находим параметр, характеризующий запаздывание:

 (2.28)

– ДПФ регулятора скорости.

Коэффициенты ДПФ имеют следующие значения (рассчитаны в MathCAD):

  

  

 

3       Оценка качества регулирования привода в

переходных и установившихся режимах

3.1 Определение реакции системы на единичное воздействие

Определение реакции системы на единичное воздействие производится с помощью программы MATLAB 6.

При использовании MATLAB 6 производится имитация или моделирование процессов, происходящих в цифровом электроприводе. Каждый элемент представлен своей передаточной функцией, а результат моделирования – переходный процесс выводится на экран осциллографа (Scope).

Рисунок 3.1 – Структурная схема ЦЭП, полученная на MATLAB.

Рисунок 3.2 – Переходной процесс, полученный на MATLAB по расчётным коэффициентам регулятора.

Коэффициент колебательности у такой системы:

.

Время переходного процесса:с или 100 мс.

Перерегулирование .

Система удовлетворяет заданным требованиям.

  3.2 Определение реакции системы на линейно-возрастающее воздействие

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

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

На выходе системы имеем переходной процесс:

t, c

W, м/с

Рисунок 3.4 – Переходной процесс при подаче линейно-нарастающего воздействия.

  3.3 Работа привода под нагрузкой

Кроме исследования работы привода без нагрузки целесообразно провести анализ его работы под действием нагрузки, которая представляет собой момент сопротивления (или статический ток ).

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

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

 (3.4)

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

 (3.5)

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

Таким образом, упрощенно структурную схему можно представить в виде, показанном на рисунке (3.6)

Рисунок 3.6-Схема привода для определения реакции его на нагрузку

Приведенная выше схема реализуется на MATLABе. В результате полученного переходного процесса можно сказать, что при нагрузке, близкой к номинальной в приводе в момент введения этой нагрузки ( определяется “дискретными часами” ) наблюдается просадка скорости, которая затем устраняется благодаря влиянию регулятора. Переходный процесс для привода при включении нагрузки представлен ниже.

Рисунок 3.7 – Структурная схема на MATLAB с набросом нагрузки.

	W, м/с

 

t, c

Рисунок 3.8 – Переходной процесс при набросе нагрузки.