Научная статья на тему 'Учет скоростей деформации в программном комплексе LS-DYNA® при тонколистовой вытяжке'

Учет скоростей деформации в программном комплексе LS-DYNA® при тонколистовой вытяжке Текст научной статьи по специальности «Технологии материалов»

CC BY
143
28
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Глобальная энергия
ВАК
Ключевые слова
ИМПУЛЬСНАЯ ВЫТЯЖКА / ТОНКОЛИСТОВОЙ МЕТАЛЛ / КОНЕЧНО-ЭЛЕМЕНТНЫЙ КОМПЛЕКС LS-DYNA / СКОРОСТИ ДЕФОРМАЦИИ / ДЕФОРМАЦИОННОЕ УПРОЧНЕНИЕ / IMPULSE DRAWING / THIN SHEET METALL / FINITE-ELEMENT SOFTWARE LS-DYNA / STRAIN RATE / STRAIN HARDENING

Аннотация научной статьи по технологиям материалов, автор научной работы — Арсентьева Ксения Сергеевна, Мамутов Вячеслав Сабайдинович

Исследовано влияние выбора модели материала и учета скоростей деформации при создании компьютерных моделей импульсного деформирования тонколистовых металлов. В конечно-элементном программном комплексе LS-DYNA разработаны компьютерные модели процессов вытяжки сферических деталей из тонколистового металла методом электрогидроимпульсной вытяжки. Компьютерные расчеты проводились для заготовок из латуни Л68 толщиной 0.24 мм и заготовок из алюминия 5754 толщиной 1 мм в матрицу диаметром 60 мм. Разработана методика уточнения коэффициента скоростного упрочнения металла на основе расчетов в комплексе LS-DYNA, с помощью которых определялся диапазон изменений скоростей эффективных деформации элементов заготовки. Проведено сравнение методов учета скоростного упрочнения металлов при задании различных моделей материала заготовки. Выполнена оценка применимости моделей материала заготовки для решения задач импульсного деформирования тонколистовых металлов.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по технологиям материалов , автор научной работы — Арсентьева Ксения Сергеевна, Мамутов Вячеслав Сабайдинович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Strain rate influence on thin sheet metal drawing simulation using LS-DYNA®

The article presents an investigation of the influence of material model selection and strain rate accounting on the calculation error during impulse sheet metal drawing simulation. Computer models of electro-hydraulic thin sheet metal drawing processes have been developed in the LS-DYNA finite-element complex. Computer simulation in this study was carried out for 0.24 mm-thick CuZn33 brass sheet blank and 1mm-thick Aluminum 5754 under impulse drawing into a die 60 mm in diameter. Technique of refining of speed hardening coefficient of metal using LS-DYNA complex has been developed. The technique is based on preliminary computer calculations, which allow defining a changes range of elements effective strain rates in different parts of workpiece. The technique is simple to execute and implement. Comparison of methods of the strain rate hardening accounting during usage of different material models of workpiece in computer simulations have been provided. The applicability of material models for high-speed computer simulation of thin sheet metal drawing has been assessed.

Текст научной работы на тему «Учет скоростей деформации в программном комплексе LS-DYNA® при тонколистовой вытяжке»

Машиностроение

DOI: 10.18721/JEST.261 03 УДК 621.7.044

К.С. Арсентьева, В.С. Мамутов

Санкт-Петербургский политехнический университет Петра Великого.

Санкт-Петербург, Россия

УЧЕТ СКОРОСТЕЙ ДЕФОРМАЦИИ В ПРОГРАММНОМ КОМПЛЕКСЕ LS-DYNA® ПРИ ТОНКОЛИСТОВОЙ ВЫТЯЖКЕ

Исследовано влияние выбора модели материала и учета скоростей деформации при создании компьютерных моделей импульсного деформирования тонколистовых металлов. В конечно-элементном программном комплексе LS-DYNA разработаны компьютерные модели процессов вытяжки сферических деталей из тонколистового металла методом электрогидроимпульсной вытяжки. Компьютерные расчеты проводились для заготовок из латуни Л68 толщиной 0.24 мм и заготовок из алюминия 5754 толщиной 1 мм в матрицу диаметром 60 мм. Разработана методика уточнения коэффициента скоростного упрочнения металла на основе расчетов в комплексе LS-DYNA, с помощью которых определялся диапазон изменений скоростей эффективных деформации элементов заготовки. Проведено сравнение методов учета скоростного упрочнения металлов при задании различных моделей материала заготовки. Выполнена оценка применимости моделей материала заготовки для решения задач импульсного деформирования тонколистовых металлов.

