Эквивалентные параметры круглой однородной мембраны Текст научной статьи по специальности «Физика»

КУСТИКА

шашг

Электронный журнал «Техническая акустика» http://www .ejta.org

2014, 8

В. И. Кезик

Федеральный медицинский биофизический центр им. А. И. Бурназяна Россия, 123182, Москва, ул. Живописная, 46, e-mail: vladimirik57@mail.ru

Эквивалентные параметры круглой однородной мембраны

Кратко рассмотрена существующая методика определения эквивалентных параметров круглой однородной мембраны. Проанализирован недавно предложенный (в 2013 г.) новый подход к проблеме определения эквивалентных параметров круглой однородной мембраны. В рамках этого подхода параметры эквивалентного поршня определяются из условия тождественности законов движения при вынужденных колебаниях «усредненной» мембраны и эквивалентного поршня того же диаметра, помещенного на место мембраны, т.е. в тех же условиях возбуждения. Параметры эквивалентного поршня оказываются при этом зависящими от частоты возбуждения (кроме массы поршня). Введение дополнительного условия равенства средней за период кинетической энергии поршня и мембраны позволяет однозначно определить массу эквивалентного поршня, которая также становится частотно зависимым параметром. Аналогичные действия с потенциальной энергией поршня и мембраны позволяют определить параметр жесткости эквивалентного поршня другим способом. Сравнение этих результатов с помощью расчета показывает их полное совпадение. Произведено сопоставление нового метода со старым методом упрощенного расчета акустических систем, включающих мембраны. Расчетные формулы проиллюстрированы графиками, построенными для реальной мембраны. Приведены примеры расчета акустических колебательных систем, включающих круглую мембрану. Учтены поправки, связанные с излучением звука мембраной, совершающей колебания под действием вынуждающей силы.

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

1. ВВЕДЕНИЕ

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

Получена 22.06.2014, опубликована 30.10.2014

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

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

Для нахождения эквивалентных параметров использовался энергетический метод [1]. Кратко опишем его.

Отклонение точек мембраны от положения равновесия под действием статического равномерно распределенного по площади мембраны давления Р описывается (для малых отклонений) дифференциальным уравнением:

Т

' ё 2п 1 ёп л

ёг2 г ёг

+ -

= ~Р, (1.1)

где п = п(г) — поперечное отклонение кольцевого элемента мембраны радиуса г в

направлении оси п, Т — напряжение в сечении мембраны, Р — величина статического давления. Сила давления действует в направлении, совпадающем с направлением оси п . Параметр Т считаем константой, не зависящей от г .

Решение этого уравнения, удовлетворяющее граничным условиям (закрепление мембраны по круговому контуру), запишется в виде:

п (г) = щ

0 Р

Г1 - 4Л

Т у

(2.2)

где п0 Р = --отклонение центра мембраны, Е00 — радиус мембраны.

Пусть теперь Р меняется по гармоническому закону. Амплитудное значение отклонения центра мембраны обозначим через п0. Если предположить, что

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

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

мембраны и поршня обеспечивается при массе поршня тр, равной одной трети полной

массы мембраны тт:

ЕШ = 2 ] п2, (3.3)

1 2 1

тР = ^ЛЛ = 3тт . (4.4)

Потенциальная энергия мембраны может быть вычислена по работе, произведенной при ее растягивании [1]. Амплитудное значение потенциальной энергии мембраны:

^ = ^Ч2. (5.5)

Приравнивая потенциальную энергию мембраны и потенциальную энергию поршня 12

(пружины), равную kpu0, определим жесткость пружины осциллятора kp :

kp = 2жТ. (6.6)

Частота основного тона мембраны определяется резонансной частотой осциллятора:

®0 =

kp I 6T

тр V Рщ^2

(7.7)

Это выражение дает значение частоты, всего на 2 % отличающееся от точного значения частоты, определяемого последовательной теорией мембран:

