Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
09 Вынужденные колебания.pdf
Скачиваний:
43
Добавлен:
11.05.2015
Размер:
277.26 Кб
Скачать

МИНИСТЕРСТВО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ

ТОМСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ СИСТЕМ УПРАВЛЕНИЯ И РАДИОЭЛЕКТРОНИКИ (ТУСУР)

Кафедра физики

ИЗУЧЕНИЕ ВЫНУЖДЕННЫХ ЭЛЕКТРОМАГНИТНЫХ КОЛЕБАНИЙ

Руководство к компьютерной лабораторной работе

2003

Министерство образования Российской Федерации

ТОМСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ СИСТЕМ УПРАВЛЕНИЯ И РАДИОЭЛЕКТРОНИКИ

(ТУСУР)

Кафедра физики

УТВЕРЖДАЮ Заведующий кафедрой физики

________________Е.М. Окс

________________

Физика

ИЗУЧЕНИЕ ВЫНУЖДЕННЫХ ЭЛЕКТРОМАГНИТНЫХ КОЛЕБАНИЙ

Руководство к компьютерной лабораторной работе

Разработчики: доценты каф. физики

___________Л.А. Троян

___________В.А. Бурдовицин

______________2003г.

Томск 2003

Введение

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

1 Краткая теория

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

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

Рассмотрим вынужденные колебания в последовательном контуре (Рис. 2.1), состоящем из R, L, C, возникающие под действием внешнего источника, э.д.с. которого изменяется по закону

E = Em cosωt

R L C

E

Рис. 2.1. Колебательный контур

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

условие

l<<λ ,

где l - размеры контура;

λ- длина волны распространяющихся электромагнитных колебаний.

Сучетом выше сказанного, на основании второго правила Кирхгофа можно запи-

сать

U L +UC +U R = E

 

 

(2.1)

Из уравнения (2.1.) получаем:

dI

 

g

 

 

 

 

L

+RI +

=

m

ω

 

 

dt

 

c

 

E

cos t

(2.2)

 

 

 

 

 

 

 

 

 

 

 

Используя соотношение

 

I =

dq

,

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

 

 

 

 

 

 

 

 

 

 

d2 q

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dq

 

 

2

(2.3)

 

 

 

 

 

 

 

dt

2

+2δ

 

 

 

 

+ω0 q=U m cosωt,

 

 

 

 

R

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

δ=

- коэффициент затухания,

(2.4)

 

 

2L

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

ω0 =

 

 

 

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

(2.5)

 

LC

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U m

=

Em

 

 

 

 

 

 

 

 

 

 

(2.6)

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (2.3) является дифференциальным уравнением вынужденных электрических колебаний в контуре. Данное уравнение является неоднородным. И как следует из теории дифференциальных уравнений решение его равно сумме общего решения, соответствую-

щего однородного уравнения

q1 (т.е. уравнения (2.3) без правой части) и частного реше-

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

q2

 

q =q +q

(2.7)

 

1

2

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

 

ω0

q=qm1 eδT cos(ω1t+ϕ) ,

(2.8)

где ω1 =

δ

 

2

2 - частота свободных затухающих колебаний.

 

Частное решение этого уравнения имеет вид

 

 

 

q2 = qm cos(ωt+ϕ),

(2.9)

где ω - частота изменения внешней э.д.с.;

ϕ - отстаивание по фазе вынужденного колебания от обусловливающей его вынуж-

дающей силы.

Слагаемое (2.8) играет существенную роль в начальной стадии процесса, при установлении колебаний.

Следовательно, вклад (2.8) в установившиеся колебания очень мал, и им можно пренебречь. Поэтому установившиеся вынужденные колебания будут определяться частным решением неоднородного уравнения

q=q2 =qm cos(ωt+ϕ)

 

 

(2.10)

Для нахождения значений qm и ϕ найдем производные q' и q''

 

q' =−ωqm sin(ωt+ϕ)=ωqm cos(ωt+ϕ+

π

)

(2.11)

2

 

 

 

 

 

q''=−ω2qm cos(ωt+ϕ)=ω2qm cos(ωt+ϕ+π)

(2.12)

 

Подставив (2.10), (2.11), (2.12) в (2.3), получим уравнение

ω2 qm cos(ωt +ϕ+π)=2δωqm cos(ωt +ϕ+π )+ ω02 qm cos(ωt +ϕ)=U m cos ωt

 

2

 

(2.13)

Из (2.13) следует, что постоянные