Ключевые слова: импульсная вытяжка, тонколистовой металл, конечно-элементный комплекс LS-DYNA, скорости деформации, деформационное упрочнение.

Ссылка при цитировании: Арсентьева К.С., Мамутов В.С. Учет скоростей деформации в программном комплексе LS-DYNA® при тонколистовой вытяжке // Научно-технические ведомости СПбПУ. Естественные и инженерные науки. 2020. Т. 26, № 1. С. 29-40. DOI: 10.18721/JEST.26103

Это статья открытого доступа, распространяемая по лицензии CC BY-NC 4.0 (https://creativecommons.Org/licenses/by-nc/4.0/)

X.S. Arsentyeva, V.S. Mamutov

Peter the Great St. Petersburg polytechnic university, St. Petersburg, Russia

STRAIN RATE INFLUENCE ON THIN SHEET METAL DRAWING SIMULATION USING LS-DYNA®

The article presents an investigation of the influence of material model selection and strain rate accounting on the calculation error during impulse sheet metal drawing simulation. Computer models of electro-hydraulic thin sheet metal drawing processes have been developed in the LS-DYNA finite-element complex. Computer simulation in this study was carried out for 0.24 mm-thick CuZn33 brass sheet blank and lmm-thick Aluminum 5754 under impulse drawing into a die 60 mm in diameter. Technique of refining of speed hardening coefficient of metal using LS-DYNA complex has been developed. The technique is based on preliminary computer calculations, which allow defining a changes range of elements effective strain rates in different parts of workpiece. The technique is simple to execute and implement.

Comparison of methods of the strain rate hardening accounting during usage of different material models of workpiece in computer simulations have been provided. The applicability of material models for high-speed computer simulation of thin sheet metal drawing has been assessed.

Keywords: impulse drawing, thin sheet metall, finite-element software LS-DYNA, strain rate, strain hardening

Citation: X.S. Arsentyeva, V.S. Mamutov, Strain rate influence on thin sheet metal drawing simulation using LS-DYNA®, St. Petersburg polytechnic university journal of engineering science and technology, 26 (01) (2020) 29-40. DOI: 10.18721/JEST.26103

This is an open access article under the CC BY-NC 4.0 license (https://creativecommons.org/licenses/by-nc/4.0/)

Введение. Для производства сложных точных деталей в мелкосерийном и индивидуальном производстве наиболее эффективны импульсные методы, в частности электрогидро-импульсная (ЭГИ) штамповка [1, 2]. Импульсный характер нагружения, за счет уменьшения пружинения, позволяет повысить точность изготовления деталей в сравнении с квазистатическими методами деформирования [3—6]. ЭГИ штамповка является относительно простой технологией, не требует сложной и дорогостоящей оснастки [7, 8], однако экспериментальные исследования процессов, на которых основывается данная технология, затруднены [9]. Аналитические методы исследования импульсных процессов также неэффективны, поскольку являются приближенными: не учитывают нелинейность физико-механических процессов и оперируют не всей системой исходных уравнений.

С развитием современных ЭВМ появилась возможность численного решения сложных задач обработки металлов давлением во всей полноте технологической и физической постановки. Развитие универсальных конечно-элементных (КЭ) комплексов, таких как LS-DYNA, позволяет моделировать процессы с учетом множества влияющих факторов [10]. Программный комплекс LS-DYNA имеет явный (explicit) решатель и изначально разрабатывался для расчетов воздействия подводных взрывов на металлические конструкции, т. е. рассчитан для решения задач взаимодействия

жидкостей и газов с деформируемой конструкцией и ориентирован на быстротекущие термомеханические процессы. Комплекс Ь8-БУКА успешно применяется в исследовании импульсных процессов, в частности для расчетов ЭГИ штамповки тонколистовых металлов [5, 11-14]. Программный комплекс позволяет решать многие технологические задачи, важные для проектирования, в том числе оценить влияние скоростей деформации на механические свойства заготовки.

Задачи данной работы: рассмотреть методы задания вида кривой деформационного упрочнения в расчетах; провести компьютерное моделирование в конечно-элементном комплексе Ь8-БУКА импульсной вытяжки-формовки на примере свободной ЭГИ вытяжки-формовки тонколистовой заготовки; оценить применение различных моделей материала и соответствующих кривых деформационного упрочнения в компьютерных расчетах импульсной вытяжки-формовки.

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

[15]. Для дальнейших расчетов используют аппроксимированные кривые деформационного упрочнения.

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

ст, = Вет,

(1)

1 Мамутов А.В. Технологии обработки материалов импульсным давлением. Технологии обработки концентрированными потоками энергии: учеб. пособие. Санкт-Петербург: Издательство Политехнического университета, 2008. 40 с

преобразуя степенной закон (1) в соответствующий:

ст, = КВет

(2)

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

Процесс ЭГИ вытяжки-формовки относится к импульсным процессам, при которых скорости деформации превышают 102 с1 и наблюдается дополнительное динамическое упрочнение металла по сравнению с квазистатическими испытаниями при скоростях деформации менее (0.1...1) с1 [16, 17]. Следовательно, для расчетов импульсного деформирования необходимо дополнительно учитывать скоростное упрочнение материала заготовки. В инженерных расчетах скоростное упрочнение рекомендуется учитывать с помощью коэффициента динамичности Кц, устанавливаемого для определенного диапазона значений скоростей и степеней деформации, характерных для того или иного процесса1. В рамках конечно-элементного программного комплекса Ь8-ОУМЛ коррекцию расчетов, связанную с влиянием скоростей, характерных для импульсных процессов, можно осуществить несколькими способами.

Один из способов — при задании модели материала МЛТ_Р0ШЯ_ЕЛЖ_Р1Л5Т1СПГ (МП) использовать коэффициент динамичности Кц,

Также в программном комплексе существует модель пластичности, которая учитывает влияние скоростей деформации и обычно используется для расчетов сверхпластичных процессов формоизменения: РЛТЕ_БЕМБ1Т1¥Е_ Р0ШЕК_ЬЛШ_Р1ЛБТ1С1ТУ (Л£РХ), соответствующая степенному закону:

ст = Ве е

(3)

где е — скорость деформации; п — коэффициент скоростного упрочнения.

Известно, что коэффициент скоростного упрочнения п для расчётов процессов импульсной штамповки определяется по следующей формуле [18]:

п = 921п

(е >

Vе о У

/1п е,

(4)

где 02 — экспериментальная константа матери -ала, определяемая по формуле:

02 = 1п К )/1п

(е ^

Vе0 У

(5)

Подставляя выражение (5) в выражение (4) получаем:

п=1п (Кл )/1п(е). (6)

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

Скорость деформаций при ЭГИ штамповке обычно находится в пределах от 102 до 104 с-1 из-за чего коэффициент п для одного и того же материала и вида технологического процесса варьируется в достаточно широком диапазоне (табл. 1). Большее значение коэффициента п соответствует меньшим скоростям деформации; меньшее значение коэффициента п соответствует большим скоростям деформации.

Таблица 1

Параметры коэффициента п при ЭГИ вытяжке-формовке

Table 1

Parameters of coefficient п during electro-hydraulic drawing

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

Методика уточнения коэффициента скоростного упрочнения металла. Методика уточнения коэффициента скоростного упрочнения n реализуется с помощью расчетов в КЭ комплексе LS-DYNA.

Предварительный компьютерный расчет проводится с целью определения диапазона скоростей эффективной деформации с учетом характера нагружения, свойств материала, а также формы и размеров заготовки, что позволит уточнить диапазон скоростей деформации, актуальных для конкретного исследуемого процесса. Полученные данные подставляются в выражение (6), что позволяет уточнить диапазон значений коэффициента скоростного упрочнения n с учетом специфики конкретного процесса.

Последовательность методики показана на примере расчета ЭГИ вытяжки-формовки тонколистовых металлов. Компьютерная модель для расчета создавалась следующим образом. Для заготовки, матрицы и прижима применялись оболочечные модели. Тип элемента для заготовки брался Thin Shell 163. Заготовка задавалась многослойной (21 слой) моментной оболочкой по типу Belytschko-Wong, который позволяет учитывать изгибную жесткость. Фактор изгиба (Shear Factor) задавался равным

0.833. Интегрирование напряжений по слоям осуществлялось трапецеидальным законом. Оснастка задавалась двухслойной оболочкой типа Belytschko. Материалы заготовок задавались моделями MPL и RSPL. Материалы прижима и матрицы предполагались абсолютно жесткими и задавались моделью RIGID с запретом перемещения и вращения.

Геометрия заготовок и оснастки показана на рис. 1. Начальное положение заготовки — плоскость XY, приложение нагрузки задавалось по оси Z. Для эффективного использования вычислительных ресурсов в расчете задавалась 1/4 заготовки и оснастки. Для линии, ограничивающей заготовку по плоскости симметрии ZX, устанавливался запрет на перемещение и скорость по оси Уи вращение по оси X. Для линии, ограничивающей заготовку по плоскости симметрии ZY, устанавливался запрет на перемещение и скорость по оси X и вращение по оси Y. При контрольном расчете с целыми заготовкой и оснасткой, расхождения в значениях величин скоростей деформации и эффективных пластических деформаций при прочих равных исходных данных не отмечено.

Характеристики материала заготовок даны в табл. 2. Задавалось полусухое трение, которое при решении задач не варьировалось, а его влияние на расчет не оценивалось. Параметры кривой деформационного упрочения B и m для обоих материалов были определены с помощью аппроксимации опытных данных степенным законом Холломана (1).

Форма импульса давления задавалась зависимостью:

p = poNp(t/Q)a exp(—bt/9) (7)

где p0 — амплитудное значение давления, 9 — характеристическое время, за которое давление уменьшается в e раз, t — время. Нормирующие величины Np, a, b определяются соотношениями: Np = exp(b — 1); a = bc; b = 1/(1 — с + с ln с), где с = tm/9 — коэффициент жесткости фронта давления; tm — время нарастания давления. В предварительном компьютерном расчете коэффициент жесткости фронта выбирался средним: с = 0.3.

Материал заготовки Коэффициент n

Латунь Л68 0.0242...0.0485

Алюминий 5754 0.0544...0.1087

Рис. 1. Конечно-элементная модель матрицы, заготовки и прижима для ЭГИ вытяжки Fig. 1. Finite-element model of a die, workpiece and holder for an electro-hydraulic drawing

Геометрические параметры и характеристики материала заготовки Geometric parameters and characteristics of the workpiece material

Таблица 2 Table 2

Материал заготовки Латунь Л68 Алюминий 5754

Диаметр заготовки, мм 110 124

Толщина заготовки, мм 0.24 1

Диаметр отверстия матрицы, мм 60 60

Радиус скругления кромки матрицы, мм 3 5

Модуль Юнга (Е), ГПа 115 71

Коэффициент Пуассона (у) 0.34 0.3

Плотность (р), кг/м3 8600 2660

Кулоновское трение покоя 0.15 0.15

движения 0.1 0.1

Параметры кривой деформационного упрочнения B 742 414.7

m 0.404 0.256

Коэффициент динамичности (Kd) 1.25 1.65

Амплитудная величина импульса давления подбиралась таким образом, чтобы максимальная деформация была ниже кривой БЬБ, соответствующей началу потери устойчивости (рис. 2).

Результаты расчета представлены ниже. При помощи программного комплекса Ь8-БУМЛ выводились графики изменения скорости эффективной деформации во времени —— для элементов в центре заготовки и дХ

ближайшей области (рис. 3).

Компьютерные расчеты показали, что максимальные скорости деформаций достигаются

в центральной части заготовок и снижаются в точках в зависимости от удаленности от центра. Для латуни Л68 максимальные скорости деформаций для выбранного процесса, а также геометрии оснастки и заготовки, достигают 6100 с-1; для алюминия — 3400 с-1.

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

\

\ \ * S V Si [\ ч у* / г

ч %

б)

0.8

-0.2 0 0.2 0.4

Вторая главная деформация

к

lo.e

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

S

CL

0

■а

01

г} I

ш и

Е

ш а.

1=

ч

ч к \ \ Ч

\ \ ч \ Ч \ . \ ч ч / И

4 V 1 ^ 1 V * л 4 и*6

-0.3 -0.2 -0.1 0 0.1 0.2 0.3 Вторая главная деформация

Рис. 2. Диаграмма предельных деформаций (FLD) заготовки, деформированной при ЭГИ вытяжке-формовке: а — для заготовки из латуни Л68; б — для заготовки из алюминия 5754

Fig. 2. Forming limit diagrams (FLD) of workpieces after electro-hydraulic drawing: a — CuZn33 brass;

b — aluminum 5754

Результаты компьютерного моделирования с заданием разных видов кривой деформационного упрочнения. Для того, чтобы определить возможность использования описанных выше степенных законов и моделей материала для расчетов задач высокоскоростного деформирования тонколистовой заготовки, были проведены сравнительные расчеты процесса свободной ЭГИ вытяжки-формовки заготовок из латуни Л68 и алюминия 5754. Для расчетов с моделью материала ЫРЬ использовались те же данные, что и в предварительном расчете (табл. 2). Для расчетов с моделью материала ЯБРЬ добавлялся коэффициент скоростного упрочнения п = 0.0257 для латуни Л68 и п = 0.0615 для алюминия 5754. При работе с моделью материала ЯБРЬ коэффициент динамичности Кл при задании характеристик материала не использовался. Форма импульса давления задавалась зависимостью (7). Для пар расчетов амплитудное значение давления задавалось одинаковым.

По результатам компьютерных расчетов определялась относительная погрешность рас-

чета эффективных пластических деформаций (8). Сравнение проводилось по результатам расчетов, полученных с использованием рассмотренных моделей материала ЫРЬ и ЯБРЬ.

б = 1£ ^1 ~ ^ 21 -100%, (8)

Е о Р

2

где 1 — соответствует эффективным пластическим деформациям заготовки, рассчитанные по степенному закону ЫРЬ; гред2 — соответствует эффективным пластическим деформациям заготовки, рассчитанные по степенному закону ЯБРЬ.

Эффективная пластическая деформация близка к интенсивности тензора накопленных деформаций (эквивалентной деформации). Точнее, эффективная деформация — это интеграл по времени интенсивности тензора скоростей деформаций. Оба варианта вычисления эффективной деформации совпадают, когда тензоры накопленных деформаций и скоростей деформаций соосны, что практически имеет место в условиях данных расчетов.

а)

" 7000,0

0.1R

■0.15R

-0.2R

го

о -&

ш ч:

ш m

о о

Ср

о ^

О

6000,0 5000,0 4000,0 3000,0 2000,0 1000,0 ,0

Я »• \ \

/ • \

1/ \\

\\

а / \ \\

¿у \\

б)

