Моделирование полимеризации – основа регулирования свойств упаковочных полимерных материалов, используемых в пищевой промышленности

 

Название работы

Моделирование полимеризации - основа регулирования свойств упаковочных полимерных материалов, используемых в пищевой промышленности

Перечень используемых сокращений

ПММА - полиметилметакрилат

ММА - метилметакрилат

ПС - полистирол

ПЭТФ - полиэтилентрефталат

ПП - полипропилен

DEHP - фталаты

ПВХ - поливинилхлорид

ПВДХ - поливинилдихлорид

МЦ - металлоцены

ГЭ - гель-эффект

полимеризация ферроцен математический моделирование

Содержание

Введение

. Литературный обзор

.1 Современные материалы для упаковки

.2 Строение и свойства полиметилметакрилата

.3 Возможность осуществления синтеза полиметилметакрилата с заданными параметрами

.4 Заключение по проведенному литературному обзору

. Экспериментальная часть

.1 Материалы

.2 Методики эксперимента

.2.1 Проведение полимеризации в присутствии ферроцена

.2.2 Построение кинетических кривых

.2.3 Определение молекулярно-массовых характеристик полимера

.2.4 Методика осуществления математического моделирования

.3 Метрологическая обработка экспериментальных данных

. Результаты эксперимента и их обсуждение

Заключение

Выводы

Список используемой литературы

Введение

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

Популярные современные упаковочные материалы имеют существенный недостаток - для придания им необходимых потребителю свойств в них добавляют различные стабилизаторы, пластификаторы, смягчители, наполнители и т.д.

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

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

Другой путь решения представляет собой изменение самого метода задания полимеру необходимых для упаковки характеристик. То есть вместо использования добавок производить синтез в таких макрокинетических условиях, при которых полимер приобретает необходимые потребителю свойства.

Единственно возможной реализацией данного подхода до настоящего времени являлся эмпирический метод - это определение оптимальных параметров для синтеза полимерных материалов с желаемыми свойствами экспериментальным путем. Данный метод требует больших финансовых вложений и временных ресурсов.

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

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

Разработка математической модели осуществлена для системы метилметакрилат (ММА) - инициатор азобисизобутиронитрил - ферроцен.

На основе созданной модели будет сформирована методика прогнозирования набора молекулярно-массовых характеристик (ММХ) получаемого ПММА. Варьированием показателей ММХ можно осуществлять изменение и контроль за физико-механическими, оптическими и электретными свойствами синтезируемого полимера.

Научная новизна данной работы заключается в том, что в мире еще не было проведено исследований в области математического моделирования процесса полимеризации ММА в присутствии ферроцена.

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

Впервые проведены теоретические и экспериментальные исследования по определению влияния температуры, начальных концентраций мономера и ферроцена на молекулярно-массовые характеристики получаемого ПММА;

Впервые разработан математический формализм, который лег в основу системы управления характеристиками ПММА, синтезируемого полимеризацией в присутствии ферроцена.

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

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

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

Выявленные в процессе работы зависимости можно использовать для моделирования полимеризации других виниловых мономеров.

Таким образом, цель курсовой работы заключается в создании методологии получения полимера с заданными молекулярно-массовыми характеристиками путем разработки модульной математической модели на примере синтеза ПММА в присутствии ферроцена.

Для достижения цели в работе были решены следующие задачи:

разработка кинетической модели процесса радикальной полимеризации ММА с учетом гель-эффекта и решение обратной кинетической задачи полимеризации ММА - восстановление констант скоростей элементарных стадий из экспериментальных данных и кинетической модели;

проверка адекватности полученной кинетической модели радикальной полимеризации ММА, проводимая сравнением экспериментальных значений молекулярно-массовых характеристик ПММА со значениями, которые были выведены теоретически;

разработка кинетической модели процесса полимеризации ММА с добавлением ферроцена и учетом гель-эффекта и решение обратной кинетической задачи полимеризации ММА при добавлении ферроцена, что позволит определить константы скоростей элементарных стадий полимеризации в координационной сфере металла;

проверка адекватности полученной модели полимеризации ММА при добавлении ферроцена на экспериментальных данных при определенной температуре;

