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

Основы_акустики_Гринченко_Вовк

.pdf
Скачиваний:
52
Добавлен:
10.04.2023
Размер:
16.81 Mб
Скачать

ωI = ω20 2δ2 , ωII = ω0 .

(2.56)

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

A(ω ) =

F0

,

A (ω ) =

F0

.

(2.57)

 

 

I

ΩR

 

υ

II

R

 

 

 

 

 

 

Здесь Ω = ω20 − δ2 собственная частота свободных колебаний сис-

темы с демпфированием.

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

Итак, для системы с демпфированием имеем три характерных частоты, а именно:

Ω = ω20 − δ2 собственных колебаний;

ωI = ω20 2δ2 резонанса по смещению;

ωII = ω0 резонанса по скорости.

Во многих практически важных случаях демпфирование в системе незначительное: δ2 << ω20 , поэтому количественные различия между

указанными характерными частотами малы. По этой причине в даль- нейшем, без специальных оговорок, мы будем говорить о резонансной частоте осциллятора, как о частоте ω0.

Для характеристики резонансных явлений важно обратить внима- ние на фазовые соотношения. Согласно формуле (2.47) для любой час- тоты ω имеем sin ϕ < 0, т.е. π < ϕ < 0. Итак, отклик всегда отстает по фазе от внешней силы. Из выражения (2.47) имеем:

tgϕ = −

2δω

.

(2.58)

2 2

 

ω

− ω

 

0

 

 

 

Если ω << ω0 , то tgϕ ≈ ϕ ≈ −2δω ω02 ;

если

ω >> ω0 , то tgϕ ≈ 2δ ω и

ϕ ≈ −π + +2δω. Таким образом, вне области резонанса колебания ос- циллятора происходят или в фазе с внешним воздействием ( ω << ω0 ), или в противофазе с ним ( ω >> ω0 ). Переход из одного состояния в

другого при изменении частоты воздействия происходит в узкой по- лосе шириной порядка 2δ вблизи резонансной частоты ω0. Точно на резонансе (ω = ω0) сдвиг фазы равняется ϕ = π/2.

51

Рис. 2.9. Нормированные час-

Рис. 2.10.Частотные зависи-

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

мости ϕ(ω) и ϕυ(ω)

смещения (кривая 1) и скорости

 

осциллятора (кривая 2)

 

Явление резонансного усиления колебаний можно осознать, если об- ратить внимание на сдвиг фазы ϕυ между внешней силой и скоростью. Согласно формуле (2.52) ϕυ = ϕ + π/2. На рис. 2.10 приведены частотные зависимости ϕ(ω) и ϕυ(ω). Если ω ω0, то между внешней силой и скоро- стью существует некоторый сдвиг фазы. Поэтому в пределах некоторой части каждого периода внешняя сила направлена против скорости, т.е. старается замедлить движение, вместо того, чтобы ускорить его. На ре- зонансе (ω = ω0) сдвиг фазы между внешней силой и скоростью ϕυ будет равен нулю (рис. 2.10), таким образом, сила всегда действует в направ- лении движения, постоянно подталкиваяего. Здесь есть наилучшая согласованность между внешней силой и внутренними свойствами ко- лебательной системы.

Пример 2.4. Осциллятор двигается по закону x(t) = x0sin(ωt), а его возбуждающая сила определяется зависимостью F(t) = F0cos(ωt). Найти коэффициент затухания осциллятора. Масса осциллятора m.

Решение. В режиме установившихся колебаний при F(t) = F0cos(ωt) смещение осциллятора определяется законом (см. (2.46)) x(t) = A(ω)cos(ωt + ϕ). По условию задачи x(t) = x0sin(ωt). Отсюда понят- но, что сдвиг фазы ϕ = π/2. Такой сдвиг фазы между смещением и силой бывает при условии резонанса, т.е. когда частота внешней си- лы равняется собственной частоте осциллятора ω = ω0. Согласно фор- муле (2.53) амплитуда смещения на резонансной частоте A = F0 /(2δωm). Отсюда коэффициент затухания δ = F0 /(2x0ωm).

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

52

