Вибрации круговой мембраны

редактировать
Один из возможных режимов вибрации идеализированной круглой головки барабана (режим u 12 {\ displaystyle u_ {12}}u_ {12 } с обозначениями ниже). Другие возможные режимы показаны внизу статьи.

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

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

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

Содержание
  • 1 Мотивация
  • 2 Проблема
  • 3 Осесимметричный случай
  • 4 Общий случай
  • 5 Анимация нескольких режимов вибрации
  • 6 См. Также
  • 7 Ссылки
Мотивация

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

Проблема

Рассмотрим открытый диск Ω {\ displaystyle \ Omega}\ Omega радиуса a {\ displaystyle a }a с центром в начале координат, который будет представлять "неподвижную" форму головки барабана. В любой момент t, {\ displaystyle t,}t, высота формы головки барабана в точке (x, y) {\ displaystyle (x, y)}(x, y) в Ω {\ displaystyle \ Omega}\ Omega , измеренное от "неподвижной" формы головки барабана, будет обозначаться u (x, y, t), {\ displaystyle u ( x, y, t),}u (x, y, t), , которые могут принимать как положительные, так и отрицательные значения. Пусть ∂ Ω {\ displaystyle \ partial \ Omega}\ partial \ Omega обозначает границу из Ω, {\ displaystyle \ Omega,}\Omega,то есть, круг радиуса a {\ displaystyle a}a с центром в исходной точке, который представляет жесткий каркас, к которому прикреплена головка барабана.

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

∂ 2 u ∂ t 2 = c 2 (∂ 2 u ∂ x 2 + ∂ 2 u ∂ y 2) для (x, y) ∈ Ω {\ displaystyle {\ frac {\ partial ^ {2} u} {\ partial t ^ {2}}} = c ^ {2} \ left ({\ frac {\ partial ^ {2} u} {\ partial x ^ {2}}} + {\ frac {\ partial ^ {2} u} {\ partial y ^ {2}}} \ right) {\ text {for}} (x, y) \ in \ Omega \,}{\ frac {\ partial ^ {2} u} {\ partial t ^ {2}}} = c ^ {2} \ left ({\ frac {\ partial ^ {2} u} {\ partial x ^ {2}}} + {\ frac {\ partial ^ {2} u} {\ partial y ^ {2}}} \ right) {\ text {for}} (x, y) \ in \ Omega \,
u = 0 на ∂ Ω. {\ displaystyle u = 0 {\ text {on}} \ partial \ Omega. \,}u = 0 {\ text {on}} \ partial \ Omega. \,

Из-за круговой геометрии Ω {\ displaystyle \ Omega}\ Omega он будет удобно использовать цилиндрические координаты, (r, θ, z). {\ displaystyle (r, \ theta, z).}(r, \ theta, z). Тогда приведенные выше уравнения записываются как

