Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

ANALYTHICAL MODEL AND METHOD OF A DOME TYPE DYNAMIC TRANSDUCERS’ OSCILLATING SYSTEM CALCULATION

Vakhitov S.Y. 1 Basharin S.A. 1 Smirnova N.A. 1 Schitov I.N. 1
1 Saint-Petersburg state Institute of cinema and television
The article describes the distributed oscillating system in the form of circular plate clamped on the edge with a rigid central part. Such oscillating systems are often found in the form of membranes in the electroacoustic equipment, for example, a dome type electrodynamical transducer, and similar in instrumentation and mechanical engineering. A method of the substitution of such plates with the simple oscillating system is proposed. Analytical expressions for calculating the amount of the equivalent flexibility and equivalent corrugated surround area and aperture as a whole were obtained. It is shown that the previously proposed formulas are not suitable for the calculation of these parameters in some dynamic converters, because, firstly, these formulae are not consistent with the classical concepts of plate vibrations, and secondly, contrary numerous experiments. Derived analytical models provide a theoretical basis for the calculation and further development of a range of instruments and devices.
analytical model
the oscillating system
dynamic converter
1. Andreeva L.E. Uprugie jelementy priborov. M.: Mashinostroenie, 1981. 392 p.
2. Vahitov S.Y. Sovremennye mikrofony. Teorija, proektirovanie. SPb.: Izd. SPbGUKiT, 2003. 396 p.
3. Vakhitov J.S. Teoreticheskie osnovy elektroakustiki i elektro-akusticheskaja apparatura. -M.: Iskusstvo, 1982. 415 p.
4. Gutin L. J. Izbrannye trudy. L.: Sudostroenie, 1977. 600 p.
5. Lemb G. Dinamicheskaja teorija zvuka. M.: Gosizdat fizmat literatury, 1960. 372 p.
6. Strutt (Raylreigh) J.W. Teorija zvuka, t.1. M.: Gosizdat tehniko-teoreticheskoj literatury, 1955. 504 p.
7. Tjurin А.М., Stashkevich А.Р., Taranov E.S. Osnovy gidroakustiki. L.: Sudostroenie, 1966. 296 p.

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

В этом случае каждый колебательный элемент можно рассматривать как простую колебательную систему (с одной эквивалентной массой, гибкостью и трением), резонансная частота которой практически совпадает с первой собственной частотой замещаемой распределенной системы. К такому замещению обычно прибегают для упрощения расчетов электроакустических аппаратов, имеющих сложную акустико-механическую систему. Как правило, применяют энергетический метод нахождения эквивалентных параметров, впервые предложенный еще Рэлеем [6], подробнее с которым можно познакомиться в отечественной литературе, например в работе [3]. Здесь мы отметим только, что основными критериями эквивалентности в этом методе являются:

– сохранение величины кинетической (Eк), потенциальной (En) и диссипативной (Eq) энергии (последнее в целом ряде случаев не является столь существенным), что позволяет определить соответственно эквивалентную массу (mэ), гибкость (сэ) и активное сопротивление (rэ);

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

Любой электроакустический аппарат рассчитывается для работы в некотором диапазоне частот, поэтому, как правило, конструкцию аппарата выбирают такой, чтобы ее упругие элементы обладали не более чем одной собственной частотой (< f01) в заданном диапазоне, но в ряде случаев, например в электроакустической аппаратуре при симметрии колебательной системы, вполне реально, что устройство работает, как простая колебательная система и в диапазоне до второй симметричной моды [2], величину которой можно вычислить по приближенной формуле

vahitov01.wmf,

где Е и ρм – модуль упругости [н/ м2] и плотность материала упругого элемента; Rк – радиус кривизны купола (рис. 1).

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

vahitov02.wmf vahitov03.wmf.

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

О вычислении полной эквивалентной массы подвижной системы

