ПРОЕКТИРОВАНИЕ
УДК 621 833.7
Г. А. Тимофеев, М. В. Самойлова
КОЭФФИЦИЕНТ ПОЛЕЗНОГО ДЕЙСТВИЯ КОМБИНИРОВАННОГО ПЛАНЕТАРНО-ВОЛНОВОГО МЕХАНИЗМА
Рассмотрены особенности определения коэффициента полезного действия комбинированного планетарно-волнового механизма и приведена методика его определения.
E-mail: [email protected]; [email protected]
Ключевые слова: зубчатые механизмы, планетарные передачи, волновые
передачи, коэффициент полезного действия, работа, мощность.
Оригинальным конструктивным решением, в котором рационально сочетаются достоинства волновых и кривошипно-планетарных механизмов, является комбинированный планетарно-волновой механизм (КПВМ) [1-3].
Конструкция этого механизма (рис. 1) обеспечивает параллельное соединение кривошип-но-планетарного механизма и волновой зубчатой передачи (ВЗП) внешнего деформирования. Входной энергетический поток при передаче и преобразовании энергии в таком механизме разбивается на несколько потоков (при двух сателлитах — на четыре), а затем вновь суммируется на выходном звене. Это позволяет существенно повысить жесткость и точность механизма без существенного увеличения габаритных размеров и массы. При параллельной работе волнового и планетарного механизмов первый обеспечивает отсутствие люфта в передаче за счет податливости гибкого колеса. При малых нагрузках, когда люфт в планетарном механизме еще не
Рис. 1. Схема КПВМ
выбран, работает только волновая передача. С увеличением нагрузки в работу включается планетарная передача. Поскольку кривошипно-планетарная передача (КПП) имеет большие модуль и диаметр зубчатых колес, то и ее жесткость больше жесткости волновой передачи, поэтому при больших нагрузках основная их часть воспринимается планетарным механизмом, а меньшая — волновым.
Механический КПД КПВМ определяется соотношением
А А — А — А
-"■их -^д ^-тр -"■у.д
П = ~А~ =-А-!
Ад Ад
где Ад — работа движущих сил и моментов за цикл установившегося движения; Ап с — работа сил и моментов полезного сопротивления; Атр — работа, затраченная на преодоление сил трения в кинематических парах (КП) и сил сопротивления среды; Ау.д — работа, затраченная на упругую деформацию звеньев.
В случае установившегося режима работы механизма КПД можно определить через средние за цикл значения соответствующих мощностей:
П =
Nn.c = N - Np - N-у.д
N N
При наличии в механизме ряда КП мощность, расходуемая на преодоление сил трения, определяется по формуле
г=к
Nт ^ ^ Nmi,i—l (Nmol + Nml2 + •••) ,
г=1
где Nmi i_1 — мощность, расходуемая на преодоление сил трения в соответствующей КП.
Исследуемый механизм (см. рис. 1) имеет 18 КП [3], в которых энергия расходуется либо на упругое деформирование гибкого колеса (создание бегущей волны деформации), либо на преодоление трения во внутреннем эвольвентном зацеплении 1-й и 2-й зон кривошипно-планетарного механизма (КП Ь и N), в двух зонах зубчатого зацепления ВЗП Т) и зацепления гибкого колеса и волновой зубчатой муфты (Р, £), а также на трение скольжения и качения в генераторе волн ВЗП, в кинематических парах деформирующие шайбы-гибкое колесо (М, К), в подшипниках деформирующих шайб (Б, С, Е, Н, Г, 7), в опорных подшипниках валов (А, В, С, Я) и на трение гибкого колеса об опорные торцевые шайбы.
Как показали экспериментальные исследования КПД ВЗП [4], основные потери мощности на трение в ВЗП приходятся на генератор волн и зацепление (рис. 2). Доля остальных потерь в общем балансе потока мощности через ВЗП незначительна (менее 1 %).
Рис. 2. Схема распределения потоков мощности в КПВМ
Мощность потерь на трение в КП равна произведению силы (или момента) трения на скорость относительного движения. Силы, действующие в КП КПВМ, определены в работе [3].
Для определения мощности потерь на трение необходимо знать относительную скорость в исследуемой паре. С внутренними поверхностями деформирующих шайб (генератора волн) контактирует наружная поверхность деформируемого гибкого колеса. Для определения относительной скорости (скорости скольжения) необходимо найти скорости точек, лежащих на наружной поверхности деформированного гибкого колеса, и скорости точек деформирующих шайб.
Точки, лежащие на наружной поверхности деформированного гибкого колеса, отстоят от срединной линии на величину Нс/2 (Нс — толщина обода гибкого колеса под зубчатым венцом), их скорости определяются согласно работе [5]. Для определения скорости точек срединной линии гибкого колеса необходимо продифференцировать по времени аналитические выражения составляющих перемещений.
Радиальные, тангенциальные и угловые перемещения находятся по следующим зависимостям (рис. 3):
wi =
Wo
vi = +
A — Б
Wo
(A cos p — Б),
A-Б
(A sin p — Б p),
di =
Wo Б
rc r A - Б
sin(p)
Рис. 3. Расчетная схема зацепления в ВЗП внешнего деформирования
— для участка I (0 ^ у ^ в), гДе гибкое колесо принудительно деформируется генератором волн по дугам постоянного радиуса;
шп —
w 0
VII =
A-Б
wo
а -
sin ß) sin ф + ^ П2 — ^J cos ф — 2sin ß — Б
A- Б
— (2+sin2 ß) cos^+^-2—ф) sin ф+^^тв—ф)
1
- Wo
flII —--л с
Гот A — Б
2 cos ф — (2 sin ß + Б) ^П2 — ф)
— для участка II < У ^ п^)' где в — угловая координата участка постоянной кривизны; у — угловая координата рассматриваемого сечения, отсчитываемая от малой оси деформации; гс.г — радиус кривизны срединной линии недеформированного гибкого колеса; — радиальное перемещение по малой оси; эд, V и 9 — радиальные, тангенциальные и угловые перемещения соответственно.
Продифференцируем по времени аналитические выражения составляющих перемещений. Поскольку изменение угла у пропорционально скорости вращения деформированной средней кривой, то можно записать у = —1
Для первого участка
Vw1 —
dwI w0 A
"тг — ~A-F7 Шж sin Шж t
dt A — Б
радиальная скорость;
вю1
V^! —
w0
dt
A- Б
шж (A cos ш * t — Б)
тангенциальная скорость;
шв1 =
Oh ~dt
Wp Б Гс.г A - Б
Шж
— угловая скорость вращения нормали.
Скорость точки срединной линии деформированного гибкого колеса равна векторной сумме радиальной и тангенциальной составляющих:
/2,2 WP vi = \/< + =
Шж
VA2 + Б2 - 2АБ cos шжг.
Запишем аналогичные выражения для второго участка:
= _ Wo
vwi1 = - A- Б
Шж
sin2 ß • cos шж£ — — шж£^ sin шж£
=
Wo
A- Б
Шж
(2 + sin2 ß) sin шж£ +
+ — cos — 2sin в — Б
Положение точки K относительно осей координат определяется зависимостями [5]
I , hc
wik = wi +
Ihc
VIK = vi + у 0-
Тогда искомые скорости определяются из следующих выражений:
VL = ^ =
dw){ vj0A
W]
dt
A- Б
шж sin шж£,
V'
V]
d^i
~dt
Wo
A-Б
Шж
A cos шж£ — 1 +
2 Гс.г
Б
для первого участка;
dWI'I
v =-
dt
Wo
A-Б
Шж
ж
sin ß cos Шж£ — ( — — Шж£ I sin Шж£
2
V' = dfa = ^ + hc d^II V]I = d£ = d£ +2 d£ '
V' = Wo-Шж
Vv]I = A - Б
(2 + sin2 в) sin шж£ + ( — — шж£ J cos шж£ — 2sin ß — Б
+
h 11 f—^[2 sinШж£ — (2 sinß + Б) \
2 Гс.г
— для второго участка.
Скорость точек деформирующей шайбы определяют, исходя из того, что они совершают поступательное движение, в котором каждая
точка движется по круговой траектории. Скорости всех точек деформирующих кольцевых шайб одинаковы и равны
Уд = aw Шд,
где aw = e — эксцентриситет кривошипных валов, равный межосевому расстоянию передачи; шд — угловая скорость деформирующей шайбы.
Для сравнения скоростей точек контакта деформирующих шайб и гибкого колеса на участке I используются проекции вектора Уд на направления w и v в каждой расчетной точке. Для точки с угловой координатой р значения этих проекций вычисляются по формулам
VW = aw Шд sin уд = aw Шд cos
Скорость скольжения в зоне контакта между гибким колесом и внутренними поверхностями деформирующей шайбы определяется следующим образом (см. рис. 3):
Уск = Vi - .
На участке I, соответствующем зоне контакта, скорость скольжения постоянна.
Мощность потерь на трение в данной кинематической паре составляет
где — усилие на деформирующей шайбе, /т — суммарный коэффициент трения скольжения и качения.
Скорость скольжения Уск = УК2-К1 в точке контакта зубьев внутреннего эвольвентного зацепления K, необходимая для определения мощности потерь на трение, определяется из следующего уравнения сложного движения:
ук2 = ук1 - ук2-к1,
где Ук2 = ш21О2К — окружная скорость точки K2 колеса 2, жестко связанного с выходным валом редуктора (ш2 = шв м); УК1 — скорость точки K1 колеса 1 с внешними зубьями (деформирующей шайбы, на внешней стороне которой нарезаны зубья).
Для определения скорости Ук2 необходимо найти расстояние от оси вращения О2 до точки К. Согласно рис. 4, на котором представлена схема внутреннего эвольвентного зацепления, это расстояние можно записать как
lü2K = \JlKР + (lO2P)2 — 2lKPlO2P • COS
Расстояние от полюса зацепления Р до точки контакта K, расположенной по середине активной линии зацепления, определяется по
формуле
1 lBi Р + lB2 Р
1кр = 2 .
Расстояние от полюса зацепления Р до точек В1 и В2 находится по формулам
_ra2 sin Yi lBi Р = -:-;
sin у
Га i sin Y2
lB2 P —
sin y
где y* = 90o — aw; y = arcsin P s'—у; = 180o — y — y;
ra2
. l°!P sin y
У2 = arcsin-; 72 = 180 — y — У2; lo,р, lo2р — радиусы на-
rai
чальных окружностей rwi, rw2.
Скорость VKl = V = е, поскольку любая точка шайбы, как отмечалось ранее, движется по круговой траектории с радиусом е = aw.
Скорость скольжения находится из векторного плана скоростей КаЪ (см. рис.4)
Уск = Укг cos Щ + Vk2 cos Щ.
Согласно схеме, угол ZO2KNl = р + y — внешний угол AO2KP, тогда из условия перпендикулярности сторон имеем
Щ = ZO2KN1 = р + y.
Угол y определяется по формуле
. (Ikp . \
Y = arcsin I --sin р I ,
VO2P J '
угол Щ2 — из соотношения
, . (Ук2 ■ ■ {I02K U2 . , \ . lo2K sin Щ Щ2= arcsin ——sin Щ = arcsin -sinЩ = arcsin
Vk- J \ шле J e\UP|
Тогда скорость скольжения зубьев в точке K запишем как Кк = шже cos Ф2 + 1o2kcos Мощность потерь на трение в исследуемом зацеплении
N = F f V
1"т ± п,/т ' ск;
где Fjj — усилие во внутреннем эвольвентном зацеплении зубьев КПП; /т — коэффициент трения скольжения.
Расчет скорости скольжения зубьев V^ = VKг-Kж в точке K в ВЗП для определения N ведется аналогично расчету внутреннего эволь-вентного зацепления и определяется из уравнения сложного движения
Укг = Укж — Укг-к ж,
где VKx = шж1ОжК — окружная скорость точки K жесткого колеса, VKr — скорость точки KT гибкого колеса, деформируемого шайбой. Расчетная схема волнового зацепления приведена на рис. 5. Число зубьев условного колеса, имеющего геометрические параметры гибкого колеса на участке постоянной кривизны 2в, составляет
2aw cos aw
zY — ¿ж +
m cos a
Далее аналогично рассчитываем следующие параметры: радиус окружности вершин зубьев условного колеса
r ау r аж + h;
радиусы начальных окружностей колес
гжт cos a
2 cos aw
rwm — lo-ж P — rwy lOryP г,шж + a
Рис. 5. Расчетная схема зацепления в ВЗП
углы
у = 90° — aw;
. loжр sin у у1 = arcsin-,
Гаж
Yi = 180° — у — yi; расстояние от полюса зацепления Р до точки В1 как
sin Yi
lB1 Р = Гаж~-;
sin у
углы
. 1огуР sin у у 2 = arcsin —--,
Гау
Y2 = 180° — у — у2; расстояние от полюса зацепления P до точки В2
sin Y21
lB2 P — ray-
sin у
расстояние от полюса зацепления до точки контакта К в среднем положении
1В1Р + 1В2 Р
lKP —
2
Из расчетной схемы с использованием теоремы косинусов определяется расстояние
1ожк = ^ 12кр + (1ожp)2 - 21rp 1ожp cos у.
Скорость скольжения Уск находится из векторного плана скоростей Kab (рис. 5) как
Уск = Уж cos ^2 + Кг cos ,
где = /ОжК^ж.
Угол /ОжKN = у + y — внешний угол треугольника АОжKP. Из условия перпендикулярности сторон следует, что угол
= /Ож KN = у + y,
где y — arcsin
KP
sin у , а угол
, . ( Vk ж . Л
— arcsi n —— si n ^ —
VK
K г
. (1Ож K ^ж . , \
— arcsi n -si n ^ — arcsi n
V ^д1 /
1ожк sin e |Up|
Скорости скольжения зубьев в точке К
Уск = ^д e cos ^2 + 1ожк^ж cos . Мощность потерь на трение в таком зацеплении запишем в виде
N = FB3n /т уск J
где FB3n — усилие в зацеплении ВЗП внешнего деформирования; /т — коэффициент трения скольжения.
Аналогично определяется мощность потерь на трение в волновой зубчатой муфте.
Для определения потерь мощности на трение в подшипнике качения (рис. 6) условно отнесем радиальную нагрузку подшипника к одному шарику [6]. Мощность в этом случае определяется по формуле
N = Na + NB J (1)
Рис. 6. Модель подшипника качения
Ыа и Мв — мощность трения в точках контакта А и В, в которой
Ыл = ^ (шг + ш);
(2)
Ыв = Гцш,
где Г — радиальная нагрузка подшипника качения; шг — угловая скорость ведущего кольца I; ш — угловая скорость шарика; ц — коэффициент трения качения.
Поскольку точка А принадлежит звену 1 и шарику и имеет одну и ту же линейную скорость, то можно записать
д
ш12 = ш^ (3)
откуда
д
ш = Ш125'
где д — диаметр цапфы вала; Б — диаметр шарика.
Подставив выражения (2) в равенство (1), учитывая зависимость (3), получим выражение для определения мощности потерь на трение в подшипниках качения:
д
Мт = Гцшг (1 + 5
Увеличение числа КП повышает суммарные потери мощности на трение. Однако в многопоточном КПВМ, где поток мощности разветвляется с последующим суммированием на выходном валу, нагрузка в зацеплениях ВЗП и КПП уменьшается и соответственно уменьшаются потери мощности. По разработанной методике определены КПД, скорости скольжения и потери мощности на трение в трех редукторах, имеющих одинаковые нагрузочные характеристики и одинаковые передаточные отношения. Исследовался КПВМ в целом и его составляющие: ВЗП с генератором волн внешнего деформирования с муфтой и КПП. При расчете использовались эмпирические зависимости определения коэффициентов трения [7-9].
Для расчета использовались следующие исходные данные: 200 — передаточное отношение; 1,5 мм — модуль КПП; 0,4 мм — модуль ВЗП внешнего деформирования; 600 Н-м — крутящий момент на выходном валу; 200 с-1 — угловая скорость на выходном валу.
Результаты расчета приведены в таблице.
Отметим, что объединение ВЗП и КПП, хоть и увеличивает число КП, но уменьшает силы в зацеплениях волновой и планетарных передач. При этом уменьшаются потери мощности на трение. Расчетный КПД комбинированного механизма получается ниже чем у КПП, но
Результаты расчета КПД
Расчетные параметры КПВМ ВЗП КПП
Сила в волновом зубчатом зацеплении, Н 166,707 4010,72 —
Скорость скольжения, м/с 0,239 0,239 —
Потери мощности в зацеплении, Вт 1,994 47,97 —
Сила во внутреннем зубчатом зацеплении, Н 2039,88 — 2128,35
Скорость скольжения, м/с 0,3183 - 0,3183
Потери мощности в зацеплении, Вт 32,46 — 33,87
КПД 0,8 0,79 0,83
Рис.7. Зависимость КПД КПВМ от момента нагрузки
выше чем у ВЗП. На рис. 7 представлена зависимость механического КПД КПВМ от момента нагрузки на выходном валу.
СПИСОК ЛИТЕРАТУРЫ
1. А.с. № 1569467 А1 СССР, МКИ F16Н 1/00.
2. Тимофеев Г. А., Самойлова М. В. Использование метода графов в структурном анализе планетарно-волнового механизма // Вестник МГТУ им. Н.Э. Баумана. Сер. "Машиностроение". - 2010. - № 2. - С. 3-14.
3. Тимофеев Г. А., Самойлова М. В. Силовой расчет комбинированного планетарно-волнового механизма с генератором волн внешнего деформирования // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. - 2010. - № 4. - С. 17-22.
4. Финогенов В. А. Исследования КПД волновых зубчатых передач: Дисс. ... канд. техн. наук (05.02.02.). - М., 1971. - 150 с.
5. Тимофеев Г. А. Разработка методов расчета и проектирования волновых зубчатых передач для приводов следящих систем: Дисс. ... д-ра техн. наук (05.02.18, 05.02.02). - М.: МГТУ им. Н.Э. Баумана, 1996. - 352 с.
6. Соловьев А. И. Коэффициент полезного действия механизмов и машин. -М: Машиностроение, 1966. - 179 с.
7. Геккер Ф. Р. Динамика машин, работающих без смазочных материалов в узлах трения. - М.: Машиностроение, 1983. - 167 с.
8. Крагельский И. В., Добычин М. Н., Комбалов В. С. Основы расчетов на трение и износ. - М: Машиностроение, 1977. - 526 с.
9. Красавин С. И. Исследование масел и разработка методики их выбора для трансмиссий транспортных машин: Автореф. дисс. ... канд. техн. наук (05.05.03.). - М., 1979. - 16 с.
Статья поступила в редакцию 12.09.2011
Геннадий Алексеевич Тимофеев родился в 1944 г., окончил МВТУ им. Н.Э.Баумана в 1969г. Д-р техн. наук, заведующий кафедрой "Теория механизмов и машин" МГТУ им. Н.Э. Баумана, заслуженный работник высшей школы РФ, лауреат премии правительства РФ в области образования за 2008 г. Автор более 190 научных и методических работ в области автоматизированного проектирования механизмов машин.
G.A.Timofeev (b.1944) graduated from the Bauman Moscow Higher Technical School in 1969. D. Sc. (Eng.), head of "Theory of Mechanisms and Machines" department of the Bauman Moscow State Technical University, laureate of Government of the Russian Federation prize in the field of education for 2008. Honoured Higher School Worker of the Russian Federation. Author of more than 190 scientific and pedagogic-methodical publications in the field of automated designing of machine mechanisms.
Марина Валерьевна Самойлова родилась в 1962 г., окончила МВТУ им. Н.Э. Баумана в 1986 г. Канд. техн. наук, доцент кафедры "Теория механизмов и машин" МГТУ им. Н.Э. Баумана. Автор более 50 научных работ в области проектирования волновых и планетарных механизмов.
M.V. Samoilova (b. 1962) graduated from the Bauman Moscow Higher Technical School in 1986. Ph. D. (Eng.), assoc. professor of "Theory of Mechanisms and Machines" department of the Bauman Moscow State Technical University. Author of more than 50 publications in the field of design of the wave and planetary mechanisms.