Научная статья на тему 'Математическое моделирование развития запроектной аварии внутри корпуса реактора на быстрых нейтронах'

Математическое моделирование развития запроектной аварии внутри корпуса реактора на быстрых нейтронах Текст научной статьи по специальности «Строительство и архитектура»

CC BY
387
81
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
РЕАКТОР / ЗАПРОЕКТНАЯ АВАРИЯ / ГИДРОДИНАМИЧЕСКОЕ ДАВЛЕНИЕ / МЕТОД КОНЕЧНЫХ ЭЛЕМЕНТОВ / ПРОЧНОСТЬ / REACTOR / BEYOND THE DESIGN BASIS ACCIDENT / HYDRAULIC PRESSURE / FINITE ELEMENT METHOD / STRENGTH

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Баженов В.Г., Жестков М.Н., Замятин В.А., Кибец А.И.

В трехмерной постановке рассматривается процесс нестационарного деформирования конструкции реактора на быстрых нейтронах с жидким теплоносителем в условиях постулируемой запроектной аварии типа ULOF (Unprotected Loss of Flow). Под этой аварией понимается расплавление активной зоны, вызванное отключением главных циркуляционных насосов первого контура с сопутствующим несрабатыванием аварийной защиты. В результате расплавления активной зоны реактора образуется область энерговыделения с повышенным уровнем давления, заполненная парами натрия. Прогрессирующее расширение области энерговыделения в теплоносителе приводит к росту уровня напряженно-деформированного состояния корпуса реактора и может привести к его разрушению. В этих условиях реакторная установка должна сохранять герметичность, обеспечивать локализацию последствий запроектной аварии внутри силового корпуса и не допускать опасного радиационного воздействия на персонал атомной станции и окружающую среду. Для описания движения теплоносителя и конструктивных элементов реактора применяется текущая лагранжевая формулировка. Уравнение движения выводится из баланса виртуальных мощностей. В качестве физических соотношений для металлов применяются уравнения теории пластического течения. В теплоносителе девиаторные компоненты напряжений полагаются равными нулю, а связь между гидростатическим давлением и плотностью берется в виде уравнения состояния квазиакустического типа. Контакт теплоносителя с конструктивными элементами реактора моделируется условиями непроникания. Решение задачи основывается на моментной схеме метода конечных элементов и явной конечно-разностной схеме интегрирования по времени типа «крест», реализованных в рамках вычислительной системы «Динамика-3». Численно исследуется изменение напряженно-деформированного состояния корпуса реактора на быстрых нейтронах в условиях запроектной аварии типа ULOF. Анализируется возможность локализации последствий запроектной аварии внутри силового корпуса реактора.

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

Похожие темы научных работ по строительству и архитектуре , автор научной работы — Баженов В.Г., Жестков М.Н., Замятин В.А., Кибец А.И.

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

Mathematical modeling of development of beyond design basis accident in integral fast-neuron reactors vessel

The process of nonstationary deformation of construction of the fast reactor with liquid metal coolant under postulated ULOF beyond design basis accident is examined. This type of accident includes core melting caused by disconnection of the main circulation pump arrangements of the primary circuit with the associated failure of the emergency protection. As a result of core melting, the area with a high energy level pressure is created. It is filled with sodium vapor. The progressive expansion of the energy area in the coolant leads to an increase of stress-strain state level of the reactor vessel and may lead to its destruction. The reactor facility must save integrity, provide localization of consequences of beyond design basis accident inside of the pressure vessel and avoid dangerous radiation effects on personnel of nuclear power station and the environment in these conditions. Current Lagrangian formulation is used for the description of coolant motion and structural elements of the reactor. The equation of motion derives from the balance of virtual capacity. Equations of the theory plastic flow are used in physical relations for metals. Deviatoric stress components are assumed to be equal to zero in the coolant; and the equation between hydrostatic pressure and density is taken as a state equation of quasiacoustic type. The contact between the coolant with structural elements of the reactor is simulated by the conditions of non-penetration. The problem solution is based on the method of moment schema of FEM and explicit finite-difference time integration scheme of the “cross” which are implemented in the computing system “Dynamics 3”. The deformation of the fast reactor vessel is investigated numerically in the ULOF beyond design basis accident. The possibility of localizing effects consequences of beyond design basis accident inside of the pressure vessel of the reactor is analyzed.

