Восстановление распределения импульсных сигналов в пространстве методами ближнепольной акустической голографии

УДК 534.6.08 OECD 01.03 А А

Восстановление распределения импульсных сигналов в пространстве методами ближнепольной акустической голографии

Костеев Д.А,1, Салин М.Б.2 1Ст. лаб.-исследователь, ИПФ РАН, Н, Новгород, ул. Ульянова, д. 46

2

Аннотация

Одно из направлений развития методов акустических измерений - это переход к измерению уровней шума на том расстоянии от объекта, где это удобно исходя из геометрии и соотношения сигнал-шум, с последующим расчетом уровня сигнала на интересующей дистанции, с применением методов акустической голографии. При этом амплитуда и фаза сигнала должна быть измерена в достаточном количестве точек. Большой интерес представляет измерение характеристик источников широкополосных сигналов. При решении подобной задачи, если речь идет об акустике, наиболее важными являются два параметра: форма временной реализации и средний уровень поля давления. Форму временной реализации нагляднее всего исследовать на примере импульсных сигналов. В настоящей работе описан соответствующий лабораторный эксперимент, проведенный в безэховой камере, методика восстановления дальнего поля, а так же результаты расчетов и сравнение с экспериментальными данными.

Ключевые слова: ближнепольная акустическая голография, дальнее поле, импульсные

сигналы.

Restoring the distribution of pulse signal in space using methods of near-ßeld

acoustic holography

Kosteev DA.1, Salin M.B.2

1 Assistant researcher, Institute of Applied Physics of the Russian Academy of Sciences,

Nizhny Novgorod, Russia

2

Nizhny Novgorod, Russia

Abstract

The development of acoustic measurement technology gives opportunity to measure noise levels at a convenient distance from the object, with respect to its the geometry and signal-to-noise ratio, and then to calculate signal levels at the desired distance. The methods of acoustic holography is widely used here. To do so, the amplitude and phase of the signal must be measured in a sufficient number of points. Of great interest is the measurement of the characteristics of broadband signal sources is of great interest. When solving such a problem, two parameters are most important: the form of temporary realization and the average level of the pressure held. The calculated waveform as a function of time is best illustrated by the example of pulse signals. In the present work, describe the corresponding laboratory experiment conducted in an anechoic chamber. The technique for reconstructing the far field is given, followed by the results of calculations and comparison with experimental data.

Key words: nearfield acoustic holography, far field, pulse signal.

*E-mail: d_kosteev@mail.ru

Введение

Применение методов ближнепольной голографии открывает возможности для измерения характеристик источников при сложной внешней обстановке. Например, ближнепольные методы находят применение, если необходимо измерить поле источника в дальней зоне, а это не представляется возможным, из-за ограниченных размеров лабораторной установки, на которой проводится эксперимент. Гидролокатор с большой апертурой можно поместить в безэховой бассейн, измерить характеристики генерируемого им сигнала в ближней зоне и рассчитать диаграмму направленности в зоне Фраунгофера [1] при том, что размеры бассейна не будут позволять проводить прямые измерения дальнего поля,

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

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

1. Экспериментальная установка

