МАШИНОСТРОЕНИЕ И МАШИНОВЕДЕНИЕ
УДК 532.517.2:534.2
Б. А. Снигерев
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ НЕИЗОТЕРМИЧЕСКОЙ ЭКСТРУЗИИ ПОЛИМЕРНОЙ ЖИДКОСТИ ИЗ СТУПЕНЧАТОЙ ФОРМУЮЩЕЙ НАСАДКИ
Аннотация. Работа посвящена моделированию медленного движения вязкоупругой жидкости со свободной поверхностью, реализующейся при входе полимерной жидкости в формующую насадку и выходе из нее. Движение жидкости описывается уравнениями сохранения массы, импульса и энергии, дополненными реологическим уравнением состояния среды Гиесекуза. На основе метода конечных элементов разработан устойчивый численный алгоритм решения задачи. Проведены численные исследования по определению формы выходной струи для различных режимов течения и формы насадки. Исследована картина распределения скорости жидкости, давления, напряжений и температуры при увеличении степени нагрева стенки насадки. Получены численные результаты зависимости эффекта разбухания полимера от параметров реологической модели и температурных факторов.
Ключевые слова: вязкоупругая жидкость, свободная поверхность.
Abstract. Numerical simulation flow of viscoelastic fluid with free surface, which is realized in entrance and output flow in extrusion die was performed. The flow of liquid is described by equations of conservation of mass, momentum and thermal energy with rheological constitutive equation of Giesekus. On basis of finite element method the stable numerical scheme was developed to solve this problem. Different numerical experiments was performed to define the configuration of outflow jet in various regimes and construction of die. The distribution of flow velocity fields, pressure and temperature are investigated on dependence of heating the walls. The ratio of extrusion in dependence of parameters the rheological model are investigated.
Key words: viscoelastic fluid, free surface.
Введение
Экструзия является одним из основных способов переработки резиновых смесей в изделия. В современных рыночных условиях на промышленных предприятиях резинотехнической, кабельной, шинной промышленности предъявляются высокие требования к качеству изделий, которое определяет во многом себестоимость продукции, а следовательно, ее конкурентоспособность. Одним из факторов, влияющих на качество изделия, является выбор технологических параметров процесса экструзии и их соблюдение. Теоретические и экспериментальные исследования процессов, проходящих в канале формующей головки экструдера при переработке резиновых смесей, в публикациях представлены весьма ограниченно [1, 2]. Для производства резино-
технических изделий высокого качества целесообразно предусмотреть контроль и автоматическое управление процессом экструзии. При разработке автоматических систем управления процессом экструзии на экструзионных машинах необходимо иметь адекватные математические модели, описывающие стационарные и нестационарные режимы ее работы. К настоящему времени число работ, посвященных исследованию гидродинамических процессов в канале формующей головки и экструзии из формующей насадки, незначительно. Работы по исследованию гидродинамических процессов в экструзионных машинах, перерабатывающих резиновые смеси, ограничены экспериментальными исследованиями, проведенными на лабораторных машинах.
1. Математическая постановка задачи
В данной работе рассматривается вытекание вязкоупругой жидкости из цилиндрической ступенчатой щелевой насадки, схема которой представлена на рис. 1 (И обозначает радиус узкого канала выходной части насадки; 4И -
радиус трубы на входе; И1 =0,5И - длина закругленной части; Ь - длина выходной части насадки). Рассматривается влияние длины выходной части насадки на степень разбухания полимерной жидкости в процессе экструзии. Выбираются значения Ь, равные Ь =0,2И,4И,10И . Анализируется также влияние реологических параметров на характер истечения из насадки.
Рис. 1. Схема вытекания вязкоупругой жидкости из ступенчатой цилиндрической насадки
Введем следующие обозначения границ: £1 - входное сечение; £2 -
выходное сечение струи; £3 - часть границы, являющейся твердой стенкой;
£4 - граница, являющаяся осью симметрии области; £5 - часть границы,
приходящаяся на свободную поверхность.
Исходной является система уравнений законов сохранения массы, импульса и энергии, для замыкания которой привлекается конститутивное реологическое соотношение Гиесекуза [2, 3]:
Эу Э V,
—L+V,—-
дґ 1 дх,
V ■)
ЭР Этц Эу,
д х,
+
д х/
= 0,
V
дг
=г/ + 2 ЛмВ,
§/_ 8t дt
У" к дг" д у
+ V —г---------------~гк,- ~гк -
д хк д х- к/ 1к д
8г" «А
V т„ 8 г/ «А V 2 - г, г, 1
г" + +^Т% =2 'Л"°->; = 2
/
ду ду/
- +- J
д х/ V J
д X ;
(
Р Ср
дТ
дТ
\
+ у ,■
V д ^ д X у
= к
д2 Т д 2 х
(ггк ' Вк/'):
(1)
где у - компоненты скорости жидкости; х-, X/ - координаты (, / = 1, 2); t - время; р - плотность; Р - давление; г-/ - девиатор напряжения; г"/ -вязкоупругая часть напряжения; Т - температура жидкости; Л = Л " + Л^ ; Л", Л^ _ динамические вязкости при нулевой скорости сдвига полимера и растворителя соответственно; А - время релаксации напряжений; к - коэффициент теплопроводности; ер - удельная теплоемкость полимера при постоянном давлении; « - безразмерный параметр модели Гиесекуза, определяющий вязкоупругие свойства полимерной жидкости.
Для зависимости вязкости и времени релаксации от температуры используется соотношение Аррениуса [1, 2]:
Л(Т) = Л0 а(Т), А(Т)=Ао а(Т), а(Т)=ехр
Е
Я (1/Т -1/Т))
где индекс 0 означает, что значение параметра берется при температуре Т0 ; Е - энергия активации; Я - универсальная газовая постоянная.
Граничные условия имеют вид
£1:у2 = у2 (х2 ), г" = г" (х2 ), Т =Т0 (х2 ) ;
£2 У = 0, ду / дх1 =0, д Т /д х1 = 0 ; £3: у = 0, Т = Т2 ;
£4:^2 = 0, г"> = 0, дТ/д х2 = 0; £5: ёх- /dt = у,
П г-/П/ - (Р-Р0)п = (1/ор )2К, п гijti = 0, кдТ/ дп=-«р (Т-Т0),
где п, t^ - компоненты единичной нормали и касательной к поверхности £5 ; Р0 - внешнее давление над свободной поверхностью; О р - коэффициент поверхностного натяжения; «р - коэффициент теплоотдачи в окружающую среду; К - кривизна свободной поверхности; Т2 - температура стенок канала; Т1 - температура жидкости на входе; Т0 - температура окружающей среды. Исходная система уравнений (1) преобразуется в безразмерный вид. Для этого в качестве определяющих параметров приняты параметры: и0 -средняя скорость жидкости на входе; Ц - характерный размер, равный полуширине канала И; давление и напряжения отнесены к значению ри0/Ц.
Определяющими являются безразмерные параметры: Яе = риЬ/Л - число Рейнольдса, We=Аи/Ь - число Вайсенберга; Ре = р Ср и Ь/к - число Пекле,
Бг=Ли /(Т1 -70) - число Бринкмана; Са =Ли0/Ор - число капиллярности; Р=Л^ / Л - параметр ретардации; Б1 = «р Ь / к - параметр Био.
Аппроксимация уравнений (1) и вычисления проводятся методом конечных элементов (МКЭ) второго порядка на нерегулярных сетках, сгущающихся в зоне истечения полимера из насадки. Для расчетов строилась последовательность сгущающихся сеток 9-узловых четырехугольных элементов (число узлов 5000, 18400). Для расчета напряжения использовались линейные четырехугольные элементы. Местоположение деформируемой свободной поверхности находится из аппроксимации кинематического условия, затем сетка конечных элементов вблизи нее перестраивается для получения решений уравнений (1), с помощью которых находится поле скоростей, давлений, напряжений и температур на новом временном слое. Стационарное положение формы выходящей струи находится методом установления эволюционной задачи с использованием традиционных для уравнений данного класса алгоритмов [4].
2. Результаты исследования
Для изотермического обтекания значения всех параметров рассматриваются при температуре Т) = 303 К. В случае неизотермического обтекания рассматривается истечение полимерной жидкости из насадки с нагретыми стенками, температура которой повышается до Т = 323, 353 К соответственно. Расчеты проводились при следующих значениях теплофизических свойств полимерной жидкости: р= 924кг/м3, Л= 104 Па• с ,
с р = 2300 Дж/(кг • К), к = 0,26 Вт/(м • К), А = 2000 К, ор = 3 10-2 н/м,
«р = 102 Вт/(м2 • К) и Я = 10-2 м, Н = 10-2 м, и0 = 4-10-3 м/сек. Время релаксации напряжений жидкости изменялось в пределах от А = 0,1 до 10 с. Получены следующие значения безразмерных критериев: Яе = 10 5, а = 0,33,
Р = 0,11, Ре = 180, Бг = 2-10-2, Са = 1,4 • 103, Б1 = 3,8, а число We изменяется в пределах от 0,1 до 5,0.
На рис. 2 представлена картина течения для насадки с размером выходного канала Ь = 0 в виде линий тока (рис. 2,а), изолиний первой разности
главных напряжений N1 (рис. 2,б) и компоненты осевой скорости (рис. 2,в)
для изотермической экструзии при We = 3,0.
На рис. 2 видно, что, как только полимерная жидкость проходит выходное сечение и попадает в зону свободного течения, давление и напряжения начинают немедленно релаксировать. Процесс релаксации происходит до тех пор, пока накопленная высокоэластичная деформация не уменьшится до значения, соответствующего эластической деформации полимера, находящегося в состоянии стационарного течения. Выделим две основные причины изменения сечения струи:
- перестройка профиля скоростей от параболического, соответствующего течению в канале, к прямоугольному, соответствующему движению струи как сплошного твердого стержня;
- релаксация высокоэластических деформаций растяжения, накопленных в материале струи за время ее прохождения по каналу. Перестройка профиля скорости приводит к возникновению продольных деформаций растяжения или сжатия.
Рис. 2. Линии тока у (а), изолинии Л (б), изолинии осевой скорости (в) для изотермический экструзии при Ше = 3,0 для насадки с длиной выходной части Ь=0
Из требования выполнения закона сохранения количества движения следует, что струя, выходящая из канала, должна сужаться. Этот эффект известен для низкомолекулярных жидкостей, не проявляющих высокоэластических свойств. Для полимерных систем реализуется упругое разбухание, обусловленное высокоэластическими деформациями. Увеличение диаметра выходной струи характеризуется параметром И^, равным отношению диаметра
выходной струи к диаметру канала и называемым степенью разбухания. На рис. 3 приведена зависимость степени разбухания И^ выходной струи полимерной жидкости от значения параметра Ше для различных насадок. Формы насадок различаются длиной узкой выходной части Ь , изменяющейся от 0 до 10 И. Результаты показывают, что степень разбухания И^ струи увеличивается при уменьшении длины выходного участка. Видно, что значения И^ зависят от кинематики предшествующего деформирования - скорости сдвига
1,36
и длины выходной части капилляра. Это можно объяснить тем, что экструзионное разбухание при переработке полимеров обусловлено высвобождением упругой энергии, запасенной при предшествующем течении в канале.
И
Рис. 3. Степень разбухания выходной струи hf в зависимости от числа We для различных насадок: 1 - Ь=10h ; 2 - Ь=^; 3 - Ь=2h; 4 - Ь=0
Вязкоупругие жидкости являются средами, обладающими наследственной памятью. Наследственная жидкость «помнит предысторию», т.е. изменение поля скоростей в предшествующие моменты времени. Длительность «памяти» характеризуется временем релаксационных процессов X, что определяется безразмерным числом Вайсенберга We. Видно, что при увеличении значения числа We степень разбухания увеличивается для всех насадок. На рис. 4 показана зависимость hf от числа We при неизотермической экструзии для ступенчатой формующей насадки с размером выходной части Ь=2h в зависимости от температурного напора между стенкой и потоком полимера. Видно, что нагрев стенки на А Т = 30 К (кривая 3) приводит к уменьшению диаметра выходной струи. Отметим при этом, что для более длинной насадки с ростом числа We и температурного напора наблюдается немонотонный характер увеличения hf. При небольших числах We (We < 1) происходит уменьшение степени разбухания.
Заключение
Приведенные результаты расчетов процесса экструзии вязкоупругой жидкости из ступенчатых формующих насадок с разной длиной выходной части показали, что степень разбухания выходной струи увеличивается для коротких насадок. Полученные результаты подтверждают вывод о том, что одним из факторов, влияющих на степень разбухания, для наследственных вязкоупругих сред является время релаксации напряжений. Неизотермич-ность потока полимерной массы, определяющаяся разностью температур стенки и основного потока, приводит к уменьшению диаметра выходной струи. Выявлено, что при увеличении температурного напора наблюдается уменьшение диаметра выходной струи с ростом числа We на начальном этапе.
Рис. 4. Степень разбухания выходной струи для неизотермической экструзии в зависимости от числа We для насадки Ь = 0 при увеличении температуры стенки £3 : 1 - ДТ = 0 К; 2 - 15 К; 3- 30 К
Данный характер немонотонного поведения степени разбухания свидетельствует о том, что температурная аномалия вязкости оказывает существенное влияние на форму выходящей струи.
Список литературы
1. Торнер, Р. В. Теоретические основы переработки полимеров / Р. В. Торнер. -М. : Химия, 1977. - 467 с.
2. Раувендаль, К. Экструзия полимеров / К. Раувендаль. - СПб. : Профессия, 2008. - 768 с.
3. Giesekus, H. A simple constitutive equation for polymeric liquids based on the concept of deformation dependent tensorial mobility / H. Giesekus // J. Non-Newtonian Fluid Mech. - 1982. - V. 11. - P. 69-109.
4. Снигерев, Б. А. Течение упруговязкой жидкости со свободной поверхностью / Б. А. Снигерев, Ф. Х. Тазюков, А. Г. Кутузов, Амер аль Раваш // Вестник Казанского технологического университета. - 2007. - № 1. - С. 85-93.
Снигерев Борис Александрович
кандидат технических наук, старший научный сотрудник, Учреждение Российской академии наук «Институт механики и машиностроения Казанского научного центра РАН»
E-mail: Snigerev@mail.knc.ru
Snigerev Boris Alexandrovich
Candidate of engineering sciences, senior staff scientist, Institute of Mechanics and Engineering, Kazan Science Center, Russian Academy of Science
УДК 532.517.2:534.2 Снигерев, Б. А.
Математическое моделирование неизотермической экструзии полимерной жидкости из ступенчатой формующей насадки / Б. А. Снигерев // Известия высших учебных заведений. Поволжский регион. Технические науки. - 2011. - № 1 (17). - С. 196-202.