источника с неограниченной мощностью. Поэтому согласно формуле (2.53) на резонансе (при ω ω0) и при условии δ 0 (осциллятор без демпфирования) амплитуда колебаний системы стремится к беско- нечности, т.е. система потребляет неограниченную энергию. Тем не менее, решение (2.51), (2.53) записано в соответствии с предположени- ем о наличии установившихся колебаний с постоянной амплитудой, что не имеет место при таких предположениях. В этом случае нужно быть более осторожными.

Рассмотрим осциллятор без демпфирования (δ 0), на который в момент времени t = 0 начинает действовать сила F(t) = F0cos(ωt), при

начальных условиях ξ(0) = 0,

 

ξ(0) = 0 . Согласно (2.47)—(2.49) при δ = 0

имеем такое выражение для общего движения системы:

ξ(t) = a1 cos(ω0t) +a2 sin(ω0t) +

F0

 

cos(ωt).

m(ω02

− ω2 )

 

 

 

 

 

 

Определяя a1 и a2 из начальных условий, получаем

 

ξ(t) =

F0

 

[cos(ωt) cos(ω0t)].

m (ω02

− ω2 )

 

 

 

 

 

 

Как видим, колебание системы состоят из двух гармонических дви- жений (вынужденного и собственного), причем амплитуды этих коле- баний одинаковы. Поэтому результирующее колебание в общем случае не будет гармоническим. Раскачивание (амплитуда колебаний) систе- мы зависит от соотношения частот ω и ω0. При ω ω0 нужно рас- крыть неопределенность типа 0/0. Таким образом, получим

 

F

 

cos(ωt)

cos(ω t)

 

F t

 

ξ(t) = lim

 

0

 

 

0

 

=

0

sin(ω t).

 

 

ω2

− ω2

2ω0m

ω→ω0

m

 

 

 

0

 

 

0

 

 

 

 

 

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

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

53

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

2.2.3. Добротность

При изучении колебательной системы без демпфирования были определены три интегральные характеристики: амплитуда, час- тота собственных колебаний, начальная фаза. Для системы с демп- фированием введем четвертую характеристику, которая служит ко- личественной мерой демпфирования, добротность колебательной системы Q . Рассмотрим вынужденные колебания системы в устано-

вившемся режиме. В соответствии с формулой (2.54) амплитуда коле- бательной скорости

Aυ (ω) =

 

ωF0

 

.

(2.59)

 

(ω02 − ω2 )2

 

m

+ (2δω)2

 

Рис. 2.11. Частотная зависимость скорости осциллятора

Амплитудно-частотная характеристика Аυ(ω) приведена на рис. 2.11. Здесь частотный интервал ω1 ω ω2 определяет частоты, на которых кинетическая энергия системы не более, чем в 2 раза мень- ше, чем на частоте ω0. Интервал Δω = ω2 − ω1 называют шириной резо-

нансной кривой. Добротность системы определяется в виде

 

Q =

ω0

 

 

 

 

.

(2.60)

 

ω

− ω

 

2

1

 

 

Выразим величину

ω = ω2 ω1

через параметры колебательной

системы. Для этого,

подставив

в уравнение

Aυ (ω0 )/ 2 = Aυ (ω)

выражение (2.54), получим соотношение ω20 − ω2 = ±2δω . Решение этих

уравнений дает такие корни: ω = −δ ± δ2 + ω20 ; ω = δ ± δ2 + ω20 . Физи- ческое содержание имеют лишь положительные корни. Поэтому пред-

54

полагаем,

что

ω1 = −δ + δ2 + ω02 ,

ω2 = δ + δ2 + ω02 .

Тогда

Δω = ω2 − ω1 = 2δ и добротность имеет такой вид:

Q =

ω

,

 

2δ =

R

(2.61)

0

 

.

 

2δ

 

 

 

m

 

Итак, добротность системы это отношение двух ее временных характеристик. Уменьшение амплитуды собственных колебаний сис- темы характеризуется постоянной времени затухания τ = 1/δ (см. (2.38)). Учитывая этот факт, из формул (2.60) и (2.61) имеем важное соотношение между шириной резонансной кривой вынужденных ко- лебаний и постоянной времени затухания свободных колебаний:

