Geum.ru
Авторефераты по всем темам  >>  Авторефераты по техническим специальностям  

УДК 519.28+621.039  На правах рукописи

Образцов Сергей Михайлович

Нейросетевые модели эволюции физико-механических свойств конструкционных материалов ядерных реакторов

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

Автореферат

диссертации на соискание ученой степени

доктора физико-математических наук

Обнинск - 2011

Работа выполнена в Государственном Научном Центре Российской Федерации - Физико-энергетический институт имени А.И. Лейпунского

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

                                                               профессор Рачков В.И.

Официальные оппоненты:

доктор физико-математических наук, профессор                        Коробейников В.В.

доктор технических наук, профессор                                Рисованый В.Д.

доктор физико-математических наук                                 Степанов В.А.

Ведущая организация:

Институт металлургии и материаловедения РАН им. А.А.аБайкова.

Защита состоится 24 июня 2011 года в 1000 часов на заседании диссертационного совета Д201.003.01 при Государственном Научном Центре Российской Федерации - Физико-энергетический институт имени А.И. Лейпунского по адресу: 249033, г. Обнинск, Калужская обл., пл. Бондаренко, 1, конференц-зал Главного корпуса

С диссертацией можно ознакомиться в библиотеке ГН - РФ ЦФизико-энергетический институт имени А.И. Лейпунского по адресу: 249033, г. Обнинск, Калужская обл., пл. Бондаренко, 1.

Автореферат разослан  л _____ апреля 2011 г.

Ученый секретарь диссертационного совета,

доктор технических наук  Прохоров Ю.А.

Общая характеристика работы

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

Для безопасной эксплуатации ядерных энергетических установок (ЯЭУ) с реакторами на тепловых нейтронах таких, как ВВЭР-1000, необходимо усовершенствование корпусных сталей, устойчивых к радиационному охрупчиванию.

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

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

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

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

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

Отсюда следует, что в настоящее время сложились условия для создания компьютерных систем, обеспечивающих решение поставленных задач. Подчеркнем, что интерес к разработке имитационных систем, предназначенных для создания конструкционных материалов с заданными свойствами, имеет глобальный характер. Например, с 1998 года реализуется проект испытаний в виртуальном реакторе (Virtual Test Reactor - VTR), в котором принимают участие ведущие страны в области ядерной энергетики, в частности, Франция, Великобритания, США и Япония. Важным преимуществом подобных систем является отсутствие негативных экологических факторов, присущих реальным реакторным экспериментам.

Одним из наиболее употребительных методов анализа является оценивание параметров нелинейных регрессионных моделей. Однако, идентификация существенно нелинейных многофакторных моделей большой размерности встречает цепь трудностей (плохая обусловленность матриц Якоби и Гессе, сильная зависимость от начального приближения и т.п), непреодолимых инструментальными способами. Это приводит к необходимости привлечения диалогового подхода Учеловек-машинаФ, в котором исследователю отводится роль диспетчера, распределяющего процесс поиска решения между различными методами. При этом возрастает вероятность принятия субъективно неверных стратегий движения к оптимуму, что может приводить к непроизводительным затратам времени и энергии.

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

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

ЦЕЛЬ РАБОТЫ И ЗАДАЧИ ИССЛЕДОВАНИЯ. Целью работы является создание системы эффективных методов и алгоритмов нелинейного оценивания на основе современных информационных технологий. В связи с этим необходимо решить следующие задачи:

  • исследовать природу закона Мейера-Нелдела;
  • на основе этого исследования построить алгоритм бутстреп-идентификации параметров экспоненциальной зависимости;
  • найти на основе нейросетевого моделирования способ решения общей задачи безусловной оптимизации и, как ее важного случая, нелинейной регрессии;
  • разработать программное обеспечение, реализующее эти подходы и тестировать их при помощи имитационных экспериментов;
  • оценить при помощи алгоритма бутстреп-идентификации экспоненциального отклика экспериментальные данные по технологии жидкометаллических теплоносителей в ядерной энергетике;
  • разработать нейросетевые модели процессов эволюции физико-механических свойств конструкционных материалов ЯЭУ;
  • разработать план реального эксперимента и получить первые образцы ферритно-мартенситных сталей с 12% содержанием хрома с повышенными прочностными характеристиками.

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

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

  • проведенном впервые исследовании математико-статистической природы закона Мейера-Нелдела;
  • построении оригинального алгоритма оценивания параметров аррениусовой зависимости при помощи бутстрепа;
  • разработке нового метода решения общей задачи безусловной оптимизации и нелинейной регрессии на основе искусственной нейронной сети;
  • создании нейросетевой модели эволюции механических свойств ферритно-мартенситной стали и расчете оптимального химического состава, обеспечивающего повышенную жаропрочность;
  • подтверждении справедливости результатов расчета в реальных экспериментах;
  • проведении имитационных экспериментов по исследованию зависимости максимумов предела прочности от температуры облучения;
  • разработке нейросетевых моделей охрупчивания корпусных сталей ВВЭР, распухания аустенитных сталей и поведения механических свойств сплава Zr+2.5%аNb.

