CC BY
43
УДК 531.717.81:514.88:51-74
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПОВЕРХНОСТИ, КОНТРОЛИРУЕМОЙ ЭКСПЕРИМЕНТАЛЬНЫМ ЛАЗЕРНЫМ КРУГЛОМЕРОМ
© 2015 О.А. Заякин
Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет)
Статья поступила в редакцию 30.07.2015
Приведена математическая модель поверхности, контролируемой экспериментальным лазерным кругломером. Прибор предназначен для контроля формы зеркально отражающей рабочей поверхности детали подшипника. Контроль проводится путем координатных измерений. Характерная особенность способа измерений - основа на обработке сигнала от зеркальной составляющей отраженного зондирующего лазерного пучка. Модель пригодна для оценки алгоритмической составляющей погрешности метода.
Ключевые слова: триангуляция, измерения геометрических величин, математическая модель, поверхность, кругломер, зеркально отраженный свет, топограмма, статистические характеристики.
Экспериментальный лазерный кругломер (двумерный лазерный триангулятор) [1] был разработан ранее в Самарском филиале Физического института РАН при участии автора. Его работа основана на методе триангуляции с использованием зеркально отраженного света зондирующего источника. Прибор предназначен для контроля круглости рабочих поверхностей внутренних колец подшипников.
Оптико-механическая схема координатных измерений приведена на рис. 1. Источник излучения 1 и каретка с фотоприемником 3 фиксируются перед началом измерений. В процессе измерений контролируемая делать 2 вращается вокруг своей оси на один полный оборот с заданным шагом. При этом контролируемая поверхность детали сканируется по ее радиальному профилю. Координаты отраженного пучка фиксируются на фотоприемнике 2 после каждого шага. Затем так же может быть просканирован другой радиальный профиль, для чего источник излучения 1 смещается на заданное расстояние вдоль оси г. Поворот каретки с фотоприемником 3 нужен только для автоматизации настройки оптической схемы перед началом цикла измерений деталей одного типоразмера. При замене однотипных деталей этого уже не требуется.
Для проверки функции преобразования и оптимизации оптической схемы нужно математическое описание контролируемой рабочей поверхности детали. Все найденные автором в литературных источниках модели [2] являются одномерными. Однако трехмерной модели не было найдено. Она была создана автором самостоятельно.
Заякин Олег Александрович, кандидат технических наук, доцент кафедра информационных систем и технологий. E-mail: [email protected]
Рис. 1. Оптико-механическая схема координатных измерений: 1 - источник излучения; 2 - контролируемая поверхность; 3 - видеокамера; АВ, С направления сканирования;
а, Ь - падающий и отраженный лучи света; у, Р - координаты точки падения отраженного луча на фотоприемнике; ^ - смещение источника излучения; Ь - радиус поворота каретки с видеокамерой
Для получения расчетной функции преобразования модель поверхности должна отвечать требованиям геометрической оптики.
Модель поверхности представлена в данной работе как жесткий сплошной экран с коэффициентом отражения, равным единице. Для анализа функции преобразования требуется еще, чтобы поверхность в модели была гладкой. Это означает, что отражение от каждого ее точечного участка описывается так же, как и отражение луча от плоскости, ориентированной перпендикулярно локальной нормали к моделируемой поверхности в ее освещенной точке. Минимальные воспроизводимые в модели детали рельефа по-
верхности ограничиваются размером светового пятна, которое перемещается по контролируемой поверхности при ее сканировании.
Уравнение восстановленной поверхности Р (г ) = 0 (см. рис. 2):
_ ^ ¿{г) = ¿¿(ГОН ^(Аг), (1)
где Г0 и г - радиус-векторы точки номинальной
и во сстановленной поверхностей; Аг=г — г0; Р0(г0) = 0 - уравнение номинальной поверхности.
Функция Аг(г0), которая выводится из
Р1(Аг) = 0 , определяет отклонения формы реальной поверхности от номинальной формы.
Для анализа удобно представить 2-е слагаемое в (1) в виде явной скалярной функции И(г0). Эта функция ставит в соответствие каждой точке номинальной поверхности значение высоты отклонения формы, отсчитанной от этой точки.
Пусть Аг(г0) не перпендикулярен номинальной поверхности. Расс мотрим точки на поверхности, определяемые г0 и г02 (см. рис. 1). Различие высот в этих двух точках не превосходит величины 0,05 рад. Тогда искомая функция
И (?0) =
(Аг • ?)
(2)
где ? = ?(г0) = ± grad Р0 (имеется ввиду градиент в физических координатах). В системе цилиндрических координат | р, р, 2 } :
grad Р0 =
дР0 1 дР0 дР0
, —х" др р др
где р, р, 2 - соответственно радиусная, аз и-мутальная и осевая координаты. Функция ? (г0) определяет перпендикуляр в каждой точке номинальной поверхности.
Функцию И(г0) назовем «функцией отклонения от номинальной формы». Отсчет Н ведем в радиальном направлении.
Представим И(г0) в виде топограммы.
Для того чтобы моделировать типичную поверхность, нами использованы статистически усредненные спектральные зависимости, характеризующие амплитуды пространственных гармоник профилей этой поверхности. Эти зависимости получают, считая профили поверхности реализациями случайного процесса. Каждая из этих зависимостей относится к множеству однотипных профилей. Она представляет собой функцию спектральной плотности мощности указанного случайного процесса. Преимущества этого подхода в том, что, во-первых, таким путем нетрудно получить замкнутые профили, а во-вторых, из спектрального представления профилей поверхности нетрудно получить ее топограмму.
д2
(3)
Рис. 2. Связь между разными представлениями рельефа поверхности: 1 - номинальная поверхность;
2 - модель реальной поверхности
Опишем, как может быть получен из указанной функции соответствующий ей модельный профиль поверхности. Основываясь на соотношениях между спектральными и статистическими характеристиками шероховатости поверхности, приведенными в [3], считаем, что для этого достаточно выполнить следующие два условия. Во-первых, отсчеты амплитуд гармоник профиля моделируемой поверхности следует расположить равномерно по оси абсцисс и совместить с графиком зависимости, представляющей собой квадратный корень из названной функции спектральной плотности мощности. График данной зависимости при этом служит огибающей дискретного спектра гармоник. И, во-вторых, начальные фазы гармоник не должны коррелировать между собой.
Мы моделируем бездефектную поверхность, то есть поверхность, удовлетворяющую требованиям к ее качеству формы. Числовые коэффициенты в формулах, описывающих такую бездефектную поверхность, получим из экспериментальных графиков.
Для моделирования трехмерной поверхности используем топограмму ее неровностей. По приведенным в научных публикациях [4] топограм-мам реальных рабочих поверхностей довольно сложно получить их математическое описание, так как информация приведена в графическом виде. Статистических зависимостей в аналитическом виде в рассмотренных публикациях не приводится. Для решения этой задачи в данной работе используем следующий подход.
Топограмма моделируемой поверхности должна учитывать статистические характеристики этой поверхности, описывающие ее профили в различных направлениях. Как правило, это два ортогональных направления в криволинейных координатах, определяемых номинальной поверхностью. Будем считать, что эти про-
фили представляют собой реализации случайных процессов, не зависящих друг от друга в том, что касается их вероятностных характеристик. Далее, для моделирования нужной нам топограммы используем понятия «двумерной пространственной частоты», широко применяемое в физической оптике и основанной на нем характеристики -функции двумерной спектральной плотности мощности случайного процесса. Взаимная статистическая независимость двух одномерных случайных процессов позволяет нам представить эту функцию в виде произведения двух функций одномерной спектральной плотности:
/ Щ,Щ )=/; (щ )/2 (Щ), (4)
где щ, (02, - пространственные частоты вдоль профилей, совпадающих с координатными линиями криволинейного базиса, соответственно, I = Яр и X (поперек номинальной поверхности - это осевое сечение).
Отсюда, с помощью преобразования Фурье функции, являющейся квадратным корнем функции (4), получим отсчет высоты микрорельефа в любой нужной точке с координатами I и х. Замкнутость профилей вдоль направляющих линий номинальной поверхности ведет к дискретности спектра частот Щ. Эти частоты должны удовлетворять условию 01 = N Я, где N - количество волн (отклонений от круглости), приходящееся на оборот поверхности. При этом функция двумерной спектральной плотности получается дискретной вдоль оси частот 01 и непрерывной вдоль оси частот 02. На графике она представляет собой набор плоских фигур, параллельных друг другу и оси 02, и перпендикулярных оси щ1. Форма этих фигур такова, что функция /(щ щ2) является для них огибающей.
Примем, что Я не зависит от X, а Я равно его среднему интегральному значению по X. Таким образом, формула для расчета топограммы моделируемой поверхности имеет вид:
н(р, х)=С £ |е(N,02)/ ()• /2(щ)х
N=1 щ2 = щ2т1п
х (^ + а2 Х~12 Щ2 )]®2, (5)
где Т]1 ^), г/2 (щ2) - функции, каждая из которых имеет значение случайной величины с равномерным законом распределения в интервале от 0 до 2ж рад. Это значение фиксировано при заданном значении аргумента. При этом выполняется условие (например, для г]1 ^)):
Г]1 (Л)=Л1 (В) , если А = В ; 7]1 (Л)ф^1 (в) , если Аф В, (6)
где А и В - отвлеченные числа. Фактически, (5) представляет собой дискретное преобразование Фурье. Пример использования этого преобразования приведен в (4).
Минимальная частота 02т1П выбирается по
возможности близкой к нулю, чтобы моделировать поверхность, не замкнутую в направлении X. Функция е (,02) имеет следующие значения:
0, если N= 0 и щ2 = 0 ;
1/2, если N= 0 и щ2 Ф0, или NФ 0 и 02 = 0;
1, если NФ0 и 02 Ф 0.
Нормирующий множитель С определим из уравнения:
^ Н = С £ |/1 (^/2 02 )Щ , (7)
где V - множество гармоник, соответствующее полосе пропускания частотного фильтра при снятии круглограммы; Е - интервал, равный всему диапазону частот щ2 модельной поверхности.
Остановимся на особенностях алгоритма, связанных с дискретностью численных расчетов. Так как моделирование поверхности с использованием формулы (2.5) проводится численно, то гармоники щ2 также, как и щ1 образуют дискретный набор:
ж М
Щ2 = С-' (8)
С1 Х0
где М =1,2,3,.... Значение С1 выбирается в пределах от четырех до пяти. При меньших значениях С1 получается, что высоты микрорельефа
в точках с координатами Х = — Х0 при одинаковых р заметно коррелируют. Это допускать не следует, если надо более близко к реальности смоделировать незамкнутый профиль. С другой стороны, если С1 будет больше этого значения, то тогда амплитуду низших гармоник уже нельзя будет считать единственным показателем их вклада в моделируемый микрорельеф. Имеются ввиду те гармоники, для которых четверть их периода больше, чем длина профиля по оси X в пределах от минус X до X. Тогда разность высот микрорельефа, в том случае, если он будет представлен только одной из этих гармоник, не достигнет значения амплитуды гармоники на всем протяжении профиля от минус Xo до Xo. И тогда вклад этих гармоник придется учитывать со своим весом, меньшим единицы, что несколько усложняет расчеты. Вместе с тем, при значениях множителя С1 из указанного диапазона получится, что во многих случаях спектры частот щ1 и щ2 будут отличаться по густоте заполнения гармониками. Это надо будет учесть при нормировании функций, описывающих огибающие спектра / (щ1) и / Щ2). В данном примере для этого использована относительно несложная оценка, о которой будет сказано ниже. Итак, указанный диапазон значений С1 можно считать оптимальным.
Соответствующая максимальной частоте
а
1тах величина Мтах = [к/4], где К - количество просканированных точек вдоль одного профиля радиального сечения, приходящееся на полный оборот вокруг поверхности. Операция [.. ] означает выделение целой части числа. Такая величина №тах выбирается потому, что при этом чувствительность высших гармоник ряда Фурье к погрешности величин в каждом из К отсчетов разлагаемого в этот ряд профиля находится в приемлемых пределах. В данном случае имеются ввиду такие пределы, при которых высшие гармоники ряда Фурье имеют физический смысл. Это ограничение на Мтах в несколько раз сокращает объем вычислений, необходимых для моделирования в данной задаче. Аналогично, используя (8), получаем
(9)
=
_ 2 Ах
где Ах - шаг дискретизации профиля поверхности вдоль оси X.
Значения К и Ах выбираются в модели так, чтобы шаг дискретизации профилей вдоль обоих направлений р и X не был больше, чем пространственное разрешение «в плане» тех деталей микрорельефа контролируемой поверхности, которые должны восстанавливаться с помощью исследуемого способа. Поясним, что речь здесь идет о размерах микронеровностей в направлениях, параллельных касательной плоскости к номинальной поверхности на ее локальном участке.
Рассмотрим для примера дорожку качения шарикоподшипника типа 201. Получим математическую модель низкочастотных составляющих микрорельефа поверхности в виде топограммы. Пусть имеем К= 600. Шаги дискретизации в направлениях р и X выберем так, чтобы на криволинейно поверхности они были по возможности одинаковыми. Итак, пусть
Ах = А/(х)=2 к-Щ)/К . Тогда, учитывая геометрические параметры номинальной поверхности: Я(х)=8,65 мм, х0 = 2,36 мм, имеем
N =150 , M = 65 • Черта над символом
max max А
означает интегральное усреднение. Координата, по которой проводится усреднение, записана в скобках.
Получим функцию двумерной спектральной плотности мощности, нужную для нахождения амплитуд пространственных гармоник. Для этого воспользуемся табл. 1. В ней приведены статистические характеристики различных случайных функций. В каждой отдельно взятой строке таблицы приведены данные, относящиеся к какой-либо одной из этих функций.
В табл. 1 приведены некоторые аналитические функции из тех, которые могут быть использованы для описания низкочастотных составляющих профилей микрорельефа рабочих поверхностей изделий машиностроения.
Сопоставим эти функции с экспериментальными спектральными зависимостями, приведенными в литературных источниках. Воспользуемся, например, графиком из [4], с. 216, где показано спектральное представление профиля шлифованной поверхности кольца подшипника качения. Из этих данных получаем следующие оценки для амплитуд H(N) гармоник профиля модельной поверхности с числом волн N на оборот этой поверхности: Я(2)=1 мкм, Н(4)=0,35 мкм, Я(16)=0,15 мкм. Амплитуда гармоники от-считывается от среднего интегрального значения радиального профиля.
Из этого сопоставления можно видеть, что ни одной из формул, приведенных в таб. 1, не удается с приемлемой точностью описать экспериментальный спектр во всем диапазоне его частот. Для повышения точности аппроксимации поступим следующим образом. Разделим весь диапазон частот этой спектральной зависимости на две части, в пределах каждой из них аппроксимируем ее отдельно, с учетом условия «сшивки» этих двух аппроксимирующих функций между собой на границе между двумя частями частотного диапазона. Границу раздела нетрудно определить по характерным различиям в поведении экспериментального графика на разных участках спек-
Таблица 1. Статистические характеристики различных случайных функций
Автоковариационная функция, аргумент т, т> 0 Функция спектральной плотности мощности, аргумент ю, 0)> 0
exp(-r/ a) 1 1 + a2 со2
exp(-T2/ a2) / exp ' о2 a2 ^ , 4 ,
1 -г/a , если т< a , 0 , если т> a sin2 Г1 1 - a ю\ 12 ) (О2
тра. Эта граница находится в пределах величин N от четырех до шести. Для упрощения расчетов примем, что это значение равно четырем, тогда оно будет согласовано с границей спектрального диапазона гранности.
Для описания спектра на участке верхних частот выберем функцию из 1-й строки таб. 1, так как из всех трех функций она наиболее близка к экспериментальному графику при N> 4. Чтобы описать спектр на другом участке, в области нижних частот, используем функцию из 3-й строки табл. 1. В пределы этого участка попадает только несколько старших гармоник. Поэтому, учитывая также «статистический разброс» экспериментальных данных, допустимо выбрать любую функцию из таблицы. При выборе мы учитываем типичный вид автоковариационной функции для низкочастотных составляющих профилей микрорельефа рассматриваемых объектов.
Используя метод наименьших квадратов, получаем для функции ¡1 (N) следующее выражение:
¿И=с,
Ф)=с1г
N
ан щ
N '
'Ш. 1
, если N= 1,2,3
(10)
1 + а,,
^ I
,если N = 4,5,6,...,150 , (11)
где С1Н =1
67740-10~6 мм2,
С1В =1, 7408-10 7 мм2, а1В =1, 403 мм
а1Н =10, 496 мм ,
Функцию ¡2 (М) опишем таким же выражением, как и функцию :
/2 М) = С 2
/2 (М)=С,
а. н М
2 2Н
2 С0 Х0
а2Н М
2 С0 Х0 1
1 + а,,
жМ
С0 Х0
если М= 1, (12)
■ , если М=2,3,4,...,65. (13)
В условиях недостатка опытных данных и отсутствия известных нам теоретических предпосылок примем, что а2н = ан и а2В = а1В .
В силу того что разность частот между соседними гармониками различается в направлениях осей р и X, будут отличаться и нормирующие множители С1Н и С2Н, а также С1В и С2В . Значения неизвестных коэффициентов С2Н и С2В найдем путем сопоставления выражений (10) - (13) друг с другом. При этом считаем, что статистическая оценка высот профилей поверхности в обоих рассматриваемых направлениях имеет одинаковые численные значения. В качестве такой статистической оценки используем
квадратичную сумму амплитуд гармоник профиля поверхности. Суммирование производим по гармоникам, которые попадают в одинаковые частотные диапазоны, занимаемые профилями модельной поверхности в направлениях р и X. Подчеркнем, что речь в данном случае идет о модельных, а не полученных экспериментально профилях поверхности. Для повышения точности оценки неизвестных значений коэффициентов используем как можно более широкий спектральный диапазон модельных профилей. Получаем С2Н = 5,9492-10 6 мм2, С2В = 3,9964-10 7 мм2.
Следующий шаг - получение выражения для огибающей функции двумерной спектральной плотности мощности случайного процесса, моделирующего топограмму интересующих нас микронеровностей. Это выражение имеет следующий вид, как следует из (4):
/(, М)=С g(, М) / () /2 (М), (14) где функция g(, М) аналогична функции g(, 02) в формуле (5).
Нормирующий множитель С, имеющий размерность [ мм 2 ], определим из уравнения (7). В левой части этого уравнения должно стоять известное значение (или оценка) амплитуды профиля реальной поверхности, с учетом ширины частотного диапазона модельной поверхности. Амплитуда профиля определяется как сумма отклонений в обе стороны от средней линии этого профиля.
Поскольку цель состоит в моделировании бездефектной поверхности, то график экспериментального профиля, уже использованный нами здесь, для этого не подходит, так как в [4], где он приведен, не содержится заключение о качестве поверхности, с которой был снят этот профиль. Вместо этого используем результаты координатных измерений колец шарикоподшипников с бездефектной поверхностью. Для подстановки числового значения Ь в (7) выберем высоту гран-ности, а не волнистости, так как использованный выше график спектра профиля поверхности из [4] полностью включает в себя диапазон частот дефектов этого вида, в то же время частотный диапазон волнистости представлен там лишь частично.
Для оценки гранности в рассматриваемой задаче мы использовали результаты измерений контактным кругломером Та1угоп^50, полученные на Самарском подшипниковом заводе СПЗ-4. Измеренная высота гранности, по этим результатам, составила 0,13 мкм.
В соответствии с [5], а также исходя из представления одномерного профиля спектром двумерных частот, получаем из (7) следующее соотношение:
(2х0,25Н!р) = С £
/1N £/2(м) gN,М)
, (15)
2
х
а
2
х
X
2
X
где
Игп= 0,13 мкм . Отсюда получаем
С = 415, Збб мм при ранее найденных зна чениях коэффициентов, входящих в функции
Ь (я) и ь (м).
С использованием полученных результатов была составлена компьютерная программа на MathCad 4.0. С ее помощью была получена модельная топограмма рассматриваемого объекта. Топограмма представлена в виде двумерного массива отсчетов высот микрорельефа поверхности. Массив имел размерность 600 отчетов, вдоль развертки кольца, на 40 отсчетов, в поперечном направлении.
Пример моделирования поверхности приведен на рис. 3 в). Для сравнения на рис. 3 а) и б) приведены микротопограммы реальных поверхностей, номинальная форма которых является плоской. Они были получены с помощью микроинтерферометра, согласно [4]. Сравнение их с результатами нашего моделирования показало, что модель достаточно верно представляет типичную обработанную металлическую поверхность детали машиностроения после ее финишной (чистовой) обработки - шлифования и полирования. Отметим, что на рис. 3 представлены бездефектные поверхности, то есть такие,
которые удовлетворяют требованиям к допустимым отклонениям своей формы.
Вместе с тем видно, что для верного представления сильно направленных следов обработки, характерного, кстати, для рабочих поверхностей деталей подшипников, модель нужно дополнить масштабированием в одном из выделенных ортогональных направлений на топограмме.
Также можно видеть на глаз некоторую «искуссвенность» модельного микрорельефа, которая заключается в ярко выраженном прямоугольном рисунке представленного фрагмента поверхности. Это можно объяснить тем, что при моделировании мы выбрали примерно одинаковый шаг дискретизации функции спектральной плотности мощности по обеим направлениям пространственных частот со1 и с2. Но если шаг для с определялся условием замкнутости радиального профиля, то для с2 такого условия нет. Такой выбор экономит ресурсы компьютера (память и время работы программы). Но по оси с2 отсчеты следует делать чаще, по крайней мере, на порядок.
Таким образом, описанная математическая модель можно признать пригодной для расчета характеристик алгоритмической составляющей погрешности метода для ее дальнейшего исполь-
/ МК.Ч
' °'90 1 1
0,60 ОМ
■ *
а)
¡а
г- . Н1- ж , /Ш, 1
У Ж 1 '
б -
■1
5 10 15 20 25 30 35 40 45 координата 1/Д1, ош. ед.
В)
Рис. 3. Микротопографические карты поверхности: - деталь из титанового сплава после вибрационного галтования в стальных шариках; стальная деталь после плоского шлифования периферией круга; а) и б) - согласно [4];
в - результат моделирования
а
зования при проверке функции преобразования и оптимизации оптической схемы экспериментального лазерного кругломера.
СПИСОК ЛИТЕРАТУРЫ
1. Лабораторный макет лазерного кругломера / В.Н. Белопухов, М.Ю. Грачев, М.В. Лесников, О.А. Заякин, Ю.О. Понявин // ПИТ-2012: сб. науч. тр. Самара: Изд. СНЦ РАН, 2012. С. 10-13.
2. Топорец, А. С. Оптика шероховатой поверхности.
Л.: Машиностроение. Ленингр. отделение, 1988. 191 с.
3. Беннет, Дж М., Маттсон М. Шероховатость поверхности и рассеяние. Вашингтон: Оптическое общество Америки, 1993. 119 с.
4. Дунин-Барковский, И.В., Карташова А.Н. Измерения и анализ шероховатости, волнистости и некругло-сти поверхности. М.: Машиностроение, 1978. 232 с.
5. Методика измерения отклонений от круглости и шероховатости поверхности деталей подшипников качения. Руководящий документ РД 37.006.106.90 Согласован. М.: НПО ВНИПП, 1990 - 40 листов.
MATHEMATICAL MODEL OF SURFACE MONITORING BY AN EXPERIMENTAL LASER ROUNDNESS METER
© 2015 O.A.Zayakin
Samara State Aerospace University named after Academician S.P. Korolyov (National Research University)
A mathematical model of surface monitoring by an experimental laser roundness measuring instrument is being discussed. The device is designed for quality conformance inspection of some mirror bearing surfaces. The procedure is conducted by means of a coordinate metrology technique based on signal processing of a mirrored probing laser beam. As a result, this model has been proved to be suitable for algoritmic component assessing of the method error.
Key words: triangulation, geometrical variables measurement, mathematical model, surface, roundness measuring insrument, mirrored light, topogram, statictical characteristics.
Oleg Zayakin, Candidate of Technics, Associate Professor at thelnformation Systems and Technology Department. E-mail: [email protected]