го

о -&

ш

:

ш m

о о

Ср g

О

4000,0 3500,0 3000,0 2500,0 2000,0 1500,0 1000,0 500,0 ,0

0,0001 0,0002 0,0003

Время, с

0,0004

0,0005

0

0.1R

■0.18R

0.25R

в • \

(Л fl 1

У \

1/

J

0,0001

0,0002

0,0003

0,0004

0,0005

Время, с

Рис. 3. График изменения скоростей деформации элементов центральной части заготовки с течением времени в зависимости от удаленности от центра (R — радиус заготовки): а) — латунь, б) — алюминий

Fig. 3. Changes of strain rates in time of workpiece elements depending on the element distance from the center (R — radius of the workpiece): a) — brass, b) — aluminum

Таблица 3

Параметры уточненного коэффициента п при ЭГИ вытяжке-формовке

Table 3

Parameters of the refined coefficient n during electro-hydraulic drawing

Материал заготовки Коэффициент n

Латунь Л68 0.0257

Алюминий 5754 0.0615

На рис. 4 показаны относительные отклонения значений эффективных пластических деформаций, рассчитанных для модели материала ИРЬ, от значений, рассчитанных для модели материала ЯБРЬ для различных элементов заготовки.

Как видно из графиков, значения эффективных пластических деформаций, рассчитанных для модели материала ИРЬ, имеют различные погрешности в зависимости от выбранного