ПРИКЛАДНАЯ ЦЕННОСТЬ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ. Внедрение разработанных автором методов, алгоритмов и программ в атомную науку и технику позволили решить ряд важных задач, касающихся технологии жидкометаллических теплоносителей и радиационного материаловедения. Уточнены параметры аррениусовой зависимости растворимости кислорода в жидком свинце и водорода в жидком натрии, жидкометаллических теплоносителях перспективных ЯЭУ. Разработка моделей эволюции механических свойств ферритно-мартенситных сталей с 12% содержанием хрома завершилась первой серией материаловедческих экспериментов, подтвердивших расчетные выводы.

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

Разработанные программы внедрены в ЦИПК и НПО Технология, что подтверждено актами внедрения.

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

Автор выносит на защиту:

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

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

  1. The second (statistical) FINSOV APL seminar, Joensuu, Finland, 1991.
  2. The third international seminar on simulation of devices and technologies, Obninsk, Russia, 1994.
  3. The international conference on APL: the language and its applications (APLТ94), Antwerp, Belgium, 1994
  4. The international conference on Applied programming languages (APLТ95),а San Antonio, Texas, United States, 1995.
  5. The APL 96 Conference on Designing the Future, Lancaster, UK, 1996.
  6. The International Joint Conference on Neural Networks, Washington, United States 1999.
  7. XV Международная конференция по физике радиационных явлений и радиационному материаловедению. Алушта, Крым, Украина, 2002.
  8. 7 Российская конференция по реакторному материаловедению. Димитровград, Россия, 2003.
  9. Семинар координационного научно-технического совета по радиационному материаловедению Главные результаты исследований ферритно-мартенситных сталей, облученных при высоких повреждающих дозах, пути их совершенствования для оболочек твэлов и чехлов ТВС быстрых реакторов. Москва, 2002.
  10. 13 Международное совещание Радиационная физика твердого тела. Севастополь, Украина, 2003.
  11. Российская научная конференция Материалы ядерной техники. Радиационная повреждаемость и свойства - теория, моделирование, эксперимент.аТуапсе, Агой, Россия, 2003.
  12. 14 Международное совещание Радиационная физика твердого тела. Севастополь, Украина, 2004.
  13. 2 Российская научная конференция Материалы ядерной техники (МАЯТ-2).аТуапсе, Агой, Россия, 2005.
  14. 18 Международное совещание Радиационная физика твердого тела. Севастополь, Украина, 2008.
  15. Vанаучно-практическая конференция материаловедческих обществ России. Цирконий: металлургия, свойства, применение. Ершово, Московская область, Россия, 2008.
  16. 20 Международное совещание Радиационная физика твердого тела. Севастополь, Украина, 2010.

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения, списка использованной литературы из 228 наименований и включает 58 рис. и 26 табл. Общий объем 184 стр.

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

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

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

.        (1)

Здесь, B= (-E/k), Е - энергия активации, эВ, k Цпостоянная Больцмана; А - предэкспоненциальный коэффициент, Т - температура, К, С - зависимая переменная.

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

Матрица J для (1) имеет вид:

,       (2)

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

При явной корреляции между столбцами в (2) имеем

,

что можно представить в виде обыкновенного дифференциального уравнения:

,        (3)

после интегрирования которого получим искомую зависимость:

,        (4)

где k1 - постоянная интегрирования, значение k2 приближенно равно середине температурного интервала.

Учитывая, что k2 выполняет роль тангенса угла наклона линии (4), автором диссертационной работы предложен принципиально новый алгоритм расчета бутстреп-оценки (1), как точки пересечения прямых (4), построенных для различных усеченных последовательностей исходной упорядоченной выборки. Для ее расчета необходимо:

  1. разделить исходный температурный интервал на несколько произвольных отрезков;
  2. провести бутстреп-процедуру, получить линейные оценки для каждой частичной выборки данных;
  3. рассчитать бутстреп - оценки для обратной задачи (4).

В итоге получается , где K1 и - множество частных оценок параметров (4), рассчитанных в пункте 2 приведенного выше алгоритма.

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

Dа=аD0exp(-E/kT),         (5)

где D0 - предэкспоненциальный множитель, Е - энергия активации самодиффузии, эВ, k - постоянная Больцмана.

На рис. 1 приведены оценки параметров (5), полученные разными авторами в разное время для - железа. Видно, что полученные независимым образом оценки можно описать одной прямой (4), а разброс значений энергии активации достигает ~1.4аэВ.

Для устранения неопределенности в оценках энергии активации автор использовал разработанный алгоритм бутстреп-идентификации экспоненты по экспериментальным данным о параметрах температурной зависимости коэффициента самодиффузии для всех аллотропических модификаций (α-, γ- и δ-) железа.

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

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

На основании пяти бутстреп - оценок для парамагнитной области были рассчитаны коэффициенты уравнения (4), которое записывается в виде:

lnаDа=а-36.05+11.06Е.         (6)

С учетом бутстреп-оценки Еа=а2.59аа0.2аэВ, разброс значений D0 определяется уравнениема(6).

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

,        (7)

где q1, q2, q3 и Е - подгоночные константы.

