CC BY
51
АНАЛИЗ ПОГРЕШНОСТЕЙ ПРИ ЦИФРОВОЙ ОБРАБОТКЕ РЕЗУЛЬТАТОВ
ИНТЕРФЕРОМЕТРИЧЕСКОГО КОНТРОЛЯ ЛОКАЛЬНЫХ ОТКЛОНЕНИЙ
НАНОМЕТРОВОГО УРОВНЯ ПОВЕРХНОСТЕЙ ОПТИЧЕСКИХ ДЕТАЛЕЙ
Д.Г. Денисов1
'Московский государственный технический университет имени Н.Э. Баумана, Москва, Россия
Аннотация
Разработан, научно обоснован и экспериментально подтверждён метод динамической интерферометрии контроля локальных отклонений нанометрового уровня поверхностей оптических деталей от заданного профиля на основе алгоритма расчёта целевой функции -спектральной плотности одномерной корреляционной функции. Представлены теоретические и экспериментальные исследования, посвящённые определению среднего квадратиче-ского отклонения локальных отклонений поверхностей оптических деталей диаметром до 100 мм и до 1000 мм, с учётом неисключённой систематической и случайной составляющих погрешностей определения целевой функции.
Ключевые слова: оптический контроль, интерферометрия, измерения поверхности, спектральная плотность мощности, методические погрешности, краевой эффект, эффект «утечки» частоты.
Цитирование: Денисов, Д.Г. Анализ погрешностей при цифровой обработке результатов интерферометрического контроля локальных отклонений нанометрового уровня поверхностей оптических деталей / Д.Г. Денисов // Компьютерная оптика. - 2017. - Т.41, № 6. -С. 820-830. - Б01: 10.18287/2412-6179-2017-41-6-820-830.
Введение
Современные оптические технологии, ориентированные на создание перспективных высокотехнологичных изделий, обладающих высокими показателями качества [1-3], на сегодняшний день невозможны без применения высокоточных методов и оптико-электронной аппаратуры аттестационного контроля.
Одним из перспективных научных направлений в области современных оптических технологий является задача создания высококачественной крупногабаритной лазерной оптики для реализации управляемого термоядерного синтеза [4]. Научно-технический задел в данном направлении существует как у отечественных («Дельфин-1» (ФИАН, Москва, СССР), «Искра-5» (РФЯЦ-ВНИИЭФ, г. Саров, РФ)), так и у зарубежных («NIF»-National Ignition Facility (США), «Мегаджоуль» (Франция), «ITER»-International Thermonuclear Experimental Reactor (США, Европа, Япония, Россия), «Nova»-Lawrence Livermore National Laboratory (США)) разработчиков, что позволяет в перспективе перейти на принципиально новый уровень получения выходной энергии - мегаджоульный. Получение такого уровня энергии требует высокого качества изготовления лазерных активных сред и разработки методов аттестационного контроля показателей их профилей [1- 6].
Для создания и успешного развития в России данного технологического направления необходимо иметь не только современную отечественную технологию производства крупногабаритных оптических элементов нового поколения, но и комплекс высокоточной оптико-электронной аппаратуры диагностики их статистических показателей качества [4-11].
1. Постановка научно-технической задачи
С целью отработки схемотехнических и конструкторских решений стендовой аппаратуры полноапер-турного (до 1000 мм) контроля локальных отклоне-
ний поверхностей крупногабаритных оптических изделий, а также оценки точностных возможностей разрабатываемого алгоритма расчёта величины среднего квадратического отклонения (СКО) в настоящей работе рассмотрен метод и его реализация на базе динамического интерферометра с диаметром выходного зрачка 100 мм и зеркально-линзового коллиматора с диаметров выходного зрачка 1000 мм. Предложенный метод основан на частотном анализе [5, 12-14] профиля контролируемой детали в заданном диапазоне масштабов локальных неоднородностей, совокупность которых описывает структуру её исследуемой поверхности.
В общем случае функцию, описывающую структуру локальных отклонений профиля оптической детали, можно представить в виде набора гармоник соответствующих пространственных частот, определяющих такие его технологические параметры, как непрямолинейность и волнистость [1- 8]. Эти локальные отклонения обусловлены свойствами обрабатывающих инструментов с погрешностью метода обработки оптических изделий [1-3]. Помимо указанных параметров формы профиля, высокочастотные гармоники характеризуют степень шероховатости, которая образуется при съёме материала в процессе формообразования детали полирующими суспензиями [1-3]. Для описания неоднородностей искажений волновых фронтов при контроле оптических элементов можно использовать международный классификатор качества изготовления активных лазерных сред технологического комплекса NIF (National Ignition Facility, США) [4, 5, 15], представленный в виде набора характеристик для четырёх пространственно-частотных диапазонов (табл. 1).
Диапазоны I-III характеризуют искажения параметров волнового фронта [7- 9] лазерного пучка, обусловленные ошибками формы исследуемых оптиче-
ских деталей, а четвёртый диапазон характеризует шероховатость профилей их поверхностей. Указанные в таблице величины искажений волнового фронта оцениваются в диапазоне I статистическими показателями PV (peak to valley - максимальное расстояние между высотой наибольшего выступа и глубиной наибольшей впадины искажений) [1- 8], RMSgrad (root mean square gradient - СКО градиента искажений) [4, 5], а в диапазонах II и III - параметром RMS [4 - 11], ха-
рактеризующим отклонения формы волнового фронта в высокочастотной области. Для оценки шероховатости в IV диапазоне используют параметр Rq [3- 6], который определяет СКО точек профиля поверхности оптической детали в пределах заданной области исследования (базовой площади). Аналитическое описание параметра Rq имеет тот же вид, что и величина RMS, но определяется на малой базовой площади [3- 6].
Табл. 1. Технологические требования к оптическим деталям при пространственно-частотном анализе неоднородностей поверхности
Характеристика Пространственно-частотный диапазон
I II III IV
Пространственная частота неоднородностей искажений волнового фронта, мм-1 0 - 0,03 0,03-0,10 0,4-1,0 10-100
Пространственный масштаб неоднородностей искажений волнового фронта, мм до 33 33-2,5 2,5-0,1 0,1-0,01
Статистические показатели искажений волнового фронта:
PV, нм 208 - - -
RMSgrad, нм/см 7 - - -
RMS, нм - 1,8 1,1 -
Rq, нм - - - 0,4
Стоит отметить, что при переходе от указанных статистических показателей искажений волнового фронта к параметрам профиля [7-11] исследуемой детали вводится нормирующий множитель 0,5, впоследствии учитываемый в алгоритме расчёта.
Важно отметить, что приведённые в табл. 1 показатели качества являются оптимизационными параметрами целевой функции аттестационного контроля профиля оптической детали, в качестве которой предлагается использовать спектральную плотность мощности сигнала (СПМ от англ. PSD (Power Spectral Density). Дословный перевод используемой целевой функции в большинстве случаев применяется при преобразовании «электрических» сигналов [12, 13]. В данной работе при обработке оптического сигнала целесообразнее говорить о спектральной плотности одномерной корреляционной функции профиля поверхности (СПКФ1) или S1D(v), нм2-мм (где v - пространственная частота, 1/мм). В качестве критериев оптимизации данной функции используются допустимые (максимальные) значения погрешностей определения приведённых в табл. 1 показателей качества.
Традиционный метод аттестационного контроля [4 - 5, 15] основан на сравнении построенной функции S1D(v) с функцией S1Diso(v), нм2-мм, определяемой согласно стандарту ISO 10110-6:2015 [15].
s1d;so (v) = 4
1
1
<v<-, (1)
v" 1000 x D 1000 x C
где S1DiSO(v) - СПКФ1, построенная в соответствии со
стандартом ISO 10110-6:2015; А - константа (1 для
стекла [4 - 6, 15], 15 для кристаллических сред [4 - 6,
15]); 1 < B< 3 - константа [15]; С - предел разрешения
прибора (мкм) [15]; D - длина, на которой производится
измерение (мкм) [15]; v - пространственная частота,
мкм-1 [15].
При превышении функции S1D(v) в указанном пространственно-частотном диапазоне функции
S1DIso(v) можно судить о наличии недопустимых локальных отклонений поверхности от заданного профиля и говорить о непригодности оптической детали для её последующего функционирования (рис. 1а, б). Детально данная методика приведена в [45, 15].
Следует отметить, что возможные пространственные неоднородности в промежуточных сечениях исходной функции S2DISO(vx, уу), нм2-мм2, обусловленные наличием локальных отклонений поверхности (рис. 1 а), могут внести существенный вклад в структуру погрешности определения результирующего значения какого-либо статистического показателя качества исследуемой формы волнового фронта. Поэтому данный метод аттестационного контроля не позволяет судить об истинном качестве оптической детали в широком диапазоне пространственных частот локальных отклонений.
В данной работе предлагается альтернативный метод, основанный на использовании алгоритма расчёта S1D(v) в результате статистического усреднения отдельно взятых направлений вдоль радиуса-вектора в информативной матрице исходной функции S2D(vx, Vy).
Для получения функции S2D(vx, Vy) можно использовать два подхода. Первый базируется на теореме Винера-Хинчина [4, 5, 14] и теории статистических измерений и позволяет определить искомую функцию как преобразование Фурье [4, 5, 12 - 14] от автокорреляционной функции распределения высот исследуемой карты профиля в зависимости от координат Н(х,у) [3 - 5, 8]:
. 4/2 а 2
S2D(vх,Пх) = — • I } АКг [н (х,у)] х
-А/2 -a/2
(2)
x exp{-2pi (vxx + v yy)} dx dy,
где а, Ь - размеры световой апертуры детали, характеризующие область исследования вдоль осей х, у соответственно; АКг[Н(х, у)] - автокорреляционная
функция или математическое ожидание произведения разностей высот профиля ^(х,у)]; г,у - разность координат высот профиля.
Локальный
поверхностный дефект
ySlD(v),
нм2-мм SlDisojv)
ч
гГТН/л Локальный . поверхностный
дефект
0,01 0,18 2,4 6 Шкала частот для максимальной области исследования 1000 мм
J
Vjb мм~'
L
Шкала частот для максимальной области исследования 1000 мм
J
б)
Рис. 1. Иллюстрация методики аттестационного контроля локальных отклонений поверхности оптического элемента технологического комплекса NIF: вид исходной функции S2Diso(vx, vy) и её сечение а) под углом 60°; б) соответствующая сечению усреднённая функция S1D(v ) и построенная согласно стандарту
ISO 10110-6:2015 функция S1Diso(v) Реализация второго метода расчёта функции S2D(vx, vy), как частный случай первого (2), предполагает нахождение искомой двумерной функции стационарного в широком смысле случайного процесса как квадрат модуля преобразования Фурье высот карты профиля оптической поверхности в зависимости от координат h(x, y), приведённой к единичной площади a*b (области исследований)
|2
где
S 2D (V x, V y ) = |3{h( x, y)}|2/(ab), (3)
a/2 4/2
3{h(x,y)} = J J h(x,y)exp(-2ni(VxX+Vyy)) dxdy -
преобразование Фурье от карты профиля оптической поверхности в зависимости от координат ^х, у).
Для дальнейших исследований будем использовать второй метод, поскольку его применение обеспечивает более высокую точность вычислений. Приведённый ал-
горитм позволяет проанализировать пространственные частоты и количественно рассчитать СКО соответствующих им неоднородностей поверхности (табл. 1) на основе теоремы Парсеваля [4, 5, 14]:
RMSS1D (V) =
2я J S1D(v)dv
(4)
где vmax, vmm - значения максимальной и минимальнои пространственных частот исследуемой области (частотный диапазон исследования) на поверхности детали (табл. 1).
Стоит отметить, что значения границ диапазонов (Vmax— Vmm) указанной частотной шкалы (4) определяются, с одной стороны, линейной величиной элемента пикселя матрицы приёмника излучения, а с другой стороны, размером световой апертурой контролируемой оптической детали. Указанные в табл. 1 требования к искажению волнового фронта лазерного пучка приводятся для световой апертуры 603 мм (диаметральный размер, соответствующий линейному размеру изделия 1000 мм). Для исследования качества поверхностей оптических деталей других диаметров необходимо ввести коэффициент масштабирования шкалы частот K = 603/D. Так, для плоской оптической детали D=100 мм пространственно-частотные диапазоны I, II и III будут составлять соответственно Avi = 0,01...0,18 мм —1; Avu = 0,18.2,4 мм —1; AvIII = 2,4.6 мм —1.Частотные шкалы диагностики целевой функции SlD(v) при расчёте СКО точек профиля поверхности оптической детали в пределах заданной области исследования диаметром 100 мм и 603 мм представлены на рис. 16. В первом случае для проведения экспериментальных исследований качества оптической поверхности используется лабораторный измерительный стенд на базе динамического интерферометра c выходной апертурой 100 мм [5], а во втором случае для расширения световой апертуры совместно с последним применяется зеркально-линзовый коллиматор (рис. 2), имеющий диаметр выходного зрачка порядка 1000 мм.
Анализируя представленную выше связь пространственных частот, в работе предлагается на первом этапе провести исследование погрешности СКО профиля поверхности оптической детали в пределах апертуры диаметром 100 мм, а на втором этапе в пределах апертуры диаметром до 1000 мм. Как можно видеть из рис. 16, первый частотный диапазон во втором случае за счёт применения зеркально-линзового коллиматора уширяется в низкочастотной области примерно в 6 раз, что даёт возможность исследовать СКО профиля крупногабаритной оптической детали на всей световой апертуре.
2. Принцип работы лабораторного измерительного стенда
C целью отработки алгоритма полноапертурного (до 1000 мм) аттестационного контроля неоднород-
V
ностей искажений волновых фронтов при контроле крупногабаритных оптических деталей в I и II диапазонах соответствующих пространственных частот (табл. 1) в результате расчёта одномерной целевой функции - S1D(v) была разработана методика и измерительный стенд на базе динамического интерферометра «IntelliumH2000» (ESDI, США). Данный интерферометр построен по оптической схеме Физо и имеет выходную апертуру диаметром 100 мм. Подробное описание функциональной схемы и принципа работы прибора представлены в [5]. Измерительный стенд состоит из динамического интерферометра и оптико-механического узла крепления плоского эталонного (контролируемого) зеркала в оправе. Отличительной особенностью используемого в измерительной схеме динамического интерферометра является то, что сдвиг фаз регистрируемых интерферо-грамм реализуется с помощью трёх поляризаторов плоскости, поляризации которых ориентированы друг относительно друга под углом 120°. Эта особенность позволяет за один кадр получать три реализации при минимальном времени экспозиции 100 мкс при помощи широкоформатных матричных приёмников излучения, что обеспечивает нечувствительность метода контроля к турбулентным потокам воздушных масс и влиянию низкочастотных вибраций, существующих в производственных условия. Стоит отметить, что в используемом приборе реализуется поляризационная калибровка [5], позволяющая учитывать инструментальную погрешность измерительного стенда. Данная составляющая общей структуры системной инструментальной погрешности влияет на результат измерения СКО точек профиля крупногабаритной контролируемой поверхности оптической детали в пределах области исследования до 1000 мм. В этой связи особое внимание в экспериментальной части работы было уделено разработке функциональной схемы стенда полноапертурного контроля совместно с сотрудниками АО «Лыткаринский завод оптического стекла».
3. Принцип работы и состав аппаратуры
полноапертурного измерительного стенда
Для реализации полноапертурного аттестационного контроля СКО искажений волнового фронта лазерного пучка, отражённого от поверхности крупногабаритной оптической детали, в работе предложена функциональная схема измерительного стенда, представленного на рис. 2.
В представленной схеме полноапертурного аттестационного контроля, аналогично функциональной схеме лабораторного стенда [5] излучение на выходе из прибора, после прохождения фокусирующего объектива 1, попадет в измерительный канал, представляющий собой зеркально-линзовый коллиматор.
В состав последнего входят плоские поворотные зеркала 2, 3, предназначенные для уменьшения габаритных характеристик стенда, коллимирующий объектив 4, фокус которого сопряжён с фокусом объектива 1, после
которого плоский волновой фронт лазерного пучка проходит оптический клин с эталонной поверхностью сравнения 5, падая на исследуемую поверхность крупногабаритной плоской контролируемой детали 6. Реализуя в схеме поляризационную калибровку [5], можно последовательно регистрировать результат суперпозиции составляющих двух пар электромагнитных волн: (5э; pд) или (pэ; 5д) соответственно. Вычитая полученные результаты в виде интерферограмм, в схеме измерительного стенда минимизируется инструментальная погрешность, обусловленная погрешностями изготовления элементов зеркально-линзового коллиматора.
4 5 6
/1 \.....
\ ______
Динамический интерферо
метр, -------
енныи
pqjepu- л \_____—с
\nocmpo- Ц. _________\>
и по . " Tg ---А
ческой W ------//
г Физо L—V
оптическои схеме Физо
-м и
Рис. 2. Функциональная схема измерительного стенда на базе динамического интерферометра «1п1еШитН2000»
и зеркально-линзового коллиматора: 1 - фокусирующий объектив; 2, 3 - плоские поворотные зеркала коллиматора, 4 - коллимирующий объектив, 5 - оптический клин с эталонной поверхностью сравнения, 6 - крупногабаритная плоская контролируемая деталь
Для отработки экспериментальных данных, полученных на основе двух модификаций измерительных стендов в лабораторных и цеховых условиях, в работе рассмотрен метод расчёта статистических показателей, характеризующих искажения волнового фронта лазерного пучка (табл. 1), возникающие при контроле плоских оптических элементов.
4. Метод определения СКО локальных отклонений в пространственно-частотных диапазонах I, II
Разработанный метод определения локальных отклонений исследуемых профилей оптических деталей основан на алгоритме расчёта целевой функции -СПКФ1 (£Ш(у)), с учётом следующих неисключён-ных систематических погрешностей: дискретизации, математического аппарата быстрого преобразования Фурье (БПФ), преобразования двумерного стационарного сигнала в одномерный, краевых эффектов БПФ, эффекта спектральной «утечки» частоты [4 - 6, 12 - 14]. Анализируемый с точки зрения погрешностей метод основан на реализации следующих этапов:
1) регистрация исходной топографической карты профиля контролируемой оптической поверхности в зависимости от координат Н(х, у);
2) расчёт информативной матрицы исходной функции 520(Ух, Уу) (рис. 1а), приведённой к еди-
ничной области (площади) исследования в соответствии с (4);
3) расчёт матрицы суммарных пространственных частот;
4) расчёт целевой функции S1D(v);
5) расчёт оптимизационных параметров целевой функции в заданных пространственно-частотных диапазонах I и II в соответствии с табл. 1:
5.1) анализ влияния краевых эффектов дискретного БПФ на определение СКО исследуемого профиля;
5.2) анализ влияния эффекта спектральной «утечки» пространственных частот на определение СКО исследуемого профиля;
6) сравнение рассчитанной функции SlD(v) в заданных пространственно-частотных диапазонах с функцией S1D/so(v), определяемой согласно стандарту ISO 10110-6:2015 [15].
На первом этапе карта профиля оптической поверхности h(x,y) контролируемой детали формируется приёмной апертурой диаметром A в плоскости матричного приёмника излучения с количеством чувствительных элементов M*N по горизонтали и вертикали соответственно, после чего реализуется операция дискретного БПФ карты профиля оптической поверхности:
M-1N-1
х , х ,, ( -2nin ) ( -2nim
К„ = L L hmn exP I —— I exp I —— fo ^ I N ) У M
(5)
где М, Ы, т, п - максимальное количество значимых точек по горизонтальной и вертикальной осям соответственно и их координаты; 1гтп - дискретное представление карты профиля оптической поверхности.
Полученный спектр 1гтп регистрируемой карты профиля формируется с учётом линейных размеров исследуемой области детали, так что нормированный спектр имеет вид:
hH = h (A/M) • ( A / N), mn mn
(6)
где A/M, A/N - шаги дискретизации с учётом линейного увеличения приёмного объектива по горизонтальной и вертикальной осям соответственно.
На втором этапе вычисляется функция S2D(vx, vy) как квадрат модуля двумерного спектра карты профиля, приведённого к единичной области (площади) исследования в соответствии с (3):
S 2D = S-1
hH mn
(7)
где = (Сх-Су)-С - площадь контролируемой (исследуемой) детали; бх = 1 / Мтах, Су = 1 / Утах - шаги дискретизации, ограничивающие максимальные пространственные частоты Мтах=М/ (2А), Утах = Ы/(2А) вдоль горизонтальной и вертикальной осей матрицы ; С - число элементов матрицы, содержащей информацию о профиле волнового фронта, отражённого от
поверхности плоского оптического элемента (информативная матрица).
При этом матрица суммарных пространственных частот определяется выражением
Vp mn = ï](Um )2 + (Vn )2
(8)
где Um = (m-M+1/2) /A, vn = (n-N+1/2) / A - векторы частотных сеток по горизонтальной и вертикальной осям соответственно.
В полученной матрице частот последовательно рассматривается каждая суммарная пространственная частота (от меньшей к большей до частоты Найкви-ста, соответствующей радиусу максимальной окружности (карты профиля оптической поверхности hmn диаметром D = 2^S/п), вписанной в матрицу фотоприёмного устройства). В итоге определяются значения векторов суммарных пространственных частот по различным направлениям vpmn и соответствующие им среднеквадратические значения векторов S1D(v), которые представляют искомую целевую функцию, приведённую к одномерному виду.
На заключительном этапе в соответствии с (4) находят СКО искажений волновых фронтов при контроле локальных отклонений изотропной поверхности [3, 6] в заданном частотном диапазоне исследования vP mm-vp max, проводя сравнительную оценку с данными, представленными в табл. 1:
rmssid (vp) =
2p
К",ЫЛ _
i S1D (Vp)dVp
Vp min=0
(9)
где vp min, vp max - радиальные значения минимальной и максимальной пространственных частот диапазона исследования искажений волнового фронта при контроле оптической детали плоской формы.
В рамках пятого этапа представленного выше алгоритма определения оптимизационных параметров целевой функции СПКФ1 проводится оценка влияния математического аппарата дискретного БПФ на искомую величину СКО. В результате проведённых исследований показано, что функция невязки, построенная на основе сравнения математической модели исходной карты профиля оптической поверхности hmn с восстановленной картой h'mn, в результате применения операторов прямого и обратного БПФ составила величину порядка 10-16. Полученная величина свидетельствует о незначительном вкладе влияния математического аппарата дискретного БПФ при обработке изотропных полированных поверхностей. Тем не менее, в задачах диагностики сложнопрофильных изделий данный эффект может быть учтён как одна из составляющих систематических видов погрешностей.
Существенным фактором, влияющим на методическую погрешность определения величины СКО во втором спектральном диапазоне при построении одномерной целевой функции S1D(v]5), является эффект спектральной «утечки» частоты (рис. 3а, б).
2
И(х),нм\ W(x)
h(x), нм
MÚ
а)
J 0\1 У \1 \1 \ S(x)=W(x)-Sn(x)
_ V
ААААА х,мкм
S„(v,f)=3(vY-v0x) S(yx)
S(Vx-Vl;г) \ Sn(\x-v0x)
S(vr-vlr) '
vlx
t
vOx
VX, 1/мкм
Z{Sn(x)}=S0(Vx) 3 {W(x)}=W(vx)
>S(vx)=fV(vx)®S0(vx) ¡
h(x), нм
KM S0M
3{So(x)}=So(vx) ^Г
3{Щх)}=Щух) [
S (vx-vlx) °
S(vx-vlx)
S(y¿ Sf>(v y)=5(vy-vOy)
HM
б)
A S(x)=W(x)-Sn(x)
А Л / \ Л Л х,мкм
У V V V \
\x. 1/mkm
Рис. 3. Иллюстрация анализа влияния эффекта спектральной «утечки» частоты на значение одномерной целевой функции БШ(vр) в заданном пространственно-частотном диапазоне: а) при использовании прямоугольной (оконной) функции выборки; б) при использовании оптимальных оконных функций с «гладкой» огибающей
Как можно видеть из рис. 3а, он проявляется в переналожении пространственных частот первого и второго диапазонов, в результате чего амплитудные значения гармоник низких пространственных частот в значительной степени превышают амплитудные значения гармоник высоких пространственных частот, что приводит к значительной погрешности вычисления одномерной целевой функции S1D(vp) и,
как следствие, измерения величины СКО (10). Данный эффект проанализирован в работе при обработке центральных сечений модельных карт h(x) = S0(x) (рис. 3а). Как известно, преобразование Фурье определено для бесконечной функции S0(x), однако при использовании реальной реализации как ограниченной (финитной) функции S(x) вследствие заданных пространственно-частотных диапазонов (табл. 1) анализируемый спектр будет равен «свёртке» спектров исходной функции с функцией выборки:
S (vx) = W (vx) ® S0(vx), где S (vx) - спектр исследуемого распределения центрального сечения модельной карты в заданном пространственно-частотном
диапазоне vmm < vx< vmax (табл. 1); S0 (vX) - спектр исходного распределения в широком пространственно-частотном диапазоне - ^ < vx< W (vx) - прямоугольная (оконная) функция выборки, в общем случае определяемая как rect(x) исследуемого масштаба неод-нородностей эквивалентного соответствующему пространственно-частотному диапазону (табл. 1).
В результате вместо ожидаемой 5-функции Дирака S0 (vx) = 5 (vx-v0X) результирующий спектр представляется в виде функции (рис. 3а)
sin с (vx — v0 x) = So(v x — v0 x). (10)
При этом амплитудные значения гармоники низкой пространственной частоты v0x первого частотно-
го диапазона в значительной степени превышают амплитудные значения гармоник высокой пространственной частоты у1х второго частотного диапазона.
Для оценки влияния данного эффекта на результирующую погрешность измерений была проведена серия расчётов СКО по 4 реализациям топографических карт искажений волнового фронта для локальных отклонений, отнесённых ко второму пространственно-частотному диапазону (табл. 2).
Табл. 2. Анализ влияния низких пространственных частот на определение СКО высокочастотных неоднородностей
№ реализации Значения параметра СКО с учётом влияния низких пространственных частот Значения параметра СКО после оконной фильтрации низких пространственных частот
1 6,045 нм 4,863 нм
2 7,149 нм 5,967 нм
3 16,176 нм 8,387 нм
4 14,545 нм 9,789 нм
Анализ представленных в табл. 2 результатов позволяет сделать вывод о необходимости предварительной пространственно-частотной фильтрации гармоник первого диапазона с целью минимизации погрешности измерения величины СКО высокочастотных пространственных неоднородностей второго диапазона (табл. 1). Проведённые в работе исследования свидетельствуют также о том, что минимизация влияния эффекта спектральной «утечки» частоты на результирующую погрешность определения величины СКО может быть устранена путём подбора оптимальных оконных фильтров, имеющих более гладкую огибающую функцию [4 - 6, 12 - 14], например, оконных фильтров Ханна или Уолша [5, 12, 13], использованных в данной работе (рис. 3б).
5. Анализ результатов аттестационного контроля СКО неоднородностей искажений волнового фронта в лабораторных и производственных условиях в пространственно-частотных диапазонах I, II
С целью апробации представленного алгоритма было разработано авторское программное обеспечение «Collapser» (рис. 4), позволяющее проводить аттестационный контроль качества оптических деталей на базе измерений, полученных на измерительных
стендах [5] (рис. 2). Разработанное программное обеспечение было протестировано на модельных топографических картах (рис. 5а, б, г), созданных в среде «МаНаЬ» с заданным СКО неоднородностей искажений волновых фронтов, а также на результатах экспериментальных измерений, полученных в лабораторных (рис. 6) и производственных (рис. 5в, г) условиях на АО «Лыткаринский завод оптического стекла» на базе измерительных стендов [5] (рис. 2) соответственно.
а)
11° Таблица СКО
ни
Диапазон частот, 1/мм Полное CKD. нм CKD над критечесой линией, нм
► 0-Q.358S 32,0499 31,9737
0,3586-0.3605 0.1257 0
0.3605-0.3665 0.2406 0,0711
0.3665-0.3686 0.1323 0 _
0,3686-0.3726 0.1923 0,0479
0.3726-0.3785 0.2405 0
0.3785-0.3826 0.1806 0,0419
0,3826-0.3866 0.1944 0
очябб-птегк Л 17В? ЛШЯ5
в)
107 106 105 ю4 103 ю2 10> 100 10-'
SlD(vr), нм2-мм
! 1
i i ^Оч I 1 I >П ' диапазон \ диапазон' диапазон \г, 1/мм
29,2
27.2
25.3
23.3
21.4
19.4
17.5
15.6
13.6
11.7 О
SHXyr), нм2-мм
диапазон i
диапазон i
V,, 1/мм
0,01
0,18
2,40 6,00
0,01 0,18 2,40 6,00 г)
Рис. 4. Внешний вид разработанного программного обеспечения «СоНарзвг»: а) таблица СКО, б) обрабатываемая карта профиля оптической поверхности и расчётная функция б2п(ух, Уу); в) целевая функция БЮ( v); г) графическое представление разброса оптимизационного параметра - СКО в пространственно-частотных диапазонах
Для минимизации краевых эффектов, обусловленных спецификой математического аппарата дискретного БПФ, были использованы непрямоугольные оконные фильтры Ханна и Уолша [5, 12, 13]. По результатам обработки построенных графических зависимостей одномерных целевых функций S1D(vp) для
пространственных частот остаточных неоднородностей (рис. 5г, 6) можно определить величину их СКО в диапазонах I, II как для деталей в диаметре менее 100 мм, так и для изделий со световой апертурой до 1000 мм в соответствии с (10).
В табл. 3 приведены сравнительные значения СКО модельных топографических карт RMSm, рассчитанные в среде «Matlab» и полученные в результате расчёта по разработанному алгоритму RMS*m. Расхождение между расчётными значениями, представленное в виде относительной среднеквадратической погрешности Dm =sj(RMSm )2 - (RMS *m)2 /RMSm, не
превышает 2 %, как следует из табл. 3.
Полученные и приведённые в табл. 3 результаты позволяют судить о корректности рассмотренных в
работе научных положений по разработке алгоритма обработки экспериментальных результатов, полученных при контроле искажений волновых фронтов для плоских оптических элементов.
Табл. 3. Результаты расчёта СКО неоднородностей для модельных топографических карт профилей оптической поверхности в пространственно-частотном диапазоне 1-11
Модель RMSm, нм RMS*m, нм Am,%
1 15,454 15,452 1,6087
2 15,476 15,479 1,9690
3 14,346 14,349 2
В табл. 4 представлены значения СКО экспериментальных данных, полученных с помощью лабораторного [5] и производственного (рис. 2) измерительных стендов, а также разработанного алгоритма - соответственно RMSe, RMS*e, а также относительные среднеквадратические погрешности
De =7(RMSe )2 - (RMS *e )2 /RMSe
Из табл. 4 можно видеть, что максимальная относительная погрешность определения величины СКО не-
однородностей искажений волновых фронтов при контроле оптических деталей со световой апертурой 100 мм составила величину не более 8 %, а при контроле крупногабаритных оптических изделий диамет-
ром до 1000 мм - не более 13 %. Различие между полученными результатами обусловлено спецификой методики и аппаратуры контроля, в частности особенностью юстировки зеркально-линзового коллиматора.
А {х,у)
а)
б)
в)
# С»1
! * * Л
ф о
107 106 ю5 104 103 ю2 10' 100 Ю-'
- ЩуД нм2-мм
1
\
2 \\ „ Г^н
диапазон диапазон \ диапазон \ >
г)
0,01
0,18
2,4
\г,1/мм
Рис. 5. Результаты апробации разработанного алгоритма расчёта СКО неоднородностей искажений волновых фронтов для топографических карт h(x, у) крупногабаритных оптических деталей диаметром до 1000 мм в I и II спектральных диапазонах: модельные топографические карты h(x, у), построенные по степенному закону и закону Лоренца, распределения неоднородностей искажений волновых фронтов и соответствующие им целевые функции б2п(ух, Уу) (а, б); экспериментально зарегистрированная топографическая карта h(x, у) и соответствующая ей двумерная целевая функция
б20(ух, Уу) (в); численные результаты построения и расчёта одномерных целевых функций БШ(ур) для модельных (1, 2)
и экспериментальной (3) топографических карт (г)
ственной частотой ур т„ в пределах заданной области (апертуры) исследования оптического элемента детали по виду характерного пика на расчётной кривой
£Ш(Ур тп) (рис. 5г, 6).
Табл. 4. Результаты расчёта СКО искажений волновых фронтов неоднородностей экспериментальных топографических карт, контролируемых плоских оптических элементов диаметрами 100 мм и 1000 мм в соответствующих пространственно-частотных диапазонах 1-11
1 >Б1Щуг), ^ нм2-мм
1000 \ ИМ$=18,57 нм ЯМ$=4,341 нм
100 Локальный '—\ поверхностный дефект
10 \ ^........5/Ояо(у)
1 0 диапазон - диапазон \ диапазон]^. ,. А ^ Ч
0,01 0,18 2,4 б
Рис. 6. Результаты расчёта разработанного алгоритма определения СКО искажений волновых фронтов при контроле локальных отклонений профиля плоской эталонной поверхности диаметром 100 мм в двух пространственно-частотных диапазонах
С использованием табл. 3, 4 можно оценить как численное значение СКО пространственных неодно-родностей волновых фронтов (10) в результате оптимизации целевой функции £Ш(ур тт), так и определённый тип локального отклонения поверхности от заданного профиля с соответствующей простран-
Реализация ЕЫБв, нм ШБ*е, нм Де, %
Световая апертура исследуемой поверхности плоской оптической детали 100 мм
I диапазон (0,01- 0,18), мм-1 18,53 18,57 6,5
II диапазон (0,18- 2,4), мм-1 4,345 4,341 7,9
Световая апертура исследуемой поверхности плоской оптической детали 1000 мм
I диапазон (0,001- 0,03), мм-1 31,76 31,77 2,5
II диапазон (0,03- 0,4), мм-1 5,21 5,26 13
Приведённые в табл. 4 данные были получены при помощи измерительных стендов диагностики неоднородностей искажений волновых фронтов для контроля плоских оптических элементов диаметром 100 мм и до 1000 мм. Фотография полноапертурного измерительного стенда представлена на рис. 7.
Заключение
Таким образом, разработанный метод динамической интерферометрии контроля неоднородностей искажений волновых фронтов для плоских оптических элементов позволил в лабораторных и производственных условиях с высокой точностью определить наличие и размер неоднородностей в частотных диапазонах I, II, для которых значения СКО составили: для диаметра 100 мм - 18,570; 4,341 нм соответственно, а для диаметра 1000 мм - 31,77; 5,26 нм соответственно.
Рис. 7. Фотография производственного измерительного стенда на базе динамического интерферометра «1п1еШишИ2000» и зеркально-линзового коллиматора для апробации разработанного алгоритма определения СКО искажений волновых фронтов при контроле локальных отклонений профиля плоской эталонной поверхности диаметром до 1000 мм
Достигнутая при оптимизации целевой функции величина максимальной неисключённой систематической погрешности алгоритма для первого и второго частотных диапазонов, в соответствии с установленными требованиями, составила величину порядка 2 %, что вполне удовлетворяет указанным точностным требованиям. Разработанный метод в перспективе позволит проводить оперативный аттестационный контроль локальных отклонений субнанометрового уровня профилей крупногабаритных оптических изделий в условиях производственных вибраций.
Благодарности
Автор выражает глубокую благодарность сотрудникам научно-производственного комплекса «Астрономическая и космическая оптика» АО «Лыткарин-ский завод оптического стекла» за помощь в проведении экспериментов и компетентных консультациях.
Отдельные параграфы работы выполнены в рамках финансирования, осуществляемого согласно назначению стипендии Президента Российской Федерации молодым учёным и аспирантам, осуществляющим перспективные научные исследования и разработки по приоритетным направлениям модернизации российской экономики на 2013-2015 годы. При-
каз № 136 Министерства образования и науки Российской Федерации от 28 февраля 2013 года.
Литература
1. Абдулкадыров, М.А. Технология изготовления высокоточных крупногабаритных облегчённых асферических зеркал с высокой стабильностью формы поверхности / М.А. Абдулкадыров, Н.С. Добриков, А.П. Патрикеев, В.Е. Патрикеев, А.П. Семёнов // Оптический журнал. - 2014. - Т. 81, № 12. - С. 6-15.
2. Абдулкадыров, М.А. Современные способы изготовления астрономических и космических зеркал / М. Абдулкадыров, А. Семёнов // Фотоника. - 2015. -№ 3. - С. 62-79.
3. Денисов, Д.Г. Разработка методов и аппаратуры лазерного интерференционного контроля формы и качества оптических поверхностей крупногабаритных зеркал на стадиях шлифования: дис. ... канд. техн. наук / Денисов Дмитрий Геннадьевич. - МГТУ им. Н.Э. Баумана, Москва, 2010.
4. Campbell, J.H. NIF optical materials and fabrication technologies: an overview / J.H. Campbell, R.A. Hawley-Fedder, Ch.J. Stolz, J.A. Menapace, M.R. Borden, P.K. Whitman, J. Yu, M.J. Runkel, M.O. Riley, M.D. Feit, R.P. Hackel // Proceedings of SPIE. - 2004. - Vol. 5341. -DOI: 10.1117/12.538462.
5. Denisov, D.G. Method for certification monitoring of surface inhomogeneities of optics based on frequency analysis of the surface profile / D.G. Denisov, N.V. Baryshnikov, Ya.V. Gla-dysheva, V.E. Karasik, A.B. Morozov, V.E. Patrikeev // Measurement Techniques. - 2017. - Vol. 60, Issue 2. -P. 121-127. - DOI: 10.1007/s11018-017-1160-0.
6. Денисов, Д.Г. Измерение параметров микронеровностей крупногабаритных шлифованных поверхностей оптических деталей при помощи лазерной интерферометрии / Д.Г. Денисов, В.Е. Карасик, В.М. Орлов // Метрология. - 2009. - № 9. - С. 15-24.
7. Nikitin, А. A device based on the Shack-Hartmann wave front sensor for testing wide aperture optics / A. Nikitin, J. Sheldakova, A. Kudryashov, G. Borsoni, D. Denisov, V. Karasik, A. Sakharov // Proceedings of SPIE. - 2016. -Vol. 9754. - 97540K. - DOI: 10.1117/12.2219282.
8. Nikitin, А. Hartmannometer versus Fizeau interferometer: advantages and drawbacks / A. Nikitin, J. Sheldakova,
A. Kudryashov, D. Denisov, V. Karasik, A. Sakharov // Proceedings of SPIE. - 2015. - Vol. 9369. - 936905. -DOI: 10.1117/12.2085263.
9. Барышников, Н.В. Метод и аппаратура аттестационного контроля радиусов кривизны сферических поверхностей оптических изделий при помощи датчика волнового фронта / Н. В. Барышников, Д. Г. Денисов, В. Е. Ка-расик, А.А. Сахаров // V Международная конференция по фотонике и информационной оптике: Сборник научных трудов. - 2016. - С. 416-417.
10. Полещук, А.Г. Лазерные интерферометры для контроля формы оптических поверхностей / А. Г. Полещук,
B.Н. Хомутов, А.Е. Маточкин, Р.К. Насыров, В.В. Чер-кашин // Фотоника. - 2016. - № 4. - С. 38-51. - DOI: 10.22184/1993-7296.2016.58.4.38.50.
11. Полещук, А.Г. Методы оперативного контроля характеристик дифракционных оптических элементов в процессе изготовления / А.Г. Полещук, В.П. Корольков, Р.К. Насыров, В.Н. Хомутов, А.С. Конченко // Компьютерная оптика. - 2016. - Т. 40, № 6. - С. 818-829. - DOI: 10.18287/2412-6179-2016-40-6-818-829.
12. Sidick, E. Power spectral density specification and analysis of large optical surfaces / E. Sidick // Proceedings of SPIE. - 2009. - Vol. 7390. - 73900L. - DOI: 10.1117/12.823844.
13. Alcock, S.G. Using the power spectral density method to characterise the surface topography of optical surfaces / S.G. Alcock, G.D. Ludbrook, T. Owen, R. Dockree // Proceedings of SPIE. - 2010. - Vol. 7801. - 780108 (8 p.). -DOI: 10.1117/12.861539.
14. Методы компьютерной оптики / А.В. Волков, Д.Л. Го-ловашкин, Л.Д. Досколович, Н.Л. Казанский, В.В. Кот-
ляр, В.С. Павельев, Р.В. Скиданов, В.А. Сойфер, В.С. Соловьев, Г.В. Успленьев, С.И. Харитонов, С.Н. Хонина; под ред. В.А. Сойфера. - Изд. 2-е, испр. - М.: Физмат-лит, 2003. - 688 с. - ISBN: 5-9221-0434-9.
15. ISO 10110-6:2015. Optics and photonics - Preparation of drawings for optical elements and systems - Part 6: Centring tolerances. - 2nd Ed. - Vernier, Geneva, Switzerland: ISO, 2015.
Сведения об авторе
Денисов Дмитрий Геннадьевич, 1981 г. рождения, окончил с отличием Московский государственный технический университет им. Н.Э. Баумана по специальности «Инженер-разработчик лазерной техники и лазерных технологий» в 2006 г., защитил кандидатскую диссертацию на тему «Разработка методов и аппаратуры лазерного интерференционного контроля формы и качества оптических поверхностей крупногабаритных зеркал на стадиях шлифования» в 2010 г., защитил диплом по специальности «Преподаватель высшей школы» в 2013 г. Работает доцентом на кафедре лазерных и оптико-электронных систем МГТУ им. Н.Э. Баумана. Область научных интересов: интерференционные методы и приборы высокоточного контроля оптических поверхностей, современные разделы оптического материаловедения, оптические материалы и технологии, методы и устройства адаптивной оптики в лазерных технологиях, системы управления параметрами лазерного излучения.
ГРНТИ 59.45.37.
Поступила в редакцию 18 июля 2017 г. Окончательный вариант - 7 ноября 2017 г.
ERROR ANALYSIS IN DIGITAL PROCESSING OF THE RESULTS OF INTERFEROMETRIC CONTROL OF NANO-SCALE LOCAL DEVIATIONS OF OPTICAL SURFACES
D.G. Denisov1
'Bauman Moscow State Technical University, Moscow, Russia Abstract
A dynamic interferometry method for controlling nano-scale local deviations of optical surfaces from a prescribed profile is developed, theoretically substantiated, and experimentally verified on the basis of an algorithm for calculating an objective function, defined as the spectral density of a one-dimensional correlation function. Theoretical and experimental approaches to determining the root-mean-square error when finding local deviations of the optical surfaces for the element diameters ranging from 100 mm to 1,000 mm are presented, while taking into account the non-excluded systematic and random error components in determining the objective function.
Keywords: optical control, interferometry, surface measurements, power spectral density, method errors, edge effect, the effect of "leakage" of frequency.
Citation: Denisov DG. Error analysis in digital processing of the results of interferometric control of nano-scale local deviations of optical surfaces. Computer Optics 2017; 41(6): 820-830. DOI: 10.18287/2412-6179-2017-41 -6-820-830.
Acknowledgments: The author expresses his deep gratitude to the employees of the Scientific and Production Complex "Astronomical and Space Optic'' of "JSC Lytkarino Factory of Optical Glass" for their assistance in conducting the experiments and competent consultations.
Separate sections of the work were carried out under the scholarship of the President of the Russian Federation, awarded to young scientists and graduate students doing research in advanced areas of modernization of the Russian economy in 2013-2015.
References
[1] Abdulkadyrov MA, Dobrikov NS, Patrikeev AP, Patrikeev VE, Semenov AP. Technology for fabricating large, high-accuracy, lightened aspheric mirrors with high stability of the surface shape. J Opt Techn 2014; 81(12): 706-713. DOI: 10.1364/JOT.81.000706.
[2] Abdulkadyrov MA. Modern ways of production of astronomical and space mirrors [In Russian]. Photonics 2015; 3: 62-79.
[3] Denisov DG. Development of laser interferometric methods and equipment to control the profile and smoothness of the surfaces of large-size optical mirrors at a polishing
stage. Thesis for the Candidate degree in Technical Sciences. Moscow, 2010.
[4] Campbell JH, Hawley-Fedder RA, Stolz ChJ, Menapace JA, Borden MR, Whitman PK, Yu J, Runkel MJ, Riley MO, Feit MD, Hackel RP. NIF optical materials and fabrication technologies: an overview. Proc SPIE 2004; 5341. DOI: 10.1117/12.538462.
[5] Denisov DG, Baryshnikov NV, Gladysheva YaV, Karasik VE, Morozov AB, Patrikeev VE. Method for certification monitoring of surface inhomogeneities of optics based on frequency analysis of the surface profile. Meas Tech 2017; 60(2): 121-127. DOI: 10.1007/s11018-017-1160-0.
[6] Denisov DG, Karasik VE, Orlov VM. Measurement of the microroughness parameters of large-size polished optics surfaces using laser interferometry [In Russian.] Metrology 2009; 9: 15-24.
[7] Nikitin A, Sheldakova J, Kudryashov A, Borsoni G, Denisov D, Karasik V, Sakharov A. A device based on the Shack-Hartmann wave front sensor for testing wide aperture optics. Proc SPIE 2016; 9754: 97540K. DOI: 10.1117/12.2219282.
[8] Nikitin A, Sheldakova J, Kudryashov A, Denisov D, Karasik V, Sakharov A. Hartmannometer versus Fizeau interferometer: advantages and drawbacks. Proc SPIE 2015; 9369: 936905. DOI: 10.1117/12.2085263.
[9] Baryshnikov NV, Denisov DG, Karasik VE, Sakharov AA. Method and equipment for attestation for control of the radius curvature of spherical surfaces of optical details by using a wave front sensor. Proc V International Conference on Photonics and Information optics 2016: 416-417.
[10] Poleshchuk AG, Khomutov VN, Matochkin AE, Nasyrov RK, Cherkashin VV. Laser interferometers for optical surfaces testing [In Russian]. Photonics 2016; 4: 38-51. DOI: 10.22184/1993-7296.2016.58.4.38.50.
Author's information
Dmitry Gennadievich Denisov, 1981, graduated with honors from the Moscow State Technical University named after N.E. Bauman in the field of engineer-developer of laser technology and laser technologies in 2006, defended his thesis on "Development of methods and equipment for laser interference control of the shape and quality of optical surfaces of large-sized mirrors at grinding stages" in 2010, defended a diploma in the specialty "Teacher of higher education" in 2013. Works as an assistant professor of Laser and Optoelectronic Systems department at the Bauman Moscow State Technical University. Scientific interests: interference methods and devices for high-precision control of optical surfaces, modern sections of optical material science, optical materials and technologies, methods and devices of adaptive optics in laser technologies, control systems for laser radiation parameters.
Received July 18, 2017. The final version - November 7, 2017.
[11] Poleshchuk AG, Korolkov VP, Nasyrov RK, Khomutov VN, Konchenko AS. Methods for on-line testing of characteristics of diffractive and conformal optical elements during the manufacturing process [In Russian]. Computer Optics 2016; 40(6): 818-829. DOI: 10.18287/2412-61792016-40-6-818-829.
[12] Sidick E. Power spectral density specification and analysis of large optical surfaces. Proc SPIE 2009; 7390: 73900L. DOI: 10.1117/12.823844.
[13] Alcock SG, Ludbrook GD, Owen T, Dockree R. Using the power spectral density method to characterise the surface topography of optical surfaces. Proc SPIE 2010; 7801: 780108. DOI: 10.1117/12.861539.
[14] Soifer VA, ed. Methods for computer design of diffractive optical elements. New York: John Wiley & Sons, Inc.; 2002. ISBN: 978-0-471-09533-0.
[15] ISO 10110-6:2015. Optics and photonics - Preparation of drawings for optical elements and systems - Part 6: Centring tolerances. 2nd Ed. Vernier, Geneva, Switzerland: ISO; 2015.