Текст научной работы на тему «Математическое моделирование развития запроектной аварии внутри корпуса реактора на быстрых нейтронах»

Баженов В.Г., Жестков М.Н., Замятин В .А., Кибец А.И. Математическое моделирование развития запроектной аварии внутри корпуса реактора на быстрых нейтронах // Вестник ПНИПУ. Механика. - 2015. - № 3. - С. 5-14. DOI: 10.15593/peim.mech/2015.3.01

Bazhenov V.G., Zhestkov M.N., Zamyatin V.A., Kibets A.I. Mathematical modeling of development of beyond design basis accident in integral fast-neuron reactors vessel. PNRPU Mechanics Bulletin. 2015. No. 3. Рр. 5-14. DOI: 10.15593/perm.mech/2015.3.01

ВЕСТНИК ПНИПУ. МЕХАНИКА № 3, 2015

PNRPU MECHANICS BULLETIN

http://vestnik.pstu.ru/mechanics/about/inf/

DOI: 10.15593/perm.mech/2015.3.01 УДК 539.3

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

В.Г. Баженов1, М.Н. Жестков1, В.А. Замятин2, А.И. Кибец1

1 Нижегородский государственный университет им. Н.И. Лобачевского, Нижний Новгород, Россия АО «ОКБМ Африкантов», Нижний Новгород, Россия

О СТАТЬЕ

АННОТАЦИЯ

Получена: 15 июня 2015 г.

Принята: 15 июля 2015 г. Опубликована: 30 сентября 2015 г.

Ключевые слова:

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

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

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

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

© ПНИПУ

© Баженов Валентин Георгиевич - доктор физико-математических наук, профессор, bazhenov@mech.unn.ru. Жестков Максим Николаевич - аспирант, e-mail: zhestkov@mail.ru.

Замятин Владимир Аркадьевич - кандидат технических наук, начальник бюро, e-mail: silaev@okbm.nnov.ru. Кибец Александр Иванович - доктор физико-математических наук, профессор, e-mail: kibec@mech.unn.ru.

Valentin G. Bazhenov - Doctor of Physical and Mathematical Sciences, Professor, bazhenov@mech.unn.ru Maxim N. Zhestkov - Postgraduate Student, e-mail: zhestkov@mail.ru

Vladimir A. Zamyatin - PhD in Technical Sciences, Head of Design Department, e-mail: silaev@okbm.nnov.ru Alexander I. Kibets - Doctor of Physical and Mathematical Sciences, Professor, e-mail: kibec@mech.unn.ru

5

Bazhenov V.G., Zhestkov M.N., Zamyatin V.A., Kibets A.I. / PNRPU Mechanics Bulletin 3 (2015) 5-14

MATHEMATICAL MODELING OF DEVELOPMENT OF BEYOND DESIGN BASIS ACCIDENT IN INTEGRAL FAST-NEURON REACTORS VESSEL

V.G. Bazhenov1, M.N. Zhestkov1, V.A. Zamyatin2, A.I. Kibets1

Lobachevsky State University of Nizhni Novgorod, Nizhni Novgorod, Russian Federation

2

JSC “Afrikantov OKBM”, Nizhni Novgorod, Russian Federation

ARTICLE INFO

Received: 15 June 2015 Accepted: 15 July 2015 Published: 30 September 2015

Keywords:

reactor, beyond the design basis accident, hydraulic pressure, finite element method, strength

ABSTRACT

The process of nonstationary deformation of construction of the fast reactor with liquid metal coolant under postulated ULOF beyond design basis accident is examined. This type of accident includes core melting caused by disconnection of the main circulation pump arrangements of the primary circuit with the associated failure of the emergency protection. As a result of core melting, the area with a high energy level pressure is created. It is filled with sodium vapor. The progressive expansion of the energy area in the coolant leads to an increase of stress-strain state level of the reactor vessel and may lead to its destruction. The reactor facility must save integrity, provide localization of consequences of beyond design basis accident inside of the pressure vessel and avoid dangerous radiation effects on personnel of nuclear power station and the environment in these conditions.

Current Lagrangian formulation is used for the description of coolant motion and structural elements of the reactor. The equation of motion derives from the balance of virtual capacity. Equations of the theory plastic flow are used in physical relations for metals. Deviatoric stress components are assumed to be equal to zero in the coolant; and the equation between hydrostatic pressure and density is taken as a state equation of quasiacoustic type. The contact between the coolant with structural elements of the reactor is simulated by the conditions of non-penetration. The problem solution is based on the method of moment schema of FEM and explicit finite-difference time integration scheme of the “cross” which are implemented in the computing system “Dynamics 3”.

