научная статья по теме МОДЕЛИРОВАНИЕ ПРОЦЕССА ДЕФОРМАЦИИ МЕТАЛЛА, ИМЕЮЩЕГО ДЕФЕКТЫ, НА ЛИТЕЙНО-КОВОЧНОМ МОДУЛЕ Машиностроение

Текст научной статьи на тему «МОДЕЛИРОВАНИЕ ПРОЦЕССА ДЕФОРМАЦИИ МЕТАЛЛА, ИМЕЮЩЕГО ДЕФЕКТЫ, НА ЛИТЕЙНО-КОВОЧНОМ МОДУЛЕ»

ПРОБЛЕМЫ МАШИНОСТРОЕНИЯ И НАДЕЖНОСТИ МАШИН

№ 2, 2014

УДК 539.3; 628

© 2014 г. Одиноков В.И., Бондаренко С.В.

МОДЕЛИРОВАНИЕ ПРОЦЕССА ДЕФОРМАЦИИ МЕТАЛЛА, ИМЕЮЩЕГО ДЕФЕКТЫ, НА ЛИТЕЙНО-КОВОЧНОМ МОДУЛЕ

Процесс получения заготовок на литейно-ковочном модуле (ЛКМ) включает заливку жидкого металла, его кристаллизацию и деформацию в заданный профиль (рис. 1). По данному процессу в Институте машиноведения и металлургии ДВО РАН были опубликованы основные результаты экспериментальных и теоретических исследований [1—4]. В зоне кристаллизующегося металла имеют место различного рода дефекты в виде пор, раковин и т.д. (рис. 2). Очень важно знать их эволюцию в процессе деформации, чтобы найти режимы работы ЛКМ для их ликвидации. С этой целью необходимо построение соответствующей математической модели.

Математическая постановка задачи. Рассмотрим процесс деформации на ЛКМ вертикального типа (рис. 1). Схема деформации металла с учетом поверхностей симметрии представлена на рис. 3. Дефект в металле моделируется наличием полости определенных геометрических размеров. Весь поворот эксцентрикового вала разобьем на т шагов. Будем исследовать процесс деформации на угле поворота 180°. Величина одно-

т

го шага Р„ = 180°/т при ^ Рт = 180°.

- т

1

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

= о, ^ = а= <0• ';.»,

I о, I *.,

1 г % _ л (1)

3

0 а = 1 = 0,5 ( и,. + I),

1 * = Т = Т(н, Е, 9), Н = (2%.)1/2,

с у = а1у (1 grad 9), (2)

дт

где г = 1, 2, 3 (суммирование по повторяющимся индексам г,у); а,у — компоненты тензора напряжения; — компоненты тензора скоростей пластических деформаций; иг — проекция скоростей перемещений по координатным осям (г = 1, 2, 3); 1 — коэффициент теплопроводности; 1 = 1(9); 9 — температура; с — удельная теплоемкость; у — удельный вес; функция Т = Т(Н, Е, 9) определяется из эксперимента; Н — интенсивность скоростей деформаций сдвига; Е — степень деформации.

%

V/VQ 1,0

120 а

Рис. 1

Рис. 4

Рис. 1. Схема литейно-ковочного модуля вертикального типа: 1 — эксцентриковые приводные валы, 2 — за-ливочно-дозирующее устройство, 3 — боковая плита, 4 — боек, 5 — эксцентриковый приводной вал, 6 — заготовка

Рис. 4. Изменение объема пустот в процессе поворота эксцентрикового вала на угол а

В области, где металл находится в жидком состоянии, имеем ньютоновскую жидкость

Ъц - = 2(3)

где р — коэффициент вязкости. Зоны твердого и жидкого металла определяются по температуре: 9 > 9кр — жидкий металл; 9 < 9кр — твердый металл, где 9кр — температура кристаллизации.

Начальные условия: т = 0 ^ а0 = 0; E0 = 0; 9 = 90 (4), где а0 — угол поворота эксцентрикового вала; E0 — начальная степень деформации металла; 90 — начальное поле температур при а0 = 0 находится из предположения, что поток металла движется стационарно через раскрытый зев наклонных бойков [4].

