Т Р А Н С П О Р Т
УДК 621.431.7
ОПРЕДЕЛЕНИЕ ПОВЕРХНОСТИ ОХЛАЖДЕНИЯ РАДИАТОРА
Докт. техн. наук, проф. ЯКУБОВИЧ А. И., асп. ТАРАСЕНКО В. Е.
Белорусский государственный аграрный технический университет
Особая роль в обеспечении заданного температурного режима дизеля при жидкостном охлаждении принадлежит радиаторам, поэтому расчет и выбор геометрических и теплотехнических параметров радиатора - крайне важный этап при проектировании всей системы охлаждения.
Конструктор может вести разработку радиатора проектируемого трактора несколькими способами, например с помощью прототипа, изменяя его геометрические параметры в зависимости от заданной мощности дизеля и компоновочных особенностей трактора или путем расчета параметров радиатора, используя математическую модель. При этом второй способ предпочтителен, так как он позволяет определять направления экономии материалов, обеспечивая лишь заранее заданный запас по тепловой эффективности. Кроме того, такой подход обеспечивает большую вероятность получения рабочего образца изделия за меньший промежуток времени и с наименьшими затратами.
Определение теплоотдачи от радиатора прежде всего предполагает тщательный расчет общей поверхности охлаждения, так как именно от этого параметра во многом зависит протекание процесса теплоотдачи. Предполагается, что общая поверхность охлаждения радиатора должна обеспечить отвод теплоты от охлаждающей жидкости системы охлаждения на всех режимах работы дизеля, включая режимы перегрузки и резкие колебания температуры окружающей среды.
Следуя разработанным О. Кришером рекомендациям [1], можно проводить расчет охлаждающей поверхности радиатора (рис. 1), используя при этом лишь некоторые его намеченные конструктивные размеры, к которым
относятся параметры охлаждающих трубок и компоновка их в сердцевине.
90-100 °С
Рис. 1. Расчетная схема радиатора системы охлаждения
Преимущества данного метода заключаются в возможности проведения расчетов по определению требуемой поверхности охлаждения радиаторов с различными исполнениями сердцевины, например при коридорном или шахматном расположении охлаждающих трубок.
При проведении расчетов количества теплоты, передаваемой поверхностью радиатора системы охлаждения потокам воздуха независимо от пространственного положения трубок сердцевины радиатора, определим такой параметр, как пористость пучка охлаждающих трубок. При коридорном расположении трубок радиатора (рис. 2), перпендикулярных к направлению набегающего потока воздуха, а также при относительном продольном шаге трубок с >1
пористость трубного пучка рассчитывается по формуле [1]
71
¥к = 1-
4с'
(1)
где с - относительный продольный шаг трубок,
В
Продольный ивг трубок по фронту,
А сИ
I
^гкУ
+
I I
Скорость потока воздуха перед фронтом радиатора,
то
Рис. 2. Расчетная схема сердцевины радиатора системы охлаждения при коридорном расположении трубок
8д
Учитывая, что с = —-, где S,, - шаг трубок 2Я д
по фронту, м; Я - радиус закругления трубки, м, получаем следующее выражение:
пЯ
¥к=1-
71
2Я
-Л-
При шахматном расположении трубок сердцевины радиатора (рис. 3), когда относительный поперечный шаг Ь < 1, формула пористости трубного пучка будет иметь вид
где Ь - относительный поперечный шаг трубок
по глубине сердцевины радиатора, Ъ - ^.
$
Учитывая, что Ь =——, где 5", - шаг трубок 2Я 1
по глубине сердцевины, м, получаем следующее выражение:
71
¥ш=1" о „
2Я 2Я
= 1-
пЯ2
Диагональный шаг трубок,
Продольный шаг трубок по фронту,
Я, =сО
Поперечный шаг трубок по глубине,
Рис. 3. Расчетная схема сердцевины радиатора системы охлаждения при шахматном расположении трубок
В случае, когда трубки расположены в сердцевине радиатора шахматным порядком с узким сечением по диагонали с! < с, пористость пучка трубок следует рассчитывать по формуле
к
4ЙГ
где ё - относительный диагональный шаг трубок радиатора, с/ =
С учетом того, что с! — где .V,, - шаг
трубок по диагонали, м, получим следующее выражение:
V =1—=
2Я
Параметром, комплексно учитывающим траекторию движения воздушных масс в сердцевине радиатора, будем считать приведенную длину обтекания V. Длину обтекания воздушным потоком радиаторной трубки, оребренную
п
пластинами охлаждения (лентами), предлагается рассчитывать по формуле
2II
(2К)2 +
л2
0,565
V
/
где Ь1 - фронтальная ширина пластины охлаждения, м; Ь2 - длина пластины, м.
Г г:— Л2
Подставив обозначение = к, получим формулу
0,565Ц Р--Д
1/ = ^(2К)2+И2.
(2)
Переход от ламинарного пограничного слоя к турбулентному, как правило, происходит при больших скоростях воздушного потока. Учитывая скорость воздуха перед фронтом радиатора, следует также прибегнуть к средней эффективной скорости потока воздуха, которая определяется по формуле
V
¥
(3)
где у0 - скорость потока воздуха перед фронтом радиатора, м/с.
Расположение охлаждающих трубок в несколько рядов влечет за собой неравномерность теплоотдачи по глубине сердцевины радиатора. Полный тепловой поток, переданный поверхностью радиатора, состоящей из п рядов трубок, будет складываться из доли теплоты, переданной в первом и последнем по глубине рядах трубок:
а также из части теплоты, переданной в (п - 1) внутренних рядах:
п
б, =(»- 1Кбщ-^1(п)) = Х'7-/7' ^/ -*т(п)\
2
где сх1 - коэффициент теплоотдачи для отдельного ряда, отнесенный к разности температур на входе, Вт/(м2°С); ^ - площадь рабочей поверхности одного ряда трубок, м2; ^ - темпе-
ратура поверхности сердцевины радиатора, °С; ^1(1) - то же воздушного теплоносителя на
входе в сердцевину радиатора, °С; Ж - расход воздушного теплоносителя, перемещаемый крыльчаткой вентилятора, м3/ч; р(Г - плотность воздушного теплоносителя, кг/м3; сРж -теплоемкость воздушного теплоносителя, кДж/(кг °С); 1(п) - температура воздушного
теплоносителя на входе во второй ряд трубок сердцевины радиатора, °С; аобщ - коэффициент
теплоотдачи для всей сердцевины, отнесенный к разности температур на входе во второй ряд, Вт/(м2 оС); а; - то же для любого ряда внутри
сердцевины, отнесенный к разности температур на входе в этот ряд, Вт/(м2 °С); - темпе-
ратура воздушного теплоносителя на входе в п-й ряд сердцевины радиатора, °С.
Учитывая, что коэффициенты теплоотдачи а отнесенные к разности температур на входе в соответствующий ряд, будут равны, представим выражение для определения безразмерного коэффициента теплоотдачи для
ряда трубок внутри сердцевины [1]
№
1
Д.,
__о ?5Ре
Ообщ .. , ' Дэкв Т,
П- 1
Г
п—1
Ыи,
0 25Ре ^зкв
и
(4)
где Р^оэкв " число Пекле, РеДзкв = ш зкв; /Д)КВ -
УшД, а
эквивалентный диаметр, м.
При коридорном расположении трубок радиатора и относительном продольном ша-
1_Л
4с
Ь'; при
4 с
ге трубок (с>1) И =
71
шахматном расположении трубок сердцевины радиатора, когда относительный поперечный шаг Ь < 1, эквивалентный диаметр рассчи-
71
4 сЪ
и.
4с
тывается по формуле 1).ЖК = — 1 -
к
В случае шахматного расположения трубок сердцевины радиатора с узким сечением по
1
1
Л
4cd
L\ а -
диагонали (с! < с) /Д)К1,,,, = —| 1 -
к
коэффициент температуропроводности, м2/с; п - количество трубных рядов по глубине.
При проведении расчетов по определению теплоотдающей способности радиатора с помощью безразмерного коэффициента теплоотдачи Ыидобщ справедливо выражение, предложенное В. Нуссельтом:
Ми _ ^общ^экв
^ и£>общ ~~ ^ >
где X - коэффициент теплопроводности, Вт/(м-°С).
Выражение коэффициента ЫиЛобщ для всех
рядов трубок сердцевины радиатора по глубине представляется зависимостью
Nu
D общ
X
--0.25РС,
Nur
L'
-х
0,25Рег
Z),
L'
(5)
В данном выражении НиЛэкв представляет собой безразмерный коэффициент теплоотдачи для любого внутреннего ряда трубок сердцевины радиатора, который определяется по рабочей диаграмме № 1 [1] на основании соотноше-
D„,
ния
а также Pe,
D„,
^ , м- ^ ^^экв ^ '
Определение безразмерного коэффициента теплоотдачи Нидобщ, а далее и коэффициента
теплоотдачи аобщ дает возможность расчета
теплоотдачи поверхности радиатора, для чего используем формулу
^Qp аобп/рад(^ Vcp)'
(6)
где F - поверхность охлаждения радиатора, 2.
м
t
Wcp
средняя температура воздушного
теплоносителя, циркулирующего через сердцевину радиатора, °С.
Тогда требуемая поверхность охлаждения радиатора
F =
рад
да
а
общ^ F
(¡F Vcp)
(7)
Алгоритм расчета поверхности охлаждения радиатора по предлагаемой методике приведен на рис. 4.
На основании представленных теоретических положений проведем расчет требуемых поверхностей охлаждения радиаторов таких тракторов, как «Беларус-1221» с дизелем Д-260.2 и «Беларус-3022» с дизелем Detroit S40E, используя при этом лишь некоторые намеченные конструктивные размеры сердцевин радиаторов рассматриваемых тракторов «Беларус» (табл. 1).
Таблица 1
Некоторые параметры сердцевин радиаторов тракторов «Беларус»
Конструктивный размер сердцевины радиатора, м Тракторы «Беларус» типоразмерного ряда
«Беларус-1221» «Беларус-3022»
Дизели типоразмерного ряда тракторов «Беларус»
Д-260.2 Detroit S40E
ъ 0,01 0,01
к 0,0015 0,0015
¿1 0,007 0,007
¿2 0,092 0,138
Температура поверхности сердцевины радиатора °С 62 69
Средняя температура воздушного теплоносителя, циркулирующего через сердцевину радиатора ^^ °С 50 61
Скорость потока воздуха перед фронтом радиатора v0, м/с 10,3 18,3
Количество трубных рядов по глубине п 4 6
Отметим, что шаг трубок по фронту 5 для
всех сердцевин радиаторов типоразмерного ряда тракторов «Беларус» равен 0,01 м; шаг трубок по глубине сердцевины 5 - 0,023 м. Радиус закругления трубки - 0,0015 м. При коридорном и шахматном расположениях трубок в сердцевине радиатора, так как с = 3,33>1 и Ь = 7,67 > 1, получаем следующее значение пористости пучка охлаждающих трубок:
п
1
1
С =
2Я
$ - продольный шаг трубок по фронту
Определение пористости трубного пучка
*■-1-й II <5 ^ - поперечный шаг трубок по глубине
4с/ 2К
Определение длины обтекания Ь'
21
(2К)2+ 0,565
- диагональный шаг трубок
Ь1 - значение фронтальной ширины пластины охлаждения
Ь1 - значение длины пластины охлаждения
Определение средней эффективной скорости потока воздуха
V
Значение скорости потока воздуха перед фронтом
радиатора соответствует средней фактической скорости потока воздуха перед фронтом радиатора
г
0,565Д 1^-Я ' V ^ V х ) 2 = й
Определение эквивалентного диаметра ДЭ1
Определение числа Пекле РеД
Ре
_ уР»,
'¿)экв
с > 1
Ь < 1
В =4£
экв.к
К
1--5-
Ас V ^ /
Г
й < с
Б =4£
экв.ш
71
1-
71
4сЪ
Г
Определение числа Нуссельта № Дэкв по рабочей
в
экв.ш
71
к
4сй?
г
диаграмме № 1 О. Кришера в зависимости от Ре
Д.
£)экв
Определение безразмерного коэффициента теплоотдачи №
Добщ
N11 =—0 25Ре зкв
Добщ ' Дэкв
и
1- 1 —
N11 г
0 25Ре ^>зкв
Определение коэффициента теплоотдачи по формуле Ии^ - ————
«--шА X
г _ ^р
Определение поверхности охлаждения радиатора -► рад ге г/ _/ ) "•общ^ '(Гер-'
V
о
V =
а
Рис. 4. Математическая модель расчета поверхности охлаждения радиатора
к(ш) 2-0,01
Дальнейшие результаты расчета сводим в табл. 2.
Таблица 2
Расчетные значения параметров
Расчетный параметр Единица измерения Тракторы «Беларус» типоразмерного ряда
«Беларус-1221» «Беларус-3022»
h - 1,673 • 10-7 3,712 • 10-7
L м 0,00471 0,00471
Vw м/с 13,473 23,937
Аэкв м 0,0153 0,0153
РеАэкв - 7921,84 14074,7
-РеАэкв ДкД - 25699,5 45660,3
№дэкв - 120 175
- 116,68 168,43
С^общ Вт/(м2-°С) 222,86 321,71
Приняв, что количество теплоты, подводимой охлаждающей жидкостью к поверхностям сердцевины радиатора, равно количеству теплоты, отводимой от радиатора потоком воздушного теплоносителя, проведем расчет требуемых поверхностей охлаждения радиаторов рассматриваемых тракторов.
Для определения количества теплоты, подводимой охлаждающей жидкостью к сердцевине радиатора, воспользуемся данными теплового баланса. Распределение теплоты по составляющим теплового баланса дизелей тракторов «Беларус», полученное по результатам экспериментальных исследований в УП «ММЗ», приведено в табл. 3. Так как данные по тепловому балансу для дизеля Detroit S40E отсутствуют, условимся считать их равными тепловыделению дизеля Д-260.16.
Таблица 3
Составляющие теплового баланса
Количество теплоты, воспринятое охлаждающей жидкостью, Вт Модель дизеля
Д-260.2С2 Д-260.1С2 Д-260.4С2 Д-260.16
Q 44000,0 54267,2 67916,1 86272,0
Используя данные табл. 2 и 3, по формуле (7) определяем требуемые поверхности охла-
ждения радиаторов типоразмерного ряда тракторов «Беларус»:
• радиатор «Беларус-1221»
F 44000,0 2
Кап = -- = 16,5 М
рад 222,86-(62-50)
• радиатор «Беларус-3022»
г 86272,0 „ - 2
=-2-= 33,5 м .
рд 321,71 (69-61)
Выполним по предлагаемой методике расчет поверхности охлаждения радиаторов типо-размерного ряда тракторов «Беларус».
Разработка охлаждающих систем тракторной техники, обеспечивающей рентабельность у производителя, невозможна без системных критериев, в основе которых лежат унификация и типизация. Цель унификации - сокращение сроков и затрат на проектирование изделий, повышение производительности и снижение затрат производства, повышение качества изделий за счет специализации и комплексной механизации производства, улучшения эксплуатации (обслуживание, ремонт) за счет взаимозаменяемости. Типизация является одним из направлений проектирования, позволяющим едиными составляющими создавать изделия с различными техническими параметрами.
При разработке охлаждающих систем тракторов нами рассматривалась унификация по «горизонтали», т. е. внутри тракторов данного тягового класса, и по «вертикали», т. е. между тракторами разных тяговых классов.
Представляется, что тракторы «Беларус» при разработке охлаждающих систем можно подразделить на четыре группы по мощности [2, 3]:
А: 43,5-66 кВт;
В: 79-88 кВт;
С: 103-132 кВт;
Б: 186-221 кВт.
Тракторы каждого диапазона мощности, как показывают расчеты, могут иметь унифицированные охлаждающие системы, включающие в себя один типоразмер водяного радиатора. Типораз-мерный ряд радиаторов тракторов «Беларус» приведен в табл. 4.
Таблица 4
Типоразмерный ряд радиаторов тракторов «Беларус»
Наименование составляющих Основные рабочие параметры по рядам
А В С Б
Радиатор (поверхность охлаждения) ^ м2 15,0
23,0
25,0
36,0
При разработке конструкции радиаторов возможны: более глубокая унификация по элементам деталей, две типовые конструкции водяных радиаторов - для четырех- и шестицилиндровых дизелей.
Надежность функционирования охлаждающих систем определяется их способностью обеспечивать заданную температуру охлаждающей жидкости и смазочного масла на всех скоростных и нагрузочных режимах работы дизеля, возможных колебаниях климатических параметров. При этом энергетические затраты дизеля на привод водяного и масляного насосов, вентилятора должны быть наименьшими.
Для оценки совершенства охлаждающей системы используется энергетический метод, характеризующийся соотношением переданного количества тепловой энергии к энергии, затраченной на передачу этой теплоты. При этом отводимый тепловой поток и затраченную энергию относят к единице поверхности охлаждения.
Исследования показали, что для оценки охлаждающей системы более рационально рассматривать радиатор комплексно с заданными теплоносителями. Предлагается эффективность системы оценивать по тепловоспринимающей и теплорассеивающей способностям радиатора, потоков жидкости и воздуха, проходящих через радиатор.
Выбор поверхности радиатора проводится по принятым расходам охлаждающей жидкости и потока воздуха. Основным условием при выборе поверхности охлаждения является обеспечение заданного температурного режима дизеля. Эффективность охлаждающей системы, т. е. выбранной поверхности радиатора и принятых расходов теплоносителей, предлагается оценивать отношением количества теплоты,
которое отводится или необходимо отвести от охлаждающей жидкости, к тепловоспринимаю-щим способностям радиатора или потока воздуха.
Системы охлаждения рекомендуется оценивать коэффициентом использования поверхности радиатора, равным отношению отводимого теплового потока для обеспечения заданного температурного режима к максимально возможной теплорассеивающей способности этого радиатора:
Предлагаемый метод позволяет качественно оценивать тепловую эффективность охлаждающей системы.
Для оценки предлагаемого типоразмерного ряда радиаторов воспользуемся количественной оценкой параметров принятых радиаторов.
На рис. 5 и 6 представлены графики типо-размерных рядов радиаторов в функции охлаждающей поверхности от мощности дизелей. Коэффициент использования поверхности радиаторов Ер с перспективными дизелями по группам равен:
А: 0,847-1,04;
В: 0,813-0,978;
С: 0,678-0,972;
Б: 0,703-0,97.
Рис. 5. Типоразмерный ряд радиаторов (типоразмер 1) по поверхности охлаждения
Рис. 6. Типоразмерный ряд радиаторов (типоразмер 2) по поверхности охлаждения (точками отмечены дизели нового типоразмерного ряда)
Расчеты показывают, что дизели группы А Д-243 и Д-245.5 не имеют запаса по охлаждающей поверхности радиаторов. Численное значение недостающей поверхности к принятой охлаждающей поверхности радиатора типоразмерного ряда составляет 2 и 4 %. Дизель группы Б Д-260.10 имеет запас поверхности охлаждения радиатора всего 3 %. При комплектовании охлаждающих систем тракторов с дизелями Д-245.5 и Д-260.10 рекомендуется использовать радиаторы этого типоразмерного ряда, но охлаждающие поверхности выполнять из более эффективных по теплопередаче материалов, например меди.
В Ы В О Д Ы
Разработан метод расчета охлаждающей поверхности радиатора, в основе которого при определении коэффициента теплопередачи приняты параметры элементов и компоновки сердцевины. Движение воздушных масс в сердце-
вине учитывается приведенной длиной обтекания охлаждающих трубок потоком воздуха. Существенным достоинством метода является учет неравномерности теплопередачи по глубине радиатора.
Предложен типоразмерный ряд радиаторов для мощностного ряда тракторов «Беларус». Оценка тепловой эффективности радиаторов предлагаемого типоразмерного ряда проведена по тепловоспринимающей и теплорассеиваю-щей способностям охлаждающих поверхностей.
Л И Т Е Р А Т У Р А
1. Конвективный тепло- и массоперенос. Единое описание для течения в каналах и внешнего обтекания тел любой формы и расположения: пер. с нем. / В. Каст [и др.]. - М.: Энергия, 1980. - 49 с.
2. Тарасенко, В. Е. Типоразмерный ряд систем охлаждения тракторов «Беларус» / В. Е. Тарасенко, А. И. Якубович // Prospective Technics and Technologies, 2006: ма-терiали Il-i Мiжнар. наук.-практ. конф. студ. i молодих учених, Mykolaiv, September 14-16, 2006 / Mykolaiv State Agrarian University; вщповщ. за випуск С. I. Пастушен-ко. - Mykolaiv, 2006. - С. 123-126.
3. Тарасенко, В. Е. Унификация систем охлаждения мощностного ряда тракторов «Беларус» / В. Е. Тарасенко, А. И. Якубович // Развитие приграничных регионов Беларуси и России на современном этапе: проблемы и перспективы: материалы науч.-практ. конф., Могилев, 12-13 окт. 2006 г. / М-во образования Респ. Беларусь, М-во образования и науки Рос. Федерации, Федеральное агентство по образованию, Бел.-Рос. ун-т; редкол. И. С. Сазонов (гл. ред.) [и др.]. - Могилев: Бел.-Рос. ун-т, 2006. -С. 22.
Поступила 24.04.2008