на основе систематизации полученных результатов осуществление создания единой методологии контроля над молекулярно-массовыми характеристиками (а как следует и свойствами) ПММА, синтезируемого в присутствии МЦ.

Автор выражает благодарность профессору кафедры ТППКМ Улитину Н.В., ассистенту кафедры ТППКМ Терещенко К.А. за консультацию по вопросам математического моделирования и профессору БашГУ Колесову С.В. за предоставленные экспериментальные данные.

1. Литературный обзор

.1 Современные материалы для упаковки

В современном мире высокой конкуренции среди производителей пищевых продуктов важным является то, как визуально выглядит готовый продукт [1]. Важность упаковки невозможно недооценить, ведь лишенный возможности анализа и дегустации потребитель оценивает продукт по внешнему виду [2, 3].

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

Одни требования предъявляются к упаковке в общем, другие - к конкретной продукции. Последние существенно различаются: к примеру, упаковка для пищевых продуктов <#"93" src="doc_zip1.jpg" />

Рис 2.1 - Структурная формула метилметакрилата

Для инициирования полимеризации в системе был использован 2,2'-азоизобутиронитрил.

МЦ, в частности, ферроцены, используемые в работе являются продуктами фирмы Aldrich, США (рисунок 2.2).

Рис 2.2 - Структурная формула ферроцена

Для осуществления модульного моделирования была использована программа Wolfram Mathematica 8. Выбор в пользу данной программы был сделан на основе простоты синтаксиса и интерфейса, функциональность, возможность реализации разнообразных численных методов и совокупность символьных, графических и численные вычисления, выполняемые в одном сеансе использования [17].

.2 Методики эксперимента

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

.2.1 Проведение полимеризации в присутствии ферроцена

Готовили раствор МЦ в мономере (ММА) с концентрацией 10 ммоль/л. Раствор помещали в ампулу из молибденового стекла и тщательно вакуумировали для удаления следов воздуха. Ампулу подвергали УФ-облучению (ртутной лампой) в течение 5 минут. Затем ампулу вскрывали, отбирали объем реакционной смеси, обеспечивающий получение концентрации [МЦ] = 1·10?3 моль/л мономера (например 0,2 мл смеси на 1,8 мл мономера) и вводили в мономер. Кинетику полимеризации исследовали дилатометрическим методом [18].

.2.2 Построение кинетических кривых

Кинетику процесса изучали дилатометрическим методом как на начальных (затворная жидкость - исходный мономер), так и на глубоких степенях превращения (затворная жидкость - глицерин) [19]. Для полимеризации в массе реакционную смесь (мономер, инициирующая система) заливали в ампулу и дегазировали раствор трехкратным повторением циклов замораживание - размораживание до остаточного давления 1,3 Па. Затем ампулу запаивали, помещали в термостат с точностью поддержания температуры ± 0.1 °С и выдерживали при выбранной температуре до достижения необходимой степени конверсии, которую рассчитывали из соотношения:

,

где V0 - первоначальный объем мономера,

DV - изменение объема мономера, - коэффициент контракции, равный:

,

где Vм и Vп - удельные объёмы мономера и полимера.

Значения k взяты из [20,21].

Начальную скорость полимеризации W0 рассчитывали по формуле:

,

где М - молекулярная масса мономера,

Dt - время, за которое объём изменился на DV,

rм - плотность мономера при данной температуре.

Кинетические кривые строили в координатах U (%) - Время (мин, час). Кинетические кривые обрабатывали в рамках пакета программ Microsoft Excel. Ошибка определения скорости полимеризации на любом участке кинетической кривой не превышает 3-5%. Ампулы после полимеризации охлаждали и вскрывали. Образовавшийся полимерный продукт растворяли в ацетоне и осаждали 10-15-кратным избытком этанола. От остатков инициатора и ферроцена, имея ввиду хорошую растворимость последних в ацетоне и хлороформе, полимеры очищали путем трехкратного переосаждения. Очищенные полимеры сушили в вакууме до постоянной массы при температуре 40 °С.

.2.3 Молекулярно-массовые характеристики полимеров

