электронное научно-техническое издание
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС 77 - 305&9. Государственная регистрация №0421100025. ISSN 1994-040S
Релаксационная стойкость винтовой цилиндрической пружины в условиях нейтронного облучения
77-30569/402924
# 04, апрель 2012 Гусев М. П., Данилов В. Л.
УДК.621
МГТУ им. Н.Э. Баумана gusev.m.p@yandex.ru vldanilov@mail.ru
Введение
Релаксация силы в сжатой пружине - это негативный фактор при эксплуатации тепловыделяющей сборки энергетического реактора, который обусловлен развитием в материале пружины необратимых деформаций ползучести под действием высокой температуры и флакса нейтронов. Релаксационная стойкость пружины - это способность материала, из которого изготовлена пружина, сопротивляться возникновению необратимых деформаций ползучести и падению напряжений в нем. Релаксационная стойкость пружины оценивается, в основном, по скорости падения силы в сжатой пружине, которую принято представлять кривыми релаксации. Основная цель разработанной математической модели - это возможность построения кривых релаксации для пружины, с различными геометрическими параметрами и работающей в различных условиях температурного и нейтронного воздействий. При этом существенным является разработка эффективного алгоритма и программы расчета, удобных в инженерной практике.
Основные соотношения математической модели
Состоянию чистой релаксации соответствует постоянства во времени деформированного состояния. Такое состояние реализуется при фиксированной высоте предварительно сжатой пружины. При этом напряженное состояние во времени изменяется и сила сжатия пружины падает. На рис. 1 представлена схема реализации условий чистой релаксации. Основные геометрические параметры деформированной пружины ясны из рисунка и остаются неизменными во времени.
Рис. 1. Схемы закрепления и нагружения пружины
В поперечном сечении витка пружины возникают внутренние силовые факторы, представленные на рис. 2. Из условий равновесия отсеченной части пружины следуют известные соотношения теории пружин [ 1 ].
М7 = т бш«н--соб«,
г 2 рв •
Му, = т соб«--бш«,
* 2
ду = р соб«,
N = Р бШ«.
(1)
где Мг - крутящий момент, Му - изгибающий момент, Qy - поперечная сила, N2 -
нормальная сила. Здесь также учтен крутящий момент т на торцах пружины, возникающий при запрете их поворота. Моменты М2 и Му, в свою очередь вызваны, изменением кручения и кривизны оси витка пружины Ак и Ах соответственно.
Рис. 1. Определение внутренних силовых факторов в пружине
Из курса дифференциальной геометрии известно, что кручение и кривизна винтовой линии пружины связаны с ее геометрическими параметрами соотношениями
sin2а 2cos2 а
к =-, Х —
(2)
Б Б
В начальный момент времени пружина находится в упругом состоянии, для которого при условии малой кривизны витка пружины справедливы соотношения
Ак — к кг. — ■
03
р
Ах — I - -1 — М.
Р Р0 Е3
(3)
У
где 3р и 3у - полярный и осевой моменты инерции сечения витка пружины;
р и р0 - измененный и начальный радиусы кривизны винтовой линии пружины.
В этом случае однозначная связь между силой поджатия пружины и изменением ее геометрических параметров будет определяться следующим соотношением
Р —
403л cos а ( sin а cos а sin а cos а 1 4Е3у sin а ( cos2 а cos2 а ^
Б
Б
Бп
Б
Б
Б,
(4)
о у
где а0 и Б0 - угол подъема наклона витка и средний диаметр пружины соответствуют
недеформированному состоянию, а а и Б - соответствуют деформированному состоянию.
Напряженное состояние пружины определяется напряжениями: касательными - тгх и т и нормальным - ог.
Компоненты напряжения тхх и связаны с полным касательным напряжением т соотношениями
\т2Х = тсоб(
\ Т • ( (5)
где угол р отвечает полярной системе координат, связанной с сечением.
Интенсивность напряжений для случая совместного кручения и изгиба имеет вид
= л/о2 + 3т2. (6)
Под действием температуры и потока нейтронов в материале пружины развиваются процессы ползучести и как следствие изменяется во времени напряженное состояние. Уравнение состояния ползучести стали с учетом нейтронного облучения имеет вид [2]
( ^ ( гл \\
Е° = а о
Ъ е се
Q,
ф
1 Ф
1 +--ехр
Ф£ I Я Т
V £ V £ а ))
(7)
где ас, пс, Ф, - экспериментальные коэффициенты, зависящие от температуры;
Ф - флакс нейтронов или поток нейтронов, н
см2 ■ с
Я - универсальная газовая постоянная,
Дж ■ К
моль
Та - температура активной зоны энергетического реактора, К,
Qф - энергия активации ползучести,
Дж
моль
Скорости линейной и угловой деформаций ползучести определяются соотношениями теории ползучести [3]
£ С £ С
£с =о, пс = 3^-т. (8)
0е 0е
Условиями релаксации являются равенства скоростей кручения и кривизны винтовой оси пружины нулю.
Алгоритм численного расчета релаксации силы сжатия пружины
Напряженное состояние витка пружины неоднородно по его сечению. Для интегрирования напряжений разбиваем сечение на малые области. На рис. 3 представлен пример дискретизации сечения.
Рис. 3. Разбиение сечения витка пружины Для определения изменения напряжений в пружине вводится соотношение между деформациями:
е' + е„ = е* = оотХ.
(9)
где е'с'1 - деформация ползучести, соответствующая номерам /, 1 области сечения; еге'1 - упругая деформация, соответствующая номерам /, 1 области сечения; е01 - начальная деформация, неизменяемая во времени.
Записывая соотношение (9) в дифференциальной форме и используя конечно-разностное представление по времени получаем следующие зависимости релаксации
о? = &к 1 ЕЦ
" =т'к1 -пк1О Мк
(10)
где Ахк - текущий шаг по времени.
Выражения (10) связывают напряжения на следующем шаге и напряжения на
предыдущем шаге. Скорости линейной и сдвиговой деформаций ползучести и пП—
определяются по напряжениям на предыдущем шаге в соответствии с зависимостями (6), (7), (8). Начальные значения напряжение определяются по упругому состоянию пружины. Интегрируя текущие значения напряжений по сечению определяем величины внутренних силовых факторов, через которые находится текущее значение силы сжатия пружины. Накопление деформаций ползучести определяется соотношениями
где а^ки у'с,]к- накопленные деформации ползучести на предыдущем шаге.
Данные соотношения позволяют определить изменения начальных размеров пружины после ее разгрузки.
На основе приведенного алгоритма расчета разработана программа в среде МайаЬ, позволяющая получать кривые релаксации силы сжатия при различных геометрических параметрах пружины и различных температурных и нейтронных воздействиях.
Результаты расчетов релаксации силы сжатия пружины
Рассматривалась задача многоэтапного поджатия пружины. Первый этап: сжатие на 15 мм, выдержка в течение 4000 ч; второй этап: дополнительное сжатие на 1 мм, общая осадка 15 +1 = 16 мм, время выдержки 4000 ч; третий этап: 16 +1 = 17 мм, 4000 ч; четвертый этап: 17 +1 = 18 мм, 6000 ч. Результаты расчета представлены на рис. 4 -нижняя кривая. По горизонтальной оси на графике отложено время, по вертикальной оси - сила сжатия пружины.
Для проверки работы алгоритма и программы расчета был выполнен расчет этой же задачи с помощью конечно-элементного комплекса Апбуб 12. На рис. 5 представлена конечно-элементная модель пружины. Условия закрепления и нагружения пружины задавались из требования тождественности начальных состояний обеих расчетных моделей. Результаты конечно-элементного расчета представлены на рис. 4 верхней кривой.
Сопоставляя результаты решения задачи по двум моделям, видно, что разработанная программа дает несколько завышенную релаксацию по сравнению с решением задачи в Апбуб 12. Однако общий вид кривых одинаковый и можно утверждать, что алгоритм разработанной программы построен и реализован верно. Завышенное значение релаксации приводит к более безопасной оценке ресурса работы пружины.
(11)
Force
-180
Relaxation, mm: D1 = 15, D2 = 1, D3 = 1, D4 = 1
1 1 1 1
у? f-
к
f
1
-My programm -
1 -Ansys data
1 1 1
-200
-220
-240
-260
-280
-300
-320
-340
2000 4000 6000 8000 10000 12000 14000 Time, h: T1 = 4Q00, T2 = 4000, T3 = 4000, T4 = 6000
16000
18000
Рис. 4. Сравнение результатов решения задачи
Рис. 5. Конечно-элементная модель пружины
Выводы
В настоящей статье изложен инженерный метод расчета релаксации силы в сжатой пружине под действием температуры и потока нейтронов. Разработанная программа в среде программирования Matlab проверялась расчетом пружины при одинаковых условиях в расчетном комплексе Ansys 12. Установлено, что алгоритм разработанной программы построен верно.
Основными достоинствами разработанной программы можно считать:
- скорость решения задачи в разработанной программе превосходит скорость решения задачи в расчетном комплексе АпБуБ 12;
- возможность варьировать геометрические параметры пружины, внешние условия работы пружины;
- возможность учета изменения физико-механических свойств материала пружины.
Библиографический список:
1. Андреева Л.Е. Упругие элементы приборов. -М.: Машиностроение, 1981.- 392 с.
2. Киселевский В.Н. Изменение механических свойств сталей и сплавов при радиационном облучении. - Киев, Наукова Думка, 1977. - 287 с.
3. Бойл Дж., Спенс Дж. Анализ напряжений в конструкциях при ползучести. - М.: Мир, 1986.- 360 с.
electronic scientific and technical periodical
SCIENCE and EDUCATION
_EL № KS 77 - 3Ü56'». .V;II421100025, ISSN 1994-jMOg_
Relaxation resistance of coiled spring under the conditions of neutron irradiation
77-30569/402924
# 04, April 2012 Gusev M.P., Danilov V.L.
Bauman Moscow State Technical University gusev.m.p@yandex.ru vldanilov@mail.ru
The authors formulate a mathematical model for relaxation resistance of a spring working under the conditions of neutron irradiation and high temperature; a method of its analysis was developed. The model and the analysis method are based on the fundamental equations of the creep theory constructed with the use of the finite-difference algorithm. This model is capable of creep and relaxation analysis of a spring based on various initial spring loads and geometrical parameters. Comparison of results of the proposed model with the well-established industrial analysis computer program "Ansys 12" is included in the article.
Publications with keywords: relaxation, analysis of springs, neutron flux, creep Publications with words: relaxation, analysis of springs, neutron flux, creep
References
1. Andreeva L.E. Uprugie elementy priborov [Elastic elements of the devices]. Moscow, Mashinostroenie, 1981. 392 p.
2. Kiselevskii V.N. Izmenenie mekhanicheskikh svoistv stalei i splavov pri radiatsionnom obluchenii [Change of mechanical properties of steels and alloys at radiation exposure]. Kiev, Naukova Dumka, 1977. 287 p.
3. Boyle J.T., Spence J. Stress Analysis for Creep, London, Butterworths Co., 1983. 283 p. (Russ. ed.: Boil Dzh., Spens Dzh. Analiz napriazhenii v konstruktsiiakhpripolzuchesti. Moscow, Mir, 1986. 360 p.).