UPGRADING OF THERMAL FRICTION CUTTING OPERATIONS DUE TO OPTIMAL STOCHASTIC DIAGNOSTICS AND STOCK REMOVAL CONTROL
Abstract and keywords
Abstract (English):
An approach that improves the quality of work on a workpiece under thermal friction cutting operations is viewed. The thermal friction cutting is presented as a synergetic object in nonequilibrium conditions, due to the influence of anisotropic properties of structural materials, that results in instability of the surface quality of products. The paper provides a mathematical description of the object, dynamic equations are constructed, and the object is represented in the form of a state space. To upgrade machinning under thermal friction cutting operations, an approach based on the synthesis of a tracking system with a Lewinberger observer is proposed, which provides operational information with the minimum allowable RMS error. The application of the classical approach to the synthesis of stochastic observers leads to significant errors due to the multiple differentiation of measurement noise. To eliminate this disadvantage, a method has been proposed that eliminates the process of multiple differentiation by rebuilding a system that does not contain instrument noise in its structure. For that purpose, a renewal equation is introduced in the state space. Based on the received operational information, a system of dynamic stabilization of the TFC process has been built. The proposed system was tested on an OPTIMUS F3 CNC machine, equipped with a set of metering equipment for measuring the relative motion of the tool, the workpiece and the temperature change in the contact area. An analysis of the results shows that for machining upgrading under thermal friction cutting operations, it is advisable to use feedback systems equipped with stochastic optimal estimation and output control units during stock removal control. Testing of the proposed control systems gives evidence of bettering the quality of the machined surface, for example, a reduction in roughness by half in terms of Rz and an increase in the stability of the microhardness index by 2 to 8 times in the TFC of 45 steel samples.

Keywords:
thermal friction cutting (TFC), contact zone, temperature fields, stochastic observer, anisotropic medium, synergetic approach
Text
Text (PDF): Read Download

Введение

 

Качество обрабатываемых изделий зависит как от качества заготовки, технологии обработки, так и от свойств технологической системы. В настоящее время для черновой и получистовой обработки заготовок стал широко применяться метод термофрикционного резания, который по экономическим критериям имеет значительные преимущества в сравнении с традиционными операциями, такими как фрезерование, обдирочное шлифование и др. [1 − 3].

Процесс термофрикционного резания (ТФР) представляет процесс формообразования путем механического удаления нагретого до определённой температуры слоя материала режущим диском. Резание материалов сопровождается высокой степенью и скоростью пластической деформации на фоне интенсивных тепловых потоков и активной технологической среды. Процесс пластической деформации и контактного взаимодействия обрабатываемого материала с материалом режущего диска генерирует ряд физико-химических явлений, таких как термоэлектрические, магнитные, диффузионные, химические и другие. Все это указывает на высокую степень неравновесности, вызванных наведением различных градиентных соотношений, которые вызывают развитие энергетических потоков различной природы, в том числе массопереноса на разных масштабных структурных уровнях [7, 8, 10].

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

При решении такой системы в качестве обобщенной координаты принимается направление, оказывающее наиболее сильное влияние на точность и качество механической обработки.

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

Сам процесс ТФР материалов можно разбить на отдельные алгоритмы, протекающие в динамическом единстве на разных масштабных структурных уровнях, при этом все они могут быть разбиты по конкурирующим парам. Такое разделение алгоритмов с философской точки зрения отвечает закону единства и борьбы противоположностей, и раскрывает их роль в обеспечении принципа наименьшего принуждения и других синергетических парадигм [1, 4, 5]. Соблюдение баланса между конкурирующими алгоритмами обеспечивает динамическую стабильность и устойчивость функционирования всей системы резания.

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

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

Анализ современных исследований в области обработки материалов резанием показывает, что одним из путей повышение качества обработки на операциях термофрикционного резания является синтез систем стохастической диагностики позволяющих получить рабочую информацию об процессе обработки  с наименьшей из возможных дисперсий, на основе которых может быть синтезирована системах управления по отклонению, которая путем обработки полученных оценок оптимальным образом осуществляет стабилизацию выходных показателей качества [6, 9].

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

 

Материалы и методы

 

