Файл: Потураев, В. Н. Резина в горном деле.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 72

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

и удовлетворением граничных условий получаем систему алгебра­ ических уравнений

 

-С ,„[ц А

 

 

 

( £^ -)]+ Я ,4 С а„ = 0,

 

п Т

( Ц/1^1 \ ! г*

( № 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* из условий

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