УДК 534.2
АНТОНОВ АЛЕКСАНДР ИВАНОВИЧ, канд. техн. наук, доцент, aiant58@yandex. ги
БАЦУНОВА АНАСТАСИЯ ВАЛЕРЬЕВНА, ассистент, nastybac@гambleг. ги
Тамбовский государственный технический университет, 392000, г. Тамбов, ул. Советская, д. 106,
ШУБИН ИГОРЬ ЛЮБИМОВИЧ, докт. техн. наук, профессор, niisf@niisf.ru
Научно-исследовательский институт строительной физики РААСН, 127238, г. Москва, Локомотивный пр., д. 21
РАСЧЕТ НЕСТАЦИОНАРНЫХ ЗВУКОВЫХ ПОЛЕЙ ПОМЕЩЕНИЙ ПРИ ЗЕРКАЛЬНО-ДИФФУЗНОЙ МОДЕЛИ ОТРАЖЕНИЯ ЗВУКА ОТ ОГРАЖДЕНИЙ*
Представлен уточненный метод расчета нестационарных звуковых полей помещений, разработанный на основе зеркально-рассеянной модели отражения звука от ограждений. Метод позволяет получить энергетические характеристики нестационарных звуковых полей помещений и дать более точную акустическую оценку помещений с позиции формирования и распространения в них отраженной звуковой энергии. Для реализации метода разработана компьютерная программа. Программа позволяет определять эквивалентные уровни звукового давления, максимальные и минимальные уровни, время реверберации, строить графики изменения в расчетных точках прямой, отраженной и суммарной звуковой энергии, а также решать ряд других исследовательских и практических задач.
Ключевые слова: нестационарные звуковые поля; зеркально-диффузная модель отражения звука; импульсный шум; защита от шума.
ALEKSANDR I. ANTONOV, PhD, A/Professor, aiant58@yandex. ru
ANASTASIYA V. BATSUNOVA, Instructor,
nastybac@rambler. ru
Tambov State Technical University,
106, Sovetskaya Str., 392000, Tambov, Russia,
IGOR'L. SHUBIN, DSc, Professor,
niisf@niisf.ru
Research Institute of Building Physics Raacs, 21, Lokomotivnyi Ave., 127238, Moscow, Russia
* Статья подготовлена в рамках выполнения НИР «Разработка методов оценки шумового режима в зданиях и на прилегающих к ним территориях для использования их при мониторинге шумового загрязнения среды и разработке мер по снижению шума в городской застройке» (код проекта 882) с финансированием из средств Минобрнауки России в рамках проектной части государственного задания в сфере научной деятельности.
© Антонов А.И., Бацунова А.В., Шубин И.Л., 2015
DIFFUSE-REFLECTION MODEL OF NON-STATIONARY ACOUSTIC FIELDS
The paper presents the estimated method of calculating non-stationary acoustic fields in rooms, based on the diffuse-reflection model of the acoustic reflection. This method allows obtaining the acoustic energy properties of non-stationary acoustic fields and provides the best estimate of room acoustics in terms of the formation and propagation of acoustic waves. The method implementation is assisted by the developed computer program which makes possible to detect the sound-pressure level, maximum and minimum levels, reverberation time, and construct diagrams of reflected and total acoustic energy at computational points.
Keywords: non-stationary acoustic field; diffuse-reflection model; impulsive noise; noise protection.
Нестационарные звуковые поля главным образом исследуют при анализе процессов реверберации, на основе которой рассчитывают средние коэффициенты звукопоглощения помещений или при оценке акустических характеристик помещений с точки зрения качественной передачи и восприятия речевой и музыкальной информации [1]. В то же время в практике борьбы с нестационарным шумом важно знать также подробные сведения о структуре шума в конкретных точках помещения (на рабочих местах), в том числе пиковые значения уровней звукового давления, соотношения уровней непостоянного и фонового постоянного шума, эквивалентные уровни шума и т. д. [2, 3]. Совершенствованию методов расчета нестационарных шумовых полей в аспекте борьбы с шумом в настоящее время не уделяется должного внимания, в частности, не разработана методика расчета нестационарных шумовых полей при сложном зеркально-диффузном характере отражения звука от ограждений.
В статье предлагается метод расчета нестационарных звуковых полей, учитывающий зеркально-диффузное отражение звука от ограждений.
Отражение звука от ограждений осуществляется по сложным пространственным зависимостям, определяемым формой поверхности, структурой материала, углом падения звука и другими параметрами. Характер отражения звука от поверхностей помещения является основным фактором, влияющим на процессы формирования отраженных звуковых полей помещений [3, 4]. Ранее в работах [5-7] установлено, что в большинстве случаев отражение звука от поверхностей носит сложный характер, при котором часть падающей на ограждение звуковой энергии отражается зеркально, а часть ее рассеивается по закону Ламберта. Такой характер отражения приводит к тому, что в структуре отраженного звукового поля помещения присутствуют две составляющие - зеркальная и рассеянная, которые различаются между собой по многим характеристикам. Они имеют разные длины свободного пробега звуковых волн [8], различные скорости затухания энергии, разные условия образования и переноса энергии и т. д. При каждом акте отражения часть энергии зеркальной составляющей переходит в рассеянную и пополняет ее. При этом обратный переход рассеянной энергии в зеркальную невозможен.
Таким образом, в случае зеркально-диффузного отражения звука от поверхности необходимо знать, кроме прямой энергии, величины зеркально и рассеяно отраженных энергий. В этом случае уровни звукового давления в расчетной точке помещения в любое время t будут определяться как
Ц = 101ё [ф^ +8? + 8р)10 ], (1)
где 10 - интенсивность на пороге слышимости; 8^, 8^, 8^- плотности звуковой энергии прямого звука, зеркальной и рассеянной составляющих отраженного звука в расчетной точке помещения в момент времени ^ с - скорость звука в воздухе.
Наиболее просто рассчитывается зеркальная составляющая звукового поля. Расчет ее при зеркально-диффузной модели отражения звука может производиться методом мнимых источников, прослеживания лучей и другими методами геометрической теории акустики [9]. При этом переход части зеркальной энергии в рассеянную при ее отражении учитывается уменьшением зеркальной энергии на сомножитель (1 - в), где в - коэффициент перехода зеркальной звуковой энергии в рассеянную. В каждом случае уменьшение акустической мощности мнимого источника за счет звукопоглощения и рассеивания происходит на величину (1 - а - в + ав). С учетом перехода зеркально отраженной энергии в рассеянную и поглощения ее на поверхности коэффициент потери звуковой энергии для зеркальной составляющей представляет собой величину, определяемую как
а ус = (1 -а)(1 -Р) = а + р-ар, (2)
где а и в - коэффициенты звукопоглощения и рассеяния звуковой энергии.
Тогда, в случае применения для расчета зеркальной составляющей метода мнимых источников, плотность зеркально отраженной звуковой энергии в момент времени t будет рассчитываться как
»р8-г,, „с - а „V, (3)
и 4кг,2
где Р - переменная во времени акустическая мощность источника звука; аус -
коэффициент, определяемый по формуле (2) (принят одинаковым для всех поверхностей); г, - расстояния от действительного (при , = 0) и мнимых (, > 0) источников до расчетной точки; N - количество учитываемых мнимых источников шума.
В помещениях других, отличных от параллелепипеда форм, выражения для расчета параметров мнимых источников имеют более сложный вид. В этих случаях, как правило, применяется метод прослеживания лучей.
Плотность рассеянно отраженной звуковой энергии в момент времени t может рассчитываться по формуле
1Шл-г/е)(1 -а, )с°*
8р = ^ (А ' л-г,/е)^ ™ (4)
5 """
где )(1 -а,)- интенсивность рассеянной звуковой энергии, излучаемой в расчетную точку помещения от поверхности ds'; г3 - расстояние от поверхности ds' до расчетной точки; 0 - угол направления от площадки на расчетную точку; а5 - коэффициент звукопоглощения поверхности ds'.
Величина ), как известно, при рассеянном отражении звука определяется интегральным уравнением Куттруфа [10], которое в случае переменной во времени акустической мощности имеет вид
где 1(^> t) - интенсивность падающей на участок йъ' звуковой энергии, определяемая вкладами отражений рассеянной энергии от всех других участков ограждений йъ; R - расстояния между элементами йъ и йъ'; и 9 2- углы между нормалями к соответствующим элементам йъ , йъ' и линией, их соединяющей; 101 - интенсивность прямого звука и звука от зеркальных отражений, участвующих в формировании отраженного рассеянного поля.
Интенсивность 101 определяется по формуле
Значения в формуле (6) соответствуют обозначениям в выражении (3).
В уравнении (5) учтена временная задержка ^ - R/c) поступления энергии, связанная с расстояниями R между элементами ограждений.
Расчет рассеянной составляющей на основе выражений (4) и (5) производится численным методом и имеет ряд трудностей. В первую очередь, это относится к тому, что источниками рассеянной энергии являются отражающие элементы помещения, излучающие в пространство помещения энергию после частичного рассеивания зеркальной составляющей поля при каждом акте отражения. Таким образом, параметры источников рассеянного звука могут быть определены только в процессе расчета зеркальной компоненты звукового поля. При этом для хранения численных значений характеристик источников рассеянного звука как нестационарных источников, а также для учета смещения во времени прихода нестационарной энергии от различных отражающих элементов требуются значительные вычислительные ресурсы.
Рассмотренная расчетная модель может быть использована при исследованиях влияния различных факторов на формирование нестационарных шумовых полей и решении практических задач борьбы с шумом. Для реализации представленной расчетной модели разработана компьютерная программа [11]. В результате вычислений по программе определяются эквивалентные уровни звукового давления для прямой, отраженной и суммарной звуковой энергии, возникающие при работе источника шума с непостоянной во времени звуковой мощностью, максимальные и минимальные уровни звукового давления, их разница, время реверберации, а также строятся графики изменения энергии во времени. Программа позволяет производить анализ влияния на распределение отраженной звуковой энергии коэффициентов а и р, устанавливать соотношение между собой энергий зеркальной и рассеянной составляющих в отраженном звуковом поле и их изменение во времени, исследовать время ревербера-
(5)
(6)
ции отдельно для зеркальной и рассеянной составляющих и в целом для всего отраженного звукового поля, а также решать ряд других задач.
Ниже в качестве примера на рисунке приведено окно программы с результатами расчета непостоянного во времени шума при действии в помещении размерами 14x4x3 м периодического источника импульсного шума. Периодичность действия источника 1,0 с, длительность импульса 0,4 с. Коэффициент звукопоглощения поверхностей а = 0,1, коэффициент рассеяния зеркально отраженной энергии в = 0,1.
* График
Координаты томки 6.00
Рисовать | КЗПср 0.10 1_ср 3,05 ВрЭйр 1,18 (11 30,54
4
3_
i
1 2 Дч \
л]
атраж зечке
|а«9
Эке. щэовми. гатимого эв 7150 Ьт1й*- и^а^лип ат и.ео
е^1гК1рмого эв |ЙС2П Вреив 1.£5
Загрузит* J
Сохранит ь дон ныв I ВыхДДиЗЬ:«! I
Результата: расчета изменения уровней звукового давления при работе в помещении нестационарного источника шума периодического действия:
1 - прямой; 2 - зеркально отраженный; 3 - рассеянно отраженный; 4 - суммарный уровень звукового давления
Из графиков видно, что рост рассеянной энергии при действии импульса источника шума отстает от роста зеркальной энергии. Скорость спадания зеркальной составляющей за счет совместного влияния поглощения звука и рассеивания значительно превышает скорость затухания рассеянной энергии. Более крутой начальный период затухания звука определяется затуханием ранних зеркальных отражений, завершающий этап затухания определяется затуханием рассеянной энергии, совпадающей с графиком затухания энергии по формуле Эйринга, используемой для оценки времени реверберации в помещениях с диффузными отраженными полями [2].
Аналогичные результаты и их анализ могут быть выполнены при других параметрах источника шума и других величинах а и р.
Таким образом, предложенная расчетная модель и разработанная на ее основе компьютерная программа представляют собой эффективный вычислительный инструмент, позволяющий анализировать процессы формирования звуковых полей в помещениях с нестационарным шумом и на этой основе производить выбор акустически эффективных и экономически оправданных строительно-акустических и других средств снижения шума.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Методика акустического проектирования одиночного пожарного оповещателя / В.Ф. Асминин, E.H. Епифанов, А.И. Антонов [и др.] // Шучный вестник ВГАСУ. Строительство и архитектура. - 2Gi3. - № 3 (3i). - С. i2i-i27.
2. Антонов, А.И. Метод расчета нестационарных шумовых полей в несоразмерных помещениях и помещениях сложных форм / А.И. Антонов, А.В. Бацунова, О.Б. Демин // Academia. Архитектура и строительство. - 2GiG. - С. i83-i85.
3. Бацунова, А.В. Расчет шумовых полей производственных помещений при работе источников шума периодического действия / А.В. Бацунова // Вопросы современной науки и практики. Университет им. В.И. Вернадского. - 2Gi5. - № 3 (57). - С. 46-52.
4. Антонов, А.И. Условия, определяющие процессы формирования шумового режима в замкнутых объемах, и их учет при оценке распределения звуковой энергии в помещениях / А.И. Антонов, А.В. Бацунова, И.Л. Шубин // Приволжский научный журнал. -ШГАСУ. - 2Gi5. - № 3. - С. 89-95.
5. Леденев, В.И. Расчет энергетических параметров шумовых полей в производственных помещениях сложной формы с технологическим оборудованием / В.И. Леденев, А.М. Макаров // кучный вестник ВГАСУ. - 2GG8. - № 2 (iG). - С. 94-iGi.
6. Антонов, А.И. Комбинированный метод расчета шумового режима в производственных зданиях теплоэлектроцентралей / А.И. Антонов, В.И. Леденев, Е.О. Соломатин // ^уч-ный вестник Воронежского государственного архитектурно-строительного университета. Строительство и архитектура. - 2Gii. - № 2. - С. i6-24.
7. Комбинированный метод расчета уровней шума в крупногабаритных газовоздушных каналах / В.П. Гусев, В.И. Леденев, М.А. Солодова [и др.] // Вестник МГСУ. - 2Gii. -№ 3 - i. - С. 33-38.
8. Антонов, А.И. Исследование длины среднего свободного пробега звуковых лучей в производственных помещениях с оборудованием / А.И. Антонов, В.И. Леденев, А.М. Макаров // Вестник Томского государственного архитектурно-строительного университета. -2GG7. - № 4. - С. ii5-i2i.
9. Осипов, Г.Л. Снижение шума в зданиях и жилых районах / Г.Л. Осипов. - М. : Стройиз-дат, i987.
iG. Kuttruff, H. Simulierte Nachhallkurven in Rechteckräumen mit diffusem Schallfeld /
H. Kuttruff // Acoustica. - i97i. - V. 25. - № 4. - P. 333-342. ii. Расчет нестационарных шумовых полей в помещениях сложных геометрических форм : свидетельство № 2Gi566i74i о регистрации программы для ЭВМ / Антонов А.И., Бацунова А.В., Леденев В.И. (РФ). - Опубл. G6.ii.2Gi5.
REFERENCES
1. Asminin V.F., Epifanov E.N., Antonov A.I., Kuznetsov S.N. Metodika akusticheskogo proek-tirovaniya odinochnogo pozharnogo opoveschatelya [Methods of acoustic design of a single fire detector]. Scientific Herald of the Voronezh State University of Architecture and Civil Engineering. Construction and Architecture. 2Gi3. No. 3. Pp. i2i-i27.(rus)
2. Antonov A.I., Batsunova A.B., Demin O.B. Metod rascheta nestatsionarnyih shumovyih poley v nesorazmernyih pomescheniyah i pomescheniyah slozhnyih form [The method of calculation of non-stationary noise fields in disproportionate and complex premises]. Academia. Arkhitektura i stroitel'stvo. 2GiG. Pp.i83-i85.(rus)
3. Batsunova A.B. Raschet shumovyih poley proizvodstvennyih pomescheniy pri rabote istochni-kov shuma periodicheskogo deystviya [The calculation of noise fields in production area under batch-type noise conditions]. Voprosy sovremennoi nauki i praktiki. Universitet im. V.I. Ver-nadskogo. 2015. No. 3. Pp. 46-52. (rus)
4. Antonov A.I., Batsunova A.B., Shubin I.L. Usloviya, opredelyayuschie protsessyi formirovani-ya shumovogo rezhima v zamknutyih ob'emah i ih uchet pri otsenke raspredeleniya zvukovoy energii v pomescheniyah [Noise conditions formed in confined space and assessment of acoustic energy distribution]. Privolzhsky Scientific Journal. 2015. No. 3. Pp. 89-95. (rus)
5. Ledenev V.I., Makarov A.M. Raschet energeticheskih parametrov shumovyih poley v proizvodstvennyih pomescheniyah slozhnoy formyi s tehnologicheskim oborudovaniem [Calculation of noise fields energy parameters in equipped industrial areas of irregular shape]. Scientific Herald of the Voronezh State University of Architecture and Civil Engineering. 2008. No. 2. Pp. 94-101. (rus)
6. Antonov A.I., Ledenev V.I., Solomatin E.O. Kombinirovannyi metod rascheta shumovogo rezhima v proizvodstvennykh zdaniyakh teploelektrotsentralei [Combined design method of noise conditions in industrial buildings of thermal power stations]. Scientific Herald of the Voronezh State University of Architecture and Civil Engineering. Construction and Architecture. 2011. No. 2. Pp. 16-24. (rus)
7. Gusev V.P., Ledenev V.I., Solodova M.A., Solomtin E.O. Kombinirovannyi metod rascheta urovnei shuma v krupnogaba-ritnykh gazovozdushnykh kanalakh [Combined design method of noise levels in large-size air-gas ducts]. Scientific and Technical Journal on Construction and Architecture. 2011. No. 3-1. Pp. 33-38. (rus)
8. Antonov A.I., Ledenev V.I., Makarov A.M. Issledovanie dliny srednego svobodnogo probega zvukovykh luchei v proizvodstvennykh pomeshcheniyakh s oborudovaniem [Investigation of mean free path of acoustic beams in equipped production areas]. Vestnik TSUAB. 2007. No. 4. Pp. 115-121. (rus)
9. Osipov G.L. Snizhenie shuma v zdaniyakh i zhilykh raionakh [Noise reduction in buildings and residential areas]. Moscow, 1987. (rus)
10. Kuttruff H. Simulierte Nachhallkurven in Rechteckräumen mit diffusem Schallfeld. Acoustica. 1971. V. 25. No. 4. Pp. 333-342.
11. Antonov A.I., Batsunova A.B., Ledenev V.I. Raschet nestatsionarnyih shumovyih poley v pomescheniyah slozhnyih geometricheskih form [Calculation of non-stationary noise fields in complex-geometry areas]. RF Certificate of State Registration of Software N 2015661741, publ. 06.11.2015. (rus)