0

0

0

элемента. Также следует отметить, что наибольшие отклонения наблюдаются в период времени до 100 мкс, который соответствует этапу начального разгона заготовки. При дальнейшем формоизменении наблюдаются отклонения в эффективных пластических деформациях менее 5 % для заготовки из латуни Л68 и менее 10 % для заготовки из алюминия 5754. К моменту остановки заготовки отличия в значениях эффективных пластических деформаций

для двух моделей материала не превышают 3 % (табл. 4), что свидетельствует о сопоставимости построенных моделей.

Также проводились сравнения по относительному профилю заготовок в момент остановки. На рис. 5 показаны кривые для заготовок из латуни Л68. Как видно из графиков, отклонения по относительному профилю не превышают 1...3 %. Для заготовок из алюминия 5754 были получены аналогичные результаты.

a)

го с

с чг

X о4

§ И

н

5 S

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Ш а. -& °

m % .о 4

3

ш

С 1_ о с

б)

fE го

5 s

ш а.

t|

О S

S S

э ¥

^ о

О g

•0

0.25R--0.5R -----0.75R---R

50% 45% 40% 35% 30% 25% 20% 15% 10% 5% 0%

100% 90% 80% 70% 60% 50% 40% 30% 20% 10% 0%

/\ , / \l // \ // \

