Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

MODELING AND STUDY OF STRESS-STRAIN STATE OF AXIAL CUTTING TOOLS WITH NUMERICAL FINITE ELEMENT METHOD IN CAD/CAM/CAE/PDM-SYSTEM SOLIDWORKS

Lukina S.V. 1 Krutyakova M.V. 1 Rybakova M.R. 1
1 Moscow State university of mechanical engineering (MAMI) (University of mechanical engineering)
Представлена методика конечно-элементного моделирования и исследования напряженно-деформированного состояния цельных и сборных осевых режущих инструментов, позволяющая прогнозировать жесткость инструментов на стадии технологической подготовки производства. Моделирование конструкций инструментов производилось созданием 3D-геометрических прототипов в CAD/CAM/CAE/PDM-системе SolidWorks. Расчетная параметризация конструкций осевых инструментов из 3D-геометрических прототипов реализована в Visual Basic 6.0. В расчетной схеме оценки напряженно-деформированного состояния учитывались силы крепления режущих элементов к корпусу инструмента, силы закрепления инструмента на станке и силы резания. Для оценки напряженно-деформированного состояния инструментов численным методом конечных элементов сформирована аналитическая расчетная модель, учитывающая деформации геометрических прототипов и жесткость стыков между ними. Разработанная модель реализована в среде Borland Delphi с использованием инструментов системы SolidWorks.
The technique of finite element modeling and investigation of the stress-strain state of the whole and partial axial cutting tools, allowing to predict the stiffness of the tools at the stage of the technological preparation of the production. The structural modeling of the tools is made by creating prototypes of 3D- geometry in CAD/CAM /CAE/PDM- system SolidWorks. The estimated parameterization of the designs axial tools of 3D-geometry prototype is implemented in Visual Basic 6.0. The fastening strength of the cutting elements to the tool body as well as the strength of the fastening of the tool in the machine and the cutting power were taken into account in the designed scheme of the estimation of the stress-stain state. To evaluate the stress- strain state of the numerical tools using finite element an analytical simulation model was formed which takes into account the deformation of geometric prototypes and stiffness of joints between them. The developed model is implemented in Borland Delphi environment using the tools of SolidWorks.
axial cutting tools
parametric model
stress- strain state
finite element method
1. Lukina S.V. Modelirovanie naprjazhenno-deformirovannogo sostojanija sbornogo rezhushhego instrumenta v processe ego izgotovlenija i jekspluatacii. Avtomatizacija i sovremennye tehnologii. 2001, no. 4, pp. 27–31.
2. Lukina S.V. Povyshenie jeffektivnosti proektirovanija sbornogo rezhushhego instrumenta na baze ustanovlennyh vzaimosvjazej konstruktorsko-tehnologicheskih i jekonomicheskih reshenij: dis. ... dokt. tehn. nauk. Moscow, 1999, pp. 180–220.
3. Lukina S.V., Krutjakova M.V., Kuprijanova O.P., Vjalkova P.F., Proektirovanie protjazhek optimalnyh konstrukcij v srede Microsoft Excel. Obrabotka metallov: tehnologija, oborudovanie, instrument, 2010, no. 1, pp. 24–26.
4. Lukina S.V., Makarov V.M., Krutjakova M.V., Issledovanie deformacij zubev kruglyh protjazhek metodom konechnyh jelementov. STIN, 2005, no.11, pp. 11–14.
5. Lukina S.V., Feofanov A.N., Krutjakova M.V., Rybakova M.R., Modelirovanie i issledovanie naprjazhenno-deformirovannogo sostojanija protjazhnogo instrumenta v CAD/CAM/CAE/PDM-sisteme SolidWorks. Tehnologija mashinostroenija, 2011, no. 12, pp. 57–59.
6. Segerlind L. Primenenie metoda konechnyh jelementov. Per. s angl. Moscow, Mir, 1979. 392 p.