Аналогично эксперименту, описанному в [4], измерения, представленные в данной работе, проводились в безэховой камере. Схема установки приведена на Рис, 1, Установка состояла из громкоговорителя без корпуса и приемного микрофона. Измерения проводились в плоскости, где была натянута система тонких металлических тросов, позволяющих закреплять микрофон в узлах прямоугольной сетки. Микрофон, с сохранением ориентации в пространстве, поочередно устанавливался в узлы сетки, динамик оставался неподвижным. За начало координат (0, 0, 0) принято положение громкоговорителя, В каждой точке производилось воспроизведение одной и той же реализации импульсного сигнала, время начала воспроизведения и время начала записи были синхронизованы. После объединения всех записей, была получена общая картина поля в данной плоскости. Были сделаны записи при двух значениях высоты закрепления громкоговорителя, которые соответствовали ближней и дальней зоне, В первом случае динамик располагался на расстоянии гпеаг = 18 см от зоны измерения над ее центром и излучал последовательность импульсов, каждый из которых имел частоту 1 кГц, длительность 4 мс, время между импульсами составляло 5 с. Оценка значения параметра Френеля ^Т = \хпеаг/Б2 составлявт: Бг = 1,53 (здесь: А - длина волны (м), Б - размер апертуры громкоговорителя (м)). Зона измерения представляла собой прямоугольник со сторонами 192 см и 131 см. Расстояние между узлами по обеим координатам было переменным и в среднем составляло 12 см. Во втором случае источник располагался на расстоянии Zfar = 215 см от зоны измерения, что соответствует зоне Фраунгофера (Бг = 18,27). На Рис. 2 построен волновой фронт по линии датчиков у = — 84 , z = znear.

Рис. 1. Схема установки

Рис. 2. Временная развертка сигнала в нескольких точках вдоль линии на границе

измерительной зоны

2. Анализ результатов

2.1. Обработка сигналов планарной антенны

Для расчета поля воспользуемся рассуждениями, приведенными в [4, 5]. Принимая во внимание то, что размер измерительного участка плоскости в ближней зоне был взят достаточно большим, чтобы поле источника спадало по краям, и на участках вне ее не возникало точек стационарной фазы, и, используя функцию Грина для мягкой границы, запишем интеграл Кирхгофа-Гельмгольца [6]:

dS

S

в данной формуле и далее по тексту: I - мнимая единица, р - комплексная амплитуда давления [обычно Па, но в данной статье используются условные единицы - В], здесь рассматривается одна частота, переход к широкому спектру будет сделан ниже, к -волновой вектор [м-1], Б - площадка в плоскости ^ = znear), в которой производились измерения давления, II - набор векторов из узлов в ближней зоне (х,у^пеаг) до необходимой точки наблюдения [м], п - нормаль к плоскости измерения. Формула (1) также называется второй интегральной формулой Гюйгенса, Она определяет значение звукового давления в полупространстве по известной величине звукового давления на плоскости.

Если расстояние до точки наблюдения много больше других линейных размеров задачи, то формулу (1) можно преобразовать следующим образом:

где r - набор векторов из начала координат в узлы измерительной зоны (x,y,znear) [м], а - угол [рад] между нормалью к плоскости измерения и R, Ввиду сделанного приближения R ^ r теперь под R можно понимать радиус-вектор из центра плоскости до необходимой точки наблюдения в дальней зоне. Такой подход для расчета поля в дальней зоне был 1Ш313Ш1 FPK (Far Plane Kirehhoff).

Выражение (2) можно обобщить на случай широкополосных сигналов, и в частности для импульсных сигналов, проведя стандартную процедуру перехода от спектрального представления к временному [7], Тогда (2) принимает вид:

где с - скорость звука [м/с], Дз - площадь "ячейки" [м2], г^ - вектор в j-yю точку измерения в ближней зоне [м].

Значения, полученные по формуле (3), сравнивались с сигналами, записанными в данной точке в дальней зоне. На Рис, 3 приведены сигналы: рассчитанные на основе ближнепольных измерений (синяя кривая) и контрольные измерения (красная).

S

Рис. 3. Сравнение результатов прямых измерений в дальней зоне и расчетов на основе ближнепольных измерений. Вверху - х = 0 см, у = 0 см, £ = Zfаr = 215 см, внизу -х = 24 см, у = 24 см, £ = Zfar Метод «РРК». По оси ординат отложено акустическое

давление в относительных единицах

