НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 • 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Метод расчета релаксации контактного взаимодействия между
тепловыделяющим элементом и ячейкой дистанционирующей
решетки
# 10, октябрь 2013
Б01: 10.7463/1013.0654596
Гусев М. П., Данилов В. Л.
УДК.621
Россия, МГТУ им. Н.Э. Баумана gusev.m.p@yandex.ru vldanilov@mail.ru
Введение
Тепловыделяющие элементы (ТВЭлы) и дистанционирующие решетки (ДР), в узлах которых закреплены ТВЭлы посредством соединения с натягом, являются основными элементами тепловыделяющей сборки (ТВС). Поскольку в активной зоне атомного реактора, которая состоит из набора ТВС, действуют интенсивные температурное и нейтронное поля, то в таких условиях происходит ослабление всех нагруженных компонентов сборки, в том числе происходит ослабление или релаксация контактного взаимодействия между ТВЭлом и ячейкой ДР. Релаксация контактного взаимодействия приводит к снижению контактных сил, действующих между ТВЭлом и ячейкой ДР с течением времени, что, в результате, приводит к проскальзыванию ТВЭла вдоль ячейки ДР. Проскальзывание большого числа ТВЭлов в топливной сборке приводит к недопустимому искривлению и снижению жесткости конструкции ТВС. Недопустимые искривление и снижение жесткости ТВС непосредственно влияют на безопасность и надежность реакторной установки, поэтому решение задачи релаксации контактного взаимодействия является актуальной задачей.
Постановка задачи
В рамках настоящей работы автором была поставлена задача разработки метода расчета релаксации контактного взаимодействия между ТВЭлом и ячейкой ДР, имеющая в себе цель определить зависимость контактной силы между ТЭлом и ячейкой ДР от времени. http://technomag.bmstu.ru/doc/654596.html 415
Разработка метода расчета релаксации контактного взаимодействия
Тепловыделяющие элементы представляют собой длинные тонкостенные стержни, в которых расположены таблетки диоксида урана. ТВЭлы объединены в пучок, посредством дистанционирующих решеток рис. 1, [1]. Объединение ТВЭлов в пучок представляет из себя тепловыделяющую сборку, которая является основным элементом активной зоны атомного реактора. Закрепление ТВЭлов в ячейках ДР осуществляется посредством соединений с натягом.
Тепловыделяющие элементы Дистанционирующие решетки
Рис. 1. Пучок тепловыделяющих элементов с дистанционирующими решетками
В результате действия высокой температуры и нейтронного облучения, возникает релаксация контактного взаимодействия между ТВЭлом и ячейкой ДР, рис. 2, [2]. Как уже отмечалось выше, релаксации контактного взаимодействия приводит к проскальзыванию ТВЭла в ячейке ДР. Эти явления снижают жесткость конструкции ТВС и приводят к недопустимому деформированию сборки [3].
Рис. 2. Схема соединения тепловыделяющего элемента и ячейки дистанционирующей
решетки
Разработка метода расчета релаксации контактного взаимодействия между ТВЭлом и ячейкой ДР требуется для оценки ресурса работоспособности ТВЭла и всей топливной сборки.
В настоящем исследовании рассматривается метод решения задачи релаксации контактного взаимодействия на основе метода конечных элементов и альтернирующего метода Шварца, [4, 5].
Метод решения задачи релаксации контактного взаимодействия состоит из следующих этапов.
На первом этапе создаются отдельные конечно-элементные модели оболочки ТВЭла и ячейки ДР. На втором этапе решается контактная задача взаимодействия оболочки ТВЭла и ячейки ДР, [4]. На третьем этапе для каждого конечного элемента двух
тел рассчитывается вектор напряжений {&а1—1У, где индекс а = 1,2 означает
принадлежность вектора напряжений к ячейке ДР или к оболочке ТВЭла, индекс 1 означает шаг по времени и е - принадлежность конкретному конечному элементу. На четвертом этапе по формуле
{¿а . -1}е=в^-1п }е,
где {¿а^У - вектор скорости деформаций ползучести для конечного элемента, В, п -экспериментальные коэффициенты, определяется вектор скоростей деформаций ползучести на промежутке времени ti — ti—1. Далее определяются деформации ползучести на текущем шаге по времени, на основе явной схемы Эйлера:
{¿а, }е={¿а,—1}е+(ь—1 ж,—1}е,
где {¿а,}е, {¿а— 1}е - вектор деформаций ползучести на текущем шаге 1 и на предыдущем шаге 1 — 1.
На пятом этапе формируется вектор узловых сил {/О,. }е для каждого конечного
элемента двух тел, эквивалентный действию деформаций ползучести {¿а. }е по формуле
{/а Уе = | [ ве ]т [ ве }е^е,
г
где [Ве ] - матрица связи деформаций с узловыми перемещениями, [Ое ] - матрица упругих постоянных и Vе - объем конечного элемента.
На шестом этапе составляется глобальный вектор узловых сил [F^, }, эквивалентный действию деформаций ползучести [Бса i}. Формируется суммарный
глобальный вектор узловых сил [F^.}, учитывающий приложенную нагрузку [Fai} и вект°р [F¿}.
Fi } = Fi } + [FI }.
На седьмом этапе решается контактная задача взаимодействия двух тел, где в качестве внешней нагрузки уже выступает суммарный глобальный вектор узловых сил
[Fai}. На основе решения контактной задачи определяются контактные усилия [FC"}
для каждого из контактирующих тел.
Далее последовательно в цикле повторяются этапы решения со второго по шестой и цикл продолжается до тех пор, пока время tt < tmax.
В результате решения задачи в каждый момент времени ' определены контактные усилия [FC"'}. На рис. 3. приведем решение задачи релаксации контактного взаимодействия при температуре 313 oC и при плотности потока нейтронов 1,73 -1018 н / м2с.
2_i_i_i_i_i_i_
О 1000 2000 3000 4000 5000 6000 7000
Время, ч
Рис. 3. Релаксация контактной силы между ТВЭлом и ячейкой ДР для одной области
контакта
Проскальзывание ТВЭла в ячейке ДР возникает в момент времени ts, когда силы трения в области контакта становятся меньше продольной внешней силы, действующей на ТВЭл.
Однако проскальзывание единичного ТВЭла в ячейке ДР не является признаком полной потери работоспособности всей топливной сборки. Полная потеря работоспособности наступает тогда, когда количество ТВЭлов, подвергнутых проскальзыванию, превышает предельно допустимое значение для топливной сборки и жесткость ТВС становится ниже критической. Время, при котором происходит проскальзывание критического количества ТВЭлов, является одной из оценок ресурса работы ТВС.
Выводы
1. Разработан метод расчета релаксации контактного взаимодействия на основе метода конечных элементов и альтернирующего метода Шварца.
2. Предложен метод определения времени проскальзывания единичного ТВЭла в ячейке дистанционирующей решетки, на основе расчета релаксации контактного взаимодействия.
3. Предложен критерий потери работоспособности топливной сборки, который заключается в определении времени проскальзывания критического количества ТВЭлов в узлах дистанционирующих решеток.
Список литературы
1) Семишкин В.П. Расчётно-экспериментальные методы обоснования поведения твэлов и ТВС ВВЭР в аварийных режимах с большой течью из первого контура РУ: дис. ... докт. техн. наук. Подольск, 2007. 360 с.
2) Троянов В.М. Расчетно-экспериментальное обоснование термомеханики активных зон реакторов типа ВВЭР: дис. ... докт. техн. наук. Обнинск, 2003. 237 с.
3) Review of Fuel Failures in Water Cooled Reactors. Vienna: International Atomic Energy Agency, 2010. 191 p.
4) Станкевич И.В., Яковлев М. Е., Си Ту Хтет. Разработка алгоритма контактного взаимодействия на основе альтернирующего метода Шварца // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2011. Спец. вып. "Прикладная математика". С. 134-141.
5) Станкевич И.В., Яковлев М.Е., Си Ту Хтет. Математическое моделирование контактного взаимодействия упругопластических сред // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2012. № 4. Режим доступа: http://technomag.edu.ru/doc/353180.html (дата обращения 26.11.2013).
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-040S
electronic scientific and technical journal
Method of relaxation analysis of contact interaction between a fuel
element and a cell of a spacer grid
# 10, October 2013
DOI: 10.7463/1013.0654596
Gusev M.P., Danilov V.L.
Bauman Moscow State Technical University, 105005, Moscow, Russian Federation
gusev.m.p@yandex.ru vldanilov@mail.ru
Development of a method for calculating relaxation of a contact between a fuel element and a spacer grid cell is a necessary condition for solving complex deformation problems of fuel assembly under intensive temperature and neutron fields. Due to relaxation of contact interaction the fuel element can slide in the spacer grid cell; this results in decreasing stiffness of overall assembly and its deformation. A method for calculating relaxation of contact interaction between a fuel element and a spacer grid cell based on the finite element method and the Schwarz alternating method was proposed in this work. Calculation of relaxation with the use of the proposed method will allow to estimate operation time of a fuel element till its sliding, and also to determine the disability criterion of fuel assembly.
Publications with keywords: FEM, fuel element, the contact interaction relaxation, fuel assembly, spacer grid
Publications with words: FEM, fuel element, the contact interaction relaxation, fuel assembly, spacer grid
References
1. Semishkin V.P. Raschetno-eksperimental'nye metody obosnovaniyapovedeniya tvelov i TVS VVER v avariynykh rezhimakh s bol'shoy tech'yu izpervogo kontura RU. Dokt.
diss. [Calculation and experimental methods of justification of behavior of fuel rods and fuel assemblies of WWER in emergency mode with big leak from the primary circuit of reactor plant. Dr. diss.]. Podol'sk, 2007. 360 p.
2. Troyanov V.M. Raschetno-eksperimental'noe obosnovanie termomekhaniki aktivnykh zon reaktorov tipa VVER. Dokt. diss. [Calculation and experimental justification of thermomechanics of active zones of WWER reactors. Dr. diss.]. Obninsk, 2003. 237 p.
3. Review of Fuel Failures in Water Cooled Reactors. Vienna, International Atomic Energy Agency, 2010. 191 p.
4. Stankevich I.V., Yakovlev M. E., Si Tu Khtet. Razrabotka algoritma kontaktnogo vzaimodeystviya na osnove al'terniruyushchego metoda Shvartsa [Development of algorithm of contact interaction on the basis of Schwarz alternating method]. VestnikMGTUim. N.E. Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2011, Spec. iss. "Prikladnaya matematika" ["Applied mathematics"], pp. 134-141.
5. Stankevich I.V., Yakovlev M.E., Si Tu Khtet. Matematicheskoe modelirovanie kontaktnogo vzaimodeystviya uprugoplasticheskikh sred [Mathematical simulation of contact interaction of elasto-plastic environments]. Nauka i obrazovanie MGTU im. N.E. Baumana [Science and Education of the Bauman MSTU], 2012, no. 4. Available at: http://technomag.edu.ru/doc/353180.html , accessed 26.11.2013.