qm и ϕ должны иметь такие значений, чтобы гар-

моническая функция U m cosωt

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

стоящих в левой части уравнения. Для нахождения значений qm и ϕ воспользуемся ме-

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

против часовой стрелки. Длина этого вектора равна амплитуде колебаний соответствующей физической величины.

Если изобразить функцию ω02 qm cos(ωt+ϕ) вектором, длина которого ω0 qm ,

направленным вправо, то функция 2δωqm cos(ωt+ϕ+π) изобразится вектором

 

 

2

длиной

2δωqm повернутым относительно вектора ω02 qm против часовой стрелки на

угол π

(Рис. 2.2), а функция

ω2 qm cos(ωt+ϕ+π) - вектором длиной ω2 qm , по-

2

 

ω02 qm на угол π . Чтобы уравнение (2.13) выполня-

вернутым относительно вектора

лось, сумма трех перечисленных векторов должна совпадать с вектором, изображающим функцию U m cosωt .

Такое совпадение возможно, как видно из рис. 2.2, если

U m2 =(ω02ω2)2 qm2 +2δ2ω2 qm2

(2.14)

2дqm

 

Um

 

 

 

 

 

 

02 –щ2) qm

щ0 2qm

щ2qm

 

 

 

 

 

 

 

 

Рис.2.2. Векторная диаграмма

Из (2.14)

найдем амплитуды

 

 

 

 

 

 

 

 

 

 

q

=

 

 

U m

 

 

 

 

 

 

 

 

 

 

 

 

(2.15)

 

 

2

 

 

 

2 2

 

 

 

 

 

 

 

 

 

 

 

 

 

2 2+2

 

 

 

 

 

 

 

 

 

 

 

 

 

m

(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

)

δ ω

 

 

 

 

 

 

 

 

 

 

 

 

ω0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.2.2 позволяет получить значения ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tgϕ=

 

2δω

 

 

 

 

 

(2.16)

 

 

 

 

 

 

 

 

 

2 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ω0

ω )

 

 

 

Учитывая, что ω02 =

1

;

δ=

 

R

;

U m =

Em

,

 

 

 

 

 

 

 

 

 

2C

 

 

 

 

2L

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

Em

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2+(ωL1

)2

 

 

Получаем

 

 

 

 

 

qm ω

 

