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

004603927

На правах рукописи

ПОЛИЩУК Игорь Николаевич

ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ ФАЗОВЫХ ПРЕВРАЩЕНИЙ ПЕРЕОХЛАЖДЕННОГО АУСТЕНИТА В СТАЛИ

05.13.18 - Математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

1 О ИЮН 2010

Тюмень-2010

004603927

Работа выполнена на кафедре информационных систем Института математики и компьютерных наук ГОУ В ПО Тюменский государственный университет

Научный руководитель: доктор технических наук, профессор

Ивашко Александр Григорьевич

Официальные оппоненты: доктор технических наук, профессор

Инютин Сергей Арнольдович доктор химических наук, профессор Андреев Олег Валерьевич

Ведущая организация: ГОУ ВПО Тюменский государственный

нефтегазовый университет

Защита диссертации состоится «25» июня 2010 г. в 1400 часов на заседании диссертационного совета Д 212.274.14 при Тюменском государственном университете по адресу 625003, г. Тюмень, ул. Перекопская, 15А, ауд. 410.

С диссертацией можно ознакомиться в библиотеке Тюменского государственного университета.

Автореферат разослан «24» мая 2010 г. Ученый секретарь

Н.Н. Бутакова

диссертационного совета " ~ ~

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

В отечественной школе фундаментальный обзор математического моделирования процесса фазовых превращений был проведен в работе Б.Я. Любова. Наиболее полная зарубежная работа, освящающая основные направления в моделировании фазовых превращений является работа Дж. Кристиана. В обоих трудах указано, что аналитические модели не учитывают явно процесс нестационарного зарождения, что противоречит физике процесса. Также отмечено, что аналитические модели не позволяют рассматривать процесс образования отдельных зерен и их роста.

Помимо аналитических моделей существует подход, заключающийся в имитации процесса зарождения и роста на основе модели И.Л. Миркина, в

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

В настоящее время проблемами моделирования фазовых превращений в нашей стране занимаются несколько школ, среди которых можно выделить школы под руководством академика В.Н. Анциферова (Пермский государственный технический университет) и профессора Ю.Г. Гуревича (Курганский машиностроительный институт), которые занимаются моделированием фазовых превращений в порошковых сталях. Группа авторов под управлением Д.А. Мирзаева (Южно-Уральский государственный университет) занимается моделированием структурных превращений в сталях. В качестве примера моделирования процесса зарождения и роста в специальных сплавах можно привести работы В.Я. Шура (Уральский государственный университет им. A.M. Горького), посвященные изучению эволюции сегнетоэлектрической доменной структуры в электрическом поле.

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

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

Для достижения поставленной цели определены следующие задачи исследования:

1) создание векторного представления структуры стали в трехмерном пространстве, которое даст возможность описать каждую структурную составляющую в любой момент времени и в любой точке пространства термообрабатываемой детали;

2) создание стохастической модели зарождения и роста структурных составляющих в результате распада аустенита при термической обработке стали;

3) разработка алгоритма определения кинетических параметров векторной имитационной модели распада аустенита по данным кинетических диаграмм полученных в изотермических условиях;

4) разработка комплекса программ, реализующих разработанные модели и алгоритмы, доступ к которым осуществляется посредством \veb-технологий;

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

Методы исследований. При решении поставленных задач используются аппарат теории массового обслуживания, методы оптимизации, математическая статистика, аналитическая геометрия, \УеЬ-технологии, СЯЮ-технологии.

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

На защиту выносятся:

1) векторная модель структуры стали, полученной в результате распада аустенита при термической обработке стали;

2) векторная имитационная модель распада аустенита;

3) модифицированный метод Хука-Дживса, нахождения кинетических параметров с учетом статистического критерия;

4) архитектура программного комплекса имитационного моделирования распада аустенита.

Научная новизна работы заключается в следующем:

1) создана векторная модель структуры стали, позволяющая описать каждую структурную составляющую в любой момент времени и в любой точке пространства;

2) разработана векторная имитационная модель распада аустенита, которая дает возможность реализовать математический аппарат теории массового обслуживания, что позволяет проводить имитационные эксперименты на основе событийного представления кинетики фазового превращения;

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

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

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

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

1) Научные исследования по тематическому плану научно-исследовательских работ №73-01 от 27.11.2006 «Разработка теоретических основ влияния пористости на кинетические параметры распада переохлажденного аустенита». Регистрационный номер НИР 01.20703455.

2) Аналитическая ведомственная целевая программа Министерства образования и науки Российской федерации «Развитие научного потенциала высшей школы (2009-2010)» - «Прогнозирование кинетики распада аустенита в порошковых сталях при непрерывном охлаждении». Регистрационный номер НИР 2.1.2/6498;

Апробация работы. Основные результаты диссертационной работы были представлены на четырех всероссийских (Пенза, Новосибирск, Томск) и двух международных конференциях: «Математические методы в технике и технологиях» (Ярославль) и «Информационно-вычислительные технологии и их приложения» (Пенза).

Публикации. Основное содержание работы отражено в 13 публикациях, из которых 2 свидетельства о государственной регистрации программы для ЭВМ и 2 статьи опубликованы в журнале из списка ВАК.

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения и списка литературы. Объем диссертации составляет 130 страниц, содержит 24 рисунка, 2 таблицы, 2 приложения. В библиографии представлено 117 наименований работ российских и зарубежных авторов.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

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

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

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

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

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

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

Рассмотренное математическое моделирование процесса в свою очередь разделяется на построение аналитических и имитационных моделей. Большинство аналитических моделей строятся на основе модели А.Н. Колмогорова и уравнения М. Авраами, которые на основе двух кинетических параметров, скорости роста и вероятности зарождения, позволяют получать долю вновь образуемой фазы во времени. Также существуют модели, которые описывают процесс превращения дифференциальным уравнением. Но получить кинетические параметры на основе экспериментальных данных в аналитических моделях не представляется

8

I

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

Хотя модели А.Н. Колмогорова и М. Авраами описывают стохастический процесс, на практике используются модели, в которых вероятность зарождения центров структурных составляющих заменяется ее усредненной характеристикой - скоростью зарождения. Таким образом, стохастическая природа процесса игнорируется. Также аналитические модели не учитывают наличие инкубационного периода, доказанного экспериментально. Согласно этим моделям превращение начинается в нулевой момент времени.

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

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

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

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

