ISSN 1995-2732 (Print), 2412-9003 (Online) УДК 621.92
DOI: 10.18503/1995-2732-2022-20-1-95-104
МОДЕЛИРОВАНИЕ ИЗМЕНЕНИЯ ГЛУБИНЫ РЕЗАНИЯ В РАЗНЫХ СЕЧЕНИЯХ ОБРАБАТЫВАЕМОЙ ПОВЕРХНОСТИ В ЦИФРОВОМ ДВОЙНИКЕ ОПЕРАЦИИ КРУГЛОГО ВРЕЗНОГО ШЛИФОВАНИЯ С ЧПУ
Акинцева А.В., Переверзев П.П.
Южно-Уральский государственный университет, Челябинск, Россия
Аннотация. Прогнозирование точности на операциях круглого шлифования с ЧПУ является сложной и до сих пор полностью не решенной задачей. Актуальность задачи обуславливается отсутствием у различных производителей цифрового инструмента, позволяющего контролировать управляющие программы на возможность обеспечения заданной точности обработки. Для решения данной задачи предложено использовать цифровой двойник операции круглого шлифования, который описывает процесс формообразования поверхности. Целью работы является разработка модели расчета глубины резания, которая позволит рассчитывать глубину резания в разных сечениях обрабатываемой поверхности при круглом врезном шлифовании вала с переменной податливостью по всей длине шлифования. Модель расчета глубины резания в разных сечениях представляет собой «цифровой двойник» процесса шлифования, который симулирует процесс съема слоя металла переменной толщины одновременно по всей длине шлифуемой поверхности на каждом обороте заготовки с учетом переменной податливости технологической системы. В результате становится возможным расчет текущих значений размеров шлифуемой поверхности, по которым можно рассчитать погрешность обработки диаметральных размеров, формы и взаимного расположения поверхностей. Для разработки модели использованы методы имитационного моделирования, аналитические методы математического моделирования, базирующиеся на фундаментальных физических законах и основных положениях механики резания и теории пластической деформации металла в зоне стружкообразования. Научная новизна приведенных в данной статье исследований заключается в разработке цифрового двойника формообразования обрабатываемой поверхности, который впервые устанавливает взаимосвязь между глубиной резания, режимами резания, упругими деформациями, силами резания и т.д. с точностью обработки. Представленную в данной статье модель расчета глубины резания возможно использовать не только для прогнозирования точности обработки в системе контроля управляющей программы ЧПУ, но и для оптимизации режимов резания проектируемой операции. Это обуславливает большую практическую значимость полученных результатов в производстве.
Ключевые слова: круглое шлифование, цифровой двойник, модель глубины резания, жесткость технологической системы, погрешность обработки.
Исследование выполнено при финансовой поддержке Министерства науки и высшего образования Российской Федерации (грант № ЕЕЫи-2020-0020).
Для цитирования
Акинцева А.В., Переверзев П.П. Моделирование изменения глубины резания в разных сечениях обрабатываемой поверхности в цифровом двойнике операции круглого врезного шлифования с ЧПУ // Вестник Магнитогорского государственного технического университета им. Г.И. Носова. 2022. Т. 20. №1. С. 95-104. https://doi.org/10.18503/1995-2732-2022-20-1-95-104
© Акинцева А.В., Переверзев П.П., 2022
Контент доступен под лицензией Creative Commons Attribution 4.0 License. The content is available under Creative Commons Attribution 4.0 License.
НОВЫЕ ТЕХНОЛОГИЧЕСКИЕ ПРОЦЕССЫ И ОБОРУДОВАНИЕ
MODELING THE CUTTING DEPTH CHANGE IN DIFFERENT SECTIONS OF THE PROCESSED SURFACE IN THE DIGITAL TWIN OF THE CIRCULAR PLUNGE CNC GRINDING OPERATION
Akintseva A.V., Pereverzev P.P.
South Ural State University, Chelyabinsk, Russia
Abstract. Predicting accuracy on circular CNC grinding operations is a challenging task, which has still not been completely solved. Relevance of the task is determined by the lack of a digital tool at various manufacturers to monitor control programs for ensuring the set processing accuracy. To accomplish this task, it is proposed to use a digital twin of the circular grinding operation, which describes the process of surface shaping. The objective of the research is to develop a model for calculating the cutting depth in different sections of the processed surface at a circular plunge grinding of the shaft with variable ease along the entire grinding length. The model for calculating the cutting depth in different sections is a «digital twin» of the grinding process, which simulates the process of removing a metal layer of variable thickness simultaneously along the entire length of the grinded surface on each revolution of the workpiece, considering the variable ease of the technological system. As a result, it becomes possible to calculate the current values of the grinded surface dimensions used to calculate the processing error of diametrical dimensions, shape and relative location of surfaces. To develop the model, the authors used simulation modeling methods with applying analytical methods of mathematical modeling based on the fundamental physical laws and fundamentals of cutting mechanics and the theory of metal flow in the chip formation zone. Scientific novelty of the studies presented in this article lies in the development of a digital twin of the processed surface shaping, which for the first time establishes the interrelation between the cutting depth, cutting modes, elastic deformations, cutting forces, etc. and processing accuracy. The model for calculating the cutting depth, presented in this paper, can be used not only to predict processing accuracy in the monitoring system of the CNC control program, but also to optimize the cutting modes of the designed operation. This determines great practical relevance of the obtained results in manufacturing.
Keywords: circular grinding, digital twin, cutting depth model, rigidity of the technological system, processing error.
The research was funded by the Ministry of Science and Higher Education of the Russian Federation (grant No. FENU-2020-0020).
For citation
Akintseva A.V., Pereverzev P.P. Modeling the Cutting Depth Change in Different Sections of the Processed Surface in the Digital Twin of the Circular Plunge CNC Grinding Operation. Vestnik Magnitogorskogo Gosudarstvennogo Tekhnicheskogo Universiteta im. G.I. Nosova [Vestnik of Nosov Magnitogorsk State Technical University]. 2022, vol. 20, no. 1, pp. 95-104. https://doi.org/10.18503/1995-2732-2022-20-1-95-104
Введение
Прогнозирование точности обработки на операциях круглого шлифования с ЧПУ, в которых программные подачи изменяются в автоматическом ступенчатом цикле, является сложной и до сих пор полностью не решенной научной задачей. Точность обработки прогнозируется по расчетным размерам радиусов прошлифованной поверхности детали, которая сформирована в результате многократного съема слоя металла на каждом обороте заготовки в течение всего цикла шлифования. Поэтому необходима модель расчета многократного съема слоя металла на каждом обороте заготовки в течение всего цикла шлифования. Цикл шлифования представляет собой ступенчатое изменение режимов резания в
зависимости от оставшейся части припуска. Применение цикла позволяет использовать весь диапазон режимов резания и сочетать при обработке заготовки различные стадии (черновую, получистовую и чистовую).
Известные на данный момент научные исследования в области прогнозирования точности обработки в области шлифования по большей части направлены на моделирование получаемого микрорельефа и волнистости обрабатываемой поверхности - шероховатости [1-8]. В работах [912] рассмотрено прогнозирование точности при шлифовании, но в данных работах нет учета силы резания, упругих деформаций технологической системы и влияния переменных технологических условий (колебание припуска в партии детали, степени затупления круга и т.д.) на точность об-
работки. По сути, происходит прогнозирование точности при обработке в идеальных условиях.
Важнейшим фактором, напрямую влияющим на точность обработки, является переменная податливость технологической системы (ТС) по длине обрабатываемой поверхности. Из-за переменной податливости ТС, под действием силы резания происходит неравномерный прогиб участков вала относительно режущего инструмента, что обуславливает разную глубину резания, а значит, разные значения радиусов и погрешность размеров, формы и расположения всей обработанной поверхности вала. Учитывая ступенчатое изменение программной подачи и нелинейное изменение фактической подачи и силы резания по длине обрабатываемой поверхности из-за переменной податливости ТС, можно сделать вывод, что глубина резания тоже будет изменяться по длине обработки на каждом обороте заготовки.
Для расчета параметров формообразования шлифуемой поверхности и прогнозирования точности обработки необходима модель одновременного расчета переменной глубины резания на всех участках обрабатываемой поверхности на каждом обороте заготовки, в зависимости от податливости каждого участка вала по всей длине шлифования. Совокупность переменных глубин резания по всей длине шлифования образует слой снимаемого металла в каждый момент времени.
Модель расчета глубины резания в разных сечениях представляет собой «цифровой двойник» процесса шлифования, который симулирует процесс съема слоя металла переменной толщины (толщина слоя равна глубине резания) одновременно по всей длине шлифуемой поверхности на каждом обороте заготовки с учетом переменной податливости ТС и переменной глубины резания (равной толщине слоя) по всей длине обрабатываемой поверхности при заданных условиях шлифования.
Цифровой двойник позволяет проводить расчет текущих значений параметров формообразования поверхности, описывающих текущие значения размеров шлифуемой поверхности в течение всего цикла шлифования, и прогнозировать размеры обработанной поверхности вала в конце цикла шлифования путем суммирования толщин слоев съема металла в разных сечениях по длине обрабатываемой поверхности. Результатом такого моделирования являются расчетные значения колебания радиусов обработанной поверхности вала в разных сечениях вала по длине обработки, по ко-
торым можно рассчитать погрешность обработки.
Таким образом, решение задачи прогнозирования точности обработки возможно при помощи цифрового двойника операции круглого врезного шлифования с ЧПУ, созданного на основе системного моделирования множества взаимосвязанных широкодиапазонных аналитических моделей процессов в единую модель расчета значений фактических подач, силы резания и точности обработки с учетом основных технологических факторов.
Полученные результаты и их обсуждение
Рассмотрим взаимосвязь глубины резания Л^ в разных сечениях (условно обозначим С1,
С2, С3, С4, С5) обрабатываемой поверхности вала с упругими деформациями ус п в каждом
сечении 5 и на каждом обороте заготовки п. Для начала рассмотрим модель процесса шлифования на первом обороте заготовки (рис. 1, п = 1).
На рис. 1 показан гибридный график, совмещающий графики подач ^ , Лп2 , Лф ,
упругих перемещений ус п и на первом обороте
п =1 заготовки радиусов заготовки Яс,п в поперечных сечениях, равномерно распределенных по длине обрабатываемой поверхности В. При построении расчетной схемы ось заготовки принята неподвижной, а упругие перемещения технологической системы (ТС) относительно оси заготовки совершает шлифовальный круг. Поэтому режущая кромка и ось шлифовального круга не параллельны оси заготовки, а из-за переменной податливости ТС по длине вала расположены под некоторым углом.
На рис. 1 по оси ординат отложены значения программной и фактической подач, а по оси абсцисс - длина шлифуемой поверхности. За время первого оборота заготовки при п = 1 режущая кромка круга подается на величину программной радиальной подачи за оборот Лп , мм/об, на г-й
ступени цикла при г = 1. Но из-за упругой деформации ТС в процессе шлифования режущая кромка круга сдвигается на графике вниз на величину упругой деформации ус п, мм. Поскольку податливость ТС по длине обрабатываемой поверхности имеет переменные значения, то ось заготовки относительно режущей кромки круга сдвигается неравномерно, в зависимости от величины податливости в поперечных сечениях вала.
НОВЫЕ ТЕХНОЛОГИЧЕСКИЕ ПРОЦЕССЫ И ОБОРУДОВАНИЕ
^Ь^фова'льныйкпуг
Рис. 1. Расчетная схема взаимосвязи программной и фактической подач с упругими деформациями и радиусами в разных сечениях обрабатываемой поверхности на первом, втором и третьем оборотах заготовки: 1 - траектория радиусов по длине обрабатываемой поверхности на n-м обороте заготовки Fig. 1. Calculation model of the interrelation between program and actual feeds and elastic deformations and radii in different sections of the processed surface on the first, second and third revolutions of the workpiece: 1 is a trajectory of radii along the length of the processed surface on the n-th revolution of the workpiece
В сечении правого торца податливость больше, чем в сечении левого торца, и упругие деформации в сечении правого торца больше, чем в левом. Поэтому режущая кромка круга показана под углом как результат неодинаковых деформаций ТС. В результате глубина резания, то есть фактическая скорость подачи А/а ,
мм/об, будет меньше, чем программная А/п^ .
На графике на рис. 1 показаны величины А/а и ус=2 п=1 для сечения С2 и п = 1. Из
фс= 2 ,п=2 '
графика видно, что сумма фактической подачи и упругой деформации усп равны программной подаче А/п2 за оборот заготовки, то
есть уравнение баланса подач и упругих деформаций за оборот заготовки можно записать в следующем виде:
„ Фс ,и + Уе,п = Ус ,п = Ус Pyn ,
(1) (2)
где - программная радиальная подача,
мм/об; А/а - фактическая радиальная подача, мм/об; усп - упругая деформация технологической системы, мм; Р - радиальная составляю-
Уп
щая силы резания на п-м обороте заготовки; ус -податливость ТС в заданном сечении вала.
Отметим, что переменные А/^ и А/ф являются скоростями подач за один оборот заготовки.
Показанная на графике подача /П , мм - это проПп
граммная подача (накопленная) за несколько п оборотов заготовки, расчитывается по формуле
Акинцева А.В., Переверзев П.П.
г N
П =ИЛпг
г=1п=1
(3)
где г - порядковый номер ступени цикла; Z - максимальное количество ступеней в цикле; п - порядковый номер оборота заготовки на г-й ступени цикла; N - максимальное количество оборотов заготовки на г-й ступеней цикла.
С учетом уравнения (3) уравнение баланса технологических перемещений для первого оборота заготовки при п = 1 запишем следующим образом (на рис. 1 размерная цепочка (баланс) выделена утолщённой линией):
'П
= Л/
П
2=1 = Лфс,п=1 + Ус ,п=1 . (4)
Из уравнения (4) выразим искомую глубину резания Лф 1 в сечении С при п = 1:
Л'
фс
= /
П„
ус,п=1 ■
(5)
Текущее значение радиуса заготовки Дс п=1
можно рассчитать как разницу между исходным радиусом заготовки Дзаг и глубиной резания
Л'
фс
(см. рис. 1), то есть
д
с, п=1
= дза
-Л',
фс
(6)
где Дс п=1 - текущее значение радиуса заготовки.
Рассмотрим взаимосвязь глубины резания Лф в разных сечениях обрабатываемой поверхности вала с упругими деформациями ус п=2
на втором обороте заготовки при п = 2. Уравнение баланса технологических перемещений для второго оборота заготовки при п = 2 запишем следующим образом (на рис. 1 размерная цепочка (баланс) выделена утолщённой прерывистой линией):
'п - =ЛП _ + Л/
п=2
п
2=1
'п
г=1
= Л',
фс ,п=2 + Ус ,п=2 + 'фс ,п =2 '
(7)
где
фс
- часть припуска (на радиус), снятого
за предыдущие (п-1) обороты заготовки, мм,
N
'фс=2,п = Х 'фс=2,п .
п=1
(8)
На первом обороте заготовки при п = 1 значение ^ =0, так как предыдущего оборота при
ф с, п=1
п = 0 не было и глубина резания на предыдущем обороте равна нулю А^ = 0 . На втором обо-
фс, п=0
роте заготовки при п = 2, величина ^ равна
фс, п=2
глубине резания л^ на предыдущем обороте
фс, п= 1
при п = 1, то есть ' = ЛЛ .
Г фс,п=2 фс,я=1
Из уравнения (7) выразим глубину резания Л'а в сечении С при п = 2 (при / - =
фс, п=2 фс ,п=2
= Л'ф ,):
тс, п=1
Л/л
фс,п=2 = 'Пп=2 - ус,п=2 - 'фс,п =2 .
(9)
Текущее значение радиуса Дс п=2 заготовки
можно рассчитать как разницу между исходным радиусом заготовки Дзаг и суммой глубин
резания на первых двух оборотах заготовки и Л^ (см. рис. 1), то есть
фс, п=1 фс, п=2 Г
д
с, п =2 Дзаг„
- Лф - Лф
фс, п=1 фс ,п=2
—д,
/фс,п= 2 Л'фс,п= 2 '
(10)
Величина снятой части припуска /- ,
г ф-сум ^п=2
мм, за два оборота заготовки рассчитывается как сумма глубин резания на всех оборотах
фс ,п
заготовки по формуле п=2
/ф сумс,п=2 = Х 'фс,п = Л'фс,п=1 - 'фс,п=2 . (11) 1
Соответственно, на первом обороте заготовки
'ф.сумс,п=1 = Л'фс,п=1 ' (12)
Текущее значение радиуса Дс п заготовки
можно рассчитать как разницу между исходным радиусом заготовки Дзаг и суммой глубин
резания (см. рис. 1), то есть
Дс,п Дзагтах /ф
тах ф.суМс,п=2 '
(13)
Аналогично получим уравнения взаимосвязи подач, упругих деформаций и радиусов заготовки в процессе круглого врезного шлифования за три оборота заготовки при п = 3. Накопленная программная подача равна трем программным подачам за оборот, то есть уравнение баланса подач и упругих деформаций можно записать в следующем виде (на рис. 1 размерная цепочка (баланс) выделена сплошной линией):
НОВЫЕ ТЕХНОЛОГИЧЕСКИЕ ПРОЦЕССЫ И ОБОРУДОВАНИЕ
/Пп=з = 3 А/П = А/фс п=з + Ус,п=3 + /фСп=з , (14)
/фс,п=3 = А/фС,п=1 +А/фС,п=2
(15)
Из уравнения (14) получим значение глуби-
ны резания:
'фс
= /
п.
ус,п=3 - /фс_п=з ■
Величина снятой части припуска /
ф.сумс
(16)
за
три оборота заготовки рассчитывается как сумма предыдущих глубин резания А/ на всех обо-
фс ,п
ротах заготовки по формуле
п=3
'ф.сум„=з = 1 А'ф„ = 1
= А/ф + А/ф + А/ф .
фс,п=1 фс,п=2 фс,п=3
(17)
Текущее значение радиуса Яс,п=з заготовки можно рассчитать как разницу между исходным радиусом заготовки Язаг и суммой глубин
резания / (см. рис. 1), то есть
КС,П = ^заттах ~ /ф.суМс,п=3 '
(18)
Таким образом, формулу глубины резания при круглом врезном шлифовании вала с переменной податливостью в заданном сечении и обороте заготовки п можно записать в следую-
щем виде:
А/
фс,п = 'П - ус,п - 'ф,
(19)
в которой /ф = 0 при п = 1.
Раскроем значение ус п путем подстановки формулы (2), получим
А/фс,п = 'П ~ГсРУп - /фс,п ■ (20)
Формула (20) представляет собой математическую модель расчета текущего значения глубины резания (фактической скорости подачи) А/ на каждом обороте п заготовки, при из-
фс ,п
вестном значении значение общей радиальной силы резания Ру и величин податливости уС в
уп
каждом сечении С вала.
Однако расчет по формуле (20) глубин резания А/ в разных сечениях вала невозможно,
фс ,п
так как неизвестно значение общей радиальной силы резания РУ , которая представляет собой
уп
сумму сил резания на всех участках вала по длине обрабатываемой поверхности.
Для установления взаимосвязи общей радиальной силы резания РУ с податливостью ТС
уп
разобьем обрабатываемую поверхность вала на элементарные участки длиной В^. На рис. 2 показана обрабатываемая поверхность вала, разбитая на четыре участка, расположенных между пятью сечениями С1-С5. Соответственно, обозначим номер каждого участка как 7 = 1, 7 = 2, 7 = 3 и 7 = 4.
При известных значениях податливости уС сечений вала можно определить податливость у7 каждого элементарного участка вала как среднее арифметическое податливости сечений, между которыми расположен участок. Тогда упругую деформацию у^п участка вала на п-м обороте
заготовки в процессе шлифования можно записать как
Угп = У1ру
Уп
(21)
где у7 - податливость 7-го элементарного участка вала, мм.
Соответственно, глубину резания на каждом участке обозначим через А/ф , а величину сня-
фг ,п
той части припуска за предыдущие обороты заготовки обозначим через / .
фг ,п
Тогда баланс технологических перемещений для каждого участка с учетом его ширины В7 (см. рис. 2) можно записать в следующем виде:
'ПпВг = (А/ф,и + Уг,п + 'ф, ■ (22)
Учитывая, что участки вала не могут перемещаться независимо друг от друга, баланс совместных технологических перемещений запишем как равенство площадей перемещений по длине обрабатываемой поверхности:
ПВ = В I (А/ф,п +Угруп + 'ф,п ) , (23)
I=1
где 7 - порядковый номер участка на обрабатываемой поверхности (I - максимальное количество участков на обрабатываемой поверхности вала).
Акинцева А.В., Переверзев П.П.
Рис. 2. Расчетная схема взаимосвязи программной и фактической подач с упругими деформациями и радиусами
в разных сечениях обрабатываемой поверхности при переменной податливости ТС по длине обрабатываемой поверхности: 1 - траектория радиусов по длине В обрабатываемой поверхности на последнем обороте заготовки; 2 - траектория радиусов по длине В обрабатываемой поверхности на предыдущем обороте заготовки Fig. 2. Calculation model of the interrelation between program and actual feeds and elastic deformations and radii in different sections of the processed surface at variable ease of the technological system along the length of the processed surface: 1 is a trajectory of radii along length B of the processed surface on the last revolution of the workpiece; 2 is a trajectory of radii along length B of the processed surface on the previous revolution of the workpiece
Количество участков вала можно рассчитать по формуле
1Пп -уш + ^1фшп + гфт1
(27)
I - — — .
(24)
После подстановки формулы (24) в уравнение (23) и проведения соответствующих преобразований получим
П„
! I
- у +УРп + Ч, 1 i-i
(25)
где ут - средняя податливость ТС на длине шлифования В, мм; Л'фт - средняя глубина резания на длине шлифования В, мм; /фт - среднее значение снятой части припуска за предыдущие обороты заготовки на длине шлифования В.
Учитывая, что среднее значение упругой деформации у равно
(П„ - Ру
i-1
- +
I I
у I у I
I !%,n +1 i-i i-i
(26)
Утп утРу„ , уравнение (27) запишем в виде
*Пп - Уш, + А(фШп + ^фшп
(28)
(29)
Поделим на количество участков I все слагаемые в уравнении (26), получим
Таким образом, формулы (27)-(29) устанавливают баланс технологических перемещений при
п
НОВЫЕ ТЕХНОЛОГИЧЕСКИЕ ПРОЦЕССЫ И ОБОРУДОВАНИЕ
общей силе резания Ру и средней податливости
уп
ут ТС, включая среднюю упругую деформацию ут и среднюю глубину резания Ыфт на каждом
обороте заготовки в цикле шлифования.
Для нахождения суммарной силы резания Руп в цикле шлифования воспользуемся полученной ранее и доказанной экспериментально формулой [13]:
Pyn = К3Афтп + К4ЛЩт,
K3 = Trdßn
V
K4 ==nß\[^ 4 3 V D + d
(30)
(31)
(32)
где Oi - интенсивность напряжении в движущемся объеме металла, характеризующая сопротивление металла пластической деформации при степени, скорости деформации и температуре, присущих процессу шлифования; sci - интенсивность скорости деформации объема металла в зоне сдвига, 1/с; п - степень затупления шлифовального круга, равная относительной опорной поверхности круга по площадкам затупления зерен; d - диаметр обрабатываемой поверхности заготовки, мм; D - диаметр заготовки, мм; В -длина обрабатываемой поверхности заготовки, мм; n - частота вращения заготовки, мм/мин;
At,
фтп
- скорость фактической подачи в задан-
ном сечении вала на п-м обороте заготовки, мм/об; Ук - скорость вращения круга, м/с.
Фактическая накопленная подача за предыдущие обороты заготовки Хфт тоже одинакова
во всех сечениях вала и рассчитывается на п-м обороте заготовки по формуле (33) как сумма глубин резания за предыдущие обороты заготовки или как часть припуска (на радиус), снятого за предыдущие (п-1) обороты заготовки:
N
Vn=I А^п
п=1
(33)
Решим уравнение (29) относительно At^m совместно с (30)-(33), получим:
At,
фтп
~ГтК4 2(1 + КзГт )
- +
V
УтК4
2(1 + КзГт )
tn " f<
фтп
1 + К3Ут
(34)
Формула (34) представляет собой математическую модель расчета скорости средней фактической подачи (глубины резания) Atфm в течение всего автоматического цикла круглого врезного шлифования с ЧПУ на каждом п-м обороте заготовки при средней податливости ТС, в зависимости от основных технологических факторов.
Расчет текущих значений глубины резания
At,
фтп
ведется в течение всего цикла шлифования
в разных поперечных сечениях вала по длине обрабатываемой поверхности с различной податливостью и при разных углах поворота радиусов профиля заготовки, имеющей исходное радиальное биение. Текущие значения радиусов обрабатываемой поверхности рассчитываются путем уменьшения исходного радиуса заготовки на величину глубины резания. Погрешность обработки в конце цикла шлифования рассчитывается по разбросу значений радиусов обработанной поверхности в разных сечениях и углах поворота радиусов обрабатываемой поверхности. Таким образом, решена задача моделирования формообразования технологического размера и его погрешности в цифровом двойнике операции круглого врезного шлифования с ЧПУ. Подобный подход является типовым при разработке цифрового двойника для других видов шлифования.
Заключение
1. Прогнозирование точности на операциях круглого шлифования с ЧПУ является сложной и до сих пор нерешенной задачей. Сложность данной задачи обусловлена наличием множества взаимосвязей между параметрами процесса резания, вызывающего теплофизические изменения, упругие деформации технологической системы и т.д. Важнейшим фактором, напрямую влияющим на точность обработки, является переменная податливость технологической системы по длине обрабатываемой поверхности. Из-за переменной податливости ТС под действием силы резания происходит неравномерный прогиб участков вала относительно шлифовального круга. Это обуславливает разную глубину резания, а
2
значит, разные значения радиусов и погрешность размеров, формы и расположения всей обработанной поверхности вала.
2. Для решения задачи прогнозирования точности обработки разработан «цифровой двойник» процесса шлифования, который симулирует процесс съема слоя металла переменной толщины одновременно по всей длине шлифуемой поверхности на каждом обороте заготовки с учетом переменной податливости ТС. Основу цифрового двойника составляет аналитическая модель расчета глубины резания, которая позволяет рассчитывать глубину резания в разных сечениях обрабатываемой поверхности при круглом врезном шлифовании вала с переменной податливостью по всей длине шлифования. Разработанная методика расчета текущих значений размеров шлифуемой поверхности за время съема припуска и ступенчатом изменении подачи позволяет прогнозировать погрешность обработки диаметральных размеров вала с переменной податливостью, отклонения формы и взаимного расположения поверхностей.
3. Прогнозирование точности обработки на основе модели расчета глубины в разных сечениях обрабатываемой поверхности в цифровом двойнике можно использовать для контроля вновь спроектированных управляющих программ для станков с ЧПУ [14] и в системе оптимизации режимов резания проектируемой операции [15].
Список литературы
1. Gong Y.D., Wang B., Wang W.S. The simulation of grinding wheels and ground surface roughness based on virtual reality technology // Journal of Materials Processing Technology. 2002. Vol. 129. Pp. 123-126.
2. Hecker R., Liang S.Y. Predictive modeling of surface roughness in grinding // I. J. of Machine Tools and Manufacture. 2003. Vol. 43(8). Pp. 755-761.
3. Zhou X., Xi F.Modeling and predicting surface roughness of the grinding process // International Journal of Machine Tools & Manufacture. 2002. Vol. 42. Pp. 969-977.
4. Kumar S., Paul S. Numerical modelling of ground surface topography: effect of traverse and helical superabra-sive grinding with touch dressing // Production Engineering. 2012. Vol. 6 (2). Pp. 199-204.
5. Широков А.В., Осипов А.П. К вопросу о прогнозировании и обеспечении параметров шероховатости шлифованной поверхности // Известия высших учебных заведений. Машиностроение. 2007. №6. C. 76-88.
6. Воронов С.А., Вэйдун Ма, Воронова И.С. Стохастическая модель процесса абразивной обработки. Кинематика плоского шлифования // Известия высших учебных заведений. Машиностроение. 2017. №11. C. 68-78.
7. Козлов А.М., Ефремов В.В. Формирование микрорельефа при обработке абразивным инструментом // Известия высших учебных заведений. Машиностроение. 2004. №1. C. 59-64.
8. Новоселов Ю.К., Богуцкий В.Б., Дзюбаба Р.Н. Прогнозирование параметров шероховатости поверхности при абразивной обработке // Известия Тульского государственного университета. Технические науки. 2017. №8. C. 262-269.
9. Voronov S.A., Veidun M. Mathematical Modeling of the Cylindrical Grinding Process // Journal of Machinery Manufacture and Reliability. 2017. Vol. 46(44). Pp. 394-403.
10. SolerYa.I., Nguyen Van Le Modeling and predicting flatness deviations under plane grinding of parts made of 30HGSA hardened steel // Вестник Иркутского государственного технического университета. 2016. № 4 (111). Pp. 54-65.
11. Колтунов И.И., Лобанов А.С. Зависимость точности шлифования от параметров наладки // Известия Московского государственного технического университета МАМИ. 2012. №1(13). C. 223-226.
12. Osipov A.A., Osipov A.P. Theoretical research of work-piece surface forming in grinding // Materials Today: Proceeding. 2020. Vol. 38. Pp. 1-3.
13. Акинцева А.В. Основные этапы моделирования послойного съема металла с применением цифрового двойника операции круглого шлифования с ЧПУ // Вестник Южно-Уральского государственного университета. Серия: Машиностроение. 2020. №3. C. 48-58.
14. Переверзев П.П. Модель формирования погрешностей внутришлифовальной обработки // СТИН. 2016. № 6. C. 25-30.
15. Akintseva A.V., Pereverzev P.P. Complex optimization of parameters for controlling the cycle of internal grinding by the method of dynamic programming // J. MATEC Web of Conferences. 2017. Vol. 129. No. 01019.
References
1. Gong Y.D., Wang B., Wang W.S. The simulation of grinding wheels and ground surface roughness based on virtual reality technology. Journal of Materials Processing Technology, 2002, vol. 129, pp. 123-126.
2. Hecker R., Liang S.Y. Predictive modeling of surface roughness in grinding. I. J. of Machine Tools and Manufacture, 2003, vol. 43(8), pp. 755-761.
3. Zhou X., Xi F. Modeling and predicting surface roughness of the grinding process. International Journal of
НОВЫЕ ТЕХНОЛОГИЧЕСКИЕ ПРОЦЕССЫ И ОБОРУДОВАНИЕ
Machine Tools & Manufacture, 2002, vol. 42, pp. 969-977.
4. Kumar S., Paul S. Numerical modelling of ground surface topography: effect of traverse and helical superabra-sive grinding with touch dressing. Production Engineering, 2012, vol. 6 (2), pp. 199-204.
5. Shirokov A.V., Osipov A.P. On the issue of predicting and ensuring the roughness parameters of the ground surface. Izvestiya vysshikh uchebnykh zavedeniy, Mashi-nostroenie [Proceedings of Higher Educational Institutions. Machine Building], 2007, no. 6, pp. 76-88. (In Russ.)
6. Voronov S.A., Weidong Ma, Voronova I.S. A stochastic model of abrasive processing. Kinematics of flat grinding. Izvestiya vysshikh uchebnykh zavedeniy, Mashi-nostroenie [Proceedings of Higher Educational Institutions. Machine Building], 2017, no. 11, pp. 68-78. (In Russ.)
7. Kozlov A.M., Efremov V.V. Formation of microrelief during processing by an abrasive tool. Izvestiya vysshikh uchebnykh zavedeniy, Mashinostroenie [Proceedings of Higher Educational Institutions. Machine Building], 2004, no. 1, pp. 59-64. (In Russ.)
8. Novoselov Yu.K., Bogutsky V.B., Dzyubaba R.N. Prediction of surface roughness parameters during abrasive processing. Izvestiya Tulskogo gosudarstvennogo uni-versiteta. Tekhnicheskie nauki [Proceedings of Tula State University. Engineering Sciences], 2017, no. 8, pp. 262-269. (In Russ.)
9. Voronov S.A., Veidun M. Mathematical modeling of the cylindrical grinding process. Journal of Machinery Manufacture and Reliability, 2017, vol. 46(44), pp. 394-403.
10. Soler Ya.I., Nguyen Van Le. Modeling and predicting flatness deviations under plane grinding of parts made of 30KhGSA hardened steel. Vestnik Irkutskogo gosudar-stvennogo tekhnicheskogo universiteta [Proceedings of Irkutsk State Technical University], 2016, no. 4 (111), pp. 54-65. (In Russ.)
11. Koltunov I.I., Lobanov A.S. Dependence between grinding accuracy and adjustment parameters. Izvestiya Mos-kovskogo gosudarstvennogo tekhnicheskogo universiteta MAMI [Proceedings of Moscow State Technical University MAMI], 2012, no. 1(13), pp. 223-226. (In Russ.)
12. Osipov A.A., Osipov A.P. Theoretical research of work-piece surface forming in grinding. Materials Today: Proceedings, 2020, vol. 38, pp. 1-3.
13. Akintseva A.V. Main stages of modeling a layer-by-layer metal removal, using a digital twin of the circular CNC grinding operation. Vestnik Yuzhno-Uralskogo gosudarstvennogo universiteta. Seriya: Mashinostroenie [Bulletin of SUSU. Series: Mechanical Engineering Industry], 2020, no. 3, pp. 48-58. (In Russ.)
14. Pereverzev P.P. A model of generating errors of internal grinding processing. STIN [Machines and Tools], 2016, no. 6, pp. 25-30. (In Russ.)
15. Akintseva A.V., Pereverzev P.P. Complex optimization of parameters for controlling the cycle of internal grinding by the method of dynamic programming. J. MATEC Web of Conferences, 2017, vol. 129, no. 01019.
Поступила 08.02.2022; принята к публикации 22.02.2022; опубликована 25.03.2022 Submitted 08/02/2022; revised 22/02/2022; published 25/03/2022
Акинцева Александра Викторовна - кандидат технических наук, доцент кафедры техники, технологии и строительства, Южно-Уральский государственный университет, Челябинск, Россия. Email: akintsevaav@susu.ru. ORCID 0000-0003-0588-1401 Переверзев Павел Петрович - доктор технических наук, доцент, профессор кафедры технологии автоматизированного машиностроения, Южно-Уральский государственный университет, Челябинск, Россия. Email: pereverzevpp@susu.ru. ORCID 0000-0002-1416-8918
Aleksandra V. Akintseva - PhD (Eng.), Associate Professor of the Department of Engineering, Technology and Construction, South Ural State University, Chelyabinsk, Russia. Email: akintsevaav@susu.ru. ORCID 0000-0003-0588-1401
Pavel P. Pereverzev - DrSc (Eng.), Associate Professor
of the Department of Automated Mechanical Engineering Technology, South Ural State University, Chelyabinsk, Russia. Email: pereverzevpp@susu.ru. ORCID 0000-0002-1416-8918