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

Российский научный центр Курчатовский институт

ИНСТИТУТ

РЕАКТОРНЫХ МАТЕРИАЛОВ И ТЕХНОЛОГИЙ

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

УДК 621.039.531

 

Сергеева Людмила Васильевна

ХАРАКТЕРНЫЕ ОСОБЕННОСТИ РАСЧЕТНОГО ОБОСНОВАНИЯ ПРОЧНОСТИ ЭЛЕМЕНТОВ КОНСТРУКЦИЙ  ЯДЕРНЫХ РЕАКТОРОВ  НА СТАДИИ ЭКСПЛУАТАЦИИ И ПРИ СОЗДАНИИ НОВЫХ УСТАНОВОК 

Специальность: 05.14.03- Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации

Автореферат

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

доктора технических наук

Москва-2007

Работа выполнена в Российском научном центре Курчатовский институт

       

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

доктор технических наук, профессор Бараненко Валерий Иванович

доктор технических наук Хмелевский Михаил Яковлевич

доктор технических наук, профессор Щепинов Валерий Павлович

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

ФГУП ОКБМ им. И.И. Африкантова (603074, г. Нижний Новгород, Бурнаковский пр.,15)

               

Защита диссертации состоится_________2007г. в ____ч.____мин.

на заседании диссертационного совета Д 520.009.06 в Российском научном

центре Курчатовский институт по адресу 123182, г. Москва, пл. Курчатова, д.1.

С диссертацией можно ознакомится в библиотеке РН - Курчатовский институт

Автореферат разослан_____________2007г.

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

доктор технических наук, профессор                       В.Г.Мадеев

Актуальность проблемы

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

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

Эти расчеты предполагают рассмотрение всего процесса деформирования кладки вплоть до стадии образования и роста трещин.

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

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

Цели и задачи работы

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

  В числе задач данной работы следует назвать:

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

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

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

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

- исследование устойчивости оболочек твэлов;

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

- проведение расчетов напряженно-деформированного состояния в трубопроводах второго контура АЭС с ВВЭР, имеющих утонения вследствие эрозионно-коррозионного износа;

- разработку математической модели и вычислительной программы нодульной коррозии канальных труб;

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

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

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

- написание программы и расчеты на прочность элементов графитовых отражателей ВТГР.

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

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

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

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

  С участием автора была получена функция повреждаемости для  материалов оболочек твэлов ВВЭР.

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

Результаты исследования прочности графитовых втулок сложного профиля вошли в УНормы расчета на прочность типовых узлов и деталей из реакторного графита уран-графитовых реакторовФ.

Программа расчета напряженно-деформированного состояния графитовых блоков отражателей ВТГР передана в отдел высокотемпературной энергетики и внедрена в практику вычислительных расчетов.

  Методика и программа для расчета площади проходного сечения сквозных трещин, как кольцевых, так и продольных, в стенке корпуса реактора типа ВВЭР были  в составе вероятностной методики и программы переданы в ОКБ Гидропресс и другие организации отрасли, а также в Китайскую народную республику.

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

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

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

- методику расчетного исследования кинетики роста  трещин в элементах конструкций активных зон с учетом воздействия внешней среды;

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

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

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

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

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

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

- проведение расчетов напряженно-деформированного состояния в трубопроводах второго контура АЭС с ВВЭР, имеющих утонения вследствие эрозионно-коррозионного износа;

- расчеты на прочность элементов графитовых отражателей ВТГР.

Апробация работы

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

  Цикл работ, выполненных Сергеевой Л.В. и в соавторстве, был отмечен почетным дипломом Академии наук СССР (1983г.), премиями им. И.В. Курчатова на конкурсе научных работ РН - Курчатовский институт (трижды).

Сообщения и публикации по теме диссертации

По теме диссертации опубликовано более 120 работ, из них более 35 печатных. 

Работы, составляющие основное содержание диссертации, опубликованы:

  • в отечественных журналах: Атомная энергия, Вопросы атомной науки и техники (ВАНТ), Вестник машиностроения, Справочник.  Инженерный журнал, в зарубежном журнале Nuclear Engineering and Design,
  • в трудах международных конференций л Proceedings International Conference on Pipeline Safety, 1997г. и 1999г.

- в материалах отраслевого семинара "Вопросы прочности и надежности элементов активных зон энергетических ядерных реакторов", г. Обнинск, 1982г.

  • в сборнике научных трудов МИФИ Проблемы материаловедения атомной  техники, Москва, Энергоатомиздат,  1989г.

- в материалах семинара "Прочность и  надежность элементов активных зон  энергетических ядерных реакторов",  Обнинск, 1991г.

Структура диссертации

Диссертация изложена на  293 страницах, включая 118 рисунков и 11 таблиц, а также список использованных источников из 238 публикаций; состоит из введения, 5-ти глав и заключения (общих выводов).

Содержание работы

Во введении сформулированы задачи настоящей работы

Глава 1.

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

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

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

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

Уравнения равновесия в полярной системе координат имеют вид:

,

  (1)

где r, z - радиальная и осевая координаты соответственно; - полярный угол,- компоненты тензора скоростей напряжений.

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

; ;;;,  (2)

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

Физические соотношения в данном случае требуют более подробного рассмотрения.

Скорость пластических деформаций в рамках рассматриваемой модели определяется по формуле:

  (3)

где Skl - компоненты девиатора напряжений (индексы i,j,k,l - сочетания из индексов r,, z); - коэффициенты анизотропии; Ф - функция трансляционного упрочнения, зависящая от работы пластического перемещения поверхности текучести после начала нагружения; - эффективное напряжение, контролирующее процесс пластических деформаций;

- частные производные функции текучести ; t-температура; q - доза накопления радиационных повреждений.

Введем обозначения ;;. Тогда ;

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

  (4)

где

, , 

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

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

. (5)

Выражение для матрицы - столбца скоростей пластических деформаций:

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

,где   и

  (6)

,,,,,,.

Остальные компоненты матрицы были получены аналогично (см. диссертацию) и в силу громоздкости выкладок здесь не приводятся.

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

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