Предлагается векторная модель структуры стали, в которой образующиеся в результате распада аустенита структурные составляющие описывается в виде объектов, для которых хранится уравнение описывающее поверхность. Так, например, в случае если структурная составляющая изоморфна, формой объекта является шар, который характеризуется координатами центра (х„ _>>„ г,) и радиусом г,.

Для определения доли объема занятого вновь образующимися объектами применяется метод Монте-Карло. В соответствии с этим методом в область

исследуемого пространства равномерно выбрасываются случайные точки с координатами (х, у, г). Если образуются объекты одного вида, то для каждой точки проверяется выполнение неравенства по всем объектам:

(х-х^-иЫ*-*,)2^. О)

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

В том случае если образуются объекты различного вида (например, в случае образования феррито-перлитной структуры из аустенита) подсчитываются точки, выбрасываемые в методе Монте-Карло, принадлежащих каждой структурной составляющей по формуле (1). В том случае, если точка принадлежит двум объектам разного вида, идентификация структурной составляющей осуществляется по следующему алгоритму. Вычисляется параметр сг,, для /-той структурной составляющей:

С; = ,У*,2*)-б(х1, У;, 2,- ), (2)

где:

(х», >'., 2.) - координаты выброшенной точки;

(х„ у,, г,) {Хр у), г,-) - координаты центров объектов для которых выполняется неравенство (1),

3{х, у, г) = [х] - х,-) • (х - х0) + (у / -У/)-{у-Уо)+ {г} - г,-) • (г - г0), (3)

П+/-Г, I \

= * + - ' '- - (4)

Уо = У, + I/ х, / V, • Ь>1 -УХ (5)

гь+1-т,

Если а, > 0, то точка принадлежит к структурной составляющей соответствующей /-му объекту.

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

г/ = г4 + / г,-, (7)

где: г/,. - начальный размер; т, - время, прошедшее с момента возникновения ¿-го объекта; / - линейная скорость роста.

Зарождение нового объекта возможно тогда, когда в систему поступает

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

1. стационарность - вероятность зарождения в течение заданного промежутка времени определяется только длиной этого промежутка и не зависит от расположения данного промежутка на временном отрезке процесса превращения;

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

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

Из чего следует, что рассматриваемый поток является простейшим, что позволяет применять модель массового обслуживания при описании процесса зарождения структурных составляющих. Существует конечное количество каналов - возможных центров зарождения. Это могут быть области с флуктуациями химического состава, дефекты кристаллической решетки, и т. п. В систему поступает поток заявок на возникновение зародышей. Центр зарождения становится зародышем новой структурной составляющей в момент, когда на свободный канал приходит заявка. При этом число занятых каналов увеличивается на единицу. Канал, обслуживающий поступившую заявку, остается занятым до самого конца фазового превращения, то есть на все время рассматриваемого процесса. Также возможный центр зарождения считается занятым до конца эксперимента, если он принадлежит новой фазе. Заявки приходят с равной вероятностью на любой канал, вне зависимости от его занятости. Таким образом, модель представляет собой многоканальную систему с отказами с непрерывным временем.

Так как поток поступления заявок на зарождение является простейшим, моменты времени т1, г2,..., тк,... поступления заявок в систему определяются следующим уравнением:

п

тп =Х (8)

/=1

где: = ~■ 1п(1 -И, ); я,-равномерно распределены на интервале (0,1).

Л

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

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

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

Описанный подход дает возможность применять методы и алгоритмы теории массового обслуживания и проводить имитационный эксперименты на основе событийного представления кинетики фазового превращения.

Ход вычислительного эксперимента для заданной температуры превращения описывается следующим алгоритмом:

ШАГ 1: Задается число итераций эксперимента (К); начальный номер итерации .7=0;

ШАГ 2: Очистка пространства превращения;

ШАГ 3: Генерация событий на зарождение и помещение их в очередь;

ШАГ 4: Генерация событий на уничтожение и помещение их в очередь;

ШАГ 5: Определение !к - ближайшего события;

ШАГ 6: Пересчет состояния системы в момент времени ^ и определение объема новых объектов методом Монте-Карло (К*);

ШАГ 7: Переход на ШАГ 3;

ШАГ 8: Если объем новых структурных составляющих соответствует заданному объему аустенита - превращение заканчивается. Для данной итерации вычислительного эксперимента массив {Ук, 4) отображается в массив (г, где г массив фиксированного процентного содержания превращенного аустенита (например, 0,1; 0,2 и т.д.);

ШАГ 9: Если Л{К-1), то ./= 7+1 и идти на ШАГ 2;

ШАГ 10: Усреднение Л - времени образования фиксированного процентного состава 2 по всем итерациям вычислительного эксперимента У.

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

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

/Ы

где: Ут, и V, - доля распавшегося аустенита в момент времени 1к при температуре Г„ полученного в имитационной модели из экспериментально построенных кинетических диаграмм;

и, - число точек экспериментальной кинетической диаграммы при температуре Тс,

А/Т) - интенсивность потоков зарождения и уничтожения;

1(Т[) - линейная скорость роста;

гь(Т) и гс(Т.) - начальный критический радиус объектов зарождения.

Для решения задачи был выбран метод Хука-Дживса, который относится к методам прямого поиска и использует только значения функции. Значение функции Ут{ при заданных ).-_, /, /'/„ г, в каждый момент времени ^ представляет собой случайную величину. В вычислительном эксперименте на шаге 10 хотя и проводится усреднение этих случайных величин, но это не дает возможность исключить стохастическую составляющую этого параметра. Для того чтобы при поиске решения в методе Хука-Дживса проверить гипотезу о том, что улучшение значения целевой функции обусловлено изменением кинетических параметров, а не связано со случайными факторами, использовался статистический критерий.

Модифицированный метод Хука-Дживса, нахождения кинетических параметров для каждой температуры Т„ с учетом статистического критерия описывается следующим образом. Выбирается начальная базисная точка ад =(/1го,1о,г/юЛсЮ,Гсп) и длина шага А,- для каждой из переменных. Производится серия компьютерных экспериментов, по результатам которой вычисляется значение АУт{а0)

АГт,(а0) = Т,), /(Т,),гь(Т,),Я/СТ,),гс(Т,), /, )-К(/(.))2 (10)

к=]