Полученные результаты демонстрируют, что данная методика, при своей простоте, позволяет с достаточной точностью восстанавливать форму сигнала и фазовый центр в дальней зоне по ближнепольным измерениям, что показано при сравнении с контрольными сигналами. Расхождение по амплитуде связано с неустойчивостью метода, т.к. интеграл в формуле (2) в общем случае может расходиться. Чтобы показать это примем участок Б бесконечно большим и подставим в качестве р - сферическую волну от точечного источника с координатами (0,0, £0). Для сокращения выкладок пусть а = 0, т.е. К сонаправлено с осью После произведенных подстановок и в пределе г ^ £0 выражение (2) можно записать в следующем виде:

(V r 1

J cos(2kr)dr + i J sin(2kr)dr

где А - длина волны [м],

Рассматривая зону интегрирования как набор колец А/4, можно заметить, что оба

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

/ А/8 3А/8+пА/4 А/4+пА/4 \

ik cos

p(R) «--

J cos(2kr)dr + £ / cos(2kr)dr + i £ / sin(2kr)dr

n A/8+nA/4 n nA/4 )

Следовательно, уравнение (2) содержит расходящийся интеграл.

Однако для источников с более высоким порядком мультиполя используемый интеграл будет сходящимся, И возвращаясь к практической задаче обработки результатов измерения, заметим, что вычисление по формуле (3) будет давать конечный результат ввиду конечного размера S и введения аподизирующего множителя, но увеличение измерительного участка плоскости может не приводить к ожидаемому результату повышения точности. Поэтому метод является работоспособным, и в ряде задач продемонстрированная точность считается достаточной,

2.2. Обработка сигналов линейной антенной

Далее рассмотрим методику, описанную в [5], и опробованную в [4], где она получила название Far Line Transfer (FLT), Данная методика позволяет восстановить поле в дальней зоне, по результатам измерения линейной антенной [7]. Здесь это имитируется выборкой для обработки данных с датчиков, расположенных на одной прямой, расположенной в плоскости источника и точки, в которой восстанавливается поле, в данном случае прямая -у = 0 z = 0,

Введем цилиндрическую систему координат с осью x, < - угол между осью z и осью приемной антенны [рад], радиус р = \Jу2 + z2 [м]. Акустическое поле может быть разложено по цилиндрическим волнам:

p

l m

где blm, clm - коэффициенты, Hl - функция Ханкеля /-го порядка, 1-го рода (задает расходящуюся волну при временной зависимости e-iUt), волновое чи ело кт =

= \Jk2 — kX2m [м-1], если k2 ^ kX, в противном случае кт = i\Jk2m — k2 и Hl

р

Основная идея предлагаемого метода заключается в устранении суммирования по l в формуле (4), Для этого необходимо указать следующее ограничение метода измерения,

1, источник должен быть вытянутым вдоль оси, параллельной оси линейной

x

2, длина антенны должна превосходить длину источника плюс два радиуса 1-й зоны Френеля,

3, требуется рассчитать поле в полу плоскости, образованной осью источника и осью антенны, здесь - полуплоскость у = у0 = 0 z > 0(< = 0),

4, источник не должен иметь слишком сложную диаграмму направленности в плоскости у, z. Разложение (4) поля источника можно ограничить по азимутальным индексам до l ^ lma^ и восстановленном распределении по x в дальней зоне можно будет ограничиться гармониками kxm ^ kx,max = л/k2 + K2miw где Kminznear > V/max + 1

Метод FLT был описан в [5]. Согласно указанной работе можно произвести расчет поля в дальней зоне в точке (R sin а, Уо, R cos а) следующим образом:

/ г>\ bm* N Ax

I^U Znear COS аа

p(a,R) =

Я V 2пс

где а - угол между направлением, перпендикулярным оси антенны, и вектором в точку в дальней зоне [рад], N - количество измерительных точек, коэффициент Ьт* находится по формулам (6):

bm* = / bmвгкшdu

1 N/2

bm ^^ h1 (xn)p(xn,y0,znear )e xm n

N

n=-N/2

где m* (индекс коэффициента bm*) выбирается так, чтобы kxm*/к = sin а. hi - функция Ханнинга, определяемая следующей формулой (7):