В результате идентификации отклика (7) расчета получены следующие значения констант q1а=а2.8, q2а=а972аК, q3а=а147.5аK, Eа=а2.63аэВ.

Для определения D0 и Е в -фазе использовалась комбинированная выборка из литературных источников, которая обрабатывалась c помощью процедуры бутстреп-идентификации экспоненты. В результате расчета были получены оценки D0а=а0.0000258ам2/с, Еа=а2.84аа0.01аэВ и связывающее их уравнение lnаDа=а-29.2+6.55Е.

Аналогично, на основе выборки данных были рассчитаны параметры коэффициента самодиффузии для Цжелеза: D0а=а0.000644ам/с, Еа=а2.63а0.004аэВ. Уравнение, связывающее D0 и Е, имеет следующий вид: lnаDа=а-21.4+5.34Е.

На рис.а2 приведены экспериментальные и вычисленные значения коэффициента самодиффузии для α-, γ-, и δ-фаз железа.

Рис.а2 - Расчетные и экспериментальные значения D

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

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

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

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

Каждый из нейронов реализует функцию вида:

,         (8)

где p - параметр, подлежащий определению.

Для безусловной минимизации целевой функции

        (9)

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

В результате проведенного исследования автором предложен нейросетевой алгоритм безусловной оптимизации и нелинейного оценивания нулевого порядка. Далее описана схема разработанной ИНС для решения (11).

Входной слой (S - слой) состоит из n узлов, соответствующих размерности искомого вектора , столько же узлов содержит второй скрытый слой (A2). Выходной слой (R) содержит один нейрон, количество узлов в первом скрытом слое (A1) назначает пользователь.

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

, ,         (10)

где W - матрица весов между соответствующими слоями, RAND(I) означает процедуру равномерного случайного выбора в интервале I, - стартовая точка. После этого начальный набор (10) корректируется как Wk+1=Wk+W, где . Здесь G обозначает нормальное распределение, 1 и 2 являются небольшими положительными числами, скаляр получатся в результате одномерной минимизации интерполирующей функции:

,        (11)

где - некоторый вектор, подлежащий определению.

Входными сигналами являются компоненты искомого вектора, а выходной нейрон рассчитывает . Алгоритм не имеет ограничений на область действия, поскольку всегда можно ввести новую функцию , где L>0 является константой. Минимизация F1 приводит к минимуму F.

Алгоритм тестировался на классических задачах Розенброка, Пауэлла и Вуда. Точное решение этих задач достигается при значениях аргумента: (1,1); (0,0,0,0); (1,1,1,1), соответственно. Каждая из этих задач решалась 65 раз, после чего рассчитывалось среднее число NAVER вычислений целевой функции и среднее значение полученной оценки . Использовалась сеть с числом нейронов во втором скрытом слое, равным 7. Условиями останова являлись уменьшение целевой функции до определённого значения или относительное замедление движения по каждой координате до 0.05%. Результаты представлены в табл.а1, где в третьем столбце приведены средние значения искомых переменных.

Таблица 1

Результаты тестирования алгоритма

Задача

Розенброка

1100

0.9995;0.9999

Вуда

54000

1.0006;1.001;0.9994;0.9989

Пауэлла

330

0.002;-0.0002;0.0005;0.0005

Отметим, что методом Нелдера-Мида не удаётся решить ни одной из приведенных тестовых задач. Метод Флетчера-Ривса (первого порядка) позволяет решить только задачу Розенброка, используя 250 вычислений цели и градиента.

Задача Вуда является наиболее трудной для решения, она имеет вид:

  + .        

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

в зависимости от номера итерации . На рис.а3 представлен график s(lgK), иллюстрирующий особенности сходимости к решению.

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

Из рис.а3 видно, что метод имеет устойчивую сходимость к решению. Повторенный несколько раз расчет коэффициентов функции Вуда показал, что траектория s(lgK) в среднем сохраняет представленную на рис.а3 форму. Это говорит о стабильности работы алгоритма.

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

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

,        (12)

где η(x,) - математическое ожидание Y при заданных независимых переменных x, ξ - случайная величина, имеющая плотность распределения вероятности f(ξ),- k-мерный вектор искомых параметров. Для оценивания используют метод наименьших квадратов в его нелинейном варианте:

.        (13)

В связи со сложностью решения (13) целесообразно применить нейросетевой метод безусловной оптимизации (10)-(11).

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

.        (14)

Решение (13) для этого отклика имеет следующий вид S(115.2; 2.31; -22.022)=7.0133.

После 60-кратного расчета оценок из стартовой точки (-10, 100, -100), было обнаружено, что метод находит решение (115.16, 2.309, -22.00) в среднем за 45730 вычислений суммы квадратов отклонений. Для сравнения эта же задача решалась широко применяющимся методом Марквардта (первого порядка), который не привел в точку оптимума из этого стартового приближения ни при каких значениях регуляризирующего параметра.

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

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