∂ 2 u ∂ t 2 = c 2 (∂ 2 u ∂ r 2 + 1 r ∂ u ∂ r + 1 r 2 ∂ 2 u ∂ θ 2) для 0 ≤ r < a, 0 ≤ θ ≤ 2 π {\displaystyle {\frac {\partial ^{2}u}{\partial t^{2}}}=c^{2}\left({\frac {\partial ^{2}u}{\partial r^{2}}}+{\frac {1}{r}}{\frac {\partial u}{\partial r}}+{\frac {1}{r^{2}}}{\frac {\partial ^{2}u}{\partial \theta ^{2}}}\right){\text{ for }}0\leq r{\ frac {\ partial ^ {2} u} {\ partial t ^ {2}}} = c ^ {2 } \ left ({\ frac {\ partial ^ {2} u} {\ partial r ^ {2}}} + {\ frac {1} {r}} {\ frac {\ partial u} {\ partial r} } + {\ frac {1} {r ^ {2}}} {\ frac {\ partial ^ {2} u} {\ partial \ theta ^ {2}}} \ right) {\ text {for}} 0 \ leq r <a, 0 \ leq \ theta \ leq 2 \ pi \,
u = 0 для r = a. {\ displaystyle u = 0 {\ text {for}} r = a. \,}u = 0 {\ text {for}} r = a. \,

Здесь c {\ displaystyle c}c - положительная константа, которая дает скорость при поперечные колебательные волны распространяются в мембране. С точки зрения физических параметров скорость волны c определяется выражением

c = N rr ∗ ρ h {\ displaystyle c = {\ sqrt {\ frac {N_ {rr} ^ {*}} {\ rho h}}}}c = {\ sqrt {{\ frac {N _ {{rr}} ^ {*} } {\ rho h}}}}

где N rr ∗ {\ displaystyle N_ {rr} ^ {*}}N _ {{rr}} ^ {* } , - это радиальная мембрана, равная на границе мембраны (r = a {\ displaystyle r = a}r = a ), h {\ displaystyle h}h - толщина мембраны, а ρ {\ displaystyle \ rho}\ rho - плотность мембраны. Если мембрана имеет равномерное натяжение, сила равномерного натяжения при заданном радиусе r {\ displaystyle r}r может быть записана как

F = r N rrr = r N θ θ r {\ displaystyle F = rN_ {rr} ^ {r} = rN _ {\ theta \ theta} ^ {r}}F = rN _ {{rr}} ^ { {r}} = rN _ {{\ theta \ theta}} ^ {{r}}

где N θ θ r = N rrr {\ displaystyle N _ {\ theta \ theta} ^ { r} = N_ {rr} ^ {r}}N _ {\ theta \ theta}} ^ {{r}} = N _ {{rr}} ^ {{r}} - результирующая мембрана в азимутальном направлении.

Осесимметричный случай

Сначала мы изучим возможные режимы колебаний круглой головки барабана, которые являются осесимметричными. Тогда функция u {\ displaystyle u}u не зависит от угла θ, {\ displaystyle \ theta,}\ theta, , и волновое уравнение упрощается до

∂ 2 u ∂ t 2 = c 2 (∂ 2 u ∂ r 2 + 1 r ∂ u ∂ r). {\ displaystyle {\ frac {\ partial ^ {2} u} {\ partial t ^ {2}}} = c ^ {2} \ left ({\ frac {\ partial ^ {2} u} {\ partial r ^ {2}}} + {\ frac {1} {r}} {\ frac {\ partial u} {\ partial r}} \ right).}{\ frac {\ partial ^ { 2} u} {\ partial t ^ {2}}} = c ^ {2} \ left ({\ frac {\ partial ^ {2} u} {\ partial r ^ {2}}} + {\ frac { 1} {r}} {\ frac {\ partial u} {\ partial r}} \ right).

Мы будем искать решения в разделенных переменных, u (r, t) = R (r) T (t). {\ displaystyle u (r, t) = R (r) T (t).}u(r,t)=R(r)T(t).Подставив это в уравнение выше и разделив обе стороны на c 2 R (r) T (t) {\ displaystyle c ^ {2} R (r) T (t)}c ^ {2} R (r) T (t) дает

