ЭНЕРГЕТИЧЕСКОЕ И ТРАНСПОРТНОЕ МАШИНОСТРОЕНИЕ
УДК 621.833
А. Л. Носко, А. П. Носко
ИССЛЕДОВАНИЕ ОХЛАЖДЕНИЯ ТОРМОЗНЫХ УСТРОЙСТВ ПОДЪЕМНО-ТРАНСПОРТНЫХ МАШИН
Рассмотрены методы определения коэффициента теплоотдачи тел из фрикционных материалов. Приведены примеры расчета коэффициента теплоотдачи тормозных накладок, выполненных из полимерных и спеченных порошковых фрикционных материалов, а также шкивов/дисков из чугуна.
При выборе тормозного устройства основными характеристиками являются долговечность и работоспособность материалов фрикционной пары тормоза. Долговечность пары трения, в первую очередь, определяется износом тормозной накладки, а работоспособность — стабильностью тормозного момента. При работе тормозного устройства фрикционная пара нагревается, что приводит к изменению указанных характеристик — интенсивность изнашивания возрастает, а стабильность тормозного момента снижается. Таким образом, при оценке долговечности и работоспособности тормозного устройства необходимо исследовать температурный режим работы фрикционной пары.
Повторно-кратковременный режим работа (а именно, в этом режиме работают тормоза большинства подъемно-транспортных машин (ПТМ), характеризуется наличием периодически повторяющихся торможений (нагрев) и пауз (охлаждение). Как показывает практика, после ряда последовательных торможений наступает установившийся тепловой режим [1], при котором количество теплоты, образующееся при трении, равно количеству теплоты, отводимому в окру ау среду. Коэффициент теплоотдачи характеризует интенсивность отвода теплот от элемента пары трения в окру ау среду (в наем случае в воздух) и является одним из параметров, определяющим нагрев материалов фрикционной пары тормоза.
Методы определения коэффициента теплоотдачи. Величина коэффициента теплоотдачи материала в окру ау среду зависит от относительной скорости среды и обработки поверхности материала.
Определение коэффициента теплоотдачи тормозных накладок. рактически все тормозные накладки Т име т форму прямоуголь-
ного параллелепипеда (пластины) или могут быть представлены в виде пластины [1,2].
В материалах, которые имеют низкую теплопроводность, например
в фрикционных полимерных материалах (ФПМ), у которых Л = 0,4 ... Вт
. 0,8
м2
◦с
распространение теплота и, соответственно, выравнива-
ние температур происходит очень медленно. Это значительно усложняет методику расчета коэффициента теплоотдачи, так как требует учета неравномерного поля температур.
Примем следующие допущения. 1. Накладка представляет собой прямоугольный параллелепипед длиной I, шириной Ь и толщиной к. 2. Материал накладки — однородный и изотропный. 3. Теплоотдача происходит через поверхность трения и боковые поверхности. При этом поверхность основания накладки считается теплоизолированной. 4. В начальный момент времени температура накладки во всех точках одинакова и равна Т0.
Распределение теплоты в накладке происходит согласно уравнению теплопроводности, которое в прямоугольной системе координат имеет вид
дТ
дт
= a
д 2T д 2T д 2T + +
(1)
дх2 ду2 дг2
где Т = Т(х, у, г, т) — температура; т — время, (х,у, г) — координаты точки накладки; а — температуропроводность материала накладки.
При температуре окружающей среды Те = 0 граничные условия принима т вид
дТ дТ
Л— = вТ, х = 0; -Л— = вТ,х = I; дх дх
дТ дТ
Л— = вТ, у = 0; —Л— = вТ,у = Ь; ду ду
дТ дТ
Л— = 0, г = 0; — Л— = вТ, г = к. дг дг
где Л — коэффициент теплопроводности; в — коэффициент теплоотдачи.
Согласно принятому четвертому допущению начальное условие можно записать в виде
Т(х, у, г, 0) = То.
В работах [3, 4] рассмотрены некоторые методы решения и примеры краевых задач остывания нагретых тел. Уравнение (1) решается методом Фурье (методом разделения переменных). Представим частное
решение уравнения (1) в виде произведения
T(x,y,z,T)= X(x) Y(y) Z(z) U(t).
Подставляя выражение (2) в уравнение (1), после разделения переменных получим
1 U _ Z" _ я
au = x + y + т =
где d = const.
Отсюда находим уравнение для функции U (т):
U' (т) + daU (т) = 0,
общее решение которого имеет следующий вид:
(3)
U (т) = C e
—dar
где С — произвольная постоянная.
С учетом граничных условий имеем три задачи Штурма-Лиувилля:
X'' + X = 0,
Ü
X'(0) - Т-X(0) = 0, А
а
X'(/) + Т-X (l) = 0, А
f Y'' + V2 Y = 0,
Y'(0) - TY(0) = 0, А
a
Y'(b) + t-Y (b) = 0, А
Z'' + w2 Z = 0,
Z'(0) = 0, fj
Z'(h) + ^Z (h) = 0, A
которые имшт следящие собственные значения:
П N А / Т2 ctg(/^n) = — ^ 1- —
,, ч А / Т2
ctg(bvm) = 2Т Vm ^ - А2"Т
ctg(h^fe) = ,
(4)
и соответствующие им собственные функции
в
Хп(х) = соэ^п х — вп), вп = аг^ ■
ут(у) = соэ(^т у — 7т), 7т = аг^ —( )
Л ^т
к ^^(г) = соэ^ г). В итоге для уравнения (3) находим собственные значения
Аптк = ^П + + Шк (6)
и собственные функции
ипшк(х, у, г) = соэ(^п х — вп) соэ(^т у — 7т) соэ^ г).
Общее решение уравнения (1) запишем в виде тройного тригонометрического ряда
Т (х,у,г,т) =
оо оо оо
= Сптк е 4аТ соэ(^п х — вп) соэ(^т У — 7т) ^(^к г).
п=0 т=0 к=0
(7)
В соответствии с теоремой Стеклова [4] коэффициент Сптк можно определить по следующей формуле:
I ь н
Т0 соэ(^п х — вп) соэ(^т у — 7т) соэ(шк г) ¿х dy Аг
с = 000_
Сптк = 1 ь н ,
соэ2(^п х — вп) соэ2(^т у — 7т) соэ2(шк г) Ах Ау Аг 000
после соответствующих преобразований окончательно получим
С =32 т ^(0,5 Уп1) соэ(0,5 — вп) х
птк 0 + Й1п(^п1) соэ(^п1 — 2 вп)
х 81п(0,5УтЬ) соэ(0,5УтЬ — 7т) 81п(^кк) (§)
^тЬ + Э1п(^тЬ) соэ(^тЬ — 27т) шкк + 0,5 эт( 2 шкк)
Если температура окружающей среды Те не равна нулю, то необходима замена переменной Т = Т — Те. В дальнейшем под Т(х, у, г, т)
будем понимать решение уравнения (1) с произвольной начальной температурой.
пределить коэффициент теплоотдачи мо но по экспериментальным данным остывания накладки: по координатам расположения термопары (х0, у0, г0); вектору температуры Т*[г] и вектору времени т*[г], г = 1, 2,... ,п, где п — количество экспериментальных значений; по теплофизическим свойствам материала накладки а, Л; по размерам накладки I, Ь, к и температуре окружающей среда Те и начальной температуре накладки Т0.
Из полученного общего решения (7) видно, что выразить коэффициент теплоотдачи в явном виде не представляется возможным. Поэтому предлагается определять коэффициент теплоотдачи путем его варьирования при неизменных значениях остальных параметров. Изменяя значение коэффициента теплоотдачи в, получаем различные значения температуры Т(х0,у0, г0,т* [г]). Сравнивая расчетные значения Т(х0,у0, г0,т* [г]) с экспериментальными значениями Т*[г], определяем в*, при котором эти значения температур будут наиболее "близки". Полученное значение в* считаем искомым. В качестве критерия "близости" температур возьмем среднее относительное отклонение расчетных значений от экспериментальных значений температур а (в), определяемое по формуле
Таким образом, задача определения коэффициента теплоотдачи сводится к задаче поиска минимума функции а (в*) ^ min.
Определение коэффициента теплоотдачи фрикционных материалов с высокой теплопроводностью .Для материалов с высокой тепло-
Вт
проводностью (для стали и чугуна Л = 40 ... 80 —◦—, для фрикцион-
м2 C
Вт
ного спеченного порошкового материала (ФСПМ) Л = 30 ... 60 ——)
методика расчета коэффициента теплоотдачи значительно упрощается, так как распространение теплота происходит настолько быстро, что температурное поле материала можно считать равномерным. Примем следующие допущения: материал тела — однородный и изотропный и температура тела во всех точках одинакова и равна Т.
При остывании тело отдает теплоту в окружающую среду. Этот процесс можно описать с помощью уравнения Ньютона
1 п
а(в) = П £
i= 1
|T(x0,y0,Zo,T*[i]) - T*[i]|
T *[i]
(9)
м2°С
dß = -es (T - Te),
(10)
где Q — внутренняя энергия тела; в — коэффициент теплоотдачи; Б — площадь поверхности охлаждения; Те — температура окружающей среды.
С другой стороны, изменение энергии тела можно представить как
dQ = cm dT,
(ii)
где с — удельная теплоемкость; т — масса тела.
Подставляя выражение (11) в уравнение (10), получим дифференциальное уравнение
ЯТ в Б
(12)
T-Te
ÖS т =--ат,
cm
решая которое при начальном условии Т(0) = То относительно в, получим
ст Т-Т
(13)
д cm To - Te
' = ST ln т^
Экспериментальные данные по охлаждению тела, коэффициент теплоотдачи которого необходимо определить, должны содержать следующие данные: вектор температуры Т*[г] и вектор времени т*[г], г = 0, 1, ..., п, где п + 1 — количество экспериментальных значений; теплоемкость материала с; массу тела т и площадь охлаждения Б; температуру окружающей среды Те и начальную температуру То.
По экспериментальным значениям времени т*[г] и температуры Т*[г] строится кривая охлаждения (рис. 1). При этом кривая разбивается на п интервалов, для каждого из которых по формуле (13) определяется соответствующее значение коэффициента теплоотдачи [г]. В качестве искомого в* возьмем среднее значение для всего интервала времени:
1
г = (14)
Определение коэффициента теплоотдачи вращающихся тормозных шкивов/дисков. Если тело неподвижно, то вся его поверхность имеет одинаковый коэффициент теплоотдачи в(т,ж,у,г) = в.
При вращении тела количество теплоты, отдаваемое в окружающую среду, увеличивается.
Рис. 1. Температурная кривая охлаждения для определения коэффициента теплоотдачи в окружающую среду
При этом различные участки поверхности тела имеют разный коэффициент теплоотдачи в = в (т,х, у, г).
Для описания процесса введем коэффициент теплоотдачи вращаю -щегося тела
в,(т) = 1 I в(т,х,у,г) 4Б. (15)
я
При этом формула (10) примет вид
^ = -вЕ 5 (Т - Те). ат
В работе [5] для определения коэффициента теплоотдачи в тормозных дисков из чугуна приводится эмпирическая формула:
в = 7,14 V0'78
где V — скорость движения среда относительно охлаждаемой поверхности.
Выразив скорость V через частоту вращения тела п, найдем коэффициент теплоотдачи участка поверхности, который находится на расстоянии Я от оси вращения:
в = 1,23 Я0'78 п0'78. (16)
Подставив выражение (16) в формулу (15), получим уравнение для определения коэффициента теплоотдачи тела:
1 23 п0'78 г
вЕ(п) = У Я0'78аБ. (17)
я
римем, что
1 23 Г Я0'78а,Б = д, (18)
1'23 I р0,78,
s J
S
причем коэффициент д полностью определяется геометрией тела. Теперь формулу (17) запишем в следующем виде:
вЕ(п)= дп0'78. (19)
Расчетно-экспериментальные исследования. В работах [1, 2] приведена информация о накладках, впускаемых отечественной промышленностью. Режимы и параметры испытаний тормозов ПТМ, результаты которых использованы в настоящей статье, подробно изложены в работах [5, 6].
Теплофизические свойства некоторых материалов пар трения тормозов ПТМ
Свойство Единица измерения Материал
ФПМ 145-40 ФМК-845 Серый чугун
Теплопроводность, Л Вт/(м • °С) 0,64 38 43... 62
Теплоемкость, с Дж/(кг • °С) 1100 630 500... 550
Плотность, р кг/м3 2500 6000 7150
Температуропроводность, а 10-6(м2/с) 0,23 10 11...17
Определение коэффициента теплоотдачи тормозной накладки из ФПМ. В качестве ФПМ возьмем материал шифра 145-40, теплофизические свойства которого представлены в таблице [1]. С помощью формулы (9) по экспериментальным данным были определены зависимости а = а (в). На рис. 2 приведены две огибающие (сверху и снизу) экспериментальные кривые. При а < 3% коэффициент теплоотдачи материала Вт
в* = 51 ... 59——.
м2-°С
Используя формулу (7), определим, как распределяется температу-
Вт
ра в накладке из ФПМ в процессе остывания, приняв в* = 54
м
2. °0
Результаты расчета приведены на рис. 3.
На рис. 4 показано распределение температуры по толщине накладки из ФПМ при остывании через определенные промежутки времени.
а
о
Щ И ФПМ ш ш
¥ 1 1
1 Ы ш V/
1 щ
50
51
54
56 в, Вт К м С)
Рис. 2. Диапазон вариации коэффициента теплоотдачи в* накладки из ФПМ 145-40 в зависимости от среднего относительного отклонения и (в) расчетных значений температуры от экспериментальных данных
Т,°С 80 60 40
го о
Vs» w \ ■ ^ 'S.
Xi
Те ^—-
200 400
600
800 t,c
Рис. 3. Графики изменения температуры накладки из ФПМ 145-40 при остывании
а
1,6 и 0,6 ОА
h,MM
Рис. 4. Распределение температуры по толщине накладки из ФПМ 145-40 (х = 0,51, у = 0,56) при остывании; 1,...,6 — т = 0; 100; 200; 300; 400; 500 с соответственно
0
ФС! IM к
\
V/J
у
в*
46
46
S0 в,Вт/(м С)
Рис. 5. Диапазон вариации коэффициента теплоотдачи 6* накладки из ФСПМ типа ФМК-845 в зависимости от среднего относительного отклонения и(&) расчетных значений температуры от экспериментальных дан-нх
Как видно, температура накладки на глубине значительно выше, чем на поверхности, что объясняется низкой теплопроводностью материала.
Определение коэффициента теплоотдачи тормозной накладки из ФСПМ. В качестве примера определим коэффициент теплоотдачи ФСПМ типа ФМК-845, теплофизические свойства которого приведены в таблице [7]. Так как рассматриваемый материал имеет высокую теплопроводность, то можно воспользоваться формулами (9) и (14). Для каждого эксперимента были определены: зависимость а = а (9) — по формуле (9), соответствующий коэффициент теплоотдачи 9* — по формуле (14). Результаты расчетов приведены на рис. 5, где показаны две огибающие (сверху и снизу) кривых а = а (9), а также указан разброс значений коэффициента теплоотдачи 9*. Объединив полученные
Вт
результаты при а < 1%, получим 9* = 46 ... 53 —5—.
м2 С
Определение коэффициента теплоотдачи тормозного шкива. При расчете коэффициента теплоотдачи кива считаем, что теплообмен осуществляется только через поверхности соприкосновения с воздуш-ной средой. Количество теплоты, проходящей через поверхности шкива, сопрягаемые с поверхностями других деталей, пренебрежимо мало, так как между этими поверхностями существует воздушный технологический зазор, который препятствует перетеканию теплоты.
Для шкивов различных диаметров О по формуле (18) был рассчитан коэффициент д (О), значения которого приведены ниже.
Диаметр шкива О, м Коэффициент д.....
0,2 0,3 0,4 0,5 0,6 0,18 0,25 0,31 0,36 0,42
Аппроксимировав полученные значения линейной зависимостью с относительной погрешностью в узлах менее 5 %, получим
©Е(Д п) = (0,07 + 0,59 О) п
0,78
(20)
Результаты расчетов по формуле (20) приведены на рис. 6. Как видно, с увеличением диаметра коэффициент теплоотдачи шкива увеличивается.
Определение коэффициента теплоотдачи тормозного диска. Примем следующие допущения: диск — идеальный цилиндр диаметром О; теплообмен осуществляется через всю поверхность, и диаметр диска намного больше его толщины.
Если при указанных допущениях по формуле (17) рассчитать коэффициент теплоотдачи диска для дисково-колодочных тормозов, то получим
6Е(Д п)=0,52 О0,78 п0,78. (21)
В работе [5] приведена экспериментальная зависимость коэффициента теплоотдачи тормозного диска диаметром 0,2 м, изготовленного из чугуна СЧ15-32, от частоты вращения, которая достаточно точно аппроксимируется формулой
6|(0,2,
п)
в^, Вт/(м2°С) 130
= 27 + 0,0165 п.
На рис. 7 приведены графики теоретической (21) и экспериментальной (22) зависимостей коэффициента теплоотдачи диска от частоты вращения. Результаты, полученные при использовании формулы (21), хорошо согласуются с экспериментальными данными для больших скоростей вращения диска. По формуле (21) для дисков разных диаметров построены графики зависимости коэффициента теплоотдачи диска от частоты вращения, приведенные на рис. 8.
Выводы. При наличии экспериментальных данных формулы (9) и (14) позволяют определить коэффициент теплоотдачи материалов тормозных накладок из ФПМ и ФСПМ,
(22)
1300 п, мин1
Рис. 6. Зависимость коэффициента теплоотдачи ©е тормозных шкивов от частоты вращения; 1,...,5 — О = 0,2; 0,3; 0,4; 0,5; 0,6 м соответственно
40 30
го ю
°500 700 900 1100 1300л,мин1
Рис. 7. С равнение теоретических и экспериментальных в| значений коэффициентов теплоотдачи тормозных дисков; В = 0,2 м; материал СЧ15-32
90 70 50 30
10 , 500 700 900 1100 1300 п, мин1
Рис. 8. Зависимость коэффициента теплоотдачи в^ тормозных дисков от частоты вращения; 1,...,5 — D = 0,2; 0,3; 0,4; 0,5; 0,6 м соответственно
elJ
Ol
а формулы (20) и (21) — коэффициент теплоотдачи тормозных шкивов/дисков, выполненных из чугуна (стали). Полученные результаты могут быть использованы в контактных температурных задачах при выборе материалов пар трения и оценке их характеристик как на стадии проектирования, так и эксплуатации колодочных, дисково-колодочных и дисковых тормозов ПТМ.
Учитывая сложность определения и отсутствие в литературе опытных данных коэффициентов теплоотдачи тормозных накладок из ФПМ и ФСПМ, точность определения коэффициента теплоотдачи для указанных материалов с использованием полученных зависимостей можно считать вполне приемлемой.
СПИСОК ЛИТЕРАТУРЫ
1. Полимеры в узлах трения машин и приборов: Справочник / А.В. Чичинадзе, А.Л. Левин, М.М. Бородулин, Е.В. Зиновьев; Под общ. ред. А.В. Чичинадзе. -2-е изд., перераб. и доп. - М.: Машиностроение, 1988. - 328 с.
2. Справочник по триботехнике. В 3-х т. Т. 3. Триботехника антифрикционных, фрикционных и сцепных устройств. Методы и средства триботехнических испытаний / Под общ. ред. М. Хебды, А.В. Чичинадзе. - М.: Машиностроение, 1992.-730 с.
3. Лыков А. В. Теория теплопроводности. - М.: Высшая школа, 1967. - 598 с.
4. Мартинсон Л. К., Малов Ю. И. Дифференциальные уравнения математической физики: Учеб. для вузов. 2-е изд. /Под ред. В.С. Зарубина, А.П. Крищенко. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2002. - 368 с.
5. Ромашко А. М. Исследование нагрева дисковых колодочных тормозов подъемно-транспортных машин. - М., 1979. - 254 с.
6. Носко А. Л. Разработка методики расчета нагрева дисково-колодочных тормозов подъемно-транспортных машин с учетом термического сопротивления контакта. - М., 1985.- 199 с.
7. Федорченко И. М., Ровинский Д. Я., Шведков Е. Л. Исследование материалов для тормозных и передаточных устройств. - Киев: Наукова думка, 1976. - 199 с.
Статья поступила в редакцию 25.02.2004
Андрей Леонидович Носко родился в 1957 г., окончил Луганский машиностроительный институт в 1979 г., канд. техн. наук, доцент кафедры "Подъемно-транспортные системы" МГТУ им. Н.Э. Баумана. Автор более 65 научных работ в области трения и износа фрикционнных пар тормозных устройств.
Andrey Leonidovich Nosko wos born in 1957, graduated from Lugansk Meckanical Engineering Institute in 1979. PhD (Eng.), assoc. professor of "Handling Transport Systems" department of the Bauman Moscow State University. Author of more than 65 publications in the field of friction and wear of brakes.
Алексей Павлович Носко родился в 1983 г., студент кафедры "Системы автоматизированного проектирования" МГТУ им. Н.Э. Баумана.
Alexey Pavlovich Nosko was born in 1983, student of the "Computer Aided Design" department of the Bauman Moscow State University.
ЖУ^МЛ "ВЕСТНИК МОСКОВСКОГО ГОСУДАРСТВЕННОГО ТЕХНИЧЕСКОГО УНИВЕРСИТЕТА имени Н.Э. БАУМАНА"
В журнале публикуются наиболее значимые результаты фундаментальных и прикладных исследований и совместных разработок, выполненных в МГТУ имени Н.Э. Баумана и других научных и промышленных организациях.
Журнал издается в трех сериях: "Приборостроение", "Машиностроение", "Естественные науки" — с периодичностью 12 номеров в год.
Подписка по каталогу "Газеты, журналы" агентства "Роспечать"
Индекс Наименование серии Объем выпуска Подписная цена (руб.)
Полугодие 3 мес. 6 мес.
72781 "Машиностроение" 2 150 300
72783 "Приборостроение" 2 150 300
79982 "Естественные науки" 2 150 300
Подписывайтесь и публикуйтесь!
Адрес редакции журнала "Вестник МГТУ имени Н.Э. Баумана": 105005, Москва,
ул. 2-я Бауманская, д. 5.
Тел.: (095) 263-62-60.
Факс: (095) 265-42-98; 263-67-07.
E-mail: [email protected], [email protected]
4 ^ у *»
ль*!