Моделирование термомеханического поведения полимерных материалов в условиях фазовых переходов

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

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

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

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

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

Научная новизна работы состоит в следующем:

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

— Экспериментально исследованы поля остаточных напряжений в круглой пластине из полиэтилена низкого давления (ПЭНД) для двух различных режимов охлаждения;

— Предложен и реализован эффективный численный алгоритм решения краевой задачи термомеханики кристаллизующейся вязкоупругой среды с учетом зависимости теплофизических и реологических свойств материала от температуры;

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

— На основе численного анализа напряженно-деформированного состояния (НДС) изделий из полимерных материалов для конкретных технологических процессов впервые установлены существенные качественные и количественные отличия в эволюции напряженного состояния в упругой и вязкоупругой кристаллизующейся механической системе.

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

Основные положения работы обсуждались на научно-технической конференции ГИТУ «Проблемы прикладной математики и механики» (Пермь, 1998 г.), на Всероссийской конференции молодых ученых «Математическое моделирование в естественных науках «(Пермь, 1998 г.), на 12-ой зимней школе по механике сплошных сред (Пермь, 1999 г.), на Всероссийской конференции молодых ученых «Математическое моделирование физико-механических процессов» (Пермь, 1999 г.), на X-ой юбилейной международной конференции «Вычислительная механика и современные прикладные программные системы» (Переяславль-Залесский, 1999 г.), на УШ-ом Всероссийском съезде по теоретической и прикладной механике (Пермь, 2001) и отражены в публикациях /66−75/.

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

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

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

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

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

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

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

— сгущения конечно-элементной сетки в области больших градиентов температуры;

— измельчение шага по времени в ходе численного решения задачи.

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

В заключении отражены основные результаты диссертационной работы и выводы.

4.5. Выводы по главе.

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

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

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

Установлено, что учет вязкоупругих свойств кристаллической фазы при определении НДС затвердевающей полимерной системы приводит к существенному снижению расчетных уровней технологических напряжений в изделиях (для рассмотренных технологических процессовдо 9,5 раз).

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

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

ЗАКЛЮЧЕНИЕ

.

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

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

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

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

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

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

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