Баев Теория колебаниы 2015
.pdfВыражение (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Žlω 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. Линейные системы с переменными параметрами
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.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 |
|||
&‡ |
ω dτ |
ω |
&τ |
|
|
|
|
|
|
|
|
Запишем выражение для полной энергии системы
M 1 , &H! H - 2 &τ
2.2.15
2.2.16
2.2.17
130