Математическая модель теплопереноса в процессе термической вулканизации покрытия из эластомера на металлической пластине Текст научной статьи по специальности «Механика и машиностроение»

ТЕХНИЧЕСКИЕ НАУКИ

DOI 10.23859/1994-0637-2018-1-82-1 УДК 662.289:678.063

Аваев Александр Алексеевич

Кандидат технических наук, доцент, Вологодский государственный университет (Вологда, Россия)

Осипов Юрий Романович

Заслуженный деятель науки Российской Федерации, доктор технических наук, профессор, Вологодский государственный университет (Вологда, Россия) E-mail: iur.osipov2011@yandex.ru

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ТЕПЛОПЕРЕНОСА В ПРОЦЕССЕ ТЕРМИЧЕСКОЙ ВУЛКАНИЗАЦИИ ПОКРЫТИЯ ИЗ ЭЛАСТОМЕРА НА МЕТАЛЛИЧЕСКОЙ ПЛАСТИНЕ

Аннотация. Данная работа представляет собой очередную попытку создать аналитическую модель переноса тепла в процессе термической вулканизации слоя эластомера, нанесенного на металлическую пластину. Данная задача связана с наличием тепловыделений и изменением тепловых параметров эластомера в течение процесса вулканизации. Несмотря на явную нелинейность этой задачи, была использована идея применения линейного уравнения теплопроводности. Теплообмен в системе «эластомер - металл» исследован применительно к процессу контактной термической вулканизации. Результаты расчетов, основанных на предлагаемой модели, иллюстрированы графически. Данный метод расчетов может быть использован в инженерной практике для совершенствования технологии термической вулканизации.

Ключевые слова: аналитическая модель, теплопроводность, слой эластомера на металлической пластине, линейное уравнение теплопроводности, граничные условия, контактная термическая вулканизация_

© Аваев А.А., Осипов Ю.Р., 2018

Avaev Alexander Alexeevich

PhD in Technical Sciences, Associate Professor, Vologda state university (Vologda, Russia)

Osipov Yury Romanovich

Hounoured worker of science, Doctor of technical sciences, professor, Vologda State University (Vologda, Russia) E-mail: iur.osipov2011@yandex.ru

MATHEMATIC MODEL OF THE HEAT TRANSFER DURING THERMAL VULCANIZATION OF ELASTOMER COAT ON THE METAL PLATE

Abstract. This work is a next attempt to create the analytic model of heat conduction in the thermal vulcanization process of elastomer coat on the metal plate. The problem is connected with heat releases and the change of elastomer heat parameters during the vulcanization process. Despite the evident non-linearity of this problem the idea of linear heat conduction equation application was used. Heat transfer in the elastomer-metal system was investigated conformably to contact thermal vulcanization process. The calculation results which based on the proposed model are illustrated by graphics. This method of calculations may be used in the engineer practice for the improvement of thermal vulcanization technology.

Keywords: analytic model, heat conduction, elastomer coat on metal plate, linear equation of heat conduction, boundary conditions, contact thermal vulcanization

Введение

Процесс термической вулканизации эластомерных покрытий на металлической пластине представляет собой сложный тепловой процесс, протекающий при изменяющихся во времени тепловых эффектах, которые приводят к изменению теплофи-зических свойств эластомера. Процесс вулканизации, кроме того, зависит от организации и режима теплообмена между объектом и внешней средой.

При вулканизации в индукционный период происходит течение эластомера, а при нагреве покрытия с большим содержанием вулканизующего агента выделяется теплота вследствие экзотермической реакции вулканизации. В ходе вулканизации необходимо обеспечить такое распределение и изменение температуры во времени, при котором для данного типа эластомера свойства готового продукта оказались бы наилучшими при оптимальной организации всего технологического процесса [1].

Выбор режима вулканизации основан на анализе температурных полей в эласто-мерном покрытии и металлической пластине. Прогнозирование оптимального режима процесса вулканизации может быть осуществлено путем его математического моделирования, при этом достоверность и надежность работы подобной модели играет первостепенную роль.

Основная часть