Для определения Mw, Mn и Mw/Mn полимеров использовали метод гель-проникающей хроматографии. Анализы проводили на жидкостном хроматографе Waters GPC 2000 System (элюент: хлороформ, скорость потока: 0.5 мл/мин). Систему колонок калибровали по полистирольным стандартам с Mw/Mn £ 1.2.

Суммарную эффективность хроматографических систем оценивали числом теоретических тарелок N (т.т.) по пику 0.3% (вес.) раствора гексана в толуоле.

Общая эффективность хроматографической системы составляла УN = 6500ч7000 т.т.

Полученная эффективность системы удовлетворяет условию:

,

где C1 и C2 - коэффициенты в уравнении VR = C1 - C2lgM,

VR - элюционный объем.

При выполнении данного условия молекулярно-массовые характеристики, найденные из хроматограмм, отличаются от результатов, полученных абсолютными методами, не более чем на 10 % [21].

Для проведения спектральных исследований была использована стандартная методика, описанная в литературе [18].

2.2.4 Методика осуществления математического моделирования

Методология синтеза и математическая модель основываются на том, что молекулярная масса и ширина ММР полимера влияют на его молекулярно-массовые, а вследствие, на их физико-механические характеристики.

Математическая модель базируется на кинетических уравнениях, составляющих процесс полимеризации ПММА. Академик В.В. Кафаров совместно с И.Н. Дороховым и Л.В. Дранишниковым в монографии «Системный анализ процессов химической технологии: процессы полимеризации» [22] сформулировали главные принципы модульного математического моделирования. Общий подход заключается в создании отдельных модулей комплексной математической модели, уточнении и отладке каждого в отдельности из этих модулей с последующим их объединением, уточнением и отладкой обобщенного математического формализма [23].

Затем на основе полученной модели проводятся численные эксперименты по оценке влияния на характеристики ПММА различных макрокинетических факторов (начальные концентрации реагентов, температура и т.п.). После чего осуществляется проведение реальных экспериментов с целью сопоставления теории и эмпирики. Если расхождение теоретически рассчитанных значений приведенной скорости полимеризации, средних молекулярно-массовых характеристик и функции ММР с этими же экспериментально исследованными характеристиками значительное, производится учет в формализме тех моментов, которыми изначально пренебрегали. Предлагаемый подход основан на совместном использовании численных и реальных экспериментов [24,25].

В комплексную математическую модель на данном этапе включены следующие модули: кинетический и прогностический. Кинетический модуль в виде системы дифференциальных уравнений, записанных на основе кинетической схемы и закона действующих масс, описывает изменение концентраций всех компонентов реакционной системы во времени. Молекулярно-массовые характеристики ПММА были получены на основе метода производящих функций из составленной системы дифференциальных уравнений [26]. Моделирование гель-эффекта было осуществлено в рамках моделей, предложенных Г.В. Жижином [27] и A.W. Hui [28]. Именно это позволило описать специфический вид и временную динамику функции ММР, наблюдаемых при получении ПММА полимеризацией в присутствии ферроцена (для данного материала это будет осуществлено впервые). Кинетические константы и их зависимости от температуры, закладываемые в данный модуль, были определены на основе литературных данных и в рамках проводимых исследований (решение обратной кинетической задачи). Фактически кинетический модуль - это математическая модель, учитывающая гель-эффект и формализующая кинетику получения ПММА полимеризацией в присутствии ферроцена. Кинетический модуль неоднократно уточнялся и отлаживался за счет сопоставления теоретически рассчитанной и экспериментально определенной приведенной скорости полимеризации, характеристик образцов ПММА с его прогнозируемыми величинами.

.3 Метрологическая обработка экспериментальных данных

Проведена оценка допущенной при измерении погрешности.

За наиболее вероятное значение величины а принимается среднее арифметическое значение результатов измерений

Расчет случайных ошибок базируются на математическом аппарате теории вероятностей с распределением Гаусса для случайных величин. Следует учитывать, что в условиях эксперимента при небольшом (n = 3¸10) числе измерений эти расчеты всегда носят оценочный характер.

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

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

Относительной погрешность i-го измерения величина

.

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

