Другие журналы

научное издание МГТУ им. Н.Э. Баумана

НАУКА и ОБРАЗОВАНИЕ

Издатель ФГБОУ ВПО "МГТУ им. Н.Э. Баумана". Эл № ФС 77 - 48211.  ISSN 1994-0408

Моделирование и расчет процессов индукционного нагрева прессового оборудования при производстве резинотехнических изделий

# 03, март 2013
DOI: 10.7463/0313.0541632
Файл статьи: Карпушкин_P.pdf (961.37Кб)
авторы: Малыгин Е. Н., Карпушкин С. В., Карпов С. В.

УДК 62-932.4

ФГБОУ ВПО "Тамбовский государственный технический университет"

malygin_en@mail.ru

karp@mail.gaps.tstu.ru

g7302@yandex.ru

 

Введение

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

Качество РТИ во многом определяется равномерностью температурного поля в объёме изделия. На равномерность поля оказывают влияние многочисленные факторы, прежде всего – параметры пресс-формы, системы её обогрева и внешнего теплообмена, а также расположение пресс-формы на поверхности плиты, теплофизические свойства материалов плиты и пресс-формы, параметры системы автоматического управления (САУ) температурой нагревательных плит.

До настоящего времени в литературе рассматривались, главным образом, задачи обеспечения температурного поля заданного профиля на рабочей поверхности нагревательных плит прессов, причем чаще всего требуется обеспечить равномерное поле с погрешностью ±3÷оС [1, 2]. Между тем, целью нагрева плит пресса является обеспечение равномерного температурного поля в вулканизируемом изделии, однако постановка и решение трёхмерных задач моделирования нестационарного нагрева системы «плита–пресс-форма–изделие» с учётом нелинейных зависимостей свойств материалов и мощности индукторов от температуры в литературе практически не рассматриваются из-за сложности математического описания (система дифференциальных уравнений в частных производных).

 

Описание: D:\Мои документы\Documents\ТГТУ\Плита\Картинки\Platen_named3.jpg

 

Рис. 1 – Конструкция плиты с индукционными нагревателями

 

В данной работе предлагается методика расчета температурных полей одиночной плиты индукционного нагрева и системы «плита индукционного нагрева–пресс-форма–изделие» с применением сравнительно нового инструмента – системы конечноэлементного анализа COMSOL Multiphysics[3]. Поскольку разработчики подобных систем не предоставляют информацию об используемом математическом описании моделируемого процесса и методике решения математических моделей, для проверки адекватности методики и результатов применения системы COMSOL Multiphysics производится сравнение результатов с данными  промышленного эксперимента по исследованию процесса нагрева одиночной индукционной плиты, а также с результатами решения разработанной математической модели процесса нагрева этой плиты методом конечных интегральных преобразований [4, 5].

Цель работы – разработка методик тепловых расчётов индукционных нагревательных плит прессов, осуществляющих процесс вулканизации РТИ в металлических пресс-формах в соответствии с требованиями технологического регламента.

Область использования результатов – автоматизированный расчёт существующего и проектирование нового прессового оборудования вулканизационных гидравлических прессов для производств РТИ.

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

1. Математическая модель процесса нагрева одиночной индукционной плиты для периода достижения рабочей температуры

Процесс нагрева одиночной индукционной плиты описывается трёхмерным нестационарным дифференциальным уравнением теплопроводности в частных производных

(1)

где T(x,y,z,τ) – температура плиты (оС) в точке ее объема с координатами (xyz) в момент времени τ;

       a = λs/(cs·ρs) – коэффициент температуропроводности материала плиты;

       λs, сs, ρs – теплопроводность (Вт/м·К), теплоемкость (Дж/кг·К) и плотность (кг/м3) материала плиты соответственно;

       – средняя температура (оС) плиты в момент времени τ;

      l, s, h – длина, ширина и высота плиты соответственно (м), т.е. x[0;l], y[0;s], z[0;h];

       – зависимость интенсивности внутренних тепловыделений от координат и средней температуры плиты (Вт/м3);

       q0j – удельная мощность j-го индуктора при средней температуре плиты, Вт/м3;

       vj – объем, в котором выделяется мощность j-го индуктора (м3);

       ni – число индукторов в плите.