Отметим здесь, что вычисление эквивалентной массы таких преобразователей, как правило, не вызывает больших трудностей. Полную эквивалентную массу подвижной системы, можно вычислить [2] по следующей формуле:

vahitov04.wmf (1)

где массу звуковой катушки можно найти по формуле

vahitov05.wmf (2)

где ρп – плотность материала провода; lп – его длина и dп – диаметр, k2 – технологический коэффициент, учитывающий массу изоляции провода и клеевого слоя, а также массу каркаса звуковой катушки (как правило, имеющего место только в случае излучателя звука). Его величина обычно лежит в пределах k2 ≈ 1,2–1,4.

pic_26.tif

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

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

vahitov06.wmf vahitov07.wmf (3)

где vahitov08.wmf vahitov09.wmf – площади купола и гофра.

Соколеблющуюся массу воздуха можно определить по известной формуле

vahitov10.wmf,

где плотность воздуха ρ0 = 1,23 кг/м3. Отметим, что масса звуковой катушки (2) обычно составляет в полной массе такой подвижной системы (1) около 70–90 %, поэтому даже в случае некоторых погрешностей при вычислении остальных составляющих полной массы обычно получаются в итоге весьма достоверные результаты. Резонансная частота ПС определяется по известной формуле

vahitov11.wmf

где с1 – эквивалентная гибкость диафрагмы. С её расчетом в подобных системах и возникают обычно сложности.

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

Определение эквивалентной гибкости и действующей площади диафрагмы

Впервые аналитическое соотношение для расчета упругости подвеса (воротника, гофра) диафрагмы динамического микрофона Sвр (рис. 2) было предложено Л.Я. Гутиным [4]. В принятых в нашей работе обозначениях (рис. 1, 2) это выражение для гибкости свр принимает вид

vahitov12.wmf (4)

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

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

pic_27.tif

Рис. 2. Физическая модель колебательной системы с жесткой центральной частью, закрепленной по периметру. Здесь ∆ – толщина упругого элемента [м]; Е и σ – модуль упругости [н/м2] и коэффициент Пуассона материала упругого элемента; а – внешний радиус закрепления [м]; b – радиус жёсткой части [м]; t = а – b [м]; γ = b/a

Известным способом замещения распределенных систем эквивалентными с сосредоточенными параметрами [3, 5] найдем эквивалентную гибкость круглой однородной пластины радиуса а и толщиной D. Смещение точек такой пластины на частотах, не превышающих первого собственного резонанса пластины (т.е. f < f01), выражается функцией

vahitov13.wmf (5)

где ξ0 – амплитуда смещения центра закрепленной по краю круглой пластины, r – радиус точки поверхности.

Считая, что площадь эквивалентного поршня Sэ = πа2, путем интегрирования этого выражения найдем объем вытесняемого пластиной воздуха V1 = ξ0·πa2/3 и величина ее эквивалентного смещения ξэ

vahitov14.wmf (6)

Согласно [6, 5] потенциальную энергию пластины в принятых нами обозначениях можно выразить как

vahitov15.wmf (7)

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

vahitov17.wmf (8)

Для гибкости, определяемой относительно центра пластины со смещением ξ0, в соответствии с [1, 5] получим

vahitov18.wmf (9)

Очевидно, что формула Гутина при ее корректности в пределе (при b1 → 0, t1 → а) в зависимости от площади приложения силы (центральная точка или равномерно распределенная по площади) и используемого смещения (ξ0 или ξэ) должна принимать вид либо (8), либо (9). Однако в указанных пределах она (4) преобразуется к виду

vahitov19.wmf (10)

Из сравнения формул (8)–(10) очевидно, что величина свр для однородной круглой пластины отличается соответственно в 16 и 16/3 раза. Это указывает на ошибочность полученного Гутиным аналитического выражения и, как следствие, неприменимости его для расчета гибкости диафрагм.