2,405 Гг~

®о — ■ (8.8)

R0 \Рт

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

Позднее стали учитывать не только равенство энергий мембраны и замещающего ее эквивалентного поршня (кинетической и потенциальной), но и равенство объемов воздуха, вытесняемого при движении мембраны и замещающего ее поршня [2]. Это нужно для сохранения воздействия на окружающую среду при переходе от мембраны к замещающей системе. Другими словами, объем воздуха, вытесняемый при колебаниях поршня, не должен отличаться от объема, вытесняемого при колебаниях мембраны [2]. (Площадь поршня по-прежнему считаем равной площади мембраны.)

Как и ранее, считаем, что параболоидальная форма смещения (1.2) сохраняется при колебаниях на частоте основного тона мембраны. Вычисляя упомянутые выше объемы и приравнивая их, получим значение амплитуды колебаний эквивалентного поршня в этом случае [2]:

ир = ^ (9.9)

При этой амплитуде колебаний поршня {ир = и0 /2) равенство кинетической энергии

мембраны и поршня обеспечивается при массе поршня, равной 4/3 полной массы мембраны:

Еш=2 [4) ир, (10.10)

4 4

тР = 3 —^2 = 3 тт. (1111)

Потенциальная энергия мембраны осталась прежней (см. (1.5)). Приравниваем потенциальную энергию мембраны и потенциальную энергию поршня, равную 2кры'2р :