Осевые режущие и комбинированные инструменты для обработки внутренних поверхностей, такие как сверла, развертки, прошивки, накатки при своей работе испытывают сложное напряженно-деформированное состояние от сил резания и сопротивления обрабатываемого материала. Одним из важнейших показателей, характеризующих долговечность работы инструмента, является его жесткость. Под действием сил резания зубья режущих инструментов упруго деформируются, в результате чего изменяются их форма и геометрические размеры. Элементы режущих кромок перемещаются, оказывая тем самым влияние на точность обработанной поверхности. Тело инструмента теряет устойчивость, следствием чего является увод оси инструмента. При обработке длинных отверстий, вследствие большой длины вылета инструмента и его низкой жесткости, из-за потери продольной устойчивости возможна поломка инструмента в процессе резания и, как следствие, отбраковка обрабатываемой детали. В этой связи задача оценки жесткости осевого режущего инструмента на начальных этапах проектирования операции металлообработки является актуальной и требует детальной проработки.

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

Расчетная параметризация конструкций осевых инструментов из 3D-геометрических прототипов реализована с помощью языка Visual Basic 6.0, позволяющего осуществлять программную интеграцию с системой solidworks (рисунок).

pic_28.tif

Расчетная параметризация конструкций осевых инструментов из 3D-геометрических прототипов

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

В качестве базовых конечных элементов для аппроксимации макроэлементов конструкций инструментов были выбраны объемные стержневые и тетраэдральные конечные элементы. Тетраэдральные конечные элементы использовались для аппроксимации режущих элементов инструмента, а стержневые – для аппроксимации корпуса инструмента. Поверхности контакта сборочных элементов моделировались множеством двумерных поверхностных стержневых и треугольных конечных элементов.

Произведенная конечно-элементная оценка напряженно-деформированного состояния параметризованной конструкции осевого инструмента учитывает деформации как внутри элементов конструкций инструментов, так и стыках между ними.

Общее разрешающее соотношение метода конечных элементов оценки объемного напряженно-деформированного состояния осевого инструмента [2, 6]:

[K]{δ} = {F}, (1)

где [K] – глобальная матрица жесткости; {δ} – вектор перемещений узлов конструкции инструмента; {F}– система внешних возмущений, отнесенных к конструкции инструмента.

Для тетраэдрального конечного элемента определен вектор-столбец перемещений узлов

Eqn42.wmf,

где u, v, w – составляющие перемещений узла тетраэдра по осям координат.

Матрица жесткости тетраэдрального конечного элемента [1]:

Eqn43.wmf (2)

где V – объем тетраэдрального элемента; [B] – матрица градиентов, выраженная через координаты узлов тетраэдра [1]; [D] – матрица упругих констант, характеризующая модуль упругости и коэффициент Пуассона материала тетраэдра [6].

Матрица жесткости геометрических прототипов элементов конструкции осевого инструмента представляет собой сумму матриц жесткости тетраэдров:

Eqn44.wmf (3)

Здесь Ne – количество тетраэдральных конечных элементов в конструкции геометрического прототипа.

Для стержневого конечного элемента определены перемещения и повороты его узлов в виде вектор-столбца:

Eqn45.wmf (4)

где Eqn46.wmf – перемещения в узлах стержневого элемента; Eqn47.wmf – углы поворота в узлах стержневого элемента относительно осей координат XYZ.

В общем случае перемещения в узлах стержневого элемента были разделены на три группы, определяющие три типа решаемых задач: изгиб в плоскостях XOY и XOZ (перемещения от изгиба внутри элемента в плоскостях XOY Eqn48.wmf и XOZ Eqn49.wmf, соответственно); осевые нагружения вдоль оси элемента ij (осевые перемещения по оси X элемента Eqn50.wmf); воздействие крутящего момента Mx (угловые повороты вокруг оси X элемента Eqn51.wmf).

Соответственно, обобщенная матрица жесткости стержневого элемента определялась суммой матриц жесткости по выделенным типам решаемых задач:

Eqn52.wmf (5)

Например, матрица жесткости стержневого элемента от сил, действующих в плоскости XOY и моментов Mzi и Mzj, определена интегралом:

Eqn53.wmf (6)

Для расчета матрицы жесткости стержневого элемента от воздействия крутящего момента Mx предварительно был определен угол закручивания элемента:

Eqn54.wmf (7)