- функция трансляционного упрочнения ползучести;- работа вязкого деформирования единичного объема, связанного со смещением поверхности потенциала ползучести в пространстве напряжений;

- координаты центра поверхности потенциала ползучести до начала нагружения;

, - коэффициенты анизотропии.

, где - мощность вязкого деформирования, связанного с формоизменением единичного объема; - эффективное напряжение, контролирующее скорость ползучести;

Вводя обозначения ;;;;; и, используя предшествующие выражения,  можно получить:

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

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

, где - матрица коэффициентов анизотропии при упругом деформировании:

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

,

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

Суммируя деформации упругости, пластичности, ползучести, радиационно-термического формоизменения, распухания и усадки, получим:

,  (7)

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

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

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

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

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

  (8)

где П - потенциальная энергия тела с трещиной; F - поверхность трещины;-вектор приращения трещины.

  После преобразований было получено

  (9)

Здесь - вектор напряжений;  [B] - матрица дифференцирования перемещений; det J - определитель матрицы Якоби; - объем элемента в локальных координатах, изменяющихся в пределах [-1, 1].

Его можно рассчитать, используя метод эквивалентного объемного интегрирования

Для случая трещины нормального отрыва

    (10)

где - поверхность малого цилиндра вокруг участка фронта трещины радиусом ε и высотой Δ; W - плотность работы напряжений на механической части деформаций, , - напряжения и перемещения в системе координат трещины; - компоненты внешней нормали к поверхности.

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

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

    (11) 

При отсутствии неупругих деформаций второе слагаемое в формуле (11) тождественно равно нулю.

 

Задаваясь функцией s  в параметрическом виде и используя интерполяцию перемещений, можно представить первое слагаемое как

где   

Площадь s-функции на

 

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

Глава 2.

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

Геометрия конечных элементов строится таким образом, что поверхности эленментов криволинейны, а поперечные сечения по толщине элементов представляют собой плоскости. Форму таких элементов можно описать с помощью задания декартовых координат точек, лежащих на поверхности.        Обозначим криволинейные координаты в срединной плоскости оболочки  ξ и η, линейную координату по толщине ζ. Предполагая, что  ξ,, η, ζ изменяются в пределах от Ч 1 до 1, зависимость

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

Здесь Ч функция формы, равная единице в i -м узле и нулю в остальных; n Ч число узлов.

 

Р и с. 1. Геометрия конечных элементов. Локальные и глобальные координаты

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

  (13)

где u, v, w Ч перемещения в направлении осей х, у, z; Ч толщина оболочки, соотнветствующая узловой точке i , ,  Ч единичные векторы векторов и .

  где

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

. Здесь , Ч локальные координаты узловых точек.

Для угловых узлов элементов второго порядка   для узлов на сторонах ; ,

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

Предполагалось, что число точек интегрирования в оболочках в обоих направленниях одинаково.

Матрица жесткости имеет вид

[K]dxdydz, где [К]= [Вт] [D] [В].        Здесь [D] Ч матрица упругости.

Матрица [В] является связующей между матрицей перемещений элемента и матрицей деформаций : . В соответствии с теорией  оболочек

