Постоянное усложнение формы деталей, и, как следствие, повышенные требования к их точности и надёжности, в машиностроении - неотъемлемая составляющая технического прогресса. Для повышения надежности требуется применение более прочных материалов, обработка которых классическими методами затруднена. К таким материалам относятся различные высокопрочные материалы, нержавеющие и жаропрочные стали, твердые и титановые сплавы, широко применяемые в авиа-, судо-, автомобиле-, двигателе- и приборостроении, инструментальном производстве [1, 2]. А для повышения точности необходимо учитывать все значимые процессы, происходящие при обработке, что, в свою очередь, позволит снизить погрешность производства деталей. Вышеназванные требования можно выполнить благодаря использованию методов обработки, основанных на электрофизикохимических воздействиях. Одним из таких методов является электроэрозионная обработка (ЭЭО).
Особенностью ЭЭО является то, что обработка ведется импульсами тока, в результате чего происходит удаление материала заготовки. Данный метод обработки позволяет получать изделия практически любой конфигурации с требуемыми параметрами качества поверхностного слоя со значительно меньшими затратами, чем при механической обработке.
Одним из наиболее распространенных методов ЭЭО является метод электроэрозионной прошивки, обеспечивающий возможность изготовления глубоких отверстий различной конфигурации. При прошивании глубоких отверстий основной задачей является поддержание процесса и его интенсификация при значительном внедрении электрода-инструмента в заготовку, когда условия эвакуации продуктов обработки из межэлектродного промежутка затруднены [3]. При прошивании отверстия частицы продуктов эрозии, двигаясь между уже обработанной поверхностью детали и образующей поверхностью электрода-инструмента, вызывают дополнительные электрические разряды, что ведет к возникновению конусности. Для решения данной проблемы используют станки с автоматическим выводом электрода-инструмента (релаксацией) из отверстия для очистки его от продуктов эрозии, а также прокачку или отсос рабочей жидкости.
Следует отметить, что ЭЭО является сложным физическим процессом, экспериментальное исследование которого в значительной степени затруднено вследствие сложности внедрения измерительных приборов непосредственно в зону обработки. Однако процесс ЭЭО можно представить в виде совокупности отдельных процессов, протекание которых подчиняется известным физическим законам, каждый из которых можно описать с помощью математических моделей. Поэтому одним из перспективных методов исследования ЭЭО является имитационное компьютерное моделирование, основанное на математических моделях физических процессов, происходящих в зоне обработки.
Имитационное компьютерное моделирование процесса электроэрозионной обработки
Этапы имитационного компьютерного моделирования процесса ЭЭО:
-
Назначаются размеры рабочей области и ячейки сетки, задаются технологические параметры процесса ЭЭО, затем на основании этих данных строится сеточная модель рабочей области.
-
Используя сеточную модель, производится расчет напряженности электрического поля в ячейках сетки в зависимости от текущего положения электродов и концентрации продуктов эрозии в рабочей жидкости.
-
Рассчитывается электрическая прочность локальных участков рабочей жидкости, в зависимости от концентрации содержащихся в них продуктов обработки, и производится проверка возможности возникновения электрического пробоя между электродом-инструментом (ЭИ) и электродом-заготовкой (ЭЗ).
-
Если для данного положения электродов и концентрации продуктов эрозии пробой межэлектродного промежутка возможен, то рассчитывается траектория образования электрической дуги.
-
В местах соприкосновения дуги с электродами формируются эрозионные лунки, и рассчитываются направления разлета продуктов эрозии, образовавшихся в лунках.
-
С учетом вновь образовавшихся продуктов обработки производится глобальный расчет движения всех продуктов эрозии в рабочей жидкости. При этом в случае использования принудительной прокачки рабочей жидкости, всем частицам продуктов эрозии придается импульс, направление которого совпадает с направлением прокачки.
-
Затем снова производится расчет напряженности электрического поля в ячейках сетки в зависимости от текущей концентрации продуктов эрозии, и процесс моделирования итерационно повторяется.
-
Если в процессе ЭЭО межэлектродный зазор увеличился настолько, что образование электрической дуги становится невозможно, то происходит опускание электрода-инструмента вниз на высоту одной ячейки сеточной модели.
-
Если на этом этапе достигнута заданная глубина прошиваемого отверстия или же достигнуто «дно» сеточной модели, то процесс моделирования прекращается.
-
Если заданная глубина прошиваемого отверстия не достигнута, то работа компьютерной модели продолжается. Следует отметить, что процесс опускания электрода-инструмента сопряжен с дополнительным приложением импульса к частичкам продуктов обработки, взвешенным в рабочей жидкости, поэтому после этого этапа необходимо снова произвести глобальный расчет движения продуктов эрозии в рабочей жидкости.
Для реализации предложенного алгоритма имитационного компьютерного моделирования процесса ЭЭО был разработан комплекс математических моделей сопряженных физических процессов при единичном электрическом разряде, включающий в себя модель распределения электрических потенциалов в зоне обработки, модель пробоя межэлектродного промежутка (МЭП), модель движения продуктов эрозии в межэлектродном зазоре.
Моделирование процесса ЭЭО производится в соответствии с алгоритмом, представленным на рис. 1.
Рис. 1. Алгоритм компьютерного моделирования процесса электроэрозионной обработки
Сеточная модель распределения электрических потенциалов в рабочей области
Для построения сеточной модели рабочей области задается диаметр электрода-инструмента Dи, глубина прошиваемого отверстия HО (рис. 2), тип рабочей жидкости (РЖ) и напряжение, подаваемое на электроды. Исходя из электрических свойств используемой РЖ и напряжения импульсов, рассчитывается номинальный (для чистой РЖ) размер МЭП hМЭП:
(1)
где Uи - напряжение импульса, В; Eпр - электрическая прочность незагрязненной РЖ, В/м.
Рис. 2. Сеточная модель рабочей области: 1 - электрод-инструмент; 2 - электрод-заготовка; 3 - рабочая жидкость
Исходя из рассчитанного hМЭП, выбирается размер ячейки сетки l таким образом, чтобы МЭП укладывался в целое число ячеек сетки и составлял не менее трех ячеек. При использовании в качестве РЖ керосина и напряжении импульсов, равном 200 В, размер ячейки сетки выбирается в интервале от 1 до 10 мкм и размер МЭП составляет от 3 до 28 ячеек.
На основании геометрических размеров электродов Dи, HО и величины МЭП задаются ширина W и высота H сеточной модели, выраженные в количестве ячеек.
Расчет электростатического поля в зоне обработки производится методом конечных разностей. При этом непосредственно для расчета электрических потенциалов ячеек используется метод резистивных сеток, где в качестве сопротивлений выступают величины, обратные диэлектрическим проницаемостям ячеек сетки. Для ячеек, соответствующих электродам, диэлектрическая проницаемость принимается равной большому конечному числу.
Пример распределения электрических потенциалов в зоне обработки до начала расчета и после завершения вычисления электростатического поля приведен на рис. 3.
Модель пробоя межэлектродного промежутка
После расчета электростатического поля происходит расчет пробоя локальных участков РЖ. Для этого вычисляется электрическая прочность ячеек МЭП. В данной модели каждая ячейка РЖ может быть загрязнена продуктами эрозии, причем концентрация загрязнения может составлять от 0 до 100%.
Рис. 3. Распределение электрических потенциалов в зоне обработки: а - до начала расчета электростатического поля; б - после вычисления электростатического поля
Принимая допущение, что продукты эрозии являются частицами проводника, и их наличие в РЖ уменьшает линейные размеры той области ячейки, которая занята диэлектриком, введем в модель линейную зависимость электрической прочности ячейки от концентрации КЗ продуктов эрозии в ней:
(2)
где Eпр.з - напряженность пробоя ячейки РЖ, загрязненной продуктами эрозии, В/м; Eпр - номинальная напряженность пробоя чистой РЖ, В/м.
Учитывая выражение (2), определяется напряженность электрического поля для ячеек РЖ.
Полученная математическая модель пробоя межэлектродного промежутка имеет вид:
где hМЭП - номинальный размер МЭП, м; hпр - размер пробиваемого МЭП, м; Uи - напряжение импульса, В; Eпр - номинальная напряженность пробоя чистой РЖ, В/м; Eпр.з - напряженность пробоя РЖ, загрязненной продуктами эрозии, В/м; КЗ - концентрация продуктов эрозии в РЖ, %.
Следует отметить, что наибольшая напряженность электрического поля наблюдается на углах ЭИ. Если для какой-либо ячейки напряженность электрического поля больше, чем напряженность пробоя, то разряд в МЭП возможен, и происходит расчет траектории образования электрической дуги.
В данной компьютерной модели эрозионные лунки в ЭЗ образуются после каждого разряда, а в ЭИ - с вероятностью, равной коэффициенту относительного износа инструмента, то есть, если в модели задан коэффициент относительного износа ЭИ, равный 6%, то эрозионная лунка будет образовываться в среднем после попадания 1/0,06 » 16 разрядов в одну и ту же ячейку ЭИ.
Модель движения продуктов эрозии в межэлектродном зазоре
Для имитации движения продуктов эрозии в РЖ в данной компьютерной модели используется клеточно-автоматная модель FHP-EP (erosion products). Модель FHP-EP базируется на клеточно-автоматной модели потока FHP-MP [4], в которой допускается более одной частицы в клетке с равными векторами скорости. Модель FHP-МP является обобщением классической модели на булевых векторах FHP (Frish, Hasslacher, Pomeau) [5].
Рис. 4. Направления движения рабочей жидкости при различных режимах принудительной прокачки: 1 - электрод-инструмент; 2 - электрод-заготовка; 3 - рабочая жидкость
Для имитации прокачки РЖ через МЭП в данной компьютерной модели также используется клеточно-автоматная модель FHP-EP. Каждой частице продуктов эрозии придается дополнительный импульс в соответствии с направлением прокачки РЖ.
После каждой операции расчета движения продуктов эрозии цикл работы компьютерной модели повторяется, то есть происходит новый расчет электростатического поля в зоне обработки с учетом появления новой эрозионной лунки и нового положения продуктов эрозии, а затем производится новый расчет пробоя РЖ.
Если при расчете пробоя выяснится, что пробой невозможен, вследствие слишком большого расстояния между электродами, то модель произведет опускание ЭИ вниз на высоту одной ячейки сеточной модели. При этом частицам продуктов эрозии, расположенным непосредственно под ЭИ, будет придан дополнительный импульс, направленный вниз, и модель выполнит соответствующую операцию расчета движения продуктов эрозии.
Эксперименты и их обсуждение
Для исследования влияния скорости потока РЖ в МЭЗ на точность формообразования при ЭЭО была проведена серия компьютерных экспериментов. В качестве изменяемого параметра использовалась скорость V потока РЖ в МЭЗ при постоянных значениях напряжения U, подаваемого на электроды. Рабочая среда - керосин, напряжение 160 В, режимы принудительной прокачки - «прокачка РЖ вправо» и «прокачка РЖ через ЭИ». Для упрощения расчета канал в инструменте принимался бесконечно узким.
Погрешность формообразования вычислялась по формуле:
где Sфакт - количество фактически удаленного материала со стенок прошиваемого отверстия; Sном - номинальное (расчетное) количество материала, которое должно быть удалено со стенок прошиваемого отверстия.
Величина Sном рассчитывается следующим образом:
где h - высота прошиваемого отверстия; d - расчетная величина пробивного зазора.
Величина Sфакт вычисляется по имитационной модели путем подсчета ячеек фактически удаленного материала.
Результаты компьютерных экспериментов представлены в таблице, графики полученных зависимостей f = D(V) представлены на рис. 5.
Результаты компьютерных экспериментов
V, м/с
U, В |
0,005 |
0,01 |
0,016 |
0,021 |
0,042 |
0,063 |
0,084 |
0,168 |
0,336 |
Погрешность формообразования Δ,% |
|||||||||
Прокачка РЖ вправо |
|||||||||
160 |
300 |
244 |
204 |
112 |
69 |
52 |
43 |
29 |
19 |
Прокачка РЖ через ЭИ |
|||||||||
160 |
419 |
403 |
362 |
195 |
110 |
86 |
70 |
45 |
23 |
Рис. 5. Зависимость погрешности формообразования D от скорости V потока рабочей жидкости в межэлектродном зазоре: 1 - режим «прокачка РЖ вправо»; 2 - режим «прокачка РЖ через ЭИ»
Анализ графиков показывает, что зависимость погрешности формообразования D от скорости V потока РЖ носит степенной характер вида D = aV-b, где a и b - эмпирические коэффициенты. Причем степенной характер этой зависимости прослеживается очень четко - среднеквадратичное отклонение от степенного закона составляет 2,5% для режима «прокачка РЖ через ЭИ» и 5,5% для режима «прокачка РЖ вправо».
Заключение
Проведенные компьютерные эксперименты подтвердили адекватность разработанной имитационной модели процесса электроэрозионной обработки. Погрешность между экспериментальными и полученными с помощью модели результатами не превысила 6%.
В результате имитационного моделирования процесса образования единичных разрядов в межэлектродном промежутке установлено, что происходит повышенный износ электрода-инструмента на углах в связи с концентрацией в этих местах силовых линий напряженности электрического поля.
При имитационном моделировании процесса образования продуктов эрозии в результате каждого единичного разряда установлено, что суммарное действие износа электрода-инструмента и продуктов эрозии изменяет величину межэлектродного зазора и, как следствие, форму и размеры обрабатываемого отверстия.
Зависимость погрешности формообразования D при ЭЭО от скорости V потока рабочей жидкости в МЭМ носит степенной характер вида D = aV-b, где a и b - эмпирические коэффициенты.
Рецензенты:
-
Черномас В.В., д.т.н., профессор, зав. лабораторией «Новые технологии в металлургии» УРАН ИММ ДВО РАН, г. Комсомольск-на-Амуре;
-
Муравьёв В.И., д.т.н., профессор кафедры «Технологии сварочного производства» ФГБОУ ВПО «КнАГТУ», г. Комсомольск-на-Амуре.
Работа поступила в редакцию 28.08.2012.