Значение первой переменной увеличивается на длину шага. Проводится новая серия экспериментов и вычисляется значение А Ут^п+Ь ] ,1о,г,,я Д^.г^). Если АУт;(Х:0+И1,10,гьоЛс10,Гсо) < АУт,(Х,0, гьо, км, гс0), то выдвигается гипотеза Но. М(АУт1(к:0+Ь1,10,гьоЛм,гса)) < М(АУт1(Х,о,1о,ПоЛм,Гсо)), для проверки которой используется /-критерий:

^2(ДЮи,(Л.0, /0,гьо,^0,гс0))+52(АКт,.(Лг0 + А„ /0,гЬ0,^0,гс0))'

где N - количество итераций вычислительного эксперимента. Наблюдаемое значение I' сравнивается со значением //.„ с V степенями свободы при уровне значимости а=0,05,

(л-2 (Д Ут, (/.„, /„, гЬ1), Лм, гс0 + 2 (Д Ут, (Лг0 + /г,, /0, гьо ,ЛМ,гс0 й)2

При /' > //.„ гипотеза принимается, в противном случае - отвергается.

В случае принятия гипотезы Нп в качестве новой базисной точки а/ берется а1=(Х2о+И1,1с,гЬо,к1ю,Гсо)\ в противном случае значение первой

переменной уменьшается на длину шага, проводится еще одна серия экспериментов и вычисляется значение &Vm,(k:n-h\, /о, >'ьо, Кю, rct>)- Если Д Vmj(X:o-hi, Iо, По, Хм, rio)< AVmifi,0, ¡о, гьо, Хм, По), то выдвигается гипотеза Н0 : reo)) < М(ЬУт-,(Х го^о.ГьоЛ'м^со)), которая проверяется по I-критерию как описано выше. В случае принятия гипотезы в качестве новой базисной точки берется <ц=(л-л-ht, /о, гьп, Хм, гсП).

Если ни один из перечисленных шагов не приводит к уменьшению значения функции, то изменению подвергается значение второй переменной и выполняется проверка соответствующих гипотез по вышеописанному алгоритму. Аналогичные шаги выполняются по всем переменным. В результате перебора по каждой из переменной алгоритм получает следующую базисную ТОЧКу Я/.

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

Кинетические параметры структурного превращения, полученные в методе Хука-Дживса, для каждой температуры использовались при построении интерполяционных функций lz(T), 1(Т), rh(T), Xj(T), rL.(T) с помощью кубических сплайнов. Полученные зависимости кинетических параметров от температуры применялись при прогнозировании кинетики превращения аустенита для заданной кинетики охлаждения.

Алгоритм имитации кинетики распада аустенита при непрерывном охлаждении аналогичен алгоритму вычислительного эксперимента в изотермических условиях, который был описан выше. При этом учитывалась зависимость кинетических параметров от температуры, а также зависимость температуры от времени T=f(т).

Размер объекта новой фазы определялся по формуле:

I

П = гь+ \l{f(r))dvt (13)

Гу

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

t>,+>"<

' ¡A{f(t))dr =-InRk+i, (14)

Р'

где: Я функции интенсивности зарождения или уничтожения.

Учитывая, что функции /(z), l(z), являются полиномиальными, вычисление интеграла проводилось методом Ныотона-Котеса.

Третья глава посвящена рассмотрению архитектурного решения и практике создания программного комплекса реализующего разработанную имитационную модель.

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

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

2. Расчет кинетических параметров моделей по данным эксперимента и их сохранение в базе данных.

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

4. Прогнозирование изотермических и термокинетических диаграмм марок сталей, для которых известны кинетические параметры моделей.

Расчет кинетических параметров производится с помощью модифицированного метода Хука-Дживса на основе экспериментальных диаграмм распада аустенита при изотермических условиях. В алгоритме реализована векторная модель распада аустенита с нестационарным зарождением. Для нахождения кинетических параметров для одной температуры необходимо около ста шагов. При каждом шаге производится усреднение полученных значений по сорока итерациям. В итерации при каждом расчете объема новых объектов методом Монте-Карло в область исследуемого пространства выбрасывается десять тысяч точек. Для каждой точки последовательно проверялась принадлежность каждому из вновь образованных объектов. Таким образом, для нахождения кинетических параметров модели необходимо большое количество вычислительных ресурсов

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

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

Модуль Web-интерфейса - это классы и Web-страницы для взаимодействия с пользователем, ввода им данных и отображения результатов экспериментов в табличном виде и в виде графиков. Классы Web-интерфейса и хранимые в нем http-страницы, публикуются через MS Internet Information Server, что дает возможность пользователям получать Web-доступ к возможностям программного комплекса. Модуль сервера базы данных служит для взаимодействия пользователя через Web-интерфейс с хранимой в базе данных информацией. С его помощью сохраняются все введенные экспериментальные кинетические диаграммы, рассчитанные кинетические параметры, а также рассчитанные на их основе термокинетические и изотермические диаграммы.

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

СЕРВЕР

.- -Сервер приложения

Модуль расчетов Сервер БД

Кднеет GÏUD WEB интерфейс

.NET Framework

БД.

MS SQL Server

\

♦Л

MS Internet Information Server (US'!

СЕРВЕР

Сервер приложения

Сервер GRID

Модуль расчетов

J

КЛИЕНТ

Internet Explorer ^jj

Рис.1. Архитектура программного комплекса

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

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

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

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

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

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

0,006 -I

0,004 ] 0,002 {

10 20 30 40 50 60

Количество итераций

Рис. 2. Дисперсия целевой функции (ЛКот;) в зависимости от количества

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

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

В ходе проведения эксперимента на основе экспериментальных диаграмм распада аустенита при изотермических условиях для стали ПК40Н2М был проведен расчет кинетических параметров для одиннадцати температур выдержки. На основе полученных кинетических параметров было проведено компьютерное моделирования распада аустенита при постоянной выдержке для тех температур, что были использованы при построении экспериметальной изотермической диаграммы. Данные полученной изотермической диаграммы были сопоставлены с экспериметальными (рис. 4). Приведенные графики и расчеты показали, что полученные в результате моделирования изотермические диаграммы идентичны диаграммам, полученным в ходе реальных экспериментов.

Количество компьютеров

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

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

-Данные эксперимента ........Результаты расчета согласно модели

Рис.4. Изотермическая диаграмма распада аустенита для стали ПК40Н2М

10 100 1000 —— Исходные данные ........Результатмодслнровання Время, сск.

Рис.5. Термокинетическая диаграмма распада аустенита для стали ПК40Н2М Таким образом, было показано, что разработанная векторная имитационная модель может быть использованы для прогнозирования структуры стали для любой кинетики охлаждения.

По результатам исследований могут быть сделаны выводы 1. Предложено векторное представление структурных составляющих стали в трехмерном пространстве. Такое представление дает возможность

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

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

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

4. Модифицирован алгоритм Хука-Дживса для получения кинетических параметров имитационной модели фазовых превращений на основе экспериментально построенных изотермических и термокинетических диаграмм.

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

6. Разработана архитектура программного комплекса на основе >УеЬ-технологий, которая обеспечивает представление результатов машинного моделирования, а также предоставляет пользователям доступ к сервисам программного комплекса. Пользователи могут проводить моделирование фазовых и структурных превращений через сеть интернет.

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

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

ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Ивашко А.Г., Полищук И.Н. Имитационное моделирование фазового превращения переохлажденного аустенита в стали // Современные математические методы и информационные технологии в образовании. Тезисы докладов. - Тюмень: Тюменский Государственный Университет, 2005.-С. 28-29.

2. Ивашко А.Г., Цыганова М.С., Полищук И.Н. Описание модели массового обслуживания для процесса зарождения зерен при переохлаждении аустенита // Проблемы информатики в образовании, управлении, экономике и технике: сборник статей VI Всероссийской научно- технической конференции. - Пенза, 2006. - С. 126- 129.

3. Ивашко А.Г., Полищук И.Н. Имитационное моделирование фазового превращения переохлажденного аустенита в стали // Математическое и информационное моделирование: Сборник научных трудов. Вып. 8. Тюмень: Издательство «Вектор Бук», 2006. - С. 92- 96.

4. Ивашко А.Г., Цыганова М.С., Полищук И.Н. Описание модели массового обслуживания для процесса зарождения зерен при переохлаждении аустенита // Наука технологии инновации // Материалы всероссийской научной конференции молодых ученых в 7- ми частях. Новосибирск: Изд-во НГТУ, 2006. Часть. 1. - С. 25- 27.

5. Ивашко А.Г., Цыганова М. С., Полищук И. Н. Применение теории массового обслуживания при моделировании фазовых превращений // Математические методы в технике и технологиях - ММТТ- 20. [Текст]: сб. трудов XX Междунар. науч. конф. В 10 т. Т. 6. Секция 12 / под общ. ред. В. С. Балакирева. - Ярославль: Изд- во Яросл. гос. техн. ун- та, 2007. - С. 233236.

6. Ивашко А.Г., Цыганова М.С., Полищук И.Н. Имитационная модель распада аустенита // Информационные технологии и телекоммуникации в образовании, экономике и управлении регионом: Материалы межрегиональной научно- практической конференции. 1-30 ноября 2006 г. Тюмень: Издательство Тюменского государственного университета, 2007. -С. 23-25.

7. Полищук И. Н. Применение теории массового обслуживания при моделировании распада аустенита // Информационно-вычислительные технологии и их приложения: сборник статей VI Международной научно-технической конференции. - Пенза: РИО ПГСХА, 2007. - С. 145- 147.

8. Ивашко А.Г., Цыганова М.С., Полищук И.Н. Имитационная модель двухфазного превращения // Молодежь и современные информационные технологии. Сборник трудов V Всероссийской научно- практической конференции студентов, аспирантов и молодых ученых. Томск, 27 февраля -1 марта 2007 г. - Томск: Изд- во ТПУ, 2007. - С. 111- 112.

9. Ивашко А.Г., Цыганова М.С., Полищук И.Н. Имитационное моделирование фазовых превращений при термической обработке порошковых сталей // Вестник Тюменского государственного университета. №5 2007 г. Тюмень: Издательство Тюменского государственного университета. 2007 г. - С. 7580.

\0.Ивашко А.Г., Цыганова М.С., Полищук И.Н. Компьютерное моделирование превращения фаз в процессе термической обработки // Вестник Тюменского государственного университета. №6. 2008 г. Тюмень: Издательство Тюменского государственного университета. 2008 г. - С. 129- 140.

11 .Ивашко А.Г., Цыганова М.С., Полищук И.Н. Модель нестационарного процесса распада переохлажденного аустенита // Молодежь и современные информационные технологии. Сборник трудов VII Всероссийской научно-практической конференции студентов, аспирантов и молодых ученых "Молодежь и современные информационные технологии". Томск, 25- 27 февраля 2009 г., ч.1. Томск: Изд- во СПБ Графике. - С. 169- 171.

12. Ивашко А.Г., Полищук И.Н. Свидетельство о государственной регистрации программы для ЭВМ №2009616252 «GRID - система прогнозирования распада переохлажденного аустенита» от 11.11.2009.

13.Ивашко А.Г., Полищук И.Н. Свидетельство о государственной регистрации программы для ЭВМ №2009616253 «Подсистема имитационного моделирования распада переохлажденного аустенита при изотермических условиях» от 11.11.2009.

Издательство «Веетор Бук» Лицензия ЛР № 066721 от 06.07.99 г.

Подписано в печать 24.05.2010. Формат 60x84/16. Бумага офсетная. Печать Riso. Печ. л. 1,0. Тираж 100 экз. Заказ 35.

Отпечатано с готового набора в типографии издательства «Вектор Бук». Лицензия ПД № 17-0003 от 06.07.2000 г.

625004, г. Тюмень, ул. Володарского, 45. Тел. (3452) 46-54-04,46-90-03.

Введение.

Глава I Состояние вопроса и задачи исследования.

Экспериментальное изучение кинетики распада аустенита.

Аналитические модели кинетики распада переохлажденного аустенита.

Имитационная модель кинетики распада переохлажденного аустенита.

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

Выводы.

Глава II Векторная модель структурных превращений.

Векторное описание структуры стали при превращении аустенита в одну структурную составляющую.

Векторное описание структуры стали при превращении аустенита в несколько структурных составляющих.

Имитационная модель структурных превращений.

Определение кинетических параметров имитационной модели на основе экспериментальных данных.

Представление кинетических параметров модели.

Выводы.

Глава III Программный комплекс имитационного моделирования распада аустенита.

Технологии и архитектура системы.

Структура базы данных.

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

Экспериментальное исследование структурных превращений в результате термической обработки заключается ■ в построении изотермических и термокинетических диаграмм, а также в изучении прокаливаемости стали [32, 33, 10].

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

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

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

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

В области моделирования кинетики фазовых превращений наиболее фундаментальные исследования были проведены в работах А. Н. Колмогорова и М. Авраами. Модель А. Н. Колмогорова и уравнение М. Авраами , описывают аналитические уравнения, где объем вновь образуемой фазы ставится в зависимость от вероятности зарождения центров, линейной скорости их роста и прошедшего времени. [27, 28].

В отечественной школе фундаментальный обзор математического моделирования процесса фазовых прёвращений был проведен в работе Б. Я. Любова. Наиболее полная зарубежная работа, освящающая все попытки моделировать кинетику фазовых превращений является работа Дж. Кристиана, которая показывает достоинства и недостатки аналитических моделей. В обоих трудах указано, что аналитические модели позволяют прогнозировать кинетику распада аустенита, но хорошо описывают процесс только в области стационарного зарождения, а в начальный период времени инкубационный период не учитывается [114]. По этому, во многих исследованиях моделирование процесса разделялся на две части — нестационарного и стационарного зарождения, что противоречит физике процесса. Также отмечено, что аналитические модели не позволяют рассматривать процесс образования отдельных зерен и их роста.

Помимо аналитических моделей существует подход, заключающийся в имитации процесса зарождения и роста на основе модели И. JI. Миркина, в которой предлагается изучать процесс кристаллизации с помощью двумерных схем [25, 50]. Модель учитывает стохастическую природу процесса, учитывает инкубационный период и позволяет получить размер зерна [25]. В тоже время стохастическая природа образования центров кристаллизации в модели не реализована в полном объеме. В данной модели пространство расплава и кристаллов представляется в виде набора элементарных фрагментов, на подобии растрового изображения в компьютерной графике, что затрудняет реализацию модели в трехмерном пространстве, а также описания процесса в непрерывном пространстве и времени.

В настоящее время проблемами моделирования фазовых превращений в нашей стране занимаются несколько школ, среди которых можно выделить следующие: В.Н.Анциферова [23, 38] Пермского государственного технического университета, Ю. Г. Гуревича [39, 42, 46] Курганского машиностроительного института занимаются моделированием фазовых превращение в порошковых сталях. Группа авторов под управлением Д. А. Мирзаева Южно-Уральского государственного университета занимается моделированием структурных превращений в сталях [115, 116]. В качестве примера моделирования процесса зарождения и роста в специальных сплавах можно привести работы В. Я. Шура Уральского государственного университета им. А. М. Горького, посвященные изучению эволюции сегнетоэлектрической доменной структуры в электрическом поле [117].

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

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

Для достижения поставленной цели определены следующие задачи исследования:

1) создание векторного представления структуры стали в трехмерном пространстве, которое даст возможность описать каждую структурную составляющую в любой момент времени и в любой точке пространства термообрабатываемой детали;

2) создание стохастической модели зарождения и роста структурных составляющих в результате распада аустенита при термической обработке стали;

3) разработка алгоритма определения кинетических параметров векторной имитационной модели распада аустенита по данным кинетических диаграмм полученных в изотермических условиях;

