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

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

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

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

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

H Q=.

2.1.54

Введя эту функцию в уравнение (2.1.52), получим

характеристическое уравнения системы относительно величины r

 

 

 

2λ ω

0,

 

 

 

 

 

2.1.55

решив которое, найдем значения величины r:

 

 

 

2.1.56

 

 

 

 

λ s l

Ž

ω

 

λ .

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение уравнения (2.1.52) имеет вид

 

 

2.1.57

 

 

 

H ž Q= ž Q= .

 

 

 

 

Если диссипативные силы не велики,

 

 

 

 

 

 

 

ω ¤ λ , то общее решение

(2.1.57) преобразуется к виду

 

 

 

 

 

 

 

 

 

H Q+G cos ω α ,

 

 

2.1.58

где

 

 

 

ω1 gω2

λ2.

 

 

 

 

2.1.59

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функция (2.1.58) описывает затухающие колебания с показателем

(или декрементом) затухания

λ

. При этом энергия системы убывает,

как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.1.60

где

o

 

 

o o Q+ G,

 

 

 

 

 

начальное значение энергии.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, колебания

 

В

случае больших

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

сил,

 

 

отсутствуют, движение монотонно затухает и

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

 

ω

 

N λ

(2.1.57), в которой

101

 

 

 

 

 

λ s

 

λ

 

ω

 

.

2.1.61

 

 

,

 

 

Ž

 

 

 

Наконец, если

 

 

 

 

,

то

 

движение тоже не

имеет

колебательного

характера, хотя

и не обязательно монотонного

 

 

ω

 

λ

2.1.62

 

 

H ž

ž Q+G.

У системы со многими степенями свободы

2.1.63

 

 

 

 

α

 

H $ 1, … , % ,

 

 

 

 

 

 

 

 

 

 

 

и уравнения движения, соответствующие функции Лагранжа (2.1.11), принимают вид

 

 

UL α

 

H Q

 

U

0 $ 1, … , % .

2.1.64

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если, как и в одномерном случае, положить

2.1.65

 

 

 

 

H 3 Q=,

то дифференциальные уравнения движения (2.1.64) сводятся3 к линейным алгебраическим уравнениям относительно амплитуд :

3

 

 

 

 

 

α

 

Q

 

0

 

$ 1, … , %

,

2.1.66

 

 

 

 

 

 

 

 

 

 

 

 

 

 

которые имеют нетривиальные решения при условии, что

2.1.67

 

 

 

&Q

α Q 0.

 

Выражение (2.1.67) является характеристическим уравнением системы относительно величины , корни которого определяют

102

собственные частоты колебаний системы и постоянные λ изменения амплитуды колебаний.

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

удобства, представим в комплексном виде

 

 

 

 

 

 

 

 

HL 2λH ω H # Q > .

 

 

2.1.68

 

Как и ранее,

будем

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

 

решение

неоднородного

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

 

U C ,

 

 

 

 

2.1.69

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где амплитуда вынужденных

колебаний

O

в

общем

случае также

комплексная величина

C b .

 

 

 

 

 

2.1.70

 

 

 

 

 

 

 

 

 

 

Подставив функцию (2.1.69) в уравнение (2.1.68), найдем

комплексную амплитуду вынужденных колебаний

 

 

 

 

C

 

ω

2

2

2lλγ .

2.1.71

 

 

 

 

γ

 

 

Ее модуль представляет так называемую амплитудно-частотную

 

системы,

 

 

 

 

связывающую

 

γ

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

 

 

 

 

амплитуду

 

вынужденныхколебаний

системы

с

амплитудой

и частотой

 

внешнего на нее воздействия

 

 

γ

 

 

 

 

 

 

 

 

 

 

b

g ω

2

 

γ

2

,

2.1.72

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

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

103

tgδ |λε|.

tgδ γ

22λγ 2

.

ω

Общее решение уравнения движения имеет вид