Еро( = лТп\ = 2 крп Р = 8 крп1. (12.12)

Отсюда находим жесткость пружины осциллятора в рассматриваемом случае: Кр = 8жТ . (13.13)

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

®0 =

к„

Р = ' "" 1 " (14.14)

тр ■■ ^

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

Мембраны имеют широкое распространение во многих отраслях техники и, являясь классическим объектом, продолжают привлекать внимание исследователей, как в качестве модельного объекта, так и в качестве конкретного технического устройства (см., например, [3, 4]). При этом мембраны, как правило, являются частью акустических колебательных систем. Разработка простого и эффективного метода расчета таких систем без ограничений по частоте колебаний представляется достаточно актуальной задачей.

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

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

2. НОВЫЙ подход к проблеме определения эквивалентных

ПАРАМЕТРОВ КРУГЛОЙ ОДНОРОДНОЙ МЕМБРАНЫ

2.1. Круглая мембрана в различных условиях возбуждения

В этом разделе дана постановка задачи, с некоторыми уточнениями по сравнению с [5].

Рассмотрим вынужденные колебания круглой закрепленной по контуру мембраны под действием приложенного к одной стороне мембраны равномерно распределенного по площади избыточного (звукового) давления, меняющегося по гармоническому закону с амплитудой р0. Другая сторона мембраны соприкасается с невозмущенной воздушной средой. (Такие условия возбуждения можно наглядно представить, если мембрану вставить в круглое отверстие тонкого бесконечно протяженного жесткого экрана, делящего пространство на две полусферы. Со стороны одной из полусфер (например, правой) на мембрану падает плоская звуковая волна, фронт которой совпадает с плоскостью экрана.) Учтем диссипативные потери (например, структурные потери), вводя в рассмотрение силы сопротивления движению элементов мембраны, которые, как обычно, будем считать пропорциональными скорости этих элементов. Однако, явления, связанные с излучением звука мембраной, совершающей вынужденные колебания в воздушной среде, пока рассматривать не будем, они будут рассмотрены в разделе 3. (Таким образом, на этом этапе мы рассматриваем колебания собственно мембраны, без учета влияния на нее прилегающей воздушной среды, то есть, как бы помещая мембрану в вакуум и воздействуя на нее некоторой силой, равномерно распределенной по площади мембраны и меняющейся по гармоническому закону.)

Такие условия возбуждения вызывают поперечные осесимметричные колебания рассматриваемой мембраны, которые описываются (для малых отклонений) дифференциальным уравнением в частных производных гиперболического типа [6]:

д 2и ды (д 2ы 1 ды л

о--У г--Т

Рт дл2 т дл

дг2 г дг

+ -

= Р0еш, (2.1)

где от - поверхностная плотность материала мембраны, и = и (г, I) - поперечное отклонение кольцевого элемента мембраны радиуса г в момент времени I (положительным считаем отклонение в сторону правой полусферы), гт - коэффициент сопротивления, отнесенный к единице площади, Т - напряжение в сечении мембраны, р0 - амплитуда поверхностной плотности вынуждающей силы, а> - круговая частота

вынуждающей силы, I = - мнимая единица. Параметры от, гт и Т считаем константами, не зависящими от г и I.

Граничное условие и условие ограниченности для функции и (г, I), справедливые при любых I, записывается в виде:

и(г, t) г=я„ = ^ I u(r, t(при г < Rо),

(2.2)

где R0 — радиус мембраны.

Для рассматриваемых нами установившихся колебаний начальные условия не требуются.

Решение уравнения (2.1) с граничными условиями (2.2) можем записать в виде (см. [5]):

и (г, I) = Р0

k 2Т

3№ )

Jо(kRо)

-1

(2.3)

где k = у—т—— (параметр k имеет смысл постоянной распространения

поперечной волны на поверхности мембраны), 30 — функции Бесселя первого рода

нулевого порядка от комплексного аргумента.

Это выражение полностью определяет установившееся движение мембраны в условиях возбуждения.

Пусть теперь на мембрану со стороны правой полусферы действует равномерно распределенное по площади мембраны избыточное давление р1, меняющееся по

гармоническому закону с амплитудой р01, а со стороны левой полусферы действует

также равномерно распределенное по площади мембраны избыточное давление р2 ,

меняющееся с той же частотой а>, но с амплитудой р02 и со сдвигом по фазе (12 по

отношению к давлению р1. Тогда правую часть уравнения (2.1) можем записать в виде

(р10 - р20вЩ2 )вш = р0вш, где р0 — комплексная амплитуда суммарного воздействия

давлений р1 и р2 на мембрану. Решение уравнения (2.1) теперь запишется в виде:

и (г, I) =

k 2т

3)

-1

(2.4)

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

В дальнейшем мы будем пользоваться более общим выражением (2.4).

2.2. Сводка необходимых формул

Далее приведем формулы, полученные в [5], которые понадобятся нам в дальнейшем. В этих формулах учтен более общий подход, даваемый выражением (2.4). Уравнение движения «усредненной» мембраны:

иш, «) = Р0

к 2Т

2 МкЮ к^ Jo(kRo)

-1

где J1 (кЯ0) — функция Бесселя первого рода первого порядка. Выражение (2.5) мы можем записать в виде итЫ(^) = итЫвш, где

(2.5)

и

к 2Т

2 МкЪ)

-1

(2.6)

— комплексная амплитуда колебаний «усредненной» мембраны. Уравнение движения поршня в тех же условиях возбуждения:

тр х + гРх + крх = р 0 ,

(2.7)

где тр — масса поршня, гр — коэффициент сопротивления движению поршня, кр —

жесткость пружины осциллятора, х — отклонение поршня от положения равновесия

(положительным считаем отклонение в сторону правой полусферы), 8р = жЯ02 —

площадь поршня, равная площади мембраны.

Решение уравнения (2.7) для установившегося режима находится методом комплексных амплитуд [7] и записывается в виде:

х(Г) =

~ о ¿а

р 0 V

кр - тра + ¡гра

Параметры эквивалентного поршня:

тР =От^р, ,

(2.8)

(2.9)

кр = тра + 8рТ Яе

2 МЯ) кЯ Jo(kRo)

-1

(2.10)

V т

гр =^^ Тт

а

2 МЯ)

кЯ Jo(kRo)

-1

(2.11)