4) разработка комплекса программ, реализующих разработанные модели и алгоритмы, доступ к которым осуществляется посредством web-технологий;

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

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

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

На защиту выносятся:

1) векторная модель структуры стали, полученной в результате распада аустенита при термической обработке стали;

2) векторная имитационная модель распада аустенита;

3) модифицированный метод Хука-Дживса, нахождения кинетических параметров с учетом статистического критерия;

4) архитектура программного комплекса имитационного моделирования распада аустенита.

Научная новизна работы заключается в следующем:

1) создана векторная модель структуры стали, полученной в результате распада аустенита при термической обработке стали;

2) разработана векторная имитационная модель распада аустенита, которая дает возможность реализовать математический аппарат теории массового обслуживания и проводить имитационные эксперименты на основе событийного представления кинетики фазового превращения;

3) обоснован метод получения термокинетических диаграмм исходя из кинетических параметров, рассчитанных на основе изотермических диаграмм.

Практическая значимость работы:

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

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

1) Научные исследования по тематическому плану научно-исследовательских работ №73-01 от 27.11.2006 «Разработка теоретических основ влияния пористости на кинетические параметры распада переохлажденного аустенита». Регистрационный номер НИР 01.20703455.

2) Аналитическая ведомственная целевая программа Министерства образования и науки Российской федерации «Развитие научного потенциала высшей школы (2009-2010)» - «Прогнозирование кинетики распада аустенита в порошковых сталях при непрерывном охлаждении». Регистрационный номер НИР 2.1.2/6498;