H Q+G cos ω α cos γ δ .

2.1.73

2.1.74

Амплитуда вынужденныхγ ω колебаний, хотя и возрастает при сближении частот и , но не обращается в бесконечность, как это было при отсутствии диссипативных сил. С помощью несложных выкладок можно убедиться в том, что амплитуда вынужденных колебаний достигает максимума (наступает резонанс) при частоте внешней силы

γres gω2 λ2.

Вблизи резонанса, где

принимают вид:

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

|ε| |γ ω| q 1,

r 2#ω√ε λ ;

2.1.75

амплитудно- и

2.1.76

2.1.77

Функция δ γ всегда отрицательна, так как представляет «запаздывание» реакции системы из-за ее инерционности на действие вынуждающей силы. Вид этой функции изображен на рис. 2.1.3.

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

104

 

δ

 

 

 

 

 

ω

 

 

 

 

 

 

γ

 

 

0

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

λ

N λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.1.3. Фазочастотная характеристика системы

 

 

Найдем энергию, поглощаемую системой в среднем за единицу

времени. Согласно выражению (1.59)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&o

 

 

 

 

u

 

 

 

 

 

 

2.1.78

 

 

 

 

 

 

& 2 ,

 

 

 

 

 

 

 

 

усредненная по периоду колебаний системы функция Релея.

где

Вuданном случае

αH

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ#

 

γ

 

 

 

γ δ

2.1.79

 

2

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

j γ λ# γ ,

 

 

2.1.80

и

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

так что вблизи резонанса

 

 

 

 

2

 

 

 

 

λ

 

 

 

 

 

&o

 

 

 

 

· ε

 

 

.

2.1.81

системы, где ¬

 

&

j ε 4

 

 

λ

 

/4#λ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На рис.2.1.4

изображена аплитудно-частотная характеристика

 

 

 

 

 

 

 

105

 

 

 

 

 

 

 

 

 

I/I0

 

1

 

 

1/2

 

λ

λ

ε

Рис. 2.1.4. Амплитудно-частотная характеристика системы

2.1.4. Критерии устойчивости систем

Критерий Рауса – Гурвица. Как было установлено выше, чтобы получить нетривиальное решение уравнений движения, должно выполняться условие (2.1.67). Это условие является алгебраическим уравнением относительно величины из представления решений уравнения движения в виде экспонент (2.1.65) и, как уже говорилось выше, называется характеристическим уравнением системы. Его корни определяют собственные частоты колебаний системы и показатели изменения амплитуд этих колебаний.

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

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

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

106

Область Область устойчивости неустойчивости

Рис. 2.1.5. Комплексная плоскость корней α lω

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

связан с числом степеней свободы системы.

 

 

 

 

Левая часть этого уравнения имеет вид полинома

 

. 2.1.82

4

 

 

 

 

 

 

 

+

 

 

 

+

l

+

\

 

 

 

 

 

 

 

 

 

 

 

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

Матрицей Гурвица полинома (2.1.82) называется матрица вида

 

 

2

,

K

L

 

 

 

,

 

, N, … ,0

 

 

 

 

 

 

2

 

K

 

 

 

,

 

,

 

, , … ,0

²

 

¯0 ,

,

, , … , 0

2.1.83

0 ,

,

, , … , 0

 

… … … … … … … … .

°

 

-0, 0 ,0 , … … … …

 

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

107

Главными диагональными минорами матрицы Гурвица назовем определители вида

 

 

\

;

 

 

 

 

 

 

2.1.84

 

\

 

 

 

 

 

,

 

 

;

 

2.1.85

 

 

³

 

,

³

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

3

 

 

0

,

2

,

4

m ;

 

 

 

 

 

 

 

, 1

 

, 3

2.1.86

 

\ m

0 ,

3

,

5

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

 

 

0

,

2

, 4

,

6, … ,0

 

$

 

 

,

 

 

,

 

,

, … ,0

 