К/Л t A

i %

к %

/ i

f \V~-_,____ ■-----

0,0001

0,0002 0,0003

Время, c

0,0004

0,0005

0.25R--0.5R -----0.75R---R

я

f \>-

K^v \ \ X

0

0,0001

0,0002 0,0003

Время, c

0,0004

0,0005

Рис. 4. Отклонения эффективных пластических деформаций для заготовки из латуни Л68 (а) и алюминия 5754 (б), рассчитанных для модели материала MPL, от значений, рассчитанных

для модели материала RSPL

Fig. 4. Deviations of effective plastic strain for brass (a) and aluminum (b) workpieces, calculated for the MPL material model, from the values calculated for the RSPL material model

0

0

Таблица 4

Погрешность эффективных пластических деформаций 8Е в момент остановки заготовки

Table 4

The Error of effective plastic strain 8E at the time of workpiece stopping

Материал 0 0,25R 0,5R 0,75R R

Латунь Л68 1,22 % 0,42% 2,71 % 0,82 % 1,47%

Алюминий 5754 0,21 % 1,06% 1,46 % 2,06 % 2,95%

Примечание: R — радиус заготовки.

0,1

0,2

0,4

0,6

0,8

о

CP

О

0 -0,1 -0,2 -0,3 -0,4 -0,5 -0,6 -0,7