T ″ (t) c 2 T (t) = 1 R (r) (R ″ (r) + 1 r R ′ (r)). {\ displaystyle {\ frac {T '' (t)} {c ^ {2} T (t)}} = {\ frac {1} {R (r)}} \ left (R '' (r) + {\ frac {1} {r}} R '(r) \ right).}{\frac {T''(t)}{c^{2}T(t)}}={\frac {1}{R(r)}}\left(R''(r)+{\frac {1}{r}}R'(r)\right).

Левая часть этого равенства не зависит от r, {\ displaystyle r,}r, и правая часть не зависит от t, {\ displaystyle t,}t, , следовательно, обе части должны быть равны некоторой константе K. {\ displaystyle K.}K. Мы получаем отдельные уравнения для T (t) {\ displaystyle T (t)}T (t) и R (r) {\ displaystyle R (r)}R (r) :

T ″ (t) = K c 2 T (t) {\ displaystyle T '' (t) = Kc ^ {2} T (t) \,}T''(t)=Kc^{2}T(t)\,
r R ″ (r) + R ′ (r) - К р R (r) = 0. {\ displaystyle rR '' (r) + R '(r) -KrR (r) = 0. \,}rR''(r)+R'(r)-KrR(r)=0.\,

Уравнение для T (t) {\ displaystyle T (t)}T (t) имеет решения, которые экспоненциально растут или затухают для K>0, {\ displaystyle K>0,}K>0, линейны или постоянны для K = 0 {\ displaystyle K = 0}K = 0 и являются периодическими для K < 0 {\displaystyle K<0}K <0 . Физически ожидается, что решение проблемы вибрирующей головки барабана будет колебательным во времени, и при этом останется только в третьем случае K < 0, {\displaystyle K<0,}K <0, , поэтому для удобства мы выбираем K = - λ 2 {\ displaystyle K = - \ lambda ^ {2}}K = - \ lambda ^ {2} . Затем T (t) {\ displaystyle T (t)}T (t) - это лин Комбинация функций синуса и косинуса,

T (t) = A cos ⁡ c λ t + B sin ⁡ c λ t. {\ displaystyle T (t) = A \ cos c \ lambda t + B \ sin c \ lambda t. \,}T (t) = A \ cos c \ lambda t + B \ sin c \ lambda t. \,

Переходя к уравнению для R (r), {\ displaystyle R (r),}R(r),с замечанием, что K = - λ 2, {\ displaystyle K = - \ lambda ^ {2},}K = - \ lambda ^ {2}, все решения этого дифференциального уравнения второго порядка являются линейной комбинацией функций Бесселя порядка 0, поскольку это частный случай дифференциального уравнения Бесселя :

R (r) = c 1 J 0 (λ r) + c 2 Y 0 (λ r). {\ displaystyle R (r) = c_ {1} J_ {0} (\ lambda r) + c_ {2} Y_ {0} (\ lambda r). \,}R (r) = c_ {1} J_ {0} (\ lambda r) + c_ {2} Y_ {0} (\ lambda r). \,

Функция Бесселя Y 0 {\ displaystyle Y_ {0}}Y_ {0} не ограничен для r → 0, {\ displaystyle r \ to 0,}r \ to 0, , что приводит к нефизическому решению вибрирующей головки барабана проблема, поэтому константа c 2 {\ displaystyle c_ {2}}c_ {2} должна быть нулевой. Мы также предположим, что c 1 = 1, {\ displaystyle c_ {1} = 1,}c_ {1} = 1, , поскольку в противном случае эта константа может быть позже преобразована в константы A {\ displaystyle A}A и B {\ displaystyle B}B из T (t). {\ displaystyle T (t).}T (t). Отсюда следует, что

R (r) = J 0 (λ r). {\ displaystyle R (r) = J_ {0} (\ lambda r).}R (r) = J_ {0} (\ lambda r).

Требование, чтобы высота u {\ displaystyle u}u была равна нулю на границе головки барабана приводит к условию

R (a) = J 0 (λ a) = 0. {\ displaystyle R (a) = J_ {0} (\ lambda a) = 0.}R (a) = J_ {0} (\ lambda a) = 0.

Функция Бесселя J 0 {\ displaystyle J_ {0}}J_ {0} имеет бесконечное количество положительных корней,

0 < α 01 < α 02 < ⋯ {\displaystyle 0<\alpha _{01}<\alpha _{02}<\cdots }0 <\ альфа _ {{01}} <\ alpha _ {{02}} <\ cdots

Мы получаем, что λ a = α 0 n, {\ displaystyle \ lambda a = \ alpha _ {0n},}\ lambda a = \ alpha _ {{0n}}, для n = 1, 2,…, {\ displaystyle n = 1,2, \ dots,}n = 1,2, \ dots, поэтому

R (r) = J 0 (α 0 nar). {\ displaystyle R (r) = J_ {0} \ left ({\ frac {\ alpha _ {0n}} {a}} r \ right).}R (r) = J_ {0} \ left ({\ frac {\ alpha _ {{0n}}}} {a}} r \ r ight).

Следовательно, осесимметричные решения u {\ displaystyle u}u проблемы вибрирующей головки барабана, которая может быть представлена ​​в виде отдельных переменных:

u 0 n (r, t) = (A cos ⁡ c λ 0 nt + B sin ⁡ c λ 0 nt) J 0 (λ 0 nr) для n = 1, 2,…, {\ displaystyle u_ {0n} (r, t) = \ left (A \ cos c \ lambda _ {0n} t + B \ sin c \ lambda _ {0n} t \ right) J_ {0} \ left (\ lambda _ {0n} r \ right) {\ text {for}} n = 1,2, \ dots, \,}u _ {{0n}} (r, t) = \ left (A \ cos c \ lambda _ {{0n}} t + B \ sin c \ lambda _ {{0n}} t \ right) J_ {0} \ left (\ lambda _ {{0n}} r \ right) {\ text {for}} n = 1,2, \ dots, \,

где λ 0 n = α 0 n / a. {\ displaystyle \ lambda _ {0n} = \ alpha _ {0n} / a.}\ lambda _ {{0n}} = \ alpha _ {{0n} } / a.

Общий случай

Общий случай, когда u {\ displaystyle u}u также может зависеть от угла θ, {\ displaystyle \ theta,}\ theta, обрабатывается аналогично. Мы предполагаем решение в разделенных переменных,

u (r, θ, t) = R (r) Θ (θ) T (t). {\ displaystyle u (r, \ theta, t) = R (r) \ Theta (\ theta) T (t). \,}u (r, \ тета, t) = R (r) \ Theta (\ theta) T (t). \,

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

T ″ (T) c 2 T (t) = R ″ (r) R (r) + R ′ (r) r R (r) + Θ ″ (θ) r 2 Θ (θ) = K {\ displaystyle {\ гидроразрыв {T '' (t)} {c ^ {2} T (t)}} = {\ frac {R '' (r)} {R (r)}} + {\ frac {R '(r) } {rR (r)}} + {\ frac {\ Theta '' (\ theta)} {r ^ {2} \ Theta (\ theta)}} = K}{\frac {T''(t)}{c^{2}T(t)}}={\frac {R''(r)}{R(r)}}+{\frac {R'(r)}{rR(r)}}+{\frac {\Theta ''(\theta)}{r^{2}\Theta (\theta)}}=K

где K {\ displaystyle K}K - постоянная величина. Как и раньше, из уравнения для T (t) {\ displaystyle T (t)}T (t) следует, что K = - λ 2 {\ displaystyle K = - \ lambda ^ {2 }}K = - \ lambda ^ {2} с λ>0 {\ displaystyle \ lambda>0}\lambda>0 и

T (t) = A cos ⁡ c λ t + B sin ⁡ c λ t. {\ displaystyle T (t) = A \ cos c \ lambda t + B \ sin c \ lambda t. \,}T (t) = A \ cos c \ lambda t + B \ sin c \ lambda t. \,

Из уравнения

R ″ (r) R (r) + R ′ (r) r R (r) + Θ ″ (θ) р 2 Θ (θ) = - λ 2 {\ displaystyle {\ frac {R '' (r)} {R (r)}} + {\ frac {R '(r)} {rR (r)}} + {\ frac {\ Theta '' (\ theta)} {r ^ {2} \ Theta (\ theta)}} = - \ lambda ^ {2}}{\frac {R''(r)}{R(r)}}+{\frac {R'(r)}{rR(r)}}+{\frac {\Theta ''(\theta)}{r^{2}\Theta (\theta)}}=-\lambda ^{2}

получаем, умножая обе стороны на r 2 {\ displaystyle r ^ {2}}r ^ {2} и разделяя переменные,

λ 2 r 2 + r 2 R ″ (r) R (r) + r R '(Г) р (г) знак равно L {\ Displaystyle \ лямбда ^ {2} г ^ {2} + {\ гидроразрыва {г ^ {2} R' '(г)} {R (г)}} + { \ frac {rR '(r)} {R (r)}} = L}\lambda ^{2}r^{2}+{\frac {r^{2}R''(r)}{R(r)}}+{\frac {rR'(r)}{R(r)}}=L

и

- Θ ″ (θ) Θ (θ) = L, {\ displa ystyle - {\ frac {\ Theta '' (\ theta)} {\ Theta (\ theta)}} = L,}-{\frac {\Theta ''(\theta)}{\Theta (\theta)}}=L,

для некоторой константы L. {\ displaystyle L.}L. Поскольку Θ (θ) {\ displaystyle \ Theta (\ theta)}\ Theta (\ theta) периодический, с периодом 2 π, {\ displaystyle 2 \ pi,}2 \ pi, θ {\ displaystyle \ theta}\ theta , будучи угловой переменной, следует, что

Θ (θ) = C cos ⁡ m θ + D sin ⁡ m θ, { \ displaystyle \ Theta (\ theta) = C \ cos m \ theta + D \ sin m \ theta, \,}\ Theta (\ theta) = C \ cos m \ theta + D \ sin m \ theta, \,

где m = 0, 1,… {\ displaystyle m = 0,1, \ точки}m = 0,1, \ dots и C {\ displaystyle C}C и D {\ displaystyle D}D - некоторые константы. Это также подразумевает L = м 2. {\ displaystyle L = m ^ {2}.}L = m ^ {2}.

Возвращаясь к уравнению для R (r), {\ displaystyle R (r),}R(r),, его решение представляет собой линейную комбинацию из функций Бесселя J m {\ displaystyle J_ {m}}J_ {m} и Y m. {\ displaystyle Y_ {m}.}Y_ {m}. Используя аргументы, аналогичные аргументам в предыдущем разделе, мы приходим к

R (r) = J m (λ mnr), {\ displaystyle R (r) = J_ {m} (\ lambda _ {mn} r), \,}R (r) = J_ {m} (\ lambda _ {{mn}} r), \, m = 0, 1,…, {\ displaystyle m = 0,1, \ dots,}m = 0,1, \ dots, n = 1, 2,…, {\ Displaystyle n = 1,2, \ dots,}n = 1,2, \ dots,

где λ mn = α mn / a, {\ displaystyle \ lambda _ {mn} = \ alpha _ {mn} / a, }\ lambda _ {{mn}} = \ alpha _ {{mn}} /a,с α mn {\ displaystyle \ alpha _ {mn}}\ alpha _ {{mn}} n {\ displaystyle n}n -й положительный корень Дж м. {\ displaystyle J_ {m}.}J_ {m}.

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

umn (r, θ, t) = (A cos ⁡ c λ mnt + В грех ⁡ с λ mnt) J м (λ mnr) (C соз ⁡ м θ + D грех ⁡ м θ) {\ displaystyle u_ {mn} (r, \ theta, t) = \ left (A \ cos c \ лямбда _ {mn} t + B \ sin c \ lambda _ {mn} t \ right) J_ {m} \ left (\ lambda _ {mn} r \ right) (C \ cos m \ theta + D \ sin m \ theta)}u _ {{mn}} (r, \ theta, t) = \ left (A \ cos c \ lambda _ { {mn}} t + B \ sin c \ lambda _ {{mn}} t \ right) J_ {m} \ left (\ lambda _ {{mn}} r \ right) (C \ cos m \ theta + D \ sin m \ theta)

для m = 0, 1,…, n = 1, 2,… {\ displaystyle m = 0,1, \ dots, n = 1,2, \ dots}m = 0,1, \ dots, n = 1,2, \ dots

Анимации нескольких режимов вибрации

Ряд режимов показан ниже вместе с их квантовыми числами. Также указаны аналогичные волновые функции атома водорода, а также соответствующие угловые частоты ω = λ mnc = α mnc / a {\ displaystyle \ omega = \ lambda _ {mn} c = \ alpha _ {mn} c / a}{\ displaystyle \ omega = \ lambda _ {mn} c = \ alpha _ {mn} c / a} .

См. Также
Список литературы
  • H. Асмар, Нахле (2005). Уравнения в частных производных с рядами Фурье и краевые задачи. Верхняя Сэдл-Ривер, Нью-Джерси: Pearson Prentice Hall. п. 198. ISBN 0-13-148096-0. Cite имеет пустой неизвестный параметр: | coauthors =()
Последняя правка сделана 2021-06-18 12:14:03
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте