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

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

CC BY
111
53
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
КРИСТАЛЛИЗАЦИЯ СПЛАВА С ДИФФУЗИОННЫМ ПЕРЕОХЛАЖДЕНИЕМ / ПАРАМЕТР ДЕНДРИТНОЙ СТРУКТУРЫ / МОДЕЛЬ ЯЧЕИСТОЙ ДВУХФАЗНОЙ ЗОНЫ / CRYSTALLIZATION OF AN ALLOY WITH DIFFUSION SUPERCOOLING / DENDRITIC STRUCTURE PARAMETER / MODEL OF A CELLULAR TWO-PHASE ZONE

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

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

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

Link of the Dendritic Structure of the Alloy Solidifying in a Two-phase Zone with Diffusion Supercooling in the Model of Cellular Two-phase Zone

The papers considers the links of the dendritic structure of metal solidifying in a two-phase zone, with a diffusion overcooling in it, through the use cellular models of two-phase zone. Considered theoretically the experimental results obtained by studying the connection of the diffusion overcooling in a two-zone crystallized alloy with a diffusion overcooling in it. The analysis is based on the use of cellular models of two-phase zone and the numerical solution of the equation. A satisfactory agreement of the change in the diffusion of overcooling in a two-phase zone, depending on its temperature (measured from the liquidus temperature of the alloy). Confirmed the theoretically derived relation between the parameter setting of the dendritic structure with a maximum diffusion overcooling of the alloy at a given cooling rate.

Текст научной работы на тему «Связь дендритной структуры сплава, затвердевающего в условиях двухфазной зоны с диффузионным переохлаждением в ней в модели ячеистой двухфазной зоны»

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

Выводы

Анализируя полученные результаты можно сделать следующие выводы:

• разработанный способ азотирования в плазме тлеющего разряда с полым катодом позволяет создавать регулярную неоднородную структуру на поверхности стали 13Х11Н2В2МФ-Ш и титанового сплава ВТ6;

• установлено, что, помимо общего упрочнения поверхности, происходит более высокое зонное упрочнение вследствие неоднородности плазмы;

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

Литература

1. Малинов Л.С., Соколов Б.К. и др. Получение высоких прочностных и пластических свойств двухфазной стали дифференцированной обработкой // МиТОМ. 1980, № 3, с. 3235.

2. Малинов Л. С. Применение дифференцированных обработок для создания в сталях и чу-гунах регулярных макроскопических градиентов структурно-фазового состояния - перспективное направление в повышении свойств // Металл и литье Украины. 2004, № 11, с. 14-18

3. Григорьянц А.Г., Сафонов А.Н. Основы лазерного термоупрочнения сталей и сплавов. М.: Высш. школа, 1988. 158с.

4. Спиридонов Н.В., Кобяков О.С., Куприянов И.Л. Плазменные и лазерные методы упрочнения деталей машин. Минск: Высш. школа, 1988. 155с.

5. Электронно-лучевая технология: Пер. с нем./ З. Шиллер, У. Гайзиг, З. Панцер: Энергия, 1980. 528с.

6. Будилов, В.В., Агзамов Р.Д., Рамазанов К.Н. Технология ионного азотирования в тлеющем разряде с полым катодом // МиТОМ. 2007. № 7. с. 25 - 29.

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

ячеистой двухфазной зоны

к.ф.-м.н. доц. Сулимцев И.И.

МГТУ «МАМИ» (495) 223-05-23, доб. 1338

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

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

Ранее [1] в связи с проблемой управления кристаллизацией поковки был рассмотрен метод изучения структурообразования в однородном (физически малом, элементарном) объеме двухфазной зоны бинарных сплавов. Полученные экспериментальные результаты пока-

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

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

Рисунок 1 - Условия затвердевания цилиндрической ячейки в элементарном объеме двухфазной зоны: X - среднее расстояние между характерными элементами структуры сплава, диаметр ячейки; - закон роста кристалла в ячейке; г - координата;

Сф - состав жидкой фазы на поверхности кристалла; С2 - состав жидкой фазы на границе ячейки; С(т,1) - состав жидкой фазы внутри ячейки; \Со - пересыщение в

ячейке; Со - исходный состав сплава

