employee from 01.01.1997 until now
Orenburg, Orenburg, Russian Federation
employee
Orenburg, Orenburg, Russian Federation
graduate student from 01.01.1920 until now
Russian Federation
UDK 62 Инженерное дело. Техника в целом. Транспорт
OKSO 15.03.05 Конструкторско-технологическое обеспечение машиностроительных производств
BBK 346 Отдельные машиностроительные и металлоперерабатывающие процессы и производства
The paper presents the result of modeling the thermal deformation behavior of the bearing system of a double-sided face grinder. Two modes of its operation is considered: idling and working strokes. Spindle bearings, heat fluxes from coolant, and heat fluxes from motors are considered as the main heat sources. The preliminary values of heat and convective fluxes are assigned according to well-known methods used in engineering calculations of machine tools. Refinement of the values of heat and convective fluxes is carried out according to the results of experimental values of temperatures and temperature displacements. In the thermal model, convective heat transfer is assigned to 450 surfaces. Heat flows are assigned to 161 surfaces. The simulation is run in Ansys for 6 hours of machine tool operation. Ansys uses a ten-node Solid227 element as a typical finite element. The problem of thermoelasticity is solved in a related setting. The error in modeling thermal processes did not exceed 0,5 degree, the error in modeling temperature displacements did not exceed 5 microns. The results of computer simulation confirmed the relative experimental positions of the left and right grinding wheels characteristic of this range of machines. At idle, the relative position of the grinding wheels represented the state «wider at the bottom», and during the working stroke, the state «at the bottom narrower» is recorded. The illustrations of the FEM model and contour plots for the temperature field and temperature displacements of the bearing system of the machine tool for two modes of its operation are presented.
face grinding machines, temperature deformations, heat resistance, accuracy improvement, techniques
Введение. Системные исследования влияния тепловых процессов на точность станков были начаты ещё в 50 годах. Однако, сегодня, несмотря на всеобщую цифровизацию экономики, по-прежнему остаются актуальными тепловые проблемы, возникающие в металлорежущих станках [1–7]. Это объясняется ростом скоростей резания и подач [8, 9]. Сегодня, благодаря усилиям отечественных и зарубежных исследователей, сформирована общая концепция исследований тепловых процессов в металлорежущих станках [1–7]. Исследованы основные источники тепловыделения в станках [10], тепловое поведение отдельных узлов станка, влияние температурных деформаций обрабатываемых деталей на их выходную точность [11]. Экспериментально определены коэффициенты теплоотдачи для деталей станка, коэффициенты теплопроводности стыков. Разработаны различные методики инженерных расчётов и компьютерного моделирования [1, 12–14].
Сегодня с более широким внедрением систем компьютерного моделирования в область тепловых исследований станков приобретают особую значимость интеллектуальные технологии компенсации температурной погрешности станков [15–23].
Рынок современных систем инженерного анализа достаточно объёмен и превышает 7 млрд. долларов США с положительной тенденцией роста более 20 % в год, начиная с 2016 года [24]. Среди наиболее известных на мировом рынке вэндоров CAE-систем: Ansys, Dassault Systems, Siemens PLM Sowtware, MSC Software, Altair Engineering, Cybernet Systems, ESI Group, Autodesk, CD-Adapco, Comsol, MAYA, PTC [25].
Несмотря на то, что сегодня вычислительная производительность компьютеров многократно превышает производительность компьютеров, на которых выполнялось моделирование первых отечественных систем инженерного анализа, реальная производительность компьютеров, используемых в научных исследованиях станков в стране и мире, а также упрощённые трибологические модели тепловых источников, ограничивают возможности реализации подробных расчётных моделей станков. Поэтому по-прежнему при разработке расчётных моделей станков актуально применение различных упрощений, главными из которых являются следующие: в качестве основных тепловых источников принимаются подшипники, тепловыделение в которых не зависит от времени. Теплоотдача учитывается в виде конвективного теплообмена. Коэффициенты теплоотдачи также принимаются постоянными. Влияние стыков в сопряжениях не учитывается.
При исследовании двусторонних торцешлифовальных станков в работе в качестве основного инструмента моделирования термодеформационного поведения несущей системы станка (расчёт нестационарных температурных полей и температурных деформаций) использован Ansys [26]. Этот инструмент широко известен в мире и России, и его позиции ещё больше укрепляются за счёт создания нового продукта Ansys Twin Builder, который используется для набирающих силу технологий создания цифровых двойников [27].
Общие положения расчета температурных полей. Анализ несущей системы торцешлифовального станка позволяет представить её деталями двух типов: коробчатыми и телами вращения.
К первому типу относятся корпуса шпиндельных бабок и станины. Ко второму типу относятся шпиндельные узлы, состоящие из шпинделя и корпуса пиноли, являющиеся телами вращения.
В качестве типового конечного элемента в Ansys использован десятиузловой элемент Solid227.
При проектировании гаммы торцешлифовальных станков их расчётные модели различаются как интенсивностью и месторасположением источников тепловыделений, так и конструктивным исполнением, определяемым способом подачи заготовки в зону обработки. Это условие определяет особенность разрабатываемой универсальной расчетной модели торцешлифовального станка, которая бы учитывала возможные особенности конструктивного исполнения, расположение тепловых источников и поверхностей теплообмена, что определяет степень влияния температурных смещений кругов на точность обработки [28–30].
Разработка расчетной схемы, учитывающей особенности конструктивного исполнения двусторонних торцешлифовальных станков с горизонтальным расположением шпинделей. Конструкцию торцешлифовальных станков данного типа (рис. 1) отличает то, что шпиндельный узел расположен в пиноли, которая обеспечивает ограниченный теплоотвод от шпиндельного узла в корпус шпиндельной бабки. Теплоотвод от шпиндельного узла в пиноль ограничен прокачиваемой смазочно-охлаждающей жидкости (СОЖ) по каналу, выполненному внутри шпинделя. Шпиндельная бабка состоит из двух частей, на нижней части которой установлена пиноль. Шпиндельные бабки (нижней частью) соответственно контактируют со станиной и имеют симметричное расположение относительно оси симметрии, проходящей через центр станины. Теплоотвод от шпиндельных бабок в станину ограничен вследствие конструктивного исполнения корпусов шпиндельных бабок и станины. Ограждение зоны шлифования (на рис. 1 не показано) контактирует со станиной и шпиндельными бабками с ограниченным теплоотводом за счет полимерных прокладок.
Рис. 1. Модель станка:
1 – станина; 2 – бабка нижняя; 3 – бабка верхняя; 4 – крышка верхняя; 5 – крышка торцовая; 6 – шкив;
7 – шпиндельный узел; 8 – тумба; 9 – редуктор; 10 – защитное устройство от СОЖ;
11 – фланец шпиндельного узла, к которому крепится шлифовальный круг; 12 – диск подачи
Температурные деформации стенок корпусов шпиндельных бабок и станины определяют смещение шпинделя в вертикальной плоскости.
Анализ конструкции, условий теплообмена, условий закрепления узлов, позволил представить термодинамическую систему торцешлифовального станка данного типа в виде расчетной модели, состоящей из следующих групп узлов:
– шпиндельный узел, включающий шпиндель, подшипники и пиноль;
– корпус шпиндельной бабки, состоящий из двух частей и направляющие;
– корпус станины с тумбой, редуктором и диском подачи.
Схема назначения граничных условий станка представлена на рис. 2 для центральной части станины и левой половины станка (для правой половины станка условия идентичны).
Рис. 2. Модель станка с граничными условиями (сечение, вид сзади):
Qпод.п., Qпод.з., Qпод.ш – тепловые потоки от подшипников; Q1=f(Vсож), Q3=f(Vсож) – тепловые потоки от СОЖ; Q2=f(n,δ1) – тепловой поток от электродвигателя и ременной передачи;
δ1 – расстояние до вращающегося шкива; δ – расстояние до вращающегося шпинделя; 1 – передний подшипник; 2 – пиноль; 3 – шпиндель; 4 – задний подшипник; 5 – направляющая; 6 – подшипник шкива
При разработке расчётной модели для определения температурного поля несущей системы станка были сделаны следующие допущения:
– основными источниками тепла в шпиндельном узле считались подшипники качения, представляемые в виде распределенных потоков тепла в местах посадки подшипников;
– мощность тепловыделения для подшипников определялась в ходе решения задачи идентификации тепловой модели, для этого использовались результаты экспериментов [29, 30], при этом распределение значений тепловых потоков от подшипника между внутренними и наружными поверхностями модели подшипника принималось пропорционально этим площадям;
– другие источники тепла (гидросистема, пары трения) также рассматривались как распределенные тепловые потоки; мощность тепловыделения, как и для подшипников, также определялась по экспериментальным данным температур и температурных перемещений;
– шпиндельный узел рассматривался в условиях действия сложного теплообмена при одновременном действии теплопроводности и конвекции;
– конвективный теплообмен для несущей системы реализован явным заданием коэффициентов теплоотдачи с учётом условий теплообмена (свободная или вынужденная конвекция; поверхности, расположенные в кольцевом зазоре или расположенные в большом пространстве и т. д.);
– все геометрические параметры элементов геометрической модели несущей системы станка назначались по рабочим чертежам;
– конструкции корпусов станины и шпиндельной бабки имеют коробчатую форму с системой внутренних перегородок, связывающих переднюю, заднюю и торцовые стенки, и обеспечивающие жесткость конструкции;
– материалы расчётной модели: чугун, конструкционная и легированная сталь, полиуретан;
– общее время моделирования работы станка разбивалось на отдельные интервалы времени, в пределах которых, все параметры модели принимались не зависящими от времени;
– в расчётной модели небольшие отверстия в стенках не учитывались;
– стенки с перепадами по высотам менее 20 мм заменялись пластинками некоторой приведенной высотой при сохранении эквивалентной жесткости.
Разработанная конечно-элементная модель станка содержит более 335000 конечных элементов и более 513000 узлов (рис. 3).
Рис. 3. Конечноэлементная модель станка
Назначение граничных условий. Анализ схемы расположения источников теплоты и их действия в термодинамической системе станка показывает, что исходными данными для расчета температурного поля торцешлифовального станка является количество теплоты, выделяемое процессом шлифования, электродвигателями (электрические и магнитные потери), ременной передачей и подшипниками.
Оценочный расчёт мощности тепловыделения для тепловых источников проводился по известным методикам [1–7].
Коэффициент теплоотдачи , характеризующий интенсивность процесса теплоотдачи, зависит от большого числа факторов. На величину
влияют геометрические размеры тела, скорость и характер движения охлаждающей среды, физические параметры используемого охладителя, режим работы. Граничные условия, характеризующие теплообмен поверхностей несущей системы станка, также назначались в соответствии с известными рекомендациями для металлорежущих станков [1–7].
Коэффициенты теплоотдачи назначались для тех поверхностей, которые наиболее существенно оказывали влияние на формирование температурного поля станка: поверхности шпиндельного узла, поверхности корпусов станины и шпиндельной бабки, обращенные к ременной передачи, а также поверхности станины и шпиндельных бабок, омываемые СОЖ.
Поверхности вращающегося шпинделя, на которых установлены подшипники, участвуют в теплообмене с коэффициентом теплоотдачи, зависящим от частоты вращения:
, (1)
где – безразмерный коэффициент теплоотдачи;
– критерий Рейнольдса для радиального движения воздушных потоков (безразмерная величина); Gr – число Грасгофа (безразмерная величина):
, (2)
где l – коэффициент теплопроводности, Вт/м2×°С; – эквивалентный диаметр, м:
, (3)
где d – расстояние между шпинделем и корпусом (пинолью), м:
, (4)
где n – вязкость охладителя, м2/с; – окружная скорость вращения шпинделя, м/с:
. (5)
где dшп – диаметр шпинделя, м; n – частота вращения шпинделя об/с.
, (6)
где b – температурный коэффициент объемного расширения; – ускорение свободного падения, м/с;
– избыточная температура, °С.
Внутренняя поверхность шпинделя, образующая отверстия для подачи СОЖ, имеет коэффициент теплоотдачи равный:
, (7)
где Pr – число Прандтля для СОЖ; th – гиперболический тангенс; – критерий Рейнольдса для осевого движения воздушных потоков:
,………………………(8)
где – осевая скорость потока охладителя, м/с.
, (9)
где – температуропроводность охладителя, м2/с;
– вязкость охладителя, м2/с.
Теплоотдача поверхностей, расположенных на небольшом расстоянии d от вращающихся деталей, зависит от вентиляционного эффекта, создаваемого вращающимися частями. К таким поверхностям относится внутренняя поверхность пиноли, обдуваемая вращающимся шпинделем.
Теплоотдача наружных поверхностей пиноли зависит от температуры окружающего воздуха в условиях отсутствия вентиляции. Коэффициент теплоотдачи наружной поверхности пиноли принят равным 5 Вт/м2×°С, при температуре окружающей среды для данной поверхности 20 °С.
Торцовые стенки (внешняя поверхность) корпуса шпиндельных бабок участвуют в теплообмене с коэффициентом теплоотдачи, зависящим от вентиляционного эффекта, создаваемого вращающимися шкивами ременной передачи . С другого торца, в центральной части станка, часть поверхности корпусов шпиндельных бабок ограничены защитным устройством от СОЖ. Теплоотдача остальных внешних стенок и внутренних поверхностей зависит от температуры окружающего воздуха в условиях отсутствия вентиляции, то есть со свободным движением воздуха. Принимаем коэффициент теплоотдачи данных поверхностей равным 5 Вт/м2×°С, температура окружающей среды также принята равной 20 °С.
Поверхности поддона, расположенного внутри станины, по которым течет СОЖ, имеют коэффициент теплоотдачи как функция от скорости потока СОЖ, определяемый по формуле (при турбулентном режиме):
. (10)
где Reохл – критерий Рейнольдса; Prохл – число Прандтля для параметров СОЖ с температурой набегающего потока; – число Прандтля для параметров СОЖ, равной температуре стенок поддона; Nu – безразмерный коэффициент теплоотдачи:
. (11)
где – теплопроводность охладителя (СОЖ), Вт/м2×°С; l – длина стенки, обтекаемой (обдуваемой) охладителем, м.
. (12)
где – средняя скорость потока охладителя, м/с.
Теплофизические параметры для СОЖ принимались из каталогов производителей СОЖ. В тепловой модели конвективный теплообмен был назначен для 450 поверхностей. Тепловые потоки (Heat Flux) были назначены для 161 поверхности.
Расчет температурных деформаций станка. Практика производства показала, что в станках данной группы температурные деформации оказывают определяющее влияние на точность обрабатываемых деталей. Поэтому расчёт температурных деформаций торцешлифовального станка обязательно проводят на стадии проектирования.
Температурное поле станка является причиной возникновения в его несущей системе температурных деформаций. Для этого в Ansys предусмотрено решение задачи термоупругости в двух постановках: несвязанной и связанной [31]. При решении задачи в несвязанной постановке для первоочередного расчёта температурного поля станка используют элементы Solid70 или Solid90. А затем для этой же сеточной модели выполняют решение задачи статики с элементами Solid185/Solid186. Использование конечных элементов типа Solid226 или Solid227 позволяет реализовать связанную задачу термоупругости. В этом случае затраты на разработку итоговой расчётной модели сокращаются. Однако, вычислительные затраты на идентификацию тепловой модели существенно возрастают, так как одновременно строится температурное поле и поле температурных деформаций (или температурных перемещений). Температурные деформации рассчитываются на той же сетке разбиения, что и расчет температурного поля. При расчете температурных деформаций приняты следующие допущения: материал конструкции работает в упругой области; нижнее основание станины жестко закреплено в восьми узлах.
На точность обрабатываемой детали, как отмечалось в [28, 29], влияет изменение взаимного положения рабочих торцов двух шлифовальных кругов относительно диска подачи заготовок по координатам Х, Y, Z (DX, DY, DZ). Перемещения по оси Y компенсируются износом шлифовальных кругов и величиной их подналадки в процессе шлифования. Перемещения по оси Х (в плоскости ХY) практически не оказывают влияния на качество обработки. Наиболее важно определить перемещения по оси Z (в плоскости YZ), оказывающие существенное влияние на точность обрабатываемой детали.
Проверка адекватности принятой расчётной модели. Для подтверждения принятых допущений в расчётной модели и назначенных граничных условий при расчете температурного поля и температурных деформаций был проведен машинный эксперимент для станка модели 3А343АДФ2 при условиях, соответствующих реальным условиям работы станка на холостом ходу и при работе под нагрузкой (при принятой мощности равной 6 кВт). Получено удовлетворительное совпадение результатов эксперимента [29, 30] и расчёта. Погрешность расчётных и экспериментальных данных не превысила 0,5 . Результаты расчёта температурного поля станка, изменение температуры характерных точек и тепловых перемещений шлифовальных кругов приведены на рис. 4–7. На рис. 5 приведены изменения температуры для точек 1, 2, 4, 5, 7, 23, в которых были установлены термодатчики [30].
Из расчётов было получено подтверждение предположения, о симметричности распределения тепловых деформаций (в Ansys – температурных перемещений) станины относительно своего центра (рис. 6). Также было подтверждено, что несимметричность температурных перемещений шлифовальных кругов возникает относительно диска подачи заготовок. Тумба, на которой установлен редуктор диска подачи и закрепленная левее центра станины, имела температурные перемещения того же порядка, что и левая часть станины. Картина температурных перемещений станка на холостом ходу и при работе под тепловой нагрузкой представлена на рис. 6а и 6б, соответственно. Погрешность расчётных и экспериментальных данных не превысила 5 мкм. Оценка адекватности модели проводилась по критерию Фишера.
Рис. 4. Результат расчета температурного поля станка
Рис. 5. Результат изменения температуры характерных точек станка
Рис. 6. Картина тепловых деформаций станка:
а – холостой ход («внизу шире»); б – работа под тепловой нагрузкой («внизу уже»)
Рис. 7. Результат расчета тепловых перемещений шлифовальных кругов
Выводы.
1. Применение известных методик построения упрощенных тепловых и деформационных моделей построения температурных полей и температурных перемещений несущей системы станка позволило провести компьютерное моделирование двустороннего торцешлифовального станка горизонтального исполнения, отличающегося повышенной теплонапряженностью и особенностями конструктивного исполнения.
2. Погрешность моделирования тепловых процессов не превысила 0,5 , погрешность моделирования температурных перемещений не превысила 5 мкм. Проверка адекватности модели была выполнена по критерию Фишера. Особенностью построения расчётной модели являлся учёт двух процессов работы станка – на холостом и рабочем ходах. Хорошее соответствие экспериментальных и расчётных результатов позволяет считать принятые допущения справедливыми и использовать эту расчётную модель торцешлифовального станка для исследования тепловых характеристик торцешлифовальных станков данной гаммы. Полученные результаты моделирования позволяют сделать заключение о возможном использовании данной упрощенной модели для последующих вычислительных исследований влияния конструктивных решений на выходную точность станка и при создании цифрового двойника этого типа станков.
1. Polyakov A.N. Computer research of thermal deformations of metal-cutting machines. Methods, models and algorithms [Komp'yuternye issledovaniya teplovyh deformacij metallorezhushchih stankov. Metody, modeli i algoritmy]. - Orenburg: GOU OGU, 2003. 382 p. (rus)
2. Ito Y. Thermal Deformation in Machine Tools. McGraw Hill Professional. 2010. 240 p.
3. Mayr J., Jedrzejewski J., Uhlmann E., Alkan Donmez M., Knapp W., Härtig F., Wendt K., Moriwaki T., Shore P., Schmitt R., Brecher C., Würz T., Wegener K. Thermal issues in machine tools. CIRP Ann Manuf Technol. 2012. Vol. 61. Is. 2. Pp. 771-791.
4. Kuznecov A.P. Thermal behavior and precision of machine tools [Teplovoe povedenie i tochnost' metallorezhushchih stankov]. M.: MGTU «Stankin», YAnus-K, 2011. 255 p. (rus)
5. Bryan J. International status of thermal error research. CIRP Ann Manuf Technol. 1990. Vol. 39. Is. 2. Pp. 645-656.
6. Kuznecov A.P. Thermal conditions of metal-cutting machines [Teplovoj rezhim metallorezhushchih stankov]. M.: MGTU «Stankin», YAnus-K, 2013. 480 p. (rus)
7. Kuznecov A.P. Thermal processes in metal cutting machines [Teplovye processy v metallorezhushchih stankah]. M.: Tekhnosfera, 2019. 488 p. (rus)
8. Bushuev V.V., Kuznetsov A.P., Khomyakov V.S., Molodtsov V.V. Precision and Efficiency of Metal-Cutting Machines. Russian Engineering Research. 2016. Vol. 36. Is. 9. Pp. 762-773.
9. Bushuev V.V., Kuznetsov A.P., Sabirov F.S., Khomyakov V.S., Molodtsov V.V. Trends in Research on Metal-Cutting Machines. Russian Engineering Research. 2016. Vol. 36. Is. 6. Pp. 488-495.
10. Fu G., Tao C., Xie Y., Lu C., Gao H. Temperature-sensitive point selection for thermal error modeling of machine tool spindle by considering heat source regions. Int J Adv Manuf Technol. 2021. Vol. 112. Pp. 2447-2460.
11. Fan J., Wang P., Tao H., Pan R. A thermal deformation prediction method for grinding machine’ spindle. Int J Adv Manuf Technol. 2021. Vol. 75.
12. Naumann A., Ruprecht D., Wensch, J. Toward transient finite element simulation of thermal deformation of machine tools in real-time. Comput Mech. 2018. Vol. 62. Pp. 929-942.
13. Ge Z., Ding X. Thermal error control method based on thermal deformation balance principle for the precision parts of machine tools. Int J Adv Manuf Technol. 2018. Vol. 97. Pp. 1253-1268.
14. Fang B., Gu T., Ye D., Luo T. An improved thermo-mechanical model for vertical machining center. Int. J. Adv. Manuf. Technol. 2016. Vol. 87. Pp. 2581-2592.
15. Polyakov A.N., Goncharov A.N. Reducing the temperature error of CNC machines based on the motion control of their working bodies [Snizhenie temperaturnoj pogreshnosti stankov s CHPU na osnove upravleniya dvizheniem ih rabochih organov]. - Saratov: Vuzovskoe obrazovanie, 2019. 199 p. (rus)
16. Li Y., Zhao W., Lan S., Ni J., Wu W., Lu B. A review on spindle thermal error compensation in machine tools. International Journal of Machine Tools and Manufacture. 2015. Vol. 95, Pp. 20-38.
17. Liu S., Lin M. Thermal-Mechanical Coupling Analysis and Experimental Study on CNC Machine Tool Feed. Mechanism. Int. J. Precis. Eng. Manuf. 2019. Vol. 20. Pp. 993-1006.
18. Zhang T., Ye W., Liang R., Lou P., Yang X. Temperature variable optimization for precision machine tool thermal error compensation on optimal threshold. Chin. J. Mech. Eng. 2013. Vol. 26. Pp. 158-165.
19. Cao H., Zhu L., Li X., Chen P., Chen Y. Thermal error compensation of dry hobbing machine tool considering workpiece thermal deformation. Int J Adv Manuf Technol. 2016. Vol. 86. Pp. 1739-1751.
20. Li T., Li F., Jiang Y., Wang H. Thermal error modeling and compensation of a heavy gantry-type machine tool and its verification in machining. Int J Adv Manuf Technol. 2017. Vol. 92. Pp. 3073-3092.
21. Ramesh R., Mannan M.A., Poo A.N. Error compensation in machine tools-a review Part II: thermal errors. Int J Mach Tool Manuf. 2000. Vol. 40. Is. 9. Pp .1257-1284.
22. Yao Xd., Du Zc., Ge Gy., Yang Jg. Dynamic temperature gradient and unfalsified control approach for machine tool thermal error compensation. J Mech Sci Technol. 2020. Vol. 34. Pp. 319-331.
23. Huang Z., Liu Y., Du L., Yang H. Thermal error analysis, modeling and compensation of five-axis machine tools. J Mech Sci Technol. 2020. Vol. 34. Pp. 4295-4305.
24. Pavlov S. High performance computing systems in 2019-2020: an overview of achievements and market analysis. Part V. Scope of PLM, including CAE and EDA [Sistemy vysokoproizvoditel'nyh vychislenij v 2019-2020 godah: obzor dostizhenij i analiz rynkov. CHast' V. Sfera PLM, vklyuchaya CAE i EDA]. CAD/CAM/CAE Observer. 2020. No. 7. Pp. 4-19.
25. Berselli G., Bilancia P. Luzi L. Project-based learning of advanced CAD/CAE tools in engineering education. Int J Interact Des Manuf. 2020. Vol. 14. Pp. 1071-1083.
26. Polyakov A.N., Dodorov A.I. Methodology for selecting solid finite element models of machine tool bearing systems during their engineering analysis [Metodika vybora tverdotel'nyh konechno-elementnyh modelej nesushchih sistem stankov pri provedenii ih inzhenernogo analiza]. Bulletin of BSTU named after V.G. Shukhov. 2019. No. 6. Pp. 102-116.
27. Park H., Easwaran A., Andalam S. Challenges in Digital Twin Development for Cyber-Physical Production Systems. In: Chamberlain R., Taha W., Törngren M. (eds) Cyber Physical Systems. Model-Based Design. CyPhy 2018, WESE 2018. Lecture Notes in Computer Science. 2019. Vol. 11615.
28. Nikitina I.P., Polyakov A.N. Problems of heat resistance of bilateral face grinding machines [Problemy teploustojchivosti dvustoronnih torceshlifoval'nyh stankov]. XIII Mezhdunarodnaya nauchno-prakticheskaya konferenciya «Tekhnika i tekhnologii: puti innovacionnogo razvitiya». YUgo-Zapadnyj gosudarstvennyj universitet. 2019. Pp. 6-12. (rus)
29. Nikitina I.P., Polyakov A.N. Experimental investigation of the temperature and accuracy characteristics of a bilateral face grinding machine [Eksperimental'noe issledovanie temperaturnyh i tochnostnyh harakteristik dvustoronnego torceshlifoval'nogo stanka]. Bulletin of BSTU named after V.G. Shukhov. 2019. No. 11. Pp. 112-120. (rus)
30. Nikitina I.P., Polyakov A.N. Features of the flow of thermal processes in double-sided face grinding machines [Osobennosti protekaniya teplovyh processov v dvustoronnih torceshlifoval'nyh stankah]. Bulletin of BSTU named after V.G. Shukhov. 2021. No. 1. Pp. 82-94. (rus)
31. Kozhevnikova M.E., Rotanova T.A., Valov A.V. Computer modeling of plane thermoelasticity problems: comparative analysis of solutions in problems and unrelated formulations [Komp'yuternoe modelirovanie ploskih zadach termouprugosti: sravnitel'nyj analiz reshenij v svyazannoj i nesvyazannoj postanovkah]. Vychislitel'naya mekhanika sploshnyh sred. 2017. Vol. 10. No. 4. Pp. 388-398. (rus)