Тысяча полезных мелочей    

Математическая модель неизотермического течения вязкой жидкости В предматричной зоне экструдера

Математическая модель неизотермического течения вязкой жидкости В предматричной зоне экструдера

Д.т.н. Остриков А.Н. к.т.н. Абрамов О.В. к.т.н. Статья в формате PDF 129 KB

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

Для исследования движения потока расплава через отверстия матрицы экструдера рассмотрим течение жидкости в ступенчатом цилиндрическом канале, который на длине l1 будет иметь диаметр d1, а на длине l2 - диаметр d2 (d1 > d2). Принимаем, что поток расплава при подходе к предматричной зоне уже установился и течение расплава можно считать стационарным.

Для каждого k - го канала можно записать систему дифференциальных уравнений неизотермического течения сплошной среды. Рассмотрим асимметричное движение среды, которое моделируется двухмерным течением. В этом случае тангенциальная составляющая скорости равна нулю. Так как течение расплава является стационарным, а сплошная среда – аномально–вязкой и несжимаемой, то пренебрегая массовыми силами, которые малы вследствие достаточно большой вязкости, можно использовать следующие дифференциальные уравнения в цилиндрической системе координат: уравнение несжимаемости, уравнения движения и уравнение энергии с использованием функции диссипации Ф, хаpaктеризующей интенсивность преобразования кинетической энергии в тепловую. При изучении течения сплошной среды будем применять реологическое уравнение, которое выражается в виде обобщенного степенного закона. Выразим обозначенные уравнения в переменных тока φ и вихря ω, для которых уравнения связи с составляющими скорости θz, θr имеют вид

  

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

В результате преобразований получим в безразмерном виде:

-уравнение для функции тока (общее уравнение движения)

где безразмерные величины

уравнение энергии

  

уравнение для вихря

граничные условия

 

Таким  образом, получена математическая модель (2 6) неизотермического течения вязкой жидкости в цилиндрическом канале. В выражениях  принято:  - число  Рейнольдса - число   Эккерта;- число Пекле

В уравнение (2) входит нелинейная функция ω содержащая члeны с производными первого и второго порядка от функции вязкости η, хаpaктеризующей    вязкостные свойства  расплава. 

Структура функции ω может быть подобрана путем моделирования течения расплава в канале использованием модели (2 -6).

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

Для моделирования процесса течения расплава вязкой жидкости разработана программа Model на языке программирования Turbo Pascal .0 в среде Windows 95.

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

Расчет проведен для расплава вязкой жидкости в канале с размерами L = 12 мм, D = 12 мм и отверстиями в матрице d = 4 мм. Общее число узлов разностной сетки области D - NzR = 1875. Число узлов области D по координате z -75, по координате R - 25. Шаг сетки по координате z  составил Δz=0,027, по координате r- ΔR = 0,0417.

Расход расплава на входе в канал Rashs = 0,0905 м3/c и через отверстие в матрице Rashs = 0,0926  м3/c, т.е. ошибка на установление процесса составляет Δ = [(0,0905 - 0,0926)/0,0905] 100 % =2,3 %.

Эпюра скоростей расплава вязкой жидкости в различных сечениях вблизи отверстия в матрице приведена на рисунке.  

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

Выводы

  1. Выявлен хаpaктер изменения поля скоростей расплава вязкой жидкости в различных сечениях вблизи отверстия в матрице экструдера.
  2. Полученная модель позволяет с достаточной точностью ( 14 %) рассчитать не только скорости течения вязкой жидкости в предматричной зоне, но и определить производительность экструдера при изменении технологических параметров процесса экструдирования.

Рис. Эпюра скоростей расплава на входе в канал матрицы в различных сечениях предматричной зоны (r/R):

1 - 2,00; 2 - 1,86; 3 - 1,59; 4 - 1,32; 5 - 1,05

 



О НЕКОТОРЫХ ОСОБЕННОСТЯХ ГОСУДАРСТВЕННОГО УЧАСТИЯ В ЭКОНОМИКЕ (В СВЕТЕ ПОЛОЖЕНИЙ ИНСТИТУЦИОНАЛЬНОЙ ЭКОНОМИКИ)

