аналитическое решение НДС длинного полого цилиндра под действием терморадиационной нагрузки Текст научной статьи по специальности «Физика»

ЕСТЕСТВЕННЫЕ И ТОЧНЫЕ НАУКИ

УДК 539.3

АНАЛИТИЧЕСКОЕ РЕШЕНИЕ НДС ДЛИННОГО ПОЛОГО ЦИЛИНДРА ПОД ДЕЙСТВИЕМ ТЕРМОРАДИАЦИОННОЙ НАГРУЗКИ

Докт. физ.-мат. наук, проф. КУЛИКОВ И. С., асп. КЛУС С. А.

Белорусский национальный технический университет

В работе предлагается аналитическое исследование напряженно-деформированного состояния (НДС) полого длинного цилиндра при неравномерном нагреве и нейтронном облучении. В основу решения данной задачи легло аналитическое решение аналогичной задачи, полученное С. П. Тимошенко и Дж. Гудьер для термических нагрузок [1].

Рассмотрим длинный полый цилиндр внутренним радиусом Я1 и внешним радиусом Я2, изготовленный из стали ОХ16Н15М3Б под действием внутренней Рв и внешней Рн распределенных нагрузок (рис. 1), поля высоких температур и нейтронного облучения. В итоге данная модель нагружения будет соответствовать модели нагружения оболочки тепловыделяющего элемента (твэла) атомного реактора. Нагрев цилиндра происходит с внутренней поверхности (в случае ядерного реактора нагрев от топливного стержня твэла), и изменение значений температуры, а, следовательно, и значений напряжений и деформаций происходит по толщине цилиндра.

Рис. 1. Сечение цилиндра, нагруженного внешним и внутренним давлениями

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

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

Соотношения между напряжениями и деформациями имеют вид [1, 2]:

8г =

) + авТ + —в5

Е ( -У(сф+с г ))- з

8ф = -1 (Сф -у(сГ +сг)) + авТ +1 в5;

" 3

))- '

=Е(с --ч

Сг + Сф

(1)

) + авТ + —в5

где 8 - полная деформация; вТ - деформации,

5

вызванные термическим расширением; в - то же радиационным распуханием; / = г, ф, г.

Так как w = 0, е г = 0, третье уравнение системы (1) примет вид

сг =у(сг +сф) + ЕI авТ +1 в5

(2)

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

вТ = аТ (г); в5 = 5 (Т (г), Ф),

(3)

где аТ(г), 5(Т(г),Ф) - заданные функции термического расширения и радиационного распухания соответственно; Ф = ф( - интегральный нейтронный поток; ф - плотность нейтронного потока; ^ - время облучения; а - коэффициент линейного расширения, который в общем случае может зависеть от Т.

Если температура Т5 на наружной поверхности цилиндра распределена равномерно, а ее значение изменяется по толщине цилиндра по некоторому закону, то температурное поле Т(г)

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

Т(г) = Т + ^ (Щ2 - г2), 4л

(4)

где Л - коэффициент теплопроводности; ду -объемное тепловыделение.

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

5(Т(г), Ф) = А ^)т -101Т(г) Т(

(г )2

(6)

где А0, т, а, Ь - некоторые константы, определяемые экспериментально.

Решая уравнение (1) относительно с г и Сф, получаем следующие соотношения:

Е

с г =■

1 -V2

Е

8г +V8ф - (1 + V) I ав

1 -V2

V8 г

-вф - (1 + V) I авТ +—в

1

1

(7)

С учетом соотношений Коши, представляющих собой зависимость деформаций от перемещений, перепишем уравнение (7) в следующем виде:

сг =-

1 -V2

Е

+ v--(1 + v) I авТ +

1

1 -V2

ди дг

ди и ( Т 1 „

V — +--(1 + v)I авТ + -в5

дг г I 3

(8)

Подставляя (8) в уравнение равновесия для тел цилиндрической геометрии

дсг

с г - Сф

дг

= 0,

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

д2ы и 1 ды . ( двТ 1 дв5

---+--= (1 + v) I а-+--

дг2 г2 г дг I дг 3 дг

(9)

Соотношение (9) можно переписать также в виде

_д_

дг

1 д(гы)

дг

=(1+1

дг 3 дг

(10)

интегрирование которого дает 1 + V

1 -V

а-

1 11

- I вТЫг н---I в5гёг

г 1 г 3

Щ Щ

Откуда главные напряжения примут вид:

+ С1г + —. (11)

г

Е 1

Сг = -

1 -V г2

Е

а

г 1 г

| вТгёг + — | в5гаТ

Щ Щ

С С2

1 + v^ 1 -2v г2

Е 1

1 -V г2

а

г 1 г |вТгйТ + ^ |в5гйТ

Я1

(12)

Т 1 5 1 Е

авТ + — в5 1+--| , .,

1 -VI 3 I 1 + vl 1 ^ г2 )

С .+ С2

Е ( Т + 1 5

с, =--1 авТ +—в5

г 1 -VI 3

2vEC!

(1 + v)(1 - 2^

где постоянные С\ и С2 определяются из граничных условий:

С г = -в С г = -Рн

при г = Щ1; при г = Щ2.

(13)

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

В качестве конкретного примера рассмотрим длинный полый цилиндр внутренним радиусом ЯI = 0,003 м и внешним радиусом Я2 = = 0,0033 м под действием внутренней Рв = = 7,5 МПа и внешней Рн = 0,1 МПа распределенных нагрузок в условиях объемных термических Т(г) и радиационных 5(Т(г), ф^) деформаций со следующими данными: Е = 2 - 105 МПа; ТЩ1) = 500 (773) °С (К); а = 18,3 - 10-6 град-1.

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

310

3,110

3,2-10"

500 0

1

О 2

0 2 ....... 3 / -500

Г -1-103

-500 * 4 *

МО3 | | 1.5-103

310-

3,1-10-

3,2-10—

в

1

2- ' " "" 2 3 0"

9

4 4

0 -

»

/

Г ' " 1

3-10-

3,1-10-

3,2-10-

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

1,5-10

Оэкв 1-10

500

0

ч Ч4 * ч V,

3 2 4

1

3-10-3 3,1-10-3 3,2-10 г

Рис. 3. Распределение эквивалентных напряжений по радиусу цилиндра при времени эксплуатации 0 (1), 1000 ч (2), 4000 ч (3) и 8000 ч (4)

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

В Ы В О Д

Полученные расчетные данные показывают важное значение применения аналитических ме-

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

Л И Т Е Р А Т У Р А

1. Тимошенко, С. П. Теория упругости / С. П. Тимошенко, Дж. Гудьер. - М.: Наука, 1975. - 576 с.

2. Куликов, И. С. Прочность элементов конструкций при облучении / И. С. Куликов, Б. Е. Тверковкин. -Минск, 1990. - 134 с.

3. Куликов, И. С. Прочность элементов конструкций при облучении / И. С. Куликов, В. Б. Нестеренко, Б. Е. Тверковкин. - Минск, 1990. - 144 с.

4. Клус, С. А. Влияние деформаций ползучести, вызванных тепловым и радиационным воздействием, на НДС тел цилиндрической геометрии / С. А. Клус, И. С. Куликов // Вестник БНТУ. - 2010. - № 2. - С. 64-66.

Поступила 17.06.2010

б

а

г

г

г