CC BY
43
УДК 535.412
И. В. Волков1, В. В. Вышинский1'2, А. О. Кисловский1
1 Московский физико-технический институт (государственный университет) 2Центральный аэрогидродинамический институт им. проф. Н. Е. Жуковского
Теоретические основы метода спекл-голографии в экспериментальной механике
Представлена математическая модель физического явления голографической и спекл-интерферометрии, объединяющая существующие теории геометрической и дифракционной оптики. Эта модель используется в экспериментальной механике при анализе пространственных перемещений и деформаций натурных конструкций. Открываются новые возможности исследования широкого спектра задач.
Ключевые слова: спекл-голография, спекл-структура, дифракционная и компьютерная оптика, дифракционные решетки, дифракционное гало, преобразование Фурье, оптическая фильтрация, оптическое дифференцирование, экспериментальная механика.
I. V. Volkov1, V. V. Vyshinsky1'2, A. O. Kislovskiy1
1 Moscow Institute of Physics and Technology (State University) 2 Central Aerohydrodynamic Institute (TsAGI)
Theoretical bases of the speckle-holography method in
experimental mechanics
The paper presents a mathematical model of the holographic and speckle-interferometry physical phenomena combining geometric and diffraction optics theories. The model is used in experimental mechanics for space displacement and deformation full-scale construction analysis. A new possibility to investigate a wide range of problems is offered.
Key words: speckleholography, specklestructure, diffracted and computer's optics, diffractive grids, diffraction halo, Fourier transformation, optical Alteration, optical differentiation, experimental mechanics.
1. Введение
Спекл-интерферометрия используется сравнительно недавно, но общепризнанной теории пока не существует. Однако результаты использования ее впечатляют и обещают широкое использование этого явления в технике. В частности, компьютерная спекл-интерферометрия [1, 2] уже используется как альтернатива мокрому голографическому бесконтактному процессу. Представленная в статье теория спекл-интерферометрии объединяет существующие модели голографической интерферометрии диффузно отражающих объектов при их деформировании. Использование представленной теории при анализе деформаций позволяет получать компоненты тензора пространственных перемещений и деформаций, используемые в механике деформированного тела. Обобщение голографических принципов в спекл-голографии и спекл-интерферометрии открывает новые возможности исследования в различных областях науки и техники.
2. Физика спекл-структур
При освещении диффузной поверхности когерентным светом лазера каждая точка ее рассеивает в окружающее пространство сложную, присущую только ей волну когерентного света. Эта волна несет информацию о структуре, материале, обработке его поверхности, шероховатости и других свойствах. В результате интерференционного сложения волн всех точек освещенной поверхности в пространстве формируется рассеянное когерентное сложное спекл-поле, несущее полную информацию о поверхности освещенного объекта. Явление останется тем же самым, если заменить глаз фотоаппаратом или видеокамерой. На изображении будут наблюдаться субъектные спеклы, которые определяются апертурой объектива. Чем больше апертура, тем тоньше структура спеклов, т.к. диаметр дифракционной картины, создаваемой объективом, убывает с увеличением его апертуры (отверстия). Чтобы получить спеклы, достаточно поместить фотопластинку на каком-либо расстоянии от объекта, и на ней будут зарегистрированы объектные спеклы. По аналогии с явлением дифракции можно назвать эти спеклы в ближней зоне спеклами Френеля или объектными спеклами, а в дальней зоне при сфокусированном изображении - спеклами Фраунгофера или субъектными спеклами. При регистрации голограммы спеклы будут наблюдаться и в восстановленном изображении, это затрудняет наблюдение объекта, его качество, однако эти спеклы можно использовать в экспериментальной механике при измерении пространственных перемещений и деформаций поверхности объекта. Природу образования спекл-структуры разные авторы представляют по-разному [2, 3]. Примеров объяснения дифракционной природы образования спекл-структуры много [5, 6, 7, 8]. Дифракция когерентного освещения на регулярной дифракционной решетке имеет строго определенную пространственную частоту, определяемую формулой Брэгга:
где Л - длина волны, Б - период решетки, а - угол дифракции.
При дифракции на шероховатостях диффузной поверхности образуется спектр пространственных частот, несущий информацию об объекте. Спекл-поле можно рассматривать как результат дифракции на совокупности решеток с изменяющимися пространственными частотами (периодами) их пространственной ориентацией и типами. Поэтому при анализе спекл-полей при оптической фильтрации может быть использован аппарат фурье-преобразования, обладающий математической полнотой функционального анализа и простотой его реализации в оптике. При использовании оптической схемы (рис. 1) [5] для регистрации спекл-голограмм сфокусированного изображения двойной апертурой пучки, опорные друг другу, создают несущую пространственную частоту по всему изображению. При оптической фильтрации дифракционных порядков спекл-голограмм по классической схеме (рис. 2) на фоне спеклов образуется регулярная решетка.
При интерференции по интенсивности в каждой точке голографического изображения регистрируется информация о пространственных частотах по всему полю объекта. Таким образом, в каждой точке геометрического изображения имеется множество опорных пространственных частот по полю объекта. Это объясняет наличие сплошного спектра в дифракционном гало в частотной плоскости при фильтрации спекл-голограмм.
При использовании дополнительного опорного пучка по полю спекл-голограммы образуется регулярная несущая решетка, которая и наблюдается на фоне спеклов. Наличие дополнительного опорного пучка увеличивает чувствительность к вибрациям, поэтому при внестендовой спекл-голографии наличие вибраций размывает эти полосы и ухудшает запись спекл-голограмм. Аналогично происходит в сендвич-голографии [8]. Данное свойство применяется в оптоэлектронике и в радиоастрономии.
Спекл-голограмма (Sg) рассматривается как свертка © аппаратной функции (А) используемой оптики и регистратора с объектной волной Бд = Б ©А. Можно представить спекл-голограмму как свертку различных видов сигналов в поле регистратора с его
(1)
/ / X
/ / / / / : А С
Рис. 1. Измерение деформаций с помощью изображающей системы двойной апертуры
X Ц и 1-2 х'
^^^^^^\ Хо ио <г
1 f г Г 1
12 3 4
Рис. 2. Классическая схема фильтрации: фурье-преобразования Ы, 3 - частотная плоскость линзы Ы, 4 - линза обратного фурье-преобразования Ь2, 5 - плоскость производной спекл-интерферограммы
аппаратной функцией. Это обобщение значительно расширяет области применения методов спекл-голографии и приводит к тензорному аппарату булевой алгебры, что может быть использовано в медицине и нейрофизиологии, а также в астрофизике.
3. Использование спеклов в экспериментальной механике
На рис. 3 представлены оптические схемы хода лучей когерентного источника при освещении диффузно отражающей поверхности объекта и получения интерференционной спекл-структуры. Слева показано образование объектной спекл-структуры, регистрируемой при получении голограмм во встречных пучках, справа - оптическая схема образования субъектной спекл-структуры, регистрируемой при получении голограмм сфокусированного изображения.
Рис. 3. Схемы хода лучей при образовании спекл-полей диффузно отражающей поверхности
Средние размеры получаемых спеклов 53 определяются первым интерференционным максимумом диска Эйри (2а), пределом разрешения оптической системы по Релею [2, 3]:
5? = 1, 22V, (2а)
\ /
5? = 1,22М-/, (2Ь)
3 Б
где Л - длина волны, Ь - расстояние между объектом и плоскостью регистрации, V -размер регистрируемого (освещенного) поля объекта, V - диаметр апертуры линзы (объектива), М - коэффициент увеличения, / - фокусное расстояние линзы (объектива).
Комплексная амплитуда рассеянного света в произвольной точке пространства равна сумме амплитуд волн по всем пространственным частотам, рассеянных каждой точкой поверхности, являющихся источником сферических волн.
Волновой фронт излучения рассеянного диффузно отражающим объектом, расположенным в плоскости (х,у), падает на регистрирующую плоскость (х',у') (рис. 3). Комплексную амплитуду а(х,у) п-го точечного излучателя, его координаты (хп,уп), п = 1...Ы, согласно принципу Гюйгенса-Френеля, можно представить в виде функции
а(хп,Уп) = 5(х — Хп,у - Уп), (3)
где 5(х — хп, у — уп) - дельта-функция Дирака, а(х, у) - непрерывная функция.
5(х — хп, у — уп) = Ц 5п(х — хп,у — уп)Ахйу = |
1 при (х = хп,у = уп); 0 при (х = хп,у = уп).
Тогда выражение для суммы пространственных волн вблизи объекта, представляющее распределение комплексных амплитуд, будет иметь вид
N
ио(х,у) = ^ 5(х — хп,у — уп )а(х,у)вхр(г^п), (4)
п=1
где фп - относительная фаза волны от п-го излучателя.
Поле в дальней зоне можно связать преобразованием Фурье:
N
и (х',у')= ! J N 5(х — хп,у — уп)а(х,у) • (хх' + уу')^ йхйу, (5)
_^^_^^ п=1
—оо — оо
где Л - длина волны когерентного излучения, Ь - расстояние между плоскостями (х,у), (х',у').
После преобразования выражение (5) приводится к виду
N 2п
и(х',у') = ^2 ап(хп,уп)ехр-[фп + (хпх' + упу')}. (6)
п=1
Выражение (6) есть сумма диффузно отражённых от объекта волн, регистрируемых спекл-голограммой во встречных пучках. Описанную этим выражением спекл-структуру рассеянного объектом излучения легко наблюдать в дальней зоне. Эта спекл-структура используется при выделении отдельных компонент тензора деформаций. Такими образом, спекл-интерферограмму можно рассматривать как суперпозицию двух световых полей, диффузно рассеянных двукратно-экспонированной спекл-голограммой. Полагая, что после экспонирования и обработки голограммы ее амплитудное пропускание пропорционально распределению интенсивности в плоскости спекл-голограммы I(х, у), запишем результат
двукратно-экспонированной голограммы объекта, испытывающего одновременно деформацию ех, £у и смещение х0 в плоскости объекта (х,у) [7, 13]:
Т(х, у) = ¿0 - к1(х, у)©[6(х, у) + 5(х - £хХ - Хо,у - £уу)}, (7)
где т(х, у) - комплексная функция амплитудного пропускания спекл-голограммы двойной экспозиции, ¿о - пропускание неэкспонированной фотопластинки, к - коэффициент пропорциональности, I(х, у) - зарегистрированное распределение интенсивности в плоскости спекл-голограммы, х,у - координаты в плоскости объекта, хо - величина смещения объекта как жесткого целого по оси х, £х,£у - относительные деформации удлинения объекта по осям х, у, © - операция свертки сигналов.
Таким образом, на спекл-голограмме регистрируются деформации и смещения объекта. При фильтрации таких голограмм наблюдается изменение пространственных частот в фурье-плоскости и формируется определённый спектр в виде дифракционного гало. Фильтрация внечастотной плоскости связана с выделением отдельных областей этого спектра.
Поле наблюдения при проведении оптической фильтрация на расстоянии I от фокальной плоскости [10] можно представить в виде
(V2 + V2)) ) • ( 1 + ехр( (^£х + V 2£у ^
цх0
где ц,и - координаты в плоскости фильтрации, f - фокусное расстояние линзы, Л - длина волны монохроматического плоского пучка.
Регистрируемое распределение интенсивности может быть получено возведением в квадрат выражения (8):
1н(V, V) ~ С082 ( (ц2£х + V2£у) + 1^хо) ).
(9)
Получаемая спекл-интерферограмма содержит компоненты деформаций и смещений. Преобразование Фурье волнового поля линзой запишется следующим образом:
и (С,п) =! J N ио(х,у)вхр(-12п • (£х + пу)^ Лхйу, (10)
— оо —оо
/■ хг уг
где С = \ и п = \— пространственные частоты волнового поля, Хf и yf - координаты
/
точки в частотной плоскости
Чем больше пространственный угол в (рис. 4) индикатрисы рассеяния, тем большие пространственные частоты записываются на спекл-голограмму.
Рис. 4. Индикатриса рассеяния - 2 в точке А диффузно отражающего объекта - 3 при освещении по стрелке - 1
При оптическом дифференцировании, осуществляемом с помощью оптической фильтрации по пространственным частотам, имеется возможность регистрировать больший диапазон плоских компонент деформаций. Диффузная часть рассеяния излучения точками поверхности используется для регистрации плоских компонент пространственных деформаций, а зеркальная составляющая, имеющая большую интенсивность, используется для регистрации нормальной или изгибной составляющей деформаций.
В силу линейности преобразования Фурье уравнение (8) можно записать в одномерном виде. И в дальнейшем использовать для анализа деформаций объекта.
В плоскости голографирования осуществляется регистрация стационарного поля суперпозиции волн с распределением амплитуд (для простоты математических выражений по одной компоненте т(х) [12]), которое соответствует пространственной структуре диффузного рассеяния объекта. После фотохимической обработки амплитудное пропускание спекл-голограммы т(х) пропорционально зарегистрированной ею интенсивности. В плоскости голограммы объектный пучок интерферирует с каждой из пространственных составляющих волн. Поэтому на спекл-голограмме регистрируется соответствующая совокупность интерференционных решеток, периоды которых определяются углами интерференции.
Для простоты сумму пространственных составляющих частот поля спекл-структуры можно записать покомпонентно:
т(х) = То - к\г(х)\2 - к £ \г(х)\2 еЩ)(Й^) - к £ \г(х)\2 еЩ)(. (11)
\ Йпт ) ^ \ Йпт )
п,т 4 ' п,т 4 '
Результат фильтрации в фурье-плоскости может быть записан в следующем виде:
и (() = Ит (х)) = То 6(() - кИ \г(х)\2) - к^Т (() ® Т*(( + (пт) - к^Т (() ® Т* (( - (пт),
^ ' ^ ' п,т п,т
(12)
где (пт = - пространственная частота спекл-структуры п-го и т-го излучателей, используемая при определении цены полосы спекл-интерферограммы поля производной перемещений по полю спекл-голограммы, Т(() - фурье-образ функции t(x), К,к - постоянные процесса регистрации (проявления), Е, ® - соответственно операторы фурье-преобразования и автокорреляции, * - знак комплексно-сопряженной функции.
Распределение амплитуд для изображения спекл-поля после фильтрации и обратного фурье-преобразования в фотоизображении можно записать в виде [9, 11]:
U (x') = k £\t(x')\2 exJi^ )X sin 9nm\ + k £\t(x')\2 exJ-i(^ )X sin 0nm), (13)
n,m ^ ' n,m ^ '
где sin Onm = \£nm.
Формула (13) для пространственных частот при фильтрации в частотной плоскости используется при расшифровке спекл-голограмм для выделения плоских компонент деформации объекта. Анализ пропускания спекл-голограмм двойной экспозиции при наличии деформации и смещения как целого был опубликован в работах [7, 10]. Впервые этот подход изложен в работах [4, 6, 10], дальнейшее его развитие отражено в работах [7, 12].
4. Измерение плоских компонент тензора деформаций
Смещению в плоскости изображений соответствует изменение пространственной частоты в фурье-плоскости, а в случае неоднородного смещения формируется определенный спектр пространственных частот, соответствующий деформациям. Способ оптического дифференцирования спекл-голограмм заключается в их оптической фильтрации. Выделяются отдельные пространственные частоты и на них строится оптическое изображение, отрожающее производную по полю перемещений:
Ч = и (14)
оиф
где ф - азимут угла фильтрации поля пространственных перемещений U по направлению плоской компоненты иф. Таким образом, при оптической фильтрации выделяются квазирегулярные структуры, соответствующие несущей пространственной частоте, по которой определяется иф:
иф = —, (15)
ф sin а
где Л - длина волны излучения, а - пространственный угол фильтрации, равный для малых углов (до 10°).
а = (16)
где а - расстояние от фильтрующего отверстия до оптической оси линзы, f - фокусное расстояние линзы.
По интерферограмме вычисляется деформация объекта. Вторая производная оптического поля (перемещения) выделяется путем двойной фильтрации оптического поля спекл-голограммы (рис. 5).
Рис. 5. Оптическая схема фильтрации: 1 - спекл-голограмма с функцией пропускания т(х), 2 -частотная плоскость линзы Ы1 д(и) = Г[Ат(х)] с выделенной пространственной частотой щ, 3 -плоскость интерферограммы т'(х') = Г[д(и)] , 4 - частотная плоскость линзы Ы3 д'(и') = Г[т'(х')], где выделяется вторая производная т(х) с пространственной частотой и'0
Аналогично выделяются остальные производные и по ним определяются главные деформации.
5. Практическое использование спекл-голографии
Возможная схема испытания натурной конструкции самолета с помощью регистрирующей голографической установки приведена на рис. 6.
В качестве тестовой задачи был выбран плоский образец с отверстием. Конфигурации полос трех ортогональных компонент перемещений и, V, и> (рис. 7) испытанного на растяжение образца полностью совпадают с аналогичными картинами полос, полученными методами Муара (геометрической и дифракционной оптики) для плоских компонент. Картина полос нормальной компоненты и> (рис. 7а) согласуется с аналогичными результатами других авторов. Главная компонента деформации приведена на рис. 7б, по которой указанным ранее способом вычисляется деформация по всему полю интерферограммы и коэффициент концентрации напряжений в упругой зоне. Погрешность измерений составила 2-5%.
На рисунке 8 приведены интерферограммы исследования влияния коррозионного повреждения (1) в образце из алюминиевого сплава.
Рис. 6. Схема испытания: 1 - лазер; 2 - спекл-голограмма; стрелками показано направление нагру-жения
Рис. 7. Интерферограммы пространственных компонент деформаций образца с отверстием а) = 0, 3 мкм; б) п\ = 3, 0 мкм; в) = 5, 2 мкм
Рис. 8. Выявление зон трещинообразования в элементе лонжерона: 1 - зона коррозионного повреждения, 2 - подповерхностная трещина, 3 - болтовое соединение лонжерона, 4 - выход трещины на поверхность лонжерона
6. Заключение и выводы
Успешное сочетание методов спекл-интерферометрии и обычной классической голоин-терферометрии, подтверждающее общность физического механизма классической гологра-
фии и оптики спеклов, позволяет существенно повысить информативность голографиче-ской интерферометрии и более конкретно говорить о связи голографии и оптики спеклов. расширяющих практическое применение спекл-голографии в различных областях науки и техники.
Разработан метод оптической фильтрации спекл-голограмм, который позволяет выделять отдельные компоненты пространственных перемещений и определять компоненты тензора деформаций.
Литература
1. Журавлев О.А., Комаров С.Ю., Попов К.Н., Прокофьев А.Б. Разработка автоматизированного метода исследования вибрационных характеристик энергоустановок // Компьютерная оптика. 2001. Вып. 21. С. 143-149.
2. Джоунс Р., Уайкс К. Голографическая и спекл- интерферометрия. М.: Мир, 1986.
3. Клименко И.С. Голография сфокусированных изображений и спекл-интерферометрия. М.: Наука, 1985.
4. Волков И.В. Исследование механики деформаций натурных образцов с концентраторами напряжений методом speckle голографии // Материалы 6-ой Всесоюзной Школы по голографии АН СССР. Л.: ЛИЯФ, 1974. С. 369-382.
5. Волков И.В. Применение метода спекл-голографии двойной апертуры для выделения плоского компонента деформации вблизи концентраторов напряжений // Ученые записки ЦАГИ. 1976. Т. VII, N 4. С.168-173.
6. Волков И.В., Клименко И.С. О некоторых особенностях получения и интерпретации спекл-интерферограмм деформируемых объектов // Журнал Технической Физики. 1980. Т. 50, вып. 5. С.1038-1048.
7. Волков И.В. Спекл-голография в экспериментальной механике. Пенза: ПГТА, 2010.
8. Abramson N. Sandwich hologram interferometry: a new dimension in holographic comparison // Applied Optics. 1974. V. 13(9). P. 2019-2025.
9. Волков И.В. Приложения голографической и спекл-интерферометрии в экспериментальной механике. Части 1, 2 // Новые промышленные технологии, Минатом России. 2006. N 3. С. 44-52; N 4. С. 45-51.
10. Волков И.В., Клименко И.С., Кварацхелия Т.Г., Голикова Н.А. К интерпретации спекл-интерферограмм деформируемых объектов // Журнал Технической Физики. 1981. T. 51, вып. 10. C.2080-2085.
11. Волков И.В., Жебракова Г.В., Короленко, Т.И. Использование когерентно-оптических методов в анализе деформаций конструкций // М.: ЦАГИ, ОНТИ. Обзор N 707. 1990.
12. Волков И.В. Внестендовая спекл-голография. Использование голографической и спекл-интерферометрии при измерении деформаций натурных конструкций // Компьютерная оптика. 2010. Т. 34, N 1. C. 82-89.
References
1. Juravlev O.A., Komarov S.J., Popov K.N., Prokofiev A.B. Development of an automated method for studying the vibration characteristics of power plants. Computer optics. 2001. V. 21. P.143-149. (in Russian).
2. Jones R.,Wykes C. Holographic and speckle interferometry. M.:Mir, 1986. (in Russian).
3. Klimenko I.S. Holography of focused-images and speckle interferometry. M.: Science, 1985. (in Russian).
4. Volkov I.V. Investigation of the deformation mechanics of full-scale samples with stress concentrators by speckle holography. Materials of 6th All-union School on holography USSR Academy of Science. L.: LINP, 1974. P.369-382. (in Russian).
5. Volkov I.V. Application of the double-aperture speckle holography method for isolating the flat strain component near the stress concentrators // Uchenye zapiski TsAGI. 1976. V. VII. N 4. P. 168-173. —(in Russian).
6. Volkov I.V., Klimenko I.S. Some features of obtaining and interpreting speckle interferograms of deformable objects. Journal of Theoretical Physics. 1980. V. 50, N 5. P.1038-1048. (in Russian).
7. Volkov I.V. Speckle holography in experimental mechanics. Penza: PGTA. 2010. (in Russian).
8. Abramson N. Sandwich hologram interferometry: a new dimension in holographic comparison. Applied Optics. 1974. V. 13(9). P. 2019-2025.
9. Volkov I.V. Applications of holographic and speckle interferometry in experimental mechanics. Part 1, 2. New industrial technologies, Minatom of Russia. 2006. N 3. C. 44-52; N 4. P. 45-51. (in Russian).
10. Volkov I.V., Klimenko I.S., Kvarazhelia T.G., Golikova N.A. Interpretation of speckle interferograms of deformable objects. Journal of Technical Physics. 1981. V. 51, N 10. C. 2080-2085. (in Russian).
11. Volkov I.V., Jebrakova G.V., Korolenko T.I. Application of coherent-optical methods in structural deformations analysis. Moscow: TsAGI, ONTI. Review N 707. 1990. (in Russian).
12. Volkov I.V. Off-stage speckle-holography. Holographic and speckle interferometry in the full-scale constructions measurments. Computer optics. 2010. V. 34, N 1. P. 82-89. (in Russian).
Поступила в редакцию 23.05.2017