The deformation of the fast reactor vessel is investigated numerically in the ULOF beyond design basis accident. The possibility of localizing effects consequences of beyond design basis accident inside of the pressure vessel of the reactor is analyzed.

© PNRPU

Оценка радиационной безопасности окружающей среды и здоровья человека является одним из ведущих факторов, определяющих перспективы развития ядерной энергетики [1-4]. В соответствии с нормативными документами [5] неотъемлемой частью обоснования безопасности любой атомной станции (АС) является анализ запроектных аварий (ЗА). В частности, энергетический сценарий развития запроектной аварии ULOF [6, 7] предусматривает обесточивание АС с потерей источников энергоснабжения, включая резервные, с одновременным несрабатыванием всех средств воздействия на реактивность. В этих условиях конструкции интегральных реакторов на быстрых нейтронах типа БН должны обеспечивать герметичность силового корпуса.

Запроектная авария характеризуется образованием в активной зоне реактора зоны энерговыделения (ЗЭ) с повышенным уровнем давления, заполненной парами натрия. От зоны энерговыделения начинают распространяться волны сжатия, способные вызвать большие смещения теплоносителя, упругопластические деформации или разрушения в конструктивных элементах реактора. Важнейшим аспектом обеспечения безопасности в рассматриваемой аварийной ситуации становится анализ динамической прочности корпуса реактора и внутриреакторного оборудования при их взаимодействии с теплоносителем. Методы и программные коды для моделирования теплофизических явлений в ядерных энергетических установках изложены в [7-11]. Ниже приводятся методика и результаты численного решения трехмерной задачи упругопластического деформирования корпуса реактора на быстрых нейтронах при гидродинамических воздействиях в условиях запроектной аварии ULOF.

6

Баженов В.Г., Жесткое М.Н., Замятин В А., Кибец А.И. / Вестник ПНИПУ. Механика 3 (2015) 5-14

Рис. 1. Расчетная схема (вид сверху)

Рис. 2. Расчетная схема (вид сбоку)

Расчетная схема реактора типа БН представлена на рис. 1, 2, где цифрами обозначены: 1 - корпус реактора; 2 - большая поворотная пробка; 3 - центральная колонна; 4 -зона энерговыделения (ЗЭ); 5 - отражатель нейтронов; 6 - опорный пояс; 7 - промежуточный теплообменник (ПТО); 8 - главный циркуляционный насос (ГЦН); 9 - элеватор. Элементы конструкции реактора выполнены из стали 08XI8H9, физико-механические свойства которой при температуре 450 °С приведены в [6]. Предполагалось, что момент времени t = 0 соответствует завершению фазы разгона реактора и образованию сферической зоны энерговыделения. Источник энерговыделения располагался в центре активной зоны. Изменение во времени давления P в зоне энерговыделения, рассчитанное с учетом теплообмена с холодными металлоконструкциями корпуса, представлено на рис. 3, где Р0 = 3 МПа - максимальное значение давления [6].

Рис. 3. Зависимость давления в теплоносителе на границе зоны энерговыделения от времени

7

Bazhenov V.G., Zhestkov M.N., Zamyatin V.A., Kibets A.I. / PNRPU Mechanics Bulletin 3 (2015) 5-14

Для описания движения конструкции реактора и теплоносителя применяется текущая лагранжевая формулировка [12-14]. Уравнение движения выводится из баланса виртуальных мощностей [12, 13]:

jo^dV+ jpU^U-dV = j PI5UIIdy + jPq5UIIdy (i, j = й), (1)

an rp гд

где Ui - компоненты вектора скорости перемещения в общей декартовой системе координат Х;; Qj, sj - компоненты тензоров напряжений Коши и скоростей деформаций (симметричной части градиента скорости перемещений); р - плотность; Ptq - контактное давление на границе Tq взаимодействия теплоносителя и конструктивных элементов реактора; P - нагрузка, распределенная на поверхности зоны энерговыделения Гp; П -исследуемая область (теплоноситель, корпус реактора и внутриреакторное оборудование); 5sj, 5 U - вариации Sj, Uj (на поверхности с заданными кинематическими граничными условиями 5Ui = 0); точка над символом означает частную производную по

