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

книги / Теория автоматического управления техническими системами

..pdf
Скачиваний:
15
Добавлен:
12.11.2023
Размер:
15.06 Mб
Скачать

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

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

3. Вводят обозначения относительных единиц и коэффи­ циентов уравнения.

Вкачестве примера рассмотрим методику составления уравнений сис­ темырегулирования угловой скорости вала электродвигателя. Для объекта регулирования исходным является уравнение баланса моментов:

(3.2)

где У—момент инерции движущихсячастей, приведенный к валу электродви­ гателя; (о—угловая скорость вращения вала; Мя—движущий (вращающий) момент на валу электродвигателя; Мс —момент сопротивления на валу; /—время.

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

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

Мя=-Мд(ю, /).

Момент сопротивления зависит от угловой скорости вала двигателя о), угла установки лопасти винта <рл и ряда других факторов (плотности воз­ духа, скорости полета и др.), изменение которых учесть трудно. Выражение для момента сопротивления имеет вид

Мс=М(со, фл, /).

Основываясь на теории двигателей, можно записать аналитические зависи­

мости полученных функций или представить их в

виде графиков ;,ит. д.

Далее, если приведенный к валу электродвигателя

момент инерции вращаю­

щихся частей

является постоянным, то, как следует из (3.2), уравнение уста­

новившегося

режима, или уравнение статики при ю=const имеет вид

МЛо=Мсо.

(3.3)

Малые отклонения моментов (приращения) от установившихся значений обо­ значимчерез ДМд и АМС.Тогда

Мд=Мдо+Д/Ид; Мс==Мсо-|-|ДМс.

Теперь, подставив в исходное уравнение режима (3.3), получим

