CC BY
57
ТЕХНИЧЕСКИЕ НАУКИ
DOI 10.23859/1994-0637-2017-1-76-1 УДК 662.289:678.063
Аваев Александр Алексеевич
Кандидат технических наук, доцент, Вологодский государственный университет (Вологда, Россия) E-mail: iur.osipov2011@yandex.ru
Осипов Юрий Романович
Доктор технических наук, профессор, Вологодский государственный университет (Вологда, Россия) E-mail: iur.osipov2011@yandex.ru
О РАСЧЕТЕ ВРЕМЕНИ ПРОГРЕВА ЭЛАСТОМЕРНОГО ПОКРЫТИЯ НА ТКАНИ
Аннотация. Работа посвящена расчету времени прогрева эластомерного покрытия на ткани. Предлагаемый метод основан на аналитической математической модели внутреннего теплопереноса в эластомер-ном покрытии. Математическая модель построена на линейном уравнении теплопроводности применительно к симметричным граничным условиям III рода. Результаты расчетов проиллюстрированы графически. Сравнение расчетов с данными опытов свидетельствует о возможности применения математической модели в инженерной практике. Представленный метод может быть использован для совершенствования технологии термической вулканизации.
Ключевые слова. Аналитическая модель, внутренний перенос тепла, эласто-мерное покрытие на ткани, линейное уравнение теплопроводности, граничные условия, термическая вулканизация.
© Аваев А.А., Осипов Ю.Р., 2017
Avaev Alexander Alexeevich
PhD (Technology), Associate professor, Vologda State University (Vologda, Russia) E-mail: iur.osipov2011@yandex.ru
Osipov Yuriy Romanovich
Doctor of technical sciences, professor, Vologda State University (Vologda, Russia) E-mail: iur.osipov2011@yandex.ru
ABOUT WARM-UP TIME CALCULATION FOR ELASTOMERIC COATING ON FABRIC
Abstract. This work is focused on the warm-up time calculation for the elastomeric coating or fabric. Suggested method is based on analytic mathematic model of inside heat transfer in the elastomeric coating. Mathematic model is built on linear equation of heat transfer for symmetric boundary conditions of III type. The calculation results are illustrated by graphs. The comparison of the calculations and the experiment data testify to possibility of analytic model application in engineering practice. The method presented in the paper may be used for the improvement of thermal vulcanization technology.
Keywords. Analytic model, inside heal transfer, elastomeric coating on fabric, lineal equation of heat transfer, boundary conditions, thermal vulcanization.
Введение
Процесс термической вулканизации эластомерных покрытий на тканевой основе представляет собой один из самых сложных тепловых процессов, протекающих при изменяющихся во времени тепловых потоках и тепловых эффектах, которые приводят к изменению теплофизических свойств эластомера. Процесс вулканизации, кро-
ме того, зависит от организации и режима теплообмена между объектом и внешней средой.
При вулканизации в индукционный период происходит течение эластомера, а при нагреве покрытия с большим содержанием вулканизующего агента выделяется теплота вследствие экзотермической реакции вулканизации. В ходе вулканизации необходимо обеспечить такое распределение и изменение температуры во времени, при котором для данного типа эластомера свойства готового продукта оказались бы наилучшими при оптимальной организации всего технологического процесса [4].
Выбор режима вулканизации основан на анализе температурных полей в эласто-мерном покрытии и на объективной оценке результатов их воздействия на процесс вулканизации. При этом важно правильно спрогнозировать период предварительного нагрева вулканизуемого объекта, предшествующего фазе течения эластомера, в частности: продолжительность этого периода и распределение температуры эластомера на момент его окончания. Возможность подобного прогнозирования может быть обеспечена путем математического моделирования этапа прогрева системы «ткань - эластомер».
Основная часть
Так как в рассматриваемом случае эластомер полностью пропитывает тканевую основу, на которую он наносится, следовательно, система «ткань - эластомер» имеет практически одинаковые теплофизические свойства [1], [2], то в качестве модели такой системы может быть предложена однородная бесконечная пластина, а ее теплообмен с внешней средой осуществляется по закону Ньютона.
В данном случае внутренний перенос тепла в системе «ткань - эластомер» описывается известным линейным уравнением теплопроводности [1]:
^М = «^М (о < X </)
2 , V- ....., (1)
дт дх
где t (х, т) — температура точки пластины толщины / с координатой х в момент времени т ; а — коэффициент температуропроводности системы «ткань - эластомер», м2/с.
При решении уравнения (1) использованы следующие краевые условия
t (х,0) = / (х); (2)
^ + (т) — t (0, т)]=0; (3)
+ а[Ч (т) — t (/,т)] = 0,
(4)
где X — коэффициент теплопроводности системы «ткань - эластомер», Вт/(моС); а — коэффициент теплоотдачи от внешней среды с температурой tc (т), в общем случае зависящей от времени т , к поверхностям х = 0 и х = / системы, Вт/(м2оС).
Условие (2) является начальным, оно соответствует распределению температуры в системе «ткань - эластомер» в начальный момент времени т = 0 и является в общем случае функцией координаты х.
Условия (3), (4) относятся к числу граничных условий, они отвечают теплообмену по закону Ньютона на поверхностях х - 0, х - I системы (граничные условия III рода).
Решение уравнения (1) при краевых условиях (2)-(4) получено с помощью конечных интегральных преобразований [3] и может быть представлено в следующем виде [2]:
< (х, ,)=, (х)+2]г .[1+(—1Г ]. , 'С08 ^ Х++ (ц),
ц 2т+2В1 + В12
(5)
где
Рт (т)= Рт (х) + Ц / (х).
Цт • | Цт • + В • ¡ЯП [ц. ' у
Рт (Т) =
| В .Цт + "2 .(^т + В 2 )• ^ Цт
Л (®).
.ёх е
' 1 ёю;
В-а-^ I — критерий Био; X - х/1 — безразмерная координата точки пластины; ю- переменная интегрирования; Цт > 0 (т -1,2,...) — корни характеристического уравнения
1 Г ц В1Л
С^ = 2 . [ В - В |
На рис. 1 представлены результаты расчета нестационарного температурного поля системы «ткань - эластомер», произведенные в соответствии с решением (5). Расчеты проводились применительно к условиям прогрева системы для эластомера марки 2566 (НК+СКБ) в аппарате непрерывного действия с псевдоожиженным слоем инертного зернистого теплоносителя при следующих параметрах [1]: I - 0,001 м;
а - 0,71 10—7 м2/с; Х-0,705 Вт/(моС); а-100 Вт/(м2оС); ГС (т) -150 оС; /(х) = 30оС.
Линии с номерами 1, 2, 3, 4 соответствуют следующим значениям текущего времени т : 1 - 0,1 с; 2 - 40 с; 3 - 80 с; 4 - 120 с. Из рис. 1 следует, что линии нестационарного температурного поля с течением времени становятся короче, превращаясь в итоге (т.е. при выходе на стационарный режим) в отрезок прямой линии, который параллелен оси абсцисс. Длина отрезка прямой при этом равна 1.
На этом принципе строится зависимость длин Ь линий поля от безразмерного
критерия времени Фурье - Ро - а .т/12, т.е. Ь - Ь(Ро). Пример такой зависимости
представлен на рис. 2.
Из рис. 2 следует, что стационарный режим нагрева системы «ткань - эластомер» наступает при значении критерия Фурье, составляющем примерно 4,0. Конечно же, прогрев рассматриваемой системы следует рассматривать не обязательно до наступления стационарного режима, а, например, до начала периода фазовых превращений, когда при дальнейшем расчете линий нестационарных температурных полей следует учитывать появление источников (стоков) тепла [1], [2].
а
t(x,r);c
160
120
80
40
L,'C 30
20
10
1,0
2,0
0 0,2 0,4 0,6 0,8 1,0 X
4,0
6,0
Рис. 1. Пример расчета температурного поля системы «ткань - эластомер» в соответствии с решением (5)
Рис. 2. Пример расчетной зависимости
длин линий нестационарного температурного поля от критерия Фурье
0
0
Предположим, что нас интересует значение критерия Фурье, соответствующее достижению максимального (среднего интегрального или минимального) значения температуры, составляющего, например, 120 оС. Это также может быть получено расчетным путем (см. рис. 3).
115
110
FOmax Fosean Fo^
1,10 1,15 1,20 1,15 1,30 Fo
Рис. 3. Пример расчетов значений Fomax, Fomeím, Fo, при заданной температуре, составляющей 120 оС
muí
ог
На рис. 3 кривыми 1-3 представлены графики зависимостей от критерия Фурье, соответственно, максимального, среднего интегрального и минимального значений температуры для линий нестационарного температурного поля. Прямая линия 4, параллельная оси абсцисс, соответствует интересующему нас значению температуры -120 оС. Абсциссы точек пересечения этой прямой с линиями 1-3 представляют собой значения критерия Фурье, соответствующие моменту достижения максимального (Ротах), среднего интегрального (Ротеа„) и минимального (Ротт) значений температуры на одной из линий температурного поля. Следует особо отметить, что все зависимости, представленные на рис. 2 и 3, могут быть автоматически рассчитаны при минимальной затрате времени. Так, в частности, расчет и построение графиков
таких зависимостей при работе с Windows-приложением MathCAD составляет всего лишь доли секунды.
Представленные расчеты апробированы экспериментально, а предлагаемый метод может быть рекомендован для применения в инженерной практике при разработке оптимальных методов организации технологического процесса термической вулканизации эластомерных покрытий на тканевой основе.
Выводы
Необходимое качество эластомерных покрытий на тканевой основе обеспечивают подбором такого сочетания параметров процесса термообработки, при котором возможна равномерность вулканизации покрытий. Она осуществляется выбором рациональных условий внутреннего и внешнего теплообмена, а также за счет конструкций вулканизационных аппаратов и варьирования характеристик эластомеров.
Особое внимание следует уделять первому этапу вулканизации - процессу прогрева эластомерного покрытия до момента начала фазы изменения агрегатного состояния эластомера, что, как правило, сопряжено с изменением его теплофизических свойств.
Литература
1. Аваев А.А. Математическое моделирование нестационарной теплопроводности при вулканизации эластомерных покрытий на тканевых и металлических подложках. Череповец, 2008. 133 с.
2. Аваев А.А., Осипов Ю.Р., Павлов В.В. Математическая модель предварительного нагрева эластомерного покрытия на тканевой основе при его вулканизации // Химическое и нефтегазовое машиностроение. 2008. №12. С. 10-11.
3. Карташов Э.М. Аналитические методы в теории теплопроводности твердых тел. М., 1985. 480 с.
4. Лукомская А.И., Баденков П.Ф., Кеперша Л.М. Тепловые основы вулканизации резиновых изделий. М., 1972. 360 с.
References
1. Avaev A.A. Matematicheskoe modelirovanie nestatsionarnoi teploprovodnosti pri vulcanizat-sii elastomernykh pokrytii na tkanevykh i metallicheskikh podlozhkakh [Mathematic models of the instationary heat conduction during the vulcanization of the elastomer coats on the fabric and metal bases]. Cherepovets, 2010, 133 p.
2. Avaev A.A., Osipov Iu.R., Pavlov V.V. Matematicheskia molel' predvaritel'nogo nagreva elastomernogo pokrytiia na tkanevoi osnove pri ego vulkanizatsii [Mathematic model of the elastomer coat on fabric base preliminary heating in it's vulcanization]. Himicheskoe i neftegazovoe ma-shinostroenie [Chemical and oil-gas machine building], 2008, vol. 12, pp. 10-11.
3. Kartashov E.M. Analiticheskie metody v teorii teploprovodnosti tvierdykh tel [Analytic methods in solid bodies heat conduction theory]. Moscow, 1985, 480 p.
4. Lukomskaia A.I., Badenkov P.F., Kepersha L.M. Teplovye osnovy vulkanizatsii rezinovykh iz-delii [Heat bases of the rubber articles vulcanization]. Moscow, 1972, 360 p.
Аваев А.А., Осипов Ю.Р. О расчете времени прогрева эластомерного покрытия на ткани // Вестник Череповецкого государственного университета. 2017. №1. С. 9-13.
For citation: Avaev A.A., Osipov Y.R. About warm-up time calculation for elastomeric coating on fabric. Bulletin of the Cherepovets State University, 2017, no. 1, pp. 9-13.