(2.17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ωC

 

 

 

 

 

 

 

 

 

 

 

tgϕ=

1

 

R

 

 

(2.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ωL

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ωC

 

 

 

 

 

 

 

Продифференцировав выражение (2.17) по t , найдем силу тока в контуре при установившихся колебаниях:

I =

dq

=−ωqm sin(ωt

+ϕ)=I m cos(ωt+ϕ+

π

),

(2.19)

 

 

 

dt

 

 

2

 

 

 

 

 

 

где I m - амплитудное значение установившегося тока в контуре

 

 

 

I m

Em

 

 

 

 

 

R2+(ωL1 ) 2

 

 

 

 

 

=ωqm=

ωC

 

 

(2.20)

 

 

 

 

 

 

В соответствии с (2.19) напряжение на резисторе

 

 

 

 

 

U R =RI m cos(ωt+ϕ+π)

 

 

(2.21)

 

 

 

2

 

 

 

Разделив выражение (2.10) на емкость, получим напряжение на конденсаторе

U C =

qm

cos(ωt+ϕ)=U cm cos(ωt+ϕ)

 

(2.22)

 

 

 

C

 

 

 

 

 

 

Здесь

qm

 

 

 

 

 

 

=

=

U m

 

=

I m

 

C

ωC R2+(ωL1

) 2

ωC

 

U cm

 

 

(2.23)

 

 

 

 

 

 

 

 

 

ωC

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

 

 

U L =L

dI

=−ωLI sin(ωt+ϕ+

π )=U Lm cos(ωt+ϕ+π)

(2.24)

dt

 

 

2

 

Сопоставляя выражения (2.21), (2.22), (2.23) видим, напряжение на емкости отстает по фа-

зе от силы тока на π

2

, а напряжение на индуктивности опережает ток на π

2

. На-

 

 

 

пряжение на активном сопротивлении изменяется в фазе с током.

 

 

Т.о., изменение I

, q , U C , U L , U R - представляют собой гармонические ко-

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

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

резонансом, а соответствующая частота резонансной частотой.

Чтобы определить резонансную частоту ωРЕЗ , любую из выше перечисленных гармонически изменяющихся величин, нужно исследовать функцию, определяющую данную величину на максимум. Так из (2.14) видно, что максимальное значение gm будет

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

ωРЕЗ :

4(ω02ω2)+8δ2ω=0

Данное уравнение имеет три решения ω=0 и ωω022δ 2 .

Решение ω=0 соответствует максимальному значению знаменателя, а отрицательное значение не имеет смысла. Т.о. резонансная частота для g и напряжение U C определя-

ются выражением:

ωgРЕЗ

=ωU C РЕЗ =

ω022δ2 =

1

R2 ω0

(2.25)

 

 

 

LC

2L2

 

На рис 2.3 приведены резонансные кривые для U С .

 

UCm

 

 

 

 

 

 

 

 

δ1

 

 

 

 

 

 

δ123

 

 

 

δ2

 

 

 

Um

 

δ3

 

 

 

 

 

 

 

 

 

δ4

 

 

 

 

 

 

ω0

 

ω

 

 

Рис. 2.3. Резонансные кривые для U С

 

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

шом затухании, когда 2δ2 >ω02 , резонанс не наблюдается, с увеличением частоты амплитуда вынужденных колебаний монотонно убывает (см. нижнюю кривую при δ4 ). При

стремлении

ω к нулю

резонансные кривые сходятся

в одной точке с ординатой

U cm =U m .

Величина U m

равна напряжению, которое возникает на конденсаторе при

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

U m . При малом затуханиии

(δ2 <<ω02 ) резонансную частоту можно брать равной ωРЕЗ =ω0 .

Резонансные кривые для силы тока приведены на рис. 2.4. Как видно из выражения

I m

 

 

δ1

 

 

 

 

 

δ1 <δ2 <δ3

 

δ3

 

δ2

 

 

 

 

ω0 ω

Рис. 2.4. Резонансные кривые для силы тока

(2.20), амплитуда тока будет иметь максимальное значение при постоянном R , если

Lω=ω1C

Следовательно, резонансная частота для силы тока совпадает с собственной частотой контура ω0

ωI РЕЗ =ω0 =

1

(2.26)

LC

 

 

При ω=0 , ток в цепи равен нулю, т.к. при постоянном напряжении ток в цепи с кон-

денсатором протекать не может.

При таком соединении элементов как рис. 2.1, когда ωРЕЗ ω0

амплитуда силы тока I =

U m

, а разность фаз

ϕ=0 , т.е. в цепи как бы нет ни емко-

 

 

R

 

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

положными по фазе. При этом напряжения U L и U C могут значительно превышать напряжение, приложенное к цепи U m .

Поэтому этот резонанс, наблюдаемый в цепи переменного тока, содержащей последовательно включенные L, C , R и источник э.д.с., называют резонансом напря-

жений.

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

Q=δT =

R 2π

=

πR

(2.27)

 

 

 

0

Lω0

 

2L

 

ω

и добротность Q . Добротность с точностью до множителя

2π равна отношению

энергии, запасенной в колебательной системе в данный момент W ЗАП к убыли этой энергии за период колебаний W

Q=2π

W ЗАП

(2.28)

W

 

 

Поскольку энергия W пропорциональна квадрату амплитуды колебаний, а ам-

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

А(t)=Ameδt , то

 

Q=

 

2π A2(t)

 

=

 

2π

 

 

 

 

 

 

 

 

 

.

 

A2(t)A2(t+T)

1e2δT

 

При малых затуханиях, когда δ2 <<ω02

 

 

 

 

 

 

 

 

 

1e2δT 2δT =2Θ

 

 

 

 

 

 

 

 

 

 

С учетом этого

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

π

 

2πL

 

Lω0

 

 

1 L

(2.29)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q=δT

=Θ

= RT =

R

=R C

 

Отсюда видно, что чем меньше затухание колебаний, тем больше добротность, а чем больше добротность, тем больше амплитуда резонансного напряжения на конденса-

торе U СРЕЗ в сравнении с амплитудой внешнего напряжения U m . Убедимся в этом, для чего найдем отношение амплитуд

U стРЕЗ

1

=

LC

 

1 L

=Q

 

U m

=ω0CR

CR

=R C

(2.30)

 

 

 

 

 

 

 

 

 

 

Таким образом, добротность контура определяет также остроту резонансных кривых. Чем больше Q , тем уже резонансная кривая. Если резонансную кривую нормировать относи-

тельно максимального значения и построить зависимость I m = f (ω) , то

I тРЕЗ