|
|
Кинотеатральный громкоговорительТЕХНИЧЕСКОЕ ЗАДАНИЕ 1. Электрическая номинальная мощность одной головки, Wэ, Вт 250 2. Номинальный частотный диапазон, fн ? fв, Гц 20 ? 320 3. Допустимая неравномерность частотной характеристики осевого давления на нижней граничной частоте, ?Nн, дБ - 12 4. Полная акустическая мощность громкоговорителя, Wак., Вт 3,5 5. Среднее стандартное звуковое давление, рст., Па 0,4 6. Общий коэффициент потерь подвижной системы головки, ? 3 7. Объём ящика акустической системы громкоговорителя, Vящ., м3 - 8. Допустимая неравномерность характеристики направленности при = 90о, ?Nд, дБ (на частоте fв) - 4 9. Сопротивление звуковой катушки, Rк, Ом 16 10. Акустическое оформление фазоинвертор СОДЕРЖАНИЕ ВВЕДЕНИЕ 4 1. ПРЕДВАРИТЕЛЬНЫЕ РАСЧЁТЫ 6 1.1. Определение акустической мощности проектируемого громкоговорителя. 6 1.2. Выбор коэффициента расстройки и соотношения гибкостей. 7 1.3. Расчёт объёма и массы звуковой катушки. 8 2. РАСЧЁТ ЗВУКОВОЙ КАТУШКИ 10 2.1. Расчёт размеров и верхней граничной частоты диафрагмы. 10 2.2. Расчёт индуктивностей звуковой катушки и размеров каркаса. 14 3. РАСЧЁТ ДИАФРАГМЫ И ПАРАМЕТРОВ ГОЛОВКИ 18 3.1. Расчёт размеров и верхней граничной частоты диафрагмы. 18 3.2. Расчёт и построение характеристики направленности. 20 3.3. Выбор и расчёт размеров головки и подвесов. 21 3.4. Расчёт акустикомеханических параметров подвижной системы. 22 3.5. Расчёт электродинамических параметров преобразователя и эффективности головки. 25 4. РАСЧЁТ МАГНИТНОЙ СИСТЕМЫ 28 4.1. Выбор типа магнитной системы и исходных данных для её расчёта. 28 4.2. Предварительный расчёт размеров магнита и магнитной системы. 29 4.3. Уточнение проводимостей и корректирование размеров магнитной системы. 31 5. РАСЧЁТ АКУСТИЧЕСКОГО ОФОРМЛЕНИЯ И АЧХ 35 5.1. Расчёт числа головок громкоговорителя. 35 5.2. Расчёт фазоинвертора. 35 5.3. Расчёт параметров и вводных данных для машинного расчёта АЧХ громкоговорителя. 37 ЗАКЛЮЧЕНИЕ 44 ЛИТЕРАТУРА 45 ПРИЛОЖЕНИЯ 46 ВВЕДЕНИЕ Кинотеатральные громкоговорители представляют собой сложные акустические системы и выполняют роль конечного элемента в цепочке аппаратуры звуковоспроизведения. Громкоговоритель представляет собой агрегат, выполняющий преобразование электрических колебаний непосредственно в акустические. При помощи кинотеатральных громкоговорителей выполняется звуковое сопровождение театрализованных мероприятий и демонстраций кинофильмов, поэтому к ним предъявляются достаточно жёсткие и высокие технические требования качества и точности звуковоспроизведения: полоса воспроизводимых звуковых частот, равномерность амплитудно-частотной характеристики, неравномерность характеристики направленности, акустической мощности и т. д. Современные требования к качеству звукового сопровождения достигаются точным инженерным расчётом. Акустическая система подбирается под конкретную конфигурацию зала и строится на основе громкоговорителей, работающих в различной полосе частот слышимого человеком звукового диапазона. Обусловлено это невозможностью качественного воспроизведения всей полосы звуковых частот одной динамической головкой или одним громкоговорителем. Среди приёмов повышения качества звуковоспроизведения следует отметить многоканальность воспроизведения (стереофония, квадрофония), разбиение звукового диапазона (в высококачественной аппаратуре как минимум 20?20000 Гц) на два (НЧ и ВЧ) или три (НЧ, СЧ и ВЧ) поддиапазона и воспроизведения каждого из них через отдельные громкоговорители, размещение громкоговорителей в зале на различных позициях и с разной направленностью на основе специальных акустических расчётов. Кинотеатральный громкоговоритель подразделяется на два громкоговорителя - низкочастотный и высокочастотный. Цель настоящего курсового проекта - выполнить в соответствии с техническим заданием инженерный расчёт низкочастотного громкоговорителя, обеспечивающего воспроизведение звуковых частот в диапазоне 20?320 Гц. Предварительный анализ технического задания позволяет сделать вывод, что данный громкоговоритель будет использоваться в качестве низкочастотного звена акустической системы для зала большой вместимости. Конструктивно проект может быть реализован в виде звуковой колонки (удлинённого ящика с рядом электродинамических головок требуемой мощности). Акустическое оформление выполняется в виде ящика - фазоинвентора, для чего в его передней стенке делается отверстие с трубкой (горлом резанатора) расчётных размеров. Это позволит повысить эффективность излучения звука на нижней граничной частоте за счёт резонансного усиления. Электродинамическая головка состоит из подвижной системы - бумажной диафрагмы конической формы (диффузора), упруго прикреплённой на двух подвесах (гофре и центрирующей шайбе) к диффузодержателю. Диафрагма и механические колебания сообщаются подклеенной к ней звуковой катушкой, находящейся в сильном магнитном поле зазора магнитной системы. При выполнении проекта использовались учебное пособие, специализированная и справочная литература. 1. ПРЕДВАРИТЕЛЬНЫЕ РАСЧЁТЫ 1.1. Определение акустической мощности проектируемого громкоговорителя. Необходимая акустическая мощность громкоговорителя (ГГ) зависит от объёма кинозала, для которого он предназначен. Каждый канал звуковоспризведения (обычно их два) должен обладать акустичской мощностью, достаточной для обеспечения в середине зала уровня звука N = 90 дБ (или в конце зала 87 дБ). Порядок расчёта акустической мощности ГГ следующий. 1.1.1. Объём зала определяется его геометрическими размерами: , где l, b, h - соответственно длина, ширина и высота зала. Если зал проектируется, то сотношение сторон принимают, как правило, наилучшим: . Откуда можно вывести следующие зависимости: . 1.1.2. Требуемый объём зала диктуется его вместимостью, то есть количеством зрительских мест. Нормативный объём на одного зрителя составляет 4,5 м3. Общий объём зала получают как результат произведения нормативного объёма на число зрительских мест. Следует отметить, что зал с числом мест до 250 относят к категории малых, от 250 до 750 - средних и более 750 - больших залов. 1.1.3. Расстояние до центра зала полагаем равным половине его длины х = l/2. Громкоговоритель в данном случае работает на полупространство (сзади находится стена) и излучает низкочастотную сферическую волну. Тогда: , . В данном случае акустическая мощность громкоговорителя WАК = 3,5 Вт обусловлена техническим заданием на проект, поэтому необходимость её вычисления отпадает. 1.2. Выбор коэффициента расстройки и соотношения гибкостей. Вблизи нижней граничной частоты громкоговорителя име.т место два резонанса: воздушного объёма ящика-резонатора (фазоинвертора) на частоте fр; подвижной системы на частоте f0. Значения резонансных частот определяются формулами: , где С0 - скорость звука, м/с; Sот - площадь отверстия в стенке ящика, м2; Vящ. - объём ящика, м3; hэф. - эффективная длина трубки, вставляемой в отверстие; , где m1 и c1 - масса и гибкость подвижной системы. Частоту fр выбираем равной нижней граничной частоте громкоговорителя fн: Частота f0 превышает значение fр пропорционально коэфициенту расстройки ?0, который принимаем по таблице 1.1. учебного пособия [4] в зависимости от величины неравномерности частотной характеристики ?N из технического задания на проект. Для ?N = - 12 дБ находим ?0 = f0/fр = 2,52. По таблице 1.1. учебного пособия определяем также и предварительное значение соотношения gс гибкостей подвижной системы с1 и воздуха в ящике св, которое будет уточнено далее в процессе расчёта и оптимизации частотной характеристики громкоговорителя. Для ?N = - 12 дБ находим gс = с1/св = 0,8. 1.3. Расчёт объёма и массы звуковой катушки. Расчёт звуковой катушки (ЗК) имеет своей целью обеспечение тепловой надёжности (отсутствия перегрева) при возможно меньшей массе, что способствует повышению отдачи ГГ. Тепловая надёжность требует, чтобы на каждый ватт электрической мощности приходилось не меньше 1 см2 охлаждающей поверхности. Охлаждающая поверхность образуется суммой внутренней и внешней поверхностей катушки, то есть удвоенной площадью её внешней поверхности. Интенсивность охлаждения зависит не только от площади, а и от толщины ЗК. Для номинальной электрической мощности одной головки Wэ = 250 Вт на основании рекомендации [4] принимаем конструкцию звуковой катушки и двух слоёв намотки (n = 1). Задавшись числом слоёв намотки, рассчитываем по эмпирическим формулам необходимые велечины. Соотношение высоты катушки hк и высоты зазора магнитной цепи hз, то есть: , где W0 = 1 Вт - опорное значение мощности. Вероятная величина удельной охлаждающей поверхности: см2/Вт. где S0 = 2,5 см2/Вт, W1 = 104 Вт. Плотность тока: , где Rк - задаваемое значение сопротивления звуковой катушки, для кинотеатральных громкоговорителей принимаем Rк = 16 Ом; R0 = 1 Ом - опорное значение сопротивления; j0 = 70 А/мм2 - опорное значение плотности тока. А/мм2. Вычисляем объём провода: , где ?пр - удельное сопротивление материала провода, Ом?м. Для меди ?пр = 1,78?10-8 Ом?м. Предварительно приводим в соответствие размерность: j = 83,990 А/мм2 = 8,399?107 А/м2. м3. Масса звуковой катушки: , где ?пр - плотность материала провода, кг/м3. Для меди ?пр = 8,9?103 кг/м3. кг. По результатам расчётов составляем таблицу 1.1. Таблица 1.1. Значения величин, выбранных и вычисленных в главе 1 ?0 gc f0, Гц j, А/мм2 Sудохл, см2/Вт n gh* Vпр, м3 mк, г 2,52 0,8 50,4 83,990 1,533 1 1,748 1,991?10-6 17,720 2. РАСЧЁТ ЗВУКОВОЙ КАТУШКИ 2.1. Расчёт размеров и намоточных данных звуковой катушки. 2.1.1. Определяем диаметр голого провода ЗК: м. Уточняем полученный диаметр d по справочной таблице стандартных размеров обмоточного провода [1] и находим его диаметр с учётом толщины изоляции dиз. Получаем: d = 0,25 мм; dиз = 0,284 мм. По уточнённому значению d находим полную длину провода ЗК: м. 2.1.2. Вычисляем толщину bк звуковой катушки, площадь боковой поверхности Sк и действительное значение удельной охлаждающей поверхности: мм; м2, где кпр - коэффициент заполнения медью объёма проводника (учитывающий объём изоляции), допускается принимать по справочной таблице или вычислять по формуле: м2/Вт. Проверяем выполнение условия: , условие не выполняется. Уменьшаем плотность тока j и производим новый расчёт: j = 70 А/мм2; м3; кг; м; d = 0,28 мм; dиз = 0,315 мм; кпр = 0,790; м; Отклонение стандартного диаметра значительно, поэтому уточняем объём и массу провода: м3; кг мм; м2; , условие выполняется. 2.1.3. Ширина зазора магнитной системы ? складывается из толщины bк намотки ЗК, толщины каркаса ?1 = 0,1?0,15 мм, ширины двух просветов ?п между внутренней поверхностью ЗК и керном, наружной поверхностью и фланцем (рис. 2.1.). Ширина просветов увеличивается с увеличением объёма намотки и может быть найдена по формуле: м. Так как Vпр ? 0,2 см3 выбираем ?п = 1,5 мм. Находим ширину зазора: мм. 2.1.4. Определяем g? = hз/?, при этом необходимо учитывать, что при значениях gh < 3 увеличивается неравномерность распределения индукции по высоте зазора, что способствует возрастанию нелинейных искажений [3, 4]. Если gh < 10, то существенно усложняется центровка ЗК, так как минимальный перекос её оси вызывает касания к керну или фланцу. Оптимальную величину g? вычисляем по эмпирической формуле: , где Sк подставляется в см2, а Wэл - в ваттах. 2.1.5. Определение диаметра керна и высоты зазора. По найденному значению g? определяем высоту зазора магнитной цепи hз (она же является толщиной переднего фланца hпф): мм. Далее вычисляем соотношение gф между суммарным магнитным током и его полезной (основной) частью ?? и ?1 по эмпирической формуле: . Позже эта величина уточняется. Находим диаметр керна: мм. 2.1.6. Определение диаметра и высоты звуковой катушки и её намоточных данных. Средний диаметр звуковой катушки Dср и её радиус ак определяем по формулам: мм; мм. Вычисляем общее число витков ? : ; Находим высоту ЗК и фактическое значение gh: мм; . Отклонение от вычисленного ранее эмпирического значения допускается в пределах +20 % и - 10 %. В данном случае: 0,9gh* = 1,537 < gh = 2,053 < 1,2gh* = 2,097, условие выполняется. Проверяем активное сопротивление ЗК: Ом. Допустимое отклонение Rк в пределах ? 2,5 % от заданного значения (16 Ом). 0,975Rк = 15,6 Ом < Rк = 16,085 Ом < 1,025Rк =16,4 Ом , условие выполняется. Площадь сечения керна: см2. Найденные и вычисленные значения сводим в итоговую таблицу 2.1. Таблица 2.1. Сводные данные по результатам расчётов в п. 2.1 d, мм dиз, мм lп, м bк, мм Sк, см ?, мм g? hз, мм gф 0,28 0,315 50,35 0,315 136,9 3,6 5,325 19,2 1,802 Dкер, мм Dср, мм ак, мм ? ?сл1 ?сл2 hк, мм gh Sкер, см2 138 141,6 70,8 125 - - 39,4 2,053 157,48 2.2. Расчёт индуктивностей ЗК и размеров каркаса. 2.2.1. Определяем основную индуктивность катушки L1. Она образуется за счёт той части витков обмотки, которые находятся внутри зазора и, следовательно, создаваемые ими переменные магнитные поля замыкаются через стальные сердечники, роль которых выполняет керн и передний фланец. Выражаем их число через общее число витков ?, зная величину hк? на которую катушка выступает из зазора: мм, . Индуктивность L1 рассчитываем по формуле для индуктивности катушки со стальным сердечником: , где ?0 = 4??10-7 - магнитная постоянная вакуума и, приближённо, воздуха. Гн 2.2.2. Находим балластную индуктивность L2 части катушки, выступающей над керном, и той части витков внутри зазора, поля которых проходят по воздушному зазору, не попадая в стальной сердечник: , где . Гн. Так как L2 = 4,763?10-4 Гн < 0,1L1 = 7,996?10-4 Гн, принимаем L2 = 0,1L1 = 7,996?10-4 Гн [4]. Полная индуктивность: Гн. 2.2.3. Определяем размеры каркаса h1, h2 и h3. Размер каркаса h2 определяем по формуле: мм. Остальные размеры принимаем: h1 = (2,5?3,5) мм = 3 мм; h3 = (1,5?2,5) мм = 2 мм. Вычерчиваем эскиз профиля ЗК с указанием всех размеров в масштабе 5: 1 на рис. 2.2. По результатам расчётов составляем таблицу 2.2. Таблица 2.2. Сводные данные по результатам расчётов в п. 2.2 hк?, мм ?1 ?2 L1, мГн L2, мГн h1, мм h2, мм h3, мм 10,1 93 32 7,996 0,476 3 9 2 Параметры зазора магнитной системы Рис. 2.1 Профиль и размеры звуковой катушки М 2 : 1 Рис. 2.2 3. РАСЧЁТ ДИАФРАГМЫ И ПАРАМЕТРОВ ГОЛОВКИ 3.1. Расчёт размеров и верхней граничной частоты диафрагмы. 3.1.1. Установочный расчёт. Вычисляем исходные значения для поискового расчёта: начальные значения радиуса диафрагмы а, средней толщины диафрагмы ? и угла раскрытия конуса диафрагмы ?0. Выбор провода из стандартного ряда обусловил увеличение массы провода по сравнению с предварительно вычисленным в главе 1. Настоящая масса провода в соответствии с принятыми намоточными данными: кг . Определяем радиус диафрагмы: , где ?0 = 1,23 кг/м3 - плотность воздуха; ?0 = 9,0?102 кг/м3 - плотность бумаги диффузора; В = ?/а - отношение толщины диафрагмы к её радиусу, принимаем равным 3,5?10-3. м. Вычисляем число радиальных узловых линий: , где ?0 - угол раскрытия конуса диафрагмы, принимаем равным 55о; g - отношение радиуса ЗК к радиусу диафрагмы; ; . Округляем значение n до ближайшего целого : . Рассчитываем частоту колебаний диафрагмы на которой происходит резкий спад её излучения: , где Ед = 2,7?109 Н/м2 - модуль упругости бумаги диффузора. Число узловых окружностей к = 2. Получаем: Вычисленное значение fп2 не попадает в допустимый интервал: 0,9fВ = 288 Гц ? fп2 = 345,651 Гц ? fВ = 320 Гц, где fВ - верхняя граничная частота по техническому заданию. Необходимо проведение поисковых расчётов. 3.1.2. Поисковые расчёты. Поисковый расчёт имеет цель получить необходимое значение fп2 (удовлетворяющее условию 0,9fВ < fп2 < fВ). Для этого вычисления производят при различных значений ?0 в диапазоне от 60о до 47,5о для фиксированных значений В, а и g. Результаты расчёта представлены в таблице 3.1. Таблица 3.1 Результаты поисковых расчётов при В = 3,5?10-3 В = 3,5?10-3, а = 165 мм, g = 4,6?10-1 ?0о 60 57,5 55 52,5 50 47,5 16 15 15 14 13 13 fп2, Гц 318,678 335,016 344,575 354,808 364,766 365,339 Полученные значения удовлетворяют условию при максимальном угле раскрытия конуса диафрагмы. Принятое значение fп2 попадает в допустимый интервал: 0,9fВ = 288 Гц ? fп2 = 318,905 Гц ? fВ = 320 Гц, при а =165 мм и ?0о = 60,0о . Достоверность значений проверена с помощью специализированной программы на языке программирования высокого уровня QuickBASIC 4.5, листинг которой приведён в приложении 1. Контрольная распечатка результатов работы программы совпадает с табличными значениями ручного расчёта. 3.1.3. Проверочный расчёт. Проверку результатов поискового расчёта производим путём определения действительной минимальности fп2. Для этого при принятых значениях В, g и ?0 вычисляем три значения частоты fn2: при , , . Результаты заносим в таблицу 3.6. Таблица 3.6 Проверка минимальности fn2 В = 5,5 ? 10-3, а = 165 мм, g =4,3?10-1 n 15 16 17 fп2, Гц 328,698 318,627 318,746 Найденное значение fn2 = 318,627 Гц действительно минимально при неизменности принятых значений остальных переменных расчётного уравнения. 3.2. Расчёт и построение характеристики направленности. Характеристика направленности выражается соотношениями: , где к = ?/с0 - волновое число; ? - угол между направлением рабочей оси диафрагмы и данным направлением; J1 - функция Бесселя первого порядка для числового значения аргумента каsin?. Расчёт производим для верхней граничной частоты, поскольку характеристика направленности любого излучателя сужается с повышением частоты. Наибольшее значение аргумента ка, соответствующее sin? = 1 и f = fв, будет равно: . Результаты вычислений характеристики направленности при значениях угла ? с шагом ?? = 10о в пределах 0 ? ? ? 90о записываем в таблицу 3.7. Таблица 3.7 Результаты расчёта характеристики направленности ?о кВаsin? R(?) 0о 0,000 1,000 10о 0,169 0,996 20о 0,334 0,986 30о 0,488 0,971 40о 0,627 0,952 50о 0,747 0,932 60о 0,845 0,913 70о 0,917 0,899 80о 0,961 0,889 90о 0,976 0,886 Неравномерность характеристики направленности: . Проверяем полученное значение на соответствие требованиям технического задания: дБ; -4,5 дБ ? ?N?=90о = - 1,051 ? -3,5 дБ. По результатам расчёта характеристики направленности строим диаграмму на рис. 3.1. 3.3. Выбор и расчёт размеров головки и подвесов. Коническая диафрагма подвешивается с помощью двух гофрированных воротников, подклеиваемых по двум контурам: основанию конуса (гофр) и несколько выше контура подклейки ЗК (центрирующая шайба). Диафрагма с подклеенной к ней ЗК образует подвижную систему (ПС) головки с общей массой m1 и общей гибкостью с1 подвесов: , где сгф - гибкость гофра; сцш - гибкость центрирующей шайбы. Поэтому: . Каждый из подвесов состоит из гофрированной части (шириной b1) и плоской части (шириной b2), зажимаемой кольцом (или кольцевыми сегментами) на соответствующей полочке диффузодержателя (ДД). Ряд размеров головки ГГ стандартизирован: габаритный размер D, превышение которого в процессе производства не допускается; диаметры d1 установочных отверстий винтов, закрепляющих прижимные кольца (или кольцевые сегменты); расстояние d2 между центрами установочных отверстий (т. е. диаметр окружности центров установочных отверстий). Примерный размер D: 2 (165 + 40) = 410 мм. По таблице стандартных размеров [4] подбираем ближайший стандартный размер D = 400 мм и ряд сопутствующих типоразмеров d1 = 5,5 мм, d2 = 375 мм. Определяем ширину плоской (зажимаемой) части гофра: 400 - 375 = 25 мм и самого гофра 0,5 ? 400 - 165 - 25 = 10 мм . Соответствующие размеры ЦШ подбираем из соотношений: (1,6 ? 2,0) ? 10 = 16 мм , 25 мм . 3.4. Расчёт акустикомеханических параметров подвижной системы. Полная масса m1 складывается из массы катушки mк, диафрагмы - mд, слабоколеблющегося с диафрагмой воздуха - mсв, эквивалентных масс гофра и ЦШ, принимаемых за 20 % от mд: . Диаграмма направленности головки громкоговорителя Рис. 3.1 Масса диафрагмы: , где ?0 - угол раскрытия конуса диффузора; ? = B ? a - толщина диффузора; ?д = 900 кг/м3 - плотность материала диффузора; S - площадь основания диффузора (эффективная площадь диафрагмы); ?0 - плотность воздуха. , кг/м3; ? = 3,5 ? 10-3 ? 16,5 ? 10-2 = 5,775 ? 10-4 мм ; ; . Гибкость общая: . Гибкость центрирующей шайбы принимаем: сцш = (4?5)сгф = 4сгф . Гибкость гофра: , откуда 1,25 ? 1,375 ?10-4 = 1,718 ? 10-4 мм/Н . Гибкость центрирующей шайбы: . Коэффициент потерь ПС: . Круговая резонансная частота головки: ?0 = 2?f0 = 2 ? 3,1415 ? 50,4 = 316,663 Гц . Активное сопротивление: . 3.5. Расчёт электродинамических параметров преобразователя и эффективности головки. К электродинамическим параметрам преобразователя относятся частота электромеханического резонанса fэм, индукция В магнитного поля в зазоре и коэфициент электромеханической связи М = Вl, где l - эффективная длина провода ЗК, то есть длина той части провода катушки, витки которой пронизываются достаточно сильным магнитным полем (находятся внутри зазора или в непосредственной близости от него). Поэтому l всегда меньше lп и их отношение gl = l/lп всегда меньше единицы. Вычисляем gl по эмпирической формуле: Эффективная длина: 50,35 / 1,709 = 29,462 м. Индукцию В вычисляем с помощью заданного в техническом задании коэффициента общих потерь подвижной системы ? = 3,0: 2 - 0,215 = 1,785 . Индукция: . Проверка: 0,6Тл ? В = 0,869 Тл ? 1,28Тл , условие выполняется. Частота электромеханического резонанса: 183,639 Гц. Результат округляем до ближайшего меньшего значения стандартного ряда частот. Частота fэм = 1,6 ? 102 Гц является опорной частотой. Вычисляем КПД на опорной частоте: Значение КПД в % : 3,278 ? 10-2 ? 100% = 3,28% . По результатам расчётов составляем таблицу 3.4. Таблица 3.4. Сводные данные по результатам расчётов в главе 3 ?0 g a, мм ?, мм D, мм d1, мм d2, мм 60о 4,3?10-1 165 5,775 ? 10-4 400 5,5 375 b1гф, мм b1цш, мм b2гф, мм mд, г m1, г с1, мм/н сгф, мм/н 10 16 25 28,250 72,540 1,375?10-4 1,718?10-4 сцш, мм/н ?1 r1, кг/с gl ?эм fэм, Гц ноп, % 6,872?10-4 0,215 4,941 1,709 1,785 200 3,28 На рис. 3.2 по вычисленным данным выполняем эскиз диафрагмы с подвесами с указанием основных размеров и построений. Диафрагма с подвесами M 1 : 2 1 - диффузор, 2 - гофр, 3 - центрирующая шайба Рис. 3.2 4. РАСЧЁТ МАГНИТНОЙ СИСТЕМЫ 4.1. Выбор типа магнитной системы и исходных данных для её расчёта. Выбор типа магнитной системы и материала магнита производим с учётом стоимости магнитного материала, мощности, размеров и стоимости проектируемой головки, технологических и конструктивных факторов. Стоимость магнитных сплавов (феррокобальтовых магнитов ФКМ) значительно выше стоимости оксидно-бариевых материалов (ОБМ). В головках сравнительно большой мощности при Wэл > 200 Вт по конструктивным соображениям целесообразно применение ФКМ, выпускаемых заводами в кольцевом варианте. Материалом для изготовления ФКМ служат специальные сплавы металлов (стали, никеля, кобальта с прочими добавками). Конструктивно кольцевой ФКМ представляет собой диск с малым диаметром и большой высотой, это обусловлено тем, что ОБМ по сравнению с ФКМ обладают меньшей остаточной индукцией и большей коэрцитивной силой. Применение ОБМ не позволит разместить звуковую катушку по высоте из-за малой высоты при неоправданно большом диметре (большем диаметра диафрагмы), применение кернового ФКМ при большой мощности конструктивно не оправдано. Таким образом, принимаем к расчёту магнитную систему на основе кольцевого ФКМ. Магнитный материал марки ЮН15ДК25БА. Строим кривую размагничивания на рис. 4.1, которая характеризует зависимость остаточной индукции данного материала (Тл) от напряжённости (коэрцитивной силы) в А/м. Магнитная энергия максимальна (при минимально возможном объёме магнита) в рабочей точке, отмеченной на кривой размагничивания. Определяем по рис. 4.1 параметры рабочей точки: В0 = 1,05 Тл; Н0 = 52 ? 103 А/м. Записываем исходные данные к расчёту: высота зазора hз = 19,2 мм; ширина зазора ? = 3,6 мм; индукция магнитного поля в зазоре В = 0,869 Тл; диаметр керна Dкер = 138 мм. 4.2. Предварительный расчёт размеров магнита и магнитной системы. Магнитный поток, создаваемый постоянным магнитом, направляется с помощью элементов магнитной системы - фланцев, керна и магнитного стакана - в зазор магнитной цепи, формируя основной поток Ф1. Но дополнительно возникает ряд паразитных потоков (потоков рассеяния) - Ф1, Ф3 и т. д., величины которых обратно пропорциональны воздушным промежуткам, через которые эти потоки проходят. Вычерчиваем расчётную схему магнитной системы на рис. 4.2, указывая зоны возникновения потоков рассеяния обозначаем порядковыми номерами 1?11. Величины паразитных потоков на этих участках оцениваются магнитными проводимостями gi, где i - порядковый номер участка. Из-за отсутствия некоторых размеров МС на начальном этапе расчёта значения некоторых проводимостей не могут быть вычисленны и достаточной степенью вероятности оцениваются по эмпирическому уравнению: , где g? = hз/? - отношение, найденное при расчёте ЗК; А - эмпирический коэффициент зависящий от типа МС. 4.2.1. Определяем высоту магнита: , где ?0 = 4??10-7 Гн/м - магнитная постоянная вакуума; Н0 - напряжённость в рабочей точке размагничивания; Q ? 0,8 - коэффициент, учитывающий влияние собственного магнитного сопротивления магнита, керна и фланцев. . Кривая размагничивания магнитного материала марки 28БА190 Рис. 4.1 Схема основного и паразитных магнитных потоков Рис. 4.2 4.2.2. Задаём расстояние ?. Для кольцевого ФКМ: . 4.2.3. Рассчитываем значения проводимостей для известных на данной стадии расчёта размеров МС: ; ; ; ; . Вероятное значение суммы недостающих проводимостей производим по эмпирической формуле, принимая значение А = 7,5 [4]: . Сумма всех проводимостей: . 4.2.4. Площадь сечения магнита: м2. 4.2.5. Внутренний диаметр магнита: м. Внешний диаметр магнита: м. При использовании ФКМ диаметры фланцев выбираются равными внешнему диаметру магнита: м. 4.3. Уточнение проводимостей и корректирование размеров МС. Определённая ранее по эмпирической формуле сумма недостающих проводимостей g? требует уточнения по теперь уже полному набору необходимых размеров МС: ; ; ; ; ; ; . Найденные значения выписываем в таблицу 4.1. и вычисляем уточнённое значение суммы . Таблица 4.1. Значения основной и паразитных магнитных проводимостей. Номера проводимостей g1 g2 g3 g4 g5 g6 g7 g8 g9 g10 g11 ?gi Численные значения, м 2,312 1,127?10-1 2,254?10-1 2,447?10-1 5,971?10-1 1,048?10-1 1,237?10-2 1,189?10-1 2,815?10-1 0 0 4,121 4.3.2. Значения индукции в зазоре: . Проверка: 1,05В = 0,912 Тл ? В* = 0,877 Тл ? 1,10В = 0,956 Тл. условие не выполнено. Корректируем значение площади магнита Sм на 5 % в необходимом направлении и пересчитываем зависящие от неё размеры до выполнения условия. Для повышения точности расчёта используем ЭВМ и вспомогательную программу расчёта (приложение 3). на основании контрольной распечатки (приложение 4) составляем из окончательных значений магнитных проводимостей таблицу 4.2, и, после расчёта толщины заднего фланца, сводную таблицу размеров магнитной системы (таблица 4.3 ). Таблица 4.2 Значения основной и паразитных магнитных проводимостей Номера проводимостей g1 g2 g3 g4 g5 g6 g7 g8 g9 g10 g11 ?gi Численные значения, м 2,308 1,127?10-1 2,254?10-1 2,463?10-1 5,957?10-1 1,048?10-1 1,232?10-2 1,120?10-1 2,831?10-1 0 0 4,2 Площадь магнита: S = 1,3 ? 10-2 м2 . Магнитная индукция в зазоре: B* = 0,920 Тл (условие выполняется). 4.3.3. Определяем толщину заднего фланца. Толщина заднего фланца вычисляется программой на основе отношение суммы всех магнитных потоков к основному: , по формуле: . Таблица 4.2 Размеры элементов магнитной системы. Dм1, мм Dм2, мм hм, мм hпф, мм hзф, мм DФл2, мм Dкер, мм ?, мм 210 250 59,8 19,2 26,2 250 138 3,6 По данным таблицы 4.2. вычерчиваем эскиз магнитной системы на рис. 4.3. Эскиз магнитной системы Рис. 4.3 5. РАСЧЁТ АКУСТИЧЕСКОГО ОФОРМЛЕНИЯ И АЧХ 5.1. Расчёт числа головок ГГ. Расчёт числа головок ГГ производим по двум показателям. По акустической мощности: , где WАК - акустическая мощность головки, Вт; WЗЛ - объём зала, м3; ноп - КПД на опорной частоте. По среднему стандартному давлению, обусловленному техническим заданием: , где - превышение роп над рср в зависимости от неравномерности АЧХ на НЧ, принимаем по справочной таблице [4] для ?Nн = 3 дБ равным 0,014; рг - звуковое давление на опорной частоте, создаваемое одной головкой на расстоянии одного метра, Па. Па; . Полученные значения Nг округляем до целого. Для дальнейшего расчёта принимаем большее число головок. Отсюда: 1. 5.2. Расчёт фазоинвертора (ящика-резонатора). 5.2.1. Определение удельного объёма. Исходные данные: m1 = 7,254?10-2 кг , с1 = 1,375?10-4 м/Н , f0 = 50,4 Гц, gс = 0,8, а = 0,165м, S = ?а2= 3,14?0,1652 = 8,553?10-2 м2, Nг = 1, fр= fн = 20 Гц. На основе найденного ранее значения св = с1/gс = Vуд/ЕвS2, определяем: м3 . Внутренний объём воздуха в ящике без учёта толщины стенок: м3 . 5.2.1. Определение конструктивных параметров фазоинвертора. Минимальная ширина передней стенки ящика: м. Минимальная высота ящика: м. Глубина ящика: м. Проверяем выполнение условия: lящ = 0,55 м > bящ = 0,48 м. Проверяем достаточность глубины ящика для размещения в нём головки с кольцевой ОБМ: lящ = 0,55 м ? 0,6D = 0,6 ? 0,4 = 0,24 м. Задаёмся значением длины горла фазоинвертора: 0,5 lящ = 0,12 ? hг = 0,15 ? 0,7lящ =0,17. Площадь отверстия: , где hэф ? 1,1hг = 1,1 ? 0,15=0,165 м. м2. Радиус отверстия: м. Проверяем выполнение условия: , . 5.2.3. Коэффициент трансформации фазоинвертора: . Базу взаимодействия излучателей определяем как расстояние между центром излучателя (или группы излучателей) и центром отверстия (рис. 5.1, 5.2): м. Вычерчиваем эквивалентную схему ГГ на рис 5.3. Параметры элементов схемы рассчитываем через резонансную частоту и характеристическое сопротивление воздушного объёма ящика: кг; ; кг/с. 5.3. Расчёт параметров и вводных данных для машинного расчёта АЧХ громкоговорителя. Вычисляем значения исходных данных, необходимых для расчёта ЧХ звукового давления на ЭВМ: ; ; ; ; Передний вид ящика М 1 : 5 Рис. 5.1 Боковой разрез ящика М 1 : 2 Рис. 5.2 Эквивалентная схема громкоговорителя Рис. 5.3 ; ; ; ; ; ; ; ; дБ. Вычисленные результаты записываем в виде таблицы 5.1. Таблица 5.1. Параметры, вводимые в ЭВМ для расчёта АЧХ. Обозначения вводимых величин В формулах hя крd gc Н ?эм ?1 Q1 Q2 ?0 кра ?эм1 G N В программе ЕY SKD GC H EM ET Q1 Q2 XO SKR XM G N Численные значения 0,1132 0,122 0,8 0,001430 1,9990 0,2151 6,280?10-2 3,741?10-3 2,52 6,098?10-2 10,0 5,149 12,0 При помощи расчёта с применением ЭВМ и специализированной программы (приложение 5) заполняем таблицу 5.2. осевых значений давлений на частотах, выбираемых через 1/3 октавы в пределах от fн до fв. В таблице приняты следующие обозначения: р? - звуковое давление, создаваемое одной головкой группы; В(?) - коэффициент повышения эффективности излучения одной головки (КПЭИ) при работе в группе, ?(?) - коэффициент осевой концентрации (КОК) группы, ?1(?) - КОК одного излучателя. Расчёт на ЭВМ позволяет заполнить девять строк таблицы. Остальные строки заполняются результатами ручного расчёта по следующим формулам: ; ; ; ; ; . Таблица 5.2 Частотная характеристика р и ?N 1 i 0 1 2 3 4 5 6 7 8 9 10 11 12 2 ?i 1,000 1,260 1,587 2,000 2,520 3,175 4,000 5,040 6,350 8,000 10,08 12,70 16,00 3 ка 0,061 0,077 0,097 0,122 0,154 0,194 0,244 0,307 0,387 0,488 0,615 0,774 0,976 4 ?из 0,118 0,016 0,010 0,010 0,013 0,018 0,026 0,039 0,059 0,091 0,139 0,211 0,315 5 Qр 0,067 0,067 0,067 0,067 0,067 0,067 0,067 0,067 0,067 0,068 0,066 0,044 0,019 6 х10 -4,139 -3,100 -2,228 -1,474 -0,800 -0,169 0,453 1,100 1,806 2,608 3,550 4,682 6,066 7 ? 3,686 5,101 3,286 2,749 2,500 2,320 2,126 1,859 1,454 0,799 -0,145 -0,401 0,370 8 ? -21,95 -3,489 -2,215 -1,426 -0,746 -0,111 0,519 1,178 1,901 2,727 3,694 4,803 6,132 9 р?, Па 0,079 0,105 0,130 0,166 0,225 0,297 0,379 0,462 0,523 0,547 0,519 0,491 0,470 10 В(?) 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 11 ?(?) 1,001 1,001 1,002 1,002 1,004 1,006 1,010 1,016 1,025 1,041 1,065 1,105 1,173 12 fi 20 25 32 40 50 64 80 101 127 160 202 254 320 13 р 0,079 0,105 0,130 0,166 0,225 0,298 0,381 0,466 0,529 0,558 0,536 0,516 0,509 14 ?N -14,07 -11,601 -9,762 -7,639 -4,990 -2,557 -0,418 1,327 2,431 2,891 2,542 2,212 2,097 По результатам расчётов строим график частотной зависимости давления р (рис 5.3) и график неравномерности АЧХ (рис. 5.4) в дБ. ЧХ осевого давления Рис. 5.3 ЧХ неравномерности давления Рис. 5.4 ЗАКЛЮЧЕНИЕ В курсовой работе выполнен проект низкочастотного кинотеатрального громкоговорителя. После предварительных расчётов были выбраны коэффициенты расстройки и соотношения гибкостей, вычислены объём провода и масса звуковой катушки. В результате расчёта звуковой катушки получены точные её размеры и намоточные данные, на основе которых в свою очередь вычислены индуктивность и размеры каркаса ЗК. Запроектирована конструкция диафрагмы электродинамической головки, обеспечивающей эффективную работу в заданном диапазоне воспроизводимых частот. Допустимая неравномерность характеристики направленности находится в пределах, обусловленных техническим заданием на проект. Определены все размеры и акустико-механические параметры подвижной системы, вычислены электродинамические показатели преобразователя и эффективности головки. Общий коэффициент потерь электродинамической головки ниже величины, обусловленной техническим заданием. В качестве оптимального типа магнитной системы выбрана кольцевая МС с ФКМ и вычислены все параметры и конструктивные размеры МС. Акустическое оформление громкоговорителя выполнено в виде ящика-резонатора (фазоинвертора) с одной электродинамической головкой - звуковой колонки. Полученный проект низкочастотного громкоговорителя полностью соответствует требованиям технического задания и обеспечивает все заданные характеристики, за исключением допустимой неравномерности частотной характеристики осевого давления на нижней граничной частоте, которая на - 2 дБ превышает заданную. Все трудоёмкие вычисления проводились и контролировались с помощью ЭВМ и специализированных программ, написанных на языке программирования высокого уровня. Это позволяет говорить о высокой точности вычислений в данном проекте. СПИСОК ЛИТЕРАТУРЫ Белорусов Н.И. Электрические кабели, провода и шнуры. Справочник. - М.: Энергоатомиздат, 1987. Вахитов Я.Ш. Теоретические основы электроакустики и электроакустическая аппаратура. - М.: Искусство, 1982. Вахитов Я.Ш. Электродинамические громкоговорители. Теория расчета. - Л.: ЛИКИ, 1983. Вахитов Я.Ш., Смирнова Н.А. Расчёт и проектирование кинотеатральных громкоговорителей. Выпуск 1. Низкочастотный громкоговоритель. Уч. пособие. - С.-П.: СПИКиТ, 1994. Вахитов Я.Ш., Смирнова Н.А. Электродинамические громкоговорители. Расчёт и проектирование. - Л.: ЛИКИ, 1985. Вахитов Я.Ш., Смирнова Н.А., Чесноков М.А. Расчёт электроакустических аппаратов и акустических систем. Методическая разработка по применению ВТ в курсовом и дипломном проектировании. Часть I. - Л.: ЛИКИ, 1988. Енохович А.С. Краткий справочник по физике. Издание 2-е, перераб. и доп. М.: Высшая школа, 1976. Приложение 1 Листинг программы расчёта на ЭВМ диафрагмы ГГ 10 CLS 20 PRINT "***РАСЧЁТ ДИАФРАГМЫ ГГ***" 30 PRINT "Введите исходные данные:" 40 INPUT "Верхняя частота ГГ fв="; FB 50 INPUT "Масса катушки mk="; MK 60 INPUT "Радиус катушки ГГ ak="; AK 70 INPUT "Плотность диффузора ГГ Pd="; PD 80 ED = 2.7 * 10 ^ 9: BP = .0035 90 F = 60 100 N = 1: FR = .0174536 * F 110 A = (MK * 2 / (3.14159 * PD * BP + 2.46)) ^ (1 / 3) 120 G = AK / A 130 U = (15.1875 * (1 + G) ^ 6 * (SIN(FR) ^ 4))/(BP ^ 2 * N ^ 4 * (1 - G) ^ 4) 140 X = .785398 * BP * SQR(ED / (3 * PD)) / (A * (1 + G) ^ 2) 150 Y = .164256 * N ^ 4 * (COS(FR) ^ 4) 160 F1 = X * SQR(Y + U) 170 IF N > 1 THEN 200 180 N = N + 1: F2 = F1: AP = A: FP = 57.29578 * FR 190 GOTO 130 200 IF F1 < F2 THEN 180 210 PRINT 220 PRINT "Bp ="; BP 230 PRINT "a ="; AP 240 PRINT "Fo ="; FP 250 PRINT "ni ="; N - 1 260 PRINT "fn2 ="; F2 270 IF F2 > .9 * FB THEN 330 280 F = F - 2: AP = 0: FP = 0: F2 = 0 290 IF F > 45 THEN 100 300 BP = BP + .0005 310 IF BP > .006 THEN 30 320 GOTO 90 330 STOP Приложение 2 Результаты расчёта на ЭВМ диафрагмы ГГ ***РАСЧЁТ ДИАФРАГМЫ ГГ*** Введите исходные данные: Верхняя частота ГГ fв=? 320 Масса катушки mk=? 0.02759 Радиус катушки ГГ ak=? 0.0708 Плотность диффузора ГГ Pd=? 900 Bp = .0035 a = .1646776 Fo = 60.00105 ni = 17 fn2 = 319.6639 Приложение 3 Листинг программы расчёта на ЭВМ магнита с кольцевой магнитной системой 10 CLS 12 LPRINT : LPRINT : LPRINT : LPRINT : LPRINT : LPRINT 15 LPRINT " РАСЧЁТ МАГНИТА С КОЛЬЦЕВОЙ МС" 20 REM INPUT "О ="; O 30 REM INPUT "hз ="; HZ 40 REM INPUT "Dкер ="; DK 50 REM INPUT "B ="; B 60 REM INPUT "Bo ="; B0 70 REM INPUT "Ho ="; H0 80 REM INPUT "Q ="; Q 90 REM INPUT "Hk ="; HK 100 REM INPUT "Go ="; G0 110 REM INPUT "A ="; A 120 REM INPUT "Z ="; Z 124 A = 10: Z = .0054 125 I = 1 130 HM = O * B / (.000001256637# * H0 * Q) 140 G1 = 3.1415 * DK * HZ / O 150 G2 = 3.1415 * DK * .25 160 G3 = 2 * G2 170 G5 = 2 * (DK + 2 * Z / 3.1415) * .4348 * LOG((2 * Z) / (3.1415 * O)) 180 GS = G1 * (A / (2 * G0 + 1) - .25) 190 GS1 = G1 + G2 + G3 + G4 + G5 + GS 200 SM = (O * B / B0) * GS1 210 DM1 = DK + 2 * Z 220 DM2 = SQR(DM1 ^ 2 + 1.2732 * SM) 230 DF2 = DM2 - .5 * HM 240 DF1 = DF2 250 G4 = DF2 * .4348 * LOG(2 + DK / O) 260 G6 = (4.398 * HM) / (.4348 * LOG(.6 * (1 + 2 * Z / DK) + .4 * DM2 / DK)) 270 G7 = (4 * Z / (3.1415 * HM)) * (DK + 2 * Z / 3.1415) 280 G8 = 1.571 * DF2 290 G9 = 3.9 * HM * DM2 / SQR(SM) 300 G10 = 5.6 * HZ * DM2 / SQR(SM) 310 G11 = (1.4 / SQR(SM)) * (DF2 ^ 2 - DF1 ^ 2) 320 GS2 = G1 + G2 + G3 + G4 + G5 + G6 + G7 + G8 + G9 + G10 + G11 327 IF I > 1 THEN GOTO 334 330 SM1 = .7854 * ((DF2 ^ 2 - DF1 ^ 2) + HM ^ 2 + 2 * (DF2 + DF1) * (HM + 2 * HZ)) 332 GOTO 340 334 SM1 = SM 340 B1 = (B0 * SM1) / (O * GS2) 350 IF B1 < 1.05 * B THEN GOTO 380 360 IF B1 > 1.1 * B THEN GOTO 390 370 GOTO 400 380 I = I + 1: SM = 1.05 * SM1 385 GOTO 220 390 I = I + 1: SM = .95 * SM1 395 GOTO 220 400 GF = GS2 / G1 410 HZF = HZ * GF * B1 / 1.2 405 LPRINT : LPRINT 420 LPRINT " ", " B ="; B1, "hзф ="; HZF 430 LPRINT " ", " O = "; O, "G1 ="; G1 440 LPRINT " ", " hз = "; HZ, "G2 ="; G2 450 LPRINT " ", " hм = "; HM, "G3 ="; G3 460 LPRINT " ", " GS1 = "; GS1, "G4 ="; G4 470 LPRINT " ", " Sм ="; SM1, "G5 ="; G5 480 LPRINT " ", "Dкер = "; DK, "G6 ="; G6 490 LPRINT " ", " DM1 = "; DM1, "G7 ="; G7 500 LPRINT " ", " DM2 ="; DM2, "G8 ="; G8 510 LPRINT " ", "Dфл1 = "; DF1, "G9 ="; G9 520 LPRINT " ", "Dфл2 ="; DF2, "G10 ="; G10 530 LPRINT " ", " hпф = "; HZ, "G11 ="; G11 550 STOP Приложение 4 Результаты расчёта на ЭВМ размеров магнита и параметров магнитной системы B = .920432 hзф = 2.623782E-02 O = .0036 G1 = 2.308531 hз = .01917 G2 = .112717 hм = 5.984379E-02 G3 = .225434 GS1 = 4.119372 G4 = .2462749 Sм = 1.299978E-02 G5 = .5957323 Dкер = .138 G6 = .1047966 DM1 = .21 G7 = .1232576 DM2 = .2462749 G8 = .1195004 Dфл1 = .2462749 G9 = .283128 Dфл2 = .2462749 G10 = 0 hпф = .01917 G11 = 0 Приложение 5 Листинг программы расчёта на ЭВМ частотной характеристики ГГ 10 CLS 20 INPUT "Nп ="; EY 30 INPUT "KpD ="; SKD 40 INPUT "gc ="; GC 50 INPUT "H ="; H 60 INPUT "Nн¬ ="; EM 70 INPUT "N1 ="; ET 80 INPUT "Q1 ="; Q1 90 INPUT "Q2 ="; Q2 100 INPUT "V0 ="; X0 110 INPUT "KpA ="; SKR 120 INPUT "Vн¬1 ="; XM 130 INPUT "G ="; G 140 INPUT "N ="; N 150 I = 0 160 V = 1.259921 ^ I 170 X = (V - 1) / V 180 Z = EY ^ 2 + X ^ 2 190 U = (GC * X0) / (V * Z) 200 KA = SKR * V 210 NZ = (H * V ^ 2 * ((V ^ 2 + (SKD + EY) ^ 2) / Z)) * (1 - KA ^ 2 / 6 + KA ^ 4 / 72 - KA ^ 6 / 1728) 220 Q = Q2 + (Q1 * XM ^ 4) / SQR(V ^ 8 - XM * V ^ 7 + XM ^ 8) 230 X1 = V / X0 - (1 + GC) * X0 / V 240 A = EM + ET - Q * V * X1 + (EY * Q + X / V) * U 250 B = X1 + Q * V * ET + (X * Q - EY / V) * U 260 P = G * SQR(NZ / (A ^ 2 + B ^ 2)) 270 IF I = 0 THEN P1 = P 280 DN = 20 * LG(P / P1) 290 PRINT 300 PRINT "i ="; I 310 PRINT "v ="; V 320 PRINT "ka ="; KA 330 PRINT "nиз ="; NZ 340 PRINT "Q ="; Q 350 PRINT "X10 ="; X1 360 PRINT "A ="; A 370 PRINT "B ="; B 380 PRINT "DN ="; DN 390 I = I + 1 400 IF I < N + 1 THEN 160 410 STOP 11 Работа на этой странице представлена для Вашего ознакомления в текстовом (сокращенном) виде. Для того, чтобы получить полностью оформленную работу в формате Word, со всеми сносками, таблицами, рисунками, графиками, приложениями и т.д., достаточно просто её СКАЧАТЬ.  |
|
Банк готовых работ
Форма заказа
Скачать работу
экономические  Copyright © refbank.ru 2005-2024
Все права на представленные на сайте материалы принадлежат refbank.ru. Перепечатка, копирование материалов без разрешения администрации сайта запрещено. |
|