Средняя относительная погрешность рассчитывается по формуле

.

Доверительный интервал. Этот интервал определяется относительно ее среднего арифметического значения, принимаемого за наилучшую оценку истинного.

Приборная погрешность является паспортной характеристикой прибора и определяется классом точности прибора.

Обозначим класс точности

аmax.,

где xiприб. - максимально возможная абсолютная приборная погрешность i-го измерения, max - величина верхнего предела шкалы измерительного прибора.

Отсюда следует, что

,

а максимальная относительная приборная погрешность i-го измерения вычисляется по формуле

(%).

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

; .

Полная погрешность. Погрешность косвенных измерений.

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

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

3. Результаты и их обсуждение

В экспериментальной части работы был проведен синтез ПММА в лабораторных условиях. Для выявления влияния ферроцена на образующийся полимер, были проведены эксперименты по полимеризации ММА при 50 °С, 60 °С и 70 °С, концентрации инициатора 0,001 и 0,003 моль/л в присутствии и отсутствии ферроцена в полимеризующейся системе. В случае наличия в системе ферроцена, инициатор и ферроцен были взяты в соотношении 1:1.

На основании полученных результатов были построены следующие графики зависимости конверсии от времени (см. рисунок 3.1-3.5).

Рисунок 3.1 - Дифференциальные кривые полимеризации ММА в присутствии [АИБН]=1·10-3 моль/л и [МеЦ]= 1·10-3моль/л при 500С.

Рисунок 3.2 - Дифференциальные кривые полимеризации ММА в присутствии [АИБН]=1·10-3 моль/л и [МеЦ]= 1·10-3моль/л при 600С.

Рисунок 3.3 - Дифференциальные кривые полимеризации ММА в присутствии [АИБН]=1·10-3 моль/л и [МеЦ]= 1·10-3моль/л при 700С.

Рисунок 3.4 - Дифференциальные кривые полимеризации ММА в присутствии [АИБН]=3·10-3 моль/л и [МеЦ]= 3·10-3моль/л при 500С.

Рисунок 3.5 - Дифференциальные кривые полимеризации ММА в присутствии [АИБН]=3·10-3 моль/л и [МеЦ]= 3·10-3моль/л при 600С.

По представленным на рис. 3.1 - 3.5 данным мы можем сделать вывод о том, что при использовании ферроцена наблюдается вырождение ГЭ увеличение степени конверсии полимера.

Перейдем к численному эксперименту курсовой работы.

Цель составления и решения прямой кинетической задачи состоит в переходе от кинетической схемы полимеризационного процесса к молекулярно-массовым характеристикам получаемого полимера. Таким образом, для осуществления решения необходимо было выполнить ряд действий.

Первым этапом к теоретическому определению молекулярно-массовых характеристик ПММА, полученном с использованием инициатора АИБН и ферроцена, было написание кинетической модели реакции. Для этого необходимо было представить сложную реакцию в качестве совокупности стадий процесса, которые содержат только исходные вещества и продукты реакции (см кинетическую схему 3.1) [29,30]. Обозначения реагирующих компонентов представлено в таблице 3.1.

Кинетическая схема 3.1 - Реакции радикальной полимеризации

1.;

.;

.;

.;

.;

Таблица 3.1 - Обозначения компонентов кинетической схемы 3.1

ОбозначениеКомпонентОбозначениеКомпонент1212I

M

R(i)

kd

kp1

ktr1

kt1

kt2Константа скорости реакции распада инициатора

Константа скорости реакции роста цепи

Константа скорости реакции передачи цепи на мономер

Константа скорости реакции рекомбинации

Константа скорости реакции диспропор-ционирования

Следующим шагом стало составление системы балансных дифференциальных уравнений (система уравнений 3.1) по каждому компоненту в камках кинетического подхода. Результатом данного действия стало математическое описание расхода каждого компонента реакции во времени.

Система уравнений 3.1 - Дифференциальные уравнения для радикальной полимеризации

1) ;

) ;

)

;

)

;

)

;

)

;

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

Введение производящей функции основывается на уравнении

,

где Ф ? производящая функция,

? концентрация полимерных цепей со случайной степенью полимеризации,