Апробация работы. Основные результаты диссертационной работы были представлены на следующих научных конференциях: Всероссийской научно-практическая конференции студентов, аспирантов и молодых ученых «Молодежь и современные информационные технологии» (Россия, Томск, февраль 2007); Научно-методические семинары кафедры информационных систем Тюменского государственного университета (2005-2010 гг.)

Публикации. Основное содержание работы отражено в 13 публикациях, из которых 2 свидетельства о государственной регистрации программы для ЭВМ и 2 статьи опубликованы в журнале из списка ВАК.

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения и списка литературы. Объем диссертации составляет 130 страниц, содержит 24 рисунка, 2 приложения. В библиографии представлено 117 наименований работ российских и зарубежных авторов.

Выводы по диссертационной работе

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

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

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

4) Модифицирован алгоритм Хука-Дживса для получения кинетических параметров имитационной модели фазовых превращений на основе экспериментально построенных изотермических и термокинетических диаграмм.

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

6) Разработана архитектура программного комплекса на основе Web-технологий, которая обеспечивает представление результатов машинного моделирования, а также предоставляет пользователям доступ к сервисам программного комплекса. Пользователи могут проводить моделирование фазовых и структурных превращений через сеть интернет.

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

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

1. Попов А.А., Попова Л.Е. Изотермические и термокинетические диаграммы распада переохлажденного аустенита: Справочник термиста. М.: Металлургия, 1965. - 500 с.

2. Кукушкин Н.Н. Исследование процессов превращения в спеченной стали при нагреве и охлаждении. Дис. . канд. техн. наук. —Л., 1974. - 178 с.

3. Гуревич Ю.Г., Рахманов В.И. Термическая обработка порошковых сталей. -М.: Металлургия, 1985. 80 с.