{ε'} = =   (14)

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

В формуле (14) деформации в направлении z' не учитываются. Аналогично выглядят выражения для напряжений. Путем преобразований матнрицу можно привести к виду где n Ч число узловых точек. Очевидно, что [В ] = [G ] [S ] [Ф ]. Матрицу [G] получим из следующего преобразования:

Для этого используется обычная операция преобразования глобальных производнных от перемещений u, v, w в локальные производные от локальных ортогональных перемещений:

 

Здесь   Ч матрица ортов по осям х',  у', z', являющаяся по существу матрицей направляющих косинусов. Для ее определения сначала нужно установить направление локальных осей. Вектор, нормальный к поверхности = const, находится как векторное произведение любых двух векторов, касательных к этой поверхности:

Матрица [G] приведена в диссертации. Матрица [S] получается из преобразования

(15)

[S] Ч представляет собой матрицу, на главной диагонали которой стоят элементы обратной матрицы от матрицы Якоби [J].

Коэффициенты матрицы Якоби получаются путем дифференцирования уравнений (12) и имеют вид

в векторе в данном случае 1 Ч номер компонента, 3 Ч номер самого вектора;

Остальные компоненты выглядят аналогично.

Далее вычисляется матрица [Ф], вводимая для получения связи между матрицей-столбцом частных производных от перемещений в декартовых координатах по криволинейным координатам и матрицей-столбцом перемещений узловых точек эленмента:

Здесь un, vn, wnЧ перемещения в направлении координатных осей х, у, z глонбальных координат; Ч углы поворота вокруг ортогональных векторов и ; nЧ номер узла.

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

Для верификации методики были выполнены следующие тесты.

  Тестирование напряжений в полом цилиндре.

Тест на растяжение пластины. При решении задачи даже только с двумя квадратичными изопараметрическими конечными эленментами перемещения верхнего края получаются равными u = 0,313.102 мм, что совнпадает с аналитическим решением.

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

Здесь h Ч толщина пластины; Ч коэффициент Пуассона; а Ч ширина пластины.

Наконец, сравнились результаты расчета зоны патрубков, полученные по разрабонтанной программе и приведенные в работе Бреббия и др. (см. /30/). Они также совпали (рис. 3). Различие в численном и аналитическом решении не превышало 5%.

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

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

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

Наконец, для тестирования программы расчетные данные сравнивались с экспериментальными данными, которые представлены как результаты исследования напряженно-деформированного состояния тройниковых соединений, полученных для  пересечения двух цилиндрических оболочек равного диаметра  (2R  =2r) и имеющих толщину (H= h). Два патрубка нагружены изгибающим моментом M, действующим в плоскости тройника. Момент действует так, что максимальными являются растягивающие напряжения.

Р и с. 2. Пример конечно-элементного разбиения пластинки в тесте

на изгиб (а) и перемещения пластины под действием распределенной нагрузки (б):-------теоретическое решение; ЧЧЧЧЧЧЧконечно-элементный расчет.

Р и с. 3. Сравнение результатов расчета для зоны патрубков, полученных по разработанной конечно-эленментной программе (-------) и в работе Бреббия  и др. (см. /30/) ( ______).

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

 

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

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

Была использована теория течения.

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

 

Рис.5

Деформации зоны патрубков при нагружении  моментами.

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

На АЭС с ВВЭР эрозионно-коррозионному износу подвержены трубопроводы второго контура, работающие в однофазных (конденсато-питательный тракт) и двухфазных (влажно-паровой тракт) потоках. На зарубежных АЭС эта проблема стоит столь же остро.

  В основном процессу эрозионно-коррозионного износа подвержено теплотехническое оборудование (64%), среди которого 25% занимают трубопроводы.

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

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

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

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

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

Рис. 6.

Рис. 7.

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

Таблица 1.

Отношение утоненной части к номинальной толщине

Отношение максимального напряжения к среднему (номинальному)

 

1

2

0,6

2,89

2,32

0,7

2,45

1,94

0,85

1,86

1,62

0,90

1,71

1,57

0,95

1,58

1,56

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

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

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

    (17)

здесь А и n - константы, зависящие от температуры и качества воды,

  (18)

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

  ;   (19)

 

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

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

    (21)

где .

Были выполнены  расчеты для исследования возможного растрескивания опускных  трубопроводов Ду300 АЭС с реакторами РБМК-1000, которые представляют собой трубы диаметром 300 мм, толщиной 15 мм, расположенные вертикально между барабан-сепараторами и всасывающим коллектором.

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

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

Таблица 2.

Глубина трещино-образно-  го дефекта

(мм)

Тол-щина
в зоне

пов-

реж-
дения

(мм)

Величина  остаточ-ных напря-жений

Значения коэффициента

  Интенсивности  напряжений 

Скорость

роста

трещины

мм/с

Время прохож-

дения
трещины насквозь  (лет)

0,45

17,25

  30.

  65.6

  67.04

44.76

4.88

  0.97

0,45

17,25

  40.

78.93

  80.65

53.84

7.07

  0.67

0,45

17,25

  50.

92.25

  94.25

62.92

9.65

  0.49

0.45

17.25

  20.

52.30

  53.44

35.67

3.10

  1.53

0.45

17,25

  0.

25.67

  26.22

17.51

7.47

  6.36

0.225

17.24

  0.

18.15

  18.26

12.38

3.62

13.1

0.675

17.25

  0.

31.43

  32.62

21.44

1.16

4.1

0.45

16.5

  30.

66.17

  67.61

45.14

4.97

1.038

0.3

16.5

  40. 

64.9

  65.64

44.27

4.68

0.984

0.3

16.5

  30.

54.03

54.64

36.8

5.4

0.682

0.3

16.5

  20.

43.16

43.64

29.44

2.07

0.435

0.45

16.5

  20.

52.86

  54.0

36.05

5.19

0.666

В представленной таблице , ,- коэффициенты интенсивности напряжений, полученные по формулам (19).

Глава 3.

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

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

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

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

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

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

В пределах шага по времени приращение полной деформации представляем в виде суммы упругой составляющей, составляющей ползучести и изменения объема:

(22)

  Приращение упругой составляющей для анизотропного материала:

(23)

и вычисляем по этой же формуле с круговой перестановкой индексов,

(24)

- значение модуля упругости в i-ом направлении, - значение коэффициента Пуассона j-ом направлении при нагружении материала в направлении i.

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

, (25)

где приращение объемной составляющей ползучести

 

Здесь- скорость изменения функции объемной ползучести

Приращение девиатора деформаций ползучести

(26)

Приращение составляющей изменения объема записываем следующим образом:

(27)

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

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

ζ

  Рис.8

                                               

        η

  ξ

 

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

Для узлов в вершинах элемента

Для узлов на серединах ребер

Матрица жесткости элемента   (28)

Выполняя перемножение по формуле (28), получим выражения для коэффициентов матрицы жесткости трехмерного квадратичного элемента:

, ,

, ,

,,

,

Здесь , , ;  m, n - верхние индексы, которые соответствуют локальным номерам узлов. Матрица деформаций может быть представлена в виде:

Здесь i,j,Еp - номера узлов.

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

  Матрица Якоби имеет вид: 

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

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

2) Были выполнены тесты  напряженно-деформированного состояния графитового  блока, обусловленные внутренним и наружным давлением на соответствующие поверхности.

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

Результаты расчета приведены на рис 9 и 10.

Рис. 9.

Рис.10.

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

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

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

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

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

В качестве примера можно привести результаты расчета напряженно-деформированного состояния графитового блока и трубы ТК реактора РБМК-1000.

а б

Рис.11. Распределение интенсивности напряжений по области определения, соответствующее эксплуатации в течение 12,7 (а) и 19,7 г. (б).

 

Рис. 12. Распределение интенсивности деформаций по области определения, соответствующее 19,7 годам эксплуатации.

Распределение интенсивности напряжений, соответствующее 12,7 и 19,7 годам работы конструкции, приведено на рис. 11. Распределение интенсивности деформаций по области определения, соответствующее 19,7 годам эксплуатации - на рис.12. Распределение радиальных составляющих перемещений по области определения, соответствующее  12,7 и 19,7 годам работы конструкции  - на рис. 13.

13. Распределение радиальных составляющих перемещений по области определения, соответствующее эксплуатации в течение 12,7 (а) и 19,7 г. (б).

Глава 4.

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

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

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

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