В данной работе рассматривается процесс так называемой контактной термической вулканизации, когда на внешних поверхностях эластомерного покрытия и металлической подложки поддерживается заданная температура. Весь процесс разбивается на несколько этапов, в течение каждого из которых изменение теплофизиче-ских свойств эластомера, а также мощности действующего источника тепла, являются незначительными, а следовательно, могут считаться постоянными. При этом распределение температуры на момент окончания очередного этапа является начальным условием при описании следующего этапа. Такой подход к моделированию процесса вулканизации позволяет использовать линейное уравнение теплопроводности, решая аналитически систему данных уравнений для эластомера и металла.

Как известно, аналитическое решение в отличие от любого его численного варианта обеспечивает надежность сходимости соответствующих рядов Фурье [2] и позволяет достаточно точно оценивать погрешности производимых вычислений.

В качестве основы модели рассматривается бесконечная двухслойная пластина общей толщины 11 +12, при этом диапазон изменения координаты х точек пластин

описывается неравенством -11 < х < 12; индекс г = 1 соответствует левой пластине, а

индекс г = 2 - правой.

Линейные уравнения теплопроводности для левой и правой пластины имеют следующий вид

д /1 (х, т) д111 (х, т) м>]

--- = а1-2— +

дх дх

(-/1 < х < 0)

СА

(1)

д/2 (х, т) д2/2 (х,т)

дх 2 дх2 с2р2

(0 < х < /2)

(2)