где Яе и Тт — действительная и мнимая части комплексного числа, стоящего в фигурных скобках.

2.3. Сопоставление амплитудно-частотных характеристик эквивалентных поршней с параметрами, определяемыми по старой и новой методике

Для иллюстрации приведенных формул и сопоставления их с используемыми ранее

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

конкретную слабо натянутую мембрану, использованную в работе [5]. В данном случае

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

меняющегося односторонне приложенного внешнего давления. Параметры мембраны:

рт =0.063 кг/м2, гт =5 кг/м2с, Т=10 Н/м, Я =0.027 м. Амплитуда вынужденных

колебаний эквивалентного поршня определяется модулем комплексной амплитуды X из уравнения (2.8):

X = ■

Р 0

0" р

кр - трю + ггрф

(2.12)

На рис. 1 сплошной линией показана зависимость амплитуды вынужденных колебаний эквивалентного поршня от частоты возбуждения /, рассчитанная с помощью формулы (2.12) после подстановки в нее параметров эквивалентного поршня, определенных по формулам (2.9) - (2.11), или, что то же самое, с помощью формулы (2.6) (это справедливо только при рассмотрении колебаний собственно мембраны). Расчет сделан для рр 0 = р0 =20 Па (120 дБ).

На этом же графике пунктирной линией показана зависимость от частоты / (рассчитанная по той же формуле (2.12)) амплитуды колебаний поршня с параметрами, определенными по формулам (1.11) и (1.13) для той же мембраны в тех же условиях

4

возбуждения. Параметр гр в (2.12) определен по формуле гр = 3гт8р .

В

0.01

1x10

1x10

1x10

1x10

- 6

200

400

600

800 1x10

Рис. 1. Зависимость амплитуды колебаний эквивалентного поршня от частоты возбуждения, рассчитанная по новой (сплошная линия) и старой (пунктирная линия) методике. По оси абсцисс - частота в Гц, по оси ординат - амплитуда колебаний в

метрах

3

4

5

0

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

Рассмотрим комплексную амплитуду иш (см. формулу (2.6)). Если < 2,

функции Бесселя J0(kR0) и J1(kR0) представимы в виде степенных рядов [1]:

Л0Ц))=1 _ «1+сиу!,

о ^ 4 64

Jl(kRo)=

kRo

1 _ + .(kRo)4

192

Подставляем эти формулы в (2.6). При выполнении преобразований (деление, умножение многочленов), отбрасываем все члены со степенями k выше 4. В результате получим:

иmid

Ро ^

Т

1

- + -

8 192

^)2

Подставляя сюда k из формулы (2.3), получим:

Umid

РХ , 4 Ро Крт „2 , 4 РоКг.

а _I-

■а .

Х =

8Т 192 Т2 192 Т2

Далее рассмотрим формулу (2.12). Представим ее в виде:

Ро^о2

( т г ^

1 _ ^ а2 + а

V К К ;

Произведя аналогичные преобразования, получим:

Х =

Ро^ , Ро^ тР 2 • Ро^Дз2 ГР

1 а _i--—-а.

kp К К

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

к Р = 8жТ.

Приравнивая коэффициенты при одинаковых степенях а и учитывая полученное значение ,.определяем , и г,: „ = 4^ . г, = 3^.

Полученные формулы совпадают с теми, что получены в разделе 1 (введение) с помощью энергетического метода (формулы (1.11) и (1.13)).

2.4. Определение массы эквивалентного поршня

Как отмечалось в [5], в рамках метода эквивалентного поршня определение массы эквивалентного поршня по формуле (2.9) является избыточным. Рассчитываемые величины смещений (или давлений) и фаз колебаний любых элементов акустических колебательных систем рассматриваемого типа от этого параметра не зависят. Это является следствием кинематического подхода к рассматриваемой проблеме. Такой

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

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

Вернемся к осциллятору. Определим среднюю за период кинетическую энергию {ЕШр) поршня осциллятора, уравнение движения которого описывается формулой

(2.7). Однако, с учетом того, что мы постулировали тождественность законов движения поршня и «усредненной» мембраны, можем воспользоваться выражением (2.5) для описания движения поршня:

х(')=кТ

2 МкЯ) кЯ Л^)