Было проведено исследование влияния учета фактора анизотропии на процесс деформирования оболочки твэла. Учет анизотропии проводился посредством метода Хилла. Расчет формоизменения топливных таблеток проводился посредством двумерной конечно-элементной методики расчета на прочность осесимметричных тел (в r-z координатах).

  Все расчетные исследования проводились в предположении, что твэлы эксплуатировались в стационарных условиях. Величина исходной овальности варьировалась от нуля до максимального значения u0, которое определялось исходя из допуска на внутренний диаметр оболочки, в частности для РБМК - Dвнутр=11,7+0,1мм и соответственно u0=0,05 мм, а для ВВЭР Dвнутр=7,85+0,05мм и  u0=0,025 мм.

  На рис. 14 приведены зависимости изменения радиуса анизотропной оболочки твэла РБМК, соответствующего двум взаимно перпендикулярным осям симметрии (1 и 2) для максимально допустимой овальности и наружного радиуса топлива. Аналогичный расчет проведен для изотропной оболочки для тех же условий эксплуатации. Этот и аналогичные расчеты показали, что неучет анизотропии приводит к существенной разнице в результатах. Из расчета с учетом анизотропии получаем, что, в рамках принятых допущений, зазор между оболочкой и топливом в твэле РБМК при максимальной исходной овальности выбирается приблизительно за 430 эффективных суток работы на постоянной мощности, в то время как из расчета без учета анизотропии получаем, что зазор перекроется приблизительно через 60 эффективных суток эксплуатации.

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

взаимно перпендикулярных осей симметрии (1 и 2) и наружного радиуса топливной таблетки от времени эксплуатации t, эффективных суток (начальная овальность uo=0,05 мм).

Таким образом, для изотропных оболочек увеличение овальности происходит существенно быстрее, время до посадки оболочки на топливо в несколько раз меньше, чем у анизотропной. Аналогичные расчеты проведены для твэлов ВВЭР. В аварийных условиях в проблеме устойчивости оболочек твэлов возникают новые аспекты. Во-первых, поведение оболочек твэлов реакторов ВВЭР и РБМК различно. В частности, при разрыве контура циркуляции реактора давление теплоносителя за 20-30 секунд падает практически до атмосферного и деформирование оболочки происходит наружу, поэтому в этом случае говорить об устойчивости оболочек нет смысла. В реакторах РБМК, напротив, при разрыве напорного коллектора, давление теплоносителя падает сравнительно медленно,  поэтому в этом случае деформирование оболочек твэлов может происходить внутрь. В этом случае, как показали расчеты, посадка оболочки на топливо для свежего твэла произойдет практически мгновенно, для твэла, который долго находился в эксплуатации, это определяется давлением ГЦД. Во- вторых, в случаях посадки оболочки на топливо существенное влияние на процесс деформирования оболочки могут оказать осевые зазоры между топливными таблетками. При достижении оболочкой температуры, превышающей 600оС, механические характеристики сплава Н-1 ухудшаются, оболочка уже не может противостоять давлению теплоносителя и начинается интенсивный процесс деформирования в осевые зазоры между таблетками.

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

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

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

  Рис.15.

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

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

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

Здесь N- концентрация зародышевых субмикротрещин:

где N0 - концентрация зародышевых трещин, соответствующая насыщению процесса их накопления; c,b - коэффициенты, зависящие от температуры и уровня приложенных напряжений.

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

 

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

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

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

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

Рис.16  Смена механизмов разрушения в процессе роста  трещины

H - толщина оболочки, I - инкубационный период, IIЦэтап разрушения, обусловленный достижением функцией повреждений предельного значения, III - этап коррозионного растрескивания, IV - этап водородного охрупчивания.

 

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

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

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

Металлографические исследования микроструктуры внутренней поверхности труб показали, что кроме нодулей 1-2 мм встречаются отдельные коррозионные очаги размером 20-80 мкм и меньше. За счет коррозионного объединения мелких мононодулей могут возникать полинодули размером 1000-2000 мкм.

Тщательное микрослойное изучение поперечных сечений (профилей) отдельных нодулей или нодулей, находящихся в составе коррозионных очагов, позволило выявить их существенные особенности:

  -  мононодуль не является двояковыпуклой линзой,

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

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

Возникновение напряжений можно объяснить радиационным распуханием карбида циркония, которое составляет 3% при флюенсе нейтронов 1021 см-2 при температуре облучения 300-700оС. Негативная роль интерметаллидов циркония объясняется их невысокой коррозионной стойкостью и последующим распуханием продуктов окисления.

Для проверки сказанного была разработана математическая модель и выполнен расчет возможного растрескивания при увеличении на 3 % объема включений второй фазы, залегающих вблизи поверхности нераспухающей матрицы, при  флюенсе  1021 см-2. Охрупчивание и растрескивание происходит при выходе на поверхность пластически деформированной зоны, в которой интенсифицируются процессы гидрирования. При этом диаметр пластической зоны при выходе на поверхность (dпл) в K раз превышает диаметр включения (dвкл). В результате металлографических исследований поверхности канальных труб, пораженных нодульной коррозией, получено, что K= dпл/ dвкл приблизительно равен 1,5-2.

Рис.17. Вид вскрытого нодуля.

Вскрытый нодуль, содержит крупное выделение второй фазы (типа представленного на рис.18).  Частицу размером приблизительно 15 мкм окружает слой гидридных пластин белого цвета (их поперечный разрез показан на рис. 18). Такие выделения идентифицированы как карбид циркония (микротвердость в исходном сплаве 2,7.104 МПа, в облученном до флюенса нейтронов 1021 см-2 3,3.104 МПа).

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

Рис. 18

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

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

  Вариантными расчетами  выявлены соотношения между размерами включения и расстоянием от его центра до поверхности трубы, соответствующие условию K= dпл/ dвкл =1,5 (см. рис. 20,б). Заштрихованная зона соответствует области диаметра включений и расстояния от него до поверхности трубы, при которых происходит разрушение перемычки между включением и поверхностью. Из расчетного графика видно, какая глубина залегания для включения заданного диаметра необходима для того, чтобы нодульная коррозия не проявилась.

 

Рис 19 Растрескивание циркониевого сплава у внешней поверхности канальной трубы вблизи отдельной карбидной частицы (х1000)

На рис.20,а показаны примеры расчета границ зон пластичности для включений размером 0,005; 0,01; 0,02 мм, находящихся на расстоянии от поверхности трубы, 0,01; 0,0125; 0,015; 0,03 мм для моментов времени 105 с: 2,763; 2,752; 2,492 и 3,312 соответственно.

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

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

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

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

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

