Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

1603

.pdf
Скачиваний:
0
Добавлен:
15.11.2022
Размер:
795.19 Кб
Скачать

Maple и MatLab для решения задач моделирования электроприводов используются явно недостаточно.

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

7.2.Программа комплексного анализа структурных схем

7.2.1.Общие сведения о программе KASPER

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

Ниже приведены основные возможности программы КАSPER:

1.Ввод структурной схемы объекта САУ в удобной для оператора форме.

2.Автоматическое преобразование схемы в таблицу соединений передаточных функций и ее просмотр.

3.Хранение отдельных передаточных функций на диске с целью доступа, сохранения, редактирования и их повторного использования.

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

5.Анализ прохождения и преобразования сигнала в системе путем установки нескольких промежуточных выходов.

6.Анализ поведения системы при ненулевых условиях путем установки дополнительных входов.

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

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

62

9. Получение частотных (амплитудной и фазовой) и временных (переходной) характеристик исследуемых объектов по результирующей передаточной функции.

10. Запись на диск после анализа результирующей передаточной функции для последующего использования (в качестве исходной передаточной функции) или для более глубокого ее анализа с выводом результатов на печать с помощью программы dinprf или аналогичных.

11.Просмотр результирующих графиков на одном экране с целью их сравнения.

12.Изменение структурной схемы или (и) исходных данных (передаточных функций).

13.Возможность использования программы в программных комплексах для анализа динамических свойств ранее введенной структурной схемы системы автоматического управления (САУ), параметры которой заложены в прикладной программе.

7.2.2. Математические методы, применяемые в программе

Преобразование структурных схем.

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

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

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

Анализ эквивалентной передаточной функции.

Анализ полученных передаточных функций возможно проводить в частотных (амплитудной АЧХ и фазовой ФЧХ характеристик) и временной (переходные характеристики) областях.

63

Расчет частотных характеристик АЧХ и ФЧХ основан на замене в передаточной функции (1.18) оператора p на jw и расчете амплитудно-фазовой характеристики.

Расчет переходного процесса заключается в обратном преобразовании Лапласа функции W(p)/p.

7.2.3. Построение программы, входные и выходные данные

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

ctc структурная схема во внутреннем формате; ott основная таблица типов;

prf передаточная функция.

Кодирование информации выполняется программой автоматически. Следует отметить, что по объему файлов можно судить об объеме схемы ctc, о числе полей таблицы ott, о порядке передаточной функции prf.

7.2.4. Описание логической структуры программы

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

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

Модули специального применения разработаны непосредственно для программы KASPER и выполняют следующие функции (табл. 7.1).

 

Таблица 7.1

Модуль

Назначение

EDITCX

Редактор структурных схем и создание таблицы расчета

RUNCX

Вычисление эквивалентной передаточной функции по таблице

 

расчета

UTISIN

Подготовка таблицы расчета для RUNCX и анализ результи-

 

рующих передаточных функций

 

Нахождение корней многочлена n-го порядка, расчет формулы

SQRPOLIN

переходного процесса, функции расчета переходного процесса и

 

частотных характеристик

На рис.7.1 показаны блоки программы в порядке проведения анализа динамических свойств объекта управления.

64

Выбор необходимого при анализе блока производится оператором через меню. Все связанные с оператором блоки имеют выход в меню.

7.2.5. Выполнение программы

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

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

KASPER

Редактор структурных схем

Анализ структуры схемы

 

 

 

Меню

 

 

 

 

 

 

 

Редактор таблицы

 

Преобразования

типов

 

структурной схемы

 

 

 

 

 

 

 

Выход

 

 

 

 

 

 

 

 

 

 

Анализ пере-

 

 

Построение

даточной

 

 

графиков

функции

 

 

расчета

 

 

 

 

 

Временной

 

Частотный

 

 

 

Рис.7.1. Функциональная схема программы KASPER

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

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

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

65

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

Для выбранных пользователем выходов может быть произведен расчет переходного процесса и частотный анализ.

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

Uo

p 1T

K +

Tp 1

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

Применение программы представлено в примере 7.1.

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

66

Uя

iя R яц

Lяц

 

diя

 

e,

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

E k w1,

 

 

 

 

 

 

 

 

 

 

M1

k iя ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dw1

(7.1)

M

 

M

 

 

M

 

 

 

J

 

 

,

1

c1

12

 

 

1

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

M

 

J

 

dw 2

