CC BY
138
УДК 536.2: 678.4.04
ИНЖЕНЕРНЫЙ АНАЛИЗ НЕСТАЦИОНАРНОЙ ТЕПЛОПРОВОДНОСТИ МНОГОСЛОЙНОЙ ПЛАСТИНЫ
Д.К. Агишева, В.М. Шаповалов
Кафедра "Технологические машины и оборудование", Волжский политехнический институт Волгоградского государственного технического университета
Представлена членом редколлегии профессором В. И. Коноваловым
Ключевые слова и фразы: метод Бубнова-Галеркина; многослойная пластина; нестационарная теплопроводность; процесс вулканизации; теплофизиче-ские свойства.
Аннотация: Предложен способ приближенного решения задачи нестационарной теплопроводности многослойной пластины. Теплофизические свойства слоев постоянны и между слоями имеет место идеальный тепловой контакт. Сопоставление полученного решения с точным показало хорошее их соответствие.
Из соображений прочности стенка надувного резинотканевого амортизатора выполняется многослойной. Она состоит из двух слоев резины толщиной 3 и 6 мм. Между слоями резины находится намотка обрезиненного корда толщиной 15 мм. Намотка корда составляет значительную часть толщины стенки, поэтому тепловой эффект вулканизации не учитываем. Между слоями материала идеальный тепловой контакт.
Теплофизические свойства резины (Х=0,163 Вт/(м-К), рс=1,66-Ш6 Дж/(м3-К)) отличаются от свойств волокна корда (Х=0,31 Вт/(м-К), рс=2,37-Ш6 Дж/(м3-К)) [1, 2], где X - коэффициент теплопроводности, рс - объемная теплоемкость. Задача состоит в определении времени прогрева стенки в процессе вулканизации. Продолжительность тепловой обработки определяет энергозатраты процесса.
В процессе вулканизации пар конденсируется и отдает тепло непосредственно внешней и внутренней стенкам амортизатора. Ввиду низкого коэффициента теплопроводности резины и высокого коэффициента теплоотдачи при конденсации число Био велико, так что можно на свободных поверхностях принять граничное условие первого рода.
Имеются точные решения задач теплопроводности многослойной стенки [3, 4], а также приближенные соотношения для эффективной теплопроводности многослойных стенок и композиционных материалов [5, 6]. Однако для инженерной физики приближенные решения сохраняют свою актуальность.
В данной работе дается приближенное решение задачи нестационарной теплопроводности многослойной пластины, пригодное для инженерной оценки, в частности, при расчете продолжительности процесса вулканизации.
В общем случае имеем пластину, состоящую из п различающихся по тепло-физическим свойствам слоев (рис. 1). Учитываем только поперечную теплопроводность. Толщина 1-го слоя - 5,, его теплофизические характеристики X, и рс,-.
n
Общая толщина пластины
PCj Pc2 Pñn
Xl
_ 5, _ _ 5, _ _ 5ÍÍ __
5
Рис. 1 Расчетная схема многослойной пластины
Вместо традиционной записи системы уравнений теплопроводности для каждого слоя и граничных условий четвертого рода (на границах соприкосновения слоев) используем одно уравнение теплопроводности, считая теплофизические свойства функциями поперечной координаты (уравнение теплопроводности с переменными коэффициентами).
Задача описывается уравнением теплопроводности с начальным и граничными условиями
, ,3T д Pc(x = ^ dt дх
\/ \дТ Х(х )— дх
(t > 0; 0 < х < 5);
(1)
где
(X, рс}
t = 0, T = T0; x = 0, T = Tc; x = 5, T = Tc,
(Xj,рCj}, 0 < x < 5j, i-J
(Xi,Pci}, X5j j=1
< x < X5j;
j=j
T0 - начальная температура пластины (однородная по толщине); Тс - температура внешних поверхностей (равная температуре конденсации пара); t - время; Т -температура; х - поперечная координата. Здесь рс (x) и X(x) являются ступенчатыми функциями, сохраняющими постоянное значение в пределах отдельного слоя.
Примем в качестве масштабов теплопроводности, объемной теплоемкости, длины и температуры Хь рсь 5, Т0 -Тс соответственно. Введем безразмерные переменные и параметры
X:
Л,
Xi
ci- =
Pc¿
5 Xj 1 pej
C учетом (2) задача (1) примет вид
с (х = А
3Fo ЗХ
а, л.
0
Т - Тс То - Тс
Fo
XJt Pcj52
Л(Х ^ V 'дХ
(Fo > 0; 0 < х < 1);
(2)
(3)
Fo = 0, 0 = 1; X = 0, 0 = 0; X = 1, 0 = 0,
где
{A,C} =
{1,1}, 0 < X <АЬ
i-1
j=1 j=1
Решение задачи (3) ведем ортогональным методом Бубнова-Галеркина. Приближенное решение, точно удовлетворяющее граничным условиям, находится в семействе функций
©n = Z ak (Fo)-Vk (X),
(4)
k=1
где в качестве базиса взята система функций
(X) = (1 - X)Хк, к = 1,2,...,п.
Ограничимся первым приближением к = 1. Составив невязку для уравнения (3) при © = ©1, потребуем ее ортогональности к базисной функции, т.е.
J L (©1 )(1 - X)XdX = 0,
0
(5)
где L (©1 ) = C (X)
3©1 д
л(х )д©1 V ' дх
1дБо дХ
Подставив выражение (4) в уравнение (5), получим
da
dFo
О J C (X )(1 - X )2 X 2dX-J dX ^л(х )Ц
(1 - X )XdX = 0.
0 0
Чтобы исключить операцию дифференцирования ступенчатой функции Л(Х) во втором слагаемом, выполним интегрирование по частям. При этом получим уравнение
dFOJe(X)(1 - X)2X2dX - X(1 - X)л())©Г +íл(X(1 - 2X) = 0. (6)
d©i
11
„ д©1 "дх
Поскольку Л и С ступенчатые функции, необходимо интервал интегрирования [0, 1] разбить на участки, соответствующие отдельным слоям пластины. Поэтому в развернутой форме уравнение (6) имеет вид
da1
dFO
1 m
J (1 - X)2 X2dX + Z
0 j=2
X:
J Cj (1 - X)2 X2dX
Xj-1
- a1 (1 - X)X(1 - 2X)
|A1
-at Z л j (1 - X )X (1 - 2X ) >
j=2 L j-1
a1
A1
X,-
J(1 - 2X)2dX + £лj J (1 - 2X)2dX
0 j=2 Xj_1
(7) 0.
J
Здесь и ниже для краткости записи обозначим X^ = ^ Дг- .
i=0
Получили дифференциальное уравнение первого порядка для функции аь Решение уравнения (7) имеет вид
О = О) • exp|-D^o |, (8)
где
.3 1 m
] =i3L(10 - 15д1 + 6Aj) + -3Y,cjx3 (10 - 15X + 6X2)
3'- 2
j=2
Xj-1
1 1 3 m Г 1 3
D = 6 - (1 - Ai) Ai (1 - 2Ai) -16 (1 - 2Ai )3 - X Л j (1 - X)X (1 - 2X) + 6 (1 - 2X)3
v
j=2
X
j-1
Для определения постоянной а0 потребуем выполнения начального условия в (3) среднеинтегрально, т.е.
Fo =
0, J©1(a = a0)dX = 1;
О = 6exp|-d^Fq |, (9)
откуда находим а0 = 6.
В частном случае трехслойной пластины (в расчетных формулах необходимо положить m = 3 , Л1 = Л3 = 1, С1 = С3 = 1) решение (8) упрощается и принимает вид
D2
D = 30 + { А3 (10 - 15Al + 6д2) - (Al + Д2)3 - 15(Al + Д2) + 6(( + Д2)2])'
D = 3 + Л23-1 {(1 - 2Ai)i + 2Ai (1 - Ai)] - (1 - 2Ai - 2А2)i + 2(Ai + A2)(1 - Ai - A2)]).
Таким образом, в результате решения удалось свести задачу нестационарной теплопроводности многослойной пластины к теплопроводности «эквивалентной» однослойной.
В случае однослойной пластны достаточно в (9) положить Л2 = 1, С2 = 1. При
D2 лп
этом показатель степени содержит множитель = 10 , что соответствует извест-
Щ
ному результату П.В. Цоя [7].
Найденное решение было сопоставлено с точным решением, полученным методом преобразования Лапласа для симметричной системы тел, состоящей из трех неограниченных пластин [8]. Численный
анализ проводился для условий: = Х.3, Si =63 = 4 мм, 62 = 15 мм, Bi = да.
Учитывались первые 20 членов ряда Фурье. Результаты представлены на рис. 2.
Сплошной линией показано точное решение, соответствующее уравнению
да 2 _1
©1 = -У-CGS(^Kfl 2X /12) ехР(-ЦиKe Kfl1FOi),
П=1 InVn
где Ka = 1,33; Ke = 1,65; Ke = 1,88; Fo1 = 4Fo(26^62 + 1)2; x = 0,0075 м (центр пластины),
-1 -1 , -1 -1 Vn = (1 + KsKeKa 2 ) Sin |ln cos V-nKeKa 2 + Ke (1 + Ke iKeKa 2 ) COS |n Sin |nKeKa 2.
0
Собственные числа цп определялись из уравнения
Ketgц. tg(^KeKfl2) = 1
методом простой итерации в универсальной математической системе Mathcad 2000.
Как видно из рисунка, хорошая сходимость приближенного и точного решений начинается с Foi = 0,1, т.е. в регулярном режиме теплообмена. Кинетика вулканизации внутренних слоев резинового изделия определяется регулярной фазой теплообмена.
Применимость математической модели определяется величиной Лi,Ct (i > 1). Точность расчетов возрастает при стремлении Л1 ,Ct к единице. Повысить точность решения можно, если взять второе приближение в (4). Однако при этом резко возрастает громоздкость расчетных формул.
В заключение отметим, что решения задачи нестационарной теплопроводности для слоисто-однородных тел в настоящее время все больше используются для изучения теплофизических свойств композиционных материалов, которые наиболее эффективно применяются в металлургии, авиационной и космической технике.
Благодарность
Авторы признательны д.т.н., профессору ВПИ ВолгГТУ И.Я. Шиповскому за предоставленные материалы и внимание к работе.
Список литературы
1. Пивень А.Н., Гречаная Н.А., Чернобыльский И.И. Теплофизические свойства полимерных материалов: Справочник. - К.: Вища школа, 1976. - 180 с.
2. Пехович А.И., Жидких В.М. Расчеты теплового режима твердых тел. - Л.: Энергия, 1958. - 304 с.
3. Беляев Н.М., Рядно А.А. Методы нестационарной теплопроводности: Учеб. пособие для вузов. - М.: Высш. школа, 1978. - 328 с.
4. Карташов Э.М. Аналитические методы в теории теплопроводности твердых тел: Учеб. пособие. - 2-е изд., доп. - М.: Высш. школа, 1985. - 480 с.
5. Лыков А.В. Тепломассообмен: Справочник. - М.: Энергия, 1978. - 480 с.
6. Дульнев Г.Н., Заричняк Ю.П. Теплопроводность смесей и композиционных материалов. - Л.: Энергия, 1974. - 264 с.
7. Цой П.В. Методы расчета задач тепломассопереноса. - М.: Энергоатом-издат, 1984. - 412 с.
8. Лыков А.В. Теория теплопроводности. - М.: Высшая школа, 1967. - 600 с.
температуры в центре трехслойной пластины от числа Фурье: сплошная линия - точное решение, пунктирная -приближенное
Engineering Analysis of Non-Steady-State Heat Conduction of Multi-Layer Plate
D.K. Agisheva, V.M. Shapovalov
Department " Technological Machines and Equipment", Volzhski Polytechnic Institute, Volgograd State Technical University
Key words and phrases: non-dimensional variables and parameters; Bubnov-Galerkin method; multi-layer plate; initial and final conditions; Non-Steady-State heat conduction; vulcanization process; thermophysical properties; Fourier number.
Abstract: The way of approximate solution of the problem of non-stationary heat conduction of multi-layer plate is offered. Thermophysical properties of layers are constant and there is ideal thermal contact between layers. Comparison of the obtained solution with the exact one showed their good correspondence.
Ingenieurische Analyse der unstationären Wärmeleitfähigkeit der vielschichtigen Platte
Zusammenfassung: Es ist die Methode der annähernden Lösung der Aufgabe der unstationären Wärmeleitfähigkeit der vielschichtigen Platte vorgeschlagen. Die wärmephysikalischen Eigenschaften von den Schichten sind ständig; zwischen den Schichten gibt es den idealen Wärmekontakt. Die Vergleichung der bekommenen Lösung mit der exakten Lösung hat ihre gute Entsprechung gezeigt.
Analyse d'ingénieur de la conductibilité calorifique non-stationnaire de la plaquette multicouche
Résumé: On a proposé le moyen de la résolution approximative du problème de la conductibilité calorifique non-stationnaire de la plaquette multicouche. Les propriétés thermophysiques des couches sont constantes et le contact calorifique idéal a lieu entre ces couches. La comparaison de la résolution reçue avec la résolution exacte a montré leur conformité parfaite.