ωτ = 2.

(2.62)

Как видим, они связаны обратно пропорциональной зависимостью. Это есть общий результат в том понимании, что указанная зависи- мость характерна для колебательных систем разной физической при- роды. Уравнение (2.62) имеет практическое значение. Часто экспери- ментально легче изучить поведение системы вблизи резонанса в ре- жиме вынужденных колебаний, чем наблюдать время затухания. В этом случае, определив ω, из формулы (2.62) легко найти τ.

Запишем одно полезное соотношение между логарифмическим декрементом затухания θ и добротностью Q:

θ = δT =

2πδ

=

2πδ

=

2π

 

.

(2.63)

Ω

ω02 − δ2

(2Q )2

 

 

 

 

1

 

При Q = 1/2 имеем θ . Это наблюдается на границе колебатель- ного режима. При Q < 1/2 собственные движения в системе не коле- бательные. В случае высокой добротности (Q >> 1) (2.63) заменяется простым приближенным соотношением

θ ≈

π

, Q

π.

(2.64)

Q

 

 

θ

 

В качестве примера приведем порядок величины добротности неко- торых типов осцилляторов:

электрический контур — 50…500; громкоговоритель (на низких частотах) — 3...10…10; рояльная или скрипичная струна — 1000; камертон — 3000; кристалл кварца — 500 000.

55

Формулы (2.53), (2.54) для амплитуды смещения и скорости осцил- лятора можно переписать, используя понятие добротности, в таком виде:

A(ω) =

 

 

 

 

 

 

 

ξст

 

 

 

 

 

 

 

 

 

 

,

 

 

 

(2.65)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω2

2

 

 

ω2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω0

 

 

 

ω0

 

 

 

 

 

 

 

 

 

 

 

A (ω) =

 

 

ωξст

 

 

 

 

 

 

 

,

 

 

ξ

ст

=

F0

,

(2.66)

 

 

 

 

 

 

 

 

 

 

 

 

 

υ

ω2

2

 

 

ω2

 

 

1

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω0

 

 

 

 

ω0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а также согласно формуле (2.58)

и соотношению

 

ϕυ

π

 

tgϕ = tg

2

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −1 tg ϕυ выражение для сдвига фаз между F(t) и ξ(t) в таком:

 

 

 

 

 

ω

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

tgϕυ = Q

 

0

 

 

 

 

, ϕ = ϕυ

 

.

 

 

(2.67)

 

ω

 

2

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.12. Частотные зависимости амплитуды смещения (а), амплитуды ско- рости (б), сдвига фазы между внешней силой и скоростью осциллятора (в)

На рис. 2.12 приведены графики зависимостей А, Аυ, ϕυ от частоты внешней силы ω при различных значениях добротности Q . Кривые на

рис. 2.12, а, б показывают, что чем выше добротность системы, тем остреерезонансная кривая, т.е. вблизи резонансной частоты ω0 пове- дение графиков определяется демпфированием. За пределами резо- нансной зоны поведение графиков совпадает, что связано с малым

56

влиянием демпфирования на колебательный процесс в окрестности этих частот. Эти характерные особенности осциллятора мы еще рас- смотрим в следующем параграфе.

Фазовые соотношения отображены на рис. 2.12, в. Как видим, сдвиг фазы между скоростью и силой определяется неравенством −π2 < ϕυ < π2; сдвиг

фазы между откликом и силой изменяется в пределах –π < ϕ < 0, т.е. отклик всегда отстает от воздействия. Еще раз укажем, что для осцилляторов с Q >> 1 (скажем, Q > 10) при ω << ω0 колебания осциллятора происходят практически в фазе с внешним воздействием, а при ω >> ω0 — в противофазе. Переход с одного состояния в другое при изменении частоты внешнего воздействия наблюдается в довольно узкой полосе вблизи резонансной частоты ω0. Точно на резонансе (ω = ω0) сдвиг фазы ϕ = –π/2, а ϕυ = 0. Укажем, что при увеличении добротности кривые на рис. 2.12, в все больше приближаются к функции-ступеньке. Для осциллятора без демпфирования (Q )