Для определения значений q0j, j=1,…,ni необходимо решать уравнение Максвелла [6, 7] и находить распределение магнитной индукции вокруг пазов индукторов, причем в [7] показано, что если требуемая температура нагрева плиты не превосходит температуру Кюри (750 оС), то задачи определения мощности, выделяемой в пазах индукторов, и распространения тепла в материале плиты можно решать независимо друг от друга. Решение уравнения Максвелла для объекта конечных размеров при граничных условиях, определяемых размерами и конфигурацией пазов индукторов, затруднительно, поэтому будем считать, что тепловыделение каждого индуктора осуществляется равномерно в объеме паза под него, причем характеристики материала этого объема, заполняемого катушкой индуктора и специальной пастой (например, органосиликатная композиция ОС-82-05 [8]), соответствуют характеристикам материала плиты.  

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

где hk – высота крышки плиты (м),

       gj, bj– глубина и ширина паза под  j-й индуктор (м);

       lj, sj – длина и ширина j-го индуктора (м);

       xcj, ycj – абсцисса и ордината геометрического центра j-го индуктора (м), причем за точку (0,0) принят левый нижний угол плиты.

Следовательно , где Qij – мощность j-го индуктора. Значения Qij, j=1,…,niпри температуре Tav будем определять по методике [7], основанной на экспериментальных исследованиях индукционного нагрева ферромагнитной стали. Эта методика предусматривает решение системы уравнений:

,                                                            (2)

где Ij– ток j-го индуктора (А),

      Uj– напряжение на j-м индукторе (В),

 –                                                          (3)

коэффициент мощности индукционного нагрева,

 –                (4)

реактивная мощность j-го индуктора (Вт),

f­– частота переменного тока в индукторах,

 –                                         (5)

магнитодвижущая сила j-го индуктора (ампер-витки),

 –                                       (6)

напряженность магнитного поля в пазу j-го индуктора (А/м),

 –                                                                  (7)

КПД индукционного нагрева,

 –                                                (8)

удельное электросопротивление материала плиты (Ом∙м) (функция получена в результате аппроксимации экспериментальной графической зависимости [7]),

 –                                     (9)

относительная магнитная проницаемость материала плиты,

 –                                         (10)

активная мощность j-го индуктора (Вт),

       di – диаметр провода индукторов (м),

       ri– удельное электросопротивление материала провода индукторов (Ом∙м),

 –                                                              (11)

число проводов в j-м индукторе.

Начальное условие для решения уравнения (1) имеет вид

,                                                             (12)

где Т0 – температура окружающего воздуха, оС.

Форма записи граничных условий определяется конструкцией пресса. При отсутствии изоляции торцов нагревательной плиты, а также контакта ее крышки с нажимной плитой, используются граничные условия (ГУ) третьего рода:

                                 (13)

Здесь ,  – комплексные коэффициенты теплоотдачи от торцов плиты,  α5 = Вн ∙ α(tст),  α6 = Вв ∙ α(tст) –коэффициенты теплоотдачи от рабочей поверхности плиты и поверхности ее крышки соответственно (Вт/м2∙К);

          α(tст) = αконв(tст) + αизл(tст) – коэффициент теплоотдачи от поверхностей плиты к окружающему воздуху конвекцией и излучением, значение которого можно определить согласно [9];

         tст – температура стенки плиты (оС);

         nкs, nкl  – число крепежных пластин в торцах плиты, соответствующих  ее ширине и длине, соответственно;

         λк – теплопроводность материала крепежных пластин плиты (Вт/м·К);

         sк, hкlк – ширина, высота сечения пластин и длина их части, выходящей за пределы плиты, соответственно (м);

         Вн, Вв – коэффициенты, учитывающие неодинаковые условия теплоотвода от нижней и верхней поверхностей плиты, соответственно.

Параметрами (независимыми переменными) математической модели (1)-(13) являются:

-  количество индукторов ni;

-  положение центров индукторов xc, yc, (м);

-  геометрия индукторов, т.е. длина lj и ширина sj каждого из них, ширина bj и глубина gjпаза (м),

-  число витков индукторов nj.

При решении модели (1)-(13) необходимо учитывать следующие ограничения:

- в момент времени th, когда температура плиты в месте установки контрольной термопары достигнет заданной, значения температуры в точках рабочей поверхности плиты T(x,y,h,τh) могут отличаться от требуемых T*(x,y,h) не более чем на заданную величину T(x,y), т.е.

                            (14)

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

