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

книги из ГПНТБ / Корытин А.М. Синтез автоматизированного электропривода на аналоговых и цифровых вычислительных машинах

.pdf
Скачиваний:
6
Добавлен:
24.10.2023
Размер:
9.72 Mб
Скачать

ный алгоритм анализа, включающий в себя (5-5), (5-7), (5-10) и "(5-11), в виде

Х"П—Х"(П_Х)

 

- f - kx(x'n

J C ' ( n - i ) )

+

 

 

 

 

m

 

 

 

klx'n + kîx"(n.l)+

Б

x*t (zi n )

 

+

 

г ѵ т т г а

 

 

( 5 4 2 >

Обобщенный

алгоритм

синтеза

включает

(5-3), (5-4),

(5-6):

 

 

 

 

 

 

 

 

m—I

 

 

 

 

 

f i

X'qn = Х' (x"n)

-\- ^

X*i

fan)

-\- T*R

3 (x"n)

* •

 

 

 

 

 

 

 

(5-13)

В (5-13) приняты те же обозначения переменных и

параметров, что и в

(4-25). Индекс «п» в двух последних

соотношениях определяет номер расчетной точки. Алго­ ритм (5-12) требует для вычислений (/и-И) подпро­ грамм определения функций по аргументу, для хранения текущей информации (m+ 4) ячейки, для хранения ко­ эффициентов 4 ячейки. Реализация (5-13) требует +2) подпрограмм для нахождения функций по аргу­ менту, (иг+ 4) ячейки для хранения текущей информа­ ции.

Располагая обобщенными алгоритмами, можно в об­ щем виде выполнить оценку погрешности, вносимую ре­ комендуемым методом синтеза. Анализ алгоритма (5-12) показывает, что он не гарантирует точности в одной из фиксированных точек шага іМ (например, в начале, кон­ це или середине интервала). Так, во втором члене числи­ теля принята переменная для начала интервала син­ теза, в знаменателе постоянная времени выбирается по той же переменной в начале №.

Один из способов повышения точности состоит в вы­ боре минимально возможного интервала синтеза. Однако в этом случае возможны излишние затраты машинного времени. Выбрать оптимальное по величине погрешности значение шага синтеза можно путем последовательного деления исходного шага St. Деление осуществляется до тех пор, пока два смежных решения не дадут погрешно­ сти меньше допустимой.

Второй способ состоит в приведении точности к концу интервала синтеза. В качестве примера рассмотрим

151

(5-12) для kx=èi = 0. Известно, что л:"71 = д;"(П _і)-ЬА*"П ,

для вырожденного соотношения (5-12)

* " « = * " ( « - , ) ( і + 7*п .я (х»( я .і)) ) +

Sx*t(zt n )'

В(5-14) второй член правой части представляет со­

бой Ave"п. Заменим в числителе

этого члена я"(П -і) через

х"п и Асс"п и выразим

приращение в виде (5-12). Тогда

у.п

Ѵп

/ I I

ЬМ.

і_

 

7 f e r + ^ .

к , . - , , ) " -

<5 -1 5 >

В (5-15)

необходимо

учитывать

знак

коэффициента

кг. Постоянная времени

Т*я.3 определяется

для

начала

интервала,

что вносит

погрешность. Для уточнения ее

следует приближенно найти величину х"п.

Здесь

может

быть рекомендован

следующий путь. Исходя из того, что

At выбиралось

для

участка

исходной характеристики,

изменяющегося

монотонно,

нельзя

ожидать

(кроме на­

чальных участков синтеза) существенных различий меж­ ду Д*"(п-і) и Ах"п. Поэтому в первом приближении мож­ но считать, что

^*и.з(л;'/(п-і)+АІХ//(п-і)) ~ Г'і : и .з (#"п) •

Если при этом иметь в виду, что характеристики ап­ паратов и машин не имеют прерывистых участков, то рекомендуемые уточнения обеспечивают вполне удовлет­ ворительную точность синтеза.

Погрешность в алгоритмы (5 - 12), (6-13) вносит за­ мена дифференциала функции ее приращением. Наибо­ лее удачно эту погрешность можно уменьшить путем вы­ бора оптимального шага синтеза при последовательном делении его исходного значения и оценке погрешности смежных решений. Следует иметь в виду то обстоятель­ ство, что синтезированные характеристики могут быть

152

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

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

5-2. ОСОБЕННОСТИ ВВЕДЕНИЯ ФУНКЦИЙ ПРИ СИНТЕЗЕ ПАРАМЕТРОВ АВТОМАТИЗИРОВАННОГО ЭЛЕКТРОПРИВОДА

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

Рис. 5-1. К введению 'функций в ЦВМ.

В ходе проектирования и теоретических исследований приходит­ ся пользоваться усредненными характеристиками электрических ма­ шин и приборов, приводимыми в каталогах. Реальные характеристи-

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

Рассмотрим функцию, заданную графически (рис. 5-1). Извест­ но, что максимальная величина аргумента, вычисленная в процессе

синтеза, не может

быть больше л-,-, a л-<.Ѵо существовать не может.

Аппроксимируем

график

отрезками

прямых и

зафиксируем

узлы

аппроксимации -по абсциссам хо, XI,

хг,

 

" ординатам j/o,

Уи </2, • • -, (/(.--м-

Введем

значения

коэффициентов

линейной

интер­

поляции на отрезках аппроксимации

h0. Л,, Л2, .

. . ,

/ І ( І - І ) . Д Л Я

каж­

дого участка получаем общий алгоритм, отличающийся только пара­

метрами узла аппроксимации. Если х лежит

в интервале Xk—-Ки+і),

то

аргумент определяется

соотношением

 

 

 

 

 

 

 

 

 

y=yii+hh(x—xk).

 

 

 

(5-16)

 

ла,

Решение сводится к логической операции нахождения

интерва­

в котором лежит х, и вычислению ординаты у по '(5-16). Усло­

вимся

хранить аргумент в ячейке а и туда же записывать

величи­

ну

у.

Хранение

узлов аппроксимации осуществляется

в

ячейках

ß,

ß + 1 ,

ß + 2 , ...,

ß+(t — 1) .

Ординаты ' узлов

хранятся

в

ячейках

у,,

Y - H ,

У+2,

у+

(і—1),

а в ячейках б, 0-И, Ô+2, ...,

Ö-H(i—I)

записаны значения коэффициентов линейной интерполяции.

 

 

 

Лучший вариант программы

описывается

оператором

 

 

 

 

 

3 f

14

it

2 f 4-1

4f 4-2 j.3

 

(5-17)

 

 

M « > P )

m [ O P + ( t - l ) M i f ( - l . 0 - с т о п .

 

Оператор построен на принципе .модификации команд. Каждый из операторов (логический и вычислительный) подвергается модифи­ кации операцией вычитания команд, уменьшающих адрес ячейки па­ мяти на единицу. Всему оператору предшествует логический, осуще­ ствляющий проверку условия х^х0. Несоблюдение соотношения приводит к останову машины. Оператор (5-17) требует по одной ячейке на каждую команду. Поэтому общее их число

 

 

m = 3j'+ко+Ал-г-ки+Ам-И,

(5-18)

где

k0 количество

команд для

выполнения логической

операции

для

трех-, двух- и

одноадресных -машин

соответственно

-равно 1,

2, 3; ел — количество

логических

команд; é n — количество

вычисли­

тельных команд; kM—количество

команд

модификации,

включаю­

щее

количество ячеек для хранения числа

модификации

(для ука­

занных машин соответственно равно 6 и 7).

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

. А „ = ^ - у ; ,

- (5-19)

•"-11 + 1 xh

154

необходимо в 'память машины ввести параметры двух узлов для вы­

числения 'Коэффициента интерполяции на

начальном участке опроса

(.ѴІ

и і/і). Тогда память машины,

занятая под знешиюю информа­

цию,

-составит 2(і+1) ячеек. Кроме

того,

потребуется дополнитель­

ная ячейка для хранения промежуточных вычислений. Учитывая, что

ko и Ііц остаются прежними, a к\

и k'M

отличаются

от введенных

в (б-і18), общее число ячеек

 

 

 

m = 2{L+.\)+ko+kn

+

k \ + k'M.

(5-20)

Совместное решение (5-18) и (5-20) относительно числа узлов аппроксимации, при которых предпочтительна программа с вычис­ лением коэффициентов линейной интерполяции, дает:

i = 3 + ( f e ' „ — k B ) + i k ' u ы ).

В полученном выражении kB и А'п , кы и А'„ для трех-, двуходиоадресных машин соответственно равны 5 и 8, 5 и 7, 12 и 16. Ре­ шение показывает, что предельными величинами числа узлов аппро­

ксимации являются

для трехадресной машины 7, для двухадресной

6, для одноадресной

10. Если число узлов превышает указанное, то

следует пользоваться программой с вычислением узлов интерполя­ ции. -Практически количество узлов аппроксимации редко превышает 5—6, тем 'более, что весьма часто встречаются симметричные функ­ ции. Это позволяет вести вычисления ординат только по одной вет­ ви, присваивая знак абсциссы вычисленной ординате (например, при вычислении -выходных координат реверсивных усилителей или элек­ трических машин). Кроме того, дополнительное число ячеек при •росте узлов аппроксимации, как это видно из (5-і18) и (5-20), уве­ личивается на единицу и не дает существенного перерасхода ячеек.

Решение по (5-19) в режиме фиксированной запятой при малой разности абсцисс функции может привести к переполнению ячеек. В этом случае рекомендуется предварительно определить величину максимального коэффициента линейной интерполяции.'По ней прини­ мается порядок числа, на единицу больший. Все ординаты узлов аппроксимированной функции перед записью по адресам у, . ..

..., у—1) уменьшаются на принятую величину (10, 100 и т, д.).

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

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

ны типа «Проминь». Общее количество ячеек,

необходимых для вы-^

числения функций,

 

Л/ф=А+Зі,

(5-21)

где k соответственно равно 12, 14 и 18 для трех-, двух- и одноад­ ресных машин.

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

155

ГЛАВА ШЕСТАЯ

СИНТЕЗ ПАРАМЕТРОВ АВТОМАТИЗИРОВАННОГО ЭЛЕКТРОПРИВОДА НА ЦВМ В ДИНАМИЧЕСКОМ РЕЖИМЕ

6-1. ВХОДНОЙ И ПРОМЕЖУТОЧНЫЙ УСИЛИТЕЛИ

Полученные выше обобщенные и частные алгоритмы синтеза и анализа элементов электропривода, а также способы вычисления функций по аргументу позволяют перейти к составлению программы синтеза системы в це­ лом. Как и при синтезе на АВМ, удобно выбрать некото­ рую обобщенную схему автоматизированного электро­ привода (см. рис. 1-2), что позволяет конкретизировать расчеты. Частные и обобщенные алгоритмы будем в дальнейшем обозначать вычислительными символами с указанием той информации, которая необходима и до­ статочна для построения программы. Такой вычислитель­ ный оператор включает информацию о постоянных рас­ четных величинах, необходимых для вычисления данной переменной, и о переменных, подлежащих предваритель­ ному расчету. Записанный в лаконичной форме вычисли­ тельный оператор позволяет удобно расположить посто­ янные в массивах и дает представление о порядке со­ ставления программ, операторных программ и программ в алгоритмических языках.

Ниже приводится весь комплекс вычислительных опе­ раторов типовых алгоритмов синтеза автоматизирован­ ного электропривода, записанных для решения задач син­ теза относительно параметров входного и промежуточно­ го усилителей. Синтез параметров автоматизированного электропривода предусматривает расчет параметров входного усилителя при наличии пяти каналов управле­ ния: задающего, отрицательной обратной связи по вы­ ходной координате, реализуемой либо в виде связи по скорости, либо связи по напряжению, обратной связи (положительной или отрицательной) по току (моменту) двигателя и автономных отрицательных обратных связей по выходному напряжению усилителя — жесткой и гиб­ кой. Последняя может охватывать два усилителя, что не отражается на структуре расчета. Такая конкретизация позволяет более строго подойти к построению общей про­ граммы синтеза или операторной программы, исключив лишние вычислительные операции. Так как более слож-

156

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

Исходной

при синтезе является

операция

вычисления

ke(ù(tn)- Формирование

текущего значения

времени при

известном шаге синтеза

осуществляется для я-й точки по

соотношению

• •

 

 

 

 

^ = і ' а ^ + ^ = 'П-, + Д^

(6-1)

Присвоим

этой

вычислительной

операции — вычисли­

тельному оператору

символ

 

 

 

 

[t(tn-u М)]п,

 

(6-2)

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

В гл. 2 указывалось, что все динамические характери­ стики могут быть приведены к трем законам управления: оптимальному, рациональному и с ограничением первой и второй производных скорости. Для первого из них вы­ числительный алгоритм

^=*!^(ь-£У (6-3)

где Ае шм — максимальная э. д. с. двигателя на заданной динамической характеристике; Г ц — время цикла.

Вычислительный оператор имеет вид:

І[АеСо(^ым; 7"ц; [tn]; 4)]„.

(6-4)

Для закона рационального управления

£(j©= ±at±âg(ùN(;tY) +bt,

(6-5)

где а и b — соответственно приведенные

ускорение' и за­

медление.

 

Знаки в (6-5) определяются режимом работы: —соот­ ветствует торможению; + соответствует разгону. Разго­ ну соответствуют условия Ав шм = 0 и 6 = 0 при изменении времени в пределах O^t^t-p, где fp =èe<aM /a. Установив­

шееся движение

при

аналоге

максимальной

скорости

/г<*шм имеет место

на

интервале

tp^ti^tp + ty,

где ty

время установившегося движения. Торможение происхо­ дит при условии а = 0 в пределах изменения времени

ір+'іу<й^Тд,

где Тц — tp + ty+tv, U = ke(ihilb — время тор­

можения.

 

157

Символ вычислительного оператора

b)]n- (6-6)

Закон движения с ограничением первой и второй про­ изводных скорости описывается уравнениями:

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

к°ш = ^Ж'

(6-?)

для участка установившегося тока

для третьего участка изменения тока

кеш = ^ L + i ^ L (ta - tJ+^L. { і _ ty,

(6-9)

где ki — темп изменения падения напряжения в якорной цепи, В/с; ti, t% — соответственно начало и конец участка установившегося тока (рис. 6-1).

Рис. 6-1. Динамические харак­ теристики привода при разго­ не с ограничением первой и второй производных скорости.

Соотношения (6-3)—^(6-9) позволяют перейти к любо­ му комбинированному закону управления.

Второй заданной характеристикой является момент статического сопротивления, представленный в виде ана­ лога— падения напряжения в якорной цепи. Ниже при­ водятся вычислительные операторы алгоритмов для раз­ личных функциональных зависимостей статического мо­ мента:

Д Л Я ПОСТОЯННОГО McRo/kM = IclRo

i[Mc(/citfo)]«;

(6-10)

158

для зависящего

от скорости

McRJkM==kew*-\-IC2R0

с(Ып,

а, /с2/?о)]„;

 

(6-М)

для зависящего от -времени McRo/kM=f3(tn)

4-/с 3 і?0

{ЛГсШп,

И». /озАо)]»;

 

(6-12)

для зависящего

от

пути

 

 

 

 

 

f п

 

 

 

McRJkM

= h

 

/ +

/„/?„

 

[ M c ( [ / 4 ] n ,

[ДА,»]*. At,

IMU

 

. (6-13)

где /оі^о постоянная

составляющая момента.

Что касается f3(t)

и f4 ( J ke&dt),

то

при

кажущемся

их многообразии

все

они

приводятся

к ступенчатому

виду либо к синусоидальной зависимости. В некоторых случаях, например в шахтных подъемниках, зависимость может быть иной. Располагая приведенными выше соот­ ношениями, можно вычислить по (5-4) либо по обобщен­ ному алгоритму (5-13) inRu с соответствующим введени­ ем переменных. Символ вычислительного оператора при этом будет иметь вид:

[ « о ( { М « , п - і ; с]п\ В; М)]п.

(6-14)

Вычислительные операции для определения э. д. с. и напряжения преобразователя совпадают. Их удобно сов­ местить. Вычислительные операторы можно записать в виде

[е„([йеш]„; {iRRo]n,n-i\

Т0, At,

AU)]n;

(6-15)

> и ([ke*U

M o k « - , ;

т0;

д £ / д ) ]

, ( б - і б )

где и"я — напряжение на выходе источника, питающего двигатель.

В (6-15) и (6-16) индексы п и (п—1) при iaRo свиде­ тельствуют о необходимости для реализации расчета рас­

полагать

двумя значениями этой переменной — в начале

и конце

интервала.

 

Следующая

в соответствии

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

электропривода

(см. рис. і-2)

операция вычисления мо­

жет выполняться по (5-9) — частному алгоритму либо по

159

обобщенному алгоритму (5-13). Для генератора т = 1 , а Г*д представляет собой отношение постоянной времени к кинетическому передаточному коэффициенту генерато­

ра (см. гл. 2) и является

постоянной

величиной. Вычис­

лительные операторы

 

 

 

 

[Un ([£H]n)]nî

)

(6-17)

[«и ([<?Іі]п, 71-1>

^*ІП

["ll]rt! ДО]

 

Необходимо иметь в виду,

что ип

= иит

поэтому при

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

* Узел модели для синтеза входного напряжения «З .п промежуточного усилителя включает три группы алгорит­ мов: вычисления выходного напряжения гибкой обратной связи, определение функций по аргументу, вычисления ис­ комой переменной. Если принять, что дифференцирую­ щее звено выполнено в виде моста, построенного на об­ мотке возбуждения генератора, то необходимо будет учитывать изменение индуктивности и, следовательно, постоянной времени цепи возбуждения в зависимости от величины и„. Так как Тп(и„) при небольшом изменении входной величины изменяется мало, то вычисления ука­ занной функции можно выполнять по начальному значе­ нию выходной переменной или конечному значению вход­ ной (5-7). Алгоритмы вычислений выходного напряже­

ния дифференцирующего звена записаны для

условий

а 1 с = 0; ß 2 c = l ; £ c n = const. Символы

вычислительных опе­

раторов

 

 

п ([ин]п)]п»

àt)}

\ (6-18)

 

 

Группа алгоритмов для определения функций по аргу­

менту

 

 

["*в.п([ии]„)]«;

[

(6-19)

i " * t . n ( l W , W ] n ;

["*r.n([«"r.n]n)]nv

Алгоритм основной вычислительной операции относи­ тельно «*з.п при m = 4 имеет следующий вычислительный оператор:

160

Соседние файлы в папке книги из ГПНТБ