Введение к работе
АКТУАЛЬНОСТЬ РАБОТЫ. Интенсивное развитие атомной энергетики требует новых конструкционных материалов с улучшенными свойствами. Так, перспектива создания быстрых ядерных реакторов, способных конкурировать на мировом энергетическом рынке (БН-1200), прямо зависит от качества конструкционных материалов, способных выдерживать напряженные условия эксплуатации и обеспечивать глубокое выгорание ядерного топлива. Прежде всего, существует острая потребность в сталях ферритно-мартенситного класса с повышенной жаропрочностью оболочек ТВЭЛ и чехлов тепловыделяющих сборок (ТВС).
Для безопасной эксплуатации ядерных энергетических установок (ЯЭУ) с реакторами на тепловых нейтронах таких, как ВВЭР-1000, необходимо усовершенствование корпусных сталей, устойчивых к радиационному охрупчиванию.
Как правило, оптимальный химический состав конструкционных материалов ЯЭУ определяют экспериментально путем изменения содержания одного из легирующих элементов при фиксировании остальных на определенных уровнях. В случае многокомпонентной системы такая стратегия приводит к недопустимым материальным и временным затратам. Часто требуется одновременно оптимизировать не одну характеристику, а комплекс служебных свойств сталей, что еще более усложняет и удорожает или делает вообще невозможным поиск перспективных материалов экспериментальным способом.
С другой стороны, недостаточная разработанность общей микроскопической теории механических свойств кристаллических тел при пластической деформации не позволяет пока определять из первых принципов эффективные технологии производства конструкционных материалов с заданными свойствами.
В этой связи в современном радиационном материаловедении возникает важнейшая проблема создания таких моделей, которые способствовали бы быстрому проектированию конструкционных материалов, удовлетворяющих требованиям безопасной и экономичной эксплуатации перспективных ЯЭУ.
Выходом из этой сложной ситуации может быть применение математико-статистических методов обработки экспериментальной информации об эволюции физико-механических свойств конструкционных материалов в процессе эксплуатации
4 ЯЭУ. На основе таких моделей становится возможным проведение вьшислительных "экспериментов", имитирующих в главных чертах реальные физико-химические процессы, и получение рекомендаций по оптимизации химических составов и условий термической обработки реакторных сталей и сплавов
Важным стимулом для такого рода исследований послужило резкое увеличение вычислительных мощностей персональных компьютеров и развитие методов математического моделирования. Необходимым условием успешного применения регрессионного анализа является наличие массивов экспериментальных данных о свойствах материалов под облучением. Такие базы данных накоплены в атомной отрасли за десятилетия эксплуатации ядерных реакторов.
Отсюда следует, что в настоящее время сложились условия для создания компьютерных систем, обеспечивающих решение поставленных задач. Подчеркнем, что интерес к разработке имитационных систем, предназначенных для создания конструкционных материалов с заданными свойствами, имеет глобальный характер. Например, с 1998 года реализуется проект испытаний в виртуальном реакторе (Virtual Test Reactor - VTR), в котором принимают участие ведущие страны в области ядерной энергетики, в частности, Франция, Великобритания, США и Япония. Важным преимуществом подобных систем является отсутствие негативных экологических факторов, присущих реальным реакторным экспериментам.
Одним из наиболее употребительных методов анализа является оценивание параметров нелинейных регрессионных моделей. Однако, идентификация существенно нелинейных многофакторных моделей большой размерности встречает цепь трудностей (плохая обусловленность матриц Якоби и Гессе, сильная зависимость от начального приближения и т.п), непреодолимых инструментальными способами. Это приводит к необходимости привлечения диалогового подхода "человек-машина", в котором исследователю отводится роль диспетчера, распределяющего процесс поиска решения между различными методами. При этом возрастает вероятность принятия субъективно неверных стратегий движения к оптимуму, что может приводить к непроизводительным затратам времени и энергии.
Таким образом, сложилась тупиковая ситуация, когда управление высокотехнологичными процессами и получение нового знания требуют максимально
5 точных моделей, а сложившийся аппарат нелинейного регрессионного анализа не может предложить эффективного способа их оценивания.
Данная работа посвящена разрешению этой коллизии путем создания методов на основе применения персональных компьютеров, с широким распространением которых открылась эра быстрых и дешевых вычислений. Особое внимание в работе уделено последовательному внедрению в нелинейное оценивание бутстреп-метода как способа извлечения более полной информации из экспериментальных выборок малого объема, и неиросетевому моделированию, позволяющему синтезировать в единый подход методы случайного и последовательного поиска оптимума.
ЦЕЛЬ РАБОТЫ И ЗАДАЧИ ИССЛЕДОВАНИЯ. Целью работы является создание системы эффективных методов и алгоритмов нелинейного оценивания на основе современных информационных технологий. В связи с этим необходимо решить следующие задачи:
исследовать природу закона Мейера-Нелдела;
на основе этого исследования построить алгоритм бутстреп-идентификации параметров экспоненциальной зависимости;
найти на основе нейросетевого моделирования способ решения общей задачи безусловной оптимизации и, как ее важного случая, нелинейной регрессии;
разработать программное обеспечение, реализующее эти подходы и тестировать их при помощи имитационных экспериментов;
оценить при помощи алгоритма бутстреп-идентификации экспоненциального отклика экспериментальные данные по технологии жидкометаллических теплоносителей в ядерной энергетике;
разработать нейросетевые модели процессов эволюции физико-механических свойств конструкционных материалов ЯЭУ;
разработать план реального эксперимента и получить первые образцы ферритно-мартенситных сталей с 12% содержанием хрома с повышенными прочностными характеристиками.
ЛИЧНЫЙ ВКЛАД АВТОРА. Все алгоритмы, методы и программы,
представленные в работе, разработаны автором лично либо при его прямом участии. НАУЧНАЯ ЗНАЧИМОСТЬ ДИССЕРТАЦИИ заключается в:
проведенном впервые исследовании математико-статистической природы закона Мейера-Нелдела;
построении оригинального алгоритма оценивания параметров аррениусовой зависимости при помощи бутстрепа;
разработке нового метода решения общей задачи безусловной оптимизации и нелинейной регрессии на основе искусственной нейронной сети;
создании нейросетевой модели эволюции механических свойств ферритно-мартенситной стали и расчете оптимального химического состава, обеспечивающего повышенную жаропрочность;
подтверждении справедливости результатов расчета в реальных экспериментах;
проведении имитационных экспериментов по исследованию зависимости максимумов предела прочности от температуры облучения;
разработке нейросетевых моделей охрупчивания корпусных сталей ВВЭР, распухания аустенитных сталей и поведения механических свойств сплава Zr+2.5%Nb.
ПРИКЛАДНАЯ ЦЕННОСТЬ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ. Внедрение
разработанных автором методов, алгоритмов и программ в атомную науку и технику позволили решить ряд важных задач, касающихся технологии жидкометаллических теплоносителей и радиационного материаловедения. Уточнены параметры аррениусовой зависимости растворимости кислорода в жидком свинце и водорода в жидком натрии, жидкометаллических теплоносителях перспективных ЯЭУ. Разработка моделей эволюции механических свойств ферритно-мартенситных сталей с 12% содержанием хрома завершилась первой серией материаловедческих экспериментов, подтвердивших расчетные выводы.
Численное исследование охрупчивания корпусов ВВЭР - явления, определяющего в значительной мере эксплуатационный ресурс тепловых реакторов, позволило наметить пути решения данной проблемы.
Разработанные программы внедрены в ЦИПК и НПО «Технология», что подтверждено актами внедрения.
ДОСТОВЕРНОСТЬ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ. Высокая эффективность разработанных автором нейросетевых алгоритмов нелинейного оценивания, бутстреп-идентификации экспоненциального отклика и предсказания механических свойств новых реакторных сталей подтверждена численными исследованиями и результатами материаловедческих экспериментов.
математико-статистическую интерпретацию природы закона Мейера-Неддела;
алгоритм бутстреп-идентификации экспоненциального отклика, основанного на выводах исследования механизма закона Мейера-Неддела;
эмпирическую зависимость коэффициента самодиффузии железа от температуры всех аллотропических модификаций;
нейросетевой алгоритм безусловной оптимизации и нелинейного оценивания;
нейросетевую модель эволюции прочностных и пластических свойств сталей ферритно-мартенситного класса с 12% содержанием хрома;
результаты экспериментального исследования образцов спрогнозированных сплавов ферритно-мартенситного класса с 12% содержанием хрома с повышенным содержанием марганца;
результаты численных экспериментов по определению положения пиков предела прочности сталей ферритно-мартенситного класса с 12% содержанием хрома в зависимости от температуры облучения.
АПРОБАЦИЯ РАБОТЫ. Результаты работы докладывались на российских и
международных научных семинарах и конференциях:
The second (statistical) FINSOV APL seminar, Joensuu, Finland, 1991.
The third international seminar on simulation of devices and technologies, Obninsk, Russia, 1994.
The international conference on APL: the language and its applications (APL'94), Antwerp, Belgium, 1994
The international conference on Applied programming languages (APL'95), San Antonio, Texas, United States, 1995.
The APL 96 Conference on Designing the Future, Lancaster, UK, 1996.
The International Joint Conference on Neural Networks, Washington, United States 1999.
XV Международная конференция по физике радиационных явлений и радиационному материаловедению. Алушта, Крым, Украина, 2002.
7 Российская конференция по реакторному материаловедению. Димитровград, Россия, 2003.
Семинар координационного научно-технического совета по радиационному материаловедению «Главные результаты исследований ферритно-мартенситных сталей, облученных при высоких повреждающих дозах, пути их совершенствования для оболочек твэлов и чехлов ТВС быстрых реакторов». Москва, 2002.
13 Международное совещание «Радиационная физика твердого тела». Севастополь, Украина, 2003.
Российская научная конференция «Материалы ядерной техники. Радиационная повреждаемость и свойства - теория, моделирование, эксперимент». Туапсе, «Агой», Россия, 2003.
14 Международное совещание «Радиационная физика твердого тела». Севастополь, Украина, 2004.
2 Российская научная конференция «Материалы ядерной техники» (МАЯТ-2). Туапсе, «Агой», Россия, 2005.
18 Международное совещание «Радиационная физика твердого тела». Севастополь, Украина, 2008.
V научно-практическая конференция материаловедческих обществ России. «Цирконий: металлургия, свойства, применение». Ершово, Московская область, Россия, 2008.
20 Международное совещание «Радиационная физика твердого тела». Севастополь, Украина, 2010.
СТРУКТУРА И ОБЪЕМ РАБОТЫ. Диссертация состоит из введения, четырех
