DOI: 10.24412/2181-144X^024-2-67-75
Каримов К.А., Эгамбердиев И.П., Яхшиев Ш.Н., Абдуллаев К.С., Мамадияров А.Ж.
РАСЧЕТ ВЫХОДНОЙ ТОЧНОСТИ ШАРИКОПОДШИПНИКА ШПИНДЕЛЬНОГО УЗЛА ТОКАРНОГО СТАНКА НТ-250 М
Каримов К.А10009-0001-4183-4955], Эгамбердиев И.П20000-0002-7549-2156], Яхшиев Ш.Н. 3[0000-0003-3259-9735], Абдуллаев к.С.4[0009-0005-0924-572Х], Мамадияров А.Ж 5[0000-0002-8658-6166]
1 Ташкентский государственный технический университет имени Ислама Каримова kamolxon.karimov@gmaii.com,
2Навоийский государственный горно-технологический университет ilhom 1977@mail.ru, 3Навоийский государственный горно-технологический университет
sheralivaxshivev1978@mail.ru,
4ПО Навоийский машиностроительный завод KS.Abdullaev@ngmk.uz,
5Навоийский государственный горно-технологический университет
amamadivarov09@gmail.com.
Аннотация. Статья посвящена анализу динамических свойств шпиндельного узла. Нами было изучено упругие свойства опор и их математическое моделирование. Рассчитаны контактные и упругие деформации в роликовых подшипниках. Также рассмотрено моделирование контактно-упругих перемещений высокоточного роликоподшипника. На основе математической модели показано, что амплитуда колебаний частоты вращения шпинделя существенно влияет на качество работы. Показано нелинейное демпфирование подшипника качения из-за нелинейности коэффициентов изгиба и периодического изменения параметров.
Ключевые слова: Металлорежущие станки, зубчатые колеса, подшипник, долговечность, прогнозирование, техническое состояние, вибродиагностика, спектральный анализ, динамические модели.
Annotatsiya: Maqola shpindel uzelining dinamik xususiyatlarini tahlil qilishga bag'ishlangan. Biz tayanchlarning elastik xususiyatlarini va ularni matematik modellashtirishni o'rganib chiqdik. Rolikli podshipniklardagi kontakt va elastik deformatsiyalar hisoblanadi. Yuqori aniqlikdagi podshipnikning kontaktelastik siljishlarini modellashtirish ham ko'rib chiqiladi. Matematik modelga asoslanib, shpindelning aylanish chastotasi tebranishlarining amplitudasi ish sifatiga sezilarli ta'sir ko'rsatishi ko'rsatilgan. Dumalash podshipnikning chiziqsiz amortizatsiyalanishi, egilish koeffitsientlarining chiziqsizligi va parametrlarning davriy o'zgarishi ko'rsatilgan.
Kalit so'zlar: metall kesish dastgohlari, tishli g'ildirak, podshipniklar, chidamlilik, baholash, texnik holat, tebranish diagnostikasi, spektral tahlil, dinamik modellar.
Annotation. The article is devoted to the analysis of the dynamic properties of the spindle assembly. We have studied the elastic properties of the supports and their mathematical modeling. Contact and elastic deformations in roller bearings are calculated. Modeling of contact-elastic displacements of a high-precision roller bearing is also considered. Based on the mathematical model, it is shown that the amplitude of the spindle rotation frequency fluctuations significantly affects the quality of work. The nonlinear damping of the rolling bearing is shown due to the nonlinearity of the bending coefficients and periodic changes in the parameters.
Key words: Metal cutting machines, gears, bearing, durability, forecasting, technical condition, vibration diagnostics, spectral analysis, dynamic models.
Введение
Точность шпиндельного узла является одним из важнейших показателей точности станка, оказывающим непосредственное влияние на параметры обрабатываемого изделия. Методы контроля точности шпиндельных узлов по ГОСТ 22267—76, а также стандарты на нормы точности и жесткости для отдельных типов станков, предусматривают комплекс проверок, результаты которых дают возможность
оценить только геометрические параметры станка. Между тем известно, что такие проверки недостаточны для обоснованного заключения о выходной точности станка. Так, преувеличении радиального биения шпинделя точность обработки может в ряде случаев не ухудшаться.
В связи с этим целесообразна разработка такой системы расчетов, которая непосредственно связывала бы погрешности элементов шпиндельного узла с погрешностями обрабатываемых поверхностей. При этом шпиндельный узел следует рассматривать как часть формообразующей системы станка, а все ее звенья и их относительные перемещения (кроме элементов шпиндельного узла) — как абсолютно точные.
Методология и методы
Здесь такая задача решается на примере анализа влияния точности шпиндельных опор на комплекс выходных параметров станка: геометрическую точность, отклонения формы, положения и размеров обработанной поверхности. Цель анализа — разработка требований к точности опорных подшипников, а также установление взаимосвязи между входными и выходными погрешностями станка для типовых схем токарной обработки. Разработана математическая модель точности шпиндельного узла, которая учитывает геометрические и статические погрешности элементов опор. В основу построения модели положены формулы, в которых все входные погрешности, кроме погрешностей шпиндельного узла, приняты равными нулю. При построении расчетной модели опоры шпиндельного узла представляются системой пружин. Тогда шесть составляющих погрешности положения системы координат шпинделя могут быть получены из шести условий статического равновесия, которые в соответствия с [1] имеют вид
СД = Р. (1)
Силы на тела качения при радиальной нагрузке подшипника распределяются неравномерно (рис.1). В восприятии нагрузки участвуют только тела качения, расположенные на дуге, не превышающей 180° (нагруженная зона. Наиболее нагруженным является шарик, расположенный по направлению действия силы на подшипник.
Задача распределения сил между телами качения является статически неопределимой.
Тела качения, расположенные симметрично относительно плоскости действия силы, нагружены одинаково. Обозначим силу на наиболее нагруженное тело качения через Fo; на тело, расположенное по отношению к плоскости нагрузки под углом у (равным угловому шагу), - через Fi, под углом 2у -через F2, под углом пу - через Fn, где n- половина тел качения в нагруженной зоне.
Предполагаем для простоты, что тела качения располагаются симметрично относительно плоскости нагузки.
Условие равновесия внутреннего кольца, нагруженного радиальной силой Fr:
Fr = F0 + 2F1 cosy + 2F2 cos2y +... + 2Fn cosny (2)
В дополнение к уравнению статики используем уравнение и перемещений. Пренебрегая изгибом колец и предполагая отсутствие радиального зазора в подшипнике, можно принять, что сближение тел качения и колец равны соответствующим проекциям полного смещения кольца öo, т.е.
= ö0 cosy; ö2 = ö0 cos2y; 8i = S0 cosiy; (3) где i- номер тела качения.
Для шарикоподшипников зависимость между б шариков и сжимающей нагрузкой F, как следует из задачи теории упругости о сжатии упругих тел,
S = cF3;
(4)
где с - коэффициент пропорциональности.
Нелинейный характер зависимости объясняется ростом площадки контакта с ростом силы.
Рисунок 1. Схема распределения сил на подшипник.
Выражая в уравнениях перемещений сближения через силы, можно написать:
F = F cos/2 y F = F cos/2 2y F = F cos/2 iy , ,■ ■ ■,
Подставив эти зависимости в уравнение равновесия, получим
í n с Л
(5)
F = F
1 + cos52 iy (6)
V 1 У
Отсюда определяем F0, причем одновременно умножаем числитель и знаменатель правой части на z и вводим обозначение
к =-(7)
1 + 2Х«
. + 2 У cos2iy 1
Тогда
F =
k ■ F
(8)
где z - общее число тел качения.
Для подшипников с числом шариков z=10.. .20 k=4.37±0.01.
В подшипниках, эксплуатируемых в обычных условиях с зазором, в работе принимают участие шарики на дуге, меньшей 1800, и наиболее нагруженный шарик сжимается с силой, большей примерно 10% чем сила по этому расчету. В соответствии с этим для однорядных шарикоподшипников принимают k=5 и Fo=5Fr/z.
В сферических двухрядных шарикоподшипниках с учетом некоторой неравномерности распределения сил между рядами силу на наиболее нагруженный шарик оценивают Fo=6Fr/z(cosa),
2
Где а-угол наклона линии контакта, z-число шариков в обоих рядах.
Распределение нагрузки между телами качения может быть несколько выравнено упругими деформациями корпусов. Отверстие должно принимать форму эллиптического цилиндра, вытянутого в направление нагрузки. Это возможно при проектировании букс железнодорожного подвижного состава.
Осевая сила при точном изготовлении и отсутствии взаимного перекоса колец распределяется между телами качении равномерно.
Резултаты
Расчет без учета влияния центробежных сил шариков
Диаметр наружной беговой дорожки
Ъ = Бв + 2й
Приведенная кривизна внутреннего контакта
2 2 =2+ъ
Наружного контакта
V-. 2 2
= 2+Ъ
Скорость качения в контактах
— Ъ + 2ё (и + и ) = — • Ъ ъв-
60 Ъ + а
(9)
(10)
(11)
(12)
Нулевое приближение для усилия в контактах, приходящегося на единицу контактной длинны шарика
Ро =
Ъ 408
61-
14,832 - lg£pe + ^рн)
>.б
Толщины смазочных слоев в контактах
к - (U + Ub )]0-75 .«0
\Ре(н)г ]015 ■[^рв(н )
B.
(н )i
К(н) = 0796^ , , = [- 10-15
[Рв(нУ ] -\^Рв(н)] [рв(н)г ]
(13)
(14)
Рисунок 2. Усилия в контактах без учета центробежных сил шариков.
61-
+ K1+ К I10
15,7 - lg((P,-, ЕРв)
(15)
Здесь s^ + h + h
л в'-1 '
суммарный натяг по шарику с учетом влияния
расклинивающего действия смазочного слоя в контактах. Расчет усилий ведется методом инерции по (15) с одновременным уточнением толщин слоев по (14) 1=1.2.3.4...
Р — Р
до Р-^ < 0,00001.
Р
1 2.1 3.2 и
Рисунок 3. График зависимости контактных деформаций шарика в зависимости от его нагружения.
Рисунок 4. График распределения контактных деформаций на наружной
обойме подшипника.
Производим статический расчёт балки с одной промежуточной опорой под действием сосредоточенных сил, действующих на наружную обойму подшипника, расположенных в вертикальной плоскости.
чЬ-7-КЙ/
503.534Экд тт -504.1368кд тт
|-7.4724кд |-325.0863кд 17.441 Зкс
Привязки У-э!4--£-^-к-^-^--^-
О 132 154 1 76 198 220 242 264 286 308 441 Пролеты ]с_¿20_±_220_^ тт
Рисунок 5. Геометрическая схема балки с установленными нагрузками.
Граничные условия
Граничные условия на каждом краю могут быть:
•жёсткая заделка;
•шарнир;
© International Journal of Advanced Technology and Natural Sciences Vol. 2(5), 2024 SJIF IF=4.023
---1 ИИЧЧМ ичкгконндл ---— я
Go о ale d — ф 71
■свободный край.
Соотношения между Прогибом, углом поворота сечения и эпюрами M и Q Из курса «Сопротивление материалов» известны очень интересные соотношения, а именно:
угол поворота сечения U есть производная от вертикального перемещения w •изгибающий момент в сечении M есть производная от угла поворота U, умноженного на EJx;
•перерезывающая сила в сечении Q есть производная от изгибающего момента
M;
•распределённая нагрузка q есть производная от перерезывающей силы Q Таким образом, имеют место равенства:
U(z) = w'(z); (16)
M(z) = EJx • U'(z); (17)
Q(z) = M'(z); (18)
q(z) = Q'(z). (19)
При построении эпюр будем руководствоваться формулами (41-44). Балка является статически неопределимой: неизвестные реакции промежуточных опор и силовые факторы на концах балки не могут быть найдены из уравнений статики. Для решения задачи применён метод начальных параметров. Дифференциальное уравнение изогнутой оси балки при этом имеет вид:
EJMz) = EJxWo + EIAz + + ЪьMk tf + Lk<zFk ^ + Ц<гtf - ^^
(20)
где:
EJxwü - прогиб в левом сечении (с точностью до множителя EJx); EJxQü - угол поворота левого сечения (также с точностью до множителя EJx); Mo и Qo - изгибающий момент и перерезывающая сила в левом сечении. Все эти параметры (они называются начальными) - неизвестны. В каждой из сумм суммирование проводится по всем силовым факторам, расположенным слева от текущего сечения. Во второй сумме (сосредоточенных усилий Fk) учитываются неизвестные реакции опор R1, R2, .... Таким образом, в уравнении (45) - n+4 неизвестных, где n - количество промежуточных опор. Если все эти неизвестные будут найдены, то будет возможно построить эпюру перемещений по формуле (45) и другие эпюры с помощью производных от (45), которые по (41-43) дают: углы поворота:
Е]хв(2) = EJA +M0z + Q0j + lXk<z Mk(z - хк) + lXk<z Fk ^ + Ъхп<2 qk - Ц<2 qk ^ (21) изгибающие моменты:
M(z) = M0 + Q0z + lXk<2 Mk + lx< Fk(z - xk) + lxn qk ^ - Ц qk ^ (22)
перерезывающие силы:
Q(z) = Qo+ lXk<z Fk + Xxn<z qk(z - x£) - Xxm<z qk(z - (23)
Для нахождения этих n+4 неизвестных имеется столько же уравнений:
• под каждой опорой перемещение равно нулю - всего n уравнений вида (24) в тех точках, где расположены опоры;
• по 2 каких-либо параметра на каждом краю балки равны нулю - всего 4 уравнения.
В зависимости от вида граничных условий будут равны нулю:
• в жёсткой заделке - перемещение и угол поворота;
• при шарнирном закреплении - перемещение и изгибающий момент;
• на свободном краю - изгибающий момент и перерезывающая сила. На Рис.1 составляется n+4 = 3+4 = 7 уравнений с 7-ю неизвестными.
Мы составили и решили эту систему уравнений методом Гаусса с выбором ведущего элемента по каждому столбцу.
После того, как составили и решили систему уравнений, описанных в предыдущем разделе, и вычислили все необходимые данные для формул (45-48), строим эпюры перемещений, углов поворота, изгибающих моментов и перерезывающих сил, при этом множитель EJx вычисляется в зависимости от текущего элемента сортамента. Зависят же от этого множителя, как видно из уравнений (45-48) только перемещения, то есть прогибы и углы поворота.
Рисунок 6. Эпюра углов поворота [град-2].
Рисунок 7. Эпюра изгибающих моментов [kg^mm].
Рисунок 8. Эпюра поперечных сил [kg].
Рисунок 9. Результаты в разных значениях
Результаты расчёта:
Характеристики элемента: Сортамент: Полоса 7х77 Масса 1 м.п. = 4,21 кд Момент инерции, Jx = 712,0000 ст4 Момент сопротивления, Wx = 102,0000 ст3 Статический момент полусечения, Sx = 58,4000 ст3 Марка стали - ШХ15 ГОСТ 801 -78 Расчётное сопротивление стали, Ry = 230 МПа Расчётное сопротивление стали сдвигу, Rs = 0,58^у = 133,40 МПа Относительный прогиб - 1/250 пролёта Модуль упругости, Е = 211000 МПа Проверка условий прочности и жесткости: Напряжения в балке: нормальное = Мтах / Wx = 0,5887 < Ry = 230 МПа касательное = Qmax•Sx / дхЧст) = 3,1417 < Rs = 0,58^у = 133,4 МПа Максимальный прогиб (с коэффициентом надёжности) = 0,0347 тт-3, что составляет 1/6373316 от максимального пролёта 221 тт.
Заключение
Упругая система шпиндельной группы включает изгибную систему заготовки, шпинделя, опор шпинделя, приспособлений для крепления заготовки (патрона, центров и т.п.) и крутильную систему (точнее изгибно-крутильную) передачи от двигателя к заготовке.
Степень влияния нагрузки от передачи на деформацию шпиндельного узла зависит от расположения передачи по отношению к опорам, от жесткости опор и их защемляющего действия, а также величины консоли, на которой рассматривается точка контакта резца и обрабатываемой заготовки.
Большая доля влияния изгибов валов в связи с их значительной длинной, необходимых для расположения и передвижения шестерен. Существенно специфической является большая доля контактных деформаций в шлицевых соединениях в общем балансе, что связанно с применением в большом количестве передвижных шестерен. Таким образом шпиндель в процессе вращения меняет свое ориентирование относительно режущего инструмента.
В работе подшипников, которые применяются в качестве опор шпинделя, важное значение имеет обеспечение надежности. Ввиду быстроходности подшипников даже незначительные ненормальности в их работе приводят к быстрому разрушению подшипника. Обеспечение надежной работы подшипников невозможно без выполнения обширного комплекса мер. В их числе оптимальное конструирование подшипника занимает ведущее место. Однако необходимо констатировать серьезное отставание в этом вопросе. Передовые зарубежные фирмы уже давно перешли на проектирование подшип -ника «по месту», т.е. на проектирование с учетом конструктивных особенностей узла, теплового режима, условий подвода смазки, удобства сборки узла и др. Мы же вынуждены практически перетяжелять конструкцию, применяя «стандартные» подшипники, которые часто в конкретных условиях работы недостаточно надежны с точки зрения оптимальных соотношений параметров желобов.
В данной работе авторами описан расчет контактных и упругих деформаций в роликовых подшипниках, а также рассчитаны контактно-упругие перемещения
высокоточного роликового подшипника марки 2-697920Л2, установленного в шпиндельном узле токарного станка модели НТ- 250И.
Решетов Д.Н. выдвинул гипотезу о том, что отверстие должно принимать форму эллиптического цилиндра, вытянутого в направлении нагрузки. Однако конкретной формы предложено не было. Нами предложена форма отверстия. Геометрия которого была спроектирована с учетом деформации подшипника. В результате удалось минимизировать величины контактных и упругих перемещений, контактных перемещений в два раза, а упругих перемещений в 3 раза. Таким образом, для деталей диаметром до 50 мм появляется возможность достижения 8-9 квалитета точности.
Список использованные литературы:
[1]. Кудинов В.А. Динамика станков. М.: Машиностроение, 1967. 359 с.
[2]. Пуш В.Э. Металлорежущие станки. - М.: Машиностроение, 1985. 256 с.
[3]. Пуш А.В., Зверев И.А. Шпиндельные узлы: Проектирование и исследование. //-М., «Станкин», -2000 г. -С. 197.
[4]. Зенкевич О. Метод конечных элементов в технике. -М.: Мир, 1985 -542 с.
[5]. Атауллаев А.О., Яхшиев Ш.Н., Мамадияров А.Ж. Методика расчета динамических характеристик шпиндельного узла Развитие науки и технологии БМТИ Узбекистан. ISSN 2181-8193-2020 г.
[6]. Sh.N. Yaxshiyev, Kh. Kh. Ashurov, A. J. Mamadiyarov. Dynamics of spindle assembly of metal-cutting machine. International Journal of Engineering and Advanced Technology (IJEAT) ISSN: 2249 - 8958, Volume-9 Issue-3, February 2020