;                                                           (15)

- значения геометрических параметров индукторов ограничены размерами плиты, правилами их взаимного расположения, а также возможностями изготовления, т.е.

; (16)

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

,                                              (17)

где χ – коэффициент заполнения паза;

- положения центров индукторов должны обеспечивать расстояния от краев плиты до индукторов не менее 5 см, т.е.

;  (18)

- суммарная мощность индукторов не должна превышать заданной максимальной мощности нагревательной плиты Qpmax Вт, т.е.

;                                                            (19)

- коэффициент мощности и КПД индукционного нагрева не должны быть меньше минимально допустимых значений (cosφmin и ηmin соответственно), т.е.

.                                               (20)

Решение модели (1)-(20) осложняется нелинейными изменениями во времени мощностей индукторов и коэффициентов теплоотдачи от поверхностей плиты. Поэтому будем считать, что коэффициенты теплоотдачи от всех поверхностей плиты и значения мощностей индукторов не зависят от температуры в пределах некоторых интервалов времени, которые названы расчётными.

Таким образом, при разработке методики решения математической модели индукционного нагрева плиты вулканизационного пресса (1)-(20) приняты следующие допущения.

1) Тепловыделение каждого индуктора осуществляется равномерно во всем объеме паза плиты, в котором размещен индуктор. Характеристики материала этого объема соответствуют характеристикам материала плиты.

2) Значения параметров α1, α2,…, α6 и Qij, j = 1,…,ni не зависят от температуры в пределах расчетных интервалов времени.

3) Значения мощности индукторов при изменении средней температуры плиты определяются по соотношениям (2)-(11).

С учётом этих допущений задача распространения тепла в плите при фиксированных значениях параметров α1, α2,…, α6 и Qij решается многократно, т.е. последовательно находятся решения для расчётных интервалов времени, в которых эти параметры можно считать константами. При этом начальным условием для всех решений, кроме первого, будет температурное поле плиты, соответствующее моменту окончания предыдущего интервала времени.

 

FIT_synopsis

 

Рис. 2 – Алгоритм решения модели (1)-(20) с помощью МКИП

 

 Алгоритм МКИП, использованный для решения модели (1)-(20), представлен на рис. 2 (подробную методику его применения см. в [10]). В качестве условий окончания предыдущего расчетного интервала времени нагрева и начала следующего использовались неравенства

 и ,

где , – значения k-го коэффициента теплоотдачи и мощности j-го индуктора для предыдущего расчетного интервала времени, ,  – их текущие значения (обоснование см. в [11]). Окончание текущего расчетного интервала времени регистрировалось при выполнении любого из этих условий.

Алгоритм метода конечных элементов (МКЭ), реализованный в системе COMSOL Multiphysics, представлен на рис. 3. Значение шага по времени Δτ = 2 с выбрано методом Рунге. Для расчета внутренних тепловыделений индукторов использовалась аналитическая зависимость Qi(τ), полученная в результате аппроксимации значений Qij, j=1,...,ni, полученных по соотношениям (2)-(11) (ее график представлен на рис. 5).

 

 

Рис. 3 – Алгоритм решения модели (1)-(20) с помощью МКЭ

 

На рис. 4 представлены результаты решения модели (1)-(20) для незакрепленной промышленной плиты размерами 500×410 мм начальной мощности 5,5 кВт с четырьмя одинаковыми индукторами (рис. 1) в процессе нагрева от начальной температуры  Т0 = 12 оС до рабочей температуры 170 оС, а также результаты промышленного эксперимента, проведенного в ОАО «АРТИ-Завод», г. Тамбов по исследованию процесса нагрева той же плитыс помощью нескольких термопар: по углам и в центре рабочей поверхности  (подробнее см. [10]). Результаты приведены для контрольной термопары, расположенной в центре плиты над индукторами. Параметрам модели (1)-(20) были присвоены постоянные значения, соответствующие конструкции исследуемой плиты.

 

 

Рис. 4 – Сравнение результатов решения модели (1)-(20) с данными эксперимента

 

рис_3

 

Рис. 5 – Сравнение результатов измерения и расчёта мощности плиты

 