? степень полимеризации,

? фиктивная переменная.

Для понимания сути метода введем понятие «момент». Статистический момент порядка ? это искусственно введенная математическая величина, которая является параметром распределения и выражается формулой

.

Следующим этапом стало выведение уравнения для статистических моментов нулевого, первого и второго порядков:

, при =0;

, при =1;

, при =2.

При подстановке к формуле 3.1 метода производящих функций, при =1 формула приобретет следующий вид:

.

При условии =1 как обязательное, возьмем первую производную от функции 3.1:

.

При взятии второй производной от функции 3.1 и подстановки в нее значение=1, уравнение приобретет следующий вид:

Таким образом, проводя параллели между уравнениями 3.3 и 3.6, 3.4 и 3.7, 3.5 и 3.8, мы можем сделать следующий вывод о связи производящей функции и ее производных со статистическими моментами:

;

;

.

В полученной кинетической системе было введено восемь производящих функций, уравнения и моменты которых представлены в табл. 3.2.

Таблица 3.2 - Производящие функции кинетической модели 1.

Производящая функцияСтатистические моменты ММР

Следующим шагом к написанию математической модели для получения молекулярно-массовых характеристик синтезируемого полимера стало приведение полученной системы кинетических уравнений к системе уравнений относительно статистических моментов ММР (см. систему уравнений 2).

Система уравнений 3.2 - Введение производящих функций

1)

)

)

)

)

Далее был проведен анализ научной литературы и были найдены значения неизвестных кинетических констант скоростей отдельных стадий реакции:

kd=1.58·1015exp(-15501/T), с-1; [31]

kp1=2.5·106exp(-23000/(R·T)), л/(моль·с); [32]

ktr1=9.958·105exp(-6197/T), л/(моль·с); [33]

kt1=(2.57·109·exp(-11700/(R·T)))/(6178·еxp(-2332/T)+1), л/(моль·с); [34] (3.24)

kt2=6178·exp(-2332/T) kt1, л/(моль·с); [34]

Вышеприведенный порядок действий описывает решение прямой задачи - нахождение молекулярно- массовых характеристик ПММА по кинетическим уравнениям. Однако, для того, чтобы модель была адекватна, необходимо учесть гель-эффект. Для описания ГЭ в случае полимеризации ММА используем зависимость, предложенную Хьюи [28], как функцию от СМ и температуры Т (К):

kt1ge = kt1 , л·моль-1·с-1

где ? степень конверсии, kt1ge ? константа скорости реакции рекомбинации c учетом ГЭ; kt1 - константа активации реакции рекомбинации в отсутствии ГЭ;

В данной модели константы А1, А2 и А3 нам не известны. В связи с этим, следующим шагом стало решение обратной задачи, то есть нахождение численного значения и температурной зависимости кинетических констант рекомбинации и диспропорционирования с учетом ГЭ. Для решения данной задачи была использована программа Wolfram Matematika, в которой был написан алгоритм, сравнивающий экспериментальные графики зависимости процента конверсии от времени (рис. 3.1 - 3.5) и графики той же зависимости, полученные теоретически, путем решения прямой задачи. Числовые значения констант, при которых происходило максимальное наложение были приняты за искомые (см. табл. 3.3).

Таблица 3.3 - Значения констант ГЭ при полимеризации ММА.

T, 0СА1А2А3505.083.33-4.85604.352.68-2.30

Следует отметить, что существует несколько наборов констант, удовлетворяющих заданным условиям. Для отбора единственно верных численных значений и проверки адекватности полученной модели, была выявлена зависимость констант от температуры:

A1 = 24.1 - 5.9010-2T, A2 = 21.6 - 5.6510-2T,

A3 = -91.3 - 0.268T.

Поучившееся при 700С наложение экспериментальной и теоретической кривой (с рассчитанными из зависимости константами) стало проверкой адекватности полученной модели (рис. 3.6).

Рисунок 6. Проверка адекватности получившейся модели.

(-- - теоретический график зависимости конверсии от времени;

? - экспериментальный график зависимости конверсии от времени)

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

Рисунок 7. Зависимость конверсии от времени ММА:

