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

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

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

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

F = −R x.

(2.25)

Здесь знак минусуказывает на тот факт, что сила направлена про- тив вектора скорости; R коэффициент демпфирования (сопротив- ления, трения), H c/м.

Если теперь учесть, что в колебательной системе (см. рис. 2.1,а) на массу действует сила сопротивления (2.25), то второй закон Ньютона будет иметь вид:

mx + Rx + Kx = 0.

(2.26)

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

ξ + 2δξ + ω02ξ = 0, 2δ = R/m, c–1.

(2.27)

Здесь ω0 собственная частота системы при отсутствии демпфиро- вания; δ коэффициент затухания.

Конкретное движение колебательной системы определяется зада- нием начальных условий (2.12). Задача нахождения решений обыч- ных дифференциальных уравнений с постоянными коэффициентами (именно таким есть уравнения (2.27)) относительно простая. Она легко сводится к задаче нахождения корней алгебраических уравнений, степень которых совпадает с порядком дифференциального уравне- ния. При этом используется тот факт, что экспоненциальная функция инвариантна относительно операции дифференцирования. Если, в

41

качестве пробного решения, в уравнение подставить искомую функ- цию ξ(t) в виде

ξ(t) = A exp(αt ),

(2.28)

(где A и α произвольные величины), то для α будем иметь уравне- ние:

 

 

α2 + 2δα + ω02 = 0,

 

 

(2.29)

которое в общем случае имеет корни

 

 

 

 

 

α1,2 = −δ ± δ2 − ω02 .

 

 

(2.30)

Понятно, что возможные три случая:

 

 

 

1) если

δ2 > ω02 ,

то

α1 < 0 ,

α2 < 0

и

решение

ξ(t) = A1 exp(α1t )+ A2 exp(α2t )

не будет описывать колебательное дви-

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

2) если δ = ω, то α1 = α2 = δ и ξ(t) = A1 exp(−δt )+ A2t exp(−δt ). В этом случае система также не будет осуществлять колебательное движение. Для акустических систем интересен третий случай, который соответ- ствует малым затуханиям, а именно:

3) если δ < ω0, то α1,2 = −δ ± i

ω02 − δ2 , и

 

 

 

 

 

 

 

2

2

 

 

i

2

2

 

(2.31)

ξ(t) = exp(−δt ) A1 exp i

ω0

− δ

t

+ A2 exp

ω0

− δ

t .

 

 

 

 

 

 

 

 

 

 

 

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

 

 

2

− δ

2

 

 

 

 

 

 

2

− δ

2

 

,

(2.32)

ξ(t) = exp(−δt ) a1cos

ω0

 

t

+ a2sin

ω0

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

− δ

2

 

 

,

 

 

 

 

(2.33)

ξ(t) = A exp(−δt )cos

ω

 

t + ϕ

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

42

где A = a2

+a2

, сos ϕ = a1/A, sin ϕ = –a1/A. Подставим (2.32)

в на-

1

2

 

 

 

 

чальные условия (2.12) и найдем константы:

 

 

 

a1 = ξ0, a2 =

υ0 + δξ0

.

(2.34)

 

 

 

 

 

 

ω2 − δ2

 

 

 

0

 

 

Рис. 2.6. Пример колебания затухающего осциллятора

Решение (2.33) описывает затухающее колебательное движение, амплитуда которого уменьшается согласно закону exp(–δt) (рис. 2.6). Понятно, что функция ξ(t), как и функция ξ(t), не является периоди- ческой, поскольку для любого T очевидно, что ξ(t) ≠ ξ(t +T ). Поэтому в

строгом понимании слова периода не существует. Тем не менее, время между двумя последовательными прохождениями системы через по- ложение равновесия (в одном направлении), или между двумя после- довательными максимальными отклонениями, есть постоянная вели- чина. Это обстоятельство позволяет ввести понятие периода для зату- хающего осцилляторного процесса:

T =

2π

,

Ω = ω02 − δ2 ,

(2.35)

 

Ω

 

 

 

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

43

 

Cn

 

=

 

A exp(−δt )cos (Ωt + ϕ)

 

= exp(δT ).

C

 

A exp −δ(t +T ) cos Ω(t +T )+ ϕ

 

 

 

 