Граничные условия (рис. 3). Полагаем, что граница исследуемой области описывается системой ортогональных поверхностей, тогда имеем

'12

3)|5. = 0 ' = 2> 7; (Ъ22 = Ъ21 = Ъ23) = 0, СТ23|= 0, ' = 4> 5> 6;

(ъ21 )|8, = , . = 4 5> 6; (ъ21 = Ъ23)^ = 0; (ъ31 = Ъ32) |58 = 0

( иск )9.

(5)

(Ъ32 )| 59 = 0; (ъ31 )| 59 = -^9т5 —;— ; (Ъ

12

3 )| 5 = 0, I = 10, 11;

(ъ22 = Ъ21 = Ъ23) |5. = 0, ' = 12, 13; (Ъ33 = Ъ31 = Ъ33 )| 5. = 0, ' =

= 14, 15;

1

= 0; и215 = -(оп)., I = 4, 5, 6; 03!^ = 0, . = 8, 9;

9= 9., . = 2, ..., 7, 9; д\8 = 0, ; = 1, 8.

Рис. 2. Дефекты в залитом в ЛКМ металле после его охлаждения в вертикальной (рис. 1) плоскости симметрии

I У

Рис. 3

Здесь т5 — предел текучести материала при сдвиге; и — нормирующая скорость, (иск), = (им — ии), — скорость скольжения металла (им) относительно скорости инструмента деформации (ии) на ,-й поверхности контакта; (и„) , — скорость перемещения инструмента по нормали на ,-й поверхности контакта; 9,- — известные из эксперимента температуры на внешних поверхностях 5,; — тепловые потоки через поверхности S¡, являющиеся поверхностями симметрии; у, — коэффициенты трения на поверхностях 5,; скорость перемещения (и„) , ( , = 4, 5, 6) определяется кинематикой движения бойков.

Для решения уравнений (1)—(3) при наличии граничных условий (4), (5) использовали численные методы, описанные в [5]. Суть методов состоит в том, что исследуемая область разбивается на ортогональные криволинейные элементы, для каждого элемента записываются в разностной форме уравнения (1)—(3), которые решаются по разработанным численным схемам [4, 5] при наличии начальных (4) и граничных условий (5), (6). Результатом решения являются поля напряжений, скоростей переме-

Рис. 5

щений и температур в исследуемой области на каждом временном шаге. При решении задачи использовался пакет программ [6, 7].

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

Параметры, определяющие геометрию системы (рис. 3): к0 = 15 мм; кх = 4 мм; Ь = 30 мм; у = 8°; R1 = 50 мм; /1 = 45 мм; /2 = 60 мм; e = 3 мм; ep = 10 мм, где e — величина эксцентрика привода наклонных бойков, ep — величина эксцентрика привода вертикальных плит, п — скорость вращения эксцентриковых валов.

При теоретических расчетах модель деформируемой среды T = T(H, E, 9) в системе (1)

построена путем аппроксимации экспериментальных данных [8] с учетом T = aS / 43 ,

H = ^л/3 , Г = Е^Ь , что справедливо для простого нагружения, при котором проводили эксперимент

Т = 1,72 -(Г)0'36 • Я0'09. (7)

Коэффициенты в (7) найдены методом наименьших квадратов. Погрешность аппроксимации не превышает 5%.

Коэффициент трения у, = 1000, , = 4, 5, 6, что соответствует практически полному прилипанию металла к инструменту, = 0 ^ ст31= 0.

Размеры пустот а1 х а2 х а3 = 3 х 2 х 7 мм задавали кратными размерам наносимой расчетной сетки. Начальные координаты пустот: 1 — х1 = 55,5 мм; х2 = 0, х3 = 0; 2 — х1 = 70,7, х2 = 6,8, х3 = 22,5; 3 — х: = 64,6, х2 = 11,7, х3 = 0; 4 — х: = 64,6, х2 = 0, х3 = 0.

