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

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

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

 

L

e(t) ~

C

R

Рис. 1.11. Последовательный колебательный контур

Величина ЭДС равна сумме напряжений на пассивных элементах схемы:

9

&$

1

; $& Q .

1.257

& $ ž

Если перейти от тока к заряду

 

 

 

$

&D

1.258

 

 

& ,

то выражение (1.257) примет вид дифференциального уравнения второго порядка

9

& D

&D

1

1.259

&

&

ž D Q .

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

магнитным полем в индуктивности

( 1 9D , 1.260 2

потенциальной энергией считать энергию, запасаемую электрическим полем в емкости:

6

1

D ,

1.261

 

81

 

функцию Релея Rel связать с резистором выражением

 

1

 

1.262

Rel 2 D ,

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

 

¡

Rel

 

D ,

1.263

5 Q

то уравнение Лагранжа сведется к уравнению (1.259).

Впредложенной электромеханической аналогии индуктивность является аналогом инерционного коэффициента (массы), величина 1/С соответствует коэффициенту упругой силы гармонического осциллятора, сопротивление R – аналог диссипативного коэффициента. Диссипативная функция Релея – это, очевидно, мощность, рассеиваемая на резисторе, порождающем аналог сил трения. В данном случае аналог обобщенной скорости– ток.

Врассмотренной схеме не учтен еще один активный элемент, а именно: источник тока i(t). Чтобы иметь возможность строить математические модели с источниками тока, рассмотрим параллельный контур, изображенный на рис. 1.12.

i(t)

~ R C L

Рис. 1.12. Параллельный контур

В контуре на рис. 1.12 складываются токи, протекающие через индуктивность, емкость и сопротивление. Если u – приложенное к

контуру напряжение, то

1

1

 

 

 

; ¢& $

1.264

ž &

¢ 9

 

 

 

 

82

 

83
# && && ~ 5¡;

или

&

 

¢

1 &¢

1

 

&$

 

 

 

¢

1.265

 

&

 

&

9

& .

 

 

 

 

 

 

 

 

 

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

(1 ž¢ , 1.266 2

аза потенциальную энергию – энергию, запасаемую магнитнымэлектрическим полем в емкости,

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

6 291 ¢ ;

функция Релея связана с резистором выражением

Rel 1 ¢ , 2

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

&&$ Rel¢ .

1.267

1.268

1.269

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

1.270

здесь £– обобщенная координата; # –~инерционный коэффициент (масса); – коэффициент силы трения; –коэффициент собственной упругой силы системы, обеспечивающий устойчивость ее движения.

Таблица 1.1

Механические

q

m

λ

k

T

 

U

 

Rel

 

5

системы

 

 

 

 

 

 

 

 

 

 

аналогия

D

 

 

ž

2 9D

 

D

 

2 D

 

1-я электро-

 

L

R

1

1

2

1

 

1

 

e(t)

механическая

 

аналогия

 

 

 

9

2 ž¢

29 ¢

 

2 ¢

 

&

2-я электро-

u

C

1

1

1

 

1

 

1

 

&$

механическая

Приведем два примера использования электромеханических аналогий.

В первом примере построим математическую модель

электрической схемы, изображенной на рис. 1.13.

 

 

 

R1

 

R2

 

 

 

e(t) ~

$ C12

 

$ L23

$2

C3

 

 

L1

 

 

R3

 

 

 

 

 

 

 

 

 

Рис.1.13. Электрическая схема

 

 

Решение

задачи.

С

помощью

первой

системы

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

кинетическую энергию

1

 

1

9

2

 

2

 

;

1.271

( 2

9 D 2

 

D D

 

 

 

 

84

 

 

 

 

 

 

 

потенциальную энергию

 

 

 

D

 

 

D

 

D ;

1.272

 

 

 

6

 

 

 

функцию Релея

 

 

 

 

12

2

 

1

 

 

 

 

 

 

 

 

 

 

D

2

D 2

D .

1.273

 

Rel 2

 

 

 

1

 

 

 

1

 

 

1

 

2 2

 

 

От этих

функций,

 

по

 

изложенной

выше

методике,

несложно

перейти к математической

модели

схемы,

 

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

1.13:

9 DL D

1

D

1

 

D

Q ;

1.274

 

 

 

 

 

 

 

 

ž

 

 

ž

 

 

 

 

 

 

 

 

9

DL 9

 

 

 

 

 

 

 

 

D

 

 

 

1

 

D

0;

1.275

DL D 1

 

 

 

 

2

 

 

2 2

 

 

 

 

 

ž

 

 

 

ž

 

 

 

 

 

9

DL 9

 

 

 

 

 

 

 

 

 

 

 

 

 

0;

1.276

 

 

DL D

 

D

 

 

 

2

2

 

 

2

 

 

 

2 2

 

 

 

ž

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

задача решена.

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

1.14.

H

 

M

R

 

 

АЭ uвх

L

C

 

L0

K

 

 

 

 

П

Рис. 1.14. Радиотехническая схема

85

Решение задачи.

Схема

состоит из

источника

питания

П

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

(АЭ),

 

M ,

регулирующего ток нагрузки

 

H . Будем считать активный

элемент

 

 

 

характеристикой

 

 

 

и

с

большим

линейным

прибором

с

 

$

 

H

 

вх

 

током

во

входным

сопротивлением, позволяющим $

пренебречь

 

•¢

 

 

 

 

 

 

входной цепи.

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

схемы:

 

&

 

D

 

 

 

 

 

 

&D

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D D 0,

 

1.277

 

 

 

9 &

 

 

 

&

 

 

ž

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

D

, D

 

 

 

 

 

 

 

 

 

 

$

 

 

 

$

 

 

 

K

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где K

H

 

заряды токов

 

K и

 

 

Н

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

 

 

Поскольку

 

 

 

 

 

 

$

 

 

 

a<

 

&

 

 

1.278

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

вх

 

 

 

&$K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

H

a<

&DK

,

 

 

 

1.279

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

то уравнение (1.277) можно переписать так:

1

 

 

 

 

 

 

&

 

 

D

 

 

 

 

a<

 

&D

 

 

 

 

1.280

 

 

 

9

&

 

 

 

 

ž !

&

 

ž D 0.

 

 

 

 

 

 

 

 

 

К

 

 

 

 

 

 

 

 

 

 

 

К

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К

 

 

 

Подставив в уравнение (1.280) функцию D< в виде

1.281

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

nKQ= ,

 

 

 

получим следующее выражение:

 

a<

 

 

1

 

 

 

 

 

 

 

 

 

 

 

K

 

9

 

 

 

 

 

 

 

 

1.282

 

 

 

 

 

n

 

 

 

 

 

ž ! ž! 0,

0,

которое справедливо в двух случаях: либо когда

nK

либо когда в

нуль обращаются скобки.

 

 

 

 

 

 

 

86

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

9

 

 

a<

1

 

 

1.283

 

ž

! ž 0.

 

 

Разрешив это уравнение

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

 

, получим

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

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

 

 

 

нарастании амплитуды колебаний, описываемых функцией (1.281):

 

 

 

N

a<

 

 

1.284

 

 

 

 

ž ,

 

 

поскольку при этом условии корни уравнения (1.283) будут иметь положительные действительные составляющие.

Область неустойчивости изображена на рис. 1.15. По оси абсцисс отложена величина h = sM/C.

R

h

Рис.1.15. Область неустойчивых колебаний (заштрихована) радиотехнической схемы

Физический смысл условия нарушения устойчивой работы схемы, изображенной на рис. 1.14, прозрачен – это малые потери в контуре.

87

Часть 2. МЕТОДЫ АНАЛИЗА МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ СИСТЕМ

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

2.1.1. Собственное движение консервативных систем

Движение линейных систем с постоянными параметрами описывается линейными обыкновенными дифференциальными уравнениями с постоянными коэффициентами.

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

устройством системы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сначала рассмотрим собственное движение одномерной

консервативной системы, функция Лагранжа которой имеет вид

 

 

 

 

9

#H

 

~H

 

,

 

 

2.1.1

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

H

отклонение системы от

ее положения в точке

H 0

;

#

инерционный коэффициент

 

(масса)

 

системы;

~

коэффициент

упругой силы.

 

 

 

 

 

 

 

 

 

 

 

 

 

Функции (2.1.1) соответствует уравнение Лагранжа

 

 

 

 

 

 

HL ω H 0,

 

 

 

 

2.1.2

где

 

ω

2

 

 

~

 

 

 

 

 

 

 

 

 

 

 

#.

 

 

 

 

2.1.3

Уравнение (2.1.2) называется уравнением гармонического осциллятора.

88

 

Его решения известны: если ω2 T 0, то

 

 

 

 

 

 

] ;

 

2.1.4

 

 

 

 

 

U cosω

 

sinω

\cos ω

 

 

если ω2 ^ 0, то

 

 

 

shω

\ch ω

 

] ;

 

 

2.1.5

 

с

с

 

 

U chω

 

 

 

 

 

 

, α

величины, определяемые начальными условиями.

здесь ,

 

,

 

 

 

Решения

(2.1.4) описывают

незатухающие

 

колебания системы

 

.

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

 

, которые она

около положения устойчивого равновесия в точке

 

совершает с амплитудой

 

 

, частотой

 

и

начальной фазой колебания

 

 

 

 

 

 

 

 

 

 

H 0

 

 

 

 

α

Полная энергия системы

 

 

~H2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

#H'

 

 

2

\

2

 

 

 

 

 

2.1.6

 

 

 

 

 

 

 

2

2

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сохраняется.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решения (2.1.5) описывают движение системы около точки

неустойчивого равновесия.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Условие

устойчивого

 

движения

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

 

системы в

данном

 

случае имеет вид

 

простого

 

неравенства:

 

 

 

когда

 

 

 

 

 

 

потенциальная энергия в точке равновесия

 

 

 

 

 

 

минимум.

 

 

 

 

имеет ω

 

¤ 0,

 

 

В линейных системах искомое решениеHx(t)

удобно представлять в

 

0

 

 

 

 

 

 

 

виде комплексных функцийH Re 3Q > ,

 

 

 

 

 

 

 

 

 

 

2.1.7

где j – мнимая единица; А – в общем случае комплексная амплитуда:

3 Q α.

2.1.8

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

89

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

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

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

 

 

 

 

 

` 2

 

 

 

 

U'U';

 

 

2.1.9

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

,

 

 

U

 

U

 

,

 

 

2.1.10

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

H 0

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

(i=1,…,n)

отклонения системы от ее положения в точке

где

 

H, i=1,…,n.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ее функции Лагранжа

 

 

 

 

U U Q

 

 

U

U

 

2.1.11

 

 

 

9 2

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

соответствует набор из n дифференциальных уравнений

2.1.12

 

 

 

 

UL Q

 

U

0 $ 1, … , % ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

описывающих движение системы.

 

 

 

 

 

 

 

 

 

 

 

 

Как известно, решениями таких систем дифференциальных

уравнений служат экспоненты

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.1.13

 

 

 

 

 

H

3 Q> $ 1, … , % .

 

 

 

 

 

 

 

 

 

 

 

 

90