ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 21.10.2024
Просмотров: 72
Скачиваний: 0
и удовлетворением граничных условий получаем систему алгебра ических уравнений
|
-С ,„[ц А |
|
|
|
( £^ -)]+ Я ,4 С а„ = 0, |
|
||||
п Т |
( Ц/1^1 \ ! г* |
7» |
( № n R 1 \ |
| |
f-1 |
п |
Т ( ^ n R l |
\ |
S^n ^П |
А |
6^ У0 \—Г ~ )+^2п^о { —r - ) + C3n— t i n h { —— |
) |
sin p„ ~ U’ |
||||||||
С,А (тг-)-я .А (-й#Ц-с. |
к в |
sin %п |
|
|
|
|||||
|
|
|
|
|
4п |
к р |
sin р„ |
|
|
|
|
|
|
С$п = |
DW*An |
|
|
|
(3.16) |
||
Здесь |
/ 0, /ц к0, к1 — функции Бесселя и Макдональда мнимого |
|||||||||
аргумента. |
|
|
|
sin (1п |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
•4 |
‘ + т ^ ) |
|
|
|
|
||
W* — усредненное |
по |
объему |
значение |
функции |
источника |
W. |
||||
Величины р„, Хп и А п определяют по таблицам |
[37]. |
|
||||||||
Пример 1. Рассмотрим теплообразование в полом цилиндри |
||||||||||
ческом |
амортизаторе |
размером |
2h — 0,175 м, |
|
2R ± = 0,07 м |
|||||
и 27?2 = 0,2 м из резины на основе НК твердостью 60—63 |
по |
|||||||||
IM-2. |
Детали такого |
типа получили |
широкое |
распространение |
в качестве элементов упругой подвески тяжелых вертикальных виброконвейеров. Температуру в амортизаторе при его гармони ческом возбуждении замеряли при следующих параметрах нагру жения: Д = 0,012 м; 6 = 0,0096 м и <а = 81,6 1/с. Реологические параметры резины по табл. 2.2 имели следующие значения: а =
— 0,59; р = 0,52; % = 0,244; G0 = 2,2 МН/м2. Используемые теплофизические характеристики определяли расчетным путем
(Н^ = |
180 1/м, Н г = 40 1/м) и заимствовали из литературы [37] |
|
(кР = |
0,293 Вт/м•°С, кв = |
0,242 Вт/м-°С). |
Расчетные значения 0Р, |
найденные для исследуемого аморти |
затора по формуле (3.13) с помощью ЭЦВМ, и экспериментальные данные, полученные статистической обработкой на основе распре деления Стьюдента с доверительной вероятностью 0,95, предста влены на рис. 54.
Несовпадение расчетных данных с экспериментальными в об ласти температур 0 > 80° С можно, по-видимому, объяснить тем, что при расчете не учитывали зависимость величины я[з от темпе ратуры. Распределение температур по высоте амортизатора при радиусе г = 0,064, полученное расчетным путем, удовлетвори тельно совпадает с экспериментальным.
С л у ч а й 2. Определим температурное поле в рассматри ваемом полом цилиндрическом амортизаторе (см. рис. 53), реологи
77
ческие характеристики которого изменялись с увеличением тем пературы. В первом приближении можно ограничиться лишь температурной зависимостью коэффициента поглощения энергии, экспериментальный график которой представлен на рис. 55 штри ховой линией [47]. Приближенно его можно представить в виде кусочно-линейной функции
+ 6р)(9р — 6р)> (3.17)
где ф0 — значение коэффициента поглощения энергии при 20° С; Ct — коэффициент пропорциональности; 0|> — значение темпе ратуры, начиная с которой зависимость ф (0р) становится суще ственно нелинейной; Н (0р, 0J) — функция Хэвисайда,
при 0р < 0р;
Я(0Р,
при 0р ss 0р.
Как и ранее, задача для вычисления температурного поля сво дится к решению уравнений теплопроводности, где функция источника тепла W определяется с учетом зависимости коэффи циента ф от температуры (3.17) и для вычисления приращения тем пературы 0Ввоздуха, заключенного в полости амортизатора (3.11).
При этом граничные условия, отражающие теплообмен с окру жающей средой, условия равенства температур и тепловых потоков на внутренней поверхности цилиндра (в интегральном смысле) остаются такими же (см. 3.12а—д).
Рис. 54. |
Распределение температуры |
Рис. 55. |
Температурная зависимость |
||
( |
в объеме образца |
коэффициента поглощения энергии |
|||
расчет; • — эксперимент) |
ПРИ со = |
80 |
1 / с ; -------------- |
кусочно |
|
|
|
линейная апроксимация; |
|||
|
|
---------------- |
|
эксперимент |
|
Будем искать решение поставленной |
задачи |
приближенно |
|||
в виде |
|
|
|
|
|
|
0р — 0р-[-0р, |
0в —QjB+ 0h> |
(3.18) |
||
где 0р и 0р удовлетворяют уравнениям |
|
|
|
||
|
у20оз + ^ | ^ 1 /(г> г) = 0) |
у20& = о |
(3.19) |
78
и граничным условиям (3.12, а—д). Величины 0), и 0д определялись уравнениями
^ + |
6*)(02, - 0*)/(г, z) —0; |
(3.20) |
|
V20b = O, |
|
причем функция источника тепла определена в области, для кото рой выполняется условие 0J, ^ 0J. Эту область приближенно (возможен и графический способ) можно представить как полый цилиндр с внутренним радиусом R x, наружным радиусом RI и высотой 2h* из условий
Q°p (RI О) = 0£; |
02, (Дх, ±h*) = Qp. |
(3.21) |
Граничные условия для уравнений (3.20) |
|
|
0J,(r, |
±h *) = 0; |
(3.22а) |
02,(Д;, z) = 0 |
(3.226) |
выберем из условия, что на торцах и внешней поверхности выде ленного цилиндра приращение температуры 0р = 0р,
— B(rdz ±h) ± H 3Q1B(r, ± h ) = 0; |
|
(3.23a) |
|
f t * |
f t * |
|
|
j* ©),(/?!, z)dz — J 0 3 (7?!, z) dz; |
|
(3.236) |
|
- f t * |
- f t * |
|
|
a e j, (R1: z) |
(Дi, z) |
dz. |
(3.23b) |
a.r |
dr |
|
|
Постоянные Csn определялись из условия минимума средне квадратичной ошибки уравнений (3.19)
Ся |
Gpco^og2 sin р„ |
(3.24) |
|
где /* — усредненное по объему значение функции / (г, z). Постоянные С1п, С2п и Ст определялись из линейной системы
3п алгебраических уравнений (3.16), получаемой при выполнении граничных условий.
В качестве приближенного решения уравнений (3.20), удовлет воряющего условиям (3.22а) и (3.23а), выберем
СО |
|
vnz |
(3.25) |
|
0ф = 2 \j>mh ( ^ |
) + D2nK0( ^ |
|||
) + />з„] cos h* ’ |
||||
П=1 |
|
|
|
|
СО |
|
|
(3.26) |
|
еЬ = 2 |
0 - Л ( - ^ ) “ |
з ^ 1 . |
я (2п — 1)
где vn:
79
Неизвестные коэффициенты D Sn определим из условия мини мума среднеквадратичной ошибки уравнений (3.20)
р. |
_ |
|
С0шб2С sin у п |
(3.27) |
|
зп ~ |
„ „ |
|
sin 2v, |
||
|
|
||||
|
|
2 л К р\ п (^Н |
|
|
|
где |
На |
h* |
|
|
|
|
|
|
|
||
|
j |
I |
(бр —0*) / (г, z ) r d r d z |
|
|
* = |
R . |
-ft*________________________ |
|
||
|
|
|
(Д*а-Д})А* |
|
Постоянные D ln, D in и D4„ определяются из граничных усло вий (3.126, г, д), что приводит к линейной системе 3п уравнений:
|
|
, |
• hnh* |
|
П lr ( |
/св sin —- — |
|
Din^i |
v'1^ 1 \ Л |
kp sin v„ |
|
U 2 nK\ |
} — u in |
Пример 2. Определим температуру нагрева в рассматриваемом цилиндрическом амортизаторе, геометрические размеры которого, а также параметры режима нагружения, реологические и тепло-
Рис. 56. Распределение температур в полом цилиндрическом амортизаторе: (------------- расчет; точки — эксперимент)
физические характеристики оставались прежними (см. пример 1 ). Зависимость ф (Т) была такой, как на рис. 55, температурной зависимостью мгновенного модуля пренебрегаем.
Параметры и 0* определялись путем обработки методом наименьших квадратов экспериментальной зависимости техни
80