Итак, пользователь задает число узлов во входном и первом скрытом слоях, соответственно, . Узлы входного слоя генерируют случайный вектор, принадлежащий равномерному распределению в интервале [0,1]. Второй скрытый слой содержит узлов, соответствующих размерности искомого вектора . Выходной слой содержит один нейрон, рассчитывающий сумму квадратов отклонений.

Автор тестировал модифицированный алгоритм при помощи отклика:

.        (15)

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

Таблица 2

Результаты оценивания (15) из стартовой точки

S,10-2

0.752

-0.152

1.09

3.31

0.897

0.755

-0.083

1.11

3.51

0.896

0.765

0.203

1.18

4.30

0.895

0.753

-0.129

1.10

3.37

0.896

Предварительные результаты позволяют сделать предположение о существовании множества локальных минимумов в случае отклика (15). Для их определения изменим схему расчета, задавая случайным образом стартовую точку из интервала [-40,40]. За ограниченное время работы программы было найдено 30 локальных минимумов. Для иллюстрации часть множества оценок представлена в табл.а3 по возрастанию суммы квадратов отклонений.

Таблица 3

Результаты оценивания (15) при случайном выборе стартового приближения

S,10-2

13.76

2.946

-1.21

2.27

0.89520

15.5

3.408

-1.21

2.62

0.89521

0.7741

0.522

1.25

5.16

0.89531

0.7751

0.568

1.26

5.28

0.89538

0.7760

0610

1.27

5.39

0.89545

0.7625

0.118

1.16

4.07

0.89546

44.1

11.04

-1.16

8.41

0.89553

84.01

21.7

-1.15

16.5

0.89567

254.4

67.2

-1.14

51.1

0.89578

Из табл.а3 видно, что при самом широком разбросе оценок сумма квадратов отклонений имеет практически одно и то же значение и для отсеивания несостоятельных оценок требуется анализ содержательной информации об объекте. Этот пример убедительно показывает, что формальная (только на основе критерия минимума невязки) идентификация приводит к многозначности решения задачи оценивания уравнения (15). Очевидно, в этом случае требуется репараметризация отклика для однозначной интерпретации решения. Поиск хорошо обусловленной физической модели может потребовать дополнительных затрат дефицитного времени, которые в случае очень большой размерности факторного пространства могут и не оправдаться.

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

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

Таблица 4

Точность приближения в зависимости от структуры ИНС

Номер по порядку

A1

A2

S,10-2

1

2

2

0.965

2

2

5

0.967

3

5

2

0.873

4

7

7

0.421

5

7

10

0.851

6

10

7

0.512

7

10

20

0.515

8

20

10

0.494

9

20

20

0.497

10

30

30

0.506

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

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

В диссертационной работе проведено исследование возможности создания алгоритма нейросетевой оптимизации с использованием бутстреп-метода. Результаты тестирования гибридного алгоритма (ИНС+бутстреп-метод) говорят о перспективности этого подхода.

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

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

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

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

В диссертационной работе описаны результаты разработки нейросетевых моделей изменения прочности и пластичности таких сталей. Однако наибольший интерес представляет модель одновременного поведения пластических и прочностных свойств. В качестве отклика были выбраны общее относительное удлинение δоа(%) и предел прочности σва(МПа).

Матрица независимых переменных состоит из 348 строк, каждая из которых содержит информацию об условиях испытания и измерения кратковременных механических свойств цилиндрических и плоских образцов сталей ЭИ-852, ЭП-450 и ЭП-823 как в исходном, так и в облученном состояниях. В число факторов, лимитирующих изменение δо и σв, включены температура облучения и доза облучения, температура испытания, полное содержание легирующих элементов и т.п.

Суммарное число столбцов матрицы независимых переменных составило 48, что соответствует числу узлов входного слоя ИНС, выходной слой включает два нейрона. Число узлов в двух скрытых слоях подбиралось методом проб и ошибок и составило 20 и 45, соответственно. На рис.а4 значками представлены экспериментальные данные испытания образцов стали ЭП-450 (закалка: 1050oC, 30амин; отпуск 720oС, 1ачас) облученной при температуре Тобл.=350С, до Dа=а0.52асна и результаты расчета (кривые) по нейросетевой модели.

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

Рис.а4 - Экспериментальные значения (значки) и расчетные зависимости (кривые) предела прочности и общего относительного удлинения облученной сталиаЭПа-а450

Так как сталь ЭП-450 применяется в качестве материала чехлов ТВС, то именно ее химический состав был взят за основу при поиске компьютерным способом оптимального сочетания легирующих элементов. Начальное содержание компонентов стали (в вес.%):

,         (16)

а исходная термообработка для этой стали: нормализация 1050С, 30амин.; отпуск 720С, 1ачас.

Целью вычислительного эксперимента было выяснение степени влияния различных комбинаций некоторых легирующих элементов на кратковременные прочностные и пластические свойства. Как показали предварительные расчеты, наибольшее воздействие на повышение σв и δо оказывают согласованные вариации содержания Mn, Mo и Nb, а также B и V.

План компьютерного эксперимента можно представить следующим образом.

1. Критериями оптимизации выбраны суммарные значения предела прочности, МПа:

        (17)

и общего относительного удлинения, %:

,         (18)

где T - температура испытания.