Г i

i / i

1

/; 2

—'---

Относительный радиус

Рис. 5. Кривые относительного профиля заготовок из латуни Л68, полученных в момент остановки:

1 - MPL; 2 - RSPL

Fig. 5. Curves of the relative profile of brass workpieces obtained at the time of stopping: 1 — MPL; 2 — RSPL Выводы

Были рассмотрены различные методы задания кривой деформационного упрочнения в компьютерных расчетах и соответствующие им модели материала заготовки. Отмечено, что для учета скоростей деформации в модели ИРЬ используется коэффициент динамичности Кй, а в модели ЯБРЬ коэффициент скоростного упрочнения п. Разработана и представлена методика уточнения коэффициента скоростного

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

В рамках КЭ комплекса Ь8-ОУМЛ разработаны компьютерные модели процессов вытяжки-формовки тонколистовых заготовок из латуни Л68 и алюминия 5754 при импульсном нагружении для моделей материала, учитывающих различные степенные законы деформационного упрочнения, и проведены соответ-

0

1

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

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

Модель материала ЯБРЬ и соответствующая ей кривая деформационного упрочнения (3) позволяет получать более точные результаты расчетов из-за применяемого вида степенного закона, однако использование данной модели связано с некоторыми трудностями. Во-первых, для использования модели требуется предварительный расчет, уточняющий коэф-

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

Проведенное исследование позволяет сказать, что применение модели материала ИРЬ и советующей кривой деформационного упрочнения (2) с коэффициентом динамичности дает погрешность расчетов 1...3 % в сравнении с моделью материала ЯБРЬ и является приемлемым для инженерных расчетов и исследования процессов импульсной вытяжки-формовки тонколистовых металлов.

СПИСОК ЛИТЕРАТУРЫ

[1] Здор Г.Н., Исаевич Л.А., Качанов И.В. Технологии высокоскоростного деформирования материалов. Минск: БНТУ, 2010. 456 с.

[2] Linnemann M., Psyk V., Djakow E., Springer R., Landgrebe D. High-Speed Incremental Forming — New Technologies For Flexible Production Of Sheet Metal Parts // Procedia Manufacturing. 2019. Vol. 27. P. 21—26.

[3] Mehdi Zohoor, Seyed Meysam Mousavi. Experimental analysis and smoothed particle hydrodynamics modeling of electrohydraulic forming of sheet metal parts // Journal of Manufacturing Processes. 2018. Vol. 35. P. 16—28.

[4] Jia Cheng, Daniel E. Green, Sergey F. Golovash-chenko. Formability enhancement of DP600 steel sheets in electro-hydraulic die forming // Journal of Materials Processing Technology. 2017. Vol. 244. P. 178—189.

[5] Christopher Maris, Amir Hassannejadasl, Daniel E. Green, Jia Cheng, Yiteng Liang. Comparison of quasi-static and electrohydraulic free forming limits for DP600 and AA5182 sheets // Journal of Materials Processing Technology. 2016. Vol. 235. P. 206—219.

