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

УДК 539.3

С.М. Шляхов, Э.Ф. Кривулина

ТЕРМОУПРУГОЕ СОСТОЯНИЕ ПРЯМОУГОЛЬНОЙ ПЛАСТИНЫ

БЕЗ ВНУТРЕННЕГО ИСТОЧНИКА ТЕПЛА ПРИ ПОРИСТОМ ОХЛАЖДЕНИИ

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

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

тепла

S-М. Shlyakhov, E.F. Krivoulina

THERMOELASTIC CONDITION OF A RECTANGULAR PLATE

WITHOUT AN INTERNAL HEAT SOURCE AT POROUS COOLING

Numerical data of the rectangular plate intense-deformed condition is presented at its porous cooling with constant porosity. The principle of freedom from communications is used, and thermophysical and mechanical characteristics of materials sensitivity to temperature are considered.

Porous cooling, rectangular plate without an internal source of heat

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

Кроме этого, принимаем температуру твердого тела и температуру жидкости или газа внутри пор одинаковыми. Это допущение также упрощает решение задачи.

Рассмотрим задачу пористого охлаждения плоской не выделяющей тепло пластины толщиной h, выполненной из пористого материала со сквозной (капиллярной) пористостью. Температурное поле внутри пластины считается одномерным. Охлаждающая жидкость, нагнетаемая сквозь пластину в положительном направлении (у), имеет при y = температуру Т0.

На верхнем торце пластинки поддерживается постоянная температура Т2. Боковые поверхности теплоизолированы. Теплофизические и механические характеристики материала зависят от местной температуры и пористости (рис. 1).

Тепловой режим в пластине - стационарный. Весовой расход охлаждающей жидкости в теле пластины G, удельная теплоемкость С, коэффициент теплопроводности X зависят от температуры Т и пористости Р [1, 2].

Перенос тепла в такой пластине можно определить как сумму двух составляющих [3]. Первая

составляющая - это теплопроводность внутри твердого тела, равная - X (1 - P) . Вторая состав-

Лу

ляющая - это теплообмен твердое тело - жидкость. За счет этой составляющей температура жидкости повышается на йТ = йд / ОС , где О = ру , V - скорость движения жидкости в теле.

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

йТ

-X'

йу

йТ й 2Т ,

--------1----------т йу

йу йу

\

= ОСйТ.

у

-

Рис. 1. Схема пористого охлаждения пластины в одномерном поле температур.

Таким образом, дифференциальное уравнение, которому удовлетворяет температура пластины, имеет вид

й2Т * йТ п ... , ОС

= 0 при 0 < у < к , где ^ =-

йу2 ^ йу ’ "^ Х(1 - Р)'

Тепловой баланс, составленный для жидкости, выглядит следующим образом: й 20 ^ й0 йу2 7 йу

Из общего решения уравнения

Т = Су + С2

получим частное решение, которое удовлетворяет граничным условиям Т = Т1 при у = 0, Т = Т2 при у = к:

й ^ йу ОС

и_у_ - ^^ = 0 при - < у < 0, где %г = —.

хf