Каждое значение σв и δо рассчитываются при фиксированных значениях температуры испытания Таi=а20,а100,а200,а300,а400,а500,а600 и 710С. Повреждающая доза полагалась равной 100асна. Принимая во внимание, что в быстрых реакторах температура облучения оболочки твэлов может достигать 710С, именно это значение было заложено в расчет суммы (17). Поскольку ферритно-мартенситные стали обнаруживают низкотемпературное радиационное охрупчивание, то температура облучения при вычислении суммы (18) принималась равной 300С.

2. В качестве базового химического состава был выбран состав стали (16).

3. Приоритетным критерием рассматривалась величина S, которая максимизировалась путем расчета по нейросетевой модели оптимальных содержаний легирующих элементов. На величину Е было наложено ограничение, согласно которому она не должна уменьшаться более чем на 5%.

4. Эксперимент проводился в два этапа. На первом этапе содержание Mn задавалось в диапазоне 0.2 6 вес.%. В процессе оптимизации содержания Mo и  Nb оставались свободными и выбирались такими, чтобы обеспечить максимальную прочность во всем диапазоне температур испытания.

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

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

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

СMo = 6.28 - 1.04СMn         (19)

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

СNb = 0.546 - 0.056 СMn + 0.007С2Mn       (20)

Для определения адекватности нейросетевой модели решающее значение имеет экспериментальное исследование влияния марганца на прочностные свойства сталей ферритно-мартенситного класса с 12% содержанием хрома. В этой связи были выплавлены опытные сплавы с различным содержанием марганца и испытаны на прочность в НИТУ МИСиС под руководством профессора БеломытцевааМ.Ю. Содержание молибдена и ниобия изменялось в соответствии с (19) и (20).

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

Таблица 5

Содержание марганца и структурные характеристики экспериментальных сплавов

№ сплава

1

2

3

4

Содержание Mn, вес%

0.265

1.3

1.98

0.30

Доля мартенсита, %

90

40

100

50

Размер зерна мартенсита, мкм

5 - 12

150

70

8 - 20

Из табл.а5 видно, что сплавы сильно различаются по содержанию мартенситной фазы, достигающей 100% в сплаве 3, и по размерам зерен мартенсита. Для унификации результатов измерений было использовано соотношение Холла-Петча, позволяющее единым образом описать предел текучести сталей с зерном различной степени дисперсности: *0,2= 0,2 d, d - размер зерна, *0,2 - константа Холла-Петча. На рис.а5 приведены значения константы Холла-Петча в зависимости от температуры испытания Т.

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

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

       (21)

       (22)

Здесь в - предел прочности, индекс i обозначает номер горба распределения; t - температуру облучения, при которой наблюдается максимум; timax среднюю температуру этого горба; γi характеризует ширину горба; А - нормировочный коэффициент; k - число горбов в распределении. На рис.а6 приведено имитированное распределение температур облучения и аппроксимации для стали ЭП-450.

Рис.а6 - Имитированное распределение температур облучения, при которых предел прочности стали ЭП-450 максимален, и кривые, описывающие это распределение: 1 - уравнение (21), 2 - уравнение (22)

На рис.а7 приведена аналогичная зависимость для стали ЭП-823.

Результаты проведенного имитационного эксперимента, а также обработка данных при помощи сумм колоколообразных функций, указывают на существование у стали ЭП-450 двух пиков предела прочности, расположенных при температурах облучения ~360 и 390С. В случае стали ЭП-823 уверенно можно говорить о наличии максимума σв при температуре облучения 390С и с меньшей вероятностью можно ожидать максимума предела прочности при температуре 360С. Сравнение гистограмм распределения пиковых температур, изображенных на рисунках 6 и 7, показывает, что различное легирование приводит к изменению расположения пиковых температур.

Рис.а7 - Имитированное распределение температур облучения максимальных значений предела прочности стали ЭП-823 и кривые, описывающие это распределение: 1 - уравнение (21), 2 - уравнение (22)

В четвертой главе представлены результаты разработки нейросетевых моделей изменения физико-механических свойств чистого циркония и циркониевого сплава Zr+2,5% Nb, корпусных сталей ВВЭР и сталей аустенитного класса.

Температура испытания оказывает сильное влияние на характер зависимости предела текучести сплава Zr+2,5% Nb от коэффициента стехиометрии гидрида циркония. На рис.а8 представлены расчетные кривые, показывающие эту взаимосвязь при различных температурах испытания.

Кроме того, в результате расчетов было обнаружено влияние содержания ниобия на расположение пика предела текучести сплава Zr+2,5% Nb на оси коэффициента стехиометрии гидрида циркония. Повышение содержания ниобия сдвигает расположение пика предела текучести вправо.

Рис.а8 - Расчетные зависимости предела текучести сплава Zr+2,5% Nb от коэффициента стехиометрии гидрида циркония при различных температурах испытания