ϕ = 0 и ϕυ = π/2 если ω < ω0 и ϕ = –π, ϕυ = –π/2 если ω > ω0.

Пример 2.5. Монокристалл сапфира в вакууме при низкой темпе- ратуре имеет добротность Q = 108…109... Частота собственных коле- баний монокристалла ω0 = 104 c–1. Оцените, во сколько раз изменится амплитуда свободных колебаний за сутки.

Решение. Добротность Q = ω0/(2δ), отсюда коэффициент затухания δ при Q = 108 равен δ = 5 10–5 с—1, а при Q = 109 имеем δ = 5 10–6 с–1. Амплитуда колебаний кристалла через сутки, т.е. за промежуток вре- мени t = 8,64 104 c, изменится в exp(δt) раз. При Q = 108 это составля- ет приблизительно 75, а при Q = 109 приблизительно 1,5 раза.

Амплитуду колебания системы на частоте ω0 называют динамиче-

ским откликом системы ξд, т.е. ξд = А(ω0). Тогда из формулы (2.65)

будем иметь соотношение:

ξд = ξстQ.

(2.68)

Таким образом, амплитуду установившегося вынужденного колебания на частоте резонанса ω0 можно определить, умножив величину пере- мещения при статической нагрузке ξст на добротность системы Q. По- нятно, что когда добротность системы Q >> 1, то амплитуда резонанс- ных колебаний может достигать значительных величин! Система, ус- тойчивая к внешнему статическому влиянию, может быть разрушена, если внешняя сила удерживает частотные компоненты, которые сов- падают с резонансной частотой системы. Именно поэтому необходи- мым есть точный расчет резонансных частот разнообразных сложных колебательных систем; например, мосты, корпусы самолетов и кораб- лей, сложные электрические схемы, мощные электроакустические преобразователи.

57

2.2.4. Комплексное механическое сопротивление

Из анализа поведения колебательной системы под влияни- ем внешней периодической силы видно, что реакция системы зависит не только от амплитуды внешней силы, но и от частоты воздействия. Как новая важная интегральная характеристика процесса вынуж- денных колебаний системы используется комплексный механический импеданс (от латинского слова impedio препятствую), или ком-

плексное механическое сопротивление. Эта величина определяется как отношение комплексных амплитуд силы и скорости в режиме ус- тановившихся колебаний системы. Если учесть (2.45) для ξ(t), то вы- ражение для комплексного механического импеданса, который обо- значается как Z, Н с м–1 = кг с–1, будет иметь вид

 

F

 

m ω2

− ω2 i2δω

 

2

2

 

 

 

 

 

 

0

 

 

ω

− ω

 

 

Z =

 

=

 

 

 

 

= m 2δ + i

0

 

.

(2.69)

ξ

 

 

 

iω

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если частота ω20 = Кm и характеристика демпфирования 2δ = R/m, то последнее выражение записываем так:

Z = R + i

K

− ωm .

(2.70)

 

ω

 

 

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

Z = Z exp(iϕυ), где

 

 

 

 

 

 

 

F

 

 

= m (2δ)2

 

ω2

2

 

 

 

ω2

− ω2

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

=

 

 

 

 

 

+

0

− ω , tgϕυ =

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

 

 

 

 

ω

 

 

 

 

2δω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

K

 

 

2

 

 

1

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

=

 

 

R

 

+

− ωm ,

tgϕυ

=

 

 

− ωm .

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

ω

 

 

(2.71)

(2.72)

Итак, модуль Z равен отношению амплитуд внешней силы и колеба- тельной скорости, а величина ϕυ определяет разность фаз между внешней силой и колебательной скоростью в системе.

Рассмотрим зависимость импеданса системы от частоты внешней силы ω:

 

 

2

 

2

 

 

 

• если ω

мала

ω0

>> ω; m

ω0

>> R , то

Z i K

. Импеданс чисто

 

 

 

ω

 

ω

 

ω

 

 

 

 

 

 

 

 

 

мнимый и положительный; он обратно пропорционален частоте ω и

58