Докритическое подрастание трещины  в результате водородного охрупчивания материала было задано в виде зависимости:

где - коэффициент интенсивности напряжений; - температура; oC.

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

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

где - приращение потенциальной энергии тела, вызванное варьированием глубины трещины; - коэффициент Пуассона; Е - модуль упругости первого рода.

 

Рис. 21. Глубина трещины в зависимости от времени, распределение интенсивности остаточных напряжений по толщине трубы.

Потенциальная энергия может быть вычислена с помощью соотношения:

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

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

,

где L - длина трещины; - номинальное окружное напряжение, - предел текучести;

- предел прочности; Rср - радиус срединной поверхности трубы; Н - толщина канальной трубы.

Формула справедлива для случая плоско-деформированного состояния и для .

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

где - толщина канальной трубы.

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

Плотность потока нейтронов принималась равной 2.1013см-2с-1. Эксплуатация трубы технологического канала рассматривалась при стационарном максимальном уровне мощности. Выход на мощность производился в течение суток. Предполагалось, что кампания составляет 15 тысяч суток.

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

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

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

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

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

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

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

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

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

    • механических свойств материала корпуса:

σт - предел текучести, МПа, σv - предел прочности, МПа, δ - равномерное удлинение, %, Е - модуль упругости, МПа, μ - коэффициент Пуассона;

    • условий нагружения, длины трещины l, мм.

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

  Так как рассматривались продольные и кольцевые трещины, размеры которых (длина 50-500 мм)  невелики по сравнению с размерами корпуса (внутренний диаметр корпуса 2075 мм), в качестве первого приближения задача решалась в плоской постановке (x-y геометрия).

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

где l - полудлина сквозного дефекта, R- - средний радиус стенки корпуса, S - толщина стенки корпуса.

 

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

  1. Сопоставление результатов расчета, полученных по двумерной методике, с результатами решения трехмерной задачи.
  2. Проверка точности расчета напряженно-деформированного состояния.
  3. Сравнение раскрытия трещины, полученного экспериментальным  и расчетным путем.

Трехмерная методика дает большие раскрытия трещины: в данном случае площади относятся как 11,5 мм2 к 13,9 мм2. 

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

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

Решение задачи рассматривалось на основе теории пластичности Хубера-Мизеса

Рассматривался диск, нагруженный внутренним давлением, равным .

Из сравнения видно, что полученные результаты и известное решение совпадают достаточно точно (рис.22).

Рис.22.

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

Эксперимент заключался в том, что образец, в центре которого было просверлено отверстие, в обе стороны от которого были сделаны два надреза, подвергался циклическому нагружению с частотой 200 циклов/мин на пульсаре ГРМ-1 до получения на концах надреза усталостных трещин длиной l=2 мм каждая. Максимальное напряжение цикла при получении трещин обычно составляло 15 кг/мм2 .

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

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

 

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

Поправка на кривизну, вводимая для продольной трещины, в случае длины трещины 50 мм, составляет лишь 1%.

Таблица 3.

Вид

Трещины

Длина

трещины, мм

Площадь

трещины, мм2

Максимальное раскрытие трещины, мм

Продольная

Трещина

50

2,79

0,0776

100

11,5

0,1598

200

51,0

0,356

300

134,0

0,616

400

286,0

0,986

500

542,0

1,498

Кольцевая

Трещина

50

1,32 (1,33)

0,0364 (0,0366)

100

5,28 (5,48)

0,0728 (0,0856)

200

21,1 (24,4)

0,1458 (0,1681)

300

47,5 (64,0)

0,218 (0,394)

400

84,5 (137,0)

0,292 (0,472)

500

132,0 (259,0)

0,364 (0,716)

 

Площади, ограничиваемые контурами трещины, представлены также в таблице 3.

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

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

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

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

Алгоритм метода можно коротко записать в виде:

, +, +,  (30)

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

Итерационный процесс продолжается до тех пор, пока не будет выполнен выбранный критерий сходимости:

(31)

Где заданная точность вычислений.

Применим обозначения; ;  ; ;;  , где ,- девиатор напряжения  и деформации; G - модуль упругости при сдвиге;  - интенсивность напряжения, соответствующая точке на диаграмме упрочнения материала (используется критерий текучести Мизеса).

В пластической области ;.

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

(32)

Это выражение позволяет определять напряжения по теории течения.

Результаты расчета показали, что переход из упругого состояния в упруго-пластическое происходит при определенном значении отношения граничного напряжения, приложенного перпендикулярно поверхности трещины, к пределу текучести. Для продольной трещины это отношение равно  / =0,4, а для кольцевой трещины -  / =0,33.

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

  При таких напряжениях рассматривать данную задачу не целесообразно, поэтому за верхнюю границу нагружения брали данное напряжение. Однако при малых пределах текучести  (392МПа) это отношение может достигать большего значения, равного 0,61.

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

1 - от = 0,33 до = 152,88 МПа, 2 - от  = 152,88 МПа до = 180,32МПа.

Результаты расчетов для кольцевых трещин представлены в таблице 4, здесь  δ=15%;  μ=0,3;  Е=1,94.104 кг/мм2 ,  l= 200мм.

Таблица 4

σт, кг/мм2

σV, кг/мм2

S, мм2

σz, кг/мм2

38

43,5

82,96

13

38

43,5

160,85

14,7

38

43,5

208,75

15,6

45

50,6

102,24

16,4

50

55,6

55,734

16,6

50

55,6

72,145

17,7

42

47,6

83,753

14,5

42

47,6

117,69

15,6

42

47,6

179,753

16,6

42

47,6

249,14

17,7

45

52

64,614

15,6

45

52

90,49

16,6

45

52

129,68

17,7

45

55

46,611

15,6

45

55

76,493

17

45

57

47,301

15,6

45

57

66,603

17

45

57

89,944

18,4

45

59

43,579

15,6

45

59

60,321

17

45

59

76,863

18,4

Аналогичные результаты, полученные для продольных трещин, приведены в диссертации.

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

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

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

В пределах упругости

За пределами упругости

Для продольной трещины

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

  для кольцевой трещины

