НЕЛИНЕЙНАЯ ЗАДАЧА ТЕПЛОПРОВОДНОСТИ ДЛЯ ТОНКОЙ ОБОЛОЧКИ Текст научной статьи по специальности «Физика»

Электронный журнал «Труды МАИ». Выпуск № 71

www.mai .ru/science/trudy/

УДК 536.2

Нелинейная задача теплопроводности для тонкой оболочки Горюнов A.B.*, Молодожникова Р.Н, Прокофьев А.И.

Московский Авиационный Институт (национальный исследовательский университет), МАИ, Волоколамское шоссе, 4, Москва, А-80, ГСП-3, 125993,

Россия *e-mail:msgor@mail. ru

Аннотация

Рассмотрена задача о нагреве произвольной тонкостенной оболочки плоскопараллельным лучистым тепловым потоком от бесконечно удаленного источника излучения. Между оболочкой и окружающей средой происходит конвективный теплообмен. Теплофизические параметры зависят от температуры. Построены различные асимптотические зависимости для расчета температурного поля оболочки. Эти решения обеспечивают необходимую для практических расчетов точность при использовании для конструкций летательных аппаратов.

Ключевые слова: нелинейная задача теплопроводности; летательные аппараты; лучистый тепловой поток; конвективный теплообмен; асимптотические решения; тонкостенная оболочка

Введение.

В конструкциях летательных аппаратов широкое применение находят тонкостенные оболочки. Одной из актуальных задач является определение несущей способности таких конструкций при различных термических и силовых нагрузках [1;2]. Первым этапом решения несвязанной

квазистатической задачи термоупругости является определение температурного поля конструкции. В случае значительного изменения температуры необходимо учитывать зависимость теплофизических параметров от температуры [3].

Постановка задачи.

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

дС 1 / д {В 1 дС \ , д {А - дС \\ , а(Р0) Л ,лх

ср — = — I — (-А —) + —(-А —)) - аВ1*(Ь - ьс) (1)

и дFо АВ\да^А даг) да2 \В да2П к ^ V су V/

при начальном

и граничном

= Ь0

условиях,где

дь дп

Р

Р=Р0

к0

ВГ

ь0яо

а0Д2

А = Г Я0'

Ч(Ро)

Fо

а°°т

Ч0(Ро) Чт

а1

а1=Т'

ь0^0

ЧтК

а 2

а2=У

а

а

а

0"

Здесь принято, что перепад температуры по толщине оболочки равен нулю, по граничному контуру Г оболочка теплоизолирована и введены

0

Г

обозначения: £0- функция распределения температуры по оболочке, -температура окружающей среды, а°0 и а% - координаты в срединной поверхности оболочки, п - нормаль к граничному контуру, т - время, ф - угол падения лучистого теплового потока на поверхность оболочки, И.0 - толщина оболочки, ц{Ро) и q!o0l - функция лучистого теплового потока и ее максимальное значение, а0 - коэффициент теплоотдачи от оболочки в окружающую среду, с0, р0, X0 и а0 - теплоемкость, плотность, теплопроводность и температуропроводность материала оболочки, с0, р0, Х0, а0 и а0 - значения соответствующих коэффициентов при начальной температуре, А и В -коэффициенты первой квадратичной формулы, Я -наименьший радиус кривизны срединной поверхности оболочки.

Будем считать, что в начальный момент времени температура оболочки равна нулю, температура окружающей среды остается нулевой в рассматриваемом промежутке времени.

Методология решения.

Известно [4,5], что при

Ро < 0,1 и ВС > 10 (2)

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

пренебречь.

При малых временах [4,5]

Ро* = ВСРо < 0,1 (3)

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

В случае выполнения условий (2), (3) и линейной зависимости

ср = с0+ схг (4)

из уравнения теплопроводности (1) получаем

с1Ь2 cosф

с01 + ~Т = — 0

I q(т)dт

Тогда решение задачи принимает вид

1 =-^--. (5)

с1

При произвольной зависимости теплофизических параметров от температуры и выполнении условий (2) из уравнения (1) методом вариации произвольной постоянной получаем формулу для расчета температурного поля оболочки

'^М-В'?^)^^^)* (6)

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

t = C-^f'■0q-(йd!; (7)

к •'0 ср

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

t = C^e-ro*f;;0q(Oeвiíd$ (8)

и

' = (9)

(последнее не учитывает конвективный теплообмен).

Если нормальная составляющая лучистого теплового потока в соответствующей части оболочки невелика (угол падения потока на поверхность оболочки достаточно тупой), то максимальные значения температур могут быть таковы, что теплофизические параметры можно считать постоянными. Тогда для расчетов температурных полей в данной зоне можно использовать решение линейной задачи (8) или (9).

Частные случаи.

Рассмотрим цилиндрическую оболочку, на которую перпендикулярно ее оси падает лучистый тепловой поток (рис. 1):

Я

Рис. 1. Нагрев цилиндрической оболочки.

Зависимость для косинуса угла падения принимает вид

cos ^ = cos ац (j — |а|) , —п<а<п,

m f1' ^ - 0' где ц =1.0' ^ < 0.

Приведем аналогичные зависимости для конической оболочки (рис. 2):

Рис. 2. Нагрев конической оболочки.

cos^ = cos — cos aq — \a\j, —n<a<n и для сферической оболочки (рис. 3):

Рис. 3. Нагрев сферической оболочки.

cos ^ = cos^ cos aq ( — — \a\l, 0 <p<n, —n<a<n.

Заключение.

Геометрические характеристики оболочек, применяемых в конструкциях летательных аппаратов, а так же условия их эксплуатации обеспечивают достаточно широкие временные промежутки, на которых построенные

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

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

1.Афанасьев П.П., Голубев И.С., Лавочкин С.Б., Новиков В.Н., Парафесь С.Г., Пестов М.Д., Туркин И.К. Под редакцией Голубева И.С., Туркина И.К. Беспилотные летательные аппараты. Основы устройства и функционирования -М.: МАИ, 2010, 654с.

2.Туркин И.К. Проектирование тонкостенных конструкций ЛА, функционирующих в экстремальных условиях. - М.: МАИ, 2000, 304с.

3.Безухов Н.И., Баженов В.Л., Гольденблат Н.И., Николаенко Н.А., Синюков А.М. Расчеты на прочность, устойчивость и колебания в условиях высоких температур (под редакцией Гольденблата Н.И.). - М.: Машиностроение, 1965, 567с.

4.3ино И.Е., Тропп Э.А. Асимптотические методы в задачах теории теплопроводности. - Л.: ЛГУ, 1978, 224с.

5.Горюнов А.В., Клименко Б.М., Румянцев Б.П., Самарин А.В Теоретико-экспериментальное исследование температурных полей в тонкостенных конструкциях при неравномерном нагреве. В сб.: Температурные задачи и устойчивость пластин и оболочек. Саратов: Издательство Саратовского Университета, 1988, с. 12-14.

6.Горшков А.Г., Горюнов А.В., Либерзон Р.Е. Односторонний нагрев цилиндрической оболочки. - Мат. методы и физ.-мех. поля, 1982, вып. 16, с. 5255.