, t , 1 Л 2п(жп - xi + Дх/2) \

hi(xn) = ^ ( 1 - cos-^-„ , л„— ) .

W

xN — x1 + Ax

N

W находится из условия нормировки ^ h2(xn) = N.

n=l

Сигнал, рассчитанный с помощью уравнения (5), представлен на Рис. 4

Рис. 4. Восстановление поля в точке (0,0,zfаr) с использованием линейной антенны методом «ПТ». По оси ординат отложено акустическое давление в относительных

единицах

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

Заключение

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

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

Авторы выражают благодарность профессору Хуанцай Лу (Ниапса1 Ьи) за полезные обсуждения.

Список литературы

1. W.-Q, Jing, Y.-B, Zhang, С.-Х. Bi, Nearfield acoustic holography-based method for measuring the sensitivity of a particle velocity sensor // Acta Acustica, 2015, Vol, 101, №,5, p. 855-858.

2. J.-M. Attendu, A. Ross. Time domain nearfield acoustical holography without wraparound error and spectral leakage for forward propagation //J. Acoust. Soc. Am. 2017. Vol. 141, №.2, p. 1039

3. Г.Ю. Годзиашвили, Определение полей отражений по измерениям в ближней зоне // Тр. X всерос, Конф, «Прикладные технологии гидроакустики и гидрофизики», Спб,: «Наука», 2010. С. 386-389.

4. М.В. Salin, D.A. Kosteev, Examples of usage of nearfield acoustic holography methods for far field estimations: Part 1, CW signals // arXiv:1812,03826 (2018) (https://arxiv.org/ftp/arxiv/papers/1812/1812,03826.pdf),

5. М.Б. Емельянов, Б.М. Салин, М.Б. Салин, А.В. Циберев. Восстановление временной зависимости и параметров сигнала широкополосных протяженных акустических источников в дальней зоне. Часть 1. Методы восстановления и технические средства // Акует. журн. 2014. Т. 60, №5. С. 608-616.

6. Е. Скучик, Основы акустики. Т. 2. Мир. М,, 1976.

7. Б.М, Салин, В.И. Турчин. Голографическое воспроизведение волновых полей с произвольной зависимостью от времени // Л куп. журн. 1992. Т. 38. JV2 1. С. 150-155.

References

1. W.-Q. Jing, Y.-B. Zhang, C.-X. Bi, Nearfield acoustic holography-based method for measuring the sensitivity of a particle velocity sensor // Acta Acustica. 2015. Vol. 101, №.5. p. 855-858.

2. J.-M. Attendu, A. Ross. Time domain nearfield acoustical holography without wraparound error and spectral leakage for forward propagation //J. Acoust. Soc. Am. 2017. Vol. 141, №.2, p. 1039

3. G.Yu Godziashvili. Determination of reflection fields from measurements in the near zone // Proceedings of the 10th All-Russian Conference "Applied technologies of hvdroacoustics and hydrophysies", Spb .: "Science", 2010. p.386-389. [in Russain]

4. M.B, Salin, D.A. Kosteev. Examples of usage of nearfield acoustic holography methods for far field estimations: Part 1. CW signals // arXiv:1812.03826 (2018) (https://arxiv.org/ftp/arxiv/papers/1812/1812.03826.pdf).

5. M.B. Emelvanov , B.M. Salin, M.B. Salin, A.V. Tsiberev, Reconstruction of the time dependence and signal parameters of far-field extended wideband sources: Part 1.

Reconstruction techniques and technical instruments // Acoustical physics 2014 Vol, 60, №5, p. 608-616.

6. E. Skudrzyk. The foundation of acoustics. Basic mathematics and basic acoustics. Springer - Verlag, Wien, New York, 1971

7. B.M. Salin, V.l. Turchin. Holographic reconstruction of wave field with arbitary time dependence // Sov, Phvs, Acoust. 1992 Vol. 38, 77.