времени t ; по повторяющимся индексам ведется суммирование. Компоненты тензора скорости деформаций определяются в метрике текущего состояния:

Sj = (U,, +Uj„)/2 (i,j = й), U,j =dU,/dXj , X, = X\ =0 + jU,dt. (2)

0

В качестве физических соотношений для металлов применяются уравнения теории пластического течения [15, 16].

t V о -V '-у ту • V - V - /о

Q'j =Qj + Q 5j, Q = -3KS , S = Sii /3, s'e = S -SV5 -Sp Sp = 0 D Q' = 2Gs'e

bij bij b bj’ bn ^Гу ^^ij’

sp = Xdf/ ',, f = Q JQJ--Ql =0; qт = qт (x^ x =

3

s Jp s JpA.

(3)

Здесь Q'j, s 'j, qv, sv - девиаторные и шаровые компоненты тензоров напряжения и скоростей деформаций; sp - скорости пластических деформаций; G, К - модули сдвига и объемного сжатия; 5iJ- - символы Кронекера; Dj - производная Яуманна [17]: djq'j =q'j-Q'kWk--QjkWik, где wj ={Uhj-Uu)/2; f- поверхность текучести Мизеса; Q.J, - динамический предел текучести; X - параметр, тождественно равный нулю при упругом деформировании и определяемый при упругопластическом деформировании из условия прохождения мгновенной поверхности текучести через конец вектора догрузки.

В настоящей работе исследуется нестационарная стадия запроектной аварии, характеризующаяся максимальным уровнем давления в зоне энерговыделения. Длительность нестационарной стадии, по данным численного решения задачи в осесимметричной постановке [6], составляет примерно 10 мс. Отклонение температуры от 450 °С незначительно. Скорость деформаций в конструктивных элементах реактора не превышает 10 с-1. Поэтому влияние изменения скорости деформаций и температуры на деформационные характеристики материала не учитывается.

8

Баженов В.Г., Жесткое М.Н., Замятин В А., Кибец А.И. / Вестник ПНИПУ. Механика 3 (2015) 5-14

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

Контакт между теплоносителем и конструктивными элементами реактора моделировался условиями непроникания по нормали и свободного скольжения вдоль касательной к поверхности контакта [18]:

= u2, = -ql, q, = qi = ° 1 = ^4. (4)

Здесь n, t), x2 - единичные векторы локального, ортогонального базиса; n - вектор нормали к поверхности контакта; т1, т2 ортогональны n; нижний индекс i означает проекцию

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

только для сжимающих сил (q'n, nl) < 0. Система уравнений (1)-(4) дополняется начальными условиями и кинематическими граничными условиями.

Решение трехмерной нелинейной задачи динамики основано на методе конечных элементов и явной конечно-разностной схеме интегрирования по времени типа «крест» [16, 18-21]. Деформируемая конструкция заменяется лагранжевой сеткой, состоящей из

8-узловых конечных элементов (КЭ). В узлах сетки определяются ускорения {U J, скорости {U} и перемещения {U) в общей системе координат {X} = {X1X2X3 }Г . Предполагается, что в конечных элементах сдвиговые и изгибные деформации малы, а смещения и углы поворота КЭ как жесткого целого произвольны. В каждом конечном элементе

вводится локальный базис {х} = {х1 х2х3}г, отслеживающий его вращение как жесткого целого [20]. Конечный элемент с помощью полилинейного изопараметрического преобразования отображается на куб -1 <5, < 1 (i = 1, 3):

X, =±x(Nt(% 15253), N = (1+^1/5k)(1+5i/52)(1+53/53)/8. (5)

k=1

Здесь х,, 5k - координаты узлов в базисах х, 5 ; Nk - функции формы. Для предотвращения развития мод нулевой энергии компоненты скорости деформаций siJ. в КЭ аппроксимируются линейными функциями

S ij = S J + S)} + S ]5 2 + S j 53. (6)

По аналогии с теорией оболочек s0 - значения компонент скорости деформаций в центре КЭ далее будут называться безмоментными составляющими, а их градиент sk = дё j / d5k = const - моментными составляющими. Чтобы не завышать сдвиговую жесткость элемента, в (6) учитываются только компоненты sk, соответствующие изгибающим и крутящим моментам в теории оболочек [21]. На основе (6) разработано семейство конечных элементов для моделирования сложных составных конструкций, включающих деформируемые среды, массивные тела и оболочки [16].

