Saratov, Saratov, Russian Federation
Saratov, Saratov, Russian Federation
VAC 05.13.18 Математическое моделирование, численные методы и комплексы программ
UDK 62 Инженерное дело. Техника в целом. Транспорт
GRNTI 27.35 Математические модели естественных наук и технических наук. Уравнения математической физики
GRNTI 55.13 Технология машиностроения
A measurement of surface roughness is commonly carried out with the aid of profilometers, but when measuring the roughness of a curved surface arise difficulties. A technical problem consists in a probe feed along the normal to a nominal surface and filtration followed to exclude a low-frequency constituent of a profile. That is why it is purposeful to use coordinate measuring machines which have a probe with a spherical tip. For the first time there is presented a simulator of a spherical probe flat contact and surface roughness. A contact model developed is based on the regulations of analytical geometry. With the use of a numerical algorithm one finds a contact point and the center of a probe circle. On the basis of the simulator there is carried out an investigation of the impact of a probe radius upon roughness parameters Ra, Rmax, Rq. As initial data were used measurement results of the profilometer and profiles modeled with the aid of Monte-Carlo statistical method. It is defined that with the increase of the radius of the spherical probe its penetrating capacity decreases. That is why height parameters of roughness are distorted. The application of the spherical probe with the radius of 5-50 mkm ensures satisfactory results.
roughness measurement, simulator, coordinate measuring machine, CMM, profilometer, spherical probe
Введение
Традиционно для измерения шероховатости используют профилометры [1; 2]. Щуп профилометра, как правило, представляет собой конус с углом 90° и радиусом при вершине 2, 5 или 10 мкм. Геометрические размеры щупа создают первый источник погрешностей при измерении шероховатости. Относительно мелкие или узкие элементы поверхности будут не полностью измерены и искажены. В этом случае щуп осуществляет механическую высокочастотную фильтрацию микронеровностей. Уменьшение угла конуса или радиуса при его вершине приводит к увеличению контактного давления, что, в свою очередь, может вызвать заклинивание щупа или царапание измеряемой поверхности.
Измерение шероховатости на сложных криволинейных поверхностях, а также труднодоступных участках, например в отверстиях малого диаметра, с помощью профилометров затруднительно [3-5]. Вместе с тем такая задача достаточно часто встречается на практике [6-11]. Тогда целесообразным становится применение координатно-измерительных машин (КИМ). Современные КИМ в большинстве случаев оснащают сферическими щупами диаметром 0,1, 2, 8 мм. С увеличением диаметра щупа увеличивается производительность контроля при прочих равных условиях. Наименьшее из известных значений радиуса сферы наконечника 0,03-0,01 мм имеется у КИМ фирмы «Лапик» - единственного отечественного производителя координатно-измерительных машин. Также КИМ имеет большую дискретность отсчета - 5 мкм, в то время как у профилометра она обычно составляет 0,5 мкм. Получение шага менее 5 мкм возможно, но требует большего времени при контактном способе. Использование сканирующих датчиков контакта позволяет многократно увеличить число измеренных точек, но при этом увеличивается погрешность. Таким образом, специфика измерения шероховатости с применением КИМ обусловливает необходимость оценки влияния радиуса щупа на параметры шероховатости. Для этого необходима математическая модель контакта сферического щупа и шероховатой поверхности.
Постановка задачи
Для аналитического описания измерения на КИМ требуется решение нескольких последовательных и взаимосвязанных задач: компенсация радиуса сферического щупа, выравнивание систем координат КИМ и детали, фильтрация результатов и исключение грубых ошибок. Компенсация радиуса сферического щупа позволяет получить координаты точек касания с поверхностью на основе координат центра щупа. Выравнивание систем координат исключает из рассмотрения номинальный профиль, что важно для криволинейных поверхностей. Фильтрация призвана отсечь шум (вибрации), а также выявить дефекты поверхности.
Одной из важных задач при измерении шероховатости является обоснование выбора метода фильтрации. Традиционно применяют три типа фильтров, имеющих одну передаточную функцию, но разные длины отсечки (ISO 16610-1:2006): 2RC - линейный Гаусса; 2RC PC - линейный Гаусса с коррекцией фазы; фильтр двойной фильтрации Гаусса - фильтр подавления впадин. Для КИМ, в отличие от профилометров, не делают аппаратную фильтрацию, поэтому данный вопрос остается актуальным. Вопросы фильтрации изложены в работах [12; 13].
Таким образом, при измерении шероховатости на КИМ актуальной будет задача определения влияния радиуса щупа на параметры шероховатости. Для этого необходимо разработать математическую модель контакта щупа и измеряемой поверхности. Указанная задача в известной литературе описана только с эвристической точки зрения [1]. Попытка аналитического описания была предложена в работе [14]. В настоящей статье представлено развитие аналитической модели и приведены результаты математического моделирования на основе реальных данных и с использованием метода Монте-Карло.
Математическая модель измерения
Математическая модель плоского контакта сферического щупа и шероховатой поверхности построена на положениях аналитической геометрии. Расчетная схема приведена на рис. 1. В качестве исходных данных выступают: исходный профиль поверхности в виде координат (xi, yi), радиус r сферического щупа, дискретность Dxj перемещения щупа вдоль измеряемого профиля.
Рис. 1. Расчетная схема измерения шероховатости
Исходные допущения модели:
- профиль шероховатой поверхности задан дискретно, так как аналитическое описание отсутствует;
- погрешность формы и расположения поверхности исключена до рассмотрения шероховатости;
- центр сферы щупа рассматривается в плоскости профиля измеряемой поверхности, поэтому его сечение представляет собой окружность с радиусом, равным радиусу сферы;
- направление подвода щупа совпадает с направлением координатной оси Y;
- радиус щупа больше, чем дискретность Dxj перемещения, поэтому контакт осуществляется только в точках xi;
- фильтрация шероховатости выполняется стандартными методами.
Аналитическое описание процесса измерения сводится к следующему. Сферический щуп (окружность в сечении) перемещается по оси из некоторого исходного положения yj, заведомо с зазором до поверхности, до касания с точкой поверхности. Касание находится в результате расчета минимального зазора между окружностью щупа и всеми точками поверхности в пределах диаметра щупа. Для этого решается параметрическое уравнение окружности, приравненное к xi, относительно угла ji. После этого центр окружности 0j смещается в направлении поверхности на величину Di min.
Схема алгоритма расчета представлена на рис. 2. Точка касания окружности щупа и профиля находится численным методом, исходя из минимума начального зазора Di min.
Рис. 2. Схема алгоритма моделирования измерения
шероховатости
Организуется два вложенных цикла: по j от 1 до n и по i от 1 до m, где m - число точек измеряемого профиля, n - число положений щупа при измерении. В общем случае числа n и m не равны, но для последующего моделирования они приняты равными. В первом цикле моделируется перемещение щупа по длине измерения по дискретно заданным координатам xj, yj. Вертикальное положение yj min соответствует минимальному расстоянию xi min между окружностью щупа и дискретно заданным профилем поверхности. Во втором (внутреннем) цикле определяются расстояния от всех точек профиля до окружности в пределах от xj – r до xj + r. Затем из этих расстояний выбирается минимальное значение xi min. Предполагается, что именно эта точка станет точкой контакта при вертикальном перемещении щупа.
Моделирование измерения шероховатости
На основе разработанных математической модели контакта щупа и шероховатой поверхности и программы расчета выполнены тестовые примеры для предварительно измеренной на профилометре поверхности. Для получения статистически достоверных данных проведено моделирование измерения 100 профилей, полученных по методу Монте-Карло по аналогии с работой [10]. Исследовались параметры шероховатости: среднеарифметическое отклонение профиля Ra, наибольшая высота профиля Rmax, среднеквадратичное отклонение профиля Rq. Усредненные данные по результатам моделирования приведены в таблице.
Таблица
Измеренные параметры шероховатости в зависимости от радиуса щупа
|
Параметры шероховатости |
Радиус щупа, мкм |
|||
|
5 |
20 |
50 |
100 |
|
|
Ra, мкм |
0,86 |
0,84 |
0,80 |
0,78 |
|
Rmax, мкм |
1,54 |
1,54 |
1,40 |
1,34 |
|
Rq, мкм2 |
0,69 |
0,68 |
0,67 |
0,63 |
Пример на рис. 3 иллюстрирует моделирование измерения шероховатости поверхности с различными радиусами сферического щупа (5, 20 и 100 мкм). Дискретность отсчета по оси X составляла 5 мкм, дискретность отсчета по оси Y - 0,1 мкм. Из профиля предварительно была исключена постоянная составляющая на основе фильтрации Гаусса.
Рис. 3. Моделирование измерения шероховатости сферическим щупом:
а - исходный профиль; б - радиус 5 мкм; в - радиус 20 мкм; г - радиус 100 мкм
Анализ результатов моделирования показал, что с увеличением радиуса сферического щупа уменьшается его проникающая способность. Поэтому происходит недоощупывание профиля, что уменьшает получаемые значения Ra, Rmax, Rq на 10-20 %. Шаговые параметры профиля, например средний шаг неровностей Sm, искажаются незначительно.
Заключение
Проведенные исследования позволяют сделать вывод, что измерение шероховатости с применением координатно-измерительных машин дает положительный результат при использовании сферических щупов диаметром 0,01-0,03 мм.
1. Uaythauz D. Metrologiya poverhnostey. Principy, promyshlennye metody i pribory. M.: Intellekt, 2009. 472 s.
2. Tabenkin A.N., Tarasov S.B., Stepanov S.N. Sherohovatost', volnistost', profil'. Mezhdunarodnyy opyt. SPb.: Izd-vo Politehn. un-ta, 2007. 136 s.
3. Fomin A.A. Obespechenie mikrogeometrii poverhnostey pri obrabotke zagotovok s neodnorodnymi svoystvami // Sborka v mashinostroenii, priborostroenii. 2012. № 12. S. 27-29.
4. Shrubchenko I.V. Opredelenie harakteristik pyatna kontakta poverhnostey kacheniya opor tehnologicheskih barabanov // Sborka v mashinostroenii, priborostroenii. 2007. № 2. S. 4-6.
5. Pechenin V.A., Bolotov M.A., Ruzanov N.V. [i dr.]. Model' sopryazheniya detaley po plosko-cilindricheskim poverhnostyam // STIN. 2017. № 3. S. 22-28.
6. Dem'yanenko E.G. Issledovanie sposoba formoobrazovaniya detaley bochkoobraznoy formy // Aviacionnaya tehnika. 2014. № 2. S. 72-76.
7. Zaharov O.V., Pogorazdov V.V. Proektirovanie formoobrazuyuschih sistem bescentrovyh superfinishnyh stankov. Saratov: SGTU, 2004. 140 s.
8. Zaharov O.V. Bescentrovoe shlifovanie konicheskih poverhnostey na stankah s prodol'noy podachey // Avtomatizaciya i sovremennye tehnologii. 2006. № 7. S. 14-16.
9. Zaharov O.V. Naladka bescentrovyh superfinishnyh stankov na osnove chislennogo modelirovaniya i optimizacii // Vestnik mashinostroeniya. 2003. № 12. S. 48-50.
10. Zaharov O.V. Upravlenie tochnost'yu bescentrovogo shlifovaniya statisticheskimi metodami // Mehatronika, avtomatizaciya, upravlenie. 2009. № 9. S. 32-35.
11. Korolev A.V., Balaev A.F., Yakovishin A.S. Obespechenie kachestva podshipnikov kacheniya na osnove primeneniya tehnologii mnogociklovogo nagruzheniya // Sborka v mashinostroenii, priborostroenii. 2017. № 7. S. 291-293.
12. Pechenin V.A., Bolotov M.A., Stepanova E.R. Ispol'zovanie bilateral'nogo fil'tra dlya ocenki i analiza otkloneniy geometrii poverhnostey // Informacionnye tehnologii i nanotehnologii: sb. tr. konf. Samara, 2016. S. 447-452.
13. Markov B.N., Melikova O.N., Shulepov A.V. Algoritm postroeniya morfologicheskogo diskovogo fil'tra dlya analiza sherohovatosti poverhnosti // Izmeritel'naya tehnika. 2017. № 5. S. 30-33.
14. Koroleva A.A., Zaharov O.V. Ocenka pronikayuschey sposobnosti sfericheskogo schupa pri izmerenii sherohovatosti // Innovacionnye tehnologii v metalloobrabotke: sb. tr. vseros. nauch.-prakt. konf. Ul'yanovsk: UGTU, 2019. S. 351-355.