О НЕКОТОРЫХ ОСОБЕННОСТЯХ ГОСУДАРСТВЕННОГО УЧАСТИЯ В ЭКОНОМИКЕ (В СВЕТЕ ПОЛОЖЕНИЙ ИНСТИТУЦИОНАЛЬНОЙ ЭКОНОМИКИ) В статье дана хаpaктеристика отдельных особенностей государственного участия в экономике. В частности, уделено внимание роли государства как гаранта существующей системы прав собственности, монополиста по производству общественных благ и хозяйствующего субъекта, стремящегося к максимизации собственных, неналоговых доходов. ...

09 08 2026 0:25:14

ДИПЛОМНОЕ ПРОЕКТИРОВАНИЕ (учебное пособие)

ДИПЛОМНОЕ ПРОЕКТИРОВАНИЕ (учебное пособие) Статья в формате PDF 97 KB...

08 08 2026 8:42:13

КОНСЕРВАТИВНОЕ ЛЕЧЕНИЕ ТРАВМ СЕЛЕЗЕНКИ

КОНСЕРВАТИВНОЕ ЛЕЧЕНИЕ ТРАВМ СЕЛЕЗЕНКИ Представлены данные литературы, посвященные изучению консервативной тактике при травматических повреждениях селезенки. Показаны показания и противопоказания и необходимые условия для проведения консервативного лечения таких повреждений. ...

05 08 2026 14:33:30

О ПЕРСПЕКТИВАХ ЭВОЛЮЦИИ НООСФЕРЫ

О ПЕРСПЕКТИВАХ ЭВОЛЮЦИИ НООСФЕРЫ Статья в формате PDF 183 KB...

28 07 2026 4:48:57

О ФИЗИКЕ СЕЙСМИЧЕСКОГО ИЗЛУЧЕНИЯ

О ФИЗИКЕ СЕЙСМИЧЕСКОГО ИЗЛУЧЕНИЯ Рассматриваются процессы формирования и распространения сейсмического излучения на основе ньютоновской механики. В источниках излучения среда приобретает механический импульс, который распространяется в виде пакета, действующего на элементы среды с силой, равной производной импульса по времени передачи. ...

27 07 2026 11:57:44

МАТЕМАТИЧЕСКИЕ МОДЕЛИ В ЗДРАВООХРАНЕНИИ

МАТЕМАТИЧЕСКИЕ МОДЕЛИ В ЗДРАВООХРАНЕНИИ Статья в формате PDF 103 KB...

20 07 2026 9:12:49

БИОЭЛЕКТРИЧЕСКАЯ АКТИВНОСТЬ СЕРДЦА ДЕВУШЕК АЛТАЙСКОЙ И РУССКОЙ НАЦИОНАЛЬНОСТЕЙ РАЗНОГО СРОКА ПРОЖИВАНИЯ НА ТЕРРИТОРИИ РЕСПУБЛИКИ АЛТАЙ

БИОЭЛЕКТРИЧЕСКАЯ АКТИВНОСТЬ СЕРДЦА ДЕВУШЕК АЛТАЙСКОЙ И РУССКОЙ НАЦИОНАЛЬНОСТЕЙ РАЗНОГО СРОКА ПРОЖИВАНИЯ НА ТЕРРИТОРИИ РЕСПУБЛИКИ АЛТАЙ Целью исследования явился анализ биоэлектрической активности сердца коренных и пришлых дeвyшек Горного Алтая алтайской и русской национальностей по данным электрокардиографии. Выявлено, что длительность интервала QT снижена во всех исследуемых группах, а интервала ТР и комплекса QRS превышает общепринятые значения. Снижение длительности всех интервалов происходит от алтаек к русским пришлым, что может свидетельствовать о более выраженной симпатикотонии среди последних. Расчетные величины показывают существенное несоответствие фактических значений должным, за исключением синусового ритма, который также снижается от алтаек к русским пришлым. Анализ амплитудно-частотных хаpaктеристик указывает на нарушения, связанные с процессами реполяризации, внутрижелудочковой и внутрипредсердной проводимости, а также гипертрофии этих отделов. О гипертрофии отделов сердца и нарушениях внутрижелудочковой проводимости в виде блокады ножек пучка Гиса свидетельствует также положение электрической оси сердца. ...

07 07 2026 20:23:26

ОСНОВНЫЕ ЭТАПЫ РАЗВИТИЯ СЛЕПОЙ КИШКИ У БЕЛОЙ КРЫСЫ

ОСНОВНЫЕ ЭТАПЫ РАЗВИТИЯ СЛЕПОЙ КИШКИ У БЕЛОЙ КРЫСЫ Статья в формате PDF 243 KB...

06 07 2026 1:47:17

Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::