4. Штейнберг С.С. Избранные статьи. -М.: Машгиз, 1950 563 с.

5. Металловедение и термическая обработка стали. Т.1. Методы испытания и исследования. /Под ред. Бернштейна М.Л., Рахштадта А.Г. -М.: Металлургия, 1983. 352с.

6. Шкляр Р.Ш., Попов А.А. Универсальный магнитометр для изучения быстро и медленно протекающих превращений в ферромагнитных материалах. // Труды Уральского политехнического института. Термическая обработка металлов. Свердловск Москва, 1954, №46. С.34-49.

7. Белоус М.В., Васильев М.А., Черепин В.Т. Магнитометр для быстропротекающих процессов // Заводская лаборатория. 1966. -№3. - С.370-371.

8. Физическое металловедение. /Под ред. Р. У. Кана, П. Хаазева. 3-е изд. Т. 1. Атомное строение металлов и сплавов. Пер. с англ. М.: Металлургия, 1987. 640 с.

9. Панынин И.Ф., Рахманов В.И. Магнитный метод исследования превращения при охлаждении стали. // Научно-исследовательская работа и учебный процесс. Курган, 1969. С. 40-46.

10. Ю.Малинкина Е.И., Ломакин В.Н. Прокаливаемость стали. — М.: Машиностроение, 1969. — 180с.

11. П.Лившиц Б.Г., Крапошин B.C., Линецкий Я.Л. Физические свойства металлов и сплавов. — М.: Металлургия, 1980. — 320с.

12. Гуревич Ю. Г., Рахманов В. И., Паныпин И. Ф. Термокинетические диаграммы превращения аустенита порошковых низколегированных сталей // Порошковые конструкционные материалы. Киев, 1980. С. 146-149.

13. Ивашко А.Г., Юшковский А.Г., Паныпин И.Ф. Магнитометр для исследования распада аустенита в компактных и порошковых материалах. Курган, 1984. -7 с. - Деп. В ВИНИТИ 1984, №3д 12396.

14. Н.Ивашко А.Г. Магнитометр для изучения фазовых превращений в порошковых сталях. Тез. докл. совещания по комплексной программе "Порошковая металлургия". Курган, 1984. С. 9—10.

15. Рахманов В.И., Паныпин И. Ф., Гуревич Ю. Г. Превращения аустенита в спеченных сталях при непрерывном охлаждении. // Исследование технологии металлических порошков и спеченных материалов. Свердловск, 1980. С. 81-85.

16. Ивашко А.Г., Гуревич Ю.Г., Юшковский А.Г., Паныпин И.Ф. Влияние углерода и способа легирования на кинетику распада переохлажденного аустенита // Технология получения изделий из порошков и исследование их свойств. Пенза, 1985. С. 33.

17. Гуревич Ю.Г., Рахманов В.И., Ивашко А.Г., Паныпин И. Ф. Магнитометрическое исследование кинетики фазовых превращений. // Тез. докл. семинара-аукциона: Порошковая металлургия и термообработка. Курган, 1989 С. 12-13.

18. Гуревич Ю.Г., Ивашко А.Г., Белканов И.А. Автоматизированный магнитометр. // Термообработка порошковых сталей: тез. докл.

19. Всесоюз. науч.-техн. конф. / Под ред. Гуревича Ю. Г. Курган, 1991. С. 13-14.

20. Гуревич Ю. Г., Рахманов В. И., Ивашко А. Г., Микуров А. И. Теория и практика термической обработки порошковых сталей. // Новые материалы и технологии в машиностроении. Матер, регион, науч.-техн. конф. Тюмень, 1997. С. 22.

21. Гуревич Ю.Г., Ивашко А.Г. Кинетика распада переохлажденного аустенита порошковых сталей. Курган: Изд-во Курган, гос. ун-та, 1998.- 153 с.

22. Гуревич Ю.Г., Ивашко А.Г., Щегловатов В.В., Колбин Д.М. Магнитометрический метод исследования распада переохлажденного аустенита. // Теория и технология производства новых конструкционных материалов: сб. научных трудов. Курган, 2000. С. 7— 13.

23. Гуревич Ю.Г., Анциферов В.Н., Буланов В .Я., Ивашко А.Г. Термокинетические и изотермические диаграммы порошковых сталей: Справочник / Под. Ред. Ю.Г. Гуревича, Екатеринбург: УрО РАН, 2001, 260 с.

24. Новиков И.И. Теория термической обработки металлов. М.: Металлургия, 1978. - 392 с.

25. Старк Б.В., Миркин И.Л., Романовский А.Н. Металловедение и термическая обработка // Труды Московского института стали, 1935, №7, с. 5-38.

26. Миркин И.Л. Фазовые превращения. // Структура и свойства сталей. — М, Оборонгиз, 1941, с. 5-158.

27. Колмогоров А. Н. К статистической теории кристаллизации металлов // Известия АН СССР. Серия математическая, 1937, № 3, с. 355-358.

28. Avraami М. "Chem. Phys.", 1939,v.7, N8 p. 1103-1112.

29. Курдюмов Г. В., Утевский Л.М., Энтин Р.И. Превращения в железе и стали. -М.: Наука, 1977. -238с.

30. Любов Б. Я. Кинетическая теория фазовых превращений. М.: Металлургия, 1969. —264с.

31. Любов Б. Я. Теория кристаллизации в больших объемах. — М.: Наука, 1975,-256с.

32. Блантер М.Е. Теория термической обработки. М.: Металлургия, 1984.-328с.

33. Блантер М.Е. Фазовые превращения при термической обработке стали. -М.: Металлургиздат, 1962. 268 с.

34. Mokarski S. Nov P/W Alloy Steels for critical components — Metals engineering quarterly, 1973, v. 13, №4, p. 21-27.

35. Математическое описание диаграммы превращения в координатах «время-температура» для изотермического превращения и непрерывного охлаждения. / ВЦП, №305 26с. - Archive fur das Eisenhiittunnvesen, 1974, v. 45, №8, p. 525-532.

36. Кисино Т. и др. Фазовые превращения, теплопроводность и упруго-пластические напряжения при закалке стали. / ВЦП, № В-45714 23 с. - Дзайре, 1979. Т. 28, №312. С.861-867.

37. Анциферов В.Н., Буланов В.Я., Богодухов С.И., Гревнов Л.М. Термохимическая обработка порошковых сталей.— Екатеренбург, УрО РАН, 1997-482с.

