Математическое моделирование профильной фрезы для обработки деталей цилиндрической передачи | Вестник Томского государственного университета. Математика и механика. 2016. № 4(42).

Математическое моделирование профильной фрезы для обработки деталей цилиндрической передачи

При обработке поверхностей деталей передаточного механизма используются различные типы фрез (сферические, тороидальные и др.). Стоимость специального профильного инструмента несколько выше, чем таких фрез. Но увеличение стоимости инструмента окупается значительным сокращением машинного времени обработки деталей. Построена математическая модель поверхности профильной фрезы (как поверхности вращения) для обработки деталей цилиндрической передачи с эксцентриково-циклоидальным зацеплением. Создана компьютерная программа для определения радиусов круговых сечений фрезы для заданного набора смещений по оси вращения.

Mathematical simulation of a profile cutter for processing parts of a cylindrical gear.pdf В [1, 2] были построены математические модели процесса формообразования деталей цилиндрического и конического передаточного механизмов с эксцентриково-циклоидальным зацеплением (ЭЦ-зацеплением [3]) с помощью сферических и тороидальных фрез. В [4] было рассмотрено построение поверхности профильной фрезы для обработки входной детали цилиндрического механизма с ЭЦ-зацеплением В данной работе моделируется построение поверхности профильной фрезы для обработки обеих деталей (входной и выходной) цилиндрического механизма с ЭЦ-зацеплением. Профильные фрезы позволяют получать поверхность впадины зубчатых колёс методом копирования «за один заход», в отличие от сферических и тороидальных фрез, обработка которыми происходит последовательным (многозаходным) удалением части металла с заготовки детали. Увеличение стоимости профильного инструмента окупается значительным сокращением машинного времени обработки деталей. Входную и выходную детали передаточного механизма принято называть «шестерней» и «колесом» соответственно. 1. Моделирование поверхности профильной фрезы для обработки зуба шестерни Выпишем константы, входящие в уравнения поверхностей деталей: Aw - межцентровое расстояние (расстояние между параллельными осями вращения деталей); е - эксцентриситет; z1 - количество зубьев шестерни; z2 - количество зубьев колеса; n = z2/zi + 1; р - радиус окружности сечения зуба шестерни; l - размер деталей по оси вращения; rc - радиус цилиндрического барабана, на котором расположены зубья шестерни; K - количество сечений поверхности фрезы; ц,- - высота подъема сечения (i = 0, 1,..., K); r,- - радиусы окружностей сечения фрезы, для различных значений ц,-; П - сдвиг оси фрезы. Поверхность зуба шестерни в цилиндрической передаче с ЭЦ-зацеплением [1] образована окружностями, расположенными в параллельных плоскостях, причем центры этих окружностей лежат на винтовой линии, принадлежащей цилиндру радиуса е (червячный элемент). Примем ось вращения колеса за ось OZ, а ось вращения шестерни направим по прямой, параллельной оси OZ, смещенной вдоль оси OX на величину Aw. Тогда параметрические уравнения поверхности зуба шестерни можно записать в виде вектор-функции двух аргументов: ' Aw -e cos и + р cos а ^ -е sin и + р sin а lr и z, (1) S (и, а) = 2п где а = 0,..,2п, и = 0,..,2n/z1. Поверхность фрезы будем строить как семейство окружностей в плоскостях, перпендикулярных оси фрезы, с центрами, лежащими на этой оси. При этом радиус каждой из этих окружностей будет определяться из условия касания ее с плоской кривой - сечением зуба плоскостью окружности. Идеальная прямая оси фрезы пересекает ось вращения шестерни под прямым углом и расположена на равных удалениях от зубьев шестерни. На практике, дабы при обработке одной боковой поверхности зуба избежать контакта с другой боковой поверхностью, вводят величину п малого параллельного сдвига оси фрезы (сохраняя перпендикулярность оси вращения шестерни). На рис. 1 изображено положение оси фрезы в плоском сечении шестерни, перпендикулярном оси ее вращения, для z1 = 4. Четыре окружности - сечения зубьев шестерни имеют радиус р, большая окружность - радиус rc. При таком расположении оси фрезы ее уравнения в плоскости сечения имеют вид + п sin | - I п 1 - cos | - I V z J + ц f \ • I п 1 sin | - I V V z1 J J Aw - rc cos | - z Osf (ц ) = (2) П + п cos rc sin При ц = 0 из (2) получаем координаты точки M0, через которую проходит плоскость первой окружности сечения поверхности фрезы. Положение плоскости последней такой окружности определяет точка пересечения M1 оси фрезы с касательной к окружности сечения зуба шестерни, перпендикулярной этой оси (см. рис. 1). Запишем уравнение этой касательной в нормальном виде: - sin | If y - р sin I п II = 0. (3) п cos | - ■ z1 х-| Aw - р cos | - I-е Рис. 1. Сечение входной детали (шестерни) цилиндрической передачи в ЭЦ-зацеплении для z1 = 4. Показано положение оси профильной фрезы и касательной к окружности сечения первого зуба шестерни, перпендикулярной оси фрезы (пунктирная линия) Fig. 1. Cross section of the cylindrical transmission inlet part (gear) in EC-gearing for z1 = 4. Position of the profile cutter axis and tangent to the section circle of the first gear tooth perpendicular to the cutter axis (dashed line) Тогда расстояние д1 от точки M0 до касательной (3) получится в результате подстановки в (3) координат точки M0: + Р - rc • д1 = е cos Таким образом, плоскости, перпендикулярные оси фрезы, в которых располагаются окружности, образующие поверхность фрезы, проходят через точки M0 и M1 оси фрезы, т.е. при изменении параметра д в (2) от 0 до дь Уравнения семейства этиих плоскостей можно записать в виде - sin| ZL Jl y - (rc + ^)sinl -П || = 0, (4) cos I f х-I Aw -(rc + ^)cosl - где д = 0,...,дь Теперь получим уравнения семейства кривых - сечений поверхности зуба шестерни плоскостями семейства (4). Для этого подставим в (4) вместо х и y две первые координаты вектор-функции (1), задающей поверхность зуба шестерни. В результате получим соотношение между параметрами и, а, и д, из которого можно получить выражение и через а и д: п+а z ( Л + rc + д р cos и(а, д) = -- +arccos (5) Z1 Sem (а, д) = (6) ' Aw - e cos^fe д)) + р cos а А -е sin^fe д)) + р sin а lr и(а, д) z1 2п Именно этих кривых должны касаться окружности, образующие поверхность профильной фрезы. Эти окружности располагаются в плоскостях семейства (4), их центры лежат на оси ферзы (2), а радиусы r меняются в зависимости от Подставляя (5) в (1), получаем искомые уравнения семейства плоских кривых - сечений зуба шестерни плоскостями семейства (4): параметра д семейства (4). Исходя из этого, запишем уравнения семейства таких окружностей (пока с произвольным радиусом г) в виде Aw -(гс + ^)cosI ~~1 + (гcosв + n)sinl - '1J (гс + д^ш| - 1 + (гcose + n)cos П Okr (в, д) = (7) г sin в В формулах (5) - (7) д = 0,...,д1. В каждом сечении д окружность (7) должна касаться кривой (6), т.е должно обращаться в нуль скалярное произведение векторов: (М0 - Sem(а, д), Sem (а, д)) = 0, (8) здесь М0 = Osf (0), а через Sem (а, д) обозначена производная вектор-функции Sem (а, д) по параметру а: С s sin(u(a, д)) и'(а, д) -р sin а ^ e cos(u(a, д)) и'(а, д) -р cos а 1г и'(а, д) z1 Sem (а, д) = v 2п . С п+а Zi р sin | -1 и'(а, д) = v 1 Для решения уравнения (8) составлена компьютерная программа, которая для каждого значения д1 д.- = - 1 1 K (K - задаваемое количество сечений поверхности фрезы) определяет начальное приближение для нахождения корня функции в левой части уравнения (8), а затем с помощью встроенной функции root пакета MathCad находит этот корень а. Далее, для определения радиуса окружности г, касающейся сечения зуба шестерни, соответствующего д., потребуем наличия общей точки у этих кривых, т.е. приравняем координаты вектор-функций (6) и (7). Равенство вторых и третьих координат дает систему двух уравнений: -e sin(u(a, д)) + р sin a - (ге + д^ш| - г cosв = п cos | - Z1 г sin в = 1г и(а, д) z1 2п из которой находим искомые радиусы окружностей для каждого значения д,: -е sin^^, д,)) + р sin а, - (rc + д, )sin| - cos | - Z1 lr и(а,-, д,) Z1 2п r = Таким образом, поверхность профильной фрезы получена как семейство окружностей вида (7), т.е. как поверхность вращения: Aw - (rc + д, )cosl ~~J + (r cosp + n)sinl - 4 J (rc + д,)sin| - 1 + (r cosp + n)cos П r sin p На рис. 2 изображена профильная фреза, смоделированная по вышеописанной схеме. Показана плоскость сечения фрезы, перпендикулярная ее оси. В этой плоскости выделена линия сечения зуба шестерни; этой линии касается окружность сечения фрезы. Fr (P, i) = Рис. 2. Профильная фреза для обработки поверхности зуба шестерни в передаче с ЭЦ-зацеплением. Показана плоскость сечения, перпендикулярная оси фрезы, и линия сечения зуба шестерни этой плоскостью Fig. 2. Profile cutter for processing the surface of the gear tooth in a transmission with EC-gearing. The cross section plane (perpendicular to the cutter axis) and line of the tooth section by this plane are shown 2. Моделирование поверхности профильной фрезы для обработки зуба колеса Поверхность выходной детали (колеса) образована вращением плоской циклоидальной кривой вокруг оси вращения колеса с одновременным смещением ее вдоль этой оси. Этой кривой является эквидистанта эпитрохоиды [5]. Зададим эпитрохоиду в виде вектор-функции аргумента т = 0,...,2ля: -e cos т + Aw cos- n -e sin т + Aw sin- Те(т) = v n У Эквидистанту эпитрохоиды (кривую, удаленную по нормалям от эпитрохоиды на величину р) запишем в виде Ne(T) (9) E (т) = Те(т) + р |Ne(T) где Aw f Л e cos т--cos - n n Aw т e sin т--sin- v n n - вектор нормали в точке эпитрохоиды. Поверхность выходной детали запишем в виде вектор-функции двух аргументов: Ne(T) = и cos-Е(т)1 + sin n -1 и - sin-Е(т)1 + cosn -1 n -1 Ir - и 2п ( U Е(т)2 ^ n -1 Ev(u, т) = Е(т)2 (10) где и = 0,...,2п - угол поворота входной детали, u/(n-1) - угол поворота эквиди-станты. Здесь обозначения E(t)j , E(т)2 означают соответствующие координаты вектор-функции (9), которые являются скалярными функциями. Ось OX лежит на оси симметрии первого профиля зуба колеса и проходит через точку Ev(0,0). Начальное положение оси фрезы - параллельно оси OX со смещением по оси OY на малую величину n. Это смещение, как и в случае шестерни (см. п.1), необходимо для того, чтобы при обработке одной боковой поверхности зуба не касаться другой боковой поверхности. Уравнение оси фрезы можно записать в виде вектор-функции Osf (ц ) = Г E (0)1 + ц n 0 Теперь запишем уравнение семейства плоскостей, перпендикулярных оси фрезы, в виде X = E (0)1 (11) и найдем кривые, получающиеся в сечениях этими плоскостями поверхности (10). Для этого приравняем правую часть (11) первой координате вектор-функции (10) и из получившегося уравнения и 1 £(х)1 + sin | - IE (т)2 = E(0)1 + ц 4n -1J Vn -1 найдем зависимость между параметрами и и т: cos E(x)t [E(т)2(E(0)t + ц) + Ed(т,ц)] E(x)j1( E (0)! + ц) - E (т)2 Ed (т,ц). ( и(т, ц) = (n -1) • arctg Ed(т, ц) = yj-E(x)!2((E(0)! + ц)2 - E(x)!2 - E(t)22) . Подставляя и(т, ц) в (10) вместо и, получаем уравнение семейства сечений зуба колеса плоскостями, перпендикулярными оси фрезы (д - параметр семейства): С E(0)1 + д 1 где , и(т, д) 1 С и(т, д) - sin |-- IE(т) + cos I-| E(т)2 Sem (т, д) = (13) n -1 n -1 и(т, д)-1г 2п Этих кривых должны касаться окружности, образующие поверхность профильной фрезы. Эти окружности располагаются в плоскостях семейства (13), их центры лежат на оси фрезы, а радиусы г меняются в зависимости от параметра д семейства (13). Исходя из этого, запишем уравнения семейства таких окружностей в виде С E(0)! +ц 1 г cos(a) +п г sin(a) Okr (а, д) = (14) Требование касания окружностей (14) с плоскими кривыми сечений (13) осуществляется как и в п.1 в два этапа. Во-первых, приравниваем координаты вектор-функций (13) и (14), что приводит к системе двух уравнений (первые координаты совпадают тождественно), из которых можно найти зависимости г и а от т для каждого значения д: Это уравнение на т при переходе к координатам имеет весьма громоздкий вид и решается с помощью специальной компьютерной программы. В результате для каждого значения высоты подъема сечения ц получаем значение тц = т(ц). Подставляя в (15) т(ц) получаем значения радиусов r(ц) окружностей сечений фрезы для каждого значения ц. Теперь из (14) получается уравнение поверхности фрезы: ' E (0)1 + | ^ Fr (а, ц) = r (|)cos(a) + n r (|)sin(a) У Если т(ц) подставить в (16), то получим параметр а(ц) точки касания окружности сечения фрезы с кривой сечения зуба колеса (для каждого значения ц). Координаты всех точек касания получатся при подстановке а(ц) в (17). Эти точки показаны на рис. 3. Рис. 3. Профильная фреза для обработки зуба колеса Fig. 3 Profile cutter for processing the wheel tooth Полученные координаты точек поверхности профильного инструмента используются для изготовления такого инструмента на специальных станках. Это с успехом используется в ЗАО «Технология маркет» (Томск, ОЭЗ) при необходимости резкого сокращения времени обработки зубчатых колёс для мелкосерийного производства.

