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

Баев Теория колебаниы 2015

.pdf
Скачиваний:
22
Добавлен:
12.11.2022
Размер:
2.85 Mб
Скачать

Выражение (2.1.117) представляет характеристическое уравнение системы с обратной: связью, . Если в него ввести выражение (2.1.116) для функции то получится рациональная функция с одинаковыми степенями полиномов в числителе и знаменателе:

v

v

+

v

+

2.1.118

4

 

+

 

+

.

 

 

 

 

 

 

 

 

Об устойчивости системы с обратной связью будем судить по

значению интеграла (2.1.93), взятого вдоль контура, изображенного на

начала

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

рис.2.1.7. Система с обратной связью будет устойчива, если при

обходе точкой

 

этого контура число оборотов вектора

 

вокруг

 

 

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

оборотов вектора

 

 

вокруг точки (1,0) равно нулю.

 

 

 

 

 

 

 

 

 

 

 

 

полуокружности

 

сколь угодно большого радиуса

Поскольку на :

 

 

 

 

 

 

 

контура интегрирования функция

 

 

, как это видно из выражения

(2.1.118), является константой, то

упомянутое выше число оборотов

4

 

 

 

 

 

 

 

 

следует

считать

 

при

перемещении

 

точки

по

 

мнимой оси

комплексной плоскости (

 

. Собственно в этом и состоит критерий

Найквиста

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Критерий Найквиста. Система с передаточной функцией

 

 

,

устойчивая при отсутствии обратной связи, сохраняет

устойчивость и

 

 

:

 

при

 

 

: lω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

при наличии обратной связи тогда и только тогда, когда годограф

Кстати,

 

ω

 

 

 

 

 

 

 

 

 

 

не обходит точку (1,0)

вектора

 

 

на комплексной плоскости (

изменяется

в пределах от

 

 

 

до

 

 

, есть не что иное,: lω

 

ω

 

изменении

 

от

 

до

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

линия, которую описывает конец вектора

 

, когда

 

 

 

 

 

 

 

 

 

 

 

 

 

плоскости

(

 

 

 

 

как граница Д-

разбиения

 

 

 

 

 

 

 

 

 

. Эту границу в критерии

комплексной

 

 

 

 

 

 

 

 

 

 

 

Найквиста называют кривой или

диаграммой Найквиста.

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

В качестве примера на критерий Найквиста решим следующую задачу: построить кривую Найквиста и найти условие устойчивости системы регулирования температуры печи, изображенной на рис. 2.1.14; считать, что печь, длиной l, одномерна и распространение в ней тепла от нагревателя описывается уравнением теплопроводности.

Регулирование осуществляется по сигналу ошибок

ε φ φ , где

φ

задаваемый уровень температуры.

 

121

 

Решение задачи. Передаточная функция системы с разомкнутой обратной связью имеет, очевидно, следующий вид – это произведение передаточной функции усилителя: с коэффициентом¸, усиления К, равной: этому коэффициенту, и передаточной функции печи

 

 

 

 

 

 

:

 

 

 

:

 

 

 

:

 

 

.

 

 

2.1.119

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

φ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К

φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Печь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε φ φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Управляющее звено

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Питание

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Функция

:

 

 

определяется

 

из

решения уравнения

 

 

 

 

 

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

 

 

 

 

 

¢

 

 

¢

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

2.1.120

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решим это уравнение методом Фурье, представив искомую

функцию u(x,t) в виде произведения двух функций:

2.1.121

 

 

 

 

 

 

 

¢ H, j H ( ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

122

 

 

 

 

 

 

 

 

 

 

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

1

&(

 

;

 

 

2.1.122

(

 

&

 

 

 

1

 

 

 

 

 

 

 

 

& j

.

 

 

2.1.123

j

&

 

 

 

 

Из решения второго уравнения

B H

 

 

 

 

 

вых

ch √

!

2.1.124

j H j

 

 

B

получаем искомую передаточную функцию печи

 

:

2

 

1

;

 

2.1.125

 

ch

 

 

 

 

H

 

 

 

 

в этом решении координата

отсчитывается от конца печи, справа

налево.

 

 

 

 

 

 

 

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

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

 

 

 

:

¸

 

 

2.1.126

ch,

 

а кривая Найквиста описывается функцией

 

 

 

 

 

 

1

 

 

 

 

где

 

 

chλ,

 

 

2.1.127

 

 

 

 

 

 

1

 

 

2.1.128

 

 

 

λ ¸.

 

 

123

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

Очевидно, что чем выше коэффициент К усиления усилителя, тем выше чувствительность узла системы, контролирующего температуру в печи и, следовательно, тем точнее можно поддерживать эту температуру. И здесь естественно возникает вопрос о том, есть ли ограничение на величину этого параметра системы.

jv

(λ) ^?

1 u chπ

Рис. 2.1.15. Д-разбиение системы регулирования работы печи

Ответ на этот вопрос дает найденное Д-разбиение системы: если

система может устойчиво работать, т.е. если

^? 0, то она будет

устойчиво работать при условии, что

1

или, что то же самое, если

λ ¤ chπ

2.1.129

¸ N chπ.

2.1.130

 

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

124

ограничение, накладываемое на коэффициент усиления К, в виде неравенства (2.1.130).

Эта задача представляет интерес еще и в связи с тем, что в отличие от того, о чем говорилось выше, в данном случае передаточная функция системы не является рациональной. Однако критерий Михайлова^ 0позволяет справиться с этой проблемой и убедиться в том,

что ? . Действительно, поскольку

где

то

где

chŽbθcosθ lshθsinθ, θ Œω2 ,

1

 

 

 

 

 

λ ω chθcosθ lshθsinθ ¢ ω lP ω ,

chθcosθ

;

¢ ω chθcosθ

 

shθsinθ

 

 

 

 

shθsinθ

 

.

P ω chθcosθ

 

shθsinθ

 

 

 

 

 

2.1.131

2.1.132

2.1.1332.1.1342.1.135

Как Pвидноω из выражений (2.1.134) и (2.1.135), корни функций u(ω) и строго чередуются на всей области^ Д-разбиения с минимально возможным значением параметра ? , и, согласно следствию из критерия Михайлова, а^ именно: критерию чередующихся корней, система в области с ? устойчива. Задача полностью решена.

125

2.2.4

Глава 2.2. Линейные системы с переменными параметрами

2.2.1. Частные случаи

Перейдем к рассмотрению линейных систем с переменными параметрами. Математической моделью таких систем являются

линейные

 

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

уравнения

 

с

 

переменными

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В качестве примера рассмотрим мятник с колеблющейся точкой

подвеса, колебания которой описываются функцией acosγt.

 

Функции Лагранжа такого маятника

 

 

 

 

 

 

 

 

 

9

#B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

φ #B γ

 

cosγ cosφ #DBcosφ 0

2.2.1

соответствует уравнение движения вида

 

 

 

 

 

 

 

 

 

 

 

D

 

γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

φL , B

 

 

B

 

cosγ - sinφ 0;

 

 

 

2.2.2

здесь

 

длина подвеса маятника;

#

его масса;

φ

угол отклонения

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

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В случае малых колебаний

 

 

 

 

 

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

линейному

дифференциальному

уравнению

 

с

переменным

 

 

 

 

 

 

|φ| q 1

 

 

 

 

 

 

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

φL

 

φ 0,

 

 

 

 

 

2.2.3

 

 

 

 

 

 

 

 

 

 

 

где

DB γB cosγ .

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

126

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

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

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

Например, применительно к уравнению (2.2.3) можно было бы

строить

¤ 0,

 

анализ его решений, исходя из следующих соображений.

Если

 

то, по аналогии с колебаниями гармонического

ω Ž

 

осциллятора, можно ввести понятие мгновенной частоты колебаний

 

и предположить, что колебания маятника, определяемые

уравнением (2.2.3), представляют в каждый момент времени

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

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

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

примера.

 

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

2.2.5

HL 1 H 0.

 

 

127

 

2.2.9

Это один из тех редких случаев, когда уравнение решается

аналитически. Если искать его решение в виде

2.2.6

H G

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

λ λ 1 0.

2.2.7

Решив это уравнение, найдем значения параметра £:

 

1

2.2.8

λ 2 1 s l√3 .

Следовательно, общее решение уравнения (2.2.5) имеет вид

H √ ,ž cos 23 ln ž sin 23 ln -

и оно изображено на рис. 2.2.1.

x

t

Рис. 2.2.1.К решению уравнения (2.2.5)

128

Из рисунка видно, что хотяH сила,0 действующая на систему, направлена к точке ее равновесия и к тому же она со временем быстро ослабевает, движение системы, тем не менее, оказывается неустойчивым. Система совершает колебания около точки

равновесия,

амплитуда

которых

неограниченно

возрастает

пропорционально

 

.

 

 

 

 

Если ввести

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

 

 

 

ω

3ln

,

2.2.10

 

 

 

 

 

2

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

может увести от истины рассуждение «по аналогии с системами,

параметры которых постоянны».

 

 

 

 

Введем

в

уравнение

 

(2.2.5)

диссипативные

силы,

пропорциональные величине 1/t:

1

 

 

 

 

 

1

 

H 0.

 

2.2.11

 

 

HL

H

 

 

 

 

 

 

 

 

 

Это уравнение также решается аналитически. Если его решение искать в виде (2.2.6), то после подстановки этой функции в уравнение

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

уравнению

относительно параметра λ: λ 1 0,

2.2.12

так что общее решение уравнения (2.2.11) имеет вид

2.2.13

H ž cos ln ž sin ln .

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

129

2.2.2. О некоторых достаточных условиях устойчивости движения

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

В первом примере рассмотрим уравнение, описывающее колебания широкого круга одномерных систем и поэтому достаточно

общее:

HL D H ω H 0.

2.2.14

 

Сделаем в этом уравнении следующую замену независимого переменного

 

 

‡ ; ω a &a

,

 

 

 

 

 

 

 

 

 

 

 

 

после чего оно примет вид

 

! &H

H 0.

&

H

1 dω D

&‡

ω

ω

 

 

 

 

 

 

 

Запишем выражение для полной энергии системы

M 1 , &H! H - 2 &τ

2.2.15

2.2.16

2.2.17

130