КАЗАКОВА Ольга Игоревна
студентка
КОЛПАКОВ Владимир Иванович
кандидат технических
наук, доцент кафедры «Технологии ракетно-космического машиностроения» (МГТУ им. Н.Э. Баумана)
e-mail: vorobushek11@gmail.com
УДК 519.6:532:539.5:621.9.04,011
Численное моделирование гидроабразивной резки листовых заготовок из алюминиевых сплавов
О.И. Казакова, В.И. Колпаков
Приведена физико-математическая постановка задачи процессов гидро- и гидроабразивной резки деталей из листовых заготовок алюминиевых сплавов. На конкретном примере рассмотрены особенности решения поставленной задачи на базе численных методов механики сплошной среды и зависимость исследуемого вида обработки от кинематических и физико-математических свойств инструмента (высокоскоростной струи) и обрабатываемой детали (мишени). Проанализированы возможности использования деформационного критерия разрушения материала детали из алюминиевого сплава АМг6.
Ключевые слова: математическая модель, гидрорезка, гидроабразивная резка, листовая заготовка, алюминиевый сплав, деформационные критерий разрушения, численное моделирование.
Numerical modeling of abrasive water Jet cutting of sheet billets from aluminum alloys
O.I. Kazakova, V.I. Kolpakov
The paper covers the physical and mathematical problem formulation of waterjet and abrasive waterjet cutting operations ofparts from aluminum alloys sheet billets. The features of the problem solution based on numerical methods of continuum mechanics and the dependence of this operation type on kinematic and physico-mechanical properties of the instrument (high-speed jet) and the workpiece (target) are shown by the particular example. Advantages of using the deformation criterion of the material failure (aluminum alloy AMg6) are discussed.
Keywords: mathematical model, water jet cutting, abrasive water jet cutting, sheet billet, aluminum alloy, deformation failure criterion, numerical modeling.
'Гккие методы обработки поверхностей заготовок, как точение, фрезерование, протягивание и шлифование, характеризуются тем, что режущие элементы инструмента в течение всего процесса удаления материала соприкасаются с обрабатываемыми поверхностями. Вследствие этого обработка материалов сопровождается нагревом и вибрацией инструмента и всей технологической системы в целом (станок, приспособление, заготовка), а также деформацией металла в зоне действия режущей кромки.
Необходимость обработки сложных фасонных поверхностей обусловила создание новых методов обработки, характеризующихся
отсутствием непосредственного механического контакта инструмента с заготовкой. В этих методах инструментом является либо электрическое поле (электрохимическая обработка), либо высокоскоростной направленный поток различных материалов на заготовку (гидроструйная, пескоструйная, дробеструйная обработки, обработка шариками и т. д.). Так, например, в настоящее время технологии гидрорезания (ГР) и гидроабразивного резания (ГАР) металлов, тугоплавких, жаропрочных и титановых сплавов, конструкционной керамики и композиционных материалов продемонстрировали высокую конкурентоспособность по сравнению с механическим, лазерным, плазменным и другими высокотехнологичными видами резания.
Физико-технологическую основу гидроструйных технологий составляет непрерывный высокоэнергетический поток жидкости, называемый также ультраструей, играющий роль основного формообразующего инструмента. Его режущая способность обусловлена составом и удельной кинетической энергией. ГР и ГАР имеют минимальные отличия. В первом случае (ГР) обработку материала обеспечивает водяная струя. Во втором случае (ГАР) струя разгоняет частицы абразива и эти частицы вместе со струей воды реализуют формообразующее воздействие на материал детали. У каждого вида резания свое назначение. Гидрорезкой без абразива режут мягкие материалы, а гидроабразивной резкой — твердые материалы, например, металлы, композиты и керамику, хотя в целом ГАР существенно эффективнее
ГР [1, 2].
Настоящая работа посвящена исследованию процессов ГР и ГАР на базе численных методов механики сплошной среды с использованием современного программного комплекса ANSYS AUTODYN. Цель работы — установление закономерностей рассматриваемого вида обработки от кинематических и физико-механических свойств инструмента (высокоскоростной струи) и обрабатываемой детали (мишени).
Поставленная задача решалась в двух стратегиях. В первой из них струя считалась подвиж-
ной, а мишень (преграда) — неподвижной, во второй, наоборот, мишень считалась подвижной, а струя неподвижной. При этом предполагалось, что вторая стратегия позволит более аккуратно изучить вопросы воздействия именно гидроабразивной струи на обрабатываемый материал. В процессе исследования использовались следующие материалы: для жидкой струи — вода, для преграды — алюминиевый сплав АМг6; для частиц абразива — кварцевый песок. Физико-механические характеристики материалов приведены в таблице. Диаметр струи был принят ( =1,0 мм, скорость струи V = 0,8 км/с, размер частиц абразива 0,4 мм, толщина преграды к = 5 мм.
Таблица
Физико-механические характеристики взаимодействующих материалов
Физико-математические Обозна- Материал
характеристики чение АМг6 8Ю2
Плотность, г/см3 Р0 2,64 232
Модуль сдвига, ГПа О 26,3 31
Предел текучести, ГПа У 0,15 0.15
Критическая интенсивность деформации * £ 0,43 0343
Модуль объемного сжатия, ГПа К 60 160
Идеализированные расчетные схемы ГР и ГАР показаны на рис. 1.
Поставленная задача решалась в двумерной осесимметричной постановке в переменных Эйлера. В этом случае система исходных соотношений имеет классический вид [3—5]:
Ф + ХрО + + Р^ = 0;
Ы дг дг г
(1)
(V да.
дг
+
дК
дг
2П + Л
+
(к да
& дг
- + Л + ^
дг
—; (2)
(е
РЩ = 0гг£' + 0*£гг + ^00 + 2° *£гг ; (3)
Р = р(р,е);
(4)
(5)
Р
Р
Рис. 1. Расчетные схемы ГР (а, в) и ГАР (б, г):
а, б — первая стратегия расчета; в, г — вторая стратегия расчета; 1 — струя; 2 — обрабатываемая деталь; 3 — частицы
абразива; Г1, Г2, Г3, Г4 — границы расчетной области
• • = ^ • = дУ7
£~ - , £ее , £77
дг
дУ,
д7
— г дvz
£ ^ =—~ ^--,£ ге = £ 7е = 0; (6)
д7 дг 7
Щгг
т
= 20
1 dp\ ППВ
£гг + Эр'лГ т
=0£*;
= 20
■ + .1 4Р) £ 77 + Эр & Г
(7)
/ = 2{Б1 + + + Б„Бв )< - У2. (8)
Здесь р — плотность; р — давление; е — удельная внутренняя энергия; I — текущее время; г, 7 — радиальная и осевая координаты; Vг, vz — компоненты вектора скорости; — метрические коэффициенты основного базиса выбранной системы координат, причем I, Ц= г, е, 7; о/Г, а77, аее, — нормальные напряжения в радиальном, осевом и тангенциальном направлениях; агг — касательные напряжения; Бц — компо-
ненты девиатора напряжений; — компоненты тензора скоростей деформаций; Б(...)/т — производная Яуманна; 0 — модуль сдвига; У— динамический предел текучести среды.
В приведенной системе уравнений в порядке следования представлены законы сохранения массы (1), импульса (2) и энергии (Э), в общем виде уравнения состояния взаимодействующих сред (4), взаимосвязь компонентов тензора полных напряжений с шаровой и де-виаторной составляющими (5), кинематические соотношения (6), физические соотношения в виде закона Гука в дифференциальной форме (7) и условия пластического течения Мизеса (8).
В начальный момент времени в соответствии с выбранной стратегией расчета задавалась либо скорость струи и частиц абразива, либо скорость преграды. Исходное состояние преграды соответствовало состоянию невозмущенной среды. В качестве граничных условий использовались следующие: в области
контактного взаимодействия струи с преградой накладывались ограничения на скорость индивидуальных точек в направлении оси Ог в соответствии с условиями непроницаемости материала
Vстр = Vпр г г '
а также на напряженное состояние, реализующегося в этих точках в соответствии с третьим законом Ньютона
Pn
Ир„ I
■а / п
где п' — вектор единичной нормали к поверхности преграды.
При формулировке граничных условий на оси симметрии (ось Oz) необходимо учитывать, что при r =0 частицы среды движутся только в осевом направлении (vz = 0), а осевые ускорения этих частиц должны быть ограничены. Из уравнений движения следует, что это может быть реализовано только при отсутствии касательных напряжений на оси симметрии (arz = = 0) [4]. При решении исходной системы уравнений на неподвижной сетке область интегрирования ограничена (см. рис. 1): снизу осью симметрии (Г1); слева, сверху, справа — открытыми поверхностями (Г2, Г3, Г4), через которые среда может вытекать или втекать.
Уравнения состояния для металлов использовали в форме линейной (linear) баро-тропной зависимости давления от плотности p = К(р/р0 — 1). Здесь р0 — начальная плот-
ность. Другие физико-механические характеристики материалов представлены в таблице. Уравнение состояния воды принималось в виде полиномиальной (polynomial) зависимости давления от плотности:
p = Л1ц + Л2ц2 + Л3ц3 + (B0 + Вгц)р0е при ц < 0; p = Тхц + Т2ц2 + B0 -р0e при ц <0,
где ц = р/р0 — 1; р, р0 = 1 г/см3 — текущая и начальная плотности соответственно. Коэффициенты уравнения состояния воды принимали следующие числовые значения: А1 D= 2,2 ГПа, А2 = 9,54 ГПа, А3 =14,57 ГПа, В0 = В1 = 0,28, Т1 = 2,2 ГПа, Т2 = 0.
В качестве критерия разрушения материала преграды использовался деформационный критерий вида [6]
£ р < £ р ,
где £,
£ P£ P — интенсивность пластиче-
ских деформаций; £ Р — пластические составляющие компонент тензора деформаций; £*Р — предельная (критическая) пластическая деформация.
На рисунках 2—5 представлены этапы ожидаемых механизмов разрушения алюминиевой преграды (АМг6) толщиной к = 5 мм, полученные расчетным путем, с использованием деформационного критерия разрушения. Вне зависимости от используемой стратегии, скорость струи или преграды составляла ¥0 = = 0,8 км/с, диаметр частиц абразива 0,4 мкм.
Рис. 2. Этапы проникновения высокоскоростной струи воды в алюминиевую преграду по первой стратегии
расчета
Рис. 3. Этапы проникновения высокоскоростной струи воды и частиц абразива в алюминиевую преграду
по первой стратегии расчета
Рис. 4. Этапы проникновения высокоскоростной струи воды в алюминиевую преграду по второй стратегии
расчета
Рис. 5. Этапы проникновения высокоскоростной струи воды и частиц абразива в алюминиевую преграду
по второй стратегии расчета
Рис. 6. Начальное и последующие конфигурации частиц абразива (позиция 1)
Сравнивая представленные рисунки необходимо отметить следующее. Несмотря на качественное совпадение регистрируемых процессов, наблюдается их количественное разночтение по времени. Это связанно со спецификой используемых расчетных стратегий и требует дополнительного объяснения. Кроме того, наблюдается более высокая скорость резки ГАР по сравнению с ГР. Отмеченная особенность наиболее ярко иллюстрируется при использовании второй расчетной стратегии, что, очевидно, связано с эффектом «размазывания» частиц абразива на эйлеровой сетке (рис. 6).
Выводы
Разработаны физико-математическая модель и стратегии расчета ГР или ГАР, основанные на использовании идельноупругопластической модели взаимодействующих сред и прямого или обращенного движения струи по отношению к обрабатываемой детали в зависимости от используемого метода обработки.
Апробированы критические условия разрушения материала детали из алюминиевого сплава АМг6 в условиях высокоскоростного воздействия гидроабразивной струи.
Показаны возможности разработанного подхода и определенные преимущества ГАР по отношению к ГР в плане скорости обработки поверхности.
Необходимо отметить, что представленная работа является лишь начальным этапом для решения более сложных задач, связанных с оптимизацией технологических параметров ГАР, в частности на характер поведения частиц абразива в фокусирующем сопле.
Литература
1. Ультраструйная технология получения микросуспензий / О.Е. Балашов, А.А. Барзов, А.Л. Галиновский и др. М.: Изд-во МГТУ им. Н.Э. Баумана, 2011. 352 с.
2. Тарасов В.А., Полухин А.Н.Оценка геометрических параметров формируемой поверхности при гидроабразивной резке // Вестник МГТУ. Сер. Машиностроение. 2012. № 1 (86). С. 107—116.
3. Физика взрыва / Под ред. Л.П. Орленко. В 2 т. Т. 2. М.: ФИЗМАТЛИТ, 2004. 656 с.
4. Бабкин А.В., Колпаков В.И., Охитин В.Н., Селиванов В.В. Численные методы в задачах физики быстропротекающих процессов. В 3 т. Т 3. М.: Изд-во МГТУ им. Н.Э. Баумана, 2006. 520 с.
5. Теоретические и экспериментальные исследования высокоскоростного взаимодействия тел / под ред. А.В. Герасимова. Томск: Изд-во Томского университета, 2007. 572 с.
6. Селиванов В.В. Механика разрушения деформируемого тела. М.: Изд-во МГТУ им. Н.Э. Баумана, 1999. 420 с.
Статья поступила в редакцию 02.05.2012