38. Анциферов В. Н., Боброва С. Н., Перельман О. М., Шацов А. А. Изотермический распад аустенита порошковой никельмолибденовой стали // МиТОМ. 1993, №8. С/18-20.

39. Гуревич Ю. Г., Ивашко А. Г. Математическая модель процесса закалки сталей с целью прогнозирования структуры. // Нелинейные процессы и проблемы самоорганизации в современном материаловедении: тез. докл. Первого Всерос. семинара. М.: МГУ, 1997. С. 86.

40. Ивашко А. Г., Богословцев JI. А., Стукало В. А. Математическая модель фазовых превращений в порошковых сталях. // Математическое и программное обеспечение научных исследований и обучения: сб. научных трудов. Курган, 1998. С. 46.

41. Гуревич Ю.Г., Ивашко А.Г., Микуров А.И., Богословцев JI.A. Вержбалович Т. А. Методика исследования фазовых превращений. // Международная науч.—техн. • конф. "От булата до современных металлов". Златоуст, 1999.

42. Гуревич Ю.Г., Ивашко А.Г., Богословцев JT.A. Определение кинетических параметров математической модели фазовых превращений. // Математическое и программное обеспечение научных исследований и обучения: сб. научных трудов. Курган, 2000.

43. Ивашко А.Г. Теоретические основы кинетики распада аустенита в порошковых сталях, разработка оптимальных параметров их термообработки: Дис. . доктора техн. наук. Курган, 1998. — 292 с.

44. Апаев Б.Л. Фазовый магнитный анализ сплавов. М.: Металлургия, 1976.-200с.

45. Ивашко А.Г., Богословцев JLA., Цыганова М.С. Статистические методы обработки магнитограмм. // Теория и технология производства новых конструкционных материалов: сб. научных трудов. Курган, 2000. С. 23-27.

46. Гуляев А.П. Металловедение. М.: Металлургия, 1977. - 647 с.

47. Банди Б. Методы оптимизации. Вводный курс. М.: "Радио и связь", 1988, 128 с.

48. Таха X. Введение в исследование операций, 7—е издание. / X. Таха: пер. с англ. М.: Издательский дом «Вильяме», 2005. — 912 е.: ил. - Парал. тит. англ.

49. Анфилатов B.C. Системный анализ в управлении: Учеб. пособие / B.C. Анфилатов, А.А. Емельянов, А.А. Кукушкин; под ред. А.А. Емельянова. М.: Изд—во: «Финансы и статистика», 2002. — 368 е.: ил. • •

50. Миллер Б.М., Панков А. Р. Теория случайных процессов в примерах и задачах. М.: ФИЗМАТЛИТ, 2002. - 320 с.

51. Блюмин С.Л. Дискретное моделирование систем автоматизации и управления: монография / С.Л. Блюмин, A.M. Корнеев; Липецкий эколого—гуманит. ин—т. Липецк: Издательство: «ЛЭГИ», 2005. -124 е.; библиогр. с. 121—123.

52. Википедия свободная энциклопедия: Petri net. Электронный ресурс. URL: http://en. wikipedia. org/wiki/Petrinet

53. Википедия свободная энциклопедия: Petri net. Электронный ресурс. URL: http://ru.wikipedia.org/wiki/CTaTHCTH4ecKHftKpHTepHfi

54. Ивашко А.Г., Полищук И.Н. Имитационное моделирование фазового превращения переохлажденного аустенита в стали. С. 92 — 96 // Математическое и информационное моделирование: Сборник научных трудов. Вып. 8. Тюмень: Издательство «Вектор Бук», 2006.

55. Полищук И.Н. Применение теории массового обслуживания при моделировании распада аустенита С. 145-147 // Информационно-вычислительные технологии и их приложения: сборник статей VI Международной научно-технической конференции. — Пенза: РИО ПГСХА, 2007.

56. Вентцель Е. С. Теория вероятностей: Учебник для студ. Вузов / Елена Сергеевна Вентцель. 10—е изд., стер. - М.: Издательский центр «Академия», 2005. - 576 с.

57. Гмурман В. Е. Теория вероятностей и математическая статистика. — М.: Высшая школа, 1977.

58. Химмельблау Д. Анализ процессов статистическими методами. М.: Мир, 1973.

59. Бусленко Н. П. Моделирование сложных систем. -М.: Наука, 1978.

60. Демидович Б. П., Марон И. А. Основы вычислительной математики. — М.: Наука, 1970.

61. Википедия свободная энциклопедия: Grid computing. Электронный ресурс.

62. URL: http://en.wikipedia.org/wiki/Gridcomputing

63. Рихтер Дж-CLR via С#. Программирование на платформе Microsoft .NET Framework 2.0 на языке С#. Мастер—класс. / Пер. с англ. — М.: Издательство «Русская Редакция» ; СПб.: Питер , 2007. — 656 стр. : ил.

64. Эспозито Д. Microsoft ASP.NET 2.0 Базовый курс. Мастер-Класс / Пер. с англ. М.: Издательство «Русская редакция»; СПб.: Питер, 2007. -688 стр.: ил.

65. Создание приложений Microsoft ASP.NET / Пер. с англ. — М.; Издательско—торговый дом «Русская Редакция», 2002. — 480 е.: ил.

66. ALGLIB PROJECT Библиотека исходных кодов математических алгоритмов Электронный ресурс. URL: http://alglib.sources.ru/interpolation/spline3.php

67. Википедия свободная энциклопедия: GNU General Public License. Электронный ресурс. URL: http://ru.wikipedia.org/wiki/GNUGeneralPublicLicense

68. Нортроп Тони, Уилдермьюс Шон, Райан Билл Основы разработки приложений на платформе Microsoft .NET Framework / Учебный курс Microsoft / Пер. с англ. М.: «Русская редакция», СПб.: «Питер», 2007. - 864 стр.: ил.

69. Арнольд P.P. Расчет и проектирование магнитных систем с постоянными магнитами. — М.; Энергия; 1969. — 184 с.

70. Ананьев Б.Л. Фазовый магнитный анализ сплавов. — М.: Металлургия, 1976. — 200с.

71. Афанасьев Ю.В. Феррозонды. — Л.: Энергия, 1969. — 168 с.