d(а

/~^ =АМЛ—АМС.

(3.2) выражение установившегося

(3.4)

Для линеаризации выражения (3.4) следует воспользоваться формулой Тей­ лора (3.1). Приращения |ДМд и ДМСопределим из выражений:

ДЛ4д=-^Г Дш+ ДМд (0.

 

 

 

Дш+

Дфл+ ДЛ*сW

 

где АМд(0

и ДМс(0 —составляющие приращений Мл и Мс, изменяющиеся

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

 

Угол установки лопастей винта <рл изменяют с помощьюисполнительного

механизм

регулятора

(серводвигателя), координату

перемещения которого

обозначают через т:

 

 

 

 

Функция <рn=f(m) задается

обычно

графически, и частную производную

(дул)1дт)

определяют

как

тангенс

угла наклона

касательной к кривой

фл—/(ш)

в точке, соотве1СТвующей

установившемуся режиму. Тогда

Полученны выражения подставляю в уравнение (3.4):

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

Для приведения дифференциального уравнения в отклонениях, выражен­ ного в абсолютных единицах, к уравнениюв относительных единицах и с безразмерными коэффициентами выбирают номинальное значение момента Мни на него делят каждый член уравнения (3.5):

где ДМ(0=АМд(0-А^с(0.

Врезультате все членыуравнения стали безразмерными. Для координат каждого приращения, входящего в полученное уравнение, выбирают соот­ ветственно некоторы постоянные значения. Так, для угловой скорости — номинальное значение ©н, для координаты серводвигателя —максимальное значение тт. Каждый член уравения, в который входит та или иная пе­ ременная, делят (умножают) на соответствующее ей выбранное постоянное значение.

После этого уравнение (3.5) будет иметь вид

Уа>н d(ù

а>„

/дМс дМД\ Аса

АМ(t)

тт дМсдул Ат

AlHtoHdt + Мн\ d©~ d(ù J

(ùu ~

Мя

~~ Мн<?Фл àm тт‘

Учитывая, что

 

 

 

 

 

d(ù

d (А©)

;

дМ дАМ

и т*д*

 

dt

— dt

d(ù~ д (А©)

 

последнее уравнение переписывают в виде

АЛЛ(О

АЛЛС

(?Фл

Ат

 

д М„

(3.6)

ЛЛ„

<?фл

^ Aw

тт'

Вводят обозначения:

Уравнение (3.6) двигателя как объекта регулирования принимает вид

нли окончательно

 

TvlU~+kc(V=fМ(0—

(3.7)

Все величины, входящие в уравнение (3.7), за исключениемt ипостоянной Тчприведенык безразмерному виду. Полученное уравнение справедливо идля других двигателей, так как их динамика во многоманалогична. Часто используют другуюформу этого уравнения (для инерционного объекта ре­ гулирования) :

T^+x2=f(t) + bxlt

где *2=Ф; *i—И; Г —постоянная

времени объекта регулирования (Т=

—Tç/Ke)’, f(t)= (1/£е)М0);

---передаточный коэффициент объекта.

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

ной форме: {Тр-fi)*2=f(0—

{T22P*+t i P ) —Xi•

T22p2-\-TlP+\)X2=f(t)—kXU

где p —символ дифференцирования (p^df(dt)).

Методика составления уравнений элементов автоматическо­

го регулятора аналогична рассмотренной ранее

методике

со-

 

ального закона изменения величи-

 

ставления уравнений объектов.

Во

 

всех

уравнениях

Т представляет

 

собой

постоянную времени

регуля­

 

тора.

случае

экспоненциального

 

В

 

(инерционного)

процесса

постоян­

 

ная времени Т может быть опреде­

Рис. 3.4. Определение постоян­

лена

известным

способом, показан­

ным на рис. 3.4.

 

физики

процес­

ной времени инерционного про­

С

точки зрения

цесса

са постоянная времени экспоненци-

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

постоянной скоростью, равной скорости ее изменения в началь­

ный момент, т. е. если бы кривая изменения х, начиная

с это­

го момента,

совпадала с

касательной

(к экспоненте),

прове­

денной в ее начальной точке.

 

 

 

 

3.3.

Свободные и вынужденные колебания

САР.

 

 

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

 

 

Система

линеаризованных

дифференциальных

уравнений,

описывающих

физические

процессы

элементов САР,

может

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

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

 

dnx

.

dn-lx .

I

_ dx I

 

 

 

an~ T t+

^ягг+ • • • +01^7+ 00*—

 

 

==

 

bm\ d^mlx+ ... + bx^7+ b0f.

 

(3.8)

Здесь f(t)—входное воздействие; x(t)—изменение выходной величины; о/, bi—постоянные величины, которые определяются физическими (техническими) свойствами системы.

Уравнение (3.8) описывает физические процессы, т. е. изме­ нение переменных х и f, в замкнутой системе регулирования. Его решение может быть:

т. е.1) общим, когда правая часть уравнения (3.8) равна нулю,

« . & + * « £ 5 - + . . . + * £ - м - (3-9)

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

Общее решение уравнения (3.9) для

случая ап=ап-\= ...

... —а^—Онаходим в виде

 

 

 

 

(3.10)

*(0=еЧ

(3.9)

позволяет

 

Решение такого вида уравнения

определять

свободное движение системы, т. е. колебания, которые

совер­

шает система, выведенная из состояния равновесия некоторым

воздействием после того, как это воздействие

стало

равным

нулю.Подставляя выражение (3.10) в (3.9), получим

 

(аДп+а„_1Г-,+1... +аД+ао)е*‘«0.

 

(3.9) при усло­

Выражение (3.10) является решением уравнения

вии, что Я—корень характеристического уравнения:

 

 

ûnA»n-f-ûtn-|Xn-1-f-... -|-fl|A.-|-ûo==0.

 

 

корней (Я|,

Так как последнее уравнение имеет п различных

Яг,..., Яп), то общее решение

уравнения (3.9)

 

может быть

представлено в виде

 

 

 

 

 

х(t)= С,еЛ*'+ Сзе**'+ ... +

,

 

 

 

 

где Ci, С2,..., Сп —постоянные, зависящие от начальных усло­ вий.

Решение уравнения (3.9) удовлетворяет условию Птx(t)=0,

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

Вычисление вынужденных колебаний. Пусть воздействие f(t) представляет собой гармоническую функцию времени, т. е. линейную комбинацию sin at и cosю/:

f(0=foCOs((ûH-<p),

где ©—угловая частота воздействия; /0—амплитуда; <р—

фаза.

Если начальная фаза фо^О, то f(t)=focosиt.

Гармоническую функцию можно рассматривать и как сумму двух экспонент:

/ (*)= Лcos = ТеУа<+4°е-У“'’ J= vr=T-

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

При f(f)=YeJ<ùt частнс^ решение уравнения (3.8) имеет вид

[о.(/о)"+а._,(/ш)"-Ч-... +0, (ja)+aa]Y(j<i,)e>°’=

 

=[ё„(/ш)"+6„_, (ja)--‘+

(/ffl)—ь0]е»‘.

 

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

 

 

Г(уш)=

bm(JOi)m+b,

 

 

чя------W ”-1+••• + *»(/©)+ *.

(ЗЛ1)

 

ал (У©)в+дл-1

+ ... +а,(/со)+а0'

Последнее выражение может быть представлено через ампли­

тудную Л(а>) и фазовую <р(ш)

частотные

характеристики си­

стемы:

 

 

 

 

У(/©)=Л(©)е™“>.

Y{j(ù)

в формулу

для *ib(0 запишем

Подставляя

выражение

Х\ь(/)

(©)

 

 

 

Если в уравнении (3.8) / (*) =

, tq

вынужденные коле­

бания

 

 

 

 

*2в(0=4°Г(-./ш) е-^,

 

 

или

 

 

 

 

х2в {*)= Т-£А (ш)

ч»с»)1.

 

 

Когда f(t)—гармоническое воздействие, т. е. ffjt)=focos©f, то

хв.(f). =xlB(t)+х2в (f)=-2 А (©)[е*“*е™»>+

 

(3.12)

+е_3“'е-*ф(в1)]=/оЛ (©)cos(©Жр (©)).

 

 

вынужденные

колебания,

Выражение

(3.12) показывает, что

вызываемые в линейной динамической системе гармоническим

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

цию времени, имеющую ту же круговую частоту ©, что и воз­

действие, но отличающуюся по амплитуде и по фазе. При этом

относительная амплитуда и фаза

вынужденных

колебаний

определяются амплитудной Л (о) и фазовой ф(©) частотными

характеристиками.

 

(3.11) в числителе и

Если в выражении для функции У(/©)

знаменателе отделить вещественную

часть

от мнимой, то

оно

может быть записано

 

 

(3.13)

r(ja)=Jtpl±^pL,

 

 

w '

c(üf)+Jd((ù) ’

 

 

v

'

a(©)=&o—&2©2+64©4—. . . i b (cd) =&i(d—bz(ùz+bs<ù6—... ; C(û)) =ûo—Û2Û)2+a4©4—... î d((ù)=ai(û—аа(о3+а5(й5—••• •

Выражение

(3.13) может быть представлено в виде

У(;m4^

g (<■>) с (<■>) + ь (ф) d (со) ,

6((0)с(<а)—g(<р) Çca)

' '

с2(а)-М2(со)

с2(©)+d2(<■>)

или

y(/0)=P(û>)+/Q(0),

где Р(©), Q((ù)—вещественная и мнимая частотные харак­ теристики системы соответственно:

_ P(û>) = (ac+bd)/(c2+d2); Q(©)= (bc—ad)/(c2+d2).

Амплитудная и фазовая частотные характеристики САР опре­ деляются выражениями:

А(©) = V~(a2+b*)/(<?+d*), ф(©)=arctg[(&c—ad)}{ac+bd)].

Частотные характеристики А(ш) и ф(щ) связаны с характери­ стиками Р(о) и Q(a>) следующими соотношениями:

Р(ш) —А (q) • cosф(©);

<2((о)=А (и) - sin ф(©); Л(ш)=/Р» + (3»;

<p(o))=arctg[SM],

На основе выражения (3.12)- можно экспериментально определить ам­

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

го к систем

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

воздействие xOx(0 =cosш/,

И=1), имеющее угловуючастоту а. Врезультате в систем возникают

переходный процесс и вынужденны колебания с частотой ©.Через неко­

торое время

произойдет затухание переходного

процесса, если система

Рис. 3.5. Экспериментальное определение частотных характеристик динамической си­ стемы(динамического звена):

а —система илизвено; б—процессына входе и выходе

устойчива, иостанутся лишь вынужденны колебания. Они будут иметь час­ тоту (о, равнуючастоте воздействия, но отличаться от входного воздейст­ вия по амплитуде ифазе (рис. 3.5). Амплитуда выходного сигнала и угол сдвига фазывыходного сигнала xDUï=A(co)-cos[(i)f+q)(co)l по отношению

к входному зависят от угловой частотыш.

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

Для экспериментального определения частотных характеристик исполь­ зуют специальнуюнизкочастотнуюаппаратуру, так как САР—обычно низ­ кочастотные системы. Всостав аппаратуры входят, например, следующие приборы(см.рис. 3.6) :

Рис. 3.6. Схема включения измерительной аппаратурыдля оп­ ределения частотных характеристик

низкочастотный генератор периодических колебаний (ИГПК) для гене­

рирования как входны колебаний синусоидальной, прямоугольной, треуголь­

ной, трапецеидальных форм, так и одиночных импульсов тех же форм;

низкочастотный фазометр-частотомер (НФ) для определения частотыи

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

ной частоты (100 кГц)

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

и за время между двумя смежными прохождениями через нуль кривых вход­

ного ивыходного напряжений при измерении фазы;

двойной пиковый вольтметр (ДПВ) для измерения амплитудына входе

ивыходе системы;

 

П1 иП2—преобразователи сигнала.

Приборырассчитанына напряжения ±100 Видиапазон частот от 0,001 до

100 Гц.

 

 

Кривые отношения амплитудывыходной переменной хВых к амплитуде

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

ставляю собой амплитуднуюи фазовуючастотные характеристики САР

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

 

Взависимости от условий эксплуатации системыиногда удобно подавать

на ее вход

колебания

прямоугольной формы (включение —выключение) и

треугольной

(равномерное открытие изакрытие регулирующего устройства)

и т. п.

 

 

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

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

Пусть воздействие f(t) удовлетворяет условиям: f(t)=0; ^<0;

J |/(0|е-е*Л<оо,

где с—абсцисса абсолютной сходимости. Тогда преобразование Лапласа для функции f(t)

Если все

члены уравнения (3.8) при нулевых

начальных

условиях умножить на е~*‘ и проинтегрировать в пределах от 0

до оо, то получим

 

(a„s,'+an_is7,_,+ ... +aiS+a0)X(s)=

(3.14)

где= (ôMs"4-6m,s—Ч-... +bls+b0)F{s)t

 

*($)= J x{t)erstdt.

 

о

 

 

Следовательно, X(s)=Y(s)F(s), где оператор

 

Y(s\_ x<s) __ àmsm+b^s"1-1+... +M+Qp

(3.15)

F (s)

ansn +an-lsn-1+ ...+a,s+ fl0

является передаточной функцией динамической системы. Согласно выражению (3.15), передаточной функцией непре­

рывной линейной стационарной динамической системы назы­ вают отношение преобразования Лапласа X(s) переменной x{t)' на выходе системы к преобразованию Лапласа F(s) воздейст­ вия f(t) на ее входе при нулевых начальных условиях.

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

Передаточные функции устойчивых динамических систем обладают следующими основными свойствами:

1) передаточная функция Y(s) представляет собой дробно­ рациональную функцию вида (3.15), причем в реальной систе­ ме порядок т числителя не превышает порядок п знаменателя; 2) все коэффициенты b0, bи ..., bm-1, bm, а0, ai,..., а„_ь а„ передаточной функции вещественны. Это следует из того, что они представляют собой функции параметров системы, т. е. мо­

гут быть только вещественными; 3) невещественные нули и полюсы передаточной функции

могут быть лишь комплексно-сопряженными;

4) все полюсы передаточной функции У(s) расположены в левой полуплоскости комплексной плоскости, что являетсяусло­ вием устойчивости системы;

5) при /ш передаточная функция (3.15) преобразуется в амплитудно-фазовую частотную характеристику системы, а при s=0—в передаточный коэффициент (для позиционных звень­

ев).

 

 

Vpfs)

r(t>

 

V

 

 

 

 

Л

 

%fs)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.7. Схема САР:

регули­

 

 

Wо(s) —передаточная функция

объекта

 

 

рования; Wр (в) —передаточная

функция

регу­

 

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

 

лятора

Общие

дифференциальные

САР.

уравнения

рассматриваемой

 

линейной САР,

приведенной на

рис. 3.7, можно записать в виде:

 

 

 

 

 

D(р)х(t)=М(р)f(t)+С(р)г(0 —для объекта регулирования,

 

 

 

 

 

 

 

 

 

(3.16)

B(p)r(t)=N{p)z(t) —для регулятора,

 

 

(3.17)

е (0=g(t)~х{t) —уравнение ошибки.

 

 

(3.18)

Применяя к уравнениям

(3.16)—(3.18) преобразование Лап­

ласа, обозначим

 

 

 

 

 

 

 

 

*(s)

x(t)e~stdt;

 

 

 

 

er3ldt;

 

 

 

G(s)=^g(t)e~s‘dt;

(s)= jj г(0e-‘W;

 

 

 

 

о

 

 

 

 

 

 

 

 

F(s)= b

 

 

 

 

 

 

 

(3.16)-

Учитывая

начальные условия MEi и Мн2, уравнения

(3.18) можно представить в виде

 

 

 

 

 

D(s X (s)=M(s)F(s)+ C(s)R{s)+Mnl(s);

(3.19)

В (s)R (s)= N(s)E(s)+ Mh2(s);

 

 

 

 

E(s)=G(s)-X(s).

(3.19) исключить функции R(s)

и E(s)

Если из уравнений

н выразить эти уравнения относительно X(s), то