для первого диапазона

в частном случае

для второго диапазона

в частном случае

где S - площадь трещины, мм2,- напряжение, нормальное к плоскости трещины, МПа, l - длина трещины, мм, - поправка на кривизну, p - давление, МПа, - равномерное удлинение,%, - окружное напряжение, МПа, - осевое напряжение, МПа.

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

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

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

В таблице 5 представлена небольшая часть (пример) результатов для критического раскрытия трещины из итогового отчета по указанной теме

                                                                                       Таблица        5                                                                        

L

N

M

σp

σq

δ

δNN

δNC

δa

1

2

3

4

5

6

7

8

9

10

5,4

мм

1

σ02h

0

27

0

1,596.10-2

1,596.10-2

3,192.10-2

1,596.10-2

2

σ02h/4

σ02h2/12

6,75

13,5

6,214.10-3

8,268.10-3

1,654.10-2

1,053.10-2

3

σ02h/4

σ02h2/6

6,75

27

8,538.10-3

1,469.10-2

2,939.10-2

2,041.10-2

4

σ02h/4

σ02h2/4

6,75

40,5

1,086.10-2

2,112.10-2

4,224.10-2

2,926.10-2

5

0

σ02h2/12

0

13,5

2,324.10-3

3,563.10-3

2,138.10-3

2,138.10-3

5a

qi=0,045 кг/мм2

2,488.10-3

3,563.10-3

7,126.10-3

3,563.10-3

6

0

σ02h/6

0

27

4,648.10-3

9,989.10-3

1,998.10-2

9,989.10-3

6a

qi=0,09 кг/мм2

4,975.10-3

9,989.10-3

1,998.10-2

9,989

10-3

7

0

σ02h2/6

0

40,5

6,972.10-3

1,642.10-2

3,283.10-2

1,642.10-2

7a

qi=0,135 кг/мм2

7,461.10-3

1,642.10-2

3,283.10-2

1,642.10-2

8

σ02h

σ02h2/6

27

27

2,020.10-2

2,881.10-2

5,762.10-2

5,120.10-2

7,2

мм

1

σ02h

0

27

0

2,488.10-2

2,769.10-2

5,537.10-2

2,769.10-2

2

σ02h/4

σ02h2/12

6,75

13,5

9,408.10-3

1,350.10-2

2,699.10-2

1,923.10-2

N- растягивающее усилие, кг/мм2, M - изгибающий момент, кг, - распределенная растягивающая нагрузка на единицу площади, кг/мм2 - распределенная изгибающая нагрузка на единицу площади, кг/мм2 qi - распределенная изгибающая нагрузка, прикладываемая по длине стержня, кг/мм2  δ - расчетное раскрытие трещины, полученное с помощью предлагаемой методики, мм,

δNN - расчетное раскрытие трещины, вычисленное по методике НИКИЭТ, с учетом поправки. δNC  - расчетное раскрытие трещины, вычисленное по методике НИКИЭТ. δa- расчетное раскрытие трещины с учетом поправочных коэффициентов; ллав означает лавинообразный рост трещины.

Коэффициенты PP, PQ, и PPQ, значения которых приведены в диссертации и отчете по теме.

  Считается, что водород мигрирует в зоны наибольших напряжений, в частности, есть мнение, что он скапливается в пластических зонах, причем скорость гидрирования связана с уровнем напряжений. Было выполнено исследование влияния поверхностных дефектов на напряженное состояние и возможность растрескивания оболочек твэлов ВВЭР и РБМК. Проведены расчеты концентрации напряжений в зоне дефектов типа царапин, имеющих место на внутренней поверхности оболочек. Результаты для оболочек твэлов ВВЭР приведены в таблице 6 (внутреннее давление газообразных продуктов деления: pвнут = 2,6 МПа) и таблице 7 (изнутри на оболочку давит растрескавшееся топливо p = 320 атм), внешнее давление теплоносителя в обоих случаях одинаковое.

  Таблица 6

l/мкм/

R/l

40

50

100

500

1/2

81,3

84,9

104

387

1/5

81,5

84,9

106

460

1/10

117

123

127

596

1/200

184

190

225

Таблица 7

l/мкм/

R/l

40

50

100

500

1/5

140

146

183

Более 400

1/10

221

234

247

Более 440

1/200

309

322

391

В таблицах полужирным курсивом (оранжевый цвет) обозначены сочетания глубины и остроты R/l дефекта (R - радиус закругления кончика трещины), при котором интенсивность напряжений больше или равна пределу текучести сплава 110, контуром  показаны сочетания размеров дефекта, при которых интенсивность напряжений находится в пределах от 2/3 до предела текучести (синий цвет).

Пределы текучести в этих случаях различные, в необлученном состоянии при температуре, соответствующей температуре эксплуатации, предел текучести равен 100 МПа, а при флюенсе Q=4,7-9,2-1020нейтр/см2 и той же температуре предел текучести равен приблизительно 300 МПа.

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

  Для оболочек твэлов РБМК-1000, работающих в пароводяной среде, целесообразно ограничить допустимую глубину дефектов в виде царапин величиной 50 мкм, так как в противном случае для наиболее вероятных значений показателей остроты трещин (R/l1/10) величина напряжений в районе рисок будет превышать уровень, близкий к , а указанный уровень соответствует началу ускоренного гидрирования циркониевой оболочки.

Для оболочек твэлов реактора ВВЭР-1000 видимо целесообразно ввести более жесткий контроль по глубине дефектов, ограничив её 40 мкм, так как даже в случае тупых  дефектов интенсивность напряжений в начальный период эксплуатации превышает 2/3 текучести.

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

Глава 5.

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

ВТГР - источники тепла с уникально высокой температурой - около 1000oC, поэтому их использование позволяет значительно расширить сферу экономически эффективного применения ядерной энергии.