Для моделирования динамки процесса ТФО рассмотрим схему, представленную на
рис. 1.

На основе принципа возможных перемещений, для системы, представленной на рис. 1, построено ее представление в виде совокупности дифференциальных уравнений, характеризующих динамику перемещений центров круга и детали в процессе круглого наружного шлифования:

где m1 , m2  – приведенные массы систем «заготовка» и «инструмент»; λ1  – коэффициент демпфирования i -го звена; ci  – параметр жесткости i -го звена; y1 , y2  – приращения координат центра масс инструмента и центра масс заготовки соответственно; ΔHΔh  – приращения расстояний от положения центра масс инструмента и заготовки, соответственно, вызванные торцевым биением инструмента и неравномерностью припуска на заготовке; L ,  ΔL  расстояние между центрами масс и его изменение.

Для упрощения записи приращений
координат и их производных в зависимостях (1), опущен предшествующий им знак приращения
Δ .

где Y0=y1y2y3y4  – вектор производных состояний системы; Y0=y1y2y3y4  – вектор (матрица-столбец), представляющая вектор состояния системы; B0=B01B02  – матрица параметров влияния торцевого биения, инструмента и неравномерности припуска заготовки; Ψ  – вектор состояний отклонений значений торцевого биения инструмента и неравномерности припуска от заданных значений Ψ=Ψ1Ψ2,  
Ψ1=[ΔR+Δr],  Ψ2=[R+r] ; U=SS  – вектор управления глубиной резания, связанный с температурой зоны контакта заготовки и инструмента;

 

A0=0100-c1+c3m1-h1+h3m1c3m1h3m10001c3m2h3m2-c2+c3m2-h2+h3m2  –

 