n +1

 

 

 

 

 

 

 

Величину θ = δT

 

называют логарифмическим декрементом затуха-

ния и определяют из формулы

 

 

 

 

 

 

 

 

 

 

θ = ln

Cn

 

= δT .

 

(2.36)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

n +1

 

 

Величина θ характеризует затухание амплитуды колебаний за один пе- риод. Важной особенностью (2.36) есть постоянство отношения Cn Cn +1 (для любого n ) — это характерная особенность принятой на-

ми модели демпфирования.

Формулу (2.36) можно использовать для экспериментального опре- деления θ. Для этого необходимо определить из эксперимента отно- шения двух соседних амплитуд колебаний. Но большая точность дос- тигается при использовании отношения двух амплитуд, отдаленных

на N периодов. При этом (2.36) приобретает вид θ =

1

ln

Cn

.

N

 

 

 

C

 

 

 

n +N

Согласно выражению (2.33) демпфирование влияет на изменение ам- плитуды колебаний осциллятора и собственную частоту системы по срав- нению с отсутствием демпфирования. Но, и это является характерной чертой для колебательных систем, которые нас интересуют, влияние демпфирования со временем на амплитуду является более весомым, чем на период колебаний. Чтобы убедиться в этом, запишем формулу (2.36) в виде

 

2π

 

 

2π

 

 

2π

 

 

 

 

ω

2

θ = δΤ = δ

 

 

 

=

 

 

=

 

 

 

 

= 2π

 

0

 

1. (2.37)

 

 

 

 

 

 

 

 

 

 

2

− δ

2

 

2

 

 

 

2

 

 

 

 

Ω

 

 

 

ω0

 

 

ω20

1

 

 

ω0

 

1

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

δ

 

 

 

ω − Ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

Тогда, например, при условии Ω/ω0 = Т0/Т = 1/1,1, т.е. когда Т отли- чается от Т0 лишь на 0,1Т0, имеем θ 2,879 и exp(δT ) 17,8. За два периода колебаний амплитуда уменьшается почти в 316 раз, т.е. ко- лебания практически прекращаются.

Обозначим буквой τ промежуток времени, за которое амплитуда колебаний уменьшается в e раз. Тогда exp(−δτ) = e1 0,368 , откуда имеем

δτ = 1, δ = 1/τ.

(2.38)

Время τ называют постоянной времени затухания. Итак, δ вели-

чина, обратная времени затухания. Пусть N число периодов, после

44

которых амплитуда уменьшается в е раз. Тогда τ = NT и θ = δT = T/τ = 1/N. Таким образом, θ величина, обратная числу пе- риодов, после которых амплитуда уменьшается в e раз. Если, напри- мер, θ = 0,01, то это означает, что амплитуда уменьшается в e раз по- сле 100 периодов.

2.2.2. Вынужденные колебания

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

ной скорости ξ(t) системы, т.е. от характера движения системы. Как

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

Следует отметить, что в механических колебательных системах не так просто с технической точки зрения влиять периодической силой непосредственно на массу, которая двигается. Значительно проще это сделать в электрических или оптических колебательных системах, на- пример, в колебательном контуре, при подключении его к внешнему источнику переменного напряжения (рис. 2.7, б). Но нетрудно понять, что поддерживать колебание системы, которая изображена, напри- мер, на рис. 2.1, а, можно, не прикладывая непосредственно внеш- нюю силу F(t) к массе тела, достаточно эту силу приложить к свобод- ному кольцу пружины, как это изображено на рис. 2.7, а.

Итак, пусть правый конец пружины на рис. 2.7, а выполняет за- данное движение в горизонтальном направлении согласно закону x(t) . В состоянии покоя пружина не деформирована и x = 0 . Тогда

уравнение движения системы записывается аналогично выражению (2.2). Однако в данном случае упругая сила пропорциональна относи-

тельному

перемещению концов

пружины

(рис. 2.7,а), а именно,

x x(t) .

Итак, уравнение

движения

будет иметь вид

mx = −Rx K (x x(t)), или

 

 

 

mx + Rx + Kx = F(t),

(2.39)

где F(t) = Kx(t) . Имеем уравнение осциллятора, в правой части которо-

го вместо нуля стоит заданная функция независимой переменной F(t). Эта функция определяет внешнее силовое влияние на систему.

45

Рис. 2.7. Примеры осцилляторов под внешним воздействием

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

Lq + Rq + q /C = U (t).

(2.39а)

После деления (2.39) или (2.39а) на коэффициент при второй про- изводной, получаем уравнение осциллятора с демпфированием при внешнем воздействии. Для уравнения (2.39) будем иметь:

ξ + 2δξ + ω02ξ =

1

F(t),

(2.40)

m

 

 

 

где 2δ = R/m, а ω20 = K /m.

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

ней силы F(t) = F0cos(ωt):

2

F0

cos(ωt),

 

ξ + 2δξ + ω ξ =

 

(2.41)

 

0

m

 

 

 

 

 

где ω частота внешней силы.

В общем случае решение этого уравнения должно удовлетворять начальным условиям (2.12). Поскольку уравнение линейное, то его решение может иметь вид суммы общего решения однородного урав- нения и некоторого частного решения неоднородного уравнения.

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

46

cos(ωt) = Re (exp(iωt )) = Re(exp(iωt )).

(2.42)

Здесь обозначение Re() указывает на то, что выделена действитель-

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

ξ + 2δξ + ω2ξ =

F0

exp(iωt ),

(2.43)

 

0

m

 

 

 

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

ся получить непосредственно из комплексного решения.

Особое внимание нужно обратить на такое обстоятельство. Как видно из (2.42), cos(ωt) есть действительная часть двух разных экспо- нент, которые различаются знаком показателя степени. Если потом, анализируя решение, выделить действительную часть, то различие между двумя возможными решениями исчезает; если рассматривать решения в комплексной форме, то можно увидеть, что они сущест- венно различаются. В акустической литературе используются обе возможные формы комплексного решения. Важно лишь при этом быть последовательным, — выбрав одну из возможных форм зависи- мости, нужно сохранять ее на протяжении всего рассмотрения. В дальнейшем будем использовать временную зависимость exp(–iωt).

Частное решение (2.43) получаем в виде:

ξ(t) = A exp(iωt ).

(2.44)

После подстановки (2.44) в уравнение (2.43) легко определить ампли- тудную характеристику А, тогда решение приобретает вид

F0

 

ξ(t) = m(ω02 − ω2 i2δω) exp(iωt ).

(2.45)

Действительная часть комплексной функции (2.45) определяется вы- ражением:

 

F0

(ω02

− ω2 )cos (ωt )+ 2δωsin(ωt )

 

 

F

cos (ωt + ϕ)

 

 

ξ(t) =

 

 

 

 

 

 

 

=

 

0

 

 

,

(2.46)

 

 

 

− ω2 )

2

+ (2δω)2

 

m (ω02

− ω2 )

2

 

 

 

 

 

 

 

 

m (ω02

 

 

 

+ (2δω)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

47

cosϕ =

ω02 − ω2

, sinϕ = −

2δω

. (2.47)

(ω02 − ω2 )2 + (2δω)2

(ω02 − ω2 )2 + (2δω)2

После получения частного решения уравнения движения общее дви- жение системы при произвольных начальных условий описывается выражением:

ξ(t) = exp(−δt ) a

cos(Ωt) +a

2

sin(Ωt) + A cos(ωt + ϕ),

(2.48)

 

1

 

 

 

 

 

где

 

 

 

 

F0

 

 

2

2

, A =

 

 

 

 

Ω = ω

− δ

 

 

 

.

(2.49)

 

 

 

0

 

m (ω02 − ω2 )2 + (2δω)2

 

 

 

 

 

 

Качественный анализ (2.48) дает достаточно простую картину по- ведения системы. Но при этом анализе формируются те фундамен- тальные понятия, которыми описываются разные ситуации в теории колебаний и акустике. Именно с точки зрения раскрытия сути этих понятий мы и выполняем дальнейший анализ. Из формулы (2.48) прежде всего, видим неравноценный вклад отдельных слагаемых. Первый из них имеет экспоненциальный множитель, который обу- словливает затухание собственных колебаний системы. Поэтому мож- но выделить два интервала времени. Сначала существенными могут быть собственные колебания системы. Процесс постепенного затуха- ния собственных колебаний называют процессом установления коле- баний в системе. Для его количественной характеристики можно принять интервал времени, за который амплитуда собственных коле- баний уменьшается на 99 %. В дальнейшем движение системы полно- стью определяется вторым слагаемым выражения (2.48). Говорят, что система находится в режиме стационарных вынужденных колебаний.

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

Рассмотрим процесс установления колебаний в случае, когда в на- чале осциллятор покоится (ξ0 = 0, υ0 = 0 в условиях (2.12)), и в момент времени t = 0 на него начинает действовать сила F(t) = F0 cos(ω0t), т.е.

имеем ω = ω0. Из формул (2.47) следует, что вынужденные колебания имеют сдвиг фазы ϕ = π/2 относительно внешней силы, а согласно выражению (2.49) амплитуда вынужденных колебаний в установив- шемся режиме A = F0 (Rω0 ). Из начальных условий (2.12) находим по-

48

стоянные: а1 = 0, a2 = −ω0 A Ω

(получите самостоятельно). Тогда ре-

шение (2.48) можно записать в виде

 

ξ(t) = A

sin(ω t) ω0 exp(−δt ) sin(Ωt) .

 

 

0

Ω

 

 

 

 

 

Поскольку коэффициент затухания δ << ω0 (что практически всегда наблюдается для осцилляторов, которые нас интересуют), то частоты

Ω = ω20 − δ2 и ω0 мало отличаются одна от одной. Поэтому погреш-

ность, которая накапливается в фазе колебаний за время установле- ния вследствие такого изменения, будет малой. Итак, положив Ω ≈ ω0, в итоге получим такой результат:

 

F0

(

 

(

 

))

0

 

ξ(t) =

Rω0

1

exp

 

−δt

 

sin(ω t),

(2.50)

 

 

 

 

 

 

 

 

типичный вид зависимости (2.50) отображен на рис. 2.8. Как видим, в начале происходит рост амплитуды колебаний, которая в установив-

шемся режиме выходит на стационарный уровень A = F0 (Rω0 ). Время наращивания колебаний можно оценить так: δ–1.

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

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

49

В соответствии с (2.48) в режиме установившихся колебаний из- менение со временем обобщенной координаты представлено зависи- мостью

ξ(t) = A cos(ωt + ϕ),

(2.51)

для скорости имеем

 

ξ(t) = −Aωsin(ωt + ϕ) = Aωcos (ωt + ϕ + π/2).

(2.52)

Поскольку колебания осуществляются с частотой внешней силы, то анализу подлежат лишь две интегральные характеристики колеба- тельного процесса: амплитуда A и фаза ϕ. Величина ϕ есть сдвиг фа- зы между внешней силой F(t) и откликом на нее осциллятора ξ(t). Итак, величина ϕ + π/2 определяет сдвиг фазы между внешней си- лой F(t) = F0 cos(ωt) и скоростью в системе ξ(t).

Согласно формулам (2.49) и (2.52) зависимость амплитуд смещения A и скорости Aυ от частоты внешнего воздействия представляется

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

A(ω) =

 

 

F0

,

 

(2.53)

m (ω02 − ω2 )2 + (2δω)2

 

Aυ(ω) = ωA(ω) =

 

ωF0

 

.

(2.54)

 

 

 

 

m

(ω02 − ω2 )2 + (2δω)2

 

Отсюда легко найти значение скорости и смещения в предельных слу- чаях малых и больших частот. Если ω << ω0, то

A(ω)

F0

=

F0

= ξст.

(2.55)

 

K

mω02

 

 

 

Отклонение системы от положения равновесия стремится к значению статического отклонения ξcт под действием заданной силы, а ско- рость к нулю. Обратим внимание на то, что в качестве масштаба, относительно которого рассматривается малая или большая частота, здесь берем величину ω0, которая является собственной частотой со- ответствующей системы без демпфирования. Если, наоборот, ω >> ω0, то амплитуды смещения и скорости стремятся к нулю.

Особый интерес имеет определение тех частот, при которых ам- плитуды смещения и скорости достигают максимальных значений. Из

уравнений dA(ω) = 0, dAυ(ω) = 0 находим, соответственно: dω dω

50