Simulation of rolling bearings for the use of finite element analysis of machine tool spindles

Full Text

Определение динамических характеристик шпиндельных узлов металлорежущих станков является важным аспектом при проектировании новых конструкций шпиндельных узлов. На этапе проектирования одним из вариантов оценки собственных частот и форм колебаний шпиндельного узла является использование пакета Ansys Workbench, основанного на методе конечных элементов. Но при этом возникает проблема моделирования упругих свойств опор качения. Как известно, упругие свойства опоры качения определяются податливостью следующих соединений: наружное кольцо – корпус, тела качения – наружное кольцо, внутреннее кольцо – тела качения, вал – внутреннее кольцо; податливостью наружного и внутреннего колец и тел качения. Моделирование опоры качения на основе вышеперечисленных податливостей при динамическом анализе в Ansys Workbench в этом случае не представляется возможным, т. к. контактные взаимодействия не учитываются. Это приводит к необходимости искать другие подходы при моделировании упругих свойств опоры качения. Одним из таких подходов является моделирование подшипника в виде сплошного кольца, обладающего упругими свойствами, идентичными упругим свойствам реального подшипника. Таким образом, цель данной работы – разработка методики моделирования опоры качения посредством целого сплошного кольца, упругие свойства которого совпадают с упругими свойствами опоры качения. Поскольку в опоре необходимо учитывать одновременно радиальную и осевую жесткость, которые, как правило, существенно отличаются, то для сплошного кольца это возможно только путем использования в модели материала с ортотропными свойствами. При использовании такого материала при конечно-элементном моделировании необходимо задавать шесть независимых характеристик: три модуля Юнга и три модуля сдвига, подбор которых методом простого перебора требует существенных затрат времени. В связи с этим на первом этапе выбора свойств материала сплошного кольца необходимо разработать модель, определяющую упругие перемещения в материале кольца в зависимости от указанных выше шести независимых характеристик. В качестве такой модели может быть использован полином, полученный при полном факторном эксперименте. Был выбран двухуровневый полный шестифакторный эксперимент с одинаковыми диапазонами варьирования факторов. Диапазоны, верхний и нижний уровни принимались одинаковыми для всех факторов и определялись в ходе предварительных численных экспериментов по минимальному значению при сходимости матрицы жесткости. Проведение полного факторного эксперимента дает возможность получить полиномы перемещений в кольце для радиального и осевого направления. Общий вид полинома: Окончательный выбор упругих характеристик материала кольца производится путем сравнения значений упругих перемещений подшипника, рассчитанных по литературным источникам или полученных экспериментально, со значениями, определенными по полиному. Производится корректировка факторов в пределах диапазона варьирования с целью совпадения (с требуемой точностью) расчетных и полиноминальных значений перемещений. Покажем на примере роликового двухрядного конического подшипника с буртом на наружном кольце 4-697713Л, нагруженного последовательно силой 100 Н в радиальном и осевом направлении, методику выбора характеристик математической модели (табл. 1). Таблица 1 Исходные данные полного факторного эксперимента Варьируемые факторы Нулевой уровень Интервал варьирования Верхний уровень Нижний уровень Модуль Юнга () Натуральное значение, Па Кодированное значение 0 1 1 -1 Модуль Юнга () Натуральное значение, Па Кодированное значение 0 1 1 -1 Модуль Юнга () Натуральное значение, Па Кодированное значение 0 1 1 -1 Модуль сдвига () Натуральное значение, Па Кодированное значение 0 1 1 -1 Модуль сдвига () Натуральное значение, Па Кодированное значение 0 1 1 -1 Модуль сдвига () Натуральное значение, Па Кодированное значение 0 1 1 -1 Данный тип подшипника используется в передней опоре шпиндельных узлов станков модели 16Б16П и 16Б16Т1 Средневолжского станкозавода. На рис. 1 представлены чертеж и объемная геометрическая модель сплошного кольца и вала. Геометрические размеры кольца соответствуют габаритным размерам реальной подшипниковой опоры данного типоразмера. Материал вала – изотропный, Па, , материал кольца – ортотропный, , модули Юнга , , и модули сдвига , , – варьируемые в пределах указанных в табл. 1 диапазонов. а б Рис. 1. Геометрическая модель Геометрическая модель разбивалась на конечные элементы размером 0,008 м (рис. 2). Рис. 2. Конечно-элементная модель опоры качения Наружной поверхности кольца запрещались перемещения по всем трем координатам. Нагрузка величиной 100 Н прикладывалась симметрично к обоим торцам вала сначала в радиальном (рис. 3, а), потом в осевом (рис. 3, б) направлении. Упругие свойства кольца оценивались по максимальным перемещениям в радиальном и осевом направлениях. а б Рис. 3. Расчетные схемы для оценки радиальной (а) и осевой (б) деформаций опоры качения По результатам численного эксперимента построены полиномы, позволяющие вычислить перемещения вала относительно закрепленной наружной поверхности кольца при действии нагрузки как в радиальном, так и в осевом направлении (т. е. упругие деформации кольца). Анализ полученных полиномов показал, что значимость учитываемых факторов существенно отличается. Для проведения ранжирования факторов в полином подставлялись кодированные значения шести факторов. Причем значения одного из факторов варьировались во всем диапазоне, а величины пяти других факторов принимались равными нулевому уровню. По результатам проведенных расчетов построены следующие графики (рис. 4). Из графиков можно сделать следующие выводы. На радиальную податливость упругого кольца оказывают влияние четыре фактора: , , ,, причем для фактора зависимость прямая, т. е. при возрастании величины перемещения в направлении оси Y также увеличиваются, для остальных трех факторов зависимость обратная – увеличение , , приводит к снижению величины перемещений. Также следует отметить, что наибольшее влияние из всех факторов оказывает , а наименьшее . В результате получена следующая последовательность факторов, влияющих на радиальную податливость, в порядке убывания их значимости: , , , . На осевую податливость упругого кольца в основном влияют два фактора: и , остальные факторы имеют слабое воздействие на перемещения в направлении оси Х. Причем значимость обоих факторов одинакова. а б Рис. 4. Графики, используемые для ранжирования факторов Таким образом, можно сделать вывод, что на упругие свойства в радиальном и осевом направлениях влияют различные модули Юнга и модули сдвига. Это позволяет, варьируя те или иные модули Юнга и модули сдвига, добиваться независимых значений податливости в радиальном и осевом направлениях. На основании сделанных выводов построены упрощенные полиномы для задания упругих свойств ортотропного материала кольца для каждого из направлений. Перемещения по оси Y: Перемещения по оси Х: По аналитическим зависимостям [1] для подшипника 4-697713Л под действием радиальной и осевой нагрузок в 100 Н получены следующие значения упругих характеристик: радиальные перемещения – м, осевые перемещения – м. Анализ матрицы эксперимента выявил, что значения перемещений в радиальном направлении, полученные по аналитическим зависимостям, лежат внутри диапазона значений аналогичных перемещений, выявленных в ходе моделирования в Ansys. Таким образом, по сокращенному полиному для упругих перемещений по оси Y можно подобрать материал с такими же упругими свойствами, как и передняя опора шпиндельного узла токарного станка 16Б16П. В то же время значения осевых перемещений лежат выше диапазона значений, полученных при численном эксперименте (максимальное значение составляет м). Поэтому для получения необходимой величины осевой податливости следует снизить верхний и нижний уровни факторов, влияющих на осевую податливость (модулей сдвига и ), с одновременным сужением диапазонов их варьирования (табл. 2). Таблица 2 Натуральные и действительные значения модулей сдвига и Варьируемые факторы Нулевой уровень Интервал варьирования Верхний уровень Нижний уровень Модуль сдвига () Натуральное значение, Па Кодированное значение 0 1 1 -1 Модуль сдвига () Натуральное значение, Па Кодированное значение 0 1 1 -1 Сделанные уточнения верхних, нижних уровней и величины диапазона варьирования модулей сдвига и позволили получить сокращенные полиномы определения осевых и радиальных перемещений в следующем виде. Перемещения по оси Y: Перемещения по оси Х: . Так как полином расчета радиальных перемещений содержит четыре неизвестных, необходимо приравнять к нулю три из них, т. е. принять значение трех факторов на нулевом уровне. Четвертый фактор, оставшийся свободным, должен быть наиболее значимым. Далее методом подбора, используя лишь одну неизвестную, находим значение перемещений, наиболее близкое к значению, найденному по аналитическим зависимостям. В итоге мы получим величину наиболее значимого фактора при определенном значении радиальных перемещений. Если при этом не удалось добиться приемлемой сходимости между значением перемещений, полученным с использованием аналитических зависимостей и при расчете полинома, то необходимо перейти к варьированию второго значимого фактора. При этом в полином подставляется уже найденное ранее предельное значение первого значимого фактора, второй фактор варьируется, два остальных остаются на нулевом уровне, т. е. равны нулю. Полученные по полиному значения сравниваются с рассчитанными по формулам. Используя описанный выше алгоритм, добиваются максимального совпадения значений, при необходимости варьируя оставшиеся два фактора. По данной методике подобраны следующие свойства материала, моделирующего подшипник 4-697713Л: Па, Па, Па, Па, Па и Па Таким образом, получены упругие характеристики материала модели опоры качения в виде сплошного кольца, которые могут быть использованы при оценке динамических характеристик шпиндельного узла.

About the authors

Alexander F Denisenko

Samara State Technical University

Email: ask@samgtu.ru
(Dr. Sci. (Techn.)), Professor. 244, Molodogvardeyskaya st., Samara, 443100

Mihail V Yakimov

Samara State Technical University

Email: ask@samgtu.ru
Assistant 244, Molodogvardeyskaya st., Samara, 443100

References

Supplementary files