Т = Т + {Т^(е"'У -1) для 0 < у < к .

Общее решение уравнения (3) имеет вид

0 = с//у + С4.

й0

йТ

Граничные условия для потока жидкости 0 = Т0 при у = -<» , X, — = X5 (1 - Р) — при у = 0.

йу

йу

Частное решение уравнения (3), удовлетворяющее этим условиям, принимает вид

0=То+(Т^ -1)е%,у для -тс- у -0.

Исключая температуру Т1 из уравнений (5) и (7), окончательно получим

Т = Т0 +

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

Для поиска нормальных температурных напряжений в пластинке исходим из посылки, что толщиной к по сравнению с размерами в плане (а, Ь) и краевыми эффектами можно пренебречь. Полагаем также, что для случая капиллярной пористости справедливы прежние уравнения МДТ, принятые для несквозной (кажущейся) пористости. В этом случае с достаточной точностью можно использовать прием, основанный на принципе освобождаемости от связей [4, 5]. Расчетная схема свободной пластины отражена на рис. 2.

При полном закреплении пластинки по краям в ней возникнут напряжения, обусловленные стесненным тепловым расширением [5],

< =<(У) =

- E

1 -Ц

Здесь а

1 Т

а

Т0 - начальная температура бруса, принимаемая чаще за ноль; Е( у) = Е (Т) • Е (Р) - переменный по толщине пластинки модуль Юнга, зависящий от пористости и температуры; ц( у) = ц(Р( у)) - переменный по толщине пластинки коэффициент Пуассона, зависящий от пористости.

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

№ =°,(;) =

a( У) E (У )Т (У) +1

^0 Ei

I —

+

1 -^ І Etht

I

-

a( У) E (У)т (У) У

I

a( У^ (У)т (У) У

І -ц

dy.

dy

Заметим, что в случае однородной пластины, т.е. при E = const, р = const (i = 1, 2, ..., n) из формул (11) вытекают известные зависимости [6]

;=1

= - аЕТ (у) + —1— и'^аЕт (у)йу

* х 1 -ц к

12 V к'2 (

+ ,3Л ч IаЕТ(у)уйу. к (1 - Ц) -к'2

На основании полученных формул проведено исследование полей температур и напряжений прямоугольной пластины, выполненной из пористого железа. Высота пластины к = 0,05 м. Пористость постоянна по высоте сечения Р=0,3. Пластинка может быть свободной или защемленной по контуру.

Взяты следующие охлаждающие агенты: вода, моторное масло, воздух (вода , моторное

масло "Сгоо , воздух і і і).

а б

Рис. 3. а - модуль Юнга (ГН/м2); б - коэффициент Пуассона

Температуры при у = -<» равны Т0 воды и моторного масла 20°С, Т0 воздуха -10°С. Температура горячей поверхности пластины Т2 во всех случаях равна 260°С. Весовой расход жидкостей и газа

п о Кг

О одинаков и равен 2------- .

с * м2

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

Влияние пористости на теплофизические характеристики - температуру и коэффициент линейного расширения - показано на рис. 4.

Рис. 5 отражает напряженное состояние при различных закреплениях пластины.

0,05 0,045 0,04 0,035 0,03 0,025 0,02 0,015 0,01 0,005 0

0 50 100 150 200 250 300

б

температура

Рис. 4. а - температура

h 0,05 -i скользящая заделка

0 -з

00 -2 00 -1 00 100

а

h 0,05 п осевые напряжения за счет теплового расширения

V \

0 -з

00 -2 00 -1 00 0

б

h 0,05 п пластина свободна от связей

0 -2

00 -1 00 100

в

Рис. 5. а - напряжения в случае жесткого закрепления по контуру

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

ЛИТЕРАТУРА

1. Крэйт Ф. Основы теплопередачи / Ф. Крэйт, У. Блэк. М.: Мир, 1983. 512 с.

2. Кашталян Ю. А. Характеристики упругости материалов при высоких температурах / Ю.А. Кашталян. Киев: Наукова думка, 1970. 112 с.

3. Шнейдер П. Инженерные проблемы теплопроводности / П. Шнейдер. М.: Иностр. лит., 1960. 479 с.

4. Тимошенко С.П. Теория упругости / С.П. Тимошенко, Дж. Гудьер; под ред. Г.С. Шапиро. 2-е изд. М.: Наука, 1979. 560 с.

5. Гейтвуд Б.Е. Температурные напряжения / Б.Е. Гейтвуд. М.: Иностр. лит., 1959. 349 с.

6. Справочник по машиностроительным материалам: в 4 т. / под ред. д-ра техн. наук проф. Г.И. Погодина-Алексеева. М.: Машгиз, 1959.

Шляхов Станислав Михайлович -

доктор физико-математических наук, профессор кафедры

«Механика деформируемого твердого тела»

Саратовского государственного технического университета имени Гагарина Ю.А.

Кривулина Эльвира Федоровна -

кандидат технических наук, доцент кафедры «Информационные технологии и прикладная математика»

Саратовского государственного аграрного университета имени Н.И. Вавилова

Статья поступила в редакцию 27.10.11, принята к опубликованию 15.11.11

Stanislav М. Shlyakhov -

Dr. Sc., Professor

Department of Deformable Solid Mechanics, Yu. Gagarin Saratov State Technical University

Elvira F. Krivoulina -

PhD, Associate Professor Department of Information Technology and Applied Mathematics,

N. Vavilov Saratov State Agrarian University

ев