где т - текущее время, с; (х, т) - температура, оС; м>1 - удельная мощность источника тепла, Вт/м3; а1 - коэффициент температуропроводности, м2/с; - удельная теплоемкость, Дж/(кг оС); р( - плотность, кг/м3.

Для решения системы уравнений (1)-(2) использованы следующие краевые условия

tj (x,0) = ajx2 +Pjx + ур

(ctj, Pj, yj = const) t2 (x,0) = a2x2 +P2x + у2; (a2, P2, у 2 = const) tj (0, t) = /2 (0,t); A gtj (0,t) a dt2 (0,t) ,

Aj - = Л2 - ;

dx dx

tj (-lj,t) = tj = const; tj (l2, t) = t2 = const,

(3)

(4)

(5)

(6)

(7)

(8)

где X1 (г = 1; 2) - коэффициент теплопроводности, Вт/(м оС).

Начальные условия (3), (4) представляют собой многочлены второй степени координаты х , которые наилучшим образом описывают линии температурного поля, а, кроме того, позволяют легко их аппроксимировать, используя метод наименьших квадратов. Как показывают расчеты, при средней абсолютной погрешности аппроксимации, составляющей около 1,0 оС, средняя абсолютная погрешность, которая получается при работе с подобными начальными условиями, составляет всего лишь 0,01 - 0,03 оС.

Граничные условия (5), (6) соответствуют идеальному тепловому контакту, а (7), (8) являются условиями I рода, что соответствует контактной схеме вулканизации эластомера.

Изображение по Лапласу решения системы уравнений (1)-(2) при данных краевых условиях может быть представлено в следующем виде

T (x, s) = Cj (s) • ch

( ^

. a,

4 1 j

- C2 (s)•sh

(

\

—x

a.

+ aj x + Pj x + у j + 2ajaj + Wj

(9)

( T ^

T2 (x, s) = Cj (s) • ch I—x + KjC2 (s)• sh

4 V a2

2 (ajaj - a2a2) + Wj - W2

. a2

4 » 2 /

ch

. a2

4 1 2

ch

1 A

Pj-P2

-sh

( T ^

x

4 . a

a2 x2 + P2 x + у 2 2a2 a2 + W2

(j0)

где = (с,.р,.) (г = 1; 2) - приведенная удельная мощность источника тепла, оС/с; 5 - аргумент по Лапласу.

При этом коэффициенты С1 (5), С2 (5) рассчитываются по формулам

С (5 )= Р11 (5) + Р12 (5 ) ; С2 (5 )= Р21 (5 ) + Р22 (5 ) ,

5-9(5) 52 ф (5)' 5-9(5) 5 2ф (5)

где

^11 (5 ) =

/2-а2/22-Р2/2-у2-(У1 -у2)сй ^2^1 "|т1Р1 -Р2 1/2

5к

а,

к ^ -/1 а1

як

Г I— А

у. а1 у

( I— А

К2. ^ У V а1 У

5к

+ к -(/1 -а^2 +Р1/1 -У1);

( I— Л

к,. -/.

5к

^12 (5) = - |[2 (а1а1 - а2а2) + - ] ск

- К1 - (2а1а1 + ^ ) ;

+ 2 а. а2 +

^ I— А

О*

у.а1 у

1"2 1 " 2 ( / I- Д

К2 — /1

^21 (5 ) =

/2 - а2/22 - Р2/2 -У2 - (У1 -У2 ) ск

( I— Л

5

Р1-Р21 /2-

А-)

5к

ки1 -/,

л

КгА~/1

ки -К

а>

ск

-( -а1/12 +Р1/1 - у1 )ск К2 ~~~/1

( I— А

уа1 у

5к

ки а

а1

Р22 (5) = -

[2 (а1а1 -а2 а2 ) + - Г2 ] ск

+ 2а2 а2 +

- (2а1а1 + РТ1 )ск

Ки -/,

а1

К и -/,

а1

ск

( I— А

а/1

Уа1 У

5к

Ки -к

а1

X

9(5) = K1ch

( I— Л

s

v < 1 У

+ ch

К

sh

f I— A

v. a1 у

i У

sh

( i— V 1

Kj-l,

v V ai у

a2 A, ' ai ^2

В данной ситуации, когда переход от изображения решения к представлению его оригинала в виде формул требует громоздких преобразований, выгоднее для каждого слагаемого под номером у вида

Ф у ( s )

sk 9(s )

(k = 1; 2)

применять формулу Ващенко-Захарченко [2]

tj (x, т) =

-—lim

(k -1)! s^0

jk-1 (

dsk

Ф;М (

9(s )

m=1 4 (Sm )

где v(s) = sk9(5); sm = -ц2m~2r, а цm (m = 1,2,3...) - положительные корни харак-

l1

теристического уравнения

sin ц • cos (К1ц) + K1 sin (К1ц)- cos ц = 0.

На рисунке представлены результаты расчетов температурного поля, проведенные по предлагаемой методике.

Г, °С 160

140

120

100

80

60

40

20

-4 - 3 - 2 -1 0 1 2 3 х ■ 103, м

Рисунок. Результаты расчета линий температурного поля эластомерного покрытия на стальной пластине

Расчеты проведены применительно для слоя эластомера марки 2566 (НК+СКБ), нанесенного на пластину из стали Ст.3. При этом толщина слоя эластомера составляла l1 = 0,004 м, а стальной пластины - l2 = 0,003 м. Температура внешних поверхностей пластин была 150 оС, а температура на момент начала прогрева системы с последующей вулканизацией слоя эластомера была 30 оС. Приведенным на рисунке линиям 1, 2, 3, 4 соответствовали значения времени т = 5, 10, 15, 20 с.

Представленные расчеты апробированы экспериментально, а предлагаемый метод может быть рекомендован для применения в инженерной практике при разработке оптимальных методов организации технологического процесса термической вулканизации эластомерных покрытий на металлических подложках.

Выводы

Необходимое качество эластомерных покрытий на металлических пластинах обеспечивают подбором такого сочетания параметров процесса термообработки, при котором достигается равномерность вулканизации покрытий, что осуществляется путем назначения рациональных условий внутреннего и внешнего теплообмена, а также за счет разработки оптимальных конструкций вулканизационных аппаратов.

Возможность проведения соответствующих расчетов, несмотря на их кажущуюся громоздкость и сложность, вполне реальна при использовании современных систем программного обеспечения и может быть рекомендована не только для проведения научных исследований, но и в инженерной практике.

Литература

1. Лукомская А.И., Баденков П.Ф., Кеперша Л.М. Тепловые основы вулканизации резиновых изделий. М., 1972. 43 с.

2. Лыков А.В. Теория теплопроводности. М., 1952. 600 с.

References

1. Lukomskaia A.I., Badenkov P.F., Kepersha L.M. Teplovyie osnovy vulkanizatsrn rezinovykh izdelii [Heat bases of the rubber articles vulcanization]. Moscow, 1972. 43 р.

2. Lykov A.V. Teoriia teploprovodnosti [Heat conduction theory]. Moscow, 1952. 600 р.

Аваев А.А., Осипов Ю.Р. Математическая модель теплопереноса в процессе термической вулканизации покрытия из эластомера на металлической пластине // Вестник Череповецкого государственного университета. 2018. №1 (82). С. 14-19.

For citation: Avaev A.A., Osipov Y.R. Mathematic model of the heat transfer during thermal vulcanization of elastomer coat on the metal plate. Bulletin of the Cherepovets State University, 2018, no. 1 (82), pp. 14-19.