9

Bazhenov V.G., Zhestkov M.N., Zamyatin V.A., Kibets A.I. / PNRPU Mechanics Bulletin 3 (2015) 5-14

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

[M ]U ] = i^!. (7)

где [M] - диагональная матрица масс; |U]; {F] - векторы, составленные из ускорений

узлов КЭ-сетки и результирующих узловых сил в общей системе координат. Система обыкновенных дифференциальных уравнений (7) интегрируется по явной конечноразностной схеме типа «крест».

Численное определение контактного давления в зонах взаимодействия деформируемых тел и статически эквивалентных ему сил в узлах КЭ-сетки осуществляется из условий непроникания и законов сохранения массы и количества движения [18]. Изложенная методика реализована в рамках вычислительной системы (ВС) «Динамика-3», аттестованной в научно-техническом центре по ядерной и радиационной безопасности (Аттестационный паспорт программного средства. Регистрационный паспорт аттестации ПС № 325 от 18.04.2013) и Госстандарте РФ (Программный продукт «Пакет прикладных программ для решения трехмерных задач нестационарного деформирования конструкций, включающих массивные тела и оболочки, «Динамика-3» (111111 «Динамика-3»): Сертификат соответствия Госстандарта России № РОСС RU.ME20.H00338). В ходе аттестации ВС «Динамика-3» были выполнены исследования точности, сходимости и устойчивости реализованных в ней математической модели и методики решения [16, 18, 20]. Результаты исследований подтвердили их эффективность в рассматриваемом классе задач.

При численном анализе напряженно-деформированного состояния реактора типа БН в условиях запроектной аварии элеватор, главные циркуляционные насосы и промежуточные теплообменники моделировались недеформируемыми цилиндрическими оболочками с одинаковыми геометрическими параметрами. С учетом циклической симметрии рассматриваемой конструкции (рис. 1) и условий нагружения в качестве расчетной области выбрана 1/20 часть реактора. Общее количество узлов дискретной модели расчетной области составило 13325. Из них 1212 узлов приходилось на корпус реактора и обечайку защиты. Применение относительно грубой дискретизации конструкции стало возможным благодаря принятой аппроксимации скорости деформаций и напряжений [16], которая является, по существу, численной реализацией шестимодального варианта теории типа Тимошенко [21] и позволяет моделировать динамику тонкостенных оболочек на сетках с одним элементом по толщине. Как показало сопоставление результатов численного решения осесиметричной задачи без учета промежуточных теплообменников и главных циркуляционных насосов реактора в трехмерной и осесимметричной постановках [6], применяемая методика [16, 18, 20] обеспечивает приемлемую точность решения на такой конечно-элементной сетке.