В диссертационной работе при помощи ИНС исследовалась зависимость радиационного распухания S аустенитных сталей от повреждающей дозы D. Экспериментальные данные, на основе которых строилась модель, характеризуются следующим рядом параметров: температура T и доза облучения, химический состав (до 19 химических элементов); температура и время аустенизации; температура и время отжига; степень холодной деформации; способ измерения S (гидростатическое взвешивание или электронная микроскопия). Дополнительно в обучающую матрицу был включен функционал, описывающий кооперативное воздействие T и D на S.

Модель представляет собой четырехслойную последовательную нейронную сеть, конфигурация которой выглядит как 29а:а11а:а3а:а1. Адекватность модели проверялась на данных, не представленных в обучающей выборке. На основе модели было проведено исследование влияния Si и Ti на дозовую и температурную зависимости распухания.

В качестве базового химического состава была выбран следующий: [С]=0.07; [Mn]=0.39; [Si]=0.4; [Cr]=15.0; [Ni]=14.75; [Ti]=0; [Nb]=0.18; [Mo]=3.0, где квадратные скобки обозначают содержание соответствующего химического элемента в вес.а%. На рис.а9 представлены температурные зависимости распухания при дозе 55 сна при различном содержании Si и Ti. Кривая 2 соответствует базовому составу; кривая 1 получена при [Si]=0.1; кривая 3 - [Si]=0.8. Кривая 4 рассчитана при УдобавленииФ в базовый химический состав 0.3 вес.а% Ti.

Рис.а9 - Температурные зависимости распухания при разных химических составах

       Из рис.а9 видно, что увеличение содержания кремния и титана уменьшает радиационное распухание.

С точки зрения безопасной эксплуатации ВВЭР интерес вызывает нейросетевой анализ влияния легирующих элементов на радиационное охрупчивание материалов корпусов. Для исследования этого явления была разработана нейросетевая модель, имеющая следующую конфигурацию: входной слой состоит из 14 узлов, выходной - из одного, ассоциативные слои включают 10 и 20 нейронов соответственно.

На основе модели были проведены расчеты влияния легирующих элементов на радиационное охрупчивание. Суть вычислительного эксперимента заключалась в расчете сдвига температуры вязкохрупкого перехода (ТВХП) как функции флюенса нейтронов при заданном химическом составе корпусной стали. Характер и степень влияния химических элементов на радиационное охрупчивание выявлялась путем поочередного увеличения вдвое их содержания относительно базового, в качестве которого был выбран основной металл корпуса 3-го блока Кольской АЭС, вес.а%: С 0,16; Si 0,23; Mn 0,41; Cr 2,63; Ni 0,15; Mo 0,66; Cu 0,09; S 0,01; P 0,011; V 0,28; As 0,013; Co 0,012.

На рис.а10 приведены расчетные кривые, полученные для базового химического состава кривая (1) и увеличенном в 2 раза содержании Cu, Ni, V, Mo, Si (2,а3,а4,а5,а6). При увеличении содержания каждого из перечисленных элементов содержание остальных оставалось неизменным. Для сравнения представлена нормативная кривая (7) ΔТк =18F1/3, (F - флюенс в единицах 1019асм-2, Е>0.5 МэВ) для основного металла корпусов ВВЭР-440.

Рис.а10 - Расчетная зависимость сдвига температуры вязкохрупкого перехода от флюенса нейтронов при варьировании содержания отдельных легирующих элементов: базовый химический состав (1), Cu (2), Ni (3), V (4), Mo (5), Si (6) и нормативная криваяа(7)

Из рис.а10 видно, что температура вязкохрупкого перехода резко возрастает в диапазоне флюенса F от 0 до 71019 см-2, затем следует почти линейный рост до ~ 651019 см-2 и начало выхода на плато в районе 1001019асм-2. Увеличение содержания никеля и меди ухудшает пластические свойства под облучением, а ванадий, молибден и кремний оказывают положительное воздействие на радиационную стойкость корпусной стали. Варьирование содержания марганца в расчетах показало, что он не оказывает заметного воздействия на величину ТВХП.

Из анализа рис.а10 следует, что существует возможность компьютерного расчета оптимального химического состава корпусной стали с помощью нейросетевой модели, разработанной с учетом эксплуатационных факторов и на основе выборки большего объема. Для грубой оптимизации химического состава автор воспользовался тем, что, исходя из расчетов, в состав конструкционного материала входят химические элементы (Si,аMo,аV), компенсирующие отрицательные последствия увеличения содержания никеля до 1а% по массе. В качестве варианта оптимизированного химического состава предложим следующий: С 0,16; Si 0,5; Mn 0,41; Cr 2,63; Ni 1,0; Mo 1,5; Cu 0,09; S 0,01; P 0,011; V 1,6; As 0,013; Co 0,012 вес.%. Для сравнения на рис.а11 изображена кривая при увеличенном содержании Ni (1 вес.%) в Убазовом Ф химическом составе.

Рис.а11 - Расчетная зависимость сдвига температуры вязкохрупкого перехода от флюенса нейтронов при оптимизированном (1), базовом химическом составе (2) и увеличенном содержании никеля (3)

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

ΔТк ∝ 14,9F 0,64.         (23)

Показатель степени при флюенсе в (23) говорит о том, что эта зависимость является более консервативной, чем нормативная.