Сопоставление расчётных и экспериментальных данных (рис. 4) относится к интервалу 0-1968 с (именно за это время температура контрольной термопары в эксперименте достигает заданного значения 170 °С). Максимальное отклонение результатов расчётов от экспериментальных значений температуры контрольной термопары не превышает 5 °С. Результаты измерения мощности плиты в ходе эксперимента в сравнении с результатами расчёта по соотношениям (2)-(11) представлены на рис. 5. Погрешность расчета мощности индукторов не превысила 3 %.

Из данных рис. 4 и 5 можно сделать вывод, что результаты решения модели (1)-(20) с помощью методики, основанной, как на МКИП, так и на МКЭ, удовлетворительно согласуются с данными эксперимента, проведенного в условиях реального производства, причем разница между результатами расчетов с применением МКИП и МКЭ не превышает 8 %. Таким образом, методика применения системы COMSOL Multiphysics для анализа нестационарного процесса индукционного нагрева адекватна и может быть использована для решения более сложных моделей.

 

2. Математическая модельпроцесса нагрева одиночной индукционной плиты для периода поддержания рабочей температуры

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

                    (21)

где Tsensor – температура контрольной термопары, оС;

      Tup,Tdown – верхний и нижний пороги срабатывания регулятора соответственно, оС.

Решение модели (1)-(21) позволяет моделировать процесс нагрева плиты не только в начальный период времени, но и в наиболее важный период регулярного (стационарного) нагрева.

Результаты решения модели (1)-(21) для вышеупомянутой промышленной плиты в системе COMSOL Multiphysicsи их сравнение с экспериментальными данными для температур контрольной (кривая 1) и центральной (кривая 2) термопар в начальный период нагрева показаны на рис. 6 (экспериментальные значения для контрольной термопары представлены в виде прямоугольных маркеров, для центральной -  в виде треугольных маркеров).

 

 

Рис. 6 – Температура контрольной (1) и центральной (2) термопар в начальный период, °С

 

Из данных рис. 6 следует, что максимальное отклонение расчётных данных от экспериментальных к моменту времени 180 с для контрольной термопары составило 6,1 °С, а для центральной термопары 5,3 °С. Средняя ошибка за весь рассматриваемый период 0-1968 с для контрольной термопары составила 6,6 %, для центральной – 6,9 %. Отметим, что по мере прогрева плиты, начиная с момента времени 1470 с, отклонение температур не превышает 2 % для контрольной термопары и 3 % для центральной. Таким образом, значения расчётных температур обеих термопар удовлетворительно согласуются с полученными в эксперименте.

На рис. 7 представлено сравнение расчётных температур контрольной термопары (кривая 2) с экспериментальными (кривая 1) в течение регулярного режима нагрева. Как видно, график расчётных температур отличается от графика экспериментальных температур только по фазе: отставание расчётной кривой от экспериментальной в среднем составляет 34 с. Для рассматриваемых интервалов времени данная ошибка является незначительной. Можно сделать вывод, что результаты численного решения модели (1)-(21) и в данном случае удовлетворительно согласуются с данными эксперимента.

 

 

Рис. 7 – Температура контрольной термопары в регулярном режиме, °С
        1 – эксперимент, 2 – расчёт

 

3. Математическая модельпроцесса нагрева системы "индукционная плита–пресс-форма–резинотехническое изделие"

Результаты решения моделей (1)-(20), (1)-(21) и их сравнения с экспериментальными данными позволяют перейти к этапу моделирования процесса нагрева многоэтажного пресса. Для этого необходимо ввести следующие допущения.

1) Между всеми элементами пресс-форм и изделиями, а также между элементами пресс-форм и нагревательными плитами имеет место идеальный тепловой контакт.

2) Тепловые эффекты вулканизации РТИ незначительны и ими можно пренебречь (справедливость этого допущения для изделий небольшого объема: прокладок, мембран, пробок и т.п. подтверждена расчетами по методике [12]).

Модель (1)-(21) необходимо дополнить уравнениями распространения тепла в  пресс-формах и изделиях

,                  (22)

,        (23)