Eqn55.wmf (8)

где Ge – модуль упругости элемента второго рода; Jp – полярный момент инерции.

Жесткость стержневого элемента определялась дифференцированием соотношения, определяющего потенциальную энергию деформирования, с учетом работы, совершаемой моментами Мxi и Мxj, прикладываемыми к узлам:

Eqn56.wmf (9)

Eqn57.wmf (10)

Матрица жесткости звена сборного осевого инструмента, расчлененного на стержневые конечные элементы, определяется суммированием матриц Eqn58.wmf по числу конечных элементов:

Eqn59.wmf (11)

Для формирования контакта между геометрическими прототипами и определения величин натягов и зазоров между соприкасающимися поверхностями элементов, имеющими реальный рельеф, в конечно-элементную модель инструмента было добавлено множество поверхностных элементов. Поскольку контакт реальных поверхностей прототипов является случайной величиной, характеризующей состояние поверхностей, положение физических точек контакта и величин зазоров, решение задачи объединения конструкции инструмента в единое целое осуществлялось множеством фиктивных контактных стержневых элементов единичной длины [1, 2, 4].

Матрица жесткости контактного стержневого элемента:

Eqn60.wmf (12)

где Eqn61.wmf – жесткость фиктивного стержневого элемента; [Т] – матрица перехода от местной системы координат поверхностного элемента XYZs к глобальной системе координат инструмента XYZ.

Eqn62.wmf (13)

Eqn63.wmf (14)

Жесткость фиктивных стержневых элементов подбиралась с учетом данных работ [1–6].

Жесткость стыка сборочных элементов определялась суммой жесткостей конечных элементов, моделирующих стык:

Eqn64.wmf (15)

Здесь Nk– количество конечных стержневых элементов, моделирующих стык.

Количество конечных стержневых элементов, моделирующих стык, принималось равным количеству узлов поверхностных линейных конечных элементов прототипа.

С учетом (3), (11), (15) глобальная матрица жесткости инструмента определялась по выражению:

Eqn65.wmf (16)

где Nc – общее количество стыков в конструкции инструмента.

Перемещения узлов объемных и поверхностных конечных элементов определялись решением управления (1) методом Гаусса.

Деформации {ε} и напряжения {σ} элементов инструмента рассчитывались по закону Гука [6]:

{ε} = [B]{U}; (17)

{σ} = [D]{ε}. (18)

Главные напряжения σ1, σ2, σ3 определялись из уравнения:

Eqn66.wmf (19)

где

Eqn67.wmf (20)

Эквивалентные напряжения рассчитывались по выражению [1]:

Eqn68.wmf (21)

Разработанная математическая модель была реализована в среде Borland Delphi с использованием инструментов системы Solid Works.

Для оценки адекватности разработанной модели была произведена серия численных экспериментов по исследованию напряженно-деформированного состояния протяжек и сборных комбинированных осевых инструментов (рисунок). Результаты расчетов позволили оценить влияние конструктивных и геометрических параметров протяжек и разработать рекомендации по их проектированию [3, 5].

Выводы:

  • в результате проведенных исследований сформирована методика конечно-элементного моделирования напряженно-деформированного состояния цельных и сборных осевых режущих инструментов, позволяющая оценить жесткость инструмента на начальных этапах проектирования операции металлообработки, обоснованно назначить конструктивные и геометрические параметры инструмента и режимы резания и тем самым повысить эффективность операции металлообработки;
  • разработанный алгоритм расчета жесткости стыков сборочных элементов осевых инструментов с учетом направленного формирования параметров поверхностных элементов позволяет аналитически моделировать поверхности контакта геометрических прототипов инструментов различного целевого назначения.

Рецензенты:

Макаров В.М., д.т.н., профессор кафедры «Станки», ФГБОУ ВПО МГТУ «Станкин», г. Москва;

Мнацаканян В.У., д.т.н., профессор, заведующая кафедрой «Технология текстильного машиностроения и конструкционных материалов», ФГБОУ ВПО «Московский государственный университет дизайна и технологии», г. Москва.

Работа поступила в редакцию 08.11.2013.