72. Википедия свободная энциклопедия: Система массового обслуживания Электронный ресурс. URL: http://ш.wikipedia.org/wiki/Cиcтeмaмaccoвoгooбcлyживaния

73. Александров П. С. Лекции по аналитической геометрии, пополненные необходимыми сведениями из алгебры с приложением собрания задач, снабженных решениями, составленного А. С. Пархоменко. — М. 1968 г., 912 стр. с илл.

74. Вентцель Е.С. Введение в исследование операций. Рисунки М. А. Герштейна. — М.: «Советское радио», 1964, 391 с.

75. Гнеденко Б. В., Хинчин А. Я. Элементарное введение в теорию вероятностей. Главная редакция физико—математической литературы изд—ва «Наука», М., 1970 г., 198 стр. с илл.

76. Калиткин Н. Н. Численные методы. Главная редакция физико— математической литературы изд—ва «Наука», М., 1978 г., 512 стр. с илл.

77. Хинчин А.Я. Работы по математической теории массового обслуживания, М., Физматгиз, 1963 г., 236 стр.

78. Булинский А. В., Ширяев А. Н. Теория случайных процессов. — М.: Физматлит, 2005. 408 с.

79. Вентцель А. Д. Курс теории случайных процессов. 2—е изд., доп. М.: Наука. Физматлит, 1996, 399 с". '

80. Миллер Р. М., Панков А. Р. Теория случайных процессов в примерах и задачах. М.: Физматлит, 2002. — 320 с.

81. Тутубалин В. Н. Теория вероятностей и случайных процессов: Учеб. Пособие. М.: Изд—во МГУ, 1992 - 400 с.

82. Каханер Д., Моулер К., Нэш С. Численные методы и программное обеспечение: Пер. с англ. Изд. Второе, стереотип. - М.: 2001. -576 е., ил.

83. Дж. Темплман, Д. Виттер. .NET Framework: Библиотека классов. Пер с англ. М.: Кудиц—Образ, 2003. - 672 с.

84. Нортроп Т., Уильдермьюс ILL, Райан Б. Основы разработки приложений на платформе Microsoft .NET Framework. Учебный курс Microsoft / Пер. с англ. М.: «Русская редакция», СПб.: «Питер», 2007. - 864 стр.: ил.

85. Майо Д. С#: Искусство программирования. Энциклопедия программиста: Пер. с англ. / Джозеф Майо СПб.: ООО «ДиаСофтЮП», 2002. - 656 с.

86. Купцевич Ю.Е. Альманах программиста, том 3. Платформа 2003: Microsoft Windows Server 2003, Microsoft Internet Information Services 6.0, Microsoft Office System. — M.: Издательско—торговый дом «Русская Редакция», 2003. — 320 е.: ил.

87. Купцевич Ю. Е. Альманах программиста, том I: Microsoft ADO.NET, Microsoft SQL Server, доступ к данным из приложений — М.: Издательско—торговый дом «Русская Редакция», 2003. — 400 е.: ил.

88. Википедия — свободная энциклопедия: Трехуровневая архитектура Электронный ресурс. URL: http://ru.wikipedia.org/wiki/Tpexypoвнeвaяapxитeктypa

89. Кремер Н.Ш. Теория вероятностей и математическая статистика: Учебник для вузов. 2—е изд., перераб. И доп. - М.: Юнити—Дана, 2004. - 573 с.

90. Корн Г., Корн Т. Справочник по математике (для научных работников и инженеров). Определения, теоремы, формулы. 6—е изд., стер., — СПб.: Издательство «Лань», 2003. 832 с.

91. Полищук И.Н., Ивашко А.Г. Свидетельство о государственной регистрации программы для ЭВМ №2009614999 «Подсистемаимитационного моделирования распада переохлажденного аустенита при изотермических условиях» от 09.10.2009.

92. Полищук И.Н., Ивашко А.Г. Свидетельство о государственной регистрации программы для ЭВМ №2009614000 «Подсистема имитационного моделирования распада аустенита при термокинетических условиях» от 09.10.2009.

93. Большая Советская Энциклопедия: Аустенит Электронный ресурс. URL: http://bse.sci—lib.com/article081760.html

94. Роджерс Д., Адаме Дж. Математические основы машинной графики: Пер. с англ. — М.: Мир, 2001. — 604 е., ил.

95. Самарский А.А. Введение в численные методы. Уч. пособие для вузов. 5—е изд. — СПб.: Издательство «Лань», 2009. — 288 с.

96. Самарский А.А., Гулин А. В. Численные методы: Учеб. Пособие для вузов. — М.: Наука. Гл. ред. физ—мат. лит., 1989. — 432 с.

97. Прокаливаемость сталей /Умэмото М. И др. Нэцу сери, 1980. Т.20. №6. С.317-322. ' ■

98. Качанов Н.Н. Прокаливаемость стали. -М: Металлургия, 1978. -192с.

99. Аксенов Г.И., Заббаров Р.И. Влияние пористости на прокаливаемость металлокерамических материалов. В кн.: Термическая и химико-термическая обработка в порошковой металлургии. Киев: ИПМ АН УССР, 1969. С. 100-112.

100. Фазовые превращения, теплопроводность и упруго-пластические напряжения при закалке стали./Кисино Т. и др. Дзайре, 1979. Т. 28, №312. С.861-867.

101. Кристиан Дж. Теория превращений в металах и сплавах. М.: Мир, 1978. 806 с.

102. Мирзаев Д.А., Окишев К.Ю., Счастливцев В.М., Мирзоев А.А., Яковлева И.Л., Карзунов С.Е. Превращение аустенита в феррит вклассическом" сплаве Fe—9 %Cr. I. Анализ литературных данных. Электронный ресурс.

103. URL: http://www.csc.as.ru/news/19982/2-3-l .pdf

104. Мирзаев Д.А., Окишев К.Ю., Счастливцев В.М., Мирзоев А.А., Яковлева И.Л., Карзунов С.Е. Превращение аустенита в феррит в "классическом" сплаве Fe-9 %Cr. II. Собственные исследования. Электронный ресурс.

105. URL: http://www.csc.as.ru/news/19982/2-3-2.pdf

106. Шур В.Я. Математическое моделирование кинетики фазовых превращений. Методическое пособие для студентов 5 курса. УрГУ. Екатеринбург 2005. Электронный ресурс.

107. URL: http://www.labfer.usu.ru/common/pdf/Phasetransition5kurs.pdf