где Tm=T(xm,ym,zm), (xm,ym,zm)Ωm – температурное поле пресс-форм, оС;

      nmolds – суммарное количество пресс-форм на всех этажах пресса;

      Trp=T(xrp,yrp,zrp), (xrp,yrp,zrp)Ωrp – температурное поле в изделиях, оС;

       λm, сm, ρm – теплопроводность (Вт/м·К), теплоемкость (Дж/кг·К) и плотность (кг/м3) материала пресс-форм, соответственно;

       λrp, сrp, ρrp – теплопроводность (Вт/м·К), теплоемкость (Дж/кг·К) и плотность (кг/м3) материала изделий, соответственно;

      Ωm, Ωrp– объём пресс-форм и РТИ, соответственно;

      nproductss – общее количество изделий во всех пресс-формах.

В качестве граничных условий для уравнений (22) и (23) используются следующие соотношения.

– Равенство температур и тепловых потоков на границе между пресс-формами и РТИ

.     (24)

Здесь Tm.b=T(xm.b,ym.b,zm.b)температурное поле элементов пресс-форм на границе с изделиями;

           Trp. =T(xrp.b,yrp.b,zrp.b) – температурное поле изделий на границе с элементами пресс-форм;

         (xm.b,ym.b,zm.b)Ωm, (xrp.b,yrp.b,zrp.b)Ωrp множество точек, принадлежащих элементам пресс-форм и изделий на границе их раздела.

– Равенство температур и тепловых потоков между рабочими поверхностями нагревательных плит и матрицами пресс-форм:

(25)

Здесь Tm.b1=T(xm.b1,ym.b1,zm.b1)температурное поле матриц пресс-форм на границе с плитами;

            Tpl.b2 =T(xpl.b2,ypl.b2,zpl.b2) – температурное поле плит на границе с матрицами пресс-форм;

          (xm.b1,ym.b1,zm.b1)Ωm, (xpl.b2,ypl.b2,zpl.b2)Ωplb2 множество точек, принадлежащих матрицам пресс-форм и нагревательных плит на границе их раздела.

– Равенство температур и тепловых потоков между крышками нагревательных плит и пуансонами пресс-форм:

.  (26)

Здесь Tm.b2=T(xm.b2,ym.b2,zm.b2)температурное поле пуансонов пресс-форм на границе с крышками;

           Tpl.b1 =T(xpl.b1,ypl.b1,zpl.b1) – температурное поле крышек плит на границе с пуансонами пресс-форм;

        (xm.b2,ym.b2,zm.b2)Ωp, (xpl.b1,ypl.b1,zpl.b1)Ωplb1 множество точек, принадлежащих пуансонам пресс-форм и крышкам нагревательных плит на границе их раздела.

– Конвективный (qm.c) и лучистый (qm.r) тепловой поток от внешних поверхностей пресс-форм в окружающий воздух: 

,                                                   (27)

.                                                    (28)

Здесь αm.c– коэффициент стеснённой теплоотдачи от поверхностей пресс-форм в окружающий воздух, Вт/м2∙К (см. [9]);

       Tm.s=T(xm.s,ym.s,zm.s) – средняя температура внешних поверхностей пресс-формы, ºК;

        εm – степень черноты поверхностей пресс-форм;

        σ = 5,67∙10-8 - постоянная Стефана-Больцмана, Вт/м24.

Таким образом, математическую модель нагрева многоэтажного пресса можно представить в виде системы уравнений (1)-(28). К числу ее параметров, кроме параметров конструкции нагревательных плит (см. раздел 1) относятся координаты расположения пресс-форм на поверхности плит. В качестве метода решения модели предлагается использовать МКЭ (система COMSOL Multiphysics).

В качестве примера использования математической модели (1)-(29) приведем результаты ее решения для промышленной плиты 432.383.35, спроектированной и применяемой на ЗАО «Завод Тамбовполимермаш» при производстве резиновых прокладок. Размеры плиты 600х600 мм, начальная мощность 6,5 кВт, значения параметров конструкции плиты введены как константы. На рис. 8а показана геометрическая модель системы, на рис. 8б – расчетное температурное поле продольного сечения прокладки в момент наибольшего перепада температур (5,35 °С – через 5326 с от начала нагрева), на рис. 8в  – графики изменения минимальной, средней и максимальной температур изделия. Изменение параметров конструкции нагревательных плит в пределах ограничений (14)-(20), расположения пресс-форм на поверхностях плит, может привести к уменьшению перепада температур в изделиях.

