Conducting a Multifactorial Experiment on the Boriding Process with Simultaneous Carbonization of Ploughshares

Research article
DOI:
https://doi.org/10.23670/IRJ.2023.130.37
Issue: № 4 (130), 2023
Suggested:
18.02.2023
Accepted:
22.03.2023
Published:
17.04.2023
760
5
XML
PDF

Abstract

The article presents the results of experiment planning for determination of optimum process values influencing hardening and increase of service life of machined parts by means of Box – Behnken plan for three factors. Variation levels of the factors under study were selected on the basis of data obtained by mathematical calculations of the course of the process and actual operating conditions of the experimental setup. The calculation of the regression equation factors of the mathematical model of the experiment was carried out by means of the computer environment "Statgraphics Plus for Windows" version 2.1. Using the experimental mathematical model, the optimum range of values of the wear index of the surface of experimental samples, estimated by the loss of thickness in millimetres, was established.

1. Введение

Актуальность темы исследований: Самой энергозатратной работой при обработке почвы является вспашка. По данным Министерства сельского хозяйства РФ общая посевная площадь на территории России увеличивается в среднем на 0,8 млн. га ежегодно на протяжении 5 лет. Ресурс современных рабочих органов почвообрабатывающих орудий в значительной степени зависит от материала заготовки и способа упрочнения. Рабочие органы подвержены, в основном, абразивному изнашиванию. Снизить интенсивность этого вида изнашивания возможно путем применения воздействий, позволяющих произвести упрочнение рабочей поверхности. В связи с этим поиск эффективных технологических процессов упрочнения является актуальной задачей, определяющей цель настоящей работы.

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

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

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

2. Основная часть

Согласно условиям образования боридов при легировании должна соблюдаться интенсивность активирования

:

img

где img – энергия активирования, затраченная на обеспечение массопереноса, Дж/м2∙с;

img – энергия, необходимая для осуществления массопереноса, Дж/м2∙с.

Ток необходимый для образования ванны расплава и передачи энергии для образования боридов

:

img

где img – необходимая энергия для протекания процесса борирования, Дж/м2∙с;

img – эффективный кпд процесса нагрева изделия электрической дугой;

img – коэффициент перевода электрохимических единиц в тепловые;

img - напряжение на электроде в процессе борирования, В.

Общая энергия, затрачиваемая на процесс борирования

,
,
:

img

где img – энергия затрачиваемая на расплавление легирующей пасты, Дж/м2∙с;

img – энергия, затрачиваемая для образования сварочной ванны, Дж/м2∙с;

img – энергия, затрачиваемая на образование боридов, Дж/м2∙с.

Необходимое количество пасты карбида бора

:

img

где img – глубина обработки, м;

img – площадь электрода, м2;

img – плотность получаемого покрытия, кг/м3;

img – требуемое содержание бора в получаемом покрытии, %;

img – коэффициент переноса состава пасты в ванну расплава

.

Расчет энергии требуемой для расплава легирующей пасты

:

img

где img – коэффициент теплопроводности пасты

, Дж/м∙с∙К;

img – толщина наносимой пасты, м;

img – температура плавления пасты, К

.

Расчет энергии, требуемой для расплава стали

:

img

где img – коэффициент теплопроводности стали, Дж/м∙с∙К;

img – температура плавления стали, К.

Расчет толщины наносимой пасты

:

img

где img - объём пасты, содержащей карбид бора, требуемый для легирования, м3.

Объём пасты, содержащей карбид бора, требуемой для легирования

:

img

где img – плотность пасты, вычесленная в результате исследований, составляет 1540 кг/м3 .

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

img

Расчет энергии, требуемой для образования боридов в ванне расплава

:

img

где img – энтальпия образования боридов в ванне расплава

, Дж/кг;

img – масса образовавшегося соединения, г.

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

img

Результаты решения данного уравнения отражены графически на рисунке 1.

Зависимость процентного содержания бора в обрабатываемом материале относительно тока

Рисунок 1 - Зависимость процентного содержания бора в обрабатываемом материале относительно тока

В качестве образцов выступали пластины, нарезанные из лемеха плуга (сталь Л53) высотой 20 мм, шириной 5 мм и толщиной 7 мм. Формирование покрытия проводилось с помощью вибродугового упрочнения. Установка включает пульт управления, вибратор с закрепленным в нем угольным электродом диаметром 10 мм, инверторный источник тока ВДУ-506.
Общий вид установки для борирования с одновременным науглероживанием

Рисунок 2 - Общий вид установки для борирования с одновременным науглероживанием

