CC BY
168
ТЕОРИЯ МАШИН И МЕХАНИЗМОВ
УДК 621.833.15 DOI: 10.17586/0021-3454-2016-59-8-619-624
ВЕРОЯТНОСТНЫЙ МЕТОД РАСЧЕТА ПАРАМЕТРОВ ТОЧНОСТИ ЭВОЛЬВЕНТНЫХ ЦИЛИНДРИЧЕСКИХ ЗУБЧАТЫХ ПЕРЕДАЧ
М. В. Абрамчук
Университет ИТМО, 197101, Санкт-Петербург, Россия E-mail: [email protected]
Рассматривается вероятностный метод расчета параметров точности эвольвент-ной цилиндрической зубчатой передачи, а именно — монтажного радиального биения зубчатой передачи и кинематической погрешности зубчатой передачи. В качестве вероятностного метода расчета применяется метод Монте-Карло.
Ключевые слова: зубчатые колеса, зубчатые передачи, точность, многозвенный механизм, метод Монте-Карло, кинематическая погрешность
Введение. Точность — один из основных показателей качества для зубчатых передач. Система нормирования параметров точности зубчатых колес и передач, существующая в России, давно нуждается в пересмотре. Прогресс в области производства передач, редукторов и мультипликаторов, отмечаемый в течение последних десятилетий, не отражен в отечественных нормативных документах, которые зачастую резко расходятся с рекомендациями ISO [1—3]. Следствием этого является неконкурентоспособность производимых зубчатых передач и многозвенных зубчатых механизмов.
На точность зубчатой передачи влияют не только погрешности зубчатых колес, но и погрешности других элементов передачи, а именно: валов, подшипников, корпуса и т.д. Вследствие влияния погрешностей монтажа степень точности зубчатого колеса отличается от значения, указанного в стандарте. Функциональные показатели передач не контролируются на предприятиях (контроль предусмотрен лишь для отсчетных передач), что обусловливает снижение качества зубчатых передач.
Следовательно, главным недостатком системы нормирования параметров точности передачи является отсутствие методик, учитывающих погрешности всех элементов зубчатой передачи. Поэтому для внесения предложений по пересмотру стандартов, в частности ГОСТ 1643-81 и ГОСТ 21098-82, необходимо разработать методы перехода от геометрических показателей точности отдельных деталей передачи к функциональным показателям зубчатых передач и многозвенных зубчатых механизмов.
Для расчета зубчатой передачи на точность в настоящей статье предлагается следующий метод: определение точностных параметров передачи по точности зубчатых колес и других элементов передачи, а также точности монтажа.
Расчет монтажного радиального биения зубчатой передачи. При расчетах на максимум-минимум учитываются только предельные отклонения параметров. В результате этого расчета определяется максимальный допуск на монтажное радиальное биение.
В качестве вероятностного метода предлагается применить метод Монте-Карло [4, 5]. Алгоритм расчета заключается в следующем: для каждого параметра случайным образом
выбирается его значение (будем считать все составляющие монтажного радиального биения распределенными, например, равномерно), затем по выведенной формуле для монтажного радиального биения вычисляется его значение. Опыт повторяется N раз, после чего определяются математическое ожидание и дисперсия или среднеквадратическое отклонение. Далее строятся гистограммы и используются общепринятые критерии для оценки соответствия результата некоторому закону распределения.
Максимальная кинематическая погрешность для цилиндрических зубчатых передач по ГОСТ 21098-82 определяется в соответствии со следующим выражением:
КГ
Ютах
= К
\2 , Е2 + ЕЕМ
2 + Е 2
(1)
где К — коэффициент фазовой компенсации; К'г1, ¥\2 — допуски на кинематические погрешности соответственно шестерни и колеса (определяются по ГОСТ 1643-81 или ГОСТ 9178-81); , Еш 2 — суммарная приведенная погрешность монтажа (по ГОСТ 21098-82).*
Допуск на кинематическую погрешность зубчатого колеса
п= Кр + //,
где Кр — допуск на накопленную погрешность шага зубчатого колеса, / грешность профиля зуба.
Суммарная приведенная погрешность монтажа рассчитывается по формуле
ег
допуск на по-
Етм ='
СОБР
+еа ,
где ег и еа — монтажные радиальное и осевое биения зубчатого колеса соответственно (показатели точности монтажа зубчатых колес); а — угол исходного профиля; в — делительный угол наклона линии зуба.
В общем случае возможно следующее расположение зубчатого колеса на валу [6, 7] относительно опор: межопорное (рис. 1, а) и консольное (рис. 1, б).
а) б)
Рис. 1
Информация о размерах а, Ь, I необходима для приведения радиального биения колец подшипников к средней плоскости венца зубчатого колеса.
Для расчета монтажного радиального биения определим допуски на погрешности, создающие первичные радиальные биения:
— ¥г — допуск на радиальное биение зубчатого венца колеса, определяется по ГОСТ 1643-81 и ГОСТ 9178-81;
— еп — зазор в посадке колеса на вал, зависящий от конструкции соединения и используемой посадки;
Здесь и далее использованы обозначения, принятые в соответствующих стандартах.
Вероятностный метод расчета точности зубчатых передач 621
— ев — допуск на радиальное биение посадочной ступени вала под зубчатое колесо относительно цапфы вала, причем для гладких валов ев = 0, определяется по ГОСТ 24643-81;
— Щ — радиальное биение внутреннего кольца собранного подшипника, определяется по ГОСТ 520-2011.
Необходимо отметить, что при наличии типовых схем зубчатых передач можно привести более полный перечень погрешностей, создающих первичные радиальные биения. В данной работе для расчета монтажного радиального биения ограничимся указанными выше погрешностями.
С учетом сказанного монтажное радиальное биение зубчатого колеса определяется следующим образом:
ег ^ + еп + ев + Ъ^у-^.
Для прямозубых зубчатых колес осевыми биениями можно пренебречь, т.е. еа =0, тогда
Е1Ы = ег ^ а .
Таким образом, алгоритм расчета монтажного радиального биения вероятностным методом будет следующим.
1. Выводится выражение для расчета и предполагается, что слагаемые в данном выражении являются случайными величинами.
2. Выбираются нижние и верхние предельные значения для каждого слагаемого:
— по таблицам соответствующих стандартов (верхний предел — заданная степень точности, нижний предел — на одну-две степени выше);
— на основе измерений параметров партии передач при производстве.
3. Задается распределение для каждого слагаемого.
4. При каждом конкретном вычислении получается случайная величина, а результатом расчета является гистограмма.
5. Осуществляется аппроксимация результата, используется критерий соответствия какому-либо закону распределения.
6. Определяются характеристики распределения (математическое ожидание, дисперсия и т.д.).
Расчет кинематической погрешности зубчатой передачи вероятностным методом.
Кинематическая погрешность — это функция, которую можно разложить на гармонические составляющие: оборотной частоты, частоты входа зубьев в зацепление и высокочастотные. Каждая гармоническая составляющая — это случайная функция с фиксированной частотой, случайной амплитудой и начальной фазой, поэтому результат измерения кинематической погрешности не будет представлять собой гармонику: рис. 2, где — наибольшая кинематическая погрешность, измеряемая в течение периода т; ф2 — угол поворота колеса в зубчатой передаче.
К
Рис. 2
Поскольку кинематическая погрешность — величина периодическая, ее период определяется по формуле
т = (2)
х
где 21 — число зубьев шестерни в передаче; х — наибольший общий делитель для числа зубьев шестерни и колеса.
Если передаточное число является целым, то х = 2\, а т = 2п. Если передача не имеет общего делителя, то т = 2л2ь Следовательно, 2 л: <т< 2л2 . Например, при г1=36 число зубьев колеса г2=42, общий делитель х=6, период т = 12п, т.е. чтобы оценить кинематическую погрешность передачи, необходимо зубчатое колесо повернуть на 6 оборотов, а шестерня сделает при этом 7.
Кинематическая погрешность зубчатой передачи является результатом влияния гармонических составляющих с различной частотой повторений к за период т (к>1 — целое число) [8].
Удвоенная амплитуда гармонической составляющей кинематической погрешности называется циклической погрешностью и обозначается/2ког, а ее допуск — /2ко (рис. 3). Удвоенная амплитуда гармонической составляющей с частотой повторений к2, равной частоте входа зубьев в зацепление, называется циклической погрешностью зубцовой частоты и обозначает-ся/-гог, а ее допуск —/=а.
Частота повторений к1 соответствует числу оборотов шестерни, к2 — числу оборотов колеса, к2 — числу входов зубьев в зацепление. Значения данных параметров для различных пар зубчатых колес приведены в таблице.
21 22 х т к1 к2 к2
48 96 48 2п 2 1 96
60 80 20 6п 4 3 240
25 30 5 10п 6 5 150
Частоты повторений определяются как
к = к2 = х х
к2=к111=к222.
(3)
Тогда, с учетом изложенного, выражение (1) для определения кинематической погрешности примет следующий вид:
+
Ъ 0 =
+
рР1 в^п (к1Ф1 +ФР1 н ) +//1 + Еш1
Рр2 вШ (к 2 Ф1 +ФР2 н ) + //2 + Е1М2 ] + /220 (к2Ф1 + Ф22о н ),
(4)
где Фр1 н, Фр 2 н, Ф220 н — начальные фазы; ф1 — угол поворота ведущего колеса (шестерни).
Вероятностный метод расчета точности зубчатых передач 623
Высокая частота повторений для погрешности профиля зуба ffr [9, 10] обусловливает ее суммирование с другими слагаемыми в формуле (4) без учета гармоник и фаз [11].
Значение кинематической погрешности, в угловых минутах, согласно ГОСТ 21098-82,
360.60^ = F0, (5)
n-1000mz d
где d — диаметр делительной окружности ведомого колеса, мм; m — модуль зацепления, мм.
Алгоритм расчета наибольшей кинематической погрешности вероятностным методом будет следующим.
1. Выводится выражение для расчета — выражение (4).
2. Принимается, что все составляющие кинематической погрешности распределены равномерно.
3. По формулам (3) рассчитываются частоты повторений k1, k2 и kz.
4. По формуле (2) рассчитывается период т.
5. Определяются случайные величины в выражении (4): амплитуды FP1, FP2 и fzzo, допуски fi, ff2 и начальные фазы.
6. Выбираются нижние и верхние предельные значения для каждого слагаемого (для начальных фаз 0.. ,2л):
— по таблицам соответствующих стандартов (верхний предел — заданная степень точности, нижний предел — на одну-две степени выше);
— на основе измерений параметров партии передач при производстве.
7. Задается распределение для каждого слагаемого.
8. Значения угла ф1 подставляются в формулу (4) с выбранным шагом, например 15°, начиная от 0 до конца периода т.
9. Определяется значение наибольшей кинематической погрешности как разность между максимальным F{or и минимальным F0 значениями, определенными по выражению (4).
10. Определяется значение F{or в угловых единицах по выражению (5).
Заключение. Использование рассмотренного вероятностного метода расчета параметров точности зубчатой передачи позволяет получить меньшие значения допусков, что более полно отражает специфику реального производства цилиндрической зубчатой передачи.
При расчете параметров точности зубчатой передачи вероятностным методом необходимо знать границы зоны рассеяния погрешностей и законы их распределения внутри зоны. Установление границ и законов распределения (либо числовых характеристик, когда законы не определены) необходимо осуществлять в налаженном производстве при стабильных числовых характеристиках распределения и достаточно большом количестве измеряемых деталей в выборке.
СПИСОК ЛИТЕРАТУРЫ
1. ISO 1328-1:1995. Cylindrical Gears — ISO System of Accuracy, Part 1: Definitions and Allowable Values of Deviations Relevant to Corresponding Flanks of Gear Teeth.
2. ISO 1328-2:1997. Cylindrical Gears — ISO System of Accuracy, Part 2: Definitions and Allowable Values of Deviations Relevant to Radial Composite Deviations and Runout Information.
3. ISO 1328-1:2013. Cylindrical Gears — ISO System of Flank Tolerance Classification, Part 1: Definitions and Allowable Values of Deviations Relevant to Flanks of Gear Teeth.
4. Соболь И. М. Метод Монте-Карло. М.: Наука, 1978. 64 с.
5. Голенко Д. И. Моделирование и статистический анализ псевдослучайных чисел на электронных вычислительных машинах М.: Наука, 1965. 228 с.
6. Тимофеев Б. П., Абрамчук М. В. Расчет погрешностей для нового базового стандарта, нормирующего точность цилиндрических зубчатых колес и передач // Сб. тез. докл. конгресса молодых ученых: Вып. 2. Проблемы механики и точности в приборостроении. СПб: НИУ ИТМО, 2013. С. 275—276.
7. Тимофеев Б. П., Абрамчук М. В., Смирнов Д. С. Расчет радиального биения колес в зубчатой передаче // Сб. докл. 11-й сессии Междунар. науч. школы, посвященной памяти В. П. Булатова: „Фундаментальные и прикладные проблемы надежности и диагностики машин и механизмов". СПб: Art-Xpress, 2013. С. 331—340.
8. Куцоконь В. А. Точность кинематических цепей приборов Л.: Машиностроение, 1980. 221 с.
9. Тимофеев Б. П., Новиков Д. В. Новые стандарты по точности зубчатых колес и передач — необходимое условие повышения конкурентоспособности отечественного редукторостроения // Вестн. машиностроения. 2013. № 5. С. 85—87.
10. Тимофеев Б. П., Новиков Д. В. Комбинирование норм точности стандартных зубчатых передач // Стандарты и качество. 2013. № 5. С. 48—52.
11. АбрамчукМ. В Совершенствование расчетов параметров точности зубчатых колес, зубчатых передач и многозвенных зубчатых механизмов: Автореф. дис. ... канд. техн. наук. СПб, 2014.
Сведения об авторе
Михаил Владимирович Абрамчук — канд. техн. наук; Университет ИТМО; кафедра мехатроники;
E-mail: [email protected]
Рекомендована кафедрой Поступила в редакцию
мехатроники Университета ИТМО 05.04.16 г.
Ссылка для цитирования: Абрамчук М. В. Вероятностный метод расчета параметров точности эвольвентных цилиндрических зубчатых передач // Изв. вузов. Приборостроение. 2016. Т. 59, № 8. С. 619—624.
PROBABILISTIC METHOD FOR CALCULATION OF PRECISION PARAMETERS OF INVOLUTE CYLINDRICAL GEAR
M. V. Abramchuk
ITMO University, 197101, St. Petersburg, Russia E-mail: [email protected]
Application of a probabilistic method for calculating the precision parameters of involute cylindrical gear is considered. The parameters in question are mounting radial runout and kinematic error of the gear, Monte Carlo method is taken for the probabilistic calculation method.
Keywords: gear wheel, gear, precision, multi-link mechanism, Monte Carlo method, kinematic
error
Data on author
Mikhail V. Abramchuk — PhD; ITMO University, Department of Mechatronics; senior lecturer;
E-mail: [email protected]
For citation: Abramchuk M. V. Probabilistic method for calculation of precision parameters of involute cylindrical gear // Izv. vuzov. Priborostroenie. 2016. Vol. 59, N 8. P. 619—624 (in Russian).
DOI: 10.17586/0021-3454-2016-59-8-619-624