На рис. 4 приведены теоретические результаты поведения дефектов (пустот) при деформации металла на ЛКМ вертикального типа. При этом V — текущий объем раковины (пустоты); V) — начальный объем раковины. Из рис. 4 видно, что пустоты наиболее интенсивно закрываются в слое (2), близком к поверхности наклонных бойков (рис. 3, поверхность 54). Кроме того, следует отметить, что результаты на рис. 4 приведены за один оборот. Каждый дефект (1, ..., 4) проходит очаг деформации, по крайней мере, за три поворота эксцентрикового вала. Следовательно, все пустоты заданных размеров (особенно у поверхности наклонных бойков) должны полностью закрыться. С этой целью был проведен эксперимент, в ходе которого была осуществлена деформация трех свинцовых полос с круглыми отверстиями диаметром 5 мм (рис. 5, а). После осуществления деформации видно, что отверстия в центральной пластине значительно уменьшились (рис. 5, б), а в пластине, примыкающей к поверхности — полностью исчезли (рис. 5, в).

Это объясняется тем, что в зоне интенсивной деформации (под поверхностями S4, S5, S6) напряжения стп, ст22, ст33 сжимающие, а гидростатическое давление ст в приграничных к бойкам слоях выше (по абсолютной величине), чем в остальном объеме.

На рис. 6 приведены графики сти, ст22, ст — в плоскости симметрии х3 = 0 при деформации полосы без дефектов. Видим, что в приконтактном к бойкам слое нормальные

Рис. 6. Графики 022, Пц, а в плоскости симметрии x3 = 0: 1 — на оси симметрии; 2 — в слое, приконтактном к наклонным бойкам

напряжения сти и ст выше (по абсолютной величине), чем в слое при x2 = 0. Это свидетельство правильности выбранного режима — прилипания металла к бойкам на поверхности Si = 0, i = 4, 5, 6.

На рис. 7 приведены графики ст22, ст в плоскости x:x2, в которой расположен дефект № 1 при повороте эксцентрикового вала на угол 120°.

Работа выполнена при финансовой поддержке РФФИ (Гранты 11-01-98500 и 12-1-0ЭММПУ-06).

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

1. Одинокое В.И., Стулов В.В. Литейно-ковочный модуль. Литье и деформация. Владивосток: Дальнаука, 1998. 150 с.

2. Одинокое В.И., Проскуряков Б.И., Черномас В.В. Непрерывный процесс кристаллизации металла при одновременном его деформировании. М.: Наука, 2006. 111 с.

3. Одинокое В.И., Проскуряков Б.И., Черномас В.В. и др. Повышение надежности литейно-ко-вочного модуля // Проблемы машиностроения и надежности машин. 2009. № 2. С. 66—69.

4. Одинокое В.И., Ловизин Н.С., Скляр С.Ю. Моделирование процесса деформации металла на литейно-ковочном модуле // Математическое моделирование. 2010. Т. 22. № 9. С. 129—145.

5. Одиноков В.И., Каплунов Б.Г., Песков А.В. и др. Математическое моделирование сложных технологических процессов. М.: Наука, 2008. 176 с.

6. Одиноков В.И., Соснин А.А. Математическое моделирование процессов деформации металла на литейно-ковочном модуле горизонтального типа. Свидетельство о государственной регистрации программ для ЭВМ № 2011615520 от 14.07.2011.

7. Одиноков В.И., Соснин А.А., Саликов С.Р. Комплексное моделирование тепловых и деформационных процессов, протекающих в кристаллизаторе установки литья и деформации металла при получении металлоизделий. Свидетельство о государственной регистрации программ для ЭВМ № 2012613947 от 10.01.2012.

8. Кроха В.А. Упрочнение металлов при холодной пластической деформации. Справочник. М.: Машиностроение, 1980. 157 с.

Комсомольск-на-Амуре Поступила в редакцию 4.III.2013

Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.

Показать целиком