Начиная с 70-х годов прошлого века, в стране были выполнены  проекты высокотемпературных гелиевых реакторов ВТГР атомных энерготехнологических станций АЭТС для химической промышленности и черной металлургии, среди которых АБТУ-50, а позднее _ проект атомной энерготехнологической станции с реактором ВГ-400 мощностью 1060 МВт для ядерно-химического комплекса по производству водорода и смесей на его основе, по выпуску аммиака и метанола, а также ряд последующих проектов этого направления.

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

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

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

  Конструкции блоков графитового отражателя, представленные на рис. 23-25, условно назовем конструкциями первого, второго и третьего типов. Характер температурного и радиационного нагружения этих конструкций приблизительно одинаков и может быть качественно охарактеризован графиками, которые представляют собой один из вариантов нагружения для конструкций первого типа (рис. 26, 27).

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

Рис. 23. Первый тип конструкций блока бокового графитового отражателя

Рис. 24. Второй тип конструкции блока  графитового отражателя

Рис. 25.  Третий тип конструкции блока  графитового отражателя

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

Рис.26. Характер распределения температур в Рис.27. Характер распределения плотностей потока

блоке  графитового отражателя  нейтронов (Е>0,1 МэВ) по блоку

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

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

При расчете симметричных вариантов конструкций блоков графитового отражателя рассматривается симметричная часть графитового блока.

Рис. 28. Зависимость глубин трещин от времени для первого варианта конструкции

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

Рис.29. Распределение интенсивности напряжений по сечению конструкции первого типа графитового блока, соответствующее 22,2 года эксплуатации

Рис.30. Распределение интенсивности напряжений по сечению конструкции  графитового блока второго типа, соответствующее 28,6 года эксплуатации

Основные выводы по диссертации

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

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

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

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

3.1. Для определения J-интеграла был использован новый эффективный метод эквивалентного объемного интегрирования.

3.2.  Была разработана  и протестирована программа, которая в качестве критерия использует предельное раскрытие трещины в вершине. С помощью этой программы получен широкий спектр данных для анализа напряженно-деформированного состояния трубопровода Ду-850 реактора ВВЭР-1000 при наличии в его стенке трещиноподобного дефекта.

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

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

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

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

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

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

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

6. Было исследовано влияние такого технологического дефекта как овальность на процесс потери устойчивости оболочек твэлов.

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

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

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

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

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

11.2. Программа в составе вероятностной модели внедрена в ОКБ Гидропресс 

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

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

14. Проведен сравнительный анализ прочности трех типов блоков  графитового отражателя проектируемого реактора ВГ-400.

14.1. Программа передана в отдел высокотемпературной энергетики.

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

Список опубликованных работ по материалам диссертации

1. Исследование растрескивания  графитовых втулок сложного профиля, Сергеева Л.В. статья в сб. "Вопросы атомной науки и техники",  серия "Атомное материалонведение", вып.2 (8), Москва, 1980 г., с.64-82.

2. Расчетное исследование кинетики роста  трещин в элементах конструкций активных зон с учетом воздействия внешней среды, Сергеева Л.В. в сб. "Вопросы атомной науки и техники",  серия "Атомное материалонведение", вып.3 (11), Москва, 1981 г., с.22-29

3. Исследование процесса растрескивания трубопроводов (оболочек) из конструкционных материалов, имеющих структурную неоднородность, Сергеева Л.В., Тутнов А.А., Тутнов И.А., там же, вып. 3 (11), Москва, 1981г., с.3-22

4. Методика расчетного исследования растрескивания  элементов конструкций активной зоны ядерных реакторов, Сергеева Л.В. "Тезисы докладов 1-го отраслевого семинара по прочности  и надежности элементов активных  зон ядерных реакнторов", сб. 1981г., с.20.

5. Методика расчетного исследования растрескивания элементов конструкций активной зоны ядерных реакторов, Сергеева Л.В."Вопросы прочности и надежности элементов активных зон энергетических ядерных реакторов",  материалы отраслевого семинара, г. Обнинск, 1982г.,  стр.59-70

6. Расчет полей напряжений в графитовых блоках отражателей ВТГР с учетом изменения теплофизических и прочностных характеристик, Сергеева Л.В., Костюк Н.Н. статья ВАНТ сер "Физ. и техн. ядерн реак." вып.6, 1985г. с.96-101

7. Методика расчета кинетики изменения напряженно-деформированного состояния оболочек  твэлов энергетических реакторов, учитывающая процесс посадки оболочек на топливный сердечник, Сергеева Л.В., Ткачев В.В., доклад в сб. "Прочность и  надежность элементов активных        зон энергетических ядерных реакторов", Москва, 1986г, с.36-37, ГК ИАЭ ЦНИИАИ

8. Конечно-элементная методика расчета напряженно-деформированного состояния оболочечных элементов конструкций ядерных реакторов сложной пространственной геометрии в условиях вязко- пластического деформирования, Сергеева Л.В., Ткачев В.В., тезисы доклада в сб. "Прочность и надежность элементов активных зон энергетических реакторов", Москва,1986г., с.39

9. К вопросу об устойчивости оболочек твэлов  энергетических реакторов, Сергеева Л.В.,  Рубцов  В.С.,  статья ВАНТ, сер "Физ. и техника ядерных реакторов",  вып.6, 1985г, с.11-16.

10. Технологический канал ядерного реактора, Маневский В.Н, Сергеева Л.В., Тутнов А.А., Тутнов И.А.,  авторское свидетельство № 1103727

11. Расчетное исследование поведения оболочек твэлов РБМК-1000 с учетом начальной овальности в стационарных условиях и в условиях аварийной ситуации, Рубцов В.С, Сергеева Л.В., статья в ВАНТ, сер. "Атомное материаловедение", вып.1(17), 1983г., с.3-9.

12. Методика расчетного исследования кинетики развития субмикротрещин в образцах из материала оболочек твэлов реакторов ВВЭР, Сергеева Л.В., Тутнов А.А,  статья в ВАНТ сер. "Атомное материаловедение", вып.1(23), 1986г., с. 33-39

13. Моделирование кинетики роста субмикротрещин в тонкостенных образцах, Сергеева Л.В., Тезисы доклада  в сб. "Прочность и надежность элементов активных зон  энернгетических ядерных реакторов",  ГК ИАЭ, ЦНИИАИ, 1986г, стр. 105-106.