матрица, характеризующая динамические свойства системы; C0=0c1m1000h1m100  – матрица управления температурой в зоне резания (за счет поперечной подачи; E0=10000100  – матрица состояния процесса измерений;
F0=f11f12f21f22  – матрица интенсивностей (погрешностей измерителей); V0=v1v2  – матрица шумов измерителей единичной интенсивности; T0=Т(tf)  – температурная матрица; Q0=-11  – матрица интенсивности совокупных измерений.

Полученные данные о значениях температуры в зоне резания отличаются от реальных значений температуры зоны контакта из-за качества процесса измерений и погрешностей измерителей, а также от случайных факторов вызванных отклонениями значений припуска и торцевого биения диска от номинальных значений. Поэтому непосредственное получение скорости изменения температуры dТ(tf)dt  дифференцированием температурной функции нецелесообразно вследствие присутствия аддитивных шумов в совокупных измерениях.

Для решения задачи оценки температуры в зоне резания и скорости ее изменения может быть использована теория стохастических наблюдателей, на основе которой может быть построена стохастический наблюдатель Льюинбергера.

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

Непосредственное применение матричных представлений (3) при моделировании нецелесообразно вследствие необходимости дифференцирования. Это не является ни особенностью матричного представления модели, ни описания динамической системы в отклонениях, а проявляется уже в исходном представлении (1) технологической системы.

Для разрешения вышеописанной задачи преобразуем вектор состояния системы к виду:

X0=D0Y0 ,                          (3)

 

где D0  – матрица интенсивности оценок параметров состояния процесса ТФО;

Тогда модифицированный вектор состояния процесса ТФР примет вид:

 

F=Y0-B02Ψ1 .                                  (4)

 

Из соотношения следует

 

Y0=F+B02Ψ1 .                                  (5)

 

Тогда вектор Y0  характеризующий скорость изменения состояния системы ТФР может быть описан выражением вида:

 

Y0=F+B02ψ1+B02ψ2  .                 (6)

 

С учетом преобразований, матричное уравнение (2) запишется:

 

Y0=A0Y0+B01Ψ1+B02Ψ2+C0U.      (7)

После подстановки выражения (F+B02Ψ1 ) в уравнение (7) получим:

 

 

Y0=A0F+(B02Ψ1)+B01Ψ1+B02Ψ2+C0U.

 

На основе выражений (6) и (7) запишем равенство:

 

F+B02ψ1+B02ψ2=A0F+B02Ψ1+B01Ψ1+B02Ψ2+C0U.               (8)

 

 

На основании полученного результата после преобразований запишем модифицированное математическое выражение характеризующее состояние процесса ТФР, не содержащее шумов измерений Ψ2 . В связи с чем исключается необходимость в дифференцировании Ψ1 :

 

 

F=A0F+[A0B02+B01-B02]Ψ1+C0U.                                                     (9)

 

 

Процедура сравнения уравнений (9) и (3) позволяет определить коэффициент, стоящий перед матрицей Ψ1 , и интерпретировать модифицированную матрицу влияния отклонения торцевого биения инструмента и неравномерности припуска заготовки:

 

B1=[A0B02+B01-B02]              (10)

Для этих условий модифицированное математическое выражение, характеризующее состояние процесса ТФР примет вид:

 

F=A0F+B1Ψ1+C0U.              (11)

 

Полученное соотношение и первое уравнение системы (2) являются эквивалентными в смысле Ляпунова.

Тогда на основании выражений (7), (4) и (11) получим соотношение для оценки
состояния
X0 :

 

X0=D0F+D0B02Ψ1              (12)

 

где F  – модифицированный вектор состояния, определяется из (11).

 

Результаты и их обсуждение

 

Рассмотрим разработанный подход
для синтеза стохастического наблю-
дателя Льюинбергера по форме фильтра
Калмана-Бюсси.

Выполним преобразование произведени матриц C0U0  системы (11) к виду

 

С0U0=С00S+С01S,

 

где С00=c1m1T , С01=0h1m100T ,

 

и перестроим модель системы (8) в форме, не требующей дифференцирования сигнала управления.

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

 

 

F=A0F+B1Ψ1+C1U;   
Z0=E0F+F0V0+DMS;                                                                           (13)
T0=Q0Z0,        

 

где F=F1F2F3F3T,  C1=0h3h2m1m20(h3+h2)h2m22-c2m2T , DM=0h2m2T .

 

Неравномерность припуска на поверхности заготовки и случайные колебания торца диска от заданных значений в процессе термофрикционного резания могут быть смоделированы в виде марковских случайных процессов второго порядка с корреляционными функциями типа:

 

R=De(-αvτ)cos⁡(βvτ) ,                    (14)

где D  − дисперсия высот неровностей рельефа заготовки, v  − скорость движения инструмента в процессе ТФР, α , β  − соответствующие коэффициенты корреляции.

Для условий ТФР плоских заготовок построен расширенный формирующий фильтр, которому в пространстве состояний соответствует система, состоящая из трех уравнений:

 

Для компактности записи модели первые два уравнения системы (16) приведем к форме:

 

X=AX+BU+EW                  (17)
Z=CX+RV

 

На основе математического описания объекта (16) может быть построен стохастический наблюдатель, позволяющий получать оптимальные в среднеквадратическом смысле оценки значений температуры в зоне резания при ТФР, путем решения матричного уравнения Риккати:

 

 

V=AV+VAT+EΨwET-VCTΨv-1CV.                                                (18)

 

Матрица коэффициентов усиления фильтра Калмана определяется следующей зависимостью

 

K=VCΨv-1                                                                                                (19)

 

Тогда алгоритм фильтрации наблюдений определяется матричными уравнениями

 

Y=AY+BU+K[Z-BU-CX                                                  (20)

 

Т(tf)Т(tf)=-101000000-1010000Y   ∧    

 

 

Соотношения (18), (19) позволяют использовать результаты измерений координат в процессе обработки детали с целью построения оценок параметров процесса обработки. Эти оценки являются оптимальными при условии гауссовости шумов измерений и возбуждений.

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

Осуществим синтез такой системы для термофрикционного резания.

Для синтеза оптимального детерминированного регулятора может быть использован подход, изложенный в работе [2]. Данная методика предполагает выбор критерия эффективности.

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

 

 

I=M12xTtQ1xt+120tf(xTQ2x+uTQ3-1)dt                                        (21)

 

 

где Q1 , Q2 , Q3 , соответствующие весовые матрицы качества управления.

Матрица коэффициентов усиления управления регулятора определяется следующей зависимостью:

На основании (17), (19) и (22), замкнутый контур регулирования, состоящий из объекта, фильтра и регулятора, для технологической системы ТФР, характеризуется обобщенной системой дифференциальных уравнений:

 

 

XX=F-CB1BHF-CB1-BH×XX+G0×W+0BΩv

 

 

Заключение

 

Испытание предложенной системы осуществили на станке модели OPTIMUM F3 CNC, оснащённом комплексом измерительной аппаратуры для измерения относительных колебаний инструмента, заготовки и температуры в зоне контакта (рис. 2).

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

В качестве образцов использовали заготовки из стали 45, на поверхности которых после разметки были проведены измерения микротвердости в 10 точках с интервалом 5 мм (табл. 2).

Анализ полученных результатов показывает, что для повышения качества обработки на операциях термофрикционного резания целесообразно использовать системы с обратной связью, оснащённые блоками стохастического оптимального оценивания и управления выходными показателями в процессе съема припуска на обработку.

Выполнение испытания предложенных систем управления свидетельствуют об улучшении качества обработанной поверхности, например снижении шероховатости в два раза по показателю Rz и повышению стабильности показателя микротвёрдости в 2 – 8 раз при ТФР образцов из стали 45.

Полученные результаты быть использованы для синтеза систем управления операциями ТФР для других материалов и могут служить основой для создания более сложных систем.

 

References

1. Zarubitsky E.U., Plakhotnik V.A., Pokintelitsa N.I. et.al. Optimization of thermal friction cutting modes // Bulletin of Mechanical Engineering. 1998, no. 9. pp. 54–56.1

2. Pokintelitsa N.I., Strelyanaya Yu.O., Bratan M.I. Plastic deformation during thermofriction processing of ma-terials. Scientific notes of the Crimean Engineering Peda-gogical University. 2022, no. 3 (77), pp. 188–191. DOIhttps://doi.org/10.34771/UZCEPU.2022.77.3.036.

3. Zarubitsky E.U. The temperature of the removable allowance under thermofriction cutting // Optimization of cutting processes of heat-resistant and extra strong materials. Ufa: UAI Publishing House, 1986, pp. 106–110.

4. Volkov O.A. Investigation of the thermal deformation effect during surface hardening of steels by thermal friction treatment // VEZHPT. 2016, no. 5 (80). URL: https://cyberleninka.ru/article/n/issledovanie-teplodeformatsionnogo-vliyaniya-pri-poverhnostnom-uprochnenii-staley-termofriktsionnoy-obrabotkoy.

5. Talantov N.V. The physical foundations of the process of cutting, tool wear and destruction. Moscow: Mashinostroenie, 1992, 240 p.

6. Kostetsky B.I., Natanson M.E., Bershadsky L.I. Mechanical and chemical processes under boundary friction. Moscow: Nauka, 1972, 170 p.

7. Balakin V.A. Friction and wear at high sliding speeds. Moscow: Mashinostroenie, 1980, 136 p.

8. Bratan S.M., Pokintelitsa N.I., Yakubov C.F., Chasovitina A.S. Method of mathematical description of tool geometry estimation taking into account thermophysical parameters // Fundamental and applied problems of engineering and technology. 2024, no. 2 (364), pp. 18–23. DOIhttps://doi.org/10.33979/2073-7408-2024-364-2-18-23.

9. Pokintelitsa N.I. Bratan M.I., Yakubov Ch.F. Methodology for ensuring the effectiveness of frictional shaping of rotations. Scientific notes of the Crimean Engineering Pedagogical University, 2023, no. 1 (79), pp. 252–260. DOIhttps://doi.org/10.34771/UZCEPU.2023.79.1.050.EDN FOCKNK.

10. Yakubov Ch.F., Pokintelitsa N.I., Bratan S.M., Chasovitina A.S. Modeling of heat transfer processes under thermal friction cutting in mild steel blanks // Science-inrensive technologies in mechanical engineering. 2024, no. 10 (160), pp. 19–26. DOIhttps://doi.org/10.30987/2223-4608-2024-19-30.

Reviews
1. Upgrading of thermal friction cutting operations due to optimal stochastic diagnostics and stock removal control Authors: Fedorov Vladimir

Login or Create
* Forgot password?