,

 

 

12

c2

2

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M12

c12 (

1

 

2 ).

 

 

 

 

 

Применив к уравнениям (7.1) преобразование Лапласа, составляется структурная схема, изображенная на рис.7.2.

u

п

 

 

 

д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

d

M

 

m

 

 

 

Н

 

1

сп

2

w

 

 

W

 

Рис.7.2. Пример ввода в ЭВМ исследуемой структурной схемы

Передаточные функции звеньев, представленных на рис.7.2, следующие:

Wп (р)

 

кп

 

 

 

 

31

 

 

тиристорный преобразователь;

 

Тп

р

1

 

0,01р

1

 

 

 

 

 

 

 

 

 

 

Wд

(р)

 

1 R а кд

 

1 0,3

0,758

 

4,4

 

двигатель;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Та

р

1

 

 

 

0,053р

1

 

 

0,053р 1

 

 

 

 

 

 

 

 

 

W1

(p)

1

 

 

 

1

 

 

 

 

первая масса;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J1p

1,48p

 

 

 

 

 

 

 

 

Wc

(p)

 

12 p

 

c12

 

 

 

1,56p

252

упругая механическая связь;

 

 

 

 

p

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W2

(p)

1

 

 

 

1

 

 

 

 

вторая масса;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J2p

 

2,5p

 

 

 

 

 

 

 

 

67

Wd (p)

1

1

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

 

 

 

кд

0,758

 

 

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

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

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

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

Примечание:

Основная таблица

 

Вход u (напряжение задания системы)

 

 

 

 

подключен к узлу У1.

ФУН

ВХОД

 

ВЫХОД

u

У0

 

У1

Передаточная функция п (преобразова-

п

У1

 

С1

тель) подключена входом к узлу У1,а выходом

д

С1

 

У2

М

У2

 

0

к сумматору С1.

 

 

У2

 

С2

Передаточная функция д (двигатель)

1

С2

 

У3

d

У3

 

-1

подключена входом к сумматору С1, а выхо-

 

W

У3

 

0

дом к узлу У2.

 

У3

 

С3

Выход М (момент двигателя) соединен с

c

С3

 

У4

m

У4

 

0

узлом У2.

 

 

У4

 

-2

Узел У2 соединен с сумматором С2.

 

У4

 

С4

2

С4

 

У5

Передаточная функция 1 (первая масса)

 

 

У5

 

-3

подключена входом к сумматору С2, а выхо-

w

У5

 

0

дом к узлу У3.

н

У0

 

-2

Н

У0

 

-4

 

 

Передаточная функция d (кд) подключена

Рис.7.3. Таблица расчета

входом к узлу У1, а выходом к сумматору С1 с отрицательным знаком.

Выход w (частота вращения вала двигателя) соединен с узлом У5.

Узел У3 соединен с сумматором С3. Передаточная функция с (жесткая связь

масс) подключена входом к сумматору С3, а выходом к узлу У4.

Выход m (момент упругой деформации) соединен с узлом У4.

68

Узел У4 соединен с сумматором С2 с отрицательным знаком. Узел У4 соединен с сумматором С4.

Передаточная функция 2 (вторая масса) подключена входом к сумматору С4, а выходом к узлу У5.

Узел У5 соединен с сумматором С3 с отрицательным знаком. Выход w (частота вращения рабочего органа) соединен с узлом У5.

Вход н (момент нагрузки первой массы) поступает на сумматор С2 с отрицательным знаком.

Вход Н (момент нагрузки второй массы) поступает на сумматор С4 с отрицательным знаком.

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

Для выхода двигателя на номинальную частоту вращения подадим сигнал задания u=7 B. Результаты исследования динамики системы по управлению приведены в табл.7.2.

 

 

Связь перед. функции

 

 

 

 

 

Таблица типов

 

 

 

 

 

 

 

 

 

 

 

 

 

D:\PRIKLPRG\BDW\PRDVIG.PRF

ФУН

ТИП

 

 

Выбрать другой файл

 

 

н

Ввод Пер Фун

 

 

Переименовать

 

 

Н

Ввод Пер Фун

 

 

Редактировать

 

 

u

Ввод Пер Фун

 

 

 

Выход

 

 

п

Ввод Пер Фун

 

 

 

 

 

 

 

д

Ввод Пер Фун

 

 

 

 

 

 

 

1

Ввод Пер Фун

 

 

 

 

 

 

 