-1

(2.13)

Дифференцированием по времени выражения (2.13) определим скорость поршня х ():

х «) =

1а рр 0

к 2Т

2 ЛТО кЯ Jо(кЯ,)

-1

(2.14)

которую запишем в виде: х(^) = Хвш, где

X =

тар 0

к 2Т

2 Jl(kR0) кЯ J0^)

-1

(2.15)

— комплексная амплитуда скорости движения поршня.

Вынужденные колебания поршня являются гармоническими, следовательно, средняя за период кинетическая энергия поршня равна половине от максимальной кинетической энергии, определяемой амплитудой скорости, которая, в свою очередь, определяется модулем комплексной амплитуды колебаний скорости XX (см. формулу (2.15)):

{ЕкТп, р^ 2

= -тр 2 р

ХТ

(2.16)

Аналогичным образом определим среднюю за период кинетическую энергию т ^

мембраны. Скорость кольцевого элемента радиуса г мембраны в каждый момент времени и(г, ^) определяется дифференцированием по времени выражения (2.4):

и(г,г) =

тр 0 к 2Т

3о(кг) 30(Щ,)

-1

(2.17)

Это выражение запишем в виде и (г, г) = и (г )вш, где

и (г ) =

тр 0

к 2Т

3о(кг)

3о(к^о)

-1

(2.18)

— комплексная амплитуда колебаний скорости кольцевого элемента мембраны радиуса г .

Каждый кольцевой элемент мембраны совершает вынужденные колебания со своей амплитудой и фазой. Скорость кольцевого элемента мембраны определяется амплитудой и фазой колебаний скорости, определяемыми комплексной амплитудой (2.18).

Определим среднюю за период кинетическую энергию (й.Ешкольцевого

элемента мембраны радиуса г и толщины йг, пользуясь теми же рассуждениями, что и при записи формулы (2.16):

(йЕШпт) = 2 йт

и (г )

(2.19)

2 2

Средняя за период кинетическая энергия всей мембраны определяется

интегралом:

(Ек1П,т) = р I \и(г)|2 гйг .

(2.20)

Массу эквивалентного поршня определим из условия равенства кинетической энергии мембраны и эквивалентного поршня, помещенного на место мембраны, то есть из равенства правых частей уравнений (2.16) и (2.20):

тр =

2Р

I - |2

\х\

I \й(г)|2 гйг .

(2.21)

Путем несложных преобразований с использованием тождества = |г2|, где и г2 — комплексные числа, выражение (2.21) можно привести к виду:

тр =

2прт\к 2| Д2

|2 31(к«0) - кД0 30(кДз)|'

||30(кг) -30(кД0)\2 гйг.

(2.22)

Таким образом, параметры поршня, эквивалентного мембране по кинематическим и энергетическим характеристикам, определяются выражениями (2.22), (2.10) и (2.11). Параметр тр в выражении (2.10) определяется выражением (2.22).

2

Д

В качестве иллюстрации на рис. 2 представлена зависимость массы эквивалентного поршня тр от частоты возбуждения / (Гц), рассчитанная по формуле (2.22) для той

же мембраны, которая использовалась при построении графиков на рис. 1.

Пунктирной линией на графике показано значение массы эквивалентного поршня, вычисленное по формуле (1.11) для этой же мембраны. Видим, что значения массы эквивалентного поршня, вычисленные по формулам (1.11) и (2.22), на низких частотах (менее 100 Гц) практически совпадают.

На рис. 3 представлена зависимость жесткости эквивалентного поршня кр от

частоты возбуждения /, рассчитанная по формуле (2.10) с учетом (2.22), для той же мембраны.

Пунктирной линией на графике показано значение жесткости эквивалентного поршня, вычисленное по формуле (1.13) для этой же мембраны. Видим, что значения жесткости эквивалентного поршня, вычисленные по формулам (1.13) и (2.10) с учетом (2.22), на низких частотах (менее 100 Гц) практически совпадают.

Вид приведенной на рис. 3 кривой существенно отличается от зависимости кр от /,

рассчитанной для случая, когда масса поршня приравнивалась массе мембраны (график этой зависимости приведен в работе [5]).

ад м

о, В

0.01

1x10

1x10

1x10

0 200 400 600 800 1x10

£ Н

Рис. 2. Зависимость массы эквивалентного поршня тр (сплошная линия) от частоты

возбуждения. По оси абсцисс — частота в Гц, по оси ординат — масса эквивалентного

поршня в килограммах

Л

1x10

1x10

1x10'

1x10

100

0 200 400 600 800 1x10

£ Н

Рис. 3. Зависимость жесткости эквивалентного поршня кр (сплошная линия) от частоты

возбуждения. По оси абсцисс — частота в Гц, по оси ординат — жесткость

эквивалентного поршня в Н/м

2.5. Потенциальная энергия мембраны

Параметры эквивалентного поршня определены. Однако для полноты картины и для проверки полученных ранее формул определим параметр кр другим способом, из

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

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

и (г, г) = Яе[и(г, г)] = Яе<

к 2Т

Л(кг)

-1

(2.23)

Потенциальная энергия, запасенная кольцевым элементом мембраны [2], может быть вычислена по формуле:

( ди Л 2

аЕРог,ш (г, г) = ^Т\ — I .

V дг

(2.24)

Потенциальная энергия всей мембраны в момент времени I определяется интегралом:

С ) = *Т]{§ I г*.

(2.25)

Функция Еро1 т (г) является периодической с частотой изменения, равной удвоенной

частоте возбуждения. Средняя за период потенциальная энергия всей мембраны определяется по формуле:

1 ТЕ

{Ерог, т ) = 1 Ерог,т (г)йг ,

Т

Е 0

1

(2.26)

где ТЕ = — период изменения потенциальной энергии мембраны в процессе 21

вынужденных колебаний, возбуждаемых с частотой / .

Определим теперь потенциальную энергию эквивалентного поршня. Вынужденные колебания эквивалентного поршня описывается формулой (2.13), которую можно представить в виде х(г) = Хв'т, где

X = -

к 2Т

2 ЛТО кД 30(кДз)

-1

(2.27)

— комплексная амплитуда колебаний эквивалентного поршня.

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

1Е ) = 1

рог, р р

(2.28)

Приравнивая правые части уравнений (2.26) и (2.28), определяем к :

ЛлТ'Е 7 (дй

К =

Т

//[§)

йг

(2.29)

IX2

Результаты расчета зависимости кр от частоты возбуждения /, сделанного по

формуле (2.29) для использованной нами ранее мембраны (см. комментарии к рис. 1) полностью совпали с результатами расчета, сделанного по формуле (2.10) (в которой тр определяется по формуле (2.22)) и представленного на рис. 3 (сплошная линия).

3. ИСПОЛЬЗОВАНИЕ ЭКВИВАЛЕНТНЫХ ПАРАМЕТРОВ КРУГЛОЙ

ОДНОРОДНОЙ МЕМБРАНЫ ДЛЯ РАСЧЕТА АКУСТИЧЕСКИХ

КОЛЕБАТЕЛЬНЫХ СИСТЕМ, СОДЕРЖАЩИХ МЕМБРАНЫ

Применение эквивалентных параметров круглой однородной мембраны для расчета акустических колебательных систем, содержащих мембраны подробно описано в работе [5], где этот метод расчета назван методом эквивалентного поршня. Здесь остановимся только на некоторых особенностях метода.

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

Если акустическая колебательная система такова, что мембрана окружена с обеих сторон полостями конечного объема V и V2, мы ограничены только тем, что должны

рассматривать вынужденные колебания элементов системы на частотах, при которых размеры этих полостей могли бы считаться малыми по сравнению с длиной звуковых волн, соответствующих этим частотам. В этом случае мембрана (или заменяющий ее эквивалентный поршень) может считаться нагруженной на упругий импеданс [8]

2 о2

ipc Sv

Z =--— . Здесь р - плотность воздуха, c - скорость звука, < - круговая частота

вынуждающей силы, - эффективный объем, определяемый из соотношения

= — + —, где V и V - объемы первой и второй полости. Эффекты, связанные с

Veff К К2

излучением звука мембраной или поршнем здесь отсутствуют. Ограничение по частотам не позволяет нам рассматривать вынужденные колебания системы на частотах, соответствующих высоким осесимметричным модам колебаний мембраны. Однако при малых объемах V и У2, малом натяжении мембраны и (или) большой поверхностной плотности материала мембраны могут быть рассмотрены относительно высокие осесимметричные моды.

Если в системе имеется открытая мембрана, одна из сторон которой соприкасается с неограниченным или, по крайней мере, обширным воздушным пространством, возникает необходимость учета присоединенной массы, величина которой зависит от конфигурации окружающего мембрану воздушного пространства, в которое происходит излучение. Нужно также учесть потери энергии на излучение, которые также зависит от конфигурации окружающего мембрану воздушного пространства, в которое происходит излучение. Эти явления описываются с помощью механического импеданса излучения [8], который будет рассмотрен ниже.

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

Теперь мы рассматриваем мембрану как акустическую колебательную систему, в полной мере взаимодействующую с прилегающей воздушной средой.

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

Активная и реактивная компоненты механического импеданса излучения круглого поршня в жестком экране (каждая сторона поршня излучает в открытое для нее полупространство) определяются с помощью функций Бесселя и Струве первого порядка [8]. Ряды, которыми могут быть представлены эти функции достаточно быстро

сходятся при 2—Я0 меньшем единицы (или, что то же самое, при X большем 4ЖЯ,, с

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

,,I (3.1)

жс

1брЯ3 (32)

Уш = —тш = —~г°. (3.2)

Здесь гг^ и уг^ — действительная и мнимая части механического импеданса излучения поршня, двусторонне излучающего в две полусферы, тгаа — присоединенная масса поршня.

Воспользуемся методом эквивалентного поршня [5], в соответствии с которым запишем уравнение движения эквивалентного поршня в рассматриваемых условиях:

(тр + тш ) * + (гр + Га* ) * + крх = р 0 8реш, (3.3)

где тр, гр, кр и 8р = жЯО — параметры эквивалентного поршня, определенные по

формулам (2.9)—(2.11), х — отклонение поршня от положения равновесия (положительным считаем отклонение в сторону правой полусферы).

Решение уравнения (3.3) для установившегося режима может быть найдено методом

комплексных амплитуд [7] и после несложных преобразований записано в виде:

—

р 0

*(*) =---

$рТк2

2МЯ) _ 1 кЯо Jo(kRo)

2 • (3.4)

_ тШ— + 1Ггас1 —

X =

Запишем эту формулу в виде х(V) = Хеш, где

Р о ЯР

Я/к2

2А(кВо) _ 1

кЯо 3о(к^о)

_ тга® + 1Гга®

(3.5)

— комплексная амплитуда колебаний эквивалентного поршня.

Амплитуда колебаний эквивалентного поршня определяется модулем комплексной амплитуды X.

На рис. 4 представлена зависимость амплитуды колебаний эквивалентного поршня |Х| от частоты возбуждения /, рассчитанная по формуле (3.5) для рассмотренной выше (см. комментарии к рис. 1) мембраны в тех же условиях возбуждения.

0.01

1x10

г^ 1x10

_ 3

В

1x10

1x10

_ 6

200

400

Г, н

600

800

Рис. 4. Зависимость амплитуды колебаний эквивалентного поршня без учета (пунктирная линия) и с учетом поправок на излучение (сплошная линия) от частоты возбуждения. По оси абсцисс — частота в Гц, по оси ординат — амплитуда колебаний

эквивалентного поршня в метрах

На рис. 4 мы видим, что учет поправок на излучение существенно сдвигает частоту первого резонанса мембраны в сторону уменьшения. Примерно такой же результат дает расчет радиационных поправок для значительно более сильно натянутых мембран. Анализ формулы (3.5) показывает, что определяющее влияние на изменение вида амплитудно-частотной характеристики поршня (поправки на излучение) оказывает величина присоединенной массы. Влияние активной компоненты импеданса излучения несущественно.

В заключение кратко рассмотрим поправки на излучение для задачи о вынужденных колебаниях мембраны барабана [5]. Будем считать в этом случае, что эквивалентный поршень, заменяющий мембрану, односторонне излучает в открытое пространство (без экрана). В рамках тех же приближений, что и в предыдущей задаче, активная и реактивная компоненты механического импеданса излучения могут быть записаны в виде [8]:

4

5

0

Г = ■

гad

Р—% 4жс

(3.6)

УгаЛ = тга— = 2РК— .

(3.7)

Здесь ггаЛ и угаЛ — действительная и мнимая части механического импеданса одностороннего излучения поршня в сферу, тгаЛ — присоединенная масса поршня. Уравнение движения эквивалентного поршня для данной задачи запишется в виде:

(

(тр + тгаЛ ) * + (Гр + ГгаЛ ) * +

к, + -2 * ^

V

V

* = Ро 8ре

(3.8)

у

где тр, гр, кр и 8р = ж Я — параметры эквивалентного поршня, определенные по

формулам (2.9)—(2.11), V — объем полости барабана.

Решение уравнения (3.8), как обычно, находим методом комплексных амплитуд [7] и после несложных преобразований комплексную амплитуду колебаний эквивалентного поршня записываем в виде:

Р о _

" " . (3 9)

X = -

ЪрТк2

2 Jl(kRo)

кЯо Jо(кЯо)

_ 1

рс2 ^

_V

2?2 р2

_ тгаЛ— + 1ГгаЛ —

Формула (3.9) определяет амплитуду и фазу вынужденных колебаний мембраны

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

4. ЗАКЛЮЧЕНИЕ

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

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

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

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

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

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

показывает их полное совпадение.

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

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

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

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

ЛИТЕРАТУРА

1. Крендалл И. Б. Акустика. М.: Книжный дом «ЛИБРОКОМ», 2009. 168 с.

2. Вахитов Я. Ш. Теоретические основы электроакустики и электроакустическая аппаратура. М.: Искусство, 1982. 415 с.

3. Канев Н. Г., Миронов М. А. Дипольный резонансный рассеиватель звука // Акуст. журн. 2003. Т. 49, № 3. С. 372-375.

4. Смирнов И. П., Бурдуковская В. Г., Кошкин А. Г., Хилько А. И. Нелинейные колебания кольцевых мембран низкочастотного акустического излучателя // Изв. вузов. Радиофизика. 2008. Т. 51, № 3. С. 199-215.

5. Кезик В. И. Вынужденные осесимметричные колебания круглой мембраны, являющейся элементом акустической колебательной системы // Электронный журнал «Техническая акустика», http://ejta.org, 2013, 9.

6. Тихонов А. Н., Самарский А. А. Уравнения математической физики. М.: Наука, 1972. 735 с.

7. Стрелков С. П. Введение в теорию колебаний. М.: Наука, 1964. 440 с.

8. Ржевкин С. Н. Курс лекций по теории звука. М.: Изд-во Московского университета, 1960. 336 с.