,2 - 500С; 3,4 - 600С, 5 - 700С; 1,3,5 - [АИБН]0=1 ммоль/л;

,4 - [АИБН]0=3 ммоль/л;

(-- - теоретический график зависимости конверсии от времени;

? - экспериментальный график зависимости конверсии от времени)

Для моделирования процесса полимеризации ММА в присутствии ферроцена необходимо дополнить существующую схему радикальной полимеризации ММА дополнительной кинетической схемой процесса (см. кинетическая схема 2). За основу была взята схема полимеризации тетрафункциональных мономеров [31].

Кинетическая схема 3.2 - Полимеризация в присутствии ферроцена

1.;

2.;

3.;

4.;

5.;

6.;

7..

Строение компонентов системы, участвующих в процессе полимеризации ПММА приведены в таблице 3.4.

Таблица 3.4 - Структурные формулы компонентов

ОбозначениеКомпонентОбозначениеКомпонентМс Расшифровка обозначений констант приводится в таблице 3.5.

Таблица 3.5 - Расшифровка обозначений констант

ОбозначениеРасшифровкаkHКонстанта скорости реакции замещения водородаkt3Константа скорости рекомбинации координации в координационной сфере металлаkc1, kc2, kc3Константы скорости реакции координации в координационной сфере металлаkp2Константа скорости реакции роста цепи в координационной сфере металлаkMetКонстанта скорости реакции замещения в координационной сфере металлаktr2Константа скорости реакции передачи на мономер в координационной сфере металла

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

По составленной кинетической схеме 2 была составлена система дифференциальных уравнений и методом производящих функций (см. таблицу 3.6) была решена сначала прямая (см. системы уравнений 3.3 - 3.4), затем, методом наложения теоретических и экспериментальных графиков, обратная задача. Следует отметить тот факт, что системами уравнений 3.3 - 3.4 модель радикальной полимеризации ПММА была дополнена и на основе объединения двух частей (радикальной и части с применением ферроцена) была составлена общая модель синтеза полиметилметакрилатного оргстекла в присутствии ферроцена.

Система уравнений 3.3 - Дифференциальные уравнения для процесса полимеризации ММА в присутствии ферроцена

1) ;

) ;

) ;

)

;

) ;

) ;

) .

Таблица 3.6 - Производящие функции кинетической схемы 2

Производящая функцияСтатистические моменты ММР12

12

Система уравнений 3.4 - Введение производящей функции в систему уравнений 3.3

1)

)

)

)

)

)

)

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

Таблица 3.7 - Кинетические константы полимеризации ММА в присутствии ферроцена

Обозначение константЗначение константkH1.34·107exp(-1453/RT)kt31.59·1010exp(-734/T)kc1, kc2, kc37.958·105exp(-2361/RT)kp22.57·1012·exp(6433/T)ktr21.58·1015exp(-2163/RT)kMet6178·exp(-6590/T)

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

Рисунок 3.8 - Проверка адекватности модели полимеризации ММА в присутствии ферроцена.

(-- - теоретический график зависимости конверсии от времени;

? - экспериментальный график зависимости конверсии от времени)

Из получившихся данных можно сделать вывод о том, что получившаяся модель адекватна и полностью описывает процесс полимеризации ММА в присутствии ферроцена.

Рисунок 3.9 - Математическая модель полимеризации ММА в программной среде Wolfram Mathematica.

В математической модели присутствует блок задания макрокинетических параметров реакции. (рисунок 3.9). Блок «вывод данных» показывает молекулярно-массовые характеристики получившегося полимера. (рисунок 3.10). Таким образом, разработанная в рамках курсового проекта модель представляет собой основу для осуществления синтеза полимерного материала с необходимыми показателями.

Рисунок 3.10 - Вывод молекулярно-массовых характеристик синтезированного ПММА

Заключение

В ходе курсовой работы была разработана математическая модель кинетики полимеризации ММА в присутствии ферроцена.

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

Адекватность созданной модели полимеризации доказывается совпадением теоретически выведенной зависимости конверсии от времени на экспериментально полученную кривую. Коэффициент корреляции при этом равен 0.98.

