
- •5 Математическое моделирование процессов переработки пластмасс
- •Оглавление
- •Список использованных обозначений
- •Введение
- •1. Основы деформации и течения расплавов полимеров
- •1.1. Понятие деформации
- •1.1.1. Деформации объемного сжатия или расширения
- •1.1.2. Упругая и пластическая деформация.
- •1.2. Понятие о скорости сдвига
- •1.2. Деформация при течении полимеров
- •1.3. Ньютоновские жидкости.
- •1.4. Вязкопластичные жидкости (тело Шведова — Бингама).
- •1.5. Дилатантные жидкости.
- •1.6. Псевдопластичные (псевдовязкие) жидкости.
- •2. Механические модели
- •2.1. Простейшие механические модели
- •2.3. Обобщенная модель Максвелла.
- •2.4. Модель Кельвина — Фойгта.
- •2.5. Обобщенная модель Кельвина — Фойгта.
- •2.6.Модель Алфрея-Гарни (Бургерса-Френкеля)
- •3. Явления переноса
- •3.1. Уравнение неразрывности
- •3.2. Уравнение движения
- •3.3. Уравнение энергии
- •3.4. Тензор напряжений
- •3.5. Тензор скоростей деформации
- •4. Общие граничащие условия и упрощающие предположения
- •5. Изотермическое установившееся течение аномально вязких жидкостей в цилиндрическом канале
- •6. Изотермическое установившееся течение аномально вязких жидкостей между вращающимися цилиндрами
- •7. Изотермическое установившееся течение аномально вязких жидкостей между параллельными пластинами
- •8. Капиллярная визкозиметрия. Входовые эффекты.
- •11. Специфические эффекты при течении полимеров
- •11.1. Развитие нормальных напряжений.
- •11.2. Эффект Барруса.
- •11.3. Течении расплавов полимеров в кольцевых каналах
- •11.4. Эффект Вайссенберга.
- •Список использованных источников
8. Капиллярная визкозиметрия. Входовые эффекты.
Рассмотрим движение расплава в цилиндрическом канале сравнительно небольшой длины с радиусом R. Расплав поступает из канала большого сечения и течет под действием перепада давления, создаваемого по длине канала (рис. ).
При переходе расплава из большого резервуара в канал небольшого сечения происходит сравнительно быстрое изменение скорости, которая стабилизируется на некотором удалении от входа. Отдельные авторы считают, что вначале профиль скорости имеет почти прямоугольную форму, однако эксперименты, выполненные на прозрачных каналах при течении с контрастным наполнителем, этого не подтверждают. Формование профиля начинается перед входом за счет образования в расплаве конического гидродинамического конфузора, как показано на рис. 2.14. Центральная струя имеет наибольшую скорость, которая уменьшается по мере удаления от оси канала и, если канал имеет прямоугольный вход, то в углах могут образоваться циркуляционные потоки (вихревые зоны).
Как уже было рассмотрено ранее, под действием напряжений сдвига в расплавах полимеров может происходить разрушение структуры и изменение конформации макромолекул, т.е. течение сопровождается развитием упругой деформации. Поскольку на входе в канал скорость изменяется по длине, т.е. имеется градиент скорости (dvz/dz)0, то упругая деформация постоянно возрастает. При увеличении скорости течения расплава полимера повышается степень разрушения структуры и уменьшается число пространственных узлов.
Изменение давления по длине канала при неустановившемся течении расплава: неустановившееся течение dp/dz const; установившееся течение dp/dz = const.
Рис. .
Поскольку все эти процессы протекают не мгновенно, то завершаются они с некоторым отставанием от темпа развития скорости. Следовательно, по мере течения расплава по длине канала изменяется напряжение сдвига и меняется доля энергии, затрачиваемая на развитие упругой деформации и разрушение структурных элементов системы. На входе в канал затрата энергии максимальна, а затем она постепенно уменьшается и при переходе к установившемуся режиму энергия расходуется лишь на развитие вязкого течения. Часть затраченной на входе энергии переходит в теплоту, а другая часть в виде потенциальной энергии (накопленная упругая деформация) сохраняется в расплаве и реализуется на выходе из канала при разбухании струи. Такое неравномерное потребление энергии вызывает непропорциональные потери давления по длине канала.
Таким образом, канал, в зависимости от режима течения, можно разбить на два участка. Первый — участок входа, когда (dp/dz)const, на котором потери давления зависят от трех составляющих:
1) неустановившееся течение в сходящемся потоке на входе в канал,
2) вязкие потери при перестройке профиля скорости потока,
3) развитие упругих деформаций.
Второй — участок с установившимся течением при (dp/dz)=const, где перепад давления обусловлен преимущественно вязкими потерями при течении расплава. Если решение для неустановившегося течения выполняется с учетом входовых потерь и проявления упругих деформаций, то в реологическое уравнение необходимо ввести дополнительные компоненты, описывающие вязкоупругие свойства и инерционные силы. То есть в уравнениях движения необходимо учесть целый ряд дополнительных членов: dvz/dz, dzz/dz и dp/dz = f(z), что очень затрудняет решение.
Для инженерных расчетов целесообразнее использовать косвенныйграфо-аналитический метод, для чего процесс течения рассматриваемкак установившийся, а непропорциональное изменение давления заменяем линейной функцией с введением поправочного входового коэффициента рвх или lвх (см. рис. ). В данном случае общий перепад давления заменяем суммой:
где dp/dz - градиент давления; l - длина канала; рвх - потери давления на входе; рl, ро - давление на входе и выходе из канала.
Таким приемом мы как бы искусственно все дополнительные затраты энергии на входном участке сосредотачиваем на входе в канал и сводим решение задачи к установившемуся режиму течения, когда
Для цилиндрического канала:
Для плоскощелевого канала:
где р = рl - р0
Приведение к установившемуся течению можно выполнить также вторым способом. Продолжим пропорциональную зависимость распределения давления влево от оси ординат до пересечения с линией давления, равного рг При этом мы искусственно удлиняем канал на величину lвх (см. рис. ), и течение заменяем установившимся. В этом случае напряжения сдвига будут равны:
Для определения скоростей сдвига пригодны уравнения (2.41) и (2.63).
Обычно для удобства расчетов вместо lвх используют входовой поправочный коэффициент т:
т = lвх /R или lвх =mR (lвх = ND – поправка Бэгли)
Поскольку входовые потери связаны с развитием вязкоупругих дополнительных деформаций и разрушением структуры, то они зависят от скорости сдвига и температуры.
Для экспериментального определения входовых потерь давления или поправочного коэффициента обычно используют каналы (сопла) различной длины, но одинакового диаметра. Измеряют последовательно зависимость давления от объемного расхода при течении в каналах различной длины и строят график (рис. 2.16). Затем находят давления на входе в канал при постоянных расходах для каждой длины (горизонтальные сечения). Определив р1, р2 и рз строят графическую зависимость изменения давления от длины канала. Если длина наиболее короткого канала больше, чем длина входового участка с неустановившимся течением, то получаем линейную зависимость (рис. 2.17). Интерполируя полученную линейную зависимость на ось ординат, находим значение давления pвх для различных расходов (скоростей сдвига). Для слишком коротких каналов линейной зависимости не получается, поэтому такой метод определения входных потерь использовать нельзя.
Если проинтерполировать полученную зависимость до пересечения с осью абсцисс, то отсекаемые отрезки будут численно равны поправочному
коэффициенту т =lвх / R.
Графическую зависимость, показанную на рис. 2.16, можно использовать также для расчета напряжения сдвига третьим способом. У каналов с одинаковыми радиусами и с постоянным объемным расходом расплава, скорость сдвига и входовые потери также будут постоянными. Поэтому для каналов различной длины, можем записать:
;
Решая совместно эти уравнения и считая, что 1 = 2 при Q1 = Q2, находим:
Зависимость давления от длины канала, используемая для нахождения входных потерь.
Рис. 2.17.
Зависимость объемного расхода от давления при различной длине канала (l2 > l1).
Рис. 2.16.
Таким образом, подставив значения давлений р1 и р2, найденные из графика рис. 2.16, мы можем найти изменение напряжения сдвига от скорости сдвига для различных расходов.
Полученные уравнения применяются для обработки данных капиллярной вискозиметрии, а также для нахождения перепада давления (обратная задача) при течении расплавов полимеров в формующих экструзионных головках или в литниковых каналах формы при литье под давлением. Применяются они в тех случаях, когда при переходе от одного канала к другому имеется резкое сужение канала или его резкий поворот.
9. Ротационная вискозиметрия
9.1. Вискозиметры конус - конус, конус - плоскость
9.2. Вискозиметры коаксиальные цилиндры
10. Практические методы количественного описания аномалии вязкого течения
Большинство расплавов обладает свойствами аномально-вязких жидкостей. Представим зависимость скорости сдвига от напряжения сдвига в обычных координатах. Если расплав обладает свойствами ньютоновской жидкости, то в этих координатах его кривая течения изображается прямой с угловым коэффициентом, равным 1/, где - ньютоновская вязкость. Если же расплав обладает свойствами аномально-вязкой жидкости, то его кривая течения выгнута по направлению к оси напряжений. Это означает, что каждой точке этой кривой соответствует свое значение эффективной вязкости, численно равное единице, деленной на угловой коэффициент прямой, проходящей через начало координат и соответствующую точку кривой.
Для полного суждения о реологических свойствах расплава надо представлять себе его кривую течения при изменении скорости сдвига в диапазоне 3—4 десятичных порядков. Поэтому экспериментальные данные представляют в логарифмических координатах. В том случае, если масштабы, выбранные для оси абсцисс и для оси ординат, одинаковы, системам, обладающим свойствами ньютоновских жидкостей, будут соответствовать прямые, наклоненные к осям под углом 45°. При этом абсолютная величина вязкости сказывается только на месте расположения прямой.
Кривые течения
1 – ньютоновская жидкость,
2 – аномально-вязкая жидкость
рис. .
Кривые течения реальных расплавов в логарифмических координатах слабо изогнуты. Обратим внимание на то, что для псевдопластичных жидкостей значение производной в любой точке этих кривых больше единицы. Это условие записывается в виде неравенства
Значение
производной позволяет легко оценивать
степень аномалии
вязкого течения. Выше мы отмечали, что
для ньютоновских жидкостей
производная равна единице. Чем больше
аномалия вязкости
расплава, тем больше значение производной.
Для расплавов термопластов
и эластомеров производная
в
пределах
одной кривой может изменяться от 1 до
5.
Слабая кривизна логарифмических кривых течения явилась причиной того, что наиболее широко для их математического описания применяется так называемый степенной закон, имеющий вид
или
где n' = 1/n.
Обычно экспериментатор располагает данными о зависимости напряжения сдвига от скорости сдвига при изменении последней в пределах 2—3 десятичных порядков. Если обследованная область достаточно далеко отстоит от области ньютоновского течения, то логарифмическая кривая течения легко аппроксимируется прямой:
Значение n, по определению, равно производной для соответствующего участка логарифмической кривой течения:
Поскольку величина n характеризует аномалию вязкости, она получила название индекса течения.
Для ньютоновской жидкости величина имеет смысл вязкости. Для аномально-вязких жидкостей величина не имеет столь четкого физического смысла. Она является своеобразным аналогом вязкости и обычно называется коэффициентом консистенции.
Известны и более сложные способы описания кривых течения. Трехпараметрическое степенное уравнение:
Здесь 0 - предельное значение вязкости, соответствующее области ньютоновского течения при минимальных напряжениях сдвига. В области малых напряжений сдвига уравнение вырождается в закон Ньютона и хорошо описывает начальный участок кривой течения.
В области высоких напряжений сдвига доминирующее значение приобретает второй член, и уравнение фактически превращается в обычный степенной закон. Наконец, в области переходных напряжений оба члена оказываются соизмеримыми, и уравнение удовлетворительно описывает участок с переменным индексом течения.
При необходимости увеличить точность аппроксимации кривой течения приходится прибегать к еще более сложным зависимостям.
В
качестве примера сошлемся на
четырехпараметрическое уравнение:
,
где 0 — ньютоновская вязкость, соответствующая начальному участку кривой течения; k1, k2 и n - эмпирические константы, определяемые обработкой экспериментальных данных.