Ключевые слова

профильная фреза, эксцентриково-циклоидальное зацепление, касание кривых, profile cutter, eccentrically cycloidal (EC) gearing, contingence of curves

Авторы

ФИООрганизацияДополнительноE-mail
Бубенчиков Алексей МихайловичТомский государственный университетдоктор физико-математических наук, профессор, заведующий кафедройbubenchikov@mail.tomsknet.ru; alexy121@mail.ru
Казакявичюс Сергей МатвеевичЗАО «Технология маркет»инженерtm@ec-gearing.ru
Костарев Артем АндреевичТомский государственный университетстудент механико-математического факультетаartem_kost@bk.ru
Щербаков Николай РомановичТомский государственный университетдоктор физико-математических наук, доцент, заведующий кафедрой геометрииnrs@math.tsu.ru
Щербакова Ирина ВладимировнаТомский государственный университеткандидат физико-математических наук, доцент кафедры общей математикиnrs@math.tsu.ru
Всего: 5

Ссылки

Камчатный С.А., Сковородин А.В., Становской А.В., Щербаков Н.Р. Математическое моделирование процесса сближения поверхностей при формообразовании деталей передаточного механизма с ЭЦ-зацеплением // Вестник ТГУ. Математика и механика. 2012. № 4(20). С. 15-24.
Камчатный С.А., Кузнецов В.М, Сковородин А.В., Щербаков Н.Р. Математическое моделирование процесса формообразования деталей конического передаточного механизма с ЭЦ-зацеплением // Вестник ТГУ. Математика и механика. 2014. 2(28). С. 5-17.
Kazakyavichyus S.M., Stanovskoy V.V., Remneva T.A., Kuznetsov W.M., Bubentchikov A.M., Shcherbakov N.R. Perfomance of eccentric-cycloid engagement with change in the interaxial distance: Modification of tooth configuration // Russian Engineering Research. 2011. V. 31. No. 3. P. 197-199.
Bubenchikov A.M., Kazakavitschyus S.M., Shcherbakov N.R. Mathematical simulation of a profile cutter as a surface of revolution // IOP Conf. Series. Materials Science and Engineering. 124 (2016) 012095. doi:10.1088/1757-899X/124/1/012095.
Савелов А.А. Плоские кривые. М.: ГИФМЛ, 1960. С. 118.
 Математическое моделирование профильной фрезы для обработки деталей цилиндрической передачи | Вестник Томского государственного университета. Математика и механика. 2016. № 4(42).

Математическое моделирование профильной фрезы для обработки деталей цилиндрической передачи | Вестник Томского государственного университета. Математика и механика. 2016. № 4(42).