В отличие от работы [2], в которой закон роста ячейки двухфазной зоны задавался в виде у(() = 41, будем считать закон роста неизвестным. В таком случае полная система уравнений, описывающих кристаллизацию ячейки должна содержать условие, информацию о котором можно получить независимым от решаемой задачи способом, например, на основе тепловых условий затвердевания (термических кривых). Это дополнительное условие можно получить следующим образом. Запишем состав примеси на фронте кристаллизации как функцию времени: С), 1 ) = Сф (() (1) .

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

2 - 102 см.

г

ции можно следующим образом выразить через известную функцию Т@) двухфазной зоны:

температуру

Сф (,) = , (2)

а

или в общем случае: Сф(,) = <р(Т(()) , (3)

где функция ф(Т) по существу представляет собой уравнение линии ликвидуса.

Заметим, что здесь температура Т@) двухфазной зоны уже не отсчитывается от ликвидуса сплава, что, впрочем, не имеет существенного значения. Связь вида (2) и является дополнительным условием, замыкающим систему уравнений, описывающих кристаллизацию ячейки в условиях проводимого эксперимента. Эта система уравнений такова: ОЛС = С,(г,,),у(()< г < А/2,, > 0

С (А/2, ,) = 0, (4)

С (г,0) = Сс,

БСГ(у((),,)+(1 - к)(у(),,) = 0,

Сф (() = (ТА - Т(())/

а.

Здесь а - наклон линии ликвидус, ТА - температура плавления чистого компонента. Система ( 4 ) принадлежит к классу обратных задач, является замкнутой и может быть решена на ЭВМ, т.е. могут быть найдены распределение концентрации примеси С(г, ,) перед фронтом, закон движения фронта у(,), сечение £ (относительная доля) жидкой фазы и зависимость от времени величины диффузионного переохлаждения ЛТа = аЛСв =а{Сф - С2),

где С2 - концентрация примеси на границе ячейки (рисунок 1).

Таким образом, как видим, при данном теоретическом рассмотрении элементарный объем кристаллизующегося в эксперименте сплава представляется в виде некоторого количества эквивалентных и одинаковым образом кристаллизующихся ячеек, связанных с характерными структурными составляющими сплава. Теоретически, следовательно, изучается процесс кристаллизации в одной, являющейся характерной для исследуемого сплава ячейке в том смысле, что наблюдаемое в опытах распределение диффузионного переохлаждения ЛТа (ЛТ) описывается значениями разности а(Сф - С2).

3

Г\

и \ и \ \ \ \ \ \

\ \ \ \ \ N \

Рисунок 2 - Пример сравнения результатов теоретического расчета (пунктир) с данными эксперимента (сплошная линия). Сплав 8п + 45 вес.% Zn, скорость охлаждения Ж = 0,20/с. Верхняя кривая получена при коэффициенте диффузии примеси в жидкой фазе Б = 1,2-10"6см2/с, нижняя - при Б = 2-10"5см2/с.

На рисунке 2 показан результат применения описанной выше теоретической схемы к

расчету распределения диффузионного переохлаждения в двухфазной зоне сплава Бп+45 вес. % 2п при скорости охлаждения 0,2 0/с. Необходимые для проведения расчета экспериментальные данные содержатся в таблице 1. Как видно из рисунка 2, теоретическое рассмотрение правильно отражает физическую сущность происходящих в условиях эксперимента процессов.

Таблица 1

Экспериментальные данные, использованные для расчета на ЭВМ, результаты которых приведены на рисунке 2 для сплава 8п + 45 вес.% Zn

Сплав Sn+45%Zn Cd+30%Zn Pb+40%Sn

Параметр

С0, %

0 С

а, — %

К

д, кал/г

2Б,

ж,

2

см с

кал г • град

45

1,8

0 17 2-10"5

(8-10"5) 0,1

30

1,5

0 16,3 4-10

-6

20

х2/т, см

с

5* Ф(5*)

ат;, 0с

теория

ат;, 0с

эксперимент

1,11-10

5,05-10"'

4,07-10"' 0,10

8,1

9,2

-3

(10-5)

0,07 8,52-10

1,64-10-'

5-10-2 0,125

7

8,6

60

2,2

0,4 9,2 10-6

(1,4-10-5)

0,05

-5

1,40-10

2,93-103,66-10-: 3,66-10-:

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

4,4 8

-5

Поскольку в [2] выражение, устанавливающее связь между АТ*; и X (максимальное диффузионное переохлаждение в двухфазной зоне и среднее расстояние между элементами дендритной структуры сплава соответственно) содержит еще и градиент О температуры в двухфазной зоне, а в проводимом эксперименте выполнялось условие О ~ 0, то непосредственное использование полученного там соотношения невозможно. Однако, используя идеи работы [2], можно получить аналогичное соотношение, но не содержащее градиента температуры.

Будем рассматривать затвердевание элементарного объема как процесс кристаллизации некоторого количества цилиндрических ячеек (связанных с дендритами), среднее расстояние между которыми есть X (рисунок 1). Ячейки будем считать одинаковыми в отношении их размеров и происходящих в них кристаллизационных процессов. Пусть I - высота ячейки. Как и в [2], полагаем, что затвердевание происходит путем радиального роста (утолщения) оси ячейки (заштриховано на рисунке 1). При этом заметим, что, как показали соответствующие расчеты, конкретная форма ячейки несущественно влияет на окончательные результаты вычислений. Обозначим через у(1;) радиус растущего в ячейке кристалла в момент времени 1 от начала кристаллизации. Найдем относительное увеличение объема твердой фазы за единицу времени для всех ячеек кристаллизующегося в эксперименте сплава.

с

Эта величина будет равна Ss■vn / V, где Бб - боковая поверхность цилиндра радиуса у(1:) и высоты I; ип - нормальная скорость роста кристалла, равная ёу(1)/ё1:; V - объем ячейки, равный п(Л/ 2)21.

_ 2луШу /йг л йу2 . .2

Следовательно, искомая величина есть —-т^ = 4-/ Л .

п(Л /2)21 йг

Если М - масса, а q - теплота кристаллизации сплава (принимается, что q - постоянная величина на интервале затвердевания), то за единицу времени выделяется количество тепла,

йу2

равное / Л2.

С другой стороны, если W - скорость охлаждения сплава, то количество тепла, отбираемого за единицу времени, равно WMх (х - средняя теплоемкость сплава). Поэтому, приравнивая два последних выражения, получаем:

Л2 = 4^ (5)

йг Жх

Закон роста у(1), фигурирующий в выражении (5), возьмем, согласно [2], в виде у (г )= л/4Б& ; здесь Б - коэффициент диффузии примеси в жидкости; б - постоянная величина. Если положить S = ЛЖ /16Бц, что соответствует условию (5), то можно воспользоваться результатами расчета, проведенного в [2]. Согласно этой работе максимальное диффузионное переохлаждение АТ* в ячейке связано с параметром б соотношением АТ* / aS = р(s), где ф - определенная функция б и к (к - коэффициент распределения), полученная в табличной форме при решении задачи на ЭВМ. Обозначив через ф-1 функцию, обратную к ф, получим для X выражение:

Л = д/16 Бд <р - (а Т* / аС 0 )/ Ж , (6)

в котором а - наклон линии ликвидуса, С0 - состав сплава.

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

Поскольку затвердевание в ячейке заканчивается в момент времени т, то Л/2 = 44Шт, откуда Л = д/16Бтр-(АТ* /аС0) . (7)

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

С целью экспериментальной проверки полученных выше результатов на поперечных шлифах исследованных сплавов для различных скоростей охлаждения элементарного объема были найдены средние расстояния между характерными элементами структуры. Исследование строения дендритных каркасов, остающихся на опорной сетке после экстрагирования междендритной жидкости, показывает, что для сплавов Бп+45 вес.% 2п, Бп+20 вес.% 2п, Сё+30 вес.% 2п первичные оси дендритов вместе с тесно прилегающими друг к другу осями высших порядков образуют плотные листочки-пластины, которые можно считать характерными элементами дендритной структуры сплава. На шлифе наблюдаются следы пересечений этих листочков-пластин плоскостью шлифа. Среднее расстояние между следами принималось за среднее расстояние между ячейками в двухфазной зоне затвердевающего элементарного объема исследуемого сплава. В случае сплава РЬ+40 вес.% Бп дендритная структура

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

Время т затвердевания ячеек двухфазной зоны в объеме кристаллизующегося сплава определялось по термической кривой Т^) непрерывно затвердевающего сплава.

Численные значения величины ф(5) можно найти либо из таблицы 1, либо графически (рисунок 3).

Рисунок 3 - Зависимость функции ф от параметра 5 (43): 1. к = 0; 2. к = 0,25; 3. к = 0,5;

4. к = 0,75; к - равновесный коэффициент распределения

В проводимых в настоящей работе расчетах значения величины 5 получались ~10-2. Ввиду малости 5 и согласно рисунку 3, а также таблице в работе [2], при к = 0 использовалось соотношение ф(5)=2,55, а при к = 0,4 ф(5) = 5. Результаты проведенных по описанному выше плану расчетов величины АТ*, экспериментальные значения АТ* для изученных сплавов, а также необходимые вспомогательные данные приведены в таблице 1. В ней помещены максимальные для каждого сплава значения величины 5, обозначаемые через 5*. Именно значениям 5* соответствуют помещенные в таблице результаты вычислений функции ф(5*). Помещенные в скобках коэффициенты диффузии приводят к тем же значениям АТ*, что и выписанные без скобок, если значение X принимать равным не величине АТ* ), соответствующей значению X при максимальном диффузионном переохлаждении в двухфазной зоне, а равным начальной (наибольшей) величине параметра X для соответствующего сплава при данной скорости охлаждения.

Физический смысл упомянутых выше условий нахождения АТ* состоит в следующем. В ячеистой теории двухфазной зоны величина X есть неизменное в процессе затвердевания расстояние между центрами кристаллизующихся ячеек. Ранее, в соответствии с положением указанной теории о максимальной величине диффузионного переохлаждения в двухфазной зоне, в качестве X принималось значение, равное X(АТ*). Однако результаты описанных опытов указывают на то, что практически конечное число дендритов образуется в самом начале (при ДТ~0,1°С) затвердевания элементарного объема двухфазной зоны. Другими словами, среднее расстояние между центрами кристаллов, которое отвечает величине X в ячеистой теории, практически не меняется в процессе эволюции элементарного объема двухфазной зоны.

Изложенные результаты позволяют сказать, что вытекающая из ячеистой модели двухфазной зоны связь между дисперсностью дендритной структуры X и скоростью охлаждения Ж имеет место при затвердевании физически малого (элементарного) объема двухфазной зоны.

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

Раздел 2. Технология машиностроения и материалы. Литература

1. Сулимцев И.И. Метод изучения структурообразования металла с целью получения информации для управления процессом затвердевания изделий, кристаллизующихся в условиях двухфазной зоны. Известия МГТУ «МАМИ», 2009, №2(8).

2. Тяжельникова И.Л., Борисов В.И., Борисов В.Т. Кристаллизация цилиндрической ячейки двузфазной зоны. Известия АН СССР, Металлы, 1970, № 5.

Метод измерения малых сегментов цилиндрических поверхностей

к.т.н. проф. Суслин В.П., к.т.н. доц. Джунковский А.В.

МГТУ «МАМИ» (495) 223-05-23 доб. 1392, labsapr@narod.ru

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

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

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

В [3, 4] рассмотрены задачи измерения малых сегментов окружностей и сфер, показана их плохая обусловленность и на практических примерах получены результаты, демонстрирующие большие ошибки измерений в таких ситуациях. Предложенная в этих работах идея использования метода регуляризации академика А.Н. Тихонова [1, 2] позволяет улучшить обусловленность задач измерения на малых областях и существенно уменьшить ошибки измерений. В данной статье приведены результаты исследований по измерениям малых сегментов цилиндрических поверхностей.

Параметрами, требующими определения при измерении цилиндра, являются: радиус R, вектор оси V=(Vx, Vy, Vz) и центральная точка C=(x0, y0, z0). Для решения задачи возможно применение разных методик. Мы используем метод наименьших квадратов, основанный на непосредственной минимизации суммы квадратов отклонений измеренных точек от цилиндрической поверхности:

n „

F(R, V,a,b,c) = X а2 ^ min, (1)

i=1 i

где: R, V, a,b,c - радиус, вектор оси и координаты центральной точки цилиндра;

А = \Pi -Pi\, i=1,2,...,n - отклонения измеренных точек от цилиндра;

Pi = (Xi, Yi, Z) - измеренная точка;

' II!

P^ = (Xi, Yi, Z.) - точка на цилиндре, соответствующая измеренной. Реализуя метод регуляризации, дополняем функцию (1) отклонениями параметров от

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