n0 , 0

, 2

, 4

 

, … , 0n

 

\

 

1

 

 

3

 

 

5

 

 

7

 

2.1.87

0 ,

1

,

3

,

5

 

, … , 0 .

 

n … … … … … … … …$. n

 

 

0, 0 ,0 , … … … …

 

Критерий Рауса – Гурвица. Система с вещественными коэффициентами в характеристическом уравнении и положительным старшим коэффициентом устойчива тогда и только тогда, когда все главные диагональные миноры матрицы Гурвица положительны.

Утверждение критерия Рауса – Гурвица приведем без доказательства, которое, как и более подробную информацию о критерии, можно найти в книге [11].

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

2 1 0.

2.1.88

108

Из требования, чтобы главные диагональные миноры матрица Гурвица

, 1, 0

@1, , 0A 2.1.89 0, , 1

были положительными, легко получаем условие устойчивого движения системы:

¤ 1.

2.1.90

Важное достоинство критерия Рауса – Гурвица состоит в том, что он связывает параметры системы, через которые выражаются коэффициенты характеристического уравнения, условиями устойчивости, записанными в аналитическом виде, как в приведенном выше простом примере (2.1.90). Это самый сильный результат, который может быть получен при анализе математических моделей, поскольку из аналитических выражений отчетливо видно, как изменение тех или иных параметров системы в тех или иных пределах влияет на устойчивость работы системы.

Критерий Михайлова. Этот критерий устойчивости непосредственно вытекает из известной формулы теории функций комплексного переменного:

 

 

1

 

4! ]

&] †,

2.1.91

 

4 ]

2πl

(C)

4 ]

где

функция, аналитическая всюду в области,

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

 

замкнутым контуром (С), за исключением тех точек, в которых она

имеет полюсы. На самом контуре (С) функции

 

достаточно быть

 

 

 

 

4 ]

нуль.

и

непрерывной, и она не должна на нем

обращаться в

 

4 ]

 

 

 

области,

 

соответственно число нулей и полюсов функции

 

в

 

 

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

109

Если функцию 4 ] записать в показательной форме

2.1.92

 

 

 

 

 

 

4 ] 3 T,

 

 

 

 

 

 

 

 

то формула (2.1.91) приводится к виду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

&ψ †.

 

 

 

 

 

 

2.1.93

 

 

 

 

 

 

(C)

 

 

 

 

 

 

Формула (2.1.93) представляет так называемый принцип

аргумента. Согласно этому принципу, когда точка

на комплексной

плоскости

 

 

 

обходит контур (С),

 

вектор,

]

изображающий

однозначную ]

аналитическую функцию

 

 

 

 

на

комплексной

плоскости (

 

 

вращаясь в положительном

направлении вокруг начала

 

 

 

4 ]

 

 

 

координат,

делает число оборотов, равное разности между числом

4 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нулей и числом

полюсов

функции

 

 

,

 

заключенных

внутри

замкнутого

контура (С), в

точках

которого

функция не

должна

 

4 ]

 

 

 

 

 

 

 

4 ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

функции

 

 

считается столько раз, какова его кратность. Сказанное

иллюстрирует рис.2.1.6.

 

 

 

 

 

 

 

 

 

 

 

 

 

jy

 

 

(z)

 

 

 

 

 

jv

 

(F)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U °

 

 

 

 

 

 

 

 

F (z)

 

 

 

 

 

 

° U

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

ψ

 

 

 

 

 

U°

 

 

 

 

 

 

 

 

 

 

 

 

 

 

°

x

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

(С)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

 

 

]( H lI

 

б

 

 

 

 

 

 

 

 

функции

4 ]

cнулями и полюсами

Рис. 2.1.6. Плоскость корней

 

 

 

 

 

 

 

 

 

 

4

] ¢ lP

 

а);

 

 

 

плоскость

вектора

 

 

 

 

 

 

 

 

 

 

 

 

(б).

110