Министерство науки и высшего образования Российской Федерации
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ
УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ
НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ ИТМО
ITMO University
ВЫПУСКНАЯ КВАЛИФИКАЦИОННАЯ РАБОТА/GRADUATION THESIS
Моделирование барабана
Автор/ Author
Шипицын Иван Антонович
Направленность (профиль) образовательной программы/Major
Математическое моделирование 2019
Квалификация/ Degree level
Магистр
Руководитель ВКР/ Thesis supervisor
Попов Игорь Юрьевич, профессор, доктор физико-математических наук, Университет
ИТМО, факультет систем управления и робототехники, профессор (квалификационная
категория "ведущий профессор")
Группа/Group
R42952
Факультет/институт/кластер/ Faculty/Institute/Cluster
факультет систем управления и робототехники
Направление подготовки/ Subject area
01.04.02 Прикладная математика и информатика
Обучающийся/Student
Документ
подписан
Шипицын Иван
Антонович
01.06.2021
(эл. подпись/ signature)
Руководитель ВКР/
Thesis supervisor
Документ
подписан
Попов Игорь
Юрьевич
01.06.2021
(эл. подпись/ signature)
Шипицын Иван
Антонович
(Фамилия И.О./ name
and surname)
Попов Игорь
Юрьевич
(Фамилия И.О./ name
and surname)
Министерство науки и высшего образования Российской Федерации
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ
УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ
НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ ИТМО
ITMO University
ЗАДАНИЕ НА ВЫПУСКНУЮ КВАЛИФИКАЦИОННУЮ РАБОТУ /
OBJECTIVES FOR A GRADUATION THESIS
Обучающийся / Student Шипицын Иван Антонович
Группа/Group R42952
Факультет/институт/кластер/ Faculty/Institute/Cluster факультет систем управления и
робототехники
Квалификация/ Degree level Магистр
Направление подготовки/ Subject area 01.04.02 Прикладная математика и информатика
Направленность (профиль) образовательной программы/Major Математическое
моделирование 2019
Специализация/ Specialization Математическое моделирование сложных систем
Тема ВКР/ Thesis topic Моделирование барабана
Руководитель ВКР/ Thesis supervisor Попов Игорь Юрьевич, профессор, доктор физикоматематических наук, Университет ИТМО, факультет систем управления и робототехники,
профессор (квалификационная категория "ведущий профессор")
Срок сдачи студентом законченной работы до / Deadline for submission of complete
thesis 20.05.2021
Техническое задание и исходные данные к работе/ Requirements and premise for the
thesis
Смоделировать динамику мембраны барабана в зависимости от характера удара
барабанной палочкой
Содержание выпускной квалификационной работы (перечень подлежащих
разработке вопросов)/ Content of the thesis (list of key issues)
1) Построение модели взаимодействия барабанной палочки и мембраны барабана
2) Численное моделирование динамики мембраны барабана
Перечень графического материала (с указанием обязательного материала) / List of
graphic materials (with a list of required material)
Исходные материалы и пособия / Source materials and publications
1) Bilbao, Stefan. (2012). Time domain simulation and sound synthesis for the snare drum. The
Journal of the Acoustical Society of America. 131. 914-25. 10.1121/1.3651240.
2) Jin, Congrui & Davoodabadi, Ali & Li, Jianlin & Wang, Yanli & Singler, Timothy. (2017).
Spherical indentation of a freestanding circular membrane revisited: Analytical solutions and
experiments. Journal of the Mechanics and Physics of Solids. 100. 10.1016/j.jmps.2017.01.005.
Дата выдачи задания/ Objectives issued on 05.05.2021
СОГЛАСОВАНО / AGREED:
Руководитель ВКР/
Thesis supervisor
Документ
подписан
Попов Игорь
Юрьевич
05.05.2021
Попов Игорь
Юрьевич
(эл. подпись)
Задание принял к
исполнению/ Objectives
assumed by
Документ
подписан
Шипицын Иван
Антонович
06.05.2021
Шипицын Иван
Антонович
(эл. подпись)
Руководитель ОП/ Head
of educational program
Документ
подписан
Попов Игорь
Юрьевич
02.06.2021
(эл. подпись)
Попов Игорь
Юрьевич
Министерство науки и высшего образования Российской Федерации
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ
УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ
НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ ИТМО
ITMO University
АННОТАЦИЯ
ВЫПУСКНОЙ КВАЛИФИКАЦИОННОЙ РАБОТЫ /
SUMMARY OF A GRADUATION THESIS
Обучающийся/ Student
Шипицын Иван Антонович
Наименование темы ВКР / Title of the thesis
Моделирование барабана
Наименование организации, где выполнена ВКР/ Name of organization
Университет ИТМО
ХАРАКТЕРИСТИКА ВЫПУСКНОЙ КВАЛИФИКАЦИОННОЙ РАБОТЫ/
DESCRIPTION OF THE GRADUATION THESIS
1. Цель исследования / Research objective
Моделирование барабана
2. Задачи, решаемые в ВКР / Research tasks
Разработка модели взаимодействия барабанной мембраны и барабанной палочки.
Численное моделирование динамики барабанной мембраны при взаимодействии с
барабанной палочкой.
3. Краткая характеристика полученных результатов / Short summary of
results/conclusions
Предложена модель модель взаимодействия барабанной мембраны и барабанной палочки.
Численно смоделирована динамика мембраны барабана при взаимодействии с барабанной
палочкой, проведены сравнительные характеристики.
4. Наличие публикаций по теме выпускной работы/ Have you produced any
publications on the topic of the thesis
5. Наличие выступлений на конференциях по теме выпускной работы/ Have you
produced any conference reports on the topic of the thesis
6. Полученные гранты, при выполнении работы/ Grants received while working on the
thesis
7. Дополнительные сведения/ Additional information
Обучающийся/Student
Документ
подписан
Шипицын Иван
Антонович
01.06.2021
(эл. подпись/ signature)
Руководитель ВКР/
Thesis supervisor
Документ
подписан
Попов Игорь
Юрьевич
01.06.2021
(эл. подпись/ signature)
Шипицын Иван
Антонович
(Фамилия И.О./ name
and surname)
Попов Игорь
Юрьевич
(Фамилия И.О./ name
and surname)
Оглавление
1.
Введение .......................................................................................................................................8
2.
Строение барабана и особенности игры.........................................................................11
Барабанная мембрана ........................................................................................................11
2.1.
2.1.1.
Идеальная круглая мембрана ....................................................................................12
2.2.
Барабанный кожух ...............................................................................................................13
2.3.
Барабанная палочка ............................................................................................................16
3.
Анализ экспериментальных данных ................................................................................19
4.
Модель взаимодействия .......................................................................................................22
5.
Колебания мембраны .............................................................................................................25
6.
Свободные колебания идеальной круглой мембраны ..............................................30
7.
Вынужденные колебания идеальной круглой мембраны ........................................40
7.1.
Параметры мембраны ........................................................................................................51
7.2.
Параметры палочки.............................................................................................................51
7.3.
Время взаимодействия ......................................................................................................52
8.
Результаты численного моделирования ........................................................................53
9.
Заключение ................................................................................................................................80
10.
Список литературы .............................................................................................................82
8
1.
Введение
Ударные инструменты являются одними из древнейших музыкальных
инструментов
в
мире.
Они
также
являются
одними
из
самых
распространенных: Музыкальный словарь Гроува [1] включает себя более
1500 вхождений для барабанов. Кроме того, все новые ударные инструменты
продолжают создаваться и разрабатываться и по сей день, как акустические,
так и электронные.
Барабаны – технически классифицируемые в семействе перкуссионных
инструментов как мембранофон – играли большую роль в развитии почти
каждой музыкальной культуры. Самые ранние барабаны представляли собой
незамысловатые конструкции из камня или дерево, водруженные над дырами
в земле. Сегодня же, самым привычным для нас барабаном является натянутые
на некоторый ограниченный объем с воздухом кожа животного или
синтетический материал. Предполагается, что такой современный вид
барабана появился примерно 6000 лет до нашей эры.
Задача моделирования барабанов вызывает живой интерес как с
практической, так и с фундаментальной стороны. Многие авторы занимались
изучением характера поведения мембраны барабана, а также моделированием
полноценных барабанов. Рэйли в 1877 году установил [2] природу
фундаментальных нот литавра, основываясь на модели идеальной круглой
мембраны. Флетчером в третьей четверти двадцатого века активно изучались
бас-барабаны [3]. В 1984 году Кристиан, Дэвис и Андерсон проанализировали
[4] влияние нагрузки со стороны воздуха на мембрану в литаврах. На
протяжении последней четверти прошлого столетия Россинг произвел [5-8]
анализ и моделирование многочисленных перкуссионных инструментов со
всего мира. Им были получены и описаны множественные краткие
9
характеристики различных инструментов. Также в 1999 году Раути была
предложена [9] и реализована модель литавров. В частности, им были
смоделированы малые барабаны с использованием модели идеально-круглой
мембраны. В 2013 году Бильбао синтезировал [10] звук малого барабана,
используя 3D-модель инструмента. Для описания мембраны была выбрана
модель идеальной круглой мембраны, возбуждаемая точечной внешней силой.
Несмотря на обширное количество исследований и предложенных
моделей, малое внимание в них уделялось взаимодействию барабана с
объектом, возбуждающи колебания. По аналогии с моделями, предложенными
для изучения взаимодействия молоточков со струнами фортепиано [12],
некоторыми исследователями [11-14] была предложена модель груз-пружина
для описания взаимодействия палочки и мембраны. В Джин [15, 16] и Бэйгли
[17] изучали поведение мембраны при контакте со сферическим и
цилиндрическим грузами.
Однако зачастую в предложенных моделях удар описывается как
действие точечной внешней силы. И хоть такая модель очень удобна и
значительно упрощает решение задачи, а также процесс вычислений, она
предназначена лишь для одиночных сильных и быстрых ударов, а также
опускает такие параметры, как, например, размер наконечника палочки.
Ударные инструменты же позволяют получать разнообразные тембры при
применении составных жестов, полноценные модели для описания которых
еще не были предложены.
В настоящей работе на основе существующих экспериментальных
данных предлагаются две модели для описания взаимодействия барабанной
палочки с мембраной барабана. На основе предложенных моделей проводится
численной моделирование динамики барабанной мембраны и анализ
изменений в динамике в зависимости от таких параметров, как тип, сила, место
удара.
10
Результаты настоящей работы были апробированы на X Конференции
Молодых Ученых. На конференции был представлен доклад, а также были
выпущены тезисы.
11
2.
Строение барабана и особенности игры
2.1. Барабанная мембрана
Основным источником возбуждения звуковых волн в ударных
инструментах является мембрана. По этой причине, вариации различных
физических параметров барабанной мембраны приводит к получению самого
разнообразного спектра звучаний. Так, многие африканские барабаны
(например, джембе) используют козлиную кожу в качестве мембраны,
натянутую на деревянное основание с помощью веревки в специальных
узорах. Такие барабаны в зависимости от предназначения, могут быть
орнаментированы различными металлическими кольцами, пластинами или
погремушками по краям мембраны.
Рис. 1. Джембе
Большинство же барабанов в современной западноевропейской
культуре в качестве мембраны, начиная с последней трети двадцатого века
[24] используется Майлар – специальная пленка, изготавливаемая на основе
полиэтилентерефталата. Данный материал в отличии от его натуральных
аналогов является однородным по своей структуре, а также у него отсутствует
чувствительность к влажности, что позволяет с большей простотой
настраивать звучание барабанных мембран из него изготовленных. Кроме
12
того, данный материал является очень тонким: так, стандарт толщины для
майларовых барабанных мембран – всего 0.19 мм.
В физико-математических моделях, как для музыкальных инструментов,
так и для акустических систем, для представления мембран наиболее
распространена модель тонкой идеальной круглой мембраны.
2.1.1.
Идеальная
Идеальная круглая мембрана
круглая
мембрана
–
модель
мембраны,
которая
характеризуется тремя параметрами: диаметром, натяжением и плотностью
материала мембраны. Модель идеальной круглой пренебрегает такими
явлениями, как неоднородность материала мембраны, сопротивление
растяжениям. Зачастую предполагается, что идеальная круглая мембрана
является жестко закрепленной по краям.
Идеальная
круглая
мембрана,
закрепленная
по
краям,
имеет
характерные колебания, описывающиеся суперпозицей колебательных мод.
Такие моды описываются парой значений ( n, m) , где число n задает
количество нодальных диаметров, а m - нодальных окружностей. На таких
нодальных линиях смещение мембраны является минимальным. На рисунке 2
представлены первые 12 колебательных мод идеальной круглой мембраны.
Известно, что частоты нормальных колебаний идеальной круглой
мембраны задаются формулой
f mn
ckmn
j
T
mn
2
2 a
(1)
где jmn - m-й ноль функции Бесселя первого рода J n ( kr ) , k – волновое число,
T – натяжение мембраны, а - поверхностная плотность мембраны. В силу
того, что функции Бесселя не являются периодическими функциями, т.е.
13
значения ее нулей не имеют некоторую фундаментальную частоту в качестве
множителя, модальные частоты идеально круглой мембраны не являются
гармонически созвучными.
В настоящей работе для моделирования мембраны барабана будет
использоваться модель идеальной круглой мембраны.
Рис. 2. Колебательные моды идеальной круглой мембраны
2.2. Барабанный кожух
Барабанный кожух является неотъемлемой частью практически всех
существующих барабанов. Размер, форма, материал и наполнение кожуха
могут радикально менять характер результирующего тембра барабана.
14
В западноевропейской культуре одной из самых распространенных
форм барабанного кожуха является цилиндр, обтянутый с одной стороны
мембраной. В таких барабанах, как том-том или бас-барабан, такой цилиндр
ограничен с двух сторон барабанными мембранами. Малые или рабочие
барабаны имеют натянутые по низу поверх второй мембраны специальные
проволоки и пружинки, придающие этому инструменту особый яркий
шипящий звук.
Подобные конфигурации производят усложненные спектральные
формы излучаемого звука в силу того, что две барабанные мембраны могут
колебаться как в одной фазе друг с другом, так и в противофазах. Подобное
поведение может результировать в «расщеплении» фундаментальных тонов
[7].
Рис. 3. Малый барабан Ludwig
Другой распространенной формой является котлообразная чаша – такую
форму принимают литавры. Несмотря на негармоничность мод колебаний
идеальной мембраны, аккуратно настроенные литавры позволяют извлекать
ярковыраженные фундаментальные ноты. В частности, известно [2], что
фундаментальный тон литавра исходит от моды (1, 1), а также обертона
порядка чистой пятой, большой седьмой и октавы над основным тоном
исходят от мод (2, 1), (3, 1) и (1, 2) соответственно. Более того, отношения
частот между модами (1, 1), (3, 1), (1, 2) в литаврах имеет вид 1 к 1.5 к 2, а
15
моды (4, 1) и (5, 1) – отношения 2.44 и 2.90 к фундаментальной моде (1, 1)
соответственно.
Рис. 4. Литавры
Кожух, как и мембрана барабана, имеет собственные моды колебания.
Несмотря на то, что по величине эти колебания малы по сравнению со
смещениями колеблющийся мембраны, сильная связь между кожухом и
мембраной может привести к значительному вкладу мод кожуха в
результирующее поведение мембраны.
Рис. 5. Рототом
16
В настоящей работе вкладом кожуха на динамику мембраны
пренебрегается, как для упрощения вычислений, так и в силу того, что
характеристики кожуха не влияют на характер взаимодействия ударным
инструментом с мембраной барабана. Такая модель может соответствовать
рототому – барабану, являющемуся мембраной на металлическом обруче.
2.3. Барабанная палочка
Для воспроизведения звука, по поверхности барабана производится
удар. Существуют различные семейства барабанов, игра на которых
подразумевает использование рук – джембе, дарбука, индийские таблы. На
некоторых баранах, таких как литавры, исполнитель чаще играет с
использованием специальных колотушек или маллетов, сделанных из дерева
или тростника со специальным наконечником, обычно деревянным, зачастую
обтянутым войлоком, пробкой или кожей. В таких стилях музыки, как джаз и
блюз, широко распространена игра специальными барабанными щеточками.
Однако в современном мире наиболее распространена игра барабанными
палочками.
Типичная барабанная палочка сделана из дерева - обычно дуба, гикори
(светлого ореха) или клена, реже могут использоваться метал, алюминий и
углеродное волокно. Состоит барабанная палочка из четырех частей:
комеля – закругленного основания со слегка увеличенным весом,
тела,
плеча – сужающейся области палочки, и
наконечника.
Наконечник может быть изготовлен как из того же материала, что и сама
палочка, так и из нейлона или ацеталя.
17
Наконечники барабанных палочек изготавливаются различных форм и
размеров. По форме, наконечники можно классифицировать по трем
основным группам:
Цилиндрические
Круглые
Овальные
Остроконечные
Различные формы наконечников позволяют исполнителю менять характер
воспроизводимого инструментом звука. Так, овальные головки позволяют
менять звучание изменяя прямо во время игры площадь поверхности
соприкосновения палочки с поверхностью барабанной мембраны за счет
регулирования угла наклона палочки к мембране. Различные формы также
помогают
исполнителю
лучше
контролировать
отскок
поверхности инструмента.
Рис. 6. Различные виды наконечников барабанных палочек
палочки
от
18
Материал изготовления, линейные размеры и вес барабанной палочки
оказывают непосредственное влияние на два параметра: сила удара и время
контакта с мембраной барабана. К сожалению, в настоящий момент не было
произведено выделенных экспериментальных исследований по зависимости
этих или иных параметров от характеристик используемой барабанной
палочки.
В настоящей работе в целях упрощения вычислений предполагается, что
удар по мембране барабана будет наносится барабанной палочкой с
наконечником идеальной круглой формы. Такие параметры палочки, как
материал палочки, материал наконечника, вес, размер и длинна не будут
непосредственно вкладываться в величину силы удара или взаимодействия
палочки с мембраной барабана.
19
3.
Анализ экспериментальных данных
Для построения модели, описывающей процесс взаимодействия
барабанной палочки с мембраной барабана, обратимся к экспериментальным
данным.
В эксперименте по описанию взаимодействия барабанщика, барабанной
палочки и мембраны барабана [18] производились измерения силы удара и
времени взаимодействия палочки с мембраной барабана при исполнении
одиночного удара.
Рис. 7. Зависимость силы давления палочки на поверхность мембраны от времени контакта в центре
мембраны с низким натяжением [18].
Из результатов, представленных на рисунке 7, можно видеть, что сила,
оказываемая барабанной палочкой, имеет промежуток нарастания (0-1 мс на
графике), величиной порядка 0.75 мс, а также промежуток спада (2-4 мс на
графике), величиной также порядка 0.75 мс. Величина силы давления палочки
на мембраны принимает максимальные значения на протяжении примерно 3.5
мс.
Также в рамках эксперимента удары были исполнены в трех
положениях, отличных от центра:
20
Рис. 8. Зависимость силы давления палочки на поверхность мембраны от времени контакта для
четырех различных областей при низком натяжении.
Из результатов, представленных на рисунке 4, можно видеть, что все три
временных параметра (время нарастания, время максимальной величины,
время спада) уменьшаются с тем, как область удара удаляется от центра
мембраны. Это связанно с тем, что в настоящих барабанных мембранах
величина натяжения меняется в большую сторону при приближении к краю
мембраны. Этой же особенностью в исследовании объясняется увеличение
силы удара в амплитуде при отдалении от центра.
Рис. 9. Сравнение времени контакта палочки с мембраной при разных величинах натяжения для
четырех различных областей.
Изучая зависимость времени взаимодействия от натяжения мембраны,
исследование обнаружило, что с увеличением натяжения мембраны,
продолжительность контакта с палочкой уменьшалась. При этом на рисунке 6
21
можно видеть, что время взаимодействия палочки с мембраной уменьшается
линейно по мере отдаления от центра мембраны.
Исследование также включало в себя использование различных
барабанных палочек. В результате было обнаружено, что использование
палочки того же материала с большим диаметром наконечника той же формы
увеличивает
амплитуду
силы
удара,
и
в
следствии
амплитуду
результирующего звука. Изменений во в
Даль [19] в своем исследовании характера звука том-тома в зависимости
от силы удара, а также амплитуды палочки, установила зависимость времени
контакта палочки с мембраной барабана от силы наносимого удара
(Результаты представлены на рисунке 10).
Рис. 10. Зависимость времени контакта палочки и мембраны от силы удара (p – piano – слабый удар,
f – forte – сильный удар).
Оба исследования заключили, что при увеличении силы удара,
наблюдаются следующие явления:
Увеличение амплитуды результирующего звука барабана
Усиление интенсивности колебательных мод, возбуждаемых при
более слабом ударе в ту же область
Уменьшение времени взаимодействия палочки и мембраны
22
4.
Модель взаимодействия
Опишем теперь модель, представляющую взаимодействие мембраны и
барабанной палочки.
Пусть идеальная круглая мембрана вступает в контакт с твердым телом
в форме полусферы, оказывающим на нее некоторое давление P0 , в
произвольной точке мембраны. Будем считать, что область контакта
мембраны и твердого тела представляет собой окружность радиуса c такую,
что c R .
Выберем теперь оси координат так, чтобы ось абсцисс лежала на
прямой, соединяющий точки центра окружностей мембраны и области
соприкосновения i , и положим центр координат O в точке центра
мембраны.
Рис. 11. Схематическое изображение идеально-круглой мембраны в контакте со сферическим
твердым телом
23
В соответствии с результатами, описанными в предыдущей главе,
взаимодействие палочки и мембраны можно представить следующими двумя
моделями.
Первая модель представляет взаимодействие палочки и мембраны как
ступенчатую функцию.
В начальный момент времени – момент контакта палочки с мембраной,
– некоторой области мембраны передается некоторый импульс P0 . На
протяжении некоторого времени T1 (порядка 1 мс), мембрана совершает
свободные колебания. По прошествии времени T1 , на мембрану в области
контакта i начинает действовать постоянная внешняя сила
F (r , ) P0
(2)
и продолжает непрерывно действовать до окончания контакта палочки и
мембраны. Будем называть такую модель «ступенчатой».
Такая модель подходит для описания ударов, в которых исполнитель
оказывает дополнительное давление на палочку, тем самым продлевая время
контакта палочки с мембраной барабана. Такая модель требует меньше
вычислений в силу отсутствия зависимости величины внешней силы от
времени.
Вторая модель представляет взаимодействие палочки и мембраны с
помощью тригонометрической функции от времени. В начальный момент
времени – момент контакта палочки с мембраной, – некоторой области
мембраны передается некоторый импульс P0 . В этот же момент на мембрану
начинает действовать внешняя сила, выразить которую можно как
F ( r , , t ) P0 sin( qt )
где
(3)
24
q
(4)
T
где T – время взаимодействия палочки и мембраны. Выбор величины T
производится на основе величины натяжения мембраны, амплитуды и силы
удара, характеристик палочки, таких как вес, материал наконечника, а также
отдаления от центра. Для каждого набора таких параметров, данная величина
может быть экспериментально определена. Например [18], для деревянной
палочки весом 53 г, с деревянным наконечником диаметра 8 мм, при
натяжении мембраны 2500 Н/м, время взаимодействия барабанной палочки и
мембраны при ударе в центральную область будет равно 5 мс, в то время как
при ударе в крайнюю область – 4 мс. Для удобства же будем считать, что при
рассмотрении
идеальной
мембраны,
время
взаимодействия
между
центральной и краевой областью одинаково. Будем называть такую модель
«тригонометрической».
Преимуществом
такой
модели
является
более
равномерное
распределение воздействия внешней силы по времени, приближенного к
соответствующим экспериментальным данным. В то же время такая модель
непригодна для моделирования более сложных жестов, во время которых
исполнитель может увеличивать продолжительность контакта палочки с
мембраной. Кроме того, данный подход увеличивает количество вычислений.
25
5.
Колебания мембраны
Мембраной в математической физике называется не оказывающая
сопротивления сдвигам и изгибам тонкая, плоская пленка.
Пусть мембрана натянута на некоторый плоский контур Г и совершает
колебания. Рассмотрим такие колебания мембраны, смещение при которых
будет перпендикулярным к плоскости этой мембраны – поперечные
колебания.
Положим, что ds - малый элемент на дуге выбранного контура,
проходящий через некоторую точку M ( x , y ) выбранной на поверхности
рассматриваемой мембраны. На такой элемент действует натяжение, положим
Tds . При этом вектор натяжения T будет лежать, в результате отсутствия
сопротивления сдвигам и сгибам, в касательной плоскости к мгновенной
поверхности мембраны, а также будет перпендикулярным к выбранному
элементу. В силу тех же причин, величина вектора натяжения T не зависит
от направления элемента поверхности ds , что означает, что вектор
T T ( x , y , z , t ) будет являться функцией лишь от x , y , t
Рассмотрим малые колебания такой мембраны с поверхностной
плотностью ( x , y ) . Пусть функция u u ( x, y , t ) определяет отклонение
точек мембраны от плоскости ( x , y ) в момент времени t . Рассмотрим нормаль
N к поверхности мембраны с набором проекций, образующих острый угол с
поверхностью u :
u u
, ,1
x y
Направляющие косинусы для такой нормали задаются формулами
26
cos
ux
1 u x2 u y2
cos
,
cos
uy
1 u x2 u y2
1
1 u x2 u y2
Рассмотрим проекцию Th ( x , y , z , t ) вектора натяжения на плоскость
( x , y ) . Заметим, что при произвольной ориентации элементарной дуги ds угол
' , образованный между вектором натяжения T и плоскостью ( x , y ) не будет
превосходить по величине угол между осью z и нормалью к поверхности
мембраны в выбранной точке. Отсюда следует, что
cos cos
1
1 u x2 u y2
(10)
В силу малости рассматриваемых колебаний будем пренебрегать квадратами
производных первого порядка функции u :
2
2
u
y 1
u
1,
x
(11)
Тогда
1
1
(12)
cos 1 ,
(13)
1 u u
2
x
2
y
Но тогда
и выражение для проекции вектора натяжения
Th ( x, y, z , t ) T cos ' T ( x, y , z , t ) ,
(14)
27
т.е. величина проекции вектора натяжения в любой момент времени равна его
абсолютной величине. Вертикальная составляющая вектора натяжения при
этом будет равна
Tu T
u
.
n
(15)
Выберем теперь на поверхности рассматриваемой мембраны такой
элемент площади, что его проекция на плоскость ( x , y ) будет представлять
собой прямоугольник, стороны которого параллельны осям координат.
Рис. 12. Проекция на плоскость ( x, y )
На такой элемент поверхности будет действовать сила натяжения
T*
Tds
(16)
ABCD
Так как проекция T * вдоль осей x , y не имеет перемещения, и потому для
соответствующих интегралов получим
C
D
Tx* T x2 , y , t dy T x1 , y , t dy
B
A
y2
y1
T x , y, t T x , y, t dy 0
2
1
(17)
28
B
D
Ty* T x, y 2 , t dy T x, y1 , t dy
A
C
x2
x1
T x, y , t T x, y , t dy 0
2
(18)
1
Тогда применяя теорему о среднем вкупе с фактом, что прямоугольник ABCD
является произвольным
T x, y1 , t T x, y2 , t
T x1 , y , t T x2 , y , t
(19)
Таким образом мы получили, что величина натяжения T не зависит от ( x , y ) .
Вычислим теперь площадь малого элемента мембрана в выбранный
момент времени:
dxdy
cos
1 u x2 u y2 dxdy d xdy
(20)
Полученный результата означает, что в процессе колебаний мембраны не
происходит растяжений, т.е. изменениями площади поверхности мембраны
можно пренебречь. Но тогда по закону Гуку можно утверждать, что величина
натяжения T не зависит от времени. Иными словами, величина натяжения T
является для рассматриваемой мембраны постоянной.
T ( x, y , z, t ) T0 const .
(21)
Получим теперь уравнение колебаний такой мембраны. Для этого будем
использовать теорему о приращении количества движения.
Пусть S1 - проекция некоторого произвольного участка мембраны на
плоскость ( x , y ) , а С1 - граница этого участка. Приравняем количественное
изменение движения импульса вертикальных составляющих внешних сил,
плотность которых равна F x , y1 , t , и сил натяжения для получения
уравнения колебаний мембраны:
29
S
ut x, y , t2 ut x, y , t1 ( x, y ) dxdy
t2
t2
u
T0
dsdt Fdxdydt
t1 С1
t1
S1
n
(22)
Потребуем, чтобы функция u u ( x, y , t ) имела непрерывные вторые
производные. Тогда применяя теорему Остроградского [21] преобразуем
контурный интеграл в поверхностный:
u
С n ds S (u xx u yy )dxdy
1
1
(23)
Тогда
t2
( x, y )u
t1
S1
tt
T0 (u xx u yy ) F ( x, y , t ) dxdydt 0
(24)
Для вычисления этого интеграла снова воспользуемся теоремой о среднем.
Учитывая, что выбор поверхности интегрирования S1 и временного отрезка
(t1 , t2 ) , можно сказать, что выражение в квадратных скобках тождественно
нулю. Тогда уравнение колебаний мембраны в дифференциальной форме
принимает вид
( x , y )utt T0 (u xx u yy ) F ( x , y , t )
(25)
Полагая, что мембрана является однородной, т.е. ( x , y ) const ,
получим
1
utt u xx u yy F ( x, y, t )
a2
a2
T0
0
(26)
(27)
30
6.
Свободные колебания идеальной круглой мембраны
Рассмотрим полученное ранее уравнение колебаний мембраны в
дифференциальном представлении
1
utt u xx u yy F ( x, y, t )
a2
(28)
Перепишем уравнение с использованием оператора Лапласа
1
utt u F ( x, y, t )
a2
(29)
Будем полагать, что плотность внешних сил F ( x, y, t ) равна нулю, а также для
мембраны заданы некоторые начальные условия
u (r , ,0) f1 ( r , )
ut ( r , ,0) f 2 ( r , )
(30)
В настоящей задаче рассматриваемая мембрана является ограниченной,
а также жестко закрепленной по краям барабанного кожуха, потому на
уравнение (29) наложим соответствующее граничное условие:
u (r0 , , t ) 0
(31)
Таким образом, мы получили уравнение о свободных колебаниях
круглой мембраны. Для удобства будем решать его в полярных координатах
(r , ) . Для этого воспользуемся тем, что оператор Лапласа в полярных
координатах имеет вид
1 u 1 2u
u
r
r r r r 2 2
Тогда уравнение (29) принимает вид:
(32)
31
1 u 1 2u 1 u 2
r
r r r r 2 2 a 2 t 2
(33)
Уравнение (31) является линейным и однородным. В силу этого, его решение
представимо в виде сумм частных решений.
Для решения уравнения (31) воспользуемся методом разделения
переменных: рассмотрим вспомогательную задачу о нахождении решения
уравнения (31), где функция u (r , , t ) выражается при помощи произведения
функции только от переменных r , и функции только от переменного
времени t :
u (r , , t ) v(r , )T (t )
(34)
Граничное условие в этом случае принимает вид
u ( R, , t ) v( R, )T (t ) 0
(35)
Равенство (33) накладывает дополнительное условие на функцию v(r , ) :
v ( R , ) 0
В противном случае, T (t ) 0 , и
(36)
u (r , , t ) 0 , соответственно, что
противоречит начальным условиям (30).
Подставим такую функцию в исходное уравнение (33):
1 v(r , )
1 2 v ( r , )
1 T 2 (t )
T (t )
v(r , ) 2
r
T (t ) 2
r r
r
r
2
a t 2
(37)
Разделим в полученном уравнении переменные:
1 v(r , ) 1 2v(r , )
r
T 2 (t ) r r
r r 2 2
const , 0 (38)
a 2T (t )
v ( r , )
Данные равенства дают следующие уравнения:
32
T (t ) a 2 T 0,
(39)
1 v ( r , ) 1 2 v(r , )
v( r , ) 0
r
r r
r r 2 2
(40)
Константа в уравнении (32) выбрана со знаком минус потому, что при
любом ее значении большем или равном нулю, уравнение (40) имеет только
тривиальное решение [23], что противоречит начальным условиям (30).
Уравнение (39) является стандартным, и имеет решение
T C1 cos(aT ) C2 sin(aT )
(41)
Так как v(r , ) необходимо являться функцией однозначной и
дифференцируемой, то в силу цикличности переменного , от функции
v ( r , )
необходимо потребовать
2 -периодичность.
Таким
образом,
уравнение (34) имеет следующий набор налагаемых на нее условий:
v ( r , ) 0
(42)
v(r , ) v( r , 2 )
(43)
v(0, )
(44)
Для решения этой задачи, снова воспользуемся методом разделения
переменных. Тогда
v(r, ) R(r )( )
(45)
Подставляя такую функцию в уравнение (34), получим
r( )
Разделим
получившееся
d dR
(r ) ( ) R (r ) r 2 0
dr dr
уравнение
на
R(t )( ) ,
соответствующие уравнения для каждой из функций:
и
(46)
рассмотрим
33
2 0
(47)
1 R (r )
2
r
( 2 ) R ( r ) 0
r r r
r
(48)
Уравнения (42-44) принимают следующий вид:
( ) ( 2 ); ( ) ( 2 );
(49)
R(r0 ) 0
(50)
Рассмотрим уравнение (47). Оно имеет нетривиальные решения,
удовлетворяющие
условия
периодичности
(49)
только
тогда,
когда
коэффициент 2 является целым числом - n 2 . В таком случае, решения
принимают вид
( ) D1n cos( n ) D2 n sin( n )
(51)
Рассмотрим теперь уравнение (48). Заметим, что точка r 0 является
для данного уравнения особой. Однако, для особых точек достаточно [20]
наложить требование ограниченности как граничное условие
R (0)
(52)
Решим задачу о собственных значениях для функции R(r ) . Введем
новую переменную x :
x r
(53)
Также обозначим исходную функцию от новой переменной функцией y ( x) :
R( r ) R(
Подставляя
данные
замены
в
x
) y ( x)
уравнение
(48)
(54)
получим
цилиндрических функций относительно функции y ( x) порядка n :
уравнение
34
d 2 y 1 dy n 2
1 y 0
dx 2 x dx x 2
(55)
Граничные условия, преобразованные для функции y ( x) , принимают вид
y ( x0 ) 0, x0 r0
y (0)
(56)
(58)
Уравнение цилиндрических функций имеет общее решение [22] вида
y ( x) d1 J n ( x ) d 2 N n ( x)
(57)
Здесь J n ( x) - функция Бесселя первого рода порядка n , N n ( x ) - функция
Неймана порядка n .
Заметим, что функции Неймана являются бесконечными в точке x 0 . Тогда,
из граничного условия (51) получаем
d2 0
(58)
Для удовлетворения граничного условия (50),
J n ( r0 ) 0
(59)
(n)
Пусть m - n -й корень уравнения J n ( x) 0 . Тогда собственное значение
принимает вид
2
n ,m
m( n )
.
r0
(60)
Такому собственному значению соответствует собственная функция
m( n )
Rn ,m (r ) y ( r ) J n
r
r
0
(61)
35
Отметим, что данные собственные функции являются ортогональными с
весом r :
r0
0
m( n1 ) m( n2 )
rJ n
r Jn
r dr 0,
r
r
0 0
n1 n2 .
(62)
Норма этих функций равна
r0
r0
0
2
m( n )
r02
rJ
r dr J n m( n ) ,
2
r0
(63)
2
m(0)
r02
rJ
r dr J1 m(0)
2
r0
(64)
2
n
и для n 0
0
2
0
Подставляя полученные результаты (51) и (61) в (40), для каждого
собственного значения n ,m существуют две собственные функции
v Jn (
c
v Jn (
s
m( n )
R
m( n )
R
r )cos(n )
(65)
r )sin(n )
(66)
c
Рассмотрим функцию v . Аналогично функциям Rn ,m (r ) , данные
собственные функции ортогональны с весом r , и норма этих функций
выражается как
36
m n1 m n2
2
1
2
0 J n1 r0 r J n2 r0 r rdr 0 cosn1 cos n2 d
0 при n1 n2 ,
0 при n n , m m ,
1
2
1
2
r02 ( n ) 2
J n m при n1 n2 n 0 и m1 m2 m
2
r 2 (n) 2
0 J 0 m 2 при n1 n2 0 и m1 m2 m
2
r0
(67)
s
Вычисления для собственной функции v аналогичны.
Выражение для нормы линейной комбинации собственных функций
m( n )
vn ,m (r , ) J n
r An ,m cos n Bn ,m sin n
r
0
(68)
можно выразить как
2 r0
vn,m (r , )drd
0 0
r02
2
n J n m( n )
2
(69)
где
2, n 0
1, n 0
n
(70)
Таким образом, решение задачи о свободных колебаниях круглой
представляется функцией
37
am( n )
am( n )
u ( r , , t ) v ( r , ) An ,m cos
t Bn ,m sin
t
r
r
n ,m 0
0
0
(71)
am( n )
a m( n )
s
vn ,m (r , ) Cn ,m cos
t Dn ,m sin
t
r
r
n ,m 0
0
0
c
n ,m
Найдем коэффициенты An ,m , Bn ,m , Cn ,m , Dn ,m . Для этого подставим
формулу (71) в начальные условия (30).
A
u (r , ,0)
n ,m 0
ut ( r , ,0)
c
n , m n ,m
v
B
n ,m 0
c
n ,m n ,m
v
Cn ,mvns,m f1 (r , )
D v
s
n ,m n ,m
(72)
am( n )
r f 2 (r , )
0
(73)
Рассмотрим равенство (67). Так как при t 0 синусы обращаются в нуль,
получим
m(0)
m( n )
f1 (r , ) A0,m J 0
r An ,m J n
r cos n
r
r
m 1
0 n1 m1
0
m( n )
Cn ,m J n
r sin n .
r
n 1 m 1
0
Такой ряд является рядом Фурье периодической функции
(74)
f1 (r , ) на
интервале (0,2 ) . Тогда множители при тригонометрических функциях
должны являться коэффициентами Фурье. Для коэффициентов Фурье же
выполняются равенства
1
2
1
2
2
0
0
m( n )
f1 (r , ) d A0,m J 0
r
r
m 1
0
m( n )
f1 (r , )cos n d An ,m J n
r
r
m 1
0
38
1
m( n )
f1 (r , )sin n d Cn ,m J n
r
r
m 1
0
2
0
В то же время, эти равенства являются разложениям некоторой функции в ряд
по функциям Бесселя первого порядка на промежутке (0, r0 ) :
m( n )
F (r ) bn ,m J n
r
r
m 1
0
(75)
Коэффициенты bn ,m такого разложения имеют вид [24]
bn , m
r0
0
m( n )
r f (r ) J n
r dr
r
0
2
r0 ( n ) 2
J n m
2
(76)
Применяя полученные соотношения, получим выражения для коэффициентов
An ,m , Cn,m
An ,m
2
r0
0
0
m( n )
f1 (r , ) J n
r cos n rdrd
r
0
2
r02 m( n )
n Jn
r0
2
r
0
(77)
Аналогично получим выражение для Сn ,m
Cn.m
2
rn
0
0
m( n )
f 2 ( r , ) J n
r sin n rdrd
r
0
2
r02 m( n )
n Jn
r0
2
r
0
Для коэффициентов аналогично Bn ,m , Cn ,m получим
(78)
39
Bn.m
Dn.m
An.m
a n ,m
Cn . m
a n ,m
(79)
(80)
40
7.
Вынужденные колебания идеальной круглой мембраны
Рассмотрим теперь задачу о вынужденных колебаниях конечной
мембраны. Эта задача соответствует второму этапу предлагаемой модели.
Данная задача описывается следующим образом: пусть на закрепленную по
краям мембрану радиуса R действует некоторая внешняя сила. Запишем это в
виде системы уравнений:
1 u 1 2u 1 u 2
f ( r , , t )
r
r r r r 2 2 a 2 t 2
(81)
u ( R , , t ) 0
(82)
u ( r , ,0) f1 ( r , )
ut ( r , ,0) f 2 ( r , )
(83)
Здесь f (r , , t ) обозначает величину внешнюю силу, рассчитанной на единицу
площади мембраны.
Воспользуемся принципом суперпозиции [20] для дифференциальных
уравнений в частных производных, и будем искать решение задачи (81)-(83) в
виде суммы решений двух подзадач: задачи о свободных колебаниях
мембраны с заданными начальными условиями, и о вынужденных колебаниях
мембраны
под
действием
внешней
силы
с
нулевыми
условиями
41
1 v (r , , t ) 1 2v (r , , t ) 1 v (r , , t ) 2
2
r
2
2
r
r
r
r
a
t 2
v ( R , , t ) 0
v ( r , ,0) f ( r , )
1
vt ( r , ,0) f 2 (r , )
(84)
1 w( r , , t ) 1 2 w( r , , t ) 1 w(r , , t ) 2
2
f ( r , , t )
r
2
r
2
a
t 2
r
r r
w( R , , t ) 0
w( r , ,0) 0
wt ( r , ,0) 0
(85)
u ( R , , t ) v ( R , , t ) w( R , , t )
(86)
Решение системы уравнений (84) уже было получено, и представлено
формулой (71).
Рассмотрим систему уравнений (85). Будем искать решения для функции
w( R , , t ) в виде ряда по собственным функциям задачи (84):
w( R , , t )
w
n, m 0
s
n ,m
w
n, m 0
c
n ,m
(t ) J n (
(t ) J n (
nm
R
nm
R
r ) cos(n )
(87)
r )sin(n )
Данное действие возможно в силу того, что набор собственных функций
Jn (
nm
R
r )cos(n ) , J n (
nm
R
r )sin(n ) является ортогональным, и тем самым
формирует базис задачи в круге. Согласно теореме Стеклова [20], всякая
непрерывная функция F ( r , )
c непрерывными первыми и вторыми
производными, удовлетворяющая граничным условиям задачи, может быть
разложена в абсолютно и равномерно сходящийся ряд
42
F (r , ) An ,m vn ,m ( r , ) Bn ,m vn ,m ( r , )
n ,m
по собственным функциям задачи о собственных значениях для круга.
Разложим также по данному базису функцию внешней силы:
f ( r , , t )
n, m 0
c
n
f (t ) J n (
nm
R
r )cos(n ) f (t ) J n (
s
m
nm
R
r )sin(n )
(88)
где
f
f
c
n, m
s
n, m
( r , , t )
( r , , t )
2
r0
0
0
2
r0
0
0
Рассмотрим слагаемое
n, m 0
m( n )
f ( r , , t ) J n
r cos n rdrd
r
0
2
r0 ( n ) 2
n J n m
2
(89)
m( n )
f ( r , , t ) J n
r sin n rdrd
r
0
2
r0 ( n ) 2
n J n m
2
(90)
wnc,m (t ) J n (
nm
R
r ) cos( n ) в выражении (87). Введем
обозначение
S ( r , ) J n (
nm
R
r )cos(n )
(91)
Подставим полученные для данного слагаемого выражения в уравнение (85).
c
2
wn ,m (t ) S ( r , )
wnc,m (t ) S ( r , )
1 n, m 0
12 n , m 0
r
r r
r
2
r
1
a2
n, m 0
wnc,m (t ) S ( r , )
t
2
n, m 0
f nc,m (t ) S ( r , )
(92)
43
Сгруппируем по производным wnc,m (t ) :
1 S ( r , ) 1 2 S ( r , )
w (t ) r r r r r 2 2
n, m 0
1
2 wn,cm (t ) S ( r , ) f nc,m (t ) S ( r , )
a n, m 0
n, m 0
c
n ,m
(93)
Заметим, что функция S (r , ) является решением для рассмотренного ранее
уравнения
1 v(r , ) 1 2v(r , )
v ( r , ) 0
r
r r
r r 2 2
(94)
Учитывая это, уравнение (93) принимает вид
n, m 0
wnc,m (t )n ,m S ( r , )
1
c
w
(
t
)
S
(
r
,
)
f nc,m (t ) S (r , )
n ,m
2
a n, m 0
n, m 0
(95)
где
n ,m
nm
r0
2
(96)
Зафиксируем переменные n, m и вынесем общие множители:
1 c
c
c
wn ,m (t )n ,m 2 wn,m (t ) f n ,m (t ) S ( r , ) 0
a
(97)
Уравнение (27) должно выполняться для произвольных r , , поэтому для его
решения будем решать обычное линейное дифференциальное уравнение
второй степени вида
1 c
wn ,m (t ) n ,m wnc,m (t ) f nc,m (t ) 0
2
a
(98)
44
Рассмотрим функцию f mc (t ) . Ее вид дан формулой (89). В настоящей задаче
подынтегральная функция f (r , , t ) представляется двумя способами: как
некоторая постоянная, и как некоторая тригонометрическая функция с
постоянным множителем. Рассмотрим оба случая по отдельности.
Пусть
f mc (t ) f mc sin(qt )
(99)
где
2
f
c
n, m
P0
1
m( n )
r1 J n r0 r cos n rdrd
r02 ( n ) 2
n J n m
2
r2
(100)
P0 const , q const .
Тем самым получим уравнение
Для
того,
чтобы
1 c
wn,m (t ) n ,m wnc ,m (t ) f nc,m sin(t ) 0
2
a
(101)
решить
сначала
уравнение
(101),
рассмотрим
дифференциальное уравнение для нахождения общего решения.
1 c
wn ,m (t ) n ,m wnc,m (t ) 0
2
a
(102)
Будем искать решение данного уравнения в виде
wnc,m (t ) e t
(103)
Подставляя такую замену в уравнение (101), получим
2
a
2
e t n ,m e t 0
(104)
45
Т.к. уравнение (104) должно выполняться для произвольного конечного , то
получаем
2
n,m 0
(105)
ia n ,m
(106)
a2
Отсюда,
И тогда общее решение принимает вид
wnc,m (t ) c1e
ia n , m
c2 e
ia n , m
(107)
Применяя формулу Эйлера и преобразовывая константы, получим
wnc,m (t ) c1 sin(a n ,m t ) c2 cos( a n ,m t )
(108)
Найдем теперь частное решение уравнения (28). Используя метод
неопределенных коэффициентов, будем искать решение уравнения
1 c
w n,m (t ) n ,m w nc,m (t ) f nc,m sin( qt )
2
a
(109)
в виде
w nc ,m (t ) b1 cos( qt ) b2 sin( qt )
(110)
Вычислим необходимые производные:
d 2 w nc ,m (t )
q 2 b1 cos( qt ) b2 sin( qt )
2
dt
(111)
c
Подставим предполагаемое решение w n ,m (t ) в уравнение (98):
q2
b1 cos(qt ) b2 sin(qt )
a2
n ,m b1 cos( qt ) b2 sin( qt ) f nc,m sin( qt )
(112)
46
Приравняем коэффициенты при функциях cos(qt ) с каждой стороны
уравнения:
q2
2 b1 n ,mb1 0
a
(113)
q2
b n ,mb2 f nc,m sin( qt )
2 2
a
(114)
Аналогично, для sin(qt )
Откуда
b1 0
(115)
a 2 f nc,m
b2 2
a n ,m q 2
(116)
и
w
c
n ,m
a 2 f nc,m
(t ) 2
sin( qt )
a n ,m q 2
(117)
Тогда решением уравнения (31) является
c
n ,m
w
a 2 f nc,m
(t ) c1 sin( a n ,m t ) c2 cos( a n ,m t ) 2
sin( qt )
a n , m q 2
(118)
Найдем теперь константы с1 , с2 из начальных условий, определенных в
(15):
a 2 f nc,m
0 c1 sin( a n ,m 0) c2 cos( a n ,m 0) 2
sin( q 0)
a n , m q 2
(119)
Откуда
c2 0
(120)
47
Вычислим производную wnc,m (t ) :
dwnc,m (t )
c1a n ,m cos( a n ,m t )
dt
a 2 f nc,m
c2 a n ,m sin( a n ,m t ) 2
cos( qt )
a n ,m q 2
(121)
0 c1a n ,m cos( a n ,m 0) c2 a n ,m sin( a n ,m 0)
(122)
a 2 qf nc,m
2
cos(q 0)
a q2
c1
aqf nc,m
( a 2 q 2 ) n , m
(123)
Окончательно, получим
c
n ,m
w
af nc,m 0 q
(t ) 2
sin( a n ,m t ) a sin( qt )
a n , m q 2 n , m
Решение для слагаемого
n, m 0
s
n ,m
w
(t ) J n (
nm
R
(124)
r )sin( n ) в (17) выводится
аналогично, и имеет вид
s
n ,m
w
af ns,m 0 q
(t ) 2
sin( a n ,m t ) a sin( qt )
a n , m q 2 n , m
(125)
где
2
f
s
n, m
P0
1
m( n )
r1 J n r0 r sin n rdrd
r02 ( n ) 2
n J n m
2
r2
(126)
Подставляя выражения (125, 126) в (77), окончательно получим решение
исходной задачи (71)
48
am( n )
am( n )
A
cos
t
B
sin
t
n ,m
n ,m
r
r
0
0
u (r , , t ) vnс,m (r , )
c
af
q
n ,m 0
n ,m
sin(a n ,m t ) a sin(qt )
2
2
a n ,m q n ,m
am( n )
a m( n )
t Dn ,m sin
t
Cn ,m cos
r
r
0
0
vns ,m (r , )
s
q
af
n , m 0
n
,
m
sin( a n ,m t ) a sin(qt )
a 2n ,m q 2 n ,m
(127)
Рассмотрим
теперь
для
слагаемого
n, m 0
c
n ,m
w
(t ) J n (
nm
R
r ) cos( n )
в
выражении (87) случай, когда функция внешней силы f (r , , t ) представима в
виде
f mc (t ) f mc
(128)
c
где f n ,m определяется формулой (90).
Решим теперь обычное линейное дифференциальное уравнение второй
степени вида
1 c
wn,m (t ) n ,m wnc,m (t ) f nc,m 0
2
a
(129)
Решением соответствующего однородного уравнения будет иметь тот же вид,
что и в первом случае:
1 c
wn ,m (t ) n ,m wnc,m (t ) 0
2
a
и его решением будет
(130)
49
wnc,m (t ) c1 sin(a n ,m t ) c2 cos( a n ,m t )
(131)
Найдем теперь частное решение уравнения (129). Используя метод
неопределенных коэффициентов, будем искать решение уравнения
1 c
w n,m (t ) n ,m w nc,m (t ) f nc,m
2
a
(132)
в виде
w nc ,m (t ) b1
(133)
Подставим предполагаемое решение в уравнение:
n ,m b1 f nc,m
(134)
Откуда
w
c
n ,m
(t )
f nc,m
(135)
n , m
Тогда решением (58) является
w (t ) c1 sin( a n ,m t ) c2 cos( a n ,m t )
c
n ,m
f nc,m
n ,m
(136)
Найдем теперь константы с1 , с2 из начальных условий, определенных в (85):
0 c1 sin( a n ,m 0) c2 cos(a n ,m 0)
f nc,m
n ,m
(137)
Откуда
c2
Вычислим необходимые производные:
f nc,m
n ,m
(138)
50
dwnc ,m (t )
c1 cos( a n ,m t ) c2 sin( a n ,m t )
dt
(139)
0 c1 cos( a n ,m 0) c2 sin(a n ,m 0)
(140)
Откуда
c1 0
(141)
Окончательно, получим
w (t )
c
n ,m
f nc,m
n ,m
Решение для слагаемого
cos(a
s
n ,m
w
(t ) J n (
f ns,m
cos(a
n, m 0
n ,m t ) 1
nm
R
(142)
r )sin( n ) в (77) выводится
аналогично, и имеет вид
w (t )
s
n ,m
n ,m
n ,m t ) 1
(143)
где f ns, m определено выражением (126).
Подставляя выражения (142, 143) в (77), окончательно получим решение
исходной задачи (71)
a m( n )
a m( n )
t Bn ,m sin
t
An ,m cos
r
r
0
0
u (r , , t ) vnс,m (r , )
с
f
n ,m 0
n ,m cos(a n ,m t ) 1
n ,m
am( n )
a m( n )
t Dn ,m sin
t
Cn ,m cos
r
r
0
0
vns ,m (r , )
s
f n ,m
n ,m 0
cos(a n ,m t ) 1
n ,m
(144)
51
7.
Выбор входных параметров
7.1. Параметры мембраны
Как было сказано ранее, модель идеальной круглой мембраны,
используемой в настоящей работе, задается тремя параметрами: радиусом,
плотностью материала и величиной натяжения. В соответствии с данными,
доступными в литературе [6,10], натяжение мембраны в расчетах принято
3200 Н/м, а плотность материала мембраны – 2690 кг/м3. Радиус мембраны
принят в соответствии со спецификацией популярного как в оркестре, так и в
обычных барабанных установках рабочего барабана Ludwig [25], имеющим
диаметр мембраны 35.56 см (14 дюймов).
7.2. Параметры палочки
В рамках настоящей работы, такие параметры палочки, как вес или
материал не фигурируют в непосредственных расчетах. Единственным
значимым параметром является ее диаметр. Наиболее распространенный
диаметр круглого наконечника палочки равен 12 мм. Стоит отметить, что в
целях упрощения вычислений смещением мембраны, вызванным контактом с
барабанной
палочкой,
пренебрегается,
однако
фактический
диаметр,
используемый в расчетах, совпадает с номинальным диаметром сферы
наконечника барабанной палочки.
В целях отыскания зависимостей динамики барабанной мембраны от
размера наконечника барабанной палочки, в ходе численного моделирования
будут рассмотрены диаметры 12 мм и 48 мм, что может соответствовать
деревянным колотушкам с круглыми наконечниками подобного размера.
52
7.3. Время взаимодействия
В соответствии с результатами, описанными в главе 3, время
взаимодействия барабанной палочки с мембраной барабана установлено 5 мс.
53
8.
Результаты численного моделирования
В ходе настоящей работы с помощью пакета Wolfram Mathematica 12.1
была разработана программа для численного моделирования динамики
барабанной мембраны при взаимодействии с барабанной палочкой.
8.1
Исходный случай
Рассмотрим сначала случай, используя ступенчатую моедль, где
параметры системы заданы в соответствии с величинами, описанными в
предыдущей главе:
Диаметр мембраны: 35.56 см
Диаметр наконечника палочки: 12 мм
Плотность мембраны: 2690 кг/м3
Натяжение мембраны: 3200 Н/м
Время взаимодействия: 5 мс
Сила удара: 100 Н
Для проведения анализа динамики мембраны в зависимости от точки
удара, а также параметров системы, будем сравнивать абсолютные значения
коэффициентов при симметричных и антисимметричных модах, данные в
формуле (71). Будем оценивать их величины в семи разных временных точках:
th – время, равное половине времени взаимодействия
t0-t5 – время, равное 0 - 500 мс после окончания
взаимодействия соответственно
Будем последовательно рассматривать удары по отдалению от центра
мембраны.
54
Рассмотрим удар в середину мембраны. На рисунках 13 и 14
представлены
графики
абсолютных
величин
коэффициентов
при
симметричных и антисимметричных модах от времени соответственно.
0 000016
0 000014
0 000012
0 00001
0 000008
0 000006
0 000004
0 000002
0
th
(0,1)
(1,1)
(2,1)
Рис. 13.
t0
(0,2)
(3,1)
t1
(1,2)
(4,1)
t2
(2,2)
(0,3)
(5,1)
t3
(3,2)
(6,1)
t4
(7, 1)
(7, 2)
t5
(7, 3)
(7, 4)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
ударе в центр мембраны от времени
0 006
0 005
0 004
0 003
0 002
0 001
0
th
(0,1)
(1,1)
t0
(2,1)
(0,2)
(3,1)
t1
(1,2)
(4,1)
t2
(2,2)
(0,3)
t3
(5,1)
(3,2)
(6,1)
t4
(7, 1)
(7, 2)
t5
(7, 3)
(7, 4)
(8,3)
55
Рис. 14.
График зависимости величин коэффициентов при антисимметричных модах
при ударе в центр мембраны от времени
Из представленных графиков видно, что при ударе в центр мембраны
возбуждение испытывают только лишь радиальные моды (0, 1), (0, 2), (0, 3) и
так далее. При этом моды высших порядков в начальный момент испытывают
большее возбуждение в сравнении с модами низших порядков, однако первые
сравнительно быстро затухают.
Зависимости величин коэффициентов для удара, смещенного на 5 см от
центра мембраны, представлены на графиках на рисунках 15 и 16:
0 00003
0 000025
0 00002
0 000015
0 00001
0 000005
0
th
(0,1)
(1,1)
Рис. 15.
(2,1)
t0
(0,2)
(3,1)
t1
(1,2)
(4,1)
t2
(2,2)
(0,3)
(5,1)
t3
(3,2)
(6,1)
t4
(7, 1)
(7, 2)
t5
(7, 3)
(7, 4)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=5 см от времени
Из графиков видно, что при смещении удара от центра возбуждается
множество новых мод, как низких, так и высоких порядков. При этом
абсолютные значения симметричных и антисимметричных коэффициентов
возрастают на порядок. Моды (1, 1), (1, 2), (3, 1) и (2, 2) выступают как моды
с наибольшей энергией на протяжении всего времени. Для симметричных мод
также больший вклад дают моды высших порядков, как, например, мода (7,4).
Радиальные моды при таком ударе практически не возбуждаются.
56
Графики на рисунках 17 и 18 показывают динамику модальных
коэффициентов при ударе, смещенном на 10 см от центра мембраны:
0 01
0 009
0 008
0 007
0 006
0 005
0 004
0 003
0 002
0 001
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(1,2)
(4,1)
(2,2)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
(7, 3)
(7, 4)
(8,3)
Рис. 16.
t2
t3
t4
t5
(0,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=5 см от времени
0 00004
0 000035
0 00003
0 000025
0 00002
0 000015
0 00001
0 000005
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
t2
(1,2)
(7, 3)
t3
(4,1)
(7, 4)
t4
(2,2)
(8,3)
t5
(0,3)
57
Рис. 17.
График зависимости величин коэффициентов при симметричных модах при
смещении d=10 см от времени
0 012
0 01
0 008
0 006
0 004
0 002
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(1,2)
(4,1)
(2,2)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
(7, 3)
(7, 4)
(8,3)
Рис. 17.
t2
t3
t4
t5
(0,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=10 см от времени
Из графиков видно, что при таком смещенном ударе мода (1, 2) заметно
ослабевает по сравнению с предыдущим случаем. При этом равномерно
возбуждаются диаметральные моды (1, 1), (2, 1), (3, 1), (4, 1), (5, 1), (6, 1), (7,
1), а также моды высших порядков (7, 4) и (8, 3). По прошествии 300 мс
высшие моды постепенно начинают затухать, и превалируют моды
диаметральные моды низких порядков.
Наконец, на графиках на рисунках 18 и 19 представлены графики
зависимости модальных коэффициентов от времени для удара, смещенного на
15 см от центра мембраны. Из графиков видно, что при таком ударе
возбуждается большое количество мод высших порядков. Первые 200 мс
наибольший вклад вносят смешанные моды (3, 2), (7, 2), (7, 3), (7, 4), а также
диаметральные моды (6, 1) и (7, 1). К моменту времени t5 = 500 мс большая
часть мод высших порядков, кроме (7, 1), ослабевает, сравниваясь по величине
58
с модами низших порядков (2,1), (3,1), (2,2). Коэффициент при моде (1, 1)
имеет наименьшее значение из всех трех случаев.
0 000035
0 00003
0 000025
0 00002
0 000015
0 00001
0 000005
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
Рис. 18.
t2
(1,2)
t3
(7, 3)
(4,1)
t4
(7, 4)
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=15 см от времени
0 01
0 009
0 008
0 007
0 006
0 005
0 004
0 003
0 002
0 001
0
(0,1)
(5,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
t2
(3,2)
(6,1)
(7, 1)
(7, 2)
(1,2)
(7, 3)
t3
(4,1)
(7, 4)
t4
(2,2)
(8,3)
t5
(0,3)
59
Рис. 19.
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=15 см от времени
Подытоживая полученные результаты, можно выделить следующие
особенности:
При
отдалении
от
центра
мембраны,
интенсивность
фундаментальных мод (1, 1) и (1, 2) ослабевает
При отдалении от центра, интенсивность радиальных мод
уменьшается, возрастают диаметральные моды низких порядков
При отдалении от центра наблюдается уменьшение в абсолютных
значениях модальных коэффициентов
При отдалении от центра повышается интенсивность смешанных
мод высших порядков
Амплитуды мод, возбуждаемых внешней силой, меньше на два
порядка, чем при возбуждении передачей импульса
8.2. Увеличенное время контакта
Рассмотрим теперь случай с увеличенным временем контакта.
Параметры системы:
Диаметр мембраны: 35.56 см
Диаметр наконечника палочки: 12 мм
Плотность мембраны: 2690 кг/м3
Натяжение мембраны: 3200 Н/м
Время взаимодействия: 15 мс
Сила удара: 100 Н
На рисунках 20-21 представлены графики абсолютных величин
коэффициентов при симметричных и антисимметричных модах от времени
при нанесении удара в центр мембраны при такой конфигурации.
60
0 00018
0 00016
0 00014
0 00012
0 0001
0 00008
0 00006
0 00004
0 00002
0
th
(0,1)
(1,1)
t0
(2,1)
Рис. 20.
(0,2)
(3,1)
t1
(1,2)
(4,1)
t2
(2,2)
(0,3)
(5,1)
t3
(3,2)
(6,1)
t4
(7, 1)
t5
(7, 2)
(7, 3)
(7, 4)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
ударе в центр мембраны, повышенное время контакта
0 025
0 02
0 015
0 01
0 005
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(1,2)
(4,1)
(2,2)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
(7, 3)
(7, 4)
(8,3)
Рис. 21.
t2
t3
t4
t5
(0,3)
График зависимости величин коэффициентов при антисимметричных модах
при ударе в центр мембраны, повышенное время контакта
61
0 00018
0 00016
0 00014
0 00012
0 0001
0 00008
0 00006
0 00004
0 00002
0
th
(0,1)
(1,1)
t0
(2,1)
Рис. 22.
(0,2)
(3,1)
t1
(1,2)
(4,1)
t2
(2,2)
(0,3)
(5,1)
t3
(3,2)
(6,1)
t4
(7, 1)
t5
(7, 2)
(7, 3)
(7, 4)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=5 см от времени, повышенное время контакта
0 025
0 02
0 015
0 01
0 005
0
(0,1)
(5,1)
Рис. 23.
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
t2
(3,2)
(6,1)
(7, 1)
(7, 2)
(1,2)
(7, 3)
t3
(4,1)
(7, 4)
t4
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=5 см от времени, повышенное время контакта
62
0 00025
0 0002
0 00015
0 0001
0 00005
0
th
(0,1)
(5,1)
Рис. 24.
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
t2
(1,2)
t3
(7, 3)
(4,1)
t4
(7, 4)
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=10 см от времени, повышенное время контакта
0 025
0 02
0 015
0 01
0 005
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(1,2)
(4,1)
(2,2)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
(7, 3)
(7, 4)
(8,3)
Рис. 25.
t2
t3
t4
t5
(0,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=10 см от времени, повышенное время контакта
63
0 00018
0 00016
0 00014
0 00012
0 0001
0 00008
0 00006
0 00004
0 00002
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(1,2)
(4,1)
(2,2)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
(7, 3)
(7, 4)
(8,3)
Рис. 26.
t2
t3
t4
t5
(0,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=15 см от времени, повышенное время контакта
0 025
0 02
0 015
0 01
0 005
0
(0,1)
(5,1)
Рис. 27.
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
t2
(3,2)
(6,1)
(7, 1)
(7, 2)
(1,2)
(7, 3)
t3
(4,1)
(7, 4)
t4
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=15 см от времени, повышенное время контакта
Графики 20-27 показывают, что увеличение времени контакта не вносит
значительных изменений во временную динамику возбуждаемых мод. При
64
этом значения коэффициентов при всех возбужденных модах увеличилось в
некоторое количество раз. Для его оценки, рассмотрим отношение
коэффициентов, например, при симметричной и антисимметричной моде
(7, 1) в случае увеличенного времени контакта и его значения в исходном
случае соответственно:
k 7 c,1
k
c
7,1
k 7 ,1s
s
k7,1
0.000169311
5.8
0.0000291585
0.020189666
2.1
0.009406022
8.3. Увеличенный радиус наконечника палочки
Рассмотрим теперь случай с увеличенным радиусом наконечника
барабанной палочки. Параметры системы:
Диаметр мембраны: 35.56 см
Диаметр наконечника палочки: 48 мм
Плотность мембраны: 2690 кг/м3
Натяжение мембраны: 3200 Н/м
Время взаимодействия: 5 мс
Сила удара: 100 Н
На рисунках 28 и 29 представлены графики абсолютных величин
коэффициентов при симметричных и антисимметричных модах от времени
при нанесении удара в центр мембраны при такой конфигурации.
65
0 000025
0 00002
0 000015
0 00001
0 000005
0
th
(0,1)
(1,1)
(2,1)
Рис. 28.
t0
(0,2)
(3,1)
t1
(1,2)
(4,1)
t2
(2,2)
(0,3)
(5,1)
t3
(3,2)
(6,1)
t4
(7, 1)
t5
(7, 2)
(7, 3)
(7, 4)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
ударе в центр мембраны, увеличенный радиус наконечника палочки
0 008
0 007
0 006
0 005
0 004
0 003
0 002
0 001
0
(0,1)
(5,1)
Рис. 29.
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
t2
(3,2)
(6,1)
(7, 1)
(7, 2)
(1,2)
(7, 3)
t3
(4,1)
(7, 4)
t4
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при антисимметричных модах
при ударе в центр мембраны, увеличенный радиус наконечника палочки
Графики 28 и 29 показывают, что увеличение радиуса наконечника палочки не
вносит значительных изменений во временную динамику возбуждаемых мод.
66
При этом значения коэффициентов при всех возбужденных модах
увеличилось на некоторую равную примерно 2*10-3 для всех модальных
коэффициентов
Зависимости величин коэффициентов для удара, смещенного на 5 см от
центра мембраны, представлены на графиках на рисунках 30 и 31:
0 00002
0 000018
0 000016
0 000014
0 000012
0 00001
0 000008
0 000006
0 000004
0 000002
0
th
(0,1)
(1,1)
Рис. 30.
(2,1)
t0
(0,2)
(3,1)
t1
(1,2)
(4,1)
t2
(2,2)
(0,3)
(5,1)
t3
(3,2)
(6,1)
t4
(7, 1)
(7, 2)
t5
(7, 3)
(7, 4)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=5 см от времени, увеличенный радиус наконечника палочки
Из графиков видно, что при таком ударе по сравнению с исходным
случаем уменьшается интенсивность мод высших порядков по отношению к
смешанным модам (1, 2), (2, 2) и (3, 2). Так же наблюдается возрастание по
величине всех модальных коэффициентов.
Из графиков, представленных на рисунках 32 и 33, видно, что для
удара, смещенного на 10 см от центра, по сравнению с исходным случаем
происходит приглушение всех возбужденных мод, а также мод высших
порядков в сравнении с модами низших порядков. Смешанная мода (7, 2)
остается при этом превалирующей.
67
0 007
0 006
0 005
0 004
0 003
0 002
0 001
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(1,2)
(4,1)
(2,2)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
(7, 3)
(7, 4)
(8,3)
Рис. 31.
t2
t3
t4
t5
(0,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=5 см от времени, увеличенный радиус наконечника
палочки
0 000025
0 00002
0 000015
0 00001
0 000005
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
Рис. 32.
t2
(1,2)
(7, 3)
t3
(4,1)
(7, 4)
t4
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=10 см от времени, увеличенный радиус наконечника палочки
68
0 009
0 008
0 007
0 006
0 005
0 004
0 003
0 002
0 001
0
(0,1)
(5,1)
Рис. 33.
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
t2
(3,2)
(6,1)
(7, 1)
(7, 2)
(1,2)
t3
(7, 3)
(4,1)
t4
(7, 4)
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=10 см от времени, увеличенный радиус наконечника
палочки
0 000025
0 00002
0 000015
0 00001
0 000005
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
Рис. 34.
t2
(1,2)
(7, 3)
t3
(4,1)
(7, 4)
t4
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=15 см от времени, увеличенный радиус наконечника палочки
69
0 006
0 005
0 004
0 003
0 002
0 001
0
th
(0,1)
(5,1)
Рис. 35.
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
t2
(3,2)
(6,1)
(7, 1)
(7, 2)
(1,2)
(7, 3)
t3
(4,1)
(7, 4)
t4
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=15 см от времени, увеличенный радиус наконечника
палочки
Аналогично, видно, что для удара, смещенного на 10 см от центра, удар,
смещенный на 15 см с увеличенным наконечником палочки, по сравнению с
исходным случаем, производит глушение всех возбужденных мод, а также мод
высших порядков в сравнении с модами низших порядков.
Таким образом, для случая с увеличенным радиусом наконечника,
можно выделить следующие особенности:
При ударе в центр амплитуда возбуждаемых мод повышается
При смещенных от центра ударах происходит приглушение мод
высших порядков и уменьшение амплитуд возбуждаемых мод по
сравнению с исходным случаем
8.4. Уменьшенная величина силы удара
Рассмотрим теперь случай с уменьшенной силой удара. Параметры
системы:
70
Диаметр мембраны: 35.56 см
Диаметр наконечника палочки: 12 мм
Плотность мембраны: 2690 кг/м3
Натяжение мембраны: 3200 Н/м
Время взаимодействия: 5 мс
Сила удара: 60 Н
0 000018
0 000016
0 000014
0 000012
0 00001
0 000008
0 000006
0 000004
0 000002
0
th
(0,1)
(1,1)
Рис. 36.
(2,1)
t0
(0,2)
(3,1)
t1
(1,2)
(4,1)
t2
(2,2)
(0,3)
(5,1)
t3
(3,2)
(6,1)
t4
(7, 1)
(7, 2)
t5
(7, 3)
(7, 4)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
ударе в центр мембраны от времени, уменьшенная величина силы
71
0 006
0 005
0 004
0 003
0 002
0 001
0
th
(0,1)
(5,1)
Рис. 37.
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
t2
(3,2)
(6,1)
(7, 1)
(7, 2)
(1,2)
t3
(7, 3)
(4,1)
t4
(7, 4)
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при антисимметричных модах
при ударе в центр мембраны от времени, уменьшенная величина силы
0 000018
0 000016
0 000014
0 000012
0 00001
0 000008
0 000006
0 000004
0 000002
0
th
(0,1)
(1,1)
Рис. 38.
(2,1)
t0
(0,2)
(3,1)
t1
(1,2)
(4,1)
t2
(2,2)
(0,3)
(5,1)
t3
(3,2)
(6,1)
t4
(7, 1)
(7, 2)
t5
(7, 3)
(7, 4)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=5 см от времени, уменьшенная величина силы
72
0 006
0 005
0 004
0 003
0 002
0 001
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(1,2)
(4,1)
(2,2)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
(7, 3)
(7, 4)
(8,3)
Рис. 39.
t2
t3
t4
t5
(0,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=5 см от времени, уменьшенная величина силы
0 000025
0 00002
0 000015
0 00001
0 000005
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
Рис. 40.
t2
(1,2)
(7, 3)
t3
(4,1)
(7, 4)
t4
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=10 см от времени, уменьшенная величина силы
73
0 008
0 007
0 006
0 005
0 004
0 003
0 002
0 001
0
(0,1)
(5,1)
Рис. 41.
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
t2
(3,2)
(6,1)
(7, 1)
(7, 2)
(1,2)
t3
(7, 3)
(4,1)
t4
(7, 4)
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=10 см от времени, уменьшенная величина силы
0 00002
0 000018
0 000016
0 000014
0 000012
0 00001
0 000008
0 000006
0 000004
0 000002
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
Рис. 42.
t2
(1,2)
(7, 3)
t3
(4,1)
(7, 4)
t4
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=15 см от времени, уменьшенная величина силы
74
0 006
0 005
0 004
0 003
0 002
0 001
0
(0,1)
th
(5,1)
Рис. 43.
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
t2
(3,2)
(6,1)
(7, 1)
(7, 2)
(1,2)
t3
(4,1)
(7, 3)
(7, 4)
t4
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=15 см от времени, уменьшенная величина силы
Из графиков, представленных на рисунках 36-43 видно, что уменьшение
силы не вносит изменений во временную динамику возбуждаемых мод. Более
того,
данные
на
представленных
графиках
соответствуют
данным,
представленным на графиках на рисунках 13-19 с точностью до множителя для
коэффициентов при симметричных и антисимметричных модах.
8.5.
Тригонометрическая модель
Рассмотрим теперь взаимодействие мембраны с барабанной палочкой,
использую
тригонометрическую
модель
с
параметрами
идентичными исходному случаю ступенчатой модели:
Диаметр мембраны: 35.56 см
Диаметр наконечника палочки: 12 мм
Плотность мембраны: 2690 кг/м3
Натяжение мембраны: 3200 Н/м
системы,
75
Время взаимодействия: 5 мс
Сила удара: 100 Н
На рисунках 44-51 представлены графики абсолютных величин
коэффициентов при симметричных и антисимметричных модах от времени
при нанесении удара в центр мембраны при такой конфигурации.
0 000018
0 000016
0 000014
0 000012
0 00001
0 000008
0 000006
0 000004
0 000002
0
th
(0,1)
(1,1)
Рис. 44.
(2,1)
t0
(0,2)
(3,1)
t1
(1,2)
(4,1)
t2
(2,2)
(0,3)
(5,1)
t3
(3,2)
(6,1)
t4
(7, 1)
(7, 2)
t5
(7, 3)
(7, 4)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
ударе в центр мембраны от времени, тригонометрическая модель
76
0 006
0 005
0 004
0 003
0 002
0 001
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(1,2)
(4,1)
(2,2)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
(7, 3)
(7, 4)
(8,3)
Рис. 45.
t2
t3
t4
t5
(0,3)
График зависимости величин коэффициентов при антисимметричных модах
при ударе в центр мембраны от времени, тригонометрическая модель
0 000009
0 000008
0 000007
0 000006
0 000005
0 000004
0 000003
0 000002
0 000001
0
th
(0,1)
(1,1)
Рис. 46.
(2,1)
t0
(0,2)
(3,1)
t1
(1,2)
(4,1)
t2
(2,2)
(0,3)
(5,1)
t3
(3,2)
(6,1)
t4
(7, 1)
(7, 2)
t5
(7, 3)
(7, 4)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=15 см от времени, уменьшенная величина силы
77
0 003
Коэффициенты при антисимметричных модах
при d = 0.1
0 0025
0 002
0 0015
0 001
0 0005
0
(0,1)
(5,1)
Рис. 47.
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
t2
(1,2)
t3
(7, 3)
(4,1)
t4
(7, 4)
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=15 см от времени, тригонометрическая модель
0 000012
Коэффициенты при симметричных модах при d
= 0.1
0 00001
0 000008
0 000006
0 000004
0 000002
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
Рис. 48.
t2
(1,2)
(7, 3)
t3
(4,1)
(7, 4)
t4
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=15 см от времени, уменьшенная величина силы
78
0 003
Коэффициенты при антисимметричных модах
при d = 0.1
0 0025
0 002
0 0015
0 001
0 0005
0
th
(0,1)
(1,1)
t0
(2,1)
(5,1)
(3,2)
(6,1)
Рис. 49.
t1
(0,2)
t2
t3
t4
(3,1)
(1,2)
(4,1)
(2,2)
(7, 1)
(7, 2)
(7, 3)
(7, 4)
(8,3)
t5
(0,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=15 см от времени, тригонометрическая модель
0 000009
0 000008
0 000007
0 000006
0 000005
0 000004
0 000003
0 000002
0 000001
0
(0,1)
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(1,2)
(4,1)
(2,2)
(5,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
(7, 3)
(7, 4)
(8,3)
Рис. 50.
t2
t3
t4
t5
(0,3)
График зависимости величин коэффициентов при симметричных модах при
смещении d=15 см от времени, уменьшенная величина силы
79
0 003
0 0025
0 002
0 0015
0 001
0 0005
0
(0,1)
(5,1)
Рис. 51.
th
(1,1)
t0
(2,1)
t1
(0,2)
(3,1)
(3,2)
(6,1)
(7, 1)
(7, 2)
t2
(1,2)
(7, 3)
t3
t4
(4,1)
(7, 4)
(2,2)
t5
(0,3)
(8,3)
График зависимости величин коэффициентов при антисимметричных модах
при смещении d=15 см от времени, тригонометрическая модель
Из графиков, представленных
на
рисунках
40-51 видно,
что
использование тригонометрической модели взаимодействия не вносит
изменений во временную динамику возбуждаемых мод. Более того, данные на
представленных графиках соответствуют данным, представленным на
графиках на рисунках 13-19 с точностью до множителя для коэффициентов
при симметричных и антисимметричных модах.
Таким образом, можно сказать, что предложенные ступенчатая и
тригонометрическая модели эквивалентны. В связи с этим, остальные случаи
конфигурации системы для тригонометрической модели в настоящей работе
рассмотрены не будут.
80
9.
Заключение
В представленной работе было изучено взаимодействие барабанной
палочки с барабанной мембраной. В частности, была исследована динамика
колебаний первых девятнадцати мод идеальной круглой мембраны при
различном удалении наносимого удара от центра мембраны, а также при
различных характеристиках удара; проведены сравнительные характеристики.
В результате моделирования были получены следующие основные
результаты:
При удалении точки удара от центра возрастает интенсивность
смешанных мод высших порядков, при этом интенсивность мод низших
порядков,
а
также
радиальных
мод,
снижается;
наблюдается
уменьшение абсолютных величин модальных коэффициентов при
отдалении от центра
Вариации таких параметров, как время контакта палочки с мембраной
или силы удара, не вносят значительных изменений во временную
динамику возбуждаемых мод, но увеличивает или уменьшает амплитуду
возбуждаемых колебаний при соответствующих изменениях данных
параметров
При взаимодействиях с длительностью меньшей 20 мс амплитуды мод,
возбуждаемых внешней силой, меньше на два порядка, чем при
возбуждении передачей импульса
Уменьшение отношения диаметра наконечника палочки к диаметру
мембраны приводит к увеличению интенсивности возбуждаемых мод
при ударе в центр мембраны, но к уменьшению при ударах, смещенных
от центра
Полученные результаты находятся в хорошем соответствии с физическими
экспериментами и с музыкальной практикой.
81
В свете полученных результатов интерес также могла бы представлять
модель мембраны с дополнительными внутренними граничными условиями в
области контакта, вызванными продолжительным давлением палочки.
Аналитическое исследование этого вопроса затрудняется существенной
сложностью получения явных выражений. В связи с этим в работе было
проведено только численное моделирование.
82
10.
Список литературы
[1] - Strickland, Edward. "Glass, Philip." The New Grove Dictionary of Music and
Musicians. Edited by Stanley Sadie and John Tyrrell. London: Macmillan, 2001. ix:
932-936
[2] - Rayleigh Lord (1945). Theory of sound, 2 nd ed, Vol. I, Macmillian, New York
(1894). Reprinted by Dover, New York
[3] - Fletcher H, Bassett IG (1978). Some experiments with the bass drum. J Acoust
Soc Am 64/6: 1570- 1576
[4] - R. S. Christian, R. E. Davis, A. Tubis, C. A. Anderson, R. I.Mills and T. D.
Rossing, “Effects of air loading on timpani membrane vibrations”, J. Acoust. Soc.
Am. 76, 1336–1345 (1984)
[5] - Rossing T, Bork I, Zhao H & Fystrom DO (1992). Acoustics of snare drums. J
Acoust Soc Am 92/1: 84-94.
[6] - T. D. Rossing, 114th meeting, Acoust. Soc. Am.; abstract in J Acoust. Soc. Am.
82, S69 (1987)
[7] - T. Rossing, “Acoustics of percussion instruments: Recent progress,” Acoust.
Sci. & Tech. 22, 177–188 (2001).
[8] - Rossing, T.D. “The Science of Sound,” Addison-Wesley, Reading,
Massachusetts (1982)
[9] - - Rhaouti, A. Chaigne, and P. Joly, “Time-domain modeling and numerical
simulation of a kettledrum,” J. Acoust. Soc. Am. 105, 3545–3562 (1999)
[10] - Bilbao, Stefan & Webb, Craig. (2013). Physical modeling of timpani drums
in 3D on GPGPUs. Journal of the Audio Engineering Society. 61. 737-748
83
[11] - Кезик В. И. “Эквивалентные параметры круглой однородной мембраны”
Техническая акустика том 14, 1819-2408 (2015).
[12] - A. Chaigne and A. Askenfelt, “Numerical simulations of struck strings. I. A
physical model for a struck string using finite difference methods,” J. Acoust. Soc.
Am. 95, 1112–1118 (1994)
[13] - K. van den Doel and D. Pai, “Modal synthesis for vibrating objects,” in Audio
Anecdotes III, edited by K. Greenebaum (A. K. Peters, Natick, Massachusetts, 2007)
[14] - Edgar Berdahl, Bill Verplank, Julius O. Smith III, and Gunter Niemeyer, “A
physically-intuitive haptic drumstick” In Proc. International Computer Music Conf,
Vol. 1, pp. 363-366 (2007)
[15] - C. Jin, A. Davoodabadi, J. Li, Y. Wang, T. Singler, Spherical indentation of a
freestanding circular membrane revisited: analytical solutions and experiments. J.
Mech. Phys. Solids 100, 2017, 85
[16] - Jin, Congrui. (2009). Theoretical study of mechanical behavior of thin circular
film adhered to a flat punch. International Journal of Mechanical Sciences - INT J
MECH SCI. 51. 481-489. 10.1016/j.ijmecsci.2009.04.003.
[17] - Begley, Matthew & Mackin, Thomas. (2004). Spherical indentation of
freestanding circular thin films in the membrane regime. Journal of The Mechanics
and Physics of Solids - J MECH PHYS SOLIDS. 52. 2005-2023.
10.1016/j.jmps.2004.03.002.
[18] - Andreas Wagner. “Analysis of Drumbeats – Interaction between Drummer,
Drumstick and Instrument”, Master’s Thesis at the Department of Speech, Music
and Hearing (TMH), 2005
[19] - S. Dahl. Spectral changes in the tom-tom related to striking force. Speech
Music and Hearing Quarterly Progress and Status Report, 38(1):059–065, 1997.
[20] - А. И. Тихонов, А. А. Самарский, Уравнение математической физики,
Наука, М., 1972
84
[21] - Смирнов В. И., Курс высшей математики, т. 2, Наука, М., 1967, стр. 196
[22] - G. N. Watson, Theory of Bessel functions, Cambridge Univ. Press,
Cambridge, 1944, стр. 106
[23]
-
Кошляков
Н.С.,
Глинер
Э.Б.,
Смирнов
М.М., Основные
дифференциальные уравнения математической физики, Гос. изд-во ФМЛ, М.,
1962
[24] - https://en.wikipedia.org/wiki/Drumhead
[25]
-
drums/supraphonic
https://www.ludwig-drums.com/en-us/ludwig/products/snare-
Отзывы:
Авторизуйтесь, чтобы оставить отзыв