Результаты численных исследований представлены на рис. 4-6 в виде графиков временной зависимости: а) контактной силы F = F /(P0SP), действующей на днище центральной колонны и поворотной пробки (SP - площадь днища центральной колонны

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

10

Баженов В.Г., Жесткое М.Н., Замятин В А., Кибец А.И. / Вестник ПНИПУ. Механика 3 (2015) 5-14

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

Рис. 4. Изменение во времени контактной Рис. 5. Изменение во времени скорости

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

колонны и поворотной пробки пробки

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

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

На рассмотренном интервале времени ГЦН и ПТО существенно снижают гидродинамическую нагрузку на корпус реактора и интенсивность напряжений в корпусе. Так, максимальное значение интенсивности напряжений, достигаемое в среднем по вертикали поперечном сечении цилиндрической обечайки, уменьшилось в три раза. Влияние ГЦН/ПТО на гидродинамическую нагрузку на днище центральной колонны несущественно. Максимальное значение интегральной гидродинамической нагрузки F(t) = Fk + Fp, действующей на днище центральной колонны Fk и поворотной пробки Fp, меняется мало. Это связано с тем, что возрастание нагрузки на днище поворотной пробки, вызванное волнами, отраженными от ГЦН/ПТО, происходит в промежуток времени, когда нагрузка, действующая на днище центральной колонны со стороны зоны энерговыделения, уже убывает.

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

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

Влияние внутрибакового оборудования на гидродинамическую нагрузку на днище центральной колонны и поворотную пробку проявляется на нестационарной стадии развития запроектной аварии (8 мс < t < 10 мс ). Поэтому исследование более длительной стадии развития запроектной аварии допустимо проводить в двумерной (осесимметричной) постановке.

11

Bazhenov V.G., Zhestkov M.N., Zamyatin V.A., Kibets A.I. / PNRPU Mechanics Bulletin 3 (2015) 5-14

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

Выводы

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

«Исследование выполнено за счет гранта Российского научного фонда (проект № 15-19-10039). Математическая модель и методика решения разработаны в Нижегородском государственном университете им. Н.И. Лобачевского. Численный анализ нестационарного деформирования реактора на быстрых нейтронах с жидким теплоносителем в условиях запроектной аварии типа ULOF выполнен совместно с сотрудниками АО „ОКБМ Африкантов“».

Библиографический список

1. Концептуальные положения стратегии развития ядерной энергетики России в 21 веке / Ав-ронин Е.Н. [и др.]; ОАО «НИКИЭТ». - М., 2012. - 62 с.

2. О стратегии ядерной энергетики России до 2050 года / П.Н. Алексеев, В.Г. Асмолов, А.Ю. Гагаринский, Н.Е. Кухаркин, Ю.М. Семченков, В.А. Сидоренко, С.А. Субботин, В.Ф. Цибульский, Я.И. Штромбах; НИЦ «Курчатовский институт». - М., 2012. - 144 с.

3. Крышев И.И., Рязанцев Е.П. Экологическая безопасность ядерно-энергетического комплекса России. - М.: ИздАТ, 2010. - 496 с.

4. Final Report of the International Mission on Remediation of Large Contaminated Areas Off-Site the Fukushima Dailchi NPP. 7-15 October 2011. - Vienna: IAEA, 2011. - 80 p.

5. Общие положения обеспечения безопасности атомных станций ОПБ-88/97, НП-001-97 (ПНАЭ Г-01-011-97) - М.: Энергоатомиздат, 1998.

6. Расчет напряженно-деформированного состояния корпусов реакторов типа БН в условиях гипотетической аварии / В.Г. Баженов [и др.] // Машиноведение. - М.: Наука, 1985. - № 3. -С. 62-68.

7. Волков А.В., Кузнецов И.А. Усовершенствованная модель кипения натрия для анализа аварий в быстром реакторе // Известия вузов. Ядерная энергетика. - 2006. - № 2. - С. 101-111.

8. Ашурко Ю.М., Волков А.В., Раскач К.Ф. Разработка программных модулей для расчета запроектных аварий в быстрых реакторах с учетом пространственно-временной кинетики // Атомная энергия. - 2013. - Т. 114, № 2. - С. 63-67.

9. Использование кода СОКРАТ-БН для расчета задач с кипением натрия и воды в элементах ЯЭУ / А.А. Бутов [и др.] // Известия высших учебных заведений. Физика. - 2012. - Т. 55, № 2/2. -С.137-141.

12

Баженов В.Г., Жесткое М.Н., Замятин В А., Кибец А.И. / Вестник ПНИПУ. Механика 3 (2015) 5-14

10. Ninokata H. A Comparative Overview of Thermal Hydraulic Characteristics of Integrated Primary System Nuclear Reactors // Journal of Nuclear Engineering and Technology. - 2006. - Vol. 38. -No. 1. - Р. 33-44.

11. Experimental verification of the fast reactor safety analysis code SIMMER-III for transient bubble behavior with condensation / K. Morita, T. Matsumoto, K. Fukuda, Y. Tobita, H. Yamano, I. Sato // Nuclear Engineering and Design. - 2008. - Vol. 238. - Р. 49-56.

12. Поздеев А.А., Трусов П.В., Няшин Ю.И. Большие упругопластические деформации: теория, алгоритмы, приложения. - М.: Наука, 1986. - 232 с.

13. Коробейников С.Н. Нелинейное деформирование твердых тел. - Новосибирск: Изд-во СО РАН, 2000. - 262 с.

14. Bathe K.-Y. Finite element procedures. - New Jersey: Upper Saddle River «Prentice Hall», 1996. - 1037p.

15. Волков И.А., Коротких Ю.Г. Уравнения состояния вязкоупругопластических сред с повреждениями. - М.: Физматлит, 2008. - 424 с.

16. Верификация конечно-элементного решения трехмерных нестационарных задач упругопластического деформирования, устойчивости и закритического поведения оболочек / А.А. Артемьева [и др.] // Вычислительная механика сплошных сред. - 2010. - Т. 3, № 2. - С. 5-14.

17. Аннин Б.Д., Коробейников С.Н. Допустимые формы упругих законов деформирования в определяющих соотношениях упругопластичности // Сиб. журн. индустр. мат. - 1998. - Т. 1, № 1. - С. 21-34.

18. Баженов В.Г., Кибец А.И., Цветкова И.Н. Численное моделирование нестационарных процессов ударного взаимодействия деформируемых элементов конструкций // Проблемы машиностроения и надежности машин. - 1995. - № 2. - С. 20-26.

19. Голованов А.И., Тюленева О.Н., Шигабутдинов А.Ф. Метод конечных элементов в статике и динамике тонкостенных конструкций. - М.: Физматлит, 2006. - 391 с.

20. Конечно-элементное моделирование упругопластического выпучивания незамкнутых сферических оболочек при сжатии / В.Г. Баженов, А.А. Артемьева, Е.Г. Гоник, А.И. Кибец, Д.В. Шошин, Т.Г. Федорова // Проблемы прочности и пластичности. - 2012. - № 74. - С. 84-91.

21. Метод конечных элементов в механике твердых тел / под ред. А.С. Сахарова и И. Альтенбаха. - Киев: Вища школа; Лейпциг: ФЕБ Фахбухферлаг, 1982. - 480 с.

22. Нормы расчета на прочность оборудования и трубопроводов атомных энергетических установок (ПНАЭ Г-7-002-86). - М.: Энергоатомиздат, 1989.

References

1. Avrorin E.N. [et al.] Kontseptual’nye polozheniia strategii razvitiia iadernoi energetiki Rossii v 21 veke [Conceptual provisions of the strategy of development of nuclear energy in Russia in the 21st century]. Moscow, 2012. 62 p.

2. Alekseev P.N., Asmolov V.G., Gagarinskii A.Iu., Kukharkin N.E., Semchenkov Iu.M., Sidorenko V.A., Subbotin S.A., Tsibulskii V.F., Shtrombakh Ia.I. O strategii iadernoi energetiki Rossii do 2050 goda [On the Strategy of Russian nuclear power by 2050]. Moscow: National Research Centre “Kurchatov Institute”, 2012. 144 p.

3. Kryshev I.I., Riazantsev E.P. Ekologicheskaia bezopasnost’ iaderno-energeticheskogo kompleksa Rossii [Environmental safety of nuclear power complex of Russia]. Moscow: IzdAT, 2010. 496 p.

4. Final Report of the International Mission on Remediation of Large Contaminated Areas Off-Site the Fukushima DaiIchi NPP. 7-15 October. Vienna: IAEA, 2011. 80 p.

5. Obshchie polozheniia obespecheniia bezopasnosti atomnykh stantsii OPB-88/97, NP-001-97 (PNAE G-01-011-97) [General safety of nuclear power plants OPB-88/97, NP-001-97 (PNAE G-01-011-97)]. Moscow: Energoatomizdat, 1998.

6. Bazhenov V.G., Zhukov V.V., Zamiatin V.A., Kochetkov A.V., Krylov S.V. Raschet napriazhenno-deformirovannogo sostoianiia korpusov reaktorov tipa BN v usloviiakh gipoteticheskoi avarii [Calculation of the stress-strain state of the buildings of fast reactors under hypothetical accident]. Mashinovedenie, 1985, no. 3, pp. 62-68.

13

Bazhenov V.G., Zhestkov M.N., Zamyatin V.A., Kibets A.I. / PNRPU Mechanics Bulletin 3 (2015) 5-14

7. Volkov A.V., Kuznetsov I.A. Usovershenstvovannaia model’ kipeniia natriia dlia analiza avarii v bystrom reaktore [Enhanced sodium boiling model for analysis of accidents in fast reactors]. Izvestiia vysshikh uchebnykh zavedenii. Iadernaia energetika, 2006, no. 2, pp. 101-111.

8. Ashurko Yu.M., Volkov A.V., Raskach K.F. Razrabotka programmnykh modulei dlia rascheta zaproektnykh avarii v bystrykh reaktorakh s uchetom prostranstvenno-vremennoi kinetiki [Development of program modules with space-time kinetics for calculating unanticipated accidents in fast reactors]. Atomic Energy, 2013, vol. 114, no. 2. pp. 77-82.

9. Butov A.A., Vozhakov I.S., Zhigach S.A., Kudashov I.G., Lezhnin S.I., Usov E.V. Ispol’zovanie koda SOKRAT-BN dlia rascheta zadach s kipeniem natriia i vody v elementakh IaEU [Use Code SOCRATES-BN for calculation tasks with boiling water and sodium in the elements of nuclear power units]. Izvestiia vysshikh uchebnykh zavedenii. Fizika, 2012, vol. 55, no. 2. pp. 137-141.

10. Ninokata H. A Comparative Overview of Thermal Hydraulic Characteristics of Integrated Primary System Nuclear Reactors. Journal of Nuclear Engineering and Technology, 2006, vol. 38, no. 1, pp.33-44.

11. Morita K., Matsumoto T., Fukuda K., Tobita Y., Yamano H., Sato I. Experimental verification of the fast reactor safety analysis code SIMMER-III for transient bubble behavior with condensation. Nuclear Engineering and Design, 2008, vol. 238, pp. 49-56.

12. Pozdeev A.A., Trusov P.V., Niashin Yu.I. Bol’shie uprugoplasticheskie deformatsii: teoriia, algoritmy, prilozheniia [The large elastoplastic deformation theory, algorithms, applications]. Moscow: Nauka, 1986. 232 p.

13. Korobeinikov S.N. Korobeynikov S.N. Nelineinoe deformirovanie tverdykh tel [Nonlinear deformation of solids]. Novosibirsk: Sibirskoe otdelenie Rossiiskoi akademii nauk, 2000. 262 p.

14. Bathe K.-Y. Finite element procedures. New Jersey: Upper Saddle River «Prentice Hall», 1996. 1037 p.

15. Volkov I.A., Korotkikh Yu.G. Uravneniia sostoianiia viazkouprugoplasticheskikh sred s povrezh-deniiami [Equations of state viscoelastoplastic media with injuries]. Moscow: Fizmatlit, 2008. 424 p.

16. Artem’eva A.A., Bazhenov V.G., Kibets A.I. [et al.] Verifikatsiia konechno-elementnogo resheniia trekhmernykh nestatsionarnykh zadach uprugoplasticheskogo deformirovaniia, ustoichivosti i zakriticheskogo povedeniia obolochek [Verification of the finite-element solution of 3d non-stationary problems of elasto-plastic deformation, stability and supercritical behavior of shells]. Vychislitelnaia mekhanika sploshnykh sred, 2010, vol. 3, no. 2, pp. 5-14.

17. Annin B.D., Korobeinikov S.N. Dopustimye formy uprugikh zakonov deformirovaniia v opredeliaiushchikh sootnosheniiakh uprugo-plastichnosti [Acceptable forms of elastic laws of deformation in the constitutive relations of elastic-plasticity]. Sibirskii zhurnal industrialnoi matematiki, 1998, vol. 1, no. 1, pp. 21-34.

18. Bazhenov V.G., Kibets A.I., Tsvetkova I.N. Chislennoe modelirovanie nestatsionarnykh protsessov udarnogo vzaimodeistviia deformiruemykh elementov konstruktsii [Numerical simulation of transient processes of shock interaction of deformable elements of constructions]. Problemy mashinostroeniia i nadezhnosti mashin, 1995, no. 2, pp. 20-26.

19. Golovanov A.I., Tiuleneva O.N., Shigabutdinov A.F. Metod konechnykh elementov v statike i dinamike tonkostennykh konstruktsii [The finite element method in statics and dynamics of thin-walled structures]. Moscow: Fizmatlit, 2006. 391 p.

20. Bazhenov V.G. Artemeva A.A. Gonik E.G. Kibets A.I. Shoshin D.V. Fedorova T.G. Konechno-elementnoe modelirovanie uprugoplasticheskogo vypuchivaniia nezamknutykh sfericheskikh obolochek pri szhatii [Finite-element modeling of elastoplastic buckling of non-closed spherical shells loaded in compression]. Problemyprochnosti iplastichnosti, 2012, no. 74, pp. 84-91.

21. Sakharov A.S., Al’tenbakh I. Metod konechnykh elementov v mekhanike tverdykh tel [The method of finite elements in the mechanics of solids]. Kiev: Vishcha shkola; Leiptsig: FEB Fakhbukhferlag, 1982. 480 p.

22. Normy rascheta na prochnost’ oborudovaniia i truboprovodov atomnykh energeticheskikh ustanovok (PNAE G-7-002-86). Moscow: Energoatomizdat, 1989.

14

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