Как видно, методика теплового расчёта систем "индукционная нагревательная плита–пресс-форма–РТИ", основанная на использовании модели (1)-(28)  и системы конечноэлементного анализа COMSOL Multiphysics,  позволяет анализировать температурное поле РТИ в процессе вулканизации, определять влияние параметров конструкции индукционных нагревательных плит на степень его равномерности, формулировать и решать задачи оптимизации параметров плиты. На основе этой методики в настоящее время разрабатывается система автоматизированного проектирования прессового оборудования для ЗАО "Завод Тамбовполимермаш" и ОАО "АРТИ-завод", г. Тамбов.

 

assembly-numbered-article

а)

1,7 – основания плит; 2,6 – крышки плит; 3,9 – индукторы; 4,5 – матрица и пуансон пресс-формы; 6 – РТИ (прокладка)

 

б)

в)

Рис 8 – Результаты решения модели (1)-(28)
а) - геометрическая модель системы; б) - температурное поле продольного сечения прокладки,  τ = 5326 с, °С ; в) – максимальная (1), средняя (2) и минимальная температура РТИ (3) °С.

  

Заключение

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

Предложена математическая модель процесса индукционного нагрева системы "индукционная нагревательная плита–пресс-форма–резинотехническое изделие", применение которой в совокупности с системой COMSOL Multiphysics позволяет формулировать и решать задачи оптимизации параметров прессового оборудования.

 

Список литературы

1. Пресс гидравлический вулканизационный 250-600 (1Э, 2Э, 4Э). Режим доступа: http://www.tambovpolimer.ru/katalog_id/33/   (дата обращения  22.12.2012).

2. Дмитриев В.Г., Мананников П.Н. Тепловая модель вулканизационного пресса //

 Горный информационно-аналитический бюллетень (научно-технический журнал) . 2002. № 9. С. 182-184.

3. Multiphysics Modeling and Simulation Software - COMSOL. Режимдоступа: http://www.comsol.com/ (датаобращения  22.12.2012).

4. Карташов Э.М. Аналитические методы в теории теплопроводности твердых тел. М.: Высшая школа, 1985. 480 с.

5. Туголуков Е.Н. Решение задач теплопроводности методом конечных интегральных преобразований: учеб. пособие. Тамбов: Изд-во Тамб. гос. техн. ун-та, 2005. 116 с.

6. Кошляков Н.С., Глинер Э.Б., Смирнов М.М. Уравнения в частных производных математической физики.  М.: Высшая школа, 1970. 712 с.

7. Кувалдин А.Б. Индукционный нагрев ферромагнитной стали. М.: Энергоатомиздат, 1988. 200 с.

8. Сайт Морозовского химического завода. Органосиликатная композиция ОС-82-05. Режим доступа: http://tdmhz.ru/archives/1217/os-82-05/ (дата обращения  22.12.2012).

9. Романков П.Г., Фролов В.Ф., Флисюк О.М., Курочкина М.И. Методы расчета процессов и аппаратов химической технологии. СПб: Химия, 1998. 496 с.

10.Малыгин Е.Н., Карпушкин С.В., Крушатин А.С. Методика теплового расчёта нагревательных плит прессов для изготовления резинотехнических изделий // Химическая промышленность сегодня. 2009.  № 11. С. 48-56.

11.Карпушкин С.В., Карпов С.В., Глебов A.O. Моделирование устройств индукционного нагрева на примере индукционных нагревательных плит вулканизационных прессов (на английском языке) // Вестник Тамбовского ГТУ. 2011. Т. 17, № 1. С. 110-120.

12. Кузнецов Г.В., Иванова Е.В. Математическое моделирование температурных полей в процессе вулканизации типичных кабельных изделий // Известия Томского политехнического университета. 2010. Т. 316, № 4. С. 38-41.

Поделиться:
 
ПОИСК
 
elibrary crossref ulrichsweb neicon rusycon
 
ЮБИЛЕИ
ФОТОРЕПОРТАЖИ
 
СОБЫТИЯ
 
НОВОСТНАЯ ЛЕНТА



Авторы
Пресс-релизы
Библиотека
Конференции
Выставки
О проекте
Rambler's Top100
Телефон: +7 (915) 336-07-65 (строго: среда; пятница c 11-00 до 17-00)
  RSS
© 2003-2024 «Наука и образование»
Перепечатка материалов журнала без согласования с редакцией запрещена
 Тел.: +7 (915) 336-07-65 (строго: среда; пятница c 11-00 до 17-00)