Сравнительные испытания борированных с одновременным науглероживанием, науглероженных, наплавленных и не упрочненных образцов производились на машине трения горизонтального типа 77 – МТ1 ГОСТ 10198 – 78. Из науглероженного лемеха были вырезаны образцы сечением 5×5×20 мм. В специально изготовленной оправке образец закреплялся в неподвижном патроне машины и прижимался к бруску из белого электрокорунда 25АСМ2К20М ГОСТ 27595 – 88 зернистостью 127. Давление образца на брусок было постоянно и составляло 0,2 МПа. Через каждые 5000 циклов (двойных ходов) производился замер высоты образца микрометром МК – 25 ГОСТ 6507 – 90 с точностью 0,01 мм и производилось взвешивание на весах ВЛР – 200 ГОСТ 24104 – 80 с точностью 0,0001 г.

Уровни варьирования исследуемых факторов были выбраны на основе данных, полученных при математических расчетах протекания процесса, фактических условий режима работы экспериментальной установки и представлены в таблице 1.

Таблица 1 - Уровни и интервалы варьирования исследуемых факторов

Наименование фактора

Кодированное обозначение фактора

Уровни варьирования

Интервал варьирования

-1

0

+1

Ток I, А

х1

140

180

220

40

Количество рабочих циклов нагрузки Nц, шт.

х2

5000

20000

35000

15000

Толщина слоя пасты hп, мм

х3

0,1

0,2

0,3

0,1

Для реализации многофакторного эксперимента был выбран оптимальный план Бокса – Бенкина для трёх факторов

. Матрица эксперимента отражена таблицей 2.

Таблица 2 - Матрица оптимального плана многофакторного эксперимента

№ п/п опыта

Факторы

Параметр оптимизации

х1

х2

х3

Н, мм

1

1

1

0

0,79

2

-1

-1

0

0,27

3

1

-1

0

0,28

4

-1

1

0

0,74

5

1

0

1

0,60

6

-1

0

-1

0,63

7

1

0

-1

0,62

8

-1

0

1

0,61

9

0

1

1

0,78

10

0

-1

-1

0,27

11

0

1

-1

0,77

12

0

-1

1

0,26

13

0

0

0

0,56

14

0

0

0

0,58

15

0

0

0

0,57

Расчёт коэффициентов уравнения регрессии математической модели эксперимента был осуществлён при помощи компьютерной среды «Statgraphics Plus for Windows» version 2.1, в результате чего было получено уравнение регрессии экспериментальной математической модели в виде полинома первой степени:

img

Для определения координат точки центра поверхности отклика, продифференцировали уравнение математической модели:

img

Решая систему линейных уравнений, были определены координаты оптимальной точки – центра поверхности отклика математической модели в кодированном виде:

img

img

img

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

img

img

img

При этом значение износа поверхности экспериментальных образцов составляет

img

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

Поверхность отклика математической модели, построенная в координатах: а – силы тока (х1) и количества рабочих циклов нагрузки (х2); б – силы тока (х1) и толщины слоя ппасты (х3); в – количества рабочих циклов нагрузки (х2) и толщины слоя пасты (х3)

Рисунок 3 - Поверхность отклика математической модели, построенная в координатах:

а – силы тока (х1) и количества рабочих циклов нагрузки (х2); б – силы тока (х1) и толщины слоя ппасты (х3); в – количества рабочих циклов нагрузки (х2) и толщины слоя пасты (х3)

Двухмерные сечения поверхности отклика математической модели изолиниями равного выхода, построенные в координатах: а – силы тока (х1) и количества рабочих циклов нагрузки (х2); б – силы тока (х1) и толщины слоя покрытия (х3); в – количества рабочих циклов нагрузки (х2) и толщины слоя покрытия (х3)

Рисунок 4 - Двухмерные сечения поверхности отклика математической модели изолиниями равного выхода, построенные в координатах:

а – силы тока (х1) и количества рабочих циклов нагрузки (х2); б – силы тока (х1) и толщины слоя покрытия (х3); в – количества рабочих циклов нагрузки (х2) и толщины слоя покрытия (х3)

Поверхность отклика экспериментальной математической модели была изучена в области «условного» оптимума с использованием элементов графоаналитического метода двухмерных сечений. Оптимальные сочетания значений исследуемых факторов по изолиниям с минимальными значениями, в соответствии с рис. 2, на каждом из двухмерных сечений, определятся следующим образом:

- двухмерное сечение а) img

- двухмерное сечение б) img

- двухмерное сечение в) img

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

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

img

img

img

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

img

img

img

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

img

3. Заключение

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

img

2. Интерпретируя математическую модель, установлено влияние исследуемых факторов на параметр оптимизации.

3. Изучена поверхность отклика экспериментальной математической модели в области «условного» оптимума с использованием элементов графоаналитического метода двухмерных сечений.

4. Определено оптимальное сочетание диапазонов значений исследуемых факторов: I∈[170;181] А; Nц∈[5750;7100] шт; hп∈[0,195;0,225] мм, при которых отклик изучаемой системы достигает своего наименьшего значения.

5. Используя экспериментальную математическую модель, установлен оптимальный диапазон значений показателя износа поверхности экспериментальных образцов, оцениваемый потерей толщины в миллиметрах: Н = 0,27 … 0,3 мм.

Article metrics

Views:760
Downloads:5
Views
Total:
Views:760