14. Методика расчета кинетики роста трещины в технологического канала реактора РБМК, Сергеева Л.В., там же, с.106

15.  Методика расчета кинетики роста трещин в трубопроводах под действием водородного охрупчивания,- Сергеева Л.В., ВАНТ, сер. "Атомное материаловедение", 1(26), 1988г., с. 29-31,  Москва, 1985г.

16. Программа TUBE-1 для расчетного моделирования напряженно-деформированного состояния и кинетики роста трещины в канальных трубах РБМК-1000,  Сергеева Л.В.,  препринт ИАЭ-4288/4         напряженно-деформированного состояния и кинетики роста трещины в канальных трубах РБМК-1000, препринт ИАЭ-4288/4

17. Математическое моделирование процесса совместного деформирования призматического  графитового блока и трубы технологического        канала с учетом анизотропии механических характеристик, Сергеева Л.В., Тутнов А.А., ВАНТ, сер. "Атомное материаловедение", 1988, вып.1(26), с. 23-29.

18. Программа ELL. Аннотация, Сергеева Л.В., ВАНТ, сер. "Физика и техника ядерных реакторов",1988г., вып.2, с.72-73.

19. Программа TUBE-1. Аннотация, Сергеева Л.В., ВАНТ, сер. "Физика и техника  ядерных реакторов", 1987г., вып. 8, с. 64-65.

20. Определяющие уравнения для трехмерного расчета  терморадиационных вязко-упругопластических напряжений и деформации в анизотропных телах с  использованием математической модели материала с одновременно смещающимися и расширяющимися поверхностями и текучести и потенциала полнзучести, Сергеева Л.В., Тутнов А.А. ВАНТ, сер. "Атомное материаловедение",  1988г., вып.1(26) с.13-23

21. Программа VG4001. Аннотация, Сергеева Л.В.,  статья ВАНТ серия, "Физика и  техника ядерных ректоров", 1988г., вып.3

22. Кинетика роста субмикротрещин в металлических фольгах, Доровской В.М, Елесин Л.А.,  Сергеева Л.В., Тутнов А.А., Проблемы материаловедения атомной  техники, сборник научных трудов МИФИ под ред. Калинина Б.А. Москва, Энергоатомиздат,  1989г. с. 43-61

23. Методика расчета напряженно-деформированного состояния графитовых блоков установки ВГ-400 с учетом их возможного растрескивания,  Сергеева Л.В., ВАНТ, сер. Материаловедение и новые материалы", вып.1(35), 1990,с. 34-40

24. Расчетное исследование термомеханического взаимодействия труб технологических каналов и элементов графитовой кладки канального реактора, Сергеева Л.В.,  Тутнов А.А., "Атомная энергия", т.68,вып.4, апрель 1990г., с. 236-241

25. Математическое моделирование процесса нодульной коррозии циркониевых труб, Сергеева Л.В., Тутнов А.А.  Материалы семинара "Прочность и  надежность элементов активных зон  энергетических ядерных реакторов", Обнинск, 1991г., с.47-48

26. Исследование прочности модифицированных графитовых блоков для строящихся РБМК, Сергеева Л.В., там же, с.57

27. Математическое моделирование роста трещин в трубопроводах в условиях коррозионно-активной среды, Сергеева Л.В., Тутнов А.А.,  статья ВАНТ, сер. "Материаловедение  и новые материалы", вып.1(35),с.49-55

28. Развитие нодульной коррозии циркониевых труб РБМК-1000 при нейтронном облучении. Карасев В.С., Ковыршин В.Г., Сергеева Л.В., Тутнов А.А.  статья ВАНТ, Физика радиационных повреждений и  радиационных материалов, вып.1(4)/2(5),1989г., с.29-37.

29. Нодульная коррозия циркониевых труб канальных реакторов, Карасев В.С.,  Ковыршин В.Г.,  Колесов В.В., Сергеева Л.В.,  Тутнов А.А., Чирко Л.И.,  Шинаков А.А., "Атомная энергия",  т.72, вып.2,  февр. 1992г., с. 124-130. 

30. Трехмерная программа расчета напряженно-деформированного состояния оболочечных конструкций сложной пространственной геометрии, Сергеева Л.В.,  Атомная энергия, т.80, вып.2, февр. 1996г., с.81-87

31. The calculated determination of stress-strain in the T-pipes at the variation of the backings, Сергеева Л.В.,  Тутнов И.А.  Proceedings International Conference on Pipeline Safety        Conference on Pipeline Safety

32. An investigation of the strength of branch areas in the pipe of nuclear power installations,  Сергеева Л.В., Nuclear Engineering and Design 196, 105-110, 2000 

33. Исследование прочности трубопроводных систем в условиях пластической деформации,  Сергеева Л.В., Тутнов И.А.,  Доклады участников третьей международной конференции "Безопасность трубопроводов" т. 2, Москва, 1999г., с. 42-50.

34. Особенности расчета напряженно-деформированного состояния трубопроводов второго контура ВВЭР, имеющих эрозионно-коррозионные утонения Сергеева Л.В., Киселев А.С.,  Атомная энергия, т. 80, 2002г

35. Исследование процесса возможного разрушения трубопроводов ядерных энергетических установок вследствие их эрозионно-коррозионного утонения, Сергеева Л.В.  журнал УВестник машиностроенияФ,  № 6, 2007г, с. 22-24.

36. Исследование возможного растрескивания корпуса реактора в рамках концепции лтечь перед разрушением,  Сергеева Л.В. УСправочник. Инженерный журнаФ № 7, 2007 г, стр.52-58

37. Методика расчетного исследования напряженно- деформированного состояния графитовых блоков  в свете обоснования продления эксплуатации графитовых кладок РБМК, Сергеева Л.В. статья в печати журнал УАтомная энергияФ.

38. Исследование напряженно-деформированного  состояния элементов графитовой кладки ядерных реакторов с учетом анизотропии графита, Сергеева Л.В.  Журнал УВестник машиностроенияФ, № 7, 2007г, с. 19-21.

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