[6] Radii a Krishna Lal, Vikas Kumar Choubey, Dwivedi J.P., Shravan Kumar. Study of factors affecting Springback in Sheet Metal Forming and Deep Drawing Process // Materials Today: Proceedings. 2018. Vol. 5, is. 2. P. 1. P. 4353—4358.

[7] Юткин Л.А. Электрогидравлический эффект и его применение в промышленности. Л.: Машиностроение, Ленингр. отд-ние, 1986. 253 с.

[8] Langstâdtler L., Pegel H., Beckschwarte B., Herrmann M., Kuhfuss B. Flexible tooling for impulse forming // Procedia Manufacturing. 2019. Vol. 27. P. 130-137.

[9] Старков Н.В., Старков И.Н. Определение параметров разрядного контура электрогидравлического пресса при штамповке листовых деталей сложной формы // Электронная обработка материалов. 2017. 53(1). P. 108-113.

[10] Hallquist J.O. LS-DYNA theoretical manual. Livermore Software Technology Corporation: Liver-more, СА, 2006. 498 p.

[11] Мамутов А.В., Мамутов В.С. Моделирование электрогидроимпульсной листовой формовки с применением комплекса LS-DYNA // Научно-технические ведомости СПбГПУ. 2014. Вып. 1 (190). С. 101-107.

[12] Mamutov V.S., Mamutov A.V., Golovaschen-ko S.F. Simulation of High-Voltage Discharge Channel in Water at Electro-Hydraulic Forming Using LS-DYNA®. 1th International LS-DYNA® Users Conference, Metal forming, Dearborn, USA, 2014, p. 9.

[13] Лепеш Г.В., Иванова Е.С., Егоров К.А. Применение CAD/CAE технологий для исследования работоспособности эластичного обтюратора в условиях импульсного нагружения высоким давлением // Технико-технологические проблемы сервиса. 2016. № 1 (35). С. 24-29.

[14] Arash Jenab, Daniel E. Green, Ahmet T. Alpas, Sergey F. Golovashchenko. Experimental and numerical

analyses of formability improvement of AA5182-O sheet during electro-hydraulic forming // Journal of Materials Processing Technology. 2018. Vol. 255. P. 914-926.

[15] Третьяков А.В., Зюзин В.И. Механические свойства металлов и сплавов при обработке давлением. М.: Металлургия, 1973. 224 с.

[16] Мазуровский Б.Я., Сизев А.Н. Электрогидравлический эффект в листовой штамповке. Киев: изд-во Наукова думка, 1983. 192 с.

[17] Jitang Fan, Minqiang Jiang. Strain hardenabil-ity of a gradient metallic alloy under high-strain-rate compressive loading // Materials & Design. 2019. Vol. 170. Article 107695.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

[18] Мамутов В.С., Мамутов А.В. Теория обработки металлов давлением. Компьютерное моделирование процессов листовой штамповки. СПб.: Изд-во Политехн. ун-та, 2006. 188 с.

СВЕДЕНИЯ ОБ АВТОРАХ

АРСЕНТЬЕВА Ксения Сергеевна — кандидат технических наук ассистент Санкт-Петербургского политехнического университета Петра Великого E-mail: xenia.ars@gmail.com

МАМУТОВ Вячеслав Сабайдинович — доктор технических наук профессор Санкт-Петербургского политехнического университета Петра Великого E-mail: mamutov_vs@spbstu.ru

Дата поступления статьи в редакцию: 15.11.2019

REFERENCES

[1] G.N. Zdor, L.A. Isayevich, I.V. Kachanov,

Tekhnologii vysokoskorostnogo deformirovaniya materialov [Technologies of high-voltage deformation of materials: monograph]. Minsk: BNTU, 2010.

[2] M. Linnemann, V. Psyk, E. Djakow, R. Springer, D. Landgrebe, High-Speed Incremental Forming — New Technologies For Flexible Production Of Sheet Metal Parts, Procedia Manufacturing, 27 (2019) 21-26.

[3] Mehdi Zohoor, Seyed Meysam Mousavi. Experimental analysis and smoothed particle hydrodynamics modeling of electrohydraulic forming of sheet metal parts, Journal of Manufacturing Processes, 35 (2018) 16-28.

[4] Jia Cheng, D.E. Green, S.F. Golovashchenko, Formability enhancement of DP600 steel sheets in electro-hydraulic die forming, Journal of Materials Processing Technology, 244 (2017) 178-189.