Анализируя графики экспериментальных кривых можно заметить, что изменение степени конверсии зависит от условия проведения эксперимента. Так, с увеличением температуры с 50 °С до 70 °С наблюдается увеличение степени конверсии. Однако, при увеличении концентрации инициатора, какого либо существенного изменения в конверсии не наблюдается. Также по экспериментальным графикам было установлено, что внесение в полимеризующуюся систему ферроценов значительно улучшают характеристики синтезируемого полимера и заметно сглаживают гель-эффект. Это говорит об эффективности использования металлоценовых соединений.

Применение ферроцена не противоречит основной идее курсового проекта, заключающейся в использовании в качестве упаковки материала без посторонних примесей. Количество вносимого в систему ферроцена имеет порядок 10-3 моль/л. Данное количество незначительно по сравнению с концентрациями вводимых стабилизаторов, пластификаторов и прочих добавок (количество порядка 10-1 моль/л каждого).

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

Выводы

1.Впервые получены константы для описания радикальной полимеризации ММА по модели Хьюи;

.Впервые была разработанная кинетическая схема для описания процесса полимеризации ММА в присутствии ферроцена;

.Впервые были получены кинетические константы для описания процесса полимеризации ММА в присутствии ферроцена;

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

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

1.Светикова, В.С. Дизайн упаковки, как инструмент воздействия на объемы продаж и управлении себестоимостью продукции / В.С. Светикова // Тара и упаковка. -2009.- №2-3. - С. 58-59

2.Розенберг, Н.А. Упаковка как фактор повышения конкурентоспособности товара / Н.А. Розенберг // Тара и упаковка. -2012.- №3. - С. 30-31

.Материалы для упаковки : учебное пособие / сост. А.А. Шагивалеев, А.А. Овчинников, А.Н. Николаев; Казан. госуд. гехн. ун-т. - Казань, 2010. - 82 с.

.Борисов, А.Ю. Уникальный проект о безопасности упаковки / А.Ю. Борисов // Тара и упаковка. -2009.- №2-3. - С. 54-56

.Тренина, В.Е. Биоразлагаемая упаковка. Мифы и реалити / В.Е. Тренина // Тара и упаковка. -2013.- №3. - С. 28-29

.Кольм, Д.А. «За» и «против» внедрения биоупаковки на современном этапе / Д.А. Кольм // Тара и упаковка. -2009.- №5. - С. 43-44

.Локс Ф. Упаковка и экология / Ф. Локс. - Москва.: Издательство МГУП, 1999. - 219 с.

.Немальцева, М.С. Полимеры в России / М. С. Немальцева // Тара и упаковка. -2010.- №6. - С. 5-6

.Потем, Т. А. Влияние полимерных упаковочных материалов на организм человека / Т.А. Потем, И.Э. Егоров // Тара и упаковка. -2011.- №2. - С. 35-37

.Жуков, А.А. Повлияет ли присоединение к ВТО на индустрию ПЭТ / А.А. Жуков // Тара и упаковка. -2012.- №2. - С. 53-54

.Рябичева, Е.Д. Америка России подарила ПВХ / Е.Д. Рябичева // Тара и упаковка. -2011.- №1. - С. 28-29

.Справочник. Тара из полимерных материалов / М.Г. Соломенко, В.Л. Шредер, В.Н. Кривошей. - Москва.: Химия, 1990. - 400 с.

.Литвининко, Г.И. О природе гель-эффекта при радикальной полимеризации стирола / Г.И. Литвиненко, М.Б. Лачинов, Е.В. Саркисова, В.А. Каминсткий // ВМС. - 1994. - том 36. - №2. - С. 327-333.

.Будтов, В.П. Основные экспериментальные закономерности радикальной полимеризации в области гель-эффекта и их теоретическая интерпретация / В.П. Будтов, Н.Г. Подосепова // ВМС. - 1988. - Том А(ХХХ). - №7. - С. 156-161.

.Головочесова, О.И. Радикальная полимеризация метилметакрилата и стирола в присутствии клатрохелатных и пивалатных комплексов железа и кобальта и гетероциклического производного ферроцена : дис. канд. хим. наук /О.И. Головочесова. - Уфа,2011. - 215 с.

