Методика моделирования процесса штамповки металла в твердожидком состоянии Текст научной статьи по специальности «Компьютерные и информационные науки»

НАУКА и ОБРАЗОВАНИЕ

Эя №<К! 77 - 30569. Государственная регистрация №0421100025.155Г11ЭЭ4-040в_

Методика моделирования процесса штамповки металла в твердожидком состоянии

# 08, август 2011 автор: Шитиков А. А.

УДК 621.7

МГТУ им Н.Э. Баумана [email protected]

1. Реологическая модель материала в твердожидком состоянии

В настоящее время исследователями предложены различные варианты реологических моделей материала в твердожидком состоянии. Все предложенные модели можно разделить на две большие группы:

- однофазные модели,

- двухфазные модели.

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

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

Агломерация - процесс образования агломератов (плотных групп) из разрозненных твердых частиц.

Деагломерация - процесс, обратный агломерации, т.е. процесс разделения группы на отдельные частицы.

Наглядно влияние этих процессов на изменение вязкости суспензии изображено на рис. 1

Рисунок 1. Влияние степени агломерации частиц s на вязкость суспензии

В двухфазных моделях рассматриваются отдельно твердая и жидкая фазы со своими законами движения, и их влияние друг на друга.

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

К недостаткам двухфазных моделей можно отнести:

- необходимость проведения большого числа экспериментальных исследований для каждого материала,

- для определения параметров модели требуются специальные реометры,

- при моделировании требуется вычисление двух полей скоростей, что приводит к значительному увеличению времени расчета. По этой причине пока отсутствуют программные пакеты производящие расчет в трехмерном пространстве.

1.1. Базовая реологическая модель

В качестве базовой для программной реализации выбрана однофазная модель, приведенная в [1]

°Е = аЬ. + °р ' Кв

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

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

- коэффициент, учитывающий вязкоупругие свойства суспензии

1.1.1. Гидродинамическое течение твердых частиц

Для описания течения твердых частиц, разделенных жидкостью в работе [1] предлагается использовать экспериментальный закон, предложенный Кригером для густых суспензий

где 1|0 - вязкость жидкости

- доля твердой фазы в суспензии

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

- показатель степени. Для глобулярных не компактных кластеров к = 2.6

Для Ньютоновских жидкостей зависимость сопротивления сдвигу от вязкости определяется законом

Т = Т] ■ у

Для использования в программных комплексах, предназначенных для моделирования процессов обработки давлением, данная зависимость должна быть приведена к форме а = /(б)

Соотношение между вязкостью жидкости при сдвиге и вязкости при удлинении определяется коэффициентом Троутона (Trouton ratio)

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

Для модели суспензии с полностью агломерированной структурой [1]

где s - коэффициент агломерации (см. п. .3) Окончательно определим как

1.1.2. Оценка влияния процесса агломерации на напряжение текучести

Для течения суспензии требуется полное либо частичное разрушение связей в агломератах твердых частиц. При этом предел текучести увеличивается на необходимое для этого напряжение &р

где Уь - объем деформируемого металла в каждой связи. Полагается [1], что эффективный объем - сфера с радиусом т.е.

п* = ■ пг

Для несжимаемых Ньютоновских жидкостей Тг = 3 [3] Учитывая, что; получим

ар ^Ь ' ' ®НТ

- объемная концентрация связей, может быть определена как общий объем пар соединенных частиц радиуса г [1]

Г _ 3 & • 5

ъ 8 ■ п г3

Су- - предел текучести твердого металла при температурах, близких к температуре плавления, может быть принят по одной из множества зависимостей, используемых для горячей объемной штамповки. Окончательно

3

Ч г

3 Гь 8 п г* ант~ 2 Лг) °нт

1.1.3. Определение степени агломерации частиц в суспензии

Изменение в структуре тиксотропной твердожидкой суспензии могут быть описаны как совокупность двух процессов:

- самопроизвольного упрочнения (загустевания) покоящейся суспензии (т.е. ее агломерации)

- разупрочнения (размягчения) суспензии при сдвиге (т.е. деагломерации)

Т.о. изменение коэффициента агломерации 5 - во времени можно представить как

где <Л (_£, Т) - функция, описывающая процесс агломерации,

Г■: г) - функция, описывающая процесс деагломерации При 5 = 0 все частицы разделены жидкостью, 5 = 1- все частицы образовали прочные связи между собой

Функция агломерации t, Г)

Агломерация - процесс образования устойчивых связей между частицами.

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

Перечислим основные параметры, определяющие процесс агломерации:

- /;- частота столкновений частиц, пропорциональна скорости деформации сдвига у

- tc - время контакта частиц при столкновении,

- (1 — 5) объемная доля свободных частиц.

- Таким образом, возможность агломерации

Для определения частоты столкновений предложена формула [1]

_8 с тс Среднее время контакта [1]

я

С У

Для образования прочной связи между двумя частицами необходима определенная взаимная ориентация их кристаллических решеток. В [1] указывается, что появление прочной связи возможно только в случае меньшеугловых границ (< 0.25 рад). Т.о. вероятность столкновения двух частиц с благоприятным взаимным положением кристаллических решеток составляет только [1]

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

В [1] приведена зависимость, для определения времени £у образования шейки радиусом ть

г2 -ть- (С±-

К

ь

где г - радиус частиц,

- тангенс угла наклона линии ликвидуса к оси концентрации на

двухкомпонентной диаграмме состояния. В общем случае т1 = —,

с1С

Сь и С5 - концентрация легирующего элемента в жидкости и твердом растворе соответственно,

- коэффициент диффузии атомов легирующего элемента, в жидкости,

Л-4 - изменение энтропии на единицу объема [2],

КТ(_С5-С,)

■ Кг (1 - О.) ■ С1

- молярный объем,

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

г

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

п

±г 8 ■ /г ■ у у

«Л =Р*-Л-^ ■(!-*) =0.02-------^---(1-*)

^ 7Г

fs

я = 0.16 (1-я)

Таким образом, <Я- ф /(у)

Функция деагломерации Ї)

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

_ ^СГІҐ СГІҐ ^

Другими словами, за время происходит изменение степени

агломерации (в начальный момент времени £0) до 5 = 0. Изменение

степени агломерации Дз за промежуток времени Ат = Г — ї0

—і? Д£

или

Аз є

Окончательно, функцию деагломерации V можно определить как

2> ' Вт ^ ^

єсги

Учитывая полученные выражения для функций <А и 2), имеем

. .. Л . . €

— = 5 = Л{я^г Т I + Т){зг£г £1 = 0.16 — ■ (1 — ,1)--------■ ^

гг £ега

Решая полученное дифференциальное уравнение, найдем

СО— ( а ) ([(( ] о + а

V 0 а — в.' а — сI

где

0.16

а

1.1.4. Вязкоупругие свойства суспензии

Предположим, что вязкоупругое поведение тиксотропной суспензии описывается моделью Максвелла [3] механический эквивалент которой приведен на рис. 2

Рисунок 2. Модель вязкоупругого тела Максвелла

а а .

Е і]

Решение данного дифференциального уравнения

)]

о

где 1М = - -постоянная времени Максвелла, определяет скорость

Е

уменьшения напряжений в зафиксированной, предварительно растянутой системе (см. рис. )

модуль упругости

Введем обозначение

Ке= 1

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

£ ■ 1

/ї '

+ 0.1-УЛ

Л? )

-2,5

ИТ

ехр

и

\£ Тмп

1.2. Расчетная реологическая модель

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

фазовый состав для условий равновесных и неравновесных фазовых переходов.

Примерами распространенных программ являются:

- MTDATA (www.mtdatasoftware.tech.officelive.com) National Physical

Laboratory (UK),

- ThermoCalc (www.thermocalc.com),

- JMatPro (www.sentesoftware.co.uk).

В данной работе для предсказания свойств материала использована демонстрационная версия программы JMatPro.

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

1) свойства, которые могут быть предсказаны в программе JMatPro:

2) параметры, которые можно определить из двухкомпонентной (псевдо двухкомпонентной) диаграммы состояния: ть, Сь, СБ.

3) параметры, которые необходимо определить экспериментально:

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

л, Л Ми-о.1*гЛ

Множитель ^1-----—-J определяет разницу между

вязкостью жидкой фазы ij0, и вязкостью тиксотропной суспензии. Однако

в программе JMatPro возможно получить значение параметра Total

Viscosity - общая вязкость (rj-p) в диапазоне температур Ts — TL, т.е. можно упростить выражение для исключив из него экспериментальный параметр f*\

В работе [1] предлагается использовать модель Зинера-Холомона

где (Т0 - предел текучести при температуре солидуса

г - = 1 ■_ - - необходима для соблюдения размерности

- энергия активации самодиффузии

- универсальная газовая постоянная

- показатель степени, зависит от свойств материала

Данные о значении параметров 0 и т приведены в литературе только для ограниченного числа легкоплавких сплавов. При моделировании штамповки из сталей, №, Co и Ti сплавов можно предположить, что

где оНТБ(Т) - параметр High Temperature Strength полученный в JMatPro, который показывает зависимость предела текучести материала во всем интервале температур вплоть до TL, при фиксированной скорости деформации 0.01 c-1 и определяется в программе по формуле

<*нг$(Т) = ес+ р

где <ж, /?, определяются по регрессионным зависимостям, полученным при обработке экспериментальных данных

Точность определения оНТ5(Т) для жаропрочных № сплавов можно оценить по графику приведенному в [6] и изображенному на рис. 3

Рисунок 3. Сравнение предела текучести № сплавов, вычисленного в JMatPro с результатами

экспериментов

При подобной замене теряется зависимость предела текучести от деформации и скорости деформации, однако, в случае использования полной версии программы ТМаРго возможно получение кривых ■7 = г т ) в широком диапазоне температур.

В [9] (п. 2.13 стр. 81) описывается схожая реологическая модель

,Т71

— Г°о

.т ■

О

Там же рекомендуется принять

Л ■ ТУ

— = о.;

1-

— £

М'

Как видно из формулы множитель —, в ней используется не в третьей

т

(как в базовой модели), а в первой степени.

Для использования в базовой реологической модели, примем с целью сохранения значения множителя

Параметр £сг1г в моделировании принят равным 0.4 . Т.к. в [1] указывается, при данной относительной деформации течение суспензии приобретает установившийся характер и твердые частицы в жидкости движутся раздельно

Сформируем на основании перечисленных допущений расчетную модель материала

Полный перечень параметров, определяемых в программе ТММРго, приведен в таблице

Таблица 1. Параметры, которые возможно определить с помощью программы JMatPro

Назначение Обозначение в формуле Обозначение в JMatPro

Параметры реологической модели fs Fraction solid (wt)

Vt Total viscosity

D Liquid diffusivity

Molar volume

Y^i Surface tension

E Young's modulus

Параметры, необходимые для решения тепловой задачи Heat capacity

Л Thermal conductivity

Предел текучести при высокой температуре для сталей, N0 ^ и ^ сплавов *ht = m High temperature strength

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

Рисунок 4. Схематизированная диаграмма состояния

С0 - равновесная концентрация суммы легирующих элементов в сплаве

С; С^ - концентрации суммы легирующих элементов в твердой и

жидкой фазе соответственно при температуре Т*

1 температура плавления чистого растворителя Как видно из рис. 7,

= = т5=Хл=

С0 С0

_ ^плА~^ _^~плА

ы £

При использовании диаграмм состояния другого вида эти выражения соответственно изменяются

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

1.3. Особенности программной реализации

Для реализации расчетной модели использованы пользовательские процедуры (User Routines) программного комплекса Deform. Изменение исходного кода процедур, и их последующая компиляция и компоновка, позволяют получить новый конечно-элементный процессор DEF_SIM.exe для расчета на однопроцессорных машинах и DEF_SIM_P4.exe и DEF_SIM_P4P.exe для параллельного расчета на многопроцессорных машинах.

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

Пользовательские процедуры делятся на процессорные и постпроцессорные.

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

Процессорные пользовательские процедуры могут использоваться

для:

- учета изменения плотности материала при штамповке порошков,

- вычисления и учета изменения объема материала,

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

- вычисления процессорных пользовательских переменных (в отличие от постпроцессорных их можно использовать в пользовательской модели материала),

- определения пользовательской модели движения инструмента (может использоваться для моделирования штамповки на нестандартном оборудовании),

- задания пользовательского критерия оценки вероятности разрушения материала,

- задания пользовательской модели износа штампов,

- задания пользовательской модели фазовых превращений в материале,

- других специфических нужд, при использовании процедуры ИБЯМБН.

Исходный текст пользовательских процедур написан на языке БОЯТКАК и для их использования необходим компилятор АЬбой БоЛгап у7.0 или более поздней версии.

При программировании рассмотренной модели материала использованы процедуры:

- ИБЯМТЯ - для определения предела текучести,

- ИБЯИРБ - для вычисления пользовательских переменных.

1.3.1. Процедура и8ЯМТЯ.

Выходные параметры процедуры:

- УБ - напряжение текучести,

- УРБ - частная производная напряжения текучести по накопленной деформации,

- Б1Р - частная производная напряжения текучести по эквивалентной скорости деформации.

Параметры УРБ и Б1Р необходимы для последующих переразбиений сетки КЭ в процессе расчета.

Расчетная модель используется в интервале температур Т5 — Ть

YPS =

до*

де

де

FIP =

да%

34

ds

■Vt +ду<Глг+5- д.

-+ S - и.

Lf2 -Т

■ 0Хр|

где сгнт = f (Є, Є,.

дх

де

~и

.(-+И-)

L4 ClT-rrf7

Параметры

U

v

и

вычисляются с помощью численного

Зе Эё

дифференцирования.

При температуре ниже используется модель материала

Данные для определения напряжений текучести &ИТ и aLT, представляют собой трехмерные массивы с максимальной размерностью (1:9, 1:7, 1:7), размещенные в текстовых файлах, при этом для определения cLZ возможно использование файлов из базы материалов DEFORM. Экстенты массива соответствуют степени деформации, скорости деформации и температуре соответственно.

Для вычисления ет при промежуточных значениях Т, Е, ё применена интерполяция:

- по £, Т - линейная,

- по е - логарифмическая, с основанием логарифма 10.

1.3.2. Процедура и8ЯиРБ

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

сетки КЭ в процессе расчета. Соответственно различают узловые и элементные пользовательские переменные.

Элементные переменные вычисляются в начале каждой итерации, узловые - только в конце сошедшегося шага, поэтому данном случае они не могут использоваться в процедуре ИБИМТЯ.

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

Таблица 2. Вычисляемые элементные пользовательские переменные

Номер переменной Обозначение в программе Обозначение в модели Обозначение в постпроцессоре

Значение на предыдущем шаге Значение, на текущем шаге

1 USRE1(1) USRE2(1) fs Solid fraction

2 USRE1(2) USRE2(2) s Agglomeration degree

3 USRE1(3) USRE2(3) Hydrodynamic flow stress

4 USRE1(4) USRE2(4) Sigma P

Номер переменной Обозначение в программе Обозначение в Обозначение в

Значение на предыдущем шаге Значение, на текущем шаге модели постпроцессоре

5 USRE1(5) USRE2(5) кв Viscoelastic coeff.

6 USRE1(6) USRE2(6) (аъ T5*T*TL \&LT> T<TS Total yield stress

7 USRE1(7) USRE2(7) T Temperature

8 USRE1(8) USRE2(8) de YPS

9 USRE1(9) USRE2(9) doj dr FIP

10 USRE1(10) USRE2(10) - Cooling rate

Номер переменной Обозначение в программе Обозначение в модели Обозначение в постпроцессоре

Значение на предыдущем шаге Значение, на текущем шаге

11 ШКШ(11) ШКБ2(11) <JHT Sigma HT

12 ШКБ1(12) ШКБ2(12) &1.Т Sigma LT

Исходные данные для вычисления 1, 2, 3 и 5 переменных являются одномерные массивы с фиксированной размерностью в 60 элементов,

представляющие собой значения функций (Т±)~~, —“тт.. V 3 - 7] т(Т;)

где Ъ = Т5 + Ть Тз - (г — 1Xг = 1,2 ...60. Данные массивы последовательно

расположены в файле SSV.txt

Значение температуры на текущем шаге для каждого элемента в программе определено как среднее арифметическое температур узлов этого элемента (учитываются только узлы с ненулевой температурой).

Для -1 шага в препроцессоре должны быть заданы начальные значения переменных и имена, под которыми они будут обозначены в постпроцессоре см. рис. 5.

Рисунок 5. Задание начальных значений элементных переменных в препроцессоре

Также должны быть заданы 7^, Ть, путь к файлам с исходными данными и название файла с данными для определения аьт см. рис. 6

Рисунок 6. Исходные данные

Все текстовые файлы с исходными данными могут быть сформированы в среде МаШсаё

2. Пример моделирования штамповки в твердожидком состоянии

На рис. 8 приведен пример моделирования штамповки ступицы тормозного диска из литейного алюминиевого сплава А356. Чертеж поковки и заготовки, схема оснастки, а также параметры технологического процесса приведены в [9].

При моделировании принято

г = 50 мкм

>-0.004б9'7’ . ^0.022 . ^0.15

,зЧ1-0.154191

О-

1Т'

- - - г- [10]

[1]

Графическая интерпретация зависимостей aLT(е,е,

и

ант(Ё,

приведена на рис. 7

Sigma LT

а)

б)

Рисунок 7. Графическая интерпретация зависемостей

А356

(ег€,Т) (а) и ент(ёгТ) (б) для сплава

Рисунок 8. Моделирование штамповки в твердожидком состоянии ступицы тормозного диска из сплава А356 (распределение твердой фазы в момент окончания заполнения штампа

и график силы деформирования)

%

36.377

29.102

21.826

14.551

7.275

0.000

Histogram

—1—1—1—1— Min: 0.498 1

Мах: 1.00

- - Аод:-0.854

5.D.: 0.158 1 1 1 1 1 1

1 г 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 1

— 1 I 1 1 I 1 1 I 1 1 I 1 1 I 1 1 1 1 шщ 1 1 1 И1

0.498

0.593

0.699

0.799

0.900

1.000

fraction

Рисунок 9. Гистограмма распределения доли твердой фазы в объеме поковки в момент

окончания заполнения штампа

Литература

1. A. Wahlen. Processing of aluminum alloys in the semi-solid state. Basic principles and constitutive models. Diss. ETH Nr. 14166, 2001

2. M. Gunduz, J.D. Hunt. The measurement of solid-liquid surface energies in the Al-Cu, Al-Si and Pb-Sn systems. - Acta metall. Vol. 33, No. 9, pp. 1651-1672, 1985

3. R.P. Chhabra, J.F. Richardson. Non-Newtonian flow in the process industries. Fundamentals and engineering applications. - Butterworth-Heinemann, 1999

4. Thixoforming. Semi-solid metal processing. Edited by G. Hirt and R. Kopp. - WILEY-VCH Verlag GmbH&Co. KGaA, Weinheim, 2009

5. Z. Guo, N. Saunders, A.P. Miodownik, J-Ph. Schille. Material properties for process simulation. Materials Science and Engineering A 499 (2009) 7-13

6. Z. Guo, N. Saunders, A.P. Miodownik, J-Ph. Schille. Quantification of high temperature strength of nikel-based superalloys, Material Science Forum, Vols. 546-549, 2007

7. М. Флемингс. Процессы затвердевания. пер. с англ. - М.: Изд-во «МИР», 1977.

8. Немнюгин М.А., Стесик О.Л. Современный Фортран. Самоучитель. -СПб.: БХВ - Петербург, 2004.

9. Семенов Б.И., Куштаров К.М. Производство изделий из металла в твердожидком состоянии. Новые промышленные технологии : учеб. пособие. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2010.

10.Хензель А., Шпиттель Т. Расчет энергосиловых параметров в процессах обработки металлов давлением: Справ. изд. Пер. с нем. М.: Металлургия, 1982.