Исследовалась также зависимость сдвига ТВХП от приращения предела текучести по данным испытания образцов корпусных сталей ВВЭР-440 и ВВЭР-1000. В результате исследования выяснилось, что влияние Ni на охрупчивание нельзя рассматривать в отрыве от содержания вредных примесей, так из литературных данных известно, что в очень чистых экспериментальных сплавах корпусной стали ВВЭР никель даже пластифицирует материал. Напротив, при суммарном содержании P и Sb выше 0.01 вес.% увеличение содержания Ni приводит к резкому охрупчиванию корпусных сталей.

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

ВЫВОДЫ

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

К основным достижениям диссертационной работы можно отнести следующее:

  1. Впервые проведено исследование, позволившее выявить математико-статистическую природу закона Мейера-Нелдела. На основе этого исследования и бутстрепа разработан принципиально новый прецизионный метод оценивания параметров экспоненциального отклика, реализованный в виде компьютерной программы. Проведены имитационные эксперименты, показавшие корректность предлагаемого метода.
  2. При помощи программного комплекса бутстреп-идентификации экспоненты уточнены оценки параметров политерм растворимости водорода в натрии и кислорода в свинце, играющие важную роль в технологии жидкометаллических теплоносителей.
  3. Применение программного комплекса бутстреп-идентификации экспоненты позволило устранить неопределенность в оценках параметров температурной зависимости самодиффузии для всех полиморфных модификаций чистого железа (754÷1788К). Предложена эмпирическая зависимость коэффициента самодиффузии в чистом железе, что крайне важно для интерпретации и моделирования радиационного повреждения реакторных сталей.
  4. Исходя из способности искусственной нейронной сети вырабатывать ассоциативные связи входного и выходного сигнала в процессе обучения, предложен новый подход к решению задачи безусловной оптимизации и нелинейного оценивания. Отличительным признаком этого подхода от известных численных процедур является синтез случайного и детерминированного методов поиска. В зависимости от способа обучения и обрабатываемой информации разработаны следующие алгоритмы (а) первого порядка, оперирующий в процессе расчета первыми производными функции отклика, (б) нулевого порядка, использующий расчет только отклика. Тестирование алгоритмов на модельных задачах безусловной оптимизации и реальных задачах нелинейной регрессии показало их высокую эффективность.
  5. Исследовано бутстреп-расширение нейросетевого подхода. Полученные результаты говорят о перспективности включения бутстрепа в процедуру обучения искусственных нейронных сетей и необходимости продолжения работы с целью повышения эффективности этого алгоритма.
  6. Впервые разработаны и тестированы нейросетевые модели изменения под облучением предела прочности и общего относительного удлинения при нейтронном облучении прочности и пластичности ферритно-мартенситных сталей с 12% содержанием хрома. С помощью разработанной обобщенной нейросетевой модели предсказан оптимальный химический состав, обеспечивающий повышенную жаропрочность и пластичность ферритно-мартенситных сталей с 12% содержанием хрома. Подтверждена важная роль марганца в химических составах сталей с заданными служебными характеристиками. Найдены регрессионные зависимости между содержаниями марганца, молибдена и ниобия.
  7. Проведена экспериментальная проверка нейросетевых прогнозов, подтвердившая их корректность.
  8. При помощи обобщенной нейросетевой модели проведено бутстреп - исследование максимумов предела прочности ферритно-мартенситных сталей с 12% содержанием хрома в зависимости от температуры облучения. Определены пики предела прочности для чистого железа и сталей марок ЭП-450, ЭП-823.
  9. Разработана и тестирована нейросетевые модели охрупчивания корпусных сталей ВВЭР - 440 и зависимости сдвига температуры вязко-хрупкого корпусных сталей ВВЭР - 440 и ВВЭР-1000 от приращения предела текучести. Компьютерные эксперименты на основе этих моделей показали положительное влияние молибдена, кремния и ванадия на устойчивость корпусных сталей к охрупчванию и резко негативное влияние суммарного содержания примесей фосфора, олова и сурьмы при повышенном содержании никеля.
  10. Проведено нейросетевое исследование зависимости предела текучести сплава Zr+2.5% Nb от коэффициента стехиометрии гидрида циркония, обнаружено, что эта зависимость изменяется при варьировании содержания ниобия.
  11. Разработан первый вариант нейросетевой модели распухания сталей аустенитного класса под воздействием нейтронного облучения. Имитационные эксперименты при помощи этой модели показали положительное влияние кремния и титана на радиационную стойкость сталей данного типа.

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

  1. Louneva I.N., ObraztsovаS.M., Shimkevich A.L. Solubility curve bootstrap>
  2. ОбразцоваС.М, ШимкевичаА.Л. Бутстреп - идентификация экспоненциальной зависимости // Заводская лаборатория. - 2000. - Т. 66. - № 1. - С. 62-64.
  3. КуприяноваВ.М, ЛуневааИ.Н., ОбразцоваС.М., СоловьеваВ.А., ШимкевичаА.Л.. Статистический анализ данных растворимости кислорода в расплаве свинца. Препринт ФЭИ-2295. ОНТИ ФЭИ. Обнинск, 1992.
  4. ОбразцоваС.М., ШимкевичаА.Л. Нелинейное оценивание растворимости кислорода в свинце // Сб.: Теплофизические исследования / Обнинск: ГН - РФ-ФЭИ. - 1999. - C. 225.
  5. ObraztsovаS.M., ShimkevichаA.L. A zero order minimization algorithm using a neural network //аProceedings of the International Joint Conference on Neural Networks, Washington. July 10-16. - 1999. - V. 1. - 1999. - P. 643-645.
  6. КоноваловаА.В., Образцов аС.М., ШимкевичаА.Л. Нейросетевой алгоритм безусловной оптимизации. Препринт ФЭИ-2684. ОНТИ ГН - РФ ФЭИ. Обнинск, 1998.
  7. ОбразцоваС.М., БиржевойаГ.А., КонобееваЮ.В., СоловьеваВ.А. и РачковаВ.И. Нейросетевая модификация стали ферритно-мартенситного класса ЭП-450 по критерию максимума прочности и пластичности // Перспективные материалы. - 2005. - №4. - C. 14-19.
  8. ОбразцоваС.М., БиржевойаГ.А., СоловьеваВ.А., РачковаВ.И. Нейросетевое моделирование прочностных свойств реакторных сталей // Труды XV Международной конференции по физике радиационных явлений и радиационному материаловедению. Алушта, 10-15 июня 2002. - C. 105-106.
  9. ОбразцоваС.М., БиржевойаГ.А., КонобееваЮ.В., СоловьеваВ.А., РачковаВ.И. Нейроносетевой анализ кратковременных механических свойств сталей ферритно-мартенситного класса // Физика и химия обработки материалов. - 2004. - № 2. - С. 92-97.
  10. ОбразцоваС.М., БиржевойаГ.А., СоловьеваВ.А. Нейросетевая модель поведения пластических свойств 12% хромистых сталей ферритно-мартенситного класса // Радиационная физика твердого тела // Труды 14 Международного совещания, Севастополь, 2004. / М. НИИ ПМТ МГИЭМ (ТУ). - 2004. - С. 620-624.
  11. ОбразцоваС.М., БиржевойаГ.А., СоловьеваВ.А., ТрояноваВ.М. РачковаВ.И. Нейросетевая оптимизация химического состава стали ЭП-450 по критерию повышения пластичности // Cб. докладов 7аРоссийской конференции по реакторному материаловедению. Димитровград, 8-12.09.03 / Димитровград. - 2004. - Т. 3. - Часть 2. - C. 36-44.
  12. ОбразцоваС.М., БиржевойаГ.А., КонобееваЮ.В., РачковаВ.И., СоловьеваВ.А. Разработка нейросетевой модели кратковременных механических свойств реакторных сталей и интерфейса пользователя. // Вопросы атомной науки и техники. Серия Материаловедение и новые материалы. - 2006. - Вып. 1. - 66. - C. 418-424.
  13. ОбразцоваС.М., БиржевойаГ.А., КонобееваЮ.В., РачковаВ.И., СоловьеваВ.А. Нейросетевая модификация стали ферритно-мартенситного класса ЭП-450 по критерию максимума прочности и пластичности // Перспективные материалы. - 2005. - № 4. - С. 14-19
  14. ОбразцоваС.М., БиржевойаГ.А., КонобееваЮ.В., РачковаВ.И., СоловьеваВ.А. Нейросетевые эксперименты по взаимовлиянию легирующих элементов на механические свойства 12% Cr ферритно-мартенситных сталей // Известия вузов. Ядерная энергетика. - 2008. - №3. - C. 119-124.
  15. ObraztsovаS.М., BirzhevoyаG.A., Konobeev Y.V., RachkovаV.I. // Bootstrap calculation of ultimate strength temperature maxima for neutron irradiated ferritic/martensitic steels // Journal of Nuclear Materials. - 2006. - 359. - P. 263-267.
  16. РачковаВ.И., ОбразцоваС.М., БиржевойаГ.А., КонобееваЮ.В., ПеченкинаВ.А. Нейросетевой анализ влияния легирующих элементов на радиационное охрупчивание материалов корпусов ВВЭР-440 // Атомная энергия. - 2006. - Т. 10. - Вып. 5, ноябрь. - C. 353-358.
  17. ОбразцоваС.М., ПеченкинаВ.А., ПримаковаН.Г., АндриановаА.Н. Нейросетевой прогноз максимума предела текучести сплава Zr+2,5% Nb в зависимости от коэффициента стехиометрии гидрида циркония // Физика и химия обработки материалов. - 2009. - № 2. - С. 33-35.
  18. ОбразцоваС.М., КонобееваЮ.В., ПеченкинаВ.А., СоловьеваВ.А. Нейросетевое исследование влияния кремния и титана на радиационное распухание аустенитной нержавеющей стали // Радиационная физика твердого тела // Труды 20 Международного совещания, Севастополь, 2010. / М. НИИ ПМТ МГИЭМ (ТУ). - 2010. - С. 63-67.
Авторефераты по всем темам  >>  Авторефераты по техническим специальностям