.Дьяконов, А.Д. Mathematica 5.1/5.2/6 Программирование и математические вычисления / А. Д. Дьяконов. - Москва.: ДМК, 2008. - 572 с.

.Насибуллин, И.И. Формирование активных центров при полимеризации метилметакрилата и стирола в присутствии металлоценовых инициирующих систем /И.И. Насибуллин. - Уфа, 2013. - 208 с.

.Гладышев, Г.П. Радикальная полимеризация при глубоких степенях превращения / Г.П. Гладышев, В.А. Попов. - Москва.: Наука, 1974. - 244 с.

19.Denisov, E.T. Oxidation and Antioxidants in Organic Chemistry and Biology / E.T. Denisov, I.B. Afanasev. - FL: CRC Press, Taylor and Francis Group, 2005. - 982 p.

20.Стыскин, Е.Л. Ициксон Л.Б., Брауде Е.В. Практическая высокоэффективная хроматография / Е.Л. Стыскин, Л.Б. Ициксон, Е.В. Брауде. - М.: Химия,1986. - 290 с.

21.Kwak, Y. А systematic Study on Activation processes in organotellurin-mediatel living radical polymerization of styrene, methyl methacrylate, methyl acrylate, and vinyl acetate / Y. Kwak, A. Goto, T. Fukuda // Macromolecules. - 2004. - V. 37. - No. 4. - Р. 121.

22.Берлин, Ал.Ал. Кинетика полимеризационных процессов / Ал.Ал. Берлин, С.А. Вольфсон, Н.С. Ениколопян.- М.: Химия, 1978. - 320 с.

.Берлин, Ал.Ал. Кинетический метод в синтезе полимеров / Ал. Ал. Берлин, С.А. Вольфсон. - М.: Химия, 1973. - 256 с.

24.Кучанов, С.И. Методы кинетических расчетов в химии полимеров / С.И. Кучанов. - М.: Химия, 1978. - 368 с.

25.Улитин Н.В. Использование кинетического подхода для моделирования радикальной полимеризации стирола, осуществляемой в присутствии дибензилтритоикарбоната / Н.В. Улитин, И.И. Насыров, Т.Р. Дебердеев, А.А. Берлин // Химическая физика, 2012. -Том 31. - № 11, С. 66-74.

26.Жижин, Г.В. Гель-эффект при радикальной полимеризации / Г. В. Жижин, И.А. Обухова // Математическое моделирование. - 1997. - Том 9. - № 11. - 3-13.

27.Hui A.W., Hamielec A.E. // J. Appl. Polym. Sci. - 1972. V. 16. P. 749.

28.Controlled/Living Radical Polymerization: Part 1 / ed. by K. Matyjaszewski. - Washington, D.C.: American Chemical Society, 2009. - 423 p.

.Kuzub L.I.,Peregudov N.I., Irzhak V.I. // Polymer Scienct. 2005. V. 47(A). № 10. Р. 1063.

30.Гайсин, Ф.Р. Моделирование методом Монте-Карло трехмерной свободно-радикальной полимеризации тетрафункциональных мономеров (на решетках разных размеров и геометрии в рамках формирования единичного трехмерного структурного элемента) / Ф.Р. Гайсин, Ю.М. Сивергин, С.М. Усманов. - Уфа: Гилем, 2009. -180 с.

31.Brandrup J., E.H. Immergut, Grulke E.A. // Polymer handbook, fourth edition.

.Nising P., Meyer T. // Ind. Eng. Chem. Res. 2004. № 43. Р. 7220-7226.

.Bamford C.H. The kinetics of vinil polymerization by radical mechanisms / C.H. Bamford, W.G. Barb, A.D. Jenkins, P.F. Onyon. - London: Butterworths scientific publication, 1958.

34.Улитин, Н.В. Формализация гель-эффекта при радикальной полимеризации метилметакрилата / Н.В. Улитин, В.А. Сидельникова, В.М. Янборисов, С.В. Колесов // Реф. в: Математическое моделирование процессов и систем.: науч.-реф. Сб. - 2013. - Вып. 5. - С.13-14.

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

Больше работ по теме: