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

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

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

возбуждения струн щипковых инструментов. Итак, начальные усло- вия (3.15) будут иметь вид

 

hx /ξ,

0 x ≤ ξ,

 

y (x,0) = Q1 (x ) =

 

 

(3.81)

h (l x )/(l − ξ), ξ ≤ x l.

 

 

 

 

y (x,0)

= Q2 (x ) = 0.

t

 

 

Рис. 3.17. Начальное отклонение струны

 

Согласно формуле (3.58)

находим коэффициенты an и bn :

 

a =

2hl2

 

sin

nπξ

,

b = 0, n =1,2,3,...

(3.82)

 

 

 

 

n

n2π2ξ(l

− ξ)

 

 

n

 

 

 

l

 

 

 

 

Подставив формулу (3.82) в выражение (3.55), получим искомое ре- шение

y (x,t ) =

 

 

2hl

2

 

1

 

 

nπξ sin

nπx

cos (ω

t ).

 

 

 

 

 

 

sin

(3.83)

 

 

 

 

 

 

 

 

 

 

 

 

π2ξ(l − ξ)n =1n2

 

 

l

 

 

n

 

 

 

 

l

 

 

 

 

 

 

 

Энергия n -ой гармоники согласно (3.80) равна

 

 

 

 

 

 

 

 

 

1

 

 

2 2

 

 

 

ρl3h2c2

 

 

2

nπξ

 

 

Hn

=

 

ρl

ωn

an

=

 

 

 

 

 

 

sin

 

 

 

 

(3.84)

4

n2π2ξ2

(l − ξ)2

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и уменьшается обратно пропорционально n2.

Если величины ξ и n такие, что sin(nπξl ) = 0, т.е. точка ξ является узлом n-ой гармоники, то an = 0 , и в полученном решении соответст-

вующая гармоника отсутствует. Например, если струна оттянута в средней точке ξ = l2, то nπξ/l = nπ/2 и для всех четных значений n

точка l/2 есть узел. Тогда в решении будут присутствовать только не-

четные гармоники. Итак, в случае ξ = l

2

решение (3.83) приобретает

вид

 

 

 

 

 

 

 

1

n

2n 1 πx

 

 

( )

sin (

 

) cos (ω

 

y (x,t ) = 8h

 

t ).

 

 

 

π2 n =1 (2n 1)2

 

l

 

2n 1

 

 

 

 

 

 

 

 

 

131

Рассмотрим второй пример. Струна возбуждается при задании начальной скорости υ0 на участке [ξ − δ,ξ + δ]. Такая ситуация соот-

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

0, 0 x < ξ − δ, Q1 (x ) = 0, Q2 (x ) = υ0, ξ − δ ≤ x ≤ ξ + δ,

0, ξ + δ < x l.

Из общих формул (3.58) находим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2l

υ

 

 

nπ(ξ − δ)

 

 

nπ(

ξ + δ)

 

 

 

 

a

= 0, b

=

 

 

 

0

 

cos

 

 

 

 

 

 

 

 

cos

 

 

 

 

 

 

,

n =1,2,...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

n

 

 

n

2

2

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

π c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, процесс колебаний описывается функцией

 

 

 

y (x,t ) =

 

4lυ

 

 

 

1

 

 

nπξ

 

 

 

nπδ

 

 

nπx

sin(ω

t ),

 

 

 

 

 

 

0

 

 

 

sin

 

 

 

sin

 

 

sin

 

 

 

(3.85)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π2c n =1n2

 

 

l

 

 

 

 

 

l

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

а энергии отдельных гармоник равны

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4ρlυ 2

2

nπξ

2

n

πδ

 

 

 

 

 

 

 

 

 

 

 

 

Hn =

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.86)

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

sin

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

n

2π2

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

В частном случае, когда начальные скорости приданы всем точкам струны (ξ = l2, δ = l2), закон движения струны имеет вид

 

4l

υ

1

 

(

2n

1 πx

 

 

y (x,t ) =

0

sin

 

)