не зависит от R. Такой импеданс называют импедансом упругого ти- па и говорят, что система управляется упругостью;

• если

ω2

− ω ≈ 0 ,

 

 

ω2

 

 

,

то

0

m

0

− ω

<< R

 

ω

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тельный и зависит только от R . Говорят, демпфированием;

Z R. Импеданс действи-

что система управляется

• если ω

 

2

 

велика ω >>

ω0

; mω >> R , то Z iωm. Импеданс чисто

 

 

ω

 

 

 

 

 

мнимый и отрицательный; он пропорционален частоте ω и не зави- сит от R . Это импеданс массового типа, и говорят, что система управляется массой.

Полученные результаты можно пояснить наглядными физически- ми соображениями. Если частота внешнего действия ω мала по срав- нению с частотой ω0, то в левой части уравнения (2.39) главную роль сыграет только слагаемое Кх и потому Kx F0cos(ωt). Внешняя сила тратится главным образом на преодоление упругой силы. Смещение совпадает по фазе с внешней силой.

Наоборот, когда частота ω значительно больше частоты ω0, то глав- ную роль играет слагаемое (mx ) и mx F0 cos (ωt ). В этом случае

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

Вобласти резонанса слагаемые (mx) и (Кх), хотя и велики каждый

вотдельности, но приблизительно равны по величине и противопо- ложны по знаку. Действительно, поскольку сдвиг фазы между внеш-

ней силой и смещением ϕ π/2 при ω ω0, колебания происходят по

закону x(t) = Asin(ωt). Если

ω

близка к

частоте

ω0 =

K m , то

слагаемые mx = −mω2A sin(ωt)

и

Kx = KA sin(ωt) в

уравнении (2.39)

компенсируют друг друга. Таким образом,

Rx F0 cos(ωt) .

Внешняя

сила направлена только на преодоление силы демпфирования. Ампли- туда скорости Aυ(ω0) = F0/R и, если демпфирование мало, то Aυ(ω0) ве- лика; скорость совпадает по фазе с внешней силой.

Проведенный анализ свойств величины Z наглядно отображается в частотных характеристиках рис. 2.12. Как видим, в окрестности час- тоты резонанса ω0 поведение осциллятора определяется демпфирова- нием, т.е. величиной добротности Q. За пределами резонансной зоны поведение графиков совпадает. Это есть принципиальный момент, ведь вне резонансной зоны влияние демпфирования на колебательное движение практически несущественно.

59

2.3. Энергетические характеристики процесса колебаний

Рассмотрим свободные колебания. Энергия осциллятора равна сумме кинетической ЕК и потенциальной ЕП энергий. Кинети-

ческая энергия EK = mξ2 2 , а потенциальная определяется работой упругих сил F = K ξ при отклонении системы от положения равнове-

ξ

сия ξ = 0 : EП = K ξdξ = K ξ2 /2. Наличие силы трения в уравнении

0

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

ная dE/dt, где E = EК + EП , должна быть отличной от нуля и отрица-

тельной. Действительно, dE /dt = ξ(mξ + K ξ), а используя

уравнение

mξ + Rξ + K ξ = 0 , получаем dE /dt = −(Rξ)ξ , где Rξ = Fтертя .

Итак, ско-

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

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

W = F v.

(2.73)

В каждом конкретном случае речь идет о произведении проекций вектора силы и обобщенной скорости на избранную координатную ось. При таком определении при вычислении значения W возможны два случая: W > 0 и W < 0, т.е. случай совпадения направления силы и скорости и случай, когда направления этих векторов противополож- ны. Физически это означает, что в первом случае энергия перетекает от источника к колебательной системе, а во втором от системы к источнику. Приведенное соотношение позволяет проанализировать поведение W со временем при произвольном законе изменения силы F(t). Детально остановимся на случае периодической внешней силы. Поскольку речь идет о вычислении квадратичных характеристик ко- лебательного смещения, то прямое использование комплексного ре- шения уже невозможно. Причина этого ясна и связана с тем, что: Re(P1P2 ) Re P1 Re P2 , где P1 и P2 некоторые комплексные величины.

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

60