d

Ввод Пер Фун

 

 

 

 

 

 

 

c

Ввод Пер Фун

 

 

Передаточная функция Wд(p)

 

 

 

 

 

 

2

Ввод Пер Фун

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Числитель

Знаменатель

 

 

 

 

 

 

b(1) = 0

a(1) = 0.053

 

 

 

 

 

 

b(0) = 4.4

a(0) = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.7.4. Формирование таблицы типов

Расчет Выход W

Рис.7.5. Окно расчета

 

 

 

Характеристики

 

 

 

Выход

Расчет

 

 

 

.

Рек.зн.

Точек

Начальн.зн.

Конечн.зн

М

[X] h(t)

0

 

0

 

3.4

80

 

[X] F(t)

0.001

0

 

49

50

W

[X] h(t)

0

 

0

 

3.4

80

 

[X] F(t)

0.001

0

 

49

50

m

[X] h(t)

0

 

0

 

3.4

80

 

[X] F(t)

0.001

69

 

49

50

 

0

 

w

[X] h(t)

0

 

0

 

3.4

80

 

[X] F(t)

0.001

0

 

49

50

Расчет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.7.6. Выбор рассчитываемых характеристик

Просмотр просчитанных графиков производится путем выбора номеров характеристик (рис.7.7).

Выбраны графики 1, 5, 6 (рис.7.7), представляющие переходную h(t), ампли- тудно-частотную A( ) и фазочастотную ( ) характеристики момента двигателя (рис.7.8).

 

 

 

 

 

Таблица 7.2

Передаточная функция

Переходный процесс

числитель

знаменатель

 

 

 

 

 

Частота вращения рабочего органа

 

 

 

a6 = 1

 

 

 

 

 

a5 = 282

164.486

+

 

 

 

a4 = 21900

1.1 10-14

e-161t +

 

 

 

a3 = 4.36 105

0.0057

e-100t +

 

b2 = 7.6

105

a2 = 7.44 106

11.07777

e-1.65t cos( 17t

29.9 ) +

b1 = 2.45

108

a1 = 8.53 107

186.421

e-1.62t

 

b0 = 1.98

1010

a0 = 1.2 108

12.337

e-15.6t

 

 

 

Частота вращения двигателя

 

 

 

a5 = 1

164.486 +

 

 

 

a4 = 120.5

1.4773

+ e-100t +

 

 

 

a3 = 2431

21.493

+ e-1.65t cos(17t

143 ) +

b2 = 1.22

106

a2 = 42830

191.33

+ e-1.62t +

 

b1 = 7.59

105

a1 = 5.23 105

45.4

+ e-15.6t

 

b0 = 1.23

108

a0 = 7.46 105

 

 

 

 

 

Момент упругой деформации

 

 

 

a5 = 1

1.41277 + e-100t +

 

 

 

a4 = 120.5

 

 

 

a3 = 2431

475.218 + e-1.65t cos(17t

234 ) +

b2 = 1.9

106

a2 = 42830

755.839 + e-1.62t +

 

b1 = 3.07

108

a1 = 5.23 105

481.8

+ e-15.6t

 

 

 

 

70

 

 

b0 = 0

a0 = 7.46 105

 

 

 

 

 

 

 

 

 

 

 

 

Момент двигателя

 

 

 

 

 

 

 

 

a5 = 1

220.052 + e-100t +

 

 

 

 

 

 

a4 = 120.5

 

 

 

 

b3

= 1.8 10

a3

= 2431

 

97.067

+ e-1.65t cos(17t

7.4 ) +

 

b2

= 3.02 106

a2

= 42830

 

1215.09 + e-1.62t +

 

 

 

 

b1

= 4.88 105

a1

= 5.23 105

 

1531.4

+ e-15.6t

 

 

 

 

b0

= 0

a0

= 7.46 105

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Результат

 

 

 

 

 

 

 

 

 

 

Выход

 

 

 

 

f(w)

 

 

 

 

 

 

 

 

 

h(t)

A(w)

 

 

 

 

График

 

 

 

M

1

5

6

 

 

 

 

 

W

2

7

8

 

 

 

 

 

 

 

 

 

 

Первый

 

1

 

 

 

 

 

m

3

9

10

 

 

Второй

 

12

 

 

 

 

 

 

 

w

4

11

12

 

 

Просмотр

Последовательно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.7.7. Выбор выводимых графиков

a)

б)

71

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]