На основе описанной методики определим коэффициент гибкости «воротника», образованного периферийной частью пластины радиуса а, центральную часть которой занимает жесткий диск радиуса b (рис. 3, б), имеющий по всей площади одинаковое смещение ξв. Следовательно, в отличие от однородной пластины (рис. 3, а), переменное смещение ξ(r) реализуется на ширине «воротника» t = a – b. А при b, стремящемся к 0, модель (б) должна стремиться к модели (а).

pic_28.tif

а

pic_29.tif

б

Рис. 3. Характер колебаний круглой пластины, однородной (а) и с жестким центром (б), под действием приложенной внешней силы

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

vahitov20.wmf (11)

Это соотношение удовлетворяет граничным условиям

vahitov21.wmf vahitov22.wmf

Очевидно, что при равенстве воздействия величина смещения жесткой центральной части ξв будет меньше смещения центра ξ0 равновеликой однородной пластины, что зависит не только от упругих характеристик участка, лежащего в пределах b < r < a, но и от радиусов защемления этого участка. Величину ξв возможно определить, произведя процедуру нахождения эквивалентной гибкости однородной пластины радиусом t = a – b, аналогичную ранее проделанной (см. выражения (5)–(9)). В результате получим

vahitov23.wmf (12)

Откуда, взяв отношение (12) к (9), найдем относительное смещение центров таких пластин при р = const, т.е.

vahitov24.wmf (13)

где r = b/a.

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

vahitov25.wmf (14)

Производя интегрирование и алгебраические преобразования с учетом (13), получим

vahitov26.wmf (15)

Разделив полученное выражение на общую площадь пластины S2 = πa2, найдем среднее по площади такой пластины смещение xэр:

vahitov27.wmf (16)

Приравнивая потенциальные энергии исследуемой пластины и равновеликой однородной [6, 5], имеем

vahitov28.wmf

Откуда, после подстановки (16) и алгебраических преобразований определим величину гибкости пластины с жестким центром

vahitov29.wmf (17)

Или в более удобном для расчетов виде:

vahitov30.wmf (18)

Разделив величину вытесняемого пластиной объема воздуха (15) на смещение его центральной части (купола) (13), получим выражение для площади эквивалентного поршня:

vahitov31.wmf (19)

Здесь первое слагаемое – площадь купола, второе – эквивалентная площадь гофра (подвеса). Действительная площадь гофра

vahitov32.wmf

Очевидно, что соотношение эквивалентной и действительной (физической) площадей гофра

vahitov33.wmf

Далее была осуществлена проверка полученных в работе формул путем сопоставления данных эксперимента и расчета по конечным формулам (18), (19), а также массы по (1)–(3). После чего рассчитывалась частота механического резонанса подвижной системы нескольких моделей динамических преобразователей (микрофонов, купольных головок громкоговорителей). Затем, в соответствии с действующими стандартами (ГОСТ Р53575-2009 и ГОСТ Р53576-2009) проводились измерения этих резонансных частот на нескольких образцах преобразователей. Расхождение расчетных и экспериментальных результатов составляло не более 10 %, что, как правило, даже меньше допустимого разброса по этому параметру изделий из одной партии.

Выводы

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

В отличие от формулы Гутина, полученные нами выражения для гибкости пластины с жестким центром приводятся к (при ρ → 0) к классическим выражениям для гибкости однородной пластины. Таким образом, последние являются частным случаем от полученных нами более общих выражений. Полученные в работе аналитические модели и расчетные соотношения позволяют создать теоретическую базу для разработки методов проектирования различных типов купольных и ряда других преобразователей.

Рецензенты:

Майзель А.Б., д.т.н., профессор, начальник отдела, ОАО «ЦКБ МТ «Рубин», г. Санкт-Петербург;

Уваров В.К., д.т.н., профессор, ФГБОУ ВО «Санкт-Петербургский государственный институт кино и телевидения», г. Санкт-Петербург.