sin(ω

t ),

 

 

 

 

 

 

 

π2c n =1 (2n 1)2

 

 

l

 

 

2n 1

 

 

 

 

 

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

(3.87)

Как ви-

дим, амплитуды гармоник An уменьшаются так: 1n2 , тогда распре- деление энергии между гармониками согласно (3.86) происходит по закону Hn 1n2 . Для этого частного случая энергия также сосредо-

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

Рассмотрим ситуацию, когда ширина молоточка 2δ стремится к нулю, т.е. струна возбуждается ударом очень острого молоточка, ко- торый дает струне импульс I в точке ξ. Это задача о влиянии на стру- ну сосредоточенного импульса. Сначала предположим, что импульс

132

равномерно распределен на

некотором малом

участке стру-

ны(ξ − δ),(ξ + δ), который окружает точку ξ, т.е. I = 2δρυ0 . Откуда

υ

=

I

.

(3.88)

 

0

 

2δρ

 

 

 

 

Далее, перейдя к пределу δ → 0, определим искомое решение. Подста- вив (3.88) в формулу (3.85), получим

y (x,t ) =

2lI 1

sin nπξ sin(nπδ l ) sin nπx sin(ω t ). (3.89)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

π2ρc n =1n2

 

δ

 

 

 

 

l

 

l

 

 

Поскольку

lim

sin(nπδ l )

= nπ

, то имеем такое решение задачи о коле-

 

 

δ→0

δ

l

 

 

 

 

 

 

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

y (x,t ) =

2I

1 sin

nπξ sin

nπx sin(ω

t ),

 

 

πρc n =1n

 

l

 

 

l

 

n

 

 

 

 

 

 

 

 

соответственно, энергии отдельных гармоник имеют вид

Hn

=

I 2

2

nπξ

 

sin

 

 

.

ρl

l

 

 

 

 

 

(3.90)

(3.91)

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

Рассмотренные примеры свидетельствуют о том, что на энергию гармоник существенно влияют значение координаты ξ, которая опре- деляет место удара молоточка, и ширина отрезка δ, на котором выпол- няется удар. Кроме того, наличие множителя sin(nπξl ) в формуле (3.86)

приводит к тому, что в случае, когда центр удара молоточка совпада- ет с узлом n-ой гармоники, то энергия соответствующей гармоники равна нулю. Эти выводы иллюстрируются рис. 3.18, на котором по- казаны относительные значения энергий гармоник Hn/Hmax при раз- ных величинах ξ и δ, здесь Нmax определяет наибольшую энергию со- ответствующей гармоники.

133

Рис. 3.18. Относительные значения энергии гармоник Hn /Hmax:

а ξ = 0,5l, δ = 0,2l ; б ξ = 0,3l, δ = 0,2l ; в ξ = 0,3l, δ = 0,05l

Указанные особенности учитывают при создании музыкальных инструментов. Наличие высоких гармоник (начиная с 7-ой) нарушает гармоничность звука и вызывает ощущение диссонанса (от латинско- го слова dissonans неблагозвучие). Наличие низких гармоник, на- оборот, вызывает чувство полноты звука. В рояле место удара моло- точка выбирают вблизи точки закрепления струны между узлами 7- ой и 8-ой гармоник, чтобы уменьшить их энергию. Регулируя ширину молоточка и его жесткость (для этого его ударную часть обтягивают кожей), стремятся увеличить энергию низких (3-ей и 4-ой) гармоник. В старых конструкциях рояля, которые имеют более резкий звук, да- же до некоторой степени звон, использовали узкие и жесткие моло-

точки [52, с. 143].

3.7. Вынужденные колебания струны

Рассмотрим струну длиной l, которая закреплена на кон- цах. Пусть функция g(x, t) характеризует внешнюю силу на единицу длины струны и не зависит от движения струны, т.е. имеем источник силы бесконечной мощности. Силы сопротивления отсутствуют (R = 0). Тогда уравнение движения струны приобретает вид (см. (3.6)):

ρ

2y

= F

2y

+ g (x,t ), или

1 2y

=

2y

+

g (x,t )

.

(3.92)

 

 

 

 

 

 

t2

x2

c2 t2

x2

 

 

 

 

 

 

F

 

Пусть внешняя сила периодическая во времени g(x, t) = g(x)exp(–iωt). Рассмотрим только вынужденные колебания; для этого необходимо найти частное решение уравнения (3.92). Колебание струны во вре-

134

мени будут отвечать действию внешней силы, поэтому частное реше- ние уравнения (3.92) будем искать в виде

y(x, t) = Y(x)exp(–iωt).

(3.93)

Подставим (3.93) в (3.92) и получим уравнение для амплитудной функции Y(x) искомого решения

d2Y

+ ω2 Y = −

g (x )

.

(3.94)

dx2

 

c2

F

 

Используя свойства нормальных колебаний струны, можно до- вольно просто найти решение уравнения (3.94). Действительно, по- скольку собственные формы колебаний струны sin(nπx/l) = sin(ωnx/c) образуют полную и ортогональную систему функций на отрезке [0, l], то мы имеем возможность, представить амплитудную функцию внешней силы g(x) рядом Фурье по собственным формам струны

{sin(ωnx/c), n = 1,2,…}в виде

g (x ) =

sin

ωn x

 

,

g

 

n

 

 

c

 

 

 

n =1

 

 

 

 

где коэффициенты ряда g

= 2 l

g

(x )sin

ωn

x

dx.

 

 

n

 

 

 

 

 

c

 

 

 

 

 

 

 

l 0

 

 

 

 

 

 

 

 

 

репишем (3.94) таким образом:

 

 

 

 

 

 

 

 

 

 

d2Y

 

ω2

 

1

 

 

 

 

 

ω

 

 

+

 

Y = −

 

 

g

n

sin

 

 

n x .

 

 

 

 

 

dx2

 

c2

 

F n =1

 

 

c

 

 

 

 

 

 

 

 

(3.95)

С учетом (3.95) пе-

(3.96)

Представим функцию Y (x ), которая определяет амплитуду колебаний

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

 

 

 

 

 

Y (x ) =

sin

ωn x

.

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n =1

 

 

 

 

 

 

 

 

Тогда уравнение (3.96) приобретает вид

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

1

 

 

 

 

 

 

y

 

ω

ωn

sin

ωn x

= −

 

g

sin

ωn x

.

 

2

 

 

 

n

 

 

c

2

 

 

 

 

 

 

n

 

 

c

 

n =1

c

 

 

 

 

c

 

 

F n =1

 

 

 

 

g c2

Откуда искомые коэффициенты yn =

n

 

. Итак, вынужденные

F (ω2 − ωn2 )

колебания струны определяются рядом

135

 

c2

g

 

ω

 

 

y (x,t ) = −

 

n

sin

n x exp(iωt ).

(3.97)

 

ωn2 − ω2

 

F

n =1

 

c

 

Решение в виде ряда (3.97) наглядно раскрывает природу явления резонанса. Согласно (3.97) для возникновения резонанса необходимо выполнение двух условий:

1) приближение частоты внешней силы ω к одной из собственных

частот струны ωn, т.е. ω ωn;

2) наличие в пространственном распределении внешней силы со- ставляющей, которая соответствует n-ой форме колебаний, т.е. gn 0. Например, когда ω ω2, то необходимо, чтобы выполнялось условие

f2 0.

Если эти два условия выполняются, то амплитуда движения во времени неограниченно возрастает (поскольку анализируем систему без демпфирования), и сумма ряда (3.97) определяется фактически одним слагаемым, для которого ω ωn. Физически это означает, что при резонансе струна приобретает форму колебаний, которая соот- ветствует гармонике с номером n.

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

3.8. Влияние демпфирования на движение струны

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

Рассматривая движение струны с учетом демпфирования, необхо- димо также обратить внимание на некоторые особенности в поведе- нии коэффициента демпфирования R. Физические процессы, кото- рые служат причиной появления демпфирования в колебательных системах, чрезвычайно разнообразны. Относительно рассматривае- мого случая движения струны можно вспомнить несколько сущест- венных механизмов рассеяния энергии. Это, прежде всего, процессы преобразования механической энергии в теплоту при деформации

136

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

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

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

g(x, t) = g(x) exp(–iωt).

(3.98)

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

(3.98):

1 2y

= −2δ(ω)y

+

2y

+

g (x )

exp(iωt ),

(3.99)

 

 

 

 

c2 t2

x2

F

t

 

 

 

 

где c2 = F/ρ, 2δ(ω) = R(ω)/F. Очевидно, что изменение во времени функции y(x,t) будет повторять зависимость от времени внешней си- лы, т.е. будет иметь вид y(x,t) = Y(x) exp(–iωt). Фазовые сдвиги при этом учитываются тем, что функция Y(x,t) комплексная. Для ее опре- деления из (3.99) получим обыкновенное дифференциальное уравне- ние с коэффициентом демпфирования, который соответствует за- данной частоте колебаний ω:

2

 

2

 

g (x )

 

 

d Y

+

ω

i2δ(ω) Y +

= 0.

(3.100)

dx2

 

F

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

137

Частное решение этого уравнения должен удовлетворять граничным условиям на концах струны. Именно поэтому при его построении це- лесообразно использовать свойства собственных форм колебаний идеальной струны sin(nπx/l), n = 1, 2, 3… Запишем ряды Фурье по собственным формам струны:

Y (x ) =

sin

nπ x

 

,

(3.101)

Y

 

n

 

 

l

 

 

 

 

n =1

 

 

 

 

 

g (x ) =

sin

nπ x

.

(3.101а)

g

 

n

 

 

l

 

 

 

 

n =1

 

 

 

 

 

Здесь gn известные при заданном распределении внешней силы ко- эффициенты:

gn

=

2 l

nπ

 

(3.102)

l

g (x )sin

l

x dx,

 

 

0

 

 

 

величины Yn должны быть определены. После подстановки формул

(3.101) и (3.101а) в уравнение (3.100) получим

 

 

2

 

n

π

2

 

gn

 

 

Y

 

ω

i2δ(ω)

 

= −

.

(3.103)

 

 

 

 

 

 

n

 

 

l

 

 

F

 

 

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это соотношение будет определять значения Yn через заданные ко- эффициенты gn разложения внешней силы в ряд Фурье. Таким обра- зом, вынужденное движение струны описывается таким выражени- ем:

 

1

 

 

 

 

 

g

 

 

 

 

 

 

nπ

 

 

y (x,t ) = −

 

 

 

 

 

 

n

 

 

 

 

 

sin

 

x exp(iωt ).

(3.104)

 

 

 

 

2

 

 

 

 

 

 

 

2

l

 

F n =1

(

ω

/c )

δ

ω

)

(n

π

 

 

 

 

 

 

 

 

i2

(

 

 

 

/l )

 

 

 

 

Здесь хорошо видно, что степень возбуждения каждой собственной формы колебаний зависит от степени близости заданного волнового числа ω/c и волнового числа собственной формы ωn/c = nπ/l. В случае резонанса, т.е. при ω/c = nπ/l, где n некоторый фиксированный номер, амплитуда колебаний струны остается конечной. Следует от- метить, что влияние демпфирования на вынужденные колебания струны будет отражаться на тех формах, для которых можно полагать

(ω/c )2 (nπ/l )2

2δ(ω).

(3.105)

Для тех собственных форм, у которых (ω/c )2 (nπ/l )2 >> 2δ(ω), на-

личие демпфирования несущественно. Следует также обратить вни-

138

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

Несколько сложнее учитывать демпфирование при свободных ко- лебаниях струны. Здесь необходимо найти решение однородного уравнения (3.99) (g(x) 0), которое удовлетворяет граничным и задан- ным начальным условиям:

y(0,t) = 0, y(l,t) = 0,

(3.106)

y(x,0) = Q1(x),

y (x,0)

= Q2 (x ).

 

t

 

 

 

 

Представим функции Q1(x) и Q2(x), которые определяют начальное возмущение, в виде рядов по собственным формам колебаний стру- ны:

Q1(x) =

 

nπ

x

 

,

qn

=

2 l

 

nπ

 

 

qn sin

 

l

 

 

l

Q1(x)sin

l

 

x dx;

 

n =1

 

 

 

 

 

 

 

 

0

 

 

 

 

 

nπ

 

 

 

 

 

 

 

2

 

l

 

nπ

 

Q2(x) =

βn sin

 

l

x

,

 

βn =

 

 

 

Q2(x)sin

 

 

x dx. (3.107)

 

 

ωnl

l

 

n =1

 

 

 

 

 

 

 

 

0

 

 

 

 

Если в этих рядах отличными от нуля являются два или более коэф- фициентов qn и βn, то это означает, что энергия начального возму- щения распределяется по нескольким собственным формам колеба- ний, и при дальнейшем движении струны будут присутствовать не- сколько нормальных колебаний. Каждое из них характеризуется соб- ственной частотой, и это нужно учитывать при рассмотрении демп- фирования, которое зависит от частоты.

Рассмотрим начальное возмущение, которое соответствует неко- торой конкретной форме sin(Nπx/l). В основу дальнейшего изложения положено предположение, что заданная форма колебаний будет оста- ваться неизменной; изменяется лишь амплитуда. Это означает, что решение однородного уравнения (3.99) для y(x,t) будем искать в виде

N

π

 

(3.108)

y (x,t ) =Y (t )sin

l

 

x .

 

 

 

 

После подстановки в выражение (3.99) при g(x) 0 решения (3.108) получим обычное дифференциальное уравнение для функции време- ни:

139

1 d2Y

+ 2δ(ω)

dY

N π 2

 

 

 

 

+

 

 

Y = 0.

(3.109)

c2

dt2

dt

l

 

 

 

 

 

Это достаточно простое дифференциальное уравнение, решение ко- торого нужно искать в виде

 

Y (t ) = AN exp(iαt ) .

 

(3.110)

Для показателя экспоненты получим равенство

 

 

N π 2

α2

+ i2αδ(ω) = 0.

 

(3.111)

 

c2

 

 

l

 

 

 

 

Решение квадратного

уравнения

α2 i2αδ(ω)c2 − ω2

= 0

несложно

найти в виде

 

 

 

N

 

 

 

 

 

 

 

 

α = iδ(ω)c

2 ±

 

ω2 − δ2 (ω)c4 .

 

(3.112)

 

 

 

 

N

 

 

Если демпфирование отсутствует, то будем иметь α = ±ωN, т.е. сис- тема осуществляет незатухающие колебания с собственной частотой ωN. В общем случае при отличном от нуля δ(ω ) колебание струны опи- сывается произведением экспонент:

2

 

±i

2

− δ

2

(ω)c

4

 

(3.113)

Y (t ) = AN exp(−δ(ω)c t )exp

ωN

 

 

t .

 

 

 

 

 

 

 

 

 

 

В случае демпфирования, меньше критического δ2 (ω)c4 < ω2N , имеем

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

допускают δ2 (ω)c4 << ω2N , и в (3.113) можно положить ω = ωN в каче-

стве первого приближения, т.е. δ(ω ) = δ(ωN). Физически это означает, что делается предположение о том, что демпфирование незначитель- но изменяет собственную частоту N-ой формы колебаний.

В случае необходимости полученный результат можно уточнить: вычисляем первое приближение для собственной частоты N-ой фор-

мы

 

 

 

 

 

ω(1)

=

ω2

− δ2 (ω

)c4 ,

(3.114)

N

 

N

N

 

 

а в окончательном результате полагаем ω = ω(N1) :

140