ПОВЫШЕНИЕ КАЧЕСТВА ОБРАБОТКИ НА ОПЕРАЦИЯХ ТЕРМОФРИКЦИОННОГО РЕЗАНИЯ ЗА СЧЕТ ОПТИМАЛЬНОЙ СТОХАСТИЧЕСКОЙ ДИАГНОСТИКИ И УПРАВЛЕНИЯ СЪЕМОМ ПРИПУСКА
Аннотация и ключевые слова
Аннотация (русский):
Рассмотрен подход, обеспечивающий повышение качества обработки заготовок на операциях термофрикционного резания. Процесс термофрикионного резания представлен как синергетический объект, находящийся в неравноместных условиях, ввиду влияние анизотропных свойств конструкционных материалов, что приводит к нестабильности качества поверхности изделий. В работе выполнено математическое описание объекта, построены динамические уравнения и объект представлен в форме пространства состояний. Для повышения качества обработки изделий на операциях термофрикицонного резания предложен подход, основанный на синтезе следящей системы с наблюдателем Льюинбергера, которая обеспечивает получение рабочей информации с минимально допустимой среднеквадратической ошибкой. Применение классического подхода к синтезу стохастических наблюдателей приводит к существенный погрешностям ввиду многократного дифференцирования шумов измерителей. Для устранения этого недостатка предложен метод, исключающий процесс многократного дифференцирования, путем перестройки системы, не содержащей в своей структуре шумов измерителей. Для этого в пространстве состояний введено уравнение восстановления. На основании полученной рабочей информации построена система динамической стабилизации процесса ТФР. Испытание предложенной системы осуществили на станке модели OPTIMUM F3 CNC, оснащённом комплексом измерительной аппаратуры для измерения относительных колебаний инструмента, заготовки и температуры в зоне контакта. Анализ полученных результатов показывает, что для повышения качества обработки на операциях термофрикционного резания целесообразно использовать системы с обратной связью, оснащённые блоками стохастического оптимального оценивания и управления выходными показателями в процессе съема припуска на обработку. Выполнение испытания предложенных систем управления свидетельствуют об улучшении качества обработанной поверхности, например снижении шероховатости в два раза по показателю Rz и повышению стабильности показателя микротвёрдости в 2 – 8 раз при ТФР образцов из стали 45.

Ключевые слова:
термофрикционое резание, зона контакта, поля температур, стохастический наблюдатель, анизотропная среда, синергетический подход
Текст
Текст (PDF): Читать Скачать

Введение

 

Качество обрабатываемых изделий зависит как от качества заготовки, технологии обработки, так и от свойств технологической системы. В настоящее время для черновой и получистовой обработки заготовок стал широко применяться метод термофрикционного резания, который по экономическим критериям имеет значительные преимущества в сравнении с традиционными операциями, такими как фрезерование, обдирочное шлифование и др. [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.

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

 

Список литературы

1. Зарубицкий Е.У., Плахотник В.А., Покинтелица Н.И. и др. Оптимизация режимов термофрикционного резания // Вестник машиностроения. 1998. № 9. С. 54–56.

2. Покинтелица Н.И., Стреляная Ю.О., Братан М.И. Пластическое деформирование при термофрикционной обработке материалов // Ученые записки Крымского инженерно-педагогического университета. 2022. № 3 (77). С. 188–191. DOIhttps://doi.org/10.34771/UZCEPU.2022.77.3.036. EDN XHTITS.

3. Зарубицкий Е.У. Температура снимаемого припуска при термофрикционном резании // Оптимизация процессов резания жаропрочных и особо прочных материалов. Уфа: Изд-во УАИ, 1986. С. 106–110.

4. Волков О.А. Исследование теплодеформационного влияния при поверхностном упрочнении сталей термофрикционной обработкой // ВЕЖПТ. 2016. № 5 (80). URL: https://cyberleninka.ru/article/n/issledovanie-teplodeformatsionnogo-vliyaniya-pri-poverhnostnom-uprochnenii-staley-termofriktsionnoy-obrabotkoy.

5. Талантов Н.В. Физические основы процесса резания, изнашивания и разрушения инструмента. М.: Машиностроение, 1992. 240 с.

6. Костецкий Б.И., Натансон М.Э., Бершадский Л.И. Механико-химические процессы при граничном трении. М.: Наука, 1972. 170 с.

7. Балакин В.А. Трение и износ при высоких скоростях скольжения. - М.: Машиностроение, 1980. 136 с.

8. Братан С.М., Покинтелица Н.И., Якубов Ч.Ф., Часовитина А.С. Методика математического описания оценки геометрии инструмента с учетом теплофизических параметров // Фундаментальные и прикладные проблемы техники и технологии. 2024. № 2 (364). С. 18−23. DOIhttps://doi.org/10.33979/2073-7408-2024-364-2-18-23.

9. Покинтелица Н.И., Братан М.И., Якубов Ч.Ф. Методология обеспечения повышения эффективности фрикционного формообразования поверхностей // Ученые записки Крымского инженерно-педагогического университета. 2023. № 1 (79). С. 252–260. DOIhttps://doi.org/10.34771/UZCEPU.2023.79.1.050. EDN FOCKNK.

10. Якубов Ч.Ф., Покинтелица Н.И., Братан С.М., Часовитина А.С. Моделирование процессов теплопередачи при термофрикционном резании в заготовках из конструкционных сталей // Наукоемкие технологии в машиностроении. 2024. № 10 (160). С. 19−26. DOIhttps://doi.org/10.30987/2223-4608-2024-19-30.

Рецензии
1. Повышение качества обработки на операциях термофрикционного резания за счет оптимальной стохастической диагностики и управления съемом припуска Авторы: Федоров Владимир Павлович

Войти или Создать
* Забыли пароль?