CC BY
85
Физика твердого тела
УДК 539.214: 621.762 А.Ф. Федотов
МОДЕЛИРОВАНИЕ И ИССЛЕДОВАНИЕ ПРОЦЕССА ДЕФОРМИРОВАНИЯ ГОРЯЧИХ ПРОДУКТОВ СВС В ОБОЛОЧКЕ ИЗ СЫПУЧЕГО МАТЕРИАЛА
Выполнено моделирование и исследование процесса деформирования горячих продуктов СВС в оболочке конечных размеров из сыпучего материала с учетом температурного режима и структурного состояния деформируемого материала. Адекватность моделей реальному физическому процессу подтверждена экспериментально. Установлено ведущее влияние структурного состояния на закономерности процесса уплотнения и формоизменения.
Среди многообразных схем горячего деформирования порошковых материалов находит применение прессование в закрытой матрице с передачей давления через сыпучую среду. Эта схема получила название квазиизостатического прессования. Квазиизостатическое прессование особенно широко применяется при получении компактных материалов из продуктов саморас-пространяющегося высокотемпературного синтеза (СВС) [1,2]. В этом случае сыпучая оболочка (песок) обеспечивает теплоизоляцию саморазогревающегося материала от холодного деформирующего инструмента и после полного сгорания шихтовой заготовки деформируется совместно с продуктами СВС. Научный и практический интерес представляет изучение закономерностей уплотнения и формоизменения синтезируемого материала и энергосиловых параметров процесса. Такую информацию можно получить при аналитическом исследовании процесса совместного деформирования продуктов синтеза и оболочки. В известных моделях [3-5] рассматривается квазиизостатическое прессование в бесконечной оболочке с заданием статических граничных условий через коэффициент неизостатичности (коэффициент бокового давления) ^ 5 отражающий реологические свойства оболочки. Эти модели не учитывают размеры и внешнее трение оболочки и не могут быть использованы для оценки количественных связей между силовыми параметрами процесса и конечными свойствами прессуемого материала.
Цель настоящей работы - моделирование процесса совместного деформирования пористого нелинейно-вязкого включения и сыпучей оболочки в замкнутом объеме конечных размеров. Реологической моделью уплотняемого включения служат горячие продукты СВС. Рассматривается осесимметричное течение при условии однородности температурного поля и изотермичное™ процесса деформирования. Моделирование кинематических и силовых параметров процесса выполняется методом верхней оценки мощности пластической деформации.
Силовому компактированию подвергаются горячие продукты синтеза, саморазогреваю-щиеся до температуры горения Тг. Технологическая задержка во времени, необходимая для сгорания всего объема шихтовой заготовки, приводит к охлаждению продуктов синтеза до температуры компактирования Тк. Соответственно исследование механического поведения продуктов СВС должно сочетаться с оценкой температуры Тк. Теплофизические и реологические свойства пористых гетерогенных материалов существенным образом зависит от параметров структурного состояния - пористости и фазового состава. Поэтому при моделировании процессов теплообмена и деформирования в первую очередь должен быть выполнен анализ макроструктурного состояния продуктов СВС.
Исходные порошковые компоненты формуются в пористую шихтовую заготовку, для которой известны высота кш , радиус , абсолютная <ЛШ и относительная рш плотности. Параметры макроструктуры продуктов СВС формируются в результате протекания трех процессов, сопровождающихся изменением объема вещества и пористого тела:
1) изменение плотности и объема конденсированной фазы в результате химического превращения;
2) выделение значительного количества газов, которые были адсорбированы или растворены в исходных компонентах;
3) твердо- и жидкофазное спекание продуктов СВС.
Примем следующую модель формирования начальной относительной плотности р/о и размеров к10, Ящ включения. При горении шихтовой заготовки, которое протекает при незначительном давлении подпрессовки, за счет внутреннего давления примесного газовыделения сохраняются размеры шихтовой заготовки. Процессом уплотнения при спекании вследствие быстротечности синтеза можно пренебречь. Тогда включение из продуктов СВС будет наследовать размеры шихтовой заготовки: кю = кш и Я к, — Кш. Начальная относительная плотность при этом будет равна:
Мш
р.»(1)
' шиСВС
где Мш , V ш - масса и объем шихтовой заготовки; й свс - теоретическая абсолютная плотность беспористых продуктов синтеза. В дальнейшем в качестве уплотняемого включения будем рассматривать продукты синтеза системы 77 - С - Ш, которыми являются карбид титана 77С и никелевая связка. Система 77 - С - № является базовой при синтезе разнообразных твердосплавных материалов на основе карбида титана. Размеры шихтовой заготовки массой 180 г. и расчетные значения начальной относительной плотности шихтовой заготовки и продуктов синтеза приведены в табл. 1.
Таблица1
Параметры плотности и размеры шихтовой заготовки и продуктов СВС
Содержание А7, % мае.
/гш, мм
Яш, мм
(1Ш, г/см
< свс,
г/см
Р ю
20
14,0
39
4,11
0,65
5,4
0,5
28
13,4
39
4,34
0,65
5,62
0,5
35
13,0
39
4,57
0,65
5,83
0,5
Продукты синтеза при температуре горения представляют собой дисперсную среду, в которой часто присутствует расплав эвтектики или легкоплавкого компонента. При расчете физических и теплофизических свойств гетерогенного материала должны учитываться пористость, состав и количество жидкой и твердой фаз. В табл. 2 приведены температуры горения исследуемых составов [6]. Величина Тг для всех составов превосходит эвтектическую температуру системы 77С - М, составляющую 1523 К [7], и продукты синтеза состоят из твердых частиц карбида 77С и карбидоникелевого расплава. Абсолютная плотность двухкомпонентных смесей (карбидо-никелевый расплав) и б?сдс (твердожидкое тело) рассчитывается по известной зависимости
(2)
Х1^2 5
где а, , <13 - плотность компонентов; лг/ - массовая концентрация 1-го компонента. Массовая концентрация жидкой фазы л:/, и никеля в расплаве х№ определяется по диаграмме состояния 77С - М [7] для соответствующей температуры.
Коэффициент теплопроводности и удельная теплоемкость Сь расплава определяются по правилу аддитивности:
^ I ~ ^А7ХЛ7 + 0 ~ ХЛ7 ) 5 СI. ~СтХш + (1 — Хт ) . (3)
Теплопроводность Хк и теплоемкость Ск компактного двухкомпонентного твердожидкого материала рассчитываются по зависимостям [8]
= ^пс(^ _ Хь) ~ 0>72|л.£ — — (4)
ск = с,,ті + Сж (/ - т1 У + 4ть (/ - ть) С^С”с ,
С/, + Сж
где гп1 - объемная доля жидкой фазы в компактном материале
т,
с1жХц,
(5)
(6)
При расчете теплофизических свойств продуктов СВС использовали справочные данные [9,10]: 3ж - 4,92 г/см3, Хпс — 46 Вт/м-К; (1щ = 8,9 г/см3; Хт = 69 Вт/м-К; С м — 735 Дж/м*К. Теплоемкость карбида 77С при температуре Тг рассчитывается по зависимости [11]
Сж = 825 + 55,8 -10~3Т -25 • 10б Т. (7)
Теплоемкость жидкого карбида 77С при Тш = 3340 А* составляет СГгСж>= 1002 Дж/кг-К. Результаты расчета теплофизических свойств сплавов с различным содержанием Ш представлены в табл. 2.
Т аблица 2
Теплофизические свойства компактных продуктов СВС
Содержание №, % мае. тг,к Хщ Ч ть Вт/м-К К, Вт/м-К Сі, Дж/кг-К г Дж/кг-К
20 2670 0,55 0,36 0,297 54,5 48,4 855,6 964
28 2510 0,61 0,46 0,382 56,1 49,6 839,5 905,8
35 2350 0,66 0,53 0,443 57,5 50,8 826 885,7
Плотность, теплопроводность и теплоемкость пористого материала включения 1 определяются по зависимостям
4=4*с0-в) 'Лх =4(1-1,50); С1=Ск(1-0) + С,0, (8)
где Су - удельная теплоемкость воздуха, равная Су = 1005 Дж/м-К; 0 - пористость включения, связанная с относительной плотностью ру зависимостью 0 = 1- р/. Теплофизические свойства пористых продуктов СВС при относительной плотности р = 0,5 представлены в табл. 3.
Таблица 3
Теплофизические свойства пористых продуктов СВС (р = 0,5)
Содержа- нием, % мае. иг, мм/с 1,0 г/см3 А,/, Вт/м-К с,, Дж/кг-К тк,к хм ті Р 50
20 15 3,1 2,7 12,1 967,5 2540 0,59 0,34 0,27 0,365
28 10 4,4 2,81 12,4 955,4 2360 0,65 0,43 0,34 0,33
35 6 7,0 2,91 12,7 945,0 2180 0,71 0,49 0,39 0,305
Для песка соответствующие свойства равны ^ ==195 г/см3; Х2 = 0,326 Вт/м-К; С2 = 795 Дж/кг-К [10].
Инициирование реакции горения производится с боковой поверхности шихтовой заготовки, имеющей форму круглой пластины. Неизохронность нагрева движущимся фронтом горения обусловливает разное время охлаждения отдельных зон включения 1. Соответственно к моменту прессования в продуктах СВС формируется неоднородное поле температур Тк и реологических свойств. Используемая схематизация теплового режима предполагает однородное распределение свойств по всему объему включения. Поэтому в дальнейшем рассматриваются осред-ненные температурно-временные характеристики, относящиеся к центральной зоне изделия. Время охлаждения этой зоны t принималось равным времени охлаждения центральной зоны круглой пластины радиусом /?/ = 39 мм:
К,
t = -L+t, и,
3 9
(9)
где tз - время задержки прессования, равное времени срабатывания исполнительных органов пресса. Для гидравлических прессов время составляет ~ 0,5 с. Скорость горения иг была определена экспериментально. Оценка температуры Тк в центральной зоне изделия проводилась из решения известной задачи об охлаждении пластины толщиной /г, в неограниченной среде с граничными условиями 4-го рода [12]:
4 1 ^ (2п -1 )к -2х (2п -1)/?, + 2х
ег/с-----77==-----+ еф-
Г- Г
Т -Т
* г с
= 1-
(10)
4у[а^ 4д/«7
где Тс - температура среды; Ке - критерий тепловой активности пластины; а/ - коэффициент температуропроводности пластаны; / - время охлаждения. Коэффициенты К£ и а/ связаны с 94
коэффициентом теплопроводности X , удельной теплоемкостью С и плотностью <2 соотношениями
к.=.1.'„1! ; А,=тт. (11)
ОД
2000
Здесь индекс 1 относится к свойствам продуктов СВС, индекс 2 - к свойствам песчаной оболочки. При тепловом расчете теплофизические свойства считались постоянными величинами. На рис. 1 представлены результаты расчета распределения температуры по высоте включения с А/ = 14 мм для сплавов с различным содержанием М. Видно, что, несмотря на существенное различие начальных температур включения и среды, распределение температуры близко к линейному и характеризуется незначительным градиентом. За температуру Тк принималось среднее по высоте значение температуры (см. табл. 3). Расчетные значения Тк для всех сплавов превышают эвтектическую температуру системы ТгС - М. составляющую 1523 К [7]. Поэтому продукты горения при деформировании находятся в твердожидком состоянии. Сопротивление сдвигу в жидкой фазе на несколько порядков ниже, чем в твердой фазе. В приближении идеальной жидкости сопротивление вязкому течению оказывает только твердая фаза, а жидкую фазу можно
рассматривать как несжимаемую пору [13]. Тогда для пористого твердожидкого тела мерой уплотнения будет относительная объемная доля твердой фазы р5, которая связана с относительной плотностью р соотношением [13]
0 1 2 3 4 5 6 7х,мм
Рис. 1. Распределение температуры по высоте включения: 1 - 20%М;
2 - 28%М; 3 - 35%М; х - расстояние от центра
р5 =р(1 -ть),
(12)
где ть - объемная доля жидкой фазы в компактном материале. С точки зрения реологии параметр р5 является эффективной плотностью твердожидкого тела. Результаты расчета начальной величины рда в сплавах с различным содержанием № представлены в табл. 3.
В нормальных условиях при отсутствии разрыхляющего действия газовыделения и начальная плотность продуктов СВС не может быть меньше насыпной плотности дисперсной среды. Насыпная плотность является структурно-чувствительной характеристикой и зависит от формы и размера частиц, состояния их поверхности, фракционного состава и плотности самого вещества. Теоретические значения насыпной плотности могут быть оценены исходя из представлений о случайных упаковках и координационных числах N. В присутствии расплава твердожидкие частицы продуктов СВС имеют форму сфер примерно одинакового диаметра. Для монодисперсных сферических частиц со случайной упаковкой при насыпной плотности координационное число N = 8, а насыпная относительная плотность составляет по данным различных авторов р„ = 0,53 - 0,6. Лучшее соответствие с экспериментом имеет значение р„ =0,53 [14]. Согласно расчетам начальная эффективная плотность ряо в продуктах СВС меньше насыпной плотности, и следует ожидать самоуплотнения материала включения. Экспериментально эффект самоуплотнения в продуктах синтеза системы 71 - С - № не наблюдается, и это означает, что за счет внутреннего давления примесных газов поддерживается взвешенное состояние твердой фазы. Начальное аэродисперсное состояние включения должно учитываться при моделировании процесса прессования.
При исследовании процесса прессования весь деформируемый осесимметричный объем разобьем на 4 блока (рис. 2): включение 1 (продукты СВС), центральный 2 и два кольцевых 3, 4 блока оболочки. В реальном объекте включение 1 со всех сторон окружено оболочкой. Поэтому в блоки 2, 4 объединены объемы оболочки, расположенные сверху и снизу включения. Скорость перемещения пуансона v0 считается известной. На границах раздела блоков скорости течения непрерывны, но их значения неизвестны. Кинематически допустимое поле скоростей
осесимметричного течения зададим из условия однородности объемной деформации е в соответствующей области:
е = ez + er + е = const,
(13)
где е2,ег,ер - компоненты тензора скоростей деформаций:
dv2 dvr vr
6z ~~dr' £r~T '
(14)
Здесь vz,vr - осевая и радиальная скорости пластического течения материала. Положим, что скорость vz зависит только
от координаты z , а скорость vr только от координаты г. Тогда условие (13) можно записать следующим образом:
Ш!
Ж
Ш
R.
R
Рис. 2. Схематизация процесса прессования
dv2
dz
= С- ^. + ^ = С
dr г 3’
(15)
где С},Сз- постоянные. Решение уравнений (15) имеет вид
у2(г) = С!г + С2; уг(г) = С3г + С4г~‘. (16)
Постоянные интегрирования находятся из кинематических краевых условий для каждого блока:
vj0) = 0; II •**» N vn(0)
vz2(fy) = vy; vz3{h) = v0; vr 2(0)
vz3(0) = 0; vz3(h) = vo\
vj0) = 0; vl4(b) = v0; vJr,)
11> Уг4\Км) = 0- (1^)
После вычисления постоянных интегрирования в (16) по краевым условиям (17) определяются компоненты тензора скоростей деформаций в каждой области деформируемого объема:
hl
Vn
Uj
— V,
**2 =■
h-hx
er2 ~~ e<p2 ~
er3=~
er4=~
u2Rj
R
м -R,
u,R,
1 + -
R
M
Rm
R,
i + -
R
г2 ,
2Л
e„, = -
u2R,
Rm2
■R,‘
l—
R
M
M
r2 у
e,„. = —
UjR,
r2 j
R
м
R,
l-
R
2\
M
Г2 у
(18)
где к, к] - текущие высоты всего деформируемого объема и включения; Л;, Ям - радиусы включения и матрицы. Текущее значение радиуса Я2 блока 2 равно радиусу включения Я] . Нетрудно убедиться, что при принятом поле скоростей скорости объемной деформации блоков не зависят от координат. Это означает, что плотность в каждом блоке во все моменты времени распределяется равномерно.
Соотношения (16) и вытекающие из них уравнения (18) справедливы в случае вязкого сжимаемого включения. При аэродисперсном состоянии вязкость продуктов СВС близка к нулю и они не оказывают сопротивления пластическому течению. Соответственно можно положить, что на границах включения и оболочки напряжения будут равны нулю:
' z 1
'г 3
О ;
0.
(19)
Из статических граничных условий (19) вытекают следующие соотношения между осевыми и радиальными компонентами скоростей течения
для блока 2
2 \хщ (к-к)
для блока 4
№оКх
«і = ~Г7------\---Т- (21)
Л[Б(1-ц)-ц]
Здесь ц - коэффициент Пуассона [15],
Зр-1
"=1^1’ (22)
В - коэффициент, учитывающий размеры включения и матрицы,
Я/ + Км
В=я1^- (23)
В итоге на начальном этапе прессования, когда продукты СВС находятся в аэродисперсном состоянии, виртуальное поле скоростей определяется только одним независимым параметром -радиальной скоростью и2.
При уплотнении включения до беспористого состояния оно становится несжимаемым:
= 0. При несжимаемом включении скорости V; и и1 связаны между собой соотношением
у1Я1
(24)
и поле скоростей определяется двумя свободными параметрами - скоростями V/ и и2 .
Действительное поле скоростей будем искать энергетическим методом. Равенство мощностей внешних и внутренних сил на виртуальном поле скоростей при заданной скорости перемещения пуансона Vo имеет вид
Ягд=1Гпл+И'трб+!Гтрп, (25)
где И'д - мощность усилия прессования (2; - суммарная мощность пластической деформа-
ции включения и оболочки; УУтрб , УУтрп - мощность удельных сил трения скольжения на боковой поверхности матрицы и пуансоне. В уравнение (25) не вошли мощность сил трения на границах включения и оболочки, а также мощность на границах блоков с разрывами скоростей. Отсутствие этих факторов диссипации энергии обусловлено следующим. Температура горения выше температуры плавления песка, и песчаная оболочка, находясь в контакте с горячими продуктами СВС, плавится. Расплав песка служит смазкой, и поэтому силы трения скольжения на границах включение-оболочка можно не учитывать. СВС-прессование производится на гидравлических прессах с невысокой скоростью перемещения деформирующего инструмента, и течение материала близко к статическому. Согласно [16] в уплотняемых пластических телах, свойства которых зависят от плотности, разрывы скоростей при квазистатическом течении невозможны.
Компоненты уравнения (25) определяются следующим образом.
Мощность усилия прессования ():
(26)
я»
где 8„ - площадь пуансона;
суммарная мощность пластической деформации
4
^=£1> (27)
(=/ у,
где Д - удельная мощность пластической деформации во включении (г = 1) и оболочке (г =2, 3, 4); ¥1 - объемы соответствующих областей;
мощность удельных сил трения скольжения на боковой поверхности матрицы площадью
&
^трб~ \DmpdS , (28)
где Отр - удельная мощность сил трения оболочки;
1Утрп = 2
мощность удельных сил трения скольжения на пуансоне
/ \
где 32, - площадь контактных поверхностей блоков 2 и 3. Множитель 2 в (29) учитывает
контактное трение оболочки на дне матрицы.
При постоянной скорости перемещения пуансона мощность прессования
Щ=<2у0. (30)
Тогда действительное поле скоростей и неизвестные параметры V; ,щ ,и2 находятся при минимизации следующего функционала:
' ( Л =Е ЩЛУ + +г
,=/ Г, 5, 45, 53
Удельная мощность пластического течения определяется реологическими свойствами уплотняемых сред. Горячее прессование продуктов СВС происходит при гомологических температурах не менее 0,7 Тт , и процесс течения твердой фазы контролируется механизмами ползучести. В общем случае течение вещества кристаллической природы осуществляется по механизму нелинейной ползучести:
(31)
У
8 = Аоп,
(32)
где 8 - скорость ползучести; <7 - напряжение; п - показатель нелинейности. Для тугоплавких соединений, в том числе и карбида 77С, п = 3. Температурная зависимость процесса течения определяется коэффициентом А :
А=~^ехр{-~й)’ <33> где Т - термодинамическая температура; II - энергия активации; к - постоянная Больцмана; С -константа. В известных работах [16-18] рассматривается мощность диссипации пористого линейно-вязкого тела. Найдем эту характеристику для пористого нелинейно-вязкого тела.
Определяющая зависимость между тензором напряжений сг,у и тензором скоростей деформаций е,, для пористого тела с произвольным законом течения твердой фазы имеет вид [17, 18]
а(и>)
(34)
где - эквивалентная скорость деформации твердой фазы; <з(м>) - эквивалентное напряжение; е - скорость изменения объема; 5„ - символ Кронекера; \|/, ф - функции относительной плотности р
у
Ф = Р •
(35)
3 1-р
Эквивалентная скорость деформации твердой фазы выражается через инварианты макроскопического тензора скоростей деформаций:
™ ~ ~г^Ц>е2 + ФН2 , (36)
л/ Р
где Н - интенсивность скоростей деформаций сдвига. Выражение (36) получено в [18] исходя из двух постулатов:
1) пористая среда диссипирует механическую энергию только за счет пластической деформации твердой фазы;
2) в процессе пластической деформации участвует и диссипирует энергию весь объем твердой фазы.
Теория, построенная на этой физической модели, адекватно описывает процессы пластического и вязкого течения малопористых спеченных материалов со связанной структурой твердой фазы. Песчаная оболочка и продукты СВС относятся к несвязанным дисперсным средам (порошки). Механизм деформирования дисперсных сред имеет две особенности. Во-первых, кроме пластического течения, происходит скольжение контактирующих частиц, и материал
уплотняется также за счет их переупаковки. При низкой плотности материала преобладающим механизмом уплотнения является скольжение, при высокой плотности - пластическое деформирование частиц. Во-вторых, деформация отдельных частиц дисперсной среды имеет неоднородный характер, и пластически деформируется не весь объем твердой фазы, а лишь его часть [19]. Соответственно определяющие уравнения теории пластичности дисперсных сред должны отражать двойственный механизм уплотнения и неоднородный характер пластической деформации частиц.
Продукты синтеза содержат жидкую фазу в виде расплава легкоплавких компонентов. В приближении идеальной жидкости расплав выполняет роль смазки, и диссипируемой при скольжении частиц энергией можно пренебречь. Согласно [19] объемная доля пластически деформируемой (или активной) твердой фазы
2 Р — Ро
ш = р --------, (37)
1 - Ро
где р0 - насыпная плотность порошкового материала. Соотношение (37) получено в [19] из вероятностного описания образования контактных связей между частицами порошка. Аналогичное выражение получено в [20] из модели уплотнения за счет роста радиусов сферических частиц вокруг фиксированных центров с дальнейшим перекрытием сфер и образованием шеек контакта.
Для песка в [21] получено соотношение для эффективного объема активной фазы, которое интегрально учитывает диссипацию энергии внутренним трением и неоднородный характер деформации частиц:
© = 1 - ехр
1,3
' 0,91 — р0Л 1,4 1п
(38)
0,91-р
Насыпная плотность песка, просеянного через сито с диаметром ячейки 1 мм, составляет р0 = 0,61. В итоге эквивалентная скорость деформации твердой фазы в дисперсных средах
и> = ——^уе2 +фН2 . (39)
д/СОр
Информация о механизме течения твердой фазы содержится в зависимости между скоростью деформации м> и напряжением &(■»>). Для линейно-вязкого тела эта зависимость имеет вид [18]:
ст(мЛ
м> = -^, (40)
2Ло
где т|0 - коэффициент сдвиговой вязкости твердой фазы; для нелинейно-вязкого тела
^ = |^,ст"(н'). (41)
Коэффициент 112 в (41) выбран из условия предельного перехода при п =1, когда ■ Ко-
эффициент А1 аналогично (33) определяет температурную зависимость течения. Мощность пластической деформации
В = а у еу ■ (42)
Из (34) следует связь между инвариантами тензоров напряжений и тензоров скоростей деформаций:
Ст(м0 а(м0 гг
с =-------и/е , Т—фН. (43)
т 2м?
После подстановки (34) в (42) с учетом (43) получим выражение мощности пластического течения для пористого нелинейно-вязкого тела (включение 1):
I
1 ( 2 V — —
д 2п (2ц>е2 + ф#2)2” . (44)
2 \А1у
Материал оболочки (несвязанный грунт) будем считать жесткопластическим телом [22]. Для жесткопластического тела напряжение а(м/) соответствует удвоенному пределу текучести при сдвиге х5 пластичных или пределу прочности хь хрупких частиц. В нашем случае частицы
песка являются хрупкими и
с(м>) = 2хь. (45)
Согласно [13] для частиц песка хь = 115 МПа. Используя (45) и (43), получим скорость диссипации механической энергии в пористом жесткопластическом теле (оболочка):
£> = л/сор-^уе2 +фЯ2 . (46)
Мощность активных сил, затрачиваемая на преодоление сил трения скольжения, равна
А»Р = , (47)
где хк - касательные напряжения на контактной поверхности. Величину х*. примем в соответствии с законом трения Прандтля [16]:
13
(48)
Скорости скольжения в соответствующих областях пуансона представляют собой радиальные компоненты скоростей течения
V
и2К,
м
а на боковой поверхности матрицы - осевую компоненту скорости
скб
и
(49)
(50)
После определения параметров действительного поля скоростей оценивалось напряженное состояние и усилие прессования (верхняя оценка). Компоненты тензора напряжений при нелинейно-вязком течении включения 1 рассчитывали по зависимости (34) с учетом соотношения (41):
1-и 2 п
С 2 \у> (_____2 . ~гг2Л
=
2\|/е + фЯ ©р
У-\ч)еЪц +Ф^
Для сыпучей оболочки после подстановки соотношения (45) в (34) будем иметь:
ст„ =
2х6Л/©р
л]2ц/е2 + фЯ2
1
(51)
(52)
Усилие прессования <2, находили из условия статического равновесия для каждого шага нагружения. В отличие от состояния пластического течения при статическом равновесии вместо закона трения Прандтля используется закон трения Кулона. Величина О равна осевому давлению на пуансоне, суммированному по его площади. При наличии внешнего трения осевое давление распределено по высоте матрицы неравномерно: максимальное давление действует на пуансоне и снижается по мере приближения к дну матрицы. Согласно [15] давление р на пуансоне примем равным
I А
р= К +/,
(53)
м
где /тр - коэффициент трения скольжения Кулона, который зависит от относительной плотности песка [21]:
= 0.63р
-0,73
(54)
Усилие £? находится суммированием давления р по площади пуансона. При этом следует учесть, что пуансон испытывает давление со стороны двух блоков - 2 и 3. Кроме того, давление на стенку матрицы оказывают также два блока 3 и 4 с радиальными напряжениями аг3 и аг4. В качестве расчетного напряжения а г примем средневзвешенное радиальное напряжение:
+ (55)
Усилие прессование будет равно:
б= Лоя|.в + Л<т„^ + ■ <5б>
Слагаемые в правой части уравнения (55) можно рассматривать как усилия, затрачиваемые соответственно на деформирование включения и центральной зоны оболочки > кольцевой ЗОНЫ оболочки Оз и на преодоление СИЛЫ трения оболочки О стенки матрицы <2дд-
е = е,+е,+&,- (57)
С увеличением плотности песка давление на стенку матрицы неограниченно возрастает и вместе с ним возрастает удельная сила трения х тр . Однако ее значение не может превысить значения, допускаемого условием текучести. Учитывая, что напряжения а г, ст ^, а 2, определяемые формулами (52), являются главными, максимальное касательное напряжение трения Прандтля будет равно [16]:
1| I
^тах = 21СТ*3-°г|- (58)
В итоге третье слагаемое в (56), представляющее собой усилие трения Qmp , расходуемое на преодоление сил внешнего трения оболочки, определяется выражением
£>тр = К, к И > (59>
где Т«р =/»рк| ПРИ /«рк1<Тш«5 Т»Р =Ттах ПРИ /шр|<*г| ^ ^тах •
Реологические свойства продуктов СВС определяются термокинетическими параметрами нелинейной ползучести карбида титана: показателем нелинейности п , предэкспонентом С и энергией активации и. Экспериментальные значения этих параметров получены в [3] при горячем изотермическом прессовании порошка карбида ТгС в жесткой пресс-форме с учетом сил внешнего трения. При выводе аналитического уравнения использовалась объемная вязкость линейного тела
р3
С = Чо-г-, (60)
1-р
где г)о - сдвиговая вязкость вещества твердой фазы. В нашем случае значение функции у (35) приводит к выражению
4 р3
<61>
з 1 р
т.е. отличается от используемого в [3] коэффициентом 4/3. Соответственно в определяющих соотношениях (44) и (51) для нелинейно-вязкого тела коэффициент А / будет связан с коэффициентом Л из (33) зависимостью
А,=[^2 А. (62)
Для карбида титана параметры уравнения (33) составляют: показатель нелинейности и = 3; предэкспонент С = 910'34 м7-Н2/с; энергия активации 11= 8,42-10"19 Дж.
Полученные реологические зависимости (44), (46), (47) в совокупности с параметрами структурного состояния материала и кинематическими соотношениями пластического течения определяют функционал (31) с точностью до свободных параметров V/, щ и щ .
Моделирование процесса квазиизостатического прессования выполняли методом ступенчатого нагружения. При этом задавали шаг изменения времени А/ , за который происходит
осадка деформируемого объема на величину ДА = у0Д/. При выборе Д/ ориентировались на технологические параметры реального процесса: скорость движения пуансона у0 = 8 мм/с; максимальное усилие прессования = 1500 кН; время прессования / = 1 с. Пробными расчетами установлено, что достаточная точность и производительность расчетов обеспечивается при Д? =0,02 с. В каждом интервале Д/ величины, зависящие от времени (относительная плотность р, высоты /г, /г/ и радиусы Л/, Я2 = /?/), считаются постоянными. Минимум функционала (31) находили методом координатного спуска с погрешностью вычисления скоростей 0,001 мм/с. Численное интегрирование при вычислении мощности и компонентов тензора напряжений в блоках 3, 4 выполняли по методу Симпсона. После минимизации функционала (31) и определения параметров V], щ, и2 вычисляются приращения осевых и радиальных размеров блоков: ДА; = У/Д/, ДК/ = м/ДГ, Д/?2 = и2Д?. Для каждого к-тото шага деформации ДА по формулам
км— А;*.; + Акй К-1к — + АЛ; ; Я2к = К2ы + АЛ2 определяются новые значения размеров бло-
ков и плотность материала в них. Плотность включения находится из закона сохранения массы:
Ь\к-\
(63)
\ к ^ 1
В начале каждого шага нагружения радиус блока 2 равен радиусу включения и граница между блоками 2 и 3 является проницаемой. При изменении радиусов Л; и Я2 на А/?; и ДЯ2 массопе-ренос через границу 2-3 составляет:
ДМ = лр(/г - \ )(я,2 - Л2 2 ). (64)
При Я]к > Н2к перетекание материала происходит из блока 3 в блок 2 и плотность р в выражении (64) равна р = рЗАг; при Я1к < Я2к поток массы направлен в противоположную сторону и р = р2к. Масса ДМ алгебраически суммируется с массами блоков и рассчитывается средняя по объему плотность материала в блоке 2 и объединенных блоках 3, 4. Такое осреднение приводит к скачкообразному изменению плотности на границе 2-3.
На рис. 3 приведены результаты расчета высоты включения А;, относительной плотности Р; и экспериментальные данные о толщине средней части синтезированных изделий в зависимости от среднего давления компактирования рК, равного отношению усилия прессования Q к площади матрицы (диаметр матрицы составлял 125 мм, начальная высота деформируемого объема А - 35 мм). Хорошее соответствие расчетных и экспериментальных значений получено для всех трех СВС-композиций. Это свидетельствует об адекватности принятых физических моделей и связанных с ними теплофизических и реологических свойств продуктов синтеза, а также достоверности модели пластического течения песчаной оболочки. Уплотнение продуктов СВС в песчаной оболочке обусловлено в основном осевой деформацией. Так, осадка включения составляет около 45%, а радиальная деформация - на порядок меньше и не превышает 5% (см. рис. 3,6). Следует отметить неоднозначное влияние содержания никеля на реологические свойства материала. С одной стороны, с увеличением содержания инертного Ж температуры горения и компактирования уменьшаются, а сопротивление уплотнению повышается. С другой стороны, с увеличением содержания Ж увеличивается объем жидкой фазы и уменьшается эффективная вязкость твердой фазы р5. Последнее приводит к снижению вязкости материала и повышению его уплотняемости. В результате оказывается, что компактируемость исследованных материалов не зависит от содержания Ж и полное уплотнение происходит при практически одном и том же давлении рк да 96 МПа.
На кривых изменения размеров и плотности можно выделить четыре характерных участка. На начальном участке уплотнение и формоизменение в зависимости от давления рк происходят по линейному закону и с высокой скоростью. Затем отмечается резкое уменьшение темпа уплотнения и формоизменения с появлением горизонтального плато, когда, несмотря на увеличения давления прессования, плотность и размеры включения остаются неизменными. За этой своеобразной площадкой текучести следует участок нелинейного уплотнения, плавно завершающийся достижением беспористого несжимаемого состояния. Заключительный четвертый участок связан с деформированием несжимаемого включения. Сложный многостадийный характер уплотнения установлен и в работе [23] при прессовании горячих продуктов синтеза системы 2г - С, также содержащих жидкую фазу.
Немонотонный характер процессов уплотнения и формоизменения следует связать с особенностями начального структурного состояния продуктов СВС. В отличие от традиционных процессов горячего прессования порошков, когда в исходном состоянии материал имеет насыпную плотность, за счет внутреннего давления примесных газов твердая фаза продуктов СВС находится во взвешенном состоянии с плотностью, меньшей, чем насыпная, и не оказывает сопротивления внешнему механическому воздействию. Закономерности уплотнения и формоизменения на этой стадии определяются исключительно реологическими и триботехническими свойствами материала оболочки. После уплотнения до насыпной плотности твердой фазы рда =0,53 и формирования пространственно устойчивой структуры наступает вторая стадия, когда уплотнение сопровождается вязким течением частиц твердой фазы. Вторая стадия начинается с площадки текучести на кривых уплотнения. При этом появление площадки текучести обусловлено степенным законом ползучести тугоплавких соединений, в частности карбида 77С. В случае нелинейного закона ползучести при п > 1 в начальный момент уплотнения
рК, МПа
рк,ШЫ
И^,ъ/ал
Р!
О 20 40 60 80 100 120
рк, МПа
В
Рк, МПа
Рис. 3. Зависимость размеров и плотности включения от давления компактирования для сплавов 77С - 20%Л7 (а, б), ТгС - 28%М (в) и 77С - 35%№ (г): 1 - высота включения А/; 2 - плотность р,;
о - эксперимент
скорость уплотнения равна или близка к нулю [24]. Появлению площадки текучести также способствует постепенное, а не скачкообразное, увеличение объема активной фазы со по зависимости (37) и мощности диссипации продуктов СВС. По мере развития нелинейно-вязкого течения происходит уплотнение материала до беспористого несжимаемого состояния, при котором согласно (12) величина р5 становится равной р$ =1- В беспористом состоянии с увеличением давления рк происходит незначительное изменение размеров включения. Так как стадия текучести и следующая за ней стадия нелинейного уплотнения являются отражением одного-физического процесса (нелинейно-вязкое течение пористого тела), то при моделировании следует рассматривать не четыре, а три стадии уплотнения и формоизменения, соответствующие определенному макроструктурному состоянию материала. Количественной мерой макроструктуры пористого твердожидкого тела служит эффективная плотность твердой фазы р5. В зависимости от величины р5 продукты СВС при деформировании могут находиться в следующих состояниях:
1) аэродисперсное взвешенное состояние, когда плотность твердой фазы меньше насыпной р5 < р$н; в этом случае диссипацией энергией в продуктах СВС можно пренебречь, и характер уплотнения и формоизменения включения полностью определяется реологическими свойствами материала оболочки;
2) пространственно устойчивая пористая структура, для которой р$н ^ Рз < 1 - Ш ; уплотнение материала с такой структурой происходит за счет вязкой деформации частиц твердой фазы и сопровождается диссипацией внешней энергии; локальная неоднородность деформации частиц дисперсной среды учитывается объемной долей пластически деформируемой фазы ш ;
3) беспористое несжимаемое структурное состояние при ря =1 - ть ; в этом состоянии в силу неизостатичности прессования в сыпучей оболочке может происходить изменение размеров включения.
В зависимости от состава и свойств продуктов СВС возможно отсутствие какого-либо вида структурного состояния и механизма деформирования. При незначительном газовыделении и развитом спекании может происходить самоуплотнение твердой фазы до насыпной плотности и формирование уже на стадии горения пространственно устойчивой структуры. При большом количестве жидкой фазы полное уплотнение может произойти при аэродисперсном состоянии. Для этого достаточно, чтобы объемная доля жидкой фазы составляла не менее Ш1 = 1 - р8н • Наконец, при малом содержании жидкой фазы и высокой вязкости твердой фазы давления прессования может быть недостаточно для достижения беспористого состояния.
Важное значение для технологического проектирования и оптимизации процесса СВС-прессования имеет информация о закономерностях уплотнения оболочки и силовых параметрах. На рис. 4 представлены кривые уплотнения материала оболочки в центральном блоке 2 и
б,кН
2 р^МПа
Рис. 4. Уплотнение зон оболочки:
1 - кольцевая зона;
2 - центральная зона
Рис. 5. Кинетика силовых параметров процесса прессования:
1-й2-е2;3-05;4-&*
кольцевом блоке 3-4 при прессовании сплава 77С 20%№. Аналогичная зависимость имеет место и при прессовании двух других сплавов. Видно, что на всех стадиях уплотнения материал оболочки в кольцевом блоке уплотняется в большей степени, чем в центральном. Соответственно для деформирования кольцевой зоны оболочки в течение всего цикла прессования затрачивается большая часть мощности и усилия пресса. Это подтверждают и результаты расчета по зависимости (61) усилий £>? и вз (рис. 5). Пока продукты СВС находятся в аэродисперсном состоянии осевое давление со стороны блока 2 на пуансон и усилие £?2 отсутствуют. Осевая реакция на пуансон со стороны включения и блока 2 появляется при уплотнении продуктов СВС до насыпной плотности. Начиная с этого момента времени происходит нарастание усилия Q2 сначала с высокой, а затем меньшей скоростью. Усилие ()з монотонно увеличивается с самого начала прессования, превосходя по величине усилие £Э2 в течение всего цикла прессования. Свой вклад в баланс усилий вносит и сила внешнего трения оболочки <2тр : согласно расчетам доля усилия ()тр в суммарном усилии прессования О, составляет 18-20%. Информация о структуре силовых характеристик показывает пути повышения эффективности прессового оборудования: увеличение в общем балансе сил доли усилия £?2. Эта задача может быть решена
двумя способами. Во-первых, кольцевая зона вокруг продуктов СВС заполняется материалом, имеющим меньшую прочность, чем песок [25]. Во-вторых, за счет увеличения высоты кольцевой зоны оболочки по сравнению с центральной можно уменьшить скорость осевой деформации и уплотнения кольцевой зоны [26].
Характеристикой напряженного состояния включения при квазиизостатическом прессовании служит коэффициент неизостатичности £,, равный отношению радиального напряжения к осевому. В моделях с бесконечной оболочкой [3-5] коэффициент £, считается постоянной величиной. В реальном процессе с конечными размерами происходит уплотнение оболочки и изме-
£
О 20 40 60 80 100 120 Рк, МПа
Рис. 6. Зависимость коэффициента неизостатичности от давления прессования
Рис. 7. Изменение радиального давления на матрицу при прессовании:
1 - в блоке 3,2 - в блоке 4
нение показателя неизостатичности. На рис. 6 показано изменение параметра £, в зависимости от давления рк при прессовании сплава 77С - 20% Ж Количественная оценка £ возможна с момента уплотнения продуктов СВС до насыпной плотности твердой фазы. Деформирование сжимаемого включения начинается с напряженного состояния, близкого к одноосному сжатию да 0). По мере уплотнения включения и оболочки коэффициент увеличивается, а неравноос-ность сжатия уменьшается. Максимальное значение показателя неизостатичности для сжимаемого включения составляет £, = 0,65. Деформирование в условиях неравноосного сжатия с малыми значениями £, обусловливает существенно меньшую радиальную деформацию включения по сравнению с осевой. При уплотнении включения до несжимаемого состояния коэффициент £ резко возрастает до значения \ - 0,825 и напряженное состояние становится близко к гидростатическому, при котором ^ = 1. В условиях почти гидростатического сжатия формоизменения несжимаемого включения не происходит (см. рис. 3).
Практическую ценность имеет информация о внутреннем давлении на стенки матрицы, необходимая для расчета на прочность инструмента. Матрица испытывает внутреннее давление со стороны блоков 3 и 4, которое равно соответствующим радиальным напряжениям. На рис. 7 приведены расчетные данные об изменении радиальных давлений рг в блоках 3 и 4 в процессе прессования. Давление со стороны блока 3 больше, чем со стороны блока 4, и его величина должна использоваться для расчета матрицы на прочность.
Таким образом, процесс квазиизостатического прессования как объект моделирования представляет собой сложную систему, и для адекватного описания процесса необходимо учитывать температуру, фазовый состав, агрегатное состояние продуктов синтеза и реологические свойства вещества твердой фазы. Особенность объекта исследования состоит в том, что в продуктах синтеза происходит самоформирование исходного структурного состояния, и это отличает СВС-прессование от прессования традиционных порошковых материалов. В целом тепловая и механическая модели позволили получить количественные оценки закономерностей уплотнения и формоизменения продуктов СВС и силовых параметров процесса. Структурное и агрегатное состояния продуктов СВС отражены в новых реологических моделях пористых вязких тел с жидкой фазой. Для описания процесса деформирования песчаной оболочки также используется новое условие пластичности, учитывающее локальную неоднородность деформации частиц дисперсной среды. Напомним, что моделировался температурный режим и процесс деформирования локального элемента изделия - его центральной зоны. Аналитическое описа-
ние процесса деформирования всего объема изделия, имеющего неоднородное температурное поле и механические свойства, может быть выполнено методом конечных элементов с использованием физических и реологических моделей настоящего исследования.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Мержанов А.Г Самораспространяющийся высокотемпературный синтез: двадцать лет поисков и находок. Черноголовка: ИСМ РАН, 1989. 92 с
2. Щербаков В.А., Грядунов А.Н., Штейнберг А.С. Макрокинетика процесса СВС-компактирования.// ИФЖ. 1992. Т. 63. № 5. С. 583-592.
3. Ковалъченко М.С. Теоретические основы горячей обработки пористых материалов. Киев: Наукова думка, 1980. 240 с.
4. Горохов В.М., Звонарев Е.В., Ковалъченко М.С. Кинетика уплотнения пористых материалов при горячем прессовании в квазиизостатических условиях // Порошковая металлургия. 1978. № 10. С. 23-27.
5. Горохов В.М., Ковалъченко М.С., Роман О.В. Уплотнение и формоизменение пористых материалов при горячем прессовании в условиях неравномерного трехосного сжатия // Порошковая металлургия. 1983. № 1. С. 8-13.
6. Рогачев А. С. и др. Безгазовое горение в системе титан-углерод-никель // ФГВ. 1988. Т. 24. № 6. С. 86-93.
7. Диаграммы состояния металлических систем. М.: Металлургия, 1964. 232 с.
8. Дулънев ГН., Заричняк Ю.П. Теплопроводность смесей и композиционных материалов. Д.: Энергия, 1974. 328 с.
9. Андриевский А.Р., Спивак ИИ Прочность тугоплавких соединений и материалов на их основе. Челябинск: Металлургия, 1989. 368 с.
10. Пехович А.И, Жидких В.М. Расчеты теплового режима твердых тел. Л.: Энергия, 1976. 352 с.
11. Кубашевский О., Олкок С.Б. Металлургическая термохимия. М.: Металлургия, 1982. 392 с.
12. Лыков А.В. Теория теплопроводности. М.: Высш. школа, 1967. 600 с.
13. Федотов А.Ф. Реологические свойства пористого вязкого тела с жидкой фазой.// Изв. вузов. Сер.- Машиностроение. 1991. № 10-12. С. 8-14.
14. Андриевский РА. Порошковое материаловедение. М.: Металлургия, 1991. 205 с.
15. Феноменологические теории прессования порошков. Киев: Наукова думка, 1985.140 с.
16. Друянов Б.А. Прикладная теория пластичности пористых тел. М.: Машиностроение, 1989. 168 с.
17. Скороход В.В., Штерн М.Б., Мартынова ИФ. Теория нелинейно-вязкого и пластического поведения пористых материалов // Порошковая металлургия. 1987. № 8. С. 23-30.
18. Скороход В.В., Олевский Е.А., Штерн М.Б. Континуальная теория спекания // Порошковая металлургия. 1993. № 1. С. 22-27.
19. Большин Ю.М., Кипарисов С.С. Основы порошковой металлургии. М.: Металлургия, 1978. 184 с.
20. Стоев П.И, Папиров ИИ, Тихинский Г.Ф., Васильев А.А. Диаграммы изостатического прессования бериллия // ФММ. Т.78. №1. С. 9-19.
21. Федотов А.Ф. Характеристики пластичности и внешнего трения сыпучих материалов теплоизолирующей оболочки для прессования продуктов СВС // Огнеупоры и техническая керамика. 1997. № 7. С. 14-17.
22. Друккер Д., Прагер В. Механика грунтов и пластический анализ или предельное проектирование // Механика. 1975. Вып. 2. С. 166-177.
23. Епишин К.Л., Питюлин А.Н., Мержанов А.Г Уплотнение материалов, образующихся при СВС // Порошковая металлургия. 1992. № 6. С. 14-19
24. Ковалъченко М.С. Особенности кинетики уплотнения порошков ковалентных веществ при горячем прессовании // Порошковая металлургия. 1990. № 5. С. 19-24.
25. Амосов А.П., Федотов А.Ф. Способ получения изделий из порошковых материалов // Патент № 2077411 (Россия). Б. И. 1997. №11.
26. Амосов А.П., Федотов А.Ф. Способ получения изделий из экзотермических порошковых материалов и устройство для его осуществления // Патент № 2060866 (Россия). Б. И. 1996. №15.
