Математическое моделирование деформирования вязкоупругих структурированных полимерных (композитных) систем тема автореферата и диссертации по механике, 01.02.04 ВАК РФ

На правах рукописи

БЕЛЯЕВА Надежда Александровна

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ДЕФОРМИРОВАНИЯ ВЯЗКОУПРУГИХ СТРУКТУРИРОВАННЫХ ПОЛИМЕРНЫХ (КОМПОЗИТНЫХ) СИСТЕМ

01.02.04 - Механика деформируемого твердого тела

Автореферат диссертации на соискание ученой степени доктора физико-математических наук

00345817Б

Черноголовка — 2008

003458176

Работа выполнена в Институте структурной макрокинетики и проблем материаловедения Российской Академии наук (ИСМАН)

Научный консультант: доктор физико-математических наук, профессор

Столин Александр Моисеевич

Официальные оппоненты: доктор физико-математических наук, старший научный сотрудник Патлажан Станислав Абрамович

доктор физико-математических наук, профессор Волынцев Анатолий Борисович

доктор технических наук, старший научный сотрудник Скульский Олег Иванович

Ведущая организация Институт прикладной механики РАН

(г. Москва)

Защита состоится 22 января 2009 г. в 10 часов

на заседании диссертационного совета Д 004.012.01 Института механики сплошных сред УрО РАН по адресу: г. Пермь, ул. Акад. Королева, 1, тел. (342)2378388, адрес ¡ШегпеЬсайта' www.icmm.ru

С диссертацией можно ознакомиться в библиотеке Института механики сплошных сред УрО РАН.

Автореферат разослан « & {» декабря 2008 г.

Ученый секретарь

диссертационного совета

Березин И.К.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

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

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

Работа выполнялась в соответствии с федеральной целевой программой «Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России» на 2007 - 2012 гг. при выполнении госконтракта 02.513.11.3377 от 26 ноября 2007 г.

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

1. Исследование процессов самоорганизации течения вязкой структурированной жидкости в области "сверханомалии вязкости".

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

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

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

5. Численный анализ качественно различных режимов твердофазной экструзии полимерных материалов. Нахождение областей реализации этих режимов в зависимости от характерных времен процесса.

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

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

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

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

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

- режим объемного отверждения цилиндрического изделия;

- режим фронтального отверждения цилиндрического изделия;

- режим объемного отверждения сферического изделия;

- режим фронтального отверждения сферического изделия.

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

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

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

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

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

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

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

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

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

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

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

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

Результаты теоретических исследований использовались в Научно-образовательном центре ТамбГТУ-ИСМАН «Твердофазные технологии» при получении изделий методами твердофазной экструзии из различных полимерных композитов. Получено заключение о практическом использовании основных положений диссертации в Институте структурной макрокинетики и проблем материаловедения РАН (г. Черноголовка) и Тамбовском государственном техническом университете (г. Тамбов).

Личный вклад соискателя. Все приведенные в диссертации результаты получены лично автором или при его непосредственном участии. Монография «Математические модели деформируемых структурированных материалов», отражающая главные направления диссертации, опубликована без соавторов. В постановке задач по отверждению цилиндрического изделия, обсуждении результатов принимали участие к.ф.-м.н. Клычников Л.В., д ф.-м н. Худяев С.И., д ф.-м н. Давтян С.П.. В постановке задачи куэттовского течения структурированной жидкости принимал участие д.ф.-м.н. Худяев С.И.; в постановке и обсуждении результатов исследования процесса твердофазной экструзии - д.ф.-м.н., профессор Сталин A.M., д.т.н. Стельмах Л.С. При подготовке статьи: Беляева Н.А, Сталин А.М., Пугачев Д.В., Стельмах Л.С. Неустойчивые режимы деформирования при твердофазной экструзии вязкоупругих структурированных систем // ДАН, 2008. Т. 420. № 6. С. 777-780. - соавтором Пугачевым Д.В. был поставлен эксперимент, результаты которого используются в названной публикации.

Апробация работы. Результаты диссертационной работы докладывались на: Всероссийской конференции «Переработка полимерных материалов в изделиях» (Ижевск, 16-19 ноября 1993 г.), 8-й Всероссийской конференции «Математические методы в химии» (Тула, 27-29 окт. 1993 г.), 2-й Международной научно-технической конференции "Актуальные проблемы фундаментальных наук" (Москва, январь 1994 г.), Симпозиуме «Математические модели деформируемого тела» (СПб, июнь 1994 г.), XXXVI Уральском семинаре «Механика и процессы управления» (Екатеринбург, 2006 г.), III Международной летней научной школе «Гидродинамика больших скоростей и численное моделирование» (Кемерово, 22-28 июня 2006 г.), Четвертой Всероссийской школе-семинаре по структурной макрокинетике для молодых ученых (Черноголовка ИСМАН , 22-25 ноября 2006 г.), конференции «Реология и физико-химическая механика гетерофазных систем» (Карачарово, 23 - 28 апреля 2007 г.), Международной конференции «XVIII сессия Международной Школы по моделям механики сплошной среды» (Саратов, 2007 г.), научной конференции-семинаре «Теория управления и математическое моделирование» (Ижевск, 2008 г.), пятой Всероссийской конференции с международным участием «Математическое моделирование и краевые задачи» (Самара, 2008 г.); научном семинаре кафедры теории упругости механико-математического факультета МГУ (Москва, 1999 г.), научном семинаре в Отделе математики Коми научного центра УрО РАН (Сыктывкар, 2004 г.), научном семинаре лаборатории гидродинамики математико-механического ф-та СПб технического университета (СПб, 2004) , научном семинаре ИСМАН РАН (г. Черноголовка, 2008 г.), научном семинаре Института механики сплошной среды (ИМСС) РАН (Пермь, 2008 г.).

Публикации. Всего по теме диссертации опубликовано 45 научных работ, из них 9 научных работ- в ведущих журналах, рекомендованных ВАК.

Структура и объем диссертации. Диссертация состоит из введения, пяти глав, заключения, библиографического списка из 191 наименования, приложения; содержит 193 страницы текста, 102 рисунка, 6 таблиц.

Основные положения, выносимые на защиту:

1. Теоретическое исследование различных типов течений структурированных текучих систем на основе взаимосвязанного рассмотрения реодинамических процессов и кинетики структурных превращений. Численная реализация структурно-однородных и пространственно структурно-неоднородных режимов деформирования. Анализ особенностей формирования диссипативных структур в условиях явления сверханомалии вязкости.

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

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

4. Анализ динамики процессов деформирования структурированных систем и объяснение характера различных критических (пороговых) явлений и колебательных режимов. Механизмы «обратной» связи неустойчивого режима деформирования соответственно в области малых и больших скоростей деформации.

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

6. Численный анализ особенностей качественно различных режимов экструзии в зависимости от характерных времен процессов уплотнения, выдавливания и структурных превращений.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

В процессе самоорганизации структурированных систем наблюдаются нелинейные эффекты динамики систем. Впервые об этом явлении было заявлено еще в 1954 году на втором Международном Конгрессе по реологии в Лондоне [1]. Здесь исследовалась неустойчивость течений Куэтга и Пуазейля, приводящая к релаксационным автоколебаниям, при этом неустойчивый режим течения возникал на падающем участке скоростной характеристики трения скольжения.

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

Взаимосвязанное рассмотрение течения и кинетики взаимного превращения структурных единиц [4-6] дает возможность выяснить интересные особенности структурированных систем. В рамках представленной в этих работах модели обосновывалась возможность существования области сверханомалии вязкости в сложных системах химической кинетики.

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

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

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

¿¡1 и цг —вязкости компонент А иАг, соответственно, а - доля А в смеси А, и А . Если = /*2> 70 вязкость ц не зависит от а - жидкость является ньютоновской; если //, > цг - псевдопластической, //, < /л2 - дилатантной.

Рассмотрим течение Куэтта структурированной жидкости в плоском зазоре. Пусть структурированная вязкая несжимаемая жидкость заполняет полосу между плоскостями ? = О и £ = Л, а ее течение вызывается движением плоскости £ = к в направлении оси т] со скоростью, не зависящей от т/ и £. Система уравнений, описывающая поведение рассматриваемой жидкости:

Здесь /I— вязкость, зависящая от степени структурных превращений а = а(^) ставляющих компонент жидкости, причем

со-

(1)

(2)

да да

(3)

(4)

На подвижной стенке рассмотрим граничные условия двух типов:

(5)

или

Здесь

Эи

К

скорость

деформации,

р —плотность,

Ф{а,у) = кг(\-а)-к1а = к1

l~a-a-±<ixp(pv(a)r+qr2)

-суммарная скорость превра-

щения структуры.

Рассмотрим однородные стационарные состояния системы (1)-(б), то есть решения вида:

Э и

а = const, — = у = const.

9S

При граничных условиях (4), (5), однородная скорость деформации имеет вид: у = и0/ h. Однородные стационарные значения степени структурных превращений а находятся из уравнения: ф(а,у) = 0 или

1-а-а^-ехр(р/1(а)у + дуг)

= 0.

Сделав в последнем уравнении замену

К 1 /А . Ч

кг А (рМг) получим реологическое уравнение кривой течения

_*(1+ЛГехр {у + 8хг)) х{1 + к)

= 5,

А + 1 + % ехр^у + S х1} Л + 1 + *'

(7)

10 х.х0х, X

Рис. 1а Реологическая кривая псевдопластической жидкости, Л = -0.7

Рис. lb. Реологическая кривая дилатантной жидкости, Л = 5

На рис. 1 а,Ь изображены реологические кривые с областями "сверханомалии вязкости" - явление падения напряжения ст от ст1(уа) до <У„(у„) при увеличении скорости деформирования у для псевдопластической жидкости {рис. 1 а), и падение скорости деформации у от уг (х„) до уи (*„) при увеличении напряжения с для дилатантной жидкости (рис. 2Ь). В областях сверханомалии происходит потеря устойчивости однородного стационарного решения, возникает пространственная неоднородность при течении жидкости. При выборе граничных условий (6), где <т0 удовлетворяет <7„ < сг0 < <Ув, мы полу-

чим три однородных стационарных решения а,, а2, а}: а, <а2 < аг и а, — а, (<т0); причем о,, а3 - устойчивые однородные состояния, а аг - неустойчивое.

Численное решение системы (1)-(б), проведенный анализ результатов, позволили построить пространственно-неоднородные решения, которые возникают в результате потери устойчивости - это монотонно-убывающая или монотонно возрастающая диссипатив-ная структура (рис. 2).

4 " Рис. 2. Пространственно-

временное распределение степени структурных превращений а(г,т): 1(т = 0)Д1)Д15)Ж37); Л = -0.7,аг =0.5125 /? = 0.001 А"=0.1,^=10"5, р, =3,9, =6Д2 =0.5;

<^5 =0-712

ТгГ

Установившееся пространственно-неоднородное решение а = ) удовлетворяет необходимому стабилизационному соотношению

1

|(а(Л)-а2)Л = 0, (8)

о

где а = аг- соответствующее неустойчивое однородное стационарное решение.

Анализ пространственно - временного распределения степени структурных превращений демонстрирует установление пространственно неоднородного течения жидкости -диссипативной структуры. График этой структуры представляет монотонно убывающую или монотонно возрастающую зависимость а = а(£). Для данной структуры выполняется стабилизирующее соотношение (8), что подтверждают результаты расчетов.

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

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

Внутренние технологические напряжения, возникающие в изделиях в ходе их изготовления, тот или иной уровень остаточных напряжений оказывают существенное влияние на прочность изделий в зависимости от условий их эксплуатации. Имеющиеся теоретические и экспериментальные исследования по объёмному (послойному) отверждению подобных изделий направлены на изучение кинетики накопления внутренних напряжений на различных этапах формирования изделий [9, 10]. В работе [11] рассмотрено образование внутренних напряжений в условиях фронтального отверждения симметричных изделий (сферы, цилиндра) при заданном заранее законе распределения температуры отверждения по радиусу изделия и в предположении, что в основном превращение отвергаемого материала происходит в узком температурном интервале вблизи температуры адиабатического разогрева.

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

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

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

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

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

На возможность одновременного протекания ре акций полимеризации и кристаллизации впервые указано в [15]. Изменение структуры полимеров и полимерных композиций, обусловленное деформациями, является в настоящее время неоспоримым экспериментально установленным фактом [16]. Однако кристаллизацию в процессе полимеризации стали изучать на основе совместного математического, физического и химического подхода только недавно. Многие вопросы этого сложного явления остаются нерешенными. В современных исследованиях для учета влияния кристаллизации на механические свойства

полимеров предлагаются различные подходы [17, 18]. рис_ з. Модель формируе-Так, в работе [17] описание свойств полимерных ма- мого изделия

териалов предлагается с помощью параллельно со-

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

В § 2 рассматривается макрокинетическая модель совмещенного процесса полимеризации и кристаллизации формования полого осесимметричного изделия (7?,<г^Л). В начальный момент мономер имеет температуру окружающей среды Т0 (объемный режим реакции). На внутренней и внешней границе образца происходит теплообмен с окружающей средой (рис. 3). В результате реакции и теплообмена в образце устанавливается неоднородное температурное поле. Определение температурного Т = Т(гД) и конверсионных полей (глубины полимеризации а — С((г, I) и кристаллизации Г] = Т](г, I)) в процессе фазового перехода производится на основе решения системы уравнений:

дТ , (ЭгТ пдт) Л да Л дт]

Л (9)

(10) (11)

г) СУ

~^к0^[-и/ЯиТ]{1~а){£> + ау, ^ = к0г ехр[-£" / ЛиТ - уТ{ !Т(ТГ - Г)](гг + т}){ат)ь - г,),

(9) - уравнение теплопроводности с нелинейным источником тепла и два кинетических уравнения относительно степени структуризации: (10) - полимеризации, (11) - кристаллизации; п - коэффициент формы (для цилиндра л=1, для сферы /1=2). Начальные и граничные условия:

,=0: а = 0; 77 = 0; =Г0; (12)

д^Щт-т0)1=о. (13)

Т =Т(г,I), ОС = а(г,(), Т] = Т](г,1)-температура, глубина полимеризации и кристаллизации, соответственно, в точке г образца в момент времени I. В уравнениях (1.1)-(1.3) исполь-

зованы обозначения: с—удельная теплоемкость материала; р — плотность; Ау - коэффициент температуропроводности; 0,р — тепловой эффект полимеризации; ()сг —тепловой эффект кристаллизации; ^Д0 -константы скоростей полимеризации и кристаллизации; и — энергия активации процесса полимеризации; Е" — энергия активации процесса кристаллизации; £х,Е2— критерии автокаталитично-сти процессов полимеризации и кристаллизации; Яи —универсальная газовая постоянная; у/ — харак-

580 т

540 -

500 -

460 -

420

380

т,к

I

; 7 11

■ : кр

-^--

50

55

60

65

70 71.5 75

Рис. 4а. Низкотемпературная (I) и высокотемпературная (II) области реакций в отсутствие кристаллизации

терная температура полимера; Т0 -начальная температура мономера и окружающей среды; й0,/г- коэффициенты теплообмена с окружающей средой; Т{ -равновесная температура плавления; 7/4 - равновесная степень кристалличности, вычисляемая по формуле

Ль = 0.52^1-(Г/7^ )4.

Влияние условий полимеризации на свойства формируемого изделия наиболее отчетливо можно проследить путем определения высокотемпературного или низкотемпературного режима проведения реакции посредством критериев Семенова

К-р=идрИ2ЦГ,)/^иТ.2-, ксг =иГ2сгКгкг{Т.)1 ,Т. - характерная температура

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

АГ* =71.5: при <71.5- низкотемпературная, >71.5 —

высокотемпературная область.

Значительное влияние кристалли-*

зации на величину Кр заметно при кс, ~103...Ю4. Так, при Ксг = 80 критическое значение

«"^=71.5, а при

Кег = 104 — Кр = 56. При переходе из области I в область П изменяется температурный режим проведения процесса, следовательно, происходит существенное изменение напряженного состояния формируемого изделия. Изменение максимальной температуры в процессе отверждения высокотемпературного и низкотемпературного режимов отверждения представлено на рис. 5. Повышение температуры на начальной стадии процесса соответствует интенсивному протеканию реакции полимеризации, последующее падение обусловлено наличием теплообмена с окружающей средой и уменьшением температурного вклада от полимеризации. Сравнение поведения кривых 1, 3 подтверждает влияние кристаллизации на температурное распределение и в низкотемпературном режиме. На кривой 4 наблюдается еще один подъем температуры, обусловленный интенсивным протеканием кристаллизации. Таким образом, в высокотемпературном режиме тепловые эффекты от полимеризации и кристаллизации могут быть разделены во времени, тогда как в низкотемпературном режиме этого сделать нельзя.

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

В § З-б рассматриваются модели определения напряженного состояния формируемого изделия: объемный (§ 3) и фронтальный (§ 4) для цилиндрического изделия, объемный

Рис. 4Ь. Низкотемпературная (I) и высокотемпературная (II) области реакций при совместном протекании реакций полимеризации и кристаллизации

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

Г=1

<Z>(£) = <T(t)+ J K(t-r)ff(r)dr,

г=0

где <р(£) — некоторая функция деформации, K(t-г) —так называемое ядро последействия.

О 150 350 500 2500 4500

Рис. 5. Зависимость максимальной температуры процесса отверждения от времени: (1)- низкотемпературный (2)- высокотемпературный (К?=1\.5) режимы без учета кристаллизации; (3)- низкотемпературный (Ки =70, ка =80), (4)-высокотемпературвый (=56, ксг =104) режимы

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

(14)

Эг г Эг г

(14) - уравнение равновесия, (15) - условие совместности деформаций; <7П, —радиальная и окружная компоненты напряжений; Е„, Е^ — радиальная и окружная компоненты деформаций. Полные компоненты деформации в (15) являются суммой вязкоупру-гой Е , температурной ЕТ и химической составляющих:

£ = £ +£' + £*, (16) где =С0(Т-Т0), а0— аналог коэффициента линейного температурного расширения материала; для простоты будем считать его постоянным; Т = Т(г, I) — текущая температура материала; £ * = £р + £а, £"р = к,ОГ -усадка материала вследствие полимеризации, £а — к2Т] -усадка вследствие кристаллизации, к,, к2 —константы.

В двумерном случае вязкоупругую компоненту деформации будем искать в тензорном виде

'1/Е -WE'

^-WE 1/Еу

f„(t-r) f„(t-г) vfpr(t-r) f„(t-r)

ojj ,т)

(17)

Т,

V — коэффициент Пуассона, Е — модуль упругости.

Применяя к соотношению (17) преобразование Лапласа, используя условие несжимаемости материала, найдем зависимости между изображениями Г^(р) и, следовательно,

соотношения между их оригиналами fjj:

ргДР) = Fyr(p) = = -F„(P).

Е Е(р)

Для стандартной модели (рис. 6) вязко-упругого тела

Л = (Е, +E1)l/3,El =Е = const; Е2= 3■ Ю10(Г„/Т-1); ц = Е2!Р,Р = Р{Пр КЛр ~ri), Д = Аехр(-£*"/RUT + Дог).

Таким образом, определены неизвестные функции ^ тензорного соотношения (17). Для определения компонент напряжений введем в рассмотрение функцию напряжений Ф = Ф(г,1):

Ф ЭФ

сг„(г,0=—. <тщ(г.0=—. (18)

Подставим (16) с учетом (17), (18) в условие (15). Принимая во внимание (14), относительно введенной функции Ф = Ф(г, I) получим уравнение

д_ Эг

r t

Эг

Эг г

Эг

Е

L V

V J r

= 0,

где в(х, t) = £•*"" + el. Интегрирование последнего уравнения приводит к уравнению

ЭФ(Г,1) Ф(г, t) Ь

Эг

= -Е

Jf(t-r)

ЭФ(г,р Ф(г,р Эг + г

dr =

'г в - в

еМ—^dr

I Г

(19)

представимому в виде системы интегро-дифференциальных уравнений 9J(r, t) J(r,t)

Эг

- = F(r,t) + ^(t),

Здесь f=f[r(t)'E, F(r,t) = -E

Ф(г, t) + jf (t - г)Ф(г,т)с!г = J(r, t).

0 ^

rf6>(r,t)-0(r,t)

(20)

dr

— функция, значения

которой определяются на основе решения макрокинетической модели (9)-(13); у/{t) -произвольная функция времени. Предположим, что на внутренней поверхности сумма нормальных напряжений равна нулю, т.е.

+ (2D

тогда, следуя (19), получим y/{t) = —F(Rt, t).

Внешнюю поверхность будем считать свободной от напряжений crj^ = 0, значит

Ф(И, t) = 0 и из (20) J(R, t) = 0, тогда первое уравнение системы (20) разрешается

1 R 1 J(r,t) = — ¡F(r,t)rdr + — F(Rvt)(R2 - г2). (22)

r\ 2 г

А второе - можно представить, учитывая f (t) = (Л- JU) exp(-fit), как

Ф(г, t) + J(A-/i) exp [-//(t - г)] Ф(г, r)dr = J(r, t). (23)

о

Записывая уравнение (23) в изображениях и разрешая полученное алгебраическое уравнение относительно Ф(г, р) — изображения функции Ф(гД), применяя обратное преобразование, находим функцию напряжений

Ф(г, t) = J(r, t) - } (Л - JU) exp [-¿(t - т)\J(r, r)dr,

(24)

где Дг, I) определяется формулой (22). Из уравнения (24), с учетом (18), получим формулы для определения компонент напряжений

(25)

(26)

Расчеты проводились для низкотемпературного и высокотемпературного режимов реакции при следующих значениях параметров задачи:

с = 0.24 -10' кал/(кг ■ град); р = 1.1-10-3 кг/см1;

Л0 = 9 -Ю"4 кал!(см ■ с ■ град)-,

£>„ = 1&кал/см*; ()кр =35кал/см1;

к0х =5-10 5с"';*01 = 2 • 104 с"';

и = 13-Ю3 кал/.моль; = 2кал/(град ■моль);

£, = 0.18; = 0.05;

= 8.8-103 кал/моль; (С = 225 К; Г„, =493^;Г0 = 373АГ; Д = 10сл(; Д, = 1сл; А0 = Л = 0.5сл<-'; ^ = 0.33; к, = к 2= -0.05; ог0 = 0.25 • 10"2 град'1; А=\ОгсП; Е"' = 5 Ю3 кал! моль; Е, = Е = 10' кал / моль; Д = 5.

£г

а

Рис. 6. Стандартное вязко-упругое тело

378

377

На рис. 7-9 представлены результаты численного анализа низкотемпературного (кп =70, ккр =80) режима объемного отверждения цилиндрического изделия. Максимум температуры (рис. 7) в процессе от- 374 верждения наблюдается внутри изделия, рост температуры происходит вследствие интенсивного протекания полимеризации. Анализ конверсионных полей (рис. 8(а, Ь)) указывает на режим параллельного протекания полимеризации и кристаллизации.

376

375

373

Рис. 7. Пространственно-временное распределение температуры Т = T(r,t), t(c) : 0(1), 660(2), 1321(3), 1981(4), 2642(5), 3303(6), 3963(7), 4624(8), 5284(9), 5945(10)

Ь)

Рис. 8(а, Ь). Пространственно-временное распределение глубины полимеризации - а) и глубины кристаллизации - Ь) в низкотемпературном режиме реакции; условия на рис. 7

На рис 9(а, Ъ) представлена динамика радиальной и окружной компонент напряжений. Радиальные напряжения в ходе отверждения являются, в основном, растягивающими, окружные - сжимающими. Наибольшее значение напряжения достигают вблизи внутренней поверхности.

Рис. 9(а, Ь). Пространственно-временное распределение радиальной и окружной компонент напряжений; условия на рис. 7

После освобождения от закрепления внутренней поверхности изделия получим распределение результирующих напряжений, характеризующих напряженное состояние образца по окончании процесса отверждения (рис. 10).

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

Далее в работе рассматривается (§ 4) модель фронтального отверждения цилиндрического изделия. Задача о фронтальном формировании (отверждении) осесимметричного изделия решается с точки зрения непрерывно наращиваемого твердого тела. Отличие от постановки задачи для фиксированного тела вытекает из того факта, что элемент растущего тела в окрестности какой-либо точки (г) деформируется совместно с другими элементами, лишь начиная с момента его отверждения I (г), совпадающего с моментом присоединения этого элемента к образованной (затвердевшей) области. Чтобы соблюсти корректность постановки рассматриваемой задачи и учесть динамический характер протекающего процесса отверждения (непрерывного наращивания), основные уравнения и граничные условия записываются в продифференцированном по времени виде.

Система уравнений макрокинетической модели (9)-(П) дополняется в рассматриваемом случае начальными и граничными условиями вида

г = 0: а = 0; ч = 0; Г|^=Г0; ГЦ =Г;

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

При определении напряженного состояния формируемого изделия принята стандартная модель вязкоупругого тела (рис. 4). Вязкоупругие компоненты радиальной и тангенциальной компонент деформации подчиняются соотношениям

¿;(г, 0 = ^(¿„(М) - у&„(х, 0) + Шо-^г.О + ^(О)а^(г, I) +

I I

+ | ¡пЦ-т)а„(т,г)йт+ |

1-(г) |'(г>

и, кг/см3

? ... 4/ 5 . . е

(<т = ст^{г)

1

Рис. 10. Результирующие напряжения в низкотемпературном режиме

C^t) = - v&a{r, t)) + ^(0)СТ„(Г, t) + fff(0)aw(r,t) -

t t + J -г)а„(т,т)йт+ J ijx-T)ajT,T)dT.

t (x)

Если на растущей поверхности Г = Г (t) задано внешнее давление p(t), а на нерас-тущей r = R, - внутреннее р4 (t), т.е. есть

*,(r'(t)) = -p(t), £Tr(R1) = -pi(t),

то радиальная и окружная компоненты напряжений определяются выражениями

г г

J"j(r,t)r2dr + c(t)R

4-dr-p(t), г

г'2 (г2 + R2) , R, (г*' + г)

Г I т

P. (t)-J Jz(r,t)r2dr

Л . А rdr

'r2(r"-R;-) r2(r*2 -R2)l

\ • N ' ' У

r-fr Л1 J r

-J \z{r,t)x4r —dr + —Jz(r,t)r2dr.

r J R|

Здесь функция X(r*t) определяется решением интегро-дифференциального уравнения

t .

Z(t, t) + f„(0)z(r, t) + f f„(r, t - r)z(r,T)dr =

= -E

dep(r,t) ¿r(r,t)-^(r,t) Эг r

в котором ^ = ^ • Е, I > г (г).

В § 5,6 представлены модели объемного и фронтального отверждения сферического изделия в условиях выше изложенных моделей. Приведены результаты численного анализа.

Система для определения радиальной СГ^гД) и тангенциальной £7п(г, О компонент напряжений в процессе объемного отверждения

Ж О + Дд^ т)%{х,г)&т = Р, (г, I),

о

Х(1,1)=^--(<Г„-(Гп), (27)

Эг г

^ ) = 0, от г

где функция Х(г'1) определяется интегральным уравнением

г -

0 +- т)х(г,т)&т = Б, (г, О, о

с принятыми обозначениями

(l-v) dr

а-10~2, кг/см2

Рис. U. Пространственно-временное распределение радиальных егг = Orr(r,t) напряжений; ?(с): 1(132), 2(264), 3(396), 4(528), 5(660), 6(792), 7(924), 8(1320)

Численное решение системы (27) проведено с начальными и граничными условиями вида:

t = 0: сг1Г(г) = О,0'п=О cr„(R1>t)=0) cr„(R,t) = 0.

Динамика распределения радиальных компонент напряжений {рис. 11) указывает на их сжимающий характер на протяжении всего процесса отверждения

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

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

В первом параграфе третьей главы проводится литературный обзор по рассматриваемой главах 3-5 проблеме.

■■ни)

о

г=-ад

Достаточно активно ведутся работы в области математического моделирования СВС-экструзии: разработка тепловых моделей СВС-экструзии и прессования [19-21], разработка неизотермических реодинамических моделей СВС-экструзии, учитывающих реологическое поведение материала и позволяющих изучить не только температурный режим, но и кинетику уплотнения, напряженно-деформированное состояние материала [22].

В работе [23] изучено распределение плотности вязкого пористого тела при одностороннем сжатии в предположении о постоянстве деформаций по всему объему (однородности деформирования); охарактеризованы область использования полученных решений и условия применимости гипотезы об однородности деформирования во взаимосвязи с гипотезой о слабой зависимости свойств тела от координат. При большой вязкости системы критерий Рейнольдса, играющий важную роль в определении характера течения, оказывается малой величиной. Малое значение этого параметра позволяет заменить уравнения движения системы условиями равновесия.

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

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

Рассмотрим (§ 2-4) процесс выдавливания вязко-упругого композитного сжимаемого материала из цилиндрической камеры в направляющий калибр под действием плунжера пресса. Начало координат г = 0 свяжем с центром отверстия основания камеры. Вся область течения разделяется на два отдельных участка: движение внутри камеры между перемещающимся поршнем г = Н(0 и выходным отверстием г = 0+ и течение внутри калибра между г = 0. (входным отверстием в калибр) и свободной поверхностью г = -Ц1), г - время. Возмущениями в обеих областях при переходе из камеры в отверстие пренебрегаем. Движение смеси в каждой из областей считаем одномерным с одной ненулевой компонентой скорости V = V .

Первоначальная длина заготовки равна На, радиус поперечного сечения г0. На дне камеры имеется круглое отверстие радиуса г}, через которое происходит выдавливание материала в направляющий цилиндрический калибр того же радиуса. Ось симметрии заготовки примем в качестве оси г, положительное направление которой противоположно направлению движения поршня {рис. 12).

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

. С ЭУ

а+—а = С—, (28)

/г Эг

Рис. 12. Модель экструзии

a = a(t). (29)

(30)

dt dz

— + У— = Вггт + ф(а'СГУ> (3!)

da т,д а _Э2а — + V— = D—1 dt dz dz

t = 0:p = pü(z),a = aa,a{Q) = a0\ (32)

= = (33)

S0 PlP{0;0 ptpm

da dz

= da

= 0; (34)

=o

Из обобщенной модели Максвелла следует, что в условиях принятого одномерного течения тензор напряжения имеет одну ненулевую компоненту <7.. В силу малости числа Рейнольдса - вязкость несжимаемой основы материала велика - уравнения движения можно заменить условиями равновесия. В рассматриваемом случае течения уравнение равновесия имеет вид Э£7г /Эг = 0. Следовательно, напряжение (Г, в камере не зависит от осевой координаты, а является функцией лишь временной координаты: О, = <7(f). Уравнение (2В) - дифференциальное уравнение состояния Максвелла, где р = ju(a) = р0 схр(ка), G = G(a) = ¡J.(d)¡tr - вязкость и модуль сдвига, соответственно, зависящие от степени структурных превращений а (концентрации межмолекулярных сшивок); (30) - уравнение неразрывности. Здесь р — относительная плотность материала, а плотность среды определяется произведением ррх, где /?, — плотность несжимаемой основы.

Обычно выдавливание материала происходит в матрицу, имеющую коническую форму, которая в данной работе не учитывается. Это допущение справедливо при относительно короткой конической части I «Н0. В случае одномерной модели движение в конической матрице характеризуется двумя характеристиками: относительным изменением плотности и гидравлическим сопротивлением, зависящим от приложенного усилия:

/3(0.,í)/p(0+,t) = В(|сг(0)|), -S,AP(0..í)V(0_1í) = /(|ff(0)|),

где — площадь сечения калибра, p(0_,t), p(0t,t) — плотность, V(0_,í), V(0+,í) — скорость материала на отверстии в калибре и камере, соответственно. Величины В и/определяются формой матрицы и реологическими свойствами материала. В случае отсутствия доуплотнения в калибре В = 1; зависимость /(|<т(0)|) для простоты аппроксимируется степенной зависимостью: /(|<7(0)|) = fcr|cr| .Параметры kr, m могут быть определены на основе опытных данных. С учетом сохранения массы материала получаем скорости (33) на отверстии в камере и калибре, соответственно; S0 — площадь сечения камеры.

Будем трактовать процесс изменения структуры в ходе деформирования, как деформацию связей молекулярно-кинетических единиц под действием приложенного механического поля интенсивностью СТ. Изменение числа межмолекулярных сшивок материала происходит в результате наложения двух процессов: разрушения и образования. По аналогии с задачами химической кинетики суммарную скорость превращения структуры по схеме Ai <-> Аг будем определять по формуле

Ф = -кг ехр(рсг)а + к2(1 - а), где a = a(z',t)~степень структурных превращений, ехр(рсг) -коэффициент скорости разрушения структуры (характеризует снижение эффективной энергии активации), к2 = const — коэффициент скорости накопления числа сшивок. Диффузионно-кинетическое уравнение с учетом, что отличны от нуля лишь осевые составляющие скорости и степени структурных превращений примет вид (31).

Введем обобщенные (лагранжевы) координаты (<7,i,): лагранжево время t, совпадает с реальным временем t, массовая координата q имеет смысл относительной массы материала, находящейся между переменным сечением г и свободной поверхностью г = -Цt). Таким образом,

q = M/S0pl=\p(Z,t)d£ + ^- } PiWl

О -1(1)

Система (28)-(36) с введением новых координат запишется G BV

ол—<j — G—p,cr = cr(t), (37)

¡1 d q

= 0, (38)

Э t д q

1Г + = Орг^ + кг\\-а- ахехр(ра)}; (39)

at dq aq

t = 0: p = p0(q), a = a0(q), cr = <70; (40)

s0 PiPiq ;0 PiPiq ;0

da

T~q

<F1 о

da 'dq

= 0, (42)

■7=9

4=0 =°- («>

В полученной системе ^ = Л, /, д* -элементарная масса, находящаяся на отверстии в рассматриваемый момент времени г.

Система (37)-(43) допускает решение в явном виде относительно скорости течения в камере V = , р — р(дт,/), выраженные через степень структуризации а — а(<7,0-

Для численного нахождения последней используется метод прогонки.

В рамках принятых допущений предполагается отсутствие сопротивления движению материала в калибре: трение о стенки калибра пренебрежимо и, в силу условия (43) на

er -02 -04 -OS -03 -1 -1 2 -1 4

100 120 140

свободной поверхности, напряжение (Т в калибре равно нулю. Это соответствует реальной ситуации, наиболее часто встречающейся на практике в СВС-экструзии тугоплавких материалов, когда процесс уплотнения в основном заканчивается в камере, и происходит выдавливание уже жесткого материала в воздух или направляющий калибр, имеющий диаметр больший, чем диаметр выдавленного стержня. В силу принятых допущений доу-плотнения материала в калибре не происходит, поэтому в области 0<q<q , соответствующей калибру, плотность будет зависеть лишь от массовой координаты д. Каждый бесконечно малый объем выдавленного материала имеет в калибре плотность и степень структуризации, приобретаемые при прохождении отверстия. Выдавленный стержень движется с переменной скоростью того элементарного объема, который в этот момент находится в отверстии. Указанная скорость определяется из закона сопротивления (41) (второе соотношение).

В четвертой главе рассмотрена модель экструзии (37)-(39), (41)-(43) с условием постоянства скорости плунжера пресса: Рассмат-

риваемые начальные условия: ? = 0:

Рис. 13. Изменение напряжения а = сг(г) в камере; V0 = -0.4165-10"4м/с; /¿0 = 9-10'Па-с; а0 = 0; Г, =0.03 с; л = 1/3; А', = 0.105 -10"2 с "

■Ра)ч1чв<

р(д,0) = р0+(рт а(?,0) = 0, ст(0) = 0.

Исследуются вопросы устойчивости течения. Неустойчивый режим течения наблюдался экспериментально при больших значениях напряжения (или при достаточно больших значениях произведения скорости деформации и времени релакса- а ции). В теоретических работах наступление неустойчивого режима течения также связывалось с достижением определенного критического значения величины напряжения. В настоящей работе приведены теоретические результаты, показывающие, что неустойчивый режим течения может наблюдаться как в области малых скоростей, так и в области больших скоростей. Устойчивый же режим течения реализуется лишь в промежуточном интервале значений скоростей. На рис. 13, 14 показано изменение напряжения в камере, степени структуризации в калибре. Из-за малости скорости

О 99995

0 99985

0 99975

О 99365

05

052

0.54

056

05В

Рис. 14. Распределение степени структуризации а = а(д) в калибре; условия на рис 13

выдавливания характерное время выдавливания велико (1600 с), и процессы вязкоупруго-

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

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

Рассмотрим систему уравнений (37)-(43) Предположим, что на плунжере пресса задано постоянное усилие: СГ = СГ0. Введем безразмерные параметры:

4> М0 '

те [0,1];

<7 = т

Р\

Здесь V' —скорость выдавливания в калибр полностью уплотненного материала, ¡х - характерное время процесса. Задача имеет три масштаба характерных времен:

где ^ — характерное время выдавливания, tc — уплотнения, /г — структурирования экс-трудируемого материала. Обозначая

*;='А=1 с44)

— безразмерные времена компактирования (уплотнения), структуризации, выдавливания, соответственно, запишем определяющую систему процесса в виде

дУ . _ .Эр 2ЭГ

гЛр-^ * = = 0,

Эг дд дд Г5

да дд

ч=я 0.

(45)

(46)

(47)

(48)

(49)

Система (45)-(50) допускает полуаналитическое решение: степень структурирования определяется численно с использованием метода прогонки. Проводится численный анализ различных режимов экструдирования материала.

...........1......•/■'-•и-р*-

■—^—-А- -1

-4 -2 0 2 4

Рис. 15 Зависимость предельной структуризации а! от параметра к, = % = 10

Рис. 16. Зависимость предельной плотности рс от параметра кс - ^(\!тс) Зависимость предельной структуризации а, от характерного параметра процесса структуризации к, = 1^(1 / Тг) представлена на рис. 15. Правая ветвь кривой соответствует полному структурированию (разрушению межмолекулярных сшивок) экструдируемого материала. Заметим, что при уменьшении параметра X предельная структуризация а, уменьшается. Левая ветвь на рис. 15 соответствует выдавливанию мало структурированного материала. В этом случае скорость структуризации значительно меньше скорости выдавливания среды. Возрастающий участок кривой соответствует переходным режи-

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

Зависимость предельной плотности рс от характерного параметра уплотнения Кс = представлена на рис. 16: правая ветвь - полное уплотнение материала, ле-

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

Рис. 17. Плоскость характерных времен (т,,^)

Рассмотренные и проанализированные в работе режимы выдавливания, определяемые соотношениями между характерными временами (44), схематично можно изобразить на плоскости (тг,тс): рис. 17.

Область 1 (Г5 « тс «1) - экструзия предельно структурированного и полностью уплотненного материала. В этом случае выдавливание осуществляется в стационарном режиме, при котором степень структуризации и плотность не меняются со временем, то да . др „

есть — = 0, — = 0, причем выдавливание происходит с постоянной скоростью

5/ 9/

V = — 1, равной скорости выдавливания полностью уплотненного материала.

Область 2 (г5 >1, тс «1) - переходные режимы структурирования: выдавливание частично структурированного материала (без структурирования), полностью уплотненного материала.

Область 3 (гс > 1, т1 «1) - переходные режимы уплотнения: выдавливание частично уплотненного материала (без уплотнения), предельно структурированного материала.

Область 4 (тг »1, тс »1) - выдавливание без структурирования и уплотнения: поршневой режим экструзии. Распределение скорости в камере указывает на ее однород-

дУ п

ность по массовой координате: —^ = 0, что следует из уравнения неразрывности в слу-

дд

чае, если — = 0. Следовательно, в данном случае реализуется квазистационарный ре-

3/

жим, при котором плотность не меняется со временем.

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

Основные результаты н выводы

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

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

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

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

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

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

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

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

1. Bulgin D., Wratten R. Measurements of steady flow characteristics of rubber at high rates of shear. Proceedings of the Second International Congress of Rheofogy, Eddited by V.G.W. Harrison. London. 1954.

2. Белкин И.М., Виноградов Г.В., Леонов А.И. Ротационные приборы. Измерение вязкости и физико-химических характеристик. М.: Машиностроение. 1967.220 с.

3. Лифшиц А.Е., Рыбников Г.Л. Диссипативные структуры и течение Куэтта неньютоновской жидкости//ДАН СССР. 1985.Т.281. № 5. С. 1088-1093.

4. Сталин А.М., Худяев С.И., Бучацкий Л.М.. К теории сверханомалии вязкости структурированных систем // ДАН СССР.1978. Т.243. № 2. С. 430-433.

5. Бучацкий Л.М., Манелис Г.Б., Столин A.M., Худяев С.И.. К теории процессов структурных превращений в текучих системах // ИФЖ. 1981. T.XLI. №6. С. 1032-1039.

6. Иванова А.Н., Тарнопольский Б.Л. Устойчивость одномерных стационарных режимов в химических реакторах и бифуркации многомерных режимов // В сб. "Матем. методы в хим. кинетике" под ред. Быкова В.И. Новосибирск. 1990. С. 74-96.

7. Столин A.M., Худяев С.И. Образование пространственно-неоднородных состояний структурированной жидкости в области сверханомалии вязкости // Докл. АН СССР. 1981. Т.260.№ 5. С. 1180-1184.

8. Худяев С.И., Ушаковский О.В. Пространственная неоднородность и автоколебания при течении структурированной жидкости // Матем. моделирование. 2002. Т. 14. № 7. С. 53-73.

9. Томашевский В.Т., Шалыгин В.Н., Яковлев В С. Моделирование условий возникновения технологических дефектов в структуре намоточных композитных полимерных материалов // Механика композитных, материалов. 1980. № 5. С. 895-900.

10. Коротков В.Н., Турусов P.A., Джавадан Э.А., Розенберг Б.А. Технологические напряжения при отверждении цилиндрических изделий из полимерных композитных материалов // Механика композитных материалов. 1986. № 1. С. 118-123.

11. Турусов P.A., Метлов В.В. Формирование напряжений при фронтальном отверждении композитов // Механика композитных материалов. 1985. № б. С. 1079-1085.

12. Клычников Л.В., Давтян С.П., Беляева H.A., Худяев С.И., Ениколопян Н.С. Об остаточных напряжениях при фронтальном отверждении сферического образца с учетом физико-механических явлений перед фронтом отверждения // Механика композитных материалов. 1985. № 3. С. 492-496.

13. Холпанов Л.П., Закиев С.Е., Помогайло А.Д. Моделирвание фронтальной полимеризации с учетом фазовых переходов во фронте реакции // ДАН, 2004. Т. 395. № 2. С. 211-214.

14. Арутюнян Н.Х., Дроздов А.Д., Наумов В.Э. Механика растущих вязко-упруго-пластичесих тел. М.: Наука 1987. 471 с.

15. Staydinger Н., Signer R., Mig G., Hegnstenberg I. Die polimer formaldehyde ein Modell der Zellulose //J. Phys. Chem. 1927. Vol 126. P. 425-429.

16. BartczakZ. //Appl. Polym. Sei. 2002. V. 86. P. 1396.

17. Ошмян В.Г., Патлажан С.A., Remond Y. Принципы структурно-механического моделирования полимеров и композитов // Высокомолек. соед., Серия А, 2006. Т. 48. № 9. С.1-12.

18 Матвеенко В.П, Шардаков И.Н., Бегишев В.П., Писцов Н.В. Simulation ofther-momechanical processes in crystallizing polymer // Int. J. Polymer Engineering and Science. 1997. V. 37. № 8. P. 1270-1279.

19. Стельмах Jl.C., Столин A.M. Математическое моделирование СВС-экструзии. 4.2. Реодинамические модели // ИФЖ. 1993. Т. 64. № 3. С. 90-94.

20. Стельмах Л.С., Столин A.M., Мержанов А.Г. Макрореологическая теория СВС-компактирования // Докл. РАН. 1995. №1. Т. 344. С. 72-77.

21. Stolin A.M., Stelmakh L.S , Zhilyaeva N.N. and Khusid B.M. The phenomenological theory of high-temperature deformation of self-propagating high-temperature synthesis (SHS) products // J. of Materials Synthesis and Processing. 1995. V. 3. №. 1. P. 19-23.

22. Стельмах Л.С., Столин A.M. О квазистационарном режиме и предельных случаях горячей экструзии порошковых материалов // Докл. РАН. 1992. Т. 322. №4. С. 732-736.

23. Бучацкий Л.М., Столин A.M., Худяев С.И. Кинетика изменения распределения плотности при горячем прессовании вязкого пористого тела // Порошковая металлургия. 1986. № 2. С. 37-42.

24. Баронин Г.С., Кербер М..Л., Радько Ю.М., Минкин Е.В. Переработка полимеров в твердой фазе. Физико-химические основы. М.: Машиностроение-1. 2002. 320 с.

25 Сверхмодульные полимеры: Пер. с англ. / Под ред. Чифери А., Уорда И., Малки-на А.Я. Л.: Химия. 272 с.

26. Славнов Е.В., Скульский О.И., Тимофеев В.М. Ориентационные эффекты в образцах после твердофазной гидроэкструзии при высокотемпературном нагружении... // Пластические массы. 1997. № 8 С.18-22.

27. Стельмах Л.С., Столин A.M., Хусид Б.М. Реодинамика выдавливания вязких сжимаемых материалов // ИФЖ. 1991. Т. 61. № 2. С. 268-276.

Список публикаций автора по теме диссертации

1. Клычников Л.В, Давтян С.П., Беляева H.A. Об остаточных напряжениях при фронтальном отверждении сферического образца с учетом физико-механических явлений перед фронтом отверждения // 1-й Всесоюз. симпозиум по макроскопической кинетике и газовой динамике: Тез. докл. Алма-Ата, 1984. С. 66-67.

2. Клычников Л.В., Давтян С.П., Беляева H.A., Худяев С.И., Ениколопян Н.С. Об остаточных напряжениях при фронтальном отверждении сферического образца с учетом физико-механических явлений перед фронтом отверждения // Механика композитных материалов. 1985. № 3. С. 492-496.

3. Клычников Л.В., Беляева H.A., Давтян С.П., Худяев С.И. Распределение технологических внутренних напряжений для цилиндрического образца с учетом кинетики отверждения // III Всесоюзный симпозиум "Теория механической переработки полимерных материалов": Тез. докл. Пермь, 1985. С. 25.

4. Беляева H.A., Клычников Л.В., Давтян С.П., Худяев С.И. Численный анализ напряженного состояния в цилиндрических изделиях, полученных фронтальным отверждением: Препринт ОИХФ АН СССР. Черноголовка, 1987. 36 с.

5. Беляева H.A. Моделирование процесса накопления внутренних напряжений при термической обработке древесины // Республиканская конф. Коми НЦ Уро АН СССР: Тез. докл. Сыктывкар, 1987. С.101.

6. Беляева H.A. Анализ напряженного состояния изделия цилиндрической формы // Республиканская конф. "Разработка и внедрение средств автоматизации в сферу науки, пр-ва, управления": Тез. докл. Сыктывкар, Коми НЦ Уро АН СССР, 1988. С.18.

7. Беляева H.A., Клычников Л.В., Давтян С.П., Худяев С.И. Образование внутренних напряжений в процессе получения изделий цилиндрической формы // Механика композитных материалов. 1988. № 6. С. 1060-1068.

8. Беляева H.A., Клычников J1.B., Давтян С.П., Худяев С.И. Анализ напряженного состояния в процессе формирования цилиндрических изделий с учетом совместного протекания полимеризации и кристаллизации // Ш Всесоюз. конф. по нелинейной теории упругости: Тез. докл. Сыктывкар, ¡989.

9. Беляева H.A. Исследование влияния неизотермических условий синтеза на некоторые механические свойства полимеров: Дис. канд. физ,-мат. наук. Черноголовка, 1990. 164 с.

10. Беляева H.A. Исследование влияния неизотермических условий синтеза на некоторые механические свойства полимеров: Автореферат дис.канд. физ -мат. наук. Черноголовка, 1990.22 с.

П. Беляева H.A. Совмещенные процессы и напряженное состояние формируемого цилиндрического изделия // XI Коми Респ. молодежная научная конференция: Тез. докл. Сыктывкар, 1990.

12. Беляева H.A., Клычников Л В., Давтян С.П., Худяев С.И. Напряженное состояние при формировании цилиндрических изделий в ходе параллельного протекания полимеризации и кристаллизации // Механика композитных материалов. 1991. № 6. С. 10911099.

13. Беляева H.A., Клычников Л.В., Давтян С.П., Худяев С И. Особенности кинетики отверждения и напряженное состояние при формировании симметричных изделий в условиях совместного протекания полимеризации и кристаллизации // Материалы Всерос. конф. Переработка полимерных материалов в изделиях. Ижевск: Ижевский механический ин-т, 16-19 ноября 1993. С.19-21.

14. Беляева H.A., Клычников Л.В. Численный анализ кинетики отверждения и напряженного состояния в цилиндрических изделиях, полученных в условиях совместного протекания полимеризации и кристаллизации // 8-я Всерос. конф. Математические методы в химии (ММХ-8): Тез. докл. Тула: Тульский гос. технич. ун-т, 27-29 окт. 1993. С.112.

15. Беляева H.A., Худяев С.И. Регулирование процесса формирования цилиндрического изделия в ходе отверждения // 2-я международная научно-техническая конференция "Актуальные проблемы фундаментальных наук": Материалы. Москва, январь 1994. Т.5.4 с.

16. Беляева H.A., Беляев Ю.Н. Исследование напряженного состояния в ходе сушки древесины // 2-я международная научно-техническая конференция "Актуальные проблемы фундаментальных наук": Материалы. Москва, январь 1994. Т. 1. С. 74-76.

17. Беляева H.A., Беляев Ю.Н. Формирование цилиндрического изделия в ходе регулируемого процесса отверждения // Симпозиум "Математические модели деформируемого тела": Материалы. СПб, июнь 1994. 8 с.

18. Беляева H.A., Беляев Ю.Н. Регулирование уровня внутренних напряжений формируемого изделия в ходе совмещенного процесса полимеризации и кристаллизации // Вестник Сыктывкарского университета. Вып. 1. 1995. С. 119-127.

19. Беляева H.A., Клычников Л.В. Метод интегрального уравнения в задаче объемного отверждения // Вестник Сыктывкарского университета. Вып.2. 1996. С.125-134.

20. Беляева H.A., Титова H.H., Юркина Е.А. Модель фронтального отверждения цилиндрического изделия с точки зрения механики непрерывно наращиваемого твердого тела // Труды Сыктывкарского лесного института. Т. 1. 1997. С. 44-52.

21. Беляева H.A., Паршукова Н.М. Термовязкоупругая модель отверждения сферического изделия // Вестник Сыктывкарского университета. Серия 1. Вып.4. 2001. С.115-124.

22. Беляева H.A., Худяев С.И., Д.Л.Горст Д.Л. Неоднородное течение Куэтга структурированной жидкости // Вестник Сыктывкарского ун-та. Сер. 1. Вып.5.2002. С. 43-48.

23. Беляева H.A., Беляев Д.Ю., Горст Д.Л. Куэтговское течение структурированной жидкости // Минобразования РФ. ОФАП. Свид. об отрасл. Регистрации разработки № 3462. 2004.

24. Беляева H.A., Беляев Д.Ю. Неоднородные структуры при течении Куэтга псевдопластической жидкости // Минобразования РФ. ОФАП. Свид. об отрасл. регистрации разработки № 3463. 2004.

25. Беляева H.A. Пространственно неоднородные течения структурированной жидкости в плоском зазоре // Нелинейные проблемы механики и физики деформируемого твердого тела: Тр. научн школы акад. В.В. Новожилова. СПб.: СПбГУ, 2004. Вып.8. С. 174-190.

26. Беляева H.A., Вербовик А.Н. Напорное течение структурированной жидкости // Нелинейные проблемы механики и физики деформируемого твердого тела: Тр. научн. школы акад. В.В. Новожилова. СПб.: СПбГУ, 2005. Вып.8. С. 151-156.

27. Беляева H.A., Размыслов Р.Ю. Сдвиговое течение структурированной жидкости // Нелинейные проблемы механики и физики деформируемого твердого тела: Тр. научн. школы акад. В.В. Новожилова СПб.: СПбГУ, 2005. Вып.8. С. 186-193.

28. Беляева H.A. Неоднородное течение структурированной жидкости // Математическое моделирование. 2006. Т. 18. №6. С. 3-14.

29. Беляева Н.А Математические модели течения структурированной жидкости // Федеральное агентство по образованию РФ. ОФАП. Свид. об отрасл. регистрации разработки № 5699. 01.03. 2006.

30. Беляева H.A. Самоорганизация куэтговского течения структурированной жидкости // Материалы III Международной летней научной школы. Кемерово, Россия. 22-28 июня 2006 г. С.321-326.

31. Беляева H.A., Столин А.М., Стельмах Л.С., Смолев Л.В. Численное моделирование твердофазной экструзии вязкоупругого структурированного материала // Четвертая Всероссийская школа-семинар по структурной макрокинетике: Тез. докл. 22-25 ноября 2006 г. Черноголовка: ИСМАН, 2006. С. 86-87.

32. Беляева H.A. Вязкоупругая модель отверждения композитного изделия // Труды XXXVI Уральского семинара "Механика и процессы управления". Екатеринбург, 2006. Т. 1.С. 22S-235.

33. Беляева H.A., Осипова В.В. Формирование цилиндрического изделия в ходе объемного отверждения // Федеральное агентство по образованию. ОФАП. Свид. об отрасл. регистрации разработки. № 7944.30 03 2007.

34. Беляева H.A., Смолев Л.В. Экструзия с заданным усилием на плунжере пресса // Федеральное агентство по образованию ОФАП. Свид. об отрасл. регистрации разработки № 7945. 30.03 2007.

35. Беляева H.A. Твердофазная экструзия с условием постоянства скорости плунжера пресса // Федеральное агентство по образованию. ОФАП. Свид. об отрасл. регистрации разработки № 7946. 30.03 2007.

36. Столин A.M., Стельмах Л С., Беляева H.A. Структурные теории в реологии // Программа и тезисы докладов конференции "Реология и физико-химическая механика гетерофазных систем". Карачарово, 23-28 апреля 2007 г.

37. Беляева H.A., Столин A.M., Стельмах U.C. Кинетика уплотнения и структуризации в твердофазной экструзии вязкоупругой среды // Инженерная физика. 2007. № 5. С. 34-41.

38. Беляева H.A., Столин A.M., Стельмах Л.С. Влияние структуризации на процесс твердофазной экструзии вязкоупругого материала // Материалы Международной конференции. XVIII сессия Международной школы по моделям механики сплошной среды. Саратов: Саратовский госуниверситет. 2007. С. 55-58.

39. Беляева H.A. Математические модели деформируемых структурированных материалов. Монография. Сыктывкар: Изд-во СыктГУ, 2008.116 с.

40. Беляева H.A. Структуризация в процессах деформирования вязкоупругих систем // Вестник Удмуртского университета. Математика. Механика Компьютерные науки. Ижевск: Удмуртский госуниверситет. 2008. Вып. 2. С. 178-179.

41. Беляева H.A., Столин A.M., Пугачев Д.В., Стельмах Л.С. Неустойчивые режимы деформирования при твердофазной экструзии вязкоупругих структурированных систем // ДАН, 2008. Т. 420. № 6. С. 777-780.

42. Беляева H.A., Столин A.M., Стельмах Л.С. Неустойчивость в твердофазной экструзии вязкоупругого структурированного материала // Труды пятой Всероссийской научной конференции с международным участием. Математическое моделирование и краевые задачи. 29-31 мая 2008 г. Самара: Самарский технический университет. Ч. 1. С. 56-58.

43. Беляева H.A., Столин А.М., Стельмах Л.С. Динамика твердофазной плунжерной экструзии вязкоупругого структурированного материала // Теоретические основы химической технологии, 2008. Ка 5. С. 579-589.

44. Belyaeva N. A., Stolin A. M., Stelmakh L. S. Dynamic of Solid-State Extrusion of Viscoelastic Cross-Linked polymeric Materials // Theoretical Foundations of Chemical Engineering, 2008. Vol. 42. No 5. P. 549-556.

45. Беляева H.A. Структурные модели процессов деформирования вязкоупругих материалов // Вестник Сыктывкарского университета. Сер 1. Вып. 8. 2008. С. 55-68.

ИПО СыктГУ. Усл.-печ.л. 2,0. Тираж 100 экз Заказ №338.

ВВЕДЕНИЕ.

ГЛАВА 1. МАТЕМАТИЧЕСКИЕ МОДЕЛИ ТЕЧЕНИЯ

СТРУКТУРИРОВАННОЙ ЖИДКОСТИ.

1. Модель структурированной жидкости.

2. Куэттовское течение в плоском зазоре.

3. Стационарные однородные режимы течения.

4. Псевдопластическая жидкость.

5. Дилатантная жидкость.

6. Сдвиговое течение.

7. Напорное течение.

ГЛАВА 2. ТЕРМОВЯЗКОУПРУГИЕ МОДЕЛИ ОТВЕРЖДЕНИЯ

СИММЕТРИЧНЫХ ИЗДЕЛИЙ.

1. Обзор литературы по проблемам получения композитных изделий методом отверждения.

2. Макрокинетическая модель совмещенного процесса.

3. Применение метода интегрального уравнения в задаче объемного отверждения цилиндрического изделия.

4. Фронтальное отверждение цилиндрического изделия.

5. Объемное отверждение сферического изделия.

6. Фронтальное формирование полого сферического изделия.

ГЛАВА 3. МОДЕЛЬ ТВЕРДОФАЗНОЙ ПЛУНЖЕРНОЙ ЭКСТРУЗИИ ВЯЗКОУПРУГИХ СТРУКТУРИРОВАННЫХ СИСТЕМ.

1. Литературный обзор по вопросам экструзии.

2. Постановка задачи.

3. Переход к обобщенным координатам.

4. Определение параметров течения в камере.

5. Параметры течения в калибре.

6. Результаты численного эксперимента.

ГЛАВА 4. ЭКСТРУЗИЯ С ЗАДАННОЙ СКОРОСТЬЮ ПЛУНЖЕРА.

1. Математическая модель.

2. Решение системы определяющих уравнений.

3. Определение параметров течения в калибре.

4. Численное решение.

5. Анализ результатов.

6. Неустойчивые режимы деформирования.

7. Новые экспериментальные данные.

ГЛАВА 5. РЕЖИМЫ ТВЕРДОФАЗНОЙ ЭКСТРУЗИИ

ВЯЗКОУПРУГИХ СТРУКТУРИРОВАННЫХ СИСТЕМ.

1. Основные уравнения.

2. Режимы экструзии.

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

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

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

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

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

Работа выполнялась в соответствии с федеральной целевой программой «Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России» на 2007 - 2012 гг. при выполнении госконтракта 02.513.11.3377 от 26 ноября 2007 г.

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

1. Исследование процессов самоорганизации течения вязкой структурированной жидкости в области "сверханомалии вязкости".

2. Построение вязкоупругой модели объемного и фронтального отверждения осесимметричных изделий (цилиндр, сфера) в ходе совместного протекания реакций полимеризации и кристаллизации. Определение температурного и конверсионных полей в процессе фазового перехода на основе макрокинетической модели совмещенного процесса полимеризации и кристаллизации.

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

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

5. Численный анализ качественно различных режимов твердофазной экструзии полимерных материалов. Нахождение областей реализации этих режимов в зависимости от характерных времен процесса.

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

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

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

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

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

• режим объемного отверждения цилиндрического изделия;

• режим фронтального отверждения цилиндрического изделия;

• режим объемного отверждения сферического изделия;

• режим фронтального отверждения сферического изделия.

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

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

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

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

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

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

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

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

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

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

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

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

Результаты теоретических исследований использовались в Научно-образовательном центре ТамбГТУ-ИСМАН «Твердофазные технологии» при получении изделий методами твердофазной экструзии из различных полимерных композитов. Получено заключение о практическом использовании основных положений диссертации в Институте структурной макрокинетики и проблем материаловедения РАН (г. Черноголовка) и Тамбовском государственном техническом университете (г. Тамбов). Личный вклад соискателя. Все приведенные в диссертации результаты получены лично автором или при его непосредственном участии. Монография «Математические модели деформируемых структурированных материалов», отражающая главные направления диссертации, опубликована без соавторов. В постановке задач по отверждению цилиндрического изделия, обсуждении результатов принимали участие д.ф.-м.н. Худяев С.И., д.ф.-м.н. Давтян С.П., к.ф.-м.н. Клычников JI.B. В постановке задачи куэттовского течения структурированной жидкости принимал участие д.ф.-м.н. Худяев С.И.; в постановке и обсуждении результатов исследования процесса твердофазной экструзии - д.ф.-м.н., профессор Столин A.M., д.т.н. Стельмах JI.C. При подготовке статьи: Беляева Н.А., Столин A.M., Пугачев Д.В., Стельмах JI.C. Неустойчивые режимы деформирования при твердофазной экструзии вязкоупругих структурированных систем // ДАН, 2008. Т. 420. № 6. С. 777-780. - соавтором Пугачевым Д.В. был поставлен эксперимент, результаты которого используются в названной публикации.

Апробация работы. Результаты диссертационной работы докладывались на: Всероссийской конференции «Переработка полимерных материалов в изделиях» (Ижевск, 16-19 ноября 1993 г.), 8-й Всероссийской конференции «Математические методы в химии» (Тула, 27-29 окт. 1993 г.), 2-й Международной научно-технической конференции "Актуальные проблемы фундаментальных наук" (Москва, январь 1994 г.), Симпозиуме

Математические модели деформируемого тела» (СПб, июнь 1994 г.), XXXVI Уральском семинаре «Механика и процессы управления» (Екатеринбург, 2006 г.), III Международной летней научной школе «Гидродинамика больших скоростей и численное моделирование» (Кемерово, 22-28 июня 2006 г.), Четвертой Всероссийской школе-семинаре по структурной макрокинетике для молодых ученых (Черноголовка ИСМАН , 22-25 ноября 2006 г.), конференции «Реология и физико-химическая механика гетерофазных систем» (Карачарово, 23 — 28 апреля 2007 г.), Международной конференции «XVIII сессия Международной Школы по моделям механики сплошной среды» (Саратов, 2007 г.), научной конференции-семинаре «Теория управления и математическое моделирование» (Ижевск, 2008 г.), пятой Всероссийской конференции с международным участием «Математическое моделирование и краевые задачи» (Самара, 2008 г.); научном семинаре кафедры теории упругости механико-математического факультета МГУ (Москва, 1999 г.), научном семинаре в Отделе математики Коми научного центра УрО РАН (Сыктывкар, 2004 г.), научном семинаре лаборатории гидродинамики математико-механического ф-та СПб технического университета (СПб, 2004) , научном семинаре ИСМАН РАН (г. Черноголовка, 2008 г.), научном семинаре Института механики сплошной среды (ИМСС) РАН (Пермь, 2008 г.). Публикации. Всего по теме диссертации опубликовано 45 научных работы, из них 9 научных работ - в ведущих журналах, рекомендованных ВАК. Структура и объем диссертации. Диссертация состоит из введения, пяти глав, заключения, библиографического списка из 191 наименования, приложения; содержит 193 страницы текста, 102 рисунка, 6 таблиц. Основные положения, выносимые на защиту:

ЗАКЛЮЧЕНИЕ

Отметим основные результаты и выводы, следующие из исследований, представленных в данной диссертационной работе:

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

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

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

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

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

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

СПИСОК ПУБЛИКАЦИЙ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ

1. Клычников Л.В., Давтян С.П., Беляева Н.А. Об остаточных напряжениях при фронтальном отверждении сферического образца с учетом физико-механических явлений перед фронтом отверждения // 1-й Всесоюз. симпозиум по макроскопической кинетике и газовой динамике: Тез. докл. Алма-Ата, 1984. С. 66-67.

2. Клычников Л.В., Давтян С.П., Беляева Н.А., Худяев С.И., Ениколопян Н.С. Об остаточных напряжениях при фронтальном отверждении сферического образца с учетом физико-механических явлений перед фронтом отверждения // Механика композитных материалов. 1985. № 3. С. 492-496.

3. Клычников Л.В., Беляева Н.А., Давтян С.П., Худяев С.И. Распределение технологических внутренних напряжений для цилиндрического образца с учетом кинетики отверждения // III Всесоюзный симпозиум "Теория механической переработки полимерных материалов": Тез. докл. Пермь, 1985. С. 25.

4. Беляева Н.А., Клычников Л.В., Давтян С.П., Худяев С.И. Численный анализ напряженного состояния в цилиндрических изделиях, полученных фронтальным отверждением: Препринт ОИХФ АН СССР. Черноголовка, 1987. 36 с.

5. Беляева Н.А. Моделирование процесса накопления внутренних напряжений при термической обработке древесины // Республиканская конф. Коми НЦ Уро АН СССР: Тез. докл. Сыктывкар, 1987. С.101.

6. Беляева Н.А. Анализ напряженного состояния изделия цилиндрической формы // Республиканская конф. "Разработка и внедрение средств автоматизации в сферу науки, пр-ва, управления": Тез. докл. Сыктывкар, Коми НЦ Уро АН СССР, 1988. С. 18.

7. Беляева Н.А., Клычников Л.В., Давтян С.П., Худяев С.И. Образование внутренних напряжений в процессе получения изделий цилиндрической формы // Механика композитных материалов. 1988. № 6. С.1060-1068.

8. Беляева Н.А., Клычников Л.В., Давтян С.П., Худяев С.И. Анализ напряженного состояния в процессе формирования цилиндрических изделий с учетом совместного протекания полимеризации и кристаллизации // III Всесоюз. конф. по нелинейной теории упругости: Тез. докл. Сыктывкар, 1989.

9. Беляева Н.А. Исследование влияния неизотермических условий синтеза на некоторые механические свойства полимеров: Дис. канд. физ.-мат. наук. Черноголовка, 1990. 164 с.

Ю.Беляева Н.А. Исследование влияния неизотермических условий синтеза на некоторые механические свойства полимеров: Автореферат дис.канд. физ.-мат. наук. Черноголовка, 1990. 22 с.

11 .Беляева Н.А. Совмещенные процессы и напряженное состояние формируемого цилиндрического изделия // XI Коми Респ. молодежная научная конференция: Тез. докл. Сыктывкар, 1990.

12.Беляева Н.А., Клычников Л.В., Давтян С.П., Худяев С.И. Напряженное состояние при формировании цилиндрических изделий в ходе параллельного протекания полимеризации и кристаллизации // Механика композитных материалов. 1991. № 6. С. 1091-1099.

13.Беляева Н.А., Клычников Л.В., Давтян С.П., Худяев С.И. Особенности кинетики отверждения и напряженное состояние при формировании симметричных изделий в условиях совместного протекания полимеризации и кристаллизации // Материалы Всерос. конф. Переработка полимерных материалов в изделиях. Ижевск: Ижевский механический ин-т, 16-19 ноября 1993. С.19-21.

Н.Беляева Н.А., Клычников Л.В. Численный анализ кинетики отверждения и напряженного состояния в цилиндрических изделиях, полученных в условиях совместного протекания полимеризации и кристаллизации // 8-я

Всерос. конф. Математические методы в химии (ММХ-8): Тез. докл. Тула: Тульский гос. технич. ун-т, 27-29 окт. 1993. С.112.

15.Беляева Н.А., Худяев С.И. Регулирование процесса формирования цилиндрического изделия в ходе отверждения // 2-я международная научно-техническая конференция "Актуальные проблемы фундаментальных наук": Материалы. Москва, январь 1994. Т.5. 4 с.

16.Беляева Н.А., Беляев Ю.Н. Исследование напряженного состояния в ходе сушки древесины // 2-я международная научно-техническая конференция "Актуальные проблемы фундаментальных наук": Материалы. Москва, январь 1994. Т. 1.С. 74-76.

17.Беляева Н.А., Беляев Ю.Н. Формирование цилиндрического изделия в ходе регулируемого процесса отверждения // Симпозиум "Математические модели деформируемого тела": Материалы. СПб, июнь 1994. 8 с.

18.Беляева Н.А., Беляев Ю.Н. Регулирование уровня внутренних напряжений формируемого изделия в ходе совмещенного процесса полимеризации и кристаллизации // Вестник Сыктывкарского университета. Вып. 1. 1995. С.119-127.

19.Беляева Н.А., Клычников JI.B. Метод интегрального уравнения в задаче объемного отверждения // Вестник Сыктывкарского университета. Вып.2. 1996. С.125-134.

20.Беляева Н.А., Титова Н.Н., Юркина Е.А. Модель фронтального отверждения цилиндрического изделия с точки зрения механики непрерывно наращиваемого твердого тела // Труды Сыктывкарского лесного института. Т. 1. 1997. С. 44-52.

21.Беляева Н.А., Паршукова Н.М. Термовязкоупругая модель отверждения сферического изделия // Вестник Сыктывкарского университета. Серия 1. Вып.4. 2001. С.115-124.

22.Беляева Н.А., Худяев С.И., Д.Л.Горст Д.Л. Неоднородное течение Куэтта структурированной жидкости // Вестник Сыктывкарского ун-та. Сер.1. Вып.5. 2002. С. 43-48.

23.Беляева Н.А., Беляев Д.Ю., Горст Д.Л. Куэттовское течение структурированной жидкости // Минобразования РФ. ОФАП. Свид. об отрасл. Регистрации разработки № 3462. 2004.

24.Беляева Н.А., Беляев Д.Ю. Неоднородные структуры при течении Куэтта псевдопластической жидкости // Минобразования РФ. ОФАП. Свид. об отрасл. регистрации разработки № 3463. 2004.

25.Беляева Н.А. Пространственно неоднородные течения структурированной жидкости в плоском зазоре // Нелинейные проблемы механики и физики деформируемого твердого тела: Тр. научн. школы акад. В.В. Новожилова. СПб.: СПбГУ, 2004. Вып.8. С. 174-190.

26.Беляева Н.А., Вербовик А.Н. Напорное течение структурированной жидкости // Нелинейные проблемы механики и физики деформируемого твердого тела: Тр. научн. школы акад. В.В. Новожилова. СПб.: СПбГУ, 2005. Вып.8. С. 151-156.

27.Беляева Н.А., Размыслов Р.Ю. Сдвиговое течение структурированной жидкости // Нелинейные проблемы механики и физики деформируемого твердого тела: Тр. научн. школы акад. В.В. Новожилова. СПб.: СПбГУ, 2005. Вып.8. С. 186-193.

28.Беляева Н.А. Неоднородное течение структурированной жидкости // Математическое моделирование. 2006. Т. 18. №6. С. 3-14.

29.Беляева Н.А. Математические модели течения структурированной жидкости // Федеральное агентство по образованию РФ. ОФАП. Свид. об отрасл. регистрации разработки № 5699. 01.03. 2006.

30.Беляева Н.А. Самоорганизация куэттовского течения структурированной жидкости // Материалы III Международной летней научной школы. Кемерово, Россия. 22-28 июня 2006 г. С.321-326.

31.Беляева Н.А., Столин A.M., Стельмах Л.С., Смолев Л.В. Численное моделирование твердофазной экструзии вязкоупругого структурированного материала // Четвертая Всероссийская школа-семинар по структурной макрокинетике: Тез. докл. 22-25 ноября 2006 г. Черноголовка: ИСМАН, 2006. С. 86-87.

32.Беляева Н.А. Вязкоупругая модель отверждения композитного изделия // Труды XXXVI Уральского семинара "Механика и процессы управления". Екатеринбург, 2006. Т. 1. С. 228-235.

33.Беляева Н.А., Осипова В.В. Формирование цилиндрического изделия в ходе объемного отверждения // Федеральное агентство по образованию. ОФАП. Свид. об отрасл. регистрации разработки. № 7944. 30.03 2007.

34.Беляева Н.А., Смолев JI.B. Экструзия с заданным усилием на плунжере пресса // Федеральное агентство по образованию. ОФАП. Свид. об отрасл. регистрации разработки № 7945. 30.03 2007.

35.Беляева Н.А. Твердофазная экструзия с условием постоянства скорости плунжера пресса // Федеральное агентство по образованию. ОФАП. Свид. об отрасл. регистрации разработки № 7946. 30.03 2007.

36.Столин A.M., Стельмах JI.C., Беляева Н.А. Структурные теории в реологии // Программа и тезисы докладов конференции "Реология и физико-химическая механика гетерофазных систем". Карачарово, 23-28 апреля 2007 г.

37.Беляева Н.А., Столин A.M., Стельмах JI.C. Кинетика уплотнения и структуризации в твердофазной экструзии вязкоупругой среды // Инженерная физика. 2007. № 5. С. 34-41.

38.Беляева Н.А., Столин A.M., Стельмах JI.C. Влияние структуризации на процесс твердофазной экструзии вязкоупругого материала // Материалы Международной конференции. XVIII сессия Международной школы по моделям механики сплошной среды. Саратов: Саратовский госуниверситет.

2007. С. 55-58.

39.Беляева Н.А. Математические модели деформируемых структурированных материалов. Монография. Сыктывкар: Изд-во СыктГУ,

2008. 116 с.

40.Беляева Н.А. Структуризация в процессах деформирования вязкоупругих систем // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. Ижевск: Удмуртский госуниверситет. 2008. Вып. 2. С. 178-179.

41.Беляева Н.А., Столин A.M., Пугачев Д.В., Стельмах JI.C. Неустойчивые режимы деформирования при твердофазной экструзии вязкоупругих структурированных систем // ДАН, 2008. Т. 420. № 6. С. 777-780.

42.Беляева Н.А., Столин A.M., Стельмах JI.C. Неустойчивость в твердофазной экструзии вязкоупругого структурированного материала // Труды пятой Всероссийской научной конференции с международным участием. Математическое моделирование и краевые задачи. 29-31 мая 2008 г. Самара: Самарский технический университет. Ч. 1. С. 56-58.

43.Беляева Н.А., Столин A.M., Стельмах JI.C. Динамика твердофазной плунжерной экструзии вязкоупругого структурированного материала // Теоретические основы химической технологии, 2008. № 5. С. 579-589.

44.Belyaeva N. A., Stolin А. М., Stelmakh L. S. Dynamic of Solid-State Extrusion of Viscoelastic Cross-Linked polymeric Materials // Theoretical Foundations of Chemical Engineering, 2008. Vol. 42. No 5. P. 549-556.

45.Беляева Н.А. Структурные модели процессов деформирования вязкоупругих материалов // Вестник Сыктывкарского университета. Сер 1. Вып. 8. 2008. С. 55-68.

1. Астарита Дж., Маруччи Дж. Основы гидромеханики неньютоновских жидкостей. М.: Мир. 1978.

2. Лойцянский Л.Г. Механика жидкости и газа. М.: Наука. 1978. 736 с.

3. Николис Д., Пригожин И. Самоорганизация в неравновесных системах. М.: Мир. 1979.

4. Bulgin D., Wratten R. Measurements of steady flow characteristics of rubber at high rates of shear. Proceedings of the Second International Congress of Rheology, Eddited by V.G.W. Harrison. London. 1954.

5. Белкин И.М., Виноградов Г.В., Леонов А.И. Ротационные приборы. Измерение вязкости и физико-химических характеристик. М.: Машиностроение. 1967. 220 с.

6. Петржик Г.Г., Трапезников А.А. Дилатансия паст двуокиси титана в растворах алкидных смол в зависимости от состава и температуры // Коллоидный журнал. 1975. Т. 37. С. 1193-1197.

7. Лифшиц А.Е., Рыбников Г.Л. Диссипативные структуры и течение Куэтта неньютоновской жидкости // ДАН СССР. 1985.Т.281. № 5. С. 1088-1093.

8. Столин A.M., Стельмах Л.С., Беляева Н.А. Структурные теории в реологии/ЯТрограмма и тезисы докладов. Конференция "Реология и физико-химическая механика гетерофазных систем". Карачарово, 23 28 апреля 2007 г. С. 35-36.

9. Eyring Н. Viscosity, Plasticity and Diffusion as Examples of Absolute Reaction Rate // J. Chem. Phys. 1936. V. 4. P. 283 291.

10. Бучацкий Л.М., Манелис Г.Б., Столин A.M., Худяев С.И. К теории процессов структурных превращений в текучих системах // ИФЖ. 1981. T.XLI. №6. С. 1032-1039.

11. З.Иванова А.Н., Тарнопольский Б. Л. Устойчивость одномерных стационарных режимов в химических реакторах и бифуркации многомерных режимов // В сб. "Матем. методы в хим. кинетике" под ред. Быкова В.И. Новосибирск. 1990. С. 74-96.

12. Столин A.M., Худяев С.И. Образование пространственно-неоднородных состояний структурированной жидкости в области сверханомалии вязкости // Докл. АН СССР.1981. Т.260. № 5. С. 1180-1184.

13. Худяев С.И., Ушаковский О.В. Пространственная неоднородность и автоколебания при течении структурированной жидкости . // Матем. моделирование. 2002. Т. 14. № 7. С. 53-73.

14. Худяев С.И. Пороговые явления в нелинейных уравнениях. М.: Физматлит.2003. 267 с.

15. Беляева Н.А., Горст Д.Л., Худяев С.И. Неоднородное течение Куэтта структурированной жидкости // Вестник СыктГУ. 2003. Вып. 5. С. 43 48.

16. Беляева Н.А., Беляев Д.Ю., Горст Д.Н. Куэттовское течение структурированной жидкости // Минобразования РФ. ОФАП. Свид. об отрасл. регистрации разработки № 3462. 2004.

17. Беляева Н.А., Беляев Д.Ю. Неоднородные структуры при течении Куэтта псевдопластической жидкости // Минобразования РФ. ОФАП. Свид. об отрасл. регистрации разработки № 3463. 2004.

18. Беляева Н.А. Пространственно-неоднородные течения структурированной жидкости в плоском зазоре // Сб. Нелинейные проблемы механики и физики деформируемого твердого тела. Санкт-Петербург. 2004. С. 174-190.

19. Беляева Н.А., Вербовик А.Н. Напорное течение структурированной жидкости // В сб.: Нелинейные проблемы механики и физикидеформируемого твердого тела. Тр. научн. школы акад. В.В. Новожилова. СПб.: СПбГУ, Вып. 8. 2005. С. 151-156.

20. Беляева Н.А., Размыслов Р.Ю. Сдвиговое течение структурированной жидкости // В сб.: Нелинейные проблемы механики и физики деформируемого твердого тела. Тр. научн. школы акад. В.В. Новожилова. СПб.: СПбГУ, Вып. 8. 2005. С. 186-193.

21. Беляева Н.А. Неоднородное течение структурированной жидкости // Математическое моделирование. Т. 18. № 6. 2006. С. 3 -14.

22. Беляева Н.А. Математические модели течения структурированной жидкости. Федеральное агентство по образованию РФ. ОФАП. Свид. об отрасл. регистрации разработки № 5699. 01.03. 2006.

23. Беляева Н.А. Самоорганизация куэттовского течения структурированной жидкости // Материалы III Международной летней научной школы. Кемерово, Россия. 22-28 июня 2006 г. С. 321-326.

24. Ландау Л.Д., Лифшиц Е.М. Гидродинамика. М.: Наука. 1988. 736 с.

25. Самарский А.А. Теория разностных схем. М: Наука. 1983. 616 с.

26. Зельдович Я.Б. К теории распространения пламени // Ж. физ. химии. 1948. Т.22. № 1. С. 27-48.

27. Зельдович Я.Б., Франк-Каменецкий Д.А. К теории равномерного распространения пламени // Докл. АН СССР. 1938. Т. 19. С. 693-697.

28. Семенов Н.Н. О некоторых проблемах химической кинетики и реакционной способности. М.: Изд-во АН СССР. 1954. 350 с.

29. Максимов Э.И. О закономерностях полимеризации в периодическом реакторе идеального смешения // ТОХТ. 1973. № 3. С. 380-385.

30. Мержанов А.Г., Абрамов Б.Г., Абрамова Л.Т. Термографический метод исследования кинетики тепловыделения // Ж. физ. химии. 1967. Т.41. № 1. С. 173-184.

31. Вольперт А.И., Худяев С.И. Анализ в классах разрывных функций и уравнения математической физики. М.: Наука. 1975. 394 с.

32. Жирков П.В. Исследование неизотермических течений реагирующих жидкостей с изменяющейся в процессе химической реакции вязкостью. Автореферат дисс. на соискание ученой степени кандидата физико-математичесих наук. Черноголовка. 1979. 16 с.

33. Боярченко В.И. Макрокинетическая теория экструзии полимерных и полимеризующихся материалов. Автореферат дисс. на соискание ученой степени доктора физико-математичесих наук. Черноголовка. 1982. 30 с.

34. Вундерлих Б. Физика макромолекул. Пер. с англ. / Под. ред. Годовского Ю.К. М.: Мир. 1979. Т. 2. 572 с.

35. Вольф Л.А., Хайтин Б.Ш. Производство поликапроамида. М.: Химия.1977. 206 с.

36. Воробьев В.А., Андрианов Р.А. Технология полимеров. М.: Высшая школа. 1980. 303 с.

37. Нелсон У.Е. Технология пластмасс на основе полиамидов. Пер. с англ. / Под. ред. Малкина А.Я. М.: Химия. 1979. 256 с.

38. Манделькерн Л. Кристаллизация полимеров. Пер. с англ. / Под. ред. Френкеля С .Я. М.: Химия. 1966. 336 с.

39. Lindenmeyer Р.Н. Crystalliztionin polymers // J. of Polymer Science. Part C. 1963. № l.P. 3-39.

40. Цахман Г. Кристаллизация и плавление полимеров // Химия и технология полимеров. 1966. № 5. С. 3-77.

41. Шарплез А. Кристаллизация полимеров. Пер. с англ. / Под. ред. Бакеева Н.Ф. М.: Мир. 1968. 200 с.

42. Avrami М. Kinetiks of phase change. I. General Theory // J. Chem. Phys. 1939. V. 7. № 3. P. 1103-1112.

43. Avrami M. Kinetiks of phase change. II. Transformation time relations for randon distribution of nuclei // J. Chem. Phys. 1950. V. 8. № 1. P. 212-224.

44. Avrami M. Granulation, phase change and microstructure. Kinetiks of phase change. Ill //J. Chem. Phys. 1941. V. 9. № 1. P. 177-184.

45. Годовский Ю.К. Теплофизические методы исследования полимеров. М.: Химия. 1976. 216 с.

46. James N.H. Crystallization kinetics and metting studies // Brit. Polym. J. 1979. V. 11. №5. P. 137-145.

47. Бегишев В.П., Болгов С.А., Малкин А.Я., Субботина Н.И. , Фролов В.Г. Тепловые эффекты при полимеризации, сопровождающейся кристаллизацией //Высокомол. соединения. 1979. Т. Б.21. № 9. С. 714-718.

48. Бегишев В.П., Кипин И.А., Малкин А.Я. Описание процесса кристаллизации полимеров с помощью макрокинетического уравнении // Высокомол. соединения. 1982. Т. Б.24. № 9. С. 656-658.

49. Бегишев В.П., Кипин И.А., Андрианова 3. С., Малкин А.Я. Кинетика неизотермического процесса кристаллизации поликапроамида // Высокомол. соединения. 1983. Т. Б.25. № 5. С. 343-346.

50. Бегишев В.П. Совмещенные процессы при синтезе кристаллизующихся полимеров. Автореферат на соискание ученой степени доктора химических наук. М. 1985.26 с.

51. Кипин И.А. Макрокинетика кристаллизации полимеров и ее приложение к моделированию процесса получения изделий из поликапроамида. Канд. дисс. на соискание учено степени кандидата физико-математических наук. Пермь. 1984. 138 с.

52. Малкин А .Я., Кипин И.А., Болгов С.А., Бегишев В.П. Альтернативный способ описания кинетики кристаллизации // Инж. физ. ж. 1984. Т. 46. № 1. С. 124-128.

53. Столин A.M., Малкин А .Я., Мержанов А.Г. Неизотермические процессы и методы исследования в химии и механике полимеров // Успехи химии. 1979. Т. 48. Вып. 8. С. 1492-1517.

54. Беляева Н.А. Исследование влияния неизотермических условий синтеза на некоторые механические свойства полимеров. Канд. дисс. на соискание ученой степени кандидата физико-математических наук. Черноголовка. 1990. 164 с.

55. Савада X. Термодинамика полимеризации. Пер. с англ. / Под. ред. Берлина А.А. М.: Химия. 1979. 312 с.

56. Keller A., Wilmouth F.M. On the morphology and origin of the fibers observed in nascent Ziegler polymethylene//Macromol. Chem. 1969.V. 121. P. 42-51.

57. Wunerlich B. Crystallizathion during polymerization // Adv. Polymer Sci. 1968/V. 5. P. 568-619.

58. Фрунзе Т.М., Шлейфман Р.В., ГодовскийЮ.К., Генин Я.В., Волкова Т.В., Котельников В.А., Курашев В.В. Некоторые особенности кристаллизации поликапроамида в процессе анионной адиабатической полимеризации // Докл.АН СССР. 1981. Т. 257. № 1. С. 162-166.

59. Болотин B.B., Болотина K.C. Расчет остаточных напряжений и деформаций в намоточных изделиях из армированных пластиков // Мех. полим. 1969. № 1. С. 134-139.

60. Тарнопольский Ю.Н., Портнов Г.Г. Изменение усилия натяжения при намотке изделия из стеклопластиков // Мех. полим. 1966. № 2. С. 278-184.

61. Николаев В.П., Инденбаум В.М. К расчету напряжений в намоточных изделих из стеклопластиков//Мех. полим. 1970. № 6. С. 1026-1030.

62. Портнов Г.Г., Горошкин В.А., Тилюк А.Г. Начальные напряжения в кольцах из стеклопластиков, изготовленных намоткой // Мех. полим. 1969. № 3. С. 505-511.

63. Благонадежин В.А., Перевозчиков В.Г. Остаточные напряжения в кольцах из стеклопластиков, полученных методом послойного отверждения // Мех. полим. 1972. № 1. С. 174-176.

64. Благонадежин В. А., Мишенков Г.В., Николаев В.П. Результаты экспериментального исследования остаточных напряжений в намоточных изделиях из стеклопластиков // Механика полимеров. 1970. № 6. С. 11161119.

65. Болотин В.В. Влияние технологических факторов на механическую надежность конструкций из композитов // Мех. полим. 1972. № 3. С. 529-540.

66. Арутюнян Х.А., Давтян С.П., Розенберг Б.А., Ениколопян Н.С. Отверждение эпоксидного олигомера ЭД-5 аминами в режиме распространения фронта реакции // Докл. АН СССР. 1975.Т. 223. С. 657-660.

67. Болотин В.В., Воронцов А.Н., Мурзыханов Р.Х. Анализ технологических напряжений в намоточных изделиях из композитов на протяжении всего процесса изготовления // Механика композитных материалов. 1980. № 3. С. 500-508.

68. Болотин В.В., Новиков Ю.Н. Механика многослойных конструкций. М. 1980. 375 с.

69. Клычников Л.В., Давтян С.П., Худяев С.И., Ениколопян Н.С. О влиянии неоднородного температурного поля на распределение остаточных напряжений при фронтальном отверждении // Механика композитных материалов. 1980. № 3. С. 509-513.

70. Клычников Л.В., Давтян С.П., Худяев С.И., Ениколопян Н.С. Трехслойная сферическая оболочка, полученная фронтальным отверждением // Механика композитных материалов. 1981. № 3. С. 494-498.

71. Кострицкий С.Н., Циркин М.З. К определению механических характеристик полимерных материалов при повышенной температуре // Заводлаб. 1980. Т. 46. № 6. С. 1087-1092.

72. Розенберг Б.А., Ениколопян Н.С. Проблемы технологической монолитности изделий из композитных материалов // Журн. Всесоюзн. хим.об-ва. 1978. Т. 23.№3. С. 298-304.

73. Томашевский В.Т., Шалыгин В.Н., Яковлев B.C. Моделирование условий возникновения технологических дефектов в структуре намоточных композитных полимерных материалов // Механика полимеров. 1980. № 5. С 895-900.

74. Аванесов K.JL, Бугаков И.И., Екельчик B.C. Расчет технологических напряжений в оболочках из стеклопластиков по теории термовязкоупругости. В кн. Исследование по теории упругости и пластичности. JI. 1976.Вып 2. С. 146-151.

75. Бахарев С.П., Миркин М.А. К вопросу о температурных напряжениях, возникающих в процессе термообработки цилиндрических оболочек из стеклопластиков//Механ. полим. 1978. № 6. С. 1118-1121.

76. Благонадежин В.А., Перевозчиков B.JL, Меркулов В.Д., Поляков B.J1. Механические свойства углепластика и остаточные напряжения в намоточных изделиях из комбинированных композитов // Механ. полим. 1975. №6. С. 996-1004.

77. Екельчик B.C., Исмайлов B.C., Кострицкий С.Н., Сборовский А.К. Описание температурной зависимости вязко-упругих свойств полимеров в динамическом режиме на основе линейной теории термовязкоупругости // Механ. полим. 1977. № 5. С. 915-918.

78. Коротков В.Н., Турусов Р.А., Розенберг Б.А. Температурные напряжения в цилиндре из композитного материала в процессе его охлаждения и хранения // Механика композитных материалов. 1983. № 2. С. 290-295.

79. Коротков В.Н., Чеканов Ю.А., Розенберг Б.А. Изотермическое неоднородное отверждение цилиндрических изделий из полимерных композитных материалов // Механика композитных материалов. 1988. № 5. С. 872-877.

80. Томашевский В.Т., Наумов В.Н., Шалыгин В.Н. Влияние неоднородности охлаждения на технологические напряжения в толстостенных цилиндрических оболочках из армированных полимеров. В кн. Полимерные материалыв машиностроении. Пермь. 1977. № 214. С. 10-18.

81. Болотин В.В., Болотина К.С. Технологические напряжения и трансверсальная прочность армированных пластиков. В сб. к 70-летию академика АН УкрССР Серенсена: Киев. 1975. С. 100-110.

82. Болотин В.В., Воронцов А.Н. Образование остаточных напряжений в изделиях из слоистых и волокнистых композитов в процессе отверждения // Механ. полим. 1976. № 5. С. 790-795.

83. Болотин В.В. К теории вязкоупругости для структурно-неустойчивых материалов // Труды Моск. энергет. ин-та. Динамика и прочность машин. 1972. Вып. 101. С. 7-14.

84. Болотин В.В., Болотина К.С. Об усадке эпоксидных связующих в процессе отверждения // Механ. полим. 1972. № 1. С. 178-181.

85. Малмейстер А.К., Тамуэн В.П., Тетере Г.А. Сопротивление жестких полимерных материалов. Рига. 1972. 498 с.

86. Болотин В.В., Болотина К.С. Результаты комплексных экспериментальных исследований кинетики отверждения эпоксидных связующих // Труды Моск. энергет. ин-та. Динамика и прочность машин. 1972. Вып. 101. С. 15-23.

87. Афанасьев Ю.А., Екельчик B.C., Иванов В.К. Регулирование остаточных напряжений в толстостенных ортотропных цилиндрических изделиях. В кн. Всесоюзн симпозиум "Новые полимерные композиционные материалы в машиностроении". Тез. докл. М. 1978. С. 25-29.

88. Коротков В.Н., Турусов Р.А., Джавадян Э.А. Розенберг Б.А. Технологические напряжения при отверждении цилиндрических изделий из полимерных композитных материалов // Механика композитных материалов. 1986. № 1. С. 118-123.

89. Арутюнян Н.Х., Дроздов А.Д., Наумов В.Э. Механика растущих вязко-упруго-пластичесих тел. М.: Наука. 1987. 471 с.

90. Турусов Р.А., Метлов В.В. Формирование напряжений при фронтальном отверждении композитов // Механика композит, материалов. 1985. № 6. С. 1079-1085.

91. Тринчер В.К. О постановке задачи определения напряженно-деформированного состояния растущего тела // Изв. АН СССР. Механика твердого тела. 1984. № 2.С. 110-124.

92. Арутюнян Н.Х., Метлов В.В. Нелинейные задачи теории ползучести наращиваемых тел, подверженных старению // Изв. АН СССР. Механика твердого тела. 1983. № 4.С. 142-152.

93. Метлов В.В., Турусов Р.А. О формировании напряженного состояния вязкоупругих тел, растущих в условиях фронтального отверждения // Изв. АН СССР. Механика твердого тела. 1986. № 6.С. 145-160.

94. Турусов Р.А., Давтян С.П., Шкадинский К.Г., Розенберг Б.А., Андреевская Г.Ф., Ениколопян Н.С. Механические явления в условиях распространения фронта отверждения // ДАН СССР. 1979. Т. 247. № 1. С. 97100.

95. Турусов Р.А., Розенберг Б.А., Ениколопян Н.С. О формировании напряжений и разрывов в процессе фронтального отверждения // ДАН СССР. 1981. Т. 260. № 1.С. 90-94.

96. Клычников JI.B., Давтян С.П., Турусов Р.А., Худяев С.И., Ениколопян Н.С. Влияние упругой оправки на распределение остаточных напряжений при фронтальном отверждении сферического образца // Механика композит, материалов. 1980. № 2. С. 300-303.

97. Чечило H.M., Ениколопян Н.С. О роли давления и начальной температуры реакционной смеси при распространении реакции полимеризации // Докл. АН СССР. 1976. Т. 230. № 3. С. 160-163.

98. Беляева Н.А., Клычников JI.B., Давтян С.П., Худяев С.И. Численный анализ напряженного состояния в цилиндрических изделиях, полученных фронтальным отверждением: Препринт ОИХФ АН СССР. Черноголовка, 1987.36 с.

99. Беляева Н.А. Моделирование процесса накопления внутренних напряжений при термической обработке древесины // Республиканская конф. Коми НЦ Уро АН СССР: Тез. докл. Сыктывкар, 1987. С.101.

100. Беляева Н.А. Анализ напряженного состояния изделия цилиндрической формы // Республиканская конф. "Разработка и внедрение средств автоматизации в сферу науки, пр-ва, управления": Тез. докл. Сыктывкар, Коми НЦ Уро АН СССР, 1988. С.18.

101. Беляева Н.А., Клычников JI.B., Давтян С.П., Худяев С.И. Образование внутренних напряжений в процессе получения изделий цилиндрической формы // Механика композитных материалов. 1988. № 6. С. 1060-1068.

102. Беляева Н.А. Исследование влияния неизотермических условий синтеза на некоторые механические свойства полимеров: Дисс. канд. физ.-мат. наук. Черноголовка, 1990. 164 с.

103. Беляева Н.А. Исследование влияния неизотермических условий синтеза на некоторые механические свойства полимеров: Автореферат дисс. канд. физ.-мат. наук. Черноголовка, 1990. 22 с.

104. Холпанов Л.П., Закиев С.Е., Помогайло А.Д. Моделирвание фронтальной полимеризации с учетом фазовых переходов во фронте реакции // ДАН, 2004. Т. 395. № 2. С. 211-214.

105. Bartczak Z. // Appl. Polym. Sci. 2002. V. 86. P. 1396.

106. Ошмян В.Г., Патлажан C.A., Remond Y. Принципы структурно-механического моделирования полимеров и композитов // Высокомолек. соед., Серия А, 2006. Т. 48. № 9. С. 1-12.

107. Матвеенко В.П., Шардаков И.Н., Бегишев В.П., Писцов Н.В. Simulation of thermomechanical processes in crystallizing polymer // Int. J. Polymer Engineering and Science. 1997. V. 37. № 8. P. 1270-1279.

108. Беляева Н.А. Совмещенные процессы и напряженное состояние формируемого цилиндрического изделия // XI Коми Респ. молодежная научная конференция: Тез. докл. Сыктывкар, 1990.

109. Беляева Н.А., Клычников Л.В., Давтян С.П., Худяев С.И. Напряженное состояние при формировании.цилиндрических изделий в ходе параллельного протекания полимеризации и кристаллизации // Механика композитных материалов. 1991. № 6. С. 1091-1099.

110. Беляева Н.А., Худяев С.И. Регулирование процесса формирования цилиндрического изделия в ходе отверждения // 2-я международная научно-техническая конференция "Актуальные проблемы фундаментальных наук": Материалы. Москва, январь 1994. Т.5. 4 с.

111. Беляева Н.А., Беляев Ю.Н. Исследование напряженного состояния в ходе сушки древесины // 2-я международная научно-техническая конференция "Актуальные проблемы фундаментальных наук": Материалы. Москва, январь 1994. Т. 1. С. 74-76.

112. Беляева Н.А., Беляев Ю.Н. Формирование цилиндрического изделия в ходе регулируемого процесса отверждения // Симпозиум "Математические модели деформируемого тела": Материалы. СПб, июнь 1994. 8 с.

113. Беляева Н.А., Беляев Ю.Н. Регулирование уровня внутренних напряжений формируемого изделия в ходе совмещенного процесса полимеризации и кристаллизации // Вестник Сыктывкарского университета. Вып. 1. 1995. С.119-127.

114. Беляева Н.А., Клычников JI.B. Метод интегрального уравнения в задаче объемного отверждения // Вестник Сыктывкарского университета. Вып.2. 1996. С.125-134.

115. Беляева Н.А., Титова Н.Н., Юркина Е.А. Модель фронтального отверждения цилиндрического изделия с точки зрения механики непрерывно наращиваемого твердого тела // Труды Сыктывкарского лесного института. Т. 1. 1997. С. 44-52.

116. Беляева Н.А., Паршукова Н.М. Термовязкоупругая модель отверждения сферического изделия // Вестник Сыктывкарского университета. Серия 1. Вып.4. 2001. С.115-124.

117. Беляева Н.А. Вязкоупругая модель отверждения композитного изделия // Труды XXXVI Уральского семинара "Механика и процессы управления". Екатеринбург, 2006. Т. 1. С. 228-235.

118. Беляева Н.А., Осипова В.В. Формирование цилиндрического изделия в ходе объемного отверждения // Федеральное агентство по образованию. ОФАП. Свид. об отрасл. регистрации разработки. № 7944. 30.03 2007.

119. Беляева Н.А. Математические модели деформируемых структурированных материалов. Монография. Сыктывкар: Изд-во СыктГУ, 2008. 116 с.

120. Беляева Н.А. Структуризация в процессах деформирования вязкоупругих систем // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. Ижевск: Удмуртский госуниверситет. 2008. Вып. 2. С. 178-179.

121. Беляева Н.А. Структурные модели процессов деформирования вязкоупругих материалов // Вестник Сыктывкарского университета. Сер 1. Вып. 8. 2008. С. 55-68.

122. Франк-Каменецкий Д.А. Диффузия и теплопередача в химической кинетике. М. 1987. 491 с.

123. Илюшин А.А. Механика сплошной среды. М.: Изд-во МГУ. 1990. 310 с.

124. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред. М.: Наука. 1982. 620 с.

125. Работнов Ю.Н. Механика деформируемого твердого тела. М.: Наука. 1979. 744 с.

126. Ван-Кревелен Д.В. Свойства и химическое строение полимеров. М. 1966.414 с.

127. Жирков П.В. Исследование неизотермических течений реагирующих жидкостей с изменяющейся в процессе химической реакции вязкостью. Автореф. дис. канд. физ.-мат. наук. Черноголовка. 1979. 16 с.

128. Тимошенко С.П., Гудьер Дж. Теория упругости. Пер. с англ. Под ред Шапиро Г.С. М.: Наука. 1975. 576 с.

129. Encyclopedia Of Polymer Science and Engineering. Second Edition. Scattering to Structural Foams. A Wiley-Interscience Publication. USA. 1985. V. 15. P. 346-362.

130. Баронин Г.С., Кербер M.JI., Минкин E.B., Радько Ю.М. Переработка полимеров в твердой фазе. Физико-химические основы. М: Машиностроение -1.2002. 320 с.

131. Merzhanov A.G., Stolin A.M., Buchatskii L.M.et al. A ceramic composite material and method of its production in Russian. European Patent №.89910469. 9 of April 15. 1991.

132. Подлесов B.B., Столин A.M., Мержанов А.Г. СВС-экструзия электродных материалов и их применение для электроискрового легирования стальных поверхностей // ИФЖ. 1993. Т. 63. № 5. С. 636-647.

133. Подлесов В.В., Радугин А.В., Столин A.M., Мержанов А.Г. Технологические основы СВС-экструзии // ИФЖ. 1992. Т. 63. № 5. С. 525537.

134. Стельмах Л.С., Столин A.M. Математическое моделирование СВС-экструзии. 4.2. Реодинамические модели // ИФЖ. 1993. Т. 64. № 3. С. 90-94.

135. Стельмах Л.С., Столин A.M., Мержанов А.Г. Макрореологическая теория СВС-компактирования // Докл. РАН. 1995. №1. Т. 344. С. 72-77.

136. Maizeliya A.V., Stolin A. M. and Stelmakh L. S. Hydrodynamic analysis of the process of SHS disintegration // Int. J. of the high-temperature synthesis. 1996. V. 5. №.2. P. 145-151.

137. Stelmakh L.S. and Stolin A. Macrorheological theory of high-temperature deformation of powder materials under SHS-compacting // J. of the European ceramic society. 17(1997). P. 453-456.

138. Стельмах Л.С., Жиляева H.H., Столин A.M. Неизотермическая .реодинамика при СВС-прессовании порошковых материалов // ИФЖ. 1991. Т. 61. № 1.С. 33-40.

139. Стельмах Л.С., Столин A.M. О квазистационарном режиме и предельных случаях горячей экструзии порошковых материалов // Докл. РАН. 1992. Т. 322. №4. С. 732-736.

140. Стельмах Л.С., Жиляева Н.Н., Столин A.M. Реодинамика и теплообмен горячего компактирования порошковых материалов // ИФЖ. 1992. Т. 63. №5. С. 612-622.

141. Stelmakh L.S, Zinenko Zh.A., Stolin A.M., Merzhanov A.G. Specific features of SHS material compacting hydrodynamic and thermal action // Int. J. of SHS. 1995. V. 4. № 1. P. 263-273.

142. Скороход В.В. Реологические основы спекания. Киев: Наук, думка. 1972. 152 с.

143. Ковальченко М.С. Теоретические основы горячей переработки пористых материалов давлением. Киев. 1980. 240 с.

144. Горохов В.М., Ковальченко- М.С., Роман О.В. Уплотнение и формоизменение пористых материалов при горячем прессовании в условиях неравномерного трехосного сжатия // Порошковая металлургия. 1983. № 1. С. 8-13.

145. Горохов В.М., Дорошкевич Е.А., Звонарев Е.В. Уплотнение и формоизменение пористых материалов при горячем прессовании в условиях неоднородного напряженно-деформированного состояния // Порошковая металлургия. 1984. № 5. с. 100-105.

146. Скотников М.В., Чулков В.П., Прилепкин Б.П., Джангирян В.Г. Распределение напряжений и скоростей при уплотнении вязких тел в замкнутых объемах // Порошковая металлургия. 1984. № 6. С. 21-27.

147. Седов Л.И. Механика сплошной среды. Т.1, 2. М.: Наука. 1983, 1984. 528 е., 560 с.

148. Бучацкий Л.М., Столин A.M., Худяев С.И. Кинетика изменения распределения плотности при горячем прессовании вязкого пористого тела // Порошковая металлургия. 1986. № 2. С. 37-42.

149. Столин A.M., Жиляева Н.Н., Хусид Б.М. Регулярный режим уплотнения горячих порошковых материалов // ИФЖ. 1990. Т. 59. № 2. С. 248-254.

150. Тадмор 3., Гогос К. Теоретические основы переработки полимеров. Пер. с англ. под ред. Торнера Р.В. М.: Химия. 1984. 632 с.

151. Смирнов-Аляев Г.А. Сопротивление материалов пластическому деформированию. Л. 1978.

152. Гессингер Г.Х. Порошковая металлургия жаропрочных сплавов. Челябинск. 1988.

153. Стельмах Л.С., Столин A.M., Мержанов А.Г. Математическое моделирование СВС-экструзии. 4.1. Тепловые модели // ИФЖ. 1993. Т. 64. № 3.С. 83-89.

154. Stelmakh L.S. and Stolin A.M. Macrorheological theory of hot compaction of composites // Mechanics of Composite Materials. 1995. V. 31. № 6. P. 840845.

155. Стельмах Л.С., Столин A.M., Хусид Б.М. Реодинамика выдавливания вязких сжимаемых материалов // ИФЖ. 1991. Т. 61. № 2. С. 268-276.

156. Баронин Г.С., Кербер М.Л., Радько Ю.М., Минкин Е.В. Переработка полимеров в твердой фазе. Физико-химические основы. М.: Машиностроение-1. 2002. 320 с.

157. Сверхмодульные полимеры: Пер. с англ. / Под ред. Чифери А., Уорда И., Малкина А .Я. Д.: Химия. 272 с.

158. Славнов Е.В., Скульский О.И., Тимофеев В.М. Ориентационные эффекты в образцах после твердофазной гидроэкструзии при высокотемпературном нагружении. // Пластические массы. 1997. № 8. С. 18-22.

159. Смолев Л.В., Беляева Н.А. Экструзия с заданным усилием на плунжере. М.: Федеральное агентство по образованию. ОФАП. Свид. об отрасл. регистрации разработки № 7945. 30.03 2007.

160. Беляева Н.А., Столин A.M., Стельмах Л.С. Динамика твердофазной плунжерной экструзии вязкоупругого структурированного материала // Теоретические основы химической технологии, 2008. № 5. 579-589.

161. Belyaeva N. A., Stolin А. М., Stelmakh L. S. Dynamic of Solid-State Extrusion of Viscoelastic Cross-Linked polymeric Materials // Theoretical Foundations of Chemical Engineering, 2008. Vol. 42. No 5. P. 549-556.

162. Успехи реологии полимеров. Под ред. Виноградова Г.В. М.: Химия. 1970. 294с.

163. Виноградов Г.В., Малкин А.Я. Реология полимеров. М.: Химия. 1977. 438 с.

164. Тагер А.А. Физикохимия полимеров. М.: Химия. 1978. 543 с.

165. Беляева Н.А., Столин A.M., Стельмах Л.С. Влияние структуризации на процесс твердофазной экструзии вязкоупругого материала. Материалы

166. Международной конференции «XVIII сессия Международной Школы по моделям механики сплошной среды». Изд-во Саратовского госуниверситета: Саратов. 2007. С. 55-58.

167. Беляева Н.А., Столин A.M., Стельмах J1.C. Кинетика уплотнения и структуризации в твердофазной экструзии вязкоупругой среды // Инженерная физика. 2007. № 5. С. 34-41.

168. Беляева Н.А. Твердофазная экструзия с условием постоянства скорости плунжера пресса. М.: Федеральное агентство по образованию. ОФАП. Свид. об отрасл. регистрации разработки № 7946. 30.03 2007.

169. Беляева Н.А., Столин A.M., Пугачев Д.В., Стельмах JI.C. Неустойчивые режимы деформирования при твердофазной экструзии вязкоупругих структурированных систем // ДАН, 2008. Т. 420. № 6. С. 777-780.