CC BY
62
Э. Д. Шакирьянов, М. Н. Фаттахов, С. М. Усманов,
Ю. М. Сивергин, Н. В. Улитин, Т. Р. Дебердеев
ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ БЛОЧНОЙ ТРЕХМЕРНОЙ
СВОБОДНО-РАДИКАЛЬНОЙ ПОЛИМЕРИЗАЦИИ ДИАЛЛИЛИЗОФТАЛАТА
Ключевые слова: диаллилизофталат, имитационное моделирование, трехмерный полимер, метод Монте-Карло.
Исследована структурная неоднородность макротела трехмерного полимера диаллилизофталата (ДАИФ) с помощью имитационного моделирования полимеризационного процесса в блоке с множественным инициированием. Получено численное распределение геометрических элементов свободного объема (ЭСО) по размерам. Представлена методика и определены прочностные характеристики поли-ДАИФ.
Keywords: diallylisophtalate, simulation modeling, three-dimensional polymer, Monte-Carlo simulation.
Structural heterogeneity of three-dimensional macro body of diallylisophtalate (DAIP) polymer is investigated by means of simulation modeling of polymerization process in the block with multiple initiation. Numerical distribution of geometrical elements of the free volume (FVE) by the sizes is received. Strengthening characteristics ofpoly-DAIP are defined and technique is provided
Введение
Задача синтеза материалов из трехмерных полимеров с заранее известными прочностными свойствами и по сей день является весьма актуальной. На данный момент наиболее целесообразным подходом представляется получение материалов с априори известными характеристиками исходя из модельных представлений об их структуре. Работ данном направлении огромное множество [1-4]. Решение этой проблемы требует не только глубоких знаний кинетики полимеризации полифункциональ-ных мономеров и олигомеров, но и проведение трудоемких и затратных экспериментов с помощью дорогостоящего лабораторного оборудования.
Одним из перспективных методов исследования данного вопроса является создание математической модели, позволяющей имитировать процесс трехмерной свободно-радикальной полимеризации полифункциональных мономеров с различными начальными условиями до предельных степеней конверсии. Свободно-радикальная полимеризация протекает по механизму классической радикальной полимеризации [5], однако при моделировании необходимо учесть тот фактор, что протекание реакционного процесса сопровождается изменением «физического» состояния реакционной массы от жидкого до стеклообразного. Исследователями кинетики процесса трехмерной полимеризации показано резкое падение скорости полимеризации с увеличением конверсии, что объяснялось низкой молекулярной подвижностью в реакционной зоне, связанной с отсутствием достаточного свободного объема [6]. В настоящей работе предложены результаты математического моделирования на основе метода Монте-Карло процессов блочной свободно-радикальной полимеризации тетрафункциональных мономеров и олигомеров.
Особенности кинетики полимеризации ДАИФ
Отличительной особенностью кинетики полимеризации ДАИФ является наличие механизма передачи реакции цепи на мономер. ДАИФ относит-
ся к группе моноаллиловые мономеров общей формулы СН2=СН-СН2-Х, которые полимеризуются по радикальному механизму с очень низкими скоростями и образуют олигомерные продукты (Рп =3 30 ) [7, 8]. Характерной чертой аллиловых мономеров является присутствие в молекулах последних СН2-групп в «-положении к двойной связи.
При атаке молекулы диаллилового мономера активным радикалом наряду с присоединением по двойной связи с высокой вероятностью может произойти отрыв атома водорода от а-СН2-группы. Это приводит к гибели исходного активного радикала и к образованию малоактивного аллильного радикала, стабилизированного внутренним напря-
жением:
—CHxCH ■
CH2 I 2
XR
CH2=CH—CH2
2 I 2
XR
—CH2CH2 ■ CH2
I
XR
CH2=CH—CH XR
СН2= ОИ-СИ ХЯ СИг еи = снхи. ,
Аллильные радикалы могут активно вступать в реакцию с двойной связью аллилового мономера только при температурах выше 300 + 400ОС или под давлением в несколько килобар (5 + 20) [8, 9]. В обычных условиях их активность близка к нулю и обычно, они рекомбинируют. Поэтому отрыв а-водорода растущим радикалом означает гибель как материальной, так и кинетической цепи, что приводит к крайне низким скоростям полимеризации диаллиловых мономеров и образованию продуктов с невысокими молекулярными массами. Такой процесс носит название деградационной передачи цепи [7-8].
deg
->RnH + А-,
/1п-+М-
где кавд - константа скорости реакции деградационной передачи цепи; А • - малоактивный аллиль-ный радикал.
Наряду с деградационной передачей цепи существует эффективная передача, после которой сохраняется возможность продолжения цепной реакции.
А • +М——
где кеП - константа скорости реакции эффективной передачи цепи
Описание модели
В работе была использована модифицированная модель [10-12] процесса трехмерной свободно-радикальной полимеризации ДАИФ в блоке до предельных степеней превращения, основанная на методе Монте-Карло. Применение данного метода основано на случайном характере протекания элементарных реакций процесса полимеризации. В основу модели реакционной системы положен виртуальный реактор, построенный на кубической решетке 100х100х100. Рост полимерной цепи происходит по ребрам, соединяющим соседние узлы, с запретом посещения ранее занятых узлов (метод исключенного объема). Модификации было подвергнуто содержимое сосредоточенного в узлах решетки микрообъема.
В раннее опубликованных работах [10-12] мономерное звено в цепи полиДАИФ занимало 3 узла, при этом предполагалось, что крайние узлы содержат реакционноспособные атомные группы, а в среднем узле сосредоточены инертные атомные группы. Однако в таком случае узлы в кубической решетке модельного реактора сильно неравноценны по сосредоточенной в них массе и занимаемому микрообъему.
В связи с этим в модели были внесены небольшие изменения: на каждое мономерное звено в полимерной цепи было отведено 14 узлов кубической решетки, согласно числу углерод-содержащих атомных групп в структурной формуле. Несмотря на то, что микрообъемы и массы углерод-содержащих атомных групп в молекуле ДАИФ различаются между собой, в этом случае дисбаланс менее ярко выражен. А увеличение количества узлов до числа уг-лерод-содержащих атомных групп, приходящихся на одну молекулу мономера, позволяет рассмотреть реакционный процесс с более высокой степенью детализации ив то же время в относительной мере несколько уравнять все узлы в решетке реактора как по сосредоточенной в них массе, так и по занимаемому микрообъему. В дальнейшем такая детализация позволила получить ряд интересных результатов.
Предполагалось, что в реакционной системе отсутствуют какие-либо примеси, способные вступать в реакцию с мономером или радикалом, а стенка модельного реактора является инертными.
Имитация роста полимерной цепи в ходе моделирования реакционного процесса предполагала заполнение узлов решетки атомными группами присоединяемой молекулы мономера. При этом помечались узлы, содержащие атомные группы подвешенных двойных связей и радикалов.
При распаде инициатора в модели считалось, что образующийся первичный радикал практически мгновенно вступает в реакцию с мономером или с не прореагировавшей двойной связью полимерной цепи, образуя полимерный радикал 13'. Таким образом, инициирование в модели лимитирова-
лось распадом инициатора и рассматривалось как один независимый элементарный процесс.
В отличие от линейной радикальной полимеризации, трехмерная полимеризация имеет ряд особенностей. В целом процесс полимеризации ДАИФ протекает по классической кинетической схеме радикальной полимеризации с учетом механизма реакции передачи цепи на мономер.
В ходе роста полимерной цепи в реакционной системе возможны не только реакции с участием молекул мономера, но и с боковыми (подвешенными) двойными связями О ! полимерных цепей.
В результате, в системе образуется так называемый «менее активный радикал» К , константа скорости реакции роста и обрыва для которого из-за стериче-ских ограничений принята на порядок меньше [12]. В зависимости от того с двойной связью какой цепи прореагировал полимерный радикал в системе могут образовываться циклы и сшивки.
Таким образом, с учетом особенностей кинетической схемы полимеризационного процесса состояние реакционной системы в модели определялось следующими элементарными реакциями, скорости представлены ниже:
И' = 2рка/,
УУ = к КМ ,
р р ’
"Р =
УУ'=0, \kR~M, (1)
р ' р ’ V >
"Р =0 = кд2,
/ =0,1 кД'Я\
УУ" =0,0\ktR R', где / , М , С>р, ^, R' - число молекул инициатора и мономера, число двойных связей, радикалов и «менее активных радикалов» в системе, ка, кр,
к^, к1 - численные значения констант распада, скорости роста, реакции с двойной связью Оро/ и обрыва соответственно; УУ. - скорость инициирования, УУ и УУр0 - скорости реакций взаимодействия полимерного радикала с мономером и двойной связью, УУ и УУр0 - скорости реакций взаимодействия «менее активного радикала» с мономером и двойной связью, УУп УУ1 и Щ - соответствующие скорости реакций обрыва; р - эффективность инициатора,
определяющая долю первичных радикалов, образующихся при распаде инициатора, обычно лежит в пределах от 0,5 до 1.
В модели механизм передачи был реализован с помощью следующих элементарных реакций:
- деградационная передача цепи на мономер
(2)
- эффективная передача цепи на мономер
(3)
где А - число аллильных радикалов в системе, И^д и ]Л/еП - скорости реакций деградационной и
эффективной передачи цепи на мономер, кдед и
ке^ - константы соответствующих элементарных
*
реакций .
Кроме этого, дополнительно было учтено участие аллильного радикала в реакциях гибели радикалов К и /?' . Для расчета модельных значений константы скорости реакции передачи цепи на мономер было принято км = кд + кд = 20 л/(мольс)
[8, 9, 13]. В модели все скорости приводятся в размерности количество реакционных актов за единицу времени (с-1).
Полимеризационный процесс в модели представляет собой набор дискретных во времени состояний реакционной системы, где интервал времени А! между дискретными состояниями определяется как:
1 (4)
М =
IX
- тип реакции.
Имитирование элементарных реакций проводится согласно условиям:
WJ, (5)
где ^ - случайные равномерно распределенные
числа в интервале [0, 1).
Генерация случайного числа проводится для каждого типа реакции. В случае, если выполняется /-ое условие (5), в системе запускается соответствующая /-ая элементарная реакция. После проведения элементарной реакции состояние реакционной системы пересчитывается, вычисляется новый шаг Д/ и процесс повторяется в цикле. Такой подход позволяет рассматривать в определенных случаях элементарные реакции в модели как набор квазине-зависимых параллельных процессов. Завершение работы алгоритма осуществляется либо по достижении определенного времени, заданной степени конверсии или по другому заданному условию. Одной из особенностей данной модели является возможность отражать реакционные процессы в системе с учетом диффузионных ограничений. Это позволяет с большей уверенностью моделировать процесс до предельный степеней конверсии.
Известно [14], что рост полимерной цепи допускает непрерывные процессы реакции диффузии и химической реакции, в результате чего общие константы реакций роста кр и обрыва ^ имеют две составляющие:
1,1,1
к к° кпП
р рк ри
(6)
Здесь значения констант, а также скоростей элементарных реакций имеют размерность с-1
1 1
—+—
(7)
где крН и кт - химический вклад, кр0 и кю -
диффузионный вклад, который согласно Хэмилеку [15] зависит от свободного объема системы и, следовательно, от степени конверсии, но не от молярной массы полимера:
'1 1
кр = С-ехР
кю=Ко-ех Р
-V,
(8)
(9)
1-А.
/ 1
где - доля свободного объем системы и - доля свободного объема при нулевой конверсии или доля свободного объема, приходящаяся на мономер, к0рО
и к^а - диффузионные вклады в константу роста и обрыва при нулевой степени конверсии, Vр и ^ -факторы перекрытия свободных объемов.
Химические вклады крК и кт определяются из условия нулевой конверсии, из которого следует кР=кр к, = к° и /■ = ^ - отсюда кр0 = к0рО,
* ~к1
'ю
'ю
, и, следовательно,
1 крн 1 Ч0" 1 к0 ро
1 кю 1 = о- 1 к . ю
В итоге для констант св
рыва получаем: и і
р Ґ
1 1 ехр V
° - к 0 ро к
к - 1
1 1 ехр Л •
1-А. / ґ
(10)
(її) и об-
(12)
рО
1-А. / ґ
(13)
ч ю ю
Согласно [16] состояние свободного объема в полимерной системе при температуре опреде-
ляется соотношением:
1 = с[ТдР - Аар (Т - ТГр)] + -1 - с)[ТдМ + Дам (т - ТдМ)] , (14)
где f р, - доля свободного объема при темпера-
туре стеклования полимера и мономера соответст-
венно, Дар = а|Р -адР, Дам = а|М -адМ - разность
между коэффициентами теплового расширения выше и ниже температуры стеклования полимера ТдР
и мономера ТдМ , с - доля полимера в реакционной
системе.
Значение 1 для мономера и полимера определялось согласно [14, 15] соотношением:
1
ДаТ =Пп
д д
(15)
При моделировании полимеризационного процесса в системе в зависимости от достигнутой степени конверсии на каждом шаге производился пересчет констант скоростей реакций роста и обрыва, что отражалось на значениях скоростей (1). Более подробно процесс моделирования кинетики полимеризации изложен в работах [10-12]. В таблице приведены исходные параметры реакционной системы с учетом описанных выше модификаций.
Таблица 1 - Параметры реакционной системы полимеризации ДАИФ
т!
/ М V, л
0,05 3617 71170 1,76е-20
т1
где —— - массовое отношение инициатора и мономера, V - эквивалентный объем виртуального реактора
Результаты численного эксперимента
Модельная система позволяет при проведении численного расчета отключать механизмы передачи цепи на мономер и диффузионного торможения. На рис.1 и 2 показаны зависимости конверсии от времени, полученные в разных режимах, и сопоставление их с данными эксперимента [13]. Кривая (1), изображенная на обоих рисунках, получена с отключенными механизмами передачи цепи и диффузионного торможения (ошибка ). Как
видно из рисунков, она качественно описывает динамику процесса и завышенные значения степени конверсии вполне очевидны. Активные полимерные радикалы, присутствующие в системе, не расходуются на образование аллильных радикалов, реакционная активность которых низка, отсутствует диффузионное торможение - в результате, процесс полимеризации не замедляется.
Рис. 1 - Расчетные конверсионные зависимости полимеризации ДАИФ: 1 - без учета передачи
цепи на мономер, 2 - с учетом передачи цепи на мономер, точками показан ход экспериментальной зависимости
Учитывая влияние реакции передачи цепи на мономер (кривая 2, рис.1), расчетная зависимость показывает хорошее совпадение с экспериментом только в начальной стадии (6^4-5%). Это можно объяснить отсутствием на этой стадии диффузионного торможения на ходе реакции. Однако при конверсиях более 30% наблюдается заметное расхож-
дение с экспериментом (б = 14-Н5%), причем в модельном процессе полимеризация протекает медленно. Однозначно объяснить такое подведение модельной кривой возможно, если провести сравнительные расчеты с экспериментальными зависимостями, полученными при различных начальных условиях полимеризации.
На рис.2 представлены расчетные конверсионные зависимости иллюстрирующие проявление диффузионного контроля на ходе моделирования реакции без учета механизма передачи. Как видно по модельным кривым 1 и 2, диффузионные эффекты начинают заметно сказываться уже в области 20% конверсии. Несмотря на сильное расхождение расчетной кривой 2 с экспериментом в начальной стадии процесса, в области высоких конверсий они сближаются.
Рис. 2 - Расчетные конверсионные зависимости полимеризации ДАИФ: 1 - без учета диффузионного эффекта, 2 - с учетом диффузионного эффекта, точками показан ход экспериментальной зависимости
Рис. 3 - Расчетные конверсионные зависимости полимеризации ДАИФ: 1 - без учета диффузионного эффекта и реакции передачи цепи, 2 - с учетом диффузионного эффекта и реакции передачи цепи, точками показан ход экспериментальной зависимости
Учет в модели обоих механизмов проиллюстрирован конверсионной зависимостью на рис.3 (кривая 2). Как видно из рисунка, полученная расчетная кривая 2 дает более качественное описание экспериментальных данных. Это позволяет использовать данную модель для исследования кинетических параметров реакционной системы при различных начальных условиях. Расхождение с экспериментом при степенях превращения более 25% вполне объяснимы несовершенством модели, ограничениями, накладываемыми кубической решеткой и неустранимыми погрешностями численного эксперимента.
Несмотря на ограниченные возможности кубической решетки, используемой в модели для имитации процессов роста и формирования полимерных цепей, ее наличие позволяет анализировать двухмерные срезы модельного реактора и получить определенные представления о структурной неоднородности полиДАИФ. Визуализация такого среза позволяет показать в реакторе незанятое пространство, часть которого может представлять собой вмороженный в матрицу полимера изолированный геометрическим микрообъемом. На рис.4. представлен один из таких срезов, где полимер обозначен темным цветом, свободное пространство - серым.
Рис. 4 - Срез модельного реактора по z = 52
Известно [17, 18], что механические свойства твердых тел определяются наличием в них свободного объема. Принято разграничивать между собой три вида свободного объема - это геометрический, флуктуационный и свободный объем, обусловленный тепловым расширением полимера [18]. Однако, экспериментальное разделение указанных видов свободного объема очень затруднено, а подчас и невозможно [17].
Классический свободный объем по Ван-дер-Ваальсу составляет около 25-30% от общего объема системы, а флуктуационный свободный объем занимает лишь 2-3% от общего объема и на порядок меньше классического свободного объема [17]. Именно флуктуационный свободный объем во многом определяет механические свойства полимерных и простых аморфных веществ.
Одним из физических параметров, который определяет механические свойства твердых тел, является коэффициент Пуассона ( ). Диапазон изме-
нения коэффициента Пуассона находится в пределах от -1 до 0,5 [19]. Причем, если ^1 = 0 , то можно утверждать, что твердое тело является пористым, оно не изменяет поперечных размеров при растяжении. Твердые тела с коэффициентом Пуассона близком к 0,5 подобны резине. Большинство других твердых тел имеют коэффициент Пуассона в пределах 0,2-
0,4. Используя связь между коэффициентом Пуассона и долей флуктуационного свободного объема, можно вычислить численное значение /У для конкретного вещества [19].
По окончании моделирования процесса блочной полимеризации ДАИФ с помощью процедуры послойного сканирования кубической решетки был проведен расчет количества и объема незанятых элементов кубической решетки - геометрических элементов свободного объема (ЭСО) «вмороженного» в матрицу полимерного тела. Суммарный объем таких ЭСО представляет собой классический свободный объем по Ван-дер-Ваальсу. Минимальный размер ЭСО составляет усредненный объем одной кубической ячейки и0, равный 17,6 А3.
На рис. 5 изображена диаграмма обнаруженных ЭСО в теле полимера. Как видно из рис. 5, в модельном теле полиДАИФ может присутствовать сравнительно большое число ЭСО, самых разнообразных по структуре и размерам, некоторые ЭСО могут достигать размеров ~ 40 ячеек, что эквивалентно расчетному объему ~ 700 А3
1> »
Общее количество микропустот в реакторе
Рис. 5 - Структурное распределение микропустот реакторе
Детальный анализ полученной диаграммы позволил построить частотную гистограмму (рис.6), иллюстрирующую состав разносортных ЭСО, заключенных внутри матрицы полимера.
Размеры микропустот в ячейках (и0=17,6 А)
Рис. 6 - Частотная шкала распределения свободного объема по количеству ячеек
В работе Сандитова Д. С. [20] отмечено, что объем среднего флуктуационного ЭСО связан с коэффициентом Пуассона // и модулем упругости следующим соотношением:
V = {'-*• К , (16)
д
где к - постоянная Больцмана, Тд - доля флуктуационного свободного объема при температуре стекло-
вания Тд (у большинства аморфных полимеров \д «0,02-5-0,03 [17]).
Из уравнения (16) найдем коэффициент Пуассона ц :
1 Vh fg E
ц =-------------------
2 6 • к Т
(17)
Значение Уь можно рассчитать из данных, полученных в ходе послойного сканирования кубической решетки виртуального реактора модели. Зная число геометрических ЭСО, значения их микрообъемов, нетрудно вычислить средний объем ЭСО. Поскольку флуктуационный свободный объем на порядок меньше классического свободного объема, то его можно рассчитать как:
■ 'ЭШ 2>
N.
-0,1,
(18)
ЭСО
где и! - объем 1-го геометрического ЭСО в кубической решетке, /ЧЭсо - число ЭСО в системе.
Подставив табличные значения модуля упругости и температуры стеклования для полиДАИФ в уравнение (17), получим /у = 0,367. Рассчитанное значение сопоставимо с коэффициентами Пуассона для полистирола и поливинилхлорида [19], что позволяет говорить о некотором сходстве прочностных свойств этих материалов.
Наличие тяжелых и крупных атомных групп в структурных звеньях полимеров обуславливает рост коэффициента Пуассона. Однако, известно, что чем выше доля флуктуационного свободного объема ^, тем выше значение коэффициента Пуассона для
аморфных полимеров.[19]. В самом деле, наличие у диаллилизофталата фенильных групп и групп СН2 в качестве боковых привесков в каждом повторяющемся звене приводит к росту нелинейности межатомных сил и флуктуационного свободного объема и, как следствие, к высокому значению /У = 0,367.
Зная коэффициент Пуассона и модуль упругости, можно рассчитать модуль сдвига [21]:
в=^-^. (19)
2(1 + ^)
Для полиДАИФ он составит значение в = 4,35 ГПа (полиэтилен - в = 0,117 ГПа).
Модуль объемного сжатия найдем по формуле [21]
К = Е ч (20)
з(1+ц)
Значение для полиДАИФ - K = 16,66 ГПа (стекло - 35-55 ГПа).
Таким образом, представленная в работе модель позволяет не только изучить кинетические характеристики полимеризационных процессов однотипных соединений, но и произвести на качественном уровне оценку механических свойств полимерных материалов.
Литература
1. Н.В. Улитин и др. Описание топологической структуры модифицированных циклокарбонатом эпоксиамин-ных систем // Вестник Казан. технол. ун-та, 2008, №5, C. 112-118.
2. Н.В. Улитин, Т.Р.Дебердеев, Р.Я.Дебердеев, Ал. Ал.Берлин, Новый подход к описанию мгновенной составляющей вязкоупругих свойств густосетчатых полимеров с применением теории наследственности, Вестник Казанского технологического университета, 2011, №14, C. 124-129.
3. Н.В. Улитин, К.А. Терещенко, Т.Р. Дебердеев, Математическое моделирование реокинетических свойств как метод оценки режима переработки изделий из эпоксидных композиционных материалов // Вестник Казан. технол. ун-та, 2012, №22, C. 17-21.
4. Н.В. Улитин и др. Кинетическая модель процесса радикальной полимеризации бутилакрилата в присутствии тритиокарбонатов // Вестник Казан. технол. ун-та, 2012, №22, C. 17-21.
5. Х.С. Багдасарян, Теория радикальной полимеризации. Наука, Москва, 1966. 299 с.
6. И.Н. Гришина, Ю.М. Сивергин, С.М. Киреева, Л.А. Жильцова, Пластические массы, 2000, №1, С. 7-9.
7. В.И. Володина, А.И. Тарасова, С.С. Спасский, Успехи химии,1970, Т. 39, №2, С. 276-303.
8. В.А. Кабанов, В.П. Зубов, Ю.Д. Семчиков, Комплексорадикальная полимеризация, Химия, Москва, 1987. 256 с.
9. А.И. Воробьева, Ю.А. Прочухан, Ю.Б. Монаков, Вы-сокомолек. соед., 2003, Т. 45, №12, С. 2118-2136.
10. Р.Р. Исмаилов, Э.Д. Шакирьянов, С.М. Усманов, Вестник УГАТУ, 2009, Т. 12, №1(30), С. 211-217.
11. М.Н. Фаттахов, Р.Р. Исмаилов, Э.Д. Шакирьянов, Ю.М. Сивергин, С.М. Усманов, Вестник Нижегородского ун-та им.Лобачевского, 2011, №3(2), С. 141-147.
12. Э. Д. Шакирьянов, Р. Р. Исмаилов, С. М. Усманов, Ю. М. Сивергин. Химическая физика, 2011, Т. 30, №9, С. 75-83.
13. I. Hace, Macromol. React. Eng., 2007, V. 1, №3, P. 391-404.
14. M. Buback, Macromol. Chem. 1990, № 191, P. 1575.
15. M. Stickler, D. Panke, A.E. Hamielec, J. Polymer Scince, 1984, V. 22, №9, P. 2243-2253.
16. F.N. Kelly, F. Bueche, J. Polym. Science, 1961 V. l, P. 549-55.
17. Д. Ферри Вязкоупругие свойства полимеров, Изд-во иностр. лит., Москва, 1963. 254 с.
18. О.Е. Ольховик, Высокомолек. соед., 1976, Т. 18(А), №5, С. 1012-1019.
19. Д.С. Сандитов, С.Ш. Сангадиев, Высокомолек. соед., 1998, Т. 40(А), №12, С.1996-2003.
20. Д.С. Сандитов, Цыдыпов Ш.Г., Баинова А.Б., Высокомолек. соед., 2004, Т. 46, №7, С. 1185-1190.
21. М.А. Колтунов, В.П. Майборода, В.Г. Зубчанинов, Прочностные расчеты изделий из полимерных материалов. Машиностроение, Москва, 1983, 239 с.
© Э. Д. Шакирьянов - к.ф.-м.н., доц., декан физико-математического факультета Бирского филиала БашГУ, eddydan@rambler.ru; М.Н. Фаттахов - асс. каф. математического моделирования и информационных систем Бирского филиала БашГУ, physlab@mail.ru; С. М. Усманов - д.ф.м.н, проф., дир. Бирского филиала БашГУ, usm@birsk.ru; Ю. М. Сивергин - д.х.н., проф. Института химической физики им. Н.Н.Семенова; Н. В. Улитин - д.х.н., проф. каф. технологии переработки полимеров и композиционных материалов КНИТТУ, n.v.ulitin@mail.ru; Т. Р. Дебердеев - д.т.н., проф. той же кафедры, deberdeev@mail.ru.
