2014
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА Математика и механика
№ 2(28)
МАТЕМАТИКА
УДК 514.8, 62.342
С.А. Камчатный, В.М. Кузнецов, А.В. Сковородин, Н.Р. Щербаков
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССА ФОРМООБРАЗОВАНИЯ ДЕТАЛЕЙ КОНИЧЕСКОГО ПЕРЕДАТОЧНОГО МЕХАНИЗМА С ЭЦ-ЗАЦЕПЛЕНИЕМ
Формообразование поверхности детали передаточного механизма в самой общей постановке является задачей сближения двух поверхностей (обрабатываемой детали и поверхности инструмента) в пространстве до состояния касания, в котором поверхности имеют общую касательную плоскость. Решение задачи аналитического описания геометрии касания поверхностей необходимо для нахождения наивыгоднейших методов и средств обработки сложных поверхностей деталей на многокоординатных станках с ЧПУ.
С этих позиций в работе построена математическая модель процесса формообразования деталей конического передаточного механизма с эксцентриково-циклоидальным зацеплением с помощью фрез двух видов. Получены системы уравнений для определения местоположения фрезы в момент касания. Найдены аналитические или численные решения этих систем. Разработан комплекс специальных компьютерных программ управления движением фрезы, обрабатывающей деталь.
Ключевые слова: формообразование, эксцентриково-циклоидальное зацепление, касание поверхностей.
В [1] была рассмотрена геометрия касания поверхностей детали и инструмента при формообразовании деталей цилиндрического передаточного механизма с эксцентриково-циклоидальным зацеплением (ЭЦ-зацеплением [2]), т.е. устройства, у которого оси вращения входной и выходной деталей параллельны. В данной работе аналогичные задачи решены для конической передачи с ЭЦ-зацеплением, у которой оси вращения деталей пересекаются (рассмотрен случай пересечения под прямым углом - см. рис.1).
Рис. 1. Входная и выходная детали конической передачи в ЭЦ-зацеплении
1. Параметрические уравнения поверхностей деталей Параметрические уравнения поверхностей деталей
X = Л(и, V),
У = /2(u, ^
г = /э(и, V)
будут записаны в виде вектор-функций двух аргументов, т.е. в виде
(/ (и,v)Л Г (u,v)= /2 (u,v)
/з(u,VI
1.1. Выходная деталь
Пусть дана сфера радиуса Я, которую можно задать вектор-функцией двух аргументов u и V, принимающих значения от 0 до 2п:
£ (и, у) =
008 и 008 у
( Я
Я 008 и БІП у Я єіи и
Л
Обозначим через х1 количество зубьев входной детали (далее - шестерни), через г2 - количество зубьев выходной детали (далее - колеса), а через п - передаточное отношение
22 П =—.
Циклоидальную кривую на сфере описывает точка окружности радиуса є, лежащей на сфере, при одновременных поворотах этой окружности вокруг своего
—а
центра на угол а и вокруг оси вращения колеса (ось 011) на угол — . Такую ок-
п
ружность можно задать поворотом радиус-вектора точки
ЧЯ2 — £2 Л 0
м =
— є
вокруг оси 0Х:
где
0кг (а) = Qx (а )М,
(1 0
Qx (а) =
0 'ї
0 соєа — 8Іп а 0 єіпа 008а
- матрица поворота вокруг оси ОХ на угол а. Тогда циклоидальная кривая на сфере (см. рис. 2) задаётся в виде вектор-функции
Бп (а) = Qz (а)Окг(а), (1.1)
где
Qz(а)=
sin
' —а і . f —
V n J —sin V:
^—аі f—1
I cos I
v n J V r
О О
—а
n
- матрица поворота вокруг оси OZ на угол
—а
Л
n
Рис. 2. Обозначено: 1 - циклоидальная кривая; 2 - окружность радиуса є,
3 - точка, описывающая циклоиду
Поверхность колеса образована движением в пространстве кривой, называемой профилем детали. В коническом ЭЦ-зацеплении профилем колеса является эквидистанта циклоидальной кривой на сфере (l.l), т.е. огибающая семейства окружностей радиуса р < є , лежащих на сфере, центры которых проектируются из
центра сферы в точки кривой (l.l). Для построения эквидистанты достаточно повернуть радиус-вектор каждой точки циклоидальной кривой в нормальных плоскостях этой кривой на угол
л/R^2
Y = arccos---------.
R
Принимая в качестве базисных векторов нормальной плоскости в точке кривой (l.l) ортогональные векторы Sn^) и векторное произведение Sn^)xSn'(а), запишем эквидистанту циклоидальной кривой на сфере в виде
~ЕГ, ч тт. ч , о • Sn(а)x Sn'(а)
Ev(а) = cos y Sn^) + R sin y-f=----=—г. (l.2)
|Sn(а) x Sn'(а)|
На рис. 3 изображен фрагмент эквидистантной кривой - профиля колеса в касании с окружностью радиуса р - профилем зуба шестерни для случая 21 = 3 .
'4-
2
Рис. 3. Обозначено: 1 - окружность радиуса е; 2 - окружности радиуса р;
3 - эквидистантная кривая
Для построения поверхности колеса нужно эквидистанту циклоидальной кри-
и Л7
вой поворачивать на угол — вокруг ОА с одновременным уменьшением радиуса
п
сферы, на которой она расположена от Я до (Я - 1г) (1г - ширина входной детали). Это приводит к записи уравнения поверхности колеса в виде вектор-функции двух аргументов:
Ек (и, а) = С (и )• Qz (и )£у(а), (1.3)
1г • х1
Я • 2п
где С (и) = 1 - с1и, с1 =
п 2п
а параметр и изменяется от 0 до —.
21
При и = 0 координатная линия поверхности Ек(и, а) лежит на сфере радиу-
о 2п Я - 1г
са Я , а при и =-------на сфере радиуса-----------.
21 Я
В [3] было показано, что точки контакта профилей деталей в реальном зубчатом ЭЦ-зацеплении лежат в достаточно узкой области в зоне перегиба циклоидальной кривой. Поэтому достаточно вытачивать только часть поверхности колеса, содержащую эту зону контакта. Эта часть поверхности определяется как лежащая по одну сторону от конуса с вершиной в центре сферы и основанием - ок-
ружностью радиуса rk на сфере. Величина rk задается из тех соображений, чтобы «обрезаемая» часть поверхности колеса не содержала точек зоны контакта. Запишем уравнения конуса, обрезающего зубья:
ЧR2 - rk2 cos v
Vr2 - rk2 sinv '
-rk у
На рис. 4 изображен фрагмент поверхности колеса и конус, «обрезающий» зуб детали.
Kon (u, v) = і
1.2. Входная деталь
Поверхность зуба входной детали получается следующим образом. Дана сфера радиуса Я и круговой конус с вершиной в центре этой сферы, пересекающий её по окружности радиуса є. Поверхность зуба образована окружностями, лежащими на концентрических сферах уменьшающихся радиусов, с центрами на конической винтовой линии, причем по мере приближения к центру сферы уменьшаются и радиусы окружностей (рис. 4). Наибольшая образующая окружность зуба лежит на сфере радиуса Я и имеет радиус р. Эту окружность можно задать в виде вектор-функции:
okr(a)=
Vr2 -e2^R2 -p2 +
ep sin a
v
R
у
Тогда уравнение поверхности зуба шестерни можно записать в виде
Е,?(г>, а) = С (и) Qx(v) окг (а). (1.4)
Уравнение конуса, на котором лежат центры окружностей, образующих зуб шестерни, имеет вид
коп (и, V) =
На рис. 5 изображён этот конус и один зуб шестерни.
и>ІЯ2 -є2 и -Є•008V
и - Є-ЄІПV
V /
ч.
Рис. 5. Поверхность одного зуба входной детали на коническом барабане
2. Определение местоположения сферической фрезы при обработке поверхностей деталей
Предлагаемый здесь алгоритм управления движением сферической фрезы подходит для обработки как входной, так и выходной детали. Поэтому в дальнейших рассуждениях будем пользоваться обозначением Е (и, а), считая Е (и, а) = Ек (и, а) (см. (1.3)), если речь идет об обработке поверхности колеса, и Е (и, а) = Еs( и, а) (см. (1.4)) - для шестерни. Фреза вытачивает деталь последовательным прохождением координатных линий поверхности детали, т. е. поверхность фрезы должна в каждый момент движения касаться наперед заданной координатной линии обрабатываемой поверхности. Обозначим через Яс1 и Яс0 (Ясг > Яс0) расстояния от центра сферы (фрезы) до оси вращения колеса в начале и в конце обработки, причем эти величины задаются так, чтобы разность Яс1 - Яс0 была несколько больше, чем ширина деталей 1г . Таким образом, в про-
цессе движения фрезы её центр находится на поверхности цилиндра переменного радиуса с осью - осью вращения колеса. Запишем уравнения этого цилиндра в виде вектор-функции:
'Rxk cos u Clk (u, v) = Rxk sin u
I v
Здесь через Rxk обозначен переменный радиус цилиндра:
( Rc - Rc0
k = RCl -( 1 K
где K - число шагов изменения радиуса цилиндра, а к - это целое число, находящееся в диапазоне от 0 до K.
Для каждого значения этого радиуса от Rc1 до Rc0 нужно найти координаты
центра сферы (фрезы), касающейся заданной линии на обрабатываемой поверхности. Линия задается фиксированным значением параметра а = а0 в уравнениях поверхности (1.3) или (1.4). Это касание обеспечивается, во-первых, наличием общей точки у сферы диаметром df и у поверхности детали, что можно в векторной форме записать в виде
df_
2
и, во-вторых, параллельностью нормалей к этим поверхностям в точке касания (общая касательная плоскость). Поскольку нормаль к сфере определяется вектором Clk (u, v) - F (и, а0), а нормаль к поверхности детали перпендикулярна касательным к двум координатным линиям на этой поверхности, то второе условие приводит к ещё двум уравнениям:
|С (u, v) - F (и, а0)]-д F (и а0) = 0; (2.2)
k
рГк (u, v) - F (и, a0)| = d-, (2.1)
до
С(и, V) - і(и, а0)]-д^(и,а0) = 0. (2.3)
л да
Таким образом, получена система (2.1) - (2.3) из трёх уравнений на три неизвестных: и, V, и. Её аналитическое решение приводит к трансцендентному уравнению на и, корни которого находятся численно с помощью встроенных программ пакета МаШСаі Благодаря специальному созданному алгоритму нахождения начальных приближений, корни находятся с высокой точностью (порядка 10-12мм). На рис. 6 изображено положение сферической фрезы с центром на заданном цилиндре при касании фрезы с поверхностью зуба колеса в точке заданной линии на этой поверхности.
3. Обработка деталей торической фрезой
При формообразовании поверхности детали малого размера, целесообразно использование не сферической фрезы, а фрезы, режущая часть поверхности которой - тор, поскольку при этом достигается меньший радиус кривизны вытачивающей поверхности фрезы.
Рис. 6. Касание сферической фрезы с поверхностью выходной детали в точке заданной линии. Показан цилиндр, на котором лежит центр фрезы
Зададим поверхность тора, полагая, что ось фрезы - это ось ОХ, в виде вектор-функции двух аргументов:
у -у оши у
где й/ - диаметр фрезы, г/ - радиус кривизны боковой поверхности тора.
3.1. Формообразование поверхности зуба колеса
При обработке поверхности зуба колеса фреза перемещается в пространстве, сохраняя положение своей оси параллельным оси вращения колеса (ось 02). При касании с поверхностью колеса в точке координатной линии, задаваемой фиксированным значением параметра а = а0 в уравнениях этой поверхности (1.3), положение фрезы определяется значениями трех смещений поверхности тора (3.1) относительно координатных осей. Обозначая величины этих смещений через Ях, Яу, Яі, запишем уравнения поверхности торической фрезы в касании с зубом колеса в виде
(3.1)
ТЯ(и, V)
/ - г/^ єіп V + г/ 008 и єіп V + Яу г/ єіп и + Яі
(3.2)
у
у
Определим дискретный набор точек на линии а = а0 значениями параметра и в уравнении поверхности колеса (1.3) по формуле
2пк
ик = — >
21К
где К - число шагов изменения параметра и, а к - целое число, изменяющееся от 0 до К. Задача определения местоположения фрезы при обработке зуба колеса сводится к тому, чтобы для каждого значения к, т.е. для каждой из К точек на линии а = а0 , найти величины смещений Ях, Яу, Яг, исходя из условий касания поверхностей фрезы (3.2) и детали (1.3). Эти условия: 1) наличие общей точки у поверхностей (3.2) и (1.3) и 2) параллельность нормалей к поверхностям в точке касания. Из первого условия, приравнивая правые части (3.2) и (1.3), получаем выражение смещений Ях, Яу, Яг через параметры и и у:
( Ях А
Ку
V К
= Ек (ц., а0) - іот (и, у).
(3.3)
Для записи соотношений, к которым приводит второе условие, обычным образом найдем вектор нормали к поверхности тора (3.1):
№ (и, у) =
діот(и, у) діот (и, у)
ди
ду
008 И БІПУ БІП И
\
Аналогично, дифференцируя по параметрам и и а вектор-функцию (1.3), найдем вектор нормали к поверхности зуба колеса:
№к (и, а) =
где обозначено:
/ (и,а)ооБ| и^ + / (и,а)БІп ^п
./1(и,а)БІп| П/,(и,а)ооБ^П /о(и,а)
/0(и, а) = ^Еу(а)0 Еу (а)0 + Еу(а) Еу (а)1) С(г)^ + ^Еу(а)0 Еу (а)1 - £У(а)1 Еу (а)01 с1
/1 (и, а) = -Еу(а)0 Еу (а)2 С(Ц^ + ^Еу(а)1 Еу (а)2 - Еу(а)2 Еу (а)1)),
/2 (и, а) = -Еу(а)1 Еу (а)2 С(Ц^ + ^Еу(а)0 Еу (а)2 - Еу(а)2 Еу (а)0)с1.
(3.4)
В соотношениях (3.4) скалярные функции Еу(а). и Еу (а). (і = 0,1,2) - это
координаты вектор-функции Еу(а) из (1.2) и ее производной по параметру а, а С (и) и с1 введены в (1.3).
Из условия 2) параллельности нормалей Ш(и,у) и №к(ик, а0) получаются два соотношения на параметры и и у, из которых удается выразить эти параметры че-
рез ик и а0. Таким образом, для каждой точки (ик) на линии а0 получаем значения параметров и и у, при подстановке которых в (3.2) получаем координаты точки касания фрезы и детали, а при подстановке этих значений в (3.3) находим величины смещений Ях, Яу, Яг, определяющих положение фрезы в пространстве в момент касания с деталью. На рис. 7 показан фрагмент поверхности зуба колеса в касании с торической фрезой в точке линии а = а0 .
Рис. 7. Фрагмент поверхности зуба колеса в касании с торической фрезой
3.2. Формообразование поверхности зуба шестерни
Запишем уравнения обрабатываемой поверхности зуба шестерни (1.4) в развернутом виде:
( окт (а)0 А
Fs(x>, а) = С (и)
окт(а)1 ооБ и - окт(а)2 БІп и окт(а)1 БІпи + окт(а)2 ооб и
(3.5)
где окг(а) ( = 0,1,2) - координаты вектор-функции окг(а), задающей наибольшую образующую окружность зуба шестерни (см. п. 1.2). Параметр а на поверхности зуба шестерни считается заданным (а = а0) - тем самым на этой поверхности выделяется винтовая линия, которой касается фреза при обработке детали. Вектор нормали в произвольной точке поверхности (3.5) имеет вид
' Я0(и,а) А
Жу(и, а) = я1(и, а)соБ и + §2 (и, а)81п и ^ (и,а)зши - g2 (и,а)соБ и
где обозначено
g0 (и,а) = окг '(а)1 (окг (а)1 С (и) - окг (а)2 с1) + окг '(а)2 (окг (а)2 С (и) + окг (а)1 с1), g1 (и,а) = -окг '(а)2 окг (а)0 с1 - окг '(а)0 (окг (а)1 С (и) - окг (а)2 с1), g2(u,а) = -окг '(а)1 окг (а)0 с1 + окг '(а)0 (окт (а)2 С (и) + окг (а)1 с1),
а окг'(а), / = 0,1,2, - координаты производной вектор-функции окг(а).
При обработке поверхности зуба шестерни ось фрезы находится в плоскости П, перпендикулярной оси вращения шестерни (ось ОХ), смещаясь вместе с этой плоскостью вдоль оси ОХ на величину Кх. Кроме того, ось фрезы поворачивается в плоскости П на угол ф, а поверхность тора смещается вдоль этой оси на величину Ш. Эти три величины определяют положение фрезы в каждый момент времени
обработки детали (т.е. для каждой точки Е^'(ик, а0) заданной линии а = а0 на обрабатываемой поверхности) и должны быть найдены из условия касания взаимодействующих поверхностей. Запишем уравнения поверхности фрезы с учетом этих смещений в виде
(1 0 0 А ( (0 А А (КхА
Tr(u, v, ф, Rx, Ht) =
tor(u,v)+
Q
V Ht УУ V
Q cos ф sin ф
VQ - sin ф cos фу
Тогда единичный вектор нормали поверхности тора в точке с криволинейным координатами (u, v) имеет вид
' cosucosv '
Nt(u,v,ф) = sinфsinu+cosфcosusinv .
V cos ф sin u-sin ф cos u sin v у
Чтобы найти точки касания поверхности детали и тора (фрезы), используем два условия: параллельность нормалей этих поверхностей в точке касания и принадлежность этой точки обеим поверхностям.
Из условия параллельности нормалей
Ns(uk, aQ) x Nt(u, v, ф) = Q
при a = aQ получаем 2 уравнения на З неизвестных, u, v, ф, из которых можно выразить u и v через ф:
sin(Uk +ф) gj (Uk, aQ) - cos(Uk +ф) g2 (Uk, aQ)
tg(u) =
tg(v) =
VgQ(uk, aQ)2 + gl(uk, aQ)2
cos(Uk +ф) g! (Uk, aQ) + sin(Uk +ф) g2 (Uk, aQ)
(З.б)
go(uk, а0)
Условие наличия общей точки у поверхностей:
Тг (и, V, ф, Кх, Ш) = -Р5(ик, а0)
приводит к ещё трём уравнениям на 3 неизвестных ф, Кх и Шt (после подстановки в эти уравнения соотношений (3.6)). Аналитическое решение этой системы приводит к трансцендентному уравнению на ф, корни которого находятся численно с помощью встроенных программ пакета МаШСа^ Благодаря специальному созданному алгоритму нахождения начальных приближений, корни находятся с высокой точностью (порядка 10-12 мм). На рис. 8 изображено положение торической
фрезы при её касании с поверхностью зуба шестерни в точке заданной линии на этой поверхности.
Рис. 8. Контактное взаимодействие поверхности зуба шестерни с торической фрезой
Решенные в данной работе задачи были использованы при разработке компьютерных программ для станка с ЧПУ, управляющих движением фрез при формообразовании деталей различных механизмов с ЭЦ-зацеплением. Эти программы успешно применяются при изготовлении в ЗАО «Технология маркет» (г. Томск) опытных образцов конических тяговых редукторов для железнодорожного транспорта и двухступенчатого редуктора [3], который прошел успешные испытания на фирме SEW (Германия).
ЛИТЕРАТУРА
1. Камчатный С.А., Сковородин А.В., Становской А.В., Щербаков Н.Р. Математическое моделирование процесса сближения поверхностей при формообразовании деталей передаточного механизма с ЭЦ-зацеплением // Вестник Томского государственного университета. Математика и механика. 2012. № 4(20). С. 15-24.
2. Патент РФ 2439401. Эксцентриково-циклоидальное зацепление зубчатых профилей (варианты) / В.В. Становской, С.М. Казакявичюс, Т.А. Ремнева, В.М. Кузнецов, А.В. Становской. Заявлено 29.01.2010; опубл. 10.01.2012, Бюлл. № 1.
3. Становской В.В., Казакявичюс С.М., Ремнева Т.А. и др. Двухступенчатый редуктор на основе эксцентриково-циклоидального зацепления (зацепление ExSyGear) // Вестник машиностроения. 2011. № 12. С. 41-44.
Статья поступила 14.03.2014 г.
Kamchatnyi S.A., Kuznetsov V.M., Skovorodin A.V., Scherbakov N.R. MATHEMATICAL SIMULATION OF THE FORMBUILDING OF THE CONIC DRIVING GEAR DETAILS WITH EC-GEARING. Formbuilding of the surface of a driving gear detail, in the most general formulation, is the task of two surfaces (the processed detail and the instrument’s surface) closing in up in the space to the point of contact where the surfaces have the common contact plane. Solving the problem of analytical description of the geometry of the surface contact is necessary for finding optimum methods and tools of processing details on multicoordinate NC units. In this
aspect, a mathematical model for the formbuilding of conic driving gear details with EC-gearing has been constructed in this work using milling cutters of two types. Systems of equations for the determination of the milling cutter location, at the point contact, have been derived. Analytical and numerical solutions of these systems are found. The complex of special computer programs to control the milling cutter moving while processing the detail was developed.
Keywords: formbuilding, eccentrically cycloidal (EC) gearing, contact of surfaces
KAMCHATNIY Sergei Alexandrovich
(M.Sc, Tomsk State University, Tomsk, Russian Federation)
E-mail: kam-serega2030@sibmail.com
KUZNETZOV Vladimir Mihailovich
(M.Sc , Tomsk State University, Tomsk, Russian Federation)
E-mail: tm@mail.tomsknet.ru
SKOVORODINAlexandr Vladimirovich
(M.Sc, Tomsk State University, Tomsk, Russian Federation)
E-mail: tm@mail.tomsknet.ru
SHCHERBAKOV Nikolay Romanovich (Doctor of Physics and Mathematics, Assoc. Prof.,
Tomsk State University, Tomsk, Russian Federation)
E-mail: nrs@math.tsu.ru
REFERENCES
1. Kamchatnyy S.A., Skovorodin A.V., Stanovskoy A.V., Shcherbakov N.R. Matematicheskoe modelirovanie protsessa sblizheniya poverkhnostey pri formoobrazovanii detaley pereda-tochnogo mekhanizma s ETs-zatsepleniem (2012) Vestnik Tomskogo gosudarstvennogo uni-versiteta. Matematika i mekhanika. No. 4(20), pp. 15-24. (in Russian)
2. Stanovskoy V.V., Kazakyavichyus S.M., Remneva T.A., Kuznetsov V.M., Stanovskoy A.V. Ekstsentrikovo-tsikloidal'noe zatseplenie zubchatykh profiley (varianty). Patent RF 2439401 / Zayavleno 29.01.2010; opubl. 10.01.2012, Byull. No 1. (in Russian)
3. Stanovskoy V.V., Kazakyavichyus S.M., Remneva T.A., Kuznetsov V.M., Bubenchikov A.M., Shcherbakov N.R., Shmidt Y. Dvukhstupenchatyy reduktor na osnove ekstsentrikovo-tsikloidal'nogo zatsepleniya (zatseplenie ExSyGear) (2011) Vestnik mashinostroeniya. No. 12, pp. 41-44. (in Russian)