[5] Maris Christopher, Amir Hassannejadasl, D.E. Green, Jia Cheng, Yiteng Liang, Comparison of quasi-static and electro hydraulic free forming limits for DP600 and AA5182 sheets, Journal of Materials Processing Technology, 235 (2016) 206-219.

[6] Radii a Krishna Lal, Vikas Kumar Choubey, J.P. Dwivedi, Shravan Kumar, Study of factors affecting Springback in Sheet Metal Forming and Deep Drawing Process, Materials Today: Proceedings, 5 (2-1) (2018) 4353-4358.

[7] L.A. Yutkin, Elektrogidravlicheskiy effekt i yego primeneniye v promyshlennosti [Electro-hydraulic effect and its application in industry]. L.: Mashinostroyeniye, Leningr. otd-niye, 1986.

[8] L. Langstadtler, H. Pegel, B. Beckschwarte, M. Herrmann, B. Kuhfuss, Flexible tooling for impulse forming, Procedia Manufacturing, 27 (2019) 130-137.

[9] N.V. Starkov, I.N. Starkov, Opredeleniye parametrov razryadnogo kontura elektrogidravlicheskogo pressa pri shtampovke listovykh detaley slozhnoy formy [Determination of the parameters of the discharge circuit of an electro-hydraulic press during stamping of sheet parts of complex shape], Elektronnaya obrabotka materialov, 53(1) (2017) 108-113.

[10] J.O. Hallquist, LS-DYNA theoretical manual. Livermore Software Technology Corporation: Livermore, SA, 2006.

[11]A.V. Mamutov, V.S. Mamutov, Modelirovaniye elektrogidroimpulsnoy listovoy formovki s primeneniyem kompleksa LS-DYNA, Nauchno-tekhnicheskiye vedomosti SPbGPU, 1 (190) (2014) 101-107.

[12] V.S. Mamutov, A.V. Mamutov, S.F. Golovaschenko, Simulation of High-Voltage Discharge Channel in Water at Electro-Hydraulic Forming Using LS-DYNA®, 1th International LS-DYNA® Users Conference, Metal forming, Dearborn, USA, (2014) 9.

[13] G.V. Lepesh, Ye.S. Ivanova, K.A. Yegorov,

Primeneniye CAD/CAE tekhnologiy dlya issledovaniya rabotosposobnosti elastichnogo obtyuratora v usloviyakh impulsnogo nagruzheniya vysokim davleniyem [Application of CAD/CAE technologies for research of operability of the elastic obturator in conditions of shock loading by the high pressure], Tekhniko-tekhnologicheskiye problemy servisa, 1 (35) (2016) 24—29.

[14] Arash Jenab, D.E. Green, A.T. Alpas, S.F. Golovashchenko, Experimental and numerical analyses of formability improvement of AA5182-O sheet during electro-hydraulic forming, Journal of Materials Processing Technology, 255 (2018) 914-926.

[15] A.V. Tretyakov, V.I. Zyuzin, Mekhanicheskiye svoystva metallov i splavov pri obrabotke davleniyem

[Mechanical properties of metals and alloys during pressure processing]. M.: Metallurgiya, 1973.

[16] B.Ya. Mazurovskiy, A.N. Sizev, Elektrogidravlicheskiy effekt v listovoy shtampovke [Electro hydraulic effect in Sheet Forming]. Kiyev: izd-vo Naukova dumka, 1983.

[17] Jitang Fan, Minqiang Jiang, Strain hardenability of a gradient metallic alloy under high-strain-rate compressive loading, Materials & Design, 170 (2019), Article 107695.

[18] V.S. Mamutov, A.V. Mamutov, Teoriya obrabotki metallov davleniyem. Kompyuternoye modelirovaniye protsessov listovoy shtampovki [Theory of metal forming. Computer modeling of sheet metal stamping processes]. SPb.: Izd-vo Politekhn. un-ta, 2006.

THE AUTHORS

ARSENTYEVA Xenia S. — Peter the Great St. Petersburg polytechnic university E-mail: xenia.ars@gmail.com

MAMUTOV Viacheslav S. — Peter the Great St. Petersburg polytechnic university E-mail: mamutov_vs@spbstu.ru

Received: 15.11.2019

© Санкт-Петербургский политехнический университет Петра Великого, 2020

i Надоели баннеры? Вы всегда можете отключить рекламу.