CC BY
55
УДК 535.31
АВТОМАТИЧЕСКАЯ СИСТЕМА КОНТРОЛЯ СОСТОЯНИЯ ЭЛЕМЕНТОВ ОСТЕКЛЕНИЯ ВОЗДУШНЫХ СУДОВ
© 2012 И.Э. Вольф1, А.Н. Малов2
1 ВУНЦ ВВС "Военно-воздушная академия имени проф. Н.Е. Жуковского и Ю.А. Гагарина", г. Воронеж 2 Иркутский государственный медицинский университет
Поступила в редакцию 10.10.2012
Автоматическая система контроля состояния элементов остекления воздушных судов разработана на основе использования методов пошагового компьютерного восстановления записанного волнового фронта изменения интенсивности лазерного излучения при прохождении его через объект диагностики и метода анализа спекл-структур регистрируемых на ПЗС-матрицу при отражении лазерного излучения от классифицируемого дефекта.
Ключевые слова: голографическое изображение, спекл-структура, дефектоскопия.
Для остекления кабин современных воздушных судов используют материалы органического и неорганического происхождения. В процессе хранения, монтажа и эксплуатации, остекление подвергается различного рода воздействиям: аэродинамическим, температурным и механическим. Определяющее влияние на разрушение его структуры оказывают большие напряжения от аэродинамических усилий и температурных перепадов во время выполнения полёта на практический потолок при максимальном ускорении и последующем торможении. В соответствии с регламентом обслуживания инженерно-техническим составом производится только субъективный оптико-визуальный контроль. Используемые методы и приборы не позволяют классифицировать дефекты остекления на соответствие техническим условиям, снижающие степень пригодности элементов остекления к дальнейшей эксплуатации, а также осуществлять мониторинг его состояния [1].
Анализ характеристик современных образцов остекления, а также технических условий применения их на воздушных судах показал, что при оценке внутренних дефектов аппаратура контроля параметров должна обеспечивать измерение геометрических размеров малоконтрастных объектов диаметром от 0,1 до 5 мм с погрешностью не более 0,01 мм и вероятностью обнаружения дефектов диаметром 0,1 мм не менее 0,95. Достижение заданных требований и автоматизацию процесса измерений, возможно за счет применения принципиально новых алгоритмов с использованием активных методов Вольф Игорь Эдуардович, старший преподаватель. E-mail: wolfigor71 @mail.ru
Малов Александр Николаевич, доктор физико-математических наук, профессор. E-mail: [email protected]
неразрушающего контроля на основе полупроводниковых лазерных модулей, средств регистрации и вычисления.
Появление быстродействующих цифровых матричных фотокамер с достаточно большим числом элементов, позволяют регистрировать массив фотоотсчетов, в реальном масштабе времени и в цифровом формате, согласованном с программным обеспечением компьютера [2]. Предлагается регистрировать голограмму участка элемента остекления на фотоматрицу с последующим компьютерным восстановлением изображений двумерных сечений голографируемого предметного поля.
Исследования были выполнены на экспериментальном стенде с использованием осевой схемы записи голограмм. Объект освещался когерентным параллельным пучком, прошедший через него свет состоит из двух компонент: 1) интенсивной однородной плоской волны, и 2) слабой рассеянной волны, состоящей из совокупности сферических волн, обусловленных дифракцией на элементах предмета. Таким образом, объект в некотором роде сам является источником опорной волны. В результате интерференции света, прямо прошедшего через объект и рассеянного им, получается распределение интенсивности, которое зависит как от амплитуды, так и от фазы рассеянной волны [3].
Распределение интенсивности I (х, у) в плоскости записи голограммы (х, у) определяется квадратом модуля суммы комплексных амплитуд объектной 0(х, у) и опорной Я(х, у) волн, а именно:
I (х, у) = 0 (х, у) + Я (х, у )2| = = Я (х, у )|2 + 0 (х, у )|2 + О (х, у) Я •( х, у) + .
+ Я(х, у)0'(х, у)
где значок * обозначает комплексно-сопряженную величину, O(х, у) = \й(x, y)\exp(/^0(x, y)) - комплексная амплитуда объектной волны с вещественной амплитудой O(x,у) и фазой <Ро(х,у), R(x у) = |R(xy)exp(/^R(x, у)) обозначает аналогичным образом комплексную амплитуду опорной волны.
Исходная объектная волна восстанавливается путем освещения голограммы опорной волной, как и для оптической схемы из рис. 1.
Мнимое изображение Д восстанавливает объект в трех измерениях. Четвертое слагаемое соответствует восстановленному действительному изображению Д объекта. Некоторые свойства действительного изображения делают его, однако, менее интересным, чем мнимое изображение.
Для пояснения математического принципа, лежащего в основе голографической реконструкции изображения, умножим записанную голограмму, т. е. интенсивность интерференционной картины I (x, у) на значения поля опорной волны R(x, у):
R(x, у)^, у)=R(x, у)| R(x, у)|2 + R(x, у) |O(x, у)|2 +
+ |R(x,у)|2 O(x,у) + R2(x,у)0*(x,у).
Первое слагаемое в правой части уравнения пропорционально значению R (x, у) опорной волны, второе слагаемое описывает пространственно изменяющееся "облако" вокруг первого порядка дифракции. Вместе они составляют нулевой порядок дифракции, или фоновую составляющую. Третье слагаемое с точностью до известного множителя |R( x, у)| представляет собой точную копию исходного поля O (x, у) = |O( x, у) exp(i^0( x, у)) и поэтому называется мнимым изображением (или просто изображением) объекта. Четвертое слагаемое относится к еще одному изображению объекта, "изображению-двойнику", которое называется действительным изображением.
В процессе записи цифровой голограммы,
исследуемый трехмерный объект с рассеивающей свет поверхностью располагается на расстоянии $ от ПЗС-матрицы. Записанная голограмма может рассматриваться как пропускающий транспарант, на котором дифрагирует опорная волна. Для получения исследуемой объектной волны в цифровой голографии применяется численный расчет распространения оптического поля в форме произведения Я(х, у)I(х, у) от плоскости голограммы до плоскости объекта (%, п).
Численное восстановление записанной цифровым образом голограммы производится в соответствии со скалярной теорией дифракции в приближении Френеля для интеграла дифракции Рэлея-Зоммерфельда. Восстановленное дифрагированное поле Q(%,ц) в плоскости изображения (%, п) на расстоянии $ от плоскости голограммы может быть представлено в параксиальном приближении следующим образом:
1 ~ ад ад
Q(£n) = ехР('' 2Лй)\I Я( х, у) I (х, У)х
X
х ехр V ~м х )2+(п-у )2 ]} ахау'
Данное уравнение служит отправной точкой для численного восстановления изображений в цифровой голографии в параксиальном приближении, при котором значения координат х и у, как и соответствующих им ^ и П, малы по сравнению с расстоянием й.
Из уравнения видно, что результирующее поле определяется двумерным Фурье-преобразованием произведения распределения интенсивности в голограмме I (х, у) комплексной амплитуды опорной волны Я (х, у) и квадратичной фазовой функции пространственного распространения волны:
w( x, у) = exp
п
i—(x2 Xd
у2)
Отметим, что интеграл является интегралом свертки первых двух сомножителей и функции.
Интенсивность 1(х, у, и фаза
р(х, у, а) восстановленного изображения могут быть получены из рассчитанного на расстоянии $ комплексного поля Q(£, п) с помощью следующих соотношений:
1 (x, у,d) = Q( x, у )|2
Рис. 1. Восстановление изображения объекта, диффузно рассеивающего свет
р( х, у, 0) = са-Щ х у)))) = arg(Q (х, у)) .
Значения фазы, р(х, у, а) полученные по этой формуле, соответствуют разрывной функции с областью изменения в пределах интервала [—п, п].
ад —ад
Для восстановления непрерывной развернутой картины значений фазы может быть применен один из широко известных алгоритмов развертывания фазы.
Распространенным методом численного восстановления записанного голографического изображения является вычисление распространения углового спектра, или так называемый "метод свертки". В этом случае выражение для результирующего поля записывается с учетом параксиального приближения в следующем виде:
I 2п Г , 2 2
Q(Уг;,vv) = ехРу\Л|_ 3 ехр(/'п\V + V))х
х З+. х, у)й( х, у)(у,,уп) ]
Здесь в явном виде используется результат применения Фурье-преобразования к квадратичному фазовому множителю w(х, у) в форме:
3+1 х, у)(у,, V ) = Ш ехр(-/п\(V + V )),
где З11 прямой (+1) и обратный (-1) непрерывные двумерные интегралы Фурье:
3й
F-1 {F {h • г} • F {g}}
Самый простой, но достаточно трудоемкий способ определения фазы восстановленного волнового поля - это просто вычислить фазу как отношение мнимой и действительной части комплексной амплитуды. Однако сделать это возможно только после второго Фурье-преобразования, т.е. после того, как мы избавимся от несущей и суперпозиции опорного и объектного пучков. Отметим, что хотя первичный массив является реальной функцией, его Фурье-образ уже функция комплексная. После второго Фурье-преобразования полученный массив и0( х, у, ) также будет комплексным. Тогда фаза в каждой точке определяется как:
ф(x, y, zt) = arctg
Im[uo( x y, z)]
f(x,у)(Уц,У„) = J J f(x,y)exp[+2п(^х+Vvy)dxdy.
Первый шаг - восстановление голограммы. Происходит моделирование опорной волны, т.е. используя математический аппарат, соответствующий освещению голограммы опорным пучком. Математический аппарат - есть воздействие на голограмму "Н" прямым и обратным преобразованиями Фурье:
Мио(х у, )] •
Здесь 2 1 обозначает, что фаза определяется на определенном расстоянии от голограммы (ПЗС-матрицы).
Другой метод расчета фазы позволяет обойтись без цифровой реконструкции (двойного Фурье-преобразования), однако в этом случае полосы пространственной частоты должны иметь расстояние d между собой больше, чем 3 пикселя, т.е. Л > 3Ах. В этом случае интенсивности на 3-х соседних пикселях можно записать в следующем виде:
11 = 210 [1 + соэф -Фя -Аф)]
12 = 210 [1 + СОБ(ф0 -фя )]
13 = 210 [1 + соб(Фо -Фя + Аф)]
где Аф =
2пАх d
Тогда, исключая из уравнений 10 и pR, получим для <р0
где r = reference, в нашем случае он равен единице (плоская волна).
Фо( x, У) = arctg
Ii -13
Л +13 - 212
Аф
tg
ехРШЛ'2+(*-+(1 -1}2АП |
^ = \ \ Л, Л, ]
Второй шаг - дифракционный интеграл от каждой точки голограммы вычисляют с заданным шагом, после чего получают ряд плоскостей - слайсы. Рассматривая каждую из полученных плоскостей, можно продвигаться вглубь голограммы, фокусируясь на необходимых расстояниях. Из каждой голограммы восстанавливают численно структуру записанного волнового фронта и находят 3D координаты частиц.
Применяя данную операцию последовательно ко всем тройкам соседних пикселей, т.е. сдвигая каждый раз зарегистрированный массив на ±Ах , получим распределение фаз для предметного пучка.
Дальнейшая обработка голограммы и извлечение информации о координатах и размерах дефектов полностью лежит на плечах компьютера или вычислительной системы. Восстановление изображения происходит по принципам компьютерной оптики. Голограмма облучается "виртуальным" опорным волновым фронтом (с параметрами реального опорного фронта). При помощи интеграла дифракции Рэлея-Зоммер-
фельда, рассчитан результат дифракции "виртуального" опорного фронта на цифровой голограмме на некотором расстоянии Z. Этим изображением является проекция в плоскости XY трехмерного распределения координат дефекта. Критерием фокусного слоя, является минимум функции зависимости усредненной интенсивности от координаты Z.
Исследуя участки остекления наиболее подверженные воздействию аэродинамических нагрузок методом пошагового восстановления компьютерного голографического изображения, становится возможным выявление внутренних дефектов на начальной стадии, при их образовании. При детальном рассмотрении полученного восстановленного изображения выявлены не только видимые дефекты поверхности, но и невидимые при визуальном осмотре, которые являются основными при дефектоскопии авиационных остеклений (фонари кабин, фотолюки, иллюминаторы).
Для оперативного определения параметров дефекта также используется метод анализа профиля рельефа и структуры дефекта в полупрозрачном объекте по рассеянному излучению. Получаемые оптические поля имеют сложную амплитудно-фазовую структуру случайного характера, образующиеся при отражении или прохождении когерентной волны через случайно неоднородные объекты, например, при отражении лазерного пучка от шероховатой поверхности. Пространственное распределение интенсивности таких полей носит случайный характер с теми или иными статистическими свойствами.
Интерференция частично когерентного света связана с проявлением корреляционных свойств временных флуктуаций поля, обусловленных процессами, развивающимися в источниках света и оптических элементах схемы.
Учитывая, что при обработке записанных спекл-изображений отраженного от объекта излучения, информативные признаки последних
не изменяются и несут туже информацию об объекте, для создания прибора целесообразно рассмотреть оптическую схему, представленную на рис. 2.
Когерентное излучение от двух полупроводниковых лазерных модулей с различной рабочей частотой, после объединения и выделения самой яркой моды, освещало исследуемый объект и отражалось . При этом часть отраженного излучения фиксировалась ПЗС-матрицей цифрового фотоаппарата.
Применительно к поставленной задаче, схема портативного устройства имеет вид [4] (рис. 2).
Когерентное излучение попадает на дефект 5 контролируемого объекта 1 и отражаясь от микроскопичных образований на всю глубину объекта, в случае если он является полупрозрачным или прозрачным, попадает на матричный фотоприемник, исключая оптические элементы.
С полученными изображениями проводим корреляционную обработку, по стандартным и оригинальным программам, разработанным в среде MathCAD.
При зондировании малоконтрастных дефектов целесообразно использовать для получения спекл-изображений две рабочих частоты лазерного излучения. Данный вариант зондирования позволяет повысить точность измерений за счет равномерного распределения интервала корреляции Дx=f(SN). Не исключая возможности получения трех спекл-изображений при различных вариантах освещения получаем большую степень свободы для принятия решения и исключения ошибки.
Способ контроля элементов остекления воздушных судов на основе выше рассмотренных методов является неразрушающим, без демонтажа исследуемого элемента и позволяет объективно произвести их оперативный контроль, оценить параметры дефекта, определить размеры площади SN (мм2) дефекта, протяженности Дh (мм) с высокой точностью и проводить мо-
Рис. 2. Оптическая схема-макет устройства для получения спекл-картин: 1 - исследуемый объект, 2, 3 - полупроводниковые лазеры, 4 - матрица цифрового фотоаппарата (камеры), 5 - дефект
ниторинг объекта. Получение результатов контроля в меньшей мере зависит от классификации специалиста, проводящего осмотр и диагностику данных элементов, и позволяют технику с низкой квалификацией выполнять данную операцию и исключить влияние субъективной составляющей на результат контроля.
С учетом полученных экспериментальных результатов, составлен алгоритм выявления, оценки, классификации и мониторинга дефектов остекления воздушных судов в процессе эксплуатации.
Полученный алгоритм позволяет максимально исключить неточности проведения контроля элементов остекления, вызванные некомпетентностью инженерно-технического состава, так как основные операции по обработке производится программными методами, с использованием портативных приборов и средств вычислительной техники. Использование данного алгоритма позволяет достаточно достоверно определить наличие малоконтрастного дефекта, произвести его оценку и классификацию, а при необходимости произвести дополнительное исследование этого дефекта с целью принятия правильного решения
по продолжению эксплуатации данного объекта, либо его ремонта или замены. Использование автоматической системы оценки, позволит производить данные операции персоналу, не имеющему высокой квалификации.
СПИСОК ЛИТЕРАТУРЫ
1. Вольф И.Э., Малов A.H., Сычевский A.B. Диагностика подповехностных дефектов прозрачных и полупрозрачных изделий // Книга "НАУЧНАЯ СЕССИЯ МИФИ-2008. Том 2. Физика плазмы и управляемый термоядерный синтез. Лазерная физика. Фотоника и информационная оптика. Физика конденсированного состояния вещества" М.: МИФИ, 2008. С.110-111.
2. Малов A.H. Лазерные методы бесконтактной диагностики машиностроительных изделий // Физика наукоемких технологий: сборник научных статей. Иркутск: ИВВАИУ. 2006. С. 41-68.
3. Burley C. The HART-11 Test: Rotor Wakes and Aeroacoustics with Higher-Harmonic Pitch Control (HHC) Inputs / Casey Burley, Thomas Brook. -Hampton: NASA Langley Research Center, 2000. 120 с.
4. Вольф И.Э., Малов A.H., Бородин A.H. Мультиспект-ральный спекл-метод исследования полупрозрачных объектов // Известия Самарского научного центра РАН. 2009. Т. 11(31). № 5. С.42-45.
THE AUTOMATIC MONITORING SYSTEM OF A CONDITION OF ELEMENTS OF A GLASS COVER OF AIRCRAFTS
© 2012 I.E. Wolf1, A.N. Malov2
1 Military Scientific-Educational Centre of Air Forces "Military-Air Academy named after prof. N.E. Zhukovsky and Jr.A. Gagarin", Voronezh 2 Irkutsk State Medical University
The automatic monitoring system of aircraft glass elements condition is developed on the basis of step-by-step computer restoration methods of the written down wavefront of laser radiation intensity change at its passage through testing object and a method of the analysis of speckle-structures registered on the CCDmatrix at reflection of laser radiation from classified defect. Keywords: the holographic image, spekl-structure, defect definition.
Igor Wolf, Senior Lecturer. E-mail: [email protected] Alexander Malov, Doctor of Physics and Mathematics, Professor. E-mail: [email protected]
